Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15

Transient flow modeling in fractured media using graphs

In this work, we describe a method to simulate transient fluid flows in fractured media using an approach based on graph theory. Our approach builds on past work where the graph-based approach was successfully used to simulate steady-state fluid flows in fractured media. We find a mean computational speedup of the order of 1400 from an ensemble of a 100 discrete fracture networks in contrast to the O(10 4 ) speedup that was obtained for steady-state flows earlier. However, the transient flows considered here involve an additional degree of complexity that was not present in the steady-state flows considered previously with a graph-based approach, that of time marching and solution of the flow equations within a time-stepping scheme. We verify our method with an analytical test case and demonstrate its use on a practical problem related to fluid flows in hydraulically fractured reservoirs. By enabling the study of transient flows, we create an opportunity for a wide set of possibilities where a steady-state approximation is not sufficient, such as the example motivated by hydraulic fracturing that we present here. This work validates the concept that graphs are able to reliably capture the topological properties of the fracture network and serve as effective surrogates in an uncertainty-quantification framework.

58 GEOSCIENCES↗

Reactive chemical transport simulations of geologic carbon sequestration: Methods and applications

Chemical reaction simulations are considerably used to quantitatively assess the long-term geologic carbon sequestration (GCS), such as CO 2 sequestration capacity estimations, leakage pathway analyses, enhanced oil recovery (EOR) efficiency studies, and risk assessments of sealing formations (caprocks), wellbores, and overlying underground water resources. All these require a deep understanding of the CO 2 - associated chemical reactions. To ensure long-term, safe CO 2 sequestration in the intended formations, modeling is the only way to plausibly assess the CO 2 flow, reaction, and transport over thousands of years. This review summarizes the multiple methodologies for describing homogeneous and heterogeneous chemical reaction patterns and multiscale application examples, the recent progress and current status of chemical reaction simulations for GCS, and the impact of such simulations on geological CO2 sequestration performance. Technical gaps and future challenges are also discussed for further study. The trends and challenges of such studies include: (1) the combination of coupled chemical, mechanical, and transport processes with calibrated experiments and associated uncertainty/risk assessments; (2) enhancement of the ability to simulate detailed geophysical and geochemical equations to mimic in situ conditions; and (3) characterization of multiscale subsurface systems with detailed conceptual models and assignment of suitable boundary conditions for field-scale sequestration fields. One major gap remaining is the current lack of accurate (and scale-justified) kinetic and equilibrium chemical reaction parameters under reservoir conditions. Advanced models that couple chemical, mechanical, and transport processes with scale-justified parameters, from lab to field-scale experiments, are required for quantitative assessments of sequestration capacity and the long-term safety of GCS projects.

58 GEOSCIENCES↗

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗

Hero Carbonsafe Phase 2 Project in the Columbia River Basalt Group: Technical Program Overview

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Hermiston is located near the center of the Columbia River Basalt Group (CRBG), which is one of the largest basalt flows in the US. Basalt CO2 storage has potential advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. associated decreases in pressure and CO2 migration risks, 3. reduced long-term monitoring requirements with respect to plume tracking, 4. widespread geographic distribution and, 5. large storage potential due to thickness, porosity, and CO2 interactions with basalt. For locations such as the Pacific Northwest (PNW), Hawaii, Iceland, India and Japan, whose localities are isolated from large sedimentary basins offering conventional saline storage options, basalt may offer the only feasible option for local CO2 storage. However, mineralization/basalt storage still has many uncertainties, as there are limited field-scale assessments of CO2 storage in basalt. There are significant uncertainties hindering the effective implementation of carbon capture utilization and storage (CCUS) in basalt. These include the lack of proven storage capacities, challenges in methodologies for modeling the area of review in igneous formations, limited understanding of mineralization kinetics and timing, and uncertainties in injectivity. Additionally, the domestic availability of specialized services and drilling expertise is constrained, and existing CCUS permitting and regulatory frameworks, originally developed for conventional saline reservoirs, may not adequately address the unique requirements of basalt systems. HERO CarbonSAFE is designed to address major research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts of the CRBG, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts.

58 GEOSCIENCES↗

Modeling CO 2 plume migration using an invasion-percolation approach that includes dissolution

