Geomechanics

Modeling Seismicity and Aseismic slip

A geological reservoir, and more generally the first 10 to 30 kilometers of the crust could be seen as a fractured and permeable rock volume (figure 1). The whole system is loaded by different mechanisms (tectonics, fluid injections and/or withdrawals, tides, remote earthquakes) that alter stresses on prexisting faults. Because of these forcings and because of the typical time and slip dependence of rock friction, slip is triggered on these faults, and emplaces either fast enough to radiate elastic waves (earthquake), or slowly so that an silent aseismic slip event occurs.

Slow slip (aseismic) events and earthquakes may also occur on the same particular fault, as shown by many seismological and geodetical studies, suggesting that contact on faults is itself heterogeneous (figure 2).

The deformations produced by slip reactare usually monitored from the earth surface by a set of geophysical sensors (seismometers, creep/sarin-meters, GPS, ...)

Figure 1 : conceptual model of a fractured and permeable geological reservoir

Figure 2 : conceptual model of a heterogeneous fault

In this context, the physical parameters controlling the triggering, the amount and the size of the earthquakes and aseismic slip events is still poorly understood. It is therefore  difficult to estimate the seismic risk associated with human operations at depth. Furthermore, the seismic activity can not be used to infer deep stress change in the crust, or fluid flows.

In this framework, our research activities are mainly dedicated to the understanding of fault (or fault network) slip reactivation under different types of loading and/or mechanical environment. We are in particular interested in the interplay between seismic and aseismic processes in the release of accumulated stress. For that, we develop mechanical models of frictional faults based on the coupling between Dieterich-Ruina rate-and-state friction law and elastic (or poroelastic) interaction. We conduct numerical and theoretical analysis of these simplified models (some examples in figures 3, 4 and 5).

Figure 4 : Map view of the propagation of a quasi-dynamic rupture (earthquake) on a set of 4 coplanar seismogenic asperities (black circles) surrounded by frictionally stable barriers. Colors represent normalized slip rate on the fault (v/vp). Results were obtained with the quasi-dynamic rate-and-state asperity model developped by (Dublanchet et al., 2013).

Figure 5 : Simulated slip accumulation on a heterogeneous rate-and-state frictional crack undergoing constant remote tectonic loading. The grey area corresponds to a major earthquake (mainshock E). Its occurrence is preceded by a sequence of small foreshocks and accelerated aseismic slip (nucleation phase). Modified from (Dublanchet, GJI, 2018).

Figure 3 : Simulated slip rate ransient on an intfinite velocity-strengthening rate-and-state fault triggered by a localized fluid injection (at x=0). Darker profiles correspond to later times after the injection. Slip rate on the fault is first excited around the injection. While decaying in amplitude, he slip rate perturbation expands along strike. Modified from (Dublanchet & Viesca, 2017).

Main references :
  • Dublanchet, P., The dynamics of earthquake precursors controlled by effective friction, Geophysical Journal International,, Vol. 112, Issue 2, 2018
  • Almakari, M., Dublanchet, P., Chauris, H., Injection-induced seismicity controlled by the pore-pressure rate, EGU General Assembly, 2018
  • Dublanchet, P., Viesca, R.C. Dynamics of Fluid Induced Aseismic Slip, 79th EAGE Conference and Exhibition, Paris, 2017
  • Dublanchet, P., Bernard, P., Favreau, P. (2013), Interactions and triggering in a 3D rate-and-state, asperity model, Journal of Geophysical Research : Solid Earth, 118(5) : 2225-2245, 2013