一维水量水质模型

第七章 一维非恒定河流和河网水量水质模型

对于中小型河流，通常其宽度及水深相对于长度数量较小，扩散质（污染物质、热量）很容易在垂向及横向上达到均匀混合，即扩散质浓度在断面上基本达到均匀状态。这种情况下，我们只需要知道扩散质在断面内的平均分配状况，就可以把握整个河道的扩散质空间分布特征，这是我们可以采用一维圣维南方程描述河流水动力特征或水量特征（水位、流量、槽蓄量等）；用一维纵向分散方程描述扩散质在时间及河流纵向上的变化状况。特别地，对于稳态水流，可以采用常规水动力学方法推算水位、断面平均流速的沿程变化；采用分段解析解法计算扩散质浓度沿纵向的变化特征。但是，在非稳态情况下（水流随时间变化或扩散质源强随时间变化）解析解法将无能为力（水流非恒定）或十分繁琐（水流稳态、源强非恒定），这时通常采用数值解法求解河道水量、水质的时间、空间分布。在模拟方法上，无论是单一河道还是由众多单一河道构成的河网，若采用空间一维手段求解，描述水流、水质空间分布规律的控制方程是相同的，只不过在具体求解方法上有所差异而已。

7.1 单一河道的控制方程 7.1.1 水量控制方程

采用一维圣维南方程组描述水流的运动，基本控制方程为：

?Q?Z?BW?q (1) ?x?tnuQ?Q?Q?Z?A?2u?gA?u2?g4/3?0 (2) ?t?x?x?xR式中t为时间坐标，x为空间坐标，Q为断面流量，Z为断面平均水位，u为断面

平均流速，n为河段的糙率，A为过流断面面积，BW为水面宽度（包括主流宽度及仅起调蓄作用的附加宽度），R为水力半径，q 为旁侧入流流量（单位河长上旁侧入流场）。此方程组属于二元一阶双曲型拟线性方程组，对于非恒定问题，现阶段尚无法直接求出其解析解，通常用有限差分法或其它数学离散方法求其数值解。在水流稳态、棱柱形河道条件下，上述控制方程组退化为水力学的谢才公式，可采用相应的方法求解水流特征。

7.1.2 扩散质输运控制方程

描述河道扩散物质运动及浓度变化规律的控制方程为：带源的一维对流分散（弥散）方程，形式如下：

2?(AC)?(QC)???c?A????AE?KAC?Sr?S (3) x???t?x?x??x?h式中，C为污染物质的断面平均浓度，Q为流量，Ex为纵向分散系数，S为单位时间内、单位河长上的污染物质排放量，K为污染物降解系数，Sr为河床底泥释

放污染物的速率。

1

此方程属于一元二阶偏微分方程，对于非恒定水流问题，微分方程位变系数的偏微分方程，现阶段尚无法直接求出其解析解，通常用有限差分法或其它数学离散方法求其数值解。在水流稳态、污染源源强恒定条件下，可按水动力特征将河道分为若干子段，在每个分段上，上述控制方程简化为常系数的常微分方程，可采用解析方法秋初起理论解。

7.2 单一河道一维水量水质模型 7.2.1 单一河道一维水量模型 （1）控制方程的离散

采用四点隐式差分格式离散方程组。如图1所示，河道被(n+1)个断面分为n个子河段，在第i个子河段M(i,i+1)上，对任一变量?取： ?(M)?(?ij??ij?1)/2 (4) ??(M)???i?1??i??xj?!j?1??(1??)???xji?1??ij? (5) ?1jj??(M)?ij?1??ij?1??i??i?1??t2?t (6) 第i个子河段2ii+1nn+1 图1 计算断面示意图 式中，上角标表示时间坐标，下脚标表示空间坐标。?为空间差商的权重系数(0???1),?=0时，此格式为显式格式，而当??0时，此格式具有隐式差分的特征。为使差分方程保持无条件稳定，必须??0.5。采用下式进行阻力项的线性化： j2??n2u?j???n2Qunuj?1j?1???Qi??4/3?Qi?1?0.5 (7) g4/3?g?4/3???R???RR?i???i?1??将式(4)-(6)代入连续方程得第i个子河段的差分方程：

