1、水文预报课程设计.水文预报课程设计报告学 院： _水利与环境学院 _专 业： _水文与水资源工程 _班 级： 200905201姓 名： _马天玉 _学 号： _20090520115_指导教师： _胡彩虹 _,.第一章 基本任务1.1 蒸发折算系数 Kc的优选根据已给数据资料及参数（本报告采用 89-92 年的历史数据），将流域作为整体：（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数 Kc；（4）8990 年的历时数据作为率定参数， 9192 年的数据作为模型检验。1.2 暴雨预报根据已给的设计暴雨资料和任务一率定的 Kc，将流域作为整体进
2、行如下计算：（1）次洪产流量计算，划分水源；（2）直接径流汇流，地下径流汇流的计算。（3）采用 2004 年暴雨数据进行预报。根据已给的资料、 参数及做过的习题， 自己编写程序， 将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较。,.第二章 基本资料2.1 流域概况白盆珠水库位于广东省东江一级支流西枝江的上游， 坝址以上集雨面积 856 km2。流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。降雨年际间变化大，年内分配不均，多年平均降雨量为 1800mm，实测年最大降雨量为 3417mm，汛期 49 月降雨量占年降雨量的
3、 81左右：径流系数 0.5 0.7 。流域内地势平缓， 土壤主要有黄壤和砂壤， 具有明显的腐殖层， 淀积层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。2流域上游有宝口水文站，流域面积 553km，占白盆珠水库坝址以上集雨面积的 64.6%。白盆珠水库有 10 年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有 7 个雨量站，其中宝口以上有 4 个。雨量站分布较均匀有10 年逐日降水资料和时段降水资料； 宝口水文站具有 10 年以上水位、流量资料；流域属山区性小流域且受到地形、 地貌等下垫面条件影响， 洪水陡涨缓落，
4、 汇流时间一般 23h，有时更短；一次洪水总历时 2 5d。图 2-1 该地区水文站分布图,.2.2 数据资料2(1) 计算流域面积为 553km。(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（
5、 i ）2004-9-23 12:001.36.29.921.617.32004-9-23 15:001.37.61620.612.62004-9-23 18:001.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
6、 15:000.925.642.739.875.42004-9-24 18:000.993.9137.612413.22004-9-24 21:000.985.390.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.
7、11.82.82.14.6,.2004-9-25 21:001.10.200.302004-9-26 0:001.200002004-9-26 3:002.100002004-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 计算参数表计算年份 参数 初始张力水蓄量Wm Um Lm Dm W WU WL WD1989 199 140 20 60 600 B
8、 C Fc IM 110 10 40 600.2 0.16 22 0.001（ 6）流域单位线3单位线过程（ m/s ）为： 0，40，80， 130，100，80， 48，20，10， 5， 0（ 7）地下径流汇流3Cg=0.978，Qg=55.3m/s,.第三章 计算公式该流域海洋性气候显著、气候温和、雨量丰沛，多年平均降雨量为 1800mm，径流系数 0.5-0.7 ，土壤主要有黄壤和砂壤，层次结构明显，透水性好，植被覆盖度高，地势平坦，由此可初步判定该流域的产流机制为蓄满产流模式。3.1 产流计算3.1.1 蒸散发计算根据流域蓄满产流特点， 蒸散发计算采用的是三层蒸散发计算模式。 三层
9、蒸发模式的具体计算如下：1) 当 WU+PEP,EU=E p,EL=0,ED=0;2)当 WU+P EP, WL CWLM, EU=WU+P,EL=(E P-EU)WL/WLM,ED=0;3)当 WU+P EP, C(E P-EU) WL CWLM, EU=WU+P,EL=C(E P-EU),ED=0;4)当 WU+P EP， WL0，则产流；否则不产流。产流时：1) 当 PE+aWMM：b+1R=PE+W-WM+WM(1-(PE+a)/WMM)2) 当 PE+a WMM：R=PE+W-WM式中：PE为扣除蒸发量后的降雨量， a 为土壤含水量 W对应的土壤水深， WM 为流域平均蓄水容量， W
10、MM为流域各地点包气带蓄水容量的最大值， b 为流域包气带蓄水容量分布的不均匀指数， R为流域产流量。3.1.3 二水源划分流域坡地上的降雨产流量因产流过程的条件和运动路径不同， 受流域的调蓄作用不同，各径流成分在流量过程线上的反应是不一样的。 在实际工作中， 常需按各种径流成分分别计算或模拟，因为要对产流量进行水源划分。直接径流和地下径流水源划分如下：1）当 PE=FC时： RS=0.0RG=R2）当 PEFC时： RG=FC*R/PERS=R-RG式中： FC为稳定下渗率， RS为直接径流， RG为地下径流。3.1.4 各层蓄水量计算降雨补充土壤含水量， 由前一天的土壤含水量推求第二天的土
11、壤含水量， 补充来源为降雨减去蒸散发减去径流量，顺序为上、下、深层依此补充。,.三层蓄水量变化的具体计算如下：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)-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 ；式中
12、： 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为单位转换系数， U流域面积 F ( km2 )3.6 * t (h)3.3 模型参数1）Kc：蒸散发能力折算系数，它主要反映流域平均高程与蒸发站高程之间差别的影响和蒸发皿散发与路面蒸散发间
13、差别的影响；2）WM:流域平均张力水容量，它表示流域蓄满的标准；3）WUM:上层张力水蓄水容量，它包括了植物截留量；4）WLM：下层张力水蓄水容量；5）b：流域包气带蓄水容量分布的不均匀指数， 在一般情况下其取值与单元流域面积有关；6）C: 深层蒸散发扩散系数，它主要取决于流域内深根植物的覆盖范围。7）IM: 不透水面积占全流域面积的比例， 它的值可由大比例尺的地形图， 通过地理信息系统现代技术量测出来，也可用历史上干旱期小洪水资料来分析。,.第四章 模型检验、结果评定及分析水文预报是一项直接服务于国家安全和国民经济建设的不可或缺的重要基本工作，是帮助人类有效地预防洪水、 减少洪灾损失， 有效
14、利用水资源的非工程措施之一。随着经济、社会发展及其全球化进程的需要， 水文预报的服务面进一步拓展，对水文预报提出了更高的要求。水文预报结果的准确率与可信程度是衡量服务质量的前提， 为了更好地为国家安全和国民经济建设服务， 必须对水文预报结果的可靠性和有效性进行评定和检验。4.1 产流模式的检验定性分析2该流域集雨面积 856km。流域地处南方，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨， 常受热带风暴影响。 降雨年际间变化大，年内分配不均，多年平均降雨量为 1800mm，实测年最大降雨量为 3417mm，汛期 49 月降雨量占年降雨量的 81左右：径流系数 0.50.7
15、。流域内地势平缓，土壤主要有黄壤和砂壤， 具有明显的腐殖层， 淀积层和母质土等层次结构， 透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。2流域上游有一水文站，控制流域面积 553km，占流域集雨面积的 64.6%。该水文站以上有 4 个雨量站。雨量站分布较均匀，有 10 年逐日降水资料和时段降水资料；该水文站具有 10 年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响， 洪水陡涨缓落， 汇流时间一般 23h，有时更短；一次洪水总历时 25d。由流域概况可以看出，该地区属于湿润多雨地区，雨量集中，地势平坦，土壤层容易蓄满， 而
16、且有场系列的降雨资料和水位流量资料， 综合从气候条件、 土壤状况、植被组成以及洪水机制看，新安江模型适用于该地区的水文预报。4.2 率定期优选蒸发折算系数 Kc,.4.2.1 计算年径流与实测年径流的比较结果及误差分析采用试算法，设置 Kc 取值在 0.9-4.9 之间，步长为 0.001 ，分别用 1987-1992年的资料数据进行计算， 得到各年全年径流量理论计算值和实测值的相对误差值随着 Kc 的变化过程如下图 1：图 1 各年全年径流量理论计算值和实测值的相对误差值随着 Kc 的变化自左至右分别为 1991、1988、1992、1989、1987、1990 年，各年相对误差最小时， K
17、c 取值如下表：表 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%以内，并使得各年的相对误差差别尽可能小。根据相对误差规律得以下三种方法
18、，同时说明：通过人为观察、比较，人工选取所给 Kc 的最优值；计算 1987-1989 年三年的相对误差累计值， 相对误差累计值在一定程度上反映了 Kc 对计算径流量造成的影响，相对误差累计值越小越好；计算 1987-1989 年三年相对误差的均方差， 它在一定程度上说明了各相对误差偏离平均相对误差的程度，也反映了 Kc 取值造成误差的稳定程度，相对误差的均方差越小越好。4.2.3 确定 Kc 的取值：缩小 Kc 的取值范围，取 Kc=0.9-1.3 ，步长为 0.001 进行计算，得到各年全年径流量理论计算值和实测值的相对误差值随着 Kc 的变化过程如下图 2.图 2 各年全年径流量理论计算
19、值和实测值的相对误差值随着 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之间进行变化绘
20、制日径流量理论计算值和实测值的对比，得到 Kc 取值为 1.300时，拟合较好，此时各年相对误差如下表3.表 3Kc=1.300 时的各年份相对误差值年份878889909192相对误0.0741-0.01690.09830.5440-0.07540.0715差4.3.2 1987-1992 年各年计算径流与实测径流的拟合结果Kc 取 1.300 时，各年逐日径流理论计算值和实测值对比图见下图 3-8.,.图 3 1987 年逐日径流理论计算值和实测值对比图图 4 1988 年逐日径流理论计算值和实测值对比图图 5 1989 年逐日径流理论计算值和实测值对比图,.图 6 1990 年逐日径流理
21、论计算值和实测值对比图,.图 7 1991 年逐日径流理论计算值和实测值对比图图 8 1992 年逐日径流理论计算值和实测值对比图以上各图是在日径流理论计算时未计算基流的情况下得到的， 可见其整个起伏趋势还是相当契合的， 但局部仍存在差异， 基本可以满足常次预报的要求。 因此确定 Kc 取 1.300 。4.4 误差来源设计的蓄满产流模型结构与流域的实际产流过程和规律不完全相符， 出现的问题以及可能误差影响因素包括：（1）1990 年数据计算出的结果误差过大，可能是因为 1990 年数据存在问,.题。（2）各年先对误差均已经很小，但不能全部达到 5%的范围之内，可能是因为调试不够，或者是因为模
22、型假设与市级的情况不尽相同。（3）逐日径流计算理论值和实测值之间存在较大的相对误差，可能是因为未对基流进行计算，模型本身精度也有限，导致这种误差 . 。（4）由于流域地理、气候、气象、水文条件上与模型假设条件存在一定程度上的差异，导致计算结果存在误差。4.5 模型的应用 - 暴雨预报对 2004 年暴雨过程进行洪水预报。运用单位线法和出流系数法分别计算直接径流出流量和地下径流出流量，两者之和即为总的流量，其中 Kc 取 1.300 。具体结果见下表 1。表 42004 年暴雨过程时间RRsRg地下径流直接径流总径流 Q月日时QgQs923123.702.730.9755.300.0055.30
23、154.303.181.1255.3510.9266.27183.642.351.2855.5734.5890.15218.446.871.5756.1170.38126.492424.9522.212.7457.96115.01172.9824334.2431.243.0060.07228.10288.17636.5333.533.0062.13454.13516.2696.043.043.0064.14781.11845.251228.7525.753.0066.11984.051050.161540.8837.883.0068.031095.581163.621892.0589.053.0069.911106.621176.532182.3379.333.0071.761495.121566.882446.0543.053.0073.562053.132126.6925354.2651.263.0075.322673.362748.68,.648.2045.203.0077.042953.153030.19914.3111.313.0078.722909.372988.10123.690.693.0080.372668.852749.23156.163.163.0081.982147.522229.50180.450.000.4580.69150
