用動態(tài)模擬方法設(shè)計(jì)Ｕ形地埋管換熱器

摘要:給出了Ｕ形地埋管換熱器的數(shù)值模型。利用實(shí)測數(shù)據(jù)對比驗(yàn)證了自主開發(fā)的模擬軟件的準(zhǔn)確性。在模擬建筑全年動態(tài)負(fù)荷的基礎(chǔ)上，采用該模擬軟件對某地源熱泵工程地埋管換熱器進(jìn)行了３０年的逐時模擬，確定了埋管方案。指出不能采用單位深度換熱量設(shè)計(jì)地源熱泵換熱器，必須進(jìn)行動態(tài)模擬。     關(guān)鍵詞：地源熱泵?。招蔚芈窆軗Q熱器　數(shù)值模型　動態(tài)模擬     ０　引言     地源熱泵空調(diào)在國內(nèi)外受到廣泛關(guān)注，發(fā)展迅速［１］。在我國地源熱泵推廣中，科學(xué)設(shè)計(jì)、系統(tǒng)優(yōu)化、規(guī)范施工及基礎(chǔ)研究方面仍需加強(qiáng)［１］。     地源熱泵系統(tǒng)設(shè)計(jì)必須進(jìn)行熱響應(yīng)實(shí)驗(yàn)、建筑動態(tài)負(fù)荷計(jì)算、熱泵動態(tài)負(fù)荷計(jì)算、系統(tǒng)方案初步擬定、地埋管換熱器水溫及土壤溫度的動態(tài)模擬等工作。地埋管換熱器設(shè)計(jì)合理與否是地源熱泵工程成功與否的關(guān)鍵，但因地埋管換熱器的非穩(wěn)態(tài)特性而使其設(shè)計(jì)復(fù)雜化。本文給出了地埋管換熱器的數(shù)值模型，在模擬建筑全年動態(tài)負(fù)荷基礎(chǔ)上，采用自主開發(fā)的Ｕ形地埋管換熱器動態(tài)數(shù)值模擬軟件，對某地源熱泵工程地埋管換熱器進(jìn)行了３０ａ的逐時模擬，確定了埋管方案。計(jì)算表明，地埋管換熱器動態(tài)數(shù)值模擬方法準(zhǔn)確性較高，快速，具有工程實(shí)用性。     １　熱響應(yīng)實(shí)驗(yàn)     部分工程簡單按經(jīng)驗(yàn)數(shù)據(jù)設(shè)計(jì)地源熱泵系統(tǒng)，但由于各種場地地質(zhì)水文條件的差異，經(jīng)驗(yàn)數(shù)據(jù)不具普適性，估算可能造成系統(tǒng)不合理。因此ＧＢ５０３６６—２００５《地源熱泵系統(tǒng)工程技術(shù)規(guī)范》［２－３］規(guī)定必須進(jìn)行熱響應(yīng)實(shí)驗(yàn)。     熱響應(yīng)實(shí)驗(yàn)的目的在于獲取試驗(yàn)孔單位深度吸放熱量、巖土的導(dǎo)熱系數(shù)及體積比熱容。必須注意，熱響應(yīng)實(shí)驗(yàn)獲得的單位深度換熱量只能作為參考進(jìn)行初步規(guī)劃，因?yàn)樵摂?shù)據(jù)是在特定溫差條件下測得的，并且熱響應(yīng)測試時間較短，不能有效反映熱量長期累積對換熱的影響。準(zhǔn)確的土壤導(dǎo)熱系數(shù)和體積比熱容才是系統(tǒng)設(shè)計(jì)模擬的基本依據(jù)。巖土熱物性參數(shù)可以在對數(shù)時間坐標(biāo)上用直線擬合方法計(jì)算，也可直接用多參數(shù)估計(jì)法計(jì)算［２－９］，其基本思想是土壤熱物性參數(shù)具代表性的取值應(yīng)該保證模型計(jì)算水溫與實(shí)測水溫的方差小。關(guān)于熱物性參數(shù)計(jì)算方法的探討，此處不予詳述。     ２　地埋管換熱器的動態(tài)模擬方法     建筑負(fù)荷計(jì)算可以采用常見的軟件。此處僅闡述地埋管換熱器的模擬方法。     ２．１　地埋管換熱器數(shù)值模型     目前國外地?zé)崮M軟件主要有ＧＬＨＥＰＲＯ，ＧＬＤ［１０］。國內(nèi)有地?zé)嶂擒浖?。在參考了現(xiàn)有計(jì)算模型［４－１４］的基礎(chǔ)上，筆者采用數(shù)值方法［１５－１７］獨(dú)立開發(fā)了Ｕ形地埋管換熱器動態(tài)模擬軟件［１８］。     該軟件采用圓柱熱源模型，鉆孔內(nèi)的傳熱簡化為穩(wěn)態(tài)傳熱，孔外傳熱認(rèn)為是非穩(wěn)態(tài)的。忽略軸向?qū)?，也忽略地面溫度波動和埋管底部傳熱的影響。離散網(wǎng)格如圖１所示。                               式（１）～（９）中?。簦?，τ為第ｉ個節(jié)點(diǎn)τ時刻的溫度，℃；Ｑ為某個時刻換熱器的熱負(fù)荷，Ｗ；Δτ為時間步長，ｓ；ρｓｃｓ為土壤的體積比熱容，Ｊ／（ｍ３·Ｋ）；Ｖｉ為控制單元的容積，ｍ３；λ為土壤的導(dǎo)熱系數(shù)，Ｗ／（ｍ·Ｋ）；Δｒ?yàn)榭臻g步長，ｍ；ｒｉ為第ｉ個節(jié)點(diǎn)對應(yīng)的半徑，ｍ；ｚ為整個換熱器豎向深度，ｍ。多孔布置時，孔間距的中心位置可以考慮為絕熱邊界，ｔＭ，τ計(jì)算式應(yīng)作相應(yīng)改變，此略。     流體節(jié)點(diǎn)與孔壁通過下式起來：                  式中?。簦鏋榱黧w平均溫度，℃；Ｒｂ為單位深度孔內(nèi)熱阻，ｍ·Ｋ／Ｗ，其計(jì)算方法見文獻(xiàn)［１１］；ｔ１為孔壁溫度，℃。     地埋管換熱器進(jìn)口溫度ｔｆ，ｉｎ和出口溫度ｔｆ，ｏｕｔ分別為：                   式（１１）～（１３）中　Δｔｆ為進(jìn)出口溫差，℃；ｍ為質(zhì)量流量，ｋｇ／ｓ；ｃｐ，ｗ為循環(huán)水的比定壓熱容，Ｊ／（ｋｇ·Ｋ）。在土壤熱物性參數(shù)、設(shè)計(jì)流量確定后，利用該模型可以設(shè)定換熱器入口溫度已知，計(jì)算出口溫度、土壤溫度及換熱量；也可以設(shè)定換熱量已知，計(jì)算進(jìn)、出口溫度和土壤溫度。     ２．２　動態(tài)模擬軟件驗(yàn)證及討論     為了驗(yàn)證自主開發(fā)軟件的準(zhǔn)確性，筆者采用２００９年４月對綿陽某場地進(jìn)行熱響應(yīng)實(shí)驗(yàn)得到的數(shù)據(jù)［１８－１９］進(jìn)行對比。方法如下：１）根據(jù)實(shí)測的進(jìn)、出水溫度和流量，用線熱源解析模型估計(jì)土壤導(dǎo)熱系數(shù)及體積比熱容。２）基于估計(jì)的導(dǎo)熱系數(shù)和體積比熱容、實(shí)測的流量，以實(shí)測入口溫度為已知條件，計(jì)算出口溫度。３）比較實(shí)測的出口溫度和模擬計(jì)算的出口溫度。在設(shè)定入口溫度的情況下，雙Ｕ形管出口溫度實(shí)測值和模擬計(jì)算值的比較如圖２～５所示。                   圖２，４表明，吸熱和放熱工況下，出口溫度模擬計(jì)算值與實(shí)測值趨勢一致，說明數(shù)值模型整體上符合物理規(guī)律。值得注意的是：熱響應(yīng)實(shí)驗(yàn)初期并不嚴(yán)格滿足恒熱流假設(shè)，恒熱流線熱源或柱熱源解析解一般用于實(shí)驗(yàn)進(jìn)行１０ｈ以后才比較準(zhǔn)確。該數(shù)值模型適用于變熱流工況，用于放熱／吸熱初期的模擬仍有較高的準(zhǔn)確度。     由圖３，５可以看出，模擬值總體偏高，但絕大部分情況下偏差在１℃以內(nèi)。偏差一方面來自于模型及算法的近似處理（比如空間步長、時間步長對精度有影響），另一方面有可能來自于實(shí)測數(shù)據(jù)本身（比如巖土熱物性參數(shù)包含有某種程度的不確定性）。筆者認(rèn)為除了空間步長、時間步長對精度有影響外，孔內(nèi)穩(wěn)態(tài)熱阻的計(jì)算也可能造成水溫的較大偏差。文獻(xiàn)［１９］只估算了導(dǎo)熱系數(shù)及體積比熱容２個參數(shù)，孔內(nèi)熱阻采用文獻(xiàn)［１１］的公式計(jì)算。筆者擬用熱響應(yīng)實(shí)驗(yàn)測試數(shù)據(jù)同時估算導(dǎo)熱系數(shù)、體積比熱容及孔內(nèi)熱阻３個參數(shù)。以此為基礎(chǔ)，用數(shù)值模型獲得的水溫模擬值可能會與實(shí)測值吻合得更好。     總體來說，自主開發(fā)的軟件具有較高的準(zhǔn)確性，可以用于工程設(shè)計(jì)。遺憾的是目前驗(yàn)證還于短期的熱響應(yīng)實(shí)驗(yàn)數(shù)據(jù)，長期模擬（１ａ以上）的準(zhǔn)確性還有待進(jìn)一步驗(yàn)證。長期模擬可以把軸向傳熱考慮進(jìn)去。該數(shù)值模型與其他計(jì)算模型的對比驗(yàn)證目前正在進(jìn)行。     除具有較高準(zhǔn)確性外，該軟件計(jì)算速度較快，動態(tài)模擬１ａ的運(yùn)行情況只需機(jī)時１ｍｉｎ，用于多種方案的分析篩選方便快捷，非常實(shí)用。     目前該軟件還未把地埋管換熱器與熱泵機(jī)組耦合起來進(jìn)行全系統(tǒng)能耗模擬，正在進(jìn)一步完善。全系統(tǒng)模擬需要給出熱泵機(jī)組的性能計(jì)算模型。     ３　地埋管換熱器設(shè)計(jì)舉例     ３．１　建筑負(fù)荷     某工程夏季空調(diào)設(shè)計(jì)總冷負(fù)荷為１?。担埃埃耄?，冬季空調(diào)設(shè)計(jì)總熱負(fù)荷為１　０００ｋＷ，全年逐時負(fù)荷如圖６所示。                    ３．２　冷熱源方案     采用２臺地埋管地源熱泵空調(diào)機(jī)組和１臺水冷螺桿式冷水機(jī)組。     熱負(fù)荷全部由２臺地源熱泵機(jī)組承擔(dān)。夏季冷負(fù)荷小于設(shè)計(jì)負(fù)荷的１／３時，開啟１臺冷水機(jī)組，以便土壤恢復(fù)冷熱平衡；冷負(fù)荷大于設(shè)計(jì)負(fù)荷的１／３時，開啟１臺冷水機(jī)組承擔(dān)設(shè)計(jì)負(fù)荷的１／３，其余部分由１臺或２臺地源熱泵機(jī)組承擔(dān)。     ３．３　地源熱泵機(jī)組負(fù)荷     根據(jù)冷熱源方案，得到熱泵機(jī)組的全年負(fù)荷，如圖７所示。                   ３．４　地埋管換熱器負(fù)荷     根據(jù)熱泵機(jī)組的ＥＥＲ和ＣＯＰ（目前沒考慮機(jī)組性能的動態(tài)變化），計(jì)算得到地埋管換熱器的全年負(fù)荷，如圖８所示。全年累計(jì)排熱量為７．０８×１０８　ｋＷ·ｈ，累計(jì)吸熱量為９．７４×１０８　ｋＷ·ｈ。全年冷熱負(fù)荷比較接近。                   ３．５　換熱器初步方案     根據(jù)當(dāng)?shù)氐刭|(zhì)條件及打孔費(fèi)用，建議鉆孔深度為７０ｍ。     圖８中，夏天地埋管換熱器承擔(dān)的大放熱負(fù)荷為９２４ｋＷ，按單位深度放熱量７０Ｗ／ｍ估計(jì)，鉆孔數(shù)為１９０個。圖８中，冬天大吸熱負(fù)荷為８０６ｋＷ，按單位深度吸熱量６０Ｗ／ｍ估計(jì)，鉆孔數(shù)為１９２個。綜合以上數(shù)據(jù)，確定初步方案：鉆２００孔。     ３．６　動態(tài)模擬及優(yōu)化     按初步方案（２００孔）全年運(yùn)行時換熱器的進(jìn)、出水溫度如圖９所示?？梢钥闯?，２００孔方案出水水溫偏高（高近４０℃，高于３２℃的時間達(dá)７４ｈ），不符合熱泵冷水機(jī)組運(yùn)行條件，說明按單位深度吸、放熱量經(jīng)驗(yàn)數(shù)據(jù)進(jìn)行設(shè)計(jì)是不可靠的。                   考慮到熱量累積的影響，打孔數(shù)增加為２６０個，全年運(yùn)行時換熱器進(jìn)、出水溫度如圖１０所示。２６０孔方案下，出水溫度超過３２℃的時間不足１８ｈ，低在５℃以上?？梢哉J(rèn)為全年運(yùn)行時該系統(tǒng)會有較高的效率。              ２６０孔方案下，孔壁溫度高時刻及低時刻土壤溫度的空間分布如圖１１所示。                   ２６０孔方案下，第３０年鉆孔周圍半徑１０ｍ內(nèi)土壤溫度分布如圖１２所示，其中孔壁及距鉆孔中心３，４ｍ處的土壤溫度全年變化如圖１３所示。從圖１１～１３可以看出：     １）土壤高溫度分布，運(yùn)行３０ａ后幾乎沒有變化?？妆诟邷囟燃s為２６℃。     ２）土壤低溫度分布，運(yùn)行３０ａ后有輕微降低。因?yàn)樵摲桨钢?，累?jì)吸熱比累計(jì)放熱大?？妆诘蜏囟燃s為１４．５℃     ３）鉆孔周圍溫度變化明顯的區(qū)域在３ｍ以內(nèi)。因此建議孔間距為６ｍ。     ４　結(jié)語     本文給出了地埋管換熱器非穩(wěn)態(tài)傳熱數(shù)值模型。該動態(tài)模型計(jì)算準(zhǔn)確、快捷，可用于工程設(shè)計(jì)。計(jì)算實(shí)例表明，不能采用單位深度換熱量設(shè)計(jì)地埋管換熱器的終方案，必須進(jìn)行地埋管換熱器動態(tài)模擬，確定換熱器方案，保證機(jī)組運(yùn)行的條件。本文只模擬了地埋管換熱器，地埋管換熱器與熱泵機(jī)組的耦合未考慮。應(yīng)對該方法進(jìn)行進(jìn)一步的完善，以便耦合模擬全系統(tǒng)的能耗。  