CiZi?CiZi?1?Qi?Qi?1?Di (8)

式中，Ci?Bwi?1/2?xi?xi(1??)jjj(Qi?Qi?1)?Ci(Zij?Zi?)?q，Di? 1i2??t??下角标i+1/2表示断面i与断面i+1河段的均值。按照同样的方法，可得动量方

程的差分方程：

EiQi?GiQi?1?FiZi?FiZi?1?Hi (9)

2

式中，Ei??xi?2uij?1/22??tn2u?g?xi??? ?4/3??2??R?ijFi?gA?Bu2??i?1/2

jGi??xi?2uij?1/22??tn2u?g?xi??? ?4/3??2??R?i?1jHi??xij1??j1??j2jjj??Qi?1/2?2ui?1/2(Qij?Qi?)?gA?Bu1i?1/2(Zi?1?Z1)

??t??对任一河段i(i?1,2,??n)，可得到方程组：

?CiZi?CiZi?1?Qi?Qi?1?Di (10) ??EiQi?GiQi?1?FiZi?FiZi?1?Hi对每一河段可列出两个线性代数方程，再加上上下游边界条件，构成完备的封闭方程组，采用追赶法可求得各个断面的水位流量。 （2）边界条件

根据上有下游边界条件类型的不同可以写成如下两种追赶形式：

·上游水位边界条件Z*1下游水位（或流量）边界条件Zn?1?Zn?1(t)（或?Z1(t)；

**，追赶形式为： Qn?1?Qn?1(t)）

?Z1?P1?R1Q1?Q?L?MQ222?1?Z2?P2?R2Q2?????Qi?Li?1?Mi?1Qi?1 (11) ?Z?P?RQi?1i?1i?1?i?1????Qn?Ln?1?Mn?1Qi?1?Z?P?RQn?1n?1n?1?n?1**??Zn?1?Zn?1(t)(or,Qn?1?Qn?1(t))式中，P,R,L,M为已知系数，依据上述方程组，可逐步由下边界水位或者流量，推算得到上游各个断面水位流量值。

·上游流量边界条件Q?Q(t)；下游水位边界条件Z1*1n?1 ?Zn?1(t)，追赶形式为：

* 3

?Q1?P1?R1Z1?Z?L?MZ222?1?Q2?P2?R2Z2?????Zi?Li?1?Mi?1Zi?1 (12) ?Q?P?RZi?1i?1i?1?i?1????Zn?Ln?1?Mn?1Zi?1?Q?P?RZn?1n?1n?1?n?1*??Zn?1?Zn?1(t)式中，P,R,L,M为已知系数，依据上述方程组，可逐步由下边界水位，推算得到上游各个断面水位流量值。

7.2.2 单一河道一维水质模型 （1）控制方程的离散与求解

对方程(3)进行离散，空间差分采用隐式迎风差分格式。顺流时(从断面i流向i+1) 有：

?(AC)(AC)i?(AC)in ??t?t?(AUC)(QC)i?(QC)i?1 ??x?xi?1?CCi?1?CiCi?Ci?1?1(AEx)?[(AEx)i?(AEx)i?1] ?x?x?xi?1?xi?xi?1?KAC?S??Ki?1Ai?1/2Ci?Si?1

得到统一形式的差分方程：

aiCi?biCi?1?ciCi?1?di (13)

式中，ai,bi,ci,di为系数，分别表示为：

bi?{[Qi?1/2,0]?(AEx)i?1/2}??t

