Finite-difference method Stokes solver (FDMSS) for 3D pore geometries: Software development, validation and case studies

Abstract Permeability is one of the fundamental properties of porous media and is required for large-scale Darcian fluid flow and mass transport models. Whilst permeability can be measured directly at a range of scales, there are increasing opportunities to evaluate permeability from pore-scale fluid flow simulations. We introduce the free software Finite-Difference Method Stokes Solver (FDMSS) that solves Stokes equation using a finite-difference method (FDM) directly on voxelized 3D pore geometries (i.e. without meshing). Based on explicit convergence studies, validation on sphere packings with analytically known permeabilities, and comparison against lattice-Boltzmann and other published FDM studies, we conclude that FDMSS provides a computationally efficient and accurate basis for single-phase pore-scale flow simulations. By implementing an efficient parallelization and code optimization scheme, permeability inferences can now be made from 3D images of up to 10 9 voxels using modern desktop computers. Case studies demonstrate the broad applicability of the FDMSS software for both natural and artificial porous media.

[1]  Dirk Mallants,et al.  Preferential flow mechanisms identified from staining experiments in forested hillslopes , 2015 .

[2]  Steven J Skerlos,et al.  Modeling of porous filter penneability via image-based stochastic reconstruction of spatial porosity correlations. , 2005, Environmental science & technology.

[3]  G. Debenest,et al.  A Comparison of Various Methods for the Numerical Evaluation of Porous Media Permeability Tensors from Pore-Scale Geometry , 2016, Mathematical Geosciences.

[4]  Kirill M. Gerke,et al.  Description and reconstruction of the soil pore space using correlation functions , 2012, Eurasian Soil Science.

[5]  Robert W. Zimmerman,et al.  Predicting the permeability of sandstone from image analysis of pore structure , 2002 .

[6]  P. Renard,et al.  Laboratory determination of the full permeability tensor , 2001 .

[7]  Qinjun Kang,et al.  Lattice Boltzmann pore-scale model for multicomponent reactive transport in porous media , 2006 .

[8]  H. Gerke,et al.  Dual‐permeability model for flow in shrinking soil with dominant horizontal deformation , 2012 .

[9]  George M. Homsy,et al.  Stokes flow through periodic arrays of spheres , 1982, Journal of Fluid Mechanics.

[10]  Jason J. Williams,et al.  Microstructural Quantification and Property Prediction Using Limited X-ray Tomography Data , 2016 .

[11]  Antti I. Koponen,et al.  Lattice-Boltzmann and finite-difference simulations for the permeability for three-dimensional porous media. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  Bruce D. Jones,et al.  Smoothed particle hydrodynamics and its applications for multiphase flow and reactive transport in porous media , 2015, Computational Geosciences.

[13]  Martys,et al.  Universal scaling of fluid permeability for sphere packings. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[14]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[15]  K. Gerke,et al.  Solution of the stokes equation in three-dimensional geometry by the finite-difference method , 2016 .

[16]  Sidney Diamond,et al.  Mercury porosimetry: An inappropriate method for the measurement of pore size distributions in cement-based materials , 2000 .

[17]  D. Assanis,et al.  Comparison of Pressure-Based and Artificial Compressibility Methods for Solving 3D Steady Incompressible Viscous Flows , 1996 .

[18]  Irina Ginzburg,et al.  Coarse- and fine-grid numerical behavior of MRT/TRT lattice-Boltzmann schemes in regular and random sphere packings , 2015, J. Comput. Phys..

[19]  Alasdair N. Houston,et al.  Adaptive-window indicator kriging: A thresholding method for computed tomography images of porous media , 2013, Comput. Geosci..

[20]  J. Bruining,et al.  Effect of fines migration on oil–water relative permeability during two-phase flow in porous media , 2016 .

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

[22]  Kevin W Eliceiri,et al.  NIH Image to ImageJ: 25 years of image analysis , 2012, Nature Methods.

[23]  F. Stillinger,et al.  Modeling heterogeneous materials via two-point correlation functions: basic principles. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  W. Curtin,et al.  Using microstructure reconstruction to model mechanical behavior in complex microstructures , 2006 .

[25]  T. Patzek,et al.  Pore space morphology analysis using maximal inscribed spheres , 2006 .

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

[27]  S. Torquato,et al.  Reconstructing random media. II. Three-dimensional media from two-dimensional cuts , 1998 .

[28]  Tae Sup Yun,et al.  Upscaling of Navier–Stokes equations in porous media: Theoretical, numerical and experimental approach , 2009 .

[29]  Kirill M. Gerke,et al.  Studying structure and determining permeability of materials based on X-Ray microtomography data (using porous ceramics as an example) , 2015, Inorganic Materials.

