青島藍(lán)色硅谷土壤熱物性測(cè)試與分析*

朱川生，葉燦滔，張冬冬，劉廣平，劉 茜，龔宇烈?

（1. 中國(guó)科學(xué)院廣州能源研究所，廣州 510640；2. 中國(guó)科學(xué)院可再生能源重點(diǎn)實(shí)驗(yàn)室，廣州 510640；3. 廣東省新能源和可再生能源研究開(kāi)發(fā)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，廣州 510640；4. 中國(guó)科學(xué)院大學(xué)，北京 100049；5. 即墨市熱電廠，山東 青島 266200）

0 引 言

我國(guó)山東省即墨市藍(lán)色硅谷區(qū)域由于特殊的地質(zhì)構(gòu)造等原因，地?zé)豳Y源較為豐富[1]。然而，由于對(duì)該區(qū)域地?zé)豳Y源的勘查程度較低，缺乏對(duì)地?zé)豳Y源的綜合評(píng)價(jià)及開(kāi)發(fā)利用規(guī)劃等原因，目前當(dāng)?shù)氐牡責(zé)豳Y源利用形式單一，優(yōu)質(zhì)地?zé)豳Y源浪費(fèi)嚴(yán)重。因此，查明該區(qū)域現(xiàn)階段可供開(kāi)發(fā)的地?zé)豳Y源量及其分布情況十分重要。

土壤源熱泵技術(shù)利用地球表面淺層地?zé)豳Y源作為冷熱源進(jìn)行能量轉(zhuǎn)換，其不但能實(shí)現(xiàn)節(jié)能和環(huán)保，而且能克服風(fēng)冷熱泵隨環(huán)境溫度變化帶來(lái)性能不穩(wěn)定的缺陷，確保建筑設(shè)備獲得穩(wěn)定和可靠的運(yùn)行負(fù)荷，提高系統(tǒng)運(yùn)行效率。目前，歐美國(guó)家正大力發(fā)展和推廣應(yīng)用土壤源熱泵技術(shù)，由于系統(tǒng)采用地下埋管換熱器，使得土壤源熱泵在設(shè)計(jì)上比空氣源和水源熱泵系統(tǒng)困難很多，巖土熱物性參數(shù)的對(duì)鉆孔的數(shù)量及深度具有顯著的影響[2]。若巖土熱物性參數(shù)測(cè)試不準(zhǔn)確，將導(dǎo)致設(shè)計(jì)負(fù)荷和實(shí)際負(fù)荷不相匹配，使埋管換熱器的準(zhǔn)確設(shè)計(jì)受影響，因此巖土熱物性是地埋管換熱系統(tǒng)設(shè)計(jì)的重要參數(shù)。當(dāng)加熱功率測(cè)量值誤差為0.1 kW時(shí)，巖土導(dǎo)熱系數(shù)計(jì)算結(jié)果偏差2.8%～4.5%；當(dāng)測(cè)量的初始溫度偏差為1℃時(shí)，計(jì)算出來(lái)的導(dǎo)熱系數(shù)變化約為12%～14.6%；當(dāng)管間距測(cè)量誤差為0.01 m時(shí)，導(dǎo)熱系數(shù)變化約為4%～8%[3]。目前，常用于確定巖土熱物性的方法有穩(wěn)態(tài)平板測(cè)試法、現(xiàn)場(chǎng)測(cè)試法、土壤類型辨別法和探針?lè)ǖ?，相關(guān)研究表明現(xiàn)場(chǎng)測(cè)試法得到的結(jié)果最接近真實(shí)情況[4]。

