Modeling induced polarization effects in helicopter time-domain electromagnetic data: Field case studies

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. B49-B61 ◽  
Author(s):  
Vladislav Kaminski ◽  
Andrea Viezzoli

Induced polarization (IP) effects are becoming more evident in time-domain helicopter airborne electromagnetic (AEM) data thanks to advances in instrumentation, mainly due to improvements in the signal-to-noise ratio and hence better data quality. Although the IP effects are often manifested as negative receiver voltage values, which are easy to detect, in some cases, IP effects can distort recovered transients in other ways so they may be less obvious and require careful data analysis and processing. These effects represent a challenge for modeling and inversion of the AEM data. For proper modeling of electromagnetic transients, the chargeability of the subsurface and other parameters describing the dispersion also need to be taken into consideration. We use the Cole-Cole model to characterize the dispersion and for modeling of the IP effects in field AEM data, collected by different airborne systems over different geologies and exploration targets, including examples from diamond, gold, and base metal exploration. We determined how multiparametric inversion techniques can simultaneously recover all four Cole-Cole parameters, including resistivity [Formula: see text], chargeability [Formula: see text], relaxation time [Formula: see text], and frequency parameter [Formula: see text]. The results obtained are in good agreement with the ancillary information available. Interpretation of the IP effects in AEM data is therefore seen by the authors as providing corrected electrical resistivity distributions, as well as additional information that could assist in mineral exploration.

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. E31-E50 ◽  
Author(s):  
Andrea Viezzoli ◽  
Vladislav Kaminski ◽  
Gianluca Fiandaca

We have developed a synthetic multiparametric modeling and inversion exercise undertaken to study the robustness of inverting airborne time-domain electromagnetic (TDEM) data to extract Cole-Cole parameters. The following issues were addressed: nonuniqueness, ill posedness, dependency on manual processing and the effect of constraints, and a priori information. We have used a 1D layered earth model approximation and lateral constraints. Synthetic simulations were performed for several models and the corresponding Cole-Cole parameters. The possibility to recover these models by means of laterally constrained multiparametric inversion was evaluated, including recovery of chargeability distributions from shallow and deep targets based on analysis of induced polarization (IP) effects, simulated in airborne TDEM data. Different scenarios were studied, including chargeable targets associated with the conductive and resistive environments. In particular, four generic models were considered for the exercise: a sulfide model, a kimberlite model, and two generic models focusing on the depth of investigation. Our study indicated that, in cases when relaxation time ([Formula: see text]) values are in the range to which the airborne electromagnetic is most sensitive (e.g., approximately 1 ms), it is possible to recover deep chargeable targets (to depths more than 130 m) in association with high electrical conductivity and in resistive environments. Furthermore, it was found that the recovery of a deep conductor, masked by a shallower chargeable target, became possible only when full Cole-Cole modeling was used in the inversion. Lateral constraints improved the recoverability of model parameters. Finally, modeling IP effects increased the accuracy of recovered electrical resistivity models.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. E75-E88 ◽  
Author(s):  
Changhong Lin ◽  
Gianluca Fiandaca ◽  
Esben Auken ◽  
Marco Antonio Couto ◽  
Anders Vest Christiansen

Recently, the interest in the induced polarization (IP) phenomenon in airborne time-domain electromagnetic (ATEM) data has increased considerably. IP may affect the ATEM data significantly and mask underlying geologic structures. To simulate 2D airborne IP data, a 2D finite-element forward-modeling algorithm has been developed with the dispersive conductivity described by the well-known Cole-Cole model. We verify our algorithm by comparison with the 1D solution of the AarhusInv code. Two-dimensional forward responses on six synthetic models, mimicking archetypal 2D conductive and chargeable anomalies, have been generated, and the results indicate that 2D IP affects the data significantly. Differences between the 2D IP responses and the 1D IP responses are evident above the 2D anomalies and at their edges. These differences are similar to what is found when comparing 2D and 1D forward responses over conductive 2D anomalies without considering IP. We evaluate an effective robust inversion scheme to recover the 2D IP parameters using the 1D laterally constrained inversion (LCI) scheme. The inversion of the synthetic data using the robust scheme indicates that not only can the IP parameters be recovered, but also the IP inversions can provide more accurate resistivity sections than a resistivity-only inversion, in terms of resistivity values and anomaly thickness/depth. The field example from Hope Bay area in Canada is even more valuable, considering that part of the profile consists of only negative data, which cannot be inverted with a resistivity-only scheme. Furthermore, the edge effects at the anomaly boundaries are less pronounced in the IP parameters than in the resistivity parameter on the synthetic models with more conductive backgrounds.


Geophysics ◽  
1996 ◽  
Vol 61 (1) ◽  
pp. 66-73 ◽  
Author(s):  
Richard S. Smith ◽  
Jan Klein

