Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geothermal Systems”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2023

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

The direct-current response of electrically conducting fractures excited by a grounded current source

Hydraulic fracture stimulation of low permeability reservoir rocks is an established and cross–cutting technology for enhancing hydrocarbon production in sedimentary formations and increasing heat exchange in crystalline geothermal systems. Whereas the primary measure of success is the ability to keep the newly generated fractures sufficiently open, long–term reservoir management requires a knowledge of the spatial extent, morphology, and distribution of the fractures — knowledge primarily informed by microseismic and ground deformation monitoring. To minimize the uncertainty associated with interpreting such data, we investigate through numerical simulation the usefulness of direct-current (DC) resistivity data for characterizing subsurface fractures with elevated electrical conductivity by considering a geophysical experiment consisting of a grounded current source deployed in a steel cased borehole. In doing so, the casing efficiently energizes the fractures with steady current. Finite element simulations of this experiment for a horizontal well intersecting a small set of vertical fractures indicate that the fractures manifest electrically in (at least) two ways: (1) a local perturbation in electric potential proximal to the fracture set, with limited farfield expression and (2) an overall reduction in the electric potential along the borehole casing due to enhanced current flow through the fractures into the surrounding formation. The change in casing potential results in a measurable effect that can be observed far from fractures themselves. Under these conditions, our results suggest that farfield, timelapse measurements of DC potentials can be interpreted by simple, linear inversion for a Coulomb charge distribution along the borehole path, including a local charge perturbation due to the fractures. As a result, this approach offers an inexpensive method for detecting and monitoring the time-evolution of electrically conducting fractures while ultimately providing an estimate of their effective conductivity — the latter providing an important measure independent of seismic methods on fracture shape, size, and hydraulic connectivity.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Project Phase 1 Report: Reducing Data Center Peak Cooling Demand and Energy Costs With Cold Underground Thermal Energy Storage (Cold UTES)

Cold Underground Thermal Energy Storage (Cold UTES) is an ultra-long duration grid energy storage technology. With Cold-UTES, low-cost grid power is converted to cold thermal energy and stored in the native subsurface rock at the point of use. Cold UTES is one approach within the general category of engineered geothermal systems. Cold UTES for peak-hour cooling of data centers (DCs) was studied for deployment in Maricopa County, Arizona and Loudoun County, Virgina using thermal storage capacities from 4 GWh-th to over 1,000 GWh-th (>1 Terawatt-hour). The two sites have different power and transmission systems, available grid energy resources, daily and seasonal load profiles, weather conditions, and grid regulatory requirements. The study results indicate high value for both locations and because of this, likely indicates value across most of the US and the world. The basis of the study was a 1,000 MW-e hourly electric use of DC computing and auxiliary loads, which was modeled as 1 GW-th of thermal load to a dry cooled heat rejection system - i.e. a cooling system that does not consume water. The electric power required for cooling the DC varies as the air temperature changes. In cold weather only the dry-coolers are used, with an electrical load for cooling load as low as 10 MW-e. In hot summer hours, chillers and dry- coolers are required, which raises the electrical load for cooling load to as much as 300 MW-e. The continuous and peak cooling electrical loads result in a grid interconnection requirement of no less than 1,300 MW-e. From both a grid and thermal design modeling perspective the 1.3 GW-e could either be a single facility or result from the total load at multiple sites.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development

