Simulating two-phase flow and geomechanical deformation in fractured karst reservoirs based on a coupled hydro-mechanical model

Abstract Two-phase flow in fractured and karstified porous media subject to coupled hydro-mechanical conditions is an important issue for oil recovery in carbonate reservoirs. However, due to the co-existence of porous media flow, fracture flow and free flow, as well as their couplings with geomechanical deformation, modeling the behavior of fractured karst systems remains challenging. In this work, a novel coupled hydro-mechanical model for simulating the complex behavior of fractured and karstified porous media is developed. Two-phase Darcy's equation is used to describe fluid flow in both matrix and fractures, while the free flow in cavities is considered based on an assumption of phase instantaneous gravity segregation. A modified Barton-Bandis's constitutive model is used to mimic the nonlinear fracture deformation. The cavity deformation is solved based on the fluid pressure on the cavity boundaries. A mixed finite volume-finite element method and a fixed-stress iterative splitting method are adopted to numerically solve the coupled system of equations. The model is then applied to a series of 2D and 3D problems to unravel the impacts of fractures and cavities on two-phase flow and geomechanical deformation in fractured karst systems. The results indicate that cavities hinder water breakthrough due to storage effects, while water may quickly migrate through highly conductive fractures. Cavities tend to dominate the flow and mechanical processes even though fractures are present as well. Significant stress concentration is observed around cavities. Furthermore, the results of 3D cases imply that phase gravity segregation in cavities leads to lower water saturation in the area above cavities and delays water breakthrough.