Next Article in Journal
A Mathematical Model for Enhancing CO2 Capture in Construction Sector Using Hydrated Lime
Next Article in Special Issue
Study of the Genesis Process and Deep Prospecting Breakthrough in the Gouli Ore Concentration of the East Kunlun Metallogenic Belt Using Audio Magnetotelluric Data
Previous Article in Journal
Correction: Leeman et al. U-Pb LA-ICP-MS Zircon Dating of Crustal Xenoliths: Evidence of the Archean Lithosphere Beneath the Snake River Plain. Minerals 2024, 14, 578
Previous Article in Special Issue
Location Prediction Study of Fluorite Ore in Shallow Cover Area: Evidence from Integrated Geophysical Surveys
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Moment Estimation from Time Domain Electromagnetic Data

by
Marc A. Vallée
* and
Mouhamed Moussaoui
Geo Data Solutions GDS Inc., 4402 Louis-B. Mayer, Laval, QC H7P 0G1, Canada
*
Author to whom correspondence should be addressed.
Minerals 2024, 14(9), 888; https://doi.org/10.3390/min14090888
Submission received: 18 July 2024 / Revised: 24 August 2024 / Accepted: 29 August 2024 / Published: 30 August 2024

Abstract

:
Moment representations have been proposed to facilitate the interpretation of geophysical time domain electromagnetic responses. We present a new methodology for estimating these moments from field data for different system waveforms when on-time and off-time measurements are available. Quadrature impulse response moments are estimated by a recursive relation involving moments of the input waveform and moments of the observed response. After adapting this method to time domain electromagnetic applications—in particular, MEGATEM and AeroTEM (AirTEM) airborne electromagnetic systems—we present the results from applying this method on synthetic and real data collected over the Reid–Mahaffy test site in northern Ontario, Canada.

1. Introduction

Time domain electromagnetic (TDEM) methods are used in exploration geophysics to search for subsurface mineralized zones that are conductive [1]. Their principle is as follows: an electric current pulse, associated with a primary magnetic field through Ampère’s law, circulates in the transmitter loop, as illustrated in Figure 1. According to Faraday’s law of induction, this primary magnetic field generates secondary currents in conductive materials in the subsurface. These secondary currents are associated with a secondary magnetic field that is detected by the receiver loop. As the pulse is of limited duration, the secondary signal can be measured while the pulse is on (on-time) and after the pulse is turned off (off-time). Typically, only off-time measurements are recorded since on-time measurements are affected by the primary field. In airborne electromagnetics (AEMs), the transmitter (source) and receiver are either located on the aircraft or are towed below the aircraft (Figure 2).
There have been many advancements since the early development of TDEM techniques to localize mineralized zones with EM measurements. Better understanding of the electromagnetic (EM) principles involved [1,2,3] and the development of computer techniques have led to improvements in various domains. Ref. [4] proposed to localize anomalies and interpret them by comparing them with the results of analog model simulations. Following advances in computer technology, Ref. [5] developed an inversion program to invert AEM data or to retrieve the parameters of the model by numerical curve-fitting using the numerical model of a plate in free space developed by [6]. The model based on the inversion of a plate was extended by [7] to consider host rock to be conductive, layered earth. Refs. [8,9] proposed an algorithm for inverting a sphere in layered earth. The thin sheets in the layered earth model of [7] were included by [10] for an initial model for the probabilistic inversion of AEM data for basement conductors.
Some of these elements were incorporated into the three-step process proposed by [11] and refined by [12], which is as follows: (1) compress the data to decay time-constant space, (2) transform the response to a conductivity–depth image, and (3) parameterize local anomalies with plate-like conductors. Similarly, Ref. [13] proposes to invert AEM data in two phases in order to determine background conductivity models with layered-earth inversions and interpret discrete conductors as electric or magnetic dipoles, while the approach in [14] approximates the subsurface current system using a 3D subsurface grid of either 3D magnetic or electric dipoles, for which locations and orientations are solved by inversion. Some of the most elaborate interpretation techniques invert the complete airborne dataset with sophisticated mathematical techniques and solve the forward model with an integral equation [15,16,17], finite volume method [18], or randomly undersampled forward model [19].
Despite these advances, elaborate EM interpretation techniques are not generalized to all AEM datasets as they are complex and costly in terms of time and computer resources and because of the enormous dataset sizes. For these reasons, simple interpretation tools have the potential of assessing data quality quickly, locating discrete conductors, and facilitating the use of more complex interpretation methods. Some of these tools include the use of parameters such as the apparent conductivity of the subsurface [20,21], the decay time constant [22,23], and inductive and resistive limits [24]. At the inductive limit, current flows do not penetrate conductive bodies, and at the resistive limit, the response varies linearly with conductivity. Resistive limit estimation was used by [21,24,25] as the first step of an EM interpretation process. Refs. [26,27] extended the concepts of these limits to the moments of the quadrature impulse response, which are time integrals of the response signal, as the zero-order and first-order moments correspond to the inductive and resistive limit, respectively. Subsequent authors developed solutions for the moments of a sphere in free-space [27], in homogeneous half-space and excited by a source located on the surface of the ground [28,29], thick horizontal layer [30], and in one-dimensional earth with finite conductance [31]. In the year 2000, CGG and Geophysical Software Solutions collaborated to develop an application called EM-Q based on the research done by [32] on the uniform field approximation of the response of a sphere, for which [8] compared the results to the solution attained in [5]. The moments were also used by [25] as an approximation for the development of a 3D EM inversion program.
However, as TDEM measurements are done with systems using different primary current waveforms, the estimation of the moments from the quadrature impulse response is sometimes problematic. Ref. [27] developed a solution by approximating the TDEM waveform as a simple current box and applied this technique to the data integral from a MEGATEM AEM system [33], for which the current waveform is a half-sine. However, this approach is not applicable to more elaborate waveforms such as those of the AeroTEM AEM system [34], for which the current waveform is a symmetrical triangular pulse. Furthermore, Ref. [35] encountered difficulties with the definition of zero time while analyzing MEGATEM data, which resulted in a noisy moment estimation. Ref. [28] only considers step-response measurements, from which it is possible to estimate moments directly. In conclusion, although there are extensive theoretical developments on the use of moments in the field of TDEM interpretation, estimation of the moments from various waveforms on synthetic or field data has been neglected, slowing the application of this method on real data.
For these reasons, we propose a method for directly estimating the quadrature impulse moments from any waveform and measured signals. This method should increase the use of the moment method and better outline its limitations. The proposed algorithm evaluates the moments of the quadrature impulse response with a recursive relation involving the moments of the pulse current derivative and the moments of the quadrature on-time and off-time responses. Consequently, it requires that on-time measurements are available. After a noise analysis, we apply this method to synthetic AirTEM results (the same waveform as AeroTEM) with a sphere model in layered earth, illustrating the conditions for which moments can be useful for outlining a discrete conductor. This is followed by application of the moment method on field data from MEGATEM and AirTEM systems collected over the Reid–Mahaffy test site, Ontario, Canada.

