青藏高原工程走廊帶凍融砂土導(dǎo)熱系數(shù)分布規(guī)律研究

金嘉林，郝建斌，劉志云，崔福慶

青藏高原工程走廊帶凍融砂土導(dǎo)熱系數(shù)分布規(guī)律研究

金嘉林，郝建斌，劉志云，崔福慶

(長(zhǎng)安大學(xué) 地質(zhì)工程與測(cè)繪學(xué)院，陜西 西安 710054)

為研究青藏高原工程走廊帶內(nèi)凍融砂土的導(dǎo)熱系數(shù)變化規(guī)律，在青藏高速(格爾木~拉薩段)現(xiàn)場(chǎng)鉆探取樣，以不同顆粒級(jí)配的粉砂、中細(xì)砂和礫砂為研究對(duì)象，采用瞬態(tài)面熱源法進(jìn)行重塑凍融砂土的導(dǎo)熱系數(shù)測(cè)試，分析3種砂土在凍融狀態(tài)下干密度、天然含水率及導(dǎo)熱系數(shù)的分布規(guī)律。測(cè)試結(jié)果表明，凍土砂土的導(dǎo)熱系數(shù)分布規(guī)律為：粉砂＜中細(xì)砂＜礫砂；融土砂土的導(dǎo)熱系數(shù)分布規(guī)律為：粉砂＜礫砂＜中細(xì)砂。根據(jù)實(shí)驗(yàn)數(shù)據(jù)擬合出砂土導(dǎo)熱系數(shù)與其干密度、天然含水率關(guān)系的經(jīng)驗(yàn)公式，為青藏高原工程走廊帶內(nèi)道路工程的建設(shè)提供一定的理論依據(jù)。

青藏高原工程走廊帶；凍融砂土；導(dǎo)熱特性

在青藏高原工程走廊帶的凍土區(qū)的工程建設(shè)中，凍土工程的溫度場(chǎng)變化引發(fā)熱力場(chǎng)中的溫度傳遞，從而影響著凍土凍融的速率變化[1]，而土體的導(dǎo)熱系數(shù)對(duì)于凍土區(qū)的寬幅路基溫度場(chǎng)穩(wěn)定起到重要的作用[2]。導(dǎo)熱系數(shù)是影響凍土溫度變化及巖土內(nèi)熱量傳遞的重要因素[3?4]，主要受溫度、土壤孔隙度、含水率等諸多因素的影響[5]。Tarnawski 等[6]對(duì)3種飽和標(biāo)準(zhǔn)砂的導(dǎo)熱系數(shù)進(jìn)行了測(cè)試，結(jié)果表明：飽和導(dǎo)熱系數(shù)隨孔隙度的增加而增加，隨溫度增加而減小。周媛媛等[7]對(duì)含水砂土的熱導(dǎo)率進(jìn)行了測(cè)定，認(rèn)為在飽和含水率以下，砂樣的熱導(dǎo)率均隨著含水率的增加而增大。目前，凍土導(dǎo)熱系數(shù)測(cè)試方法主要有穩(wěn)態(tài)法、瞬態(tài)法和正規(guī)狀態(tài)法[8]。其中瞬態(tài)法以其測(cè)量時(shí)間短的優(yōu)勢(shì)近年來(lái)得到廣泛應(yīng)用，如周媛媛等[7]利用熱線法測(cè)試了含水砂土導(dǎo)熱系數(shù)隨含水率變化的規(guī)律，李棟偉等[9]測(cè)試了人工凍黏土的導(dǎo)熱系數(shù)。Tarnawski等[6]使用非穩(wěn)態(tài)探針測(cè)定3種飽和標(biāo)準(zhǔn)砂的導(dǎo)熱系數(shù)；Manohar 等[10]利用熱探針?lè)y(cè)定了相同密度粗(細(xì))砂的導(dǎo)熱系數(shù)；HU等[11]基于瞬態(tài)平面法原理，采用DRE-III分別測(cè)定了長(zhǎng)江地區(qū)的黏土和砂質(zhì)土的導(dǎo)熱系數(shù)。青藏高原工程走廊帶總體線路較長(zhǎng)、地質(zhì)構(gòu)造復(fù)雜且土體的物理力學(xué)性質(zhì)不穩(wěn)定[12]，目前尚缺乏科學(xué)系統(tǒng)的實(shí)驗(yàn)測(cè)試結(jié)果。本文依托青藏高速初步勘察項(xiàng)目，對(duì)現(xiàn)場(chǎng)鉆探取得的砂土樣本進(jìn)行導(dǎo)熱系數(shù)測(cè)試。針對(duì)此次測(cè)試樣品數(shù)量龐大，本文采用基于瞬態(tài)面熱源技術(shù)的Hot DiskTPS1500熱參數(shù)分析儀進(jìn)行測(cè)試。

