空調(diào)水系統(tǒng)運(yùn)行參數(shù)關(guān)聯(lián)準(zhǔn)則研究

王玉杰 鄧亞宏 伍學(xué)智 徐 侃 徐新華

（1.華中科技大學(xué)建筑環(huán)境與能源應(yīng)用工程系 武漢 430074；2.中建三局第二建設(shè)工程有限責(zé)任公司 武漢 430000）

0 引言

2030年碳達(dá)峰和2060年碳中和是我國重大發(fā)展戰(zhàn)略。全球建筑能源消耗占總能源消耗的三分之一以上，并且是二氧化碳排放量的重要組成部分[1]。2018年，中國民用建筑建造能耗為5.2 億噸標(biāo)準(zhǔn)煤當(dāng)量，與民用建筑建造相關(guān)的碳排放總量約為18 億噸CO2[2]。在建筑能源消耗中，用于空調(diào)系統(tǒng)的運(yùn)行通常占比30%～40%，空調(diào)系統(tǒng)的運(yùn)行節(jié)能是建筑節(jié)能工作的一個重點[3,4]。在中央空調(diào)系統(tǒng)設(shè)計時，由于設(shè)計工況和余量的考慮，設(shè)計人員往往會選擇型號偏大的制冷主機(jī)、水泵和冷卻塔等，導(dǎo)致空調(diào)系統(tǒng)長期處于最大負(fù)荷的65%以下運(yùn)行[5,6]，系統(tǒng)中的水泵風(fēng)機(jī)等仍然以工頻或者固定頻率運(yùn)行，造成大量的能源浪費(fèi)[7]。

數(shù)據(jù)挖掘是一種新興的人工智能技術(shù)，它提供了處理海量且復(fù)雜數(shù)據(jù)的新方法[8]。數(shù)據(jù)挖掘是一個多學(xué)科的學(xué)科，集成了統(tǒng)計、機(jī)器學(xué)習(xí)、人工智能和高性能計算等技術(shù)。數(shù)據(jù)挖掘技術(shù)在建筑領(lǐng)域的預(yù)測[9,10]、故障檢測[10,11]和控制優(yōu)化[11,12]有很多研究與應(yīng)用。徐欣與田喆[7]利用多元線性回歸的方法預(yù)測日平均負(fù)荷。實際案例表明，采用該方法的預(yù)測精度可控制在15%以內(nèi)。Fan 與Xiao 等[13]研究了深度學(xué)習(xí)的方法在預(yù)測建筑物未來24 小時內(nèi)冷負(fù)荷的應(yīng)用，結(jié)果表明，使用無監(jiān)督深度學(xué)習(xí)提取的特征變量作為負(fù)荷預(yù)測的輸入時，可以明顯提高預(yù)測的準(zhǔn)確度。閆軍威等[14]根據(jù)空調(diào)實際運(yùn)行數(shù)據(jù)，利用最小二乘方法對各模型的參數(shù)進(jìn)行辨識，以冷源系統(tǒng)運(yùn)行設(shè)備總能耗最小為目標(biāo)，探討基于遺傳算法的中央空調(diào)冷源系統(tǒng)最佳運(yùn)行參數(shù)值，并將優(yōu)化的運(yùn)行參數(shù)用于控制進(jìn)行實驗驗證。數(shù)據(jù)挖掘可分為監(jiān)督分析和無監(jiān)督分析[15]。監(jiān)督分析需要在很多先驗知識的基礎(chǔ)上進(jìn)行研究分析。例如，在開發(fā)冷水機(jī)組能耗預(yù)測模型時，預(yù)先選擇了模型輸入，比如冷凍水的進(jìn)水溫度和回水溫度以及冷卻水的進(jìn)水溫度和回水溫度，領(lǐng)域的專業(yè)知識告訴我們，這些變量是影響冷水機(jī)組能耗的重要變量[16]。無監(jiān)督分析可以在沒有先驗知識或者先驗知識有限的情況下發(fā)現(xiàn)新穎的知識。無監(jiān)督分析專注于探索數(shù)據(jù)的內(nèi)在結(jié)構(gòu)、關(guān)聯(lián)和模式，因此具有發(fā)現(xiàn)潛在有用但先前未知的知識的能力。更重要的是，實施無監(jiān)督分析的成功較少取決于領(lǐng)域的專業(yè)知識且不受高質(zhì)量標(biāo)記培訓(xùn)數(shù)據(jù)的可用性的影響。

