羅曉曉，俞祁浩，馬勤國(guó)，張建明

?

基于附面層理論的路基熱邊界模型

羅曉曉1, 2，俞祁浩1，馬勤國(guó)3，張建明1

(1. 中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院 凍土工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，甘肅 蘭州，730000；2. 中國(guó)科學(xué)院大學(xué) 工程科學(xué)學(xué)院，北京，100049；3. 華南理工大學(xué) 亞熱帶建筑科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，廣東 廣州，510641)

基于“附面層原理”，建立考慮氣象條件和路面結(jié)構(gòu)熱物理因素影響下的路面結(jié)構(gòu)熱交換計(jì)算模型。選用合適的判據(jù)給出熱邊界厚度的確定方法，然后，通過(guò)數(shù)值模擬和灰關(guān)聯(lián)分析，考察太陽(yáng)輻射強(qiáng)度、風(fēng)速、路面結(jié)構(gòu)等因素對(duì)熱邊界底部溫度增量的影響。最后，對(duì)各影響因素及其交互作用下的熱邊界底部溫度增量進(jìn)行二次正交回歸分析，建立熱邊界底部溫度計(jì)算模型，并以青藏高原氣象數(shù)據(jù)和青藏高速公路路基試驗(yàn)段實(shí)測(cè)的溫度驗(yàn)證模型的可靠性。研究結(jié)果表明：熱邊界厚度與路基等效熱擴(kuò)散率線性正相關(guān)；路面吸收率對(duì)熱邊界底部溫度增量的影響最大，太陽(yáng)輻射強(qiáng)度的影響次之，對(duì)流換熱系數(shù)的影響中等，等效熱擴(kuò)散率的影響較小，發(fā)射率的影響最?。辉撃Ｐ涂煽啃愿?，可以用來(lái)預(yù)測(cè)路基熱邊界。

附面層；熱邊界；敏感性分析；正交分析法；預(yù)測(cè)模型

與天然地表相比，瀝青路面具有強(qiáng)烈的吸熱作用和較弱的蒸發(fā)作用[1]。瀝青路面因顯著的溫度敏感性常常會(huì)產(chǎn)生車(chē)轍、推移、擁包等現(xiàn)象。在多年凍土區(qū)，瀝青路面的這種強(qiáng)吸熱性會(huì)造成路基下方多年凍土融化、上限下移，進(jìn)而引發(fā)路基下沉、路面開(kāi)裂、凹陷等病害，嚴(yán)重影響行車(chē)安全[2?4]，可見(jiàn)，溫度對(duì)路基熱穩(wěn)定性有顯著的影響。路基溫度場(chǎng)模擬是研究路基熱穩(wěn)定性的基礎(chǔ)，而路基熱邊界條件的確定是路基溫度場(chǎng)分析的關(guān)鍵[5]。路基熱邊界是路表與外界環(huán)境能量交換的結(jié)果，許多學(xué)者采用能量平衡法進(jìn)行估算[5?7]。但受太陽(yáng)輻射強(qiáng)度、風(fēng)速、大氣溫度等氣象條件和路面吸收率、反射率、導(dǎo)熱系數(shù)、比熱容等路面結(jié)構(gòu)熱物理參數(shù)的綜合影響，對(duì)路基溫度場(chǎng)進(jìn)行長(zhǎng)期模擬時(shí)，要精確確定路基熱邊界相當(dāng)困難。目前，對(duì)寒區(qū)工程熱邊界的處理方法主要如下3種：一是現(xiàn)場(chǎng)實(shí)測(cè)溫度，如凍土工程國(guó)家重點(diǎn)實(shí)驗(yàn)室在青藏鐵路沿線布設(shè)37個(gè)測(cè)溫孔，其中多年凍土區(qū)30個(gè)，融區(qū)6個(gè)，季節(jié)凍土區(qū)1個(gè)[8]。這種方法受限于測(cè)溫孔以地溫監(jiān)測(cè)為主，往往不能確定路表的溫度并且測(cè)溫傳感器容易失效。二是理論分析加數(shù)值模擬方法，這種分析方法根據(jù)傳熱學(xué)原理，將熱邊界作為隱含量，視大氣、土體為整體系統(tǒng)，來(lái)分析路基溫度場(chǎng)[9?12]。這種方法需要詳細(xì)的氣象、地質(zhì)資料，在工程應(yīng)用時(shí)需確定的參數(shù)較多，會(huì)有諸多不便。三是統(tǒng)計(jì)分析研究，即利用氣象統(tǒng)計(jì)資料，分析建立路基表面溫度與大氣溫度等因素間的關(guān)系，以此來(lái)研究溫度場(chǎng)[13?15]。這種方法需要對(duì)大量的氣象數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，并且統(tǒng)計(jì)結(jié)果帶有一定的區(qū)域局限性，同時(shí)，這種分析方法不考慮氣象條件和路面結(jié)構(gòu)熱物理參數(shù)各因素間的交互作用，所得經(jīng)驗(yàn)公式誤差相對(duì)偏大。朱林楠[16]經(jīng)過(guò)長(zhǎng)期的氣候及地溫觀測(cè)，將流體力學(xué)中的附面層概念引入巖土工程，提出了路基溫度場(chǎng)計(jì)算的“附面層理論”。該理論定義路表上下一定范圍內(nèi)溫度、濕度受大氣及太陽(yáng)輻射作用強(qiáng)烈的物理層為“附面層”，附面層底部溫度受太陽(yáng)輻射、地面紊流熱交換等因素影響較小，以致可以忽略，進(jìn)而建立了附面層溫度與大氣溫度之間的經(jīng)驗(yàn)關(guān)系，即附面層底部溫度等于大氣平均溫度加上1個(gè)溫度增值。附面層理論可用于有效地分析路基溫度 場(chǎng)[17?19]，但在應(yīng)用時(shí)存在如下問(wèn)題，1) 對(duì)附面層厚度的確定沒(méi)有統(tǒng)一標(biāo)準(zhǔn)；2) 附面層底部溫度均以現(xiàn)場(chǎng)實(shí)測(cè)地溫為依據(jù)，通過(guò)函數(shù)擬合得到，這在溫度場(chǎng)分析時(shí)極為不便；3) 附面層理論雖然對(duì)附面層厚度和附面層底部溫度的確定給出了相應(yīng)的算法，但對(duì)其中參數(shù)的確定缺少理論分析和科學(xué)依據(jù)，實(shí)際應(yīng)用時(shí)無(wú)法直接采用模型計(jì)算。為此，本文作者基于附面層理論，綜合考慮影響路面溫度場(chǎng)的氣象條件和路面結(jié)構(gòu)熱物理參數(shù)各因素，建立路面結(jié)構(gòu)熱交換計(jì)算模型；以附面層厚度作為路基熱邊界厚度，以附面層底部溫度作為熱邊界底部溫度，通過(guò)數(shù)值模擬和灰關(guān)聯(lián)分析，考察太陽(yáng)輻射強(qiáng)度、風(fēng)速、路面結(jié)構(gòu)等因素對(duì)熱邊界底部溫度增量的影響，經(jīng)理論分析，選用合適的判據(jù)給出熱邊界厚度的確定方法，并利用二次正交回歸對(duì)太陽(yáng)輻射強(qiáng)度、風(fēng)速、大氣溫度等氣象條件和路面吸收率、反射率、導(dǎo)熱系數(shù)、比熱容等路面結(jié)構(gòu)熱物理參數(shù)交互變化時(shí)的熱邊界底部溫度增量進(jìn)行研究，構(gòu)建可以有效預(yù)測(cè)路基熱邊界條件的模型，為路基熱邊界的確定提供合理依據(jù)。

