，，，，

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

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

1 研究區(qū)概況

本文選取太河流域?yàn)檠芯繀^(qū)。太河流域位于山東省淄博市淄川區(qū)，東西部地勢(shì)較高，中部地勢(shì)較低，南部是東西走向的淄博冠魯山。流域面積約780 km2，長(zhǎng)度約40.5 km，寬度約19.3 km,干流長(zhǎng)度約61.1 km。

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

圖1 太河流域

2 理論與方法

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

Xt=Nt+Pt+St

(1)

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

2.1 趨勢(shì)成分分析

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

2.1.1滑動(dòng)平均法識(shí)別趨勢(shì)成分[20]

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

(2)

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

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

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

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

2.2 跳躍成分分析

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

2.2.1有序聚類法識(shí)別突變點(diǎn)[22]

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

(3)

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

2.2.2滑動(dòng)游程法檢驗(yàn)顯著性[23]

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

(4)

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

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

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

3.1 資料的三性審查

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

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

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

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

3.2 趨勢(shì)成分分析

3.2.1滑動(dòng)平均法識(shí)別趨勢(shì)成分

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

圖3 滑動(dòng)平均法回歸趨勢(shì)線

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

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

3.3 跳躍成分分析

3.3.1有序聚類法識(shí)別突變點(diǎn)

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

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

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

3.3.2滑動(dòng)游程法檢驗(yàn)顯著性

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

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

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

3.4 結(jié)果分析

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

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

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

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

4 結(jié)論

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

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

c) 受自然因素與人為因素的影響，對(duì)水文時(shí)間序列的預(yù)測(cè)具有一定的不確定影響。基于此，盡管在2002年太河流域出現(xiàn)了特枯水年現(xiàn)象，但未來(lái)太河流域僅有一定的概率出現(xiàn)枯水年現(xiàn)象。