Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMIC”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Seismic Waves, 4th order accurate

SW4 is a program for simulating seismic wave propagation on parallel computers. SW4 colves the seismic wave equations in Cartesian corrdinates. It is therefore appropriate for regional simulations, where the curvature of the earth can be neglected. SW4 implements a free surface boundary condition on a realistic topography, absorbing super-grid conditions on the far-field boundaries, and a kinematic source model consisting of point force and/or point moment tensor source terms. SW4 supports a fully 3-D heterogeneous material model that can be specified in several formats. SW4 can output synthetic seismograms in an ASCII test format, or in the SAC finary format. It can also present simulation information as GMT scripts, whixh can be used to create annotated maps. Furthermore, SW4 can output the solution as well as the material model along 2-D grid planes.

Source record↗
DOE OSTI2023

Seismic Response Processing Module

The seismic response processing module is a Java based library that provides support for removing instrument response signals from seismic recordings. The response module provides support for EvalResp, PZF, PAZ, FAP, PAZFIR, PAZFAP, FIRFAP, and CSS response types and parallel processing support for transfer operations. Additionally, the response module uses the JSR-363 units of measurement specification to allow for transfer to and from a wide range of units types to represent the time-series.

Dodge, Douglas↗
DOE OSTI2024

Bayesian OED for Seismic Monitoring

SAND2024-13870O The Bayesian OED (Optimal Experiment Design) for Seismic Monitoring code provides the tools to analyze and optimize seismic monitoring networks using Bayesian OED. This method designs a utility function for experiments (network designs) using network analysis and network optimization. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Catanach, Thomas↗
DOE OSTI2016

A source-synchronous filter for uncorrelated receiver traces from a swept-frequency seismic source

We have developed a novel algorithm to reduce noise in signals obtained from swept-frequency sources by removing out-of-band external noise sources and distortion caused from unwanted harmonics. The algorithm is designed to condition nonstationary signals for which traditional frequency-domain methods for removing noise have been less effective. The source synchronous filter (SSF) is a time-varying narrow band filter, which is synchronized with the frequency of the source signal at all times. Because the bandwidth of the filter needs to account for the source-to-receiver propagation delay and the sweep rate, SSF works best with slow sweep rates and moveout-adjusted waveforms to compensate for source-receiver delays. The SSF algorithm was applied to data collected during a field test at the University of California Santa Barbara’s Garner Valley downhole array site in Southern California. At the site, a 45 kN shaker was mounted on top of a one-story structure and swept from 0 to 10 Hz and back over 60 s (producing useful seismic waves greater than 1.6 Hz). The seismic data were captured with small accelerometer and geophone arrays and with a distributed acoustic sensing array, which is a fiber-optic-based technique for the monitoring of elastic waves. The result of the application of SSF on the field data is a set of undistorted and uncorrelated traces that can be used in different applications, such as measuring phase velocities of surface waves or applying convolution operations with the encoder source function to obtain traveltimes. Lastly, the results from the SSF were used with a visual phase alignment tool to facilitate developing dispersion curves and as a prefilter to improve the interpretation of the data.

58 GEOSCIENCES↗
DOE OSTI2020

Gas hydrate quantification in Walker Ridge block 313, Gulf of Mexico, from full-waveform inversion of ocean-bottom seismic data

The Gulf of Mexico (GOM) Joint Industry Project Leg 2 logging-while-drilling data in Walker Ridge lease block 313 (WR313) in the GOM detected gas hydrate in coarse- and fine-grained sediments at sites WR313-G and WR313-H. The coarse-grained units are thin ([Formula: see text]) and highly saturated, whereas the fine-grained unit is thick (approximately 200 m) with low saturation and fracture-filling gas hydrate. Unlike its coarse-grained counterpart, the seismic character of the fine-grained unit does not clearly indicate the presence of gas hydrate, which would likely have remained undiscovered in the absence of drilling. In this paper, through frequency-domain acoustic full-waveform inversion (FWI) of ocean-bottom seismometer data along a 2D multichannel seismic transect near sites WR313-G and WR313-H, we detect and quantify gas hydrate in the fine-grained unit. Key results are as follows: First, the base of the gas hydrate stability zone, which is not obvious in the reflection profile, can be discerned in the FWI results. Second, the gas hydrate in the fine-grained unit is mainly confined to the area between two sets of opposite-dipping normal faults implying that the fault architecture may be partially responsible for this gas hydrate accumulation and distribution.