1 實(shí)驗(yàn)概況

1.1 實(shí)驗(yàn)設(shè)備簡(jiǎn)介

實(shí)驗(yàn)采用Hot DiskTPS1500熱參數(shù)分析儀進(jìn)行凍融重塑砂土的導(dǎo)熱系數(shù)測(cè)試，另配備有冰柜、密封的樣品架以及有機(jī)玻璃材質(zhì)的實(shí)驗(yàn)盒等輔助儀器。圖1為樣品架及Kapton探頭。

1.2 測(cè)試樣品來(lái)源及其物性參數(shù)

本次實(shí)驗(yàn)依托青藏高速(格爾木~拉薩段)工程地質(zhì)勘察項(xiàng)目，測(cè)試區(qū)段為青藏高原工程走廊帶內(nèi)北起西大灘(即走廊帶多年凍土北界)，南至唐古拉山，共437 km。最大鉆孔深度40 m，每孔沿不同深度進(jìn)行取樣，取樣深度分布見表1。試驗(yàn)中，試樣的土質(zhì)類別、干密度、天然含水率等參數(shù)依照《公路土工試驗(yàn)規(guī)程》(JTG E40—2007)進(jìn)行測(cè)定。試驗(yàn)對(duì)同一砂土樣進(jìn)行凍融測(cè)試，共測(cè)試土樣220組，其中共測(cè)試融土樣133組、凍土樣205組。

圖1 測(cè)試樣品架

表1 青藏高原工程走廊帶內(nèi)凍(融)砂土測(cè)試樣本深度分布

1.2.1 砂土土性分布

對(duì)鉆探獲取的砂土土樣進(jìn)行土性統(tǒng)計(jì)，按照砂土粒徑級(jí)配分為粉砂、中細(xì)砂以及礫砂3類，各類土樣數(shù)目如圖2所示。由圖2可知，走廊帶內(nèi)砂土主要以礫砂和中細(xì)砂為主，粉砂相對(duì)較少。

圖2 砂土分類

1.2.2 干密度與天然含水率分布特征

對(duì)3類土的干密度進(jìn)行統(tǒng)計(jì)，得到其概率分布如圖3所示。由圖可知，礫砂土和中細(xì)砂的干密度分布特征基本接近，礫砂土、中細(xì)砂和粉砂的主要分布區(qū)間(累積分布概率20%~80%)分別為1.51~ 1.90，1.60~1.93和1.45~1.71 g/cm3，均值分別為1.71，1.75和1.59 g/cm3。走廊帶內(nèi)砂土的干密度分布規(guī)律為：粉砂＜礫砂＜中細(xì)砂。

對(duì)粉砂、中細(xì)砂和礫砂的天然含水率進(jìn)行統(tǒng)計(jì)，得到其概率分布如圖4所示。由圖可知，中細(xì)砂和礫砂的天然含水率分布較為接近，粉砂、中細(xì)砂和礫砂的天然含水率主要分布區(qū)間(累積分布概率20%~80%)分別為17.1%~25.6%，11.5%~26.1%和8.5%~33.6%，均值分別為24.09%，19.44%和23.5%。走廊帶內(nèi)砂土天然含水率分布規(guī)律為：粉砂＞礫砂＞中細(xì)砂。

