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 127 records · Page 7

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗

EGS Collab Experiment 1: DNA tracer data on transport through porous media

This submission contains DNA tracer data that supports the analysis and conclusions of the publication, "DNA tracer transport through porous media -The effect of DNA length and adsorption." https://doi.org/10.1029/2020WR028382. This experiment used DNA as an artificial reservoir tracer. Groundwater tracing is an effective way to identify fluid flow pathways and estimate hydrogeologic properties, which are important premises for building reliable hydrological models for transport predictions or contamination mitigations. The objective of this study was to understand the effect of DNA length (i.e., number of base pairs for dsDNA) and adsorption on DNA transport.

15 GEOTHERMAL ENERGY↗

A comparison of economic evaluation models as applied to geothermal energy technology

Several cost estimation and financial cash flow models have been applied to a series of geothermal case studies. In order to draw conclusions about relative performance and applicability of these models to geothermal projects, the consistency of results was assessed. The model outputs of principal interest in this study were net present value, internal rate of return, or levelized breakeven price. The models used were VENVAL, a venture analysis model; the Geothermal Probabilistic Cost Model (GPC Model); the Alternative Power Systems Economic Analysis Model (APSEAM); the Geothermal Loan Guarantee Cash Flow Model (GCFM); and the GEOCOST and GEOCITY geothermal models. The case studies to which the models were applied include a geothermal reservoir at Heber, CA; a geothermal eletric power plant to be located at the Heber site; an alcohol fuels production facility to be built at Raft River, ID; and a direct-use, district heating system in Susanville, CA.

Ziman, G. M.↗

mosartwmpy sample input data; 1980 - 1985

Sample input data spanning the years 1980-1985 for running the `mosartwmpy` water routing and management model: https://github.com/IMMM-SFA/mosartwmpy. More information is available in the README.md. Changelog: v0.0.7: new reservoir files with corrected locations based on an exhaustive review by Dan Broman for the 9505 project; corresponding long term mean demand and flow files updated based on VIC4 simulations using Daymet forcing. v0.0.6: updates reservoir parameters to indicate which reservoirs should follow generic operating rules and which should follow ISTARF data-driven operating rules. v0.0.5: adds supporting files for running the Farmer Agent Based Model of adaptive water demand.

Rexer, Emily↗

Repository of HydroSMADE: Hydropower Site-level Monthly Availability Data Ensemble for 1950-2100 at Existing and Potential Global Sites

This repository presents HydroSMADE—Hydropower Site-level Monthly Availability Data Ensemble, a new open dataset that provides monthly hydropower availability for 1,593 existing and 124,333 potential sites worldwide over the period 1950–2100. The dataset is generated by using a global hydrologic model (Xanthos) with explicit representation of hydropower operation. Specifically, HydroSMADE distinguishes between storage and diversion sites, applies optimized operating rules, and incorporates site-specific characteristics such as generation capacity, maximum turbine flow, and reservoir storage. Driven by bias-corrected meteorological inputs, the data is provided for 30 alternative future scenarios. The scenarios consist of the full factorial combination of three standard CMIP6 atmospheric forcing pathways (SSP1-2.6, SSP3-7.0, and SSP5-8.5) and ten CMIP6 General Circulation Models (GCMs): GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, EC-Earth3, CanESM5, MIROC6, CNRM-ESM2-1, UKESM1-0-LL, and CNRM-CM6-1. The repository contains a total of 122 files: a text file (readme.txt) containing a brief description of the included data, a CSV file containing site attributes, and the remaining 120 files (in CSV) containing site-level monthly hydropower availability. Example Jupyter Notebooks to explore the HydroSMADE dataset are available on GitHub at https://github.com/kamal0013/HydroSMADE More details on the methods and technical validation of HydroSMADE are available in the following paper by the same authors: Chowdhury, A. K., Abeshu, G. W., Zhao, M., Wild, T. B., Hassan, N., Ying, Z., Kim, G. J., Matthew, B., Jonathan, L., & Li, H.-Y. (Submitted). Hydropower Site-level Monthly Availability Data Ensemble for 1950-2100 at Existing and Potential Global Sites.

Existing and Potential Sites↗

A combined experimental and modelling study of granite hydrothermal alteration

Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.

58 GEOSCIENCES↗

Deposition, Diagenesis, and Sequence Stratigraphy of the Pennsylvanian Morrowan and Atokan Intervals at Farnsworth Unit

