暴雨强度公式编制方法研究 (1).doc

《暴雨强度公式编制方法研究 (1).doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、城市暴雨强度公式编制方法研究顾骏强陈海燕徐集云（浙江省气候中心杭州310017）摘要根据暴雨强度统计资料，利用概率分布摸式对不同历时暴雨强度的频数分布进行分析，通过比较选择适合公式编制的概率模式，并通过解非线性超定方程组改进公式参数估计的精度。关键词暴雨强度公式概率模式参数估计引言城市短历时暴雨强度公式是设计城市雨水排水管道及小汇水流域排水管渠的主要依据。虽然国家有关规范和规范性手册对暴雨强度公式的编制与应用作了许多规定，但是在具体工作中，往往出现有悖于规范的作法，这给市政工程建设造成极大的隐患。随着社会的进步，法

2、规的健全，科技的进步，要求市政工程的规划与设计建立在科学的基础上。因此，不仅每座城市都应拥有暴雨强度公式，而且要求公式的精度更高，误差更小。这就需要对公式编制涉及到的一系列问题作研究，以便提高公式编制精度，改进公式编制方法。我国现行排水设计规范中的暴雨强度公式形式为：i＝(1)式中A1、C、b、n为地方性参数，T为重现期，t为暴雨历时(min)，i为暴雨强度(mm/min)。暴雨历时一般采用5、10、15、20、30、45、60、90、120min共9个历时。重现期一般取0.25、0.33、0.5、1、2、3、5、

3、10、20、50、100年。编制暴雨强度公式的过程是根据当地自记雨量资料，采用数理统计方法确定暴雨强度抽样分布模式，计算各历时不同重现期理论暴雨强度值，然后通过适当的方法估计公式中的参数，使公式能反映当地的暴雨规律。国内对暴雨强度抽样频数分布拟合，较多地采用指数分布和皮尔逊（Pearson）Ⅲ型分布模式。关于公式参数估计，由经验值与计算相结合以及图解与计算相结合的方法，逐步被数值求解法取代。1暴雨强度抽样分布的拟合1.1概率分布模式1.1.1耿贝尔(Gumbel)分布耿贝尔分布的概率密度函数和分布函数分别为：f(x

4、)=（2）Ｆ(x)＝(3)式中y为抽样xi(i=1,2,……,N)的函数，N为抽样数。耿贝尔推得：y＝a×(xi－b)，其中a>0为尺度参数，b是分布密度的众数。通过积分，可以求得y的数学期望值和均方差分别为：my=0.577216(4)σy=π/(5)因此，由y=a×(xi－b)，可以得到：a=σx(6)b=mx－my(7)这里mx和σx分别为xi的数学期望和均方差。分布函数参数确定以后，通过下式求对应于重现期T的分位数：xT=b－ln(－ln(1－))（8)1.1.2指数分布指数分布的概率密度函数和分布函数分别

5、为：f(x)＝(9)Ｆ(x)＝＝1－(10)式中xi(i=1,2,……,N)为抽样资料，N为抽样数。a表示分布曲线离散程度的参数，r为分布曲线的下限。分布函数超过概率：Ｐ＝1－Ｆ(x)＝(11)故对应于重现期T的分位数：xT=r－lnT(12)令1/a=α，利用抽样的经验概率分布和最小二乘法，可以计算出系数α、r的估计，并通过(12)式计算各重现期的分位数。1.1.3皮尔逊(Pearson)Ⅲ型分布皮尔逊Ⅲ型分布密度函数和分布函数为：ｆ(x)＝(13)Ｆ(x)＝(14)式中抽样变量xi≥0，参数α、β≥0分别表示概

6、率曲线的形状和尺度参数，a为概率曲线起点与序列零点的距离。用皮尔逊型Ⅲ分布拟合暴雨强度抽样资料，首先要对参数α、β、a作估计。这里采用矩法和极大似然法。根据矩法得到皮尔逊Ⅲ型的数学期望、方差和偏度系数分别为：＝a+攭σ2＝(15)Cs＝2/由此得到：α=β=(16)a=－由于高阶矩误差较大，因此，矩法对皮尔逊Ⅲ型参数的估计被认为误差偏大。为了提高估计精度，可用权函数对参数进行修正，即：=α=β=(17)a=－其中和为样本的均值和均方差，而E=(xi－)ψ(xi)G=(xi－)2ψ(xi)(18)ψ(x)=ψ(x)是

7、正态概率分布密度函数。参数估计加权后，减少了由于抽样频数分布尾部取值较少而造成的误差，提高了计算精度，这就是所谓的“权函数法”。由极大似然法函数得到以下似然方程组：lnβ－ψ(α)+ln(xi－a)=0α=β(－a)(19)β=(α－1)上式消去β得到：α=ψ(α)=ln+ln(xi－a)(20)由于ψ(α)=－0.577216+(α－1)(i(i－1+α))－1，式中s为一正整数。因此，通过迭代可求得参数α、β的估计。皮尔逊Ⅲ型概率分布模式参数确定后，可根据达到某一精度要求的拟合曲线进行各概率或重现期对应的分位数

8、的估算。可以证明皮尔逊Ⅲ型分布曲线某一概率或重现期对应的分位数只依赖于参数α或偏度系数Cs。因此，给定一个参数α或Cs值后，就可通过曲线离均系数表，查算出概率P或重现期T对应的分位数。若对Cs作适当调整，使拟合误差达到最小，这便是所谓的“适线法”。1.1.4韦伯(Weibull)分布三参数韦伯分布的概率密度函数和分布函数为：f(x)=x>a(21)F(x)=