李傳憲，逯雯雯，石亞男，杜世聰，鄭琬郁，李鵬宇

管道運(yùn)輸具有輸量大、成本低、受環(huán)境影響小的優(yōu)點(diǎn)，在石油天然氣行業(yè)有著越來越重要的作用。管道腐蝕、施工破壞，以及人為破壞等因素使管道頻繁發(fā)生泄漏，造成了嚴(yán)重的人員傷亡、巨大的財(cái)產(chǎn)損失和環(huán)境破壞[1]，因此提高輸油管道泄漏檢測(cè)的準(zhǔn)確性十分必要。目前常用的輸油管道泄漏檢測(cè)方法主要有光纖檢測(cè)法、實(shí)時(shí)模型法及負(fù)壓波法[2]，其中負(fù)壓波法具有精度高、投資少、原理簡(jiǎn)單的優(yōu)點(diǎn)，應(yīng)用最為廣泛。在管道泄漏檢測(cè)中，現(xiàn)場(chǎng)往往存在各種不同工況操作引起的干擾以及傳感器本身的噪聲干擾等，都會(huì)導(dǎo)致負(fù)壓波的信噪比較低，進(jìn)而影響管道泄漏檢測(cè)的準(zhǔn)確性。

針對(duì)管道泄漏引起的負(fù)壓波信號(hào)的去噪研究，小波分析[3]和EMD分解[4]都是常用的方法，然而小波分析不具有自適應(yīng)性，需要選擇合適的小波基函數(shù)和閾值函數(shù)[5]，EMD分解方法雖然具有自適應(yīng)性，但是分解過程會(huì)造成模態(tài)混疊和端點(diǎn)效應(yīng)[6]，針對(duì)這一問題，本文提出了改進(jìn)的添加成對(duì)白噪聲的完全集合經(jīng)驗(yàn)?zāi)B(tài)分解算法（改進(jìn)的CEEMDAN）對(duì)負(fù)壓波信號(hào)進(jìn)行預(yù)處理，該方法不僅具有自適應(yīng)性，克服了模態(tài)混疊，而且消除了輔助白噪聲，提高了運(yùn)算效率，使噪聲信號(hào)和泄漏信號(hào)可以更好的分離開，提高了負(fù)壓波信號(hào)的去噪效果。針對(duì)管道不同工況下信號(hào)的特征提取，能量特征[7?9]和峭度特征[10]常用于正常運(yùn)行、泄漏以及調(diào)泵等差別較大工況的特征描述，而不適用于差別較小工況的特征描述，為了全面表征管道不同工況的壓力信號(hào)特征，本文提出了基于熵的特征向量提取方法，計(jì)算負(fù)壓波信號(hào)的能量熵、峭度熵以及排列熵，從能量、峭度、時(shí)間序列復(fù)雜性三個(gè)角度全面提取不同工況下的信號(hào)特征，并第一次將排列熵應(yīng)用到管道不同工況的識(shí)別中。

1 EMD、EEMD與 CEEMDAN方法介紹及仿真模擬

1.1 EMD及EEMD方法

經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）是希爾伯特黃變換(HHT)方法的關(guān)鍵部分，它依據(jù)信號(hào)自身的時(shí)間尺度特征進(jìn)行分解，而不依賴于任何基函數(shù)，具有自適應(yīng)性[11]。EMD的實(shí)質(zhì)是按照時(shí)間尺度由小到大的順序把不同頻率的分量逐步從原始信號(hào)中分離出來的過程，得到不同頻率的分量稱為固有模態(tài)分量(IMF)[7]。

原始信號(hào)x(t)經(jīng)EMD分解后可表示為:

式中，IMFi(t)為包含不同時(shí)間尺度的固有模態(tài)分量，r(t)為殘余分量，N為固有模態(tài)分量的個(gè)數(shù)。

