1. 引言
乌江流域位于长江南岸,是长江上游南岸最大的入库支流,是三峡水库重要的水沙来源。乌江于重庆涪陵入汇长江,属于三峡水库常年回水区末端,乌江水沙直接进入三峡水库常年回水区内。因此,乌江水沙条件变化对三峡库区水沙及冲淤特性有着直接影响 [1] [2] [3] [4]。关于乌江流域水沙特性变化,陈松生 [5]、吴晓玲 [6] 等基于实测资料,对乌江流域水沙变化特征进行研究并开展了归因分析;邵俊 [7]、郭文献 [8] 等从多个方法入手,分析了乌江流域径流变化特征及归因识别;石国珏、许全喜等 [9] - [15] 从长江上游总体出发,分析了长江上游来沙变化及对三峡入库水沙的影响。以上研究成果对揭示乌江流域水沙变化及其对三峡入库泥沙的影响起到了一定的支撑作用。然而,已有研究成果多基于对实测资料成果的表观分析,对资料序列规律的深度挖掘还不够。本文采用多种统计理论,对三峡库区入库水沙控制站之一的乌江武隆站多年水沙资料开展深层规律分析。本文所采用的实测水沙资料均来自于长江水利委员会水文局。
2. 武隆站水沙变化特征
2.1. 水沙总体变化情况
武隆水文站位于乌江口上游约70 km处,设立于1951年,隶属于长江水利委员会水文局。武隆站径流量资料系列起讫年份为1952~2021年,输沙量资料系列起讫年份为1955~2021年。同时,为深入分析乌江水沙补给对三峡入库水沙特性的影响,将时间段划分为1991~2002年、2003~2012年、2013~2021年三个时段。
实测资料表明,乌江武隆站多年平均径流量和输沙量分别为485.8亿m3、2070万t,年径流量、输沙量最大值分别出现在1954年、1977年,最小值分别出现在2006、2013年。
根据对三峡水库不同蓄水时期乌江来水来沙变化情况的分析,武隆站1991~2002年多年平均径流量和输沙量分别为532亿m3、2040万t,与1990年以前相比,径流量增加39亿m3,增幅为8%,而输沙量减小930万t,减幅达31%;三峡水库蓄水运用后,2003~2012年武隆站年径流量有所减小,与1991~2002年均值相比减幅约为21%,而输沙量再次减少1470万t,减幅高达72%;金沙江下游溪洛渡、向家坝水库相继蓄水运用后,三峡入库水沙特性进一步发生显著变化,2013~2021年武隆站多年平均径流量和输沙量分别为495亿m3、298万t,与2003~2012年相比,径流量增幅为17%,而输沙量呈进一步减小态势,减幅为48%。总体来看,乌江武隆站年径流量呈现波动变化,在不同时期有增有减,而输沙量呈现持续减小态势,特别是2003年以来尤为明显。综上分析,对于径流量变化而言,可能存在着年际间的周期性变化规律,对于输沙量变化而言,年际间的跃变特征更为显著。武隆站历年来水来沙变化过程见图1。
![](//html.hanspub.org/file/3-2411172x8_hanspub.png?20220901081840730)
Figure 1. Annual flow and sediment transport process of Wulong station
图1. 武隆站年径流量与年输沙量变化过程
从不同时期乌江水沙补给对三峡入库水沙条件的影响程度来看(表1),乌江来水和三峡入库总径流量基本保持稳定,占比波动变化;而输沙量占比表现为持续减小,1990年以前,乌江来沙占三峡入库总沙量的6.1%;1991~2002年,乌江入库沙量占比减小至5.8%;2003年至今,占比又进一步减小至3.0%。值得注意的是,在三峡入库沙量整体大幅减少的情况下,乌江来沙占比仍在不断减小,说明乌江流域输沙量减少的情况较三峡入库总沙量减少更为突出。
![](Images/Table_Tmp.jpg)
Table 1. The ratio of flow and sediment between Wujiang and TGR at different time periods
表1. 不同时段乌江水沙补给占三峡入库水沙总量比值
注:在三峡建库前后,均采用朱沱站 + 北碚站 + 武隆站水沙之和作为三峡入库水沙总量。
2.2. 径流量周期性变化特征
从实测资料分析来看,武隆站多年径流过程呈现波动变化特征,具有一定的周期性。本文采用小波分析理论,对武隆站径流量的年际间周期性变化特征展开研究。
从武隆站径流量小波等值线(图2)可以看出,径流量从上至下存在着24~26a,15~16a,10~12a和6~7a (a为无量纲时间尺度)四类尺度的周期变化规律,且各尺度的周期变化波动性明显。从尺度较大的24~26a分析来看,武隆站径流量多年过程经历了4次丰枯交替,其中在1964~1969年、1980~1984年、1995~2001年、2009~2016年4个时期小波系数为正值,表示为历史多水期,其余各时间段为少水期,呈现显著的丰枯交替规律。而在其他较小尺度下,相关性相对较弱。同时从小波方差图(图3)也可以看出,1952年以来共出现了4个峰值,对应径流量变化的主周期分别为24a,16a,10a和6a,其中尤其以24a为第一主周期的变化最为显著。由此可见,武隆站径流量在一定程度上存在着丰枯交替的周期性变化规律。
2.3. 输沙量跃变特征
跃变是水文系统所具有的非线性特殊表现形式,指时间序列在某时刻发生急剧变化的一种形式,表现为跳跃点前后平均值发生陡升或是陡降。跳跃分析的本质是有序聚类分析,首先找出跳跃点,然后进行检验 [14]。一般把由子序列均值作为检测跃变的基本量而得到的跃变叫做第I类跃变,由子序列方差作为检测跃变的基本量而得到的跃变叫做第II类跃变。
受多种因素影响,长江上游各支流上泥沙跃变情况较为普遍且突出。根据武隆站1955~2021年共67年输沙
![](//html.hanspub.org/file/3-2411172x9_hanspub.png?20220901081840730)
Figure 2. Wavelet analysis of annual runoff at Wulong station
图2. 武隆站年径流量小波分析
![](//html.hanspub.org/file/3-2411172x10_hanspub.png?20220901081840730)
Figure 3. The main cycle of runoff changes at Wulong station
图3. 武隆站径流量变化主周期
![](//html.hanspub.org/file/3-2411172x11_hanspub.png?20220901081840730)
Figure 4. Jump phenomenon of sediment transport at Wulong station
图4. 武隆站年输沙量跃变图
量序列计算结果来看(图4),武隆站多年输沙过程经历了一次第I类跃变,两次第II类跃变。其中1983年为一级跳跃点,输沙量明显骤降,输沙量减小幅度达到了64.5% (见表2)。1966年、2003年为次级跳跃点,其中1966年为跃升点,输沙量明显增加,跳跃幅度为57.5%;2003年为跃降点,输沙量减小幅度为79.3%。
![](Images/Table_Tmp.jpg)
Table 2. The sediment transport in different periods at Wulong station
表2. 武隆站各分期年平均输沙量
3. 主要影响因素及后续影响分析
已有研究成果表明,影响长江上游地区流域产输沙的主要因素包括流域下垫面条件、降水、人类活动等3个主要方面 [10]。具体到乌江流域,主要影响因素可归纳为大型水利枢纽建设和土地利用方式改变两个方面。
3.1. 大型水利枢纽建设
乌江流域水能资源丰富,水利枢纽建设密度较大。乌江流域来沙减少的最主要原因就是水库拦沙,尤其是与大型水利枢纽建设有密切关联 [5] [6]。从乌江流域水库总库容变化来看(图5),1960年,支流猫跳河上的红枫湖水库建成,水库库容为7.53亿m3,乌江流域水库总库容迎来第一次显著增大;1982年,乌江干流乌江渡水库建成,这是乌江流域首座大Ⅰ型水库,水库库容23亿m3;2002年,支流六冲河上洪家渡水库建成,水库库容49.47亿m3;2009年,乌江干流构皮滩水库建成,水库库容64.54亿m3,乌江流域水库总库容增长迅速。
大型水利枢纽修建,极大程度上导致乌江河道泥沙输移受阻,进入三峡库区内的泥沙不断减少。实测资料表明,在这些大型水库蓄水运用后,下游输沙量锐减,拦沙作用显著 [16] [17]。如洪家渡水电站于2004年4月正式蓄水,蓄水前下游多年平均输沙量为611万t,而2004年实测输沙量仅为11万t,基本拦截了支流六冲河的全部沙量,对乌江流域输沙影响极为显著。
![](//html.hanspub.org/file/3-2411172x12_hanspub.png?20220901081840730)
Figure 5. Accumulative diagram of reservoir storage change in the Wujiang River basin since 1956
图5. 1956年以来乌江流域水库库容变化累计图
大型水利枢纽修建对水沙输移的影响可通过径流量–输沙量双累积曲线进行反映。从武隆站年径流量–输沙量双累积曲线来看(图6),1954年以来,武隆站径流–输沙关系主要经历了3次明显变化,分别出现在1966年、1983年、2003年,其中尤以1983年和2003年最为突出,双累积曲线斜率明显减小,反映出同径流量下输沙量大幅减少。这与2.3节得到的武隆站输沙量跃变节点基本一致,也与上述大型水库修建时间节点基本吻合,反映出大型水利枢纽修建对乌江水沙输移的直接影响。但2003年之后,尽管乌江流域水库建设仍在继续,但双累积曲线再未出现过显著转折,可见水库修建对河道输沙的影响已逐渐达到上限,输沙量进一步减少的空间十分有限。
![](//html.hanspub.org/file/3-2411172x13_hanspub.png?20220901081840730)
Figure 6. Double mass curve between annual runoff and sediment transport at Wulong station
图6. 武隆站年径流量–输沙量双累积曲线
3.2. 土地利用方式改变
土地利用方式改变对流域产输沙具有一定的影响,尤其是退耕还林(草)及水土保持工程,对治理流域水土流失具有重要意义。但土地利用方式改变对流域产输沙的影响效果相对滞后,与水沙变化关键节点的呼应较差,其中间往往存在一定的时间差。
乌江上游毕节地区是乌江流域水土流失较为严重的地区,自1989年开展“长江上游水土保持重点防治区综合治理工程”以来,水土流失情况逐渐好转,初步估计,1989~2004年毕节地区水土保持措施年均减蚀量在1000万t左右,极大程度从源头上减少了乌江入河泥沙总量,这也是2003年以来武隆站泥沙输移量锐减的主要原因之一。
从乌江流域土地利用方式变化来看,长期以来,乌江流域各土地类型变化特征具体表现为“四增两减”,即林地、建设用地、水域与未利用地的面积增加;耕地、草地的面积减少。根据2017年调查结果 [18],目前贵州乌江流域林地面积为35,741.05 km2,占流域总面积的比重为53.43%;其次为耕地,面积为21,056.33 km2,占流域总面积的比重为31.48%;草地面积为8344.43 km2,占流域总面积的比重为12.47%。林地、耕地与草地构成了贵州乌江流域土地利用的主体,三者面积之和占流域总面积的97.38%,其中林地和草地占比为65.9%。总体来看,乌江流域土地利用方式整体上是朝着流域产沙减少的方向发展的,这是未来乌江流域产沙继续维持较低水平的重要前提。
4. 水沙变化趋势分析
4.1. 径流量趋势研究
通过前文分析可知,武隆站径流量年际间呈现着周期性变化特征,存在着丰枯交替的变化规律。基于小波方差检验的结果,得到乌江武隆站径流演变的第一和第二主周期小波系数(图7)。在不同的时间尺度下,流域径流存在着平均周期及丰水期~枯水期变化特征。如在16a特征时间尺度上(图7(a)),流域的平均变化周期为10年左右,大约6个周期的丰水期~枯水期变化;在24a特征时间尺度上(图7(b)),乌江武隆站年径流量变化的平均周期为17年左右,大约经历了3个丰水期~枯水期的循环阶段。根据较大尺度(24a)的小波系数变化趋势,初步可以看出乌江武隆站年径流过程自2016年以来进入偏枯期,未来2~3年内武隆站来流将持续偏枯态势,后将由偏枯期转为偏丰期。
(a) 16a特征尺度
(b) 24a特征尺度
Figure 7. Wavelet real part process line of annual runoff change at Wulong station
图7. 武隆站年径流量变化小波实部过程线
4.2. 输沙量趋势研究
Mann-Kendall法是一种非参数统计检验方法,在水文学中常常被用来分析水文时间序列的变化趋势和突变年份。本节即基于Mann-Kendall秩相关检验法的统计理论,对乌江流域泥沙趋势性变化展开分析。
利用趋势分析,可判断时间序列中是否具有上升或是下降的趋势。其基本原理是:对于时间序列
,先确定其对偶值,记为P (即当
,且
出现的次数);若序列按顺序前进的值全部大于前一个值,即具有上升趋势,则
;若序列倒过来,即为下降趋势,则P = 0。因此,对于无趋势的序列
;当P < E(P)表示序列可能是有下降趋势的;P > E(P)时序列可能是上升趋势。此检验采用统计量:
(1)
当N增加时,U很快收敛于标准化正态分布。设置信水平α = 0.05,查正态分布表得临界值U0.05/2 = 1.96。
根据上述趋势分析方法对乌江武隆站的输沙量变化趋势进行分析,分析结果表明,在0.05显著水平下,乌江武隆站输沙量呈明显减小趋势。具体来看(图8),1968年以前,乌江武隆站年输沙量M-K值基本在0值以下,表明该时段年输沙量呈下降趋势,1968~1990年,年输沙量表现为上升趋势,而1991年以来,武隆站年输沙量再次呈下降趋势。2003年,UF和UB曲线相交,且交点在置信水平区间[−1.96, 1.96]内,说明2003年武隆站输沙量呈现突变性下降状态;此后M-K值突破了显著水平的临界值下限,年输沙量下降趋势显著。但考虑到当前乌江流域水利规划实施已较为充分,水利工程修建对泥沙输移的影响已接近饱和,尤其是受干流几座大型水库的拦截,武隆站的年输沙量已接近最低值,未来进一步减少的空间十分有限。因此,未来乌江流域进入三峡库区的泥沙量将维持当前的较低水平。
![](//html.hanspub.org/file/3-2411172x22_hanspub.png?20220901081840730)
Figure 8. M-K statistic series of annual sediment transport at Wulong station from 1950 to 2020
图8. 1955~2020年武隆站年输沙量M-K统计量序列
5. 结论
1) 实测资料表明,乌江武隆站多年径流过程(1952~2021年)呈现波动变化特征,具有一定的周期性,大尺度上存在着丰枯交替的规律;输沙量(1955~2021年)则表现为年际间跃变特征显著,其中1983年为一级跳跃点,输沙量明显骤降,减小幅度达到了64.5%。1966年、2003年为次级跳跃点,其中1966年为跃升点,跳跃幅度为57.5%;2003年为跃降点,减小幅度为79.3%。
2) 大型水利枢纽修建极大程度上阻碍了乌江河道泥沙输移,尤其是1982年乌江渡水库和2002年洪家渡水库的修建,对乌江流域输沙减少影响最为显著。但2003年之后,水库修建对河道输沙的影响已逐渐达到上限,加上乌江流域土地利用方式总体朝着流域产沙减少方向发展,未来乌江流域产输沙将继续维持当前较低水平。
3) 水沙变化趋势分析表明,武隆站年径流量呈现周期性变化规律,在24a特征时间尺度上,乌江武隆站年径流量变化的平均周期为17年左右,当前乌江正处于水量偏枯期,未来2~3年仍将维持偏枯状态,之后将转入丰水期;武隆站2003年以来年输沙量下降趋势明显,但今后输沙量进一步减少的空间有限,未来仍将维持当前较低水平。
基金项目
研究内容得到国家自然科学基金联合基金项目(U2040218)资助。