Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reservoir flow model”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Calibration of Axisymmetric and Quasi-1D Solvers for High Enthalpy Nozzles

The proposed paper will present a numerical investigation of the flow characteristics and boundary layer development in the nozzles of high enthalpy shock tunnel facilities used for hypersonic propulsion testing. The computed flow will be validated against existing experimental data. Pitot pressure data obtained at the entrance of the test cabin will be used to validate the numerical simulations. It is necessary to accurately model the facility nozzles in order to characterize the test article flow conditions. Initially the axisymmetric nozzle flow will be computed using a Navier Stokes solver for a range of reservoir conditions. The calculated solutions will be compared and calibrated against available experimental data from the DLR HEG piston-driven shock tunnel and the 16-inch shock tunnel at NASA Ames Research Center. The Reynolds number is assumed to be high enough at the throat that the boundary layer flow is assumed turbulent at this point downstream. The real gas affects will be examined. In high Mach number facilities the boundary layer is thick. Attempts will be made to correlate the boundary layer displacement thickness. The displacement thickness correlation will be used to calibrate the quasi-1D codes NENZF and LSENS in order to provide fast and efficient tools of characterizing the facility nozzles. The calibrated quasi-1D codes will be implemented to study the effects of chemistry and the flow condition variations at the test section due to small variations in the driver gas conditions.

Papadopoulos, P. E.↗
DOE OSTI2021

Geomechanical and flow implications with continued bioconversion of coal to methane: Experiments and modeling

Microbial conversion of coal to methane is a promising technology during transition of coal-based energy from conventional coal mining to natural gas recovery from coal. Significant research advances have been made towards engineering ideal microbial communities and nutrients for bio-stimulation of coal. However, actual field applications require geomechanical and flow behavior characterization of coal during the bioconversion process as well as gas production over the life of the created “biogenic gas reservoirs”. This work presents the results of an experimental investigation to estimate gas/methane production using bioconversion of coal and analyses of variations in bulk modulus, strain and permeability. Bulk modulus of coal demonstrated time-dependent behavior with continued bioconversion of solid coal. The phenomenon was modeled using mass balancing in a closed environment and the logistic equation based biogenic gas production with time. The modeled gas production and changes in modulus of coal showed excellent agreement with the corresponding experimental results. Next, numerical simulation of biogenic conversion of coal in a constant stress in-situ condition, replicating the Huff ‘n Puff method under field conditions, was carried out. The results showed that bioconversion can lead to decrease in effective stress, increase in permeability and decrease in the modulus of coal with time. Repeated recharge of coal with nutrients and recovering the produced gas, that is, sequential cycles of Huff n’ Puff, showed further increase in permeability and decrease in coal strength, potentially leading to coal failure in-situ, further increasing the permeability, thus enhancing the prospect of field application of the technology.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2000

The Sugar Model: Catalytic Flow Reactor Dynamics of Pyruvaldehyde Synthesis from Triose Catalyzed by Poly-L-Lysine Contained in a Dialyzer

The formation of pyruvaldehyde from triose sugars was catalyzed by poly-L-lysine contained in a small dialyzer (100 MWCO) suspended in a much larger triose substrate reservoir. The polylysine confined in the dialyzer functioned as a catalytic flow reactor that constantly brought in triose from the substrate reservoir by diffusion to offset the drop in triose concentration within the reactor caused by its conversion to pyruvaldehyde. A 400 mM solution of poly-L-lysine contained in a 0.35 ml dialyzer placed in a 120 ml solution of triose substrate (pH 5.5, 40 C) generated pyruvaldehyde 11 -times faster than an a control reaction without the catalytic dialyzer. However, since the catalytic dialyzer's volume was 343-times smaller than the control reaction, the synthetic intensity (rate/volume) of pyruvaldehyde synthesis within the catalytic dialyzer was 3400-times greater than that of the control reaction and substrate solution. A similar result was obtained using a dialyzer with a 500 MWCO value. Acting as a catalytic flow reactor the polylysine catalytic dialyzer synthesized about 3.5 molecules of pyruvaldehyde per lysine residue in 7 days -- an amount of triose equal to twice the weight of the catalyst. At 7 days the catalytic activity of polylysine was 16% of its initial value, a result indicating catalyst-poisoning caused by reaction of pyruvaldehyde with the e-amino groups of polylysine. The dialyzer method of catalyst containment was selected it provides a simple, flexible, and easily manipulated experimental system for studying the dynamics and evolutionary development of confined autocatalytic processes related to the origin of life under anaerobic conditions.

Weber, Arthur L.↗
NASA NTRS2001

