[发明专利]针对VTI介质逆时偏移方法、设备及介质

权利要求书

1.一种针对VTI介质逆时偏移方法，其特征在于，包括：

步骤S1，获取炮集数据；

步骤S2，在所述炮集数据中提取震源坐标，并生成震源子波；

步骤S3，根据正演方案，将震源子波加载至波场中，得到震源波场，通过qP波波动方程对所述震源波场沿时间正向延拓，共进行Nc次延拓，Nc为t₀至t_Nc时间段内设定的延拓次数，并在延拓过程中选择若干部分波场进行存储；

在所述qP波波动方程中引入稳定性因子σ；

利用稳定性因子公式实现对横波沿着对称轴方向的传播速度V_sz的计算，其中，通过数值模拟，所述稳定性因子σ为0.75，ε和δ为Thomsen参数，V_pz为纵波沿着对称轴方向的传播速度；

所述步骤S3中，根据正演方案，通过qP波波动方程对所述震源子波沿时间正向延拓，具体为：

根据所述纵波和横波沿着对称轴方向的传播速度以及Thomsen参数作为基于qP波波动方程的VTI介质逆时偏移成像方法的模型参数；

将地震子波作为激发震源生成地震波场，利用有限差分数值解法，对所述震源波场沿时间正向Nc次延拓，Nc为t₀至t_Nc时间段内设定的延拓次数；

其中，qP波波动方程如以下公式所示，

式中，V_pz和V_sz分别为纵波和横波沿着对称轴方向的传播速度，ε和δ为Thomsen参数，V_px为纵波沿着垂直于对称轴方向的传播速度，V_pn为动校正速度，p为声压，q为辅助变量；

步骤S4，根据所述正演方案，通过所述qP波波动方程对所述若干部分波场沿时间逆向延拓，共进行Nc次延拓，实现对所述震源波场的重建；

假设所述震源波场的正向延拓的总时间为t_Nc，在所述震源波场中确定N个检查点；

在所述震源波场正向延拓的过程中，保存每个检查点时刻及前一时刻的p和q全部波场；

将所述每个检查点时刻及前一时刻的p和q全部波场作为若干部分波场进行存储，其中，p为声压，q为辅助变量；

所述步骤S4，具体为：

在对所述检查点N-1至检查点N时段内波场进行重建时，读取所述检查点N-1时刻及前一时刻的p和q全部波场信息作为初值进行c次正向延拓，每次延拓记录内边界层信息，保存至图形处理器的显存p和q缓冲区，所述显存p和q缓冲区为两个检查点之间p和q内边界大小；

将t_Nc和t_Nc-1两时刻的p和q全部波场作为反向延拓的初值，每延拓一次将所述保存至显存上p和q缓冲区对应时刻内边界信息载入重建波场中，反向延拓c-1次得到重建t_(N-1)c时刻波场；

重复上述步骤重建t_Nc至t₀波场，得到震源波场；

步骤S5，在所述炮集数据中提取地震数据，并对所述地震数据进行插值，得到插值后的地震数据；

步骤S6，根据所述正演方案，通过所述qP波波动方程对所述插值后的地震数据沿时间逆向延拓，共进行Nc次延拓，得到检波点波场；

步骤S7，将所述震源波场和所述检波点波场进行互相关成像，得到一炮偏移剖面，重复上述步骤，得到多炮偏移剖面；

步骤S8，将所述多炮偏移剖面进行叠加，得到最终成像剖面。

2.根据权利要求1所述的方法，其特征在于，

所述步骤S2之后，还包括：

对所述震源子波进行插值，得到插值后的震源子波。

3.根据权利要求1所述的方法，其特征在于，

基于qP波波动方程进行数值模拟计算具体为：

利用高阶有限差分数值方法求解qP波波动方程，采用PML完全匹配层边界条件吸收边界反射，得到模拟地震波场。

4.根据权利要求1～3任一项所述的方法，其特征在于，

在进行所述正向延拓和逆向延拓时，通过图形处理器进行计算。