1 路面結(jié)構(gòu)熱交換計(jì)算模型

由于瀝青材料的小孔隙度阻礙路基面層下方水分的蒸發(fā)過(guò)程，路基表層不能形成有效的湍流散熱[20]。降雨、降雪對(duì)地表能量傳輸?shù)挠绊懼饕w現(xiàn)在地表反照率的變化方面，反照率的變化主要受地面持水性和入滲特性的影響，與天然地表相比，瀝青路面具有低孔隙度和低滲透性，其反照率是天然地表的30%左 右[21]。瀝青路面的強(qiáng)吸熱作用會(huì)使路面積雪保存時(shí)間較短。鑒于上述原因，瀝青路面結(jié)構(gòu)傳熱分析時(shí)，忽略水汽蒸發(fā)耗熱和降雨、降雪的影響。路面與周?chē)h(huán)境熱交換過(guò)程主要包括熱傳導(dǎo)、對(duì)流和輻射。

1.1 輻射換熱

路表與環(huán)境的輻射換熱過(guò)程包括路表吸收的太陽(yáng)輻射S和路表反射的長(zhǎng)波輻射F?？紤]到大氣對(duì)太陽(yáng)輻射的反射作用，路表實(shí)際吸收的太陽(yáng)輻射S表示為

式中：為參數(shù)，=12/0，0為實(shí)際有效日照時(shí)間；0為太陽(yáng)日輻射強(qiáng)度最大值，0=0.131d，d為太陽(yáng)日輻射總強(qiáng)度，為角頻率，=2π/24；為時(shí)間。

路表向外發(fā)射的長(zhǎng)波輻射F可以表示為[23]

式中：為路面發(fā)射率；為玻爾茲曼常數(shù)；為路表溫度；k為絕對(duì)零度；a為大氣溫度。大氣溫度的日變化可以表示為

1.2 對(duì)流換熱

空氣對(duì)流帶給路表的熱量可以表示為

式中：為對(duì)流換熱系數(shù)，W/(m2·℃)。按照經(jīng)驗(yàn)公式，其與風(fēng)速的關(guān)系為[23]

式中：w為日平均風(fēng)速，m/s。

1.3 熱傳導(dǎo)

路基內(nèi)的熱傳導(dǎo)問(wèn)題實(shí)際是路基面層、基層和土基等各層組成的層狀結(jié)構(gòu)的溫度場(chǎng)問(wèn)題。假設(shè)路基層間接觸良好，熱阻為0，路基內(nèi)熱傳導(dǎo)方程表示為

根據(jù)能量守恒原理，該問(wèn)題的有限元計(jì)算公式 為[24]

1.4 路面溫度場(chǎng)有限元模型

以某路段實(shí)際路面結(jié)構(gòu)和材料情況為路面溫度場(chǎng)的計(jì)算模型，路基結(jié)構(gòu)從上至下依次為改性瀝青、重交瀝青、水泥穩(wěn)定碎石基層、石灰土和土基層，各層厚度如圖1所示。路基各層材料熱物理參數(shù)與氣象數(shù)據(jù)皆采用文獻(xiàn)[25]中實(shí)測(cè)值，如表1所示，當(dāng)?shù)靥?yáng)輻射強(qiáng)度、風(fēng)速及氣溫等氣象數(shù)據(jù)如表2所示。采用ANSYS有限元軟件建立該路段路基傳熱模型，模型中瀝青路面太陽(yáng)輻射吸收率取0.9，路面反射率取0.81，下邊界熱通量溫度梯度取0.03℃/m[25]，即熱流取 37.5 mW/m2。

圖1 路面結(jié)構(gòu)示意圖

表1 路面結(jié)構(gòu)各層材料熱物性參數(shù)

表2 當(dāng)?shù)貧庀髤?shù)

1.5 計(jì)算結(jié)果驗(yàn)證

