Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

The Energy & Environmental Research Center (EERC), in partnership with the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the EERC’s State Energy Research Center (SERC), MPLX Operations LLC, and TC Energy Development Holdings Inc. (a subsidiary of TC Energy Corporation), studied the potential for subsurface hydrogen storage and recovery in the Williston Basin of western North Dakota. The project’s goal was to evaluate the feasibility of large-scale, secure geologic H 2 storage to support future hydrogen generation, storage, and use. This work included laboratory testing, H 2 –rock–fluid exposure experiments, literature reviews on H 2 embrittlement, and reservoir modeling and simulations. The study included an assessment of storage potential across three types of storage reservoirs using both reservoir simulation and DOE’s web-based tool SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration – Hydrogen Estimator for Logistical Planning), as well as investigation of potential H 2 production and markets for commercial-scale deployment. Building on prior EERC gas storage research, three storage options were selected for detailed evaluation: the Broom Creek Formation (a clastic saline reservoir), the Dickinson Lodgepole Mounds (DLM) complex (carbonate mud mound structures) of the Lodgepole Formation (an active oil and gas producing reservoir), and the Dunham Salt Interval of the Piper Formation (to be used for engineered salt cavern development). These targets were prioritized based on prior EERC research using datasets related to seal capacity, reservoir quality, mechanical integrity, and injectivity. Exposure tests on Broom Creek and DLM samples showed mineral dissolution and precipitation that increased brine salinity and altered reservoir rock surfaces. Although these results provide useful insight, they are limited by small sample sizes and short-term (30-day) exposure, requiring further study to assess long-term storage integrity. Salt formations were not tested because of their known nonreactivity and established mechanical stability. Results of reservoir simulations performed for a single site demonstrated that the Broom Creek Formation may be capable of receiving up to 42,000 tonnes of injected H 2 over 7 months via one well. H 2 recovery took place over 5 months, resulting in approximately 26,000 tonnes (~62% without cushion gas [CG]). This work suggests water production may be important and subsequent cycles of injection and production may perform more efficiently; however, significant site-specific work in the future is needed to assess actual reservoir performance of injection and withdrawal of H 2 storage. For oil reservoir potential, a multiple-well model was used to simulate injection of approximately 32,000 tonnes of H 2 into a single wellbore while simultaneously producing in place reservoir fluids from four offset wells to maintain reservoir pressure. The simulation results suggested a high recovery (~98%); in addition, cost advantages through existing infrastructure could be realized. Challenges in this reservoir include vi managing gas purity and leakage risks. In both scenarios, production of H 2 takes place in a single-well scenario with 10 cycles (7 months of injection and 5 months of production) over 10 years. Finally, the use of engineered caverns in the Dunhan Salt was evaluated, and the results suggest that while they have a smaller capacity (<1000 tonnes per cavern), they exhibit nearly complete gas recovery (>99%), fast response times, and low purity risk. While caverns in North Dakota may be smaller in capacity, fields can be developed in galleries to accommodate the volumetric needs and rapid turnaround times necessary to meet market demands. Geographic limitations and thin salt intervals in North Dakota may represent less total storage potential than salt domes elsewhere, but significant opportunities exist to expand this market for gas storage in North Dakota. A basinwide assessment was performed to estimate a first-of-its-kind value for H 2 storage on a large scale. DOE’s SHASTA-HELP, combined with EERC simulation work, was used to perform the assessment. Estimated H 2 storage potential varied widely for each formation type. The Broom Creek saline formation was estimated to have a storage potential of approximately 1.7–90.5 million tonnes (MMt). The DLM oil reservoirs were estimated to have 0.07–0.19 MMt of capacity. Notably, each of these estimates relies on significant assumptions regarding reservoir thickness, porosity, permeability, and CG needed for operation. Much research is needed to understand the true site-specific storage resource potential of each formation. Using the Dunham Salt Interval for cavern development may result in as much as 4.79 MMt (up to 2.87 MMt working gas) of H 2 storage potential. An important note for consideration is that the values presented here need significantly more geological characterization and engineering assessments prior to gaining confidence in performance. This will be a focal point for future research and development needs. The basinwide evaluation also indicated that North Dakota has significant H 2 generation potential, with estimates up to a possible 13 MMt annually, suggesting a substantial opportunity for H 2 market development and thus the need for commercial-scale H 2 storage to facilitate growth and resilience. Key Recommendations 1. Conduct detailed site characterization (3D seismic, well logs, core sampling) to reduce geologic uncertainty. 2. Perform techno-economic analyses incorporating market, regulatory, and incentive frameworks. 3. Investigate long-term interactions among H 2 , CGs, well materials, and formations to assess risks. 4. Develop pilot- and field-scale demonstrations to validate models and establish best practices.