Airborne induced‐polarization (IP) measurements can be obtained with standard time‐domain airborne electromagnetic (EM) equipment, but only in the limited circumstances when the ground is sufficiently resistive that the normal EM response is small and when the polarizability of the ground is sufficiently large that the IP response can dominate the EM response. Further, the dispersion in conductivity must be within the bandwidth of the EM system. One example of what is hypothesized to be IP effects are the negative transients observed on a GEOTEM® survey in the high arctic of Canada. The dispersion in conductivity required to explain the data is very large, but is not inconsistent with some laboratory measurements. Whether the dispersion is caused by an electrolytic or dielectric polarization is not clear from the limited ground follow‐up, but in either case the polarization can be considered to be induced by eddy currents associated with the EM response of the ground. If IP effects are the cause of the negative transients in the GEOTEM data, then the data can be used to estimate the polarizabilities in the area.


Geophysics ◽  
1981 ◽  
Vol 46 (6) ◽  
pp. 932-933 ◽  
Author(s):  
T. Lee

Recently Pelton et al. (1978) used a Cole‐Cole relaxation model to simulate the transient voltages that are observed during an induced‐polarization survey. These authors took the impedance of the equivalent circuit Z(ω) to be [Formula: see text]They then gave the expression for the transient voltage [Formula: see text] as [Formula: see text]In equation (2), [Formula: see text] was misprinted as [Formula: see text]. In these equations, [Formula: see text] and [Formula: see text], [Formula: see text] and τ are constants to be determined for the given model. [Formula: see text] is the height of the step current that will flow in the transmitter. A disadvantage of equation (2) is that it is only slowly convergent for large t/τ. Pelton et al. (1978) used a τ which ranged from [Formula: see text] to [Formula: see text]. The purpose of this note is to provide an alternative expression for [Formula: see text] that is valid only at the later stages but which does not have this disadvantage. The trivial case of c = 1.0 is ignored.


Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. E213-E225 ◽  
Author(s):  
Gianluca Fiandaca ◽  
Esben Auken ◽  
Anders Vest Christiansen ◽  
Aurélie Gazoty

Time-domain-induced polarization has significantly broadened its field of reference during the last decade, from mineral exploration to environmental geophysics, e.g., for clay and peat identification and landfill characterization. Though, insufficient modeling tools have hitherto limited the use of time-domain-induced polarization for wider purposes. For these reasons, a new forward code and inversion algorithm have been developed using the full-time decay of the induced polarization response, together with an accurate description of the transmitter waveform and of the receiver transfer function, to reconstruct the distribution of the Cole-Cole parameters of the earth. The accurate modeling of the transmitter waveform had a strong influence on the forward response, and we showed that the difference between a solution using a step response and a solution using the accurate modeling often is above 100%. Furthermore, the presence of low-pass filters in time-domain-induced polarization instruments affects the early times of the acquired decays (typically up to 100 ms) and has to be modeled in the forward response to avoid significant loss of resolution. The developed forward code has been implemented in a 1D laterally constrained inversion algorithm that extracts the spectral content of the induced polarization phenomenon in terms of the Cole-Cole parameters. Synthetic examples and field examples from Denmark showed a significant improvement in the resolution of the parameters that control the induced polarization response when compared to traditional integral chargeability inversion. The quality of the inversion results has been assessed by a complete uncertainty analysis of the model parameters; furthermore, borehole information confirm the outcomes of the field interpretations. With this new accurate code in situ time-domain-induced polarization measurements give access to new applications in environmental and hydrogeophysical investigations, e.g., accurate landfill delineation or on the relation between Cole-Cole and hydraulic parameters.


2021 ◽  
Author(s):  
Alexis Neven ◽  
Pradip Kumar Maurya ◽  
Anders Vest Christiansen ◽  
Philippe Renard

Abstract. Quaternary deposits are complex and heterogeneous. They contain some of the most abundant and extensively used aquifers. In order to improve the knowledge of the spatial heterogeneity of such deposits, we acquired a large (more than 1400 hectares) and dense (20 m spacing) Time Domain ElectroMagnetic (TDEM) dataset in the upper Aare Valley, Switzerland. TDEM is a fast and reliable method to measure the magnetic field directly related to the resistivity of the underground. In this paper, we present the inverted resistivity models derived from this acquisition, and all the necessary data in order to perform different inversions on the processed data (https://doi.org/10.5281/ZENODO.4269887 (Neven et al., 2020)). The depth of investigation ranges between 40 to 120 m depth, with an average data residual contained in the standard deviation of the data. These data can be used for many different purposes: from sedimentological interpretation of quaternary environments in alpine environments, geological and hydrogeological modeling, to benchmarking geophysical inversion techniques.


Sign in / Sign up

Export Citation Format

Share Document