孫明遠， 郝學(xué)軍

(北京建筑大學(xué) 環(huán)境與能源工程學(xué)院， 北京 100044)

1 概述

國內(nèi)外學(xué)者對燃氣調(diào)壓器故障診斷的方法有了階段性進展[1-2]，但在燃氣調(diào)壓器故障預(yù)警方面的研究仍相對匱乏。高-中壓燃氣調(diào)壓器運行參數(shù)類型眾多，多變量之間的相關(guān)系數(shù)與典型故障之間存在一定的關(guān)聯(lián)度，準(zhǔn)確找出其關(guān)聯(lián)對實現(xiàn)高-中壓燃氣調(diào)壓器故障診斷與安全預(yù)警至關(guān)重要。因此，本研究采用小波包算法與Apriori算法相結(jié)合的方法實現(xiàn)高-中壓燃氣調(diào)壓器的故障診斷預(yù)警。

2 研究概況

本文的研究對象是高-中壓燃氣調(diào)壓器，型號為FL-BP的間接作用式軸流式調(diào)壓器。

調(diào)壓站調(diào)壓工藝流程見圖1。為了方便開展對調(diào)壓器故障診斷預(yù)警的研究，除去調(diào)壓站內(nèi)已有可獲取的監(jiān)測數(shù)據(jù)類型外，通過外加傳感器自行對監(jiān)測數(shù)據(jù)類型擴展到了8種，監(jiān)測的變量類型見表1。表1中調(diào)壓支路運行調(diào)壓器進口壓力p1、調(diào)壓支路運行調(diào)壓器出口壓力p2和調(diào)壓支路工況流速v這3個量的傳感器監(jiān)測點分別見圖1中點1、點2和點3。

圖1 調(diào)壓工藝流程

表1 監(jiān)測的變量類型

運行調(diào)壓器有4種壓力：運行調(diào)壓器的進口壓力p1；運行調(diào)壓器的出口壓力p2；由進口壓力p1經(jīng)過指揮器過濾、穩(wěn)壓后得到的指揮器穩(wěn)定壓力ps；負載壓力p3，又稱為中間壓力，指由指揮器穩(wěn)定壓力ps和出口壓力p2經(jīng)調(diào)壓彈簧相互作用調(diào)定的一個壓力。

通過外加傳感器可以測得調(diào)壓支路工況流速v和運行調(diào)壓器閥位相對開度K2。

通過大量前期調(diào)研與數(shù)據(jù)分析可知，調(diào)壓器往往在冬供期間(11月15日至次年3月15日)因用氣量大容易在運行過程中出現(xiàn)異常，而非冬供期間天然氣用量小，故障概率隨之減小甚至為0。因此，本文選取某市4個進口壓力為高壓B級、出口壓力為中壓A級的調(diào)壓站，在2017—2018年和2018—2019年兩個冬供期間部分日期的歷史運行數(shù)據(jù)作為原始數(shù)據(jù)，共挑選了147 d，挑選的數(shù)據(jù)包括調(diào)壓器正常運行的數(shù)據(jù)和出現(xiàn)各類故障的數(shù)據(jù)。燃氣調(diào)壓器典型故障分為4類：用氣低谷出口壓力偏高、用氣高峰出口壓力偏低、主閥喘振和指揮器喘振。根據(jù)不同種類故障的嚴(yán)重程度，分為預(yù)警和報警兩種級別。預(yù)警表示調(diào)壓器暫時沒有發(fā)生故障，但未來大概率會發(fā)生故障，應(yīng)當(dāng)予以提示檢修；報警表示調(diào)壓器已發(fā)生故障，應(yīng)立即對其進行檢修。

本研究以1 d的數(shù)據(jù)為1組，1組中有出口壓力p2、負載壓力p3、工況流速v和運行調(diào)壓器閥位相對開度K2這4種變量，每1種變量以1 d 中00: 00—23: 59為周期，每1 min記錄1個數(shù)據(jù)。每組共有5 760個數(shù)據(jù)，總共147組。數(shù)據(jù)統(tǒng)計見表2。

表2 數(shù)據(jù)統(tǒng)計

