1955－2009年澧水流域徑流泥沙演變及影響因素

劉碧維，宋 楠，耿勝慧，陳 茜

（1.湖南省水利水電勘測(cè)設(shè)計(jì)研究總院，湖南 長(zhǎng)沙410007；2.湖南省水利水電職業(yè)技術(shù)學(xué)院，湖南 長(zhǎng)沙410131）

澧水是湖南省四大河流之一，流域內(nèi)暴雨徑流量大，荒山跡地多，地質(zhì)構(gòu)造運(yùn)動(dòng)活躍，為湖南省水土流失嚴(yán)重地區(qū)之一［1］。隨著人類活動(dòng)不斷增強(qiáng)，洞庭湖各主要支流的徑流量和輸沙量發(fā)生了顯著變化，引起了國(guó)內(nèi)學(xué)者的關(guān)注。他們針對(duì)洞庭湖流域的徑流泥沙變化特征、演變規(guī)律和驅(qū)動(dòng)因素開展了一系列研究［2－7］。然而，以往的研究主要集中在洞庭湖全流域上，以澧水單個(gè)流域水沙演變?yōu)閷?duì)象的研究相對(duì)較少。因此本研究采用Pearson相關(guān)系數(shù)法對(duì)澧水年徑流量及年輸沙量的演變規(guī)律進(jìn)行分析，研究氣候變化（主要為降水）和人為活動(dòng)（主要是水土保持綜合治理）兩大因素對(duì)澧水年徑流量和年輸沙量的影響，以期為深入認(rèn)識(shí)澧水乃至洞庭湖流域水沙演變規(guī)律提供思路和參考。

1 流域概況

澧水發(fā)源于湖南省桑植縣杉木界（北源），流經(jīng)桑植、張家界、慈利、石門、澧縣、津市等縣市，于小渡口注入西洞庭湖，全長(zhǎng)390km，流域面積18 583km2。該流域主要支流有溇水、渫水、道水和涔水。澧水流域?qū)賮啛釒Ъ撅L(fēng)濕潤(rùn)氣候區(qū)，氣候溫和，雨量豐沛，四季分明，光照充足；年平均溫度16－18℃，多年平均降水量1 200～1 600mm。流域內(nèi)由于大氣環(huán)流和復(fù)雜地形的影響，降水量年內(nèi)變化大，且集中程度高，4—9月降水量占年總量的72%～78%。位于桑植境內(nèi)的八大公山年均降水量可達(dá)2 300mm，是湖南省三大暴雨區(qū)之一。流域內(nèi)植被主要分布在渠水、婁水、澧水北源源頭區(qū)，廣大中下游區(qū)植被稀少，生態(tài)效益低下［8］。流域范圍內(nèi)山高坡陡，河流縱坡較大，斷層發(fā)育，河谷深切，總體地勢(shì)西北高，東南低［1］。

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

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

年徑流量和年輸沙量選用澧水流域把口站三江口水文站1955—2009年同步資料；年降水量選用1955—2009年五道水、涼水口、南岔、大庸、溪口、三江口、石門7個(gè)雨量站逐年降水?dāng)?shù)據(jù)（圖1）。由于雨量站站網(wǎng)密度水平較低，本研究采用泰森多邊形法計(jì)算澧水流域的面平均降雨量［9］，資料均來(lái)自湖南省水文局。

圖1 研究區(qū)雨量站示意圖

2.2 研究方法

運(yùn)用Pearson相關(guān)系數(shù)法來(lái)檢驗(yàn)降水、徑流、輸沙與年份的相關(guān)關(guān)系。Pearson相關(guān)系數(shù)描述了兩定距變量間聯(lián)系的緊密程度，Pearson相關(guān)系數(shù)r是反映兩隨機(jī)變量線性相關(guān)性的統(tǒng)計(jì)量，r＞0表明兩變量正相關(guān)，r＜0表明兩變量負(fù)相關(guān)，r的絕對(duì)值越大表明相關(guān)性越強(qiáng)。其表達(dá)式為：

