[发明专利]一种基于局部方差分析的自适应曲波阈值地震去噪方法

权利要求书

1.一种基于局部方差分析的自适应曲波阈值地震去噪方法，其特征在于，所述方法包括以下步骤:

步骤1：对地震数据进行曲波变换获得曲波变换系数；

步骤2：利用最精细尺度下曲波系数的标准差来确定随机噪声在时间域的标准差，并以此计算曲波域中每一尺度及角度下的噪声标准差；

步骤 3：利用含噪信号在每一尺度及角度下的噪声标准差，结合局部方差分析，推算有效信号的标准差；

步骤 4：利用上述获得的噪声与有效信号的标准差建立阈值，进行软阈值处理；

步骤 5：进行曲波逆变换，获得时间域的去噪结果。

2.根据权利要求1所述的基于局部方差分析的自适应曲波阈值地震去噪方法，其特征在于，步骤1中所述的曲波变换采用USFFT或Wrapping方法。

3.根据权利要求1所述的基于局部方差分析的自适应曲波阈值地震去噪方法，其特征在于，步骤3中所述的有效信号的标准差的计算公式为：

式中代表含噪信号的总体方差， j,l 分别代表尺度与角度，代表输入含噪信号的总体标准差，代表噪声标准差。

4.根据权利要求3所述的基于局部方差分析的自适应曲波阈值地震去噪方法，其特征在于，在估算有效信号在曲波域的标准差时，引入一种局部方差处理方法，对指定尺度与角度下的曲波系数矩阵，逐点计算含噪信号的总体方差，即

上式表示计算以目标点(m,n)为中心的一个局部窗口内的方差作为该目标点的方差，进而利用权利要求3所述的公式能产生更准确的有效信号在曲波域的标准差。

5.根据权利要求1所述的基于局部方差分析的自适应曲波阈值地震去噪方法，其特征在于，步骤4中所述的阈值设为：

其中与分别代表尺度和角度，为调节因子，与分别代表随机噪声和有效信号在不同尺度和角度下的标准差。