冰川冻土, 2022, 44(2): 458-469 doi: 10.7522/j.issn.1000-0240.2022.0051

寒区工程与灾害

青藏工程走廊冻融土热扩散系数特性与预测模型研究

刘志云,1, 钟振涛1, 崔福庆,1, 陈建兵2, 彭惠2

1.长安大学 地质工程与测绘学院,陕西 西安 710054

2.中交第一公路勘察设计研究院有限公司 高寒高海拔地区 道路工程安全与健康国家重点实验室,陕西 西安 710075

Study on characteristic and prediction model of permafrost thermal diffusivity of Qinghai-Tibet engineering corridor

LIU Zhiyun,1, ZHONG Zhentao1, CUI Fuqing,1, CHEN Jianbing2, PENG Hui2

1.College of Geology Engineering and Geomatics,Chang’an University,Xi’an 710054,China

2.State Key Laboratory of Road Engineering Safety and Health in Cold and High-Altitude Regions,CCCC First Highway Consultants Co. Ltd. ,Xi’an 710075,China

通讯作者: 崔福庆,讲师,主要从事冻土传热传质研究. E-mail: cfq731@chd.edu.cn

收稿日期: 2021-05-31   修回日期: 2021-08-29  

基金资助: 国家自然科学基金项目.  51574037.  41502292
中国交建科技研发项目.  2020-ZJKJ-PTJS04.  2020-ZJKJ-QNCX09.  2020-ZJKJ-PTJS12

Received: 2021-05-31   Revised: 2021-08-29  

作者简介 About authors

刘志云,副教授,主要从事冻土工程热灾害防治研究.E-mail:dcdgx33@chd.edu.cn , E-mail:dcdgx33@chd.edu.cn

摘要

热扩散系数是多年冻土对外界热扰动敏感程度的重要影响参数之一,也是寒区工程设计与建设的关键基础数据。基于瞬态平面热源法导热系数测试结果和质量加权法计算获取的比热容理论值,计算获得青藏工程走廊西大滩—唐古拉山沿线典型类别土样热扩散系数,分析对比了走廊带内冻融土热扩散系数的分布特征和参数影响规律,提出了基于经验拟合公式法和RBF神经网络方法的冻融土热扩散系数预测模型,并比较了不同预测模型的预测效果。研究结果表明:青藏工程走廊带内土的热扩散系数与粒径整体呈正相关性,融土热扩散系数按黏性土、粉土、全风化岩类、砂土及碎石土依次增大,冻土热扩散系数按黏性土、全风化岩类、粉土、碎石土及砂土依次增大;热扩散系数与容重及天然含水率相关性随土类及冻融状态差异明显,冻、融土热扩散系数呈显著正线性关系;以融土热扩散系数为拟合参数的冻土热扩散系数三元预测模型的预测精度明显高于二元经验公式;RBF神经网络模型在冻、融土热扩散系数预测中均具有最优的预测精度,为最佳预测模型。

关键词: 青藏工程走廊 ; 多年冻土 ; 热扩散系数 ; 瞬态平面热源法 ; 预测模型

Abstract

Thermal diffusivity is a prime influencing factor of permafrost thermal response sensitivity to external heat disturbances, and it is also the key basic data for engineering design and construction in cold regions. In present work, thermal diffusivity of typical soil samples along the Xidatan-Tanggula Mountain of the Qinghai-Tibet engineering corridor were calculated using the thermal conductivity test results of the transient plane heat source method and the theoretical value of specific heat calculated by mass weighting method. Then, the distribution characteristics and the parameter influence law of thermal diffusivity of frozen and unfrozen soil in the corridor zone had been analyzed and compared. Finally, the prediction models of thermal diffusivity of frozen and unfrozen soil based on empirical fitting formula method and RBF neural network method had been developed and compared. The research results showed that: (1) The thermal diffusivity of the soil in the Qinghai-Tibet engineering corridor is positively correlated with the particle size. The increasing order of thermal diffusivity value of unfrozen soil is cohesive soil, silt, fully weathered rock, sandy soil and gravel soil, and order is cohesive soil, fully weathered rock, silt, gravel soil and sandy soil for frozen soil. (2) The correlation among thermal diffusion coefficient, volume weight and natural moisture content varies with soil types and freezing/thawing states, and the thermal diffusivity of frozen and unfrozen soil has a significant positive linear relationship. (3) Using thermal diffusivity of unfrozen soil as fitting parameter, the prediction accuracy of ternary fitting prediction model is significantly higher than that of binary fitting model. (4) The RBF neural network model has the greatest prediction accuracy for both frozen and unfrozen soils, which is the best soil thermal diffusivity prediction model.

Keywords: Qinghai-Tibet engineering corridor ; permafrost ; thermal diffusivity ; transient plane source method ; prediction model

