貴州大龍經濟開發區極端氣象參數研究

李 霄，帥士章 ，張東海，朱 軍

(貴州省氣候中心，貴州 貴陽 550002)

0 引言

在全球氣候變暖的背景下，各類極端天氣氣候事件不斷增多，極端強降水事件、極端高溫事件更加嚴重和頻繁。研究表明全國范圍內平均降水、降水強度、極端強降水和連續性強降水呈增強趨勢[1-3]。張嬌艷等[4]在研究全球升溫1.5 ℃和2.0 ℃情景下貴州省極端降水的變化特征時，發現貴州在RCP2.6 和RCP4.5 情景下各極端降水指數總體上均呈現增加趨勢。袁文德等[5]、戴聲佩等[6]發現西南地區、華南地區極端高溫事件和極端低溫事件分別呈上升和下降趨勢。朱大運等[7]同樣發現貴州極端氣溫暖系列指數呈上升趨勢、冷系列指數呈下降趨勢。與極端強降水事件、極端高溫事件呈上升趨勢不同的是，有研究表明我國平均風速、最大風速和極大風速整體以減少趨勢為主，但不同區域風速的波動性存在差異，如西南區東部平均風速變化趨勢較不明顯[8-10]，江瀅等[11]研究指出21 世紀中國區域近地層年平均風速呈減小的趨勢。隨著我國城鎮化建設，預計未來城市極端高溫、極端降水事件將更為頻發，氣候變化風險進一步加大[12]。因此，在極端天氣氣候事件頻發的今天，氣候可行性論證是科學應對氣候變化、做好防災減災工作的具體行動和措施，是在預設的安全系數下為工程項目算好經濟帳的基礎性工作，氣候可行性論證工作正在重大工程及區域經濟發展中發揮著重要的作用[13-14]。

貴州大龍經濟開發區位于貴州省銅仁市玉屏縣，是1999年經貴州省人民政府批準設立、2006年經國家發改委審核通過的省級經濟開發區，2014年成為國家級循環化改造示范試點園區。本文以貴州大龍經濟開發區為例，討論在氣候變化背景下，影響開發區安全生產的主要極端氣象參數重現期。選取距離開發區最近、地形地貌相似的玉屏、新晃國家氣象觀測站資料，對影響開發區工程項目建設和安全生產較大的極端氣溫、極端強降水和極值風速進行分析，旨在為開發區企業科學認識當地氣象災害現狀、防范氣象災害提供科技支撐。

1 資料和方法

1.1 資料來源

本文使用的資料為玉屏、新晃氣象觀測站1959—2021年歷年極端最高、最低氣溫資料，歷年1 d、3 d、5 d、10 d最大降水量資料，1980—2021年歷年最大風速、極大風速資料。資料來自全國綜合氣象信息共享平臺(CIMISS)及大數據云平臺“天擎”。

1.2 研究方法

在進行極端氣象參數計算時，選用國內外應用比較廣泛的皮爾遜Ⅲ(Pearson-Ⅲ)和耿貝爾(Gumbel)分布函數，對玉屏、新晃氣象觀測站極端氣溫、降水量和風速數據進行擬合，推算不同重現期的設計值。

在處理風速資料時，采用風速廓線的冪指數公式對觀測高度進行修正，同時對10 min平均年最大風速缺測數據采用定時觀測2 min的年最大風速進行插補。

1.2.1 皮爾遜Ⅲ(Pearson-Ⅲ)分布函數 Pearson-Ⅲ型分布的概率密度函數f(x)和累積分布函數F(x)分別見下式：

(1)

(2)

式中F(x)為α的伽瑪函數，α為形狀參數，β為比例參數(尺度參數)，ɑ0為位置參數。皮爾遜Ⅲ型分布其他常用參數有x均值、Cv離(變)差系數、Cs偏態系數、σ標準差，參數間換算關系：

(3)

(4)

(5)

(6)

(7)

(8)

σ=x·Cv

(9)

已知累積概率P(xp)，求極值xp時，可根據分布函數進行變量轉換并按伽瑪函數積分求得，也可查皮爾遜Ⅲ型分布的離均系數表求得。

1.2.2 耿貝爾(Gumbel)分布函數 Gumbel分布的概率密度函數和累積分布函數分別見下式：

