Bayesian inference of spectral induced polarization parameters for laboratory complex resistivity measurements of rocks and soils

Abstract Spectral induced polarization (SIP) measurements are now widely used to infer mineralogical or hydrogeological properties from the low-frequency electrical properties of the subsurface in both mineral exploration and environmental sciences. We present an open-source program that performs fast multi-model inversion of laboratory complex resistivity measurements using Markov-chain Monte Carlo simulation. Using this stochastic method, SIP parameters and their uncertainties may be obtained from the Cole-Cole and Dias models, or from the Debye and Warburg decomposition approaches. The program is tested on synthetic and laboratory data to show that the posterior distribution of a multiple Cole-Cole model is multimodal in particular cases. The Warburg and Debye decomposition approaches yield unique solutions in all cases. It is shown that an adaptive Metropolis algorithm performs faster and is less dependent on the initial parameter values than the Metropolis-Hastings step method when inverting SIP data through the decomposition schemes. There are no advantages in using an adaptive step method for well-defined Cole-Cole inversion. Finally, the influence of measurement noise on the recovered relaxation time distribution is explored. We provide the geophysics community with a open-source platform that can serve as a base for further developments in stochastic SIP data inversion and that may be used to perform parameter analysis with various SIP models.

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

[2]  Lee Slater,et al.  On the relationship between induced polarization and surface area in metal-sand and clay-sand mixtures , 2006 .

[3]  A. Tarasov,et al.  Time domain spectral induced polarization of disseminated electronic conductors: Laboratory data analysis through the Debye decomposition approach , 2013 .

[4]  James R. Clark,et al.  Constraints on the Genesis of the Archean Oxidized, Intrusion-Related Canadian Malartic Gold Deposit, Quebec, Canada , 2014 .

[5]  Andreas Kemna,et al.  A comparison between Gauss-Newton and Markov-chain Monte Carlo–based methods for inverting spectral induced-polarization data for Cole-Cole parameters , 2008 .

[6]  Nicolas Florsch,et al.  Direct estimation of the distribution of relaxation times from induced-polarization spectra using a Fourier transform analysis , 2012 .

[7]  Nicolas Florsch,et al.  Determination of permeability from spectral induced polarization in granular media , 2010 .

[8]  Charles L. Bérubé,et al.  Structural setting for Canadian Malartic style of gold mineralization in the Pontiac Subprovince, south of the Cadillac Larder Lake Deformation Zone, Québec, Canada , 2017 .

[9]  C. Torres‐Verdín,et al.  Interfacial polarization of disseminated conductive minerals in absence of redox-active species — Part 1: Mechanistic model and validation , 2016 .

[10]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[11]  David Huard,et al.  PyMC: Bayesian Stochastic Modelling in Python. , 2010, Journal of statistical software.

[12]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[13]  A. Gelman,et al.  Weak convergence and optimal scaling of random walk Metropolis algorithms , 1997 .

[14]  H. Haario,et al.  An adaptive Metropolis algorithm , 2001 .

[15]  Stephen P. Brooks,et al.  Markov chain Monte Carlo method and its application , 1998 .

[16]  B. Carlin,et al.  Markov Chain Monte Carlo conver-gence diagnostics: a comparative review , 1996 .

[17]  Andrey Tarasov,et al.  Induced polarization of disseminated electronically conductive minerals: a semi-empirical model , 2015 .

[18]  Lee Slater,et al.  Spectral induced polarization (SIP) response of mine tailings. , 2015, Journal of contaminant hydrology.

[19]  J. Rosenthal,et al.  Coupling and Ergodicity of Adaptive Markov Chain Monte Carlo Algorithms , 2007, Journal of Applied Probability.

[20]  Nicolas Florsch,et al.  Induced polarization response of porous media with metallic particles — Part 1: A theory for disseminated semiconductors , 2015 .

[21]  David P. Lesmes,et al.  Dielectric spectroscopy of sedimentary rocks , 2001 .

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