實(shí)際檢測(cè)到的信號(hào)具有復(fù)雜性，當(dāng)信號(hào)中出現(xiàn)間斷或跳躍性變化時(shí)，傳統(tǒng)EMD分解無法將具有不同特征時(shí)間尺度的成分完全分離開來，會(huì)使不同時(shí)間尺度的信號(hào)出現(xiàn)在同一分量中，產(chǎn)生模態(tài)混疊現(xiàn)象，進(jìn)而導(dǎo)致各IMF失去原來的物理意義。集合經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）利用白噪聲功率譜密度均勻分布的特點(diǎn)，引入輔助白噪聲，在一定程度上削弱了模態(tài)混疊現(xiàn)象[12]。

EEMD的步驟為：

（1）在原始信號(hào)x(t)中添加標(biāo)準(zhǔn)差為kε的零均值白噪聲 n(t)，得到混合信號(hào) x（t）+kεn（t），其中，ε為原始信號(hào)的標(biāo)準(zhǔn)差，k為白噪聲幅度系數(shù)。一般k取0.01～0.50較合適，經(jīng)過大量實(shí)驗(yàn)，本文取其為0.20。

（2）利用EMD對(duì)上述混合信號(hào)進(jìn)行分解。

（3）每次加入不同的n(t)，重復(fù)步驟（1）、（2）。

（4）最后計(jì)算同一特征時(shí)間尺度的所有分量的均值- -- --------IMFi(t)，并將其作為最終的IMFi(t)。

1.2 CEEMDAN及其改進(jìn)算法

EEMD通過多次加入白噪聲求均值的方法，在一定程度上減少了模態(tài)混疊現(xiàn)象的產(chǎn)生，但是加入的白噪聲所帶來的誤差不會(huì)因?yàn)榍缶档倪^程完全消除，并且隨著重復(fù)加入白噪聲次數(shù)的增多，導(dǎo)致運(yùn)算量大大增加。針對(duì)這個(gè)問題，引入添加自適應(yīng)白噪聲的完全集合經(jīng)驗(yàn)?zāi)B(tài)分解（CEEMDAN）[13?14]。

利用EMD進(jìn)行分解得到的第i個(gè)模態(tài)分量用imfi()表示，利用CEEMDAN方法分解所得到的第i個(gè)模態(tài)分量記為- -- -- ------IMFi(t)，那么 CEEMDAN方法步驟如下：

（1）CEEMDAN方法第一個(gè)固有模態(tài)分量的獲得與EEMD算法相同，即對(duì)混合信號(hào)x（t）+kεnj（t），j=1，2，…，M進(jìn)行分解，并保留M次分解的第一個(gè)固有模態(tài)分量，將所得到的分量取平均值，作為CEEMDAN方法的第一個(gè)固有模態(tài)分量。

（2）白噪聲經(jīng)EMD分解得到的模態(tài)分量表示為imfi(kεnj(t))，將其與上一步計(jì)算得到的殘余信號(hào)組合得到混合信號(hào)為 ri(t)+imfi(kεnj(t))，i=1。利用EEMD對(duì)得到的混合信號(hào)進(jìn)行分解，保留每次分解的第一個(gè)固有模態(tài)分量，取平均值可以得到CEEMDAN方法的第二個(gè)固有模態(tài)分量。

（3）繼續(xù)利用步驟（2）的算法對(duì)殘余分量進(jìn)行分解，可得第n+1個(gè)固有模態(tài)分量及其殘余分量為：

（4）當(dāng)殘余分量滿足終止條件時(shí)，算法終止。

在CEEMDAN的算法步驟中可以發(fā)現(xiàn)，多次利用了EEMD方法計(jì)算每階模態(tài)分量。在求得每階固有模態(tài)分量時(shí)仍然需要重復(fù)大量運(yùn)算以消除添加白噪聲的影響，并且添加的白噪聲很難被完全消除。針對(duì)這種現(xiàn)象，提出添加成對(duì)白噪聲的方法，也就是將所加的每個(gè)噪聲都采用正負(fù)成對(duì)的形式添加到信號(hào)中，如式（6）所示：

式中，S是原始信號(hào)；N是添加的白噪聲；M1是原始信號(hào)和正噪聲形成的混合信號(hào)；M2是原始信號(hào)和負(fù)噪聲形成的混合信號(hào)。

