[发明专利]基于3D-Shearlet变换和尺度自适应的地震随机噪声压制方法

权利要求书

1.一种基于3D-Shearlet变换和尺度自适应的地震随机噪声压制方法，其特征在于，包括如下步骤：

步骤1、读取多炮二维地震数据；对多炮二维地震数据进行3D-Shearlet变换，获得多炮二维地震数据在3D-Shearlet域各尺度和方向上的变换系数C(j，l，k)；

步骤2、设定去除随机噪声的自适应阈值函数，根据不同尺度和方向下的3D-Shearlet变换系数，计算尺度自适应阈值；

步骤3、对3D-Shearlet变换系数C(j，l，k)进行阈值化处理，保留对应尺度和方向大于自适应阈值的变换系数，设置对应尺度和方向小于自适应阈值的变换系数为零，进而得到阈值后3D-Shearlet变换系数；

步骤4、利用3D-Shearlet反变换，对阈值后3D-Shearlet变换系数进行3D-Shearlet反变换处理，获取随机噪声压制后的结果。

2.根据权利要求1所述的基于3D-Shearlet变换和尺度自适应的地震随机噪声压制方法，其特征在于，步骤1中，所用3D-Shearlet变换系数算法如下：

式中，C(j，l，k)为3D-Shearlet变换系数；3DSH{·}表示3D-Shearlet正变换；·，·表示内积计算；f为待分解信号；

表示为：

式中，j，l，k分别表示尺度参数、方向参数和平移参数；A表示各向异性矩阵；B表示剪切矩阵；F为尺度方向函数；d＝1，2，3；

其中W表示如下：

式中，为Meyer小波的尺度函数。

3.根据权利要求1所述的基于3D-Shearlet变换和尺度自适应的地震随机噪声压制方法，其特征在于，步骤2中，自适应阈值函数T表达如下：

T(j，l，k)＝ε·H·λ·σ (4)

式中，H为常数变量；

λ表示某一尺度和方向下3D-Shearlet变换系数的均方根；

σ表示噪声标准差；

ε为尺度自适应因子，表示为：

ε＝2^-(J-j)/2 (5)

式中，J为尺度值。

4.根据权利要求1所述的基于3D-Shearlet变换和尺度自适应的地震随机噪声压制方法，其特征在于，步骤3中，阈值化处理表示如下：

式中，表示阈值后3D-Shearlet变换系数。

5.根据权利要求1所述的基于3D-Shearlet变换和尺度自适应的地震随机噪声压制方法，其特征在于，步骤4中，3D-Shearlet反变换表达如下：

式中，表示随机噪声压制后的地震数据；

3DSH^-1{·}表示3D-Shearlet反变换。