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

说明书

技术领域

本发明涉及石油勘探反射波地震数据处理过程中的叠前偏移成像技术，具体的说是一种在复杂地质构造下，基于速度分区和线性相位移波场插值的高精度的叠前深度偏移方法。

背景技术

在石油勘探反射波地震数据处理过程中对复杂地质构造成像时，波动方程叠前深度偏移方法得到了越来越广泛的应用，它可以自然地解决波在复杂地层中传播的多路径问题，明显改善复杂构造的成像质量。裂步傅立叶方法(SSF)和傅立叶有限差分法(FFD)是当前波动方程叠前深度偏移的典型算法。SSF方法可以适应缓慢的横向速度变化，但当横向剧烈变速时SSF方法达不到成像精度要求。FFD方法可以适应比较剧烈的横向变速，在处理三维问题时采用多方向分裂的方法，但是分裂方向越细计算量就越大，而且需要大内存，同时具有频散和角度限制问题，该方法目前在生产中应用的成本较高。

发明内容

本发明是提供一种适应复杂地质构造且成本低的高精度叠前深度偏移方法。

本发明的高精度叠前深度偏移方法，具体步骤包括：

1)野外采集的地震数据，通过常规预处理流程，得到统一基准面上的叠前炮集数据；

2)根据地震资料的复杂程度确定延拓步长；

步骤2)所述的延拓步长是通常延拓步长的5到10倍。

步骤2)所述的延拓步长优选为通常步长的5倍。

3)根据速度场横向变化的剧烈程度将速度场分为n个小速度区；

步骤3)所述的将速度场分为n个小速度区是根据已知的每个延拓层的最大速度v_max和最小速度v_min，求出每个分区的临界速度，有n+1个临界速度就有n个小速度区；

步骤3)所述的将速度场分为n个小速度区是将临界速度按从小到大的顺序排列，则最小的临界速度v₁＝v_min，最大的临界速度v_n+1＝v_max，其它的临界速度：v_i＝v_min+(v_max-v_min)/n，i＝2，3，…，n，则第i(i＝1，2，…，n-1)个分区的速度值区间为[v_i，v_i+1)，第n个速度分区的速度值区间为[v_n，v_n+1]，根据n个速度值区间，计算出相应的n个网格点的集合，V_i(i＝1，2，…，n-1)为第i个网格点的集合，则V_i中的元素为该延拓层中所有速度值大于等于v_i且小于v_i+1的网格点，而V_n中的元素为该延拓层中所有速度值大于等于v_n且小于等于v_n+1的网格点；

4)在每一个小速度区内，采用裂步傅立叶方法依据给定的大步长进行波场延拓。

5)对分区延拓后得到的每个速度区的波场，按照实际对应的网格点进行合并，并进行光滑滤波处理，得到延拓后的整体波场，

延拓后的整体波场可由下式计算：

P(x,z+Δz,ω)=Σk=1nηk(x)Pk(x,z+Δz,ω)]]>