澧水流域径流泥沙变化及驱动力研究

刘碧维，宋 楠，耿胜慧，陈 茜

（1.湖南省水利水电勘测设计研究总院，湖南长沙410007；2.湖南省水利水电职业技术学院，湖南 长沙410131）

澧水是湖南省四大河流之一，流域内暴雨径流量大，荒山迹地多，地质构造运动活跃，为湖南省水土流失严重地区之一[1]。随着人类活动不断增强，洞庭湖各主要支流的径流量和输沙量发生了显著变化，引起了国内学者的关注，他们针对洞庭湖流域的径流泥沙变化特征、演变规律和驱动因素开展了一系列研究[2-4]。然而，以往的研究主要集中在洞庭湖全流域上，以澧水单个流域水沙演变为研究对象的相对较少，且采用的统计方法也较为传统。因此本研究采用Pearson相关系数法对澧水年径流量及年输沙量的演变规律进行分析，研究气候变化（主要为降水）和人为活动（主要是水土保持综合治理）两大因素对澧水年径流量和年输沙量的影响，以期为深入认识澧水乃至洞庭湖流域水沙演变规律提供思路和参考。

1 流域概况

澧水发源于湖南省桑植县杉木界（北源），流经桑植、张家界、慈利、石门、澧县、津市等县市，于小渡口注入西洞庭湖，全长390 km，流域面积18583 km2。主要支流有溇水、渫水、道水和涔水。流域属亚热带季风湿润气候区，气候温和，雨量丰沛，四季分明，光照充足；年平均温度16～18℃，多年平均降水量在1200～1600mm，流域内由于大气环流和复杂地形的影响，降水量年内变化大，且集中程度高，4～9月降水量占年总量的72%～78%，位于桑植境内的八大公山年均降水量可达2 300 mm，是湖南省三大暴雨区之一[5]。

2 数据来源及研究方法

2.1 数据来源

年径流量和年输沙量选用澧水流域把口站三江口水文站1955～2009年同步资料，年降水量选用1955～2009年五道水、凉水口、南岔、大庸、溪口、三江口、石门等7个雨量站逐年降水数据，资料均来自湖南省水文局。

2.2 研究方法

图1 研究区雨量站示意图Fig.1 Schematic diagram of rainfall station in study area

运用Pearson相关系数法来检验降水、径流、输沙与年份的相关关系。Pearson相关描述了两定距变量间联系的紧密程度，Pearson相关系数r是反映两随机变量线性相关性的统计量，r>0表明两变量正相关，r<0表明两变量负相关，r的绝对值越大表明相关性越强。其表达式为：

对于样本容量为n的两样本序列，在已知自由度v=n-2和给定显著性水平α的情况下，可通过查相关系数表的方法对r进行显著性检验，若|r|大于相关系数检验值则拒绝原假设H0即认为两变量线性相关是显著的，否则就接受原假设H0即认为两变量线性相关性不显著[6]。

3 结果分析

3.1 近55年澧水径流泥沙演变规律

3.1.1 澧水年净流量和输沙量序列的趋势性分析 根据澧水流域把口站三江口水文站1955～2009年输沙量和径流量数据，点绘其年径流量、年输沙量及相应5年滑动平均值序列（见图2）。

图2 1955～2009年澧水流域年径流量和年输沙量年际波动Fig.2 Interannual fluctuation of annual runoff and annual discharge in Lishui River from 1955 to 2009

由图2可以看出，澧水流域年径流量和年输沙量年际波动很大，但呈现基本一致的变化规律，即大水大沙。但由于其随年际波动变化现象明显，从中并不能直观地看出其变化规律，为了探寻年径流量和年输沙量随年份的长期变化规律，笔者借助于Pearson相关系数法对年输沙量和径流量进行趋势显著性检验。经计算，澧水1955～2009年年输沙量和径流量的Pearson相关系数r分别为-0.35和-0.09。取显著水平α=0.01，则自由度v为53时相应的检验临界值tα＝0.345。由于输沙量|r|>t α，因此，澧水流域近55年来年输沙量在0.01显著性水平上呈现明显下降趋势，而年径流量没有明显的变化趋势。

3.1.2 澧水年径流量和输沙量序列的阶段性分析 降水-输沙双累积曲线常被用于检测输沙在年际间波动变化规律。即如果它们之间斜率变大，即同样降水量所对应的输沙量增多，也就是说区域输沙增强，而斜率变小即意味着同样降水所造成的流域输沙量在减少[7]。通过建立澧水流域累积降水量和累积输沙量关系（见图3）表明，澧水流域双累积曲线在1983年之前基本呈直线，但之后曲线斜率明显减少，说明澧水流域从1983年开始，其输沙强度开始下降，根据曲线斜率的变化可将澧水流域年降水-输沙相关关系分成1955～1983年和1984～2009年2个时段。

3.2 澧水径流泥沙演变的驱动力分析

影响流域径流量和输沙量的因素可分为自然和人为因素，其中又以气候（主要是降水）和人类活动（如水利水保工程）对径流、泥沙的变化影响较大[8]。下面重点分析降水和水利水保工程对澧水径流泥沙演变的影响。

