Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation
Liu Yun1, Wang Xu-Ben2, and Wang Yun1
1. State Key Laboratory of Ore Deposit Geochemistry, Institute of Geochemistry Chinese Academy of Sciences, Gui Yang 550002, China.
2. School of Geophysics, Chengdu University of Technology, Chengdu 610059, China.
Abstract To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time-domain transient electromagnetic secondary field of the line source based on the DuFort-Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force (Emf), and determined the forward time step. By extending upward the earth–air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the flat layered model with high precision—the maximum relative error is less than 0.01% between our method and the analytical method—and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.
The research is supported by the National High Technology Research and Development Program (863 Program) (2009AA06Z108).
Cite this article:
LIU Yun,WANG Xu-Ben,WANG Bin. Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation[J]. APPLIED GEOPHYSICS, 2013, 10(2): 134-144.
[1]
Adhidjaja, J. I., Hohmann, G. W., and Oristaglio, M. L., 1985, Two-dimensional transient electromagnetic responses: Geophysics, 50, 2849 - 2861.
[2]
Birtwistle, G. M., 1968, The explicit solution of the equation of heat conduction: Comput. J., 11, 317.
[3]
DuFort, E. C., and Frankel, S. P., 1953, Stabilitv conditions in the numerical treatment of parabolic differential equations: Math. Tables and Other Aids to Comput, 7, 135 - 152.
[4]
Guptasarma, D., 1982, Computation of the time-domain response of a polarizable ground: Geophysics, 47(11), 1574 - 1576.
[5]
He, G. Y., and Gao, Y. L., 2002, Visual Fortran common collection numerical algorithms: Science Press, Beijing.
[6]
Liu, Y., and Wang, X. B., 2010, Complex terrain effect and studying on 2.5 D TEM: Proceedings of the 4th International Conference on Environmental and Engineering Geophysics, Chengdu, China. 452 - 456.
[7]
Nabighian, M. N., (Translator Zhao J. X.), 1992, Electromagnetic Methods in Applied Geophysics, Volume 1: Theory: Geological Publishing House Beijing.
[8]
Niu, Z. L., 1992, The Theory of Time-Domain Electromagnetic Methods: Central South University of Technology Press, Changsha.
[9]
Oristaglio, M. L., 1982, Diffusion of electromagnetic fields into the earth from a line source of current: Geophysics, 47, 1585 - 1592.
[10]
Oristaglio, M. L., and Hohmann, G. W., 1984, Diffusion of electromagnetic fields in two-dimensional earth: A finite difference approach. Geophysics, 49, 870 - 894.
[11]
Piao, H. R., 1990, Theory of electromagnetic: Geology Press, Beijing.
[12]
Ruan, B. Y., and Xu, S. Z., 1996, Finite-difference calculations of transient time domain electromagnetic field of a two-dimensional earth for a line source of current excitation: Geological Journal of China Universities (in Chinese), 2(4), 437 - 447.
[13]
Ruan, B. Y., 1995, Asymptotic Inversion Schemes for Electrical Sounding Data and Numerical Method for 2D Time Domain Electromagnetic Response: Ph.D. Thesis, Qingdao, Ocean University of Qingdao, Qingdao.
[14]
Ruan, B. Y., 1996, The application of Guptasarma method on TEM forward problem: Journal of Guilin Institute of Technology (in Chinese), 16(2), 167 - 170.
[15]
Stoyer, C. H., and Greenfield, R. J., 1976, Numerical solutions of the response of a two-dimensional earth to an oscillating magnetic dipole source: Geophysics, 41(3), 519 - 530.
[16]
Wang, H. J., and Luo Y. Z., 2003, Algorithm of a 2.5 dimensional finite element method for transient electromagnetic with a central loop: Chinese J. Geophysics (in Chinese), 46(6), 855 - 862.
[17]
Wang, T., and Hohmann, G. W., A finite-difference time-domain solution for three-dimensional electromagnetic modeling: Geophysics, 1993, 58(6), 797 ? 809.
[18]
Xiong, B., 2004, The research on large-loop transient electromagnetic 2.5D forward algorithm and some key data technologies: Ph.D. Thesis (in Chinese), China University of Geosciences (Wuhan).
[19]
Xiong, B., and Luo, Y. Z., 2006, Finite element modeling of 2·5DTEMwith block homogeneous conductivity: Chinese J. Geophys (in Chinese), 49(2), 590 - 597.
[20]
Yan, S., Chen, M. S., and Fu, J. M., 2002, Direct time-domain numerical analysis of transient electromagnetic field: Chinese J. Geophysics (in Chinese), 45(2), 275 - 284.
2013, Vol. 10 