03 NATURAL GAS↗

Development Fiber Optic Distributed System for Direct Detection of Subsurface Gases Leakages

Carbon, natural gas, and hydrogen gas storage is an emerging solution to safeguard us against pollution, support goals of negative carbon emission, and protect sources of renewable energy. Properly constructed storage wells provide a virtually impervious barrier to any unintended subsurface transmission. The ability to ensure the long-term integrity of such wells is vital to the success of any storage operation and be successful in the public eyes. Therefore, robust monitoring of any gas migration into the subsurface is highly sought. A fiber-optic distributed chemical sensor (DCS) enables monitoring of long-term well integrity along its depth, ensuring the success of any storage operation and bolsters public acceptance of the safety of the reservoir via leak early detection. The same technique can be applied to gas monitoring in pipeline networks and nuclear stockpile monitoring applications. Fiber based Raman spectroscopy enables DCS, as optical fibers can be deployed in virtually any environment and relay spectroscopic information over long distances back to the user. Hollow core fibers (HCF) make excellent DCSs as the air core of the fiber allows gas from the environment to diffuse into the core, which interacts with the laser signal that is carried in the air core. This work builds upon the previous LDRD project, Fiber Optic System for Direct Detection of Carbon Dioxide Leakage in Carbon Storage Wells (21-FS-003), in which the feasibility of Raman spectroscopy detection of Carbon Dioxide (CO2) in HCF detection was demonstrated. We mitigated the risk of this DCS technology by establishing and completing five objectives. The first objective was to model and optically characterize HCF uptake of CO2, establishing the relationship between HCF length, gas diffusion time, detectable gas concentration, and measured Raman intensity. In objective two, we developed a fiber core drilling recipe to enable additional diffusion ports in the fiber core and established a method for maintaining fiber strength and integrity post drilling. Objective three characterized the drilled fibers against the undrilled fibers, establishing the differences in the gas mechanics and optical properties and provided parameters to iterate the drilling process. In objective four, a fusion splicing technique was developed to join the HCF to conventional single-mode fibers, localizing the gas detection point at the drilled HCF hole, emulating a DCS. Lastly, objective five was the testing of the sensor in Edgar Mines at Colorado School of Mines on a CO2 pipeline with a simulated leak, to showcase the ability to detect CO2 leaks. This capstone result showed CO2 leak detection in < 10 minutes, raising the technology readiness level of HCF segments as deployable DCS.

organic↗

Self-Healing Cements with Improved Toughness at Casing and Formation Interfaces for Subsurface Applications (CRADA 530)

Over the three-year CRADA, Pacific Northwest National Laboratory (PNNL) advanced a new class of Molecular Velcro™ self-healing wellbore cements designed to improve the long-term integrity of wells used for geothermal energy production and carbon storage. Conventional wellbore cement can crack due to repeated mechanical loading, temperature swings, and exposure to aggressive subsurface fluids. Once cracks form, they can create leakage pathways, drive expensive remediation (“workovers”), and in extreme cases shorten well life. This project addressed that problem by developing polymer-modified cement systems that can autonomously re-seal cracks and better tolerate the extreme conditions experienced in subsurface wells, while remaining compatible with practical cementing operations.

15 GEOTHERMAL ENERGY↗

Subsurface hydrological controls on the short-term effects of hurricanes on nitrate–nitrogen runoff loading: a case study of Hurricane Ida using the Energy Exascale Earth System Model (E3SM) Land Model (v2.1)

When the nutrient level in the soil surpasses vegetation demand, nutrient losses due to surface runoff and subsurface leaching are the major reasons for the deterioration of water quality. The lower Mississippi River basin (LMRB) is one of the sub-basins that deliver the highest nitrogen loads to the Gulf of Mexico. Potential changes in episodic events induced by hurricanes may exacerbate water quality issue in the future. However, uncertainties in modeling the hydrologic response to hurricanes may limit the modeling of nutrient losses during such events. Using a machine learning approach, we calibrated the land component of the Energy Exascale Earth System Model (E3SM), or ELM, version 2.1, based on the water table depth (WTD) of a calibrated 3D subsurface hydrology model. While the overall performance of the calibrated ELM is satisfactory, some discrepancies in WTD remain in slope areas with low precipitation due to the missing lateral flow process in ELM. Simulations including biogeochemistry performed using ELM with and without model calibration showed important influences of soil hydrology, precipitation intensity, and runoff parameterization on the magnitude of nitrogen runoff loss and the leaching pathway. Despite such sensitivities, both ELM simulations produced reduced WTD and increased runoff and accelerated nitrate–nitrogen runoff loading during Hurricane Ida in August 2021, consistent with the observations. With observations suggesting more pronounced effects of Hurricane Ida on nitrogen runoff than the simulations, we identified factors for model improvement to provide a useful tool for studying hurricane-induced nutrient losses in the LMRB region.

