The direct flow simulation using high‐resolution micro‐computed tomographic (μ‐CT) images of porous rock can be used to help understand the flow characteristics at the pore‐scale and to estimate fluid properties; however, segmentation of pore space in grayscale 3‐D μ‐CT images, a necessary step in this process, is challenging because of issues related to the image resolution and pore‐filling matrix in the gray‐level images. We present a novel process for determining the voxel porosity and permeability of the gray‐level regime and evaluating the bulk permeability using the Brinkman force lattice Boltzmann method. In this study, μ‐CT images of Berea sandstone are acquired with two spatial resolutions. After the pore size distribution curve is experimentally obtained, the “apparent pore” voxels and “gray pore” voxels are determined based on the designated gray‐level values in the gray‐level (CT number) histogram, the cumulative and fractional pore volumes, and the linear relationship between CT number and individual voxel porosity. The results show that the boundary between apparent and gray pores in terms of size determines the volumetric fraction and specific surface area. Additionally, the permeability computed by considering the gray pore regime based on the proposed method is more similar to the experimentally measured value than the results of other segmentation methods because the gray pore domain is fully incorporated into the flow domain while preserving the pore connectivity. The proposed sequence of pore segmentation presents a method for handling gray‐level pore space without compromising pore connectivity.
We propose a 3D-shape preferred orientation (SPO) measurement method of rigid grains using synchrotron micro-computational tomography (μ-CT). The method includes oriented sampling, 3D μ-CT imaging, image filtering, ellipsoid fitting, and SPO measurement. After CT imaging, all processes are computerized, and the directions of thousands of rigid grains in 3D-space can be automatically measured. This method is optimized for estimating the orientation of the silt-sized rigid grains in fault gouge, which indicates P-shear direction in a fault system. This allows us to successfully deduce fault motion sense and quantify fault movement. Because this method requires a small amount of sample, it can be applied as an alternative to study fault systems, where the shear sense indicators are not distinct in the outcrop and the fault gouge is poorly developed. We applied the newly developed 3D-SPO method for a fault system in the Yangsan fault, one of the major faults in the southeastern Korean Peninsula, and observed the P-shear direction successfully.
The hydraulic permeability is a key parameter for simulating the flow-related phenomenon so that its accurate estimation is crucial in both experimental and numerical simulation studies. 3D pore structure can be readily taken by X-ray computed tomography (CT) and it often serves as a flow domain for pore-scale simulation. However, one encounters the challenges in segmenting the authentic pore structure owing to the finite size of image resolution and segmentation methods. Therefore, the loss of structural information in pore space seems unavoidable to result in the unreliable estimation of permeability. In this study, we propose a novel framework to overcome these limitations by using a flexible ternary segmentation scheme. Given the pore size distribution curve and porosity, three phases of pore, solid, and gray regions are segmented by considering the partial volume effect which holds the composition information of unresolved objects. The resolved objects such as solid and pore phases are taken to equivalently solve Stokes equation while the fluid flow through unresolved objects is simultaneously solved by Stokes-Brinkmann equation. The proposed numerical scheme to obtain the permeability is applied to Indiana limestone and Navajo sandstone. The results show that the computed hydraulic permeability is similar to the experimentally obtained value without being affected by image resolution. This approach has advantages of achieving consistent permeability values, less influenced by segmentation methods.
Fault gouges play an important role in the shear deformation of fault zones, by causing weakness and frictional instability in structures. Previous studies have investigated the evolution of shear deformation of fault zones by observing experiments using remolded and synthetic gouge specimens at a micro-scale. However, how the spatial configuration of the rock constituents accounts for the 3D anisotropy of intact structures of fault gouges, particularly at the core-scale, is not well understood. We obtained 3D μ-CT images of directionally cored gouge specimens and performed statistical analysis to quantify the major orientation of the internal structures. Direct shear tests were conducted to investigate the relationship between the distribution of the internal structures and geomechanical behavior. The results show that the undisturbed fault gouge has a clear anisotropy parallel to the fault plane even at the core-scale. Moreover, the direct shear test results show that the frictional resistance of a fault gouge has anisotropy related to the fault plane. The simple, yet robust method proposed in this study confirms that the core-scale structural anisotropy is correlated to the anisotropic shear resistance.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.