Geological carbon dioxide (CO 2 ) sequestration has been proposed as a viable technique to decrease effective emissions of CO 2 into the atmosphere. However, the security of this sequestration is tied to our understanding of the long-term migration of CO 2 in subsurface. The dissolution of CO 2 in the reservoir brine is one of the main long-term trapping mechanisms. However, the assumptions used in large-scale reservoir simulations usually lead to an overestimation of the dissolution volume. We propose a modified approach based on the macroscopic invasion-percolation (MIP) theory that allows the dissolution of CO 2 into brine. We used a high-resolution geological model to compare the Darcy-, modified MIP-, and classic MIP-based simulation results. We observed a significant shrinkage in the nonaqueous plume volume when dissolution is considered during the MIP simulation. In the case of Darcy-based simulation, the plume was completely trapped inside the reservoir with limited migration even after a thousand-year simulation. On the other hand, the majority of the plume migrated out of the simulated reservoir in the case of MIP. Our approach provides more realistic estimation of the dissolution volume and nonaqueous plume extent while leveraging the computational efficiency enjoyed by MIP.

54 ENVIRONMENTAL SCIENCES↗

MRCI Subtask 2.4/2.5: Regional/Subregional Analysis and Risk Assessment Final Technical Summary Report

The objective of the Midwest Regional Carbon Initiative (MRCI) project is to implement a collaborative Regional Initiative (RI) to accelerate the deployment of carbon capture, and storage (CCS) in the Midwest-Northeastern quadrant of the United States. This report is a Technical Summary report describing work performed on Tasks 2.4 (Conducting Regional/Subregional Analysis) and 2.5 (Assessing and Managing Risk) during the MRCI project. In Task 2.4, detailed numerical reservoir simulation models were developed for selected carbon storage (CS) systems identified under Task 2.1 (Battelle, 2021a) in the MRCI study area. The objective of Task 2.4 is to demonstrate a dynamic modeling methodology for evaluating the suitability of the (selected) CS systems in the MRCI region for hosting a commercial-scale storage project. In this study, CO2 injectivity was evaluated for different CS systems with an annual injection rate of 1 million metric tonnes (MMT) of CO2 considered as the minimum requirement for a commercial scale project. The objective of Task 2.5 is to assess key risks associated with storing CO2 in the different CS systems across the MRCI region and to demonstrate a method(s) for assessing these risks that may be used by developers of future CO2 storage projects in the region. The risk analysis was limited to evaluating two types of leakage risks (i.e., wellbore leakage, flow across unfractured caprock) at three modeled sites considered in Task 2.4.

MRCI,Report,Summary,Technical Challenges,dynamic m↗

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗

toughio v1.2.2

TOUGH is a widely used general purpose numerical simulator designed for fluid and heat flows of multiphase, multicomponent fluid mixtures in porous and fractured media, which has been applied to many real-world problems such as underground geological storage, geothermal reservoir engineering or nuclear waste disposal, to name a few. `toughio` is a lightweight, object-oriented and vectorized Python library that aims to provide user-friendly routines to facilitate pre- and post-processing of a TOUGH simulation. Although its primary target is the latest version TOUGH3, `toughio` can also be used to read/write TOUGH2 input/output files. Currently, to the best of our knowledge, only PyTOUGH offers an exhaustive number of features to carry out a complete TOUGH simulation using a scripting language. `toughio` and PyTOUGH share the same objectives, yet with different approaches and philosophies. On the one hand, a PyTOUGH mesh is represented as a MULGRAPH geometry where elements can be unstructured horizontally but only layered vertically. It also checks that the model parameters are valid to ensure that a TOUGH simulation runs smoothly. On the other hand, although it provides basic meshing features, `toughio` mostly relies on common third-party softwares (e.g. Abaqus, FLAC3D, Gmsh, LaGriT) to generate the mesh by importing and converting it to a TOUGH mesh, which also conveniently facilitates the coupling of TOUGH with any other simulator that also supports the same mesh formats. Besides, `toughio` is designed to be just a TOUGH input and output file reader/writer. It assumes that the user knows what he is doing and does not verify the validity of the model parameters (`toughio` does verify the functional arguments).

Luu, Keurfon↗

Fourier-MIONet: Fourier-enhanced multiple-input neural operators for multiphase modeling of geological carbon sequestration

