Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMICITY”

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.

At least 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Frequency-Domain-Based Algorithm for Detecting Microseismicity Using Dense Surface Seismic Arrays

We propose a new frequency-domain-based algorithm for detecting small-magnitude seismic events using dense surface seismic arrays. Our proposed method takes advantage of the high energy carried by S waves, and approximate known source locations, which are used to rotate the horizontal components to obtain the maximum amplitude. By surrounding the known source area with surface geophones, we achieve a favorable geometry for locating the detected seismic events with the backprojection method. To test our new detection method, we used a dense circular array, consisting of 151 5 Hz three-component geophones, over a 5 km aperture that was in operation at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) in southcentral Utah. We apply the new detection method during a small-scale test injection phase at FORGE, and during an aftershock sequence of an Mw 4.1 earthquake located ~30 km north of the geophone array, within the Black Rock volcanic field. We are able to detect and locate microseismic events (Mw<0) during injections, despite the high level of anthropogenic activity, and several aftershocks that are missing from the regional catalog. By comparing our method with known algorithms that operate both in the time and frequency domain, we show that our proposed method performs better in the case of the FORGE injection monitoring, and equally well for the off-array aftershock sequence. Our new method has the potential to improve microseismic event detections even in extremely noisy environments, and the proposed location scheme serves as a direct discriminant between true and false detections.

Geochemistry & Geophysics↗
DOE OSTI2024

Seismic Tomography 2024

Seismic tomography is the most abundant source of information about the internal structure of the Earth at scales ranging from a few meters to thousands of kilometers. It constrains the properties of active volcanoes, earthquake fault zones, deep reservoirs and storage sites, glaciers and ice sheets, or the entire globe. It contributes to outstanding societal problems related to natural hazards, resource exploration, underground storage, and many more. The recent advances in seismic tomography are being translated to nondestructive testing, medical ultrasound, and helioseismology. Nearly 50 yr after its first successful applications, this article offers a snapshot of modern seismic tomography. Focused on major challenges and particularly promising research directions, it is intended to guide both Earth science professionals and early-career scientists. The individual contributions by the coauthors provide diverse perspectives on topics that may at first seem disconnected but are closely tied together by a few coherent threads: multiparameter inversion for properties related to dynamic processes, data quality, and geographic coverage, uncertainty quantification that is useful for geologic interpretation, new formulations of tomographic inverse problems that address concrete geologic questions more directly, and the presentation and quantitative comparison of tomographic models. Finally, it remains to be seen which of these problems will be considered solved, solved to some extent, or practically unsolvable over the next decade.

58 GEOSCIENCES↗
DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗
DOE OSTI2017

Identifying and Correcting Timing Errors at Seismic Stations in and around Iran

A fundamental component of seismic research is the use of phase arrival times, which are central to event location, Earth model development, and phase identification, as well as derived products. Hence, the accuracy of arrival times is crucial. However, errors in the timing of seismic waveforms and the arrival times based on them may go unidentified by the end user, particularly when seismic data are shared between different organizations. Here, we present a method used to analyze travel-time residuals for stations in and around Iran to identify time periods that are likely to contain station timing problems. For the 14 stations with the strongest evidence of timing errors lasting one month or longer, timing corrections are proposed to address the problematic time periods. Finally, two additional stations are identified with incorrect locations in the International Registry of Seismograph Stations, and one is found to have erroneously reported arrival times in 2011.

58 GEOSCIENCES↗
DOE OSTI2018

Interactive Visualization of Complex Seismic Data and Models Using Bokeh

Visualizing multidimensional data and models becomes more challenging as the volume and resolution of seismic data and models increase. But thanks to the development of powerful and accessible computer systems, a model web browser can be used to visualize complex scientific data and models dynamically. In this paper, we present four examples of seismic model visualization using an open-source Python package Bokeh. One example is a visualization of a surface-wave dispersion data set, another presents a view of three-component seismograms, and two illustrate methods to explore a 3D seismic-velocity model. Unlike other 3D visualization packages, our visualization approach has a minimum requirement on users and is relatively easy to develop, provided you have reasonable programming skills. Finally, utilizing familiar web browsing interfaces, the dynamic tools provide us an effective and efficient approach to explore large data sets and models.

58 GEOSCIENCES↗
DOE OSTI2021

