珊瑚礁礁坪宽度对波浪传播变形及增水影响的实验研究

姚宇, 何天城, 唐政江, 蒋昌波, 3



珊瑚礁礁坪宽度对波浪传播变形及增水影响的实验研究

姚宇1, 2, 何天城2, 唐政江2, 蒋昌波2, 3

1. 河海大学, 海岸灾害及防护教育部重点实验室, 江苏 南京 210098; 2. 长沙理工大学水利工程学院, 湖南 长沙 410004; 3. 水沙科学与水灾害防治湖南省重点实验室, 湖南 长沙 410004

通过在波浪水槽中进行一系列物理模型实验, 研究珊瑚礁礁坪宽度变化对珊瑚礁海岸附近波浪传播变形及礁坪上波浪增水的影响。物理实验采用理想化的珊瑚礁模型, 测试了3种礁坪宽度下的一系列不规则波工况。实验结果分析表明: 波浪沿礁传播过程中, 短波持续衰减; 低频长波波高沿礁逐渐增大, 直到海岸线附近达到最大; 随着礁坪宽度的增加, 海岸线附近的短波波高呈下降趋势, 低频长波波高的变化规律不显著; 礁坪上的波浪增水受礁坪宽度变化的影响不明显; 通过对海岸线附近的波浪进行频谱分析发现, 礁坪上低频长波的运动存在着一阶共振模式, 且共振放大效应强度受礁坪水深、入射波峰周期和礁坪宽度共同影响。

低频长波; 波浪增水; 礁坪宽度; 礁坪共振

近年来, 我国在南中国海的珊瑚岛礁周围进行了填礁造岛工程, 并修建了诸如灯塔、码头、机场、油气开采平台等设施。此类工程设施常处于远海地区, 海况恶劣, 周围海底地形变化剧烈(姚宇等, 2017)。如何保障已建成的岛礁建筑物在风暴潮等极端气候下的安全, 降低因洪涝灾害造成的人员和财产损失引起了学者们的广泛关注。典型的珊瑚礁地形主要由礁前斜坡与礁坪组成, 礁前斜坡与礁坪的连接处称为礁缘。在常浪条件下, 波浪作用于礁前斜坡, 通过礁缘处的破碎和礁坪上的摩擦损耗后, 到达礁后岸滩用于波浪爬高的能量几乎可以忽略不计。波浪爬高是珊瑚礁海岸地区洪涝灾害发生的主要因素, 主要由短波(频率区段为0.04~0.4Hz)、低频长波(频率区段为0.001~0.04Hz)和波浪增水共同组成(Merrifield et al, 2014)。其中低频长波成分在某些情况下特别是由风暴潮引起的高潮位时, 易于在半封闭礁坪上发生共振现象引起海岸线附近波浪的放大效应, 加剧波浪海岸爬高即发生洪涝灾害的风险(Péquignet et al, 2009)。

物理模型实验是国内外研究波浪沿礁演变规律的手段之一。例如Gourlay (1996)进行了一系列的模型实验, 较为全面地研究了规则波在概化珊瑚礁物理模型上的传播, 增水和波生流问题。梅弢等(2013)通过实验模拟了常年平均波浪和重现期为50a两种波浪条件, 研究了波浪在珊瑚礁上的传播规律。丁军等(2015)通过港池实验模拟某典型岛礁附近的真实地形, 研究了岛礁地形上波浪的传播变形规律。姚宇等(2017)通过开展波流水槽实验, 研究了潮汐流影响下规则波在珊瑚礁海岸地形上的传播变形和增水规律。上述物理模型实验研究均未考虑礁形变化等因素对波浪在珊瑚礁地形上传播变形的影响。

近年来, 此领域的学者们主要通过数值模拟对礁形变化影响下波浪的传播变形规律进行研究, 例如Shimozono等(2015)以超强台风“海燕”过境为背景, 基于Boussinesq方程对太平洋Eastern Samar岛珊瑚礁海岸的波浪爬高进行了研究, 结果显示低频长波和短波是岸滩爬高的主要成分且二者所占的比重受到礁坪宽度和岸滩坡度的影响。Quataert等(2015)在马绍尔群岛Roi-Namur礁采集现场观测数据, 基于XBeach模型研究了礁前斜坡坡度、礁坪宽度和礁面糙率等礁形变化对海岸洪涝灾害的影响。