?xi?1(AEx)i?1/2(AEx)i?1/2??t???t?KiVi?t ?xi?xi?1ai?Vi?[Qi?1/2,0]?t?[?Qi?1/2,0]?t?(AEx)i?1/2}??t ?xici?{?[Qi?1/2,0]?

4

di?Vi?Cin?[Qi,0]?Si??xi?t?[?Qi,0]Si?1??xi?1?t 方程（13）两边同时除以ai得到：

nCi?EiCi?1?FiCi?1?Gi (14)

在顺流情况下，各河段差分方程可写成：

?C2?E2C1?F2C3?G2?C?EC?FC?G32343?3??? (15) ?Ci?EiCi?1?FiCi?1?Gi????Cn?1?En?1Cn?2?Fn?1Cn?Gn?1??Cn?EnCn?1?FnCn?1?Gn对首断面给定第一类边界条件，对末断面给定第二类边界条件，可得到如下封闭的方程组：

?C1?C1*(t)??C2?E2C1?F2C3?G2?C3?E3C2?F3C4?G3???? (16) ?Ci?EiCi?1?FiCi?1?Gi????Cn?1?En?1Cn?2?Fn?1Cn?Gn?1??Cn?EnCn?1?FnCn?1?Gn?Cn?1?En?1Cn?Gn?1?

对方程组（16）采用追赶法可容易求得n,n?1,n?2,?,3,2等断面的扩散质的浓度。 （2）参数确定

·纵向分散系数E的确定

X

EX与水流流速、水面宽度成正比，与水深成反比，常采用下面的经验公式：

Ex??C0?2q

式中，C0?c/g是无尺度谢才系数，c为谢才系数，? =B/h为宽深比，q为单宽流量，? ?0.01?为经验常数。

·降解系数K的确定

5

可采用监测资料对降解系数进行率定，或根据经验得到。

7.2.2 应用实例［］

三峡大坝位于宜昌县三斗坪中堡岛,葛洲坝位于南津关下游的宜昌市境内.两坝间水域处于鄂西山区向平原的过渡地带,周围地形地貌呈西北高东南低之势.三峡大坝至葛洲坝河段长38Km,两坝间江面宽210m至1500m,大部分处于西陵峡谷中。该实例建立了两坝间水量水质模型,分别运用大坝一期围堰及二期围堰施工期间的同步水文水质实测资料对模型进行了率定和验证,取得了较好效果.此模型可用以预测大坝施工期间及投入运行后两坝间水流及水质特性的变化. （１）水量模型率定

利用1996年三斗坪、白庙子及黄陵庙等断面水文观测资料率定糙率，得到各子河段的糙率，率定结果显示糙率分布取值范围为0.036至0.050。1996年白庙子及黄陵庙等断面实测水位过程线及计算水位过程线图略。 （２）水量模型验证

采用太平溪断面1998年实测流量作模型验证的上边界条件，采用葛洲坝坝前断面1998年实测水位作模型验证的下边界条件，对1998年两坝间水流进行模拟，通过对各水文观测断面的水文要素的观测值及计算值进行比较可见，吻合程度较好，因此此水量模型可用来模拟两坝间的一维水流情况，1998年白庙子及黄陵庙等断面实测水位过程线及计算水位过程线见图２。

H(m)74.0072.0070.0068.0066.0064.0006012018024030036068.0066.0064.00060120180240300360三斗坪断面H(m)72.0070.00黄陵庙断面T(daY)T(daY)

图2模型验证各断面实测及计算水位过程线

(３) CODMn模型验证

高锰酸盐指数的降解系数由实测资料，取经验值，不采用模型率定。采用1998年7~12月的CODMn浓度监测资料对CODMn模型参数可靠性进行验证，模拟时CODMn污染源分点源和面源两种情况。部份断面CODMn水质因子实测及计算浓度值见图3。

6

三斗坪断面COD(mg/L)3.00实测值计算值COD(mg/L)3.00白庙子断面实测值计算值2.002.001.000306090120T(day)1501801.0003060

COD(mg/L)90120T(day)150180

COD(mg/L)3.00黄陵庙断面实测值计算值南津关断面实测值计算值3.002.002.001.000306090120T(day)1501801.00

0306090120T(day)150180

图3 模型验证CODMn浓度实测及计算值

7.3 平原河网水量水质模型

河网地区是中国社会、经济、文化的发达地区，在国家经济发展中占有举足轻重的地位。随着地区经济的进一步发展、居民物质精神生活水准的进一步提高，水资源问题日益突出。如中国著名的长江三角洲、珠江三角洲地区，已成为水质型缺水区，水资源、水环境问题已成为制约经济社会发展、事关区域可持续发展的重要因素。生产实践的迫切需要使得河网地区的水环境保护研究呈不断深入、系统之势。作为基础，河网水流水质模拟方法是进行区域环境规划、环境管理等的必备工具，在水环境问题研究中占有重要地位。因此，该领域的研究一直是环境科学研究人员、环境管理决策部门十分关注的重要问题。在物理模型、数学模型这两大手段中，因为数学模型具有经济、快捷、实用等优点，加之水网地区河道密布等客观条件的限制，现阶段只能采用数值方法模拟水网地区的水流运动及污染物输运规律。但由于理论、技术及各种客观条件的制约，目前水流特别是水质数值模拟精度不是十分令人满意。

7.3.1 平原河网水量模型

按河网水流的控制方程及对河网的概化处理方式不同，河网地区水流数值模拟方法可分为两大类：第一类为常用的一维圣维南方程组数值解法；第二类为所谓的“组合单元解法”。，其中，一维圣维南方程组数值解法又可分为直接解法和间接解法两种。在直接解法中，较有代表性的有文献[G.Noseda, Mathematical Model of Unsteady Flow in Open Channels Networks, Proceedings of the International

Symposium on Unsteady Flow in Open Channels,1976]，[J.J.R. Williams, T.R.E. Chidley, Nonliear Analysis of Unsteady Flow in Open Channel Networks, Proceedings of the International Symposium on Unsteady Flow in Open Channels,1976]提出的方法，该类方

法将计算断面交替取为水位和流量断面，对河网在所有计算断面上统一对一维圣维南方程组差分离散并求解。但该方法未知数数量较多，在河网规模较大的情况下，因为河道的交叉衔接，形成的矩阵是一个不规则、不对称的大型稀疏矩阵。为减少存贮，J.J.Dronkers 于1976年提出间接解法[5] 的思想，以后又有许多学者对其作了进一步完善。为提高计算效率，需缩小矩阵规模，中山大学数力系1976年提出了河网非恒定流隐式方程组稀疏矩阵解法[6李岳生

7

，该方法从河网矩阵本身的特点出发，能够

有效地节省存储并提高计算速度，但矩阵中需包含所有断面的未知数，方程规模仍然较大，其实际使用也受到限制。而间接解法是将断面未知数往交汊点集中，待求出交汊点未知数后，再求解各单一河道未知变量，计算效率较高。间接解法的思想首先由荷兰水力学专家Dronkers于1976年提出，以后又有许多学者对此方法进行进一步的完善，相继提出了河网非恒定流的二级解法[徐正凡，明渠非恒定流[M].武汉水利电力学院，1983年.]、三级解法[8张二骏，河网非恒定流的三级联合解法，华东水利

][9]

学院学报，1982年第1卷：1~12和四级解法吴寿红，河网非恒定流四级解法，水利学报，1985年，第8期：42~50。此类间接解法中，以三级解法最为常用[][][]。其基本求解思路可概括为“单一河道—交汊点—单一河道”，即（1）先将单一河道划分为若干子河段，在计算断面上对一维Saint-Venant方程组进行有限差分运算，得到各单一河道差分方程组，进行消元计算，得到单一河道首、末断面间流量与水位的相互关系；（2）根据河道交汊点水量守恒方程，得到并求解交汊点水位方程组，得到所有交汊点水位；（3）根据交汊点水位返回各单一河道，求得各计算断面水位、流量值。

“组合单元解法”[10]由法国水力学专家Jean A. Cunge于1975年首次提出，国内也有研究者采用此方法进行了水网地区的水力模拟[11][12]。此方法的基本思想是：将河网地区水力特性相似、水位变幅不大的水体概化成单元。取单元中心的水位为代表水位，采用谢才公式模拟单元间流量交换，根据水量守恒建立每一单元的微分形式的水量守恒方程，离散并得到以单元水位为自变量的代数方程，辅以边界条件，可求得各单元水位、单元间流量。

在上述两类方法中，组合单元解法对河道进行了简单概化，以单元为计算单位，计算相对简单，但模拟精度相对较低，仅实用于大尺度水域的水力模拟。而Saint-Venat方程组数值解法可以精确计算每一条河道的水流状况，所以成为目前河网水力模拟的主流方法，其中又以三级解法最为常用[13][14][15][16]。 当采用三级解法模拟长系列大范围河网水流特性时，需对河道进行概化处理，即将等级较小的河道概化为一条“概化”河道，要求该概化河道与被概化河道的过流能力、调蓄量相当。同时，为考虑降雨对河网水力特性的影响，还必须对河道包围的陆域面积进行产汇流计算，将产流量以包围陆域的河道长度为权重分配到周围河道[17]

。由于三级解法以单一河道为模拟对象，同时考虑了降雨等的影响，计算精度较高，可以满足水利、航运、环境保护等的需要。

等，河网不恒定流隐式方程组稀疏矩阵解法，中山大学学报（自然科学版），1979年3月

]

（1）控制方程的离散与求解

采用四点隐式差分格式离散方程，得子河段差分方程：

CiZi?CiZi?1?Qi?Qi?1?Di (17)

EiQi?GiQi?1?FiZi?FiZi?1?Hi (18)

式中Ci,Di,Ei,Fi,Gi,Hi由时段初值及河道特征求得。

（2）节点连接条件

水流运动在河网各节点上应满足质量守恒及能量守恒，即满足以下两个连接条件： · 质量守恒条件

进出某一节点的流量与该节点内水量蓄量的增减相平衡，定量表示为：

?Qi?1mij?1j?(?K??K)/?t (19)

式中K为节点编号，m为流入(流出)节点K的河道数，?K为节点蓄量，Q为流量(流入为

8

正；流出为负)。若节点汇合区容积与子河段容积相比可忽略不计，则此节点称为无调节节点，方程(19)可简化为：

?Qi?1mi?0

· 能量守恒条件

不计节点处能量损失有：

Zi?Zj i=1,2,?m, j?1,2,?m (20)

(3) 边界条件

有三种类型的边界条件： 水位边界条件

即在边界河道上给定水位随时间的变化过程： Z?Z(t)； 流量边界条件

即在边界河道上给定流量随时间的变化过程：Q?Q(t) ；

水位流量关系

当边界河道上有水工建筑物(如水闸、堰、堤坝等)时，通常给定水位流量关系：Q?Q(Z)。

（4） 方程的求解

利用消元法，方程(3)、(4)经递推运算，写成如下形式：

Qi??i??iZi??iZn?1 (21)

Qi?1??i?1??i?1Zi?1??i?1Z1 (22)

式中，?i,?i,?i,?i,?i,?i为追赶系数，可逐步递推求得。利用公式(7)依次由末断面向首断面递推，可将各断面流量表为该断面水位及末断面水位的函数。利用公式(22)依次由首断面向末断面递推，可将各断面流量表为该断面水位及首断面水位的函数。特殊的，单一河道首末断面分别与节点相连，分别对应如下追赶方程：

Q1??1??1Z1??1Zn?1 (23) Qn?1??n?1??n?1Zn?1??n?1Z1 (24)

将节点各支流相应的(23)或(24)代入公式(19)，并将节点各相邻断面水位统一表示成节点水位，得节点方程：

fi(Zi,Z邻1,Z邻2,?Z邻KI)?0 i=1,2,?,M

式中，M为节点总数，KI为第i个节点相连河道数。

当某单一河道为边界河道时，分三种类型分别有：对于第一类边界条件，式(24)中Z1已知；对于第二类边界条件，式(23)中Q1已知；对于第三类边界条件，已知线性关系Q1=Q(Z1)，代入(23)消去Q1。无论何种边界条件，最终可增加一

9

水位方程，封闭节点方程组fi(Zi,Z邻1,Z邻2,?Z邻KI)?0，可解得各节点水位，再返回单一河道方程，最终求得各断面水位及流量。

7.3.2 平原河网水质模型

(1) 单一河道控制方程的离散

采用如下网格对方程(3)进行离散，空间差分采用隐式迎风差分格式，由于河网地区水流流向可能谁随时间变化，根据流向的不同，分顺流及逆流两种状况给出差分形式：

Ci?1 Ci Ci?1

1 2 3 ? i?1 i i?1 ? n?1 n ?xi?1 ?xi

图4水质模拟网格布置示意图 顺流时(从断面i流向i+1) 有： ?(AC)(AC)i?(AC)ni?t??t ?(AUC)(QC)i?(QC)i?1?x??x i?1?(AE?C1Ci?1?CiCi?Ci?1?xx?x)??x[(AEx)i?(AEx)i?1] i?1?xi?xi?1?KAC?S??Ki?1Ai?1/2Ci?Si?1

逆流时(从断面 i+1流向i) 有：

?(AC)(AC)i?(AC)ni?t??t ?(AUC)(QC)i?1?(QC)?x?i?x i?(AE?C1C?CiC?Ci?1xx)??x[(AEi?1x)i?(AEi?x?x)i?1] i?xi?xi?1?KAC?S??KiAi?1/2Ci?Si

为考虑流向顺逆变化的影响，引入流向调节因子rc及rd，对于断面i表示为： 10

为考虑流向顺逆变化的影响，引入流向调节因子rc及rd，对于断面i表示为：

?Qw?(Qi?Qi?1)/2Q?e?(Qi?Qi?1)/2 ?rcw?(Qw?Qw)/2Qrce?(Qe?Q?we)/2Qe?rdw?(Qw?Qw)/2Qr?(QQ

?wdee?e)/2Qe?rc?rd?0,（当Qw,Qe?0时）得到任意流向下，统一形式的差分方程：

?iCi?1??iCi??iCi?1?Zi 式中，?i,?i,?i,Zi为系数，对于一般断面(i=2,3,…,n-1)分别表示为： ?i??[rcwDww?rdwDpw?Fcw]??t/V

?i??[rcwDww?rdwDpw?rceDpp?rdeDep?Fcp?Fdp]??t/V?rcwKi?1?rdeKi]??t?10.

