2020 年流域性大洪水下长江下游干流河道演变分析

2023-09-04 12:44罗红雨
水利水运工程学报 2023年4期
关键词:三峡水库大通河床

曹 双,程 龙,李 洁,罗红雨

(长江水利委员会 水文局长江下游水文水资源勘测局,江苏 南京 210011)

长江下游起于湖口,流经湖北、江西、安徽、江苏和上海五省市,干流全长约938 km,属于分汊河型,河道平面呈藕节状,宽窄相间。其中湖口至江阴段按照节点、束窄段或特定的支流河口可分成15 个河段[1],自上而下依次为上下三号、马垱、东流、官洲、安庆、太子矶、贵池、大通、铜陵、黑沙洲、芜裕、马鞍山、南京、镇扬、扬中河段(图1)。两岸入汇的主要支流有江西省鄱阳湖水系,安徽省华阳河、皖河、裕溪河、滁河、青弋江、水阳江、后河、秋浦河等,江苏省淮河入江水道及大运河、秦淮河等。河段中最大的江心洲是扬中的太平洲,河床最低处在江阴鹅鼻嘴。

图1 长江下游湖口至江阴河段河势Fig. 1 River map of Hukou to Jiangyin reach of the lower Yangtze River

长江下游干流多年来陆续实施了河道治理、航道整治、崩岸应急治理等工程,加之天然节点控制及人工强化的河道边界限制,河势基本稳定,但局部河段河势仍有较大变化[2-4]。2003 年6 月三峡水库开始蓄水以来,大量泥沙在库区淤积[5-6],清水下泄改变了长江中下游河道的水沙条件[7-8],河势变化特性及河道冲淤也产生了新的变化[9-13]。李义天等[14]曾通过与其他水库下游冲淤规律的对比研究,预测三峡水库下游的冲刷发展过程将比已有的计算成果缓慢;姚仕明等[15]研究认为,虽然长江中下游河段的河(航)道整治工程有效控制了总体河势,但三峡水库运行后河道仍向冲刷方向发展,并通过数学模型计算和物理模型试验预测2008—2022 年荆江河段将持续冲刷;许全喜等[16]对宜昌至大通干流河道冲淤规律分析发现,三峡水库蓄水运行后,长江中下游河道由三峡工程建设前“上冲下淤”和“滩淤槽冲”的特征转变为近乎全程的“滩槽普冲”;朱玲玲等[17]通过原型观测数据分析认为,三峡大坝下游河道进入了强冲刷状态。

2020 年长江流域性大洪水发生后,许全喜等[18]对2020 年长江中下游(宜昌至长江口)的河道总体冲淤变化分析发现,上游来沙和三峡水库出库泥沙明显增多、中下游干流洪水顶托严重、汛期水位快涨和汛后水位慢退是宜昌至湖口河段河床冲刷量偏小的主要原因;湖口至江阴段水面比降大,加之水流含沙量低,使冲刷量偏大。另有多名学者研究了2020 年大洪水对个别洲滩、分汊河段的影响[19-20],得到了一些局域性河道演变特征。本文在以往研究成果的基础上,对2020 年流域性大洪水作用后长江下游干流河道水沙变化、河道冲淤及河势演变特征等进行较为全面及深入的研究和探讨。

1 材料与方法

研究河段内设有九江、大通两个常年水文站,水文资料采用两站的长系列水位、流量和泥沙观测数据,以及沿程各主要汊道测验数据等。河道地形采用1998 年以来长江下游水文水资源勘测局测绘的河道地形图及断面实测资料,其中包含了2020 年11 月实测大断面资料、2021 年3 月1∶10 000 长程地形资料。

本次冲淤分析采用断面法计算冲淤量。根据河段内主要水文站的水位流量关系及河段地形情况确定水面线,分枯水、 平均 、平滩 和洪水4 个水位进行计算,对应大通站流量分别为15 000、30 000、45 000 和60 000 m3/s。固定断面平均间隔约1.0 km,弯曲段及变化较大部位的间隔缩小至500 m 左右。河道断面一般按1∶5 000 进行数据提取,即横向平均点距35~40 m,陡坎及局部地形变化剧烈区域横向点距加密至约20 m。采用截锥公式计算相邻两断面间的河道槽蓄量。

式中:Vi(zj)为在水位zj条件下,断面i和断面i+1 间的河道槽蓄量;Ai为断面i的过水断面面积;Li为断面间距。

