The development of a novel method for modelling fluid flow and heat transfer in naturally fractured geothermal reservoirs represents a significant advancement in geothermal energy research. This Study presents a hybrid approach, which combines discrete fracture and single continuum techniques, to effectively capture the complex interactions between fluid flow and heat transfer in geothermal fractured reservoirs. In addition, the incorporation of the local thermal nonequilibrium method for simulating heat transmission accounts for the disparities in temperature between the rock matrix and the fluid, providing a more realistic representation of heat transfer processes. The study also presents a fully coupled thermo-poro-elastic framework that integrates fluid flow and heat transfer to comprehensively evaluate reservoir responses to injection/production scenarios. This coupled approach allows for the prediction of changes in reservoir properties, such as permeability and porosity, under varying fluid pressure and temperature conditions. The application of the proposed model to evaluate a geothermal reservoir’s long-term response to injection/production scenarios provides valuable insights into the reservoir’s behaviour and potential energy production capacity. The sensitivity analysis further enhances the model’s utility by identifying the key reservoir parameters that significantly influence the thermal depletion of the reservoir. Overall, this novel modelling approach holds promise for improving the understanding and management of naturally fractured geothermal reservoirs, contributing to the optimization of geothermal energy extraction strategies.