一维水量水质模型.docx

一维水量水质模型.docx

《一维水量水质模型.docx》由会员分享，可在线阅读，更多相关《一维水量水质模型.docx（17页珍藏版）》请在冰豆网上搜索。

一维水量水质模型

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

对于中小型河流，通常其宽度及水深相对于长度数量较小，扩散质（污染物质、热量）很容易在垂向及横向上达到均匀混合，即扩散质浓度在断面上基本达到均匀状态。

这种情况下，我们只需要知道扩散质在断面的平均分配状况，就可以把握整个河道的扩散质空间分布特征，这是我们可以采用一维圣维南方程描述河流水动力特征或水量特征（水位、流量、槽蓄量等）；用一维纵向分散方程描述扩散质在时间及河流纵向上的变化状况。

特别地，对于稳态水流，可以采用常规水动力学方法推算水位、断面平均流速的沿程变化；采用分段解析解法计算扩散质浓度沿纵向的变化特征。

但是，在非稳态情况下（水流随时间变化或扩散质源强随时间变化）解析解法将无能为力（水流非恒定）或十分繁琐（水流稳态、源强非恒定），这时通常采用数值解法求解河道水量、水质的时间、空间分布。

在模拟方法上，无论是单一河道还是由众多单一河道构成的河网，若采用空间一维手段求解，描述水流、水质空间分布规律的控制方程是相同的，只不过在具体求解方法上有所差异而已。

7.1单一河道的控制方程

7.1.1水量控制方程

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

（1）

（2）

式中t为时间坐标，x为空间坐标，Q为断面流量，Z为断面平均水位，u为断面平均流速，n为河段的糙率，A为过流断面面积，BW为水面宽度（包括主流宽度及仅起调蓄作用的附加宽度），R为水力半径，q为旁侧入流流量（单位河长上旁侧入流场）。

此方程组属于二元一阶双曲型拟线性方程组，对于非恒定问题，现阶段尚无法直接求出其解析解，通常用有限差分法或其它数学离散方法求其数值解。

在水流稳态、棱柱形河道条件下，上述控制方程组退化为水力学的谢才公式，可采用相应的方法求解水流特征。

7.1.2扩散质输运控制方程

描述河道扩散物质运动及浓度变化规律的控制方程为：

带源的一维对流分散（弥散）方程，形式如下：

（3）

式中，C为污染物质的断面平均浓度，Q为流量，为纵向分散系数，S为单位时间、单位河长上的污染物质排放量，K为污染物降解系数，Sr为河床底泥释放污染物的速率。

此方程属于一元二阶偏微分方程，对于非恒定水流问题，微分方程位变系数的偏微分方程，现阶段尚无法直接求出其解析解，通常用有限差分法或其它数学离散方法求其数值解。

在水流稳态、污染源源强恒定条件下，可按水动力特征将河道分为若干子段，在每个分段上，上述控制方程简化为常系数的常微分方程，可采用解析方法秋初起理论解。

7.2单一河道一维水量水质模型

7.2.1单一河道一维水量模型

（1）控制方程的离散

采用四点隐式差分格式离散方程组。

如图1所示，河道被（n+1）个断面分为n个子河段，在第i个子河段M（i,i+1）上，对任一变量取：

（4）

（5）

（6）

图1计算断面示意图

式中，上角标表示时间坐标，下脚标表示空间坐标。

为空间差商的权重系数（）,=0时，此格式为显式格式，而当时，此格式具有隐式差分的特征。

为使差分方程保持无条件稳定，必须。

采用下式进行阻力项的线性化：

（7）

将式（4）-（6）代入连续方程得第i个子河段的差分方程：

（8）

式中，

，

下角标i+1/2表示断面i与断面i+1河段的均值。

按照同样的方法，可得动量方程的差分方程：

（9）

式中，

对任一河段

，可得到方程组：

（10）

对每一河段可列出两个线性代数方程，再加上上下游边界条件，构成完备的封闭方程组，采用追赶法可求得各个断面的水位流量。

（2）边界条件

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

·上游水位边界条件

；下游水位（或流量）边界条件

（或

），追赶形式为：

（11）

式中，

为已知系数，依据上述方程组，可逐步由下边界水位或者流量，推算得到上游各个断面水位流量值。

·上游流量边界条件

；下游水位边界条件

，追赶形式为：

（12）

式中，

为已知系数，依据上述方程组，可逐步由下边界水位，推算得到上游各个断面水位流量值。

7.2.2单一河道一维水质模型

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

对方程（3）进行离散，空间差分采用隐式迎风差分格式。

顺流时（从断面i流向i+1）有：

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

（13）

式中，

为系数，分别表示为：

方程（13）两边同时除以

得到：

（14）

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

（15）

对首断面给定第一类边界条件，对末断面给定第二类边界条件，可得到如下封闭的方程组：

（16）

对方程组（16）采用追赶法可容易求得

等断面的扩散质的浓度。

（2）参数确定

·纵向分散系数EX的确定

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

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

·降解系数K的确定

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

7.2.2应用实例［］

三峡大坝位于县三斗坪中堡岛,葛洲坝位于南津关下游的市境.两坝间水域处于鄂西山区向平原的过渡地带,周围地形地貌呈西北高东南低之势.三峡大坝至葛洲坝河段长38Km,两坝间江面宽210m至1500m,大部分处于西陵峡谷中。

该实例建立了两坝间水量水质模型,分别运用大坝一期围堰及二期围堰施工期间的同步水文水质实测资料对模型进行了率定和验证,取得了较好效果.此模型可用以预测大坝施工期间及投入运行后两坝间水流及水质特性的变化.

