顧笑偉，王智偉，王占偉，何所畏，閆增峰

?

基于密度權(quán)重支持向量數(shù)據(jù)描述的冷水機(jī)組故障檢測

顧笑偉1，王智偉1，王占偉1，何所畏1，閆增峰2

（1西安建筑科技大學(xué)環(huán)境學(xué)院，陜西西安 710055；2西安建筑科技大學(xué)建筑學(xué)院，陜西西安 710055）

虛警率(FAR)是評價冷水機(jī)組故障檢測性能的關(guān)鍵指標(biāo)，用戶無法接受過高的FAR。為了降低支持向量數(shù)據(jù)描述(SVDD)在冷水機(jī)組故障檢測時的FAR，將密度權(quán)重集成到SVDD中，提出了一種基于密度權(quán)重支持向量數(shù)據(jù)描述(DW-SVDD)的冷水機(jī)組故障檢測方法，該方法考慮了樣本數(shù)據(jù)在真實空間中的密度分布情況。使用ASHRAE RP-1043冷水機(jī)組實驗數(shù)據(jù)對提出的方法進(jìn)行驗證，并將檢測結(jié)果與傳統(tǒng)SVDD的冷水機(jī)組故障檢測方法進(jìn)行比較。結(jié)果表明，提出的方法將FAR從10.5%降低到7%，同比下降超過了30%，同時對4個劣化等級下的7種典型冷水機(jī)組故障有著優(yōu)良的檢測性能。

支持向量數(shù)據(jù)描述；算法；集成；冷水機(jī)組；故障檢測；模型

引 言

冷水機(jī)組是暖通空調(diào)系統(tǒng)中主要耗能部件，冷水機(jī)組運(yùn)行能耗占整個空調(diào)系統(tǒng)能耗的40%～50%[1]。且冷水機(jī)組運(yùn)行一段時間后常常會發(fā)生故障，造成的能源損失占商業(yè)建筑總能耗的15%～30%[2-3]。據(jù)Comstock等[4]統(tǒng)計，42%的維修服務(wù)和26%的維修費(fèi)用是由冷水機(jī)組故障引起的。將故障檢測應(yīng)用于冷水機(jī)組，及時發(fā)現(xiàn)故障并予以排除，對于維持室內(nèi)環(huán)境舒適度，減少設(shè)備損耗，節(jié)約能源具有重要意義[5]。

近些年，基于模式識別的故障檢測方法由于其算法成熟、檢測正確率高、魯棒性好等優(yōu)點[6]，受到了研究人員的廣泛關(guān)注。Han等[7-12]將支持向量機(jī)(SVM)應(yīng)用到冷水機(jī)組故障檢測與診斷中，但該方法在缺少故障數(shù)據(jù)時無法進(jìn)行檢測[13]。由于成本原因，大多數(shù)用戶或廠家不愿進(jìn)行故障模擬實驗，在現(xiàn)場運(yùn)行的冷水機(jī)組中難以獲取大量的故障數(shù)據(jù)，而正常運(yùn)行數(shù)據(jù)容易獲得。因此，冷水機(jī)組故障檢測是典型的一類劃分問題。Zhao等[6,13]將支持向量數(shù)據(jù)描述(SVDD)用于冷水機(jī)組故障檢測與診斷，將故障檢測與診斷看做一類劃分問題。李冠男等提出了一種基于主元分析(PCA)與SVDD的冷水機(jī)組故障檢測方法[14-15]，并將SVDD用于冷水機(jī)組傳感器故障檢測[16-17]。研究結(jié)果表明：將冷水機(jī)組故障檢測看作一類劃分問題，并將SVDD用于冷水機(jī)組故障檢測是有效的。然而SVDD只考慮數(shù)據(jù)點與超球體之間的核距離，未考慮數(shù)據(jù)的密度分布情況，使得到的超球體忽略那些高密度區(qū)域[18-19]。SVDD的這一缺點會導(dǎo)致冷水機(jī)組故障檢測性能下降，主要體現(xiàn)在檢測過程中的虛警率(FAR)較高。FAR是評價故障檢測性能的關(guān)鍵指標(biāo)，對于不會影響生命安全或造成重大經(jīng)濟(jì)損失的HVAC系統(tǒng)，用戶無法接受過高的FAR[2]。