取計(jì)算結(jié)果中重交瀝青AC-20上層點(diǎn)和下層點(diǎn)(圖1)溫度計(jì)算值與實(shí)測(cè)值進(jìn)行比較，結(jié)果如圖2所示。從圖2可以看出，點(diǎn)和點(diǎn)溫度的計(jì)算值較實(shí)測(cè)值均有不同程度的滯后，在升溫階段，滯后性更明顯，而在降溫階段滯后性較好，這主要是由于模型計(jì)算時(shí)，沒(méi)有考慮路基內(nèi)的蒸發(fā)換熱過(guò)程。模型溫度計(jì)算值與實(shí)測(cè)值總體變化趨勢(shì)一致，尤其是下層點(diǎn)的計(jì)算值與實(shí)測(cè)值很吻合，證明參數(shù)的選取和模型的建立是正確的。

(a) 中面層上部A點(diǎn)溫度；(b) 中面層下部B點(diǎn)溫度

2 路基溫度場(chǎng)熱邊界厚度的確定

假設(shè)路基為半無(wú)限體，在周期性變化的熱邊界條件作用下，路基溫度場(chǎng)的解析解為[26]

式中：為路基熱傳導(dǎo)衰減指數(shù)。

由于路基的熱傳導(dǎo)實(shí)質(zhì)是路基面層、基層和土基等各層組成的層狀結(jié)構(gòu)的溫度場(chǎng)問(wèn)題，則溫度振幅經(jīng)過(guò)多層路基材料的總衰減度可以表示為各層路基結(jié)構(gòu)衰減度疊加的形式：

根據(jù)熱?質(zhì)交換特點(diǎn)和溫度變化連續(xù)性，路基表面至淺層地表劃分為對(duì)流換熱區(qū)、輻射換熱區(qū)和熱傳導(dǎo)換熱區(qū)，其中，輻射換熱區(qū)分為上附面層(與空氣銜接)和下附面層(與土介質(zhì)銜接)[16]。在路基某一深度處，大氣溫度和太陽(yáng)輻射強(qiáng)度的日波動(dòng)對(duì)地溫振幅的影響很小，可以忽略，則定義該深度為附面層底部，其厚度為f，以此作為熱邊界厚度。熱邊界底部溫度變化可以表示為

式中：T為熱邊界底部溫度相對(duì)大氣平均溫度的增值；f為熱邊界溫度的振幅；為初相。

考慮到熱邊界底部，太陽(yáng)輻射對(duì)地溫的影響較小，本文采用文獻(xiàn)[27]中標(biāo)準(zhǔn)對(duì)熱邊界的位置進(jìn)行確定：

通過(guò)式(14)確定熱邊界底部的熱傳導(dǎo)衰減指數(shù)f，然后假定熱邊界底部位于多層路面結(jié)構(gòu)體系的第層，并利用式(15)確定第層的位置：

當(dāng)熱邊界位于路面結(jié)構(gòu)體系的第層時(shí)，根據(jù)路基結(jié)構(gòu)中熱傳導(dǎo)的衰減，采用如下方法確定熱邊界厚度：

熱邊界厚度因路基填料不同而異，瀝青面層的等效熱擴(kuò)散率為4.3×10?7~8.7×10?7m2/s[28]，分別計(jì)算不同等效熱擴(kuò)散率瀝青路面的熱邊界厚度，如圖3所示。由圖3可以看出：熱邊界厚度與等效熱擴(kuò)散率具有很好的正相關(guān)性，并且熱邊界厚度的理論解析值和數(shù)值計(jì)算值吻合較好，說(shuō)明通過(guò)式(16)計(jì)算熱邊界厚度是可行的。

圖3 瀝青路面等效熱擴(kuò)散率與熱邊界厚度hf之間的關(guān)系

Fig. 3 Relationship between equivalence thermal diffusivity and thermal boundary layer thickness

3 路面溫度場(chǎng)影響因素敏感性分析

根據(jù)前述分析，給出了確定路基熱邊界位置的方法，而輻射強(qiáng)度、大氣溫度、風(fēng)速等外部氣象因素和路面吸收率、反射率、導(dǎo)熱系數(shù)、比熱容等路面結(jié)構(gòu)熱物理參數(shù)會(huì)對(duì)熱邊界底部溫度增量T產(chǎn)生影響。采用灰關(guān)聯(lián)分析法，以灰關(guān)聯(lián)度作為衡量標(biāo)準(zhǔn)，分析影響熱邊界底部溫度增量的各因素的敏感性。

3.1 灰關(guān)聯(lián)分析原理

以熱邊界底部溫度增量的各影響因素為比較列，記作

以各影響因素相應(yīng)的溫度增量為參考列，記作

其中：與列每項(xiàng)各有若干取值；X[x1，x2，…，x]；Y[y1，y2，…，y]。采用初值化、均值化、區(qū)間相對(duì)值化和歸一化等方法對(duì)X和Y的數(shù)值進(jìn)行歸一化處理：

采用下式得到差異序列矩陣，在差異序列矩陣中找出最大值max和最小值min：

關(guān)聯(lián)系數(shù)V表示比較因素和參考因素的相關(guān)性：

式中：為分辨系數(shù)，用來(lái)提高關(guān)聯(lián)系數(shù)間差異的顯著性，∈[0,1]，一般取0.5。

對(duì)關(guān)聯(lián)系數(shù)求平均值得到關(guān)聯(lián)度A，影響因素的關(guān)聯(lián)度越大，說(shuō)明該因素對(duì)目標(biāo)影響越大，其敏感性越大。

3.2 各因素敏感性分析

