张娜，张庆河，邹国良，蒋学炼

(1.天津城建大学 天津市软土特性与工程环境重点实验室，天津 300384； 2.天津大学 水利工程仿真与安全国家重点实验室，天津 300072； 3.南京水利科学研究院，南京 江苏 210024)

越浪堤后透射波周期特性研究

张娜1,2，张庆河2，邹国良3，蒋学炼1

(1.天津城建大学 天津市软土特性与工程环境重点实验室，天津 300384； 2.天津大学 水利工程仿真与安全国家重点实验室，天津 300072； 3.南京水利科学研究院，南京 江苏 210024)

为了研究不规则波越过光滑堤后的波周期变化，本文基于近岸波浪传播变形非静压模型分析了相对淹没深度、深水波陡、相对堤顶宽度和坡度对越浪堤波周期变化的影响。根据数值模拟的结果建立了光滑越浪堤波周期透射系数的预测公式。结果表明：堤后平均波周期受高频谐波的影响小于入射波周期；堤顶相对淹没深度以及深水波陡是波周期透射系数的主要影响因子，且与其线性相关；相对堤顶宽度的影响次之，且出水堤时对于波周期透射系数变化的影响大于淹没堤；斜坡坡度对波周期透射系数的影响最不明显。将预测公式得到的预测值与实测值相比，确定性系数高达0.95，表明其良好的预测效果。

越浪堤； 光滑堤； 非静压波浪模型； 波周期透射系数； 预测公式

越浪堤包括允许少量越浪的出水堤，也包括淹没在水中的潜堤[1-2]。修建越浪堤能对近岸工程区起到消浪掩护的作用，因此在实际工程中有着较为广泛的应用。在实际工程中，堤后波高和周期是确定堤顶高程的重要因素[2-3]。然而，已有的大多数研究主要是针对波高。除了波高，堤后透射波周期的估算在设计中也是非常必要的[4]。事实上，堤后平均波周期受高频谐波生成的影响与堤前平均波周期相比会发生变化，因此寻求一种实用有效的方法估算越浪堤后的波周期具有重要价值。

经验公式法是估算光滑堤后波参数的简单有效的方法。已有研究集中于在大量模型实验数据的基础上拟合得到波高透射系数[5-6]。Van der Meer等提出了一种波浪越过光滑结构物后的波谱预测方法[2]。这种方法假定频率低于1.5倍主频的波谱与入射波谱的波谱类型相似，且这部分低频谱能量占到了总能量的60%；而频率在1.5～3.5倍主频之间的波能量则被假定为均匀分布。Carevic等研究发现高频区的透射波能量随堤顶相对淹没深度的增加线性递减，从而进一步改进了该波谱预测模型[7]。Zhang等基于非静压数值波浪水槽模拟结果建立了光滑潜堤有效波高与平均波周期透射系数的预测公式，提高了预测精度，但该公式中没有包含允许越浪的出水堤[8]。

为此，本文将进一步利用近岸波浪传播变形非静压波浪模型对光滑越浪堤(包括淹没堤和允许越浪的出水堤)进行模拟，并在大量模拟数据分析的基础上针对平均波周期的透射系数(堤后透射平均波周期与堤前入射平均波周期的比值)建立经验预测公式[9-13]。

1 数值模拟方法

非静压波浪模型的控制方程为包含非静压项的浅水方程，其笛卡尔坐标系统的垂向二维控制方程为

(1)

(2)

(3)

式中：t表示时间，x为水平坐标，u和w分别为x方向和z方向上的流速，ζ为自由表面，ρ是密度常数，g是重力加速度，q是非静水压力，τij是水平湍流应力。

湍流应力可表示为

(4)

式中vt是由于波浪破碎引起的水平涡粘系数。

自由表面位移可通过对式(1)进行沿水深方向的积分并结合自由表面运动学条件获得，最终的自由表面方程为

(5)

式中d是静水深。

底部摩擦力采用曼宁粗糙度系数表示

(6)

式中：cf是无量纲的底部摩擦系数；n是曼宁粗糙度系数，m-1/3·s；总水深h=ζ+d。