為了降低FAR，本文將密度權(quán)重支持向量數(shù)據(jù)描述(DW-SVDD)集成到冷水機(jī)組故障檢測中。在SVDD中集成密度權(quán)重思想首次被Cha等[20]提出，相比于SVDD，DW-SVDD考慮了數(shù)據(jù)的密度分布情況，使得形成的超球體優(yōu)先描述高密度區(qū)域。而冷水機(jī)組運(yùn)行數(shù)據(jù)存在密度分布不均勻的特性，為此，本文提出了一種基于DW-SVDD的冷水機(jī)組故障檢測方法。采用ASHRAE RP-1043冷水機(jī)組實驗數(shù)據(jù)驗證所提方法的故障檢測性能，并將結(jié)果與基于SVDD的故障檢測方法進(jìn)行比較。

1 SVDD及DW-SVDD算法

1.1 SVDD算法

SVDD是由Tax等[21-23]提出的一種分類算法，其基本思想是將目標(biāo)數(shù)據(jù)映射到高維空間，在高維空間構(gòu)建出一個以為中心為半徑且體積最小的超球體，使其能夠包含絕大多數(shù)目標(biāo)數(shù)據(jù)。假設(shè)有一組目標(biāo)數(shù)據(jù){∈,=1,2,…,}，則SVDD的目標(biāo)函數(shù)為式(1)

其中，為懲罰因子，控制著超球體的體積與分類錯誤率之間的松弛；為松弛變量，它容許一些訓(xùn)練樣本被錯分；為將目標(biāo)數(shù)據(jù)點映射到高維特征空間的非線性映射。

該優(yōu)化問題通過引入拉格朗日乘子求解，同時采用核函數(shù)計算內(nèi)積，得到其對偶形式見式(2)

其中，(,)為滿足mercer條件的核函數(shù)，為拉格朗日乘子。本文采用高斯核函數(shù)，它只有一個可變參數(shù)。根據(jù)的大小，目標(biāo)數(shù)據(jù)分為3類：（1）=0，目標(biāo)數(shù)據(jù)處在超球體內(nèi)部；（2）0<<，目標(biāo)數(shù)據(jù)處在超球面上，稱為無界支持向量；（3）=，目標(biāo)數(shù)據(jù)在超球體外部，它們的松弛變量不為零，稱為有界支持向量。超球體半徑可將任意無界支持向量代入式(3)求出。對于一組觀測數(shù)據(jù)，其到超球體中心的距離見式(4)。

(3)

在進(jìn)行故障檢測時，SVDD反映在二維空間為圖1(a)，根據(jù)圖1(a)，目標(biāo)數(shù)據(jù)為正常運(yùn)行數(shù)據(jù)，SVDD判別邊界即故障檢測邊界，邊界內(nèi)為正常運(yùn)行區(qū)域，邊界外為故障區(qū)域。通過比較觀測數(shù)據(jù)到超球體中心距離()與超球體半徑的大小來判斷是否發(fā)生故障。若()≤，觀測數(shù)據(jù)在超球體內(nèi)，判定為正常；反之，判定為故障。

1.2 DW-SVDD算法

SVDD只考慮數(shù)據(jù)點與超球體之間的距離，未考慮數(shù)據(jù)的密度分布情況，超球體可能會忽略高密度區(qū)域而無法準(zhǔn)確反映目標(biāo)數(shù)據(jù)特征，導(dǎo)致分類性能下降。為了解決這一問題，Cha等[20]在SVDD中引入密度權(quán)重思想，提出了DW-SVDD算法。

密度權(quán)重基于近鄰距離建立，是一種相對密度，它通過比較每個數(shù)據(jù)點的近鄰距離和數(shù)據(jù)集的最大近鄰距離來反映，定義式為

式中，()為的密度權(quán)重，為的近鄰點，(,)為與之間的歐氏距離。位于較高密度區(qū)域的數(shù)據(jù)點距離其臨近數(shù)據(jù)點的更近，因而分配較大的密度權(quán)重。反之，較低密度區(qū)域的數(shù)據(jù)點遠(yuǎn)離其臨近數(shù)據(jù)點，為其分配的密度權(quán)重較小。通過密度權(quán)重()乘以松弛變量，將密度權(quán)重引入SVDD中，目標(biāo)函數(shù)(1)變?yōu)槭?6)

