(1.郴州水文水資源勘測(cè)局，湖南 郴州 423000；2.湖南農(nóng)業(yè)大學(xué) 生命科學(xué)技術(shù)學(xué)院，湖南 長(zhǎng)沙 410128； 3.湖南省水文水資源勘測(cè)局，湖南 長(zhǎng)沙 410000)

資水是湖南第三大河，西以雪峰山與沅江分界，東隔衡山與湘江毗鄰，左源為發(fā)源于湖南省城步縣北茅坪黃馬界的赧水，右源為發(fā)源于廣西資源縣的夫夷水，兩河在邵陽(yáng)縣雙江口匯合后始稱資水。資水流經(jīng)邵陽(yáng)市納邵水，新邵以下納石馬江、大洋江、渠江，安化以下納數(shù)溪、洢溪、沂溪等支流，于益陽(yáng)市甘溪港注入洞庭湖。資水流域面積28 211 km2，全長(zhǎng)653 km，平均坡降1.44%。資水流域地勢(shì)西南高東北低，自西南蜿蜒流向東北，上下游海拔相差較大；河流大部分流經(jīng)丘陵與山谷，河谷兩岸山脈逼近，雨水集流快，河流水位陡漲陡落，具有山溪性河流的特征。

河川徑流受氣候、地貌、植被等自然條件以及人類活動(dòng)的綜合影響。隨著全球變暖和人類活動(dòng)影響的加劇，徑流發(fā)生了顯著的時(shí)空變化，直接影響了流域水資源的配置、開(kāi)發(fā)與利用，以及河流生態(tài)系統(tǒng)的物理、化學(xué)和生物過(guò)程[1-2]。為了實(shí)現(xiàn)水資源的優(yōu)化管理和可持續(xù)利用，研究水資源在氣候變化和人類活動(dòng)影響下的時(shí)空變化規(guī)律具有重大意義[3]。

本文運(yùn)用M-K趨勢(shì)檢驗(yàn)和雙累積曲線法討論資水水沙演變的趨勢(shì)性和突變性，并對(duì)成因進(jìn)行分析，以期為資水流域水資源管理、資水尾閭濕地生態(tài)保護(hù)提供參考。

1 數(shù)據(jù)來(lái)源和研究方法

1.1 數(shù)據(jù)來(lái)源

桃江站為資水入洞庭湖控制站，集水面積26 748 km2，控制資水流域面積的 94.8%。研究的主要數(shù)據(jù)包括資水1956～2016年的年降雨量、桃江站1956～2016年逐月徑流量和1959～2016年輸沙量。

1.2 研究方法

(1)采用Mann-Kendall非參數(shù)秩次相關(guān)檢驗(yàn)法(簡(jiǎn)稱“M-K法”)[4]，檢驗(yàn)徑流量和輸沙量時(shí)間序列的趨勢(shì)性和突變性。當(dāng)統(tǒng)計(jì)變量Z為正值時(shí)表示增加趨勢(shì)，為負(fù)值時(shí)表示減少趨勢(shì)，絕對(duì)值大于等于1.28,1.64,2.32時(shí)分別表示通過(guò)90%，95%，99%的檢驗(yàn)。用M-K法檢測(cè)突變時(shí)，UFK為標(biāo)準(zhǔn)正態(tài)分布，若UFK的絕對(duì)值大于Uα/2時(shí)，表示序列趨勢(shì)變化明顯。當(dāng)UFk與UBk在臨界線之間有交點(diǎn)時(shí)，交點(diǎn)所對(duì)應(yīng)時(shí)間為突變開(kāi)始的時(shí)間。

(2)采用雙累積曲線法檢驗(yàn)輸沙量和徑流量變化的趨勢(shì)轉(zhuǎn)折點(diǎn)。該方法是檢驗(yàn)兩個(gè)參數(shù)間關(guān)系一致性及其變化的常用方法。雙累積曲線是在直角坐標(biāo)系中繪制的同期內(nèi)一個(gè)變量的連續(xù)累積值與另一個(gè)變量連續(xù)累積值的關(guān)系線，它可用于水文氣象要素一致性的檢驗(yàn)、缺值的插補(bǔ)或資料校正，以及水文氣象要素的趨勢(shì)性變化及其強(qiáng)度的分析[5]。

