有限差分正演：复杂介质地震反演的信任基石

简介：这份PPT面向地震勘探领域的地球物理师、地质科研人员及高校相关专业学习者，系统讲解Tesseral 2D模型正演软件的功能与特点，帮助读者理解如何借助正演模拟应对复杂介质下常规处理假设失效、解释结果多解等实际问题。资源包共1个PPT文件，约2MB，以图文并茂的幻灯片形式组织内容，便于课堂讲解与自学浏览。目前已有226人学习关注。内容围绕Model builder、Computation Engine、Viewer与Processing Block四大功能模块展开，涵盖LAS测井曲线建模、基于有限差分的P波与SV波合成记录计算、波场快照可视化，以及速度分析、CMP抽道叠加与叠前叠后偏移处理等环节；同时介绍梯度、声学、弹性及各向异性介质模拟与AVO效应分析，并结合溶洞、砂泥互层等模型正演示例，说明其在采集参数设计、处理流程查证与解释成果验证中的实用价值。

1. 正演不是“画个图看看”：为什么复杂介质下反演必须先有正演兜底

地震勘探拿到的是一张地下反射界面图，但真正决定井位的，往往不是这张图本身，而是“这张图在复杂介质里到底可不可信”。常规处理流程建立在一堆假设上——水平层状、均匀介质、射线路径近似，一旦遇到横向变速、薄互层、各向异性或者溶洞，这些假设就集体失效，反演结果出现多解性，解释出来的构造和实际钻井对不上。Tesseral 2D 这类模型正演软件解决的正是这个信任问题：先给定一个已知的地下模型，用有限差分把波场正着推一遍，看合成记录长什么样，再拿它去校准采集参数、处理流程和解释结论。换句话说，正演是给反演准备的一颗“后悔药”，让你在花钱打井之前，先在计算机里把各种可能性跑一遍。搞懂它的底层机制，比会点按钮重要得多。

2. 有限差分正演引擎：波场是怎么一步步被“推”出来的

2.1 从波动方程到差分网格：为什么必须离散化

地震波在介质中的传播由弹性波方程描述，方程本身是连续的偏微分方程，计算机没法直接求解连续域，必须把空间和时间切成网格。Tesseral 2D 的 Computation Engine 核心就是这套有限差分格式：把二维空间剖分成规则网格，每个网格点存速度、密度等介质参数，时间也切成固定步长，然后按差分格式逐时间步推进波场。常见的实现思路是交错网格有限差分，速度和应力分量放在不同网格位置上，这样能更稳定地处理自由表面和介质界面。

这里有个关键取舍：网格越细、时间步越小，模拟精度越高，但计算量按网格数的平方甚至立方增长。所以实际使用中，网格间距通常取目标波长的一小部分，比如每波长若干个网格点，时间步长则受 CFL 稳定性条件约束——步长太大，数值解会直接发散，波场变成一团噪声。这也是为什么复杂模型跑起来慢，不是软件写得差，是物理规律决定的。

# 简化的二维声波交错网格有限差分推进（示意，非Tesseral源码）
# p: 压力波场, vx/vz: 质点速度分量, dt: 时间步长, dx/dz: 网格间距
for it in range(nt):
    # 更新速度场：压力梯度驱动速度
    vx[1:-1, :] -= dt / rho[1:-1, :] * (p[1:-1, :] - p[:-2, :]) / dx
    vz[:, 1:-1] -= dt / rho[:, 1:-1] * (p[:, 1:-1] - p[:, :-2]) / dz
    # 更新压力场：速度散度驱动压力
    p[1:-1, 1:-1] -= dt * bulk[1:-1, 1:-1] * (
        (vx[2:, 1:-1] - vx[1:-1, 1:-1]) / dx +
        (vz[1:-1, 2:] - vz[1:-1, 1:-1]) / dz
    )
    # 每若干步输出波场快照，供Viewer可视化
    if it % snapshot_interval == 0:
        save_snapshot(p)

上面这段伪代码把机制说清楚了：每一时间步先由压力梯度更新速度，再由速度散度更新压力，交替推进。 rho 是密度， bulk 是体积模量，两者共同决定波速。 snapshot_interval 控制波场快照的输出频率，这正是 Tesseral 里“波场快照技术”的底层来源——不是事后画图，而是在推进过程中按固定时刻把整个波场存下来。原文给出的 t1=0.350 s 到 t5=1.420 s 那组快照，就是按这个节奏截出来的，用来追踪波前在介质中的传播路径。

2.2 P 波和 SV 波：弹性介质里为什么要分开算

