Magnetotelluric Impedance Estimation Method

A magnetotelluric impedance estimation method is provided, belonging to the field of exploration. The method specifically includes the following steps: Step 1, performing Fourier transform on electromagnetic time series data of a plurality of electromagnetic field components to obtain a frequency spectrum array corresponding to each of electromagnetic field components; Step 2, identifying an interference frequency of each of electromagnetic field components according to the frequency spectrum array; Step 3, merging the interference frequency and a deviation value; Step 4, attenuating the frequency spectrum array using the merged interference frequency and deviation value; Step 5, recovering the electromagnetic time series data using the attenuated frequency spectrum array, so as to obtain a transform result; and Step 6, performing magnetotelluric impedance estimation according to the transform result. By the scheme disclosed, the adaptability and accuracy of magnetotelluric impedance estimation are improved.