Search NASA⌕ Search

SEARCH · Search NASA

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 343 records · Page 19

Permeability Prediction Using Vision Transformers

Accurate permeability predictions remain pivotal for understanding fluid flow in porous media, influencing crucial operations across petroleum engineering, hydrogeology, and related fields. Traditional approaches, while robust, often grapple with the inherent heterogeneity of reservoir rocks. With the advent of deep learning, convolutional neural networks (CNNs) have emerged as potent tools in image-based permeability estimation, capitalizing on micro-CT scans and digital rock imagery. This paper introduces a novel paradigm, employing vision transformers (ViTs)—a recent advancement in computer vision—for this crucial task. ViTs, which segment images into fixed-sized patches and process them through transformer architectures, present a promising alternative to CNNs. We present a methodology for implementing ViTs for permeability prediction, its results on diverse rock samples, and a comparison against conventional CNNs. The prediction results suggest that, with adequate training data, ViTs can match or surpass the predictive accuracy of CNNs, especially in rocks exhibiting significant heterogeneity. This study underscores the potential of ViTs as an innovative tool in permeability prediction, paving the way for further research and integration into mainstream reservoir characterization workflows.

58 GEOSCIENCES↗

A reactive tracer method for predicting EGS reservoir geometry and thermal lifetime: development and field validation

The project summarized here (DOE Award No. DE-EE0006764) was intended to develop a methodology for predicting advective heat transfer in fracture-dominated crystalline rock. Our goal was to determine if a combination of inert and reactive tracers could adequately constrain the effective heat transfer surface area between an injection-production well pair. Our approach consisted of: 1. developing a novel computational framework; 2. Performing heat and tracer experiments at meso-scale; and 3. Comparing predictions of advective heat transfer to the “true” thermal breakthrough measured at the Altona site. Below is a summary of project activities/findings, a summary of project tasks, and a conclusion

15 GEOTHERMAL ENERGY↗

BiGEST

Natural products have provided a rich reservoir of beneficial compounds in public health including antibiotics, therapeutics, and immunosuppressants. These natural products are synthesized by enzymes encoded by Biosynthetic Gene Clusters (BGCs), clusters of co-localized biosynthetic genes. Computational detection of BGCs has become a crucial step in natural product discovery. While this process has been facilitated in bacterial and fungal organisms thanks to the currently available tools (e.g., antiSMASH), a large spectrum of eukaryotic organisms have been neglected by these existing tools due to the scarcity and incompleteness of genome annotation resources. Here, we introduce Biosynthetic Gene cluster Extensive Search Tool (BiGEST) to provide an extensive annotation-free search for BGCs in diverse eukaryotic organisms. As a result, BiGEST uncovers eukaryotic BGCs that could be undetected by other BGC detection tools.

Adriani, Lisa↗

Legacy Well Leakage Risk Analysis at the Farnsworth Unit Site

This paper summarizes the results of the risk analysis and characterization of the CO 2 and brine leakage potential of Farnsworth Unit (FWU) site wells. The study is part of the U.S. DOE’s National Risk Assessment Partnership (NRAP) program, which aims to quantitatively evaluate long-term environmental risks under conditions of significant geologic uncertainty and variability. To achieve this, NRAP utilizes risk assessment and computational tools specifically designed to quantify uncertainties and calculate the risk associated with geologic carbon dioxide (CO 2 ) sequestration. For this study, we have developed a workflow that utilizes physics-based reservoir simulation results as input to perform leakage calculations using NRAP Tools, specifically NRAP-IAM-CS and RROM-Gen. These tools enable us to conduct leakage risk analysis based on ECLIPSE reservoir simulation results and to characterize wellbore leakage at the Farnsworth Unit Site. We analyze the risk of leakage from both individual wells and the entire field under various wellbore integrity distribution scenarios. The results of the risk analysis for the leakage potential of FWU wells indicate that, when compared to the total amount of CO 2 injected, the highest cemented well integrity distribution scenario (FutureGen high flow rate) exhibits approximately 0.01% cumulative CO 2 leakage for a 25-year CO 2 injection duration at the end of a 50-year post-injection monitoring period. In contrast, the highest possible leakage scenario (open well) shows approximately 0.1% cumulative CO 2 leakage over the same time frame.

