SEARCH · Engineering Papers

Results for “velocity model building”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗

DOE OSTI2023

SMART-CS Initiative: Developing Velocity and Saturation Models With Machine Learning

Objectives - Build velocity model: Develop a data driven approach for velocity prediction. Assess saturation models: Use ML to develop a data and a physical model-driven approach to estimate CO2 saturation (SCO2).

James, Maureen↗

DOE OSTI2022

Image-domain DAS 3D VSP elastic transmission tomography

SUMMARY Full-wavefield elastic imaging of active-source seismic data acquired by downhole receivers commonly offers higher-resolution subsurface images in the vicinity of a borehole compared to conventional surface seismic data sets, which can lack higher-frequency wavefield components due to longer travel paths and increased attenuation. An increasingly used approach for downhole acquisition is vertical seismic profiling (VSP), which has become more attractive when coupled with distributed acoustic sensing (DAS) using optical fibres installed in wells. The main difficulty for generating high-quality images with full-wavefield imaging tools for DAS VSP data, though, is the need for an accurate velocity model. To build plausible velocity models using active-source DAS VSP data, we adopt a 3-D image-domain elastic transmission tomography technique, originally developed for surface-recorded passive (microseismic) data, by exchanging the source and receiver positions (i.e. reciprocity) to mimic a passive-seismic surface monitoring scenario. The inversion approach exploits various images for each source constructed through time-reverse imaging (TRI) of downgoing P- and S-wave first-arrival waveforms. The TRI process uses the kinetic term of the (extended) PS energy imaging condition that exhibits sufficient sensitivity to velocity model errors. The method automatically updates the P- and S-wave velocity models to optimize image focusing via adjoint-state inversion. We illustrate the efficacy of the adopted elastic inversion technique using an active-source DAS 3-D VSP field data set acquired in the North Slope of Alaska. The numerical experiments demonstrate that the inverted elastic velocity models can be further used in full-wavefield acoustic/elastic imaging algorithms to obtain accurate subsurface images.

Oren, Can (ORCID:000000019534265X)↗

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↗

DOE OSTI2023

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↗

DOE OSTI2021

Progressive transfer learning for low-frequency data prediction in full waveform inversion

To effectively overcome the cycle-skipping issue in full waveform inversion (FWI), we developed a deep neural network (DNN) approach to predict the absent low-frequency components by exploiting the hidden physical relation connecting the low- and the high-frequency data. To efficiently solve this challenging nonlinear regression problem, two novel strategies were proposed to design the DNN architecture and to optimize the learning process: (1) dual data feed structure; (2) progressive transfer learning. With the dual data feed structure, not only the high-frequency data, but also the corresponding beat tone data are fed into the DNN to relieve the burden of feature extraction. The second strategy, progressive transfer learning, enables us to train the DNN using a single evolving training dataset. Within the framework of the progressive transfer learning, the training dataset continuously evolves in an iterative manner by gradually retrieving the subsurface information through the physics-based inversion module, progressively enhancing the prediction accuracy of the DNN and propelling the inversion process out of the local minima. Here, the synthetic numerical experiments suggest that, without any a priori geological information, the low-frequency data predicted by the progressive transfer learning are sufficiently accurate for an FWI engine to produce reliable subsurface velocity models free of cycle-skipping artifacts.

02 PETROLEUM↗

DOE OSTI2021

A physics-informed data-driven low order model for the wind velocity deficit at the wake of isolated buildings

Detailed simulations of the Atmospheric boundary layer (ABL) provides information on wind velocity, direction and turbulence for optimizing the design of wind turbines and wind farms and their operations. Distributed wind turbines have been considered for use in urban or suburban areas for covering energy needs using a sustainable energy solution. Simulations of high Reynolds number turbulent ABL flows past obstacles can be used for the prediction of the velocity deficit and turbulence characteristics in the leeward side of building-like obstacles, allowing for the estimation of power losses and fatigue loads in wind turbines. However, fully three-dimensional simulations of high Reynolds number ABL flows using Large Eddy Simulation (LES) or Reynolds Average Navier-Stokes Simulations (RANS) continue to require large computational/ execution times. As a result, these 3D simulations may not be ideal tools for the rapid assessment of wind energy potential. In the present study, 3D numerical results from wall resolved RANS simulations were used for the development of a low order model that allows the rapid assessment of wind velocity and turbulence in the leeward side of an isolated building.

17 WIND ENERGY↗

DOE OSTI2023

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↗

DOE OSTI2024

Inverse-Velocity Transformation Wall Model for Reacting Turbulent Hypersonic Boundary Layers

