The thermo-hydrological (TH) coupling model constitutes the foundational framework for investigating the temperature distribution of surrounding rock in cold region tunnels. In this study, a fully coupled TH model is proposed that takes into account multiple physical phenomena during the freezing process of surrounding rock. Firstly, the model was established based on thermodynamics, seepage theory, and ice–water phase change theory, which accounted for unfrozen water, latent heat of phase change, ice impedance, and convective heat transfer. The model was successfully verified by comparing its results to field data. Next, the sensitivity of surrounding rock temperature to environmental, thermodynamic, seepage, and coupling parameters in the fully coupled TH model was systematically studied using a numerical analysis method. The results show that the annual temperature amplitude and thermal conductivity represent the main factors affecting the surrounding rock temperature at a radial depth of 0 m, while the initial temperature and porosity are the key factors at a radial depth of 5 m. Permeability has the least influence on the surrounding rock temperature, but the temperature field will experience sudden changes if its value exceeds its value exceeds 1 × 10−12 m2. Finally, using the proposed numerical model, the thickness of insulation layer was simulated, and the degree of influence of the parameters on the thickness of insulation layer was analyzed. This study reveals that the annual temperature amplitude has the greatest influence on the calculation of insulation layer thickness, with its normalized sensitivity factor being approximately 50%. These findings not only expand the methodology for exploring the laws of TH coupling but also provide a theoretical foundation for improving the parameter calibration efficiency and calculation accuracy of the fully coupled TH model, and they have significant reference value.
Wu et al. (Fri,) studied this question.