ISRN Geophysics

Volume 2013 (2013), Article ID 720839, 4 pages

http://dx.doi.org/10.1155/2013/720839

Research Article

Error Analysis in Measured Conductivity under Low Induction Number Approximation for Electromagnetic Methods

George Caminha-Maciel1,2 and Irineu Figueiredo1,3

¹Observatório Nacional (ON), MCT, Rua General José Cristino 77, 20921-400 Rio de Janeiro, RJ, Brazil

²Universidade Federal do Pampa (UNIPAMPA), Avenida Pedro Anunciação s/n°, 96570-000 Caçapava do Sul, RS, Brazil

³Universidade Estadual do Rio de Janeiro (UERJ), Rua São Francisco Xavier 524, 20550-900 Rio de Janeiro, RJ, Brazil

Received 28 August 2013; Accepted 5 November 2013

Academic Editors: G. Chisham, A. De Santis, G. Mele, and A. Tzanis

Copyright © 2013 George Caminha-Maciel and Irineu Figueiredo. 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.

We present an analysis of the error involved in the so-called low induction number approximation in the electromagnetic methods. In particular, we focus on the EM34 equipment settings and field configurations, widely used for geophysical prospecting of laterally electrical conductivity anomalies and shallow targets. We show the theoretical error for the conductivity in both vertical and horizontal dipole coil configurations within the low induction number regime and up to the maximum measuring limit of the equipment. A linear relationship may be adjusted until slightly beyond the point where the conductivity limit for low induction number () is reached. The equations for the linear fit of the relative error in the low induction number regime are also given.

1. Introduction

The induction method consists basically in determining subsurface rock conductivities with the help of electromagnetic fields generated by a coil at the Earth’s surface and by catching the response to this field from the conducting media under surface by using a reception coil [1–3].

From the Maxwell equations, in particular the Faraday induction law applied to an infinite homogenous half-plane, the subsurface rock conductivity can be estimated through the ratio between the magnetic field measured in the receiving coil and the magnetic field produced at the transmission coil with both at surface. Then we can take laterally distributed measurements along a transect for identifying conductivity-related anomalies. We can also get information on the vertical conductivity structure by varying the coil’s dipole configurations (vertical dipole or horizontal dipole) as well as by increasing the instrument height. This information is very useful in several geophysical problems as, for example, water prospecting or mapping pollution plumes.

The basic model for both configurations is described in Figure 1 where a transmission coil Tx with a given alternate electric current at a given frequency is located on the terrain (assumed to be an uniform semiplane) and a receiving coil Rx is located at a short distance from Tx. The time variation of the magnetic field , called primary magnetic field, produced by the electric current in the transmission coil generates a small alternate current in the soil. This electric current, on its turn, produces a magnetic field , called secondary field, which can be measured at the receiving coil together with the primary field.

Figure 1: Representation of the transmission (Tx) and reception (Rx) coils for both the vertical and the horizontal dipole configurations.

In general, the secondary magnetic field is a complicated function of the distance between coils , of the magnetic permeability of the medium —which will be considered exactly the same as for vacuum (A·m)—of the angular frequency of the oscillatory electric current , and of the soil conductivity . Under well-established conditions, technically known as “operation with low values of induction number,” the secondary magnetic field becomes a simple function of these variables when we consider only the quadrature component:

This result can be used when coils are in the vertical or horizontal dipole orientations. Actually, these two situations are described by different sets of equations for the secondary field; however, these equations give the same result for the component in quadrature of the ratio between the primary and secondary fields under low induction number conditions.

In order to understand the induction number we need to define the “electromagnetic skin depth” :

which is equal to the distance in a conducting medium an electromagnetic front wave has to travel in order to reduce the amplitude to of its value outside the conducting medium.

The induction number describes the distance between the transmission coil and the reception coil in units of skin depth :

The recent literature contains some debate on which values of (and consequently ) are valid for low induction number approximations [4, 5]. It also contains suggested corrections for nonzero height dependence and nonlinear departure for higher values of conductivity even above the low induction number regime [4].

In this work we introduce an error analysis for the induction method taking as a starting point the set of equations given by McNeill [1] for the field ratios on both coil dipole configurations. We also consider zero elevation of the coils and, as in Beamish [4], we discard medium magnetic and dielectric components as we are working in frequencies lower than 15 kHz. Those equations represent the ratio between primary and secondary fields as the response of a homogeneous earth. We calculate the relative difference between this conductivity value and that obtained through the low induction number approximation.

2. Methods

The equations for the primary and secondary field ratios for vertical and horizontal dipole configurations as given by McNeill [1] are

where

 , ,  = frequency (Hz),  = permeability of free space, .

Expressions (4) and (5) are multivalued functions of the conductivity, frequency, and intercoil spacing, and as functions of the complex exponential they have several branches in the complex plane. Therefore, as a matter of simplification, the induction number, a nondimensional parameter defined as the ratio between intercoil spacing and skin depth, is introduced, which will act as a conformal map simplifying equations.

To do this we begin by writing as a function of the induction number :

Remembering that

we have

By substituting in (4) the following expression for the vertical dipole configuration comes out:

