应用气象学报  2007, 18 (2): 237-241   PDF    
引用本文
陈正洪, 王海军, 张小丽. 水文学中雨强公式参数求解的一种最优化方法[J]. 应用气象学报, 2007, 18(2): 237-241.
Chen Zhenghong, Wang Haijun, Zhang Xiaoli. An Optimized Method for Estimating Parameters of the Rainstorm Intensity Formula[J]. Journal of Applied Meteorological Science, 2007, 18(2): 237-241  
水文学中雨强公式参数求解的一种最优化方法
陈正洪1,2, 王海军3, 张小丽4     
1. 武汉区域气候中心, 武汉 430074;
2. 中国气象局武汉暴雨研究所, 武汉 430074;
3. 湖北省气象信息与技术保障中心, 武汉 430074;
4. 深圳市气象台, 深圳 518008
2006-07-24 收到, 2006-11-28 收到修改稿.
资助项目: 深圳市气象局科研开发项目 (2003-3-03)、湖北省气象局科技基金 (2004Y07) 和中国气象局气候变化专项 (CCSF2006-31) 共同资助    
摘要: 提出了一种客观的、最优化的暴雨强度公式参数估算方法:先将公式线性化, 确定出未知参数b, C取值范围, 给定一个b值 (分公式)、b, C组合 (总公式), 再对雨强-历时-重现期 (i-t-T) 三联表数据进行最小二乘法拟合可得到参数A, n, 以总误差最小为控制条件, 理论上可得到最优的一组参数估算值。并以深圳、武汉两市为例, 进行暴雨强度公式参数估算, 精度高于国家标准要求, 且明显优于对比方法。该法已被编制成计算机软件, 只要输入原始资料就可以很快输出结果, 包括曲线型估计、参数估算、误差分析、图表, 使用极其方便, 可向全国各地推广应用。
关键词: 暴雨强度公式    线性化    最优化    误差控制    
An Optimized Method for Estimating Parameters of the Rainstorm Intensity Formula
Chen Zhenghong1,2, Wang Haijun3, Zhang Xiaoli4     
1. Wuhan Regional Climate Center, Wuhan 430074;
2. Wuhan Institute of Heavy Rain, CMA, Wuhan 430074;
3. Hubei Meteorological Information and Technological Support Center, Wuhan 430074;
4. Shenzhen Meteorological Observatory, Shenzhen 518008
Abstract: In hydrology, the form of the rainstorm intensity formula is known as the nonlinear one with excess factors, and its parameters are not easy to be solved with normal methods, which makes the method design of parameters estimate and the elimination of the error from parameters estimate the most crucial. For a long time, many methods to estimate the parameters of the formula have been put forward by hydrological and meteorological experts that have exerted a profound influence on urban water drainage design. But all the methods above have certain error from objective error by personal judging, approximate supposed error, or missing error from skip-over searching. An objective and optimized method for estimating the parameters of the rainstorm intensity formula is put forward that can avoid the above errors skillfully. First, the non-linear formula is linearized and the scope of the parameters b and C is decided. Second, the parameters A and n can be determined based on known rainstorm intensity, duration and return period (i-t-T) tri-relation table by least regression method after the value of b (single period formula) or a combination of the value of b and C(multi-period formula) is given. Then the formula with least error is the optimum one. This method is used in Shenzhen and Wuhan separately, and the accuracy of the formula can meet the requirement of the national standard and is superior to the compared methods. A software has been designed according to the new method and can be used easily and popularized nationwide. If the original data are put in, the estimating result of the curve pattern, formula parameters, error, figures and table can be obtained quickly.
Key words: rainstorm intensity formula     linearize     optimize     error control    
引言

根据国家《给水排水设计手册》[1]和《室外排水设计规范 (GBJ14-87) 》[2]规定, 暴雨强度是设计水库大坝高度, 确定公路、铁路涵洞直径以及城市雨污分流管道系统的关键技术参数, 也是设计防洪及水利工程设施中的重要指标, 而这些工程排水的可靠与否和采用的暴雨强度公式的准确性和精度有直接关系。

暴雨强度公式是反映一定重现期、历时下的平均暴雨强度, 有许多种经验公式, 在我国一般采用如下公式[1-2] :

(1)

