Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic State Estimation”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of Power-Maximizing Region 2 Controllers for Wind and Marine Turbines

Wind and marine energy are rapidly growing and complementary technologies that share some techniques for simplified modeling and control, particularly in below-rated flow speeds. A turbine operator has several choices of controller for maximizing power in Region 2. The simple and ubiquitous KΩ 2 control law is often effective but limited in its flexibility. Alternative controllers use reference tracking to split the control objectives into a low-bandwidth optimal tip-speed ratio tracking loop to maximize steady-state power and a higher-bandwidth proportional-integral control loop to reject inflow turbulence. Several options exist for identifying the slowly varying optimal set point during operation, based on estimating the inflow velocity or filtering the power or torque signals. This study compares the trade-offs between performance and other design priorities for a few choices of reference-tracking controller in the literature for reference wind and marine turbines. Analysis is performed in the frequency domain using the linearization of each controller, and the impact of turbulent disturbances on the closed-loop system is described. The controllers are simulated in OpenFAST to analyze their performance with higher-order nonlinear turbine dynamics.

17 WIND ENERGY↗
DOE OSTI2019

Statistically Rigorous Uncertainty Quantification for Physical Parameter Model Calibration with Functional Output

In experiments conducted on the Z-machine at Sandia National Laboratories, dynamic material properties cannot be analyzed using traditional analytic methods, necessitating solving an inverse problem. Bayesian model calibration is a statistical framework for solving an inverse problem to estimate parameters input into a computational model in the presence of multiple uncertainties. Disentangling input parameter uncertainty and model misspecification is often poorly identified problem. When using computational models for physical parameter estimation, the issue of parameter identifiability must be carefully considered to obtain accurate and precise estimates of physical parameters. Additionally, in dynamic material properties applications, the experimental output is a function, velocity over time. While we can sample an arbitrarily large number of points from the measured velocity, these curves only contain a finite amount of information about the calibration parameters. In this report, we propose modifications to the Bayesian model calibration framework to simplify and improve the estimation of physical parameters with functional outputs. Specifically, we propose scaling the likelihood function by an effective sample size rather than modeling the discrepancy function; and modularizing input nuisance parameters with weakly identified parameters. We evaluate the performance of these proposed methods using a statistical simulation study and then apply these methods to estimate parameters of the tantalum equation of state. We conclude that these proposed methods can provide simple, fast, and statistically valid alternatives to the full Bayesian model calibration procedure; and that these methods can be used to estimate parameters of the equation of state for tantalum.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Microstructural Changes to Thermally Sprayed Materials Subjected to Dynamic Compression

Dynamic compression of materials can induce a variety of microstructural changes. As thermally-sprayed materials have highly complex microstructures, the expected pressure at which changes occur cannot be predicted a priori. In addition, typical in-situ measurements such as velocimetry are unable to adequately diagnose microstructural changes such as failure or pore collapse. Quasi-isentropic compression experiments with sample recovery were conducted to examine microstructural changes in thermally sprayed tantalum and tantalum-niobium blends up to 8 GPa pressure. Spall fracture was observed in all tests, and post-shot pore volume decreased relative to the initial state. The blended material exhibited larger spall planes with fracture occurring at interphase boundaries. An estimate of the pressure at which pore collapse is complete was determined to be ~26 GPa for pure tantalum and ~19 GPa for the tantalumniobium blend under these loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Estimates of Decadal Climate Predictability From an Interactive Ensemble Model

Decadal climate predictability has received considerable scientific interest in recent years, yet the limits and mechanisms for decadal predictability are currently not well known. It is widely accepted that noise due to internal atmospheric dynamics at the air-sea interface influences predictability. The purpose of this paper is to use the interactive ensemble (IE) coupling strategy to quantify how internal atmospheric noise at the air-sea interface impacts decadal predictability. The IE technique can significantly reduce internal atmospheric noise and has proven useful in assessing seasonal-to-interannual variability and predictability. Here we focus on decadal timescales and apply the nonlinear local Lyapunov exponent method to the Community Climate System Model comparing control simulations with IE simulations. This is the first time the nonlinear local Lyapunov exponent has been applied to the state-of-the-art coupled models. The global patterns of decadal predictability are discussed from the perspective of internal atmospheric noise and ocean dynamics.

58 GEOSCIENCES↗
DOE OSTI2025

Energy impacts of nationwide window upgrades in commercial buildings

