李 揚

（山西省水利水電科學(xué)研究院 山西太原 030002）

0 引言

“十年九旱、旱澇交錯”是山西省的主要氣候特點。太原市是山西省省會，位于山西省中部盆地，面積6 988km2，屬溫帶半干旱大陸性季風(fēng)氣候。全市1951～2010年平均降水量為441.37mm，6～9月降水量占全年的70%以上，其中7、8月降水量占全年的40%，年內(nèi)及年際分布極不均勻，旱澇災(zāi)害時有發(fā)生。

標準化降水指數(shù)（Standardized Precipitation Index，SPI）是基于當(dāng)?shù)赜炅康拈L期氣候狀況，表征某時段某一降水量出現(xiàn)概率的指標之一，由McKee最早提出，可用實測降水序列直接計算。SPI值為正表示降水量偏多，SPI值為負表示降水量偏少，因此可根據(jù)SPI值的大小來劃分旱澇等級并評估旱澇程度。由于SPI具有計算簡便、輸入?yún)?shù)單一（僅降水量）、適用于不同時間尺度等優(yōu)點，在全球70多個國家的研究和工程實踐中得到了廣泛的應(yīng)用。

1 標準化降水指數(shù)SPI

我國通常采用分布Γ來描述降水量的變化。標準化降水指標（SPI）是指計算出某一時段內(nèi)降水量的分布概率后，對其進行標準正態(tài)化處理，最終用標準化降水累計頻率分布來劃分旱澇等級。其計算步驟如下：

假設(shè)某一時段降水量為隨機變量x，其Γ分布的概率密度函數(shù)為：

式中，α＞0為形狀參數(shù)，β＞0為尺度參數(shù)，Γ（）α為gamma函數(shù)。參數(shù)可用極大似然法估計，其估計值表達式為：

由于Γ分布未定義x=0的情形，而實測降水量可能為0，降水量為0時的概率用下式估計：

式中，n為降水量序列的長度，m為實測降水量x=0的樣本數(shù)。

對Γ分布進行標準正態(tài)化處理，即將式（5）或式（6）代入標準化正態(tài)分布函數(shù)，得：

式中的Z即為SPI。求解式（7），可得SPI值的計算公式：

依據(jù)國家氣象干旱等級標準（GB/T20481-2006）及Lloyd-Hughes和 Saunders的文獻[1]，確定以 SPI值劃分的旱澇等級共9級，如表1所示。

表1 SPI旱澇等級劃分

2 SPI與降水量關(guān)系

選擇1、3、6、9、12和24個月共 6 個時間尺度，分別繪制各時間尺度的SPI值（分別記作SPI-1，…，SPI-24）與相應(yīng)時段降水量的散點圖，如圖1所示?？梢钥闯?，各時間尺度的SPI均隨降水量的增加而呈增大趨勢。當(dāng)時間尺度小于1年（12個月）時，SPI與降水量的關(guān)系曲線不是單值的，而是一族單調(diào)遞增曲線，曲線的條數(shù)與年內(nèi)時段數(shù)一致。如：時間尺度為1個月時，年內(nèi)時段數(shù)為12，則相應(yīng)的SPI-降水量曲線為一族12條單調(diào)遞增曲線；時間尺度為3個月時，年內(nèi)時段數(shù)為4，則相應(yīng)的SPI-降水量曲線為一族4條單調(diào)遞增曲線。每一時段內(nèi)，SPI和降水量的關(guān)系為單調(diào)遞增，SPI-降水量關(guān)系按時段的不同分成幾部分，各部分曲線位置存在明顯差異，說明在時間尺度小于1年時，SPI能夠較好地反映年內(nèi)不同時段的旱澇狀況，而不受相鄰其他時段旱澇情況的影響。

圖1 不同時間尺度SPI與降水量的關(guān)系

3 基于SPI的太原市旱澇情況分析

3.1 基于SPI的各月旱澇情況

計算太原站1951～2010年各月的SPI，參照表1標準分類，得到太原站各月旱澇次數(shù)如表2所示。

由表2可知，太原屬于旱澇災(zāi)害頻發(fā)地區(qū)，在月尺度下，60年間太原地區(qū)共發(fā)生旱災(zāi)206次（其中特旱8次，重旱30次，輕中旱168次）；澇災(zāi)224次（其中特澇16次，重澇37次，輕中澇171次）；從旱澇發(fā)生的次數(shù)來看，太原的旱澇發(fā)生次數(shù)相差不大，且略偏澇。經(jīng)進一步分析發(fā)現(xiàn)：在降水量相對較小的非汛期，太原發(fā)生澇災(zāi)的次數(shù)相當(dāng)可觀，甚至超過了汛期；而旱災(zāi)的發(fā)生則多集中于降水量相對較大的汛期，尤其重旱和特旱的發(fā)生次數(shù)集中于4～10月，該時段正值作物生長期，旱災(zāi)對農(nóng)業(yè)生產(chǎn)危害較大。

