The stochastic finite element equations for random temperature are obtained using the first-order per-turbation technique taking into account the random thermal properties and boundary condition, based on heat transfer variational principle. The local average method for 2-D is used to discretize random fields. Then, the random temperature fields of embankment in cold regions are investigated on condi-tion that the thermal properties and boundary condition are taken as random fields, respectively, by using the program, which is written by the methods. The expected value of temperature field and the standard deviation of the temperature field of embankment in cold regions are obtained and analyzed.