Stress Chatter via Fluid Flow and Fault Slip in a Hydraulic Fracturing‐Induced Earthquake Sequence in the Montney Formation, British Columbia

2020 ◽  
Vol 47 (14) ◽  
Author(s):  
A. F. Peña Castro ◽  
M. P. Roth ◽  
A. Verdecchia ◽  
J. Onwuemeka ◽  
Y. Liu ◽  
...  
2021 ◽  
Vol 2021 ◽  
pp. 1-21
Author(s):  
Zhiyong Niu ◽  
Shiquan Wang ◽  
Hongrui Ma ◽  
Songbao Feng ◽  
Hengjie Luan ◽  
...  

Hot dry rock (HDR) geothermal energy has many advantages, such as being renewable, clean, widely distributed, and without time and weather limitations. Hydraulic fracturing is usually needed for the exploitation of HDR geothermal energy. It has many hidden faults in the reservoir/caprock sequences. Injecting fluid into underground formations during hydraulic fracturing often induces fault slip and leads to earthquakes. Therefore, to well understand the induced fault slip and earthquakes is important for the applications and development of HDR geothermal exploitation. In this study, we investigated the hazardous injection area of the induced earthquakes during hydraulic fracturing. The study was based on a hydraulic fracturing test in Qiabuqia geothermal field in China. According to the field, a fault-surrounding rock-fracturing region system was developed to study the influences of fluid injection on the stability of the specific fault. A total of 60 hydraulic fracturing regions and 180 numerical experiments were designed. The results revealed that the hazardous injection regions that threaten the fault’s stability were near to the fault and concentrated on the following four areas: (a) above the top of the fault in underlying strata; (b) above the top of the hanging wall of the fault in underlying strata; (c) near to the fault planes in both footwall and hanging wall; (d) at the bottom of the footwall of the fault in underlying strata. The hazardous injection area can be controlled effectively by adjusting the injection pressure.


2020 ◽  
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 (


Sign in / Sign up

Export Citation Format

Share Document