Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir computing”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Novel Proppant Logging Technique for Infill Drilling of Unconventional Shale Wells

Summary During the development of an unconventional play, wells are drilled and completed in batches, and depending on the development plans, current and expected energy market trends, as well as other developmental considerations, new wells are drilled and hydraulically fractured later near existing producing laterals. This creates challenges in terms of optimizing resource recovery and reducing interwell communication. A novel approach is proposed that utilizes systematic composite sampling and analysis of drilling mud returns to look for and quantitatively identify sand particles. The workflow involves cleaning, drying, and segregation of samples into sizes of interest to us (size distribution of pumped proppant in offset parent wells). These samples are imaged at a very high resolution and analyzed for grains using characteristic optical imaging properties to classify proppant sand particles using computer vision algorithms. Further analysis, such as elemental compositional analysis, is used to validate the results from the imaging workflow. We present a case study from the Permian Basin, where a new child well was used as a test case to prove this technology at the Hydraulic Fracturing Test Site (HFTS-2) in Delaware Basin. We introduce new proppant parameters that help identify sustained proppant zones vs. localized propped fractures. We have used additional diagnostics and data collected at the test site to validate observations from the proppant log and have successfully interpreted significantly propped vs. unpropped zones. A key finding from this test has been the significant proppant transport distances observed away from parent wells. Observable proppant was found at a lateral distance of approximately 425 m for one set of parent wells and more than 915 m for another set of parent wells. While a major limitation of this technique is the sampling rate, given adequate sampling, the proposed technology represents a systematic and one-of-a-kind interpretation of spatial proppant distribution while drilling infill wells. It provides us with unique opportunities to better understand the current state of the reservoir being targeted, including zones that are likely highly drained relative to others, and how the planned hydraulic fracturing of child wells can be improved.

Energy & Fuels↗
DOE OSTI2022

Simulink Modeling and Dynamic Study of Fixed-Speed, Variable-Speed, and Ternary Pumped Storage Hydropower

Pumped Storage Hydropower (PSH) is one of the most popular energy storage technologies in the world. It uses an upper reservoir to store water which can be later used during high-demand. In the United States, most of the energy storage capability actually corresponds to PSH. Moreover, PSH also brings multiple benefits to grid operation. This report presents the Simulink models of three common PSH technologies: Fixed-Speed (FS), Variable-Speed (VS), and Ternary (T)-PSH. These models are available to the general public on this GitHub repository, which contains the MATLAB model initialization files, the Simulink model files, and supplementary MATLAB code used to obtain the figures in this work. For each PSH model, an introductory description of the model components and other relevant functionalities are provided. For further information regarding the models and the initialization parameters, the reader is referred to the shared files in the repository. This report also presents the dynamic behavior of each model. The response of such models to a load event is analyzed and matched with each model's features. A custom IEEE 39 bus case is employed for the FS and T-PSH simulations, while the VS-PSH is simulated on a simplified three-bus test system due to the computational complexity of the model. For the T-PSH, the steady-state and the switching between several operating modes are also studied in this work.

13 HYDRO ENERGY↗
DOE OSTI2026

Illinois Storage Corridor, CarbonSAFE Phase III: UIC Class VI Permitting Plan

