史啟朋，吳曉華，姜福紅，李霖，宋帥良，吳延明，孟甲

(山東省魯南地質(zhì)工程勘察院，山東 兗州 272100)

0 引言

近年來隨著地源熱泵技術的成熟及國家對環(huán)境保護工作的重視，淺層地溫能這種清潔環(huán)保的可再生能源在現(xiàn)代城市建筑中得到了廣泛的應用[1-2]，2007年底全國利用淺層地溫能供暖制冷的建筑面積8000萬m2，2012年4月底超過2.4億m2，2014年底達到3.6億m2。淺層地溫能開發(fā)利用在快速發(fā)展的同時，也存在著工程設計參數(shù)不足和監(jiān)督管理機制不健全等制約因素[3-4]。

為研究巖土體熱物性參數(shù)特征，李少華等[5]對嘉興市一試驗孔分別進行連續(xù)15d和6.25d的巖土熱響應試驗，研究了測試時間、舍棄初始時間、加熱功率對埋管深度范圍內(nèi)巖土體導熱系數(shù)測試的影響；畢文明等[6]在室內(nèi)搭建了一套地埋管換熱器綜合微縮實驗臺，開展了干砂、飽和砂及飽和砂在不同滲流速度下?lián)Q熱參數(shù)研究；韓斯東等[7]引入分布式光纖測溫系統(tǒng)，進行了試驗場區(qū)巖土體分層熱物性測試，取得了一定的研究成果。這些研究僅局限于特定的巖土層，目前缺少不同水文地質(zhì)單元，不同巖土體熱物性能參數(shù)的對比試驗研究。

因此，進行不同水文地質(zhì)單元，不同巖土體熱響對比試驗，求取熱物性參數(shù)[8-9]，對指導類似地區(qū)淺層地溫能開發(fā)利用工程應用，合理利用淺層地溫能資源具有重要意義。

1 地熱地質(zhì)概況

1.1 區(qū)域地熱地質(zhì)條件

魯西地區(qū)主要分布黃河沖洪積平原、山前沖洪積平原和中低山丘陵等水文地質(zhì)單元[10-11]，各水文地質(zhì)單元地熱地質(zhì)條件如下：

(1)黃河沖洪積平原。 位于京杭運河以西地區(qū)，地貌類型為沖洪積平原及沖湖積平原，第四系較發(fā)育，由東向西厚度逐漸增加，最大厚度大于300m，下伏新近紀碎屑巖。恒溫帶埋深20～25m，溫度15.5～16.0℃。大地熱流值為56.0～79.0mW/m2，地溫梯度為2.5～4.0℃/hm，埋深100m處平均地溫17.0～19.0℃。含水巖組顆粒較細，巖性為粉土、粉質(zhì)粘土、細砂—粉細砂。含水層厚度15.25～52.09m，滲透系數(shù)2～3m/d。

(2)山前沖洪積平原。 位于京杭運河以東的山前平原地區(qū)，地貌類型為沖洪積平原，第四系平均厚度在150m左右，下伏新近紀、石炭-二疊紀和奧陶紀地層。恒溫帶埋深為23～30m，溫度16.3℃左右。大地熱流值為52.0～68.0mW/m2，地溫梯度為1.20～3.0℃/100m，埋深100m處平均地溫16～18℃。含水巖組顆粒較粗，巖性主要為粉質(zhì)粘土、粘土、中粗砂。含水層平均厚度32m左右，井孔單位涌水量大于500m3/d·m。

(3)中低山丘陵區(qū)。 位于泰安—萊蕪—棗莊—臨沂一帶基巖山區(qū)，地貌類型以丘陵山地為主，其次為山間平原與山前沖積平原。地層巖性以寒武—奧陶紀碳酸鹽巖及太古宇花崗巖、變質(zhì)巖為主，區(qū)內(nèi)第四系不發(fā)育，厚度普遍小于30m。恒溫帶埋深受地層巖性的影響，碳酸鹽巖發(fā)育區(qū)約25m，太古宇花

崗巖、變質(zhì)巖發(fā)育區(qū)約22m，溫度12.0～13.5℃。大地熱流值為35.0～60.0mw/m2，地溫梯度為2.5～4.3℃/hm，埋深100m處平均地溫15～18℃。地下水類型以碳酸鹽巖巖溶裂隙水及花崗巖、變質(zhì)巖裂隙水為主，在山間平原及山前沖積平原區(qū)分布有松散巖類孔隙水。

1.2 試驗場地水文地質(zhì)條件

