Time domain electromagnetic‐induced polarisation: extracting more induced polarisation information from grounded source time domain electromagnetic data

Electrical induced polarization (EIP) surveys are used to detect chargeable materials in the earth. For interpretation of time domain EIP data a common procedure is to first invert the DC data (electric current on-time) to recover conductivity and then invert the IP data (current off-time) to recover chargeability. This DC-IP procedure assumes that the off-time data are free of secondary EM induction effects. To comply with this, early time data are often discarded, or not recorded. For mid-time data, an EM-decoupling technique, which removes EM induction in the observations, needs to be implemented. Usually responses from a half-space or a layered earth are subtracted. Recent capability in 3D time domain EM (TEM) forward modelling and inversion allows to revisit these procedures. In a TEM-IP survey, a high sampling rate allows early time channels of the EM data to be recorded. The recovery of chargeability then follows a 3-step workflow: (a) invert early time channel TEM data to recover the 3D conductivity; (b) use that conductivity to compute the TEM response at later time channels and subtract this fundamental response from the observations to extract the IP responses, and (c) invert the IP responses to recover a 3D chargeability. This workflow effectively removes EM induction effects in the observations and produces better chargeability and conductivity models compared to conventional approaches. In a synthetic example involving a gradient array we show that the conductivity structure obtained from the early time channel data, which are usually discarded, is superior to that obtained from the steady state DC voltages. This adds a further reason to collect these EM data. This article is protected by copyright. All rights reserved

[1]  Douglas W. Oldenburg,et al.  On recovering distributed IP information from inductive source time domain electromagnetic data , 2016 .

[2]  Douglas W. Oldenburg,et al.  Three-dimensional modeling of IP effects in time-domain electromagnetic data , 2014 .

[3]  Gregory A. Newman,et al.  Induced-polarization effects in time-domain electromagnetic measurements , 1989 .

[4]  James R. Wait,et al.  On electromagnetic coupling “removal” from induced polarization surveys , 1986 .

[5]  Stanley H. Ward,et al.  Mineral discrimination and removal of inductive coupling with multifrequency IP , 1978 .

[6]  Richard S. Smith,et al.  THE TIME-DOMAIN ELECTROMAGNETIC RESPONSE OF POLARIZABLE BODIES: AN APPROXIMATE CONVOLUTION ALGORITHM1 , 1988 .

[7]  P. Veeken,et al.  Benefits of the induced polarization geoelectric method to hydrocarbon exploration , 2009 .

[8]  Esben Auken,et al.  Time-domain-induced polarization: Full-decay forward modeling and 1D laterally constrained inversion of Cole-Cole parameters , 2012 .

[9]  J. Macnae,et al.  Induced polarization in airborne EM , 2012 .

[10]  Michael Commer,et al.  A parallel finite-difference approach for 3D transient electromagnetic modeling with galvanic sources , 2004 .

[11]  D. Mao,et al.  Induced polarization response of porous media with metallic particles — Part 2: Comparison with a broad database of experimental data , 2015 .

[12]  Andrey Tarasov,et al.  On the use of the Cole–Cole equations in spectral induced polarization , 2013 .

[13]  E. Haber,et al.  Inversion of 3D electromagnetic data in frequency and time domain using an inexact all-at-once approach , 2004 .

[14]  Yuval,et al.  Computation of Cole-Cole parameters from IP data , 1997 .

[15]  Lindsey J. Heagy,et al.  SimPEG: An open source framework for simulation and gradient based parameter estimation in geophysical applications , 2015, Comput. Geosci..

[16]  A. Hördt,et al.  Efficient spectral IP-modelling in the time domain , 2006 .

[17]  E. O. Kudryavceva,et al.  Geoelectric modelling with separation between electromagnetic and induced polarization field components , 2009 .

[18]  D. Oldenburg,et al.  Inversion of geophysical data over a copper gold porphyry deposit; a case history for Mt. Milligan , 1997 .

[19]  G. W. Hohmann,et al.  Transient electromagnetic responses of surficial, polarizable patches , 1990 .

[20]  A. Christiansen,et al.  Resolving spectral information from time domain induced polarization data through 2-D inversion , 2013 .

[21]  Douglas W. Oldenburg,et al.  Three dimensional inversion of multisource time domain electromagnetic data , 2013 .

[22]  Douglas W. Oldenburg,et al.  Electromagnetic coupling in frequency-domain induced polarization data: A method for removal , 2001 .

[23]  J. Wynn,et al.  EM coupling, its intrinsic value, its removal and the cultural coupling problem , 1975 .

[24]  D. Oldenburg,et al.  Inversion of induced polarization data , 1994 .

[25]  David L. Alumbaugh,et al.  3D time-domain simulation of electromagnetic diffusion phenomena: A finite-element electric-field approach , 2010 .

[26]  Douglas W. Oldenburg,et al.  3-D inversion of magnetic data , 1996 .

[27]  H. Seigel THE MAGNETIC INDUCED POLARIZATION (MIP) METHOD , 1974 .