湖泊科学   2021, Vol. 33 Issue (1): 289-298.  DOI: 10.18307/2021.0121.
0

研究论文——水文与水资源保护

引用本文 [复制中英文]

贲鹏, 虞邦义, 张辉, 胡勇, 洪泽湖水沙变化趋势和冲淤时空分布及驱动因素. 湖泊科学, 2021, 33(1): 289-298. DOI: 10.18307/2021.0121.
[复制中文]
Ben Peng, Yu Bangyi, Zhang Hui, Hu Yong. Spatiotemporal runoff and sediment variation, deposition-erosion characteristics and their driving factors in Lake Hongze. Journal of Lake Sciences, 2021, 33(1): 289-298. DOI: 10.18307/2021.0121.
[复制英文]

基金项目

国家重点研发计划项目(2017YFC0405602)和国家自然科学基金项目(51809002)联合资助

通信作者

贲鹏, E-mail:benpeng11172219@126.com

文章历史

2020-03-10 收稿
2020-06-15 收修改稿

码上扫一扫

洪泽湖水沙变化趋势和冲淤时空分布及驱动因素
贲鹏1,2 , 虞邦义1,2 , 张辉1,2 , 胡勇1,2     
(1: 安徽省(水利部淮河水利委员会)水利科学研究院, 蚌埠 233000)
(2: 水利水资源安徽省重点实验室, 蚌埠 233000)
摘要:为了研究洪泽湖水沙特性、变化趋势与冲淤时空分布规律,运用累积距平法、Mann-Kendall趋势与突变检验以及R/S分析法等方法,分析了洪泽湖1950—2016年的水沙特征;采用地理信息技术,基于1992年和2016年实测地形,对湖区泥沙冲淤空间分布进行了定量计算与分析.结果表明,入湖径流量无明显增加或减少的趋势,输沙量和含沙量呈明显减小趋势,1990年以后含沙量基本稳定在0.2 kg/m3以下;淮河干流(包括池河)入湖水量和沙量约占入湖总量的89.6%,三河闸出湖水沙占总出湖量的60%.淮干入湖口和溧河洼为主要淤积区域,淤积量分别为2300×104和1900×104 m3,平均淤积厚度分别为0.35和0.25 m;其他区域自然冲淤基本平衡.上游水库和河道闸坝的拦沙作用、农业种植结构变化和水土保持、大规模人工采砂等是入湖沙量减少的主要影响因素;湖区水动力特性是泥沙自然淤积主导因素,而湖区库容变化的主因则是人工采砂、围湖造田和围网养殖,且人类活动的影响远大于自然冲淤.
关键词洪泽湖    水沙变化    湖盆演变    冲淤分布    驱动因素    
Spatiotemporal runoff and sediment variation, deposition-erosion characteristics and their driving factors in Lake Hongze
Ben Peng1,2 , Yu Bangyi1,2 , Zhang Hui1,2 , Hu Yong1,2     
(1: Anhui & Huaihe River Institute of Hydraulic Research, Bengbu 233000, P. R. China)
(2: Key Laboratory of Water Conservancy and Water Resources of Anhui Province, Bengbu 233000, P. R. China)
Abstract: In order to study the water and sediment variation trend, changes in sediment deposition and erosion and spatial distribution of Lake Hongze, the methods of cumulative anomaly, Mann-Kendall trend and mutation test, and R/S analysis were used to analyze the water and sediment characteristics of Lake Hongze from 1950 to 2016. Using geographic information technology, based on the terrain of the lake area in 1992 and 2016, the spatial distribution of sediment erosion and siltation was quantitatively calculated and analyzed. The results show that, on the premise of no increasing or decreasing of the flow rate, the sediment transport and concentration show an obvious decreasing trend. After 1990, the sediment concentration has stabilized below 0.2 kg/m3; The amount of water and sediment flowing into the mainstream of the Huaihe River (including the Chihe River) accounts for about 89.6% of the total amount of water flowing into the lake, and the outflow of Sanhe Sluice accounts for 60% of the total outflow. The main sedimentary areas of the Lake Huaihe inlet and Lihewa are the sedimentation volumes of 2300×104 and 1900×104 m3, and the average sedimentation thickness is 0.35 and 0.25 m, respectively; Basic balance of natural erosion and deposition in other areas. The effects of the upstream reservoir sediment blocking, agricultural planting changing underlying surface, and artificial sand mining are the key reasons for the continuous decline of sediment concentration entering the lake. The hydrodynamic characteristics of the lake area are the dominant factors for the natural sedimentation of the sediment, while the main reasons for the change in the storage capacity of the lake area are artificial sand mining, lake farming, and seine culture, and the impact of human activities is far greater than that of natural erosion and deposition.
Keywords: Lake Hongze    water and sediment variation trend    lake basin evolution    erosion and deposition distribution    driving factors    

