[发明专利]一种对地震资料剔除野值的滤波方法

说明书

技术领域

本发明属于地震地质领域中对地震资料进行信号处理的领域，尤其涉及一种对信号中野值信号进行有效剔除的方法。

背景技术

在地震资料中，野值信号的频率和相位是不确定的，它可能是傅立叶(周期波动)信号，也可能是小波(非周期波动)信号。因此，针对频率和相位的各种滤波方法对剔除野值来说都是收效甚微的。野值信号的唯一特征是其相对振幅极大，但是，野值的振幅门槛同样无法确定。单独采用平滑滤波(如中值滤波)可有效剔野值，然而，对其它傅立叶信号来说，却可能出现频率或相位畸变。现有技术在实数道内自动剔除野值是非常困难的，目前，各地震资料处理系统都有自动剔除野值的功能，但其效果均不理想，均无法有效剔除野值同时不产生信号的频率或相位畸变。

发明内容

本发明为了解决现有技术中存在的技术问题，研发了一种复数道剔除野值的方法，该方法能有效剔除野值，并保留有效波信号及野值信号的频率和相位，地震道信号无畸变现象出现。

本发明所采用的方法为针对地震资料信号中存在的野值信号进行剔除和滤波，并保留有效波信号及野值信号的频率和相位；本发明采用地震复数道滤波的方法，包括对地震资料信号进行信号变换和分解步骤以及平滑滤波步骤，并将平滑滤波后的瞬时包络与原始信号波相位函数重构得到剔除野值波的输出信号。

具体的说，本发明方法包括如下步骤，

①.地震资料信号输入；

②.变换分解输入信号：采用希尔伯特变换分解方法

③.瞬时包络平滑滤波步骤；

④.将滤波后的瞬时包络与原始相位函数重构并输出信号：

即重构步骤为滤波后的瞬时包络与原始相位函数cosθ(t)的积：

yr(t)=cosθ(t)·a(t)‾]]>

其中，y_r(t)为已剔除野值的且无相位和频率畸变的最终输出。

其中，所述步骤二中采用希尔伯特变换分解方法为

实函数x_r(t)的希尔伯特变换为x_h(t)；

xh(t)=1π∫-∞+∞xr(τ)t-τdτ=xr(t)*1πt]]>

式中：t为时间序列，*为卷积；