The application of the Hydraulic Fracturing (HF) technology to exploit geothermal energy from Hot Dry Rocks (HDR) reservoirs is addressed. HF is achieved by extensively pumping geothermal fluids to already existing fractured HDR reservoirs of low permeability. High fluid pressures are expected to drive cracks to evolve and connect. The newly created burgeoning hydraulic conduits should supposedly enhance the permeability of the existing HDR reservoirs. The flow rate/pressure values at which geothermal fluids should be pumped, as well as the pumping schedule to initiate HF, depend primarily on the existing geostatic conditions (geostatic stresses, initial HDR pressure and temperature) as well as on HDR fracture properties (initial mean fracture length, mean fracture aperture, density and orientational distribution of fractures). While these components, in addition to their effects on borehole stability, are scrutinized in this research, focus is on the evolution during circulation processes of the fracture properties. A fracturing model that is capable of tracking fracture evolution in all possible spatial orientations is used to obtain the time course of the anisotropic permeability tensor. This evolving property is integrated into a domestic finite element code which is developed to solve thermo–poroelastic BVPs: emphasis is laid on the efficiency of the doublet flow technique where a fluid gains thermal energy by circulating through the HDR reservoir from the injection well to the production well. The spurious oscillations in the hyperbolic solutions of the approximated finite element approach that are commensal with the phenomenon of forced heat convection are healed/mitigated through several stabilization approaches.Ce travail concerne l'utilisation de la technique de Fracturation Hydraulique (FH) pour exploiter l'énergie géothermique des réservoirs profonds de roches sèches chaudes (HDR). La fracturation hydraulique est réalisée par injection de fluides géothermiques dans des réservoirs partiellement fracturés de faible perméabilité. Les fluides à haute pression sont destinés à faire évoluer les fissures et leur connectivité. Les valeurs de débit/pression auxquelles les fluides géothermiques doivent être pompés, ainsi que le calendrier de pompage pour initier la fracturation hydraulique, dépendent principalement des conditions géostatiques existantes (contraintes géostatiques, pression fluide et température initiales de l'HDR) ainsi que des propriétés des fissures de l'HDR (longueur, épaisseur, densité et distribution directionnelle initiales moyennes de fissures). Tous ces éléments, en sus de leurs effets sur la stabilité des forages, sont analysés dans cette recherche. Des modèles de fracturation, qui sont capables de suivre l'évolution des fissures dans toutes les orientations spatiales possibles, sont utilisés pour obtenir le tenseur anisotrope de perméabilité. Ces modèles sont intégrés dans un code domestique d'éléments finis qui est développé pour résoudre des problèmes aux limites thermo-poroélastiques. Pour supprimer/diminuer les oscillations qui accompagnent les solutions paraboliques et/ou hyperboliques lors de la convection forcée, plusieurs techniques de stabilisation ont dû être implémentées
