2023年可持续水资源管理P65P85.doc

可持续水资源管理—理论 方法 应用〔P65—P85〕 2．2 水质模型 建立反映水环境污染客观规律的数学模型是对是对水环境进行预测评价、规划和科学管理的根底，也是建立社会资源管理的根底，也是建立可持续水资源管理量化模型的根底模型之一，本文首先简介水环境的模型化过程，然后重点描述各类水体的水质模型。 2．2．1水环境模型化过程 所谓水环境的模型化，就是把与水环境有关的变量之间的关系及其影响规律抽象成一系列反映水环境性能和机理的数学方程式形式，并利用数学手段对其进行求解，已解决实际的水环境问题。它的建立是在对水环境进行反复地考察研究，通过实验或实地监测取得大里的有关信息和数据，进而对所研究的系统行为和变化规律有了深刻认识的根底上，经过简化和数学演绎而得出的一系列数学表达式.因此，建立水环境的数学棋型一直是从事水资源和水环境研究的一项重要根底工作。 水环境的数学模型，主要是指各类水体的水质模型、水质评价模型和水环境规划模型等。一般而言，建立的水环境模型应具有如下特性。 ①现实性即所建一立的模型在一定程度上要能够反映和符合系统客观实体的情况和规律，应把系统本质的东西和系统内外部不阅组分之间的主要关系反映出来，把非本质或次要的东西省略，而又不影响对现实模拟的真实程度。 ②准确性 由于任何一个水环境系统都是作常复杂的，其自然规律目前人类还没有完全掌握，因此在实际建模过程中不得不忽一些难以用数学公式表达的关系，但这种忽略要建立在满足解决问题所必需的精度要求之上。在确认模型参数的过程中，也要致力于保证参数的真实性和可靠性，从而保证模型的准确性。 ③简单性 在保证模型的现实性和准确性的根底上，要尽可能地使所建立的模型简单化，以节省建模和运用的费用和时间，并增强其可应用性。如果一个简单的模型能使实际问题得以满意的解决，那它就优先于一个复杂的模型。 在水环境评价和规划中!模型化也起着十分重要的作用。其原因在于，在没有或难以建立实体系统以供观测和试验的条件下，为了求得系统的参数和确定各种制约条件下系统的变化趋势，模型方法确实是一种有效的科学手段。目前!水环境评价和规划工作，几乎都是借助模型来完成的。一般做法是利用建立的模型!给出不同的输入参数和根底数据，观察系统的输出，然后根据最后的结果评价环境影响，并选择最正确的系统设计方案。 2.2.1.2 水环境模型的分类 水环境的数学模型可从不同的方面进行分类。 按照水环境模型研究对象的水体特征，可以分为河流河口水质模型、湖泊水库水质模型、海洋水质模型和地下水水质模型等。河流模型的开发已经比拟成熟，能够较好地反映实际情况。湖、海模型都比拟复杂，可靠性较差。 按照水环境模型中的水质组分，可以分为单组分、耦合组分、多组分和水生生态模型等，其BOD-DO耦合模型是应用最多的一种水质模型。 按照水环境数学模型与时间的关系，可以分为稳态模型和非稳态模型两类。稳态是指污染物的浓度不随时间变化，然而天然河流实际上并不存在这种理想状态，但我们可以把河流某一段时间内相对稳定的流量和排污情况作为稳态来处理。动态是指污染物的浓度随时间变化而变化，它可用来估算暴雨径流和污染事故泄漏等瞬时变化的情况。 按照模型的空间维数，可以分为零维模型、一维模型、二维模型和多维模型。在实际应用中!选用几维模型应根据研究区域的空间范围及其水体中污染物的混合情况而定。例如，如果仅是对区域性水质进行粗略的模拟计算!采用零维模型即可；如果对一个较长的河段或河流进行水质模拟计算，可以采用一维模型；如果要研究河流局部河段范围内的水质状况，排污口或入海口附近的污染物分布以及污染带等类似间题时，就要考虑采用二维模型或三维模型。 按照所建模型的物质的反响特性，水环境模型可以分为保守性物质模型和非保守性物质模型。前者指像重金属、农药等类的物质，其在水环境中不(或不易)发生生物氧化分解；后者指可〔易〕降解的一类有机污染物。 2.2.1.2水环境系统建模过程 建立水环境模型的过程一般可分为如下几个步骤。 ①分析模型的使用目的和要求，确定模型的功能。在建立模型之前，首先必须明确以下几个问题：为什么样的对象建立模型需要建立什么样的模型？如何保证使建立的模型符合使用要求 ②进行科学的系统分析，并根据目的和要求从时间和空间等角度确定研究系统的边界条件，确认系统的输入、输出和状态等变量及其性质。 ③在对系统的变化机制有了初步分析、推断的根底上，对系统进行调研、监测，取得建模所必需的数据和资料。 ④识别模型的结构，在理解系统各要素之间的作用关系的前提下，建立起相应的函数关系。 ⑤根据监测得到的数据，使用适当的数学技术估算出模型中的各个参数值。 ⑥将几组在参数估值时未使用过的数据代入模型，并计算模型的精度，了解是否符合模型的精度要求和使用要求。假设符合要求!那么建模过程结束。假设不符合要求，那么重新确定模型的结构，重复以上④～⑥过程。 2.2.2 河流水质模型 地表水的水质模型可以分为河流水质模型/湖泊和水库水质模型/海洋水质模型/非点源水质模型等,其中每种类型的模型又都可按照不同的特点分为各种子类别,例如河流水质模型又可分为零维、一维、二维、多维和BOD-DO耦合模型等。下面就逐一进行介绍。 针对水体中不同的混合特点，可以从空间维数来区分以下几种根本方程。 图2-12 连续完全混合反响器 2.2.2.1零维水质模型 当把一个水体〔如一个河段〕看作是一个完全混合反响器〔见图2-12〕并假设水流进入该系统后即完全分散到整个系统，其中各水团是完全混合均匀的。对这种连续流完全混合反响器的理想情况，可列出以下的质量平衡关系： 积累率=进入量-输出量+衰减量 即 式中，V为反响器内水的体积,L3；C0，C分别为入流与反响器内污染物的浓度，ML-3；S为反响器内其他来源和漏失的污染物量总和，MT-1；r〔C〕为反响器内过程的反响速率，ML-3T-1； Q为流经反响器的流量。 如果系统内S=0，那么上式可变为 如果研究的污染物在反响器内的反响符合一级动力学反响，那么r〔C〕=-K1C，式〔2-57〕可进一步改写为 式中，K1为一级反响速率常数，T-1。 式〔2-59〕就是常用的零维水质模型根本方程.当假设河段中没有源和漏失量〔S=0〕，污染物反响属于一级动力学反响，且河流处于稳定条件下，即dC/dt=0，式〔2-59〕转化为： 对于任一具有以上条件的河段，可以将其分割为i个子河段〔见图2-13〕.对于其中第i个河段，又可以分为m个单元〔1,2，…，m〕。对于每个单元，都可以看作是一个完全混合反响器。因此，可以把这m个单元看作m个串联的完全混合反响器。并由此可以得到以下各式： 图2-13 河流零维模型概念 起始断面浓度为 式中，C0为第i河段起始断面浓度，[ML-3]；Cj为第i河段第j子段的污染物浓度，[ML-3]；U为河段平均流速，[LT-1]；△x为单元长度，[L]；Qi-1为上游流入第i河段的流量，[ML-1]；qx-1为第i河段起始端旁侧入流流量，[ML-1]；Ci-1为上游来水的污染物浓度，[ML-3]；Cx-1为第i河段起始端旁侧入流污染物浓度，[ML-3]。 利用上式可以计算出各河段中每个微小单元的浓度Cj。零维河流水质模型适用于均匀河段，并且要求取的△x足够小，否那么误差较大。 2.2.2.2 一维水质模型 〔1〕水量平衡方程 考察河流中长度为△x的微小河段的水量平衡关系〔见图2-14〕。图中：Q(x)，Q(x+△x)为流入和流出量，[L-3T-1]；q、qb分别为侧向输入和底部渗出流量，[L-3T-1 L-1]；pR、pE分别为单位水面的降水量与蒸发量，[L-3T-1 L-2]；A为河床断面，[L2]；B为河面宽，[L]。那么△x时间间隔内微小河段中的水质量平衡方程为： 式中△m为△t时间内微元河段中水质量的增量，[M]；ρ为水的密度，[ML-3]，ρ≈1. 用△t、 △x除方程两边： 图2-14微小河段的水量平衡关系 图2-15河流中的一个假想微元 当△t→0，△x→0时： 通常可以把pR、pE和qb忽略，那么产生了水平衡方程式： 〔2〕水质平衡方程 首先假定水流运动是推流，即所有水片以同一速度v运动〔见图2-15〕.那么污染物质量变化可以写成： 式中S1、S1’分别为单位时间、单位长度上侧向和底部的源和漏，[ML-1T-1]；Ss为单位时间、单位面积上的源和漏，[ML-2T-1]；Sv为单位时间、单位体积内的源和漏，[ML-3T-1]对式中的源取正值，漏取负值。 用△x、 △t除方程两边，并使△x、 △t趋于0，得 或 该式即为推流的水质平衡方程，其中AS= S1+ S1’+BSs+ASv (3)根本方程 当考虑分子扩散/湍流扩散和弥散过程时,推流流动水质平衡方程又可变化为以下3种形式。 ①增加分子扩散〔点的瞬时值〕用流场中水质点的瞬时值C和u来表达其浓度和流速，那么式〔2-71〕变为： ②增加湍流扩散〔点的时平均〕采用流场中点的平均值、表达，此时，C=+C’，u=+u’(C’和u’ 为因湍流扩散所引起的偏差量) 。把此转换关系和Q=uA代入式〔2-72〕，经整理可得： ③增加弥散〔断面平均值〕此时应采用河流断面的平均、来表达。其转换关系为C=+Cx，u=+ux(Cx和ux为弥散引起的偏差量)。代入式〔2-73〕并经推导整理可得到如下的一维河流水质微分方程： 或写成 此时，必须注意式〔2-75〕中的水流流速和水中污染物的浓度都是指断面的平均值!而非任意点的瞬时值或平均值，该式即为一维河流水质微分方程。 实际应用中，可以忽略分子扩散作用，因此一维河流水质根本方程即为： 对于均匀河段!断面面积且为常数，此时一维河流水质根本方程为： 对于有弥散现象的河流，一般x方向的弥散系数Dx比扩散系数Ex大得多〔Dx≈10～103，Ex≈10-2～1〕，因此，往往可忽略扩散项，那么式〔2-77〕可改写成一维均匀河流水质根本方程： 其中 如忽略河流底部渗流项S1’，并用平均水深H代替A/B，那么 对于河底无渗漏、忽略面源的侧向输入、污染物为一级衰减反响的一维均匀河流水质模型的根本方程为： 式中，L为t时存有的污染物浓度，[ML-3]；K1为反响速率常数，[T-1]。 〔4〕根本方程的解析解 一维均匀河流水质模型的根本方程的通式为： 按照河流水文和排污状况，先确定求解条件，求解条件分类如下： ①稳态解 稳态是指均匀河段、定常排污的条件，即污染物输入量W 、断面面积A、流速u和Dx恒定，均不随时间变化。因此，有，可得如下常微分方程： 给定边界条件：x=0，C=C0；x→∞，C=0.可用解特征多项式的方法对方程求解，得到稳态解： ②忽略弥散当河流的弥散系数很小时，一般可以忽略，即D=0，于是一维河流水质方程变为： 该偏微分方程可化为由下式给出的两个常微分方程： 式中，x〔t〕为该方程的特征线。对上式积分得到： 由上式可看出!只要知道初始浓度C=C0，即可求得水团流到下游某个地点的时间与浓度。 图2-16 瞬时排污的情况 ③瞬时突然排污条件时的解首先构造一个δ(t)函数〔见图2-16〕.设河段起始端〔t=0，x=0〕，突然瞬时排放质量为W 、流量为Q的平面污染源。其排污状况可以用一个以短时间△t为等比指标的δ(t)函数来描述。当等量的平面污染源W/D[ML-3T]，用不同的时间间隔〔如△T、△T/2、△T/4、…〕释放，那么排放的浓度可表