RS2和RS3試驗孔位于黃河沖洪積平原區(qū)，地層巖性主要為粉土、粉質(zhì)粘土、粉細砂等；RS4孔位于山前沖洪積平原，地層巖性以粘土、粉質(zhì)粘土、中砂為主；RS6、RS7和K13孔位于中低山丘陵區(qū)，RS6孔地層巖性為花崗巖；RS7和K13地層巖性為粉質(zhì)粘土、灰?guī)r，RS7孔位于羊莊盆地巖溶水徑流排泄區(qū)，裂隙巖溶發(fā)育，地下水流徑流速度較快*山東省魯南地質(zhì)工程勘察院，山東省淺層地熱能調(diào)查評價(魯西)報告，2012年。*山東省魯南地質(zhì)工程勘察院，臨沂、棗莊、濟寧、菏澤、泰安市淺層地溫能調(diào)查評價報告，2016年。[12]，各試驗孔水文地質(zhì)條件見表1。

表1 試驗孔地層水文地質(zhì)條件

2 試驗設備及場地工程布置

2.1 試驗設備

該次熱響應試驗采用天津地熱勘查開發(fā)設計院生產(chǎn)的FTPT1-1型地層熱響應測試儀進行。該儀器由主機循環(huán)泵、電加熱器、水箱、溫度和流量測試組件、控制與記錄組件、電源控制、管道等部件組成，可完成無負荷循環(huán)試驗、加熱恒功率等試驗。

2.2 場地工程布置

RS2試驗孔位于菏澤市鄆城縣鄆城鎮(zhèn)，孔徑180mm，下管類型PE32雙U，下置深度120.9m，回填料為中粗砂，場地淺層地下水位埋深3.762m。RS3試驗孔位于菏澤市成武縣成武鎮(zhèn)，孔徑180mm，下管類型PE32雙U，下置水深度120.5m，場地淺層地下水位埋深2.609m。RS2和RS3試驗孔目的層位為黃河沖洪積平原粉土、粉質(zhì)粘土及粉細砂等。

RS4試驗孔位于濟寧兗州興隆莊鎮(zhèn)，孔徑180mm，下管類型PE32雙U，下置深度100.5m，回填料為中粗砂，場地淺層地下水位埋深5.622m。試驗孔目的層位為汶泗河沖洪積平原粉質(zhì)粘土、粘土及中粗砂等。

RS6試驗孔位于濟寧泗水泗張鎮(zhèn)，孔徑180mm，下管類型PE32雙U，下置深度70.0m，回填料為中粗砂，場地淺層地下水位埋深13.410m。試驗孔目的層位為花崗巖。

RS7試驗孔位于棗莊滕州羊莊鎮(zhèn)，孔徑180mm，下管類型PE32雙U，下置深度100.0m，回填料為中粗砂，場地淺層地下水位埋深9.504m。試驗孔目的層位為灰?guī)r、泥灰?guī)r。K13試驗孔位于濟南市高新區(qū)港溝街道辦事處,孔徑180mm，下管類型PE32雙U，下置深度120.5m，回填料為中粗砂，試驗孔目的層位為灰?guī)r、泥灰?guī)r。

3 工藝流程及參數(shù)的影響因素分析

3.1 試驗工藝流程

(1)該次熱響應試驗是在試驗孔鉆探成孔，PE32雙U地埋管打壓測試，下管、回填礫料，再次試壓和管口密封靜止48h后進行。以RS4孔的巖性地層為例，試驗工藝流程詳見圖1。

(2)試驗時將熱響應測試儀主機與地埋管相連接，組成一個閉合回路，向回路內(nèi)注入干凈的自來水，接入三項五線380V電源，開啟主機水泵讓系統(tǒng)循環(huán)，檢查系統(tǒng)是否漏水。

圖1 熱響應試驗工藝流程圖

(3)進行加熱4kW測試，加熱負荷和流量基本保持恒定(波動范圍在±5%)，換熱管內(nèi)傳熱介質(zhì)流速不應低于0.2m/s，實時記錄回路中傳熱介質(zhì)的流量和溫度，數(shù)據(jù)采集頻率為1次/min。地埋管換熱器出水溫度1h內(nèi)變化小于1℃后，觀測時間不少于24h后試驗結束。

(4)測試結束后，先關閉電加熱器和水泵，再關閉系統(tǒng)，斷開電源，將系統(tǒng)與地埋管連接斷開，排干凈系統(tǒng)內(nèi)的水，封閉地埋管的管口。

3.2 熱物性參數(shù)的影響因素分析

不同巖土體熱物性參數(shù)大小及影響因素與測試設備儀表的精度、測試方法、數(shù)據(jù)分析模型、回填材料、巖土體類型和地下水徑流條件等相關。