Geochemistry & Geophysics↗
DOE OSTI2022

Seismic monitoring of well integrity

Seismic tube waves, produced by flexure of the well boundary, pressure changes in the fluid in the well, and deformation of the material immediately surrounding the well, are particularly sensitive to variations in the state of the well. In this work, we evaluate a direct approach for generating and observing tube waves as a means of detecting well damage. While we find that it can be difficult to reliably excite observable tube waves without a very strong surface source, time-frequency techniques can be employed to increase the detectability of tube wave reflections. New technologies, particularly distributed acoustic sensing, hold great promise for evaluating well integrity by monitoring tube waves, temperature changes, and seismic noise due to well deformation and fluid leakage.

58 GEOSCIENCES↗
DOE OSTI2023

Utah FORGE 5-2419: Seismicity-Permeability Relationships Probed via Nonlinear Acoustic Imaging - Workshop Presentation

This is a presentation on the Seismicity-Permeability Relationships Probed via Nonlinear Acoustic Imaging project by Pennsylvania State University, presented by Derek Elsworth. The project's objectives were to explore controls and acoustic signatures of aseismic through seismic evolution of friction-stability-permeability relationships on fractures in shear-reactivation, link this reactivation to key features of the pre-existing stress state, and upscale these indexes to reservoir scale to drive successful reservoir stimulation. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

Significantly Improving Regional Seismic Amplitude Tomography at Higher Frequencies by Determining S -Wave Bandwidth

Characterizing regional seismic signals continues to be a difficult problem due to their variability. Calibration of these signals is very important to many aspects of monitoring underground nuclear explosions, including detecting seismic signals, discriminating explosions from earthquakes, and reliably estimating magnitude and yield. Amplitude tomography, which simultaneously inverts for source, propagation, and site effects, is a leading method of calibrating these signals. A major issue in amplitude tomography is the data quality of the input amplitude measurements. Pre-event and prephase signal-to-noise ratio (SNR) tests are typically used but can frequently include bad signals and exclude good signals. The deficiencies of SNR criteria, which are demonstrated here, lead to large calibration errors. To ameliorate these issues, we introduce a semi-automated approach to assess the bandwidth of a spectrum where it behaves physically. We determine the maximum frequency (denoted as F max ) where it deviates from this behavior due to inflections at which noise or spurious signals start to bias the spectra away from the expected decay. We compare two amplitude tomography runs using the SNR and new F max criteria and show significant improvements to the stability and accuracy of the tomography output for frequency bands higher than 2 Hz by using our assessments of valid S-wave bandwidth. We compare Q estimates, P/S residuals, and some detailed results to explain the improvements. Lastly, for frequency bands higher than 4 Hz, needed for effective P/S discrimination of explosions from earthquakes, the new bandwidth criteria sufficiently fix the instabilities and errors so that the residuals and calibration terms are useful for application.

58 GEOSCIENCES↗
DOE OSTI2018

Dynamic Tuning of Seismic Signal Detector Trigger Levels for Local Networks

The quality of automatic signal detections from sensor networks depends on individual detector trigger levels (TLs) from each sensor. The largely manual process of identifying effective TLs is painstaking and does not guarantee optimal configuration settings, yet achieving superior automatic detection of signals and ultimately, events, is closely related to these parameters. In this paper, we present a Dynamic Detector Tuning (DDT) system that automatically adjusts effective TL settings for signal detectors to the current state of the environment by leveraging cooperation within a local neighborhood of network sensors. After a stabilization period, the DDT algorithm can adapt in near-real time to changing conditions and automatically tune a signal detector to identify (detect) signals from only events of interest. Our current work focuses on reducing false signal detections early in the seismic signal processing pipeline, which leads to fewer false events and has a significant impact on reducing analyst time and effort. This system provides an important new method to automatically tune detector TLs for a network of sensors and is applicable to both existing sensor performance boosting and new sensor deployment. Finally, with ground truth on detections from a local neighborhood of seismic sensors within a network monitoring the Mount Erebus volcano in Antarctica, we show that DDT reduces the number of false detections by 18% and the number of missed detections by 11% when compared with optimal fixed TLs for all sensors.