影響熱邊界底部溫度增量的各因素的合理取值如表3所示。分別計(jì)算路面吸收率、發(fā)射率、太陽(yáng)輻射強(qiáng)度、對(duì)流換熱系數(shù)和等效熱擴(kuò)散率對(duì)熱邊界底部溫度增量T的影響，計(jì)算結(jié)果如圖4所示。由圖4可以看出：吸收率、對(duì)流換熱系數(shù)和太陽(yáng)輻射強(qiáng)度對(duì)熱邊界底部溫度增量的影響較大。吸收率與溫度增量成線性正相關(guān)，這與吸收率本身的物理含義一致，當(dāng)吸收率從0.50變到0.91時(shí)，溫度增量變化將近1.5 ℃。發(fā)射率與溫度增量成線性負(fù)相關(guān)，這是因?yàn)榘l(fā)射率越大，路面向外的輻射能力越強(qiáng)，釋放的熱能就越多。當(dāng)發(fā)射率在合理范圍內(nèi)變化時(shí)，溫度增量減小0.355 ℃，這比吸收率對(duì)溫度增量的影響要小很多.這是因?yàn)槁访鎸?duì)太陽(yáng)輻射的吸收主要是對(duì)短波輻射的吸收，而路面本身對(duì)外輻射屬于長(zhǎng)波輻射，在同等條件下，短波輻射能量大于長(zhǎng)波輻射能量。太陽(yáng)輻射強(qiáng)度與溫度增量的關(guān)系成線性正相關(guān)，當(dāng)該因素在合理范圍內(nèi)變化時(shí)，溫度增量變化很明顯，這是由于太陽(yáng)輻射強(qiáng)度直接決定進(jìn)入路基的熱量，也直接影響路基溫度場(chǎng)。對(duì)流換熱系數(shù)與溫度增量呈冪函數(shù)關(guān)系，總體來(lái)說(shuō)，溫度增量隨換熱系數(shù)的增大逐漸減小，遞減速率由大變小，在導(dǎo)熱系數(shù)合理取值范圍內(nèi)溫度增量減小3.457 ℃。對(duì)流換熱系數(shù)越大說(shuō)明單位面積路面與環(huán)境交換熱能的能力越強(qiáng)，而這種換熱能力又受風(fēng)速和路面粗糙度等的影響，風(fēng)速會(huì)受近地表氣壓影響，隨著換熱的進(jìn)行，近地表氣壓發(fā)生變化，換熱系數(shù)也會(huì)相應(yīng)變化。等效熱擴(kuò)散率是對(duì)路面結(jié)構(gòu)的綜合反映，溫度增量隨等效熱擴(kuò)散率增大而線性減小，取值在合理范圍內(nèi)變化時(shí)，溫度減小0.384 ℃。對(duì)圖4中給出的不同影響因素對(duì)應(yīng)的熱邊界底部溫度增量T進(jìn)行統(tǒng)計(jì)回歸，得到的回歸方程見(jiàn)表4。

表3 各影響因素取值范圍

1—發(fā)射率；2—等效熱擴(kuò)散率；3—太陽(yáng)有效輻射強(qiáng)度；4—吸收率；5—對(duì)流換熱系數(shù)。

根據(jù)敏感性分析方法，將不同的路面吸收率、發(fā)射率、太陽(yáng)輻射強(qiáng)度、對(duì)流換熱系數(shù)和等效熱擴(kuò)散率以及計(jì)算得到的熱邊界底部溫度增量T代入式(17)和(18)中，得

表4 熱邊界底部溫度增量與各影響因素回歸方程

式中：1，2，3，4和5分別為路面吸收率、路面發(fā)射率、太陽(yáng)輻射強(qiáng)度、對(duì)流換熱系數(shù)和等效熱擴(kuò)散率，1，2，3，4和5為該因素對(duì)應(yīng)的熱邊界底部溫度增量值。

對(duì)式(23)和(24)進(jìn)行矩陣歸一化、差異矩陣計(jì)算等，最終求得關(guān)聯(lián)度：

由關(guān)聯(lián)度可以看出，影響熱邊界底部溫度增量的5個(gè)因素的敏感性由大到小依次為路面吸收率、太陽(yáng)輻射強(qiáng)度、對(duì)流換熱系數(shù)、等效熱擴(kuò)散率和發(fā)射率。

4 溫度增量的二次正交回歸分析

為建立有效預(yù)測(cè)熱邊界條件的模型，采用二次回歸正交分析法對(duì)熱邊界底部溫度增量進(jìn)行研究，以此建立有效預(yù)測(cè)熱邊界條件的模型。

4.1 因素水平編碼

溫度增量在路面吸收率1、發(fā)射率2、太陽(yáng)輻射強(qiáng)度3、對(duì)流換熱系數(shù)4和等效熱擴(kuò)散率5交互作用下變化規(guī)律的研究采用二次正交回歸分析方法，設(shè)置因素水平，見(jiàn)表5。

表5 因素水平表

表5中x代表試驗(yàn)中各因素的實(shí)際取值，x(=1,2,3,4,5)可按式(26)換算成與其相應(yīng)的編碼值z：

式中：[1i,2i]為第個(gè)因素的變化范圍；0i為第個(gè)因素的上下水平均值；Δ為第個(gè)因素的變化間距；為星號(hào)臂長(zhǎng)度，依據(jù)二次正交設(shè)計(jì)確定，文中取1.547。

4.2 回歸正交試驗(yàn)結(jié)果

根據(jù)二次回歸正交設(shè)計(jì)原理[35]，采用L16(215)正交表進(jìn)行正交回歸分析，因素水平編碼值的回歸方程如式(27)，二次回歸方差分析表見(jiàn)表6。

式中：1, 2, …,m，為因素?cái)?shù)；1, 2,…,s；s為樣本總數(shù)。

表6 方差分析表 I

注：0.05(1,6)=5.99；0.01(1,6)=13.7；**表示回歸系數(shù)達(dá)到非常顯著水平；*表示回歸系數(shù)達(dá)到顯著水平；表中省去了未達(dá)到顯著水平的方差分析值。

