回填空氣間隙對地埋管巖土溫度恢復性能的影響

王 勇，金逸韜

（重慶大學 三峽庫區生態環境教育部重點實驗室，重慶400045）

地源熱泵是隨全球能源環境問題的可持續發展而興起的一種節能環保的地熱利用技術［1］。在地源熱泵工程中，地埋換熱器周圍巖土溫度的恢復一直是研究、設計人員非常關心的問題，相關學者對此進行了大量的研究［2］。劉俊等［3］提出地溫的恢復特性主要取決于巖土熱物性、管群布置、系統啟停比、冷熱負荷強度和冷熱負荷不平衡率等因素。Shang等［4］研究了間歇運行情況下，影響地溫恢復的因素：巖土及回填材料的熱物性、多孔性對巖土溫度的恢復有較大影響，當多孔性指數下降時，巖土溫度恢復所需的時間減少。趙利君等［5］通過對間歇運行巖土溫度在不同深度及半徑下的恢復率的研究，發現由于巖土對傳熱的衰減作用，離孔壁距離越遠，溫度受到間歇運行工況的影響就越小；巖土恢復能力隨著運行周期的增加而不斷降低。高源等［6］的實驗結果表明，同一埋管深度處，機組的運行時間越長，巖土短期內可恢復的可用穩定溫度越高，且需要的穩定恢復時間也越長。陳穎等［7］通過實驗發現熱泵機組開機12h后，巖土溫度趨于穩定；停機12h后地下換熱器處的巖土溫度可恢復約95%，但基本得以恢復需要停機81h。

由于地埋管的鉆孔造價較高，鉆孔的深度、直徑、管材及鉆井周圍巖土的熱物性都難以隨意更改，近年來巖土源熱泵的換熱性能研究除了較多的集中在回填材料的熱物性上以外，越來越重視回填密實度對鉆孔換熱器的換熱性能的影響。由于人工回填施工困難，會導致回填不密實，增加了傳熱總熱阻，極大的影響了地埋管換熱性能［8－9］。Zhang等［10］建立了一個由泥土、水、空氣混合而成的巖土顆粒模型進行研究。研究表明：巖土的多孔性對巖土的傳熱、蓄熱性能有較大的影響。隨著多孔性數值的增加，巖土的導熱系數和比熱值會下降。朱清宇等［11］發現人工回填事實上更廣泛地應用于各個土壤源熱泵項目之中，但回填密實性的檢測無法進行，且回填密實性在多大程度上影響U形管的換熱效果也未得到合理的解釋。在實際工程中，由于施工人員與施工單位的素質問題，非常不重視回填這道工序，人為的質量問題更加導致了回填的不密實。

筆者在已有研究成果的基礎上，針對工程上易出現的埋管空氣間隙，通過理論計算得到該條件下導致的地溫分布情況，研究其對地埋管巖土溫度恢復性能的影響。

1 地下換熱器與巖土的換熱耦合計算模型

地埋管與巖土之間的換熱是一個不穩定的導熱過程，其傳熱過程復雜且影響因素很多，如管內水流動、材料物性參數、巖土熱物性參數、回填密實度以及地下水分遷移等問題。為了便于計算分析，作出以下簡化：

生產率因素，基于異質性企業貿易原理可得，采用合理措施確保企業生產率呈現大幅度提升，與國外消費者所要求的高質量、高品質、以及高安全度與環保度的要求相符，進入壁壘的能力將得到提升，背負出口初始沉淀成本的能力就越強。

巖土平均過余溫度的計算對象和計算方法已經在前面闡述，根據式（14），可得到圖3中的計算結果。

2）由于巖土的溫度變化范圍不大，假定巖土的導熱系數、比熱、密度等物性參數不隨溫度的變化而變化，且是均勻一致的；

3）地埋管同截面具有相同的溫度和流速；

4）無地下水流動換熱，忽略巖土的濕遷移。

1.1 數學模型的建立

1.1.1 管內流體 對于管內流動，流動為不可壓縮流體，采用標準k－>ε模型，近壁區利用壁面函數法求解。不考慮源項時，連續性方程、動量方程以及運輸方程見式（1）～（5）［12－14］：

1.2.5 巖土及回填上表面 巖土及回填上表面設為壁面，選擇對流換熱邊界條件，考慮到巖土及回填上表面直接與空氣接觸，傳熱量與室外空氣溫度、風速、輻射等因素有關，定義為第三類邊界條件，邊界面周圍空氣溫度設為31℃（夏季室外計算平均溫度），經計算，邊界面與空氣之間的表面傳熱系數h＝1.2W／（m2·K）［16］。

