[发明专利]对大构造倾角的倒转盐体叠前深度偏移构造成像方法

权利要求书

1.对大构造倾角的倒转盐体叠前深度偏移构造成像方法，其特征在于采用以下步骤：

1)激发并采集地震波，取得成像空间对应的深度域速度场；

2)根据速度场划分矩形网格；

3)在已经划分的矩形网格上进行射线追踪，而且对射线角度大于90度的射线继续追踪；记录射线在空间位置走过的时间，并将射线没有到达的地方进行插值得到充满空间的旅行时表；

4)应用旅行时对叠前炮集数据进行叠前深度偏移处理，得到共成像点道集；

5)将共成像点道集数据进行叠加得到该工区深度域的构造成像结果，并将结果显示为地层剖面图像。

2.根据权利要求1所述的方法，其特征在于步骤3)所述的射线追踪按照以下公式计算：

x(τ+Δτ)=x(τ)+Δτ6(Kx1+2Kx2+2Kx3+Kx4)]]>

z(τ+Δτ)=z(τ)+Δτ6(Kz1+2Kz2+2Kz3+Kz4)]]>

Px(τ+Δτ)=Px(τ)+Δτ6(Lx1+2Lx2+2Lx3+Lx4)]]>

Pz(τ+Δτ)=Pz(τ)+Δτ6(Lz1+2Lz2+2Lz3+Lz4)]]>

K_x1＝f(τ，x(τ)，P_x(τ))  K_z1＝f(τ，x(τ)，P_z(τ))

L_x1＝f(τ，x(τ)，P_x(τ))  L_z1＝f(τ，x(τ)，P_z(τ))

Kx2=f(τ+Δτ2,x(τ)+12ΔτKx1,Px(τ)+12ΔτLx1)]]>

Kz2=f(τ+Δτ2,z(τ)+12ΔτKz1,Pz(τ)+12ΔτLz1)]]>

Lx2=g(τ+Δτ2,x(τ)+12ΔτKx1,Px(τ)+12ΔτLx1)]]>

Lz2=g(τ+Δτ2,z(τ)+12ΔτKz1,Pz(τ)+12ΔτLz1)]]>

Kx3=f(τ+Δτ2,x(τ)+12ΔτKx2,Px(τ)+12ΔτLx2)]]>

Kz3=f(τ+Δτ2,z(τ)+12ΔτKz2,Pz(τ)+12ΔτLz2)]]>

Lx3=g(τ+Δτ2,x(τ)+12ΔτKx2,Px(τ)+12ΔτLx2)]]>

Lz3=g(τ+Δτ2,z(τ)+12ΔτKz2,Pz(τ)+12ΔτLz2)]]>

K_x4＝f(τ+Δτ，x(τ)+ΔτK_x3，P_x(τ)+ΔτL_x3)

K_z4＝f(τ+Δτ，z(τ)+ΔτK_z3，P_z(τ)+ΔτL_z3)

L_x4＝g(τ+Δτ，x(τ)+ΔτK_x3，P_x(τ)+ΔτL_x3)

L_z4＝g(τ+Δτ，z(τ)+ΔτK_z3，P_z(τ)+ΔτL_z3)

其中，其中x、z是位置坐标分量，v是速度矢量，关于x、z的函数。P_x、P_z是慢度矢量在X、Z方向的分量，也是关于X、Z的函数，τ为当前时间，总是从零时间开始，Δτ为时间步长，τ+Δτ为下一时刻，其函数的表达如下：f(a，b，c)＝a+b²c，