1951-2011年湘江徑流泥沙演變過程與特征分析

胡光偉， 毛德華， 李正最， 許 瀅

(1.湖南師范大學 資源與環境科學學院， 湖南 長沙 410081；2.湖南省水文水資源勘測局， 湖南 長沙 410007； 3.湖南工業大學， 湖南 株洲 412008)

湘江是洞庭湖流域最大的支流，是湖南省的母親河，對洞庭湖流域的水資源安全有較大貢獻。隨著人類活動的不斷增強，洞庭湖各主要支流徑流量和輸沙量均呈現出顯著變化，在國內學術界引起廣泛關注，學者們針對洞庭湖流域的徑流泥沙變化特征、演變規律和驅動因素開展了一系列研究[1-8]。但單獨對湘江流域徑流泥沙演變過程的系統研究較少，且采用的是傳統統計學方法。鑒于徑流量和輸沙量受多種影響因素制約，是非線性極強的隨機變量，利用傳統統計學方法很難準確地揭示出其內在變化規律，因此本文采用較為精確的時間序列分析方法中的Mann—Kendall趨勢分析法和均值差異T檢驗法對湘江年徑流量和年輸沙量的演變規律進行分析，為深入研究湘江流域及洞庭湖區水沙演變規律提供參考。

1 流域概況

湘江是湖南省最大的河流，為長江7大支流之一，流域內人口將近5 000萬，占湖南省總人口的60%左右。它發源于廣西壯族自治區東北部興安、靈川、灌陽、全州等縣境內的海洋山，上游稱海洋河，在湖南省永州市區與瀟水匯合，稱湘江，向東流經永州、衡陽、株洲、湘潭、長沙，至湘陰縣蘆林潭入洞庭湖后歸長江。全長856 km，湘江干流在湖南省內全長670 km，流域面積94 660 km2，其中湖南省內面積為85 383 km2，占總面積的90.2%，沿途接納大小支流2 157條，主要支流有瀟水、耒水、洣水、舂陵水、蒸水、漣水等，多年平均徑流量9.39×1010m3，實測最大流量20 800 m3/s。湘江流域屬于亞熱帶季風濕潤氣候，雨量豐沛，但年內分配不均，降水多集中在春夏季，多年平均降水量1 300～1 500 mm，降雨多集中在4—6月，占全年的40%～45%；7—9月干旱少雨，降水量約占年降水量的18%；1—2月最少，僅占全年的8%。湘江流域的洪水主要由氣旋雨形成，每年的4—9月為汛期，年最大洪水和洪峰水位多發生于每年的4—8月，其中5，6月出現次數最多，次洪歷時10 d左右，并且與暴雨發生時間相對應。

2 數據來源與研究方法

2.1 數據來源

采用湘江入洞庭湖控制站的湘潭水文站1951—2011年的年徑流量和年輸沙量數據以及湘潭站1956—2011年逐月徑流量和輸沙量數據資料，數據來源于湖南省水文水資源勘測局、《長江泥沙公報2000—2011》[9]和中國氣象科學數據共享服務網。

2.2 研究方法

(1)Mann—Kendall(M—K)非參數檢驗法。用于對長序列數據的趨勢突變進行檢驗分析，具體計算參照文獻[10—12]。M—K法因為不要求原始數據服從特定概率分布，只滿足時間序列隨機獨立即可，因此受到水文學者的廣泛認可。設定時間序列x1，x2，x3，…，xn，定義統計量：

(1)

式中：mi——第i個樣本xi大于xj(1≤j≤i)的累計數。

(2)小波分析。小波分析的基本思想是用一簇小波函數系來表示或逼近某一信號或函數，小波函數是指具有震蕩性、能夠迅速衰減到零的一類函數，即ψ(t)∈L2(R)且滿足：

(2)

式中：ψa,b(t)——基小波函數，它可通過尺度的伸縮和時間軸上的平移構成一簇函數系。

(3)

式中：ψa,b(t)——子小波，對于給定能量有限信號f(t)∈L2(R)，其連續小波變換(continue wavelet transform， CWT)為：

(4)

式中：Wf(a,b)——小波變換系數；f(t)——一個信號或平方可積函數；a——伸縮尺度；b——平移參數。

(3)均值差異T檢驗法。在徑流泥沙長時間序列中，若出現從某變點起前期及后期平均徑流量、輸沙量均呈現出明顯增加或減少的現象，稱為跳躍點或突變點。定義樣本長度為N的序列的突變指數AIj為：