河段槽蓄量V由各断面间槽蓄量累加得到,两测次之间的槽蓄量之差即为冲淤量 ∆V。冲淤强度 ∆G为单位河长年平均冲淤量,即∆G=∆V/(Ln),其中:L为河段长度,n为两测次间隔的年数。采用套图比较法分析河道地形和断面的变化规律。

2 三峡水库蓄水以来长江下游水沙特性变化

2.1 年际变化

以三峡水库蓄水时间(2003 年)为界,分蓄水前及蓄水后两个时段进行水沙特性空间差异性统计分析。对比三峡水库蓄水前后数据可知,九江水文站径流量减小4.6%,输沙量减小67.8%(图2(a));大通水文站径流量减小2.5%,输沙量减小69.5%(图2(b))。三峡水库蓄水后,下游径流量微减,而输沙量锐减,输沙量减幅显著大于径流量的减幅,水沙关系发生重大改变。

图2 九江站、大通站年径流量及输沙量变化Fig. 2 Changes in annual runoff and sediment runoff at Jiujiang and Datong hydrological stations

2.2 年内变化

大通水文站是长江下游流量控制站,其下游干流区间入江流量仅约占其流量的3%,大通水文站的流量、泥沙特征基本代表长江下游来水、来沙特征。根据1950 年以来大通站的水沙资料,以2003 年为界分别统计流量和输沙率的年内差异(图3)。

图3 大通站月均流量及输沙率变化Fig. 3 Variation of monthly average discharge and sediment transport rate at Datong hydrological station

由图3 可见,长江下游水沙年内分配不均,存在明显的汛期和枯水期。三峡水库蓄水以后,汛期洪峰有所坦化,枯水期流量略有增大,常年流量不低于10 000 m3/s;沙峰的坦化更为明显,三峡水库蓄水后的汛期输沙率远小于蓄水前和1950 年以来的平均值,长江下游汛期的不饱和输沙现象更为突出。三峡工程对长江下游输沙年内分配的影响远大于对径流的影响。

2.3 2020 年流域性大洪水特性

2020 年,长江发生了自1949 年以来仅次于1954、1998 年的流域性大洪水,长江干流发生5 次编号洪水。其中,长江上游发生特大洪水,寸滩站洪峰水位居实测记录第2 位,三峡水库出现建库以来最大入库流量;长江中下游干流监利至大通段洪峰水位列有实测记录以来的第2~5 位,马鞍山至镇江段潮位达历史最高[21-22]。大通站于7 月13 日达到最高水位16.24 m,仅低于1954 和1998 年;同时达到了最大洪峰流量83 800 m3/s,该值仅次于1954 和1999年,位列历史第3;最大1 日洪量和最大3 日洪量介于1954 和1998 年之间,但最大7 日和最大15 日洪量均低于1954 和1998 年。大通站2020 年汛期日均流量、水位过程见图4。

图4 2020 年长江洪水大通站水位、流量过程Fig. 4 Process of water level and flow at Datong hydrological station during the 2020 Yangtze River flood

长江下游河段水情还与鄱阳湖来水密切相关。2020 年,鄱阳湖虽发生流域性超历史大洪水,湖口站洪峰水位居实测记录第2 位,鄱阳“五河”最大合成流量43 200 m3/s,但湖口站年径流量(1 547 亿m3)却维持在多年均值附近。这主要是由于长江水位持续居高,顶托鄱阳湖形成持续高水位,造成湖口出流不畅,进一步延缓了长江洪峰过后的退水过程。

2020 年主汛期,长江流域具有入梅早、出梅晚、梅雨量大,暴雨强度大、极端性强,暴雨阶段性分布明显、各阶段强雨区重叠度高等特征。在长江中下游超长超强梅雨和四川盆地强降雨影响下,2020 年长江流域性大洪水主要有以下6 个特点:(1)上游来水早,洪水发生范围广;(2)干流区间洪水突出;(3)中下游干流水位涨势猛;(4)中下游洪峰水位高、高水位持续时间长,南京站出现最高水位10.39 m,超历史最高水位0.17 m,超警戒水位48 天,超保证水位6 天;(5)上游洪水峰高量大,大通站最大30 天洪量接近1998 年水平;(6)水库群联合调度影响显著[21,23]。

3 河道冲淤变化特性

3.1 湖口至大通河段