Farnsworth Field Unit (FWU), a mature oilfield currently undergoing CO2-enhanced oil recovery (EOR) in the northeastern Texas panhandle, is the study area for an extensive project undertaken by the Southwest Regional Partnership on Carbon Sequestration (SWP). SWP is characterizing the field and monitoring and modeling injection and fluid flow processes with the intent of verifying storage of CO2 in a timeframe of 100–1000 years. Collection of a large set of data including logs, core, and 3D geophysical data has allowed us to build a detailed reservoir model that is well-grounded in observations from the field. This paper presents a geological description of the rocks comprising the reservoir that is a target for both oil production and CO2 storage, as well as the overlying units that make up the primary and secondary seals. Core descriptions and petrographic analyses were used to determine depositional setting, general lithofacies, and a diagenetic sequence for reservoir and caprock at FWU. The reservoir is in the Pennsylvanian-aged Morrow B sandstone, an incised valley fluvial deposit that is encased within marine shales. The Morrow B exhibits several lithofacies with distinct appearance as well as petrophysical characteristics. The lithofacies are typical of incised valley fluvial sequences and vary from a relatively coarse conglomerate base to an upper fine sandstone that grades into the overlying marine-dominated shales and mudstone/limestone cyclical sequences of the Thirteen Finger limestone. Observations ranging from field scale (seismic surveys, well logs) to microscopic (mercury porosimetry, petrographic microscopy, microprobe and isotope data) provide a rich set of data on which we have built our geological and reservoir models.

04 OIL SHALES AND TAR SANDS↗

Aromatic amino acid metabolism and active transport regulation are implicated in microbial persistence in fractured shale reservoirs

Abstract Hydraulic fracturing has unlocked vast amounts of hydrocarbons trapped within unconventional shale formations. This large-scale engineering approach inadvertently introduces microorganisms into the hydrocarbon reservoir, allowing them to inhabit a new physical space and thrive in the unique biogeochemical resources present in the environment. Advancing our fundamental understanding of microbial growth and physiology in this extreme subsurface environment is critical to improving biofouling control efficacy and maximizing opportunities for beneficial natural resource exploitation. Here, we used metaproteomics and exometabolomics to investigate the biochemical mechanisms underpinning the adaptation of model bacterium Halanaerobium congolense WG10 and mixed microbial consortia enriched from shale-produced fluids to hypersalinity and very low reservoir flow rates (metabolic stress). We also queried the metabolic foundation for biofilm formation in this system, a major impediment to subsurface energy exploration. For the first time, we report that H. congolense WG10 accumulates tyrosine for osmoprotection, an indication of the flexible robustness of stress tolerance that enables its long-term persistence in fractured shale environments. We also identified aromatic amino acid synthesis and cell wall maintenance as critical to biofilm formation. Finally, regulation of transmembrane transport is key to metabolic stress adaptation in shale bacteria under very low well flow rates. These results provide unique insights that enable better management of hydraulically fractured shale systems, for more efficient and sustainable energy extraction.

04 OIL SHALES AND TAR SANDS↗

An Integrated Hydroclimatic Assessment of Future Reservoir and Hydropower Operations in the U.S.

The engineering of rivers by dams is a formative feature of human-nature systems and the interconnectivity of water, energy, and the climate. Sufficient and broad-based representations of dams in large-scale hydrological models prove essential to mapping their extensive regulation of river flow and biogeochemistry and gauging climate-linked provisions, including freshwater supply and hydropower. We present an integrated modeling framework to investigate future streamflow and hydropower generation in the Contiguous U.S. (1990–2075), leveraging an ensemble of six downscaled and bias-corrected General Circulation Models (GCMs) from the high-end SSP585 scenario of the CMIP6. To achieve this, we develop a reservoir operations and parameterization scheme for 1,384 dams in a high-resolution river network, including simulated hydropower generation for 326 dams. For the GCM ensemble mean, we simulate a widespread increase in regulated streamflow into the late-century (11% annual and 17% in winter for the dam median) with region-specific changes in summer streamflow that feature prominent declines in the Northwest (−7%). Mediation by reservoirs is shown to dampen intra-annual streamflow changes, delivering additional summer releases that partially mitigate declining flows. Total hydropower generation is projected to increase modestly (+3%), with boosted generation in the winter (+9%) and spring (+5%) offsetting declined summer generation (−3.4%), suggesting strong adaptation potential for hydropower in the future energy portfolio. Further analysis reveals that the choice of GCM, particularly in western regions, has significant bearing on projected streamflow and hydropower changes.

13 HYDRO ENERGY↗

