Investigating the Effect of Fracture Properties on Peclet Number of Enhanced Geothermal Systems

Understanding heat transfer mechanisms in subsurface environments is crucial for deep geothermal energy exploitation. Despite multiple studies on heat transfer in porous media, the combination of different heat transfer mechanisms with fracture networks leads to uncertainty in the temperature distribution in geothermal sites. More specifically, heat convection and conduction are two major mechanisms responsible for heat transfer in porous media. Conduction occurs through rock matrix and is more dominant than convection in intact rocks due to their low porosity. However, fracture networks in rocks increase fluid transfer in certain directions, making convection more important in heat transfer. At a certain point, it becomes difficult to identify which mechanism is more prevalent. Besides this difficulty, achieving a balance between these two key mechanisms and determining the optimal Peclet number is very crucial for the geothermal energy industry to extract adequate volume of water at the desired temperature, which is essential for smooth operation of geothermal systems. In this study, a two dimensional coupled thermo-hydraulic model in COMSOL is developed to simulate heat transfer mechanisms in enhanced geothermal systems on field scale. Characteristic time variables based on hydraulic and thermal diffusivities are defined to monitor heat distribution through the domain caused by each mechanism. Finally, a sensitivity analysis is performed to identify how the temperature is affected by rock (thermal and hydraulic) parameters and fracture patterns. The simulation results indicate that the Peclet number is highly dependent on fracture network.