计算湖口至大通河段1998—2020 年间各时段的冲淤量、冲淤强度见表1。可见湖口至大通段河床累计为冲刷。三峡水库蓄水前,除枯水期外该河段河床总体表现为淤积;三峡水库蓄水后,该河段河床转淤为冲,其后呈单向冲刷态势。具体特征为:(1)2001—2006 年三峡水库蓄水前后,河床呈槽淤滩冲,且冲刷量、冲刷强度均随水位抬高而增大,说明该时段冲刷部位主要在河滩;(2)2006—2011 年上游水库群运用前,河床表现为滩槽普冲,但冲刷相对较缓;(3)2011—2016 年上游水库相继投入运行,冲刷逐渐增大,这期间各计算水位下河槽冲刷量较为接近,说明冲刷主要发生在枯水河槽;(4)2016—2020 年上游水库群实现初期蓄水后,冲刷有所减缓,但整体冲刷强度仍高于2011 年以前水平。

表1 湖口至大通段冲淤量和冲淤强度统计Tab. 1 Statistics of erosion/deposition amount and intensity from Hukou to Datong

2019—2020 年河床冲刷幅度显著增大,冲刷量占2016—2020 年5 年冲刷量的50%左右,冲刷强度是2016—2020 年平均值的2.0~2.3 倍。2001—2020 年间,该段河床平滩河槽冲刷总量达49 682 万m3,其中2020 年洪水作用后的冲刷量为6 501 万m3,占比13.2%,年均冲刷强度是2001 年以来均值的2.5 倍。

3.2 大通至江阴河段

大通至江阴河段1998—2020 年间各时段的冲淤量、冲淤强度见表2。1998—2020 年大通至江阴段河床呈单向冲刷趋势。三峡水库蓄水前的1998—2001 年,该河段河床冲刷相对较小,冲刷量随计算水位升高而减小,说明该河段河滩总体呈小幅淤积态势;三峡水库蓄水后,该河段河床保持单向冲刷态势。具体特征为:(1)2001—2006 年三峡水库蓄水前后,河床滩槽普冲,冲刷相对较缓,冲刷量及冲刷强度随水位抬高而增大;(2)2006—2011 年,冲刷强度显著增大,约为前一时段的2.5 倍;(3)2011—2016 年上游水库群相继投入运行,冲刷又出现一定程度的减小,但仍约为2001—2006 年的2 倍;(4)2016—2020 年上游水库群实现初期蓄水后,河床冲刷再次加剧。

表2 大通至江阴段冲淤量和冲淤强度统计Tab. 2 Statistics of erosion/deposition amount and intensity from Datong to Jiangyin

由表2 可见,与1998—2020 年冲淤相比,2019—2020 年河床冲刷强度显著增大,2019—2020 年河床冲刷量约占2016—2020 年冲刷量的60%,冲刷强度是2016—2020 年平均值的2.3~2.8 倍;1998—2020 年间,该河段河床平滩河槽冲刷总量达124 101 万m3,其中2020 年洪水作用下的冲刷量为21 478 万m3,占比17.3%,年均冲刷强度是1998 年以来均值的3.8 倍。

2020 年大洪水作用下,大通至江阴河段洪水河槽冲刷强度达54.31 万m3/(km·a),是上游湖口至大通河段的1.8 倍,其余各水位河槽冲刷强度也是上游河段的1.6~2.2 倍。湖口至大通河段与大通至江阴河段的来水来沙条件相近,两者冲刷强度不同主要是受河床边界条件的影响,长江下游干流床沙总体趋势为沿程变细,故而同等条件下越向下游越易冲刷。已有研究[24]表明,随着三峡水库的蓄水运用,坝下游河道发生沿程冲刷,并逐步向下游发展。上段河床经过多年剧烈冲刷后,冲刷强度已有所减缓,而下段则正处于剧烈冲刷时期,在大洪水的作用下冲刷强度增大更为显著。

3.3 南京河段

2019—2020 年,长江下游除南京段出现微淤[19](图5,淤积强度为12.6 万m3/(km·a)),其他河段以冲刷为主。长江湖口至南京段没有大的支流入汇,该河段流量差仅约2%。在经历湖口以下多个河段约430 km河长连续冲刷后,江水含沙量应有明显增大。但事实上,通过对比大通水文站及南京水文实验站的含沙量资料发现,南京站含沙量反而小于上游的大通站(图6),鉴于两测站流量差甚小,可推知有大量泥沙在南京河段落淤,由表3 统计的数据可见,从新济洲汊道段开始出现了较为显著的淤积。

表3 2019—2020 年南京河段冲淤量统计Tab. 3 Statistics of erosion/deposition amount in Nanjing reach during 2019-2020

