蔡蕓 姚木申 盧寶 趙蘇民

摘要 針對不同基巖地層條件，開展地埋管鉆孔的現(xiàn)場熱響應(yīng)試驗與地下?lián)Q熱實驗，分析相關(guān)影響因素。結(jié)果表明，在目前地埋管鉆孔深度范圍內(nèi)，初始平均地溫均高于年平均氣溫，且二者呈正相關(guān)關(guān)系;就基巖地區(qū)熱響應(yīng)試驗而言，巖性差異在很大程度上決定了測試結(jié)果，其中巖石熱導(dǎo)率和熱擴散率大小順序為沉積巖（灰?guī)r）> 變質(zhì)巖 >火山巖（安山巖、玄武巖）。地埋管換熱量不僅受地層巖性、初始地溫、管內(nèi)流體溫度等影響，而且還受地下水徑流作用影響。地下水流動性越強，對地埋管換熱的強化效應(yīng)也越強。

關(guān) 鍵 詞 淺層地?zé)崮?熱響應(yīng)試驗;熱物性;換熱量

中圖分類號 TK529? ? ?文獻(xiàn)標(biāo)志碼 A

Abstract Thermal response tests （TRTs） and heat transfer tests of ground heat exchangers （GHEs） are performed under different rocky stratigraphic conditions， and the influencing factors are analyzed. Results show that within the current depth range of GHEs， the initial average ground temperature is always higher than the annual average ambient temperature， and there is a positive relationship between the two temperatures. The results of TRTs are obviously affected by the lithologic differences， and the thermal conductivity together with the thermal diffusivity of rocks in the work are as follows： sedimentary rock （limestone） > metamorphic rock > volcanic rock （andesite and basalt）. The heat transfer rate of GHEs is affected not only by stratigraphic and lithologic conditions， the initial ground temperature and the fluid temperature， but also by the groundwater flow conditions. The more rapid the groundwater flows， the stronger the enhanced effect on heat transfer becomes.

Key words shallow geothermal energy; thermal response test; thermal properties; heat transfer rate

0 引言

巖土體熱物性特征是地埋管地源熱泵系統(tǒng)設(shè)計的主要參考因素之一。在系統(tǒng)設(shè)計時，若對地埋管鉆孔換熱能力不了解，盲目地增加或者減少鉆孔數(shù)量，極易造成機組選型不合理，成本增加，大馬拉小車，或者導(dǎo)致用戶供熱制冷量不足，熱泵機組運行不穩(wěn)定，從而影響系統(tǒng)的使用壽命。巖土體的熱物性特征受氣候條件、地質(zhì)構(gòu)造、巖性、水文地質(zhì)條件等綜合影響，不同條件下表現(xiàn)特征各異。目前，現(xiàn)場熱響應(yīng)試驗是獲得巖土體熱物性特征及地埋管換熱能力的主要技術(shù)手段，前人在此方面開展了較多理論與實驗工作[1-2]，但主要集中在第四系發(fā)育較厚的平原地區(qū)，較少涉及基巖地層。周世玉等[3]分析了重慶市基巖地質(zhì)條件下測試時間對巖土熱物性結(jié)果的影響，建議測試時間宜大于60 h（單U）和50 h（雙U）。王松濤[4]等在山東威海環(huán)翠區(qū)開展了地埋管換熱試驗，對各類型地埋管換熱量進(jìn)行了對比分析，其中地下12 m以淺為第四系松散沉積層，以下為元古代巖漿巖地層（以花崗巖為主）?；谏鲜霰尘?，本文更廣泛地在沉積巖、火山巖和變質(zhì)巖等基巖地層上，開展地埋管換熱孔的現(xiàn)場熱響應(yīng)試驗和換熱特性實驗，對不同類型巖土體熱物性參數(shù)進(jìn)行計算對比，并分析其主要影響因素，旨在對含有基巖地層條件下的地埋管地源熱泵系統(tǒng)設(shè)計與優(yōu)化提供一定的參考依據(jù)。

1 試驗概況

1.1 試驗鉆孔及地埋管

選取位于四川理塘、遼寧錦州、山東濟南和山東濰坊的4眼地埋管換熱鉆孔作為研究對象，分別命名為K1（100°17′E，29°58′N）、K2（121°04′″E，41°12′4N）、K3（117°10′E，36°38′N）和K4（119°05′E，36°37′N），海拔高度分別為4 080 m、180 m、159 m和51 m。

