1951—2016年臨沂市降水變化特征分析

丁文超 趙欣 王凱

摘要 利用1951—2016年臨沂市年均降水量數據，采用數理統(tǒng)計和Morlet小波分析方法，研究了臨沂市降水變化多時間尺度的周期性變化規(guī)律，并對臨沂市未來降水量變化趨勢進行了預測。結果表明，1951—2016年臨沂市年均降水量最低值為529.5 mm，最高值為1 415.6 mm，平均值為864.3 mm，中位數為846.5 mm，符合正態(tài)分布；臨沂市年均降水量在15年、28年、40年左右周期震蕩劇烈，有明顯的規(guī)律；降水偏少的趨勢將結束，未來幾年降水量將逐漸偏多；40年時間尺度的小波系數的模方值和方差系數最大，能量最強、周期最顯著，為降水量變化的第1主周期。

關鍵詞 降水量；數理統(tǒng)計；小波分析；變化特征；山東臨沂；1951—2016年

中圖分類號 S161.6 文獻標識碼 A 文章編號 1007-5739（2018）10-0205-02

Analysis on Precipitation Change Characteristics in Linyi City from 1951 to 2016

DING Wen-chao 1 ZHAO Xin 1 WANG Kai 2 XU Zhi-hui 1 ZHANG Ying 1 WU Di 1 WANG Chang-wen 1 *

（1 College of Resources and Environment，Linyi University，Linyi Shandong 276000； 2 Linyi Meteorological Bureau）

Abstract Based on the annual average precipitation data in Linyi City from 1951 to 2016，the mathematical statistics and Morlet wavelet analysis method were used to study the periodic change pattern of precipitation with multiple time scales，and forecast the precipitation trend in the future of Linyi City.The results of average annual precipitation from 1951 to 2016 in Linyi City was as follows：the lowest value was 529.5 mm，the maximum value was 1 415.6 mm，the average value was 864.3 mm，the median value was 846.5 mm，which conformed to the normal distribution. The average annual precipitation in Linyi City vibrated in the period of 15 years，28 years，40 years，having obvious regularity. The trend of less precipitation will be end and there will be more precipitation in the next few years. Modulus of wavelet coefficient and variance coefficient of 40 years scale were maximum，which had strongest energy and the most obvious cycle，and became the first main precipitation period.

Key words precipitation；mathematical statistics；wavelet analysis；change characteristic；Linyi Shandong；1951- 2016

Morlet提出的小波分析具有良好的時域、頻域多分辨分析能力，逐漸地被應用于各領域，且取得了較好的效果。近年來，在降水量的時間序列中引入小波分析，在對研究降水量的規(guī)律性變化方面取得了卓有成效的進展[1-2]。劉忠陽等[3]利用Morlet小波分析方法對鄭州市近54年降水的小波變化特征進了分析；邱海軍等[4]利用小波分析對近60年來中國洪澇災情變化趨勢持續(xù)性和周期性進行了研究；楊 義等[5]以甘肅省白銀市近50年氣候變化為例，采用Mann-Kendall法和小波分析研究了中小尺度多年氣候變化特征。

臨沂市屬溫帶季風區(qū)大陸性氣候，降水主要集中在夏季。農業(yè)生產仍以大氣降水為主，尤其北部沂蒙山區(qū)，對降水的蓄積成為農業(yè)灌溉的主要方式。此外，全球氣候變暖對臨沂區(qū)域氣象條件、水資源、農業(yè)生產和生態(tài)環(huán)境產生深刻的影響。因此，掌握降水量的變化規(guī)律，對合理支配臨沂市自然降水資源十分必要。本文對1951—2016年臨沂市年均降水量進行統(tǒng)計學研究，利用Morlet連續(xù)小波分析年均降水量變化的特點，揭示不同層次上的降水結構、異常變化規(guī)律及突變特征，并預測未來降水量的變化趨勢，服務農業(yè)生產。

1 資料與方法

1.1 資料來源

本文選取1951—2016年臨沂市年平均降水量資料，數據資料由臨沂市氣象局提供。

1.2 分析方法

對近65年臨沂市年平均降水量進行統(tǒng)計分析，獲取算數平均值、中位數等統(tǒng)計特征值；對降水量以100 mm為檔統(tǒng)計次數，研究分布規(guī)律。

采用Matlab對近65年臨沂市年平均降水量進行小波分析，具體程序參考《Matlab小波分析工具箱原理與應用》[6]。Morlet連續(xù)小波變換圖反映降水在不同時間尺度上的周期震蕩，信號的強弱通過小波系數的大小來表示。小波光譜圖用顏色的深度代表偏離平均值的程度，紅色代表降水偏少，藍色代表降水偏多，小波系數為零則對應著突變點。小波方差圖反映了能量隨時間的分布，可以確定1個時間序列中各種尺度擾動的相對強度，對應峰值處的尺度稱為該序列的主要時間尺度，用以反映時間序列的主要周期[7-8]。

