水文预报课程设计报告
水文预报课程设计
《水文预报》课程设计报告学院：_____水利与环境学院_____专业：____水文与水资源工程____班级： 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 二．流域产流方案...............................................
水文预报课程设计报告书
《水文预报课程设计》姓名：学号：学院：水利与环境学院班级：指导老师：时间：2013.1.23《水文预报课程设计》说明书1、设计目的1、流域水文模型的用途：洪水预报方案是现代实时洪水预报调度系统的核心部分，是提高预报精度和增长预见期的关键技术。
对水资源可持续利用：流域水文模型是水资源评价、开发、利用和管理的理论基础。
对水环境和生态系统保护：流域水文模型是构建面污染模型和生态评价模型的主要平台。
流域水文模型还是分析研究气候变化和人类活动对洪水、水资源和水环境影响的有效工具。
本次课程设计的目的是通过一个具体的降雨～径流预报方案的制作，使学生了解生产单位对预报任务的要求。
2、通过课程设计，要求掌握如下内容：1)流域综合退水曲线、地下水退水曲线的制作以及次洪分割方法；2)熟悉降雨～径流相关图编制的完整过程；3)新安江两水源模型结构及产流参数率定方法；4)流域经验单位线的推求方法；5)洪水预报方案精度评定方法；6)利用预报方案进行实行洪水预报方法；7)利用马斯京根分段连续演算法进行长河段洪水演进预报。
2、设计基本资料该流域集水面积1884.6km2，干流河长约273km。
流域气候温湿，年降雨量在1700毫左右。
地下水位较高，且随季节变幅小，因次，一般情况下，土壤含水量较大。
根据流域自然地理条件情况和气候条件，以及洪水流量过程线分析，可知流域产流规律符合湿润地区的蓄满产流模型特征。
采用降雨径流相关图制作流域蓄满产流方案，用二层蒸发模型计算蒸发，水源划分考虑两水源划分。
3、课程设计资料1. 一场历史洪水的流量过程，相应的各雨量站（3个雨量站）时段雨量与权重（时段长为3h ）2. 洪水的前期日降雨量，日最大蒸发量资料（历史洪水与实时预报洪水）3. XX 场洪水的退税过程4. 部分场次洪水降雨、径流特征值成果表5. 干流河段的马斯京根参数及分段数（3段）4、产流计算4.1 面平均雨量计算 4.11 计算方法泰森多边形法i i P P α⨯=∑(1)其中i α为i 站的面积权重，i P 为i 站的实测降雨 4.12 各测站基本资料如表1、表2所示表1 各测站降雨量 单位:mmA 站降雨B 站降雨C 站降雨 18.18181818 8.510638298 12.90322581 90.90909091 42.55319149 64.51612903 168.1818182 78.72340426 119.3548387 45.45454545 21.27659574 32.25806452 36.3636363617.021276625.80645161表2 各测站权重A 站权重B 站权重C 站权重 0.220.470.314.13 计算步骤将表1、表2中数据代入公式(1)中并使其结果除以3即得面平均降雨量如表3所示表3 面平均雨量 单位:mm流域面雨量4.00 20.00 37.00 10.00 8.00再由实测洪水过程可画出流量过程线如图1图1 实测洪水流量过程线4.2 次洪分割4.21 计算流域平均消退系数Cgg由最小二乘法公式 ∑∑==∧=n i i ni i i g QQ QC 121,12,1, (2)算得Cg=0.837(从流量开始消退时算起)同理可算得地下水消退系数Cgg=0.746(从退水开始算起)4.22 划分水源由地下水退水过程线不变的规律，结合实测洪水过程(无后续洪水)和实测洪水过程可画出图2图2 分割复式洪水再由洪峰(如此题中的683.91)与Cg 的乘积的递推，算出一系列的Q 值，用径流深计算公式可累加出一系列的R 值，即得蓄泻关系曲线，如图3图3 蓄泻关系曲线用图3可查出此场洪水对应的初和末e e R R 值，则由公式初末e 12210/26.3R R A Q Q Q t R e n i i -+⎪⎭⎫ ⎝⎛++∆=∑-= (2)算得mm R 637.390=再算得mm R S 814.34= mm R G 832.4= 另附以下几张图表表4 地面与地下流量过程 单位:mmQ Qg Qs 1.7571 1.7571 0 1.9518 4.2332 0.8572 17.3571 6.7094 10.6477 117.0003 9.1856 107.8146 361.8473 11.6612 350.1856 620.1595 14.1379 606.0216 683.9124 16.6140 667.2984 666.2792 19.0902 647.1887 558.2343 21.5663 536.6681 430.3991 24.0425 406.3565 315.2504 26.5187 288.7317 235.7673 28.9948 206.7724 168.2515 31.4710 136.7805 112.7568 33.9471 78.8096 62.1171 36.4233 25.6938 38.899438.8994图5 流域平均退水曲线图6 地下水退水曲线表5 蓄泻关系曲线表单位:m3/s Q R0.00068891 0.000002360.00072262 0.000007210.00075798 0.000012290.00079508 0.000017620.00083399 0.000023210.0008748 0.000029070.00091762 0.000035220.00096252 0.000041680.001009630.00004844图7 相邻时段流量关系图4.3 流域初始土壤蓄水量的计算 4.31 前期影响雨量Pa 的计算WMEpk -=1 (3) )(,1,t a t t a P P k P +=+ (4) 4.32 计算实例已知资料见表6表6 Pa 基本数据 单位:mm日降雨量 Pa 蒸发能力 WM Pa,t+1 1.4 60 9.7 80 53.95525 0 47.41318 5.9 46.84895 1.9 42.83814 0.738.259144.4 降雨径流相关图制作及误差检验 4.41 计算原理：（1）流域蓄满产流模型计算公式)()1(1WMM PE a WMMa PE WM WM W PE Rb <=++-+-+=+ (5) )(WMM PE a WMW PE R >+-+=(6)（2）蒸发计算（两层蒸发模式）PL U L P U P U E E E E E E E W P =+===⇒>+E 0(7)LU LMU L U U P U E E W E E W E E W P +===⇒<=+E W )-(E L P (8)表7 降雨径流相关表单位:mm W0 0 10 20 30 38.25914 40 50 60 63.1 70 80 P \ a 0 10.1638 20.695677 31.6747 41.15998 43.21656 55.506258 68.879817 73.3425 84.11941 106.40 0 0 0 0 -7.1E-15 -7.1E-15 0 0 -7.1E-15 7.11E-15 02.7875 0.01212 0.10381 0.2060671 0.32223 0.43239 0.457702 0.6221349 0.835937 0.919062 1.158914 2.787510 0.15846 0.49578 0.873091 1.30347 1.713695 1.808288 2.4264865 3.2432283 3.566412 4.528131 1015.75 0.39826 0.94089 1.549526 2.24635 2.913684 3.068082 4.0830814 5.4464816 5.996761 7.704987 15.7520 0.64862 1.34905 2.1364135 3.04064 3.910063 4.111797 5.444911 7.263802 8.013052 10.48246 2021.575 0.75764 1.51796 2.3733851 3.35697 4.304224 4.524275 5.9815807 7.9834773 8.816008 11.67635 21.57530 1.4958 2.59071 3.8290544 5.26353 6.659053 6.985688 9.1807722 12.35853 13.81252 20 3040 2.73081 4.25936 6.0016437 8.04337 10.06335 10.54252 13.86103 20 23.1 30 4050 4.39211 6.40489 8.7230367 11.4852 14.29346 14.97618 20.138794 30 33.1 40 5052.5 4.87954 7.02318 9.4999044 12.467 15.51396 16.26221 22.5 32.5 35.6 42.5 52.557.3625 5.91606 8.32778 11.134169 14.5397 18.13038 19.04014 27.3625 37.3625 40.4625 47.3625 57.362560 6.52935 9.09493 12.093885 15.7648 19.71786 20.75181 30 40 43.1 50 6066.15 8.10818 11.0603 14.556409 18.9586 24.40914 26.15 36.15 46.15 49.25 56.15 66.1570 9.20978 12.4279 16.280051 21.2713 28.25914 30 40 50 53.1 60 7072.9375 10.1134 13.5498 17.705627 23.2865 31.19664 32.9375 42.9375 52.9375 56.0375 62.9375 72.937581.25 12.9982 17.151 22.425269 31.25 39.50914 41.25 51.25 61.25 64.35 71.25 81.2580 12.5311 16.5645 21.633086 30 38.25914 40 50 60 63.1 70 8090 16.6527 21.8387 30 40 48.25914 50 60 70 73.1 80 90100 21.9032 30 40 50 58.25914 60 70 80 83.1 90 100 110.000 30.000 40.000 50.000 60.000 68.259 70.000 80.000 90.000 93.100 100.000 110.000 120.000 40.000 50.000 60.000 70.000 78.259 80.000 90.000 100.000 103.100 110.000 120 130.000 50.000 60.000 70.000 80.000 88.259 90.000 100.000 110.000 113.100 120.000 130.000 140.000 60.000 70.000 80.000 90.000 98.259 100.000 110.000 120.000 123.100 130.000 140.000.专业.整理..专业.整理.图8 降雨径流相关图4.4.2 误差检验用其他场次PE～R关系查降雨径流相关图，校核是否合格并算出合格率与等级，见表8表8 径流深校核单位：mmPa 有效降雨PE 次洪径流深演算径流深许可误差实际误差合格否合格率等级15.2 45.7 10.36973824 6.452868934 3 3.916869311 否0.88 甲级48.8 58.2 29.95190859 27 5.990381719 2.951908594 是43.9 25.2 2.406251416 6.160410485 3 -3.75415907 否46 81.5 46.96073863 47.5 9.392147726 -0.53926137 是40.1 113.3 71.16978066 71.55914295 14.23395613 -0.3893623 是45.3 66 31.99972577 31.3 6.399945153 0.699725765 是35.7 64.7 27.56503271 20.70111833 5.513006542 6.863914379 否10 71.9 10.51539884 13.14539778 3 -2.62999893 是29.8 76.2 31.84636564 26 6.369273129 5.846365644 是44.4 114.2 77.20418345 78.6 15.44083669 -1.39581655 是44.9 81.9 44.8442484 46.8 8.96884968 -1.9557516 是51 65.9 38.72889177 36.9 7.745778353 1.828891765 是80 101.7 103.827372 101.7 20 2.127372048 是10.7 37.6 4.762956997 3.923568862 3 0.839388135 是52.4 100.5 68.49763184 72.9 13.69952637 -4.40236816 是42.7 74.7 35.32604807 37.4 7.065209614 -2.07395193 是39.6 150 106.9960346 109.6 20 -2.60396537 是27.8 60.5 13.56898305 15.10650626 2.71379661 -1.53752321 是38 94.9 50.55374053 52.9 10.11074811 -2.34625947 是52.4 85.1 58.77736748 57.5 11.7554735 1.277367482 是.专业.整理.47.9 46.7 19.71769713 16.68307565 3.943539427 3.034621487 是11.5 109.4 40.1479979 40.9 8.029599579 -0.7520021 是56.4 10 4.5795851 2.916876906 3 1.662708194 是44.7 18.9 6.657111634 4.380469246 3 2.276642387 是.专业.整理.5、汇流计算5.1、fc 的试算 5.1.1试算方法iiFCPE CPE R f RG i ∑>=(9) ()iiFCPE C i PE R f PE RS i ∑<=-=(10) 5.1.2 试算实例 (1) 基本数据如表9表9 单位:mm再假设hmm f C /79.2<，则有h mm f C /79.2959.1)17145.03863.02061.01551.0/(823.4<=++++=故h mm f C /959.1= 划分时段净雨如表10所示表10 划分地面及地下净流深 单位:mmRg Rs 0.30 0.13 0.40 3.47 0.76 13.07 1.40 4.88 1.964.835.2 计算单位线(试错法) 5.2.1 演算实例如表11所示表11 科林法试算表单位:m3/s再用调整的单位线计算地面流量，比较预测与实际值间的误差大小。
三峡大学水文预报课程设计
水文预报课程设计说明书姓名：班级：2013学号：2013学院：水利与环境学院指导老师：2016年1月目录第1章设计任务简介 (1)1.1设计目的.............................................................错误！未定义书签。
1.2设计基本资料 (1)1.3课程设计资料 (1)第2章流域产流方案的编制 (2)2.1预报方案制作 (2)2.1.1流域平均退水曲线 (2)2.1.2径流分割与次洪径流深计算 (4)2.1.3制作降雨径流关系图 (5)2.1.4净雨预报检验 (7)2.1.5二层蒸发模式计算 (8)2.1.6稳定下渗率计算 (8)2.1.7地面径流单位线计算 (9)2.1.8汇流计算 (11)2.1.9次洪模拟预报精度评价 (13)2.2实时洪水预报 (14)2.2.1划分净雨 (14)2.2.2汇流计算 (14)2.2.3河道洪水演算 (16)第3章设计体会 (18)3.1重难点归纳与解决方法 (18)3.2体会 (18)第1章设计任务简介1.1 设计目的1、流域水文模型的用途：洪水预报方案是现代实时洪水预报调度系统的核心部分，是提高预报精度和增长预见期的关键技术。
对水资源可持续利用：流域水文模型是水资源评价、开发、利用和管理的理论基础。
对水环境和生态系统保护：流域水文模型是构建面污染模型和生态评价模型的主要平台。
流域水文模型还是分析研究气候变化和人类活动对洪水、水资源和水环境影响的有效工具。
本次课程设计的目的是通过一个具体的降雨～径流预报方案的制作，使学生了解生产单位对预报任务的要求。
2、通过课程设计，要求掌握如下内容：（1）流域综合退水曲线、地下水退水曲线的制作以及次洪分割方法；（2）掌握次洪径流深及地面、地下流量分割方法；（3）掌握降雨～径流相关图（APN模型）编制的完整过程；（4）新安江两水源模型结构及产流参数率定方法；（5）流域经验单位线的推求方法；（6）洪水预报方案精度评定方法；（7）利用预报方案进行实时洪水预报方法；（8）利用马斯京根分段连续演算法进行长河段洪水演进预报。
水文预报课程设计报告
水文预报课程设计报告学院：资源与环境学院专业：水文与水资源工程年级:2012级姓名：刘燨元学号：2012215334指导教师：余倩2015 年6 月9 日水文预报（hydrologic forecasting）就是据已知的信息对未来一定时期内的水温状态作出定性或定量的预测。
已知信息，广义上指对预报水文状态有影响的一切信息，最常用的是水文与气象要素信息，如降水、蒸发、流量、水位、气温和含沙量等观测信息。
预报的水文状态变量可以是任一水文要素也可以是水文特征量，不同的状态量预报要求的已知信息不同、预报方法不同、预见期也不同。
目前通常预报的水文要素有流量、水位、冰情和旱情等.水文预报方法以水文基本规律、水文模型研究为基础，结合生产实际问题的需要，构成具体的预报方法或预报方案，服务于生产实际.一般水文预报研究的重点和关键有两部分：①共性规律研究，即具有一定普遍性的水文基本规律模拟方法和流域水文模型研究；②个性问题研究，对反映具体问题的特性、方法进行了解，构成具有解决各种具体实际问题的、具有较高预报精度的预报方案。
1、新安江水文模型简介 (4)1。
1 概述 (4)1。
2 新安江模型的基本原理 (4)1.3 新安江模型的结构 (5)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 概述流域水文模型可分为物理模型、概念性模型和系统模型。
在水文预报中，概念性模型和系统模型应用较多,此处主要介绍概念性流域水文模型.概念性流域水文模型属于数学模型，它与物理模型相比，具有许多优点：一是它的所有条件均可由原型观测资料直接给出，不受比尺的限制，即数学模型无相似律问题；二是它的边界条件及其它条件可严格控制，也可随时按实际需要改变；三是它的通用型较强，只要研制出一种适用的应用软件，就可用来解决不同的实际问题；四是它具有理想的抗干扰性能，只要条件不变，重复模拟可以得到相同的结果，不会因人、因地而异；五是它的研制费用相对较低.因此，流域水文模型的研制和应用受到水文学家和水文工作者的普普遍重视。
栾川站水文预报课程设计
栾川站水文预报课程设计一、课程目标知识目标：1. 让学生掌握栾川站所在流域的基本水文特征，包括流域面积、气候类型、降水分布等；2. 使学生了解水文预报的基本原理，如降水-径流过程、水位-流量关系等；3. 帮助学生掌握栾川站历史水文数据的基本分析方法，并运用数据进行短期和长期水文预报。
技能目标：1. 培养学生运用地理信息系统（GIS）和相关软件进行水文数据分析和处理的能力；2. 提高学生根据实际天气和流域情况，独立进行水文预报的技能；3. 培养学生撰写水文预报报告和口头表达能力。
情感态度价值观目标：1. 培养学生对水文预报工作的兴趣和热情，增强环保意识和水资源保护观念；2. 培养学生严谨的科学态度和团队协作精神，养成分析问题、解决问题的良好习惯；3. 增强学生对我国水文事业的认识，激发学生为国家和地方水文预报工作贡献力量的决心。
课程性质：本课程为实践性较强的学科课程，旨在通过实际案例分析，使学生掌握水文预报的基本知识和技能。
学生特点：学生为高中年级，具备一定的地理、数学和科学知识基础，具有较强的学习能力和好奇心。
教学要求：结合学生特点，课程设计应注重理论与实践相结合，充分调动学生的积极性和参与度，提高学生的实际操作能力。
通过课程学习，使学生在知识、技能和情感态度价值观方面均取得具体、可衡量的学习成果。
二、教学内容1. 流域水文特征分析：介绍栾川站所在流域的自然地理概况，包括流域面积、地形地貌、气候类型、降水分布等，关联教材中流域水文特征相关章节。
教学大纲：引导学生通过实地考察和查阅资料，分析流域水文特征对水文预报的影响。
2. 水文预报原理：讲解降水-径流过程、水位-流量关系等基本原理，结合教材中相关理论章节。
教学大纲：通过案例分析，使学生理解并掌握水文预报的基本原理。
3. 水文数据分析：教授学生如何收集、整理和分析栾川站历史水文数据，关联教材中水文数据处理相关章节。
教学大纲：指导学生运用GIS和相关软件进行水文数据分析，培养学生实际操作能力。
水文预报课程设计VB报告
水文预报课程设计VB报告一、课程目标水文预报课程设计VB报告旨在达成以下课程目标：1. 知识目标：（1）使学生掌握水文预报的基本概念、原理和方法；（2）培养学生运用VB编程语言解决水文预报问题的能力；（3）让学生了解我国水文预报的现状和发展趋势。
2. 技能目标：（1）培养学生收集、整理和分析水文数据的能力；（2）提高学生运用VB编程语言编写水文预报程序的能力；（3）锻炼学生通过实际案例，解决实际水文预报问题的能力。
3. 情感态度价值观目标：（1）激发学生对水文预报学科的兴趣，培养其科学探索精神；（2）引导学生关注水文预报在国民经济和社会发展中的作用，增强社会责任感；（3）培养学生团结协作、积极进取的精神风貌。
本课程针对高中年级学生，充分考虑了学生的认知特点、学习兴趣和实际需求。
在教学过程中，注重理论与实践相结合，以实际案例为引导，让学生在掌握基本知识的同时，提高解决实际问题的能力。
通过本课程的学习，使学生能够独立完成水文预报相关任务，为未来进一步学习和从事相关工作奠定坚实基础。
课程目标的设定，旨在为教学设计和评估提供明确的方向，确保教学效果的达成。
二、教学内容水文预报课程设计VB报告的教学内容如下：1. 水文预报基本概念与原理- 水文预报的定义与分类- 水文预报的基本原理- 水文预报的主要方法2. VB编程基础- VB编程环境与基本操作- 数据类型、变量与常量- 控制结构（顺序、选择、循环）- 过程与函数3. 水文数据处理与分析- 水文数据收集与整理- 水文数据分析方法- VB编程实现水文数据分析4. 水文预报模型- 水文预报模型概述- 常见水文预报模型介绍- VB编程实现水文预报模型5. 案例分析与实际应用- 案例分析：某地区洪水预报- VB编程解决实际水文预报问题- 水文预报结果评价与分析教学内容依据课程目标，结合教材章节进行组织。
在教学过程中，按照以下进度安排教学内容：1. 第1-2周：水文预报基本概念与原理2. 第3-4周：VB编程基础3. 第5-6周：水文数据处理与分析4. 第7-8周：水文预报模型5. 第9-10周：案例分析与实际应用三、教学方法水文预报课程设计VB报告的教学方法如下：1. 讲授法：教师通过生动的语言、形象的表达，为学生讲解水文预报的基本概念、原理和VB编程基础等理论知识。
水文预报课程设计
《水文预报》课程设计学号：2013100540姓名：吴旺军指导教师：王永文专业班级：13水文与水资源工程01班南昌工程学院水文与水资源工程2016年05月目录一、背景资料二、课程设计目的和任务三、预报一(降雨径流法)三、预报二(相应水位（流量）法)四、预报三(合成流量法)五、总结一、背景资料水文预报是指根据前期或现时的水文气象资料，对某一水体、某一地区或某一水文站在未来一定时间内的水文情况作出定性或定量的预测。
对防洪、抗旱、水资源合理利用和国防事业中有重要意义。
对未来水文情况作出科学预测并发布预报的技术与作业。
水文预报是水文学为经济和社会服务的重要方面，特别是对灾害性水文现象做出预报，对综合利用大型水利枢纽做出短期、中期和长期的预报，作用很大。
中国已开展预报服务的项目有:洪水水位与流量、枯水水位与流量、含沙量、各种冰情、水质等。
二、课程设计目的和任务熟悉水文预报的基本方法及原理，在实践中提高我们的动手能力，理论结合实际，充分的理解水文预报的精髓并更好的掌握这门知识。
三、预报一：降雨径流预报（1）、预报原理按降雨径流行成过程的原理，利用流域内的降雨资料通过单位线预报出流域出口断面的洪水过程线，主要的就是求P、Pa、R还有实际流量过程线和出口断面的流量过程线。
（2）、预报的过程基本资料：以下的图表的右侧是老师给我们资料所在区域的基本情况，左边为给的七个站的站名与站码、单位线等；除此之外，还给出了从1981到1995年的流量、每个站码的时段雨量与日雨量等资料：图一：站点分布（3）、求总雨量P步骤一：选取一场合适的洪水，并将其绘制在坐标纸上，其总的流量绘制出如下图所示：选取场次洪水的依据为1所选的场次洪水所包含的时间所对应的时段雨量资料齐全，或者资料缺乏的站不多于两个站。
2同时也拥有日雨量资料且场次洪水洪量足够大，特别要注意的是选一场洪水，不能选两场洪水的中间一部分或者是一场洪水的一部分。
以下就是我选的洪水资料。
水文预报课设报告
当 W0=0 时，WM=P-E-R。因此，选用几次 W0 很小，P 很大， 估计雨止时土壤含水量达到最大值的实测洪水资料，及可分析出 WM。 在南方湿润地区，若用三层蒸发模式计算蒸发，WM 可采用 140mm (WUM=20mm, WLM=60mm, WDM=60mm)左右。 本次课程 设计选用 WM=140mm。 b、b 的确定 b 是个地形地质因素，表示流域蓄水量分布的不均匀性。当全 流域蓄满后，这个因素不再起作用，b 可取 0.3~0.4 左右。本次课 程设计选用 b=0.37。 c、蒸发能力 Ep 的推求 流域内若有蒸发皿资料，则可以直接用蒸发皿实测蒸发值乘
(3) 计算 W0 根据水量平衡原理，W0 的计算公式为：
W0，t+1 = W0,t + Pt − Et − R t (3 − 7)
Rt 可由产流计算公式计算；Et 采用三层蒸发计算模型计算；W0 的 起算值，可以以久旱以后（W0=0）或大雨后（W0=WM）开始。
表 1: 密塞流域降水径流计算表（1986 年 3 月至 10 月）单位：(mm)(WM=140 K=0.86 b=0.37)
三、课程设计的原理和内容 .............................. 3
(一) 流域产流方式的初步论证 ................................ 3 (二) 径流预报方案.......................................... 3 1、用分析法计算径流量 .................................. 3 2、由实测资料计算次洪径流量 R 实 ........................ 10 3、成果分析........................................... 12 （三）直接径流汇流过程预报方案 ...........................ห้องสมุดไป่ตู้ 13 1、选用资料........................................... 13 2、从实测流量过程线上分割地下径流 ..................... 14 3、确定单位线时段长∆t ................................. 18 4、从实测降雨过程求净雨过程 ........................... 18 5、单位线分析......................................... 27 6、分析由两次洪水所得出的单位线，确定单位线预报方案 ... 29 （四）地下径流汇流过程预报方案。 .......................... 29 （五）降雨径流过程预报方案的校核预报 ...................... 30 四、设计的体会与建议.......................................... 34
周传争水文预报课程设计报告剖析
水文预报课程设计（文天学院11级）班级学号姓名周成绩河海大院2014年1月4日一基本任务任务一：根据已给数据资料（本文选取1995、1996年的历史数据）、参数，编写程序将流域作为整体：（1）进行日模型产流量计算；（2）将计算年径流与实测年径流进行比较，优选蒸发折算系数Kc；任务二根据已给设计暴雨资料、参数,编写程序，将流域作为整体：(1)进行次洪产流量、划分水源(2)直接径流汇流、地下径流汇流计算；（3）绘出直接径流过程、地下径流过程、总的流量过程(绘制成图)。
二基本资料（1）研究流域概况呈村流域控制面积为290 km2，地势南高北低，相对高差较大，平均海拔高程为583 m，流域河道平均坡度为0.95%，最大汇流路径长度为36 km。
流域内植被良好，雨量充沛，多年平均降雨量约为2100 mm，流域降水在年内年际分配极不均匀，为典型的湿润流域。
该流域植被类型主要包括常绿针叶林、落叶阔叶林、混合林、森林地、林地草原、牧草地与作物地，土壤类型主要为黏壤土。
流域内地势平缓，土壤主要有黄壤和砂壤，具有明显的腐殖层，淀积层和母质土等层次结构，透水性好。
台地、丘陵多生长松、杉、樟等高大乔木；平原则以种植农作物和经济作物为主，植被良好。
三基本数据km。
(1) 计算的呈村流域面积2902(2) 流域内有10个雨量站：呈村、董坑呜、棣甸、樟源口、汪村、用功城、左龙、冯村、田里、大连,各站面积比例分别为：0.1, 0.1, 0.10, 0.1，0.1,0.1, 0.1, 0.1, 0.1, 0.1(3) 资料放在/周传争Data/文件夹下，以年份命名的文件，数据格式为：T（i）、Q(i)、E(i)、P1(i)、 P2(i)、P3(i)、P4(i)、P5(i)、P6(i)、P7(i)、P8(i)、P9(i)、P10(i)。
(4)暴雨资料放在/周传争Data/文件夹下pdata.csv(5)单位线数据在/周传争Data/文件夹下unitgraph.csv(5)计算参数如下表表一任务一计算参数表二任务二参数表表三各参数物理意义四计算公式1.蒸散发计算根据流域特点，蒸散发计算采用的是三层蒸散发计算模式。
水文预课程设计报告指示书
.水文预报课程设计指示书水文预报课程设计指示书题目:制作江西省上犹江水库干流入库站的降雨～径流预报方案1方案采用的技术途径1.1蒸散发方案采用一层蒸散发模式。
1.2产流方案根据流域的自然地理情况和气候条件,以及洪水流量过程线的分析,可知流域产流规律符合湿润地区的蓄满产流特征。
采用蓄满产流的降雨径流相关图形式表达方案。
用一层蒸发模型计算蒸发。
计算时段为3小时。
蓄满产流方案可由流域蓄水容量曲线表达，曲线共有两个参数：[1]流域蓄水容量WM；[2]流域蓄水容量曲线指数b。
根据“水文学原理”或“水文预报”知识，流域蓄水容量曲线的参数确定后，可将流域蓄水容量曲线转化为蓄满产流的降雨径流相关图。
因此，蓄满产流的产流方案也可用蓄满产流的降雨径流相关图表达。
1.3水源划分采用变动稳定下渗率Fc～R关系作两水源划分。
1.4汇流方案⑴分型经验单位线作直接径流汇流方案。
⑵采用矩形入流的马斯京根线性水库演算作地下径流汇流方案。
1.5预报模式预报模式见图1。
图1 预报模式示意图2流域概况上犹江发源于湖南省汝城县诸广山的东南麓，干流称为古亭水。
上犹江水库位于江西省上犹县，水库建于古亭水之上。
水库入库站——麟潭站控制面积1067km2，干流河长93km。
流域地貌属低山丘陵区，以低山分布为主，相对高差多在500m左右。
上游部分地区分布有原始森林, 森林植被以松树、杉树、竹类为主，灌木次之；山间盆地种植农作物，流域植被率在80%以上。
土壤多为红色砂壤土,间有亚粘土层。
山坡残积坡积一般厚度为1～2m,最厚者约4～5m。
在山麓坡积层与基岩接触面上，或河流下切至接触部位时，常有泉水出露、涌水量较大者可达每秒数升左右。
流域气候温暖湿润，年降雨量为1700mm左右。
汛期4～9月降雨量约占全年降雨量的73%左右；冬季有降雪，但量不大。
地下水位一般较高,且季节性变幅较小,因此,一般情况下，土壤含水量较大。
洪水流量过程线极不对称，涨洪历时仅数小时至十多小时,而洪水退水历时可达数日至十余日。
- 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第七章：总结和心得水文预报是一项直接服务于国家安全和国民经济建设的不可或缺的重要基本工作，是帮助人类有效地预防洪水、减少洪灾损失，有效利用水资源的非工程措施之一。