图5 2019—2020 年长江下游河段平滩河槽冲淤强度Fig. 5 Variation of mean erosion/deposition intensity under bankfull water level along the lower Yangtze River during 2019-2020

图6 2019—2020 年大通站和南京站含沙量变化Fig. 6 Variation of sand content at Datong and Nanjing hydrological stations during 2019-2020

初步分析得出南京河段冲淤异常的成因主要包括两个方面:一是由于2020 年大洪水期间叠加了高潮位,在南京河段形成了长时间与流量不匹配的高水位,使水流流速下降明显;二是长江南京河段先后进行了9 次较大规模的河道整治工程,护岸总长约110 km[20]。整治工程的实施,使水流顶冲岸段的岸线得到有效保护,基本控制了河势,河床冲刷被限制在江心较窄范围内。南京河段冲淤异常的详细成因及其影响机理尚需进一步研究。

4 河势演变特征

20 世纪50 年代以来,长江下游河道在自然因素及人工干预下演变,且人工干预的影响近年来明显增强。具体表现为:总体河势基本稳定,两岸岸线相对稳定,崩岸频次及总长有所减小,河道冲淤基本限制在远离河岸的江心河槽部位;三峡水库蓄水运行以来,长江下游河道几乎呈单向冲刷态势,部分河段冲刷幅度较大;人工干预未改变河道演变的基本规律。2020 年长江流域性大洪水作用虽未改变总体基本稳定的态势,但仍造成了河床冲刷明显加剧、局部河段河势有所调整、个别河段河势剧烈变化等影响。

4.1 主流走势与分流格局

三峡水库蓄水以来,长江下游主流走势及各洲滩分流格局总体相对稳定,但局部河段随洲头低滩冲刷变化,分流区及进口段主流走势有所调整;有些支汊分流比呈微增趋势,个别河段出现主支汊易位趋势。

湖口至江阴段20 多个汊道大部分保持了相对稳定的态势,如南京河段梅子洲汊道,多年来主支汊分流比稳定在19∶1 左右;八卦洲汊道在经历了多轮河道整治工程后,左汊缓慢衰退的趋势得到进一步遏制,近5 年汛期分流比基本维持在15%以上,2020 年大洪水后,分流格局未发生大的变化(图7(a))。

图7 部分汊道分流比变化Fig. 7 Variation of diversion ratio of selected branching channels

2020 年大洪水冲刷后,太平洲、落成洲等原已表现出支汊发展势头的汊道,其支汊发展趋势进一步巩固。太平洲进口段主流进一步右偏,使得落成洲右汊及太平洲右汊更易获得进流,支汊冲刷在大水作用后明显加速。2021 年7 月,太平洲右汊实测分流比达13.7%(图7(b)),为多年来实测最大值,落成洲右汊分流比则达24.1%(图7(c)),超越了航道部门控制目标(20%)和整治前水平。

铜陵河段成德洲长期以左汊为主汊,1959—1998 年左汊分流比为56.5%~66.3%,1998 年以后有所下降,但仍维持主汊地位。2016 年9 月和2022 年2 月两次实测右汊分流比超过左汊,达到53.4%(图7(d))。可见在右汊冲刷发展的汊道演变趋势下,经过2016、2018、2020 年连续多次大水冲刷后,成德洲已呈现出主支汊易位状况。

4.2 岸线与滩槽格局

长江下游河道经过多年治理,岸线、滩槽格局相对稳定,大部分河段的变化仅限于局部小幅度调整。但个别河段有崩岸发生,洲滩相对薄弱,江心洲依附的边滩有多处明显冲刷。

扬中河段是长江下游崩岸较为频发的河段之一。2017 年11 月在太平洲左缘指南村附近江岸发生较大尺度灾害性窝崩。经应急抢险和系统治理,2020 年大洪水后,虽在崩窝外侧河床形成−60 m 冲刷坑,但未发生新的窝崩险情。在系统加固治理工程保护下,2020 年大洪水后未发生类似灾害性崩岸,但仍有局部高滩陡坎条崩。如落成洲右汊丰乐桥、雷公岛一带沿岸流增强,丰乐桥出现条崩,右汊分流比骤然回调。

