基于时变窗长时频峰值滤波的地震信号去噪研究.ppt

基于时变窗长时频峰值滤波的地震信 号去噪研究 地震勘探是利用人工激发的地震波，探测地下地质构造，来查 明石油、煤炭、天然气等的分布和储量的一种方法。在野外采 集的地震数据中，由于随机干扰和相干噪声的存在，使得一些 有效波被淹没在噪声中，采集到的地震勘探资料常常无法使用， 所以提高地震资料的信噪比是地震资料处理的首要任务。 地震信号中的随机噪声通常是由各种不可预知的因素综合作用 而形成的，没有一定的规律，频带相对较宽，几乎覆盖整个有 效波的频带。 2004 年 B.Boashash 提出时频峰值滤波（ TFPF ）算 法，它是提取有效信号、压制随机噪声的一种有效算法，近几 年已经应用于地震勘探资料处理中。 时频峰值滤波是一种基于瞬时频率估计的信号增强算法，该算 法首先通过频率调制将含噪信号调制成解析信号的瞬时频率， 然后取解析信号的 Wigner-Ville 分布的峰值来进行瞬时频率估 计，从而恢复有效信号。针对实际地震资料的非线性特性，本 文选用加窗的 Wigner-Ville 分布即伪 Wigner-Ville 分布（ PWVD ） 来实现 TFPF ，从而提高被处理信号的线性度，减小瞬时频率估 计的偏差。 对于 TFPF 来说，时频分布（ TFD ）中的窗长是一个非常重要的参 数，窗长大小的选取影响着 TFPF 的滤波效果。不同频率的信号 需要不同的最优窗长，最优窗长的选取需要很好的平衡有效信 号的恢复和随机噪声的压制。 窗取得越短，窗内的瞬时频率（ IF ）越近似线性，这样可以很 好的恢复有效信号，但是其在频域对噪声压制的能力会下降。 窗取得越长，窗内噪声的统计特性就越接近于真实高斯白噪声 的零均值特性，这样就会增强对噪声的压制能力，但是窗内瞬 时频率的线性度降低，导致瞬时频率估计的偏差增大，造成有 效信号幅值的衰减。 尤其是在信号的波峰和波谷的位置，解析信号的瞬时频率达到 最大值，信号变化最剧烈，在很难保证瞬时频率线性条件的情 况下，选用传统的固定窗长时频峰值滤波很难同时达到有效信 号幅值的良好恢复和随机噪声的有效压制。对于本文研究的地 震信号来说，每一道不同时刻地震信号的非线性程度和随机噪 声的强度是不同的，所以当采用传统的固定窗长时频峰值滤波 进行信号处理时，在有效压制随机噪声的同时容易引起有效信 号幅值的衰减，从而会丢失部分有效信号同相轴。 本文在窗长度的选取上综合考虑信号的非线性程度和噪声强度 两方面的因素，对基于 PWVD 的瞬时频率估计方法进行改进，根 据瞬时频率偏差和方差的近似公式并利用瞬时频率置信区间交 叉准则，提出时变窗长时频峰值滤波算法，并将其应用于地震 勘探资料的处理中。为了验证本文改进算法的有效性和实用性， 进行了大量的仿真实验，将该算法分别应用于模拟地震记录和 实际地震资料的处理。