在入流处采用域内源造波，并在入流和出流边界的位置设置了海绵层用于消波[12]。具体的边界条件、方程离散及数值计算方法详见文献[9]。

2 模型验证

非静压波浪模型已经被证实能够准确模拟不规则波与光滑越浪堤相互作用过程中的破碎、越浪、增减水等传播变形现象[10]。本文在此基础上进一步利用物理模型实验对不规则波作用下的光滑越浪堤后波参数进行验证。

2.1物理模型实验布置

验证数据来自Carevic 等[7]开展的关于光滑潜堤堤后波参数的物理模型实验。实验水槽宽1 m，高1.1 m。潜堤建在平底上，由光滑木板护面。堤顶宽度0.2 m，堤身高度0.35 m，前后坡坡度均为1∶2(图1)。采用具有主动吸能装置的活塞式造波机造波，入射不规则波谱为JONSWAP谱 (γ=3.3,σa=0.07,σb=0.09)。

图1 光滑潜堤波浪水槽实验布置图Fig.1 Experimental setup of the wave flume with a smooth submerged breakwater

2.2数学模型参数

与Carevic等[7]的物理模型实验水槽基本一致，数值波浪水槽的长度是60 m。网格步长0.02 m，垂向平均划分为2层。初始时间步长0.005 s，模拟时长20 min。波高和波谱基于300个波进行统计分析，底摩阻系数cf取为0。水槽两端各设置了海绵层用来消除波浪反射的影响，海绵层宽度至少为每组波长的3倍。

2.3数学模型验证

采用2.1节的一组实验用于数值模拟验证。入射有效波高(Hm0-i)为0.105 m，入射平均波周期(T0,2-i)为0.92 s，入射谱峰波周期(Tp-i)为1.10 s，水深是0.45 m。与图1一样，G1和G4用于计算入射和透射波参数。入射波参数是基于空槽率而定，使得G1点测得的波谱与理论入射波谱一致，最大相对误差不超过2%。透射波参数是在水槽中设置了光滑堤之后测量得到的。堤前入射波的理论谱和模拟谱的对比如图2(a)，堤后实测波谱和模拟谱的对比如图2(b)所示。结果表明，模拟谱与实测谱(理论谱)吻合较好。模拟的透射谱和实测谱的fp分别为0.91 Hz和0.93 Hz。高频部分的谱能量分布在1.5fp～3.5fp。模拟的透射有效波高(Hm0-t)、透射平均周期(T0,2-t)以及谱峰周期(Tp-t)与实测值的对比见表1。模拟值与实测值的最大相对误差不超过3.53%。其中波高误差产生的原因是物理模型实验中，无论堤表面再光滑仍然存在一定摩阻，所以物理模型的数据(实测值)与模拟值相比略偏小是正常的。误差结果表明建立的非静压数值波浪水槽能够准确地模拟不规则波越过光滑潜堤之后的波参数。

图2 数值模拟谱与入射理论谱(G1)以及堤后实测谱(G4)的对比Fig.2 Comparison of theoretical and simulated incident spectra (G1) and the measured and simulated transmitted spectra (G4)

Table1Thecomparisonbetweenthecalculatedandmeasuredtransmittedwaveparameters

波参数Hm0-t/mT0,2-t/sTp-t/s实测值0.0810.851.07模拟值0.0830.821.10相对误差/%2.47-3.532.80

3 越浪堤波浪周期变化影响因素分析

除了G1和G4用于分析入射和透射波参数，数

图3 光滑越浪堤数值水槽测点布置图Fig.3 Numerical flume layout of the smooth low-crested breakwater

3.1相对淹没深度Rc/Hm0-i的影响