目前文献中尚未发现有开展礁坪宽度变化对波浪传播变形影响的物理模型实验工作。本文通过波浪水槽物理模型实验, 研究珊瑚礁礁坪宽度变化对波浪传播变形和礁坪上波浪增水的影响, 并通过波谱分析了海岸线附近低频长波能量的变化规律。研究结果将为评估风暴潮极端波浪条件下珊瑚礁对岛礁工程建设的保护作用提供一定的理论参考。

1 实验设置

图1为物理模型实验设置示意图。实验在长沙理工大学水利实验中心长40m, 宽0.5m, 高0.8m的波浪水槽中进行。水槽左端为一台推板式造波机, 在距造波机27.3m处采用坡度为1 : 6的斜面模拟礁前斜坡, 斜坡后接一定长度的水平平台模拟礁坪, 礁坪距水槽底高度为0.35m, 礁坪后接坡度为1 : 3.3的斜面模拟礁后岸滩。整个礁体物理模型由PVC材料制作, 并通过固定支架悬挂于水槽壁上, PVC板宽度与水槽宽度一致。使用玻璃胶填充模型与水槽壁之间、礁坪相接处及礁前斜坡与水槽底部之间的缝隙。

图1 实验布置图

图中0为远海水深;r为礁坪水深

Fig. 1 Experimental setup.0donates the offshore water depth;rdonatesthe reef-flat water depth

实验采用8个电容式浪高仪(G1~G8)测量自由液面高程的变化, 布置方式如图1所示。在外海侧设置浪高仪G1、G2和G3用以测量分离入射波和反射波, G4被放置于礁坪斜坡中部测量波浪的浅水变形。G5、G7和G8分别布置于礁缘附近(礁缘向岸侧0.2m)、礁坪中部和海岸线附近(岸滩斜坡坡脚向海侧0.2m), 用来监测从礁缘向海岸线波浪的传播变形过程。G6放置于G5和G7浪高仪的正中间位置。在实验过程中, 所有浪高仪通过数据采集系统同步采集, 采样频率为50Hz, 采样时长为自造波机启动后的1200s。

实验中的物理模型尺寸参照Becker等(2014)对马绍尔群岛进行现场观测的资料, 按照弗如德相似准则确立, 几何比尺为1 : 20。实验测试了4个深水有效波高S0(0.04m、0.06m、0.08m、0.10m)、2个礁坪水深r(0.05m、0.10m)和3个峰周期p(1.0s、1.5s、2.0s)组合下的一系列不规则波工况。在上述所有波况下, 分别测试3个礁坪宽度r(4.8m、8m、11.2m)对波浪传播变形的影响, 对应的原型礁坪宽度尺寸分别为96m、160m、224m。为了消除初始数据不稳定性的影响, 本文选取实验测量的200~ 1200s时间段的波浪数据进行分析。不规则波的波浪频谱根据浪高仪所测得的自由液面时间序列通过快速傅立叶变化(FFT)的方法分析得到。对于短波波高与低频长波波高的计算, 由于物理实验与现场观测有所区别, 因此本文采用Buckley等(2016)对低频长波和短波频段的划分方法, 将入射波的谱峰频率的二分之一作为低频长波和短波频段划分的界限。短波波高和低频长波波高的具体计算公式如下:

式中:SS为短波波高(单位: m);IG为低频长波波高(单位: m);p是入射波的谱峰频率(单位: Hz);()为波浪谱(单位: m2×s)。

2 结果分析

2.1 礁坪宽度对沿礁波高变化的影响

2.1.1 波高的沿礁变化

图2展示了典型波浪工况S0=0.08m,p=1.5s,r=0.05m时的短波波高(SS)和低频长波波高(IG)在不同礁坪宽度(r)情况下的沿礁变化规律。图2表明: 波浪以群波的包络形式由远海向近岸传播时, 远海处(G1~G3)存在少量的低频长波, 这可能是由于在包络线中束缚着一定的低频长波(van Dongerenet al, 2013)。到达礁前斜坡附近(G4)时, 由于浅水作用导致SS值增大, 随后波浪在礁缘附近(G5)发生破碎,SS值迅速减小, 而IG值则增加显著, 这是由于群波效应造成的破碎点的移动产生的(Pomeroy et al, 2012)。波浪到达礁坪中部(G7)时, 破碎过程结束, 继续向海岸线附近(G8)传播时, 由于礁面摩擦导致SS值沿礁持续减小并在G8时达到最小值, 而IG值沿礁持续增大并在G8达到最大, 这是由于低频长波在礁坪上存在一阶共振放大效应(详见2.4节)。对于总波高T, 波浪在破碎前(G1~G4), 其成分主要是SS; 随着波浪从G5传播到G7,SS成分逐渐减小,IG成分逐渐增大; 当波浪传播至G8, 其成分主要为IG。

