Modeling of heat and solute transport in a fracture-matrix mine thermal energy storage system and energy storage performance evaluation

Repurposing groundwater-filled mine cavities for thermal energy storage has demonstrated promising potential to buffer the imbalance of energy supply and demand. Fractured formations are widespread in old mines due to previous excavation activities, which dominates the performance and efficiency of the mine thermal energy storage (MTES) system. Therefore, understanding the mechanisms of fluid flow, heat, and solute transport in fractured reservoirs in the MTES system is essential for its application and the assessment of environmental impact. In this study, fluid flow, heat, and solute transport in a 3D fracture-matrix MTES system are modeled simultaneously based on site-specific data in Freiberg, Germany. The simulations are conducted with a coupled Hydro-Thermo-Component process newly developed in the open-source software OpenGeoSys (OGS). To stabilize the simulation in the fracture-matrix hybrid system, a flux-corrected transport scheme is specifically implemented and applied in the model. Thus, heat and mass exchange between the embedded fractures and the matrix in the formation of MTES can be accurately simulated on the basis of the conforming mesh. In a hypothetical short-term scenario of a 15-day heating and cooling cycle of the MTES operation, the solute is transported faster than the heat in the fractured formation, indicating that disturbance of mineral compositions in the original mine water can affect a larger region than heat. The embedded fracture network in the surrounding formation is also found to increase the thermal energy storage capacity by 44%, while decreasing the recovery ratio by 14% compared to those of the unfractured formation. More energy will be transported and stored in the more fractured formation from heated mine water, but less energy can be recovered. The behavior of the MTES system strongly depends on local geology and hydraulic conditions and therefore needs to be evaluated in a site- and application-specific manner. This study provides preliminary insights into the performance and efficiency, as well as the environmental impact of the MTES operation in the short term. •A coupled HTC model with fracture effect stabilized by FCT for MTES is implemented.•The heat and mass transport characteristics in DFM formation are analyzed.•The effects of embedded DFN and groundwater on the MTES performance are studied.•DFN-embedded formation in MTES promotes heat to be stored but hinders recovery.