洪泽湖位于淮河中游和下游的结合部,承泄淮河上、中游15.8×104 km2流域面积的来水来沙,是一座具有防洪、蓄水、水产、航运和生态等综合利用功能的大型平原水库,也是南水北调东线工程的重要输水通道.洪泽湖属浅水湖,湖盆呈浅碟形,岸坡平缓,湖底较平坦,湖底高程一般在10~11 m之间.

洪泽湖是淮河行洪、蓄洪的重要载体,淮河中游洪涝与洪泽湖的形成与扩大有直接的联系.淮河干流河道与洪泽湖是一对矛盾的统一体,表现出复杂的互馈关系,淮河干流水沙条件决定了洪泽湖冲淤演变规律,而洪泽湖是淮河干流的侵蚀基准面,较高的湖区水位导致淮河洪水下泄不畅,成为淮河严重洪涝灾害的主要原因之一.洪泽湖来水来沙变化趋势和冲淤时空分布规律研究主要涉及以下问题,在自然和人类活动影响下,洪泽湖水沙和冲淤如何变化?如何量化这些变化?是否具有趋势性?未来的情势如何?变化的原因是什么?这些变化与河床和湖盆演变关系如何?维持河床与湖盆形态稳定性条件是什么?这些的基础科学研究对新时期解决淮河中游洪涝灾害,探讨淮河与洪泽湖的关系、湖区综合治理、水资源保护利用等关键问题的意义重大.

目前,关于湖泊水沙变化规律及冲淤特性研究较多,多集中在鄱阳湖和洞庭湖,而洪泽湖相关研究较少.代稳等[1]对洞庭湖降水量、径流量和输沙量之间的相互关系进行分析,确定影响水沙变化的驱动因素;罗蔚等[2]分析了鄱阳湖水沙变化特征并初步探讨了水沙变化规律;邓吉秋等[3]基于GIS对目平湖泥沙冲淤变化进行空间定量分析;高俊峰等[4]研究表明,洞庭湖近24 a来总的趋势是淤积的,局部有冲刷,但总体上淤积量大于冲刷量.针对淮河与洪泽湖冲淤演变研究主要集中在降雨径流变化特性和沙量平衡冲淤计算方面,孙鹏等[5]分析了淮河流域中上游径流年内分配、年际变化、径流趋势、突变特征及周期变化等径流过程变化特征;刘源等[6]分析了洪泽湖入湖水沙的周期性变化;虞邦义等[7]基于实测水沙资料,得到洪泽湖区泥沙冲淤情况.目前,基于洪泽湖长系列水沙资料的径流、泥沙变化规律探究较少,也尚未见运用GIS技术对洪泽湖泥沙冲淤时空分布进行定量研究.

本文采用长序列水沙资料和最新地形数据,通过累积距平法、Mann-Kendall(M-K)、R/S分析法等解析洪泽湖水沙的组成、时空演变规律,并预测其变化趋势;采用GIS数据处理方法分析洪泽湖泥沙冲淤量和空间分布以及入湖水沙变化和冲淤分布的主要驱动因素,揭示自然因素和人类活动耦合作用下洪泽湖冲淤格局的时空演变规律及响应关系.

1 水沙时空分布规律 1.1 河流水系情况

洪泽湖入湖河流主要有淮河干流、池河、怀洪新河、新汴河、濉河、老濉河和徐洪河等,大多分布于湖区的西北部;出湖河道主要有入江水道、入海水道、分淮入沂、苏北灌溉总渠,均位于湖区东南部.洪泽湖水系概化见图 1.

图 1 洪泽湖水系概化图 Fig.1 Water system generalization of Lake Hongze
1.2 入湖水沙历史变化趋势

