暴雨强度的总公式与分公式，并对其计算精度进行分析比较

摘要：本文采用规范规定的年多个样法进行雨样选取处理，通过P-Ⅲ、Weibull分布模型对雨样进行频率分析优选出误差最小的理论分布推求出P-i-t表，应用遗传算法、拟牛顿算法、高斯牛顿法^[1]推求暴雨强度的总公式与分公式，并对其计算精度进行分析比较。结果表明：不管总公式还是分公式其计算精度均能满足《室外排水设计规范》^[3]的要求，但分公式的计算精度明显高于总公式。
关键词：  暴雨强度    暴雨公式    精度比较
1暴雨公式推求的步骤
  城市暴雨强度公式作为城市排水设计的基础公式，其正确与否将直接关系到城市基础设施建设的科学性。推求暴雨强度公式的工作的程序按时间先后可以分为三部分：第一部分为选取雨样；第二部分为频率分析；第三部分为推求公式。
1.1 选取雨样
  城市暴雨资料的收集、暴雨强度资料的选样与统计方法及与之相关的频率分布线型选择是城市暴雨强度公式制定过程中的前段工作，也是极其关健的环节，因为它直接影响编制暴雨强度公式所需的P-i-t经验数据表的质量对暴雨公式的精度有相当大的影响。选样是从现有的记录中合理选择若干个数值以组成一个样本，来作为频率分析的依据，因此暴雨强度资料的选样工作极其重要，在选样过程中要充分做到每个单元具有一致性和独立性，所组成的样本具有代表性和足够的可行性，从而在此基础上认真分析研究，再选择与之适应的频率分布线。
  根据《室外排水设计规范》(GB50014-2006) ^[2]的规定，本文采用年多个样法选样，完整收集1953年～2006年间某市的降雨量，对每年分别挑选并读取8场最大的暴雨中每5分钟、10分钟、15分钟、20分钟、30分钟、45分钟、60分钟、90分钟、120分钟九个时段的降雨量数据进行整理，然后统一排序，从大到小取资料年数4倍的最大值作为统计的基础。
1.2 频率分析
  频率分布线型的选择对城市暴雨强度公式的精确制定起到保证作用。因为它直接关系到编制公式所需的重现期-暴雨强度-历时(即P-i-t)经验数据表的可靠性。本文采用矩法、拟牛顿法、适线法、极大似然法、遗传算法对P-Ⅲ分布模型^[3]，以及采用最小二乘法、遗传算法、拟牛顿法对Weibull分布模型进行拟合分析优选出误差最小的理论分布。拟合分析的结果见表1：
表1 频率分析误差比较表

根据表1可以看出相对于Weibull最小二乘推求的P-i-t数据的误差最小，总绝对方差为0.0308，总相对方差为2.95%，均满足于规范规定的精度要求。因此将相对于Weibull最小二乘推求的P-i-t数据作为暴雨强度公式推求的基础数据。
1.3 推求公式
   利用以上频率分析所得的最优P-i-t数据，分别采用遗传算法、拟牛顿算法、高斯牛顿法推求暴雨强度公式。

   暴雨强度公式的推求根据我国的排水规范，在推求单一重现期暴雨强度公式(以下简称分公式)时采用的公式形式为：

   而在推求各重现期统一公式(以下简称总公式)时，采用的公式形式为：

 式中：
i —— 降雨强度(mm/min)；
P —— 重现期(a)；
t —— 降雨历时(min)；
A₁、C、b、n——地方参数。
   采用遗传算法、拟牛顿算法、高斯牛顿法推求暴雨强度分公式的平均绝对方差和平均相对方差见表2。表3为采用以上三种方法推求暴雨强度总公式的参数成果。
表2 分公式平均误差比较

从表2中可以看出用拟牛顿法推求的暴雨强度分公式的平均误差最小

   平均绝对方差为0.0180 mm/min，平均相对方差为1.04%，均满足于规范规定的精度要求。因此将用拟牛顿法推求的暴雨强度分公式的参数成果作为研究对象具体数据见表3。同理，根据表4中的参数可知，采用遗传算法推求暴雨强度总公式的误差最小。总绝对方差为0.0448 mm/min，总相对方差为3.26%，均满足于规范规定的精度要求。因此将用遗传算法推求的暴雨强度总公式的参数成果作为比较对象。

表 3 拟牛顿法推求分公式参数成果