2 結果與分析

2.1 數理統(tǒng)計分析

1951—2016年臨沂市年均降水量的年際變化如圖1所示?？梢钥闯觯?5年，臨沂市年際變化較大，年均降水量最低值為529.5 mm（1988年），最高值為1 415.6 mm（1960年），平均值為864.3 mm，中位數為846.5 mm。

對1951—2016年臨沂市降水量進行分組，分析其分布規(guī)律。從圖2可以看出，年降水量基本符合正態(tài)分布，正態(tài)分布的相關系數R2=0.902 0，擬合結果如下：

y=y0+Ae■（1）

式中，y0=0.722 5，A=14.779 8，μ=3.854 3，σ=1.614 5。當x=3.854 3（對應降水量為835.4 mm）時，y值最大，為15.502 3（即降雨次數最多，概率為23.49%）。95%置信區(qū)間為0.436 9～7.271 2，對應降水量為493.7～1 177.1 mm。

2.2 Morlet連續(xù)小波分析

利用Matlab對1951—2016年臨沂市年平均降水量進行Morlet連續(xù)小波分析。不同周期上年均降水量的小波系數分布如圖3所示。臨沂市年平均降水量在7年、15年、28年、40年左右周期振蕩，有明顯的規(guī)律。在7年尺度上，存在明顯的小周期振蕩規(guī)律；在15年尺度上，存在“少—多”交替的8個階段，2016年降水偏少的等值線中心已經閉合，但多值區(qū)域尚未結束；在28年尺度上，存在“少—多”交替的3個階段，2016年少值區(qū)域已經出現(xiàn)，但中心尚未閉合；在40年尺度上，存在“少—多”交替的2個階段，規(guī)律更加清晰，2016年降水少區(qū)等值線中心已經閉合，但多值區(qū)域尚未結束；從15年、28年、40年3個尺度來看，臨沂市年均降水偏少的趨勢將結束，未來幾年降水量將增加。

2.3 小波方差

從圖4可以看出，1951—2016年臨沂地區(qū)降水量存在3個較為明顯的峰值，依次對應著15年、28年、40年。其中，最大的峰值對應著40年左右的時間尺度，說明在該尺度附近的周期振蕩最強，為降水量變化的第1主周期；28年時間尺度對應著第2峰值，為降水量變化的第2主周期，第3峰值對應著15年時間尺度。上述3個周期的波動控制著整個時間域內降水量的變化特征。

根據小波方差檢驗的結果，繪制出15年、28年、40年尺度的小波實部過程線。由圖5可以看出，在15年特征的時間尺度上，降水量變化的平均周期約為10年，大約經歷了6個“少—多”轉換期；在28年特征尺度上，降水量的平均變化周期為20年左右，大約經歷了3個“少—多”轉換期；而在40年特征尺度上，降水量的平均變化周期為27年左右，大約經歷了2個“少—多”轉換期。

3 結論

1951—2016年臨沂市年均降水量最低值出現(xiàn)在1988年，為529.5 mm；最高值出現(xiàn)在1960年，為1 415.6 mm；平均值為864.3 mm，中位數為846.5 mm，降水頻率符合正態(tài)分布。臨沂市年均降水量在15年、28年和40年左右周期震蕩劇烈，有明顯的規(guī)律；降水偏少的趨勢將結束，未來幾年降水量將逐漸偏多。40年時間尺度的小波系數的模方值和方差系數最大，能量最強、周期最顯著，為降水量變化的第1主周期。

4 參考文獻

[1] 王文圣，丁晶，向紅蓮.小波分析在水文學中的應用研究及展望[J].水科學進展，2002，13（4）：515-520.

[2] 馬柱國，華麗娟，任小波.中國近代北方極端干濕事件的演變規(guī)律[J].地理學報，2003，58（增刊1）：69-74.

[3] 劉忠陽，王勇，丁園圓，等.鄭州近54年降水變化的多時間尺度分析[J].氣象科技，2005（增刊1）：125-128.

[4] 邱海軍，曹明明，胡勝，等.近60a來中國洪澇災情變化趨勢持續(xù)性和周期性研究[J].地球與環(huán)境，2014，42（1）：17-24.

[5] 楊義，舒和平，馬金珠，等.基于Mann-Kendall法和小波分析中小尺度多年氣候變化特征研究：以甘肅省白銀市近50年氣候變化為例[J].干旱區(qū)資源與環(huán)境，2017，31（5）：126-131.

[6] 董長虹.Matlab小波分析工具箱原理與應用[M].北京：國防工業(yè)出版社，2004.

[7] 李遠平，楊太保.柴達木盆地近50年來年氣溫、降水的小波分析[J].干旱區(qū)地理（漢文版），2007，30（5）：708-713.

[8] 王濤，霍彥峰，羅艷.近300a來天山中西部降水與太陽活動的小波分析[J].干旱區(qū)研究，2016，33（4）：708-717.