54 ENVIRONMENTAL SCIENCES↗

Machine learning deciphers CO 2 sequestration and subsurface flowpaths from stream chemistry

Endmember mixing analysis (EMMA) is often used by hydrogeochemists to interpret the sources of stream solutes, but variations in stream concentrations and discharges remain difficult to explain. We discovered that machine learning can be used to highlight patterns in stream chemistry that reveal information about sources of solutes and subsurface groundwater flowpaths. The investigation has implications, in turn, for the balance of CO 2 in the atmosphere. For example, CO 2 -driven weathering of silicate minerals removes carbon from the atmosphere over ~106-year timescales. Weathering of another common mineral, pyrite, releases sulfuric acid that in turn causes dissolution of carbonates. In that process, however, CO 2 is released instead of sequestered from the atmosphere. Thus, understanding long-term global CO 2 sequestration by weathering requires quantification of CO 2 - versus H 2 SO 4 -driven reactions. Most researchers estimate such weathering fluxes from stream chemistry, but interpreting the reactant minerals and acids dissolved in streams has been fraught with difficulty. We apply a machine-learning technique to EMMA in three watersheds to determine the extent of mineral dissolution by each acid, without pre-defining the endmembers. The results show that the watersheds continuously or intermittently sequester CO 2 , but the extent of CO 2 drawdown is diminished in areas heavily affected by acid rain. Prior to applying the new algorithm, CO 2 drawdown was overestimated. The new technique, which elucidates the importance of different subsurface flowpaths and long-timescale changes in the watersheds, should have utility as a new EMMA for investigating water resources worldwide.

54 ENVIRONMENTAL SCIENCES↗

State of Stress Uncertainty Quantification and Geomechanical Risk Analysis for Subsurface Engineering

All subsurface stress measurement methods available for deep subsurface applications are indirect and involve a significant amount of uncertainty. Quantifying these uncertainties is critical for applications where the state of stress is a component of a risk analysis with significant safety, environmental, and economic considerations. A publicly available and open source tool called the State of Stress Analysis Tool (SOSAT) has been developed to aid in performing a Bayesian uncertainty quantification and geomechanical risk analysis using many types of commonly available data. A new feature of this tool is using data from the FutureGen 2.0 site. The new feature of the tool was uses information about the absence or presence of borehole breakouts to constrain the state of stress. This new feature allows the user to specify probability distributions for the relevant parameters such as the drilling mud pressure, temperature, and the formation strength properties. The tool computes the posterior joint probability distribution of the horizontal principal stresses. Using random samples from this posterior distribution the tool then computes the probability of activating a critically oriented fault at a range of pore pressures. This latter feature can be used to aid in the assessment of the risk of induced seismicity and leakage in geologic carbon sequestration and enhanced geothermal systems.

Burghardt, Jeffrey A.↗

Methods and systems for obtaining reconstructed low-frequency seismic data for determining a subsurface feature

A computer-implemented method for obtaining reconstructed seismic data for determining a subsurface feature, includes: determining an initial training velocity model, training a machine learning model based on first training seismic data and second training seismic data generated from the training velocity model, the first training seismic data corresponding to one or more first frequencies, the second training seismic data corresponding to one or more second frequencies lower than the one or more first frequencies, obtaining, based on measured seismic data and the machine learning model, reconstructed seismic data corresponding to the one or more second frequencies, generating a velocity model based on the measured seismic data, the reconstructed seismic data, and a full waveform inversion (FWI), and when the generated velocity model does not satisfy a preset condition, updating the training velocity model based on the generated velocity model, to obtain updated reconstructed seismic data for determining a subsurface feature.

58 GEOSCIENCES↗

Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects

