Direct simulation of non-additive properties on unstructured grids

Abstract Uncertainties related to permeability heterogeneity can be estimated using geostatistical simulation methods. Usually, these methods are applied on regular grids with cells of constant size, whereas unstructured grids are more flexible to honor geological structures and offer local refinements for fluid-flow simulations. However, cells of different sizes require to account for the support dependency of permeability statistics (support effect). This paper presents a novel workflow based on the power averaging technique. The averaging exponent ω is estimated using a response surface calibrated from numerical upscaling experiments. Using spectral turning bands, permeability is simulated on points in each unstructured cell, and later averaged with a local value of ω that depends on the cell size and shape. The method is illustrated on a synthetic case. The simulation of a tracer experiment is used to compare this novel geostatistical simulation method with a conventional approach based on a fine scale Cartesian grid. The results show the consistency of both the simulated permeability fields and the tracer breakthrough curves. The computational cost is much lower than the conventional approach based on a pressure-solver upscaling.

[1]  Julián M. Ortiz,et al.  Two approaches to direct block-support conditional co-simulation , 2011, Comput. Geosci..

[2]  A. Mantoglou,et al.  The Turning Bands Method for simulation of random fields using line generation by a spectral method , 1982 .

[3]  Maxine Akhurst,et al.  Risk Assessment-Led Characterisation of the SiteChar UK North Sea Site for the Geological Storage of CO2 , 2015 .

[4]  L. Durlofsky,et al.  Unstructured 3D gridding and upscaling for coarse modelling of geometrically complex reservoirs , 2005, Petroleum Geoscience.

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

[6]  José Elguero,et al.  STATISTICAL ANALYSIS OF SOLVENT SCALES. PART 1 , 1997 .

[7]  Vincent C. Tidwell,et al.  Laboratory method for investigating permeability upscaling , 1997 .

[8]  Hans Wackernagel,et al.  Change-of-Support Models on Irregular Grids for Geostatistical Simulation , 2016, Mathematical Geosciences.

[9]  Christian Lantuéjoul,et al.  Geostatistical Simulation: Models and Algorithms , 2001 .

[10]  M. Shinozuka,et al.  Digital simulation of random processes and its applications , 1972 .

[11]  C. Deutsch,et al.  Upscaling permeability to unstructured grids using the multipoint flux approximation , 2012 .

[12]  Dongxiao Zhang,et al.  Efficient analytical upscaling method for elliptic equations in three-dimensional heterogeneous anisotropic media , 2020 .

[13]  Benoit Noetinger,et al.  Influence of the Permeability Anisotropy Ratio on Large-Scale Properties of Heterogeneous Reservoirs , 1993 .

[14]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.

[15]  J. Jaime Gómez-Hernández,et al.  Stochastic analysis of three-dimensional hydraulic conductivity upscaling in a heterogeneous tropical soil , 2018, Computers and Geotechnics.

[16]  B. Lévy,et al.  Voronoi grids conforming to 3D structural features , 2014, Computational Geosciences.

[17]  Whitney Trainor-Guitton,et al.  Uncertainty and risk evaluation during the exploration stage of geothermal development: A review , 2019, Geothermics.

[18]  Benoît Nœtinger,et al.  Use of the Fourier-Laplace transform and of diagrammatical methods to interpret pumping tests in heterogeneous reservoirs , 1998 .

[19]  Jean-Marc Gratien,et al.  Reservoir Simulator Runtime Enhancement Based on a Posteriori Error Estimation Techniques , 2016 .

[20]  Benoit Noetinger,et al.  Calculation of Internodal Transmissivities in Finite Difference Models of Flow in Heterogeneous Porous Media , 1995 .

[21]  Martin J. Mlacnik,et al.  Locally Streamline-Pressure-Potential-Based PEBI Grids , 2003 .

[22]  Toby Aiken,et al.  Geological storage of CO2 in saline aquifers—A review of the experience from existing storage operations , 2010 .

