Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “borehole pressure”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bryan Mound Abandoned Cavern 3 Stability Analysis - 2021 Review

The U.S. Strategic Petroleum Reserve is moving towards employing an expanded enhanced monitoring program. In doing so it has become apparent that there is a need for a better project wide understanding of the current state of Bryan Mound abandoned Cavern 3 stability. Cavern 3 has been inaccessible since 1988 when it was plugged and abandoned and thus this comprehensive report is structured by focusing on 1) a summarization of what can be discerned from historical records prior to 1988 and 2) a presentation and discussion of our current understanding of Cavern 3 based solely on surface monitoring and geomechanical analyses. Historical literature state the cavern was deemed unsuitable for oil storage, as it could not be definitively determined if fluid pressure could be maintained in the borehole. Current surface monitoring indicates the largest surface subsidence rates are occurring above Cavern 3. The subsidence rates are linear with no evidence of acceleration. Cavern collapse could occur if there is insufficient pressure holding up the roof. Next steps are to implement a microseismic system that will lend to a better understanding of cavern stability, as well as provide an improved early warning system for loss of integrity.

02 PETROLEUM↗
DOE OSTI2020

Distributed Acoustic Sensing as a Distributed Hydraulic Sensor in Fractured Bedrock

Distributed acoustic sensing (DAS) was originally intended to measure oscillatory strain at frequencies of 1 Hz or more on a fiber optic cable. Recently, measurements at much lower frequencies have opened the possibility of using DAS as a dynamic strain sensor in boreholes. A fiber optic cable mechanically coupled to a geologic formation will strain in response to hydraulic stresses in pores and fractures. A DAS interrogator can measure dynamic strain in the borehole, which can be related to fluid pressure through the mechanical compliance properties of the formation. Because DAS makes distributed measurements, it is capable of both locating hydraulically active features and quantifying the fluid pressure in the formation. We present field experiments in which a fiber optic cable was mechanically coupled to two crystalline rock boreholes. The formation was stressed hydraulically at another well using alternating injection and pumping. The DAS instrument measured oscillating strain at the location of a fracture zone known to be hydraulically active. Rock displacements of less than 1 nm were measured. Laboratory experiments confirm that displacement is measured correctly. These results suggest that fiber optic cable embedded in geologic formations may be used to map hydraulic connections in three-dimensional fracture networks. A great advantage of this approach is that strain, an indirect measure of hydraulic stress, can be measured without beforehand knowledge of flowing fractures that intersect boreholes. The technology has obvious applications in water resources, geothermal energy, CO 2 sequestration, and remediation of groundwater in fractured bedrock.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Recent Developments at the Raft River Geothermal Field

The Raft River geothermal field, located in Cassia County in southwestern Idaho, is the site of a Department of Energy Enhanced Geothermal System project. U.S. Geothermal, Inc. currently produces about 11 MWe from Precambrian metamorphic rocks. These lie beneath ~5,000 ft of Quaternary and Tertiary volcaniclastic and volcanic deposits. Maximum temperatures range from 271°F to 300°F. Well RRG-9 ST1, the well targeted for stimulation is located approximately 1 mile south of the main bore field. The open hole section of the well, from 5,551 to 5,900 ft MD, consists of Precambrian Elba Quartzite, the stimulation target, granite and minor diabase. Prior to setting the casing acoustic, gamma ray, and density logs were run. After completing the well, a step rate/step down test was conducted. The maximum injection rate achieved was 18 bpm at a wellhead pressure of 1,150 psig. A borehole televiewer run in the open hole section showed evidence of more than eighty fractures. The majority of these fractures trend from N20⁰W to N20°E and dip from 40° to 60°W. Permeable fractures were encountered in the Elba Quartzite at 5,640-5,660 ft MD. Analysis of the injection test indicates that the minimum in-situ principal stress in this zone is 3,050-3,200 psi, corresponding to a fracture gradient of 0.59-0.62 psi/ft. A discrete fracture network model was developed using measured and inferred fracture orientations, distributions and dimensions. A three-phase stimulation program is proposed for RRG-9 ST-1. During the first two stages, water at 140°F, and later 40°F, will be injected to pre-condition and thermally fracture the reservoir. Here, the third stage will consist of a high rate, large volume conventional hydraulic stimulation.

Bradford, Jacob↗
DOE OSTI2025

