[发明专利]基于小波分解的自适应中值强噪声干扰地震数据去噪方法

权利要求书

1.基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，该基于小波分解的自适应中值强噪声干扰地震数据去噪方法包括：

步骤1，对地震数据进行小波分解，对分解结果进行初始中值滤波器处理；

步骤2，通过以离散余弦变换为基的滤波器对数据进行卷积处理，估计得到的响应数据的方差和峰度值，进而估计噪声方差；

步骤3，在参考中值滤波器窗口的基础上，根据噪声水平估计的不同进行滤波器窗口长度的合理取选，来定义中值滤波器窗口长度；

步骤4，使用自适应中值滤波器对数据进行滤波处理，根据不同区域噪声水平的大小自适应的调节中值滤波器的窗口长度；

步骤5，对中值滤波处理后的数据进行重构，得到去噪后的地震数据。

2.根据权利要求1所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，该基于小波分解的自适应中值强噪声干扰地震数据去噪方法还包括，在步骤1之前，输入地震数据，进行数据分析。

3.根据权利要求1所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，在步骤1中，对地震数据进行小波分解，将地震数据变换到小波域进行处理，得到去除了大部分混合噪声的地震数据；对处理结果用初始固定窗口长度的中值滤波器进行滤波处理。

4.根据权利要求3所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，在步骤1中，使用滤波窗长度较大的参考中值滤波器，对经过小波分解的子频带进行粗略的滤波处理，初始滤波器长度的选择要使得大部分混合噪声被去除，不用管有效信号的微小损失，初始滤波器长度在7到11之间选取。

5.根据权利要求1所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，在步骤2中，通过离散余弦变换为基的多个滤波器对数据进行卷积，生成响应数据；估计得到各个响应数据的方差和峰度值；对得到的响应数据的方差和峰度值，忽略离散余弦分量，来估计噪声方差。

6.根据权利要求5所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，在步骤2中，对地震数据进行分块，N表示分块的大小，用N×N个以离散余弦变换为基的滤波器对分块数据进行卷积生成响应数据，对在范围(2…N²)内的每个相应数据求得其方差和峰度值，以此来得到经过常规中值初始滤波器处理后的数据的局部噪声方差估计，从而得到不同位置的噪声水平估计。

7.根据权利要求1所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，在步骤3中，针对地震数据中不同区域的不同噪声方差，在定义的参考中值滤波器窗口的基础上，利用估计噪声方差定义自适应滤波器窗口长度，参考滤波器窗口将会根据噪声水平估计的不同进行取。

8.根据权利要求1所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，该基于小波分解的自适应中值强噪声干扰地震数据去噪方法还包括，在步骤3之后，计算自适应滤波器的映射。

9.根据权利要求1所述的基于小波分解的自适应中值强噪声干扰地震数据去噪方法，其特征在于，在步骤4中，使用自适应中值滤波器对数据进行滤波处理，为了适应不同区域，采用不同的窗口长度，有效信号位置应用小的时变滤波器窗口，以消除噪声时保护有效信号，而在噪声位置应用较大的时变滤波器窗口以用来衰减噪声能量。