淮河干流是洪泽湖最大的来水来沙河道,吴家渡水文站是距离洪泽湖最近的长系列水文资料观测站点,其水沙系列可以反映洪泽湖入湖水沙特性及变化规律.由图 2可知,吴家渡站近66 a来,多年平均径流量为266.10×108 m3,各年代平均径流量整体上呈现交替状态,并无明显增加或减少的趋势;多年平均输沙量为827.02×104 t,各年代平均输沙量随时间推移均呈明显减小趋势,以22.882×104 t/a的速率减少,且输沙量在1985年后整体上均低于近66 a平均值;多年平均含沙量为0.28 kg/m3,各年代平均含沙量随时间推移均呈明显减小趋势,以0.0074 kg/(m3·a)的速率减少,1990年以后含沙量基本稳定在0.2 kg/m3以下.

图 2 吴家渡站径流量(a)、输沙量(b)和含沙量(c)的逐年变化 Fig.2 Changes in the annual runoff (a), sediment discharge (b) and sediment concentration (c) at the Wujiadu Station
1.3 入湖泥沙变异性 1.3.1 分析方法

1) 累积距平法

累积距平法是由曲线直观判断离散数据点变化趋势的一种非线性统计方法,根据曲线上下起伏特性,可以判断序列要素长期显著的变化趋势及持续性变化,同时也还可以判断出发生突变的大致时间[8].

对于某一序列x,其某一时刻t的累积距平表示为:

$ {\hat X_t} = \mathop \sum \limits_{i = 1}^t \left( {{X_i} - \bar X} \right){\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {t = 1, 2, 3, \cdots , n} \right) $ (1)

式中,X为序列的平均值,$\bar X = 1/t\left( {\mathop \sum \limits_{i = 1}^t {X_i}} \right)$t为时间.

2) M-K分析法

M-K分析法适用于不符合正态分布的水沙数据,不受少数异常值的干扰,计算简便[9].

对于具有ri个样本量的时间序列x,构造一秩序列:

$ {S_k} = \mathop \sum \limits_{i = 1}^k {r_i}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {k = 2, 3, 4, \cdots , n} \right) $ (2)

式中,Sk为第i时刻数值大于第j时刻数值个数的累计数,ri为第i个样本xi大于yi(1≤ji)的累计数.

Sk标准化得到如下统计量:

$ \begin{array}{l} U{F_k} = \left( {{S_k} - E\left( {{S_k}} \right)} \right)/\sqrt {{\mathop{\rm var}} \left( {{S_k}} \right)} \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {k = 1, 2, 3, 4, \cdots , n} \right) \end{array} $ (3)

式中,UF1=0,E(Sk)和var(Sk)分别为Sk的均值和方差.给定显著性水平α,若$\left| {U{F_k}} \right| > {U_{a/2}}$, 则表明时间序列存在明显的变化趋势.

将时间序列x按逆序重复上述步骤,同时使:

$ U{B_k} = - U{F_k}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {k = n, n - 1, n - 2, \cdots , 1} \right) $ (4)

通过对UBkUFk曲线的分析,得到突变时间与突变区域.

1.3.2 变异性规律

采用累积距平法对吴家渡站年均来沙量进行多沙少沙变化规律分析,年输沙量与年含沙量多沙少沙变化过程基本相同,1980s前为多沙期,随后转入少沙期(图 3).采用M-K法对来沙序列进行突变检验,来沙量的UFk值整体都为负值,表明来沙量一直处于减少的趋势,减少趋势显著,突变发生在1980s(图 4).

图 3 来沙量累积距平曲线 Fig.3 Cumulative anomaly curve of incoming
图 4 来沙量M-K法统计量曲线 Fig.4 Statistic curve of M-K method in sand volume
1.4 入湖泥沙未来变化趋势 1.4.1 分析方法

R/S分析法最先是由英国水文家Hurst提出的一种基于长程相关思想的时间序列统计方法,已发展为一种广泛应用的趋势分析方法[10],其基本原理如下:

考虑一个时间序列{ξ(t)},t=1,2,….对于任意正整数τ≥1,定义均值序列为:

$ {\left\langle \xi \right\rangle _t} = \frac{1}{\tau }\mathop \sum \limits_{i = 1}^\tau \xi \left( t \right){\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {\tau = 1, 2, 3, 4, \cdots , n} \right) $ (5)

累积离差X(tτ)为:

$ X\left( {t, \tau } \right) = \mathop \sum \limits_{u = 1}^t \left( {\xi \left( u \right) - {{\left\langle \xi \right\rangle }_t}} \right){\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {1 \le t \le \tau } \right) $ (6)

极差R(τ)为:

