Impulse interference processing for MT data based on a new adaptive wavelet threshold de-noising method

Wavelet de-noising method is often used in the processing of magnetotelluric (MT) signal, and the wavelet hard and soft threshold de-noising method are the most popular although there is much room for improvement in the threshold function selection and threshold determination. A new adaptive wavelet threshold de-noising method was proposed by selecting a new threshold function and presenting an adaptive method for obtaining optimal threshold based on the multi-resolution Stein unbiased risk estimation. New threshold function and an adaptive method to determine threshold were discussed, and the principle and implementation of the algorithm were given. The simulated signal and the measured MT data contaminated by impulse interference were analyzed, and the obtained results were compared with those of the conventional wavelet hard and soft threshold de-noising methods. The results show that the proposed method overcomes the defects of the traditional wavelet soft and hard threshold due to a new threshold function, and a new method to determine the threshold of each layer is applied and provides an adaptive method for filtering MT data in the wavelet domain that requires a minimum of human intervention. The presented de-noising method is very suitable for suppressing the impulse interference for MT data and can get higher signal-to-noise ratio than the traditional wavelet threshold de-noising methods. After de-noising, the accurate data is loaded for further impedance estimation and geological interpretation.