表1給出了試驗孔地層巖性特征及其地埋管的基本情況[5-8]。其中，K1地處甘孜藏族自治州西南部、橫斷山脈北段，試驗孔巖性上部為上更新統(tǒng)冰水堆積層（以砂、沙礫、泥礫和黏土為主），其下為三疊系上統(tǒng)侏倭組基巖（石英砂巖與板巖不等厚互層）。K2地處錦州市小凌河下游的錦凌水庫，屬于山地丘陵區(qū)，試驗孔巖性上部以第四系粉質(zhì)黏土為主，下部為義山組安山巖。K3地處濟南市歷城區(qū)，屬于山前沖積平原區(qū)，試驗孔巖性上部為第四系松散堆積層，下部為灰?guī)r地層。K4地處濰西南中低山丘陵區(qū)，試驗巖性上部多為第四系粉質(zhì)黏土，下部為玄武巖和砂礫巖。

1.2 熱響應(yīng)試驗

采用FTPT-11型現(xiàn)場熱響應(yīng)測試儀進(jìn)行（圖1）。測試儀主要由熱泵系統(tǒng)、電加熱器、補水箱、循環(huán)泵、循環(huán)管路、溫度和流量檢測元件等組成。其中，地埋管供、回水管及循環(huán)管路均布置有Pt1000型溫度傳感器，精度為±0.1 ℃，地埋管供水管處安裝電磁流量計，精度為±0.01 m3/h，用于測量管內(nèi)循環(huán)流量。地埋管排熱量和取熱量測試采用恒溫法，巖土熱物性測試采用恒熱流法，加熱功率范圍為0～8 kW，可任意調(diào)節(jié)，符合國家標(biāo)準(zhǔn)《地源熱泵系統(tǒng)工程技術(shù)規(guī)范》（GB50366—2009）中的相關(guān)技術(shù)要求。

現(xiàn)場熱響應(yīng)試驗的具體工作流程如下：1）首先進(jìn)行初始地溫空載循環(huán)測試，當(dāng)?shù)芈窆苓M(jìn)出口溫度12 h內(nèi)變化幅度小于0.5 ℃，且溫差小于0.2 ℃時，獲得巖土體初始平均溫度;2）采用恒熱流加熱工況，保證加熱功率和流量波動范圍在±1 %以內(nèi)，每間隔1 min記錄回路中的流量和進(jìn)出口溫度，直至溫度達(dá)到基本穩(wěn)定（持續(xù)變化小于0.5 ℃/d）;3）進(jìn)行夏季排熱工況試驗，排熱溫度設(shè)置在30～33 ℃左右;4）進(jìn)行地溫恢復(fù);5）進(jìn)行冬季取熱工況試驗，取熱溫度設(shè)置在5～7 ℃左右。

1.3 數(shù)據(jù)處理

上述4組熱響應(yīng)試驗的結(jié)果處理均采用了國標(biāo)推薦的線熱源理論模型。該模型中，在持續(xù)的恒定排熱熱流作用下，地埋管的供回水平均溫度隨加熱時間的升溫變化可表述為[2]

式（1）可進(jìn)一步簡化為

式中：[k]為[Tf-ln（t）]曲線的線性段變化斜率;[m]為直線截距;可分別計算如下：

式中：[Tf]為埋管內(nèi)流體平均溫度，取進(jìn)出口溫度的算術(shù)平均值;[Qheat]為加熱功率;[H]為垂直方向上地埋管深度;[λ]為巖土體平均熱導(dǎo)率;[a]為巖土體熱擴散系數(shù);[t]為持續(xù)測試時間;[r]為試驗鉆孔半徑;[γ]為歐拉常數(shù);[Rb]為鉆孔熱阻;[T0]為巖土體遠(yuǎn)處未受擾動的溫度，即初始平均地溫。

巖土體熱擴散系數(shù)[a]計算公式為

式中：[ρ]為巖土體平均密度;[c]為巖土體比熱。

2 結(jié)果與討論

2.1 巖土初始平均溫度

