趙志勇，國洪琴

(遼寧省葫蘆島水文局，遼寧葫蘆島 125000)

太子河流域降水演變特征淺析

趙志勇，國洪琴

(遼寧省葫蘆島水文局，遼寧葫蘆島 125000)

文章以太子河流域本溪水文站1952—2014年降水數(shù)據(jù)為依據(jù)，綜合應(yīng)用氣候傾向率分析、模比系數(shù)差積曲線對太子河流域降水量多年變化特征進行了研究，分析出本溪水文站降水量變化趨勢；汛前、汛期降水量演變特征，為本溪水文站乃至太子河流域的中長期降雨預報奠定基礎(chǔ)。

降水演變規(guī)律；降雨預報；本溪水文站；太子河流域

1 資料與方法

以厄爾尼諾現(xiàn)象為主的氣候變化對人類的生存環(huán)境和社會經(jīng)濟產(chǎn)生了深遠的影響，其中，對降水量的影響尤甚。遼河支流太子河流域是東北地區(qū)重要的經(jīng)濟增長地區(qū)，其降水量的變化趨勢也是水利科學研究的重點方向之一。

1.1 流域氣候概況及資料來源

太子河發(fā)源于本溪縣東草帽頂子山，本溪境內(nèi)河長168 km，流域面積4 428 km2，河道比降約80%左右。太子河流域多年平均降水量為824 mm，年內(nèi)降水量分布不均，多集中在汛期6—9月份。太子河流域多年平均水資源量為14.48億 m3，開發(fā)利用率為22.75%，主要集中在本溪市太子河下游地區(qū)。

本溪水文站是太子河控制站，設(shè)立于1933年11月，位于本溪市平山區(qū)體育街，集水面積4 324 km2，至河口距離252 km。

文章分析資料來源于太子河流域本溪水文站1952—2014年逐月降水資料，采用降水傾向率法和累積距平法進行分析，并采用Mann—Kenddall法進行趨勢檢驗。

1.2 降水傾向率法

降水量的長期趨勢變化，用線性回歸方程表示：式中：P為降水量，mm；a1=dP/dt為線性趨勢；a0為常數(shù)。

1.3 模比系數(shù)差積曲線法

模比系數(shù)差積曲線可以反映一段時期內(nèi)降雨的變化趨勢。其公式為：

式中：LPi為距平累計值；xi為第 i年的時間序列值為所計算時間序列的多年平均值。

1.4 Mann—Kenddall檢驗

Mann—Kenddall的檢驗方法是非參數(shù)方法。非參數(shù)檢驗方法亦稱無分布檢驗，其優(yōu)點是不需要樣本遵從一定的分布，也不受少數(shù)異常值的干擾，更適用于類型變量和順序變量，計算也比較簡便。

由于最初由曼（H.B.Mann）和肯德爾（M.G. Kendall）提出了原理并發(fā)展了這一方法，故稱其為曼—肯德爾（Man—Kendall）法。但是，當時這一方法僅用于檢測序列的變化趨勢。后來經(jīng)其他人進一步完善和改進，才形成目前的計算格式。下文采用該方法進行檢驗。

方法概述：

對于具有n個樣本量的時間序列X，構(gòu)造一秩序列：

可見，秩序列Sk是第i時刻數(shù)值大于j時刻數(shù)值個數(shù)的累計數(shù)。在時間序列隨機獨立的假定下，定義統(tǒng)計量：

式中：UF1=0，E Sk()，Var Sk()是累計數(shù)Sk的均值和方差，在x1，x2，…，xn相互獨立，且有相同連續(xù)分布時，它們可由式（5）算出：

2 降水量變化趨勢分析

2.1 降水傾向率法

2.1.1 年降水趨勢分析

根據(jù)降水傾向率法對本溪站降水量的變化趨勢進行了分析，見圖1（a）。降水傾向率b（mm/10 a）取10 a，當b＞0時，說明降水隨時間t呈上升的趨勢；當b＜0時，說明降水隨時間t呈下降的趨勢。

圖1 本溪水文站年降水傾向率圖

