Abstract — Integral expressions of electromagnetic (EM) field components due to a large circular loop source carrying an alternating (ac) current and placed on or above the surface of a layered earth model are transformed to such suitable forms that facilitate numerical computation of field response in quasi-static as well as non-quasi-static regions. The improper integrals occurring in expressions of EM field components are evaluated by converting these integrals into the convergent integrals using the process of subtraction or addition of an integral expression inside the integral sign and subsequently adjusting it or its equivalent analytic expression outside the integral sign. The adjusted integral expressions, in turn, are evaluated using the functional relationships described in this paper. The computation method based on this formulation takes into consideration the effects of both conduction as well as displacement currents, and is well suitable for any position of the source loop either in the air or on the surface of the model, in contrary to the earlier methods which face convergence problem. Moreover, the formulation is equally efficient for computing the EM response at any arbitrary receiver position either inside or outside the source loop. For illustrating the accuracy and applicability of the method and studying the nature of EM response of a loop source over a layered earth model, we have applied it for the computation of amplitude and phase of Hz field over the various 2-layer and 3-layer models. Results show their characteristic variations, and depict good resolution for the subsurface layering. The results are in agreement with those of the published results for the quasi-static region, and are new extension of quasi-static variation in the non-quasi-static region. The agreement of computed results with published results demonstrates the accuracy of the method. Moreover, this is the initial presentation of numerical results for an arbitrary in-loop point (other than the center point) inside a large circular loop source.
[1]
M. Poddar,et al.
A rectangular loop source of current on a two-layered earth
,
1982
.
[2]
G. W. Hohmann,et al.
4. Electromagnetic Theory for Geophysical Applications
,
1987
.
[3]
H. Frank Morrison,et al.
ELECTROMAGNETIC FIELDS ABOUT A LOOP SOURCE OF CURRENT
,
1970
.
[4]
D. F. Hays,et al.
Table of Integrals, Series, and Products
,
1966
.
[5]
G. Keller,et al.
Frequency and transient soundings
,
1983
.
[6]
Irene A. Stegun,et al.
Handbook of Mathematical Functions.
,
1966
.
[7]
I. S. Gradshteyn,et al.
Table of Integrals, Series, and Products
,
1976
.
[8]
William H. Press,et al.
Numerical Recipes: FORTRAN
,
1988
.
[9]
H. J. Hagger,et al.
Electromagnetic Waves in Stratified Media
,
1996
.
[10]
Walter L. Anderson,et al.
Computer program; numerical integration of related Hankel transforms of orders O and 1 by adaptive digital filtering
,
1979
.
[11]
D. Guptasarma,et al.
New digital linear filters for Hankel J0 and J1 transforms
,
1997
.
[12]
L. Milne‐Thomson.
A Treatise on the Theory of Bessel Functions
,
1945,
Nature.
[13]
Ronald Cools,et al.
A Matlab Implementation of an Algorithm for Computing Integrals of Products of Bessel Functions
,
2006,
ICMS.
[14]
H. F. Morrison,et al.
QUANTITATIVE INTERPRETATION OF TRANSIENT ELECTROMAGNETIC FIELDS OVER A LAYERED HALF SPACE
,
1969
.
[15]
John Knight,et al.
Transient electromagnetic calculations using the Gaver-Stehfest inverse Laplace transform method
,
1982
.
[16]
A. A. Kaufman,et al.
The theory of inductive prospecting
,
2001
.