Engineering topics

Smith, J. Torquil

Publications and source records attributed to Smith, J. Torquil.

7 records

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↗

DOE OSTI2021

Adjoint-Based Inversion of Geodetic Data for Sources of Deformation and Strain

An adjoint-based formulation leads to a particularly efficient approach for inverting geodetic measurements for the source of the deformation. Specifically, the quantities necessary to iteratively improve the fit to the observations can be computed with just three forward calculations, one to obtain the current residuals, another to solve the adjoint problem, and a third to compute the step length. An inversion algorithm utilizing the adjoint-based gradient is applied to a set of Interferometric Synthetic Aperture Radar (InSAR) data gathered between 2016 and 2018 over the Tulare Basin in California's Central Valley. Because the measured deformation is due to groundwater withdrawal, a penalty function is included in the inversion to avoid placing aquifer volume change in locations that are far from any documented wells. The solution of the inverse problem provides estimates of aquifer compaction that provide a match to the observed range changes while honoring the well data. The solution indicates an average aquifer volume loss of 2.17 km 3 /year over the two year period from January 2016 to January 2018, encompassing one drought year (2016) and one wet year (2017). Finally, this magnitude of lost volume is compatible with the 3.1 km 3 /year decrease in water volume for the entire Central Valley, estimated from GRACE satellite gravity data.

58 GEOSCIENCES↗

DOE OSTI2021

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗

DOE OSTI2015

Newberry Well 55-29 Stimulation Data 2014

The Newberry Volcano EGS Demonstration in central Oregon, a 5 year project begun in 2010, tests recent technological advances designed to reduce the cost of power generated by EGS in a hot, dry well (NWG 55-29) drilled in 2008. First, the stimulation pumps used were designed to run for weeks and deliver large volumes of water at moderate well-head pressure. Second, to stimulate multiple zones, AltaRock developed thermo-degradable zonal isolation materials (TZIMs) to seal off fractures in a geothermal well to stimulate secondary and tertiary fracture zones. The TZIMs degrade within weeks, resulting in an optimized injection/ production profile of the entire well. Third, the project followed a project-specific Induced Seismicity Mitigation Plan (ISMP) to evaluate, monitor for, and mitigate felt induced seismicity. An initial stimulation was conducted in 2012 and continued for 7 weeks, with over 41,000 m3 of water injected. Further analysis indicated a shallow casing leak and an unstable formation in the open hole. The well was repaired with a shallow casing tieback and perforated liner in the open hole and re-stimulated in 2014. The second stimulation started September 23rd, 2014 and continued for 3 weeks with over 9,500 m3 of water injected. The well was treated with several batches of newly tested TZIM diverter materials and a newly designed Diverter Injection Vessel Assembly (DIVA), which was the main modification to the original injection system design used in 2012. A second round of stimulation that included two perforation shots and additional batches of TZIM was conducted on November 11th, 2014 for 9 days with an additional 4,000 m3 of water injected. The stimulations resulted in a 3-4 fold increase in injectivity, and PTS data indicates partial blocking and creation of flow zones near the bottom of the well. This submission includes all of the files and reports associated with the stimulation, pressure testing, and monitoring included in the scope of the project.

15 GEOTHERMAL ENERGY↗

DOE OSTI2013

dwholeprinceresp.f90

Program dwholeprinceresp.f90 computes responses at receivers of a detection/discrimination system in the time domain to an induced magnetic dipole polarizability object specified in the frequency domain, in a whole space of a given conductivity. Magnetic fields of a set of loop transmitters are propagated through the conductive medium (e.g. seawater) to an object whose magnetic dipole moment polarizabilities are specified in the frequency domain, the resultant induced magnetic moments calculated, the magnetic fields of the induced moments the propagated to a set of user specified receivers (approximated as point dipole receivers), the response at the receivers is convolved with the transmitter waveform, and the result inverse Fourier transformed to get the time domain response at the receivers. The time domain response of the receivers to the primary (direct) transmitter magnetic field is computed separately in a similar manner. Transmitter loops may be approximated with any number of line segments, and any number of receivers may be specified. Square pulse and half sine pulse transmitter wave forms of arbitrary duration are modeled. The waveform may be uni-modal or bi-model (alternating sign). The polariability object may be positioned and oriented arbitrarily with respect to the system.

Smith, J. Torquil↗

DOE OSTI2012

Borehole EM Monitoring at Aquistore: Final Report to the Carbon Capture Project (CCP)

Geologic carbon sequestration (GCS) is a technology whose goal is to prevent atmospheric release of greenhouse gases via injection of carbon dioxide (CO 2 ) into an underground reservoir for long term storage. GCS is typically part of a program of carbon capture and storage (CCS) that captures CO 2 from point sources such as power plants, transports the CO 2 to a storage site, and operates an injection facility. One recent CCS pilot project is the Aquistore CO 2 sequestration project, near Estevan, Saskatchewan, Canada. The Aquistore project is managed by the Petroleum Technology Research Centre (PTRC) and will be one of the first integrated CCS projects storing CO 2 in a deep saline aquifer from a coal fired power plant (PTRC, 2011). Aquistore is expected to store 500,000 tons of CO 2 during its lifetime (Ministry of Environment, 2012). Assuring the long-term, safe storage of CO 2 requires the development of effective monitoring strategies. As part of the geophysical monitoring effort at Aquistore, there were initial plans for deployment of borehole electrodes for electrical or electromagnetic measurements to monitor CO 2 within the reservoir. The injected CO 2 displaces saline brine in the reservoir, and because CO 2 has a high resistivity compared to brine, the overall resistivity of the formation increases and can be monitored by measuring electric or magnetic fields. Previous work by Lawrence Berkeley National Laboratory (LBNL) had indicated that borehole-to-surface electromagnetic monitoring, using an electric dipole source near the bottom of a well penetrating the reservoir, could detect the resistivity change induced by GCS. To assess the potential application of electromagnetic monitoring at Aquistore, Lawrence Berkeley National Laboratory and Multi-Phase Technologies collaborated on a two-part study including (1) numerical forward modeling of a time-lapse, controlled-source electromagnetic (CSEM) monitoring survey, and (2) an initial engineering study of instrumentation and proposed design for a borehole electric dipole source and electrode sensors.

54 ENVIRONMENTAL SCIENCES↗