Permeability scaling relationships of volcanic tuff from core to field scale measurements

A recent chemical explosive test in P-Tunnel at the Nevada National Security Site, Nevada, USA, was conducted to better understand how signals propagate from explosions in the subsurface. A primary signal of interest is the migration of gases that can be used to differentiate chemical from nuclear explosions. Gas migration is highly dependent on the rock permeability which is notoriously difficult to determine experimentally in the field due to a potentially large dependence on the scale over which measurements are made. Here, we present pre-explosion permeability estimates to characterize the geologic units surrounding the recent test. Permeability measurements were made at three scales of increasing size: core samples (≈2 cm), borehole packer system tests (≈1 m), and a pre-shot cavity pressurization test (> 10 m) across ten tuff units. Permeability estimates based on core measurements showed little difference from borehole packer tests. However, permeability in most rock units calibrated from cavity pressurization tests resulted in higher permeability estimates by up to two orders of magnitude. Here, we demonstrate that the scale of the measurement significantly impacts the characterization efforts of hydraulic properties in volcanic tuff, and that local-scale measurements (< 10 m scale) do not incorporate enough heterogeneity to accurately predict field-scale flow and mass transport.

Environmental sciences↗
DOE OSTI2020

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

UT-GOM2-2 Drilling Fluid (Technical Note)

It is proposed that a salt-saturated, water-based mud might improve borehole stability for UT-GOM2-2 relative to a fresh-water-based mud. However, the primary objective for UT-GOM2-2 is to safely and successfully acquire uncompromised hydrate-bearing sediment cores. This can only be done if conditions stay within the hydrate stability zone, and well away from the hydrate stability boundary. The addition of salt to the drilling mud, shifts the hydrate stability boundary closer to estimated conditions. This shift, by some estimates of pressure and temperature, where conditions now fall outside of the hydrate stability zone, may result in borehole enlargement and the release of free gas into the borehole. This shift by all estimates, even when coring conditions stay inside the hydrate stability zone, shrinks the window between the estimated conditions and the hydrate stability boundary. This shift is likely to compromise the hydrate-bearing cores while they are being recovered from the bottom-hole to the rig floor. In this report we present two examples: 1) Assuming temperature and pressure from in-situ estimates: Hydrate within the target reservoirs (Orange and Blue) at Walker Ridge Block 313 (WR 313) will be outside the hydrate stability zone in the presence of a drilling mud with 9.5 wt.% salinity (the salinity of the proposed 10.5 ppg salt-based mud) at in situ pressure and temperature. Thus, the hydrate will be unstable. A 10.5 ppg salt-based mud may result in dissociation of the hydrate into its components: water and gas. 10.5 ppg salt-based mud may enlarge the borehole, release free gas into the borehole, and compromise the cores. 2) Assuming temperatures equal to the measured LWD borehole temperatures at this location: Hydrate within the target reservoirs will be stable with a salinity of 9.5 wt.% (the salinity of the proposed 10.5 ppg salt-based mud), but possibly unstable (just at the methane hydrate stability boundary) for a salinity of 17.2 wt.% (the salinity of the proposed 13.5 ppg salt-based kill mud). if the borehole has the same temperature as recorded during previous LWD drilling at this location. A 13.5 ppg salt-based kill mud may destabilize the borehole. More importantly, a 10.5 ppg salt-based mud does not provide a large enough window between the estimated conditions and the hydrate stability boundary. Core temperatures and pressures fluctuate from the bottom-hole conditions as the core is brought up from the bottom-hole to the rig floor. These fluctuations are likely to cause the hydrate in the cores to touch or cross the stability boundary. Therefore, cores captured using a 10.5 ppg salt-based mud are likely to be compromised during core recovery to the rig floor.

03 NATURAL GAS↗
DOE OSTI2021

Time-reversed nonlinear acoustic downhole pore pressure measurements

Apparatus (10) and methods for measurement of pore pressure in rock formations through a metal borehole casing (32) after a well is cased and cemented, are described. Such measurements may be accomplished by using the Dynamic Acoustic Elasticity (DAE) method for characterizing nonlinear parameters by perturbing a selected rock formation region with a High Amplitude, Low Frequency (HALF) acoustic strain wave, and probing this region using a Low Amplitude, High Frequency (LAHF) acoustic wave (18), (22). Time reversal techniques (36) may be employed for focusing acoustic energy into the formation in the vicinity of the pipe or open hole. The change in wave speed of the probe pulses as the HALF induced strain wave oscillation propagates through the formation, as a function of the induced strain, may be used to determine the nonlinear elastic parameters α, β, δ, and A of the pore pressure, from which the pore pressure may be determined in the region of the HALF wave.

