[发明专利]基于经验周期模态单元滤波的地震信号消噪方法

说明书

一、技术领域：

本发明涉及信号消除噪音的方法，特别是一种地震信号的消噪的方法。

二、背景技术：

地震信号被广泛应用于地质勘探，地震预报等应用领域。由于地震信号采集环境的不同和采集仪器中存在各种干扰，所采集到的地震波形中通常包含了很多噪声信号，严重影响了对地震信号的进一步应用。要提高地震资料的信噪比，就需要去除地震波形中的随机噪声和想干噪声。目前常用的用于消除信号噪声的方法有：基于IIR的滤波方法、FIR滤波方法、基于小波阈值滤波方法等。但用上述方法时需要反复修改参数或基函数，并需要不断观察消噪的效果是否理想。

三、发明内容：

本发明的目的是克服上述不足问题，提供一种地震信号的消噪的方法。该方法采用自适应消噪方法，抛弃了传统阈值消噪方法中以单个信号点为阈值处理对象的思想，根据经验周期模态单元滤波的原理，该方法不需要先验基函数和不断改变和修改参数。方法简单、可行。

本发明为实现上述目的所采用的技术方案是：基于经验周期模态单元滤波的地震信号消噪方法，首先采集地震波形图和频谱图，利用MATLAB软件中的imread命令进行读取处理，初始预定信号模拟单元的经验周期范围，确定两个经验周期的阈值；然后根据确定的阈值和模态单元的处理系数确定去除噪声后的波形形状，最后利用MATLAB软件中的plot命令处理结果显示出来。

所述确定阈值和模态单元的方法为：先将信号模态单元的经验周期范围初始化为[thpi^-，thpi⁺]，把imf_i中第k个模态单元记为m_i^k，其中i＝1，2，...，I。m_i^k的经验周期阈值处理方法为：

mik(thpi-,thpi+)=χ(thpi-,thpi+)mik,]]>其中

χ(thpi-,thpi+)=1thpi-≤pik≤thpi+0otger]]>