王茂枚，朱 昊，赵 钢，徐 毅，蔡 军

（江苏省水利科学研究院，南京 210017）

0 引 言

长江河势急剧变化区通常指水深、流急、水动力条件极其复杂，造成河床、洲滩摆动性大，江岸冲淤变动频繁的河段［1］。据有关资料，自20 世纪90年代以来，长江中下游河道已发生数百起崩岸现象，最严重时一年内发生数十起，长江中下游干流河道两岸岸线总长4 250 km，崩岸段的长度就达到1 520 km，占比达35.7%。新中国成立之后，随着水利事业的蓬勃发展，长江保护和治理受到高度重视，持续开展堤岸防护、河势控导以及应急抢护等工作，长江下游实施了大量河道整治工程，大大减少了崩岸现象的发生，大幅度提高了江岸的抗洪能力。以张家港老海坝段为例，根据2014年以来逐月河势监测资料表明，冲淤最大深度的月际变化超过10 m，河势变化极为剧烈［2］。水下抛石护岸能适应各种岸坡地形需要，而且造价低，施工方便，在长江中下游护岸工程中被广泛采用，也是长江河势急剧变化区治理的主要护岸型式［3，4］，研究抛石护岸工程的岸坡防护效果和对河势的控制效果对维持岸线稳定、保障防洪安全具有重要意义［5，6］。

国内外对抛石护岸的研究，主要从抛石护岸的稳定性、防冲促淤效果和河势控制效果三个方面进行考虑。抛石护岸的稳定性受水流条件影响显著，相关研究主要包括抛石粒径［7］、抛石厚度［8］、抛投落距［9］、施工技术［10］等。一些学者通过水槽实验方法，分析了块石的起动流速以及不同铺石厚度和面积对河床的防冲促淤作用，对抛石护岸的破坏机理进行了深入研究［11-13］。在河势影响方面，常用的研究方法是根据历年水文泥沙及水下地形等基础资料［14，15］，分析抛石护岸工程后河势的变化。Islam［16］基于卫星遥感图像对印度孟加拉河上的Farakka 大坝实施前后引起的河道弯曲处的曲率变化进行了研究，其研究成果在工程后期得到验证。赵钢［17，18］等基于点云数据分析对抛石护岸效果的抛投准确性、抛投均匀性、抛石增厚值等方面开展了水下抛石工艺抛投效果分析。罗青［19，20］等基于GIS 空间分析技术从抛石施工质量、施工过程及施工成本等方面对水下抛石效果进行评价研究。随着数值模拟技术的发展，其应用价值大幅提升，李大鸣［21］等提出对网格进行错位计算的方法，改进了相同层网格计算误差较大的缺点，并结合水平有限元和垂向有限差分的分层方法建立了三维水流泥沙数学模型，在海河下游进行了模拟，模拟精度较好。姜果等［22］采用平面二维水沙数值模型分析了抛石工程对水位、流速的影响，为河道治理提供了理论依据，然而其未考虑感潮河段涨潮和落潮流情形下抛石护岸工程影响的差异。

抛石护岸作为国内外应用最广泛的护岸型式，目前国内外研究主要集中于抛石施工工艺、施工方法、工程设计和破坏机理等方面，以往关于河势急剧变化区的抛石护岸防护效果的研究较少，且主要根据实测资料的结果进行河势演变分析，方法较为单一。此外，由于水下地形监测需要较大人力、物力的投入，监测时间受到较大的限制，目前仍难以做到实时水下地形监测。针对以上研究中存在的不足，本文以长江老海坝河段为例，基于Delft3D 开源软件建立三维水流数值模型，对抛石护岸实施后工程河段水流条件变化进行研究，对比分析抛石护岸工程实施后不同涨、落潮情形下河势急剧变化区水流结构的变化，为类似区域的长期治理提供参考依据。

1 研究区域概况

1.1 河道概况