對(duì)于樣本容量為n的兩樣本序列，在已知自由度ν＝n－2和給定顯著性水平α的情況下，可通過查相關(guān)系數(shù)表的方法對(duì)r進(jìn)行顯著性檢驗(yàn)，若｜r｜大于相關(guān)系數(shù)檢驗(yàn)值，則拒絕原假設(shè)H0，即認(rèn)為兩變量線性相關(guān)是顯著的；否則就接受原假設(shè)H0，即認(rèn)為兩變量線性相關(guān)性不顯著［10］。

3 結(jié)果分析

3.1 徑流泥沙變化情況

3.1.1 年徑流量和輸沙量序列的趨勢(shì)性分析 根據(jù)澧水流域把口站三江口水文站1955—2009年輸沙量和徑流量數(shù)據(jù)，點(diǎn)繪其年徑流量、年輸沙量及相應(yīng)5a滑動(dòng)平均值序列（圖1）。

由圖2可以看出，澧水流域年徑流量和年輸沙量年際波動(dòng)很大，但呈現(xiàn)基本一致的變化規(guī)律，即大水大沙。但由于其隨年際波動(dòng)變化現(xiàn)象明顯，從中并不能直觀地看出其變化規(guī)律。

為了探尋年徑流量和年輸沙量隨年份的長(zhǎng)期變化規(guī)律，本文借助于Pearson相關(guān)系數(shù)法對(duì)年輸沙量和徑流量進(jìn)行趨勢(shì)顯著性檢驗(yàn)。經(jīng)計(jì)算，澧水1955—2009年年輸沙量和徑流量的Pearson相關(guān)系數(shù)r分別為－0.35和－0.09。取顯著水平α＝0.01，則自由度ν為53時(shí)相應(yīng)的檢驗(yàn)臨界值tα＝0.345。由于輸沙量｜r｜＞tα，因此，澧水流域近55a來(lái)年輸沙量在0.01顯著性水平上呈現(xiàn)明顯下降趨勢(shì)，而年徑流量沒有明顯的變化趨勢(shì)。

圖2 1955－2009年澧水流域年徑流量和年輸沙量年際變化

3.1.2 年徑流量和輸沙量序列的階段性分析 降水—輸沙雙累積曲線常被用于檢測(cè)輸沙在年際間波動(dòng)變化規(guī)律。即如果它們之間斜率變大，即同樣降水量所對(duì)應(yīng)的輸沙量增多，也就是說(shuō)區(qū)域輸沙增強(qiáng)，而斜率變小即意味著同樣降水所造成的流域輸沙量在減少［11］。繪制澧水流域累積降水量和累積輸沙量關(guān)系曲線（圖3）可以發(fā)現(xiàn)，澧水流域雙累積曲線在1983年之前基本呈直線，但之后曲線斜率明顯減小，說(shuō)明澧水流域從1983年開始，其輸沙強(qiáng)度開始下降。根據(jù)曲線斜率的變化可將澧水流域年降水—輸沙相關(guān)關(guān)系分成1955—1983年和1984—2009年2個(gè)時(shí)段。

降水—徑流雙累積曲線常被用于檢測(cè)徑流量在年際間波動(dòng)變化規(guī)律。即如果它們之間斜率變大，即同樣降水量所對(duì)應(yīng)的徑流量增多，也就是說(shuō)區(qū)域徑流增強(qiáng)，而斜率變小即意味著同樣降水所造成的流域徑流量在減少。繪制澧水流域累積降水量和累積徑流量關(guān)系曲線（圖4）可以發(fā)現(xiàn)，澧水流域雙累積曲線在近55a來(lái)基本呈直線，說(shuō)明澧水流域近55a來(lái)徑流強(qiáng)度維持不變。

圖3 澧水流域年降水量和年輸沙量雙累積曲線

圖4 澧水流域年降水量和年徑流量雙累積曲線

3.2 澧水徑流泥沙演變的驅(qū)動(dòng)力分析

