Modeling Fine‐Scale Geological Heterogeneity—Examples of Sand Lenses in Tills

Sand lenses at various spatial scales are recognized to add heterogeneity to glacial sediments. They have high hydraulic conductivities relative to the surrounding till matrix and may affect the advective transport of water and contaminants in clayey till settings. Sand lenses were investigated on till outcrops producing binary images of geological cross-sections capturing the size, shape and distribution of individual features. Sand lenses occur as elongated, anisotropic geobodies that vary in size and extent. Besides, sand lenses show strong non-stationary patterns on section images that hamper subsequent simulation. Transition probability (TP) and multiple-point statistics (MPS) were employed to simulate sand lens heterogeneity. We used one cross-section to parameterize the spatial correlation and a second, parallel section as a reference: it allowed testing the quality of the simulations as a function of the amount of conditioning data under realistic conditions. The performance of the simulations was evaluated on the faithful reproduction of the specific geological structure caused by sand lenses. Multiple-point statistics offer a better reproduction of sand lens geometry. However, two-dimensional training images acquired by outcrop mapping are of limited use to generate three-dimensional realizations with MPS. One can use a technique that consists in splitting the 3D domain into a set of slices in various directions that are sequentially simulated and reassembled into a 3D block. The identification of flow paths through a network of elongated sand lenses and the impact on the equivalent permeability in tills are essential to perform solute transport modeling in the low-permeability sediments.

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

[2]  P. Renard,et al.  Calculating equivalent permeability: a review , 1997 .

[3]  P. Huggenberger,et al.  Ground-probing radar as a tool for heterogeneity estimation in gravel deposits: advances in data-processing and facies analysis , 1994 .

[4]  Deli Chen,et al.  Comparison of Sequential Indicator Simulation and Transition Probability Indicator Simulation Used to Model Clay Content in Microscale Surface Soil , 2009 .

[5]  M. Houmark‐Nielsen Signature and timing of the Kattegat Ice Stream: onset of the Last Glacial Maximum sequence at the southwestern margin of the Scandinavian Ice Sheet , 2003 .

[6]  P. Jussel,et al.  Transport modeling in heterogeneous aquifers: 1. Statistical description and numerical generation of gravel deposits , 1994 .

[7]  R. Sidle,et al.  Transport and biodegradation of creosote compounds in clayey till, a field experiment , 2000 .

[8]  W. Cowan,et al.  Moraine formation in northwestern Ontario: product of subglacial fluvial and glaciolacustrine sedimentation , 1990 .

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

[10]  R. Bossi,et al.  Preferential flow and pesticide transport in a clay‐rich till: Field, laboratory, and modeling analysis , 2002 .

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

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

[13]  L. Murdoch,et al.  Field characterization of vertical bromide transport in a fractured glacial till , 2007 .

[14]  濱田 吉郎,et al.  信頼性制御と Frobenius Norm 最適化 , 1997 .

[15]  Tavi Murray,et al.  Characterisation of glacial sediments using geophysical methods for groundwater source protection , 2005 .

[16]  Andrew Richard Gardiner,et al.  Best practice stochastic facies modeling from a channel-fill turbidite sandstone analog (the Quarry outcrop, Eocene Ainsa basin, northeast Spain) , 2006 .

[17]  L. McKay,et al.  Distribution, Origin, and Hydraulic Influence of Fractures in a Clay-rich Glacial Deposit , 1995 .

[18]  S. Kullback,et al.  Information Theory and Statistics , 1959 .

[19]  Philippe Renard,et al.  Connectivity metrics for subsurface flow and transport , 2013 .

[20]  R. Reese Geostatistics for Environmental Scientists , 2001 .

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

[22]  Социальная педагогика,et al.  Society for Sedimentary Geology , 2010 .

[23]  P. Witherspoon,et al.  Porous media equivalents for networks of discontinuous fractures , 1982 .

[24]  P. Goovaerts Geostatistics in soil science: state-of-the-art and perspectives , 1999 .

[25]  A. Dreimanis Small to medium-sized glacitectonic structures in till and in its substratum and their comparison with mass movement structures , 1993 .

[26]  Graham E. Fogg,et al.  Multi-scale alluvial fan heterogeneity modeled with transition probability geostatistics in a sequence stratigraphic framework , 1999 .

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

