[发明专利]基于波形形态特征稀疏化建模的谐波噪声压制方法

摘要

本发明公开一种基于波形形态特征稀疏化建模的谐波噪声压制方法，包括：1)、根据可控震源滑动扫描采集地震记录中谐波噪声波形形态特征构造Chirplet变换，并与连续小波变换构成超完备字典；2)、Chirplet正、反变换的快速实现；3)、基于相关后数据时频分布特征确定Chirplet变换参数；4)、根据参考扫描信号的起始频率确定谐波噪声的滤波截止频率，实现有效信号与谐波噪声保真分离。本发明解决了可控震源滑动扫描采集地震数据中谐波噪声干扰的问题，达到了提高了地震资料信噪比的目的。本发明方法根据谐波噪声的时频分布特征确定选择Chirplet变换，保证了稀疏性；根据数据驱动自动确定Chirplet变换的参数，具有很强的适应性，且单道计算，便于并行处理。

主权项

1.基于波形形态特征稀疏化建模的谐波噪声压制方法，其特征在于，包括以下步骤：1)、根据可控震源滑动扫描采集地震记录中谐波噪声的波形形态特征构造Chirplet变换，并与连续小波变换构成超完备字典；2)、Chirplet正、反变换的快速实现；3)、基于相关后数据时频分布特征确定Chirplet变换参数；4)、根据参考扫描信号的起始频率确定谐波噪声滤波截止频率，实现有效信号与谐波噪声保真分离；步骤2)具体包括：步骤101：从原始数据中读出单道时域信号，记为x；步骤102：计算得到信号x的Chirplet正变换的结果，记为y；步骤103：计算得到y的Chirplet反变换结果,记为步骤201：构造一个线性升频信号模型，记为s，其长度与信号x相同；步骤202：计算得到信号s的Chirplet正变换结果，记为s0；步骤203：计算得到s₀的Chirplet反变换结果，记为步骤204：根据模型反变换结果计算修正系数，即其中||.||₂表示二范数；步骤104：将修正系数r作用到单道时域信号x的Chirplet反变换得到修正的Chirplet反变换结果，记为则