The Illinois Storage Corridor (ISC) project evaluated two distinct sites to determine the feasibility of commercial-scale CO₂ storage at each. The project leveraged the region's exceptional geological characteristics, particularly the well-characterized Cambro-Ordovician Storage Complex, to enable permanent geological storage of more than 50 million tonnes of CO₂ over 30 years. The two storage sites are located at One Earth Energy (OEE) facility in northcentral Illinois and Prairie State Generating Company (PSGC) in southcentral Illinois. Once operational, these facilities will combine to capture and store more than 6.5 million tonnes of CO₂ per year, positioning the ISC among the largest carbon storage regions globally. Preliminary homogeneous dynamic modeling based on regional and site-specific reservoir characteristics indicates promising injection capabilities at both locations. For the OEE site, modeling predicts a maximum allowable injection rate of 3.7 MTPA, with a baseline scenario of 1.7 MTPA over 30 years producing a CO₂ plume radius of 1.4 miles at end of injection. For the PSGC site, incorporating recent well data, modeling indicates a single-well maximum injection rate of 2.1 MTPA, with a plume radius of 4.2 miles at end of injection for the 60 MT over 30 years scenario. The primary objective of this CarbonSAFE Phase III project is to develop and submit Class VI Underground Injection Control (UIC) permit applications to the U.S. Environmental Protection Agency Region 5. The permitting plan outlines comprehensive site characterization, Area of Review delineation, monitoring programs, well construction designs, financial responsibility provisions, and post-injection site care procedures necessary to demonstrate safe, permanent CO₂ storage protective of underground sources of drinking water. Three UIC Class VI permit applications for the OEE site were submitted to EPA in October 2022 and are progressing through technical review, with final permit decision projected by June 2026. Through five rounds of Requests for Additional Information and responses, the applications have been refined to address computational modeling, area of review delineation, well integrity, monitoring protocols, and financial assurance requirements. For the PSGC site, finalized characterization and permitting documentation was delivered directly to the facility in July 2023 due to business constraints precluding formal federal regulatory submission. This comprehensive permitting effort builds upon extensive prior subsurface evaluations and demonstration projects that have confirmed the feasibility of widespread commercial-scale carbon storage in the region.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Using Knowledge-Guided Machine Learning To Assess Patterns of Areal Change in Waterbodies across the Contiguous United States

Lake and reservoir surface areas are an important proxy for freshwater availability. Advancements in machine learning (ML) techniques and increased accessibility of remote sensing data products have enabled the analysis of waterbody surface area dynamics on broad spatial scales. However, interpreting the ML results remains a challenge. While ML provides important tools for identifying patterns, the resultant models do not include mechanisms. Thus, the “black-box” nature of ML techniques often lacks ecological meaning. Using ML, we characterized temporal patterns in lake and reservoir surface area change from 1984 to 2016 for 103,930 waterbodies in the contiguous United States. We then employed knowledge-guided machine learning (KGML) to classify all waterbodies into seven ecologically interpretable groups representing distinct patterns of surface area change over time. Many waterbodies were classified as having “no change” (43%), whereas the remaining 57% of waterbodies fell into other groups representing both linear and nonlinear patterns. This analysis demonstrates the potential of KGML not only for identifying ecologically relevant patterns of change across time but also for unraveling complex processes that underpin those changes.

58 GEOSCIENCES↗
DOE OSTI2021

Pore characterization of the Marcellus Shale by nitrogen adsorption and prediction of its gas storage capacity

In a shale gas reservoir, pore characterization is an important factor used to determine gas storage capacity. However, the nanometer (nm)-scale pore system in shale is difficult to explore by traditional optical, scanning electron microscopy, or even nuclear magnetic resonance well logging. We have investigated the pore structure and storage capacity of the Marcellus Shale through integration of petrophysical analysis from laboratory and well-logging data and nitrogen adsorption. The isotherm of Marcellus Shale is a composite isotherm, which has features of type I, type II, and type IV isotherms with type H4 of the hysteresis loop, suggesting slit-like pores developed in the Marcellus Shale. Quantitative analysis of pore volumes from the nitrogen adsorption indicates that density porosity may be more properly used to approximate the shale porosity and estimate the shale gas volume. In addition, the specific surface area, micropore, and mesopore volumes have a positive relationship with the kerogen volume and total organic content (TOC). By using the Langmuir and Brunauer-Emmet-Teller models, the simulated result indicates that the higher adsorbed quantity of the Marcellus Shale could be the result of the increase of micropore volume contributed, by the increase of kerogen or TOC content. The proposed equations rapidly compute TOC, a key parameter to predict gas storage capacity in overmature shale such as the Marcellus Shale.

Geochemistry & Geophysics↗
DOE OSTI2021

A Systematic Interpretation of Subsurface Proppant Concentration from Drilling Mud Returns: Case Study from Hydraulic Fracturing Test Site (HFTS-2) in Delaware Basin