[30]  Andrew G. Stack,et al.  The dynamic nature of crystal growth in pores , 2016, Scientific Reports.

[31]  Mark L. Rivers,et al.  Using X-ray computed tomography in hydrology: systems, resolutions, and limitations , 2002 .

[32]  M. Holt,et al.  Proceedings of the Second International Conference on Numerical Methods in Fluid Dynamics , 1971 .

[33]  Quantifying fluorescent tracer distribution in allophanic soils to image solute transport , 2007 .

[34]  M. Blunt,et al.  Pore-scale imaging and modelling , 2013 .

[35]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[36]  R. Helmig.  Multiphase Flow and Transport Processes in the Subsurface: A Contribution to the Modeling of Hydrosystems , 2011 .

[37]  Martin J Blunt,et al.  Pore-network extraction from micro-computerized-tomography images. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  A. Acrivos,et al.  Slow flow through a periodic array of spheres , 1982 .

[39]  C. Torres‐Verdín,et al.  Finite-Difference Approximation for Fluid-Flow Simulation and Calculation of Permeability in Porous Media , 2012, Transport in Porous Media.

[40]  R. Temam,et al.  Navier-Stokes equations: theory and numerical analysis: R. Teman North-Holland, Amsterdam and New York. 1977. 454 pp. US $45.00 , 1978 .

[41]  Faisal Khan,et al.  3D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver , 2011, Journal of Soils and Sediments.

[42]  Kirill M. Gerke,et al.  Improving stochastic reconstructions by weighting correlation functions in an objective function , 2015 .

[43]  P. Delmas,et al.  Application of smoothed particle hydrodynamics (SPH) and pore morphologic model to predict saturated water conductivity from X-ray CT imaging in a silty loam Cambisol , 2015, Geoderma.

[44]  Kirill M. Gerke,et al.  A new way to parameterize hydraulic conductances of pore elements: A step towards creating pore-networks without pore shape simplifications , 2017 .

[45]  Pejman Tahmasebi,et al.  Reconstruction of three-dimensional porous media using a single thin section. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[46]  V. Hejtmánek,et al.  On the measurement of transport parameters of porous solids in permeation and Wicke-Kallenbach cells , 2014 .

[47]  Ali Q. Raeini,et al.  Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method , 2012, J. Comput. Phys..

[48]  Kejian Wu,et al.  SHIFT: an implementation for lattice Boltzmann simulation in low-porosity porous media. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[49]  Charles J Werth,et al.  Evaluation of the effects of porous media structure on mixing-controlled reactions using pore-scale modeling and micromodel experiments. , 2008, Environmental science & technology.

[50]  C. Torres‐Verdín,et al.  A Generalized Finite-Difference Diffusive-Advective (FDDA) Model for Gas Flow in Micro- and Nano-Porous Media , 2009 .

[51]  Kirill M. Gerke,et al.  Tomographic method of studying soil pore space: Current perspectives and results for some Russian soils , 2012, Eurasian Soil Science.

[52]  Dirk Mallants,et al.  Improving pattern reconstruction using directional correlation functions , 2014 .

[53]  M. Blunt,et al.  Numerical study of the effects of particle shape and polydispersity on permeability. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[54]  Pierre M. Adler,et al.  Flow in simulated porous media , 1990 .

[55]  Martin J. Blunt,et al.  Simulation of Flow and Dispersion on Pore-Space Images , 2012 .

[56]  Farzam Javadpour,et al.  Numerical Simulation of Shale-Gas Production: From Pore-Scale Modeling of Slip-Flow, Knudsen Diffusion, and Langmuir Desorption to Reservoir Modeling of Compressible Fluid , 2011 .

[57]  Katja Bachmeier,et al.  Numerical Heat Transfer And Fluid Flow , 2016 .

[58]  P. Adler,et al.  Boundary flow condition analysis for the three-dimensional lattice Boltzmann model , 1994 .

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

[60]  Philippe C. Baveye,et al.  The Operational Significance of the Continuum Hypothesis in the Theory of Water Movement Through Soils and Aquifers , 1984 .

[61]  Adrian Sheppard,et al.  Techniques for image enhancement and segmentation of tomographic images of porous materials , 2004 .

[62]  K. Mecke,et al.  Reconstructing complex materials via effective grain shapes. , 2003, Physical review letters.

[63]  T. Taylor,et al.  Computational methods for fluid flow , 1982 .

[64]  Veerle Cnudde,et al.  Prediction and Evaluation of Time-Dependent Effective Self-diffusivity of Water and Other Effective Transport Properties Associated with Reconstructed Porous Solids , 2015, Transport in Porous Media.

