水文预报课程设计(河海)

1、水文预报课程设计（河海）第一章基本任务本次课程设计从2013年7月8日至2012年7月12日，主要任务是对广东省东江一级支流西枝江白盆珠水库的上游宝口流域编制预报方案与产汇流计算。其基本任务为：任务一：根据已给的资料、参数与做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较；每人 计算两年。任务二：根据已给设计暴雨资料、参数与做过的习题，自己编写程序，将流域作为整体进行次洪产流量、划分水源、直接径流汇流、地下径流汇流计算；绘出直接径流过程、地下径流过程、总的流量过程。第二章基本资料2.1流域概况白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积

2、856 2。流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。降雨年际间变化大，年内分配不均，多年平均降雨量为1800，实测年最大降雨量为 3417，汛期49月降雨量占年降雨量的 81 %左右；径流系数0.5-0.7 o流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积水文预报课程设计（河海）层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大 乔木；平原则以种植农作物和经济作物为主，植被良好。流域上游有宝口水文站，流域面积为5532，占白盆珠水库坝址以上集雨面积的64.6 %。白盆珠水库有6年逐日入库流量资料、逐

3、日蒸发资料 和时段入库流量资料：流域内有 7个雨量站，其中宝口以上有 4个。雨量 站分布较均匀，有6年逐日降水资料和时段降水资料；宝口水文站具有6年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面 条件影响，洪水陡涨缓落，汇流时间一般2 3小时（h），有时更短；一次洪水总历时25天（d ）。2.2基本资料任务一参数与相关实测资料计算流域为白盆珠水库上游的宝口流域，其流域面积为553 2。该流域内有四个雨量站：禾多布、马山、高潭、宝口，其权重系数分别为：0.33、0.14、0.33、0.20。本次课设提供的资料为1987年至1992年该流域径流量、蒸发量与各雨量站的日实测资料，其中

4、我计算年份为1987年与1988年，所使用的相关参数值如下表所示：()()()BCW()()()()14020600.20.160.001110104060其中蒸散发折算系数需编程优选，优选范围为0.90-1.30。优选原则为：计算的2年内每年的年径流相对误差尽可能不超过5%。水文预报课程设计（河海）222任务二参数与相关实测资料日蒸散发降雨量（mm）月日时(mm)木多布马山咼潭宝口923121.36.29.921.617.3151.37.61620.612.6181.36.26.414.915.9211.38.817.229.418.524241.22534.835.324.630.929.

5、929.243.937.860.938.624.846.93390.96.97.56.112.3120.928.329.934.228.5150.925.642.739.875.4180.993.9137.612413.2210.985.390.88575.925240.851.547.449.238.531.139.870.342.197.761.143.247.361.545.991.120.513.315.813.1121.110.581.83.3151.17.48.47.610.9181.11.82.82.14.6211.10.200.3026241.2000032.1000062.1

6、000092.10000122.10000152.10000182.100002120000任务二产流参数为任务一中日模优选的值，为10.03h，产流参数用日模优选的值，其他参数取任务 1中各组的参数，初始张力水蓄量取各组的 容量值。任务二汇流参数为0.968，o=55.3m 3，单位线序号1至11对应的Q（m3）分别为：0、40、80、130、100、80、48、20、10、5、02.3产流方式论证宝口流域地处我国南方湿润地区，气候暖热，雨量充沛，多年平均降 雨量为1800>1000 ，年径流系数在0.5 0.7之间，大于0.4 ;流域内土 质疏松，植被良好，不易超渗；一次洪水的流量过

7、程陡涨缓落，持续时间 25d左右。从流域的气象条件，下垫面条件和流量过程的分析知，该流 域降雨径流关系具有蓄满产流的特点，可以按蓄满产流建立产流量预报方水文预报课程设计 （河海 ）第三章 计算公式3.1 降雨量计算由已知资料知，该地区雨量站分布均匀，且宝口以上有四个雨量站： 禾多布、马山、高潭、宝口，其权重系数分别为： 0.33 、0.14 、0.33 、0.20 ， 所以宝口流域平均降雨量 0.33*P1+0.14*P2+0.33*P3+0.20*P4 。3.2 蒸散发量计算蒸散发计算采用三层蒸发计算模式，即：上层蒸发量：下层蒸发量： *深层蒸发量： *总蒸发量：式中：为流域蒸发能力（） ；

8、为下层土壤含水量（） ；为下层土壤含水 量（）； C 为蒸发扩散系数。三层蒸发模式按照先上层后下层的次序，具体计算为：1）当 时， 002）当，* 时， （）*03）当, C*（）C* 时， *（）04）当, C*（） 时，*（）3.3 降雨产流量计算水文预报课程设计(河海)1) a + 冬时，+ W ( 1 )12) a+时，一(一W)3.4水源划分计算通过稳渗率可划分产流中的直接径流和地下径流。次洪的各水源分量为:总地下径流：rRGfciriiPEiiPEi FCPEi FC总直接径流：rRS(PEi fc)十iPEiPEi FC3.5汇流计算根据流域净雨和流域径流单位线，采用卷积的差分形

