Research Article| May 01, 2018 Interevent Triggering in Microseismicity Induced by Hydraulic Fracturing Samira Maghsoudi; Samira Maghsoudi aComplexity Science Group, Department of Physics and Astronomy, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4, davidsen@phas.ucalgary.cacAlso at Department of Geoscience, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4. Search for other works by this author on: GSW Google Scholar Jordi Baró; Jordi Baró aComplexity Science Group, Department of Physics and Astronomy, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4, davidsen@phas.ucalgary.ca Search for other works by this author on: GSW Google Scholar Alana Kent; Alana Kent bDepartment of Geoscience, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4 Search for other works by this author on: GSW Google Scholar David Eaton; David Eaton bDepartment of Geoscience, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4 Search for other works by this author on: GSW Google Scholar Jörn Davidsen Jörn Davidsen aComplexity Science Group, Department of Physics and Astronomy, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4, davidsen@phas.ucalgary.ca Search for other works by this author on: GSW Google Scholar Bulletin of the Seismological Society of America (2018) 108 (3A): 1133–1146. https://doi.org/10.1785/0120170368 Article history first online: 01 May 2018 Cite View This Citation Add to Citation Manager Share Icon Share Facebook Twitter LinkedIn Email Tools Icon Tools Get Permissions Search Site Citation Samira Maghsoudi, Jordi Baró, Alana Kent, David Eaton, Jörn Davidsen; Interevent Triggering in Microseismicity Induced by Hydraulic Fracturing. Bulletin of the Seismological Society of America 2018;; 108 (3A): 1133–1146. doi: https://doi.org/10.1785/0120170368 Download citation file: Ris (Zotero) Refmanager EasyBib Bookends Mendeley Papers EndNote RefWorks BibTex toolbar search Search nav search search input Search input auto suggest search filter All ContentBy SocietyBulletin of the Seismological Society of America Search Advanced Search Abstract Permeability enhancing treatments such as hydraulic fracturing (HF) induce microseismic events with typical magnitudes in the −3.0 to −0.5 range, although significantly larger induced earthquakes up to 4.7 in moment magnitude have been reported. Diffusion of pore pressure away from the hydraulic fracture system is thought to be a primary controlling mechanism. Understanding other primary or secondary triggering mechanisms during HF is expected to furnish new insights regarding stress, strength of faults, and rupture initiation and propagation. Using novel methods from statistical seismology, we present evidence for the existence of event–event triggering cascades in microseismicity induced by HF. Although background seismicity dominates, we find that these triggering cascades exhibit features that also characterize tectonic aftershock sequences such as the empirical Omori–Utsu relation and the productivity relation. This suggests that the underlying physical earthquake–earthquake triggering mechanisms are similar in both cases, as also observed for other earthquake swarms. The presence of triggering cascades is of direct relevance for optimizing the effectiveness of the stimulation. You do not currently have access to this article.
Abstract The shape of the spatial aftershock decay is sensitive to the triggering mechanism and thus particularly useful for discriminating between static and dynamic stress triggering. For California seismicity, it has been recently recognized that its form is more complicated than typically assumed consisting of three different regimes with transitions at the scale of the rupture length and the thickness of the crust. The intermediate distance range is characterized by a relative small decay exponent of 1.35 previously declared to relate to dynamic stress triggering. We perform comprehensive simulations of a simple clock‐advance model, in which the number of aftershocks is just proportional to the Coulomb‐stress change, to test whether the empirical result can be explained by static stress triggering. Similarly to the observations, the results show three scaling regimes. For simulations adapted to the depths and focal mechanisms observed in California, we find a remarkable agreement with the observation over the whole distance range for a fault distribution with fractal dimension of 1.8, which is shown to be in good agreement with an independent analysis of California seismicity.
Abstract. Crucial to the development of earthquake forecasting schemes is the manifestation of spatiotemporal correlations between earthquakes as highlighted, for example, by the notion of aftershocks. Here, we present an analysis of the statistical relation between subsequent magnitudes for a recently proposed self-similar aftershock rates model of seismicity, whose main distinguishing feature is that of interdependence between trigger and triggered events in terms of a time-varying frequency magnitude distribution. By means of a particular statistical measure, we study the level of magnitude correlations under specific types of time conditioning, explain their provenance within the model framework and show that the type of null model chosen in the analysis plays a pivotal role in the type and strength of observed correlations. Specifically, we show that while the variations in the magnitude distribution can give rise to large trivial correlations between subsequent magnitudes, the non-trivial magnitude correlations are rather minimal. Simulations mimicking Southern California show that these non-trivial correlations cannot be observed at the 3σ-level at the current level of completeness. We conclude that only the time variations in the frequency-magnitude distribution might lead to significant improvements in earthquake forecasting.
The formation and stability of social hierarchies is a question of general relevance. Here, we propose a simple generalized theoretical model for establishing social hierarchy via pair-wise interactions between individuals and investigate its stability. In each interaction or fight, the probability of "winning" depends solely on the relative societal status of the participants, and the winner has a gain of status whereas there is an equal loss to the loser. The interactions are characterized by two parameters. The first parameter represents how much can be lost, and the second parameter represents the degree to which even a small difference of status can guarantee a win for the higher-status individual. Depending on the parameters, the resulting status distributions reach either a continuous unimodal form or lead to a totalitarian end state with one high-status individual and all other individuals having status approaching zero. However, we find that in the latter case long-lived intermediary distributions often exist, which can give the illusion of a stable society. As we show, our model allows us to make predictions consistent with animal interaction data and their evolution over a number of years. Moreover, by implementing a simple, but realistic rule that restricts interactions to sufficiently similar-status individuals, the stable or long-lived distributions acquire high-status structure corresponding to a distinct high-status class. Using household income as a proxy for societal status in human societies, we find agreement over their entire range from the low-to-middle-status parts to the characteristic high-status "tail". We discuss how the model provides a conceptual framework for understanding the origin of social hierarchy and the factors which lead to the preservation or deterioration of the societal structure.
Earth and Space Science Open Archive This preprint has been submitted to and is under consideration at Journal of Geophysical Research - Solid Earth. ESSOAr is a venue for early communication or feedback before peer review. Data may be preliminary.Learn more about preprints preprintOpen AccessYou are viewing the latest version by default [v1]Limited earthquake interaction during a geothermal hydraulic stimulation in Helsinki, FinlandAuthorsGrzegorzKwiatekiDPatriciaMartínez-GarzónJörnDavidseniDPeter EricMaliniDAinoKarjalaineniDMarcoBohnhoffiDGeorgDresenSee all authors Grzegorz KwiatekiDCorresponding Author• Submitting AuthorHelmholtz Centre Potsdam GFZ German Research Centre for GeosciencesiDhttps://orcid.org/0000-0003-1076-615Xview email addressThe email was not providedcopy email addressPatricia Martínez-GarzónHelmholtz Centre Potsdam GFZ German Research for Geosciencesview email addressThe email was not providedcopy email addressJörn DavidseniDUniversity of CalgaryiDhttps://orcid.org/0000-0002-8103-4442view email addressThe email was not providedcopy email addressPeter Eric MaliniDAdvanced Seismic Instrumentation & ResearchiDhttps://orcid.org/0000-0002-5315-9731view email addressThe email was not providedcopy email addressAino KarjalaineniDSt1 OyiDhttps://orcid.org/0000-0002-3279-2239view email addressThe email was not providedcopy email addressMarco BohnhoffiDHelmholtz Centre Potsdam GFZ German Research Centre for GeosciencesiDhttps://orcid.org/0000-0001-7383-635Xview email addressThe email was not providedcopy email addressGeorg DresenGFZ Potsdamview email addressThe email was not providedcopy email address
Short-term forecasting of estimated maximum magnitude ([Formula: see text]) is crucial to mitigate risks of induced seismicity during fluid stimulation. Most previous methods require real-time injection data, which are not always available. This study proposes two deep learning (DL) approaches, along with two data-partitioning methods, that rely solely on preceding patterns of seismicity. The first approach forecasts [Formula: see text] directly using DL; the second incorporates physical constraints by using DL to forecast seismicity rate, which is then used to estimate [Formula: see text]. These approaches are tested using a hydraulic-fracture monitoring dataset from western Canada. We find that direct DL learns from previous seismicity patterns to provide an accurate forecast, albeit with a time lag that limits its practical utility. The physics-informed approach accurately forecasts changes in seismicity rate, but sometimes under- (or over-) estimates [Formula: see text]. We propose that significant exceedance of [Formula: see text] may herald the onset of runaway fault rupture.
Abstract Previous studies of injection-induced earthquake sequences have shown that the maximum magnitude ( M max ) of injection-induced seismicity increases with the net injected volume ( V ); however, different proposed seismic-hazard paradigms predict significantly different values of M max . Using injection and seismicity data from two project areas in northeastern British Columbia, Canada, where hydraulic fracturing induced seismicity was observed, we test the predictive power and robustness of three existing and one novel method to estimate M max . Due to their vastly different values of seismogenic index (Σ), these two project areas represent end-member cases of seismogenic response. Our novel method progressively adjusts the M max forecast under the assumption that each recorded event embodies an incremental release of fluid-induced stress. The results indicate that our method typically provides the lowest upper bound of the tested methods and it is less sensitive to site-specific calibration parameters such as Σ. This makes the novel method appealing for operational earthquake forecasting schemes as a real-time mitigation strategy to manage the risks of induced seismicity.
ABSTRACT Hydraulic fracturing operations are generally associated with induced seismic events with magnitudes that are typically lower than those of tectonic-based earthquakes. Although fluid migration and the resulting pore-pressure increase appear to be the primary driving mechanisms, one must also understand fault activation and the secondary interevent triggering process(es) leading to the generation of aftershock sequences. This requires distinguishing between seismic events that are a direct consequence of fluid injections and those that predominantly arise because of secondary mechanisms such as static or dynamic stress triggering caused by preceding earthquakes. Focusing on a case study in the Montney area in western Alberta, Canada, where both surface array and borehole geophone data are available, we (1) establish the importance of secondary mechanisms and aftershock sequences in this area and (2) show that the statistical properties of the aftershock sequences can allow one to establish whether significant fault activation is at play and whether the activated fault zone was critically prestressed. Specifically, we find evidence for aftershock sequences using different statistical methodologies, although the presence of aftershocks in the studied hydraulic-fracturing-induced context is significantly less dominant compared with tectonic seismicity. The identified aftershock sequences are characterized by long-range spatial aftershock zones, setting them apart from other cases in fluid-driven settings but comparable to tectonic aftershock sequences and indicative of fault activation. Despite this observation, we find overall a strong dominance of smaller triggers for aftershock sequences and an aftershock production associated with large-magnitude events that is significantly lower than what has been observed in tectonic settings. This suggests that the activated fault zone as a whole was not critically prestressed, and the associated initial seismic hazard was significantly lower compared with tectonic cases.
We show that seismic waiting time distributions in California and Iceland have many features in common as, for example, a power‐law decay with exponent α ≈ 1.1 for intermediate and with exponent γ ≈ 0.6 for short waiting times. While the transition point between these two regimes scales proportionally with the size of the considered area, the full distribution is not universal and depends in a non‐trivial way on the geological area under consideration and its size. This is due to the spatial distribution of epicenters which does not form a simple mono‐fractal. Yet, the dependence of the waiting time distributions on the threshold magnitude seems to be universal.