SEARCH · Search NASA

Results for “forward rate”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Southern Hemisphere In Situ Observations of OH, HO2, CIO and BrO from the ER-2 Aircraft for the 1994 ASHOE Mission

A summary of the first order scientific conclusions that emerged from the research done under this grant are as follows: (1) For the first time, the concentration of the key hydrogen and halogen radicals OH, H02, ClO and BrO were determined on a global scale extending from the arctic circle to the antarctic circle, over the altitude domain of the ER-2. That domain extends from 15-20 km altitude, covering a critical part of the lower stratosphere; (2) Simultaneous, in situ measurements of the concentrations of OH, H02, ClO, BrO, NO and NO2 demonstrate the predominance of odd-hydrogen and halogen free radical catalysis in determining the rate of removal of ozone in the lower stratosphere over the complete ASHOE mission. This extends to the global scale the "first look" data obtained during the NASA Stratospheric Photochemistry and Dynamics Experiment (SPADE), executed out of Ames Research Center in June 1993. This represents a major rearrangement of our understanding with respect to the hierarchy of dominant catalytic cycles controlling ozone loss in the lower stratosphere. For the past twenty years, it has been assumed that nitrogen radicals dominate the destruction rate of ozone in the lower stratosphere; (3) Throughout the altitude and latitude range covered by ASHOE, it was determined that a single catalytic cycle, HO2 + O3 yields OH + 2O2, accounted for one half of the total O3 removal in this region of the atmosphere. Halogen radical catalytic cycles were found to account for one third of the ozone loss, and nitrogen radicals were found to account for 20% of the loss; (4) Simultaneous observations of the full complement of radicals, tracers, ozone, and water vapor during ASHOE demonstrated quantitatively the coupling that exists between the rate limiting radicals and other reactive species in the photochemical reaction network. Specifically, the concentrations of ClO and HO2 are inversely correlated with the concentration of NOx. This carries the implication that the NOx effluent from the proposed High Speed Civil Transport may be less destructive to stratospheric ozone than had previously been thought. ASHOE brought this conclusion forward for the first time on a global basis; and (5) The density of BrO was measured on a global scale during ASHOE in the lower stratosphere. It was found that bromine is responsible for 55-65% of the local rate of catalytic destruction of ozone by reactions involving bromine and chlorine. Normalizing calculated loss rates to total available inorganic bromine and chlorine explicitly demonstrates that bromine is 60-80 times more efficient than chlorine in removing ozone in the lower stratosphere. An inferred value of total inorganic bromine is in excellent agreement with measurements of their source species, organic bromine compounds in the troposphere.

Anderson, James G.↗

NASA NTRS2010

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗

DOE OSTI2021

Life Cycle Assessment and Techno-Economic Assessment of Lithium Recovery from Geothermal Brine

Lithium-ion batteries (LIB) play an essential role in the electrification of the transportation sector, and battery demand for lithium compounds will see a significant increase in the coming decades. This has raised concerns on the supply of lithium, and as a result, technologies are being developed to process unconventional lithium sources. One promising technology is to extract lithium from geothermal brine using lithium-aluminum-layered double hydroxide chloride (LDH) sorbent and forward osmosis. A combined life cycle assessment (LCA) and techno-economic assessment (TEA) is conducted to evaluate the environmental and economic performance of this technology. It is assumed that the lithium extraction unit is an add-on to a 50 MW geothermal power plant located in California. The analysis is based on lab-scale experimental data and stoichiometry while considering the economy of scale for an industrial system. LCA results suggest that, compared with conventional LiOH and Li 2 CO 3 production pathways, the new technology achieves 1–95% reduction in environmental impacts. Even higher reduction can be achieved for LiOH produced via electrolysis. This add-on unit for lithium extraction could achieve a payback period of less than 1 year and reach net present values of $454M and $315M and internal rates of return of 792 and 1130% for LiOH and Li 2 CO 3 production pathways, respectively. The favorable environmental and economic performance suggests that LDH sorption coupled with forward osmosis has great potential to enable the domestic production of battery lithium compounds and that further development should be carried out.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Commercial, industrial, and institutional discount rate estimation for efficiency standards analysis: Sector-level data 1998–2022

