近60年來西江流域年輸沙量變化特征分析

,,,

(中水珠江規劃勘測設計有限公司，廣東廣州510610)

受氣候變化和人類活動影響，水文序列往往容易發生變異而導致其原變化趨勢發生改變。目前，在水文序列趨勢變化和突變分析上廣泛采用的方法是Mann-Kendall法[1]。此方法是一種非參數檢驗方法，可檢驗水文時間序列的變化趨勢和變化程度，即可揭示已發生的水文時間序列是上升還是下降趨勢和說明上升或下降趨勢變化是否顯著，但無法獲取水文序列未來的趨勢變化。大量的研究發現，水文時間序列具有統計自相似性，未來與過去存在相同或相反的變化特征，即具有分形特征[2-4]。R/S法是基于分形理論的提出的一種統計方法[5-6]，可用于分析獲取水文時間序列的分形特征，但無法確定水文序列未來趨勢變化是上升或下降。

鑒于R/S法和Mann-Kendall法分別在分析水文時間序列過去、未來趨勢變化上各有特點，本文嘗試將R/S法和Mann-Kendall法相結合，以西江梧州站1960—2015年年輸沙量時間序列為研究對象，力求揭示西江中上游年輸沙量時間序列的過去及未來的趨勢變化。

1 研究方法與資料

1.1 研究方法

水文時間序列趨勢變化分析常用分析方法有累積距平法、線性傾向估計法、滑動t檢驗法、游程檢驗法、Mann-Kendall法和R/S法等[7-8]。本文綜合運用R/S法和Mann-Kendall法分析水文時間序列過去及未來的趨勢變化特征。

1.1.1R/S法

R/S分析法，也稱重標極差分析法(Rescaled Range Analysis)，是由水文學家Hurst經大量實證研究提出的一種非參數方法，可用于揭示非線性系統如水文時間序列長期記憶性[7]。其原理和方法是對給定某一時間序列，t=1，2，…，可構建其極差和標準差：

(1)