(5)

(6)

式中：m1，m2——前后兩段序列樣本長度；Sp——聯合樣本方差。

(7)

式中：S1，S2——前后兩時段序列的方差。該統計量t服從自由度為m1+m2-2的t分布，當給出一定的顯著性水平ɑ，若t>tɑ，則在ɑ的顯著性水平上，基準點兩側m1年和m2年的均值有顯著差異，即在基準點處發生了突變，并取這一區間中最大突變指數AIj的年份作為突變的發生時間。

3 結果與分析

3.1 年徑流量和年輸沙量演變的趨勢性

圖1為1951—2011年湘江湘潭站年徑流量、年輸沙量變化過程和累積年徑流量、累積年輸沙量變化過程線。從圖1可以看出，湘江年徑流量和年輸沙量變化規律基本一致，即大水大沙，小水小沙；徑流量累積變化過程線基本呈現直線狀態，表明湘江徑流量隨時間沒有發生顯著趨勢性變化；輸沙量累積變化過程呈上凸形態，1980s之前基本上為直線，表明年輸沙量波動不大，平均年輸沙量1.14×107t，1980s以后輸沙量累積過程線開始偏離原增長趨勢，年輸沙量有減少趨勢，年平均輸沙量也相應減少至6.64×106t，減少幅度達41.91%。

圖1 湘潭站年徑流量和年輸沙量(a)及累積徑流量和累積輸沙量(b)變化過程

湘江多年平均徑流量為6.54×1010m3，其中汛期徑流量為3.89×1010m3，占全年的59.51%。年徑流量最大值1.04×1011m3(1994年)，年徑流量最小值2.81×1010m3(1963年)。

從圖2中可以看出，湘江徑流量圍繞多年平均值震蕩，以直線擬合時不滿足R2檢驗。湘江多年平均輸沙量為9.39×106t，年輸沙量最大值2.95×107t(1954年)，年輸沙量最小值1.27×106t(2011年)。湘江輸沙量擬合的趨勢線直線斜率為負，表示有較為明顯的衰減趨勢。

圖2 湘潭站徑流量和輸沙量線性擬合圖

3.2 突變分析

根據M—K檢驗結果(圖3)，湘江徑流量在1951—1955和1994—2011年徑流量處于增長的波動變化，1956—1993年徑流量處于波動減少變化，但并未突破95%置信水平。年徑流量UF和UB曲線相交于1955，1980，1983，1987等，且均落入0.05顯著性水平閥值內，即湘江徑流量出現突變可能發生在1955，1980，1983，1987年。因UF曲線未突破0.05顯著性水平閥值，表明徑流量沒有出現顯著性變化。湘江輸沙量在1951—1955，1961—1962年和1976—1987年呈現增加波動變化，其余時間均呈現輸沙量減少的波動變化。UF和UB曲線在信度區間內相交于1997年，說明在這一年輸沙量可能發生了突變現象，并且UF曲線在1999年突破95%的置信水平，表明輸沙量呈現出顯著減少的趨勢。

圖3 湘江湘潭站徑流量、輸沙量M-K統計值

為了驗證M—K突變分析結果的準確性，下面通過均值差異T檢驗方法檢驗可能發生突變年份的科學性。取顯著性水平ɑ為0.01，臨界值tɑ為2.704，根據均值差異T檢驗結果(表1)。當M=5，30，33，38時，突變指數AIj分別為0.242，0.085，0.012和0.152，相對應的徑流統計量t為1.032，-1.664，-0.096tɑ，意味著湘江輸沙量在第47 a附近發生了明顯的趨勢變化，即在1997年發生了一次泥沙突變。

表1 湘潭站年徑流量和輸沙量均值差異T檢驗統計量

3.3 周期分析

