，，，，

(山東農(nóng)業(yè)大學(xué)水利土木工程學(xué)院，山東泰安271000)

受流域氣候因素、下墊面因素以及人類活動因素的影響，降水量等水文時間序列具有隨機性、非線性、多時間尺度變化等復(fù)雜的特性[1]，因而對于準(zhǔn)確的掌握其變化規(guī)律并做出合理的預(yù)測是眾多研究的困擾。目前，學(xué)界對于水文時間序列的研究有序列相關(guān)性分析方法、水文頻率分析方法、模糊分析方法等科學(xué)的方法[2-17]，但較多學(xué)者仍然通過序列相關(guān)性分析方法進(jìn)行研究。該方法基于水文時間序列的平穩(wěn)性假設(shè)，是一種研究序列變化趨勢和規(guī)律并對未來的情況進(jìn)行判斷和預(yù)測的有效的研究方法。蔡繼等曾在相關(guān)研究中[4]運用該方法對深圳水庫的降水量序列進(jìn)行了研究，結(jié)果表明深圳水庫年和汛期降水量呈增加趨勢且存在6、10、17 a的主周期，非汛期降水量呈微弱減少趨勢；袁滿等在相關(guān)研究中[5]運用該方法對沱江流域三皇廟站1941—2008年平均徑流序列進(jìn)行研究，結(jié)果表明顯著突變點出現(xiàn)在2001年。由此可見該方法具有科學(xué)性、有效性。本文采用序列相關(guān)性分析方法對太河流域的降水量序列進(jìn)行研究。

1 研究區(qū)概況

本文選取太河流域為研究區(qū)。太河流域位于山東省淄博市淄川區(qū)，東西部地勢較高，中部地勢較低，南部是東西走向的淄博冠魯山。流域面積約780 km2，長度約40.5 km，寬度約19.3 km,干流長度約61.1 km。

太河流域?qū)儆谂瘻貛О霛駶櫦撅L(fēng)性氣候，春季多風(fēng)，夏季多雨，秋季多旱，冬季少雪，具有春旱、夏澇、晚秋又旱的氣候特點。流域多年平均降水量約697.2 mm，受氣候與地形影響，流域降水量時空分布不均，季節(jié)性變化明顯，汛期降雨約占全年總降水量的73.5%，且雨季主要集中在7—8月份，大暴雨洪水發(fā)生概率較大[18-19]。流域內(nèi)設(shè)有太河水庫、鎮(zhèn)后、郝峪等水文站，設(shè)有田莊、峨莊、李家莊、郭莊、趙莊、中郝峪、燕峪山頂、夏莊、南博山、南邢、石馬等雨量站，其中，田莊、趙莊雨量站位于流域內(nèi)太河干流上，峨莊、郭莊雨量站位于太河二級支流上，其余雨量站皆位于太河三、四級支流上，見圖1。

圖1 太河流域

2 理論與方法

水文時間序列是在水文要素隨時間變化過程中記錄產(chǎn)生的，降水量就是一種典型的水文時間序列。水文時間序列用符號Xt表示，其線性表示方法為：

Xt=Nt+Pt+St

(1)

式中Nt——水文序列的非周期成分(含趨勢、跳躍、突變成分)；Pt——水文序列的周期成分(含簡單周期、復(fù)合周期及近似周期)；St——隨機成分(含平穩(wěn)、非平穩(wěn)狀態(tài))。其中，趨勢成分多指水文時間序列在較長的一段時期內(nèi)具有一定的規(guī)則的變化，而跳躍成分多指水文時間序列從一種狀態(tài)急劇變化到另一種狀態(tài)的形式[20]。

2.1 趨勢成分分析

對水文時間序列的趨勢成分進(jìn)行分析時，先對序列趨勢成分的存在性進(jìn)行識別，若序列存在趨勢成分，則須對趨勢成分的顯著性進(jìn)行檢驗；否則，無須進(jìn)行顯著性檢驗。

2.1.1滑動平均法識別趨勢成分[20]

設(shè)有某水文序列x1,x2,x3,……,xn,對該序列的前期值和后期值取平均，得到新的序列yt,使得原序列光滑化，其相應(yīng)的數(shù)學(xué)表達(dá)式為：

(2)

若序列中存在趨勢成分，選定k值，利用yt即可將趨勢清晰地顯示出來。

2.1.2Kendall秩次相關(guān)法檢驗顯著性[21]

設(shè)有某水文序列x1,x2,x3,……,xn,確定序列中所有對偶值(xi,xj)(j>i)中xi

