Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Leakage From Coexisting Geologic Forcing and Injection‐Induced Pressurization: A Semi‐Analytical Solution for Multilayered Aquifers With Multiple Wells

Abstract Abnormal fluid pressures (above or below hydrostatic pressure) can develop and persist in sedimentary basins. The common occurrence of abnormal pressures may cause challenges for project permitting of geological carbon sequestration (GCS), particularly in reservoirs with pre‐injection overpressure. The leaky wells that may exist in some sedimentary basins can provide flow paths between deep brine aquifers and shallower freshwater aquifers. Pre‐injection relative overpressures can cause brine leakage through leaky wells even before any injection occurs. The tendency for flow through leaky wells is coupled with the process of pressure dissipation that occurs through aquitards. Specifically, with non‐zero permeability, aquitards can dissipate pressure over large areal extents and thereby reduce leakage rates through leaky wells. This study presents development of a semi‐analytical solution for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. The geologic forcing that causes abnormal pressures in the multilayered system can coexist with any number of injection, extraction, and leaky wells that also affect fluid pressure. The semi‐analytical model is applied to explore how leakage through leaky wells varies as functions of pressurization rate, along with aquitard and leaky well properties in an overpressured multilayered system. The results show that although injection‐induced pressures can dissipate rapidly through suitably permeable aquitards, coexisting geologic forcing may create sustained rates of brine leakage into freshwater aquifers through leaky wells. In GCS, a very low‐permeability aquitard with high capillary entry pressure to free‐phase CO 2 is desired to serve as the caprock to prevent leakage of CO 2 from the storage reservoir. Nevertheless, the results from this study show that the brine leakage impact to shallow freshwater aquifers through leaky wells might substantially decrease with increasing aquitard permeability values, as long as small aquitard permeability and high capillary entry pressure serve to prevent CO 2 leakage.

Cihan, A.↗
DOE OSTI2020

Influence of Clay Wettability Alteration on Relative Permeability

Abstract Understanding the wettability of porous materials is important to model fluid flow in the subsurface. One of the critical factors that influences wetting in real reservoirs is the composition of geologic materials. The wetting properties for clay minerals can have a particularly strong impact on flow and transport. In this work, we analyze the chemical composition of a Mt. Simon sandstone core to resolve the microscopic structure of clay regions and assess how alterations to the local wetting properties influence multiphase transport based on core flooding experiments and relative permeability simulations. We show that whichever fluid has greater affinity toward clay minerals will tend to accumulate within these high surface area regions, leading to dramatic shifts in the relative permeability. This work establishes the essential importance of the mineral composition and associated wetting properties in the modeling of flow and transport in reservoir‐scale systems.

Fan, Ming↗
DOE OSTI2021

Validation of Fault Displacements from Dynamic Rupture Simulations against the Observations from the 1992 Landers Earthquake

ABSTRACT Coseismic fault displacements in large earthquakes have caused significant damage to structures and lifelines on and near fault lines. Coseismic displacements represent a real threat, especially to distributed infrastructure systems. For infrastructure systems that can not avoid active faults, engineering displacement demands are defined using probabilistic fault-displacement hazard analyses (PFDHA). However, PFDHA models are sparse and poorly constrained partly due to the scarcity of detailed fault-displacement observations. Advancements in dynamic rupture simulation methods make them an attractive approach to address this important issue. Because fault displacements can be simulated for various geologic conditions as constrained by current knowledge about earthquake processes, they can be used to supplement the observation datasets. In addition to providing on-fault displacements, when used with appropriate constitutive models for the bulk medium, they can capture off-fault distributed inelastic deformations as well. For viable extrapolation, simulations must first be validated against data. In this article, we summarize the calibration and validation of the dynamic rupture model against the observations of the well-documented 1992 Landers earthquake. We defined a preferred model that reproduces several first-order fault-displacement metrics such as the on-fault partition of the total displacement, the mean fault-zone width, and the location of the peak displacement. Simulated ground motions consistent with the observations ensure that all physics important to modeling have been properly parameterized. For the extrapolation, we generated a suite of dynamic rupture models to quantify expected fault-displacement metrics, their intercorrelations, and magnitude dependencies, which are in part supported by the Landers and other recent earthquakes. Our validation and extrapolation exercise paves the way for using dynamic rupture modeling to quantitatively address fault-displacement hazard on a broader scale. The results are promising and are expected to be useful to inform PFDHA model development.

Geochemistry & Geophysics↗
DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2021

Improving ChemCam LIBS long-distance elemental compositions using empirical abundance trends

