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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗
NASA NTRS1972

Q and structure.

Different values of seismic Q may be obtained from measurements of different portions of a seismic record, indicating a separation of the effects on energy propagation of scattering and real loss parameters. In cases in which these parameters can be separated, they can then be used in seismic modeling. As an example, a model is presented in which scattering occurs to a depth of about 10 km below the volcanic ridge in the Tonga-Fiji region. The case of the lunar maria is discussed in terms of the relation of measured seismic Q, scattering parameters, and direct ray propagation parameters. The results indicate that a fairly simple jointed bedrock model is compatible with observed lunar seismic data.

Warren, N.↗
DOE OSTI2020

Composition of Magma and Characteristics of the Hydrothermal System of Newberry Volcano, Oregon, From Magnetotellurics

Abstract We use 3‐D magnetotellurics to improve our understanding of the structure and magma composition of Newberry Volcano in Oregon, USA. Newberry is a broad shield volcano with a summit caldera and strongly bimodal magmatism. Newberry has long been the subject of geothermal exploration research, but that work has focused on the volcano's west flank, leaving the caldera largely unstudied with geophysical methods until recently. Our modeling shows a relatively resistive magma reservoir of approximately 50 Ωm. Our work builds upon recent seismic models and petrological analysis to interpret Newberry's magma reservoir as a dry rhyolite with 8–11% partial melt, which matches the seismically determined melt fraction. Finding the conditions within the magma reservoir that allow the resistivity and seismic analysis to come to the same melt fraction helps us narrow down the magma temperature and composition. From this we infer a dry rhyolitic magma at 850 °C. We also image a prominent vertical conductive anomaly along the south rim below the vent that produced the most recent eruption. The anomaly extends from magma reservoir depths of 3 km to 1 km below the surface where it disperses into the caldera fill. We interpret this as the main conduit for magmatic fluids to reach Newberry's hydrothermal system. Other features that our model shows include higher conductivity along the caldera rim and a resistive older pluton on the west flank.

Bowles‐Martinez, E.↗
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2011

Further Constraints and Uncertainties on the Deep Seismic Structure of the Moon

The Apollo Passive Seismic Experiment (APSE) consisted of four 3-component seismometers deployed between 1969 and 1972, that continuously recorded lunar ground motion until late 1977. The APSE data provide a unique opportunity for investigating the interior of a planet other than Earth, generating the most direct constraints on the elastic structure, and hence the thermal and compositional evolution of the Moon. Owing to the lack of far side moonquakes, past seismic models of the lunar interior were unable to constrain the lowermost 500 km of the interior. Recently, array methodologies aimed at detecting deep lunar seismic reflections found evidence for a lunar core, providing an elastic model of the deepest lunar interior consistent with geodetic parameters. Here we study the uncertainties in these models associated with the double array stacking of deep moonquakes for imaging deep reflectors in the Moon. We investigate the dependency of the array stacking results on a suite of parameters, including amplitude normalization assumptions, polarization filters, assumed velocity structure, and seismic phases that interfere with our desired target phases. These efforts are facilitated by the generation of synthetic seismograms at high frequencies (approx. 1Hz), allowing us to directly study the trade-offs between different parameters. We also investigate expected amplitudes of deep reflections relative to direct P and S arrivals, including predictions from arbitrarily oriented focal mechanisms in our synthetics. Results from separate versus combined station stacking help to establish the robustness of stacks. Synthetics for every path geometry of data were processed identically to that done with data. Different experiments were aimed at examining various processing assumptions, such as adding random noise to synthetics and mixing 3 components to some degree. The principal stacked energy peaks put forth in recent work persist, but their amplitude (which maps into reflector impedance contrast) and timing (which maps into reflector depth) depend on factors that are not well constrained -- most notably, the velocity structure of the overlying lunar interior. Thus, while evidence for the lunar core remains strong, the depths of imaged reflectors have associated uncertainties that will require new seismic data and observations to constrain. These results strongly advocate further investigations on the Moon to better resolve the interior (e.g., Selene missions), for the Moon apparently has a rich history of construction and evolution that is inextricably tied to that of Earth.

Lin, Pei-Ying Patty↗
DOE OSTI2024

Report on the sensitivity kernel construction and updating the SPiRaL model with waveform data

The WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes research leading towards the prediction of entire seismic and acoustic waveforms produced by natural and manmade events, including explosions. A key aspect of this research is the development of Earth (seismic) models at multiple scales including crustal, regional, and global. The work described here pertains to the global-scale seismic tomography effort led by LLNL.

58 GEOSCIENCES↗
DOE OSTI2023

Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data

Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.

58 GEOSCIENCES↗
DOE OSTI2021

Seismic Visibility of Eclogite in the Earth’s Upper Mantle—Implications From High Pressure‐Temperature Single‐Crystal Elastic Properties of Omphacite

Abstract Identifying and locating the geochemical and geophysical heterogeneities in the Earth’s interior is one of the most important and challenging tasks for the deep Earth scientists. Subducted oceanic crust metamorphizes into the dense eclogite in the upper mantle and is considered as a major cause of geochemical and geophysical heterogeneities in the deep Earth. In order to detect eclogitic materials inside the Earth, precise measurements of the high pressure‐temperature single‐crystal elasticity of major minerals in eclogite are thus exceedingly important. Omphacite, a Na,Al‐bearing clinopyroxene, constitutes up to 75 vol% of eclogite. In the present study, we performed the first high pressure‐temperature single‐crystal elasticity measurements of omphacite using Brillouin spectroscopy. Utilizing the finite‐strain approach, we obtained the following thermoelastic parameters for omphacite: K S0 ’ = 4.5(1), G 0 ’ = 1.53(5), ∂K S0 /∂ T = −0.029(5) GPa/K, ∂G 0 /∂ T = −0.013(5) GPa/K, with K S0 = 123(3) GPa, G 0 = 74(2) GPa, and ρ 0 = 3.34(1) g/cm 3 . We found that the seismic velocities of undeformed eclogite are similar to pyrolite at the depths of 200–300 and 410–500 km, thus eclogite is seismically invisible at these depths. Combined with the lattice‐preferred orientations of the omphacite in naturally deformed eclogites, we also modeled seismic anisotropy of eclogite at various pressure‐temperature conditions. A 10 km thick subducted eclogitic crust can result in ∼0.2 s shear wave splitting in the Earth’s upper mantle.

Hao, Ming↗
DOE OSTI2025

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗
NASA NTRS1991

A scale-invariant cellular-automata model for distributed seismicity

In the standard cellular-automata model for a fault an element of stress is randomly added to a grid of boxes until a box has four elements, these are then redistributed to the adjacent boxes on the grid. The redistribution can result in one or more of these boxes having four or more elements in which case further redistributions are required. On the average added elements are lost from the edges of the grid. The model is modified so that the boxes have a scale-invariant distribution of sizes. The objective is to model a scale-invariant distribution of fault sizes. When a redistribution from a box occurs it is equivalent to a characteristic earthquake on the fault. A redistribution from a small box (a foreshock) can trigger an instability in a large box (the main shock). A redistribution from a large box always triggers many instabilities in the smaller boxes (aftershocks). The frequency-size statistics for both main shocks and aftershocks satisfy the Gutenberg-Richter relation with b = 0.835 for main shocks and b = 0.635 for aftershocks. Model foreshocks occur 28 percent of the time.

Barriere, Benoit↗
DOE OSTI2021

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗
DOE OSTI2025

Prediction of Complex Observed Shear Wave Splitting Patterns at Ryukyu Subduction Zone Using a Strong Intra-Slab Anisotropy Model

Complex shear wave splitting (SWS) patterns in subduction zones are often interpreted geodynamically as resulting from complex mantle flow; however, this may not always be necessary. We analyzed 7,093 high-quality SWS measurements from teleseismic S waves recorded by Hi-net stations across the Ryukyu arc in Japan. Our findings show a systematic rotation of the fast S polarization from trench-parallel to trench-perpendicular depending on the earthquake backazimuth. For the same earthquake, the measured splitting patterns also vary spatially across the southwest Japan. Using full-wave seismic modeling, we showed that a dipping slab with ~30% shear anisotropy of the tilted transverse isotropy (TTI) type, with a symmetry axis perpendicular to the slab interface, can predict the observed delay times and polarization rotation. Our results highlight the importance of considering dipping anisotropic slabs in interpreting SWS at subduction zones.

58 GEOSCIENCES↗
DOE OSTI2022

Report on Depth-Dependent Q from Frequency-Dependent Lg Q