產(chǎn)品咨詢請北京鴻鷗儀器（bjhocy），產(chǎn)品搜索：地源熱泵測溫，地埋管測溫
　　
　　關(guān)鍵詞：地源熱泵地埋管溫度測量系統(tǒng)實(shí)現(xiàn)實(shí)時溫度在線監(jiān)測/地源熱泵換熱井實(shí)時溫度電腦監(jiān)測系統(tǒng)/GPRS式豎直地埋管地源熱泵溫度監(jiān)控系統(tǒng)/地源熱泵溫度場測控系統(tǒng)/地埋管測溫/地源熱泵溫度監(jiān)控/地源熱泵測溫
　　
　　遠(yuǎn)程全自動地溫監(jiān)測系統(tǒng)/鐵路凍土地溫監(jiān)測系統(tǒng)/地溫監(jiān)測系統(tǒng)/城市地溫監(jiān)測自動化系統(tǒng)/礦井深部地溫/地源熱泵監(jiān)測研究/地源熱泵溫度測量系統(tǒng)/淺層地?zé)釡y溫/深水測溫儀/深井測溫儀/深水測溫儀/深井測溫儀
　　
　　推薦產(chǎn)品如下：
　　
　　地源熱泵溫度監(jiān)控系統(tǒng)/地源熱泵測溫/淺層地溫能動態(tài)監(jiān)測系統(tǒng)