圖3 砂土干密度概率分布

圖4 砂土天然含水率概率分布

1.3 試樣制備及測(cè)試

對(duì)于未飽和土，將其烘干并進(jìn)行粉碎和過(guò)篩后，按既定含水量(天然含水率)進(jìn)行配制，混合均勻，然后用塑料薄膜遮蓋悶樣8~12 h后，裝入制樣模具(圖5)，采用靜壓法進(jìn)行制樣(圖6)；對(duì)于飽冰和含土冰層土樣，則需加入破碎后的冰沙來(lái)進(jìn)行制樣。為保護(hù)測(cè)試傳感器且保證測(cè)試數(shù)據(jù)的可靠性，試樣取出后，再將試樣表面進(jìn)行磨平處理。

圖5 制樣模具

圖6 制樣機(jī)

試驗(yàn)中，將含水率大于15%的砂土試樣稱為大含水量試樣，僅測(cè)定凍土的導(dǎo)熱系數(shù)；其余為小含水量試樣，先測(cè)定融土的導(dǎo)熱系數(shù)，再測(cè)定凍土的導(dǎo)熱系數(shù)。融土樣的測(cè)試溫度控制在20 ℃左右；對(duì)于凍土樣，先將融土樣置于冰柜內(nèi)冷凍24 h后進(jìn)行測(cè)試，其中柜內(nèi)溫度范圍控制在?10 ℃左右。

2 凍(融)砂土導(dǎo)熱系數(shù)概率分布及影響因素分析

2.1 凍(融)砂土導(dǎo)熱系數(shù)概率分布統(tǒng)計(jì)

對(duì)3類砂土的凍土導(dǎo)熱系數(shù)進(jìn)行概率分布統(tǒng)計(jì)，得到3種砂土的導(dǎo)熱系數(shù)概率分布規(guī)律，如圖7所示。從圖中可看出，礫砂和中細(xì)砂導(dǎo)熱系數(shù)概率分布比較接近，但總體上礫砂導(dǎo)熱系數(shù)高于中細(xì)砂導(dǎo)熱系數(shù)，粉砂導(dǎo)熱系數(shù)在累積分布概率區(qū)間20%~80%之間導(dǎo)熱系數(shù)相對(duì)最小。粉砂、中細(xì)砂和礫砂的凍土導(dǎo)熱系數(shù)主要分布區(qū)間(累積分布概率20%~80%)分別為1.96~2.48，1.75~2.73和1.74~2.78 W/(m?K)，均值分別為2.19，2.28 和2.36 W/(m?K)?？梢?，砂土凍土導(dǎo)熱系數(shù)隨其粒徑的增大而增大，青藏高原走廊帶內(nèi)砂土凍土的導(dǎo)熱系數(shù)分布基本排序?yàn)椋悍凵埃贾屑?xì)砂＜礫砂。

圖7 凍土導(dǎo)熱系數(shù)概率分布

