3、冲刷深度预测：经验公式法、数值模拟与物理模型试验

说到冲刷深度预测，我这些年打交道最多的就是这三种方法。经验公式、数值模拟、物理模型试验，各有各的脾气。说白了，没有哪个是万能的，关键看你的项目阶段和精度要求。

我个人习惯是：初步设计用经验公式快速估算，详细设计用数值模拟优化，关键工程或复杂海况下，必须上物理模型试验验证。嗯，咱们一个一个聊。

3.1 经验公式法：快速估算的利器

经验公式，说白了就是前人用大量试验和实测数据“凑”出来的关系式。优点是快，缺点嘛——适用范围有限。我在项目中最常用的是 HEC-18 和 DNV 公式。

3.1.1 HEC-18 公式

这个公式来自美国联邦公路局，主要用于桥墩冲刷。我最早接触它是在一个跨海大桥的初步设计阶段，当时时间紧，就用它快速算了一版。

核心表达式：

ys / a = 2.0 * K1 * K2 * K3 * (y0 / a)^0.35 * Fr^0.43

其中：

• ys：冲刷深度（m）
• a：桩径（m）
• y0：水深（m）
• Fr：弗劳德数
• K1, K2, K3：修正系数（桩型、来流角度、床沙粒径）

3.1.2 DNV 公式

DNV（挪威船级社）的公式在海洋工程里更常见。它考虑了波浪和流的联合作用，这一点比 HEC-18 更贴近海上风电场景。

DNV 推荐的形式：

S / D = 1.3 * (1 - exp(-0.03 * (KC - 6)))   （当 KC > 6 时）

其中：

• S：平衡冲刷深度
• D：桩径
• KC：Keulegan-Carpenter 数（反映波浪作用强弱）

3.2 数值模拟方法：精细化分析

数值模拟，说白了就是用计算机“算”出流场和冲刷坑的演化过程。我常用的工具是 FLOW-3D 和 OpenFOAM。

为什么要用数值模拟？因为经验公式只能给一个最终深度，而数值模拟能告诉你：

• 冲刷坑的形态（是V形还是U形？）
• 随时间的变化过程（什么时候达到平衡？）
• 不同防护方案的效果对比

3.2.1 基本流程

1. 建立几何模型：桩基、海床、边界条件
2. 设置网格：近壁面加密，远场稀疏
3. 选择湍流模型：我一般用 k-ε 或 SST k-ω
4. 泥沙输运模型：考虑悬移质和推移质
5. 动网格或地形更新：每步更新床面高程

这里给一个简单的 FLOW-3D 设置片段（示意）：

# 泥沙参数设置
sediment_density = 2650  ! kg/m3
median_diameter = 0.0002 ! 200 μm
critical_shields = 0.03  ! 临界希尔兹数
bed_load_coeff = 5.7     ! 推移质系数

3.2.2 数值模拟的局限性

说实话，数值模拟也不是万能的。它有几个硬伤：

• 计算时间长：一个 3D 模型跑 24 小时是常事
• 参数敏感：泥沙参数稍微调一下，结果差一倍
• 无法模拟极端工况：比如台风浪下的冲刷，数值模型容易发散

3.3 物理模型试验：最可靠的验证手段

物理模型试验，说白了就是造一个缩尺的水池，放个模型桩进去，用水泵造流、造波机造浪，直接看冲刷效果。这是我最信任的方法，没有之一。

为什么？因为物理模型能真实反映水动力-泥沙-结构三者的耦合作用。数值模拟里那些简化假设，在物理模型里都不存在。

3.3.1 相似准则

做物理模型，最核心的是相似准则。我常用的有：

相似准则	控制参数	适用场景
弗劳德相似	Fr 数	重力主导的流动（波浪、明渠流）
雷诺相似	Re 数	粘性力主导（但缩尺后很难满足）
希尔兹相似	θ 数	泥沙起动条件

3.3.2 试验步骤

1. 模型制作：桩基用有机玻璃或金属，床沙用天然沙或轻质沙
2. 率定：先不放桩，率定流速和波高
3. 放桩试验：记录冲刷深度随时间的变化
4. 数据采集：用声呐或激光扫描冲刷坑形态
5. 结果分析：换算到原型尺度

3.3.3 物理模型的坑

嗯，这里我要多说几句。物理模型看着简单，其实坑很多：

• 缩尺效应：小尺度下泥沙的粘性行为和大尺度不一样
• 时间缩尺：冲刷过程的时间缩尺很难准确
• 边界效应：水池壁面会影响流场

3.4 三种方法的对比与选择

我把这三种方法的特点整理成了一张图，方便你理解它们的关系：

你看，这三种方法不是互相替代的关系，而是递进的关系。我个人的工作流是：先用经验公式算个大概，心里有个数；然后用数值模拟细化，看看冲刷坑形态和时变过程；最后如果项目重要或者海况复杂，一定上物理模型。

说到最后，我想强调一点：不管用哪种方法，都要留有余量。海洋环境太复杂了，你永远不知道下一个台风会带来什么。我一般会在预测值上加 20% 的安全裕度，这是用真金白银换来的教训。