[发明专利]一种消除弹性逆时偏移串音干扰的方法、装置及计算机存储介质

权利要求书

1.一种消除弹性逆时偏移串音干扰的方法，其特征在于，包括以下步骤：

a建立地球物理介质模型；

b给出源小波，利用检波器接收地球物理介质模型中的多分量记录；

c利用1/ω²滤波器对源小波和多分量记录进行滤波，基于滤波后的源小波构造震源波场u，通过散度旋度计算得到分离的源波场P(p_src)和S(s_src)；

d基于滤波后的多分量记录求解双向弹性波动方程，构造接收器波场和额外的接收器波场，得到分离后的接收器波场P波(p_rec)和S波(s_rec)以及额外的接收器波场H_t(p_rec)和H_t(s_rec)；

e基于分离后的源波场P(p_src)和S(s_src)、接收器波场P波(p_rec)和S波(s_rec)以及额外的接收器波场H_t(p_rec)和H_t(s_rec)进行希尔伯特变换，应用成像条件，并对所有分量求和，得到PP、PS、SP和SS图像；

其中，PP是指纵波，即P波入射，P波反射；PS是指转换横波，即P波入射，S波反射；SP是指转换纵波，即S波入射，P波反射；SS是指横波，即S波入射，S波反射；

所述步骤d包括如下步骤：

d1以滤波后的多分量记录为边界条件，求解双向弹性波动方程，构造接收器波场，得到分离后的接收器波场P波(p_rec)和S波(s_rec)；

d2基于滤波后的多分量记录进行希尔伯特变换并以其为边界条件，求解双向弹性波动方程，构造额外的接收器波场，得到分离后的额外接收器波场分别为H_t(p_rec)和H_t(s_rec)；

所述步骤e包括如下具体步骤：

对于扩展的复杂函数ε，希尔伯特算子H既能作用于时间变量t也能作用于空间坐标x或z，在此对深度z进行定义：

ε_z[f(z)]＝f(z)+iH_z[f(z)], (13)

公式(13)中用于表示源端和接收端波场：

其中，ε_t[p_src(x,t)]为对于p_src(x,t)关于时间变量t的函数计算，ε_t[p_rec(x,t)]是p_rec(x,t)关于时间变量t的函数计算，ε_z{ε_t[p_rec(x,t)]}是ε_t[p_rec(x,t)]关于深度z的函数计算，ε_z[p_src(x,t)]是p_src(x,t)关于深度z的函数计算，ε_z{ε_t[s_rec(x,t)]是ε_t[s_rec(x,t)]关于深度z的函数计算；

弹性逆时偏移中PP和PS的最终成像条件为：

其中，H_z(p_src)为对于p_src中空间变量z的希尔伯特变换，H_z(p_rec)为对于P波的接收端中空间变量z的希尔伯特变换，H_t(s_rec)为S波的接收端中时间变量t的希尔伯特变换，H_t(p_rec)为P波的接收端中时间变量t的希尔伯特变换，H_z[H_t(p_rec)]为对于H_t(p_rec)中空间变量z的希尔伯特变换，H_z[H_t(s_rec)]为对于H_t(s_rec)中空间变量z的希尔伯特变换；

忽略被积函数中的空间和时间变量，上述公式(15)的PP和PS方程推广到SP和SS的成像条件：

其中，s_src为S波的震源端信息，p_rec为P波的接收端信息，s_rec为S波的接收端信息，H_z(s_src)为对于s_src空间变量z的希尔伯特变换，H_z(p_rec)为对于P波的接收端中空间变量z的希尔伯特变换，H_t(s_rec)为S波的接收端中时间变量t的希尔伯特变换，H_t(p_rec)为P波的接收端中时间变量t的希尔伯特变换，H_z[H_t(p_rec)]为对于H_t(p_rec)中空间变量z的希尔伯特变换，H_z[H_t(s_rec)]为对于H_t(s_rec)中空间变量z的希尔伯特变换。