Critical Reynolds number for nonlinear flow through rough-walled fractures: The role of shear processes

This paper experimentally investigates the role of shear processes on the variation of critical Reynolds number and nonlinear flow through rough-walled rock fractures. A quantitative criterion was developed to quantify the onset of nonlinear flow by comprehensive combination of Forchheimer's law and Reynolds number. At each shear displacement, several high-precision water flow tests were carried out with different hydraulic gradients then the critical Reynolds number was determined based on the developed criterion. The results show that (i) the Forchheimer's law was fitted very well to experimental results of nonlinear fluid flow through rough-walled fractures, (ii) the coefficients of viscous and inertial pressure drops experience 4 and 7 orders of magnitude reduction during shear displacement, respectively, and (iii) the critical Reynolds number varies from 0.001 to 25 and experiences 4 orders of magnitude enlargement by increasing shear displacement from 0 to 20 mm. These findings may prove useful in proper understanding of fluid flow through rock fractures, or inclusions in computational studies of large-scale nonlinear flow in fractured rocks

Quantifying fracture geometry with X-ray tomography: Technique of Iterative Local Thresholding (TILT) for 3D image segmentation

This paper presents a new method—the Technique of Iterative Local Thresholding (TILT)—for processing 3D X-ray computed tomography (xCT) images for visualization and quantification of rock fractures. The TILT method includes the following advancements. First, custom masks are generated by a fracture-dilation procedure, which significantly amplifies the fracture signal on the intensity histogram used for local thresholding. Second, TILT is particularly well suited for fracture characterization in granular rocks because the multi-scale Hessian fracture (MHF) filter has been incorporated to distinguish fractures from pores in the rock matrix. Third, TILT wraps the thresholding and fracture isolation steps in an optimized iterative routine for binary segmentation, minimizing human intervention and enabling automated processing of large 3D datasets. As an illustrative example, we applied TILT to 3D xCT images of reacted and unreacted fractured limestone cores. Other segmentation methods were also applied to provide insights regarding variability in image processing. The results show that TILT significantly enhanced separability of grayscale intensities, outperformed the other methods in automation, and was successful in isolating fractures from the porous rock matrix. Because the other methods are more likely to misclassify fracture edges as void and/or have limited capacity in distinguishing fractures from pores, those methods estimated larger fracture volumes (up to 80 %), surface areas (up to 60 %), and roughness (up to a factor of 2). These differences in fracture geometry would lead to significant disparities in hydraulic permeability predictions, as determined by 2D flow simulations