$ R\left( \tau \right) = {\max _{1 \le t \le \tau }}X\left( {t, \tau } \right) - {\min _{1 \le t \le \tau }}X\left( {t, \tau } \right){\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {\tau = 1, 2, 3, 4, \cdots , n} \right) $ (7)

标准差S(τ)为:

$ S\left( \tau \right) = {\left\{ {\frac{1}{\tau }\mathop \sum \limits_{i = 1}^\tau {{\left[ {\mathop \sum \limits_{u = 1}^t \left( {\xi \left( t \right) - {{\left\langle \xi \right\rangle }_t}} \right.} \right]}^2}} \right\}^{1/2}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {\tau = 1, 2, 3, 4, \cdots , n} \right) $ (8)

研究发现:

$ R\left( \tau \right)/S\left( \tau \right) = K{\left( {\tau /2} \right)^H} $ (9)

式中,K为比例系数,H为Hurst指数.

1.4.2 变化趋势分析

采用R/S分析法对吴家渡站未来一段时间内来沙趋势进行判断.由图 5可知,来沙量时间序列R/S法统计量H回归拟合程度较好,输沙量、含沙量与历时的相关系数R2分别可达到0.9645和0.9628.输沙量和含沙量时间序列的Hurst指数均大于0.5,表明未来一段时间输沙过程具有持续性,且与历时有正相关性,即未来变化趋势将和输沙量历史总体趋势保持相同.

图 5 吴家渡站来沙量R/S法统计量曲线 Fig.5 R/S statistic curve of incoming sand volume at the Wujiadu Station
1.5 入出湖水沙空间分布

淮河干流小柳巷站是距离洪泽湖最近的水文站,于1982年建站,能够更好地反映近期淮河干流入洪泽湖水沙量,以其近34 a的水沙系列分析入、出湖水沙空间分布, 对于个别年份缺测资料进行了相关插补[13].由表 1可知,1983-2016年淮河(淮河干流和池河)来水来沙量均占入湖总量的89.6 %,平均含沙量为0.19 kg/m3;怀洪新河来水量占5.7 %,来沙量仅占1.5 %,平均含沙量为0.05 kg/m3,怀洪新河下游地区地势低洼,湖泊、河道众多,泥沙沉积后入湖量较小;溧河洼西北支流来水量占3.8 %,来沙量占6.5 %,平均含沙量为0.33 kg/m3,含沙量较大.由表 2可知,三河闸是水沙出湖的主要通道,约60 %水沙量经入江水道排入长江,未考虑蒸发和取用水.洪泽湖水沙同源,输沙量的大小受控于来水量大小,符合大水大沙的特点.

表 1 1983-2016年洪泽湖入湖水沙空间分布 Tab. 1 Spatial distribution of water and sediment in Lake Hongze from 1983 to 2016
表 2 1983-2016年洪泽湖出湖水沙空间分布 Tab. 2 Spatial distribution of water and sediment out of Lake Hongze from 1983 to 2016
2 输沙平衡冲淤分析 2.1 近期冲淤年际分布规律

根据洪泽湖入湖和出湖控制站计算出湖区年际泥沙冲淤量,冲淤年际分布见图 6.洪泽湖1983-2016年总体趋势是淤积,总来沙量为19726.8×104 t,出沙量为9522.5×104 t,淤积量为10204.8×104 t,淤积率为51.7 %,多年平均淤积量为300.1×104 t,其中2007年淤积量最大,为1044.4×104 t.系列资料34 a中2001、2009、2011、2012年发生微冲,冲刷量都小于50×104 t;1994和2014年输沙基本平衡;1992、1993、1995、1997和1999年共5 a发生微淤,淤积量小于100×104 t,其余23 a洪泽湖都发生淤积,淤积量大于100×104 t.

图 6 洪泽湖年际冲淤分布 Fig.6 Interannual scouring and deposition of Lake Hongze
2.2 长系列冲淤量及变化趋势

文献[7]推算了1950-1959年、1966-1976两个时段的洪泽湖冲淤量;文献[12]计算了1960-1965年、1972-1982年两个时段的洪泽湖冲淤量;结合近期水沙资料,进一步推算出1950-2015年洪泽湖总淤积量为30174×104 t,年均淤积457×104 t.从表 3可以看出,洪泽湖淤积量呈减少趋势,洪泽湖淤积量减少与淮河干支流来沙量变化趋势一致.

表 3 洪泽湖各年代平均淤积量 Tab. 3 Average sedimentation of Lake Hongze in each age
3 湖盆冲淤分布规律

以1992年(1 :10000地形图,1985国家高程)和2016年(1 :10000地形图,1985国家高程)实测地形图为研究的基础数据,基于GIS的泥沙冲淤变化空间定量分析和地理空间数据,运用数字化、GIS的数据处理和空间分析功能,实现对不同时间泥沙淤积状态的计算与对比,从而对泥沙冲淤的空间分布变化进行计算与分析.

3.1 冲淤高程分布特征

本次仅分析洪泽湖湖区13.5 m高程以下湖盆的冲淤变化,不包含周边滞洪圩区及湖泊洼地.湖区1992年高程见图 7a,2016年高程见图 7b,1992和2016年库容变化见表 4,可以看出,整个湖区库容增加0.91×108 m3,包括采砂量和航道疏浚量.9.5 m高程以下湖区变化最大,库容增加0.92×108 m3,主要变现为局部深坑;10.5 m高程以下库容增加1.29×108m3,9.5~10.5 m高程区域库容增加0.37×108 m3,主要变现为局部深坑和航道变化;11.5 m高程以下库容增加0.71×108 m3,10.5~11.5 m高程区域库容减少0.58×108 m3,说明该区域发生了淤积;11.5~13.5 m之间湖区库容变化较小,冲淤幅度较小.

图 7 洪泽湖湖区1992年(a)和2016年(b)高程 Fig.7 Elevations of Lake Hongze in 1992 (a) and 2016 (b)
表 4 各高程湖区库容及变化量(×108 m3) Tab. 4 Storage capacity and erosion and variation of each elevation lake area (×108 m3)

2016年湖底出现较多深坑,局部区域下切近6 m,初步分析,湖区深坑不是自然冲淤的结果,而是人工采砂造成的[13].为了反映湖区真实冲淤量,需要对局部DEM地形进行修正,扣除明显的采砂坑和航道疏浚量后,湖区表现为淤积,淤积量为4200×104 m3,平均淤积厚度为0.03 m.航道整治方量修正计算方法,以2016年地形航道边界为基准,将航道区域冲淤量设置为0 m3;人工采砂量周边修正计算方法,将采砂坑周边局部静水或流速较小且冲刷厚度大于0.3 m区域的冲刷量设置为0 m3.

3.2 冲淤空间分布特征

湖区不同区域的形态和水流特性相差较大,冲淤变化亦不同,根据入、出流位置和湖区形态将湖盆分为5个区域,分别为淮干入湖口区、溧河洼区、成子湖区、三河闸区和中心湖区(图 8).淤积区域主要集中在淮干入湖口、溧河洼以及临淮镇和半城镇局部区域;成子湖和三河闸区冲淤幅度较小,局部回淤;中心湖区冲淤总体表现为冲刷,冲刷量较小.淮干入湖口区淤积量约为2300×104 m3,局部淤积厚度达1.0 m,平均淤积厚度为0.35 m.溧河洼主要淤积在怀洪新河入湖口以上湖区,淤积量约为1900×104 m3,平均淤积厚度为0.25 m;怀洪新河入湖口以下湖区未淤积,有微冲.临淮镇和半城镇等局部区域淤积约为900×104 m3.中心湖区总体表现为冲刷,西部冲淤基本平衡,东部湖区湖底下切,下切深度为0.1~0.2 m,库容约增加1600×104 m3,由于区域内采砂严重,库容增加并非自然冲刷的结果.

图 8 湖区冲淤空间分布(+为淤积;-为冲刷) Fig.8 Spatial distribution of erosion and deposition in the lake area (+ : deposition; -: scour)
3.3 地形法与输沙法对比

按地形法计算,1992-2016年间,扣除人工采砂和航道疏浚量后,湖区表现为淤积,淤积区域为淮干入湖口和溧河洼区西北部,淤积量为4200×104 m3,这两个区域是湖区沙源所在.按输沙法计算,1992 2016年间,湖区共淤积5600×104 t,相当于淤积量为3700×104 m3.两种计算方法结果基本吻合,也可以说明洪泽湖来沙主要淤积在淮干入湖口和溧河洼区域,其他湖区冲淤基本平衡.

4 水沙变化和湖盆演变驱动因素 4.1 水沙变化趋势原因分析

出湖水沙特性受入湖水沙、降雨、蒸发、湖区冲淤、水资源利用等因素影响,其核心因素是入湖水沙条件,这里仅分析入湖水沙变化原因.

1) 洪泽湖入湖水量变化的主导因素为天然降雨.1950s以来,洪泽湖流域面雨量呈波动增加趋势,但是总体变化趋势不明显[11, 14-15],流域用水总量虽然呈现缓慢增长趋势,但增长缓慢,并且增加量主要源于地下水[18].面雨量和用水量均没有发生明显变化,洪泽湖入湖水量也无明显变化趋势.

