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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Geologic Framework Model for the Dry Alluvium Geology (DAG) Experiment Testbed Yucca Flat, Nevada National Security Site

Geologic framework models (GFMs) provide a methodology for integrating geology into, and thus geologically informing, other modeling and simulation activities. GFMs provide a threedimensional (3-D), geology-based digital framework for developing and parametrizing meshes and evaluating model and simulation results. This report describes a 3-D GFM constructed for the Source Physics Experiment Phase II Dry Alluvium Geology test series located in Yucca Flat at the Nevada National Security Site. The geology in the Yucca Flat region is complex and diverse, which creates challenges to modeling seismic wave propagation from SPE tests. The Yucca Flat GFM helps address these challenges by providing the 3-D distribution of relevant geologic features and physical properties necessary to more effectively model seismic wave propagation. The GFM includes 7 model layers and 48 faults that cut and offset the layers. An appendix is included that provides quantitative data on physical properties for each model layer.

58 GEOSCIENCES↗
DOE OSTI2020

Geologic framework model for the Dry Alluvium Geology (DAG) experiment testbed, Yucca Flat, Nevada National Security Site

Geologic framework models (GFMs) provide a methodology for integrating geology into, and thus geologically informing, other modeling and simulation activities. GFMs provide a three-dimensional (3-D), geology-based digital framework for developing and parametrizing meshes and evaluating model and simulation results. This report describes a 3-D GFM constructed for the Source Physics Experiment Phase II Dry Alluvium Geology test series located in Yucca Flat at the Nevada National Security Site. The geology in the Yucca Flat region is complex and diverse, which creates challenges to modeling seismic wave propagation from SPE tests. The Yucca Flat GFM helps address these challenges by providing the 3-D distribution of relevant geologic features and physical properties necessary to more effectively model seismic wave propagation. The GFM includes 7 model layers and 48 faults that cut and offset the layers. An appendix is included that provides quantitative data on physical properties for each model layer.

58 GEOSCIENCES↗
DOE OSTI2021

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1988

Lunar composition - A geophysical and petrological synthesis

Lunar compositional constraints are derived on the basis of geophysical data (in particular, the lunar seismic model as revised by Nakamura) and petrological arguments. Only in the case of extreme assumptions can critical aspects of bulk lunar composition be demonstrated to be equivalent to the present-day terrestrial mantle; specifically, the moon has an Mg number that is too low and an alumina abundance that is too high. Over a broad range of crustal densities the presence of a metallic core at least 150 km in radius is necessary to reconcile the upper mantle lunar seismic model with mass and moment of inertia constraints. A significant seismic discontinuity at 500 km depth may mark the lowest extent of early lunar differentiation, possibly representing a transition between highly fractionated upper mantle and less fractionated, perhaps even primordial, middle/lower mantle.

Mueller, Steve↗
DOE OSTI2026

Seismic velocity modeling and earthquake relocations for the southern Nevada National Security Site (NNSS), with a focus on Rock Valley

The Rock Valley Direct Comparison (RV/DC) is the third phase of the Source Physics Experiment (SPE; Snelson et al., 2014), a project aimed at improving seismic discrimination between explosive sources and natural earthquakes. Phases I and II of SPE focused on studying the generation of shear waves in endmember geologies (granite and alluvium, respectively), while Phase III aims to detonate a chemical explosion co-located with the source region of the 1993 Rock Valley earthquake sequence, located in the Rock Valley Fault Zone (RVFZ), in the southern portion of the Nevada Test Site (now the Nevada National Security Site, or NNSS).

58 GEOSCIENCES↗
DOE OSTI2021

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗
DOE OSTI2021

Density of Fe-Ni-C Liquids at High Pressures and Implications for Liquid Cores of Earth and the Moon

The presence of light elements in the metallic cores of the Earth, the Moon, and other rocky planetary bodies has been widely proposed. Carbon is among the top candidates in light of its high cosmic abundance, siderophile nature, and ubiquity in iron meteorites. It is, however, still controversial whether carbon-rich core compositional models can account for the seismic velocity observations within the Earth and lunar cores. In this paper we report the density and elasticity of Fe 90 Ni 10 -3 wt.% C and Fe 90 Ni 10 -5 wt.% C liquid alloys using synchrotron-based X-ray absorption experiments and first-principles molecular dynamics simulations. Our results show that alloying of 3 wt.% and 5 wt.% C lowers the density of Fe 90 Ni 10 liquid by ~2.9-3.1% at 2 GPa, and ~3.4-3.6% at 9 GPa. More intriguingly, our experiments and simulations both demonstrate that the bulk moduli of the Fe-Ni-C liquids are similar to or slightly higher than those of Fe-Ni liquids. Thus, the calculated compressional velocity ($v_p$) of Fe-Ni-C liquids are higher than that of pure Fe-Ni alloy, promoting carbon as a possible candidate to explain the elevated v p in the Earth's outer core. However, the values and slopes of both density and v p of the studied two Fe-Ni-C liquids do not match the outer core seismic models, suggesting that carbon may not be the sole principal light element in Earth's outer core. The high $v_p$ of Fe-Ni-C liquids does not match the presumptive $v_p$ of the lunar outer core well, indicating that carbon is less likely to be its dominant light element.

58 GEOSCIENCES↗
DOE OSTI2024

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗
DOE OSTI2022

Distinguishing Unique Earthquakes with Overlapping Signals in Oklahoma