The sugar model: catalytic flow reactor dynamics of pyruvaldehyde synthesis from triose catalyzed by poly-l-lysine contained in a dialyzer

The formation of pyruvaldehyde from triose sugars was catalyzed by poly-l-lysine contained in a small dialyzer with a 100 molecular weight cut off (100 MWCO) suspended in a much larger triose substrate reservoir at pH 5.5 and 40 degrees C. The polylysine confined in the dialyzer functioned as a catalytic flow reactor that constantly brought in triose from the substrate reservoir by diffusion to offset the drop in triose concentration within the reactor caused by its conversion to pyruvaldehyde. The catalytic polylysine solution (400 mM, 0.35 mL) within the dialyzer generated pyruvaldehyde with a synthetic intensity (rate/volume) that was 3400 times greater than that of the triose substrate solution (12 mM, 120 mL) outside the dialyzer. Under the given conditions the final yield of pyruvaldehyde was greater than twice the weight of the polylysine catalyst. During the reaction the polylysine catalyst was poisoned presumably by reaction of its amino groups with aldehyde reactants and products. Similar results were obtained using a dialyzer with a 500 MWCO. The dialyzer method of catalyst containment was selected because it provides a simple and easily manipulated experimental system for studying the dynamics and evolutionary development of confined autocatalytic processes related to the origin of life under anaerobic conditions.

Biogenesis↗
DOE OSTI2021

An Integrated Multiscale Modeling Framework for Unconventional Stimulation and Production (Final Report)

The production of oil and gas from unconventional reservoirs largely depends upon two main features operating at different scales: (1) the establishment of a reservoir scale stimulated fracture network that effectively communicates with the rock volume, enhancing permeability and transport to the wellbore and (2) the coupled multi-phase flow, chemical and mechanical processes affecting the migration of hydrocarbons from the low permeability country rock adjacent to the stimulated fracture network. To date, there has been no simulation framework that allows seamless and integrated prediction of these features across spatial scales extending from the pore structure of the reservoir rock to the volume of the reservoir. In addition, there has been a lack of suitable field measurements to test such models, as stimulation and production data are often proprietary and not freely available to national laboratories and academic institutions. New multi-scale simulation capabilities are needed that are validated against suitable field-based research experiments on hydraulic fracturing and shale production.

02 PETROLEUM↗
DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Primary Electric Propulsion Technology Study

An investigation of the 30-cm engineering-model-thruster technology with emphasis placed on the development of models for understanding and predicting the operational characteristics and wear-out mechanisms of the thruster as a function of operating or design parameters is presented. The task studies include: (1) the wear mechanisms and wear rates that determine the useful lifetime of the thruster discharge chamber; (2) cathode lifetime as determined by the depletion of barium from the barium-aluminate-impregnated-porous-tungsten insert that serves as a barium reservoir; (3) accelerator-grid-system technology; (4) a verification of the high-voltage propellant-flow-electrical-isolator design developed under NASA contract NAS3-20395 for operation at 10-kV applied voltage and 10-A equivalent propellant flow with mercury and argon propellants. A model was formulated for predicting performance.

Poeschel, R. L.↗
DOE OSTI2020

Real Time 3D Observations of Portland Cement Carbonation at CO2 Storage Conditions

Depleted oil reservoirs are considered a viable solution to the global challenge of CO 2 storage. A key concern is whether the wells can be suitably sealed with cement to hinder the escape of CO 2 . Under reservoir conditions, CO 2 is in its supercritical state, and the high pressures and temperatures involved make real-time microscopic observations of cement degradation experimentally challenging. Here, we present an in situ 3D dynamic X-ray micro computed tomography (μ-CT) study of well cement carbonation at realistic reservoir stress, pore-pressure, and temperature conditions. The high-resolution time-lapse 3D images allow monitoring the progress of reaction fronts in Portland cement, including density changes, sample deformation, and mineral precipitation and dissolution. By switching between flow and nonflow conditions of CO 2 -saturated water through cement, we were able to delineate regimes dominated by calcium carbonate precipitation and dissolution. For the first time, we demonstrate experimentally the impact of the flow history on CO 2 leakage risk for cement plugging. In-situ μ-CT experiments combined with geochemical modeling provide unique insight into the interactions between CO 2 and cement, potentially helping in assessing the risks of CO 2 storage in geological reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

Thermal extremes in regulated river systems under climate change: an application to the southeastern U.S. rivers