This report presents comprehensive estimates of the energy impacts of nationwide commercial building window upgrades in the United States, using a conservative approach. Windows play a substantial role in determining building energy use and occupant experience. Estimates point to commercial building windows impacting loads that represent more than 6 quads (approximately 6%) of annual primary energy use in the U.S. (Harris, 2022). Beyond heating and cooling loads, windows also have effects on lighting and occupant comfort. The fastest route to improving the energy efficiency of windows in U.S. buildings is upgrading or replacing windows in existing buildings. This is due to poor performance of windows in older existing buildings compared to most new construction, low levels of window replacement, and long window service life compared to energy-using building components. Nationwide window upgrades were considered using the following technologies: • Secondary glazing systems • Double pane (clear and tinted) • Triple pane (clear and tinted) • Electrochromic glazing Nationwide upgrades provide on the order of 4%–6% site energy savings in typical buildings, or up to 26% in buildings with the highest savings potential. Electrochromic windows, with their ability to adapt dynamically to environmental conditions, can provide additional benefits, ranging from median savings of 7.2% in buildings with window to wall ratio (WWR) greater than 10% and up to 28% for some buildings. Savings increase substantially for buildings with higher WWR. This study’s approach focused on isolating the direct energy benefits from improvement in window performance, and does not take into account the following additional benefits from window retrofits, which are likely to be substantial: • Managing peak demand and enabling HVAC equipment downsizing. • Energy savings from customizing upgrades to building type and climate. • Energy savings and comfort improvements resulting from post-retrofit reductions in air leakage. • Non-energy benefits, such as occupant comfort and resilience during extreme weather.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

First-Principles Approach for Assessing Cold Electron Injection Efficiency of Dye-Sensitized Solar Cell: Elucidation of Mechanism of Charge Injection and Recombination

The adequacy of the inclusion of spacer units in the metal-free D-π–A organic dyes pertaining to the augmentation of dye-sensitized solar cell (DSSC) efficiency has been examined through the excited-state simulations of the charge injection and recombination processes at the dye–semiconductor interface. Within the framework of the time-dependent density functional theory, the proposed computational studies focus on the precise evaluation of pivotal factors controlling the rates of photoinduced charge-transfer and energy-transfer processes, including electronic coupling, reorganization energy, and threshold energy barrier in the semiclassical Marcus formalism. The estimation of the fluorescent state appears to be the crucial step while explaining the ultrafast electron injection process and the charge recombination at the Marcus inverted region, as revealed by the obtained results. The retardation of charge recombination is facilitated by the insertion of a thiophene moiety between the π-bridge and the acceptor units. The estimated cold electron injection efficiencies deploying the Onsager–Braun theory, that rely on the computations of cold electron injection lifetime and cold electron lifetime, show a linear correlation with the experimental photovoltaic parameters of the DSSC comprising short-circuit current density, open-circuit voltage, and power conversion efficiency. The outcomes of the present investigation establish a basis for unraveling the mechanism of intricate dynamical processes upon photoexcitation of the sensitizers, as well as devising plausible routes for functional DSSC materials.

14 SOLAR ENERGY↗
DOE OSTI2025

Solving reaction dynamics with quantum computing algorithms

The description of quantum many-body dynamics is extremely challenging on classical computers, as it can involve many degrees of freedom. However, the time evolution of quantum states is a natural application for quantum computers that are designed to efficiently perform unitary transformations. Here, in this paper, we study quantum algorithms for response functions, relevant for describing different reactions governed by linear response. We focus on nuclear-physics applications and consider a qubit-efficient mapping on the lattice, which can efficiently represent the large volumes required for realistic scattering simulations. For the case of a contact interaction, we develop an algorithm for time evolution based on the Trotter approximation that scales logarithmically with the lattice size and is combined with quantum phase estimation. We eventually focus on the nuclear two-body system and a typical response function relevant for electron scattering as an example. We also investigate ground-state preparation and examine the total circuit depth required for a realistic calculation and the hardware noise level required to interpret the signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Steps to achieving high-resolution NMR spectroscopy on solutions at GPa pressure