式 (1) 中i, t, T均为变量; i为暴雨强度 (单位: mm/min); T为重现期 (单位:a), 取值范围为0.25~100 a; t为降雨历时 (单位:min), 取值范围为1~120 min。重现期越长、历时越短, 暴雨强度就越大。而A1, C, b, n是与地方暴雨特性有关且需求解的参数 (b, n亦称气候参数): A1雨力参数, 即重现期为1 a时的1 min设计降雨量 (单位:mm), C为雨力变动参数, b为降雨历时修正参数, 即对暴雨强度公式两边求对数后能使曲线化成直线所加的一个时间常数 (单位:min), n为暴雨衰减指数, 与重现期有关。

(2)

式 (2) 中, A为雨力参数, 即不同重现期下的1 min设计降雨量 (单位: mm)。

可见暴雨强度公式为已知关系式的超定非线性方程, 总公式、分公式各有4个和3个参数, 常规方法将无能为力, 所以参数估计方法设计和减少估算误差尤为关键。长期以来水文气象工作者不懈探索, 提出许多参数估算方法, 在城市排水设计中发挥了重要的作用。但这些方法仍然存在较大的人为判断误差 (经验法、图解法), 一定的近似假设误差 (麦夸尔特法, 遗传法, 加速遗传法, Marqardt-Hartley法, 解超定非线性方程组), 以及跳跃搜索中遗漏误差 (二分搜索法或黄金分割法)。本文在以上工作基础上提出了一种客观、对参数进行全组合、可找到最小误差的参数估算方法, 可巧妙避开以上3类问题的出现。虽然占用机时较多, 如对总公式, 若参数均保留3位小数, 有近亿种组合, 在PIV/2.8 G微机上约需机时10 min左右, 但在计算速度不断加快的今天已不是问题。

1 方法与误差分析

首先对式 (1) 和 (2) 进行线性化处理, 再推导出未知参数b与C的取值范围, 给定一个b值 (分公式) b与C组合 (总公式), 对已求出的雨强-历时-重现期 (i-t-T) 三联表数据, 进行最小二乘法拟合便可得到参数A, n, 以总误差最小为控制条件, 理论上可得到最优的一组参数估算值。采用均方根误差 (σ) 或相对均方根误差 (f) 并参照国家标准[2]进行误差评判。

1.1 确定b和C的范围

根据大量文献, 从我国200多个暴雨强度公式[1, 3-27]中可见, b, C是有一定范围的, 最大值分别为46.4(江苏无锡), 1.537(河南济源)。另外1 + ClgT必须大于0, 而T=0.25时, lgT=-0.602, 所以必须有C < 1.666。

1.2 对分公式线性化及最优化目标控制

对分公式, b取 (0, 50.000), 以0.001为间隔, 共50000种情景, 对分公式两端求对数:

(3)
(4)

通过最小二乘法求出b0, b1, 从而可求出A, n以及i′ (拟合值), 同时求出均方根误差:

(5)

以其为目标函数, 取使σ最小的一组参数, 同时计算出相对均方根误差。

同理, 也可求出i与i′的相关系数r, 取使r最大的一组参数。

1.3 对总公式线性化及最优化目标控制

对总公式, b取 (0, 50.000), C取 (0, 1.666), 以0.01为间隔, 共5000×166=830000种组合, b取 (0, 50.000), C取 (0, 1.666), 以0.001为间隔, 共50000×1666=83300000种组合对总公式两端求对数:

(6)
(7)

通过最小二乘法求出b0, b1, 从而可求出A1, n以及i′, 同时求出总的误差:

(8)

以其为目标函数, 取使σ最小的一组参数。其中也可m0=8, 即前8个重现期 (0.25 a, …, 10 a)。并将总公式分解到分公式 (代入不同的T便可), 算出分公式均方根误差和相对均方根误差。

(9)

同理, 也可求出i与i′的相关系数r, 取使r最大的一组参数。

1.4 误差标准

按照国家标准[2]规定, 暴雨强度公式参数估算误差, 以均方差σ≤0.05 mm为主要衡量指标, 对于深圳、武汉这样降雨强度大的地方, 也可采用另一规定:平均相对误差f≤5%, 或适当放宽条件为均方差在多数重现期下 (0.25 a, …10 a, 共8个重现期) 不高于0.05 mm, 或平均相对误差不高于5%。

1.5 试验资料

利用指数分布, 根据深圳市1954—2003年50年、武汉市1961—1995年35年降水资料序列求出两地各自的雨强-历时-重现期 (i-t-T) 三联表数据, 见表 1①和表 2。可见深圳的雨强比武汉大。