2) 流域大型水库和河道闸坝的蓄水拦沙是入湖沙量减少的重要原因[16].淮河上游山区和淮北平原水系是淮河流域的主要产流产沙区,1970s建成了五岳、泼河和鲇鱼山水库,1987年板桥水库复建,1993年石漫滩水库复建,随着水库库容的增加,拦沙作用更为明显,使得上游泥沙淤积在水库中,淮河中游河道含沙量减小. 1990s,淮河水系水库和闸坝建设基本完成,区域来沙也基本稳定下来.

3) 农业种植结构改变和水土保持改变了流域下垫面[17].淮北平原是我国重要的粮食生产基地,淮北支流洪汝河、沙颍河来沙量对淮河干流影响较大.1980s开始,为了增加粮食产量,采用精耕细作、涵养水分和养分、增加复种指数等措施,这些技术措施正是水土保持的第一道防线.同时,种植结构也发生一些变化,如旱地作物向水田作物转变、粮食作物向经济作物转变等,地面裸土面积减小,改变了淮北地区的下垫面.同时,治淮骨干工程要求进一步加大水土保持工程建设,水土流失得到了有效的遏制,使得淮河干流含沙量继续减小.

4) 人工采砂直接影响入湖沙量,而且改变河床形态.1980s中后期,淮河干流及主要支流河道开始大规模采砂,根据淮干历次地形及历史水沙资料计算分析,仅淮河干流王家坝至浮山段河道人工采砂量约为2×108 m3,无序采砂导致河床出现大量不连续的深坑,河道成为沉砂池.