采用數(shù)據(jù)挖掘技術(shù)從大量的空調(diào)系統(tǒng)日常運(yùn)行數(shù)據(jù)中挖掘能夠指導(dǎo)空調(diào)系統(tǒng)高效運(yùn)行的參數(shù)，能夠為空調(diào)系統(tǒng)帶來可觀的節(jié)能效果。對中央空調(diào)水系統(tǒng)進(jìn)行參數(shù)優(yōu)化控制研究，提出在不同的邊界條件下，適合中央空調(diào)水系統(tǒng)高效穩(wěn)定運(yùn)行的狀態(tài)參數(shù)，對于空調(diào)系統(tǒng)高效運(yùn)行、降低運(yùn)行能耗和延長使用壽命等有著重要的意義。本文提出采用關(guān)聯(lián)規(guī)則挖掘的方法，以實際空調(diào)水系統(tǒng)為對象，進(jìn)行數(shù)據(jù)挖掘，獲得空調(diào)系統(tǒng)的運(yùn)行參數(shù)與相關(guān)邊界條件的關(guān)聯(lián)準(zhǔn)則，并在模擬平臺上進(jìn)行應(yīng)用分析，節(jié)能效果明顯。

1 關(guān)聯(lián)規(guī)則數(shù)據(jù)挖掘介紹

關(guān)聯(lián)規(guī)則挖掘是數(shù)據(jù)挖掘中常用的一種方法，是發(fā)現(xiàn)數(shù)據(jù)中頻繁模式的一種重要手段[17]。它在諸如購物籃、計算機(jī)網(wǎng)絡(luò)、推薦系統(tǒng)和醫(yī)療保健等眾多應(yīng)用中都取得了巨大的成功。關(guān)聯(lián)規(guī)則的基本概念如下：

交易數(shù)據(jù)庫（D）由n個不同的項目和m個交易組成。是一組n個項目，是一組有限的交易。事務(wù)t是集合I中的一組項目（t?I）?；谥С趾椭眯哦榷攘?，提取關(guān)聯(lián)規(guī)則X→Y。規(guī)則X→Y的支持度（用α表示）是同時包含X和Y的交易百分比，其中X和Y分別是規(guī)則的左側(cè)和右側(cè)項集。支持度計算為|X∪Y|/m，其中m是數(shù)據(jù)庫中的事務(wù)數(shù)，而|X∪Y|是規(guī)則的頻率。置信度是包含X的事務(wù)也包含Y事務(wù)的比例，用β表示。置信度的計算方式為|X∪Y|/|X|，其中|X∪Y|和|X|分別是規(guī)則和項目集X的支持計數(shù)。如果某項集的支持高于最小支持閾值，則該項集很頻繁。當(dāng)規(guī)則的置信度高于由βmin決定的最低置信度閾值時，規(guī)則就很強(qiáng)。

Apriori 算法是常用的一種關(guān)聯(lián)準(zhǔn)則的數(shù)據(jù)挖掘方法[13]。該算法使用頻繁項集性質(zhì)。大致思想就是逐層搜索，首先，它會掃描數(shù)據(jù)庫里面的所有數(shù)據(jù)，計算出每一種類型的數(shù)據(jù)個數(shù)，保留那些滿足最小支持度的數(shù)據(jù)類別。其次，找到第一個頻繁項集，然后從剩余項集里找出第二個頻繁項集。以此類推，不斷尋找，直到找不出來新的頻繁項集為止。提升度是一個表示規(guī)則挖掘是由有效的參數(shù)，定義為：

lift(X→Y)=P(Y/X)/P(X)

該定義表示在含有X的條件下，同時含有Y的概率，與Y同時發(fā)生的概率之比。滿足最小支持度和最小置信度的規(guī)則叫做強(qiáng)關(guān)聯(lián)準(zhǔn)則，但是強(qiáng)關(guān)聯(lián)準(zhǔn)則也是有有效和無效之分的。如果提升度lift(X→Y)＞1，那么規(guī)則就是有效的規(guī)則；如果提升度lift(X→Y)＜1，那么規(guī)則就是無效的規(guī)則；特別的，當(dāng)lift=1 的時候，代表X和Y是相互獨(dú)立的。

2 空調(diào)水系統(tǒng)及模擬分析

2.1 空調(diào)系統(tǒng)介紹

