Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir simulation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of axial tensile well deformation during reservoir compaction in offshore unconsolidated methane hydrate-bearing formation

Sand production encountered in the 2013 offshore field gas production tests at the Nankai Trough, Japan, could be attributed to well failure during reservoir compaction. In this study, well integrity under various reservoir compaction patterns for the Nankai Trough case is examined using a well-formation finite element model. Here, the modelling details include the inclusion of a cement sheath as well as the modelling of construction processes (such as cement shrinkage). Well elongation in the overburden layer becomes significant when the reservoir subsidence is localized near the wellbore under large depressurization. Results show that the maximum plastic deviatoric strain level in the cement could reach 0.7% when the maximum reservoir subsidence reaches 0.85 m and cement shrinkage is limited. When cement shrinkage rises to 0.75%, the maximum plastic deviatoric strain increases to 2.4% as the cement accumulates additional plastic strain during shrinkage due to its deformation being constrained by the casing. In order to prevent the cement from failure, it might be effective to hold the pressure drawdown at a low level (e.g., several MPa) until the hydrate dissociation front advances to a certain radius from the well (e.g., a couple of tens of metres).

58 GEOSCIENCES↗
DOE OSTI2022

Coupling subsurface and above-surface models for optimizing the design of borefields and district heating and cooling systems in the presence of varying water-table depth. In: Proceedings, 46th Workshop on Geothermal Reservoir Engineering

Dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes, without flow of groundwater. On the other hand, detailed simulation tools for subsurface heat and mass transfer exist, but these fall short in simulating above-surface energy systems. To support the design and operation of such systems, we have developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the ground. For the first, we use the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, we use the TOUGH simulator. TOUGH can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. In previous work, we described the coupling of these software packages, including how time-dependent boundary conditions for the borehole walls are synchronized for use in Modelica and TOUGH. We verified that the coupled Modelica/TOUGH code produced consistent results with the original Modelica code for an idealized problem in which heat transfer was purely by conduction in a uniform geologic medium. Here, we examine less idealized problems for which TOUGH’s advanced capabilities for modeling fluid flow are required. The first problem has a shallow vadose zone, and the second problem has a thicker vadose zone with a water-table depth that varies in time, which requires a fine vertical grid discretization for the TOUGH model.

Doughtry, Christine↗
DOE OSTI2026

Energy to Communities Program Short Introduction: Evaluation of Enhanced Geothermal System for Cornell District Heating [Slides]

Through the U.S. Department of Energy's Energy to Communities (E2C) program, NLR, other national laboratory experts, and select organizations provide Expert Match - free, short-term technical assistance to address near-term energy challenges and questions. Expert Match is for community stakeholders who have decision-making power or influence in their community but need access to additional energy expertise to inform key upcoming decisions. This Expert Match request supported Cornell University (Itaca, NY) community efforts in developing Enhanced Geothermal System (EGS) reservoir performance, by simulating fracture flow to estimate pressure drop, flow distribution, and heat extraction efficiency in fractured rock. Thermal drawdown modeling was used to assess long-term resource sustainability, while system optimization studies refined key parameters such as well configuration and flow rate to improve overall reservoir efficiency.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Creep Properties of Shale and Predicted Impact on Proppant Embedment for the Caney Shale, Oklahoma

The Caney shale is an emerging hydrocarbon play located in southwest Oklahoma, USA. Within the Caney shale exist facies which were initially dubbed “reservoir” and “ductile” based on evaluation of well logging data. While past work has shown the distinction of “brittle” and “ductile” is not mechanically justifiable according to formal definitions, here the current work shows some important differences between nominally ductile and reservoir zones. First, the “ductile” zones are more clay rich and have textural differences which can be expected to lead to differences in mechanical properties. One important impact of these differences is observed in triaxial creep experiments showing the “ductile” zones are more prone to creep deformation. Numerical simulations predict the “reservoir” zones will experience very little proppant embedment due to creep deformation of hydraulic fractures around proppant particles. On the other hand, “ductile” zones can be expected to undergo creep-driven proppant embedment leading to loss of fracture aperture ranging up to 100% loss, depending upon the spatial density of the proppant distribution. Hence, this research shows the identification of nominally “ductile” zones from well logs, while a misnomer, can be useful in finding clay-rich, creep-prone zones which will be the most prone to proppant embedment and hence vulnerable to greater production decline over time.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling study of deep direct use geothermal on the West Virginia university campus-morgantown, WV

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use (DDU) system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. Consequently, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided only the single permeability model is needed; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the prediction, therefore, heterogeneity is included in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis will be used for further economic studies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2024

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗
DOE OSTI2023

