[发明专利]一种高精度的叠前深度偏移方法

摘要

本发明涉及石油勘探反射波地震数据处理过程中的高精度的叠前深度偏移方法，包括以下步骤：采集得到统一基准面上的叠前炮集数据，确定延拓步长，将速度场分为小速度区，在每一小速度区内用裂步傅立叶方法依据大步长进行波场延拓，5)对分区延拓后得到的每个速度区的波场，按照实际对应的网格点进行合并，并进行光滑滤波处理，得到延拓后的整体波场，对两个相邻波场进行线性相位移插值，得到中间每一个以小延拓步长进行延拓的波场。本发明依据互相关成像原理进行成像，得到叠前深度偏移结果，并将偏移结果显示为地层剖面图像，对复杂地质构造成像质量明显提高且成本低。

主权项

一种高精度的叠前深度偏移方法，其特征在于包括以下步骤：1)野外采集的地震数据，通过常规预处理流程，得到统一基准面上的叠前炮集数据；2)根据地震资料的复杂程度确定延拓步长；3)根据速度场横向变化的剧烈程度将速度场分为n个小速度区；4)在每一个小速度区内，采用裂步傅立叶方法依据给定的大步长进行波场延拓；5)对分区延拓后得到的每个速度区的波场，按照实际对应的网格点进行合并，并进行光滑滤波处理，得到延拓后的整体波场，延拓后的整体波场由下式计算： P ( x , z + Δz , ω ) = Σ k = 1 n η k ( x ) P k ( x , z + Δz , ω ) 其中P(x，z+Δz，ω)为延拓后得到的横向位置x深度位置z+Δz处的波场，ω为圆频率，n为步骤(3)中确定的分区的个数，Pk(x，z+Δz，ω)为第k(k＝1，2，…，n)个分区从深度z延拓到深度z+Δz处的波场，ηk(x)是关于自变量x的函数，算式如下： η k ( x ) = 1 , if x ∈ V k 0 , if x ∉ V k , k = 1,2 , . . . , n 其中Vk(k＝1，2，…，n)为第k个分区的网格点的集合；6)对两个相邻波场进行线性相位移插值，得到中间每一个以小延拓步长进行延拓的波场；7)依据互相关成像原理进行成像，得到叠前深度偏移结果，并将偏移结果显示为地层剖面图像。