成都龙泉驿区暴雨变化特征

Characteristics of Rainstorm Variation in Longquanyi District of Chengdu

摘要: 利用龙泉驿区气象局1980~2019年近40年逐日降水观测资料与2009~2019年19个自动站的逐时降水观测资料，采用多项式插值和反距离权重插值等方法，分析了龙泉驿区暴雨的时空变化特征，得到以下结果：1) 近40年来暴雨量、暴雨日数、暴雨强度存在明显的年际变化，其旬、候变化呈现单峰型。2) 近40年来夏季暴雨量、暴雨日数在多年来暴雨中占比最高，暴雨主要集中在7、8月份。3) 在空间上，暴雨量、暴雨日数整体呈现南北高中西低，暴雨强度整体呈现东高西低。

Abstract: Based on the daily precipitation observation data of Longquanyi Meteorological Bureau from 1980 to 2019 and hourly precipitation observation data of 19 automatic stations from 2009 to 2019, the methods of polynomial interpolation and inverse distance weight interpolation are adopted, and the following results were obtained: 1) For 40 years, the precipitation of rainstorm, the days of rainstorm and the intensity of rainstorm have obvious inter-annual variation, and its ten-day and pentad variation presents a single peak pattern. 2) For 40 years, the precipitation of rainstorm and the days of rainstorm in summer accounted for the highest proportion of the rainstorm in many years, and the rainstorm mainly concentrated in July and August. 3) From the spatial point of view, the precipitation and days of rainstorm were higher in the north and south, lower in the middle and west, while the intensity of rainstorm was higher in the east and lower in the west.

1. 引言

根据联合国政府间的气候变化委员会在第五次科学评估的发现，在以往的几十年中，极端天气如强降水、热浪等呈增多趋势，未来极端天气出现的次数将会更多 [1]。中国处于亚欧大陆东部，暴雨天气在中国十分常见，其原因与夏季东亚季风有关，并且许多暴雨过程都存在季风水汽输送带，因夏季风从海上带来的水汽充足很容易形成洪涝灾害 [2]。中国的暴雨降水量自1951年开始逐年增加，降水的范围也在扩大 [3]，到了1990年代暴雨天气有所减少，在21世纪又开始增多，这种现象在南方地区表现明显 [4]。

在四川省，暴雨洪涝是发生频率最高、造成危害最重的自然灾害之一，其特点有范围广、强度大等，影响四川暴雨的系统有青藏高原低压、西南涡、西太平洋副热带高压、低空急流等 [5]。对于一些暴雨个例，来自孟加拉湾、阿拉伯海、西太平洋以及高原地区的水汽输送对其影响十分明显，由于水汽来源多样，处于盆地的四川暴雨过程机制复杂 [6] [7]。近几十年来，川渝地区暴雨主要发生在四川盆地，川西高原地区暴雨较少，四川盆地年暴雨发生频数大于2次 [8]。成都龙泉驿区地处成都市东部，被国务院命名为“中国水蜜桃之乡”，其境内包含国家经济技术、生态园、汽车基地等产业，在四川省内其发展潜力巨大 [9]。以全球气候变暖为背景，成都市龙泉驿境内暴雨灾害事件增多 [10]，据统计，暴雨洪灾造成的人口死亡、经济损失的数量在全球各种灾害排序中居首位，因此暴雨天气受各国普遍关注 [11]。在暴雨的变化特征的研究中，鲍名等人结合暴雨与洪涝的关系以及暴雨的气候背景，分析了暴雨日数的气候分布与暴雨日数及其距平的逐年变化 [12]。黄明策通过常规的统计分析方法，分析了广西暴雨量、暴雨日数的空间分布与时间分布 [13]。林爱兰等人通过新指标下的广东暴雨日和持续性暴雨的定义，分析了广东暴雨的变化特征及其环流形势 [14]。周长艳等人运用不同级别的降水事件频率和强度，揭示了四川暴雨以及其他降水的变化特征 [15]。在前人对暴雨的研究中，其研究的区域大多以省或国家为代表，并且对于暴雨的旬、候变化分析较少，而本文研究的龙泉驿区暴雨其空间范围小，且对暴雨量、日数、强度都做出了分析，在暴雨的时间演变部分细化到了旬、候变化，通过1980~2018年龙泉国家站的逐日降水量资料以及2009~2019年龙泉驿区内19个自动站的逐时降水量资料，分析暴雨量、暴雨日数、暴雨强度数据，以此对龙泉驿区暴雨的时空变化特征进行概括、总结。

