Search NASA⌕ Search

Engineering topics

Preston, Leiph

Publications and source records attributed to Preston, Leiph.

At least 19 records

Deep Learning for Full Waveform Inversion of Elastic Active-Source Seismic Data to Estimate P-Wave Velocity Models

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗

Modeling of Seismic Waves Through Geologic Metamaterials

This project conducted a modeling study on seismic invisibility cloaks that render geologic targets invisible to seismic waves, using the concept of seismic metamaterials. We present a parametric numerical study on the behaviors of seismic waves through cloaks with different design parameters as well as degrees of geologic heterogeneity. In addition, a seismic cloaking strategy is proposed for a future field-scale experiment at a real-world test bed. This feasibility study will guide future field experiment designs and ultimately allow us to conduct systematic field-scale tests employing Sandia’s existing resources and field expertise. The ultimate goal is to develop methods and design parameters of seismic invisibility cloaks to protect against natural and man-made seismic waves. Seismic cloaking has potential applications in several areas of national security, energy, and natural hazard reduction.

58 GEOSCIENCES↗

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.

58 GEOSCIENCES↗

NetCDFFile v.1.0

SAND2024-02969O The NetCDFFile provides a generic interface to NetCDF and/or PnetCDF library calls, allowing the user to write generic code to read and write from NetCDF files. The NetCDFFile library requires a C++20 compliant compiler as well as NetCDF and PnetCDF libraries installed on the machine. 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.

Preston, Leiph↗

Rock Valley Direct Comparison Relocation Working Group Location Results and Recommendations

Work accomplished: Collected and compared historic data for the 1993 Rock Valley earthquake sequence; Compared preliminary and prior location work from different location algorithms, phase pick sets, station constellations, and velocity models; Selected a common set of stations that could be used across all location methods for consistency; Reviewed 8 different sets of phase picks and converged on a single, reviewed set of picks for all common stations; Evaluated four pre-existing regional velocity models and incorporated new and preliminary results for five new velocity models that provide information on the very shallow (< 2km) structure near station RTPP; Compared location results from different methods while using the common sets of picks, stations, and velocity models

58 GEOSCIENCES↗

Rock Valley Direct Comparison Relocation Working Group Location Results and Recommendations (Version 4.0)

The Source Physics Experiments (SPE) were designed to improve our physics-based understanding of explosion sources for the purposes of nuclear test monitoring. Phase I consisted of 6 chemical explosions in the Climax Stock Granite of the Nevada National Security Site (NNSS), while Phase II consisted of 4 explosions in a contrasting dry alluvium geology (DAG) in Yucca Flat, providing essential data in various media and emplacement conditions to further modeling efforts. For Phase III, the Rock Valley Direct Comparison (RVDC) seeks to directly compare earthquake and explosion source types. An unusually shallow series of events in 1993 along the Rock Valley Fault Zone in the southeastern portion of the NNSS has been targeted for this direct comparison. Depth ranges for the events, previously estimated to be less than 3 km, is achievable by modern drilling techniques and accessibility to the epicentral locations would require minimal improvements to the infrastructure. The events providing this unique opportunity for direct comparison are the focus of this report.

58 GEOSCIENCES↗

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