Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗
DOE OSTI2025

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗
DOE OSTI2020

New Opportunities to Study Earthquake Precursors

The topic of earthquake prediction has a long history, littered with failed attempts. Part of the challenge is that possible precursory signals are usually reported after the event, and the systematic relationships between potential precursors and main events, should they exist, are unclear. Furthermore, several recent studies have shown the potential of new approaches to simultaneously detect earthquake foreshocks and slow-slip phenomena through ground deformation, seismic, and gravitational transients—weeks to months before large subduction zone earthquakes. The entire international community of earthquake researchers should be engaged in deploying instrumentation, sharing data in real time, and improving physical models to resolve the extent to which slow slip events and earthquake swarms enhance the likelihood (or not) for later, larger earthquakes.

58 GEOSCIENCES↗
DOE OSTI2026

Elastic-wave sensitivity-guided adaptive seismic survey design for cost-effective monitoring of geological carbon storage

Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments

58 GEOSCIENCES↗
DOE OSTI2022

Spatiotemporal Graph Convolutional Networks for Earthquake Source Characterization

Abstract Accurate earthquake location and magnitude estimation play critical roles in seismology. Recent deep learning frameworks have produced encouraging results on various seismological tasks (e.g., earthquake detection, phase picking, seismic classification, and earthquake early warning). Many existing machine learning earthquake location methods utilize waveform information from a single station. However, multiple stations contain more complete information for earthquake source characterization. Inspired by recent successes in applying graph neural networks (GNNs) in graph‐structured data, we develop a Spatiotemporal Graph Neural Network (STGNN) for estimating earthquake locations and magnitudes. Our graph neural network leverages geographical and waveform information from multiple stations to construct graphs automatically and dynamically by adaptive message passing based on graphs' edges. Using a recent graph neural network and a fully convolutional neural network as baselines, we apply STGNN to earthquakes recorded by the Southern California Seismic Network from 2000 to 2019 and earthquakes collected in Oklahoma from 2014 to 2015. STGNN yields more accurate earthquake locations than those obtained by the baseline models and performs comparably in terms of depth and magnitude prediction, though the ability to predict depth and magnitude remains weak for all tested models. Our work demonstrates the potential of using GNNs and multiple stations for better automatic estimation of earthquake epicenters.

58 GEOSCIENCES↗
DOE OSTI2023

Interseismic Strain Accumulation between the Colorado Plateau and the Eastern California Shear Zone: Implications for the Seismic Hazard near Las Vegas, Nevada

ABSTRACT Crustal deformation in the central Basin and Range between the Colorado plateau and the Eastern California Shear Zone is active but slow, making it a challenge to assess how strain is distributed and crustal motion transferred. However, knowledge of strain rates is very important, particularly for addressing the seismic hazard for both the Las Vegas urban area and the site of the proposed Yucca Mountain nuclear waste repository, in southern Nevada. Global Positioning System (GPS) data provide important constraints, particularly now that the GPS network in the area has substantially expanded in recent years. However, because deformation is slow, it is important to mitigate any transient tectonic and nontectonic signals to obtain the most accurate long-term interseismic motion and robust estimation of strain rates. We use data from all GPS stations in the region including both long-running continuous and semicontinuous stations. We model and remove postseismic displacements at these stations using source parameters for 41 events, dating back to the 1700 Cascadia megathrust earthquake, which contribute significantly to the deformation field within the central Basin and Range. We also remove correlated noise from the time series with the common-mode component imaging technique. We find that removal of both the postseismic transients and common-mode noise substantially reduces the uncertainties and spatial variation in the velocities. We find east–west extension across the Las Vegas Valley of 0.5–0.6 mm/yr. The interseismic strain rate field, calculated with the final velocities, reveals higher strain rates through southern Nevada than in previous studies, with rates within Las Vegas Valley of 8.5±2.4×10−9 yr−1. Our results also confirm shear along the Pahranagat shear zone, but the estimated amplitude is strongly affected by postseismic relaxation.

Geochemistry & Geophysics↗
DOE OSTI2024

LYNM-PE1 Seismic Parameters from Borehole Log, Laboratory, and Tabletop Measurements

The goal of this work is to provide a database of quality-checked seismic parameters that can be integrated with the Geologic Framework Model (GFM) for the LYNM-PE1 (Low Yield Nuclear Monitoring – Physical Experiment 1) testbed. We integrated data from geophysical borehole logs, tabletop measurements on collected core, and laboratory measurements. We reviewed for internal consistency among each measurement type, documented the caveats of measurement conditions, and integrated lithologic logs to check the validity of outlier values. The resulting consolidated parameter tables can be used as inputs for modeling and analysis codes and are designed to interface with the GFM, which is being actively developed.