3 小波包頻帶能量值

3.1 小波包算法理論

在信號處理過程中，小波包算法會對信號進行頻帶分解，將信號中的低頻諧波部分連續(xù)分解為高頻信號和低頻信號，同時還把上層高頻諧波部分再一次分解為高頻信號和低頻信號，從而把頻段分解得更加精細。以3層分解為例，小波包分解結(jié)構(gòu)見圖2。圖2中，S是原始信號序列，L是低頻信號序列，H是高頻信號序列，最末尾數(shù)字是信號分解層數(shù)。

圖2 小波包分解結(jié)構(gòu)

小波包具有能量守恒的物理特點[3]，因此，當(dāng)機械故障發(fā)生時，某些頻帶的能量信號必然會發(fā)生變化，通過分析這些頻帶能量值的變化可以對調(diào)壓器的不同故障類型進行有效識別。在調(diào)壓器實際運行工況下，由于壓力信號是不斷波動變化的，因此，對采集到的調(diào)壓器出口壓力p2以及負載壓力p3這2種時域信號(隨時間變化的動態(tài)信號)，應(yīng)用小波包算法分解為頻帶能量值這種頻域信號(隨頻率變化的動態(tài)信號)，從而實現(xiàn)故障診斷預(yù)警。

根據(jù)圖1可以看出，每個壓力信號經(jīng)3層分解后最終將得到8個頻帶特征能量值，記為Ei(i=1～8)。分解出的E1是低頻能量值，其代表系統(tǒng)的低頻波動強度，波動幅度越大信號就越不穩(wěn)定。因此，E1越大，調(diào)壓器低頻故障(調(diào)壓器出口壓力偏高或者偏低)就越嚴(yán)重。剩余的7個能量值則是高頻能量值，其代表系統(tǒng)的高頻波動強度，頻率越高信號波動就越劇烈。因此，E2～E8越大，調(diào)壓器高頻故障(調(diào)壓器喘振)就越嚴(yán)重。

3.2 頻帶能量值求解

每次在MATLAB中應(yīng)用小波包算法，只調(diào)用1組原始數(shù)據(jù)中的1 440個出口壓力數(shù)據(jù)所組成的序列和1 440個負載壓力數(shù)據(jù)所組成的序列，147組原始數(shù)據(jù)共調(diào)用147次，選擇db6小波函數(shù)[4]對其進行分解，求出壓力信號的小波包頻帶能量值的分布情況，結(jié)果如下。

① 用氣低谷出口壓力偏高故障

預(yù)警的條件：p2的E1在(1.62,1.85]區(qū)間內(nèi)；

報警的條件：p2的E1在(1.85, +∞)區(qū)間內(nèi)。

② 用氣高峰出口壓力偏低故障

預(yù)警的條件：p2的E1在(1.62,1.8]區(qū)間內(nèi)；

報警的條件：p2的E1在(1.8, +∞)區(qū)間內(nèi)。

③ 主閥喘振故障

預(yù)警的條件：p2的E2～E8同時都在(0.52,0.86]區(qū)間內(nèi)，且p3的E2～E8同時都在(0.4,0.95]區(qū)間內(nèi)；

報警的條件：p2的E2～E8同時都在(0.86,+∞)區(qū)間內(nèi)，且p3的E2～E8同時都在(0.4,0.95]區(qū)間內(nèi)。

④ 指揮器喘振故障

預(yù)警的條件：p3的E2～E8同時都在(1.47,1.93]區(qū)間內(nèi)；

報警的條件和預(yù)警的條件完全相同。

由上述結(jié)果可以看出，用氣低谷出口壓力偏高故障和用氣高峰出口壓力偏低故障的頻帶能量值分布基本相同，很難精準(zhǔn)判斷；指揮器喘振故障的預(yù)警與報警判斷中頻帶能量值分布更是完全相同。單靠小波包算法無法實現(xiàn)研究目標(biāo)，因此，需要引入多變量強關(guān)聯(lián)規(guī)則完善模型。

4 高-中壓燃氣調(diào)壓器多變量關(guān)聯(lián)分析

4.1 Apriori算法理論

