水文预报课程设计.docx

1、水文预报课程设计 水文预报课程设计报告 学 院：_水利与环境学院_专 业：_水文与水资源工程_班 级： 200905201 姓 名：_马天玉_学 号：_20090520115_ 指导教师：_胡彩虹_ 第一章 基本任务1.1蒸发折算系数Kc的优选根据已给数据资料及参数（本报告采用89-92年的历史数据），将流域作为整体：（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数Kc；（4）8990年的历时数据作为率定参数，9192年的数据作为模型检验。1.2暴雨预报根据已给的设计暴雨资料和任务一率定的Kc，将流域作为整体进行如下计算：（1）次洪产流量计算，

2、划分水源；（2）直接径流汇流，地下径流汇流的计算。（3）采用2004年暴雨数据进行预报。 根据已给的资料、参数及做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较。第二章 基本资料2.1流域概况白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期49月降雨量占年降雨量的81左右：径流系数0.50.7。流域内地势平缓，土壤主要有黄

3、壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的64.6%。白盆珠水库有10年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。雨量站分布较均匀有10年逐日降水资料和时段降水资料；宝口水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般23h，有时更短；一次洪水总历时25d。图2-1 该地区水文站分布图2.2数据资料(1)

4、计算流域面积为553km2。(2) 流域内有四个雨量站，权重系数分别为0.33、0.14、0.33、0.20。(3) 日产流模型数据，具体见资料文件名：87-92data.xls，数据格式为： T(i) Q(i) E(i) P1(i) P2(i) P3(i) P4(i)(4)暴雨预报的数据，见表2-1表2-1 2004年暴雨过程数据表时间蒸散发(mm)降雨量（mm）禾多布马山高潭宝口T（i）E（i）P1（i）P2（i）P3（i）P4（i）2004-9-23 12:001.36.29.921.617.32004-9-23 15:001.37.61620.612.62004-9-23 18:001

5、.36.26.414.915.92004-9-23 21:001.38.817.229.418.52004-9-24 0:001.22534.835.324.62004-9-24 3:000.929.929.243.937.82004-9-24 6:000.938.624.846.9332004-9-24 9:000.96.97.56.112.32004-9-24 12:000.928.329.934.228.52004-9-24 15:000.925.642.739.875.42004-9-24 18:000.993.9137.612413.22004-9-24 21:000.985.390

6、.88575.92004-9-25 0:000.851.547.749.238.52004-9-25 3:001.139.870.342.197.72004-9-25 6:001.143.247.361.545.92004-9-25 9:001.120.513.315.813.12004-9-25 12:001.110.581.83.32004-9-25 15:001.17.48.47.610.92004-9-25 18:001.11.82.82.14.62004-9-25 21:001.10.200.302004-9-26 0:001.200002004-9-26 3:002.1000020

7、04-9-26 6:002.100002004-9-26 9:002.100002004-9-26 12:002.100002004-9-26 15:002.100002004-9-26 18:002.100002004-9-26 21:0020000（5）计算参数数据，见表2-2 表2-2 计算参数表计算年份参数初始张力水蓄量19891990WmUmLmDmWWUWLWD140206060110104060BCFcIM0.20.16220.001（6）流域单位线单位线过程（m3/s）为：0，40，80，130，100，80，48，20，10，5，0（7）地下径流汇流Cg=0.978，Qg=5

8、5.3m3/s第三章 计算公式 该流域海洋性气候显著、气候温和、雨量丰沛，多年平均降雨量为1800mm，径流系数0.5-0.7，土壤主要有黄壤和砂壤，层次结构明显，透水性好，植被覆盖度高，地势平坦，由此可初步判定该流域的产流机制为蓄满产流模式。3.1产流计算 3.1.1蒸散发计算根据流域蓄满产流特点，蒸散发计算采用的是三层蒸散发计算模式。三层蒸发模式的具体计算如下：1)当WU+PEP, EU=Ep,EL=0,ED=0;2)当WU+PEP, WLCWLM, EU=WU+P,EL=(EP-EU)WL/WLM,ED=0;3)当WU+PEP, C(EP-EU) WLCWLM, EU=WU+P,EL=C

9、(EP-EU),ED=0;4)当WU+PEP，WL0，则产流；否则不产流。产流时： 1)当PE+aWMM： R=PE+W-WM+WM(1-(PE+a)/WMM)b+12)当PE+aWMM： R=PE+W-WM式中：PE为扣除蒸发量后的降雨量，a为土壤含水量W对应的土壤水深，WM为流域平均蓄水容量，WMM为流域各地点包气带蓄水容量的最大值，b为流域包气带蓄水容量分布的不均匀指数，R为流域产流量。3.1.3二水源划分流域坡地上的降雨产流量因产流过程的条件和运动路径不同，受流域的调蓄作用不同，各径流成分在流量过程线上的反应是不一样的。在实际工作中，常需按各种径流成分分别计算或模拟，因为要对产流量进行