動量方程：

王祥一哆嗦，感覺天都要塌了，自己夸下海口的這批玉器要是賣不出價，就真的沒臉回去見家鄉父老了。他又轉念一想，老道士久經江湖，說不定對錢數的認識和自己不一樣。自己只要10萬就滿足了，說不定……想到這里，王祥馬上追問老道：“道長，你覺得我這批玉器能值多少錢呢？”

能量方程：

Given the rarity of this disease there is no data on the incidence, relative risks, or life time risks of SBA and at present no guidelines exist for surveillance. Some authors do recommend upper endoscopy every 3-5 years from age 15, and repeated annually if polyps are diagnosed[42].

湍動能方程：

耗散率方程：

式中湍動粘度μt和由于平均速度梯度引起的湍動能k的產生項Gk的表達式見式（6）和式（7）：

式（1）～（7）中各常數的取值為：

1.1.2 固體部分 無內熱源非穩態導熱方程見式（11）和式（12）：

式中：kp為管壁傳熱系數；tf為流體溫度；tp為管壁溫度。其中在制冷工況下巖土對流體為冷卻作用，對流換熱系數h采用式（9）和迪圖斯－貝爾特公式（10）求解：

(1) 外圍地區車站差異分析：從天津的城市開發情況來看，居住人口仍然集中在城市中心區，外圍車站所處的快速環路以外區域的土地利用強度不高、人口較為分散，因此按照600 m半徑計算的各類指標均較少。而且城市外圍地區車站的接駁條件較好，實際服務范圍應大于600 m，車站越靠近外環線，服務范圍就越大，因此測算的評價指標低于客流表現。

對管內流體和管壁換熱見式（8）：

式中，k為傳熱系數，W／（m·K）；ρ為密度。

不考慮接觸熱阻，則不同材料接觸面處溫度相等：t1＝t2，見式（13）：

式中：λf為流體導熱系數；r為管半徑。

1.2 邊界條件和初始條件設置

圖6顯示了計算機仿真和實驗室測試結果。如圖所示,濾波器測量結果與計算機仿真結果吻合良好。所設計的濾波器具有4個帶通頻道,中心頻率分別為1.27 GHz、2.2 GHz、2.63 GHz、3.4 GHz,相應的3 dB相對帶寬分別為9.84%、13.25%、2.66%、2.35%,最小插入損耗為分別為1.5 dB、1.1 dB、1.7 dB、1.9 dB,濾波器通帶范圍內回波損耗低于-13 dB。另外,本文所設計的四頻帶通濾波器具有多個傳輸零點,進一步提高了濾波器的頻率選擇性。

1.2.2 埋管出口 埋管的出口各參數都由埋管進口參數及管內流動換熱情況確定，設為壓力出口。

1.2.3 U 型管管壁

U型管管壁設為壁面［15］，對于流動方程，是固定、無滑移壁面，對于邊界節點速度為零；對于能量方程，選擇耦合的傳熱條件，是位于管內流體和管外回填材料這兩個區域間的壁面。

1.2.4 回填側壁 回填側壁設為壁面，是耦合的傳熱壁面，位于回填材料和巖土這兩個區域間。

連續性方程：

1.2.6 遠邊界巖土及底部巖土 遠邊界巖土及底部巖土設為壁面，選擇給定的壁面溫度20℃。

1.3 回填條件的設置

根據回填物與孔壁和換熱器管壁之間不同的回填密實度，在系統連續運行9h及恢復15h的運行工況下，建立了3種不同回填工況：1）回填物與孔壁及換熱器管壁間均回填密實，即無空氣層工況，簡稱工況1；2）回填物與孔壁存在空氣層工況，簡稱工況2；3）回填物與換熱器管壁存在空氣層工況，簡稱工況3。具體設置見表1。工況3中，為了便于模型的建立和網格劃分，將流體和回填材料之間3mm壁厚的PE管和1mm厚的空氣薄層組合為一種材料，簡稱“綜合材料”，根據兩種材料組合前后傳熱熱阻、密度、比熱計算得到綜合材料的熱物性參數。

表1 3種不同回填工況空氣層的設置

1.4 計算條件