Geologic carbon sequestration (GCS) is a safety-critical technology that aims to reduce the amount of carbon dioxide in the atmosphere, which also places high demands on reliability. Multiphase flow in porous media is essential to understand CO 2 migration and pressure fields in the subsurface associated with GCS. However, numerical simulation for such problems in 4D is computationally challenging and expensive, due to the multiphysics and multiscale nature of the highly nonlinear governing partial differential equations (PDEs). It prevents us from considering multiple subsurface scenarios and conducting real-time optimization. Here, we develop a Fourier-enhanced multiple-input neural operator (Fourier-MIONet) to learn the solution operator of the problem of multiphase flow in porous media. Fourier-MIONet utilizes the recently developed framework of the multiple-input deep neural operators (MIONet) and incorporates the Fourier neural operator (FNO) in the network architecture. Once Fourier-MIONet is trained, it can predict the evolution of saturation and pressure of the multiphase flow under various reservoir conditions, such as permeability and porosity heterogeneity, anisotropy, injection configurations, and multiphase flow properties. Compared to the enhanced FNO (U-FNO), the proposed Fourier-MIONet has 90% fewer unknown parameters, and it can be trained in significantly less time (about 3.5 times faster) with much lower CPU memory (<15%) and GPU memory (<35%) requirements, to achieve similar prediction accuracy. In addition to the lower computational cost, Fourier-MIONet can be trained with only 6 snapshots of time to predict the PDE solutions for 30 years. Furthermore, we observed that Fourier-MIONet can maintain good accuracy when predicting out-of-distribution (OOD) data. The excellent generalizability of Fourier-MIONet is enabled by its adherence to the physical principle that the solution to a PDE is continuous over time. Furthermore, the developed Fourier-MIONet makes it possible to solve the long-time evolution of geological carbon sequestration in a large-scale three-dimensional space accurately and efficiently.

97 MATHEMATICS AND COMPUTING↗

Rheology and stability of nanoparticle-stabilized CO 2 foam under reservoir conditions

This work investigates the rheology and stability of nanoparticle-stabilized CO 2 foams under reservoir conditions (high temperature and high pressure) for fracturing applications. The effects of different parameters on foam apparent viscosity and foam stability were experimentally investigated including the effects of nanoparticle concentration, salinity, foam quality (Γ), shear rate and temperature. The power law model was applied to calculate foam apparent viscosity owing to its pseudo-plastic behavior, and the changes of foam heights over time were used to evaluate foam stability. Results showed that the CO 2 foam apparent viscosity featured a mountain-shaped curve versus Γ, with the peak apparent viscosity obtained at 70% foam quality. The increase of salinity (up to 11%) in the continuous phase improved both foam stability and rheology. Higher nanoparticle concentration could contribute to better foam stability, but there was a threshold concentration, above which the foam apparent viscosity remained stabilized. The CO 2 foam stabilized by nanoparticles displayed a shear-thickening behavior as the experiment flow rate increased from 6 mL/min to 18 mL/min. Further studies showed that as the total flow rate increased, the CO 2 foam became finer-textured with better stability. Elevated temperatures could undermine foam apparent viscosity and long-term stability. The results of this study could provide guidelines as to the design of foam fracturing systems for field applications.

03 NATURAL GAS↗

Techno-economic feasibility of geothermal energy production using inactive oil and gas wells for district heating and cooling systems in Tuttle, Oklahoma

Oil and gas wells have been drilled in the United States and often encounter hot co-produced water possibly suitable for geothermal direct-use applications. Here this study evaluated geothermal resources, heating and cooling demand, and techno-economic potential of four oil and gas wells repurposed for geothermal energy production to serve two public schools and 250 nearby houses in Tuttle, Oklahoma, USA. Subsurface geology in the study area primarily consists of Permian to Mississippian sandstone and limestone formations approximately from 2 km to 3.5 km depth. With a geothermal gradient of 25.3 °C/km, the temperature of geothermal energy production was estimated as 68.2 °C at 2.1 km. The potential of calcite scaling within boreholes and pipes was minimal. Using the characterized reservoir properties, techno-economic analyses were conducted for evaluating levelized costs of geothermal energy production for direct-use heating and cooling and district heating systems with five different production scenarios. Heating and cooling demand in two schools and 250 houses and additional 0.3 MW and 0.6 MW base heating loads for potential geothermal direct-use applications in the study area were also incorporated in the modeling. Results indicated that the levelized cost of heating of geothermal district heating systems utilizing existing wells were significantly lower than those of conventional geothermal energy systems that include well drillings. The geothermal energy production with doublet and quartet configurations was generally sufficient to supply total annual heating demand in the district, while a peaking boiler was used to supply peak loads. Doublet configurations showed higher production temperature with a lower thermal drawdown than the quartet configurations. The doublet system at 2.1 km production depth with 9.3 kg/s flow rate represented the optimal levelized cost of $71/MWh and 91.2% utilization factor. These results imply that the Tuttle geothermal district heating system is techno-economically feasible and comparable to geothermal district heating systems in the United States and the regional natural gas rate for residentials by recycling inactive oil and gas wells with no drillings.