The ChemCam instrument on the Curiosity rover provides chemical compositions of Martian rocks and soils using remote laser-induced breakdown spectroscopy (LIBS). The elemental calibration is stable as a function of distance for Ti, Fe, Mg, and Ca. The calibration shows small, systematically increasing abundance trends as a function of distance for Al, Na, K, and to some extent, Si. The distance effect is known to be due to a dependence with distance on the relative strengths of atomic transition lines. Emission lines representing transitions from relatively low energy levels remain intense at longer distances while emission lines representing transitions from higher energy levels decrease in intensity more rapidly as a function of distance. The multivariate algorithms used to determine elemental compositions rely on a large number of emission lines in many cases, so rather than trying to correct all emission lines, a study was made of the predicted compositions as a function of distance, in order to determine an empirical correction. Abundance trends can be well approximated by a linear trend with distance within the ranges of abundances and distances observed up to ~6 m. Data from 11 distinct geological members and data groups of the Murray formation in Gale crater, Mars, were used to form the model, selecting the members and data groups yielding the best statistics. The model was tested using data from several targets observed from two different distances, and using data from the Kimberley formation, the composition of which is significantly different from the Murray formation, showing that the model works on other compositions beyond those used to build the model. For long-distance observations up to ~6 m, corrections can be made back to an equivalent composition at the median distance of ChemCam observations (2.6 m). Finally, the model has been validated up to 6.2 m, although ChemCam is able to observe bedrock targets to >7 m, and iron meteorites to distances of >9 m.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

DEEPEN Data Catalog for Magmatic Geothermal Systems in the United States

This data catalog contains information related to the Training Site Analysis for the Geothermica project "DE-risking Exploration of geothermal Plays in magmatic ENvironments (DEEPEN)." The DEEPEN project aims to reduce exploration risk for geothermal fluids in magmatic systems by developing improved an improved framework for interpretation of exploration data using the Play Fairway Analysis (PFA) methodology. The Training Site Analysis performed for DEEPEN leverages existing datasets to develop a customized PFA approach to exploration for multiple geothermal resource types in magmatic systems (conventional hydrothermal resources, supercritical fluid and superheated steam resources, and superhot EGS resources). This data catalog contains links to publicly available data files related to 8 training sites in the United States. US training sites are: the Cascades/Aleutians PFA project; the Hawaii PFA project, the Oregon Cascades PFA project, the Snake River Plain, Idaho PFA project, the Washington State PFA project, Newberry Volcano, Coso Geothermal Field, and the Geysers Geothermal field. This database contains an overview of these training sites, data sources, and links to publicly available exploration datasets. For the five PFA projects, details on exploration data related to PFA components (heat, fluid, permeability, sometimes seal) are provided, including a summary of data weighting methodologies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Evaluation of Containment and Geomechanical Risks at Integrated Mid-Content Stacked Carbon Storage Hub Sites

The U.S. Department of Energy is leading the development of CCS through collaborative projects and programs to not only develop and demonstrate the technology but to also quantify the risks and provide tools that site operators can use to develop a site-specific understanding of their risk to enable successful project operations. This report summarizes results and findings from the application of risk-assessment tools to an industrial-scale Geologic Carbon Storage (GCS) project. Specifically, we applied the National Risk Assessment Partnership’s (NRAP) Open Integrated Assessment Model (NRAP-Open-IAM) and State of Stress Assessment Tool (SOSAT) to two candidate sites being considered for storage by the Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-Hub) project team. These sites, Sleepy Hollow Field in Nebraska and Patterson Field in Kansas, have historical oil and gas production and thus are attractive candidates to store the 50 million metric tons (Mt), which is the CarbonSAFE objective. Because these sites have historical operations, they have a significant number of existing wells that pose a risk for well leakage. Additionally, the storage formations will undergo significant pore pressure perturbation (i.e., increase due to CO 2 injection). Hence, we have selected the two NRAP tools best suited to study the risks associated with well leakage and geomechanical risks. The objective of the study was not only to assess the risk at the site, but to also improve the NRAP tools through application using real site data on an ongoing project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Insert Modeling in UNF ST&DARDS

The Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS) is a software tool that integrates a used nuclear fuel (UNF) or spent nuclear fuel (SNF) relational database and key analysis capabilities to simplify and automate numerous UNF management and fuel cycle–related activities. UNF-ST&DARDS is being developed for the US Department of Energy’s Office of Nuclear Energy Spent Fuel and Waste Disposition program. UNF-ST&DARDS provides an integrated framework that uses advanced modeling and simulation to predict the behavior of SNF over the timescales associated with permanent disposal in a geologic repository. After leaving the spent fuel pool, SNF is transferred to dry storage in a dual-purpose canister (DPC). DPCs are considered “dual purpose” because they are designed for both storage and transportation, removing the need to transfer the fuel to a separate transportation cask. However, much research has been conducted investigating the feasibility of directly disposing of DPCs in geologic repositories. Direct disposal of DPCs could reduce worker exposure during repackaging, reducing the amount of low-level waste from the discarded DPCs and potentially saving billions of dollars. Therefore, direct disposal of as-loaded DPCs is desirable if it can be done safely.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Geologic Seawater Air Conditioning (GeoSWAC) System: Resource Assessment and Techno-Economic Evaluation in Puerto Rico

Puerto Rico's hot, humid climate drives a high and persistent demand for cooling, straining an aging and fuel-dependent energy grid while increasing peak electricity loads. Much of this challenge stems from inefficient air-conditioning systems operating in buildings without passive cooling design, making cooling both costly and vulnerable to disruption - especially during hurricanes. To address these issues, researchers evaluated a new alternative called Geologic Seawater Air Conditioning (GeoSWAC), which uses inland wells connected to naturally cold, deep seawater. By avoiding long offshore pipelines and energy-intensive refrigeration cycles, GeoSWAC can cut electricity use by 75-90%, reduce environmental impacts, and operate more reliably during power outages. A case study at the University of Puerto Rico's Rio Piedras campus found that GeoSWAC could deliver significantly lower levelized costs of cooling, reduced operational expenses, and improved water conservation compared to the campus's existing chilled-water plant. The system offers key advantages - including access to a constant cold heat sink, low pumping requirements, and enhanced resilience - for coastal regions such as Puerto Rico, the broader Caribbean, and Florida. However, its effectiveness depends heavily on site-specific geological conditions, such as subsurface connectivity and aquifer characteristics, which require detailed investigation. Future research should refine hydrogeological models, conduct performance and economic analyses, and address regulatory and permitting frameworks to support broader adoption of this promising technology.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Growth and Stability of Stratiform Carrollite (CuCo2S4) in the Tenke-Fungurume Ore District, Central African Copperbelt

Abstract Carrollite (CuCo2S4) is the main ore mineral of the critical battery metal cobalt, yet, surprisingly, detailed characterization of its thermodynamic properties and of its trace element contents remain sparse in the scientific literature. To fill these knowledge gaps, we generated the first thermodynamic data set for carrollite using well-characterized samples obtained from the stratiform sulfide mineralization style at Fungurume 88 deposit in the Democratic Republic of Congo. Detailed mineralogical investigation of these same samples suggests that carrollite has limited ability to incorporate other ‘sweetener’ elements (except Cu and Ni) by substitution mechanisms. Modeling using CHNOSZ 2.0.0. software and an integrated Cu-Co-S-O thermodynamic data set provides insights into the stability of carrollite relative to other Co-bearing minerals under physicochemical conditions relevant to the Central African Copperbelt. Further exploration of these models indicates that dropping redox potential in the mineralizing system is insufficient to explain the mineral-scale and deposit-scale zonation patterns commonly observed in Cu-Co deposits of the Central African Copperbelt. Instead, factors such as increasing pH, decreasing T, decreasing sulfur activity, or some combination of changing physicochemical parameters may need to be invoked to explain these zonation trends. At the Fungurume 88 deposit, the absence of linnaeite (Co3S4) suggests that ore precipitation occurred at temperatures cooler than 211 °C. From our pH-logfO2 diagrams, and relative to higher temperatures, such a lower T regime is marked by a greater degree of offset between the modeled Cu and Co solubility contours, thus facilitating a more differentiated and zoned distribution of these two metals at the deposit scale. Our models are offered to the geoscience community as exceptionally useful tools for modeling and visualizing the Cu-Co-S-O systematics and may be applied to a variety of geological and geometallurgical questions in the Central African Copperbelt and other deposits elsewhere.

Mineralogy↗
DOE OSTI2021

Distributed Thermal Response Multi-Source Modeling to Evaluate Heterogeneous Subsurface Properties