54 ENVIRONMENTAL SCIENCES↗

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗

Computer Program For Variable-Conductance Heat Pipes

VCHPDA provides accurate mathematical models of transient as well as steady-state performance of variable-conductance heat pipes over wide range of operating conditions. Applies to heat pipes with either cold, wicked or hot, nonwicked gas reservoirs and uses ideal-gas law and "flat-front" (negligible vapor diffusion) gas theory. Calculates length of gas-blocked region and temperature of vapor in active portion of heat pipe by solving set of nonlinear equations for conservation of energy and mass. Written in FORTRAN 77.

Antoniuk, D.↗

A Model Tree Generator (MTG) Framework for Simulating Hydrologic Systems: Application to Reservoir Routing

Data-driven algorithms have been widely used as effective tools to mimic hydrologic systems. Unlike black-box models, decision tree algorithms offer transparent representations of systems and reveal useful information about the underlying process. A popular class of decision tree models is model tree (MT), which is designed for predicting continuous variables. Most MT algorithms employ an exhaustive search mechanism and a pre-defined splitting criterion to generate a piecewise linear model. However, this approach is computationally intensive, and the selection of the splitting criterion can significantly affect the performance of the generated model. These drawbacks can limit the application of MTs to large datasets. To overcome these shortcomings, a new flexible Model Tree Generator (MTG) framework is introduced here. MTG is equipped with several modules to provide a flexible, efficient, and effective tool for generating MTs. The application of the algorithm is demonstrated through simulation of controlled discharge from several reservoirs across the Contiguous United States (CONUS).

54 ENVIRONMENTAL SCIENCES↗

Tracer-Guided Characterization of Dominant Pore Networks and Implications for Permeability and Wettability in Shale

Pore network characterization is an important aspect in unconventional reservoir evaluation. While application of the technique of scanning electron microscope (SEM) brings substantial advances in pore characterization in shale, understanding the connected pore network that dominates flow in shale samples is limited by using SEM alone because of small fields of view and lack of views of connectivity in 3D. In this research, a technique integrating tracer imbibition, micro-computed tomography (CT) imaging, and SEM imaging was developed to provide a solution for multiscale imaging in shale. Tracer imbibition indicates pore connectivity; micro-CT imaging after tracer imbibition thus provides an overview of the connected pore network at the millimeter scale. With guidance from micro-CT images after tracer imbibition, a more accurate and detailed characterization of pore systems and related mineralogy can be conducted using higher-resolution SEM. The method was applied to five samples from Wolfcamp and Eagle Ford Formations. Results reveal the effectiveness of the integrated method by showing different patterns of distribution of the dominant pore network and different controlling mineralogy. Dominant porosity, estimated from grayscale analyses, displays a good correlation with permeability. This result indicates that dominant porosity is more relevant to permeability than is total porosity. Results from imbibition tests are also compared with that from contact angle measurement, and important implications on wettability can be obtained. The integrated method thus has the capacity to link the dominant pore network and wettability with microscale to submicroscale mineralogy, which can help better understand the pore systems and fluid flow in shale. Plain Language Summary Pore network characterization and wettability are two important aspects in unconventional reservoir evaluation and are crucial to understanding the nature of shale oil or gas storage and production. Here, a new technique integrating tracer imbibition, micro-computed tomography imaging, and scanning electron microscope imaging was developed for a more representative and effective pore characterization in shale. This integrated technique can provide a direct visualization of the dominant pore network in a millimeter-scale sample and, at the same time, can provide a guided pore network characterization under nanoscale resolution. This method thus can generate results that are more representative and more accurate in a more efficient way compared with the traditional scanning electron microscope method. The new technique can also reveal the local wettability in microscale. High-resolution analysis on the tracer imbibition or fluid displacement and mineralogy can shed light on the local wettability and the controlling mineralogy in shale.