對(duì)3類砂土的融土導(dǎo)熱系數(shù)進(jìn)行概率分布統(tǒng)計(jì)，結(jié)果如圖8所示。粉砂、中細(xì)砂和礫砂主要分布區(qū)間(累積分布概率20%~80%)分別為1.44~1.98，1.35~2.34和1.32~2.17 W/(m?K)，均值分別為1.69，1.90和1.79 W/(m?K)。認(rèn)為走廊帶內(nèi)砂土的融土導(dǎo)熱系數(shù)與土的成分有顯著關(guān)系，砂土融土導(dǎo)熱系數(shù)分布基本排序?yàn)椋悍凵埃嫉[砂＜中細(xì)砂。

2.2 凍(融)砂土導(dǎo)熱系數(shù)影響因素分析

天然含水率和干密度是影響導(dǎo)熱系數(shù)的2個(gè)重要影響因素。土的持水性受干密度的影響，干密度越大其持水性越弱，導(dǎo)致天然含水率通常與干密度呈負(fù)相關(guān)關(guān)系。為了研究單一因素對(duì)導(dǎo)熱系數(shù)的影響及相關(guān)程度，分析其與導(dǎo)熱系數(shù)的密切程度，對(duì)走廊帶內(nèi)凍融砂土導(dǎo)熱系數(shù)與干密度和天然含水率進(jìn)行統(tǒng)計(jì)，得出3類砂土導(dǎo)熱系數(shù)與天然含水率、干密度之間的偏相關(guān)性，結(jié)果如表2所示。

圖8 融土導(dǎo)熱系數(shù)概率分布

表2 凍融砂土導(dǎo)熱系數(shù)影響因素及其偏相關(guān)性

由表2可知，在控制其中一個(gè)影響因素(干密度或天然含水率)的情況下，另一影響因素與3種凍融砂土導(dǎo)熱系數(shù)均有顯著的正相關(guān)關(guān)系，其顯著性(雙側(cè))均在0.01以上，說(shuō)明隨著干密度和天然含水率的增大，其導(dǎo)熱系數(shù)相應(yīng)增大。這是由于干密度增大時(shí)，土越密實(shí)，單位體積土顆粒中礦物骨架數(shù)量增加，孔隙減少導(dǎo)致土顆粒之間緊密接觸，使土體的導(dǎo)熱性能增強(qiáng)；當(dāng)天然含水率增大時(shí)，由于水的導(dǎo)熱系數(shù)較空氣的導(dǎo)熱系數(shù)大，相應(yīng)地提高了土體的導(dǎo)熱系數(shù)。

3 凍(融)砂土導(dǎo)熱系數(shù)擬合

3.1 凍(融)砂土導(dǎo)熱系數(shù)擬合

對(duì)走廊帶內(nèi)凍融砂土與干密度、天然含水率之間的擬合關(guān)系式進(jìn)行曲線估計(jì)后發(fā)現(xiàn)，基本都呈現(xiàn)為指數(shù)函數(shù)形式，得出導(dǎo)熱系數(shù)與含水量、干密度相關(guān)性的擬合公式為：

式中：λ和λ分別為融土和凍土導(dǎo)熱系數(shù)，W/(m?K)；ρ為干密度，g/m3；為天然含水率；，和為待定的擬合系數(shù)。

利用以上擬合關(guān)系式對(duì)3類凍融砂土的導(dǎo)熱系數(shù)分別與含水量和干密度之間的關(guān)系進(jìn)行擬合，分別得出3類凍融砂土，和各參數(shù)值，見表3所示。

表3 凍融砂土導(dǎo)熱系數(shù)擬合結(jié)果

3.2 凍(融)砂土導(dǎo)熱系數(shù)擬合誤差分析

3類砂土的擬合效果如圖9~11所示。由圖可看出，粉砂凍融土導(dǎo)熱系數(shù)的擬合效果稍低于中細(xì)砂和礫砂的擬合效果，各類凍融土導(dǎo)熱系數(shù)的擬合誤差在各區(qū)間內(nèi)比較均勻，且在可接受范圍內(nèi)，可見擬合效果較好。

(a) 融土；(b) 凍土

為了進(jìn)一步對(duì)擬合參數(shù)值進(jìn)行驗(yàn)證，將3類砂土樣本中的含水量和干密度值分別代入式(4)和式(5)，通過(guò)與實(shí)測(cè)值比較，得出導(dǎo)熱系數(shù)相對(duì)誤差值，對(duì)其作相對(duì)誤差概率分布，如圖12~13所示。

由圖9~10可知，粉砂、中細(xì)砂和礫砂的融土導(dǎo)熱系數(shù)相對(duì)誤差在±15%以內(nèi)的累計(jì)誤差概率分別為11%~95%，5%~95%和14%~86%之間，即分別有84%，90%和72%的樣本相對(duì)誤差在±15%以內(nèi)；而凍土導(dǎo)熱系數(shù)相對(duì)誤差在±15%以內(nèi)的累計(jì)誤差概率分別為11%~95%，21%~91%和12%~95%之間，且分別有84%，70%和83%的樣本相對(duì)誤差在±15%以內(nèi)，則3類砂土導(dǎo)熱系數(shù)的擬合效果較為理想。

(a) 融土；(b) 凍土

(a) 融土；(b) 凍土

圖12 融土導(dǎo)熱系數(shù)誤差概率分布

圖13 凍土導(dǎo)熱系數(shù)誤差概率分布

4 結(jié)論

1) 采用瞬態(tài)平面法對(duì)青藏高原工程走廊帶內(nèi)重塑凍融砂土的導(dǎo)熱系數(shù)進(jìn)行了測(cè)定，結(jié)果表明青藏工程走廊帶內(nèi)粉砂、中細(xì)砂及礫砂的干密度與其含水率呈負(fù)相關(guān)關(guān)系；砂土干密度和含水率與其導(dǎo)熱系數(shù)均有顯著的正相關(guān)關(guān)系。

