Lattice Boltzmann Simulations of Fluid Flow in Continental Carbonate Reservoir Rocks and in Upscaled Rock Models Generated with Multiple-Point Geostatistics

Microcomputed tomography (μCT) and Lattice Boltzmann Method (LBM) simulations were applied to continental carbonates to quantify fluid flow. Fluid flow characteristics in these complex carbonates with multiscale pore networks are unique and the applied method allows studying their heterogeneity and anisotropy. 3D pore network models were introduced to single-phase flow simulations in Palabos, a software tool for particle-based modelling of classic computational fluid dynamics. In addition, permeability simulations were also performed on rock models generated with multiple-point geostatistics (MPS). This allowed assessing the applicability of MPS in upscaling high-resolution porosity patterns into large rock models that exceed the volume limitations of the μCT. Porosity and tortuosity control fluid flow in these porous media. Micro- and mesopores influence flow properties at larger scales in continental carbonates. Upscaling with MPS is therefore necessary to overcome volume-resolution problems of CT scanning equipment. The presented LBM-MPS workflow is applicable to other lithologies, comprising different pore types, shapes, and pore networks altogether. The lack of straightforward porosity-permeability relationships in complex carbonates highlights the necessity for a 3D approach. 3D fluid flow studies provide the best understanding of flow through porous media, which is of crucial importance in reservoir modelling.

[1]  J. Soete,et al.  Shrub morpho-types as indicator for the water flow energy - Tivoli travertine case (Central Italy) , 2017 .

[2]  A. Foubert,et al.  Multiscale approach to (micro)porosity quantification in continental spring carbonate facies: Case study from the Cakmak quarry (Denizli, Turkey) , 2016 .

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

[4]  Alessandro Comunian,et al.  Building a training image with Digital Outcrop Models , 2015 .

[5]  Alain Burgisser,et al.  Characterization of juvenile pyroclasts from the Kos Plateau Tuff (Aegean Arc): insights into the eruptive dynamics of a large rhyolitic eruption , 2009 .

[6]  Xin Wang,et al.  The construction of carbonate digital rock with hybrid superposition method , 2013 .

[7]  S. Kele,et al.  Stable isotope geochemical study of Pamukkale travertines: New evidences of low-temperature non-equilibrium calcite-water fractionation , 2011 .

[8]  T. Ford,et al.  A review of tufa and travertine deposits of the world , 1996 .

[9]  Yu-Wei Chang,et al.  Comparison of the Quaternary travertine sites in the Denizli extensional basin based on their depositional and geochemical data , 2013 .

[10]  A. Foubert,et al.  Comparative study of the Pleistocene Cakmak quarry (Denizli Basin, Turkey) and modern Mammoth Hot Springs deposits (Yellowstone National Park, USA) , 2017 .

[11]  S. Bakke,et al.  Pore network modelling on carbonate: a comparative study of different micro-CT network extraction methods , 2008 .

[12]  D. Bourlès,et al.  Dating the Homo erectus bearing travertine from Kocabaş (Denizli, Turkey) at at least 1.1 Ma , 2014 .

[13]  Daniel H. Rothman,et al.  Lattice-Boltzmann simulations of flow through Fontainebleau sandstone , 1995 .

[14]  Kenneth Stuart Sorbie,et al.  Representation of multiscale heterogeneity via multiscale pore networks , 2013 .

[15]  Veerle Cnudde,et al.  Multi-scale, micro-computed tomography-based pore network models to simulate drainage in heterogeneous rocks , 2015 .

[16]  Orestis Malaspinas,et al.  Synchrotron X-ray microtomography and lattice Boltzmann simulations of gas flow through volcanic pumices , 2010 .

[17]  D. Wildenschild,et al.  X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems , 2013 .

[18]  D. D. Marais,et al.  Depositional facies and aqueous-solid geochemistry of travertine-depositing hot springs (Angel Terrace, Mammoth Hot Springs, Yellowstone National Park, U.S.A.). , 2001, Journal of sedimentary research. Section A, Sedimentary petrology and processes : an international journal of SEPM.

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

[20]  J. Soete,et al.  Shrub and pore type classification: Petrography of travertine shrubs from the Ballık-Belevi area (Denizli, SW Turkey) , 2017 .

[21]  V. Cnudde,et al.  Software tools for quantification of X-ray microtomography at the UGCT , 2007 .

[22]  Wei Shyy,et al.  Lattice Boltzmann Method for 3-D Flows with Curved Boundary , 2000 .

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

[24]  Thomas Tran,et al.  Improving variogram reproduction on dense simulation grids , 1994 .

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

