Effect of X-ray CT resolution on the quality of permeability computation for granular soils: definition of a criterion based on morphological properties
Miriam Patricia Ortega Ramírez A , Laurent Oxarango![https://orcid.org/0000-0003-1062-4303](/media/client/orcid_16x16.png)
![https://orcid.org/0000-0001-9668-5822](/media/client/orcid_16x16.png)
A Université Grenoble Alpes, CNRS, IRD, Grenoble-INP, IGE, F-38000 Grenoble, France.
B Instituto de Ciencias Aplicadas y Tecnologia, Universidad Nacional Autónoma de Mexico, Mexico, D.F., Mexico.
C Corresponding author. Email: laurent.oxarango@univ-grenoble-alpes.fr
Soil Research 57(6) 589-600 https://doi.org/10.1071/SR18189
Submitted: 16 July 2018 Accepted: 24 May 2019 Published: 16 July 2019
Abstract
In this study, the quality of soil permeability estimation based on computational fluid dynamics is discussed. Two types of three-dimensional geometries were considered: an image of Fontainebleau sand obtained from X-ray computed micro-tomography and a virtual pack of spheres. Numerical methods such as finite difference or lattice Boltzmann can conveniently use the image voxels as computational mesh elements. In this framework, the image resolution is directly associated with quality of the numerical computation. A higher resolution should promote both a better morphological description and discretisation. However, increasing the resolution may prevent the studied volume from being representative. Here, each sample was scaled and analysed at five resolutions. The dependence of soil properties with respect to the image resolution is discussed. As resolution decreased, the permeability and specific surface values tended to diverge from the reference value. This deterioration could be attributed to the shift of the pore size distribution towards badly resolved pores in the voxelised geometry. As long as granular soils are investigated, the volume fraction of pores smaller than six voxels in diameter should not exceed 50% to ensure the validity of permeability computation. In addition, based on an analysis of flow distribution, the volume fraction of pores smaller than four voxels should not exceed 25% in order to limit the flow rate occurring in badly discretised pores under 10%. For the Fontainebleau sand and virtual pack of spheres, the maximum voxel size meeting this criterion corresponded to 1/14 and 1/20 of the mean grain size respectively.
Additional keywords: digital rock physics, image resolution, soil morphology, soil permeability.
References
Ahmadi MA (2015) Connectionist approach estimates gas–oil relative permeability in petroleum reservoirs: application to reservoir simulation. Fuel 140, 429–439.| Connectionist approach estimates gas–oil relative permeability in petroleum reservoirs: application to reservoir simulation.Crossref | GoogleScholarGoogle Scholar |
Al-Raoush RI, Willson CS (2005) Extraction of physically realistic pore network properties from three-dimensional synchrotron X-ray microtomography images of unconsolidated porous media systems. Journal of Hydrology 300, 44–64.
| Extraction of physically realistic pore network properties from three-dimensional synchrotron X-ray microtomography images of unconsolidated porous media systems.Crossref | GoogleScholarGoogle Scholar |
Al-Raoush RI, Thompson K, Willson CS (2003) Comparison of network generation techniques for unconsolidated porous media. Soil Science Society of America Journal 67, 1687–1700.
| Comparison of network generation techniques for unconsolidated porous media.Crossref | GoogleScholarGoogle Scholar |
Bear J (2013). ‘Dynamics of fluids in porous media.’ (Dover publications Inc.: New York).
Borujeni AT, Lane NM, Thompson K, Tyagi M (2013) Effects of image resolution and numerical resolution on computed permeability of consolidated packing using LB and FEM pore-scale simulations. Computers & Fluids 88, 753–763.
| Effects of image resolution and numerical resolution on computed permeability of consolidated packing using LB and FEM pore-scale simulations.Crossref | GoogleScholarGoogle Scholar |
Calonne N, Geindreau C, Flin F, Morin S, Lesaffre B, Rolland Du Roscoat S, Charrier P (2012) 3-D image-based numerical computations of snow permeability: links to specific surface area, density, and microstructural anisotropy. The Cryosphere 6, 939–951.
| 3-D image-based numerical computations of snow permeability: links to specific surface area, density, and microstructural anisotropy.Crossref | GoogleScholarGoogle Scholar |
Chatzigeorgiou G, Picandet V, Khelidj A, Pijaudier-Cabot G (2005) Coupling between progressive damage and permeability of concrete: analysis with a discrete model. International Journal for Numerical and Analytical Methods in Geomechanics 29, 1005–1018.
| Coupling between progressive damage and permeability of concrete: analysis with a discrete model.Crossref | GoogleScholarGoogle Scholar |
Chen G, Weil RR, Hill RL (2014) Effects of compaction and cover crops on soil least limiting water range and air permeability. Soil & Tillage Research 136, 61–69.
| Effects of compaction and cover crops on soil least limiting water range and air permeability.Crossref | GoogleScholarGoogle Scholar |
Cnudde V, Boone MN (2013) High-resolution X-ray computed tomography in geosciences: a review of the current technology and applications. Earth-Science Reviews 123, 1–17.
| High-resolution X-ray computed tomography in geosciences: a review of the current technology and applications.Crossref | GoogleScholarGoogle Scholar |
Cooper D, Turinsky A, Sensen C, Hallgrimsson B (2007) Effect of voxel size on 3D micro-CT analysis of cortical bone porosity. Calcified Tissue International 80, 211–219.
| Effect of voxel size on 3D micro-CT analysis of cortical bone porosity.Crossref | GoogleScholarGoogle Scholar | 17340226PubMed |
Danahy EE, Agaian SS, Panetta KA (2007) Algorithms for the resizing of binary and grayscale images using a logical transform. Proceedings Imaging Processing: Algorithms and Systems 6497, 64970Z
| Algorithms for the resizing of binary and grayscale images using a logical transform.Crossref | GoogleScholarGoogle Scholar |
Doughty DA, Tomutsa L (1996) Multinuclear NMR microscopy of two-phase fluid systems in porous rock. Magnetic Resonance Imaging 14, 869–873.
| Multinuclear NMR microscopy of two-phase fluid systems in porous rock.Crossref | GoogleScholarGoogle Scholar | 8970097PubMed |
Fonseca J, O’Sullivan C, Coop MR, Lee PD (2012) Non-invasive characterization of particle morphology of natural sands. Soil and Foundation 52, 712–722.
| Non-invasive characterization of particle morphology of natural sands.Crossref | GoogleScholarGoogle Scholar |
Fredrich JT, DiGiovanni AA, Noble DR (2006) Predicting macroscopic transport properties using microscopic image data. Journal of Geophysical Research. Solid Earth 111, B03201
| Predicting macroscopic transport properties using microscopic image data.Crossref | GoogleScholarGoogle Scholar |
Grevera GJ, Udupa JK (1996) Shape-based interpolation of multidimensional grey-level images. IEEE Transactions on Medical Imaging 15, 881–892.
| Shape-based interpolation of multidimensional grey-level images.Crossref | GoogleScholarGoogle Scholar | 18215967PubMed |
Guibert R, Nazarova M, Horgue P, Hamon G, Creux P, Debenest G (2015) Computational permeability determination from pore-scale imaging: sample size, mesh and method sensitivities. Transport in Porous Media 107, 641–656.
| Computational permeability determination from pore-scale imaging: sample size, mesh and method sensitivities.Crossref | GoogleScholarGoogle Scholar |
Happel J, Brenner H (2012) ‘Low Reynolds number hydrodynamics: with special applications to particulate media.’ (Springer: Netherlands)
Houston AN, Otten W, Falconer R, Monga O, Baveye PC, Hapca SM (2017) Quantification of the pore size distribution of soils: assessment of existing software using tomographic and synthetic 3D images. Geoderma 299, 73–82.
| Quantification of the pore size distribution of soils: assessment of existing software using tomographic and synthetic 3D images.Crossref | GoogleScholarGoogle Scholar |
Juanes R, Spiteri EJ, Orr FM, Blunt MJ (2006) Impact of relative permeability hysteresis on geological CO2 storage. Water Resources Research 42, W12418.
Kanit T, Forest S, Galliet I, Mounoury V, Jeulin D (2003) Determination of the size of the representative volume element for random composites: statistical and numerical approach. International Journal of Solids and Structures 40, 3647–3679.
| Determination of the size of the representative volume element for random composites: statistical and numerical approach.Crossref | GoogleScholarGoogle Scholar |
Khan F, Enzmann F, Kersten M, Wiegmann A, Steiner K (2012) 3D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver. Journal of Soils and Sediments 12, 86–96.
| 3D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver.Crossref | GoogleScholarGoogle Scholar |
Kumar S, Anderson SH, Udawatta RP, Gantzer CJ (2010) CT-measured macropores as affected by agroforestry and grass buffers for grazed pasture systems. Agroforestry Systems 79, 59–65.
| CT-measured macropores as affected by agroforestry and grass buffers for grazed pasture systems.Crossref | GoogleScholarGoogle Scholar |
Legland D, Kiêu K, Devaux MF (2011) Computation of Minkowski measures on 2D and 3D binary images. Image Analysis & Stereology 26, 83–92.
| Computation of Minkowski measures on 2D and 3D binary images.Crossref | GoogleScholarGoogle Scholar |
Lehmann G, Legland D (2012) Efficient n-dimensional surface estimation using Crofton formula and run-length encoding. The Insight Journal 2, 1–11.
Lehmann TM, Gonner C, Spitzer K (1999) Survey: interpolation methods in medical image processing. IEEE Transactions on Medical Imaging 18, 1049–1075.
| Survey: interpolation methods in medical image processing.Crossref | GoogleScholarGoogle Scholar | 10661324PubMed |
Lehmann P, Wyss P, Flisch A, Lehmann E, Vontobel P, Krafczyk M, Kaestner A, Beckmann F, Gygi A, Flühler H (2006) Tomographical imaging and mathematical description of porous media used for the prediction of fluid distribution. Vadose Zone Journal 5, 80–97.
| Tomographical imaging and mathematical description of porous media used for the prediction of fluid distribution.Crossref | GoogleScholarGoogle Scholar |
Liang H, Song Y, Liu Y, Yang M, Huang X (2010) Study of the permeability characteristics of porous media with methane hydrate by pore network model. Journal of Natural Gas Chemistry 19, 255–260.
| Study of the permeability characteristics of porous media with methane hydrate by pore network model.Crossref | GoogleScholarGoogle Scholar |
Mostaghimi P, Blunt MJ, Bijeljic B (2013) Computations of absolute permeability on micro-CT images. Mathematical Geosciences 45, 103–125.
| Computations of absolute permeability on micro-CT images.Crossref | GoogleScholarGoogle Scholar |
Muñoz-Ortega FJ, San José Martínez F, Caniego Monreal FJ (2015) Volume, surface, connectivity and size distribution of soil pore space in CT images: comparison of samples at different depths from nearby natural and tillage areas. Pure and Applied Geophysics 172, 167–179.
| Volume, surface, connectivity and size distribution of soil pore space in CT images: comparison of samples at different depths from nearby natural and tillage areas.Crossref | GoogleScholarGoogle Scholar |
Nakaya S, Yohmei T, Koike A, Hirayama T, Yoden T, Nishigaki M (2002) Determination of anisotropy of spatial correlation structure in a three-dimensional permeability field accompanied by shallow faults. Water Resources Research 38(8), 35–1–35–14.
Parker JA, Kenyon RV, Troxel DE (1983) Comparison of interpolating methods for image resampling. IEEE Transactions on Medical Imaging 2, 31–39.
| Comparison of interpolating methods for image resampling.Crossref | GoogleScholarGoogle Scholar |
Peng Z, Duwig C, Gaudet JP, Gastelum Strozzi A, Charrier P, Denis H (2015) Visualization and characterization of heterogeneous water flow in double-porosity media by means of X-ray computed tomography. Transport in Porous Media 110, 543–564.
| Visualization and characterization of heterogeneous water flow in double-porosity media by means of X-ray computed tomography.Crossref | GoogleScholarGoogle Scholar |
Reboul N, Vincens E, Cambou B (2008) A statistical analysis of void size distribution in a simulated narrowly graded packing of spheres. Granular Matter 10, 457–468.
| A statistical analysis of void size distribution in a simulated narrowly graded packing of spheres.Crossref | GoogleScholarGoogle Scholar |
Ridler TW, Calvard S (1978) Picture thresholding using an iterative selection method. IEEE Transactions on Systems, Man, and Cybernetics 8, 630–632.
| Picture thresholding using an iterative selection method.Crossref | GoogleScholarGoogle Scholar |
Rozenbaum O, Rolland du Roscoat S (2014) Representative elementary volume assessment of three-dimensional X-ray microtomography images of heterogeneous materials: application to limestones. Physical Review. E 89, 053304
| Representative elementary volume assessment of three-dimensional X-ray microtomography images of heterogeneous materials: application to limestones.Crossref | GoogleScholarGoogle Scholar |
Rutka V (2008) A staggered grid-based explicit jump immersed interface method for two-dimensional Stokes flows. International Journal for Numerical Methods in Fluids 57, 1527–1543.
| A staggered grid-based explicit jump immersed interface method for two-dimensional Stokes flows.Crossref | GoogleScholarGoogle Scholar |
Shah SM, Gray F, Crawshaw JP, Boek ES (2016) Micro-computed tomography porescale study of flow in porous media: effect of voxel resolution. Advances in Water Resources 95, 276–287.
| Micro-computed tomography porescale study of flow in porous media: effect of voxel resolution.Crossref | GoogleScholarGoogle Scholar |
Vogel HJ, Roth K (2001) Quantitative morphology and network representation of soil pore structure. Advances in Water Resources 24, 233–242.
| Quantitative morphology and network representation of soil pore structure.Crossref | GoogleScholarGoogle Scholar |
Vogel HJ, Weller U, Schlüter S (2010) Quantification of soil structure based on Minkowski functions. Computers & Geosciences 36, 1236–1245.
| Quantification of soil structure based on Minkowski functions.Crossref | GoogleScholarGoogle Scholar |
Wiegmann A (2007) ‘Computation of the permeability of porous materials from their microstructure by FFF-Stokes.’ (Bericht des Fraunhofer ITWM: Kaiserslautern, Germany)
Wiegmann A, Bube KP (2000) The explicit-jump immersed interface method: finite difference methods for PDE with piecewise smooth solutions. SIAM Journal on Numerical Analysis 37, 827–862.
| The explicit-jump immersed interface method: finite difference methods for PDE with piecewise smooth solutions.Crossref | GoogleScholarGoogle Scholar |
Yoon H, Dewers TA (2013) Nanopore structures, statistically representative elementary volumes, and transport properties of chalk. Geophysical Research Letters 40, 4294–4298.
| Nanopore structures, statistically representative elementary volumes, and transport properties of chalk.Crossref | GoogleScholarGoogle Scholar |