58 GEOSCIENCES↗
DOE OSTI2018

Depth Discrimination Using Rg-to-Sg Spectral Amplitude Ratios for Seismic Events in Utah Recorded at Local Distances

Short-period fundamental-mode Rayleigh waves (Rg) are commonly observed on seismograms of anthropogenic seismic events and shallow, naturally occurring tectonic earthquakes (TEs) recorded at local distances. In the Utah region, strong Rg waves traveling with an average group velocity of about 1.8 km/s are observed at ~1 Hz on waveforms from shallow events ( depth<10 km ) recorded at distances up to about 150 km. At these distances, Sg waves, which are direct shear waves traveling in the upper crust, are generally the dominant signals for TEs. Here in this study, we leverage the well-known notion that Rg amplitude decreases dramatically with increasing event depth to propose a new depth discriminant based on Rg-to-Sg spectral amplitude ratios. The approach is successfully used to discriminate shallow events (both earthquakes and anthropogenic events) from deeper TEs in the Utah region recorded at local distances ( <150 km ) by the University of Utah Seismographic Stations (UUSS) regional seismic network. Using Mood’s median test, we obtained probabilities of nearly zero that the median Rg-to-Sg spectral amplitude ratios are the same between shallow events on the one hand (including both shallow TEs and anthropogenic events), and deeper earthquakes on the other, suggesting that there is a statistically significant difference in the estimated Rg-to-Sg ratios between the two populations. We also observed consistent disparities between the different types of shallow events (e.g., mining blasts vs. mining-induced earthquakes), implying that it may be possible to separate the subpopulations that make up this group. Lastly, this suggests that using local distance Rg-to-Sg spectral amplitude ratios one can not only discriminate shallow events from deeper events but may also be able to discriminate among different populations of shallow events.

58 GEOSCIENCES↗
DOE OSTI2019

Classification of Local Seismic Events in the Utah Region: A Comparison of Amplitude Ratio Methods with a Spectrogram-Based Machine Learning Approach

Abstract The capability to discriminate low-magnitude earthquakes from low-yield anthropogenic sources, both detectable only at local distances, is of increasing interest to the event monitoring community. Here, we used a dataset of seismic events in Utah recorded during a 14-day period (1–14 January 2011) by the University of Utah Seismic Stations network to perform a comparative study of event classification at local scale using amplitude ratio (AR) methods and a machine learning (ML) approach. The event catalog consists of 7377 events with magnitudes M C ranging from –2 and lower up to 5.8. Events were subdivided into six populations based on location and source type: tectonic earthquakes (TEs), mining-induced events (MIEs), and mining blasts from four known mines (WMB, SMB, LMB, and CQB). The AR approach jointly exploits Pg-to-Sg phase ARs and Rg-to-Sg spectral ARs in multivariate quadratic discriminant functions and was able to classify 370 events with high signal quality from the three groups with sufficient size (TE, MIE, and SMB). For that subset of the events, the method achieved success rates between about 80% and 90%. The ML technique used trained convolutional neural network (CNN) models to classify the populations. The CNN approach was able to classify the subset of events with accuracies between about 91% and 98%. Because the neural network approach does not have a minimum signal quality requirement, we applied it to the entire event catalog, including the abundant extremely low-magnitude events, and achieved accuracies of about 94%–100%. Furthermore, we compare the AR and ML methodologies using a broad set of criteria and conclude that a major advantage to ML methods is their robustness to low signal-to-noise ratio data, allowing them to classify significantly smaller events.

58 GEOSCIENCES↗
DOE OSTI2021

Deep Learning Denoising Applied to Regional Distance Seismic Data in Utah