關(guān)聯(lián)規(guī)則是一種用于探索數(shù)據(jù)趨勢規(guī)律性的方法，Apriori算法[5]是用于尋找關(guān)聯(lián)規(guī)則的算法之一。要歸納關(guān)聯(lián)規(guī)則，首先需要建立數(shù)據(jù)實例項目集，定義項目集如式(1)。

G= {g1,g2, …,gm, …,gn}

(1)

式中G——項目集

gm——項目集G的第m個項目，m=1,2,…,n

n——項目集G中的項目總數(shù)量

G中某些項目所構(gòu)成的集合為事務(wù)，定義事務(wù)如式(2)。

D= {g1,g2, …,gu, …,gw},w≤n,D?G

(2)

式中D——事務(wù)

gu——事務(wù)D的第u個項目，u=1,2,…,w

w——事務(wù)D中的項目總數(shù)量

項集是包含0個或多個項的集合，表示為Wn，即項集Wn為包含了n個項的集合；在事務(wù)集合中對于某個指定的項集被包含的事務(wù)數(shù)量為支持度計數(shù)；項集X和Y中所含的項在事務(wù)集中的頻繁程度(百分比)表示為支持度，項集Y在包含關(guān)聯(lián)規(guī)則項集X的事務(wù)中出現(xiàn)的頻繁程度(百分比)表示為置信度，其中X和Y為項目集G中的某2個項集，且X∩Y=?。 分別記min_sup和min_conf為最小支持度與最小置信度，即用戶所能接受關(guān)聯(lián)結(jié)果的程度；支持度不小于最小支持度閾值的所有項集表示為頻繁項集[6]。

支持度表示規(guī)則的有用性，而置信度表示規(guī)則的確定性。關(guān)聯(lián)規(guī)則挖掘旨在發(fā)現(xiàn)所有高于用戶指定的最小支持度min_sup和最小置信度min_conf的強關(guān)聯(lián)規(guī)則。Apriori算法的主要思想是逐層迭代，通過連接步與剪枝步產(chǎn)生頻繁項集進而產(chǎn)生關(guān)聯(lián)規(guī)則[7]。

4.2 關(guān)聯(lián)規(guī)則求解

由表1變量類型可以看出，在這8種變量中，有3種變量是開關(guān)量，剩下的5種變量為模擬量，由于開關(guān)量用0和1表示關(guān)和開，不滿足關(guān)聯(lián)分析的條件，因此，關(guān)聯(lián)規(guī)則研究采用出口壓力p2、負載壓力p3、支路流速v和閥位相對開度K2這4個模擬量作為研究對象，應(yīng)用Apriori算法來計算多變量間相關(guān)系數(shù)與典型故障的強關(guān)聯(lián)規(guī)則。

4.2.1相關(guān)系數(shù)求解

在求解關(guān)聯(lián)規(guī)則前首先求出原始數(shù)據(jù)的相關(guān)系數(shù)。有學(xué)者研究表明[8]：出口壓力p2與工況流速v、工況流速v與閥位相對開度K2、出口壓力p2與負載壓力p3這3組變量存在一定相關(guān)性，因此，可采用皮爾遜(Pearson)相關(guān)系數(shù)法分別計算出這3組變量的相關(guān)系數(shù)。為方便表達，x和y均代表p2、p3、v和K2這4個變量中的某1個變量，計算公式見式(3)。表2中9種數(shù)據(jù)類型的147組原始數(shù)據(jù)經(jīng)過計算得到441個相關(guān)系數(shù)。

(3)

式中ρ(x,y)——變量x和y的相關(guān)系數(shù)

Cov(x,y)——變量x和y的協(xié)方差

σ(x)——變量x的標(biāo)準(zhǔn)差

σ(y)——變量y的標(biāo)準(zhǔn)差

4.2.2事務(wù)數(shù)據(jù)轉(zhuǎn)換

