Power-line interference suppression of MT data based on frequency domain sparse decomposition  

基于频域稀疏分解的大地电磁工频干扰压制（英文）

作　　者：TANG Jing-tian LI Guang ZHOU Cong LI Jin LIU Xiao-qiong ZHU Hui-jie 汤井田;李广;周聪;李晋;刘晓琼;朱会杰

机构地区：[1]School of Geosciences and Info-Physics,Central South University,Changsha 410083,China [2]Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring,Ministry of Education(Central South University),Changsha 410083,China [3]College of Information Science and Engineering,Hunan Normal University,Changsha 410081,China [4]The First Engineering Scientific Research Institute of General Armaments Department,Wuxi 214035,China

出　　处：《Journal of Central South University》2018年第9期2150-2163,共14页中南大学学报（英文版）

基　　金：Project(2014AA06A602)supported by the National High-Tech Research and Development Program of China;Projects(41404111,41304098)supported by the National Natural Science Foundation of China;Project(2015JJ3088)supported by the Natural Science Foundation of Hunan Province,China

摘　　要：Power-line interference is one of the most common noises in magnetotelluric(MT)data.It usually causes distortion at the fundamental frequency and its odd harmonics,and may also affect other frequency bands.Although trap circuits are designed to suppress such noise in most of the modern acquisition devices,strong interferences are still found in MT data,and the power-line interference will fluctuate with the changing of load current.The fixed trap circuits often fail to deal with it.This paper proposes an alternative scheme for power-line interference removal based on frequency-domain sparse decomposition.Firstly,the fast Fourier transform of the acquired MT signal is performed.Subsequently,a redundant dictionary is designed to match with the power-line interference which is insensitive to the useful signal.Power-line interference is separated by using the dictionary and a signal reconstruction algorithm of compressive sensing called improved orthogonal matching pursuit(IOMP).Finally,the frequency domain data are switched back to the time domain by the inverse fast Fourier transform.Simulation experiments and real data examples from Lu-Zong ore district illustrate that this scheme can effectively suppress the power-line interference and significantly improve data quality.Compared with time domain sparse decomposition,this scheme takes less time consumption and acquires better results.工频干扰是指由电网产生的基频及其奇次谐波干扰,它是大地电磁信号采集过程中最为普遍的干扰之一。尽管大部分采集设备都设计有抑制工频干扰的陷波电路,但由于电网中电流的实际频率会随着负载的变化有所波动,而陷波器的中心频率是固定的,因此在实际采集时,大地电磁信号依然受到工频噪声的严重影响。实践经验表明,当受工频干扰影响时,远参考法时常难以奏效;工频干扰奇次谐波的幅值随着频率的增大而骤减,在时间域难以准确识别,因此时间域编辑法效果不佳;此外,由于干扰源是固定的,工频干扰通常存在于整个采集过程中,通过数据段筛选也无法去除噪声。本文基于频域稀疏分解,先对采集的大地电磁信号进行傅里叶变换。然后设计与干扰信号相匹配而对有用信号不敏感的冗余字典原子,结合IOMP算法分离出频域信号中的工频干扰成分。最后将处理后的频域信号进行傅里叶逆变换。仿真实验及案例分析表明,所述方法能够在较好地保留有用信号的前提下有效压制工频干扰,在大大降低时间消耗的基础上取得比时域稀疏分解更好的效果,改善大地电磁数据质量。

关 键 词：sparse representation magnetotelluric signal processing power-line noise improved orthogonal matching pursuit redundant dictionary 

分 类 号：P631.325[天文地球—地质矿产勘探]