目前用于計(jì)算周圍巖土與地埋管換熱器換熱的模型主要有線熱源模型和柱熱源模型，而基于這兩種模型又有不同的數(shù)據(jù)處理方法，分別為最小二乘法和參數(shù)估計(jì)法。AUSTIN[5]和顏亮等[6]的研究發(fā)現(xiàn)雖然基于線熱源模型比柱熱源模型的計(jì)算結(jié)果更加準(zhǔn)確，但是線熱源模型不能提供巖土的容積比熱容。盧軍等[7]進(jìn)行了定熱流熱響應(yīng)實(shí)驗(yàn)，采用線熱源模型和柱熱源模型得到的巖土熱物性參數(shù)有較大差距，導(dǎo)熱系數(shù)相對(duì)差距在8.8%～10.4%。于明志等[8]和楊衛(wèi)波等[9]的實(shí)驗(yàn)研究表明，對(duì)于采用線熱源模型數(shù)據(jù)擬合法來(lái)說(shuō)，在測(cè)試時(shí)間約為30 h后，導(dǎo)熱系數(shù)開(kāi)始穩(wěn)定，表明延長(zhǎng)測(cè)試的時(shí)間可以有效提高計(jì)算結(jié)果的穩(wěn)定性。章云根等[10]搭建巖土熱響應(yīng)實(shí)驗(yàn)平臺(tái)，并基于線熱源模型和非線性最小二乘法建立計(jì)算模擬程序，得到結(jié)果的誤差在8%以內(nèi)。常桂欽等[11]基于線熱源模型采用參數(shù)估計(jì)法對(duì)地埋管理論值和模擬值進(jìn)行了優(yōu)化，使模擬值盡量接近真實(shí)情況。清華大學(xué)的桂樹(shù)強(qiáng)等[12]對(duì)比了埋管換熱器的實(shí)驗(yàn)值和數(shù)值模擬計(jì)算值，結(jié)果顯示吻合度較好，因此數(shù)值模擬作為更簡(jiǎn)便易行的研究工具，可有效提高研究效率。

本文采用現(xiàn)場(chǎng)測(cè)試方法搭建了地埋管換熱測(cè)試系統(tǒng)，對(duì)青島藍(lán)色硅谷區(qū)域地埋管換熱器的換熱性能以及土壤的熱物理性能進(jìn)行測(cè)試，同時(shí)針對(duì)垂直地埋管的地源熱泵系統(tǒng)，對(duì)現(xiàn)場(chǎng)地質(zhì)情況進(jìn)行了解，基于線熱源理論測(cè)得測(cè)試孔的熱響應(yīng)曲線，計(jì)算得到巖土熱物性參數(shù)，并結(jié)合柱熱源理論的模擬計(jì)算結(jié)果進(jìn)行對(duì)比分析，以期為土壤換熱器提供更經(jīng)濟(jì)的設(shè)計(jì)，并且在測(cè)試孔的鉆井過(guò)程中，了解到項(xiàng)目所在地的地質(zhì)情況，以此確定青島藍(lán)色硅谷是否可以適用土壤源熱泵技術(shù)，為將來(lái)施工機(jī)械使用以及工期安排提供參考。

1 巖土熱響應(yīng)試驗(yàn)

巖土熱物性現(xiàn)場(chǎng)測(cè)試法又稱為熱響應(yīng)測(cè)試法，為了更準(zhǔn)確地確定地埋管地源熱泵系統(tǒng)鉆井的深度與數(shù)量，必須通過(guò)巖土熱響應(yīng)試驗(yàn)獲取土壤的導(dǎo)熱系數(shù)等熱物性參數(shù)。根據(jù)測(cè)試儀采集到測(cè)量溫度、水流量、電功率等相關(guān)參數(shù)并根據(jù)建立的地埋管與周圍巖土換熱模型進(jìn)行計(jì)算。

