Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geothermal fields”

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 253 records · Page 14

GOOML - Real World Applications of Machine Learning in Geothermal Operations

GOOML (Geothermal Operational Optimization with Machine Learning) is a machine-learning based framework that enables geothermal power plant operators to explore optimization opportunities for their assets in an efficient and robust digital environment. Backed by real-world data sources, thermodynamic constraints and steamfield intelligence, the GOOML environment provides new tools to explore how to best operate steamfields as well as test new scenarios and configurations prior to implementation in the field. To prove the effectiveness of GOOML, we have undertaken optimization experiments using reinforcement learning (RL) to generate operational suggestions using a balance of mass-take targets, sustainability considerations and net generation. Our experiments use the GOOML construct to explore different field parameters and perform multiple reinforcement learning experiments. Like a comprehensive laboratory workbench, we can change out components of a steamfield to perform testing under a variety of conditions (restrict mass, increase pressure, reroute steam, etc.). This flexibility allows us to explore conditions that would require significant infrastructure changes in a real-world setting at a fraction of the cost and time in a digital environment. The results highlight the benefits of using digital twins and advanced data analytics for the geothermal industry.

forecasting↗

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

All-digital Sensor System for Distributed Downhole Pressure Monitoring in Unconventional Fields

This project developed and validated (through field tests) a new low-cost all-digital pressure sensing technology for in situ distributed downhole pressure monitoring in unconventional oil and gas (UOG) fields. The all-digital sensing technology uses a built-in non-electric analog-to-digital converter (ADC) to transform the pressure information into a combination of binary (ON/OFF) states. As such, the system does not need downhole electronics for signal conditioning and telemetry. The all-digital sensors can be remotely logged over a long distance, and many sensors can be multiplexed for distributed sensing. Based on a review of unconventional wells in the Lower 48 states, the specification of the sensor is to measure pressure up to 69 MPa (10,000 psi) and temperature up to 250°C. A sensor with a helical bourdon sensing element and a digital signal decoder of 50 mm diameter and 109 mm length was constructed. The helical bourdon sensing element was made of 304L stainless steel and filled with motor oil. The digital converter was made up of 8 digital reading pads constructed of high-temperature epoxy with conductive inserts made of stainless steel. The sensor had a linear response to pressure with an accuracy of 0.14 MPa (20 psi). To withstand the high pressure, the sensor was enclosed in a stainless-steel pressure housing with a wall thickness of 5.5 mm, a diameter of 73 mm, and a length of 724 mm. In the laboratory tests, the sensor exhibited no temperature-related effects on the results. The sensor did not show drift over a 14-day test period at elevated pressure. A field test was conducted where the sensor was deployed in a test wellbore at the Quest drilling test facility to a depth of 0 feet over three weeks. The sensor was attached to the production rods, along with a downhole reference sensor of PPS27 type, which is a permanent downhole monitoring system. During the testing phase, the test well annular blow-out preventer was closed, and the well was pressurized at the surface to 11 MPa (1600 psi). The sensor read the elevated bottom hole pressure of 1500 psi. A multiplexing unit was created for the sensor to deploy multiple sensors on one data transmission line in a distributed approach. The multiplexing unit was tested in a simulated environment of 3048 m (10,000 ft) with five sensors distributed. The sensors were pressurized at different intervals. The multiplexed sensors recorded the correct pressure, and the multiplexing did not interfere with the readings. The proposed concept of an all-digital pressure sensor for harsh downhole environments was designed, manufactured, and tested in the laboratory and tested at the field to a up to 69 MPa and 250°C. This technology has high-temperature tolerance and has potential in downhole areas outside oil and gas, such as carbon capture and storage (CCS) and geothermal wells. The sensor concept has been proven in this project, but to create a commercially viable product, manufacturing a sensor with a smaller diameter needs to be performed.

02 PETROLEUM↗

Grid-Scale, Long-Term Energy Storage: Repurposing Hydrocarbon Reservoirs, Resources, and Infrastructure to Store CO 2 and Heat (Mid-Project Progress Report)