（１）水量模型率定

利用1996年三斗坪、白庙子及黄陵庙等断面水文观测资料率定糙率，得到各子河段的糙率，率定结果显示糙率分布取值围为0.036至0.050。

1996年白庙子及黄陵庙等断面实测水位过程线及计算水位过程线图略。

（２）水量模型验证

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

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

（３）CODMn模型验证

高锰酸盐指数的降解系数由实测资料，取经验值，不采用模型率定。

采用1998年7~12月的CODMn浓度监测资料对CODMn模型参数可靠性进行验证，模拟时CODMn污染源分点源和面源两种情况。

部份断面CODMn水质因子实测及计算浓度值见图3。

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

7.3平原河网水量水质模型

河网地区是中国社会、经济、文化的发达地区，在国家经济发展中占有举足轻重的地位。

随着地区经济的进一步发展、居民物质精神生活水准的进一步提高，水资源问题日益突出。

如中国著名的长江三角洲、珠江三角洲地区，已成为水质型缺水区，水资源、水环境问题已成为制约经济社会发展、事关区域可持续发展的重要因素。

生产实践的迫切需要使得河网地区的水环境保护研究呈不断深入、系统之势。

作为基础，河网水流水质模拟方法是进行区域环境规划、环境管理等的必备工具，在水环境问题研究中占有重要地位。

因此，该领域的研究一直是环境科学研究人员、环境管理决策部门十分关注的重要问题。

在物理模型、数学模型这两大手段中，因为数学模型具有经济、快捷、实用等优点，加之水网地区河道密布等客观条件的限制，现阶段只能采用数值方法模拟水网地区的水流运动及污染物输运规律。

但由于理论、技术及各种客观条件的制约，目前水流特别是水质数值模拟精度不是十分令人满意。

7.3.1平原河网水量模型

按河网水流的控制方程及对河网的概化处理方式不同，河网地区水流数值模拟方法可分为两大类：

第一类为常用的一维圣维南方程组数值解法；第二类为所谓的“组合单元解法”。

，其中，一维圣维南方程组数值解法又可分为直接解法和间接解法两种。

在直接解法中，较有代表性的有文献[G.Noseda,MathematicalModelofUnsteadyFlowinOpenChannelsNetworks,ProceedingsoftheInternationalSymposiumonUnsteadyFlowinOpenChannels,1976]，[J.J.R.Williams,T.R.E.Chidley,NonliearAnalysisofUnsteadyFlowinOpenChannelNetworks,ProceedingsoftheInternationalSymposiumonUnsteadyFlowinOpenChannels,1976]提出的方法，该类方法将计算断面交替取为水位和流量断面，对河网在所有计算断面上统一对一维圣维南方程组差分离散并求解。

但该方法未知数数量较多，在河网规模较大的情况下，因为河道的交叉衔接，形成的矩阵是一个不规则、不对称的大型稀疏矩阵。

为减少存贮，J.J.Dronkers于1976年提出间接解法[5]的思想，以后又有许多学者对其作了进一步完善。

为提高计算效率，需缩小矩阵规模，大学数力系1976年提出了河网非恒定流隐式方程组稀疏矩阵解法[6岳生等，河网不恒定流隐式方程组稀疏矩阵解法，大学学报（自然科学版），1979年3月]，该方法从河网矩阵本身的特点出发，能够有效地节省存储并提高计算速度，但矩阵中需包含所有断面的未知数，方程规模仍然较大，其实际使用也受到限制。

而间接解法是将断面未知数往交汊点集中，待求出交汊点未知数后，再求解各单一河道未知变量，计算效率较高。

间接解法的思想首先由荷兰水力学专家Dronkers于1976年提出，以后又有许多学者对此方法进行进一步的完善，相继提出了河网非恒定流的二级解法[徐正凡，明渠非恒定流[M].水利电力学院，1983年.]、三级解法[8二骏，河网非恒定流的三级联合解法，华东水利学院学报，1982年第1卷：

1~12]和四级解法[9吴寿红，河网非恒定流四级解法，水利学报，1985年，第8期：

42~50]。

此类间接解法中，以三级解法最为常用[][][]。

其基本求解思路可概括为“单一河道—交汊点—单一河道”，即

（1）先将单一河道划分为若干子河段，在计算断面上对一维Saint-Venant方程组进行有限差分运算，得到各单一河道差分方程组，进行消元计算，得到单一河道首、末断面间流量与水位的相互关系；

（2）根据河道交汊点水量守恒方程，得到并求解交汊点水位方程组，得到所有交汊点水位；（3）根据交汊点水位返回各单一河道，求得各计算断面水位、流量值。

“组合单元解法”[10]由法国水力学专家JeanA.Cunge于1975年首次提出，国也有研究者采用此方法进行了水网地区的水力模拟[11][12]。

此方法的基本思想是：

将河网地区水力特性相似、水位变幅不大的水体概化成单元。

取单元中心的水位为代表水位，采用谢才公式模拟单元间流量交换，根据水量守恒建立每一单元的微分形式的水量守恒方程，离散并得到以单元水位为自变量的代数方程，辅以边界条件，可求得各单元水位、单元间流量。

在上述两类方法中，组合单元解法对河道进行了简单概化，以单元为计算单位，计算相对简单，但模拟精度相对较低，仅实用于大尺度水域的水力模拟。

而Saint-Venat方程组数值解法可以精确计算每一条河道的水流状况，所以成为目前河网水力模拟的主流方法，其中又以三级解法最为常用[13][14][15][16]。

当采用三级解法模拟长系列大围河网水流特性时，需对河道进行概化处理，即将等级较小的河道概化为一条“概化”河道，要