2. 资料与方法

龙泉驿区位于四川省成都市平原东部，由图1可知：龙泉驿处于北纬30˚27'52~30˚43'23与东经104˚08'19~104˚27'09之间，境内东西长29.8 km、南北宽28.75 km，最高海拔为1050 m左右，最低海拔450 m左右。其总面积达556.98 km²，境内龙泉山位于东南部，呈东北–西南走向，平坝、丘陵、低山等地貌在龙泉都有表现，各地貌面积占比平坝最大、低山与丘陵次之。

2.1. 资料来源与地形概况

本文选用1980~2018年龙泉驿区气象局的逐日降水量观测资料，以及如图1所示的黄土、柏合、西河、万兴、大面、茶店、洪安、洛带、桃花故里、七里、红花、照壁、天公山、梨园村、同安、十陵、经开区、大佛村、成龙一号共19个自动站的2009~2019年逐时降水量观测资料，得到暴雨量、暴雨日数、暴雨强度的时间演变和空间演变的数据。

Figure 1. Terrain and automatic station distribution in Longquanyi District

图1. 龙泉驿区地形和自动站分布

2.2. 研究方法

2.2.1. 多项式插值

通过沈剑华 [16] 编著的数值计算基础得知，在许多实际的问题中我们可以用函数 $y=f\left(x\right)$ 来表示变量间的一些内在规律的数量关系。但是在科学研究或工程技术中，对于函数 $f\left(x\right)$ 有时只能通过实验或观测等手段得到它有限个不同点上的函数值，即表现为 $y=f\left(x_i\right)\left(i=0,1,2,\cdots ,n\right)$ 的情况，有时虽然函数有明确的表达式但是过于复杂，这些情况都不便于我们计算、使用。

为了研究上述函数的性质，我们希望做出一个能反映原函数特性，且便于计算的简单函数 $P\left(x\right)$ 来近似代替 $f\left(x\right)$，对于函数 $P\left(x\right)$ 要满足 $P\left(x_i\right)=y_i\left(i=0,1,2,\cdots ,n\right)$ 的条件，称这类问题为插值问题， $P\left(x\right)$ 为插值函数。设函数 $y=f\left(x\right)$ 在区间 $\left[a,b\right]$ 上有定义，并且它在 $n+1$ 个不同点 $a<=x_0<x_1<x_2<\cdots <x_n<=b$ 上的函数值 $y_0,y_1,\cdots ,y_n$ 都已知，若存在一个次数不超过 $n$ 次的多项式(其中 $a_i$ 为实数)：

$P_n\left(x\right)=a_0+a_{1}x+a_2{x}^2+\cdots +a_n{x}^n$ (1)

满足条件：

$P_n\left(x_i\right)=y_i,\left(i=0,1,2,\cdots ,n\right)$ (2)

称 $P_n\left(x\right)$ 为函数 $f\left(x\right)$ 的 $n$ 次代数插值多项式。本文通过1980~2018年近40年的龙泉驿区气象局逐日降水量观测资料，得到暴雨的时间演变数据，利用绘图工具绘制暴雨量、暴雨日数、暴雨强度的时间演变图形，并用多项式插值来叠加曲线，分析暴雨随时间变化的特征，最后得到龙泉驿区暴雨的时间演变概况。

2.2.2. 反距离权重插值

