Question

三维地震勘探资料处理

Answer 1

二维处理所包括的主要内容，在三维地震勘探资料处理的流程中一般也是需要的。专门用于三维处理的三维偏移以及成果显示是流程的重要环节，其中包括三维速度分析、三维剩余静校正、三维叠加、宽线处理、三维偏移等方面。下面仅对三维处理中最重要的速度分析和偏移归位加以讨论。

（一）三维速度分析

1.叠加速度的方向性变化

对于平倾斜界面而言，叠加速度随方位的变化可以方便地导出。众所周知，此时共中心点时距曲线方程仍可写成双曲线方程形式，只须将所用速度改为等效速度v_φ。因此，速度分析求出的叠加速度v_st就是等效速度：

反射波地震勘探原理和资料解释

式中，v_a为均方根速度v_rms或称沿界面走向的速度，φ为反射界面视倾角。由视倾角与真倾角ψ及方位角α^*的关系：

sinφ＝sinψcosα^*

可以得到：

反射波地震勘探原理和资料解释

这是一个椭圆，主轴沿界面倾向和走向。

当反射界面不是一个平界面时，在一阶近似的意义上叠加速度的方向变化可用此椭圆关系表示。

若覆盖地层比较单纯，则沿界面走向为最小叠加速度，沿界面倾向为最大叠加速度。当上覆地层复杂时，v_a不再能认为是均方根速度，叠加速度方位变化椭圆的主轴也不一定沿界面倾向和走向，整个上覆地层对椭圆形状和主轴都有影响。此时，式（7-5-2）需要化为

反射波地震勘探原理和资料解释

式中，e为椭圆的离心率：

反射波地震勘探原理和资料解释

v_max、v_min分别是最大和最小叠加速度，α^*是椭圆主轴方位角，

是求叠加速度方向的方位角。此式中v_max、v_min和α^*均为未知数。因此，欲确定叠加速度方位变化椭圆，必须至少求出三个方位处的叠加速度。

2.三维叠加速度的求取

由前述第五章可知，在二维地震勘探中叠加速度取决于地层速度和界面倾角，共反射点时距曲线符合双曲线规律，故通过速度谱分析就可以求得叠加速度。在三维地层模型中，反射界面的倾角的走向可以是任意的。这样，反射同相轴的时差就不仅受地层速度的倾角的控制，而且还受炮点-检波点方位与地层走向方位之间的相对方位角的控制。因此，对三维面积观测的地震资料就不能用前述的求取二维地震资料的方法来求取二维叠加速度，而必须进行三维速度分析以求取比较可靠的三维叠加速度。目前，求取三维叠加速度的方法很多，这里介绍一个简单易行的扇形分析技术。

考虑一个CDP选排中像蜘蛛网似的炮检矢量分布，图7-5-10。将它划分成若干个扇形。划分的原则是不能太小，太小会增加工作量，且造成一扇形内道数太少，影响结果精度；但也不能太大，太大反映叠加速度方位变化就不灵敏，误差也比较大。另外，尽量使每个扇形内保持有数量相近的道数。划分了扇形之后，将一个扇形内的所有道组成一个虚二维共反射点道集，用标准的二维速度分析方法计算速度，结果置于扇形的中心方位。然后用最小平方拟合技术求出叠加速度方位椭圆（求最大叠加速度、最小叠加速度和椭圆主轴方位角）。为了完成椭圆的计算至少需要划分三对扇形。用虚二维CDP道集计算叠加速度时可适当考虑加权问题。

图7-5-10 扇形分析示意图

在速度横向变化不大的地区，常常忽略方位对叠加速度的影响。为了保证采用统一速度校正对叠加结果的影响为最小，在数据采集的设计阶段对最大炮检距要加以限制。实际中一般应使最大炮检距S_max小于或等于目的层深度z。另外，处理中还需采用平滑的方法。因为对各方位叠加速度取平均是求取统一速度的常用方法，但这样求出的速度总是不准确的，势必造成速度曲面（即同一界面叠加速度空间变化曲面）出现抖动现象。对速度曲面进行二维平滑处理有助于求出较准确的叠加速度，消除抖动，当然即使用这种准确的叠加速度对道集进行校正叠加，仍会引起共反射点的分散，使时间剖面上同相轴的连续性、光滑性受到影响。这一问题可在自动剩余静校正中解决或直接对同相轴进行光滑处理。

（二）三维偏移归位处理

1.三维偏移归位处理概述