f(x)=ɑe-ɑ(x-b)-e-ɑ(x-b)(-∞0,-∞

(10)

F(x)=e-e-α(x-b)

(11)

式中ɑ為尺度參數，b為位置參數。可設轉換變量y=ɑ(x-b)以方便計算。已知累積概率P(xp)，按下式求極端事件的極值。

(12)

1.2.3 測風數據修正 玉屏、新晃站測風觀測儀器經歷了幾次變換，需要將不同觀測高度的觀測數據修正到標準高度，即地面上方10 m處高度。本文采用風速廓線的冪指數公式來做觀測高度修正[12]。即：

(13)

式中：h為參照高度，取10 m；z為測量高度，單位為 m；Vh、Vz為相應的風速；α為冪指數(風速切變指數)，取決于表面粗糙度等因素，開發區位于房屋比較稀疏的鄉鎮，屬B類地區，冪指數(風速切變指數)取0.15。

玉屏站1998—2004年10 min平均年最大風速缺測，用定時觀測2 min年最大風速進行插補，定時觀測2 min年最大風速與10 min平均年最大風速換算關系[13-14]見下式：

V10 min=1.14V2 min(定時)

(14)

2 極端氣象參數計算

2.1 極端氣溫

統計玉屏、新晃氣象觀測站1959—2021年歷年極端最高、最低氣溫(圖1)。由圖1看出，極端最高氣溫出現在新晃站1971年的40.9 ℃(圖1a)，極端最低氣溫出現在玉屏站1977年的-10.7 ℃(圖1 b)。玉屏、新晃2站歷年極端最高、最低氣溫變化基本相似，總體上，新晃極端最高氣溫比玉屏偏高，但20世紀80年代中期及之后，新晃極端最低氣溫比玉屏偏低。

圖1 玉屏、新晃歷年極端最高氣溫(a)和極端最低氣溫(b)變化圖Fig.1 Time series of extreme maximum temperature (a) and minimum temperature (b) in Yuping and Xinhuang

利用玉屏、新晃氣象觀測站的最高氣溫資料系列，采用Pearson-Ⅲ和Gumbel分布進行擬合，計算出2個氣象觀測站最高氣溫百年一遇值，結果見表1，分布曲線擬合見圖2～3。由圖可以看出，Pearson-Ⅲ 和Gumbel分布對玉屏和新晃最高氣溫擬合效果較好，高重現期或低重現期2種分布函數擬合結果均大于2站實測值。由表1看出，2站最高氣溫百年一遇值Pearson-Ⅲ 和Gumbel分布計算結果均大于各站實測值，擬合結果新晃站Gumbel分布最高為42.7 ℃，基于企業工程項目建設和安全生產考慮，開發區區域最高氣溫百年一遇設計值采用42.7 ℃。

表1 玉屏、新晃極端最高氣溫百年一遇設計值(單位：℃)Tab.1 Design value of extreme maximum temperature at 100-year return period in Yuping and Xinhuang (Unit: ℃)

圖2 玉屏最高氣溫Pearson-Ⅲ(a)和Gumbel(b)分布曲線擬合Fig.2 Pearson-Ⅲ (a) and Gumbel (b) distribution curve fitting of the maximum temperature in Yuping

圖3 新晃最高氣溫Pearson-Ⅲ(a)和Gumbel(b)分布曲線擬合Fig.3 Pearson-Ⅲ (a) and Gumbel (b) distribution curve fitting of the maximum temperature in Xinhuang

利用玉屏、新晃氣象觀測站的最低氣溫資料系列，采用Pearson-Ⅲ 和Gumbel分布進行擬合，計算出2個氣象觀測站最低氣溫百年一遇值，結果見表2，分布曲線擬合見圖4。由圖4可以看出，Pearson-Ⅲ 和Gumbel分布對玉屏和新晃最高氣溫擬合效果較好，尤其是低重現期部分，Gumbel分布函數對2站在高重現期部分的擬合結果大于2站實測值。由表2看出，玉屏、新晃氣象觀測站最低氣溫百年一遇值Pearson-Ⅲ 計算結果均低于實測值，但Gumbel分布計算結果均高于實測值；2個氣象觀測站最低氣溫百年一遇值Pearson-Ⅲ 和Gumbel分布計算結果均以玉屏站為最低。因此，基于企業工程項目建設和安全生產考慮，開發區區域最低氣溫百年一遇設計值采用-10.9 ℃。

表2 玉屏、新晃最低氣溫百年一遇設計值(單位：℃)Tab.2 Design value of Extreme minimum temperature at 100-year return period in Yuping and Xinhuang (Unit: ℃)

圖4 玉屏(a、b)新晃(c、d)最低氣溫Pearson-Ⅲ(左)和Gumbel(右)分布曲線擬合Fig.4 Pearson-Ⅲ (left) and Gumbel (right) distribution curve fitting of the minimum temperature in Yuping(a、b) Xinhuang(c、d)

2.2 極端強降水

統計玉屏、新晃氣象觀測站1959—2021年歷年1 d、3 d、5 d、10 d最大降水量(圖5)，由圖5看出，1 d、10 d極端最大降水量出現在1960年的玉屏，分別為226.2 mm、316.7 mm(圖5a)，3 d、5 d極端最大降水量出現在2007年的新晃，分別為245.3 mm、279.9 mm(圖5b)。

利用玉屏、新晃氣象觀測站的1 d、3 d、5 d、10 d最大降水量資料系列，采用Pearson-Ⅲ和Gumbel分布進行擬合，計算2個氣象觀測站不同日數最大降水量百年一遇值，結果見表3，分布曲線擬合見圖6(部分分布曲線擬合圖略)。由圖可以看出，除1 d最大降水量分布擬合效果不理想外，其余不同日數最大降水量分布函數對2站的擬合效果均較好，高重現期部分擬合值和實測值相當。由表3可知，2站不同日數最大降水量Pearson-Ⅲ 和Gumbel分布百年一遇擬合計算結果均大于實測值，基于企業工程項目建設和安全生產考慮，開發區區域1 d、3 d、5 d、10 d最大降水量百年一遇設計值分別采用248.0 mm、261.8 mm、300.8 mm、327.7 mm。

圖5 玉屏(a)、新晃(b)歷年不同日數最大降水量變化圖Fig.5 Time series of Maximum precipitation in different durations in Yuping (a) and Xinhuang (b)

表3 玉屏、新晃不同日數最大降水量百年一遇設計值(單位：mm)Tab.3 Design value of maximum precipitation at 100-year return period at different durations in Yuping and Xinhuang(unit: mm)

圖6 玉屏1 d(a)、10 d(b)和新晃3 d(c)、5 d(d)最大降水量Pearson-Ⅲ 分布曲線擬合Fig.6 Pearson-Ⅲ distribution curve fitting of the maximum precipitation of 1-day (a)、10-day (b) in Yuping and 3-day(c)、5-day(d) in Xinhuang

2.3 極端風速

2.3.1 測風資料 《地面氣象觀測規范》(GB/T 35227-2017)對最大風速和極大風速作如下規定：最大風速“從10 min滑動平均風速值中挑取，并記錄相應的風向和時間”。極大風速“從3 s滑動平均風速值中挑取，并記錄相應的風向和時間”。

據統計(表略)，新晃觀測站極大風速為39.6 m·s-1，出現在2002年4月1日；最大風速為25.3 m·s-1，出現在2015年6月8日；均為建國以來開發區區域內出現的極大和最大風速正式記錄，是強對流天氣過程導致(風雹)。

由于玉屏氣象觀測站1998—2004年風速缺測，故采用2 min年最大風速進行插補，統計得出玉屏、新晃站10 min平均最大風速資料(表4)。由表4看出，新晃平均最大風速、平均極大風速均大于玉屏，表明新晃的風速比玉屏的大。經測風高度修正后得出玉屏、新晃站年最大風速變化情況(圖7)。由圖看出，玉屏、新晃10 min年平均最大風速變化趨勢基本一致，新晃最大值25.3 m·s-1出現在2015年，玉屏最大值23.1 m·s-1出現在2014年，新晃最小值8.1 m·s-1出現在2018年，玉屏最小值4.3 m·s-1出現在1987年。由圖7還可以看出，2站10 min年平均最大風速最大值基本出現在20世紀90年代初期、21世紀初期和10年代中前期，而最小值基本出現在20世紀80年代中后期、90年代中后期和21世紀 10年代初期。總體來說，近40 a新晃最大風速、極大風速均比玉屏偏大。

表4 玉屏、新晃極值風速資料統計表Tab.4 Information about extreme wind speed in Yuping and Xinhuang

圖7 玉屏、新晃歷年10 min平均最大風速變化圖Fig.7 Time series of average maximum wind speed of 10 minutes in Yuping and Xinhuang

2.3.2 最大風速統計分析 利用新晃、玉屏氣象觀測站30 a以上10 min平均年最大風速資料序列，采用Pearson-Ⅲ和Gumbel分布計算頻率為3.33%、2.0%、1%、0.5%和0.1%的最大風速，結果見表5，分布曲線擬合見圖8。玉屏、新晃站實測最大風速分別為23.1 m·s-1和25.3 m·s-1。由圖8可以看出，Pearson-Ⅲ和Gumbel分布對玉屏和新晃最大風速擬合效果較好，2種分布函數百年一遇擬合結果均大于2站實測值。

從表5可知，各設計頻率(重現期)風速，新晃站均大于玉屏站，新晃站Pearson-Ⅲ擬合值大于Gumbel擬合值。基于企業工程項目建設和安全生產考慮，推薦開發區采用新晃站Pearson-Ⅲ計算成果，該站10 m高度百年一遇最大風速擬合值為26.7 m·s-1，大于實測值。按照《建筑結構荷載規范》(GB50009-2012)，當R=50時，風壓=0.3 kN·m-2，新晃累年氣壓974.9 hPa，氣溫17.5 ℃，水氣壓16.5 hPa，可計算得到五十年一遇重現期風速為22.7 m·s-1，小于推薦值。R=100時，風壓=0.35 kN·m-2，新晃累年氣壓974.9 hPa，氣溫17.5 ℃，水氣壓16.5 hPa，可計算得到百年一遇重現期風速為24.5 m·s-1，小于推薦值。

2.3.3 最大風速與極大風速關系分析 新晃站2001年開始有年極大風速觀測資料，為使分析結果更具代表性，收集該站10 min平均最大風速出現時間，挑選同處1次大風過程且同小時出現的極大風速資料，進行高度訂正后，最大風速與極大風速的歷年變化見圖9。由圖可以看出，新晃歷年10 min平均最大風速與極大風速變化趨勢基本一致，最大風速最大值為2015年的25.3 m·s-1，最小值為2018年的8.1 m·s-1；極大風速最大值為2002年的38.9 m·s-1，最小值為2001年的12.3 m·s-1。

表5 玉屏、新晃10 min平均最大風速各頻率擬合值(單位：m·s-1)Tab.5 Fitting values of frequency of average maximum wind speed at 10 minute in Yuping and Xinhuang(unit: m·s-1)

有極大風速觀測的近21 a，新晃瞬時極大風速與10 min平均最大風速同時次的平均比值為1.6。

2.3.4 各高度各頻率最大風速與極大風速計算 根據新晃站采用Pearson-Ⅲ分布計算的頻率為3.33%、2.0%、1%、0.5%和0.1%最大風速結果，計算各離地高度不同重現期最大風速；同時利用極大風速與最大風速比值計算得到的不同頻率極大風速結果，計算各離地高度不同重現期極大風速。各高度的推算采用風廓線冪次律公式，冪指數值取0.15，結果見表6。由表6可知，開發區離地10 m高度百年一遇最大風速為26.7 m·s-1，百年一遇極大風速為42.7 m·s-1。

圖8 玉屏(a、b)、新晃(c、d)最大風速Pearson-Ⅲ(左)和Gumbel(右)分布曲線擬合Fig.8 Pearson-Ⅲ (left) and Gumbel distribution (right) curve fitting of the maximum wind speed in Yuping(a、b)and Xinhuang(c、d)

3 結論與討論

通過對大龍經濟開發區距離最近、地形地貌相似的玉屏、新晃國家氣象觀測站極端氣溫、極端強降水和極值風進行分析，采用Pearson-Ⅲ和Gumbel分布擬合，計算百年一遇設計值，得出如下結論：

①大龍經濟開發區區域極端最高氣溫為40.9 ℃，極端最低氣溫為-10.7 ℃；1 d、3 d、5 d、10 d極端最大降水量分別為226.2 mm、245.3 mm、279.9 mm 、316.7 mm；最大風速為25.3 m·s-1，極大風速為39.6 m·s-1。

表6 開發區最大風速、極大風速各高度各頻率風速Tab.6 Maximum wind speed and maximum instantaneous wind speed of different frequency at different altitudes in developing area

②大龍經濟開發區區域最高氣溫百年一遇設計值為42.7 ℃，最低氣溫百年一遇設計值為-10.9 ℃；1 d、3 d、5 d、10 d最大降水量百年一遇設計值分別為248.0 mm、261.8 mm、300.8 mm、327.7 mm；離地10 m高度百年一遇最大風速為26.7 m·s-1，百年一遇極大風速為42.7 m·s-1。

本文針對大龍經濟開發區的極值氣溫、極端強降水和極值風不同重現期進行了分析，可為開發區工程項目設計和安全生產提供科學的氣象設計參數。但貴州地處低緯高原，境內山高峪深、下墊面崎嶇不平，造成氣溫、降水空間分布復雜多變，風速局地差異明顯，只選取玉屏、新晃2個站點作為大龍經濟開發區極端氣象參數研究還存在一定的局限性，有待進一步研究。