This is the conference paper accompanying an oral presentation “Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects” at the 17th International Conference on Greenhouse Gas Control Technologies GHGT-17 held in Calgary, Canada, October 20-24 , 2024. Carbon capture and storage (CCS) technology is critical for mitigating climate change but requires effective subsurface reservoir management to ensure safe containment of injected CO2. Accurate predictions of reservoir pressure and saturation are essential for assessing long-term CCS performance. Traditional numerical simulations, while effective, are computationally intensive, time-consuming, and constrained by data discretization. Previous work has shown the effectiveness of MeshGraphNets (MGN), a graph-based machine learning framework, as an innovative alternative for predicting reservoir behavior. MGN leverages graph neural networks (GNNs) and mesh representations to model complex geological formations, offering superior adaptability across different discretizations and reservoir configurations. Classic MGN implementations utilize an autoregressive technique to predict future behavior based on current predictions, but this technique is hampered by error accumulation over time. To enhance the model accuracy in time-series predictions, this study implemented a multi-step rollout strategy that integrates autoregressive predictions during training to stabilize prediction of saturation over time. Using the Illinois Basin – Decatur Project (IBDP) dataset, comprising 100 simulations of CO2 injection, pressure, and saturation changes, the framework demonstrated its ability to learn spatial dependencies and temporal dynamics. With inputs including permeabilities, porosities, and injection rates, MGN accurately predicted CO2 plume evolution over time, even with limited training data. Moreover, the addition of a multi-step rollout procedure during training improved the ability of MGN to predict stably over time by ~15%. This research positions MGN, enhanced with multi-step rollout capabilities, as a robust and efficient tool for CCS applications. It advances the field by enabling precise, computationally efficient predictions of reservoir behavior, providing a foundation for the broader adoption of machine learning frameworks in CCS and other geoscience domains.

Holcomb, Paul↗

Editorial: Subsurface microbiology within hydrocarbon resources or stored gases

A Research Topic on the microbiology of hydrocarbon and gas storage reservoirs has far reaching industrial applications. In recent decades, there has been a growing interest in understanding microbial communities in subsurface energy reservoirs, such as coal, oil, and shale beds. This area of research has broadened to include gas storage reservoirs for hydrogen and CO 2 . Scientists are beginning to unravel the unexpected impact microorganisms have on these systems, through changing the fluid geochemistry, the gas content, and even the permeability. By recognizing the influence of these tiny organisms on our engineered environments, we can develop better risk assessments, target mitigation strategies, expand energy production, and refine operational guidance, ultimately contributing to a more sustainable energy future.

59 BASIC BIOLOGICAL SCIENCES↗

Subsurface mafic and ultramafic rock mapping and analysis for carbon mineralization in the United States (SubMAP-CO 2 )

Achieving durable, gigaton-scale carbon dioxide (CO₂) removal will require the development of storage strategies that are both permanent and scalable. Among the available options, carbon mineralization in mafic and ultramafic rocks offers a uniquely robust pathway for long-term CO₂ sequestration, converting CO₂ into thermodynamically stable carbonate minerals. While surface exposures of these lithologies have been extensively studied, particularly in ophiolitic and volcanic settings, the potential for carbon mineralization in the subsurface remains comparatively underexplored at national scale. This gap limits the ability to assess storage capacity, accessibility, and deployment pathways for mineralization-based carbon management in the United States.

08 HYDROGEN↗

Degradation Diagnostics from the Subsurface of Lithium-Ion Battery Electrodes

Despite the long-established rocking-chair theory of lithium-ion batteries (LIBs), developing novel characterization methodology with higher spatiotemporal resolution facilitates a better understanding of the solid electrolyte interphase studies to shape the reaction mechanisms. In this work, we develop a Xenon ion plasma focused ion beam (Xe + PFIB)-based characterization technique to probe the cross-sectional interface of both ternary cathode and graphite anode electrodes, with the focus on revealing the chemical composition and distribution underneath the electrode surface by in-depth analysis of secondary ions. Particularly, the lithium fluoride is detected in the pristine cathode prior to contact with the electrolyte, reflecting that the electrode degradation is in the form of the loss of lithium inventory during electrode preparation. This degradation is related to the hydrolysis of the cathode material and the decomposition of the PVDF binder. Here, through the quantitative analysis of the transition-metal degradation products, manganese is found to be the dominant element in the newly formed inactive fluoride deposition on the cathode, while no transition metal signal can be found inside the anode electrode. These insights at high resolution implemented via a PFIB-based characterization technique not only enrich the understanding of the degradation mechanism in the LIBs but also identify and enable a high-sensitivity methodology to obtain the chemical survey at the subsurface, which will help remove the capacity-fade observed in most LIBs.

25 ENERGY STORAGE↗

The subsurface structure of abraded Al–Zn–Mg–Cu alloy

