□ 韩永林（浙江省水利水电建筑监理公司）

0 引言

新泓口围堤方案为北直堤长1.70km，南直堤长1.10km，顺堤（东直堤）长3.60km，围海面积537hm2。二期管道镇海登陆点位于顺堤北部，距离北直堤约300m。围海工程于2008年至2009年建成，其建设对海堤附近的海底地形产生了影响。文章根据实测的地形和水文泥沙资料，分析该工程海域的水文泥沙和地形特征；采用潮流数学模型，结合泥沙冲淤半理论半经验公式，计算分析该工程建设对潮流和泥沙淤积的影响。

1 潮流数学模型

二维潮流的基本方程包括水连续性方程和动量守恒方程，即

式中:t为计算时间；x、y为计算区域坐标；H为海面到海底的距离；Zb为海底高程；U、V分别为垂线平均流速在x、y方向的分量；f为地转参数，g为重力加速度；τsx、τsy为海面风应力分量；τbx、τby为海底摩阻分量，ρ为海水密度；ve为水平涡动粘性系数。上述参数具体计算见文献[1]。

上述方程式是一个初边值问题，求解需要适当的初始条件和边界条件。

计算边界分为两类，一类是岸边界，另一类是水边界。初始条件与潮位和流速的初始分布有关。采用有限差分法离散基本方程。通过离散得到的线性代数方程组，用高斯消去法求解。

2 泥沙淤积计算

根据潮汐水流中的悬沙运动方程，可得海岸工程引起的床面冲淤强度，即

式中:p为年冲淤强度，m/a；α为沉降机率；ω为泥沙沉速，m/s；γd为泥沙干容重，kg/m3；t为一年中的冲淤历时，s；S1、S2为工程前、后的水流挟沙力，kg/m3。其中计算参考文献[3]和文献[4]

3 数值模拟分析

3.1 工程简介

本路由总体上呈NE－SW穿越杭州湾口灰鳖洋海域，路由海域的海底地形特征为西南部高，平缓东倾，东北部低且起伏。……