2) 在凍結(jié)狀態(tài)下，3類砂土的導(dǎo)熱系數(shù)大小關(guān)系為：粉砂＜中細(xì)砂＜礫砂；融化狀態(tài)下，3類砂土的導(dǎo)熱系數(shù)大小關(guān)系為：粉砂＜礫砂＜中細(xì)砂；相同條件下，凍土導(dǎo)熱系數(shù)大于融土導(dǎo)熱系數(shù)。

3) 根據(jù)試驗(yàn)數(shù)據(jù)，對(duì)青藏高原走廊帶內(nèi)3類凍融砂土的導(dǎo)熱系數(shù)公式進(jìn)行了擬合，擬合結(jié)果較為理想。

[1] 姜雄. 多年凍土區(qū)高溫凍土導(dǎo)熱系數(shù)試驗(yàn)研究[D]. 徐州: 中國(guó)礦業(yè)大學(xué), 2015. JIANG Xiong. Experimental study on thermal conductivity of high temperature frozen soil in permafrost regions[D]. Xuzhou: China University of Mining and Technology, 2015.

[2] 劉戈, 汪雙杰, 袁堃, 等. 復(fù)合措施在多年凍土區(qū)寬幅路基建設(shè)中的適用性[J]. 公路, 2016, 61(3): 12?17. LIU Ge, WANG Shuangjie, YUAN Wei, et al. Applicability of composite measures in wide-span subgrade construction in permafrost regions[J]. Highway, 2016, 61(3): 12?17.

[3] TIAN Z, LU Y, Horton R, et al. A simplified devries-based model to estimate thermal conductivity of unfrozen and frozen soil[J]. European Journal of Soil Science, 2016, 67(5):564?572.

[4] 朱志武, 寧建國(guó), 馬巍. 土體凍融過(guò)程中水、熱、力三場(chǎng)耦合本構(gòu)問(wèn)題及數(shù)值分析[J]. 工程力學(xué), 2007, 24(5): 138?144. ZHU Zhiwu, NING Jianguo, MA Wei. Three-field coupling constitutive problems of water, heat and force in soil freezing and thawing process and numerical analysis[J]. Engineering Mechanics, 2007, 24(5): 138? 144.

[5] 杜宜臻, 李韌, 吳通華, 等. 土壤熱導(dǎo)率的研究現(xiàn)狀及其進(jìn)展[J]. 冰川凍土, 2015, 37(4): 1067?1074. DU Yizhen, LI Ren, WU Tonghua, et al. Research status and progress of soil thermal conductivity[J]. Journal of Glaciology and Geocryology, 2015, 37(4): 1067?1074.