利用Mann—Kenddall法對本溪水文站多年平均降水量進行趨勢檢驗，給定顯著性水平α= 0.05，即U0.05=±1.96，趨勢檢驗結(jié)果表明：本溪站多年平均降水量1960—1985年有明顯下降的趨勢，1985—1995年有緩慢上升的趨勢，1995—2014年又有下降的趨勢。10年滑動平均曲線表示本溪站的降水傾向率為9.461 mm/10 a，R2=0.16。

2.1.2 汛前和汛期降水變化趨勢

現(xiàn)將1月1日至5月31日劃為汛前，6月1日至9月30日劃為汛期。利用Mann—Kenddall法對本溪水文站汛前和汛期降水量進行趨勢檢驗，給定顯著性水平α=0.05，即U0.05=±1.96。降水變化趨勢檢驗結(jié)果表明：

1）汛前：本溪站汛前1973—1990年有緩慢上升的趨勢，2000—2006年有下降的趨勢，其余年份變動不明顯。根據(jù)UF和UB曲線交點的位置，確定本溪水文站降水量減少是一突變現(xiàn)象，具體是從1990年開始的。10年滑動平均曲線表示本溪站汛前的降水傾向率為0.135 mm/10 a，R2=0.000 3。

2）汛期：由于本溪水文站汛期的降水量占全年降水量的比重較大，因此汛期降水量的變化趨勢與多年平均降水量的變化趨勢相似，見圖1（b）。由此可見，本溪站汛前降水量，1960—1985年有明顯下降的趨勢 ，1985—1995年有緩慢上升的趨勢，1995—2014年又有下降的趨勢。10年滑動平均曲線表示本溪站的降水傾向率為9.441 mm/10 a，R2=0.116 8。

2.2 年降水模比系數(shù)差積曲線法

2.2.1 年降水量周期性分析

對1952—2014年逐年降水量做模比系數(shù)差積曲線，由圖2（a）可以看出，年降水量變化大體可分為兩個時期，在20世紀50年代至70年代中期，總體上是增加的趨勢，70年代中期以后呈現(xiàn)周期震蕩減少的趨勢。

2.2.2 汛前和汛期降水變化趨勢

圖2 本溪水文站模比系數(shù)差積曲線圖

對本溪站1952—2014年汛前和汛期降水量做模比系數(shù)差積曲線分析，汛前降水量變化大體可分為3個時期，從20世紀50年代至60年代初呈現(xiàn)周期震蕩下降的趨勢，60年代至80年代中期呈現(xiàn)周期震蕩上升的趨勢，90年代中期以后呈現(xiàn)周期震蕩下降的趨勢。

汛期降水量變化大體可分為3個時期，從20世紀50年代至60年代中期呈現(xiàn)周期震蕩上升的趨勢，80年代中期至90年代中期又呈震蕩上升的趨勢，90年代中期以后呈現(xiàn)周期震蕩下降的趨勢。

3 結(jié) 論

以降水傾向率和模比系數(shù)差積曲線方法分析多年降水變化趨勢，并用Mann—Kenddall檢驗方法對分析結(jié)果進行趨勢檢驗，是我國近年來主要采用的方法。該方法經(jīng)過多次實踐，理論成熟、精度較高。

根據(jù)上述方法得出的降水量多年變化趨勢，在太子河流域中長期降水量預報中得到應(yīng)用，效果較好，為今后的防汛抗旱、減災避險等提供了數(shù)據(jù)基礎(chǔ)。

［1］欒兆擎，章光新，鄧偉，等.三江平原50年來氣溫及降水變化研究［J］.干旱區(qū)資源與環(huán)境，2007，21（11）：39-43.

［2］張雪花，步新影，孔繁德.氣候變化對秦皇島地區(qū)降水量影響分析［J］環(huán)境科學與管理，2007，32（8）：62-64.

［3］張劍明，黎祖賢，章新平.長沙近50年來降水的多時間尺度分析［J］.水文，2007，27（4）：78-80.

P337

A

1002－0624（2016）12－0033－03

2016-07-24