scholarly journals Maturity of nearby faults influences seismic hazard from hydraulic fracturing

2018 ◽  
Vol 115 (8) ◽  
pp. E1720-E1729 ◽  
Author(s):  
Maria Kozłowska ◽  
Michael R. Brudzinski ◽  
Paul Friberg ◽  
Robert J. Skoumal ◽  
Nicholas D. Baxter ◽  
...  

Understanding the causes of human-induced earthquakes is paramount to reducing societal risk. We investigated five cases of seismicity associated with hydraulic fracturing (HF) in Ohio since 2013 that, because of their isolation from other injection activities, provide an ideal setting for studying the relations between high-pressure injection and earthquakes. Our analysis revealed two distinct groups: (i) deeper earthquakes in the Precambrian basement, with larger magnitudes (M > 2), b-values < 1, and many post–shut-in earthquakes, versus (ii) shallower earthquakes in Paleozoic rocks ∼400 m below HF, with smaller magnitudes (M < 1), b-values > 1.5, and few post–shut-in earthquakes. Based on geologic history, laboratory experiments, and fault modeling, we interpret the deep seismicity as slip on more mature faults in older crystalline rocks and the shallow seismicity as slip on immature faults in younger sedimentary rocks. This suggests that HF inducing deeper seismicity may pose higher seismic hazards. Wells inducing deeper seismicity produced more water than wells with shallow seismicity, indicating more extensive hydrologic connections outside the target formation, consistent with pore pressure diffusion influencing seismicity. However, for both groups, the 2 to 3 h between onset of HF and seismicity is too short for typical fluid pressure diffusion rates across distances of ∼1 km and argues for poroelastic stress transfer also having a primary influence on seismicity.

2020 ◽  
Vol 91 (6) ◽  
pp. 3171-3181 ◽  
Author(s):  
Maomao Wang ◽  
Hongfeng Yang ◽  
Lihua Fang ◽  
Libo Han ◽  
Dong Jia ◽  
...  

Abstract Human activity-induced earthquakes are emerging as a global issue, and revealing its underlying mechanisms is essential for earthquake hazard mitigation and energy development. We investigated the relationship between the seismotectonic model and seismic sequences from moderate Mw 4.3 and Mw 5.2 earthquakes that occurred in February and September 2019, respectively, in the Weiyuan anticline of Sichuan basin, China. We found that the Mw 5.2 earthquake ruptured a back thrust of structural wedges and released most aftershocks near the wedge tip. However, the two foreshocks of the Mw 4.3 earthquake sequence occurred in hydrofractured Silurian shale at depth of 2.5–3 km, and the mainshock ruptured the overlying oblique tear fault at a depth of ∼1  km. Hydraulic fracturing in the sedimentary cover of this block may induce earthquakes through fluid pressure diffusion in the Silurian shale and through poroelastic effects on back thrusts within structural wedges, respectively. We assessed the hazard potential of four seismic sources in the Weiyuan block and suggest it is critical to conduct a coupled flow-geomechanics assessment and management on induced seismicity and related cascading effects in the densely inhabited and seismically active Sichuan basin.


2020 ◽  
Vol 91 (6) ◽  
pp. 3182-3194
Author(s):  
Hongfeng Yang ◽  
Pengcheng Zhou ◽  
Nan Fang ◽  
Gaohua Zhu ◽  
Wenbin Xu ◽  
...  

Abstract Earthquakes rarely occur at extremely shallow depths, for example, less than 2 km. Even for induced earthquakes that are typically shallower than tectonic events, only very small ones have been reported in such depths. The ML 4.9 earthquake (Mw 4.3) that struck the Rongxian County, Sichuan, China on 25 February 2019 was an extremely shallow event. Seismological and geodetic data constrained the mainshock depth at ∼1  km with a thrust-faulting mechanism, consistent with the Molin fault orienting northwest. Two foreshocks with magnitudes larger than 4 occurred on an unmapped fault striking northeast, right next to an injection well where hydraulic fracturing (HF) was conducted. The focal depths of the two foreshocks were at ∼2.7  km, coinciding with the depth of HF. Coulomb failure stresses of the two foreshocks on the Molin fault was ∼3  kPa, smaller than typical static triggering threshold (10 kPa), and thus their triggering effects were mild. As the fault was hydraulically sealed from HF, we suggested that the ML 4.9 earthquake was possibly triggered by nearby HF activities through poroelastic stress transfer. Such findings held significant implications for shale gas development by considering seismic hazard associated with shallow faults.


