Three-dimensional inversion of CSAMT data in the presence of topography
Changhong Lin 1 2 Handong Tan 1 Wangwang Wang 1 Tuo Tong 1 Miao Peng 1 Mao Wang 1 Weihua Zeng 11 School of Geophysics and Information Technology, China University of Geosciences (Beijing), Beijing 100083, China.
2 Corresponding author. Email: linchh@cugb.edu.cn
Exploration Geophysics 49(3) 253-267 https://doi.org/10.1071/EG16067
Submitted: 15 June 2016 Accepted: 6 April 2017 Published: 23 May 2017
Abstract
3D controlled-source audio frequency magnetotelluric (CSAMT) responses can be distorted strongly by topography and should be accounted for in data inversion and interpretation. In this paper we present a scheme to incorporate topographic distortions into the inversion instead of correcting them. This approach has been verified by comparing the modelling results with 2D FEM CSAMT solutions and synthetic inversion examples. Compared with the responses generated from a half-space model with flat surface, it is found that not only the topography in the survey area but also that at the source position may strongly distort the CSAMT responses. The field example indicates that results with topography are much better than those without considering topography to map the distribution of coal seam underground, which also illustrates the effectiveness of our approach.
Key words: 3D inversion, CSAMT, topography.
References
Chave, A. D., 1983, Numerical integration of related Hankel transforms by quadrature and continued fraction expansion: Geophysics, 48, 1671–1686| Numerical integration of related Hankel transforms by quadrature and continued fraction expansion:Crossref | GoogleScholarGoogle Scholar |
Commer, M., and Newman, G., 2008, New advances in three-dimensional controlled-source electromagnetic inversion: Geophysical Journal International, 172, 513–535
| New advances in three-dimensional controlled-source electromagnetic inversion:Crossref | GoogleScholarGoogle Scholar |
Key, K., and Ovall, J., 2011, A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic modeling: Geophysical Journal International, 186, 137–154
| A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic modeling:Crossref | GoogleScholarGoogle Scholar |
Li, S. H., Booker, J. R., and Aprea, C., 2008, Inversion of magnetotelluric data in the presence of strong bathymetry/topography: Geophysical Prospecting, 56, 259–268
| Inversion of magnetotelluric data in the presence of strong bathymetry/topography:Crossref | GoogleScholarGoogle Scholar |
Lin, C. H., Tan, H. D., and Tong, T., 2011, Three-dimensional conjugate gradient inversion of magnetotelluric full information data: Applied Geophysics, 8, 1–10
| Three-dimensional conjugate gradient inversion of magnetotelluric full information data:Crossref | GoogleScholarGoogle Scholar |
Lu, X. Y., Unsworth, M., and Booker, J., 1999, Rapid relaxation inversion of CSAMT data: Geophysical Journal International, 138, 381–392
| Rapid relaxation inversion of CSAMT data:Crossref | GoogleScholarGoogle Scholar |
Mackie, R. L., and Madden, T. R., 1993, Three-dimensional magnetotelluric inversion using conjugate gradients: Geophysical Journal International, 115, 215–229
| Three-dimensional magnetotelluric inversion using conjugate gradients:Crossref | GoogleScholarGoogle Scholar |
Mackie, R., Watts, M. D., and Rodi, W., 2007, Joint 3D inversion of marine CSEM and MT data: 77th Annual International Meeting, SEG, Expanded Abstracts, 574–578.
Madden, T. R., and Mackie, R. L., 1989, Three-dimensional magnetotelluric modeling and inversion: Proceedings of the IEEE, 77, 318–333
| Three-dimensional magnetotelluric modeling and inversion:Crossref | GoogleScholarGoogle Scholar |
Nam, M. J., Kim, H. J., Song, Y., Lee, T. J., Son, J. S., and Suh, J. H., 2007, 3D magnetotelluric modeling including surface topography: Geophysical Prospecting, 55, 277–287
| 3D magnetotelluric modeling including surface topography:Crossref | GoogleScholarGoogle Scholar |
Newman, G. A., and Alumbaugh, D. L., 1995, Frequency-domain modeling of airborne electromagnetic responses using staggered finite differences: Geophysical Prospecting, 43, 1021–1042
| Frequency-domain modeling of airborne electromagnetic responses using staggered finite differences:Crossref | GoogleScholarGoogle Scholar |
Newman, G. A., and Alumbaugh, D. L., 2000, Three-dimensional magnetotelluric inversion using non-linear conjugate gradients: Geophysical Journal International, 140, 410–424
| Three-dimensional magnetotelluric inversion using non-linear conjugate gradients:Crossref | GoogleScholarGoogle Scholar |
Newman, G. A., and Boggs, P. T., 2004, Solution accelerators for large-scale three-dimensional electromagnetic inverse problems: Inverse Problems, 20, S151–S170
| Solution accelerators for large-scale three-dimensional electromagnetic inverse problems:Crossref | GoogleScholarGoogle Scholar |
Nocedal, J., and Wright, S. J., 2006, Numerical optimization: Springer-Verlag.
Penz, S., Chauris, H., Donno, D., and Mehl, C., 2013, Resistivity modeling with topography: Geophysical Journal International, 194, 1486–1497
| Resistivity modeling with topography:Crossref | GoogleScholarGoogle Scholar |
Rodi, W., and Mackie, R. L., 2001, Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion: Geophysics, 66, 174–187
| Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion:Crossref | GoogleScholarGoogle Scholar |
Schwarzbach, C., Börner, R.-U., and Spitzer, K., 2011, Three-dimensional adaptive higher order finite element simulation for geo-electromagnetics – a marine CSEM example: Geophysical Journal International, 187, 63–74
| Three-dimensional adaptive higher order finite element simulation for geo-electromagnetics – a marine CSEM example:Crossref | GoogleScholarGoogle Scholar |
Smith, J. T., 1996, Conservative modeling of 3D electromagnetic fields part II: biconjugate gradient solution and an accelerator: Geophysics, 61, 1319–1324
| Conservative modeling of 3D electromagnetic fields part II: biconjugate gradient solution and an accelerator:Crossref | GoogleScholarGoogle Scholar |
Streich, R., 2009, 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: direct solution and optimization for high accuracy: Geophysics, 74, F95–F105
| 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: direct solution and optimization for high accuracy:Crossref | GoogleScholarGoogle Scholar |
Tan, H. D., Yu, Q. F., Booker, J., and Wei, W. B., 2003, Three-dimensional magnetotelluric modeling using the staggered-grid finite difference method: Chinese Journal of Geophysics, 46, 1011–1020
| Three-dimensional magnetotelluric modeling using the staggered-grid finite difference method:Crossref | GoogleScholarGoogle Scholar |
Tikhonov, A. N., and Arsenin, V. Y., 1977, Solutions of ill-posed problems: V. H. Winston and Sons.
Unsworth, M. J., Bryan, J. T., and Alan, D. C., 1993, Electromagnetic induction by a finite electric dipole source over a 2-D earth: Geophysics, 58, 198–214
| Electromagnetic induction by a finite electric dipole source over a 2-D earth:Crossref | GoogleScholarGoogle Scholar |
Wannamaker, P. E., Stodt, J. A., and Rijo, L., 1986, Two-dimensional topographic responses in magnetotellurics modeled using finite elements: Geophysics, 51, 2131–2144
| Two-dimensional topographic responses in magnetotellurics modeled using finite elements:Crossref | GoogleScholarGoogle Scholar |
Ward, S. H., and Hohmann, G. W., 1988, Geophysical electromagnetic theory, in M. N. Nabighian, ed., Electromagnetic methods in applied geophysics, volume 1, Theory: Society of Exploration Geophysicists.
Wu, X. P., and Xu, G. M., 2000, Study on 3-D resistivity inversion using conjugate gradient method: Chinese Journal of Geophysics, 43, 450–458
| Study on 3-D resistivity inversion using conjugate gradient method:Crossref | GoogleScholarGoogle Scholar |
Zhdanov, M. S., Lee, S. K., and Yoshioka, K., 2006, Integral equation method for 3D modeling of electromagnetic fields in complex structures within homogeneous background conductivity: Geophysics, 71, G333–G345
| Integral equation method for 3D modeling of electromagnetic fields in complex structures within homogeneous background conductivity:Crossref | GoogleScholarGoogle Scholar |