[28]  R. Lark,et al.  Geostatistics for Environmental Scientists , 2001 .

[29]  R. Sidle,et al.  Mass transport and scale-dependent hydraulic tests in a heterogeneous glacial till–sandy aquifer system , 2001 .

[30]  Richard R. Jones,et al.  Characterization of fluvial architectural elements using a three-dimensional outcrop data set: Escanilla braided system, South-Central Pyrenees, Spain , 2007 .

[31]  R. Gillham,et al.  Field experiments in a fractured clay till: 1. Hydraulic conductivity and fracture aperture , 1993 .

[32]  R. D. Barker,et al.  The use of electrical resistivity tomography in deriving local-scale models of recharge through superficial deposits , 2009 .

[33]  G. Pedersen,et al.  Lacustrine sedimentation and tectonics: an example from the Weichselian at Lønstrup Klint, Denmark , 1997 .

[34]  Steven F. Carle,et al.  Geologic heterogeneity and a comparison of two geostatistical models: Sequential Gaussian and transition probability-based geostatistical simulation , 2007 .

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

[36]  N. Eyles The role of meltwater in glacial processes , 2006 .

[37]  D. Seifert,et al.  Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time. , 2004, Journal of Contaminant Hydrology.

[38]  Solomon Kullback,et al.  Information Theory and Statistics , 1970, The Mathematical Gazette.

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

[40]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[41]  Steven F. Carle,et al.  CONDITIONAL SIMULATION OF HYDROFACIES ARCHITECTURE: A TRANSITION PROBABILITY/MARKOV APPROACH1 , 1998 .

[42]  Marcel G. Schaap,et al.  Percolation Theory for Flow in Porous Media , 2006 .

[43]  Steven F. Carle,et al.  Connected-network paradigm for the alluvial aquifer system , 2000 .

[44]  U. Asprion,et al.  Aquifer architecture analysis using ground-penetrating radar: Triassic and Quaternary examples (S. Germany) , 1997 .

[45]  Steven F. Carle,et al.  Geologically based model of heterogeneous hydraulic conductivity in an alluvial setting , 1998 .

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

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

[48]  N. Robinson,et al.  Impact of permeable conduits on solute transport in aquitards: Mathematical models and their application , 2007 .

[49]  Katherine A. Klise,et al.  Exploring solute transport and streamline connectivity using lidar‐based outcrop images and geostatistical representations of heterogeneity , 2009 .

[50]  P. Bertran,et al.  Facies and microfacies of slope deposits , 1999 .

[51]  Denis Allard On the Connectivity of Two Random Set Models: The Truncated Gaussian and the Boolean , 1993 .

[52]  Timothy T. Eaton,et al.  On the importance of geological heterogeneity for flow simulation , 2006 .

[53]  P. Binning,et al.  Modeling multi-component transport and enhanced anaerobic dechlorination processes in a single fracture-clay matrix system. , 2010, Journal of contaminant hydrology.

[54]  J. Sminchak,et al.  Indicator geostatistical analysis of sand interconnections within a till , 1996 .

[55]  P. L. Bjerg,et al.  Characterization of sand lenses embedded in tills , 2012 .

[56]  F. G. Bastante,et al.  Comparison of indicator kriging, conditional indicator simulation and multiple-point statistics used to model slate deposits , 2008 .

[57]  G. Fogg,et al.  Transition probability-based indicator geostatistics , 1996 .

[58]  G. Fraser,et al.  Late Wisconsinan Proglacial Sedimentation Along the West Chicago Moraine in Northeastern Illinois , 1982 .

[59]  E. Phillips Micromorphology of a debris flow deposit: evidence of basal shearing, hydrofracturing, liquefaction and rotational deformation during emplacement , 2006 .

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

[61]  L. Wassenaar,et al.  Characterizing the hydrogeology of a complex clay-rich aquitard system using detailed vertical profiles of the stable isotopes of water , 2004 .

[62]  Robert W. Ritzi,et al.  Behavior of indicator variograms and transition probabilities in relation to the variance in lengths of hydrofacies , 2000 .

[63]  J. Boyce,et al.  Evaluation of heterogeneity and field-scale groundwater flow regime in a leaky till aquitard , 2001 .

[64]  P. Gravesen,et al.  Fractures and biopores in weichselian clayey till aquitards at Flakkebjerg, Denmark , 1999 .