PDF (6168KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘志云, 钟振涛, 崔福庆, 陈建兵, 彭惠. 青藏工程走廊冻融土热扩散系数特性与预测模型研究[J]. 冰川冻土, 2022, 44(2): 458-469 doi:10.7522/j.issn.1000-0240.2022.0051

LIU Zhiyun, ZHONG Zhentao, CUI Fuqing, CHEN Jianbing, PENG Hui. Study on characteristic and prediction model of permafrost thermal diffusivity of Qinghai-Tibet engineering corridor[J]. Journal of Glaciology and Geocryology, 2022, 44(2): 458-469 doi:10.7522/j.issn.1000-0240.2022.0051

0 引言

青藏工程走廊作为最重要的陆路进出藏通道,穿越了约550 km多年冻土区,特殊的工程地质条件引发了冻胀、融沉、开裂、拥包、塌陷及不均匀变形等病害1-2。在未来的青藏高速公路工程建设中,高速公路黑色沥青路面特有的强烈吸热、聚热效应将给路基下伏冻土带来更大的热冲击3,诱发更为严重的工程病害4-5。如何保障这一特殊工程场景下公路工程的安全性和稳定性,是寒区冻土工程研究重要课题之一6。热扩散系数是表征土体内温度扩散速率的物理量,反映了传热过程中导热系数和比热容两大基础参数协同作用的效果7。同时,冻土的热扩散系数一定程度上也是冻土本身对外部热量输入响应敏感性的重要标尺,是描述多年冻土特征的关键参数之一8。因此,揭示多年冻土热扩散系数的变化规律、提出适用的预测模型对于未来青藏高速的设计施工和病害防治具有重要意义。

热扩散系数主要获取方法有理论计算和试验测试,其中理论计算是基于土壤为常热扩散系数的半无界介质假设,以一维热传导方程为理论基础获得计算公式,主要包括振幅法、相位法、对数法、反正切法、数值法和谐波法等9-10。许多研究者针对不同计算方法进行了对比分析,Horton等11从计算所需数据条件和计算结果方面对比评价了六种热扩散系数计算方法,发现具有显式方程的振幅、相位、反正切和对数方法所需地温数据少,但计算结果偏差大,而以隐式求解的数值法和谐波法通常以大量实测地温为计算基础,结果更为可靠。缪育聪等12根据现场观测地温数据计算了地表浅层5~20 cm土壤层的热扩散系数,对比发现充分利用实测地温信息的谐波法为最优估算方法。韩炳宏等13对青海南部高寒草地土壤的热扩散系数进行了计算,研究表明热传导对流法计算结果及其拟合效果最好,同时发现除干土层外随土壤深度加深,热扩散系数逐渐递减。同时,不少学者根据土壤的传热特征对计算方法进行了修正,Gao等14-16综合考虑了土壤热传输过程中的热传导和多孔热对流方式,给出了一维热传导对流方程,采用谐波法和拉普拉斯变换方法推导了热传导对流方程的解析解及热扩散系数的计算公式,并基于实测地温数据验证了改进方法的准确性和优越性。原黎明等17基于改进的热传导对流方法对青藏高原中部活动层表层下5~20 cm深度的土壤热扩散系数进行了估算,结果发现5~10 cm处土壤热扩散系数显著大于10~20 cm深度,且融化季对应热扩散系数显著高于冻结季。章永辉等18基于青藏高原理塘区野外观测地温数据,研究发现利用耦合热传导-对流法计算所得热扩散系数模拟地温时的准确性最高。在试验方面,甄作林等19利用瞬态平面热源法测试了原状和重塑砂土的热扩散系数,分析发现热扩散系数与干密度成正比关系,而含水率低时热扩散系数与含水率为正相关性,达到一定含水率时热扩散系数趋于稳定或与其呈负相关性。

关于热扩散系数影响因素方面,研究主要集中于土壤质地、内部结构、含水率、温度及深度等。王可里等20通过对青藏高原那曲地区非均质土壤热扩散系数计算发现,热扩散系数分布具有明显的深度和季节性变化特征,且冷暖季采用不同热扩散系数可获得较好地温模拟效果。周亚等21根据青藏高原观测站长期监测所得0.8 m和3.2 m深度土壤温度,给出了青藏高原不同地区深层土壤热扩散系数的年际和季节性变化特征。刘经星等22利用数值差分求解方法对不同质地条件土壤的热扩散系数进行了估算,结果发现土壤质地轻且土熵条件好对应热扩散系数小。Roxy23、安可栋等24根据观测数据研究了土壤水分和热扩散系数之间关系,研究表明热扩散系数随土壤水分增加先增大后减小。邸佳颖等25通过室内试验对比研究了原状土和装填土的热特性差异,结果表明装填土在处于中等含水率时热扩散系数较原状土有所增大,而接近饱和时热扩散系数趋于一致。马欣等26研究发现土壤热扩散系数与深度为非线性正相关性,且降水增多会降低热扩散系数。在计算模型方面,Arkhangelskaya等27以土壤质地、容重和有机碳百分含量作为影响变量,建立了不同土类任意含水率状态下的热扩散系数回归模型。

但应看到的是,目前冻土热扩散系数的研究报道还相对较少,缺乏针对青藏工程走廊带典型类别土样热扩散系数的大规模测试与预测模型研究。为了进一步揭示青藏工程走廊带沿线多年冻土的热物理性质,本文基于青藏公路改扩建(格尔木—拉萨段)工程地质勘察项目912组土样导热系数测试结果和质量加权法计算获取的比热容理论值,计算获得西大滩—唐古拉山沿线典型类别土样热扩散系数,分析了冻融土热扩散系数的分布特征和参数影响规律,提出了基于经验拟合公式法和RBF神经网络方法的冻融土热扩散系数预测模型,并对不同预测模型的预测效果进行了对比分析,以期为青藏高速公路的热工设计提供数据参考。

1 热扩散系数计算方法

热扩散系数定义为导热系数与体积比热的比值,是反映物体增温快慢的综合物性参数,计算公式如下:

α=λC

式中:α为土壤热扩散系数;λ为土壤导热系数,W·m-1·K-1C为体积比热,kJ·m-3·K-1

本文青藏工程走廊重塑冻融土导热系数采用Hot Disk TPS1500热常数分析仪以瞬态平面热源法进行测试,其中测试区段北起西大滩(即走廊带多年冻土北界)南至唐古拉山共437 km,对应青藏公路里程桩号为K2870~K3307。钻孔取样深度自地表向下最深达40 m,共计获取土样912组,其中测试融土样633组、冻土样858组。

对钻探获取的土样进行土性统计,将土样数超过15组的列为典型土样,按照土的工程分类,可将走廊带内所取土样分为黏性土、粉土、全风化岩类、砂土和碎石土共计5大类12个亚类,各土类统计如图1所示。由图可知,走廊带内主要以黏性土为主,砂土和碎石土次之,且二者数量相近。

图1

图1   青藏工程走廊导热系数测试区段及土性分布

Fig.1   Thermal conductivity test section and soil property distribution of Qinghai-Tibet engineering corridor


导热系数测试主要操作包括:(1)试样制备:土样烘干筛分、增湿至天然含水率、静置闷样及加压制样;(2)测试预处理:平整试样表面和预先冻结冻土试样;(3)物性测试:确定测试时间及功率、重复测试和数据误差对比。其中试样制备依据《土工试验方法标准》(GB/T50123—2019)规范内容进行,具体测试系统如图2所示。

图2

图2   导热系数测试系统

Fig.2   Thermal conductivity test system


体积比热按照文献[28]所给公式进行计算,该研究认为土比热是多相组分的质量加权平均结果,冻融土比热因固态冰的存在具有明显差异,具体计算式如下:

Cu=Csu+WCw1+Wρu
Cf=Csf+W-WuCi+WuCw1+Wρf

式中:CuCf分别为融土和冻土的体积热容量,kJ·m-3·K-1CsuCsfCwCi分别为融土骨架、冻土骨架、水及冰的质量比热,kJ·kg-1·K-1WWu分别为总含水量和未冻水含量;ρuρf分别为融土和冻土的天然密度,kg·m-3。其中不同土类冻融土骨架及水分质量比热取值见表1所示。

表 1   不同土类土骨架及水质量比热取值

Table 1  Mass specific heat of soil skeleton and water for different soil types

状态不同土类比热/(kJ·kg-1·K-1水比热/(kJ·kg-1·K-1
黏性土、粉质黏性土粉土、全风化千枚岩、 全风化泥岩、全风化板岩粉砂、中细砂、砾砂圆砾、角砾、碎石土
融化0.840.840.840.794.18
冻结0.770.750.730.712.09

新窗口打开| 下载CSV


2 冻融土热扩散系数分布特征及影响因素分析

2.1 冻融土热扩散系数分布特征

对五类土热扩散系数进行统计,得到其概率分布如图3所示。由图3可知,融土热扩散系数主要分布区间整体依黏性土、粉土、全风化岩类、砂土和碎石土偏右分布,表明融土热扩散系数主要分布范围值依次增大。可以看出,土的成分对融土的热扩散系数有显著影响,粒径越小,土的持水性越强,则结合水的比例越高,降低了水分在矿物骨架之间的联系作用,从而导致其热扩散系数降低。为排除由取样随机性带来的统计误差,取累积概率分布区间20%~80%热扩散系数为对比值,统计得黏性土、粉土、全风化岩类、砂土和碎石土融土热扩散系数分布区间为0.441×10-6~0.595×10-6、0.485×10-6~0.694×10-6、0.486×10-6~0.677×10-6、0.589×10-6~0.829×10-6和0.653×10-6~0.858×10-6 m2·s-1,均值分别为0.507、0.565、0.586、0.699和0.759×10-6 m2·s-1。可以看到,走廊带内各土类融土热扩散系数与粒径呈正相关性,分布基本排序为黏性土、粉土、全风化岩类、砂土及碎石土依次增大。

图3

图3   土类热扩散系数概率分布

Fig.3   Probability distribution of thermal diffusivity of soils


同理,图3中各土类冻土热扩散系数主要分布区间整体按黏性土、全风化岩类、粉土、碎石土和砂土依次增大,其累积概率分布区间20%~80%对应的热扩散系数区间为0.742×10-6~1.057×10-6、0.826×10-6~1.134×10-6、0.810×10-6~1.213×10-6、0.996×10-6~1.374×10-6和0.998×10-6~1.410×10-6 m2·s-1,均值分别为0.882×10-6、0.976×10-6、0.985×10-6、1.172×10-6和1.194×10-6 m2·s-1。对比发现,各土类冻融土热扩散系数分布规律不尽相同,可能与水分冻结引发土体结构变化有关,总体来说,走廊带内土类冻土热扩散系数分布基本排序为黏性土、全风化岩类、粉土、碎石土及砂土依次增大。

对比各土类冻土与融土热扩散系数发现,冻土热扩散系数均显著大于融土,其中黏性土、粉土、全风化岩类、砂土和碎石土均值相对增幅分别为73.96%、74.34%、66.55%、70.82%和54.41%,这是因为冻结状态下土体内水变为冰后具有更强的导热能力和较小的比热。同时,还可发现粒径较小土类整体具有较大相对增幅,可能与土体冻结过程土骨架结构连接性增强相关,粒径较小土类表现出更大增幅的导温性。

2.2 土体天然含水率、干密度对热扩散系数的影响

对走廊带内冻融土热扩散系数与天然含水率、干密度之间的偏相关性进行统计分析,结果如表2所示。由表可知,黏性土、粉土和碎石土冻融土热扩散系数与干密度呈显著的正相关关系,这是因为含水率不变的情况下,其质量比热也不变,而导热系数随干密度变大而增大,故随着干密度变大其冻融土热扩散系数相应增大。但也存在全风化岩类冻融土及砂土冻土热扩散系数与干密度为负相关性情况,由图4可知,在自然环境状况下不同干密度条件对应热扩散系数分布随机且离散,并不具有较强的正相关性。

表2   热扩散系数与影响因素偏相关分析

Table 2  Partial correlation analysis between thermal diffusivity and influencing factors

土性变量相关参数含水率干密度融土热扩散系数冻土热扩散系数
黏性土含水率相关性1.000-0.629-0.187-0.005
显著性0.0000.0400.934
干密度相关性-0.6291.0000.2070.118
显著性0.0000.0010.034
融土热扩散系数相关性-0.1870.2071.0000.803
显著性0.040.0010.000
粉土含水率相关性1.000-0.363-0.4750.002
显著性0.0080.0000.989
干密度相关性-0.3631.0000.0570.188
显著性0.0080.6890.118
融土热扩散系数相关性-0.4750.0571.0000.843
显著性0.0000.6890.000
全风化岩类含水率相关性1.000-0.842-0.259-0.210
显著性0.0000.0270.044
干密度相关性-0.8421.000-0.062-0.280
显著性0.0000.6030.007
融土热扩散系数相关性-0.259-0.0621.0000.785
显著性0.0270.6030.000
砂土含水率相关性1.000-0.567-0.1610.063
显著性0.0000.0650.368
干密度相关性-0.5671.0000.200-0.020
显著性0.0000.0210.778
融土热扩散系数相关性-0.1610.2001.0000.898
显著性0.0650.0210.000
碎石土含水率相关性1.000-0.610-0.0990.073
显著性0.0000.2620.355
干密度相关性-0.6101.0000.1690.149
显著性0.0000.0540.057
融土热扩散系数相关性-0.0990.1691.0000.821
显著性0.2620.0540.000

新窗口打开| 下载CSV


图4

图4   全风化岩类融土热扩散系数与干密度及天然含水率分布关系

Fig.4   Relationship between thermal diffusivity and dry density and natural water content distribution of unfrozen soil of fully weathered rocks


同时,可发现各土类冻融土热扩散系数与含水率间关系差异显著,大部分土类热扩散系数与含水率呈负相关性,但粉土、砂土及碎石土冻土热扩散系数与含水率呈正相关性,同样如图3所示,这可能与天然条件下土类处于不同赋水状态时,热扩散系数随含水率变化规律不同相关。此外,还可发现土体天然含水率与其干密度有强烈的负相关性,冻融土热扩散系数之间具有显著正线性相关性。

3 冻融土热扩散系数预测

3.1 冻融土热扩散系数预测模型

3.1.1 二元拟合

对走廊带内冻融土热扩散系数与干密度、天然含水率之间的拟合关系进行曲线估计后发现,基本都呈为多项式函数形式,得到走廊带内冻融土热扩散系数一般拟合公式如下:

α=a1+b1×ω+c1×ω2+d1×ω3+e1×ρd+f1×ρd2

式中:α为土体热扩散系数,10-6 m2·s-1ω为天然含水率,%;ρd为干密度,g·m-3a1b1c1d1e1f1均为待定的拟合系数,各土类拟合值见表3所示。

表 3   热扩散系数二元拟合相关参数

Table 3  Related parameters of binary fitting for thermal diffusivity

土类土壤状态参数
a1b1c1d1e1f1
黏性土融土0.581-0.005-1.591×10-43.395×10-6-0.0040.023
冻土1.450-0.0111.051×10-4-3.487×10-7-0.4860.157
粉土融土-0.1270.008-8.024×10-41.135×10-50.851-0.237
冻土0.4040.036-8.191×10-45.180×10-6-0.1630.161
全风化岩类融土2.7470.036-2.312×10-33.091×10-5-2.3280.602
冻土2.1810.044-1.603×10-31.286×10-5-1.2930.258
砂土融土0.8920.056-2.895×10-33.630×10-5-0.7920.297
冻土1.7520.041-1.163×10-38.211×10-6-1.1870.382
碎石土融土-1.1010.063-5.069×10-31.221×10-41.486-0.330
冻土1.8720.031-1.007×10-36.090×10-6-0.8720.211

新窗口打开| 下载CSV


3.1.2 基于融土热扩散系数的冻土三元拟合

青藏高原走廊带内土类分布多样且分散,各土类矿物成分及组分粒径等分布界限并不一定明显,土类间存在相互混杂现象,所以简单将其分类后采用干密度和含水率对其热扩散系数进行预测会产生一定偏差。此外,相对于冻土而言,融土试样制备及测试更为简单,且融土热扩散系数结果隐含有决定土性成分和粒径等信息,前述偏相关分析表明二者具有极强的正线性相关性,同时考虑干密度、含水率对冻土热扩散系数预测结果进行修正,得到走廊带内冻土热扩散系数的三元拟合公式如下:

αf=a2+b2×αu+c2×lnω+d2×lnρd

式中:αfαu分别为冻土和融土热扩散系数,10-6 m2·s-1a2b2c2d2均为待定的拟合系数,各土类拟合结果见表4所示。

表 4   热扩散系数三元拟合相关参数

Table 4  Related parameters of ternary fitting for thermal diffusivity

土类参数
a2b2c2d2
黏性土-0.6641.7690.240-0.055
粉土-1.3101.9650.3830.222
全风化岩类-0.1481.4790.201-0.399
砂土-1.6771.6860.2610.038
碎石土-0.7991.5040.3100.234

新窗口打开| 下载CSV


3.1.3 径向基函数(RBF)神经网络

径向基函数(radical basis function)神经网络具有计算量小、学习收敛速度快及泛化能力强的特点,其较强的非线性函数逼近能力在诸多领域得以广泛应用29-30。将干密度、天然含水率及融土热扩散系数(冻土热扩散系数三元回归)作为RBF神经网络的输入层,热扩散系数作为输出层,建立以高斯函数为隐含层激活函数的前馈型网络模型。同时,将走廊带内冻融土热扩散系数测试结果按9:1随机分块,其中90%样本数据用于神经网络模型的训练样本,剩余10%用作验证神经网络的预测能力。根据预测结果的反馈,对偏差较大的样本进行剔除,从而有效提高预测精度。图5为全风化岩类融土热扩散系数的RBF神经网络模型预测结果,可以看到预测值与实测值较为一致(R2=0.78,数据剔除比例为9%)。

图5

图5   全风化岩类融土热扩散系数RBF神经网络预测结果

Fig.5   Prediction results of thermal diffusivity based on RBF neural network of unfrozen soil of fully weathered rocks


3.2 融土热扩散系数预测模型对比

图6为各土类融土热扩散系数两种预测模型的预测结果,可以看到两种预测模型所得预测值与实测结果整体较为吻合,大多数样本点分布在相对误差为10%范围内,证明了预测方法的有效性和工程应用价值。此外,对比两种方法预测值与实测值间差异分布,可发现RBF神经网络模型整体预测效果更佳。

图6

图6   各土类融土热扩散系数预测结果对比

Fig.6   Comparison of prediction results of thermal diffusivity of unfrozen soil of different soil types


3.3 冻土热扩散系数预测模型对比

图7为各土类冻土热扩散系数四种预测模型的预测结果,同样可看到四种预测模型具有良好的预测性能。同时,对比四种方法预测值分布范围及与实测值间差异可发现,将融土热扩散系数纳入回归模型的三元预测方法较二元预测方法具有更宽的样本分布范围,且整体预测精度更高。

图7

图7   各土类冻土热扩散系数预测结果对比

Fig.7   Comparison of prediction results of thermal diffusivity of frozen soil of different soil types


3.4 走廊带内冻融土热扩散系数预测方法对比及应用总结

表5表6分别为冻融土热扩散系数预测模型所得结果的相关系数R2和相对误差在10%内占比,由表可知融土热扩散系数两种预测方法中RBF神经网络在剔除样本比例少的情况下,各土类预测精度均显著高于二元拟合方法,其中粉土融土预测效果最好,相关系数R2和相对误差10%内占比分别为0.90和82%。同时,对比冻土热扩散系数预测模型可发现,三元模型预测效果明显优于二元模型,进一步说明融土热扩散系数包含有重要的土壤成分信息,可有效弥补不完全参数回归模型的缺陷。此外,三元RBF神经网络方法要优于三元拟合方法,说明RBF神经网络方法可高效捕捉热扩散系数与影响因素间关系特征,进而有效提升预测精度,其中全风化岩类预测效果最好,相关系数R2和相对误差10%内占比分别为0.96和90.77%。

表 5   热扩散系数预测模型相关系数R2

Table 5  Correlation coefficient R2 of thermal diffusivity prediction models

预测模型剔除比例/%R2
黏性土粉土全风化岩类砂土碎石土
融土冻土融土冻土融土冻土融土冻土融土冻土
二元经验公式拟合15~200.450.250.370.320.520.500.510.250.250.26
二元RBF神经网络5~100.650.360.900.840.780.710.670.550.510.66
三元经验公式拟合3~50.750.800.800.870.84
三元RBF神经网络2~30.860.900.960.910.93

新窗口打开| 下载CSV


表 6   热扩散系数预测模型相对误差在±10%内占比

Table 6  The relative error of the thermal diffusivity prediction models is within ±10%

预测模型拟合参数P±10%(%)
黏性土粉土全风化岩类砂土碎石土
融土冻土融土冻土融土冻土融土冻土融土冻土
二元经验公式拟合ωρd55.3445.9551.1152.6357.1453.1653.3856.2557.1460.31
二元RBF神经网络ωρd65.6052.0582.0070.1586.5772.6269.1761.0868.6467.11
三元经验公式拟合ωρdαu66.9866.6776.5670.1880.17
三元RBF神经网络ωρdαu81.3684.3190.7784.4893.28

新窗口打开| 下载CSV


依据前述分析可知,青藏工程走廊带各土类冻融土热扩散系数分布离散,不同预测方法应用总结如下:

(1)二元拟合方法针对于热扩散系数主体分布区间,其应用范围较窄且误差较大,但获取形式简单(只需干密度和含水率参数),可满足一般工程估算的相应要求;

(2)三元拟合方法局限于冻土热扩散系数,且不可避免高成本和冗长的实验测试过程,但具有适用样本分布范围广、预测性能好及应用简单的优点;

(3)RBF神经网络需以大量实测样本数据和数学算法为基础,其开放性不强且应用相对复杂,但具有应用范围广、预测精度高的优势,总体来看该方法为热扩散系数最佳预测方法。

4 结论

本文针对青藏工程走廊典型土类的冻融土热扩散系数进行了计算,分析了热扩散系数分布特征和参数影响规律,提出了不同预测模型并对比了预测效果,得出主要结论如下:

(1)走廊带土类以黏性土、砂土和碎石土居多,热扩散系数与土类粒径整体呈正相关性,具体为融土热扩散系数按黏性土、粉土、全风化岩类、砂土及碎石土依次增大,冻土热扩散系数按黏性土、全风化岩类、粉土、碎石土及砂土依次增大。

(2)热扩散系数与干密度及天然含水率相关性随土类及冻融状态差异明显,冻、融土热扩散系数呈显著正线性关系。

(3)融土热扩散系数预测模型中,二元RBF神经网络预测效果显著优于二元经验公式;冻土热扩散系数预测模型中,基于融土热扩散系数的三元预测模型的预测精度明显高于二元模型,且三元RBF神经网络模型预测效果更佳。

(4)综合热扩散系数不同预测模型预测效果和误差分析可得,二元拟合方法形式简单,可满足一般工程估算需求;三元拟合方法需以融土实验结果为基础,适用样本分布范围广且预测精度高;RBF神经网络模型具有最优的预测精度且应用范围最广,为最佳预测模型。

参考文献

Wang ShuangjieJin LongMu Keet al.

Distresses and countermeasures of highway subgrade in plateau permafrost regions

[J]. Strategic Study of CAE, 2017196): 140-146.

[本文引用: 1]

汪双杰金龙穆柯.

高原冻土区公路路基病害及工程对策

[J]. 中国工程科学, 2017196): 140-146.

[本文引用: 1]

Wang TaoZhou GuoqingWang Jianzhouet al.

Statistical characteristics and probabilistic analysis of uncertain settlement of subgrade in permafrost regions

[J]. Cold Regions Science and Technology, 2020175103079.

[本文引用: 1]

Wang QingzhiFang JianhongChao Ganget al.

Influence of pavement type on thermal state of block-stone subgrade in permafrost region

[J]. Journal of Glaciology and Geocryology, 2019415): 1087-1097.