通过已知样品点采集的数据来推算未知点的数据称为空间插值，这种方法常用于气象领域，反距离权重法作为一种常用于气象要素的空间插值方法，其原理简单、使用方便，在气象等领域被广泛使用。反距离权重插值其插值结果好坏与搜索半径以及样品点数量等因素有关系，搜索半径越大差值结果越光滑，效果越好。若样品点数量少，则反距离权重插值结果随机性大，会出现较大差异 [17] [18]。反距离权重插值法是由Tobler定理为基础提出的空间插值方法。通过推算未测量点附近各点的测量值的加权平均来进行插值是其原理，在空间上越相近的事物其相似性越大，在附近已测点取得的权重越大，对于反距离权重插值法，其表达式为：

$Z_o=\left({\displaystyle {\sum }_{i=1}^n\frac{z_i}{d_i^r}}\right)/\left({\displaystyle {\sum }_{i=1}^n\frac{1}{d_i^r}}\right)$ (3)

式子中 $z_o$ 为 $O$ 点估算值， $z_i$ 控制点 $i$ 的值， $d_i$ 为控制点 $i$ 与点 $O$ 的距离， $n$ 为估算时所用控制点数量， $r$ 为指定的幂次数 [19]。

本文通过2009~2019年近10年的龙泉驿区内19个自动站点逐时降水量观测资料，得到暴雨的空间演变数据，利用GIS等软件的反距离权重插值，绘制出暴雨量、暴雨日数、暴雨强度的空间分布图，通过分析不同时间下暴雨的空间变化特征，最后得到龙泉驿区暴雨的空间演变概况。

3. 龙泉驿区暴雨时间演变特征

3.1. 暴雨量时间演变特征

根据图2(a)~(c)龙泉驿区春、夏、秋各季暴雨量的时间演变图来看：春季暴雨量最高值为1997年的91.0 mm，最低值为1983年的50.2 mm，平均春季暴雨量为5.2 mm。39年间只有1983、1986、1997年这3年存在春季暴雨，因此春季暴雨的变化趋势参考性不大；夏季暴雨量最高值为1984年的364.0 mm，最低值为1986年的51.8 mm，平均夏季暴雨量为113.4 mm。夏季暴雨量在1980~1989年左右呈上升趋势，在1990~2002年左右呈下降趋势，2003~2018年夏季暴雨量波动较平稳；秋季暴雨量最高值为2001年的116.6 mm，最低值为2012年的50.9 mm，平均秋季暴雨量为24.8 mm。39年间秋季暴雨量呈先增再减趋势；冬季无暴雨，故不讨论。从图2(d)龙泉驿区年暴雨量的时间演变曲线来看：龙泉驿区暴雨量最高值为1984年的364.0 mm，最低值为1993年的59.3 mm，平均暴雨量为143.4 mm。差值曲线上，暴雨量在1985~1990年左右呈上升趋势在1991~2002年左右呈下降趋势，2005~2015年暴雨量升降趋势较为平稳。从图2(e)~(f)暴雨量的旬、候变化来看：8月中旬的暴雨量最高，8月第4候的暴雨量最高，暴雨量的旬、候变化都呈现单峰型，且只在4~10月有暴雨量其他月份无暴雨量。

综合来看，暴雨量存在明显的年际变化，且暴雨量的旬、候变化都呈现单峰型，夏季暴雨量约占总暴雨量的79%，而秋季暴雨约占17%，春季暴雨约占4%，冬季无暴雨出现。近40年间1~12月中7、8月份暴雨量占比最高。

Figure 2. Temporal variation of rainstorm precipitation in Longquanyi District. (a)~(c): Rainstorm precipitation in spring, summer and autumn (Unit: mm), (d): Annual rainstorm precipitation (Unit: mm), (e): Ten-day rainstorm precipitation (Unit: mm), (f): Pentad rainstorm precipitation (Unit: mm)

图2. 龙泉驿区暴雨量的时间演变。(a)~(c)：春、夏、秋季暴雨量(单位：mm)，(d)：年暴雨量(单位：mm)，(e)：旬暴雨量(单位：mm)，(f)：候暴雨量(单位：mm)

3.2. 暴雨日数时间演变特征

