From characterisation of pore-structures to simulations of pore-scale fluid flow and the upscaling of permeability using microtomography: A case study of heterogeneous carbonates

Abstract We propose a workflow to provide detailed characterisation of microstructures, simulations of permeability at micro-scale, and to analyse the upscaling of permeability. The workflow is fully tested using two very different microtomographic datasets of carbonate samples with strong heterogeneity. The characterisation of microstructures includes not only routine parameters but also anisotropy of the pore-structure. The critical pore diameter is also obtained through morphological implementations. The method of stochastic analysis of extended local porosity theory is used to determine the sizes of representative volume elements (RVE). We demonstrate criteria for determining the size of the RVE and show a sample that satisfies the criterion. We also discuss a sample that does not suffice the RVE criterion, because of the strong anisotropy. Permeabilities are computed using Lattice-Boltzmann (LB) simulations on the RVE and compared with laboratory measurements. Techniques and procedures are used to extract scaling parameters for both of the samples including: 1) percolation threshold — by using a shrinking/expanding algorithm; 2) crossover length — by analysing the mass density; 3) the critical exponent of correlation length — by using the finite-size scaling scheme; and 4) the critical exponent of permeability — by running LB simulations on a series of derivative models close to the percolation threshold. Our results of critical exponents are different from the previous studies and the mass density distributions are irregular. This pilot study provides new information on the relationships between microstructures and permeability of natural rocks with complex microstructures. The study also reveals the scaling parameters. Our results clearly put into question whether natural rocks can be idealised by classical theoretical solutions. A robust workflow for embracing a computational approach of microstructures may provide a solution.

[1]  Alison Ord,et al.  Theoretical and numerical analyses of chemical‐dissolution front instability in fluid‐saturated porous rocks , 2008 .

[2]  B. Hobbs,et al.  The influence of faulting on host-rock permeability, fluid flow and ore genesis of gold deposits: a theoretical 2D numerical model , 2003 .

[3]  Yanhua Zhang,et al.  Numerical modelling of orogenic processes and gold mineralisation in the southeastern part of the Yilgarn Craton, Western Australia , 2002 .

[4]  B. Hobbs,et al.  Theoretical analyses of nonaqueous phase liquid dissolution‐induced instability in two‐dimensional fluid‐saturated porous media , 2010 .

[5]  Xiaodong Jia,et al.  Combining X-ray microtomography with computer simulation for analysis of granular and porous materials , 2010 .

[6]  Alison Ord,et al.  Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity , 2003 .

[7]  Chongbin Zhao,et al.  Fluids in geological processes — The present state and future outlook , 2010 .

[8]  Ali Karrech,et al.  Multiscale coupling and multiphysics approaches in earth sciences: Theory , 2013 .

[9]  Christoph H. Arns,et al.  Computation of linear elastic properties from microtomographic images: Methodology and agreement between theory and experiment , 2002 .

[10]  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.

[11]  Alison Ord,et al.  Convective and Advective Heat Transfer in Geological Systems , 2008 .

[12]  Liang-ming Liu,et al.  Numerical Modeling of Coupled Geodynamical Processes and Its Role in Facilitating Predictive Ore Discovery: An Example from Tongling, China , 2005 .

[13]  Chongbin Zhao,et al.  Dynamic and Transient Infinite Elements: Theory and Geophysical, Geotechnical and Geoenvironmental Applications , 2016 .

[14]  Paul W. Cleary,et al.  Comparison of permeability of model porous media between SPH and LB , 2012 .

[15]  D. Stauffer,et al.  Self similarity and correlations in percolation , 1983 .

[16]  B. Hobbs,et al.  Potential effects of upward throughflow on thermal structure models within the continental lithospheric mantle-crust , 2008 .

[17]  S. Yun,et al.  Reaction path modeling of hydrogeochemical evolution of groundwater in granitic bedrocks, South Korea , 2012 .

[18]  Feng,et al.  Differences between lattice and continuum percolation transport exponents. , 1985, Physical review letters.

[19]  Philippe Montarnal,et al.  Impact of spatial heterogeneities on oxygen consumption in sediments: Experimental observations and 2D numerical modeling , 2012 .

[20]  A. Katz,et al.  Prediction of rock electrical conductivity from mercury injection measurements , 1987 .

[21]  Jie Liu,et al.  Multiscale coupling and multiphysics approaches in earth sciences : applications , 2013 .

[22]  Christoph H. Arns,et al.  Pore Scale Characterization of Carbonates Using X-Ray Microtomography , 2005 .

[23]  Chongbin Zhao,et al.  Structural controls on fluid flow and related mineralization in the Xiangshan uranium deposit, Southern China , 2006 .

[24]  A. Bannari,et al.  Assessment of soil contamination around an abandoned mine in a semi-arid environment using geochemistry and geostatistics: Pre-work of geochemical process modeling with numerical models , 2013 .

[25]  Salih Muhammad Awadh,et al.  Geochemical exploration using surveys of spring water, hydrocarbon and gas seepage, and geobotany for determining the surface extension of Abu-Jir Fault Zone in Iraq: A new way for determining geometrical shapes of computational simulation models , 2013 .