試驗(yàn)的過(guò)程為在已鉆好的測(cè)試孔中埋設(shè)PE管，并參照設(shè)計(jì)技術(shù)規(guī)范[13]進(jìn)行回填作業(yè)，地埋管內(nèi)充滿水，待測(cè)試孔溫度場(chǎng)經(jīng)過(guò)48 h恢復(fù)至初始狀態(tài)后，循環(huán)系統(tǒng)開(kāi)始工作，使管內(nèi)的水處于循環(huán)流動(dòng)狀態(tài)，排氣工作完成后，測(cè)試埋管的狀態(tài)基本穩(wěn)定，測(cè)試流體的進(jìn)出口水溫趨于恒定，這時(shí)可以認(rèn)為以環(huán)路進(jìn)出水平均溫度作為埋管換熱器埋深范圍內(nèi)的土壤初始平均溫度。之后開(kāi)啟電加熱系統(tǒng)，使循環(huán)水在地層散失的熱量得到補(bǔ)充，通過(guò)恒定功率加熱，以維持地下放熱率的恒定，測(cè)量系統(tǒng)運(yùn)行，其主要是通過(guò)兩個(gè)溫度傳感器和一個(gè)流量計(jì)來(lái)測(cè)試進(jìn)出口水溫和循環(huán)水流量，通過(guò)加熱功率、埋管出口及進(jìn)口水溫和流量參數(shù)，可計(jì)算得出該地源熱泵系統(tǒng)的土壤綜合導(dǎo)熱系數(shù)?，F(xiàn)場(chǎng)測(cè)試系統(tǒng)裝置示意圖如圖 1所示。溫度傳感器、流量計(jì)、電加熱器循環(huán)水泵均集中置于熱物性測(cè)試箱體中，設(shè)定加熱功率后電加熱器以恒定功率對(duì)地埋管內(nèi)的流體進(jìn)行加熱，流量計(jì)、溫度傳感器等每隔一段時(shí)間采集一次數(shù)據(jù)，采集到的數(shù)據(jù)均儲(chǔ)存于自動(dòng)記錄儀中，各測(cè)試參數(shù)變量及其精度如表1所示。

圖1 熱物性測(cè)試系統(tǒng)裝置示意圖Fig. 1 Thermal property test system installation

表1 現(xiàn)場(chǎng)試驗(yàn)測(cè)試參數(shù)及精度表Table 1 Field experiment test parameters and accuracy table

試驗(yàn)地點(diǎn)位于青島即墨市某地?zé)豳Y源勘探項(xiàng)目所在地。在進(jìn)行熱物性測(cè)試之前，對(duì)該區(qū)域的地質(zhì)進(jìn)行勘探，該區(qū)域的地質(zhì)構(gòu)造如表2所示。

表2 測(cè)試地點(diǎn)豎直方向地質(zhì)構(gòu)造Table 2 Geological structure in the vertical direction of the test area

該鉆孔的深度為100 m，采用的是雙U型埋管形式。地埋管采用的是高密度聚乙烯管材（HDPE管），埋管內(nèi)徑為26 mm，外徑為32 mm。放管結(jié)束后進(jìn)行原漿回填作業(yè)。

2 地埋管傳熱模型及算法

2.1 線熱源模型

1948年，INGERSOLL等[14]在前人研究的基礎(chǔ)上提出了線熱源模型，該模型將地埋管換熱器與周圍接觸的土壤之間的傳熱看作為恒定線熱源傳熱，測(cè)試井看作為無(wú)限長(zhǎng)圓柱體的一維非穩(wěn)態(tài)導(dǎo)熱問(wèn)題。為簡(jiǎn)化模型，引進(jìn)如下假設(shè)：（1）鉆孔周圍土壤的熱物特性是均勻的；（2）埋管與周圍土壤的換熱認(rèn)為是鉆孔中心的一根線熱源與周圍土壤進(jìn)行換熱，沿長(zhǎng)度方向傳熱量忽略不計(jì)；（3）地埋管與周圍巖土的換熱強(qiáng)度維持不變（試驗(yàn)時(shí)通過(guò)控制加熱功率實(shí)現(xiàn)）。無(wú)限長(zhǎng)線熱源模型既可以用于定熱流條件的熱響應(yīng)實(shí)驗(yàn)，也可用于定埋管出口溫度條件下的熱響應(yīng)實(shí)驗(yàn)。本文試驗(yàn)采用恒熱流條件的熱響應(yīng)實(shí)驗(yàn)，利用數(shù)據(jù)擬合[9]對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析。其溫度解析式為：

