[发明专利]一种P-SV转换波叠前逆时深度偏移的方法

权利要求书

1.一种P-SV转换波叠前逆时深度偏移的方法，其特征在于，该方法包括以下步骤：

分别进行炮点波场的正向延拓和检波点波场的逆时延拓；

将炮点波场和检波点波场分别进行分解，形成上、下、左、右行波；

将炮点波场和检波点波场中传播方向相反的波互相关成像；以及

完成转换波逆时偏移成像。

2.根据权利要求1所述的方法，其特征在于，利用纵波速度实现炮点纵波波场的正向延拓；利用横波速度实现检波点转换波波场的逆时延拓。

3.根据权利要求2所述的方法，其特征在于，在交错网格空间中用纵波速度对声波方程(1)进行高阶有限差分法离散，进行炮点纵波波场的正向延拓和检波点转换波波场的逆时延拓：

∂p∂t=-K(∂vx∂x+∂vz∂z)∂vx∂t=-1ρ∂p∂x∂vz∂t=-1ρ∂p∂z---(1)]]>

其中K＝ρv²，ρ为密度，v为地震波传播速度，t为时间，x、z为直角坐标系的两个方向，p为位移，v_x、v_z分别为质点在x和z方向的振动速度。

炮点纵波波场正向延拓的差分格式为：

pi,jn=pi,jn-1+ρ(i,j)K(i,j)ΔtΔxΣk=1NDk(N)[vxn-12(i+12+k,j)-vxn-12(i+12-k,j)]+ρ(i,j)K(i,j)ΔtΔzΣk=1NDk(N)[vzn-12(i,j+12+k)-vzn-12(i,j+12-k)]vxn+12(i+12,j)=vxn-12(i+12,j)+Δtρ(i,j)ΔxΣk=1NDk(N)[pn(i+k,j)-pn(i-k+1,j)]vzn+12(i,j+12)=vzn-12(i,j+12)+Δtρ(i,j)ΔzΣk=1NDk(N)[pn(i,j+k)-pn(i,j-k+1)]---(2)]]>

其中：v_p为纵波层速度，Δx、Δz分别为x和z方向的空间离散步长，Δt为时间离散步长，N为差分阶数的一半，i、j为空间离散点序号，n为时间离散点序号，为差分系数；

检波点转换波波场的逆时延拓差分格式：

pi,jn-1=pi,jn-ρ(i,j)K(i,j)ΔtΔxΣk=1NDk(N)[vxn-12(i+12+k,j)-vxn-12(i+12-k,j)]-ρ(i,j)K(i,j)ΔtΔzΣk=1NDk(N)[vzn-12(i,j+12+k)-vzn-12(i,j+12-k)]vxn-12(i+12,j)=vxn+12(i+12,j)-Δtρ(i,j)ΔxΣk=1NDk(N)[pn(i+k,j)-pn(i-k+1,j)]vzn-12(i,j+12)=vzn+12(i,j+12)-Δtρ(i,j)ΔzΣk=1NDk(N)[pn(i,j+k)-pn(i,j-k+1)]---(3)]]>

其中：v_s为横波层速度。