This work extends the hydro-mechanical phase-field fracture model to non-isothermal conditions with micromechanics based poroelasticity, which degrades Biot's coefficient not only with the phase-field variable (damage) but also with the energy decomposition scheme. Furthermore, we propose a new approach to update porosity solely determined by the strain change rather than damage evolution as in the existing models. As such, these poroelastic behaviors of Biot's coefficient and the porosity dictate Biot's modulus and the thermal expansion coefficient. For numerical implementation, we employ an isotropic diffusion method to stabilize the advection-dominated heat flux and adapt the fixed stress split method to account for the thermal stress. We verify our model against a series of analytical solutions such as Terzaghi's consolidation, thermal consolidation, and the plane strain hydraulic fracture propagation, known as the KGD fracture. Finally, numerical experiments demonstrate the effectiveness of the stabilization method and intricate thermo-hydro-mechanical interactions during hydraulic fracturing with and without a pre-existing weak interface.