[发明专利]基于全波波场VSP/RVSP地震资料的逆时偏移方法

权利要求书

1.一种基于全波波场VSP/RVSP地震资料的逆时偏移方法，其特征在于通过以下过程来实现：

(1)VSP/RVSP地震资料的观测系统变换；

(2)有限差分算法求解波动方程，具体如下：

1)时间-空间域的波动方程，经过差分后，利用有限差分法进行数值求解，其中，任意阶差分精度的全波波动方程如下：

2)对上式进行傅里叶变换，得到频空域的形式：

3)确定波场外推计算的时间步长，即需要满足的稳定性条件：

Dt<h/(2.0*Vmax)

其中，h表示网格剖分后最大的网格间距，Vmax表示速度模型中最大速度值，Dt表示波场外推计算的最大时间步长；

(3)依据速度场分布情况对VSP/RVSP逆时偏移的计算区域进行变网格剖分；

(4)利用高阶有限差分算法对网格化的计算区域进行波动方程求解，实现对震源波场的正向延拓和检波波场的反向延拓；

(5)对延拓后的波场应用成像条件，该成像条件是考虑照明度的影响，根据下式

对震源波场求取其能量分布，以此对成像的相干值进行规则化；

(6)在所有数据处理完后，进行最终的偏移叠加，并利用拉普拉斯滤波或带通滤波方法做叠后去噪，得出最终的VSP/RVSP成像剖面。

2.根据权利要求1所述的基于全波波场VSP/RVSP地震资料的逆时偏移方法，其特征在于所述的过程(1)具体为：在逆时偏移前需要对数据进行道头编辑；基于射线路径可逆原理，进行炮检置换和道集重排，得到适用于逆时偏移的共炮点道集。

3.根据权利要求1所述的基于全波波场VSP/RVSP地震资料的逆时偏移方法，其特征在于所述的过程(3)中根据地下地质构造的形态，偏移速度场的分布情况，合理进行计算区域的网格化，保证网格尺寸满足以下公式，从而波场计算不会产生频散：

h<Vmin/(10.0*Fmax)

其中，h表示网格剖分后最大的网格间距；Vmin表示速度模型中最小速度值；Fmax表示地震子波的最大频率。