基于暴雨衰减特性的上海市长历时综合暴雨公式

贾卫红,徐卫忠,李琼芳,周正模

(1. 上海市水务规划设计研究院,上海 200233；2. 上海市气候中心,上海 200030；3. 河海大学水文水资源学院,江苏 南京 210098；4. 长江保护与绿色发展研究院,江苏 南京 210098)

城市洪涝问题历来是一个非常突出的问题[1],不仅影响城市居住环境质量,而且事关城市社会经济的可持续发展[2- 3],是城市发展面临的主要问题,也是水利、市政等相关部门面临的巨大挑战。针对城市洪涝问题,水利部门和市政部门关注的侧重点有所不同,水利部门关注区域性洪涝灾害防御,市政部门关注城市管网排水[4- 5]。两部门采用的设计暴雨量确定方法也不同,市政部门通常利用3 h以内的年最大暴雨所编制的短历时暴雨强度公式来确定设计暴雨量[6],水利部门通常对年最大24 h暴雨进行频率分析来确定设计暴雨量。随着城市化进程的加剧,城市管网排水与区域洪涝灾害防御的相互影响日益凸显,协同考虑排水、治涝与防洪的相互影响,已成为系统治水的必由之路,也是水务一体化的必然之举。由此可见,编制适用于不同历时的综合暴雨公式对城市实现系统治水至关重要。然而,基于水利部门与市政部门现有的设计暴雨核定方法很难获得3～24 h的综合暴雨公式。

目前,国内外对短历时暴雨公式的研究不胜枚举[7- 10],但对长历时暴雨公式的研究屈指可数。国际上更多聚焦暴雨强度- 历时- 频率(IDF)关系对气候变化的响应研究[11- 13],很少涉足长历时暴雨公式的编制；国内的相关研究主要集中在利用频率曲线推算设计暴雨量,以及在不考虑暴雨衰减特性的条件下编制暴雨公式。如严方家等[14]、王捷音[15]、郑长陵和史桂菊[16]采用年最大值选样方法,计算了基于P- Ⅲ频率分布曲线的不同重现期不同长历时的设计暴雨量,未进一步建立长历时暴雨公式,限制了成果的直接推广应用；袁婷婷等[17]、唐颖等[18]基于P- Ⅲ分布线型对年最大值法得到的5～1 440 min不同历时的实测暴雨样本进行频率分析,借助市政部门推荐的短历时暴雨强度公式推导出长历时暴雨公式,未考虑暴雨衰减规律在不同历时之间的差异性,精度有待提高。因此,从暴雨衰减特性出发,编制适用于不同重现期不同历时的综合暴雨公式迫在眉睫。

本文选取易受风暴潮影响、暴雨内涝风险大的上海市作为研究对象,在揭示其暴雨衰减规律的基础上建立1～24 h长历时范围的综合暴雨公式。

1 上海市暴雨特点

上海市地处长江三角洲最前缘、长江流域最下游,境内地势低平、河网密布,为典型的平原感潮河网地区。上海市属于亚热带季风气候区,雨量充沛,年平均降雨量1 191.3 mm,汛期6—9月降雨量占全年的60%以上。上海地区年均暴雨量约占年均降雨量的50%,静止锋暴雨、暖区暴雨、低压暴雨、台风暴雨为主要暴雨类型,约占暴雨总数的83%。静止锋暴雨和台风暴雨的雨时较长、降雨范围广、总雨量大,而暖区暴雨和低压暴雨虽然雨强大,但雨时短、降雨范围小。相关研究[19]表明,长历时、大范围的大暴雨容易造成区域性洪涝灾害,但暴雨历时一般在24 h以内；短历时强暴雨容易造成局部地区严重积水,但历时短于1 h的暴雨造成的危害和影响一般不大。加强历时为1～24 h暴雨特点的研究对上海市防灾减灾至关重要。

上海市各区的暴雨频次和强度虽略有差异(总体上中心城区略高于郊区),但各区气候气象、地形地貌基本一致,其暴雨特性的差异并不十分显著[20],且郊区的城市化程度也在不断提高,各区分别编制暴雨公式意义不大。此外,以徐家汇站为代表的中心城区各历时年最大雨量,除了极个别历时略低于其他区,其余历时均略高于其他区。因此,选取徐家汇站资料来编制全市统一暴雨公式是偏安全的。

2 上海市暴雨衰减规律

暴雨是比较复杂的自然现象,但也有一定的规律性。暴雨强度随历时延长而衰减的规律是暴雨的重要规律之一,一般称之为暴雨强度衰减规律,其具体表现为:暴雨历时越短,雨强衰减程度越小；暴雨历时越长,雨强衰减程度越大。为了厘清暴雨强度与历时的确切关系,国内外长期水文气象研究成果建议在双对数纸上点绘不同重现期的暴雨强度—历时关系图,从而得到一组斜率呈现阶段性变化的折线,不同历时区间的直线斜率就代表该区间暴雨强度衰减指数。绝大部分地区同一重现期的暴雨强度—历时双对数折线会在t=1 h 处出现明显的转折点,部分地区也会在t=6 h处存在转折点[21],因此,可以假设t<1 h的暴雨衰减指数为n1,1 h≤t<6 h时为n2,t≥6 h时为n3,依据不同时间段的直线斜率确定其暴雨强度衰减指数。