4.2 湖盆演变原因分析

2016年较1992年洪泽湖13.5 m高程以下库容增加9100×104 m3,而大部分区域表现为淤积或者冲淤基本平衡,仅局部区域出现深坑,湖区不同区域冲淤特性差别较大,导致这种现象的主要原因如下:

1) 湖区水动力特性是泥沙自然淤积分布的重要影响因素[18].湖区泥沙淤积区域为淮河干流入湖口和溧河洼北部,是湖区最大的两个来水来沙区.淮河干流入湖口水面突然开阔,入湖河段河底倒比降,洪水期河道段流速为1~2 m/s,进入湖区扩散段后流速降至0.5 m/s以下,水动力条件减弱,挟沙能力明显下降,导致河口泥沙累积性淤积,特别是入湖河口淤积浅滩发育.

2) 围湖造田和围网养殖是洪泽湖局部淤积的重要原因.洪泽湖的圈圩垦殖和围网养殖始于1960s,在高程12.5 m以下的滩地进行圈地种植,主要分布在溧河洼湖区、淮干入湖口、临淮镇和半城镇等西北湖区,这些区域正是湖区的淤积位置.圈圩垦殖和围网养殖不仅占用洪泽湖库容,而且严重阻水,成为泥沙的沉积区,这也是淮干入湖口和溧河洼湖区泥沙淤积的重要原因.1980s以来,洪泽湖西岸临淮镇附近的湖泊变化最大,大部分湖泊己经被变成了围网养殖区[19],从航拍图片(图 9)来看,近岸的湖面被隔堤和渔网切割成一个个小块,向内湖延伸近千米.

图 9 溧河洼围网养殖航拍图 Fig.9 Aerial photography of seine culture in Lihewa

3) 人工采砂是影响湖区库容的主要因素之一.中心湖区东部出现大量深坑,且湖底出现大面积下切,此区域并非主流区,不会出现大范围冲刷,原因是人工采砂[13].2005年,洪泽湖开始出现采砂;2006年,采砂区域向洪泽湖水域纵深发展;2010年,采砂区域扩大至泗洪、泗阳等多水域;2015年,采砂区域扩大整个湖区,采砂船已增加到600余条(图 10).随着采砂区域不断扩大和采砂船数量不断增加,大规模采砂活动对湖区库容产生较大影响.1992-2016年,洪泽湖采砂量超1×108 m3,除采砂区域,期间洪泽湖表现为淤积或冲淤基本平衡状态,人类活动的影响远大于自然冲淤.

