一种消除逆时偏移低频假像的方法

摘要

本发明在于提供一种在进行滤波的同时，尽可能地保持地震波动力学特征的消除逆时偏移低频假像的方法。采用深度域层速度模型和正演炮集数据，首先计算补偿后的震源子波，在深度域进行单炮波动方程正演，存储所有时刻单炮源波场；然后用单炮记录作为初始边界，沿时间轴逆向进行波动方程正演得到接收波场，在波场逆推时读取相应时刻单炮源波场，计算补偿后的单炮逆时偏移成像，得到低频假像消除后的单炮逆时偏移成像。在构造成像和波形相对保持、沿层地震属性空变稳定性方面，本发明提出的逆时偏移低频假像的消除方法要好于Laplace算子直接滤波方法。

主权项

1.一种消除逆时偏移低频假像的方法，特点是采用以下步骤实现：1）采集地震勘探数据，处理得到深度域层速度模型和正演炮集数据；2）利用公式计算补偿后的震源子波；其中S°(t)为补偿后的震源子波，S(t)为震源子波，*为褶积，u(t)为阶跃函数，t为时间；3）利用补偿后的震源子波作为震源，利用下式在深度域进行单炮波动方程正演：并存储所有时刻补偿后的单炮源波场；其中：v为地震波速度，P_s°(X,t)为补偿后的单炮源波场，Δ为Laplace算子，δ(X-X_s)为单位脉冲函数，X=(x,z)为波场所在二维空间位置坐标，X_s=(x_s,z_s)为震源所在二维空间位置坐标，x和x_s分别为波场和震源横坐标，z和z_s分别为波场和震源纵坐标；4）利用单炮地震记录作为初始边界条件，用下式从最大记录时刻开始，沿时间轴逆向进行波动方程正演，得到接收波场：$\frac{1}{v^2}\frac{{\partial }^2{P}_g(X,t)}{\partial t^2}-\Delta P_g(X,t)=G(X_g,t);$其中：P_g(X,t)为接收波场，G(X_g,t)为单炮地震记录，X _g=(x_g,z_g)为检波器所在二维空间位置坐标，x_g为检波器横坐标，z_g为检波器纵坐标；5）在接收波场逆推时，读取相应时刻补偿后的单炮源波场，利用以下公式计算补偿后的单炮逆时偏移成像：I°(X)＝∫P_s°(X,t)P_g(X,t)dt；6）利用公式I_SF(X)=S□F{I°(X)}进行逆时偏移低频假像消除，得到低频假像消除后的单炮逆时偏移成像；其中：S为满足二阶导数连续的平滑算子以保证成像波形的连续性（所述的平滑算子采用三次样条函数插值平滑算子），F=v²Δ为滤波算子；7）将所有单炮逆时偏移成像进行叠加，得到消除低频假像的叠加逆时偏移成像。