[发明专利]三维倾角域稳相叠前深度偏移方法和装置

说明书

本发明提供一种三维倾角域稳相叠前深度偏移方法和装置，涉及地震勘探技术领域，该方法包括获取目标工区的叠前地震数据；根据该叠前地震数据计算走时表以及目标线的深度域倾角道集；该走时表包括炮点走时表和检波点走时表；该深度域倾角道集包括X方向深度域倾角道集和Y方向深度域倾角道集；该目标线为表征目标工区构造形态的测线；根据深度域倾角道集计算目标工区全部成像点的偏移孔径；根据上述走时表和偏移孔径计算目标工区的各个成像点对应不同偏移距的偏移幅值。本发明实施例提供的一种三维倾角域稳相叠前深度偏移方法和装置，可以缓解射线角度分裂导致的菲涅耳带求取困难的问题，增加倾角道集计算的鲁棒性，并获得更高信噪比的偏移结果。

技术领域

本发明涉及地震勘探技术领域，尤其是涉及一种三维倾角域稳相叠前深度偏移方法和装置。

背景技术

在现行叠前偏移方法的实践中，选择偏移孔径是一个重要工作。较小的偏移孔径可减少偏移计算量，但存在着不能对陡倾角构造正确成像的风险；过大的偏移孔径又带来了偏移噪音和较大的偏移计算量。由于偏移算法的实现方式限制，在一个偏移作业中，通常总是选取一个统一的偏移孔径；这样，由于构造倾角是随着空间变化而变化的，对某一成像点合适的孔径，对其他成像点就可能是过大或过小。此外，由于偏移前人们很难对拟成像构造的倾角有一个准确的估计，偏移孔径的选取只能采用保守的方法即选取较大的孔径，从而导致不能很好地压制偏移噪音。

另外，在叠前深度偏移流程中发展稳相偏移方法的主要困难在于倾角的求取。与叠前时间偏移中不同，严格来讲，深度偏移中需要计算瞬时入射射线和瞬时反射射线所夹的法线夹角所决定的平面的瞬时倾角。而在深度偏移中，走时表的求取利用的是射线追踪方法，当速度存在较大横向差异的时候射线发生弯曲，此时需要求出周围点的走时，然后构造出走时梯度方向即射线方向，继而再求取反射面角度。利用这种方法求取出的倾角的连续性较差，对走时表的光滑程度要求很高，稍微不光滑的走时面都会导致射线角度的分裂，即相邻位置处求取的倾角差别很大，非常不利于菲涅耳带的拾取和偏移孔径的筛选。而过于平滑走时面又使得叠前深度偏移方法处理速度横向变化的能力降低。

发明内容

有鉴于此，本发明的目的在于提供一种三维倾角域稳相叠前深度偏移方法和装置，可以缓解射线角度分裂导致的菲涅耳带求取困难的问题，增加倾角道集计算的鲁棒性，并获得更高信噪比的偏移结果，提高计算效率。

第一方面，本发明实施例提供了一种三维倾角域稳相叠前深度偏移方法，包括：获取目标工区的叠前地震数据；根据该叠前地震数据计算走时表以及目标线的深度域倾角道集；该走时表包括炮点走时表和检波点走时表；该深度域倾角道集包括X方向深度域倾角道集和Y方向深度域倾角道集；该目标线为表征目标工区构造形态的测线；根据该深度域倾角道集计算目标工区全部成像点的偏移孔径；根据上述走时表和偏移孔径计算目标工区的各个成像点对应不同偏移距的偏移幅值。

结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，上述根据该叠前地震数据计算目标线的深度域倾角道集的步骤，包括：对于该叠前地震数据中的每一地震道，利用深度域伪倾角计算公式计算该目标线每个CDP的各个成像点的X方向深度域伪倾角和Y方向深度域伪倾角；计算该叠前地震数据中每个地震道的偏移结果；根据上述X方向深度域伪倾角和Y方向深度域伪倾角的大小对该偏移结果进行分选和叠加，利用深度域倾角道集计算公式计算得到该目标线每个CDP的深度域倾角道集。