58 GEOSCIENCES↗

How does the connectivity of open-framework conglomerates within multi-scale hierarchical fluvial architecture affect oil-sweep efficiency in waterflooding?

Understanding multi-phase fluid flow and transport processes within aquifers, candidate reservoirs for CO 2 sequestration, and petroleum reservoirs requires understanding a diverse set of geologic properties of the aquifer or reservoir, over a wide range of spatial and temporal scales. We focus on multiphase flow dynamics with wetting (e.g., water) and non-wetting (e.g., gas or oil) fluids, with one invading another. This problem is of general interest in a number of fields and is illustrated here by considering the sweep efficiency of oil during a waterflood. Using a relatively fine-resolution grid throughout a relatively large domain in these simulations and probing the results with advanced scientific visualization tools (Reservoir Visualization Analysis [RVA]/ ParaView software) promote a better understanding of how smaller-scale features affect the aggregate behavior at larger scales. We studied the effects on oil-sweep efficiency of the proportion, hierarchical organization, and connectivity of high-permeability open-framework conglomerate (OFC) cross-sets within the multi-scale stratal architecture found in fluvial deposits. We further analyzed oil production rate, water breakthrough time, and spatial and temporal distribution of residual oil saturation. As expected, the effective permeability of the reservoir exhibits large-scale anisotropy created by the organization of OFC cross-sets within unit bars, and the organization of unit bars within compound- bars. As a result, oil-sweep efficiency critically depends on the direction of the pressure gradient. However, contrary to expectations, the total amount of trapped oil due to the effect of capillary trapping does not depend on the magnitude of the pressure gradient within the examined range. Hence the pressure difference between production and injection wells does not affect sweep efficiency; although the spatial distribution of oil remaining in the reservoir depends on this value. Whether or not clusters of connected OFC span the domain affects only the absolute rate of oil production—not sweep efficiency.

58 GEOSCIENCES↗

CO 2 storage in the Paluxy formation at the Kemper County CO 2 storage complex: Pore network properties and simulated reactive permeability evolution

The Paluxy formation is being considered as a prospective CO 2 reservoir at the Kemper County CO 2 Storage Complex. Herein, the pore and pore-throat size distributions and connectivity of the Paluxy formation is evaluated through analysis of 3D X-ray Computed Tomography images. In spite of resolution limitations that constrain the pore-throat sizes detectable by imaging, the permeability contributing pore-throats are successfully characterized through 3D imaging analysis. Image-obtained pore and pore-throat size distributions and pore connectivity are then utilized to construct pore network models and simulate permeability. After CO 2 is injected, it will dissolve into formation brine and create conditions favorable for dissolution of primary minerals and precipitation of secondary minerals. These reactions will alter the porosity and permeability of the system to varying degrees depending on the spatial location of reactions. In this research, the possible porosity-permeability evolution is simulated using pore network models considering mineral reactions occurring uniformly and non-uniformly throughout the network. For a given change in porosity, there is a large range of possible permeability outcomes. Depending on the extent and spatial location of mineral reactions, permeability may decrease by more than one order of magnitude as minerals precipitate. During dissolution, simulated permeability increases as much as 500%.

03 NATURAL GAS↗

Hydropower Scheduling Toolchains: Comparing Experiences in Brazil, Norway, and USA and Implications for Synergistic Research

While hydropower scheduling is a well-defined problem, there are institutional differences that need to be identified to promote constructive and synergistic research. We study how established toolchains of computer models are organized to assist operational hydropower scheduling in Brazil, Norway, and the United States’ Colorado River System (CRS). These three systems have vast hydropower resources, with numerous, geographically widespread, and complex reservoir systems. Although the underlying objective of hydropower scheduling is essentially the same, the systems are operated in different market contexts and with different alternative uses of water, where the stakeholders’ objectives clearly differ. This in turn leads to different approaches when it comes to the scope, organization, and use of models for operational hydropower scheduling and the information flow between the models. Here, we describe these hydropower scheduling toolchains, identify the similarities and differences, and shed light on the original ideas that motivated their creation. We then discuss the need to improve and extend the current toolchains and the opportunities to synergistic research that embrace those contextual differences.

