We present a mathematical model for determining the temperature field around radioactive waste containers in very deep geological boreholes. The mathematical and physical ingredients of the model, underlying assumptions and approximations, and method of solution are discussed, together with results for two realistic disposal scenarios.