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

（1.湖南省水利水電勘測設計研究總院，湖南 長沙410007；2.湖南省水利水電職業(yè)技術學院，湖南 長沙410131）

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

1 流域概況

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

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

2.1 數(shù)據(jù)來源

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

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

2.2 研究方法

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

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

3 結(jié)果分析

3.1 徑流泥沙變化情況

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

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

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

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

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

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

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

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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