測試設備中的溫度傳感器、流量計、功率計等精度和可靠性是影響試驗數(shù)據(jù)準確性的重要因素；穩(wěn)定熱流或穩(wěn)定工況測試方法，獲取的試驗數(shù)據(jù)差別較大；線熱源理論模型或圓柱熱源理論模型分析數(shù)據(jù)，得到的結果各有不同；回填材料對地埋管換熱器換熱性能也有著顯著的影響，回填材料是多孔介質(zhì)，孔隙之間空氣的導熱率最小，因此，回填材料顆粒愈大，孔隙度愈小，導熱率愈大。反之，回填料顆粒愈小，孔隙愈大，導熱率愈小。

另外，巖石地層的導熱性越強，所處地下水徑流循環(huán)速度越快，試驗孔巖土體熱物性參數(shù)值就越大。

該次開展的6組熱響應試驗采用的測試儀器，測試方法及數(shù)據(jù)分析模型相同，回填料均為中粗砂，顆粒大，導熱率大，且各孔回填料粒徑基本相同，可有效避免或減少因測試設備、測試方法、數(shù)據(jù)分析模型、回填材料的不同對試驗結果帶來的影響，得到巖土體真實的熱物性參數(shù)。

4 試驗方法及結果分析

4.1 試驗方法

國內(nèi)外針對巖土體熱物性測試和地埋管換熱器換熱特性測試做了大量理論研究及試驗工作，目前常用的測試方法是由Mogensen提出的熱響應方法(TRT)。這種方法是對測試孔中的U型管施加一個恒定的熱流，通過獲取熱流的加熱功率，U型管進出口流體溫度及流量等參數(shù)，計算巖土體熱物性參數(shù)。地埋管換熱器的熱響應特性試驗在理論上可以歸結為在一定熱流邊界條件下的非穩(wěn)態(tài)傳熱問題。其數(shù)學解析主要有2種模型：①基于線熱源理論的線模型；②基于圓柱熱源理論的柱模型。該次試驗所用設備設計和計算方法均采用Hart和Couvillion建立的線熱源理論數(shù)學模型[13-15]，根據(jù)該理論模型，地埋管換熱器進出水口平均溫度值計算公式為：

(1)

式中：Tf—隨時間變化的地埋管換熱器進出水口平均溫度(k)；Q—單孔換熱量(W)；H—垂直埋管深度(m)；λ—巖土體平均熱導率[W/(m·k)]；a—熱擴散系數(shù)(m2/s)；t—測試時間(s)；r—鉆孔半徑(m)；γ—歐拉常數(shù)(取0.5772)；Rb—鉆孔熱阻(m·k/W)；T0—地層初始溫度(k)。

公式(1)可寫為線性形式，即：

Tf=klnt+m

(2)

(3)

(4)

單孔換熱量計算公式：

Q=G×△t/0.86

(5)

式中：Q—地埋管每小時產(chǎn)生的熱量kW；G—地埋管循環(huán)水流量(m3/h)； △t—地埋管進出水口溫差(℃)。

根據(jù)各試驗孔的測試數(shù)據(jù)，繪制Tf隨lnt變化曲線，求取其斜率(k)，由(2)～(5)式可計算巖土體的平均熱導率λ和每延米換熱量。

4.2 試驗結果

熱響應試驗加熱功率均為4kW，試驗孔地埋管換熱器進出水口溫度隨時間變化情況見圖2。

圖2 試驗孔進出水口溫度隨時間變化曲線圖

圖3 試驗孔擬合曲線圖

根據(jù)加熱4kW測試數(shù)據(jù)繪制Tf(進出水口平均溫度)隨lnt的變化曲線，并選擇測試換熱穩(wěn)定狀態(tài)下的數(shù)據(jù)進行線性擬合，擬合圖及擬合方程式見圖3。每延米換熱量(W/m)和平均熱導率(W/m·k)計算結果見表2。

通過對試驗孔地質(zhì)條件和試驗測試結果綜合分析，初步試驗結果表明：每延米換熱量(W/m)和平均熱導率(W/m·k)花崗巖>灰?guī)r>粉質(zhì)粘土、中粗砂>粉土、粉細砂。在巖溶地下水富水性強，徑流速度快的地段灰?guī)r平均熱導率>花崗巖。熱導率現(xiàn)場試驗結果與同類巖石樣品試驗室測定結果大小規(guī)律基本一致。

4.3 試驗結果分析

根據(jù)該次試驗測試結果和試驗孔處水文地質(zhì)條件，試驗孔每延米換熱量和平均熱導率存在如下特點：

表2 每延米換熱量、平均熱導率計算結果

