An official website of the United States government
Pore-scale modeling is essential in understanding and predicting flow and transport properties of rocks. Generally, pore-scale modeling is dependent on imaging technologies such as Micro Computed Tomography (micro-CT), which provides visual confirmation into the pore microstructures of rocks at a representative scale. However, this technique is limited in the ability to provide high resolution images showing the pore-throats connecting pore bodies. Pore scale simulations of flow and transport properties of rocks are generally done on a single 3D pore microstructure image. As such, the simulated properties are only representative of the simulated pore-scale rock volume. These are the technological and computational limitations which we address here by using a stochastic pore-scale simulation approach. This approach consists of constructing hundreds of 3D pore microstructures of the same pore size distribution and overall porosity but different pore connectivity. The construction of the 3D pore microstructures incorporates the use of Mercury Injection Capillary Pressure (MICP) data to account for pore throat size distribution, and micro-CT images to account for pore body size distribution. The approach requires a small micro-CT image volume (7–19 mm3) to reveal key pore microstructural features that control flow and transport properties of highly heterogeneous rocks at the core-scale. Four carbonate rock samples were used to test the proposed approach. Permeability calculations from the introduced approach were validated by comparing laboratory measured permeability of rock cores and permeability estimated using five well-known core-scale empirical model equations. The results show that accounting for the stochastic connectivity of pores results in a probabilistic distribution of flow properties which can be used to upscale pore-scale simulated flow properties to the core-scale. The use of the introduced stochastic pore-scale simulation approach is more beneficial when there is a higher degree of heterogeneity in pore size distribution. This is shown to be the case with permeability and hydraulic tortuosity which are key controls of flow and transport processes in rocks.
more »
« less
Multiscale Digital Porous Rock Reconstruction Using Template Matching
Abstract 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):
- 1724469
- PAR ID:
- 10455753
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 55
- Issue:
- 8
- ISSN:
- 0043-1397
- Page Range / eLocation ID:
- p. 6911-6922
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Effect of X‐Ray Computed Tomography Imaging Parameters on Quantification of Petrophysical PropertiesAbstract 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
-
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
-
AbstractElectrokinetic in-situ recovery is an alternative to conventional mining, relying on the application of an electric potential to enhance the subsurface flow of ions. Understanding the pore-scale flow and ion transport under electric potential is essential for petrophysical properties estimation and flow behavior characterization. The governing physics of electrokinetic transport is electromigration and electroosmotic flow, which depend on the electric potential gradient, mineral occurrence, domain morphology (tortuosity and porosity, grain size and distribution, etc.), and electrolyte properties (local pH distribution and lixiviant type and concentration, etc.). Herein, mineral occurrence and its associated zeta potential are investigated for EK transport. The new Ek model which is designed to solve the EK flow in complex porous media in a highly parallelizable manner includes three coupled equations: (1) Poisson equation, (2) Nernst–Planck equation, and (3) Navier–Stokes equation. These equations were solved using the lattice Boltzmann method within X-ray computed microtomography images. The proposed model is validated against COMSOL multiphysics in a two-dimensional microchannel in terms of fluid flow behavior when the electrical double layer is both resolvable and unresolvable. A more complex chalcopyrite-silica system is then obtained by micro-CT scanning to evaluate the model performance. The effects of mineral occurrence, zeta potential, and electric potential on the three-dimensional chalcopyrite-silica system were evaluated. Although the positive zeta potential of chalcopyrite can induce a flow of ferric ion counter to the direction of electromigration, the net effect is dependent on the occurrence of chalcopyrite. However, the ion flux induced by electromigration was the dominant transport mechanism, whereas advection induced by electroosmosis made a lower contribution. Overall, a pore-scale EK model is proposed for direct simulation on pore-scale images. The proposed model can be coupled with other geochemical models for full physicochemical transport simulations. Meanwhile, electrokinetic transport shows promise as a human-controllable technique because the electromigration of ions and the applied electric potential can be easily controlled externally. Graphical abstractmore » « less
-
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
