Importance of shear dilation to two-phase flow in naturally fractured geological media: A numerical study using zero-thickness interface elements

In this paper, a hydro-mechanical coupling model is proposed to investigate the geomechanical effects on two-phase flow in fractured rocks. The discrete fracture model (DFM) and modified Barton-Bandis's model are used to model two-phase fluid flow and mechanical behavior of fractures, respectively. The coupled problem is solved by a mixed finite volume/finite element method, together with the fixed-stress split algorithm. Our model is able to capture the complex coupled behavior induced by the interaction between pore pressure and in-situ stress loadings. Then the proposed model is applied to investigate the two-phase flow behaviors in the fracture networks with different fracture configurations under various in-situ stress conditions. The results illustrate significant influences of shear dilation on the hydraulic properties of fractured rock. With the increase of differential stress, fracture aperture gradually increases due to shear dilation, which accordingly enhances the equivalent permeability of the fractured medium. The degree of shear dilation is highly dependent on the fracture orientation and length distribution. Consequently, channelized flow is formed through the dilated fractures, which results in early water breakthrough and reduces the water sweep efficiency. This study highlights the necessity of considering hydro-mechanical coupling effect for accurate prediction of saturation distributions in fractured rocks.