(6)

對式(6)引入拉格朗日乘子求解，同時采用核函數(shù)計算內(nèi)積，優(yōu)化問題變?yōu)槭?7)。求解出式(7)中的拉格朗日乘子后，DW-SVDD的超球體半徑′及觀測數(shù)據(jù)到超球體中心′的距離()′的求解方法同SVDD，見式(3)、式(4)。

DW-SVDD在二維空間進(jìn)行故障檢測時的示意圖見圖1(b)，通過比較觀測數(shù)據(jù)到超球體中心距離′()與超球體半徑′的大小來判斷是否發(fā)生故障，具體判定方法同SVDD。比較圖1(a)、(b)可以看出，由于正常運(yùn)行數(shù)據(jù)密度分布不均勻，在參數(shù)及目標(biāo)數(shù)據(jù)不變的情況下，DW-SVDD的超球體中心′、半徑′相比SVDD均發(fā)生變化。為了使目標(biāo)函數(shù)最小化，高密度區(qū)域具有較高密度權(quán)重的數(shù)據(jù)點被錯分時的懲罰增加，容忍度下降，使超球體向這些區(qū)域轉(zhuǎn)移。相反，低密度區(qū)域數(shù)據(jù)點被錯分時的懲罰減小，容忍度提高，使超球體偏離低密度區(qū)域。

圖1 SVDD及DW-SVDD在二維空間故障檢測示意圖

Fig.1 SVDD and DW-SVDD map in two dimensions space for fault detection

2 基于DW-SVDD的冷水機(jī)組故障檢測流程

圖2為基于DW-SVDD的冷水機(jī)組故障檢測流程圖，該流程分為離線模型訓(xùn)練及在線故障檢測兩部分。離線模型訓(xùn)練階段采用機(jī)組歷史記錄的正常工作運(yùn)行數(shù)據(jù)。數(shù)據(jù)的預(yù)處理，包括穩(wěn)態(tài)檢測、特征選擇及數(shù)據(jù)標(biāo)準(zhǔn)化。其中穩(wěn)態(tài)檢測是為了剔除瞬態(tài)數(shù)據(jù)，篩選出穩(wěn)態(tài)數(shù)據(jù)。特征選擇要結(jié)合現(xiàn)場情況，選出對故障敏感的特征。圖2中，表示第個近鄰，值尋優(yōu)方法將在3.3節(jié)詳細(xì)敘述。根據(jù)確定的值計算-NN距離及每個樣本點的密度權(quán)重，再經(jīng)過網(wǎng)格搜索及五折交叉驗證確定懲罰因子及核參數(shù)(,)[24]，并建立DW-SVDD模型。

在線故障檢測階段對于在線采集到的機(jī)組實時運(yùn)行數(shù)據(jù)進(jìn)行與模型訓(xùn)練相同的數(shù)據(jù)預(yù)處理。基于訓(xùn)練好的DW-SVDD模型，計算觀測數(shù)據(jù)到超球體中心′的距離′并同超球體半徑′進(jìn)行比較。若′>′，觀測數(shù)據(jù)在超球面之外，判定為故障；反之，檢測為正常。

3 實驗驗證與分析

3.1 實驗數(shù)據(jù)

本文數(shù)據(jù)來源于ASHRAE RP-1043冷水機(jī)組故障實驗數(shù)據(jù)[4,25-26]。該實驗采用一臺90冷噸離心式水冷冷水機(jī)組，蒸發(fā)器、冷凝器均為殼管式換熱器，管程為水，制冷劑為R134a。通過特殊設(shè)計的試驗臺模擬了7類典型故障，獲得了27個工況下的64個參數(shù)的數(shù)據(jù)，每類故障由低到高分4個劣化等級，詳細(xì)信息見表1。同時還獲得了機(jī)組正常運(yùn)行時27工況下的正常樣本數(shù)據(jù)。此處，數(shù)據(jù)采集時間間隔為10 s。

表1 冷水機(jī)組故障類型及故障劣化等級

3.2 數(shù)據(jù)預(yù)處理

