基于地表数据分离的层间多次波压制.pdfVIP

石 油 物 探 第53卷 形态往往比自由表面多次波更加复杂E 。此外，由 方法是基于数据矩阵的，因此可适用于二维和三维 于在地层界面之间经历了多次反射，层间多次波与 地震数据。 一 次波的速度差异较小，基于多次波周期性与可分 在海洋地震勘探中，海水表面的下行反射作用 离性的滤波方法很难取得压制效果。因此，波动方 会产生强烈的自由表面多次波，因此，地震数据可 程预测相减法成为衰减层间多次波的首选。 表示为一次波P。(含有层间多次波)与自由表面多 反馈迭代法是一种基于波动方程的多次波预 次波M。之和，即 测相减方法，它以描述地震数据正演过程的反馈模 P—Po+ (1) 型为理论基础，利用地震波在自由表面、层界面间 Verschuur等 提出了衰减自由表面多次波 的传播规律来预测多次波 ]，能够适应复杂的地下 的SRME方法，该方法通过引入表面算子A，得到 结构。针对层间多次波，反馈迭代法一般采用层剥 了完全数据驱动的自由表面多次波预测公式： 离的方式，按照反射界面由浅至深的顺序，递归消 一 PoAP (2) 除各界面的下行反射效应。实际地下介质中，能够 产生明显层间多次波的强反射界面数量有限，因此 其中，A S_。R ，s表示震源子波，R为海水表面 可以在叠加剖面上人工选取基准面，把强反射界面 的下行反射系数矩阵，一般R— I。 划分到不同的反射层，从而以反射层为基本物理单 图1描述了公式(2)所示的自由表面多次波预 元预测层间多次波，然后通过自适应相减实现层间 测过程，由炮点 激发，检波点 接收的自由表面多 多次波的压制。 次波，可通过红线所示的共炮点道集P与蓝线所示 目前，反馈迭代层间多次波压制方法主要分为 的共检波点道集 ，沿地表进行一致性褶积后叠加 两种：基于延拓数据的共聚焦点(CFP)方法 与 得到，表面算子A的作用相当于对预测得到的自由 表面多次波进行震源子波反褶积。实际应用中子波 基于地表数据分离的Jakubowicz方法_1。11]。CFP 方法需要将地表波场反向延拓至地下界面，构建炮 未知，预测得到的多次波与真实地震记录存在振幅 点在地表，检波点在地下的CFP道集，以及炮点和 或相位上的差异，因此，必须从地震记录中自适应减 检波点分别位于地下的网格点道集，计算量较大， 去，以达到衰减自由表面多次波的目的。 对地下速度模型有一定程度的依赖。相比之下，基 于地表数据分离的方法直接利用地表观测数据中 分离出的三部分波场数据，进行地表一致性相关与 褶积运算来预测层间多次波，从而避免了地表数据 的反向延拓，具有计算效率高，不依赖地下速度的 优势。 我们首先阐述基于地表数据分离的层间多次 波压制方法原理，并通过模拟算例展示该方法压制 层间多次波的具体步骤；然后针对地下介质中存在 图1 自由表面多次波预测示意图解 多个强反射界面的地震数据，重点论述逐层压制层 间多次波的递归实现过程；最后用实际地震资料验 按照界面由浅至深的顺序逐步消除多次波的 证了方