巖土體初始平均溫度（或初始平均地溫）是淺層地?zé)崮苜Y源勘查評價與開發(fā)應(yīng)用的關(guān)鍵參數(shù)之一，往往受地質(zhì)構(gòu)造、地層巖性和層序、地形地貌、地下水、太陽輻射、氣候條件等參數(shù)影響。圖2給出了4眼試驗孔的初始平均地溫與年平均氣溫的變化關(guān)系?？梢钥闯?，初始平均地溫均高于年平均氣溫，且就整體而言，在地埋管鉆孔深度范圍內(nèi)，二者之間呈正相關(guān)變化。Wang等[9]對河北平原區(qū)域內(nèi)的5眼100～120 m深地埋管試驗孔進(jìn)行了初始平均地溫測試，結(jié)果滿足：[T0=0.928 8Ta+2.92]，其中[Ta]為年平均氣溫。但本文測試結(jié)果中，只有K3滿足Wang模型，其主要原因在于該鉆孔地處山前沖積平原區(qū)，與文獻(xiàn)[9]中試驗孔具有類似的地貌特征。K2和K4與Mohamed模型吻合良好[（T0=0.951Ta+4.514][10]，且二者均屬于山地丘陵地貌，這也暗示該模型與Wang模型的適用范圍有所差異。K1既不符合Wang模型，也不符合Mohamed模型，其主要原因在于該鉆孔地處高海拔山區(qū)。綜合對比可知，關(guān)于“初始平均地溫通常比年平均氣溫高2～3 ℃”的這一習(xí)慣說法，實質(zhì)上是對于平原區(qū)而言的，并不適合于山地丘陵區(qū)和高海拔山區(qū)，不能一概而論。例如，對于K3，二者溫度差為1.9 ℃;對于K2和K4，二者溫度差增加至3.9 ℃;對于K1，二者溫度差高達(dá)8.6 ℃。此外，根據(jù)文獻(xiàn)[11]，青島市某100 m深地埋管鉆孔（沿海丘陵區(qū)，地下3.5 m以深為花崗巖地層），[T0]= 16.3 ℃，比[Ta]偏高4.0 ℃，符合Mohamed模型。根據(jù)文獻(xiàn)[5]，山東威海沿海丘陵區(qū)地埋管鉆孔[T0]= 16.3 ± 0.3 ℃，比[Ta]偏高4.2 ℃，符合Mohamed模型。根據(jù)筆者單位以往測試結(jié)果，北京延慶區(qū)某100 m深地埋管鉆孔（115°48.74′E，40°30.06′N），[T0]= 13.0 ℃，與[Ta]= 8.7 ℃相比，偏高4.3 ℃，符合Mohamed模型。

另一方面，初始平均地溫還與地質(zhì)構(gòu)造及其活動關(guān)系密切。例如，文獻(xiàn)[12]結(jié)果表明，雖然同屬于河北平原區(qū)，曹妃甸某地埋管鉆孔[T0]= 14.5 ℃，比[Ta]偏高3.1 ℃，更接近Wang模型，而衡水某地埋管鉆孔，由于地處地溫場異常區(qū)，[T0]= 16.9 ℃，比[Ta]偏高4.3 ℃，更接近Mohamed模型。由此可見，初始平均地溫與年平均氣溫之間的函數(shù)關(guān)系是因地而異的，需要區(qū)別分析。由于高海拔山區(qū)的測試點較少且公開文獻(xiàn)報導(dǎo)欠缺，本文尚無法給出適合的預(yù)測模型，這還有待于以后進(jìn)一步數(shù)據(jù)積累與深入研究。

2.2 巖土熱物性參數(shù)

表2給出了4眼試驗孔的巖土熱物性參數(shù)測試結(jié)果。結(jié)果表明，K3的平均熱導(dǎo)率和熱擴散率最高，大于3 W/（m·K），其次為K1，介于2～3 W/（m·K）之間，而K2和K4相對較低，均小于2 W/（m·K），造成此差異的主要影響因素在于基巖巖性。如前所述，K3為沉積巖（灰?guī)r），K1為變質(zhì)巖，而K2和K4為火山巖（安山巖、玄武巖）。其中，較致密的灰?guī)r熱導(dǎo)率通常在2.5～3.8 W/（m·K），略低于白云巖[13]，而K3熱導(dǎo)率為3.62 W/（m·K），顯然屬于高熱導(dǎo)率灰?guī)r類型。

