Compared to traditional tunnel construction, artificial freezing involves two distinct stages: freezing and thawing. However, this process can pose a risk to the surrounding environment if it is not possible to accurately analyze the frost heave and thawing settlement during the freezing process. The thermomechanical coupled mathematical model of formation frost heave was established by considering boundary conditions, such as formation temperature and convective heat transfer, and introducing the parameters of instantaneous volumetric strain and denaturation characteristic coefficient. In addition, the numerical analysis method for the whole construction process of the tunnel by the horizontal freezing method was established by programming the user subroutine that takes into account the characteristic coefficients of frozen soil relying on the secondary development technology of ABAQUS (version 2022). Then, the method was applied to the horizontal freezing engineering of a double-line tunnel, and the distribution laws of the freezing temperature field and frost heave displacement field were obtained and compared with the field measurement results. The numerical analysis method for determining the deformation law of the ground surface has been shown to be reliable through comparison with field measurements. This method can serve as a reference for designing effective ground surface heaving control schemes during the freezing construction period of tunnels in complex environments.