Bayesian inversion of geoelectrical resistivity data

Summary. Enormous quantities of geoelectrical data are produced daily and often used for large scale reservoir modelling. To interpret these data requires reliable and efficient inversion methods which adequately incorporate prior information and use realistically complex modelling structures. We use models based on random coloured polygonal graphs as a powerful and flexible modelling framework for the layered composition of the Earth and we contrast our approach with earlier methods based on smooth Gaussian fields. We demonstrate how the reconstruction algorithm may be efficiently implemented through the use of multigrid Metropolis-coupled Markov chain Monte Carlo methods and illustrate the method on a set of field data.

[1]  L. S. Edwards,et al.  A modified pseudosection for resistivity and IP , 1977 .

[2]  W. G. Brown.  Enumeration of Triangulations of the Disk , 1964 .

[3]  R. Barker,et al.  Rapid least-squared inversion of apparent resisitivity pseudosections by a quasi-Newton method , 1996 .

[4]  B. Jacobsen,et al.  Multichannel deconvolution, MCD, in geophysics and helioseismology , 1999 .

[5]  C. Geyer.  Markov Chain Monte Carlo Maximum Likelihood , 1991 .

[6]  B. Jacobsen,et al.  Rapid inversion of 2-D geoelectrical data by multichannel deconvolution , 2001 .

[7]  David E. Boerner,et al.  Fréchet derivatives and single scattering theory , 1989 .

[8]  R. Barker,et al.  Least-squares deconvolution of apparent resistivity pseudosections , 1995 .

[9]  A. Kirsch.  An Introduction to the Mathematical Theory of Inverse Problems , 1996, Applied Mathematical Sciences.

[10]  D. Surgailis,et al.  Markov fields with polygonal realizations , 1989 .

[11]  Goodman,et al.  Multigrid Monte Carlo method. Conceptual foundations. , 1989, Physical review. D, Particles and fields.

[12]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[13]  R. Barker,et al.  Rapid least-squares inversion of apparent resistivity pseudosections , 1994 .

[14]  P. Clifford,et al.  Point-based polygonal models for random graphs , 1993, Advances in Applied Probability.

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

[16]  Judea Pearl,et al.  Probabilistic reasoning in intelligent systems - networks of plausible inference , 1991, Morgan Kaufmann series in representation and reasoning.

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

[18]  C. Geyer,et al.  Annealing Markov chain Monte Carlo with applications to ancestral inference , 1995 .

[19]  Dave Higdon,et al.  A Bayesian approach to characterizing uncertainty in inverse problems using coarse and fine-scale information , 2002, IEEE Trans. Signal Process..

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

[21]  E. Somersalo,et al.  Statistical inversion and Monte Carlo sampling methods in electrical impedance tomography , 2000 .

[22]  Geoff K. Nicholls,et al.  Bayesian image analysis with Markov chain Monte Carlo and coloured continuum triangulation models , 1998 .

[23]  W. Daily,et al.  Electrical resistance tomography experiments at the Oregon Graduate Institute , 1995 .

[24]  Stephen P. Brooks,et al.  A Bayesian approach to crack detection in electrically conducting media , 2001 .

[25]  M. Dabas,et al.  Recent developments in shallow‐depth electrical and electrostatic prospecting using mobile arrays , 1998 .

[26]  P. Green.  Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , 1995 .

[27]  Douglas W. Oldenburg,et al.  Approximate inverse mappings in DC resistivity problems , 1992 .

[28]  K. Sørensen.  Pulled Array Continuous Electrical Profiling , 1996 .

[29]  P. McGillivray,et al.  Forward modeling and inversion of DC resistivity and MMR data , 1992 .

[30]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

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

[32]  Klaus Mosegaard,et al.  Resolution analysis of general inverse problems through inverse Monte Carlo sampling , 1998 .