9、式算出流域出口的流量过程。其计算公式：直接径流过程：(i)(i)*U地下径流过程：(i)*(1)+(1)*(i)*U总的流量过程：Q(i)(i)(i)第四章计算结果4.1任务一结果"w舉肛Ln15CJ35 Im.Q上上 02器89 ? £ 1 11O 069 3 O 9 -1JO I f n3 1 B _& ? B -u 叩 g g 9 g g 勺 UP.M 9SKuaeml.OS 0 L.Q3 0 L04 P0E foTFTH n nilU匸忖侶I*饰 41Q)=1 27SLT0住匸苗5&O-JQ 1S123.23经过程序优选，0.98年份实测R（）计算R

10、（）绝对误差（）相对误差（%）19871080110019.31.791988824813-10.2-1.234.2任务二结果t(i)P(i)R(i)(i)(i)(i)(i)Q(i)2004/9/23 12:001412.72.710069.969.914.12.2004/9/23 15:00192.91010.884.194.92004/9/23 18:00119.709.733.297.313018.17.2004/9/23 21:00747.41058.311116929.28.2004/9/24 0:007518.51094.312321835.2004/9/24 3:0036125.

11、11018413632038.37.2004/9/24 6:003527.5103811485292004/9/24 9:007.86.906.964515479930.29.2004/9/24 12:005619.61079916696542.41.2004/9/24 15:006731.710874177105093.92.2004/9/24 18:008982.910864188105084.83.2004/9/24 21:001373.310122019814202004/9/25 0:0047.646.836.810177020819802004/9/25 3:0056.455.34

12、5.310237021825902004/9/25 6:0050.449.339.310265022728702004/9/25 9:0016.515.45.410260023628402004/9/25 12:005.84.804.8236023726002004/9/25 15:008.37.207.2186024121002004/9/25 18:002.61.501.5125023614802004/9/25 21:000.200077422810002004/9/26 0:0000003962216172004/9/26 3:0000001682143822004/9/26 6:00

13、000072.82072802004/9/26 9:00000025.02002252004/9/26 12:0000002.701941972004/9/26 15:0000000.001881882004/9/26 18:0000000.001821822004/9/26 21:0000000.0017617623 0:002324 0:002425 0:002526 0:002627 0:002712:00 12:00 12:0012:0012:00第五章误差统计与分析5.1精度评定从计算结果可见，年产流量绝对误差均小于100，所产流量的相对误差均小于5%。精度统计表明，率定的模型参数是

14、基本合理的。但由于课设时间限制与任务要求，每位同学只用两年的资料来率定，不满足水文情报预报规范中规定：洪水预报方案要求使用样本数量不少于10年的水文气象资料，其中应包括大、中、小水各种代表性年份，并保证有足够代表性的场次洪水资料。显然0.98的结果还是存在一定问题的。5.2误差分析影响流域降雨径流过程的因素很多，利用蓄满产流新安江模型的结构 与参数能够很好反映湿润地区降雨径流过程的主要规律和特点，因而能获 得较好的精度。但是模型本身以与模型计算中有很多概化，会造成误差。 造成宝口流域产汇流计算方案误差来源主要有以下几个方面：1. 资料代表性的影响。我们每位同学只用两年实测资料来率定，资料 必不

15、能满足时期要求与代表性要求。实际操作中，应要求有 12 年以上的 连续的对未来有代表性的实测资料，其中 10 年为率定期， 2 年为检验期， 且这种资料应具备丰、平、枯水年代表性，资料系列前后一致，受人类影 响较小。所以仅用两年的资料率定势必造成一定的误差，而且误差可能还 很大。2. 量测误差。实测的降雨、蒸发、径流量等水文气象信息与河流、湖 泊、地形等下垫面信息是研制预报模型或编制洪水预报方案或进行作业预 报的主要依据，在现有站网、仪器设备、观测技术条件下，各种信息的时 空变化是难以准确反映的，加上受自然因素等客观条件影响，势必造成各 种信息的量测误差。3. 模型结构误差。在该蓄满产流新安江

16、模型中，有很多将非线性现象 概化为线性现象或者将某些随机因子近似作为确定因子描述等都会带来 误差。4. 模型参数误差。模型参数是根据输入，通过模型计算输出，再将输 出过程与实测过程进行比较，用系统识别的方法作优化调试的，上述所率 定出的模型参数可能不是最优。5. 人类活动的影响。 随着社会经济的快速发展， 人类活动的影响加剧， 流域内可能新建了一些大中型水库或其他工程措施。影响了该地区原有的 产流特点，如汛期提前。5.3 实时校正 模型计算值与实测值直接总是存在一定的误差。造成两者间误差的因水文预报课程设计 (河海 )素很多，若针对每一个单一因素是难于描述或预见的，一般采用实时校正模型来解决。

