Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir simulation”

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 505 records · Page 28

Water and Climate Impacts on ERCOT Long-Term Systems Assessment

Water and climate change pose many potential challenges to the electric power system. Substantial water is withdrawn every day to support thermoelectric power generating unit operations, and changes to water supply have the potential to affect generation dispatch. Climate change can accelerate growing demand for electricity, which can necessitate additional generating capacity, often in locations with limited water supply. Drought conditions also threaten thermoelectric power plant operations due to streamflow and reservoir levels dropping below intake structures, or water temperatures exceeding a power plants' permitted operating conditions. Here we explore how future climate change might influence decisions related to electricity capacity expansion planning in Electric Reliability Council of Texas (ERCOT) using a multi-model framework. Specifically, water resource modeling is used to simulate climate impacts on the future water supply for thermoelectric and hydropower generation for four future climate projections. Separately, temperature impacts on electricity load are evaluated for these scenarios. These climate impacts are applied to five alternative electricity futures in an electricity capacity expansion model that projects future generation and transmission capacity additions in ERCOT. Results indicate that climate has a measurable influence on future generation and transmission capacity needs, with temperature-driven increases in peak and average load resulting in 5-15 GW additional generating capacity and up to 1 GW additional transmission capacity. Additional capacity is a diverse mix of PV, natural gas, and wind, depending on the makeup of economic, policy and technology assumptions. Climate impacts increase total system costs 2-5%, while the marginal cost of energy and emissions are not affected substantially by climate change effects.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Observations and Modeling of Fiber-Optics Strain on Hydraulic Fracture Height Growth in HFTS-2

Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced dataset to allow us to observe and understand fracture geometries rigorously. Low frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. Furthermore, the simulation results improved the understating of the HFTS-2 LF-DAS data and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

58 GEOSCIENCES↗

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗

CO 2 geological storage: Critical insights on plume dynamics and storage efficiency during long-term injection and post-injection periods

In the quest for net-zero greenhouse gas emission, CO 2 capture and storage technologies play a prominent role. Estimates suggest a potential of sequestering up to several thousand metric gigatons of CO 2 . However, some of the main challenges in its development are an assessment of reservoir-specific storage capacities due to variabilities and heterogeneities in the underlying properties (structure, formations, facies, petrophysics, reservoir architecture, and others) and understanding the migration of CO 2 in the subsurface. How much CO 2 can be stored in a reservoir, and how quickly can it be injected? What will be the dominant storing (trapping) mechanisms? These questions, among others, remain as challenges to the commercial deployment of this technology. This article critically investigates CO 2 plume characteristics and determines the evolving contributions of different trapping mechanisms during a 100-year injection period. Here, simulation results on the CO 2 plume dynamics indicate that lateral propagation is much larger than vertical propagation. The ratio of average lateral to vertical dispersion ranged between 6.3 and 18.8 for the cases investigated. On average, for all the scenarios investigated, the free phase (mobile supercritical CO 2 ) is the dominant trapping mechanism (in terms of the amount of CO 2 stored) during the injection period, accounting up to 45% of the CO 2 injected. Both residual (30%) and solubility (25%) trapping storage ratios decline with time during injection, and no significant mineral trapping appeared during the 100-year injection period. However, the behavior of CaCO 3 suggests a predominant dissolution of CaCO 3 , especially in the vicinity of the injection well.

03 NATURAL GAS↗

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗

Utah FORGE 5-2428: Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior - Workshop Presentation

This is a presentation on the Fracture Permeability Impact on Seismic Slip Behavior project by Lawrence Livermore National Laboratory, presented by Dr. Kayla A. Kroll. The project's objective is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Utah FORGE 5-2557: Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes - Workshop Presentation

This is a presentation on the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University, presented by Distinguished Professor of Physics & Astronomy, Laura J. Pyrak-Nolte. The project's objective was to develop and validate a macroscopic model that accounts for local deformation/frictional behavior, seismic/aseismic behavior, chemical reactions, and determine the adequacy of classic Coulomb failure vs. rate-and-state friction. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

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

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

58 GEOSCIENCES↗

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf galaxy evolution