Recent theoretical advances and molecular-dynamic estimates of the dielectric constant of water have extended the HKF model for aqueous solution chemistry up to 6.0 gigapascals (GPa), which are conditions well beyond the capabilities of conventional NMR spectroscopy (Pan et al., 2013; Sverjensky et al., 2014, see also Wasserman and others, 1995). These developments provide strong motivation to design a simple NMR probe that allows experiments on solutions, at high resolution, to pressures of a few GPa (Pautler et al., 2014; Ochoa et al., 2015; 2016). Here we describe the performance of a compact NMR probe that can reach several GPa pressures. The probe is made by placing a solenoid microcoil inside of a standard piston-cylinder device used in solid-state physics. High pressures are achieved in the sample by applying force to a coaxial piston. Early designs of the probe, although useful, were limited in sample size to ~10-15μL. Here we describe modifications that allow greatly improved resolution and sensitivity, including 1H-1H NMR correlation spectroscopy, on solutions at 2.8 GPa pressure. Sample sizes can be expanded if, instead of a standard NMR spectrometer that is built around a superconducting magnet, one employs a magnetic-resonance imaging (MRI) system that is built around a permanent magnet. The MRI systems can apply larger magnetic field gradients than conventional spectrometers, and thus have more robust shimming capability, which is needed because of the juxtaposition of different alloys in the pressure cell. More importantly, these systems can have magnetic fields that are oriented perpendicular to the magnet bore, which allows rotation of the coil within the axis of the pressure cell to increase sample volumes and pressure. Sensitivity is improved by replacing the traditional, but reduced-volume, closed sample container mounted in the center of the NMR detection coil with an open NMR coil that dangles freely in solution, thus making the sample solution itself the pressure-transmission fluid. The efficacy of these modifications are demonstrated by measuring 1 H NMR spectra for ethyl alcohol and ethyl alcohol/methyl alcohol mixtures at pressures up to 2.8 GPa. As a result, further developments are discussed that will allow geochemists to acquire aqueous solution NMR spectra at higher 3-4 GPa pressures.

58 GEOSCIENCES↗
DOE OSTI2025

Can Neutron Star Tidal Effects Obscure Deviations from General Relativity?

Abstract One of the main goals of gravitational-wave astrophysics is to study gravity in the strong-field regime and constrain deviations from general relativity (GR). Any such deviation affects not only binary dynamics and gravitational-wave emission but also the structure and tidal properties of compact objects. In the case of neutron stars, masses, radii, and tidal deformabilities can all differ significantly between different theories of gravity. Currently, the measurement uncertainties in neutron star radii and tidal deformabilities are quite large. However, much less is known about how the large uncertainty in the nuclear equation of state (EOS) might affect tests of GR using binary neutron star mergers. Conversely, using the wrong theory of gravity might lead to incorrect constraints on the nuclear EOS. Here, we study this problem within scalar–tensor (ST) theory. We apply the recently derived ℓ = 2 tidal Love numbers in this theory to parameter estimation of GW170817. Correspondingly, we test if physics beyond GR could bias measurements of the nuclear EOS and neutron star radii. We find that parameter inference for both the GR and ST cases returns consistent component masses and tidal deformabilities. The radius and the EOS posteriors, however, differ between the two theories, but neither is excluded by current observational limits. This indicates that measurements of the nuclear EOS may be biased and that deviations from GR could go undetected when analyzing current binary neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Vitrification of Hanford Tank Waste 241-AP-107 in a Continuous Laboratory-Scale Melter

Low-activity waste (LAW) stored in underground tanks on the Hanford Site in Washington State is planned to be filtered for solids removal and processed through ion exchange columns for cesium removal. These pretreatment steps will allow the waste to be transferred to the Hanford Tank Waste Treatment and Immobilization Plant’s LAW Facility for immobilization into glass. The liquid waste will be combined with glass-forming chemicals (GFCs) to form a waste feed slurry that can be fed to electric melters for vitrification. The process of continuously converting the aqueous feed slurry into a melt is dynamic and includes multiple reactions, degassing, and dissolution processes that depend on heat from the melt below. In this conversion process, waste components are partitioned into one of two streams: glass and off-gas. Washington River Protection Solutions has requested processing information and chemical information associated with these waste products for actual waste from Hanford tank 241-AP-107 (referred to herein as AP-107). To acquire this type of information, a small-scale melter system was designed that would not require high volumes of input waste or the large resource commitment of a full-scale melter system, while also providing dynamic information that would be difficult to determine from batch reactions in a crucible system. A continuous laboratory-scale melter (CLSM) has been designed to operate with a continuous feeding process while periodically pouring glass product and collecting off-gas. The CLSM vessel has been sized to collect the relevant process and chemical information from obtainable volumes of AP-107 waste samples. A total 8.6 L of actual AP-107 tank waste were received after filtering for solids removal and ion exchange for cesium removal. This volume of AP-107 waste was mixed with GFCs to form an estimated 11.8 L of melter feed slurry. A mass of 7.01 kg of glass product were poured from the CLSM vessel during 10.07 hours of charging the AP-107 melter feed, indicating that about 15.0 kg of melter feed was vitrified into glass. The CLSM system was also designed with the capability to fully divert the flow of off-gas produced in the CLSM vessel to a sampling line that could capture the volatile species of interest, such as technetium-99 ( 99 Tc). This novel sampling system avoided the difficulties of slipstream sampling and could be activated once the feeding reached a steady state within the CLSM vessel. Off-gas product samples captured via high-efficiency particulate air (HEPA) filters using this method, as well as selected glass product samples were sent to an analytical lab for chemical analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
DOE OSTI2026