①陈正洪, 张海军, 王小丽, 等.深圳市新一代暴雨强度公式的编制与创新研究———技术报告. 2005.

表 1 深圳市不同历时 (t)、重现期 (T) 对应的雨强 (指数分布, 单位: mm/min) Table 1 Rainstorm intensity with duration (t) and return period (T) in Shenzhen (exponential distribution, unit: mm/min)

表 2 武汉市不同历时 (t)、重现期 (T) 对应的雨强 (指数分布, 单位: mm/ min) Table 2 Rainstorm intensity with duration (t) and return period (T) in Wuhan (exponential distribution, unit: mm/min)

2 试验结果分析 2.1 分公式

根据深圳、武汉的i-t-T资料, 采用最优法求出两地各自的暴雨强度分公式各参数和误差, 见表 3、表 4。

表 3 最优法所求深圳市暴雨强度分公式参数和误差一览表 (指数分布) Table 3 The value of the parameters and error of the rainstorm intensity formula for single return period calculated by optimized method in Shenzhen (exponential distribution)

表 4 最优法所求武汉市暴雨强度分公式参数和误差一览表 (指数分布) Table 4 The value of the parameters and error of the rainstorm intensity formula for single return period calculated by optimized method in Wuhan (exponential distribution)

可见, 两地重现期为20年以下的各分公式的绝对误差均在0.05 mm以下, 最小只有0.0122 mm; 重现期为50年、100年两种情况也只略大于临界值, 最大也只有0.0667 mm。两地相对误差则全部在2.62%以下。至于绝对误差、相对误差平均值, 在深圳分别为0.0324 mm, 1.717%, 在武汉分别为0.0334 mm, 2.284%, 结果相当理想。另外还发现, 两地分公式的绝对误差、相对误差均随重现期增加而增大。深圳市不同重现期下暴雨强度随历时变化曲线见图 1(武汉市图略)。

图 1. 深圳市不同重现期的暴雨强度随历时变化曲线 (指数分布, 从上到下对应重现期为100 a, 50 a, ..., 0.25 a) Fig 1. The curve of rainstorm intensity with duration in different return period in Shenzhen (exponential distribution, the return period: 100 a, 50 a, ..., 0.25 a from top to bottom)

2.2 总公式

采用最优法求出了两地各自的暴雨强度总公式各参数,

(10)
(11)

以上两公式误差见表 5。

表 5 最优法所求深圳、武汉暴雨强度总公式误差一览表 (指数分布) Table 5 The value of the error of the rainstorm intensity formula for any return period calculated by optimized method in Shenzhen and Wuhan (exponential distribution)

由表 5可见, 在深圳, 重现期为10年以下公式绝对误差不高于0.05 mm标准的仅4次 (4/8), 但相对误差不高于5%标准的则有6次 (6/8);在武汉, 重现期10年以下公式绝对误差则全部达标 (8/8), 相对误差达标的也有6次 (6/8)。深圳绝对误差平均不达标, 但平均相对误差可达标; 在武汉, 平均绝对误差、相对误差均达标, 完全可以满足国家标准的要求。

2.3 与其他方法计算结果比较

同时计算出二分搜索法 (黄金分割法) [16, 19]推算出总公式参数后的回代误差 (表 6), 并与表 5进行对比。可见, 二分搜索法在两地求出的暴雨强度公式的绝对误差均不能达标; 至于相对误差, 在深圳重现期10年以下公式的相对误差个数 (5/8) 及平均值刚好可达标, 在武汉虽然重现期10年以下公式的相对误差个数不达标 (4/8), 但所有重现期下公式的相对误差个数可达标 (7/11), 平均值仍不能达标。无疑, 二分搜索法整体效果仍不及最优法。

表 6 二分搜索法 (黄金分割法) 所求深圳、武汉暴雨强度总公式误差一览表 (指数分布) Table 6 The value of the error of the rainstorm intensity formula for any return period calculated by two-part searching method in Shenzhen and Wuhan (exponential distribution)

3 结语

1) 在前人工作基础上, 考虑计算条件的极大改善, 提出了一种简单、有效的暴雨强度公式参数的求解方法, 较好地解决了暴雨强度公式参数的估算问题。

2) 以深圳和武汉多年短历时暴雨样本资料为例, 应用该算法求取了两地暴雨强度公式参数, 结果表明两地公式精度均高于国家规范要求, 尤其是较好地控制了大雨地区总公式的误差。