從圖4a,4b徑流量和輸沙量小波分析可以得出，徑流量變化過程中存在20～25，5～7，10～13和2～3 a 4類尺度的周期變化規律，其中在20～25 a尺度上出現了豐—枯交替的3次震蕩，而且該時間尺度的周期變化在整個分析時段表現出非常穩定的特性，具有全局性特征，相反，小尺度的周期變化波動性較強。由徑流量小波方差圖4b得出，湘江徑流存在20,7,10和3 a的周期性。而且20～25 a周期尺度上的方差最大，表明該時間尺度周期性變化最明顯，其他尺度上的周期性變化較小。輸沙量變化(圖4c， 4d)，湘潭站的輸沙量周期性與徑流量的周期較為一致，輸沙量有20～25，6～8，3～5和11～14 a 4類尺度的周期變化規律，其中20～25 a尺度上出現了豐—枯交替的3次振蕩，而且這個時間尺度的周期變化在整個統計時段表現出非常穩定的特性，具有全局性特征，而小尺度的周期變化波動性較強。從小波方差圖4d看，輸沙量存在20,7,3和13 a的周期性規律，其中第一主周期為20 a，其對應的小波方差最大。

3.4 徑流泥沙演變的階段性

徑流泥沙變化往往出現豐枯年交替的現象，一段時期表現出徑流豐沛、輸沙量大，一段時期表現出徑流較枯、輸沙量小。年徑流量和輸沙量的累積距平百分數能準確直觀地反映徑流量和泥沙量的階段性變化特征，計算公式為：

(8)

圖4 湘潭站徑流量和輸沙量小波分析等值線圖和小波方差

根據湘潭站年徑流泥沙序列階段性特征統計(表2)和圖5計算結果，可以把湘江徑流量序列變化分為4 個階段：1951—1954,1992—2002年為多水期，1955—1991,2003—2011年為少水期，這4個階段的多年平均徑流量分別為7.82×1010,6.15×1010,7.87×1010和5.94×1010m3；同樣湘江輸沙量序列可以分為2個階段：1951—1985年為多沙期，1986—2011年為少沙期，這2個階段的多年平均輸沙量為1.143×107和6.64×106t。湘江輸沙量演變的階段性與湘江上游興建水庫密切相關，湖南省最大的東江水庫曾于1958年動工，1978年4月復工，1980年11月成功截流合龍，1986年8月下閘蓄水，1987年11月投產發電，1992年樞紐全面竣工。進入20世紀，湘江輸沙量明顯減少的另一重要原因是湖南省水土保持工程的積極作用，說明水土保持工程措施發揮了不可替代的作用。

表2 湘潭站年徑流泥沙序列階段性特征統計

圖5 湘潭站徑流泥沙累積距平百分數

3.5 徑流量和輸沙量演變集中度與集中期

集中度與集中期分析法是最早應用于氣象領域內分析氣象因素在年內分配的一種向量方法，1982年湯奇成等[13]把集中度與集中期分析法引入到河流月徑流年內分配規律中來，并成功分析中國河川徑流年內分配規律，之后被應用到不同的流域中[14-17]。集中度反映了徑流泥沙年內的非均勻分布特征，集中期則反映了全年徑流量和輸沙量集中出現的時期(月份)，具體計算方法參照文獻[18]，計算結果見表3。

由表3可知，湘江徑流泥沙年內分配具有較為明顯的不均勻性，輸沙量不均勻性尤其明顯。徑流量集中度為38.2%，輸沙量集中度為57.2%，輸沙量年內分配不均勻性更顯著，輸沙量合成向量的方向為159.1°，也即輸沙量主要集中在每年的6月，能夠反映湘江輸沙量最大值出現的月份，與實際輸沙量年內統計數據一致。從相對變化幅度來看，輸沙量相對變化幅度為21.9，明顯大于徑流量，輸沙量絕對變化幅度為1.73×106t，也比徑流量絕對變化幅度(9.11×109m3)大。

表3 湘潭站徑流量和輸沙量集中度與集中期計算結果

4 湘江徑流泥沙演變的影響因素分析

本文對湘江干流湘潭水文控制站1951—2011年的水沙序列變化特征進行全面的分析，并根據已有的研究成果與湘江流域已注冊的水庫信息，對湘江流域60 a來水沙變化特征及其影響因素做了探討。

4.1 湘江流域氣候變化特征對湘江徑流泥沙演變過程的影響

