Seismic Monitoring for High-Precision Delineation of Fault Geometry and Stress; Case Study for the West Texas Subscription Array

Author(s):  
Sepideh Karimi ◽  
Adam Baig ◽  
Aaron Booterbaugh ◽  
Yoones Vaezi ◽  
Mark Stacey ◽  
...  

<p>Seismicity potentially induced through wastewater disposal, hydraulic fracture completion, or other industrial operations, has been a cause for increasing public concern over the last decade. Monitoring for this activity has focussed on the problems of location and characterization, often to a relatively rough degree of precision. Regulations typically spell out responses for operators should an event exceed a magnitude threshold within a specified distance of their facilities. While this type of monitoring is critical for ensuring injections be conducted effectively while minimizing potential damage from shaking and public alarm, it often leaves many unanswered questions in terms of the underlying processes.</p><p>Understanding these questions entails that we demand more out of the seismic networks, essentially upgrading the data products to a “next generation” level.  The data from the network needs to be used to provide a detailed understanding of critical geological structures and geomechanics of the study area. This goal is facilitated through both a densification of hardware and a higher order of event processing. High-precision locations delivered through relative relocation methodologies delineate slipping fault structures, often resolving previously unknown features. Moment tensor inversion processing also helps reveal the orientations of faults and provides information on stress in the region. The resolution of these structures provides critical insight into understanding how a field is reacting.</p><p>We illustrate the application of this “next-generation” seismicity monitoring system to the Delaware Basin in West Texas, where we have deployed a network of 25 broadband seismometers complementing monitoring from TexNet and other networks. Despite being an exceptionally challenging recording environment, by aggregating all of these data we obtain a high-resolution catalog of earthquake hypocenters delineating a number of fault features. Inverting the stresses from the moment tensors of the highest-quality events shows a dominantly normal stress regime and tangibly resolves a rotation of axes transitioning across the basin. We illustrate both the logistical and processing requirements necessary for timely delivery of results highlighting the dynamics of seismicity in an active study area.</p>

2019 ◽  
Vol 219 (2) ◽  
pp. 1447-1462 ◽  
Author(s):  
Alexandre P Plourde ◽  
Michael G Bostock

SUMMARY We introduce a new relative moment tensor (MT) inversion method for clusters of nearby earthquakes. The method extends previous work by introducing constraints from S-waves that do not require modal decomposition and by employing principal component analysis to produce robust estimates of excitation. At each receiver, P and S waves from each event are independently aligned and decomposed into principal components. P-wave constraints on MTs are obtained from a ratio of coefficients corresponding to the first principal component, equivalent to a relative amplitude. For S waves we produce constraints on MTs involving three events, where one event is described as a linear combination of the other two, and coefficients are derived from the first two principal components. Nonlinear optimization is applied to efficiently find best-fitting tensile-earthquake and double-couple solutions for relative MT systems. Using synthetic data, we demonstrate the effectiveness of the P and S constraints both individually and in combination. We then apply the relative MT inversion to a set of 16 earthquakes from southern Alaska, at ∼125 km depth within the subducted Yakutat terrane. Most events are compatible with a stress tensor dominated by downdip tension, however, we observe several pairs of earthquakes with nearly antiparallel slip implying that the stress regime is heterogeneous and/or faults are extremely weak. The location of these events near the abrupt downdip termination of seismicity and the low-velocity zone suggest that they are caused by weakening via grain-size and volume reduction associated with eclogitization of the lower crustal gabbro layer.


2021 ◽  
Author(s):  
Savvaidis Alexandros ◽  
Roselli Pamela

<p>In the scope to investigate the possible interactions between injected fluids, subsurface geology, stress field and triggering earthquakes, we investigate seismic source parameters related to the seismicity in West Texas (USA). The analysis of seismic moment tensor is an excellent tool to understand earthquake source process kinematics; moreover, changes in the fluid volume during faulting leads to existence of non-double-couple (NDC) components (Frohlich, 1994; Julian et al., 1998; Miller et al., 1998). The NDC percentage in the source constitutes the sum of absolute ISO and CLVD components so that %NDC= % ISO + %CLVD and %ISO+%CLVD+%DC=100%. It is currently known that the presence of NDC implies more complex sources (mixed shear-tensile earthquakes) correlated to fluid injections, geothermal systems and volcano-seismology where induced and triggered seismicity is observed.</p><p>With this hypothesis, we analyze the micro-earthquakes (M <2 .7) recorded by the Texas Seismological Network (TexNet) and a temporary network constituted by 40 seismic stations (equipped by either broadband or 3 component geophones). Our study area is characterized by Northwest-Southeast faults that follow the local stress/field (SH<sub>max</sub>) and the geological characteristic of the shallow basin structure of the study area. After a selection based on signal-to-noise ratio, we filter (1-50 Hz) the seismograms and estimate P-wave pulse polarities and the first P-wave ground displacement pulse in time domain. Then, we perform the full moment tensor analysis by using hybridMT technique (Andersen, 2001; Kwiatek et al., 2016) with a detailed 1D velocity model. The key parameter is the polarity/area of the first P-wave ground displacement pulse in time domain. Uncertainties of estimated moment tensors are expressed by normalized root-mean-square (RMS errors) between theoretical and estimated amplitudes (Vavricuk et al., 2014). We also evaluate the quality of the seismic moment tensors by bootstrap and resampling. In our preliminary results we obtain NDC percentage (in terms of %ISO and %CLVD components), Mw, seismic moment, P, T and B axes orientation for each source inverted.</p>