As

then

and doing , it follows that

By taking only the quadrature of we have

Similarly, for the horizontal dipole configuration we obtain

The imaginary part (quadrature term) is related to the conductivity measurements under low induction number conditions. The magnitude of the secondary magnetic field is now directly proportional to the conductivity and its phase leads the primary field to 90°. Under low induction number () the expressions (4) and (5) can be expanded in a power series, and considering only the first terms the quadrature part (neglecting the in-phase term) can be written as [1]

This is the low induction number approximation key expression which leads to the conductivity value from the equipment readings.

The induction number is a function of the material conductivity for a given frequency and the separation between coils . Furthermore, is a function of the conductivity and of the product . For the EM34-3 and EM34-X (Geonics Ltd.) equipment settings all different possibilities of frequency and intercoil spacing give the same product (800 mS1/2 * m1/2). Since the low induction number condition () implies that we should have , then, considering the EM34-3 instrument that works in the frequencies 6400, 1600, and 400 Hz with distances between coils of 10, 20, and 40 m, respectively, it comes out that the value for the soil conductivity measured by this instrument should be quite lower than 400 mS/m.

We define the error associated with the low induction number approximation as the theoretical deviation from the exact model which is the electromagnetic response for a homogenous half-plane as follows:

From the induction number definition we have

and from (15)

In this way, the error associated with the approximation is

Substituting those values in expression (19) we obtain the relative theoretical error for the approximation. Figure 2 represents the plot of the error for the whole conductivity range (and a little further) mostly used in practical applications and shows the deviation from the linear response. We add two vertical black lines marking the limits of interest: the low induction number limit () and the usual limit of the equipment (1000 mS/m).

Figure 2: Relative percent error (%) as a function of conductivity (mS/m). It also shows error isolines for the error levels of 100%, 50%, and 25%. The vertical black lines mark two main limits: induction number equals to one (, for which  mS/m) and the EM34 instrument limit of 1000 mS/m.

From the studied lower limit through the whole low induction number regime (considering its limit as ) the relationship is clearly linear, but it also can easily be extended to 100 mS/m or even to 400 (mS)/m (), as it is shown in Figure 2.

Figure 3 shows the absolute apparent conductivity as a function of true conductivity until the linear response limit is surpassed. McNeill [1] shows this plot for a narrower range of conductivity values, and therefore the departure from linearity was not clear. The same linear behaviour is seen in the apparent conductivity (Figure 3) for vertical and horizontal dipole configurations, as seen by the fitted straight line to the logarithmic scaled error curves on the low induction number regime given by

where is given in percent (%) and in mS/m.

Figure 3: Apparent conductivity versus true conductivity in mS/m until slightly above the EM34 instrument limit. The vertical black lines mark two main limits: induction number equals to one (, for which  mS/m) and the EM34 instrument limit of 1000 mS/m. This curve shows a slight departure from the linear behaviour after the low induction number limit (~40 mS/m) until the limit for (~400 mS/m) and a strong departure from the linear behaviour after that.

3. Conclusions

It is well known that the vertical dipole configuration is more sensitive to conductivity anomalies in deep, while the horizontal one is more sensitive to near surface conductivity variations [1]. However, if we consider only the measured conductivity values, the error is lower for the horizontal than for the vertical dipole configuration as noted in Figure 2. For both configurations conductivity converges to the same values when conductivity tends to zero.

An approximately linear tendency can be observed over the whole low induction number range and up to 400 mS/m. Above this value there are strong deviations from linearity indicating that we are close to singularities of the complex function (not only the quadrature term). Despite a correction procedure suggested by Beamish [4], we believe that under statistical fluctuations of the measured signal the response could depart from the expected model leading to very high errors in the estimated conductivity due to the proximity of the function singularity.

It is worth noting that the determination of these bounds depends on the configurations (operation frequency and intercoil spacing) of the considered instrument.

Acknowledgments

This work was supported by the Brazilian agency CNPq through a postdoctoral scholarship to George Caminha-Maciel. Thanks are due to E. La Terra for the valuable discussions and M. Ernesto for helping in improving the paper. Any use of product names is for descriptive purposes only and does not imply any endorsement.

References

J. D. McNeill, “Electromagnetic terrain conductivity measurements at low induction numbers,” Technical Note TN-6, Geonics Ltd., Mississauga, Canada, 1980.

G. V. Keller and F. C. Firschknecht, Electrical Methods in Geophysical Prospecting, Pergamon, New York, NY, USA, 1966.

J. R. Waitt, “A note on the electromagnetic response of a stratified Earth,” Geophysics, vol. 27, pp. 382–385, 1962.

D. Beamish, “Low induction number, ground conductivity meters: A correction procedure in the absence of magnetic effects,” Journal of Applied Geophysics, vol. 75, no. 2, pp. 244–253, 2011.

J. B. Callegary, T. P. A. Ferré, and R. W. Groom, “Vertical spatial sensitivity and exploration depth of low-induction-number electromagnetic-induction instruments,” Vadose Zone Journal, vol. 6, no. 1, pp. 158–167, 2007.