Case Studies of Geothermal System Response to Perturbations in Groundwater Flow and Thermal Regimes.

Global demands for energy-efficient heating and cooling systems coupled with rising commitments toward net zero emissions is resulting in wide deployment of shallow geothermal systems, typically installed to a depth of 100 to 200 m, and in the continued growth of the global ground source heat pump (GSHP) market. Ground coupled heat pump (GCHP) systems take up to 85% of the global GSHP market. With increasing deployment of GCHP systems in urban areas coping with limited regulations, there is growing potential and risk for these systems to impact the subsurface thermal regime and to interact with each other or with nearby heat-sensitive subsurface infrastructure. In this paper, we present three numerical modeling case studies, from the UK and Canada, which examine GCHP systems' response to perturbation of the wider hydrogeological and thermal regimes. The studies demonstrate how GCHP systems can be impacted by external influences and perturbations arising from subsurface activities that change the thermal and hydraulic regimes in the area surrounding these systems. Additional subsurface heat loads near existing schemes are found to have varied impacts on system efficiency with reduction ranging from <1% to 8%, while changes in groundwater flow rates (due to a nearby groundwater abstraction) reduced the effective thermal conductivity at the study site by 13%. The findings support the argument in favor of regulation of GCHP systems or, to a minimum, their registration with records of locations and approximate heat pump capacity-even though these systems do not abstract/inject groundwater.

Borehole Heat Exchangers-Addressing the Application Gap with Groundwater Science.

Hydrogeologists and mechanical engineers approach the design of geoexchange systems, and the associated borehole heat exchanger (BHE) fields, in different ways, each focusing on their knowledge areas. While these differences have created a strong research base, with well-published innovations and designs that collectively allow for sustainable systems, industry has not embraced these recent advancements. Despite abundant research demonstrating how complex shallow groundwater flow and temperature conditions can influence BHE design and operation, the low-temperature geothermal industry remains largely fixed on simple analytical codes and assumed uniform ground conditions. Geoexchange system inefficiencies become masked via reduced heat pump performance and increased electricity consumption. Similarly, interactions between BHE fields and infrastructure in urban areas are slow to manifest and are often unrealized due to a lack of field temperature data. While regulations that include hydrogeological factors have been developed in some jurisdictions, they are largely voluntary or rudimentary and can be unapplied in industry. Addressing this application gap may be unreasonable as designing and installing thermally efficient geoexchange systems can put them out of the cost envelope of competing heating and cooling systems. Perhaps for hydrogeologists, the way forward lies in the use of BHE's to facilitate contaminated sites remediation, an area we are familiar with, and one that allows for innovative technologies to reduce cost envelopes. Following that path, hydrogeologists can help improve system efficiencies while fully considering the dynamic nature of advective and thermal transport by groundwater.

Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow.

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.

Impact of Groundwater Flow and Energy Load on Multiple Borehole Heat Exchangers.

The effect of array configuration, that is, number, layout, and spacing, on the performance of multiple borehole heat exchangers (BHEs) is generally known under the assumption of fully conductive transport. The effect of groundwater flow on BHE performance is also well established, but most commonly for single BHEs. In multiple-BHE systems the effect of groundwater advection can be more complicated due to the induced thermal interference between the boreholes. To ascertain the influence of groundwater flow and borehole arrangement, this study investigates single- and multi-BHE systems of various configurations. Moreover, the influence of energy load balance is also examined. The results from corresponding cases with and without groundwater flow as well as balanced and unbalanced energy loads are cross-compared. The groundwater flux value, 10(-7) m/s, is chosen based on the findings of previous studies on groundwater flow interaction with BHEs and thermal response tests. It is observed that multi-BHE systems with balanced loads are less sensitive to array configuration attributes and groundwater flow, in the long-term. Conversely, multi-BHE systems with unbalanced loads are influenced by borehole array configuration as well as groundwater flow; these effects become more pronounced with time, unlike when the load is balanced. Groundwater flow has more influence on stabilizing loop temperatures, compared to array characteristics. Although borehole thermal energy storage (BTES) systems have a balanced energy load function, preliminary investigation on their efficiency shows a negative impact by groundwater which is due to their dependency on high temperature gradients between the boreholes and surroundings.

