Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inversion”

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

Supporting Data: Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

Supporting Computational Data from the manuscript titled "ITI and Uranium-Arene s-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial" by Benedikt Kestel, Daniel Pividori, Sabyasachi Roy Chowdhury, Frank W. Heinemann, Bess Vlaisavljevich, and Karsten Meyer published in J. Am. Chem. Soc. (DOI: 10.1021/jacs.5c05029).

Chowdhury, Sabyasachi Roy [Department of Chemistry↗
DOE OSTI2024

ML Classifier Fusion for Three Data Streams with Quality Inversely Proportional to Time Resolution

We consider a monitoring scenario of phenomenon using three different streams of measurements whose quality is proportional to their constant inter-arrival times. Each measurement of a stream needs to be binary-classified to reflect the state of interest of the phenomenon. A set of classifiers is separately trained and fused for each stream at its time resolution using measurements collected under known states. We present a machine learning method to fuse the outputs of these fusers to provide a final classification at the finest time resolution. We show that this fused-fusers method provides decisions with likely superior classification probability compared to the best individual classifiers and fused-classifiers. We derive generalization equations that guarantee a superior classification probability of fused-fusers with a confidence probability specified by the classifiers’ generalization equations. We apply these results to study a practical problem of classifying Pu/Np target dissolution events at a radiochemical processing facility using gamma spectral measurements of effluent flows.

Rao, Nageswara↗
DOE OSTI2024

Phase-inversion polymer composite material, fabricating methods and applications of same

A composite film usable as a separator of an electrochemical device includes hBN nanosheets and at least one polymer. The hBN nanosheets are uniformly dispersed within a matrix of said least one polymer to achieve a highly porous microstructure. Said at least one polymer comprises one or more electrically insulating and electrochemically inert polymers.

Hersam, Mark C.↗
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↗
DOE OSTI2025

Permafrost, Peatland, and Cropland Regions Are Key to Reconciling North American Carbon Sink Estimates

Persistent discrepancies between bottom-up, terrestrial biosphere models (TBMs), and top-down, atmospheric inversions, have made it difficult to quantify the magnitude of the North American terrestrial carbon sink. Previous studies have compared aggregated continent-scale estimates of carbon fluxes from TBMs and inversions for all of North America, but this provides limited insights into finer-scale mismatches that contribute to the overall discrepancies. Here we evaluate agreement between TBM and inversion carbon flux estimates at 1° × 1° resolution to provide more direct insights into where models disagree and what underlying factors drive discrepancies. We find that the additional carbon uptake estimated by inversions, in just 16% of the area of North America, is large enough to account for the discrepancy between TBMs and inversions across the whole continent. The majority of these differences occur in permafrost, peatland, and cropland regions. In these regions, we find a higher likelihood of potential biases in the weaker sink estimates from TBMs, suggesting that the stronger sink implied by inversions is more likely to be realistic. However, the current observational coverage is insufficient for fully assessing the causes of discrepancies or the magnitude of biases in either approach. Encouragingly, improved representation of agricultural processes in a TBM led to better agreement with inversions in croplands. Efforts to accurately model cropland dynamics will help improve agreement between TBMs and inversions. Overall, this work presents a clear path for reconciling the discrepancies between inversion and TBM estimates of the North American carbon sink that have persisted for two decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Global Methane Budget 2000–2020

Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. CH4 is the second most important human-influenced greenhouse gas in terms of climate forcing after carbon dioxide (CO2), and both emissions and atmospheric concentrations of CH4 have continued to increase since 2007 after a temporary pause. The relative importance of CH4 emissions compared to those of CO2 for temperature change is related to its shorter atmospheric lifetime, stronger radiative effect, and acceleration in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in quantifying the factors responsible for the observed atmospheric growth rate arise from diverse, geographically overlapping CH4 sources and from the uncertain magnitude and temporal change in the destruction of CH4 by short-lived and highly variable hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to improve, synthesise, and update the global CH4 budget regularly and to stimulate new research on the methane cycle. Following Saunois et al. (2016, 2020), we present here the third version of the living review paper dedicated to the decadal CH4 budget, integrating results of top-down CH4 emission estimates (based on in situ and Greenhouse Gases Observing SATellite (GOSAT) atmospheric observations and an ensemble of atmospheric inverse-model results) and bottom-up estimates (based on process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). We present a budget for the most recent 2010–2019 calendar decade (the latest period for which full data sets are available), for the previous decade of 2000–2009 and for the year 2020. The revision of the bottom-up budget in this 2025 edition benefits from important progress in estimating inland freshwater emissions, with better counting of emissions from lakes and ponds, reservoirs, and streams and rivers. This budget also reduces double counting across freshwater and wetland emissions and, for the first time, includes an estimate of the potential double counting that may exist (average of 23 Tg CH4 yr−1). Bottom-up approaches show that the combined wetland and inland freshwater emissions average 248 [159–369] Tg CH4 yr−1 for the 2010–2019 decade. Natural fluxes are perturbed by human activities through climate, eutrophication, and land use. In this budget, we also estimate, for the first time, this anthropogenic component contributing to wetland and inland freshwater emissions. Newly available gridded products also allowed us to derive an almost complete latitudinal and regional budget based on bottom-up approaches. For the 2010–2019 decade, global CH4 emissions are estimated by atmospheric inversions (top-down) to be 575 Tg CH4 yr−1 (range 553–586, corresponding to the minimum and maximum estimates of the model ensemble). Of this amount, 369 Tg CH4 yr−1 or ∼ 65 % is attributed to direct anthropogenic sources in the fossil, agriculture, and waste and anthropogenic biomass burning (range 350–391 Tg CH4 yr−1 or 63 %–68 %). For the 2000–2009 period, the atmospheric inversions give a slightly lower total emission than for 2010–2019, by 32 Tg CH4 yr−1 (range 9–40). The 2020 emission rate is the highest of the period and reaches 608 Tg CH4 yr−1 (range 581–627), which is 12 % higher than the average emissions in the 2000s. Since 2012, global direct anthropogenic CH4 emission trends have been tracking scenarios that assume no or minimal climate mitigation policies proposed by the Intergovernmental Panel on Climate Change (shared socio-economic pathways SSP5 and SSP3). Bottom-up methods suggest 16 % (94 Tg CH4 yr−1) larger global emissions (669 Tg CH4 yr−1, range 512–849) than top-down inversion methods for the 2010–2019 period. The discrepancy between the bottom-up and the top-down budgets has been greatly reduced compared to the previous differences (167 and 156 Tg CH4 yr−1 in Saunois et al. (2016, 2020) respectively), and for the first time uncertainties in bottom-up and top-down budgets overlap. Although differences have been reduced between inversions and bottom-up, the most important source of uncertainty in the global CH4 budget is still attributable to natural emissions, especially those from wetlands and inland freshwaters. The tropospheric loss of methane, as the main contributor to methane lifetime, has been estimated at 563 [510–663] Tg CH4 yr−1 based on chemistry–climate models. These values are slightly larger than for 2000–2009 due to the impact of the rise in atmospheric methane and remaining large uncertainty (∼ 25 %). The total sink of CH4 is estimated at 633 [507–796] Tg CH4 yr−1 by the bottom-up approaches and at 554 [550–567] Tg CH4 yr−1 by top-down approaches. However, most of the top-down models use the same OH distribution, which introduces less uncertainty to the global budget than is likely justified. For 2010–2019, agriculture and waste contributed an estimated 228 [213–242] Tg CH4 yr−1 in the top-down budget and 211 [195–231] Tg CH4 yr−1 in the bottom-up budget. Fossil fuel emissions contributed 115 [100–124] Tg CH4 yr−1 in the top-down budget and 120 [117–125] Tg CH4 yr−1 in the bottom-up budget. Biomass and biofuel burning contributed 27 [26–27] Tg CH4 yr−1 in the top-down budget and 28 [21–39] Tg CH4 yr−1 in the bottom-up budget. We identify five major priorities for improving the CH4 budget: (i) producing a global, high-resolution map of water-saturated soils and inundated areas emitting CH4 based on a robust classification of different types of emitting ecosystems; (ii) further development of process-based models for inland-water emissions; (iii) intensification of CH4 observations at local (e.g. FLUXNET-CH4 measurements, urban-scale monitoring, satellite imagery with pointing capabilities) to regional scales (surface networks and global remote sensing measurements from satellites) to constrain both bottom-up models and atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) integration of 3D variational inversion systems using isotopic and/or co-emitted species such as ethane as well as information in the bottom-up inventories on anthropogenic super-emitters detected by remote sensing (mainly oil and gas sector but also coal, agriculture, and landfills) to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GKQ9-2RHT (Martinez et al., 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Information theory optimization of signals from small-angle scattering measurements

Small-angle X-ray scattering (SAXS) of particles in solution informs on the conformational states and assemblies of biological macromolecules (bioSAXS) outside of cryo- and solid-state conditions. In bioSAXS, the SAXS measurement under dilute conditions is resolution limited, and through an inverse Fourier transform, the measured SAXS intensities directly relate to the physical space occupied by the particles via the P (r)-distribution. Yet, this inverse transform of SAXS data has been historically cast as an ill-posed, ill-conditioned problem requiring an indirect approach. Here, we show that through the applications of matrix and information theories, the inverse transform of SAXS intensity data is a well-conditioned problem. The so-called ill-conditioning of the inverse problem is directly related to the Shannon number. By exploiting the oversampling enabled by modern detectors, a direct inverse Fourier transform of the SAXS data is possible, provided the recovered information does not exceed the Shannon number. The Shannon limit corresponds to the maximum number of significant singular values that can be recovered in a SAXS experiment, suggesting this relationship is a fundamental property of band-limited inverse integral transform problems. This correspondence reduces the complexity of the inverse problem to the Shannon limit and maximum dimension. We propose a hybrid scoring function using an information theory framework that assesses both the quality of the model-data fit as well as the quality of the recovered P (r)-distribution. The hybrid score utilizes the Akaike information criteria and Durbin-Watson statistic that considers parameter-model complexity, i.e., degrees of freedom, and the randomness of the model-data residuals. The described tests and findings extend the boundaries for bioSAXS by completing the information theory formalism initiated by Peter B. Moore to enable a quantitative measure of resolution in SAXS, robustly determine maximum dimension, and more precisely define the best parameter model appropriately representing the observed scattering data.

Rambo, Robert P. [Science and Technology Facilitie↗
DOE OSTI2024

Neural Networks to Find the Optimal Forcing for Offsetting the Anthropogenic Climate Change Effects

Abstract Of great relevance to climate engineering is the systematic relationship between the radiative forcing to the climate system and the response of the system, a relationship often represented by the linear response function (LRF) of the system. However, estimating the LRF often becomes an ill-posed inverse problem due to high-dimensionality and nonunique relationships between the forcing and response. Recent advances in machine learning make it possible to address the ill-posed inverse problem through regularization and sparse system fitting. Here, we develop a convolutional neural network (CNN) for regularized inversion. The CNN is trained using the surface temperature responses from a set of Green’s function perturbation experiments as imagery input data together with data sample densification. The resulting CNN model can infer the forcing pattern responsible for the temperature response from out-of-sample forcing scenarios. This promising proof of concept suggests a possible strategy for estimating the optimal forcing to negate certain undesirable effects of climate change. The limited success of this effort underscores the challenges of solving an inverse problem for a climate system with inherent nonlinearity. Significance Statement Predicting the climate response for a given climate forcing is a direct problem, while inferring the forcing for a given desired climate response is often an inverse, ill-posed, problem, posing a new challenge to the climate community. This study makes the first attempt to infer the radiative forcing for a given target pattern of global surface temperature response using a deep learning approach. The resulting deeply trained convolutional neural network inversion model shows promise in capturing the forcing pattern corresponding to a given surface temperature response, with a significant implication on the design of an optimal solar radiation management strategy for curbing global warming. This study also highlights the technical challenges that future research should prioritize in seeking feasible solutions to the inverse climate problem.

Ren, Huiying↗
DOE OSTI2025

Quantifying uncertainty in regional-scale seismic moment tensors

We examine the ability of three different inversion methods: 1) first motion (FM) inversions, 2) amplitude inversions, and 3) full waveform (FW) time variable moment tensor (TVMT) inversions, to recover an accurate source mechanism for simulated data for an earthquake as well as an explosion. The ability of inversion methods, such as the ones described above, to recover accurate models representative of the data depends on both a priori information as well as the quality of data being inverted. Therefore, we examine the effect that station geometry, geologic model, and noise level has on the inversion and estimated source mechanism. We find that FM data can provide more robust solutions than amplitude data and are not as sensitive to inaccurate earth models, especially in low-noise cases, and overall FW inversions are the most accurate out of the methods examined, but can still be biased by inaccurate earth models.