[23]  G. Matheron.  Forecasting Block Grade Distributions: The Transfer Functions , 1976 .

[24]  P. Indelman,et al.  Effective permittivity of log-normal isotropic random media , 1995 .

[25]  R. Ababou,et al.  Implementation of the three‐dimensional turning bands random field generator , 1989 .

[26]  Carl W. Gable,et al.  Large-scale numerical simulation of groundwater flow and solute transport in discretely-fractured crystalline bedrock , 2011 .

[27]  A. De Wit,et al.  Correlation structure dependence of the effective permeability of heterogeneous porous media , 1995 .

[28]  Benoit Noetinger,et al.  Multiscale Description and Upscaling of Fluid Flow in Subsurface Reservoirs , 2004 .

[29]  X. Emery,et al.  On a continuous spectral algorithm for simulating non-stationary Gaussian random fields , 2018, Stochastic Environmental Research and Risk Assessment.

[30]  Benoît N∄tinger,et al.  The effective permeability of a heterogeneous porous medium , 1994 .

[31]  A. Boschan,et al.  On the multiscale characterization of effective hydraulic conductivity in random heterogeneous media: A historical survey and some new perspectives , 2020 .

[32]  I. A. Antonov,et al.  An economic method of computing LPτ-sequences , 1979 .

[33]  M. Claeys-Bruno,et al.  Construction of space-filling designs using WSP algorithm for high dimensional spaces , 2012 .

[34]  C. R. Dietrich,et al.  A Simple and Efficient Space Domain Implementation of the Turning Bands Method , 1995 .

[35]  Ivar Aavatsmark,et al.  Discretization on Unstructured Grids for Inhomogeneous, Anisotropic Media. Part I: Derivation of the Methods , 1998, SIAM J. Sci. Comput..

[36]  Guillaume Caumon,et al.  High Resolution Geostatistics on Coarse Unstructured Flow Grids , 2005 .

[37]  A. Boschan,et al.  Scale Dependence of Effective Hydraulic Conductivity Distributions in 3D Heterogeneous Media: A Numerical Study , 2012, Transport in Porous Media.

[38]  Philippe Renard,et al.  A fast algorithm for the estimation of the equivalent hydraulic conductivity of heterogeneous media , 2000 .

[39]  Gedeon Dagan,et al.  Higher-order correction of effective permeability of heterogeneous isotropic formations of lognormal conductivity distribution , 1993 .

[40]  G. Marsily,et al.  Some current methods to represent the heterogeneity of natural media in hydrogeology , 1998 .

[41]  Clayton V. Deutsch,et al.  Geostatistical Assignment of Reservoir Properties on Unstructured Grids , 2002 .

[42]  Clayton V. Deutsch,et al.  Power Averaging for Block Effective Permeability , 1986 .

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

[44]  A. J. Desbarats,et al.  Spatial averaging of hydraulic conductivity in three-dimensional heterogeneous porous media , 1992 .

[45]  Benoit Noetinger,et al.  Permeability Averaging for Well Tests in 3D Stochastic Reservoir Models , 1996 .

[46]  Christophe Preux,et al.  About the Use of Quality Indicators to Reduce Information Loss When Performing Upscaling , 2016 .

[47]  Rachid Ababou,et al.  Random Porous Media Flow on Large 3-D Grids: Numerics, Performance, and Application to Homogenization , 1996 .

[48]  Darius Mottaghy,et al.  Reducing temperature uncertainties by stochastic geothermal reservoir modelling , 2010 .

[49]  X. Emery,et al.  An improved spectral turning-bands algorithm for simulating stationary vector Gaussian random fields , 2016, Stochastic Environmental Research and Risk Assessment.

[50]  Christian Lantuéjoul,et al.  TBSIM: A computer program for conditional simulation of three-dimensional Gaussian random fields via the turning bands method , 2006, Comput. Geosci..