48 research outputs found
Recommended from our members
Waveform-relocated earthquake catalog for Oklahoma and southern Kansas illuminates the regional fault network
For much of Oklahoma, augmentation of the seismic network with new public stations in the activated areas has followed rather than preceded the spread of seismicity across the state, and consequently the network geometry is often unfavorable for resolving the underlying fault structures. With this study, we reanalyze the existing earthquake catalog with additional data from two industry-operated networks for the period May 2013 to November 2016. These networks include 40 seismic stations and cover seismically active north-central Oklahoma with a station spacing on the order of 25 km. Relative locations obtained from waveform cross correlation reveal a striking pattern of seismicity, illuminating many previously unmapped faults. Absolute depths are usually well constrained to within 1 km. Relative locations provide about one order of magnitude better precision for resolving the structure of seismicity clusters. Relocated epicenters tend to cluster in linear trends of less than 1 km to more than 20 km in length. In areas with stations closer than about 10 km, we can resolve fault planes by strike and dip. These are generally in agreement with surface-wave-derived moment-Tensor solutions
Recommended from our members
A Systematic Assessment of the Spatiotemporal Evolution of Fault Activation Through Induced Seismicity in Oklahoma and Southern Kansas
Much of Oklahoma and southern Kansas has seen widespread seismic activity in the last decade that is attributed to large-scale wastewater disposal into the Arbuckle group. Using a waveform-relocated earthquake catalog, we perform a systematic study of the activity on several hundreds of identified faults. We use 93 sequences with at least 30 events for a detailed analysis of their spatiotemporal evolution. For most awakened faults, seismicity tends to initiate at shallower depth and migrates deeper along the faults as the sequence proceeds. No major sequence starts with the largest earthquake, and many sequences initiate months before they rise to peak activity. We study temporal clustering as a means to quantify earthquake interactions. Some sequences show no temporal clustering similar to Poissonian background seismicity but at much higher rate than the natural background. Other sequences exhibit strong temporal clustering akin to main shock-aftershock sequences. We conclude that once initiated by anthropogenic forcing, portions of the activated faults in the Oklahoma/Kansas area are close enough to failure to continue failing through earthquake interactions. In many sequences, including those with the largest earthquakes, seismicity continues within the previously activated region rather than by growing the activated area. Therefore, monitoring seismicity with a low magnitude threshold and high location precision has the potential to detect minor activity as it initiates failure on specific faults and thus provide time to take actions to mitigate the occurrence of potentially damaging earthquakes
Recommended from our members
Waveform-relocated earthquake catalog for Oklahoma and southern Kansas illuminates the regional fault network
For much of Oklahoma, augmentation of the seismic network with new public stations in the activated areas has followed rather than preceded the spread of seismicity across the state, and consequently the network geometry is often unfavorable for resolving the underlying fault structures. With this study, we reanalyze the existing earthquake catalog with additional data from two industry-operated networks for the period May 2013 to November 2016. These networks include 40 seismic stations and cover seismically active north-central Oklahoma with a station spacing on the order of 25 km. Relative locations obtained from waveform cross correlation reveal a striking pattern of seismicity, illuminating many previously unmapped faults. Absolute depths are usually well constrained to within 1 km. Relative locations provide about one order of magnitude better precision for resolving the structure of seismicity clusters. Relocated epicenters tend to cluster in linear trends of less than 1 km to more than 20 km in length. In areas with stations closer than about 10 km, we can resolve fault planes by strike and dip. These are generally in agreement with surface-wave-derived moment-Tensor solutions
Recommended from our members
Quantifying the heterogeneity of the tectonic stress field using borehole data
The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with subvertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures, and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km
Recommended from our members
Quantifying the heterogeneity of the tectonic stress field using borehole data
The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with subvertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures, and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km
Recommended from our members
How faults wake up: The Guthrie-Langston, Oklahoma earthquakes
Large-scale wastewater disposal has led to a fast-paced reawakening of faults in the Oklahoma/Kansas region. High-resolution earthquake relocations show that the inventory of ancient basement faults in the study region differs from results of seismic surveys and geologic mapping focused on the sedimentary cover. We analyze the evolution of seismic activity in the Guthrie-Langston sequence in central Oklahoma in greater detail. Here, seismic activity has reactivated a network of at least 12 subvertical faults in an area less than 10 km across. Recorded activity began in late 2013, peaked about six months later, and includes two M 4 earthquakes. These earthquakes characteristically occur at about 4 km depth below the top of the basement and do not reach the sedimentary cover. The sequence shows a radial growth pattern despite being no closer than 10 km to significant wastewater disposal activity. Hydrologic modeling suggests that activity initiated with a time lag of several years relative to early injection activity. Once initiated, earthquake interactions contribute to the propagation of seismicity along the reactivated faults. As a result, the spatiotemporal evolution of the seismicity mimics a diffusive pattern that is typically thought to be associated with injection activity. Analysis of the fault slip potential shows that most faults are critically stressed in the contemporary stress field. Activity on some faults, for which we find low slip probability, suggests a significant contribution of geomechanical heterogeneities to the reawakening of these ancient basement faults
How faults wake up: The Guthrie-Langston, Oklahoma earthquakes
Large-scale wastewater disposal has led to a fast-paced reawakening of faults in the Oklahoma/Kansas region. High-resolution earthquake relocations show that the inventory of ancient basement faults in the study region differs from results of seismic surveys and geologic mapping focused on the sedimentary cover. We analyze the evolution of seismic activity in the Guthrie-Langston sequence in central Oklahoma in greater detail. Here, seismic activity has reactivated a network of at least 12 subvertical faults in an area less than 10 km across. Recorded activity began in late 2013, peaked about six months later, and includes two M 4 earthquakes. These earthquakes characteristically occur at about 4 km depth below the top of the basement and do not reach the sedimentary cover. The sequence shows a radial growth pattern despite being no closer than 10 km to significant wastewater disposal activity. Hydrologic modeling suggests that activity initiated with a time lag of several years relative to early injection activity. Once initiated, earthquake interactions contribute to the propagation of seismicity along the reactivated faults. As a result, the spatiotemporal evolution of the seismicity mimics a diffusive pattern that is typically thought to be associated with injection activity. Analysis of the fault slip potential shows that most faults are critically stressed in the contemporary stress field. Activity on some faults, for which we find low slip probability, suggests a significant contribution of geomechanical heterogeneities to the reawakening of these ancient basement faults
Recommended from our members
Analysis of borehole breakout development using continuum damage mechanics
Damage distribution and evolution have a significant effect on borehole stress concentrations. To model the complex fracturing process and inelastic deformation in the development of the borehole breakout, we implement a continuum damage mechanics (CDM) concept that takes tensile and compressive failure mechanisms into account. The proposed approach explicitly models the dissipative behavior of the material due to cracking and its evolution, which leads to an inhomogeneous redistribution of material properties and stresses in the vicinity of the borehole wall. We apply a constitutive plastic model for Berea sandstone and compare our numerical results to laboratory experiments performed on Tablerock sandstone. We are able to reproduce several characteristics of the failure process during the breakout development as observed in experimental tests, e.g. localized crack distribution in the vicinity of the borehole wall, damage evolution, which exhibits a widening process in the beginning followed by subsequent growth in depth, and shear fracturing-dominated breakout growth in sandstone. A comparison of our results with laboratory experiments performed on a range of stress conditions shows a good agreement of the size of borehole breakouts. The importance of the constitutive damage law in defining the failure mechanisms of the damaging processes is discussed. We show that the depth and the width of breakouts are not independent of each other and no single linear relation can be found between the size of breakouts and the magnitude of the applied stress. Consequently, only one far field principal stress component can be estimated from breakout geometry, if the other two principal stresses are known and sufficient data on the plastic parameters are available
Recommended from our members
Coupled Hydromechanical Modeling of Induced Seismicity From CO2 Injection in the Illinois Basin
Injection of CO2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate