Based on the ＇double-temperature＇ concept proposed by Aifantis and Beskos [1], a thermoelastic analytical model is presented for the study of fractured geothermal reservoirs. The fluid flow is not explicitly formula fed as a result of analogous behavior between the heat and the fluid conduction. Therefore, the thermal flux transported in the porous media may be considered as an average quantity split between the fluid and the rock matrix. This equivalent thermal transport model may be envisioned as a phenomenological simplification of reality, because the heat in general is transported by the moving fluid in a geothermal reservoir. Using a displacement potential, the effect of reversible volumetric expansion on the change of thermal flux has been partially decoupled in the temperature field. The full coupling is restored in the solid equilibrium by placing the strain due to thermal stress into the deformation assemblage. The thermoelastic solutions are obtained in radial flow coordinates, imitating an extraction scenario from a centrally located well. Parametric studies determine the key factors such as thermal conductivity ratio and fracture spacing, critical to the reservoir performance. The double-porosity phenomenon of heat flow, similar to the fluid flow in a fractured porous medium, has been identified.