Effect of the Injection Scenario on the Rate and Magnitude Content of Injection-Induced Seismicity

in collaboration with Pierre Dublanchet and Herve Chauris (MINES ParisTech)

Injection of fluids into underground formations reactivates preexisting faults and modifies the seismic hazard, as demonstrated by the 2011 Mw 5.7 and the 2016 Mw 5.8 earthquakes in Oklahoma. Currently, the effect of injection remains poorly understood. We model the seismicity triggered by a fluid flowing inside a Dietrich-Ruina heterogeneous 2-D fault, which can generate irregular sequences of events with magnitudes obeying Gutenberg Richter distribution. We consider a punctual injection scenario where injection pressure increases at a constant rate until a maximum pressure is reached and kept constant. We show that such a fluid injection leads to a sharp increase in the seismicity rate, which correlates with the time series of the pore pressure rate, for a wide range of injection pressure. Increasing the final pressure leads to an increase in the amplitude and the duration of the seismicity rate perturbation but also to a decrease in the frequency of large-magnitude events. The maximum seismicity rate during the sequence also increases with the injection pressure rate, as long as a pressure-rate threshold is not exceeded. Beyond it, the effect of increasing the injection rate is to make large-magnitude earthquakes more frequent. While the total number of induced earthquakes is essentially controlled by the maximum pressure, the total seismic moment liberated increases with both the maximum pressure and the pressure rate. The comparison of our model to Dietrich’s (1994, https://doi.org/10.1029/93JB02581) model shows the important trade-off existing between seismicity rate perturbations and magnitude content variations of fluid induced seismicity.

Time series of the seismicity rate (black curve), mean pore pressure, pore pressure at injection point, and mean pore pressure rate, on a semilogarithmic scale. The black horizontal dashed line represents the initial seismicity rate S0 (before injection). Dotted lines represent the point (tmax,Smax), where the seismicity rate reaches its maximum. The time axis is normalized from the injection start time ti.
Seismicity rate increase Smax∕S0 (ratio of maximum seismicity rate over seismicity rate before injection) as a function of the injection pore pressure rate 𝛽. Continuous lines represent the results of the numerical model of the heterogeneous fault; dashed lines represent the seismicity rate increase predicted by analytical approximation based on Dietrich's 1994 model

Corresponding publication: Almakari at el, 2019