水文预报课程设计河海.docx
- 文档编号:11791757
- 上传时间:2023-04-01
- 格式:DOCX
- 页数:25
- 大小:152.63KB
水文预报课程设计河海.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-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 水文 预报 课程设计
![提示](https://static.bdocx.com/images/bang_tan.gif)