The capillary pressure vs. saturation curve for a fractured rock mass: fracture and matrix contributions

Alejandro Cardona; Qi Liu; J Carlos Santamarina

doi:10.1038/s41598-023-38737-y

The capillary pressure vs. saturation curve for a fractured rock mass: fracture and matrix contributions

Sci Rep. 2023 Jul 25;13(1):12044. doi: 10.1038/s41598-023-38737-y.

Authors

Alejandro Cardona^{1

2}, Qi Liu³, J Carlos Santamarina^{2

4}

Affiliations

¹ Institute for Geophysics, Jackson School of Geosciences, The University of Texas at Austin, 10601 Exploration Way, Austin, TX, 78758, USA.
² King Abdullah University of Science and Technology KAUST, Thuwal, 23955-6900, Saudi Arabia.
³ King Abdullah University of Science and Technology KAUST, Thuwal, 23955-6900, Saudi Arabia. qi.liu@kaust.edu.sa.
⁴ School of Civil and Environmental Engineering, Georgia Institute of Technology, 790 Atlantic Drive, N.W., Atlanta, GA, 30332, USA.

Abstract

The fractal topography of fracture surfaces challenges the upscaling of laboratory test results to the field scale, therefore the study of rock masses often requires numerical experimentation. We generate digital fracture analogues and model invasion percolation to investigate the capillarity-saturation P_c-S_w fracture response to changes in boundary conditions. Results show that aperture is Gaussian-distributed and the coefficient of variation is scale-independent. The aperture contraction during normal stress increments causes higher capillary pressures and steeper P_c-S_w curves, while shear displacement results in invasion anisotropy. The three-parameter van Genutchen model adequately fits the fracture capillary response in all cases; the capillary entry value decreases with fracture size, yet the fracture P_c-S_w curve normalized by the entry value is size-independent. Finally, we combine the fracture and matrix response to infer the rock mass response. Fracture spacing, aperture statistics and matrix porosity determine the rock mass capillarity-saturation P_c-S_w curve. Fractures without gouge control the entry pressure whereas the matrix regulates the residual saturation at high capillary pressure P_c.