其中，Tf=(Tg,in+Tg,out)/2，Tg,in和Tg,out分別為地埋管的進(jìn)、出口溫度的測(cè)量值，℃；π為圓周率；Tsur為t時(shí)刻土壤源處邊界的初始溫度，℃；Q為地埋管加熱功率，W；λ為地埋管換熱器的綜合導(dǎo)熱系數(shù)，W/(m·K)；H為埋管深度，m；t為測(cè)試持續(xù)的時(shí)間，s；a為土壤的熱擴(kuò)散率，m2/s；Rb為土壤的平均導(dǎo)熱熱阻，K/(W·m)；γ為歐拉系數(shù)，取為0.577 2；rb為測(cè)試孔的外徑，m。

式（1）可以簡(jiǎn)化為：

式中，m和k都是定值。其計(jì)算式分別為：

其中，鉆孔熱阻根據(jù)公式（3）可知k是導(dǎo)熱系數(shù)的一個(gè)參數(shù)，從而可以得到導(dǎo)熱系數(shù)為：

通過(guò)實(shí)驗(yàn)測(cè)量得到的加熱功率及不同時(shí)刻地埋管內(nèi)流體平均溫度Tf值，利用斜率法，繪制Tf隨ln(t)變化的關(guān)系曲線，采用最小二乘法擬合得到式（2）的斜率k和截距m，從而可根據(jù)得出的k值計(jì)算出地埋管換熱器的綜合導(dǎo)熱系數(shù)λ值。

2.2 柱熱源模型

柱熱源模型將地埋管用當(dāng)量直徑等價(jià)于有限半徑垂直圓柱管，假定土壤各向同性，鉆孔井壁與周圍土壤之間的換熱在恒熱流條件下進(jìn)行。TRNSYS軟件中的地埋管換熱器模塊的數(shù)學(xué)模型建模過(guò)程對(duì)其進(jìn)行了一些假設(shè)，其布置形式及其溫度場(chǎng)、蓄熱區(qū)域的熱容積、埋管附近土壤的熱物特性等都是以換熱區(qū)域?yàn)檩S的中心對(duì)稱，呈圓柱形分布。地埋管換熱器以圓柱為軸對(duì)稱的形式分布于圓柱形區(qū)域內(nèi)，圓柱中心線作為換熱過(guò)程的內(nèi)邊界，沒(méi)有熱量通過(guò)該區(qū)域。外邊界通過(guò)不影響蓄熱體內(nèi)換熱過(guò)程的基準(zhǔn)來(lái)進(jìn)行確定。根據(jù)CARSLAW等給出的解析法[15]，本模型可給出其精確解，假設(shè)沿深度方向單位鉆孔熱阻恒定不變且為Rb，該模型得出的恒定熱流下埋管內(nèi)流體的平均溫度解析解為：

地埋管周圍土壤的平均溫度為：

式中，r1為地埋管的管間當(dāng)量半徑，即視為絕熱邊界，在該半徑范圍之外幾乎沒(méi)有溫度擾動(dòng)；q為恒定加熱量；C為土壤的比熱容。分析式（1）和式（6）可以發(fā)現(xiàn)，兩個(gè)簡(jiǎn)化模型中都有三個(gè)未知數(shù)土壤的綜合導(dǎo)熱系數(shù)λ、地埋管外側(cè)表面單位管長(zhǎng)的熱量傳到附近土壤時(shí)的熱阻Rb和土壤的熱擴(kuò)散率a。利用傳熱反問(wèn)題求解，結(jié)合最優(yōu)化方法同時(shí)確定Rb和λ，該問(wèn)題就變成為Rb和λ的雙參數(shù)估計(jì)問(wèn)題。

參數(shù)估計(jì)法的基本思想是通過(guò)不斷調(diào)整待求參數(shù)，當(dāng)實(shí)驗(yàn)記錄水溫對(duì)應(yīng)時(shí)刻的理論計(jì)算溫度方差求和最小時(shí)的參數(shù)值為最優(yōu)計(jì)算結(jié)果。即保證目標(biāo)函數(shù)f最小，其計(jì)算式為：

式中，f為目標(biāo)函數(shù)，Tcal,i為i時(shí)刻的計(jì)算水溫，℃；Texp,i為i時(shí)刻的實(shí)驗(yàn)測(cè)試水溫，℃；N為按照時(shí)間順序?qū)嶒?yàn)測(cè)試的數(shù)據(jù)組數(shù)。

本實(shí)驗(yàn)采用的是通過(guò)控制現(xiàn)場(chǎng)測(cè)試的恒定加熱功率來(lái)滿足常熱流邊界條件。因此通過(guò)式（1）和式（6）計(jì)算得到埋管內(nèi)流體的平均溫度與實(shí)際測(cè)試得到的數(shù)據(jù)進(jìn)行對(duì)比，利用參數(shù)估計(jì)法，對(duì)熱物性參數(shù)不斷進(jìn)行調(diào)整，直至得到土壤熱物性參數(shù)的優(yōu)化值。

3 數(shù)據(jù)計(jì)算與結(jié)果分析

3.1 實(shí)測(cè)數(shù)據(jù)處理與分析

基于恒熱流測(cè)試原理，利用巖土熱響應(yīng)測(cè)試儀對(duì)青島即墨市的地埋管測(cè)試孔的巖土熱物性進(jìn)行現(xiàn)場(chǎng)測(cè)試。鉆孔孔徑為150 mm，測(cè)得土壤的初始溫度為16.4℃。設(shè)置恒定加熱功率為7 kW，測(cè)試48 h的地埋管進(jìn)出口平均水溫如圖2所示。測(cè)試裝置初期處于試運(yùn)行階段，需要進(jìn)行排氣和測(cè)試原始土壤溫度等操作，因此從圖中可以看出測(cè)試的前0.5 h的數(shù)據(jù)變化較大，此時(shí)的測(cè)試數(shù)據(jù)與實(shí)際測(cè)試井的傳熱并不相符。為了排除這些不良影響，本文將測(cè)試的前0.5 h數(shù)據(jù)予以舍棄。

圖2 現(xiàn)場(chǎng)測(cè)試地埋管流體進(jìn)出口水溫變化曲線圖Fig. 2 Inlet and outlet water temperature curve of on-site testing of buried pipe fluid

根據(jù)測(cè)試得到的48 h數(shù)據(jù)，舍棄測(cè)試前期的不穩(wěn)定數(shù)據(jù)，利用ORIGIN軟件對(duì)平均溫度曲線圖分析后，得出平均溫度時(shí)間對(duì)數(shù)曲線圖及相應(yīng)公式，其計(jì)算決定系數(shù)（即擬合優(yōu)度）R2為0.995 4，說(shuō)明其擬合的回歸線對(duì)測(cè)量值擬合的程度很好。

由圖3對(duì)數(shù)曲線圖得出下式：

對(duì)應(yīng)公式（2）及2.1所述相關(guān)原理，得出k=2.446 7。

將公式（5）代入相關(guān)數(shù)據(jù)得到土壤綜合導(dǎo)熱系數(shù)：

圖3 測(cè)試井環(huán)路平均溫度回歸曲線Fig. 3 Average temperature regression curve of the test well