图2 不同礁坪宽度()下短波波高(HSS)、低频长波波高(HIG)和总波高(HT)的沿礁变化(波况: HS0=0.08m, Tp=1.5s, hr=0.05m)

2.1.2 海岸线附近波高的变化

为研究礁坪宽度的变化对海岸线附近(G8)波高变化的影响, 图3展示了S0=0.08m时G8的短波波高(SS)和低频长波波高(IG)随礁坪宽度(r)的变化规律。图3a表明:SS随礁坪水深的减小而减小, 这是由于水深较小时将会增大波浪破碎强度;SS随入射波峰周期的增大而增大; 同时SS随着礁坪宽度的增大而减小, 这是由于礁坪宽度的增大导致礁面摩阻损耗增加。图3b表明:IG随礁坪水深的增大而减小, 这可能是由于水深的增大限制波浪在礁缘附近的破碎强度, 从而减少了上述由于破碎点移动产生的IG; 同时入射波峰周期越大, G8附近的IG值越大, 这是由于峰周期的增大减小了频谱的带宽, 增加了包络线中低频长波的能量(Yao et al, 2016);IG随着礁坪宽度的增长变化规律不明确, 这与礁坪宽度的变化改变了珊瑚礁地形固有共振频率, 因此在不同的波浪条件下会产生不同强度的共振响应有关。

图3 海岸线附近(G8)短波波高(HSS)和低频长波波高(HIG)随礁坪宽度(Wr)的变化(HS0=0.08m)

图中圆形表示p=1.0s; 方形表示p=1.5s; 三角形表示p=2.0s; 实心表示r=0.05m; 空心表示r=0.1m

Fig. 3 Variations of swell height (SS) and infragravity wave height (IG) with reef-flat width (r) near the coastline at G8 (S0=0.08m). Circles denotep=1.5s; Squares denotep=1.5s; Triangles denotep=2.0s; Solid markers denoter=0.05m; Open markers denoter=0.1m

2.2 礁坪波浪增水的变化规律

图4 礁坪波浪增水()随入射有效波高(a)和礁坪宽度(b)的变化

图中圆形表示p=1.0s; 方形表示p=1.5s; 三角形表示p=2.0s; 实心表示r=0.05m; 空心表示r=0.1m

2.3 波浪频谱分析

为了揭示波浪能量的沿礁变化, 图5展示了典型入射波况S0=0.08m,p=1.5s,r=0.05m时, 3种礁坪宽度(r)影响下4个代表性测点G1 (远海处)、G5 (礁缘附近)、G7 (礁坪中部)和G8 (海岸线附近)的波浪频谱图。观察图5波浪频谱图, 可以发现: G1附近存在少量的低频长波能量(= 0~ 0.33Hz), 这是波群向岸传播时, 存在少量束缚在包络线中的低频长波(van Dongerenet al, 2013)。在G5附近, 波浪发生破碎, 短波区间(> 0.33Hz)能量衰减剧烈, 同时与图2对应, G5处低频长波区间能量有较大的增长, 这同样是由于群波效应引起的破碎点移动产生的低频长波(Pomeroy et al, 2012)。波浪在礁坪上继续向岸传播时, 短波能量由于礁坪摩擦继续减少, 低频长波能量由于共振放大效应沿礁增大。波浪传播到达G8时, 短波能量已大幅削弱, 而低频长波能量则达到最大值。

对比3种不同礁坪宽度情况, 发现海岸线G8附近波浪在短波频段的能量随着礁坪宽度的增大而减小; 而在礁坪宽度r=8m时低频长波频段的能量明显高于r为4.8m和11.2m时的能量, 具体原因分析见2.4节。

图5 不同礁坪宽度(Wr)条件下波浪谱的沿礁变化(HS0=0.08m, Tp=1.5s, hr=0.05m)

Fig. 5 Variation of wave spectra with different reef-flat widths (r) across the reef profile (S0=0.08m,p=1.5s,r=0.05m).The dashed line denotes the splitting frequency between the swell wave and the infragravity wave frequency bands

2.4 有关礁坪共振的讨论

礁缘到海岸线附近低频长波能量的放大可以利用礁坪共振理论进行分析。珊瑚岛礁地形近似于半开放盆地港湾, 理论上礁坪共振周期可参照开放盆地的共振模式(Nwogu et al, 2010)进行分析, 其形式如下:

