河流水质模型

河流水质模型第1页，课件共75页，创作于2023年2月例如排出水与河水之间因流速分布和紊动作用而产生的质量交换；排出水与河水的温差所产生的热交换；排出水与河水的密度差而产生的浮力作用；排出水与河水的动量交换等。对于浮力为中性的非射流排放情形，竖向混合区的长度与水深成正比，大致为排放处水深的几十倍到一百倍，其距离较短。对于浮力出流和射流问题，也有许多估算其混合区范围的公式。李光炽水质模型第2页，课件共75页，创作于2023年2月(2)从竖向充分混合起至河流横向开始充分混合为止。天然河流的河床一般是宽浅型的，宽深比大于10。达到横向混合所需要的河段长度比达到竖向混合所需要的河段长度大得多，河越宽则所需距离就越大，可达几公里、几十公里，对于大河甚至达上百公里。

李光炽水质模型第3页，课件共75页，创作于2023年2月(3)从横断面上开始充分混合以后的阶段。在这个阶段，河流断面上各点水质浓度的偏差远比各横断面间的断面平均浓度偏差小。因此，一般只需考虑断面平均浓度沿河流纵向的变化就可以了。从排放口至第三混合阶段开始之间的距离L可按下式估算：

李光炽水质模型第4页，课件共75页，创作于2023年2月排放口混合距离L估算式：

式中B—河流平均宽度；u—河流平均流速；u*—摩阻流速；g—重力加速度；I—河流水力坡降；H—

平均水深。

李光炽水质模型第5页，课件共75页，创作于2023年2月根据河流中的各混合阶段，可采用不同的河流水质模型。例如当研究河流水质规划时，河段长度比横向尺度大得多，可以主要考虑第三混合阶段，采用一维模型或零维模型。而当研究污染物排放口与取水口相互位置的关系，而且河面相对较宽阔时，宜应采用二维模型或三维模型。李光炽水质模型第6页，课件共75页，创作于2023年2月5.2

均匀混合水质模型(零维水质模型)

均匀混合水质模型是把一个水体，如一个河段，看作是一个均匀混合的反应器。假定入流进入反应器后立即均匀分散，各水团完全均匀混合。其质量平衡方程为李光炽水质模型式中V—反应器内水的体积；Q—反应器入流及出流流量；C0、C－入流及反应器内污染物浓度；r(C)－反应器内过程的反应速率，与C有关；S－除含r(C)项以外的源和汇强度。第7页，课件共75页，创作于2023年2月如果S＝0，可简化为当所研究的水质组分在反应器内的变化过程符合一级反应动力学规律时，可表示为

式中

―

一级反应速率。

李光炽水质模型第8页，课件共75页，创作于2023年2月对于稳态情形，，则有

由此可解得C为

或

其中t=V／Q，称滞留时间。李光炽水质模型第9页，课件共75页，创作于2023年2月李光炽水质模型第10页，课件共75页，创作于2023年2月将零维模型应用于实际河流的稳态水质模型或预测时，首先需要将河流分成若干河段，每一河段再划分为长度为的若干微段，每一个微段即视为一个完全混合的反应器

式中C0

—第i河段起始断面污染物浓度；m—第i河段内的微段数。李光炽水质模型第11页，课件共75页，创作于2023年2月每一河段起始断面的污染物浓度C0，可按质量平衡方法求得式中—上游流入i河段的流量；

—河段起始端旁侧入流(支流或排污口)的流量；

—上游来水的污染物浓度；

—

旁侧入流的污染物浓度。均匀混合模型适用于均匀河段，要求足够小，否则会造成较大误差。

李光炽水质模型第12页，课件共75页，创作于2023年2月5.3

一维BOD-DO水质模型BOD-DO模型的基本假定是：(1)BOD的降解符合一级动力学反应规律；即在任何时候反应速率都和剩余的有机物数量成正比。以L表示BOD浓度，则。(2)水体中溶解氧DO的减少只是由于BOD降解所引起的，而且与BOD的降解有相同的速率。(3)水体中复氧的速率与氧亏量成正比。氧亏是指溶解氧浓度O与可能达到的饱和溶解氧浓度Os之差。李光炽水质模型第13页，课件共75页，创作于2023年2月一维BOD-DO水质方程(Streeter-Phelps模型)式中L—BOD浓度；O—溶解氧浓度，ux