The purpose of this project is to create a viable business case for CO 2 capture, use, and storage (CCUS) by enabling grid-scale, long-term energy storage in mature hydrocarbon reservoirs and saline aquifers. Our CCUS approach, called the Earth Battery, is based on patents issued to Lawrence Livermore National Laboratory on Multi-Fluid Geo-Energy and Earth Battery Systems and to the University of Minnesota and TerraCOH, Inc. on CO 2 Plume Geothermal (CPG). Excess energy can be used to compress and store CO 2 and hot brine underground in permeable formations, including: (1) saline aquifers suitable for geologic CO 2 storage (GCS) and (2) hydrocarbon reservoirs, the focus of this project featuring a pilot test in a mature oil and natural gas (NG) field. Excess energy can come in the form of excess electricity from the grid or excess NG that does not have a market. It can also come in the form of excess heat from other sources, such as from the heat of air or CO 2 compression or from a nearby base-load power plant or solar thermal energy (STE) facility, when heat is not needed to generate electricity from those sources.

15 GEOTHERMAL ENERGY↗

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗

Pilgrim Hot Springs: GEOPHIRES Inputs and Outputs for Direct-Use Geothermal District Heating and Cooling

This dataset includes files for a techno-economic analysis conducted using the GEOPHIRES simulator to examine the feasibility of expanding a larger district heating site in a remote location: Pilgrim Hot Springs, Alaska. Files included here are GEOPHIRES inputs and outputs for five different scenarios with varying demand, cycle, and system design characteristics to analyze. Also included is the link to the GEOPHIRES GitHub, as well as a link to the dataset that contains the energy modelling used to determine the heating demand for the district. For a list of the differences between scenarios, see the included "Input Overview.txt" file. Fields included in the input files are: subsurface technical parameters, surface technical parameters, financial parameters, capital and O&M parameters, as well as simulation parameters. The output files are case reports that summarize all equipment, reservoir characteristics, costs, and heating profiles.

15 GEOTHERMAL ENERGY↗

Detecting fractures and monitoring hydraulic fracturing processes at the first EGS Collab testbed using borehole DAS ambient noise

Enhanced geothermal systems (EGS) require cost-effective monitoring of fracture networks. We validate the capability of using borehole distributed acoustic sensing (DAS) ambient noise for fracture monitoring using core photos and core logs. The EGS Collab project has conducted 10 m scale field experiments of hydraulic fracture stimulation using 50–60 m deep experimental wells at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The first EGS Collab testbed is located at 1616.67 m (4850 ft) depth at SURF and consists of one injection well, one production well, and six monitoring wells. All wells are drilled subhorizontally from an access tunnel called a drift. The project uses a single continuous fiber-optic cable installed sequentially in the six monitoring wells to record DAS data for monitoring hydraulic fracturing during stimulation. We analyze 60 s time records of the borehole DAS ambient noise data and compute the noise root-mean-square (rms) amplitude on each channel (points along the fiber cable) to obtain DAS ambient noise rms amplitude depth profiles along the monitoring wellbore. Our noise rms amplitude profiles indicate amplitude peaks at distinct depths. We compare the DAS noise rms amplitude profiles with borehole core photos and core logs and find that the DAS noise rms amplitude peaks correspond to the locations of fractures or lithologic changes indicated in the core photos or core logs. We then compute the hourly DAS noise rms amplitude profiles in two monitoring wells during three stimulation cycles in 72 h and find that the DAS noise rms amplitude profiles vary with time, indicating the fracture opening/growth or closing during the hydraulic stimulation. Our results demonstrate that borehole DAS passive ambient noise can be used to detect fractures and monitor fracturing processes in EGS reservoirs.

58 GEOSCIENCES↗

Selecting Appropriate Model Complexity: An Example of Tracer Inversion for Thermal Prediction in Enhanced Geothermal Systems

Abstract A major challenge in the inversion of subsurface parameters is the ill‐posedness issue caused by the inherent subsurface complexities and the generally spatially sparse data. Appropriate simplifications of inversion models are thus necessary to make the inversion process tractable and meanwhile preserve the predictive ability of the inversion results. In this study, we investigate the effect of model complexity on fracture aperture inversion and thermal performance prediction in a field‐scale EGS model. Principal component analysis was used to map the aperture field to a low‐dimensional latent space. The complexity of the inversion model was quantitatively represented by the percentage of total variance in the original aperture fields preserved by the latent space. Tracer, pressure and flow rate data were used to invert for fracture aperture through an ensemble‐based inversion method, and the inferred aperture field was used to predict thermal performance. With an over‐simplified aperture model, ensemble collapse occurred. The inverted aperture models failed to resolve necessary flow and transport features, leading to a biased thermal performance prediction. A complex aperture model involved excessive features and was prone to overinterpreting the inversion data. Both the tracer/pressure/flow rate data reproduction and thermal prediction showed significant uncertainties, making it difficult to properly estimate long‐term thermal performance. Fortunately, our results indicate that there exists an appropriate model complexity which can simultaneously match inversion data and predict thermal performance with an acceptable uncertainty. The quality of the fit of tracer data appears to be a useful indicator of such an appropriate model complexity.