为了进一步研究礁坪水深(r)、波浪峰周期(p)和礁坪宽度(r)对海岸线附近波浪低频长波频段频谱以及可能存在的共振的影响, 图6展示了2种礁坪水深、3种波浪峰周期和3种礁坪宽度情况下G8位置(海岸线附近)低频长波频段的波浪频谱图。

观察图6可发现波浪频谱中的峰值频率基本集中在一阶礁坪共振频率附近, 可见珊瑚礁礁坪上低频长波的运动存在着一阶共振模式。对比图6a和6b可知:r=0.05m时较r=0.1m到达海岸线的低频长波具有更大的能量, 这与Quataert等(2015)现场观测发现的礁坪水深越大, 海岸低频长波能量越大的结论并不一致, 这是因为此工况(S0=0.08m,p=1.5s,r=8m)在r=0.05m时更接近于一阶理论共振频率从而在海岸线激发了更大的低频长波共振放大效应。

对比图6a、6c和6d可知: G8附近低频长波的能量随着波浪峰周期的增大而增大, 这是由于峰周期的增大减小了频谱的带宽, 增加了包络线中低频长波的能量所造成(Yao et al, 2016)。对比图6a、6e和6f可知:r=8m时低频长波能量明显大于r为11.2m和4.8m时的能量值, 且可以发现r=8m时低频长波频段的峰值频率最接近一阶礁坪理论共振频率, 表明在入射波条件相同时,r=8m相较于r为4.8m和11.2m, 更有利于激发珊瑚礁坪的一阶共振模式。这与图3b中对应工况(S0=0.08m,p=1.5s,r=0.05m)在不同礁坪宽度情况下, 低频长波波高出现的先增大后减小规律相对应。通过上述分析, 可知共振放大效应强度由礁坪水深、入射波峰周期和礁坪宽度共同决定。

图6 海岸线附近(G8)低频波浪频谱(f = 0~0.33Hz)随礁坪水深(hr), 波浪周期(Tp)和礁坪宽度(Wr)的变化

a.S0=0.08m,p=1.5s,r=0.05m,r=8m; b.S0=0.08m,p=1.5s,r=0.1m,r=8m; c.S0=0.08m,p=1.0s,r=0.05m,r=8m; d.S0=0.08m,p=2.0s,r=0.05m,r=8m; e.S0=0.08m,p=1.5s,r=0.05m,r=4.8m; f.S0=0.08m,p=1.5s,r=0.05m,r=11.2m。各子图中左虚线为一阶共振频率, 右边虚线为二阶共振频率

3 结论

本文通过物理模型实验研究礁坪宽度变化影响下珊瑚礁海岸附近波浪传播变形。主要结论如下: 1) 不规则波通过礁坪向岸传播过程中, 短波波高由于波浪破碎以及礁面摩阻损耗逐渐衰减, 低频长波则由波浪破碎点的移动产生并在礁坪上发生共振效应逐渐增大; 2) 礁坪宽度的增大能显著降低波浪在岸线附近的短波波高, 而低频长波随礁坪宽度变化的规律不显著; 3) 礁坪宽度变化时礁坪波浪增水值趋于一致; 4) 波谱分析进一步表明, 珊瑚礁坪上低频长波的运动存在着一阶共振模式, 且共振放大效应强度受礁坪水深、入射波峰周期和礁坪宽度共同影响。

丁军, 田超, 王志东, 等, 2015. 近岛礁波浪传播变形模型试验研究[J]. 水动力学研究与进展, 30(2): 194–200. DING JUN, TIAN CHAO, WANG ZHIDONG, et al, 2015. Experimental research on wave deformation near the typical island[J]. Chinese Journal of Hydrodynamics, 30(2): 194–200 (in Chinese with English abstract).

梅弢, 高峰, 2013. 波浪在珊瑚礁坪上传播的水槽试验研究[J]. 水道港口, 34(1): 13–18. MEI TAO, GAO FENG, 2013. Flume experiment research on law of wave propagation in reef flat[J]. Journal of Waterway and Harbor, 34(1): 13–18 (in Chinese with English abstract).

姚宇, 唐政江, 杜睿超, 等, 2017. 潮汐流影响下珊瑚岛礁附近波浪传播变形和增水试验[J]. 水科学进展, 28(4): 614–621. YAO YU, TANG ZHENGJIANG, DU RUICHAO, et al, 2017. Laboratory study of wave transformation and wave-induced setup over reef islands under the effect of tidal current[J]. Advances in Water Science, 28(4): 614–621 (in Chinese with English abstract).

