Three-dimensional magnetotelluric inversion using a heterogeneous smoothness-constrained least-squares method
Tateyuki Negi 1 4 Hideki Mizunaga 2 Koichi Asamori 3 Koji Umeda 31 Nittetsu Mining Consultants Co. Ltd, 2-3, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan.
2 Department of Earth Resources Engineering, Kyushu University, 744 Motooka Nishi-ku, Fukuoka 819-0395, Japan.
3 Tono Geoscientific Research Unit, Geological Isolation Research and Development Directorate, Japan Atomic Energy Agency, 959-31 Jorinji, Izumi-cho, Toki-shi, Gifu 509-5102, Japan.
4 Corresponding author. Email: negi@nmconsults.co.jp
Exploration Geophysics 44(3) 145-155 https://doi.org/10.1071/EG13026
Submitted: 4 March 2013 Accepted: 12 March 2013 Published: 16 April 2013
Abstract
This paper presents a fast algorithm for electromagnetic data inversion to three-dimensional (3D) resistivity models. The algorithm is distinctive for the level of accuracy it attains while bypassing the sensitivity matrix update. A common sensitivity matrix for homogeneous half-space is used in all iterations. Instead of updating the sensitivity matrix, the smoothness filter coefficients at each model element are updated, based on the spatial variations in resistivity in the model derived from the latest iteration. This substitution is expected not only to reduce the computation time required for large-scale inversions, such as those for 3D surveys, but also to allow the resolution of sharp boundaries in resistivity structures. Our algorithm was applied to 3D magnetotelluric inversion in order to confirm its effectiveness. Using synthetic examples under several conditions, we demonstrated that the method can reduce the number of forward calculations required to reduce data misfits to noise level, and that the method is robust for constructing target models even with sharp boundaries without generating fatally false resistivity structures or boundaries under noisy conditions.
Key words: inversion, magnetotelluric, smoothness, three-dimensional.
References
Akaike, H., 1980, Likelihood and the Bayes procedure (with discussion), in J. M. Bernardo, M. M. De Groot, D. V. Lindley, and A. F. M. Smith, eds., Bayesian statistics: University Press (Valencia, Spain), 143–166.Avdeev, D. B., Kuvshinov, A. V., Pankratov, O. V., and Newman, G. A., 1997, High-performance three-dimensional electromagnetic modeling using modified Neumann seriesa: wide-band numerical solution and examples Journal of Geomagnetism and Geoelectricity, 49, 1519–1539
| High-performance three-dimensional electromagnetic modeling using modified Neumann seriesa: wide-band numerical solution and examplesCrossref | GoogleScholarGoogle Scholar |
Bjorck, A., 1967, Solving linear least squares problems by Gram-Schmidt orthogonalization BIT Numerical Mathematics, 7, 1–21
| Solving linear least squares problems by Gram-Schmidt orthogonalizationCrossref | GoogleScholarGoogle Scholar |
de Groot-Hedlin, C., and Constable, S. C., 1990, Occam’s inversion to generate smooth, two-dimensional models from magnetotelluric data Geophysics, 55, 1613–1624
| Occam’s inversion to generate smooth, two-dimensional models from magnetotelluric dataCrossref | GoogleScholarGoogle Scholar |
de Groot-Hedlin, C., and Constable, S. C., 2004, Inversion of magnetotelluric data for 2D structure with sharp resistivity contrasts Geophysics, 69, 78–86
| Inversion of magnetotelluric data for 2D structure with sharp resistivity contrastsCrossref | GoogleScholarGoogle Scholar |
Fujino, S., 2002, GPBiCG (m, ℓ): a hybrid of BiCGSTAB and GPBiCG methods with efficiency and robustness Applied Numerical Mathematics, 41, 107–117
| GPBiCG (m, ℓ): a hybrid of BiCGSTAB and GPBiCG methods with efficiency and robustnessCrossref | GoogleScholarGoogle Scholar |
Kimura, T., Goto, T., Kasaya, T., Okamoto, T., Mikada, H., Sanada, Y., Watanabe, T., and Ashida, Y., 2010, Two-dimensional inversion of magnetotelluric data with sharp structural boundary Butsuri Tansa, 63, 185–196
Loke, M. H., and Barker, R. D., 1996, Rapid least-squares inversion of apparent resistivity pseudosections by a quasi-Newton method Geophysical Prospecting, 44, 131–152
| Rapid least-squares inversion of apparent resistivity pseudosections by a quasi-Newton methodCrossref | GoogleScholarGoogle Scholar |
Mackie, R. L., Smith, J. T., and Madden, T. R., 1994, Three-dimensional electromagnetic modeling using finite difference equations: the magnetotelluric example Radio Science, 29, 923–935
| Three-dimensional electromagnetic modeling using finite difference equations: the magnetotelluric exampleCrossref | GoogleScholarGoogle Scholar |
Meijerink, J. A., and Van der Vorst, H. A., 1981, Guidelines for the usage of incomplete decompositions in the solving sets of linear equations as they occur in practical problems Journal of Computational Physics, 44, 134–155
| Guidelines for the usage of incomplete decompositions in the solving sets of linear equations as they occur in practical problemsCrossref | GoogleScholarGoogle Scholar |
Ogawa, Y., and Uchida, T., 1996, A two-dimensional magnetotelluric inversion assuming Gaussian static shift Geophysical Journal International, 126, 69–76
| A two-dimensional magnetotelluric inversion assuming Gaussian static shiftCrossref | GoogleScholarGoogle Scholar |
Sasaki, Y., 1989, Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data Geophysics, 54, 254–262
| Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity dataCrossref | GoogleScholarGoogle Scholar |
Sasaki, Y., 1994, 3-D resistivity inversion using the finite element method Geophysics, 59, 1839–1848
| 3-D resistivity inversion using the finite element methodCrossref | GoogleScholarGoogle Scholar |
Sasaki, Y., 2001, Full 3-D inversion of electromagnetic data on PC Journal of Applied Geophysics, 46, 45–54
| Full 3-D inversion of electromagnetic data on PCCrossref | GoogleScholarGoogle Scholar |
Sasaki, Y., 2004, Three-dimensional inversion of static-shifted magnetotelluric data Earth, Planets, and Space, 56, 239–248
Siripunvaraporn, W., and Egbert, G., 2009, WSINV3DMT: vertical magnetic field transfer function inversion and parallel implementation Physics of the Earth and Planetary Interiors, 173, 317–329
| WSINV3DMT: vertical magnetic field transfer function inversion and parallel implementationCrossref | GoogleScholarGoogle Scholar |
Siripunvaraporn, W., Egbert, G., Lenbury, Y., and Uyeshima, M., 2005, Three-dimensional magnetotelluric inversion: data space method Physics of the Earth and Planetary Interiors, 150, 3–14
| Three-dimensional magnetotelluric inversion: data space methodCrossref | GoogleScholarGoogle Scholar |
Smith, J. T., 1996, Conservative modeling of 3-D electromagnetic fields, Part II: biconjugate gradient solution and an accelerator Geophysics, 61, 1319–1324
| Conservative modeling of 3-D electromagnetic fields, Part II: biconjugate gradient solution and an acceleratorCrossref | GoogleScholarGoogle Scholar |
Smith, T., Hoversten, M., Gasperikova, E., and Morrison, F., 1999, Sharp boundary inversion of 2D magnetotelluric data Geophysical Prospecting, 47, 469–486
| Sharp boundary inversion of 2D magnetotelluric dataCrossref | GoogleScholarGoogle Scholar |
Uchida, T., 1993, Smoothness-constrained 2D inversion for DC resistivity data by ABIC minimization method Butsuri Tansa, 46, 105–119
Uchida, T., and Sasaki, Y., 2006, Stable 3D inversion of MT data and its application to geothermal exploration Exploration Geophysics, 37, 223–230
| Stable 3D inversion of MT data and its application to geothermal explorationCrossref | GoogleScholarGoogle Scholar |
Wannamaker, P. E., 1991, Advances in three-dimensional magnetotelluric modeling using integral equations Geophysics, 56, 1716–1728
| Advances in three-dimensional magnetotelluric modeling using integral equationsCrossref | GoogleScholarGoogle Scholar |