Numerical Modeling of Fracture Network Evolution in Organic‐Rich Shale With Rapid Internal Fluid Generation

Abstract

When low‐permeability and organic‐rich rocks such as shale experience sufficient heating, chemical reactions including shale dehydration and maturation of organic matter lead to internal fluid generation. This may cause substantial pore fluid overpressure and fracturing. In the vicinity of igneous intrusions emplaced in organic‐rich shales, temperatures of several hundred degrees accelerate these processes and lead to intense fracturing. The resulting fracture network provides hydraulic pathways, which allow fluid expulsion and affect hydrothermal fluid flow patterns. However, the evolution of these complex fracture networks and controls on geometry and connectivity are poorly understood. Here, we perform a numerical modeling study based on the extended finite element method to investigate coupled hydromechanical fracture network evolution due to fast internal fluid generation. We quantify the evolution of different initial fracture networks under varying external stresses by analyzing parameters including fracture length, opening, connectivity, and propagation angles. The results indicate a three‐phase process including (1) individual growth, (2) interaction, and (3) expulsion phase. Magnitude of external stress anisotropy and degree of fracture alignment with the largest principal stress correlate with increased fracture opening. We additionally find that although the external stress field controls the overall fracture orientation distribution, local stress interactions may cause significant deviations of fracture paths and control the coalescence characteristics of fractures. Establishing high connectivity in cases with horizontally aligned initial fractures requires stress anisotropy with σV > σH, while the initial orientation distribution is critical for connectivity if stresses are nearly isotropic.