58 GEOSCIENCES↗
DOE OSTI2025

Machine learning enabled discovery of superhard and ultrahard carbon polymorphs

The demand for multifunctional materials has motivated the move from near-equilibrium materials to metastable i.e. out-of-equilibrium phases that can meet several desired target properties. The search for such metastable phases with exotic properties is non-trivial and often serendipitous. Inverse design approaches based on evolutionary search have been powerful tools, but such traditional searches have focused on identifying primarily stable and metastable materials with the lowest enthalpy. The inverse design of materials, with a focus on a desired property such as, for example, hardness is a challenging task because of the expensive computational cost involved in sampling multiple structures. The recent advances in machine learning have brought new powerful AI techniques to the forefront which can potentially revolutionize the inverse design and discovery of materials, especially metastable phases capable of meeting multifunctionality. Here, in this work, we develop and apply an automated reinforcement learning workflow for inverse design that integrates first principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the superhard and ultrahard metastable phases of Carbon. We demonstrate an automatic machine learning based inverse design workflow to map new undiscovered metastable states ranging from near equilibrium to those far-from-equilibrium that satisfy multiple property objectives, specifically bulk moduli, shear moduli and hardness. We create a comprehensive library of carbon stable and metastable phases with varying hardness and subsequently shortlist 10 top performing candidate carbon structures, including two newly reported phases, based on their hardness and characterize their temperature dependent mechanical properties. A neural network model is built using featurization of allotropes of carbon to predict the quasi-harmonic Gibbs free energies. The Gibbs free energies of the top performing phases are analyzed to get an estimate of the experimental synthesizability of these superhard and ultrahard carbon phases. In general, we show using machine learning based inverse design approaches how hitherto inaccessible metastable states can be identified and potentially synthesized to meet the demand for multifunctional materials.