[23]  Nicolas Florsch,et al.  Spectral induced polarization porosimetry , 2014 .

[24]  Nicolas Florsch,et al.  Bayesian inference of the Cole–Cole parameters from time‐ and frequency‐domain induced polarization , 2007 .

[25]  David P. Lesmes,et al.  Inversion for dielectric relaxation spectra , 1994 .

[26]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

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

[28]  André Revil,et al.  Is it the grain size or the characteristic pore size that controls the induced polarization relaxation time of clean sands and sandstones? , 2012 .

[29]  Nicolas Florsch,et al.  Inversion of generalized relaxation time distributions with optimized damping parameter , 2014 .

[30]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[31]  J. Schön,et al.  Technical Note: A Relation between the Quadrature Component of Electrical Conductivity and the Specific Surface Area of Sedimentary Rocks , 1991 .

[32]  J. Wong,et al.  An electrochemical model of the induced‐polarization phenomenon in disseminated sulfide ores , 1979 .

[33]  Markku Peltoniemi,et al.  Spectral IP studies of Finnish ore prospects , 1992 .

[34]  Andreas Weller,et al.  A new approach to fitting induced-polarization spectra , 2008 .

[35]  Andrew Binley,et al.  Textural controls on low-frequency electrical spectra of porous media , 2010 .

[36]  H. Vinegar,et al.  Induced polarization of shaly sands , 1984 .

[37]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[38]  Lee Slater,et al.  On the estimation of specific surface per unit pore volume from induced polarization: A robust empirical relation fits multiple data sets , 2010 .

[39]  Andreas Kemna,et al.  Estimating the spatiotemporal distribution of geochemical parameters associated with biostimulation using spectral induced polarization data and hierarchical Bayesian models , 2012 .

[40]  A. Ghorbani,et al.  Complex conductivity of water-saturated packs of glass beads. , 2008, Journal of colloid and interface science.

[41]  Andreas Kemna,et al.  Relationship between low-frequency electrical properties and hydraulic permeability of low-permeability sandstones , 2010 .

[42]  C. Dias,et al.  Developments In A Model To Describe Low-Frequency Electrical Polarization Of Rocks - (I): Theory , 1999 .

[43]  A. Kemna,et al.  Debye decomposition of time-lapse spectral induced polarisation data , 2016, Comput. Geosci..

[44]  André Revil,et al.  Connecting complex conductivity spectra to mercury porosimetry of sedimentary rocks , 2016 .

[45]  A. Tarantola,et al.  Generalized Nonlinear Inverse Problems Solved Using the Least Squares Criterion (Paper 1R1855) , 1982 .

[46]  J. Rosenthal,et al.  Optimal scaling for various Metropolis-Hastings algorithms , 2001 .

[47]  Ian M. Johnson,et al.  Spectral induced-polarization parameters as determined through time-domain measurements , 1984 .

[48]  Andrew Binley,et al.  Laboratory SIP signatures associated with oxidation of disseminated metal sulfides. , 2013, Journal of contaminant hydrology.

[49]  Andreas Weller,et al.  Estimating permeability of sandstone samples by nuclear magnetic resonance and spectral-induced polarization , 2010 .

[50]  C. Torres‐Verdín,et al.  Interfacial polarization of disseminated conductive minerals in absence of redox-active species — Part 2: Effective electrical conductivity and dielectric permittivity , 2016 .

[51]  Heikki Haario,et al.  Adaptive proposal distribution for random walk Metropolis algorithm , 1999, Comput. Stat..

[52]  Andrew Binley,et al.  Markov-chain Monte Carlo estimation of distributed Debye relaxations in spectral induced polarization , 2012 .

[53]  Heikki Vanhala,et al.  MAPPING OIL‐CONTAMINATED SAND AND TILL WITH THE SPECTRAL INDUCED POLARIZATION (SIP) METHOD , 1997 .

[54]  Andrew Binley,et al.  Relationship between spectral induced polarization and hydraulic properties of saturated and unsaturated sandstone , 2005 .