A particle-tracking formulation of advective-diffusive heat transport in deformable fracture networks

Modeling heat transfer in complex heterogeneous fractured system is key for geothermal energy applications. Discrete fracture network (DFN) modeling is the classical framework to reproduce the advective part of the transport, which is determined by the fracture connectivity and heterogeneity. This approach in general sacrifices the representation of the rock matrix, disregarding both its diffusive heat exchange with the fractures and the effects of its thermo-mechanical deformation on the fracture aperture. Here we propose a new semi-analytic formulation that can be implemented in a DFN simulator with particle tracking approach. The contribution of the rock matrix in terms of diffusive heat exchange and thermal contraction/expansion is analytically evaluated, which respectively impact the advective heat transfer and the fracture aperture variation. The method is proved to be accurate and robust. Results from simulations of cold fluid injection show that rock contraction affects the transmissivity, which accelerates the advective transport resulting in a faster recovery of cold fluid at the outlet. The methodology enables investigating the reservoir behavior and optimizing the geothermal performance while keeping the computational effort within reasonable values. This allows exploring the uncertainty in cases when the in situ characterization is poor, which is the spirit of the DFN modeling.

Automated summary

Modeling heat transfer in complex heterogeneous fractured systems is crucial for geothermal energy applications. A new semi-analytic formulation has been proposed to evaluate the contribution of the rock matrix in diffusive heat exchange and thermal contraction/expansion, impacting the advective heat transfer and fracture aperture variation. This method has been proven accurate and robust, allowing for studying reservoir behavior and optimizing geothermal performance effectively.