参考文献
[1] 北京市市政设计院. 给水排水设计手册 (第五册). 北京: 中国建筑工业出版社, 1985: 48-87.
[2] 上海市建设委员会. 中华人民共和国国家标准-室外排水设计规范 (GBJ14-87, 1997年版). 北京: 中国计划出版社, 1998: 4-5.
[3] 毛慧琴, 宋丽莉, 杜尧东. 珠江三角洲地区城市暴雨强度公式研究. 自然灾害学报, 2003, 12, (2): 341–345.
[4] 戴慎志, 陈践. 城市给水排水工程规划. 合肥: 安徽科技出版, 2001: 1-10.
[5] 邓培德. 城市暴雨公式统计中若干问题. 中国给水排水, 1992, 18, (2): 45–48.
[6] 邓培德. 暴雨选样与频率分布模型及其应用. 给水排水, 1996, 22, (2): 5–9.
[7] 任伯帜, 许仕荣, 王涛. 编制现代城市暴雨强度公式的统计方法研究. 湖南城建高等专科学校学报, 2001, 10, (2): 31–33;74.
[8] 周胜昔. 浙江省城市暴雨强度公式研究课题介绍. 浙江建筑, 1996, (6): 13–26.
[9] 夏宗尧. 编制暴雨强度公式中应用P—Ⅲ曲线与指数曲线的比较. 中国给水排水, 1992, 18, (2): 32–38.
[10] 季日臣, 郭晓东, 刘有录. 编制兰州市暴雨强度公式中频率曲线的比较. 兰州铁道学院学报 (自然科学版), 2002, 21, (1): 64–66.
[11] 顾骏强, 陈海燕. 瑞安市暴雨强度概率分布公式参数估计研究. 应用气象学报, 2000, 11, (3): 355–363.
[12] 周玉文, 周胜昔. 极大似然法求皮尔逊Ⅲ型分布参数. 给水排水, 1997, 23, (6): 19–21.
[13] 常福宣, 丁晶, 姚健. 降雨随历时变化标度性质的探讨. 长江流域资源与环境, 2002, 11, (1): 79–82.
[14] 周桂明. 计算机搜索法推求暴雨强度公式参数. 杭州大学学报, 1998, 25, (4): 85–89.
[15] 植石群, 宋而莉, 罗金玲, 等. 暴雨强度计算系统及其应用. 气象, 2000, 26, (6): 30–33.
[16] 乔华, 张理, 高俊发. 西安市暴雨强度公式的推导与研究. 西北建筑工程学院学报, 1996: 65–71.
[17] 季日臣, 郭晓东, 刘有录. 兰州市暴雨强度关系的研究. 兰州铁道学院学报 (自然科学版), 2002, 21, (6): 65–68.
[18] 杨开, 程晓如. 暴雨强度公式中系数B统计计算法一例. 人民长江, 1996, 27, (3): 16;22.
[19] 张爱英, 钱喜镇, 王栋成, 等. 短历时暴雨强度公式研究及应用软件介绍. 山东气象, 2003, (1): 32–33.
[20] 王世刚. 城市暴雨公式参数优化计算程序. 中国给水排水, 1987, 3, (4): 50–52.
[21] 赵建国. 迭代法推求暴雨强度公式参数. 给水排水, 1997, 23, (12): 9–12.
[22] 顾骏强, 徐集云, 陈海燕, 等. 暴雨强度公式参数估计及其应用. 南京气象学院学报, 2000, 23, (1): 63–67.
[23] 李树平, 刘遂庆, 黄延林. 用麦夸尔特法推求暴雨强度公式参数. 给水排水, 1999, 25, (2): 26–28.
[24] 任伯帜. 带因子—迭代法求解暴雨强度公式参数. 中国给水排水, 2002, 18, (2): 40–42.
[25] 任伯帜, 许仕荣. 基于Marqardt—Hartley法及其在求解城市暴雨强度公式参数中的应用研究. 给水排水, 2002, 29, (3): 96–100.
[26] 杨晓华, 金菊良, 张国桃. 加速遗传算法及其在暴雨强度公式参数优化中的应用. 自然灾害学报, 1998, 7, (3): 71–75.
[27] 朱颖元. 暴雨强度公式参数率定方法. 中国给水排水, 1999, 15, (7): 32–33.
水文学中雨强公式参数求解的一种最优化方法
陈正洪, 王海军, 张小丽