3. 地震数据预处理与反演准备

叠后地震反演的质量，在很大程度上并不取决于反演算法本身，而取决于进入反演之前的数据准备是否扎实。地震数据加载、极性判断、井震标定、子波提取与质量控制这五项工作，构成了反演准备的完整链条：加载保证数据"进得来、对得上"，极性判断保证反射系数与地震响应"符号一致"，井震标定建立时间域地震与深度域测井之间的对应关系，子波提取给出反演所需的算子，质量控制则贯穿始终，决定这些中间成果能否被信任。本节按这一逻辑展开，重点讲清每一步"做什么、怎么做、怎么判断做对了"。

一、地震数据加载：先解决坐标、采样与道头

加载不是简单地把 SEG-Y 读进来。反演对输入的要求比常规解释更严格，需要重点核对三件事。其一是采样率与时间范围：反演通常要求地震数据与测井曲线在相同的时间采样间隔上工作，常见为 1 ms 或 2 ms，若地震为 2 ms 而期望 1 ms，需要明确是重采样还是保持原采样，避免后续子波与反射系数序列长度不匹配。其二是道头与坐标：检查 CDP 号、Inline/Crossline、炮检距、面元坐标是否完整，井位坐标能否正确投影到地震网格上，这是井震标定能"找对道"的前提。其三是振幅与单位：确认数据是相对振幅还是经过增益处理，是否存在异常道、空道、极性反转道。

一个实用的加载检查清单如下（示例参数，实际以工区为准）：

二、极性判断：SEG 标准与"同相轴符号"

极性是反演准备中最容易被忽视、后果却最严重的一环。业界约定：SEG 正极性表示地下波阻抗增大（反射系数为正）对应地震波形的正瓣（波峰）；SEG 负极性则相反，波阻抗增大对应波谷。若地震数据与测井反射系数序列极性相反，标定相关系数会明显偏低，强行反演会得到与地质认识相悖的阻抗趋势。

判断方法通常有两条：一是用已知的强反射界面（如海底、区域不整合、厚层灰岩顶）对照，看波阻抗增大方向对应波峰还是波谷；二是用合成记录与井旁道做正、负两种极性的相关对比，相关系数高者即为数据真实极性。工程上建议在工区内选取 3～5 口分布不同的井重复验证，避免单井偶然性。极性一旦确定，应在整个反演流程中固定，不得中途翻转。

三、井震标定：建立时间—深度—反射系数的对应

井震标定是把深度域的声波（AC）与密度（DEN）曲线转换为时间域的波阻抗与反射系数，再与井旁地震道对比，通过调整时深关系使合成记录与地震道最佳吻合。关键步骤包括：曲线环境校正与去异常值、由 AC 与 DEN 计算波阻抗 Z = ρ·V、由相邻采样点计算反射系数 R = (Z₂−Z₁)/(Z₂+Z₁)、给定初始子波做合成记录、用相关分析微调时深漂移。

标定的质量判据不是"看起来像"，而是可量化的：井旁道与合成记录的相关系数一般要求达到 0.7 以上，目标层段内主要同相轴应逐一对应，时深漂移量应平滑、无突变。若某一段相关系数骤降，优先怀疑曲线异常、极性错误或该段地震资料信噪比低，而不是反复拉伸时深表去"凑"。

四、子波提取：相位、频带与长度

子波是反演算子，其相位误差会直接转化为阻抗的漂移。提取方式分两类：统计性子波由井旁地震道自相关估算，适用于井少或井曲线质量一般的地区；确定性子波由标定后的井反射系数与井旁道联合求取，精度更高但依赖标定质量。无论哪种，都应关注三个参数：相位（优先零相位，若数据为混合相位需在反演中保留）、频带（应与地震有效频带一致，常见主频 25～40 Hz、频带约 8～60 Hz，示例值）、长度（一般 80～120 ms，过长会引入噪声，过短则旁瓣截断）。

判断子波是否可用，可看其自相关旁瓣是否对称衰减、振幅谱是否与地震振幅谱形态接近、多井提取的子波是否彼此相似。若不同井提取的子波相位差异很大，说明标定或极性存在问题，应先回到前两步排查，而不是任选一个子波强行反演。

五、质量控制：把"能反演"变成"敢反演"

质控应贯穿全流程并形成可核查的记录。核心指标包括：井旁道与合成记录相关系数（示例要求 ≥0.7）、时深漂移的平滑性、子波在多井间的一致性、地震数据在目标层段的信噪比与频带宽度。建议在正式反演前，用 2～3 口参与标定的井做"盲测"式检查——即用其余井建立的子波与趋势去预测该井阻抗，观察吻合程度，这比单纯看相关系数更能暴露系统性问题。

完成上述五项工作后，应形成一套可追溯的反演输入包：统一采样率的地震数据、明确的极性标识、各井的时深关系与反射系数序列、选定的子波及其参数、以及质控指标汇总表。这套输入包的质量，基本决定了后续约束稀疏脉冲反演或地质统计学反演结果的上限。