通過添加成對(duì)白噪聲進(jìn)行輔助的方法，每添加一次白噪聲所得到的固有模態(tài)分量中，其中一部分包含了正噪聲的殘余量，另一部分包含了負(fù)噪聲的殘余量，在利用均值求最終的固有模態(tài)分量時(shí)，大部分噪聲會(huì)正負(fù)抵消，這樣可以在運(yùn)算量較少的情況下，較精確的去除添加輔助噪聲的影響。

1.3 仿真模擬

以模擬信號(hào)為例對(duì)以上幾種方法進(jìn)行驗(yàn)證，對(duì)余弦信號(hào)x(t)=2cos(20πt)加入間歇性高頻振蕩信號(hào)，采樣點(diǎn)數(shù)為2 000，得到如圖1所示波形。

圖1 合成模擬信號(hào)Fig.1 The synthetic analog signals

分別利用EMD、EEMD、CEEMDAN方法和改進(jìn)的CEEMDAN方法對(duì)上述合成模擬信號(hào)進(jìn)行分解，其中EEMD、CEEMDAN和改進(jìn)的CEEMDAN三種方法都加入幅度系數(shù)為0.20的白噪聲運(yùn)算100次，得到如圖2所示的4種方法分解結(jié)果。從圖2（a）中可以發(fā)現(xiàn)，經(jīng)典EMD方法無法完全把高頻分量和低頻分量分離開，出現(xiàn)了嚴(yán)重的模態(tài)混疊現(xiàn)象。而圖 2（b）、（c）、（d）基本上都把高頻分量和低頻分量分離開了，說明EEMD、CEEMDAN和改進(jìn)的CEEMDAN三種方法基本上消除了模態(tài)混疊現(xiàn)象。分別將EEMD、CEEMDAN和改進(jìn)的CEEMDAN分解后的各個(gè)IMF分量重構(gòu)，并計(jì)算原始信號(hào)能量與重構(gòu)后的信號(hào)能量，Eo=4.040 7×103， EEEMD=4.043 8×103， ECEEMD=4.040 7 × 103，E改CEEMD=4.040 7× 103。發(fā)現(xiàn)EEMD的重構(gòu)信號(hào)引入了很多噪聲，重構(gòu)信號(hào)的能量與原始信號(hào)的能量差別較大，影響后續(xù)處理，而CEEMDAN和改進(jìn)的CEEMDAN幾乎不引入噪聲，重構(gòu)信號(hào)的能量與原始信號(hào)的能量一樣，保證了算法的完備性。

圖2 四種方法對(duì)圖1復(fù)合信號(hào)的分解結(jié)果Fig.2 The decomposing results of Fig.1 composite signal by four methods

對(duì)比圖2（c）、（d）可以發(fā)現(xiàn)，CEEMDAN分解得到的前4個(gè)IMF分量以噪聲為主，無法識(shí)別出其中的高頻分量，而經(jīng)過改進(jìn)的CEEMDAN分解后，前4個(gè)IMF分量包含的噪聲比較少，從第4個(gè)IMF分量開始已經(jīng)明顯包含高頻振蕩分量，并且IMF 5、IMF6分量的光滑性更好，說明加入成對(duì)噪聲的方法，既可以克服模態(tài)混疊，又能消除大量噪聲。對(duì)比CEEMDAN及改進(jìn)的CEEMDAN方法計(jì)算每階IMF分量的迭代次數(shù)，如圖3所示，圖中藍(lán)色盒子里的紅色短線代表迭代次數(shù)的中位數(shù)，而紅色的“+”表示偏離中位數(shù)較大的點(diǎn)，其位置越靠上代表迭代次數(shù)越多，運(yùn)算時(shí)間越長(zhǎng)，可以發(fā)現(xiàn)改進(jìn)的CEEMDAN方法的迭代次數(shù)明顯降低，運(yùn)算效率大幅提高。綜合以上分析，改進(jìn)的CEEMDAN方法，在克服模態(tài)混疊，消除輔助白噪聲和提高運(yùn)算效率三個(gè)方面明顯好于其他三種方法，因此本文選擇改進(jìn)的CEEMDAN方法進(jìn)行后續(xù)處理。