?i??[rceDpp?rdeDep?Fde]??t/V Zi??Cni?[rcwSi?1?1A?rdeSi?1i?1/2A]??t

i?1/2对于首断面(i=1)，逆流时有：

??1?0?????t1?[rdeDep?Fdp]?x?rdeKi?t?rdeiAi?1/2 ???t??1??[rdedep?F

de]??xiAi?1/2???Z1?Cn1?r?tdeS1Ai?1/2

对于末断面(i=n)，顺流时有：

???n??[r?tcwDww?Fcw]?An?1/2??t ???n?[rcwDww?Fcp]?x?rcwKi?1?t?rcwi?1Ai?1/2

???n?0???Zn?Cnn?r?tcwSn?1Ai?1/2

式中，

(25)

11

Ai?1/2?(Ai?1?Ai)/2, V1?Ai?1/2??xi?1, V?rcwV1?rdeV2

??Dww?(AEx)i?1/?xi?1, ???Dpw?(AEx)i/?xi?1,Ai?1/2?(Ai?1?Ai)/2V2?Ai?1/2??xi

Dpp?(AEx)i/?xiDep?(AEx)i?1/?xi

??Fcw?(Qi?1?Qi?1)/2, ???Fdp?(Qi?Qi?1)/2,Fcp?(Qi?Qi)/2Fde?(Qi?1?Qi?1)/2

