3. 薄互层地震正演模拟方法

要理解薄互层为什么会让地震波形“变样”，最直接的办法不是去猜，而是用正演把已知的地质模型“合成”成地震记录，再观察波形随参数如何变化。正演是薄互层解释的“标尺”：只有先知道不同厚度、速度、密度组合会生成什么样的波形，才能在实际资料中反过来判断厚度和岩性。本节围绕三类常用正演手段展开，并重点说明薄互层参数对波形的影响规律。

三类正演手段的定位与适用场景

褶积模型是最基础、最快的方法。它把地下视为水平层状介质，先由波阻抗差异算出反射系数序列，再与地震子波做褶积，得到合成道。其核心假设是：地震波垂直入射、层间无多次波、子波不随深度变化。对于薄互层，褶积模型能清晰展示“反射系数对”与子波的干涉关系，是理解调谐效应最直观的工具。但它无法处理层间多次波和透射损失，薄层数量多、阻抗差异大时误差明显。

反射率法在频率域严格求解层状介质的弹性波传播，能同时包含一次波、层间多次波和透射效应，适合薄互层这种“层多、层薄、多次波不可忽略”的情形。它计算量大于褶积模型，但精度高，常用于检验褶积模型结果的可靠性。使用时需给定每层的厚度、纵波速度、横波速度和密度，并指定子波与入射角范围。

有限差分法直接在时间—空间域求解波动方程，能模拟复杂构造、横向变速和非水平界面，是验证薄互层横向变化、绕射和层间干涉的有力手段。代价是计算量大、需要设置合适网格与边界条件，且频散和稳定性对参数敏感。实际工作中常形成“褶积快速扫描—反射率法校验—有限差分验证复杂情形”的组合策略。

薄互层参数对波形的影响

薄互层波形是顶底反射相互干涉的结果，其形态主要由三个参数控制。第一是厚度：当单层厚度小于调谐厚度（约四分之一主波长）时，顶底反射无法分开，合成波表现为一个振幅增强的复合波；厚度继续减小，复合波振幅先增后减，视厚度却持续变大。第二是速度与密度：二者共同决定波阻抗，波阻抗差越大，反射系数越强，复合波振幅越高；若相邻薄层阻抗相近，反射相互抵消，波形变弱甚至“消失”。第三是薄层组合方式：多层薄互层叠加时，各层反射的极性和时移共同决定复合波是增强还是削弱，这正是薄互层解释容易产生“假厚层”或“假弱反射”的原因。

下面给出一个可直接复现的褶积模型示例（参数为教学假设，非实测）。设主频 35 Hz 的雷克子波，采样间隔 1 ms，取三层薄互层，厚度分别设为 4 m、6 m、8 m、10 m 四组，速度与密度保持不变，观察复合波振幅与视厚度的变化。

# 褶积模型示例（教学假设参数）
dt = 0.001          # 采样间隔 1 ms
f0 = 35             # 子波主频 Hz
v  = [2800, 3200, 2900]   # 各层纵波速度 m/s
rho= [2.35, 2.50, 2.40]   # 各层密度 g/cm3
h_list = [4, 6, 8, 10]    # 薄层厚度 m（四组对比）

# 1) 由波阻抗计算反射系数
Z = [vi*ri for vi, ri in zip(v, rho)]
rc = [(Z[i+1]-Z[i])/(Z[i+1]+Z[i]) for i in range(len(Z)-1)]

# 2) 按厚度换算双程旅行时，放置反射系数
# 3) 生成雷克子波并与反射系数序列褶积
# 4) 读取复合波峰值振幅与视厚度（波峰-波谷时间差）

按上述流程计算，可观察到：厚度从 10 m 减到 6 m 时，复合波峰值振幅逐渐增大，在约 6～7 m 附近达到最大（接近调谐厚度）；厚度继续减到 4 m，振幅开始回落，但波峰与波谷的时间差反而继续增大，形成“视厚度变大”的假象。若把速度差调小（如把中间层速度改为 3000 m/s），整体振幅明显下降，说明波阻抗差是振幅的主控因素。实际解释中，应先用井旁道标定“振幅—厚度”和“视厚度—厚度”两条关系曲线，再用于横向外推。

选择正演方法时，可遵循一条实用原则：先用褶积模型快速扫描厚度与阻抗组合，建立波形变化规律；再用反射率法校验多次波和透射损失对薄互层振幅的影响；对横向变化明显或构造复杂的靶区，用有限差分做最终验证。三种方法互相印证，才能把“波形为什么这样”讲清楚，也才能把正演结论可靠地用到厚度识别中。