防淤积型透空式防波堤水动力特性分析.pdf

第 44 卷第 8 期2023 年 8 月哈 尔 滨 工 程 大 学 学 报Journal of Harbin Engineering UniversityVol.44.8Aug.2023防淤积型透空式防波堤水动力特性分析于通顺1,唐俊辉2,孟晓宇1(1.中国海洋大学 工程学院,山东 青岛 266000;2.中国电力工程顾问集团东北电力设计院有限公司,吉林 长春 130021)摘 要:为探究新型防淤积型透空式防波堤的水动力性能,本文基于标准的 k-紊流模型构建了规则波与透空式防波堤相互作用的数值模型,开展了物理模型试验。本文利用数值模拟方法探究了规则波作用下透空式防波堤的消波性能、荷载特征等水动力特性。研究发现:挡板入水深度增大 50%,透射系数减小 9%28%;入水深度增大100%,透射系数减小 42%58%。防淤积型透空式防波堤可通过调节上拦沙板的下落高度保证潮位变化过程中的挡浪效果。上拦沙板的无量纲水平波浪力和上托力随波陡的增大而减小,在大波作用下,上拦沙板整体处于挡浪板后时其承受的水平波浪力更大。防淤积型透空式防波堤的特殊结构及其良好的消波效果,决定了其有成为理想的海岸防护结构的巨大潜力。关键词:透空式防波堤;规则波;数值模拟;消波性能;透射系数;荷载特征;波浪力DOI:10.11990/jheu.202106027网络出版地址:https:/ 文献标志码:A 文章编号:1006-7043(2023)08-1295-09Analysis of the hydrodynamic characteristics of an antideposition permeable breakwater YU Tongshun1,TANG Junhui2,MENG Xiaoyu1(1.College of Engineering,Ocean University of China,Qingdao 266000,China;2.Northeast Electric Power Design Institute Co.,Ltd.of China Power Engineering Consulting Group,Changchun 130021,China)Abstract:A numerical model of the interaction between regular waves and permeable breakwaters was constructed on the basis of the standard k-turbulence model to explore the hydrodynamic performance of a new antideposition permeable breakwater.Physical model tests were conducted to verify the accuracy of the numerical model.Numeri-cal simulation was used to explore the hydrodynamic characteristics,such as wave dissipation performance and load characteristics,of the permeable breakwater under regular waves.When the water entry depth of the baffle in-creased by 50%,the transmission coefficient decreased by 9%28%.When the water entry depth increased by 100%,the transmission coefficient decreased by 42%58%.The antideposition permeable breakwater can ensure the wave-blocking effect during the change in tidal levels by adjusting the drop height of the upper sand-board for sediment trapping.The dimensionless horizontal wave force and uplift force of the upper sand baffle,which de-creased with the increase in wave steepness,can ensure the wave-blocking effect during the change in tidal levels by adjusting the drop height of the upper board for sediment trapping.Under the action of large waves,the upper sand-blocking plate bore a great horizontal wave load when it was entirely behind the wave baffle.The special struc-ture of antideposition permeable breakwater and its good wave dissipation effect determine its great potential to be-come an ideal coastal protection structure.Keywords:permeable breakwater;regular wave;numerical simulation;wave attenuation performance;transmis-sion coefficient;load characteristics;wave force收稿日期:2021-06-08.网络出版日期:2023-05-06.基金项目:国家自然科学基金项目(51739010).作者简介:于通顺,男,教授,博士生导师.通信作者:于通顺,E-mail:tshyu707 .随着我国海洋生态文明建设的不断推进,对防波堤的环保属性提出了更高要求。传统堤身封闭式的防波堤阻断了港池内外水体交换,易引发港池内水质劣化1,已不适应“蓝色海湾、银色沙滩”的建设需要。为解决传统防波堤的环保问题,透空式防波堤2-3应运而生,其水动力性能受到国内外研究哈 尔 滨 工 程 大 学 学 报第 44 卷人员的广泛关注,对于竖直板式透空堤的理论推导做了大量的工作。Ursell4根据修正的贝赛尔函数推导了有限潜深的竖直薄板在波浪作用下透射系数的解析解。Weigel5提出了透空式防波堤的理念,在忽略了波浪的绕射的前提下通过微幅波理论推导了竖直挡板式防波堤的反射、透射系数的解析解。邱大洪等6将波浪的绕射现象考虑在内,推导在任意水深下单一薄板式防波堤的透射、反射系数解析解,为后来有关板式防波堤的研究奠定了重要的理论基础。琚烈红等7在忽略波浪与防波堤相互作用过程中能量损耗的前提下引入了经验系数,基于波浪理论和能量守恒定律探究单层挡板透空式防波堤的透射系数解析解。朱大同8在考虑流体粘性情况下对单层竖直挡板防波堤的消波特性和受力特性展开研究。9在 Wiegel 公式5基础上引入局部损失阻力系数,得到透射系数解析解,并给定了相应的适用条件,目前我国防波堤设计与施工规范中采用的公式与文献9推导一致。防波堤水动力性能的模拟研究是透空式防波堤研发的重要内容。殷福安10采用物理模型断面试验方法给出了特定波浪条件下单层竖直挡板式防波堤透射系数的公式。Neelamani 等11研究了由互相平行的 2 个竖直板组成的挡浪结构的水动力特性,分析 2 个板中间水体的波动,并基于结构物的消波效果,发现双侧挡板对不规则波的衰减作用较规则波更明显。文献12-13研究双侧挡板式透空式防波堤波浪要素、结构尺寸等对消波特性的影响。针对波浪作用过程中水质点做椭圆形运动的特点及透空式防波堤只对特定参数范围内的波浪有明显消波效果的特征,周效国等14设计一种多层直立挡板式防波堤,各挡板设置不同表面开孔率和底高程,其中下部挡板可以前后运动,对不同水位、波高和周期的波浪逐层消减。严以新等15对竖直多层挡板式挡浪结构进行物理模型试验,根据结果给出了设置 4层挡板以及将挡板进行前排稀疏后排加密排列的建议。Rao16将挡板斜向布置,对倾斜角度不同的板式结构进行了全面的物理模型试验研究。目前竖直挡板式防波堤的研究比较成熟,在国内已经成功得到应用17-19。然而港区频繁涨落潮过程中港池外泥沙会随水流大量进入港池内部,造成港池内泥沙淤积。针对这一问题,于通顺等20在竖直挡板式防波堤结构基础上,提出了一种防淤积型透空式防波堤。考虑到海水中泥沙主要集中在下层水体中21,该防波堤在竖向桩柱底端设置下拦沙板,挡板下表面与海床紧密接触,在上部挡浪板后侧设置上拦沙板,挡板可沿桩柱表面滑道上下移动。当海水中泥沙含量较低时,上拦沙板全部隐藏在挡浪板后,当泥沙含量超过设定的阈值时,上拦沙板会沿着滑道下降,与下拦沙板对接在一起,将港池内外过水通道完全阻隔,避免泥沙进入港池内。本文针对防淤积型透空式防波堤的水动力性能,采用数值模拟的方法开展研究。建立了波浪与防淤积型透空式防波堤相互作用的数值模型,开展物理模型试验以验证数值模型的准确性,基于数值模型考察了透空式防波堤的消波效果与受力特征。1 水动力分析数值理论 本文基于流体连续性方程、动量方程和标准的k-紊流模型构建波浪与防淤积型透空式防波堤相互作用的数值模型。连续性方程:x(uAx)+y(vAy)+z(wAz)=0(1)动量方程:ut+1VF(uAxux+vAyuy+wAzuz)=-1px+Gx+fx(2)vt+1VF(uAxvx+vAyvy+wAzvz)=-1py+Gy+fy(3)wt+1VF(uAxwx+vAywy+wAzwz)=-1pz+Gz+fz(4)式中:u、v、w 是 x、y、z 方向的流速;Ax、Ay、Az是 x、y、z方向控制方程的各个波浪面积分项;VF表示控制方程的可流动体积分数;为介质密度;Gx、Gy、Gz与fx、fy、fz为流体在三维空间中的重力以及粘滞力加速度;p 为压力。标准 k-紊流模型为:kt+ukx+vky=+tk()2kx2+2ky2()+1ktxkx+tyky()+2tux()2+y()2+tuy+x()2-(5)t+ux+vy=+t()2x2+2y2()+1txx+tyy()+2C1ktux()2+y()2+C1ktuy+x()2-C22k(6)6921第 8 期于通顺,等:防淤积型透空式防波堤水动力特性分析式中:k 为湍动能;为耗散率;为运动粘滞系数,t=Cuk2/为紊动粘性系数,标准的 k-模型中的参数选择为 C1=1.43,C2=1.92,k=1.0,=1.322。本文采用基于界面重构技术的自由表面处理方法(volume of fluid method,VOF)23捕捉自由液面,根据不同时刻流体在网格单元中所占体积函数 F的不同追踪自由液面的变化过程,可较好地捕捉自由液面破碎、冲击、翻滚等现象。2 数值模型的建立2.1 数值水槽建立 所建数值水槽长宽高分别为 18、0.6、0.7 m。造波机设置在水槽 x 轴方向 0 m 位置,末端设置消波网,模型中心位于水槽 x 轴方向 7 m 处,模型前设置的 2 个监测点分别在模型前 2 m 和 1.5 m 处,模型后设置的 3 个监测点分别在模型后 0.75、1 和 2 m 处,用以监测模型前后自由液面的变化情况。数值模型及监测点的布置情况如图 1 所示。图 1 波浪与防波堤相互作用数值模型示意Fig.1Numerical model diagram of interaction between wave and breakwater2.2 计算网格划分 将整个计算区域沿 x 轴方向分为 3 部分,第 1部分 0 6.7 m,第 2 部分 6.7 7.3 m,第 3 部分7.318 m,并对第 2 部分防波堤结构所在区域和 z方向波高附近网格进行加密,整个计算区域划分为15 块,如图 2 所示。其中,块 1、块 2 和块 3 的 x 方向网格尺寸为 0.02 m,块 4、块 5 和块 6 的 x 方向网格尺寸为 0.01 m,块 7、块 8 和块 9 的 x 方向网格尺寸为 0.005 m,计算区域 y 方向网格尺寸为 0.02 m,块 1、块 4 和块 7 的 z 方向网格尺寸为 0.01 m,块 3、块 6 和块 9 的 z 方向网格尺寸为 0.015 m。图 2 防波堤网格划分示意Fig.2 Schematic diagram of breakwater grid division2.3 模型尺寸及工况设计 防波堤设计与施工规范24规定,对基本不越浪的直立式防波堤,堤顶高程宜定在设计高水位以上不小于 1.25 倍设计波高处,当透空式防波堤采用单侧挡浪结构时,挡板入水深度与水深之比宜取0.30.5。基于广东海安新港荔枝湾码头25所在位置的水利条件和桩基透空式防波堤结构的建造尺寸,设计数值模拟的水利条件及防波堤的模型尺寸。模拟工况采用规则波,根据上拦沙板下落距离 D 的不同将模型分为 3 种结构形式。试验在水深 0.4 m与 0.45 m 的 2 种条件下进行,波高范围为 0.06 0.14 m,波浪周期为 1.01.4 s,上拦沙板下落范围为 00.10 m,组合不同的波浪参数共计 64 组工况。防波堤模型挡浪板的高度为 0.4 m,当水深为 0.4 m时,挡板初始入水深度为 0.1 m,上、下拦沙板的高度均为 0.15 m,挡板厚度为 0.02 m,方形桩柱边长为 0.06 m。2.4 网格独立性和模型准确性验证 为了更精确地获取波高数据,需对 z 方向波高附近网格进行加密,因此首先要进行网格独立性验证。在波高范围内设置 3 种网格数量方案,分别为20、40 和 60。在三维水槽模型中对水深 0.4 m,波高 0.08 m,周期 1.2 s 波浪工况进行计算,比较 3 种网格尺寸下相同位置波面的历时曲线,结果如图 3所示。比较后发现,3 种方案在 x-z 方向上的网格最大纵横比都为2,网格数量为40 与60 的2 种方案的计算结果基本吻合,仅在峰值和谷值处存在细微差别,网格数量 20 的方案与其他 2 种方案的计算结果差距较大。表 1 给出了不同网格设置方案下计算域的网格数量及求解器运行时间,综合考虑计算精度和时间成本,在波高范围内 z 方向划分的网格数量为 40,因此块 2、块 5 和块 8 在 z 方向网格尺寸为0.005 m。图 3 不同网格尺寸的波面历时曲线Fig.3 Wave duration curves of different mesh sizes提取 x=7 m 处的波高历时曲线,并将其与理论波形作对比,如图 4 所示。可以看出数值模拟的波浪与理论波形相比存在一定的变形,波峰处略尖陡,波谷处略平坦,在波面下降阶段,数值模拟结果相对理论值略小,但整体而言与理论吻合较好,表明本文所建立的数值水槽能产生持续稳定的规则波浪,数值水槽 z 方向液面附近网格划分合理。7921哈 尔 滨 工 程 大 学 学 报第 44 卷表 1 三维数值水槽流体域计算精度与时间成本Table 1Calculation accuracy and time cost of fluid do-main in three-dimensional numerical flume数量(波高范围)网格大小/cm计算域中的网格总数求解器运行时间/min200.011 152 000160400.0052 999 800424600.0036 624 200936图 4 波面历时曲线计算值与理论值比较Fig.4 Comparison between numerical and theoretical re-sults of wave surface为验证数值模型建立的准确性,在中国海洋大学工程水动力学实验室的波流水槽中,设计并进行了防淤积型透空式防波堤水动力特性的物理模型试验。依据重力相似准则,将防波堤的模型比尺定为1 10 进行试验,试验模型如图 5 所示。采用电容式波高仪进行波高数据采集。透空式防波堤的透射系数为:Kt=Ht/H(7)式中:Ht为透射波高;H 为入射波高。比较水深 h=0.4 m,波高 H=0.08 m,周期 T=1.0、1.2、1.4 和 1.6 s 的 4 组工况的数模结果和物理模型试验透射系数,如图 6 所示。可以看出:在 4组工况条件下,物理模型试验透射系数与数值模拟结果差距分别为 8.1%、8.1%、5.6%和 6.0%,因此可以认为建立的波浪与防波堤相互作用的模型是准确的。3 数值结果分析与讨论 建立模型时,在防淤积型透空式防波堤上拦沙板底面和前后表面布置多个压力测点,如图 7(a)和图 8(a)(c)所示。假设测点处的压强可以代表该点周围面积压强的平均值,则该测点附近所承受的波浪力可以表示为测点压强与面积的乘积,方向垂直作用于受力表面,对上拦沙板各表面所受的波浪力进行矢量合成,得到作用在挡板上总的波浪压力F,并用 gH3对其无量纲化,其中密度 取值1 000 kg/m3,重力加速度 g 取值 9.81 m/s2。图 5 波浪与防波堤相互作用模型试验示意Fig.5Model test diagram of interaction between wave and breakwater图 6 物理模型试验透射系数和数值结果对比Fig.6 Comparison of transmission coefficients of physical model tests and numerical results图 7 防波堤上拦沙板底面压力测点Fig.7 Pressure measuring point at the bottom of the upper board for sediment trapping on breakwater 设定防淤积型透空式防波堤模型 2 个竖向支撑桩柱之间一段完整防波堤的宽度为 40 cm,利用上拦沙板表面监测点压强计算挡板受力时,首先根据结构对称性选取上拦沙板各点的压强值,如图 7(b)所示。3.1 防波堤消波效果 防淤积型透空式防波堤透射系数随波陡的变化如图 9 所示,采用多项式回归方法对同一工况下数据点进行拟合,获得透射系数随波陡变化的二阶曲线。可以看出:透射系数随着波陡的增大8921第 8 期于通顺,等:防淤积型透空式防波堤水动力特性分析而减小,且减小速率逐渐变大。上拦沙板下落距离 D 越大,即挡板的相对入水深度越大,透射系数越小,且随着下落距离的增大,透射系数的减小速率加快。挡板的初始入水深度为 0.1 m,当上拦沙板下落 0.05 和 0.10 m 时,入水深度分别增大了50%和 100%。经比较,在 T=1.0 s 时,上拦沙板下落距离 D=0.05 m 和 D=0.10 m 比 D=0 m 时波浪的透射系数最高分别减少 28%和 58%,在 T=1.2 s 时,波浪的透射系数最高分别减少 15%和45%,在 T=1.4 s 时,波浪的透射系数最高分别减少 9%和 42%,在 T=1.6 s 时,波浪的透射系数最高分别减少 11%和 47%。图 8 防波堤上拦沙板侧面压力测点布置示意Fig.8 Schematic diagram of pressure measuring point arrangement on the side of the upper board for sediment trap-ping on breakwater图 9 透射系数随波陡变化关系Fig.9 Relation of transmission coefficient with wave steepness 防淤积型透空式防波堤透射系数随 H/gT2的变化如图 10 所示,可以看出:拟合的函数曲线可以9921哈 尔 滨 工 程 大 学 学 报第 44 卷较准确的预测新型透空式防波堤的透射系数。透射系数预测公式采用二阶多项式函数的形式,使公式透射系数的预测值与模拟结果之间的误差更小,可以更准确预测其他工况条件下防波堤的透射系数。图 10 透射系数随 H/gT2变化关系Fig.10 Relation of transmission coefficient with H/gT2不同水深条件,上拦沙板下落不同距离时,防淤积型透空式防波堤透射系数随周期的变化如图 11所示,上拦沙板没有下落时,水深 h 由 0.45 m 降为0.4 m,防波堤透射系数增大,将上拦沙板的位置降低 0.05 m 后,透射系数有所减小,但仍大于水深 h=0.45 m 时,当上拦沙板降低 0.10 m 后,防波堤透射系数急剧减小,消波效果明显优于水深 h=0.45 m时。由此可见,在防波堤所处位置水位变化较大时,可以通过降低上拦沙板使防波堤消波效果重新满足港池内的设计要求,当水位较高时将上拦沙板升起,使挡板开度大小充分满足水体交换所需的过流面积。图 11 透射系数随周期变化关系Fig.11 Relation of transmission coefficient with periodKoraim26研究了 3 种竖直挡板透空式防波堤的透射系数,在相对入水深度 d/h=0.25 和 d/h=0.5 时,3 种模型与本文模型透射系数与周期的曲线关系如图 12(a)和 12(b)所示。由图 12(a)可见 4种结构形式的防波堤透射系数都随周期增大而增大,但 3 种传统透空式防波堤的透射系数随周期的增长的度更快,相同工况下,防淤积型透空式防波堤的透射系数小于 3 种传统透空式防波堤。由图 12(b)可见,除了周期 T 在 1.11.3 s 时,其他工况下防淤积型透空式防波堤的透射系数小于以上 3 种透空式防波堤。还可看出,防淤积型透空式防波堤透射系数随周期增大趋于稳定,而 3 种传统透空式防波堤的透射系数随周期的增大持续增大,因此可以认为入射波浪的周期越大,采用防淤积型透空式防波堤进行消波的优势越明显。图 12 本文防波堤和传统透空式防波堤透射系数比较Fig.12Comparison of transmission coefficients between breakwater and traditional permeable breakwa-ter in this paper图 12(c)给出相对入水深度 d/h=0.4 和 d/h=0.6 时,本文防波堤与琚烈红7研究的防波堤透射0031第 8 期于通顺,等:防淤积型透空式防波堤水动力特性分析系数与周期的关系。可以看出,除个别工况外,防淤积型透空式防波堤的透射系数小于传统透空式防波堤的透射系数,并且防淤积型透空式防波堤透射系数随周期增大逐渐趋于稳定,而传统透空式防波堤的透射系数随周期增大而增大。因此可以认为在竖直挡板式防波堤基础上增加下拦沙板具有更好的消波效果,尤其在入射波浪周期较大时,这种优势更加明显。3.2 防波堤上拦沙板水平波浪力 防淤积型透空式防波堤上拦沙板水平波浪力无量纲参数 F/gH3随波陡 H/L 的变化如图 13 所示,可以看出:F/gH3随 H/L 增大而减小,但其减小速率逐渐变慢,最后趋于稳定。对比挡板位于 3 种不同位置时的波浪力曲线可以发现,在波陡较小时,上拦沙板下落距离 D 越大,水平波浪力也越大,但上拦沙板下落距离 D=0.10 m 的曲线纵坐标 F/gH3下降速率更快,当波陡增大到一定数值后,上拦沙板下落距离D=0.10 m 时的水平波浪力小于D=0.05 m时,甚至小于 D=0 m 时。这是由于上拦沙板下落后,挡板前侧面直接承受波浪荷载的面积增大,大波高波浪作用在上拦沙板前侧面的波浪力增大,挡板两侧压力差减小,故作用在上拦沙板的合力减小。图 13 水平波浪力无量纲参数 F/gH3随波陡变化关系Fig.13 Relationship between dimensionless number of horizontal wave force F/gH3 and wave steepness3.3 防波堤上拦沙板底面上托力 防淤积型透空式防波堤上拦沙板底面上托力无量纲参数 F/gH3随波陡 H/L 的变化如图 14 所示,可以看出:F/gH3随 H/L 的增大持续减小,但F/gH3减小速率逐渐变慢,趋于稳定后曲线末端略有提升。从数据点的分布来看,压力数值在持续减小,因此曲线末端的提升应与二阶函数的拟合形式有关。模拟过程中上拦沙板处于完全静止状态,从各曲线图可以得出,上拦沙板下落距离 D越大,上托力逐渐变大。随着水深增加,波浪的静水压力部分增大,而动水压力减小。挡板底面压力主要受静水压强的影响,上拦沙板底面位于水下位置越深时,静水压强的影响越大,底面上托力也越大。1031哈 尔 滨 工 程 大 学 学 报第 44 卷图 14 底面上托力无量纲参数 F/gH3随波陡变化关系Fig.14 Relationship between dimensionless number of the upper support force F/gH3 and wave steepness4 结论 1)防淤积型透空式防波堤的透射系数随着波陡 H/L 的增大而减小,具有较好的消波效果,且更适用于波浪周期较大的海域。2)工作海域水位降低时,新型透空式防波堤透射系数增大,可以通过降落上拦沙板降低透射系数,挡板入水深度增大 50%时,透射系数减小 9%28%,入水深度增大 100%,透射系数减小了 42%58%,可使防波堤的消波效果满足设计需求。3)上拦沙板水平波浪力的无量纲参数 F/gH3随波陡 H/L 的增大而减小,在波陡较小时,上拦沙板下落距离 D 越大波浪力也越大。4)上拦沙板底面上托力的无量纲参数 F/gH3随着波陡 H/L 的增大而减小,上托力主要受静水压强的影响,上拦沙板下落距离 D 越大,挡板所受上托力越大。参考文献:1 AHMED K.A comparative study between some different types of permeable breakwaters according to wave energy dissipationJ.Ain shams engineering journal,2022,13(4):101646.2 LI Yanting,WANG Dengting,SUN Tianting,et al.Fully coupled simulation of interactions among waves,permeable breakwaters and seabeds based on N-S equationsJ.Chi-na ocean engineering,2021,35(1):26-35.3 KORAIM A S.Hydrodynamic characteristics of double per-meable breakwater under regular wavesJ.Marine struc-tures,2011,24(4):503-527.4 URSELL F.The effect of a fixed vertical barrier on surface waves in deep waterJ.Mathematical proceedings of the Cambridge philosophical society,1947,43(3):374-382.5 WEIGEL R L.Closely spaced piles as a breakwaterJ.The dock and harbour authority,1961,42(491):150-165.6 邱大洪,王学庚.深水薄板式防波堤的理论分析J.水运工程,1986(4):8-12.QIU Dahong,WANG Xuegeng.Theoretical analysis of deep water thin plate breakwaterJ.Port&waterway en-gineering,1986(4):8-12.7 琚烈红,杨正己.设有挡浪板透空堤波浪透射系数实验研究J.水运工程,2008(4):19-22.JU Liehong,YANG Zhengji.Experimental research on co-efficient of wave transmission through single immersed ver-tical barrier of open-type breakwaterJ.Port&waterway engineering,2008(4):19-22.8 朱大同.直立单层挡板透空堤的水动力学特性C/中国海洋学会海洋工程分会.第十八届中国海洋(岸)工程学术讨论会论文集(下).北京,2017:894900.ZHU Datong.Hydrodynamic characteristics of vertical monolayer baffle permeable dike C/Marine Engineering 2031第 8 期于通顺,等:防淤积型透空式防波堤水动力特性分析Branch of China Oceanic Society.Proceedings of the 18th China Offshore Engineering Symposium(II).Beijing,2017:894-900.9(苏)斯米尔诺夫(.)等著,吴德镇译.港口与港口建筑物M.北京:人民交通出版社,1984.Port and port buildings M.Beijing:China Communications Press,1984.10 殷福安.挡板式透空堤透浪特性研究D.南京:河海大学,2005:18-35.YIN Fuan.Study on wave permeability characteristics of baffle-type permeable breakwater D.Nanjing:Hohai University,2005:18-35.11 NEELAMANI S.Wave interaction with partially immersed twin vertical barriersJ.Ocean engineering,2002,29(2):215-238.12 范骏,王宇楠,杨斯汉,等.双挡板透空堤透浪与反射系数实验研究J.海洋工程,2011,29(4):60-67.FAN Jun,WANG Yunan,YANG Sihan,et al.Experi-mental research on reflection coefficient and transmitted coefficient of double barriers penetrated breakwaterJ.The ocean engineering,2011,29(4):60-67.13 黄璐.T 型透空式防波堤消波性能分析D.大连:大连理工大学,2013:22-26.HUANG Lu.Analysis of wave-absorbing performance of T-shaped permeable breakwaterD.Dalian:Dalian Uni-versity of Technology,2013:22-26.14 周效国,李雨,江沭淮,等.多层直立开孔挡板透空式防波堤消浪性能试验研究J.水运工程,2014(1):31-35.ZHOU Xiaoguo,LI Yu,JIANG Shuhuai,et al.Experi-ment study on wave dissipation performance of permeable breakwater of multi-layer upright wave board with holesJ.Port&waterway engineering,2014(1):31-35.15 严以新,郑金海,曾小川,等.多层挡板桩基透空式防波堤消浪特性试验研究J.海洋工程,1998,16(1):68-75.YAN Yixin,ZHENG Jinhai,ZENG Xiaochuan,et al.Characteristics of wave dissipation for pile foundation tier retainer breakwatersJ.The ocean engineering,1998,16(1):68-75.16 RAO S,SHIRLAL K G,VARGHESE R V,et al.Physi-cal model studies on wave transmission of a submerged in-clined plate breakwaterJ.Ocean engineering,2009,36(15/16):1199-1207.17 谢怀东.黄岐渔港新型透空式防波堤设计研究J.河海大学学报(自然科学版),1997,25(4):42-46.XIE Huaidong.Study and design of new type perforated breakwater for Huangqi fishery harbour J.Journal of Hohai University(natural sciences),1997,25(4):42-46.18 左其华,窦希萍.第十三届中国海洋(岸)工程学术讨论会论文集M.北京:海洋出版社,2007:311-315.ZUO Qihua,DOU Xiping.Proceedings of the 13th China Ocean(Offshore)Engineering SymposiumM.Beijing:Ocean Press,2007:311-315.19 左其华,窦希萍.第十四届中国海洋(岸)工程学术讨论会论文集M.北京:海洋出版社,2009:704-707.ZUO Qihua,DOU Xiping.Proceedings of the 14th China Ocean(Coastal)Engineering Symposium M.Beijing:Ocean Publishing House,2009:704-707.20 于通顺.一种兼具拦沙功能的透水型装配式防波堤及其施工方法P.中国:201811199519,2018-10-15.YU Tongshun.Water-permeable assembled breakwater with sand blocking functionP.China:201811199519,2018-10-15.21 庞重光,韩丹岫,赵恩宝.黄东海悬浮泥沙浓度的垂向分布特征及其回归模型J.泥沙研究,2009(5):69-75.PANG Chongguang,HAN Danxiu,ZHAO Enbao.Char-acteristics and regression equations for vertical distribution of suspended sediment concentration in the Yellow Sea and East China SeaJ.Journal of sediment