将差分方程(25)用向量形式表示为：

AcCc?Bc (26)

式中Ac为一?(n?1)?n?阶的三对角的系数矩阵，Cc为此河道断面平均浓度的n维列向量，Bc为已知的(n-1)维列向量。(26)式中含(n-1)个方程，n个未知数，方程组不闭合，需引入节点方程及边界条件。 （2）河网节点方程

对于污染物充分混合的节点，从时刻n?t到(n?1)?t，根据节点有无调蓄作用，分别给出相应的节点方程。

若节点有m条单一河道，其中流入节点的河道m1条，流出节点的河道m2条，

m1及m2随流场变化而变化(m1+m2=m)。

若节点本身具有调蓄作用，有：

Cout,i?CN， i=m1?1,m1?2,?,m (27)

?Qi?1m1in,iCin,i??i?m1?1?Qmout,iCout,i (28)

1jj?K?t?NCN??NCNe?SN?t??式中，CN 为节点的浓度，Cout,i为流出节点的第i条河道与该节点相邻断面的污染物浓度，Cin,i为流入节点的第i条河道与该节点相邻断面的污染物浓度，Q分别为与之相应的流量，?N为节点的蓄量，SN为节点的污染源加入项，方程(28)右端括号中末项为时段初节点污染物降解至时段末的残留量。(27)式称为充分混

