Optimal X-ray micro-CT image based methods for porosity and permeability quantification in heterogeneous sandstones
Optimal X-ray micro-CT image based methods for porosity and permeability quantification in heterogeneous sandstones
3-D X-ray micro-CT (XCT) is a non-destructive 3-D imaging method, increasingly used for a wide range of applications in Earth Science. An optimal XCT image-processing workflow is derived here for accurate quantification of porosity and absolute permeability of heterogeneous sandstone samples using an assessment of key image acquisition and processing parameters: image resolution, segmentation method, representative elementary volume (REV) size and fluid-simulation method. XCT image-based calculations obtained for heterogeneous sandstones are compared to two homogeneous standards (Berea sandstone and a sphere pack), as well as to the results from physical laboratory measurements. An optimal XCT methodology obtains porosity and permeability results within ±2 per cent and vary by one order of magnitude around the direct physical measurements, respectively, achieved by incorporating the clay fraction and cement matrix (porous, impermeable components) to the pore-phase for porosity calculations and into the solid-phase for permeability calculations. Two stokes-flow finite element modelling (FEM) simulation methods, using a voxelized grid (Avizo) and tetrahedral mesh (Comsol) produce comparable results, and similarly show that a lower resolution scan (∼5 μm) is unable to resolve the smallest intergranular pores, causing an underestimation of porosity by ∼3.5 per cent. Downsampling the image-resolution post-segmentation (numerical coarsening) and pore network modelling both allow achieving of a REV size, whilst significantly reducing fluid simulation memory requirements. For the heterogeneous sandstones, REV size for permeability (≥1 mm3) is larger than for porosity (≥0.5 mm3) due to tortuosity of the fluid paths. This highlights that porosity should not be used as a reference REV for permeability calculations. The findings suggest that distinct image processing workflows for porosity and permeability would significantly enhance the accurate quantification of the two properties from XCT.
Core, Microstructure, Numerical modelling, X-ray micro-CT, image processing, permeability and porosity
1210-1229
Callow, Ben
15166203-d3e6-4b28-8369-e99e1bd00240
Falcon-Suarez, Ismael
f5cdbc61-326b-424d-a90f-593a8698a4d2
Marin-Moreno, Hector
e466cafd-bd5c-47a1-8522-e6938e7086a4
Bull, Jonathan M.
974037fd-544b-458f-98cc-ce8eca89e3c8
Ahmed, Sharif
37570e92-ba6b-4e03-9144-c70fa7722c51
Callow, Ben
15166203-d3e6-4b28-8369-e99e1bd00240
Falcon-Suarez, Ismael
f5cdbc61-326b-424d-a90f-593a8698a4d2
Marin-Moreno, Hector
e466cafd-bd5c-47a1-8522-e6938e7086a4
Bull, Jonathan M.
974037fd-544b-458f-98cc-ce8eca89e3c8
Ahmed, Sharif
37570e92-ba6b-4e03-9144-c70fa7722c51
Callow, Ben, Falcon-Suarez, Ismael, Marin-Moreno, Hector, Bull, Jonathan M. and Ahmed, Sharif
(2020)
Optimal X-ray micro-CT image based methods for porosity and permeability quantification in heterogeneous sandstones.
Geophysical Journal International, 223 (2), .
(doi:10.1093/gji/ggaa321).
Abstract
3-D X-ray micro-CT (XCT) is a non-destructive 3-D imaging method, increasingly used for a wide range of applications in Earth Science. An optimal XCT image-processing workflow is derived here for accurate quantification of porosity and absolute permeability of heterogeneous sandstone samples using an assessment of key image acquisition and processing parameters: image resolution, segmentation method, representative elementary volume (REV) size and fluid-simulation method. XCT image-based calculations obtained for heterogeneous sandstones are compared to two homogeneous standards (Berea sandstone and a sphere pack), as well as to the results from physical laboratory measurements. An optimal XCT methodology obtains porosity and permeability results within ±2 per cent and vary by one order of magnitude around the direct physical measurements, respectively, achieved by incorporating the clay fraction and cement matrix (porous, impermeable components) to the pore-phase for porosity calculations and into the solid-phase for permeability calculations. Two stokes-flow finite element modelling (FEM) simulation methods, using a voxelized grid (Avizo) and tetrahedral mesh (Comsol) produce comparable results, and similarly show that a lower resolution scan (∼5 μm) is unable to resolve the smallest intergranular pores, causing an underestimation of porosity by ∼3.5 per cent. Downsampling the image-resolution post-segmentation (numerical coarsening) and pore network modelling both allow achieving of a REV size, whilst significantly reducing fluid simulation memory requirements. For the heterogeneous sandstones, REV size for permeability (≥1 mm3) is larger than for porosity (≥0.5 mm3) due to tortuosity of the fluid paths. This highlights that porosity should not be used as a reference REV for permeability calculations. The findings suggest that distinct image processing workflows for porosity and permeability would significantly enhance the accurate quantification of the two properties from XCT.
Text
Callow_et_al_2020_GJI-S-19-1155.R1
- Accepted Manuscript
More information
Accepted/In Press date: 25 June 2020
e-pub ahead of print date: 27 June 2020
Additional Information:
Funding Information: we acknowledge funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 654462 - STEMM-CCS. We acknowledge Diamond Light Source for time on Beamline I13-2 under Proposal MT18758. We give our thanks to Hans Deyhle, Laurence North, Christina Reinhard, Shashi Marathe, Andrew Bodey and the Diamond Support Scientists at I13-2 for their support during the experiment. We are also very grateful for the advice and guidance with image reconstruction provided by Kaz Wanelik. We acknowledge the use of the IRIDIS High Performance Computing Facility, and associated support services at the University of Southampton, in the completion of this work. We would like to thank our colleagues at the ì-VIS X-Ray Imaging Centre, University of Southampton for their helpful advice and support. We would also like to thank Andrew Hurst and Antonio Grippa (University of Aberdeen, Sand Injection Research Group) for their assistance and advice with fieldwork and rock sampling in California. We would like to thank the FEI Visualization Sciences Group for providing the use of the Avizo 9.3.0 software for XCT image processing, Simpleware for the use of their ScanIP meshing software and Comsol for use of theirMultiphysics software. Finally, we would like to sincerely thank Samuel Jackson and two anonymous reviewers for their constructive comments, which greatly improved the manuscript. Author contribution statement: B.C., I.F.S.: Conceptualization, Formal Analysis, Validation. B.C.: Data curation, Methodology, Software, Visualization. B.C., J.M.B., S.A.: Funding acquisition. B.C., I.F.S., H.M.M., S.A.: Investigation. B.C., I.F.S., H.M.M., J.M.B., S.A: Writing. The data set associated with this work is available online at https://doi.org/10.5258/SOTON/D1441. Here we provide the raw and segmented XCT tomographic image data, and the associated image processing files.
Keywords:
Core, Microstructure, Numerical modelling, X-ray micro-CT, image processing, permeability and porosity
Identifiers
Local EPrints ID: 441855
URI: http://eprints.soton.ac.uk/id/eprint/441855
ISSN: 0956-540X
PURE UUID: 632b5eaa-372f-4141-945c-a13fcf5a3b6e
Catalogue record
Date deposited: 30 Jun 2020 16:31
Last modified: 03 Jul 2025 02:37
Export record
Altmetrics
Contributors
Author:
Ben Callow
Author:
Ismael Falcon-Suarez
Author:
Hector Marin-Moreno
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics