A fully coupled model for predicting geomechanical and multiphase flow behaviour in fractured rocks

Geomechanical and multiphase flow characteristics are essential in recovering oil from naturally fractured rocks during hydrocarbon production because of changes in pore pressure and tension within the rock. It is a wellestablished fact that the geomechanical and multiphase flow characteristics of fractured rocks are interdependent on each other. Evaluation of these characteristics, for hydrocarbons displaced by water in fractured rocks under external stress loading, is severely lacking in published literature. This study aims to develop a novel numerical framework for a fully coupled model of fractured rocks, taking into consideration the pore pressure and porous media discontinuity at the fracture-matrix interface, along with an expanded Darcy's equation. The fully coupled Finite Element Method (FEM) and Computational Fluid Dynamics (CFD) model developed in this study is shown to accurately predict geomechanical and multiphase flow behaviour at the fracture-matrix interface. The results show that as external stress loading on the fractured rock increases, the porosity and permeability of the rock matrix decrease, capillary pressure at the fracture-matrix interface decreases, and the relative permeability curves shift to the right, indicating a water-soaked fracture-matrix interface. The findings of this study can be used to develop innovative strategies for enhanced oil recovery from fractured rocks.