2021 ◽  
Vol 9 ◽  
Author(s):  
Tom Kettlety ◽  
James P. Verdon

We investigate the physical mechanisms governing the activation of faults during hydraulic fracturing. Recent studies have debated the varying importance of different fault reactivation mechanisms in different settings. Pore pressure increase caused by injection is generally considered to be the primary driver of induced seismicity. However, in very tight reservoir rocks, unless a fracture network exists to act as a hydraulic conduit, the rate of diffusion may be too low to explain the spatio-temporal evolution of some microseismic sequences. Thus, elastic and poroelastic stress transfer and aseismic slip have been invoked to explain observations of events occurring beyond the expected distance of a reasonable diffusive front. In this study we use the high quality microseismic data acquired during hydraulic fracturing at the Preston New Road (PNR) wells, Lancashire, UK, to examine fault triggering mechanisms. Injection through both wells generated felt induced seismicity—an ML 1.6 during PNR-1z injection in 2018 and an ML 2.9 during PNR-2 in 2019—and the microseismic observations show that each operation activated different faults with different orientations. Previous studies have already shown that PNR-1z seismicity was triggered by a combination of both direct hydraulic effects and elastic stress transfer generated by hydraulic fracture opening. Here we perform a similar analysis of the PNR-2 seismicity, finding that the PNR-2 fault triggering was mostly likely dominated by the diffusion of increased fluid pressure through a secondary zone of hydraulic fractures. However, elastic stress transfer caused by hydraulic fracture opening would have also acted to promote slip. It is significant that no microseismicity was observed on the previously activated fault during PNR-2 operations. This dataset therefore provides a unique opportunity to estimate the minimum perturbation required to activate the fault. As it appears that there was no hydraulic connection between them during each stimulation, any perturbation caused to the PNR-1z fault by PNR-2 stimulation must be through elastic or poroelastic stress transfer. As such, by computing the stress transfer created by PNR-2 stimulation onto the PNR-1z fault, we are able to approximate the minimum bound for the required stress perturbation: in excess of 0.1 MPa, orders of magnitude larger than stated estimates of a generalized triggering threshold.


2021 ◽  
Author(s):  
Rebecca O. Salvage ◽  
David W. Eaton

&lt;p&gt;The global pandemic of COVID-19 furnished an opportunity to study seismicity in the Kiskatinaw area of British Columbia, noted for hydraulic-fracturing induced seismicity, during a period of anthropogenic quiescence. A total of 389 events were detected from April to August 2020, encompassing a period with no hydraulic-fracturing operations during a government-imposed lockdown. During this time period, observed seismicity had a maximum magnitude of M&lt;sub&gt;L&lt;/sub&gt; 1.2 and lacked temporal clustering that is often characteristic of hydraulic-fracturing induced sequences. Instead, seismicity was persistent over the lockdown period, similar to swarm-like seismicity with no apparent foreshock-aftershock type sequences. Hypocenters occurred within a corridor orientated NW-SE, just as seismicity had done in previous years in the area, with focal depths near the target Montney formation or shallower (&lt;2.5 km). Based on the Gutenberg-Richter relationship, we estimate that a maximum of 21% of the detected events during lockdown may be attributable to natural seismicity, with a further 8% possibly due to dynamic triggering of seismicity from teleseismic events. The remaining ~70% cannot be attributed to direct pore pressure increases induced by fluid injection, and therefore is inferred to represent latent seismicity i.e. seismicity that occurs after an unusually long delay following primary activation processes, with no obvious triggering mechanism. We can exclude pore-pressure diffusion from the most recent fluid injection, as is there is no clear pattern of temporal or spatial seismicity migration. If elevated pore pressure from previous injections became trapped in the subsurface, this could explain the localization of seismicity within an operational corridor, but it does not explain the latency of seismicity on a timescale of months. However, aseismic creep on weak surfaces such as faults, in response to tectonic stresses, in addition to trapped elevation pore-pressure could play a role in stress re-loading to sustain the observed pattern of seismicity.&lt;/p&gt;


