环境系统分析PPT第4讲

1、主讲主讲: 李明俊李明俊 教授教授2006.5.8南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.第二章第二章 环境问题的模型化环境问题的模型化 一、模型与模拟一、模型与模拟 1、模型、模型 对真实系统的描述且是一种抽象。对真实系统的描述且是一种抽象。 形象模型（放大、缩小）形象模型（放大、缩小） 模型模型 抽象模型（符号、图表）抽象模型（符号、图表）南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. 模拟模型模拟模型（如电路如电路 系统模拟力学系统等系统模拟力学系统等）抽象模型抽象模型（符号、图表符号、图表） 概念模型概念

2、模型（质点、理质点、理 想气体、原子结构想气体、原子结构） 数学模型数学模型（用数学语用数学语 言表达言表达）2、模拟、模拟 模型建立后，在人为控制的条件下，通模型建立后，在人为控制的条件下，通过改变特定参数来观察分析模型的响应，预过改变特定参数来观察分析模型的响应，预测真实环境条件下实际系统的行为，这一过测真实环境条件下实际系统的行为，这一过程称为模拟。程称为模拟。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.主要采用计算机模拟，分为两大类。主要采用计算机模拟，分为两大类。 连续模拟（连续模拟（用状态连续变化的模型用状态连续变化的模型） 模拟模拟 离散

3、模拟（离散模拟（用状态不连续变化的模型，用状态不连续变化的模型， 或将连续模型离散化后进行离散模拟或将连续模型离散化后进行离散模拟） 河流水质模型可描述状态的连续变化，河流水质模型可描述状态的连续变化，但由于监测的不连续性常将问题离散化，进但由于监测的不连续性常将问题离散化，进行离散模拟。行离散模拟。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.3 3、数学模型概述（着重于系统问题的数学模、数学模型概述（着重于系统问题的数学模型）型）定义：应用数学语言和方法来描述环境污定义：应用数学语言和方法来描述环境污染过程中的物理，化学，生物化学，生物染过程中的物理

4、，化学，生物化学，生物生态以及社会等方面的内在规律和相互关生态以及社会等方面的内在规律和相互关系的数学方程。系的数学方程。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.建立过程建立过程: 它是建立在对环境系统进行反复的观它是建立在对环境系统进行反复的观察研究，通过实验或现场监测取得了察研究，通过实验或现场监测取得了大量的大量的有关信息和数据，进而对所研究的有关信息和数据，进而对所研究的系统行为系统行为动态动态、过程本质过程本质和和变化规律变化规律有了较深刻认识有了较深刻认识的基础上，经过简化和数学演绎而得出的一的基础上，经过简化和数学演绎而得出的一些数学

5、表达关系，这些表达关系描述了环境些数学表达关系，这些表达关系描述了环境系统中各变量及其参数间的关系系统中各变量及其参数间的关系。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. 应用应用在环境问题中的应用分三类：在环境问题中的应用分三类： A A、环境规划和管理中的应用环境规划和管理中的应用如：制定环境质量标准和排放标准。如：制定环境质量标准和排放标准。 确定减少排放污染物的数量。确定减少排放污染物的数量。 对不同治理方案的经济性，有效性的分析。对不同治理方案的经济性，有效性的分析。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006

6、.4.18.B B、环境影响评价和环境质量预测中的应用环境影响评价和环境质量预测中的应用 如：新、改扩工程的环境影响评价（污染因如：新、改扩工程的环境影响评价（污染因子识别，污染程度预测模式等）。子识别，污染程度预测模式等）。 对城市结构变化、人口增长、污染物对城市结构变化、人口增长、污染物增长、能源结构改变和经济发展等造成环增长、能源结构改变和经济发展等造成环境质量变化进行预测。境质量变化进行预测。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.C C、污染物治理和给水、排水、水资源利用等污染物治理和给水、排水、水资源利用等方面应用方面应用 如：城市污水处

7、理流程优化、污染治理最佳如：城市污水处理流程优化、污染治理最佳运行控制、给排水管网系统的优化、多目标运行控制、给排水管网系统的优化、多目标水资源开发等。水资源开发等。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.系统模型化的基本理论（系统模型化的基本理论（4 4种）种）a a“黑箱黑箱”理论理论 对环境系统内部结构和行为不清楚，通对环境系统内部结构和行为不清楚，通过实验或现场实测方法得到系统输入过实验或现场实测方法得到系统输入输出输出规律，常用传递函数来描述系统。规律，常用传递函数来描述系统。b b“白箱白箱”理论理论 对环境系统内部结构和行为已掌握清对环

