水科学进展Vol-应用基础与工程科学学报.DOC

流域水量调控非恒定流数值计算模型( 冶运涛1,2，蒋云钟1，雷晓辉1，假冬冬2，尚毅梓3 （） 摘要：。 关键词：中图分类号：TV 122 文献标志码：A 文章编号：1001-6791 就数值计算模型建立而言，水工建筑物的修建改变了水体运动动量守恒特性，导致了不能直接搬用圣维南方程进行计算，而且给数值模拟带来了难度，其原因有二：一是调控使水工建筑物过流能力模拟计算与上下游河道的过流能力不一致，易造成模拟计算的不稳定；二是过流计算格式与上下游河道水流计算格式不一致，易造成计算震荡。因此，需要探讨描述水工建筑物过流能力和调度方程的计算格式。 .1 基本方程 用圣维南方程组描述渐变性一维非恒定明渠水流，包括连续方程（质量守恒方程）和运动方程（动量守恒方程）。河网中某一河道水流控制方程如下： 连续方程： （1） 运动方程： （2） 式中，Q为流量；A为过流断面面积；q表示侧向入流；(为断面动量修正系数；g为重力加速度；Z为水位；Sf为阻力坡降；x为沿着河道纵向距离；t为时间坐标。 摩擦阻力坡降Sf利用曼宁公式计算： （3） 式中，n为曼宁糙率系数；R为水力半径（过流面积与湿周比值）。 1.2 边界条件 对河网中存在入流或出流的外河道的开边界，一般为缓流，每个开边界处提供1个边界条件，可以是流量过程、水位过程、流量-水位关系中任意一种。采用普林斯满格式计算急流或者跨临界流时，则需要修正离散格式近似处理[3]。对单独的河道内，可以采用圣维南方程组描述；但是对于河道交界处的汊点或者存在水工建筑物处，需要等价于圣维南方程组的两个方程描述，即质量守恒方程和能量守恒方程。假定汊点处水量不发生变化，连续方程写为： （4） 假定忽略汊点处的损失和流速水头差，汊点处的能量守恒方程可以近似为： （5） 式中，Zi和Zo分别表示入流水位和出流水位。 1.3 离散求解 采用Preissmann四点隐式差分格式离散控制方程 （6） 如果河网中节点数目为N个，未知数为2N个，离散圣维南方程组得到2(N-1)个代数方程，由边界条件补充2个代数方程，共形成2N个非线性代数方程组。 采用三级解法求解，该算法将河网分成微段、河段、节点三级，逐级处理，再联合运算。 经过离散后，消元处理，忽略时间上标，得到如下差分方程。 （7） 式中：αj、βj、γj、θj、ηj、ξj （8） 2 特殊问题处理 2.1 河网初始水面和流量的赋值 不合理的初始自由水面和流量赋值可能不利于数值计算稳定性。通常采用以下三种方式进行初始水面和流量赋值：①在水文站点密集情况下，利用实测数据并进行水面插值方式确定初始自由水面和流量；②按河道大小分配流量，采用恒定流模型给出自由水面分布；③按河道过水面积大小确定流量，初始自由水面则采用较高水平水位，而后逐渐降低水位在河段末端达到实测值，经过迭代计算就能使初值分布较为合理。本文采用第三种方法，自由水面和流量的分配误差会随着迭代最终趋于合理的值，适合于复杂河网非恒定流的数值模拟。 2.2 河道的干湿边界处理 水位涨落、河道地形冲淤变化以及闸坝调控等情形下，部分河道就会出现干湿交替情况，会影响数值计算的稳定性。本文采用河道断面窄缝法解决干湿交替的问题，该方法就是在河道每个断面中挖出宽度约为1mm左右和深度约为20m的窄深“缝隙”，当实际河道状态由湿变干时，有微量水流从中通过，对计算结果影响较小，能够满足实际要求。使用该方法简化了干湿交替计算的具体过程，保证计算的稳定性。。 3 水量调控计算模式 3.1 闸堰调控计算模式 闸门是按照某种规则实行闸门开度的变化而调整自然水系或人工水系中水位或流量变化。由于天然来水和用水需求的不确定性，闸门调控运行有完全关闭或开启，或者处于一定开度等多种工况，这就导致闸下过流可能出现闸孔自由出流、闸孔淹没出流、堰流自由出流、堰流淹没出流等多种流态,且各种流态都属于急变流，给数值模拟带来了难度，因而要求数值计算模型能够描述各个流态以及流态之间转化过程的复杂关系。 闸孔出流与堰流是密切相关的，在一定边界条件驱动下，根据闸门开度与上下游水位差之间的关系判定是属于闸孔出流还是属于堰流；进而判定闸孔出流是属于淹没出流还是自由出流，堰流是属于自由出流还是淹没出流。在数值计算时，首先要判定闸门是否处于关闭状态，若处于关闭状态，则需要特殊处理，否则就根据上述思路进行计算。将堰所在的河段作为单独的计算河段，为了便于融合到整个河网计算模式中，就要将过堰流量非线性表达转化为上下游水位的线性化关系。 当闸门关闭时，过流量恒为零，此时，双向追赶系数均等于零。 当闸门开启时，假定闸底坎宽顶堰，闸门开度为G0，闸门前堰上水头为H，当G0/H≤0.65时，按闸孔出流公式计算，当G0/H＞0.65时，按堰流出流