表3給出了K1鉆孔巖芯實驗室測試結(jié)果。可以看出，在垂直鉆孔方向上，隨著深度增加，巖石孔隙率呈逐漸降低趨勢，同時密度、熱導(dǎo)率和熱擴散率呈逐漸增加趨勢。巖樣的熱導(dǎo)率變化范圍1.56～3.08 W/（m·K），而平均熱導(dǎo)率測試值為2.14 W/（m·K），二者相對一致。在火山巖中，安山巖和玄武巖均屬于低熱導(dǎo)率類型。以玄武巖為例[13]，熱導(dǎo)率多在1.6～2.4 W/（m·K）。特別地，對于多孔狀玄武巖，由于氣孔發(fā)育較多，熱導(dǎo)率甚至可低至1.0 W/（m·K）以下。K4平均熱導(dǎo)率為1.40 W/（m·K），與上述變化范圍下限接近。表4給出了K2鉆孔巖芯實驗室測試結(jié)果。可以看出，安山巖樣的熱導(dǎo)率較低，變化范圍為1.26～1.60 W/（m·K），和K4差別不大。綜上可知，對于基巖地區(qū)地埋管鉆孔熱響應(yīng)試驗而言，巖性差異在很大程度上決定了測試結(jié)果，是影響巖土體熱物性特征的關(guān)鍵因素之一。因此，在開展現(xiàn)場熱響應(yīng)試驗的同時，適當(dāng)采集巖芯樣品進(jìn)行實驗室分析具有重要的實際指導(dǎo)意義，有助于形成更為完整合理的數(shù)據(jù)解譯。

2.3 單位延米換熱量

圖3給出了試驗孔的單位延米換熱量測試結(jié)果，其中K4鉆孔沒有開展換熱實驗，K3鉆孔沒有開展夏季排熱工況實驗?？梢钥闯?，就K1和K2鉆孔而言，二者地溫接近，但后者夏季換熱量略高于前者，冬季排熱量略低于前者。例如，K1和K2的冬季取熱量分別為24.6 W/m和22.3 W/m（供回水溫度分別為5 ℃和8 ℃），夏季排熱量分別為54.7 W/m和59.9 W/m（供回水溫度分別為30 ℃和25 ℃）。從鉆孔巖性上看，K1平均熱導(dǎo)率高于K2，但從水文地質(zhì)條件上看，K1上部第四系厚度更大，且下部巖性為高壓下形成的變質(zhì)砂巖和板巖，透水性較差。相比之下，K2上部第四系厚度僅6 m左右，下部安山巖多為破碎狀，大氣降水入滲補給作用強烈，地下水徑流作用強烈，有利于熱的傳遞和散失，因而相同工況下的夏季排熱量更高一些。對于K3，由于初始地溫（15.7 ℃）明顯高于K1（12.3 ℃）和K2（11.9 ℃），且鉆孔巖性為熱導(dǎo)率很高的灰?guī)r（表2），因此其單位延米換熱量整體上高于K1和K2。例如，K3的冬季取熱量達(dá)到46.6 W/m（供回水溫度分別為5.9 ℃和9 ℃），比K1和K2分別偏高1.9倍和1. 1倍;夏季排熱量的預(yù)測值為66.8 W/m，比K1和K2分別偏高22%和12%。通常而言，地下水流動性越強，對地埋管換熱的強化效應(yīng)越強。王松濤等[4]測試結(jié)果表明，當(dāng)初始地溫16.7 ℃時，DN32雙U形地埋管的夏季排熱量達(dá)到105～128 W/m，冬季取熱量達(dá)到75～93 W/m，與此同時測試熱導(dǎo)率達(dá)到11.5 W/（m·K），而巖芯測試熱導(dǎo)率僅為2. 65 W/（m·K）（以黑云母花崗巖、角閃巖為主），這表明地下水流動不僅會使地埋管單位延米換熱量增加，還會影響到熱響應(yīng)試驗結(jié)果。從此角度判斷，K2熱響應(yīng)試驗結(jié)果為1.67 W/（m·K），而巖芯測試熱導(dǎo)率為1.26～1.60 W/（m·K），因此其地下水流動的強化傳熱效應(yīng)遠(yuǎn)弱于文獻(xiàn)[4]。

3 結(jié)論

