本文首先介绍了考虑相变的数理控制方程的推导过程,其中对相变的处理方法进行了详细的理论说明及公式推导,然后对多年冻土区的天然温度场进行了数值模拟,包括模型的建立、材料参数的确定、初始值和边界条件的确定、子步的确定等,最后将模拟的多年冻土区天然温度场与实测值进行了对比,验证了模拟的正确性,为今后多年冻土区温度场的相关研究奠定了基础。