[发明专利]一种基于双树复小波域的地震信号降噪方法及系统

权利要求书

1.一种基于双树复小波域的地震信号降噪方法，包括：

步骤a：建立含噪地震信号模型，对含噪地震信号进行双树复小波域分解变换，得到每一尺度每个方向的含噪小波系数；

步骤b：对同一方向实部小波系数或虚部小波系数与其对应的模建立双变量模型得到阈值函数，利用阈值函数对双树复小波域中的含噪小波系数进行双变量阈值处理，得到不含噪小波系数的估计值；

步骤c：利用双树复小波逆变换对不含噪小波系数的估计值进行重构得到降噪后的原始地震信号。

2.根据权利要求1所述的基于双树复小波域的地震信号降噪方法，其特征在于，在所述步骤a中，所述含噪地震信号模型为：

g＝x+v

其中，x为原始地震信号，v是服从分布为的高斯白噪声，g为含噪地震信号。

3.根据权利要求1所述的基于双树复小波域的地震信号降噪方法，其特征在于，在所述步骤a中，所述对含噪地震信号进行双树复小波域分解变换还包括：在双树复小波域中，含噪地震信号模型为：

y＝w+n

其中，y为含噪小波系数，w为原始地震信号小波系数，n为噪声的小波系数。

4.根据权利要求1所述的基于双树复小波域的地震信号降噪方法，其特征在于，在所述步骤a中，所述对含噪地震信号进行双树复小波域分解变换具体为：采用双树算法对地震信号进行二维双树复小波分解变换，二维双树复小波变换通过一维双树复小波的张量积得到；变换过程中，各尺度分别生成2个低频子带系数和方向分别为±15°、±45°和±75°的6个不同高频子带系数，且每个方向有实部小波系数与虚部小波系数两个部分。

5.根据权利要求1所述的基于双树复小波域的地震信号降噪方法，其特征在于，在所述步骤b中，所述双变量模型具体为：

y1=w1+n1m=u+n0]]>

其中，含噪地震信号小波系数的模y₁是实部小波系数，y₂为对应的虚部小波系数，u为原信号小波系数的模，n₀为噪声小波系数的模。