湘江徑流主要來源于降水，年內分配不均勻，3—7月徑流量占全年的66.6%，其中5月最大，占全年的17.3%；8月至翌年2月徑流量占33.4%，其中1月最小，僅占全年的3.3%。湘江枯水徑流一年出現2次，第1次是10月至翌年2月的冬季枯水，這5個月平均徑流只占年徑流量的21.2%，湘潭站歷年實測徑流量最小值100 m3/s(1996年10月6日)；第2次是夏季枯水。湘江泥沙主要來自降水特別是暴雨對流域表土的侵蝕，汛期河流含沙量最大，且含沙量、侵蝕模數、水沙比自上而下不斷增大，水土流失現象也從上而下逐漸加劇。湘江多年平均侵蝕模數在100～600 t/km2。以老埠頭站以上來水的含沙量最大為92%，衡陽站至湘潭站區間來水的含沙量又比衡陽以上大51%。

總體上看，湘江流域年徑流量存在波動變化，但均無顯著變化特征，流域內年徑流量演變最大的影響因素是降水量的多少。李景保[19]和張劍明等[20]的研究表明，湘江的水位變化與流域的降水保持一致。毛德華[21]對湖南四水流域年最大流量的研究表明，自1985—1995年以來，四水流域年最大流量正處于一個較高的時期，而且這種狀態還將延續若干年，這與四水流域1996年和1988年出現了較大的洪峰流量是一致的。通過與降水周期比較后發現，湘江流域的徑流變化對降水有較強的依賴性，與降水的變化特征保持較為一致的趨勢。影響河流徑流的因素很多，流域內徑流變化主要受降水多年變化的影響，還同時受人類生產活動所引起的下墊面條件的改變，包括流域內水庫建設、居民和工農業用水量變化等綜合因素的影響。截至2001年湖南省興建了大中小型水庫13 318座[22]，其中大型水庫就達19座。湖南水庫總庫容達3.69×1010m3，包括興利庫容2.48×1010m3，占四水多年平均徑流量的15%，可使全省增加有效灌溉面積達1.55×106hm2。

4.2 水利設施建設和水土保持工程措施對湘江徑流泥沙演變的影響

湘江年徑流量在長時間序列1951—2011年變化趨勢不明顯，是因為湘江水系發育，流域支流眾多，降水豐沛，流域生產生活用水量并未對流域內徑流量產生根本性影響。2010年湘江流域工農業生產和城鄉居民生活用水總量為1.73×1010m3，并且每年都以7.9%的速度增加，其中，湘江流域規模以上企業用水總量就達5.10×1010m3。根據湘江湘潭站實測資料統計，湘潭站徑流量在時間序列上波動變化幅度相對較小，這表明湘江流域產水量大，水系發達，地表蒸發與植物蒸騰損失量小，而水庫的大量修建以及工農業生產和城鄉居民生活用水量的增加未對湘江流域的徑流產生根本性的影響。湘江干流的輸沙在時間序列上顯著減少，主要原因是產沙量的變化受人為和自然因素的雙重影響：一方面，水庫的蓄水作用攔截了大量泥沙淤積庫中，流域支流已建大型水庫10座，總庫容達1.19×1010m3，其中9座按省防汛辦計劃預留防洪庫容1.06×109m3，已建中型水庫中有116座留有防洪庫容4.00×108m3，這些水庫的建設對流域泥沙輸移產生重要影響。根據湖南省水利水電課題組的研究成果，湖南省四水流域水庫在20世紀60—80年代共截留泥沙達9.55×108t。另一方面，自20世紀50年代以來，湖南省水土流失綜合治理面積達22 296.4 km2，共營造水土保持林10 035.3 km2，對1 100條小流域開展綜合治理，建成一批水土保持示范工程，初步建立起覆蓋全省的水土保持監測網絡。在水利工程和水土保持措施綜合作用影響下，湘江泥沙的輸移量顯著減少。

5 結 論

(1)根據徑流量和輸沙量的趨勢性分析，湘江徑流量沒有發生明顯變化，而輸沙量下降趨勢顯著，輸沙量在1980s以后開始表現出明顯的下降，1985年之前年平均輸沙量為1.143×107t，而1985年之后多年平均輸沙量為6.64×106t，減少幅度達41.91%。

(2)通過Mann—Kendall法和均值差異T檢驗法對湘江湘潭站1951—2011年年徑流量和年輸沙量序列變化規律進行了分析，結果表明，湘江干流60 a徑流量在1951—1955和1994—2011年兩個時段有上升趨勢，1956—1993年呈現下降趨勢，但是變化并不顯著；湘江輸沙量在1997年發生了突變，呈現顯著下降趨勢，突變點相對應的統計量t值為3.520。