原始數(shù)據(jù)中摻雜著瞬態(tài)數(shù)據(jù)，由于這些數(shù)據(jù)是機(jī)組在各工況之間動態(tài)運(yùn)行時測得的，對于檢測機(jī)組當(dāng)前的運(yùn)行狀態(tài)會產(chǎn)生干擾，因此必須對原始數(shù)據(jù)進(jìn)行穩(wěn)態(tài)檢測。本文選用Glass等[27-28]提出的方法進(jìn)行穩(wěn)態(tài)檢測。經(jīng)過穩(wěn)態(tài)篩選，從27個工況下的正常樣本中均勻地隨機(jī)選出600個穩(wěn)態(tài)樣本，不同工況樣本混合后，將其隨機(jī)劃分為200個樣本的訓(xùn)練集和400個樣本的測試集；再從每類故障的4個劣化等級中各均勻地隨機(jī)選取27個工況下的400個穩(wěn)態(tài)樣本，混合后，加入測試集。原始實驗數(shù)據(jù)中共有64個參數(shù)變量，其中48個為傳感器直接測得，16個為軟件實時計算得出。然而現(xiàn)場冷水機(jī)組并不具備如此多的傳感器，且有些傳感器對故障反應(yīng)不敏感[29]。本文主要從以下兩個方面進(jìn)行特征選擇：（1）特征參數(shù)在現(xiàn)場易獲得，且數(shù)量不宜過多；（2）特征參數(shù)對故障反應(yīng)敏感，以提高檢測正確率。Zhao等[30]對14個現(xiàn)場冷水機(jī)組的在線測量系統(tǒng)進(jìn)行了調(diào)研，由此確定了現(xiàn)場容易獲得的測量參數(shù)。本文根據(jù)這一結(jié)果，選擇了易在現(xiàn)場獲得的16個特征，且Zhao等[4]已經(jīng)驗證了這16個特征對于冷水機(jī)組故障的敏感性，見表2。

表2 冷水機(jī)組故障檢測的特征選擇

3.3 參數(shù)尋優(yōu)

在冷水機(jī)組故障檢測過程（圖2）中，共有3個參數(shù)需要優(yōu)化，近鄰中的值，懲罰因子，以及核參數(shù)。首先是要確定值，目前未有統(tǒng)一的值確定方法。本文中令分別等于1到20，迭代計算出訓(xùn)練集各數(shù)據(jù)點的密度權(quán)重，并分別訓(xùn)練DW-SVDD模型，從中選出訓(xùn)練集FAR最低時所對應(yīng)的值。本文經(jīng)過尋優(yōu)得到=9。

對于(,)，網(wǎng)格搜索范圍初步定為[2-4,24]和[2-4,24]，圖3為使用五折交叉驗證進(jìn)行網(wǎng)格搜索的結(jié)果。從圖中可以看出，五折交叉驗證正確率大于90%的區(qū)域并不是很大，且當(dāng)>20時，交叉驗證正確率基本穩(wěn)定不變，而隨著的增加，交叉驗證正確率提升很明顯。懲罰因子控制著超球體體積與分類錯誤率之間的松弛，值越小，超球體體積越小，在保證正確率的前提下，盡可能取小值?？刂瞥蝮w形狀，值越大，決策邊界越松弛同樣盡量取小值。本文中，訓(xùn)練集的交叉驗證正確率控制在90%左右。進(jìn)一步縮小尋優(yōu)范圍，在[2-1.5,2-0.5]和[21.5,21.75]范圍內(nèi)尋優(yōu)，找出訓(xùn)練集FAR最低時(、)的值，并訓(xùn)練DW-SVDD模型。

3.4 故障檢測結(jié)果與分析

3.4.1 冷水機(jī)組運(yùn)行數(shù)據(jù)密度分布 根據(jù)3.3節(jié)尋優(yōu)確定的=9計算訓(xùn)練集各樣本點的密度權(quán)重見圖4，由1.2節(jié)可知，密度權(quán)重是一種相對密度，能夠反映數(shù)據(jù)的密度分布情況。高密度區(qū)域數(shù)據(jù)分布密集具有較高的密度權(quán)重，低密度區(qū)域數(shù)據(jù)分布稀疏具有較低的密度權(quán)重。由圖4可以看出各樣本的密度權(quán)重均不相同，最大權(quán)重接近1，最小為零，且整體出現(xiàn)了明顯的波動，這表明冷水機(jī)組運(yùn)行數(shù)據(jù)密度分布不均勻。如采用SVDD，各數(shù)據(jù)點等權(quán)對待，無法反映數(shù)據(jù)的密度分布情況。而采用DW-SVDD時，根據(jù)數(shù)據(jù)密度分布情況為各數(shù)據(jù)點分配密度權(quán)重，使最終得到的超球體向高密度區(qū)域偏移的同時，偏離低密度區(qū)域，從而提升冷水機(jī)組故障檢測性能。