This DOE GTO funded project was a collaborative effort between Array Information Technology and the University of California at Berkeley. The project goal is to improve technology to assess in-situ stress magnitude and orientation, kinematic fracture parameters, rupture size, rupture orientation as well as temporal changes and volumetric distribution of the injected fluid during Enhanced Geothermal System (EGS) resource development. The project leveraged high-frequency seismic data recorded by the LBNL 34-station permanent geophone network and seismic broadband data recorded by a temporary 33-station seismometer network that operated in The Geysers during the injection phase of the DOE GTO funded EGS demonstration project at the well Prati-32. The operation of the broadband network with high station density and good azimuthal coverage to monitor the injection phase produced an unprecedented dataset that is typically unavailable for EGS operations. The two datasets offer the advantages of studying the parameters under investigation over a broader frequency band and the higher sensitivity of the broadband seismometers resulted in a richer dataset with a lower magnitude of completeness. Moment tensor (MT) inversion have been employed to develop a catalog of earthquake source solutions including moment tensor and finite source rupture area. This resulted in a 168-event waveform-based seismic moment tensor catalog for events 0.7≤Mw≤3.9 based on a semi-automated moment tensor analysis technique. The MT solutions were used for in-situ stress estimation during the EGS injection phases. We found an approximate 15-degree counterclockwise rotation of the least compressive stress sigma-3, and a rotation of the maximum compressive stress sigma-1 toward the vertical as the injected volume of water increased. We developed a rupture-area magnitude scaling relationship of The Geysers earthquakes obtained from finite-source inversion for fault slip that reveals a very high correlation to published results. Based on the rupture area-magnitude relationship we mapped Mw from the LBNL earthquake catalog to estimate the activated fracture area, its location and orientation within the injection volume of Prati-32. The double-difference Wadati (DDW) technique was applied to high-accuracy P- and S-wave differential travel times derived from waveform cross correlation to image fluid saturation in the subsurface based on Vp/Vs-ratio. The recently developed technique has been applied to estimate temporal changes of fluid distribution in fault zones and in volcanic regions but not yet in engineered geo-reservoirs. The improvement of the technology includes high-resolution waveform cross correlation applied to the high-quality datasets with lower magnitude of completeness resulting in a higher number of recorded events. The improvement allowed to estimate both, the spatial and temporal distribution of injected fluids and the state of the injected fluid during reservoir stimulation. Results of the DDW analysis revealed the development of an injection-derived steam plume during the early injection phase, particularly in the northern part of the Prati-32 study area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Integration of microseismic monitoring data into coupled flow and geomechanical models with ensemble Kalman filter

Abstract Hydraulic stimulation of low‐permeability rocks in enhanced geothermal systems, shale resources, and CO 2 storage aquifers can trigger microseismic events, also known as microearthquakes (MEQs). The distribution of microseismic source locations in the reservoir may reveal important information about the distribution of hydraulic and geomechanical rock properties. In this paper, we present a framework for conditioning heterogeneous rock permeability and geomechanical property distributions on microseismic data. To simulate the multiphysics processes in these systems, we combine a fully coupled flow and geomechanical model with the Mohr‐Coulomb type rock failure criterion. The resulting multiphysics simulation constitutes the forecast model that relates microseismic source locations to reservoir rock properties. We adopt this forward model in an ensemble Kalman filter (EnKF) data assimilation framework to jointly estimate reservoir permeability and geomechanical property distributions from injection‐induced microseismic response measurements. We show that integration of a large number of spatially correlated microseismic data with practical ensemble sizes can lead to severe underestimation of ensemble spread, and eventually ensemble collapse. To mitigate the variance underestimation issue, two low‐rank data representation schemes are presented and discussed. In the first approach, microseismic data are projected onto a low‐dimensional subspace defined by the left singular vectors of the perturbed observation matrix. The second method uses a coarser grid for representing the microseismic data. A series of numerical experiments is presented to evaluate the performance of the proposed methods and to illustrate their applicability for assimilating microseismic data into coupled flow and geomechanical forward models to estimate multiphysics rock properties.

Tarrahi, Mohammadali↗
DOE OSTI2016

An analytical thermohydraulic model for discretely fractured geothermal reservoirs

In discretely fractured reservoirs such as those found in Enhanced/Engineered Geothermal Systems (EGS), knowledge of the fracture network is important in understanding the thermal hydraulics, i.e., how the fluid flows and the resulting temporal evolution of the subsurface temperature. Here, the purpose of this study was to develop an analytical model of the fluid flow and heat transport in a discretely fractured network that can be used for a wide range of modeling applications and serve as an alternative analysis tool to more computationally intensive numerical codes. Given the connectivity and structure of a fracture network, the flow in the system was solved using a linear system of algebraic equations for the pressure at the nodes of the network. With the flow determined, the temperature in the fracture was solved by coupling convective heat transport in the fracture with one-dimensional heat conduction perpendicular to the fracture, employing the Green’s function derived solution for a single discrete fracture. The predicted temperatures along the fracture surfaces from the analytical solution were compared to numerical simulations using the TOUGH2 reservoir code. Through two case studies, we showed the capabilities of the analytical model and explored the effect of uncertainty in the fracture apertures and network structure on thermal performance. While both sources of uncertainty independently produce large variations in production temperature, uncertainty in the network structure, whenever present, had a predominant influence on thermal performance.

58 GEOSCIENCES↗
DOE OSTI2024

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗
DOE OSTI2014

Wellbore cement fracture evolution at the cement–basalt caprock interface during geologic carbon sequestration