[6] Tarnawski V R, Momose T, Leong W H. Thermal conductivity of standard sands ii. saturated conditions[J]. International Journal of Thermophysics, 2011, 32(5): 984?1005.

[7] 周媛媛, 江海峰, 李輝, 等. 含水砂土導(dǎo)熱系數(shù)實(shí)驗(yàn)研究[J]. 工程熱物理學(xué)報(bào), 2015, 36(2): 265?268. ZHOU Yuanyuan, JIANG Haifeng, LI Hui, et al. Experimental study on thermal conductivity of water-containing sand[J]. Journal of Engineering Thermophysics, 2015, 36(2): 265?268.

[8] 陳之祥, 李順群, 夏錦紅, 等. 凍土導(dǎo)熱系數(shù)測(cè)試和計(jì)算現(xiàn)狀分析[J]. 建筑科學(xué)與工程學(xué)報(bào), 2019, 36(2): 101?115. CHEN Zhixiang, LI Shunqun, XIA Jinhong, et al. Analysis of the current situation of thermal conductivity test and calculation of frozen soil[J]. Journal of Architecture and Civil Engineering, 2019, 36(2): 101? 115.

[9] 李棟偉, 劉娟, 蔡強(qiáng), 等. 人工凍融黏土導(dǎo)熱系數(shù)試驗(yàn)研究[J]. 地下空間與工程學(xué)報(bào), 2017, 13(6): 1452?1456. LI Dongwei, LIU Juan, CAI Qiang, et al. Experimental study on thermal conductivity of artificial freeze-thaw clay[J]. Chinese Journal of Underground Space and Engineering, 2017, 13(6): 1452?1456.

[10] Manohar K, Shirish R. Apparent thermal conductivity of sand[J]. International Journal of Heat Exchangers, 2005, 6:250?259.

[11] HU Q, LI J, LIN J, et al. Experimental study of soil thermal properties based on the transient method[C]// Proceedings of the 2017 6th International Conference on Energy, Environment and Sustainable Development. Atlantis Press, 2017: 372?380.

[12] 趙濤. 青藏高原凍土工程走廊場(chǎng)地地震動(dòng)特征與路基動(dòng)力響應(yīng)分析[D]. 蘭州: 蘭州交通大學(xué), 2017. ZHAO Tao. Ground motion characteristics and subgrade dynamic response analysis of the Qinghai-Tibet Plateau permafrost engineering corridor[D]. Lanzhou: Lanzhou Jiaotong University, 2017.

Distribution law of thermal conductivity offrozen-thawed sandy soil in Qinghai-Tibet plateau engineering corridor

JIN Jialin, HAO Jianbin, LIU Zhiyun, CUI Fuqing

(School of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China)

In order to study the variation law of thermal conductivity of frozen-thawed sandy soil in the Qinghai-Tibet Plateau Engineering Corridor, drilling and sampling were carried out in the Qinghai-Tibet Expressway (Germu-Lhasa section), with silt, medium-fine sand and gravel sand of different grain compositions as the research object. The thermal conductivity of reconstituted frozen-thawed sand was tested by Transient Plane Source Method. The distribution of dry density, natural moisture content and thermal conductivity of three kinds of sand under freezing and thawing conditions was analyzed. The test results show that the distribution law of thermal conductivity of frozen sand soil is: silt

Qinghai-Tibet Plateau engineering corridor; frozen-thawed sandy soil; thermal conduction characteristic

U419

A

1672 ? 7029(2020)03 ? 0608 ? 07

10.19713/j.cnki.43?1423/u.T20190537

2019?06?17

國(guó)家自然科學(xué)基金資助項(xiàng)目(51574037，41502292，41472266)

劉志云(1979?)，男，江西吉安人，副教授，博士，從事寒區(qū)工程熱災(zāi)害與防治方面的研究；E?mail：dcdgx33@chd.edu.cn

(編輯 涂鵬)