2020 ◽  
Vol 91 (6) ◽  
pp. 3550-3562
Author(s):  
Qipeng Bai ◽  
Sidao Ni ◽  
Risheng Chu ◽  
Zhe Jia

Abstract Earthquake moment tensors and focal depths are crucial to assessing seismic hazards and studying active tectonic and volcanic processes. Although less powerful than strong earthquakes (M 7+), moderately strong earthquakes (M 5–6.5) occur more frequently and extensively, which can cause severe damages in populated areas. The inversion of moment tensors is usually affected by insufficient local waveform data (epicentral distance <5°) in sparse seismic networks. It would be necessary to combine local and teleseismic data (epicentral distance 30°–90°) for a joint inversion. In this study, we present the generalized cut-and-paste joint (gCAPjoint) algorithm to jointly invert full moment tensor and centroid depth with local and teleseismic broadband waveforms. To demonstrate the effectiveness and explore the limitations of this algorithm, we perform case studies on three earthquakes with different tectonic settings and source properties. Comparison of our results with global centroid moment tensor and other catalog solutions illustrates that both non-double-couple compositions of the focal mechanisms and centroid depths can be reliably recovered for very shallow (<10  km) earthquakes and intermediate-depth events with this software package.


2020 ◽  
Vol 110 (5) ◽  
pp. 2225-2241 ◽  
Author(s):  
Alexandros Savvaidis ◽  
Anthony Lomax ◽  
Caroline Breton

ABSTRACT Most current seismicity in the southern U.S. midcontinent is related to oil and gas operations (O&G Ops). In Texas, although recorded earthquakes are of low-to-moderate magnitude, the rate of seismicity has been increasing since 2009. Because of the newly developed Texas Seismological Network, in most parts of Texas, recent seismicity is reported on a daily basis with a magnitude of completeness of ML 1.5. Also, funded research has allowed the collection of O&G Op information that can be associated with seismicity. Although in the Dallas–Fort Worth area, recent seismicity has been associated mostly with saltwater disposal (SWD), in the South Delaware Basin, West Texas, both hydraulic fracturing (HF) and SWD have been found to be causal factors. We have begun to establish an O&G Op database using four different sources—IHS, FracFocus, B3, and the Railroad Commission of Texas—with which we can associate recent seismicity to HF and SWD. Our approach is based on time and epicentral location of seismic events and time, location of HF, and SWD. Most seismicity occurs in areas of dense HF and SWD-well activity overlapping in time, making association of seismicity with a specific well type impossible. However, through examination of clustered seismicity in space and time, along with isolated clusters of spatiotemporal association between seismicity and O&G Ops, we are able to show that a causation between HF and seismicity may be favored over causation with SWD wells in areas of spatially isolated earthquake clusters (Toyah South, Reeves West, Jeff Davis Northeast, and Jeff Davis East). Causality between SWD and seismicity may be inferred for isolated cases in Reeves South and Grisham West.


1996 ◽  
Vol 86 (5) ◽  
pp. 1255-1269 ◽  
Author(s):  
Michael E. Pasyanos ◽  
Douglas S. Dreger ◽  
Barbara Romanowicz

Abstract Recent advances in broadband station coverage, continuous telemetry systems, moment-tensor procedures, and computer data-processing methods have given us the opportunity to automate the two regional moment-tensor methods employed at the UC Berkeley Seismographic Station for events in northern and central California. Preliminary solutions are available within minutes after an event has occurred and are subsequently human reviewed. We compare the solutions of the two methods to each other, as well as the automatic and revised solutions of each individual method. Efforts are being made to establish robust criteria for determining accurate solutions with human review and to fully automate the moment-tensor procedures into the already-existing automated earthquake-location system.


1983 ◽  
Vol 73 (2) ◽  
pp. 419-434
Author(s):  
Jeffery S. Barker ◽  
Charles A. Langston

abstract Teleseismic P-wave first motions for the M ≧ 6 earthquakes near Mammoth Lakes, California, are inconsistent with the vertical strike-slip mechanisms determined from local and regional P-wave first motions. Combining these data sets allows three possible mechanisms: a north-striking, east-dipping strike-slip fault; a NE-striking oblique fault; and a NNW-striking normal fault. Inversion of long-period teleseismic P and SH waves for the events of 25 May 1980 (1633 UTC) and 27 May 1980 (1450 UTC) yields moment tensors with large non-double-couple components. The moment tensor for the first event may be decomposed into a major double couple with strike = 18°, dip = 61°, and rake = −15°, and a minor double couple with strike = 303°, dip = 43°, and rake = 224°. A similar decomposition for the last event yields strike = 25°, dip = 65°, rake = −6°, and strike = 312°, dip = 37°, and rake = 232°. Although the inversions were performed on only a few teleseismic body waves, the radiation patterns of the moment tensors are consistent with most of the P-wave first motion polarities at local, regional, and teleseismic distances. The stress axes inferred from the moment tensors are consistent with N65°E extension determined by geodetic measurements by Savage et al. (1981). Seismic moments computed from the moment tensors are 1.87 × 1025 dyne-cm for the 25 May 1980 (1633 UTC) event and 1.03 × 1025 dyne-cm for the 27 May 1980 (1450 UTC) event. The non-double-couple aspect of the moment tensors and the inability to obtain a convergent solution for the 25 May 1980 (1944 UTC) event may indicate that the assumptions of a point source and plane-layered structure implicit in the moment tensor inversion are not entirely valid for the Mammoth Lakes earthquakes.


Sign in / Sign up

Export Citation Format

Share Document