滩槽变化多为依附在江心洲的边滩局部调整。马垱河段棉船洲洲头左缘河滩崩坍,左汊进口段河槽右移;安庆河段江心洲左缘下段2.6 km 范围河滩崩坍(图8(a)),中汊出口段深槽右摆约100 m,略向顺直发展;太子矶河段铜板洲洲头右缘边滩发生趋势性冲退,2016—2021 年累计冲退480 m(图8(b));贵池河段长沙洲右缘中下段约3.8 km 范围内边滩大幅冲退超300 m(图8(c));铜陵河段汀家洲左缘铜陵沙发生撇岸切滩,凸岸冲刷(图8(d)),2016—2021 年最大累计冲退约340 m;芜裕河段新沙洲由于未守护,累计趋势性冲退1.4 km,洲尾下延,2016—2021 年冲退近400 m,面积缩小约1/4。

图8 部分洲滩岸线变化Fig. 8 Shoreline changes on selected beaches

上述滩槽局部调整多发生在未经守护或守护相对薄弱的江心洲,且对上游来水的水沙特性较为敏感。河滩崩坍或冲退虽为2016—2021 年间形成,但较三峡水库蓄水以来平均的冲刷变化有明显加剧,综合本文前述冲刷强度的变化数据可推知,2020 年大洪水的造床作用对其影响最为显著。

4.3 河道断面

近年来,长江下游河道总体上以冲刷为主,纵断面进一步下切。以马垱河段棉船洲右汊为例,2001—2021 年深泓纵剖面以冲刷下切为主(图9),20 年间纵剖面最大冲深约14.5 m,平均冲深约2.0 m。

图9 马垱河段棉船洲右汊深泓纵剖面变化Fig. 9 Changes in the longitudinal profile of the right branch of Mianchuan Island in Madang reach

河道束窄段横断面形态有进一步窄深化的趋势,例如官洲河段进口段(图10(a)),受吉阳矶节点控制河道束窄,主流自右向左过渡,2001 年以来左缘稳定,右缘淤积束窄,深槽略有冲深;南京河段龙潭河口至三江口段亦为束窄段(图10(b)),主流贴右岸下行,近年来左缘淤积,深槽略有冲深,河相系数由2.3减小到1.6。鉴于河道束窄段河势通常较展宽段和分汊段更为稳定,其变化趋势可在一定程度上反映河段总体的演变趋势,长江下游河段束窄段多有窄深化的发展趋势,也印证了河段总体上呈冲刷态势。

图10 长江下游部分束窄段河道横断面变化Fig. 10 Changes in river cross-sections in some narrow sections of the lower Yangtze River

5 结 语

2020 年流域性大洪水后长江下游总体河势基本稳定,大多数河段冲刷加剧,且呈现滩槽普冲态势,各水位河槽冲刷强度为2016—2020 年平均值的2.0~2.8 倍,大通至江阴河段冲刷强度是上游湖口至大通河段的1.8 倍,单河段以芜裕河段冲刷强度最大。南京河段不冲反淤,一方面是由于高潮位顶托降低了流速;另一方面与相对较为完善的河道守护限制了冲淤变化范围有关。南京河段的淤积成因与机理尚需深入研究。各汊道主流线走势整体稳定,分汊及汇流区仍有一定摆幅;整体分流格局相对稳定,但铜陵河段成德洲汊道有主支汊易位趋势,多次测得右汊分流比超过左汊;扬中河段支汊分流比持续微增,2020 年大洪水冲刷使其支汊发展趋势进一步巩固,虽不危及主支汊格局,但趋势很难扭转。

长江下游两岸岸线基本稳定,局部虽有条崩发生,但基本无指南村类似的大型、灾害性窝崩发生,由于2020 年高水位持续时间长,不少江心洲高滩、陡坎长期浸泡,未护段及薄弱段条崩频发。江心洲洲滩冲刷明显,洲滩面积萎缩,以江西、安徽两省尤为明显,主要原因是洲头守护工程相对于江苏薄弱,铜陵河段的铜陵沙一带出现撇岸切滩现象,滩槽冲淤仍较为频繁。长江下游河道演变规律虽未改变,但人工干预日益增强,双重影响下长江下游水沙特性及河床冲淤、河势变化将更加复杂,后续应进一步加强原型观测及分析研究工作。

猜你喜欢
三峡水库大通河床
距离完美只差一点 上汽大通MAXUS G90
三峡水库迎战长江2020年第2号洪水
寻找你心中的完美之车
——上汽大通D90……虞展
上汽大通:C2B个性化定制未来
上海大通的黄金十年
ArcGIS在河床冲淤量分析中的应用
三峡水库澎溪河水-气界面CO2、CH4扩散通量昼夜动态初探
缱绻的回眸诗歌
河床下伏采空区岩移规律及治理方案研究
三峡水库综合管理创新研究