Field Laboratory for Emerging Stacked Unconventional Plays in Central Appalachia

The goal of the Field Laboratory for Emerging Stacked Unconventional Plays (ESUP) in Central Appalachia project was to investigate and characterize the resource potential for multi-play production of emerging unconventional reservoirs in Central Appalachia. Project activities included drilling, logging and coring of a vertical characterization well drilled to basement to approximately 15,000 feet in depth. The data from drilling and well logs as well recovered core samples were used by the research team to support characterization of the geology and potential pay zones within the Nora Gas Field of southwestern Virginia and the greater Central Appalachian Basin. This project was led by Dr. Nino Ripepi of the Virginia Center for Coal and Energy Research (VCCER) in close collaboration with EnerVest Operating LLC (EnerVest). VCCER is housed in the Department of Mining and Minerals Engineering at Virginia Tech and was created by an Act of the Virginia General Assembly on March 30, 1977, as an interdisciplinary study, research, information, and resource facility for the Commonwealth of Virginia. EnerVest is a top tier, low-cost oil and natural gas company with a long history of traditional operating relationships with institutional investors as well as public and private companies. This report summarizes two major research activities: 1) Log and Core Analysis and 2) Numerical Modeling. Core analysis comprises the characterization efforts of the ESUP Field Laboratory and includes a discussion of the well logs run in the deep characterization wells, an initial analysis of those, the number, size and location of cores recovered, core testing and analysis done to date as well as plans for future analysis. Numerical modeling is comprised of seven major modeling efforts, including (1) Multiphysics shale transport modeling work, (2) the developed in-house compositional simulator studies, (3) reservoir modeling studies using GEM software, (4) fracture modeling work using EFRAC3D software and (5) numerical modeling efforts with ABAQUS software, (6) Numerical modeling efforts using FLAC3D and PFC3D, and (7) Automatic Machine Learning (AutoML) Studies.

02 PETROLEUM↗
DOE OSTI2021

Marcellus Shale Energy and Environment Laboratory (MSEEL) (Final Report)

The objective of the Marcellus Shale Energy and Environment Laboratory (MSEEL) was to provide a long-term field site to develop and validate new knowledge and technology to improve recovery efficiency and minimize environmental implications of unconventional resource development. MSEEL initiated in October 2014 and completed in September 2021. Total project value was $\$29,765,067$, support from the US Department of energy totaled $\$16,608,355$ with a cost share of $\$13,156,712$ primarily from Northeast Natural Energy. This report in a departure from previous reports summarizes the overarching results and outlines the approach taken. We cover two recent efforts in machine learning and reservoir characterization and simulation. Numerous quarterly reports, public presentations and numerous external publications cover specific results by subtopic and in detail. Publications are listed in the Appendix.

03 NATURAL GAS↗
DOE OSTI2025

Foundational Dataset for Developing Large-Sample Stream Temperature Models in the Conterminous United States

This dataset provides inputs, evaluation results, and trained weights from a large-sample Long Short-Term Memory (LSTM) model designed to predict daily stream temperatures across unregulated river reaches in the conterminous United States (CONUS). It includes dynamic meteorological and hydrologic forcings, static physiographic attributes, and model outputs from cross-validation experiments spanning 300 basins. It supports reproducible modeling, direct application for new basins, and provides data suitable for integration with reservoir and river simulations under current and future climates. It contains two .zip files described below · RQ-AI_runs.zip: Model outputs from 10-fold cross-validation experiments, including observed and predicted daily stream temperatures, along with test performance metrics for water years 2017–2019. Two versions are included: 1. Model trained and validated using subbasin-area weighted dynamic features. 2. Model trained and validated using whole-basin area weighted dynamic features. · RQ-AI_inputs.zip: Collection of all formatted dynamic and static predictor datasets (meteorological, hydrologic, and physiographic features) used in model training and analysis. Detailed instructions and data structure is held at the following GitLab repository: https://code.ornl.gov/tempwise/training.