A Data-Driven Framework for Automated Detection of Aircraft-Generated Signals in Seismic Array Data Using Machine Learning

Abstract Ground motions associated with aircraft overflights can cover a significant portion of the seismic data collected by shallowly emplaced seismometers, such as new nodal and Distributed Acoustic Sensing systems. This article describes the first published framework for automated detection of aircraft on single channel and multichannel seismic data. The seismic data are converted to spectrograms in a sliding time window and classified as aircraft or nonaircraft in each window using a deep convolutional neural network trained with analyst-labeled data. A majority voting scheme is used to convert the output from the sequence of sliding time windows onto a decision time sequence for each channel and to combine the binary classifications on the decision time sequences across multiple channels. Precision, recall, and F-score are used to quantify the detection performance of the algorithm on nodal data using fourfold time-series cross validation. By applying our framework to data from the Sage Brush Flats nodal array in Southern California, we provide a benchmark performance and demonstrate the advantage of using an array of sensors.

Geochemistry & Geophysics↗
DOE OSTI2022

The Imperial Valley Dark Fiber Project: Toward Seismic Studies Using DAS and Telecom Infrastructure for Geothermal Applications

We report that the Imperial Valley is a seismically active basin occupying the southern end of the Salton trough, an area of rapid extension, high heat flow, and abundant geothermal resources. This report describes an ongoing large-scale distributed acoustic sensing (DAS) recording study acquiring high-density seismic data on an array between Calipatria and Imperial, California. This 27 km array, operating on dark fiber since 9 November 2020, has recorded a wealth of local seismic events as well as ambient noise. The goal of the broader Imperial Valley Dark Fiber project is to evaluate passive DAS as a tool for geothermal exploration and monitoring. This report is intended to provide installation information, noise characteristics, and metadata for future studies utilizing the data set. Because of the relatively small number of basin-scale DAS studies that have been conducted to date, we also provide a range of lessons learned during the deployment to assist future researchers exploring this acquisition strategy.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Tracking Seismic Velocity Perturbations at Ridgecrest Using Ballistic Correlation Functions

Abstract We present results based on data of a dense nodal array composed of 147 stations, deployed in 2022 near the epicenter of the 2019 Mw 7.1 Ridgecrest earthquake to investigate characteristics of the seismic wavefields. Through array analyses, we identified two primary components. First, we observed far-field P waves dominating the 0.5–1.2 Hz frequency range, which are likely primarily generated by wind-driven oceanic swell activity. Second, we detected near-field body waves resulting from anthropogenic activities in the frequency range 2–8 Hz. We examined noise correlation functions derived from data of the dense deployment and regional stations to explore fault-zone seismic velocity changes using ballistic arrivals, with a focus on velocity perturbation shortly before and after the Ridgecrest earthquake sequence. Our findings exhibit distinct behavior compared to results obtained through standard coda-wave interferometry. Particularly, we observed a decrease in P-wave travel time on certain station pairs prior to the 2019 earthquake sequence. Supported by detailed investigation of the local seismic wavefields, we interpret the decreasing P-wave travel time as likely caused by a velocity increase away from the fault, possibly related to fluid migration. However, additional information is necessary to verify this hypothesis.

Geochemistry & Geophysics↗
DOE OSTI2025

Evaluation of Seismic Artificial Intelligence with Uncertainty

Artificial intelligence has transformed the seismic community with deep learning models (DLMs) that are trained to complete specific tasks within workflows. However, there is still a lack of robust evaluation frameworks for evaluating and comparing DLMs. Here, we address this gap by designing an evaluation framework that jointly incorporates two crucial aspects: performance uncertainty and learning efficiency. To target these aspects, we meticulously construct the training, validation, and test splits using a clustering method tailored to seismic data and enact an expansive training design to segregate performance uncertainty arising from stochastic training processes and random data sampling. The framework’s ability to guard against misleading declarations of model superiority is demonstrated through the evaluation of PhaseNet (Zhu and Beroza, 2018), a popular seismic phase picking DLM, under three training approaches. Our framework helps practitioners choose the best model for their problem and set performance expectations by explicitly analyzing model performance with uncertainty at varying budgets of training data.

58 GEOSCIENCES↗
DOE OSTI2011

Comparative in situ X-ray Diffraction Study of San Carlos Olivine: Influence of Water on the 410 km Seismic Velocity Jump in Earth’s Mantle