以巖土源熱泵單U型換熱器為模型，埋管深80m，鉆孔直徑130mm，流量0.82m3／h，輸入巖土的冷負荷為3.347kW（恒定）。各材料的熱物性參數如表2所示：

1.2.1 埋管進口 埋管的進口設為速度進口［15］，速度值設定為0.43m／s（根據流量及內管徑計算得出），溫度由前一次計算的埋管出口溫度及冷凝器散熱量決定，通過UDF函數輸入，初始溫度設為21℃。

安徽省政府近日審議通過《安徽省淮河流域綜合治理規劃綱要》，根據規劃，到2020年前，安徽省將投資1929億元啟動新一輪淮河治理，標志著安徽治淮工作從防汛抗旱步入減災與防污、供水并重的綜合治理階段。

表2 各材料熱物性參數

整個模型為一直徑3m，高100m的圓柱體，采用ANSYS建模，網格劃分如圖1所示。

圖1 模型網格劃分圖

1.5 巖土平均過余溫度及其計算方法

研究巖土溫度的恢復情況，目的是考查下一個運行周期前的巖土溫度，使得下一個運行周期盡量不受到上一個運行周期負荷痕跡的影響。筆者主要研究回填空氣間隙對埋管周圍巖土溫度恢復性能的影響，基于以下2點，選擇研究溫度恢復的主要對象：

1）近管壁處的巖土溫度是影響地埋管下一個運行周期埋管換熱性能的主要因素。

2）假設3種不同的回填工況，其中工況2和工況3的空氣層均位于回填區域內。

(2)高磁、激電兩種地球物理方法是尋找隱伏、半隱伏接觸交代矽卡巖礦體的常用方法，但在物探異常解釋過程中，需對獲取的疊加異常進行處理。通過向上、向下延拓，正確區分深、淺部異常，綜合各類異常特征，進行綜合評價和解釋，推測成礦富集的有利地段。

故選擇回填孔壁及其以內的回填區域作為計算巖土平均過余溫度的對象。

計算方法如式（14）、（15）所示：

1）對于某一深度z處的回填區域截面，平均過余溫度為：

2）對于整個回填區域體，平均過余溫度為：

根據情景C、D，把種植業結構調整與高效灌溉技術結合起來設定情景E：2020年糧、經、草種植比例為41. 3∶45. 3∶13. 4，微灌面積達70%，2025年糧、經、草比例達41∶40∶19，微灌面積達100%。

總的來說，2000年后揚州市水環境總體生態安全預警狀態逐年好轉，警情有所緩解，從評價結果來看，揚州市水環境生態安全的狀態評價值有下降趨勢，經濟社會發展以及自然水環境條件的下降是揚州市水環境生態安全的主要警源。雖然揚州市采取了積極的措施，但是目前還不能完全恢復水環境的生態功能，受外界干擾后容易惡化。

式中，σz為深度為z米處回填區域截面；txy為σz區域內各單元面積的溫度值；V為整個回填區域體積；txyz為V區域體內各單元體積的溫度值，T0為巖土初始溫度值。

2 模型計算及分析

同種回填工況下，同半徑處不同深度的巖土恢復初始溫度的最大溫差出現在出水管壁處，為1.2～1.4℃。考慮到豎向的巖土初始過余溫差不是很大，在最初共24h的運行與恢復時間段內，同半徑處不同深度的巖土溫度恢復情況基本相同，選取了深度為25m的典型平面作為半徑方向溫度恢復的研究對象［4，17］。深度25m 處，3種回填工況在恢復期的巖土過余溫度如表3、圖2所示。

把試卷中三種題型實際難度做為自變量(第一大題為x1，第二大題為x2，第三大題為x3)，整卷的實際難度做為因變量y，通過標準化轉化，然后建立多元回歸分析模型，基于表1數據可以得到三個自變量的回歸系數(P值)為0.331(P<0.001)、0.239(P<0.001)和0.62(P<0.001).

圖2 深度25m處不同回填工況巖土溫度分布

從表3中可以發現，停機時刻，3種回填工況在深度25m處的熱擴散影響半徑均為r＝0.1～0.5m，近管壁處（管壁至r＝0.1m）的巖土過余溫度均在停機時刻達到了最大值。與工況1相比，工況2在回填孔壁及其以內巖土區域的溫度上升了近3.5℃，而r＝0.1m處的巖土溫度略有下降。這表明位于孔壁處的1mm厚空氣層把由內而外擴散的熱量阻隔在鉆孔內，使得回填孔壁以內的巖土溫度上升，不能有效的將輸入熱量向孔壁外的大地進行傳熱。類似的，與工況1相比，由于空氣層存在于進、出水管壁處，工況3在此處的溫度上升了近8.5℃，換熱器管壁與回填物之間存在的空氣層導致換熱器無法有效的與大地進行換熱。