A thorough assessment of thermal properties in heterogeneous subsurface is necessary in design of low-temperature borehole heat exchangers (BHEs). For this study, a distributed thermal response test (DTRT), which combines distributed temperature sensing (DTS) with a conventional thermal response test (TRT), was conducted in a U-bend geothermal loop installed in an open borehole at the University of Illinois at Urbana-Champaign to estimate thermal properties by analyzing the thermal response of different geologic materials while applying a constant heat input rate. Fiber-optic cables in the DTRT were deployed both inside the U-bend geothermal loop and in the center of the borehole to improve the accuracy of calculated heat-loss rates and borehole temperature profile measurements. To assess the subsurface thermal conductivity during the heating phase of the DTRT, a single-source model and a multi-source model, both based on the infinite line source method, were developed using the borehole temperature data and temperatures inside and along the outside of the loop, separately. The two models returned similar thermal conductivity values. The multi-source modeling has the advantage of predicting the thermal conductivity of heterogeneous geologic materials from borehole temperature profiles during the DTRT heating phase. Additionally, based on the distributed thermal conductivity measured in the borehole, estimates were made for both radial thermal impacts and the rate of heat loss in the BHE.

58 GEOSCIENCES↗
DOE OSTI2025

Combining Eddy Covariance Towers, Field Measurements, and the MEMS 2 Ecosystem Model Improves Confidence in the Climate Impacts of Bioenergy With Carbon Capture and Storage

ABSTRACT Carbon dioxide removal technologies such as bioenergy with carbon capture and storage (BECCS) are required if the effects of climate change are to be reversed over the next century. However, BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty. We place these net ecosystem carbon balance forcings in the context of associated radiative forcings from changes to soil N 2 O and CH 4 fluxes, land surface albedo, embedded fossil fuel use, and geologically stored carbon. Results from long term eddy covariance measurements, a soil and plant carbon inventory, and the MEMS 2 process‐based ecosystem model all agree that establishing perennials such as switchgrass or mixed prairie on former cropland resulted in net negative radiative forcing (i.e., global cooling) of −26.5 to −39.6 fW m −2 over 100 years. Establishing these perennials on former grassland sites had similar climate mitigation impacts of −19.3 to −42.5 fW m −2 . However, the largest climate mitigation came from establishing corn for BECCS on former cropland or grassland, with radiative forcings from −38.4 to −50.5 fW m −2 , due to its higher plant productivity and therefore more geologically stored carbon. Our results highlight the strengths and limitations of each method for quantifying the field scale climate impacts of BECCS and show that utilizing multiple methods can increase confidence in the final radiative forcing estimates.

