暴雨强度公式参数估计及其应用.
- 格式:doc
- 大小:66.50 KB
- 文档页数:12
第23卷第1期2000年3月南京气象学院学报Jou rnal of N an jing In stitu te of M eteo ro logy. 23N o . 1V o lM ar . 2000文章编号:100022022(2000 20120063205暴雨强度公式参数估计及其应用3顾骏强徐集云陈海燕黄建萍(浙江省气候中心, 杭州310017摘要:, 计, 提高了公式参数估计的精度。
关键词:非线性模型; 暴雨强度公式; 中图分类号:P 333. 2A, 它直接影响到排水工程的投资预算[1]和可靠性(GBJ 14287 》规定, 在进行城市排水管网设计时, 雨水管。
城市排水工程的可靠性与采用的暴雨强度公式有直接的关系, 而排水工程直接影响到城市防灾减灾的功能和城市环境。
因此, 编制各地暴雨强度公式具有重要的现实意义。
国家规范给定的暴雨强度公式形式为一非线性模型, 过去对公式中的参数通过经验值、解析图解或摆试法确定, 然后利用线性最小二乘法确定公式中的其他参数[2]。
文献[3, 4]利用最优化方法解非线性超定方程组, 也有用加速遗传算法[5]对暴雨强度公式参数估计进行研究, 但具体应用的对象均为暴雨强度公式的参数估计。
本文利用最优化问题的求解方案, 较好地解决非线性模型的参数求解问题, 实现的模式参数求解方法具有广泛的适用性。
在编程中, 充分考虑计算机数值计算能力, 有效地控制了迭代误差, 提高了暴雨强度公式参数估计精度。
1非线性模型参数估计方法暴雨强度公式是一非线性模型, 对于非线性系统数值求解通常采用迭代法。
非线性方程的一般形式为:给定m 个函数f i (i =1, 2, …, m , f i 包含n 个未知参数变量, 即f 1(x 1, x 2, …, x n =0……n ≤m 。
f m (x 1, x 2, …, x n =(1参数的估计问题, 就是设法确定n 个参数x 1, x 2, …, x n , 使实际上的特性函数f 在Ξ的某一区[Α, Β]上逼近规定的特性函数Ω(Ξ , 并令其达到指定的精确度Ε, 即(2 m ax f (x 1, x 2, …, x n ; Ξ -Ω(Ξ ≤Ε。
Α≤Ξ≤Β3收稿日期:1999206204; 改回日期:1999211230基金项目:浙江省重大科研项目(991103364 资助作者简介:顾骏强, 男, 上海市人, 1962年12月生, 高级工程师。
主要研究方向:应用气象, 自然灾害学64南京气象学院学报23卷为了方便求解, 将其离散化; 在区间[Α, Β]上取m 个Ξ值, Ξ1, Ξ2, …, Ξm 。
令Ωi =Ω(Ξi , 则(1 式可近似表示成f (x 1, x 2, …, x n ; Ξi -Ωi ≤Ε; i =1, 2, …, m 。
令m(3 (4Υ(x 1, x 2, …, x n =m in ∑f 2i (x 1, x 2, …, x n 。
i =1(5因此, 问题化成求参数x 1, x 2, …, x n , 使函数Υ为极小。
其中Υi (x 1, x 2, …, x n =f (x 1, x 2, …, x n ;Ξi -Ωi , i =1, 2, …, m , 故非线性方程组的求解问题即化成一个无约束最优化问题。
由于(5 式是参数变量的二次函数, 故其最小值总存在。
用F (x 表示以f 1(x , f 2(x , …, f m (x 为分量的m 维向量, 则Υ(x =F (x T F (x (2因为=2∑f i (x , j =1, , x j x j i =1m(6f i (x x i i =1m(x =…f i (x ∑x n i =1m=2A T F (x 。
(7式中A =[a ij ],a ij =, i =1, 2, …, m ; j =1, 2, …, n 。
x j设x 0是初始值, 令x =x 0+∃x 。
(8 (9 (10将函数f i (x 在点x 0展开为泰勒(T aylo r 级数, 并略去二次项以上的展开项, 得Tf i (x ≈f i (x 0 +f i (x 0 (x -x 0 ; i =1, 2, …, m ,即F (x ≈f 0+A 0∃x 。
其中f 0=f (x 0 , A 0=[];i =1, 2, …, m ; j =1, 2, …, n 。
x j根据最小二乘法, 求∃x 使‖A 0x +f 0‖=极小。
因为T T‖A 0∃x +f 0‖2=‖f 0‖2+2f T 0A 0∃x +∃x A 0A 0∃x 。
(11 (12 (13所以(7 式可化成A 0A 0∃x +A 0f 0=0,TT则∃x =--1T(A T A 0f 0。
0A 0在求出∃x 后, 不直接用x 0+∃x 作为解的第一次近似, 而把∃x 当作寻查方向。
从x 0点出发, 沿该方向用一维寻查方法求函数Υ(x 的极小点, 即(14 Υ(x 0+t ∃x =极小。
设此时t 为t 0, 则x 1=x 0+t ∃x 作为第一近似, 然后用x 1代替x 0重复前述步骤, 直到得出满足精度要求的解为止。
1期顾骏强等:暴雨强度公式参数估计及其应用65在算法的编程中, 对于差商求偏导数, 利用计算机语言中具有高位精度的扩展型实数作处理, 以使误差控制在10-7~10-5范围内, 保证数值迭代的计算精度。
2暴雨强度公式参数估计根据我国现行室外排水设计规范(GBJ 14287 [1], 暴雨强度公式的形式为I =(t +B N(15式中I 为暴雨强度(mm 编m in , t 为降水历时(m in , T 为重现期, A 1、C 、B 、N 为地方性参数。
制暴雨强度公式就是根据当地降雨的自记雨量记录推求(15 中的参数。
当确定某一重现期时, 即(15 式的分子部分为一参数, 假设为A , 此时公式称做暴雨强度分公式; 相应地(15 式也可称做暴雨强度总公式。
根据前面的讨论, 解问题化为多元函数的极小化问题, NL S 。
根据式(15 , 可将(1 式中的f i 具体表示为:f i (A 1, C , B , N =I i -; I =1, 2, …, m 。
(t i (16m 为暴雨强度样本数, 即n =3, f i 中包含三个参数。
通过前述[2~6], 在文献[5, 6], 这里利用NL S 数值求解方法对该例子进行计算, 。
表1给出暴雨强度公式各种参数估计方法对同一算例计算结果的绝对均方差, 可见, NL S法对参数估计计算精度较高, 表明该方法对公式的拟合效果有明显提高。
另外, 根据文献[4]提供的数据, 利用NL S 方法对上海地区的暴雨强度总公式参数作估计, 并与文献中给出的总公式编制误差做比较(表2 。
从表2的绝对均方差不难看出, N SL 法的计算误差比文献介绍计算结果偏小17. 5%, 进一步说明本文实现对暴雨强度公式参数估计方法, 其计算精度比过去使用的方法有明显提高。
表1各参数估计方法拟合效果T able 1 F itting effects of different param eter esti m ati on m ethods 重现期(a0. 250. 330. 501. 002. 005. 0010. 0020. 0050. 00100. 00传统法0. 0130. 0130. 0130. 0160. 0190. 0310. 0400. 0480. 0570. 0720. 0322优选法[6]0. 00750. 00730. 00750. 01200. 01900. 03000. 03900. 04600. 05500. 06900. 0292A GA 法[5]0. 00770. 00820. 00850. 01180. 01910. 03020. 05650. 04570. 05470. 06720. 0309NL S 法0. 00670. 00820. 00820. 01370. 01970. 02710. 03880. 04420. 04430. 06730. 0278平均66南京气象学院学报表2上海暴雨强度总公式参数的两种估计方法比较T able 2Comparison betw een tw o param eters esti m ati on m ethods into rrential rain intensity general fo r m ula in ShanghaiA123卷C B N Ρ0. 06640. 0547文献[4]N SL 法9. 45009. 72700. 71890. 67795. 54005. 96540. 65140. 6521注:Ρ表示绝对均方差3公式参数估计的应用根据绍兴市1963~1995年逐日降水自记记录, 按5、10、15、20、、90、120m in 9个降水历时, 8, 每个历时都可得到264个样本, 通过计算机处理, , 并取每个历时的前132, 0. 25a 暴雨强度理论值; 同时, 这种。
利用指数分布函数对各历时抽样频数分布进行拟合, (表3 。
表3绍兴各历时不同重现期的暴雨强度T able 3To rrential rainfall intensity of different durati onsand return peri ods in Shaoxing重现期(a0. 250. 330. 501. 002. 003. 005. 0010. 0020. 0050. 00100. 005m in 1. 3631. 5081. 7232. 0832. 4432. 6542. 9193. 2803. 6404. 1164. 47610m in 1. 1091. 2161. 3751. 6421. 9082. 0642. 2602. 5262. 7933. 1453. 41115m in 0. 8950. 9961. 1471. 3991. 6511. 7991. 9852. 2372. 4892. 8223. 07420m in 0. 7550. 8450. 9801. 2051. 4301. 5621. 7281. 9532. 1782. 4762. 70130m in 0. 5730. 6480. 7600. 9471. 1341. 2341. 3811. 5681. 7552. 0032. 19045m in 0. 4260. 4860. 5740. 7220. 8700. 9571. 0661. 2141. 3621. 5581. 70660m in 0. 3440. 3930. 4670. 5890. 7110. 7830. 8730. 9951. 1181. 2791. 402mm m in90m in 120m in 0. 2610. 2960. 3480. 4350. 5220. 5730. 6370. 7240. 8110. 9261. 0120. 2140. 2410. 2830. 3510. 4200. 4600. 5110. 5790. 6480. 7390. 808根据表3中的数据, 利用前述的非线性模型参数估计方法, 可求出各指定重现期下的暴雨强度分公式参数, 及不定重现期时的暴雨强度总公式参数。