水库调洪演算试算法
一、水库调洪计算的任务
入库洪水流经水库时，水库容积对洪水的拦蓄、滞留作用，以及泄水建筑物
对出库流量的制约或控制作用，将使出库洪水过程产生变形。与入库洪水过程相比，出库洪水的洪峰流量显著减小，洪水过程历时大大延长。这种入库洪水流经水库产生的上述洪水变形，称为水库洪水调节。水库调洪计算的目的是在已拟定泄洪建筑物及已确定防洪限制水位（或其他的起调水位）的条件下，用给出的入库洪水过程、泄洪建筑物的泄洪能力曲线及库容曲线等基本资料，按规定的防洪调度规则，推求水库的泄流过程、水库水位过程及相应的最高调洪水位和最大下泄流量。
若水库不承担下游防洪任务，那么水库调洪计算的任务是研究和选择能确保水工建筑物安全的调洪方式，并配合泄洪建筑物的形式、尺寸和高程的选择，最终确定水库的设计洪水位、校核洪水位、调洪库容及二种情况下相应的最大泄流量。若水库担负下游防洪任务，首先应根据下游防洪保护对象的防洪标准、下游河道安全泄量、坝址至防洪点控制断面之间的区间入流情况，配合泄洪建筑物形式和规模，合理拟定水库的泄流方式，确定水库的防洪库容及其相应的防洪高水位；其次，根据下游防洪对泄洪方式的要求，进一步拟定为保证水工建筑物安全的泄洪方式，经调洪计算，确定水库的设计洪水位与校核洪水位及相应的调洪库容。
二、水库调洪计算基本公式
洪水进入水库后形成的洪水波运动，其水力学性质属于明渠渐变不恒定流。
常用的调洪计算方法，往往忽略库区回水水面比降对蓄水容积的影响，只按水平面的近似情况考虑水库的蓄水容积（即静库容）。水库调洪计算的基本公式是水量平衡方程式：
2(Qt1?Qt?1)?t?12(qt?qt?1)?t?Vt?1?Vt （3-1）
式中： ?t——计算时段长度（s）；
Qt,Qt?1——t时段初、末的入库流量（m3/s）；
qt,qt?1——t时段初、末的出库流量（m/s）；
3
Vt,Vt?1——t时段初、末水库蓄水量（m）。
3
当已知水库入库洪水过程线时，Qt,Qt?1均为已知；Vt,qt则是计算时段t开始的初始条件。于是，式中仅Vt?1,qt?1为未知数。必须配合水库泄流方程q=f（V）与上式联立求解Vt?1,qt?1的值。当水库同时为兴利用水而泄放流量时，水库泄流量应计入这部分兴利泄流量。假设暂不计及自水库取水的兴利部门泄向下游的流量，若泄洪建筑物为无闸门表面溢洪道，则下泄流量q的计算公式为: q1??mBh12gh1 （3-2）
式中： ε 侧收缩系数； m 流量系数； B 溢洪道宽； h1 堰上水头。
若为孔口出流，则泄流公式为：
q2???2gh2 （3-3）
式中： μ 孔口出流系数；
ω 孔口出流面积； h2 孔口中心水头。
由式（3-2）或（3-3）所反映泄流量q与泄洪建筑物水头h的函数关系可转换为泄流量q与库水位Z的关系曲线q ＝f（Z）。借助于水库容积特性V ＝f（Z）,
可进一步求出水库下泄流量q与蓄水容积V的关系，即
q＝f（V） （3-4）
说明如何进行一次洪水的水库调洪计算。图中Q～t为入库洪水过程线；q～t为水库调洪计算需要推求的出库流量过程线。设?t为计算过程的面临时段，由入库洪水资料可知时段初、末的流量Qt,Qt?1的数值，Vt，qt为该时段已知的初始条件。图中阴影线的面积表示该时段水库蓄水量的增量ΔV，即ΔV=Vt+1-Vt。利用式（3-1），（3-4）可求解时段末的水库蓄水量Vt+1和相应的出库流量qt?1。前一个时段的Vt?1,qt?1求出后，其值即成为后一时段的Vt,qt值，使计算有可能逐时段地连续进行下去。 必须指出，上述水库调洪计算中采用的泄流函数式
q =f（v）是基于泄洪设施为自由溢流的条件建立的。所谓自由溢流是指泄洪设施不设闸门，或虽设有闸门，但闸门达到的开度不对水流形成制约的情况。
三、水库调洪演算试算法
水库调洪演算就是联解式(3-1)和(3-4) 。常用的算法有试算法(迭代法)和图解法。试算法可达到对计算结果高精度的要求，但以往靠人工计算时，此法计算工作量大；图解法是为了避免繁琐的试算工作而发展起来的，它实用于人工操作，可大大减轻试算法的人工计算工作量。随着计算机科学技术的迅速发展，上述水库调洪计算的试算法很适合编制电算程序，即在计算机上进行迭代计算，不必再提倡采用图解法来完成调洪计算。在进行迭代计算时,可先假定计算时段末的出库流量的qt?1值, 求出式中待定的时段末水库蓄水量Vt?1的值；也可先假定Vt?1的值, 求出式中待定的qt?1值。最后，在迭代过程中算出满足精度的解。下述迭代算法(以先假定qt?1的值为例)的步骤可以作为编制水库调洪计算软件的程序流程。
（1）初步假设计算时段末的出库流量qt?1的值，代入式(3-1)，可初步求出