图4给出了Rc/Hm0-i与Kp的关系曲线。图4(a)显示了当Rc/Hm0-i从0.5变到-1.5时，Kp<1，即越过光滑堤之后的透射平均波周期要小于入射平均波周期。这主要是由于波浪越过光滑堤之后产生高频谐波，部分能量从低频转移到高频。为进一步说明高频谐波的产生原理，采用Briganti等[14]定义的能量分布参数e.d.p.来表示f>1.5fp的高频波能在入射波谱和透射波谱能量之间的转换。当e.d.p.值大于0时，说明透射谱的高频部分能量(f>1.5fp)向低频部分发生了转移。图4(b)表明，Rc/Hm0-i从-1.5变化到0.5，对于Sop=0.01，e.d.p.值从0.39变为4.44，说明随着水深变浅，波谱能量逐渐从低频部分转移到高频部分，导致平均波周期的减小。这与Carevic等[7]的结论一致。

表2 数值模拟研究的入射波要素

注：入射波谱采用JONSWAP谱 (γ=3.3,σa=0.07,σb= 0.09)，h是总水深。

图4 Rc/Hm0-i与Kp及e.d.p的关系(B/Hm0-i=2.0，斜坡坡度为1∶2)Fig.4 Kp and e.d.p in relation to Rc/Hm0-i (B/Hm0-i is 2.0, slope is 1∶2)

图4也表明，随着Rc/Hm0-i从0.5变化到-1.5时，对于Sop分别为0.04、0.03、0.02、0.01，Kp分别增加37%、38%、41%和48%。因此Rc/Hm0-i是影响Kp的一个重要参数。Kp与Rc/Hm0-i之间存在较强的负相关，相关系数为-0.79。

3.2深水波陡Sop的影响

图5显示了Sop-Kp的关系曲线。由图5可知，Kp随着Sop的增大而增大。在波高和水深相同的情况下，Sop越大则对应的波长L越小，根据Ursell数[15](Ur=HL2/h3)，L越小则波浪的非线性越弱，堤后谐波分离越不明显，高频谐波对平均周期的影响也就越小，因此Kp也就越大。

当Sop从0.01变到0.06时，对于Rc/Hm0-i=0.5，0.0，-0.5，-1.0，-1.5时，Kp分别增加了42%、37%、34%、29%和27%。表明Sop也是影响Kp的一个重要因素。Kp与Sop之间存在较强的正相关，相关系数为0.72。

3.3相对堤顶宽度B/Hm0-i的影响

图6显示了B/Hm0-i与Kp的关系曲线。结果表明，对于Rc/Hm0-i=0.5，当B/Hm0-i从1变到8，对于Sop分别为0.01、0.02、0.03和0.04，Kp分别减小了17.0%、15.0%、13.6%和12.5%。对于Rc/Hm0-i=-0.5，当B/Hm0-i从1变到8，对于Sop分别为0.01、0.02、0.03和0.04，Kp分别减小了10.8%、9.8%、9.1%和6.0%。对比结果表明在出水堤时，B/Hm0-i对于Kp变化的影响要大于淹没堤。

图5 Sop-Kp的关系曲线 (B/Hm0-i=2.0，斜坡坡度为1∶2)Fig.5 Kp in relation to Sop (B/Hm0-i is 2.0, slope is 1∶2)

3.4坡度的影响

由图7可知，Kp随着坡度的增加而减小。当坡度从1变到3时，Kp仅减小了5%，表明相对于Rc/Hm0-i和Sop，坡度对于Kp变化的影响较小。这是由于当波浪沿着缓坡越过光滑堤的过程中，相比于陡坡有更多的波能从低频转移到高频。图7中的T0,2空间演化可较好的解释Kp随坡度的增加而减小的现象。

图6 B/Hm0-i-Kp的关系曲线 (斜坡坡度为1∶2)Fig.6 Kp in relation to B/Hm0-i (slope is 1∶2)

图7 不同斜坡坡度的T0,2空间演化(Rc/Hm0-i=-0.5，B/Hm0-i=2.0，Sop=0.04)Fig.7 The spatial evolution of T0,2 (b) (Rc/Hm0-i is -0.5, B/Hm0-i is 2.0 and Sop is 0.04)

4 周期变化经验公式的建立

根据第3节的影响因素分析可知，Rc/Hm0-i和Sop是Kp的主要影响因子。B/Hm0-i的影响次之，且出水堤时对于Kp的影响大于淹没堤。前坡坡度对Kp的影响最不明显。为此，基于Rc/Hm0-i、Sop和B/Hm0-i，采用最小二乘法拟合建立光滑越浪堤

