Articles written in Journal of Earth System Science
Volume 107 Issue 2 June 1998 pp 149-154
This paper reviews the validity of earlier models obtained after quantitative interpretation of GDS data and presents a fresh model using the inversion scheme EM2INV. The 2-D inversion of data is more objective than the earlier interpretation performed by using trial and error method. The inversion results indicate that the present model differs from the earlier ones. The reason could be that available GDS data are sufficient only for deriving the horizontal variation of subsurface resistivity. In order to study the vertical resistivity variation additional MT sounding data would be required. It would therefore be desirable to carry out MT survey in the specified area. A more comprehensive/appropriate model could be derived from joint inversion of GDS and MT data.
Volume 108 Issue 4 December 1999 pp 233-253
The paper presents an efficient finite difference based 2D-inversion algorithm, EM2INV, for geoelectromagnetic 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 points different from grid points,
• Gaussian elimination forward matrix solver, that enables reuse of already decomposed coefficient matrix,
• super-block notion that reduces the number of blocks with unknown resistivities and, in turn, the size of Jacobian matrix and
• bi-conjugate gradient matrix solver for inverse problem which circumvents the need of explicit Jacobian matrix computation.
The algorithm is tested rigorously by setting up exercises of diverse nature and of practical significance. The stability of the algorithm is established 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 vertical position of the target and Geomagnetic Depth Sounding data deciphers the horizontal variations in a better way. The conductive and resistive bodies are better resolved by inversion of E- and B-polarization data respectively. The results of multi-frequency inversion imply that the increase in the number of frequencies does not necessarily enhance the inversion quality especially when the spread of observation points is sufficiently large to sense the target. The study of a minimum number of observation points highlights the importance of single point inversion that furnishes useful information about the inhomogeneity.