由表7可知：0.05(7,19)=2.54＜216.192，回歸方程達(dá)到顯著水平，方程有意義。由表7可以看出各回歸系數(shù)均達(dá)到顯著水平，為不顯著因素對(duì)其余7個(gè)顯著因素所產(chǎn)生的影響，構(gòu)造新的統(tǒng)計(jì)量：

式中：1和2分別為修正前后因素水平編碼值回歸方程的相關(guān)系數(shù)；1和2分別為修正前后因素水平編碼值回歸方程的自變量個(gè)數(shù)；值服從自由度為1?2和s?1?1個(gè)自由度的分布。根據(jù)前文中的計(jì)算值得=0.38<0.05(13,5)=3.03，則在給定的顯著水平=0.05下，所刪除的未達(dá)到顯著性的因素對(duì)方程不造成顯著性影響。

表7 方差分析表 II

注：0.05(1,19)=4.38；0.01(1,19)=8.18；**表示回歸系數(shù)達(dá)到非常顯著水平；*表示回歸系數(shù)達(dá)到顯著水平。

由式(29)可得因素實(shí)際值溫度增量方程：

式中：i(=1，2，3，4，5)分別代表實(shí)際路面吸收率、路面發(fā)射率、太陽(yáng)輻射強(qiáng)度、對(duì)流換熱系數(shù)和等效熱擴(kuò)散率的取值。

5 熱邊界條件預(yù)測(cè)模型

5.1 熱邊界條件預(yù)測(cè)模型

以式(13)作為路基溫度場(chǎng)長(zhǎng)期模擬的熱邊界溫度，考慮到熱邊界底部溫度振幅可以忽略，另外，大氣平均氣溫為年波動(dòng)函數(shù)，一般以周期為1 a的正弦函數(shù)表示，路基溫度場(chǎng)長(zhǎng)期模擬時(shí)的熱邊界條件預(yù)測(cè)模型為

式中：為坡面系數(shù)，其值可根據(jù)路面實(shí)際情況參考文獻(xiàn)[36]得到；e為太陽(yáng)日輻射總強(qiáng)度；a和b分別為太陽(yáng)輻射強(qiáng)度的年平均值和年振幅，其值按實(shí)測(cè)太陽(yáng)輻射強(qiáng)度擬合；2為太陽(yáng)輻射強(qiáng)度函數(shù)的初相位；4(對(duì)流換熱系數(shù))可通過(guò)現(xiàn)場(chǎng)實(shí)時(shí)監(jiān)測(cè)的風(fēng)速，按式(7)計(jì)算得到。風(fēng)速由于受大氣環(huán)流和下墊面異常的影響較大，其季節(jié)波動(dòng)性很強(qiáng)，區(qū)域差異較大，隨著全球氣候變暖，我國(guó)年平均風(fēng)速呈現(xiàn)不同程度的降低趨 勢(shì)[37?38]，為此采用如下方法處理風(fēng)速：

式中：a為年平均風(fēng)速；為路基溫度場(chǎng)長(zhǎng)期模擬年限；為年平均風(fēng)速每10年的降低速率，其值根據(jù)路基所處風(fēng)速區(qū)采用文獻(xiàn)[38]得到；b為風(fēng)速年振幅；3為風(fēng)速函數(shù)的初相位。

5.2 預(yù)測(cè)模型驗(yàn)證

為了驗(yàn)證預(yù)測(cè)模型的適用性，取青藏高原北麓河地區(qū)路基試驗(yàn)段現(xiàn)場(chǎng)監(jiān)測(cè)溫度和氣象數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證。日平均氣溫、太陽(yáng)有效日輻射及按實(shí)測(cè)風(fēng)速計(jì)算得到熱交換系數(shù)的實(shí)測(cè)值和擬合值分別如圖5~7所示。熱邊界條件預(yù)測(cè)所需參數(shù)見(jiàn)表8。

將各因素取值和氣溫平均值代入式(31)和(32)計(jì)算，可得模型預(yù)測(cè)值，如圖8所示。從圖8可以看出模型預(yù)測(cè)值與實(shí)測(cè)值較吻合，說(shuō)明熱邊界預(yù)測(cè)模型具有適用性。

圖5 太陽(yáng)有效日輻射實(shí)測(cè)值與擬合值

圖6 對(duì)流換熱系數(shù)實(shí)測(cè)值與擬合值

圖7 日平均氣溫實(shí)測(cè)值與擬合值

表8 預(yù)測(cè)模型參數(shù)取值

圖8 熱邊界底部溫度實(shí)測(cè)值與計(jì)算值

6 結(jié)論

1) 路基溫度場(chǎng)的熱傳遞過(guò)程包括輻射、對(duì)流和熱傳導(dǎo)，是外部因素(太陽(yáng)輻射強(qiáng)度、氣溫、對(duì)流換熱系數(shù))和內(nèi)部因素(路面吸收率、發(fā)射率、路面結(jié)構(gòu))綜合作用的結(jié)果，通過(guò)理論分析，給出了路基熱邊界條件的確定方法。熱邊界厚度與路基等效熱擴(kuò)散率線性正相關(guān)。熱邊界底部溫度增量與路面吸收率、太陽(yáng)輻射強(qiáng)度、等效熱擴(kuò)散率線性正相關(guān)；與路面發(fā)射率線性負(fù)相關(guān)；與對(duì)流換熱系數(shù)成負(fù)指數(shù)相關(guān)。

2) 通過(guò)灰關(guān)聯(lián)方法分析熱邊界底部溫度增量各影響因素的敏感性，表明路面吸收率對(duì)熱邊界底部溫度增量的影響最大，太陽(yáng)輻射強(qiáng)度的影響次之，對(duì)流換熱系數(shù)的影響中等，等效熱擴(kuò)散率的影響較小，發(fā)射率的影響最小。減小吸收率可以顯著降低熱邊界底部溫度增量，增大發(fā)射率也可以起到降低溫度增量的作用，但影響小于降低吸收率。

