四川省水库洪水计算和调洪演算
小（一）型水库洪水计算程序：根据四川省一九七九年版《水文手册》方
册》方法计算洪峰流量和洪水总量。
小（一）型水库洪水计算程序：根
序：根据四川省一九七九年版《水文手册》方法计算洪峰流量和洪水总量。
算，不要改动；兰色格内数据可不改动；绿色格内可不输入数据；紫色格内根据暴雨历时确定是否输入数据；白色格内要输洪 峰 流 量 计 算
小（一）型水库洪水计算程序：根据四川省一九七九年版《水文手册》方法
》方法计算洪峰流量和洪水总量。
小时时
算
小（一）型水库洪水计
值
洪水总量计算
洪水计算程序：根据四川省一九七九年版《水文手册》方法计算洪峰流量和洪水总量。
算机计算，不要改动；兰色格内数据可不改动；绿色格内可不输入数据；紫色格内根据暴雨历时确定是否输入数据；白色格内
分区计算小（一）型水库洪水计算程序：根据四川省一九七九年版《水文手册》方
调 洪 演 算
演 算不输入数据；紫色格内根据暴雨历时确定是否输入数据；白色格内要输入所需数据。
册》方法计算洪峰流量和洪水总量。
小（一）型水库洪水计算
说明：表中黄色格
调 洪
调 洪 演 算
水计算程序：根据四川省一九七九年版《水文手册》方法计算洪峰流量和洪水总量。
黄色格内数据由计算机计算，不要改动；兰色格内数据可不改动；绿色格内可不输入数据；紫色格内根据暴雨历时确定。
时确定是否输入数据；白色格内要输入所需数据。
洪 演 算调 洪。
调洪演算
1.5 调洪演算调洪演算的基本原理是水量平衡，其方程为121221--22Q Q q qt t V V ++∆∆= 式中： Q 1、Q 2分别为计算时段Δt 始、末入库流量； q 1、q 2分别为计算时段Δt 始、末出库流量； V 1、V 2分别为计算时段Δt 始、末水库库容； Δt 为计算时段。
采用列表试算法，计算工作量较大，这里采用半图解法（单辅助线法）。
将水量平衡方程变形得：2212111222V q Q Q V q q t t +⎛⎫+=-++ ⎪∆∆⎝⎭式中右边为已知项，左边为未知项。
我们可以先确定q 与2V q t ⎛⎫+ ⎪∆⎝⎭之间的关系，绘制2q V q t ⎛⎫+ ⎪∆⎝⎭～的辅助曲线。
方法为由已知的q 查上游水位与泄流量关系曲线得上游水位H 上，在查水位库容关系曲线得相应的库容，Δt 为计算时段，在这里为24h ，进而求得对应的2V q t ⎛⎫+ ⎪∆⎝⎭。
从第一时段开始，由入库洪水过程和起始条件就可以知道Q 1、Q 2、q 1、V 1，由上式求得222V q t ⎛⎫+ ⎪∆⎝⎭，然后由2q V q t ⎛⎫+ ⎪∆⎝⎭～的辅助曲线查的对应的q 值即为q 2，然后按此方法依次计算q 。
计算过程如下，先确定q 与2V q t ⎛⎫+ ⎪∆⎝⎭之间的关系。
表1.19A q 与2V q t ⎛⎫+ ⎪∆⎝⎭关系表绘制2q V q t ⎛⎫+ ⎪∆⎝⎭～的辅助曲线:图1.7A 2q V q t ⎛⎫+ ⎪∆⎝⎭～的辅助曲线然后进行调洪演算，过程如下：表1.20A 调洪演算过程表图1.8A 调洪演算过程曲线调洪演算后的最大泄流量为两线的交点，表中计算的q max=4819m3/s，对应的Q=4800 m3/s，不相等，但很接近，则q max比4819m3/s稍微大些，参照图得q max=4825m3/s。
“水库调洪演算的原理和方法”教案讲义
2、蓄泄方程
qf(V) 或q ~ V
Z(m)
0
Z~q洞
Z~q溢
H H0
q溢M1BH3/2
q洞M2H01/2
溢洪 q(m3/s) 道
q~H Z ~V q~f(V)
泄洪 洞
方程不是显式，无法直接求解
三、水库调洪计算方法
1、列表试算法
（1）由Z~V曲线、泄流计算公式推求q~V曲线。
（2）推求水库下泄流量过程线q~t。
• 二、无闸溢洪道水库的防洪计算 • (一)下游无防洪要求时 • (1)拟定方案；(2)调洪计算；(3)计算坝高；
(4)计算各方案的投资(5)选定方案
• (二)下游有防洪要求
• 在防洪水利计算中需要考虑下游安全泄量的要求， 分别对枢纽标准与下游防洪标准的洪水进行调洪计 要。具体的步骤是：
• (1)假定不同的溢洪道宽度B方案。
• 水库防洪计算的内容，主要包括以下几点： • (1)根据库区地形、地质等条件，分析洪水特性及
灾害情况，考虑兴利库容与调洪(防洪)库容结合的 可能和程度，拟定苦干个泄流建筑物形式、位置、 尺寸，以及汛期运用方式的方案。 • (2)对各方案进行调洪计算，求得每个方案相应于 各种设计洪水的最大下泄流量、调洪库容和最高洪 水位。 • (3)计算各方案的大坝造价、淹没损失、泄流建筑 物投资、下游堤防造价、水库防洪效益等经济指标、 进行技术经济分析与比较，选择最佳方案。
（4）推求Vm、V调、Zm。
V调=Vm-V汛限
q~V
qmax
Vm
Z~V
Zm
列表试算法小结：
洪水过程线 起始条件
q~V曲线的绘制←Z~V q~t曲线的推求←试算
qm的推求←绘图 Zm、Vm推求←查图（表）
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 泄流孔口底槛高程。
水库兴利调节及调洪演算
⽔库兴利调节及调洪演算《⽔资源规划及利⽤》课程设计说明书在过去的⼀周⾥，我们进⾏了《⽔资源规划及利⽤》的课程设计，我们做的是梅⼭⽔库的三个典型年的兴利演算及其发电量的计算还有68年梅⼭⽔库的⼀次调洪演算。
现在，我们基本上已经设计完毕，通过课程设计，进⼀步加强了我们对所学内容的理解⽔平和应⽤能⼒，培养了我们分析问题与解决问题的能⼒。
下⾯，就是我们在课程设计的过程。
概述1．梅⼭⽔库概况：梅⼭⽔库位于淮河⽀流史河上游的安徽省⾦寨县境内，东与淠河西源为邻，西与灌河隔岭为界，南源于⼤别⼭北麓，北距史河⼊淮⼝130km。
⽔库流域南北长约70km，东西宽约40km，流域⾯积1970km2。
梅⼭⽔库按500年⼀遇洪⽔设计，5000年⼀遇洪⽔校核，设计洪⽔位137.66m，校核洪⽔位139.93m，正常蓄⽔位128.0m，汛限⽔位125.27m，死⽔位94.00m，总库容22.64亿m3，兴利库容9.57亿m3，死库容1.26亿m3，为年调节⽔库。
梅⼭⽔库现有⽔电站装机容量为4万kW，4台发电机组，单机最⼤过⽔流量29.8m3/s，电站主要结合灌溉供⽔或利⽤泄洪弃⽔发电，原则上不单独为发电⽬的⽽放⽔⼊横排头⽔库。
现状情况下多年平均发电量为9925万kW·h。
2．设计内容：①.熟悉资料，绘出相关曲线；②.根据梅⼭⽔库的1969-2008年的⼀年中各旬的⼊库流量资料，定出对应设计保证率为10%，50%，90%的设计典型年；③.分别对梅⼭⽔库的三个典型年的⼀、⼆种⽅案进⾏兴利调节和⽔能计算，求出各种⽅案的年平均发电量，并且⽐较各⽅案的利弊。
④.⽤第⼀种⽅案分别对梅⼭⽔库的1969年的⼀次洪⽔进⾏调洪演算。
3．设计⽅案：兴利⽅案：①⽅案⼀：正常蓄⽔位：128 m，汛限⽔位：125.27 m。
②⽅案⼆：正常蓄⽔位：128 m，汛限⽔位：124.57 m。
调洪⽅案：汛限⽔位125.27m，防洪⾼⽔位130.2，控泄流量75m3/s 4．梅⼭⽔库兴利调度原则：⾮汛期为5-10⽉，汛期为6-9⽉，其中9⽉为过渡期，即从9⽉份开始⽔库开始慢慢蓄⽔，并在9⽉末时保证⽔库⽔位在死⽔位和正常⾼⽔位之间，10-5⽉份⽔库处于⾮汛期，这段时间保证⽔库⽔位在死⽔位和正常⾼蓄⽔位之间，并保证⾮汛期末时⽔库⽔位在死⽔位与汛限⽔位之间。
借助Excel电子表格实现水库泄流全过程调洪演算
借助Excel电子表格实现水库泄流全过程调洪演算【内容提要】借助Excel电子表格实现水库泄流全过程调洪演算，以井研县毛坝水库枢纽病害整治工程初步设计中，其上游水库红岩水库p=0.1%校核洪水过程的调洪演算为算例，探索出两种不同的计算方法：一种是利用Excel电子表格自身的函数和公式运算，实现半自动化的水库调洪演算；二是以Microsoft Excel2010以上版本为依托，通过VBA语言编程来实现全自动化的水库调洪演算。
两种方法的计算结果，与本人当年使用日本夏普便携式PC-1500计算机通过Basic语言编程的水库多功能调洪演算程序的计算结果完全一致。
作者推荐第二种方法，它能在大约1至2秒之内实现多行多列上百数据数百次试算的入库和出库洪水流量演算全过程。
由于它的快速和智能化设计，给水库洪水预报盈得时间，是编制已成水库防洪预案和水利水电工程设计及病害整治工作的有力助手。
【正文】现以井研县毛坝水库1995年工程整治中红岩水库调洪演算为算例，展现更先进的计算方法。
本算例在计算入库洪水过程中，所涉及暴雨洪水参数来源有二：一是四川省水利电力局水文总站1979年10月编制的《四川省水文手册》；二是四川省水利电力厅1984年6月编制的《四川省中小流域暴雨洪水计算手册》。
虽然后者的暴雨等值线图已经被近20多年来新增实测暴雨资料的补充分析成果所代替。
四川省水文水资源勘测局已增补到2000年的暴雨资料，于2010年12月出版了新的《四川省暴雨统计参数等值线图》，按理应以新的暴雨资料为依据，但为了保持该工程整治设计的原貌，为了与当时本人使用PC-1500计算机和原来的调洪演算程序的运算结果进行比较，仍然使用《四川省中小流域暴雨洪水计算手册》所附暴雨等值线图计算该水库设计洪水过程线。
一、水库的基本情况毛坝水库位于岷江左支泥溪河上游，井研县天云乡境内。
水库大坝的地理坐标为E103°57′24″，N29°45′37″。
暴雨产流计算(推理公式-四川省)演示教学
1
0.938
0.835
0.77
1
0.928
0.832
0.771
1
0.972
0.909
0.861
1
0.815
0.539
0.423
1
0.962
0.875
0.81
根据实测资料的对应分析，μ值随集水面积的减小而增大。在四川其平均变化关系为：μ＝k·F－0.19，k值变化归
计算编号
地区
流域地形地貌
1
青衣江～鹿头山暴雨区 相对高差在200m以上，地势较陡，切割较深，植被较好，有部分荒山或坡地
9
#NAME?
10
#NAME?
11
#NAME?
12
#NAME?
13
#NAME?
分区名称
分区范围
盆地腹部丘陵区
岷江、沱江中 游及涪江、嘉
岷沱江下游平行岭谷区 岷江、沱江下游及平行岭谷地区
长江南岸区
高县至綦江长江南岸地区
赤水河古蔺区
赤水河古蔺地区
乌江下游及巫山区 乌江下游及巫山地区
沅江区
秀山地区
大巴山暴雨区
1 0.95 0.8 0.6 0.4 0.2 0.1 0.05
0
0.76 0.52
0.79 0.56
单峰
Ⅰ 0 0.2 0.26 0.34 0.4 0.5 0.62 0.75 0.9 0.98 1.09 1.24 1.44 1.79 2.1 2.38 3.1
四川省分区不同流域面积F（km2）综合暴雨24h面深折减系数αt
分区
上限面积km2
25
100
Ⅰ1
240
1
0.986
暴雨产流计算(推理公式_四川省)
Ⅰ
地区
Ⅰ1
洪水过程线
雨型
3
计算参数
编号
5
东部地区 双峰 Ⅱ 金沙江
时段
10′ 1h 3h 6h 24h
点雨量均值 Ht(mm) 18 46 60 86 110
变差系数 Cv
0.35 0.42 0.45 0.53 0.55
偏差系数 Cs/Cv 3.5 3.5 3.5 3.5 3.5
模比系数 KP
1.67 1.82 1.88 2.05 2.09
Ⅰ
0.94
1
Ⅱ
1
1
Ⅲ
1
1
Ⅳ
1
1
Ⅴ
1
1
雨型
历时T
单峰洪水
44.03
双峰洪水
61.92
洪水过程线计算参数
东部地区
单峰
川西南地区
双峰
设计洪水流量 WP(万m3) #NAME?
Q0(m3/s)
确定以日还是以时计算暴雨
1.83
日
2.58
日
计算参数
是
否
值m=0.2J－1/3，m范围为m=0.3～2.5。
影响因素
25
0.811
110
1
0.63
70
1
0.734
70
1
0.473
40
1
6h面深折减系数αt
100
300
0.966
0.901
0.924
0.802
0.932
0.822
0.888
0.7
0.851 0.887 0.795 0.811 0.868 0.763 0.738 0.84 0.4 0.78
水库调洪演算的原理和方法
V t
q 2
f2 (Z )
V t
q 2
f3(Z) q
f3 (Z )
f1(Z )
f2 (Z )
Q (m3 / s), V q (m3 / s), V q (m3 / s)
t 2
t 2
调洪计算半图解法的双辅助线
水利水能规划
水库调洪计算的半图解法
V调=Vm-V汛限
Vm
Z~V
Zm
【例 题】
水利水能规划
【补偿调节】
水库
QB=q+Q区
Q
A
QB
qB
坝址
Q区 6h
防
洪
保
护 区
河流
水利水能规划
水利水能规划
水利水能规划
水利水能规划
水利水能规划
水利水能规划
水利水能规划
水利水能规划
水利水能规划
【补偿调节】
水库
QB=q+Q区
Q
A
QB
qB
坝址
水利水能规划
水库调洪计算的半图解法
由上节知道列表试算法麻烦工作量大，故人们比较喜欢 用半图解法。
Q1
Q2 2
Δt
q1
q2 2
Δt
V2
V1
Q Q1 Q2 2
等式两边同时除以△t,并移项
Q V1 q1 V2 q2 t 2 t 2
第十四章 水库防洪计算
水利水能规划
• 三、有闸溢洪道水库的防洪计算
水利水能规划
水利水能规划
• 四、具有非常泄洪设施水库的防洪计算
调洪计算方法
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。
第十三章 水库洪水调节及计算
简化三角形解析法、简化三角形图解法
30
第十一章 水库洪水调节及计算
三、水库调洪计算的基本方法
原理
流量
入库洪水总量W 1 2 QmT
调洪库容：Vm W 1 2 qmT qm QmT (1 ) Qm
1 2
Qm
A
Vm
O
qm
C
B
qm Vm W (1 ) Qm
时间
T
Vm 或 qm Qm (1 ) W
推求水库下泄流量过程线 q ~ t ，最大下泄量 qm ，调 洪库容 V洪 和水库最高洪水位 Z洪 。
常用的方法有列表试算法、半图解法和简化三角形法
等。
18
第十一章 水库洪水调节及计算
三、水库调洪计算的基本方法
1.列表试算法
步骤： (1)计算并绘制水库的蓄泄曲线q ~ V；
(2) 分析确定调洪的起始条件：起调水位和起调库
下游有防洪任务时，水库的防洪库容和泄洪建筑 物尺寸应在确保水库主要建筑物(大坝等)防洪安全的 前提下，同时满足下游防洪要求。( 有闸门控制) 方案选择 根据计算成果，通过技术和经济比较，选出最优 方案。
15
第十一章 水库洪水调节及计算
二、水库调洪计算的原理
水库调洪计算的基本原理是逐时段地联立求解水量 平衡方程和水库的蓄泄方程。
一、水库的调洪作用
有闸门控制时， 水库的调洪作用
q兴
胸墙 闸门
q允
Z堰
溢洪道
12
第十一章 水库洪水调节及计算
一、概述
2.水库调洪计算的任务与过程
规划设计阶段
管理运行阶段
由
13
第十一章 水库洪水调节及计算