Underlying each of the U.S. Department of Energy’s (DOE’s) federal appliance and equipment energy conservation standards are a set of complex analyses of the projected costs and benefits of regulation. Any new or amended standard must be designed to achieve significant additional energy conservation, provided that it is technologically feasible and economically justified (42 U.S.C. 6295(o)(2)(A)). DOE determines economic justification based on whether the benefits exceed the burdens, considering a variety of factors, including the economic impact of the standard on consumers of the product and the savings in lifetime operating cost compared to any increase in price or maintenance expenses (42 U.S.C. 6295(o)(2)(B)). As part of this determination, DOE conducts a life-cycle cost (LCC) analysis, which models the combined impact of appliance first cost and operating cost changes on a representative commercial building sample to identify the fraction of customers achieving LCC savings or incurring net cost at the considered efficiency levels. Thus, the commercial discount rate value(s) used to calculate the present value of energy cost savings within the LCC model implicitly plays a role in estimating the economic impact of potential standard levels. This report provides an in-depth discussion of the commercial discount rate estimation process. It is an update to previous reports on estimating commercial discount rates from firm-level and sector-level financial data (e.g., Fujita, 2021, 2016). Major topics covered in this report include the following: -Discount rate estimation methods and rationale -Data sources used and data limitations -Discount rate distributions for use in standards analysis -Discount rate estimation methods and distributions specific to the small business subgroup analysis A version of this analysis was most recently released in 2022. Going forward, this report will be updated as data allow and analyses necessitate.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1991

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

NASA NTRS2012

Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow

We present preliminary development of an approach for simulating high Reynolds number steady compressible flow in two space dimensions using a Cartesian cut-cell finite volume method. We consider both laminar and turbulent flow with both low and high cell Reynolds numbers near the wall. The approach solves the full Navier-Stokes equations in all cells, and uses a wall model to address the resolution requirements near boundaries and to mitigate mesh irregularities in cut cells. We present a quadratic wall model for low cell Reynolds numbers. At high cell Reynolds numbers, the quadratic is replaced with a newly developed analytic wall model stemming from solution of a limiting form of the Spalart-Allmaras turbulence model which features a forward evaluation for flow velocity and exactly matches characteristics of the SA turbulence model in the field. We develop multigrid operators which attain convergence rates similar to inviscid multigrid. Investigations focus on preliminary verification and validation of the method. Flows over flat plates and compressible airfoils show good agreement with both theoretical results and experimental data. Mesh convergence studies on sub- and transonic airfoil flows show convergence of surface pressures with wall spacings as large as approx.0.1% chord. With the current analytic wall model, one or two additional refinements near the wall are required to obtain mesh converged values of skin friction.

Berger, Marsha↗

NASA NTRSDate not supplied

Decision Aid for Conjunction Risk Mitigation by Differential Drag

In the previous five years, the rate of conjunctions that the NASA Conjunction Assessment Risk Analysis (CARA) team processed and analyzed has more than tripled. (NASA CARA, 2024) New missions in the early development phases are now required to plan for dealing with conjunctions under the present space environment, and also projecting forward into a future likely with even further increased utilization of the space environment. Some missions are investigating the possibility of using differential drag to remediate conjunctions without expending limited fuel or for missions without propulsive capabilities. The NASA CARA team studied the historical record of conjunctions to evaluate the circumstances under which differential drag may be successfully applied and have developed a series of tables to use as a decision aid for missions considering differential drag. CARA records all conjunctions of their protected payloads, with historical records starting in 2005 (with significant conjunction events starting to occur on or after 2013). From this record, approximately 7,300 had a probability of collision (Pc) greater than 1 in 10,000, the nominal requirement for a risk mitigation maneuver (RMM) to be made per NASA Procedural Requirements (NPR) 8079.1 (NASA, 2023), at 3 days prior to the time of closest approach, the time analyzed for differential drag efficacy. CARA’s maneuver trade space tool (MTS) was used to prop-agate the primary satellite forward from that decision point with varying degrees of increase to ballistic coefficient (BC), and then recalculate the Pc to evaluate whether or not the conjunction was mitigated (Pc < 3E-6). These results were then binned and sorted along several axes, including altitude, amount of delta-BC, and (pre-maneuver) rate of energy dissipation, to identify underlying patterns. The altitude plot is shown in Figure 1. We found that differential drag was most successful for satellites with perigees below 560 km, and which could adopt an average delta-BC of 2 or greater (that is, increasing their ballistic coefficient by a factor of 3). However, this is a difficult threshold for a mission to clear; very few spacecraft are capable of adopting a high-drag configuration for 72 hours continuously. Planet’s Dove spacecraft use differential drag to remediate conjunctions (Griffith, et al., 2021), and they have a maximum delta-BC factor of 9, but in practice (with mission and charging constraints) they achieve a time-averaged delta-BC that is closer to 2 (Foster, et al., 2017). In contrast, a NASA mission in development reached out to CARA to evaluate the utility of differential drag to remediate potential conjunctions for a spacecraft which is capable of a similar high delta-BC, but with operational constraints that limit their time-averaged delta-BC to less than 1. CARA has developed tables that can be used as decision aids to advise missions-in-development about the best way to utilize their differential drag capabilities. For missions below 560 km with the operational flexibility to devote multiple days to holding a high-drag configuration (or a sufficiently high drag ratio to compensate for limitations on that time), they are able to successfully remediate high-risk conjunctions. Conversely, missions that do not meet these ex-acting criteria – most missions – can instead be advised to use on-board propulsion systems to perform RMMs or to turn their minimum-area face towards the approach vector, thereby reducing Pc at the moment of conjunction. (NASA, 2023)