The present study builds on prior work by taking advantage of the novel framework proposed by Griffin et al. (hereafter referred to as the GFM) as a baseline. The model is progressively extended to multicomponent reacting mixtures, accounting for differential diffusion and finite-rate chemistry in a similar fashion to Di Renzo & Urzay and Di Renzo et al. The accuracy of the present approach, as well the prior model of Di Renzo & Urzay, is assessed in an a priori sense for the first time in a turbulent reacting boundary layer, using the boundary-layer data of Williams et al. Five species are included in the present analysis, i.e., N s = 5, namely N 2 , O 2 , NO, N and O, a neutral mixture most representative of dissociation/recombination phenomena for temperatures below 6000 K. A schematic of a flow over a wedge representative of the described configuration is presented in Figure 1. The brief is organized as follows: In Section 2, the wall-model equations and the computational framework are presented. In Section 3, the a priori results of the proposed model are described and compared to the extended EWM. Finally, in Section 4, some conclusions are offered.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A dynamic model of wind turbine yaw for active farm control

This paper presents a graph-based dynamic yaw model to predict the dynamic response of the hub-height velocities and the power of a wind farm to a change in yaw. The model builds on previous work where the turbines define the nodes of the graph and the edges represent the interactions between turbines. Advances associated with the dynamic yaw model include a novel analytical description of the deformation of wind turbine wakes under yaw to represent the velocity deficits and a more accurate representation of the interturbine travel time of wakes. The accuracy of the model is improved by coupling it with time- and space-dependent estimates of the wind farm inflow based on real-time data from the wind farm. The model is validated both statically and dynamically using large-eddy simulations. An application of the model is presented that incorporates the model into an optimal control loop to control the farm power output.

17 WIND ENERGY↗

DOE OSTI2021

Investigation of main bearing operating conditions in a three-Point mount wind turbine drivetrain

Wear-related failures of spherical roller bearings in the main bearing position of three-point mount wind turbines have been higher than expected and can contribute to higher than anticipated operation-and-maintenance costs. In this paper, the operational conditions of such a main bearing—including measured axial displacement and velocity subject to the estimated axial loads—are described for an instrumented, commercial wind turbine. The field measurements suggest a maximum axial speed between the bearing rings that is less than 2 millimeters per second. It is estimated that the axial speed between the ring and rollers is approximately 25% of this value. When compared to a speed of rolling from approximately 200–352 millimeters per second, the measured axial sliding is therefore significantly less than 1% of the speed of rolling. Previous numerical studies of lubricant film formation in rolling contacts have shown that the effect of axial sliding starts to be noticeable only when this ratio exceeds 10%; therefore, the axial velocity represents only a small disturbance to the nominal pure rolling case, and the influence on oil film building can be neglected. A simple analytic model of the main bearing motion was also developed and demonstrated similar displacement and velocity characteristics.

17 WIND ENERGY↗

DOE OSTI2021

Investigation of Main Bearing Operating Conditions in a Three-Point Mount Wind Turbine Drivetrain

Wear related failures of spherical roller bearings in the main bearing position of three-point mount wind turbines have been higher than expected and can contribute to higher than anticipated operation-and-maintenance costs. In this paper, the operational conditions of such a main bearing, including measured axial displacement and velocity subject to the estimated axial loads are described for an instrumented, commercial wind turbine. The field measurements suggest a maximum axial speed between the bearing rings that is less than 2 millimeters per second. It is estimated that the axial speed between the ring and rollers is approximately 25% of this value. When compared to a speed of rolling from approximately 200-352 millimeters per second, the measured axial sliding is therefore significantly less than 1% of the speed of rolling. Previous numerical studies of lubricant film formation in rolling contacts have shown that the effect of axial sliding starts to be noticeable only when this ratio exceeds 10%; therefore, the axial velocity represents only a small disturbance to the nominal pure rolling case, and the influence on oil film building can be neglected. A simple analytic model of the main bearing motion was also developed and demonstrated similar displacement and velocity characteristics.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

DOE OSTI2024

Reviving MeV-GeV indirect detection with inelastic dark matter

Thermal relic dark matter below ∼ 10 GeV is excluded by cosmic microwave background data if its annihilation to visible particles is unsuppressed near the epoch of recombination. Usual model-building measures to avoid this bound involve kinematically suppressing the annihilation rate in the low-velocity limit, thereby yielding dim prospects for indirect detection signatures at late times. In this work, we investigate a class of cosmologically viable sub-GeV thermal relics with late-time annihilation rates that are detectable with existing and proposed telescopes across a wide range of parameter space. We study a representative model of inelastic dark matter featuring a stable state χ 1 and a slightly heavier excited state χ 2 whose abundance is thermally depleted before recombination. Since the kinetic energy of dark matter in the Milky Way is much larger than it is during recombination, χ 1 χ 1 → χ 2 χ 2 upscattering can efficiently regenerate a cosmologically long-lived Galactic population of χ 2 , whose subsequent coannihilations with χ 1 give rise to observable gamma-rays in the ∼ 1 MeV − 100 MeV energy range. We find that proposed MeV gamma-ray telescopes, such as e-ASTROGAM, AMEGO, and MAST, would be sensitive to much of the thermal relic parameter space in this class of models and thereby enable both discovery and model discrimination in the event of a signal at accelerator or direct detection experiments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Qualification of ANSI/HPS N13.1-2011 Mixing Criteria by Computational Fluid Dynamics Modeling for the 3430 Building Fan Addition and Increased Ventilation Capacity

