SEARCH · Search NASA

Results for “Reservoirs”

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 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricityand river water management for industry and agriculture. Although the hydrological and geophysical impacts ofwater relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack ofin situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial waterstorage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satelliteobservations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The studysite is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia.Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that wasclearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period(from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200mmequivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, aforward model is developed to determine the increased water surface level corresponding to GRACE observations,estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannotcapture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the landsubsidence was calculated using the disk load model constructed based on the GRACE-derived lake level andLandsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the westerncoast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, whichcauses the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation ofthe surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidencefeature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE inconjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smallerspatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-basedalgorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysicalchanges from local to regional scales.

Tangdamrongsub, Natthachet↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS2026

Synchronized Eruptions on Io: Possible Evidence of Interconnected Subsurface Magma Reservoirs

On 27 December 2024, Juno's JIRAM (Jovian InfraRed Auroral Mapper) instrument observed an unprecedented volcanic event in Io's southern hemisphere, covering a vast region of ∼65,000 km2, near 73°S, 140°E. Within the imaged region, only one hot spot was previously known (Pfd454). This feature was earlier estimated to cover an area of 300 km2 with a total power output of 34 GW. JIRAM results show that the region produces a power output of 140–260 TW, over 1,000 times higher than earlier estimates and likely exceeding the brightest eruption ever recorded on Io, that of Surt in 2001 (∼80 TW). Three adjacent hot spots also exhibited dramatic power increases: P139, PV18, and an unnamed feature south of the main one that surged to ∼1 TW, placing all of them among the top 10 most powerful hot spots observed on Io. A temperature analysis of the features supports the simultaneous onset of these brightenings and suggests a single eruptive event propagating beneath the surface across hundreds of kilometers; this is the first time this has been observed on Io. This in turn would imply a connection among the hotspots' magma reservoirs, while other nearby hotspots that have been known to be active in the recent past, such as Kurdalagon Patera, appear unaffected. The simultaneity supports models of massive, interconnected magma reservoirs. The topology of these regional magma systems may resemble that of a large-scale sponge, in which the massive reservoirs are the pores, interconnected through a largely solid outer shell.

A Mura↗

DOE OSTI2021

Method for optimization of huff-n-puff gas injections in hydrocarbon reservoirs

Methods for optimization of liquid oil production by huff-n-puff in shale reservoirs to achieve an improved (and optimal) oil recovery factor. The process determines and utilizes the optimum huff and puff times, number of cycles and soaking time under practical operation and reservoir conditions. The huff time in the process is a period so long that the pressure near the wellbore reaches the set maximum injection pressure during the huff period. The puff time in the process is the time required for the pressure near the wellbore to reach the set minimum production pressure during the puff period. Soaking is typically not necessary during the huff-n-puff gas injection in shale oil reservoirs. The number of huff-n-puff cycles is determinable by the time in which the economic rate cut-off is reached.

Sheng, James J.↗

DOE OSTI2023

Method for optimization of huff-n-puff gas injection in hydrocarbon reservoirs

Sheng, James J.↗

DOE OSTI2025

Initial Conditions of Planet Formation: Time Constraints from Small Bodies and the Lifetime of Reservoirs in the Solar Protoplanetary Disk

This review explores the timescales of the initial phase of planet formation, from nebular dust (CAIs and chondrules) to planetesimal accretion and differentiation, using evidence from meteorite research. Aluminium-Mg systematics of CAIs indicate either an extended period of CAI formation (∼0.3 Ma) or an initial 26 Al heterogeneity, with evidence supporting a homogeneous 26 Al abundance in the protoplanetary disk. Thermal and aqueous alteration on the parent body can disturb the U-Pb and Al-Mg chronometers in chondrules. Focusing on relatively robust isochron data from plagioclase of pristine (types ≤3.05) chondrites indicates a shift in chondrule formation locations, moving from the inner to the outer disk over time. Ages of basaltic achondrites show that silicate differentiation on small bodies was well underway within the first few million years (Ma) of our solar system. Their age record, however, reveals inconsistencies between different chronometers, partly caused by (i) secondary disturbances due to thermal metamorphism, aqueous alteration, or impacts, (ii) the presence of xenolithic minerals, and (iii) potentially variable initial 26 Al abundances due to disturbances at the mineral scale. Nucleosynthetic isotope data indicate that parent bodies of iron and stony meteorites formed in two distinct regions within the protoplanetary disk: the inner, non-carbonaceous (NC) and the outer, carbonaceous (CC) region. Based on Hf-W chronometry it has been demonstrated that NC and CC parent bodies of magmatic iron meteorites segregated their cores within ∼1–3 Ma after CAI formation, implying that parent body accretion occurred within <1 Ma in both reservoirs. Combining accretion ages with nucleosynthetic data further reveals that, at first order, NC and CC reservoirs in the solar protoplanetary disk were established within 1 Ma and existed over several Ma with limited exchange between them. In the CR chondrite accretion region of the disk, planetary bodies formed over at least 3 Ma, while in most other regions, formation spanned at least 1 Ma, with minimal changes in nucleosynthetic isotope compositions. Aerodynamical size sorting of dust likely introduced or amplified some of these variations.

26Al distribution↗

DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗

DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗

DOE OSTI2020

CO 2 plume evolution in a depleted natural gas reservoir: Modeling of conformance uncertainty reduction over time

Uncertainty in the long-term fate of CO 2 injected for geologic carbon sequestration (GCS) is a significant barrier to the adoption of GCS as a greenhouse-gas emission-mitigation for industry and regulatory agencies alike. Here, we present a modeling study that demonstrates that the uncertainty in forecasts of GCS site performance decreases over time as monitoring data are used to update operational models. We consider a case study of GCS in a depleted natural gas reservoir, with CO 2 injection occurring over 20 years, with a 50-year post-injection site care period. We constructed a detailed model to generate the actual model output, which is considered synthetic observation data. A series of simpler operational models based on limited data and assumptions about how an operator would model such a site are then run and compared against actual model output at specific monitoring points after one year, two years, etc. The operational model is updated and improved using the synthetic observation data from the actual model at the same time intervals. Model parameter values and model features needed to be updated over time to improve matches to the actual model. These kinds of model adjustments would be a normal part of reservoir engineering and site management at GCS sites. Uncertainty in two key measures related to site performance decreases with time: extent of the CO 2 plume up-dip migration, and radial extent of the pressure pulse. This conclusion should help allay the concerns of industry and regulators about uncertainty in long-term fate of CO 2 at GCS sites.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Gigaton commercial-scale carbon storage and mineralization potential in stacked Columbia River basalt reservoirs

This work presents a detailed supercritical CO 2 storage resource estimation for the stacked basalt reservoirs in the Grande Ronde Basalt of the Columbia River Basalt Group in eastern Washington and Oregon. The assessment aims to derisk the commercialization potential of geologic carbon storage in basalt by leveraging both structural and mineralization trapping of CO 2 in basalt. The structural closures formed by anticlinal ridges and synclinal valleys in Yakima Fold Belt are excellent physical traps to accommodate injected supercritical CO 2 . Rigorous hydraulic testing, well logs and simulation results from the Wallula Basalt Pilot #1 well showed the occurrence of 17 suitable permeable injection zones (up to 2,496 mD) intercalated with dense seals (~2.6E-10 mD) in the Grand Ronde Basalt. In addition, geochemical studies showed fast reactions between supercritical CO 2 and dissolved basalt minerals to form stable carbonates. In conclusion, our calculation indicates up to 40 gigatons (P90) of mineralization storage resources exist in the Grande Ronde Basalt reservoirs.

58 GEOSCIENCES↗

DOE OSTI2020

Controls on reservoir properties in organic-matter-rich shales: Insights from MICP analysis

Mercury injection capillary pressure (MICP) analysis was used to study organic-matter-rich Devonian shales varying in maturity from marginally mature to postmature to discuss controls on drainage and imbibition in shale reservoirs. Six samples come from the New Albany Shale; their total organic carbon (TOC) ranged from 3.3–8.4 % and vitrinite reflectance (R o ) from 0.54–1.42 %. One Marcellus Shale sample with R o 2.50 % was used to extend the maturity range to the advanced dry gas window. Three coal samples of varying rank were also included for comparison. Porosity of the shales ranges from 1.5–5.0 %, and MICP-derived permeability ranges from 0.17–9.9 mD. Injection capillary pressure curves have a range of entry pressures from 0.055 to0.421 MPa, and no relationship of entry pressure with maturity was revealed. Maturity, however, has an influence on pore size distribution, with the higher-maturity samples having significantly lower volumes of pores with throats less than 10 nm in diameter. Withdrawal capillary pressure curves document a wide range of withdrawal efficiencies from 55 % to almost 100 %, and pores having throats between 3–50 nm in diameter play a major role in the ability of mercury to be released from the pore system. The obtained data provide valuable information about shale properties and their behavior not only in the air-mercury system but also under reservoir conditions.

02 PETROLEUM↗

DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗

DOE OSTI2021

Immunological Correlates of the HIV-1 Replication-Competent Reservoir Size

Understanding what shapes the latent human immunodeficiency virus type 1 (HIV-1) reservoir is critical for developing strategies for cure. Here, we measured frequency of persistent HIV-1 infection after 5 years of suppressive antiretroviral therapy initiated during chronic infection. Pretreatment CD8 + T-cell activation, nadir CD4 count, and CD4:CD8 ratio predicted reservoir size.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Utah FORGE: EGS Reservoir Produced Fluids Geochemistry 2022-2024

This dataset contains geochemical analyses of produced fluids from the Utah FORGE site, specifically from wells 16A(78)-32, 16B(78)-32, and 58-32, collected during various stimulation, flowback, and circulation tests conducted between 2022 and 2024. The data contains element concentrations, pH, and dissolved gas compositions. Geochemical analyses for 2022 and 2023 were performed at the Brigham Young University geochemistry laboratory, while the 2024 results were obtained from Thermochem. Dry gas samples were collected using a mini-separator attached to the single-phase production line between the wellhead and separator, where fluid was flashed to atmospheric pressure. Gas concentrations were recalculated to a single phase reservoir liquid based on heat and mass balance expressions. Additional contextual information, including interpretations of geochemical trends and reservoir behavior, is available in the included report from Simmons et al. (2025), which was presented at the Stanford Geothermal Workshop in February 2025.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Development of tools for coupled InSAR and Seismicity monitoring of EGS reservoir development and management (Final Report)

The report details a series of software tools and workflows to analyze surface deformation data obtained from InSAR, deformation evidence by earthquake, geological constraints, and reservoir management to define and monitor the geothermal reservoir.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Assessing the Solar Photovoltaic (PV) Potential in Puerto Rican Brownfields and Reservoirs: Detailed Results and Methodology Annex [Slides]

This is the technical annex for "Assessing the Solar Photovoltaic (PV) Potential in Puerto Rican Brownfields and Reservoirs." This study by the National Renewable Energy Laboratory (NREL) evaluates the potential for solar photovoltaic (PV) development on brownfields and reservoirs in Puerto Rico. Redeveloping these sites offers significant benefits, such as expanding energy infrastructure, supporting environmental remediation, and enhancing community revitalization. The study identifies up to 3.3 gigawatts (GW) of solar PV capacity, including 213 megawatts (MW) on closed landfills, 1-2.5 GW on contaminated sites, 78 MW on decommissioned plants, 21-50 MW on transmission line rights-of-way, and 636 MW on waterbodies for floating PV. This research highlights the strategic importance of investments and regulatory support to achieve Puerto Rico's goal of 100% clean energy by 2050.

14 SOLAR ENERGY↗

DOE OSTI2025

Reservoir Dynamics in Proposed Operational Scenarios Where CO 2 -EOR Fields are Transitioned From CO 2 -Flood Enhanced Oil Recovery to Dedicated Carbon Storage: A Field Case Study

This study models the transition of CO 2 -enhanced oil recovery (CO 2 -EOR) fields to dedicated carbon storage using a generalized reservoir model informed by real field data, SACROC. Through scenario-based simulations, the work evaluates how reservoir depletion levels, boundary conditions, fluid properties, and domain size affect pressure buildup, CO 2 plume migration, and long-term containment performance. The results inform practices for repurposing oil fields into reliable CO 2 storage sites.

02 PETROLEUM↗