(3)由小波分析可知，湘江徑流量和輸沙量存在20 a的主要周期，并且，湘江徑流量和輸沙量還存在7 a的第二主周期。從第一、第二主周期來看，湘江徑流量和輸沙量的周期性變化規律非常吻合，可以進一步說明徑流量和輸沙量的相關關系非常顯著，徑流量是引起輸沙量變化的最重要因素。

(4)根據徑流泥沙階段性分析，可以把湘江徑流量序列變化分為4 個階段：1951—1954和1992—2002年為多水期，1955—1991和2003—2011年為少水期；同樣湘江輸沙量序列可以分為2個階段：1951—1985年為多沙期，1986—2011年為少沙期。而集中度與集中期計算結果表明，湘江徑流量和輸沙量年內分配不均勻性較為明顯，尤其以輸沙量的不均勻性更為顯著，輸沙量集中度為52.7%，主要集中在每年的6月，且輸沙量絕對變化幅度達1.731×106t。

(5)對湘江水沙演變過程的影響因素分析可知，60 a來湘江年徑流量與年降水量的變化具有一致性，輸沙量的變化主要受人為因素的影響，水利設施建設和水土保持工程措施是導致流域泥沙量減少的重要因素。

[參考文獻]

[1]高俊峰,張琛,姜加虎,等.洞庭湖的沖淤變化和空間分布[J].地理學報,2001,56(3):271-276.

[2]張曉陽,杜耘,蔡述明.洞庭湖演變趨勢分析[J].長江流域資源與環境,1995,4(1):64-68.

[3]郭小虎,姚仕明,晏黎明.荊江三口分流分沙及洞庭湖出口水沙輸移的變化規律[J].長江科學院院報,2011,28(8):80-86.

[4]李景保,張磊,王建,等.水沙過程變化下洞庭湖區的生態效應分析[J].熱帶地理,2012,32(1):16-21.

[5]盧承志.洞庭湖治理與江湖關系[J].湖南水利水電,2001(1):13-15.

[6]李景保,王克林,秦建新,等.洞庭湖年徑流泥沙的演變特征及其動因[J].地理學報,2005,60(3):503-510.

[7]劉卡波,叢振濤,欒震宇.長江向洞庭湖分水演變規律研究[J].水力發電學報,2011,30(5):16-19.

[8]林承坤.洞庭湖水沙特征與湖泊沉積[J].地理科學,1987,7(1):10-17.

[9]水利部長江水利委員會.長江泥沙公報2000-2011[M].武漢:長江出版社,2000-2011.

[10]覃愛基,陳雪英,鄭艷霞.宜昌徑流時間序列的統計分析[J].水文,1993(5):15-21.

[11]劉宇峰,孫虎,原志華.近60年來汾河入黃河水沙演變特征及驅動因素[J].山地學報,2010,28(6):668-673.

[12]符淙斌,王強.氣候突變的定義和檢驗方法[J].大氣科學,1992,16(4):482-493.

[13]湯奇成,程天文,李秀云.中國河川月徑流的集中度和集中期的初步研究[J].地理學報,1982,37(4):383-393.

[14]郭巧玲,楊云松,暢祥生,等.1957—2008年黑河流域徑流年內分配變化[J].地理科學進展,2011,30(5):550-556.

[15]丁瑾佳,許有鵬,潘光波.杭嘉湖地區城市發展對降水影響的分析[J].地理科學,2010,30(6):886-891.

[16]張秋孌,路紫.旅游網站信息流距離衰減的集中度研究[J].地理科學,2011,37(7):886-890.

[17]白紅英,侯欽磊,馬新萍.50年來秦嶺金錢河流域水文特征及其對降水變化的響應[J].地理科學,2012,32(10):1230-1235.

[18]劉賢趙,李嘉竹,宿慶,等.基于集中度與集中期的徑流年內分配研究[J].地理科學,2007,27(6):791-795.

[19]李景保,吳桂生,劉曉清.近數十年湘江流域河流水文變化規律[J].熱帶地理,1997,17(3):295-302.

[20]張劍明,黎祖賢,章新平.1960—2005 年湖南省降水的變化[J].氣候變化研究進展,2008,4(2):101-105.

[21]毛德華.洪災綜合風險分析的理論方法與應用研究[M].北京:中國水利水電出版社,2009.

[22]葉澤鋼.湖南省大型水庫特點及防洪格局[J].水資源研究,2003,24(1):20-23.