Additional ventilation capacity has been designed for the 3430 Building filtered exhaust stack system. The updated system will increase the number of fans from two to three and include new ductwork with a larger diameter to integrate the new fan into the existing stack. Stack operations will involve running various fan combinations at any given time. The air monitoring system of the existing two-fan stack previously was found to comply with the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard. Full-scale, three-dimensional computational fluid dynamics (CFD) modeling was used to evaluate the modified three-fan system for compliance with the ANSI/HPS N13.1-2011 standard, which essentially is equivalent to the ANSI/HPS N13.1-1999 standard. The four mixing criteria evaluated are 1) flow angle, 2) velocity, 3) gas tracer, and 4) particle tracer. Benchmarking of the CFD modeling methodology showed good agreement with previous testing used to qualify the stack, and modeling of the existing two-fan system showed good agreement with test data collected from the 3430 Building stack. Modeling was performed to develop a suitable three-fan design. Initial modeling of the three-fan design and basic ductwork showed that flow angles and velocity uniformity were acceptable; however, the gas tracer and particle tracer mixing results were not acceptable. To meet ANSI/HPS N13.1-2011 criteria, an air blender was added to the stack design.

42 ENGINEERING↗

DOE OSTI2022

Qualification of ANSI/HPS N13.1-2011 Mixing Criteria by Computational Fluid Dynamics Modeling for the 3430 Building Fan Addition and Increased Ventilation Capacity

42 ENGINEERING↗

DOE OSTI2024

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↗

DOE OSTI2020

Assessment of the 3420 Building Filtered Exhaust Stack Sampling Probe Location: Stack Verification Following Fan and Air Blender Additions

The Pacific Northwest National Laboratory (PNNL) 3420 Building, which is part of the Physical Sciences Facility (PSF), houses radiological capabilities that results in the requirement that emissions monitoring must be conducted for potential radionuclides in the exhaust air discharge of this building. The air monitoring system is required to conform to Title 40 of the Code of Federal Regulations part 61 (40 CFR 61) Subpart H, which in turns requires a sampling probe in the exhaust stream to conform to the criteria of American National Standards Institute / Health Physics Society (ANSI/HPS) N13.1-2011, Sampling and Monitoring Releases of Airborne Radioactive Substances from the Stack and Ducts of Nuclear Facilities. To support the air emissions permit for the 3420 Building on the PNNL campus, stack testing that used computational fluid dynamics (CFD) modeling as a surrogate stack and verification tests of velocity uniformity and flow angle on the retrofitted facility stack was performed. The ANSI/HPS N13.1-2011 criteria for the air monitoring probe location are that velocity uniformity, gaseous tracer uniformity, and particulate tracer uniformity must be less than or equal to 20%COV. Furthermore, no point in the sampling location may have a gaseous tracer concentration that varies from the mean concentration by more than 30%. Additionally, the flow angle at the sampling location must not be more than 20°. The CFD modeling of the stack, as reported by Recknagle et al. (2018) demonstrated that the stack meets the criteria at the probe location. The velocity uniformity and flow angle results from the 3420 stack verification tests, performed in October 2020, demonstrated that the CFD model results may be used to support the qualification of the stack sampling location. The measured velocity uniformity verification test result was 1.4%COV. This value is well within the uniformity criterion, which is that the velocity uniformity be =20%COV. Additionally, this value is well within the criterion that the actual stack measurement must be within 5% of the surrogate stack result of 2.1%COV when all four fans were operating. Additionally, the measured average flow angle at the 3420 stack monitor location was 15.5 degrees. Although this is higher than expected based on the CFD model, the result is =20 degrees, so the criterion is met. Based on these stack verification test results, the reconfigured 3420 Building filtered exhaust stack meets the qualification criteria given in the ANSI/HPS N13.1-2011 standard. Further changes to the system configuration or operating conditions that are outside the bounds described in this and the CFD report (Recknagle et al., 2018) may require additional tests and additional analysis to determine compliance with the standard.

3420 Building↗

DOE OSTI2023

Assessment of the 3430 Building Filtered Exhaust Stack Sampling Probe Location: Stack Verification Following Fan and Air Blender Additions

The velocity uniformity and flow angle results from the 3430 stack verification tests, performed in April 2023, demonstrated that the CFD model results may be used to support the qualification of the stack sampling location. The measured velocity uniformity verification test result was 2.1 %COV. This value is well within the uniformity criterion, which is that the velocity uniformity be =20 %COV. Additionally, this value is well within the criterion that the actual stack measurement must be within 5% of the surrogate stack (i.e., CFD modeled stack); in this case the CFD modeled average result of 2.85 %COV for the nominal operating range of 22,800 cfm to 62,400 cfm. Additionally, the measured average flow angle at the 3430 stack monitor location was 5.6 degrees. The result is =20 degrees, so the criterion is met.

3430 Building↗