1. 课程导入与波浪谱概念建立

海洋工程、船舶耐波性、海岸动力分析以及海洋能评估,都离不开对海浪的定量描述。然而真实海面是随机起伏的三维过程,任意一点的波面高度随时间变化,既没有严格周期,也没有固定振幅。若只记录一条波高时程曲线,我们很难直接回答“这片海域的能量主要集中在哪些频率上”“哪种波长的波浪对平台影响最大”这类问题。波浪谱正是为解决这些问题而引入的数学工具:它把随机波面从时域转换到频域,用能量分布的方式刻画海况。

学习目标

完成本章学习后,应当能够:说明波浪谱的物理意义,理解谱密度与波面方差、有效波高之间的关系;区分频率谱与波数谱、一维谱与方向谱的适用场景;读懂常见谱图(谱密度曲线、方向分布极坐标图)的坐标含义;并建立对后续谱估计、谱拟合与可视化流程的整体认识。本章不要求立即动手计算,但需要把“谱是什么、为什么这样定义、图上每个量代表什么”这三件事想清楚,因为后续所有谱分析代码都建立在这些概念之上。

波浪谱的物理意义

把随机波面 η(t) 看作由无数不同频率、不同相位的正弦波叠加而成,是谱概念最直观的出发点。设第 i 个成分的角频率为 ωi,振幅为 ai,则波面可写成各成分之和。单个成分携带的能量正比于振幅平方,于是把所有成分按频率排列,就得到能量随频率的分布。定义谱密度 S(ω) 使得 S(ω)dω 表示频率区间 [ω, ω+dω] 内波面成分的方差贡献,则波面总方差(即波面高度平方的期望)满足:

σ² = ∫₀^∞ S(ω) dω
H_s ≈ 4σ = 4·sqrt(∫₀^∞ S(ω) dω)

这里 σ 是波面标准差,Hs 为有效波高。这个关系是波浪谱最核心的实用价值:只要拿到谱曲线,积分一次就能得到有效波高,无需逐波统计。需要强调,谱密度不是“某个频率上真实存在一个振幅为 S(ω) 的波”,而是能量在频率上的分布密度,单位通常为 m²·s/rad 或 m²·s。谱曲线下的面积才有物理量纲意义,单点纵坐标不能直接当作波高读。

另一个常被忽略的要点是谱矩。定义 n 阶谱矩 mn = ∫ωnS(ω)dω,则 m0 就是方差,平均周期可用 m0/m1 或 sqrt(m0/m2) 等形式表达。谱形越宽,能量分散在越多频率上,波面越不规则;谱形越窄,越接近单一周期的规则波。这一“宽窄”特征在后续判断海况成熟度时非常有用。

频域与波数域表示

频率谱描述的是固定点上波面随时间的起伏,自变量是角频率 ω 或频率 f,二者关系为 ω=2πf。工程上常用 f(Hz)作图,因为周期 T=1/f 更直观;理论推导多用 ω。绘图时必须标明用的是哪一种,否则横轴数值会相差 2π 倍,这是初学者最容易出错的地方之一。

波数谱则把自变量换成波数 k=2π/L(L 为波长),描述的是空间上的起伏结构。在深水条件下色散关系为 ω²=gk,频率与波数一一对应,因此频域与波数域可以互相转换。波数域在分析波浪对结构物的空间作用、折射绕射以及遥感反演时更自然;而单点观测(如浮标)天然给出的是频率谱。理解两者的对应关系,才能在不同数据源之间正确换算。

真实海面并非单向传播,能量还随传播方向分布,于是有方向谱 S(f,θ)。它在一维频率谱基础上增加方向维,满足 ∫S(f,θ)dθ = S(f)。方向分布通常用集中度参数或 cos2s(θ/2) 型分布描述。方向谱的可视化常用极坐标图或方向-频率二维色图,这是本课程数据可视化部分的重要内容。

波浪谱:从随机波面到能量分布 时域:随机波面 η(t) 单点记录,难以直接读出能量结构 傅里叶 变换 频域:谱密度 S(ω) 面积 = 方差 σ²,H_s ≈ 4σ ω²=gk 深水色散 波数域:S(k) 空间结构,用于绕射与遥感反演 方向谱 S(f, θ):能量随频率与传播方向联合分布 主波向 θ ∫ S(f,θ) dθ = S(f) 方向集中度反映海况的成熟程度 可视化:极坐标图 / 方向-频率色图

课程整体框架与数据案例

本课程按“概念—估计—拟合—可视化—应用”的主线组织。前半部分建立谱的定义与常用谱型(如 Pierson-Moskowitz 谱、JONSWAP 谱),说明它们各自适用的海况条件与参数含义;中段进入实测数据处理,包括去趋势、去均值、窗函数选择、快速傅里叶变换与谱平滑,重点讲清楚为什么原始周期图需要平滑、平滑会带来什么代价;随后是谱拟合与参数反演,用观测谱估计有效波高、谱峰周期和谱宽参数;最后进入可视化与工程应用,包括谱图规范绘制、方向谱呈现以及谱在结构响应估算中的用法。

贯穿全课程的数据案例设定为一段教学示例:某近海观测点以 1 Hz 采样、连续记录 1024 个波面高度样本,采样时长约 17 分钟。示例中假设海况为单一峰值的风浪,谱峰周期约 8 s,有效波高约 2.0 m。这一组参数仅用于演示计算流程与图形判读,并非实测结果。后续章节会基于这组设定演示:如何从 1024 点序列得到频率分辨率约 0.00098 Hz 的原始谱、如何用分段平均降低谱估计方差、以及如何从平滑谱积分反推有效波高并与设定值对照。读者应重点关注“输入参数—处理步骤—谱形变化—结果校验”之间的对应关系,而不是记住某个具体数值。

学习建议:每遇到一张谱图,先确认三件事——横轴是 f 还是 ω、纵轴是谱密度还是对数谱密度、曲线下面积对应什么物理量。这三项确认清楚,后续读图与计算就不容易混淆。