Geothermal waste heat utilization from in situ thermal bitumen recovery operations.

In situ thermal methods for bitumen extraction introduce a tremendous amount of energy into the reservoirs raising ambient temperatures of 13 °C to as high as 200 °C at the steam chamber edge and 50 °C along the reservoir edge. In essence these operations have unintentionally acted as underground thermal energy storage systems which can be recovered after completion of bitumen extraction activities. Groundwater flow and heat transport models of the Cold Lake, Alberta, reservoir, coupled with a borehole heat exchanger (BHE) model, allowed for investigating the use of closed-loop geothermal systems for energy recovery. Three types of BHEs (single U-tube, double U-tube, coaxial) were tested and analyzed by comparing outlet temperatures and corresponding heat extraction rates. Initial one year continuous operation simulations show that the double U-tube configuration had the best performance producing an average temperature difference of 5.7 °C, and an average heat extraction of 41 W/m. Given the top of the reservoir is at a depth of 400 m, polyethylene piping provided for larger extraction gains over more thermally conductive steel piping. Thirty year operation simulations illustrate that allowing 6 month cyclic recovery periods only increases the loop temperature gain by a factor of 1.2 over continuous operation. Due to the wide spacing of existing boreholes and reservoir depth, only a small fraction of the energy is efficiently recovered. Drilling additional boreholes between existing wells would increase energy extraction. In areas with shallower bitumen deposits such as the Athabasca region, i.e. 65 to 115 m deep, BHE efficiencies should be larger.

Radar determination of the spatial structure of hydraulic conductivity.

Spatial variability of hydraulic conductivity exerts a predominant control on the flow of fluid through porous media. Heterogeneities influence advective pathways, hydrodynamic dispersion, and density-dependent dispersion; they are, therefore, a key concern for studies of ground water resource development, contaminant transport, and reservoir engineering. Ground-penetrating radar contributes to the remote, geophysical characterization of the macroscale variability of natural porous media. On a controlled excavation of a glacial-fluvial sand and gravel deposit in the Fanshawe Delta area (Ontario, Canada), the hydraulic conductivity field of a 45 x 3 m vertical exposure was characterized using constant-head permeameter measurements performed on undisturbed horizontal sediment cores. Ground-penetrating radar data were collected along the excavation face in the form of both reflection and common midpoint surveys. Comparison of geostatistical analyses of the permeameter measurements and the radar data suggests thatthe horizontal correlation structure of radar stack velocity can be used to directly infer the horizontal correlation structure of hydraulic conductivity. The averaging nature of the common midpoint survey is manifest in the vertical correlation structure of stack velocity, making it less useful. Radar reflection data do not exhibit a spatial structure similar to that of hydraulic conductivity possibly because reflections are a result of material property contrasts rather than the material properties themselves.

Errors with small volume elastic seepage meter bags.

The use of small volume elastic collection bags (condoms) has become popular in seepage meter studies in recent years, despite minimal field or laboratory validation of their use and, specifically, the impact of their elasticity on seepage measurements. A laboratory study was initiated after field results using small elastic collection bags produced seepage data that did not correlate with hydrometric data. The laboratory data demonstrate that condoms undergo significant mechanical relaxation during seepage measurement times typically observed in field settings. Unlike conventional nonelastic collection bags, which mechanically relax over several minutes, the condoms suffered from a slow mechanical relaxation or equilibration. Over nine hours, condoms gained 43 mL of water, approximately 50% of maximum workable volume (between mechanical relaxation effect and elastic limit), under stagnant flow conditions. This long-term equilibration invalidates simple subtraction of equilibration volumes from collection volumes as a correction technique. Previously published studies using flexible small-volume elastic measurement bags (condoms) have not reported a mechanical relaxation effect. Overall, because the condom's small workable volume and inherent variability, we would not recommend any small-volume elastic measurement bags for quantitative seepage measurements.