[发明专利]一种基于同相轴形态差异的面波提纯方法及装置

权利要求书

1.一种基于同相轴形态差异的面波提纯方法，其特征在于，所述方法包括，

S10：对地震反射波数据y_r进行初始化；

S20：根据原始地震数据y及所述反射波数据y_r，求取面波在频域高分辨率LRT变换域的稀疏表示系数，并重构面波数据y_g，包括：

根据原始地震数据y以及反射波数据y_r计算获得第一残差数据Δy₁，其中，Δy₁＝y-y_r；

对第一残差数据Δy₁进行频域高分辨率LRT变换获得第一变换系数z_g；

对获得的第一变换系数z_g利用第一阈值σ_g进行硬阈值处理得到面波在频域高分辨率LRT变换域的稀疏表示系数z_g'；

对面波稀疏表示系数z_g'进行频域反LRT变换，重构面波数据y_g；

S30：根据所述原始地震数据y及重构后的所述面波数据y_g，求取反射波在时域高分辨率HRT变换域的稀疏表示系数，并重构反射波数据y_r；

迭代进行步骤S20、步骤S30以获得面波提纯结果。

2.如权利要求1所述的基于同相轴形态差异的面波提纯方法，其特征在于，所述根据所述原始地震数据y及重构后的面波数据y_g，求取反射波在时域高分辨率HRT变换域的稀疏表示系数，并重构反射波数据y_r的步骤包括：

根据原始地震数据y以及重构后的面波数据y_g计算获得第二残差数据Δy₂，其中Δy₂＝y-y_g；

对第二残差数据Δy₂进行时域高分辨率HRT变换得到第二变换系数z_r；

对获得的第二变换系数z_r利用第二阈值σ_r进行硬阈值处理得到反射波在时域高分辨率HRT变换域的稀疏表示系数z_r'；

对反射波稀疏表示系数z_r'进行时域反HRT变换，重构反射波数据y_r。

3.如权利要求2所述的基于同相轴形态差异的面波提纯方法，其特征在于，设定一循环变量k、给定迭代次数K，循环变量k的初始值k＝0，在所述步骤S20之前，所述方法还包括：

令循环变量k＝k+1。

4.如权利要求3所述的基于同相轴形态差异的面波提纯方法，其特征在于，在每次迭代过程中，更新第一阈值σ_g、第二阈值σ_r的值，具体包括：

其中，k为循环变量，K为迭代次数。

5.如权利要求1所述的基于同相轴形态差异的面波提纯方法，其特征在于，所述对地震反射波数据y_r进行初始化包括：

令反射波数据初始值为y_r＝0。

6.如权利要求5所述的基于同相轴形态差异的面波提纯方法，其特征在于，所述迭代进行步骤S20、步骤S30以获得面波提纯结果的步骤包括：

当循环变量k小于或等于迭代次数K时，重复执行步骤S20～S30以获得面波提纯结果。

7.如权利要求5所述的基于同相轴形态差异的面波提纯方法，其特征在于，所述方法还包括：

当循环变量k大于迭代次数K时，停止迭代，输出面波提纯数据。

8.如权利要求1所述的基于同相轴形态差异的面波提纯方法，其特征在于，在对面波数据y_g、反射波数据y_r进行初始化之前，所述方法还包括：

获取原始地震数据y。

9.一种基于同相轴形态差异的面波提纯装置，其特征在于，所述装置用于执行如权利要求1～8任意一项所述的基于同相轴形态差异的面波提纯方法，所述装置包括：

初始化单元，用于对反射波数据y_r进行初始化；

第一重构单元，用于根据原始地震数据y及所述反射波数据y_r，求取面波在频域高分辨率LRT变换域的稀疏表示系数，并重构面波数据y_g；

第二重构单元，用于根据所述原始地震数据y及重构后的所述面波数据y_g，求取反射波在时域高分辨率HRT变换域的稀疏表示系数，并重构反射波数据y_r；

所述第一重构单元及所述第二重构单元还用于进行多次迭代重构以获得面波提纯结果。