Inversion of generalized relaxation time distributions with optimized damping parameter

Retrieving the Relaxation Time Distribution (RDT), the Grains Size Distribution (GSD) or the Pore Size Distribution (PSD) from low-frequency impedance spectra is a major goal in geophysics. The ``Generalized RTD'' generalizes parametric models like Cole-Cole and many others, but remains tricky to invert since this inverse problem is ill-posed. We propose to use generalized relaxation basis function (for instance by decomposing the spectra on basis of generalized Cole-Cole relaxation elements instead of the classical Debye basis) and to use the L-curve approach to optimize the damping parameter required to get smooth and realistic inverse solutions. We apply our algorithm to three examples, one synthetic and two real data sets, and the program includes the possibility of converting the RTD into GSD or PSD by choosing the value of the constant connecting the relaxation time to the characteristic polarization size of interest. A high frequencies (typically above 1 kHz), a dielectric term in taken into account in the model. The code is provided as an open Matlab source as a supplementary file associated with this paper. (C) 2014 Elsevier B.V. All rights reserved.

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

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

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

[4]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

[5]  Michael J. Tompkins,et al.  On the topography of the cost functional in linear and nonlinear inverse problems , 2012 .

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

[7]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

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

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

[10]  D. Miranda,et al.  Determination of Cole–Cole parameters using only the real part of electrical impedivity measurements , 2008, Physiological measurement.

[11]  S. Havriliak,et al.  A complex plane analysis of α‐dispersions in some polymer systems , 2007 .

[12]  E. Barsoukov,et al.  Impedance spectroscopy : theory, experiment, and applications , 2005 .

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

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

[15]  Nicolas Florsch,et al.  Quantification of slag heap volumes and masses through the use of induced polarization: application to the Castel-Minier site , 2011 .

[16]  Andrey Tarasov,et al.  Relaxation time distribution from time domain induced polarization measurements , 2007 .

[17]  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 .

[18]  M. Tong,et al.  Determining capillary-pressure curve, pore-size distribution, and permeability from induced polarization of shaley sand , 2006 .

[19]  E. Warburg,et al.  Ueber das Verhalten sogenannter unpolarisirbarer Elektroden gegen Wechselstrom , 1899 .

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

[21]  G. Schwarz.  A THEORY OF THE LOW-FREQUENCY DIELECTRIC DISPERSION OF COLLOIDAL PARTICLES IN ELECTROLYTE SOLUTION1,2 , 1962 .

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

[23]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[24]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

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

[26]  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 .

[27]  A. Revil.  On charge accumulation in heterogeneous porous rocks under the influence of an external electric field , 2013 .

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

[29]  Vivek Agarwal,et al.  Total Variation Regularization and L-curve method for the selection of regularization parameter , 2003 .