Seismic waveform data are generally contaminated by noise from various sources. Suppressing this noise effectively so that the remaining signal of interest can be successfully exploited remains a fundamental problem for the seismological community. To date, the most common noise suppression methods have been based on frequency filtering. These methods, however, are less effective when the signal of interest and noise share similar frequency bands. Inspired by source separation studies in the field of music information retrieval (Jansson et al., 2017) and a recent study in seismology (Zhu et al., 2019), we implemented a seismic denoising method that uses a trained deep convolutional neural network (CNN) model to decompose an input waveform into a signal of interest and noise. In our approach, the CNN provides a signal mask and a noise mask for an input signal. The short-time Fourier transform (STFT) of the estimated signal is obtained by multiplying the signal mask with the STFT of the input signal. To build and test the denoiser, we used carefully compiled signal and noise datasets of seismograms recorded by the University of Utah Seismograph Stations network. Results of test runs involving more than 9000 constructed waveforms suggest that on average the denoiser improves the signal-to-noise ratios (SNRs) by ~5 db, and that most of the recovered signal waveforms have high similarity with respect to the target waveforms (average correlation coefficient of ~0.80) and suffer little distortion. Application to real data suggests that our denoiser achieves on average a factor of up to ~2-5 improvement in SNR over band-pass filtering and can suppress many types of noise that band-pass filtering cannot. For individual waveforms, the improvement can be as high as ~15 db.

58 GEOSCIENCES↗
DOE OSTI2022

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗
DOE OSTI2023

Testing and Design of Discriminants for Local Seismic Events Recorded during the Redmond Salt Mine Monitoring Experiment

The Redmond Salt Mine (RSM) Monitoring Experiment in Utah was designed to record seismoacoustic data at distances less than 50 km for algorithm testing and development. During the experiment from October 2017 to July 2019, six broadband seismic stations were operating at a time, with three of them having fixed locations for the duration, whereas the three other stations were moved to different locations every one-and-half to two-and-half months. RSM operations consist of nighttime underground blasting several times per week. The RSM is located in proximity to a belt of active seismicity, allowing direct comparison of natural and anthropogenic sources. Using the recorded data set, we built 1373 events with local magnitude (M L ) of -2.4 and lower to 3.3. For 75 blasts (RMEs) from the Redmond Salt Mine and 206 tectonic earthquakes (EQs), both M L and the coda duration magnitude (MC) are well constrained. We used these events to test and design discriminants that separate the RMEs from the EQs and are effective at local distances. The discriminants consist of M L -M C , low-frequency Sg to high-frequency Sg, Pg/Sg phase-amplitude ratios, and Rg/Sg spectral amplitude ratios, as well as different combinations of two or more of these classifiers. The areas under the receiver operating characteristic curves (AUCs) of 0.92–1.0 for M L -M C , low-frequency Sg to high-frequency Sg, and Rg/Sg indicate that these discriminants are very effective. Conversely, the AUC of only 0.57 for Pg/Sg suggests that this discriminant is only slightly better than a random classifier. Among the effective classifiers, Rg/Sg, shows the lowest likelihood of misclassification (4.3%) for the populations. In conclusion, results of joint discriminant analyses suggest that even the arguably ineffective single classifier, like Pg/Sg in this case, can provide some value when used in combination with others.

58 GEOSCIENCES↗
DOE OSTI2016

Iterative Strategies for Aftershock Classification in Automatic Seismic Processing Pipelines