3) 基于二次正交回歸分析理論建立了熱邊界條件與太陽(yáng)輻射強(qiáng)度、氣溫、對(duì)流換熱系數(shù)、路面吸收率、發(fā)射率、等效熱擴(kuò)散率的函數(shù)關(guān)系。通過(guò)現(xiàn)場(chǎng)氣象和地溫監(jiān)測(cè)資料對(duì)該模型適用性進(jìn)行了驗(yàn)證，為預(yù)測(cè)路基熱邊界條件提供了技術(shù)支撐。

[1] 王紹令, 趙林, 李述訓(xùn), 等. 青藏公路多年凍土段瀝青路面熱量平衡及路基穩(wěn)定性研究[J]. 冰川凍土, 2001, 23(2): 111?118. WANG Shaolin, ZHAO Lin, LI Shuxun, et al. Study on thermal balance of asphalt pavement and roadbed stability in permafrost regions of the Qinghai-Tibet highway[J]. Journal of Glaciology and Geocryology, 2001, 23(2): 111?118.

[2] 竇建明, 胡長(zhǎng)順, 何子文, 等. 青藏公路多年凍土段路基病害分布規(guī)律[J]. 冰川凍土, 2002, 24(6): 780?783. DOU Mingjian, HU Changshun, HE Ziwen, et al. Distributing regularities of subgrade diseases in permafrost section of the Qinghai-Tibet highway[J]. Journal of Glaciology and Geocryology, 2002, 24(6): 780?783.

[3] 張魯新, 原思成, 楊永平. 青藏鐵路多年凍土區(qū)路基變形裂縫發(fā)生機(jī)理及其防治[J]. 第四紀(jì)研究, 2003, 23(6): 604?610. ZHANG Luxin, YUAN Sicheng, YANG Yongping. Mechanism and prevention of deformation cracks of embankments in the permafrost region along Qinghai-Xizang railway[J]. Quaternary Sciences, 2003, 23(6): 604?610.

[4] 俞祁浩, 谷偉, 錢(qián)進(jìn), 等. 多年凍土區(qū)高等級(jí)公路建設(shè)面臨問(wèn)題分析[J]. 公路, 2010, 11(11): 74?80. YU Qihao, GU Wei, QIAN Jin, et al. Problem analysis of high grade highway construction in permafrost regions[J]. Highway, 2010 , 11(11): 74?80.

[5] ZHANG Mingyi, LAI Yuanming, GAO Zhihua, et al. Influence of boundary conditions on the cooling effect of crushed-rock embankment in permafrost regions of Qinghai-Tibetan Plateau[J]. Cold Regions Science and Technology, 2006, 44(3): 225?239.

[6] ZHANG T, OSTERKAMP T E, STAMNES K. Influence of the depth hoar layer of the seasonal snow cover on the ground thermal regime[J]. Water Resources Research, 1996, 32(7): 2075?2086.

[7] ZHANG T, STAMNES K. Impact of climatic factors on the active layer and permafrost at Barrow, Alaska[J]. Permafrost and Periglacial Processes, 1998, 9(3): 229?246.

[8] 穆彥虎. 青藏鐵路凍土區(qū)路基溫度和變形動(dòng)態(tài)變化過(guò)程研究[D]. 蘭州: 中國(guó)科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所, 2012: 24?34. MU Yanhu. Analyses on dynamic variations of embankment thermal regime and deformation along the Qinghai-Tibet Railway in permafrost regions[D]. Lanzhou: The Chinese Academy of Sciences. Cold and Arid Regions Environmental and Engineering Research Institute, 2012: 24?34.

[9] 王可麗, 程國(guó)棟.青藏鐵路沿線地表和路基表面熱力學(xué)模式(Ⅰ): 物理過(guò)程與實(shí)驗(yàn)方案[J]. 冰川凍土, 2002, 24(6): 759?764. WANG Keli, CHENG Guodong. Thermodynamic model of the ground surface and the roadbed surface along the Qinghai-Tibetan Railway (Ⅰ): physical process and experimental scheme[J]. Journal of Glaciology and Geocryology, 2002, 24(6): 759?764.

[10] 王可麗, 程國(guó)棟, 江灝, 等. 青藏鐵路沿線地表和路基表面熱力學(xué)模式(Ⅱ): 無(wú)云大氣條件下模擬試驗(yàn)結(jié)果分析[J]. 冰川凍土, 2004, 26(2): 171?176. WANG Keli, CHENG Guodong, JIANG Hao, et al. Thermodynamic model of the ground surface and embankment surface along the Qinghai—Tibet Railway (Ⅱ): results in the cloud-free condition[J]. Journal of Glaciology and Geocryology, 2004, 26(2): 171?176.

[11] 江灝, 吳青柏, 王可麗, 等. 青藏鐵路沿線地表和路基表面熱力學(xué)模式(Ⅲ): 參數(shù)化方案[J]. 冰川凍土, 2005, 27(5): 680?685. JIANG Hao, WU Qingbai, WANG Keli, et al. Thermodynamic model of the ground surface and embankment surface along the Qinghai—Tibet Railway (Ⅲ): parameterization schemes[J]. Journal of Glaciology and Geocryology, 2005, 27(5): 680?685.

[12] LIU Zhiyun, CHEN Jianbing, JIN Long, et al. Roadbed temperature study based on earth-atmosphere coupled system in permafrost regions of the Qinghai—Tibet Plateau[J]. Cold Regions Science and Technology, 2013, 86(2): 167?176.

