On structure-based priors in Bayesian geophysical inversion

Geophysical Journal International - Tập 208 Số 3 - Trang 1342-1358 - 2017
Giulia De Pasquale1, Niklas Linde1
11 Applied and Environmental Geophysics Group, Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland

Tóm tắt

AbstractBayesian methods are extensively used to analyse geophysical data sets. A critical and somewhat overlooked component of high-dimensional Bayesian inversion is the definition of the prior probability density function that describes the joint probability of model parameters before considering available data sets. If insufficient prior information is available about model parameter correlations, then it is tempting to assume that model parameters are uncorrelated. When working with a spatially gridded model representation, this overparametrization leads to posterior realizations with far too much variability to be deemed realistic from a geological perspective. In this study, we introduce a new approach for structure-based prior sampling with Markov chain Monte Carlo that is suitable when only limited prior information is available. We evaluate our method using model structure measures related to standard roughness and damping metrics for l1- and l2-norms. We show that our structure-based prior approach is able to adequately sample the chosen prior distribution of model structure. The usefulness and applicability of the methodology is demonstrated on synthetic and field-based crosshole ground penetrating radar data. We find that our method provides posterior model realizations and statistics that are significantly more satisfactory than those based on underlying assumptions of uncorrelated model parameters or on explicit penalties on model structure within an empirical Bayes framework.

Từ khóa


Tài liệu tham khảo

Bodin, 2009, Seismic tomography with the reversible jump algorithm, Geophys. J. Int., 178, 1411, 10.1111/j.1365-246X.2009.04226.x

Caers, 2004, Multiple-point geostatistics: a quantitative vehicle for integrating geologic analogs into multiple reservoir models, AAPG Mem., 80, 383

Casella, 1985, An introduction to empirical Bayes data analysis, Am. Statisticians, 39, 83, 10.1080/00031305.1985.10479400

Coscia, 2011, 3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water, Geophysics, 76, G49, 10.1190/1.3553003

Cowles, 1996, Markov chain Monte Carlo convergence diagnostics: a comparative review, Am. Statistician Assoc., 91, 883, 10.1080/01621459.1996.10476956

de Groot-Hedlin, 1990, Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data, Geophysics, 55, 1613, 10.1190/1.1442813

Doetsch, 2010, Zonation for 3D aquifer characterization based on joint inversions of multi-method crosshole geophysical data, Geophysics, 75, G53, 10.1190/1.3496476

Engel, 1994, An iterative bandwidth selector for kernel estimation of densities and their derivatives, Nonparametric Stat., 4, 21, 10.1080/10485259408832598

Fraser, 1958, Statistics: An Introduction, 10.1037/14003-000

Gamerman, 2006, Markov Chain Monte Carlo, Stochastic Simulation for Bayesian Inference, 10.1201/9781482296426

Gelman, 1992, Inference from iterative simulation using multiple sequences, Stat. Sci., 7, 457, 10.1214/ss/1177011136

Gelman, 1995, Bayesian Dataset Analysis, 10.1201/9780429258411

Gilks, 1996, Markov Chain Monte Carlo in Practice

Hansen, 2016, Probabilistic integration of geo-information, Integrated Imaging of the Earth: Theory and Applications, 10.1002/9781118929063.ch6

Hastings, 1970, Monte Carlo sampling method using Markov chains and their applications, Biometrika, 57, 97, 10.1093/biomet/57.1.97

Hawkins, 2015, Geophysical imaging using trans-dimensional trees, Geophys. J. Int., 203, 972, 10.1093/gji/ggv326

Huysmans, 2009, Application of multiple-point geostatistics on modeling groundwater flow and transport in cross-bedded aquifer (Belgium), Hydrogeol. J., 17, 1901, 10.1007/s10040-009-0495-2

Jeffreys, 1939, Theory of Probability

Kass, 1995, Bayes factors, J. Am. Stat. Assoc., 90, 773, 10.1080/01621459.1995.10476572

Klotzsche, 2013, 3D characterization of high-permeability zones in a gravel aquifer using 2D crosshole GPR full-waveform inversion and waveguide detection, Geophys. J. Int., 195, 932, 10.1093/gji/ggt275