—断面平均流速，DL

—纵向弥散系数；K1

—耗氧系数：K2

—复氧系数；Os

—

饱和溶解氧浓度。李光炽水质模型第14页，课件共75页，创作于2023年2月对于稳态情形

李光炽水质模型边界条件

第15页，课件共75页，创作于2023年2月解为

李光炽水质模型第16页，课件共75页，创作于2023年2月当忽略弥散项时有如下形式的解：

李光炽水质模型表示氧亏浓度

第17页，课件共75页，创作于2023年2月

BOD和DO的沿程变化李光炽水质模型第18页，课件共75页，创作于2023年2月溶解氧的沿程变化曲线表明，河段内可能出现量低溶解氧浓度值，相应的氧亏值称临界氧亏Dc，相应的距离称为临界距离xc。

李光炽水质模型式中f为水体的自净系数。

第19页，课件共75页，创作于2023年2月在实际应用中，Streeter—Phelps基本模型暴露出它的某些局限性，因面出现了许多修正模型，主要有以下三种。(一)BOD沉淀的修正式(Thomas修正式)当废、污水很少处理就排放时，在排放口下游经常可以看到某些物质的沉淀，从面导致了BOD的减少。为此，可在BOD方程中引入沉淀系数K3，将K3与K1相加；但在DO方程中仍然保留一个Kl，这是因为BOD的这一部分减少并不是降解所致，与DO的减少无关。修正的模型可表示为李光炽水质模型第20页，课件共75页，创作于2023年2月以上修正式也称Thomas修正式。其中K3为负值时表示BOD沉淀物的再悬浮。K3的量值一般在－0.36～0.36d－1之间变化。李光炽水质模型第21页，课件共75页，创作于2023年2月当边界条件为，，时，其解为李光炽水质模型第22页，课件共75页，创作于2023年2月(二)地表汇流、底泥耗氧，光合作用增氧修正式(Dobbins-Camp修正式)Dobbins和Camp在Streeter—Phelps方程的基础上考虑了如下的修正：①由于地表面源污染的汇入或(和)底泥中有机物重新悬浮，引起河流BOD的变化，其变化速率以常数R表示，附加在B0D方程右侧。②由于底泥有机物降解耗氧，水生植物光合作用增氧和呼吸耗氧等综合作用，引起溶解氧浓度发生变化，其变化速率以常数P表示，附加在DO方程的右侧，修正式如下所示：李光炽水质模型第23页，课件共75页，创作于2023年2月李光炽水质模型当边界条件为，，时，其解为第24页，课件共75页，创作于2023年2月李光炽水质模型第25页，课件共75页，创作于2023年2月上述式称Dobbins-Camp修正式。如果R＝O，P=O，则成为Thomas修正式；如果只R=0，P=O，K3=O，则成为Streeter-Phelps基本模型。

同样可求解溶解氧的临界距离xc：