[13] 張文綱, 李述訓(xùn), 吳通華, 等. 青藏高原地氣溫差變化分析[J]. 地理學(xué)報(bào), 2006, 61(9): 899?910. ZHANG Wengang, LI Shuxun, WU Tonghua, et al. Changes of the differences between ground and air temperature over the Qinghai—Xizang Plateau[J]. Acta Geographica Sinica, 2006, 61(9): 899?910.

[14] 李述訓(xùn), 吳通華. 青藏高原地氣溫度之間的關(guān)系[J]. 冰川凍土, 2005, 27(5): 627?632. LI Shuxun, WU Tonghua. The relationship between air temperature and ground temperature in the Tibetan Plateau[J]. Journal of Glaciology and Geocryology, 2005, 27(5): 627?632.

[15] KARUNARATNE K C, BURN C R. Relations between air and surface temperature in discontinuous permafrost terrain near Mayo, Yukon Territory [J]. Canadian Journal of Earth Sciences, 2004, 41(12): 1437?1451.

[16] 朱林楠. 高原凍土區(qū)不同下墊面的附面層研究[J].冰川凍土, 1988, 10(1): 8?14. ZHU Linnan. Study of the adherent layer on different types of ground in permafrost regions on the Qinghai-Xizang Plateau[J]. Journal of Glaciology and Geocryology, 1988, 10(1): 8?14.

[17] ZHANG Mingyi, LAI Yuanming, DONG Yuanhong. Numerical study on temperature characteristics of expressway embankment with crushed-rock revetment and ventilated ducts in warm permafrost regions[J]. Cold Regions Science & Technology, 2009, 59(1): 19?24.

[18] 趙世運(yùn). 青藏鐵路凍土區(qū)“冷卻地基”設(shè)計(jì)理論的研究和實(shí) 踐[D]. 蘭州: 中國(guó)科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所, 2007: 173?209. ZHAO Shiyun. The study and practice of cooling embankment in permafrost regions of Qinghai-Tibet Railway[D]. Lanzhou: The Chinese Academy of Sciences. Cold and Arid Regions Environmental and Engineering Research Institute, 2007: 173?209.

[19] 孫志忠. 青藏鐵路多年凍土區(qū)塊石護(hù)坡路基試驗(yàn)研究[D]. 蘭州: 中國(guó)科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所, 2006: 90?105. SUN Zhizhong. Study on the crushed-rock revetment embankment in permafrost regions of Qinghai—Tibet Railway[D]. Lanzhou: The Chinese Academy of Sciences. Cold and Arid Regions Environmental and Engineering Research Institute, 2006: 90?105.

[20] ASAEDA T, CA V T, WAKE A. Heat storage of pavement and its effect on the lower atmosphere[J]. Atmospheric Environment, 1996, 30(3): 413?427.

[21] 張中瓊. 多年凍土區(qū)瀝青路面熱效應(yīng)機(jī)理研究[D]. 蘭州: 中國(guó)科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所, 2012: 25?30. ZHANG Zhongqiong. Study on the mechanism of asphalt pavement’s thermal effects in permafrost regions[D]. Lanzhou: The Chinese Academy of Sciences. Cold and Arid Regions Environmental and Engineering Research Institute, 2012: 25?30.

[22] 嚴(yán)作人. 層狀路面體系的溫度場(chǎng)分析[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 1984, 16(3): 79?88. YAN Zuoren. Analysis of the temperature field in layered pavement system[J]. Journal of Tongji University(Natural Science), 1984, 16(3): 79?88.

[23] 王孫富. 瀝青路面結(jié)構(gòu)溫度場(chǎng)與溫度應(yīng)力的數(shù)值模擬分析[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué)交通科學(xué)與工程學(xué)院, 2010: 15?16. WANG Sunfu. Simulation analysis of temperature field and thermal stress for bituminous pavement structure[D]. Harbin: Harbin Institute of Technology. School of Transportation Science and Engineering, 2010: 15?16.

[24] 張朝暉. ANSYS熱分析教程與實(shí)例解析[M]. 北京: 中國(guó)鐵道出版社, 2007: 51?52. ZHANG Chaohui. ANSYS thermal analysis tutorial and case analysis[M]. Beijing: China Railway Press, 2007: 51?52.

[25] 馮德成, 胡偉超, 于飛, 等. 瀝青路面材料熱物性參數(shù)對(duì)溫度場(chǎng)的影響及敏感性分析[J]. 公路交通科技, 2011, 28(11): 12?19. FENG Decheng, HU Weichao, YU Fei, et al. Impact of asphalt pavement thermophysical property on temperature field and sensitivity analysis[J]. Journal of Highway and Transportation Research and Development, 2011, 28(11): 12?19.

[26] 趙鎮(zhèn)南. 傳熱學(xué)[M]. 2版. 北京: 高等教育出版社, 2007: 142?145. ZHAO Zhennan. Heat transfer[M]. 2nd ed. Beijing: Higher Education Press, 2007: 142?145.

[27] 白青波, 李旭, 田亞護(hù). 路基溫度場(chǎng)長(zhǎng)期模擬中的地表熱邊界條件研究[J]. 巖土工程學(xué)報(bào), 2015, 37(6): 1142?1149. BAI Qinbo, LI Xu, TIAN Yahu. Upper boundary conditions in long-term thermal simulation of subgrade[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(6): 1142?1149.

[28] 張慧彧, 鄒玲, 紀(jì)小平. 瀝青混合料導(dǎo)熱系數(shù)的試驗(yàn)研究[J]. 公路, 2011(10): 50?51. ZHANG Huiyu, ZOU Lin, JI Xiaoping. Experimental study on the heat conduction of asphalt mixtures[J]. Highway, 2011(10): 50?51.

[29] 陳明宇. 瀝青路面太陽(yáng)能集熱性能研究[D]. 武漢: 武漢理工大學(xué)材料科學(xué)與工程學(xué)院, 2010: 22?23. CHEN Mingyu. Research on solar energy collection using asphalt pavement[D]. Wuhan: Wuhan University of Technology. School of Materials Science and Engineering, 2010: 22?23.

[30] 張?chǎng)? 瀝青路面熱反射與熱阻技術(shù)降溫機(jī)理與應(yīng)用研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué)交通科學(xué)與工程學(xué)院, 2011: 38?40. ZHANG Xin. Research on heat reflection and thermal resistance technology of asphalt pavement cooling mechanism and its application[D]. Harbin: Harbin Institute of Technology. School of Transportation Science and Engineering, 2011: 38?40.