表3 深度25m處恢復工況下不同半徑處巖土過余溫度 ℃

巖土溫度恢復的最初3h內，溫度下降較快，尤其是工況2的回填孔壁內區域以及工況3的管壁區域。這是因為在恢復初期，溫度較高的巖土區域存在較大的溫度梯度，使得傳熱速度加快。隨著恢復的進行，3種回填工況巖土溫度恢復的速度大幅減小，并在恢復9h后基本趨于一致。

相反的，停機后3種回填工況在r＝0.5m處的巖土溫度卻略微上升，近管壁處的熱量以不穩定狀態從靠近埋管的區域轉移到了靠近遠邊界的巖土，形成了蓄熱。而r＝1.0m、1.5m處的巖土溫度仍然保持初始溫度，說明在運行和恢復總共24h內，r＝0.5m以外的巖土幾乎沒有受到熱擴散的影響。

1）埋管內液體的流速在徑向上均勻一致（忽略重力對流速的影響）；

圖3 深度25m處不同回填工況巖土平均過余溫度

從工況2與工況1的對比分析可知，由于工況2存在孔壁空氣層，導致熱量積累集中在孔洞中，使得孔洞內的初始平均過余溫度高于工況1（兩者分別為8.97、5.39℃）。恢復開始后，工況2的平均過余溫度始終高于工況1，在恢復前6h內，溫度下降速度較快，9h后兩者的平均過余溫度基本趨于一致（0.80℃左右）。

工況3與工況1相比，由于熱源熱量集中在進、出水管壁附近，相同情況下通過管壁傳導到孔洞回填材料中的熱量比工況1少，因此，孔洞內的初始平均過余溫度低于工況1（兩者分別為5.15、5.39℃）。但在恢復開始后，由于積聚在管壁附近的熱量得到迅速擴散，使得工況3的平均過余溫度下降的速度比工況1慢。從圖3中可以發現，從恢復后0～3h間的某個時刻起，工況3的平均過余溫度略高于工況1，6h后兩者的平均過余溫度基本趨于一致。

草兒的到來，使牧兒終于結束了二十多年的單身生活，使得幾間空落落的破舊的屋子里有了色彩，如一朵毫無生機的花兒，重放光彩。草兒買回了油光紙，把墻壁，頂棚糊了，房子就亮堂起來。草兒把牧兒的衣服洗得干干凈凈，疊放得整整齊齊。被子里外都洗了，疊得有棱有角。鞋子洗干凈了，整整齊齊地放在床下。家具、鋤頭也整整齊齊地放好了。往日散亂的一切整齊了，往日臟的一切干凈了。屋子里似變了一個環境，讓牧兒感到無比歡喜，他臉上洋溢著幸福，充滿了朝氣。

工況2與工況3相比，前者的空氣層包圍了整個平均過余溫度的計算區域；而后者的空氣層只是計算區域中靠近管壁處的一部分，由于空氣層的熱阻效應，介于空氣層與孔壁之間的其他計算區域吸收的熱量較工況2少。故工況2的平均過余溫度始終高于工況3。

由式（15）計算得出了豎向80m回填孔洞體積內的巖土平均過余溫度，如表4所示。通過與深度25m處回填孔洞平面的巖土平均過余溫度比較發現，兩者的數值相差在±2%以內，說明在土壤溫度恢復的15h內，深度25m處的回填土壤溫度變化情況與整個回填體積內的土壤溫度變化情況相接近。

表4 深度25m處回填孔洞平面（＊）與豎向80m回填孔洞體積（＊＊）的巖土平均過余溫度 ℃

3 計算模型與實際運行工況的對比

3.1 實際工程情況

根據重慶某公司的巖土源熱泵地下換熱器熱響應實驗，該埋管深50m，鉆孔直徑130mm，冷負荷為恒定的2.092kW，流量0.82m3／h，地埋管初始進水溫度為21℃。經測試，巖土初始地溫為20℃左右，地面10m以下地溫豎向分布基本一致。巖土及材料的熱物性參數與筆者計算情況相接近。根據施工的情況，其回填狀況與筆者計算假設中的工況2接近。實驗按照上述負荷、流量的單U埋管連續運行9h，停機恢復15h。

