本研究的主要作者是来自法国波尔多大学UMR 5805 EPOC实验室的Kévin Martins和Philippe Bonneton，以及来自格勒诺布尔阿尔卑斯大学UMR 5519 LEGI实验室的Hervé Michallet。这项题为“缓坡实验室海滩上传播和破碎的非线性波浪的弥散特性”的研究发表于Elsevier旗下的期刊《Coastal Engineering》2021年第167卷。
本研究属于海岸工程与近岸水动力学领域，核心关注点是风生短波（短波）在近岸缓坡海滩上传播、演变直至破碎过程中的非线性特性及其对波浪弥散关系的影响。线性波浪理论在描述波浪传播的许多物理过程（如折射、浅水变形）方面表现良好，是一种广泛应用的一阶近似。然而，当波浪传播至近岸浅水区，随着水深变浅、地形坡度影响以及非线性效应增强，波浪形态会发生显著变化，从近对称变为水平方向上的偏斜（skewed），并在破碎前及破波带内呈现出垂直方向上的不对称（asymmetric）。这些非线性变化不仅影响海面高程，也影响近底波动速度，从而对海滩的短期和长期地貌演变起到关键作用。
这种变化主要由弱非线性三波相互作用（triad interactions）驱动，即两个主频率分量（f1, f2）通过和或差相互作用，激发一个次生频率分量f’。这个次生分量与主分量“捆绑”（bound）在一起，不遵循经典的线性波浪弥散关系。与深水区需要长距离传播才能生效的四波共振相互作用不同，近岸区的三波相互作用是非共振或近共振的，能量转换效率非常高。这导致波浪场中同时存在自由波分量和强迫波分量，对于一个给定的频率，不再存在唯一的波数与之对应。当强迫波分量在频谱的某些区域占主导时，会导致从线性弥散关系预测得到的单值优势波数和相速度谱出现显著偏差。
尽管过去已有基于现场观测的研究证实了这些非线性弥散特性的偏离，但现场研究通常存在一些局限：如通常使用水下压力或流速计进行测量，需要基于线性弥散关系进行水深衰减校正，这本身就可能在高频处引入误差；空间分辨率不足；以及观测距离有限等。因此，高频处非线性效应在何处开始占主导，即波浪场空间结构随距离变化的详细特征，仍然不甚明确。
本研究旨在利用一个高分辨率的实验室数据集“GLBEX”，来深入探究单向不规则波浪在缓坡海滩上传播和破碎过程中的弥散特性。该数据集直接测量了海面高程，空间分辨率极高（在传播的任何阶段均有数个测点对应一个波长），且波浪为单向传播，排除了方向性的影响。这使得本研究能够以前所未有的精度，量化非线性能量传递对波浪弥散特性的影响，并明确识别波浪传播过程中的不同阶段。
1. 实验数据集 研究采用了“缓坡海滩实验”（GLBEX）项目在荷兰代尔夫特Scheldegoot波浪水槽中采集的数据。水槽长110米，带有1：80的混凝土缓坡。实验使用了21个电容式和电阻式波高仪（采样频率128 Hz），以高空间分辨率沿水槽布置，直接测量自由水面高程。通过移动测点并重复实验，实现了对整个水槽高密度的空间覆盖。研究选取了A系列（A1, A2, A3）三个70分钟的不规则波试验进行分析，这些试验涵盖了从温和到剧烈、从宽谱到窄谱的不同波浪条件（详见附表1）。
2. 分析方法 本研究主要采用了两种频谱分析方法：交叉谱分析和双谱分析。 * 交叉谱分析：用于计算优势波数谱κ(f)和相速度谱c(f)。该方法通过分析相邻两个波高仪信号的相位差，得到每个频率对应的波数（κ(f)=相位差/测点间距）和相速度（c(f)=2πf*测点间距/相位差）。这里的κ(f)代表在能量平均意义上的“单值”优势波数，它反映了混合了自由波和强迫波的波浪场的总体传播特性。为确保结果的可靠性，研究采用了多个不同间距的测点组合进行估计，并对相干性高于一定阈值（coherence^2 > 0.25）的结果进行集合平均。 * 双谱分析：用于量化非线性三波相互作用的强度和非线性能量传递。双谱B(f1, f2)表征了频率f1、f2和f1+f2三个分量之间的相位耦合强度。双相位（biphase）包含了相对相位的信息，决定了近共振分量间能量传递的量。双谱分析可以帮助理解高频强迫能量的来源和强度。
3. 非线性能量传递的量化模型 为了量化高频处强迫能量的相对重要性，研究采用了基于Boussinesq近似的谱能量平衡模型。该模型描述了波浪能量通量谱F(f)在跨岸方向的变化，由非线性能量交换源项Snl(f)和耗散项Sdis(f)决定。本研究在忽略破波耗散（仅考虑底部摩擦）的前提下，直接利用实验测量得到的双谱来计算源项Snl，然后通过空间积分估算从起始点到目标点之间通过非线性耦合获得（或损失）的能量通量。该积分值与目标点的总能量通量之比，被用作该点处强迫能量相对比例的近似估计。这种方法避免了从起始点向前预测双谱的不确定性。
通过对A2和A3试验数据的交叉谱分析，研究人员根据局部厄塞尔数（Ursell number, Ur，定义为非线性参数ε与弥散参数μ之比）的变化，清晰地识别出波浪传播的四个不同阶段：
阶段I：线性传播阶段（深水区，Ur ~ 0.01-0.1） 主要结果：在深度较大的区域，波浪场主要由自由波分量主导。测量的优势波数κ(f)和相速度c(f)与线性弥散关系的预测高度吻合，偏差极小。 逻辑关系与贡献：此阶段确立了基线，表明在非线性效应可忽略的区域，线性理论是有效的。
阶段II：弱非线性浅水变形阶段（远离平均破碎点，Ur ~ 0.3） 主要结果：这是本研究揭示的关键阶段。 * 对于主频分量（约0.6fp至1.5fp），其波数和相速度仍遵循线性弥散关系，空间结构可由线性理论很好描述。 * 但对于高频分量（>1.5fp），波数κ(f)和相速度c(f)谱开始出现显著偏离。在窄谱试验A3中，2fp和3fp处的相速度与峰值频率fp处的相速度几乎相等；而在宽谱试验A2中，这些高频的相速度值则介于fp处的相速度和线性预测值之间。 * 双谱分析表明，该阶段存在明显的非线性能量向高频（特别是2fp和3fp）的传递，且相互作用主要表现为实部（斯托克斯型非线性）。 * 低频（次重力波频段）的波分量仍与短波波群“捆绑”，传播速度接近波群的传播速度。
逻辑关系与贡献：此结果表明，非线性效应在距离破碎点尚远的区域就已对高频分量产生重要影响，且频谱宽度（宽带vs窄带）显著影响了非线性能量传递的效率和方式，从而导致了不同的弥散特性。
阶段III：破碎点附近的非线性浅水变形阶段（Ur ~ 0.7） 主要结果： * 高频处的波数和相速度谱变得对频率的依赖性更弱，因为大部分高频能量被“捆绑”在主频分量上。 * 非线性的振幅效应（amplitude dispersion）对主频分量本身的影响开始显现，导致κ(fp)和c(fp)也开始偏离线性预测值。 * 双谱中虚部（与不对称性相关）开始增强，但波形的水平偏斜（skewness）仍是主导特征。
逻辑关系与贡献：此阶段标志着非线性效应从仅影响高频扩展到也开始影响主频，弥散特性进一步向非色散状态过渡，为进入破波带做准备。
阶段IV：破波带阶段（内破波带，Ur ~ 3.2） 主要结果： * 波数κ(f)和相速度c(f)谱在整个短波频段内几乎不随频率变化（即非色散）。 * 所有波分量都以略大于经典浅水波速（√gh）的速度传播。研究发现，关系c = √[gh(1+ε)]能很好地描述观测到的相速度，其中ε=Hrms/h。 * 双谱显示实部和虚部均较弱且接近，双相位接近-π/2，反映了破碎波浪前倾（俯冲）的不对称形态。
逻辑关系与贡献：此阶段验证了在强烈非线性、深度受限的破波带内，波浪的传播特性主要由振幅效应控制，可用一个修正的浅水波速来描述。
4. 非线性能量传递对κ和c的影响机制 通过Boussinesq模型量化非线性能量传递，研究发现： * 窄谱条件（A3）促进了更高效的非线性能量向2fp和3fp的传递。在阶段II和III，分别有约70%/80%和>90%的2fp和3fp能量是强迫的。这使得其相速度与fp处的相速度锁定一致。 * 宽谱条件（A2）下，非线性能量交换更为复杂，向高频传递的能量相对比例较低（阶段II/III时2fp约为20%/30%，3fp约为30%/50%），其相速度值因此介于强迫波和自由波的特征速度之间。 * 优势波数κ和相速度c的偏离程度，直接与特定频率处强迫能量所占的相对比例相关。强迫能量比例越高，偏离线性弥散关系的程度就越大。
结论： 本研究利用高分辨率实验室数据，系统地揭示了单向不规则波浪在缓坡海滩上从深水传播到破波带全过程中弥散特性的演变规律。研究识别出四个具有不同特征的传播阶段，并定量证明了非线性三波相互作用产生的强迫能量是导致高频波浪弥散特性显著偏离线性理论预测的根本原因。频谱宽度通过影响非线性能量传递的效率和分布，决定了这种偏离的具体模式。在破波带，波浪呈现非色散特性，其相速度可由一个包含非线性参数的修正浅水公式准确描述。
科学价值与应用价值： * 科学价值：深化了对近岸波浪非线性变换物理过程的理解，特别是明确了非线性效应开始显著影响波浪弥散特性的具体位置（对应厄塞尔数约0.3的区域，此时κph约0.6-0.7），并定量揭示了强迫能量比例与弥散特性偏离之间的直接关系。 * 应用价值： 1. 对基于视频/雷达的水深反演算法（如cBathy）的启示：此类算法通常基于线性弥散关系。研究表明，在非线性显著区域（远在平均破碎点之前），若使用高频（>1.5fp）数据反演，由于强迫波相速度的偏离，会导致反演出的相速度被高估20-100%，进而引起水深反演误差。因此，在非线性预期较强的区域，应优先使用峰值频率附近的数据对。 2. 对水下压力数据重建海面高程的启示：这是近岸现场研究的常用方法。线性校正公式会因高估高频波数（κ）而导致高频能量“爆炸式”放大。本研究结果强调了在浅水区（μ ≲ 0.3）需要使用弱色散非线性校正方法，而在更深水域，则需要波浪场空间结构信息才能进行精确的非线性完全色散重建。
附录部分提供了对所用分析方法有效性和鲁棒性的验证：附录A通过合成数据验证了交叉谱分析给出的确实是能量平均意义上的“优势”波数和相速度；附录B通过对比双谱计算与统计方法得到的偏度和不对称性，验证了双谱计算的正确性；附录C讨论并展示了不同双相干（bicoherence）定义的结果，进一步支持了关于非线性耦合强度的分析。这些内容增强了整个研究方法论的严谨性和可信度。
（附表1：试验波浪参数） | 试验 | 有效波高 Hs [m] | 峰值频率 fp [Hz] | 谱峰增强因子 γ | | :— | :— | :— | :— | | A1 | 0.10 | 0.633 | 3.3 (宽谱) | | A2 | 0.20 | 0.444 | 3.3 (宽谱) | | A3 | 0.10 | 0.444 | 20 (窄谱) |