Kuruwita, 2010, Density estimation using asymmetric kernels Bayes bandwidth with censored data, J. Stat. Plan. Inference, 140, 1765, 10.1016/j.jspi.2010.01.001

Linde, 2013, Distributed soil moisture from crosshole ground-penetrating radar travel times using stochastic inversion, Vadose Zone J., 12, 10.2136/vzj2012.0101

Linde, 2006, Improved hydrogeophysical characterization using joint inversion of crosshole electrical resistance and ground penetrating radar traveltime data, Water Resour. Res., 42, W12404, 10.1029/2006WR005131

Linde, 2015, Geological realism in hydrogeological and geophysical inverse modeling: a review, Adv. Water Resour., 86, 86, 10.1016/j.advwatres.2015.09.019

Lochbühler, 2014, Probabilistic electrical resistivity tomography of a CO2 sequestration analog, J. Appl. Geophys., 107, 80, 10.1016/j.jappgeo.2014.05.013

Lochbühler, 2015, Summary statistics from training images as prior information in probabilistic inversion, Geophys. J. Int., 201, 157, 10.1093/gji/ggv008

Malinverno, 2002, Parsimonious Bayesian Markov chain Monte Carlo inversion in a nonlinear geophysical problem, Geophys. J. Int., 151, 675, 10.1046/j.1365-246X.2002.01847.x

Malinverno, 2004, Expanded uncertainty quantification in inverse problems: hierarchical Bayes and empirical Bayes, Geophysics, 69, 1005, 10.1190/1.1778243

Menke, 1989, International Geophysics Series Volume 45, Geophysical Data Analysis: Discrete Inverse Theory

Metropolis, 1953, Equations of state calculations by fast computing machines, J. Chem. Phys., 21, 1087, 10.1063/1.1699114

Mosegaard, 1995, Monte Carlo sampling of solutions to inverse problems, J. geophys. Res., 1001, 12 431, 10.1029/94JB03097

Pedretti, 2013, An automatic locally-adaptive method to estimate heavily-tailed breakthrough curves from particle distributions, Adv. Water Resour., 59, 52, 10.1016/j.advwatres.2013.05.006

Podvin, 1991, Finite difference computation of travel times in very contrasted velocity models: a massively parallel approach and its associated tools, Geophys. J. Int., 105, 271, 10.1111/j.1365-246X.1991.tb03461.x

Robert, 2010, Introducing Monte Carlo Methods with R, 10.1007/978-1-4419-1576-4

Rosas-Carbajal, 2014, Two-dimensional probabilistic inversion of plane-wave electromagnetic data: methodology, model constraints and joint inversion with electrical resistivity data, Geophys. J. Int., 196, 1508, 10.1093/gji/ggt482

Sambridge, 2002, Monte Carlo methods in geophysical inverse problems, Rev. Geophys., 40, 3.1, 10.1029/2000RG000089

Sambridge, 2006, Trans-dimensional inverse problems, model comparison and the evidence, Geophys. J. Int., 167, 528, 10.1111/j.1365-246X.2006.03155.x

Shimazaki, 2010, Kernel Bandwidth optimization in spike rate estimation, J. Comput. Neurosci., 29, 171, 10.1007/s10827-009-0180-4

Silverman, 1986, Density Estimation for Statistics and Data Analysis, 10.1007/978-1-4899-3324-9

Strebelle, 2002, Conditional simulation of complex geological structures using multiple-point statistics, Math. Geol., 34, 1, 10.1023/A:1014009426274

Tarantola, 2005, Inverse Problem Theory and Model Parameter Estimation, 10.1137/1.9780898717921

Tarantola, 1982, Inverse problem = quest for information, Geophysics, 50, 150

Yang, 2013, Searching for efficient Markov chain Monte Carlo proposal kernels, Proc. Natl. Acad. Sci. USA, 110, 19 370, 10.1073/pnas.1311790110

Thông tin
Thông tin xuất bản

Geophysical Journal International

Tập 208 Số 3

1342-1358

Thông tin tác giả