2 資水水沙變化特征

2.1 年徑流量變化

桃江站1956～2016年的多年平均流量為 223.7×108m3，采用M-K法檢驗(yàn)，Z=0.70，無(wú)明顯變化趨勢(shì)。逐年徑流量變化呈波動(dòng)狀態(tài)(見(jiàn)圖1)，最小值(135.5×108m3)出現(xiàn)在1963年，最大值(358.9×108m3)出現(xiàn)在1994年。

圖1 1956～2016年資水徑流量的變化

圖2 資水不同年代徑流量年內(nèi)分配的變化

2.2 月徑流量變化

資水桃江站徑流量年內(nèi)分配不均，每年3月至8月為豐水期，9月至次年的2月為枯水期，徑流量最大值出現(xiàn)在5月，最小值出現(xiàn)在12月。不同年代的徑流量年內(nèi)分配變化見(jiàn)圖2。月徑流量最大值(89.16×108m3)出現(xiàn)在1998年6月，最小值(1.551×108m3)出現(xiàn)在1956年12月。1950s，1960s，1970s和1980s中連續(xù)最大月徑流量出現(xiàn)在4～6月，3個(gè)月的總徑流量分別占4個(gè)不同年代全年徑流量的55.4%，45.4%，45.5%和37.2%。1990s，2000s和2010s中連續(xù)最大月徑流量出現(xiàn)在5～7月，3個(gè)月的總徑流量分別占3個(gè)年代全年徑流量的41.8%，42.3%和53.5%。

M-K法檢驗(yàn)結(jié)果(見(jiàn)圖3)表明：①1月(Z=2.49)和7月(Z=2.39)通過(guò)置信水平檢驗(yàn)，月徑流量有明顯上升趨勢(shì)(見(jiàn)圖3(a))，突變時(shí)間均發(fā)生在1965年。1956～1964年1月平均徑流量為6.666×108m3, 1965～2016年1月平均徑流量為10.58×108m3，增加58.7%。②1956～1964年7月平均徑流量為14.59×108m3,1965～2016年7月平均徑流量為 28.36×108m3，增加94.4%(見(jiàn)圖3(b) )。③5月Z值為-1.86，呈顯著下降趨勢(shì)，突變時(shí)間發(fā)生在1959年。1956～1958年5月平均徑流量為 59.64×108m3,1959～2016年5月平均徑流量為33.42×108m3，減小44.0%(見(jiàn)圖3(c))。其他各月的Z的絕對(duì)值均小于1.64，無(wú)趨勢(shì)變化。

圖3 桃江站1，7月和5月徑流量M-K法檢驗(yàn)統(tǒng)計(jì)值

2.3 年輸沙量變化

桃江站1959～2016年的多年平均輸沙量為141×104t，總體呈波動(dòng)下降趨勢(shì)(見(jiàn)圖4)。年輸沙量最大值(583×104t)出現(xiàn)在1969年，最小值(2.47×104t)出現(xiàn)在2009年。

圖4 資水桃江站年輸沙量變化

M-K法檢驗(yàn)表明：Z值為-4.15，1959～2016年輸沙量呈明顯下降趨勢(shì)，突變時(shí)間為1999年(見(jiàn)圖5)。1956～1998年年平均輸沙量180×104t， 1999～2016年的年平均輸沙量53.8×104t，減少了 70.1%。

3 成因分析

3.1 徑流量變化