為了便于Apriori算法程序后續(xù)進行讀取和處理，對相關(guān)系數(shù)和數(shù)據(jù)類型進行符號化處理，進而轉(zhuǎn)換為事務(wù)數(shù)據(jù)：任意變量x和y的相關(guān)系數(shù)ρ(x,y)絕對值的取值范圍為[0,1]，用|ρ(x,y)|來表示，|ρ(x,y)|越大，相關(guān)程度越高，反之，則相關(guān)程度越低。將相關(guān)系數(shù)絕對值以0.05為單位長度劃分成20個單位區(qū)間，并對單位區(qū)間進行符號化處理，相關(guān)系數(shù)符號見表3。例如，數(shù)據(jù)中|ρ(p2,v)|為0.03，則該數(shù)值所對應(yīng)表3中的符號為A1。再把數(shù)據(jù)類型進行符號化處理，見表4。

表3 相關(guān)系數(shù)符號

續(xù)表3

表4 數(shù)據(jù)類型符號

至此，441個相關(guān)系數(shù)及其對應(yīng)的數(shù)據(jù)類型都能依據(jù)表3和表4用符號表示，從而全部轉(zhuǎn)換成441個事務(wù)數(shù)據(jù)。基于147組數(shù)據(jù)，可以得到147條事務(wù)數(shù)據(jù)記錄，挑選部分事務(wù)數(shù)據(jù)記錄見表5。以序號2為例進行說明：當(dāng)|ρ(p2,v)|=0.14，|ρ(v,K2)|=0.65，|ρ(p2,p3)|=0.43時，調(diào)壓器用氣低谷出口壓力偏高故障級別為預(yù)警，轉(zhuǎn)換成事務(wù)數(shù)據(jù)記錄為：A3，B13，C9，F(xiàn)2。

表5 部分事務(wù)數(shù)據(jù)記錄

4.2.3應(yīng)用Apriori算法求得強關(guān)聯(lián)規(guī)則

將147條事務(wù)數(shù)據(jù)記錄輸入到用MATLAB語言編寫的Apriori算法程序中，得到294條調(diào)壓器多變量在不同故障類型下的關(guān)聯(lián)規(guī)則。最小支持度和最小置信度直接影響關(guān)聯(lián)規(guī)則的生成，經(jīng)過反復(fù)實驗對比后，本研究設(shè)定最小支持度為30%，最小置信度為50%。根據(jù)計算結(jié)果，本研究所需的具有實際意義的強關(guān)聯(lián)規(guī)則為滿足最小支持度與最小置信度的關(guān)聯(lián)規(guī)則，經(jīng)篩選后得到的高-中壓調(diào)壓器多變量強關(guān)聯(lián)規(guī)則見表6。

表6 高-中壓調(diào)壓器多變量強關(guān)聯(lián)規(guī)則

對表6中的強關(guān)聯(lián)規(guī)則進行說明：強關(guān)聯(lián)規(guī)則前件表示事件發(fā)生的條件，強關(guān)聯(lián)規(guī)則后件表示事件發(fā)生的結(jié)果。以表6中第①條強關(guān)聯(lián)規(guī)則為例，強關(guān)聯(lián)規(guī)則前件“A1，A2，A3，A4，A5，A6，A7，A8，A9，A10，A11，A12”對照表3，其表達的物理意義為：|ρ(p2,v)|在[0,0.6]范圍內(nèi)；強關(guān)聯(lián)規(guī)則后件“F2”對照表4，表達的物理意義為：用氣低谷出口壓力偏高預(yù)警。因此，第①條強關(guān)聯(lián)規(guī)則為：當(dāng)|ρ(p2,v)|在[0,0.6]范圍內(nèi)時，調(diào)壓器用氣低谷出口壓力偏高故障級別為預(yù)警，故障級別為預(yù)警表示調(diào)壓器極有可能發(fā)生故障。其他強關(guān)聯(lián)規(guī)則表示的含義以此類推。

5 燃氣調(diào)壓器故障診斷與預(yù)警模型

以采用小波包算法對調(diào)壓器出口壓力及負載壓力進行分解得到的頻帶能量值作為基礎(chǔ)判據(jù)，以Apriori算法求解出的強關(guān)聯(lián)規(guī)則作為關(guān)鍵判據(jù)，建立高-中壓燃氣調(diào)壓器故障診斷與預(yù)警模型。