BECKER J M, MERRIFIELD M A, FORD M, 2014. Water level effects on breaking wave setup for Pacific Island fringing reefs[J]. Journal of Geophysical Research: Oceans, 119(2): 914–932.

BUCKLEY M L, LOWE R J, HANSEN J E, et al, 2016. Wave setup over a fringing reef with large bottom roughness[J]. Journal of Physical Oceanography, 46(8): 2317–2333.

GOURLAY M R, 1996. Wave set-up on coral reefs. 1. Set-up and wave-generated flow on an idealised two dimensional horizontal reef[J]. Coastal Engineering, 27(3–4): 161–193.

MERRIFIELD M A, BECKER J M, FORD M, et al, 2014. Observations and estimates of wave-driven water level extremes at the Marshall Islands[J]. Geophysical Research Letters, 41(20): 7245–7253.

NWOGU O, DEMIRBILEK Z, 2010. Infragravity wave motions and runup over shallow fringing reefs[J]. Journal of Waterway, Port, Coast, and Ocean Engineering, 136(6): 295–305.

PÉQUIGNET A C N, BECKER J M, MERRIFIELD M A, et al, 2009. Forcing of resonant modes on a fringing reef during tropical storm Man-Yi[J]. Geophysical Research Letters, 36(3): L03607.

POMEROY A, LOWE R, SYMONDS G, et al, 2012. The dynamics of infragravity wave transformation over a fringing reef[J]. Journal of Geophysical Research: Oceans, 117(C11): C11022.

QUATAERT E, STORLAZZI C, VAN ROOIJEN A, et al, 2015. The influence of coral reefs and climate change on wave driven flooding of tropical coastlines[J]. Geophysical Research Letters,42(15): 6407–6415.

SHIMOZONO T, TAJIMA Y, KENNEDY A B, et al, 2015. Combined infragravity wave and sea-swell runup over fringing reefs by super typhoon Haiyan[J]. Journal of Geophysical Research: Oceans, 120(6): 4463–4486.

VAN DONGEREN A, LOWE R, POMEROY A, et al, 2013. Numerical modeling of low-frequency wave dynamics over a fringing coral reef[J]. Coastal Engineering, 73: 178–190.

YAO YU, BECKER J M, FORD M R, et al, 2016. Modeling wave processes over fringing reefs with an excavation pit[J]. Coastal Engineering, 109: 9–19.

Laboratory study on the effect of varying reef-flat width on wave transformation and wave-induced setup

YAO Yu1, 2, HE Tiancheng2, TANG Zhengjiang2, JIANG Changbo2, 3

1. Key Laboratory of Coastal Disasters and Defense of Ministry of Education,Hohai University, Nanjing 210098, China; 2. School of Hydraulic Engineering, Changsha University of Science & Technology, Changsha 410004, China; 3. Key Laboratory of Water-Sediment Sciences and Water Disaster Prevention of Hunan Province, Changsha 410004, China

A series of laboratory experiments were carried out in a wave flume to study the influences of different reef-flat widths on wave transformation over reef coasts as well as wave-induced setup on reef flats. A series of incident irregular wave conditions were tested based on the idealized reef model with three reef-flat widths. Results show that as waves propagate across the reef, the swell height decays all the way, whereas the infragravity wave height increases gradually along the reef and reaches the maximum at the shoreline. Near the shoreline, the swell height decreases with increasing reef-flat width, but the variation of infragravity wave height with the reef-flat width is not obvious. The effect of varying reef-flat width on wave-induced setup over the reef flat is insignificant. The first resonant mode is identified on reef flat for the infragravity waves, and the level of resonant amplification near the shoreline is determined by the combined effect of reef-flat water depth, incident peak wave period and reef-flat width.

infragravity wave; wave-induced setup; reef-flat width; reef-flat resonance

2018-07-29;

2018-11-05. Editor: YIN Bo

National Natural Science Foundation of China (51679014);Hunan Science and Technology Plan Program (2017RS3035); Open Foundation of Key Laboratory of Coastal Disasters and Defense of Ministry of Education (201602)

YAO Yu. E-mail: yaoyu821101@163.com

P731.22; P731.3; P754.2; P737.11; P714.13

A

1009-5470(2019)02-0013-07

2018-07-29;

2018-11-05。殷波编辑

国家自然科学基金(51679014); 湖南省科技计划(2017RS3035); 河海大学海岸灾害及防护教育部重点实验室开放研究基金(201602)

姚宇(1982—), 男, 湖南省湘潭市人, 副教授, 从事近海水动力学研究。E-mail: yaoyu821101@l163.com

10.11978/2018079