中咀湾波浪传播的数学模型.docxVIP

中咀湾波浪传播的数学模型 1 防波堤堤前波和工程海域港内波浪要素分析 大陈岛中滩湾位于12152e12153，2828n2829。三面环山，西山居首阶段。这是一个独特的地形，可以防止船只。为充分发挥其优势,为当地渔民就近避风创造良好条件,台州市椒江区海洋与渔业局拟在该处新建大陈岛渔船避风港,在港区入口处建设防波堤,在中咀湾建造消浪设施,并沿海湾岸线修建护岸,完善该处的避风锚泊条件,使该处形成船只台风期间锚泊避风,平时可修整网具等具备多功能的场所。为了给工程设计提供科学依据,需要对防波堤堤前波高和工程海域港内波浪要素进行研究。 本文依据大陈海洋监测站多年风速、风向的实测资料用皮尔逊Ⅲ型曲线得到3个主风向不同重现期的设计风速,进而根据浅水风浪成长公式推算计算区域边界处的波浪要素,再用波浪数学模型计算极端高水位和设计高水位时3种波况分别在3组重现期时工程海域的波浪要素。详细过程介绍如下: 2 域中咀湾口门处推导 依据大陈海洋监测站26a风向和3个风向年最大风速实测资料进行频率分析,用皮尔逊Ⅲ型曲线推求3个风向(NW-NNW,W-WNW,SW-WSW)50a、25a和2a重现期的设计风速,见表1。 外海波浪要素根据工程海域情况选用《工程环境海洋学》中浅水风浪成长公式,自椒江口门向大陈工程海域中咀湾口门处推算求得,公式如下: gHsU2=5.5×10?3(gFU2)0.35tanh[30(gdU2)0.8(gFU2)0.35](1)gTsU=0.55(gFU2)0.233tanh2/3[30(gdU2)0.8(gFU2)0.35](2)gΗsU2=5.5×10-3(gFU2)0.35tanh[30(gdU2)0.8(gFU2)0.35](1)gΤsU=0.55(gFU2)0.233tanh2/3[30(gdU2)0.8(gFU2)0.35](2) 式中:d为计算水域的水深;计算时,有效波高Hs及风区长度F的单位为m,有效周期Ts的单位为s;U为海面上10m高处的风速,单位为m/s。 NW-NNW方向风区长度约为55km,W-WNW方向风区长度约为40km,SW-WSW方向风区长度约为90km,工程水域的水深范围为7～13m。由公式(1)和公式(2)推算出口门边界处的波浪要素,如表2所示,据调查与当地实际情况基本相符。 3 波浪传播的影响因素 早先的研究已经证明了用数学模型研究波浪的优越性,外海波浪在由深水向近岸或岛屿传播过程中以及当边界条件发生突变时,因折射、绕射、反射、浅水变形、波能损耗和波形破碎等波浪近岸变形的影响,波浪要素不断发生变化。为了考虑波浪传播中的诸多因素,需选用合理的波浪数学模型,计算波浪要素。目前世界上用数学模型研究外海波浪的方法主要有缓坡方程、Boussinesq方程、能量守恒方程等方法,本文采用考虑陡变地形和能量耗散效应的改进型缓坡方程为控制方程,结合固边界的反射边界条件,构成波浪传播变形的联合折射、绕射和反射的数学模型。本文所使用的模型已被一些学者用许多地形实测数据验证过,模型的计算结果与试验结果吻合很好,说明本模型能较好的模拟波浪在近岸的传播变形,反映了波浪在传播过程中的折射和绕射效应,适合近海地区的波浪传播模拟。 4 工程及港内波高分析 中咀湾的走向大致呈NE-SW,水域相对狭小,西深东浅。此海域西面毗邻大陆,所以外海混合浪主要从北向、东向和南向传来,由于中咀湾位于上大陈岛西南部的特殊地理位置,外海混合浪基本被上、下大陈二岛屏蔽,对湾内水域影响有限;而对中咀湾影响较大的波向主要有3个:NW-NNW,W-WNW,SW-WSW,这3个方向的波浪主要是当地局部风场产生的局地风浪。 计算区域以及防波堤位置如图1所示,坐标为北京54坐标,西南边界为开边界,距防波堤1.4km左右。网格空间步长为10m×10m,网格数为28624,底摩擦系数取0.01。防波堤长度为170m,实心式防波堤为斜坡式结构,透空式防波堤的透射系数为0.3～0.4,计算时各水位均统一到大陈海洋监测站平均海平面。本文对工程海域分别进行了3个重现期(50a、25a、2a)、3个波向(NW-NNW向、W-WNW向、SW-WSW向)和2种设计水位(极端高水位3.556m、设计高水位2.256m)组合条件下的波浪折射、绕射及反射所引起的波浪变形的计算。下面只给出NW-NNW向重现期为50a极端高水位情况下工程前以及工程后实心式防波堤和透空式防波堤三种情况下的工程区域局部波高分布图(见图2)。 中咀湾中受防波堤影响的水域面积达246 123m2,根据椒江区海洋与渔业局的要求,船舶停泊面积定义为波高1m的工程水域面积。从实心堤工程前与工程后的波高分布图可知,工程前波浪未受阻挡,一直可传播到中咀湾湾顶,而工程后,虽然波浪在堤外反射,能量累积,部分区域甚至出现立波