The aim of this study is generation and validation of a proppant log using analysis of drilling mud returns for child wells. Proppant log provides qualitative as well as quantitative insights into spatial distribution of proppant sand particles from prior stimulation of parent wells. While the basic methodology was developed and formalized during analysis of material collected from through fracture cores at Hydraulic Fracturing Test Site in Midland Basin (HFTS – 1), the test wells at HFTS – 2 in the neighboring Delaware Basin allowed the opportunity to validate the workflow on actual mud return samples from subsurface. As a child well is being drilled, periodic mud return samples are collected at the rig site and preserved for analysis. The workflow involves systematic cleaning of the samples including various steps such as washing, drying and segregation of samples into relevant size fractions of interest (< Mesh 20) based on specifications of pumped sand during stimulation of the parent well. Clean samples are imaged using high resolution transparency scanning. Scan images are then systematically analyzed for particles of interest using computer vision techniques. Sample counts are further validated using elemental analysis of smaller sub-samples at various depths of interest. This step is necessary to isolate proppant versus other naturally occurring minerals such as sulphates and carbonates which show similar optical properties. We successfully correlated proppant distribution against the existing parent well and validated propped versus relatively un-propped zones for a child well at the test site. The advantage of testing the proppant log concept at the HFTS – 2 site is the plethora of additional diagnostic data that is available to validate our primary observations. We can correlate spatial proppant distribution against variability in stimulation response based on independent observations such as image logs, microseismic attributes as well as DAS response, all of which tend to corroborate one another. One of our significant successes was being able to describe varying degrees of impact of the parent well along the lateral length of a stimulated child well. Our workflow represents a systematic and one-of-a-kind interpretation of spatial proppant distribution while drilling child wells. This provides unique opportunities to better understand the current state of the Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S031R003/2477415/urtec-2021-5189-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5189 2 reservoir being targeted including zones which are likely more drained relative to others and how the planned completion of the child well can be improved. Lastly, this log can be useful is validating optimal well spacing in relatively new fields under development.

58 GEOSCIENCES↗
DOE OSTI2025

Heterobimetallic multi-site concerted proton electron transfer (MS-CPET) promotes coordination-induced O–H bond weakening

Coordination-induced bond weakening of X–H bonds (X = O, N, C) has been observed in a number of low-valent transition metal compounds. However, the impact of an appended electron reservoir on the bond dissociation free energy of the O–H bond (BDFE O–H ) of a substrate bound to a d 0 metal is poorly understood. To gain insight into the ability of separated deprotonation and oxidation sites to decrease the BDFE O–H during proton-coupled electron transfer (PCET) reactions, a bimetallic system in which the sites of proton and electron loss are two distinct metal sites is described. Herein, the interconversion of tris(phosphinoamide) Zr/Co complexes HO–Zr(MesNP i Pr 2 ) 3 CoCN t Bu and O$≡$Zr(MesNP i Pr 2 ) 3 CoCN t Bu via hydrogen atom addition/abstraction was studied. Since the Zr center remains in the d 0 Zr IV state throughout these transformations, the electron transfer process is mediated by the appended redox-active Co 0/I center. A series of open-circuit potential (OCP) measurements on the HO–Zr(MesNP i Pr 2 ) 3 CoCN t Bu and O$≡$Zr(MesNP i Pr 2 ) 3 CoCN t Bu complexes was performed, from which the BDFE O–H was found to be 64 ± 1 kcal mol −1 . The BDFE O–H value was further verified through a series of stoichiometric H atom transfer reactions, stoichiometric protonation/deprotonation reactions, and computational studies.

Chemistry↗
DOE OSTI2023

Ab initio quantum scattering calculations and a new potential energy surface for the HCl( X 1Σ+)–O2(X3Σg−) system: Collision-induced line shape parameters for O2-perturbed R(0) 0–0 line in H35Cl