High river temperatures, or “thermal extremes”, can cause fish mortality and thermoelectric powerplant derating. Under climate change, projected higher air temperature and stronger surface energy fluxes will lead to increased water temperatures, exacerbating thermal extremes. However, cold hypolimnetic releases from thermally stratified reservoirs can depress tailwater temperatures and therefore alleviate thermal extremes. Thermal extremes are more harmful when they coincide with low flows, which we refer to as “hydrologic hot-dry events”. To assess multisectoral impacts of climate change over large regions, we evaluate thermal events according to three impact attributes: duration (D), intensity (I), and severity (S). We apply an established model framework to simulate streamflow and stream temperature over the southeastern US regulated river system. We quantify climate change impacts (by the 2080s under RCP8.5) by comparing historical and future periods and quantify regulation impacts by comparing unregulated and regulated scenarios. We find that climate change will exacerbate thermal extremes (all three metrics) in both unregulated and regulated scenarios, albeit less in the regulated scenario. Thermal mitigation from reservoir regulation will be stronger under climate change, decreasing the three metrics compared to the unregulated case. Even so, thermal extremes in the regulated scenario will still be more severe under climate change, and only 12.2%, 19.7%, and 26.0% of D, I, and S can be mitigated by reservoirs. Despite stronger reservoir stratification, the number of regulated river segments that experience simultaneous high temperature and low flow events (hydrologic hot-dry events) will increase by 21.4% by the 2080s under RCP8.5. These events will have a median annual duration of 10.3 days/year, over 10 times the historical value.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Maximizing oil production from water alternating gas (CO 2 ) injection into residual oil zones: The impact of oil saturation and heterogeneity

Residual oil zones (ROZs) are widespread reservoirs, characterized by oil at residual saturation, either underlying oil fields (brownfield) or lateral (greenfield) to such fields. These reservoirs have the potential to produce volumes of oil sufficiently significant to make appreciable impacts on the US’s oil reserves and associated incidental CO 2 sequestration. The objective of this study is to improve our understanding the impact of heterogeneous and low oil saturations, in brownfield ROZs, on the effectiveness of water alternating gas (WAG) injection strategies. ROZs occur in the Permian Basin and elsewhere, and operators are using CO 2 injection for enhanced oil recovery (EOR) in these zones. The consensus model for the formation of ROZs is that they were formed by the effect of faster regional aquifer flow, acting over millions of years. Both the magnitude of oil saturation and the spatial distribution of oil differ from water-flooded main pay zones (MPZs). To explore the most effective injection strategies, we conducted simulations of CO 2 injection into synthetic geologic reservoirs. These simulations focused on injection into reservoirs subject to either man-made waterflooding or long-term natural waterflooding. By exploring the impact of varying: oil saturation; well patterns; reservoir heterogeneity; and permeability anisotropy, we attempt to quantify the factors that most influence the effectiveness of WAG injection. WAG ratios (the ratio of injected water and CO 2 , in reservoir volumes) of interest are those that either minimize the net CO 2 utilization ratios or maximize oil production rates. In general, the most effective WAG ratios for ROZs, are consistently less than those observed undergoing CO 2 injection in the same geologic reservoir models after traditional (man-made) waterflooding. This work demonstrates that most favorable WAG ratios for oil production in ROZs are different from those in traditional MPZs because of oil saturation differences. Therefore, CO 2 injection into both zones or directly copying WAG injection designs from MPZs to ROZs might not maximize oil production.

04 OIL SHALES AND TAR SANDS↗
NASA NTRS1974

Martian climate - An empirical test of possible gross variations

There appears to be evidence for a cyclic behavior of the Martian climate in which the surface pressure periodically reaches values compatible with the flow of water in equatorial regions on the planet. A relatively simple test of such hypotheses is pointed out. The premise on which cyclic models are based is that a substantial reservoir of volatils exist in frozen form at one or both poles. The proposed test involves a determination of the relative abundances of neon and argon isotopes. The required measurements may be made after the soft landing next February of Soviet spacecraft presently en route to the planet.

Owen, T.↗
NASA NTRS2011

Enceladus and Europa: How Does Hydrothermal Activity Begin at the Surface?

The question of how the surface hydrothermal activity (e.g., eruptive plumes and heat flow) is initiated can be addressed within the frame-work of our "Perrier Ocean" model. This model delivers the necessary heat and chemicals to support the heat flow and plumes observed by Cassini in Enceladus' South Polar Region. The model employs closed-loop circulation of water from a sub-surface ocean. The ocean is the main reservoir of heat and chemicals, including dissolved gases. As ocean water moves up toward the surface, pressure is re-duced and gases exsolve forming bubbles. This bub-bly mixture is less dense than the icy crust and the buoyant ocean-water mixture rises toward the surface. Near the surface, heat and chemicals, including some volatiles, are delivered to the chambers in which plumes form and also to shallow reservoirs that keep the surface ice "warm". (Plume operations, per se, are as described by Schmidt et al. and Postberg et al. and are adopted by us.) After transferring heat, the water cools, bubbles contract and dissolve, and the mixture is now relatively dense. It descends through cracks in the crust and returns to the ocean. Once the closed-loop circulation has started it is self-sustaining. Loss of water via the erupting plumes is relatively negligible compared to the amount needed to maintain the heat flow.We note that the activity described herein for the the "Perrier-Ocean" model could, a priori, apply to all small icy bodies that sheltered an interior ocean at some point in their history.

