Fracture in geological subsurface is a major concern in many underground geoengineering applications, including CO2 sequestration, fossil energy recovery, and underground hydrology engineering. Due to the high cost of experimental measurements associated with geomaterial failure, many geoengineering studies primarily rely on computational modeling and simulations. However, the modeling of two-phase flow-driven fracture in geological media is a critical challenge due to the complex coupling physics of two-phase flow, porous skeleton and fracture, and there is still a lack of effective computational tools and numerical models. In this work, we present a novel scheme to model two-phase flow-driven fracture in geological porous media using the phase-field method. First, we derived the two-phase flow governing equations based on the mass continuity equations and movement balance equations of each fluid phase in fractured and unbroken porous media. After coupling the governing equations with porous media elasticity and phase-field equations, we construct a partial differential equation (PDE) system to model the two-phase flow-driven fracture in porous media. This modeling scheme provides an effective way to study two-phase flow-driven fracture in porous media. The authors would like to acknowledge the DOE-Center for Frontiers of Subsurface Energy Security (CFSES) for providing support on this research.

This content is only available via PDF.
You can access this article if you purchase or spend a download.