Abstract During routine operations monitoring Oklahoma earthquakes, we found that certain earthquakes occurred closely both in space and time and had overlapping phases at the recording stations. Through further scrutiny and analysis, we determined that rather than being distinctly different earthquakes, some of the earthquakes exhibited multiphase arrivals and longer than expected coda due to unique ray paths that encounter impedance contrasts such as at the sedimentary rock-basement. Of course, some of these events truly were distinct events, which we term overlapping earthquakes, for which perceived coda duration overlaps and obscures the phase arrivals of the second event due to the source proximity in both time and space. We detail our classification scheme to separate the local earthquakes in Oklahoma as single, overlapping earthquakes, or those associated with multiphase arrivals. We forward model seismic wave propagation in a 2D crustal model and develop a methodology that utilizes waveform correlation to distinguish phases from overlapping earthquakes to those from crustal reverberations. Duration analysis shows a more elongated duration, qualitatively similar to the duration produced by overlapping earthquakes, at the sites where multiphase arrivals are observed.

Geochemistry & Geophysics↗
NASA NTRS1976

Seismic structure of the moon - A summary of current status

A seismic model of the moon, refined from previous models with the aid of new data and analyses for seismic events recorded by four Apollo lunar seismic stations and a gravimeter, is outlined. The current model has five zones. The crust is 55-60 km thick. The surface is covered by regolith of varying thicknesses and low seismic velocity. The seismic Q's are extremely high. The upper mantle, 250 km thick, is characterized by seismic velocities close to those of the earth's upper mantle, and shear wave Q of about 4000. Tidally controlled deep moonquakes are concentrated near the boundary of the middle mantle, 600-700 km thick. A lower mantle extending to an undetermined depth is characterized by high attenuation of shear waves. The existence of a low-velocity core is only tentatively proposed.

Nakamura, Y.↗
DOE OSTI2020

Updates to the Regional Seismic Travel Time (RSTT) Model: 1. Tomography

Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).

58 GEOSCIENCES↗
DOE OSTI2021

Wellbore integrity assessment with casing-based advanced sensing

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage 1 and waste disposal. After installation, well casing and cement are subject to mechanical stress due to near-well pressure changes and fluid induced corrosion. This is exacerbated for geothermal wells where produced fluid is at high temperature and corrosive. The current state-of-the-art technologies for wellbore integrity assessments are an array of cased hole logging tools. Wireline deployed acoustic, electromagnetic and mechanical tools are all available to inspect steel casing corrosion and casing-cement bond and these tools can provide high-resolution assessment of borehole conditions. They are intrusive, however, in terms of borehole preparation and interruption to the normal operation of the wells, and not suitable for high temperature or highly deviated well deployments. In addition, these measurements are performed infrequently due to high cost, and are therefore incapable of providing frequent data to better predict borehole degradation trajectory, which can help provide early warning of potential borehole failures. For this project we are developing a suite of novel, non-invasive, casing based tools for wellbore integrity assessment, combining fast/low cost screening with higher-precision investigation. Our approach is based on monitoring the response of the casing when energized at the wellhead, thereby interrogating the casing without well intervention. Lab, field and numerical approaches are used in our study. During the early stage of the research, we focus on numerical simulations, which have shown the sensitivity of the low frequency electromagnetic (EM) signals to changes in borehole depths and have successfully tested the concept at a field site with different length well casings. Initial seismic modeling efforts have also demonstrated our capability to simulate seismic tube wave and seismic field alterations due to borehole breakage and associated fluid leakage. Further numerical, laboratory and field experiments are underway for additional technology sensitivity analysis, particularly the transient EM/Seismic reflectometry methods, data acquisition optimization, and numerical simulation improvements.

Wilt, Michael↗
DOE OSTI2022

Report on the LLNL Global Seismic Waveform Tomography Modeling Project

Earth models have important applications for seismic event monitoring including the location and characterization of potential underground nuclear tests. Seismic waves generated by earthquakes and man-made events are altered by variations in Earth’s subsurface properties; and these distortions effect the arrival times of energy packets (phases) and obscure the characteristics of the original energy source (event mechanism). Three-dimensional models of Earth’s seismic properties can help explain and predict the distortions to the seismic wavefield and reveal the properties of a source that generated the waves, including the event location using the model-predicted timing of the phases and source mechanism using the modelpredicted waveform characteristics.

58 GEOSCIENCES↗
NASA NTRS2019

Estimation of the Seismic Moment Release Rate of Mars from InSight Seismic Data

Seismicity models for Mars usually estimate the long-term average annual seismic moment rate, and also the average annual event rate. This holds for estimations based on geological evidence (Golombek et al., 1992, Golombek, 2002, Taylor et al., 2013) as well as for models based on thermal evolution and cooling of the Martian interior (Phillips, 1991, Knapmeyer et al., 2006, Plesa et al., 2018). All studies are compatible with the conclusion based on the non-observation of any unambiguous event by Viking (Anderson et al., 1977, Goins & Lazarewicz, 1979) that Martian seismicity lies somewhere between that of the Moon and that of the Earth. We developed tools to derive reasonable estimations of the annual seismic moment rate from a number of events as small as one, provided that the observed events are beyond the global completeness threshold for observable events. Numerical tests as well as evaluation of terrestrial data shows the feasibility of the approach.

Knapmeyer, Martin↗
DOE OSTI2023

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