基于递归滤波的工频噪声压制方法研究.pdfVIP

第 51卷第 4期 石 油 物 探 Vo1．51，No．4 2012年 7月 GEOPHYSICALPROSPECTING F0Rl PETR0LEUM Ju1．，2012 文章编号：1000—1441(2012)04—0362—05 基于递归滤波的工频噪声压制方法研究 康怡亭 ，王 鑫 ，赵彦青 (1．中国石油化工股份有限公司胜利油田分公司，山东东营257000；2．上海锦迪软件开发有限公 司，上海 200092) 摘要：在老油区地震勘探中，工频噪声干扰严重影响了地震资料的品质。压制工频噪声的常规方法是陷波处理， 该方法在压制噪声的同时对有效信号损伤很大。另一类基于正弦或余弦函数逼近的工频噪声压制方法，通常假 设工频噪声的振幅、频率和相位保持固定不变，在实际资料处理时其适用性受到很大限制 ，并且计算量较大，算 法效率不高。在分析已有方法的基础上，提出了一种基于递归滤波的工频噪声压制方法，利用工频噪声的周期 性 ，导出递归滤波器，并对递归滤波算子作适当修改 ，保证了滤波器的稳定性和物理可实现性。理论记录试算与 实际资料处理结果表明，该方法能够有效压制工频噪声，且不损伤有效信号。 关键词 ：工频噪声；递归滤波；噪声压制 DOI：10．3969／j．issn．1000—1441．2012．04．007 中图分类号：P631．4 文献标识码：A 老油区工业化程度高、高压线密集，采集的地 量较大，算法效率不高。近年来，还有一种新兴的 震记录中常常存在强工频噪声，极大地影响了地震 压噪方法——独立分量分析 (ICA)方法[]，该方 资料的品质。这种强工频干扰规律性较强，主频约 法最大的问题是噪声分离后尺度与顺序不确定，无 50Hz，是典型的加性噪声，常与有效信号叠合在 法获得独立信号的真实幅值[1，从而无法完全去 一 起，在很大程度上影响了地震资料的叠加和偏移 除噪声。 成像的质量_1]。 针对工频噪声的特点，且考虑到现有方法的缺 目前压制工频噪声的方法有很多，最常见的是 陷，我们提出了一种基于递归滤波的工频噪声压制 陷波滤波法。由于在频率域工频噪声与有效信号 方法，可以充分适应不同时间段工频噪声的变化， 相叠合，陷波滤波在去除工频噪声的同时也损伤了 以提高去噪的精度。 同频率段的有效信号 ]。凌云等[4最早提出了利 用频率域中值滤波压制工频噪声的方法。该方法 1 方法原理 能够压制工频噪声主要能量，但从频谱上看，工频 噪声依然存在。李文艳等[5]提出了基于互相关原 对于具有一定频率的周期性噪声，设第 1个周 理的50Hz干扰波剔除方法，但重构的噪声主频 期内噪声信号为 ，l()，噪声的周期为T，经过周期 为50Hz，与实际工频噪声主频(49~51Hz)不符。 T时间延迟后，噪声为 厂2(￡)，经过 个周期后噪声 高少武等[6]分析了频率域滤波的缺陷，并提出基于 为 (￡)，即 余弦函数的自适应方法，利用余弦函数逼近 50Hz fl( T ㈩ 噪声，计算其振幅、频率及初始相位，达到压制与去 除噪声的目的。刘洋口]提出了应用共轭梯度法估 f0 0 t≤ T 算工频噪声的振幅、频率及初始相位，利用正弦波 ()=== fl(￡一T) T t≤ 2T (2) 逼近噪声从而消除噪声。胡伟等[8]提出了改进的 l0 2T 余弦逼近法，采用精细扫描、标准余弦信号超采样 ；