Innovative numerical procedure for simulating borehole heat exchangers operation and interpreting thermal response test through MODFLOW-USG code

In recent years, among renewable energies, the geothermal resource exploitation shows a constant growth; specifically, in countries engaged in CO2 emissions reduction and dependent on energy from abroad, the lowtemperature geothermal energy (geo-exchange) for air conditioning of buildings represents a cost-effective and green solution. In closed-loop systems borehole heat exchangers (BHE) are coupled with ground-source heat pumps (GSHP) constituting the key component of the heating ventilation air-conditioning (HVAC) system. Therefore, the design of the BHE and the correct interpretation of in situ Thermal Response Tests (TRT) are essential to supply the building energy demand. To support the design, Modflow-USG Connected Linear Network (CLN) and Drain Return Flow (DRT) packages are adapted and improved to reproduce the operation of one or more BHE in aquifers and to analyze the TRT. The improvements are compared with a previously developed numerical model and two different analytical solutions (infinite line source and moving line source) by imposing a constant heat rate injection into the aquifer. The results show good agreement between the new approach and previous ones (discrepancy lower than 2% for models with highly refined grid), but the new approach is much more accurate and expeditious in both implementation and execution, also allowing for an easy numerical simulation of multiple BHE.