[发明专利]挤压短时傅里叶变换的地震勘探信号随机噪声压制方法

说明书

本发明公开了一种挤压短时傅里叶变换的地震勘探信号随机噪声压制方法，首先对高维地震数据体中的每一道信号进行短时傅里叶变换，然后对大于设定阈值1的短时傅里叶变换系数进行挤压操作得到挤压短时傅里叶变换系数，接着利用相邻道挤压短时傅里叶变换系数的位置为约束，对当前道挤压短时傅里叶变换的系数进行阈值处理，最后将处理完的短时傅里叶变换系数进行重构得到压制随机噪声后的地震勘探信号，为后续地震数据处理也解释提供高信噪比的基础数据。

技术领域

本发明属于勘探地球物理技术领域，具体涉及一种利用挤压短时傅里叶变换的地震勘探信号随机噪声压制方法。

背景技术

随机噪声主要来源于地面微震、仪器以及由于介质不均匀性造成的散射波。它在地震记录上表现为杂乱无章的振动，并且频带很宽。同时由于复杂的地表和地形条件，如山地、黄土塬、崎岖海底、沙漠、戈壁和多期叠合构造等因素，对深层地震记录产生各种随机噪声干扰。随机噪声干扰会降低地震记录的信噪比及质量，对后续的地震信号处理及解释环节造成很大的影响。因此为了更加充分和准确的应用采集到的地震数据，需要对勘探地震数据中的随机噪声进行压制以提高地震信号的信噪比，进而保证后续处理解释环节的可靠性。

目前已有的地震勘探信号随机噪声压制方法包含傅里叶域滤波类方法，频率-波数域滤波类方法，时频域变换滤波类方法和其它二维小波变换类滤波方法，针对傅里叶域滤波类方法，宽频带地震信号与随机噪声在频率域相互交叠，并不能将两者完全分离；同时随机噪声在时域的分布并不完全均匀，直接在频率域压制噪声并不完全且对信号损伤较大；针对频率-波数域滤波类方法，宽频带地震信号与随机噪声在频率-波数域相互交叠，相较于傅里叶域滤波方法交叠减少，但也不能完全将两者完全分离；同时随机噪声在时间-空间域的分布并不完全均匀，直接在频率-波数域压制噪声并不完全且对信号损伤较大。针对时频域变换滤波类方法，计算效率较低且压制噪声的效果依赖于时频变换工具，时频变换的分辨率及聚集度对滤波影响很大。针对其它二维小波变换类滤波方法，计算量大，且各种变换的分辨率受不确定性原理的约束都不高，因此信号的系数在变换域并不聚焦。

发明内容

本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种挤压短时傅里叶变换的地震勘探信号随机噪声压制方法，对高维地震数据体中的每一道信号进行短时傅里叶变换，对大于某个阈值的短时傅里叶变换系数进行挤压操作得到挤压短时傅里叶变换系数，然后利用相邻道挤压短时傅里叶变换系数的位置和大小为约束，同时利用阈值对当前道挤压短时傅里叶变换的系数进行阈值等处理，将处理完的短时傅里叶变换系数进行重构得到压制随机噪声后的信号，为后续地震数据处理也解释提供高信噪比的基础数据。

本发明采用以下技术方案：

挤压短时傅里叶变换的地震勘探信号随机噪声压制方法，先对高维地震数据体中的每一道信号进行短时傅里叶变换；然后对大于设定阈值1的短时傅里叶变换系数进行挤压操作得到挤压短时傅里叶变换系数；再利用相邻道挤压短时傅里叶变换系数的位置和大小为约束，对当前道挤压短时傅里叶变换的系数进行阈值处理；最后将处理完的短时傅里叶变换系数进行重构得到压制随机噪声后的信号，为后续地震数据处理也解释提供高信噪比的基础数据，重复以上步骤直到二维数据体中所有的Nx道地震信号全部被压制随机噪声。

具体的，对高维地震数据体中的每一道信号进行短时傅里叶变换具体为：

S1、采集二维含噪地震数据体，选取第m道信号x_m(t)作为当前道；

S2、对选取的第m道地震信号x_m(t)做短时傅里叶变换，得到x_m(t)的短时傅里叶变换为STFT_g(x_m；t,ω)。

进一步的，步骤S2中，取实数窗函数为g(t)，对第m道地震信号x_m(t)做短时傅里叶变换，得到窗函数为g(t)时x_m(t)的短时傅里叶变换为STFT_g(x_m；t,ω)为：