水文学中雨强公式参数求解的一种最优化方法
An Optimized Method for Estimating Parameters of the Rainstorm Intensity Formula
-
摘要: 提出了一种客观的、最优化的暴雨强度公式参数估算方法:先将公式线性化, 确定出未知参数b, C取值范围, 给定一个b值 (分公式)、b, C组合 (总公式), 再对雨强-历时-重现期 (i-t-T) 三联表数据进行最小二乘法拟合可得到参数A, n, 以总误差最小为控制条件, 理论上可得到最优的一组参数估算值。并以深圳、武汉两市为例, 进行暴雨强度公式参数估算, 精度高于国家标准要求, 且明显优于对比方法。该法已被编制成计算机软件, 只要输入原始资料就可以很快输出结果, 包括曲线型估计、参数估算、误差分析、图表, 使用极其方便, 可向全国各地推广应用。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 深圳市不同历时 (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)
表 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)
表 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)
表 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)
-
[1] 北京市市政设计院.给水排水设计手册 (第五册).北京:中国建筑工业出版社, 1985:48-87. [2] 上海市建设委员会.中华人民共和国国家标准-室外排水设计规范 (GBJ14-87, 1997年版).北京:中国计划出版社, 1998:4-5;160-161. [3] 毛慧琴, 宋丽莉, 杜尧东.珠江三角洲地区城市暴雨强度公式研究.自然灾害学报, 2003, 12(2):341-345. [4] 戴慎志, 陈践.城市给水排水工程规划.合肥:安徽科技出版, 2001:1-10. [5] 邓培德.城市暴雨公式统计中若干问题.中国给水排水, 1992, 18(2):45-48. http://www.cnki.com.cn/Article/CJFDTOTAL-GSPS199203011.htm [6] 邓培德.暴雨选样与频率分布模型及其应用.给水排水, 1996, 22(2):5-9. http://www.cnki.com.cn/Article/CJFDTOTAL-JZJS199602000.htm [7] 任伯帜, 许仕荣, 王涛.编制现代城市暴雨强度公式的统计方法研究.湖南城建高等专科学校学报, 2001, 10(2):31-33;74. http://www.cnki.com.cn/Article/CJFDTOTAL-HNCG200102011.htm [8] 周胜昔.浙江省城市暴雨强度公式研究课题介绍.浙江建筑, 1996, (6):13-26. http://www.cnki.com.cn/Article/CJFDTOTAL-ZJZU199606008.htm [9] 夏宗尧.编制暴雨强度公式中应用P—Ⅲ曲线与指数曲线的比较.中国给水排水, 1992, 18(2):32-38. http://www.cnki.com.cn/Article/CJFDTOTAL-GSPS199003006.htm [10] 季日臣, 郭晓东, 刘有录.编制兰州市暴雨强度公式中频率曲线的比较.兰州铁道学院学报 (自然科学版), 2002, 21(1):64-66. http://www.cnki.com.cn/Article/CJFDTOTAL-LZTX200201017.htm [11] 顾骏强, 陈海燕.瑞安市暴雨强度概率分布公式参数估计研究.应用气象学报, 2000, 11(3):355-363. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20000352&flag=1 [12] 周玉文, 周胜昔.极大似然法求皮尔逊Ⅲ型分布参数.给水排水, 1997, 23(6):19-21. http://www.cnki.com.cn/Article/CJFDTOTAL-JZJS706.004.htm [13] 常福宣, 丁晶, 姚健.降雨随历时变化标度性质的探讨.长江流域资源与环境, 2002, 11(1):79-82. http://www.cnki.com.cn/Article/CJFDTOTAL-CJLY200201017.htm [14] 周桂明.计算机搜索法推求暴雨强度公式参数.杭州大学学报, 1998, 25(4):85-89. http://www.cnki.com.cn/Article/CJFDTOTAL-HZDX804.015.htm [15] 植石群, 宋而莉, 罗金玲, 等.暴雨强度计算系统及其应用.气象, 2000, 26(6):30-33. http://www.cnki.com.cn/Article/CJFDTOTAL-QXXX200006006.htm [16] 乔华, 张理, 高俊发.西安市暴雨强度公式的推导与研究.西北建筑工程学院学报, 1996, 65-71. http://www.cnki.com.cn/Article/CJFDTOTAL-XBJG604.011.htm [17] 季日臣, 郭晓东, 刘有录.兰州市暴雨强度关系的研究.兰州铁道学院学报 (自然科学版), 2002, 21(6):65-68. http://www.cnki.com.cn/Article/CJFDTOTAL-LZTX200206017.htm [18] 杨开, 程晓如.暴雨强度公式中系数B统计计算法一例.人民长江, 1996, 27(3):16;22. http://www.cnki.com.cn/Article/CJFDTOTAL-RIVE603.009.htm [19] 张爱英, 钱喜镇, 王栋成, 等.短历时暴雨强度公式研究及应用软件介绍.山东气象, 2003, (1):32-33. http://www.cnki.com.cn/Article/CJFDTOTAL-SDQX200301011.htm [20] 王世刚.城市暴雨公式参数优化计算程序.中国给水排水, 1987, 3(4):50-52. http://www.cnki.com.cn/Article/CJFDTOTAL-GSPS198704016.htm [21] 赵建国.迭代法推求暴雨强度公式参数.给水排水, 1997, 23(12):9-12. http://www.cnki.com.cn/Article/CJFDTOTAL-JZJS712.001.htm [22] 顾骏强, 徐集云, 陈海燕, 等.暴雨强度公式参数估计及其应用.南京气象学院学报, 2000, 23(1):63-67. http://www.cnki.com.cn/Article/CJFDTOTAL-NJQX200001009.htm [23] 李树平, 刘遂庆, 黄延林.用麦夸尔特法推求暴雨强度公式参数.给水排水, 1999, 25(2):26-28. http://www.cnki.com.cn/Article/CJFDTOTAL-JZJS902.008.htm [24] 任伯帜.带因子—迭代法求解暴雨强度公式参数.中国给水排水, 2002, 18(2):40-42. http://www.cnki.com.cn/Article/CJFDTOTAL-GSPS200202027.htm [25] 任伯帜, 许仕荣.基于Marqardt—Hartley法及其在求解城市暴雨强度公式参数中的应用研究.给水排水, 2002, 29(3):96-100. http://www.cnki.com.cn/Article/CJFDTOTAL-HNDX200203020.htm [26] 杨晓华, 金菊良, 张国桃.加速遗传算法及其在暴雨强度公式参数优化中的应用.自然灾害学报, 1998, 7(3):71-75. http://www.cnki.com.cn/Article/CJFDTOTAL-ZRZH803.011.htm [27] 朱颖元.暴雨强度公式参数率定方法.中国给水排水, 1999, 15(7):32-33. http://www.cnki.com.cn/Article/CJFDTOTAL-GSPS907.010.htm