A comparative study of the equation of states of hydrous (0.4 wt% H{sub 2}O) and anhydrous San Carlos olivine (<30 ppm H2O) was conducted using synchrotron X-rays up to 11 GPa in a diamond anvil cell (DAC) at ambient temperature. Both samples were loaded in the same high-pressure chamber of the DAC to eliminate the possible pressure difference in different experiments. The obtained compression data were fitted to the third-order Birch-Murnaghan equation of state, yielding a bulk modulus K{sub 0} = 123(3) GPa for hydrous olivine and K{sub 0} = 130(4) GPa for anhydrous olivine as K{sub 0}' is fixed at 4.6. Therefore, 0.4 wt% H{sub 2}2O in olivine results in a 5% reduction in bulk modulus. Previous studies reported bulk modulus reduction by water in olivine's high-pressure polymorph (wadsleyite), to which the transformation from olivine gives rise to the seismic discontinuity at 410 km depth. The new data results in a reduction in the magnitude of the discontinuity by 50% in v{sub P} and 30% in v{sub S} (for 1:5 water partitioning between olivine and wadsleyite) with respect to anhydrous mantle. Previous knowledge of the influence of water on this phase transition has been in opposition to a large amount of water [e.g., 200 ppm by Wood (1995)] existing at 410 km depth. Calculation of the seismic velocities based on newly available elasticity data of the hydrous phases indicates that the presence of water is favorable for the mineral composition model (pyrolite) and seismic observations in terms of the magnitude of the 410 km discontinuity.

58 GEOSCIENCES↗
DOE OSTI2010

ADVANCED SEISMIC BASE ISOLATION METHODS FOR MODULAR REACTORS

Advanced technologies for structural design and construction have the potential for major impact not only on nuclear power plant construction time and cost, but also on the design process and on the safety, security and reliability of next generation of nuclear power plants. In future Generation IV (Gen IV) reactors, structural and seismic design should be much more closely integrated with the design of nuclear and industrial safety systems, physical security systems, and international safeguards systems. Overall reliability will be increased, through the use of replaceable and modular equipment, and through design to facilitate on-line monitoring, in-service inspection, maintenance, replacement, and decommissioning. Economics will also receive high design priority, through integrated engineering efforts to optimize building arrangements to minimize building heights and footprints. Finally, the licensing approach will be transformed by becoming increasingly performance based and technology neutral, using best-estimate simulation methods with uncertainty and margin quantification. In this context, two structural engineering technologies, seismic base isolation and modular steel-plate/concrete composite structural walls, are investigated. These technologies have major potential to (1) enable standardized reactor designs to be deployed across a wider range of sites, (2) reduce the impact of uncertainties related to site-specific seismic conditions, and (3) alleviate reactor equipment qualification requirements. For Gen IV reactors the potential for deliberate crashes of large aircraft must also be considered in design. This report concludes that base-isolated structures should be decoupled from the reactor external event exclusion system. As an example, a scoping analysis is performed for a rectangular, decoupled external event shell designed as a grillage. This report also reviews modular construction technology, particularly steel-plate/concrete construction using factory prefabricated structural modules, for application to external event shell and base isolated structures.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2010

OSI Passive Seismic Experiment at the Former Nevada Test Site

