Stabilized formation temperatures were estimated at different depths in 40 wells from the Los Humeros geothermal field, Mexico, using the Horner and the spherical radial flow (SRF) methods. The results showed that the Horner method underestimates formation temperatures, while the SRF method gives temperatures that are closer to the true formation temperatures. This was supported by numerical simulation of a combined circulation and shut-in period in several wells, and results for well H-26 are presented. Numerical reproduction of logged temperature is more feasible if an initial temperature profile based on the SRF method is employed instead of using an initial temperature profile based on the Horner method.