2019 ◽  
Author(s):  
Sergey Turuntaev ◽  
Evgeny Zenchenko ◽  
Maria Trimonova ◽  
Petr Zenchenko ◽  
Nikolay Baryshnikov ◽  
...  

2021 ◽  
Author(s):  
Georgios Michas ◽  
Vasilis Kapetanidis ◽  
George Kaviris ◽  
Filippos Vallianatos

&lt;p&gt;Earthquake diffusion is frequently observed in the spatiotemporal evolution of seismic clusters and regional seismicity, a characteristic that is attributed to a triggering mechanism, such as fluid flow, aseismic creep and/or stress transfer effects. In this work, we study the earthquake diffusion properties in the Western Gulf of Corinth (central Greece), an area that presents high extension rates, moderate to large magnitude earthquakes, intense microseismicity and frequent seismic swarms. We focus on the period 2013&amp;#8211;2014 that is characterized by intense background microseismic activity along with significant seismic sequences. More specifically, the latter include the 2013 Helike swarm, the 2014 seismic sequence between Nafpaktos and Psathopyrgos, which culminated with an Mw 4.9 event on 21 September 2014, as well as moderate magnitude events that were followed by aftershock sequences. In the herein analysis, we employ a relocated earthquake catalogue of ~9000 events which delineates the activated areas during the study period in high-resolution. We consider the most significant seismic sequences and calculate their respective spatial correlation histograms and the evolution of the mean squared distance of the hypocenters with time, in order to study the earthquake diffusion rates and possible variations that might be related to the triggering mechanisms of seismicity. Our results demonstrate a weak earthquake diffusion process, analogous to subdiffusion within a stochastic framework, for the seismic sequences under consideration, providing further evidence for slow earthquake diffusion in regional and global seismicity. In addition, the earthquake diffusion rates exhibit variations that can be associated with the triggering mechanism. In particular, seismic sequences which are related with pore-fluid pressure diffusion present considerably higher diffusion rates than mainshock/aftershock sequences or the background activity. Such results may provide novel constraints on the triggering mechanisms of clustered seismic activity based on the study of the earthquake diffusion rates.&amp;#160;&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Acknowledgements&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;We would like to thank the personnel of the Hellenic Unified Seismological Network (http://eida.gein.noa.gr/) and the Corinth Rift Laboratory Network (https://doi.org/10.15778/RESIF.CL) for the installation and operation of the stations used in the current article. The present research is co-financed by Greece and the European Union (European Social Fund- ESF) through the Operational Programme &amp;#171;Human Resources Development, Education and Lifelong Learning 2014-2020&amp;#187; in the context of the project &amp;#8220;The role of fluids in the seismicity of the Western Gulf of Corinth (Greece)&amp;#8221; (MIS 5048127).&lt;/p&gt;


SPE Journal ◽  
2021 ◽  
pp. 1-12
Author(s):  
Gang Hui ◽  
Shengnan Chen ◽  
Zhangxin Chen ◽  
Fei Gu ◽  
Mathab Ghoroori ◽  
...  

Summary The relationships among formation properties, fracturing operations, and induced earthquakes nucleated at distinctive moments and positions remain unclear. In this study, a complete data set on formations, seismicity, and fracturing treatments is collected in Fox Creek, Alberta, Canada. The data set is then used to characterize the induced seismicity and evaluate its susceptibility toward fracturing stimulations via integration of geology, geomechanics, and hydrology. Five mechanisms are identified to account for spatiotemporal activation of the nearby faults in Fox Creek, where all major events [with a moment magnitude (Mw) greater than 2.5] are caused by the increase in pore pressure and poroelastic stress during the fracturing operation. In addition, an integrated geological index (IGI) and a combined geomechanical index (CGI) are first proposed to indicate seismicity susceptibility, which is consistent with the spatial distribution of induced earthquakes. Finally, mitigation strategy results suggest that enlarging a hydraulic fracture-fault distance and decreasing a fracturing job size can reduce the risk of potential seismic activities.