[65]  Farzam Javadpour,et al.  Multiscale, Multiphysics Network Modeling of Shale Matrix Gas Flows , 2013, Transport in Porous Media.

[66]  Dorthe Wildenschild,et al.  Efficiently engineering pore-scale processes: The role of force dominance and topology during nonwetting phase trapping in porous media , 2015 .

[67]  K. Sorbie,et al.  Efficient flow and transport simulations in reconstructed 3D pore geometries , 2010 .

[68]  O. Dikinya,et al.  Decrease in hydraulic conductivity and particle release associated with self-filtration in saturated soil columns , 2008 .

[69]  Min Liu,et al.  Impact of mineralogical heterogeneity on reactive transport modelling , 2017, Comput. Geosci..

[70]  David F. Fouhey,et al.  Representative elementary volume estimation for porosity, moisture saturation, and air‐water interfacial areas in unsaturated porous media: Data quality implications , 2011 .

[71]  B. Ghanbarian,et al.  Permeability in Two‐Component Porous Media: Effective‐Medium Approximation Compared with Lattice‐Boltzmann Simulations , 2016 .

[72]  O. Vizika,et al.  Macroscopic Conductivity of Vugular Porous Media , 2002 .

[73]  John R. Williams,et al.  Characterizing flow in oil reservoir rock using SPH: absolute permeability , 2016 .

[74]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[75]  Martin J. Blunt,et al.  Computations of Absolute Permeability on Micro-CT Images , 2012, Mathematical Geosciences.

[76]  Yang Jiao,et al.  Modeling and characterizing anisotropic inclusion orientation in heterogeneous material via directional cluster functions and stochastic microstructure reconstruction , 2014 .

[77]  C. Noiriel,et al.  Effects of pore-scale precipitation on permeability and flow , 2016 .

[78]  Keijo Mattila,et al.  Using microtomography, image analysis and flow simulations to characterize soil surface seals , 2012, Comput. Geosci..

[79]  Dirk Mallants,et al.  Universal Spatial Correlation Functions for Describing and Reconstructing Soil Microstructure , 2015, PloS one.

[80]  Jan Havelka,et al.  Compression and reconstruction of random microstructures using accelerated lineal path function , 2016, 1601.04359.

[81]  Wolfgang Fink,et al.  Three‐Dimensional Multiphase Segmentation of X‐Ray CT Data of Porous Materials Using a Bayesian Markov Random Field Framework , 2012 .

[82]  Kenneth Stuart Sorbie,et al.  Efficient extraction of networks from three‐dimensional porous media , 2007 .

[83]  Timothy D. Scheibe,et al.  Pore‐scale and multiscale numerical simulation of flow and transport in a laboratory‐scale column , 2015 .

[84]  W. Brent Lindquist,et al.  Image Thresholding by Indicator Kriging , 1999, IEEE Trans. Pattern Anal. Mach. Intell..

[85]  Hans-Jörg Vogel,et al.  Segmentation of X-ray microtomography images of soil using gradient masks , 2010, Comput. Geosci..

[86]  J. Thovert,et al.  Grain reconstruction of porous media: application to a Bentheim sandstone. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[87]  Stig Bakke,et al.  Reconstruction of Berea sandstone and pore-scale modelling of wettability effects , 2003 .

[88]  Lateef T. Akanji,et al.  Finite Element-Based Characterization of Pore-Scale Geometry and Its Impact on Fluid Flow , 2010 .

[89]  Boming Yu,et al.  Developing a new form of permeability and Kozeny-Carman constant for homogeneous porous media by means of fractal geometry , 2008 .

[90]  Yongming Liu,et al.  Numerical investigation of microstructure effect on mechanical properties of bi-continuous and particulate reinforced composite materials , 2016 .

[91]  Wei Wang,et al.  Observer-dependent variability of the thresholding step in the quantitative analysis of soil images and X-ray microtomography data , 2010 .

[92]  M. J. Blunt,et al.  Time-of-Flight Distributions and Breakthrough Curves in Heterogeneous Porous Media Using a Pore-Scale Streamline Tracing Algorithm , 2015, Transport in Porous Media.

[93]  Veerle Cnudde,et al.  Recent progress in X-ray CT as a geosciences tool , 2006 .

[94]  R. E. Larson,et al.  A periodic grain consolidation model of porous media , 1989 .

[95]  Sam L. Butler,et al.  Numerical modeling of fluid and electrical currents through geometries based on synchrotron X-ray tomographic images of reservoir rocks using Avizo and COMSOL , 2014, Comput. Geosci..