3.4.2 對正常樣本檢測結(jié)果 圖5為采用SVDD及DW-SVDD進(jìn)行故障檢測時訓(xùn)練集和測試集的FAR?？梢钥闯?，DW-SVDD訓(xùn)練集、測試集的FAR均低于SVDD。其中，訓(xùn)練集的FAR降低最多，相比SVDD從6%變?yōu)?%；測試集的FAR相比SVDD從10.5%變?yōu)?%，同比下降超過30%?？梢?，DW- SVDD的故障檢測FAR相比SVDD有明顯改善。

采用SVDD和DW-SVDD對正常樣本測試集進(jìn)行故障檢測時各樣本的結(jié)果分別見圖6(a)和(b)。圖中橫坐標(biāo)為測試集樣本編號，縱坐標(biāo)為測試樣本到超球體中心距離的平方，虛線表示超球體半徑的平方值，虛線上方視為故障，下方視為正常。對比圖6(a)和(b)發(fā)現(xiàn)，SVDD進(jìn)行故障檢測時，有很大一部分錯分?jǐn)?shù)據(jù)點到中心的距離與超球體半徑大小相當(dāng)?？梢?，這些錯分?jǐn)?shù)據(jù)點距離超球體很近，且多分布在高密度區(qū)域，一旦被錯分會對分類性能造成較大影響，導(dǎo)致FAR明顯偏大。而DW-SVDD進(jìn)行故障檢測時，絕大多數(shù)錯分?jǐn)?shù)據(jù)點到超球體中心的距離明顯大于其半徑。說明這些錯分?jǐn)?shù)據(jù)點遠(yuǎn)離超球體，且處在數(shù)據(jù)密度分布較低的區(qū)域，即便被錯分也不會對分類性能造成太大影響。同時，由圖6(a)和(b)可以發(fā)現(xiàn)，SVDD和DW-SVDD的超球體半徑及各樣本點到中心的距離并未有明顯差異，這說明DW-SVDD的超球體中心′、半徑′相比SVDD僅發(fā)生微小變化，卻使超球體向高密度區(qū)域偏移的同時，偏離低密度區(qū)域，使得錯分?jǐn)?shù)據(jù)減少，F(xiàn)AR顯著降低。

3.4.3 對故障樣本檢測結(jié)果 表3為SVDD、DW-SVDD對于7類故障在4劣化等級下的故障檢測正確率?？梢钥闯?，對任一故障的檢測正確率隨著劣化等級的增加而提高。表3顯示，與SVDD相比，對大多數(shù)故障，DW-SVDD保持了同樣高的檢測正確率，甚至提高了某些故障的檢測正確率，如RedEW。同時顯示，DW-SVDD對于CdFoul、Ncg、ExOil這3類故障檢測正確率最高，幾乎都能達(dá)到100%。RefOver、RefLeak、RedCW在劣化等級為SL-1時的故障檢測正確率分別為82.5%、60.8%、57.3%，在劣化等級為SL-2時分別為88.8%、55.8%、79.3%，但在劣化等級超過SL-3時，正確率也都接近100%。對RedEW的檢測正確率較低，但在各劣化等級下，正確率均超過50%，這對于報警一個故障依然是足夠的。

表3 兩種不同方法對故障數(shù)據(jù)的故障檢測結(jié)果

