Numerical Modeling of a Soil‐Borehole Thermal Energy Storage System

Core Ideas Borehole thermal energy storage is studied with a 3D transient fluid flow and heat transfer model. BTES heat extraction efficiency increases with decreasing soil thermal conductivity. BTES efficiency decreases with convective heat losses associated with high soil permeability. Borehole thermal energy storage (BTES) in soils combined with solar thermal energy harvesting is a renewable energy system for the heating of buildings. The first community‐scale BTES system in North America was installed in 2007 at the Drake Landing Solar Community (DLSC) in Okotoks, AB, Canada, and has since supplied >90% of the thermal energy for heating 52 homes. A challenge facing BTES system technology is the relatively low efficiency of heat extraction. To better understand the fluid flow and heat transport processes in soils and to improve BTES efficiency of heat extraction for future applications, a three‐dimensional transient coupled fluid flow and heat transfer model was established using TOUGH2. Measured time‐dependent injection temperatures and fluid circulation rates at DLSC were used as model inputs. The simulations were calibrated using measured soil temperature time series. The simulated and measured temperatures agreed well with a subsurface having an intrinsic permeability of 1.5 × 10 −14 m 2 , thermal conductivity of 2.0 W m −1 °C −1 , and a volumetric heat capacity of 2.3 MJ m −3 °C −1 . The calibrated model served as the basis for a sensitivity analysis of soil thermal and hydrological parameters on BTES system heat extraction efficiency. Sensitivity analysis results suggest that: (i) BTES heat extraction efficiency increases with decreasing soil thermal conductivity; (ii) BTES efficiency decreases with background groundwater flow; (iii) BTES heat extraction efficiency decreases with convective heat losses associated with high soil permeability values; and (iv) unsaturated soils show higher overall heat extraction efficiency due to convection onset at higher intrinsic permeability values.