測(cè)試時(shí)長(zhǎng)對(duì)土壤的綜合導(dǎo)熱系數(shù)的影響情況如圖4所示。 分析圖4可發(fā)現(xiàn)，隨有效時(shí)間t的后移，土壤綜合導(dǎo)熱系數(shù)下降，在開(kāi)始的12 h，導(dǎo)熱系數(shù)下降非常快，在試驗(yàn)過(guò)程中有兩個(gè)換熱過(guò)程，即電加熱器對(duì)流體進(jìn)行加熱，流體通過(guò)壁面與土壤進(jìn)行換熱。由于通過(guò)恒定功率對(duì)埋管內(nèi)的流體進(jìn)行加熱，埋管內(nèi)的流體邊加熱邊與地埋管周圍的土壤進(jìn)行換熱，前12 h流體平均溫度升高較快，之后升高速度緩慢，從式（1）進(jìn)行分析，其右邊第一項(xiàng)與時(shí)間呈對(duì)數(shù)關(guān)系，第二項(xiàng)為與加熱量和埋管相關(guān)設(shè)置參數(shù)有關(guān)的量。埋管周圍土壤的平均溫度的變化趨勢(shì)為線性變化，可通過(guò)式（7）進(jìn)行表示。由于傳熱過(guò)程是連續(xù)變化的，將每個(gè)點(diǎn)的熱通量視為瞬時(shí)且連續(xù)的變化情況，換熱器進(jìn)口的載熱流體溫度會(huì)隨著管內(nèi)流體流速的變化而變化，前期處于非穩(wěn)定狀態(tài)時(shí)，其流體溫度也會(huì)隨著時(shí)間的遞進(jìn)而發(fā)生變化，可見(jiàn)流體溫度是與時(shí)間有關(guān)的函數(shù)。一段時(shí)間過(guò)后，電加熱器的加熱功率等于管內(nèi)流體通過(guò)壁面與土壤的換熱功率，此時(shí)達(dá)到了穩(wěn)定狀態(tài)，因此流體進(jìn)出口溫度不變，導(dǎo)熱系數(shù)也逐漸趨于穩(wěn)定。說(shuō)明測(cè)試時(shí)間越長(zhǎng)，得到的計(jì)算結(jié)果越準(zhǔn)確，前期受實(shí)驗(yàn)條件和傳熱穩(wěn)定性的影響，計(jì)算結(jié)果會(huì)出現(xiàn)不穩(wěn)定波動(dòng)現(xiàn)象，因此，為了提高計(jì)算結(jié)果的準(zhǔn)確度，可適當(dāng)延長(zhǎng)試驗(yàn)測(cè)試時(shí)間，一般情況下測(cè)試48 h左右即可滿足實(shí)際工程設(shè)計(jì)的要求。

圖4 測(cè)試時(shí)間對(duì)導(dǎo)熱系數(shù)的影響Fig. 4 Effect of test time on thermal conductivity

整個(gè)熱響應(yīng)測(cè)試加熱過(guò)程中每延米換熱量隨時(shí)間變化曲線見(jiàn)圖5。由于初期運(yùn)行時(shí)未開(kāi)啟加熱器加熱，且處于排氣階段，導(dǎo)致其在測(cè)試初期的波動(dòng)較大。在1 h 之后每延米換熱量逐漸趨于平穩(wěn)并接近于70 W/m。為驗(yàn)證試驗(yàn)的準(zhǔn)確性，對(duì)測(cè)試誤差進(jìn)行分析。測(cè)試口換熱趨于穩(wěn)定時(shí)，流體進(jìn)出口平均溫差為4.3℃，流量為1.4 m3/h，根據(jù)得：

單位井深換熱量為：

計(jì)算得到的單位井深換熱量約為 69.89 W/m，因此其準(zhǔn)確率為69.89/70 = 99.8%，說(shuō)明試驗(yàn)結(jié)果與理論值幾乎相同，試驗(yàn)結(jié)果良好，其中的微小誤差是由實(shí)驗(yàn)過(guò)程中加熱功率、測(cè)量的溫度和流量等存在誤差導(dǎo)致的。

圖5 每延米換熱量隨時(shí)間變化曲線Fig. 5 Heat transfer rate per meter over time

測(cè)試孔測(cè)試結(jié)果如表3所示，夏季工況進(jìn)出口水溫約為35.1～30.8℃，溫差為4.3℃。

表3 地源熱泵測(cè)試結(jié)果Table 3 Ground source heat pump test results

3.2 柱熱源模擬分析