[31] 翁篤鳴. 中國(guó)太陽(yáng)直接輻射的氣候計(jì)算及其分布特征[J]. 太陽(yáng)能學(xué)報(bào), 1986, 7(2): 3?12. WENG Duming. Climatical method for direct solar radiation calculation and its distribution over China[J]. Acta Energiae Solaris Sinica, 1986, 7(2): 121?130.

[32] 馬有哲, 劉小寧, 許松. 中國(guó)太陽(yáng)輻射數(shù)據(jù)集及其質(zhì)量檢驗(yàn)分析[J]. 氣象科技, 1998(2): 53?56. MA Youzhe, LIU Xiaoning, XU Song. Analysis of data set and quality inspection radiation in China[J]. Meteorological Science and Technology, 1998(2): 53?56.

[33] 吳其重, 王自發(fā), 崔應(yīng)杰. 我國(guó)近20年太陽(yáng)輻射時(shí)空分布狀況模式評(píng)估[J]. 應(yīng)用氣象學(xué)報(bào), 2010, 21(3): 343?351. WU Qizhong, WANG Zifa, CUI Yingjie. Evaluating the solar radiation resources of China in recent 20 years by meteorological model[J]. Journal of Applied Meteorological Science, 2010, 21(3): 343?351.

[34] SHARPLES S. Full sealed measurements of convective energy losses from exterior building surfaces[J]. Building Environment, 1984, 19: 31?38.

[35] WILLIAM D B. 理解回歸假設(shè)[M]. 余珊珊, 譯. 上海格致出版社, 上海出版社, 2012: 20?50. WILLIAM D B. Understanding regression assumptions[M]. YU Shanshan, trans. Shanghai: Truth and Wisdom Press, Shanghai People’s Publishing House, 2012: 20?50.

[36] 王鐵行. 多年凍土地區(qū)路基計(jì)算原理及臨界高度研究[D]. 西安: 長(zhǎng)安大學(xué)公路工程學(xué)院, 2001: 19?21. WANG Tiexing. Research on calculation principle and critical height of subgrade in permafrost regions[D]. Xi’an: Chang’ an University. School of Highway Engineering, 2001: 19?21.

[37] 陳練. 氣候變暖背景下中國(guó)風(fēng)速(能)變化及其影響因子研究[D]. 南京: 南京信息工程大學(xué)大氣科學(xué)學(xué)院, 2013: 14?15.CHEN Lian. Changes and their impact factors of wind speed (energy) over china under the background of climate warming[D]. Nanjing: Nanjing University of Information Science and Technology. College of Atmospheric Science, 2013: 14?15.

[38] 江瀅, 羅勇, 趙宗慈, 等. 近50年中國(guó)風(fēng)速變化及原因分析[C]// 中國(guó)氣象學(xué)會(huì)2007年年會(huì)氣候變化分會(huì)場(chǎng)論文集. 北京, 2007: 80?89. JIANG Ying, LUO Yong, ZHAO Zongci, et al. Analysis of wind speed changes and causes in China in the past 50 years[C]// Proceedings of the 2007 Annual Conference of the Chinese Meteorological Society on Climate Change. Beijing, 2007: 80?89.

A prediction model on thermal boundary condition of embankment based on boundary layer theory

LUO Xiaoxiao1, 2, YU Qihao1, MA Qinguo3, ZHANG Jianming1

(1. State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China;2. School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China;3. State Key Laboratory of Subtropical Building Science, South China University of Technology, Guangzhou 510641, China)

Based on the boundary-layer theory, a model for heat exchange of embankment structure was established firstly considering all influence factors including meteorological conditions and thermal physical parameters of pavement structure to determine thermal boundary condition. And then the method to determine the thermal boundary condition was proposed. Secondly, sensitivity analysis of the influencing factors which affect the temperature increment at the thermal boundary was implemented using the method of grey relation. Finally, considering the influences of all factors and their interaction, the quadratic regression orthogonality analysis method was used to study the temperature increment at the thermal boundary. A prediction model on thermal boundary condition was established and verified by the measured temperature of Qinghai—Tibet Highway. The results show that the thickness of thermal boundary has a positive relationship with the equivalence thermal diffusivity of embankment. The road surface absorptivity has the largest influence on the thermal boundary, followed by the solar radiation, with the convective heat transfer coefficient in the middle, the equivalent thermal diffusivity being the next factor, and the road surface emissivity influences the smallest. The model is reliable and can be used to predict the thermal boundary of embankment.

boundary layer; thermal boundary condition; sensitivity analysis; orthogonality analysis method; prediction model

U416.1+68

A

1672?7207(2019)03?0658?11

10.11817/j.issn.1672-7207.2019.03.020

2018?04?09；

2018?08?11

國(guó)家科技支撐計(jì)劃項(xiàng)目(2014BAG05B03)；國(guó)家自然科學(xué)基金資助項(xiàng)目(41401088) (Project(2014BAG05B03) supported by the National Science and Technology Pillar Program of China; Project(41401088) supported by the National Natural Science Foundation of China)

俞祁浩，博士，研究員，從事凍土路基穩(wěn)定性研究；E-mail：yuqh@lzb.ac.cn

(編輯 趙俊)