从图3(a)~(c)龙泉驿区春、夏、秋各季暴雨日数的时间演变曲线来看：春季暴雨日数最高、最低值都为1 d，平均春季暴雨日数为0.1 d。由于春季暴雨日数较少，不讨论其变化趋势；夏季暴雨日数最高值为1984年的5 d，最低值为1 d，平均夏季暴雨日数为1.7 d。在1980~1989年左右夏季暴雨日数呈上升趋势，1990~2003年左右呈下降趋势，2004年~2018年暴雨日数波动不明显；秋季暴雨日数最高值为2 d，最低值为1 d，平均秋季暴雨日数为0.4 d。插值曲线显示，39年来秋季暴雨日数变化波动不大。冬季无暴雨，故不讨论。从图3(d)龙泉驿区年暴雨日数的时间演变曲线来看：暴雨日数最高值为1984、2018年的5 d，最低值为1 d，平均暴雨日数为2.2 d。在1985~1989年左右暴雨日数呈上升趋势，1990~2000年左右呈下降趋势，2001~2015年左右暴雨日数波动平稳。从图3(e)、图3(f)暴雨日数的旬、候变化来看：8月中旬的暴雨日数最高，8月第4候的暴雨日数最高，暴雨日数的旬、候变化都呈现单峰型，且只在4~10月有暴雨日数其他月份无暴雨日数。

综合来看，暴雨日数存在明显的年际变化，且暴雨日数的旬、候变化都呈现单峰型，夏季暴雨日数约占总暴雨降水日数的77%，秋季暴雨日数约占18%，春季暴雨日数约占5%，而冬季无暴雨出现。近40年间1~12月份中只有4~10月份存在暴雨日数，其中7、8月份暴雨日数占比最高。

Figure 3. Temporal variation of rainstorm days in Longquanyi District. (a)~(c): Number of rainstorm days in spring, summer and autumn (Unit: d), (e): Number of annual rainstorm days (Unit: d), (e): Number of ten-day rainstorm days (Unit: d), (f): Number of pentad rainstorm days (Unit: d)

图3. 龙泉驿区暴雨日数的时间演变。(a)~(c)：春、夏、秋季暴雨日数(单位：d)，(d)：年暴雨日数(单位：d)，(e)：旬暴雨日数(单位：d)，(f)：候暴雨日数(单位：d)

3.3. 暴雨强度时间演变特征

从图4(a)~(c)龙泉驿区季暴雨强度的时间演变曲线来看：春季暴雨强度最高值为1997年的91.0 mm/d，最低值为1983年的50.2 mm/d，平均春季暴雨强度为68.2 mm/d。由于39年间春季暴雨数量少，因此不分析其降水强度的升降趋势；夏季暴雨强度最高值为1987年的91.4 mm/d，最低值为1986年的51.8 mm/d，平均夏季暴雨强度为65.2 mm/d。差值曲线显示，在1983、2002年左右夏季暴雨强度存在波动谷值，1991、2011年左右夏季强度存在波动峰值；秋季暴雨强度最高值为1995年的96.9 mm/d，最低值为2015年的50.0 mm/d，平均秋季暴雨降水强度65.6 mm/d。差值曲线显示，39年来秋季暴雨强度呈多次波动状态。冬季无暴雨，故不讨论。从图4(d)龙泉驿区年暴雨强度的时间演变曲线来看：近40年来，龙泉驿区年暴雨强度最高值为1987年的91.4 mm/d，最低值为1994年的53.8 mm/d，平均暴雨强度为57.8 mm/d。差值曲线显示，1980~1994年左右暴雨强度呈下降趋势，1985~1993年左右暴雨强度呈上升趋势，1994~2018年左右暴雨强度主要呈下降趋势。从图4(e)、图4(f)暴雨强度的旬、候变化来看：4月下旬暴雨降水强度最高，5月下旬暴雨降水强度最低，6月第2候的暴雨降水强度最高，5月第5候的暴雨降水强度最低。

综合来看，暴雨强度存在明显的年际变化，且暴雨强度的旬、候变化都呈现单峰型，龙泉驿区暴雨主要集中在7、8月份，近40年间平均暴雨强度为66.3 mm/d。春、夏、秋季平均暴雨强度都超过了60.0 mm/d。