2. Theory

2.1. Moment Estimation

The quadrature impulse response of a TDEM system can be developed as follows [27]. The EM step response of the subsurface can be written as
s ( t ) = B u ( t ) h ( t )
where t is time, B is the response at t = 0 , u ( t ) is the unit step-on function, and h ( t ) is a dimensionless function characterizing the secondary decay associated with currents in the ground. The total impulse response r ( t ) is the time derivative of s ( t ) :
r ( t ) = B δ ( t ) h ( t ) + u ( t ) h ( t ) t
Analogous to frequency-domain systems, Ref. [26] names the first and second terms of this equation the in-phase and quadrature components, respectively, and focuses on the second term. This is represented as
i ( t ) = B u ( t ) h ( t ) t
The moments of the quadrature impulse response i ( t ) are defined as
I n = 0 t n i ( t ) d t
In practice, the primary waveform is different from a step response and is convolved with the impulse response. Furthermore, the receiver coil detects the derivative of the magnetic field. This can be written as
y ( t ) = x ( t ) i ( t )
where x ( t ) represents the waveform derivative, * is the convolution operator, and y ( t ) is the measured secondary field, from which the on-time primary field and the secondary in-phase field have been removed. This equation represents the result y ( t ) of a causal convolution of the system waveform x ( t ) and the system transfer function i ( t ) , where x ( t ) = i ( t ) = 0 for t < 0
y ( t ) = 0 i ( h ) x ( t h ) d h = 0 i ( t h ) x ( h ) d h
Following [36,37], we develop the moments of the convolution response as
Y n = 0 y ( t ) t n d t = 0 0 i ( h ) x ( t h ) t n d h d t
If we replace u = t h , then t = u + h and d u = d t ; then, the expression becomes
Y n = h h 0 i ( h ) x ( u ) ( u + h ) n d h d u
Since x ( u ) = 0 for u < 0 and h , we can re-write the equation as
Y n = 0 0 i ( h ) x ( u ) ( u + h ) n d h d u
Developing the power of a binomial expansion, it is possible to separate the two integrals in the following way:
Y n = 0 0 i ( h ) x ( u ) k = 0 n n k u n k h k d h d u Y n = k = 0 n n k 0 x ( u ) u n k d u 0 i ( h ) h k d h
where
n k = n ! k ! ( n k ) !
If we define the input and impulse response moments as
X n = 0 x ( t ) t n d t
and
I n = 0 i ( t ) t n d t ,
we can write
Y n = k = 0 n n k X n k I k .
For X 0 0 , if we develop
Y n = k = 0 n 1 n k X n k I k + X 0 I n ,
we can evaluate I n from X k and Y k , k [ 0 , n ] , as
I n = Y n k = 0 n 1 n k X n k I k X 0 .
However, in the case of some AEM systems, the current waveform or derivative is symmetric, with positive and negative lobes, such as in MEGATEM and AeroTEM, and X 0 = 0 . Then we need to use the following relation:
Y n = k = 0 n 2 n k X n k I k + n n 1 X 1 I n 1
from which we deduce
I n 1 = Y n k = 0 n 2 n k X n k I k n n 1 X 1
or
I n = Y n + 1 k = 0 n 1 n + 1 k X n + 1 k I k n + 1 n X 1
The next subsections demonstrate how this relation can be applied to the MEGATEM and AeroTEM systems.