Kp的透射系数经验公式[16]。

根据上一节的分析可知，Kp与Rc/Hm0-i和Sop都具有较好的线性相关关系，因此Kp的待拟合公式可表示为

Kp=a1Rc/Hm0-i+a2Sop+a3B/Hm0-i+a4

(7)

图8 Kp的拟合值与模拟值对比以及文献[7]和文献[2]的实测值与Kp预测值的对比Fig.8 The comparison of Kp based on simulations and calculated values and measurements from Ref.[7] and Ref.[2] with calculated values

5 结论

1)根据数值模拟的结果分析，波浪越过光滑越浪堤后的平均波周期会小于入射值，这主要是由于波浪越堤过程中产生的高频谐波引起的。Rc/Hm0-i和Sop是Kp的主要影响因子，且与其线性相关。B/Hm0-i的影响次之，且出水堤时对于Kp的影响大于淹没堤。前坡坡度对Kp的影响最不明显。

该成果可为越浪堤工程设计提供依据，并为进一步研究带有护面的越浪堤波参数演化奠定基础。

[1] 冯卫兵, 陈凯华, 薛波, 等. 不同护面块体潜堤的消浪特性试验研究[J]. 水运工程, 2014(3): 35-40.

FENG Weibing, CHEN Kaihua, XUE Bo, et al. Lab research on wave dissipation characteristics of different armor blocks used in submerged breakwater[J]. Port & waterway engineering, 2014(3): 35-40.

[2] VAN DER MEER J W, REGELING H J, DE WAAL J P. Wave transmission: spectral changes and its effects on run-up and overtopping[C]// Coastal Engineering Conference， ASCE, 2000: 2156-2168.

[3] VAN DER MEER J W, BRIGANTI R, ZANUTTIGH B, et al. Wave transmission and reflection at low-crested structures: Design formulae, oblique wave attack and spectral change[J]. Coastal engineering, 2005, 52(10): 915-929.

[4] 梁书秀, 常艳玲, 孙昭晨. 深水破碎波周期变化的试验研究[J]. 水运工程, 2016(6): 12-17.

LIANG Shuxiu, CHANG Yanling, SUN Zhaochen. Experimental study on evolution characteristics of wave period of deep-water breaking waves[J]. Port & waterway engineering, 2016(6): 12-17.

[5] 冯卫兵, 王明明, 崔川川. 潜堤透浪系数的试验研究[J]. 水运工程, 2012(9): 1-6.

FENG Weibing, WANG Mingming, CUI Chuanchuan. Experimental study of wave transmission coefficient over a submerged breakwater[J]. Port & waterway engineering, 2012(9): 1-6.

[6] VAN DER MEER J W, WANG B, WOLTERS A, et al. Oblique wave transmission over low-crested structures[C]// Proceedings Coastal Structures, USA, 2003: 567-579.

[7] CAREVIC D, LONCAR G, PRSIC M. Wave parameters after smooth submerged breakwater[J]. Coastal engineering, 2013, 79: 32-41.

[8] ZHANG N, ZHANG Q, ZOU G, et al. Estimation of the transmission coefficients of wave height and period after smooth submerged breakwater using a non-hydrostatic wave model[J]. Ocean engineering, 2016, 122: 202-214.

[9] ZIJLEMA M, STELLING G, SMIT P. SWASH: An operational public domain code for simulating wave fields and rapidly varied flows in coastal waters[J]. Coastal engineering, 2011, 58(10): 992-1012.

[10] 张娜, 邹国良. 斜坡上波浪破碎与越浪非静压数值模拟[J]. 海洋工程, 2015, 33(2): 32-41.

ZHANG Na, ZOU Guoliang. Numerical simulation of wave transformation and overtopping based on non-hydrostatic model[J]. The ocean engineering, 2015, 33(2): 32-41.

[11] 邹国良, 张庆河. 非静压方程与波作用谱模型的波浪传播嵌套模拟[J]. 哈尔滨工程大学学报, 2014, 35(1): 126-131.

