SEARCH · Engineering Papers
Results for “fault location”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves
Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.
End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods
Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.
Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation
At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.
Graph-Learning-Assisted State and Event Tracking for Solar-Penetrated Power Grids with Heterogeneous Data Sources
Unlike transmission systems, distribution systems do not typically contain sufficient metering to enable real-time state estimation. The lack of sufficient real-time measurements prohibits accurate and timely monitoring of the state of distribution systems. As a result, control and optimal operation of distribution systems, especially those containing large numbers of renewable generation units are not possible without proper data and information about the current state of the system. The main motivation of this project is to address this shortcoming by developing an approach which provides “predicted” real-time measurements so that they can be used to execute a distribution system state estimator. Thus, the objective of the project is to make the distribution systems fully observable, such that the hosting capacity for solar generation can be accurately estimated, and unnecessary solar curtailments can be avoided. In order to accomplish this goal, the project investigated the use of a grid-model-informed machine learning (ML) tool which integrates heterogeneous data streams obtained from AMI meters, SCADA as well as PMU measurements and created synchronous measurement snapshots for the state estimator (SE); and developed a hybrid robust SE which provides not only accurate state estimates but also real-time feedback for the ML model refinement.
Towards Automated/Semiautomated Extraction of Faults from Lidar Data
The Pajarito fault system is a complex zone of deformation and a seismically active region nestled within the Rio Grande rift in north-central New Mexico. Numerous laterally discontinuous faults and associated folds and fractures interact in a manner that has important implications for seismic hazards and risk mitigation. Previous efforts have established a foundation for the location of lineaments and structures in the Pajarito fault system; however, ensuring the completeness of the current lineament mapping is required for identifying areas for field validation, evaluating the potential for future seismic activity, and better understanding fault interaction. Assistance with this fault-mapping task via automated or semiautomated techniques as applied to lidar data over a large area of interest is highly desirable. A proof-of-concept processing flow which transforms lidar point-cloud data into a raster of surficial fault candidates is described and illustrated herein. We report these initial results hold great promise toward achieving our ultimate goal.
In-situ experiment reveals CO2 enriched fluid migration in faulted caprock
Abstract The sealing characteristics of the geological formation located above a CO 2 storage reservoir, the so-called caprock, are essential to ensure efficient geological carbon storage. If CO 2 were to leak through the caprock, temporal changes in fluid geochemistry can reveal fundamental information on migration mechanisms and induced fluid–rock interactions. Here, we present the results from a unique in-situ injection experiment, where CO 2 -enriched fluid was continuously injected in a faulted caprock analogue. Our results show that the CO 2 migration follows complex pathways within the fault structure. The joint analysis of noble gases, ion concentrations and carbon isotopes allow us to quantify mixing between injected CO 2 -enriched fluid and resident formation water and to describe the temporal evolution of water–rock interaction processes. The results presented here are a crucial complement to the geophysical monitoring at the fracture scale highlighting a unique migration of CO 2 in fault zones.
Fault MLReal: a fault delineation study for the Decatur $CO_{2}$ field data using neural network predicted passive seismic locations
Explore the source record for details and available documents.
Precursory Locking Precedes Slip Events on Laboratory Fault
Earthquakes occur through repeated cycles of locking and slip. What is the relationship between transiently locked patches and the size, location, and timing of the subsequent slip events? Here we use an experimental fault made of transparent rubber to image the slip history of over 18,000 events. The soft elastic system allows for contained ruptures and thus a natural evolution of stress heterogeneity over multiple earthquake cycles. We observe a precursory locking phase preceding slip nucleation for events of all sizes. The integrated locking distance increases linearly with the moment of the eventual event. The most intense (largest slip deficit) locking events are spatially and temporally correlated with the hypocenters of the largest rupture events. These results show that the transient locking is statistically, but not deterministically, predictive of the location, timing, and size of the eventual slip event.
IoT-Enabled Traveling Wave Microgrid Protection
Traveling Wave Protection based on the Internet of Things (TWP-IoT) is developed to enable ultra-resilient microgrids. TWP-IoT adopts a directional traveling wave approach that utilizes the voltage and current modal components obtained from Discrete Hilbert Transform to identify the fault features via ultra-lightweight IoT hardware. Furthermore, the main contributions of this work include (1) a novel TWP-IoT for microgrid and distribution networks that can achieve the time-of-arrival fault detection and localization with composite wave impedance and conductance using Discrete Hilbert Transform, (2) a cost-effective smart IoT prototype with functionalities outperforming commercial relays, and (3) a real-time hardware-in-the-loop prototype based on RTDS that verifies the robustness and efficacy of TWP-IoT subject to a diverse set of working conditions. TWP-IoT is found to have excellent performance under a wide spectrum of internal and external fault conditions across a myriad of microgrid locations, which is unattainable by today’s TWP products.
How Good Is Your Location? Comparing and Understanding the Uncertainties in Location for the 1993 Rock Valley Sequence
Accurate event locations are important for many endeavors in seismology, and understanding the factors that contribute to uncertainties in those locations is complex. In this article, we present a case study that takes an in-depth look at the accuracy and precision possible for locating nine shallow earthquakes in the Rock Valley fault zone in southern Nevada. These events are targeted by the Rock Valley Direct Comparison phase of the Source Physics Experiment, as candidates for the colocation of a chemical explosion with an earthquake hypocenter to directly compare earthquake and explosion sources. For this comparison, it is necessary to determine earthquake hypocenters as accurately as possible so that different source types have nearly identical locations. Our investigations include uncertainty analysis from different sets of phase arrivals, stations, velocity models, and location algorithms. For a common set of phase arrivals and stations, we find that epicentral locations from different combinations of velocity models and algorithms are within 600 m of one another in most cases. Event depths exhibit greater uncertainties, but focusing on the S-P times at the nearest station allows for estimates within approximately 500 m.
Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective
Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.
End-to-end microgrid protection using distributed data-driven methods
This paper introduces an end-to-end microgrid protection framework that offers real-time system monitoring, fault-related decision making, and circuit breaker control. This is achieved through the design of distributed data-driven techniques based on the support vector machine method, where each relay is responsible for distributed data collection, fault detection, fault localization, and fault isolation. Local communication is established among neighboring relays, fostering cooperative fault localization and isolation. This decentralized design not only reduces the computational and communication requirements but also enables the adaptability of each relay under varying operational dynamics. The proposed end-to-end protection framework was validated using MATLAB/Simulink simulations on a 100% renewable microgrid, achieving an accuracy of 93.1% with response time of 0.0523 s, in protecting against a range of fault scenarios that are characterized by various types, locations, impedances, load conditions, photovoltaic power levels, and microgrid operating modes.
Temporal and spatial characterization of a thermogenic, fault-controlled gas hydrate system, Woolsey Mound, Gulf of Mexico
Woolsey Mound, located at Mississippi Canyon Lease Block 118 (MC118), is the site of the Gulf of Mexico hydrate research consortium’s seafloor observatory, where gas hydrates outcrop at the seafloor. The presence of gas hydrates in the mound is confirmed directly by coring and indirectly by 3D seismic reflection data. Craters, pockmarks, chemosynthetic communities, and authigenic carbonates populate the seafloor at Woolsey Mound. Each crater is characterized by a network of shallow crestal faults that connect the hydrate mound to the underlying allochthonous salt body. We characterize the temporal and spatial evolution of gas hydrates at Woolsey Mound under natural perturbations using four collocated 3D seismic reflection data sets that span over 14 years. Data acquisition differences embedded in the data sets arising from variation in geometry, sample rate, and phase are minimized using the “cross-equalization” method. Our results indicate that hydrate formation and dissociation vary temporally and spatially in close connection to the shallow crestal faults. Evidence of gas hydrate dissociation is observed over a period of three years (2000–2003), where major dissociation occurred along the southern portion of the crestal fault in the southeast crater. The dissociation is less prominent in the southwest crater. Evidence of methane venting is observed between 2000 and 2010, which is mostly concentrated in the southeast crater. The residual amplitude anomalies observed between 2000 and 2014 in the mound are mostly positive, implying that the methane venting had increased significantly. The positive anomalies are correlated with the methane seepage recorded in 2011. Our results indicate the evolution of a fault-controlled gas hydrate system in the northern Gulf of Mexico, which would aid in assessing its impact on the seafloor.
Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations
SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.
Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks
Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.
Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers
A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.
Performance assessment of near-fault buildings subjected to physics-based simulated earthquake ground motions with fling step
The effects of the co-seismic static offset (known as fling step) and associated velocity pulses on civil structures have been difficult to study because the static offset is typically removed during the processing of earthquake ground motion records. Simulated ground motions contain fling features and require no processing; therefore, they create new opportunities for representing fling features in seismic hazard analysis and assessing their influence on the seismic demands on near-fault structures. We use physics-based fault rupture simulations to study the characteristics of ground motions with fling step and the sensitivity of the near-fault structural demands to strong fling features. We uncover that simulated ground motions with a large fling step tend to have higher spectral intensity than those without a fling step at the same rupture distance, especially at periods longer than 2 s. As a result, the structural demands on flexible buildings tend to be the most sensitive to the fling features. Statistical analysis suggests that the ground motion spectral shape (represented by spectral accelerations at multiple periods) is—in most cases—a sufficient predictor of the structural demands on near-fault low-rise and mid-rise buildings at locations that are susceptible to strong fling effects. Finally, ground motion record selection experiments reveal that representing the spectral shape features at periods that are most relevant to a given structure may be an effective strategy to reduce the bias in the estimated demands on near-fault long-period structures when the available database of records is considered deficient in fling features.