[发明专利]基于均炮源波场照明的逆时偏移空间振幅补偿方法

说明书

技术领域

本发明涉及勘探和油藏地球物理技术，具体是一种基于均炮源波场照明的逆时偏移空间振幅补偿方法。

背景技术

在勘探和油藏地球物理中，地震偏移成像是其中最关键的一步，成像的好坏直接影响到后续的地震和地质解释，从而影响油气储层的识别。地震偏移成像可以分为几大类：（1）基于射线理论的Kirchhoff积分法偏移；（2）基于单程波动理论的F-K偏移和有限差分偏移；（3）基于双程波动理论的逆时偏移。（1）和（2）属于常规地震偏移，对横向速度剧烈变化、高陡倾角等复杂地质构造，成像能力有限。而逆时偏移是基于双程波动理论，理论上它能解决常规地震偏移所面临的问题，当然它也会带来很多新的问题。逆时偏移是利用源波场和接收波场零延迟互相关来进行成像（Claerbout,1971;Biondi和Shan,2002）。对于简单地质构造，它能很好地成像；对于逆掩断层等复杂地质构造，由于逆掩断层下盘地震波场照明较差，从而导致成像振幅很弱。为了对这些弱振幅进行补偿，一些学者进行了研究和尝试。Kaelin和Guitton（2006）通过对每炮互相关成像除以源波场或接收波场照明来改进成像条件，从而进行振幅补偿。Chattopadhyay和McMechan（2008）对几种不同成像条件的振幅特性进行分析，并与理论反射系数进行了比较，通过源波场照明补偿的互相关成像有正确的角度特性、尺度因子、符号和单位。一般采用基于单炮源波场照明的振幅补偿方法来进行弱振幅补偿，这种方法比较容易实现，也不会增加大的计算和存储量，但是它会破坏地震波的动力学特征，从而影响后续的油藏精细描述。

发明内容

本发明目地是克服基于单炮源波场照明的逆时偏移振幅补偿方法在地震波动力学特征保持方面存在的问题，提供一种在进行振幅补偿的同时，尽可能地保持地震波动力学特征的基于均炮源波场照明的逆时偏移空间振幅补偿方法。

本发明通过如下步骤实现：

1）采集地震勘探数据，处理得到深度域层速度模型和正演炮集数据；

2）利用下式在深度域进行单炮波动方程正演：

1v2∂2Ps(X,t)∂t2-ΔPs(X,t)=δ(X-Xs)S(t),]]>

并存储所有时刻的单炮源波场；同时利用下式计算单炮源波场照明：

I_Source(X)=∫P_s(X,t)P_s(X,t)dt；

其中：v为地震波速度，P_s(X,t)为单炮源波场，Δ为Laplace算子，δ(X-X_s)为单位脉冲函数，S(t)为震源子波，X=(x,z)为波场所在二维空间位置坐标，X_s=(x_s,z_s)为震源所在二维空间位置坐标，x和x_s分别为波场和震源横坐标，z和z_s分别为波场和震源纵坐标，t为时间；

3）利用单炮地震记录作为初始边界条件，用下式从最大记录时刻开始，沿时间轴逆向进行波动方程正演，得到接收波场：