DW-SVDD在目標(biāo)函數(shù)式(7)中引入了正常運(yùn)行數(shù)據(jù)的密度權(quán)重()，且正常運(yùn)行數(shù)據(jù)密度分布不均勻，使得訓(xùn)練好的DW-SVDD模型能夠反映正常運(yùn)行數(shù)據(jù)密度分布情況，但并未反映故障運(yùn)行數(shù)據(jù)的密度分布情況，而通常故障運(yùn)行數(shù)據(jù)分布有別于正常運(yùn)行數(shù)據(jù)。以劣化等級為SL-1時的ExOil為例，采用SVDD和DW-SVDD對其進(jìn)行故障檢測時各樣本的結(jié)果分別見圖7(a)和(b)。由圖7不難看出，DW-SVDD對故障檢測時，與SVDD差異不大。這是因為DW-SVDD的超球體中心′、半徑′相比SVDD僅發(fā)生微小變化。除極個別的樣本外，絕大多數(shù)故障樣本點到中心的距離同超球體半徑有著明顯的差異?？梢酝茰y，超球面附近故障樣本的分布極少。因此，與SVDD相比，對大多數(shù)故障，DW-SVDD保持了同樣高的檢測正確率。

需要指出的是，本文在進(jìn)行模型訓(xùn)練、構(gòu)建超球體的過程中，將27個工況下正常數(shù)據(jù)混合進(jìn)行模型訓(xùn)練，并未對27個工況的正常數(shù)據(jù)分別進(jìn)行模型訓(xùn)練。雖然不同工況下冷水機(jī)組運(yùn)行數(shù)據(jù)分布存在一定差異，針對不同工況構(gòu)建不同超球體，或許可以更好地反映各個工況的數(shù)據(jù)分布，但檢測計算量、計算時間等都會有較大幅度增加。事實上，本文中的超球體是在高維空間進(jìn)行構(gòu)建，且考慮了數(shù)據(jù)的密度分布情況，能夠較好地包含各工況下的正常數(shù)據(jù)，從而將各工況下正常數(shù)據(jù)與故障數(shù)據(jù)區(qū)分開，并取得較滿意的結(jié)果。

4 結(jié) 論

為了克服SVDD沒有考慮數(shù)據(jù)密度分布的缺點，降低故障檢測FAR，本文在SVDD中集成密度權(quán)重的思想，提出了一種基于DW-SVDD的冷水機(jī)組故障檢測方法，并采用RP-1043實驗數(shù)據(jù)驗證了該方法的有效性，發(fā)現(xiàn)冷水機(jī)組運(yùn)行數(shù)據(jù)存在密度分布不均勻的情況，隨后將該方法故障檢測性能同SVDD故障檢測方法進(jìn)行了比較，得出如下結(jié)論。

（1）同SVDD故障檢測方法相比，本文提出的DW-SVDD冷水機(jī)組故障檢測方法顯著降低了FAR，對測試集的FAR由10.5%下降到7%，同比下降了30%以上；

（2）與SVDD相比，對大多數(shù)故障，DW-SVDD保持了同樣高的檢測正確率，甚至提高了某些故障的檢測正確率，如RedEW。對于CdFoul、Ncg、ExOil的檢測正確率幾乎都能達(dá)到100%；RefOver、RefLeak、RedCW在劣化等級超過SL-3時，正確率也都接近100%；RedEW較難檢測，但各劣化等級的檢測正確率均高于50%?？梢姡珼W-SVDD對4個劣化等級下的7種典型冷水機(jī)組故障有著優(yōu)良的檢測性能。

符 號 說 明

a——SVDD超球體中心 a′——DW-SVDD超球體中心 C——懲罰因子 d(xi,xki)——xi與xki之間的歐氏距離 D(z)——觀測數(shù)據(jù)z到SVDD超球體中心的歐氏距離 D′(z)——觀測數(shù)據(jù)z到DW-SVDD超球體中心的歐氏距離 g——高斯核參數(shù)， k——k近鄰中的k值 K(xi,xj)——核函數(shù) max——極大值優(yōu)化問題 min——極小值優(yōu)化問題 R——SVDD超球體半徑 Rd——d維特征空間 R′——DW-SVDD超球體半徑 s.t.——約束條件 xi——目標(biāo)數(shù)據(jù) xki——xki為xi的k近鄰點 z——觀測數(shù)據(jù) ai——拉格朗日乘子 xi——松弛變量 r(xi)——xi的密度權(quán)重 s——高斯核參數(shù) j——非線性映射 下角標(biāo) i, j——數(shù)據(jù)集矩陣的行序號 train set——訓(xùn)練集

References

[1] DOE. Buildings Energy Data Book[R]. Washington DC: U.S. Department of Energy, 2007.