Composite Portland cement-basalt caprock cores with fractures, as well as neat Portland cement columns, were prepared to understand the geochemical and geomechanical effects on the integrity of wellbores with defects during geologic carbon sequestration. The samples were reacted with CO2-saturated groundwater at 50 ºC and 10 MPa for 3 months under static conditions, while one cement-basalt core was subjected to mechanical stress at 2.7 MPa before the CO2 reaction. Micro-XRD and SEM-EDS data collected along the cement-basalt interface after 3-month reaction with CO2-saturated groundwater indicate that carbonation of cement matrix was extensive with the precipitation of calcite, aragonite, and vaterite, whereas the alteration of basalt caprock was minor. X-ray microtomography (XMT) provided three-dimensional (3-D) visualization of the opening and interconnection of cement fractures due to mechanical stress. Computational fluid dynamics (CFD) modeling further revealed that this stress led to the increase in fluid flow and hence permeability. After the CO2-reaction, XMT images displayed that calcium carbonate precipitation occurred extensively within the fractures in the cement matrix, but only partially along the fracture located at the cement-basalt interface. The 3-D visualization and CFD modeling also showed that the precipitation of calcium carbonate within the cement fractures after the CO2-reaction resulted in the disconnection of cement fractures and permeability decrease. The permeability calculated based on CFD modeling was in agreement with the experimentally determined permeability. This study demonstrates that XMT imaging coupled with CFD modeling represent a powerful tool to visualize and quantify fracture evolution and permeability change in geologic materials and to predict their behavior during geologic carbon sequestration or hydraulic fracturing for shale gas production and enhanced geothermal systems.

basalt↗
DOE OSTI2016

An overview of TOUGH-based geomechanics models

After the initial development of the first TOUGH-based geomechanics model 15 years ago based on linking TOUGH2 multiphase flow simulator to the FLAC3D geomechanics simulator, at least 15 additional TOUGH-based geomechanics models have appeared in the literature. This development has been fueled by a growing demand and interest for modeling coupled multiphase flow and geomechanical processes related to a number of geoengineering applications, such as in geologic CO 2 sequestration, enhanced geothermal systems, unconventional hydrocarbon production, and most recently, related to reservoir stimulation and injection-induced seismicity. This paper provides a short overview of these TOUGH-based geomechanics models, focusing on some of the most frequently applied to a diverse set of problems associated with geomechanics and its couplings to hydraulic, thermal and chemical processes.

58 GEOSCIENCES↗
DOE OSTI2017

Simulations of CO 2 injection into fractures and faults for improving their geophysical characterization at EGS sites

Here, we propose the use of CO 2 in push-pull well tests to improve geophysical identification and characterization of fractures and faults at enhanced geothermal system (EGS) sites. Using TOUGH2/ECO2N, we carried out numerical experiments of push-pull injection-production cycling of CO 2 into idealized vertical fractures and faults to produce pressure-saturation-temperature conditions that can be analyzed for their geophysical response. Our results show that there is a strong difference between injection and production mainly because of CO 2 buoyancy. While the CO 2 -plume grows laterally and upward during injection, not all CO 2 is recovered during the subsequent production phase. Even under the best conditions for recovery, at least 10% of the volume of the pores still remains filled with CO 2 . To improve EGS characterization, comparisons can be made of active seismic methods carried out before and after (time lapse mode) CO 2 injection into the fracture or fault. Here we find that across the CO 2 saturation range, C11 (the normal stiffness in the horizontal direction perpendicular to the fracture plane) varies between maximum and minimum values by about 15%. It reaches a maximum at around 6% gas saturation, decreasing exponentially to a minimum at higher saturations. Finally, our results suggest that CO 2 injection can be effectively used to infiltrate fault and fracture zones reaching about optimal saturation values in order to enhance seismic imaging at EGS sites.

58 GEOSCIENCES↗
DOE OSTI2024

Seismic monitoring of EGS fracture stimulations at Utah FORGE (Part 1): Time-lapse variations of b-values and Shear-wave Splitting Rates of induced microearthquakes

