全波形反演：把地下速度模型，用地震波一点点磨准

地震勘探里有个绕不开的东西，叫速度模型——地下每一层岩石里地震波跑多快，都得先摸清。它是后面成像、归位、看构造的地基：速度算歪一点，井位就可能打偏一截。可地下几千米，谁也没法拿尺子去量。那这张"速度地图"到底怎么来的？ 有一种如今被视为最高精度的办法，专门干这件事，它叫全波形反演。

一句话先给答案

全波形反演（FWI，Full Waveform Inversion）干的活儿，是拿一张假想的速度模型去"正演"模拟出理论上该收到的地震波形，再和野外实测的地震记录逐点对比、算出差异（残差），然后反复修改速度模型，让模拟波形一步步逼近真实波形——差异磨到足够小，那张模型就足够准了。 说白了，它不是一次算出答案，而是"猜一版、对一遍、改一版"地反复迭代，直到理论和现实对上号。之所以叫"全波形"，是因为它不像老办法只挑走时或振幅这几个特征，而是把地震波形里的全部信息——连相位都不放过——都拿来拟合。

它是怎么"对答案"的

!OGC-B-033 封面 <!-- 【配图占位符·主视觉】金#C6A24C + 靛#3B3A52；抽象科技型：左侧一版粗糙的地下速度模型，中间正演模拟出的合成地震波形与实测波形叠合比对、标出残差，右侧经反复迭代后收敛出的高分辨率速度模型；金色引线示意"正演→残差→迭代→高精度模型"，不出现地震灾害/预测意象 -->

拆开看，其实就是一个闭环在转圈。第一步正演：给定一版速度模型，用波动方程算出地震波在这个"假想地下"里怎么传，会在地面收到什么样的波形。第二步比对：把这份合成波形和真实观测记录摆在一起，两者的差就是残差。第三步反演更新：顺着残差往回推，算出该把模型的哪块速度调高、哪块调低，改出一版更好的模型。然后拿新模型再正演、再比对、再更新……一圈圈迭代下去，残差越缩越小，速度模型也越磨越贴近真实的地下。 有点像调一台老收音机——你不知道台在哪，只能一点点拧旋钮，听着杂音一次次变小，直到声音清亮，才知道对上了频率。

为什么它能比老办法更准

关键就在"全波形"这三个字。传统的层析成像这类方法，往往只用地震波的走时信息，等于只看波"什么时候到"，好算、稳，但地下的细节被大量丢掉了。全波形反演不一样——它把波形的振幅、相位、每一处起伏都当成约束，信息用得满，反演出的速度模型分辨率自然高得多。正因为吃进了全部波形信息，它才有本事刻画复杂构造、盐下、以及那些小尺度的速度变化，把老办法糊成一团的地方看清楚。地下越复杂、越难啃，它的优势越明显。至于这套"正演—比对—迭代"的链条在实际工区怎么搭、初始模型怎么准备，我们整理了一份报告，放在文末。

高精度不是白来的：它有三个软肋

好东西都有代价，全波形反演的代价还不小。

其一，算力惊人。 每迭代一轮都要在整个模型上做一次正演模拟，反复几十上百轮，数据量和计算量堆起来极其庞大，通常得靠超级计算机才跑得动，慢、也贵。

其二，怕初始模型不准。 这是它最要命的软肋。反演本质是在"摸黑找最低点"，如果一开始给的初始模型离真实地下太远，它很容易走错方向、掉进一个假的最优解爬不出来——这在行业里叫局部最优。

其三，周波跳跃（cycle skipping）。 这是初始模型不准最典型的翻车方式：当模拟波形和实测波形错开超过半个周期，程序会把模拟波的这一个波峰，错认成实测波的下一个波峰去对齐，于是模型越改越偏，一路错到底。要躲开它，就得有足够好的初始模型、足够低频的数据，把起点摆到离真身够近的地方。所以全波形反演对数据质量的要求，也比一般方法苛刻得多。

说到底

全波形反演的本分，就是拿地震波形的全部信息反复迭代拟合，把一张粗糙的初始速度模型，一遍遍磨成高分辨率的地下速度模型。 它精度高，是因为一点信息都舍不得丢；它算力大、怕起点不准，也正是因为要拟合的东西太满、要走的路太长。从只看走时的老办法，到连相位都不放过的全波形，物探这些年在"把地下看真切"这条路上较的劲，说到底是同一件事：让屏幕上那张速度模型，尽可能贴近真实的地下。

一张能把地下速度看准的模型，值得被讲清楚。

想把全波形反演这套"给地下速度建模"的门道彻底看明白？ 我们做了份《全波形反演高精度速度建模与复杂构造成像研究》——公众号后台回复「地震」，拿去看。