Markov Random Field Models for High-Dimensional Parameters in Simulations of Fluid Flow in Porous Media
- 1 August 2002
- journal article
- Published by Informa UK Limited in Technometrics
- Vol. 44 (3), 230-241
- https://doi.org/10.1198/004017002188618419
Abstract
We give an approach for using flow information from a system of wells to characterize hydrologic properties of an aquifer. In particular, we consider experiments where an impulse of tracer fluid is injected along with the water at the input wells and its concentration is recorded over time at the uptake wells. We focus on characterizing the spatially varying permeability field, which is a key attribute of the aquifer for determining flow paths and rates for a given flow experiment. As is standard for estimation from such flow data, we use complicated subsurface flow code that simulates the fluid flow through the aquifer for a particular well configuration and aquifer specification, in particular the permeability field over a grid. The solution to this ill-posed problem requires that some regularity conditions be imposed on the permeability field. Typically, this regularity is accomplished by specifying a stationary Gaussian process model for the permeability field. Here we use an intrinsically stationary Markov random field, which compares favorably to Gaussian process models and offers some additional flexibility and computational advantages. Our interest in quantifying uncertainty leads us to take a Bayesian approach, using Markov chain Monte Carlo for exploring the high-dimensional posterior distribution. We demonstrate our approach with several examples. We also note that the methodology is general and is not specific to hydrology applications.Keywords
This publication has 14 references indexed in Scilit:
- Bayesian modelling of inseparable space‐time variation in disease riskStatistics in Medicine, 2000
- Partitioning Tracers for Measuring Residual NAPL: Field-Scale Test ResultsJournal of Environmental Engineering, 1998
- Optimal estimation of residual non–aqueous phase liquid saturations using partitioning tracer concentration dataWater Resources Research, 1997
- Markov chain Monte Carlo methods for conditioning a permeability field to pressure dataMathematical Geology, 1997
- Fully Bayesian estimation of Gibbs hyperparameters for emission computed tomography dataIEEE Transactions on Medical Imaging, 1997
- Partitioning Tracer Test for Detection, Estimation, and Remediation Performance Assessment of Subsurface Nonaqueous Phase LiquidsWater Resources Research, 1995
- Bayesian Computation and Stochastic SystemsStatistical Science, 1995
- Simulated Tempering: A New Monte Carlo SchemeEurophysics Letters, 1992
- Bayesian image restoration, with two applications in spatial statisticsAnnals of the Institute of Statistical Mathematics, 1991
- A statistical approach to the inverse problem of aquifer hydrology: 1. TheoryWater Resources Research, 1979