10、水源划分。直接径流和地下径流水源划分如下：1）当PE=FC时：RS=0.0 RG=R2）当PEFC时：RG=FC*R/PE RS=R-RG式中：FC为稳定下渗率，RS为直接径流，RG为地下径流。3.1.4各层蓄水量计算降雨补充土壤含水量，由前一天的土壤含水量推求第二天的土壤含水量，补充来源为降雨减去蒸散发减去径流量，顺序为上、下、深层依此补充。三层蓄水量变化的具体计算如下：1）WUi+Pi-EUi-RiUM，WLi-ELi+(WUi+Pi-EUi-Ri-UM)UM，WLi-ELi+(WUi+Pi-EUi-Ri-UM)LM， WDi-EDi+WLi-ELi+(WUi+Pi-EUi-Ri-UM)-

11、LMUM，WLi-ELi+(WUi+Pi-EUi-Ri-UM)LM， WDi-EDi+WLi-ELi+(WUi+Pi-EUi-Ri-UM)-LMDM时， WUi+1=UM; WLi+1=LM; WDi+1=DM; Wi+1=WUi+1+WLi+1+WDi+1；式中：i表示第i天。3.2汇流计算根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。直接径流汇流可根据该流域的时段单位线推求，地下径流汇流由线性水库演算法推求。1）直接径流汇流计算公式：QS(i)=RS(i)UH；式中：UH为该流域的单位线2）地下径流汇流计算公式：QG(i)=CGQG(i-1)+(1-CG)RG(i

12、)UQ(i)=QS(i)+QG(i)式中：U为单位转换系数，3.3模型参数1）Kc：蒸散发能力折算系数，它主要反映流域平均高程与蒸发站高程之间差别的影响和蒸发皿散发与路面蒸散发间差别的影响；2）WM:流域平均张力水容量，它表示流域蓄满的标准；3）WUM:上层张力水蓄水容量，它包括了植物截留量；4）WLM：下层张力水蓄水容量；5）b：流域包气带蓄水容量分布的不均匀指数，在一般情况下其取值与单元流域面积有关；6）C:深层蒸散发扩散系数，它主要取决于流域内深根植物的覆盖范围。7）IM:不透水面积占全流域面积的比例，它的值可由大比例尺的地形图，通过地理信息系统现代技术量测出来，也可用历史上干旱期小洪水

13、资料来分析。第四章 模型检验、结果评定及分析水文预报是一项直接服务于国家安全和国民经济建设的不可或缺的重要基本工作，是帮助人类有效地预防洪水、减少洪灾损失，有效利用水资源的非工程措施之一。随着经济、社会发展及其全球化进程的需要，水文预报的服务面进一步拓展，对水文预报提出了更高的要求。水文预报结果的准确率与可信程度是衡量服务质量的前提，为了更好地为国家安全和国民经济建设服务，必须对水文预报结果的可靠性和有效性进行评定和检验。4.1产流模式的检验 定性分析该流域集雨面积856km2。流域地处南方，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。降雨年际间变化大，

14、年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期49月降雨量占年降雨量的81左右：径流系数0.50.7。流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。流域上游有一水文站，控制流域面积553km2，占流域集雨面积的64.6%。该水文站以上有4个雨量站。雨量站分布较均匀，有10年逐日降水资料和时段降水资料；该水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般23h，有时更短；一

15、次洪水总历时25d。由流域概况可以看出，该地区属于湿润多雨地区，雨量集中，地势平坦，土壤层容易蓄满，而且有场系列的降雨资料和水位流量资料，综合从气候条件、土壤状况、植被组成以及洪水机制看，新安江模型适用于该地区的水文预报。4.2率定期优选蒸发折算系数Kc4.2.1计算年径流与实测年径流的比较结果及误差分析采用试算法，设置Kc取值在0.9-4.9之间，步长为0.001，分别用1987-1992年的资料数据进行计算，得到各年全年径流量理论计算值和实测值的相对误差值随着Kc的变化过程如下图1：图1 各年全年径流量理论计算值和实测值的相对误差值随着Kc的变化自左至右分别为1991、1988、1992、

16、1989、1987、1990年，各年相对误差最小时，Kc取值如下表：表1 各年Kc最佳取值年份Kc最佳取值相对误差19871.4830.0001333619881.2270.00029593519891.4235.32E-0519902.5650.00013631819911.0820.00010660919921.4180.000149347可以看出，1990年的资料不可取，因此舍弃不用。4.2.2 蒸发折算系数Kc值的优选方法在流域整体的计算径流量与实测径流量的相对误差满足5%以内的前提下，尽可能使1987-1989年连续三年的相对误差均在10%以内，并使得各年的相对误差差别尽可能小。根据