The remote sensing of abundance and properties of HCl—the main atmospheric reservoir of Cl atoms that directly participate in ozone depletion—is important for monitoring the partitioning of chlorine between “ozone-depleting” and “reservoir” species. Such remote studies require knowledge of the shapes of molecular resonances of HCl, which are perturbed by collisions with the molecules of the surrounding air. In this work, we report the first fully quantum calculations of collisional perturbations of the shape of a pure rotational line in H35Cl perturbed by an air-relevant molecule [as the first model system we choose the R(0) line in HCl perturbed by O2]. The calculations are performed on our new highly accurate HCl(X1Σ+)–O2(X3Σg−) potential energy surface. In addition to pressure broadening and shift, we also determine their speed dependencies and the complex Dicke parameter. This gives important input to the community discussion on the physical meaning of the complex Dicke parameter and its relevance for atmospheric spectra (previously, the complex Dicke parameter for such systems was mainly determined from phenomenological fits to experimental spectra and the physical meaning of its value in that context is questionable). We also calculate the temperature dependence of the line shape parameters and obtain agreement with the available experimental data. We estimate the total combined uncertainties of our calculations at 2% relative root-mean-square error in the simulated line shape at 296 K. This result constitutes an important step toward computational population of spectroscopic databases with accurate ab initio line shape parameters for molecular systems of terrestrial atmospheric importance.

Chemistry↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
DOE OSTI2025

Laboratory evaluation of cyclic underground hydrogen storage in the Temblor sandstone of the San Joaquin Basin, California

Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.

08 HYDROGEN↗
DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗
DOE OSTI2020

Sodium Ion Expansion Power Block for Distributed CSP

The Sodium Ion Expansion Power Block for Distributed CSP was a three-plus-one-year effort under the Concentrating Solar Power: Advanced Projects Offering Low LCOE Opportunities (CSP: APOLLO) funding program within the U.S. Department of Energy Solar Energy Technologies Office. The primary objective of this project is to develop a dual-stage modular sodium thermal electrochemical converter (Na-TEC) heat engine power block, which can be potentially integrated with either a small-scale dish solar or large-scale heliostats and parabolic trough CSP. Na-TEC is a heat engine that generates electricity through the isothermal expansion of sodium ions. The Na-TEC is a closed system that can theoretically achieve conversion efficiencies above 45% when operating between thermal reservoirs at 1150 K and 550 K. However, thermal designs have confined previous single-stage devices to thermal efficiencies below 20%. To mitigate some of these limitations, we consider dividing the isothermal expansion into two stages; one at the evaporator temperature (1150 K) and another at an intermediate temperature (650 K –1050 K). This dual-stage Na-TEC takes advantage of regeneration and reheating, and could be amenable to better thermal management. In light of this, we first designed and developed a thermo-electrochemical model, and thermodynamically demonstrated how the dual-stage device can improve the efficiency by up to 8% points over the best performing single-stage device. We also established an application regime map for the single- and dual-stage Na-TEC in terms of the power density and the total thermal parasitic loss. Moreover, a thermal design of an axisymmetric dual-stage Na-TEC is developed to guide the scale-up and fabrication of sub-components of prototype module. A reduced-order finite-element model is used in conjunction with a Na-TEC thermodynamic model that was developed to determine the total parasitic heat loss of this dual-stage design. A number of simplifications are applied in the reduced-order model to decrease the computational time while maintaining acceptable accuracy. According to this analysis, a maximum efficiency of 29% and a maximum power output of 125 W can be achieved. Ultimately, we were able to demonstrate thermal efficiency improvements of the Na-TEC heat engine from 19% up to 40.3%, in a dual-stage (non-optimized) prototype module that we designed, fabricated, and tested with high temperature stage at 923 K. Furthermore, a cost-performance analysis for this improved dual-stage design was carried out for distributed-CSP systems. A high-level techno-economic analysis (TEA) explores four scenarios where a Na-TEC is used as the heat engine for a distributed-CSP system. Overnight capital cost and levelized cost of electricity (LCOE) are estimated for a system lifetime of 30 years, revealing that overnight capital costs in a range from $3.57 to $17.71 per We are feasible, which equate to LCOEs from 6.9 to 17.2 cents/kWh e -1 . This analysis makes a significant contribution by concurrently quantifying the efficiency and unit costs for a range of multistage configurations, and demonstrating that a Na-TEC may be a promising alternative to Stirling engines for distributed-CSP systems at residential scale of 1–5 kW e .

14 SOLAR ENERGY↗
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗
DOE OSTI2023

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Water Network Tool for Resilience (WNTR). User Manual, Version 0.2.3

