On structure-based priors in Bayesian geophysical inversion
Tóm tắt
Bayesian 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
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
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
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
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
Strebelle, 2002, Conditional simulation of complex geological structures using multiple-point statistics, Math. Geol., 34, 1, 10.1023/A:1014009426274
Tarantola, 1982, Inverse problem = quest for information, Geophysics, 50, 150