Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic modeling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Subducted slabs and the geoid - Constraints on mantle rheology and flow

The total geoid anomaly which is the result of a given density contrast in a convecting viscous earth is affected by the mass anomalies associated with the flow induced deformation of the upper surface and internal compositional boundaries, as well as by the density contrast itself is discussed. If the internal density contrasts can be estimated, the depth and variation of viscosity with depth of the convecting system can be constrained. The observed long wavelength geoid is highly correlated with that predicted by a density model for seismically active subducted slabs. The amplitude of the correlation is explained if the density contrasts associated with subduction extend into the lower mantle or if subducted slabs exceeding 350 km in thickness are piled up over horizontal distances of thousands of km at the base of the upper mantle. Mantle wide convection in a mantle that has a viscosity increasing with depth provides the explanation of the long-wavelength geoid anomalies over subduction zones. Previously announced in STAR as N83-22874

Hager, B. H.↗
NASA NTRS1993

Optimism Experiment and Development of Space-qualified Seismometers in France

The OPTIMISM experiment will put two magnetometers and two seismometers on the Martian floor in 1995, within the framework of the Mars '94 mission. The seismometers are put within the two small surface stations. The seismometer sensitivity will be better than 10 exp -9 g at 1 Hz, 2 orders of magnitude higher than the Viking seismometer sensitivity. A priori waveform modeling for seismic signals on Mars shows that it will be sufficient to detect quakes with a seismic moment greater than 10 exp 15 Nm everywhere on Mars. Such events, according to the hypothesis of a thermoelastic cooling of the Martian lithosphere, are expected to occur at a rate close to one per week and may therefore be observed within the l-year lifetime of the experiment. Other aspects of the experiment are discussed.

Lognonne, P.↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Modeling and Analysis for Spent Nuclear Fuel Seismic Testing

Note: This is the final public release version of PNNL-31671 Draft, which was previously released to the sponsor for review. The intent is the public release version will be PNNL-31671. There are no significant changes from the earlier draft version. This unlimited distribution report is the deliverable for M3SF-21PN010202014: Modeling and Analysis for Spent Nuclear Fuel Seismic Testing. This report summarizes the modeling, analysis, and test plan support completed by Pacific Northwest National Laboratory for the spent nuclear fuel dry storage system seismic test plan through May of 2021. Test plan preparation is planned to continue until the seismic test is completed in July of 2022. This report covers preliminary structural dynamic model development and computer-aided design of test hardware. Based on preliminary modeling, the strongest earthquakes under consideration in this test program provide mechanical loading on the spent nuclear fuel that is comparable to the 30 cm cask drop scenario, although the loads do not appear to be strong enough to cause significant permanent deformation of the fuel assembly spacer grids. The potential for grid deformation during the test will be assessed when the proposed shake table motion becomes available. The weakest earthquakes considered in this test are expected to be comparable to the mechanical loads witnessed in the multimodal transportation test of 2017. The horizontal canister system is predicted to provide nearly-uniform loading condition on the fuel assemblies it contains, while the vertical cask system is predicted to cause non-uniform dynamic loads on the fuel assemblies. In the horizontal case, gravity keeps the fuel assemblies in contact with one basket wall surface unless the loads are strong enough to cause a separation. In the vertical case, the fuel assemblies are relatively long and slender, and seismic motion in the anticipated test range is predicted to cause the assemblies to lean, tilt, and impact the basket walls throughout the seismic event. These gap closures are a nonlinear force transmission condition, so the vertical system is expected to have more variation and variability than the horizontal system. While the horizontal system is expected to have a relatively more linear response than the vertical system, there is still the potential for nonlinear behavior in the horizontal system because the fuel assemblies are free to slide, bounce, and impact the basket walls if the seismic loads are strong enough. One major conclusion of this study is that the use of mixed fuel assemblies in the canister will be acceptable. There are differences in overall system response when the mass, center of gravity, or gaps are changed, but the changes in response are within the bounds of a system that contains completely homogeneous fuel assemblies. One important observation is that the loading for each individual fuel assembly within the vertical canister is expected to be different from that of the others because of nonlinearities in the system. The horizontal canister system is expected to have a more uniform and more predictable response than the vertical canister, making variations in fuel assembly characteristics easier to account for. One important recommendation that comes from this modeling work is to repeat some of the strongest shake tests at different angles, particularly for the vertical cask system. Modeling predicts that strong seismic motion will cause the fuel assemblies to close initial gaps and impact the fuel basket walls when the canister is in the vertical cask configuration. Each shake test will impose a pre-defined three-dimensional time history on the cask system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Methods and systems for obtaining reconstructed low-frequency seismic data for determining a subsurface feature