假設(shè)原序列無趨勢(H0)，給定顯著性水平α，將計算求得的U值與Uα/2值進(jìn)行比較，若U

2.2 跳躍成分分析

對水文時間序列的跳躍成分進(jìn)行分析時，先對序列跳躍成分的存在性進(jìn)行識別，即識別序列有無突變點，若序列存在突變點，則須對突變點的顯著性進(jìn)行檢驗；否則，無須進(jìn)行顯著性檢驗。

2.2.1有序聚類法識別突變點[22]

設(shè)有某水文序列x1,x2,x3,……,xn,假設(shè)該序列的突變點為τ，構(gòu)造目標(biāo)函數(shù)：

(3)

其中，Vτ、Vn-τ為突變點前后兩部分的離差平方和。當(dāng)S取極小值時對應(yīng)的τ即為突變點。

2.2.2滑動游程法檢驗顯著性[23]

突變點τ確定后，將水文序列x1,x2,x3,……,xn突變點τ前后兩部分分別用不同的字母表示，將原序列值從小到大排序并用相應(yīng)的符號代替，統(tǒng)計連續(xù)出現(xiàn)相同字母序列的個數(shù)即游程記為k，前面部分的序列個數(shù)為n1，后面部分的序列個數(shù)為n2，構(gòu)造如下統(tǒng)計量：

(4)

給定顯著性水平α，將計算求得的U值與Uα/2值進(jìn)行比較，若U

3 數(shù)據(jù)與分析

選取太河流域內(nèi)太河水庫、鎮(zhèn)后、郝峪3處水文站以及田莊、峨莊、李家莊、南邢、石馬5處雨量站1979—2012年共計34 a實測降水量資料。

3.1 資料的三性審查

分別選用郝峪和李家莊2個站為代表站、石馬和峨莊為其相應(yīng)的鄰近站進(jìn)行可靠性分析。由以上2個站的逐年降水量相關(guān)性分析可以判斷兩個代表站與相應(yīng)的鄰近站逐年降水量的差值均較小，因此該流域降水量資料的可靠性較好；采用單累計模比過程線法分析降雨資料的一致性，郝峪、李家莊2個代表站的逐年降水量模比系數(shù)過程線的總趨勢均呈單一直線關(guān)系，相關(guān)系數(shù)分別為0.9944、0.9973，相關(guān)性較好，因此該流域的降水量資料的一致性較好；采用差積曲線法、累計平均值模比系數(shù)過程線法、長短系列相對誤差分析法求得逐年降水量資料的代表期和穩(wěn)定期為34 a。綜上所述，太河流域1979—2012年降水量資料的可靠性、一致性與代表性均可以保證。

借助GIS系統(tǒng)對實測的8個站的年降水量資料采用泰森多邊形法得到太河流域的年降水量資料序列。對該序列進(jìn)行趨勢和跳躍成分分析。該序列過程線見圖2。

圖2 太河流域1979—2012年逐年降水量過程線

由圖2，太河流域1979—2012年的降水量序列中，1989年出現(xiàn)降水量的最小值為252.4 mm，1990年出現(xiàn)降水量的最大值為808.1 mm。該水文時間序列的最大波動幅度出現(xiàn)在1989—1990年之間，極差為555.7 mm。從過程線的整體角度來看，該水文時間序列整體趨勢不明顯且于1989—2006年之間出現(xiàn)較為明顯的波動，僅通過過程線無法科學(xué)合理的確定該序列的具體組成成分，遂使用序列相關(guān)性分析方法對該序列的趨勢和跳躍成分進(jìn)行識別和顯著性檢驗。

3.2 趨勢成分分析

3.2.1滑動平均法識別趨勢成分

對太河流域該水文時間序列使用滑動平均法進(jìn)行趨勢擬合，采用滑動平均法得到結(jié)果見圖3。圖3中的滑動平均過程線顯示，太河流域的降雨量在1987—1997年、2002—2007年出現(xiàn)較明顯的波動，但在2007—2012年呈現(xiàn)較為平穩(wěn)的上升趨勢；總的來看，該水文時間序列具有上升的趨勢，綜合趨勢線走勢判斷該序列具有上升趨勢。

圖3 滑動平均法回歸趨勢線

3.2.2Kendall秩次相關(guān)法檢驗顯著性

用Kendall秩次相關(guān)檢驗法對該水文時間序列的上升趨勢進(jìn)行趨勢成分顯著性檢驗。取顯著性水平α=5%，假設(shè)原序列上升趨勢成分不顯著，計算得U=2.179，由顯著性水平α=5%查正態(tài)分布表得Uα/2=1.96。故U>Uα/2，則拒絕原假設(shè)，該序列上升趨勢具有顯著性。