Balasubramanian, Karthik [Univ. of Illinois, Chica↗
DOE OSTI2026

Dependence of wind-farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration

This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .

17 WIND ENERGY↗
DOE OSTI2024

EPCAPE-PT-LANL Measurements: Photoacoustic Extinctometer

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Photoacoustic Extinctiometer at 870 nm (Droplet Measurement Technology) Header: - AbsorptionNoise_870nm[1/Mm]: Noise level associated with the absorption measurement at 870 nm, expressed in inverse megameters (Mm⁻¹). - Absorption_870nm[1/Mm]: Absorption coefficient of aerosols measured at 870 nm, expressed in inverse megameters (Mm⁻¹). - Extinction_870nm[1/Mm]: Extinction coefficient of aerosols measured at 870 nm, expressing the sum of scattering and absorption by aerosols, in inverse megameters (Mm⁻¹). - Scattering_870nm[1/Mm]: Scattering coefficient of aerosols measured at 870 nm, expressed in inverse megameters (Mm⁻¹). - BackgroundAbsorption_870nm[1/Mm]: Background absorption measurement at 870 nm, used to correct the primary absorption data, expressed in inverse megameters (Mm⁻¹). - BackgroundScattering_870nm[1/Mm]: Background scattering measurement at 870 nm, used to correct the primary scattering data, expressed in inverse megameters (Mm⁻¹). - Temperature[C]: The ambient temperature at the time of the measurement, expressed in degrees Celsius. - RelativeHumidity[%]: The relative humidity at the time of the measurement, expressed as a percentage. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
DOE OSTI2024

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. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) 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 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and 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 closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. 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. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