Adaptive Sampling-Based Bi-Fidelity Stochastic Trust Region Method for Stochastic Derivative-Free Optimization

Bi-fidelity stochastic optimization has gained increasing attention as an efficient approach to reduce computational costs by leveraging a low-fidelity (LF) model to optimize an expensive high-fidelity (HF) objective. In this paper, we propose ASTRO-BFDF, an adaptive sampling trust-region method specifically designed for unconstrained bi-fidelity stochastic derivative-free optimization problems. In ASTRO-BFDF, the LF function serves two purposes: (i) to identify better iterates for the HF function when the optimization process indicates a high correlation between them and (ii) to reduce the variance of the HF function estimates using bi-fidelity Monte Carlo (BFMC). The algorithm dynamically determines sample sizes while adaptively choosing between crude Monte Carlo and BFMC to balance the trade-off between optimization and sampling errors. We prove that the iterates generated by ASTRO-BFDF converge to a first-order stationary point almost surely. Additionally, we demonstrate the effectiveness of the proposed algorithm through numerical experiments on synthetic benchmarks and simulation optimization problems involving discrete event systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Static and dynamic magnetic properties of the spiin-$\frac{5}{2}$triangle lattice antiferromagnet Na 3 Fe(PO 4 ) 2 studied by 31 P NMR

We report 31 P nuclear magnetic resonance (NMR) measurements have been carried out to investigate the magnetic properties and spin dynamics of Fe 3+ (S = 5/2) spins in the two-dimensional triangular lattice (TL) compound Na 3 Fe(PO 4 ) 2 . The temperature (T) dependence of nuclear spin-lattice relaxation rates (1/T 1 ) shows a clear peak around Néel temperature, T N = 10.9 K, corresponding to an antiferromagnetic (AFM) transition. From the temperature dependence of NMR shift (K) above T N , an exchange coupling between Fe 3+ spins was estimated to be J/k B ≃ 1.9 K using the spin-5/2 Heisenberg isotropic-TL model. The temperature dependence of 1/T 1 T divided by the magnetic susceptibility (χ), 1/T 1 T χ , above T N proves the AFM nature of spin fluctuations below ~50 K in the paramagnetic state. In the magnetically ordered state below T N , the characteristic rectangular shape of the NMR spectra is observed, indicative of a commensurate AFM state in its ground state. The strong temperature dependence of 1/T 1 in the AFM state is well explained by the two-magnon (Raman) process of the spin waves in a 3D antiferromagnet with a spin-anisotropy energy gap of 5.7 K. The temperature dependence of sublattice magnetization is also well reproduced by the spin waves. Those results indicate that the magnetically ordered state of Na 3 Fe(PO 4 ) 2 is a conventional 3D AFM state, and no obvious spin frustration effects were detected in its ground state

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Identified charged hadron production in Au+Au collisions at $\sqrt{s_{NN}}$ = 54.4 GeV with the STAR detector

Here, we present results on the production of 𝜋±, 𝐾±, 𝑝, and $\bar{𝑝}$ in Au+Au collisions at $\sqrt{𝑠_{N⁢N}}$ = 54.4~GeV using the STAR detector at RHIC, at mid-rapidity (|𝑦|< 0.1). Invariant yields of these particles as a function of transverse momentum are shown. We determine bulk properties such as integrated particle yields (𝑑 ⁢𝑁 /𝑑 ⁢𝑦 ), mean transverse momentum (⟨𝑝 𝑇 ⟩), particle ratios, which provide insight into the particle production mechanisms. Additionally, the kinetic freeze-out parameters (𝑇 kin and ⟨𝛽 𝑇 ⟩), which provide information about the dynamics of the system at the time of freeze-out, are obtained. The Bjorken energy density (𝜖 𝐵⁢𝐽 ), which gives an estimate of the energy density in the central rapidity region of the collision zone at the formation time 𝜏, is calculated and presented as a function of multiplicity for various energies. The results are compared with those from the models such as A Multi-Phase Transport (AMPT) and Heavy Ion Jet INteraction Generator (HIJING) for further insights.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Quantification of MagLIF Morphology using the Mallat Scattering Transformation