[本文引用: 1]

王青志房建宏晁刚.

路面类型对多年冻土区片块石路基热状态的影响

[J]. 冰川冻土, 2019415): 1087-1097.

[本文引用: 1]

Zhu DongpengDong YuanhongLiu Geet al.

Study of the influence of large-width asphalt-concrete pavement on the thermal characteristics of underlying permafrost

[J]. Journal of Glaciology and Geocryology, 2014364): 845-853.

[本文引用: 1]

朱东鹏董元宏刘戈.

宽幅沥青路面热效应对其下部土体热状态的影响

[J]. 冰川冻土, 2014364): 845-853.

[本文引用: 1]

Liu ZhiyunCui FuqingChen Jianbinet al.

Study on the permafrost heat transfer mechanism and reasonable interval of separate embankment for the Qinghai-Tibet expressway

[J]. Cold Regions Science and Technology, 2019170102952.

[本文引用: 1]

Sun ZheZhao LinHu Guojieet al.

Influence of lower boundary conditions on the numerical simulation of permafrost temperature field changes

[J]. Journal of Glaciology and Geocryology, 2021432): 357-369.

[本文引用: 1]

孙哲赵林胡国杰.

下边界条件对多年冻土温度场变化数值模拟的影响

[J]. 冰川冻土, 2021432): 357-369.

