水文预报课程设计Word下载.docx
《水文预报课程设计Word下载.docx》由会员分享，可在线阅读，更多相关《水文预报课程设计Word下载.docx（27页珍藏版）》请在冰豆网上搜索。
陶湾。
伊河洪水为暴雨所形成，受降雨特性影响，洪水过程陡涨陡落，洪水多发生在7～10月份，特大洪水发生在7～8月份。
流域内水系和站网分布见下图。
流域内水系和站网分布图
四、设计步骤及技术要求
（一）编制产流量预报方案
1．绘制汛期栾川站流量过程线和相应的降雨量过程线
根据1964年和1978年栾川站洪水水文要素摘录表可以绘制出汛期栾川站流量过程线，见图如下。
图1流量过程线
2．计算流域平均次降雨量P
由各站累积雨量摘录表（表1各站累积雨量摘录表如下）能够算出1964年和1978年流域平均次降雨量P1964年=，P1978年=。
由各站累积雨量摘录表知道各个时刻的降雨量，然后把各个时刻的雨量相加的该时间段的降雨总量，用降雨总量除以该段时间的天数得到该流域的平均次降雨量。
3．分析栾川站流量过程的退水规律，制作退水方案
根据已知资料，算出栾川站流量的退水流量过程，其详细计算成果表见表2所示。
4．划分洪水，计算各次洪水的实测径流深R
水源的划分：
蓄满产流条件下，产流量是产流面积上的降雨量，能形成地下径流的也只是产流面积上由稳定下渗率fc形成的水量。
因此水源公式为
当h＞fc时，rg=
①
当h≤fc时,
②
一次洪水地下水和直接径流总产量分别为
③
④
式中①、②就是划分时段和次洪水源的公式。
由以上公式可以计算出表3和表4的结果，其成果见附表中的表3及表4所示。
5.初定蓄满产流模型参数
蓄满产流模型参数主要有蓄水容量曲线参数b，Wm，蒸发模型参数Kc，C分水原参数fc共五个。
其中对产流量预报影响最大的是Kc。
模型参数在预报中起重要作用，故在建立产流模型前要确定其初值，供模型的调试分析作用。
（1）流域稳定下渗率fc的推求
Fc值可用实测降雨径流资料反推而得，方法较多，常用试算法。
对一次洪水，首先用上一步所述方法求出实测径流深
，然后根据①、②，假设不同的fc，算出
，直至由降雨资料计算得
与实测值相等，则fc值求定。
为了使fc值有一定的精度，选择洪水时应选用次洪产流量计算误差较小和能达全流域产流的洪水，并分析多次，如差别不大，可取平均值。
有时，需对fc值作古算，精度要求不高时，可选择雨初基本能达全流域产流而雨强又较大的洪水，算出
，并统计出降雨历时T，则
由多次洪水分析的fc值可能会差别较大，应该具体分析，不能一概简单地取平均值。
如因未计时段内下渗历时，会使fc值受时段均化而变小，时段越长，这种现象越明显，但fc有自身规律，尚需研究。
（2）蒸发模型参数Kc的推求
通过分析流域的降雨径流资料，可得到较为准确的初值，应用水量平衡原理，进行流域蒸发接近于Ep时期的水平衡计算，令湿润期为T，在T时段内水量平衡方程为
式中
为T时段内流域蒸发量，
R为总降雨量和由其形成的产流量，均为实测，
为时段始末土壤蓄水量之差。
如在汛期中某湿润期时段始末都能蓄满，则
=0，湿润期的蒸发也接近于Ep，并假定其相等。
另外，统计出T时期蒸发器观测的水面蒸发量
按Kc的意义，可用下式估算出
⑤
估算的Kc值中综合了多种误差，故必须多分析几次洪水，求其平均值才比较可靠。
另外，也可用多年平均值，据式⑤算出Kc值，但只能作为Kc的下限值，因为实际蒸发量E总是小于Ep的。
（3）
的推求
一般天然流域，
值的确定，通常由实测降雨径流资料分析而得。
是流域缺水容量，故应尽量选取前期很干旱（
≈0，缺水量最大），本次降雨又较大，能达到全流域产流（缺水量全部能不足）的洪水，其降雨损失量即
值，可多分析几次洪水取其较大者作为
初值。
有时不易找到满足上述条件的洪水资料，则可找出长期干旱以后的多次洪水，取至达全流域蓄满时作为计算时段，进行水量平衡计算，用下式定出
。
但此时由于时间较长，故需估算出蒸发期
值，可参考同期
确定。
（4）b的推求
B值在部分流域面积产流时起作用，全流域产流时，由于
为1.故b值大小与产流量无关。
在湿润地区，b只影响小洪水或大洪水的前期产流量（不影响总产流量），故对R影响不大；
但干旱地区则局部产流较多而全流域蓄满情况较少，故b值影响较大，根据经验，一般取值为~，山区大于平坦地区，b值一般随流域面积增大而增大。
（5）C的推求
C值的确定，由于研究不够，常取经验数据，约~.一般作物多的地区旱季深根植物散发量较大，C值也较大。
6．应用蓄满产流模型计算各次洪水的径流深R
（1）确定雨量站，计算流域平均雨量P；
（2）分析流域退水曲线，划分Q~t，计算次洪实测产流量R0；
（3）初定参数b，
，Kc，C,fc；
（4）计算次洪Rc，根据蒸散发计算模型，计算逐日E，并计算出
，由
结合P~R关系（或查图）算出逐日Rc，据此统计出次洪水径流深Rc值，即R值。
由以上四个步骤即可计算出各次洪水的径流深R，其计算成果见下表4所示。
7.对方案进行精度评定
模型的误差有三个方向。
一是原始资料的观测误差，二是模型计算方法和参数确定方面的误差，三是模型的理论误差。
原始资料的误差，除了较明显的外，一般难以理解，但是具有这样的概念，有助于分析判断。
计算方法的误差涉及面比较广，主要是降雨径流关系各要素P，R0，W0的计算及其中的取定，常见一下几个方面。
（1）流域平均降雨量的计算
很多流域表明，面雨量的计算误差常成为点据误差的主要因素。
首先是测站的控制性问题，在雨量分布不均匀的地区及季节，现有站网往往难以控制面雨量，这个问题实际上很难完全解决。
其次，面雨量的各种计算方法基本依据是点雨量对其附近的面积往往有代表性，但是，由于雨型空间多变，加上测站控制性问题，采用何种方法最为合适或如何结合雨型采用不同的方法，须对多年计算结果进行对比分析，了解计算误差后，才能做出判断和改进。
一般流域上游降雨较大，而站网稀疏，易漏掉暴雨而使算得的流域平均雨量偏小，这种情况可通过洪水资料的分析比较来确定。
分块计算的流域，其子流域的平均雨量也存在上述类似的情况。
（2）次洪R0的计算
误差来自采用统一的地下水退水曲线，没有考虑不同退水规律的差异，地下水退水本身也有季节性变化。
此外，有些小洪水，虽流量较小，却并非全是地下水，还包括直接径流，往往出现小洪水比大洪水退水快，孤立洪水比连续洪水退水快的现象。
因此，在退水曲线变化较大或退水分割误差较大时，要考虑实际情况，取不同的退水曲线或改进计算方法。
（3）蒸散发量计算
影响蒸散发计算的参数主要有两个即Kc和Wm。
①Kc
Kc决定了Ep，常对蒸散发能力计算起主要作用。
如取得不准，将使点据出现系统偏差。
如Kc取得过小，Ep也小，则W0偏大，算的R0偏大，点据将偏离理论线上方，反之则偏于下方。
若取得合适，则无系统偏离。
有些地区Kc还成季节性变化，与植物散发有关，一般是夏季大于冬季。
②Wm的影响不如Wm明显，它通过蒸发计算来影响蒸发量。
若Wm取得过小，又采用一层蒸发模型计算时，会使损失量（Wm-W0）偏小，而使计算的产流量偏大，主要影响久旱后的洪水。
一般采用分层模型时，Wm较大，上述情况不易出现。
③b
B值只影响局部流域产流时的产流量，b值越大，R越大，对于湿润地区，大部分洪水为全流域蓄满，但对半湿润地区，则常有局部产流情况，需要注意的是b值对R的影响。
（4）模型简化生产的误差
①将流域作为整体计算
如将降雨径流各项要素都取流域平均值，不考虑流域上的不均性，则流域越大，不均匀性程度越剧烈，将影响计算成果。
首先是p在面上分布的不均匀性。
虽不均匀，但在最小面积上雨末也能蓄满，因全流域都蓄满，损失达最大值，为一常数，则对R计算无影响。
因此，只有在p不均匀而形成部分流域产流时，用平均p计算R与实际才会不同，其点据将偏于理论线下方，即计算径流深偏小。
当p在面上分布不均但最终达到全流域产流时，其产流量与p分布不均匀时相同，但产流量时程分配是不同的。
其次，在降雨分布明显不均匀的流域，往往引起流域缺水量也不均匀，其影响与p分布不均是相似的。
为消除误差，宜于各雨量站分别计算R，再求流域平均R值，这种分站计算法可以大致考虑到p分布不均和下垫面W不均的影响，提高精度。
②模型主要考虑下渗损失量
模型着重考虑下渗损失量，当其他损失项目突出时，会产生误差，常见的是中小水库蓄水，稻田蓄水等人类活动。
这些因素往往在短暂时间内较大地影响流域蓄水量。
如久旱之后，由于上列原因产生径流量比正常情况下小。
如衢县流域，根据调查，其影响量可达20~30mm。
显然，也会因后期水库泄水、农田放水而造成相反的误差。
这在半湿润地区影响尤为明显。
模型的理论误差，即发生超渗产流情况。
湿润地区在久旱后如遇短历时暴雨，可能发生超渗产流现象，其中R0大于用蓄满产流计算出的Rc值，点据偏于理论线下方。
但根据经验，这种情况并不多见，南方地区并非每年都能发生这种情况。
以上分析了蓄满产流模型误差来源及主要影响因素，但对于有误差点据，首先应找出可能产生这种误差的原因，然后逐个分析，排除无关因素，找出主因。
如Rc偏大，则可能是R0偏小，p偏小，Kc偏小，b偏大等。
然后逐个分析其可能性和论据。
其次，分析时要注意是否有系统性，要对所有资料作分析，因为有些因素是有系统影响的，如Kc值，对所有蒸发计算都有影响，b值则对局部产流的洪水有影响。
明确了误差的性质和找到原因后，需修改参数或计算方法，对全部资料重新计算，一般需反复多次才能使误差最小，最后确定方案和模型参数。
综合以上分析，根据所得的各次洪水径流深对方案进行精度评定，见下表6。
8．确定蓄满产流模型参数
（1）确定计算时段ΔT
由已知图表知：
△T=2h
（2）按ΔT摘录流量过程
经上一步的计算结果知△T=2h，所以按△T=2h时段摘录流量过程见下表7。
（3）由蓄满产流模型计算时段径流量
由蓄满产流模型计算时段径流量表见表8。
（4）对时段径流量进行修正
对时段径流量修正表见表9。
】
附表
表1各站累积雨量摘录表
日期
时间
站名
年
月
日
起
迄
P
∑P
1964
7
26
9:
00
16:
11:
40
18:
20:
4:
22:
6:
8:
27
10:
12:
14:
平均次降雨量
1978
4
23:
0:
2:
1:
55
5:
1
56
7:
表2退水资料摘录表
时段Δt=6小时
640727
750808
780805
790712
810715
830908
80
60
43
48
2
36
3
38
18
30
24
13
5
21
6
17
8
20
9
10
19
11
16
12
15
14
22
23
25
28
29
表3退水流量相应径流深计算表
Qt（m3/s）
ΣQt（m3/s）
Re（mm）
表4实测次洪径流深计算表
1964年
流量
时段长
平均流量
Δw
ΣΔw
126
172
322
320
314
305
13:
15:
221
17:
182
99
1978年
93
127
196
243
3:
235
178
123
107
39
表5产流计算表
p
E水
Em
Eu
El
E
P-E
Wu
WL
WO
R
31
升级会员 