季 云曹 弋

(1.常州信息職業(yè)技術(shù)學(xué)院信息中心，江蘇 常州 213164；2.大連交通大學(xué)交通運(yùn)輸工程學(xué)院，遼寧 大連 116028)

有載調(diào)壓變壓器在電力系統(tǒng)中扮演重要角色，其可以連通電網(wǎng)、調(diào)整電網(wǎng)電壓使其處于穩(wěn)定狀態(tài)，并且可以使電網(wǎng)在較高的電壓中正常運(yùn)行，為電力企業(yè)節(jié)約成本，提高經(jīng)濟(jì)效益[1]。隨著有載調(diào)壓變壓器在電力系統(tǒng)中的應(yīng)用越來(lái)越廣泛，有載調(diào)壓變壓器故障也頻頻出現(xiàn)。有載調(diào)壓變壓器應(yīng)滿足使用需求，對(duì)其實(shí)施有效的檢測(cè)是保障電網(wǎng)正常運(yùn)行的基礎(chǔ)[2]。

相關(guān)領(lǐng)域的大量研究人員對(duì)變壓器故障檢測(cè)進(jìn)行研究，并取得了一定的研究成果。文獻(xiàn)[3]應(yīng)用無(wú)源射頻識(shí)別傳感器標(biāo)簽，對(duì)測(cè)量變壓器振蕩信號(hào)進(jìn)行采集，利用深度學(xué)習(xí)算法，提取測(cè)量信號(hào)特征，通過(guò)加權(quán)貝葉斯分類模型，實(shí)現(xiàn)故障診斷。該算法具有較高的故障診斷正確率，但存在故障檢測(cè)效率低的問(wèn)題。文獻(xiàn)[4]通過(guò)對(duì)故障數(shù)據(jù)進(jìn)行歸一化處理，采用深度信念網(wǎng)絡(luò)，提取故障數(shù)據(jù)特征，運(yùn)用改進(jìn)的模糊C均值聚類，聚類提取的特征，完成對(duì)變壓器故障的診斷。該算法具有能夠有效縮短故障診斷時(shí)間，但存在檢測(cè)結(jié)果誤差較大的問(wèn)題。文獻(xiàn)[5]使用指示量疊加電抗方法來(lái)定位變壓器故障元件。采用用于故障定位的校準(zhǔn)因子，并且可以提供故障電容器元件數(shù)量的實(shí)時(shí)報(bào)告，但算法較復(fù)雜，計(jì)算時(shí)間需要進(jìn)一步提高。基于上述問(wèn)題的存在，本文提出基于奇異值分解的有載調(diào)壓變壓器故障智能檢測(cè)算法。通過(guò)建立Hankel 矩陣，對(duì)設(shè)備信號(hào)進(jìn)行奇異值分解，采用局域波分解方法，分解有載調(diào)壓變壓器設(shè)備中存在的電流反饋信號(hào)，計(jì)算有載調(diào)壓變壓器設(shè)備信號(hào)的瞬時(shí)頻率和幅值，實(shí)現(xiàn)有載調(diào)壓變壓器故障智能檢測(cè)。所提算法的故障檢測(cè)精度和效率較高，檢測(cè)結(jié)果誤差較小，在該領(lǐng)域具有一定意義。

1 電流反饋信號(hào)去噪

因?yàn)橛休d調(diào)壓變壓器設(shè)備故障信號(hào)的獲取往往是直接在設(shè)備周圍的空氣介質(zhì)中獲得，所以信號(hào)的采集實(shí)際上是由多種聲源疊加而成。為保障變壓器故障檢測(cè)的精確度，需要對(duì)電流反饋信號(hào)進(jìn)行去噪處理。有載調(diào)壓變壓器設(shè)備故障智能檢測(cè)方法采用奇異值分解方法，去除有載調(diào)壓變壓器設(shè)備信號(hào)中存在的噪聲。首先建立Hankel 矩陣[6]，對(duì)設(shè)備信號(hào)進(jìn)行奇異值分解，將較小的奇異值設(shè)置為零，得到去噪后有載調(diào)壓變壓器設(shè)備的電流反饋信號(hào)。

假設(shè)x(n)代表有載調(diào)壓變壓器設(shè)備的電流反饋信號(hào)，反饋信號(hào)中包含噪聲信號(hào)u(n)和實(shí)際設(shè)備信號(hào)s(n)。

