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

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.

Ways to sparse representation:An overview

Many algorithms have been proposed to find sparse representations over redundant dictionaries or transforms. This paper gives an overview of these algorithms by classifying them into three categories： greedy pursuit algorithms, lp norm regularization based algorithms, and iterative shrinkage algorithms. We summarize their pros and cons as well as their connections. Based on recent evidence, we conclude that the algorithms of the three categories share the same root： lp norm regularized inverse problem. Finally, several topics that deserve further investigation are also discussed.

Ways to Sparse Representation: A Comparative Study

Many algorithms have been proposed to achieve sparse representation over redundant dictionaries or transforms. A comprehensive understanding of these algorithms is needed when choosing and designing algorithms for particular applications. This research studies a representative algorithm for each category, matching pursuit （MP）, basis pursuit （BP）, and noise shaping （NS）, in terms of their sparsifying capability and computational complexity. Experiments show that NS has the best performance in terms of sparsifying ca- pability with the least computational complexity. BP has good sparsifying capability, but is computationally expensive. MP has relatively poor sparsifying capability and the computations are heavily dependent on the problem scale and signal complexity. Their performance differences are also evaluated for three typical ap- plications of time-frequency analyses, signal denoising, and image coding. NS has good performance for time-frequency analyses and image coding with far fewer computations. However, NS does not perform well for signal denoising. This study provides guidelines for choosing an algorithm for a given problem and for designing or improving algorithms for sparse representation.