Search NASA⌕ Search

SEARCH · Search NASA

Results for “Velocity Model”

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

Seismic Waveform Inversion Capability on Resource-Constrained Edge Devices

Seismic full wave inversion (FWI) is a widely used non-linear seismic imaging method used to reconstruct subsurface velocity images, however it is time consuming, has high computational cost and depend heavily on human interaction. Recently, deep learning has accelerated it’s use in several data-driven techniques, however most deep learning techniques suffer from overfitting and stability issues. In this work, we propose an edge computing-based data-driven inversion technique based on supervised deep convolutional neural network to accurately reconstruct the subsurface velocities. Deep learning based data-driven technique depends mostly on bulk data training. In this work, we train our deep convolutional neural network (DCN) (UNet and InversionNet) on the raw seismic data and their corresponding velocity models during the training phase to learn the non-linear mapping between the seismic data and velocity models. The trained network is then used to estimate the velocity models from new input seismic data during the prediction phase. The prediction phase is performed on a resource-constrained edge device such as Raspberry Pi. Raspberry Pi provides real-time and on-device computational power to execute the inference process. In addition, we demonstrate robustness of our models to perform inversion in the presence on noise by performing both noise-aware and no-noise training and feeding the resulting trained models with noise at different signal-to-noise (SNR) ratio values. We make great efforts to achieve very feasible inference times on the Raspberry Pi for both models. Specifically, the inference times per prediction for UNet and InversionNet models on Raspberry Pi were 22 and 4 s respectively whilst inference times for both models on the GPU were 2 and 18 s which are very comparable. Finally, we have designed a user-friendly interactive graphical user interface (GUI) to automate the model execution and inversion process on the Raspberry Pi.

Manu, Daniel (ORCID:0000000154982677)↗

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗

No-net-rotation model of current plate velocities incorporating plate motion model NUVEL-1

NNR-NUVEL1 is presented which is a model of plate velocities relative to the unique reference frame defined by requiring no-net-rotation of the lithosphere while constraining relative plate velocities to equal those in global plate motion model NUVEL-1 (DeMets et al., 1990). In NNR-NUVEL1, the Pacific plate rotates in a right-handed sense relative to the no-net-rotation reference frame at 0.67 deg/m.y. about 63 deg S, 107 deg E. At Hawaii the Pacific plate moves relative to the no-net-rotation reference frame at 70 mm/yr, which is 25 mm/yr slower than the Pacific plate moves relative to the hotspots. Differences between NNR-NUVEL1 and HS2-NUVEL1 are described. The no-net-rotation reference frame differs significantly from the hotspot reference frame. If the difference between reference frames is caused by motion of the hotspots relative to a mean-mantle reference frame, then hotspots beneath the Pacific plate move with coherent motion towards the east-southeast. Alternatively, the difference between reference frames can show that the uniform drag, no-net-torque reference frame, which is kinematically equivalent to the no-net-rotation reference frame, is based on a dynamically incorrect premise.

Argus, Donald F.↗

Probability density and exceedance rate functions of locally Gaussian turbulence

A locally Gaussian model of turbulence velocities is postulated which consists of the superposition of a slowly varying strictly Gaussian component representing slow temporal changes in the mean wind speed and a more rapidly varying locally Gaussian turbulence component possessing a temporally fluctuating local variance. Series expansions of the probability density and exceedance rate functions of the turbulence velocity model, based on Taylor's series, are derived. Comparisons of the resulting two-term approximations with measured probability density and exceedance rate functions of atmospheric turbulence velocity records show encouraging agreement, thereby confirming the consistency of the measured records with the locally Gaussian model. Explicit formulas are derived for computing all required expansion coefficients from measured turbulence records.

Mark, W. D.↗

Resolution of Lava Tubes with Ground Penetrating Radar: The TubeX project