On-site inspection (OSI) is one of the four verification provisions of the Comprehensive Nuclear Test Ban Treaty (CTBT). Under the provisions of the CTBT, once the Treaty has entered into force, any signatory party can request an on-site inspection, which can then be carried out after approval (by majority voting) of the Executive Council. Once an OSI is approved, a team of 40 inspectors will be assembled to carry out an inspection to ''clarify whether a nuclear weapon test explosion or any other nuclear explosion has been carried out in violation of Article I''. One challenging aspect of carrying out an on-site inspection (OSI) in the case of a purported underground nuclear explosion is to detect and locate the underground effects of an explosion, which may include an explosion cavity, a zone of damaged rock, and/or a rubble zone associated with an underground collapsed cavity. The CTBT (Protocol, Section II part D, paragraph 69) prescribes several types of geophysical investigations that can be carried out for this purpose. One of the methods allowed by the CTBT for geophysical investigation is referred to in the Treaty Protocol as ''resonance seismometry''. This method, which was proposed and strongly promoted by Russia during the Treaty negotiations, is not described in the Treaty. Some clarification about the nature of the resonance method can be gained from OSI workshop presentations by Russian experts in the late 1990s. Our understanding is that resonance seismometry is a passive method that relies on seismic reverberations set up in an underground cavity by the passage of waves from regional and teleseismic sources. Only a few examples of the use of this method for detection of underground cavities have been presented, and those were done in cases where the existence and precise location of an underground cavity was known. As is the case with many of the geophysical methods allowed during an OSI under the Treaty, how resonance seismology really works and its effectiveness for OSI purposes has yet to be determined. For this experiment, we took a broad approach to the definition of ''resonance seismometry''; stretching it to include any means that employs passive seismic methods to infer the character of underground materials. In recent years there have been a number of advances in the use of correlation and noise analysis methods in seismology to obtain information about the subsurface. Our objective in this experiment was to use noise analysis and correlation analysis to evaluate these techniques for detecting and characterizing the underground damage zone from a nuclear explosion. The site that was chosen for the experiment was the Mackerel test in Area 4 of the former Nevada Test Site (now named the Nevada National Security Site, or NNSS). Mackerel was an underground nuclear test of less than 20 kT conducted in February of 1964 (DOENV-209-REV 15). The reason we chose this site is because there was a known apical cavity occurring at about 50 m depth above a rubble zone, and that the site had been investigated by the US Geological Survey with active seismic methods in 1965 (Watkins et al., 1967). Note that the time delay between detonation of the explosion (1964) and the time of the present survey (2010) is nearly 46 years - this would not be typical of an expected OSI under the CTBT.

58 GEOSCIENCES↗
DOE OSTI2011

Application of the Split Hopkinson Resonant Bar Test for Seismic Property Characterization of Hydrate-bearing Sand Undergoing Water Saturation

Conventional resonant bar tests allow the measurement of seismic properties of rocks and sediments at low frequencies (several kilohertz). However, the tests require a long, slender sample which is often difficult to obtain from the deep subsurface and weak and fractured formations. We present an alternative low-frequency measurement technique to the conventional resonant bar tests. This technique involves a jacketed core sample placed between a pair of long, metal extension rods with attached seismic source and receiver—the same geometry as the split Hopkinson pressure bar test for large-strain, dynamic impact experiments. Because of the added length and mass to the sample, the resonance frequency of the entire system can be lowered significantly, compared to the sample alone. The proposed “Split Hopkinson Resonant Bar (SHRB)” test is applied in two steps. In the first step, extension and torsion-mode resonance frequencies and attenuation of the system are measured. Then, numerical inversions for the compressional and shear wave velocities and attenuation are performed. We initially applied the SHRB test to synthetic materials (plastics) for testing its accuracy, then used it for measuring the seismic velocities and attenuation of a rock core containing supercritical CO 2 , and a sediment core while methane hydrate formed in the pore space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Second Generation Toolset for Calculation of Induced Seismicity Risk Profiles

This report describes development and demonstration of Generation 2 of the NRAP toolset for probabilistic analysis of risks from earthquakes that could potentially be induced by CO 2 injection. The capability to carry out hazard estimation at multiple sites and over multiple time periods of interested has been added to the integrated toolset. The capability to carry out hazard runs that utilize multiple simulations, based on random sampling of input parameter distributions, to bound hazard curve uncertainties has also been added. During Generation 2, an innovative method of constructing “nuisance fragility” functions for seismic ground motion has been developed. These functions are combined with seismic hazard curves to estimate of the risk of nuisance from ground shaking. The results described in this report demonstrate the functionality of the toolset based on its application to hypothetical scenarios involving one or two faults. Additional computer program development, including code parallelization, will be needed to bring the toolset up to full functionality. Alternative approaches to implementing the simulation-based approach to probabilistic induced seismicity hazard analysis are also outlined.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2014

Towards the Understanding of Induced Seismicity in Enhanced Geothermal Systems