[本文引用: 1]

Yang ShimingTao wenquan. Heat transfer[M]. 4th ed. Higher Education Press2006.

[本文引用: 1]

杨世铭陶文铨. 传热学[M]. 4版. 高等教育出版社2006.

[本文引用: 1]

Guo DonglinYang MeixueWang Huijun.

Characteristics of land surface heat and water exchange under different soil freeze/thaw conditions over the central Tibetan Plateau

[J]. Hydrological Processes, 20112516): 2531-2541.

[本文引用: 1]

Li HuoqingWu XinpingMamtimin Aliet al.

Comparison of the different methods to estimate soil thermal diffusivity, temperature and heat flux in sub-surface of Taklimakan Desert

[J]. Chinese Journal of Soil Science, 2016474): 805-813.

[本文引用: 1]

李火青吴新萍买买提艾力·买买提依明.

塔克拉玛干沙漠地表浅层土壤热扩散系数、温度和热通量计算方法的比较研究

[J]. 土壤通报, 2016474): 805-813.

[本文引用: 1]

Wang LinlinGao ZhiqiuHorton R.

Comparison of six algorithms to determine the soil apparent thermal diffusivity at a site in the Loess Plateau of China

[J]. Soil Science, 20101752): 51-60.

[本文引用: 1]

Horton RWierenga P JNielsen D R.

