水库调洪计算试算法

水库调洪计算试算法
水库调洪演算试算法一、水库调洪计算的任务
入库洪水流经水库时，水库容积对洪水的拦蓄、滞留作用，以及泄水建筑物对出库流量的制约或控制作用，将使出库洪水过程产生变形。

与入库洪水过程相比，出库洪水的洪峰流量显著减小，洪水过程历时大大延长。

这种入库洪水流经水库产生的上述洪水变形，称为水库洪水调节。

水库调洪计算的目的是在已拟定泄洪建筑物及已确定防洪限制水位(或其他的起调水位)的条件下，用给出的入库洪水过程、泄洪建筑物的泄洪能力曲线及库容曲线等基本资料，按规定的防洪调度规则，推求水库的泄流过程、水库水位过程及相应的最高调洪水位和最大下泄流量。

若水库不承担下游防洪任务，那么水库调洪计算的任务是研究和选择能确保水工建筑物安全的调洪方式，并配合泄洪建筑物的形式、尺寸和高程的选择，最终确定水库的设计洪水位、校核洪水位、调洪库容及二种情况下相应的最大泄流量。

若水库担负下游防洪任务，首先应根据下游防洪保护对象的防洪标准、下游河道安全泄量、坝址至防洪点控制断面之间的区间入流情况，配合泄洪建筑物形式和规模，合理拟定水库的泄流方式，确定水库的防洪库容及其相应的防洪高水位;其次，根据下游防洪对泄洪方式的要求，进一步拟定为保证水工建筑物安全的泄洪方式，经调洪计算，确定水库的设计洪水位与校核洪水位及相应的调洪库容。

二、水库调洪计算基本公式
洪水进入水库后形成的洪水波运动，其水力学性质属于明渠渐变不恒定流。

常用的调洪计算方法，往往忽略库区回水水面比降对蓄水容积的影响，只按水平面的近似情况考虑水库的蓄水容积(即静库容)。

水库调洪计算的基本公式是水量平衡方程式:
11(Q，Q),t,(q，q),t,V,V (3-1) tt，1tt，1t，1t22
,t式中: ——计算时段长度(s);
3Q,Q ——t时段初、末的入库流量(m/s); tt，1
3q,q ——t时段初、末的出库流量(m/s); tt，1
3V,V ——t时段初、末水库蓄水量(m)。

tt，1
当已知水库入库洪水过程线时，均为已知;则是计算时段t开Q,QV,qtt，1tt 始的初始条件。

于是，式中仅为未知数。

必须配合水库泄流方程q=f(V)V,qt，1t，1
与上式联立求解的值。

当水库同时为兴利用水而泄放流量时，水库泄流V,qt，1t，1
量应计入这部分兴利泄流量。

假设暂不计及自水库取水的兴利部门泄向下游的流量，若泄洪建筑物为无闸门表面溢洪道，则下泄流量q的计算公式为: q,,mBh2gh (3-2) 111
式中: ε 侧收缩系数;
m 流量系数;
B 溢洪道宽;
h堰上水头。

1
若为孔口出流，则泄流公式为:
q,,,2gh (3-3) 22
式中: μ 孔口出流系数;
ω 孔口出流面积;
h 孔口中心水头。

2
由式(3-2)或(3-3)所反映泄流量q与泄洪建筑物水头h的函数关系可转换为泄流量q与库水位Z的关系曲线q ,f(Z)。

借助于水库容积特性V ,f(Z), 可进一步求出水库下泄流量q与蓄水容积V的关系，即
q,f(V) (3-4)
说明如何进行一次洪水的水库调洪计算。

图中Q,t为入库洪水过程线;
,tq,t为水库调洪计算需要推求的出库流量过程线。

设为计算过程的面临时段，由入库洪水资料可知时段初、末的流量的数值，V，q为该时段已知的初
Q,Qtttt，1
始条件。

图中阴影线的面积表示该时段水库蓄水量的增量ΔV，即ΔV=V-V。

t+1t
q利用式(3-1)，(3-4)可求解时段末的水库蓄水量V和相应的出库流量。

t+1t，1前一个时段的求出后，其值即成为后一时段的值，使计算有可能逐
V,qV,qt，1t，1tt
时段地连续进行下去。

必须指出，上述水库调洪计算中采用的泄流函数式q =f(v)是基于泄洪设施为自由溢流的条件建立的。

所谓自由溢流是指泄洪设施不设闸门，或虽设有闸门，但闸门达到的开度不对水流形成制约的情况。

三、水库调洪演算试算法
水库调洪演算就是联解式(3-1)和(3-4) 。

常用的算法有试算法(迭代法)和图解法。

试算法可达到对计算结果高精度的要求，但以往靠人工计算时，此法计算工作量大;图解法是为了避免繁琐的试算工作而发展起来的，它实用于人工操作，可大大减轻试算法的人工计算工作量。

随着计算机科学技术的迅速发展，上述水库调洪计算的试算法很适合编制电算程序，即在计算机上进行迭代计算，不必再提倡采用图解法来完成调洪计算。

在进行迭代计算时,可先假定计
qV算时段末的出库流量的值, 求出式中待定的时段末水库蓄水量的值;也t，1t，1
qV可先假定的值, 求出式中待定的值。

最后，在迭代过程中算出满足精度t，1t，1
q的解。

下述迭代算法(以先假定的值为例)的步骤可以作为编制水库调洪计
t，1
算软件的程序流程。

