水文预报课程设计文档格式.docx

1、蒸散发（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:7.61620.612.62004-9-23 18:6.414.915.92004-9-23 21:8.817.229.418.52004-9-24 0:1.22534.835.324.62004-9-24 3:0.929.929.243.937.82004-9-24 6:38.624.846.9332004-9-24 9:6.97.56.112.32004-9-24 12:28.334.228.5200

2、4-9-24 15:25.642.739.875.42004-9-24 18:93.9137.612413.22004-9-24 21:85.390.88575.92004-9-25 0:0.851.547.749.238.52004-9-25 3:1.170.342.197.72004-9-25 6:43.247.361.545.92004-9-25 9:20.513.315.813.12004-9-25 12:10.581.83.32004-9-25 15:7.48.410.92004-9-25 18:2.82.14.62004-9-25 21:0.20.32004-9-26 0:2004

3、-9-26 3:2004-9-26 6:2004-9-26 9:2004-9-26 12:2004-9-26 15:2004-9-26 18:2004-9-26 21:2（5）计算参数数据，见表2-2 表2-2 计算参数表计算年份参数初始张力水蓄量19891990WmUmLmDmWWUWLWD14020601101040BCFcIM0.16220.001（6）流域单位线单位线过程（m3/s）为：0，40，80，130，100，80，48，20，10，5，0（7）地下径流汇流Cg=0.978，Qg=55.3m3/s第三章 计算公式 该流域海洋性气候显着、气候温和、雨量丰沛，多年平均降雨量为180

4、0mm，径流系数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） WLC EU=WU+P,EL=C（EP-EU）,ED=0;4）当WU+PEP，WL0，则产流；否则不产流。产流时： 1）当PE+aWMM

5、： R=PE+W-WM+WM（1-（PE+a）/WMM）b+12）当PE+aWMM： R=PE+W-WMPE为扣除蒸发量后的降雨量，a为土壤含水量W对应的土壤水深，WM为流域平均蓄水容量，WMM为流域各地点包气带蓄水容量的最大值，b为流域包气带蓄水容量分布的不均匀指数，R为流域产流量。水源划分流域坡地上的降雨产流量因产流过程的条件和运动路径不同，受流域的调蓄作用不同，各径流成分在流量过程线上的反应是不一样的。在实际工作中，常需按各种径流成分分别计算或模拟，因为要对产流量进行水源划分。直接径流和地下径流水源划分如下：1）当PE=FC时：RS=0.0 RG=R2）当PEFC时：RG=FC*R/PE

6、 RS=R-RGFC为稳定下渗率，RS为直接径流，RG为地下径流。降雨补充土壤含水量，由前一天的土壤含水量推求第二天的土壤含水量，补充来源为降雨减去蒸散发减去径流量，顺序为上、下、深层依此补充。三层蓄水量变化的具体计算如下：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）-LMLM， WDi-EDi+WLi-ELi+（WUi+Pi-EUi-Ri-UM）-LMDM时， WDi+1=DM; Wi+1=WUi+1+WLi+1+W

7、Di+1；i表示第i天。3.2汇流计算根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。直接径流汇流可根据该流域的时段单位线推求，地下径流汇流由线性水库演算法推求。1）直接径流汇流计算公式：QS（i）=RS（i）UH；UH为该流域的单位线2）地下径流汇流计算公式：QG（i）=CGQG（i-1）+（1-CG）RG（i）UQ（i）=QS（i）+QG（i）U为单位转换系数，3.3模型参数1）Kc：蒸散发能力折算系数，它主要反映流域平均高程与蒸发站高程之间差别的影响和蒸发皿散发与路面蒸散发间差别的影响；2）WM:流域平均张力水容量，它表示流域蓄满的标准；3）WUM:上层张力水蓄水

8、容量，它包括了植物截留量；4）WLM：下层张力水蓄水容量；5）b：流域包气带蓄水容量分布的不均匀指数，在一般情况下其取值与单元流域面积有关；6）C:深层蒸散发扩散系数，它主要取决于流域内深根植物的覆盖范围。7）IM:不透水面积占全流域面积的比例，它的值可由大比例尺的地形图，通过地理信息系统现代技术量测出来，也可用历史上干旱期小洪水资料来分析。第四章 模型检验、结果评定及分析水文预报是一项直接服务于国家安全和国民经济建设的不可或缺的重要基本工作，是帮助人类有效地预防洪水、减少洪灾损失，有效利用水资源的非工程措施之一。随着经济、社会发展及其全球化进程的需要，水文预报的服务面进一步拓展，对水文预报提

