Digital rock core pore analysis combines high-resolution X-ray CT imaging, advanced image processing, and pore-scale numerical simulation to quantitatively characterize the three-dimensional pore structure of reservoir and geotechnical rocks. This article explains the technical workflow, sample requirements, measurable parameters, applicable standards, and typical applications of digital core pore analysis. It serves as a practical reference for petroleum engineers, geoscientists, and laboratory professionals seeking non-destructive, voxel-level insight into rock porosity and permeability.
Technical Overview
Digital rock core pore analysis, also known as digital rock physics (DRP) pore characterization, is a laboratory technique that converts a physical rock sample into a high-resolution three-dimensional digital volume and then quantifies its pore system through image analysis and numerical simulation. A cylindrical core plug or rock fragment is scanned with X-ray computed tomography (micro-CT or nano-CT) or imaged with focused ion beam scanning electron microscopy (FIB-SEM). The resulting gray-scale tomogram is reconstructed, filtered, and segmented into solid grains, pore space, and mineral phases.
From the segmented 3D volume, a pore network model is extracted that represents the rock as a lattice of pore bodies connected by throats. Quantitative descriptors such as porosity, pore size distribution, throat radius distribution, coordination number, and tortuosity are computed directly from the digital volume. Fluid flow properties, including absolute permeability and relative permeability, can be estimated using Lattice Boltzmann methods, direct Navier-Stokes solvers, or pore network simulation, providing a non-destructive complement to conventional core flooding and mercury intrusion testing.
The complete workflow typically includes sample preparation, imaging, image reconstruction and denoising, phase segmentation (thresholding, watershed, or machine-learning based), pore network extraction, petrophysical property calculation, and validation against physical measurements performed by accredited third-party laboratories following ISO, ASTM, or GB standards.
Testing Samples
Conventional Sandstone and Conglomerate Core Plugs
Well-cemented sandstones are the most common samples for micro-CT digital pore analysis. Mini-core plugs with diameters of 2 to 10 mm are drilled parallel or perpendicular to bedding, cleaned of hydrocarbons and salts, and dried before scanning. At voxel resolutions of 1 to 10 micrometers, intergranular pores, pore-filling cements, and micro-fractures can all be resolved with high fidelity.
Carbonate Rocks
Limestones and dolomites often contain vugs, molds, and microporosity spanning several orders of magnitude in scale. Depending on heterogeneity, carbonate samples may require multi-scale imaging: whole-core medical or industrial CT for vug distribution, micro-CT for mesopores, and FIB-SEM or argon-ion milled SEM imaging for micrite-hosted micropores. Larger or irregularly shaped fragments can be scanned in multiple stages and registered into a composite volume.
Shale and Tight Reservoir Rocks
Nanoscale organic-matter pores and clay-hosted micropores in shale demand resolutions from a few nanometers to a few hundred nanometers. FIB-SEM serial sectioning and broad ion beam (BIB-SEM) techniques are used to build representative 3D volumes, while nano-CT captures larger-scale microfractures and pyrite or organic nodule distributions. Sample preparation includes low-damage polishing with ion milling to expose undamaged pore geometry.
Unconsolidated Sediments and Drill Cuttings
Weakly consolidated sands and cuttings can be analyzed by carefully mounting them in low-attenuation sleeves or by scanning frozen or resin-stabilized material. Drill cuttings, when full-diameter cores are unavailable, provide small but useful volumes for porosity typing, grain characterization, and screening-level digital permeability estimation.
Testing Items
Digital core pore analysis delivers a comprehensive suite of quantitative pore-structure and petrophysical parameters, including but not limited to:
- Total, effective (connected), and isolated porosity, expressed as voxel-based volume fractions
- Pore body size distribution and pore throat radius distribution
- Pore-throat aspect ratio and shape factors
- Coordination number and pore network connectivity statistics
- Absolute permeability (single-phase) computed by Lattice Boltzmann or direct Stokes flow simulation
- Relative permeability and capillary pressure curves derived from two-phase pore network or direct simulations
- Tortuosity of connected flow paths and formation factor estimation
- Specific surface area of pore-solid interfaces
- Fractal dimension of the pore space
- Mineral phase volume fractions, including clay, carbonate cement, and organic matter segmentation
- Micro-fracture aperture, orientation, and density
- Fluid saturation distribution and residual fluid configuration after drainage or imbibition simulation
Testing Methods
The core of the method is non-destructive X-ray micro-CT imaging. The rock sample is rotated 360 degrees while a microfocus X-ray source illuminates it; transmitted projections are reconstructed via filtered back-projection or iterative algorithms into a 3D gray-scale volume. Laboratory micro-CT systems typically achieve spatial resolutions from approximately 0.5 to 50 micrometers depending on sample size, while synchrotron CT and nano-CT extend resolution below 400 nanometers. For nanoscale pore systems, FIB-SEM volumes are constructed by alternating ion-beam milling and SEM imaging in a fully automated serial-sectioning routine.
Image post-processing applies edge-preserving filters (for example, non-local means or anisotropic diffusion) followed by segmentation. Machine-learning classifiers trained on manually labeled slices are increasingly used to improve phase separation in complex lithologies. Pore network models are then extracted using maximal-ball or medial-axis algorithms, converting the continuous pore space into a graph of pore bodies and throats from which geometrical and topological properties are computed.
Numerical simulation methods include single-phase Lattice Boltzmann and finite-volume Stokes solvers for permeability, two-phase pore network invasion algorithms for capillary pressure and relative permeability, and percolation analysis for connectivity assessment. Digital results should be benchmarked against independent physical measurements, and commonly referenced standards and practices include:
- ISO 15901 series for pore size distribution and porosity determination by mercury intrusion and gas adsorption, used for cross-validation
- ASTM D4404 for pore volume and pore volume distribution of rock by mercury intrusion porosimetry
- ASTM D4543 for preparation of rock core specimens and precision statements for physical benchmarking
- GB/T 29172 covering core analysis practices from which sampling and conditioning procedures are adapted
- ISO 15901-2 gas adsorption methods for nanopore volume and surface area validation in shale and tight rocks
Applications
Digital pore analysis supports a wide range of energy, environmental, and geotechnical workflows:
- Reservoir characterization: rock typing, flow-unit classification, and calibration of wireline log porosity-permeability relationships
- Unconventional reservoirs: nanoscale porosity quantification in shale gas, tight oil, and coalbed methane formations
- Enhanced oil recovery: screening of wettability alteration, gas injection, and chemical flooding through pore-scale two-phase simulation
- CO2 geological storage: caprock seal integrity, residual trapping, and injectivity assessment at the pore scale
- Geothermal reservoir evaluation: fracture-matrix interaction and thermo-hydraulic property upscaling
- Nuclear waste and geotechnical engineering: long-term transport properties of host rocks and clay barriers
- Formation damage diagnosis: identification of fines migration, clay swelling, and scale precipitation pathways
- Digital twin development: generation of training datasets for machine-learning prediction of rock properties from logs or images
Because the method is non-destructive, the same sample can be imaged repeatedly, for example before and after reactive fluid exposure, enabling time-lapse monitoring of dissolution, precipitation, and mechanical compaction within the pore network.
FAQ
What sample size and resolution are required?
Sample diameter and voxel resolution are linked: smaller samples yield finer voxels. A 5 mm mini-core typically achieves 2 to 5 micrometer resolution, suitable for sandstones, whereas shale nanopores require FIB-SEM volumes of roughly 10 to 30 micrometers edge length at 5 to 20 nanometer resolution. The laboratory selects the combination that best resolves the target pore size range while maintaining a representative elementary volume.
How representative is a small digital volume?
Representativeness is assessed through the representative elementary volume (REV) analysis: porosity and connectivity are computed on progressively larger sub-volumes until values stabilize. For heterogeneous carbonates, multi-scale imaging and domain decomposition are used so that both vuggy macroporosity and interparticle microporosity are captured.
How accurate is digital permeability compared with physical measurement?
For clean sandstones imaged at adequate resolution, simulated permeability commonly agrees with steady-state core flooding within 10 to 30 percent. Discrepancies in carbonates and shales arise from unresolved microporosity, so results are routinely reconciled with mercury intrusion (ASTM D4404, ISO 15901-1), gas permeability, and gas adsorption (ISO 15901-2) data from accredited third-party laboratories.
Can digital analysis replace conventional core analysis?
Digital pore analysis is best regarded as a complementary, high-information-density extension of conventional core analysis. It provides 3D pore geometry and simulation capability that physical tests cannot, while conventional plug measurements remain essential for validation and for properties governed by larger-scale heterogeneity.
What is a typical project timeline?
A standard workflow covering sample preparation, micro-CT scanning, segmentation, pore network extraction, and permeability simulation generally takes one to three weeks per sample group. Projects involving FIB-SEM volumes, multi-scale registration, or history-matched two-phase simulation may require additional time for image processing and computational runs.