[发明专利]一种地震波分阶逆时偏移加权叠加成像方法

权利要求书

1.一种地震波分阶逆时偏移加权叠加成像方法，其特征在于它具体包括以下步骤：

(1)基于SRME方法，剔除地震数据中的所有多次波数据，得到纯一次波数据；

(2)基于纯一次波数据，应用SRME方法预测获得纯一阶多次波数据，然后依次预测得出所有N阶多次波数据；

(3)应用逆时偏移成像方法分别对一次波数据和所有N阶多次波数据，根据基于正时和逆时波场能量的归一化互相关条件公式进行成像，得到一次波和N-1阶多次波的N个逆时偏移成像剖面；

a、以炮点处的震源作为正时扰动，进行时间正向传播模拟得到正时波场U₀(x,z,t)(其中x、z、t分别表示水平方向空间坐标、垂向空间坐标、时间坐标，以下同)，同时以检波点处的纯一次波数据作为逆时扰动，进行时间逆向传播模拟得到逆时波场R₀(x,z,t)，并按下式进行成像：

S0(x,z)=ΣtU0(x,z,t)*R0(x,z,t)ΣtU0(x,z,t)*R0(x,z,t),]]>由此得到一次波逆时偏移剖面；

b、以检波点处的纯一次波数据作为正时扰动，进行时间正向传播模拟得到正时波场U₁(x,z,t)，同时以检波点处的纯一阶多次波数据作为逆时扰动，进行时间逆向传播模拟得到逆时波场R₁(x,z,t)，并按下式进行成像：

S1(x,z)=ΣtU1(x,z,t)*R1(x,z,t)ΣtU1(x,z,t)*R1(x,z,t),]]>由此得到一阶多次波逆时偏移剖面；

c、以此类推，以检波点处的纯N-1阶多次波数据作为正时扰动，进行时间正向传播模拟得到正时波场U_N-1(x,z,t)，同时以检波点处的纯N阶多次波数据作为逆时扰动，进行时间逆向传播模拟得到逆时波场R_N-1(x,z,t)，并按下式进行成像:

SN-1(x,z)=ΣtUN-1(x,z,t)*RN-1(x,z,t)ΣtUN-1(x,z,t)*RN-1(x,z,t),]]>由此可得N-1阶多次波成像剖面；

(4)将N个剖面进行加权叠加，得到最终的地震波分阶逆时偏移加权叠加成像剖面。

2.根据权利要求1所述的方法，其特征在于所述步骤(4)为将成像剖面由浅至深划分为M层，根据各个剖面在每层的成像精度确定该剖面每个点的权重系数W_i(x,z)；将各个剖面按如下公式进行加权叠加，得到最终的成像剖面，所述公式为