[发明专利]变窗口加权地震数据尖峰噪声压制方法

摘要

本发明涉及一种变窗口加权地震数据尖峰噪声压制方法。是针对尖峰噪声具有能量强、出现时间随机和频带宽等特点，噪声源复杂和多变导致多个尖峰噪声具有不同的持续时间和频带宽度，采取固定窗长的尖峰噪声压制方法往往不能良好的压制，且自适应性能较差等诸多问题。经试验，本发明的方法针对受不同持续时间的强尖峰噪声干扰的地震数据具有良好的压制效果，不仅有效地压制了尖峰噪声并且根据检测的尖峰噪声变化窗长自主估计权值得以保留了信号的特征。与现有的中值替代方法相比，对识别尖峰噪声的分辨率更高，充分的保留了信号能量，同时适用于不同的环境情况下，实时分析尖峰噪声的特征并进行压制，有效保证采集的数据质量。

主权项

1.一种变窗口加权地震数据尖峰噪声压制方法，其特征在于，包括以下步骤：/na、读取一道含多尖峰噪声干扰的地震信号x(n)(n＝1,2,...,N)，N为地震信号总采样点数，通过观察分析含尖峰噪声的时域波形，确定尖峰噪声的最短持续时间t_p，信号采样率为Fs；/nb、计算地震信号x(n)的特征函数/n /n利用递归法计算地震信号振幅的短时变化量/n /n长时变化量/n /n式中L_l和L_s为整数，且满足0.01Fs≤L_s≤9Fs·t_p且L_l/L_s＝2；定义/n /nc、定义尖峰噪声的检测阈值/n /nd、满足ψ(n)＞η所有连续采样点区间E＝{E₁,E₂,...,E_K}，K为超过阈值的连续采样点区间的个数，第k个连续采样点区间为和为该连续采样点区间起始采样点序号和终止采样点序号,k∈[1,K]；/ne、对每一个E_k判断其长度是否满足条件/n /n若满足条件，该区间则被认为存在尖峰噪声，否则视为不存在尖峰噪声；/nf、经过步骤e判断得出尖峰噪声的连续采样点区间E′＝{E′₁,E′₂,...,E′_J}，J为尖峰噪声个数，为第j个尖峰噪声段，和为该尖峰噪声所在连续采样点区间的起始采样点序号和终止采样点序号,j∈[1,J]；定义第j个尖峰噪声的持续采样点数的一半为/n /n 为向上取整运算；/ng、对于每一个E′_j的定义两个权值/n /n式中/nh、对于每一个E′_j的采样点，在其两侧各取L_j窗长的信号/n /n令/n /n式(10)中Υ()为计算中值的函数；x′(n)为压制尖峰噪声后的信号。/n