Bayesian inference of the Cole–Cole parameters from time‐ and frequency‐domain induced polarization

The inversion of induced-polarization parameters is important in the characterization of the frequency electrical response of porous rocks. A Bayesian approach is developed to invert these parameters assuming the electrical response is described by a Cole–Cole model in the time or frequency domain. We show that the Bayesian approach provides a better analysis of the uncertainty associated with the parameters of the Cole–Cole model compared with more conventional methods based on the minimization of a cost function using the least-squares criterion. This is due to the strong non-linearity of the inverse problem and non-uniqueness of the solution in the time domain. The Bayesian approach consists of propagating the information provided by the measurements through the model and combining this information with a priori knowledge of the data. Our analysis demonstrates that the uncertainty in estimating the Cole–Cole model parameters from induced-polarization data is much higher for measurements performed in the time domain than in the frequency domain. Our conclusion is that it is very difficult, if not impossible, to retrieve the correct value of the Cole–Cole parameters from time-domain induced-polarization data using standard least-squares methods. In contrast, the Cole–Cole parameters can be more correctly inverted in the frequency domain. These results are also valid for other models describing the induced-polarization spectral response, such as the Cole–Davidson or power law models.

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

[2]  Tamaz Chelidze,et al.  Electrical spectroscopy of porous rocks: a review—I. Theoretical models , 1999 .

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

[4]  Andreas Kemna,et al.  Field Applications of Complex Resistivity Tomography , 1999 .

[5]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[6]  R. Hilfer Analytical representations for relaxation functions of glasses , 2002 .

[7]  Daizhan Cheng,et al.  On the adequacy of identified cole--cole models , 2003 .

[8]  J. Schön,et al.  Physical Properties of Rocks: Fundamentals and Principles of Petrophysics , 1996 .

[9]  L. Cathles,et al.  Electrical conductivity in shaly sands with geophysical applications , 1998 .

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

[11]  C. Dias,et al.  Analytical model for a polarizable medium at radio and lower frequencies , 1972 .

[12]  Kurt M. Strack,et al.  Society of Exploration Geophysicists , 2007 .

[13]  A. Niemeijer,et al.  Compaction of quartz sands by pressure solution using a Cole-Cole distribution of relaxation times , 2006 .

[14]  M. Taherian,et al.  Measurement of dielectric response of water‐saturated rocks , 1990 .

[15]  H. Seigel,et al.  Some case histories of source discrimination using time-domain spectral IP , 1997 .

[16]  A. Tarantola,et al.  Inverse problems = Quest for information , 1982 .

[17]  J. Wait A PHENOMENOLOGICAL THEORY OF OVERVOLTAGE FOR METALLIC PARTICLES , 1959 .

[18]  L. M. M.-T. Theory of Probability , 1929, Nature.

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

[20]  J. Hinderer,et al.  The search for weak harmonic signals in a spectrum with application to gravity data , 1995 .

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

[22]  K. Cole,et al.  Dispersion and Absorption in Dielectrics I. Alternating Current Characteristics , 1941 .

[23]  D. Oldenburg,et al.  Regularized in version of spectral IP parameters from complex resistivity data , 1998 .

[24]  Daizhan Cheng,et al.  Direct inversion of the apparent complex-resistivity spectrum , 2001 .

[25]  David P. Lesmes,et al.  Electrical‐hydraulic relationships observed for unconsolidated sediments , 2002 .

[26]  S. R. Jaggar,et al.  Forward and Inverse Cole-Cole modelling in the analysis of frequency domain electrical impedance data , 1988 .

[27]  J. Tombs,et al.  The feasibility of making spectral IP measurements in the time domain , 1981 .

[28]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

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

[30]  W. R. Sill,et al.  Electrical properties of artificial clay-bearing sandstone , 1982 .

[31]  Albert Tarantola,et al.  Monte Carlo sampling of solutions to inverse problems , 1995 .

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

[33]  J. Hinderer,et al.  Bayesian estimation of the free core nutation parameters from the analysis of precise tidal gravity data , 2000 .

[34]  H. Robain,et al.  A rapid electrical sounding method The «three-point» method: a Bayesian approach , 2001 .

[35]  H. Seigel Mathematical formulation and type curves for induced polarization , 1959 .

[36]  C. Schlumberger,et al.  Étude sur la prospection électrique du sous-sol , 1920 .

[37]  A. Tarantola Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[38]  Roy H Wilkens,et al.  Physical Properties of Rocks: Fundamentals and Principles of Petrophysics , 1997 .

[39]  Andreas Kemna,et al.  Crosshole IP imaging for engineering and environmental applications , 2004 .

[40]  Yanzhong Luo,et al.  Theory and Application of Spectral Induced Polarization , 1998 .