有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)中存在的噪聲采用奇異值分解(Singular Value Decomposition，SVD)方法去除[7]。重新組織x(n)的結(jié)構(gòu)，得到矩陣A。如果矩陣A是由噪聲和信號(hào)構(gòu)成，則有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)中噪聲和信號(hào)的能力集中，通過(guò)矩陣A的奇異值A(chǔ)1，…，Ai，…，Ar反映[8]。依據(jù)遞減的順序?qū)ζ娈愔颠M(jìn)行排序，去除有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)中存在的噪聲信號(hào)[9]。

設(shè)N代表有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)的長(zhǎng)度，通過(guò)x(n)構(gòu)造Hankel 矩陣X1，則有:

式中:M表示矩陣列數(shù)；K表示矩陣行數(shù)；奇異值分解矩陣X1，將較小奇異值的值設(shè)為零，進(jìn)行奇異值分解反變換，得到有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)q(n)為:

通過(guò)調(diào)節(jié)電流反饋信號(hào)中的噪聲信號(hào)，提高故障檢測(cè)精度。如果有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)中存在的噪聲信號(hào)具有一定周期性[10]，則需分割x(n)，將其分割成m段，每段長(zhǎng)l，得到數(shù)據(jù)矩陣X2為:

根據(jù)分割周期性噪聲的數(shù)據(jù)矩陣X2，得到去噪后的有載調(diào)壓變壓器設(shè)備電流反饋信號(hào)q′(n)為:

以此，完成有載調(diào)壓變壓器電流反饋信號(hào)去噪，為后續(xù)故障智能檢測(cè)提供支持。

2 有載調(diào)壓變壓器故障智能檢測(cè)

2.1 信號(hào)分解

有載調(diào)壓變壓器設(shè)備系統(tǒng)的故障狀態(tài)受外界環(huán)境等因素的干擾，在任意時(shí)刻有載調(diào)壓變壓器設(shè)備均存在若干個(gè)振蕩模式信號(hào)[11]。含有多個(gè)信號(hào)的單信道接收信號(hào)為多分量信號(hào)，當(dāng)有載調(diào)壓變壓器設(shè)備信號(hào)為多分量信號(hào)時(shí)，多個(gè)時(shí)頻重疊信號(hào)會(huì)影響有載調(diào)壓變壓器設(shè)備故障檢測(cè)結(jié)果，則需要分解有載調(diào)壓變壓器設(shè)備信號(hào)，將其分解為若干個(gè)信號(hào)分量[12]，本文通過(guò)局域波分解方法分解設(shè)備信號(hào)。在局域波分解中，所有信號(hào)均由基本模式分量構(gòu)成。

采用局域波分解法求解均值，由局部極值和最小值的包絡(luò)表征實(shí)現(xiàn)。當(dāng)找到有載調(diào)壓變壓器設(shè)備信號(hào)中的全部局部極值點(diǎn)后，用三次樣條連接設(shè)備信號(hào)的局部極值，構(gòu)成上包絡(luò)，通過(guò)同樣的方法得到下包絡(luò)[13]。

如果b1表示上下包絡(luò)線均值，h1表示信號(hào)q′(n)與均值b1之間的差值，那么:

差值h1在理想狀態(tài)下為基本模式分量，在實(shí)際檢測(cè)過(guò)程中，包絡(luò)均值與真實(shí)局部均值之間存在誤差。為了提高有載調(diào)壓變壓器設(shè)備信號(hào)分解的精準(zhǔn)度，需要進(jìn)行過(guò)濾處理，去除包絡(luò)均值的疊加波，提高波形的對(duì)稱性[14]。

式中:分量o1中存在的信號(hào)周期較長(zhǎng)。該過(guò)程可以對(duì)有載調(diào)壓變壓器設(shè)備信號(hào)中的所有分量oj進(jìn)行處理，滿足設(shè)定的條件后停止。最后把有載調(diào)壓變壓器信號(hào)分解成剩余的on和n兩個(gè)局域波分量。

設(shè)X(t)代表分解后的有載調(diào)壓變壓器設(shè)備信號(hào)，其可通過(guò)式(8)進(jìn)行分解，即:

2.2 故障信號(hào)檢測(cè)

根據(jù)信號(hào)分解的結(jié)果.檢測(cè)故障信號(hào)。利用尖峰能量法[15]，對(duì)分解后的有載調(diào)壓變壓器設(shè)備信號(hào)進(jìn)行故障檢測(cè)。計(jì)算有載調(diào)壓變壓器設(shè)備信號(hào)的瞬時(shí)頻率和幅值，構(gòu)建故障檢測(cè)模型，得到有載調(diào)壓變壓器設(shè)備故障的具體類型，進(jìn)行有載調(diào)壓變壓器設(shè)備故障的智能檢測(cè)。

