Magnetotelluric (MT) sounding data can easily be damaged by various types of noise, especially in industrial areas, where the quality of measured data is poor. Most traditional de-noising methods are ineffective to the low signal-to-noise ratio of data. To solve the above problem, we propose the use of a de-noising method for the detection of noise in MT data based on discrete wavelet transform and singular value decomposition (SVD), with multiscale dispersion entropy and phase space reconstruction carried out for pretreatment. No “over processing” takes place in the proposed method. Compared with wavelet transform and SVD decomposition in synthetic tests, the proposed method removes the profile of noise more completely, including large-scale noise and impulse noise. For high levels or low levels of noise, the proposed method can increase the signal-to-noise ratio of data more obviously. Moreover, application to the field MT data can prove the performance of the proposed method. The proposed method is a feasible method for the elimination of various noise types and can improve MT data with high noise levels, obtaining a recovery in the response. It can improve abrupt points and distortion in MT response curves more effectively than the robust method can.

We have developed a mineral exploration method for the joint inversion of 2D gravity gradiometry and magnetotelluric (MT) data based on data-space and normalized cross-gradient constraints. To accurately explore the underground structure of complex mineral deposits and solve the problems such as the non-uniqueness and inconsistency of the single parameter inversion model, it is now common practice to perform collocated MT and gravity surveys that complement each other in the search. Although conventional joint inversion of MT and gravity using model-space can be diagnostic, we posit that better results can be derived from the joint inversion of the MT and gravity gradiometry data using data-space. Gravity gradiometry data contains more abundant component information than traditional gravity data and can be used to classify the spatial structure and location of underground structures and field sources more accurately and finely, and the data-space method consumes less memory and has a shorter computation time for our particular inversion iteration algorithm. We verify our proposed method with synthetic models. The experimental results prove that our proposed method leads to models with remarkable structural resemblance and improved estimates of electrical resistivity and density and requires shorter computation time and less memory. We also apply the method to field data to test its potential use for subsurface lithofacies discrimination or structural classification. Our results suggest that the imaging method leads to improved characterization of geological targets, which is more conducive to geological interpretation and the exploration of mineral resources.

We present a joint 2D inversion approach for magnetotelluric (MT) and gravity data with elastic-net regularization and cross-gradient constraints. We describe the main features of the approach and verify the inversion results against a synthetic model. The results indicate that the best fit solution using the L2 is overly smooth, while the best fit solution for the L1 norm is too sparse. However, the elastic-net regularization method, a convex combination term of L2 norm and L1 norm, can not only enforce the stability to preserve local smoothness, but can also enforce the sparsity to preserve sharp boundaries. Cross-gradient constraints lead to models with close structural resemblance and improve the estimates of the resistivity and density of the synthetic dataset. We apply the novel approach to field datasets from a copper mining area in the northeast of China. Our results show that the method can generate much more detail and a sharper boundary as well as better depth resolution. Relative to the existing solution, the large area divergence phenomenon under the anomalous bodies is eliminated, and the fine anomalous bodies boundary appeared in the smooth region. This method can provide important technical support for detecting deep concealed deposits.