一种解决逆时偏移波形不连续性的方法

摘要

本发明是能够提高油藏精细描述的有效解决逆时偏移波形不连续性的方法。采用深度域层速度模型和正演炮集数据，首先选取震源子波，在深度域进行单炮波动方程正演，存储所有时刻单炮源波场；然后用单炮地震记录作为初始边界条件，从最大记录时刻开始沿时间轴逆向进行波动方程正演得到接收波场，利用本发明提出的平滑成像条件，得到波形连续后的单炮逆时偏移成像，将所有单炮逆时偏移成像进行叠加得到逆时偏移成像。本发明在波形连续性、沿层地震瞬时频率和相位空变稳定性上，平滑成像条件要好于常规成像条件。

主权项

一种解决逆时偏移波形不连续性的方法，特点是通过以下步骤实现：1）采集地震数据，处理得到深度域层速度模型和正演炮集数据；2）利用下式在深度域进行单炮波动方程正演： 1 v 2 ∂ 2 P s ( X , t ) ∂ t 2 - Δ P s ( X , t ) = δ ( X - X s ) S ( t ) , 并存储所有时刻的单炮源波场；其中：v为地震波速度，Ps(X,t)为单炮源波场，Δ为Laplace算子，δ(X‑Xs)为单位脉冲函数，S(t)为震源子波，t为时间；X=(x,z)为波场所在二维空间位置坐标，Xs=(xs,zs)为震源所在二维空间位置坐标，x和xs分别为波场和震源横坐标，z和zs分别为波场和震源纵坐标；3）用单炮地震记录作为初始边界条件，利用下式从最大记录时刻开始，沿时间轴逆向进行波动方程正演，得到接收波场： 1 v 2 ∂ 2 P g ( X , t ) ∂ t 2 - Δ P g ( X , t ) = G ( X g , t ) ; 其中：Pg(X,t)为接收波场，G(Xg,t)为单炮地震记录，X g=(xg,zg)为检波器所在二维空间位置坐标，xg为检波器横坐标，zg为检波器纵坐标；4）在接收波场逆推时，读取相同时刻的单炮源波场，利用下式得到波形连续后的单炮逆时偏移成像：IS(X)=S{∫Ps(X,t)Pg(X,t)dt}；其中：S为满足二阶导数连续的平滑算子；5）将所有单炮逆时偏移成像进行叠加，得到波形连续后的叠加逆时偏移成像。