圖3 CEEMDAN方法改進(jìn)前后迭代次數(shù)Fig.3 Iteration times of CEEMDAN method and impr oved CEEM DAN method

2 現(xiàn)場(chǎng)數(shù)據(jù)中的應(yīng)用

實(shí)驗(yàn)以西部某輸油管線為研究對(duì)象，輸油管線長(zhǎng)52.6 km，管徑557 mm，設(shè)計(jì)壓力8 MPa。壓力傳感器安裝在管線兩端。在實(shí)驗(yàn)中，采取人工下油的方式模擬泄漏，將管道首末兩端采集到的負(fù)壓波信號(hào)作為檢測(cè)對(duì)象。為了區(qū)別泄漏與其他正常工況操作引起的負(fù)壓波，分別對(duì)切罐、切泵等工況引起的負(fù)壓波進(jìn)行采集，進(jìn)行分析處理。對(duì)不同工況下的數(shù)據(jù)進(jìn)行統(tǒng)一，下文所有用到的負(fù)壓波信號(hào)數(shù)據(jù)均采用了式（7）進(jìn)行標(biāo)準(zhǔn)化處理：

式中，x(t)代表原始信號(hào)；x′(t)代表標(biāo)準(zhǔn)化以后的信號(hào)；min(x(t))代表原始信號(hào)的最小值；max(x(t))代表原始信號(hào)的最大值。不同工況的時(shí)域波形如圖4所示。

圖4 泄漏、切罐、切泵三種不同工況下的時(shí)域波形Fig.4 Waveforms of leakage,cutting tanks and cutting pumps

2.1 信號(hào)預(yù)處理

除了不同的工況操作會(huì)對(duì)泄漏的識(shí)別產(chǎn)生干擾以外，現(xiàn)場(chǎng)環(huán)境噪聲以及獲取信號(hào)的軟硬件設(shè)備的噪聲都會(huì)導(dǎo)致采集到的負(fù)壓波信噪比較低，影響泄漏工況的識(shí)別精度。因此需要對(duì)采集到的負(fù)壓波信號(hào)進(jìn)行去噪預(yù)處理，以提高工況識(shí)別準(zhǔn)確性。

改進(jìn)的CEEMDAN可以有效抑制模態(tài)混疊現(xiàn)象，并且可以最大程度地降低外加白噪聲對(duì)信號(hào)分解的影響，但是實(shí)際管道泄漏信號(hào)中包含了大量的噪聲，只有提取主要包含泄漏信號(hào)的分量即有效IMF分量，去除噪聲影響，才能最大程度得到有意義的特征值。管道一旦發(fā)生泄漏，其產(chǎn)生的負(fù)壓波會(huì)同時(shí)向管道上下游傳播，本文提出利用管道上下游的兩路信號(hào)進(jìn)行相關(guān)分析的方法來提取有效IMF分量。在管道上下游采集到的負(fù)壓波信號(hào)中，既包含泄漏信號(hào)又包含噪聲信號(hào)，但是只有泄漏信號(hào)具有相關(guān)性，而噪聲信號(hào)互不相關(guān)。利用這一規(guī)律，將一路信號(hào)分解得到的IMF分量分別與另一路信號(hào)進(jìn)行相關(guān)分析就可以提取含有泄漏信號(hào)的有效分量[15]。

以一次實(shí)驗(yàn)獲得的泄漏負(fù)壓波信號(hào)為例進(jìn)行分析，圖5所示分別是管道上下游采集到的壓力信號(hào)x1（n）和x2(n)。圖6所示的IMF分量圖是利用改進(jìn)的CEEMDAN對(duì)上游負(fù)壓波信號(hào)x1(n)進(jìn)行分解得到的，可以發(fā)現(xiàn)前5個(gè)IMF分量以噪聲為主，第6個(gè)IMF分量及以后包含的噪聲明顯減少，但是哪些是有效IMF分量不能進(jìn)行明確判斷，因此需要進(jìn)一步分析。