A computer-implemented method for obtaining reconstructed seismic data for determining a subsurface feature, includes: determining an initial training velocity model, training a machine learning model based on first training seismic data and second training seismic data generated from the training velocity model, the first training seismic data corresponding to one or more first frequencies, the second training seismic data corresponding to one or more second frequencies lower than the one or more first frequencies, obtaining, based on measured seismic data and the machine learning model, reconstructed seismic data corresponding to the one or more second frequencies, generating a velocity model based on the measured seismic data, the reconstructed seismic data, and a full waveform inversion (FWI), and when the generated velocity model does not satisfy a preset condition, updating the training velocity model based on the generated velocity model, to obtain updated reconstructed seismic data for determining a subsurface feature.

58 GEOSCIENCES↗
DOE OSTI2020

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗
DOE OSTI2020

Regional Seismic Travel Time v.3.2.0

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗
DOE OSTI2026

Yield Estimates of three Historical Atmospheric Nuclear Explosions at Lop Nor, from Reduced Order Models of Regionally Recorded Rayleigh Waveforms and One Seismic Station: A Brief Communication

We use a reduced order model (ROM) for Rayleigh waveforms sourced by large atmospheric explosions along with data collected from a single seismic station (TLG) to estimate the yields of three atmospheric nuclear tests with order-megaton (MT) yields. These historical tests were conducted by China between 1973 and 1980 (CHIC 15, CHIC 16, and CHIC 26). We use our ROM to estimate yields of 2.34 MT (CHIC 15), 0.51 MT (CHIC 16), and 1.15 (CHIC 26) in the 20.5 s period band that are consistent with publicly accepted values of 3 MT (CHIC 15), 0.50 MT (CHIC 16), and 0.730 MT (CHIC 26).

58 GEOSCIENCES↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2021

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗
NASA NTRS1986

Finite difference simulations of seismic scattering Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity

Synthetic seismographs that were obtained by the finite difference method are presently applied to the study of elastic and acoustic wave scattering in two-dimensional media with random spatial variations in seismic velocity. The seismograms are analyzed to determine the variation in travel times and waveforms across arrays of receivers. The random media with Gaussian and exponential correlation functions considered differ in the spectral falloff of their velocity fluctuations at wavelengths smaller than 2pi times the correlation distance. It is found that alternative models of crustal heterogeneity can be tested by improved measurements of the frequency dependence of the crustal Q at frequencies greater than about 1 Hz, assuming that scattering is responsible for most of the attenuation at such frequencies.

Frankel, A.↗
DOE OSTI2022

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↗
DOE OSTI2020

Urban Seismic Site Characterization by Fiber-Optic Seismology

ccurate ground motion prediction requires detailed site effect assessment, but in urban areas where such assessments are most important, geotechnical surveys are difficult to perform, limiting their availability. Distributed acoustic sensing (DAS) offers an appealing alternative by repurposing existing fiber-optic cables, normally employed for telecommunication, as an array of seismic sensors. We present a proof-of-concept demonstration by using DAS to produce high-resolution maps of the shallow subsurface with the Stanford DAS array, California. We describe new methods and their assumptions to assess H/V spectral ratio—a technique widely used to estimate the natural frequency of the soil—and to extract Rayleigh wave dispersion curves from ambient seismic field. These measurements are jointly inverted to provide models of shallow seismic velocities and sediment thicknesses above bedrock in central campus. The good agreement with an independent survey validates the methodology and demonstrates the power of DAS for microzonation.

Geochemistry & Geophysics↗
DOE OSTI2025

Seismicity-constrained fault detection and characterization with a multitask machine learning model

Geological fault detection and characterization are crucial for understanding subsurface dynamics across scales. While methods for fault delineation based on either seismicity location analysis or seismic image reflector discontinuity are well-established, a systematic approach that integrates both data types remains absent. We develop a novel machine learning model that unifies seismic reflector images and seismicity location information to automatically identify geological faults and characterize their geometrical properties. The model encodes a seismic image and a seismicity location image separately, and fuses the encoded features with a spatial-channel attention fusion module to improve the learning of important features in both inputs. We design an automated strategy to generate high-quality synthetic training data and labels. To improve the realism of the seismicity location image, we include random seismicity noise and missing seismicity location associated with some of the faults. We validate the model’s efficacy and accuracy using synthetic data examples and two field data examples. Moreover, we show that fine-tuning the trained model with a small, domain-specific dataset enhances its fidelity for field data applications. The results demonstrate that integrating seismicity location and seismic images into a unified framework allows the end-to-end neural network to achieve higher fidelity and accuracy in delineating subsurface faults and their geometrical properties compared with image-only fault detection methods. Our approach offers an adaptive data-driven tool for geological fault characterization and seismic hazard mitigation, bridging the gap between seismicity location and image-based fault detection methods.

58 GEOSCIENCES↗