Using Multiple-Point Geostatistics for Tracer Test Modeling in a Clay-Drape Environment with Spatially Variable Conductivity and Sorption Coefficient

This study investigates the effect of fine-scale clay drapes on tracer transport. A tracer test was performed in a sandbar deposit consisting of cross-bedded sandy units intercalated with many fine-scale clay drapes. The heterogeneous spatial distribution of the clay drapes causes a spatially variable hydraulic conductivity and sorption coefficient. A fluorescent tracer (sodium naphthionate) was injected in two injection wells and ground water was sampled and analyzed from five pumping wells. To determine (1) whether the fine-scale clay drapes have a significant effect on the measured concentrations and (2) whether application of multiple-point geostatistics can improve interpretation of tracer tests in media with complex geological heterogeneity, this tracer test is analyzed with a local three-dimensional ground-water flow and transport model in which fine-scale sedimentary heterogeneity is modeled using multiple-point geostatistics. To reduce memory needs and calculation time for the multiple-point geostatistical simulation step, this study uses the technique of direct multiple-point geostatistical simulation of edge properties. Instead of simulating pixel values, model cell edge properties indicating the presence of irregularly shaped surfaces are simulated using multiple-point geostatistical simulations. Results of a sensitivity analysis show under which conditions clay drapes have a significant effect on the concentration distribution. Calibration of the model against measured concentrations from the tracer tests reduces the uncertainty on the clay-drape parameters. The calibrated model shows which features of the breakthrough curves can be attributed to the geological heterogeneity of the aquifer and which features are caused by other processes.

[1]  J. Caers,et al.  Conditional Simulation with Patterns , 2007 .

[2]  H. Vereecken,et al.  Analysis of solute transport in a heterogeneous aquifer: the Krauthausen field experiment , 2000 .

[3]  A. Dassargues,et al.  Delineation of groundwater protection zones based on tracer tests and transport modeling in alluvial sediments , 1998 .

[4]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[5]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[6]  R. Zhang,et al.  Applied Contaminant Transport Modeling: Theory and Practice , 1991 .

[7]  Zhiming Lu,et al.  Coupling Basin‐ and Site‐Scale Inverse Models of the Espanola Aquifer , 2003, Ground water.

[8]  D. Mikeš Sampling procedure for small-scale heterogeneities (crossbedding) for reservoir modelling , 2006 .

[9]  R. M. Srivastava,et al.  Multivariate Geostatistics: Beyond Bivariate Moments , 1993 .

[10]  Hydrogeological modeling of radionuclide transport in low permeability media: a comparison between Boom Clay and Ypresian Clay , 2006 .

[11]  A. Dassargues,et al.  Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium) , 2009 .

[12]  Alexandre Boucher,et al.  A SGeMS code for pattern simulation of continuous and categorical variables: FILTERSIM , 2008, Comput. Geosci..

[13]  Fred J. Molz,et al.  Analysis and Interpretation of Single-Well Tracer Tests in Stratified Aquifers , 1985 .

[14]  Georg Teutsch,et al.  Forced and natural gradient tracer tests in a highly heterogeneous porous aquifer: instrumentation and measurements , 1994 .

[15]  P. A. Hsieh,et al.  Documentation of a computer program to simulate horizontal-flow barriers using the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model , 1993 .

[16]  Kenneth Stuart Sorbie,et al.  Immiscible flow behaviour in laminated and cross-bedded sandstones , 1993 .

[17]  C. Leibundgut,et al.  Tracers in Hydrology , 2009 .

[18]  A. Dassargues,et al.  Relationship between sedimentary features and permeability at different scales in the Brussels Sands , 2012 .

[19]  Steven M. Gorelick,et al.  Geological modeling of submeter scale heterogeneity and its influence on tracer transport in a fluvial aquifer , 2010 .

[20]  Alain Dassargues,et al.  Direct Multiple-Point Geostatistical Simulation of Edge Properties for Modeling Thin Irregularly Shaped Surfaces , 2011 .

[21]  Orhun Aydin,et al.  Image transforms for determining fit-for-purpose complexity of geostatistical models in flow modeling , 2013, Computational Geosciences.

[22]  Alexandre Boucher,et al.  Applied Geostatistics with SGeMS: A User's Guide , 2009 .

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

[24]  E. Sudicky,et al.  Spatial Variability of Strontium Distribution Coefficients and Their Correlation With Hydraulic Conductivity in the Canadian Forces Base Borden Aquifer , 1991 .