Falvo, Grant [Department of Plant, Soil and Microb↗
DOE OSTI2020

Characterizing Baselines and Change in Gas Hydrate Systems using EM Methods

The objective of this project was to advance our understanding of gas hydrate systems in nature by characterizing their electrical properties in the field and in the laboratory. In the laboratory measurements, methane hydrate was synthesized from pure water ice and flash frozen seawater, with varying amounts of sand or silt added. Electrical conductivity was determined by impedance spectroscopy, using equivalent circuit modeling to separate the effects of electrodes and to gain insight into conduction mechanisms. Silt and sand increase the conductivity of pure hydrate, inferred to be contaminant NaCl contributing to conduction in hydrate, to a peak conductivity in agreement with peak resistivities observed in well logs through massive hydrate (3,000--10,000 Ωm). The addition of silt and sand lowers the conductivity of hydrate synthesized from seawater, by an amount consistent with Archie's Law. All samples were characterized using cryogenic scanning electron microscopy and energy dispersive spectroscopy, which shows good connectivity of salt and brine phases. Electrical conductivity measurements of pure hydrate and hydrate mixed with silt during pressure-induced dissociation supports previous conclusions that sediment increases dissociation rate. In order to characterize gas hydrate systems in the field, we collected 360 line kilometers of controlled-source electromagnetic data on Walker Ridge 313, Orca Basin (WR100), Mad Dog (GC781), and Green Canyon 955 in the Gulf of Mexico, all areas with known or seismically inferred gas hydrate deposits and which have be drilled or targeted for future drilling. We deep-towed an EM transmitter that generates an alternating electric field which propagates through the seafloor geology. Data were recorded on 6 receivers towed behind the transmitter at distances between 550 and 1550 m. In the presence of conductive geology, the electric fields will be attenuated, and conversely, in resistive geology the fields will be preserved. Our data were inverted using a 2D inversion method that first optimizes the model-data misfit, then finds the smoothest model fitting the data. This ensures that resistivity structures present in the final model are likely necessary. At each of the proposed drilling sites we found increased resistivity, interpreted as increased hydrate concentrations. However, not only were the primary drilling sites not always more resistive than the alternate sites, at WR313 the strongest resistors were not at the locations targeted for drilling.

03 NATURAL GAS↗
DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗
DOE OSTI2025

Toward Drilling the Perfect Geothermal Well: An International Research Coordination Network for Geothermal Drilling Optimization Supported by Deep Machine Learning and Cloud Based Data Aggregation

The EDGE project, supported by the U.S. Department of Energy Geothermal Technologies Office under award DE-EE0008793, established a data-driven framework for improving the efficiency, cost-effectiveness, and reliability of geothermal well drilling. The project focused on developing scalable data infrastructure, advanced machine learning and probabilistic models, and integrated analytics tools to support continuous drilling optimization. A central objective was to reduce geothermal drilling costs by up to seventy percent while minimizing the risk of well failure through predictive diagnostics and adaptive planning. Over the project period, a comprehensive data repository was designed and deployed, incorporating records from over one hundred geothermal wells across varied geological settings. This repository supported both structured and unstructured data and adhered to FAIR data principles, enabling provenance tracking, quality control, and standardized metadata. The project introduced automated ingestion pipelines and a cloud-hosted platform that facilitated access to raw, processed, and derived datasets. This infrastructure served as the foundation for model development and analysis. Machine learning workflows were developed to predict key drilling metrics including rate of penetration, non-productive time, and total drilling costs. Self-organizing maps and dimensionality reduction methods were used to uncover operational patterns and outliers, while supervised learning algorithms such as random forests and deep neural networks were applied to forecast performance outcomes. The models were validated on heterogeneous datasets from both U.S. and Icelandic fields, demonstrating variable but significant predictive accuracy. The results indicated that finer temporal resolution, inclusion of lithological data, and consistency in operational annotations could substantially improve model performance. The project also implemented process mining techniques to reconstruct state-transition models from drilling event logs. These models enabled the identification of deviations from optimal workflows and provided insights into recurring failure modes. Analysis of non-productive time highlighted the impact of equipment failures, geological challenges, and human factors, offering opportunities for targeted mitigation strategies. The EDGE Dashboard was developed as a web-based expert system integrating data visualization, model outputs, and user-driven queries. It provided an accessible interface for operators to explore historical data, evaluate predicted outcomes, and compare drilling scenarios. Initial feedback from project partners suggested that the dashboard could serve as a foundation for more advanced advisory and optimization tools. Overall, the EDGE project demonstrated the feasibility and value of applying modern data science techniques to geothermal drilling. It delivered a set of interoperable tools and models that can support more efficient, lower-risk well development. The findings point toward a viable path for transitioning from advisory analytics to semi-autonomous drilling systems, contingent on continued collaboration, expanded datasets, and field validation. The project results have immediate relevance for drilling operations, data management practices, and future geothermal R&D efforts aimed at achieving reliable, cost-competitive geothermal energy at scale.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Geologic hydrogen: From natural occurrences to anthropogenic generation – A review of fundamentals, potential, challenges and prospects

Growing demand for hydrogen is exposing the environmental and economic limits of reforming-based and carbon-managed supply chains, while the scale-up of electrolytic capacity remains capital-constrained. Geologic hydrogen, defined as molecular H₂ generated and stored within the Earth's crust offers a complementary, potentially lower-cost resource, yet exploration is still ad hoc. This review (1) revisits a global inventory of confirmed hydrogen seeps and subsurface occurrences; (2) analyzes the controlling reactions, migration pathways, and trapping conditions governing these occurrences; (3) proposes a process-based geologic hydrogen system concept analogous to, yet distinct from, the petroleum system; and (4) evaluates potential geologic hydrogen systems within the United States as a representative case study. Here, we contrast natural systems powered by serpentinization, mantle degassing or radiolysis with anthropogenic systems that stimulate the same reactions or convert in-situ hydrocarbons. Stable hydrogen accumulations require generation rates that outpace combined physical, chemical and microbial losses; the Bourakébougou field (Mali) exemplifies a self-recharging, free-gas reservoir sustained by meteoric-water serpentinization beneath an efficient caprock. Prospective geologic hydrogen resources are likely to occur in regions where iron-rich lithologies, deep-seated faults, and low-permeability sealing formations coexist. Applying this principle, we highlight three promising hydrogen play types in U.S. geological terrains: ophiolite belts (Appalachian and Californian regions), the Midcontinent Rift and the Lake Superior banded‑iron formations. Multiphysics numerical models and positive-unlabeled machine-learning workflows help to accelerate play screening and de-risk future production; yet, reaction kinetics, stimulation strategies, and full techno-economic and life-cycle assessments remain pivotal knowledge gaps.

Anthropogenic hydrogen generation↗
DOE OSTI2022

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