影響流域徑流量和輸沙量的因素可分為自然和人為因素，其中又以氣候（主要是降水）和人類活動(dòng)（如水利水保工程）對(duì)徑流、泥沙的變化影響較大［12］。3.2.1 降水對(duì)澧水徑流泥沙演變的影響 降水是流域產(chǎn)流產(chǎn)沙的動(dòng)力條件，其大小與時(shí)空分布決定了流域徑流量和輸沙量的大小。在一定地表?xiàng)l件下，降水量越大流域徑流量和輸沙量也越大。也就是說(shuō)，流域徑流量和輸沙量的變化與降水量的變化趨勢(shì)基本一致，只有在局部地區(qū)受水土保持?jǐn)r沙和水利工程蓄水等人為因素的影響時(shí)才會(huì)呈現(xiàn)出不一致的變化趨勢(shì)。

收集1955—2009年五道水、涼水口、南岔、大庸、溪口、三江口、石門7個(gè)雨量站點(diǎn)逐年降水?dāng)?shù)據(jù)，采用泰森多邊形面積加權(quán)后得到澧水流域年降水量年際變化情況（圖5）。由圖2、圖4可以看出，近55a來(lái)澧水年徑流量與年降水量變化的一致性較好，而年輸沙量與年降水量變化的一致性較差。進(jìn)一步借助Pearson相關(guān)系數(shù)法對(duì)年降水量進(jìn)行分析，結(jié)果表明，過去55a澧水流域年降水量并沒有明顯變化，不存在明顯趨勢(shì)（r＝－0.024，p＞0.1，n＝55）。同時(shí)澧水輸沙量呈明顯減少趨勢(shì)，故可以認(rèn)為研究期內(nèi)降水量對(duì)澧水流域輸沙量影響不顯著。3.2.2 水利水保工程對(duì)澧水徑流泥沙演變的影響澧水年徑流量在時(shí)間序列上（1955—2009年）下降趨勢(shì)不明顯，這表明澧水水系發(fā)育較好，產(chǎn)水量大，因生產(chǎn)生活發(fā)展引起的工農(nóng)業(yè)、居民用水量增加尚未對(duì)澧水河流水文特征及其系統(tǒng)功能產(chǎn)生根本性影響，入湖徑流量未發(fā)現(xiàn)顯著變化；澧水入湖輸沙量在時(shí)間序列上（1955—2009年）呈明顯減少趨勢(shì)，在降水量沒有增加的情況下，其減少的主要原因是水土保持綜合治理和水利工程攔沙等人為因素的作用。對(duì)水利水保工程的實(shí)地調(diào)查印證了上述結(jié)論。從1980—2000年，澧水流域已相繼完成“七五”期間擬定的八條小流域共403km2的治理工作。其中，生物措施治理了79.6km2，工程措施控制泥石流面積0.48km2，坡改梯2.3km2，封山育林工作也得到了基本落實(shí)；完成了“長(zhǎng)防”工程共2.92×105hm2的營(yíng)林任務(wù)，全面改造、改種1.20×105hm2的油桐林；水土流失區(qū)的輕度流失和劇烈流失的面積得到了基本控制、改善和治理，其中，強(qiáng)度流失區(qū)治理了40%，中度流失區(qū)治理了50%。

圖5 1955－2009年澧水流域年降水量年際變化

為發(fā)展農(nóng)田水利，澧水流域修建了一些大、中、小型水庫(kù)等蓄水工程，截止2009年12月底，澧水流域已建大、中、小型水庫(kù)9 000余座（大型水庫(kù)2座，中型水庫(kù)22座），其中澧水流域最大的兩座水庫(kù)江埡水庫(kù)（庫(kù)容1.74×109m3）和皂市水庫(kù)（庫(kù)容1.44×109m3）分別于1999年5月和2008年4月建成蓄水，對(duì)減少澧水入湖泥沙量發(fā)揮著重要作用。