图 10 2014年洪泽湖局部水域采砂航拍图 Fig.10 Aerial photo of sand mining in Lake Hongze in 2014
5 结语与探讨

1) 基于洪泽湖1950-2016年的实测水沙资料以及1992年和2016年地形资料,对其水沙变化趋势和湖区冲淤分布进行了定量计算与定性分析.结果表明,洪泽湖入湖水量无明显增加或减少的趋势,输沙量和含沙量呈明显减小趋势,淮河干流(包括池河)入湖水量和沙量约占入湖总量的89.6 %;(2)淮河干流入湖口和溧河洼上部为主要淤积区域,平均淤积厚度0.35和0.25 m,其他区域自然冲淤基本平衡;(3)水库和道闸坝的拦沙作用,农业种植结构变化,水土保持、大规模采砂等因素导致入湖沙量减少;湖区水动力特性是泥沙自然淤积的主导因素,而湖区库容变化的主因则是人工采砂、围湖造田和围网养殖,且人类活动的影响远大于自然冲淤.

2) 1992-2016年间,在自然和人类活动作用下,淮干入湖口和溧河洼区域出现淤积,但淤积厚度不大;湖区冲淤基本平衡,淮干入湖口至三河闸、淮干入湖口至二河闸和怀洪新河入湖口至二河闸等以及湖区主要流路也没有出现淤积现象.根据入湖水沙变化规律,在流域下垫面不发生较大变化条件下,入湖含沙量将维持在较低水平,淤积量也将处于较低水平.洪泽湖综合整治涉及湖区疏浚,河床稳定性以及入海水道、入江水道、冯铁营引河和老河道如何衔接等问题,水沙变化趋势和冲淤时空分布规律可以为湖区综合治理提供参考.

6 参考文献