3.2 計算結果對比

假設一個該實驗工況對應的數值計算條件，則此數值計算條件與數值計算條件相比［18］：1）兩者是同類現象；2）同屬非穩態對流換熱，各物理量隨時間變化趨勢相同，其他幾何條件、邊界條件及物理條件都分別成比例，故單值性條件相似已得到滿足；3）兩者同名的已定準則Nu＝f（Re，Pr）相等。故可得兩個現象的流動及換熱相似，可將該實驗的實測結果與計算結果相比較。

學生在掌握有一定寫作素材后，緊跟著需要學習有關作文方面的技巧。例如：作文開頭方式可以為欲揚先抑、開門見山、排比反復、設問開篇、倒敘開頭等；結尾技巧有凸顯主旨式、首尾呼應式、升華主題式等。技巧很多，也很零散，學生在寫作時不需要全部應用，只需要做到學以致用，在有些時候能為文章添加點綴。最后，學生要懂得謀篇布局，審題、立意、選材，一篇好作文即可完成。

實驗及數值計算的地埋管進、出水水溫如圖4所示。在實際運行的前4h，系統處于向穩定運行的過渡期，在運行4h后，實驗結果與工況2的數值計算結果相接近，兩者的誤差小于±7%。實驗進、出水最終溫度為32.5、30.5℃，回填工況2的進、出水溫度分別為34.3、30.8℃，兩者吻合的較好。

圖4 實驗數據與計算結果的埋管進、出水溫度比較

4 結論

1）在計算條件下，回填孔壁處的空氣間隙對土壤溫度恢復性能影響較大，降低了地埋管的換熱性能，與回填密實工況相比，回填孔壁及其以內巖土區域的初始過余溫度上升了約3.5℃，巖土平均過余溫度在恢復初始時上升了3.6℃，在恢復期的前6h明顯高于回填密實工況。與回填密實工況相比，進、出水管壁處的空氣間隙使得管壁處巖土的初始過余溫度上升了約8.5℃，管壁空氣間隙對運行時地埋管的進、出水溫度影響較大。建議采用間歇運行的模式來保證運行后較低的土壤過余溫度。

當然，產假制度也有其消極影響。一方面，對于企業來說，員工休產假會打亂企業平常的節奏，企業還要承擔部分因員工產假而導致的工作調動、繳納社保、工資等費用。另一方面，對于女職工自身而言，她們會因為休產假而長期脫離工作崗位，可能存在著影響個人晉升與職業長期發展的擔憂。公司和員工在產假這個問題上進行博弈，由于雇主擁有權力和信息方面的絕對優勢，可以通過釋放一些模棱兩可的信息和隱晦的威脅，從而對員工產生壓力并且在實際操作中變相降低她們的收入等。

2）數值計算結果與實驗結果相比較，實驗中單U地埋管的換熱性能與數值計算回填孔壁處存在空氣間隙的狀況下的換熱性能相接近，說明該數值計算方法可以應用于不同回填工況下的地埋管傳熱計算。

3）由于巖土溫度向巖土初始溫度的恢復速率不一致，前期恢復時間較短，而到接近初始地溫階段，恢復時間較長，故需要合理的恢復時間來保證地埋管的換熱能力。

［1］田慧峰，曹偉武.地埋管長度計算中關鍵參數的計算方法研究［J］.土木建筑與環境工程，2009，31（1）：110－>124.TIAN Huifeng，CAO Weiwu.ThEkey parametersmethod of calculating geothermal heat exchanger length［J］.Journal of Civil，Architectural ＆ Environmental Engineering，2009，31（1）：110－>124.

［2］Bhuttamm A B，Hayat N，Bashirm H，et al.CFD applications in various heat exchangers design：a review［J］.Applied Thermal Engineering，2012，32：1－>12.

［3］劉俊，張旭，高軍，等.地源熱泵土壤溫度恢復特性研究［J］.暖通空調，2008，38（11）：147－>150.LIU Jun，ZHANG Xun，GAO Jun，et al.Research of soil temperaturErestoration characteristics of groundsourcEheat pumPsystems［J］.Heating Ventilating ＆Air Conditioning，2008，38（11）：147－>150.

