水库调洪计算的数值解法及其程序
- 格式:pdf
- 大小:363.53 KB
- 文档页数:6
洪水调节计算方法
洪水调节计算可是个很重要的事儿呢。
有一种简单的方法叫静库容法。
就好比把水库当成一个大盆子,洪水来了就往盆子里装。
这个方法假设水库的水面是水平的,只考虑水库的库容和进出库的水量关系。
计算的时候,就看进来多少水,出去多少水,然后算出水库里水位的变化。
比如说,入库洪水流量是一个数值,出库流量如果是个固定的值,那两者的差值就是水库里水量增加或者减少的速度啦,根据这个就能算出不同时间水库的水位情况呢。
还有动库容法哦。
这就比静库容法复杂一点啦。
因为它考虑到了水库里的水不是像静库容法里那样完全水平的,水面是有坡度的。
就像你看河流里的水,它不是平平整整的。
动库容法在计算的时候,要考虑到这个坡度对库容的影响。
这就需要更多的测量数据啦,像不同位置的水深啊之类的。
另外,还有水量平衡法。
这个就像是在算家里的收支账一样。
把水库看成一个整体,进来的水就是收入,出去的水就是支出。
在一段时间内,入库水量减去出库水量,就等于水库里水量的变化量。
如果入库的洪水忽大忽小,出库流量也跟着调整,那通过这个方法就能很好地算出每个时段水库的水量情况,从而知道水位的变化情况。
不管是哪种方法,都是科学家们经过很多研究才得出来的,它们在防洪减灾方面可起着超级大的作用哦。
C-2 水库调洪演算的数值解程序作者 张校正(新疆水利厅 )一、程序功能已知水库的水位--水面面积关系,洪水量过程线,对于每一种调洪方案(包括泄流条件、调洪方式、泄水建筑物参数)由调洪起始水位依次计算,直至洪水过程结束,计算机输出各时段末之水位、泄洪洞流量、溢洪道流量、水库出库总流量等。
并用彩色曲线绘制洪水过程线、泄洪过程线和水库水位变化线。
二、算法简介1,水库水量平衡分方程的数值解:水库水量平衡微分方程:q Q dt dZ f -=式中: f=f(z) 水库水面面积,是水位z 的函数;Z=Z(t) 水位,是时间t 的函数;Q=Q(t) 入库流量,是时间t 的函数;Q=q(z) 出库流量,是水位z 的函数。
将上式移项,并定义调洪函数)()()(),(z f Z q t Q Z t F -=则得 ⎪⎩⎪⎨⎧==00)(),(Z t Z Z t F dt dZ 这是一个一阶常微分方程的初值问题。
应用定步长的龙格-库塔方法求解。
其公式为:)22(6143211K K K K Z Z n n ++++=-式中: )()()(),(111111------⨯=⨯=n n n n n Z f Z q t Q T Z t F T K)21()2()2()2,2(11111112K Z f K Z q T t Q T K Z T t F T K n n n n n ++-+⨯=++⨯=----- )2()2()2()2,2(212112113K Z f K Z q T t Q T K Z T t F T K n n n n n ++-+⨯=++⨯=----- )()()(),(3131314K Z f K Z q t Q T K Z t F T K n n n n n ++-⨯=+⨯=---T 为洪水流量时段间隔;n=1,2,……,J2,泄流量公式:当泄水建筑物为深孔时,)(2111111A C Z g B A M q --=式中:M 1 流量系数;A 1 泄流孔口高;B 1 泄流孔口宽;Z 水位;C 1 泄流孔口底槛高程。
调洪演算计算说明书一、 相关资料中包水利枢纽工程是三等工程,溢洪道设计洪水标准为五十年一遇(P=2%)至一百年一遇(P=1%),校核洪水标准为千年一遇(P=0.1%).二、基本原理1.泄水建筑物尺寸:溢洪道堰顶高程519m ,采用3孔86m m ⨯(宽⨯高)的弧形门控制。
由2/302q H g m nb ⋅=ε (其中侧收缩系数ε=0.92,n 为所开孔数, 流量系数m=0.48,单孔堰顶宽度b=8m ,g=9.812/m s ,堰顶水头0H =水位Z-堰顶高程,。
不计流速水头。
) 计算出下泄流量2.设计洪水来临时,用左右2孔泄洪;校核洪水来临时,用3孔泄洪。
3.基本计算公式为:()()()t V V q q Q Q ∆-=+-+/2/2/122121式中: Q 1, Q 2--分别为计算时段初、末的入库流量,m 3/s ; v 1,v 2--分别为计算时段初、末水库的蓄水量,m 3 ; q 1,q 2--分别为计算时段初、末的下泄流量,m 3/s ; t ∆--计算时段,一般取1小时。
4.下游安全泄量及起调水位该水利枢纽没有下游防洪要求,一般在洪水来临时,水库将预泄库水至水库防洪限制水位,以便有足够的库容蓄洪或滞洪。
防洪限制水位是水库在汛期允许兴利蓄水的上限水位,则调洪计算从水位525.3m 起调。
5.水库运行方式根据题目分析,本工程采用3孔溢洪道泄洪,设计洪水来临时,用左右2孔泄洪;校核洪水来临时,用3孔泄洪。
在洪水期间洪水来临时,先用闸门控制下泄流量q 并使其等于洪水来水量Q,使水库水位保持在防洪限制水位525.3m不变;当洪水来水量Q继续增大时,闸门逐渐打开;当闸门达到全开后,就不再用闸门控制,下泄流量q随水库水位z 的升高而增大,流态为自由流态,情况与无闸门控制一样。
6.计算方法:先决定开始计算时刻和此时的q1、V1,然后假定下泄流量q2值,再由计算V2值,再查q-V表得出q2’值,水量平衡方程()()()t-+2/2/=+/VV-qqQ∆Q211122比较q2和q2’,若二者基本相等,则假定正确,否则重新试算,直到大致相等为止,依次计算下去。
2.4.2 调洪计算方法水库调洪是在水量平衡和动力平衡的支配下进行的,本次计算单辅助线法计算。
水量平衡的数学表达式为:221Q Q +t ∆ -221q q + t ∆=V 2-V 1式中:Q 1,Q 2——时段初、末入库流量,m 3/s ;q 1,q 2——时段初、末出库流量,m 3/s ;V1,V2——时段初、末水库蓄水量,m 3;t ∆——计算时段,t ∆=1h=3600s 。
将水量平衡方程进行变换得到:)(22)2(1112221q t V q Q Q q t V +∆+-+=+∆ 建立q ~2q t V +∆函数关系曲线,绘出q ~2q t V +∆辅助线,连续求出水库的下泄流量过程。
2.4.3 调洪演算成果按照不同频率入库设计洪水过程线,逐时段查算辅助曲线,确定水库出库流量过程。
根据上述入库设计洪水过程线、库容曲线、起调水位进行调洪演算。
本次调洪演算成果见表2-9。
调洪演算成果表2-92.5 坝顶高程计算水库主坝为浆砌石坝,坝顶超高计算公式采用《砌石坝设计规范》(SL25-2006)中公式进行计算:c z b h h H H ++∆=式中:H ∆——坝顶超高,m ;H b ——波浪高,m ;H z ——风浪中心线至正常蓄水位或校核洪水位的高差,m ;H c ——安全超高,5级坝,设计情况A=0.3m ,校校情况A=0.2m 。
根据当地提供的风速风向资料,水库水面以上10m 高度处,年最大平均风速为16m/s 。
根据《砌石坝设计》(SL25-2006)及《水利水电等级划分及洪水标准》(SL252-2000)有关规定,永久建筑物级别为5级。
根据《砌石坝设计规范》(SL25-2006)波高、波长按官厅公式(C.4.1-1)和(C.4.1-2)计算: )(11.4.)(0076.03/12020121-=C v gD v v gh o b)(21.4.)(33.015/42020157-=C v gD v v gLm o式中:H b ——波高(当2502020-=v gD时,为累积频率5%的波高h s%;当当100025020-=v gD 时,为累积频率10%的波高h 10%),m ;L m ——平均波长,m ;v 0——计算最大风速(设计情况采用多年平均年最大风速的1.5倍 ,校核情况采用多年平均最大风速),m/s ; D ——风区长度,m ;g ——重力加速度,9.81m/s 2。
调洪演算计算说明书一、 相关资料中包水利枢纽工程是三等工程,溢洪道设计洪水标准为五十年一遇(P=2%)至一百年一遇(P=1%),校核洪水标准为千年一遇(P=0.1%).二、基本原理1.泄水建筑物尺寸:溢洪道堰顶高程519m ,采用3孔86m m ⨯(宽⨯高)的弧形门控制。
由2/302q H g m nb ⋅=ε (其中侧收缩系数ε=0.92,n 为所开孔数, 流量系数m=0.48,单孔堰顶宽度b=8m ,g=9.812/m s ,堰顶水头0H =水位Z-堰顶高程,。
不计流速水头。
) 计算出下泄流量2.设计洪水来临时,用左右2孔泄洪;校核洪水来临时,用3孔泄洪。
3.基本计算公式为:()()()t V V q q Q Q ∆-=+-+/2/2/122121式中: Q 1, Q 2--分别为计算时段初、末的入库流量,m 3/s ; v 1,v 2--分别为计算时段初、末水库的蓄水量,m 3 ; q 1,q 2--分别为计算时段初、末的下泄流量,m 3/s ; t ∆--计算时段,一般取1小时。
4.下游安全泄量及起调水位该水利枢纽没有下游防洪要求,一般在洪水来临时,水库将预泄库水至水库防洪限制水位,以便有足够的库容蓄洪或滞洪。
防洪限制水位是水库在汛期允许兴利蓄水的上限水位,则调洪计算从水位525.3m 起调。
5.水库运行方式根据题目分析,本工程采用3孔溢洪道泄洪,设计洪水来临时,用左右2孔泄洪;校核洪水来临时,用3孔泄洪。
在洪水期间洪水来临时,先用闸门控制下泄流量q 并使其等于洪水来水量Q,使水库水位保持在防洪限制水位525.3m不变;当洪水来水量Q继续增大时,闸门逐渐打开;当闸门达到全开后,就不再用闸门控制,下泄流量q随水库水位z 的升高而增大,流态为自由流态,情况与无闸门控制一样。
6.计算方法:先决定开始计算时刻和此时的q1、V1,然后假定下泄流量q2值,再由计算V2值,再查q-V表得出q2’值,水量平衡方程()()()t-+2/2/=+/VV-qqQ∆Q211122比较q2和q2’,若二者基本相等,则假定正确,否则重新试算,直到大致相等为止,依次计算下去。
第三章调洪计算调洪计算目的水库调洪计算的目的是在已拟定泄洪建筑物及已确定防洪限制水位(或其他的起调水位)的条件下,用给出的入库洪水过程、泄洪建筑物的泄洪能力曲线及库容曲线等基本资料,按规定的防洪调度规则,推求水库的泄流过程、水库水位过程及相应的最高调洪水位和最大下泄流量。
调洪演算的原理水库调洪计算的基本公式是水量平衡方程式:t t t t t t V V t q q t Q Q -=∆+-∆++++1121121)()( (3-1)式中t ∆—计算时段长度,s ;1,+t t Q Q —t 时段初、末的入库流量,m 3/s ; 1,+t t q q —t 时段初、末的出库流量,m 3/s ; 1,+t t V V —t 时段初、末水库蓄水量,m 3。
水库泄流方程 :q =f (V ) (3-2)用已知(设计或预报)的入库洪水过程线Q ~t ,由起调水位开始,逐时段连续求解(3-1)和(3-2)组成的方程组,从而求得水库出流过程q ~t ,这就是调洪演算的基本原理。
这里采用单辅助线半图解法,联解(2-1)和(2-2)两个方程,将(3-1)改写为:(V t /△t+q t /2 )+Q-q t = (V t+1/△t)+(q t+1/2 ) (3-3)式中Q —计算时段平均入流量,Q =(Q t + Q t+1)/2;其他同(3-1)也就是说,可以事先绘制q ~(V/△t)+(q/2 )的关系曲线,即调洪演算工作曲线,因式3-3)的左端各项为已知数,故式(3-3)右端项也可求出,然后根据(V t+1/△t)+(q t+1/2 )的值,通过工作曲线q~(V/△t)+(q/2 )可查出q t+1的值。
因第一时段的V2、q2就是第二时段的V1、q1,于是可重复以上步骤连续进行计算,直到求出结果。
调洪计算结果整理调洪演算基本资料水库特征水位:正常蓄水位1856m,汛期限制水位1854m,死水位1852m积石峡入库洪水过程线见下表:表2-1积石峡入库洪水过程线调洪计算过程及结果方案一:1. 拟定泄水建筑物型式、尺寸及堰顶(或底坎)高程:左岸溢洪道: 单孔溢洪道,B=,H=18m,堰顶高程为1833m。
第三章调洪计算3.1调洪计算目的水库调洪计算的目的是在已拟定泄洪建筑物及已确定防洪限制水位(或其他的起调水位)的条件下,用给出的入库洪水过程、泄洪建筑物的泄洪能力曲线及库容曲线等基本资料,按规定的防洪调度规则,推求水库的泄流过程、水库水位过程及相应的最高调洪水位和最大下泄流量。
3.2调洪演算的原理水库调洪计算的基本公式是水量平衡方程式:t t t t t t V V t q q t Q Q -=∆+-∆++++1121121)()( (3-1)式中t ∆—计算时段长度,s ;1,+t t Q Q —t 时段初、末的入库流量,m 3/s ; 1,+t t q q —t 时段初、末的出库流量,m 3/s ; 1,+t t V V —t 时段初、末水库蓄水量,m 3。
水库泄流方程 :q =f (V ) (3-2)用已知(设计或预报)的入库洪水过程线Q ~t ,由起调水位开始,逐时段连续求解(3-1)和(3-2)组成的方程组,从而求得水库出流过程q ~t ,这就是调洪演算的基本原理。
这里采用单辅助线半图解法,联解(2-1)和(2-2)两个方程,将(3-1)改写为:(V t/△t+q t/2 )+Q-q t= (V t+1/△t)+(q t+1/2 ) (3-3)式中Q—计算时段平均入流量,Q=(Q t + Q t+1)/2;其他同(3-1)也就是说,可以事先绘制q~(V/△t)+(q/2 )的关系曲线,即调洪演算工作曲线,因式3-3)的左端各项为已知数,故式(3-3)右端项也可求出,然后根据(V t+1/△t)+(q t+1/2 )的值,通过工作曲线q~(V/△t)+(q/2 )可查出q t+1的值。
因第一时段的V2、q2就是第二时段的V1、q1,于是可重复以上步骤连续进行计算,直到求出结果。
3.3调洪计算结果整理3.3.1调洪演算基本资料水库特征水位:正常蓄水位1856m,汛期限制水位1854m,死水位1852m 积石峡入库洪水过程线见下表:表2-1积石峡入库洪水过程线3.3.2调洪计算过程及结果 方案一:1. 拟定泄水建筑物型式、尺寸及堰顶(或底坎)高程:左岸溢洪道: 单孔溢洪道, B=16.5m,H=18m,堰顶高程为1833m 。