分別計(jì)算各個(gè)IMF分量與x1(n)和x2(n)的相關(guān)系數(shù)，得到表1，表中r表示原始信號(hào)經(jīng)過分解得到的最終殘余分量。從表1可以發(fā)現(xiàn)，IMF 1和IMF2與x1(n)的相關(guān)系數(shù)明顯偏低，而IMF3-IMF9與x1(n)的相關(guān)系數(shù)差別不大，這可能是因?yàn)镮MF3-IMF9中包含了與x1(n)相關(guān)性較大的噪聲部分，因此無法辨別哪些是相關(guān)分量，而IMF1-IMF5與x2(n)的相關(guān)系數(shù)較其他明顯偏低，并且與圖7所示的信息相吻合，因此前5個(gè)IMF分量以噪聲為主，是應(yīng)當(dāng)去除的不相關(guān)分量，而IMF6-IMF11是應(yīng)當(dāng)與殘余分量一起參與重構(gòu)的有效分量。

圖5 管道上下游采集到的負(fù)壓波信號(hào)Fig.5 Negative pressure wave signals collected fr om upstream and downstream pipes

圖6 上游負(fù)壓波信號(hào)的IMF分量圖Fig.6 IMF components’diagrams of upstream negative pr essur e wave signal

2.2 特征向量提取

準(zhǔn)確的特征向量提取可以大大提高輸油管道泄漏檢測(cè)的準(zhǔn)確性。本文提出了基于熵的特征提取，Shannon提出的信息熵[16]是衡量系統(tǒng)隨機(jī)性的指標(biāo)，信息熵越大代表系統(tǒng)的隨機(jī)性越大。系統(tǒng)X包括多個(gè)事件 X={x1，x2，…，xn}，并且事件發(fā)生的概率為 P={p1，p2，…，pn}。那么事件所攜帶的信息量：I(xi)=-pilg pi，全部事件的信息量累加為系

表1 各個(gè)IMF分量分別與x 1(n)和x2(n)的相關(guān)系數(shù)Table 1 Correlation coefficients between IMF components and x1(n)/x2(n)

當(dāng)管道發(fā)生泄漏時(shí)，其負(fù)壓波信號(hào)所攜帶的能量與其他工況相比有很大的差異性，并且在不同的頻率段中能量分布情況會(huì)發(fā)生變化。因此，本文利用能量熵對(duì)管道泄漏工況進(jìn)行識(shí)別。經(jīng)過2.1信號(hào)預(yù)處理，已經(jīng)獲得了有效IMF分量，然后對(duì)其提取能量熵，其步驟如下：

（1）分別計(jì)算各個(gè)有效IMF分量的能量：

式中，N為信號(hào)的采樣長(zhǎng)度；xi為信號(hào)的幅值。

（2）計(jì)算有效IMF分量的能量總和：

式中，m為有效IMF的個(gè)數(shù)。

（3）用單個(gè)IMF的能量占總能量的比例作為其概率，按照Shannon信息熵的形式計(jì)算每個(gè)樣本的能量熵：

圖7所示為泄漏、切泵、切罐的22個(gè)樣本的能量熵。

圖7 泄漏、切泵和切罐的22個(gè)樣本的能量熵Fig.7 Energy entropy of 22 samples of leakage,cutting pumps and cutting tanks

由圖7可以發(fā)現(xiàn)，切泵產(chǎn)生的負(fù)壓波的能量熵小于泄漏和切罐產(chǎn)生的負(fù)壓波的能量熵，但是泄漏和切罐的能量熵沒有明顯區(qū)別。因此，能量熵可以一定程度上描述三種工況的負(fù)壓波的特征，但是不夠全面。

管道發(fā)生泄漏時(shí)，產(chǎn)生具有顯著沖擊特點(diǎn)的負(fù)壓波。峭度可以描述信號(hào)的沖擊特性，但是峭度值一般只能用來判斷有沒有泄漏發(fā)生，而不能用來區(qū)分其他工況。因此本文仿照能量熵提出了峭度熵的概念，既可以突出峭度對(duì)沖擊信號(hào)靈敏度高的特點(diǎn)，又能考察不同頻率下的沖擊特性。峭度熵的求解過程如下：

（1）分別計(jì)算各個(gè)有效IMF分量的峭度：

