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

权利要求书

1.一种三维倾角域稳相叠前深度偏移方法，其特征在于，包括：

获取目标工区的叠前地震数据；

根据所述叠前地震数据计算走时表；并对于所述叠前地震数据中的每一地震道，利用深度域伪倾角计算公式计算目标线每个CDP的各个成像点的X方向深度域伪倾角和Y方向深度域伪倾角；计算所述叠前地震数据中每个地震道的偏移结果；根据所述X方向深度域伪倾角和所述Y方向深度域伪倾角的大小对所述偏移结果进行分选和叠加，利用深度域倾角道集计算公式计算得到所述目标线每个CDP的深度域倾角道集；所述走时表包括炮点走时表和检波点走时表；所述深度域倾角道集包括X方向深度域倾角道集和Y方向深度域倾角道集；所述目标线为表征所述目标工区构造形态的测线；

根据所述深度域倾角道集计算所述目标工区全部成像点的偏移孔径；

根据所述走时表和所述偏移孔径计算所述目标工区的各个成像点对应不同偏移距的偏移幅值。

2.根据权利要求1所述的三维倾角域稳相叠前深度偏移方法，其特征在于，所述深度域伪倾角计算公式、所述深度域倾角道集计算公式为：

其中，θ_x为X方向深度域伪倾角，θ_y为Y方向深度域伪倾角，x、y是成像点的横向坐标，z是成像点的成像深度，τ_s为炮点到成像点的走时，τ_g为检波点到成像点的走时，r_s为炮点到成像点的距离，r_g为检波点到成像点的距离，I(x,y,z,θ_x)为X方向深度域倾角道集，I(x,y,z,θ_y)为Y方向深度域倾角道集，N_trace是所述叠前地震数据包含的地震道总数目，f′_m(τ_s+τ_g；x_s,y_s,x_g,y_g)为地震道的一阶导数，x_s、y_s为炮点横向坐标，x_g、y_g为检波点横向坐标。

3.根据权利要求1所述的三维倾角域稳相叠前深度偏移方法，其特征在于，所述根据所述深度域倾角道集计算所述目标工区全部成像点的偏移孔径的步骤，包括：

根据所述深度域倾角道集计算目标CDP的稳相点；所述目标CDP为表征所述目标工区目标线构造形态的CDP；

根据所述深度域倾角道集和所述稳相点计算所述目标CDP的成像点的偏移孔径；

根据所述目标CDP的成像点的偏移孔径，通过插值和平滑得到所述目标工区全部成像点的偏移孔径。

4.根据权利要求3所述的三维倾角域稳相叠前深度偏移方法，其特征在于，所述根据所述深度域倾角道集计算目标CDP的稳相点的步骤，包括：

沿角度方向分别叠加X方向深度域倾角道集、叠加Y方向深度域倾角道集得到对应方向的偏移剖面；

根据所述偏移剖面，利用第一关系式沿深度方向分别对X方向、Y方向的深度域倾角道集进行振幅值系数度量，得到对应方向的反射同相轴的深度系列Z_j；其中，j＝1,2...J,J为反射同相轴的个数；

利用第二关系式分别检验X方向、Y方向的深度域倾角道集中所述深度系列Z_j的反射同相轴的同相性，当所述反射同相轴满足所述第二关系式时，记录所述反射同相轴的稳相点，分别得到目标CDP在对应方向的稳相点集合。