赫斯特指數H(0

1.1.2Mann-Kendall法

Mann-Kendall法是基于秩的非參數檢驗法[9-10]，檢驗的統計量Z計算如下式：

(2)

當統計量Z大于0時，時間序列呈上升趨勢，小于0時呈下降趨勢。給定置信水平α，若|Z|≥Z1-α/2，則時間序列數據存在顯著的上升或下降趨勢。當α取5%時，|Z1-0.05/2|=1.96。

當進一步進行時間序列突變分析時，其統計量計算式如下：

(3)

UFk服從標準正態分布，給定顯著性水平α，如果|UFk|≥Uα/2，則表明序列趨勢變化顯著。將時間序列x按逆序排列，再按照上式計算，同時使：

(4)

綜合分析UFk和UBk可進一步獲取時間序列的趨勢變化、突變點以及突變區域。若UFk值大于0，表明時間序列呈上升趨勢；小于0則表明時間序列呈下降趨勢。當它們超過臨界值時，表明時間序列上升或下降趨勢顯著。超過臨界值的范圍確定為出現突變的時間區域。如果UFk和UBk兩條曲線出現交點，且交點在臨界直線之間，則交點對應的時刻便為突變開始的時間。

高職項目式教學要打破傳統的專業(學科)教學模式，按照高等職業教育的本質要求，從職業崗位的實際出發，依據工作的典型任務，選擇和安排學習項目，使得課程的學習內容與工作任務一致、學習過程與工作過程一致[6]。

1.2 研究資料

本文使用的年徑流量和輸沙量數據由水利部珠江水利委員會水文局提供，數據可靠。選取西江干流重要控制站梧州站為研究，資料年限接近60 a(1960—2015年)，系列較長，代表性好。

2 研究區概況

西江是珠江流域的主干流，從源頭至思賢滘全長2 075 km，平均坡降0.58‰，集水面積353 120 km2，占珠江流域總面積的77.8%。西江自馬雄山源頭下流，自西向東各河段依次稱為南盤江、紅水河、黔江、潯江和西江，總稱西江，西江流至思賢滘后匯入珠江三角洲河網區。西江為少沙河流，多年平均含沙量在0.032～0.570 kg/m3之間(梧州站)。雖然含沙量較少，但由于年徑流量大，因此輸沙量較大，多年平均徑流量和輸沙量分別為2 017億m3和5 412萬t(梧州站，見表1)。同時，輸沙年際變化大，年輸沙量變差系數Cv值為0.60，極值比高達33.6。

梧州水文站位于西江下游，西江干流與支流桂江匯合口以下約3 km處，集水面積327 006 km2，占西江流域總集水面積的92.6%?？梢?，梧州站是西江干流的重要控制站，分析研究該站的泥沙變化特征，對了解西江流域中上游水庫建設后清水下泄的對泥沙變化的影響有著十分重要意義。西江流域水系見圖1。

圖1 西江流域示意

項目時段多年平均Cv最大值最小值極值比輸沙量1960—2000年6 705萬t0.4214 037萬t1 691萬t8.32001—2015年1 879萬t0.645 341萬t418萬t12.81960—2015年5 412萬t0.6014 037萬t418萬t33.6徑流量1960—2015年2 017萬m30.192 961萬m31 024萬m32.9

3 結果分析

西江梧州站1960—2015年年輸沙量的M-K突變檢驗曲線見圖2，年徑流量M-K突變檢驗曲線見圖3，逐年輸沙量及徑流量變化過程見圖4。西江梧州站1960—2015年年輸沙量和徑流量趨勢特征見表2。

3.1 年輸沙量趨勢分析

采用Mann-Kendall法對梧州站年輸沙量進行檢驗可得，梧州站1960—2015年輸沙量統計量Z=-5.22，|Z|>|Z0.05/2|=1.96，因此，梧州站輸沙量序列呈顯著下降趨勢(表2)。

3.2 年輸沙量突變分析

3.2.1突變點

由圖2可知，自1995年以來，西江中上游年輸沙量呈減少趨勢，且2001—2015年的減少趨勢大大超過α=0.05的臨界線。根據UF、UB曲線的交點位置，可以確定西江中上游輸沙量的減少為突變現象，突變點為2000年。

表2 西江梧州站1960—2015年年輸沙量和徑流量趨勢特征

突變點前后輸沙差異較大。由表1可知，1960—2000年多年平均輸沙量為6 705萬t，2001—2015年多年平均輸沙量為1 879萬t，下降幅度達72.0%；變差系數Cv值由0.42增大為0.64，極值比由8.3增大為12.8。

3.2.2原因分析

河流輸沙量變化受自然因素、徑流變化和人類活動等影響，其中人類活動主要有水庫大壩修建、水土保持措施等。限于資料，本文僅從徑流變化及流域大型建設攔沙兩方面進行突變原因分析。

圖2 西江梧州站年輸沙量M-K突變檢驗

a) 徑流變化。對梧州站1960—2015年徑流量進行Mann-Kendall檢驗分析，力求獲取該站輸沙量下降趨勢及突變是否與徑流量變化有關。梧州站1960—2015年徑流量Mann-Kendall檢驗統計量Z=-0.77(表1)，徑流量呈下降趨勢。梧州站年徑流量UF值在α=0.05的臨界線之間波動(圖3)，未超過顯著檢驗，說明徑流量下降趨勢不明顯，這與圖4線性趨勢分析結果一致。另外，UF和UB曲線基本在臨界線內波動，雖有幾個交點，但變化趨勢不顯著，且年徑流量Cv較小(Cv=0.19，表1)，不存在突變現象。由圖4可知，梧州站1960—2015年的逐年輸沙量和徑流量均呈鋸齒形震蕩，但輸沙量總體呈下降趨勢，而徑流量變化趨勢不明顯。從20世紀90年代中期開始，輸沙量呈減少趨勢，但徑流量仍圍繞多年平均值(2 017億m3，表1)波動，即徑流量并未與輸沙量序列同步波動。