[2] KATIPAMULA S, BRAMBLEY M R. Review article: methods for fault detection, diagnostics, and prognostics for building systems—a review(Ⅰ)[J]. HVAC&R Research, 2005, 11(1): 3-25.

[3] KATIPAMULA S, BRAMBLEY M R. Review article: methods for fault detection, diagnostics, and prognostics for building systems—a review(Ⅱ)[J]. HVAC&R Research, 2005, 11(2): 169-187.

[4] COMSTOCK M C, BRAUN J E, BERNHARD R. Development of Analysis Tools for the Evaluation of Fault Detection and Diagnostics in Chillers[M]. Purdue University, 1999: 25.

[5] LI H R, BRAUN J E. Economic evaluation of benefits associated with automated fault detection and diagnosis in rooftop air conditioners[J]. ASHRAE Transactions, 2007, 113(2): 200-210.

[6] ZHAO Y, WANG S W, XIAO F. Pattern recognition-based chillers fault detection method using support vector data description (SVDD)[J]. Applied Energy, 2013, 112: 1041-1048.

[7] HAN H, CAO Z K, GU B,. PCA-SVM-based automated fault detection and diagnosis (AFDD) for vapor-compression refrigeration systems[J]. HVAC&R Research, 2010, 16(3): 295-313.

[8] HAN H, GU B, KANG J,. Study on a hybrid SVM model for chiller FDD applications[J]. Applied Thermal Engineering, 2011, 31(4): 582-592.

[9] HAN H, GU B, HONG Y C,. Automated FDD of multiple-simultaneous faults (MSF) and the application to building chillers[J]. Energy and Buildings, 2011, 43(9): 2524-2532.

[10] 韓華, 谷波, 康嘉. 基于遺傳算法和支持矢量機(jī)參數(shù)優(yōu)化的制冷機(jī)組故障檢測與診斷研究[J]. 機(jī)械工程學(xué)報, 2011, 47(16): 120-126. HAN H, GU B, KANG J. Chiller fault detection and diagnosis based on genetic algorithm and support vector machine with parameter tuning[J]. Journal of Mechanical Engineering, 2011, 47(16):120-126.

[11] 谷波, 韓華, 洪迎春, 等. 基于 SVM 的制冷系統(tǒng)多故障并發(fā)檢測與診斷[J]. 化工學(xué)報, 2012, 62(S2): 112-119. GU B, HAN H, HONG Y C,. SVM-based FDD of multiple-simultaneous faults for chillers[J]. CIESC Journal, 2012, 62(S2): 112-119.

[12] 韓華, 谷波, 任能. 主元分析和支持向量機(jī)在制冷系統(tǒng)故障診斷中的應(yīng)用[J]. 上海交通大學(xué)學(xué)報, 2011, 45(9): 1355-1361. HAN H, GU B, REN N. Fault diagnosis for refrigeration systems based on principal component analysis and support vector machine[J]. Journal of Shanghai Jiaotong University, 2011, 45(9): 1355-1361.

[13] ZHAO Y, XIAO F, WEN J,. A robust pattern recognition-based fault detection and diagnosis (FDD) method for chillers[J]. HVAC&R Research, 2014, 20(7): 798-809.

[14] LI G N, HU Y P, CHEN H X,. An improved fault detection method for incipient centrifugal chiller faults using the PCA-R-SVDD algorithm[J]. Energy and Buildings, 2016, 116: 104-113.

