1、ICS07.060CCSA 47东 北 三 省 一 区 区 域 协 同 地 方 标 准DB21/T 20019-2025DB22/T 20019-2025DB23/T 20019-2025DB15/T 20019-2025城市暴雨强度公式编制技术规范Specification for establishment of urban rainstorm intensity formula2025-09-23 发布2025-12-24 实施辽 宁 省 市 场 监 督 管 理 局吉 林 省 市 场 监 督 管 理 局黑 龙 江 省 市 场 监 督 管 理 局内蒙古自治区市场监督管理局发布DB21/T
2、20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-2025I目次前言.III1范围.12规范性引用文件.13术语和定义.14样本选取和资料处理.25暴雨强度公式参数推算方法.36长历时重现期雨量计算方法.47暴雨强度公式适用性分析.58常用查算图表编制.5附录 A（规范性）暴雨强度公式参数求解方法.7附录 B（资料性）常用查算图表格式.10DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-2025III前言为推进东北三省一区经济社会高质量发展
3、辽宁省市场监督管理局、吉林省市场监督管理厅、黑龙江省市场监督管理局、内蒙古自治区市场监督管理局共同组织制定本区域标准，在东北三省一区内适用。本文件按照GB/T 1.12020标准化工作导则第1部分：标准化文件的结构和起草规则的规定起草。本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由内蒙古自治区气象标准化技术委员会（SAM/TC 23）提出。本文件由辽宁省气象局、吉林省气象局、黑龙江省气象局、内蒙古自治区气象标准化技术委员会（SAM/TC 23）归口。本文件起草单位：内蒙古自治区气候中心、吉林省气候中心、辽宁省气候中心、黑龙江省气候中心。本文件主要起草人：马玉峰、
4、刘诗梦、王曼霏、孙玉、赵艳丽、屠志鑫、东方、姜忠宝、朱玲、徐永清、刘柏鑫、顾正强、于秀晶、李金泽、徐艺璇、王婷、范彦斌。DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-20251城市暴雨强度公式编制技术规范1范围本文件规定了短历时城市暴雨强度公式编制的样本选取和资料处理、频率分布曲线拟合、公式参数推算、长历时重现期雨量计算、公式适用性分析以及相关图表编制等方面的技术要求。本文件适用于城镇、工业园区、经济开发区的暴雨强度公式编制工作，其他地区亦可参照执行。2规范性引用文件下列文件中的内容通过文中的规范性引用而构
5、成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T 35228地面气象观测规范 降水量GB 50014室外排水设计标准QX/T 118气象观测资料质量控制 地面QX/T 341降雨过程强度等级3术语和定义下列术语和定义适用于本文件。3.1降雨量rainfall某一时段内从天空降到地面上的雨，未经蒸发、渗透、流失而在水平面上积聚的深度。注：以毫米（mm）为单位。来源：QX/T 341-2016,2.1,有修改3.2降雨强度rainfall intensity某一历时内单位时间（每分钟或每小时）的
6、降雨量。注：以毫米/分钟（mm/min）或毫米/小时（mm/h）为单位。3.3降雨历时rainfall duration降雨过程中任意连续时段的持续时间。注：以分钟（min）为单位。3.4暴雨重现期rainstorm return period某一强度的暴雨重复出现的统计平均时间间隔。注：以年（a）为单位。3.5暴雨强度公式rainstorm intensity formulaDB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-20252表征暴雨强度、重现期和降雨历时三者之间函数关系的数学表达式，反映一定频率的暴
7、雨在规定时段的最不利时程分配的平均强度。3.6短历时降雨short duration rainfall任意连续时段持续时间不大于180 min降雨。3.7长历时降雨long duration rainfall任意连续时段持续时间大于180 min且不大于1440 min的降雨。4样本选取和资料处理4.1站点选择4.1.1应满足区域代表性、历史数据连续性和一致性要求。4.1.2优先选择国家气象站作为代表站，不应采用数个站点的数据混合样本。4.1.3选择站点在资料年限内发生迁址、雨量记录仪更换时，应对降雨资料进行一致性检验，并对降雨资料的代表性和一致性进行论证和说明。4.1.4宜使用最近 30 年
8、以上的降雨资料。4.1.5应至少选择一个代表站，若城区规模大或狭长分散、海拔高差较大、降雨分布不均的城市，宜选择多个代表站分区域编制暴雨强度公式。4.2原始数据4.2.1原始降雨资料采用逐分钟记录的基础数据，包括以自记纸形式记录的逐分钟降雨资料和自动气象站逐分钟降雨资料。4.2.2应按照 GB/T 35228 和 QX/T 118 的规定，对原始资料进行数据质量检查、审核，其中历史降雨自记记录资料应采用信息化提取并经省级以上气象主管部门审核后的标准分钟降水数据集资料。4.2.3由自记纸记录和自动记录混合形成的原始降雨基础数据，应进行数据的一致性分析。4.2.4应收集并给出基础资料涉及的各个观测
9、站的降雨观测历史沿革。4.3统计样本4.3.1短历时降雨应采用 5 min、10 min、15 min、20 min、30 min、45 min、60 min、90 min、120min、150 min、180 min 共 11 个历时的样本。4.3.2长历时降雨应在短历时降雨 11 个历时基础上，增加 240 min、360 min、540 min、720 min 和1440 min 降雨历时的样本。4.3.3统计样本选取，应采用逐分钟滑动统计法，选取各降雨历时雨量逐年最大值的降雨场次，记录选定降雨场次过程开始时间及逐分钟降雨量，作为各降雨历时的有效降雨资料样本。同场次降雨过程中，同一种历时
10、降雨不可以交叉，不受日、月界的限制（但不跨年）。4.3.4应按照 GB 50014 的规定，采用年最大值法进行暴雨样本选样，即选取降雨资料样本中各降雨历时雨量的逐年最大值，作为有效暴雨资料样本。注：自记雨量记录在10年至20年的地区，暴雨样本选样方法宜采用年多个样法。4.4特殊问题处理4.4.1宜选择 1440 min 降雨量大于或等于 10.0 mm 为指标，筛选出场次降雨过程，纳入统计样本遴选范围。DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-20253注：各地降雨量年际变化大，尤其部分干旱地区在少雨年
11、份难以选出足够样本，一些短历时（如5min）强降水过程可能丢失，可根据当地降雨特点，选择1440 min降雨量大于或等于5.0 mm为降雨过程选取指标。4.4.2选取样本时，偶遇特大值时，应进行专门分析确认。宜采用样本均值加上 3 倍样本标准差作为判定标准，大于判定标准可确定为特大值。4.4.3不宜舍弃特大值，可采用皮尔逊型分布函数进行人工适线，妥善处理离群数据。确需舍弃特大值，宜单独进行说明，论证舍弃特大值可能导致的不确定性。5暴雨强度公式参数推算方法5.1暴雨强度公式的表达形式暴雨强度公式的表达形式见公式 1。=1671 1+log10+(1)式中：q 暴雨强度，单位为升/(秒公顷)L/(
12、shm2)；A1 雨力参数；C 雨力变动参数；P 暴雨重现期，单位为年（a）；t 降雨历时，单位为分钟（min）；b 降雨历时修正参数，单位为分钟（min）；n 暴雨衰减指数。注1：本文件中暴雨强度公式指按照GB 50014的规定，采用5 min、10 min、15 min、20 min、30 min、45 min、60min、90 min、120 min、150 min、180 min共11个降雨历时统计样本的暴雨强度公式。注2：室外排水设计采用的暴雨强度q（单位为L/(shm2)）与气象行业常用的降雨强度i（单位为mm/min）的换算关系为：q=167i。5.2暴雨强度重现期和频率计算暴雨
13、强度重现期应按2年、3年、5年、10年、20年、30年、50年、100年统计，相对应的频率为50%、33.3%、20%、10%、5%、3.3%、2%、1%。5.3理论频率分布曲线拟合5.3.1选取的各历时统计样本，应用理论频率曲线进行趋势性拟合调整。理论频率曲线类型一般采用皮尔逊型分布曲线、耿贝尔分布曲线和指数分布曲线。根据确定的频率分布曲线，得出降雨强度、重现期和降雨历时三者的关系，即 i-P-t 关系表。注：本文件所列举的皮尔逊型分布曲线、耿贝尔分布曲线和指数分布曲线的函数表达式可参照相关文献，不作具体规定。5.3.2应根据当地降雨特点，进行多种频率分布函数的拟合试验，从中选取拟合效果较好
14、的理论频率曲线函数。5.3.3应结合当地统计分析降雨频率的经验，在控制拟合精度的同时，保证频率曲线的总体协调。5.3.4宜对拟合确定的频率分布曲线适线外延至 1%频率（100 年重现期）。DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-202545.4暴雨强度公式参数求解方法5.4.1暴雨强度公式可以采用符合数理统计理论基础的数学优化方法求解参数，包括求解非线性方程的方法或最优化方法率定暴雨强度公式参数。推荐采用数值逼近法和高斯-牛顿法，按照附录 A。5.4.2将求得的各参数最优解代入公式(1)，即得到当地的
15、暴雨强度公式。5.4.3在编制重现期 2 年至 100 年范围的暴雨强度公式时，应重点保证重现期 2 年至 20 年区间的拟合精度。5.5精度检验5.5.1应按照 GB 50014 的要求，计算暴雨强度公式的平均均方根误差和相对均方根误差。计算重现期在 2 年至 20 年时，平均均方根误差不宜大于 0.05 mm/min，或相对均方根误差不宜大于 5%。5.5.2平均均方根误差和相对均方根误差分别由公式（2）和公式（3）计算。=1=12(1)=1=12 100%(2)式中：平均均方根误差，单位为毫米/分钟（mm/min）；相对均方根误差，单位为%；理论降雨量，单位为毫米（mm）；Ri-P-t曲
16、线确定的降雨量，单位为毫米（mm）；t降雨历时，单位为分钟（min）；i第i个样本；n样本总数。6长历时重现期雨量计算方法6.1频率曲线拟合方法和误差检验6.1.1频率分布曲线拟合应基于选取的统计样本，采用经验频率曲线或理论频率曲线进行趋势性拟合调整，宜选择使用的理论频率曲线有皮尔逊-型分布、耿贝尔型分布和指数型分布，通过多种频率分布函数拟合试验，选取拟合误差较小的理论频率曲线。注：当多种频率曲线拟合误差均满足精度要求时，一般推荐皮尔逊-型分布函数。6.1.2用理论频率曲线拟合频率分布时，应结合当地统计分析降雨频率的经验，宜在控制拟合精度的同时，保证频率曲线的总体协调。6.1.3宜根据不同资料
17、时段样本的频率曲线，拟合 2 年至 100 年误差结果，重点保证重现期 20 年至100 年区间的拟合精度，并利用各历时的样本与曲线拟合误差结果，结合降雨年代特征等分析结果，选取适宜的资料时段。6.1.4频率曲线拟合误差计算按公式（2）和公式（3）。6.2重现期雨量计算6.2.1应根据频率分布曲线，推算不同重现期下的长历时降雨量。6.2.2宜通过频率分布曲线计算的降雨量与当地已发布的暴雨强度公式计算结果进行对比，分析不同DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-20255重现期下降雨量在城市防涝中的适宜
18、性。6.2.3宜对比分析频率分布曲线计算的降雨量与当地暴雨洪水图集、水文图集中长历时重现期雨量结果，分析不同重现期雨量在城市排水规划中的适用性。7暴雨强度公式适用性分析7.1时间分布特征分析7.1.1绘制各历时暴雨（纳入了暴雨公式编制的样本）出现日（次）数的年际变化图，分析各历时暴雨的逐年或年代变化特征，以描述各历时暴雨年际变化特征。7.1.2可以把各年降雨数据入选各历时有效暴雨资料样本的占比作为评价指标，分析城市暴雨的时间变化特征。7.1.3可以各历时降雨量的前 10 个最大值在各年代的分布特征，分析城市强暴雨的时间变化特征。7.1.4根据当地暴雨变化特征，采用代表性气象站点不同时期基础资料
19、编制不同时段暴雨强度公式，并对计算结果进行差异性分析，论证降雨资料年限选择的合理性。7.2空间分布特征分析7.2.1对于地形地貌差异较大或涵盖范围大的城市，应补充利用本地加密观测或高精度再分析降雨资料，分析本地强降雨的区域分布特征，为编制的暴雨强度公式的适用范围提供依据。注：应用高精度再分析降雨资料时，宜进行精度检验和降尺度处理。7.2.2可对比代表性站点与加密观测站点或高精度再分析降雨资料，绘制部分典型历时的暴雨事件发生次数、最大降雨量分布图，以及代表性站点与加密观测站点特征参数的比值或差值等，分析暴雨的区域分布特征。7.2.3城区规模大或狭长分散、海拔高差较大、降雨分布不均的城市，可开展
20、不同分区暴雨强度公式的差异性分析，确定暴雨分区的合理性。7.3社会影响分析7.3.1依据空间分布特征分析结论，说明是否需要划分暴雨分区、站点代表性以及新编暴雨强度公式的适应地域范围。7.3.2依据气候变化、城市化发展和强暴雨时间分布特征，说明样本年限确定的合理性和适用性，并提出当地暴雨强度公式修编年限。7.3.3将同一气候区域内的城市已经批准实行的暴雨强度公式作为参考，评估当地新编暴雨强度公式的区域适用性。7.3.4新编暴雨强度公式与现用暴雨强度公式的对比分析，应注意代表性站点选择、样本年限、选样方法、线型拟合、公式参数计算方法等产生的差异，并给出分析结论。7.4分析结论和建议7.4.1提出推
21、荐编制成果（含公式、图表），并说明其适用范围和条件。7.4.2宜每 5 年对短历时暴雨强度公式和长历时重现期雨量进行误差复核或修编。7.4.3若在最近时段内出现了超过历史极值的暴雨，应对短历时暴雨强度公式和长历时重现期雨量进行复核修正。注：最近时段具体指暴雨强度公式通过审定并公布后至最近发生的超过历史极值暴雨年份的这段时期。8常用查算图表编制DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-202568.1常用图表和格式8.1.1常用查算图表包括：暴雨强度查算表、暴雨强度曲线图和不同历时的暴雨强度频率曲线图。8
22、1.2常用查算图表格式见附录 B。8.2图表编制说明8.2.1采用数值逼近法可拟合出暴雨强度总公式、单一重现期暴雨强度公式和重现期区间参数公式，三者计算得出常用查算图表的计算精度和数值略有差异，在编制暴雨强度公式及图表时，应规定适用范围和条件，并标注说明。8.2.2编制暴雨强度查算图表时，应明确标注变量单位，并注意不同单位之间的换算关系。DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-20257AA附录A（规范性）暴雨强度公式参数求解方法A.1数据逼近法A.1.1单一重现期公式参数拟合暴雨强度公式的表达形式
23、见A.1：=1671 1+log10+(A.1)令=1671 1+log10 ，则得到一个简化的表达式，即为单一重现期公式见 A.2：=+(A.2)式中：A 雨力参数，即不同重现期下1min的设计降水量（mm）。对其两边取对数，并令：y=lnq，b0=ln A，b1=-n，x=ln(t+b)，则公式进一步简化为一个一元线性方程形式：=0+1(A.3)采用最小二乘法，可求出（A.3）式中的b0和b1，则可求出A、n。由于（A.2）式中的 b 也是未知数，采用“数值逼近法”来处理：先给定一个 b 值，采用最小二乘法进行计算，得出相应的 A、n 值，同时求出其均方根误差，不断调整 b 值，直至使其值
24、达到最小时，得到最为合理的 A、b、n 值。同理，以此方法，可将 11 个降雨历时的单一重现期暴雨强度公式逐个推算出来。A.1.2区间参数公式拟合为计算两个单一重现期之间的暴雨强度，引入重现期区间参数见公式A.4：=1+2 ln +(A.4)式中：y A、b、n参数中的任一个；P 重现期；C 常数。首先把2年至100年划分为二个区间：2年至10年；：10年至100年。将A、b、n 代入区间参数公式见A.5、A.6、A.7：=1+2 ln +(A.5)=1+2 ln +(A.6)=1+2 ln +(A.7)DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20
25、019-2025;DB15/T 20019-20258上面三式中，A、b、n 和 P 是已知数，A1、A2、CA、b1、b2、Cb 及 n1、n2、Cn都是未知数。根据求得的单一重现期 P 下的 A、b、n 值，利用二分搜索法和最小二乘法，可解得未知数 A1、A2、CA、b1、b2、Cb及 n1、n2、Cn，从而得到 2 年至 10 年和 10 年至 100 年二个区间的 A、b、n 值，即可得到 2 年至 100 年任意区间单一重现期暴雨强度公式。A.1.3总公式参数拟合对公式（A.1）的两端求对数见公式A.8：ln =ln 1671+ln 1+log10 ln +(A.8)设 =ln ln
26、 1+log10 ，0=ln 1671，b1=-n,=ln +，则上式可见公式A.9：=0+1(A.9)以最小二乘法求出b0、b1，从而可求出A1、n以及q(拟合值)，同时求出总公式的均方根误差见A.10：=10=101=1 2(A.10)式中：均方根误差，单位毫米/分钟（mm/min）；q实际降雨强度，单位为毫米/分钟（mm/min）；q拟合降雨强度，单位为毫米/分钟（mm/min）；m短历时取值为 10，长历时取值为 16；m0 8 个重现期；第个降雨历时；第个重现期。取使最小的一组参数A1、b、n，即为最佳拟合参数。A.2高斯-牛顿法高斯-牛顿法的基本思想是使用泰勒级数展开式去近似地代替
27、非线性回归模型，然后通过多次迭代，多次修正回归系数，使回归系数不断逼近非线性回归模型的最佳回归系数，最后使原模型的残差平方和达到最小。初始值的选择有三种方法：a）根据以往的经验选定初始值；b）用分段法求出初始值；c）对于可线性化的非线性回归模型，通过线性变换，用最小平方法求出初始值。泰勒级数展开式。设一般的非线性回归模型见公式A.11、A.12、A.13：=,+(A.11)=1,2,(A.12)=1,2,(A.13)式中：f 一般函数；x 可以是单个自变量，也可以是 r 个自变量；DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB1
28、5/T 20019-20259p 维参数向量；随机误差项，且 0,2。设对 y 和 x 通过 m 次观测，得到 m 组数据：1,2,;i=1,2,m的两个角标，第一个代表观测序号，i=1m；第二个代表自变量序号j=1r。求“最小二乘”拟合曲线，就是求的估计值，使得公式（A.14）为最小。=1 2=1 ,2(A.14)对于非线性模型，无法直接求“最小二乘”解，若把它对待定参数=1,2,在 0=1 0,2 0,0处展开成只包括一次项的泰勒级数，从而使非线性模型线性化。为方便起见，以下用 代替,，导出公式（A.14）中 为最小的参数递推公式，并写成矩阵形式见公式 A.15、A.16：+1=+1 (A
29、15)式中：k 递推次数。=1 11 21 2 12 22 1 2=(A.16)=1,2,，=1 ,2 ,。用公式（A.15）从 0 开始，一步步递推下去，直到 收敛稳定，即 +1 的值小于或等于预先指定的小正数，从而得到的估计值。DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-202510BB附录B（资料性）常用查算图表格式B.1暴雨强度查算表暴雨强度查算样表见表B.1。表 B.1 暴雨强度查算表（样表）P=2(a)t：minq：L/(shm2)tqtqtqtqtqtqtq1335.02527127.25
30、05382.9257962.76510551.03913143.29515737.7642312.74128124.5535481.8798062.20210650.68313243.04815837.5823293.52029121.9805580.8638161.65010750.33313342.80515937.4024276.75530119.5235679.8748261.10910849.98913442.56516037.2255261.98931117.1745778.9128360.57810949.65013542.32816137.0496248.87732114.92
31、65877.9758460.05811049.31513642.09316236.8757237.14633112.7725977.0638559.54711148.98613741.86216336.7028226.58434110.7066076.1758659.04611248.66213841.63316436.5329217.01935108.7236175.3098758.55411348.34213941.40816536.36310208.31236106.8186274.4648858.07211448.02714041.18416636.19611200.34937104.
32、9866373.6418957.59811547.71714140.96416736.03112193.03638103.2226472.8389057.13311647.41114240.74616835.86713186.29439101.5246572.0549156.67611747.10914340.53116935.70514180.0574099.8866671.2899256.22711846.81214440.31817035.54515174.2694198.3076770.5419355.78611946.51814540.10817135.38616168.881429
33、6.7826869.8119455.35312046.22914639.90017235.22917163.8514395.3106969.0989554.92712145.94414739.69517335.07318159.1444493.8867068.4009654.50912245.66314839.49217434.91919154.7304592.5097167.7199754.09712345.38614939.29117534.76720150.5804691.1777267.0529853.69312445.11215039.09317634.61621146.671478
34、9.8867366.4009953.29512544.84215138.89717734.46622142.9824888.6367465.76110052.90312644.57515238.70317834.31823139.4944987.4247565.13710152.51812744.31315338.51117934.17124136.1915086.2497664.52510252.14012844.05315438.32118034.02625133.0585185.1087763.92610351.76712943.79715538.13326130.0825284.0007863.34010451.40013043.54415637.947B.2暴雨强度曲线图暴雨强度曲线样图见图B.1。DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-202511图 B.1 暴雨强度曲线图（样图）DB21/T 20019-2025;DB22/T 20019-2025;DB23/T 20019-2025;DB15/T 20019-202512B.3不同历时的暴雨强度频率曲线图不同历时的暴雨强度频率曲线样图见图B.2。图 B.2 不同历时的暴雨强度频率曲线图（样图）