本文以中央空調(diào)水系統(tǒng)為研究對象，如圖1所示。該系統(tǒng)為某地下空間提供冷量。該系統(tǒng)通過冷卻塔制取冷卻水，冷凍水由制冷機(jī)制取，末端為空氣處理機(jī)組，有兩臺冷凍水泵和兩臺冷卻水泵。制冷機(jī)最小運(yùn)行負(fù)載為15%。系統(tǒng)中的水泵與風(fēng)機(jī)全部采用變頻控制。制冷機(jī)額定制冷量為534.3kW，蒸發(fā)器額定流量為92m3/h，冷凝器額定流量為110m3/h，額定功率為104.9kW。冷凍水泵額定流量為171m3/h，功率為22kW，揚(yáng)程為32.7mH2O。冷卻水泵額定流量為163.5m3/h，功率為18.5kW，揚(yáng)程為29.3m。冷卻塔額定冷卻水流量150m3/h，風(fēng)機(jī)功率5.5kW，冷卻水進(jìn)出水的設(shè)計溫度37/32℃。水系統(tǒng)設(shè)置A、B 支路為大空調(diào)機(jī)組提供冷凍水，C、D 為小空調(diào)機(jī)組提供冷凍水。

圖1 空調(diào)水系統(tǒng)示意圖Fig.1 Schematica of air conditioning water system

2.2 系統(tǒng)模擬及標(biāo)定

本文在TRNSYS 軟件中對建筑的負(fù)荷與空調(diào)系統(tǒng)進(jìn)行了建模，通過制冷機(jī)實際運(yùn)行參數(shù)對制冷機(jī)模型進(jìn)行驗證，采用實際測量數(shù)據(jù)對冷卻塔、空氣處理機(jī)組（AHU）和水泵等進(jìn)行模型驗證。本文采用的冷機(jī)模型[18]考慮了機(jī)組內(nèi)部參數(shù)，能夠較為準(zhǔn)確地模擬機(jī)組的運(yùn)行。冷卻塔模型采用傳熱單元法（NTU）根據(jù)進(jìn)口參數(shù)計算出口空氣及冷卻水出水狀態(tài)。該方法先假定出口空氣濕球溫度，利用熱交換器換熱原理采用迭達(dá)的方法進(jìn)行出口狀態(tài)的計算，在此過程中，假定劉易斯數(shù)為1，同時忽略冷卻塔側(cè)壁的傳熱。本文所用的AHU 模型是TRNSYS 中自帶的組件Type32。圖2 是制冷機(jī)的模型預(yù)測與實際測量功率對比，誤差在±5%之間。圖3 是冷卻塔預(yù)測模型與實際測量出水溫度對比，誤差范圍在+2%～-5%之間。

圖2 模型與實際測量的制冷機(jī)功率對比Fig.2 Comparison of the model predicted chiller power and the measurement

圖3 冷卻塔模型預(yù)測與實際測量的冷水出水溫度對比Fig.3 Comparison of the cooling water outlet temperatures of the cooling tower model prediction and the measurement

2.3 挖掘數(shù)據(jù)樣本

將空調(diào)水系統(tǒng)模型與建筑物負(fù)荷模型耦合在一起，空調(diào)水系統(tǒng)模型主要包括制冷機(jī)、冷水泵、冷卻塔集管網(wǎng)等。在模型中，制冷機(jī)輸出的冷凍水通過冷凍水泵輸送至AHU 中經(jīng)過換熱后返回到制冷機(jī)。從AHU 中出來的低溫空氣送至室內(nèi)空間，用PID 控制器控制冷凍水泵的頻率來控制冷凍水流量，進(jìn)而控制室內(nèi)空間的回風(fēng)溫度，使得回風(fēng)溫度維持在28℃左右，送風(fēng)量恒定。在模擬平臺運(yùn)行過程中，分別給定不同的制冷機(jī)冷凍水出水溫度設(shè)定值，從7℃至18℃，以1℃為間隔。在空調(diào)系統(tǒng)的冷卻水側(cè)，制冷機(jī)輸出的冷卻水經(jīng)過冷卻塔散熱返回到制冷機(jī)中，分別給定不同的冷卻水泵運(yùn)行頻率，從24Hz 至50Hz，每2Hz 為間隔。通過上述的運(yùn)行策略，可以得到大量的空調(diào)水系統(tǒng)運(yùn)行數(shù)據(jù)，主要關(guān)注數(shù)據(jù)為制冷機(jī)冷凍水出水溫度、冷凍水流量和回水溫度，冷卻水流量、冷卻水供回水溫度，各個設(shè)備的能耗與系統(tǒng)能效等。圖4 是TRNSYS 中的建筑與空調(diào)系統(tǒng)的耦合模擬平臺模型。圖中上半部分是建筑物的負(fù)荷模型，下半部分是空調(diào)水系統(tǒng)模型，左側(cè)是冷卻側(cè)，右側(cè)是冷凍側(cè)，中間是制冷機(jī)。