Characterizing geologic formations near geothermal wells is crucial for efficient extraction of subsurface geothermal energy resources. To monitor effective fracture stimulation of enhanced geothermal systems (EGS), we analyze a catalog of induced microearthquakes to characterize seismicity induced by Phase 2C fracture stimulations of the Utah FORGE project and derive the b-values for 424 induced microearthquakes in three stages of fracture stimulation. Here, we find the time-lapse variations of these b-values. We apply eigenvector rotation and waveform correlation to 3C geophone data to separate fast S 1 and slow S 2 waves, determine their directions, measure their lag time, and calculate the S-wave splitting rate (SSR) for six induced microearthquakes in Stage 1 and 21 in Stage 3. Both SSR values and b-values increase with time. The increase in these values may be related to heavily fractured zones created by fracture stimulations during Stage 2 injections. Large and highly variable SSR values in Stage 3 may reflect heterogeneities in the fracture zones. Time-lapse variations of SSR values and b-values may be one of the most useful approaches to monitoring and characterizing the effectiveness and efficiency of various types of fracture stimulations.

58 GEOSCIENCES↗
DOE OSTI2015

Uncultivated thermophiles: current status and spotlight on ‘Aigarchaeota’

Meta-analysis of cultivation-independent sequence data shows that geothermal systems host an abundance of novel organisms, representing a vast unexplored phylogenetic and functional diversity among yet-uncultivated thermophiles. There are a number of thermophiles that were recently interrogated using metagenomic and/or single-cell genomic approaches, including members of taxonomic groups that inhabit both thermal and non-thermal environments, such as 'Acetothermia' (OP1) and 'Atribacteria' (OP9/JS1), as well as the exclusively thermophilic lineages 'Korarchaeota', 'Calescamantes' (EM19), 'Fervidibacteria' (OctSpA1-106), and 'Aigarchaeota' (HWCG-I). The 'Aigarchaeota', a sister lineage to the Thaumarchaeota, likely includes both hyperthermophiles and moderate thermophiles. Furthermore, they inhabit terrestrial, marine, and subsurface thermal environments and comprise at least nine genus-level lineages, several of which are globally distributed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Ecology of methyl-coenzyme M reductase encoding Thermoproteota

The recent demonstration that members of at least three classes of archaea affiliated with the Thermoproteota superphylum are involved in the production of the climate-active gas methane has sparked discussions about how well we understand the diversity of methanogens. Here, in this work, we show that members of all three of these lineages, as well as several other, yet uncultured and physiologically uncharacterized groups within the Thermoproteota that encode the key enzyme of anaerobic methane cycling, methyl-coenzyme M reductase (MCR), are widely distributed in anoxic ecosystems. We postulate that the taxonomic, metabolic, and ecological diversity of methanogenic and MCR-encoding Thermoproteota are poorly understood, and that the contribution of methylotrophic and thermoproteotal methanogenesis to methane production is largely unknown. We hypothesize that thermoproteotal methanogens could contribute, potentially substantially, to methane emissions in many anoxic environments that harbor methylated precursors, including wetlands, sediments, peat, rice paddies, wastewater sludge, and geothermal systems. We highlight the necessity to experimentally test the (eco)physiology of these widely distributed archaea using both culture-dependent ( in vitro ) and culture-independent ( in situ ) approaches to assess their potential contribution to methane emissions. Last, we stress the importance of remaining agnostic about the physiology of MCR-encoding Thermoproteota in the absence of experimental data because most of these archaea also carry the genetic potential to grow non-methanogenically.

MCR↗
DOE OSTI2023

Modeling aqueous association constants and mineral solubilities at subcritical and supercritical temperatures

The need for sustainable power generation has increased interest in the use of hydrothermal fluids for industrial applications. New high-enthalpy geothermal systems and biowaste-to-fuel processes are two relevant examples that employ supercritical fluids which require an in-depth understanding of complex chemical reactions occurring near the supercritical temperature of water (374 °C). As these processes operate in thermodynamic regimes that are not currently covered by a standard molar Gibbs energy of formation model, only empirical fits for single reaction systems are available which limit the use of multi-component phase equilibria calculations that are standard practice for less extreme environments. Here, in this paper, we advance a standard molar Gibbs energy of formation model able to operate in these otherwise inaccessible thermodynamic states to include species needed for key mineral solubility systems and ion association reactions. This work extends a model based on molecular statistical thermodynamics (MST) into four new systems (Na 3 PO 4 -H 2 O, LiOH-H 2 O, KOH-H 2 O, and BaSO 4 -H 2 O) by extending the model to cover 10 new species. For each of these systems, model predictions were consistently within the experimental uncertainties for the new systems covered. A breakdown of MST contributions to the model revealed that electrostatic and hard sphere contributions were key to reproducing density dependencies of standard molar Gibbs energy of formation values around the critical point of water.

42 ENGINEERING↗
DOE OSTI2025

