An alternative method of the numerical solution of the heat equation is presented. A three-dimensional modelling approach is used in this paper rather than one or two dimensional models in order to account for the lateral heterogeneous thermal conductivity coefficients or the heat sources or sinks commonly encountered in geothermal reservoir engineering studies.
