余曉平，陳杰，周仁君，劉雨

（重慶科技學院 建筑工程學院，重慶 401331）

0 引言

地源熱泵系統(tǒng)是以巖土體、地下水或地表水為低溫熱源，由熱泵機組、室外換熱系統(tǒng)和建筑物內(nèi)供熱空調(diào)系統(tǒng)組成的系統(tǒng)。 地源熱泵技術屬于可再生能源利用技術，近年來得到大力推廣。在地源熱泵技術的應用過程中，獲取巖石、土壤的熱物性參數(shù)對地源熱泵系統(tǒng)設計起著至關重要的作用， 其值一般通過熱響應測試試驗獲取。

工程現(xiàn)場的巖土熱響應試驗，測試環(huán)境較復雜，測試周期長[1]，環(huán)境影響因素較多，具體現(xiàn)場測試環(huán)境如圖1、圖2 所示，在該環(huán)境下主要還存在以下問題：

圖1 測試現(xiàn)場環(huán)境

圖2 測試現(xiàn)場設備

（1） 測試地較偏遠，工程現(xiàn)場的電壓波動較大， 影響設備的工作性能[2]，導致地埋管內(nèi)溫度、流量數(shù)據(jù)異常，影響測試結(jié)果的準確性；

（2） 人工進行熱響應測試穩(wěn)態(tài)判定易出錯，影響測試結(jié)果準確性；

（3） 測試設備的控制手段落后，沒有集成的系統(tǒng)控制終端[3]，測試數(shù)據(jù)的后期處理重復工作量大[4]。

因此，筆者研發(fā)一種適用于工程現(xiàn)場環(huán)境需求的熱響應測試系統(tǒng)，既可提高測試自動化程度，又可用于學校虛擬仿真實驗，為建環(huán)專業(yè)學生提供開放實驗訓練。 本文通過測試系統(tǒng)軟件模塊搭建和軟件編程， 研究了系統(tǒng)如何實現(xiàn)數(shù)據(jù)采集、數(shù)據(jù)處理和設備控制，分析了工程環(huán)境對測試的影響因素，可為同類型的工程測試系統(tǒng)開發(fā)和虛擬仿真實驗臺建設提供參考。

1 基于MCGS 熱響應測試系統(tǒng)搭建

本系統(tǒng)的搭建流程主要分為窗口界面設計、通過Modbus 模塊實現(xiàn)數(shù)據(jù)采集和設備控制以及Matlab 處理測試數(shù)據(jù)幾個步驟，具體如圖3 所示。

圖3 系統(tǒng)搭建流程圖

1.1 窗口界面設計

本系統(tǒng)基于監(jiān)視與控制通用系統(tǒng)Monitor and Control Generated System（簡稱MCGS）研發(fā)，根據(jù)系統(tǒng)需要實現(xiàn)數(shù)據(jù)監(jiān)測和設備控制功能，設置4 個運行窗口，分別是系統(tǒng)基本測試參數(shù)界面、測試參數(shù)顯示界面、控制面板界面和數(shù)據(jù)歷史曲線界面，具體如圖4—圖7 所示。

圖4 系統(tǒng)基本測試參數(shù)界面

圖5 測試參數(shù)顯示界面

圖6 控制面板界面

圖7 數(shù)據(jù)歷史界面

1.2 通過Modbus 模塊實現(xiàn)數(shù)據(jù)采集及設備控制

Modbus 通訊協(xié)議類似于人與人之間交流的語言， 它是用于電子控制器上的一種通用語言，能夠?qū)崿F(xiàn)設備與設備之間的數(shù)據(jù)傳輸。 本系統(tǒng)主要采用支持Modbus 通訊協(xié)議的熱電偶采集器和多通道采集器采集溫度、流量、電壓和電流數(shù)據(jù)，采用支持Modbus 通訊協(xié)議的工業(yè)繼電器控制板對水泵以及加熱器進行啟?？刂?。RS485 接口是Modbus 通訊協(xié)議的標準接口，在工控機上一般沒有此接口，因此需要RS485 轉(zhuǎn)USB 裝置實現(xiàn)工控機與Modbus工業(yè)模塊的數(shù)據(jù)傳輸。 在建立工控機與各個Modbus 工業(yè)模塊通訊時，對不同的采集數(shù)據(jù)以及控制變量，需在MCGS 軟件中設置相應的數(shù)據(jù)變量和數(shù)據(jù)傳輸通道， 并在軟件和Modbus 工業(yè)模塊中設置對應的地址、波特率以及校驗方式，才能進行正常的數(shù)據(jù)傳輸。 本系統(tǒng)的控制原理如圖8 所示。

圖8 系統(tǒng)控制原理圖

1.3 利用Matlab 處理熱響應測試數(shù)據(jù)