1）在目前地埋管鉆孔深度范圍內(nèi)，初始平均地溫均高于年平均氣溫，且二者呈正相關(guān)關(guān)系，其具體函數(shù)關(guān)系因地而異。其中，K2和K4滿足Mohamed模型，K3滿足Wang模型，而K1不滿足現(xiàn)有模型。

2）對于基巖地區(qū)地埋管鉆孔熱響應(yīng)試驗而言，巖性差異在很大程度上決定了測試結(jié)果，是影響巖土體熱物性特征的關(guān)鍵因素之一。適當(dāng)采集巖芯樣品進(jìn)行實驗室分析，有助于形成更為完整合理的數(shù)據(jù)解譯。就本文結(jié)果而言，巖石熱導(dǎo)率和熱擴散率大小順序為沉積巖（灰?guī)r）> 變質(zhì)巖 > 火山巖（安山巖、玄武巖）。

3）地埋管換熱量不僅受地層巖性、初始地溫、管內(nèi)流體溫度等影響，而且還受地下水徑流作用影響。地下水流動性越強，對地埋管換熱的強化效應(yīng)越強，同時還會影響熱響應(yīng)測試結(jié)果。此外，上部蓋層厚度較薄的地區(qū)，巖土體受氣候、水文循環(huán)等因素影響較大一些。

參考文獻(xiàn)：

[1]? ? 崔俊奎，王震，任常在. 不同巖土熱物性參數(shù)測試方法的比對試驗研究[J]. 工程勘察，2013，41（11）：36-39.

[2]? ? WANG H J，LIU L B，QI C Y. Comparisons of test methods to determine the ground thermal conductivity for geothermal applications [J]. Geothermal Resource Council Transactions，2010，34（1）：969-973.

[3]? ? 周世玉，崔文智，廖全，等. 基巖地質(zhì)條件熱響應(yīng)測試時間優(yōu)化[J]. 暖通空調(diào)，2017，47（1）：129-133.

[4]? ? 王松濤，吳振，王華軍，等基巖地區(qū)地埋管換熱器傳熱特性的測試研究[J]. 河北工業(yè)大學(xué)學(xué)報，2012，41（6）：69-73.

[5]? ? 盧寶，蔡蕓. 四川省甘孜州理塘縣城區(qū)地?zé)峁┡椖繜犴憫?yīng)試驗總結(jié)報告[R]. 天津：天津地?zé)峥辈殚_發(fā)設(shè)計院，2017.

[6]? ? 姚木申，盧寶. 遼寧省錦州市錦凌水庫項目辦公樓地源熱泵工程現(xiàn)場熱響應(yīng)試驗總結(jié)報告[R]. 天津：天津地?zé)峥辈殚_發(fā)設(shè)計院，2017.

[7]? ? 姚木申，楊玉新. 濟南市淺層地?zé)崮苷{(diào)查評價熱響應(yīng)試驗總結(jié)報告[R]. 天津：天津地?zé)峥辈殚_發(fā)設(shè)計院，2012.

[8]? ? 蔡蕓，姚木申. 山東省淺層地?zé)崮苷{(diào)查評價（濰坊）現(xiàn)場熱響應(yīng)試驗報告[R]. 天津：天津地?zé)峥辈殚_發(fā)設(shè)計院，2016.

[9]? ? WANG H J，F(xiàn)ENG L F，WANG Y H，et al. Test investigations of recovery process of the ground temperature of boreholes after being disturbed[J]. Electronic Journal of Geotechnical Engineering，2015，20（10）：4249-4262.

[10]? OUZZANE M，ESLAMI-NEJAD P，BADACHE M，et al. New correlations for the prediction of the undisturbed ground temperature[J]. Geothermics，2015，53：379-384.

[11]? 杜付昌，饒慧鋒，賀富嬌，等. 青島旅游學(xué)校地源熱泵系統(tǒng)工程的設(shè)計與實施[C]//第5屆中國地源熱泵行業(yè)高層論壇論文集，天津，2013.

[12]? 王華軍，齊承英，王華軍，等. 地下熱響應(yīng)實驗中土壤初始溫度的探討[J]. 暖通空調(diào)，2010，40（1）：95-98.

[13]? ROBERTSON E C. Thermal conductivity of rocks [R]. US Geological Survey，1979.

[責(zé)任編輯 田 豐]