The paper is concerned with mathematical modeling of heat and mass transfer in subsurface geothermal systems. Heat exchanger is driven as a vertical channel in an impermeable rock massive. The temperature of rock is driven in accord with a usual thermal conductance equation on a cylinder coordinates basis. The boundary conditions between the rock and heat exchanger are given by the equality of a thermal flow.