本文研究区域为长江澄通河段（图1），上起江阴市鹅鼻嘴，下迄常熟市徐六泾，总长度约96.8 km，由福姜沙汊道、如皋沙群段、通州沙汊道等组成，属于长江下游感潮河段，江中沙洲众多，在径流和潮汐的双重作用下，河势演变剧烈。福姜沙河道上段顺直单一，下段为向南弯曲的稳定性较好的双分汊河道；如皋沙群段主要由民主沙、长青沙和横港沙等组成，如皋中汊及浏海沙水道上段的两股水流经民主沙汇合后进入浏海沙水道下段（长青沙南汊）；通州沙水道西起十三圩，下至徐六泾，随着通州沙整治工程的实施，目前向通州沙东水道和西水道分流的稳定双分汊河道发展。老海坝所在如皋沙群段位于长江口澄通河段中上部，处于河口潮汐与径流的交汇区，长约28 km，西起青龙港、东至天生港，该河段水流流路分歧多变，属于典型的多汊型河段。河床对水流约束的边界条件宽，使得水流随着流量的变化而动荡，岸线坍塌和淤积反反复复，变化较大。

抛石工程区域位于福姜沙南汊浏海沙水道下段右岸（见图1 方框）。20 世纪70年代初，主流进入浏海沙水道，直接顶冲右江岸老海坝一带，造成老海坝江岸全线崩塌。九龙港以下岸段因在顶冲区以下，受水流直接冲刷作用较弱，河岸与上游相比较为稳定。随后，如皋中汊分流开始增加，到70年代末中汊分流比已接近10%，浏海沙水道主流顶冲点由以往老海坝上段下移到九龙港～十一圩，因此九龙港以下近岸河床冲刷强度明显增强，江岸堤防安全受到直接影响。1990年以来，如皋中汊的发展趋于缓和，分流比渐渐趋于稳定，但是九龙港～十一圩主流顶冲区域的河床冲刷强度依然较大，河床被普遍刷深，深泓右移，前沿由-30 m 发展到-50 m 的贯通深潭，多处出现-60 m深潭。2014年逐月河势监测资料表明，张家港市老海坝河段河势变化极为剧烈，个别位置冲淤厚度的月际变化超过10 m，岸坡已处于临界状态。

图1 研究河段河势图Fig.1 River regime map of study reach

1.2 工程概况

张家港老海坝节点综合整治工程位于长江下游澄通河段浏海沙水道右岸，隶属长江澄通河段中部如皋沙群段，工程范围为张家港市一干河～九龙港～二干河以下1 260 m 之间的岸线（图2），里程总长约7 250 m，沿岸线走势长度约6 800 m。其中平顺护岸工程抛石分为堤防外坡的散抛石（外围堤外由码头平台及栈桥围成的内侧水域抛石）与码头外部前沿深水区抛石（散抛大块石、散抛小粒径块石防崩层）。工程分两期实施，一期工程抛石量约为96.4 万m2，二期工程抛石总量约147.8 万m2，散抛大粒径块石的粒径范围为0.4～0.65 m，散抛小粒径块石的粒径范围为0.16～0.4 m。水下抛石主要功能是提高岸坡的抗冲刷能力，维持岸坡稳定，为保证工程效果，设计抛石宽度为40～120 m，设计最大抛石厚度2.5 m。

图2 老海坝抛石护岸工程范围示意图Fig.2 Scope of Laohaiba riprap revetment project

2 数学模型的建立

2.1 模型计算范围及网格划分

本文主要基于软件Delft3D 进行研究，其水动力模型控制方程为基于浅水假设以及Boussines 假设的不可压缩流Navier-Stoke 方程［23］，通过合适的边界条件与控制方程组成定解条件，方程在空间上采用有限差分法进行离散、在时间上采用欧拉向前差分法离散。求解变量在网格中交错布置，采用ADI 交替隐式方法进行求解。

模型计算区域网格如图3所示，上游边界为江阴，下游边界至徐六泾，横向网格673 个，纵向网格168，对老海坝抛石工程区进行加密，最小网格尺寸为10 m 左右。计算时间步长为0.3 min，以确保柯朗数（Courant number）在计算过程中不超过1.0。模型采用2014年实测地形资料作为模型初始地形。平面坐标系采用北京1954 投影坐标，地形高程及计算潮位均采用1985国家高程基准。

