Toward understanding focal mechanism of hydraulic fracturing induced earthquakes using constrained inversion: method and synthetic tests

2015 ◽  
Author(s):  
Shuhei Iida ◽  
Ahyi Kim
Geophysics ◽  
2020 ◽  
pp. 1-74
Author(s):  
Han Li ◽  
Xu Chang ◽  
Xiao-Bi Xie ◽  
Yibo Wang

Through the study of microseismic focal mechanisms, information such as fracture orientation, event magnitude, and in-situ stress status can be quantitatively obtained, thus, providing a reliable basis for unconventional oil and gas exploration. Most source inversion methods assume that the medium is isotropic. However, hydraulic fracturing is usually conducted in sedimentary rocks, which often exhibit strong anisotropy. Neglecting this anisotropy may cause errors in focal mechanism inversion results. We propose a microseismic focal mechanism inversion method that considers velocity anisotropy in a vertically transverse isotropic (VTI) medium. To generate synthetic data, we adopt the moment-tensor model to represent microearthquake sources. We use a staggered-grid finite-difference (SGFD) method to calculate synthetic seismograms in anisotropic media. We perform seismic moment-tensor (SMT) inversion with only P-waves by matching synthetic and observed waveforms. Both synthetic and field datasets are used to test the inversion method. For the field dataset, we investigate the inversion stability using randomly selected partial datasets in the calculation. We pay special attention to analyze the sensitivity of the inversion. We test and evaluate the impact of noise in the data and errors in the model parameters ( VP0, ε, and δ) on the SMT inversion using synthetic datasets. The results indicate that for a surface acquisition system, the proposed method can tolerate moderate noise in the data, and deviations in the anisotropy parameters can cause errors in the SMT inversion, especially for dip-slip events and the inverted percentages of non-double-couple components. According to our study, including anisotropy in the model is important to obtain reliable non-double-couple components of moment tensors for hydraulic fracturing induced microearthquakes.


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.


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.


2017 ◽  
Author(s):  
Hongyu Zhai* ◽  
Xu Chang ◽  
Yibo Wang ◽  
Ziqiu Xue ◽  
Yi Zhang ◽  
...  

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. KS169-KS182 ◽  
Author(s):  
Xiong Zhang ◽  
Jie Zhang

Similar to a web search engine, we have developed a microseismic search engine that can estimate an event location and the focal mechanism in less than a second to monitor the hydraulic fracturing process. The method was extended from a real-time earthquake monitoring approach for seismological applications. We first calculate the full waveforms of all possible microseismic events over a 3D grid with a known velocity model for a given acquisition geometry to create a database. We then index and rank all of the seismic waveforms in the database by following the characteristics of the phase and amplitude of the waveform through a computer fast search technology, specifically, the multiple randomized k-dimensional tree method. When a microseismic event occurs, the approximate best matches to the entry waveform are found immediately by comparing the characteristic features between the input data and the database. The method returns not just one but a series of solutions, similar to a web search engine. Thus, we can obtain a solution space that delineates the resolution and confidence level of the results. Also similar to a web search engine, the microseismic search engine does not require any input parameter or processing experience; thus, the solutions are the same for any user. Numerical tests suggest that the waveform search approach is insensitive to random and correlated noises. However, if the correlation values between the input data and best matches in the database are too low, suggesting unreliable results, the solution may be rejected automatically by applying a preset threshold. We have applied the method to real data, and found great potential for the routine real-time monitoring of microseismic events during hydraulic fracturing.


2021 ◽  
Vol 9 ◽  
Author(s):  
Min Zhao ◽  
Feng Long ◽  
Guixi Yi ◽  
MingJian Liang ◽  
Jiangtao Xie ◽  
...  

The 3 February 2020 MS 5.1 Qingbaijiang earthquake, southwestern China, is the closest recorded MS ≥ 5.0 event to downtown Chengdu City to date, with an epicentral distance of only 38 km. Here we analyze seismic data from the Sichuan and Chengdu regional seismic networks, and employ a multi-stage location method to relocate the earthquakes that have occurred along the central and northern segments of the Longquanshan fault zone since 2009, including the MS 5.1 Qingbaijiang earthquake sequence, to investigate the seismogenic structure of the region. The relocation results indicate that the seismicity along the central and northern segments of the Longquanshan fault zone has occurred mainly along the eastern branch since 2009, with the hypocentral distribution along a vertical cross-section illustrating a steep, NW-dipping parallel imbricate structure. The terminating depth of the eastern branch is about 12 km. The distribution of the MS 5.1 Qingbaijiang earthquake sequence is along the NE–SW-striking Longquanshan fault zone. The aftershock focal depths are in the 3–6 km range, with the mainshock located at 104.475°E, 30.73°N. Its initial rupture depth of 5.2 km indicates that the earthquake occurred above the shallow decollement layer of the upper crust in this region. The hypocentral distribution along the long axis of the aftershock area highlights that this earthquake sequence occurred along a fault dipping at 56° to the NW. Our surface projection of the inferred fault plane places it near the eastern branch of the Longquanshan fault zone. We infer the MS 5.1 mainshock to be a thrust faulting event based on the focal mechanism solution via the cut-and-paste waveform inversion method, with strike/dip/rake parameters of 22°/36°/91° and 200°/54°/89° obtained for nodal planes I and II, respectively. We identify that the seismogenic fault of the MS 5.1 Qingbaijiang earthquake lies along the eastern branch of the Longquanshan fault zone, and nodal plane II represents the coseismic rupture plane, based on a joint analysis of the event relocation results, mainshock focal mechanism, and regional geological information. Our study provides vital information for assessing the seismic hazard of the Longquanshan fault zone near Chengdu City.


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.


Sign in / Sign up

Export Citation Format

Share Document