The morphology of the stagnated plasma resulting from Magnetized Liner Inertial Fusion (MagLIF) is measured by imaging the self-emission x-rays coming from the multi-keV plasma, and the evolution of the imploding liner is measured by radiographs. Equivalent diagnostic response can be derived from integrated rad-MHD simulations from programs such as Hydra and Gorgon. There have been only limited quantitative ways to compare the image morphology, that is the texture, of simulations and experiments. We have developed a metric of image morphology based on the Mallat Scattering Transformation (MST), a transformation that has proved to be effective at distinguishing textures, sounds, and written characters. This metric has demonstrated excellent performance in classifying ensembles of synthetic stagnation images. We use this metric to quantitatively compare simulations to experimental images, cross experimental images, and to estimate the parameters of the images with uncertainty via a linear regression of the synthetic images to the parameter used to generate them. This coordinate space has proved very adept at doing a sophisticated relative back-ground subtraction in the MST space. This was needed to compare the experimental self emission images to the rad-MHD simulation images. We have also developed theory that connects the transformation to the causal dynamics of physical systems. This has been done from the classical kinetic perspective and from the field theory perspective, where the MST is the generalized Green's function, or S-matrix of the field theory in the scale basis. From both perspectives the first order MST is the current state of the system, and the second order MST are the transition rates from one state to another. An efficient, GPU accelerated, Python implementation of the MST was developed. Future applications are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SPRUCE Root Production Assessed with Manual Minirhizotrons Resolved to Plant Functional Type, 2015-2021

This dataset contains raw root length and diameter for individual roots and estimated root population production measurements from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experimental site within the Marcell Experimental Forest in northern Minnesota, USA. Measurements started at the beginning of whole ecosystem warming manipulations in 2015 through 2021 (2015-05-26 to 2021-09-01). Root morphology and estimated production were quantified throughout the peat profile with manual minirhizotrons deployed within SPRUCE plots. Images were processed using commercial software to quantify the length and diameter of individual roots. Roots were visually assigned to a plant functional type (PFT) of either (ericaceous) shrub, herb (sedges and Maianthemum trifolium), or tree (Larix laricina, Picea mariana) based on expert opinion. The biomass of individual roots was estimated using PFT-specific allometric equations (Iversen et al., 2018). Production per day was estimated as the length of new roots produced between imaging sessions, divided by the number of days between imaging sessions. These values were placed on a m2 aboveground area basis and scaled to a standard depth of 1m (roots are not evenly distributed with depth, do not interpret value as being on a m3 basis). Maximum and average (weighted by production length) depth of each PFT were also estimated within each minirhizotron tube. Annual production was interpolated as the average of four methods to scale these data (see Weber et al, 2026). Standing crop of roots was estimated for each tube as the maximum visible amount (both length and mass) of roots of that PFT for that year. These data expand the ability of researchers to accurately estimate the belowground dynamics of peatland vegetation, as well as the role that fine roots may play in impacting the fluxes of carbon within peatlands. This dataset contains three data files in comma-separate values (*.csv) format. This dataset contains one data file in comma-separate values (.csv) format. Additional metadata are provided: three data dictionaries and a file-level metadata file in comma-separate values (.csv) format and a user guide in PDF (*.pdf) format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Imaging of a fluid injection process using geophysical data — A didactic example

In many subsurface industrial applications, fluids are injected into or withdrawn from a geologic formation. It is of practical interest to quantify precisely where, when, and by how much the injected fluid alters the state of the subsurface. Routine geophysical monitoring of such processes attempts to image the way that geophysical properties, such as seismic velocities or electrical conductivity, change through time and space and to then make qualitative inferences as to where the injected fluid has migrated. The more rigorous formulation of the time-lapse geophysical inverse problem forecasts how the subsurface evolves during the course of a fluid-injection application. Using time-lapse geophysical signals as the data to be matched, the model unknowns to be estimated are the multiphysics forward-modeling parameters controlling the fluid-injection process. Properly reproducing the geophysical signature of the flow process, subsequent simulations can predict the fluid migration and alteration in the subsurface. The dynamic nature of fluid-injection processes renders imaging problems more complex than conventional geophysical imaging for static targets. This work intents to clarify the related hydrogeophysical parameter estimation concepts.

58 GEOSCIENCES↗