[发明专利]一种高精度自动识别和压制地震数据单频干扰的方法

摘要

本发明提供一种高精度自动识别和压制地震数据单频干扰的方法，包括输入原始地震数据，对深部地震信号进行求导运算得到地震信号导数的频谱，并对其进行叠加处理；根据定义的单频干扰判别准则在频谱上识别单频干扰，并将其频率作为初始频率；重构单频干扰信号，对构建的目标函数求导得到单频干扰振幅；计算由目标函数转换的归一化互相关函数值，将互相关函数的极大值对应的频率确定为单频干扰的最佳频率；将单频干扰的频率及其振幅代入正余弦函数公式，得到重构的单频干扰信号,从原始地震数据中减去重构的单频干扰信号，完成单频干扰的压制工作。本发明提供的方法能够自动识别和压制地震数据中的单频干扰，并能取得保真去噪的效果。

主权项

1.一种高精度自动识别和压制地震数据单频干扰的方法，其特征在于，包括如下步骤：S1，输入原始地震数据，选取2s以下数据进行n阶求导，利用傅氏变换将求导信号转换到频率域，得到原来单频干扰振幅4倍以上的频谱；S2，对得到的频谱进行叠加处理，消除锯齿状振幅起伏现象对单频干扰识别的影响；S3，根据单频干扰的凸起频率振幅特征定义单频干扰判别准则，根据判别准则在频谱上一次性识别出所有单频干扰，并将单频干扰对应的频率作为后续计算的初始频率；S4，利用正余弦函数逼近法重构单频干扰信号，将原始地震信号与重构单频干扰信号差的平方求和定义为目标函数，对目标函数求导得到单频干扰的振幅；S5，将定义的目标函数转换为原始地震信号与重构单频干扰的归一化互相关函数，采用频率折半逼近算法进行频率扫描，得到多个互相关函数值，将其极大值对应的频率视为单频干扰的最佳频率；S6，将单频干扰频率及用其计算的振幅代入正余弦函数公式，得到重构的单频干扰信号；S7，从原始地震数据中减去重构的单频干扰信号，获得无单频干扰的地震数据，完成消除单频干扰处理工作。