15 GEOTHERMAL ENERGY↗

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗

Machine Learning‐Assisted Microearthquake Location Workflow for Monitoring the Newberry Enhanced Geothermal System

Abstract Enhanced geothermal systems (EGS) offer a sustainable energy source but face challenges in accurately locating microearthquakes induced during reservoir stimulation. Locating these microearthquakes provides reliable feedback on the stimulation progress. Current deep learning methods for locating earthquakes require extensive data sets for training, which is problematic as detected microearthquakes are often limited. To address the scarcity of training data, we propose a practical workflow using probabilistic multilayer perceptron (PMLP) which predicts microearthquake locations from cross‐correlation time lags in waveforms. Utilizing a 3D velocity model of Newberry site derived from ambient noise interferometry, we generate numerous synthetic microearthquakes and 3D acoustic waveforms for PMLP training. Accurate synthetic tests prompt us to apply the trained network to the 2012 and 2014 stimulation field waveforms. To enhance the accuracy of source localization, we carefully handpick the P‐arrival times. Predictions on the 2012 stimulation data set show major microseismic activity at depths of 0.5–1.2 km, correlating with a known casing leakage scenario. In the 2014 data set, the majority of predictions concentrate at 2.0–2.9 km depths, consistent with results obtained from conventional physics‐based inversion, and align with the presence of natural fractures from 2.0 to 2.7 km. We validate our findings by comparing the synthetic and field picks, demonstrating a satisfactory match for the first arrivals. By combining the benefits of quick inference speeds and accurate location predictions, we demonstrate the feasibility of using realistic synthetic data set to locate microseismicity for EGS monitoring.

15 GEOTHERMAL ENERGY↗

The Lawrence Berkeley Laboratory geothermal program in northern Nevada

The Lawrence Berkeley Laboratory's geothermal program began with consideration of regions where fluids in the temperature range of 150 to 230 C may be economically accessible. Three valleys, located in an area of high regional heat flow in north central Nevada, were selected for geological, geophysical, and geochemical field studies. The objective of these ongoing field activities is to select a site for a 10-MW demonstration plant. Field activities (which started in September 1973) are described. A parallel effort has been directed toward the conceptual design of a 10-MW isobutane binary plant which is planned for construction at the selected site. Design details of the plant are described. Project schedule with milestones is shown together with a cost summary of the project.

Mirk, K. F.↗

A Frequency-Domain-Based Algorithm for Detecting Microseismicity Using Dense Surface Seismic Arrays

We propose a new frequency-domain-based algorithm for detecting small-magnitude seismic events using dense surface seismic arrays. Our proposed method takes advantage of the high energy carried by S waves, and approximate known source locations, which are used to rotate the horizontal components to obtain the maximum amplitude. By surrounding the known source area with surface geophones, we achieve a favorable geometry for locating the detected seismic events with the backprojection method. To test our new detection method, we used a dense circular array, consisting of 151 5 Hz three-component geophones, over a 5 km aperture that was in operation at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) in southcentral Utah. We apply the new detection method during a small-scale test injection phase at FORGE, and during an aftershock sequence of an Mw 4.1 earthquake located ~30 km north of the geophone array, within the Black Rock volcanic field. We are able to detect and locate microseismic events (Mw<0) during injections, despite the high level of anthropogenic activity, and several aftershocks that are missing from the regional catalog. By comparing our method with known algorithms that operate both in the time and frequency domain, we show that our proposed method performs better in the case of the FORGE injection monitoring, and equally well for the off-array aftershock sequence. Our new method has the potential to improve microseismic event detections even in extremely noisy environments, and the proposed location scheme serves as a direct discriminant between true and false detections.