声学介质只需要一个压力场，但真实地下是弹性介质，纵波（P）和横波（SV）同时存在，传播速度不同，遇到界面时的反射透射行为也不同。Tesseral 的 Computation Engine 明确支持 P 波和 SV 波的分别计算，这意味着底层用的是弹性波方程而不是简单的声波方程，网格上要同时维护多个应力分量和速度分量。多分量采集和极化分析之所以能做，前提就是正演阶段把两个波场都算出来了。

常见的实现思路是把弹性波方程写成速度-应力一阶方程组，交错网格上速度和正应力、切应力分别布置，边界处理用吸收边界或者阻尼层来压制人工反射。这里最容易翻车的地方是边界条件：如果边界没处理好，模型边缘会产生虚假反射，混进有效波场里，看起来像地下真有个界面。所以正演结果里那些“莫名其妙”的同相轴，先怀疑边界，再怀疑模型。

2.3 介质类型：梯度、各向异性、吸收是怎么进方程的

Tesseral 支持的介质类型不是摆设，每一种都对应方程里参数的改变。梯度介质是把速度写成空间位置的函数，网格上每个点的 bulk 和 rho 随深度或横向变化；各向异性介质则需要把体积模量从标量换成刚度张量，方程形式更复杂，计算量也更大；吸收介质引入衰减机制，通常用品质因子或者粘弹性模型来近似，让波在传播过程中能量逐渐耗散。原文提到“考虑波的能量衰减、初至到达时间场”，初至时间场一般是在推进过程中记录每个网格点首次达到阈值的时刻，用来做走时分析。

这些介质类型的选择直接决定了正演能不能反映真实地质体的响应。比如薄互层模型，如果只用声学介质，就丢掉了横波信息，AVO 效应根本算不准；溶洞模型如果忽略吸收，波场能量分布会偏乐观。所以建模时选哪种介质，不是随手勾选，而是要先想清楚你要验证的是什么物理现象。

3. 建模与可视化：LAS 曲线怎么变成可计算的剖面

3.1 Model builder：从测井曲线到密度-速度模型

Model builder 的职责是把地质剖面的几何结构和物性参数变成 Computation Engine 能吃的二维数组。原文特别提到可以用 LAS 格式的测井曲线文件输入 P 波和 SV 波的速度、密度值，这一点很关键：LAS 是测井行业的标准格式，里面按深度记录了声波时差、密度等曲线，Model builder 要做的是把这些一维深度序列插值、映射到二维网格上，再配合用户画的地层界面，构建出完整的密度-速度模型。

常见的实现思路是分层赋值加横向插值：先定义地层框架，每层指定从哪条曲线取值，然后在层内沿横向做插值或外推，遇到断层、尖灭这类复杂结构再手工调整。这里有个血泪经验——LAS 曲线的深度单位和模型网格的深度单位必须对齐，采样率也要匹配，否则插值出来的速度场会出现阶梯状伪影，正演记录上表现为莫名其妙的周期性反射。建模阶段偷的懒，后面处理阶段都要加倍还回来。

3.2 Viewer：波场快照和合成记录到底在看什么

Viewer 不是简单的看图工具，它承担的是“把计算结果翻译成可判断的信息”这个任务。波场快照让你看到某一时刻波前在介质中的位置和形态，合成炮集记录让你看到检波器实际会收到什么。原文说 Viewer 可以“对比不同的波场分量，提取信号振幅值、相关属性值”，这意味着它支持分量切换和振幅提取，做 AVO 分析时可以直接从合成记录上读振幅随偏移距的变化。

波场快照的浏览方式值得单独说：按固定时刻截取的快照序列，本质上是在时间轴上采样波场，把传播过程拆成若干帧。判断一个反射到底来自哪个界面，靠的就是追踪波前从激发点出发、遇到界面、产生反射、再传到接收点的完整路径。如果只看最终炮集，很容易把多次波或者绕射波误认成一次反射；有了快照，波场的“来龙去脉”就清楚了。这也是原文强调“准确识别勘探目标所产生的波场、追踪其在介质中的传播过程”的原因。

3.3 Processing Block：合成数据为什么要再走一遍处理流程

Processing Block 的存在容易被忽略，但它恰恰是正演闭环的关键一环。它把 Computation Engine 生成的合成炮集当成真实数据来处理：速度分析、CMP 抽道、叠加、叠前叠后时间偏移、叠前叠后深度偏移，一整套流程走下来，看处理结果能不能恢复出当初建模时的地质结构。如果能恢复，说明这套处理流程和参数在类似地质条件下是可信的；如果恢复不出来，说明流程里有环节不适合这种介质。

