水文预报课程设计

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

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

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

4、：()()()BCW()()()()14020600.20.160.001110104060其中蒸散发折算系数需编程优选，优选范围为0.90-1.30。优选原则为：计算的2年内每年的年径流相对误差尽可能不超过5%2.2.2任务二参数及相关实测资料0蒸散发降雨量(mm)3(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.929.243.937.860.938.624.846.93390.96.97.56.11

5、2.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.1000092.10000122.10000152.10000182.10000212000

6、0任务二产流参数为任务一中日模优选的值，为10.03h,产流参数用日模优选的值，其他参数取任务1中各组的参数，初始张力水蓄量取各组的容量值。任务二汇流参数为0.968,0=55.3m3,单位线序号1至11对应的Q(m3)分别为:0、40、80、130、100、80、48、20、10、5、02.3产流方式论证宝口流域地处我国南方湿润地区，气候暖热，雨量充沛，多年平均降雨量为18001000,年径流系数在0.50.7之间，大丁0.4;流域内土质疏松，植被良好，不易超渗；一次洪水的流量过程陡涨缓落，持续时间25d左右。从流域的气象条件，下垫面条件和流量过程的分析知，该流域降雨径流关系具有蓄满产流的特

7、点，可以按蓄满产流建立产流量预报方案。第三章计算公式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当*时,()*0当,C*()C*时，*()0当

8、,时，一(W3.4水源划分计算通过稳渗率可划分产流中的直接径流和地下径流。次洪的各水源分量为：总地下径流：RGfcLniPEiiPEiFCPEiFC总直接径流：RS(PEifQ-iPEiPEiFC3.5汇流计算根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。其计算公式：直接径流过程：(i)(i)*U地下径流过程：(i)*(1)+(1)*(i)*U总的流量过程：Q(i)(i)(i)第四章计算结果4.1任务一结果经过程序优选，0.98年份实测R()计算R()绝对误差()相对误差(%19871080110019.31.791988824813-10.2-1.234.2任务二结

9、果t(i)P(i)R(i)(i)(i)(i)(i)Q(i)2004/9/2312:001412.7r2.7101069.969.92004/9/2315:0014.112.92.91010.884.194.92004/9/2318:001119.7r09.7133.2197.311302004/9/2321:0018.717.47.41058.31111692004/9/240:0029.728.5r18.51094.3112312182004/9/243:003635.125.1101841363202004/9/246:0038.337.527.5103811485292004/9/249

10、:007.86.906.96451547992004/9/2412:0030.529.619.6107991669652004/9/2415:0042.641.731.71087417710502004/9/2418:00:93.892.9P82.9108649/2421:0084.183.373.310122019814202004/9/250:0047.646.836.810177020819802004/9/253:0056.455.3F45.31012370j21825902004/9/256:0050.449.339.310265022728702004/9

11、/259:0016.515.45.410260023628402004/9/2512:005.84.8r04.8123601237126002004/9/2515:008.37.2p07厂186024121002004/9/2518:002.61.501.5125023614802004/9/2521:000.20r007741228110002004/9/260:0000P00396r2216172004/9/263:0000001682143822004/9/266:00000072.82072802004/9/269:0000r00二25.0n2002252004/9/2612:0000

12、002.701941972004/9/2615:0000r000.001881882004/9/2618:00了00P00了0.00n182n1822004/9/2621:0000000.00176176第五章误差统计与分析5.1精度评定从计算结果可见，年产流量绝对误差均小丁100,所产流量的相对误差均小丁5%精度统计表明，率定的模型参数是基本合理的。但由丁课设时间限制与任务要求，每位同学只用两年的资料来率定，不满足水文情报预报规范中规定：洪水预报方案要求使用样本数量不少于10年的水文气象资料，其中应包括大、中、小水各种代表性年份，并保证有足够代表性的场次洪水资料。显然0.98的结果还是存在一

13、定问题的。5.2误差分析影响流域降雨径流过程的因素很多,利用蓄满产流新安7工模型的结构与参数能够很好反映湿润地区降雨径流过程的主要规律和特点，因而能获得较好的精度。但是模型本身以及模型计算中有很多概化，会造成误差。造成宝口流域产汇流计算方案误差来源主要有以下几个方面：1. 资料代表性的影响。我们每位同学只用两年实测资料来率定，资料必不能满足时期要求及代表性要求。实际操作中，应要求有12年以上的连续的对未来有代表性的实测资料，其中10年为率定期，2年为检验期，且这种资料应具备丰、平、枯水年代表性，资料系列前后一致，受人类影响较小。所以仅用两年的资料率定势必造成一定的误差，而且误差可能还很大。2.

14、 量测误差。实测的降雨、蒸发、径流量等水文气象信息及河流、湖泊、地形等下垫面信息是研制预报模型或编制洪水预报方案或进行作业预报的主要依据，在现有站网、仪器设备、观测技术条件下，各种信息的时空变化是难以准确反映的，加上受自然因素等客观条件影响，势必造成各种信息的量测误差。3. 模型结构误差。在该蓄满产流新安江模型中，有很多将非线性现象概化为线性现象或者将某些随机因子近似作为确定因子描述等都会带来误差。4. 模型参数误差。模型参数是根据输入，通过模型计算输出，再将输出过程与实测过程进行比较，用系统识别的方法作优化调试的，上述所率定出的模型参数可能不是最优。5. 人类活动的影响。随着社会经济的快速发

15、展，人类活动的影响加剧，流域内可能新建了一些大中型水库或其他工程措施。影响了该地区原有的产流特点，如汛期提前。5.3实时校正模型计算值与实测值直接总是存在一定的误差。造成两者问误差的因素很多，若针对每一个单一因素是难于描述或预见的，一般采用实时校正模型来解决实时校正模型常用的有卡尔曼滤波、自回归模型等。第六章编写程序1()22()321()122()m=040=0.9+0.01*m(2),(2),(2),(3)W0,0,0,00=10:0=11:0=12W0=0+0+0i=12=1986+i用于定位401000=365=366读入数据&.
t()()()P1()P2()P3()P4()j=

16、1365#1,t(j),(j),(j),P1(j),P2(j),P3(j),P4(j)#1参数赋值=2:=3:C=4=5:B=6:=7P()，()，()E(),(),(),()W(),(),(),()R()j=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)(j)=P(j)-(j)(j)0(j)=(j)蒸发(j)=0(j)=0A,=*(1+B)/(1-)(W(j-1)-)0.01A=A=*(1-(1-W(j-1)/)A(1/(1+B)A=A+(j)A第一层(j)=(j)=W(j-1

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

18、j)+(j)+(j)W0=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)*10+0.5)/10(i)=(i)-(i)*10+0.5)/10(i)=(i)-(i)/(i)*100+0.5)/100(i)=0.05=1987;(i);(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=11:0=12

19、W0=0+0+0i=12=1986+i用于定位401000=365=366读入数据#1&t()()()P1()P2()P3()P4()j=1365#1,t(j),(j),(j),P1(j),P2(j),P3(j),P4(j)#1=2:=3:C=4=5:B=6:=7P()，()，()E(),(),(),()W(),(),(),()R()j=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)(j)=P(j)-(j)(j)0(j)=(j)蒸发(j)=0(j)=0=*(1+B)/(1-)(W

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

21、(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)W0=W():0=():0=():0=()输出文件j=1#2,t(j),P(j),(j),(j),(j),(j),(j),E(j),(j),(j),(j),W(j),R(j)#2年统计j=1-1(i)=(i)+(j)*3.6*24/553年实测径流深j=1(i)=(i)+R(j)年计算径流深(i)=(i),0.00)(i)=(i),0.00)(i)=(i)-(i)(i)=(i)-(i)/

22、(i)*100(i)=(i),0.00)(i)=(i),0.00)”实测径流(&i&)=&(i)”计算径流(&i&)=&(i)”绝对误差(&i&)=&(i)”相对误差(&i&)=&(i)i31()任务二读入数据&.#1(1)=+1t()()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)=+(3)+(4):(0)=(2):(0)=(3):(0)=(4)j=1P(j)=(0.33*P1(j

23、)+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由上至下，逐层填蓄(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)=(

24、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)E(j)=(j)+(j)+(j)W(j)=(j)+(j)+(j)j=1#2,t(j),P(j),(j),(j),(j),(j),(j),E(j),(j),(j),(j),W(j),R(j)#2(),()(),()=10#i=1()(i)(i)=*(R(i)-(i)*)/(i)(i)0(i)=)(i)=(R(i)-(i)*)(i)=(i)*10+0.5)/