A Dynamic Poro-Visco-Elasto-Plastic Earthquake Simulator with Spontaneous Dilatant Coseismic Rupture: Modelling Fault Deformations at Reservoir-Induced Seismicity

Compared to other kinds of fluid-related seismicity, reservoir-induced seismicity (RIS) is usually characterized by higher magnitudes. Seismic and water level monitoring and statistical modelling, however, do not provide comprehensive understanding of the RIS mechanism and controls. This study presents a novel finite element method-based 2D poro-visco-elasto-plastic fully dynamic earthquake model, specifically applicable to RIS simulations. In the first stage of the simulations, Drucker-Prager plasticity is used to generate a normal fault in the Earth’s upper crust with enhanced porosity, over a long time-scale of millions of years. In the second stage of the simulations, RIS is modelled under typical reservoir impoundment dynamics, producing four seismic sequences, triggered by pore pressure increase at the fault at shallow depth below the reservoir. This pressurization is released by aftershocks in every seismic cluster, accompanied by permeability hikes at the fault and associated with fault “valving” behaviour. A dynamic coseismic rupture phase driven by wave-mediated stress transfers coupled with rate-and-state dependent friction coefficient weakening is modelled, along with interseismic deformations. Coseismic crack opening in a dilatant regime, inducing porosity and permeability hikes especially emphasized at the fault, is implemented. The model component verifications demonstrate convincing agreement with theoretical predictions. The model allows investigation of spatio-temporal RIS characteristics and their controls. It may contribute to earthquake prediction in situ and facilitate earthquake mitigation policies.