Nitrogen fixation in Pampinifervens , a new species-rich genus of Aquificaceae that inhabits a wide pH range in terrestrial hot springs

Members of the Aquificota play important roles in hydrogen, sulfur, and carbon cycles in geothermal systems as ubiquitous and abundant chemolithoautotrophs; however, their roles in the nitrogen cycle are poorly defined. Here, we show that Aquificaceae isolate T-2 T from Tengchong, China, fixes nitrogen based on growth with dinitrogen as the sole nitrogen source and incorporation of 15 N atoms when grown with 15 N 2 . We further show evidence that suggests its nitrogen fixation (nif) genes are transcribed in situ in a nearby spring. We propose the name Pampinifervens diazotrophicum gen. et sp. nov., with strain T-2 T (= JCM 35475 T = DSM 116324 T ) as the nomenclatural type. Its close relative designated Pampinifervens florentissimum sp. nov. T-8 T (= CGMC 1.5214 T = JCM 33569 T ) does not fix nitrogen but respires nitrate. Analysis of 104 metagenome-assembled genomes from 84 metagenomes allowed us to circumscribe the genus Pampinifervens and identify at least 16 sympatric Pampinifervens species that dominate the Aquificaceae populations in diverse hot springs in southwest China, with distinct species globally. Analysis of these Pampinifervens genomes revealed variable presence of nif and denitrification genes and frequent gene loss and horizontal transfer. Genes for pH homeostasis were also present in Pampinifervens genomes, with a proliferation in species from Tengchong, consistent with the wide pH growth range of both Pampinifervens isolates and the prevalence of Pampinifervens in both alkaline and acidic springs. The dynamic evolution of nif and denitrification genes, along with pH homeostasis genes, may facilitate the diversification of Pampinifervens into diverse springs with different nitrogen availability.

Aquificae↗
DOE OSTI2016

Numerical Simulation of Permeability Change in Wellbore Cement Fractures after Geomechanical Stress and Geochemical Reactions Using X-ray Computed Tomography Imaging

X-ray microtomography (XMT) imaging combined with a three-dimensional (3D) computational fluid dynamics (CFD) modeling technique was used to study the effect of geochemical and geomechanical processes on fracture properties in composite Portland cement–basalt caprock core samples. The effect of fluid properties and flow conditions on fracture permeability was numerically studied by using fluids with varying physical properties and simulating different pressure conditions. CFD revealed that the application of geomechanical stress led to increased fluid flow, which resulted in increased fracture permeability. After CO2-reaction, XMT images displayed preferential precipitation of calcium carbonate within the fractures in the cement matrix and less precipitation in fractures located at the cement–basalt interface. CFD predicted changes in flow characteristics and differences in absolute values of flow properties due to different pressure gradients. CFD was able to highlight the profound effect of fluid properties on flow characteristics and hydraulic properties of fractures. This study demonstrates the applicability of XMT imaging and CFD as powerful tools for characterizing the hydraulic properties of fractures in a number of applications like geologic carbon sequestration and storage, hydraulic fracturing for shale gas production, and enhanced geothermal systems.

CO2 storage↗
DOE OSTI2017

Partitioning of HCl Between Concentrated Brines and Supercritical CO 2 Under Geologic Reservoir Conditions

HCl is generated in subsurface brines exposed to high CO 2 pressures via reactions such as CO 2 + NaCl + H 2 O ↔ NaHCO 3 + HCl. The extent of partitioning of HCl between a concentrated chloride brine and supercritical CO 2 (scCO 2 ) has not been measured before and it cannot be estimated with the accuracy sufficient for predictive modeling of subsurface processes. The partitioning of HCl between 4.92 mol/kg aqueous solutions of NaCl and scCO 2 was measured at 100 and 150 °C and at pressures between 9 and 16.2 MPa. At P = 15 MPa, the concentrations of HCl in the scCO 2 phase were 7 × 10 –4 and 4 × 10 –3 mol/L at 100 and 150 °C, an increase by more than 5 and 3 orders of magnitude, respectively, relative to the concentrations in water vapor in the absence of CO 2 . Furthermore, these levels of HCl may accelerate reactions with reservoir rocks, natural seals (shale caprocks) and downhole materials, potentially enhancing or impairing operation and long-term performance of subsurface technologies such as geologic carbon sequestration/storage, enhanced geothermal systems using compressed CO 2 , and scCO 2 extraction of hydrocarbons, and other resources in the presence of chloride brines.

58 GEOSCIENCES↗