58 GEOSCIENCES↗
DOE OSTI2024

Seismic Signal Detection on International Monitoring System 3-Component Stations using PhaseNet

In this report we discuss training a deep learning seismic signal detection model on 3-component stations from the International Monitoring System (IMS) using the PhaseNet architecture. Using 14 years of associated signals from the International Data Centre’s (IDC) Late Event Bulletin (LEB), we auto-curated training data consisting of signal windows containing associated arrivals, and noise windows that contain no LEB-associated signals. We trained several models using different waveform window durations (30 seconds and 100 seconds), with and without bandpass filtering. We evaluated the effectiveness of our models using associated signals from the Unconstrained Global Event Bulletin (UGEB) and found that several of our models outperformed the signal detections from the IDC’s Selected Event List 3 (SEL3) arrival table. The SEL3 bulletin evaluated on the UGEB dataset with 100-second waveform windows registered a precision and recall of .15 and .48, respectively, versus .19 and .59 for our filtered-data model. For the 30-second waveform window dataset, the SEL3 bulletin achieved a precision and recall of .31 and .47, respectively, versus .32 and .60 for our filtered-data model. Finally, our models detected signals from all source-to-receiver distances, suggesting it is feasible to use a single PhaseNet model for the IMS network.

58 GEOSCIENCES↗
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2025

Reducing the Parameter Dependency of Phase-Picking Neural Networks with Dice Loss

Training a neural network for picking seismic phase arrivals has been commonly posed as a segmentation problem. It is a highly imbalanced segmentation problem in the sense that the background vastly dominates the foreground because we are trying to pick the optimal single sample point that represents the arrival of a seismic phase in a many seconds long time window. Here, we test the Dice loss, which is a preferred loss function for highly imbalanced image segmentation problems. We show that phase-picking neural networks trained on the Dice loss behave in a binary fashion for which the prediction output is almost always either nearly 1 or nearly 0. This feature removes the strong dependence of data processing workflows on the prediction score threshold, which is an otherwise critical parameter to determine when using neural networks trained on the cross-entropy loss. When strategically used, models trained on the Dice loss can reduce the parameter dependency of machine learning-based seismic monitoring.

58 GEOSCIENCES↗
DOE OSTI2025

Quantifying uncertainty in regional-scale seismic moment tensors

We examine the ability of three different inversion methods: 1) first motion (FM) inversions, 2) amplitude inversions, and 3) full waveform (FW) time variable moment tensor (TVMT) inversions, to recover an accurate source mechanism for simulated data for an earthquake as well as an explosion. The ability of inversion methods, such as the ones described above, to recover accurate models representative of the data depends on both a priori information as well as the quality of data being inverted. Therefore, we examine the effect that station geometry, geologic model, and noise level has on the inversion and estimated source mechanism. We find that FM data can provide more robust solutions than amplitude data and are not as sensitive to inaccurate earth models, especially in low-noise cases, and overall FW inversions are the most accurate out of the methods examined, but can still be biased by inaccurate earth models.

58 GEOSCIENCES↗
DOE OSTI2020

The Complicated Link Between Material Properties and Microfracture Density for an Underground Explosion in Granite

Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (V p ), and S wave velocity (V s ), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.

58 GEOSCIENCES↗
DOE OSTI2024

Safety Basis Compliant Flexible Operations [Capstone Project]

Compliant nuclear operations must be weighed against the radiological risk to the worker, public, and the environment. This project set out to look at one facility at Los Alamos National Laboratory (LANL or the Laboratory) and compare the risk with the operational restrictions associated with the Technical Safety Requirements (TSRs). The hope was to identify controls that could be reduced or simplified while maintaining compliance with 10 CFR 830, Nuclear Safety Management. This would result in increasing facility operational flexibility and efficiency. Controls were identified to be simplified or eliminated based on risk. This identification was initially envisioned from the perspective of reanalysis of the accident and hazard analysis consequences. Additionally, it was envisioned that fire testing could show that the consequences of a fire are not as dramatic as analyzed. Specific facility controls were identified that cause operational inefficiencies. Based on data collection, two controls were identified for further evaluation. The first control was combustible loading within the waste storage buildings, and the second control is related to the seismic power cutoff switch. In evaluating the combustible loading control, the thought is that the fire suppression system (FSS) provides adequate protection and the combustible loading controls do not need elevation to a TSR control. The FSS, along with the Fire Protection Program, would provide protection to prevent a fire from impacting large amounts of stored waste material. The seismic switch was more of a challenge. When reviewing the facility operation, it was determined that the facility does not operate anywhere near the established material-at-risk (MAR – radioactive material) limits. By adjusting these limits, a reduced consequence could be achieved and hence, a significant risk reduction achieved. A new operating mode would be used to achieve this reduction. A reduced MAR limit was successfully used as a test case for snow removal during this project. This project expands on that trial case as related to combustible loading and the seismic power cutoff switch TSRs. This will increase the efficiency of operations with reduced TSR controls while reducing overall facility risk. This is a win-win situation for Triad National Security, LLC (Triad) and the U.S. Department of Energy (DOE). Due to schedule and budgetary constraints, not all envisioned work scope was accomplished. However, it is maintained in this report for future efforts.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Seismic Ambient Noise Analyses Reveal Changing Temperature and Water Signals to 10s of Meters Depth in the Critical Zone

