Combined boundary element and finite volume methods for modelling fluid-induced seismicity in fault networks within low-permeability rocks

SUMMARY To better understand the mechanics of injection-induced seismicity, we developed a 2-D numerical code to simulate both seismic and aseismic slip on non-planar faults and fault networks driven by fluid diffusion along permeable faults, in an impervious host rock. Our approach integrates a boundary element method to model fault slip governed by rate-and-state friction with a finite-volume method to simulate fluid diffusion along fault networks. We demonstrate the capabilities of the method with two illustrative examples: (1) Fluid injection inducing slow slip on a primary rough, rate-strengthening fault, which subsequently triggers microseismicity on nearby secondary, smaller faults, and (2) fluid injection on a single fault in a network of intersecting faults, leading to fluid diffusion and reactivation of slip throughout the network. This work highlights the importance of distinguishing between mechanical and hydrological processes in the analysis of induced seismicity, providing a powerful tool for improving our understanding of fault behaviour in response to fluid injection, in particular when a network of geometrically complex faults is involved.