[3] J Noorishad,TSANG C-F,WITHERSPOON P-A.Coupled thermal-hydraulic-mechanical phenomena in saturated fractured porous rocks:Numerical approach[J].Journal of Geophysical Research:Solid Earth,1984,89(B12)：10365-10373.

[17] YUE Zu-run,TAI Bo-wen,SUN Tie-cheng.Analysis of temperature field characteristics based on subgrade site measurements of Harbin-Qiqihar High-speed Railway in a deep seasonal frozen soil region[J].Sciences in Cold and Arid Regions,2015,7(5):7.

[18] BOWEN Tai,LIU Jian-kun,YUE Zu-run,et al.Effect of sunny-shady slopes and strike on thermal regime of subgrade along a high-speed railway in cold regions,China[J].Engineering Geology,2018,232：182-191.

[19] BOWEN Tai,LIU Jian-kun,WANG Teng-fei,et al.Numerical modelling of anti-frost heave measures of high-speed railway subgrade in cold regions[J].Cold Regions Science and Technology,2017,141：28-35.

[24] ZENG Yu-chao,HE Bin,TANG Lian-sheng,et al.Numerical simulation of temperature field and pressure field of the fracture system at Zhangzhou geothermal field[J].Environmental Earth Sciences,2020,79(3)：1-13.

[25] G-S Taylor,LUTHIN J-N.A model for coupled heat and moisture transfer during soil freezing[J].Revue Canadienne De Géotechnique,1978,15(4):548-555.

[27] M-Th, van GENUCHTEN.A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1[J].Soil Science Society of America Journal,1980,44(5)：892-898.

[28] ZHANG Ming-li,WEN Zhi,XUE Ke,et al.A coupled model for liquid water,water vapor and heat transport of saturated-unsaturated soil in cold regions:model formulation and verification[J].Environmental Earth Sciences,2016,75(8)：1-19.