如果τ[X(t)]表示信號(hào)X(t)能量信號(hào)，那么:

式中:β表示振動(dòng)參數(shù)；t表示信號(hào)振動(dòng)時(shí)間；ξ表示信號(hào)振動(dòng)頻率；?代表標(biāo)準(zhǔn)因子。

將式(10)代到式(9)中得到:

設(shè)E代表有載調(diào)壓變壓器設(shè)備信號(hào)的能量常數(shù)，計(jì)算公式如下:

能量常數(shù)與能量信號(hào)之間差值為b/2，因此可以采用能量算子描述有載調(diào)壓變壓器設(shè)備信號(hào)能量的變化情況。

假設(shè)x′(t)代表有載調(diào)壓變壓器設(shè)備信號(hào)的調(diào)頻調(diào)幅信號(hào)，則有:

式中:a(t)代表帶限信號(hào)，根據(jù)公式(13)得到的調(diào)頻調(diào)幅信號(hào)，計(jì)算有載調(diào)壓變壓器設(shè)備信號(hào)的能量信號(hào)τ[x′(t)]，即:

同理，得到下式:

結(jié)合式(15)和式(16)可計(jì)算有載調(diào)壓變壓器設(shè)備信號(hào)的瞬時(shí)頻率η和幅值f，即:

根據(jù)有載調(diào)壓變壓器設(shè)備信號(hào)的瞬時(shí)頻率和幅值構(gòu)建故障檢測(cè)模型V(t)為:

通過(guò)故障檢測(cè)模型可檢測(cè)有載調(diào)壓變壓器設(shè)備的具體故障類型，根據(jù)檢測(cè)結(jié)果，采取對(duì)應(yīng)的措施修理有載調(diào)壓變壓器設(shè)備。

3 實(shí)驗(yàn)分析

為了驗(yàn)證基于奇異值分解的有載調(diào)壓變壓器故障智能檢測(cè)算法的有效性。分別采用所提算法、文獻(xiàn)[3]算法、文獻(xiàn)[4]算法和文獻(xiàn)[5]算法檢測(cè)有載調(diào)壓變壓器設(shè)備故障，對(duì)比4 種不同方法的故障檢測(cè)時(shí)間，標(biāo)準(zhǔn)因子的計(jì)算準(zhǔn)確度以及故障檢測(cè)準(zhǔn)確性。

3.1 實(shí)驗(yàn)環(huán)境

在MATLAB 仿真平臺(tái)上測(cè)試有載調(diào)壓變壓器設(shè)備故障。實(shí)驗(yàn)采用高頻變壓檢測(cè)機(jī)進(jìn)行檢測(cè)，實(shí)驗(yàn)環(huán)境如圖1 所示。

圖1 實(shí)驗(yàn)環(huán)境

3.2 實(shí)驗(yàn)參數(shù)

實(shí)驗(yàn)數(shù)據(jù)由某電科院提供的120 條樣本數(shù)據(jù)，樣本數(shù)據(jù)中包括有載調(diào)壓變壓器正常狀態(tài)(S0)和5 種有載調(diào)壓變壓器故障類型，分別為低溫過(guò)熱(S1)、中溫過(guò)熱(S2)、高溫過(guò)熱(S3)、低能放電(S4)、高能放電(S5)以及局部放電(S6)故障。從樣本數(shù)據(jù)中隨機(jī)抽取100 條為實(shí)驗(yàn)的訓(xùn)練樣本，20 條作為測(cè)試樣本。實(shí)驗(yàn)樣本數(shù)據(jù)故障類型如表1 所示。

表1 實(shí)驗(yàn)樣本數(shù)據(jù)故障類型

3.3 實(shí)驗(yàn)結(jié)果分析

3.3.1 不同方法故障檢測(cè)時(shí)間對(duì)比

為了驗(yàn)證所提算法的有載調(diào)壓變壓器設(shè)備故障檢測(cè)時(shí)間，實(shí)驗(yàn)對(duì)比了4 種方法的檢測(cè)樣本數(shù)據(jù)故障用時(shí)，實(shí)驗(yàn)結(jié)果如圖2 所示。

圖2 不同方法的故障檢測(cè)時(shí)間對(duì)比

