水文预报课设报告
桥涵水文课程设计
桥涵水文课程设计一、教学目标本课程的教学目标是让学生掌握桥涵水文的基本知识、技能和应用方法。
通过本课程的学习，学生应能理解桥涵水文的基本概念、原理和方法，掌握水文测验、水文计算和水文预报的基本技能，并能够运用所学知识解决实际工程问题。
具体来说，知识目标包括：1.掌握桥涵水文的基本概念、原理和方法。
2.了解水文测验、水文计算和水文预报的基本原理和方法。
3.熟悉桥涵工程中水文问题的处理方法和技巧。
技能目标包括：1.能够进行水文测验，包括降雨量、径流量、流速、水位等的观测和计算。
2.能够进行水文计算，包括洪水频率计算、径流系数计算、设计洪水计算等。
3.能够进行水文预报，包括短期洪水预报、长期径流预报等。
情感态度价值观目标包括：1.培养学生的科学精神和创新意识，使其能够积极探求新的知识和方法。
2.培养学生的团队合作意识和能力，使其能够与他人合作完成工程项目。
3.培养学生的社会责任感和职业道德，使其能够考虑到工程对环境和社会的影响，并遵循相关伦理和法律规范。
二、教学内容本课程的教学内容主要包括桥涵水文的基本知识、水文测验、水文计算和水文预报等方面的内容。
具体来说，教学内容的安排如下：1.桥涵水文的基本知识：包括桥涵水文的概念、原理和方法，水文循环的基本过程等。
2.水文测验：包括降雨量、径流量、流速、水位等的观测和计算方法，以及相关仪器的使用和维护。
3.水文计算：包括洪水频率计算、径流系数计算、设计洪水计算等方法，以及相关公式的应用和计算技巧。
4.水文预报：包括短期洪水预报、长期径流预报等方法，以及预报模型的建立和应用。
教材的章节安排如下：1.第一章：桥涵水文的基本知识2.第二章：水文测验3.第三章：水文计算4.第四章：水文预报三、教学方法为了激发学生的学习兴趣和主动性，本课程将采用多种教学方法，包括讲授法、讨论法、案例分析法和实验法等。
讲授法将用于传授桥涵水文的基本知识和原理，通过清晰的讲解和示例，帮助学生理解和掌握相关概念和方法。
水文预报课程设计
《水文预报》课程设计报告学院：_____水利与环境学院_____专业：____水文与水资源工程____班级： 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. 利用统计方法和机器学习算法，对水文数据进行处理和分析；4. 提出改进水文预报方法的建议，提高预报准确性。
三、研究内容本研究将围绕以下几个方面展开：1. 水文要素的时空变化规律分析通过对历史水文数据的收集和整理，分析水文要素（如降雨量、径流量、水位等）的时空变化规律。
利用统计方法和时序分析技术，揭示水文要素的周期性和趋势性变化。
2. 水文要素与环境因素的关联性研究探索水文要素与气象、地质、人类活动等因素之间的关联性。
通过分析多个变量之间的相关性，建立水文要素与环境因素的数学模型，为水文预报提供更准确的输入数据。
3. 水文数据处理和分析针对水文数据中的噪声、缺失值等问题，采用数据插值和异常值处理方法，提高数据质量。
利用统计方法、机器学习算法等技术，对水文数据进行挖掘和分析，提取有用的信息。
4. 改进水文预报方法根据前期研究的结果，提出改进水文预报方法的建议。
可以尝试引入新的模型和算法，结合多个模型进行集成预报，提高预报准确性和可靠性。
四、研究方法本研究将采用以下方法进行实施：1. 数据收集和整理收集历史水文数据、气象数据、地质数据等相关数据，进行数据清洗和整理，建立完整的数据集。
2. 数据分析和建模利用统计方法和时序分析技术，对水文要素的时空变化规律进行分析。
通过回归分析、相关性分析等方法，研究水文要素与环境因素的关联性。
matlab水文预报课程设计
matlab水文预报课程设计一、课程目标知识目标：1. 学生能理解水文预报的基本概念，掌握利用Matlab进行水文预报的基本原理；2. 学生能掌握Matlab编程基础，运用相关函数和工具箱进行水文数据的处理和分析；3. 学生能掌握建立和优化水文预报模型的方法，并运用模型进行洪水预测。
技能目标：1. 学生能运用Matlab软件收集和整理水文数据，进行数据预处理；2. 学生能编写简单的Matlab程序，实现对水文数据的可视化分析；3. 学生能运用Matlab建立水文预报模型，对实际案例进行模拟和预测。
情感态度价值观目标：1. 学生培养对水文预报的兴趣，增强环保意识和水资源管理观念；2. 学生通过团队协作解决问题，培养合作精神和沟通能力；3. 学生在解决实际问题的过程中，培养勇于探索、积极创新的精神。
课程性质：本课程为应用性实践课程，以Matlab软件为工具，结合水文预报理论，培养学生实际操作和解决问题的能力。
学生特点：学生已具备一定的水文预报知识和Matlab基础，具有较强的学习兴趣和动手能力。
教学要求：教师需引导学生将理论知识与实际应用相结合，注重培养学生的实际操作能力和创新思维，提高学生的水文预报技能。
在教学过程中，关注学生的学习进度和反馈，及时调整教学方法和内容。
通过课程学习，使学生能够独立完成水文预报模型的建立和优化，为实际工作打下坚实基础。
二、教学内容1. 水文预报基本概念：介绍水文预报的定义、意义、方法及在水资源管理中的应用。
教材章节：第一章 水文预报概述2. Matlab软件入门：讲解Matlab软件的基本操作、数据类型、矩阵运算和编程基础。
教材章节：第二章 Matlab基础3. 水文数据处理：学习利用Matlab进行水文数据的收集、预处理、可视化分析等操作。
教材章节：第三章 水文数据处理4. 水文预报模型：介绍常见的水文预报模型及其原理，如ARIMA模型、SVM 模型等。
教材章节：第四章 水文预报模型5. Matlab在水文预报中的应用：讲解如何利用Matlab软件建立和优化水文预报模型。
水文预报课程设计洪水预报
水文预报课程设计洪水预报一、课程目标知识目标：1. 学生能理解水文预报的基本概念，掌握洪水形成的原因及其发展过程。
2. 学生能够掌握洪水预报的主要方法及其适用条件，如降雨径流模型、统计模型等。
3. 学生能够了解我国洪水预报的现状及发展趋势，了解相关法规政策。
技能目标：1. 学生能够运用所学知识，分析特定流域的洪水形成过程，进行简单的洪水预报。
2. 学生能够运用水文预报软件，进行数据收集、处理和分析，提高解决实际问题的能力。
3. 学生能够通过小组合作，进行洪水预报案例的研究，提高沟通协作能力。
情感态度价值观目标：1. 学生能够认识到水文预报在防洪减灾中的重要作用，增强社会责任感和使命感。
2. 学生在学习过程中，培养严谨的科学态度，树立正确的价值观。
3. 学生通过了解我国在水文预报方面的成就，增强民族自豪感，激发为国家和人民服务的情怀。
本课程针对高中年级学生，结合学科特点，注重理论知识与实践操作的相结合，培养学生的实际应用能力。
课程设计以学生为中心，充分考虑学生的认知水平、兴趣和需求，采用案例教学、小组合作等方法，激发学生的学习兴趣，提高教学效果。
通过本课程的学习，使学生能够掌握洪水预报的基本知识和技能，提高防洪减灾意识，为我国水文预报事业贡献力量。
二、教学内容1. 水文预报基本概念：洪水定义、洪水分类、洪水周期与洪水频率。
2. 洪水形成原因及发展过程：降水过程、流域特性、径流形成与汇集。
3. 洪水预报方法：- 降雨径流模型：水箱模型、单位线法、流域水文模型。
- 统计模型：时间序列分析、回归分析、人工神经网络。
4. 洪水预报软件应用：数据收集、处理、分析和预报结果输出。
5. 我国洪水预报现状与发展趋势：技术进展、政策法规、防洪减灾体系。
6. 实践案例：分析特定流域洪水预报实例，掌握预报流程和操作方法。
教学内容依据课程目标，参照教材相关章节进行组织。
教学大纲安排如下：第一周：水文预报基本概念及洪水形成原因。
《工程水文学》(第4版)第5章水文预报
主 讲：张 峰 手 机：13856993385 邮 箱：zf_hhu@
安徽水电学院水文与水资源专业教研室 二Ο 一五年二月一日
工程水文学 (Engineering Hydrology)
第5章 水文预报
二Ο一五年二月一日
一、概述
1、降雨径流预报
几种具体确定预见期的方法： 1）对于源头流域可把主要降雨结束到 预报断面洪峰出现这个时间差作为洪水预见 期。 2）对于区间流域洪水预报或河段洪水 预报，当区间来水对预报断面洪峰影响不大 时，洪水预见期就等于上下游断面间水流的 传播时间；如果暴雨中心集中在区间流域， 那么预见期就接近于区间洪水主要降雨结束 到下游预报断面洪峰出现这个时差。假如降 雨空间分布比较均匀，上断面和区间都形成 了有影响的洪水，其预见期常取河段传播时 间和区间流域水流平均汇集时间的最小值。
水利部水文局水文情报预报中心
二、短期洪水预报
短期洪水预报包括降雨径流预报、河段洪水预报以及考虑实时修正的实时洪水 预报。
降雨径流预报是指按降雨径流形成过程的原理，利用流域内的降雨资料预报出 流域出口断面的洪水过程。
河段洪水预报是指以河槽洪水波运动理论为基础，由河道上游断面的水位、流 量过程预报下游断面的水位和流量过程。
2、河段洪水预报
（2）合成流量法 在有支流河段，若支流来水量大，干支流洪水之间干扰影响不可忽略，可采
用合成流量法。所谓合成流量是指同时到达下游断面的各上游站相应流量之和。 以下式为根据建立预报方案称为合成流量法。
合成流量法的关键是Ʈ值的确定。
n
Qi f( I ) i,t i i 13、时洪水预报1、水文预报
水文预报是指根据已知信息对未来一定时期内的水文情势作出定性或定量预 报。
水文预报课程设计教学体会与问题探索
水文预报课程设计教学体会与问题探索水文预报是水文学的一个重要分支，它主要是通过气象、降雨模型等手段对未来水文情况进行预测和预报，对水资源管理、防洪减灾等方面具有重要的应用价值。
水文预报课程是水文学专业的重要课程之一，通过学习这门课程，学生能够掌握水文预报的基本理论和方法，提高水文预报能力，为将来从事水资源管理、水灾防治等工作奠定扎实的理论基础。
在教学实践中，我对水文预报课程设计和教学体会进行了总结与问题探索，希望能够对今后的教学工作有所借鉴和提升。
一、教学体会在水文预报课程的教学中，我深刻体会到了以下几点：1. 理论和实践相结合水文预报课程中涉及到很多理论知识，如降雨模型、径流模型、水文预报方法等。
在教学过程中，我注重理论与实践相结合，引导学生通过实际案例分析和水文数据模拟，加深对理论知识的理解和掌握。
通过实践操作，学生们能够更加深入地理解水文预报理论知识，提高解决实际问题的能力。
2. 知识串联与综合应用在教学中，我将水文预报课程与气象学、地质学、水资源管理等相关学科进行了深入的串联，引导学生全面、综合地理解水文预报的理论和方法。
通过案例分析和课程设计，让学生感受到水文预报在实际工程中的应用，培养学生的综合应用能力。
3. 实践能力培养水文预报课程注重培养学生的实践能力，我通过课程设计和实验操作，引导学生进行水文数据的分析和预报模拟，让他们亲自动手，掌握水文预报的基本方法和技能。
我还组织实地考察和调研活动，让学生亲自走出教室，了解真实的水文工程和案例，培养学生的实践操作和创新能力。
二、问题探索在教学实践中，我也遇到了一些问题和困惑，需要进一步探索和改进：1. 教学手段更新随着科技的发展，教学手段和教学资料也需要与时俱进。
在水文预报课程的教学中，我发现学生对新型的数据处理软件、水文分析工具等存在一定的陌生和不熟悉。
我需要进一步学习和掌握新的教学手段和软件工具，引入和应用于水文预报课程的教学中，提高教学效果。
水文与水资源水文预报课设报告
《水文预报》课程综合实训任务书一、实训任务本次实训的任务主要包括以下两个方面：（一）编制伊河栾川以上流域降雨径流量预报方案（二）编制伊河栾川以上流域降雨径流过程预报方案二、实训目的水文预报课程是水文水资源专业的核心技能课程，通过水文预报综合实训，能进一步掌握水文预报的基本理论和基本方法，培养依据设计任务和所给资料，运用所学知识和规范，编制水文预报方案的专业技能，同时提高技术报告的编写能力和计算机的应用能力。
在业务组织能力和实际工作能力方面得到全面锻炼，为毕业后直接顶岗参与水文预报工作奠定基础。
三、流域概况（一）流域自然地理概况：伊河栾川水文站建于东径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所示，其中，桓仁水电站作为龙头电站，具有不完全年调节能力，回龙和太平哨为日调节电站，桓仁水电站运行决定着整个库群的发电效益。
桓仁水电站位于浑江中游。
浑江发源于长白山系龙岗山脉的老爷岭南麓，是鸭绿江右侧的最大支流，全长432km，自东北向西南流经吉林、辽宁两省。
桓仁水电站坝址以上河流长247km，控制流域面积为10400km2，流域内共有10个雨量站，分别为三岔子、八道江、八里哨、通化、三棵树、八道沟、东村、业主沟、花甸、桓仁。
浑江流域地处我国东北暴雨中心北部边缘，多山地，山势陡峭，植被良好，属温带季风型大陆性气候，多年平均年降水量860mm，多年平均径流系数为0.52。
70%的雨量集中在6～9月间，大洪水发生在7月下旬～8月中旬。
由于流域内地形起伏较大，且河道坡降大、槽蓄能力小，易在暴雨过后形成陡涨陡落型洪水，3天洪量占一次洪水总量的80%。
冬季一般从11月份开始到次年3月末或4月初结束，期间主要以降雪为主，在3月末或4月初开始融雪，河流开始解冻，形成桃汛。
桓仁水电站以发电为主，兼顾防洪灌溉等综合效益。
拦河坝按千年一遇洪水设计，万年一遇洪水校核。
下游桓仁镇的防洪标准为五十年一遇，库区耕地补偿与居民迁移，分别为十年一遇与三十年一遇洪水标准，相应坝前水位为303米与304.7米。
课程设计任务：已知桓仁流域场次洪水“19770731”的降雨径流信息如附件1所示，求：（1）洪水“19770731”的次洪径流深。
桓仁流域的退水曲线如附件2所示，基流为60m3/s。
（2）求该场次洪水的预报流量过程，要求产流模型使用三层蒸散发模型（模型参数组合不少于4种），汇流模型使用单位线汇流。
所需单位线见附件1。
（3）绘制所有结果的降雨流量过程图（包括实际流量和预报流量），并计算确定性系数、径流深的相对误差和绝对误差、洪峰的相对误差和绝对误差、峰现时间误差。
水文预报课程设计
退水曲线成果图
将半透明纸上的退水曲线去拟合自己所画的Q-T 过程线中的退水过程，找到一条拟 合度最好的，将这条退水曲线绘制在Q-T过程线上。
流域径流深---R的计算
从开始涨洪时间点开始，计算整个涨洪过程的洪量----W.
注意：整个洪水过程要计算完毕，不能遗漏，计算的洪流应该是绘制过退水 曲线之后的Q-T过程线围成的闭合图形的面积。 计算过程中要注意时间单位的转化。1.3代表的是一点半
时段单位线修匀
修匀要求：
1,修匀后计算峰值与实测值差值在允许范围内。 2,修匀后的单位线产生的径流要仍为10mm。
修匀方法：找到最敏感和最迟钝的单位线上的数据。
误差分析
实测误差。实测雨量误差。 模型误差。概化模型。 误差大----批判---误差产生原因分析
回流检验 到峰时间 洪峰流量 ，时间误差，径流深 确定性系数（查手册）
P-T过程线的绘制
1，降雨数据整理计算。 根据各个测站实测资料，依据各测站数据权重大小计算得平 均时段雨量
2,注意是先降雨后产流。降雨时间提前1-2天。 3,根据所得平均时段降雨将降雨过程线（柱状图）与Q-T曲线绘制 在一张坐标纸上。
成果图
退水曲线的绘制
1,数据选择。 每个小组选择一次退水过程---所设计年的一次退水。一般选取一年内退水退到流量 最小值的一次退水。画圆圈部分。
注意：绘制退水过程线时，坐标纸横纵坐标与Q-T过程线保持 一致。 将12条退水曲线绘制到半透明纸上。 绘制方法：保持横坐标对齐，左右移动半透明纸，使得每条退 水曲线的尾部重合。 遇到的问题： 1，所绘制退水曲线有交叉---选取退水时段长，流量最小的那条保 留。 2，某些退水曲线没有重合部分，形成原因，退水结束时流量过大。
水文预报课程设计报告材料
水文预报课程设计报告学院：资源与环境学院专业：水文与水资源工程年级：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 概述流域水文模型可分为物理模型、概念性模型和系统模型。
在水文预报中，概念性模型和系统模型应用较多，此处主要介绍概念性流域水文模型。
概念性流域水文模型属于数学模型，它与物理模型相比，具有许多优点：一是它的所有条件均可由原型观测资料直接给出，不受比尺的限制，即数学模型无相似律问题；二是它的边界条件及其它条件可严格控制，也可随时按实际需要改变；三是它的通用型较强，只要研制出一种适用的应用软件，就可用来解决不同的实际问题；四是它具有理想的抗干扰性能，只要条件不变，重复模拟可以得到相同的结果，不会因人、因地而异；五是它的研制费用相对较低。
水文预报课程设计报告
水文预报课程设计报告学院：资源与环境学院专业：水文与水资源工程年级：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 概述流域水文模型可分为物理模型、概念性模型和系统模型。
在水文预报中，概念性模型和系统模型应用较多,此处主要介绍概念性流域水文模型。
概念性流域水文模型属于数学模型,它与物理模型相比，具有许多优点：一是它的所有条件均可由原型观测资料直接给出，不受比尺的限制，即数学模型无相似律问题；二是它的边界条件及其它条件可严格控制，也可随时按实际需要改变；三是它的通用型较强,只要研制出一种适用的应用软件,就可用来解决不同的实际问题；四是它具有理想的抗干扰性能,只要条件不变,重复模拟可以得到相同的结果，不会因人、因地而异；五是它的研制费用相对较低。
水文预报课程设计报告书
《水文预报课程设计》姓名：学号：学院：水利与环境学院班级：指导老师：时间：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再用调整的单位线计算地面流量，比较预测与实际值间的误差大小。
水文预报教学设计
水文预报教学设计一、教学目标1.掌握水文预报的基本概念和原理；2.熟悉水文预报的流程和方法；3.理解水文预报在防洪减灾中的作用；4.能够灵活运用水文预报技术进行实际操作。
二、教学内容1. 水文预报的基本概念和原理1.1 水文预报定义1.2 水文预报的基本原理1.3 水文预报的分类及特点2. 水文预报的流程和方法2.1 数据采集与处理2.2 大坝坝前预报2.3 主坝坝后预报2.4 超前预报2.5 定量预报3. 水文预报在防洪减灾中的作用3.1 洪水预警预报3.2 洪水保护预报3.3 水文预报在防洪中的应用4. 水文预报技术的实际操作4.1 软件的操作方法4.2 数据处理方法4.3 编制预报报表三、教学方式本次教学采用讲授、案例分析、实践操作相结合的方式进行。
四、教学方法1. 讲授讲授主要通过课件、PPT等经典的教学方式进行，针对每个教学内容提供详细且系统的介绍和分析，帮助学生掌握水文预报的基本概念、原理及流程、方法和应用等内容。
2. 案例分析通过真实的案例，让学生更加深入地了解水文预报的应用价值及实践过程，帮助学生在学术理论知识上，更加贴近实践。
3. 实践操作通过实验室等方式，让学生实际操作水文预报技术，理解技术的实际操作方法及完全体验技术的精髓之处，帮助学生提高相关技能，为懂得掌握水文预报技术打下坚实的实践基础。
五、教学质量评估针对本次教学，采用多种评估方式，其中包括：课堂答问、实验操作等，提高教学质量和效果。
六、教学时间本次教学时间为两天，每天共18学时，以坚果2名为主。
七、教学材料本次教学提供相关的教材、PPT、计算机备课等材料，以帮助学生更好地掌握相关知识和技能，达到教学目标。
水文预报课程设计
水文预报课程设计
水文预报是水文学中的一个重要分支，涉及到气象学、地质学、工程学等多个领域。
水文预报的目标是对未来一定时间内水质量的预报，帮助人们做出决策，保障水资源的充分利用。
水文预报课程设计可以包括以下步骤:
1. 确定课程的目标和范围：明确本课程的教学目标和教学内容，包括水文预报的基本概念、计算方法、数据处理、水文预报模型的应用等。
2. 收集和整理相关资料：通过学习国内外相关教材、文献和网站，收集和整理水文预报的相关资料，包括水文预报的基本概念、计算方法、数据处理、水文预报模型的应用等。
3. 编写课程教案：根据课程的目标和范围，编写课程教案，包括教学内容、教学方法、教学评价等。
4. 实施教学：按照课程教案进行教学，可以采用讲授、讨论、案例分析等形式进行，帮助学生掌握水文预报的基本概念、计算方法、数据处理、水文预报模型的应用等。
5. 教学评价：对教学效果进行评价，包括学生的反馈、教学内容的满意度、教学方法的满意度等，根据评价结果进行改进和完善。
6. 完成课程设计：完成水文预报课程设计，包括课程教案的编写、教学评价的制定、教学内容的优化等。
以上是水文预报课程设计的一般流程，具体实施时需要根据具体情况进行调整和完善。
水文预报课设报告
当 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
课程设计报告
水文预报课程设计报告课程设计说明书：（一）马斯京根流量演算过程（1）K、x的确定和流量演算方程的确定①K、x值试算表马斯京根法基本原理：马斯京根法基于下面的槽蓄方程：'=+-=（1-1）W K[xI(1x)Q]KQ式中，Q’为示储流量，单位为m3/s, Q’的计算式子如下：'(1x)Q=+-（1-2）Q xIK为蓄量流量关系曲线坡度，K=△W/△Q'； x为流量比重系数。
根据所给资料，利用公式（1-2）可试算得到x的试算表如表1：表1 X值试算表单位：（m3/s）表中，（3）、（4）列为已知，（5）列 iO I q I I-=⨯∑∑∑ （1-3）（6）=（4）+（5）,（7）=（6）-（4）,（8）列中，1112221[(I Q )(I Q )]2W q q ∆=+-++- ,(9)列等于第（8）列依次累计相加，（10）到（13）列为x 的试算值，先分别假定x=0.1，x=0.2，x=0.3，x=0.35，然后根据公式（1-2）计算可得。
（2）试算过程中不同的x 值所对应的Q ’—W 曲线根据四个假定的x 值算出的Q ’和W 可绘制Q ’—W 关系曲线，如图1（a,b,c,d ）所示。
图1（a ），x=0.1图1（a）Q'—W关系曲线-50005000100001500020000W（m 3/s）图1（b）Q'—W关系曲线-50005000100001500020000W（m3/s）Q '(m 3/s )图1（c ），x=0.3图1（c）Q'—W关系曲线-50005000100001500020000W(m 3/s)图1（d）Q'—W关系曲线-50005000100001500020000W（m 3/s）Q (m 3/s )x 能代表水面曲线的形状，反映楔形的大小，还反映河段的调蓄能力和河道的水力特性，x 值会随水位的变化而变化，经试算和分析可知，当x=0.35时，计算分段流量演算时较符合题意，所以取x=0.35。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基本任务任务一：根据已给数据资料及参数（本报告采用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）产流量计算：根据流域特点，产流量计算系根据蓄满产流理论得出的。
蓄满产流，即任一地点上，土壤含水量达田间持水量前，降雨量全部补充土壤含水量，不产流；当土壤蓄满后，其后续降雨量全部产生径流。
a=WMM(1-(1-W0/WM)^(1/(b+1))土壤三层含水量WU WL WD 计算1.PE>0，则产流，即R>0；此时土壤中上中下层含水量会发生变化，如果R ΔW产流，则PE-R后的值就给三层土壤来分，上层的土壤先补充水，如果上层WU>WUM，中层开始补充水，如果WL>WLM，下层开始补充水，如果三层都达到最大值，则土壤水分保持最大值不变。
2.否则不产流，R=0。
即P-E<=0 此时土壤中水分含量随着蒸发减少，上层土壤先开始蒸发，然后中层，最后下层。
产流时，当PE+a<WMM：R=PE+W-WM+WM*(1-(PE+a)/WMM)^(b+1)当PE+a<WMM：R=PE+W-WM（3）水源划分：流域坡地上的降雨产流量因产流过程的条件和运动路径不同，受流域的调蓄作用不同，各径流成分在流量过程线上的反应是不一样的。
各水源分量如下：1）当PE〈=FC时：RD=0.0RG=R2）当PE〉FC时：RG=FC*R/PERD=R-RG（4）汇流计算：根据流域净雨和流域径流单位线，采用卷积的差分形式算出流域出口的流量过程。
地上产生的径流用单位线推流，地下部分用线性水库。
计算公式：QS(I)=RS(I)*UQG(I)=CG*RG(I)*U+(1-CG)*QG(I-1)Q(I)=QS(I)+QG(I)上式中RS即为地表径流Rd.基本数据任务一：1.计算流域为白盆珠水库上游的宝口流域，流域面积553 km2。
2.流域内有四个雨量站：禾多布、马山、高潭、宝口，权中重分别为0.55、0.14、0.33、0.20。
3. 资料文件名：87-92data csv，数据格式为：T(i) Q(i) E(i) P1(i) P2(i) P3(i)P4(i)4.计算参数：任务二：1.计算参数：WM=140,UM=20,LM=60,DM=60IM=0.002,FC=11.0W=110,WU=10,WL=40,WD=60b=0.3,C=0.16 Cg=0.978 Qg=55.3 2.计算基本数据如下表：计算结果任务一：kc优选计算表1991年KC 优选1992年KC 优选1991，1992两年KC 优选由于个人认为KC优选应该以单独一年为单位，因此本程序只是对单独一年的KC 优选做了程序，如果以两年为单位使用范围会很小。
因此，两年共同的KC优选，在本程序中的解决方法是，将单独一年KC从0.9到1.3间隔0.01的值全部带入计算并得出所有R与实测R计算相对误差的值输出并可以导入EXCEL，因此可以把两年所有KC和相对误差导入EXCEL。
通过绘制图表来反映共同KC时，相对误差之间的关系。
以选出尽可能好的KC值。
两年KC优选的图表：可以看出在上图中两年在共同小于1%时有个交点，即KC=1.07 此时对于两年径流的计算，为共同最优KC。
任务二：次洪计算结果此次选出的KC=1.07 为下部汇流计算的基础。
最后绘制时间与降雨，地表流量，地下流量，总流量的次洪流量过程线：计算程序及说明程序代码：（1）全自动单年KC优选结果自动导入EXCEL：Dim i As Integer, P1(366) As Single, Eo(366) As Single, P2(366) As Single, P3(366) As Single, P4(366) As SingleDim Ep(366) As Single, ro As Single, Q(366) As SinglePrivate Sub Command1_Click()Dim ExcelApp As Excel.ApplicationPrint MsgBox("资料导入过称有点慢，请耐心等待30秒左右，此过程根据个人RP不同快慢不同o(∩_∩)o ...", 64, "温馨提示")Set ExcelApp = CreateObject("excel.application")If Combo1.Text = "1991" ThenExcelApp.Workbooks.Open (App.Path & "\1991.csv")Else: ExcelApp.Workbooks.Open (App.Path & "\1992.csv")End IfWith MSFlexGrid1.Rows = 367.Cols = 20.TextMatrix(0, 0) = "T".TextMatrix(0, 1) = "Q".TextMatrix(0, 2) = "Eo".TextMatrix(0, 3) = "P1".TextMatrix(0, 4) = "P2".TextMatrix(0, 5) = "P3".TextMatrix(0, 6) = "P4".TextMatrix(0, 7) = "P".TextMatrix(0, 8) = "WU".TextMatrix(0, 9) = "WL".TextMatrix(0, 10) = "WD".TextMatrix(0, 11) = "W".TextMatrix(0, 12) = "EU".TextMatrix(0, 13) = "EL".TextMatrix(0, 14) = "ED".TextMatrix(0, 15) = "E".TextMatrix(0, 16) = "PE".TextMatrix(0, 17) = "R"i = 0v = 0For R = 1 To .Rows - 1For c = 1 To .Cols.TextMatrix(R, c - 1) = ExcelApp.Sheets(1).Cells(R, c)Nexti = i + 1Q(i) = Val(.TextMatrix(R, 1))Eo(i) = Val(.TextMatrix(R, 2))P1(i) = Val(.TextMatrix(R, 3))P2(i) = Val(.TextMatrix(R, 4))P3(i) = Val(.TextMatrix(R, 5))P4(i) = Val(.TextMatrix(R, 6))v = v + Q(i)Nextro = (v * 3.6 * 24) / 553End WithExcelApp.QuitEnd SubPrivate Sub Command2_Click() ‘开始各参数计算Dim P(366) As Single, Wum As Integer, Wlm As Integer, Wdm As Integer, b As Single, c As Single, kc As Single, fc As SingleDim Eu(366) As Single, El(366) As Single, Ed(366) As Single, E(366) As Single, PE(366) As Single, Wu(367) As Single, Wl(367) As SingleDim Wd(367) As Single, W(367) As Single, Wm As Single, Wmm As Single, a(366) As Double,Dim sum(41) As Single, j As Integer, z As IntegerPrint MsgBox("此次计算是对你选择的年份进行的产流以及所有KC进行径流误差试算，结果将在下面的图表中显示", 64, "温馨提示")Wum = Val(Text1.Text)Wlm = Val(Text2.Text)Wdm = Val(Text3.Text)b = Val(Text4.Text)c = Val(Text5.Text)fc = Val(Text7.Text)Wu(1) = 10Wl(1) = 40Wd(1) = 60Wm = 140Wmm = Wm * (1 + b)sum(1) = 0z = 1j = 1For kc = 0.9 To 1.3 Step 0.01For i = 1 To MSFlexGrid1.Rows - 1Ep(i) = kc * Eo(i)P(i) = 0.33 * P1(i) + 0.14 * P2(i) + 0.33 * P3(i) + 0.2 * P4(i)If Wu(i) + P(i) >= Ep(i) Then ‘三层蒸散发判别计算Eu(i) = Ep(i)ElseEu(i) = Wu(i) + P(i)End IfIf Wu(i) + P(i) >= Ep(i) ThenEl(i) = 0 And Ed(i) = 0ElseIf Wu(i) + P(i) < Ep(i) And Wl(i) >= Wlm * c ThenEl(i) = (Ep(i) - Eu(i)) * Wl(i) / WlmEd(i) = 0ElseIf Wu(i) + P(i) < Ep(i) And Wl(i) >= (Ep(i) - Eu(i)) * c And Wl(i) < Wlm * c Then El(i) = (Ep(i) - Eu(i)) * cEd(i) = 0ElseIf Wu(i) + P(i) < Ep(i) And Wl(i) < (Ep(i) - Eu(i)) * c ThenEl(i) = Wl(i)Ed(i) = (Ep(i) - Eu(i)) * c - El(i)End IfE(i) = Eu(i) + El(i) + Ed(i)PE(i) = P(i) - E(i)W(i) = Wu(i) + Wd(i) + Wl(i)a(i) = Wmm * (1 - (1 - W(i) / Wm) ^ (1 / (b + 1)))If PE(i) <= 0 ThenR(i) = 0ElseIf PE(i) > 0 And a(i) + PE(i) <= Wmm ThenR(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / Wmm) ^ (b + 1)R(i) = PE(i) + W(i) - WmEnd IfIf Wu(i) + PE(i) - R(i) <= 0 And Wl(i) < 0 Then ‘三层土壤含水量判别Wu(i + 1) = 0Wl(i + 1) = 0Wd(i + 1) = Wd(i) - Ed(i)ElseIf Wu(i) + PE(i) - R(i) <= 0 And Wl(i) > 0 And Wl(i) <= Wlm ThenWu(i + 1) = 0Wl(i + 1) = Wl(i) - El(i)Wd(i + 1) = Wd(i)ElseIf Wu(i) + PE(i) - R(i) > 0 And Wu(i) + PE(i) - R(i) <= Wum ThenWu(i + 1) = Wu(i) + PE(i) - R(i)Wl(i + 1) = Wl(i)Wd(i + 1) = Wd(i)ElseIf Wu(i) + PE(i) - R(i) >= Wum And Wl(i) + Wu(i) - Wum + PE(i) - R(i) < Wlm Then Wu(i + 1) = WumWl(i + 1) = PE(i) + Wu(i) - Wum - R(i) + Wl(i)Wd(i + 1) = Wd(i)ElseIf Wu(i) + PE(i) - R(i) >= Wum And Wu(i) - Wum + PE(i) - R(i) + Wl(i) >= Wlm And Wu(i) - Wum + Wl(i) - Wlm + PE(i) - R(i) < Wdm ThenWu(i + 1) = WumWl(i + 1) = WlmWd(i + 1) = Wd(i) + Wu(i) - Wum + Wl(i) - Wlm + PE(i) - R(i)ElseIf Wu(i) = Wum And Wl(i) + PE(i) - Wlm - R(i) >= Wlm And Wd(i) < Wdm ThenWu(i + 1) = WumWl(i + 1) = WlmWd(i + 1) = PE(i) + Wl(i) - Wlm - R(i) + Wd(i)ElseWu(i + 1) = WumWl(i + 1) = WlmWd(i + 1) = WdmEnd IfEu(i) = Format(Eu(i), "0.0")El(i) = Format(El(i), "0.0")Ed(i) = Format(Ed(i), "0.0")E(i) = Format(E(i), "0.0")PE(i) = Format(PE(i), "0.0")Wu(i) = Format(Wu(i), "0.0")Wl(i) = Format(Wl(i), "0.0")Wd(i) = Format(Wd(i), "0.0")W(i) = Format(W(i), "0.0")R(i) = Format(R(i), "0.0")MSFlexGrid1.TextMatrix(i, 7) = P(i)MSFlexGrid1.TextMatrix(i, 8) = Wu(i)MSFlexGrid1.TextMatrix(i, 10) = Wd(i)MSFlexGrid1.TextMatrix(i, 11) = W(i)MSFlexGrid1.TextMatrix(i, 12) = Eu(i)MSFlexGrid1.TextMatrix(i, 13) = El(i)MSFlexGrid1.TextMatrix(i, 14) = Ed(i)MSFlexGrid1.TextMatrix(i, 15) = E(i)MSFlexGrid1.TextMatrix(i, 16) = PE(i)MSFlexGrid1.TextMatrix(i, 17) = R(i)sum(j) = sum(j) + R(i)Next iMSFlexGrid2.TextMatrix(0, 0) = "kc"MSFlexGrid2.TextMatrix(0, 1) = "实测Ro"MSFlexGrid2.TextMatrix(0, 2) = "计算R"MSFlexGrid2.TextMatrix(0, 3) = "绝对误差"MSFlexGrid2.TextMatrix(0, 4) = "相对误差%"MSFlexGrid2.TextMatrix(z, 0) = Val(kc)MSFlexGrid2.TextMatrix(z, 1) = Val(ro)MSFlexGrid2.TextMatrix(z, 2) = Val(sum(j))MSFlexGrid2.TextMatrix(z, 3) = Val(sum(j) - ro)MSFlexGrid2.TextMatrix(z, 4) = Abs(Val(((sum(j) - ro) / ro) * 100))j = j + 1z = z + 1Next kcEnd SubPrivate Sub Command3_Click()Dim z As Integer, c As Double, T As Double, s As DoubleDim a As Stringa = MsgBox("此按钮进行的是单独一年的最优KC，如果要进行多年KC优选，请点击“导入EXCEL”将多年的误差在EXCEL进行绘图分析，优选！", 0 + 64, "温馨提示")If a = 1 ThenT = 5For z = 1 To 40If MSFlexGrid2.TextMatrix(z, 4) < T ThenT = MSFlexGrid2.TextMatrix(z, 4)s = zc = 0.89 + 0.01 * sEnd IfNext zText9.Text = Str(T)Text10.Text = Str(c)Private Sub Command4_Click()Dim a As StringDim xlApp As Excel.ApplicationDim xlBook As Excel.WorkbookDim xlsheet As Excel.Worksheeta = MsgBox("如果你要继续，请认真阅读下面内容，如果你点击“确定”会把你选择年份的所有KC和误差结果导入EXCEL", 1 + 48, "提示信息")If a = 1 ThenSet xlApp = CreateObject("Excel.Application")xlApp.Visible = TrueSet xlBook = xlApp.Workbooks.AddSet xlsheet = xlBook.Worksheets(1)For i = 1 To 42For j = 1 To 5xlsheet.Cells(i, j) = MSFlexGrid2.TextMatrix(i - 1, j - 1)NextNextElseExit SubEnd IfEnd SubPrivate Sub Command5_Click()Dim a As Stringa = MsgBox("感谢你的使用，此程序所有解释权归水文二班殷钊所有，谢谢！", 1 + 64, "温馨提示")If a = 1 ThenUnload MeElseExit SubEnd IfEnd SubPrivate Sub Command6_Click()任务一.Hide任务二.ShowEnd SubPrivate Sub Form_Activate()Dim a As Stringa = MsgBox("如果你使用的电脑EXCEL版本不是2010可能会导致无法运行下列代码，请运行exe 程序！本程序自动试算所有可能KC并自动优选出最优值。