Figure 4. Temporal variation of rainstorm intensity in Longquanyi District. (a)~(c): Rainstorm intensity in spring, summer and autumn (Unit: mm/d), (d): Annual rainstorm intensity (Unit: mm/d), (e): Ten-day rainstorm intensity (Unit: mm/d), (f): Pentad rainstorm intensity (Unit: mm/d)

图4. 龙泉驿区暴雨强度的时间演变。(a)~(c)：春、夏、秋季暴雨强度(单位：mm/d)，(d)：年暴雨强度(单位：mm/d)，(e)：旬暴雨强度(单位：mm/d)，(f)：候暴雨强度(单位：mm/d)

4. 龙泉驿区暴雨空间演变特征

4.1. 暴雨量空间演变特征

从图5(a)~(c)龙泉驿区春、夏、秋季总暴雨量的空间分布来看：春季暴雨量呈现三点分布，这三点为西北部的西河和龙泉山北部的梨园村南部的红花，这三处附近暴雨量在49.61~74.39 mm之间，其余地区降水量基本在0.00~12.40 mm间，无暴雨；夏季暴雨量呈现南北高中西低，在龙泉北部黄土、洛带和龙泉山南部的红花附近暴雨量最高达到1393.38 mm，南北大部分地区暴雨量在830.57 mm以上，龙泉西部十陵暴雨量最低，在267.77~455.37 mm间，龙泉中西部暴雨量基本在830.58 mm以下；秋季暴雨量呈现东北高西南低，龙泉北部洪安、洛带附近暴雨量最高达到302.07 mm，龙泉东北半区暴雨量基本在151.09 mm以上，西南部最低暴雨量在0.10~50.43 mm间，龙泉西南半区暴雨量基本在151.10 mm以下。冬季无暴雨，故不讨论。从图5(d)龙泉驿区总暴雨量的空间分布来看：龙泉驿区总暴雨量呈现南北高中西低，近10年间龙泉驿区北部黄土、洛带、洪安和龙泉山南部红花的暴雨量最高，在1371.31~1563.94 mm之间，中西部地区暴雨量偏低，西部十陵暴雨两最低达到408.10 mm，南北大部分地区暴雨量在986.02 mm以上，中西大部分地区暴雨量在986.03 mm以下。

综合来看，暴雨量在龙泉驿区的分布特点为北部、东南部多，西部和中部少。在北部黄土、洪安、洛带以及东南部的七里、红花、茶店附近暴雨量较高，西部十陵中部成龙一号附近暴雨量较低。

(a) (b) (c) (d)

Figure 5. Spatial distribution of rainstorm precipitation in Longquanyi District. (a)~(c): Total rainstorm precipitation in spring, summer and autumn (Unit: mm); (d): Total rainstorm precipitation (Unit: mm)

图5. 龙泉驿区暴雨量的空间分布。(a)~(c)：春、夏、秋季总暴雨量(单位：mm)，(d)：总暴雨量(单位：mm)

4.2. 暴雨日数空间演变特征

从图6(a)~(c)龙泉驿区春、夏、秋季总暴雨日数的空间分布来看：春季暴雨日数主要分布在龙泉西北部的西河以及龙泉山北部的梨园村和南部的七里、红花，暴雨日数最高达到1 d，其余地区基本上无暴雨存在；夏季暴雨日数呈现南北高中西低，在龙泉北部的黄土和龙泉山南部的红花附近暴雨日数最高达到22 d，南北大部分地区暴雨日数在13 d以上，中西大部分地区暴雨日数在14 d以下，龙泉西部的十陵暴雨日数最低达到4 d；秋季暴雨日数呈现东北高西南低，在龙泉北部的洛带附近暴雨日数最高达到5 d，东北半区的大部分地区暴雨日数在3 d以上，西南半区暴雨日数基本在3 d以下，西南部的柏合等地无暴暴雨；冬季无暴雨，故不讨论。从图6(d)龙泉驿区总暴雨日数的空间分布来看：龙泉驿区暴雨日数整体呈现南北高中西低，暴雨日数最高的地区为龙泉北部的黄土、洛带附近，其暴雨日数在23~25 d之间，暴雨日数最低的地区为龙泉西部的十陵附近，在6~9 d之间，南北大部分地区暴雨日数在16 d以上，中心部地区暴雨日数基本在17 d以下。