3.2.1 降水对澧水径流泥沙演变的影响 降水是流域产流产沙的动力条件，其大小与时空分布决定了流域径流量和输沙量的大小。在一定地表条件下，降水量越大流域径流量和输沙量也越大。也就是说，流域径流量和输沙量的变化与降水量的变化趋势基本一致，只有在局部地区受水土保持拦沙和水利工程蓄水等人为因素的影响时才会呈现出不一致的变化趋势。

图3 澧水流域年降水量和年输沙量双累积曲线Fig.3 The cumulative curve of annual precipitation and annual sediment discharge in Lishui River

收集1955～2009年五道水、凉水口、南岔、大庸、溪口、三江口、石门等7个雨量站点逐年降水数据，采用泰森多边形面积加权后得到澧水流域年降水量年际波动（图4）。由图2、图4可以看出：近55年来澧水年径流量与年降水量变化的一致性较好，而年输沙量与年降水量变化的一致性较差。进一步借助Pearson相关系数法对年降水量进行分析，结果表明，过去55年澧水流域年降水量并没有明显变化，不存在明显趋势（r=-0.024，p>0.1，n=55），同时澧水输沙量呈明显减少趋势，故可以认为研究期内降水量对澧水输沙量影响不显著。

图4 1955～2009年澧水流域年降水量年际波动Fig.4 Interannual fluctuation of annual precipitation in Lishui River from 1955 to 2009

3.2.2 水利水保工程对澧水径流泥沙演变的影响 澧水年径流量在时间序列上（1955～2009年）下降趋势不明显，这表明澧水水系发育，产水量大，因生产生活发展引起的工农业、居民用水量增加尚未对澧水河流水文特征及其系统功能产生根本性影响，入湖径流量未发现显著变化；澧水入湖输沙量在时间序列上（1955～2009年）呈明显减少趋势，在降水量没有增加的情况下，其减少的主要原因是受到水土保持综合治理和水利工程拦沙等人为因素的影响。对水利水保工程的实地调查印证了上述结论。从1980～2000年，澧水流域已相继完成“七五”期间拟定的八条小流域共403 km2的治理工作。其中，生物措施治理了79.6 km2，工程措施控制泥石流面积0.48 km2，坡改梯2.3 km2，封山育林工作也得到了基本落实；完成了“长防”工程共29.2万hm2的营林任务，全面改造、改种12万hm2的油桐林；水土流失区的轻度流失和剧烈流失的面积得到了基本控制、改善和治理，其中，强流失区治理40%，中度流失区治理50%。

为发展农田水利，澧水流域修建了一些大、中、小型水库等蓄水工程，截止2009年12月底，澧水流域已建大、中、小型水库9 000余座（大型水库2座，中型水库22座），其中澧水流域最大的两座水库江垭水库（库容17.4亿m3）和皂市水库（库容14.4亿m3）分别于1999年5月和2008年4月建成蓄水，对减少澧水入湖泥沙量发挥了重要作用。

在水土保持综合治理和水利工程建设等人为驱动力的综合影响下，澧水入湖泥沙量每年相应减少约459.03×107kg。

4 结论

（1）近55年来澧水年径流量没有发生明显变化，而年输沙量呈显著下降趋势。根据阶段性分析结果，认为澧水泥沙的演变过程（1955～2009年）可分为2个阶段，即1955～1983年和1984～2009年2个阶段。

（2）对澧水水沙演变的驱动力分析结果表明，近55年来澧水年径流量与年降水量变化的一致性较好，而年输沙量与年降水量变化的一致性较差。由于澧水年降水量无明显变化趋势，故可认为其年输沙量显著减少的主要原因是人为作用，水土保持综合治理工程的实施和水利工程拦沙导致澧水入湖泥沙量每年相应减少约459.03×107kg。

（3）1983年澧水年输沙量发生突变的主要原因是澧水流域自20世纪80年代初开展了以小流域综合治理为模式的大规模水土流失治理，水土保持治理初显成效，加之一些水利工程的兴建有效拦截了泥沙下泄。

[1]李景保.澧水流域物质侵蚀强度及其迁移特征[J].水土保持学报，1990，4(2)：62-69.

[2]覃红燕，谢永宏，邹冬生.湖南四水入洞庭湖水沙演变及成因分析[J].地理科学，2012，32(5)：609-615.

[3]李景保，王克林，秦建新，等.洞庭湖年径流泥沙的演变特征及其动因[J].地理学报，2005，60(3)：503-510.

[4]王铁生.洞庭湖水系河流泥沙量的初步研究[J].水土保持通报，1982，(4)：48-56.

[5]莫靖龙，,夏卫生，罗 轶.澧水流域水土流失原因分析及保持措施初探[J].安徽农学通报，2009，15(5)：75-76.

[6]魏凤英.现代气候统计诊断与预测技术[M].北京：气象出版社，1999.69-72.

[7]郑海金，方少文，杨 洁，等.近40年赣江年径流泥沙变化及影响因素分析[J].水土保持学报，2012，26(1)：28-31.

[8]丁文峰，张平仓，任洪玉.近50年来嘉陵江流域径流泥沙演变规律及驱动因素定量分析[J].长江科学院院报，2008，25(3)：23-27.