2.2. MEGATEM

MEGATEM [33] is a fixed-wing system with a transmitter loop wound around the nose, wingtips, and tail of a DeHavilland Dash-7. The transmitter loop covers an area of 406 m2, while the receiver comprises three magnetic coils that are towed 50 m below and 131 m behind the aircraft. As stated before, the current pulse is a half-sine. The primary field extraction [38], also called stripping, is done using a reference waveform that is measured at altitude and by minimizing the functional Φ , where
Φ = n T ( t i ) α R ( t i ) = n Q ( t i )
T ( t i ) is the total field response, R ( t i ) is the reference waveform, and α is an unknown factor that is estimated from the minimization. As [38] points out, this procedure estimates the on-time in-phase and quadrature responses. Accordingly, this approach provides the information required for moment estimation using the proposed method.

2.3. AeroTEM

The AeroTEM system [34] is a helicopter-towed EM system with coincident transmitter and receivers characterized by a triangular current pulse (Figure 3) for which the derivative is a constant with alternating polarity. The voltage sensed by the receiver, which is surrounded by a bucking coil to cancel the primary field, is a series of step-responses, i.e., a first step of positive polarity, a second step of double-negative polarity, and a last step of positive polarity. The on-time field is measured as a sum of two voltages separated by half the period.
The response derivative can be seen as the sum of three step-response of alternating polarities and is expressed as
z ( t ) = B u ( t ) h ( t ) + 2 B u ( t P / 2 ) h ( t P / 2 ) B u ( t P ) h ( t P )
where P is the pulse width, and B, u ( t ) , and h ( t ) are defined as in Equation (1). This can be decomposed according to time:
z ( t ) = B h ( t ) 0 < t < P / 2 = B h ( t ) 2 h ( t P / 2 ) P / 2 < t < P = B h ( t ) 2 h ( t P / 2 ) + h ( t P ) P < t
The in-phase x ( t ) and quadrature y ( t ) responses are, respectively,
x ( t ) = B h ( 0 ) 0 < t < P / 2 = B h ( 0 ) P / 2 < t < P = 0 P < t
and
y ( t ) = B h ( t ) h ( 0 ) 0 < t < P / 2 = B h ( t ) 2 h ( t P / 2 ) + h ( 0 ) P / 2 < t < P = B h ( t ) 2 h ( t P / 2 ) + h ( t P ) P < t
The on-time sum measurement typical of the AeroTEM system is
d ( t ) = z ( t ) + z ( t + P / 2 ) = B h ( t ) B h ( t + P / 2 ) 2 h ( t ) = B h ( t + P / 2 ) h ( t )
If we assume that h ( t + P / 2 ) < < h ( t ) , we can approximate that d ( t ) B h ( t ) . Using this approximation, y ( t ) can be estimated from Equation (23) as
y ( t ) d ( 0 ) d ( t ) y ( t + P / 2 ) 2 d ( t ) d ( 0 )

2.4. Incomplete Moments

In practice, response moments are estimated between the early and late times, as represented by the following equation [28]:
Y n = 0 t 1 y ( t ) t n d t + t 1 t n y ( t ) t n d t + t n y ( t ) t n d t
where [28] labels the first and last terms the “head” and the “tail”, respectively, and the middle term is the incomplete moment. As demonstrated by [28], if the waveform is a step, the tail and head can be estimated from the early and late-time apparent resistivity. However, for a complex waveform, this method cannot be applied, and the solution we develop in this paper can only estimate the incomplete moments, which are more complex than for a step response.
The error in this estimation can be assessed by comparing the synthetic moment estimates with the analytic moment estimates provided by [30]. For this purpose, we present results for the AirTEM system [39] (shown in Figure 2), which is characterized by a triangular waveform similar to that of the AeroTEM system. This system employs a horizontal loop transmitter with a diameter of 4.25 m and a vertical magnetic coil receiver (dB/dt) located on the side of the transmitter loop. These are towed by a helicopter (Figure 2), and the EM system loop operates approximately 40 m above the ground. The transmitter waveform is a 1.85 ms bipolar triangular pulse with a frequency of 90 Hz. The time derivative of the magnetic field is sampled with 40 off-time channels between 0.02 and 3.3 ms and on-time differences covering the first half of each on-time pulse. The time domain synthetic half-space response is estimated using the program Airbeo [7] with Hankel filters [40] and frequency-to-time domain transformations [41]. The first two moments are estimated from the half-space response using the method proposed in this paper and are compared to the analytic solutions from [30], as seen in Figure 4.

2.5. Normalization

