Straightforward inversion scheme (SIS) for one-dimensional magnetotelluric data

Abstract

This paper presents a Straightforward Inversion Scheme (SIS) for interpreting one-dimensional magnetotelluric sounding data. The basic steps of SIS are (i) parameterization of the layered model such that the layer thickness, expressed in units of its skin depth, is a constant (α); (ii) expansion of the reflection function at each interface as a power series in parameter u = exp(-2(1 +j)α√f);(iii) development of a recurrence relation between the coefficients of the same powers ofu in the power series of reflection functions of any two successive layers; (iv) estimation of the impedance power series coefficients using regressed minimum norm estimator; and (v) evaluation of layer resistivities and thicknesses using the inverse recurrence relation. The power of SIS is established by inverting four synthetic data sets and two field data sets. The effect of noise is extensively studied on a synthetic data set, deliberately corrupted with increasing levels of Gaussian random noise up to 25%. It is found that the scheme can retrieve broad features of the true model even with noise levels as high as 25%. On the basis of findings of different experiments conducted on SIS, it is concluded that SIS is an efficient, robust algorithm with high resolving power. Further, being linear, it is non-iterative and it dispenses with the requirement of having to choose an initial guess model.