15 GEOTHERMAL ENERGY↗

Filtration of Hanford Tank 241-AN-107 Supernatant at 16 °C

Approximately 9 liters of supernatant from Hanford waste tank 241-AN-107 was delivered by Washington River Protection Solutions to the Radiochemical Processing Laboratory (RPL) at Pacific Northwest National Laboratory. The thirty-six AN-107 sample bottles consisted of six sets of six samples, with each set pulled from a unique tank sampling level. Prior to testing, samples from each level were composited to provide nominally level-independent feed for dead end filtration and ion exchange testing. The composited 241-AN-107 supernatant was chilled to 16 °C for 1 week prior to testing. Filtration testing was then conducted using a backpulse dead-end filter (BDEF) system equipped with a feed vessel and a Mott inline filter Model 6610 (Media Grade 5) in the hot cells of the RPL. The purpose of this testing is to a) demonstrate dead-end filtration (DEF) of AN-107 feed at reduced temperature to obtain prototypic tank side cesium removal (TSCR) flux rates and identify issues that may impact filtration after dilution to 5.5M Na, and b) provide feed for a follow on ion exchange unit operation. The feed was filtered through the BDEF system at a targeted flux of 0.065 gpm/ft 2 . During filtration the differential pressure required to effect filtration at 0.065 gpm/ft 2 was slow to increase for most of the filtration campaign. After all the feed bottles had been pumped into the slurry reservoir, the bottoms of the bottles were added to the reservoir and transmembrane pressure (TMP) reached 2.0 psid (the TSCR action limit). The prototypic filter cleaning process was unable to effectively restore filter performance, and cleaning with oxalic acid was required before flow through the filter could be restored. This indicates that the Media Grade 5 filter may require an alternative cleaning protocol when processing AN-107 supernatant. After completing filtration of the AN-107 feed, the filter was cleaned. Solids concentrated from the backpulse solutions were composed of natrophosphate, Mn-Fe phases, and fluoro-natrophosphate that occurred as particle agglomerates. The individual particles were in some cases 100s of micrometers across which is consistent with prior observations from AN-107 supernate waste characterizations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of hydrogen fueling, recycling, and confinement at Wendelstein 7-X via a single-reservoir particle balance

Abstract A single-reservoir particle balance for the main plasma species hydrogen has been established for W7-X. This has enabled the quantitative characterization of the particle sources in the standard island divertor configuration for the first time. Findings from attached scenarios with two different island sizes with a boronized wall and turbo molecular pumping are presented. Fueling efficiencies, particle flows and source locations were measured and used to infer the total particle confinement time τ p . Perturbative gas injection experiments served to measure the effective particle confinement time τ p * . Combining both confinement times provides access to the global recycling coefficient R ¯ . Hydrogen particle inventories have been addressed and the knowledge of particle sources and sinks reveals the core fueling distribution and provides insight into the capability of the magnetic islands to control exhaust features. Measurements of hydrogen fueling efficiencies were sensitive to the precise fueling location and measured between 12% and 31% with the recycling fueling at the strike line modeled at only 6%, due to much higher densities. 15% of the total 5.2 × 10 22 a/s recycling flow ionizes far away from the recycling surfaces in the main chamber. It was shown that 60% of recycled particles ionize above the horizontal and 18% above the vertical divertor target, while the remainder of the recycling flow ionizes above the baffle (7%). Combining these source terms with their individual fueling efficiencies resolves the core fueling distribution. Due to the higher fueling efficiency in the main chamber, up to 51% of the total 5.1 × 10 21 s −1 core fueling particles are entering the confined plasma from the main chamber. τ p values in the range of 260 ms were extracted for these discharges. Together with τ p , the global recycling coefficient R ¯ was resolved for every τ p * measurement and a typical value close to unity was obtained. An increase of the island size, resulted in no change of τ p , but doubled τ p * , indicating the feasibility of the control coils as an actuator to control exhaust features without affecting core confinement properties.