Surface abrasion has been shown to introduce heavy shear stress and produce a few hundred nanometers thick altered surface layer (ASL) on Al alloys. Such ASLs exhibit microstructural changes at room temperature and different corrosion resistance than the underlying substrate. Here, we report how sample dimensions and temperature affect the deformed subsurface microstructure evolution on an abraded Al–Zn–Mg–Cu alloy using transmission electron microscopy (TEM) and atom probe tomography. The ASL microstructure evolves differently in the bulk sample than in the thin TEM specimen. In the freshly abraded ASL, the original η'/η precipitates and grain boundaries were deformed toward the abrasion direction. After 7-day natural aging of the freshly abraded bulk sample, well-defined subgrains formed locally in the ASL, and unusual Al 2 Cu (θ) and Mg x Zn y phases precipitated at the subgrain boundaries. However, during natural aging of this TEM specimen taken from the abraded and aged bulk samples, much faster changes occurred, including the continued formation of subgrains and θ particles in the entire ASL, and the diffusion of Zn atoms out of Mg x Zn y particles to form a Zn particle. On the other hand, no obvious microstructural changes occurred in the ASL during 3-month natural aging of the TEM specimen immediately taken from the freshly abraded sample. Finally, this work advances the understanding of how Al alloy microstructures respond during and following any shear deformation, and complex and unexpected effects of sample dimensions on the behavior.

36 MATERIALS SCIENCE↗

Accessing the Subsurface Biosphere Within Rocks Undergoing Active Low-Temperature Serpentinization in the Samail Ophiolite (Oman Drilling Project)

The Oman Drilling Project established an “Active Alteration” multi-borehole observatory in peridotites undergoing low-temperature serpentinization in the Samail Ophiolite. The highly serpentinized rocks are in contact with strongly reducing fluids. Distinct hydrological regimes, governed by differences in rock porosity and fracture density, give rise to steep redox (Eh +200 to –750 mV) and pH (pH range 8.5–11.2) gradients within the 300–400 m deep boreholes. The serpentinites and fluids host an active subsurface ecosystem. Microbial cell abundances in serpentinite vary at least six orders of magnitude, from ≤3.5 × 10 1 to 2.9 × 10 7 cells/g. Low levels of biological sulfate reduction (2–1,000 fmol/cm 3 /day) can be detected in rock cores, particularly in rocks in contact with reduced groundwaters with pH < 10.5. Thermodesulfovibrio is the predominant sulfate reducer identified via metagenomic sequencing of adjacent groundwater communities. Here, we infer that transport and reaction of microbially generated sulfide with the serpentine and brucite assemblages gives rise to optical darkening and sulfide overprinting, including the formation of tochilinite-vallerite group minerals, potentially serving as an indicator that this system is inhabited by microbial life. Olivine mesh-cores replaced with ferroan brucite and minor awaruite, abundant veins containing hydroandradite garnet and polyhedral serpentine, and late-stage carbonate veins are suggested as targets for future spatially resolved life-detection investigations. The high-quality whole-round core samples that have been preserved can be further probed to define how life distributes itself and functions within a system where chemical disequilibria are sustained by low-temperature water/rock interaction, and how biosignatures of in situ microbial activity are generated.

58 GEOSCIENCES↗

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗

Model data for Flood Frequency Analysis using Stochastic Storm Transposition and an Integrated Surface-Subsurface Hydrological Model

This archived provides scripts and input files used for the implementation of a novel approach to conduct process-based Flood Frequency Analysis using a Stochastic Storm Transposition (SST) and an Integrated Surface-Subsurface Hydrological Model (ISSHM). As a proof-of-concept, this study uses the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct flood frequency analysis by simulating the flood response to 5,000 annual synthetic storm events in a ~2000 km2 Southeast Texas watershed.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include TXT, CSV, DAT, SBATCH, SHP, TIF, NetCDF, and HDF5 files, which can be read through Python scripts. The input files for the ATS model and RainyDay model have .XML and .SST extensions, respectively, and can be edited in any commonly used text editors.This archive contains:* Scripts and data files essential for generating the ATS model input. It uses the Watershed Workflow package to produce both mesh and ATS input files. * Jupyter notebooks designated for the ATS model evaluation, covering both long-term simulations and 40 rainfall-runoff events.* Input files required to simulate SST storm events using RainyDay.

54 ENVIRONMENTAL SCIENCES↗

Fractional Modeling in Action: A Survey of Nonlocal Models for Subsurface Transport, Turbulent Flows, and Anomalous Materials

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas – subsurface transport, turbulence, and anomalous materials – in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

97 MATHEMATICS AND COMPUTING↗