對(duì)資水年降雨量與年徑流量進(jìn)行了線性回歸分析[6]，二者高度相關(guān)(R2=0.732 8，見(jiàn)圖6)。用M-K法對(duì)1956～2016年年降雨量進(jìn)行分析，Z值為 0.89，年降雨量無(wú)明顯變化趨勢(shì)。分析結(jié)果表明，資水流域年降雨量在近60多年里保持相對(duì)穩(wěn)定，逐月徑流量的突變可能是人為調(diào)控的結(jié)果。眾多研究表明，水利工程的建設(shè)極大改變了流域徑流在上下游間和季節(jié)間的既有分配格局，導(dǎo)致自然水文情勢(shì)的變化[7-8]。資水流域現(xiàn)已建5座大型水庫(kù)：柘溪(Ⅰ級(jí)，)、筱溪(Ⅱ級(jí))、馬跡塘(Ⅱ級(jí))、修山(Ⅱ級(jí))和六都寨(Ⅱ級(jí))，均以發(fā)電為主，兼顧防洪、灌溉、航運(yùn)及其他綜合利用。其中，柘溪水庫(kù)是資水流域最大的水庫(kù)，總庫(kù)容是其他大型水庫(kù)總庫(kù)容之和的 6.9倍，占多年年均徑流量的16.0%。柘溪水庫(kù)修建于1960s初期，與逐月徑流量1月，5月，7月發(fā)生的突變時(shí)間吻合，水庫(kù)蓄水削峰防洪和水力發(fā)電等人為分配調(diào)節(jié)導(dǎo)致月徑流量發(fā)生顯著變化。

圖6 資水1956～2016年降雨量與年徑流量關(guān)系

3.2 輸沙量

1959～2016年徑流量和輸沙量雙累積曲線檢驗(yàn)(見(jiàn)圖7)表明， 轉(zhuǎn)折期出現(xiàn)在1999年， 與M-K法檢驗(yàn)結(jié)果一致。 采用雙累積曲線趨勢(shì)方程對(duì)資水桃江站1999～2016年輸沙量進(jìn)行估算， 資水入洞庭湖輸沙總量人為減少 2.124×107t， 年均減少118×104t， 占多年平均年輸沙量的83.7%。 導(dǎo)致輸沙量變化的主要原因包括徑流變化[9]、 水庫(kù)修建[10]和流域生態(tài)環(huán)境(水土保持工程)改善等。 柘溪水庫(kù)修建于20世紀(jì)60年代初， 該研究所用輸沙量資料始于1959年， 難以檢驗(yàn)柘溪水庫(kù)對(duì)資水流域年輸沙量的實(shí)際影響。 1956～2016年流域年降雨量和年徑流量的趨勢(shì)檢驗(yàn)表明， 二者皆沒(méi)有顯著趨勢(shì)性變化， 因此基本可以判定水庫(kù)修建和水土保持工程建設(shè)是導(dǎo)致資水流域年輸沙量減少的主要原因。 資水流域建成大中型水庫(kù)40余座， 總庫(kù)容19×108m3， 水庫(kù)建設(shè)是導(dǎo)致資水流域輸沙量變化的重要驅(qū)動(dòng)因素[11]。 1980s以來(lái)， 湖南省廣泛開(kāi)展綠化造林、 封山育林和生態(tài)公益林建設(shè)等[12]， 對(duì)改善流域生態(tài)環(huán)境， 減少資水流域產(chǎn)沙過(guò)程也發(fā)揮了重要作用。

圖7 資水1959～2016年年徑流量與年輸沙量雙累積曲線變化

4 結(jié) 論

基于1956～2016年資水桃江站長(zhǎng)時(shí)間水文系列資料分析結(jié)果表明：

(1)資水流域1956～2016年年徑流量總體呈波動(dòng)狀態(tài)，無(wú)明顯變化趨勢(shì)；1月與7月徑流量在1965年后有顯著增加趨勢(shì)；5月徑流量在1959年后有顯著下降趨勢(shì)；月徑流量格局的變化主要受柘溪水庫(kù)建設(shè)的影響。

(2)資水流域1959～2016年輸沙量總體呈明顯下降趨勢(shì)，突變時(shí)間發(fā)生在1999年，年輸沙量發(fā)生變化的主要原因是流域梯級(jí)水庫(kù)的修建和流域生態(tài)環(huán)境改善措施的實(shí)施。