ZOU Guoliang, ZHANG Qinghe. Simulation of wave transformation by nesting the non-hydrostatic equation and wave action spectrum model[J]. Journal of Harbin Engineering University, 2014, 35(1): 126-131.

[12] 邹国良, 张庆河. 非静压波浪模型无反射造波[J]. 海洋工程, 2012, 30(4): 55-61.

ZOU Guoliang, ZHANG Qinghe. Wave generation without re-reflection for non-hydrostatic wave model[J]. The ocean engineering, 2012, 30(4): 55-61.

[13] 王岗, 郑金海. 非静压波浪模型新理论[C]// 第十七届中国海洋(岸)工程学术讨论会论文集(下).南京，2015: 763-772.

WANG Gang, ZHENG Jinhai. New theory of non-hydrostatic wave model[C]// Proceedings of the 17th China Marine (Offshore) Engineering Symposium. Nanning, 2015: 763-772.

[14] BRIGANTI R, VAN DER MEER J W, BUCCINO M, et al. Wave transmission behind low crested structures[C]// Proceedings of the 3rd Coastal Structures Conference，ASCE, 2003: 580-592.

[15] URSELL F. Short surface waves due to an oscillating immersed body[C]// Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences，The Royal Society, 1953, 220(1140): 90-103.

[16] GHASEMI S E, HATAMI M, AHANGAR G H R M, et al. Electrohydrodynamic flow analysis in a circular cylindrical conduit using least square method[J]. Journal of electrostatics, 2014, 72(1): 47-52.

[17] KURZ-KIM J R, LORETAN M. On the properties of the coefficient of determination in regression models with infinite variance variables[J]. Journal of econometrics, 2014, 181(1): 15-24.

本文引用格式：

张娜，张庆河，邹国良，等. 越浪堤后透射波周期特性研究[J]. 哈尔滨工程大学学报， 2017， 38(11): 1715 -1720.

ZHANG Na, ZHANG Qinghe, ZOU Guoliang, et al. Study on transmitted mean wave period after low-crested breakwater[J]. Journal of Harbin Engineering University， 2017， 38(11): 1715 -1720.

Studyontransmittedmeanwaveperiodafterlow-crestedbreakwater

ZHANG Na1,2, ZHANG Qinghe2, ZOU Guoliang3, JIANG Xuelian1

(1.Tianjin Key Laboratory of Soft Soil Characteristics & Engineering Environment, Tianjin Chengjian University, Tianjin 300384, China; 2.State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China; 3.Nanjing Hydraulic Research Institute, Nanjing 210024, China)

In this paper, to study the change of mean wave period after the smooth low-crested breakwater under irregular waves, the influence of the relative submergences, wave steepness, the relative crest width and slope on wave period variation of low-crested breakwater was discussed based on the non-hydrostatic model simulating waves till shore. Based on the numerical results, an empirical formula of transmission coefficient of mean wave period was proposed. According to the analysis of the numerical results, the transmitted mean wave period is less than the incident mean wave period because of the generation of high-frequency harmonics. Moreover, the relative submergences and wave steepness are the main factors linearly influencing the transmission coefficient of mean wave period. The relative crest width also has considerable effect on the transmission coefficient of mean wave period especially when the structure is non-immersed. However, the seaside slope of the structure shows little effect on the transmission coefficient of mean wave period. The coefficients of determination between measurements and values calculated using the empirical formula is 0.95, representing good predicting performance for mean wave period over a smooth low-crested breakwater.

low-crested breakwater; smooth breakwater; non-hydrostatic wave model; transmission coefficient of mean wave period; empirical formula

10.11990/jheu.201608058

http://www.cnki.net/kcms/detail/23.1390.u.20170427.1321.012.html

TV139.2

A

1006-7043(2017)11-1715-06

2016-08-29.

网络出版日期：2017-04-27.

国家自然科学基金项目(51509177，51609152)；天津市自然科学基金项目(14JCYBJC22100)；天津市建交委科技计划(2013-8).

张娜(1978-), 女, 副教授；

张庆河(1966-), 男, 教授，博士生导师.

张庆河，E-mail: qhzhang@tju.edu.cn.