三维偏移归位是三维资料处理中的核心部分，它集中体现了三维勘探的优点。前述的二维偏移处理，只能使存在于二维（x，z）剖面内的反射同相轴归位，绕射波、回转波等收敛。实际中地下地质体是三维的，即使在地面沿一条测线观测时，只要满足斯奈尔定律，在包括测线的各个射线平面内的反射都会反映到该测线剖面中来，即水平叠加剖面上必然会混杂着由测线位置的法平面之外传来的侧面波。因此，只进行所谓“平面型”的二维偏移处理不可能使这些侧面波归位和收敛，这些侧面波会构成假象掩盖真实界面。故只有进行“立体型”的三维偏移处理才可能使侧面波归位，去掉假象。目前常用的三维偏移方法主要有全三维偏移和两步法偏移两种，其原理见图7-5-11。图中为地下P点的三维共深度点或中心点道集所形成的旋转双曲面，全三维偏移就是把双曲面上所有的点一次归位到P点上去。两步法偏移是先对某一方向做好第一步的二维偏移，再在其垂直方向做好第二步的二维偏移从而得到最终成果。例如，先在y方向把A、B、C、D、E、F、G……各点上资料偏移归位到D点，见图7-5-11，并以此类推；再在x方向把第一步偏移归位好的D_-1、D₀、D₁、D₂、D₃……各点的资料偏移归位到P点。由于三维勘探得到的数据量十分大，全三维偏移（或称一步法偏移）工作量相当大，因此目前多用某些效率较高的近似方法，就是两步法和分裂法。

图7-5-11 全三维偏移和两步法偏移原理图

2.两步法和分裂法偏移

（1）两步法偏移。两步法的概念建立在认为一个三维问题可以完全分离成若干个二维问题解决的基础上。在具体实现时，一个三维偏移用互相垂直方向上依次进行的两次二维偏移代替，先沿x方向对所有的三维资料进行一次二维偏移，然后将偏移后的各道按y方向进行一次二维偏移就完成了近似的三维偏移工作。

其优点是极大地减少了计算工作量，而且现行的所有二维偏移程序在用于实现三维偏移时均可不作任何修改（仅在其中加上若干组织工作程序）。

（2）分裂法偏移。比两步法更精确的方法是分裂法。它是将一个三维偏移方程分裂为两个二维偏移方程，即将三维偏移化为两个二维偏移。但这两个二维偏移方程互相耦合，不是独立的，每延拓一次要交替使用两个方程。因此，分裂法与两步法不同之处在于它要交替地计算x方向和y方向上的延拓剖面，相当于一种递归运算，是一种递归偏移。

由此可见，虽然分裂法比两步法精确，但工作量大一些。图7-5-12给出了三维偏移和二维偏移的对比图。图7-5-12（b）是二维偏移剖面，可以在2.0s左右看到有明显的侧面波同相轴和反射同相轴斜交；图7-5-12（a）是经过三维偏移处理的剖面，侧面波被去掉，记录形态更清楚。

图7-5-12 三维偏移（a）和二维偏移（b）剖面比较

3.一步法与两步法的区别

用绕射扫描求和思想很容易说明一步法和两步法的区别。见图7-5-13（a），三维偏移要将绕射双曲面上任一点E（x_e，y_e，t）处的采样值归位于双曲面顶点D（x_d，y_d，t₀）处。一步法偏移使用直接归位法，即直接将E点处的样值取至D点与其他归位样值叠加。两步法先将E点的样值归位于沿x方向切割的双曲面顶点A（x_α，y_α，t₀），待此方向二维偏移做完后，再将这些双曲线顶点处的归位值沿y方向归位于双曲线顶点D。从形式上看二者毫无区别，实际上在速度纵向变化时二者结果是不同的。

图7-5-13 一步法和两步法的区别

做一步法偏移时，计算绕射双曲面所用的速度是D点处的速度v（t₀）［图7-5-13（b）］。两步法沿x方向做第一次二维偏移时，计算绕射双曲线所用的速度也应当是v（t₀），但实际使用的却是A点处速度v（T₀）图7-5-13（b）（因t₀在何处此时并不知道）。当v（T₀）＝v（t₀）时，当然没有问题。只要存在简单的纵向速度变化，就必然出现误差。此时并没有按切割绕射双曲面得到的双曲线进行归位，而是按另一条用v（T₀）算出的双曲线（即图7-5-12（b）中的虚线）进行归位。归位到A点的能量并不完全是D点绕射波的能量，结果有问题。误差的大小取决于速度、速度梯度、旅行时、倾角、方位角等因素。不过，虽然两步法有误差，但通常情况下误差不很大，且它的效率十分高，故仍经常使用。