[发明专利]一种平面波叠前深度偏移方法

权利要求书

1.一种平面波叠前深度偏移方法，特征是采用以下步骤：

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

2)根据地震资料的复杂程度确定稀疏频率采样系数α和延拓大步长NΔz，所述Δz为通常延拓步长，N为延拓大步长相对于通常延拓步长的倍数；

3)对频率进行稀疏采样；

利用步骤2中给出的采样系数α，按如下公式对频率进行重采样：

f_m＝α^-mf₁，m＝0，1，…，M

其中f₁为偏移的最小频率，f_n为偏移的最大频率，M为重采样后的频率个数，

4)对重采样后得到的每一个频率按照平面波叠前深度偏移方法以NΔz的大步长进行波场延拓，得到下一个延拓步的波场；

5)对相邻两个大延拓步之间的波场进行线性相位移插值，得到这两个延拓步中间每一个以小延拓步长Δz进行延拓的波场；

6)依据互相关成像原理进行成像，得到最终的叠前深度偏移结果；

7)将偏移结果显示为地层剖面图像。

2.如权利要求1所述的方法，特征是步骤2所述的稀疏频率采样系数α是为重采样后的频率个数为重采样前的一半至三分之一。

3.如权利要求1所述的方法，特征是步骤2所述的N为5至10。

4.如权利要求1或3所述的方法，特征是步骤2所述的N在地震资料复杂时为5。

5.如权利要求1所述的方法，特征是步骤5小延拓步长Δz进行延拓的波场计算是，

深度位置z₁到深度z₁+NΔz之间以小步长step1＝Δz延拓得到的波场可由如下公式计算：

P(x,z1+nΔz,ω)=(N-nN)Pd(x,z1+nΔz,ω)+(nN)Pu(x,z1+nΔz,ω)]]>

其中，P(x，z₁+nΔz，ω)为横向位置x深度位置z₁+nΔz(n＝1，2，…，N-1)处的近似波场，ω为圆频率；P_d(x，z₁+nΔz，ω)为横向位置x处用相位移方法由深度z₁以步长step1向下延拓到深度z₁+nΔz(n＝1，2，…，N-1)处的波场，P_u(x，z₁+nΔz，ω)为横向位置x处用相位移方法由深度z₁+NΔz向上延拓到度z₁+nΔz(n＝1，2，…，N-1)的波场，其计算公式分别为：

Pd(x,z1+nΔz,ω)=P1(x,z1,ω)exp[iωv(x)nΔz]]]>

Pu(x,z1+nΔz,ω)=P2(x,z1+NΔz,ω)exp[-iωv(x)(N-n)Δz]]]>

其中，P₁(x，z₁，ω)为横向位置x深度位置z₁处的波场，P₂(x，z₁+NΔz，ω)为横向位置x深度位置z₁+NΔz处的波场(由波场P₁(x，z₁，ω)以大步长step2＝NΔz延拓得到)，i为虚数，v(x)为横向位置x处的速度。