Evaluation of methods for determining the apparent thermal diffusivity of soil near the surface

[J]. Soil Science Society of America Journal, 1983471): 25-32.

[本文引用: 1]

Miao YucongLiu ShuhuaShihua et al.

A comparative study of computing methods of soil thermal diffusivity, temperature and heat flux

[J]. Chinese Journal of Geophysics, 2012552): 441-451.

[本文引用: 1]

缪育聪刘树华吕世华.

土壤热扩散率及其温度、热通量计算方法的比较研究

[J]. 地球物理学报, 2012552): 441-451.

[本文引用: 1]

Han BinghongZhou BingrongSun Yinget al.

The characteristics of soil temperature variations and heat exchange in freeze-thaw period over southern alpine grasslands in Qinghai Province

[J]. Climate Change Research, 2018145): 456-464.

[本文引用: 1]

韩炳宏周秉荣孙瑛.

青海南部冻融区高寒草地土壤温度变化及热量传输特征

[J]. 气候变化研究进展, 2018145): 456-464.

[本文引用: 1]

Gao ZhiqiuFan XingangBian Lingen.

An analytical solution to one-dimensional thermal conduction-convection in soil

[J]. Soil Science, 20031682): 99-107.

[本文引用: 1]

Gao ZhiqiuBian LingenZhang Yabinet al.

