skip to main content

Title: Multiscale Digital Porous Rock Reconstruction Using Template Matching

Rocks are heterogeneous multiscale porous media: two rock samples with identical bulk properties can vary widely in microstructure. The advent of digital rock technology and modern 3‐D printing provides new opportunities to replicate rocks. However, the inherent trade‐off between imaging resolution and sample size limits the scales over which microstructure and macrostructure can be identified and related to each other. Here, we develop a multiscale digital rock construction strategy by combining X‐ray computed microtomography and focused‐ion beam (FIB)‐scanning electron microscope (SEM) images, and we apply the technique to a tight sandstone. The computed tomography (CT) scanning images characterize macroscale pore structures, while the FIB‐SEM images capture microscale pore textures. The FIB‐SEM images are then coupled to CT images via a template‐matching algorithm and superposition. Bulk properties, including porosity and pore and throat size distribution, can be recovered with this approach. Permeability prediction with a pore network model for the largest connected pore network are 3 orders and 1 order of magnitude greater than the bulk rock measured value using the CT‐only and the SEM‐CT coupled images, respectively.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Page Range / eLocation ID:
p. 6911-6922
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Three‐dimensional (3D) X‐ray computed tomography (X‐ray CT) imaging has emerged as a nondestructive means of microstructural characterization. However, obtaining and processing high‐quality and high‐resolution images is time‐consuming and often requires high‐performance computing, particularly with a high number of projections. This work evaluates the effect of 3D X‐ray CT imaging parameters on pore connectivity and surface area quantification in sandstone samples of varied composition. Samples from Bentheimer and Torrey Buff formations are imaged via 3D X‐ray CT imaging at resolutions ranging from 1.25 to 15 μm, bin sizes 1, 2, and 4, and number of projections from 400 to 4,500. Collected images are processed and analyzed using ImageJ and MATLAB to discern petrophysical properties and the results are compared with each other and Mercury Intrusion Capillary Pressure (MICP) results. Overall, little variation in bulk porosity with changing scanning parameters is observed. However, for low resolution and projection numbers, connected porosity is lower compared to bulk porosity due to a failure to capture microfeatures. Overall, mineral surface area is observed to decrease with increasing bin size, voxel size, and projection numbers, except an observed increase with projection numbers for Torrey Buff. The Torrey Buff samples contains comparatively more clays and even the highest resolution (1.25 μm) fails to separate the micrograins, which is reflected in the pore size distribution. Identifying these variations are helpful as discrepancies in imaged pore connectivity and surface area can largely impact assessments of fluid flow and transport in reactive transport simulations informed by this data.

    more » « less
  2. Focused ion beam (FIB) – scanning electron microscopy (SEM) allowed the characterization of the microstructure of two solid oxide fuel cells prepared at different sintering temperatures. 3D volume reconstruction showed that a relatively low sintering temperature significantly and positively affected distribution, volume and particle size of yttria-stabilized zirconia, nickel, and pore phases inside the anode, as well as the extent of the important triple-phase boundary interface. The poor performance of the T1 sample sintered at a higher temperature is explained by the poorly connected pore network and very low-density triple-phase boundary. The pore space inside the T1 anode was unable to ensure continuous hydrogen flow from the inlet to the outlet and thus exhibited very low gas permeability. In contrast, the T2 sample sintered at a lower temperature had approximately equal amounts of YSZ and nickel and larger pores, which allowed formation of significantly more TPB electrochemical reaction sites. The higher power density of the T2 cell was also the result of its robust pore network capable of transporting hydrogen throughout the anode. The methodology used in this paper eliminates the need for employing hypothetical structures and provides accurate estimates of the investigated parameters by evaluating microstructures that were successfully reconstructed using high-resolution microscopy techniques. 
    more » « less
  3. Abstract

    A typical ground investigation for characterizing geotechnical properties of soil requires sampling soils to test in a laboratory. Laboratory X-ray computed tomography (CT) has been used to non-destructively observe soils and characterize their properties using image processing, numerical analysis, or three-dimensional (3D) printing techniques based on scanned images; however, if it becomes possible to scan the soils in the ground, it may enable the characterization without sampling them. In this study, an in-situ X-ray CT scanning system comprising a drilling machine with an integrated CT scanner was developed. A model test was conducted on gravel soil to verify if the equipment can drill and scan the soil underground. Moreover, image processing was performed on acquired 3D CT images to verify the image quality; the particle morphology (particle size and shape characteristics) was compared with the results obtained for projected particles captured in a two-dimensional (2D) manner by a digital camera. The equipment successfully drilled to a target depth of 800 mm, and the soil was scanned at depths of 700, 750, and 800 mm. Image processing results showed a reasonable agreement between the 3D and 2D particle morphology images, and confirmed the feasibility of the in-situ X-ray CT scanning system.

    more » « less

    The chemo-mechanical loading of rocks causes the dissolution and precipitation of multiple phases in the rock. This dissolution and precipitation of load-bearing mineral phases lead to the stress redistribution in neighboring phases, which in turn results in deformational changes of the sample composite. The aim of this study is to investigate the link between microstructural evolution and creep behavior of shale rocks subjected to chemo-mechanical loading through modeling time-dependent deformation induced by the dissolution-precipitation process. The model couples the microstructural evolution of the shale rocks with the stress/strain fields inside the material as a function of time. The modeling effort is supplemented with an experimental study where shale rocks were exposed to CO2-rich brine under high temperature and pressure conditions. 3D snapshots of the sample microstructure were generated using segmented micro-CT images of the shale sample. The time-evolving microstructures were then integrated with the Finite element-based mechanical model to simulate the creep induced by dissolution and precipitation processes independent of the intrinsic viscoelasticity/viscoplasticity of the mineral phases. After computation of the time-dependent viscoelastic properties of the shale composite, the combined microstructure model and finite element model were utilized to predict the time-dependent stress and strain fields in different zones of reacted shale.


    Determination of viscous behavior of shale rocks is key in wide range of applications such as stability of reservoirs, stability of geo-structures subjected to environmental forcing, underground storage of hazardous materials and hydraulic fracturing. Short-term creep strains in hydraulic fracturing can change stress fields and in turn can impact the hydraulic fracturing procedures(H. Sone & Zoback, 2010; Hiroki Sone & Zoback, 2013). While long-term creep strains can hamper the reservoir performance due to the reduction in permeability of the reservoir by closing of fractures and fissures(Du, Hu, Meegoda, & Zhang, 2018; Rybacki, Meier, & Dresen, 2016; Sharma, Prakash, & Abedi, 2019; Hiroki Sone & Zoback, 2014). Owing to these significance of creep strain, it is important to understand the viscoelastic/viscoplastic behavior of shales.

    more » « less
  5. Abstract

    Accurate prediction of physical alterations in carbonate reservoirs under dissolution is critical for development of subsurface energy technologies. The impact of mineral dissolution on flow characteristics depends on the connectivity and tortuosity of the pore network. Persistent homology is a tool from algebraic topology that describes the size and connectivity of topological features. When applied to 3D X‐ray computed tomography (XCT) imagery of rock cores, it provides a novel metric of pore network heterogeneity. Prior works have demonstrated the efficacy of persistent homology in predicting flow properties in numerical simulations of flow through porous media. Its ability to combine size, spatial distribution, and connectivity information make it a promising tool for understanding reactive transport in complex pore networks, yet limited work has been done to apply persistence analysis to experimental studies on natural rocks. In this study, three limestone cores were imaged by XCT before and after acid‐driven dissolution flow‐through experiments. Each XCT scan was analyzed using persistent homology. In all three rocks, permeability increase was driven by the growth of large, connected pore bodies. The two most homogenous samples saw an increased effect nearer to the flow inlet, suggesting emerging preferential flow paths as the reaction front progresses. The most heterogeneous sample showed an increase in along‐core homogeneity during reaction. Variability of persistence showed moderate positive correlation with pore body size increase. Persistence heterogeneity analysis could be used to anticipate where greatest pore size evolution may occur in a reservoir targeted for subsurface development, improving confidence in project viability.

    more » « less