1. 引言
大型水利工程的建造一方面可以充分利用流域水能资源,减轻下游防洪的压力,另一方面却改变了局部空间尺度上的水循环演变规律,对下游的水文情势造成了一定的影响 [1] [2]。近几十年来,随着经济社会发展,长江流域水资源开发利用程度逐步加大,溪洛渡、向家坝和三峡水库是长江干流控制性工程,调蓄能力巨大,对长江中下游河道及两湖水文情势造成了一定的影响 [3] [4] [5]。另外,随着近年来气候变化造成的局部地区干旱灾害严重,长江流域多次遭遇上下游同枯的流域性枯水水情,蓄水期上游梯级水库群联合蓄水可能会对长江中下游的供水产生一定的影响 [6] [7]。因此,为了充分发挥上游梯级水库群的综合利用效益,需全面了解上游梯级水库群蓄水期蓄水对长江中下游水文情势影响程度。
长江上游梯级水库群,尤其是三峡水库运行对下游河段的水文情势造成的影响备受学术界的关注。对于梯级水库群调度运行对下游水文情势影响研究主要分为三类:第一类是基于观测水文资料的研究,主要采用统计分析方法,对比三峡运行前后干流主要控制站点的水文情势变化情况,Zhang等 [8] 采用统计检验法分析了长江中下游的水沙变化情况,得出三峡工程修建使得中下游水文序列产生突变。李长春等 [9] 分析了蓄水期长江干流和两湖地区水文情势特征,认为少数年份存在蓄水期长江干流来水与两湖本身来水遭遇枯水的情况,使得三峡水库的蓄水形势更为严峻。第二类是采用Range of Variability Approach (RVA)指标分析方法,通过定义能够反映水沙时空变化的分析指标,定量评估三峡大坝对长江中下游水文情势的影响。班璇等 [10] 通过定义32个具有生态学意义的指标,评估了三峡蓄水后中游水沙时空分布特征。第三类方法为数值模拟法,通过构建水动力学模型,定量评估三峡蓄水对于下游的影响。王俊等 [11] 结合2009年三峡水库蓄水情况,揭示了三峡工程运行使得蓄水期下游水位下降了2~3 m。目前,数值模拟研究多集中在通过还原流量计算三峡工程对于下游水位影响,研究多集中在三峡工程对于长江中下游河段水文情势的影响,2014年后,溪洛渡、向家坝水库先后蓄水,长江中下游水文情势将进一步变化,三个水库联合蓄水对于中下游影响程度目前还未有研究。
本研究拟构建长江中下游一维水动力演算模型,利用水库运行资料对模型上边界宜昌站的流量进行还原,选取溪洛渡、向家坝两座水库正常蓄水时期(2014~2015年)作为分析时段,定量评估溪洛渡、向家坝、三峡水库联合蓄水对中下游干流水位的影响,并重点分析溪洛渡、向家坝的影响,为减轻上游梯级水库群蓄水对长江中下游生产、生活、生态、航运等用水的影响提供可靠的技术支撑。
2. 计算方法及验证
2.1. 模型方法
研究模拟范围包括长江干流宜昌至大通河段,支流包括清江、湘江、资水、沅江、澧水、洞庭湖区、汉江、鄱阳湖区、饶河、信江、抚河、赣江、修水、昌江等一级支流,三口松滋河、虎渡河、藕池河等分流口及其二级支流,河网水系概化见图1。干流河道长1095 km,采用2012年的实测河道断面资料,一共设置了744个断面,断面平均间距为2 km。根据河网概化结果,模型边界一共有25个,其中宜昌水文站为上边界,大通站为下边界,区间有洞庭湖四水、鄱阳湖五河、清江、汉江、鄂东等23个控制站作为区间点源入流。
采用长江中下游一维水动力学模型,分析梯级水库蓄水对长江中下游水位的影响。其数学表达式为。
式中:x、t分别为计算点空间和时间的坐标,A为过水断面面积,Q为过流流量,h为水位,q为旁侧入流流量,R为水力半径,α为动量校正系数,g为重力加速度。
2.2. 模型率定
利用长江中下游干支流河道断面数据及主要控制站2003~2015年实测水位、流量资料,采用2003~2007年的水文资料进行率定,根据各河段的河道形态及水位流量实测资料,针对不同水位分三段对曼宁糙率系数进行率定 [12] [13],例如对于宜昌水文站,水位低于40 m时曼宁系数n为0.028,在40 m与49 m之间时曼宁系数n为0.025,高于49 m时曼宁系数n为0.020,其他河段的率定结果见表1。
![](Images/Table_Tmp.jpg)
Table 1. Manning roughness coefficients in the coupled model
表1. 曼宁糙率系数率定成果表
2.3. 模型验证
选取2003~2007年作为模型参数率定期,2008~2015年作为模型参数的验证期。根据率定期选择的参数,模拟2008~2015年干流各个站点的逐日流量过程,并与实测流量进行对比,结果见图2。由图可以看出,各个站点模拟流量过程与实测过程拟合程度较好。干流各站的Nash系数和平均相对误差统计结果见表2,可以看出,各站的Nash系数均在0.91以上,平均相对误差均在0.1%之内,拟合程度较好。根据表2模型模拟效果评价指标等级,所有站点Nash系数和平均相对误差都达到优秀级别。本模型拟合效果较好,可用于梯级水库群调度对下游水文情势的影响的评估 [14]。
![](//html.hanspub.org/file/8-2410979x12_hanspub.png)
Figure 1. Conceptual sketch of network of waterway in the middle and lower reaches of the Yangtze River
图1. 长江中下游河网水系概化示意图
![](Images/Table_Tmp.jpg)
Table 2. Checkout of water level fitting degree of each station in the middle and lower reaches of the Yangtze River
表2. 长江中下游各站水位拟合程度检验
3. 溪洛渡和向家坝水库蓄水对三峡入库流量影响
根据溪洛渡和向家坝水库2014~2015年的实际运行资料,采用水量平衡法分别还原得到各水库的6h入库流量,考虑水流传播时间,采用马斯京根演算方法将入库流量演算到下游控制点与区间水量叠加,推求得到屏山站2014~2015年的还原后6 h流量过程。将屏山站还原后的流量过程演算到三峡水库入库点清溪场,与屏山~清溪场区间水量过程叠加,得到考虑溪洛渡、向家坝水库还原的三峡天然入库流量过程。
将还原得到的三峡水库的入库流量与实测的三峡入库流量进行对比,如表3。由表可以看出,溪洛渡、向家坝9月份蓄水一定程度上减少了三峡水库的入库流量。2014年9月,两座水库蓄水使得三峡入库各旬平均流量均减少1300 m3/s,2015年9月两座水库蓄水使得三峡入库各旬平均流量分别减少1800 m3/s、1500 m3/s以及1700 m3/s。10月份,两座水库对三峡入库流量影响较小,仅2015年下旬三峡入库的旬平均流量增加1000 m3/s,其他时间段旬平均流量变化不明显。11月份,两座水库运行一定程度增加了三峡水库的入库流量,表现为2014年三峡11月上、中旬平均流量均增加了400 m3/s,2015年三峡11月上、中旬流量分别增加1100 m3/s和600 m3/s,而下旬流量减少500 m3/s。
![](Images/Table_Tmp.jpg)
Table 3. Analysis of impounding impact of the Xiluodu and Xiangjiaba reservoirs on TGR (unit: m3/s)
表3. 溪洛渡和向家坝水库蓄水对三峡入库流量影响分析表(单位:m3/s)
注:① 表示考虑上游溪洛渡和向家坝水库蓄水影响的三峡入库流量;② 未考虑上游溪洛渡和向家坝水库蓄水影响的三峡入库流量。
4. 水库联合蓄水对干流水位的影响
4.1. 研究工况
本次共设置3种情景,各情景区别主要是模型上边界宜昌站的流量过程,其他边界条件不变,得到3种情景下的中下游各站水文过程。将各站考虑水库联合蓄水的还原水位(流量)与实测水位(流量)进行比较,以此分析水库联合蓄水对中下游水位(流量)的影响。将各站考虑水库联合蓄水的天然水位(流量)与仅考虑三峡蓄水天然水位(流量)进行比较,以此分析溪洛渡、向家坝蓄水对中下游水位(流量)的影响,具体的技术路线见图3。
![](//html.hanspub.org/file/8-2410979x19_hanspub.png)
Figure 3. Mind map of the analysis of techniques
图3. 分析技术路线图
4.2. 对宜昌站水位影响
将考虑水库联合蓄水的还原水位、仅考虑三峡蓄水的还原水位以及实测水位进行对比,如表4。
可以看出,2014年9月,宜昌站各旬平均水位分别下降1.09 m、0.43 m以及0.76 m,其中溪洛渡、向家坝水库蓄水使宜昌站各旬平均水位分别下降0.35 m、0.38 m以及0.35 m,三峡水库蓄水使宜昌站各旬平均水位分别下降0.74 m、0.05 m以及0.41 m,溪洛渡和向家坝联合蓄水对宜昌站各旬水位下降的贡献率分别为32.1%、88.4%以及46.1%。2015年9月,宜昌站各旬平均水位分别下降0.84 m、1.60 m以及1.71 m,其中溪洛渡、向家坝水库蓄水使宜昌站各旬平均水位分别下降0.58 m、0.44 m以及0.53 m,三峡水库蓄水使宜昌站各旬平均水位分别下降0.26 m、1.17 m以及1.18 m,溪洛渡和向家坝联合蓄水对于宜昌站各旬水位下降的贡献率分别为69.0%、27.5%以及31.0%。
2014年10月,宜昌站各旬平均水位分别下降1.33 m、0.51 m以及0.15 m,其中溪洛渡、向家坝水库维持在正常蓄水位附近运行,对宜昌站水位影响较小,仅2015年10月下旬使宜昌站水位上升0.40 m。宜昌站10月份水位主要受三峡水库蓄水的影响。
2014年11月,宜昌站上旬平均水位上升0.72 m,中旬和下旬水位下降0.05 m以及0.16 m,其中溪洛渡、向家坝水库运行使宜昌站11月各旬平均水位略有上升。2015年11月,宜昌站上旬和中旬水位上升0.87以及0.99m,下旬水位下降0.22 m,其中溪洛渡、向家坝水库运行使宜昌站11月上中旬平均水位上升0.50 m、0.31 m,下旬水位下降0.27 m,三峡水库运行使宜昌站11月各旬平均水位分别上升0.37 m、0.68 m以及0.05 m。11月份宜昌站水位受三库蓄水影响程度较小。
可以看出,溪洛渡、向家坝两座水库运行对宜昌站的水位的影响主要集中在9月份,两座水库运行对宜昌站10月份和11月份的水位影响较小。三峡水库运行对宜昌站影响主要集中在9月下旬和10月份,三峡水库对宜昌站11月份水位影响较小。
![](Images/Table_Tmp.jpg)
Table 4. Impact of the Xiluodu and Xiangjiaba reservoirs impoundment on water level at the Yichang hydrological station (unit: m)
表4. 溪洛渡和向家坝水库蓄水对宜昌站水位影响分析表(单位:m)
4.3. 对下游控制站水位的影响
根据分析,溪洛渡、向家坝两座水库运行对长江中下游干流的水位和流量影响主要集中在9月份,两座水库运行对干流10月份和11月份的水位变化影响较小。以2014~2015年为例,进一步分析9月份各旬溪洛渡、向家坝水库,以及三峡水库对长江中下游干流水位的影响影响程度,见图4。
由图4可以看出,2014年9月上旬,溪洛渡、向家坝水库对于下游各个站点水位影响幅度要小于三峡水库,两座水库蓄水使中下游干流各站水位下降的幅度在0.00~0.35 m,三峡工程蓄水使干流各站水位下降幅度在0.25~0.74 m;2014年9月中旬,对宜昌、枝城、沙市以及螺山站,溪洛渡和向家坝水库对于这些站点水位影响幅度要大于三峡水库,对于汉口、九江以及大通站点,溪洛渡和向家坝水库对于这些站点水位影响幅度要小于三峡水库;溪洛渡、向家坝两座水库蓄水使中下游干流各站水位下降的幅度在0.01~0.38 m,三峡工程蓄水使干流各站水位下降幅度在0.05~0.26 m。2014年9月下旬,溪洛渡和向家坝水库对于各个站点水位影响幅度要小于三峡水库;长江干流各站水位变化主要受三峡影响,溪洛渡、向家坝两座水库蓄水使中下游干流各站水位下降的幅度在0.04~0.35 m,三峡工程蓄水使干流各站水位下降幅度在0.23~0.41 m。
2015年9月上旬,溪洛渡、向家坝水库对于下游各个站点水位影响幅度要大于三峡水库,两座水库蓄水使中下游干流各站水位下降的幅度在0.21~0.56 m,三峡工程蓄水使干流各站水位下降幅度在0.07~0.26 m;2015年9月中下旬,溪洛渡和向家坝水库对于这些站点水位影响幅度要小于三峡水库;2015年9月中旬,溪洛渡、向家坝两座水库蓄水使中下游干流各站水位下降的幅度在0.27~0.44 m,三峡工程蓄水使干流各站水位下降幅度在0.27~1.17 m。2015年9月下旬,溪洛渡、向家坝两座水库蓄水使中下游干流各站水位下降的幅度在0.32~0.53 m,三峡工程蓄水使干流各站水位下降幅度在0.60~1.18 m。
5. 结论
1) 基于mike模型构建了长江中游宜昌至大通的一维水动力演算模型,分析了溪洛渡、向家坝以及三峡水库联合蓄水对于下游河道水位的影响,模型模拟精度较高,得到了实测资料良好验证。
2) 根据分析,三个水库联合蓄水一定程度上降低了蓄水期下游各站点水位。溪洛渡、向家坝两座水库运行对长江中下游干流的水位的影响主要集中在9月份,两座水库运行对干流10月份和11月份的水位影响较小。以宜昌站为例,两座水库运行对宜昌站9月份旬水位的贡献率在27.5%~88.4%,三峡水库运行对长江中下游干流水位影响主要集中在9月下旬和10月份,三峡水库对干流11月份水位影响较小。
3) 以2014~2015年为例,分析了水库了联合蓄水对于中下游各个站点水位的影响。2014年9月溪洛渡、向家坝水库对中下游干流水位影响最大的是中旬,溪洛渡、向家坝两座水库蓄水使中下游干流各站水位下降的幅度在0.01~0.38 m,三峡工程蓄水使干流各站水位下降幅度在0.04~0.26 m。2015年9月溪洛渡、向家坝水库对中下游干流水位影响最大的是上旬,溪洛渡、向家坝两座水库蓄水使中下游干流各站水位下降的幅度在0.21~0.58 m,三峡工程蓄水使干流各站水位下降幅度在0.07~0.26 m。
基金项目
中国长江三峡集团有限公司资助(0704169);国家重点研发计划项目(2017YFC0405302)。
参考文献