http://www.paper.edu.cn
河网水动力模型及水质模拟研究展望
陈新永,朱静
(河海大学,环境科学与工程学院, 江苏南京 210024)
(Chenxinyong2007@163.com;biye2004@eyou.com)
摘要:论述了河网地区 4 类水动力模型及其基本思想,对水质模拟的方法理论、
发展过程和生产应用进行简单分析。结果表明,三级联解法是河网地区水质计算
的主流方法,组合单元解法有其良好的应用条件。在水动力模型研究发展基础之
上,指出减小模拟假设误差、改进和设计新的计算方法、数值模拟可视化将成为
河网水质模拟的主要发展方向。
关键词:河网;水动力模型;水质模拟
1. 前言
平原河网地区湖泊密布、河网发达,是我国社会、经济和文化较发达的地区。
然而,随着城市人口、经济的发展,各种工业废水和生活污水大量排入河网,致使
水环境质量日益恶化,成为困扰城市环境和经济可持续发展的一个主要问题。
所有这些对平原河网的灌溉、排涝、防洪、规划和污染控制等提出了一系列新的
研究课题。为此,人们不仅日益重视该地区水资源保护的研究,同时还采取相应的
水流、水质模拟方法对其进行环境规划和管理.但由于理论、技术及各种客观条
件的制约,目前水流,特别是水质数值模拟精度不十分令人满意.由于河网模拟区
域范围很大,大多数情况下只能采用数值方法进行模拟,其核心问题是河网数学
模型的建立及求解,其中水动力模型是其他模型(如水质模型、水环境容量模型等)
的基础。
河网水动力数学模型大体可以分为节点-河道模型、单元划分模型、混合模
型以及人工神经网络模型 4 类[1]。平原河网不同于单一河流的特点,在于河网错
综复杂性以及由此带来的方程组离散和求解上的困难,这是多年来人们研究河网
问题的一大难点,有关文献[2 ,3]给出了方程组的离散及求解推导的详细过程。本
文仅对各类河网水动力模型进行综述分析,并在此基础之上论述河网水质模型的
1
研究与进展,为河网网地区水流水质计算提供科学参考。
http://www.paper.edu.cn
2. 河网水动力数学模型
2.1 节点-河道模型
该模型的基本思想[4]是:将河网中的每一河道视为单一河道,其控制方程均
为一维 Saint-Venant 方程组;河道连接处称为节点(汊点),每个节点处均应满足
水流连续性方程和能量守恒方程。求解由边界条件、Saint-Venant 方程组和汊
点衔接方程联立闭合方程组,即可得到各河段内部断面的未知水力要素。
河网非恒定流问题水动力模型的控制方程为 Saint-Venant 方程组:
A
∂
t
∂
Q
∂
t
∂
+
+
Q
∂
x
∂
∂
x
∂
=
q
α
Q
2
A
+
gA
Z
∂
x
∂
+
gA
=
0
QQ
K
2
式中:x 为沿水流方向的距离;t 为时间;A 为过水断面面积;Q 为流量;
Z 为水位;q 为单位河段长度的旁侧入流流量;K 为流量模数;α为动量校
正系数;g 为重力加速度。
节点-河道模型原则上可以求解任何水网的水力参数,不同的求解方法各有
优缺点。目前常用的求解方法有:特征线法[5],有限差分法[6],直接解法[7],间
接解法[8-12],矩阵标示法[13],有限体积法[14]等。
2.2 单元划分模型
在大型复杂河网计算中,尤其在湖泊、水库较多的情况下,单元划分模型更
具优点,其基本思想[15]是:将水力特性相似、水位变化不大的某一片水体概化为
一个单元,单元间流量交换的媒介是连接河道,其本身无调蓄作用。假定单元间存
在两种连接方式:在无水工建筑物或障碍物(不存在局部水头损失)的情况下认为
是河型连接;堰型连接存在局部水头损失,又可分为自由出流和淹没出流两种形
式。
单元划分模型的主要控制方程为,对任意单元 i 建立时段(
t)1n(,tn
∆
∆+
)内
2
的质量平衡方程:
As
i
i
dZ
dt
=
)t(P
i
∑+
k
)Z,Z(Q
k
k, i
i
http://www.paper.edu.cn
式中, iAs 为水位为 iZ 时单元 i 的水面面积; )t(Pi 为 i 单元的降雨强度;
)Z,Z(Q
k
k,i
i
为从任意相邻单元 k 流入或流出单元 i 的流量。
取单元几何中心的水位为单元代表水位,给出水位与水面面积关系。将计算
河网分解为一定数量的单元,再进行分组,然后确定各单元间的连接类型。对每个
单元给出微分形式的质量平衡方程,经有限差分法离散后得到以单元水位为基本
未知量的方程组,进而求解各单元的代表水位和单元间流量。
2.3 混合模型
姚琪等人认为,运河河网地区地势平坦,区内无长大的天然河流;大多数河流
坡降平缓,流量很小;农灌渠道不计其数,再加上泵站、水闸、船闸等水利控制工
程,使河网的水力学描述更加复杂,因而在建模工作中完全如实地模拟如此庞大
复杂的水系几乎是不可能的。节点-河道模型和单元划分模型都不能很好地适应
运河水网的特性,前者失之过繁,后者失之过简。将节点-河道模型和单元划分模
型中与平原河网特性相适应的优点综合起来,并避免其不相适应的缺点,构成新
的数学模型,即混合模型。
建立混合模型的基本思想[16]是:将平原河网的水域区分为骨干河道和成片
水域两类,对骨干河道采用节点-河道模型;对成片水域采用单元划分的方法将其
划分为单元,再引入当量河宽的概念,把成片水域的调蓄作用概化为骨干河道的
滩地,将其纳入节点-河道模型一并计算。
2.4 人工神经网络模型
人工神经网络是采用计算机来模拟生物神经网络的结构和功能,由大量神经
元构成的并行分布式系统。人工神经网络与平原河网在结构上有许多相似之处,
两者都是由各个内部单元通过并联或串联形成一个相互制约的整体网络结构,通
过调整系统内部各个“神经元”之间的相互作用可以达到系统输入变量和输出变
3
http://www.paper.edu.cn
量之间的最优化或平衡,因此,人工神经网络理论可用于复杂河网水动力模型的
数值模拟[17]。
实际运用人工神经网络时,由于天然河网十分复杂,河流湖泊众多,首先应
对河网进行概化[18]。李荣等人从河网水流运动连续方程和漕蓄方程出发,对河网
作了概化[19]。经过上述概化,整个河网由河网水源的输入、河网内部相互串联和
并联而成的水库、河网输出三部分构成,如图 1 所示。第一层和最后一层节点的
输入输出为简单的线性关系,利用节点水量平衡方程和整个河网水量守恒关系保
证各个节点和整个河网水量守恒。
模型输入向量
中间层
模型输出向量
图 1 河网概化的神经网络模型
神经网络模型与水动力模型的有机结合可以减少水动力模型计算量、加快计
算速度;同时水动力模型为神经网络模型提供参考,在重要河段或区域由于资料
缺乏时得到广泛应用。神经网络模型是“黑箱”模型[19],反映网络内部节点之间
关系参数的物理意义模糊,就河网水动力模拟而言,它实际上反映的是河网各种
因素(如河网结构、河道糙率、断面形状、调蓄容量等)的综合影响,这就给模型
的验证带来困难,同时有可能影响计算结果的精度。
3.分级解法
分级解法是目前河网水流求解最主要的方法,该方法首先由荷兰学者提出.
4
http://www.paper.edu.cn
分级解法的基本思想是[20],先将未知数集中到汊点上,待汊点未知数求出后,再
将各河段作为单一河段求解.分级解法按方程组的连接形式,又可以分为:二级解
法[21]、三级解法[22,23]四级解法[24]、汊点分组解法[25]和树型河网分组解法[26].先就
其中几种解法进行简介。
(1)二级解法:二级解法的基本思想是将所有的边界方程和河段方程一起构
成一个封闭的方程组,也即是二级连接方程组,求解这样的方程组,便可以求得河
网各河端未知数,然后再利用微段方程,求出全部内部计算断面的未知数.
(2)三级解法:三级解法是在二级解法的基础上提出来的.它的基本思想是在
二级解法的基础上,将所得到的河段方程自相消元,可以得到一对以水位或流量
为隐函数的方程组,方程组经转化变形可以直接代入相应的汊点和边点方程,消
去其中的流量改正值,则剩余含有 2Nr 个未知的水位改正值变量的 2Nr 个方程,
求解连接矩阵得到各汊点上各断面的水位改正值,回代各河段方程得到汊点各断
面的流量改正值,再回代节点的流量和水位值的递推公式,得到所有各断面上的
水位和流量改正值,从而可以得出各断面的水位和流量。
(3)四级解法:四级解法是在三级解法的基础上,进一步从三级连接方程组中
分离出外边界方程和汊点能量衔接方程,最后由剩下的方程构成四级连接方程组.
然后,将第一步所得结果和汊点能量方程代入汊点水量平衡方程中,消去除各汊
点第一断面水位改正值以外的全部未知数,这样形成的方程组即是所谓的四级连
接方程组,四级连接方程组共含有 Nj 个方程,Nj 个未知数可以求解。四级连接方
程组一经解出,将其乘上相应的系数,便可以得到三级解,再按三级解法相同的方
法求出各断面的水位和流量值。
(4)汊点分组解法:汊点分组解法首先由我国学者李义天提出[25]。该方法首
先将河网中的汊点分为 NG 组,用 ng 表示汊点组的序号(1≤ng≤NG),分组时,要求
每一汊点组最多与两个汊点组相连每组汊点中的河段,可以是连接本组汊点的河
段,也可以是一端连接本组汊点,另一端连接前一组或后一组汊点的河段,对每一
组汊点均可以形成分组后的汊点方程组。联立求解方程组,并回代求解可得河段
内部各断面的水位和流量增量。
在直接解法中,系数矩阵阶数为河网中选取的计算断面数的两倍,这对于大
型河网的求解几乎是不可能的。二级解法的系数矩阵的阶数为 4Nr,三级解法系
5
http://www.paper.edu.cn
数矩阵为 2Nr,四级解法系数矩阵的阶数为 Nj,汊点分组解法系数矩阵的阶数仅
与分组以后汊点数最大值相等,因而使大型河网问题在一般计算机上也能实现。
树型河网分组解法不涉及解大型矩阵,但该方法不适用于环形河网。因此汊点分
组解法是目前最优的解决河网问题的方法,值得我们借鉴使用。
4. 河网水质模拟研究现状
目前,河网水质模型及模拟方法都比较成熟。在水质模型中,美国的
QUAL-II(E)模型、传统 CSTR(Continuously Stirred Tank Reactor)模型等被广
泛用于河流水质模拟。WASP 模型及 Carlos 的河网水质模型[33]则只适用于小型河
网。河海大学的 H w q now 模型能反映水体在不同溶解氧状态下的水质变化,且
适用于大型河网[34]。组合单元水质模型也是一种有效的河网模型[35]。清华大学
的水文-水质模型能应用于环状河网[29]。
在水动力模型的基础之上,根据水质控制方程及河道概化方式,河网地区水
质模拟方法可分为两类[31]:第一类为常用一维纵向分散方程求解法;第二类为
组合单元解法。
第一类水质模拟中,三级联解法是水质模拟最为常用,其求解思路[33]:(a)
对每条河道的扩散质一维对流纵向分散方程进行有限差分离散,经消元计算得到
各单一河道出流断面质量浓度与入流断面质量浓度间的线性关系;(b)在交叉口
引入均匀混合假设——流出交叉口的断面质量浓度等于流入交叉口各断面质量
浓度的平均值,导得以各交叉口节点质量浓度为未知变量的方程组,求解得各交
叉口质量浓度;(c)返回单一河道计算各断面质量浓度.此方法被广泛应用于河网
地区环境模拟和环境规划[25]研究中,并已成为河网地区水质计算的主流方法。
第二类水质模拟中,其组合单元解法求解思路[15]:按水力水质特性相近的原
则,将河网划分成若干单元,要求划分后的单元内部,水位流速变化均匀、水质浓
度无悬殊变化。单元间污染物质的交换通过随水体对流及交界面上的扩散而进
行。单元污染物浓度以形心处浓度代替,它随时间的变化,取决于交界面上的对流
扩散输运、单元内污染物的降解及源项的加入。此方法由于单元概化使精度受损,
仅适用于大尺度水环境规划,但该方法以各单元物理量为直接求解对象,为环境
6
http://www.paper.edu.cn
容量计算、水质规划及反问题[32,33]研究提供了方便。
河网水质模拟实际上是一个复杂的过程,不单单与污染质有关,更与其本身
性质息息相关,如油污染质不同于一般污染质,它不易溶于水,易挥发等特性,
使之难于用一般水质模拟方法进行模拟分析;另外,像 COD、BOD、有机物等污
染物质,本身还有反应变化,这就增加了数值模拟的难度,无法对其进行精确模
拟。因此,河网水质数学模拟仍存在大量问题,有待进一步研究
5. 结语
目前,无论是水力模拟还是水质模拟,都基本上能满足水利、航运、环保等
生产实践的需要。随着水动力模型的不断发展,水质模型也在发展。模型主要取
决于以下几个因素:工程实际的需要[36-38](如水环境质量、水污染控制等),计算
数学理论和方法的成熟和应用[39-40] (如数学物理方程、离散方法、偏微分方程组
的求解方法等),计算机软硬件技术的发展[41] (如内存,CPU,GIS,多媒体技术等),
以及用于率定模型数据的量测精度等。
水质模拟的主要控制方程为对流扩散方程,寻求一种更能真实有效反映平原
河网水力水质特性的数学模型和高效的求解方法仍将是今后河网水质模拟的主
要发展方向:
a.减小模拟假设误差。从河网水力水质特性出发[33],使模拟假设有效地反
映具体河流实际情况(如水位、流量、污染物浓度等);减少被概化河流个数,
选择水质状况、水环境功能相近的河流分组概化。
b.改进和设计新的计算方法。随着计算机软硬件技术的飞速发展,内存、存
储空间和 CPU 运行速度将允许在数值离散上采用很小的时间和空间步长,这对于
提高计算精度来说,比采用较为复杂的高阶算法更为有效。
c.数值模拟可视化。计算机性能的提高为复杂河网的数值模拟提供了广阔的
前景,借助计算机的图形处理能力和多媒体技术,实现计算可视化[42-43],必将大
大提高河网计算数据的处理速度和质量。
7
参考文献:
http://www.paper.edu.cn
[1],程开宁.河网水流模拟计算数学模型探讨[J].东北水利水电,2005 年第 3 卷.
[2],汪德爟,计算水力学理论与运用[M],河海大学出版社,1989 年 9 月.
[3],金忠青,王玲玲.N-S 方程的数值解及紊流模型[M],河海大学出版社,2005 年 8 月.
[4],卢士强,徐祖信.平原河网水动力模型及求解方法探讨[J].水资源保护 2003 年第 3 期.
[5] , 顾 元 木 炎 , 王 尚 毅 . 具 有 岔 道 支 流 的 明 渠 不 恒 定 流 计 算 [J]. 天 津 大 学 学
报,1995(1):11~17.
[6] , 谭 维 炎 , 胡 四 一 , 王 银 堂 , 等 . 长 江 中 游 洞 庭 湖 防 洪 系 统 水 流 模 拟 [J]. 水 科 学 进
展,1996(12):337~345.
[7],李岳生.河网不恒定流隐式方程组稀疏矩阵解法[J].中山大学学报 1977(3):7-13.
[8],Dronkers J J.河流近海区和外海的潮汐计算[J].水利水运科技进展,1976(增刊):24~
67.
[9],张二俊,张东生,李挺.河网非恒定流三级联合算法[J]华东水利学院学报,1982(1):1~
13.
[10],吴寿红.河网非恒定流四级解算法[J].水利学报,1985 (8):42~50.
[11],李义天.河网非恒定流隐式方程组的汊点分组解法[J]水利学报,1997(3):49~57.
[12],侯玉,卓建民,郑国权.河网非恒定流汊点分组解法[J]水科学进展,1999(3):49~52.
[13],李光炽,王船海.大型河网水流模拟的矩阵标识法[J]河海大学学报,1995(1):37~42.
[14] , 曾 凡 棠 , 黄 水 祥 . 珠 江 三 角 洲 潮 汐 河 网 水 环 境 数 学 模 型 评 述 [J]. 海 洋 环 境 科
学,2000(11):47~50.
[15],韩龙喜,张书农,金忠青.复杂河网非恒定流计算模型——单元划分法[J].水利学
报,1994(2):53~56.
[16],姚琪等.运河河网水量数学模型的研究和应用[J].河海大学学报,1991(7):9~17.
[17],DibikeY B,Olomatine D,Abbott M B.On the encapsulation of numerical- hydraulic
models
Hydraulic
Research,1999,37(2):141~167.
[18],吴寿红.河网非恒定流四级解算法[J].水利学报,1985(8):42~50.
[19] , 李 荣 等 . 河 网 水 情 预 报 的 神 经 网 络 模 型 及 应 用 [J]. 应 用 基 础 与 工 程 科 学 学
报,2000(2):179~186.
[20] 白 玉 川 , 万 春 艳 , 黄 本 胜 , 等 . 河 网 非 恒 定 流 数 值 模 拟 的 研 究 进 展 [J]. 水 利 学
报,2000(12):43~47.
[21] DronkersJJ.河流近海岸区和外海的潮汐计算[J].水利水运科技情报,1976 年增刊.
[22] 张二骏,等.河网非恒定流三级联合算法[J].华东水利学院学报,1982,(1).
[23] Schulze K W.Finite Element Analysis of Long Waves in open Channel Systems.Finit
eElement Method in Flow Problems. Chapter5 ,Ed. By Oden J T etal., The Univ. Of
Alabama Press, 1974,PP.295 303.
network[J].
artificial
Journal
of
in
neural
[24] 吴寿红.河网非恒定流四级解算法[J],水利学报,1985,(8).
[25] 李义天.河网非恒定流隐式方程组的汊点分组解法[J],水利学报,1997,(8).
[26]Naidu B J, Murty BhallamudiS, and Narasimhan S, GVF Computation in tree type
channel Networks.[J] J.HydrEngrg,1997,(8):700 708.
[27]TucciCEM,ChenYH.Unsteady water quality model for river network[J]. Proc
ASCE,1981,107(WR2):477-493.
8