The moments of a wire loop were developed in [27], which illustrates the distribution of the moments for a simple electromagnetic model. For a wire loop with a time constant τ w = L / R , where L and R are the inductance and resistance, respectively, of the loop, the moments are
M n = B n ! τ w n
where B is the inductive limit and n is the moment order. This implies that time normalization will have an impact on the variation with the moment order. Accordingly, smaller time constants can be enhanced by scaling the moments. The moments of a scaled integral are developed as follows. If we define the scaled version of Y n as Y ν n , where
Y ν n = 0 y ( t / ν ) t n d t ,
if t / ν = τ , this becomes
Y ν n = 0 y ( τ ) ( ν τ ) n ν d τ = ν n + 1 0 y ( τ ) τ n d τ = ν n + 1 Y n
Scaling corresponds to multiplication of the moments by the power of a time factor. The impulse response moment from Equation (16) becomes
I ν n = Y ν n + 1 k = 0 n 1 n + 1 k X ν n + 1 k I ν k n + 1 n X ν 1
Using a recurrent relation, it can be shown that
I ν n = ν n I n
This allows enhancing the moments of higher orders, if required.

2.6. Moment Estimation from Noisy Data

In practice, moments are estimated from noisy data, and an analysis of the proposed method requires an evaluation of the impact of noisy data on moment estimation. First, we must estimate the impacts of noisy data on the estimations of waveforms and output moments. These are estimated by numerical integration of the waveform and output signals, which can be written as
X n = 0 l X ( t l ) t l n Δ t l Y n = 0 m Y ( t m ) t m n Δ t m
where X ( t ) and Y ( t ) are the waveform and output signals, respectively, and t l , t l are time samples. Assuming independent measurements and neglecting correlations, the variances of these signals are estimated using the variance formula [42]:
s f = f x 2 s x 2 + f y 2 s y 2 + f z 2 s z 2 + . . .
where s f is the variance of the function f = f ( x , y , z , . . . ) , and s x , s y , s z are the variances of the parameters x , y , z . Developing the partial derivatives of Equation (33), we write
X n ( t l ) t l = n t l n 1 X ( t l ) Δ t l Y n ( t m ) t m = n t m n 1 Y ( t m ) Δ t m
Then, Equation (34) can be developed as
s X n = 0 l n t l n 1 X ( t l ) Δ t l 2 s l 2 s Y n = 0 m n t m n 1 Y ( t m ) Δ t m 2 s m 2
where s X n and s Y n are the moment waveform and output variances respectively. Equation (36) shows that these moment variances increase with the moment order.
The impulse moment estimate variance can be developed in the following way. If we simplify Equation (19) for n = 1 , we have
I 0 = Y 1 X 1
From Equation (34), we write
s I 0 = I 0 Y 1 2 s Y 1 2 + I 0 X 1 2 s X 1 2 s I 0 = s Y 1 2 ( X 1 ) 2 + ( Y 1 ) 2 s X 1 2 ( X 1 ) 4
In general, from Equation (19), the variance becomes
s I n = I n Y n + 1 2 s Y n 2 + k = 0 n 1 I n I k 2 s I k 2 + k = 1 n I n X k 2 s X k 2
Without developing this expression in detail, it is easy to again conclude that the moment variances increase with the moment order.

3. Results

3.1. Synthetic Models

We apply this method on synthetic results from numerical models of a localized target in a conductive environment. We use the model of [9] to represent the sphere response to an AirTEM system in layered earth. The sphere has a radius of 50 m and a conductivity of 10 S/m, and its center is located 100 m below the surface of the ground, while the AirTEM system is at an altitude of 40 m. We start with a sphere in a half-space with 0.1 mS/m conductivity, for which the response is shown in Figure 5. The sphere has a distinctive moment response. The peak moments of a similar sphere with varying conductivity are displayed in Figure 6. This figure shows that the sphere’s peak moments are negligible for conductivity below 1 S/m, even if the sphere has a meaningful off-time response, and the moment reaches a maximum value at about 8 S/m and decreases for conductivities above this value.
We now examine the sphere moments below a conductive overburden of 25 m thickness with varying conductivity. Figure 7 shows the peak nomogram for the overburden with varying conductivity (“anomalous” represents the sphere-only response). At low overburden conductivity, the response comes mainly from the sphere. As the overburden conductivity increases, the sphere peak becomes negative, while the overburden response increases. Figure 8 show the corresponding moments. At low overburden conductivity, the sphere moments dominate. However, above 0.1 S/m, the sphere moments decrease as the overburden moments increase strongly. Figure 8 suggests that we can still observe the sphere response if the overburden is 0.1 S/m (100 mS/m). This is confirmed by the subsequent model in Figure 9, wherein the top layer has a conductivity of 100 mS/m and 25 m thickness and is above a half-space with 0.1 mS/m conductivity. Even if the sphere’s off-time responses are barely visible, the moment responses are clearly visible.

3.2. Field Data

Moment estimation was tested on field data collected with the MEGATEM and the AirTEM systems over the Reid–Mahaffy test site in northern Ontario, Canada [43]. This site is characterized by the presence of several conductors located in a resistive basement below a variably conductive overburden.