q (1)初步假设计算时段末的出库流量的值，代入式(3-1)，可初步求出t，1 式中待定的时段末水库蓄水量的值。

Vt，1
(2)利用关系，用初求的值，按插值法求出对应的出库流量q。

Vq,f(V)t，1
q (3)检验步骤(1)所假设的时段末的出库流量步骤(2)得到的出库t，1
q,q,流量q的相符合情况。

若设定的允许误差为，，则满足计算精度,,t，1 q要求，结束该时段计算，时段末出库流量及水库蓄水量即为计算的结果。

Vt，1t，1
qq否则，重新假设=(+q)/2，返回步骤(1)进行下一轮迭代计算。

t，1t，1 以上仅以某一计算时段为例，说明水库调洪计算的原理和方法。

对于一场入库洪水的调洪计算，必须从洪水起涨开始，依时序逐时段进行。

第一个计算时段(t,1)可将起调水位(规划设计中对一定设计标准的洪水的调洪计算一般采用防洪限制水位作为起调水位)及其相应的泄水建筑物的泄流能力作为计算初始条件，即已知该时段初的出库流量q和水库蓄水量V，通过调洪计算求11
出时段末的出库流量q和水库蓄水量V。

接着进行第二时段的调洪计算，此时22
q，V已成为第二时段的初始条件，可按同样的方法进行此时段的调洪计算。

22 循此执行逐时段调洪计算，直到水库水位消落至防洪限制水位(或根据要求只
推算到出现水库最高调洪水位)。

现将具体的演算过程用一例子加以说明。

某水库的泄洪建筑物形式和尺寸已定，溢洪堰设有闸门控制。

水库的运行方式是在洪水来临时，先用闸门控制，使水库泄流量等于入库流量，水库保持汛期防洪
限制水位(38m)不变。

随着入库流量继续增大，闸门逐渐开启直至达到全部开启，水库泄流q随库水位的升高而加大，闸门全部开启后的流态为自由泄流。

已知堰顶高程为36m，水库容积曲线V =f(Z)，并根据泄洪建筑物形式和尺
寸，算出水位和下泄量关系曲线q =f(Z)，见表3-4。

计算过程见表3-5。

并按下
列步骤计算。

33表3-4 单位:水位Z(m)，库容V(万m)，泄流量q(m/s) Z 36. 36.5 37
37.5 38 38.5 39 39.5 40 40.5 V 4 330 4 800 5 310 5 860 6 450 7 080 7 760 8 540 9 420 10 250 q 0 22.5 55 105 173.9 267.2 378.3 501.9 638.9 786.1
(1) 将已知入库洪水流量过程线列入表中的第(1)、(2)栏，取计算时段Δ
t =3h =10 800 s;起始库水位为Z=38.0 m，在图中可查出闸门全开时限
3相应的q =173.9 m/s。

3 (2) 在第18小时以前，入库流量Q均小于173.9 m/s，水库按q =Q泄流。

3水库不蓄水，无需进行调洪计算。

从第18小时起，Q开始大于173.9 m/s，
以
3，而初始的第18小时为开始调洪计算的时刻，此时初始的q即为173.9 m/s1 3V为6 450万m。

然后，按水量平衡方程进行计算，将计算结果列入表3-5中相1
应时段的各栏，并点绘在图3-4上。

(3)由表3-5可见，在第36小时，水库水位Z=40.51m、水库蓄水量 V =10 232
333万m、Q =900 m/s、q =781 m/s;而在第39 小时，Z =40.51 m，V =10280
万333m，Q =760 m/s，q =790 m/s。

按前述水库调洪原理，当q 出现时，一定max是q =Q ，此时Z，V均达最大值。

显然，q 将出现在第36小时与第39小时max
之间，在表中并末算出。

通过进一步试算，在第38小时16分钟处，可
33得出q = Q =795 m/s，Z= 40.51 m，V=10 290万m。

max max max
了解以上试算过程后，如果借助计算机将会很快得出计算结果。

必须注意到前面介绍水库调洪计算时，采用了泄水建筑物泄流能力曲线来反映水库出流量与水库蓄水量的函数关系，即q,f(V)。

工程实践中，对于存在闸门开度控制较复杂的调洪情况，可以根据防洪要求，从拟定水库泄洪方式入手，研究确定一种合理的开闸程序，包括启用闸门和变动开度的操作过程，以实现所拟定泄洪方式的泄流过程。

表3-5 调洪计算列表试算法
时间入库洪水时段平均下泄流量时段平均下时段内水水库存水库水位
泄流量qt(h) 流量Q 入库流量 q 库存水量水量 V Z
3333 Q(m/s) /s) (m/s) 变化ΔV ) (m) (m(万m
33(m/s) (万m)
(1) (2) (3) (4) (5) (6) (7) (8) 18 174 257 173.9 180.5 83 6450 38.0 21 340 187 6533 38.1
595 224.5 400
24 850 262 6933 38.4
1385 343.5 1125
27 1920 425 8058 39.2
1685 522.5 1256
30 1450 620 9314 39.9
1280 677.0 651
33 1110 734 9965 40.3
1005 757.5 267
36 900 781 10232 40.5
830 785.5 48
39 760 790 10280 40.51
685 781.0 -104
42 610 772 10176 40.4
535 751.5 -234
45 460 731 9942 40.3
410 702.5 -316
48 360 674 9626 40.1
325 645.5 -346
51 290 617 9280 39.9
qQ注: 表中数字下有横线者为初始已知值; ΔV=( -)Δt。