地震勘探中的全波形反演方法
近几十年来，地球科学产业，特别是石油天然气勘探、地热能和地下地质研究领域，地震成像技术取得了飞速发展。过去，地下构造解释主要依赖于相对简单的反射映射和速度分析，而如今，更先进的方法正日益占据主导地位。全波形反演（FWI）被认为是构建高分辨率地下速度模型的“黄金标准”。该方法利用地震波中包含的全部信息，而不仅仅是走时或有限的反射振幅。因此，FWI 是提高地下模型精度和降低解释不确定性的有力工具。
全波形反演的基本概念
简而言之，全波形反演（FWI）是一种非线性反演过程，旨在找到地下参数（例如，P波速度、S波速度、密度或其他弹性参数）的模型，使得由该模型计算得到的合成地震数据与现场观测到的地震数据相匹配。与仅利用部分波信息（例如，层析成像中拾取初至）的传统方法不同，FWI 力求匹配完整的波形：相位、振幅以及波事件之间的干涉。
全波形反演（FWI）通常被建模为一个优化问题：最小化观测数据与合成数据之间的误差。这种误差可以用目标函数来表示，例如两个数据集差异的最小二乘范数。最小化过程是迭代的，逐步更新模型，直到误差达到最小值或满足收敛准则。
为什么全波形反演在地震勘探中如此重要？
地下成像很大程度上取决于速度模型的质量。在反射地震中，速度模型决定了偏移精度、反射层位置、图像清晰度，甚至油气圈闭的解释。如果速度模型不正确，构造图像可能会发生偏移、倾斜，甚至出现伪影，最终误导勘探决策。
全波形反演 (FWI) 通过生成更符合波传播物理规律的精细模型来帮助解决这个问题。与走时层析成像相比，FWI 具有更高的分辨率，可以捕捉到细微的横向和纵向速度变化，包括以下复杂区域：
– 高速对比的盐体
– 近地表复杂性（风化层、空心碳酸盐岩）
复杂的断层和褶皱构造
– 薄储层或中等尺度非均质性
其结果是得到了更清晰的地下模型，有助于迁移过程（RTM/LSRTM），并提高了地质解释的可靠性。
FWI 的工作原理和阶段
一般来说，全波形反演（FWI）工作流程包括以下几个主要阶段：
1. 初始模型的准备
全波形反演(FWI)高度依赖于初始模型的质量。该模型可以来源于层析成像、地质解释、偏移速度模型或多种方法的组合。如果初始模型与实际情况偏差过大，则可能导致反演无法收敛。
2. 正演建模（波浪模拟）
基于当前模型，利用波动方程（声学/弹性）模拟地震波。输出结果为合成地震图。
3. 不匹配计算
拟合误差计算为合成数据与观测数据之间的差异。比较可以在时域、频域甚至特定属性域中进行。
4. 梯度计算（伴随状态法）
为了更新模型，需要了解参数变化的方向以减少拟合误差。伴随状态方法即使对于非常大的参数数量也能高效地计算梯度。
5. 更新模型（迭代优化）
使用优化算法（例如最速下降法、共轭梯度法或拟牛顿法（L-BFGS））更新模型。此过程重复多次，直至结果令人满意。
在实际应用中，全波形反演通常采用多尺度策略：反演从低频开始（以捕捉较大的速度变化趋势），然后逐步增加高频分量以提高细节分辨率。这种策略对于避免陷入局部解陷阱至关重要。
主要挑战：非线性与跳周期
全波形反演（FWI）面临的最大挑战之一是非线性，这会导致目标函数曲面存在许多局部最小值。最著名的问题是周期跳跃，即合成值和观测值之间的差异超过半个波周期，导致反演“卡住”，最终使模型得到错误的解。
周期跳跃通常发生在以下情况下：
初始模型与实际模型相差太远
低频信号不可用或质量差。
– 采集几何形状不足（缺少偏移/方位角）
信噪比低或数据处理不当
为了降低这种风险，采用了多种方法，例如使用低频数据、选择替代目标函数（基于包络或走时误差）、正则化以及整合地质信息和约束。
全波形反演类型：声学型、弹性型和各向异性型
FWI 有多种类型。选择哪种类型取决于媒介的复杂程度和研究目标：
– 声学全波形反演：假设介质为声学介质（不含剪切波）。这种方法更简单，更常用作初始方法，尤其适用于以纵波为主的海洋环境。
– 弹性全波形反演：考虑P波和S波以及模式转换。适用于复杂的陆地和近地表数据，但计算成本更高，且对噪声更敏感。
– 各向异性全波形反演：考虑了各向异性效应（VTI、TTI），这些效应在层状沉积物或具有特定地质构造的区域中非常重要。这提高了精度，但也增加了需要反演的参数数量。
在现代勘探中，趋势是朝着越来越现实的反演方向发展：弹性反演和各向异性反演，尤其是在最终目标是储层表征和岩石性质的情况下。
数据和计算需求
全波形反演（FWI）需要高质量的数据、良好的偏移距覆盖范围和严格的处理控制。数据必须经过去噪、源校正、去鬼波（针对海洋数据）、静校正（针对陆地数据）以及振幅/相位调整等步骤，以确保合成数据与观测数据之间具有可比性。
从计算角度来看，全波形反演（FWI）的计算成本非常高，因为正向和伴随建模需要跨越多个信号源并进行多次迭代。因此，FWI 的实现几乎总是利用集群或 GPU 上的并行计算。随着高性能计算（HPC）技术的进步，这种方法在工业应用中正变得越来越实用。
应用及对勘探决策的影响
全波形反演已被证明可以提高速度模型的质量，进而对以下方面产生直接影响：
– 更集中的迁移和更精确的反射器定位
– 改进了构造解释（断层、盐丘侧翼、河道）
降低钻井目标深度的不确定性
– 更好地与井数据、垂直地震剖面和区域地质数据整合
对于勘探公司而言，深度和图像精度方面的微小改进可以带来巨大的成本节约，并降低钻井干孔的风险。
关闭
全波形反演（FWI）是现代勘探地震工程领域最重要的突破之一，它充分利用了全面的地震波信息。FWI采用非线性反演方法和精细的波物理建模，能够生成高分辨率的地下模型，这对于地质成像和解释至关重要。尽管面临诸如周期跳跃、高质量数据需求以及计算开销巨大等挑战，但算法的开发、多尺度策略的进步以及并行计算技术的进步，使得FWI技术日趋成熟并得到广泛应用。在勘探领域，该方法不仅提高了地下图像的质量，而且为更准确、高效、低风险的决策提供了关键基础。