[26]  David Hodgetts,et al.  Laser scanning and digital outcrop geology in the petroleum industry: A review , 2013 .

[27]  J. Soete,et al.  Geobody architecture, genesis and petrophysical characteristics of the Budakalász travertines, Buda Hills (Hungary) , 2017 .

[28]  Veerle Cnudde,et al.  HECTOR: A 240kV micro-CT setup optimized for research , 2013 .

[29]  A. Foubert,et al.  Acoustic properties in travertines and their relation to porosity and pore types , 2015 .

[30]  A. Reed,et al.  Advances in Computed Tomography for Geomaterials: GeoX 2010 , 2010 .

[31]  Cédric Gommes,et al.  Practical methods for measuring the tortuosity of porous materials from binary or gray‐tone tomographic reconstructions , 2009 .

[32]  P. Pedley,et al.  An Introduction to Fluid Dynamics , 1968 .

[33]  Tuanfeng Zhang,et al.  MPS-Driven Digital Rock Modeling and Upscaling , 2015, Mathematical Geosciences.

[34]  P. Ronchi,et al.  Continental carbonates as a hydrocarbon reservoir, an analog case study from the travertine of Saturnia, Italy , 2015 .

[35]  Wim Degruyter,et al.  Controls on magma permeability in the volcanic conduit during the climactic phase of the Kos Plateau Tuff eruption (Aegean Arc) , 2009 .

[36]  S. Claes.  Pore classification and upscaling strategy in travertine reservoir rocks , 2015 .

[37]  Philippe Renard,et al.  A practical guide to performing multiple-point statistical simulations with the Direct Sampling algorithm , 2013, Comput. Geosci..

[38]  Mark Bentley,et al.  Reservoir Model Design: A Practitioner's Guide , 2015 .

[39]  V. Wright.  Lacustrine carbonates in rift settings: the interaction of volcanic and microbial processes on carbonate deposition , 2012 .

[40]  S. Schröder,et al.  Algal–microbial carbonates of the Namibe Basin (Albian, Angola): implications for microbial carbonate mound development in the South Atlantic , 2015 .

[41]  RW3D.m: three-dimensional random walk program for the calculation of the diffusivities in porous media , 2002 .

[42]  R. Swennen,et al.  Geobody architecture of continental carbonates: “Gazda” travertine quarry (Süttő, Gerecse Hills, Hungary) , 2017 .

[43]  Philippe Renard,et al.  Constraining distance-based multipoint simulations to proportions and trends , 2015, Environ. Model. Softw..

[44]  M. Blunt,et al.  Prediction of permeability for porous media reconstructed using multiple-point statistics. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[45]  R. Riding,et al.  Hot‐spring travertine facies and sequences, Late Pleistocene, Rapolano Terme, Italy , 1998 .

[46]  Veerle Cnudde,et al.  A three-dimensional classification for mathematical pore shape description in complex carbonate reservoir rocks , 2016, Mathematical Geosciences.

[47]  Y. Nakashima,et al.  DMAP.m: A Mathematica® program for three-dimensional mapping of tortuosity and porosity of porous media. , 2004 .

[48]  Veerle Cnudde,et al.  High-resolution X-ray computed tomography in geosciences: A review of the current technology and applications , 2013 .

[49]  Rudy Swennen,et al.  Application of microfocus computed tomography in carbonate reservoir characterization: Possibilities and limitations , 2008 .

[50]  Dirk Mallants,et al.  Universal Stochastic Multiscale Image Fusion: An Example Application for Shale Rock , 2015, Scientific Reports.

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

[52]  Jonas Latt.  Choice of units in lattice Boltzmann simulations , 2008 .

[53]  A. Saller,et al.  Presalt stratigraphy and depositional systems in the Kwanza Basin, offshore Angola , 2016 .

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

[55]  Abdellah Kharicha,et al.  3D simulation of interdendritic flow through a Al-18wt.%Cu structure captured with X-ray microtomography , 2012 .

[56]  Andre G. Journel,et al.  Geostatistics: Roadblocks and Challenges , 1993 .

[57]  W. Janna,et al.  Introduction to Fluid Mechanics , 2012 .

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

[59]  J. Soete,et al.  Sedimentology, three‐dimensional geobody reconstruction and carbon dioxide origin of Pleistocene travertine deposits in the Ballık area (south‐west Turkey) , 2015 .

[60]  Samir Khatib,et al.  Études stratigraphique, sédimentologique et paléomagnétique des travertins de Kocabas ¸, Bassin de Denizli, Anatolie, Turquie, contenant des restes fossiles quaternaires , 2014 .