水文预报课程设计报告
水文预报课程设计
《水文预报》课程设计报告学院：_____水利与环境学院_____专业：____水文与水资源工程____班级： 200905201姓名：________马天玉__________学号：______20090520115___指导教师：________胡彩虹________第一章基本任务1.1蒸发折算系数Kc的优选根据已给数据资料及参数（本报告采用89-92年的历史数据），将流域作为整体：（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数Kc；（4）89～90年的历时数据作为率定参数，91～92年的数据作为模型检验。
1.2暴雨预报根据已给的设计暴雨资料和任务一率定的Kc，将流域作为整体进行如下计算：（1）次洪产流量计算，划分水源；（2）直接径流汇流，地下径流汇流的计算。
（3）采用2004年暴雨数据进行预报。
根据已给的资料、参数及做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较。
第二章基本资料2.1流域概况白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。
流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。
暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。
降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期4～9月降雨量占年降雨量的81％左右：径流系数0.5～0.7。
流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。
台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。
流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的64.6%。
白盆珠水库有10年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。
雨量站分布较均匀．有10年逐日降水资料和时段降水资料；宝口水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般2—3h，有时更短；一次洪水总历时2～5d。
水文预报课程设计洪水预报
水文预报课程设计洪水预报一、课程目标知识目标：1. 学生能理解水文预报的基本概念，掌握洪水形成的原因及其发展过程。
2. 学生能够掌握洪水预报的主要方法及其适用条件，如降雨径流模型、统计模型等。
3. 学生能够了解我国洪水预报的现状及发展趋势，了解相关法规政策。
技能目标：1. 学生能够运用所学知识，分析特定流域的洪水形成过程，进行简单的洪水预报。
2. 学生能够运用水文预报软件，进行数据收集、处理和分析，提高解决实际问题的能力。
3. 学生能够通过小组合作，进行洪水预报案例的研究，提高沟通协作能力。
情感态度价值观目标：1. 学生能够认识到水文预报在防洪减灾中的重要作用，增强社会责任感和使命感。
2. 学生在学习过程中，培养严谨的科学态度，树立正确的价值观。
3. 学生通过了解我国在水文预报方面的成就，增强民族自豪感，激发为国家和人民服务的情怀。
本课程针对高中年级学生，结合学科特点，注重理论知识与实践操作的相结合，培养学生的实际应用能力。
课程设计以学生为中心，充分考虑学生的认知水平、兴趣和需求，采用案例教学、小组合作等方法，激发学生的学习兴趣，提高教学效果。
通过本课程的学习，使学生能够掌握洪水预报的基本知识和技能，提高防洪减灾意识，为我国水文预报事业贡献力量。
二、教学内容1. 水文预报基本概念：洪水定义、洪水分类、洪水周期与洪水频率。
2. 洪水形成原因及发展过程：降水过程、流域特性、径流形成与汇集。
3. 洪水预报方法：- 降雨径流模型：水箱模型、单位线法、流域水文模型。
- 统计模型：时间序列分析、回归分析、人工神经网络。
4. 洪水预报软件应用：数据收集、处理、分析和预报结果输出。
5. 我国洪水预报现状与发展趋势：技术进展、政策法规、防洪减灾体系。
6. 实践案例：分析特定流域洪水预报实例，掌握预报流程和操作方法。
教学内容依据课程目标，参照教材相关章节进行组织。
教学大纲安排如下：第一周：水文预报基本概念及洪水形成原因。
水文预报课程设计
水文预报课程设计《水文预报》课程设计报告学院：_____水利与环境学院_____专业：____水文与水资源工程____班级：姓名：________马天玉__________学号：_________指导教师：________胡彩虹________第一章基本任务1.1蒸发折算系数Kc的优选根据已给数据资料及参数（本报告采用89-92年的历史数据），将流域作为整体：（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数Kc；（4）89～90年的历时数据作为率定参数，91～92年的数据作为模型检验。
1.2暴雨预报根据已给的设计暴雨资料和任务一率定的Kc，将流域作为整体进行如下计算：（1）次洪产流量计算，划分水源；（2）直接径流汇流，地下径流汇流的计算。
（3）采用2004年暴雨数据进行预报。
根据已给的资料、参数及做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较。
第二章基本资料2.1流域概况白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。
流域地处粤东沿海的西部，海洋性气候显着，气候温和，雨量丰沛。
暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。
降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期4～9月降雨量占年降雨量的81％左右：径流系数0.5～0.7。
流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。
台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。
流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的64.6%。
白盆珠水库有10年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。
雨量站分布较均匀．有10年逐日降水资料和时段降水资料；宝口水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般2—3h，有时更短；一次洪水总历时2～5d。
水文预报课设报告
基本任务任务一：根据已给数据资料及参数（本报告采用89-90年的历史数据），将流域作为整体（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数Kc。
任务二：根据已给的设计暴雨资料（宝口流域23日至26日暴雨）和任务一率定的Kc，将流域作为整体进行（1）次洪产流量，划分水源.直接径流汇流，地下径流汇流计算；（2）绘出直接径流过程，地下径流过程，总的流量过程。
基本资料白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。
流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。
暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。
降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期4—9月降雨量占年降雨量的81％左右：径流系数0.5一0.7。
流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。
台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。
流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的64.6%。
白盆珠水库有10年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。
雨量站分布较均匀．有10年逐日降水资料和时段降水资料；宝口水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般2—3h，有时更短；一次洪水总历时2~5d。
计算公式（1）蒸散发计算根据流域特点，蒸散发计算采用的是三层蒸散发计算模式。
三层蒸发模式的具体计算如下：1)当WU+P>=E P,EU=E p,EL=0,ED=0;2)当WU+P<E P, WL>=C.WLM,EU=WU+P,EL=(E P-EU)*WL/WLM,ED=0;3)当WU+P<E P, C.(E P-EU)<=WL<C.WLM,EU=WU+P,EL=C*(E P-EU),ED=0;4)当WU+P<E P, WL<C*(E P-EU),EU=WU+P,EL=WL,ED=C*(E P-EU)-EL.（2）产流量计算：根据流域特点，产流量计算系根据蓄满产流理论得出的。
水文与水资源水文预报课设报告
《水文预报》课程综合实训任务书一、实训任务本次实训的任务主要包括以下两个方面：（一）编制伊河栾川以上流域降雨径流量预报方案（二）编制伊河栾川以上流域降雨径流过程预报方案二、实训目的水文预报课程是水文水资源专业的核心技能课程，通过水文预报综合实训，能进一步掌握水文预报的基本理论和基本方法，培养依据设计任务和所给资料，运用所学知识和规范，编制水文预报方案的专业技能，同时提高技术报告的编写能力和计算机的应用能力。
在业务组织能力和实际工作能力方面得到全面锻炼，为毕业后直接顶岗参与水文预报工作奠定基础。
三、流域概况（一）流域自然地理概况：伊河栾川水文站建于东径110036'北纬33047‘，是伊河上游的第一个水文站，集水面积340平方公里，河长36.9公里，流域内属石山、林区、纵坡较陡、森林茂密，植被率在50%左右，土层覆盖较薄，表层土为壤土，下层为沙卵石，下渗能力较大，河道比降陡，河网密度大，河系呈扇形。
伊河流域的气候变化，完全受季风的支配，每年的9月份以后季风逐渐南退，地方冷空气不断南下，则为北方冷高压所控制造成季风寒冷，雨雪稀少的干燥气候，6月份以后，夏季风开始活跃，南方暖湿空气不断袭来，形成了夏季炎热，而潮湿的多雨气候。
流域内各站多年平均降水量700～800毫米，6～10月为汛期，汛期雨量占全年的60%以上，7～8月为暴雨期，暴雨历时短，强度大，一般降雨历时一至两天，超过三天以上很少出现，暴雨中心多在栾川、陶湾。
伊河洪水为暴雨所形成，受降雨特性影响，洪水过程陡涨陡落，洪水多发生在 7～10月份，特大洪水发生在7～8月份。
流域内水系和站网分布见下图。
（二）流域内水系和站网分布：四、采用资料(一) 编制方案采用的资料有：1．流域出口断面逐时流量；2．流域内雨量站逐日降水量；3．流域内雨量站逐时降水量；4．实测蒸发皿日蒸发量。
五、产流方案的编制12.00 15.3 4 15.55 62.2020.00 12.4 8 13.85 110.8023.00 11.7 3 12.05 36.151964年1978年1966年R′75.25 R′0.00 R′0.00 R初13.62 R初 5.02 R初0.16 R末8.54 R末 3.63 R末 5.15 R 70.17 R -1.39 R 4.99产流量计算表中用到了以下几个公式：①时，蒸发在上层进行，② 0<时，上、下层都有蒸发③=0时，蒸发在下层进行，Wm=100mmWum=20mmWlm=80mmb=0.25Ke=0.96表格如下图所示六、汇流方案的编制汇流方案的步骤：1、选择洪水。
水文预报课程设计报告书
水文预报课程设计报告书一、设计目的一、流域水文模型的用途：洪水预报方案是现代实时洪水预报调度系统的核心部份，是提高预报精度和增加预见期的关键技术。
对水资源可持续利用：流域水文模型是水资源评判、开发、利用和治理的理论基础。
对水环境和生态系统爱惜：流域水文模型是构建面污染模型和生态评判模型的要紧平台。
流域水文模型仍是分析研究气候转变和人类活动对洪水、水资源和水环境阻碍的有效工具。
本次课程设计的目的是通过一个具体的降雨～径流预报方案的制作，使学生了解生产单位对预报任务的要求。
二、通过课程设计，要求把握如下内容：（1）流域综合退水曲线、地下水退水曲线的制作和次洪分割方式；（2）熟悉降雨～径流预报方案编制的完整进程。
（3）新安江两水源模型结构及产流参数率定方式。
（4）流域综合单位线的推求方式。
（5）洪水预报方案精度评定方式（6）利用水文模型进行实时洪水预报方式。
二、设计大体资料古亭水是上饶江的干流，源于湖南省汝城县广山的东南麓，麟潭站以上集水面积为1067km2，干流河长约93 km，地貌属于低山丘陵区，以低山散布较广，相对高差多在500 km 左右，最上游尚有部份原始丛林，以松、衫、竹类为主，灌木次之。
山涧盆地种植农作物，植被总覆盖面积占全流域的80%以上，土壤多为红色沙壤土，间有亚粘土层。
山坡残积坡坡积一样厚度为1~2米，最厚约为4~5米，在山麓坡积层与基岩接触面上，或河流下切至接触部位时，常有泉水流出，涌水量较大每秒可达数升。
本流域气候温湿，年降雨量在1700毫米左右，汛期为4~9月，降雨量占全年的73%左右，冬季有降雪，量很少。
地下水位较高，且季节变幅小，因此，一样情形下，土壤含水量较大。
洪水流量进程线极不对称，涨洪历时仅数小时至十多小时，而退水历时竟达十余日，洪水结尾的底水明显举高，说明洪水期潜水、壤中流补给十分丰硕。
说明本流域以蓄满产流为主。
水源划分暂考虑为两水源。
三、课程设计资料1、6770一、73628号洪水的流量进程，相应的时段雨量表2、6770一、73628号洪水的前期日降雨量。
水文预报课程设计
《水文预报》课程设计报告学院专业班级姓名学号负责教师2019 年 12 月 29日资料：1.沅水下游沅陵至王家河河段，河长112公里，沅陵以上流域面积为76400平方公里，王家河以上流域面积为80500平方公里，区间面积为4100平方公，约占总面积的5%，流域内雨量站网的分布见图(一)。
2.两次洪水(一九六九年五月二十五日至三十一日，一九七〇年九月二十三日至二十六日)的上下游实测流量过程，区间面积上的降雨过程和小河流量过程，分别见表(一)、(二)、(三)。
工作步骤：1．用试算法分析七〇年九月份洪水的x ，K 值。
a ．根据实测流量过程点绘上下游流量过程线。
（见图1）b.根据河段水量平衡方程计算区q ，计算时段长根据报汛要求选用3t =∆小时。
（见表1）c.计算各时段内的河槽蓄量W ∆及W ，假定不同的x 值，据公式x 10x 10+-'）（，计算出相应的0'，当其一值能使0'~W 绳套曲线基本合一时，x 即为所求，则K=W ∆/0'∆。
（见图2、3、4）d.点绘0'~W 关系曲线并分析其特点。
（见图2、3、4）2.用表（四）所写出的参数，对六九年五月份共水进行河段洪水演算。
a ．用整段演算方法进行校核预报。
(1) 根据已知的10000=M Q 时的河段演算参数x ，K 值及t ∆计算0C ，1C ，2C 值并填入表(四)中。
（见表2）(2) 根据已知的0C ，1C ，2C 值建立演算方程1211202O C I C I C O ++=，并列表进行流量演算。
（见表3）b ．取t ∆=E K ，用分段演算法进行校核预报。
(1) 据公式N K =E K ，22x 1N 21x )(--=及10000=M Q 时的x ,K 值，计算E X 、E K 值，并填入表（四）中。
（见表2）(2) 计算单元间段10000=M Q 时的0C ，1C ，2C 值填入表（四）中。
水文预报课程设计报告材料
水文预报课程设计报告学院：资源与环境学院专业：水文与水资源工程年级：2012级姓名：燨元学号：2012215334指导教师：余倩2015 年 6 月 9 日水文预报（hydrologic forecasting)就是据已知的信息对未来一定时期的水温状态作出定性或定量的预测。
已知信息，广义上指对预报水文状态有影响的一切信息，最常用的是水文与气象要素信息，如降水、蒸发、流量、水位、气温和含沙量等观测信息。
预报的水文状态变量可以是任一水文要素也可以是水文特征量，不同的状态量预报要求的已知信息不同、预报方法不同、预见期也不同。
目前通常预报的水文要素有流量、水位、冰情和旱情等。
水文预报方法以水文基本规律、水文模型研究为基础，结合生产实际问题的需要，构成具体的预报方法或预报方案，服务于生产实际。
一般水文预报研究的重点和关键有两部分：①共性规律研究，即具有一定普遍性的水文基本规律模拟方法和流域水文模型研究；②个性问题研究，对反映具体问题的特性、方法进行了解，构成具有解决各种具体实际问题的、具有较高预报精度的预报方案。
1、新安江水文模型简介 (4)1.1 概述 (4)1.2 新安江模型的基本原理 (5)1.3 新安江模型的结构 (6)1.4新安江模型的参数 (8)1.4.1参数的物理意义 (8)1.4.2 模型参数率定 (9)2、新安江日模、次模调参成果图 (11)2.1 日模模拟结果 (11)2.2次模模拟结果 (14)3、新安江水库日模型、次模型模拟结果及精度统计表 (18)3.1 新安江水库日模型 (18)4、心得体会 (18)导师评语1、新安江水文模型简介1.1 概述流域水文模型可分为物理模型、概念性模型和系统模型。
在水文预报中，概念性模型和系统模型应用较多，此处主要介绍概念性流域水文模型。
概念性流域水文模型属于数学模型，它与物理模型相比，具有许多优点：一是它的所有条件均可由原型观测资料直接给出，不受比尺的限制，即数学模型无相似律问题；二是它的边界条件及其它条件可严格控制，也可随时按实际需要改变；三是它的通用型较强，只要研制出一种适用的应用软件，就可用来解决不同的实际问题；四是它具有理想的抗干扰性能，只要条件不变，重复模拟可以得到相同的结果，不会因人、因地而异；五是它的研制费用相对较低。
水文预报课程设计报告
水文预报课程设计澴河流域洪水预报方案编制姓名：韩夕林学号：2014********组号：15目录一、设计目的 (2)二、课程设计任务 (2)任务一、A 站降雨径流相关图制作与预报 (2)1 流域自然地理特征 (2)2 水文气象特征 (2)3 预报站基本情况 (2)4 水利工程概况 (2)5 预报方案编制说明 (2)5.1 流域特征提取 (2)5.2 水文数据整理 (2)5.2.1 雨量站权重 (2)5.2.2 原始流量资料格式处理 (2)5.2.3 历年最枯流量平均值 (2)5.2.4 流域平均蓄水容量WM (2)5.2.5 流域月蒸发能力Em (2)5.2.6 消退系数K (2)5.2.7 退水曲线 (2)5.2.8 蓄泄关系曲线 (2)5.3 产流分析 (2)5.3.1 产流模式 (2)5.3.2 洪峰流量频率分析 (2)5.3.3 产流量计算 (2)5.4 水源划分 (2)5.4.1 两场洪水分析Fc (2)5.4.2 Fc分析 (2)5.5 汇流分析 (2)5.5.1 单位线 (2)5.5.2 线性水库 (2)5.5.3 水位流量关系图 (2)6 方案的检验、评定 (2)6.1 产流方案检验评定 (2)6.2 汇流方案检验评定 (2)7 预报方案的应用和注意事项 (2)8 方案作业预报应用 (2)任务二、A 站与H 站相应水位流量法预报 (2)1 预报站基本情况 (2)2 水利工程概况 (2)3 预报方案编制说明 (2)3.1 基本资料引用 (2)3.2 采用方法及有关处理技术 (2)4 方案的检验、评定 (2)5 方案作业预报应用 (2)三、课程设计的总结与心得 (2)附录 (33)一、设计目的流域洪水预报方案的用途：洪水预报方案是现代实时洪水预报调度系统的核心部分，是提高预报精度和增长预见期的关键技术。
对水资源可持续利用：流域水文模型是水资源评价、开发、利用和管理的理论基础。
对水环境和生态系统保护：流域水文模型是构建面污染模型和生态评价模型的主要平台。
水文预报实习报告
（2014 级）
学生姓名： 学 号：
韦昭华 2014301580040 陈杰
指导老师： 专
业：水文与水资源工程
武汉大学水利水电学院 水文水资.......................................................................................................................................... 2 1.1 流域概况............................................................................................................................. 2 1.2 流域水文资料..................................................................................................................... 2 1 . 2 . 1 老灌河流域水文资料的整理................................................................................ 2 1.3 任务要求............................................................................................................................. 3 二．流域产流方案...............................................
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《水文预报》课程设计报告书**大学*******学院目录第一章基本任务 (2)第二章：基本资料 (3)1、流域概况 (3)2、基本资料 (4)3.计算参数 (4)第三章计算公式 (5)1、蒸散发计算 (5)2、产流量计算 (5)3、水源划分 (5)4、汇流计算 (6)第四章：基本数据 (7)第五章：结果评定及分析 (9)第六章：计算程序及说明 (12)第七章：总结和心得 (23)第一章基本任务任务一：根据已给数据资料及参数（本报告采用89-92年的历史数据），将流域作为整体：（1）进行日模型产流量计算；（2）比较计算年径流与实测年径流；（3）通过误差分析，优选蒸发折算系数Kc；（4）89～90年的历时数据作为率定参数，91～92年的数据作为模型检验。
任务二：根据已给的设计暴雨资料和任务一率定的Kc，将流域作为整体进行如下计算：（1）次洪产流量计算，划分水源；（2）直接径流汇流，地下径流汇流的计算。
（3）采用2004年暴雨数据进行预报。
根据已给的资料、参数及做过的习题，自己编写程序，将流域作为整体进行产流量计算；将计算年径流与实测年径流进行比较。
第二章：基本资料1、流域概况白盆珠水库位于广东省东江一级支流西枝江的上游，坝址以上集雨面积856 km2。
流域地处粤东沿海的西部，海洋性气候显著，气候温和，雨量丰沛。
暴雨成因主要是锋面雨和台风雨，常受热带风暴影响。
降雨年际间变化大，年内分配不均，多年平均降雨量为1800mm，实测年最大降雨量为3417mm，汛期4～9月降雨量占年降雨量的81％左右：径流系数0.5～0.7。
流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。
台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。
流域上游有宝口水文站，流域面积553km2，占白盆珠水库坝址以上集雨面积的64.6%。
白盆珠水库有10年逐日入库流量资料、逐日蒸发资料和时段入库流量资料：流域内有7个雨量站，其中宝口以上有4个。
雨量站分布较均匀．有10年逐日降水资料和时段降水资料；宝口水文站具有10年以上水位、流量资料；流域属山区性小流域且受到地形、地貌等下垫面条件影响，洪水陡涨缓落，汇流时间一般2—3h，有时更短；一次洪水总历时2～5d。
图2-1 该地区水文站分布图2、基本资料(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)3.计算参数计算参数见表2-1。
第三章计算公式1、蒸散发计算根据流域特点，蒸散发计算采用的是三层蒸散发计算模式。
三层蒸发模式的具体计算如下：1)当WU+P≥E P,EU=E p,EL=0,ED=0;2)当WU+P<E P, WL≥C×WLM,EU=WU+P,EL=(E P-EU)×WL/WLM,ED=0;3)当WU+P<E P, C.(E P-EU)≤WL<C.WLM,EU=WU+P,EL=C×(E P-EU),ED=0;4)当WU+P<E P，WL<C×(E P-EU)，EU=WU+P,EL=WL,ED=C×(E P-EU)-EL.2、产流量计算根据流域特点，产流量计算系根据蓄满产流理论得出的。
蓄满产流，即任一地点上，土壤含水量达田间持水量前，降雨量全部补充土壤含水量，不产流；当土壤蓄满后，其后续降雨量全部产生径流。
a=WMM×(1-(1-W0/WM)^(1/(b+1))PE>0，则产流；否则不产流。
产流时，当PE+a≤WMM；R=PE+W-WM+WM×(1-(PE+aR ΔW)/WMM)^(b+1)当PE+a＞WMM；R=PE+W-WM3、水源划分流域坡地上的降雨产流量因产流过程的条件和运动路径不同，受流域的调蓄作用不同，各径流成分在流量过程线上的反应是不一样的。
各水源分量如下： 1）当PE ≤FC 时：RS=0.0 RG=R 2） 当PE ＞FC 时：RG=FC*R/PE RS=R-RG 4、汇流计算根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。
计算公式： QS(I)=RS(I)×UQG(i)=CG ×QG(i-1)+(1-CG)×RG(i)×U Q(I)=QS(I)+QG(I)式中，U 为单位转换系数，)(*6.3)(2h t km F U ∆=流域面积第四章：基本数据任务一:1.计算流域为白盆珠水库上游的宝口流域，流域面积553 km2。
2.流域内有四个雨量站：禾多布、马山、高潭、宝口，权中重分别为0.55、0.14、0.33、0.20。
3. 资料数据格式为：T(i) Q(i) E(i) P1(i) P2(i) P3(i) P4(i)4.计算参数：1、计算参数：WM=140,UM=20,LM=60,DM=60，IM=0.002,FC=11.0W=140,WU=20,WL=60,WD=60b=0.3,C=0.16，Cg=0.978,Qg=55.32、暴雨过程：1989～1992年数据见文件89-92data.xls。
2004年暴雨过程数据见表4-1。
表4-1 2004年暴雨过程数据表2、单位线单位线过程（m3/s）为：0，40，80，130，100，80，48，20，10，5，0闰年按天数366，非闰年天数365。
流域蒸散发折算系数Kc优选原则：计算的2年资料的Kc应相同并使得2年内每年的年径流相对误差尽可能不超过5％。
第五章：结果评定及分析水文预报结果的准确率与可信程度是衡量服务质量的前提，为了更好地为国家安全和国民经济建设服务，必须对水文预报结果的可靠性和有效性进行评定和检验。
然而，设计的模型结构与流域的实际产流过程和规律也不完全相符，预报误差是客观存在的，故应采用水文预报中结果评定方法进行评定模型，并且分析结果中误差产生的主要原因。
1、任务一：日模型计算结果图5-1 当Kc=1.42时日模型计算结果根据1989年数据资料，当流域蒸散发折算系数Kc=1.42时，1989年数据计算与实际的相对误差最小，为0.167%。
但是此时计算得出1990年数据的相对误差达约44.8%，故不适宜取1.42。
图5-2 Kc=2.56时日模型计算结果同理，当流域蒸散发折算系数Kc=2.56时，1990年数据计算的相对误差最小，为0.148%，但此时1989年计算的相对误差约39%，故也不适宜取2.56。
图5-3 Kc=1.954时1989年、1990年日模型计算结果对比两年的资料来率定参数Kc，误差较大。
虽然误差比较大，但经过计算机程序优选，计算得出使两者误差最接近时的流域蒸散发折算系数为Kc=1.954，结果显示出1989年和1990年数据的误差均约为21.6%，选择这样的流域蒸散发折算系数能使之具有较好的稳定性。
由图5-1、图5-2可知，当Kc增加或者减小时，其中一年误差虽然在减小，但是另一年的误差却会随之增加，故取Kc=1.954时，使结果具有良好的稳定性。
当取流域蒸散发折算系数Kc=1.954时，运用1991年和1992年资料进行检验，分别计算出的相对误差为34.8%和19.4%，可知造成此误差的主要原因有：1、资料数据较少，不能适用于整体。
2、模型的构造与实际有一定的不符，造成误差。
但根据1989年和1990年资料进行率定时，取Kc=1.954时的稳定性较好，故最终选择1.954作为流域蒸散发折算系数。
任务二：洪水预报表5-1 次洪计算成果表6 0 0 0 65.4566 182.0063 247.46299 0 0 0 21.9134 178.0022 199.915612 0 0 0 1.6558 174.0861 175.741915 0 0 0 0 170.2562 170.256218 0 0 0 0 166.5106 166.510621 0 0 0 0 162.8474 162.8474运用单位线法和出流系数法分别计算直接径流出流量和地下径流出流量，两者之和即为总的流量，具体结果见上表5-1。
将直接径流、地下径流出流过程和总流量过程绘制出洪水流量过程线，见下图5-4。
图5-4 洪水流量过程线第六章：计算程序及说明Dim P(1 To 366) As Single '流域的平均降雨量Dim i As Integer, P1(1 To 366) As Single, P2(1 To 366) As Single, P3(1 To 366) As Single, P4(1 To 366) As SingleDim E0(1 To 366) As Single '流域的水面蒸发Dim Ep(1 To 366) As Single '流域的蒸发能力Dim E(1 To 366) As Single '流域总的蒸发量Dim EU(1 To 366) As Single, EL(1 To 366) As Single, ED(1 To 366) As SingleDim Q(1 To 366) As Single '流域实测径流Dim R(1 To 366) As Single '流域计算径流Dim RS(1 To 366) As Single, RG(1 To 366) As Single, PE(1 To 366) As SingleDim W(1 To 367) As Single '流域总的蓄水量Dim WU(1 To 367) As Single, WL(1 To 367) As Single, WD(1 To 367) As Single, WMM As Single, a(1 To 367) As SingleConst h1 = 0.33, h2 = 0.14, h3 = 0.33, h4 = 0.2 '各雨量站的权重Const Wm = 140, Um = 20, Lm = 60, Dm = 60Const B = 0.3, C = 0.16, IM = 0.002Private Sub Command1_Click() '任务一运用程序优选KcDim sumQ As Single, sumR As Single, sumR89(2000) As Single, sumQ89(2000) As Single, sumR90(2000) As Single, sumQ90(2000) As SingleDim JD89(2000) As Single, XD89(2000) As Single, j As Integer, JD90(2000) As Single, XD90(2000) As SingleDim JDB(2000) As Single, XDB(2000) As Single, Y As Single, minj As Integer, Min As Single, Kc(2000) As SingleConst Fc = 24For j = 1 To 2000 '运用1989年资料率定Kc(j) = 0.9 + 0.001 * jOpen "C:\Documents and Settings\Administrator\桌面\水文预报\1989年资料.txt" For Input As #1 For i = 1 To 365 '流域平均降雨量计算Input #1, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i)P(i) = 0.33 * P1(i) + 0.14 * P2(i) + 0.33 * P3(i) + 0.2 * P4(i)sumQ89(j) = sumQ89(j) + Q(i) * 24 * 3.6 / 553Next iClose #1W(1) = 110: WU(1) = 10: WL(1) = 40: WD(1) = 60 '流域三层蒸发计算W(1) = WU(1) + WL(1) + WD(1)WMM = Wm * (1 + B)a(1) = WMM * (1 - (1 - (W(1) / Wm)) ^ (1 / (1 + B)))For i = 1 To 365Ep(i) = E0(i) * Kc(j)Next iFor i = 1 To 365If WU(i) + P(i) >= Ep(i) ThenEU(i) = Ep(i)EL(i) = 0ED(i) = 0End IfIf WU(i) + P(i) < Ep(i) ThenIf WL(i) >= C * Lm ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * WL(i) / LmED(i) = 0ElseIf WL(i) < C * Lm And WL(i) >= C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * CED(i) = 0ElseIf WL(i) < C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = WL(i)ED(i) = (Ep(i) - EU(i)) * C - EL(i)End IfEnd IfE(i) = EU(i) + EL(i) + ED(i)PE(i) = P(i) - E(i) '流域产流计算If PE(i) > 0 Then '当产流时If PE(i) + a(i) < WMM ThenR(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)W(i + 1) = W(i) + PE(i) - R(i)a(i + 1) = PE(i) + a(i)ElseIf PE(i) + a(i) >= WMM ThenR(i) = PE(i) + W(i) - WmW(i + 1) = Wma(i + 1) = WMMEnd IfEnd IfIf WU(i) + P(i) - EU(i) - R(i) <= Um ThenWU(i + 1) = WU(i) + P(i) - EU(i) - R(i)WL(i + 1) = WL(i) - EL(i)WD(i + 1) = WD(i) - ED(i)ElseWU(i + 1) = UmIf WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) <= Lm ThenWL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)WD(i + 1) = WD(i) - ED(i)ElseWL(i + 1) = LmIf WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <= Dm Then WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm ElseWD(i + 1) = DmEnd IfEnd IfEnd IfIf PE(i) <= 0 Then '当不产流时R(i) = 0W(i + 1) = W(i) + PE(i)a(i + 1) = WMM * (1 - (1 - W(i + 1) / Wm) ^ (1 / (1 + B)))End IfNext iFor i = 1 To 365sumR89(j) = sumR89(j) + R(i)Next iNext jFor j = 1 To 2000 '运用1990年资料率定Kc(j) = 0.9 + 0.001 * jOpen "C:\Documents and Settings\Administrator\桌面\水文预报\1990年资料.txt" For Input As #2 For i = 1 To 365Input #2, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i)P(i) = 0.33 * P1(i) + 0.14 * P2(i) + 0.33 * P3(i) + 0.2 * P4(i)sumQ90(j) = sumQ90(j) + Q(i) * 24 * 3.6 / 553Next iClose #2W(1) = 110: WU(1) = 10: WL(1) = 40: WD(1) = 60W(1) = WU(1) + WL(1) + WD(1)WMM = Wm * (1 + B)a(1) = WMM * (1 - (1 - (W(1) / Wm)) ^ (1 / (1 + B)))For i = 1 To 365Ep(i) = E0(i) * Kc(j)Next iFor i = 1 To 365If WU(i) + P(i) >= Ep(i) ThenEU(i) = Ep(i):EL(i) = 0:ED(i) = 0End IfIf WU(i) + P(i) < Ep(i) ThenIf WL(i) >= C * Lm ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * WL(i) / LmED(i) = 0ElseIf WL(i) < C * Lm And WL(i) >= C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * CED(i) = 0ElseIf WL(i) < C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = WL(i)ED(i) = (Ep(i) - EU(i)) * C - EL(i)End IfEnd IfE(i) = EU(i) + EL(i) + ED(i)PE(i) = P(i) - E(i)If PE(i) > 0 ThenIf a(i) + PE(i) < WMM ThenR(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)a(i + 1) = PE(i) + a(i)W(i + 1) = W(i) + PE(i) - R(i)ElseR(i) = PE(i) + W(i) - Wma(i + 1) = WMMW(i + 1) = WmEnd IfElseR(i) = 0W(i + 1) = W(i) + PE(i)a(i + 1) = WMM * (1 - (1 - W(i + 1) / Wm) ^ (1 / (1 + B)))End IfIf WU(i) + P(i) - EU(i) - R(i) <= Um ThenWU(i + 1) = WU(i) + P(i) - EU(i) - R(i)WL(i + 1) = WL(i) - EL(i)WD(i + 1) = WD(i) - ED(i)ElseWU(i + 1) = UmIf WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) <= Lm ThenWL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)WD(i + 1) = WD(i) - ED(i)ElseWL(i + 1) = LmIf WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <= Dm Then WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm ElseWD(i + 1) = DmEnd IfEnd IfEnd IfNext iFor i = 1 To 365sumR90(j) = sumR90(j) + R(i)Next iNext jFor j = 1 To 2000JD89(j) = sumR89(j) - sumQ89(j)XD89(j) = (sumR89(j) - sumQ89(j)) / sumQ89(j)JD90(j) = sumR90(j) - sumQ90(j)XD90(j) = (sumR90(j) - sumQ90(j)) / sumQ90(j)Next jFor j = 1 To 2000JDB(j) = Abs(XD90(j)) - Abs(XD89(j))XDB(j) = (Abs(XD90(j)) - Abs(XD89(j))) / Abs(XD89(j))Next jMin = 1For j = 1 To 2000 '运用尽量相近原则优选KcIf Abs(XDB(j)) < Min ThenMin = Abs(XDB(j))minj = jEnd IfNext jLabel2.Caption = Kc(minj) ‘输出计算结果JD89(minj) = sumQ89(minj) - sumR89(minj) '绝对误差XD89(minj) = (sumQ89(minj) - sumR89(minj)) / sumQ89(minj) '相对误差Text1.Text = sumQ89(minj)Text2.Text = sumR89(minj)Text3.Text = JD89(minj)Text4.Text = XD89(minj)JD90(minj) = sumQ90(minj) - sumR90(minj) '绝对误差XD90(minj) = (sumQ90(minj) - sumR90(minj)) / sumQ90(minj) '相对误差Text5.Text = sumQ90(minj)Text6.Text = sumR90(minj)Text7.Text = JD90(minj)Text8.Text = XD90(minj)End SubPrivate Sub Command2_Click()'任务二次洪流量计算Dim Qg(1 To 28) As Single, Qs(1 To 28) As Single, UH(1 To 28) As IntegerConst Fc = 11, Cg = 0.978, Qgchu = 55.3Kc = Val(Label2.Caption)Open "C:\Documents and Settings\Administrator\桌面\水文预报\暴雨资料.txt" For Input As #3 For i = 1 To 28Input #3, E0(i), P1(i), P2(i), P3(i), P4(i)P(i) = h1 * P1(i) + h2 * P2(i) + h3 * P3(i) + h4 * P4(i)Ep(i) = E0(i) * KcNext iClose #3W(1) = 140: WU(1) = 20: WL(1) = 60: WD(1) = 60W(1) = WU(1) + WL(1) + WD(1)WMM = Wm * (1 + B)a(1) = WMM * (1 - (1 - (W(1) / Wm)) ^ (1 / (1 + B)))For i = 1 To 28If WU(i) + P(i) >= Ep(i) ThenEU(i) = Ep(i):EL(i) = 0:ED(i) = 0End IfIf WU(i) + P(i) < Ep(i) ThenIf WL(i) >= C * Lm ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * WL(i) / LmED(i) = 0ElseIf WL(i) < C * Lm And WL(i) >= C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * CED(i) = 0ElseIf WL(i) < C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = WL(i)ED(i) = (Ep(i) - EU(i)) * C - EL(i)End IfEnd IfE(i) = EU(i) + EL(i) + ED(i)PE(i) = P(i) - E(i) '流域产流计算If PE(i) > 0 ThenIf PE(i) + a(i) < WMM ThenR(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)W(i + 1) = W(i) + PE(i) - R(i)a(i + 1) = PE(i) + a(i)ElseIf PE(i) + a(i) >= WMM ThenR(i) = PE(i) + W(i) - WmW(i + 1) = Wma(i + 1) = WMMEnd IfEnd IfIf WU(i) + P(i) - EU(i) - R(i) <= Um ThenWU(i + 1) = WU(i) + P(i) - EU(i) - R(i)WL(i + 1) = WL(i) - EL(i)WD(i + 1) = WD(i) - ED(i)ElseWU(i + 1) = UmIf WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) <= Lm ThenWL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)WD(i + 1) = WD(i) - ED(i)ElseWL(i + 1) = LmIf WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <= Dm Then WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm ElseWD(i + 1) = DmEnd IfEnd IfEnd IfIf PE(i) <= 0 ThenR(i) = 0W(i + 1) = W(i) + PE(i)a(i + 1) = WMM * (1 - (1 - W(i + 1) / Wm) ^ (1 / (1 + B)))End IfNext iFor i = 1 To 28 '水源划分If PE(i) > 0 ThenIf PE(i) <= Fc ThenRG(i) = R(i)RS(i) = 0ElseRG(i) = Fc * R(i) / PE(i)RS(i) = R(i) - RG(i)End IfEnd IfIf PE(i) <= 0 ThenR(i) = 0RG(i) = 0RS(i) = 0End IfNext i'出流系数法推求地下径流Qg(1) = Cg * Qgchu + (1 - Cg) * RG(1) * 553 / (3 * 3.6)For i = 2 To 28Qg(i) = Cg * Qg(i - 1) + (1 - Cg) * RG(i) * 553 / (3 * 3.6)Next i'单位线推求直接径流Open "C:\Documents and Settings\Administrator\桌面\水文预报\单位线.txt" For Input As #4 For i = 1 To 11Input #4, UH(i)Next iClose #4For i = 1 To 28For j = 1 To 28If 1 <= i - j + 1 And i - j + 1 <= 28 ThenQs(i) = Qs(i) + RS(j) / 10 * UH(i - j + 1)End IfNext jNext iFor i = 1 To 28 '总的流量Q(i) = Qs(i) + Qg(i)Next i'次洪计算结果输出Open "C:\Documents and Settings\Administrator\桌面\水文预报\次洪流量过程.txt" For Output As #5For i = 1 To 28Print #5, "直接径流:" & Qs(i); "地下径流:" & Qg(i); "次洪总流量:" & Q(i)Next iClose #5End SubPrivate Sub Command3_Click() ‘退出EndEnd SubPrivate Sub Command4_Click() ‘人工优选和检验Dim R(1 To 366) As Single, RS(1 To 366) As Single, RG(1 To 366) As SingleDim PE(1 To 366) As Single, W(1 To 366) As Single, WU(1 To 366) As Single, WL(1 To 366) As Single, WD(1 To 366) As SingleDim WMM As Single, a(1 To 366) As SingleDim sumQ As Single, sumR As Single, sumQ90 As Single, sumR90 As SingleKc = Val(Text9.Text)If Option1.Value = True ThenOpen "C:\Documents and Settings\Administrator\桌面\水文预报\1989年资料.txt" For Input As #1For i = 1 To 365Input #1, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i)P(i) = h1 * P1(i) + h2 * P2(i) + h3 * P3(i) + h4 * P4(i)sumQ = sumQ + Q(i) * 24 * 3600 * 1000 / 553000000Text1.Text = sumQNext iClose #1End IfIf Option2.Value = True ThenOpen "C:\Documents and Settings\Administrator\桌面\水文预报\1990年资料.txt" For Input As #2For i = 1 To 365Input #2, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i)P(i) = h1 * P1(i) + h2 * P2(i) + h3 * P3(i) + h4 * P4(i)sumQ90 = sumQ90 + Q(i) * 24 * 3600 * 1000 / 553000000Text5.Text = sumQ90Next iClose #2End IfIf Option3.Value = True ThenOpen "C:\Documents and Settings\Administrator\桌面\水文预报\1991年资料.txt" For Input As #6For i = 1 To 365 '运用1991年资料检验Input #6, Q(i), E0(i), P1(i), P2(i), P3(i), P4(i)P(i) = h1 * P1(i) + h2 * P2(i) + h3 * P3(i) + h4 * P4(i)sumQ = sumQ + Q(i) * 24 * 3600 * 1000 / 553000000Text1.Text = sumQNext iClose #6End IfW(1) = 110: WU(1) = 10: WL(1) = 40: WD(1) = 60W(1) = WU(1) + WL(1) + WD(1)WMM = Wm * (1 + B)a(1) = WMM * (1 - (1 - (W(1) / Wm)) ^ (1 / (1 + B)))For i = 1 To 365Ep(i) = E0(i) * KcNext iFor i = 1 To 365If WU(i) + P(i) >= Ep(i) ThenEU(i) = Ep(i): EL(i) = 0: ED(i) = 0End IfIf WU(i) + P(i) < Ep(i) ThenIf WL(i) >= C * Lm ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * WL(i) / LmED(i) = 0ElseIf WL(i) < C * Lm And WL(i) >= C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = (Ep(i) - EU(i)) * CED(i) = 0ElseIf WL(i) < C * (Ep(i) - EU(i)) ThenEU(i) = WU(i) + P(i)EL(i) = WL(i)ED(i) = (Ep(i) - EU(i)) * C - EL(i)End IfEnd IfE(i) = EU(i) + EL(i) + ED(i)PE(i) = P(i) - E(i)If PE(i) > 0 ThenIf PE(i) + a(i) < WMM ThenR(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)W(i + 1) = W(i) + PE(i) - R(i)a(i + 1) = PE(i) + a(i)ElseIf PE(i) + a(i) >= WMM ThenR(i) = PE(i) + W(i) - WmW(i + 1) = Wma(i + 1) = WMMEnd IfEnd IfIf WU(i) + P(i) - EU(i) - R(i) <= Um ThenWU(i + 1) = WU(i) + P(i) - EU(i) - R(i)WL(i + 1) = WL(i) - EL(i)WD(i + 1) = WD(i) - ED(i)ElseWU(i + 1) = UmIf WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) <= Lm ThenWL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)WD(i + 1) = WD(i) - ED(i)ElseWL(i + 1) = LmIf WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <= Dm Then WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm ElseWD(i + 1) = DmEnd IfEnd IfEnd IfIf PE(i) <= 0 ThenR(i) = 0W(i + 1) = W(i) + PE(i)a(i + 1) = WMM * (1 - (1 - W(i + 1) / Wm) ^ (1 / (1 + B)))End IfNext iIf Option1.Value = True Or Option3.Value = True ThenFor i = 1 To 365sumR = sumR + R(i)Next iText2.Text = sumRText3.Text = sumQ - sumRText4.Text = (sumQ - sumR) / sumQEnd IfIf Option2.V alue = True ThenFor i = 1 To 365sumR90 = sumR90 + R(i)Next iText6.Text = sumR90Text7.Text = sumQ90 - sumR90Text8.Text = (sumQ90 - sumR90) / sumQ90End IfEnd Sub第七章：总结和心得水文预报是一项直接服务于国家安全和国民经济建设的不可或缺的重要基本工作，是帮助人类有效地预防洪水、减少洪灾损失，有效利用水资源的非工程措施之一。