Study on analytical resolution to soil thermal conductive equation and soil thermal diffusivity over Nagqu area

[J]. Acta Meteorologica Sinica, 2002603): 352-360.

高志球卞林根张雅斌.

土壤热传导方程解析解和那曲地区土壤热扩散率研究

[J]. 气象学报, 2002603): 352-360.

Dai ChengyingGao ZhiqiuWang Linlinet al.

Intercomparison between two soil temperature algorithms

[J]. Chinese Journal of Atmospheric Sciences, 2009331): 135-144.

[本文引用: 1]

代成颖高志球王琳琳.

两种土壤温度算法的对比分析

[J]. 大气科学, 2009331): 135-144.

[本文引用: 1]

Yuan LimingZhao LinHu Guojieet al.

Hydro-thermal dynamic and soil thermal diffusivity characteristics of typical active layer on the central Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2020422): 378-389.

[本文引用: 1]

原黎明赵林胡国杰.

青藏高原中部典型下垫面活动层水热动态及其热扩散率研究

[J]. 冰川冻土, 2020422): 378-389.

[本文引用: 1]

Zhang YonghuiGao ZhiqiuTong Bing.

Estimating the soil thermal parameters and modeling the soil temperature in Litang of Qinghai-Xizang Plateau

[J]. Plateau Meteorology, 2017361): 33-44.