因此，SPI能考慮到降水在時間分配上的差異，較為準確地反映了較小時間尺度下時段的實際旱澇情況，這也是太原地區(qū)在月尺度下澇略多于旱的原因，如在降水量偏少的冬春季，較小的降水也可能引發(fā)洪澇，而在降水量偏多的夏秋季，即使降水量較大也可能發(fā)生干旱。

3.2 基于SPI的各季度旱澇情況

計算太原站時間尺度為3個月的SPI如表3所示，用以分析太原地區(qū)年內(nèi)各季度旱澇發(fā)生情況。表3顯示，太原地區(qū)一年四季都有旱澇發(fā)生，各季度旱澇的總次數(shù)相差不大，但程度有所不同。第1季度旱澇多為交替出現(xiàn)，20世紀90年代以來中度以上旱澇次數(shù)呈增加趨勢；第2季度連年旱澇現(xiàn)象頻繁，20世紀80年代以前中度以上旱澇相對較多，此后澇的程度略有減輕，但中旱及重旱仍相對嚴重；第3季度是中旱及重旱多發(fā)期，此時太原地區(qū)高溫少雨，是造成干旱的主要原因；第4季度總體偏旱，但輕澇頻發(fā)，周晉紅等[2]對近年來太原市降水趨勢的研究結(jié)果表明，太原市春季、夏季、秋季降水呈減少趨勢，冬季降水呈增加趨勢，是導(dǎo)致第4季度輕澇多發(fā)的原因。

表2 太原站多年各月旱澇情況表

表3 太原站多年各季度旱澇情況表

3.3 基于SPI的各時間尺度旱澇情況

為反映太原地區(qū)不同時段的旱澇情況及其特征，通過分析1、3、6、9、12、24個月共 6 種不同時間尺度SPI值的變化趨勢來考察太原站1951～2010年共60年來的總體旱澇特征。

分析發(fā)現(xiàn)，時間尺度為6個月以下時，SPI的波動趨勢與相應(yīng)時間尺度降水量的波動趨勢不完全一致，受季節(jié)、氣候及降水特征的影響，存在著“降水量較大時澇災(zāi)不一定嚴重”、“澇災(zāi)嚴重時的降水量相對不大”等現(xiàn)象；當(dāng)時間尺度為6個月及以上時，SPI的變化趨勢與相應(yīng)時間尺度降水量的變化基本一致。

其中，時間尺度為12個月即1年時，SPI的變化趨勢表明：20世紀50年代中期、60年代、70年代初和90年代中期為太原地區(qū)旱澇災(zāi)害頻發(fā)時期，1956年、1964年、1966年、1969年和1973年都出現(xiàn)了SPI值接近甚至大于2的重澇和特澇，1957年、1965年、1972年、1986年、1997年和2005年都出現(xiàn)了SPI值小于-1.5的重旱和特旱，其余年份旱澇交替出現(xiàn)且變幅不大。這一結(jié)果與王濤等[3]采用Z指數(shù)分析太原站多年旱澇情況所得結(jié)論相同，說明當(dāng)時間尺度為1年時，SPI與Z指數(shù)對于旱澇情況的描述具有較好的一致性。史料顯示，1972年我國華北遭遇幾十年不遇的嚴重干旱，其中太原屬于特旱區(qū)，這一記載也與圖中所示1972年SPI值為-2.36的特旱等級吻合，亦可說明SPI對旱澇情況的描述較為準確可靠。

此外，較大的旱澇災(zāi)害往往在相鄰年份交替出現(xiàn)，如1956年和1957年先重澇后重旱，1964、1965和1966年則為特澇→重旱→重澇，1972年和1973年先特旱后重澇，因此，在實際中發(fā)生較大的旱災(zāi)時，需防范次年發(fā)生大型洪澇災(zāi)害的可能，而發(fā)生較大洪澇的年份需防范次年發(fā)生大旱的可能。

4 結(jié)論

本文采用太原站1951-2010年月降水序列，計算了多時間尺度下的SPI，并對SPI與降水量之間的關(guān)系及太原市旱澇規(guī)律進行分析，得出以下結(jié)論：

1）年內(nèi)時段數(shù)決定了SPI-降水量散點圖的分布及趨勢。

2）SPI可以較為準確地描述時段降水特征，反映降水在時間分配上的差異，從而在一定程度上避免由于降水時程分配不均對旱澇等級劃分結(jié)果產(chǎn)生的影響，對不同時間尺度的降水特征及旱澇分布描述更為客觀。

3）SPI能夠較好地反映太原市實際旱澇特征，是旱澇監(jiān)測的有效手段。

[1]Lloyd-Hughes B and SaundersmA.A drought climatology for Europe[J].International Journal of Climatology，2002，22：1571-1592.

[2]周晉紅，趙彩萍，常呈祥.太原市降水趨勢的初步研究[J].山西氣象，2003（2）：18-19.

[3]王 濤，朱敏嘉，崔愛萍，王計川.山西近45年降水變化對旱澇的影響分析[J].山西農(nóng)業(yè)科學(xué)，2010，38（12）：61-63.