We explore the role of earthquake interactions during an injection‐induced seismic sequence. We propose a model, which considers both a transient pressure and static stress redistribution due to event interactions as triggering mechanisms. By calibrating the model against observations at the Enhanced Geothermal System of Basel, Switzerland, we are able to reproduce the time behavior of the seismicity rate. We observe that considering earthquake interactions in the modeling leads to a larger number of expected seismic events (24% more) if compared to a pressure‐induced seismicity only. The increase of the number of events is particularly evident after the end of the injection. We conclude that implementing a model for estimating the static stress changes due to mutual event interactions increases significantly the understanding of the process and the behavior of induced seismicity.

No abstract

The 2012 Emilia Romagna (Italy) seismic sequence has been extensively studied given the occurrence of two mainshocks, both temporally and spatially close to each other. The recent literature accounts for several fault models, obtained with different inversion methods and different datasets. Several authors investigated the possibility that the second event was triggered by the first mainshock with elusive results. In this work, we consider all the available InSAR and GPS datasets and two planar fault geometries, which are based on both seismological and geological constraints. We account for a layered, elastic half-space hosting the dislocation and compare the slip distribution resulting from the inversion and the related changes in Coulomb Failure Function (CFF) obtained with both a homogeneous and layered half-space. Finally, we focus on the interaction between the two main events, discriminating the contributions of coseismic and early postseismic slip of the mainshock on the generation of the second event and discuss the spatio-temporal distribution of the seismic sequence. When accounting for both InSAR and GPS geodetic data we are able to reproduce a detailed coseismic slip distribution for the two mainshocks that is in accordance with the overall aftershock seismicity distribution. Furthermore, we see that an elastic medium with depth dependent rigidity better accounts for the lack of the shallow seismicity, amplifying, with respect to the homogeneous case, the mechanical interaction of the two mainshocks.

Understanding the injection-induced triggering mechanism is a fundamental step towards controlling the seismicity generated by deep underground exploitation. Here we propose a modeling approach based on coupling the TOUGH2 simulator with a geomechanical-stochastic model. The hydromechanical-stochastic model provides a good representation of different mechanisms influencing each other during and after the injection phase. Each mechanism affects the induced seismicity in a different way and at different times during the reservoir stimulation, confirming that a complex interaction is in place, and that more sophisticated and physics-based approaches coupled with statistical model are required to explain such a complex interaction. In addition to previous statistical and hybrid models, our approach accounts for a full 3D formulation of both stresses and fluid flow, further including all the TOUGH2 capabilities. Furthermore, it includes interactions between triggered seismic events through calculation of static stress transfer. In this work, we present the main capabilities of TOUGH2-SEED and apply the model to the Basel EGS case, successfully reproducing the injection pressure as well as the evolution of the seismicity.