Abstract:Summary Earthquakes can be triggered by fluid injection into underground formations. Fluid injection can cause large changes in the underground volume that exert stresses on nearby preexisting faults, leading to seismic activity. Assuming an increase in underground development activities in the future, our understanding of the mechanism underlying induced seismicity must be improved, and methods must be developed to properly assess the risk of seismic events. The objective of this study is to develop a seismicity prediction model that calculates the magnitude and timing of triggered earthquakes or seismic events occurring during various subsurface fluid injection activities. We developed an injection-induced seismicity analysis model that predicts the dynamic earthquake nucleation caused by changes in stress and pore pressure that occur during various subsurface activities. The governing equations consisting of the dynamic motion of the poroelastic spring-slider system, rate and state friction laws and pore pressure diffusion equation were solved using the embedded semi-implicit Runge-Kutta (SIRK) method. The dynamic sliding model was also incorporated into the finite element method (FEM) model, considering the variations in the stresses and pore pressures in the formation. A field case study was also conducted to compare the model results with typical microseismicity responses observed from hydraulic fracturing treatments in shale fields. Contrary to the popular understanding derived from Amonton’s law, the dynamic friction model revealed that a large normal stress on the fault leads to rapid sliding. A larger normal stress accumulates a large amount of elastic energy until it slips owing to fluid injection, nucleating large seismic waves. The poroelastic spring-slider model estimated reasonable microseismic magnitudes for hydraulic fracturing treatment but overestimated the time required to trigger a microseismic event under field conditions. To improve the analysis results, the poroelastic spring-slider model was coupled with a linear elastic FEM that considered the complex interplay of stress changes from hydraulic fracturing and the associated pore pressure variation in the formation. Compared with the field data, the coupled simulation model estimated a reasonable timing for the induced microseismic events when the increasing pore pressure during hydraulic fracturing penetrated deep into the formation. These findings suggest the existence of permeable natural fractures in the formation, which intensify early frictional sliding during treatment. The seismicity prediction model presented in this study simulates the magnitude and timing of seismic nucleation, helping to manage and mitigate the environmental impacts of induced seismicity during various subsurface development activities, such as oil and gas extraction, hydraulic fracturing, geothermal, and carbon dioxide sequestration. Moreover, the case study results imply that the time series of seismic events predicted by the model can be used to understand the possible fracture geometry and extent of fluid invasion for field applications. Key Terms and Phrases induced seismicity, subsurface fluid injection, rate- and state-dependent friction law, embedded semi-implicit Runge-Kutta method, finite element analysis

Adjoint-based uncertainty quantification for inhomogeneous friction on a slow-slipping fault

Seismic Versus Aseismic Slip: Probing Mechanical Properties of the Northeast Japan Subduction Zone

Uncertainty quantification for fault slip inversion

Numerical experiment toward simultaneous estimation of the spatiotemporal evolution of interseismic fault slips and block motions in southwest Japan

Adjoint-based inversion for stress and frictional parameters in earthquake modeling

Data assimilation for fault slip monitoring and short-term prediction of spatio-temporal evolution of slow slip events: application to the 2010 long-term slow slip event in the Bungo Channel, Japan

Uncertainty Quantification in Friction Model for Earthquakes using Bayesian inference

A Two-Step Inversion for Fault Frictional Properties Using a Temporally Varying Afterslip Model and Its Application to the 2019 Ridgecrest Earthquake

A model for possible crustal deformation prior to a coming large interplate earthquake in the Tokai district, central Japan

Analysis of Fluid-Injection-Induced Seismicity Using a Dynamic Sliding Model Incorporating the Rate- and State-Dependent Friction Law

Spatially variable fault friction derived from dynamic modeling of aseismic afterslip due to the 2004 Parkfield earthquake

Estimating Fault Friction from Seismic Signals in the Laboratory

Patterns of Causative Faults of Normal Earthquakes in the Fluid‐Rich Outer Rise of Northeastern Japan, Constrained With 3D Teleseismic Waveform Modeling

Fluid-driven slow slip and earthquake nucleation on a slip-weakening circular fault

Slip rate deficit and earthquake potential on shallow megathrusts

Estimating the Physical State of a Laboratory Slow Slipping Fault from Seismic Signals

Fluids control along-strike variations in the Alaska megathrust slip

Verification of Theoretical Evaluation of Near-Fault Seismic Ground Motions by Difference of Sub-faults Division Based on Thin Layer Method

Shallow slip amplification and enhanced tsunami hazard unravelled by dynamic simulations of mega-thrust earthquakes

Stress‐driven Relaxation of Heterogeneous Upper Mantle and Time‐dependent Afterslip Following the 2011 Tohoku Earthquake

Segmentation of Shallow Slow Slip Events at the Hikurangi Subduction Zone Explained by Along-Strike Changes in Fault Geometry and Plate Convergence Rates