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

摘要

本发明提供的针对VTI介质逆时偏移方法、设备及介质，方法为：在炮集数据中提取震源坐标，并生成震源子波；根据正演方案，通过qP波波动方程对震源子波沿时间正向延拓，在延拓过程中选择若干部分波场进行存储；通过qP波波动方程利用若干部分波场沿时间逆向延拓，得到震源波场；在炮集数据中提取地震数据，通过qP波波动方程对插值后的地震数据沿时间逆向延拓，得到检波点波场；将震源波场和检波点波场进行互相关成像，得到一炮偏移剖面，重复上述步骤，得到多炮偏移剖面；将多个炮偏移剖面进行叠加，得到最终成像剖面。本发明在大幅减小计算量的情况下，能够获得相同运动学特征的纵波分量，使成像剖面的反射同相轴准确归位，提升复杂介质的成像精度。

主权项

1.一种针对VTI介质逆时偏移方法，其特征在于，包括：/n步骤S1，获取炮集数据；/n步骤S2，在所述炮集数据中提取震源坐标，并生成震源子波；/n步骤S3，根据正演方案，将震源子波加载至波场中，得到震源波场，通过qP波波动方程对所述震源波场沿时间正向延拓，共进行Nc次延拓，Nc为t₀至t_Nc时间段内设定的延拓次数，并在延拓过程中选择若干部分波场进行存储；/n在所述qP波波动方程中引入稳定性因子σ；/n利用稳定性因子公式实现对横波沿着对称轴方向的传播速度V_sz的计算，其中，通过数值模拟，所述稳定性因子σ为0.75，ε和δ为Thomsen参数，V_pz为纵波沿着对称轴方向的传播速度；/n所述步骤S3中，根据正演方案，通过qP波波动方程对所述震源子波沿时间正向延拓，具体为：/n根据所述纵波和横波沿着对称轴方向的传播速度以及Thomsen参数作为基于qP波波动方程的VTI介质逆时偏移成像方法的模型参数；/n将地震子波作为激发震源生成地震波场，利用有限差分数值解法，对所述震源波场沿时间正向Nc次延拓，Nc为t₀至t_Nc时间段内设定的延拓次数；/n其中，qP波波动方程如以下公式所示，/n /n /n /n /n式中，V_pz和V_sz分别为纵波和横波沿着对称轴方向的传播速度，ε和δ为Thomsen参数，V_px为纵波沿着垂直于对称轴方向的传播速度，V_pn为动校正速度，p为声压，q为辅助变量；/n步骤S4，根据所述正演方案，通过所述qP波波动方程对所述若干部分波场沿时间逆向延拓，共进行Nc次延拓，实现对所述震源波场的重建；/n假设所述震源波场的正向延拓的总时间为t_Nc，在所述震源波场中确定N个检查点；/n在所述震源波场正向延拓的过程中，保存每个检查点时刻及前一时刻的p和q全部波场；/n将所述每个检查点时刻及前一时刻的p和q全部波场作为若干部分波场进行存储，其中，p为声压，q为辅助变量；/n所述步骤S4，具体为：/n在对所述检查点N-1至检查点N时段内波场进行重建时，读取所述检查点N-1时刻及前一时刻的p和q全部波场信息作为初值进行c次正向延拓，每次延拓记录内边界层信息，保存至图形处理器的显存p和q缓冲区，所述显存p和q缓冲区为两个检查点之间p和q内边界大小；/n将t_Nc和t_Nc-1两时刻的p和q全部波场作为反向延拓的初值，每延拓一次将所述保存至显存上p和q缓冲区对应时刻内边界信息载入重建波场中，反向延拓c-1次得到重建t_(N-1)c时刻波场；/n重复上述步骤重建t_Nc至t₀波场，得到震源波场；/n步骤S5，在所述炮集数据中提取地震数据，并对所述地震数据进行插值，得到插值后的地震数据；/n步骤S6，根据所述正演方案，通过所述qP波波动方程对所述插值后的地震数据沿时间逆向延拓，共进行Nc次延拓，得到检波点波场；/n步骤S7，将所述震源波场和所述检波点波场进行互相关成像，得到一炮偏移剖面，重复上述步骤，得到多炮偏移剖面；/n步骤S8，将所述多炮偏移剖面进行叠加，得到最终成像剖面。/n