［4］Shang Y，Li S F，Li H J.Analysis of geo－>temperaturErecovery under intermittent operation of ground－>sourcEheat pump［J］.Energy and Buildings，2011，43：935－>943.

［5］趙利君，周亞素，張行洋，等.地源熱泵地埋管周圍土壤溫度恢復的模擬［J］.建筑節能，2010，38（10）：34－>38.ZHAO Lijun，ZHOU Yasu，ZHANG Xingyang，et al.Simulation on soil temperaturErecovery for ground sourcEheat pump［J］.Construction Conserves Energy，2010，38（10）：34－>38.

［6］高源，謝俊，董鵬博，等.土壤源熱泵間歇運行特性及地溫恢復規律研究［J］.節能，2010，341（12）：33－>36.GAO Yuan，XIEJun，DONG Pengbo，et al.ThEexperimental study on intermittentmovement and ground temperaturErestorativEcharacteristics of ground sourcEheat pump［J］.Energy Conservation，2010，341（12）：33－>36.

［7］陳穎，楊敏，史保新.土壤源熱泵制冷間歇工況土壤溫度響應實驗研究［J］.太陽能學報，2009，30（10）：1193－>1197.CHEN Ying，YANGmin，SHI Baoxin.Experimental investigation on soil pemperaturErestorativEcharacteristics for soil sourcEheat pumPin intermittent cooling operation ［J］.Acta Energy Solaris Sinica，2009，30（10）：1193－>1197.

［8］孟召賢，余躍進，郭友明.淺析回填材料在地源熱泵應用中的影響［J］.建筑節能，2010，38（3）：69－>71.MENG Zhaoxian， YU Yuejin， GUO Youming.InfluencEof thEbackfillmaterials on thEground－>sourceheat pump（Gshp）［J］.Construction Conserves Energy，2010，38（3）：69－>71.

［9］Zeng H Y，Diao N R，Fang Z H.Heat transfer analysis of boreholEin vertical ground heat exchangers［J］.International Journal of Heat andmass Transfer.2003，46：4467－>4481.

［10］Zhang H F，GEX S，YEH，et al.Heat conduction and heat storagEcharacteristics of soils ［J］. Applied Thermal Engineering，2007，27：369－>373.

［11］朱清宇，肖龍.U形管換熱器的回填工藝［J］.供熱制冷，2009，11：52.ZHU Qingyu，XIAO Long.ThEadd－>backprocess of U－>tubEground heat exchangE［J］. Heating ＆Refrigeration，2009，11：52.

［12］Wang Y，Wong kkL，Liu Q H，et al.Improvement of energy efficiency for an open－>looPsurfacEwater sourcEheat pumPsystem via optimal design of waterintake［J］.Energy and Buildings，2012，51：93－>100.

［13］Choi J C，LeES R，LeED S.Numerical simulation of vertical ground heat exchangers：intermittent operation in unsaturated soil conditions ［J］.Computers and Geotechnics，2011，38：949－>958.

［14］Florides G A，Christodoulides P，Pouloupatis P.An analysis of heat flow through a boreholEheat exchanger validatedmodel［J］.Applied Energy，2012，92：523－>533.

［15］王勇，唐曦.衛生熱水蓄熱方式對地理管換熱性能的影響分析［J］.重慶大學學報，2012，35（3）：140－>146.WANG Yong，TANG Xi.InfluencEof heat storagemodEof domestic hot water on thEperformancEof ground heat exchanger in ground－>sourcEheat pamPsystem［J］.Journal of Chongqing University，2012，35（3）：140－>146.

［16］唐曦，王勇.水平埋管地下巖土傳熱模型中上墊面邊界條件的確定方法研究［J］.制冷與空調，2011，107（25）：197－>201.TANG Xi，WANG Yong.Method of ground surfacEboundary condition of soil heat transfermodel around horizontal ground heat exchangers［J］.Refrigeration and Air Conditioning，2011，107（25）：197－>201.

［17］李新國，趙軍，周倩.U型垂直埋管換熱器管群周圍土壤溫度數值模擬［J］.太陽能學報，2004，25（5）：703－>707.LI Xinguo，ZHAO Jun，ZHOU Qian.Numerical simulation on thEground temperaturEfield around U－>pipEunderground heat exchangers［J］.Acta Energy Solaris Sinica，2004，25（5）：703－>707.

［18］張靖周，常海萍.傳熱學［M］.北京：科學出版社，2009.