Gomez-Velez, Jesus [Oak Ridge National Laboratory ↗
NASA NTRS1984

Vectorization of a penalty function algorithm for well scheduling

In petroleum engineering, the oil production profiles of a reservoir can be simulated by using a finite gridded model. This profile is affected by the number and choice of wells which in turn is a result of various production limits and constraints including, for example, the economic minimum well spacing, the number of drilling rigs available and the time required to drill and complete a well. After a well is available it may be shut in because of excessive water or gas productions. In order to optimize the field performance a penalty function algorithm was developed for scheduling wells. For an example with some 343 wells and 15 different constraints, the scheduling routine vectorized for the CYBER 205 averaged 560 times faster performance than the scalar version.

Absar, I.↗
NASA NTRS2016

Validation of a Cephalad Fluid Shift Countermeasure

INTRODUCTION: This project will provide critical data required to objectively determine how an optimized thigh cuff could be incorporated into the NASA integrated physiological countermeasure suite. This project will determine if thigh cuffs used during simulated spaceflight impact intracranial pressure (ICP), ocular structure and function, and intraocular pressure (IOP) using state of-the-art techniques. Additionally, some of the same methods, hardware, and protocols will be employed in the present investigation to enable direct comparisons to the International Space Station (ISS) "Fluid Shifts" experiment with Chibis-Lower Body Negative Pressure (LBNP). This study will determine the temporal physiological responses of thigh cuff application and removal on ocular and cerebral variables (including invasive ICP) in a microgravity analog. Furthermore, this proposed study will determine tissue pressure distribution applied by thigh cuffs in order to improve comfort, mobility, and efficacy of the countermeasure. Our specific aim is to determine the efficacy of a novel thigh cuff device to mitigate cephalad fluid shifts. We hypothesize that a thigh cuff countermeasure employed in a microgravity analog will temporarily reverse or attenuate ocular and cerebral-volume-pressure variables, approaching normal Earth-based seated posture, the most frequent posture assumed in daily life. In addition, we hypothesize that the magnitude of fluid and pressure redistribution using a thigh cuff countermeasure may require a longer exposure time than that of Chibis-LBNP (using ground-based data from our "Fluid Shifts" project). This project directly addresses Critical Path Roadmap Risks and Questions regarding "Risk of Spaceflight-Induced Intracranial Hypertension/Vision Alterations," and IRP Gap VIIP13: We need to identify preventative and treatment countermeasures to mitigate changes in ocular structure and function and intracranial pressure during spaceflight. METHODS: Noninvasive measures and tissue pressure distributions beneath thigh cuffs The objectives of this study are to: 1) determine the distribution of skin surface pressures beneath the advanced thigh cuff in ten subjects, 2) calibrate the built-in pressure measurement system of the advanced thigh cuff using an industry standard device, and 3) collect subjective feedback and data on the new cuff design to allow for further adjustments prior to invasive studies. A Tekscan Industrial Sensing (I-Scan) system will measure the pressure distribution of the advanced thigh cuff against the skin. In addition, we will measure blood pooling in the thigh and record the circumference of the thigh using Hokanson strain gauge plethysmography. The advanced thigh cuff will be adjusted to obtain a skin contact pressure of 30-50 mmHg as visualized on the Tekscan system. The built-in advanced thigh cuff pressure monitor will be recorded simultaneously to allow direct comparison to the Tekscan measurements. The volunteer will then remove the thigh cuff and remain at rest for five minutes with no legging applied. The thigh cuff will be donned again and pressure measurements will be taken in the same manner for up to 10 repetitions to show reproducibility of pressure after donning. At the conclusion of the study, subjects will be asked to flex their knee, stand, walk, and sit with the thigh cuff activated. During each of these maneuvers the subject will rate their pain/comfort using a modified Borg scale. Effect of thigh cuffs on ICP during simulated microgravity Ommaya reservoir patients will be recruited from the John Wayne Cancer Institute. Ommaya reservoirs provide safe and direct access for the measurement of ICP. Subjects will be instrumented for continuous blood pressure, ECG, and invasive ICP measures. The subjects will be positioned in the upright sitting posture for a 10-minute stabilization period. After the 10-minute stabilization period, imaging measures [ICP, Optical Coherence Tomography, IOP, ocular and vascular ultrasound] will be performed. Following baseline seated measures, the subject will be positioned randomly in the supine, 15deg head-down-tilt, and 15deg head-down-tilt with thigh cuffs and measures repeated. DISCUSSION: Tests to down-select thigh cuff designs will occur in early 2016. Invasive ICP and noninvasive eye imaging tests will begin in spring 2016. Supported by NSBRI through NCC 9-58.

Macias, B.↗
NASA NTRS2019

Measuring the D/H Ratios of Exoplanets and Brown Dwarfs