(1)第四系山前沖積平原粉質(zhì)粘土、中粗砂地層每延米換熱量和平均熱導率大于黃河沖洪積平原粉土、粉細砂地層。

(2)低山丘陵區(qū)基巖每延米換熱量和平均熱導率均大于第四紀松散層，花崗巖大于灰?guī)r。

(3)相同水文地質(zhì)單元，地下水富水性強、徑流速度快的地段，鉆孔每延米換熱量和平均熱導率大。

5 結論及建議

(1)魯西黃河沖積平原粉土、粉質(zhì)粘土、細砂地層，山前沖洪積平原粉質(zhì)粘土、粘土、中細砂地層，中低山丘陵灰?guī)r和花崗巖地層，平均每延米換熱量分別為36.47W/m，37.70W/m，41.57和62.14W/m;平均熱導率分別為1.49W/m·k,1.93W/m·k,3.18W/m·k和2.65W/m·k。

(2)每延米換熱量和平均熱導率第四系山前沖洪積平原粗顆粒的粉質(zhì)粘土、中粗砂大于黃河沖洪積平原細顆粒的粉土、粉細砂；花崗巖大于灰?guī)r，灰?guī)r大于第四系松散層。熱導率現(xiàn)場試驗結果與同類巖石樣品試驗室測定結果大小規(guī)律基本一致，現(xiàn)場熱響應試驗測試結果可靠性強。

(3)相同水文地質(zhì)單元，地下水富水性強、徑流速度快的試驗孔每延米換熱量和平均熱導率大。

魯西地區(qū)不同巖性地層熱響應試驗結果，可為該區(qū)淺層地溫能開發(fā)利用工程設計、淺層地溫能調(diào)查評價和開發(fā)利用提供基礎數(shù)據(jù)。

[1] 周華慧，王景剛，佘軍.地源熱泵巖土熱響應現(xiàn)場測試及分析[J].制冷與空調(diào)，2012，12(2)：57-60.

[2] 胡波，黎逢良，譚秀全.泗河沖洪積扇松散巖土體地溫場自恢復能力研究[J].山東國土資源，2017，33(7)：55-60.

[3] 王秉沈，田廷山，趙繼昌，等.我國地溫資源開發(fā)與地源熱泵技術應用、發(fā)展及存在問題[J].地熱能，2009(1)：23-26.

[4] 劉春雷，王貴玲，王婉麗，等.基于現(xiàn)場熱響應測試方法的地下巖土熱物性分析[J].吉林大學學報(地球科學版)，2014，44(5)：1602-1608.

[5] 李少華，秦祥熙，牛定輝，等.現(xiàn)場熱響應試驗測定導熱系數(shù)的影響因素分析[J].暖通空調(diào)，2015，45(12)：49-52.

[6] 畢文明，岳麗燕，韓再生，等.地埋管換熱性能綜合微縮試驗研究[J].水文地質(zhì)工程地質(zhì)，2014，41(1)：144-148.

[7] 韓斯，金光，畢文明，等.某地源熱泵場地淺層土壤分層熱物性響應試驗[J].水文地質(zhì)工程地質(zhì)，2015，42(4)：144-149.

[8] 中華人民共和國建設部.地源熱泵系統(tǒng)工程技術規(guī)范(2009年版)(GB/T50366-2005)[S].北京：中國建筑工業(yè)出版社，2009.

[9] 史啟朋.汶泗河沖積扇松散層現(xiàn)場熱響應試驗數(shù)據(jù)分析及應用[J].山東國土資源，2014，30(2)：54-57.

[10] 朱友強，劉桂儀，佟光玉，等.山東省淺層地溫能開發(fā)與應用實例[A]//地溫能資源與地源熱泵技術應用論文集[C].北京：地質(zhì)出版社，2011：3-8.

[11] 周亞醒.魯西地區(qū)淺層地溫能地熱地質(zhì)條件與開發(fā)利用適性評價[J].山東國土資源，2015，31(8)：45-49.

[12] 李傳謨,康鳳新.巖溶水資源及增源增采模型[M].濟南：山東科學技術出版社，1999：22-33.

[13] 刁乃仁.地熱換熱器的傳熱問題研究及其工程應用[D].北京：清華大學，2005：16-69.

[14] 畢文明，郭艷春，王琳.地埋管地源熱泵系統(tǒng)巖土體熱物性參數(shù)研究[A]//地溫能資源與地源熱泵技術應用論文集[C].北京：地質(zhì)出版社，2009：105-111.

[15] 于明志，彭曉峰，方肇洪，等.基于線熱源模型的地下巖土熱物性測試方法[J].太陽能學報，2006，27(3)：279-283.