Accelerating Multiphase Simulations With Denoising Diffusion Model Driven Initializations

This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics‐based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical properties of interest. Pore‐scale simulations enhance our understanding of applications such as assessing hydrogen and storage efficiency in underground reservoirs. Nevertheless, they are computationally expensive and the presence of non‐unique solutions can require multiple simulations within a single geometry. To overcome the computational cost hurdle, we propose a method that couples generative diffusion models and physics‐based simulations. While training the data‐driven model, we simultaneously generate initial conditions and perform physics‐based simulations using these. This integrated approach enables us to receive real‐time feedback on a single compute node equipped with both CPUs and GPUs. By efficiently managing these processes within a single compute node, we can continuously monitor performance and halt training once the model meets the specified criteria. To test our model, we generate realizations in a real Berea sandstone fracture which shows that our technique is up to 4.4 times faster than commonly used flow simulation initializations.

36 MATERIALS SCIENCE↗

Impact of Rheology Models on Horizontal Well Polymer Flooding in a Heavy Oil Reservoir on Alaska North Slope: A Simulation Study

Polymer rheology can have either a positive or a negative effect on polymer flooding performance under varied circumstances. Many researchers have studied the effect of polymer rheology in a vertical well, but no field scale studies have been conducted to investigate whether polymer rheology is beneficial to polymer flooding in heavy oil reservoirs developed by horizontal wells. In this paper, we conducted a numerical simulation study to examine the effect of HPAM polymer rheology on a polymer flooding pilot, which is the first-ever project conducted on a heavy oil reservoir from Alaska North Slope (ANS) developed by horizontal wells. Three rheology types were considered in the study including the apparent viscosity measured during coreflooding of using a HPAM polymer, the bulk viscosity measured with a viscometer, and a Newtonian flow model. The results suggest that using the bulk viscosity in simulation underestimates the conformance control and the water-oil-ratio reduction capability of the HPAM polymer solution. When the apparent viscosity is used, the incremental oil and sweep were largely increased, and the optimal recovery period of polymer flooding was extended greatly, especially for the heterogeneous formations. Therefore, the rheology type of polymer plays a significant role in the incremental oil recovery and injection profile of the horizontal well system given the pilot test conditions. This study has provided practical guidance to field operators for the ongoing polymer flooding pilot on ANS and will also provide valuable information for other polymer projects conducted in similar conditions.

Leng, Jianqiao↗

History Matching and Prediction of a Polymer Flood Pilot in Heavy Oil Reservoir on Alaska North Slope

The first-ever polymer flood pilot to enhance heavy oil recovery on Alaska North Slope is ongoing. After more than 3 years of polymer injection, significant benefit has been observed from the decrease in water cut from 65% to less than 15% in the project producers. The primary objective of this study is to develop a robust history-matched reservoir simulation model capable of predicting future polymer flood performance. In this work, the reservoir simulation model has been developed based on the geological model and available reservoir and fluid data. In particular, four high transmissibility strips were introduced to connect the injector-producer well pairs, simulating short-circuiting flow behavior that can be explained by viscous fingering and reproducing the water cut history. The strip transmissibilities were manually tuned to improve the history matching results during the waterflooding and polymer flooding periods, respectively. It has been found that higher strip transmissibilities match the sharp water cut increase very well in the waterflooding period. Then the strip transmissibilities need to be reduced with time to match the significant water cut reduction. The viscous fingering effect in the reservoir during waterflooding and the restoration of injection conformance during polymer flooding have been effectively represented. Based on the validated simulation model, numerical simulation tests have been conducted to investigate the oil recovery performance under different development strategies, with consideration for sensitivity to polymer parameter uncertainties. The oil recovery factor with polymer flooding can reach about 39% in 30 years, twice as much as forecasted with continued waterflooding. Besides, the updated reservoir model has been successfully employed to forecast polymer utilization, a valuable parameter to evaluate the pilot test’s economic efficiency. All the investigated development strategies indicate polymer utilization lower than 3.5 lbs/bbl in 30 years, which is less than that of the same polymer used in a polymer pilot in Argentina.

Wang, Xindan↗

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗

Heavy Oil Polymer EOR in the Challenging Alaskan Arctic - It Works!

