建桥对钱塘江涌潮影响的数值分析.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 建桥对钱塘江涌潮影响的数值分析 鲁海燕 潘存鸿 曾剑 （浙江省水利河口研究院 杭州, 310020） 摘要：本研究应用基于Boltzmann 方程的KFVS （Kinetic Flux Vector Splitting ）格式，建 立了富阳-澉浦的平面二维涌潮数学模型，在实测资料验证的基础上，计算分析了江东大桥建 设对涌潮高度、涌潮流速、涌潮传播速度等的影响。计算结果表明，建桥对桥位下游的涌潮 几乎没有影响，上游侧离桥位越近，影响越大，且对主槽内涌潮的影响程度大于边滩。桥位 近区涌潮高度减小的绝对值一般在10cm 以内，流速减幅最大达10%，桥轴线断面桥跨之间涌 潮流速增加0%~5% ，对涌潮传播速度的影响较小。 关键词：涌潮；KFVS 格式；桥梁；数值模拟 1 前言 钱塘江河口桥梁工程在满足南北两岸往来交通需求，加快实现长江三角洲地区经济交通 一体化发展中起着至关重要的作用，目前钱塘江河口闻家堰以下已建和在建跨江桥梁有8 座， 根据浙江省交通规划，沿江还将兴建 3 座。建桥后桥墩占据河道过水面积，并增加了水流局 部阻力，从而对涌潮产生一定程度的影响。钱塘江涌潮是我国乃至世界得天独厚的自然奇观 和极具特色的旅游资源，目前，有关主管部门已明确指出，尽管涌潮破坏力极大，但是兴建 涉水工程削减、消灭涌潮，也是一种不可挽回的资源损失。 [1] 曾剑等 应用钱塘江河口定床实体模型，模拟了涌潮的形成与传播过程，研究了单座桥梁 工程对涌潮的形成、形态与高度的影响，并根据模型试验数据建立了涌潮高度减小率的经验 [2] 公式。孙志林等 采用Delft3D 数学模型就单桥、多桥不同组合方案对钱塘江河口流场进行系 列数值试验，探讨了潮流对多桥的复杂响应。潘存鸿等[3]应用基于有结构网格的Godunov 格式 建立了二维涌潮数学模型，预测了萧山通道（钱江十桥）桥位断面的设计流速和设计流量， 但应用数学模型研究桥梁工程对涌潮的影响未见公开文献报道。 本研究应用基于Boltzmann 方程的KFVS （Kinetic Flux Vector Splitting）格式数学模 型[4]研究了江东大桥（钱塘江九桥）建成后对涌潮的影响，工程位于钱塘江河口段的中游，大 桥桥位如图1，跨江大桥长2 100m，有两个主跨和一个副跨，总共20 个桥墩，主跨为256m， 副跨为 160m，其余跨径为75m，85m，91.5m 不等。主墩承台高程为-1m，其他为0m，桥梁的 承台、桥墩柱的阻水面积约占1%频率高潮位下桥位断面过水面积的4.3%。该河段受径流和潮 汐的双重作用，河床冲淤变化剧烈，水动力条件非常复杂，桥位处多年平均潮差为1.55m，附 近河段涌潮可能最大高度约3.0m。 - 234 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 图1 钱塘江河口形势及潮位站、流速测点位置图 2 数学模型及其验证 平面二维涌潮数学模型的范围为富阳~澉浦，计算域内共布置42789 个三角形单元，21940 个节点，在桥附近局部加密，最小空间步长为2m，网格布置见图2。 图2 大范围网格图 模型对2000 年9 月、2005 年6 月、2006 年4 月及2007 年 10 月共4 次水文实测资料进 行了验证，其中 2007 年 10 月在江东大桥附近设有流速点，沿程潮位站及各年流速测点位置 如图1。验证计算所采用的地形及上、下游水边界与所验证的水文实测时期对应一致，计算的 时间步长为 0.1s，糙率系数按涨、落潮分别给定，上游糙率大于下游糙率，涨潮糙率取值 0.004~0.015，落潮糙率取值0.006~0.020。 图3 为2007 年 10 月仓前、盐官两站的潮位过程验证图，图4 为位于桥位附近的测点流 速、流向过程验证图。从图中可看出，潮