Enceladus↗
DOE OSTI2023

User’s Manual for Seal_Flux: A Seal Barrier Reduced-Order Model (Update)

This report provides a brief description on the use of the Seal_Flux computer program developed as part of the effort to quantify the risk of geologic storage of carbon dioxide (CO 2 ) under the U.S. Department of Energy’s (DOE) National Risk Assessment Partnership (NRAP). The Seal_Flux code simulates the flow of CO 2 through a low permeability rock horizon or seal formation overlying the storage reservoir into which CO 2 is injected. A two-phase, relative permeability approach with Darcy’s law is used for one-dimensional (1D) flow computations of CO 2 through the horizon in the vertical direction. The code also allows the simulation of time-dependent processes that can influence such flow. However, as part of its design, Seal_Flux is what can be termed a “reduced-order model” (ROM) and is not intended as a full-functioning flow code. The theory and simulation in the code is streamlined and directed towards the implementation of Monte Carlo risk analyses of CO 2 transport or as termed in this context as “leakage.” While presented in this report as a stand-alone tool, the Seal_Flux code is intended to function in the future as one of several models as part of an integrated, systems-level model of CO 2 storage performance. Finally, the code is written in Python 3.10 to provide an open framework for further development by others and to assist in linking the code with other modules in an integrated assessment model.

58 GEOSCIENCES↗
DOE OSTI2025

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗
DOE OSTI2020

GEOPHIRES Simulations for Deep Direct Use (DDU) Projects

This folder contains the GEOPHIRES codes and input files for running the base case scenarios for the six deep direct-use (DDU) projects. The six DDU projects took place during 2017-2020 and were funded by the U.S. Department of Energy Geothermal Technologies Office. They investigated the potential of geothermal deep direct-use at six locations across the country. The projects were conducted by Cornell University, West Virginia University (WVU), University of Illinois (U of IL), Sandia National Laboratory (SNL), Portland State University (PSU), and National Renewable Energy Laboratory (NREL). Four projects (Cornell, WVU, U of IL, SNL) investigated geothermal for direct heating of a local campus or community, the project by PSU considered seasonal subsurface storage of solar heating, and the NREL project investigated geothermal heating for turbine inlet cooling using absorption chillers. To allow comparison of techno-economic results across the six DDU projects, GEOPHIRES simulations were set up and conducted for each project. The GEOPHIRES code was modified for each project to simulate the local application and incorporate project-specific assumptions and results such as reservoir production temperature or financing conditions. The base case input file is included which simulates the base case conditions assumed by each project team. The levelized cost of heat (LCOH) is calculated and matches the base case LCOH reported by the project teams.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗
DOE OSTI2023

Efficient Optimization of Energy Recovery From Geothermal Reservoirs With Recurrent Neural Network Predictive Models

Improving the long-term energy production performance of geothermal reservoirs can be accomplished by optimizing field development and management plans. Reliable prediction models, however, are needed to evaluate and optimize the performance of the underlying reservoirs under various operation and development strategies. In traditional frameworks, physics-based simulation models are used to predict the energy production performance of geothermal reservoirs. However, detailed simulation models are not trivial to construct, require a reliable description of the reservoir conditions and properties, and entail high computational complexity. Data-driven predictive models can offer an efficient alternative for use in optimization workflows. This paper presents an optimization framework for net power generation in geothermal reservoirs using a variant of the recurrent neural network (RNN) as a data-driven predictive model. The RNN architecture is developed and trained to replace the simulation model for computationally efficient prediction of the objective function and its gradients with respect to the well control variables. The net power generation performance of the field is optimized by automatically adjusting the mass flow rate of production and injection wells over 12 years, using a gradient-based local search algorithm. Two field-scale examples are presented to investigate the performance of the developed data-driven prediction and optimization framework. Furthermore, the prediction and optimization results from the RNN model are evaluated through comparison with the results obtained by using a numerical simulation model of a real geothermal reservoir.

15 GEOTHERMAL ENERGY↗