基于独立分量分析工频干扰消除技术.docVIP

基于独立分量分析工频干扰消除技术 　　（中国地质大学 地球物理与信息技术学院,北京 100083） 　　摘 要：简要介绍了ICA的基本原理和快速算法，在分析地震信号和工频干扰特点的基础上，利用ICA技术来消除地震记录中的工频干扰，并与常规方法进行比较。研究结果表明ICA在有效消除工频干扰的同时，能够保护有效信号，并且在提高资料的信噪比方面更有优势，具有良好的应用前景。 　　关键词：独立分量分析；盲源分离；快速独立分量分析算法；工频干扰 　　中图分类号：TP391 文献标志码：A 　　 文章编号：1001?B3695(2009)01?B0227?B03 　　 　　Power interference removal method based on independent component analysis 　　WEI Wei,LIU Xuewei 　　（ School of Geophysics Information Technology， China University of Geosciences， Beijing 100083， China） 　　Abstract:This paper studied the fundamentals and fast algorithm of ICA, tried to remove the power interference in seismic data on the base of analyzing the features of seismic data and the power interference. Compared with the results of the conventional methods, ICA has good effect on removing the power interference and can protect the available signals effectively. Results show that ICA has better effect on improving the signal to noise ratio and has a good application perspective. 　　Key words:independent component analysis（ICA）; blind source separation（BSS）; FastICA algorithm; power interference 　　在地震勘探数据采集时，如果有工业电传输线穿过勘探区域，则在采集的地震数据中将包含工频干扰，这势必会影响地震勘探的精度。目前压制工频干扰主要采用两种方法：a）陷波滤波法。它是在频率域中实现的，设计多吸收点的、具有一定阻带宽度的陷波滤波器组，用于消除工频干扰的基波和谐波成分。如果信号频谱与工频干扰的频谱有混叠，则陷波滤波器在滤除工频干扰的同时也会造成目标信号的损失。b）自适应滤波法。该方法的前提是目标信号和干扰波不相关。在此基础上，利用自适应滤波算法(最小均方或最小二乘法)，自动调整滤波器系数，以跟踪输入过程的变化，并实现工频干扰的自适应抵消。这种方法在时间域进行，其本质与陷波滤波是一致的[1]。?? 　　ICA是近期发展起来的一种非常有效的盲源分离技术。其处理的对象是一组相互统计独立的信源经线性组合而产生的混合信号，最终从混合信号中提取出各独立的信号分量。目前ICA 在生物医学信号处理、混合语音信号分离、图像消噪等方面已取得了较好的应用效果[2~5]。本文在对ICA的基本原理和快速算法的研究基础上，分析了地震数据和工频干扰的特点，深入研究了基于ICA的工频干扰消除技术，并与常规方法进行比较，结果表明ICA技术可以有效消除工频干扰，并且不伤害有效信号在提高资料的信噪比方面更有优势，结果令人满意。?? 　　 1 ICA基本原理?? 　　下面以两个信号的盲源分离为例介绍ICA的基本原理。在图1中,s1(t)、s2(t)是两个非高斯分布的独立源；x1(t)、x2(t)是两个观察信号, 混合矩阵A是未知的模型系数。?? 　　x1(t)x2(t)=a11 a12 　　a21 a22×s1(t)s2(t)(1) 　　或用矩阵形式描述为?? 　　X=A×S（2） 　　在独立源信号S和混合矩阵A都是未知的情况下，希望能寻找到一个分离矩阵W=A-1从观测信号中进行源信号的分离,即?? 　　y=A-1×X=W×X=S(3) 　　其中：X、S分别为观测信号和源信号向量；A为混合矩阵。?? 　　ICA方法可以有效解决上述盲源分离问题，其基本思想就是中心极限定理