综合来看，暴雨日数在龙泉驿区的分布特点为北部、东南部多，西部和中部少。在北部的黄土、洛带附近以及龙泉山南部的红花、七里、茶店附近较高，在西部的十陵和中部成龙一号附近较低。

(a) (b) (c) (d)

Figure 6. Spatial distribution of rainstorm days in Longquanyi District. (a)~(c): Total days of rainstorm in spring, summer and autumn (Unit: d), (d): Total days of rainstorm (Unit: d)

图6. 龙泉驿区暴雨日数的空间分布。(a)~(c)：春、夏、秋季总暴雨日数(单位：d)，(d)：总暴雨日数(单位：d)

4.3. 暴雨强度空间演变特征

从图7(a)~(c)龙泉驿区春、夏、秋季暴雨强度的空间分布来看：春季暴雨强度在西北部的西河以及龙泉山北部的梨园村、南部的红花附近最高，在62.00~74.39 mm/d之间，其余地区基本无暴雨；夏季暴雨强度整体呈现东高西低，龙泉山南部的桃花故里附近暴雨降水强度最高，在71.18~73.45 mm/d之间，东部大部分地区暴雨强度在66.62 mm/d以上，西北部西河的暴雨强度最低，在59.80~62.08 mm/d之间，西部大部分地区暴雨强度在66.63 mm/d以下；秋季暴雨强度整体呈现东北高西南低，在东北大部分地区暴雨强度在63.97~76.75 mm/d之间，西南部柏合附近降水强度最低在0.05~12.84 mm/d之间，其余地区暴雨强度基本在51.19~63.96 mm/d间。冬季无暴雨，故不讨论。从图7(d)龙泉驿区暴雨强度的空间分布来看：龙泉驿区暴雨强度整体呈现东高西低，在北部万兴和龙泉山南部大佛村、桃花故里、茶店附近暴雨强度最高，在68.08~69.68 mm/d之间，龙泉山一带暴雨强度在66.48 mm/d以上，龙泉西北西河附近暴雨强度最低，在60.05~61.66 mm/d之间，西部大部分地区暴雨强度在64.87 mm/d以下。

综合来看，5~9月暴雨强度在龙泉山北部的梨园村和南部的茶店、桃花故里附近最高。4月和10月暴雨强度在北部的西河、黄土、万兴附近最高，其余地区基本无暴雨。

(a) (b) (c) (d)

Figure 7. Spatial distribution of rainstorm intensity in Longquanyi District. (a)~(c): Rainstorm intensity in spring, summer and autumn (Unit: mm/d), (d): Rainstorm intensity (Unit: mm/d)

图7. 龙泉驿区暴雨强度的空间分布。(a)~(c)：春、夏、秋季暴雨强度(单位：mm/d)，(d)：暴雨强度(单位：mm/d)

5. 结论

综合上述对龙泉驿区暴雨的时间和空间演变的分析，得到以下结论：

1) 近40年间，暴雨量、日数、强度存在明显的年际变化，其旬、候变化表现为单峰型，在四季中，夏季暴雨量、日数最高，春秋两季较低。近40年间的各月中，4~10月存在暴雨，且暴雨主要集中在7、8月份。

2) 近10年间龙泉驿区暴雨量、日数呈现南北高中西低，强度呈现东高西低。龙泉驿区北部和龙泉山南部暴雨量、日数较高。北部和龙泉山南部附近暴雨强度较高，西北部的暴雨强度较低。

本文主要从暴雨量、日数、强度三个方面来分析龙泉驿区的暴雨时空变化特征，探究了龙泉驿区暴雨的时空变化规律。在暴雨的形成、发展等方面还未进行探讨，对于这些方面还应采用更加多样的分析方法进行详细研究，这是我下一步的研究目标。