Goodman, Harvey E.↗
DOE OSTI2021

Thermal Stimulation and Injectivity Testing at Raft River, ID EGS Site

The Raft River geothermal field is the site of an innovative Department of Energy Enhanced Geothermal System (EGS) project to determine the viability of using combined thermal and hydraulic stimulation techniques to improve energy production. Well RRG-9 is currently undergoing a stimulation program using injectate from the US Geothermal Raft River Power Plant and cold water from a cooling tower make-up water well. The stimulation began on 13 June 2013 with injection from the power plant at a temperature of about 39 °C and a pressure of 275 psig. Next, two positive displacement plunger type pumps were used to increase the injection pressure and flow rate for about one month. The highest rate achieved was 258 gpm at a pressure of 741 psig. During this time, fluid from the cooler water well was injected for about 2 weeks at various pressures. Then, the pumps were removed and plant injection resumed on 25 September. Plant injection will continue until the spring of 2014, when a high pressure hydraulic stimulation will be conducted. A series of seismic monitoring stations deployed around the well are providing data on seismic events occurring at the site. Over the past year, 51 microseismic events have been recorded, all less than Magnitude 1. During injection, several diagnostic tests were conducted to gain a better understanding of the well and reservoir. A step-rate test was performed on 22 August to measure the in-situ stress and aid in modeling in-situ fractures. A tracer was injected into the well on 9 September. No tracer was detected in adjacent production wells after several months. A second borehole televiewer survey was conducted for comparison to pre-stimulation images. Here, a third borehole televiewer survey is planned after the high pressure stimulation. Injection test data is evaluated in real time. A modified Hall plot analysis indicates the effective permeability is increasing. The injectivity index supports the results of the modified Hall plot analysis. As the thermal stimulation has continued, the injectivity index has consistently followed an upward trend from 0.1 gpm/psi to 0.53 gpm/psi.

enhanced geothermal systems↗
DOE OSTI2025

Insights into the hydro-mechanical behavior of a decimeter-scale fracture using the mini-SIMFIP probe

ABSTRACT: Understanding hydro-mechanical couplings in fractured rocks is essential for predicting the rock mass response during high-pressure fluid injection, including the stimulation of enhanced geothermal systems. However, fluid-driven fracture dislocations are challenging to measure due to the need for local displacement data at high fluid pressures. In this study, fluid-driven displacement across a decimeter-scale laboratory fracture was investigated using the mini-SIMFIP (step rate injection method for fracture in-situ properties) probe, which is a smaller version of the SIMFIP tool (Guglielmi et al., 2014). The mini-SIMFIP probe is able to resolve the full 3D displacement vector of a fracture. The probe was installed in one of two boreholes across a decimeter-scale saw-cut granite fracture. Two pressure step injection tests were conducted under the same isotropic stress conditions. By varying injection between the boreholes, we estimated the aperture profile across the fracture. The results show a consistent pressure-dependent opening as long as steady-state flow was maintained. At a certain pressure step, steady-state conditions were no longer achievable and the fracture opened rapidly. The pressure-opening relationship diverged between the two tests and indicated a homogenization of the aperture profile and an increase in the non-linearity of the flow regime.

Osten, J↗
DOE OSTI2020

Temperature response and brine availability to heated boreholes in bedded salt

