Morison 式波浪荷载模拟器 返回
海洋工程·波浪荷载

Morison 式波浪荷载模拟器 — 海洋构造物的惯性力·阻力

为夹套平台或海上风电单桩基础等细长海洋构造物,用 Morison 式(惯性力+阻力)进行波浪荷载的实时计算工具。改变有义波高、周期、水深、构件直径、C_D/C_M,可立即获得波长、水质点速度、KC 数、全波浪荷载和倾覆力矩。

参数设置
有义波高 H_s
m
由谱定义的代表波高(最大波的约 1.6 倍)
峰值周期 T_p
s
JONSWAP 频谱的峰值周期
水深 d
m
构件直径 D
m
构造形状
代表性 C_D / C_M 范围因形状而异
阻力系数 C_D
粗糙 1.0-1.2/光滑 0.65/海洋生物附着时增加
惯性系数 C_M
理想势流为 2.0,光滑面为 1.6
没水高度 z
m
波浪荷载积分的有效水中长度
计算结果
波长 λ (m)
最大水质点速度 u_max (m/s)
阻力/单位长 (kN/m)
惯性力/单位长 (kN/m)
全波浪荷载 (kN)
KC 数
海面与没水构件 — 水质点轨迹·力矢量

传播中的线性波与没水圆柱(夹套腿)。蓝箭头=惯性力 F_i,红箭头=阻力 F_d,白点=水质点轨迹。KC 数表示主控模式。

波周期内的荷载时程 — F_drag·F_inertia·F_total
KC 数 vs 构件直径 D — 域判定
理论·主要公式

$$F = \rho C_M V \dot u + \tfrac{1}{2}\rho C_D A |u|u,\qquad KC = \frac{u_{max} T}{D}$$

第一项=惯性力(与流体加速度成正比),第二项=阻力(与速度平方成正比,符号由 |u|u 保持)。ρ=海水密度 1025 kg/m³,C_M=惯性系数(粗糙圆柱 2.0),C_D=阻力系数(粗糙圆柱 1.0-1.2,光滑 0.65),V=单位长排除体积,A=单位长投影面积。

$$u_{max}=\omega\,\tfrac{H_s}{2},\qquad \dot u_{max}=\omega^{2}\,\tfrac{H_s}{2},\qquad \lambda=\tfrac{2\pi}{k},\quad k=\tfrac{\omega^{2}}{g}\ (\text{deep water})$$

线性波理论(Airy 波)给出的水面(z=0)处水质点速度和加速度的最大值。深海近似 k=ω²/g,ω=2π/T。Stokes 5 阶等高阶理论会使峰值增加 10-30%,但与 KC、C_D 选取的不确定性相比属次要。

$$F_{total}/L = \sqrt{F_{drag}^{2}+F_{inertia}^{2}}\quad (\text{quadrature sum at peak})$$

惯性力与阻力相位相差 90°,故其峰值的平方和平方根接近时程中的最大荷载值。严格上应按时程评估合成波形。

基于 Morison 式的海洋构造物波浪荷载计算