分析圖2 可以看出，隨著故障樣本數(shù)據(jù)的增加，不同算法的故障檢測(cè)時(shí)間隨之增大。當(dāng)故障樣本數(shù)據(jù)達(dá)到120 條時(shí)，文獻(xiàn)[3]算法的故障檢測(cè)時(shí)間為58 s，文獻(xiàn)[4]算法的故障檢測(cè)時(shí)間為47 s，文獻(xiàn)[5]算法的故障檢測(cè)時(shí)間為45 s，而所提算法的故障檢測(cè)時(shí)間僅為20 s。由此可知，所提算法的故障檢測(cè)時(shí)間較短。

3.3.2 不同方法的檢測(cè)結(jié)果誤差對(duì)比

為了驗(yàn)證所提算法的可靠性，實(shí)驗(yàn)對(duì)比了4 種方法的誤差平方和，誤差平方和的值越小，表明故障檢測(cè)結(jié)果誤差越小。得到不同方法的檢測(cè)結(jié)果誤差對(duì)比結(jié)果如圖3 所示。

圖3 不同方法的檢測(cè)結(jié)果誤差對(duì)比

分析圖3 可以看出，隨著故障樣本數(shù)據(jù)中訓(xùn)練樣本數(shù)量的增加，不同算法的誤差平方和隨之增大。當(dāng)故障樣本數(shù)據(jù)中訓(xùn)練樣本數(shù)量為100 條，所提算法的誤差平方和為28，文獻(xiàn)[3]算法的誤差平方和為51，文獻(xiàn)[4]算法的誤差平方和為53，文獻(xiàn)[5]算法的誤差平方和為48。由此可知，所提算法的誤差平方和較小，表明所提算法的故障檢測(cè)結(jié)果誤差較小。因?yàn)楸疚脑跈z測(cè)變壓器故障前，對(duì)電流反饋信號(hào)進(jìn)行了去噪處理，利用構(gòu)建的Hankel 矩陣，對(duì)設(shè)備信號(hào)進(jìn)行了奇異值分解，從而得到了去噪后的電流反饋信號(hào)，減小了檢測(cè)誤差。

3.3.3 不同方法的故障檢測(cè)精度對(duì)比

為了進(jìn)一步驗(yàn)證所提算法的有效性，選取20 條測(cè)試樣本，分別采用4 種方法對(duì)實(shí)驗(yàn)中有載調(diào)壓變壓器設(shè)備存在的6 種故障進(jìn)行檢測(cè)，得到不同方法的故障檢測(cè)結(jié)果如表2 所示。

表2 不同方法的故障檢測(cè)精度對(duì)比

分析表2 中的數(shù)據(jù)可知，4 種方法對(duì)每種故障類型檢測(cè)的準(zhǔn)確率有所差異。所提算法檢測(cè)6 種故障的準(zhǔn)確率為95%，而文獻(xiàn)[3]算法、文獻(xiàn)[4]算法及文獻(xiàn)[5]算法檢測(cè)6 種故障的準(zhǔn)確率分別為80%、75%和65%。由此可知，所提算法的故障檢測(cè)精度較高，能夠準(zhǔn)確檢測(cè)有載調(diào)壓變壓器設(shè)備故障類型。因?yàn)楸疚姆椒ㄏ葘?duì)電流反饋信號(hào)進(jìn)行了奇異值分解，去除了信號(hào)噪聲。在多分量信號(hào)分解的基礎(chǔ)上，利用尖峰能量法監(jiān)測(cè)設(shè)備故障信號(hào)，從而保證了故障信號(hào)檢測(cè)的準(zhǔn)確性。

4 結(jié)束語(yǔ)

本文提出的基于奇異值分解的有載調(diào)壓變壓器故障智能檢測(cè)算法，通過(guò)對(duì)有載調(diào)壓變壓器設(shè)備信號(hào)的提取、分解等，將奇異值分解結(jié)合到所提算法中，實(shí)現(xiàn)了有載調(diào)壓變壓器故障智能檢測(cè)。實(shí)驗(yàn)結(jié)果表明，所提算法能夠有效縮短故障檢測(cè)時(shí)間，提高故障檢測(cè)精度，減小檢測(cè)結(jié)果誤差，為有載調(diào)壓變壓器故障的檢測(cè)提供了理論基礎(chǔ)。雖然所提算法在現(xiàn)階段取得了一定成果，但還存在較多不足。有載調(diào)壓變壓器故障類型較多，本文實(shí)驗(yàn)僅考慮了其中的幾種。因此，在今后的工作中將檢測(cè)更多類型的故障，以提高所提算法的普遍適用性。