13 HYDRO ENERGY↗

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗

Co-Optimization of CO 2 -EOR and Storage Processes in Mature Oil Reservoirs

This article presents an optimization methodology for CO 2 enhanced oil recovery in partially depleted reservoirs. A field-scale compositional reservoir flow model was developed for assessing the performance history of an active CO 2 flood and for optimizing both oil production and CO 2 storage in the Farnsworth Unit (FWU), Ochiltree County, Texas. A geological framework model constructed from geophysical, geological, and engineering data acquired from the FWU was the basis for all reservoir simulations and the optimization method. An equation of state was calibrated with laboratory fluid analyses and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). Initial history calibrations of primary, secondary and tertiary recovery were conducted as the basis for the study. After a good match was achieved, an optimization approach consisting of a proxy or surrogate model was constructed with a polynomial response surface method (PRSM). The PRSM utilized an objective function that maximized both oil recovery and CO 2 storage. Experimental design was used to link uncertain parameters to the objective function. Control variables considered in this study included: water alternating gas cycle and ratio, production rates and bottom-hole pressure of injectors and producers. Other key parameters considered in the modeling process were CO 2 purchase, gas recycle and addition of infill wells and/or patterns. The PRSM proxy model was ‘trained’ or calibrated with a series of training simulations. This involved an iterative process until the surrogate model reached a specific validation criterion. A sensitivity analysis was first conducted to ascertain which of these control variables to retain in the surrogate model. A genetic algorithm with a mixed-integer capability optimization approach was employed to determine the optimum developmental strategy to maximize both oil recovery and CO 2 storage. The proxy model reduced the computational cost significantly. The validation criteria of the reduced order model ensured accuracy in the dynamic modeling results. The prediction outcome suggested robustness and reliability of the genetic algorithm for optimizing both oil recovery and CO 2 storage. The reservoir modeling approach used in this study illustrates an improved approach to optimizing oil production and CO 2 storage within partially depleted oil reservoirs such as FWU. Lastly, this study may serve as a benchmark for potential CO 2 –EOR projects in the Anadarko basin and/or geologically similar basins throughout the world.

02 PETROLEUM↗

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

Geophysical and geochemical constraints favoring the capture hypothesis

The early fractionation of the Moon as a result of melting the outer layer of the Moon is discussed. This melting can be attributed to the accretional heating of the moon. Computations of the early thermal evolution of the Earth and Moon based on the assumption that these bodies accreted independently were carried out. The results are consistent with many observations. The concept of geochemical reservoirs was also applied to the Moon.

Turcotte, D. L.↗

StrmtbGen Fortran Program, Version 2

This zip file contains files storing Fortran source code, compiled files, an executable file for Windows computers, and example input and output files for StrmtbGen. StrmtbGen is a Fortran program that generates streamlines and stream tubes. StrmtbGen and StrmtbFlow comprise the FE/NETL CO2 Prophet Model, which is a streamline/stream tube reservoir simulator for modeling CO2 EOR. StrmtbGen's user’s manual as well as StrmtbFlow and its user's manual and a report on the mathematical basis for the FE/NETL CO2 Prophet Model are available on NETL's website under the Collection Name: FE/NETL CO2 Prophet Model.

Morgan, David↗

Impacts of Small-Scale Heterogeneity on Large-Scale Subsurface Flow

Discussion of NETL's CT scanning and core characterization capabilities, followed by a rapid review of several recent publications that utilize these capabilities to answer reservoir questions at the micro scale with relevance to field scale operations.

core characterization↗