International Journal of Geophysics
Volume 2018 (2018), Article ID 9891548, 7 pages
https://doi.org/10.1155/2018/9891548

Research Article

TEM Response of a Large Loop Source over the Multilayer Earth Models

A. K. Tiwari, S. P. Maurya, and N. P. Singh

Department of Geophysics, Faculty of Science, Banaras Hindu University, Varanasi 221005, India

Correspondence should be addressed to N. P. Singh; singhnpbhu@yahoo.co.in

Received 8 March 2018; Accepted 11 April 2018; Published 8 May 2018

Academic Editor: Rudolf A. Treumann

Copyright © 2018 A. K. Tiwari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The general expression of TEM response of large loop source over the layered earth models is not available in the literature for arbitrary source-receiver positions, except for the case of central loop and coincident loop configurations over the homogeneous earth model. In the present study, an attempt is made to present the TEM response of a large loop source over the layered earth model for arbitrary receiver positions. The frequency domain responses of large loop source over the layer earth model for arbitrary receiver positions are converted into the impulse (time derivative of magnetic field) TEM response using Fourier cosine or sine transform. These impulse TEM responses in turn are converted into voltage responses for arbitrary receiver positions, namely, central loop, arbitrary in-loop, and offset-loop TEM responses over the layered earth models. For checking the accuracy of the method, results are compared with the results obtained using analytical expression over a homogeneous earth model. The complete matching of both of the results suggests that the present computational technique is capable of computing TEM response of large loop source over the homogeneous earth model with high accuracy. Thereafter, the technique is applied for computation of TEM response of a large loop source over the layered earth (2-layer, 3-layer, and 4-layer) models for the central loop, in-loop, and offset-loop configurations and the results are presented in voltage decay form. The results depict their characteristic variations. These results would be useful for modeling and inversion of large loop TEM data over the layer earth models for all the possible configurations resulting from a large loop source.

1. Introduction

A large horizontal loop on or above the earth is one of the most widely used sources in TEM and airborne TEM (ATEM) methods. Using a large loop source, one can measure at the center of the loop, at any arbitrary point inside the loop, coincident loop, and at arbitrary offset-loop points. A number of studies on methods of computation of TEM response and modeling can be found in the literature [1–13]. All these studies are based on considering transmitter as a large circular loop and receiver either at the center of the loop and/or coincident with the loop, over the simple layer/homogenous earth models, and majority of them are based on the computation of TEM responses from their frequency responses only for the central and/or coincident loop configurations, because of the computational intricacies associated with the frequency domain response computation of large loop TEM source for arbitrary receiver positions (except for central and coincident loop configurations) over the layer earth models which involve product of two Basel functions that make the integral unstable and nonconvergent. Singh and Mogi [14, 15] for the first time presented a new computational method for computation of EM responses of large loop source over the layer earth models for arbitrary receiver positions. Thereafter, a number of researches appeared in the literature [16–24]. Thus, there was scarcity of EM modeling due to large loop sources over the layer earth models even in frequency domain, which resulted in a lack of studies in TEM response of large loop sources over the layer earth models for arbitrary receiver positions, except for the case of central loop and coincident loop configurations that were too over the homogeneous half-space because of nonavailability of analytical expression for arbitrary receiver positions, and sophisticated computational methods for computing the frequency domain responses from which TEM responses are usually derived. Moreover, during the process of converting the frequency domain EM response into time domain, even a small error of less than 1% in frequency response can produce a considerable difference in time domain EM response of large loop source over the layered earth models for all the configurations. Therefore, with the view of overcoming this drawback and filling the gap in the existing literature for TEM response of large loop sources, in this study, an attempt is made to present a reliable method for computation of TEM response of a large loop source over the layer earth models for any arbitrary receiver positions using the frequency response computation method described in [14, 15]. The computed TEM voltage response is compared with the available TEM responses over the homogeneous earth model for checking the accuracy and reliability of the method.

2. Theoretical Background

The plan view of large loop TEM method with central loop, in-loop, and offset-loop configurations over a layer earth is shown in Figure 1. The large loop presents a source loop and the small loop at the center of the source loop () represents receiver position corresponding to the central loop configuration. Similarly, the large loop presents a source loop and small loop at arbitrary position inside the source loop () represents receiver position corresponding to the in-loop configuration and small loop at arbitrary point outside the source loop () represents receiver position corresponding to arbitrary offset-loop configuration.

Figure 1: Schematic diagram showing all the possible configurations of large loop TEM system over a 3-layer earth model.