式中，xi為信號(hào)的幅值，為信號(hào)的均值，σ為信號(hào)的標(biāo)準(zhǔn)差，N為信號(hào)的采樣長(zhǎng)度。

（2）計(jì)算有效IMF分量的峭度總和：

式中，m為有效IMF的個(gè)數(shù)。

（3）將單個(gè)IMF的峭度占總峭度的比例作為其概率，按照Shannon信息熵的形式計(jì)算每個(gè)樣本的峭度熵：

圖8為泄漏、切泵、切罐的22個(gè)樣本的峭度熵。

圖8 泄漏、切泵和切罐的22個(gè)樣本的峭度熵Fig.8 Kurtosis entropy of 22 samples of leakage,cutting pumps and cutting tanks

由圖8可以發(fā)現(xiàn)，切泵產(chǎn)生的負(fù)壓波的峭度熵最小，泄漏產(chǎn)生的負(fù)壓波的峭度熵最大，而切罐產(chǎn)生的負(fù)壓波的峭度熵處于中間大小，因此峭度熵可以在一定程度上用于三種工況的模式識(shí)別。

排列熵[17?18]（PE）不同于上文所述的能量熵和峭度熵，它是時(shí)間序列的復(fù)雜度的評(píng)價(jià)指標(biāo)，其值越大，時(shí)間序列的復(fù)雜度越高，而且它對(duì)時(shí)間序列的突變現(xiàn)象具有很高的敏感性。輸油管道發(fā)生泄漏或者進(jìn)行其他工況操作時(shí)，壓力波往往會(huì)發(fā)生突變現(xiàn)象，因此可以用排列熵表征不同工況的特征。其求解過程如下：

（1）將2.1中得到的有效IMF分量與殘余分量一起重構(gòu)得到新序列{s(n)，n=1，2，…，N}。

（2）對(duì)新序列進(jìn)行相空間重構(gòu)時(shí)，采用互信息法確定最優(yōu)時(shí)延τ為17，采用偽近鄰法[19]確定最優(yōu)嵌入維數(shù)m為2，得到重構(gòu)序列Xj=[s(j)，s(j+τ)，…，s(j+(m-1)τ)]， 其 中 下 標(biāo) j=1，2，…，N-(m-1)τ。

（3）以重構(gòu)序列Xj為行向量構(gòu)造矩陣X，并對(duì)矩陣X的每一行向量進(jìn)行升序排列，得到s[j+(j1-1)τ]≤ s[j+(j2-1)τ]≤ … ≤ s[j+(jm-1)τ]。其中 j1，j2，…，jm表示重新排列后各個(gè)元素所在列的下標(biāo)，是一組符號(hào)序列。

（4）矩陣X的每一行向量，經(jīng)過重新排列后都可以得到一組符號(hào)序列 S（g）=（j1，j2，…，jm）。其中g(shù)=1，2，…，l，l≤ m!。統(tǒng)計(jì)每一組符號(hào)序列出現(xiàn)的次數(shù)，進(jìn)而可以計(jì)算出每一種符號(hào)序列出現(xiàn)的概

（5）仿照Shannon熵的形式，定義時(shí)間序列{s(n)，n=1，2，…，N}的 PE為：

圖9為泄漏、切泵、切罐的22個(gè)樣本的排列熵，可以發(fā)現(xiàn)，泄漏產(chǎn)生的負(fù)壓波的排列熵明顯小于切泵和切罐產(chǎn)生的負(fù)壓波的排列熵，而切泵和切罐產(chǎn)生的負(fù)壓波的排列熵幾乎重合。這說明，排列熵可以用于泄漏工況的識(shí)別，而對(duì)于其他工況的識(shí)別效果不好。

2.3 工況識(shí)別