巖土的熱物性參數(shù)一般采用基于線熱源模型的斜率法求得[5]，此過程需要對大量的溫度數(shù)據(jù)進行處理，再經(jīng)過復雜的數(shù)學公式計算得到結(jié)果。 MCGS 本身不具備處理大量數(shù)據(jù)以及進行復雜數(shù)學計算的能力， 因此本文提出采用Matlab進行巖土熱物性參數(shù)的計算， 再將計算結(jié)果返回到MCGS中，為熱響應測試復雜的數(shù)據(jù)處理問題提供一種解決方案。

Matlab 的數(shù)學函數(shù)庫是數(shù)學算法的一個巨大集合[6]。 在Matlab 中以線熱源模型為基礎進行數(shù)學模型的建立，首先根據(jù)線熱源模型將測試過程中的進出水平均溫度Tf與時間對數(shù)ln(τ)運用polyfit 函數(shù)擬合成y=ax+b 的線性關系式，相關代碼如下：

k=60; %溫度采集時間間隔，S

T=k∶k∶n*k; %設置溫度采集時間范圍，n 表示測試數(shù)據(jù)的組數(shù)

y1=polyfit(log(T),tf,1); %對溫度采集的時間對數(shù)以及進出水平均溫度擬合

a=y1(∶,1); %得到擬合函數(shù)的斜率

b=y1(∶,2); %得到擬合函數(shù)的截距

得到斜率a 以及截距b，再經(jīng)過線熱源模型的理論公式（1）

轉(zhuǎn)換得到導熱系數(shù)以及容積比熱容的具體計算表達式，見公式（2）及公式（3）[7-8]。

式中，Q——地埋管換熱量，W；H——地埋管深度，m；rb——鉆孔半徑，m；γ——歐拉常數(shù)，取0.5；Rb——管內(nèi)流體到鉆孔壁熱阻，m.℃/w；TO——初始地溫，℃。

Rb需要導入不同溫度水的熱物性參數(shù)表， 根據(jù)不同的水溫選取參數(shù)值計算。根據(jù)式（1）—式（3）的計算表達式，編寫程序最后得到基于線熱源的數(shù)學計算模型。 通過OPC[9]建立MCGS 與Matlab 之間的連接，Matlab 通過OPC 讀取MCGS 的基本測試數(shù)據(jù)進行計算， 計算完成后再將熱物性參數(shù)的值通過OPC 傳回MCGS。

采用Matlab 建立線熱源模型處理實驗數(shù)據(jù)不僅能夠快速得到結(jié)果，同時Matlab 自帶的擬合工具箱還能夠繪制出擬合圖像以便于分析，通過調(diào)整擬合度來提高擬合結(jié)果的準確性。

2 影響測試結(jié)果的因素分析

2.1 設備運行參數(shù)影響

工程現(xiàn)場熱響應測試常因電壓不穩(wěn)定導致測試中斷，因地埋管水溫上升導致地埋管流速變化，從而對測試結(jié)果產(chǎn)生影響。 根據(jù)《地源熱泵系統(tǒng)工程技術規(guī)范》[10]對地埋管換熱器內(nèi)流速和溫度設定的要求，當實測數(shù)據(jù)超出設定范圍時，系統(tǒng)應能夠?qū)Ξ惓?shù)據(jù)作出及時反饋并控制設備以及時調(diào)節(jié)。

通過對系統(tǒng)設定報警顯示， 能夠有效解決上述問題。 在MCGS 軟件中， 對需要設定范圍的數(shù)據(jù)設置上下限報警， 如圖9所示； 在控制策略中編寫判定條件以及反饋控制的腳本程序，如圖10 所示。 當流速超出規(guī)定范圍時，系統(tǒng)界面會進行報警顯示，并自動調(diào)節(jié)水泵的變頻器，使流速達到規(guī)定范圍，與傳統(tǒng)的PID控制相比方法更加簡單，方便測試人員的日常使用。

圖9 異常報警設置

圖10 設置判定條件及反饋控制

2.2 穩(wěn)態(tài)判定影響

熱響應測試包括初溫測試以及巖土換熱兩個階段。 巖土初溫測試一般采用無功循環(huán)法，通過地埋管內(nèi)水循環(huán)，記錄循環(huán)流體進出口的溫度隨時間的變化， 當出口溫度達到穩(wěn)定時地埋管的進出口平均溫度值即是巖土的初始溫度。 巖土換熱階段持續(xù)時間較長，一般在48h 以上，只有當測試達到換熱平衡時才能進行巖土換熱系數(shù)以及容積比熱容的計算，《地源熱泵系統(tǒng)工程技術規(guī)范》[10]中對換熱平衡時的條件作了詳細規(guī)定。 熱響應測試初溫測試階段以及巖土換熱階段，都存在一個穩(wěn)態(tài)判定過程。目前在工程現(xiàn)場， 此工作基本都是由人工完成， 存在因測試數(shù)據(jù)量大，主觀判定容易出錯的問題。