Geochemistry & Geophysics↗

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗

Shallow Geothermal Resources for Cooling Applications at the University of Hawai'i

Drilling activities account for 30% to 57% of the cost to develop and install a geothermal plant. Therefore, an accurate representation of the cost to drill a well is paramount in techno-economic analysis to determine the feasibility of a geothermal power project. In 2022, the National Renewable Energy Laboratory (NREL) endeavored to revise the U.S. Department of Energy (DOE) GeoVision baseline drilling cost curves due to extensive improvement in drilling rates at the Utah Frontier Observatory Research in Geothermal Energy (FORGE) demonstration site. That effort did not culminate in the recommendation of new curves because the actual project costs did not match the reported performance improvements and were at or above the GeoVision baseline. The need for another iteration of this analysis has arisen from industry record drilling performance reported by recent commercial field-scale and demonstration projects, including Fervo Energy's Cape Station, the Utah FORGE 16B(78)-32 demonstration and the Geysers Power Company's GDC-36 demonstration. Therefore, in this work, we have estimated the resulting industry average rate of penetration (ROP) and bit life and applied these parameters as inputs to the Well Cost Simplified model used in the GeoVision analysis. The resulting revised cost curves show a significant decline from the GeoVision baseline. For vertical wells, the magnitude of cost decline ranges between 12% and 24% while for deviated wells, cost reductions between 18% and 26% are estimated. The revised cost curves are in good agreement with actual cost data and therefore, quantify the economic impact of the utilization of (and advances in) polycrystalline diamond compact (PDC) bit technology and the application of physics-based methodologies that optimize mechanical specific energy.

building cooling↗

2025 Geothermal Drilling Cost Curves Update: Preprint

Drilling activities account for 30% to 57% of the cost to develop and install a geothermal plant. Therefore, an accurate representation of the cost to drill a well is paramount in techno-economic analysis to determine the feasibility of a geothermal power project. In 2022, the National Renewable Energy Laboratory (NREL) endeavored to revise the U.S. Department of Energy (DOE) GeoVision baseline drilling cost curves due to extensive improvement in drilling rates at the Utah Frontier Observatory Research in Geothermal Energy (FORGE) demonstration site. That effort did not culminate in the recommendation of new curves because the actual project costs did not match the reported performance improvements and were at or above the GeoVision baseline. The need for another iteration of this analysis has arisen from industry record drilling performance reported by recent commercial field-scale and demonstration projects, including Fervo Energy’s Cape Station, the Utah FORGE 16B(78)-32 demonstration and the Geysers Power Company’s GDC-36 demonstration. Therefore, in this work, we have estimated the resulting industry average rate of penetration (ROP) and bit life and applied these parameters as inputs to the Well Cost Simplified model used in the GeoVision analysis. The resulting revised cost curves show a significant decline from the GeoVision baseline. For vertical wells, the magnitude of cost decline ranges between 12% and 24% while for deviated wells, cost reductions between 18% and 26% are estimated. The revised cost curves are in good agreement with actual cost data and therefore, quantify the economic impact of the utilization of (and advances in) polycrystalline diamond compact (PDC) bit technology and the application of physics-based methodologies that optimize mechanical specific energy.

15 GEOTHERMAL ENERGY↗

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗

GOOML Big Kahuna Forecast Modeling and Genetic Optimization Files

This submission includes example files associated with the Geothermal Operational Optimization using Machine Learning (GOOML) Big Kahuna fictional power plant, which uses synthetic data to model a fictional power plant. A forecast was produced using the GOOML data model framework and fictional input data, and a genetic optimization is included which determines optimal flash plant parameters. The inputs and outputs associated with the forecast and genetic optimization are included. The input and output files consist of data, configuration files, and plots. A link to the Physics-Guided Neural Networks (phygnn) GitHub repository is also included, which augments a traditional neural network loss function with a generic loss term that can be used to guide the neural network to learn physical or theoretical constraints. phygnn is used by the GOOML framework to help integrate its machine learning models into the relevant physics and engineering applications. Note that the data included in this submission are intended to provide a demonstration of GOOML's capabilities. Additional files that have not been released to the public are needed for users to run these models and reproduce these results. Units can be found in the readme data resource.

15 GEOTHERMAL ENERGY↗