3.2.1. MEGATEM

MEGATEM data collected over the Reid–Mahaffy test site in 2002 are available from [44]. In particular, the moments of Line 15, shown in Figure 10a, were analyzed by [27] based on B field estimation. For comparison, moments estimated using our proposed method are presented in Figure 10b. The various bedrock conductors can be clearly identified. Another profile of Line 3 in Figure 11 shows three distinct conductors whose moment signatures are remarkably significant.

3.2.2. AirTEM

AirTEM data were collected for testing purposes over the Reid–Mahaffy test site in 2019. The profiles for Line 15 are displayed in Figure 12a. AirTEM is a helicopter TDEM system with a much smaller loop flown closer to the ground. Consequently, the bedrock conductors are less apparent and are partly masked by the overburden response. Despite this fact, some deep conductors to the north of the profile are visible and can be observed in the moment response shown in Figure 12b. Line 3 conductor responses are shown in Figure 13. Clear moment responses are observed over these distinct conductors.

4. Discussion

Although not exhaustive, the results presented in this study allow a discussion of the strengths and limitations of the moment method. The methodology applied in this paper has been used successfully to interpret synthetic and field data. An important requirement of the method is the availability of on-time quadrature measurements. For the MEGATEM system, there is an established procedure for this purpose. For AeroTEM and AirTEM, some approximation is required. For other systems, this could be a major limitation to the implementation of this method.
The synthetic models show that for small background conductivity and a good conductor (in our models, more than a few S/m), the conductor moments are easily detected. However, as shown in Figure 6, a bad conductor (in our models, less than 1 S/m) will not be detected. Furthermore, when the conductivity of the background is significant, even a good conductor response is not detectable. Another important consideration when applying this technique is to consider the noise present in the data. The noise analysis shows that noise increases with the moment order, which may explain why high-order MEGATEM moments are less noisy than equivalent AirTEM moments, as the MEGATEM data, originating from a much more powerful system, are less noisy than the AirTEM data.
The application of this method on datasets from different systems over the same area allowed for a pertinent comparison of the performances of the systems in question. In the case presented, the Reid–Mahaffy area is covered with significant conductive overburden, and the AirTEM system, with smaller coincident loops, has some difficulty penetrating the overburden. Consequently, the discrete conductor moments are scattered. On the other hand, the MEGATEM system—operating at the same frequency but with much higher power, a larger transmitter loop, and greater transmitter–receiver separation—detects distinct conductors that are deeper below the overburden easily, comparatively speaking. These conductors are clearly defined by moments.
As mentioned in Section 2.4, the method presented in this paper provides an estimate of incomplete moments. This must be considered if results are compared with analytic models, and in that case [28], total moment relations must be applied. However, for numerical models wherein system time windows are included, observed moments can be directly related to model moments.

5. Conclusions

The method presented in this paper extends the application of moment estimation to arbitrary EM system waveforms and generalizes the use of moments; it can provide data quality assessments, conductor locations, and preliminary interpretations. A limitation of this approach is the availability of on-time measurements. For systems that can resolve this limitation, such as MEGATEM and AeroTEM (AirTEM), we show that this method successfully estimates moments and allows a comparison of system performance over a given area. We advocate that this method could be extended to other EM systems—either airborne, ground, or borehole—when on-time measurements are available. Our model also shows that for a conductor to be detected, the conductivity range must be appropriate: good conductors are easier to detect as long as there is contrast with the surrounding environment.
This algorithm can be easily implemented for data post-processing and can be used to glean information about the EM data quality quickly. As such, it is seen a supplement to simplified interpretation methods that only provide single parameters like decay time constants or apparent conductivities.
As this method relies on on-time and off-time measurements, improvements in data collection, like reducing noise and improving calibration methods, should improve the moment estimates. This is particularly true for on-time measurements: these are more difficult to collect, which explains why not all systems measure the full waveform.

Author Contributions

Conceptualization, M.A.V. and M.M.; methodology, M.A.V.; software, M.A.V.; validation, M.A.V. and M.M.; investigation, M.A.V.; resources, M.A.V.; writing—original draft preparation, M.A.V.; writing—review and editing, M.A.V. and M.M.; visualization, M.A.V.; supervision, M.M.; project administration, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

MEGATEM data from Reid–Mahaffy from the given reference are available from the Ontario Geology website. AirTEM data are available on request from the first author.

Acknowledgments

The authors acknowledge Khorram Khan for revising the English.

Conflicts of Interest

The authors are both working for Geo Data Solutions GDS Inc., which commercializes the AirTEM system.