This DOE funded project was a collaborative effort between Array Information Technology (AIT), the University of California at Berkeley (UCB), the Helmholtz Centre Potsdam - German Research Center for Geosciences (GFZ) and the Lawrence Berkeley National Laboratory (LBNL). It was also part of the European research project “GEISER”, an international collaboration with 11 European partners from six countries including universities, research centers and industry, with the goal to address and mitigate the problems associated with induced seismicity in Enhanced Geothermal Systems (EGS). The goal of the current project was to develop a combination of techniques, which evaluate the relationship between enhanced geothermal operations and the induced stress changes and associated earthquakes throughout the reservoir and the surrounding country rock. The project addressed the following questions: how enhanced geothermal activity changes the local and regional stress field; whether these activities can induce medium sized seismicity M > 3; (if so) how these events are correlated to geothermal activity in space and time; what is the largest possible event and strongest ground motion, and hence the potential hazard associated with these activities. The development of appropriate technology to thoroughly investigate and address these questions required a number of datasets to provide the different physical measurements distributed in space and time. Because such a dataset did not yet exist for an EGS system in the United State, we used current and past data from The Geysers geothermal field in northern California, which has been in operation since the 1960s. The research addressed the need to understand the causal mechanisms of induced seismicity, and demonstrated the advantage of imaging the physical properties and temporal changes of the reservoir. The work helped to model the relationship between injection and production and medium sized magnitude events that have jeopardized, and in some cases suspended, the generation of energy from EGS systems worldwide.

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

Seismic Attenuation Inversion with t* Using tstarTomog.

Seismic attenuation is defined as the loss of the seismic wave amplitude as the wave propagates excluding losses strictly due to geometric spreading. Information gleaned from seismic waves can be utilized to solve for the attenuation properties of the earth. One method of solving for earth attenuation properties is called t*. This report will start by introducing the basic theory behind t* and delve into inverse theory as it pertains to how the algorithm called tstarTomog inverts for attenuation properties using t* observations. This report also describes how to use the tstarTomog package to go from observed data to a 3-D model of attenuation structure in the earth.

Preston, Leiph↗
DOE OSTI2014

Development of Seismic Isolation Systems Using Periodic Materials

Advanced fast nuclear power plants and small modular fast reactors are composed of thin-walled structures such as pipes; as a result, they do not have sufficient inherent strength to resist seismic loads. Seismic isolation, therefore, is an effective solution for mitigating earthquake hazards for these types of structures. Base isolation, on which numerous studies have been conducted, is a well-defined structure protection system against earthquakes. In conventional isolators, such as high-damping rubber bearings, lead-rubber bearings, and friction pendulum bearings, large relative displacements occur between upper structures and foundations. Only isolation in a horizontal direction is provided; these features are not desirable for the piping systems. The concept of periodic materials, based on the theory of solid-state physics, can be applied to earthquake engineering. The periodic material is a material that possesses distinct characteristics that prevent waves with certain frequencies from being transmitted through it; therefore, this material can be used in structural foundations to block unwanted seismic waves with certain frequencies. The frequency band of periodic material that can filter out waves is called the band gap, and the structural foundation made of periodic material is referred to as the periodic foundation. The design of a nuclear power plant, therefore, can be unified around the desirable feature of a periodic foundation, while the continuous maintenance of the structure is not needed. In this research project, three different types of periodic foundations were studied: one-dimensional, two-dimensional, and three-dimensional. The basic theories of periodic foundations are introduced first to find the band gaps; then the finite element methods are used, to perform parametric analysis, and obtain attenuation zones; finally, experimental programs are conducted, and the test data are analyzed to verify the theory. This procedure shows that the periodic foundation is a promising and effective way to mitigate structural damage caused by earthquake excitation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Permafrost Active Layer Seismic Interferometry Experiment (PALSIE).

We present findings from a novel field experiment conducted at Poker Flat Research Range in Fairbanks, Alaska that was designed to monitor changes in active layer thickness in real time. Results are derived primarily from seismic data streaming from seven Nanometric Trillium Posthole seismometers directly buried in the upper section of the permafrost. The data were evaluated using two analysis methods: Horizontal to Vertical Spectral Ratio (HVSR) and ambient noise seismic interferometry. Results from the HVSR conclusively illustrated the method's effectiveness at determining the active layer's thickness with a single station. Investigations with the multi-station method (ambient noise seismic interferometry) are continuing at the University of Florida and have not yet conclusively determined active layer thickness changes. Further work continues with the Bureau of Land Management (BLM) to determine if the ground based measurements can constrain satellite imagery, which provide measurements on a much larger spatial scale.

54 ENVIRONMENTAL SCIENCES↗