The Water Network Tool for Resilience (WNTR, pronounced winter) is a Python package designed to simulate and analyze resilience of water distribution networks. Here, a network refers to the collection of pipes, pumps, valves, junctions, tanks, and reservoirs that make up a water distribution system. WNTR has an application programming interface (API) that is flexible and allows for changes to the network structure and operations, along with simulation of disruptive incidents and recovery actions. WNTR is based upon EPANET, which is a tool to simulate the movement and fate of drinking water constituents within distribution systems. Users are encouraged to be familiar with the use of EPANET and/or should have background knowledge in hydraulics and pressurized pipe network modeling before using WNTR. EPANET has a graphical user interface that might be a useful tool to facilitate the visualization of the network and the associated analysis results. Information on EPANET can be found at https://www.epa.gov/water-research/epanet. WNTR is compatible with EPANET 2.00.12 [Ross00]. In addition, users should have experience using Python, including the installation of additional Python packages. General information on Python can be found at https://www.python.org/.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hydropower representation in water and energy system models: a review of divergences and call for reconciliation

Abstract Reservoir-based hydropower systems represent key interactions between water and energy systems and are being transformed under policy initiatives driven by increasing water and energy demand, the desire to reduce environmental impacts, and interacting effects of climate change. Such policies are often guided by complex system models, whereby divergence in system representations can potentially translate to incompatible planning outcomes, thereby undermining any planning that may rely on them. We review different approaches and assumptions in hydropower representation in water and energy systems. While the models and issues are relevant globally, the review focuses on applications in California given its extensive development of energy and water models for policy planning, but discusses the extent to which these observations apply to other regions. Structurally, both water-driven and energy-driven management models are similar. However, in energy models, hydropower is often represented as a single-priority output. Water management models typically allocate water for competing priorities, which are generally uninformed by dynamic electricity load demand, and often result in a lower priority for hydropower. In water models, constraints are increasingly resolved for non-energy components (e.g. inflow hydrology and non-energy water demand); few analogues exist for energy models. These limitations may result in inadequate representations of each respective sector, and vastly different planning outcomes for the same facilities between the two different sectors. These divergent modeling approaches manifest themselves in California where poorly reconciled outcomes may affect decisions in hydropower licensing, electricity grid flexibility and decarbonization, and planning for environmental water. Fully integrated water-energy models are computationally intensive and specific to certain regions, but better representation of each domain in respective efforts would help reconcile divergences in planning and management efforts related to hydropower across energy and water systems.

Rheinheimer, David E. (ORCID:0000000315259069)↗
DOE OSTI2024

Advancing spatiotemporal forecasts of CO 2 plume migration using deep learning networks with transfer learning and interpretation analysis

Accurate and timely forecasts of CO 2 plume distribution throughout the injection and post-injection phases are crucial for detecting plume migration, assessing leakage risks, and supporting operational decisions in geologic carbon storage (GCS). Current convolutional neural network-based approaches primarily focus on spatial information and overlook temporal dependencies in plume distributions, thus limiting their ability to capture dynamic movement effects and provide accurate predictions of plume migration. In this work, we propose two deep learning models, Auto-Encoder (AE)-LSTM and Encoder-Decoder (ED)-ConvLSTM, each uniquely designed to capture both spatial and temporal features. We apply the proposed methods to forecast the dynamic distribution of CO 2 plumes based on 108 reservoir simulations over a 30-year injection and a 30-year post-injection period. The results indicate that the ED-ConvLSTM model outperforms the AE-LSTM model in accurately predicting the spatiotemporal dynamics of CO 2 plume migration, achieving R 2 values above 0.99. To provide a deeper understanding of these model predictions, we employ a gradient-based explanation method on the trained models. This approach provides insights into the influence of input variables on plume migration forecasts and uncovers the underlying prediction mechanisms of the proposed models. Furthermore, we introduce a transfer learning technique, enabling fast and accurate plume migration forecasting in the post-injection phase by leveraging the trained model during the injection phase. This reduces the necessity for extensive data collection or re-training. In conclusion, the methods proposed in our work enhances the performance and interpretability of CO 2 plume migration forecasts, thereby facilitating informed decision-making throughout the entire lifecycle of GCS applications.

58 GEOSCIENCES↗