ABSTRACT We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR), respectively, to study 78 radio active galactic nuclei (AGN) in nearby (z < 0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages, and accretion rates are ∼1035 W, ∼5 Myr, and ∼10−3.4 M⊙ yr−1, respectively. The median mechanical energy output of these sources is ∼100 times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability, and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf galaxy evolution, as is the case in massive galaxies.

Davis, F. (ORCID:0000000285224980)↗

Distributed electric field sensing using fibre optics in borehole environments

In the past decade, rapid advances in distributed optical fibre sensing technologies have made it possible to record various geophysical data (e.g. strain, temperature and pressure) continuously in both time and space along the fibre, providing an unprecedented quantity and spatial density of data compared to traditional geophysical measurements as well as reducing data acquisition cost. To date, no distributed fibre-based electromagnetic field sensing system has been implemented although electromagnetic sensing could have a broad range of applications to geophysical imaging and monitoring in borehole environments. The goal of this paper is to provide a theoretical feasibility study regarding the design and use of an electromagnetic sensing optical fibre for geophysical applications. First, we present the sensitivity analysis of a ‘hypothetical’ optical fibre coated with polyvinylidene fluoride, a polymer that provides relatively high piezoelectric properties, yet unlike ceramics, is flexible. Next, using a two-dimensional electromagnetic modelling algorithm, we simulate the earth electric-field-to-fibre-strain transfer function and estimate the theoretical sensitivity of the optical fibre to electric fields. Given the state-of-the-art distributed acoustic sensing strain sensitivities in the picometres strain range, our numerical modelling analysis suggests that a perfectly coupled polyvinylidene fluoride–coated optical fibre can measure electric field values in the mV/m to V/m amplitude range. We then apply a cylindrically symmetric modelling algorithm to simulate numerical models demonstrating the applicability of such a fibre in an oilfield environment. Scenarios investigated employ an electric field source and suggest that the measurements can be used to distinguish the oil versus water ratio with a fibre mounted inside a producing steel cased oil well as well as distinguishing between brine and hydrocarbon filled reservoir zones with a fibre located outside of the casing.

58 GEOSCIENCES↗

CO 2 /Brine/Rock Interactions in the Cedar Keys-Lawson Formation

It has been reported that among the various geologic storage options, deep saline aquifers have the largest estimated capacity for CO 2 storage. Obtaining knowledge of possible artificially geochemically induced changes to the permeability and porosity of host CO 2 storage deposits will enable us to gain insight on long-term reservoir behavior under CO 2 storage conditions. An experimental study of the interaction of CO 2 /brine/rock on saline formations was conducted in a static system under CO 2 storage conditions. Chemical interactions in the Cedar Keys-Lawson Formation carbonate during exposure to CO 2 and brine under sequestration conditions were studied. Samples were exposed to the simulated in-situ reaction conditions for one and six months. The samples were exposed to a model brine at 55 °C and CO 2 pressure of 23.8 MPa (3,500 psig). Computed tomography (CT), x-ray diffraction (XRD), scanning electron microscopy (SEM)-energy dispersive x-ray spectroscopy (EDS), brine composition, core porosity, and core permeability analyses were conducted prior to and after the exposure experiments. Preliminary permeability measurements obtained from the core samples showed changes after they were exposed to CO 2 -saturated brine for one and six months. This observation suggests that mineral dissolution and mineral precipitation could occur in the host deposit altering its characteristics for CO 2 storage over time. The 3D images of the pore space clearly illustrate the degree of dissolution that occurred during exposure. It is noted that the dissolution that occurred during the six-month exposure could have enhanced the connectivity between voids. This may contribute the increase of permeability after the CO 2 /brine exposure. In addition, the primary minerals comprising the core are dolomite and gypsum. Both minerals could dissolve in the CO 2 /brine environment resulting in the increase of permeability after the six-month exposure.

58 GEOSCIENCES↗

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

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

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Controlling solute channel formation using magnetic fields

Solute channel formation introduces compositional and microstructural variations in a range of processes, from metallic alloy solidification, to salt fingers in ocean and water reservoir flows. Applying an external magnetic field interacts with thermoelectric currents at solid/liquid interfaces generating additional flow fields. This thermoelectric (TE) magnetohydrodynamic (TEMHD) effect can impact on solute channel formation, via a mechanism recently drawing increasing attention. To investigate this phenomenon, we combined in situ synchrotron X-ray imaging and Parallel-Cellular-Automata-Lattice-Boltzmann based numerical simulations to study the characteristics of flow and solute transport under TEMHD. Observations suggest the macroscopic TEMHD flow appearing ahead of the solidification front, coupled with the microscopic TEMHD flow arising within the mushy zone are the primary mechanisms controlling plume migration and channel bias. Two TE regimes were revealed, each with distinctive mechanisms that dominate the flow. Further, we show that grain orientation modifies solute flow through anisotropic permeability. These insights led to a proposed strategy for producing solute channel-free solidification using a time-modulated magnetic field.

36 MATERIALS SCIENCE↗

Nonlinear quantum behavior of ultrashort-pulse optical parametric oscillators

The quantum features of ultrashort-pulse optical parametric oscillators (OPOs) are investigated theoretically in the nonlinear regime near and above threshold. Viewing the pulsed OPO as a multimode open quantum system, we rigorously derive a general input-output model that features nonlinear coupling among many cavity (i.e., system) signal modes and a broadband single-pass (i.e., reservoir) pump field. Under appropriate assumptions, our model produces a Lindblad master equation with multimode nonlinear Lindblad operators describing two-photon dissipation and a multimode four-wave-mixing Hamiltonian describing a broadband, dispersive optical cascade, which we show is required to preserve causality. To simplify the multimode complexity of the model, we employ a supermode decomposition to perform numerical simulations in the regime where the pulsed supermodes experience strong single-photon nonlinearity. We find that the quantum nonlinear dynamics induces pump depletion as well as corrections to the below-threshold squeezing spectrum predicted by linearized models. We also observe the formation of non-Gaussian states with Wigner-function negativity and show that the multimode interactions with the pump, both dissipative and dispersive, can act as effective decoherence channels. Finally, we briefly discuss some experimental considerations for potentially observing such quantum nonlinear phenomena with ultrashort-pulse OPOs on nonlinear nanophotonic platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗

Evaluation and Optimization of Well Completion Options for the Utah FORGE Site

Orientation and completion for well pairs that have been subjected to multi-zonal stimulation play a critical role in the long-term performance of an Enhanced Geothermal Reservoir. Enhanced geothermal systems often rely on preferential flow along fractures between well injection and production locations. Modeling this preferential flow using discrete fracture networks (DNF) relies on stochastic realizations of the DFN based on geological sampling. Here we present the development of a stochastic optimization methodology to determine well completion options in a discrete fracture network based on using parallel subset simulation. Stochastic optimization will provide insight into regions where placements of the injection and production wells are optimal. An example optimization of well-pair location optimization based on a deterministic-stochastic DFN model representing FORGE follows a discussion of the theory.

15 GEOTHERMAL ENERGY↗

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R 2 ) and a low mean squared error. Specifically, the achieved R 2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. In conclusion, the present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.

58 GEOSCIENCES↗

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