References

  1. Grant, F.S.; West, G.F. Interpretation Theory in Applied Geophysics; McGraw-Hill Book Company: New York, NY, USA, 1965. [Google Scholar]
  2. Ward, S.H.; Hohmann, G.W. Electromagnetic theory for geophysical applications. In Electromagnetic Methods in Applied Geophysics: Volume I, Theory; Nabighian, M.N., Ed.; Society of Exploration Geophysicists: Tulsa, OK, USA, 1987; pp. 131–311. [Google Scholar]
  3. Kaufman, A.A.; Alekseev, D.; Oristaglio, M. Principles of Electromagnetic Methods in Surface Geophysics; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar]
  4. Palacky, G.J. Interpretation of INPUT AEM Measurements in Areas of Conductive Overburden. Geophysics 1975, 40, 490–502. [Google Scholar] [CrossRef]
  5. Keating, P.B.; Crossley, D.J. The inversion of time-domain airborne electromagnetic data using the plate model. Geophysics 1990, 55, 705–711. [Google Scholar] [CrossRef]
  6. Annan, A.P. The Equivalent Source Method for Electromagnetic Scattering Analysis and Its Geophysical Application. Ph.D. Thesis, Memorial University of Newfoundland, St. John’s, NL, Canada, 1974. [Google Scholar]
  7. Raiche, A.; Sugeng, F.; Wilson, G. Practical 3D AEM inversion—The P223F software suite. ASEG Ext. Abstr. 2007, 1, 1–5. [Google Scholar]
  8. Vallée, M.A. New developments in AEM discrete conductor modelling and inversion. Explor. Geophys. 2015, 46, 97–111. [Google Scholar] [CrossRef]
  9. Vallée, M.A.; Moussaoui, M. Modelling the electromagnetic response of a sphere located in a layered earth. Explor. Geophys. 2023, 54, 362–375. [Google Scholar] [CrossRef]
  10. Hauser, J.; Gunning, J.; Annetts, D. Probabilistic inversion of airborne electromagnetic data for basement conductors. Geophysics 2016, 81, E389–E400. [Google Scholar] [CrossRef]
  11. Macnae, J.; King, A.; Stolz, N.; Osmakoff, A.; Blaha, A. Fast AEM data processing and inversion. Explor. Geophys. 1998, 29, 163–169. [Google Scholar] [CrossRef]
  12. Macnae, J. Automated airborne EM anomaly picking and 3D model fitting. ASEG Ext. Abstr. 2015, 1, 1–4. [Google Scholar] [CrossRef]
  13. Sattel, D.; Reid, J. Modelling of Airborne Em Anomalies with Magnetic and Electric Dipoles Buried in a Layered Earth. Explor. Geophys. 2006, 37, 254–260. [Google Scholar] [CrossRef]
  14. Kolaj, M.; Smith, R.S. Inductive electromagnetic data interpretation using a 3D distribution of 3D magnetic or electric dipoles. Geophysics 2017, 82, E187–E195. [Google Scholar] [CrossRef]
  15. Cox, L.H.; Wilson, G.A.; Zhdanov, M.S. 3D inversion of airborne electromagnetic data. Geophysics 2012, 77, WB59–WB69. [Google Scholar] [CrossRef]
  16. Cox, L.H.; Zhdanov, M.S.; Pitcher, D.H.; Niemi, J. Three-Dimensional Inversion of Induced Polarization Effects in Airborne Time Domain Electromagnetic Data Using the GEMTIP Model. Minerals 2023, 13, 779. [Google Scholar] [CrossRef]
  17. Cox, L.H.; Zhdanov, M.S.; Prikhodko, A. Inversion for 3D Conductivity and Chargeability Models Using EM Data Acquired by the New Airborne TargetEM System in Ontario, Canada. Minerals 2024, 14, 237. [Google Scholar] [CrossRef]
  18. Oldenburg, D.W.; Haber, E.; Shekhtman, R. Three dimensional inversion of multisource time domain electromagnetic data. Geophysics 2013, 78, E47–E57. [Google Scholar] [CrossRef]
  19. Ren, X.; Lai, M.; Wang, L.; Yin, C.; Liu, Y.; Su, Y.; Zhang, B.; Ben, F.; Huang, W. A fast 3-D inversion for airborne EM data using pre-conditioned stochastic gradient descent. Geophys. J. Int. 2023, 234, 737–754. [Google Scholar] [CrossRef]
  20. Fraser, D.C. Dighem resistivity techniques in airborne electro-magnetic mapping. In Airborne Resistivity Mapping; Palacky, G.J., Ed.; Geol. Surv. Canada Paper 86-22; Geological Survey of Canada: Ottawa, Canada, 1986; pp. 63–69. [Google Scholar]
  21. Annan, A.P.; Smith, R.S.; Lemieux, J.; O’Connell, M.D.; Pedersen, R.N. Resistive-limit, time-domain AEM apparent conductivity. Geophysics 1996, 61, 93–99. [Google Scholar] [CrossRef]
  22. Nelson, P.H.; Morris, D.B. Theoretical response of a time-domain, airborne, electromagnetic system. Geophysics 1969, 34, 729–738. [Google Scholar] [CrossRef]
  23. Martinez, J.M.; Smith, R.; Diaz Vazquez, D. On the time decay constant of AEM systems: A semi-heuristic algorithm to validate calculations. Explor. Geophys. 2020, 51, 94–107. [Google Scholar] [CrossRef]
  24. Macnae, J.; King, A.; Stolz, N.; Klinkert, P. 3-D EM Inversion to the Limit. In Three-Dimensional Electromagnetics; Oristaglio, M., Spies, B., Eds.; Society of Exploration Geophysicists: Tulsa, OK, USA, 1999; pp. 489–501. [Google Scholar]
  25. Schaa, R.; Fullagar, P.K. Rapid, approximate 3D inversion of transient electromagnetic (TEM) data. SEG Tech. Program Expand. Abstr. 2010, 650–654. [Google Scholar]
  26. Smith, R.S.; Lee, T.J. The impulse-response moments of a conductive sphere in a uniform field, a versatile and efficient electromagnetic model. Explor. Geophys. 2001, 32, 113–118. [Google Scholar] [CrossRef]
  27. Smith, R.S.; Lee, T.J. The moments of the impulse response: A new paradigm for the interpretation of transient electromagnetic data. Geophysics 2002, 67, 1095–1103. [Google Scholar] [CrossRef]
  28. Schaa, R.; Fullagar, P.K. Vertical and horizontal resistive limit formulas for a rectangular-loop source on a conductive half-space. Geophysics 2012, 77, E91–E99. [Google Scholar] [CrossRef]
  29. Fullagar, P.K.; Schaa, R. Analytic formulas for complete and incomplete first-order TEM moments below ground in a conductive half-space. Geophysics 2020, 86, E1–E11. [Google Scholar] [CrossRef]
  30. Smith, R.S.; Lee, T.J. Using the moments of a thick layer to map conductance and conductivity from airborne electromagnetic data. J. Appl. Geoph. 2002, 49, 173–183. [Google Scholar] [CrossRef]
  31. Lee, T.J.; Smith, R.S. Multiple-order moments of the transient electromagnetic response of a one-dimensional earth with finite conductance—Theory. Explor. Geophys. 2021, 52, 1–15. [Google Scholar] [CrossRef]
  32. Smith, R.; Hyde, C.; Lee, T.; Almond, R. Impulsive moments at work. ASEG Ext. Abstr. 2003, 1, 1–7. [Google Scholar] [CrossRef]
  33. Smith, R.S.; Fountain, D.; Allard, M. The MEGATEM fixed-wing transient EM system applied to mineral exploration: A discovery case history. First Break 2003, 21, 71–75. [Google Scholar] [CrossRef]
  34. Balch, S.J.; Boyko, W.P.; Paterson, N.R. The AeroTEM airborne electromagnetic system. Lead. Edge 2012, 22, 562–566. [Google Scholar] [CrossRef]
  35. Smith, R.S.; Lee, T.J. Multiple-order moments of the transient electromagnetic response of a one-dimensional earth with finite conductance—The Gaussian variation applied to a field example. Explor. Geophys. 2022, 53, 262–274. [Google Scholar] [CrossRef]
  36. Nash, J.E. Systematic Determination of Unit Hydrograph Parameters. J. Geoph. Res. 1959, 64, 111–115. [Google Scholar] [CrossRef]
  37. Koutsoyiannis, D.; Xanthopoulos, T. On the Parametric Approach to Unit Hydrograph Identification. Water Resour. Manag. 1989, 3, 107–128. [Google Scholar] [CrossRef]
  38. Smith, R.S. On removing the primary field from fixed-wing time-domain airborne electromagnetic data: Some consequences for quantitative modelling, estimating bird position and detecting perfect conductors. Geophys. Prospect. 2001, 49, 405–416. [Google Scholar] [CrossRef]
  39. Balch, S.J.; Samson, C.; Sanchez, J.L. Correcting for HTEM altitude variations by transmitter mutual inductance normalization. SEG Tech. Program Expand. Abstr. 2016, 2139–2143. [Google Scholar]
  40. Christensen, N.B. Optimized fast Hankel transform filters. Geophys. Prospect. 1990, 38, 545–568. [Google Scholar] [CrossRef]
  41. Raiche, A. Modelling the time-domain response of AEM systems. Explor. Geophys. 1998, 29, 103–106. [Google Scholar] [CrossRef]
  42. Ku, H.H. Notes on the use of propagation of error formulas. J. Res. Natl. Bur. Stand. 1966, 70C, 262. [Google Scholar] [CrossRef]
  43. Witherly, K.; Irvine, R.; Godbout, M. Reid Mahaffy Test Site, Ontario Canada: An example of benchmarking in airborne geophysics. SEG Tech. Program Expand. Abstr. 2004, 1202–1204. [Google Scholar]
  44. Ontario Geological Survey. Ontario Airborne Geophysical Surveys, Magnetic and Electromagnetic Data, Timmins Area MEGATEM®; Ontario Geological Survey, Geophysical Data Set 1041; Ontario Geological Survey: Sudbury, ON, Canada, 2002. [Google Scholar]
