EM2INV - A finite difference based algorithm for two-dimensional inversionof geoelectromagnetic data

The paper presents an efficient finite difference based 2D-inversion geaele ctromagnetic data. The special features of the algorithm are. optimal grid generation based on grid design thumb rules, finite domain boundary conditions, interpolation matrix that permits generation of response at observation poi nts different from grid points, Gaussian elimination forward matrix solver, that enables reuse of already d ecomposed coefficient matrix, super-black notion that reduces the number of blocks with unknown resistivi ties and, in turn, the size of Jacobian matrix and bi-conjugate gradient matrix solver for inverse problem which circumvents t he need of explicit Jacobian matrix computation. The algorithm is tested rigorously by setting up exercises of diverse natur e and of practical significance. The stability of the algorithm is establis hed by inverting the synthetic response corrupted with Gaussian noise. The inversion experiments are aimed at studying relative performance of response functions; inversion quality of E and B-polarization data, efficacy of single and multi-frequency data inversion, minimum number of frequencies and observation points needed for successful data inversion. It has been observed that the Magneto-telluric data deciphers better the ve rtical position of the target and Geomagnetic Depth Sounding data deciphers the horizontal variations in a better way. The conductive and resistive bo dies are better resolved by inversion of E- and B-polarization data respect ively. The results of multi-frequency inversion imply that the increase in the number of frequencies does not necessarily enhance the inversion qualit y especially when the spread of observation points is sufficiently large to sense the target. The study of a minimum number of observation points high lights the importance of single point inversion that furnishes useful infor mation about the inhomogeneity.