2021 ◽  
Author(s):  
Samuel Chapman ◽  
Jan V. M. Borgomano ◽  
Beatriz Quintal ◽  
Sally M. Benson ◽  
Jerome Fortin

&lt;p&gt;Monitoring of the subsurface with seismic methods can be improved by better understanding the attenuation of seismic waves due to fluid pressure diffusion (FPD). In porous rocks saturated with multiple fluid phases the attenuation of seismic waves by FPD is sensitive to the mesoscopic scale distribution of the respective fluids. The relationship between fluid distribution and seismic wave attenuation could be used, for example, to assess the effectiveness of residual trapping of carbon dioxide (CO2) in the subsurface. Determining such relationships requires validating models of FPD with accurate laboratory measurements of seismic wave attenuation and modulus dispersion over a broad frequency range, and, in addition, characterising the fluid distribution during experiments. To address this challenge, experiments were performed on a Berea sandstone sample in which the exsolution of CO2 from water in the pore space of the sample was induced by a reduction in pore pressure. The fluid distribution was determined with X-ray computed tomography (CT) in a first set of experiments. The CO2 exosolved predominantly near the outlet, resulting in a heterogeneous fluid distribution along the sample length. In a second set of experiments, at similar pressure and temperature conditions, the forced oscillation method was used to measure the attenuation and modulus dispersion in the partially saturated sample over a broad frequency range (0.1 - 1000 Hz). Significant P-wave attenuation and dispersion was observed, while S-wave attenuation and dispersion were negligible. These observations suggest that the dominant mechanism of attenuation and dispersion was FPD. The attenuation and dispersion by FPD was subsequently modelled by solving Biot&amp;#8217;s quasi-static equations of poroelasticity with the finite element method. The fluid saturation distribution determined from the X-ray CT was used in combination with a Reuss average to define a single phase effective fluid bulk modulus. The numerical solutions agree well with the attenuation and modulus dispersion measured in the laboratory, supporting the interpretation that attenuation and dispersion was due to FPD occurring in the heterogenous distribution of the coexisting fluids. The numerical simulations have the advantage that the models can easily be improved by including sub-core scale porosity and permeability distributions, which can also be determined using X-ray CT. In the future this could allow for conducting experiments on heterogenous samples.&lt;/p&gt;


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. KS105-KS118 ◽  
Author(s):  
Himanshu Barthwal ◽  
Mirko van der Baan

Hydraulic fracturing in low-permeability hydrocarbon reservoirs creates/reactivates a fracture network leading to microseismic events. We have developed a simplified model of the evolution of the microseismic cloud based on the opening of a planar fracture cavity and its effect on elastic stresses and pore pressure diffusion during fluid injection in hydraulic fracturing treatments. Using a material balance equation, we compute the crack tip propagation over time assuming that the hydraulic fracture is shaped as a single penny-shaped cavity. Results indicate that in low-permeability formations, the crack tip propagates much faster than the pore pressure diffusion front thereby triggering the microseismic events farthest from the injection domain at any given time during fluid injection. We use the crack tip propagation to explain the triggering front observed in distance versus time plots of published microseismic data examples from hydraulic fracturing treatments of low-permeability hydrocarbon reservoirs. We conclude that attributing the location of the microseismic triggering front purely to pore pressure diffusion from the injection point may lead to incorrect estimates of the hydraulic diffusivity by multiple orders of magnitude for low-permeability formations. Moreover, the opening of the fracture cavity creates stress shadow zones perpendicular to the principal fracture walls in which microseismic triggering due to the elastic stress perturbations is suppressed. Microseismic triggering in this stress shadow region may be attributed mainly to pore pressure diffusion. We use the width, instead of the longest size, of the microseismic cloud to obtain an enhanced diffusivity measure, which may be useful for subsequent production simulations.


Sign in / Sign up

Export Citation Format

Share Document