2D joint inversion of CSAMT and magnetic data based on cross-gradient theory
Wang Kun-Peng1, Tan Han-Dong1,2, and Wang Tao3
1. School of Geophysics and information Technology, China University of Geosciences(Beijing), Beijing 100083, China.
2. Key Laboratory of Geo-detection (China University of Geosciences), Ministry of Education, Beijing 100083, China.
3. Key Laboratory of Electromagnetic Radiation and Sensing Technology, Chinese Academy of Sciences, Beijing 100190, China.
Abstract：
A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric (CSAMT) method is developed to invert data in the entire region (near, transition, and far) and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.
. 2D joint inversion of CSAMT and magnetic data based on cross-gradient theory[J]. APPLIED GEOPHYSICS, 2017, 14(2): 279-290.
[1]
Colombo, D., and Stefano, M. D., 2007, Geophysical modeling via simultaneous joint inversion of seismic, gravity, and electromagnetic data: Application to prestack depth imaging: The Leading Edge, 26(3), 326−331.
[2]
Constable, S. C., Parker, R. L., and Constable, C. G., 1987, Occam’s inversion: a practical algorithm for generating smooth models from electromagnetic sounding data: Geophysics, 52(3), 289−300.
[3]
DeGroot-Hedlin, C., and Constable, S., 1990, Occam’s inversion to generate smooth, two-dimensional models from magnetotelluric data: Geophysics, 55(12), 1613−1624.
[4]
Fregoso E. and Gallardo L. A., 2009, Cross-gradients joint 3D inversion with applications to gravity and magnetic data: Geophysics, 74(4), L31−L42.
[5]
Gallardo, L. A., and Meju, M. A., 2003, Characterization of heterogeneous near-surface materials by joint 2D inversion of DC resistivity and seismic data: Geophysical Research Letters, 30(13), 1658−1661.
[6]
Guan, Z. N., 2005, Geomagnetic Field and Magnetic Exploration: Geological Publishing House, Beijing, 107−110.
[7]
Jegen, M. D., Hobbs, R. W., Tarits, P., et al, 2009, Joint inversion of marine magnetotelluric and gravity data incorporating seismic constraints: Preliminary results of sub-basalt imaging off the Faroe Shelf: Earth and Planetary Science Letters, 282(1), 47−55.
[8]
Key, K., 2016, MARE2DEM: a 2-D inversion code for controlled-source electromagnetic and magnetotelluric data: Geophysical Journal International, 207(1), ggw290.
[9]
Li, Y. G., and Oldenburg, D. W., 1996, 3-D inversion of magnetic data: Geophysics, 61(2), 394−408.
[10]
Li, Y. G., and Oldenburg, D. W., 2003, Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method: Geophys. J. Int., 152, 251−265.
[11]
Lin, C. H., Tan, H. D., Shu Q., et al., 2012, Three-dimensional conjugate gradient inversion of CSAMT data: Chinese J. Geophys. (in Chinese), 55(11), 3829−3838.
[12]
McGillivray, P. R., Oldenburg, D. W., Ellis, R. G., et al., 1994, Calculation of sensitivities for the frequency-domain electromagnetic problem: Geophys. J. Int., 116, 1−4.
[13]
Mitsuhata, Y., Uchida, T., and Amano, H., 2002, 2.5-D inversion of frequency-domain electromagnetic data generated by a grounded-wire source: Geophysics, 67(6), 1753−1768.
[14]
Peng, M., 2012, Joint Inversion of Magnetotelluric and Teleseismic data: PhD Thesis, China University of Geosciences, Beijing.
[15]
Rodi, W., and Mackie, R. L., 2001, Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion: Geophysics, 66(1), 174−187.
[16]
Sasaki, Y, 1989, Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data: Geophysics, 54(2), 254−262
[17]
Siripunvaraporn, W., and Egbert, G., 2000, An efficient data-subspace inversion method for 2-D magnetotelluric data: Geophysics, 65(3), 791−803.
[18]
Vozoff, K., and Jupp, D. L. B., 1977, Effective search for a buried layer: An approach to experimental design in geophysics: Exploration Geophysics, 8(1), 6−15.
[19]
Yang, B., 2012, Three Dimensional Marine Controlled Source Electromagnetic Data Forward Modeling and Inversion with Topography: PhD Thesis, China University of Geosciences.
[20]
Zhang, B., 2012, Research of 2D CSAMT Forward and Inversion: MSc. Thesis, China University of Geosciences, Beijing.
[21]
Zhou, L. F., 2012, Two Dimensional Joint Inversion of MT and Seismic data: MSc. Thesis, China University of Geosciences, Beijing.
2017, Vol. 14 