MCGS 提供的事件策略功能，可以通過設定條件來完成系統(tǒng)的自動判穩(wěn)。 根據(jù)《地源熱泵系統(tǒng)工程技術規(guī)范》中提到的要求，在MCGS 中編寫判定初溫測試以及巖土換熱階段達到穩(wěn)態(tài)的程序。 以下是本系統(tǒng)的初溫穩(wěn)態(tài)判定程序：

If (地源側(cè)進口溫度-地源側(cè)出口溫度)〈0.1 Then

初始地溫判定計數(shù)器=初始地溫判定計數(shù)器+1

If 初始地溫判定計數(shù)器〉43200 Then %判定進出水溫差連續(xù)12h 變化不超過0.1℃

初始地溫=（地源側(cè)進口溫度+地源側(cè)出口溫度）/2

EndIf

EndIf

在不同測試階段， 當運行參數(shù)滿足各階段程序設定的條件時，系統(tǒng)就會自動得出土壤初溫和計算巖土換熱系數(shù)以及容積比熱容，實現(xiàn)從測試到計算輸出的自動化。

2.3 系統(tǒng)自動測試的控制流程

通過集成的控制終端能夠?qū)崿F(xiàn)熱響應測試從傳統(tǒng)的手動控制到設備自動控制的優(yōu)化，提升了設備的自動化程度，減少了工作人員的工作量。 本系統(tǒng)的自動控制流程如圖11 所示。

圖11 系統(tǒng)自動控制流程

3 案例應用分析

表1 為某項目豎直雙U 地埋管測試試驗孔埋管換熱器的基本參數(shù)。

表1 某項目試驗孔基本測試參數(shù)

將實測參數(shù)數(shù)據(jù)輸入到軟件中，得到時間與進出水平均溫度曲線以及擬合的線性曲線，如圖12 所示。

圖12 進出水平均溫度與時間對數(shù)曲線

通過對不同時間的數(shù)據(jù)段進行擬合，得到測試不同時長的計算結(jié)果，具體計算結(jié)果如表2 所示。

表2 系統(tǒng)測試與人工測試計算結(jié)果

從系統(tǒng)測試的結(jié)果看， 測試時長20h 與測試時長49.3h 導熱系數(shù)相差4%,容積比熱容相差6.3%，差異并不明顯，且導熱系數(shù)隨時間增大而增大，容積比熱容隨時間增大而減小。 理論上隨著測試時間增大，在恒熱流條件下進出水溫差也隨之增大，換熱量增大，所以導熱系數(shù)也隨之增大，當換熱量增大時，由式（2）和式（3）可知，容積比熱容減小，可見該軟件計算結(jié)果符合理論公式變化的規(guī)律，具有工程應用的可靠性。在49.3h 系統(tǒng)測試的結(jié)果與人工測試的結(jié)果對比中，導熱系數(shù)相差5.1%，容積比熱容相差9%，可見采用該系統(tǒng)測試與人工測試，結(jié)果具有一致性，系統(tǒng)能夠代替人工測試的結(jié)果，減少測試人員的工作量。

4 結(jié)論

（1） 建立了一種適用于工程現(xiàn)場的熱響應測試系統(tǒng)，能夠?qū)犴憫獪y試設備集成控制， 對測試數(shù)據(jù)實時記錄并自動計算，說明了該系統(tǒng)的搭建路徑以及實驗數(shù)據(jù)的處理方法。

（2） 分析了工程現(xiàn)場測試環(huán)境影響測試結(jié)果的因素，給出了通過設置運行策略來降低外界影響的解決方案。

（3） 通過案例驗證，系統(tǒng)測試時長在20h 以及49.3h 計算得到的導熱系數(shù)以及容積比熱容分別相差4%和6.3%，對于工程現(xiàn)場熱響應測試在不要求準確精度的情況下，合理的穩(wěn)態(tài)判定能夠有效縮短測試時間。 對49.3h 系統(tǒng)測試結(jié)果與人工測試結(jié)果進行對比可知，導熱系數(shù)相差5.1%，容積比熱容相差9%，結(jié)果具有一致性，可見系統(tǒng)測試能夠代替人工測試的結(jié)果，減少測試人員的工作量。

本文通過搭建熱響應測試系統(tǒng)，實現(xiàn)了設備各模塊之間的集成化和自動化，與人工測試相比降低了環(huán)境因素對測試結(jié)果的影響程度。 本系統(tǒng)后續(xù)進一步的優(yōu)化方向?qū)⑹羌友b無線通訊模塊，實現(xiàn)設備的遠程控制以及數(shù)據(jù)的遠傳。