三峡工程运行后洞庭湖区水沙过程变异特征及其影响因素.docxVIP

三峡工程运行后洞庭湖区水沙过程变异特征及其影响因素 由于湖泊和河流关系的互动，湖泊水沙过程的变化一直是水保护工作者和科学家的中心问题。洞庭湖作为长江中游重要的通江型湖泊,承载着分泄荆江松滋、太平、藕池三口来水来沙,汇合湘、资、沅、澧四水,再由城陵矶汇入长江,江湖关系非常复杂(图1)。近100a来在人类活动和自然的双重作用下,洞庭湖不断萎缩,防洪形势严峻。特别是三峡工程建成运行以来,洞庭湖面临荆江三口新的来水来沙条件,江湖关系更加复杂。前人在长江与洞庭湖的江湖关系有了很多研究成果,特别是近10a来,许多学者从人类活动(围垦、三峡工程建设等)对洞庭湖的影响方面做了大量探索[8,9,10,11,12,13,14,15,16,17,18]。目前对洞庭湖的研究多集中在泥沙淤积计算、洪涝与防洪、水情、水环境及水位的影响、生态系统服务功能、湿地及水生生物等方面,系统全面研究在自然和人类活动双重影响下湖泊系统对水沙变异的响应,仍是洞庭湖研究中的一个前沿领域。 1 市场水沙变异响应与分析 系列水沙数据来源于湖南省水文水资源勘测局和《长江泥沙公报2000—2011》。文中选取洞庭湖湘、资、沅、澧四水、荆江三口和出口城陵矶等主要控制站(湘潭、桃江、桃源、石门、新江口、沙道观、弥陀寺、康家岗、管家铺、城陵矶)1951—2011年的实测水沙数据作为洞庭湖水沙过程变异分析的基础。基于以上数据,文章主要采用水文描述法、变差系数法、双累积曲线法、数学计量等方法探讨在三峡运行影响下洞庭湖水沙变异特性与规律性,分析人类活动影响下洞庭湖水沙变异响应,以及江湖水沙交互关系。变差系数可以反映一个流域径流泥沙的年内分配不均匀程度,Cv值越大,就表明各月径流量、输沙量相差越悬殊,表示年径流泥沙的年内丰枯变化剧烈,即年内分配不均匀,反之亦然。其计算公式为: 2 分析湖泊流域水沙特征 2.1 入湖泥沙与年输沙量变化规律 洞庭湖水沙主要来自荆江三口、湘资沅澧四水以及湖区区间,水沙组合复杂多变,影响因素众多。根据1951—2011年洞庭湖区实测水沙资料,三口、四水多年平均入湖径流量为2 573亿m3/a,城陵矶出湖多年平均径流量2 885亿m3/a,其中三口来水907亿m3/a,占出湖总水量的31.4%,四水来水1 666亿m3/a,占出湖总水量的57.8%,区间补给量为312亿m3/a,占出湖水量的10.8%。三口、四水多年平均输沙量为14 361万t/a,其中三口入湖沙量11 642万t/a,占入湖总沙量的81.1%,四水入湖沙量为2 720万t/a,占入湖总沙量的18.9%。城陵矶出湖多年平均输沙量3 890万t/a,湖内多年淤积沙量10 471万t/a,淤积率72.9%。 图2为洞庭湖三口、四水入湖和城陵矶出湖年径流量与年输沙量双累积曲线,从图2可以看出,入湖和出湖水沙60a以来均已发生明显的变化,洞庭湖水沙关系表现为以泥沙衰减为主要特征的阶段性变化,特别是荆江三口的入湖沙量和城陵矶出湖沙量衰减十分显著。从双累积曲线斜率分析,入湖水沙双累积曲线斜率K由第一时段的8.005 7衰减到第五时段的0.945 8,衰减率为88.2%;而出湖水沙双累积曲线斜率K由第一时段的1.824衰减到第五时段的0.721 9,衰减率为60.4%;特别要指出的是,第四阶段入湖水沙曲线斜率K由4.259锐减到0.945 8,说明2003年洞庭湖的入湖水沙与三峡投入运行关系密切,而出湖的水沙曲线走势变化并不强烈。可见,长江干流的水利工程对洞庭湖入湖、出湖水沙的影响不同步,入湖泥沙的衰减变化速率远大于出湖泥沙的衰减速率,且入湖水沙的变化受三峡运行的影响比较明显。 三峡蓄水运用前期(1994—2002年)四水总量为1 908亿m3/a,三峡蓄水后四水总量为1 555亿m3/a,衰减率为18.5%;而三口多年平均水量为907亿m3/a,三峡蓄水运用前期水量为616亿m3/a,蓄水后水量为500亿m3/a,与多年平均值相比,两个时期的水量衰减率分别为32.1%和44.9%。下荆江裁弯前后三口和四水入湖水量分别为886亿m3/a和834亿m3/a、1 729亿m3/a和1 699亿m3/a,两者裁弯前后水量变化并不大,可以看出三峡蓄水运行后,三口四水来水量均发生了变化,三口来水量的锐减应该归因于三峡水库的蓄水作用。来沙量方面,三口入湖沙量多年平均值为11 642万t/a,占入湖沙量的81.1%;下荆江裁弯前1951—1966年,三口来沙量20 236万t/a,占入湖沙量的85.3%,随着下荆江裁弯的实施,三口来沙量的比例降为77.6%,至三峡蓄水前三口入湖沙量的比例稳定在75%~81%。三峡水库蓄水后由于水库的拦沙作用,荆江泥沙量减少,三口入湖泥沙大幅下降,三口在三峡运行期(2003—2010)的来沙量仅为