8、境系统内部结构和行为已掌握清楚，应用已知的科学知识进行描述。楚，应用已知的科学知识进行描述。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.c“灰箱灰箱”理论理论 介于上述两者之间，对系统内部结构和介于上述两者之间，对系统内部结构和行为的主要部分清楚，其他部分不清楚，是行为的主要部分清楚，其他部分不清楚，是“白箱白箱”与与“黑箱黑箱”理论相结合的一种方法。理论相结合的一种方法。d数理统计和分析的理论数理统计和分析的理论 对属于对属于“黑箱黑箱”但又不能进行实验的但又不能进行实验的系统，采用数理统计和分析的方法建立系统系统，采用数理统计和分析的方法建立系统模

9、型的理论和方法。在环境系统预测模型中模型的理论和方法。在环境系统预测模型中应用居多。应用居多。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.4、数学模型分类、数学模型分类 a按模型变量的随机性分：按模型变量的随机性分：n确定性模型：输入、输出均确定确定性模型：输入、输出均确定n非确定性模型：输入是随机的，解是不非确定性模型：输入是随机的，解是不稳定的，不具有唯一性。稳定的，不具有唯一性。b按随时间变化规律分：按随时间变化规律分：n稳态模型：系统内物质量不随时间而变稳态模型：系统内物质量不随时间而变n动态模型：系统内的物持量随时间而变。动态模型：系统内的物

10、持量随时间而变。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. c. c. 按空间维数分：按空间维数分：n一维模型，仅一个方向上有梯度。一维模型，仅一个方向上有梯度。n二维模型，二个方向上有梯度。二维模型，二个方向上有梯度。n三维模型，三个方向上有梯度。三维模型，三个方向上有梯度。 d d按物质的输移特性分：按物质的输移特性分：n平流模型，可忽略扩散项时。平流模型，可忽略扩散项时。n扩散模型，可忽略平流项时。扩散模型，可忽略平流项时。n平流扩散模型，两项均不可忽略时。平流扩散模型，两项均不可忽略时。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分

11、析课件2006.4.18.e按反应动力学的性质可分为：按反应动力学的性质可分为：n纯输移模型纯输移模型n纯反应模型纯反应模型n输移及反应模型输移及反应模型n生态模型生态模型 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.5. 数学模型的建立步骤数学模型的建立步骤（5步）观测数观测数据组据组模型结模型结构选择构选择模型模型应用应用观测数据组观测数据组参数参数估计估计检 验 与 验检 验 与 验证证图图2-2南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.参数估计：如弥散系数参数估计：如弥散系数E E、生物降解系数生物降解系

12、数k k1 1、大气复氧系数大气复氧系数k k2 2等的估值。等的估值。参数估计数值准确与否，关系到模型的实参数估计数值准确与否，关系到模型的实际应用中能否正确模拟实际情况，是项极际应用中能否正确模拟实际情况，是项极其重要的工作。其重要的工作。实际上，水质模拟过程，也就是对系统模实际上，水质模拟过程，也就是对系统模型进行识别，对模型参数进行估计，再用型进行识别，对模型参数进行估计，再用实际观测值进行检验、调整的反复试验过实际观测值进行检验、调整的反复试验过程。程。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.模型检验：检验与实际情况的吻合程度。模型检验：检

13、验与实际情况的吻合程度。 原因：建模时作过一些假定。原始数原因：建模时作过一些假定。原始数据误差可能使参数估计产生误差。据误差可能使参数估计产生误差。灵敏度分析：模型参数变动时造成的影响。灵敏度分析：模型参数变动时造成的影响。首先变动一个参数，其余参数保持不变，首先变动一个参数，其余参数保持不变，然后检查目标函数的变化程度，若变化不然后检查目标函数的变化程度，若变化不大，说明目标函数对这个参数不敏感，对大，说明目标函数对这个参数不敏感，对这个参数的估计可不要求很准确，若特别这个参数的估计可不要求很准确，若特别不敏感，说明这个参数在该模型中是多余不敏感，说明这个参数在该模型中是多余的，可剔除。的

14、，可剔除。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.二、环境系统现象模型推导基本原理二、环境系统现象模型推导基本原理 所谓现象模型我们指的是环境现象的所谓现象模型我们指的是环境现象的物理模型（客观），而不是经验模型物理模型（客观），而不是经验模型（非物理模型）。为了建立环境系统的（非物理模型）。为了建立环境系统的现象模型，需要以下两方面的知识：现象模型，需要以下两方面的知识：作为输送污染物的介质（如水、空气等）作为输送污染物的介质（如水、空气等）的流体的流动状态的描述。的流体的流动状态的描述。在污染物被输送过程中，污染物所发生在污染物被输送过程中，污染

15、物所发生的质与量的变化的描述。的质与量的变化的描述。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. 利用以上两方面信息，根据物质和能利用以上两方面信息，根据物质和能量的守恒原理，对某一污染物在流体介质中量的守恒原理，对某一污染物在流体介质中的浓度、流体的质量、动量或热量进行衡算，的浓度、流体的质量、动量或热量进行衡算，即用下列基本法则建立数学表达式：即用下列基本法则建立数学表达式： 积蓄量积蓄量= =（流入量流出量）（流入量流出量）+ +（产生量消灭量）（产生量消灭量） （式（式2-12-1）南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课

