Issues in characterizing heterogeneity and connectivity in non-multiGaussian media

The performances of kriging, stochastic simulations and sequential self-calibration inversion are assessed when characterizing a non-multiGaussian synthetic 2D braided channel aquifer. The comparison is based on a series of criteria such as the reproduction of the original reference transmissivity or head fields, but also in terms of accuracy of flow and transport (capture zone) forecasts when the flow conditions are modified. We observe that the errors remain large even for a dense data network. In addition some unexpected behaviours are observed when large transmissivity datasets are used. In particular, we observe an increase of the bias with the number of transmissivity data and an increasing uncertainty with the number of head data. This is interpreted as a consequence of the use of an inadequate multiGaussian stochastic model that is not able to reproduce the connectivity of the original field.

[1]  M. Mosley.  Analysis of the effect of changing discharge on channel morphology and instream uses in a Braided River, Ohau River, New Zealand , 1982 .

[2]  Xian-Huan Wen,et al.  Significance of conditioning to piezometric head data for predictions of mass transport in groundwater modeling , 1996 .

[3]  P. Renard,et al.  Dealing with spatial heterogeneity , 2005 .

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

[5]  J. Gómez-Hernández,et al.  To be or not to be multi-Gaussian? A reflection on stochastic hydrogeology , 1998 .

[6]  J. Bear,et al.  On the movement of water bodies injected into aquifers , 1965 .

[7]  J. Valstar,et al.  Inverse modeling of multimodal conductivity distributions , 2006 .

[8]  Sebastien Strebelle,et al.  Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics , 2002 .

[9]  Charles F. Harvey,et al.  When good statistical models of aquifer heterogeneity go bad: A comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields , 2003 .

[10]  L. Hu,et al.  Gradual deformation and iterative calibration of truncated Gaussian simulations , 2001, Petroleum Geoscience.

[11]  S. Redner,et al.  Introduction To Percolation Theory , 2018 .

[12]  J. Caers,et al.  The Probability Perturbation Method: A New Look at Bayesian Inverse Modeling , 2006 .

[13]  Jef Caers,et al.  Comparing the Gradual Deformation with the Probability Perturbation Method for Solving Inverse Problems , 2007 .

[14]  Matthijs van Leeuwen,et al.  Stochastic determination of well capture zones conditioned on regular grids of transmissivity measurements , 2000 .

[15]  J. Gómez-Hernández,et al.  Simulation of Non-Gaussian Transmissivity Fields Honoring Piezometric Data and Integrating Soft and Secondary Information , 1999 .

[16]  Ghislain de Marsily,et al.  Three‐dimensional interference test interpretation in a fractured aquifer using the Pilot Point Inverse Method , 2001 .

[17]  Fabien Cornaton,et al.  Deterministic models of groundwater age, life expectancy and transit time distributions in advective-dispersive systems , 2003 .

[18]  J. W. Mercer,et al.  Contaminant transport in groundwater. , 1992 .

[19]  S. Gorelick,et al.  Heterogeneity in Sedimentary Deposits: A Review of Structure‐Imitating, Process‐Imitating, and Descriptive Approaches , 1996 .

[20]  D. Rosgen.  A classification of natural rivers , 1994 .

[21]  Timothy Scheibe,et al.  An Evaluation of Conditioning Data for Solute Transport Prediction , 2003, Ground water.

[22]  L. Feyen,et al.  Quantifying geological uncertainty for flow and transport modeling in multi-modal heterogeneous formations , 2006 .

[23]  Andres Alcolea,et al.  Pilot points method incorporating prior information for solving the groundwater flow inverse problem , 2006 .

[24]  J. Jaime Gómez-Hernández,et al.  ISIM3D: and ANSI-C three-dimensional multiple indicator conditional simulation program , 1990 .

[25]  O. Banton,et al.  A Comparison of Two Stochastic Inverse Methods in a Field‐Scale Application , 2003, Ground water.

[26]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[27]  M. G. Anderson.  Encyclopedia of hydrological sciences. , 2005 .

[28]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[29]  A. Journel,et al.  The Necessity of a Multiple-Point Prior Model , 2007 .

[30]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric head data—3. Application to the Culebra formation at the Waste Isolation Pilot Plan (WIPP), New Mexico, USA , 1998 .

[31]  Jesús Carrera,et al.  Scale effects in transmissivity , 1996 .

[32]  A. Journel,et al.  Entropy and spatial disorder , 1993 .

[33]  Daniel C. Yoder,et al.  Groundwater Monitoring Network Design Using Minimum Well Density , 1993 .

[34]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

[35]  Giuseppe Gambolati,et al.  Groundwater contour mapping in Venice by stochastic interpolators: 1. Theory , 1979 .

[36]  Xian-Huan Wen,et al.  Numerical modeling of macrodispersion in heterogeneous media: a comparison of multi-Gaussian and non-multi-Gaussian models , 1998 .