这背后的逻辑是：正演给你一个已知答案，处理流程是解题过程，最后对答案。常规处理在复杂介质下失效，往往不是某个模块写错了，而是假设条件不成立，而通过合成数据走一遍，能定位到到底是速度分析不准、偏移算法不适配，还是叠加方式有问题。原文说 Processing Block“可以非常方便地对复杂介质估算 AVO 效应”，也是同一个道理——合成记录上振幅随偏移距的变化是已知模型产生的，拿处理后的结果去对比，就能检验 AVO 分析流程的可靠性。

4. 避坑与排查：正演跑出来的结果不对劲，先查这几处

4.1 波场快照边缘出现强反射：边界条件没压住

现象：波场快照在模型左右或底部边缘出现明显的弧形同相轴，像是地下有个强界面。原因：吸收边界或者阻尼层参数设置不足，人工边界产生了虚假反射。解决：加宽吸收层、调整阻尼系数，或者把模型范围扩大，让有效波场在到达边界前已经衰减掉。判断方法很简单，把模型横向加宽一倍再跑，如果那个“界面”位置跟着边界移动，就是边界问题。

4.2 合成记录上出现周期性高频噪声：网格频散

现象：炮集记录上出现一串周期性的高频振荡，尤其在速度变化剧烈的区域。原因：网格间距相对波长太大，有限差分格式产生数值频散，波速随频率变化。解决：减小网格间距，或者提高时间采样率满足 CFL 条件；也可以在主频较高的模型里适当降低震源主频。注意，频散和真实的高频反射在形态上不同，频散通常表现为波前后面的“尾巴”，跟着波前一起走。

4.3 LAS 曲线导入后速度场呈阶梯状：深度采样不匹配

现象：Model builder 里生成的速度模型沿深度方向出现台阶状跳变，正演记录上有对应深度的异常反射。原因：LAS 曲线的采样间隔和模型网格的深度步长不一致，插值时没有做重采样。解决：导入前先把曲线重采样到模型网格的深度间隔，检查深度起止范围是否对齐，单位是否统一。这个坑很常见，尤其是不同来源的测井曲线混用时。

4.4 各向异性模型跑得极慢甚至内存溢出：参数维度没控制好

现象：切换到各向异性介质后计算时间暴涨，或者直接报内存错误。原因：各向异性需要存储刚度张量的多个分量，网格上每个点的参数量成倍增加，如果网格还按各向同性时的密度划分，内存和计算量都会超预期。解决：适当放宽网格间距，或者缩小模型范围分段计算；确认刚度张量参数是否只填了必要的分量，避免冗余存储。

4.5 处理后的剖面恢复不出原模型：偏移速度场用错了

现象：Processing Block 走完偏移，剖面和建模时的地层结构对不上，界面位置偏移或者形态扭曲。原因：偏移用的速度场和正演时的真实速度场不一致，或者偏移算法假设的介质类型和模型不符。解决：先用正演的速度场做偏移，确认能恢复后再逐步替换成处理流程估计的速度场，对比差异定位问题。这一步本质上是在检验速度分析环节的精度。

5. 设计取舍：为什么正演软件的价值在“可控的已知”而不是“逼真的未知”

Tesseral 2D 这类工具最核心的设计取舍，是把复杂度控制在二维、可控、可验证的范围内。它不追求三维全弹性、全各向异性的极致逼真，而是用有限差分在二维剖面上把 P 波、SV 波、梯度、各向异性、吸收这些关键因素分别打开，让你一次只验证一个变量。这种“分而治之”的思路，比一上来就堆全物理模型更实用——因为正演的目的是建立因果关系，不是复刻现实。现实中你永远不知道真实地下是什么，但在正演里你知道，这就是它最大的价值。

从实现角度看，四个功能包的划分也体现了这个取舍：Model builder 负责输入可控，Computation Engine 负责物理可控，Viewer 负责观察可控，Processing Block 负责验证可控。每一环都把“已知”牢牢握在手里，才能在处理流程里定位问题。常见的实现思路里，正演软件最容易走偏的方向是追求视觉上的逼真，波场快照画得漂亮但物理参数经不起推敲；真正有用的正演，是参数可追溯、边界可检查、结果可复现。

我自己的习惯是，每次拿到一个新的地质模型，先跑一个最简单的均匀介质或者单界面模型，确认波场快照、初至时间、反射形态都符合解析解或者手算预期，再逐步加复杂度。这个“从已知到未知”的验证顺序，比直接上复杂模型然后对着异常结果猜原因要省时间得多。正演软件给的是可控的已知，用好它的前提是你自己先守住这个顺序。希望帮到你。

本文还有配套的精品资源，点击获取