b) 流域大型水庫攔截泥沙。西江干流泥沙主要來源于黔江，而黔江泥沙主要來源為紅水河[11]。流域各梯級水庫庫容復核顯示，天生橋一級水庫年均淤積量1 081萬m3(1985—2005年)，龍灘水庫年均淤積量2 384萬m3(2007—2016年)，巖灘水庫年均淤積量1 322萬m3(1993—2016年)，大壩攔沙影響明顯。可見，西江中上游年徑流量呈不顯著下降趨勢，輸沙量呈顯著下降趨勢且發生變異，但輸沙量的減少與變異與徑流量變化關系不大，主要受流域大型水利樞紐工程建設攔沙影響。西江中上游20世紀90年代至21世紀初期興建了多座大型水利樞紐工程，如天生橋一級水庫1997年底下閘蓄水，龍灘水庫2006年9月底下閘蓄水，巖灘水庫1992年3月下閘蓄水，光照電站2007年12月31日下閘蓄水，百色水庫2005年8月26日下閘蓄水。

圖3 西江梧州站年徑流量M-K突變檢驗

圖4 西江梧州站1960—2015年年輸沙量與徑流量年際變化過程

3.3 年輸沙量未來趨勢分析

在應用R/S和Mann-Kendall法前，首先進行方法的可行性檢驗。對Mann-Kendall法揭示出水文時間序列的變異點前的序列進行分析，對變異點后序列進行檢驗，然后再結合整個序列的趨勢變化來預測其未來的趨勢特征。

3.3.1方法檢驗

西江梧州站年輸沙量變異點為2000年。采用Mann-Kendall法分析梧州站1960—2000年和2001—2015年輸沙量的趨勢特征，西江中上游1960—2000年Mann-Kendall檢驗統計量Z=-1.70，|Z|<|Z0.05/2|=1.96，輸沙量呈下降趨勢，但不顯著；2001—2015年Mann-Kendall檢驗統計量Z=-1.48，|Z|<|Z0.05/2|=1.96，輸沙量呈不顯著的下降趨勢。用R/S法計算梧州站1960—2000年序列H=0.8274>0.5，說明梧州站1960—2000年序列趨勢變化具有持續性，即2001—2015年的趨勢變化與1960—2000年相同，即如果1960—2000年呈上升趨勢，那么2001—2015年呈上升趨勢；相反，則有下降趨勢。根據Mann-Kendall法的分析結果，梧州站2001—2015年輸沙量序列呈不顯著的下降趨勢，與R/S法對1960—2000年序列的分析結果吻合。綜上所述，R/S法和Mann-Kendall法可應用于水文時間序列趨勢分析，且能很好地解釋水文時間序列未來的趨勢變化。分析結果見表1。

3.3.2分析結果

采用已檢驗可行的R/S法和Mann-Kendall法對西江梧州站1960—2015年的輸沙量序列未來的變化趨勢進行分析。Mann-Kendall法分析結果顯示，西江梧州站1960—2015年年輸沙量序列呈顯著下降趨勢，見表1。經R/S檢驗，西江梧州站1960—2015年年輸沙量序列極差與標準差關系為：

R(τ)/S(τ) =(0.3875τ)0.9232

(5)

可見，H=0.9232>0.5，說明西江中上游未來輸沙量的變化趨勢與1960—2015年一致，即西江中上游1960—2015年的變化趨勢具有持續性。

綜合R/S法和Mann-Kendall法的分析結果可得，西江流域中上游地區輸沙未來呈下降趨勢。西江梧州站年輸沙量這種趨勢變化特征可以為西江流域應對清水下泄對沿江防洪工程的影響提供科學的依據。

4 結語

本文對西江梧州站年輸沙量進行統計分析，并將R/S法和Mann-Kendall法應用于其變化特征分析，結果表明：西江中上游輸沙量在近60 a來年際變化大，呈顯著下降趨勢，且在2000年發生突變，其顯著下降趨勢與突變原因主要受流域大型水利樞紐工程的建設攔沙影響；西江中上游年輸沙量下降趨勢具有持續性。