We report aftershock sequences following very large earthquakes present enormous challenges to near-real-time generation of seismic bulletins. The increase in analyst resources needed to relocate an inflated number of events is compounded by failures of phase-association algorithms and a significant deterioration in the quality of underlying, fully automatic event bulletins. Current processing pipelines were designed a generation ago, and, due to computational limitations of the time, are usually limited to single passes over the raw data. With current processing capability, multiple passes over the data are feasible. Processing the raw data at each station currently generates parametric data streams that are then scanned by a phase-association algorithm to form event hypotheses. We consider the scenario in which a large earthquake has occurred and propose to define a region of likely aftershock activity in which events are detected and accurately located, using a separate specially targeted semiautomatic process. This effort may focus on so-called pattern detectors, but here we demonstrate a more general grid-search algorithm that may cover wider source regions without requiring waveform similarity. Given many well-located aftershocks within our source region, we may remove all associated phases from the original detection lists prior to a new iteration of the phase-association algorithm. We provide a proof-of-concept example for the 2015 Gorkha sequence, Nepal, recorded on seismic arrays of the International Monitoring System. Even with very conservative conditions for defining event hypotheses within the aftershock source region, we can automatically remove about half of the original detections that could have been generated by Nepal earthquakes and reduce the likelihood of false associations and spurious event hypotheses. Lastly, further reductions in the number of detections in the parametric data streams are likely, using correlation and subspace detectors and/or empirical matched field processing.

58 GEOSCIENCES↗
DOE OSTI2020

Effects of Shallow-Velocity Reductions on 3D Propagation of Seismic Waves

Abstract We perform 3D simulations of seismic wavefields to clarify effects of strong reductions of shallow velocities on long-period seismic waves. The simulations use a reference Community Velocity Model of southern California and a modified version with strong velocity reductions in the top 500 m of the Los Angeles basin. Differences between wavefields generated by 10 earthquakes in the reference and perturbed models are analyzed. Velocity changes are estimated by measuring relative time shifts between reference and perturbed seismograms using wavelet cross-correlation spectra. The results indicate that strong localized temporal velocity drops near the surface, such as those observed during strong ground motions, may generate regional perturbations of wavefields at periods up to 20 s. These perturbations may be misinterpreted as generated by temporal changes at seismogenic depths. The results also have important implications for waveform tomography studies.

Geochemistry & Geophysics↗
DOE OSTI2021

Pycheron: A Python-Based Seismic Waveform Data Quality Control Software Package

Supplementing an existing high-quality seismic monitoring network with openly available station data could improve coverage and decrease magnitudes of completeness; however, this can present challenges when varying levels of data quality exist. Without discerning the quality of openly available data, using it poses significant data management, analysis, and interpretation issues. Incorporating additional stations without properly identifying and mitigating data quality problems can degrade overall monitoring capability. If openly available stations are to be used routinely, a robust, automated data quality assessment for a wide range of quality control (QC) issues is essential. To meet this need, we developed Pycheron, a Python-based library for QC of seismic waveform data. Pycheron was initially based on the Incorporated Research Institutions for Seismology’s Modular Utility for STAtistical kNowledge Gathering but has been expanded to include more functionality. Pycheron can be implemented at the beginning of a data processing pipeline or can process stand-alone data sets. Its objectives are to (1) identify specific QC issues; (2) automatically assess data quality and instrumentation health; (3) serve as a basic service that all data processing builds on by alerting downstream processing algorithms to any quality degradation; and (4) improve our ability to process orders of magnitudes more data through performance optimizations. This article provides an overview of Pycheron, its features, basic workflow, and an example application using a synthetic QC data set.

58 GEOSCIENCES↗
DOE OSTI2021

Analysis of Seismic Signals Generated by Vehicle Traffic with Application to Derivation of Subsurface Q-Values

Correct identification and modeling of anthropogenic sources of ground motion are of considerable importance for many studies, including detection of small earthquakes and imaging seismic properties below the surface. Here, to understand signals generated by common vehicle traffic, we use seismic data recorded by closely spaced geophones normal to roads at two sites on San Jacinto fault zone. To quantify the spatiotemporal and frequency variations of the recorded ground motions, we develop a simple analytical solution accounting for propagation and attenuation of surface waves. The model reproduces well-observed bell-shaped spectrograms of car signals recorded by geophones close to roads, and it can be used to estimate frequency-dependent Q-values of the subsurface materials. The data analysis indicates Q-values of 3–40, for frequencies up to 150 Hz for road-receiver paths at the two examined sites. The derived Q-values are consistent with attenuation factors of surface waves previously obtained with other methods. The analytical results and analysis procedure provide a highly efficient method for deriving Q-values of shallow subsurface materials.

58 GEOSCIENCES↗