[本文引用: 1]

章永辉高志球童兵.

青藏高原理塘地区土壤热参数的确定及其土壤温度模拟试验

[J]. 高原气象, 2017361): 33-44.

[本文引用: 1]

Zhen ZuolinZhu JianghongZhang Huyuanet al.

Study on the measurement and prediction of thermal properties for sand

[J]. Chinese Journal of Underground Space and Engineering, 2018146): 1577-1586.

[本文引用: 1]

甄作林朱江鸿张虎元.

砂土导热性能测试与预测研究

[J]. 地下空间与工程学报, 2018146): 1577-1586.

[本文引用: 1]

Wang KeliCheng GuodongJiang Changchanget al.

Variation of thermal diffusivity and temperature simulation of soils of vertical heterogeneity in Nagqu prefecture in the Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2007293): 470-474.

[本文引用: 1]

王可丽程国栋江畅畅.

青藏高原那曲地区非均质土壤导温率的变化及土壤温度数值模拟

[J]. 冰川冻土, 2007293): 470-474.

[本文引用: 1]

Zhou YaGao XiaoqingLi Zhengchaoet al.

Spatio-temporal distribution of thermal diffusivity in deep soil in Qinghai-Tibetan Plateau

[J]. Acta Pedologica Sinica, 2018552): 351-359.