[1]
Dai W, Lv DQ, Li JB et al. Evolution characteristics and driving factors of runoff and sediment changes in Dongting Lake During 1951-2014. Journal of Soil and Water Conservation,, 2017, 31(2): 142-150, 176. [代稳, 吕殿青, 李景保等. 1951-2014年洞庭湖水沙阶段性演变特征及驱动因素分析. 水土保持学报, 2017, 31(2): 142-150, 176.]
[2]
Luo W, Zhang X, Deng ZM et al. Runoff and sediment load transport and cause analysis in Poyang Lake basin over the period 1956-2008. Advances in Water Science, 2014, 25(5): 658-667. [罗蔚, 张翔, 邓志民等. 1956-2008年鄱阳湖流域水沙输移趋势及成因分析. 水科学进展, 2014, 25(5): 658-667.]
[3]
Deng JQ. GIS-based spatial quantitative analysis on change of sediment deposition and erosion in Muping Lake. Journal of Central South University:Science and Technology, 2007(5): 1000-1006. [邓吉秋. 基于GIS的目平湖泥沙冲淤变化空间定量分析. 中南大学学报:自然科学版, 2007(5): 1000-1006. DOI:10.3969/j.issn.1672-7207.2007.05.036]
[4]
Gao JF, Zhang C, Jiang JH et al. Changes in sediment deposition and erosion and their spatial distribution in the Dongting Lake. Acta Geographica Sinica, 2001(3): 269-277. [高俊峰, 张琛, 姜加虎等. 洞庭湖的冲淤变化和空间分布. 地理学报, 2001(3): 269-277. DOI:10.3321/j.issn:0375-5444.2001.03.003]
[5]
Sun P, Sun YY, Zhang Q et al. Temporal and spatial variation characteristics of runoff processes and its causes in Huaihe Basin. J Lake Sci, 2018, 30(2): 497-508. [孙鹏, 孙玉燕, 张强等. 淮河流域径流过程变化时空特征及成因. 湖泊科学, 2018, 30(2): 497-508. DOI:10.18307/2018.0221]
[6]
Liu Y, Xu GB, Duan Y et al. Multi-time scale wavelet analysis of streamflow and sediment sequences into Hongze Lake. Water Resources and Hydropower Engineering, 2020, 51(2): 128-135. [刘源, 徐国宾, 段宇等. 洪泽湖入湖水沙序列的多时间尺度小波分析. 水利水电技术, 2020, 51(2): 128-135.]
[7]
Yu BY, Yu YS. Sediment deposition in Hongze Lake. Journal of Sediment Research, 2010(6): 38-43. [虞邦义, 郁玉锁. 洪泽湖泥沙淤积分析. 泥沙研究, 2010(6): 38-43.]
[8]
Liu Y, Wang Y, Yang XG. Trends in reference crop evapotranspiration and possible climatic factors in the North China Plain. Acta Ecologica Sinica, 2010, 30(4): 923-932. [刘园, 王颖, 杨晓光. 华北平原参考作物蒸散量变化特征及气候影响因素. 生态学报, 2010, 30(4): 923-932.]
[9]
Yao WY, Ran DC, Chen JN. Recent changes in runoff and sediment regimes and future projections in the Yellow River basin. Advances in Water Science, 2013, 24(5): 607-616. [姚文艺, 冉大川, 陈江南. 黄河流域近期水沙变化及其趋势预测. 水科学进展, 2013, 24(5): 607-616.]
[10]
Mandelbrot BB, Wallis JR. Some long-Run properties of geophysical records. Water Resources Research, 1969, 5(2): 321-340. DOI:10.1029/WR005i002p00321
[11]
Wang KQ, Zeng Y, Xie ZQ et al. Change trend of temperature and precipitation in Huaihe River Basin from 1961 to 2008. Scientia Meteorologica Sinica, 2012, 32(6): 671-677. [王珂清, 曾燕, 谢志清等. 1961-2008年淮河流域气温和降水变化趋势. 气象科学, 2012, 32(6): 671-677. DOI:10.3969/2012jms.0061]
[12]
Zhu SQ, Dou HS eds. Hongze Lake-Water resources and aquatic living resources. Hefei: China Science and Technology Press, 1993. [朱松泉, 窦鸿身. 洪泽湖——水资源和水生物资源. 合肥: 中国科学技术出版社, 1993.]
[13]
Yan DY. Analysis on sand mining management in Hongze Lake. Jiangsu Science & Technology Information, 2015(31): 49-50. [严登余. 洪泽湖采砂管理分析. 江苏科技信息, 2015(31): 49-50.]
[14]
Ye ZW, Li ZH. Characteristics and trend of areal precipitation in the Hongze Lake basin since 1951. Resources and Environment in the Yangtze Basin, 2010, 19(12): 1392-1396. [叶正伟, 李宗花. 1951年来洪泽湖流域面雨量变化特征与趋势分析. 长江流域资源与环境, 2010, 19(12): 1392-1396.]
[15]
Wang H, Mei M, Zhan TS et al. Current situation and prospect of water resources development and utilization in Huaihe River area. Water Resources Planning and Design, 2005(1): 10-29-31. [王浩, 梅梅, 詹同涛等. 淮河区水资源开发利用现状及展望. 水利规划与设计, 2005(1): 10-29-31.]
[16]
Shi HL, Hu CH, Wang YG et al. Variation trend and cause of runoff and sediment load variations in Huaihe River. Journal of Hydraulic Engineering, 2012, 43(5): 571-579. [史红玲, 胡春宏, 王延贵等. 淮河流域水沙变化趋势及其成因分析. 水利学报, 2012, 43(5): 571-579.]
[17]
Mao SM, Xu H. Discussion on the paper "Variation trend and cause of runoff and sediment load variations in Huaihe River". Journal of Hydraulic Engineering, 2013, 44(7): 880-882. [毛世民, 许浒. 对"淮河流域水沙变化趋势及其成因分析"一文的商榷. 水利学报, 2013, 44(7): 880-882.]
[18]
Yu BY, Ni J, Yang XJ et al. Research on hydrodynamic numerical model for main stream of Huaihe River from Fushan to outlet of Hongze Lake. Water Resources and Hydropower Engineering, 2011, 42(8): 38-42. [虞邦义, 倪晋, 杨兴菊等. 淮河干流浮山至洪泽湖出口段水动力数学模型研究. 水利水电技术, 2011, 42(8): 38-42.]
[19]
Ruan RZ, Xia S, Chen Y et al. Change of lake nearby Linhuai Town in West Bank of Hongze Lake during 1979-2006. Wetland Science, 2012, 10(3): 344-349. [阮仁宗, 夏双, 陈远等. 1979-2006年洪泽湖西岸临淮镇附近湖泊变化研究. 湿地科学, 2012, 10(3): 344-349.]