A BME solution of the inverse problem for saturated groundwater flow

Abstract.In most real-world hydrogeologic situations, natural heterogeneity and measurement errors introduce major sources of uncertainty in the solution of the inverse problem. The Bayesian Maximum Entropy (BME) method of modern geostatistics offers an efficient solution to the inverse problem by first assimilating various physical knowledge bases (hydrologic laws, water table elevation data, uncertain hydraulic resistivity measurements, etc.) and then producing robust estimates of the subsurface variables across space. We present specific methods for implementing the BME conceptual framework to solve an inverse problem involving Darcy’s law for subsurface flow. We illustrate one of these methods in the case of a synthetic one-dimensional case study concerned with the estimation of hydraulic resistivity conditioned on soft data and hydraulic head measurements. The BME framework processes the physical knowledge contained in Darcy’s law and generates accurate estimates of hydraulic resistivity across space. The optimal distribution of hard and soft data needed to minimize the associated estimation error at a specified sampling cost is determined.

[1]  George Christakos,et al.  Modern Spatiotemporal Geostatistics , 2000 .

[2]  G. Christakos A Bayesian/maximum-entropy view to the spatial estimation problem , 1990 .

[3]  Marc L. Serre,et al.  Modern geostatistics: computational BME analysis in the light of uncertain physical knowledge – the Equus Beds study , 1999 .

[4]  Alberto Guadagnini,et al.  Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains: 2. Computational examples , 1999 .

[5]  S. P. Neuman,et al.  Prediction of steady state flow in nonuniform geologic media by conditional moments: Exact nonlocal , 1993 .

[6]  J. Bear Dynamics of Fluids in Porous Media , 1975 .

[7]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

[8]  D. D'Or,et al.  Application of the BME approach to soil texture mapping , 2001 .

[9]  George Christakos,et al.  Efficient BME estimation of subsurface hydraulic properties using measurements of water table elevation in unidirectional flow , 2003 .

[10]  George Christakos,et al.  On the assimilation of uncertain physical knowledge bases: Bayesian and non-Bayesian techniques. , 2002 .

[11]  Alberto Guadagnini,et al.  Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion , 2003 .

[12]  George Christakos,et al.  Random Field Models in Earth Sciences , 1992 .