水库调洪半图解法 matlab编程

1、作者：双辅助线法双辅助线法2 水库调洪计算的直接目的，在水库调洪计算的直接目的，在于求出水库逐时段的蓄水、泄水变于求出水库逐时段的蓄水、泄水变化过程，从而获得调节该次洪水后化过程，从而获得调节该次洪水后的水库最高洪水位和最大下泄流量，的水库最高洪水位和最大下泄流量，供进一步防洪计算分析之用。供进一步防洪计算分析之用。3122121VVt2qqt2QQ222211qtVqtVQ等式两边同时除以等式两边同时除以t, t,并并移移项项2tV2tV21qq、均可与水库水位均可与水库水位Z建立函数关系建立函数关系221QQQ 双辅助线法的基本原理仍然是逐时段连续求解水库的水量平衡方程和蓄泄方程，但为了避

2、免试算，需对这两个方程作适当的变换。 因此，可根据选定的计算时段t t，已知的水库容积关系曲线，以及根据水力学公式算出的水位下泄流量关系曲线，事先计算并绘制曲线组：2)(1qtVZf2)(2qtVZfqZf)(3)/(2),/(2),/(333smqtVsmqtVsmQ)(mZ调洪计算半图解法的双辅助线调洪计算半图解法的双辅助线)(1Zf)(2Zf)(3Zf c.根据根据Z2值，利用水库曲线即可求出值，利用水库曲线即可求出V2。 d.将点代表的将点代表的2值作为下一时段的值作为下一时段的1值值a.根据已知的根据已知的Qt过程线、过程线、ZV曲线、曲线、Z限限、计算时、计算时段段t，确定调洪计算

3、的起始时段，并划分各计算，确定调洪计算的起始时段，并划分各计算时段。算出时段。算出平均入库流量平均入库流量Q1以及定出第一时段以及定出第一时段初初始的始的Z1、q1、V1各值。各值。 b.利用辅助线在图上求解得出利用辅助线在图上求解得出Z2。 e.求出该时段的求出该时段的2、q2、V2值。如此逐时段进行值。如此逐时段进行计算，即可得到下泄流量过程线。计算，即可得到下泄流量过程线。ABCABC)2()2(2211qtVqtVQQABCABC(2)在图的水位坐标轴上在图的水位坐标轴上量取第一时段的量取第一时段的Z1，得，得a点点。作水平线。作水平线ac交曲交曲线线A于于b点，并点，并使使bc= 。

4、因曲线因曲线A是（是（V/t-q/2)=f1(z），），a点代表点代表Z1，ab=(V1/t-q1/2)，ac= + (V1/t-q1/2)，按式按式(3-5)，即等于，即等于（V2/t+q2/2）。）。QQ从从c c点作垂线交曲线点作垂线交曲线B B于于d d点，过点，过d d点作水平线点作水平线dede交水位坐标轴于交水位坐标轴于e e，显然显然de=acde=ac= =（V V2 2/t t+ +q q2 2/2/2）。因曲线）。因曲线B B是是 （V V/t t- -q q/2)=/2)=f f2 2(Z(Z），），d d点在曲线点在曲线B B上，上，e e就应代表就应代表Z Z2 2

5、，从，从e e点可读出点可读出Z Z2 2值值。ABCABCABCABC4) de交曲线交曲线C于于f点，过点，过f 点作垂线交点作垂线交q坐标轴于坐标轴于g点。因曲线点。因曲线C是是q=f3(Z),e代表代表Z2,而而ef是是q2，即从，即从g点可以读出点可以读出q2的值。的值。(5) 根据根据Z2值，利用水库水位容积关系曲线可值，利用水库水位容积关系曲线可求出求出V2的值的值。(6) 将将e点所代表的点所代表的Z值作为第二时段的值作为第二时段的Z1，按，按上述方法依次进行图解计算，又可求出第二时上述方法依次进行图解计算，又可求出第二时段的段的Z2、Q2、V2值。按逐时段进行计算，将结值。按

6、逐时段进行计算，将结果列表，即完成全部计算。果列表，即完成全部计算。水库调洪半图解法的实例水库调洪半图解法的实例 1.2 1.2 水库调洪计算的半图解法水库调洪计算的半图解法 解解 根据已知的入库洪水流量过程线、水库水位容积关系曲根据已知的入库洪水流量过程线、水库水位容积关系曲线、汛期防洪限制水位、计算时段等，确定调洪计算的起线、汛期防洪限制水位、计算时段等，确定调洪计算的起始时段始时段t，并划分各计算时段。算出各时段的，并划分各计算时段。算出各时段的平均入库流平均入库流量量 ，并定出第一时段初始的，并定出第一时段初始的Z1、q1、V1值。值。Q（V/t-q/2)=f1(z）和和(V/t+q/