水文预报课程设计河海.docx
《水文预报课程设计河海.docx》由会员分享，可在线阅读，更多相关《水文预报课程设计河海.docx（25页珍藏版）》请在冰点文库上搜索。
水文预报课程设计河海
第一章基本任务
本次课程设计从2013年7月8日至2012年7月12日，主要任务是对广东省东江一级支流西枝江白盆珠水库的上游宝口流域编制预报方案与产汇流计算。
其基本任务为：
任务一：
根据已给的资料、参数及做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较；每人计算两年。
任务二：
根据已给设计暴雨资料、参数及做过的习题,自己编写程序，将流域作为整体进行次洪产流量、划分水源、直接径流汇流、地下径流汇流计算；绘出直接径流过程、地下径流过程、总的流量过程。
第二章基本资料
2.1流域概况
白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积8562。
流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。
暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。
降雨年际间变化大，年内分配不均，多年平均降雨量为1800，实测年最大降雨量为3417，汛期4—9月降雨量占年降雨量的81％左右；径流系数0.5-0.7。
流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。
台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。
流域上游有宝口水文站，流域面积为5532，占白盆珠水库坝址以上集雨面积的64.6％。
白盆珠水库有6年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：
流域内有7个雨量站，其中宝口以上有4个。
雨量站分布较均匀，有6年逐日降水资料和时段降水资料；宝口水文站具有6年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般2—3小时（h），有时更短；一次洪水总历时2~5天（d）。
2.2基本资料
2.2.1任务一参数及相关实测资料
计算流域为白盆珠水库上游的宝口流域，其流域面积为5532。
该流域内有四个雨量站：
禾多布、马山、高潭、宝口，其权重系数分别为：
0.33、0.14、0.33、0.20。
本次课设提供的资料为1987年至1992年该流域径流量、蒸发量及各雨量站的日实测资料，其中我计算年份为1987年与1988年，所使用的相关参数值如下表所示：
（）
（）
（）
B
C
W（）
（）
（）
（）
140
20
60
0.2
0.16
0.001
110
10
40
60
其中蒸散发折算系数需编程优选，优选范围为0.90-1.30。
优选原则为：
计算的2年内每年的年径流相对误差尽可能不超过5%。
2.2.2任务二参数及相关实测资料
任务二产流参数为任务一中日模优选的值，为10.03h，产流参数用日模优选的值，其他参数取任务1中各组的参数，初始张力水蓄量取各组的容量值。
任务二汇流参数为0.968，0=55.3m3，单位线序号1至11对应的Q（m3）分别为：
0、40、80、130、100、80、48、20、10、5、0
2.3产流方式论证
宝口流域地处我国南方湿润地区，气候暖热，雨量充沛，多年平均降雨量为1800>1000，年径流系数在0.5—0.7之间，大于0.4；流域内土质疏松，植被良好，不易超渗；一次洪水的流量过程陡涨缓落，持续时间2~5d左右。
从流域的气象条件，下垫面条件和流量过程的分析知，该流域降雨径流关系具有蓄满产流的特点，可以按蓄满产流建立产流量预报方案。
第三章计算公式
3.1降雨量计算
由已知资料知，该地区雨量站分布均匀，且宝口以上有四个雨量站：
禾多布、马山、高潭、宝口，其权重系数分别为：
0.33、0.14、0.33、0.20，所以宝口流域平均降雨量0.33*P1+0.14*P2+0.33*P3+0.20*P4。
3.2蒸散发量计算
蒸散发计算采用三层蒸发计算模式，即：
上层蒸发量：
下层蒸发量：
*
深层蒸发量：
*
总蒸发量：
式中：
为流域蒸发能力（）；为下层土壤含水量（）；为下层土壤含水量（）；C为蒸发扩散系数。
三层蒸发模式按照先上层后下层的次序，具体计算为：
1）当>时，00
2）当<，>*时，（）*0
3）当<,C*（）<4）当<,3.3降雨产流量计算
1）a＋≤时，＋W－＋（1－
）1
2）a＋＞时，－（－W）
3.4水源划分计算
通过稳渗率可划分产流中的直接径流和地下径流。
次洪的各水源分量为：
3.5汇流计算
根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。
其计算公式：
直接径流过程：
（i）（i）*U
地下径流过程：
（i）*
（1）+
（1）*（i）*U
总的流量过程：
Q（i）（i）（i）
第四章计算结果
4.1任务一结果
经过程序优选，0.98
年份
实测R（）
计算R（）
绝对误差（）
相对误差（%）
1987
1080
1100
19.3
1.79
1988
824
813
-10.2
-1.23
4.2任务二结果
t（i）
P（i）
R（i）
（i）
（i）
（i）
（i）
Q（i）
2004/9/2312:
00
14
12.7
2.7
10
0
69.9
69.9
2004/9/2315:
00
14.1
12.9
2.9
10
10.8
84.1
94.9
2004/9/2318:
00
11
9.7
0
9.7
33.2
97.3
130
2004/9/2321:
00
18.7
17.4
7.4
10
58.3
111
169
2004/9/240:
00
29.7
28.5
18.5
10
94.3
123
218
2004/9/243:
00
36
35.1
25.1
10
184
136
320
2004/9/246:
00
38.3
37.5
27.5
10
381
148
529
2004/9/249:
00
7.8
6.9
0
6.9
645
154
799
2004/9/2412:
00
30.5
29.6
19.6
10
799
166
965
2004/9/2415:
00
42.6
41.7
31.7
10
874
177
1050
2004/9/2418:
00
93.8
92.9
82.9
10
864
188
1050
2004/9/2421:
00
84.1
83.3
73.3
10
1220
198
1420
2004/9/250:
00
47.6
46.8
36.8
10
1770
208
1980
2004/9/253:
00
56.4
55.3
45.3
10
2370
218
2590
2004/9/256:
00
50.4
49.3
39.3
10
2650
227
2870
2004/9/259:
00
16.5
15.4
5.4
10
2600
236
2840
2004/9/2512:
00
5.8
4.8
0
4.8
2360
237
2600
2004/9/2515:
00
8.3
7.2
0
7.2
1860
241
2100
2004/9/2518:
00
2.6
1.5
0
1.5
1250
236
1480
2004/9/2521:
00
0.2
0
0
0
774
228
1000
2004/9/260:
00
0
0
0
0
396
221
617
2004/9/263:
00
0
0
0
0
168
214
382
2004/9/266:
00
0
0
0
0
72.8
207
280
2004/9/269:
00
0
0
0
0
25.0
200
225
2004/9/2612:
00
0
0
0
0
2.70
194
197
2004/9/2615:
00
0
0
0
0
0.00
188
188
2004/9/2618:
00
0
0
0
0
0.00
182
182
2004/9/2621:
00
0
0
0
0
0.00
176
176
第五章误差统计与分析
5.1精度评定
从计算结果可见，年产流量绝对误差均小于100，所产流量的相对误差均小于5%。
精度统计表明，率定的模型参数是基本合理的。
但由于课设时间限制与任务要求，每位同学只用两年的资料来率定，不满足《水文情报预报规范》中规定：
洪水预报方案要求使用样本数量不少于10年的水文气象资料，其中应包括大、中、小水各种代表性年份，并保证有足够代表性的场次洪水资料。
显然0.98的结果还是存在一定问题的。
5.2误差分析
影响流域降雨径流过程的因素很多，利用蓄满产流新安江模型的结构与参数能够很好反映湿润地区降雨径流过程的主要规律和特点，因而能获得较好的精度。
但是模型本身以及模型计算中有很多概化，会造成误差。
造成宝口流域产汇流计算方案误差来源主要有以下几个方面：
1.资料代表性的影响。
我们每位同学只用两年实测资料来率定，资料必不能满足时期要求及代表性要求。
实际操作中，应要求有12年以上的连续的对未来有代表性的实测资料，其中10年为率定期，2年为检验期，且这种资料应具备丰、平、枯水年代表性，资料系列前后一致，受人类影响较小。
所以仅用两年的资料率定势必造成一定的误差，而且误差可能还很大。
2.量测误差。
实测的降雨、蒸发、径流量等水文气象信息及河流、湖泊、地形等下垫面信息是研制预报模型或编制洪水预报方案或进行作业预报的主要依据，在现有站网、仪器设备、观测技术条件下，各种信息的时空变化是难以准确反映的，加上受自然因素等客观条件影响，势必造成各种信息的量测误差。
3.模型结构误差。
在该蓄满产流新安江模型中，有很多将非线性现象概化为线性现象或者将某些随机因子近似作为确定因子描述等都会带来误差。
4.模型参数误差。
模型参数是根据输入，通过模型计算输出，再将输出过程与实测过程进行比较，用系统识别的方法作优化调试的，上述所率定出的模型参数可能不是最优。
5.人类活动的影响。
随着社会经济的快速发展，人类活动的影响加剧，流域内可能新建了一些大中型水库或其他工程措施。
影响了该地区原有的产流特点，如汛期提前。
5.3实时校正
模型计算值与实测值直接总是存在一定的误差。
造成两者间误差的因素很多，若针对每一个单一因素是难于描述或预见的，一般采用实时校正模型来解决。
实时校正模型常用的有卡尔曼滤波、自回归模型等。
第六章编写程序
1
1（）
2
2（）
3
2
1（）
1
2
2（）
m=040
=0.9+0.01*m
（2）,
（2）,
（2）,（3）
W0,0,0,0
0=10:
0=11:
0=12
W0=0+0+0
i=12
=1986+i'用于定位
4<>0100<>0
=365
=366
'读入数据
&"\..\\"&&""#1
t（）
（）
（）
P1（）
P2（）
P3（）
P4（）
j=1365
#1,t（j）,（j）,（j）,P1（j）,P2（j）,P3（j）,P4（j）
j
#1
'参数赋值
=2:
=3:
C=4
=5:
B=6:
=7
P（）,（）,（）
E（）,（）,（）,（）
W（）,（）,（）,（）
R（）
j=1
W（0）=W0:
（0）=0:
（0）=0:
（0）=0
P（j）=0.33*P1（j）+0.14*P2（j）+0.33*P3（j）+0.2*P4（j）
（j）=*（j）
（j）=P（j）-（j）
（j）>0
（j）=（j）'蒸发
（j）=0
（j）=0
A,
=*（1+B）/（1-）
（W（j-1）-）<0.01
A=
A=*（1-（1-W（j-1）/）^（1/（1+B）））
A=A+（j）
A<=
R（j）=（j）+W（j-1）-+*（1-A/）^（1+B）
R（j）=（j）+W（j-1）-
（（j-1）+（j）-R（j））>'第一层
（j）=
（（j-1）+（j）-R（j）-）+（j-1）>'第二层
（j）=
（j）=W（j-1）+（j）-R（j）-（j）-（j）
（j）=（（j-1）+（j）-R（j）-）+（j-1）
（j）=（j-1）
（j）=（j-1）+（j）-R（j）
（j）=（j-1）
（j）=（j-1）
'不产流，按三层蒸发算
（j-1）+P（j）>=（j）
（j）=（j）
（j）=0
（j）=0
（j）=（j-1）+P（j）
（j-1）>*C
（j）=（（j）-（j））*（（j-1）/）
（j）=0
（j-1）>=C*（（j）-（j））
（j）=C*（（j）-（j））
（j）=0
（j）=（j-1）
（j）=C*（（j）-（j））-（j）
（j）=（j-1）+P（j）-（j）
（j）=（j-1）-（j）
（j）=（j-1）-（j）
E（j）=（j）+（j）+（j）
W（j）=（j）+（j）+（j）
j
W0=W（）:
0=（）:
0=（）:
0=（）
j=1-1
（i）=（i）+（j）*3.6*24/553'年实测径流深
j
j=1
（i）=（i）+R（j）'年计算径流深
j
（i）=（（i）*10+0.5）/10
（i）=（（i）*10+0.5）/10
（i）=（（（i）-（i））*10+0.5）/10
（i）=（（（i）-（i））/（i）*100+0.5）/100
（i）<=0.05
=1987
1;（i）
2;（i）
i
（i-2）<=0.05（i-1）<=0.05
（i）=（（i-2）+（i-1））/2'（i）为两年的平均值3;（i）
m
1=（"请输入您选择的"）'使用者选用自己认为合适的
=1
（2）,
（2）,
（2）,（3）
0=10:
0=11:
0=12
W0=0+0+0
i=12
=1986+i'用于定位
4<>0100<>0
=365
=366
'读入数据
&"\..\\"&&""#1
t（）
（）
（）
P1（）
P2（）
P3（）
P4（）
j=1365
#1,t（j）,（j）,（j）,P1（j）,P2（j）,P3（j）,P4（j）
j
#1
=2:
=3:
C=4
=5:
B=6:
=7
P（）,（）,（）
E（）,（）,（）,（）
W（）,（）,（）,（）
R（）
j=1
W（0）=W0:
（0）=0:
（0）=0:
（0）=0
P（j）=0.33*P1（j）+0.14*P2（j）+0.33*P3（j）+0.2*P4（j）
（j）=*（j）
（j）=P（j）-（j）
（j）>0
（j）=（j）'蒸发
（j）=0
（j）=0
=*（1+B）/（1-）
（W（j-1）-）<0.01
A=
A=*（1-（1-W（j-1）/）^（1/（1+B）））
A=A+（j）
A<=
R（j）=（j）+W（j-1）-+*（1-A/）^（1+B）
R（j）=（j）+W（j-1）-
（（j-1）+（j）-R（j））>'第一层
（j）=
（（j-1）+（j）-R（j）-）+（j-1）>'第二层
（j）=
（j）=W（j-1）+（j）-R（j）-（j）-（j）
（j）=（（j-1）+（j）-R（j）-）+（j-1）
（j）=（j-1）
（j）=（j-1）+（j）-R（j）
（j）=（j-1）
（j）=（j-1）
'不产流，按三层蒸发算
（j-1）+P（j）>=（j）
（j）=（j）
（j）=0
（j）=0
（j）=（j-1）+P（j）
（j-1）>*C
（j）=（（j）-（j））*（（j-1）/）
（j）=0
（j-1）>=C*（（j）-（j））
（j）=C*（（j）-（j））
（j）=0
（j）=（j-1）
（j）=C*（（j）-（j））-（j）
（j）=（j-1）+P（j）-（j）
（j）=（j-1）-（j）
（j）=（j-1）-（j）
E（j）=（j）+（j）+（j）
W（j）=（j）+（j）+（j）
j
W0=W（）:
0=（）:
0=（）:
0=（）
'输出文件
&"\"&&""#2
j=1
#2,t（j）,P（j）,（j）,（j）,（j）,（j）,（j）,E（j）,（j）,（j）,（j）,W（j）,R（j）
j
#2
'年统计
j=1-1
（i）=（i）+（j）*3.6*24/553'年实测径流深
j
j=1
（i）=（i）+R（j）'年计算径流深
j
（i）=（（i）,"0.00"）
（i）=（（i）,"0.00"）
（i）=（i）-（i）
（i）=（（i）-（i））/（i）*100
（i）=（（i）,"0.00"）
（i）=（（i）,"0.00"）
"实测径流（"&i&"）="&（i）
"计算径流（"&i&"）="&（i）
"绝对误差（"&i&"）="&（i）
"相对误差（"&i&"）="&（i）
i
3
1（）
'任务二
'读入数据
&"\..\\"#1
（1）
=+1
t（）
（）
P1（）
P2（）
P3（）
P4（）
#1,t（）,（）,P1（）,P2（）,P3（）,P4（）
#1
'参数赋值
=1:
=20:
=60:
C=0.16
=140:
B=0.2:
=0.001
P（）,（）,（）
E（）,（）,（）,（）
W（）,（）,（）,（）
R（）
'初始状态
W（0）=
（2）+（3）+（4）:
（0）=
（2）:
（0）=（3）:
（0）=（4）
j=1
P（j）=（（0.33*P1（j）+0.14*P2（j）+0.33*P3（j）+0.2*P4（j））*10+0.5）/10
（j）=*（j）
（j）=P（j）-（j）
（j）>0'产流计算
（j）=（j）
（j）=0
（j）=0
A,
=（1+B）*/（1-）
（W（j-1）-）<0.1
A=
A=*（1-（1-W（j-1）/）^（1/（1+B）））
A=A+（j）
A<=
R（j）=（j）+W（j-1）-+*（1-A/）^（1+B）
R（j）=（j）+W（j-1）-
R（j）=（R（j）*10+0.5）/10
（（j-1）+（j）-R（j））>'由上至下，逐层填蓄
（j）=
（（（j-1）+（j）-R（j）-）+（j-1））>
（j）=
（j）=W（j-1）+（j）-R（j）-（j）-（j）
（j）=（（j-1）+（j）-R（j）-）+（j-1）
（j）=（j-1）
（j）=（j-1）+（j）-R（j）
（j）=（j-1）
（j）=（j-1）
（j-1）+P（j）>=（j）
（j）=（j）
（j）=0
（j）=0
（j）=（j-1）+P（j）
（j-1）/>C
（j）=（（j）-（j））*（j-1）/
（j）=0
（j-1）>=C*（（j）-（j））
（j）=C*（（j）-（j））
（j）=0
（j）=（j-1）
（j）=C*（（j）-（j））-（j）
（j）=（j-
升级会员 