🙋
「Morison 式」在海上风电和石油平台设计中经常听到,到底是在算什么呢?
🎓
简单说就是「海水中的细长柱子,波浪施加的力有多大」的计算公式。1950 年 Morison、O'Brien、Johnson、Schaaf 四人在墨西哥湾石油平台腿柱数据基础上,大胆地把流体阻力理论和附加质量理论直接相加,得出了 F = ρ·C_M·V·du/dt + ½·ρ·C_D·A·|u|·u。第一项是「流体加速时被拖拽的惯性力」,第二项是「流速快时后面产生漩涡的阻力」。
🙋
原来如此,流体的两种效应相加。但我注意到左边滑块把「构件直径 D」拉到 10 m 时,KC 数变得很小,显示「惯性支配」,这是什么意思?
🎓
好问题!KC = u_max·T/D 表示「一个波周期内水质点的运动距离是构件直径的多少倍」。D 越大 KC 越小。比如 KC=2 时,水质点只运动了柱子直径的 2 倍,柱子后面的漩涡还没来得及充分发展,所以「加速水流的惯性」比「漩涡的阻力」更重要。反过来细管子(立管、KC=50)水质点就运动了 50 个直径,漩涡充分产生,阻力就主导了。
🙋
那夹套平台的腿柱(D=2m 左右)是不是正好在中间呢?
🎓
对的,默认条件下 KC≈9.4,正好在过渡域 5 < KC < 20 内。下面的「KC 数 vs 构件直径」图表可以看到夹套腿就在这个带状区间内。这个过渡域其实最麻烦,因为 C_D 和 C_M 会随 KC 剧烈变化。Sarpkaya 的摆荡水槽实验表明 KC=10 附近 C_D 能跳到 1.4,所以 API RP 2A 在粗糙带海洋生物的情况下推荐 C_D=1.05、C_M=1.2 这样比较保守的值。
🙋
提到海洋生物了,这在设计里要考虑多少呢?
🎓
影响很大。北海夹套平台在运行 20 年后,贻贝和水螅虫可能长到 5~10 cm 厚,实际直径增加、表面粗糙度暴增,C_D 从光滑的 0.65 增加到 1.2 左右,翻了接近一倍。设计实务中会明确写「海洋生物层厚 50 mm、密度 1325 kg/m³、C_D=1.2」这样的条件来评估寿命中后期的恶化荷载。日本沿岸藤壶生长快,海上风电单桩的过渡段(TP)需要定期清洁防附着。
🙋
最后,计算出来的「全波浪荷载 293 kN」对一根夹套腿来说大不大?
🎓
没水长 50 m 的单根腿柱来说这是典型数量级。4 腿夹套的话,垂直于波的两根腿同时承受最大荷载,底部剪力大约 500~600 kN。再乘以离底部几十米高的臂长,就是 10~15 MN·m 的倾覆力矩。真正的北海夹套设计中,100 年重现波(H_s≈14 m, T_p≈15 s)会产生底部剪力 50~100 MN、倾覆力矩 3~5 GN·m,这才是设计控制荷载。本工具只是定阶段用,真正设计要用 SACS、USFOS 这样的软件做时程非线性分析。

常见问题

Morison 式是针对「相对波长构造物直径足够小(约 D < λ/5)的细长构件」的标准经验公式。具体包括固定式夹套平台的腿柱、水平斜撑、导管、立管、海底管道、海上风电单桩基础等。而对于大直径重力式结构物(GBS)或半潜式浮体的浮筒等,当直径与波长同阶时,绕射效应不能忽略,应使用绕射理论(MacCamy-Fuchs 等)或 3D 绕射分析。
对于粗糙圆柱(夹套腿),设计实务中通常采用 C_D ≈ 1.0~1.2、C_M ≈ 2.0,这是 API RP 2A WSD 和 DNV-OS-J101 的共同建议。光滑表面时分别约为 C_D ≈ 0.65、C_M ≈ 1.6。这些值是 Reynolds 数 Re 和 Keulegan-Carpenter 数 KC、表面粗糙度的函数,由 Sarpkaya、Isaacson 等的实验数据确定。海洋生物附着物(藤壶、贻贝)增长时,C_D 可增加到 2 倍,需要单独考虑寿命中期的恶化荷载。
KC = u_max·T/D 表示「一个波周期内水质点运动距离是构件直径的多少倍」,反映阻力和惯性力哪个起主导作用。KC < 5 时水质点振幅较小,流体对构件的「摇晃」中不会形成发达的剥离涡,惯性力(C_M 项)起主导(大直径圆柱、单桩)。KC > 20~30 时有足够的振幅使涡脱落发育,阻力(C_D 项)起主导(立管、管线)。5 < KC < 20 是过渡区,两者同阶,是设计中最需注意的区域。
判断标准是「波长 λ 与构件直径 D 的比值」。当 D/λ < 0.2(约 1/5)时,Morison 式(细长物体近似)足够。例如本工具默认条件(H_s=6 m, T_p=12 s)中 λ≈225 m,D=2 m 的夹套腿 D/λ≈0.009 很小,Morison 式适用。当 D/λ > 0.2 的大直径结构(GBS 重力式基座、储油罐、8~10 m 级巨大单桩)时,构造物对波场本身的绕射、反射效应显著,需用 MacCamy-Fuchs 解析解或边界元法(BEM、WAMIT/HydroD)进行绕射分析。

实际应用

固定式夹套平台(墨西哥湾、北海):壳牌、埃克森美孚、挪威油气公司运营的石油天然气生产平台大多采用钢制夹套结构,腿柱(直径 1~2 m)和水平、斜撑是 Morison 式的主要应用对象。针对 100 年重现波(北海中部 H_s≈14 m, T_p≈15 s)进行 SACS 或 USFOS 时程非线性分析,确定底部剪力和倾覆力矩,设计夹套腿和桩基。运行超过 50 年的老化平台会经历比设计时更恶劣的波浪条件,需重新用 Morison 式评估是否适合延寿。

