建桥对望虞河行洪及引水能力影响的计算分析

建桥对望虞河行洪及引水能力影响的计算分析

曲红玲，高祥宇，夏益民

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

摘要：沿江开发高等级公路跨越望虞河，通过平面二维水流数学模型计算分析拟建大桥对望虞河行洪和引水时流量、水位、流速的影响。结果表明：大桥修建后对河道行洪及引水能力的影响较小，桥墩周围水位略有变化，河道内流速影响范围不大。 关键词：行洪；引水；跨河大桥；数值模拟

望虞河是太湖流域骨干泄洪引水通道，对流域防洪、水资源配置和水环境保护有举足轻重的作用。新一轮太湖流域防洪规划，将进一步扩大望虞河行洪和引水能力，建成太湖行洪的高速通道和引江“清水走廊”[1]。沿江开发高等级公路在常熟海虞镇境内跨越望虞河，大桥对其行洪和引水能力影响研究十分必要。 1 工程概况

1.1建设项目概况[2]

沿江开发高等级公路望虞河大桥位于常熟市海虞镇境内，大桥距下游常熟枢纽节制闸约1 265 m ，距上游虞王桥约1 250 m 。桥轴线与水流方向交角为90°。图1为望虞河现状示意图和拟建大桥工程位置。

图1 望虞河现状示意图和拟建大桥工程位置

拟建大桥桥型为7×30+(55+85+55)+6×30预应力混凝土连续箱梁结构，主桥下部结构采用薄壁式墩，钻孔灌注桩基础。桥梁全长592.2 m 。大桥主墩左右各1排2个，顺流排列，长12 m 、宽2.5 m ，两端迎水面和背水面均为半圆形。

1.2工程所在河道情况[3]

望虞河通航等级五级，通航要求净宽50 m ，净高5 m 。桥位处河口宽约150 m ，河道断面标准为：底高程-3.0 m （吴淞高程）、底宽82 m 、坡比1:3~1:5（即高程-1.0 m 以上1:3，以下1:5）；东西岸堤防标准为：顶高程分别为6.0 m 、5.5 m ，顶宽分别为11.5 m 、4 m ，坡比1:2，青坎高程4.5 m ，青坎宽度分别为5 m 、2 m 。 2 平面二维水流数学模型

2.1控制方程及其离散求解

在流体力学中，浅水流动是对实流动的一种简化和概化的数学模型。本模型采用基于圣维南方程基本假设的浅水方程组：

()()()()E U G U U U F U S U t t x y

∂∂∂∂+∇=++=∂∂∂∂ （1） 式中：h U hu hv ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，F Ei Gj =+v v ，222hu gh E hu huv ⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎣⎦，222hv G huv gh hv ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥+⎣⎦，0()()ox fx oy fy S gh S S gh S S ⎡⎤⎢⎥=−⎢⎥⎢⎥−⎣⎦，h 为水深；u 为x 方向流速；

联系方式：xygao@

v 为y 方向流速；ox S 为x 方向底坡；oy S 为y 方向底坡；x

方向摩阻坡降f x S

=y 方向摩阻坡降3

4222h v u v n S fy +=；n 为曼宁粗糙系数。

[4] 对控制体积分方程（1），应用Gauss-Green 公式，化为沿其周界的线积分，得：

()x y s U d En Gn ds Sd t ΩΩ

∂Ω=++Ω∂∫∫∫ （2） 对于m 边凸多边形，

（2）式等号右边第一项可离散成m 项之和，在数值上等于被积函数在控制体各边上的法向值与该边长度的乘积，假定水力要素在各控制体内均匀分布，则式（2）可化为：

1()m i i i n n

i U A E G L AS t =∂=++∂∑ （3） 利用欧拉方程的旋转不变性，将式（3）变为：

111()()()m n n i i i A U U t T F U L AS θ+−=⎡⎤−=Δ+⎢⎥⎣⎦∑

（4）

式（4）左边表示控制体内守恒变量在t Δ内的变化，右边第一项表示沿第i 边法向输出的平均通量乘以相应边长，第二项表示控制体内源项（入流及外力）在t Δ内的作用；这反映了守恒物理量的守恒原理：守恒物理量在控制体内随时间的变化量等于各边法向数值通量的时间变化量和源项的时间变化量。式（4）中的输出通量()F U 采用通量差分裂格式（FDS ）求出[5,6]。