图3 模型网格示意图Fig.3 Model grid diagram

2.2 模型验证

模型水动力验证采用2014年7月27日-7月31日的实测潮位、流速和流向资料（测站位置见图4），潮位验证选择如皋港、太子圩港、天生港、五干河、七干河、望虞河、营船港、汇丰码头等8 个潮位站的实测潮位过程，流速验证采用计算区域内的四条垂线的实测平均流速过程。验证结果如图5所示。模拟计算的水位、流速、流向的幅值和相位与实测值吻合较好，因此模拟结果可信，模型在模拟水流运动时是较为合理和准确的。

图4 研究区域水尺、测流点、采样点布置Fig.4 Layout of water ruler and flow measuring point in study area

图5 监测点潮位、流速及流向验证Fig.5 Layout of water ruler and flow measuring point in study area

2.3 抛石区糙率修正

由于河床演变在汛期的演变速率明显高于非汛期。为综合分析抛石护岸工程实施对研究河段水流条件的影响，选择汛期实测潮位模拟典型水动力边界下河段水动力过程。抛石护岸工程通过加密局部地形网格（图3）后抬高地形以及抛石区修正局部糙率的方法进行模拟。护底加糙后糙率与抛石前糙率关系采用以下公式［24］：

式中：n为加糙后的糙率系数；n0为加糙前的糙率系数；D为加糙后抛石中值粒径；D0为加糙前床沙粒径；k、k0为经验系数。

3 计算结果

3.1 对潮位的影响

为详细分析工程对水流特性产生的影响，在不同河段内选取多个采样点，位置如图4所示，其中P1-P13位于浏海沙水道，P14、P15 位于如皋中汊，P16～P18 位于如皋左汊。潮位变化主要体现的是抛石护岸工程对河段防洪、排涝的影响，因此重点对洪季条件下工程对潮位的影响进行分析。工程河段受非正规半日浅海潮影响，水位呈周期性涨落变化。为分析工程潮位变化对河段防洪、排涝的影响，选取模拟期内涨憩（潮水涨到最高时）和落憩（潮水落到最低时）的两个时间点进行工程前后的水位变化分析。

工程后涨憩、落憩时刻水位变化分布如图6 所示。由图可见，涨憩时刻抛石护岸工程实施后，离工程区越远水位变化幅度越小，整体呈上游水位壅高下游水位下降的规律。落憩时刻与涨憩时刻水位变化规律相似且变化规律更为明显，以一干河到十一圩港为分界，上游水位壅高下游水位降低。总体来看工程实施对河段的水位影响较小，不论是涨憩还是落憩时刻水位变幅大多在0.01 m以内。表1给出了工程河段内监测点水位的变化结果，由表可见，所有监测点水位变化值都在0.01 m 以内，抛石工程区刚好为水位壅高与水位降低的分界点，因此上游抛石区P7 水位壅高，下游抛石区P12 水位下降，而P10 则在涨憩时刻水位壅高，落憩时刻水位下降。其余监测点也符合上游测点水位壅高，下游测点水位下降的规律。总体上，落憩时刻水位变化幅度略高于涨憩时刻水位变幅，这可能是由于涨潮时，水流在径流和潮汐的共同作用下应抛石区地形拔高而产生的影响有些许抵消作用，而落潮时，河道内水流受径流作用，呈现较为规律的水位变化情势。

表1 监测点潮位变化统计 mTab.1 Statistics of tidal level changes at monitoring points

图6 工程后水位变化分布图Fig.6 Distribution map of water level change after engineering

3.2 对流速的影响