海上风电单桩基础(北海、日本沿岸):Hornsea、Walney Extension 等全球大型海上风电场主要采用直径 7~10 m 的钢管单桩基础。直径大导致 KC < 5 的惯性支配域,使用「修正 Morison 式」加上 MacCamy-Fuchs 绕射修正系数评估。日本千叶、秋田等商业海上风电也按 DNV-OS-J101/IEC 61400-3 用同样方法设计。FAST、Bladed 等软件分析风与水的动力耦合作用(aero-hydro-elastic coupling)。

浮体式海上风电和浮式生产平台的系泊系统:Hywind Scotland(全球首座商业浮体式风电场)和日本福岛浮体式实证采用圆筒型或半潜式浮体搭载风机。浮体的系泊缆线(链条直径 100~150 mm、聚酯绳)也要用 Morison 式评估波浪和潮流荷载,用 OrcaFlex 软件做系泊张力动力分析。Vestas、三菱-维斯塔斯等从实证向量产过渡,Morison 计算是系泊索疲劳寿命评估的基础数据。

海底管道和立管设计:大水深油气田的海底管线(直径 0.3~0.6 m)和刚性立管在底层流和内部波作用下长期承受 Morison 荷载。构件细、KC 值大,阻力主导,涡激振动(VIV)是长期疲劳破坏的主要模式。Shear7、VIVA 等专用软件基于 Morison 做模态分析,通过加装防护条(strake)或调整张力来应对。巴西 Petrobras 深海油田是典型案例。

常见误区和注意点

最大陷阱是「把 C_D 和 C_M 当成常数固定」。实际上它们随 Reynolds 数 Re、KC 数和表面粗糙度 k/D 大幅变化。Sarpkaya 的经典摆荡水槽实验显示 KC=10 附近 C_D 从 0.6 能跳到 1.4,KC=30 时又平均落到 0.8,行为很复杂。北海夹套再评估案例中出现过「设计时 C_D=0.7,实测改为 C_D=1.1 后荷载增加 60%」的情况。设计规范(API RP 2A、DNV-OS-J101)推荐的保守值最稳妥,特殊条件下应做水槽试验或 CFD 单独验证。

其次,「只考虑水面波」的问题。本工具用线性波(Airy 波)只算水面 z=0 的最大值,但实际上水质点速度沿水深指数衰减(深海 exp(-kz))。50 m 深夹套腿在水面最大速度 u_max=1.57 m/s,到 30 m 深衰减到 0.43 m/s(=1.57·exp(-0.028·30))。设计实务要沿水深方向积分分布荷载,本工具的「全波浪荷载 = 峰值荷载 × 没水长」只是上界估算,SACS/USFOS 自动进行水深积分。

最后,「Morison 式深海浅海不分」的错误。本工具用深海近似 k=ω²/g,但浅海(d/λ < 0.05)差很大。浅海中 λ ≈ T·√(g·d),波长变成水深函数。例如水深 d=15 m、T_p=10 s 的浅海条件,深海近似给出 λ=156 m 而实际 λ≈115 m(误差 30%),u_max 也偏高。浅海和中间水深设计要用通用分散关系 ω² = g·k·tanh(k·d),用 Newton-Raphson 法求解。沿岸防波堤、海上风电近岸场址必须用浅海效应,否则荷载严重高估。

使用指南

  1. 输入有义波高(Hs)和周期(Tp)。北海油田典型值为 Hs=6m、Tp=12s。
  2. 设置水深和构件直径(如单桩外径)。海上风电通常水深 30~50m、直径 8~12m。
  3. 调整阻力系数 CD(0.6~1.2)和惯性系数 CM(1.5~2.0),实时检查波长λ、水质点速度 u_max、单位长阻力和惯性力。
  4. 若 KC 数 < 5 说明惯性力主导,> 15 说明阻力主导,5~15 是过渡区。

具体计算例

夹套平台(水深 40m、构件直径 1.2m、Hs=5m、Tp=10s):波长λ≈156m,最大水质点速度 u_max≈1.85m/s,KC≈18.5(阻力型)。取 CD=1.0、CM=1.6,单位长阻力=14.2kN/m,惯性力=6.8kN/m。全波浪荷载由有效长度积分,乘以安全系数 1.35 得设计荷载(约 152kN)。

实务注意事项