The relative abundance of deuterium and hydrogen is a potent tracer of planet formation and evolution. Jupiter and Saturn have protosolar atmospheric D/H ratios, a relic of substantial gas accretion from the nebula, while the atmospheres of Neptune and Uranus are enhanced in D by accretion of ices into their envelopes. For terrestrial planets, D/H ratios are used to determine the mechanisms of volatile delivery and subsequent atmosphere loss over the lifetime of the planet. Planets and brown dwarfs more massive than ~13 M(J) quickly fuse their initial D reservoir. Here, we simulate spectra for giant exoplanets and brown dwarfs (2 M(Neptune) to ~10 M(Jupiter)) from T(eff) = 200–1800 K including both CH3D and HDO to determine the observability of these dominant deuterium isotopologues in mid-infrared thermal emission spectra. Colder objects have stronger molecular features in their spectra, due to the temperature dependence of molecular cross sections. CH3D is easier to observe than HDO at all temperatures considered, due to the strength of its absorption bands and locations of features at wavelengths with few other strong absorption features. We predict that for nearby cool brown dwarfs, the CH3D abundance can be measured at high signal to noise with the James Webb Space Telescope (JWST); for objects from 200 to 800 K closer than 10 pc, a protosolar D/H ratio would be readily observable in 2.5 hr. Moderately young Jupiter-mass planets (100–300 Myr) and young Neptunes (10 Myr) may be discovered with JWST and provide the best targets for detecting deuterium on an exoplanet in the coming decade. Future telescope designs should consider the importance of isotopes for understanding the formation and evolution of planetary atmospheres.

Caroline V. Morley↗
DOE OSTI2022

Effect of Initial Water Saturation on The Performance of Fracturing Fluids With and Without Polyallylamine under Simulated EGS Conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO2-reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO2, can be used as a less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO2 fracturing fluids outperformed other fluids such as water, CO2, and CO2/water in generating large fractures with less fluid consumed. The objective of this work is to study the effect of initial water saturation of rock on the performance of StimuFrac fluid in ½ foot cubic rock samples and under representative EGS pressure/temperature conditions using cyclic and constant flow rate injection strategies. The fracturing results are compared with results using different fracturing fluids in terms of controlling fracture propagation rates, fracture hydraulic conductivity, breakdown pressures and volumes of fluids required. Methods/Procedures/Process: In all tests, water was initially injected into the rock to increase the water saturation before the fracturing processes to simulate actual geothermal reservoir conditions. For the cyclic injection, one complete cycle consisted of (1) a PAA slug (or water slug) injection followed by (2) CO2 injection to initiate the fracture. In the second step of the first cycle, the pressure of CO2 is increased until a maximum pressure is reached (fracture is initiated at this moment), and then the injection of CO2 is allowed for another 30 seconds to propagate the fracture. Then, the two-step cycle of PAA followed by CO2 injection (up to 2-4 mL/min) was continued. For the constant flow rate injection strategies, the initial water saturation was increased by injecting water at 1000 psi and 200°C for three days. After that, an initial slug of water, CO2, or PAA was injected and then fracturing was initiated and propagated by injecting CO2 at a constant flow rate. Applications/Significance/Novelty: The results of this study suggest that water saturation, especially near the wellbore region, will significantly affect the fracturing fluid transmission into the rock porous media by changing the relative permeability of CO2 or water, thus affecting the fracture initiation and propagation. In this study, fracturing with cyclic injection or constant flow rate injection methods were performed using three different kinds of fluids systems. These fluids are water, CO2, or CO2 with PAA. Splitting the rock samples in half after fracturing reveals that the fracture propagation is significantly limited under these high water saturated conditions compared to dry initial conditions: The fractures propagate less than 1/3 length of the distance from the wellbore to rock surface, and in some cases no fracture is generated. This may be caused by the fact that leak-off is dominating the fracturing process and the injected fluid flow rate is not high enough to overcome the leak-off even under high flow rate injection conditions. Additionally, CO2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate pressure gradients away from the near-wellbore region.

Jian, Guoqing↗
DOE OSTI2020

Optimal pumping schedule with high-viscosity gel for uniform distribution of proppant in unconventional reservoirs