选取工程实施后河段涨、落急时刻垂向平均流速进行分析，工程实施后涨急、落急情形下流速变化分布如图7所示。模拟结果表明，工程实施后涨急情形下，流速变化主要发生在抛石工程区所在河段，即浏海沙水道从上游一干河至下游西界港所在区域，最大变幅在0.2 m/s 以内，其余河段流速变化均小于0.01 m/s。流速增大和减少区域无明显规律，究其原因，抛石护岸工程不仅增加了抛石区的河底高程，还改变了河道及岸坡的糙率，加上涨潮流与径流相互作用导致该处流速流向较为复杂，总体表现为抛石区及岸坡流速减小，抛石区外侧流速增大，深槽内流速又减少的规律。落急情形下工程实施后，流速变化区域与涨急时刻类似，主要局限抛石区及下游10 km 范围河段内，流速变化大多在0.1 m/s以内，离工程区越远流速变化越小。流速变化规律较为明显，抛石区及其下游带状区域流速减小，而抛石区外侧河道内则流速增大，这是由于抛石施工增加了河底高程，且块石增大了河道的糙率，从而增大了水流经过的阻力，使得一部分水流从左侧河道通过，增大了浏海沙水道左侧河道的流速。

图7 工程后流速变化分布图Fig.7 Distribution map of velocity change after engineering

为了更加直观的了解工程附近的流速变化情况，表2 给出了监测点工程前后的涨急和落急的流速大小，并给出了其差值大小及变化率。由表可知，涨急情形下，工程区以外各测点流速变化很小，变幅小于0.005 m/s，工程区内流速变化在-0.15～0.2 m/s，最大变幅发生在P10 测点，流速增幅达43.52%。落急时刻，非工程区流速变化较涨急时刻有所增大，整体流速变化在0.02 m/s以内，浏海沙水道及如皋中汊流速以减少为主，如皋左汊流速有所增加。工程区内流速变化较大，流速变化范围为-0.3～0.1 m/s，最大变幅发生在P12测点，流速减幅达31.51%。综上，工程实施后，不论是涨急时刻还是落急时刻，流速变化主要集中于抛石工程区域，流速最大变幅在0.3 m/s 左右，主要呈现抛石区流速减少深槽流速增大的规律。

3.3 对分流比的影响

对于分汊河道来说，分流比是反映汊道兴衰变化的重要水动力学指标，工程实施后河段的汊道形态、汊道阻力、河槽容积等均有所改变，影响汊道分流比，进而影响本河段的河势变化，为分析工程的实施对河道分流比的影响，在工程附近河道布置断面如图8。

图8 断面布置图Fig.8 Section layout diagram

工程前后各断面分流比统计见表3，民主沙将福姜沙水道分为如皋中汊（断面1）和浏海沙水道（断面2），其中浏海沙水道落潮时分流比占比超70%，落潮时刻占比大于涨潮时刻，工程前后两河道的涨落潮分流比有所变化。如皋中汊分流比在涨急时增加0.84%，在落急时增加0.02%，浏海沙水道分流比变化则相反。长青沙的存在使如皋中汊的水流一部分从如皋左汊流过，两股水流合并为浏海沙水道（断面4），与如皋左汊（断面3）共同流入通州沙水道，其中浏海沙水道落潮流时水流占比超98%，落急时刻占比大于涨急时刻，工程后断面4分流比略微减小，断面3略微增大。总体而言，工程的实施阻碍了浏海沙水道的水流，使得断面2 和断面4 分流比有所减小，断面1 和断面3分流比增大。同时，涨潮流所受工程影响大于落潮时刻，特别是断面1分流比变化最大，变幅达0.84%。

表3 工程后各断面分流比变化 %Tab.3 Change of diversion ratio of each section after engineering

4 结 论

基于Delft3D 采用三维水流数学模型针对长江下游河势急剧变化区抛石护岸工程对水流条件影响开展研究。分别针对落憩时刻和涨憩时刻工程实施后潮位变化进行分析，结果表明，落憩时呈现工程区上游壅高下游降低的规律变化，涨憩时由于潮汐的影响潮位呈不对称波动式变化；针对工程后涨急、落急时刻流速及分流比变化可知，流速变化主要集中于抛石工程区，呈现抛石区流速减少深槽流速增大的规律，工程实施后浏海沙水道分流比减小且涨急时刻所受影响大于落急时刻。所得成果为分析和预测抛石护岸工程实施效果及未来抛石方案的制定提供参考依据，也对保障河道防洪安全和堤防安全具有重要意义。 □