Abstract There is a growing need for disposal of high‐level nuclear waste. To reduce uncertainty associated with brine availability to repository excavations in salt formations, a collaboration between Sandia, Los Alamos, and Lawrence Berkeley National Laboratories is performing a series of borehole‐scale coupled process tests. Here, we report on the first round of the Brine Availability Test in Salt (BATS) project, a “shakedown” experiment called Phase 1s. Experimental testing included placing a resistive heater, a 260‐W radiative heater, and a 750‐W radiative heater within previously drilled horizontal boreholes at the Waste Isolation Pilot Plant (WIPP) while monitoring temperature and water inflow. The experiments successfully achieved the targeted temperature of 120 °C when using the 750‐W radiative heater. Simulations using FEHM (Finite Element Heat and Mass transfer code) and TOUGH‐FLAC (Transport Of Unsaturated Groundwater and Heat–Fast Lagrangian Analysis of Continua) were able to accurately predict the coupled thermo–hydro–mechanical–chemical response of salt, matching the observed temperature and brine production. Due to the extremely low permeability of salt, these systems take many years to reach steady state when perturbed by mining activities. Long‐term numerical simulations are used to develop the initial pressure and saturation conditions. The inclusion of a damaged rock zone with higher permeability around the borehole also affects the saturation and pressure distributions and plays an important role in dissipating the potential for thermal pressurization. Knowledge gained from this round of experimentation and simulation will be used to conduct the next BATS project experiment in new boreholes at WIPP.

58 GEOSCIENCES↗
DOE OSTI2021

Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation

At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

58 GEOSCIENCES↗
DOE OSTI2026

Quarterly Research Performance Progress Report

The goal of this project is to develop a high temperature, high-pressure well integrity evaluation device that can operate in enhanced geothermal system (EGS) boreholes without active cooling or substantial mitigation of borehole conditions. The product must provide consistent high resolution data regarding the efficacy of wellbore casing and cement for long-term operation. This project aims to deliver a ruggedized tool that can withstand a corrosive geothermal environment with pressure ranges up to 15,000 psi and temperatures up to 250°C for a 24-hour operating period.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
DOE OSTI2021

In-situ observation of pre-, co- and post-seismic shear slip at 1.5 km depth

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and transient afterslip. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. < p >The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment.

Schoenball, Martin↗
DOE OSTI2024

The EGS Collab project: Outcomes and lessons learned from hydraulic fracture stimulations in crystalline rock at 1.25 and 1.5 km depth

With the goal of better understanding stimulation in crystalline rock for improving enhanced geothermal systems (EGS), the EGS Collab Project performed a series of stimulations and flow tests at 1.25 and 1.5 km depths. The tests were performed in two well-instrumented testbeds in the Sanford Underground Research Facility in Lead, South Dakota, United States. The testbed for Experiment 1 at 1.5 km depth contained two open wells for injection and production and six instrumented monitoring wells surrounding the targeted stimulation zone. Four multi-step stimulation tests targeting hydraulic fracturing and nearly year-long ambient temperature and chilled water flow tests were performed in Experiment 1. The testbed for Experiments 2 and 3 was at 1.25 km depth and contained five open wells in an outwardly fanning five-spot pattern and two fans of well-instrumented monitoring wells surrounding the targeted stimulation zone. Experiment 2 targeted shear stimulation, and Experiment 3 targeted low-flow, high-flow, and oscillating pressure stimulation strategies. Hydraulic fracturing was successful in Experiments 1 and 3 in generating a connected system wherein injected water could be collected. However, the resulting flow was distributed dynamically, and not entirely collected at the anticipated production well. Thermal breakthrough was not observed in the production well, but that could have been masked by the Joule-Thomson effect. Shear stimulation in Experiment 2 did not occur - despite attempting to pressurize the fractures most likely to shear - because of the inability to inject water into a mostly-healed fracture, and the low shear-to-normal stress ratio. The EGS Collab experiments are described to provide a background for lessons learned on topics including induced seismicity, the correlation between seismicity and permeability, distributed and dynamic flow systems, thermoelastic and pressure effects, shear stimulation, local geology, thermal breakthrough, monitoring stimulation, grouting boreholes, modeling, and system management.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2021

Downhole nonlinear acoustics measurements in rock formations using dynamic acoustic elasticity and time reversal

Apparatus and methods for measurement of pore pressure in rock formations through an open, or cemented and/or cased, borehole are described. Such measurements are achieved using the Dynamic Acoustic Elasticity (DAE) method for characterizing nonlinear parameters by perturbing a selected rock formation volume with a High Amplitude, Low Frequency (HALF) acoustic strain wave, and probing this volume using a Low Amplitude, High Frequency (LAHF) acoustic wave. Time reversal techniques may be employed for focusing acoustic energy Into the formation in the vicinity of the pipe or open hole.

Goodman, Harvey Edwin↗
DOE OSTI2022

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗
DOE OSTI2024

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↗