As a phase that travels through the crustal waveguide, Lg samples the velocity and attenuation structure of the crust. We have developed methods (Pasyanos et al., 2009a; Pasyanos et al., 2009b) to estimate Qp and Qs of the crust and upper mantle from the amplitudes of regional phases, and have applied it to a number of regions around the world, and over a broad frequency band (0.5 - 10 Hz). In prior attenuation modeling, Lg phase amplitudes are able to effectively distinguish the high Q (low attenuation) of old stable cratonic crust from the low Q (high attenuation) of more recently active tectonic regions. For a variety of purposes (e.g. waveform modeling, seismic hazard, etc.), it is useful to have more information on the depth-dependence of the Q structure in the crust. One obvious candidate for doing this is using the frequency-dependence information of Lg Q, in the manner of surface waves. Unfortunately, however, Lg samples the whole of the crust without much differentiation with respect to frequency that is characteristic of fundamental mode surface waves. Lg is generally thought to be the sum of supercritically reflected S-waves trapped in the crustal waveguide (e.g. Xie and Lay, 1994). The phase can also be modeled as a superposition of surface waves. At lower frequencies, it is primarily composed of fundamental and low-order surface waves, while at higher frequencies, it is comprised purely of many high-order surface waves. The low-order surface waves that comprise Lg in low frequency bands (< 1 Hz) are preferentially sensitive to shallow earth structure, such as sedimentary basins and upper crust, while high-order surface waves have a more homogenized sensitivity to the whole crust. Furthermore, at the lower frequencies, we might expect to see some differentiation in the sensitivity between early Lg (herein referred to as Lg1) and late Lg (referred to as Lg2) where the latter is comprised of lower-order surface waves and, hence, shallower structure. Therefore, we might expect to see the signature of shallow features, primarily the lower Q of sediments as compared to crystalline crust, at lower frequencies and in later arriving Lg amplitudes.

58 GEOSCIENCES↗
DOE OSTI2023

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Strength, plasticity, and spin transition of Fe-N compounds in planetary cores

Elastic and plastic properties of Fe-light element alloys and compounds are needed to determine the compositions and dynamics of planetary cores. Elastic strength and plastic deformation mechanisms and their relationship to electronic properties of ε-Fe 7 N 3 and γ’-Fe 4 N mixture were investigated by x-ray diffraction and x-ray emission spectroscopy in the diamond anvil cell from 1 bar up to 60 GPa. X-ray diffraction shows that ε-Fe 7 N 3 reaches a pressure of 15–20 GPa before undergoing bulk plasticity at a differential stress of 4.4–10.4 GPa. ε-Fe 7 N 3 is stronger than γ’-Fe 4 N and hcp-Fe which achieve a flow stress of 1.5–3.6 GPa at 10–15 GPa and 2–3 GPa at ~20 GPa, respectively. X-ray emission spectroscopy shows that a decrease in electronic spin moment begins before and completes after plastic flow onset for each nitride, suggesting that pressure-driven changes in electronic arrangement do not trigger a plastic response although they may modify the strength and plastic behavior of Fe-N compounds. Plastic deformation in ε-Fe 7 N 3 and hcp-Fe results in a preferred orientation of (0001) normal to maximum compression, while γ’-Fe 4 N develops a maximum in the (110). Furthermore, these observations may be combined with measurements of elasticity to model seismic properties of cores of small planetary bodies such as Mars, Mercury, and the Moon.

58 GEOSCIENCES↗
DOE OSTI2023

Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai

Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga-Hunga Ha'apai eruption, four remarkably similar seismic subevents within a 5-min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (~25 s each). Each subevent begins with an unusual negative P-wave polarity which is inferred, using full-wave seismic modeling, to be caused by an upward single-force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights.

58 GEOSCIENCES↗
DOE OSTI2023

Data-driven time-scale separation of ODE right-hand sides using dynamic mode decomposition and time delay embedding

Multi-physics simulation often involves multiple different scales. The ARKODE ODE solver package in the SUNDIALS library addresses multi-scale problems with a multi-rate time-integrator that can work with a right-hand side that has fast scale and slow scale components. In this report, we use dynamic mode decomposition and time delay embedding to extract the fast and slow components of the right-hand sides of a simple ODE from data. We then use the extracted components to solve the ODE with ARKODE. Finally, to move towards a real-world use case, we attempt to extract fast and slow scale dynamics from synthetic seismic modeling data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