A comparison between Gauss-Newton and Markov-chain Monte Carlo–based methods for inverting spectral induced-polarization data for Cole-Cole parameters

Abstract

We develop a Bayesian model to invert spectral induced polarization (SIP) data for ColeCole parameters using Markov chain Monte Carlo (MCMC) sampling methods. We compare the performance of the MCMC based stochastic method with an iterative Gauss-Newton based deterministic method for Cole-Cole parameter estimation through inversion of synthetic and laboratory SIP data. The Gauss-Newton based method can provide an optimal solution for given objective functions under constraints, but the obtained optimal solution generally depends on the choice of initial values and the estimated uncertainty information is often inaccurate or insufficient. In contrast, the MCMC based inversion method provides extensive global information on unknown parameters, such as the marginal probability distribution functions, from which we can obtain better estimates and tighter uncertainty bounds of the parameters than with the deterministic method. Additionally, the results obtained with the MCMC method are independent of the choice of initial values. Because the MCMC based method does not explicitly offer single optimal solution for given objective functions, the deterministic and stochastic methods can complement each other. For example, the stochastic method can first be used to obtain the means of the unknown parameters by starting from an arbitrary set of initial values and the deterministic method can then be initiated using the means as starting values to obtain the optimal estimates of the Cole-Cole parameters.3

Calcite precipitation dominates the electrical signatures of zero valent iron columns under simulated field conditions

“…We inverted the electrical data using a deterministic algorithm based on a leastsquares approach with Marquardt regularization (Kemna, 2000). Although this inversion method offers reasonable parameter estimates satisfying the purpose of this study, recently developed stochastic inversion algorithms (using a Bayesian framework) have shown improved parameter estimates and tighter uncertainty bounds of parameters (Ghorbani et al, 2007;Chen et al, 2008).…”

Impact of changes in grain size and pore space on the hydraulic conductivity and spectral induced polarization response of sand

“…The inversion procedure is based on the work of Mosegaard and Tarantola (1995). For a detailed description of the application to SIP data, and a comparison to a deterministic approach of fitting the Cole-Cole model, please see Chen et al (2008). We considered uniform prior distributions for all parameters except for the low-frequency asymptote of the resistivity ρ 0 , which we set to the resistivity value observed at the lowest measurement frequency of 1 mHz.…”

Application of time domain induced polarization to the mapping of lithotypes in a landfill site

“…The induced polarization method has been used increasingly in environmental investigations, because IP measurements are very sensitive to the low frequency capacitive properties of rocks and soils (Chen et al, 2008). These properties are associated with diffusion-controlled polarization processes that occur at the mineral-fluid interface (Slater and Lesmes, 2002).…”