5.1 用氣低谷出口壓力偏高故障模型

用氣低谷出口壓力偏高故障模型見式(4)～(7)。

p2:E1∈(1.62,1.85]

(4)

|ρ(p2,v)|∈[0,0.6]

(5)

p2:E1∈(1.85, +∞)

(6)

|ρ(p2,v)|∈[0,0.5]

(7)

預(yù)警的條件：同時滿足式(4)和式(5)；

報警的條件：同時滿足式(6)和式(7)。

5.2 用氣高峰出口壓力偏低故障模型

用氣高峰出口壓力偏低故障模型見式(5)、(8)～(10)。

p2:E1∈(1.62,1.8]

(8)

p2:E1∈(1.8, +∞)

(9)

|ρ(p2,v)|∈[0,0.45]

(10)

預(yù)警的條件：同時滿足式(5)和式(8)；

報警的條件：同時滿足式(9)和式(10)。

5.3 主閥喘振故障模型

主閥喘振故障模型見式(11) ～ (15)。

p2:Ei (i=2～8)∈(0.52, 0.86]

(11)

p3:Ei (i=2～8)∈(0.4,0.95]

(12)

|ρ(v,K2)|∈[0.56,1]

(13)

p2:Ei (i=2～8)∈(0.86, +∞)

(14)

|ρ(v,K2)|∈[0.66,1]

(15)

預(yù)警的條件：同時滿足式(11) ～ (13)；

報警的條件：同時滿足式(12)、(14)、(15)。

5.4 指揮器喘振故障模型

指揮器喘振故障模型見式(16) ～ (18)。

p3:Ei (i=2～8)∈(1.47,1.93]

(16)

|ρ(p2,p3)|∈[0.51,1]

(17)

|ρ(p2,p3)|∈[0.56,1]

(18)

預(yù)警的條件：同時滿足式(16)和式(17)；

報警的條件：同時滿足式(16)和式(18)。

從上述數(shù)學(xué)模型中可以看出，在小波包頻帶能量值判據(jù)的基礎(chǔ)上，加入多變量相關(guān)系數(shù)與典型故障的強關(guān)聯(lián)規(guī)則后便能有效區(qū)分各類故障以及預(yù)警與報警，彌補了單靠小波包算法無法精確區(qū)分部分類型故障以及無法有效判斷預(yù)警與報警的缺陷，提升了模型的精準(zhǔn)度。

6 模型驗證

首先選取某燃氣公司檢修中心整理的2020年調(diào)壓器檢修記錄中專家技師判定的7組故障歷史數(shù)據(jù)和14組異常歷史數(shù)據(jù)，共21組數(shù)據(jù)。之后把這21組數(shù)據(jù)經(jīng)小波包與Apriori算法分析后代入到高-中壓燃氣調(diào)壓器故障診斷與預(yù)警模型中，該模型對7組故障數(shù)據(jù)進行了報警判斷，報警所判斷的故障與專家技師給出的判定完全吻合。14組異常數(shù)據(jù)模型判斷為預(yù)警，然后經(jīng)檢修中心技師確認，該14組異常數(shù)據(jù)所對應(yīng)的14個調(diào)壓器經(jīng)檢修均發(fā)現(xiàn)不同程度零部件老化破損，如不及時更換零部件，調(diào)壓器將會發(fā)生故障。

因此，此次驗證證明了該模型具備對高-中壓燃氣調(diào)壓器進行故障診斷及預(yù)警的能力。

7 結(jié)論

① 采用小波包算法對調(diào)壓器出口壓力及負載壓力進行分解得到的頻帶能量值可作為基礎(chǔ)判據(jù)。

② Apriori算法得出的多變量相關(guān)系數(shù)與故障類型之間的強關(guān)聯(lián)規(guī)則可以作為關(guān)鍵判據(jù)，提升高-中壓燃氣調(diào)壓器故障診斷與預(yù)警模型的精準(zhǔn)度。

③ 應(yīng)用小波包算法和Apriori算法相結(jié)合的方法建立的高-中壓燃氣調(diào)壓器故障診斷與預(yù)警模型經(jīng)實例驗證可行。