[发明专利]一种自适应抗频散的分块时空双变逆时偏移方法

权利要求书

1.一种自适应抗频散的分块时空双变逆时偏移方法，其特征在于，包括以下步骤：

1)对用于采集地面地震的观测系统覆盖的计算区域做变网格划分后得到空间网格；

2)采用变网格方法在对空间网格做可变处理的基础上，并考虑时间步长的变化；同时，利用波场滤波和三次样条平滑解决变网格方法在大采样时间下的稳定性问题；

3)采用Lanczos滤波方法处理在变网格过渡区域中边点的波场值；

4)对观测系统覆盖的计算区域利用分块时空双变算法，根据逆时偏移的流程进行震源波场的正向外推并存储震源波场；

5)对观测系统覆盖的计算区域利用相同的分块时空双变算法，进行检波波场的反向外推并存储检波波场；

6)应用互相关成像条件，对外推后的震源波场和检波波场进行相关成像，并对成像结果进行常规叠加和拉普拉斯滤波等叠后处理得到最终的逆时偏移结果。

2.如权利要求1所述的一种自适应抗频散的分块时空双变逆时偏移方法，其特征在于，在上述步骤2)中，采用式(1)中的各向同性非均匀介质的2D声波一阶速度—应力方程进行变网格下的波场计算：

ρ∂vx∂t=∂p∂xρ∂vz∂t=∂p∂z∂p∂t=ρVp2(∂vx∂x+∂vz∂x)---(1)]]>

其中，ν_x,ν_z分别表示横向和纵向的质点速度；p表示应力向量；ρ表示正演所用模型的密度；V_p表示正演所用模型的速度；t表示时间；x和z分别表示空间的横向和纵向坐标。