《2022年水文预报课程设计 .pdf》由会员分享，可在线阅读，更多相关《2022年水文预报课程设计 .pdf（30页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、1 水文预报课程设计报告学院：_水利与环境学院 _ 专业：_水文与水资源工程 _ 班级： 200905201 姓名：_马天玉 _ 学号：_20090520115_ 指导教师： _胡彩虹 _ 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 30 页1 第一章 基本任务根据已给数据资料及参数本报告采用89-92 年的历史数据，将流域作为整体：1进行日模型产流量计算；2比较计算年径流与实测年径流；3通过误差分析，优选蒸发折算系数Kc；48990年的历时数据作为率定参数，9192 年的数据作为模型检验。根据已给的设计暴雨资料和任务一率定的Kc，
2、将流域作为整体进行如下计算：1次洪产流量计算，划分水源；2直接径流汇流，地下径流汇流的计算。3采用 2004 年暴雨数据进行预报。根据已给的资料、 参数及做过的习题， 自己编写程序， 将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 30 页2 第二章 基本资料白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。降雨年际间变化大，年内分配不均，
3、多年平均降雨量为1800mm ，实测年最大降雨量为3417mm ，汛期 49 月降雨量占年降雨量的 81左右：径流系数0.5 0.7 。流域内地势平缓， 土壤主要有黄壤和砂壤， 具有明显的腐殖层， 淀积层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的 64.6%。白盆珠水库有10 年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7 个雨量站，其中宝口以上有4 个。雨量站分布较均匀有10 年逐日降水资料和时段降水资料； 宝口水文站具有10年以上水
4、位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响， 洪水陡涨缓落， 汇流时间一般 23h，有时更短；一次洪水总历时25d。图 2-1 该地区水文站分布图精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 30 页3 资料(1) 计算流域面积为 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) 暴雨预报的数据，
5、见表2-1 表 2-1 2004 年暴雨过程数据表时间蒸散发(mm) 降雨量 mm 禾多布马山高潭宝口Ti Ei P1i P2i P3i P4i 2004-9-23 12:00 2004-9-23 15:00 16 2004-9-23 18:00 2004-9-23 21:00 2004-9-24 0:00 25 2004-9-24 3:00 2004-9-24 6:00 33 2004-9-24 9:00 2004-9-24 12:00 2004-9-24 15:00 2004-9-24 18:00 124 2004-9-24 21:00 85 2004-9-25 0:00 2004-9-2
6、5 3:00 2004-9-25 6:00 2004-9-25 9:00 2004-9-25 12:00 8 2004-9-25 15:00 2004-9-25 18:00 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 30 页4 2004-9-25 21:00 0 0 2004-9-26 0:00 0 0 0 0 2004-9-26 3:00 0 0 0 0 2004-9-26 6:00 0 0 0 0 2004-9-26 9:00 0 0 0 0 2004-9-26 12:00 0 0 0 0 2004-9-26 15:00 0
7、0 0 0 2004-9-26 18:00 0 0 0 0 2004-9-26 21:00 2 0 0 0 0 5计算参数数据，见表2-2 表 2-2 计算参数表计算年份参数初始张力水蓄量19891990 Wm Um Lm Dm W WU WL WD 140 20 60 60 110 10 40 60 B C FcIM 0.2 220.001 6流域单位线单位线过程 m3/s 为： 0，40，80，130，100，80，48，20，10，5，0 7地下径流汇流3/s 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 30 页5 第三章 计
8、算公式该流域海洋性气候显著、气候温和、雨量丰沛，多年平均降雨量为1800mm ，径流系数 0.5-0.7 ，土壤主要有黄壤和砂壤，层次结构明显，透水性好，植被覆盖度高，地势平坦，由此可初步判定该流域的产流机制为蓄满产流模式。蒸散发计算根据流域蓄满产流特点， 蒸散发计算采用的是三层蒸散发计算模式。三层蒸发模式的具体计算如下：1) 当 WU+P EP, EU=Ep,EL=0,ED=0; 2) 当 WU+P EP, WL C WLM, EU=WU+P,EL=(EP-EU) WL/WLM,ED=0; 3) 当 WU+P EP, C(EP-EU) WL C WLM, EU=WU+P,EL=C(EP-EU
9、),ED=0; 4) 当 WU+P EP，WL 0 ，则产流；否则不产流。产流时： 1)当 PE+a WMM： R=PE+W-WM+WM (1-(PE+a)/WMM)b+12)当 PE+a WMM： R=PE+W-WM 式中：PE为扣除蒸发量后的降雨量， a 为土壤含水量 W对应的土壤水深， WM为流域平均蓄水容量， WMM 为流域各地点包气带蓄水容量的最大值，b 为流域包气带蓄水容量分布的不均匀指数，R为流域产流量。水源划分流域坡地上的降雨产流量因产流过程的条件和运动路径不同，受流域的调蓄作用不同， 各径流成分在流量过程线上的反应是不一样的。在实际工作中， 常需按各种径流成分分别计算或模拟，
10、因为要对产流量进行水源划分。直接径流和地下径流水源划分如下：1当 PE=FC 时：RS RG=R 2当 PE FC时：RG=FC*R/PE RS=R-RG 式中： FC为稳定下渗率， RS为直接径流， RG为地下径流。降雨补充土壤含水量， 由前一天的土壤含水量推求第二天的土壤含水量，补充来源为降雨减去蒸散发减去径流量，顺序为上、下、深层依此补充。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 30 页7 三层蓄水量变化的具体计算如下：1WUi+Pi-EUi-RiUM，WLi-ELi+(WUi+Pi-EUi-Ri-UM)UM， WLi-E
11、Li+(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；式中： i 表示第 i 天。汇流计算根据流域净雨和流域径流单位线， 采用卷积的差分形式算出流域出口的流量过程。直接径流汇流可根据该流域的时段单位线推求，地下径流汇流由线性水库演算法推求。1直接径流汇流计算公式：QS(i)
12、=RS(i )U H；式中： UH为该流域的单位线2地下径流汇流计算公式：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 30 页8 QG(i)=CG QG(i-1)+(1- CG) RG(i) UQ(i)=QS(i)+QG(i) 式中： U为单位转换系数，)(*6 .3)(2htkmFU流域面积1Kc：蒸散发能力折算系数，它主要反映流域平均高程与蒸发站高程之间差异的影响和蒸发皿散发与路面蒸散发间差异的影响；2WM: 流域平均张力水容量，它表示流域蓄满的标准；3WUM: 上层张力水蓄水容量，它包括了植物截留量；4WLM ：下层张力水蓄水
13、容量；5b：流域包气带蓄水容量分布的不均匀指数，在一般情况下其取值与单元流域面积有关；6C:深层蒸散发扩散系数，它主要取决于流域内深根植物的覆盖范围。7IM: 不透水面积占全流域面积的比例，它的值可由大比例尺的地形图，通过地理信息系统现代技术量测出来，也可用历史上干旱期小洪水资料来分析。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 30 页9 第四章 模型检验、结果评定及分析水文预报是一项直接服务于国家安全和国民经济建设的不可或缺的重要基本工作，是帮助人类有效地预防洪水、减少洪灾损失， 有效利用水资源的非工程措施之一。 随着经济、 社
14、会发展及其全球化进程的需要，水文预报的服务面进一步拓展，对水文预报提出了更高的要求。水文预报结果的准确率与可信程度是衡量服务质量的前提，为了更好地为国家安全和国民经济建设服务， 必须对水文预报结果的可靠性和有效性进行评定和检验。产流模式的检验定性分析该流域集雨面积 856km2。流域地处南方，海洋性气候显著，气候温和，雨量丰沛。暴雨成因主要是锋面雨和台风雨， 常受热带风暴影响。 降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm ，实测年最大降雨量为3417mm ，汛期49 月降雨量占年降雨量的81左右：径流系数0.50.7 。流域内地势平缓，土壤主要有黄壤和砂壤， 具有明显的腐殖层
15、， 淀积层和母质土等层次结构，透水性好。台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。流域上游有一水文站，控制流域面积553km2，占流域集雨面积的64.6%。该水文站以上有 4 个雨量站。雨量站分布较均匀，有10 年逐日降水资料和时段降水资料；该水文站具有10 年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响， 洪水陡涨缓落， 汇流时间一般 23h，有时更短；一次洪水总历时 25d。由流域概况可以看出，该地区属于湿润多雨地区，雨量集中，地势平坦，土壤层容易蓄满， 而且有场系列的降雨资料和水位流量资料，综合从气候条件、 土壤状况、植被
16、组成以及洪水机制看，新安江模型适用于该地区的水文预报。4.2 率定期优选蒸发折算系数Kc 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 30 页10 采用试算法，设置 Kc取值在 0.9-4.9之间， 步长为 0.001， 分别用 1987-1992年的资料数据进行计算， 得到各年全年径流量理论计算值和实测值的相对误差值随着 Kc 的变化过程如下列图1：图 1 各年全年径流量理论计算值和实测值的相对误差值随着Kc的变化自左至右分别为 1991、1988、1992、1989、1987、1990 年，各年相对误差最小时， Kc取值如下表
17、：表 1 各年 Kc最正确取值年份Kc 最正确取值相对误差1987 1988 1989 1990 1991 1992 可以看出， 1990年的资料不可取，因此舍弃不用。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 30 页11 4. 蒸发折算系数 Kc 值的优选方法在流域整体的计算径流量与实测径流量的相对误差满足5% 以内的前提下，尽可能使 1987-1989 年连续三年的相对误差均在10% 以内，并使得各年的相对误差差异尽可能小。根据相对误差规律得以下三种方法，同时说明：通过人为观察、比较，人工选取所给Kc 的最优值；计算 198
18、7-1989 年三年的相对误差累计值， 相对误差累计值在一定程度上反映了 Kc 对计算径流量造成的影响，相对误差累计值越小越好；计算 1987-1989 年三年相对误差的均方差， 它在一定程度上说明了各相对误差偏离平均相对误差的程度，也反映了Kc取值造成误差的稳定程度，相对误差的均方差越小越好。确定 Kc 的取值：缩小 Kc 的取值范围，取 Kc=0.9-1.3 ，步长为 0.001 进行计算，得到各年全年径流量理论计算值和实测值的相对误差值随着Kc 的变化过程如下列图2. 图 2 各年全年径流量理论计算值和实测值的相对误差值随着Kc的变化表 2 各年 Kc最正确取值精选学习资料 - - -
19、- - - - - - 名师归纳总结 - - - - - - -第 12 页，共 30 页12 年份Kc 最正确取值相对误差1987 1988 1989 1991 1992 因为 1990 年数据舍弃，原定4 年率定，改为采用 1987-1989 年资料进行率定，得到 KC取 1.263-1.430之间数据时，误差均在5% 以内，其中， Kc 取值为1.3450 时，误差最小。模型检验及评价4.3.1 模型检验结果通过对 Kc 取值在 1.263-1.430之间进行变化绘制日径流量理论计算值和实测值的比照，得到 Kc 取值为 1.300 时，拟合较好，此时各年相对误差如下表3. 年份87 88
20、 89 90 91 92 相对误差0.0741 -0.0169 0.0983 0.5440 -0.0754 0.0715 4.3.2 1987-1992年各年计算径流与实测径流的拟合结果Kc取 1.300 时，各年逐日径流理论计算值和实测值比照图见下列图3-8. 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页，共 30 页13 图 3 1987年逐日径流理论计算值和实测值比照图图 4 1988年逐日径流理论计算值和实测值比照图图 5 1989年逐日径流理论计算值和实测值比照图精选学习资料 - - - - - - - - - 名师归纳总结
21、- - - - - - -第 14 页，共 30 页14 图 6 1990年逐日径流理论计算值和实测值比照图精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页，共 30 页15 图 7 1991年逐日径流理论计算值和实测值比照图图 8 1992年逐日径流理论计算值和实测值比照图以上各图是在日径流理论计算时未计算基流的情况下得到的，可见其整个起伏趋势还是相当契合的， 但局部仍存在差异， 基本可以满足常次预报的要求。因此确定 Kc取 1.300。误差来源设计的蓄满产流模型结构与流域的实际产流过程和规律不完全相符，出现的问题以及可能误差影响因素包括
22、：11990年数据计算出的结果误差过大，可能是因为1990年数据存在问精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 16 页，共 30 页16 题。2各年先对误差均已经很小，但不能全部到达5% 的范围之内，可能是因为调试不够，或者是因为模型假设与市级的情况不尽相同。3逐日径流计算理论值和实测值之间存在较大的相对误差，可能是因为未对基流进行计算，模型本身精度也有限，导致这种误差. 。4由于流域地理、气候、气象、水文条件上与模型假设条件存在一定程度上的差异，导致计算结果存在误差。对 2004 年暴雨过程进行洪水预报。运用单位线法和出流系数法分别计算直
23、接径流出流量和地下径流出流量，两者之和即为总的流量，其中Kc 取 1.300。具体结果见下表1。表 4 2004 年暴雨过程时间月日时R Rs Rg 地下径流Qg 直接径流Qs 总径流 Q 9 23 12 3.70 2.73 0.97 55.30 0.00 55.30 15 4.30 3.18 1.12 55.35 10.92 66.27 18 3.64 2.35 1.28 55.57 34.58 90.15 21 8.44 6.87 1.57 56.11 70.38 126.49 24 24.95 22.21 2.74 57.96 115.01 172.98 24 3 34.24 31.24
24、 3.00 60.07 228.10 288.17 6 36.53 33.53 3.00 62.13 454.13 516.26 9 6.04 3.04 3.00 64.14 781.11 845.25 12 28.75 25.75 3.00 66.11 984.05 1050.16 15 40.88 37.88 3.00 68.03 1095.58 1163.62 18 92.05 89.05 3.00 69.91 1106.62 1176.53 21 82.33 79.33 3.00 71.76 1495.12 1566.88 24 46.05 43.05 3.00 73.56 2053.
25、13 2126.69 25 3 54.26 51.26 3.00 75.32 2673.36 2748.68 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 17 页，共 30 页17 6 48.20 45.20 3.00 77.04 2953.15 3030.19 9 14.31 11.31 3.00 78.72 2909.37 2988.10 12 3.69 0.69 3.00 80.37 2668.85 2749.23 15 6.16 3.16 3.00 81.98 2147.52 2229.50 18 0.45 0.00 0.45 80.6
26、9 1500.58 1581.27 21 0.00 0.00 0.00 78.91 964.96 1043.87 24 0.00 0.00 0.00 77.17 540.63 617.81 26 3 0.00 0.00 0.00 75.48 254.57 330.05 6 0.00 0.00 0.00 73.82 122.02 195.84 9 0.00 0.00 0.00 72.19 50.44 122.64 12 0.00 0.00 0.00 70.60 12.66 83.26 15 0.00 0.00 0.00 69.05 3.50 72.55 18 0.00 0.00 0.00 67.
27、53 1.58 69.11 21 0.00 0.00 0.00 66.05 0.00 66.05 24 64.59 0.00 64.59 27 3 63.17 0.00 63.17 6 61.78 0.00 61.78 9 60.42 0.00 60.42 12 59.09 0.00 59.09 15 57.79 0.00 57.79 18 56.52 0.00 56.52 21 55.28 0.00 55.28 24 54.06 0.00 54.06 28 3 52.87 0.00 52.87 将直接径流、地下径流出流过程和总流量过程绘制出洪水流量过程线，见下列图 9. 图 9 洪水流量过程
28、线精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页，共 30 页18 第五章 总结和心得此次课程设计，做了很久，期间碰到诸多棘手的问题，终于一一解决，得以完成，虽然还有些不尽人意， 但总体上还是较好得完成了此次课程设计的各项要求。通过此次课程设计，收获良多。首先是对于水文预报这门课程，为了完成这次课设，对课本进行了深入系统的复习，尤其是与新安江模型相关的内容，使得对这门课程的掌握更加扎实牢固，理解也更加深入；其二，本次课设采用的是matlab 软件进行数据处理，这也是基于数据较多，计算分析比较复杂的缘故，利用原有的一点儿编程基础， 在这期间
29、进一步不断地学习， 对该软件的使用能力有了很大程度的增强， 这对以后的学习工作都是大有裨益的。这次课设的意义就是在于不断逼迫自己去学习更多的新东西，并把以前学过的东西进行整合贯穿，到达提升自身水平的效果； 其三，在做课程设计的过程中， 老师对我们进行了多次耐心认真的辅导，同学们之间也不断进行互相的交流，也得到了很多的乐趣。总的来说，这次课程设计受益匪浅。 这是毕业设计之前的最后一次课设，也是由书本上知识运用到实践中的一次尝试，让我们意识到自身知识的匮乏和有限，长叹“书到用时方恨少” ，同时这也是一个有趣的过程。试想，通过自己的不懈努力终于把它攻克， 把成果展现在自己眼前时的喜悦，又有什么能比得
30、上？感谢老师，感谢此次课设。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 19 页，共 30 页19 附件：年连续六年，前三年率定期和后两年检验期的计算径流量与实测径流量绝对误差和相对误差计算程序：%*% % 1987-1989 年计算结果% % % %*% clc clear ZL=load(1987 到 1989 资料.txt);%ZL 表示导入的数据WM=140;UM=20;LM=60;DM=60;B=0.3;C=0.16; IM=0.002; %流域平均张力水容量WMmm),上层张力水容量 UM mm),下土层张力水容量LMmm),深层张
31、力水容量 DM(mm), 张力水蓄水容量曲线方次B，深层蒸散发折算系数 C，%不透水面积占全流域面积的比例IM Q=ZL(:,3); %日径流量实测值Qm3/s 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 20 页，共 30 页20 E0=ZL(:,4); %日蒸发量蒸发皿实测值 mmP1=ZL(:,5); P2=ZL(:,6); P3=ZL(:,7); P4=ZL(:,8); %P1,P2,P3,P4分别表示流域四个地区的日降雨量实测值mmsumQ=0; %sumQ 表示年径流量实测值 mmsumR=zeros(4000 ，1); h=len
32、gth(ZL);% 求出数据的天数R=zeros(h,1) ; %R 表示日径流量的理论计算值EU=zeros(1,h); for i = 1:h %流域平均降雨量计算P(i) = 0.33 * P1(i) + 0.14 * P2(i) + 0.33 * P3(i) + 0.2 * P4(i); %逐日降雨深的计算sumQ = sumQ + Q(i) * 24 * 3.6 / 553; %实际测得年径流量的计算end %求逐日降雨深和实测全年径流量W(1) = 110; WU(1) = 10; WL(1) = 40; WD(1) = 60; %初始的土壤总，上，下，深层土壤含水量%流域三层蒸发
33、计算WMM = WM * (1 + B); %初始土壤含水量确实定a(1) = WMM * (1 - (1 - (W(1) / WM) (1 / (1 + B); for j=1:4000 for i = 1:h EP(i) = E0(i) * Kc(j); %流域逐日蒸散发深if WU(i) + P(i) = EP(i) EU(i) = EP(i); EL(i) = 0; ED(i) = 0; end if WU(i) + P(i) = C * LM EU(i) = WU(i) + P(i); EL(i) = (EP(i) - EU(i) * WL(i) / LM; ED(i) = 0; e
34、lseif WL(i) = C * (EP(i) - EU(i) EU(i) = WU(i) + P(i); EL(i) = (EP(i) - EU(i) * C; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 21 页，共 30 页21 ED(i) = 0; elseif WL(i) 0 %当产流时if PE(i) + a(i) = WMM R(i) = PE(i) + W(i) - WM; W(i + 1) = WM; a(i + 1) = WMM; end end %产流计算完毕if WU(i) + P(i) - EU(i) - R(i)
35、= UM WU(i + 1) = WU(i) + P(i) - EU(i) - R(i); WL(i + 1) = WL(i) - EL(i); WD(i + 1) = WD(i) - ED(i); else WU(i + 1) = UM; if WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) = LM WL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM); WD(i + 1) = WD(i) - ED(i); else WL(i + 1) = LM; if WD(i)
36、 - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) - LM = DM WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) - LM; else WD(i + 1) = DM; end end 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 22 页，共 30 页22 end %此处计算的是降雨补充土壤含水量，由前一天的土壤含水量推求第二天的土壤含水量，补充来源为降雨减去蒸
37、散发减去径流量，顺序为上、下、深依此补充。if PE(i) = EP(i) EU(i) = EP(i); EL(i) = 0; ED(i) = 0; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 24 页，共 30 页24 end if WU(i) + P(i) = C * LM EU(i) = WU(i) + P(i); EL(i) = (EP(i) - EU(i) * WL(i) / LM; ED(i) = 0; elseif WL(i) = C * (EP(i) - EU(i) EU(i) = WU(i) + P(i); EL(i) = (
38、EP(i) - EU(i) * C; ED(i) = 0; elseif WL(i) 0 %当产流时if PE(i) + a(i) = WMM R(i) = PE(i) + W(i) - WM; W(i + 1) = WM; a(i + 1) = WMM; end end %产流计算完毕if WU(i) + P(i) - EU(i) - R(i) = UM WU(i + 1) = WU(i) + P(i) - EU(i) - R(i); WL(i + 1) = WL(i) - EL(i); WD(i + 1) = WD(i) - ED(i); else WU(i + 1) = UM; if W
39、L(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) = LM WL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM); WD(i + 1) = WD(i) - ED(i); 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 25 页，共 30 页25 else WL(i + 1) = LM; if WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) -
40、LM = DM WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) - LM; else WD(i + 1) = DM; end end end %此处计算的是降雨补充土壤含水量，由前一天的土壤含水量推求第二天的土壤含水量，补充来源为降雨减去蒸散发减去径流量，顺序为上、下、深依此补充。if PE(i) = EP(i) EU(i) = EP(i); EL(i) = 0; ED(i) = 0; end if WU(i) + P(i) = C * LM 精选学习资料 - - - - - - -
41、- - 名师归纳总结 - - - - - - -第 27 页，共 30 页27 EU(i) = WU(i) + P(i); EL(i) = (EP(i) - EU(i) * WL(i) / LM; ED(i) = 0; elseif WL(i) = C * (EP(i) - EU(i) EU(i) = WU(i) + P(i); EL(i) = (EP(i) - EU(i) * C; ED(i) = 0; elseif WL(i) 0 %当产流时if PE(i) + a(i) = WMM R(i) = PE(i) + W(i) - WM; W(i + 1) = WM; a(i + 1) = W
42、MM; end end %产流计算完毕if WU(i) + P(i) - EU(i) - R(i) = UM WU(i + 1) = WU(i) + P(i) - EU(i) - R(i); WL(i + 1) = WL(i) - EL(i); WD(i + 1) = WD(i) - ED(i); else WU(i + 1) = UM; if WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) = LM WL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM); WD(i +
43、 1) = WD(i) - ED(i); else WL(i + 1) = LM; if WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) - LM = DM 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 28 页，共 30 页28 WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - UM) - LM; else WD(i + 1) = DM; end end end
44、 %此处计算的是降雨补充土壤含水量，由前一天的土壤含水量推求第二天的土壤含水量，补充来源为降雨减去蒸散发减去径流量，顺序为上、下、深依此补充。if PE(i) 0 if PE(i) = Fc/8 % Fc稳定下渗率RG(i) = R(i);% RG表示地下径流深RS(i) = 0; %RS表示地表径流深else RG(i) = Fc /8* R(i) / PE(i); RS(i) = R(i) - RG(i); end end if PE(i) = 0 R(i) = 0; RG(i) = 0; RS(i) = 0; end end dwx=0,40,80,130,100,80,48,20,10
45、,5,0; n=length(dwx); %出流系数法推求地下径流Qg(1)=55.3; for i = 2:m Qg(i) = Cg * Qg(i - 1) + (1 - Cg) * RG(i) * 553 / (3 * 3.6);%Qg是算出来的逐日地下径流过程精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 29 页，共 30 页29 end for i = m+1:(m+n-1) Qg(i) = Cg * Qg(i - 1); end %单位线法求地表径流for i=1:m for j=1:n qs(i,j)=RS(i)*dwx(j)/10; end end Qs=zeros(1,m+n-1); for k=1:(m+n-1) for i=1:m for j=1:n if k=i+j-1 Qs(k)=Qs(k)+qs(i,j);%Qs 是算出来的逐日地表径流过程end end end end %把地表径流和地下径流相加推求地下径流for i=1:(m+n-1) Q(i)=Qg(i)+Qs(i);%Q 表示逐日总的径流量end hsgc=Qg Qs Q aaa=R Rs Rg 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 30 页，共 30 页