圖4 TRNSYS 中的建筑物與空調(diào)系統(tǒng)模型Fig.4 Models of the building and the air conditioning system in the TRNSYS platform

3 系統(tǒng)運(yùn)行參數(shù)關(guān)聯(lián)規(guī)則挖掘

本文選取了與空調(diào)水系統(tǒng)運(yùn)行有關(guān)的水系統(tǒng)運(yùn)行參數(shù)，包括冷凍水供水溫度、冷凍水溫差、冷卻水回水溫度、冷卻水溫差，以及室外空氣溫度、濕度與系統(tǒng)需求冷量為影響因素。將模擬得到的數(shù)據(jù)采用Apriori 算法進(jìn)行，獲得在不同的室外空氣溫濕度和系統(tǒng)冷量需求的邊界條件下，能夠讓空調(diào)水系統(tǒng)高效運(yùn)行的參數(shù)，即（冷凍水供水溫度、冷凍水溫差、冷卻水回水溫度、冷卻水溫差）→（水系統(tǒng)COP）樣式的關(guān)聯(lián)準(zhǔn)則，表1 是部分工況下的關(guān)聯(lián)準(zhǔn)則。

表1 關(guān)聯(lián)準(zhǔn)則表Table 1 Associated rules

4 關(guān)聯(lián)規(guī)則應(yīng)用及結(jié)果分析

本節(jié)對挖掘得到的關(guān)聯(lián)準(zhǔn)則的有效性和在中央空調(diào)水系統(tǒng)運(yùn)用的節(jié)能性進(jìn)行分析。首先按照實際的運(yùn)行模式進(jìn)行模擬，得到整個制冷季按照既有模式運(yùn)行的數(shù)據(jù)。既有模式即為冷凍水泵、冷卻水泵和冷卻塔風(fēng)機(jī)等工頻運(yùn)行，主機(jī)供水溫度恒定，不隨冷量需求變化而調(diào)整。關(guān)聯(lián)準(zhǔn)則運(yùn)行模式則是通過讀取建筑物冷量需求、空氣溫度與濕度數(shù)據(jù)，模型平臺中的關(guān)聯(lián)準(zhǔn)則應(yīng)用策略模塊（見圖5 左下角）經(jīng)過判斷可以尋找該邊界條件下的最優(yōu)參數(shù)值（即目標(biāo)值），即冷凍水供水溫度設(shè)定值、冷卻水供回水溫差設(shè)定值、冷卻水回水溫度設(shè)定值、冷凍水供回水溫差。通過PID 控制器控制冷卻水泵、冷凍水泵的運(yùn)行頻率進(jìn)而控制冷卻水供回水溫差、冷凍水供回水溫差，通過PID 控制器控制冷卻塔風(fēng)機(jī)轉(zhuǎn)速控制冷卻水回水溫度，使得系統(tǒng)按照關(guān)聯(lián)準(zhǔn)則表中的最優(yōu)參數(shù)來運(yùn)行。將模型中按照關(guān)聯(lián)準(zhǔn)則運(yùn)行得到的各項數(shù)據(jù)與實際測量（即既有模式的運(yùn)行數(shù)據(jù)）的數(shù)據(jù)進(jìn)行對比，分析關(guān)聯(lián)準(zhǔn)則的有效性和節(jié)能潛力。圖5 為關(guān)聯(lián)準(zhǔn)則的優(yōu)化應(yīng)用流程圖。

圖5 參數(shù)優(yōu)化運(yùn)用流程圖Fig.5 Implementation flowchart of parameter optimization

