[发明专利]一种绕射波与反射波分离成像方法

说明书

技术领域

本发明属于地震资料成像处理领域，具体涉及一种绕射波与反射波分离成像方法，用于绕射波信号提取和小地质异常体或断层清晰成像。

背景技术

由于地震记录中来自小地质异常体的绕射波信号能量较弱，要对其精确成像通常需要在成像前进行波场分离或绕射波提取。目前已有的平面波域解构滤波绕射波与反射波分离成像方法是一种有效的方法，其包含τ-p变换、共p平面波数据分选、平面波域解构滤波、平面波叠前偏移等步骤。然而要使平面波叠前偏移与传统炮域叠前深度偏移的成像效果相当，需要在τ-p变换时分解出大量的平面波数据，这往往又为后面的偏移带来巨大的计算量，甚至使得平面波叠前偏移的效率反而比传统炮域偏移的效率还低。为此，可以将波场分离后的平面波数据通过τ-p反变换重新获得叠前道集，再采用传统炮域叠前偏移实现地震成像。常规二维τ-p正反变换所用信息仅仅是偏移距的值，故当两道的偏移距值相近而炮检位置差异很大时，两道地震数据在τ-p正反变换中将混叠在一起，出现τ-p反变换数据信号失真的问题，影响绕射波与反射波分离成像的效果。

发明内容

本发明针对叠前偏移成像中绕射波信号弱、成像不清晰的问题，提供一种绕射波与反射波分离成像方法，根据地震数据中的炮检坐标信息进行三维时间域τ-p正反变换，将平面波解构滤波后的平面波域绕射波数据完整地转换为叠前道集，进而对绕射波信息进行炮域叠前偏移成像，从而避免二维τ-p正反变换信号失真问题。

本发明是通过以下技术方案实现的：

一种绕射波与反射波分离成像方法，包括以下步骤：

(1)对每一条线内的每一个道集进行三维τ-p正变换得到平面波数据d(p_y，p_x，y，x，τ)；

(2)对由步骤(1)输出的平面波数据d(p_y，p_x，y，x，τ)采用平面波域解构滤波进行处理，得到主要含绕射波能量的平面波数据体d_r(p_y，p_x，y，x，τ)；

(3)从步骤(2)得到的主要含绕射波能量的平面波数据体d_r(p_y，p_x，y，x，τ)中抽取每一个道集对应的三维τ-p域数据，并从原道集中获取每个道集的炮检点坐标，然后对所述三维τ-p域数据做三维τ-p反变换得到波场分离后的叠前道集；

(4)对由步骤(3)输出的波场分离后的叠前道集进行炮域叠前深度偏移获得绕射波分离成像结果(即分离掉了反射波而只包括绕射波的成像结果)。

所述步骤(1)具体如下：

(11)采用下面的三维τ-p正变换公式对每一条线内的每一个道集进行三维τ-p正变换得到τ-p域数据ψ：

ψ(py,px,τ)=Σiφi(τ+pxhix+pyhiy)]]>

其中，φ_i表示第i道数据，h_ix、h_iy是第i道数据的两个方向的偏移距分量，此两个偏移距分量是由第i道数据的炮检点坐标计算得到的，p_x、p_y表示两个方向的射线参数；