2.2定解条件

2.2.1 初始条件

计算时取初始流速00u =和00v =，初始水位0z 为一给定数。模型经过一定时间的运行，初始条件的影响将会消除。

2.2.2 边界条件

模型开边界用水位过程或流量过程控制；闭边界采用不可入条件，即取法向流速0n v =，n 为边界的外法线方向。

2.3计算范围、网格及地形

2.4糙率系数和时间步长

本次计算糙率采用曼宁公式，参考《水工设计手册》推荐的取值和有关望虞河类似大桥的洪水影响评

价成果，并考虑本河段的实际情况，糙率系数取值为0.025。

根据模型计算稳定条件及计算精度的要求，大模型的时间步长取为1 s，小模型的计算步长取为0.01 s。

2.5计算条件

本次计算选择1999年型（相当于100年一遇）特大洪水、1954年型（相当于50年一遇）设计洪水和

引水时望虞河的最大流量及相应水位作为模型边界控制条件。根据常熟枢纽调度情况，本次计算所采用的

具体水文条件见下表1。

表1 计算采用水文条件表

水文类型上边界条件（望亭立交闸下水位）下边界条件（常熟枢纽情况）备注

多年平均高水位4.00 m

百年一遇洪水 5.09 m

多年平均低水位1.56 m

上下边界都是水位控制

闸上实测水位过程线

50年一遇洪水 4.84 m 多年平均低水位1.56 m

引水 3.0 m 350 m3/s 上边界水位控制下边界流量控制表中上边界水文条件是指望亭立交处的水文条件，下边界水文条件是指常熟枢纽处水文条件。闸上实

测水位过程线是1999年7月13日5时52分~12时33分之间开闸泄洪期间水位过程线。

大模型计算结果采用1999年7月13号实测琳桥水位验证，当望亭立交闸下水位为5.09 m时，琳桥实

测水位4.73 m，计算水位为4.77 m，验证结果表明模型计算结果与实际情况基本相符。

3 拟建工程对流量、水位和流速的影响

3.1拟建大桥对河道过水面积变化分析

拟建大桥在望虞河河道内布置了两组薄壁桥墩（共4个），每个薄壁桥墩立柱顺水流方向12.0 m，垂

直水流方向2.5 m。根据二维水流模型计算结果，本工程布置在河道内的桥墩阻水面积比均低于3.5%，阻

水面积比较小。望虞河行洪和引水情况下桥墩阻水面积统计见表2。

表2 常熟西段望虞河大桥桥墩阻水面积统计表

水文条件桥址水位/m 河槽过水面积/m2阻水面积/m2阻水面积比/%

下边界高水位 4.02 804.78 26.25 3.26 百年

一遇下边界低水位 1.77 500.21 15.13 3.02

50年

下边界低水位 1.68 488.78 14.90 3.05 一遇

供水 3.77 773.76 25.05 3.24 3.2拟建大桥对望虞河行洪及引水流量影响分析

常熟西段望虞河大桥建设后，由于河道内桥墩的建设，使河道桥址处过水断面减小，加上桥墩阻水作用，会对望虞河过水流量有所影响，为了比较流量变化情况，在桥址所在断面计算工程前后的流量情况。计算结果表明：百年一遇洪水下边界高水位和低水位的过流量分别为516 m3/s和525 m3/s，50年一遇洪水下边界低水位的过流量为511 m3/s。引水时，断面过流量为350 m3/s。建桥后，断面过流量略有减小，百年一遇洪水下边界高水位和低水位的过流量分别减小1.3 m3/s和2.3 m3/s，50年一遇洪水下边界低水位的过流量减小2.1 m3/s。引水时，在工程前后上下游边界水位不变的情况下，断面过流量减小1 m3/s。

由此可见，常熟西段望虞河大桥建设后，对望虞河行洪和引水能力有所影响，但影响程度较小。

3.3拟建大桥对望虞河行洪和引水水位影响分析

常熟西段望虞河大桥建设后，由于桥墩阻水，使望虞河在行洪时桥上游产生壅水，桥下游水位有所降低；引水时，桥下游水位壅高，上游水位降低。为分析拟建大桥对望虞河行洪和引水水位影响，在河道内

桥上下游及桥孔之间共布置了29个采样点（见图5）进行工程前后水位变化（图6~9）比较。