覃红燕
, 谢永宏, 邹冬生
QIN Hong-yan, XIE Yong-hong, ZOU Dong-sheng
中图分类号: P349
文献标识码: A
文章编号: 1000-0690（2012）05-0609-07
通讯作者:
收稿日期: 2011-05-25
修回日期: 2011-07-22
网络出版日期: 2012-05-20
版权声明: 2012 《地理科学》编辑部 本文是开放获取期刊文献，在以下情况下可以自由使用：学术研究、学术交流、科研教学等，但不允许用于商业目的.
基金资助:
作者简介:
作者简介：覃红燕（1975-）,女,湖南石门人,博士研究生,主要从事生态水文研究。E-mail: qhy@hunau.net
展开
摘要
以湖南四水入洞庭湖的代表性控制站1951~2009年长系列水文数据为基础,运用Mann- Kendall趋势突变检验等方法探讨水沙演变过程并分析了其成因。结果表明：四水年径流量变化比较复杂,存在多个上升-下降过程,但总体上无明显上升或下降趋势;年输沙量总体上呈明显下降趋势,且突变时间存在差异,湘水是1996年,资水是1973~1974年,沅水是1997~1998年,澧水是1998年。年径流量的减少、植被覆盖增加以及大中型水库的建设等综合作用是2001年后年输沙量较大幅度减少的主要原因。
关键词：
Abstract
On the basis of hydrological data from 1951 to 2009 of four representative gauging stations (Xiangtan, Taojiang, Taoyuan, Shimen), the changes of runoff and sediment discharge into Dongting Lake from the Four Rivers (Xiang River, Zi River, Yuan River and Li River) in Hunan Province were investigated using Mann-Kendall trend test, moving average value and double mass analysis. The results showed that annual runoff was highest in Xiang River (657×108m3), intermediate in Yuan River (640×108m3) and Zi River (229×108m3), and lowest in Li River (147×108m3). However, the changes of annual runoff of the Four Rivers were insignificant (P > 0.05), with several up-down processes. Annual sediment discharge was highest in Yuan River (1 051×104t), intermediate in Xiang River (959×104t) and Li River (382×104t), and lowest in Zi River (199×104 t).The annual sediment discharge from the Four Rivers showed a significantly decreasing trend generally (P < 0.001). Changes of annual runoff and annual sediment discharge were consistent in Xiang River before 1980s, in Zi River after 1960s, in Yuan River before 1980s and in Li River before 1990s. Moreover, the time of abrupt change in annual sediment discharge was different among the Four Rivers, with Xiang River in 1996, Zi River in 1973-1974, Yuan River in 1997-1998 and Li River in 1998, and the time of abrupt beyond 95% confidence interval was 1999, 1963, 1988, 2001 in Xiang River, Zi River, Yuan River and Li River, respectively. Double mass analysis between annual runoff and annual sediment discharge showed that Xiang River had one turning point in 1991, Zi River three points in 1961, 1981 and 1997, Yuan River two points in 1979 and 1997, Li River one point in 1996. Turning points of annual sediment discharge in Four Rivers were closely related to the construction of large-scale reservoirs. The reduction of annual runoff, increase of vegetation cover and construction of medium- and large-scale reservoirs were the main reasons accounting for the great reduction of annual sediment discharges of Four Rivers after 2001.
Keywords：
洞庭湖位于湖南省东北部,长江中游荆江南岸,为中国第二大淡水湖,是长江中游最重要的行蓄洪湖泊之一,发挥着巨大的调蓄洪水和保护生物多样性等重要生态服务功能。洞庭湖总面积为2 625 km2,其多年平均地表径流达3 126×108m3,约为黄河的6倍,鄱阳湖的3倍,太湖的10倍。入湖径流按来源通常分为长江四口、湖南四水及区间三部分,其中来自四口的为1 119×108m3,占37.1％;来自四水的为1 647×108m3,占54.6％,在洞庭湖水资源组成中占有重要地位。近年来,由于气候异常、人类高强度干扰等多重因素的影响,湖泊泥沙淤积严重,湿地面积不断萎缩,洪涝灾害加剧,生物多样性受损,生物灾害频发等多种生态环境问题日趋严峻[1],而最近几年几乎年年都发生的水资源短缺事件尤为引人关注,并引发人们对湿地生态环境加剧恶化的担忧[2]。因此,阐明流域水沙演变趋势是湿地生态环境保护的重要基础性工作之一。
流域的水沙变化一直是河流地貌、水利工程等领域倍受关注的科学问题[3,4]。水文特征如水沙的流动对很多因素都很敏感,包括水利工程（如水库、大坝、水管理措施）、土地利用变化以及气候变化等[5~7]。目前,国内外对于河流水沙变化的影响因素已有大量研究,尤其是人类活动对河流水沙变化的影响方面较为广泛。这些研究对于揭示地表水减少过程、气候变化的水沙响应以及人类活动的干扰过程与调控机制有着重要意义[8,9]。然而,以往的研究主要集中在土地利用、水库拦沙、水土保持、植被覆盖等方面[10~12],且主要集中在单个河流上,以全流域为基础、以河流进入湖泊的水沙演变为研究对象的相对较少。基于此,本研究的主要目的是通过移动平均值法、Mann-Kendall趋势突变检验、双累积曲线等方法探讨湖南四水入洞庭湖水沙演变的规律性、突变性及其成因,以期为充分合理利用四水水资源,保护洞庭湖湿地生态环境提供依据。
长系列水沙数据来源于湖南省水文局和《长江泥沙公报2003~2009》[13]。本研究以湘潭水文站（湘水）、桃江水文站（资水）、桃源水文站（沅水）、石门水文站（澧水）1951~2009年的实测数据作为湖南四水入洞庭湖水沙演变分析的基础（图1）。主要指标包括年径流量和年输沙量。
图1 湖南四水代表性控制站位置
Fig.1 Location of representative hydrological stations of Four Rivers in Hunan Province
主要研究方法如下：
1） 采用移动平均值法判定水沙演变趋势,期数设定为5。
2） 采用非参数Mann-Kendall （以下简称M-K）法检验趋势突变点,具体计算方法参见[8,14]。M-K法以时间序列平稳为前提, 且时间序列是随机独立的,其概率分布等同。设定原假设H0：在时间序列没有变化的情况下,设此序列为x1、x2、x3…xm, 用mi表示第i个样本xi大于xj ( 1≤j≤i)的累计数, 则可以定义一统计量:
3） 采用双累积曲线法检验年径流量和年输沙量累积变化的转折点。同时,通过转折点的获取估计了不同时段的减沙量,即通过前一时段的线性模拟方程按后一时段最后1 a 的累积径流量换算得到的后一时段最后1 a的理论累积输沙量。减少输入总量为后一时段最后1 a的理论累积输沙量与实际累积输沙量的差值,平均减少量为减少输入总量与后一时段年数的比值。
湖南四水入洞庭湖代表性水文站在1951~2009年的总径流量多年平均为1 673× 108m3。年径流量以湘水最高（657×108m3,占四水总径流的39.3%）,其次为沅水（640×108m3,占38.3%）,再次为资水（229×108m3,占13.7%）,最少的为澧水（147×108m3,占8.8%）。从移动平均值来看（图2）,总体上呈波动状态,除澧水外,湘、资、沅均体现为三个波峰（1950 s初期、1960 s末至1970 s初期、1990 s）和三个波谷（1950 s末至1960 s初、1980 s、2001~2009）,且三个水系的最大径流量都在1990 s。相比之下,澧水波动更为频繁但幅度相对较小。
湖南四水入洞庭湖代表性水文站在1951~2009年总输沙量多年平均2 591×104t。与年径流量不同的是,年输沙量以沅水最高（1 051×104t,占四水总输沙量的40.6%）,其次为湘水（959×104t,占37.0%）,再次为澧水（382×104t,占14.7%）,最少的为资水（199×104t,占7.7 %）。从不同年代来看,四水的总输沙量呈不断下降趋势,如1951~1960年为3 543×104t,1981~1990年为2 280×104t,到2001~2009年仅为876×104t,下降75.3%。从移动平均值来看（图2）,总体上呈不断下降趋势,其中以沅水下降最多,2001~2009年比1951~1960年下降90.6%;其次为资水,下降88.3%;再次为澧水,下降73.5%;最少的为湘水,下降53.9%。从快速下降的年代来看,存在明显差异,湘水在1980 s,资水在1960 s,沅水在1980 s和2001~2009年,澧水在1990 s后期。从径流量和输沙量的变化趋势看,湘水在1980 s以前、资水在1960 s以后、沅水在1980 s年前、澧水在1990 s以前呈明显的一致性变化趋势。
图2 湖南四水入洞庭湖代表性水文站1951~2009年的年径流量和年输沙量移动均值
Fig.2 Moving average values of annual runoff and sediment discharge in representative hydrological stations of Four Rivers of Hunan Province in 1951-2009
从M-K检验结果可以看出（图3）,四水入洞庭湖代表性水文站年径流量曲线变化较为复杂,存在多个上升-下降过程,表现为波动状态。以湘水湘潭站为例,在经历了1954~1958年的下降后,1958~1963年有所上升,1963~1967年又有所下降,1967~1983年为波浪式上升,1983~1990年又有所下降,1993~2002年间的上升过程相对较大,此后又呈下降趋势。虽然四水的UFk和UFb都存在多个交叉点,但M-K曲线一直处于置信区间内波动（׀Z׀ <1.96, P上或P下>0.05）,说明变化并不显著。相比之下,澧水的波动较小,而湘水、资水和沅水的波动幅度相对较大。可以认为,四水的年径流量在近59 a里未表现明显的上升或下降的趋势性变化。
从M-K检验结果可以看出（图3）,四水入洞庭湖代表性站点年输沙量总体上存在极其显著的下降趋势（Z<-3. 29, P下<0.001）,但曲线变化过程较为复杂,存在多个上升-下降过程。以湘水湘潭站为例,1956~1965、1971~1986和1998~2009年有三个明显的下降过程,并在1987~1998年有一个上升过程。湘水、资水和沅水的波动幅度相对较大,而澧水的波动较小,在1951~1996年间UFk为小幅波动,其后呈明显下降趋势。从突变时间来看,也存在差异,湘水为1996年,1999年后突破95%置信区间;资水的突变时间为1973~1974年间,在1963年即突破95%置信区间;沅水的突变时间为1997~1998年间,1988年突破95%置信区间;澧水的突变时间为1998年,此后几乎呈直线下降,突破95%置信区间时间为2001年。此外,从四水年径流量和年输沙量的M-K曲线变化过程来看,湘水在1996年前基本一致,而在1997年后不一致。可以认为,四水年输沙量在近59 a里呈极显著的下降趋势,但突变时间存在极大差异。
图3 1951~2009年湖南四水代表性水文站年径流量（a, b, c, d）和年输沙量（e, f, g, h）的M-K检验统计值
Fig.3 Statistical values of M-K test for annual runoff (a, b, c, d) and sediment discharge (e, f, g, h) of the representative hydrological stations of Four Rivers in Hunan Province in 1951-2009
尽管湖南四水入洞庭湖年径流量存在波动,但都无明显的上升或下降趋势,与其它河流的研究结论是一致的[11]。流域内年径流量的变化可能主要受降水量丰缺的影响[2]。李景保等对湘江流域数十年河流水文变化规律的研究表明,湘江及其支流水位变化规律与流域降水变化完全一致[15]。四水在1990 s的最大径流量与同期的最大降水量是非常吻合的[16]。影响入湖径流因素众多,包括流域内水库、居民和工农业用水量变化等。至2008年,湖南省14个地州市共建有大中小型水库11 466座,约占全国水库总数1/9,总库容为392.6×108m3[17],占1951~2009年湖南四水多年平均径流23.5%。同时,随工农业发展和居民生活水平提高,用水量呈持续增长状态,如2009年湖南规模工业企业取水总量68.02×108m3,比2005年增长3.07%（http://www.hunan. gov. cn / zwgk / tjsj / sjfxjd / 201005/t20100505_206096.html）。因此,湖南四水入湖径流无明显的上升或下降趋势说明湖南四水流域产水量大,水系发达,水利工程蓄水和工农业用水量的增加未能对四水河流整体水文特征及其系统功能产生根本性的影响。
湖南四水年输沙量变化的原因可分两个方面,一是自然因素（如降雨）导致的年径流量变化;二是人类活动干扰。泥沙移动是伴随水文而发生的自然现象,即年径流量对年输沙量的调控作用明显。尽管湖南四水入洞庭湖的年径流量和年输沙量变化较复杂,均表现出年径流与年输沙量变化趋势高度一致现象,如1980 s以前湘水湘潭站、资水桃江站以及澧水石门站及1960 s以后沅水桃源站。
湘水湘潭站双累积曲线只有1991年一个转折点（图4a）。经计算,1991~2009年共减少泥沙输入11 095×104t,年平均为584×104t。在四水流域,470余座大中型水库的建设是影响年输沙量的主要因素之一[18]。在湘水流域,全省最大的东江水库（库容为91.48×108m3,占多年平均径流的13.9%）于1987年开始二次建设,于1991年开始蓄水发电,1992年全面竣工。建成运行后的东江水库对输沙量的减少作用非常明显,致使1993~2002年间湘潭站年输沙量比多年平均少17.9%,双累积曲线也在蓄水的当年发生趋势性变化。
资水桃江站年径流量和年输沙量的双累积曲线有三个转折点,即1961、1981和1997年（图4b）。经计算,从1961年到1981年共减少泥沙输入总量为5 639×104t,年平均为282×104t,减沙效应最为明显;1982~1998年共减少泥沙输入总量为942×104 t,年平均为55×104t;1999~2009年共减少泥沙输入总量为1 105×104t,年平均为111×104t。资水干流的柘溪水库库容为35.65×108m3,占多年平均径流的15.6%,于1958年开始建设,1961年2月开始蓄水,1962年开始发电,这一工程的建设是导致桃江站双累积曲线最明显变化的最主要原因。在资水流域,陆续建成的中型水库有40座,总库容达19×108m3,占多年平均径流的8.3%,大量中型水库的综合作用是1982和1998年转折点形成主要原因。
沅水桃源站年径流量和年输沙量的双累积曲线有两个转折点,即1979和1997年（图4c）。经计算,1980~1997年共减少泥沙输入总量为11 175×104t,年平均为 621×104t;1998~2009年共减少泥沙输入总量为7 439×104t,年平均为620×104t。沅江支流的凤滩水库（库容为16.76×108m3,占多年平均径流的2.6%）于1970年开始建设,于1978年建成蓄水发电,这一工程是1979年转折点形成的主要原因。此外,沅江干流的五强溪水库（库容为42.9×108m3,占多年平均径流的6.7%）于1989年开始建设,于1994年建成蓄水发电,这一水库的建成是1997年转折点形成的重要原因。同时,1998年以后陆续建成的凌津滩水库（库容为6.34×108m3）、碗米坡水库（库容为3.78×108m3）、洪江水库（库容为3.2×108m3）和蟒塘溪水库（库容为1.53×108m3）对输沙量的下降也有一定作用,体现了水库群的综合减沙效应,对1997年转折点的形成有一定的推动作用。
图4 1951~2009年湖南四水年径流量和年输沙量双累积曲线变化
Fig.4 Double mass analysis between annual runoff and sediment discharge of Four Rivers of Hunan Province in 1951-2009
澧水石门站年径流量和年输沙量的双累积曲线仅在1996年形成转折点（图4d）。经计算,1997~2009年共减少泥沙输入总量为4 616×104t,年平均为355×104t。澧水支流的江垭水库（库容为17.41×108m3,占多年平均径流的11.9 %）于1993年开始建设,于1997年建成,是形成这一转折点的主要原因。
此外,1980 s以来开展的植树造林活动使湖南省的森林面积不断增加,森林覆盖率由1983年的40.2%提高为2007年的51.6%[19],森林的大面积增加有助于提高四水流域的水土保持功能,减少水土流失,降低四水的输沙量。因此,径流量偏小、森林面积的增加以及大中型水库的综合作用可能是导致2001~2009年湖南四水入洞庭湖年输沙量下降的主要原因。
1） 59 a来湖南四水入洞庭湖年径流量存在多个上升-下降过程,但无明显的增加或减少趋势。虽然四水年输沙量也存在多个上升-下降过程,但总体上呈下降趋势。
2） 年输沙量突变分析表明,59 a来湘、资、沅、澧四水的突变时间分别为1996、1973~1974、1997~1998、1998年,突破95%置信区间的时间分别为1999、1963、1988、2001年,其中资水和沅水的突变趋势比湘水和澧水的更强。
3） 双累积曲线法消除径流量的影响后表明,大中型水库的建设是年输沙量减少的主要原因。此外,水土保持工程的实施、大中型水库建设以及径流量的减少等综合作用是2001年后输沙量减少的主要原因。
The authors have declared that no competing interests exist.
| [1] |
洞庭湖生态环境的演变、问题及保护措施 [J]. |
| [2] |
The hydrological effect between Jingjiang River and Dongting Lake during the initial period of Three Gorges Project operation [J]. |
| [3] |
南水北调细线一期工程调水区径流量与影响因子关系分析——以达曲为例 [J]. |
| [4] |
Sediment concentration changes in runoff pathways from a forest road network and the resultant spatial pattern of catchment connectivity [J]. |
| [5] |
塔里木河源流区气候变化和年径流量关系初探 [J]. |
| [6] |
近45年来降水变化和人类活动对潮河流域年径流量的影响 [J].
|
| [7] |
土地利用/覆被变化下的暴雨径流过程模拟分析——以太湖上游西苕溪流域为例 [J]. |
| [8] |
近60年来汾河入黄河水沙演变特征及驱动因素 [J]. |
| [9] |
黑河流域近60年来径流量变化及影响因素 [J]. |
| [10] |
Sediment and runoff changes in the Yangtze River basin during past 50 years [J]. |
| [11] |
鄱阳湖流域水沙时空演变特征及其机理 [J]. |
| [12] |
长江上游重点产沙区产沙量对人类活动的响应 [J]. |
| [13] |
|
| [14] |
气候突变的定义和检验方法 [J].
|
| [15] |
近数十年湘江流域河流水文变化规律 [J]. |
| [16] |
1960~2005年湖南省降水的变化 [J]. |
| [17] |
|
| [18] |
湖南省大型水库特点及防洪格局 [J]. |
| [19] |
湖南省农业洪旱灾害演变趋势和成因及防灾减灾措施 [J]. |
/
| 〈 |
|
〉 |