The critical zone sustains terrestrial life, but we have few tools to explore it efficiently beyond the first few meters of the subsurface. Using analyses of high-frequency ambient seismic noise from densely spaced seismometers deployed in the forested Shale Hills subcatchment of the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), we show that temporal changes in seismic velocities at depths from ~1 m to tens of m can be detected. Furthermore, these changes are driven by variations at the land surface. The Moving-Window Cross-Spectral (MWCS) method was employed to measure seismic-velocity changes in coda waves at hourly resolution in 10 different frequency bands. We observed a diurnal signal, a seasonal signal, and a meteorological-event-based signal. These signals were compared to time-series measurements of precipitation, well water levels, soil moisture, soil temperature, air temperature, latent heat flux, and air pressure in the heavily instrumented catchment. Most of the velocity changes can be explained by variations in temperature that result in thermoelastic strains that propagate to depth. But some double minima in seismic velocity time-series observed after large rain events were attributed in part to the effects of water infiltration. These results show that high-frequency ambient noise data may in some locations be used to detect changes in the critical zone from ~1 to ~100 m or greater depth with hourly resolution. But interpretation of such data requires multiple environmental data sets to deconvolve the complex interrelationships among thermoelastic and hydrological effects in the subsurface critical zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

WUS256: An Adjoint Waveform Tomography Model of the Crust and Upper Mantle of the Western United States for Improved Waveform Simulations

Abstract We report a new model (WUS256) of radially anisotropic seismic wavespeeds of the crust and upper mantle of the western United States (WUS) obtained from adjoint waveform tomography for the purpose of improving synthetic waveform fits to observed data. WUS256 is based on inversion of over 94,000 waveforms from 72 earthquakes recorded by nearly 3,400 stations. We started with the SPiRaL global model (Simmons et al., 2021, https://doi.org/10.1093/gji/ggab277 ) and waveforms in the period band of 50–120 s. We followed a conservative multiscale inversion approach with eight stages and 256 total inversion iterations which enabled monotonic misfit reduction to 20‐s minimum‐period waves. WUS256 relied on time‐frequency (TF) phase misfits and a trust region limited memory Broyden–Fletcher–Goldfarb–Shanno (L‐BFGS) optimization. Hessian‐vector products were used to qualitatively assess model resolution. Results indicate that WUS256 has good coverage of the continental regions to depths of about 150 km and is able to resolve features on lateral scales of about 200 km. We quantify waveform fits by the reduction in TF and normalized amplitude difference misfits between WUS256 and the SPiRaL starting model. WUS256 significantly improves waveform fits with misfit reduction 64% for both inversion and validation data sets compared to the SPiRaL starting model and shows even better fits compared to other models. Waveform fits illustrate that WUS256 reproduces body‐waves, fundamental mode surface waves as well as late arriving dispersed and/or scattered short period surface waves. The improvement in waveform fit indicates that WUS256 can be used to reproduce path effects on regional complete waveforms and moment tensor inversions.

58 GEOSCIENCES↗
DOE OSTI2023

Daily and Seasonal Variations of Shallow Seismic Velocities in Southern California From Joint Analysis of H/V Ratios and Autocorrelations of Seismic Waveforms

Establishing a baseline of ongoing secular velocity variations at the subsurface can improve the accuracy of detecting and interpreting short-term velocity changes, and advance the understanding of observed seismic motions and the behavior of subsurface materials. Toward these goals, we develop and apply a deconvolved autocorrelation (DA) method to estimate regional daily and seasonal changes of seismic velocities in southern California. The DA method combines advantages of traditional autocorrelation and Horizontal-to-Vertical Spectral Ratio, and is used to analyze over 10 years of data recorded by 50 stations. The results indicate widespread daily and seasonal changes of up to 10% and 4%, respectively, in the top tens of meters of the crust. The thickness of the surface layer, distance from the coast, and topographic variations are important factors controlling the amplitudes of the resolved velocity variations. The results suggest that changes of soil moisture and thermoelastic strain are likely dominant factors affecting the daily and seasonal variations, respectively. The developed DA method can improve the accuracy and robustness of estimated changes of subsurface materials at other locations.

58 GEOSCIENCES↗