This article presents a solution to the geothermal problem of transient heat production from hot dry rocks using a horizontal well. Dimensionless forms of the governing equations are derived, including conduction in the rock, convection between the wellbore rock and the fluid, and advection and conduction in the fluid along the well. Ten model material and geometric parameters are reduced to three dimensionless parameters: α=4LDSt is the ratio of the rate of heat storage in the fluid to the rate the heat convection to the fluid from the rock, β=1Pe is the ratio of the rate of conductive versus advective heat transfer in the fluid, and γ=2LDBi is the ratio of heat convection to the rock from the fluid to the rate of heat depletion in the rock. An axisymmetric finite element method (FEM) program is developed and yields the solution for the temperature of the rock mass and fluid over time. For physically meaningful inputs, analysis indicates that the effect of β is negligible, that combinations of very large α and very low γ (and vice versa) do not occur, and that all other things being equal, increasing α or decreasing γ leads to higher fluid outlet temperatures. System response at different times and dimensionless temperature are captured in contour plots for values of α and γ spanning practical injection rates, well geometries and fluid and rock properties. Power law response surface models are fitted using model outputs at discrete intervals and provide a means to accurately and rapidly compute the temperature-time histories of practical geothermal systems.