16、件2006.4.18. 在对象系统的介质流体中，取一均匀的在对象系统的介质流体中，取一均匀的微小单元微小单元V，并设系统的状态变量为并设系统的状态变量为（可以是流体的动量、流体的质量、流体可以是流体的动量、流体的质量、流体中污染物浓度、流体温度、热量等）由（中污染物浓度、流体温度、热量等）由（2-1）可得：可得： V=（AiJiA0J0）t+（BD）Vt (式式2-2)南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.式中：式中： t时间间隔中单位体积中时间间隔中单位体积中的增量的增量 （）。）。 Ji 、 J0 分别表示单位时间、单位截面上分别表示单位时间、

17、单位截面上 流入、流出的流入、流出的通量向量（通量向量（L3/L2T）。）。 Ai、A0分别表示流入和流出的截面积。分别表示流入和流出的截面积。 （L2） B、D分别是单位时间、单位体积中分别是单位时间、单位体积中产产生和消灭速率（生和消灭速率（/T）南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.n式（式（2-3）中状态变量）中状态变量赋予不同意义时可赋予不同意义时可得下列四种很有用处的平衡方程式：得下列四种很有用处的平衡方程式： 的意义的意义 所得的平衡方程式所得的平衡方

18、程式流体的动量流体的动量 运动方程式运动方程式 v x、y、z、t流体的质量流体的质量 连续性方程式连续性方程式 x、y、z、t流体中污染物浓度流体中污染物浓度 扩散方程式扩散方程式 c x、y、z、t流体温度或热量流体温度或热量 热扩散方程式热扩散方程式 T x、y、z、t南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.（场：速度场、密度场、浓度场、温度场等）（场：速度场、密度场、浓度场、温度场等） 以上方程式是我们研究污染物在大气和以上方程式是我们研究污染物在大气和水体中迁移转化规律时的基本方程式。水体中迁移转化规律时的基本方程式。南昌航空工业学院环境系

19、统分析课件南昌航空工业学院环境系统分析课件2006.4.18.三、环境系统现象模型的建立（三、环境系统现象模型的建立（B B、D D、的方程）的方程） 在研究环境系统中污染物的时空分在研究环境系统中污染物的时空分布时，除上述的一些方程式以外，还必布时，除上述的一些方程式以外，还必须收集到有关污染源产生强度以及污染须收集到有关污染源产生强度以及污染物在输送过程中所发生变化（如生物化物在输送过程中所发生变化（如生物化学反应等）的速率等信息，为此，就需学反应等）的速率等信息，为此，就需要建立要建立污染源模型污染源模型（B B）和和反应速度方程反应速度方程（D D），），以及描述环境系统的状态变量以及

20、描述环境系统的状态变量（）和参数（如温度、密度、压力、）和参数（如温度、密度、压力、等）之间关系的等）之间关系的“状态方程状态方程”。南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.建立这些方程时的基础信息关系可参见下图：建立这些方程时的基础信息关系可参见下图： 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. 此处此处“状态方程状态方程”就是把作为介质的就是把作为介质的流体温度、密度、压力等的热力学关系表流体温度、密度、压力等的热力学关系表示出来，以对温度在垂直方向的不同来控示出来，以对温度在垂直方向的不同来控制密度成层

21、的形成。制密度成层的形成。 “ “变化速度式变化速度式”是把水中的有机污染物是把水中的有机污染物质的分解和沉积，大气中的光化学反应等质的分解和沉积，大气中的光化学反应等由于化学、物理和生物学的各种原因而形由于化学、物理和生物学的各种原因而形成的发生和消失速度表示出来。成的发生和消失速度表示出来。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. 在适当的边界条件和初始条件下求解式在适当的边界条件和初始条件下求解式（2-3）便可得）便可得在时间上和空间上的分布信在时间上和空间上的分布信息，但是对于和息，但是对于和V、D、B，一起求解时往往一起求解时往往很困难，

22、既使不考虑很困难，既使不考虑B、D只考虑只考虑V时也不容时也不容易求解，因此，常必须考虑把流动特性易求解，因此，常必须考虑把流动特性V及及变化速度式变化速度式D、B按照实用的目的来进行简按照实用的目的来进行简化。化。 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.（1）流动特性的简化，）流动特性的简化，a弥散扩散方程式弥散扩散方程式 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18. 南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.南昌航空工业学院环境系统分析课件南昌航空工业学院环境系统分析课件2006.4.18.b湍流扩散方程的降维湍流扩散方程的降维 经简化后的式（经简化后的式（2-4）仍是一个关于）仍是一个关于x、y、z、t四个独立变量的偏微分方程，四个独立变量的偏微分方程，求解需知道求解需知道v和和E在在x、y、z方向上的分量，方向上的分量，取得这些数据是很困难的，所以在实际取得这些数据是很困难的，所以在实际