为了探明上海市暴雨衰减规律,选取上海市徐家汇站1949—2013年实测降雨资料,统计得到历时分别为5 min、10 min、15 min、20 min、30 min、45 min、60 min、90 min、120 min、180 min、240 min、360 min、720 min、1 440 min的年最大降雨量的逐年系列,并通过了一致性检验,不同历时前几位暴雨量与均值关系分析表明,前3位暴雨量均未超过均值的3倍,且基于上海市暴雨综合统计参数估算得到的各历时前几位暴雨量的重现期均未超过100 a,表明不需要进行特大值处理。选用P- Ⅲ型分布曲线,依据“均值不变”、“频率曲线不交叉”等原则,采用目估适线法进行频率分析(见图1和表1)。由表1可知,变差系数(CV)在0.25～0.50之间变化,没有出现大于0.6的异常情况,且总体上随历时的增加而增加；CS/CV值取3.50(CS为偏差系数),符合中国暴雨统计参数的分布特征与变化规律[21]。

由P- Ⅲ曲线拟合结果计算得到重现期分别为2 a、5 a、10 a、20 a、50 a、100 a的不同历时的设计暴雨量(见表2)。在双对数坐标里点绘暴雨强度与历时关系图,据此可得到不同重现期下不同历时区间的暴雨强度衰减指数(见表3)。

上海市暴雨强度与历时的双对数关系在t=1 h处发生了明显的转折,而在t=6 h处不明显。此外,不同重现期的10 min～1 h历时区间的暴雨强度衰减指数n1差异显著,且呈现随重现期增大而减小的趋势,在坐标中表现为互不平行的直线簇；而不同重现期的1～24 h历时区间的暴雨强度衰减指数n2在0.73～0.75内变化,n3在0.72～0.76内变化,n2与n3的变化区间基本重合,均值都为0.74,在坐标中均呈现出几乎互相平行的直线簇,绝对值差值在±0.02内,变幅约±2.7%。因此,可以确定历时1～24 h的暴雨具有基本一致的衰减规律,其衰减指数为0.74。

图1 徐家汇站1949—2013年不同历时年最大降雨量P- Ⅲ分布曲线

表1 徐家汇站1949—2013年部分历时年最大降雨量P- Ⅲ分布曲线拟合参数

表2 徐家汇站部分历时不同重现期设计暴雨量统计表 mm

表3 上海市不同重现期不同时段暴雨衰减指数统计

3 上海市长历时暴雨公式

3.1 单一重现期暴雨公式

水利部门常用的长历时暴雨强度公式的一般形式[22]如下:

(1)

式中:i为暴雨强度,mm/h；A为雨力,代表t=1 h的平均雨强,mm/h；n为暴雨强度衰减指数；t为暴雨历时,h。

根据上海市暴雨衰减规律研究成果,n=0.74,2 a、5 a、10 a、20 a、50 a、100 a重现期的A值分别为40.60 mm/h、59.75 mm/h、73.16 mm/h、86.11 mm/h、102.82 mm/h、115.24 mm/h,代入式(1)可得到不同重现期的暴雨强度公式。

3.2 综合暴雨公式

由式(1)得出的单一重现期暴雨强度公式未包含暴雨重现期这一重要变量,不同重现期对应不同的暴雨公式,不便于实际应用。因此,借鉴《城市暴雨强度公式编制和设计暴雨雨型确定技术导则》[23]推荐的雨力计算公式,在公式(1)中引入暴雨重现期变量,雨力计算公式如下:

A=A1+ClgT

(2)

式中:A1、C为待求参数；T为重现期,a。

根据不同重现期T对应的A值,采用最小二乘法拟合得到参数A1=28.621、C=43.707,确定性系数0.999 1,拟合结果如图2所示。

图2 上海市暴雨重现期与雨力关系曲线

据此,不同的单一重现期暴雨公式可以统一用综合暴雨公式表达:

(3)

上式等号两侧同乘以历时t,可得到暴雨量计算公式:

H=(28.621+43.707lgT)t0.26

(4)

式中:H为暴雨量,mm。

3.3 暴雨公式精度检验

参考《城市暴雨强度公式编制和设计暴雨雨型确定技术导则》[23],当重现期为2～20 a时,在一般强度的地方,平均绝对均方根误差不宜大于0.05 mm/min；在较大强度的地方,平均相对均方根误差不宜大于5%。本文分别采用平均相对均方根误差和平均绝对均方根误差2项指标来检验推导的长历时综合暴雨公式的精度。经检验[23],上海市长历时暴雨公式的平均相对均方根误差为1.9%,平均绝对均方根误差为0.009 mm/min,均满足暴雨公式编制精度要求。

3.4 暴雨重现期公式

为了便于在实际工作中及时准确地确定实际场次暴雨的重现期,将式(4)推导为暴雨重现期计算公式:

(5)

由式(5)可以确定上海市1～24 h任意历时场次暴雨的重现期。

4 结 论

采用上海市徐家汇站1949—2013年实测降雨资料,确定暴雨强度衰减指数,结合雨力公式得到上海市长历时综合暴雨公式。主要结论如下:

(1) 不同重现期暴雨呈现的衰减规律一致,历时1～24 h的暴雨强度衰减指数均约为0.74。

(2) 根据暴雨强度衰减指数和1 h设计暴雨量,可以编制多个单一重现期的暴雨公式。

(3) 推求得到的长历时综合暴雨公式,实现了1～24 h任意历时、2～100 a任意重现期的设计暴雨计算公式的统一；由长历时综合暴雨公式推导的暴雨重现期计算公式,可以估算1～24 h任意历时场次暴雨的重现期。

(4) 编制的长历时综合暴雨公式满足暴雨公式编制的精度要求,平均相对和平均绝对均方根误差分别仅为1.9%和0.009 mm/min。

本文研制的长历时综合暴雨公式经实践检验是可行的,已纳入上海市《治涝标准:DB31/T1121—2018》,对其他城市编制长历时综合暴雨公式具有参考价值。