河道清淤工程对防洪的累积影响研究.pdf

2023.10孙光宝1，刘红奎2，李宗骏1（1.长江水利委员会水文局荆江水文水资源勘测局，湖北 荆州 434099；2.湖北省襄阳市水文水资源勘测局，湖北 襄阳 441003）揖摘要铱 为研究多个取水口清淤对防洪的定量影响，文章以长江柳林洲新厂河段的（颜家台泵站、耀兴电灌站、普济中心水厂）3 个取水口为研究对象，构建二维水动力学模型，分别采用 2019 年 7 月 21 日和 7 月 24 日的流速和水位资料进行模型率定和验证。结果表明：模型模拟结果良好，断面流速计算值和实测值偏差在依0.1 m/s 以内。清淤疏浚工程实施后，对局部行洪有利，对附近水域的水位、流速以及护岸工程没有明显的影响。对比颜家台泵站下移前后的水文特征，发现多个清淤工程的叠加影响大于单个工程，但小于其线性叠加值。揖关键词铱 防洪影响；长江取水口；河道清淤；二维水动力学模型揖中图分类号铱 TV87揖文献标志码铱 A揖文章编号铱 1009-6159渊2023冤-10-0064-05本文研究河段上起沙市盐卡港（荆 46 断面）下迄新厂（荆 82 断面），长度约 66 km，整体位于上荆江下段沙市至藕池口，流经荆州市沙市区，公安县和江陵县，为微弯分汊型河道。拟实施清淤项目位于长江柳林洲-新厂河段左岸，河段为沙质河床，具有二元结构特征，主流线摆动幅度较大，河岸主要由厚层黏土组成，且受到两岸护岸工程保护，多年来工程河段岸线较为稳定1。自三峡工程建成以后，荆江水沙条件发生了较大的变化2，加之长江航道部门实施筑坝拦洪，主泓南移工程后，北岸泥沙淤塞严重等因素影响，导致研究区内颜家台泵站、耀兴电灌、普济中心水厂取水口无法正常取水。因辖区内种养殖结构发生变化及 荆州市四湖总干渠污染防治三年行动计划（20182020 年）的要求，为确保四湖总干渠流域内河水水质达到芋类标准，同时兼顾农业种植结构调整后农业常态化用水需求，需对灌区内水系常年进行生态补水。因此，需对颜家台泵站、耀兴电灌、普济中心水厂取水口开展清淤。1水动力模型构建与验证在水动力模型构建中，首先要考虑平面尺度与垂直尺度的相对关系，本文所选取的长江柳林州新厂河段为典型的大江大河天然河道，平面尺度远大于垂直尺度，故可构建水动力二维模型进行相关计算求解3。为准确评估取水口清淤对河道演变的实际影响，基于研究区域河床演变趋势初步判断结果，以实例进行验证，通过对比计算结果和实测成果，验证模型的适应性，用于计算和分析工程实施后对河道水位、流速的实际影响。本文基于 Mike21 构建二维水动力模型，分析工程河段的流速变化情况。本次选用 Mike21 模型中的 Hydrodynamic（HD）水动力模块，利用Saint-Venant 方程组描述二维非恒定流运动规律，由质量守恒的连续性方程和能量守恒的动量方程组成，可模拟湖泊，河道的水流运动规律，利用 ADI 法求解。河道地形大部采用 2018 年 10 月 1颐10 000地形图，研究区域局部地形采用 2019 年 7 月1颐2 000 地形图。其中耀兴电灌站段 1颐2 000 地形图 2.3 km2，颜家台闸段 1颐2 000 地形图 3.2 km2，普济中心水厂取水口段 1颐2 000 地形图 3.8 km2。选取郝穴站基本水尺断面和新厂水位站基本水尺断面为水位观测断面；选取荆 53、荆 58 左汊、荆58 右汊为流量测验断面。分别选取 2019 年 7 月24 日（26 600 m3/s）和 7 月 21 日的（21 000 m3/s）的收稿日期：2023-05-31作者简介：孙光宝（1986），男，工程师山东水利64窑窑2023.10表 3 工程影响计算水流条件表工况流量/（m3 s-1）下边界水位/m备注1998 年洪水53 70039.24多年平均12 10029.5820032018 年设计流量9 40027.30颜家台闸设计水位 27.57 m枯水期6 00026.2320032018 年表 2 汊道分流比实测值与计算值对比统计表断面名施测时间郝穴站水位/m流量/（m3 s-1）实测分流比计算分流比荆 537 月 24 日12：3235.9826 600荆 58左汊7 月 24 日17：0835.97730027.428.1荆 58右汊7 月 24 日18：4635.9718 70070.371.9图 1 26 000 m3/s 流量下流速率定结果渊荆 53冤资料进行率定和验证。经试算率定并参照经验值得本河段主槽的糙率取值范围为 0.0170.020，滩地的糙率取值范围为 0.0200.025。表 1 为各断面水位计算值与 2019 年 7 月24 日实测值的差值统计表。结果表明两者基本吻合，最大误差在 0.01 m 以内。结果表明：所建模型可很好地描述各断面情况。以图 1 为例，计算河段率定流场图及断面流速分布率定及验证结果。可见二维数学模型计算所得流场变化平顺，滩、槽水流运动区分明显；断面流速计算值和实测值吻合较好，流速偏差在依0.1 m/s 以内。2019 年 7 月 24 日 实 测 长 江 干 流 流 量26 600 m3/s，突起洲左汊分流量 7 300 m3/s，突起洲右汊分流量 18 700 m3/s，差值 600 m3/s，占干流流量的 2.3%，在测验误差允许范围内。实测和计算的汊道分流比结果见表 2。由表 圆 可知：干流与汊道流量测验时段内郝穴水位站水位变幅 0.01 m，水位过程与上游沙市站、下游新厂站相应，水位变幅基本一致，符合恒定流模型要求。综合对比和分析结果，经率定后平面二维数学模型可以较准确地模拟本河段的水流运动特性，计算结果和实测成果误差较小，吻合较好，可用于计算和分析取水口清淤工程实施后对研究区域水位和流速的影响。2清淤方案对实际工程的影响为了反映清淤疏浚对防洪、通航等的影响，基于 3 个水文站（沙市水文站、郝穴水位站、新厂水文站）的资料选取了 1998 年洪水、多年平均、枯水期、设计流量共四级流量进行工程影响分析计算（见表 3）。考虑到 1998 年洪水沙市站最高水位达 45.22 m（冻结），故选取 1998 年洪水期实测最高水位和最大流量作防洪设计水位和流量。计算河段下边界在新厂站基本水尺断面下游约3.4 km，下边界控制水位用新厂站实测水位、河段比降推算。表 3 中下边界设计水位 27.30 m 按郝穴、新厂两站枯季比降推出。颜家台闸站原设计水位为27.57 m，原设计流量 5 500 m3/s，在设计初期是正确的。但河道经几十年的演变，特别是三峡水库蓄水后河道冲刷下切，枯季水库调蓄运用下泄水量加大，枯水期同水位条件下河道流量加大较多，再用 5 500 m3/s 作为设计流量不符合实际。因此，根据 2018 年枯季郝穴站水位与沙市站流量关系推求出颜家台闸设计流量为 9 400 m3/s，枯季流量为 20032018 年沙市站年最小流量平均值取整，下边界水位为新厂站同期年最低水位平均值。当建筑物尺寸相对网格尺寸较小时，假定河底高程增加值所阻挡的流量与工程阻挡的流量表 1 各断面水位计算值与实测值对比统计表 m断面名实测计算差值荆 5337.3537.346-0.004荆 58 左汊37.0537.047-0.003荆 58 右汊37.137.098-0.002郝穴站基本水尺35.9835.978-0.002孙光宝，刘红奎，李宗骏：河道清淤工程对防洪的累积影响研究65窑窑2023.10表 5 不同计算方案下取水区域水位流速最大影响值清淤区域洪水组次流量/（m3 s-1）流速最大增加值/（m s-1）流速最大减小值/（m s-1）耀兴电灌98 年洪水53 7000.070.12多年平均12 1000.110.14设计流量9 4000.120.16枯水期6 000/颜家台闸98 年洪水53 7000.130.16多年平均12 1000.160.21设计流量9 4000.190.25枯水期6 0000.210.24普济中心水厂98 年洪水53 7000.120.16多年平均12 1000.170.24设计流量9 4000.210.25枯水期6 0000.230.3表 4 不同工况下清淤区域水位最大影响值清淤区域洪水组次流量/（m3 s-1）水位最大壅高/cm水位最大降低/cm耀兴电灌98 年洪水53 7001.140.76多年平均12 1002.41.77设计流量9 4001.893.26枯水期6 000/颜家台闸98 年洪水53 7005.273.86多年平均12 1004.623.05设计流量9 4001.61.15枯水期6 0000.520.88普济中心水厂98 年洪水53 7006.854.91多年平均12 1002.562.51设计流量9 4001.51.77枯水期6 0001.331.5相同，通过增加工程所在网格节点的河底高程来反映工程的阻水影响；当建筑物尺寸大于或与网格尺寸相当时，可直接根据建筑物高度来修改相应网格节点的河底高程4-5。假定垂线流速沿水深呈指数分布，则垂线上某点流速 u 可表示为6-7：u=u0伊（y/h）m（1）式中：y 为该点离床面距离，m；u0为 y=h 处水流表面流速，m/s；h 为水深，m；指数 m 可取 1/6。令 a为工程阻水增加的河底高程，b1、b2分别为工程阻水宽度和网格宽度，则工程阻挡的流量Q1和河底高程增加阻挡的流量 Q2分别为：a=h伊（b1/b2）1/（2+m）（2）Q1=b1u0h0（yh）mdy=b1u0hm+1Q2=b2u0h0（yh）mdy=b2u0hm+1扇墒设设设设设设设缮设设设设设设设（3）由 Q1=Q2可解出河底高程增加值。2.3.1原始清淤方案下模拟计算利用构建的水动力平面二维模型计算了1998 年洪水、多年平均、枯水期、设计流量共四级流量工况下，3 个水利工程取水口清淤引起的水位和流速变化，清淤引起突起洲分流比变化，以及对岸坡稳定影响的变化。1）水位变化分析。在设计流量和多年平均流量工况下，3 个取水口清淤疏浚前后的水位变化见表 4。由表 4 可知，中洪水流量下清淤疏浚后，清淤区上游水位略有降低，清淤区下游局部水位略升高；枯水设计流量下疏浚区水位变化规律基本一致，但因工程前不过水，清淤疏浚后过水，受计算方法影响，清淤区显示为水位降低。由模型计算结果可知，耀兴电灌站清淤疏浚后该区域水位变化较小，水位变化在 0.2 cm，影响范围主要位于工程区域上游约 1 km 范围内；其中在枯水期文村夹左汊基本不过流，此时耀兴电灌处无流量，因此无水位、流速变化。颜家台闸清淤疏浚后，在设计流量条件下清淤疏浚区尾端水位壅高，最大壅高值约为 5.27 cm；清淤疏浚区上游局部水位下降，最大水位下降约3.86 cm。普济中心水厂清淤疏浚后，水位最大壅高值约为 6.85 cm，最大水位下降约 4.91 cm。两者在设计流量下水位变化影响区域较小，水位变化在 0.2 cm 的影响范围局限在清淤疏浚区上游600 m 至下游 400 m 内。两者各自多年平均流量、设计流量和枯水期流量下的水位变化规律与设计流量条件下的变化规律基本一致。2）流速变化分析。工程前后采样点流速变化等值线及流速变化值见表 5。根据模型计算结果可知，耀兴电灌由于上下游筑有潜丁坝等航道整治建筑物，施工后该区域水位流速变化较小，流速变化在 0.01 m/s 影响范围主要位于工程区域上游约 1 km 范围内。颜家台闸清淤疏浚后，在设计流量条件下，清淤疏浚区上游和下游局部流速增大，流速最大增加值约为 0.13 m/s，流速最大减小值 0.16 m/s。普济中心水厂流速最大增加值约为 0.12 m/s，流速最大减小值 0.16 m/s。这两者多年平均流量下的流速变孙光宝，刘红奎，李宗骏：河道清淤工程对防洪的累积影响研究66窑窑2023.10表 8 不同计算方案下取水区域水位流速最大影响值清淤区域流量/（m3 s-1）流速最大增加值（m s-1）流速最大减小值（m s-1）颜家台泵站下移 3 km12 1000.160.22表 7 清淤区域水位最大影响值清淤区域洪水组次流量/（m3 s-1）水位最大壅高/cm水位最大降低/cm颜家台泵站（假定下移 3 km 后）多年平均12 1003.872.54表 6 不同计算方案下突起洲分流比变化化规律与其各自设计洪流量条件下的变化规律基本一致，但流速变化值略有所变化。枯水期流量下都因施工前河床不过水、工程后河床过水，导致施工后清淤疏浚区域及附近区域流速变化较大。3）突起洲分流比变化。各级流量下工程前后突起洲左右汊分流流量及分流比变化情况如表 6所示。各级流量下，工程前后突起洲左右汊分流情况无明显变化。入口流量为 6 000 m3/s 时，突起洲左汊保持不过流状态，主要从右汊过流。随着入口流量增加，左汊分流比逐渐增加，右汊的分流比逐渐减小，工程过后，突起洲左右汊流量分流比变化较小，其中 12 100 m3/s 时，工程对分流比影响最大，但此时分流量的变化只有个位数、分流比变化幅度仍较小。4）清淤对岸坡稳定影响分析。三大取水口距离上下游护岸工程至少有 100 m 以上的距离，因此三大取水口清淤作业对护岸工程、取水工程和航道工程安全无不利影响，且清淤区底部处在近期床沙活动层上部，清淤不会对岸坡稳定产生不利影响。2.3.2河流模拟计算由前述的模拟结果可知，相距较远的 3 个取水口清淤所引起的研究河段水文特征值的变化基本不会相互影响。为进一步研究距离较近的取水口清淤对所在河段水文特征值所造成的累积影响效应，假定颜家台泵站沿长江河道下移 3 km，下移后颜家台泵站与普济水厂取水口两个清淤区之间相距不足 1 km。在此假设背景下，选取多年平均流量工况重新对颜家台泵站清淤后的河流进行模拟计算，为进一步探讨两个距离较近的取水口清淤所造成的叠加效应提供对比方案。1）水位变化分析。多年平均流量下清淤的水位变化见表 7。由表 7 可知：多年平均流量下，清淤疏浚工程实施后，清淤区上游水位略有降低，清淤区下游局部水位略升高。由模型计算结果可知，颜家台闸下移 3 km实施清淤疏浚后，在多年平均流量条件下，本次拟清淤疏浚区尾端水位壅高，最大壅高值约为3.87 cm；清淤疏浚区上游侧局部水位下降，最大水位下降约 2.54 cm。2）流速变化分析。工程前后采样点流速变化值见表 8。根据数学模型计算结果可知，颜家台泵站下移 3 km 后实施清淤疏浚，在多年平均流量下，工程河段流速变化增加值最大为 0.16 m/s，流速减小值最大 0.22 m/s，流速变化 0.1 m/s 范围局限在清淤疏浚区上游 800 m 至下游 400 m，流速变化最大影响范围均在清淤疏浚范围附近。2援猿援猿颜家台泵站下移后产生的叠加效应为了进一步研究距离较近的水利工程取水口清淤对所在河段水文特征值所造成的叠加影响效应，在上节假定的基础下，仍选取多年平均流量的工况对颜家台泵站清淤后的河流模拟计算，不同的是本节以下移的颜家台泵站和普济水厂取水口同时清淤所产生的相互影响，来研究两个距离较近的取水口清淤所造成的叠加效应。1）水位影响叠加效应。上述两处水利工程取水口实施清淤后引起的水位变化见表 9。由表 怨可知：多年平均流量下，清淤疏浚工程实施后，清淤区上游水位略有降低，清淤区下游局部水位略升高。由模型计算结果可知，在多年平均流量条件入口流量工程前工程后左汊流量/（m3 s-1）分流比/%右汊流量/（m3 s-1）分流比/%左汊流量/（m3 s-1）分流比/%右汊流量/（m3 s-1）分流比/%6 000006 0000006 00009 400404.474.308 995.5395.70405.234.318 993.9395.6812 1001 639.4813.5510 460.5286.451 648.2813.6210 441.4986.2953 70020 831.6838.7932 868.3261.2120 833.6438.8032 863.1661.20孙光宝，刘红奎，李宗骏：河道清淤工程对防洪的累积影响研究67窑窑2023.10清淤区域洪水组次流量/（m3 s-1）流速最大增加值/（m s-1）流速最大减小值/（m s-1）颜家台泵站（下移）多年平均12 1000.160.22普济水厂多年平均12 1000.170.24颜家台泵站（下移）+普济水厂多年平均12 1000.190.28表 10 不同计算方案下取水区域流速最大影响值下，上述两个取水口同步清淤，造成水位最大壅高值 4.15 cm，大于任一取水口单独清淤造成的水位壅高值（3.87 cm；2.56 cm），小于两个取水口清淤造成的水位壅高值之和（4.15约3.87+2.56）。与此同时，两个取水口同步清淤，造成淤疏浚区上游侧局部水位下降，最大水位下降约 2.68 cm，大于任一取水口单独清淤造成的水位壅高值（2.54 cm；2.51 cm），小于两个取水口清淤造成的水位下降值之和（2.68约2.54+2.51）。由上述河流模拟计算分析可以看出，颜家台泵站（下移）和普济水厂清淤对清淤区在上下游引起的水位壅高和降低产生了叠加效应。2）流速变化分析。上述两处取水口清淤后的引起的流速变化见表 10。由表 10 可知：清淤疏浚工程实施后，工程河段流速变化增加值最大为0.19 cm/s，流速减小值最大 0.28 cm/s。由模型计算结果可知，在多年平均流量条件下，清淤疏浚区上游和下游局部流速增大，颜家台泵站（下移）和普济水厂取水口同步实施清淤后造成上游和下游局部流速最大增加值为 0.19 m/s，大于任一取水口单独清淤造成的流速最大增加值（0.16 m/s；0.17 m/s），小于两个取水口清淤造成的流速增加值之和（0.19约0.16+0.17）。上述两个取水口同步清淤，造成清淤疏浚前沿区域流速减小，最大流速减小值为 0.28 m/s，减小值大于任一取水口单独清淤造成的水位壅高值（0.22 m/s；0.24 m/s），小于两个取水口清淤造成的水位下降值之和（0.28约0.22+0.24）。由上述河流模拟计算分析可以看出，颜家台泵站（下移）和普济水厂清淤对清淤区上游和下游局部地区以及清淤疏浚前沿区域引起的流速增大或减小产生了叠加效应。3结 论1）清淤工程实施后，局部过水断面面积扩大，无阻洪建筑物构建，对局部行洪有利。清淤疏浚作业后，工程区域上下游附近水域的水位与流速变化幅度较小，且影响范围局限在疏浚区及其上下游局部河段范围。清淤施工不涉及堤身护护岸，对堤防不产生直接影响，不会引起崩岸或切滩等剧烈的河势变化，对河势稳定没有很大的影响，开挖区至护岸脚为近期堆积的中细沙，清淤对岸坡稳定也不产生不利影响。2）构建了二维水动力数学模型，利用 2019年 7 月 24 日（26 600 m3/s）和 2019 年 7 月 21 日（21 000 m3/s）流速和水位资料进行模型率定和验证，得出的结果表明断面流速计算值和实测值吻合较好，流速偏差在依0.1 m/s 以内。3）利用构建的二维水动力模型，得出结果表明，相邻取水口同步清淤对河道水位和流速的影响大于其中任一取水口单独清淤所产生的影响，但小于两个取水口单独清淤所产生影响的线性之和。参考文献1王光谦.长江中下游特大洪水分析J.群言，1998，（10）：15-17.2高耶.三峡工程运行后荆江三口与洞庭湖的水沙变化D.长沙：湖南农业大学，2019.3何杰，徐志扬，辛文杰.浅水方程 Roe 型格式的平衡性J.河海大学学报（自然科学版），2009，37（4）：450-456.4雷云龙.三峡蓄水后太平口水道滩槽演变机理及趋势研究D.杭州：浙江大学，2019.5赵建锋，李国斌，尚倩倩.长江下游江心洲河段航道整治工程防洪评价J.人民长江，2013，44（S1）：73-75，154.6张细兵，余新明，金琨.桥渡壅水对河道水位流场影响二维数值模拟J.人民长江，2003，（4）：23-24，40-48.7石波堂，唐昌焕，黎前查.浅析荆岳长江公路大桥对河势及防洪的影响J.中国农村水利水电，2007，（7）：99-101.（责任编辑 赵其芬）表 9 颜家台泵站渊下移冤和普济水厂同时清淤水位最大影响值清淤区域洪水组次流量/（m3 s-1）水位最大壅高/cm水位最大降低/cm颜家台泵站（下移）多年平均12 1003.872.54普济水厂多年平均12 1002.562.51颜家台泵站（下移）+普济水厂多年平均12 1004.152.68孙光宝，刘红奎，李宗骏：河道清淤工程对防洪的累积影响研究68窑窑