[26]  Klaus Regenauer-Lieb,et al.  Application of percolation theory to microtomography of structured media: percolation threshold, critical exponents, and upscaling. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  B. Hobbs,et al.  Fault-related dilation, permeability enhancement, fluid flow and mineral precipitation patterns: numerical models , 2008 .

[28]  Mark L. Rivers,et al.  An introduction to the application of X-ray microtomography to the three-dimensional study of igneous rocks , 2012 .

[29]  G. G. Pereira,et al.  NUMERICAL PORE-SCALE MODELING OF THREE-PHASE FLUID FLOW : COMPARISON BETWEEN SIMULATION AND EXPERIMENT , 1999 .

[30]  Feng,et al.  Transport properties of continuum systems near the percolation threshold. , 1987, Physical review. B, Condensed matter.

[31]  Zhenhua Chai,et al.  Non-Darcy flow in disordered porous media: A lattice Boltzmann study , 2010 .

[32]  Joel Koplik,et al.  Conductivity and permeability of rocks , 1984 .

[33]  P. Gow,et al.  Copper‐gold mineralisation in New Guinea: Numerical modelling of collision, fluid flow and intrusion‐related hydrothermal systems , 2002 .

[34]  Alison Ord,et al.  Fundamentals of Computational Geoscience: Numerical Methods and Algorithms , 2009 .

[35]  Chongbin Zhao,et al.  Geodynamic constraints on orebody localization in the Anqing orefield, China: Computational modeling and facilitating predictive exploration of deep deposits , 2011 .

[36]  Muhammad Sahimi,et al.  Non-linear and non-local transport processes in heterogeneous media: from long-range correlated percolation to fracture and materials breakdown , 1998 .

[37]  Chongbin Zhao,et al.  Finite element modeling of pore-fluid flow in the Dachang ore district, Guangxi, China: Implications for hydrothermal mineralization , 2011 .

[38]  R. Ketcham.  Three-dimensional grain fabric measurements using high-resolution X-ray computed tomography , 2005 .

[39]  J. Almeida,et al.  Managing borehole samples of unequal lengths to construct a high-resolution mining model of mineral grades zoned by geological units , 2013 .

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

[41]  Yanhua Zhang,et al.  Geodynamic modelling of the Century deposit, Mt Isa Province, Queensland , 2002 .

[42]  C. L. Farmer.  Upscaling: a review , 2002 .

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

[44]  Xiao‐Ming Li,et al.  Apatite fission track age of Mesozoic sandstones from Beipiao basin, eastern China: Implications for basin provenance and tectonic evolution , 2003 .

[45]  Andrew Squelch,et al.  Improved estimates of percolation and anisotropic permeability from 3‐D X‐ray microtomography using stochastic analyses and visualization , 2009 .

[46]  J. Boon.  The Lattice Boltzmann Equation for Fluid Dynamics and Beyond , 2003 .

[47]  D. d'Humières,et al.  Multiple–relaxation–time lattice Boltzmann models in three dimensions , 2002, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[48]  Andre Peters,et al.  Prediction of capillary hysteresis in a porous material using lattice-Boltzmann methods and comparison to experimental data and a morphological pore network model , 2008 .

[49]  Cass T. Miller,et al.  An evaluation of lattice Boltzmann schemes for porous medium flow simulation , 2006 .

[50]  Jens Harting,et al.  Quantitative analysis of numerical estimates for the permeability of porous media from lattice-Boltzmann simulations , 2010, 1005.2530.

[51]  W. V. Pinczewski,et al.  3D imaging and flow characterization of the pore space of carbonate core samples , 2006 .

[52]  Matthew D. Jackson,et al.  Detailed physics, predictive capabilities and macroscopic consequences for pore-network models of multiphase flow. , 2002 .

[53]  S. Baxter,et al.  A framework for stochastic mechanics , 2006 .

[54]  Chongbin Zhao,et al.  Coupled geodynamics in the formation of Cu skarn deposits in the Tongling–Anqing district, China: Computational modeling and implications for exploration , 2010 .

[55]  R. Hilfer.  Review on Scale Dependent Characterization of the Microstructure of Porous Media , 2001, cond-mat/0105458.

[56]  Jan D. Miller,et al.  Saturated flow characteristics in column leaching as described by LB simulation , 2005 .

[57]  Yanhua Zhang,et al.  Application of coupled deformation, fluid flow, thermal and chemical modelling to predictive mineral exploration , 2000 .

[58]  Huilin Xing,et al.  Three-dimensional finite element simulation of large-scale nonlinear contact friction problems in deformable rocks , 2008 .

[59]  Paul W. Cleary,et al.  SPH modelling of fluid at the grain level in a porous medium , 2011 .

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

[61]  Chongbin Zhao,et al.  Numerical analysis and simulation experiment of lithospheric thermal structures in the South China Sea and the Western Pacific , 2009 .

[62]  Hilfer,et al.  Local-porosity theory for flow in porous media. , 1992, Physical review. B, Condensed matter.