17、实时校正模型常用的有卡尔曼滤波、自回归模型等。第六章 编写程序11( )22()321()122()m = 0 40= 0.9 + 0.01 * m(2), (2), (2), (3)W0 , 0 , 0 , 00 = 10: 0 = 11: 0 = 12W0 = 0 + 0 + 0i = 1 2 = 1986 + i' 用于定位4 <> 0 100 <> 0= 365= 366'读入数据& "." && "" #1t()()()P1()P2()P3()P4()j = 1 365#1, t(

18、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 = 1W(0) = W0: (0) = 0: (0) = 0: (0) = 0P3(j) + 0.2P(j) = 0.33 * P1( j) + 0.14 * P2(j) + 0.33* P4(j)(j) = * (j)(j) = P(j) - (j)(j) > 0'蒸发(j) = (j)(j) = 0(j) = 0A ,=

19、 * (1 + B) / (1 - ) (W(j - 1) - ) < 0.01A =* (1 - (1 - W(j - 1) / ) A (1 / (1 + B)A = A + (j)A <=+* (1 - A / ) A (1R(j) = (j) + W(j - 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) =

20、 (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) =

21、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)jW0 = W(): 0 = (): 0 = (): 0 = ()j = 1- 1'年实测径流'年计算(i) = (i) + (j) * 3.6 * 24 / 553 深j = 1(i) = (i) + R(j)径流深(i) = (i) * 10 + 0.5) / 10(i) = (i)

22、 * 10 + 0.5) / 10(i) = (i) - (i) * 10 + 0.5) / 10(i) = (i) - (i) / (i) * 100 + 0.5) / 100(i) <= 0.05= 19871 ; (i)水文预报课程设计 (河海 )2 ; (i)i(i - 2) <= 0.05(i - 1) <= 0.05(i) = (i - 2) + (i - 1) / 2'(i) 为两年的平均值 3 ;(i)m1 = (" 请输入您选择的= 1") ' 使用者选用自己认为合适的(2), (2), (2), (3)0 = 10: 0

23、 = 11: 0 = 12W0 = 0 + 0 + 0i = 1 2= 1986 + i'用于定位4 <> 0100 <> 0= 365= 366'读入数据& "." && "" #1t()()()P1()P2()P3()P4()j = 1 365#1, t(j), (j), (j), P1(j), P2(j), P3(j), P4(j)j#1= 2: = 3: C = 4= 5: B = 6: = 7P(), (), ()E(), (), (), ()W(), (), (), ()R()j

24、 = 1W(0) = W0: (0) = 0: (0) = 0: (0) = 0P(j) = 0.33 * P1(j) + 0.14 * P2(j) + 0.33 * P3(j) + 0.2* P4(j)(j) = * (j)+ B)R(j) = (j) + W(j - 1) -(j - 1) + (j) - R(j) >(j) =(j - 1) + (j) - R(j) - ) + (j - 1) >'第一层(j) = P(j) - (j)(j) > 0(j) = (j)'蒸发(j) = 0(j) = 0= * (1 + B) / (1 - )(W(j -

25、1) - ) < 0.01A =A =* (1 - (1 - W(j - 1) / ) A (1 / (1 + B)A = A + (j)A <=第二层(j) =R(j) = (j) + W(j - 1) -+* (1 - A / ) A (1(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

26、) >= (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) = 0R(j)j#2'年统计j = 1- 1(i) = (i) + (j) * 3.6 * 24 / 553'年实测径流深水文预报课程设计 (河海 )(j) = (j - 1)(j) = C * (j) - (j) - (j)(j) = (j - 1) + P(j) - (

27、j)(j) = (j - 1) - (j)(j) = (j - 1) - (j)E(j) = ( j) + (j) + (j)W(j) = (j) + (j) + (j)jW0 = W(): 0 = (): 0 = (): 0 = ()'输出文件& "" && "" #2j = 1#2, t(j), P(j), (j), (j), (j), (j), (j), E(j), (j), (j), (j), W(j),j = 1'年计算(i) = (i) + R(j)径流深j(i) = (i), "0.00&q

28、uot;)(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 & ")=" & (

29、i)"相对误差 (" & i & ")=" & (i) i31()'任务二'读入数据t()()P1()P2()P3()P4()#1, t(), (), P1(), P2(), P3(), P4()#1'参数赋值= 1: = 20: = 60: C = 0.16= 140: B = 0.2:= 0.001P(), (), ()E(), (), (), ()W(), (), (), ()R()'初始状态W(0) = (2) + (3) + (4): (0) = (2): (0) = (3): (0) =

30、 (4)j = 1P(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) = 0A ,= (1 + B) * / (1 - )(W(j - 1) - ) < 0.1A =A =* (1 - (1 - W(j - 1) / ) A (1 / (1 + B)A = A + (j)A <=B)R(j) = (j) + W(j - 1) -+* (1 - A

31、 / ) A (1 +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)

32、 = (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& ""#2j = 1j), W(j), R(j)#2, t(j), P( j), ( j), (j), ( j), (j), (j), E(j), (j), (j), ( j#2() , ()(), ()= 10#i = 1 ()(i) &