2.2 中提出基于熵的特征提取，分別對(duì)能量熵、峭度熵、排列熵進(jìn)行了研究，將能量熵、峭度熵、排列熵作為特征向量，可以全面刻畫不同工況產(chǎn)生的負(fù)壓波的特性，利于后續(xù)的模式識(shí)別，分別以三種工況下的一組信號(hào)為例展示特征向量如表2所示。由于支持向量機(jī)（SVM）適用于小樣本數(shù)據(jù)的分類[20?21]，而現(xiàn)場(chǎng)實(shí)際輸油管道能提供的泄漏樣本很少，因此選用支持向量機(jī)具有較好的分類精度。構(gòu)造兩個(gè)支持向量機(jī)分類模型，分別用于區(qū)分泄漏和切罐，泄漏和切泵，將其命名為SVM 1和SVM 2。SVM 1輸出值為1表示泄漏，-1表示切罐；SVM 2輸出值為1表示泄漏，-1表示切泵。其中每種工況的22個(gè)樣本用于訓(xùn)練模型，另外22個(gè)樣本用于測(cè)試模型的分類效果。

圖9 泄漏、切泵、切罐的22個(gè)樣本的排列熵Fig.9 Permutation entropy of 22 samples for leakage,cutting pumps and cutting tanks

表2 改進(jìn)的CEEMDDAN分解三種工況的特征向量Table 2 The eigenvectors of three working conditions decomposed by improved CEEMDDAN method

利用改進(jìn)的CEEMDAN和EEMD方法分別對(duì)原始負(fù)壓波信號(hào)進(jìn)行分解，然后求其有效IMF分量的能量熵、峭度熵以及排列熵，利用支持向量機(jī)對(duì)三種工況進(jìn)行分類，得到表3。

表3 改進(jìn)的CEEMDAN?SVM（C?S）和EEMD?SVM（E?S）診斷結(jié)果Table 3 Diagnostic results of improved CEEMDAN?SVM(C?S)and EEMD?SVM(E?S)

由表3可知，改進(jìn)的CEEMDAN得到的分類結(jié)果均為100%，而EEMD對(duì)于泄漏和切泵的分類準(zhǔn)確率較高，但是對(duì)泄漏和切罐的分類準(zhǔn)確率較低，說明本文改進(jìn)的CEEMDAN方法效果較好。

為了驗(yàn)證所選特征向量的有效性，在利用改進(jìn)的CEEMDAN分解后，分別單獨(dú)利用能量熵和峭度熵，能量熵和排列熵，峭度熵和排列熵對(duì)三種工況進(jìn)行分類得到表4。由表4可知，能量熵和峭度熵對(duì)泄漏和切泵的識(shí)別效果較好，對(duì)泄漏和切罐的識(shí)別效果較差，說明其對(duì)不同工況的適應(yīng)性較差；能量熵和排列熵組合的工況識(shí)別結(jié)果均較差，峭度熵和排列熵組合的工況識(shí)別效果較好，但仍然不如能量熵、峭度熵和排列熵組合的診斷效果好，說明本文提出的特征向量可以全面描述不同工況的負(fù)壓波特性，具有較強(qiáng)的適應(yīng)性。

3 結(jié) 論

分別對(duì)EMD、EEMD、CEEMDAN和改進(jìn)的CEEMDAN方法進(jìn)行仿真模擬，并對(duì)現(xiàn)場(chǎng)實(shí)際輸油管道不同工況下的負(fù)壓波信號(hào)進(jìn)行預(yù)處理，選用不同的特征向量進(jìn)行分類。得到以下結(jié)論：

（1）通過對(duì)4種方法的仿真模擬研究發(fā)現(xiàn)，改進(jìn)的CEEMDAN方法能有效抑制模態(tài)混疊，消除噪聲，信號(hào)所包含的信息可以比較真實(shí)的通過IMF的特征來反映。將改進(jìn)的CEEMDAN應(yīng)用于輸油管道泄漏檢測(cè)中，發(fā)現(xiàn)基于改進(jìn)的CEEMDAN的管道泄漏識(shí)別方法比基于EEMD的泄漏識(shí)別方法準(zhǔn)確率更高。

（2）提出基于熵的特征向量，通過對(duì)不同特征向量的模式識(shí)別效果進(jìn)行對(duì)比，發(fā)現(xiàn)能量熵、峭度熵、排列熵作為不同工況的特征向量可以全面地反映不同工況的特征信息，比任何單獨(dú)兩種熵組成的特征向量效果更好，對(duì)不同工況具有更強(qiáng)的適應(yīng)性。