前言 水是人类生活与生产劳动所必需的物质。水文学是研究 自然界各水体的运动、变化和分布规律的一门学科。随着工、 农业的迅速发展和人口的急剧增长，人们对水的需求日益增 加，因此对水文学的要求也就更高。 讨算机是门新生科学，其FORTRAN语言是目前国内 外科技计算中使用最广泛的一种计算语言。计算机技术已广 泛用于水文领域.在解决水文实际问题和进行水文规律研究 上都起着重要的作用，大大促进了水文学学科的发展。 为了使计算机技术更多更好地在水文领域发挥作用，便 于实现联机预报，提高洪水预报的精度和增长洪水预报的有 效预见期.为水利水电工程的设计、调度管理提供方便，特编 著《常用水文预报算法和计算程序》一书。 该书主要介绍一些国内常用的水文预报算法，包播流域 降雨径流模型、纳须时段单位线、河道汇流模型、水库调洪计 算等，并提供了它们相应的FORTRAN语言子程序。为了使 程序具有灵活和通用性，在程序中没有使用任何库软件，语言 尽可能简化和通俗。可用FORTRAN66或FORTRAN77在 各类小型机和微机上编译运行。 本书共分十章.第一章至第五章分别介绍了国内常用的 几种流域降雨径流模型它们是河海大学的新安江（三水源） .1· 模型与陕北模、美国的萨克拉门托模型、日本的水箱模型， 以及降雨径流经验相关（P～P。~R）模型。以上五个模型的柜 应程字都是一个单元面积上的降雨径流计算模型，输入是单 元面积上的平均降雨量和蒸发能力序列，输出除了P~P～ R是单元面积上的净雨过程外，其它都是单元出口处的流量 过程。各程字都是遵照源文编制的，并宜尽可能使之容易调 用。不过，新安江（三水源）模型的三水源划分有别于参考文献 【1]这是因为参考文献（1】中水源划分公式有误，本书作了改 正。球北模型中采问了流域下渗能力分布曲线，使之可以用于 较大的流域，并具有实价值。第六章是水库调洪计算。它是 由水库入库流量序列，根据初始库水位和水库在一次洪水期 间的实际或计划泄水建筑物启闭情况，进行水库调洪计算，求 出该调度方式下的水库出库流量、库水位、库蓄水量序列。该 程序适用于水库的实时调度。第七章是马斯京根法。主要介 绍马斯京根法（以下简称马法）在国内几种不同的应用形式， 即基本马法、分段冯法、变参数马法、非线性马法、分层马法， 以及各方法的特点和相互关系，并分别给出了各种形式的计 算程序。第八章和第九章分别是纳须单位线法和反馈模拟实 时校正及它们柜应的计算程序。第十章给出了直线插值、抛物 线插值、水位与流量的换算、单位线推流计算、时段单位线的 转换、马斯京根演算参数K与及汇流系数的计算、误差统 计与分析、等时段流量（或水位、雨量）序列的生成与打印过程 线、径流量与输沙量的计算、流量序列的特征值统计等13种 其它常用的水文计算子程序。 为了便于使，前九章除了包括方法介绍、计算程字及说 ·2· 明外，还分别给出了计算实例。 河海大学教授芮孝芳、水利部水利信息中心高级工程师 刘金清、黄河水利委员会教授级高工辛国荣对该书进行审阅， 避免了筑漏：冯立亚司志在本书的编写中做了大量的工作；本 书的出版得到黄污水利委员会及委河务局等单位领导和周围 同志们的大力支持，在此一并表示表心感谢。 由于作者永证所限·不妥之处在所难免·愿请读者批评指 正。 翟家瑞 1994年 6月 ·3· 目 前 言 第一亚 新安江（三水源）模型 (1) 第一节模型结构 (1) 、蒸散发计算 (1) 二、产流量计算… (4) 三、水源划分· (5) 四、汇流计算· (8) 五、不同时段长度的转化 (12) 第二节模型计算程序， (13) 、计算程序 (13) 二、程序说明 (18) 第三节‧实例…· (20) 、主程序 (20) 二、实例说明 (21) 第二章萨克拉门托模型 . (23) 第一节模型结构 (23) 一、产流量计算· (23) 二、蒸散发计算 (26) 三、水分交换计算 (27) 四、下渗计算 (27) 五、镶型参数的意义 (29) 六、几个主要参数对洪水的影响 (30) 第二节模型计算程序 (30) .1: (31) 、计算程序 二、程序说明 (38) 第三节‧实例… (39) 一、主程序 (39) 二、实例说明 (40) 第三章水箱模型… (42) 第节模型结构 (42) 第节模型计算程序 (45) 一、计算程序 (46) 二、程序说明 (47) 第三节‧实例…· (49) 、主程序 (49) 二、实例说明 (50) 第四章 陕北模型…· (52) 第节模型结构· (52) 一、点降雨径流模型 (52) 二、降雨和下垫面空间分布不均匀的考 (55) 三、陕北模型推流步骤 (56) 第二节模型计算程序， (58) 、计算程序 (58) 二、程序说明 (61) 第三节‧实例… (62) 、主程序 (62) 二、实例说明 (63) 第五童降雨径流相关围法 (65) 第一节‧降雨径流相关图 . (65) 、概述 (65) 二、P,的计算 (66) 三、相关图推流计算 (66) 第二节 P～P,～R 计算程序 (68) 、P～P,～R程序 (68) .2: (31) 、计算程序 二、程序说明 (38) 第三节‧实例… (39) 一、主程序 (39) 二、实例说明 (40) 第三章水箱模型… (42) 第节模型结构 (42) 第节模型计算程序 (45) 一、计算程序 (46) 二、程序说明 (47) 第三节‧实例…· (49) 、主程序 (49) 二、实例说明 (50) 第四章 陕北模型…· (52) 第节模型结构· (52) 一、点降雨径流模型 (52) 二、降雨和下垫面空间分布不均匀的考 (55) 三、陕北模型推流步骤 (56) 第二节模型计算程序， (58) 、计算程序 (58) 二、程序说明 (61) 第三节‧实例… (62) 、主程序 (62) 二、实例说明 (63) 第五童降雨径流相关围法 (65) 第一节‧降雨径流相关图 . (65) 、概述 (65) 二、P,的计算 (66) 三、相关图推流计算 (66) 第二节 P～P,～R 计算程序 (68) 、P～P,～R程序 (68) .2: (31) 、计算程序 二、程序说明 (38) 第三节‧实例… (39) 一、主程序 (39) 二、实例说明 (40) 第三章水箱模型… (42) 第节模型结构 (42) 第节模型计算程序 (45) 一、计算程序 (46) 二、程序说明 (47) 第三节‧实例…· (49) 、主程序 (49) 二、实例说明 (50) 第四章 陕北模型…· (52) 第节模型结构· (52) 一、点降雨径流模型 (52) 二、降雨和下垫面空间分布不均匀的考 (55) 三、陕北模型推流步骤 (56) 第二节模型计算程序， (58) 、计算程序 (58) 二、程序说明 (61) 第三节‧实例… (62) 、主程序 (62) 二、实例说明 (63) 第五童降雨径流相关围法 (65) 第一节‧降雨径流相关图 . (65) 、概述 (65) 二、P,的计算 (66) 三、相关图推流计算 (66) 第二节 P～P,～R 计算程序 (68) 、P～P,～R程序 (68) .2: