Analysis of hydrothermal coupling response characteristics of frozen soil subgrade of high-speed railway in cold region
-
摘要:
冻土水热耦合问题的控制方程具有强耦合特性,使得相应的数值计算存在一定挑战,进而影响其在工程实践中的应用。根据能量守恒、质量守恒原理及土体冻结曲线,给出了考虑相变效应的冻土水热耦合理论模型,而后数学推导得到解耦的理论模型方程以便优化数值求解。基于COMSOL平台二次开发实现了冻土水热耦合过程的数值建模。使用兰州—乌鲁木齐(兰新)客运专线路基的实测数据进行数值计算的验证,并在拟合的地表边界条件下开展了该冻土路基水热耦合的数值分析。分析表明:(1)不同深度观测点对应的温度和含水率数值解与实测值具有较好的一致性,从而验证了所解耦的冻土水热耦合理论模型的可靠性。(2)土层对温度和含水率周期性变化时的幅值均有“削峰”作用,且不同深度观测点的温度和含水率正弦变化曲线均有一定的相位滞后现象;其中,温度幅值削峰和相位滞后是热传导过程的能量耗散引起,而含水率曲线的类似现象则可能是冰水相变改变土层渗透性的缘故。(3)近地表附近温度等值线较密,而远地表土层中温度等值线较疏,表明路基表层更易受外界温度波动影响,夏季时温度自上而下逐渐降低,而冬季时温度自上而下逐渐升高。(4)路基断面中含水率随着深度的增加而增大,约在含水泥粗粒土材料和填料的界面附近达到峰值,体现了材料界面对水分迁移的影响效应,而后随深度增加含水率逐渐降低。研究成果可为寒区路基等工程的建造提供一定的理论支撑。
Abstract:Coupled heat and water transfer in frozen soil presents certain challenges for numerical computation due to the strong coupling nature of its control equations, thereby affecting their application in engineering practice. Based on the principles of energy conservation, mass conservation, and the soil freezing curve, a theoretical model for the coupled heat and water transfer in frozen soil, considering phase change effects, was proposed. Subsequently, decoupled control equations were derived through mathematical deduction to optimize numerical solutions. Numerical modeling of the coupled water and heat transfer process in frozen soil was implemented through secondary development on the COMSOL platform. The measured data from the subgrade of the Lanxin Passenger Dedicated Line were used for numerical validation, and numerical analysis was conducted under fitted surface boundary conditions. The analysis indicates the numerical solutions of temperature and moisture content at different depth exhibit good agreement with the measured values, thereby validating the reliability of the decoupled theoretical model for coupled heat and water transfer in frozen soil. The soil layer has a “peak damping” effect on the amplitude of periodic variations in temperature and moisture content, and the sine wave curves of temperature and moisture content at different depth present certain phase lag phenomena. Among these, temperature amplitude attenuation and phase lag are caused by energy dissipation in the heat conduction process, while similar phenomena in moisture content curves may be due to phase changes between ice and water altering soil permeability. The temperature contour lines near the surface are denser, while those in the deeper soil layers are sparser. The surface layer of the roadbed is more susceptible to external temperature fluctuations. During summer, the temperature gradually decreases from top to bottom, whereas in winter, the temperature gradually increases from top to bottom. The moisture content in the subgrade section increases with depth, reaching a peak near the interface between the cemented coarse-grained soil material and the fill material, highlighting the influence of material interfaces on moisture migration, and then gradually decreases with increasing depth. The research findings provide technical support for the construction of engineering projects such as roadbeds in cold regions.
-
-
表 1 主要物性参数
Table 1. Physics parameters
物性参数 取值 含水泥粗粒土 填料 路基填土 水 冰 密度/(kg·m−3) 2 300 2 060 1 900 1 000 918 体积热容/(kJ·m−3·K−1) 920 860 1 000 4 180 1 874 导热系数/(W·m−1·K−1) 1.51 1.41 1.00 0.58 2.22 渗透系数/(m·s−1) 1×10−7 1×10−4 1×10−6 — — 表 2 计算参数
Table 2. Calculation parameters
参数 a0 l m θs θr Tf 数值 2 0.5 0.5 0.48 0.05 − -
[1] 郑拓. 我国高速铁路与经济发展研究[J]. 铁道学报,2020,42(7):34 − 41. [ZHENG Tuo. Research on development of China’s high-speed railway and economy[J]. Journal of the China Railway Society,2020,42(7):34 − 41. (in Chinese with English abstract)]
ZHENG Tuo. Research on development of China’s high-speed railway and economy[J]. Journal of the China Railway Society, 2020, 42(7): 34 − 41. (in Chinese with English abstract)
[2] PHILIP J R,DE VRIES D A. Moisture movement in porous materials under temperature gradients[J]. Transactions,American Geophysical Union,1957,38(2):222 − 232. doi: 10.1029/TR038i002p00222
[3] HARLAN R L. Analysis of coupled heat-fluid transport in partially frozen soil[J]. Water Resources Research,1973,9(5):1314 − 1323. doi: 10.1029/WR009i005p01314
[4] HANSSON K,SIMUNEK J,MIZOGUCHI M,et al. Water flow and heat transport in frozen soil:Numerical solution and freeze-thaw applications[J]. Vadose Zone Journal,2004,3(2):693 − 704.
[5] 毛雪松. 多年冻土地区路基水热力场耦合效应研究[D]. 西安:长安大学,2004. [MAO Xuesong. Study on coupling effect of hydrothermal field of subgrade in permafrost region[D]. Xi’an:Chang’an University,2004. (in Chinese with English abstract)]
MAO Xuesong. Study on coupling effect of hydrothermal field of subgrade in permafrost region[D]. Xi’an: Chang’an University, 2004. (in Chinese with English abstract)
[6] 周家作,李东庆,房建宏,等. 开放系统下饱和正冻土热质迁移的数值分析[J]. 冰川冻土,2011,33(4):791 − 795. [ZHOU Jiazuo,LI Dongqing,FANG Jianhong,et al. Numerical analysis of heat and mass transfers in saturated freezing soil in an open system[J]. Journal of Glaciology and Geocryology,2011,33(4):791 − 795. (in Chinese with English abstract)]
ZHOU Jiazuo, LI Dongqing, FANG Jianhong, et al. Numerical analysis of heat and mass transfers in saturated freezing soil in an open system[J]. Journal of Glaciology and Geocryology, 2011, 33(4): 791 − 795. (in Chinese with English abstract)
[7] 常启昕,孙自永,马瑞,等. 冻土区地下水流过程及其与地表水转化关系研究进展[J]. 水利水电科技进展,2016,36(5):87 − 94. [CHANG Qixin,SUN Ziyong,MA Rui,et al. A review of groundwater flow and its interaction with surface water in permafrost region[J]. Advances in Science and Technology of Water Resources,2016,36(5):87 − 94. (in Chinese with English abstract)] doi: 10.3880/j.issn.1006-7647.2016.05.016
CHANG Qixin, SUN Ziyong, MA Rui, et al. A review of groundwater flow and its interaction with surface water in permafrost region[J]. Advances in Science and Technology of Water Resources, 2016, 36(5): 87 − 94. (in Chinese with English abstract) doi: 10.3880/j.issn.1006-7647.2016.05.016
[8] 梁昆淼. 数学物理方法[M]. 4版. 北京:高等教育出版社,2010. [LIANG Kunmiao. Methods of mathematical physics[M]. 4th ed. Beijing:Higher Education Press,2010. (in Chinese)]
LIANG Kunmiao. Methods of mathematical physics[M]. 4th ed. Beijing: Higher Education Press, 2010. (in Chinese)
[9] ASMAR N H. 偏微分方程教程[M]. 陈祖墀,宣本金,译. 2版. 北京:机械工业出版社,2006. [ASMAR N H. Partial differential equations with Fourier series and boundary value problems[M]. CHEN Zuchi,XUAN Benjin,trans. 2nd ed. Beijng:China Machine Press,2006. (in Chinese)]
ASMAR N H. Partial differential equations with Fourier series and boundary value problems[M]. CHEN Zuchi, XUAN Benjin, trans. 2nd ed. Beijng: China Machine Press, 2006. (in Chinese)
[10] 尚松浩,雷志栋,杨诗秀. 冻结条件下土壤水热耦合迁移数值模拟的改进[J]. 清华大学学报(自然科学版),1997,37(8):62 − 64. [SHANG Songhao,LEI Zhidong,YANG Shixiu. Numerical simulation improvement of coupled moisture and heat transfer during soil freezing[J]. Journal of Tsinghua University (Science and Technology),1997,37(8):62 − 64. (in Chinese with English abstract)]
SHANG Songhao, LEI Zhidong, YANG Shixiu. Numerical simulation improvement of coupled moisture and heat transfer during soil freezing[J]. Journal of Tsinghua University (Science and Technology), 1997, 37(8): 62 − 64. (in Chinese with English abstract)
[11] 白青波,李旭,田亚护,等. 冻土水热耦合方程及数值模拟研究[J]. 岩土工程学报,2015,37(增刊2):131 − 136. [BAI Qingbo,LI Xu,TIAN Yahu,et al. Equations and numerical simulation for coupled water and heat transfer in frozen soil[J]. Chinese Journal of Geotechnical Engineering,2015,37(Sup2):131 − 136. (in Chinese with English abstract)]
BAI Qingbo, LI Xu, TIAN Yahu, et al. Equations and numerical simulation for coupled water and heat transfer in frozen soil[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(Sup2): 131 − 136. (in Chinese with English abstract)
[12] GHIAS M S,THERRIEN R,MOLSON J,et al. Numerical simulations of shallow groundwater flow and heat transport in continuous permafrost setting under impact of climate warming[J]. Canadian Geotechnical Journal,2019,56(3):436 − 448. doi: 10.1139/cgj-2017-0182
[13] 毛卫南,刘建坤. 冻土斜坡路基温度水分分布特性分析[J]. 铁道学报,2015,37(6):88 − 94. [MAO Weinan,LIU Jiankun. Analysis of temperature and water content distribution of slope embankment in permafrost regions[J]. Journal of the China Railway Society,2015,37(6):88 − 94. (in Chinese with English abstract)] doi: 10.3969/j.issn.1001-8360.2015.06.013
MAO Weinan, LIU Jiankun. Analysis of temperature and water content distribution of slope embankment in permafrost regions[J]. Journal of the China Railway Society, 2015, 37(6): 88 − 94. (in Chinese with English abstract) doi: 10.3969/j.issn.1001-8360.2015.06.013
[14] 雷志栋,尚松浩,杨诗秀,等. 地下水浅埋条件下越冬期土壤水热迁移的数值模拟[J]. 冰川冻土,1998,20(1):51 − 54. [LEI Zhidong,SHANG Songhao,YANG Shixiu,et al. Numerical simulation on simultaneous soil moisture and heat transfer under shallow ground water table in winter[J]. Journal of Glaciology and Geocryology,1998,20(1):51 − 54. (in Chinese with English abstract)]
LEI Zhidong, SHANG Songhao, YANG Shixiu, et al. Numerical simulation on simultaneous soil moisture and heat transfer under shallow ground water table in winter[J]. Journal of Glaciology and Geocryology, 1998, 20(1): 51 − 54. (in Chinese with English abstract)
[15] BAI Bing,ZHOU Rui,YANG Guangchang,et al. The constitutive behavior and dissociation effect of hydrate-bearing sediment within a granular thermodynamic framework[J]. Ocean Engineering,2023,268:113408. doi: 10.1016/j.oceaneng.2022.113408
[16] 杨世铭,陶文铨. 传热学[M]. 4版. 北京:高等教育出版社,2006. [YANG Shiming,TAO Wenquan. Heat transfer[M]. 4th ed. Beijing:Higher Education Press,2006. (in Chinese)]
YANG Shiming, TAO Wenquan. Heat transfer[M]. 4th ed. Beijing: Higher Education Press, 2006. (in Chinese)
[17] 徐学祖,王家澄,张立新. 冻土物理学[M]. 北京:科学出版社,2001:74 − 97. [XU Xuezu,WANG Jiacheng,ZHANG Lixin. Physics of frozen soils[M]. Beijing:Science Press,2001:74 − 97. (in Chinese)]
XU Xuezu, WANG Jiacheng, ZHANG Lixin. Physics of frozen soils[M]. Beijing: Science Press, 2001: 74 − 97. (in Chinese)
[18] 卢宁,WILLIAM J L. 非饱和土力学[M]. 韦昌富,侯龙,简文星,译. 北京:高等教育出版社,2012. [LU Ning,WILLIAM J L. Unsaturated soil mechanics[M]. WEI Changfu,HOU Long,JIAN Wenxing,trans. Beijing:Higher Education Press,2012. (in Chinese)]
LU Ning, WILLIAM J L. Unsaturated soil mechanics[M]. WEI Changfu, HOU Long, JIAN Wenxing, trans. Beijing: Higher Education Press, 2012. (in Chinese)
[19] ZHAO Junlin,ZHANG Pei,YANG Xiao,et al. On the uniaxial compression strength of frozen gravelly soils[J]. Cold Regions Science and Technology,2020,171:102965. doi: 10.1016/j.coldregions.2019.102965
[20] 孟生勇,江兴元,杨义,等. 降雨诱发堆积体滑坡水土响应与稳定性时空演化试验研究[J]. 水文地质工程地质,2023,50(1):104 − 112. [MENG Shengyong,JIANG Xingyuan,YANG Yi,et al. An experimental study of spatial-temporal evolution of water-soil response and stability of a rainfall-induced accumulation landslide[J]. Hydrogeology & Engineering Geology,2023,50(1):104 − 112. (in Chinese with English abstract)]
MENG Shengyong, JIANG Xingyuan, YANG Yi, et al. An experimental study of spatial-temporal evolution of water-soil response and stability of a rainfall-induced accumulation landslide[J]. Hydrogeology & Engineering Geology, 2023, 50(1): 104 − 112. (in Chinese with English abstract)
[21] 刘青灵,简文彬,许旭堂,等. 基于可靠度方法的全基质吸力段土-水特征模型研究[J]. 水文地质工程地质,2022,49(1):92 − 100. [LIU Qingling,JIAN Wenbin,XU Xutang,et al. A study of the soil-water reliability model in the whole matric suction range[J]. Hydrogeology & Engineering Geology,2022,49(1):92 − 100. (in Chinese with English abstract)]
LIU Qingling, JIAN Wenbin, XU Xutang, et al. A study of the soil-water reliability model in the whole matric suction range[J]. Hydrogeology & Engineering Geology, 2022, 49(1): 92 − 100. (in Chinese with English abstract)
[22] 李强,李同录,李华,等. 毛细水作用下非饱和土压缩过程的微观非连续变形数值分析[J]. 水文地质工程地质,2022,49(4):135 − 143. [LI Qiang,LI Tonglu,LI Hua,et al. Numerical analysis of evolution of the unsaturated soil micro-structure with capillary action during compression[J]. Hydrogeology & Engineering Geology,2022,49(4):135 − 143. (in Chinese with English abstract)]
LI Qiang, LI Tonglu, LI Hua, et al. Numerical analysis of evolution of the unsaturated soil micro-structure with capillary action during compression[J]. Hydrogeology & Engineering Geology, 2022, 49(4): 135 − 143. (in Chinese with English abstract)
[23] TAYLOR G S,LUTHIN J N. A model for coupled heat and moisture transfer during soil freezing[J]. Canadian Geotechnical Journal,1978,15(4):548 − 555. doi: 10.1139/t78-058
[24] BURT T P,WILLIAMS P J. Hydraulic conductivity in frozen soils[J]. Earth Surface Processes,1976,1(4):349 − 360. doi: 10.1002/esp.3290010404
[25] 徐学祖,邓友生. 冻土中水分迁移的实验研究[M]. 北京:科学出版社,1991. [XU Xuezu,DENG Yousheng. Experimental study on water migration in frozen soil[M]. Beijing:Science Press,1991. (in Chinese)]
XU Xuezu, DENG Yousheng. Experimental study on water migration in frozen soil[M]. Beijing: Science Press, 1991. (in Chinese)
[26] 邵珠杰. 高海拔季节冻土区高速铁路路基水-热-冻胀变形特征研究 ——以兰新客专民乐段为例[J]. 冰川冻土,2018,40(3):588 − 597. [SHAO Zhujie. The characteristics of high-speed railway subgrade’s temperature,moisture and frost heave deformation in high altitude and seasonal frozen region:Taking the Minle section of Lanzhou-Xinjiang passenger railway line as an example[J]. Journal of Glaciology and Geocryology,2018,40(3):588 − 597. (in Chinese with English abstract)]
SHAO Zhujie. The characteristics of high-speed railway subgrade’s temperature, moisture and frost heave deformation in high altitude and seasonal frozen region: Taking the Minle section of Lanzhou-Xinjiang passenger railway line as an example[J]. Journal of Glaciology and Geocryology, 2018, 40(3): 588 − 597. (in Chinese with English abstract)
[27] BAI Bing,YANG Guangchang,LI Tao,et al. A thermodynamic constitutive model with temperature effect based on particle rearrangement for geomaterials[J]. Mechanics of Materials,2019,139:103180. doi: 10.1016/j.mechmat.2019.103180
[28] 曹元兵,盛煜,吴吉春,等. 上边界条件对多年冻土地温场数值模拟结果的影响分析[J]. 冰川冻土,2014,36(4):802 − 810. [CAO Yuanbing,SHENG Yu,WU Jichun,et al. Influence of upper boundary conditions on simulated ground temperature field in permafrost regions[J]. Journal of Glaciology and Geocryology,2014,36(4):802 − 810. (in Chinese with English abstract)]
CAO Yuanbing, SHENG Yu, WU Jichun, et al. Influence of upper boundary conditions on simulated ground temperature field in permafrost regions[J]. Journal of Glaciology and Geocryology, 2014, 36(4): 802 − 810. (in Chinese with English abstract)
-