[25]  P. Renard,et al.  Three-dimensional high resolution fluvio-glacial aquifer analog - Part 2: Geostatistical modeling , 2011 .

[26]  Philippe Renard,et al.  3D multiple-point statistics simulation using 2D training images , 2012, Comput. Geosci..

[27]  Gregoire Mariethoz,et al.  The Direct Sampling method to perform multiple‐point geostatistical simulations , 2010 .

[28]  R. Houthuys Vergelijkende studie van de afzettingsstruktuur van getijdenzanden uit het eoceen en van de huidige vlaamse banken , 1990 .

[29]  Jef Caers,et al.  Geological modelling and history matching of multi-scale flow barriers in channelized reservoirs: methodology and application , 2011 .

[30]  Jef Caers,et al.  Multiple-point Geostatistics: A Quantitative Vehicle for Integrating Geologic Analogs into Multiple Reservoir Models , 2004 .

[31]  R. Dimitrakopoulos,et al.  High-order Stochastic Simulation of Complex Spatially Distributed Natural Phenomena , 2010 .

[32]  M. Hill,et al.  Development and evaluation of a local grid refinement method for block-centered finite-difference groundwater models using shared nodes , 2002 .

[33]  Philippe Renard,et al.  Three-dimensional high resolution fluvio-glacial aquifer analog: Part 1: Field study , 2011 .

[34]  L. Y. Hu,et al.  Multiple‐point geostatistics for modeling subsurface heterogeneity: A comprehensive review , 2008 .

[35]  Gq,et al.  Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution , 2002 .

[36]  G. Mariéthoz,et al.  An Improved Parallel Multiple-point Algorithm Using a List Approach , 2011 .

[37]  Wolfgang Kinzelbach,et al.  3D-Groundwater Modeling with PMWIN , 2003 .

[38]  P. Adler,et al.  Multiple-Point Statistics for Modeling Facies Heterogeneities in a Porous Medium: The Komadugu-Yobe Alluvium, Lake Chad Basin , 2011 .

[39]  Okke Batelaan,et al.  GIS-based recharge estimation by coupling surface–subsurface water balances , 2007 .

[40]  P. Corbett,et al.  Detailed analysis of probe permeameter and interval pressure transient test permeability measurements in a heterogeneous reservoir , 2002, Petroleum Geoscience.

[41]  P. Renard,et al.  Comparison of three geostatistical methods for hydrofacies simulation: a test on alluvial sediments , 2012, Hydrogeology Journal.

[42]  Steffen Mehl,et al.  Comparison of Local Grid Refinement Methods for MODFLOW , 2006, Ground water.

[43]  Yuhong Liu,et al.  Using the Snesim program for multiple-point statistical simulation , 2006, Comput. Geosci..

[44]  Alain Dassargues,et al.  Modeling the effect of clay drapes on pumping test response in a cross-bedded aquifer using multiple-point geostatistics , 2012 .

[45]  R. Houthuys A sedimentary model of the Brussels Sands, Eocene, Belgium , 2010 .

[46]  Alain Dassargues,et al.  Relating small-scale sedimentary structures and permeability in a cross-bedded aquifer , 2008 .

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

[48]  Paul Switzer,et al.  Filter-Based Classification of Training Image Patterns for Spatial Simulation , 2006 .

[49]  A. Dassargues,et al.  Application of Multiple-Point Geostatistics on Modelling Groundwater Flow and Transport in a Cross-Bedded Aquifer , 2010 .

[50]  H. Reineck,et al.  Depositional Sedimentary Environments, With Reference to Terrigenous Clastics , 1973 .

[51]  Louis J. Durlofsky,et al.  A multi-resolution workflow to generate high-resolution models constrained to dynamic data , 2011 .

[52]  Brian J. Willis,et al.  Quantitative Outcrop Data for Flow Simulation , 2000 .

[53]  L. Stright Modeling, Upscaling, and History Matching Thin, Irregularly-Shaped Flow Barriers: A Comprehensive Approach for Predicting Reservoir Connectivity , 2006 .

[54]  S. Brouyère,et al.  Integrating geophysical and tracer test data for accurate solute transport modelling in heterogeneous porous media , 2002 .

[55]  A. Journel,et al.  Reservoir Modeling Using Multiple-Point Statistics , 2001 .