Some issues on 2.5-D transient electromagnetic forward

Abstract

Numerical calculation for two integral transforms in 2.5-D transient electromagnetic forward is a difficult and key task, namely, the inverse Fourier transform and the inverse Laplace transform. Some effective algorithms for them were described. Based on the known algorithms in DC resistivity on wave-number distribution and selection, we proposed a principle on how to choose the least wave-number concerning the central-loop transient electromagnetic method. First, observe the behavior of transformation function curve with regard to wave-number in Fourier domain. In the light of its asymptote, ascertain the coverage scope of wave-number. Compared with analytic solution, the least wave-number in Fourier domain can be derived. Furthermore, the Laplace numerical inversion algorithm which needs only a few Laplace variables in pure real domain was also introduced here. The procedure was applied to forward modeling on transient electromagnetic field of a vertical magnetic dipole over uniform half-space to demonstrate them effectiveness and general applicability.

Keywords

transient electromagnetic method / 2.5-D forward modeling / inverse fourier transform / inverse laplace transform / wave-number

References

[1] LuoYan-zhong, ZhangGui-qingApplication of Electronic Computer in Electrical Prospecting [M], 1987, Wuhan, Wuhan College of Geology Press(in Chinese)

[2] WangHua-junMethods Study on 2.5-D Forward and Inversion for Transient ElectromAgnetic with Central Loop [D], 2001, Wuhan, China University of Geosciences(in Chinese)

[3] MengYong-liang, LuoYan-zhongLuoYan-zhong. Two dimensional direct finite element algorithm of electrical dipolar source CSAMT method [A]. Contributions to the Exploration Geophysics and Geochemistry [C], 1996, Beijing, Geological Press: 103-114(in Chinese)

[4] WangHua-jun, LuoYan-zhong. Algorithm of a 2.5 dimensional finite element method for transient electromagnetic with a central loop [J]. Chinese Journal of Geophysics, 2003, 46(6): 855-862(in Chinese)

[5] LuoYan-zhong, MengYong-liang. Some problems on resistivity modeling for two-dimensional structures by the finite element method [J]. Chinese Journal of Geophysics, 1986, 29(6): 613-621(in Chinese)

[6] XuShi-zheFinite Element in Geophysics [M], 1994, Beijing, Science Press(in Chinese)

[7] DuanBen-chun, ZhuYong-sheng. The application of optimized wave number to direct calculation of point source two-dimensional finite element [J]. Geophysical & Geochemical Exploration, 1994, 18(3): 186-191(in Chinese)

[8] LeppinM. Electromagnetic modeling of 3-D sources over 2-D inhomogeneities in the time domain [J]. Geophysics, 1992, 57(8): 994-1003

[9] UnsworthM J, TravisB J, ChaveA D. Electromagnetic induction by a finite electric dipole source over a 2-D earth [J]. Geophysics, 1993, 58(2): 198-214

[10] EverettM E, EdwardsR N. Transient marine electromagnetics: The 2.5-D forward problem [J]. Geophys J Int, 1992, 113: 545-561

[11] MitsuhataY, UchidaT, AmanoH. 2.5-D inversion of frequency-domain electromagnetic data generated by a grounded-wire source [J]. Geophysics, 2002, 67(6): 1753-1768

[12] YanShu, ChenMing-sheng, FuJun-mei. Direct time-domain numerical analysis of transient electromagnetic fields [J]. Chinese Journal of Geophysics, 2002, 45(2): 275-284 in Chinese)

[13] LuoYan-zhong, ZhangShen-gye, WangWeiping. A research on one-dimension forward for aerial electromagnetic method in time domain [J]. Chinese Journal of Geophysics, 2003, 46(5): 719-724(in Chinese)

[14] StehfestH. Numerical inversion of Laplace transforms [J]. Communications of the ACM, 1970, 13(1): 47-49

[15] CrumpK S. Numerical inversion of Laplace transforms using a Fourier series approximation [J]. J ACM, 1976, 23(1): 89-96

[16] GuptasarmaD. Computation of time-domain response of a polarizable ground [J]. Geophysics, 1982, 47(11): 1574-1576

[17] AndersonW L. Fast Hankel transforms using related and lagged convolutions [J]. ACM Trans on Math Software, 1982, 8(4): 344-368

[18] KnightJ H, RaicheA P. Transient electromagnetic calculation using Gaver-Stehfest inverse Laplace transform method [J]. Geophysics, 1982, 47(1): 47-50

[19] HeGuang-yu, WangWei-hong. An application of an accurate and efficient numerical inversion of the Laplace transform in modern well-testing analysis [J]. Acta Petrolei Sinica, 1995, 16(1): 96-104(in Chinese)

[20] HeGuang-yu, LiLun. An application of AWG method of numerical inversion of Laplace transform for flow in a fluid finite-conductivity vertical fractures [J]. Petroleum Exploration and Development, 1995, 22(6): 47-50(in Chinese)

[21] RuanBai-yao. An application of Guptasarma algorithm in transient electromagnetic forward [J]. Journal of Guilin Institute of Technology, 1996, 16(2): 167-170(in Chinese)

[22] LuoYan-zhong, ChangYan-jun. A rapid algorithm for G-S transform [J]. Chinese Journal of Geophysics, 2000, 43(5): 684-690 in Chinese)

[23] KaufmanA A, KellerG VFrequency and Transient Soundings [M], 1987, Beijing, Geological Press(in Chinese)