本研究對空調(diào)制冷季進(jìn)行了優(yōu)化運(yùn)行模擬分析，選擇中等負(fù)荷日的數(shù)據(jù)進(jìn)行展示分析。圖6 與圖7 分別是某中等負(fù)荷日的室外空氣溫度和相對濕度與空調(diào)系統(tǒng)的需求冷量。圖8 是采用關(guān)聯(lián)準(zhǔn)則前后（即既有模式與優(yōu)化模式）的冷凍水供水溫度。優(yōu)化后冷凍水供水溫度大部分運(yùn)行時間內(nèi)高于既有模式運(yùn)行的供水溫度，與關(guān)聯(lián)準(zhǔn)則的目標(biāo)值一致。圖9 是優(yōu)化前后冷凍水供回水溫差對比圖，從圖中可以看出優(yōu)化后冷凍水溫差大于優(yōu)化前的溫差，與關(guān)聯(lián)準(zhǔn)則的目標(biāo)值基本一致。圖10 是冷卻水回水溫度對比圖，從圖中可以看出優(yōu)化模式的回水溫度比既有模式的回水溫度低。圖11 是優(yōu)化前后水系統(tǒng)COP 對比圖，圖中的“目標(biāo)COP”曲線是在該負(fù)荷日運(yùn)行工況下關(guān)聯(lián)準(zhǔn)則中對應(yīng)的COP數(shù)值。從圖中進(jìn)一步可以看出，在一部分的運(yùn)行時間內(nèi)，優(yōu)化后比優(yōu)化前的水系統(tǒng)COP 提升了1.5左右（水系統(tǒng)的能耗包括冷水機(jī)組、冷凍泵、冷卻泵、冷卻塔風(fēng)機(jī)的電耗，不包括末端的能耗），并且在全天運(yùn)行過程中，優(yōu)化模式的空調(diào)水系統(tǒng)COP比關(guān)聯(lián)規(guī)則對應(yīng)的COP 曲線略低，但很接近，說明了關(guān)聯(lián)準(zhǔn)則的有效性。在該負(fù)荷日下，系統(tǒng)既有運(yùn)行模式的水系統(tǒng)能耗為1869.6kWh，采用關(guān)聯(lián)準(zhǔn)則的優(yōu)化運(yùn)行模式的水系統(tǒng)能耗為1705.7kWh，節(jié)能率為8.8%，采用關(guān)聯(lián)準(zhǔn)則進(jìn)行空調(diào)系統(tǒng)的優(yōu)化運(yùn)行有明顯的節(jié)能效果。

圖6 室外空氣溫度和相對濕度Fig.6 Outdoor air temperature and relative humidity

圖7 建筑逐時冷量需求Fig.7 Hourly cooling demand of the building

圖8 既有模式與優(yōu)化模式的冷凍水供水溫度Fig.8 Chilled water supply temperature profiles of the existed operation mode and the optimal operation mode

圖9 既有模式與優(yōu)化模式的冷凍水溫差Fig.9 Chilled water temperature difference profiles of the existed operation mode and the optimal mode

圖10 既有模式與優(yōu)化模式的冷卻水回水溫度Fig.10 Cooling water temperature profiles of the existed operation mode and the optimal mode

圖11 既有模式與優(yōu)化模式的水系統(tǒng)COPFig.11 Water system COP profiles of the existed operation mode and the optimal mode

5 結(jié)論

本文利用關(guān)聯(lián)準(zhǔn)則數(shù)據(jù)挖掘方法，通過對某空調(diào)系統(tǒng)制冷季的運(yùn)行數(shù)據(jù)的挖掘，獲得了在不同工況條件（系統(tǒng)冷量需求和室外空氣溫濕度）下的空調(diào)水系統(tǒng)高效運(yùn)行的參數(shù)，即根據(jù)不同的系統(tǒng)冷量需求、室外空氣溫度和空氣濕度，得到保證系統(tǒng)高效運(yùn)行的冷凍水供水溫度設(shè)定值、冷凍水溫差、冷卻水溫差和冷卻水回水溫度等。進(jìn)一步對關(guān)聯(lián)準(zhǔn)則進(jìn)行了應(yīng)用分析。在中等負(fù)荷日下，優(yōu)化后冷凍水供水溫度大部分運(yùn)行時間內(nèi)高于既有模式運(yùn)行的供水溫度，優(yōu)化后冷凍水溫差大于優(yōu)化前的溫差，優(yōu)化模式的冷卻水回水溫度比既有模式的回水溫度低，與關(guān)聯(lián)準(zhǔn)則的目標(biāo)值基本一致。優(yōu)化后的水系統(tǒng)COP 在大部分運(yùn)行時間內(nèi)比既有模式高1.5左右。在該負(fù)荷日下，既有運(yùn)行模式的能耗為1869.6kWh，采用關(guān)聯(lián)準(zhǔn)則的優(yōu)化模式的能耗為1705.7kWh，節(jié)能率為8.8%。