[15] 李冠男, 胡云鵬, 陳煥新, 等. 基于 PCA-SVDD 的冷水機(jī)組故障檢測方法[J]. 華中科技大學(xué)學(xué)報(自然科學(xué)版), 2015, 43(8): 119-122. LI G N, HU Y P, CHEN H X,. PCA-SVDD-based chiller fault detection method[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2015, 43(8):119-122.

[16] 李冠男, 胡云鵬, 陳煥新, 等. 基于 SVDD 的冷水機(jī)組傳感器故障檢測及效率分析[J]. 化工學(xué)報, 2015, 66(5): 1815-1820. LI G N, HU Y P, CHEN H X,. SVDD-based chiller sensor fault detection method and its detection efficiency[J]. CIESC Journal, 2015, 66(5): 1815-1820.

[17] LI G N, HU Y P, CHEN H X,. A sensor fault detection and diagnosis strategy for screw chiller system using support vector data description-based D-statistic and DV-contribution plots[J]. Energy and Buildings, 2016, 133: 230-245.

[18] LEE K Y, KIM D W, LEE D,. Improving support vector data description using local density degree[J]. Pattern Recognition, 2005, 38(10): 1768-1771.

[19] LEE K Y, KIM D W, LEE K H,. Density-induced support vector data description[J]. IEEE Transactions on Neural Networks, 2007, 18(1): 284-289.

[20] CHA M, KIM J S, BAEK J G. Density weighted support vector data description[J]. Expert Systems with Applications, 2014, 41(7): 3343-3350.

[21] TAX D M J, DUIN R P W. Support vector data description[J]. Machine Learning, 2004, 54(1): 45-66.

[22] TAX D M J, DUIN R P W. Support vector domain description[J]. Pattern Recognition Letters, 1999, 20(11): 1191-1199.

[23] TAX D M J, YPMA A, DUIN R P W. Pump failure detection using support vector data descriptions[C]//HAND D J, KOK J N, BERTHOLD M R. International Symposium on Intelligent Data Analysis. Berlin Germany: Springer-Verlag, 1999: 415-425.

[24] HSU C W, CHANG C C, LIN C J. A practical guide to support vector classification[EB/OL].[2016-05-19]. http://www.csie.ntu.edu.tw/～cjlin/papers/guide/guide.pdf.

[25] COMSTOCK M C, BRAUN J E, GROLL E A. A survey of common faults for chillers/discussion[J]. Ashrae Transactions, 2002, 108: 819.

[26] COMSTOCK M C, BRAUN J E, BERNHARD R. Experimental Data from Fault Detection and Diagnostic Studies on a Centrifugal Chiller[M]. Purdue University, 1999:16-66.

[27] GLASS A S, GRUBER P, ROOS M,. Qualitative model-based fault detection in air-handling units[J]. IEEE Control Systems, 1995, 15(4): 11-22.

[28] HAN H, GU B, WANG T,. Important sensors for chiller fault detection and diagnosis (FDD) from the perspective of feature selection and machine learning[J]. International Journal of Refrigeration, 2011, 34(2): 586-599.

[29] DASH M, LIU H. Feature selection for classification[J]. Intelligent Data Analysis, 1997, 1(3): 131-156.

[30] ZHAO X Z, YANG M, LI H R. Field implementation and evaluation of a decoupling-based fault detection and diagnostic method for chillers[J]. Energy and Buildings, 2014, 72: 419-430.

Chiller fault detection by density weighted support vector data description

GU Xiaowei1, WANG Zhiwei1, WANG Zhanwei1, HE Suowei1, YAN Zengfeng2

(1School of Environmental, Xi’an University of Architecture and Technology, Xi’an 710055, Shaanxi, China;2School of Architecture, Xi’an University of Architecture and Technology, Xi’an 710055, Shaanxi, China)

False alarm rate (FAR) is a key indicator to evaluate performance of chiller fault detection methods, since customers cannot accept high FAR. In order to reduce FAR of support vector data description (SVDD)-based chiller fault detection, a density weighted support vector data description (DW-SVDD)-based chiller fault detection method was proposed by integration of density weight into SVDD with a consideration of density distribution of sample data in real space. The proposed method was validated with experimental data of RP-1043 ASHRAE and detection results were compared to those of traditional SVDD chiller fault detection methods. The results showed that the new method could reduce FAR from 10.5% to 7%, which was lowered about 30%, and had excellent detection performance for 7 typical chiller faults at 4 severity levels.

support vector data description; algorithm; integration; chiller; fault detection; model

10.11949/j.issn.0438-1157.20161077

TB 65

A

0438—1157（2017）03—1099—10

“十二五”國家科技支撐計劃項目（2011BAJ03B06）。

2016-08-01收到初稿，2016-11-04收到修改稿。

聯(lián)系人：王智偉。第一作者：顧笑偉（1991—），男，碩士研究生。

2016-08-01.

Prof.WANG Zhiwei, wzhiwei-atu@163.com

supported by “Twelfth Five-Year” National Key Technology Research and Development Program of China(2011BAJ03B06).