C2水库调洪演算的数值解程序.docVIP

Ｃ－２ 水库调洪演算的数值解程序 作者 张校正（新疆水利厅 ） 一、程序功能 已知水库的水位--水面面积关系，洪水量过程线，对于每一种调洪方案（包括泄流条件、调洪方式、泄水建筑物参数）由调洪起始水位依次计算，直至洪水过程结束，计算机输出各时段末之水位、泄洪洞流量、溢洪道流量、水库出库总流量等。并用彩色曲线绘制洪水过程线、泄洪过程线和水库水位变化线。 二、算法简介 1，水库水量平衡分方程的数值解： 水库水量平衡微分方程： 式中: f=f(z) 水库水面面积，是水位z的函数； Z=Z(t) 水位，是时间t的函数； Q=Q(t) 入库流量，是时间t的函数； Q=q(z) 出库流量，是水位z的函数。 将上式移项，并定义调洪函数 则得 这是一个一阶常微分方程的初值问题。应用定步长的龙格－库塔方法求解。其公式为： 式中： T 为洪水流量时段间隔； n=1,2,……,J 2，泄流量公式： 当泄水建筑物为深孔时， 式中：M1 流量系数； A1 泄流孔口高； B1 泄流孔口宽； Z 水位； C1 泄流孔口底槛高程。 当泄水建筑物为堰流时： 式中：M2 流量系数； B2 堰顶宽度； C2 堰顶高程。 3，水位－ 水面面积关系，以数值方式输入，程序中用线性插值以水位求面积值。 三、程序性能 1，可以设置10个泄洪洞,10个溢洪道。 2，可以设置变宽变高的泄流孔(随水位不同,孔宽孔高可以不同)。 3，可以设置变宽的溢洪道(随水位不同,溢流宽度可以不同)。 4，可以用其它泄流方式,即用水位～泄流量曲线调洪。 5，以上各种泄流方式可以兼而有之,也可只有部分。 6，设置调洪起始水位，设置水库防洪下限水位，设置下游安全限制泄量，设置电站、灌溉等常流量。 四、使用方法： 1，输入数据准备： 数据按下列顺序: 工程名C$, 洪水频率D$, 水库水位～水面面积关系曲线结点数 K， K个 水库水位ZP(米)～水面面积F(万平方米)， 洪水过程时段数 J，时段间隔 T(秒)， 洪水过程 QL(0～J)(立方米/秒) 防洪下限水位 ZM (米)，调洪起始水位 Z0 (米),泄洪起始流量 Q(0) (立方米/秒)，电站常流量 W (立方米/秒)，下游安全限制泄量 QA (立方米/秒) 泄洪洞个数GO(若泄洪洞个数为0,以下泄洪洞参数不要) 第一个泄洪洞:进口宽 B1(米), 进口高 A1(米),进口高程 C1(米), 流量系数 M11， 闸门宽B11(米), 闸门高A11(米),闸门底高程 C11(米), 限泄流量QX(立方米/秒)， 第二个泄洪洞:进口宽 B1(米), 进口高 A1(米),进口高程 C1(米), 流量系数 M11， 闸门宽B11(米), 闸门高A11(米),闸门底高程 C11(米), 限泄流量QX(立方米/秒)， ...... 溢洪道个数HO(若溢洪道个数为0,以下溢洪道参数不要) 第一个溢洪道:流量系数 M2, 底宽 B2(米), 底高程 C2(米) 第二个溢洪道:流量系数 M2, 底宽 B2(米), 底高程 C2(米) ...... 变宽变高的泄流孔不同状况个数KK1(若KK1为0,以下参数不要) 第一种情况:适用的下限水位ZZ1(米), 适用的上限水位ZZ2(米), 闸门总宽度BB(米), 闸门开启高度AA(米), 流量系数MM, 闸底高程VV(米)， 第二种情况:适用的下限水位ZZ1(米), 适用的上限水位ZZ2(米), 闸门总宽度BB(米), 闸门开启高度AA(米), 流量系数MM, 闸底高程VV(米)， 变宽的溢洪道不同状况个数KK0(若KK0为0,以下参数不要)， 第一种情况:适用的下限水位ZZ3(米), 适用的上限水位ZZ4(米), 流量系数MM, 闸孔宽度DD(米), 闸底高程EE(米)， 第二种情况:适用的下限水位ZZ3(米), 适用的上限水