risk mitigation↗

DOE OSTI2023

High event rate analysis technique for the dual-axis duo-lateral position-sensitive silicon detectors of FAUST

The dual-axis duo-lateral (DADL) position-sensitive silicon detector was developed to obtain precise position and energy information for detected charged particles. The Forward Array Using Silicon Technology (FAUST) is currently equipped with 68 DADL detectors backed by CsI(Tl) scintillators for the study of charged particle correlations in heavy-ion collisions where precise position and energy information is essential. When conventional signal processing electronics were used for the DADL detectors, a position dependence of the measured energy as well as distortions in the calculated particle positions were observed. In previous work, waveforms from the detector after preamplification were studied to better understand the features that give rise to these distortions; therein, a waveform analysis technique was developed to improve the energy resolution and linearity in position reconstruction. However, the reading and writing of waveforms for an entire detector array limits data collection rates and adds significant burden in data storage and analysis speed. In this work, the integrators of a Struck SIS3316 ADC were utilized to process 228 Th source data to develop and optimize a new analysis method that captures the benefits of the waveform analysis technique while circumventing the waveform writing requirement. This integrator method – capable of 59 keV (FWHM) energy resolution – was used in the collection of 35 MeV/nucleon 28 Si + 12 C collision data using FAUST to investigate exotic decays of highly excited highly deformed nuclei. In this data, a position resolution of 0.4 mm (FWHM) was obtained for 25 MeV α-particles; for α-particles near this energy that originate from 8 Be ground state decays, a 8 Be ground state width of 30 keV (FWHM) was obtained. The impact of the energy-dependent DADL position resolution emergent from electronic noise on the quality of excited state measurement was modeled and compared to the experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1972

A high speed sequential decoder

The performance and theory of operation for the High Speed Hard Decision Sequential Decoder are delineated. The decoder is a forward error correction system which is capable of accepting data from binary-phase-shift-keyed and quadriphase-shift-keyed modems at input data rates up to 30 megabits per second. Test results show that the decoder is capable of maintaining a composite error rate of 0.00001 at an input E sub b/N sub o of 5.6 db. This performance has been obtained with minimum circuit complexity.

Lum, H., Jr.↗

NASA NTRS1991

Sinc-Galerkin estimation of diffusivity in parabolic problems

A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion coefficients in linear partial differential equations is presented. Because the parameter recovery problems are inherently ill-posed, an output error criterion in conjunction with Tikhonov regularization is used to formulate them as infinite-dimensional minimization problems. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which displays an exponential convergence rate and is valid on the infinite time interval. The minimization problems are then solved via a quasi-Newton/trust region algorithm. The L-curve technique for determining an approximate value of the regularization parameter is briefly discussed, and numerical examples are given which show the applicability of the method both for problems with noise-free data as well as for those whose data contains white noise.

Smith, Ralph C.↗

NASA NTRS1990

Challenges in the 1990's for astronaut training simulators

New challenges for the simulation community at the Johnson Space Center both in near and long terms are considered. In the near term, the challenges of supporting an increasing flight rate, maintaining operations while replacing obsolete subsystems, and incorporating forthcoming changes to the Space Shuttle are discussed, and focus is placed on a change of forward flight-deck instruments from electro-mechanical devices to electronic displays. Training astronauts for complex concurrent missions involving multiple spacecraft and geographically dispersed ground facilities is considered to be foremost of the long-term challenges, in addition to the tasks of improving the simulator reliability and the operational efficiency of the facilities.

Brown, Patrick M.↗

NASA NTRS2022

A Generalized Dual-Frequency Ratio (DFR) Approach for Rain Retrievals

The dual-frequency ratio of radar reflectivity factors (DFR) has been shown to be a useful quantity as it is independent of the number concentration of the particle size distribution and primarily a function of the mass-weighted particle diameter, Dm. A drawback of DFR-related methods for rain estimation, however, is the non-unique relationship between Dm and DFR. At Ku- and Ka-band frequencies, two solutions for Dm exist when DFR is less than zero. This ambiguity generates multiple solutions for the range profiles of the particle size parameters. We investigate characteristics of these solutions for both the initial-value (forward) and final-value (backward) forms of the equations. To choose one of among many possible range profiles of Dm, number concentration and rain rate, R, independently measured path attenuations are used. For the backward approach the possibility exists of dispensing with externally measured path attenuations by achieving consistency between the input and output path attenuations. The methods are tested by means of a simulation based on disdrometer-measured raindrop size distributions and the results are compared with a simplified version of the operational R-Dm method.

