3D controlled source electromagnetic inversion with topography in the frequency domain
Zhu Cheng1, Li Tonglin2, Yang Haibin3, Liu Yongliang2, Wu Liang4, Hu Yingcai5
1. Minmetals Exploration and Development Co. LTD, Beijing 100010, China;
2. School of Geo-Exploration Science & Technology, Jilin University, Changchun, Jilin 130000, China;
3. Department of Geosciences and Natural Resource Management, University of Copenhagen, Denmark, DK1350, Denmark;
4. Changjiang Engineering Geophysical Prospecting Company, Wuhan, Hubei 430010, China;
5. Beijing Research Institute of Uranium Geology, Beijing 100029, China
Abstract:Seismic exploration is sometimes unsuccessful in some areas with difficult geological and topographical conditions such as complex mountainous area. Experiments show that the joint application with seismic and non-seismic methods especially electromagnetic (EM) is an effective way to ensure high quality prospecting. This paper discusses the theory and algorithm of three dimensional forward modeling and inversion of controlled-source electromagnetic method (CSEM) with topography, which is tested to be accurate in theoretical data inversion. The paper consists of two parts:3D forward modeling and non-linear conjugate gradient (NLCG) inversion. For 3D forward modeling, the finite difference based on scattered fields is implemented, and the obtained apparent resistivity in far field is similar with that of magnetolluric (MT) method on a typical model. Therefore the forward modeling algorithm is suitable and accurate. On the other hand, NLCG is used in the 3D inversion and initial solution optimization is adopted, which can efficiently reduce computation time consumption. Numerical experiments illustrate that the distribution of anomalous bodies can be well recovered by the proposed inversion and the initial solution optimization can dramatically reduce computation time consumption.
朱成, 李桐林, 杨海斌, 刘永亮, 伍亮, 胡英才. 带地形频率域可控源电磁法三维反演研究[J]. 石油地球物理勘探, 2016, 51(5): 1031-1039.
Zhu Cheng, Li Tonglin, Yang Haibin, Liu Yongliang, Wu Liang, Hu Yingcai. 3D controlled source electromagnetic inversion with topography in the frequency domain. OGP, 2016, 51(5): 1031-1039.
Zhao Bangliu, He Zhanxiang, Wen Baihong. Non-seismic technique for direct oil-gas detection and its exploration practice. China Petroleum Exploration, 2005, 10(6):29-37.
Jia Jindou, He Zhanxiang, Kong Fanshu. The application and effects of comprehensive non-seismic geophysical and geochemical exploration techniques. OGP,1998, 33(5):625-631.
Dong Hao.Study of Three-dimensional Magnetotelluric Inversion Including Surface Topography Based on Finite-difference Method[D]. China University of Geosciences (Beijing), Beijing, 2013.
[8]
Wannamaker P E, Stodt J A, Rijo L. Two-dimen-sional topographic responses in magnetotellurics modeled using finite elements. Geophysics, 1986,51(11):2131-2144.
[9]
Mitsuhata Y, Uchida T. 3D magnetotelluric modeling using the T-Ω finite-element method. Geophysics, 2004, 69(1):108-119.
Lei Da.Studies and application of 2-D CSAMT modeling and inversion with a dipole source and topography. Chinese Journal of Geophysics, 2010, 53(4):982-993.
Yang Bo, Xu Yixian, He Zhanxiang et al. 3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method. Chinese Journal of Geophysics, 2012, 55(4):1390-1399.
Dong Hao, Wei Wenbo, Ye Gaofeng et al. Study of three-dimensional magnetotelluric inversion including surface topography based on finite-difference method. Chinese Journal of Geophysics, 2014, 57(3):939-952.
Gu Guanwen, Wu Wenli, Li Tonglin. Modeling for the effect of magnetotelluric 3D topography based on the vector finite-element method. Journal of Jilin University(Earth Science Edition), 2014, 44(5):1678-1686.
[17]
Newman G A, Alumbaugh D L. Frequency-domain modeling of airborne electromagnetic responses using staggered finite-differences. Geophysical Prospecting, 1995, 43(8):1021-1042.
[18]
Weiss C J, Constable S. Mapping thin resistors and hydrocarbons with marine EM methods, Part Ⅱ-Modeling and analysis in 3D. Geophysics,2006,71(6):G321-G332.
[19]
Barrett R, Berry M W, Chan T F et al. Templates for the solution of linear systems:building blocks for iterative methods. SIAM, 2001, 20-24.
[20]
Freund R W, Nachtigal N M. An implementation of the QMR method based on coupled two-term recurrences. SIAM Journal on Scientific Computing, 1994, 15(2):313-337.
Weng Aihua, Liu Yunhe, Jia Dingyu et al. Three-dimensional controlled source electromagnetic inversion using non-linear conjugate gradients. Chinese Journal of Geophysics, 2012, 55(10):3506-3515.
[22]
Newman G A, Alumbaugh D L. Electromagnetic modeling of subsurface 3D structures.//International Geoscience and Remote Sensing Symposium. IEEE, 1996, 1941-1944.
[23]
Smith J T. Conservative modeling of 3-D electromagnetic fields,Part 2:Biconjugate gradient solution and an accelerator. Geophysics, 1996, 61(5):1319-1324.
Zhu Cheng, Li Tonglin, Fan Cuisong. Research on divergence correction method in 3D numerical modeling of 3D controlled source electromagnetic fields. Global Geology, 2014, 33(4):916-924.
[25]
Newman G A, Alumbaugh D L. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients. Geophysical Journal International, 2002, 140(2):410-424.
[26]
Commer M, Newman G A, Williams K H et al. 3D induced-polarization data inversion for complex resistivity. Geophysics, 2011, 76(3):F157-F171.
[27]
Kelbert A, Egbert G D, Schultz A. Non-linear conjugate gradient inversion for global EM induction:resolution studies. Geophysical Journal International, 2008,173(2):365-381.
[28]
Newman G A, Boggs P T. Solution accelerators for large-scale three-dimensional electromagnetic inverse problems. Inverse Problems, 2004, 20(6):S151-S170.
[29]
Nocedal J, Wright S J. Numerical Optimization.Springer, Berlin, 2006, 56-60.