Remote sensing surveys of the Moon and Mars show evidence of lava tubes, which are potential safe havens for human crews and their equipment. Ground penetrating radar (GPR) can be used to map tubes because the void/rock interface at tube ceilings and floors strongly reflects radar pulses. We have tested the capacity of GPR to sense lava tube geometry at Lava Beds National Monument (LBNM) in California, USA. GPR and detailed LiDAR data are presented for two tubes: Skull Cave, with a few meters of overburden, diameter ~10-20 m, and a rubbly floor; and Valentine Cave, with similarly thin overburden, diameter ~1-3 meters, and a flatter smoother floor. On both caves GPR clearly resolves the ceiling and permits good estimates of the cave width as validated with LiDAR data. Where GPR fails, the primary cause is inferred to be strong out-of-plane effects due to complex 3D geometries. Recovery of the floor position requires migrating the GPR data with a 2D velocity model, as signal velocity is faster in void space. We find that floor position is recoverable in caves whose voids are taller than the radar wavelength (~3 m in this study). Forward modeling assuming planetary parameters suggests GPR should be similarly successful on the Moon or Mars.

S. Esmaeili↗

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗

Mantle Plume Magmatism in Elysium Planitia as Constrained by InSight Seismic Observations

Based on cratering statistics, there has been active volcanism in the last 100 Ma and possibly the last 1 to 10 Ma in both central Elysium Planitia and in the Cerberus Fossae region of south-eastern Elysium. Seismic measurements by the InSight mission show a strong concentration of seismicity in Cerberus Fossae which are interpreted as due to magma transport along dike systems. Together, these observations strongly support the existence of geologically recent and possibly of currently active magmatism in Elysium. InSight seismic observations were also used to construct seismic velocity models for the martian mantle. These models are effectively a regional aver-age of the seismic velocity structure for the Elysium region, where both the InSight lander and most of the measured seismic events are located. These seismic velocity models have been interpreted in terms of the lithospheric thickness and mantle potential temperature. Khan et al. found a lithospheric thickness of 400 to 600 km and a mantle potential temperature of 1325 to 1425 °C. Durán et al. found a lithospheric thickness of ~450 km and a mantle potential tempera-ture of 1375-1475 °C. Drilleau et al. found a thermal lithosphere of 420-660 km and a mantle potential temperature of 1380-1560 °C. In this study, we combine these InSight seismic observations with finite element mantle plume magma production models to place new constraints on the properties of the mantle beneath Elysium Planitia.

Walter S Kiefer↗

Multivalued velocity field model of turbulence

The considered study is based on new theoretical concepts regarding a post-instability model of a fluid discussed by Zak (1980). The model permits the completion of the governing equations of turbulence by introducing multivalued fields of velocities. Attention is given to the mechanism of energy dissipation, the characteristic wave propagation, a simplified model, the formation of turbulence around stagnation points, the formulation of boundary conditions, and the mechanism of turbulence formation. The mechanism of turbulence formation can be understood as propagation of initial discontinuities from the boundaries into a flow with the characteristic velocity which is defined by the normal (to the boundary) velocity components. These components emerge at the boundary as a result of jumps in the tangential components due to the continuity equation.

Zak, M.↗

Seismic identification of geothermal prospecting in Harrat Rahat, Northern Arabian Shield

Harrat Rahat volcanic field lies along the Makkah-Madinah-Nafud active volcanic line. It has extensive linear weak zones of wrench faults that could provide ideal conduits for upwelling magma and hydrothermal fluids as well. Here, passive seismic data recorded by seismic stations installed within and around Harrat Rahat has been processed for the detection of microseismic events. Moreover, ambient noise was also gathered and analyzed to construct a subsurface velocity model. This model indicates low velocity zone in Harrat Rahat, and two areas with coincident basement faults and possibly with enhanced subsurface temperatures were identified. It is concluded that Harrat Rahat has high geothermal potential and promising for geothermal renewable power.

58 GEOSCIENCES↗

Velocity structure and evolution of the moon

Seismic data from the Apollo Passive Seismic Network stations are analyzed to determine the velocity structure and to infer the composition and physical properties of the lunar interior. Data from artificial impacts (SIBV booster and LM-ascent stage) cover a distance range of 9 to 1750 km. Travel times and amplitudes, as well as theoretical seismograms, are used to derive a velocity model for the outer 150 km of the moon. The P-wave velocity model confirms an earlier report of a lunar crust in the eastern part of Oceanus Procellarum. The crust is about 60 km thick and may consist of two layers in the mare regions. Possible values for the P-wave velocity in the uppermost mantle are between 7.6 and 9.0 km/sec. The 9 km/sec velocity represents either a localized heterogeneous unit, or a thin layer less than about 40 km in thickness. The elastic properties of the deep interior, as inferred from the seismograms of natural events (meteoroid impacts and moonquakes) occurring at great distances, indicate that there is an increase in attenuation and a possible decrease of velocity at depths below about 1000 km.