在水土保持綜合治理和水利工程建設(shè)等人為驅(qū)動(dòng)力的綜合作用下，澧水入湖泥沙量每年相應(yīng)減少約4.59×106t。

4 結(jié)論

（1）近55a來(lái)澧水年徑流量沒有發(fā)生明顯變化，而年輸沙量呈顯著下降趨勢(shì)。根據(jù)階段性分析結(jié)果，認(rèn)為澧水泥沙的演變過程（1955—2009年）可分為2個(gè)階段，即1955—1983年和1984—2009年2個(gè)階段。

（2）對(duì)澧水水沙演變的驅(qū)動(dòng)力分析結(jié)果表明，近55a來(lái)澧水年徑流量與年降水量變化的一致性較好，而年輸沙量與年降水量變化的一致性較差。由于澧水年降水量無(wú)明顯變化趨勢(shì)，故可認(rèn)為其年輸沙量顯著減少的主要原因是人為作用，水土保持綜合治理工程的實(shí)施和水利工程攔沙導(dǎo)致澧水入湖泥沙量每年相應(yīng)減少約4.59×106t。

（3）1983年澧水年輸沙量發(fā)生突變的主要原因是澧水流域自20世紀(jì)80年代初開展了以小流域綜合治理為模式的大規(guī)模水土流失治理，水土保持治理初顯成效，加之一些水利工程的興建有效攔截了部分泥沙。

（4）在前人對(duì)洞庭湖全流域水沙變化研究的基礎(chǔ)上，本文以澧水單個(gè)流域?yàn)檠芯繉?duì)象，揭示了澧水流域徑流泥沙的演變規(guī)律，并對(duì)其變化的影響因素進(jìn)行了分析。由于澧水流域地形復(fù)雜，引起水沙變化的原因眾多，今后需對(duì)其進(jìn)行進(jìn)一步論證。

［1］ 李景保.澧水流域物質(zhì)侵蝕強(qiáng)度及其遷移特征［J］.水土保持學(xué)報(bào)，1990，4（2）：62-69.

［2］ 覃紅燕，謝永宏，鄒冬生.湖南四水入洞庭湖水沙演變及成因分析［J］.地理科學(xué)，2012，32（5）：609-615.

［3］ 李景保，王克林，秦建新，等.洞庭湖年徑流泥沙的演變特征及其動(dòng)因［J］.地理學(xué)報(bào)，2005，60（3）：503-510.

［4］ 王鐵生.洞庭湖水系河流泥沙量的初步研究［J］.水土保持通報(bào)，1982，2（4）：48-56.

［5］ 張祥志.洞庭湖水沙特性和泥沙淤積分析［J］.華東師范大學(xué)學(xué)報(bào)（自然科學(xué)版），1996，42（1）：63-69.

［6］ 鄒文發(fā).洞庭湖泥沙沉積與土壤侵蝕［J］.中國(guó)水土保持，1992，13（6）：18-21.

［7］ 黎昔春，張水云.洞庭湖的泥沙輸移特性［J］.泥沙研究，2003，48（4）：73-76.

［8］ 莫靖龍，夏衛(wèi)生，羅軼.澧水流域水土流失原因分析及保持措施初探［J］.安徽農(nóng)學(xué)通報(bào)，2009，15（5）：75-76.

［9］ 楊瑋，韋宏鵠.利用泰森多邊形法分析汾河流域中段降水量［J］.山西水利，2006，22（2）：72-73.

［10］ 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)［M］.北京：氣象出版社，1999：69-72.

［11］ 鄭海金，方少文，楊潔，等.近40年贛江年徑流泥沙變化及影響因素分析［J］.水土保持學(xué)報(bào)，2012，26（1）：28-31.

［12］ 丁文峰，張平倉(cāng)，任洪玉.近50年來(lái)嘉陵江流域徑流泥沙演變規(guī)律及驅(qū)動(dòng)因素定量分析［J］.長(zhǎng)江科學(xué)院院報(bào)，2008，25（3）：23-27.