[本文引用: 1]

周亚高晓清李振朝.

青藏高原深层土壤热扩散率的时空分布特征

[J]. 土壤学报, 2018552): 351-359.

[本文引用: 1]

Liu JingxingLi ChuntaoLiang Yuxianget al.

Influence of soil texture on the thermal diffusivity

[J]. Journal of Desert Research, 2010303): 577-581.

[本文引用: 1]

刘经星李春桃梁玉祥.

质地对土壤热扩散率的影响

[J]. 中国沙漠, 2010303): 577-581.

[本文引用: 1]

Roxy M SSumithranand V BRenuka G.

Variability of soil moisture and its relationship with surface albedo and soil thermal diffusivity at Astronomical Observatory, Thiruvananthapuram, south Kerala

[J]. Journal of Earth System Science, 20101194): 507-517.

[本文引用: 1]

An KedongWang WenkeWang Zhoufenget al.

Soil thermal parameters of aeolian sand and its relationship with water content

[J]. Hydrogeology & Engineering Geology, 2015425): 129-133.

[本文引用: 1]

安可栋王文科王周峰.

风积沙热参数计算及其与含水率关系研究

[J]. 水文地质工程地质, 2015425): 129-133.

[本文引用: 1]

Di JiayingLiu XiaonaRen Tusheng.

Comparative study on thermal properties of intact and repacked soil samples

[J]. Transactions of the Chinese Society of Agricultural Engineering, 20122821): 74-79.

[本文引用: 1]

邸佳颖刘晓娜任图生.

原状土与装填土热特性的比较

[J]. 农业工程学报, 20122821): 74-79.

[本文引用: 1]

Ma XinZhang TangtangChen Jinlei.

Observations of soil thermal properties in a typical mesa over Chinese loess plateau

[J]. Plateau Meteorology, 2019383): 507-517.

[本文引用: 1]

马欣张堂堂陈金雷.

黄土高原典型塬区土壤热性质变化特征研究

[J]. 高原气象, 2019383): 507-517.

[本文引用: 1]

Arkhangelskaya TLukyashchenko K.

Estimating soil thermal diffusivity at different water contents from easily available data on soil texture, bulk density, and organic carbon content

[J]. Biosystems Engineering, 201816883-95.

[本文引用: 1]

Xu XiaozuWang JiachengZhang Lixin. Permafrost physics[M]. BeijingScience Press2010.

[本文引用: 1]

徐斅祖王家澄张立新. 冻土物理学[M]. 北京科学出版社2010.

[本文引用: 1]

Cong Shuang.

The function analysis and application study of radial basis function network

[J]. Computer Engineering and Applications, 2002383): 85-87.

[本文引用: 1]

丛爽.

径向基函数网络的功能分析与应用的研究

[J]. 计算机工程与应用, 2002383): 85-87.

[本文引用: 1]

Wei MinYu Le’an.

A RBF neural network with optimum learning rates and its application

[J]. Journal of Management Sciences in China, 2012154): 50-57.

[本文引用: 1]

卫敏余乐安.

具有最优学习率的RBF神经网络及其应用

[J]. 管理科学学报, 2012154): 50-57.

[本文引用: 1]

/