Latent seismicity driven by aseismic creep and enhanced pore-fluid pressure in NE British Columbia

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

<p>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<sub>L</sub> 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 (<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.</p>

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.


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

Abstract. Recent seismicity in Alberta and British Columbia has been attributed to ongoing oil and gas development in the area, due to its temporal and spatial correlation. Prior to such development, the area was seismically quiescent. Here, we show evidence that latent seismicity may occur in areas where previous operations may have occurred, even during a shutdown in operations. The global pandemic of COVID-19 furnished the unique opportunity to study seismicity during a period of anthropogenic quiescence. A total of 389 events were detected within the Kiskatinaw area of British Columbia from April to August 2020, which encompasses a period with no hydraulic fracturing operations during a government imposed lockdown. Apart from a reduction in seismicity rate, the general characteristics of the observed seismicity were similar to the preceding time period of active operations. During the shutdown, observed event magnitudes fell between ML −1 and ML 1.2, but lacked temporal clustering that is often characteristic of hydraulic-fracturing induced sequences. Hypocenters occurred in a corridor orientated NW-SE, just as seismicity had done in previous years in the area, and locate at depths associated with the target Montney formation or shallower (


2021 ◽  
Vol 40 (6) ◽  
pp. 413-417
Author(s):  
Chunfang Meng ◽  
Michael Fehler

As fluids are injected into a reservoir, the pore fluid pressure changes in space and time. These changes induce a mechanical response to the reservoir fractures, which in turn induces changes in stress and deformation to the surrounding rock. The changes in stress and associated deformation comprise the geomechanical response of the reservoir to the injection. This response can result in slip along faults and potentially the loss of fluid containment within a reservoir as a result of cap-rock failure. It is important to recognize that the slip along faults does not occur only due to the changes in pore pressure at the fault location; it can also be a response to poroelastic changes in stress located away from the region where pore pressure itself changes. Our goal here is to briefly describe some of the concepts of geomechanics and the coupled flow-geomechanical response of the reservoir to fluid injection. We will illustrate some of the concepts with modeling examples that help build our intuition for understanding and predicting possible responses of reservoirs to injection. It is essential to understand and apply these concepts to properly use geomechanical modeling to design geophysical acquisition geometries and to properly interpret the geophysical data acquired during fluid injection.


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.


Solid Earth ◽  
2021 ◽  
Vol 12 (3) ◽  
pp. 765-783
Author(s):  
Rebecca O. Salvage ◽  
David W. Eaton

Abstract. Recent seismicity in Alberta and north-east British Columbia has been attributed to ongoing oil and gas development in the area, due to its temporal and spatial correlation. Prior to such development, the area was seismically quiescent. Here, we show evidence that latent seismicity may occur in areas where previous operations have occurred, even during a shutdown in operations. The global COVID-19 pandemic furnished the unique opportunity to study seismicity during a long period of anthropogenic quiescence. Within the Kiskatinaw area of British Columbia, 389 events were detected from April to August 2020, which encompasses a period with very little hydraulic fracturing operations. This reduction in operations was the result of a government-imposed lockdown severely restricting the movement of people as well as a downturn in the economic market causing industry stock prices to collapse. Except for a reduction in the seismicity rate and a lack of temporal clustering that is often characteristic of hydraulic fracturing induced sequences, the general characteristics of the observed seismicity were similar to the preceding time period of active operations. During the period of relative quiescence, event magnitudes were observed between ML −0.7 and ML 1.2, which is consistent with previous event magnitudes in the area. Hypocentres occurred in a corridor orientated NW–SE, just as seismicity had done in previous years, and were located at depths associated with the target Montney formation or shallower (<2.5 km). A maximum of 21 % of the detected events during lockdown may be attributable to natural seismicity, with a further 8 % potentially attributed to dynamic triggering of seismicity from teleseismic events and 6 % related to ongoing saltwater disposal and a single operational well pad. However, this leaves ∼65 % of the seismicity detected during lockdown being unattributable to primary activation mechanisms. This seismicity is unlikely to be the result of direct pore pressure increases (as very little direct injection of fluids was occurring at the time) and we see no patterns of temporal or spatial migration in the seismicity as would be expected from direct pore pressure increases. Instead, we suggest that this latent seismicity may be generated by aseismic slip as fluids (resulting from previous hydraulic fracturing injection) become trapped within permeable formations at depth, keeping pore pressures in the area elevated and consequently allowing the generation of seismicity. Alternatively, this seismicity may be the result of fault and fracture weakening in response to previous fluid injection. This is the first time that this latent seismicity has been observed in this area of British Columbia and, as such, this may now represent the new normal background seismicity rate within the Kiskatinaw area.


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 ◽  
Author(s):  
Alessandro Verdecchia ◽  
Bei Wang ◽  
Yajing Liu ◽  
Rebecca Harrington ◽  
Marco Roth ◽  
...  

&lt;p&gt;The Dawson-Septimus area near the towns of Dawson Creek and Fort St. John, British Columbia, Canada has experienced a drastic increase in seismicity in the last ~ 6 years, from no earthquakes reported by Natural Resources Canada (NRCan) prior to 2013 to a total of ~ 200 cataloged events in 2013 &amp;#8211; 2019. The increase follows the extensive horizontal drilling and multistage hydraulic fracturing activity that started to extract shale gas from the unconventional siltstone resource of the Montney Formation. In addition to hydraulic fracturing, ongoing wastewater disposal in the permeable sandstones and carbonates located stratigraphically above and below the Montney formation may also be contributing to elevated seismicity in the region. Earthquakes occur in close spatial and temporal proximity to hydraulic fracturing wells, at distances up to ~ 10 km. The expected diffusion time scales in the low-diffusivity siltstone rock units and the temporal and spatial scale of seismic activity beg questions about the possible processes controlling the location and timing of earthquakes.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;Here, we investigate the causative mechanisms for two of the largest events in the Montney Basin, British Columbia: the August 2015 M4.6 earthquake near Fort St. John, and the November 2018 M4.5 earthquake near Dawson Creek. Both events are thought to have occurred within the crystalline basement, ~2 km below the injected shale units (Montney formation). &amp;#160;We use a finite-element 3D poroelastic model to calculate the coupled evolution of elastic stress and pore pressure due to injection at several hydraulic fracturing stages. Initially, we consider a simple layered model with differing hydraulic parameters based on lithology. Subsequently, also considering the seismicity distribution for each sequence, we introduce hypothetic hydraulic conduits connecting the injection intervals with the crystalline basement, where the respective mainshock occurred. We test a range of permeability values (10&lt;sup&gt;-15&lt;/sup&gt; m&lt;sup&gt;2&lt;/sup&gt;&amp;#8211; 10&lt;sup&gt;-12&lt;/sup&gt; m&lt;sup&gt;2&lt;/sup&gt;) commonly implemented for fault zones.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;Our results show that, for both cases, the poroelastic stress perturbation may be not sufficient to trigger events in the basement. Instead, a scenario with a high-permeability (10&lt;sup&gt;-13&lt;/sup&gt; m&lt;sup&gt;2&lt;/sup&gt;&amp;#8211; 10&lt;sup&gt;-12&lt;/sup&gt; m&lt;sup&gt;2&lt;/sup&gt;) conduits connecting the Montney formation to the fault responsible for the mainshock could better explain the relationship between the hydraulic stimulation and the timing of the two M &gt; 4 earthquakes. For the 2018 M4.5 event, aftershock distribution can be mainly attributed to earthquake-earthquake interaction via static Coulomb stress transfer from the mainshock slip. In addition to the modeling of single well/event sequences, future work will include the long-term poroelastic effect due to multiple disposal wells located in the region.&lt;/p&gt;


2015 ◽  
Vol 52 (2) ◽  
pp. 112-122 ◽  
Author(s):  
Amir Mansour Farahbod ◽  
Honn Kao ◽  
Dan M. Walker ◽  
John F. Cassidy

We systematically re-analyzed historical seismograms to verify the existence of background seismicity in the Horn River Basin of northeast British Columbia before the start of regional shale gas development. We also carefully relocated local earthquakes that occurred between December 2006 and December 2011 to delineate their spatiotemporal relationship with hydraulic fracturing (HF) operations in the region. Scattered seismic events were detected in the Horn River Basin throughout the study periods. The located seismicity within 100 km of the Fort Nelson seismic station had a clearly increasing trend, specifically in the Etsho area where most local HF operations were performed. The number of events was increased from 24 in 2002–2003 (prior to HF operations) to 131 in 2011 (peak period of HF operations). In addition, maximum magnitude of the events was shifted from ML 2.9 to ML 3.6 as the scale of HF operation expanded from 2006–2007 to 2011. Based on our relocated earthquake catalog, the overall b value is estimated at 1.21, which is higher than the average of tectonic/natural earthquakes of ∼1.0. Our observations highly support the likelihood of a physical relationship between HF operation and induced seismicity in the Horn River Basin. Unfortunately, due to the sparse station density in the region, depth resolution is poor for the vast majority of events in our study area. As new seismograph stations are established in northeast British Columbia, both epicentral mislocation and depth uncertainty for future events are expected to improve significantly.


Sign in / Sign up

Export Citation Format

Share Document