Toksoz, M. N.↗

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗

Compensating attenuation effects in full-waveform inversion with dissipation-dispersion decoupling

Seismic attenuation poses challenges to velocity model building from acoustic/elastic full-waveform inversion (FWI). In particular, when constructing the FWI gradient for velocity inversion, recent studies have indicated that the high-attenuation structure could distort the gradient by damping the amplitudes and shifting the kinematic phases, resulting in an imbalanced update and thus unreliable velocity model. These Q effects are particularly significant in the reflection acquisition geometry due to the “double-damping” issue. Here, we develop a Q-compensated FWI algorithm for constructing a Q-free FWI gradient. By using a recently developed viscoacoustic wave propagator, this compensation can be done conveniently by keeping the dispersion term and flipping the dissipation term in the wave equation when we simulate the forward and time-reversed adjoint wavefields. The resultant gradient obtained by interacting these two wavefields has correct kinematics and Q-free amplitudes. This Q-compensated FWI can balance the update between Q- and no- Q-areas, which we determine using synthetic examples. In addition, we illustrate how to take advantage of the dissipation-dispersion decoupling to determine the anomalous Q value in the Q model building workflow via Q-compensated reverse time migration.

Geochemistry & Geophysics↗

Shear wave velocities in the earth's mantle.

Direct measurement of the travel time gradient for S waves together with travel time data are used to derive a shear velocity model for the earth's mantle. In order to satisfy the data it is necessary to discard the usual assumption of lateral homogeneity below shallow depths. A shear velocity differential is proposed for a region between western North America and areas of the Pacific Ocean. Distinctive features of the velocity model for the upper mantle beneath western North America are a low-velocity zone centered at 100 km depth and zones of high velocity gradient beginning at 400, 650, and 900 km.

Robinson, R.↗

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↗

Preliminary estimates of Gulf Stream characteristics from TOPEX data and a precise gravimetric geoid

TOPEX sea surface height data has been used, with a gravimetric geoid, to calculate sea surface topography across the Gulf Stream. This topography was initially computed for nine tracks on cycles 21 to 29. Due to inaccurate geoid undulations on one track, results for eight tracks are reported. The sea surface topography estimates were used to calculate parameters that describe Gulf Stream characteristics from two models of the Gulf Stream. One model was based on a Gaussian representation of the velocity while the other was a hyperbolic representation of velocity or the sea surface topography. The parameters of the Gaussian velocity model fit were a width parameter, a maximum velocity value, and the location of the maximum velocity. The parameters of the hyperbolic sea surface topography model were the width, the height jump, position, and sea surface topography at the center of the stream. Both models were used for the eight tracks and nine cycles studied. Comparisons were made between the width parameters, the maximum velocities, and the height jumps. Some of the parameter estimates were found to be highly (0.9) correlated when the hyperbolic sea surface topography fit was carried out, but such correlations were reduced for either the Gaussian velocity fits or the hyperbolic velocity model fit. A comparison of the parameters derived from 1-year TOPEX data showed good agreement with values derived by Kelly (1991) using 2.5 years of Geosat data near 38 deg N, 66 deg W longitude. Accuracy of the geoid undulations used in the calculations was of order of +/- 16 cm with the accuracy of a geoid undulation difference equal to +/- 15 cm over a 100-km line in areas with good terrestrial data coverage. This paper demonstrates that our knowledge or geoid undulations and undulation differences, in a portion of the Gulf Stream region, is sufficiently accurate to determine characteristics of the jet when used with TOPEX altimeter data. The method used here has not been shown to be more accurate than methods that average altimeter data to form a reference surface used in analysis to obtain the Gulf Stream characteristics. However, the results show the geoid approach may be used in areas where lack of current meandering reduces the accuracy of the average surface procedure.

Rapp, Richard H.↗