3. 合成地震记录的基本原理与流程

合成地震记录是连接地质模型与地震响应的桥梁。它的核心思想可以用一句话概括：把地下界面的反射系数序列与地震子波做卷积，得到地震道。这一节围绕这个卷积模型展开，讲清楚它为什么成立、需要哪些输入、完整流程如何走通，以及模型成立所依赖的假设条件——这些假设恰恰是实际工作中误差的主要来源。

一、卷积模型的数学表达与物理含义

在一维水平层状介质、垂直入射（自激自收）的简化条件下，地震道 s(t) 可以写成反射系数序列 r(t) 与子波 w(t) 的卷积，再叠加噪声 n(t)：

s(t) = w(t) * r(t) + n(t)
     = Σ w(τ)·r(t-τ) + n(t)

离散形式：s[i] = Σ_k w[k]·r[i-k] + n[i]

物理含义很直观：每一个反射界面产生一个幅度正比于反射系数的子波，整道记录就是所有界面产生的子波在时间轴上叠加的结果。反射系数序列由波阻抗差决定：

R = (ρ₂V₂ - ρ₁V₁) / (ρ₂V₂ + ρ₁V₁)

其中 ρ 为密度，V 为纵波速度，ρV 即波阻抗。界面波阻抗差越大，反射系数绝对值越大；波阻抗由小变大时 R 为正，反之为负。这个正负号约定在后续与子波极性匹配时必须保持一致，否则合成记录与实测道会出现整体极性反转。

二、从反射系数序列到合成记录的完整流程

实际制作一条合成记录，通常按以下链条推进，每一步的输出都是下一步的输入：

1. 建立波阻抗模型：从测井曲线（声波时差 AC、密度 DEN）出发，计算纵波速度与波阻抗曲线。若密度曲线缺失，可用 Gardner 公式 ρ ≈ 0.23·V^0.25（英制单位经验式，需注意单位换算）近似补齐。
2. 时深转换：测井曲线是深度域的，而地震记录是时间域的。用声波时差积分得到单程旅行时 t = ∫(1/V)dz，把波阻抗曲线由深度域采样到时间域，采样间隔通常取 1 ms 或 2 ms。
3. 计算反射系数序列：对时间域波阻抗逐点求相邻样点的反射系数，得到与地震采样率一致的 r(t)。
4. 提取或给定子波：可从井旁地震道提取（统计性子波），也可用理论子波（如 Ricker 子波、带通子波）。子波的主频、相位、长度直接决定合成记录的形态。
5. 卷积合成：将 r(t) 与 w(t) 卷积，得到初始合成道。
6. 标定与匹配：把合成道与井旁实际地震道对比，通过整体时移、拉伸压缩（drift correction）和相位旋转，使主要反射同相轴对齐，从而建立时深关系、标定层位。

三、卷积模型的假设条件及其现实偏差

上面的公式看似简单，但它建立在一组明确的假设之上。理解这些假设，才能判断合成记录在什么情况下可信、在什么情况下会失真。

• 水平层状介质假设：模型假定反射界面是水平的、横向不变的。实际地层存在倾角、断层、尖灭和横向变速，这时一维卷积无法表达绕射、侧向反射和横向振幅变化，合成记录只能近似对应井旁一小段范围。
• 垂直入射假设：公式中的反射系数是垂直入射（零炮检距）反射系数。实际地震数据有炮检距，存在 AVO 效应，反射振幅随入射角变化，用零炮检距子波合成会与远道不匹配。
• 子波不变假设：模型认为子波在传播过程中波形和相位不变。实际地层对高频吸收衰减、速度频散会使子波随深度变宽、主频降低，深层合成记录若仍用浅层子波，同相轴会偏"胖"。
• 线性叠加假设：要求反射系数较小、层间多次反射可忽略。薄层调谐、强波阻抗界面产生的层间多次波会破坏线性关系，使合成振幅偏离实际。
• 噪声可加假设：把噪声当作与信号独立相加，便于分析，但实际噪声（如面波、多次波）往往与信号相关，不能简单相加处理。

四、一个可操作的示例参数

假设某井目的层段深度 2000～2100 m，声波时差由 250 μs/m 渐变到 220 μs/m，密度约 2.4 g/cm³。取 V=1/AC，则速度约由 4000 m/s 增至 4545 m/s，波阻抗由约 9.6×10⁶ 增至 10.9×10⁶ kg/(m²·s)，界面反射系数约为 (10.9−9.6)/(10.9+9.6)≈0.063，属于中等强度反射。选用主频 35 Hz 的零相位 Ricker 子波、采样间隔 1 ms、子波长度 120 ms 做卷积，可得到一组正极性、波峰对应波阻抗增大界面的合成道。随后与井旁道对比，若整体时移约 +8 ms，则对合成道做整体时移校正；若浅深段时移量不一致，则分段拉伸压缩，并检查校正量是否超过半个子波周期——超过则说明速度模型存在系统误差，应回到时深转换环节复核。