Figure 1. Schematic of a transmitter–receiver target system.
Figure 1. Schematic of a transmitter–receiver target system.
Minerals 14 00888 g001
Figure 2. Example of a helicopter TDEM system (AirTEM) towed below the aircraft. The receiver is located on the same bird as the transmitter.
Figure 2. Example of a helicopter TDEM system (AirTEM) towed below the aircraft. The receiver is located on the same bird as the transmitter.
Minerals 14 00888 g002
Figure 3. Waveform and derivative waveform of the AeroTEM system.
Figure 3. Waveform and derivative waveform of the AeroTEM system.
Minerals 14 00888 g003
Figure 4. Comparison of synthetic and analytic half-space low-order moments.
Figure 4. Comparison of synthetic and analytic half-space low-order moments.
Minerals 14 00888 g004
Figure 5. (a) On-time and off-time responses and (b) moments of a sphere with 10 S/m conductivity and a 50 m radius located in a conductive half-space of 0.1 mS/m 100 m below the surface.
Figure 5. (a) On-time and off-time responses and (b) moments of a sphere with 10 S/m conductivity and a 50 m radius located in a conductive half-space of 0.1 mS/m 100 m below the surface.
Minerals 14 00888 g005
Figure 6. (a) On-time and off-time peak response nomogram and (b) associated moments for a sphere model with varying conductivity and 50 m radius located in a conductive half-space of 0.1 mS/m 100 m below the surface.
Figure 6. (a) On-time and off-time peak response nomogram and (b) associated moments for a sphere model with varying conductivity and 50 m radius located in a conductive half-space of 0.1 mS/m 100 m below the surface.
Minerals 14 00888 g006
Figure 7. (a) Total and (b) anomalous on-time and off-time peak responses of a sphere model of 10 S/m located in two-layered earth with varying conductivity and 25 m thickness above a half-space with 0.1 mS/m conductivity.
Figure 7. (a) Total and (b) anomalous on-time and off-time peak responses of a sphere model of 10 S/m located in two-layered earth with varying conductivity and 25 m thickness above a half-space with 0.1 mS/m conductivity.
Minerals 14 00888 g007
Figure 8. (a) Total and (b) anomalous peak response moments of a sphere model of 10 S/m located in two-layered earth with varying conductivity and 25 m thickness above a half-space with 0.1 mS/m conductivity.
Figure 8. (a) Total and (b) anomalous peak response moments of a sphere model of 10 S/m located in two-layered earth with varying conductivity and 25 m thickness above a half-space with 0.1 mS/m conductivity.
Minerals 14 00888 g008
Figure 9. (a) On-time and off-time responses and (b) moments of a sphere model of 10 S/m, 50 m radius, and 100 m depth located in two-layered earth with 100 mS/m conductivity and 25 m thickness above a half-space with 0.1 mS/m conductivity.
Figure 9. (a) On-time and off-time responses and (b) moments of a sphere model of 10 S/m, 50 m radius, and 100 m depth located in two-layered earth with 100 mS/m conductivity and 25 m thickness above a half-space with 0.1 mS/m conductivity.
Minerals 14 00888 g009
Figure 10. MEGATEM X component (a) on-time and off-time responses and (b) moments estimated from L15 of the Reid–Mahaffy test site.
Figure 10. MEGATEM X component (a) on-time and off-time responses and (b) moments estimated from L15 of the Reid–Mahaffy test site.
Minerals 14 00888 g010
Figure 11. MEGATEM X component (a) on-time and off-time responses and (b) moments estimated from L3 of the Reid–Mahaffy test site.
Figure 11. MEGATEM X component (a) on-time and off-time responses and (b) moments estimated from L3 of the Reid–Mahaffy test site.
Minerals 14 00888 g011
Figure 12. AirTEM (a) on-time and off-time responses and (b) estimated moments from L15 of the Reid–Mahaffy test site.
Figure 12. AirTEM (a) on-time and off-time responses and (b) estimated moments from L15 of the Reid–Mahaffy test site.
Minerals 14 00888 g012
Figure 13. AirTEM (a) on-time and off-time responses and (b) estimated moments from L3 of the Reid–Mahaffy test site.
Figure 13. AirTEM (a) on-time and off-time responses and (b) estimated moments from L3 of the Reid–Mahaffy test site.
Minerals 14 00888 g013
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vallée, M.A.; Moussaoui, M. Moment Estimation from Time Domain Electromagnetic Data. Minerals 2024, 14, 888. https://doi.org/10.3390/min14090888

AMA Style

Vallée MA, Moussaoui M. Moment Estimation from Time Domain Electromagnetic Data. Minerals. 2024; 14(9):888. https://doi.org/10.3390/min14090888

Chicago/Turabian Style

Vallée, Marc A., and Mouhamed Moussaoui. 2024. "Moment Estimation from Time Domain Electromagnetic Data" Minerals 14, no. 9: 888. https://doi.org/10.3390/min14090888

APA Style

Vallée, M. A., & Moussaoui, M. (2024). Moment Estimation from Time Domain Electromagnetic Data. Minerals, 14(9), 888. https://doi.org/10.3390/min14090888

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop