Exploring possible causative mechanisms for earthquakes triggered by hydraulic fracturing: examples from the Montney Basin, BC, Canada.

Author(s):  
Alessandro Verdecchia ◽  
Bei Wang ◽  
Yajing Liu ◽  
Rebecca Harrington ◽  
Marco Roth ◽  
...  

<p>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 – 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.</p><p> </p><p>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).  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<sup>-15</sup> m<sup>2</sup>– 10<sup>-12</sup> m<sup>2</sup>) commonly implemented for fault zones.</p><p> </p><p>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<sup>-13</sup> m<sup>2</sup>– 10<sup>-12</sup> m<sup>2</sup>) 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 > 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.</p>

2020 ◽  
Vol 91 (5) ◽  
pp. 2744-2756 ◽  
Author(s):  
Marco P. Roth ◽  
Alessandro Verdecchia ◽  
Rebecca M. Harrington ◽  
Yajing Liu

Abstract The number of earthquakes in the western Canada sedimentary basin (WCSB) has increased drastically in the last decade related to unconventional energy production. The majority of reported earthquakes are correlated spatially and temporally with hydraulic fracturing (HF) well stimulation. In this study, we use waveform data from a new deployment of 15 broadband seismic stations in a spatial area of roughly 60×70km2, covering parts of the Montney Formation, to study the relationship between earthquakes and HF operations in the Dawson-Septimus area, British Columbia, Canada, where the two largest HF-related earthquakes in WCSB to date, an Mw 4.6 on 17 August 2015 and an ML 4.5 on 30 November 2018, have occurred. We use an automated short-term average/long-term average algorithm and the SeisComP3-software to detect and locate 5757 local earthquakes between 1 July 2017 and 30 April 2019. Using two clustering techniques and double-difference relocations of the initial catalog, we define event families that are spatially associated with specific wells, and exhibit temporal migration along a horizontal well bore and/or multiple fractures close to wells. Relocated clusters align in two dominant orientations: one roughly perpendicular to the maximum horizontal regional stress direction (SH) and several conjugate structures at low angles to SH. Comparing the two predominant seismicity lineations to regional earthquake focal mechanisms suggests that deformation occurs via thrust faulting with fault strike oriented perpendicular to SH and via strike-slip faulting with strike azimuth at low angles to SH. Local scale seismicity patterns exhibit clustering around individual HF wells, whereas regional scale patterns form lineations consistent with deformation on faults optimally oriented in the regional stress field.


2021 ◽  
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>


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 (


2019 ◽  
Vol 3 (1) ◽  
pp. 1-14
Author(s):  
Miriam R. Aczel ◽  
Karen E. Makuch

High-volume hydraulic fracturing combined with horizontal drilling has “revolutionized” the United States’ oil and gas industry by allowing extraction of previously inaccessible oil and gas trapped in shale rock [1]. Although the United States has extracted shale gas in different states for several decades, the United Kingdom is in the early stages of developing its domestic shale gas resources, in the hopes of replicating the United States’ commercial success with the technologies [2, 3]. However, the extraction of shale gas using hydraulic fracturing and horizontal drilling poses potential risks to the environment and natural resources, human health, and communities and local livelihoods. Risks include contamination of water resources, air pollution, and induced seismic activity near shale gas operation sites. This paper examines the regulation of potential induced seismic activity in Oklahoma, USA, and Lancashire, UK, and concludes with recommendations for strengthening these protections.


Author(s):  
Bin Chen ◽  
Beatriz Ramos Barboza ◽  
Yanan Sun ◽  
Jie Bai ◽  
Hywel R Thomas ◽  
...  

AbstractAlong with horizontal drilling techniques, multi-stage hydraulic fracturing has improved shale gas production significantly in past decades. In order to understand the mechanism of hydraulic fracturing and improve treatment designs, it is critical to conduct modelling to predict stimulated fractures. In this paper, related physical processes in hydraulic fracturing are firstly discussed and their effects on hydraulic fracturing processes are analysed. Then historical and state of the art numerical models for hydraulic fracturing are reviewed, to highlight the pros and cons of different numerical methods. Next, commercially available software for hydraulic fracturing design are discussed and key features are summarised. Finally, we draw conclusions from the previous discussions in relation to physics, method and applications and provide recommendations for further research.


Sign in / Sign up

Export Citation Format

Share Document