12

合假定，(28)式反映节点内污染物的质量守恒。

当节点调蓄可忽略时，称为无调蓄节点，有：

?N?0

（3）模型的动态耦合数值计算

在流场已知时，由(26)-(28)式加上边界条件及初始条件，构成封闭的线性代数方程组，即可求解。但对于大型复杂河网来说，其规模巨大，较为经济的做法是：将方程(26)作适当的运算，得单一河道首末断面浓度关系方程，并与充分混合假定(27)及边界条件，一并代入方程组(28)，形成封闭的节点方程组。

由于方程组(26)系数矩阵为三对角矩阵，通过消元法，可得到任一单一河道首末断面浓度关系方程：

Ci,ni?fiCi,1，i = 1, 2, ??????M (29)

式中M为河道总数，ni为第i条河道断面数，Ci,1,Ci,ni 分别为第i条河道首断面及末断面浓度，f i为线性函数关系，据上式，所有与节点相邻的河道首末断面处的浓度皆可表示为节点浓度的线性函数。将方程(29)及河网入流边界条件等代入方程(28)，并考虑方程(28)，即可得到整个河网各节点浓度的代数方程组： FjC1,C2,?Ci,?CK?bj,????j?1,2,?,K (30)

式中，K为河网节点总数，方程组(7)还可用矩阵表示为：