TRNSYS軟件對(duì)地埋管換熱器模塊的數(shù)學(xué)模型建模過(guò)程進(jìn)行了一些假設(shè)：其溫度場(chǎng)、蓄熱區(qū)域的熱容積以及埋管附近土壤的熱物特性等均以換熱區(qū)域的中心軸呈圓柱形中心對(duì)稱分布。地埋管換熱器以圓柱為軸對(duì)稱的形式分布于圓柱形區(qū)域內(nèi)，圓柱中心線作為換熱過(guò)程的內(nèi)邊界，沒(méi)有熱量通過(guò)該區(qū)域。外邊界通過(guò)不影響蓄熱體內(nèi)換熱過(guò)程的基準(zhǔn)來(lái)進(jìn)行確定。對(duì)該模塊的相關(guān)參數(shù)進(jìn)行設(shè)置，其主要的設(shè)置參數(shù)如表4所示。

表4 TRNSYS中地埋管換熱器模塊設(shè)置參數(shù)表Table 4 The ground heat exchanger module settings parameter

基于參數(shù)估計(jì)法，利用48 h的測(cè)試數(shù)據(jù)，采用柱熱源模型對(duì)該地埋管換熱器進(jìn)行模擬研究，得到的地埋管進(jìn)出口水溫及平均水溫變化情況如圖 6所示。從圖中可以看出，采用柱熱源模型模擬得到的流體進(jìn)出口溫度與實(shí)驗(yàn)得到的流體進(jìn)出口溫度的變化趨勢(shì)相同，且模擬夏季工況進(jìn)、出水溫度為35.07～30.4℃，與實(shí)驗(yàn)工況進(jìn)、出水溫度 35.1～30.8℃具有很好的吻合度，進(jìn)出口溫差在4.5℃左右。這說(shuō)明柱熱源模型是可靠的，可用于土壤源熱泵的熱物性分析，指導(dǎo)實(shí)際工程。

圖6 柱熱源法模擬地埋管進(jìn)出口水溫變化曲線圖Fig. 6 Inlet and outlet water temperature curve of column heat source simulation of ground source heat pump

4 結(jié) 論

在設(shè)計(jì)土壤源熱泵系統(tǒng)時(shí)，巖土熱物性參數(shù)的確定至關(guān)重要，本文采用現(xiàn)場(chǎng)測(cè)試法得到的實(shí)測(cè)數(shù)據(jù)，并基于線熱源模型計(jì)算了巖土的熱物性參數(shù)，同時(shí)基于參數(shù)估計(jì)法模擬在柱熱源模型下的進(jìn)出口水溫變化，并與實(shí)測(cè)值進(jìn)行比較，計(jì)算結(jié)果表明：

（1）試驗(yàn)測(cè)試地區(qū)土壤平均初始溫度為16.4℃，地溫適宜，利用48 h的試驗(yàn)測(cè)試數(shù)據(jù)，算出土壤導(dǎo)熱系數(shù)為2.278 W/(m·K)，導(dǎo)熱系數(shù)較高，地層熱傳導(dǎo)條件優(yōu)良，可以保證土壤源熱泵的換熱能力，說(shuō)明藍(lán)色硅谷區(qū)域適宜使用地源熱泵技術(shù)。

（2）測(cè)試時(shí)間對(duì)計(jì)算結(jié)果的影響較大，測(cè)試時(shí)間越長(zhǎng)，計(jì)算的導(dǎo)熱系數(shù)越準(zhǔn)確，本試驗(yàn)在測(cè)試36 h后開(kāi)始趨于穩(wěn)定，為了提高結(jié)果的準(zhǔn)確度，通常測(cè)試時(shí)長(zhǎng)為48 h左右可以滿足實(shí)際工程需求。

（3）采用柱熱源模型模擬得到的流體進(jìn)出口溫度與實(shí)驗(yàn)得到的流體進(jìn)出口溫度的變化趨勢(shì)相同，且模擬夏季工況進(jìn)、出水溫度為 35.07℃及30.4℃，與實(shí)驗(yàn)工況進(jìn)、出水溫度 35.1℃及 30.8℃具有很好的吻合度。這說(shuō)明采用柱熱源模型進(jìn)行模擬較為可靠，該研究可用于指導(dǎo)土壤源熱泵的熱物性分析以及藍(lán)色硅谷區(qū)域地?zé)豳Y源的綜合開(kāi)發(fā)。