3.3 跳躍成分分析

3.3.1有序聚類法識別突變點

對該水文時間序列采用有序聚類法對序列中的突變點進(jìn)行識別，經(jīng)處理得到有序聚類法統(tǒng)計曲線見圖4。由圖4，該水文時間序列的有序聚類法目標(biāo)函數(shù)S值在1996—2005年出現(xiàn)較明顯的波動，且波動幅度遠(yuǎn)大于整體的平均波動幅度，但在2005—

圖4 有序聚類法統(tǒng)計曲線

2011年呈現(xiàn)較平穩(wěn)的波動；綜合整體來看，S值在2002年有S=Smin，則突變點出現(xiàn)在2002年。

3.3.2滑動游程法檢驗顯著性

用滑動游程法對序列的跳躍成分進(jìn)行顯著性檢驗。經(jīng)處理得到滑動游程統(tǒng)計曲線見圖5。

圖5 滑動游程法統(tǒng)計曲線

由圖5，該水文時間序列的滑動游程法統(tǒng)計量U值在1985年有U=Umax=0.0013，在2009年有U=Umin=-0.0013即U值均滿足-0.0013

3.4 結(jié)果分析

本次研究表明，太河流域1979—2012年降水量序列具有顯著上升的趨勢成分與不顯著的跳躍成分。綜合考慮太河流域氣候因素、地質(zhì)因素以及人為因素等其他因素的影響，分析該水文時間序列的趨勢和跳躍成分的成因分析如下。

a) 自20世紀(jì)末，全球氣候變暖現(xiàn)象日趨明顯，由其引起的冰川消融、海平面上升、降水量增多等將進(jìn)一步導(dǎo)致洪、澇、旱等自然災(zāi)害[24]。相關(guān)研究[25-26]表明，全球氣候變暖主要是通過夏季風(fēng)影響我國的降雨分布。全球溫度升高，引起亞洲東部地區(qū)夏季風(fēng)強度增大，導(dǎo)致海陸溫差加大，使夏季風(fēng)邊界向北向內(nèi)陸伸展，使得中國受影響地區(qū)降水量增多。1979—2012年，隨著全球經(jīng)濟(jì)的不斷發(fā)展，全球氣候變暖趨勢不斷加強，太河流域受氣候變暖因素的影響，降水量整體呈現(xiàn)上升的趨勢。

b) 太河水庫是太河流域內(nèi)的一座大型水庫，較大的庫區(qū)面積對于氣候必定產(chǎn)生一定的調(diào)節(jié)作用，進(jìn)而對降雨的產(chǎn)生具有一定的影響。

c) 跳躍突變點的出現(xiàn)多與水文時間序列極值點的出現(xiàn)有關(guān)，受水文時間序列極值預(yù)報精度普遍較差的影響[27]，其誘發(fā)的自然災(zāi)害越來越嚴(yán)重。據(jù)有關(guān)資料[28]，2002年山東省全年平均降水量比歷年偏少，全省出現(xiàn)罕見的四季連旱現(xiàn)象，旱災(zāi)影響范圍較大，全省因干旱造成的直接經(jīng)濟(jì)損失達(dá)260億元，其中，淄博市即太河流域所處地級市春旱、夏旱、秋旱受災(zāi)耕地面積分別占全部耕地面積的66%、83%、61%，受災(zāi)害影響較大，因而出現(xiàn)跳躍突變點。

4 結(jié)論

a) 在未來的時期，太河流域的年降水量具有較大的概率出現(xiàn)逐年上升的趨勢。降水量的逐年增加，使得地表徑流量逐年增加，對地表的沖刷作用也越明顯[29]，故太河流域的下墊面條件會受到一定程度的破壞。

b) 降水量的增加，使得太河流域的水資源量有所增加。當(dāng)?shù)乜梢酝ㄟ^太河水庫對這部分水資源加以調(diào)節(jié)，實現(xiàn)有效利用，更好地為農(nóng)工業(yè)發(fā)展服務(wù)。從而可以有效地緩解由于社會發(fā)展、人口增長帶來的水資源短缺問題，實現(xiàn)水資源的供需平衡。

c) 受自然因素與人為因素的影響，對水文時間序列的預(yù)測具有一定的不確定影響?；诖?，盡管在2002年太河流域出現(xiàn)了特枯水年現(xiàn)象，但未來太河流域僅有一定的概率出現(xiàn)枯水年現(xiàn)象。