ANCN?BN (31)

式中，AN为系数矩阵，CN为节点浓度列向量，BN为包括污染源在内的已知项的列向量，求解上述方程组可得各节点浓度列向量。再将与单一河道相应的已知方程(27)代入方程(26)，也即在方程组(26)中增加一类似于方程(27)的反映入流断面与相连节点浓度关系的代数方程，即可封闭河道方程(26)，从而求得该河道各断面浓度。

7.3.3 算例（南通河网水量水质计算）

南通市境内水系以如泰运河为界，分属淮河和长江两大水系。长江水系流域面积5974.5km2，包括通扬运河、焦港、如海运河、通启运河、九圩港、通吕运河；淮河水系流域面积为2524.1km2，主要河道有通扬运河、通榆运河、拼茶运河、如泰运河等。水流流向受自然因素和人为因素的影响，流向多变。该地区地势平坦，河道纵横交织成网，南与长江相连，东与黄海相接，水流既受地形的支配，又受长江、黄海及内部闸门人为控制，水流留向顺逆不定。

根据南通河网的水文情势、水力特性及河道大小等因素，对其进行概化。概化后的河网如图5所示，共有92个节点，145个内部河道，37条边界河道，其中26个水位边界，11个流量边界，594个计算断面，因缺乏流量资料，选用1988年5-9月货隆、搬泾、金沙、丁堰等水位资料率定模型。用1994年元月水位观测资料进行模型验证，图6为模型率定情况下，观测水位及计算水位过程线。

p I C T U

13

图5 南通河网概化图 R E P I C P

图6 模型率定：实测及计算水位过程线（实线：实测；虚线：计算）

利用1994年1月水质监测资料率定模型，5-9月水质监测资料校核模型。部分断面校核成果如图7所示。由图可知，计算精度尚可。

图 7 模型验证：COD MN实测及计算比较

14

本文来源：https://www.bwwdw.com/article/piu2.html