The frequency domain expressions of EM field components at a point on or above the surface of an -layered earth due to a finite horizontal circular loop of radius , carrying a current and placed at the height above the surface of layered earth model, can be found in [9]. The expression of field component at a measurement point on the surface of -layered earth (i.e., at ) can be written as

where with (intrinsic admittance of free space) and   (surface admittance at = 0).

For an -layer case, the surface admittances are given by the recurrence relation

and with , = , and .

Here, is source-receiver offset (measured from the center of the loop). For calculation purposes, is used in its exponential form for stability reasons [5].

Thereafter, the time derivative of vertical magnetic field is obtained by transforming the frequency domain solution of vertical magnetic field computed using [14] into the time domain solution using the Fourier cosine and/or sine transform as given in [25].

where and are the real and imaginary parts of the vertical magnetic field over a layered earth model in frequency domain. The components of vectors and are the resistivities and thickness of different layers of the layered earth model, and is the angular frequency.

In general, the data collected from a large loop TEM system consist of vertical voltage measurements made at various time intervals after the current in the transmitter is turned off. These voltage measurements are related to the time derivatives of vertical magnetic field ( in accordance with the following relation:

where is the area-turns product of the receiver coil. These voltage data can be further transformed into the apparent resistivity because sometimes it is preferable to use apparent resistivity transformation to have a direct relation with the geoelectrical section, as well as an initial estimate of the layer resistivities, which are often required in nonlinear inversion for interpretation of TEM data.

Therefore, starting with the computation of field (see (1)), using the method and algorithm described in [14, 15], we have computed the time derivative of vertical magnetic field using the Fourier cosine and sine transforms (see (3)) [8, 25–27]. Thereafter, the transient voltage response is computed using (4). Computations are performed for the homogeneous, 2-layer, 3-layer, and 4-layer resistive and conductive earth models and the results are presented in the following section.

3. Results and Discussions

3.1. Check for Validity of the Method

For checking the validity and accuracy of the method, we applied it for the computation of TEM (impulse) response of a large circular loop source over the surface of a homogeneous earth model for central loop configurations, and the result is compared with the published results for TEM response of a large loop source generated using the central loop analytical expression for impulse response [8, 9]. The analytic expression for the impulse response of the field at the center of a large loop source of radius (a) over the homogeneous earth model of conductivity σ can be written as [8, 9]

where , erf indicates the error function, and means the delay time after the current is turned off.

3.2. Illustration of TEM Response over Layer Earth Models

For illustrating the accuracy, applicability, and efficiency of the program for generating voltage response due to large loop TEM methods over the layer earth models for arbitrary receiver positions, we applied it for computing the TEM response over the homogeneous earth model of conductivity 0.01 S/m and compared the results with the published results [8, 9] (Figure 2). From Figure 2, it is evident that the computed results match very well the published results computed using the analytical expressions [9], thereby indicating the accuracy and validity of the method for computation of TEM response of large loop source over the homogeneous earth models.

Figure 2: Comparison of results computed using the present method for voltage response of vertical magnetic field at the center of a large circular loop source of radius, = 50 m, placed over the surface of a homogeneous earth model of conductivity ( = 0.01 S/m) with the results generated using the analytical expression given by ward and Hohmann (1988) and Singh et al. (2009).

Further, the method is applied for computation of TEM responses of a large loop source over the 2-layer, 3-layer, and 4-layer earth models for arbitrary receiver positions, that is, central loop, in-loop, and offset-loop configurations, for different loop sizes, and the results are presented in Figures 3–5.

(a)

(b)

(c)

Figure 3: Transient voltage response due to large circular loop source of radii 100 m and 200 m over the 2-layer earth model (as shown in the inset of the figure) for (a) central loop configuration, , (b) arbitrary in-loop configuration, , and (c) arbitrary offset-loop configuration, .

(a)

(b)

(c)

Figure 4: Transient voltage response of large circular loop source of radii 100 m and 200 m over the 3-layer earth model (as shown in the inset of the figure) for (a) central loop configuration, , (b) arbitrary in-loop configuration, , and (c) arbitrary offset-loop configuration, .

(a)

(b)

(c)

Figure 5: Transient voltage response of large circular loop source of radii 100 m and 200 m over the 4-layer earth model (as shown in the inset of the figure) for (a) central loop configuration, , (b) arbitrary in-loop configuration, , and (c) arbitrary offset-loop configuration, .

Figure 3 presents TEM responses of a large loop source over a two-layer earth model (as shown in the inset of the figure) for central loop, in-loop, and offset-loop configurations, respectively, for source loop radii 100 m and 200 m. The computed results depict characteristic features of TEM response over the two-layer earth models.

Figure 4 presents TEM response of a large loop source over the three-layer earth model (as shown in the inset of the figure) for loop sizes of radii 100 m and 200 m, for central loop, in-loop, and offset-loop configurations. The computed results depict characteristic response over the three-layer earth models.

Figure 5 depicts TEM response of a large loop source over the four-layer earth model (as shown in the inset of the figure) for source loop sizes of radii 100 m and 200 m, for central loop, in-loop, and arbitrary offset-loop configurations, respectively. The results depict characteristic response of large loop sources over the four-layer earth models.

From these figures (Figures 3–5), it is clear that, with the increase in loop size, the TEM response shows smooth and well-defined characteristics. Moreover, it is also noticed that the central loop TEM responses are more regular as compared to the offset-loop and arbitrary in-loop responses.

4. Conclusion

The present article describes a simple and sophisticated method for computation of TEM response of a large loop source over layered earth model for arbitrary receiver positions, that is, at the center of the loop, at arbitrary in-loop, and at arbitrary offset-loop points. The method is based on conversion of frequency domain results computed using EMLCLTR program of Singh and Mogi [15] into the time derivative of the vertical magnetic field, that is, impulsive TEM responses. The method is simple and suitable for computation of TEM response of a large loop source at any arbitrary receiver locations, that is, central loop, in-loop, and offset-loop configurations. During its frequency domain computation, it has the option of computing frequency domain responses with or without the inclusion of displacement current factor and hence is a more reliable and accurate one.

For illustrating the nature and characteristics of large loop TEM responses over the layer earth models at arbitrary receiver positions, results are presented for the TEM response at source-receiver offsets ,  , and pertaining to the central loop, in-loop, and offset-loop configurations over 2-layer, 3-layer, and 4-layer earth models. The results depict their characteristic variations. From the results, it is noticed that the voltage response curves are more regular for central loop, followed by offset loop, and the least for the in-loop configurations. Further, it is also observed that, with the increase in source loop size, the voltage responses become smoother.

This study would enable the prospect, development, and use of loop-in-loop method (in-loop method) and loop-offset-loop method (offset-loop method) along with the well-developed central loop and coincident loop methods using large loop sources, and thus it would enhance the applicability and cost-effectiveness of the large loop source transient electromagnetic method for exploration and applied geophysics applications.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

One of the authors (A. K. Tiwari) is thankful to the Department of Geophysics, Banaras Hindu University, Varanasi, for providing the opportunity and facility for completion of this work.

References

A. V. Christiansen and N. B. Christensen, “A quantitative appraisal of airborne and ground-based transient electromagnetic (TEM) measurements in Denmark,” Geophysics, vol. 68, no. 2, pp. 523–534, 2003.

H. Huang and G. J. Palacky, “Damped least-squares inversion of time-domain airborne EM data based on singular value decomposition,” Geophysical Prospecting, vol. 39, no. 6, pp. 827–844, 1991.

T. Ingeman-Nielsen and F. Baumgartner, “CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys,” Computers & Geosciences, vol. 32, no. 9, pp. 1411–1419, 2006.

H. Jang, H. Jang, K. H. Lee, and H. J. Kim, “Step-off, vertical electromagnetic responses of a deep resistivity layer buried in marine sediments,” Journal of Geophysics and Engineering, vol. 10, no. 2, Article ID 025011, 2013.

J. H. Knight and A. P. Raiche, “Transient electromagnetic calculations using the Gaver- Stehfest inverse Laplace transform method.,” Geophysics, vol. 47, no. 1, pp. 47–50, 1982.

H. F. Morrison, R. J. Phillips, and D. P. Obrien, “Quantitative interpretation of transient electromagnetic fields over a layered half space,” Geophysical Prospecting, vol. 17, no. 1, pp. 82–101, 1969.

A. P. Raiche, “Transient electromagnetic field computations for polygonal loops on layered earths,” Geophysics, vol. 52, no. 6, pp. 785–793, 1987.

N. P. Singh, M. Utsugi, and T. Kagiyama, “TEM response of a large loop source over a homogeneous earth model: A generalized expression for arbitrary source-receiver offsets,” Pure and Applied Geophysics, vol. 166, no. 12, pp. 2037–2058, 2009.

S. H. Ward and G. W. Hohmann, “Electromagnetic theory for geophysical applications,” in Electromagnetic Methods in Applied Geophysics, M. N. Nabighian, Ed., vol. 1, pp. 131–311, Society of Exploration Geophysicists, Tulsa, Oklahoma, 1988.

G.-Q. Xue, C.-Y. Bai, and X. Li, “Extracting the virtual reflected wavelet from TEM data based on regularizing method,” Pure and Applied Geophysics, vol. 169, no. 7, pp. 1269–1282, 2012.

G. Q. Xue, S. Yan, and N. N. Zhou, “Theoretical study on the errors caused by dipole hypothesis of large-loop TEM response,” Diqiu Wuli Xuebao, vol. 54, no. 9, pp. 2389–2396, 2011.

M. S. Zhdanov, “Electromagnetic geophysics: Notes from the past and the road ahead,” Geophysics, vol. 75, no. 5, pp. 75A49–75A66, 2010.

M. S. Zhdanov, D. A. Pavlov, and R. G. Ellis, “Localized S-inversion of time-domain electromagnetic data,” Geophysics, vol. 67, no. 4, pp. 1115–1125, 2002.

N. P. Singh and T. Mogi, “EMLCLLER-a program for computing the EM response of a large loop source over a layered earth model,” Computers & Geosciences, vol. 29, no. 10, pp. 1301–1307, 2003.

N. P. Singh and T. Mogi, “Electromagnetic response of a large circular loop source on a layered Earth: a new computation method,” Pure and Applied Geophysics, vol. 162, no. 1, pp. 181–200, 2005.

S. Cristina and M. Parise, “Fast calculation of theoretical response curves for induction depth sounding,” in Proceedings of the European Microwave Week 2009, EuMW 2009: Science, Progress and Quality at Radiofrequencies - 39th European Microwave Conference, (EuMC '09), pp. 1567–1570, October 2009.

M. Jamie, “Note on: ‘EMLCLLER—A program for computing the EM response of a large loop source over a layered earth model’ by N.P. Singh and T. Mogi, Computers & Geosciences 29 (2003) 1301–1307,” Computers & Geosciences, vol. 96, pp. 236–246, 2016.

M. Parise, “Fast computation of the forward solution in controlled-source electromagnetic sounding problems,” Progress in Electromagnetics Research, vol. 111, pp. 119–139, 2011.

M. Parise, “Full-wave analytical explicit expressions for the surface fields of an electrically large horizontal circular loop antenna placed on a layered ground,” IET Microwaves, Antennas & Propagation, vol. 11, no. 6, pp. 929–934, 2017.

Y. Qi, L. Huang, X. Wu, G. Fang, and G. Yu, “Effect of loop geometry on TEM response over layered earth,” Pure and Applied Geophysics, vol. 171, no. 9, pp. 2407–2415, 2014.

J. T. Ratnanather, J. H. Kim, S. Zhang, A. M. Davis, and S. K. Lucas, “Algorithm 935: IIPBF, a MATLAB toolbox for infinite integral of products of two Bessel functions,” ACM Transactions on Mathematical Software (TOMS), vol. 40, no. 2, p. 14, 2014.

N. P. Singh and T. Mogi, “EMDPLER: A F77 program for modeling the EM response of dipolar sources over the non-magnetic layer earth models,” Computers & Geosciences, vol. 36, no. 4, pp. 430–440, 2010.

J. Deun and R. Cools, “Note on electromagnetic response of a large circular loop source on a layered earth: A new computation method,” in Pure and Applied Geophysics, N. P. Singh and T. Mogi, Eds., vol. 164, pp. 1107–1111, 2007.

G.-Q. Xue, X. Li, L. J. Gelius, Z.-P. Qi, N.-N. Zhou, and W.-Y. Chen, “A new apparent resistivity formula for in-loop fast sounding TEM theory and application,” Journal of Environmental & Engineering Geophysics, vol. 20, no. 2, pp. 107–118, 2015.

G. A. Newman, G. W. Hohmann, and W. L. Anderson, “Transient electromagnetic response of a three-dimensional body in a layered earth,” Geophysics, vol. 51, no. 8, pp. 1608–1627, 1986.

W. L. Anderson, “Numerical integration of related Hankel transforms of orders 0 and 1 by adaptive digital filtering,” Geophysics, vol. 44, no. 7, pp. 1287–1305, 1979.

W. L. Anderson, “Nonlinear least-squares inversion of transient soundings for a central induction loop system (Program NLSTCI),” U.S. Geological Survey 82, 1982.