Under co-sponsorship of the US Department of Energy and Hilcorp Alaska LLC the first ever polymer field pilot commenced on 8/28/2018 in the Schrader Bluff heavy oil reservoir at the Milne Point Field on Alaska North Slope (ANS). The primary objective of the pilot is to prove the efficacy of polymer Enhanced Oil Recovery (EOR) to unlock the vast heavy oil resources on ANS. More than two and half years after startup, the polymer injection, supporting laboratory experiments and simulation studies steadily continue. The pilot started injecting hydrolyzed polyacrylamide (HPAM), at a concentration of 1,750 ppm to achieve a target viscosity of 45 cP, into the two horizontal injectors in the flood pattern. Production is monitored in the two horizontal producers. Based on laboratory measurements of polymer viscosity at reservoir conditions, the team decided to reduce polymer concentration to 1,200 ppm since July 2020 in an effort to control injection pressure and optimize polymer utilization. Quality control (QC) on the field ensures uniform polymer solution properties. Representative rock and fluid systems and test conditions are utilized in the corefloods on polymer retention, effect of injection water salinity, polymer loading, and their combinations on oil recovery. A history matched reservoir simulation model for forecasting oil recovery was developed on the basis of all the available field data. Field concerns related to the post-polymer breakthrough impact on flow assurance is addressed via specialized laboratory tests. Notwithstanding early operational disruptions and hydration issues, continuous polymer injection in both injectors has been achieved. To date, 950,000 lbs of polymer or 2 million barrels of polymer solution, equating to 13% of total pore volume (PV), has been placed in the flood pattern, serving as an effective indicator of adequate polymer injectivity. So far, the success of polymer EOR is evident from drastically reduced water cut in the producers, an estimated incremental 1,000 bopd over waterflood, and a favorable polymer utilization of 1.7 lbs/barrel of incremental oil. Polymer breakthrough was observed 26 months after the start of polymer injection. Main observations from corefloods are unusually high polymer retention values in some cases and a positive response to low salinity water. Although the heterogeneity in the flood pattern and exceptionally low water cut pose some challenges, persistent novel and justifiable simulation approaches have resulted in a robust history matched model. Experimental results on produced fluids treatment provide operational guidance for improved oil-water separation and mitigation of heater tube fouling. The ongoing success of the pilot is a key indicator of bringing the team ever closer to meeting the project’s primary objective. The encouraging results of the pilot is one of the drivers that has provided the impetus to apply polymer EOR throughout the Milne Point Field, which would increase oil recovery and extend the economic life of the Trans Alaska Pipeline System. The scientific knowledge, including the many lessons learned from this pilot also has referential value for other potential heavy oil EOR projects throughout the world. The collected data, operational lessons learned, and the overall success of the pilot are summarized in the paper.

None, None↗

Fundamentals of chemistry modeling applicable to a vectorized particle simulation

This paper describes the fundamentals of extending the vectorized particle simulation method derived by Baganoff and McDonald (1990), McDonald and Baganoff (1988), and McDonald (1989) for modeling chemically reacting flows. Details of reaction mechanics per reaction are presented, with particular attention given to the quantum nature of the vibrational mode. The models of reactive flows developed here were verified through a simulation of a superheated diatomic gas relaxing thermochemically to equilibrium in a reservoir.

Haas, Brian L.↗

A geospatial model for the analysis of time-dependent land subsidence induced by reservoir depletion

Land subsidence due to fluid depletion is an outcome of physical processes operating across a wide range of time and length scales. Although geomechanical models are crucial to simulate reservoir compaction and predict its long-term fate, their use across large regions often bears prohibitive computational costs. To overcome this obstacle, this paper proposes a simplified modelling framework consisting of (i) a near-field numerical solver simulating the coupling between fluid flow and rock deformation with reference to a simplified one-dimensional geometry and (ii) a far-field geospatial algorithm mapping ground settlements across a region through the superposition of poroelastic computations at multiple wells. The model computes the delay between depletion history and reservoir compaction in proximity of a producing well by assuming basal depletion of a fluid-saturated deformable disk, while the Geertsma solution of nucleus of strain is used to extrapolate the impact of such time-varying reservoir compaction around the well. This approach has been used to back-analyze the spatio-temporal progression of subsidence at the Groningen gas field. The results are compared against measurements collected over 50 year of production at 25 benchmark locations scattered over an area of 900 km 2 . It is shown that coupled simulations based on average values of rock compressibility and permeability lead to nonlinear trends of subsidence evolution in good agreement with field measurements, while uncoupled analyses overpredict settlements by more than 70%. Lastly, synthetic forecasts based on different rates of depletion were provided. The results suggest that slower depletion rates lead to lower subsidence at a given time, and that residual subsidence may continue to develop for several decades after interruption of production activities.

58 GEOSCIENCES↗

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