9、出了更高的要求。水文预报结果的准确率与可信程度是衡量服务质量的前提，为了更好地为国家安全和国民经济建设服务，必须对水文预报结果的可靠性和有效性进行评定和检验。4.1产流模式的检验 定性分析该流域集雨面积856km2。流域地处南方，海洋性气候显着，气候温和，雨量丰沛。降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期49月降雨量占年降雨量的81左右：径流系数0.50.7。流域上游有一水文站，控制流域面积553km2，占流域集雨面积的64.6%。该水文站以上有4个雨量站。雨量站分布较均匀，有10年逐日降水资料和时段降水资料；该水文站具有10年以上水位

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

11、对误差最小时，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%以内，并使得各年的相对误差差别尽可能小。根据相对误差规律得以下三种方法，同时说明：

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

13、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时，拟合较好，此时各年相对误差如下表3.表3 Kc=1.300时的各年份相对误差值8788899091920.0741 -0.0169 0.0983 0.5440 -0.0754 0.0715 4.3.2 1987-1992年各年计算径流与实测径流的拟合结果Kc

14、取1.300时，各年逐日径流理论计算值和实测值对比图见下图3-8. 图3 1987年逐日径流理论计算值和实测值对比图图4 1988年逐日径流理论计算值和实测值对比图图5 1989年逐日径流理论计算值和实测值对比图图6 1990年逐日径流理论计算值和实测值对比图图7 1991年逐日径流理论计算值和实测值对比图图8 1992年逐日径流理论计算值和实测值对比图以上各图是在日径流理论计算时未计算基流的情况下得到的，可见其整个起伏趋势还是相当契合的，但局部仍存在差异，基本可以满足常次预报的要求。因此确定Kc取1.300。4.4误差来源设计的蓄满产流模型结构与流域的实际产流过程和规律不完全相符， 出现的问

15、题以及可能误差影响因素包括：（1）1990年数据计算出的结果误差过大，可能是因为1990年数据存在问题。（2）各年先对误差均已经很小，但不能全部达到5%的范围之内，可能是因为调试不够，或者是因为模型假设与市级的情况不尽相同。（3）逐日径流计算理论值和实测值之间存在较大的相对误差，可能是因为未对基流进行计算，模型本身精度也有限，导致这种误差.。（4）由于流域地理、气候、气象、水文条件上与模型假设条件存在一定程度上的差异，导致计算结果存在误差。4.5模型的应用-暴雨预报对2004年暴雨过程进行洪水预报。运用单位线法和出流系数法分别计算直接径流出流量和地下径流出流量，两者之和即为总的流量，其中Kc取

16、1.300。具体结果见下表1。表4 2004年暴雨过程月日时RRsRg地下径流Qg直接径流Qs总径流Q923123.70 2.73 0.97 55.30 0.00 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.87 1.57 56.11 70.38 126.49 2424.95 22.21 2.74 57.96 115.01 172.98 334.24 31.24 3.00 60.07 228.10 288.17 636.53 33.53 62.13 454.13 516.26 6.

17、04 3.04 64.14 781.11 845.25 28.75 25.75 66.11 984.05 1050.16 40.88 37.88 68.03 1095.58 1163.62 92.05 89.05 69.91 1106.62 1176.53 82.33 79.33 71.76 1495.12 1566.88 46.05 43.05 73.56 2053.13 2126.69 54.26 51.26 75.32 2673.36 2748.68 48.20 45.20 77.04 2953.15 3030.19 14.31 11.31 78.72 2909.37 2988.10 3

18、.69 0.69 80.37 2668.85 2749.23 6.16 3.16 81.98 2147.52 2229.50 0.45 80.69 1500.58 1581.27 78.91 964.96 1043.87 77.17 540.63 617.81 2675.48 254.57 330.05 73.82 122.02 195.84 72.19 50.44 122.64 70.60 12.66 83.26 69.05 3.50 72.55 67.53 1.58 69.11 66.05 64.59 2763.17 61.78 60.42 59.09 57.79 56.52 55.28 54.06 2852.87 将直接径流、地下径流出流过程和总流量过程绘制出洪水流量过程线，见下图9.图9 洪水流量过程线第五章 总结和心得此次课程设计，做了很久，期间碰到诸多棘手的问题，终于一一解决，得以完成，虽然还有些不尽人