Robert Meneghini↗

NASA NTRS1982

Development of an LSI maximum-likelihood convolutional decoder for advanced forward error correction capability on the NASA 30/20 GHz program

The particular coding option identified as providing the best level of coding gain performance in an LSI-efficient implementation was the optimal constraint length five, rate one-half convolutional code. To determine the specific set of design parameters which optimally matches this decoder to the LSI constraints, a breadboard MCD (maximum-likelihood convolutional decoder) was fabricated and used to generate detailed performance trade-off data. The extensive performance testing data gathered during this design tradeoff study are summarized, and the functional and physical MCD chip characteristics are presented.

Clark, R. T.↗

NASA NTRS1999

The Yearly Variation in Fall-Winter Arctic Winter Vortex Descent

Using the change in HALOE methane profiles from early September to late March, we have estimated the minimum amount of diabatic descent within the polar which takes place during Arctic winter. The year to year variations are a result in the year to year variations in stratospheric wave activity which (1) modify the temperature of the vortex and thus the cooling rate; (2) reduce the apparent descent by mixing high amounts of methane into the vortex. The peak descent amounts from HALOE methane vary from l0km -14km near the arrival altitude of 25 km. Using a diabatic trajectory calculation, we compare forward and backward trajectories over the course of the winter using UKMO assimilated stratospheric data. The forward calculation agrees fairly well with the observed descent. The backward calculation appears to be unable to produce the observed amount of descent, but this is only an apparent effect due to the density decrease in parcels with altitude. Finally we show the results for unmixed descent experiments - where the parcels are fixed in latitude and longitude and allowed to descend based on the local cooling rate. Unmixed descent is found to always exceed mixed descent, because when normal parcel motion is included, the path average cooling is always less than the cooling at a fixed polar point.

Schoeberl, Mark R.↗

NASA NTRS1978

Effect of tropospheric aerosols upon atmospheric infrared cooling rates

An investigation has been made of the impact of wind-blown dust particles upon local climate of arid regions. The case of Northwest India is specifically considered, where a dense layer of dust persists for several months during the summer. In order to examine the effect of this dust layer on the infrared radiative flux and cooling rates, a method is presented for calculating the IR flux within a dusty atmosphere which allows the use of gaseous band models and is applicable in the limit of small single scattering albedo and pronounced forward scattering. The participating components of the atmosphere are assumed to be water vapor and spherical quartz particles only. The atmospheric window is partially filled by including the water vapor continuum bands for which empirically obtained transmission functions have been used. It is shown that radically different conclusions may be drawn on dust effects if the continuum absorption is not considered. The radiative transfer model, when applied to a dusty atmosphere, indicates that there is a moderate enhancement in the atmospheric greenhouse and a 10% increase in the mean IR radiative cooling rate, relative to the dust free case, within the lower troposphere. These results have been compared with previous work by other authors in the context of the possibility of dust layers inhibiting local precipitation.

Harshvardhan, MR.↗

NASA NTRS2012

Motion-Base Simulator Evaluation of an Aircraft Using an External Vision System

Twelve air transport-rated pilots participated as subjects in a motion-base simulation experiment to evaluate the use of eXternal Vision Systems (XVS) as enabling technologies for future supersonic aircraft without forward facing windows. Three head-up flight display concepts were evaluated -a monochromatic, collimated Head-up Display (HUD) and a color, non-collimated XVS display with a field-of-view (FOV) equal to and also, one significantly larger than the collimated HUD. Approach, landing, departure, and surface operations were conducted. Additionally, the apparent angle-of-attack (AOA) was varied (high/low) to investigate the vertical field-of-view display requirements and peripheral, side window visibility was experimentally varied. The data showed that lateral approach tracking performance and lateral landing position were excellent regardless of AOA, display FOV, display collimation or whether peripheral cues were present. However, the data showed glide slope approach tracking appears to be affected by display size (i.e., FOV) and collimation. The monochrome, collimated HUD and color, uncollimated XVS with Full FOV display had (statistically equivalent) glide path performance improvements over the XVS with HUD FOV display. Approach path performance results indicated that collimation may not be a requirement for an XVS display if the XVS display is large enough and employs color. Subjective assessments of mental workload and situation awareness also indicated that an uncollimated XVS display may be feasible. Motion cueing appears to have improved localizer tracking and touchdown sink rate across all displays.

Kramer, Lynda J.↗

DOE OSTI2026

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