[发明专利]基于极化滤波的多波多分量地震资料的去噪方法

权利要求书

1.一种基于极化滤波的多波多分量地震资料的去噪方法，其特征在于， 该方法包括以下步骤：

获取多波多分量地震资料的纵波资料和转换波资料；

在预设的频带范围内对获取的纵波资料和转换波资料进行低通滤波，以 获得纵波初始数据和转换波初始数据；

A、根据所述多波多分量地震资料的纵波初始数据和转换波初始数据， 采用视速度差异、F-K滤波或者τ-p变换获得面波初始水平参量和初始垂直 参量；

B、根据所述面波初始水平参量和初始垂直参量，利用极化滤波迭代方 法计算面波目标水平参量和面波目标垂直参量，所述极化滤波迭代方法包括：

B1、利用极化椭圆参数求解方法求解得到多个椭圆；

B2、对所述求解的多个椭圆进行全局优选得到所需的椭圆(P)，并将椭 圆(P)的长轴(a)赋值给面波初始垂直参量及短轴的负值(-b)赋值给面 波初始水平参量；

B3、根据所述面波初始水平参量和面波初始垂直参量分别计算转换波初 始数据与面波初始水平参量差值(R₀(N)-R_p′(N))的振幅包络(E_rp(N))以及纵 波初始数据与面波初始垂直参量差值(Z₀(N)-Z_p′(N))的振幅包络(E_zp(N))；

B4、根据所述振幅包络E_rp(N)和E_zp(N)利用互相关公式来分别计算此时 面波水平参量的振幅包络(E′_r)和面波垂直参量振幅包络(E′_z)；

B5、根据所述E′_z和E′_r来计算当所述(R₀(N)-R_p′(N))与E′_r的模量值 (||(R₀(N)-R_p(N))，E′_r||)以及(Z₀(N)-Z_p′(N))与E′_z的模量值(||(Z₀(N)-Z_p(N))，E′_z||) 为最小值时的面波水平参量和面波垂直参量，当所述两个模量值都达到最小 时计算得到的所述面波水平参量和面波垂直参量即为面波目标水平参量和目 标垂直参量，否则将求得的所述面波水平参量和面波垂直参量分别赋值给面 波初始水平参量和面波初始垂直参量后依次重复步骤B3、B4和B5；

C、从所述纵波初始数据中减去所述面波目标垂直参量，及从所述转换 波初始数据中减去所述面波目标水平参量，以得到去除噪音的数据。

2.根据权利要求1所述的方法，其特征在于，所述极化椭圆参数求解方 法利用最小平方二乘法实现。