recycling↗

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

Utilizing Temperature and Brine Inflow Measurements to Constrain Reservoir Parameters During a Salt Heater Test

Brine availability in salt has multiple implications for the safety and design of a nuclear waste storage facility. Brine availability includes both the distribution and transport of brine through a damaged zone around boreholes or drifts excavated into the salt. Coupled thermal, hydrological, mechanical, and chemical processes taking place within heated bedded salt are complex; as part of DECOVALEX 2023 Task E this study takes a parsimonious modeling approach utilizing analytical and numerical one-dimensional simulations to match field measurements of temperature and brine inflow around a heater. The one-dimensional modeling results presented arrive at best-fit thermal conductivity of intact salt, and the permeability and porosity of damaged salt of 5.74 W/m∙K, 10 –17 m 2 , and ≈ 0.02, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Shifts in hydropower operation to balance wind and solar will modify effects on aquatic biota

To avoid negative consequences to freshwater biota from climate change, society must complete the transition from fossil to renewable electricity sources. However, temporal patterns in hydropower generation (and flow releases that affect aquatic biota) may change with increased wind and solar penetration. We used power cost modeling to characterize current and future within-day and seasonal patterns in hydropower generation across the Eastern Interconnection in a wet and a dry year. Compared to the baseline, future hydropower generation across the grid decreased during the day and increased before dawn and after dusk. At a project level, such a pattern would suggest 'double peaking' operation (up- and down-ramping before dawn and after dusk, with lower releases midday). Variation in generation was higher in wet years than dry years, foreshadowing possible flow constraints on hydropower flexibility. At the grid scale, projected ramping rates were higher in all seasons. A review of the ecological literature suggests that these changes would shift the timing of invertebrate drift and elevate the risk of nest scouring during up-ramping and the risk of stranding or dewatering during down ramping. Thermal conditions may be moderated by increased ramping. Strategies for adapting to future shifts in the renewable portfolio range from re-regulation in reservoir cascades to providing flow refuge (structures and vegetation) below individual projects. Coordinated basin-scale operation can distribute peaking operation to maintain grid support while restricting local ramping at critical ecological times. In addition, research to design hybrid renewable systems that add battery storage is needed to understand how we can mitigate future risks to aquatic communities while promoting the use of renewable energy. This study, which is among the first to examine ecological side-effects of the shift to renewable energy in freshwater ecosystems, lays out a path toward understanding and navigating changes to flow regimes under the energy transition.

54 ENVIRONMENTAL SCIENCES↗

Floating solar panels on reservoirs impact phytoplankton populations: A modelling experiment

Floating solar photovoltaic (FPV) deployments are increasing globally as the switch to renewable energy intensifies, representing a considerable water surface transformation. FPV installations can potentially impact aquatic ecosystem function, either positively or negatively. However, these impacts are poorly resolved given the challenges of collecting empirical data for field or modelling experiments. In particular, there is limited evidence on the response of phytoplankton to changes in water body thermal dynamics and light climate with FPV. Given the importance of understanding phytoplankton biomass and species composition for managing ecosystem services, we use an uncertainty estimation approach to simulate the effect of FPV coverage and array siting location on a UK reservoir. FPV coverage was modified in 10% increments from a baseline with 0% coverage to 100% coverage for three different FPV array siting locations based on reservoir circulation patterns. Results showed that FPV coverage significantly impacted thermal properties, resulting in highly variable impacts on phytoplankton biomass and species composition. The impacts on phytoplankton were often dependent on array siting location as well as surface coverage. Changes to phytoplankton species composition were offset by the decrease in phytoplankton biomass associated with increasing FPV coverage. We identified that similar phytoplankton biomass reductions could be achieved with less FPV coverage by deploying the FPV array on the water body's faster-flowing area than the central or slower flowing areas. The difference in response dependent on siting location could be used to tailor phytoplankton management in water bodies. Simulation of water body-FPV interactions efficiently using an uncertainty approach is an essential tool to rapidly develop understanding and ultimately inform FPV developers and water body managers looking to minimise negative impacts and maximise co-benefits.

14 SOLAR ENERGY↗

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