长江口北支咸潮倒灌工程方案的水流数值模拟分析.docxVIP

长江口北支咸潮倒灌工程方案的水流数值模拟分析 长江口北支是长江口三级支流和第四海水支流的初级支流。从西到东，在海口和兴港之间蜿蜒而行，全长约83公里。北支的上、下半部分向东弯曲，湾顶在大红港和大新河之间，下段呈喇叭形。中下段河槽宽浅,江中暗沙较多,受河口涨落潮动力条件、上游河道演变和人类活动的影响,自20世纪50年代末起,北支由落潮槽演变为涨潮槽。近几十年来,河道呈淤积萎缩趋势,河槽日益淤浅,容积不断减小,大量水、沙、盐倒灌南支,造成南支水质变咸和河势不稳定。因此,有必要采取工程措施来控制北支咸潮倒灌,消除或减轻对长江口淡水资源开发利用的不良影响,以提高上海市对长江口淡水资源开发利用的能力。 本文应用美国陆军工程兵团编制的数学模型软件SMS8.0,针对北支下段不同束窄程度的大、中、小3个控制工程方案进行模拟计算,分析不同工程方案实施前后潮量、流速及潮位等变化,比较不同方案的实施效果,为长江口北支整治提供科学依据。 1 数学模型 1.1 地球角速度的变化 二维水流数学模型采用的基本方程是沿深度积分的平面二维质量与动量守恒方程: 式中:h为水深;u,v为直角坐标系下的速度;x,y,t为直角坐标及时间;ρ为水的密度;E为涡动黏滞系数;g为重力加速度;a为底部高程;n为曼宁糙率系数;ω为地球角速度;φ为纬度。 1.2 边形、三角形单元 上述方程组采用伽辽金加权余量有限元法进行数值求解。单元的形式可为一维线性单元、二维四边形或三角形单元,适用于岸线曲折的边界。速度形函数为二次,深度形函数为线形。空间积分采用高斯积分,时间离散利用非线性有限差分。在每一时间步上,变量表达式为 式中:a,b,c是常数,c=1.5。方程采用全隐式求解,解法为Newton-Raphson非线性迭代。 1.3 水下地形模型 本模型的计算范围上边界取至长江下游的江阴,下边界取至东经123°,北边界取至吕四港附近(北纬32°15′),南边界取至芦潮港附近(北纬30°53′)。计算网格采用六节点二次三角形单元和八节点二次四边形单元相结合,空间步长100～3 000m,时间步长取15min。模型采用1998～2002年长江口实测水下地形资料。 模型中较为重要的参数曼宁系数和涡动黏滞系数由以下形式给出。 a. 曼宁系数。模型中曼宁系数n考虑有、无植被影响,根据各单元的平均水深按公式(5)自动赋值并实时调整。 式中:h为单元的平均水深;n0为无植被区最大的曼宁系数;n′为有植被区的曼宁系数;h0为植被影响糙率的水域水深;γ为深度系数。n取值范围为0.010～0.018。 b. 涡动黏滞系数。模型中涡动黏滞系数E根据各单元的平均速度和水流方向上单元长度自动赋值并实时调整,计算公式为 式中:P为Peclet常数;ρ1为液体密度;u为单元平均速度;dx为主流方向的单元长度。 1.4 水文资料及验证 本次计算模型的验证采用2002年9月洪季大潮的同步实测水文资料进行率定,采用2002年3月枯季的水文资料进行验证计算。率定和验证的结果表明,模型所采用的计算参数合理,结果合理、可靠,可以较真实地反映长江口地区的流场特性,可以模拟控制工程方案。 2 工程方案和影响分析 2.1 底段圈围 根据长江口综合整治研究规划,本文模拟了北支束窄控制工程方案,有大、中、小3个方案,如图1所示。大束窄方案:封堵新村沙右汊,北支下段沿新隆沙暗沙群-2 m线圈围,圈围崇明东滩,适当圈围连兴港下游滩地,口门宽度由11 km缩窄为5.5 km;中束窄方案:不封堵新村沙右汊,北支下段沿省市边界线适当后退圈围,不圈围崇明东滩,适当圈围连兴港下游滩地,口门宽度缩窄为6 km;小束窄方案:封堵新村沙右汊,北支下段沿新隆沙暗沙群-2 m线圈围,圈围崇明东滩,适当圈围连兴港下游滩地,口门宽度缩窄为9.2 km。 2.2 计算水文条件 上游江阴边界选用2002年3月枯季和9月洪季的实测水位过程,外海边界采用同期预报潮位过程。 2.3 不同项目的比较和影响分析 2.3.1 枯季涨潮量变化 北支是涨潮流占优的河道,造成北支咸潮倒灌南支一般有两个主要因素:一是上游来流量小,二是北支中下段呈喇叭口状,易形成涌潮且潮差较大。北支整治束窄方案的直接目的就是减小河道放宽率,控制进潮量,以达到减小潮差、减轻咸潮倒灌南支的目的。本文统计了北支进口断面连兴港断面的洪枯季涨潮量,大、中、小束窄方案在连兴港断面洪季涨潮量分别减少了37.3%,30.7%和27.8%;枯季分别减小了36%,29%和26.1%。可见,束窄方案较大幅度地减少了北支涨潮量,可以从潮量上控制北支咸潮倒灌南支。 束窄方案在减少外海咸潮来量的同时,北支的涨潮流动力作用也有所减弱,这可从崇头断面涨潮量和南北支会潮点的下移清晰地反映出。崇头断面统计显示,枯季大、中、小束窄方案涨潮量分