17、相对误差规律得以下三种方法，同时说明：通过人为观察、比较，人工选取所给Kc的最优值；计算1987-1989年三年的相对误差累计值，相对误差累计值在一定程度上反映了Kc对计算径流量造成的影响，相对误差累计值越小越好；计算1987-1989年三年相对误差的均方差，它在一定程度上说明了各相对误差偏离平均相对误差的程度，也反映了Kc取值造成误差的稳定程度，相对误差的均方差越小越好。4.2.3确定Kc的取值：缩小Kc的取值范围，取Kc=0.9-1.3，步长为0.001进行计算，得到各年全年径流量理论计算值和实测值的相对误差值随着Kc的变化过程如下图2.图2 各年全年径流量理论计算值和实测值的相对误差值随

18、着Kc的变化表2 各年Kc最佳取值年份Kc最佳取值相对误差19871.30.07413011719881.2270.00029593519891.30.06879280519911.0820.00010660919921.30.049825263因为1990年数据舍弃，原定4年率定，改为采用1987-1989年资料进行率定，得到KC取1.263-1.430之间数据时，误差均在5%以内，其中，Kc取值为1.3450时，误差最小。4.3 模型检验及评价4.3.1 模型检验结果通过对Kc取值在1.263-1.430之间进行变化绘制日径流量理论计算值和实测值的对比，得到Kc取值为1.300时，拟合较好

19、，此时各年相对误差如下表3.表3 Kc=1.300时的各年份相对误差值年份878889909192相对误差0.0741 -0.0169 0.0983 0.5440 -0.0754 0.0715 4.3.2 1987-1992年各年计算径流与实测径流的拟合结果Kc取1.300时，各年逐日径流理论计算值和实测值对比图见下图3-8. 图3 1987年逐日径流理论计算值和实测值对比图图4 1988年逐日径流理论计算值和实测值对比图图5 1989年逐日径流理论计算值和实测值对比图图6 1990年逐日径流理论计算值和实测值对比图图7 1991年逐日径流理论计算值和实测值对比图图8 1992年逐日径流理论计

20、算值和实测值对比图以上各图是在日径流理论计算时未计算基流的情况下得到的，可见其整个起伏趋势还是相当契合的，但局部仍存在差异，基本可以满足常次预报的要求。因此确定Kc取1.300。4.4误差来源设计的蓄满产流模型结构与流域的实际产流过程和规律不完全相符， 出现的问题以及可能误差影响因素包括：（1）1990年数据计算出的结果误差过大，可能是因为1990年数据存在问题。（2）各年先对误差均已经很小，但不能全部达到5%的范围之内，可能是因为调试不够，或者是因为模型假设与市级的情况不尽相同。（3）逐日径流计算理论值和实测值之间存在较大的相对误差，可能是因为未对基流进行计算，模型本身精度也有限，导致这种误

21、差.。（4）由于流域地理、气候、气象、水文条件上与模型假设条件存在一定程度上的差异，导致计算结果存在误差。4.5模型的应用-暴雨预报对2004年暴雨过程进行洪水预报。运用单位线法和出流系数法分别计算直接径流出流量和地下径流出流量，两者之和即为总的流量，其中Kc取1.300。具体结果见下表1。表4 2004年暴雨过程时间月日时RRsRg地下径流Qg直接径流Qs总径流Q923123.70 2.73 0.97 55.30 0.00 55.30 154.30 3.18 1.12 55.35 10.92 66.27 183.64 2.35 1.28 55.57 34.58 90.15 218.44 6.

22、87 1.57 56.11 70.38 126.49 2424.95 22.21 2.74 57.96 115.01 172.98 24334.24 31.24 3.00 60.07 228.10 288.17 636.53 33.53 3.00 62.13 454.13 516.26 96.04 3.04 3.00 64.14 781.11 845.25 1228.75 25.75 3.00 66.11 984.05 1050.16 1540.88 37.88 3.00 68.03 1095.58 1163.62 1892.05 89.05 3.00 69.91 1106.62 1176.5

23、3 2182.33 79.33 3.00 71.76 1495.12 1566.88 2446.05 43.05 3.00 73.56 2053.13 2126.69 25354.26 51.26 3.00 75.32 2673.36 2748.68 648.20 45.20 3.00 77.04 2953.15 3030.19 914.31 11.31 3.00 78.72 2909.37 2988.10 123.69 0.69 3.00 80.37 2668.85 2749.23 156.16 3.16 3.00 81.98 2147.52 2229.50 180.45 0.00 0.45

24、 80.69 1500.58 1581.27 210.00 0.00 0.00 78.91 964.96 1043.87 240.00 0.00 0.00 77.17 540.63 617.81 2630.00 0.00 0.00 75.48 254.57 330.05 60.00 0.00 0.00 73.82 122.02 195.84 90.00 0.00 0.00 72.19 50.44 122.64 120.00 0.00 0.00 70.60 12.66 83.26 150.00 0.00 0.00 69.05 3.50 72.55 180.00 0.00 0.00 67.53 1.58 69.11 210.00 0.00 0.00 66.05 0.00 66.05 2464.59 0.00 64.59 27363.17 0.00 63.17 661.78 0.00 61.78 960.42 0.00 60.42 1259.09 0.00 59.09 1557.79 0.00 57.79 1856.52 0.00 56.