Static behaviour of induced seismicity

16Citations
Citations of this article
26Readers
Mendeley users who have this article in their library.

Abstract

The standard paradigm to describe seismicity induced by fluid injection is to apply non-linear diffusion dynamics in a poroelastic medium. I show that the spatio-temporal behaviour and rate evolution of induced seismicity can, instead, be expressed by geometric operations on a static stress field produced by volume change at depth. I obtain laws similar in form to the ones derived from poroelasticity while requiring a lower description length. Although fluid flow is known to occur in the ground, it is not pertinent to the geometrical description of the spatio-temporal patterns of induced seismicity. The proposed model is equivalent to the static stress model for tectonic foreshocks generated by the Non-Critical Precursory Accelerating Seismicity Theory. This study hence verifies the explanatory power of this theory outside of its original scope and provides an alternative physical approach to poroelasticity for the modelling of induced seismicity. The applicability of the proposed geometrical approach is illustrated for the case of the 2006, Basel enhanced geothermal system stimulation experiment. Applicability to more problematic cases where the stress field may be spatially heterogeneous is also discussed..

Figures

  • Figure 1. Seismicity spatio-temporal behaviour described by the N-C PAST static stress model (tectonic case; Mignan, 2012): (a) spatio-temporal evolution of the stress field σ(r, t) generated by constant stress accumulation τ̇ on a fault located at r = 0 (Eq. 1). Background, quiescence, and activation seismicity regimes are described by densities of events δb0, δbm, and δbp for |σ | ≤ σ ∗ 0 ±1σ∗, σ < σ∗ 0 −1σ∗, and σ > σ∗ 0 +1σ∗, respectively; (b) temporal evolution of quiescence and activation envelopes r∗(t), with σ ( r∗ ) = σ∗ 0 ±1σ∗ (Eq. 2); (c) rate of events µ(t) in a disc of constant radius max ( r∗ ) (Eq. 3); (d) cumulative number of events N(t) (Eq. 4) of power law form (Eq. 5), with t0 = 0, tmid = 1, tf = 2, h= 1, τ̇ = 0.1, σ∗ 0 = 0, 1σ∗ = 10−2, δbm = 0.001, δb0 = 0.1, δbp = 1, n= 3, k = π , d = 2, 1t = 0.01.
  • Figure 2. 2006 Basel EGS stimulation experiment data with activation and quiescence envelope fits: (a) flow rate Q(t) (digitised from Häring et al., 2008); (b) spatio-temporal distribution of relocated induced seismicity (Kraft and Deichmann, 2014) with r the distance from the borehole. The activation and quiescence envelopes r∗ A (t) and r∗ Q (t) are defined from Eq. (13) with parameters 1̂σ∗ = 0.007 day−1 (dark curves) and 1t = 1/4 day. The light curves represent the range 1̂σ∗ ∈ [10−3,10−1] day−1 in 0.1 increments on the log10 scale. Points represent the induced earthquakes; which colour indicates how they are declared. (c) Score S = (wA+wQ)/2, with wA and wQ being the ratio of events of distance r ≤ r∗ A and r ≥ r∗ Q in the injection and bleeding-off phases, respectively. The vertical line represents 1̂σ∗ = 0.007 day−1.
  • Figure 3. Induced seismicity production time series, observed and predicted: (a) histogram of the observed 6 h induced seismicity rate µ(t) with fit based on Eq. (15) with MLE parameters δbp = 4.68× 10−7 eventm−3 day−1 (production parameter) and τ = 1.18 day (diffusion parameter); (b) cumulative number of induced earthquakes N(t) with fit based on Eq. (4) with µ(t) of Eq. (15).
  • Figure 5. Sketch on how anisotropy and other types of heterogeneities can be implemented in the geometrical approach by adding a historical tectonic static stress field (ad hoc parameter values used for sake of simplicity). Here a past overloading field (σ > σ∗ 0 +1σ∗) on a nearby fault would have been “planed” to the threshold σ∗ 0 +1σ∗ by temporal diffusion (Eq. 14), leaving only a “ghost” of that historical static stress field (for the homogeneous case, see Fig. 1a).
  • Figure 4. Description length defined as the count of physical steps required to describe induced seismicity, in poroelasticity and in the newly proposed geometrical approach. In the latter, Biot’s theory is entirely bypassed.

References Powered by Scopus

General theory of three-dimensional consolidation

8716Citations
N/AReaders
Get full text

Power-law distributions in empirical data

6676Citations
N/AReaders
Get full text

Injection-induced earthquakes

2005Citations
N/AReaders
Get full text

Cited by Powered by Scopus

Induced earthquake magnitudes are as large as (statistically) expected

216Citations
N/AReaders
Get full text

Induced seismicity closed-form traffic light system for actuarial decision-making during deep fluid injections

83Citations
N/AReaders
Get full text

Hierarchical Bayesian Modeling of Fluid-Induced Seismicity

43Citations
N/AReaders
Get full text

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Cite

CITATION STYLE

APA

Mignan, A. (2016). Static behaviour of induced seismicity. Nonlinear Processes in Geophysics, 23(2), 107–113. https://doi.org/10.5194/npg-23-107-2016

Readers over time

‘15‘16‘17‘18‘19‘20‘21‘22‘23‘24‘2502468

Readers' Seniority

Tooltip

PhD / Post grad / Masters / Doc 12

60%

Researcher 7

35%

Professor / Associate Prof. 1

5%

Readers' Discipline

Tooltip

Earth and Planetary Sciences 11

65%

Engineering 3

18%

Energy 2

12%

Social Sciences 1

6%

Save time finding and organizing research with Mendeley

Sign up for free
0