李光炽水质模型当给定临界氧亏Dc的值时，可以确定BOD的允许排放浓度L0：第26页，课件共75页，创作于2023年2月（三）硝化阶段耗氧的修正式(O'Connor修正式)

一般有机物的氧化过程包括碳化和硝化两个阶段。O‘Connor假定总的BOD是碳化需氧和硝化需氧量这两部分之和，即，相应阶段的反应速率分别记为和，则修正模型可表示为李光炽水质模型第27页，课件共75页，创作于2023年2月李光炽水质模型模型的解为

第28页，课件共75页，创作于2023年2月硝化需氧量LN0可用水体中可氧化的氮的总量来估计。一般污水中可氧化的氮主要以氨氮形式出现，当它被氧化时，一个氮原子完全氧化需两个氧分子(4个氧原子)，所以，当氨氮浓度给定时，LN0为李光炽水质模型第29页，课件共75页，创作于2023年2月5.4

一维河流水质模拟对于一条实际河流应用水质模型进行水质模拟和水质预测时，首先应对河流水量水质的基本数据和资料进行如下分析和整理。(1)河流的径流量。河流的径流量对稀释作用和自净能力有重要影响。天然河流的径流量一般变化幅度较大，而污水流量变化幅度小，因而在河道的枯水季节，河流的水质状况较差。李光炽水质模型第30页，课件共75页，创作于2023年2月进行水质模拟，应根据摸拟计算的目的，确定河流计算流量。对于河流水污染控制规划，可选用十年最枯月平均流量作为计算流量。(2)河流流速、断面平均水深、水面宽度、河床纵坡、糙率等。(3)河流支流的上述资料及汇入主流的位置。(4)沿河流的污染源分布，污水流量，主要污染物。(5)河流水质，如COD、BOD5、DO等的监测资料。李光炽水质模型第31页，课件共75页，创作于2023年2月对河流分段。分段的原则是使得每一河段大体适用一维河流水质模型应用的条件，即主要使得平均流速和有关参数保持恒定。根据河流自然特征与沿程流量输入输出状况，分为n个河段，共含n+1个断面，依次由上游向下游编号为0，1，…，i，i+1，…，n，其中编号0代表上游输入断面，n代表下游输出断面。

李光炽水质模型第32页，课件共75页，创作于2023年2月一维河流水质模拟第33页，课件共75页，创作于2023年2月Qi

——在断面i处注入河流的污水流量；Q1i

——由上游流到断面i的河流流量；Q2i

——由断面i向下游流出的河流流量；Q3i

——在断面i处引走的流量；Li、Oi

——在断面i处注入河流的污水BOD5浓度与DO浓度；李光炽水质模型第34页，课件共75页，创作于2023年2月L1i、O1i

——由上游流到断面i的河水BOD5浓度与DO浓度；L2i、O2i

——由i断面向下游输出的河水BOD5浓度与DO浓度；K1i、K2i、K3i

——分别是由i-1断面至i断面间的BOD5衰减速率常数，复氧速率常数与沉淀或再悬浮速率常数；ti——

河水由i-1断面至i断面的流动时间。

李光炽水质模型第35页，课件共75页，创作于2023年2月水量关系：

水质关系：

李光炽水质模型第36页，课件共75页，创作于2023年2月i-1断面至i断面之间的BOD5关系与复氧关系：

若令

李光炽水质模型第37页，课件共75页，创作于2023年2月可以简写成:令李光炽水质模型第38页，课件共75页，创作于2023年2月可以写出任—个断面处BOD5值与上游各断面及汇入本断面的污水BOD5值之间的递推关系式李光炽水质模型............第39页，课件共75页，创作于2023年2月用一个矩阵方程来表示：李光炽水质模型其中第40页，课件共75页，创作于2023年2月A和B是两个n×n阶矩阵，A非奇异，A的逆矩阵A－1存在，则可以推导出：李光炽水质模型是给出i＝1断面处BOD5值的n维列向量，其中第41页，课件共75页，创作于2023年2月矩阵方程给出了河流每一个断面向下游输出BOD5值()与各个断面输入河流的BOD5值()之间的关系，在水质模拟时，是一组已知量，是需要求解的量。而在进行水污染控制系统规划时，是一组已知的河流BOD5的约束量，则是需要确定的量。

李光炽水质模型第42页，课件共75页，创作于2023年2月DO的模拟李光炽水质模型若令得第43页，课件共75页，创作于2023年2月以矩阵方程形式表示如下：李光炽水质模型其中第44页，课件共75页，创作于2023年2月C和D是两个n×n阶矩阵，C非奇异，其逆矩阵存在，所以李光炽水质模型是由河流各断面往下游输出的DO值组成的n维列向量。是各断面输入河流的污水的DO浓度组成的n维列向量，通常这是一组已知的量。第45页，课件共75页，创作于2023年2月李光炽水质模型及都是表征起始条件影响的n维列向量其中将的表达式代入得

第46页，课件共75页，创作于2023年2月令

则可写成

李光炽水质模型第47页，课件共75页，创作于2023年2月U和V是两个由给定数据计算的n×n阶下三角矩阵，和是两个由给定数据计算的n维向量。每输入一组BOD5()值，就可以获得一组相应的河流BOD5值和DO值(和)，由于U和V反映了这种输入、输出的因果关系，故称U为河流BOD稳态响应矩阵，称V为河流DO稳态响应矩阵。李光炽水质模型第48页，课件共75页，创作于2023年2月5.5QUAL—Ⅱ河流水质综合模型一、简介QUAL—Ⅱ是一个具有多种用途的河流水质模型，它能按照使用者的要求，以各种组合方式描述以下十三种水质参数(或称为水质变量)：(1)溶解氧(DO)；(2)生化需氧量(BOD)；(3)水温(T)；(4)叶绿素—α——藻类；(5)氨氮；(6)亚硝酸氮；(7)硝酸氮；(8)可溶性磷；(9)大肠杆菌；(10)任选的一种可降解物质；(11)三种任选的不降解物质。

李光炽水质模型第49页，课件共75页，创作于2023年2月该模型假设：在河流里物质的主要迁移方式是对流和弥散，而且认为这种迁移只发生在河流或水道的纵轴方向上，所以该模型是一个一维模型。这个模型可以描述同时有多个排放口并有支流流入和取水流出的河流系统。在水力学方面，QUAL—Ⅱ只限于描述流量不随时间变化时的水质情况。在其它方面QUAL—Ⅱ既可作为稳态模型也可作为动态模型，所谓动态模型是反映水质随气象条件和排放量的日变化而变化的动态过程的模型。李光炽水质模型第50页，课件共75页，创作于2023年2月在水力学参数随时间的变化是在足够缓慢的情况下仍可应用QUAL—Ⅱ模型，这时只要把时间分成一定的间隔并认为在一个时间间隔内水力学参数是常数即可。QUAL—Ⅱ可用于水质变化对河流的影响。该模型能够研究由于藻类生长和呼吸过程引起的溶解氧日变化，还能被用于研究污染物的瞬时排放对水质的影响，如有关污染源的事故性排放、季节性排放或周期性排放对水质的影响。李光炽水质模型第51页，课件共75页，创作于2023年2月二、各水质变量之间的相互关系

李光炽水质模型大气溶解氧BOD氨氮亚硝酸氮叶绿素α（藻类）硝酸氮正磷酸盐11121684271514511310936第52页，课件共75页，创作于2023年2月1．复氧作用；2．河底生物(包括底泥)的耗氧；3．碳化合物BOD耗氧；4．光合作用产氧；5．氨氮氧化耗氧；6．亚硝酸氮氧化耗氧；7．碳化合物BOD的沉淀；8．浮游植物对硝酸氮的吸收；李光炽水质模型第53页，课件共75页，创作于2023年2月9．浮游植物对磷(磷酸盐磷)的吸收；10．浮游植物呼吸产生磷(磷酸盐)；11．浮游植物的死亡和沉淀；12．浮游植物呼吸产生氨氮；13．底泥释放氨氮；14．氨氮转化为亚硝酸氮；15．亚硝酸氮转化为硝酸氮；16．底泥释放磷。

李光炽水质模型第54页，课件共75页，创作于2023年2月三、河流系统的概化

QUAL—Ⅱ模型把河流系统表达为一系列河段构成的网络，用节点把这些河段联系在一起，同时假定在同一河段里水力学参数保持不变。每一河段又被分成许多小节(小段)。把所有的节分为七种类型：1．源头节(主流和支流的第一节)；2．正常节；3．支流入口的上游节；4．支流入口节；5．河系的末节；6．含有点源的节(点源指除支流外的污染源)；7．有出流的节(如取水灌溉等等)。李光炽水质模型第55页，课件共75页，创作于2023年2月李光炽水质模型第56页，课件共75页，创作于2023年2月四、模型方程

QUAL—Ⅱ的基本方程是一个对流—弥散质量迁移方程，它能描述任一水质变量的时间与空间变化情况。在方程里除平移和弥散项外还包括由化学、物理和生物作用引起的源漏项(包括支流和排放口的影响)。对于任意的水质变量C，这个方程可以写成如下形式：

李光炽水质模型第57页，课件共75页，创作于2023年2月QUAL—Ⅱ模型假设流量处于稳定状态，即=0，=0，于是变成如下形式：

式中E——河流纵向弥散系数，m2／s；Sint——水质变量C的内部的源和汇(如化学反应等)，kg／(s.m3)；Sext——外部的源和漏(如支流的影响等)，kg／(s.m3)。李光炽水质模型第58页，课件共75页，创作于2023年2月五、模型的水力学部分

QUAL—Ⅱ模型假设水力学系统处于稳定状态，即=0，从计算单元的水量平衡关系可得：

式中（qx）i——第i个计算单元的外部入流与出流的总和。当已知流量Q之后，流速u和水深H就可用经验公式得到：，。可以用经验公式的方法来确定纵向离散系数。李光炽水质模型第59页，课件共75页，创作于2023年2月六、模型的源汇项

（1）叶绿素—α(浮游藻类)叶绿素—α的浓度与浮游藻类的生物物质量的浓度成正比，为了建立关于叶绿素—α的模型，用下面的简单关系将藻类的生物物质转换为叶绿素—α的量：式中Cca——叶绿素—α的浓度；CA——藻类生物质量的浓度；α0——转换系数。李光炽水质模型第60页，课件共75页，创作于2023年2月描述藻类(叶绿素—α)生长与产量的微分方程，可由下面的关系得到：式中μ——藻类比生长率，它随温度而变化；ρA——藻类呼吸速率常数，随温度变化；σ1——藻类沉淀速率常数；H——平均水深。李光炽水质模型第61页，课件共75页，创作于2023年2月(2)氮的循环在QUAL—Ⅱ模型里考虑了三种形态的氮：氨氮(CN1)、亚硝酸氮(CN2)和硝酸氮(CN3)。氨氮：

式中α1——藻类生物量中氨氮的比例；σ3——水底生物的氨氮释放速率；A——平均横截面积；KN1——氨氮氧化速率常数；李光炽水质模型第62页，课件共75页，创作于2023年2月亚硝酸氮：式中KN2——亚硝酸氮氧化的速率常数；其它符号意义同前。硝酸氮：式中符号意义同前。李光炽水质模型第63页，课件共75页，创作于2023年2月(3)磷循环在QUAL—Ⅱ里关于磷循环的模型不象氮循环模型那样复杂，该模型只考虑了溶解性磷和藻类的相互关系，以及底泥释放磷的项，模型方程如下：式中Cph——磷酸盐(换算成磷)的浓度；α2——在藻类生物质量中磷所占的比例；σ2——底泥释放磷的速率；

李光炽水质模型第64页，课件共75页，创作于2023年2月(4)碳化DOD碳化BOD的变化速率按—级反应来考虑，可得到下面微分方程：式中K1——碳化BOD的降解速率常数，与温度有关；K3——由于沉淀作用而引起的碳化BOD消耗速率常数。李光炽水质模型第65页，课件共75页，创作于2023年2月(5)溶解氧(DO)在QUAL—Ⅱ模型中描述溶解氧变化速度的微分方程形式如下：

式中C——DO浓度；CS——DO的饱和浓度；李光炽水质模型第66页，课件共75页，创作于2023年2月α3——单位藻类的光合作用产氧率；α4——单位藻类的呼吸作用的耗氧率；α5——单位氨氮氧化时的耗氧率；α6——单位亚硝酸氮氧化时的耗氧率；K2——复氧系数；K4——水底耗氧常数。其它符号意义同前。李光炽水质模型第67页，课件共75页，创作于2023年2月(6)大肠杆菌大肠杆