Hydrocarbon recovery from a hydraulic-fractured unconventional reservoir can be enhanced by replacing slickwater with non-Newtonian high-viscosity gel as a fracturing fluid. However, in the literature, the role of gel rheology on hydrocarbon production from a fractured reservoir is not very well understood. This lack of understanding makes it difficult for the industry to choose appropriate viscosity parameters of fracturing gels. Motivated by these limitations, first, we have developed a high-fidelity model of non-Newtonian fluid flow to simulate fracture propagation, fluid leak-off, and gel flowback processes to predict the oil production from a hydraulic fractured reservoir. Simulation results of the high-fidelity model show how fracture geometry, formation damage, and gel cleanup process varies with the fracturing fluid rheology. Therefore, the developed model is used to understand the effect of gel rheology on hydrocarbon production via a sensitivity analysis, which provides a set of rheological parameters that maximize the cumulative amount of oil production. Finally, utilizing this optimal high-viscosity gel, a model predictive controller is designed to obtain an optimal pumping schedule necessary for producing a target fracture geometry and proppant concentration. The closed-loop simulation results demonstrate that the obtained pumping schedule successfully achieved the desired fracture geometry and proppant concentration at the end of pumping, leading to the maximum oil production from an unconventional reservoir.

13 HYDRO ENERGY↗
DOE OSTI2024

Multi-fidelity Fourier neural operator for fast modeling of large-scale geological carbon storage

Deep learning-based surrogate models have been widely applied in geological carbon storage (GCS) problems to accelerate the prediction of reservoir pressure and CO2 plume migration. Large amounts of data from physics-based numerical simulators are required to train a model to accurately predict the complex physical behaviors associated with this process. In practice, the available training data are always limited in large-scale 3D problems due to the high computational cost. Therefore, we propose to use a multi-fidelity Fourier neural operator (FNO) to solve large-scale GCS problems with more affordable multi-fidelity training datasets. FNO has a desirable grid-invariant property, which simplifies the transfer learning procedure between datasets with different discretization. Here, we first test the model efficacy on a GCS reservoir model being discretized into 110 k grid cells. The multi-fidelity model can predict with accuracy comparable to a high-fidelity model trained with the same amount of high-fidelity data with 81% less data generation costs. We further test the generalizability of the multi-fidelity model on a same reservoir model with a finer discretization of 1 million grid cells. This case was made more challenging by employing high-fidelity and low-fidelity datasets generated by different geostatistical models and reservoir simulators. We observe that the multi-fidelity FNO model can predict pressure fields with reasonable accuracy even when the high-fidelity data are extremely limited. The findings of this study can help for better understanding of the transferability of multi-fidelity deep learning surrogate models.

58 GEOSCIENCES↗
NASA NTRS2012

The Connection between Model Performance on the CCMVal Transport Diagnostics and Simulated Sensitivity of Ozone to Chlorine Change

The chemistry climate models that contributed simulations for past and future ozone evolution to the 2010 Scientific Assessment of Ozone Depletion were subject to extensive evaluation by the SPARC (Stratospheric Processes and their Role in Climate) CCMVal (Chemistry-Climate Model Validation) activity. The sensitivity of ozone to changes in composition and climate varies among the models, but the relationship between these variations and the model evaluations of CCMVal is not obvious. We have learned that the transport evaluation can be used to interpret the comparisons between observed and simulated columns of chlorine reservoirs, hydrochloric acid (HCl) and chlorine nitrate (ClONO2); these comparisons were part of the CCMVal evaluation of chemistry. The simulations with best performance on the transport diagnostics most faithfully reproduce the evolution and seasonal variation of the chlorine reservoirs as observed at NDACC (Network for Detection of Atmospheric Composition Change) stations (NyAlesund 78.9N, Kiruna 67.8N, Harestua 60.2N, Jungfraujoch 46.6N, Toronto 43.6N, Kitt Peak 31.9N, Izana 28.3N, Mauna Loa 19.5N, Lauder 45S and Arrival Heights 77.8S). In the simulations, the HCl in the lower stratosphere depends on total inorganic chlorine (Cly) and partitioning between HCl and ClON02. Total inorganic chlorine depends on the fractional release of chlorine from source gases, and ratio of ClON02 to HCl is inversely dependent on methane and varies quadratically with ozone. Simulated HCl from various models may agree with observations even though Cly is in error, partitioning is in error, or both. Simulated ozone sensitivity to chlorine is shown to be greater for models that produce larger values of chlorine nitrate for background chlorine levels, and vice versa. Comparisons with the NDACC data show why the models with 'best' transport have similar sensitivity to chlorine change. The realistic evolution of the simulated HCl and ClONO2 columns suggests realistic levels of Cly in the lower atmosphere. In addition, the wide range values for the sensitivity of ozone to chlorine obtained from the CCMVal simulations is explained by the wide range in lower atmospheric columns of ClONO2 and the concomitant wide range of levels for chlorine monoxide.

Douglass, Anne↗