陳朝陽(yáng)，楊文榮，張雨蒙，石小暉

（河北工業(yè)大學(xué)省部共建電工裝備可靠性與智能化國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300130）

0 引 言

為保障國(guó)民經(jīng)濟(jì)提供可靠的能源供應(yīng)，亟需制定動(dòng)態(tài)檢修策略，進(jìn)而量化故障帶來(lái)的風(fēng)險(xiǎn)。變壓器機(jī)械故障屬于潛在性故障，初期大多是由機(jī)械結(jié)構(gòu)問題所引起，其中繞組故障所占比重較大［1］。配電變壓器繞組從外觀結(jié)構(gòu)上有圓形、長(zhǎng)圓形、橢圓形等形狀，長(zhǎng)圓形繞組設(shè)計(jì)具有優(yōu)良的節(jié)材效果，廣泛應(yīng)用于中小型配電變壓器設(shè)計(jì)中［2］。長(zhǎng)圓形繞組的配電變壓器向大容量方向發(fā)展，在配電網(wǎng)中大力推廣，長(zhǎng)圓形繞組因機(jī)械性能較差，外部短路事故時(shí)有發(fā)生，影響其使用壽命［3］，因此，對(duì)長(zhǎng)圓形繞組進(jìn)行早期機(jī)械故障診斷具有重要意義。

目前，針對(duì)長(zhǎng)圓形繞組的機(jī)械故障診斷研究較少，主要研究圓形繞組松動(dòng)、形變等故障最多［4，5］。針對(duì)圓形繞組機(jī)械故障診斷，Wu Y Y等人［6］采用麻雀搜索算法（sparrow search algorithm，SSA）與深度信念網(wǎng)絡(luò)（deep belief network，DBN）相結(jié)合的方法對(duì)振動(dòng)信號(hào)的特征向量進(jìn)行分類，可以有效提高分類算法的分類精度5 %～20 %。Zhou L J 等人［7］提出基于二叉樹支持向量機(jī)（binary tree support vector machine，BT-SVM）的診斷模型，可以解決小樣本學(xué)習(xí)的問題，可用于電力變壓器的機(jī)械故障診斷。但上述方法存在特征提取量較大、診斷過程較為復(fù)雜的問題，不利于實(shí)時(shí)在線故障診斷。因此，設(shè)計(jì)一種算法流程簡(jiǎn)單、樣本點(diǎn)少、診斷速度快的故障診斷方法是有必要的。

針對(duì)長(zhǎng)圓形繞組更容易產(chǎn)生位移形變的問題，本文首先建立了長(zhǎng)圓形繞組輻向和軸向數(shù)學(xué)模型，分析其不同方向上的振動(dòng)特性；其次，建立長(zhǎng)圓形繞組的有限元仿真模型，分析其不同位置的振動(dòng)特性；最后，設(shè)計(jì)繞組松動(dòng)和翹曲故障實(shí)驗(yàn)，針對(duì)繞組薄弱位置進(jìn)行小波包能量特征提取，采用灰狼優(yōu)化（grey wolf optimization，GWO）算法優(yōu)化SVM參數(shù)的方法對(duì)長(zhǎng)圓形繞組不同機(jī)械狀態(tài)進(jìn)行故障診斷。

1 變壓器長(zhǎng)圓形繞組振動(dòng)數(shù)學(xué)模型分析

長(zhǎng)圓形高壓和低壓繞組輻向電磁力如圖1所示。

圖1 輻向電磁力示意

由圖1可發(fā)現(xiàn)，長(zhǎng)圓形繞組的輻向電磁力Fx比輻向電磁力Fy強(qiáng)，F(xiàn)x易使高、低壓繞組導(dǎo)線被拉長(zhǎng)，進(jìn)而發(fā)生永久形變；此外，匝間絕緣也被拉長(zhǎng)，導(dǎo)致匝間絕緣破裂，進(jìn)而形成匝間短路。

對(duì)于工頻電流為i＝Imcosωt的長(zhǎng)圓形繞組，將其振動(dòng)過程分解為輻向和軸向

式中Fr和Fz分別為輻向和軸向電磁力，kz和kr分別為繞組軸向和輻向坐標(biāo)的相關(guān)系數(shù)矩陣，Im為繞組電流的有效值。

1.1 繞組輻向振動(dòng)數(shù)學(xué)模型

建立繞組輻向振動(dòng)的動(dòng)力學(xué)方程為

式中M為質(zhì)量矩陣，Cr為輻向阻尼矩陣，μ為摩擦系數(shù)，zr為各節(jié)點(diǎn)輻向位移列向量。

式（2）的齊次通解為

式中D1，D2為常數(shù)，與繞組振動(dòng)初始條件有關(guān)。

式（2）有兩項(xiàng)特解，第一項(xiàng)特解為

第二項(xiàng)特解為

繞組輻向振動(dòng)的動(dòng)力學(xué)方程通解為

在負(fù)載電流不變的條件下，繞組輻向振動(dòng)位移與繞組質(zhì)量、輻向阻尼、絕緣樹脂與墊塊的摩擦系數(shù)、繞組預(yù)緊力及繞組位置有關(guān)，主要振動(dòng)頻率為基頻（100 Hz）及其倍頻。

1.2 繞組軸向振動(dòng)數(shù)學(xué)模型

建立繞組軸向振動(dòng)的動(dòng)力學(xué)方程為

式中 zz為各節(jié)點(diǎn)軸向位移列向量。Cz為軸向阻尼矩陣，軸向剛度矩陣K。

式（7）的齊次通解為

式中A1和θ取決于方程初始條件的積分常數(shù)。

式（7）有兩項(xiàng)特解，第一項(xiàng)特解為

第二項(xiàng)特解為

繞組軸向振動(dòng)的動(dòng)力學(xué)方程通解為

在負(fù)載電流不變的條件下，繞組軸向振動(dòng)位移與繞組質(zhì)量、剛度、軸向阻尼、預(yù)緊力及繞組位置有關(guān)，主要振動(dòng)頻率為基頻100 Hz。

當(dāng)繞組發(fā)生結(jié)構(gòu)變形時(shí)，主要改變100 Hz倍頻的振動(dòng)；當(dāng)繞組發(fā)生松動(dòng)時(shí)，主要影響基頻100 Hz的振動(dòng)，即機(jī)械狀態(tài)改變引起振動(dòng)特征的改變。因此，通過監(jiān)測(cè)變壓器的振動(dòng)信號(hào)可以評(píng)估繞組的機(jī)械狀態(tài)。

2 長(zhǎng)圓形繞組機(jī)械故障有限元仿真

在中性點(diǎn)接地方式中，單相繞組故障為變壓器故障的主要類型，占變壓器繞組故障數(shù)量的50%以上。因此，以一臺(tái)SC10-30／10.5干式變壓器為研究對(duì)象，繞組為長(zhǎng)圓形結(jié)構(gòu)，設(shè)計(jì)單相繞組松動(dòng)和翹曲故障。

建立三維模型時(shí)，沿疊片鐵心中心平面取求解域總體的1／2，變壓器長(zhǎng)圓形繞組有限元仿真參數(shù)設(shè)置：額定電流為1.65／43.3A，匝數(shù)比為525／20，相對(duì)介電常數(shù)為3.56，密度為4 450 kg／m3，彈性模量為11 ×103MPa，泊松比為0.35。

數(shù)值計(jì)算了時(shí)間同為0.326 s的A相繞組正常、松動(dòng)和翹曲狀態(tài)下的位移分布，如圖2所示。

圖2 繞組振動(dòng)位移分布

從圖2（a）中可以看出，繞組yoz平面中間位置處（R點(diǎn)）的振動(dòng)位移最大，最大振動(dòng)位移為4. 82 ×10-6mm，表明該處繞組易發(fā)生輻向形變，即該處為繞組薄弱點(diǎn)；由圖2（b）可看出，繞組松動(dòng)時(shí)振動(dòng)位移分布發(fā)生改變，其最大振動(dòng)位移向繞組端部的松動(dòng)處移動(dòng)，最大振動(dòng)位移為7.87 ×10-6mm；圖2（c）為繞組翹曲時(shí)的振動(dòng)位移分布，翹曲后的振動(dòng)位移分布不發(fā)生改變，但振動(dòng)幅值減小，最大振動(dòng)位移為1.36 ×10-6mm。由于繞組發(fā)生松動(dòng)與輻向形變具有累積特性，即當(dāng)繞組存在松動(dòng)和較小形變時(shí)會(huì)導(dǎo)致其抗短路能力下降，經(jīng)過長(zhǎng)時(shí)間運(yùn)行的累積從而更容易發(fā)生嚴(yán)重的機(jī)械故障，因此有必要對(duì)較小的機(jī)械故障進(jìn)行狀態(tài)監(jiān)測(cè)。

3 基于GWO-SVM的長(zhǎng)圓形繞組機(jī)械故障診斷

3.1 振動(dòng)信號(hào)采集

定制SC10-30／10.5干式變壓器的試驗(yàn)樣機(jī)，將繞組外側(cè)的絕緣板去掉，繞組線圈裸露在外。根據(jù)實(shí)際變壓器故障類型及繞組的受力分析，制造了變壓器繞組松動(dòng)和翹曲故障，其結(jié)構(gòu)如圖3所示。

圖3 繞組松動(dòng)與翹曲故障

測(cè)試探頭為356A16型壓電式三軸加速度傳感器，其精度為0.001 m／s2，通過振動(dòng)噪聲分析系統(tǒng)SQuadriga Ⅱ處理采集到的加速度振動(dòng)信號(hào)。由第2 節(jié)仿真分析可知，可以從繞組R點(diǎn)的振動(dòng)特性中提取變壓器機(jī)械故障的特征量，因此將測(cè)試探頭布置在繞組薄弱R點(diǎn)，試驗(yàn)測(cè)試如圖4所示。

圖4 變壓器機(jī)械故障試驗(yàn)測(cè)試

試驗(yàn)時(shí)，用銅導(dǎo)線短接變壓器的低壓繞組，在高壓側(cè)施加激勵(lì)電壓，此時(shí)鐵心振動(dòng)忽略不計(jì)，加速度傳感器檢測(cè)到的振動(dòng)信號(hào)可認(rèn)為全部來(lái)源于繞組。負(fù)載電流Im在（1 ±20%）IN波動(dòng)時(shí)，采集變壓器繞組正常運(yùn)行情況以及繞組松動(dòng)和翹曲故障狀態(tài)下的振動(dòng)信號(hào)各20組。

3.2 小波包能量特征提取

按照能量劃分表示的小波包分解結(jié)果稱為小波包能量譜。小波包能量譜是基于小波包分解系數(shù)求解各個(gè)節(jié)點(diǎn)下的信號(hào)能量，將求解出的能量值按照順序排列組成特征向量，可供變壓器繞組故障識(shí)別使用。各節(jié)點(diǎn)能量計(jì)算公式為

式中M為第j個(gè)頻帶的長(zhǎng)度，dj，i（n）為小波包分解后節(jié)點(diǎn)（j，i）對(duì)應(yīng)的第n個(gè)小波包系數(shù)。

歸一化處理，得到各節(jié)點(diǎn)能量占比D

經(jīng)小波包分解后的第3層能量占比構(gòu)建特征向量T

Symlets小波具備更好的對(duì)稱性，在一定程度上減少信號(hào)重構(gòu)時(shí)產(chǎn)生的相位失真。因此，對(duì)變壓器繞組0～1 000 Hz范圍內(nèi)的振動(dòng)信號(hào)進(jìn)行Symlets三層小波包求解，基于分解系數(shù)求解8個(gè)節(jié)點(diǎn)下的信號(hào)能量，進(jìn)而求解各節(jié)點(diǎn)能量占比。對(duì)變壓器繞組正常、松動(dòng)和翹曲3種不同機(jī)械狀態(tài)下的振動(dòng)信號(hào)進(jìn)行特征提取，得到小波包分解能量譜，如圖5所示。

圖5 繞組小波包分解能量譜

由圖5可知，故障前后的振動(dòng)信號(hào)能量特征變化明顯。繞組正常時(shí)，振動(dòng)能量主要為D30，能量占比高于70%，D31～D37段能量占比在1%～10%不等，其中D33能量占比約為10%；繞組松動(dòng)時(shí)，振動(dòng)能量主要集中在D30～D33頻率區(qū)間，D30占總能量的85%以上，D32次之，能量占比約為10%，相比于繞組正常，繞組松動(dòng)故障導(dǎo)致D30頻率分量占比更大，能量占比向低頻段集中；繞組翹曲時(shí)，振動(dòng)能量占比改變較大，D30能量占比下降至35%左右，相應(yīng)的D34～D37頻率區(qū)間能量占比達(dá)到60%以上。繞組翹曲位置的漏磁場(chǎng)較小且形狀發(fā)生改變，所受電磁力減小，導(dǎo)致獲得的D34～D37振動(dòng)能量大幅增加。

3.3 GWO-SVM算法介紹

2014年，Mirjalili S等人［8］受自然中灰狼群體捕食獵物啟發(fā)，首次提出了GWO算法，該算法具有收斂快，調(diào)參較少和程序簡(jiǎn)潔等特點(diǎn)。

該算法將灰狼群體按等級(jí)依次劃分為α，β，γ和ω。在狼群中，其他灰狼必須聽從和執(zhí)行α狼指令，α狼被稱為支配狼；β狼協(xié)助α狼做出正確的決策，并聽從于α狼；γ狼支配剩余層級(jí)的狼，并聽從α狼和β狼；ω狼等級(jí)最低，服從并跟隨前3等級(jí)的狼對(duì)獵物進(jìn)行追蹤和圍捕，同時(shí)平衡狼群內(nèi)部關(guān)系。

式中Xp（t）為獵物的位置向量，即α，β，γ三個(gè)等級(jí)的最優(yōu)解；X（t）為灰狼的位置向量，即所有可能解。D為灰狼與獵物的距離，X（t＋1）為位置更新公式，A，C為t時(shí)刻位置系數(shù)。

在狩獵過程中，由α，β，γ狼主導(dǎo)尋找獵物，并要求狼群根據(jù)最優(yōu)解更新各自的位置，逐漸向獵物靠近，最終得到的最優(yōu)解。

不同的機(jī)械狀態(tài)對(duì)應(yīng)不同的振動(dòng)特性，可從繞組正常、松動(dòng)和翹曲故障樣本數(shù)據(jù)中各隨機(jī)抽取10 組樣本作為訓(xùn)練樣本，其余的數(shù)據(jù)均作為測(cè)試樣本，進(jìn)行繞組故障診斷。SVM能夠很好地處理小樣本問題，更加適用于變壓器繞組機(jī)械故障的分類問題。本文構(gòu)建了GWO-SVM 的識(shí)別與診斷模型，具體流程如圖6所示。

圖6 GWO-SVM故障診斷流程

3.4 GWO-SVM故障診斷與結(jié)果分析

采用GWO算法優(yōu)化SVM的懲罰參數(shù)c和核函數(shù)參數(shù)g，以獲得模型較高的準(zhǔn)確率。利用GWO算法優(yōu)化SVM獲得的最佳適應(yīng)度曲線如圖7所示。

圖7 GWO尋優(yōu)參數(shù)適應(yīng)度曲線

由圖7可以看出，在適應(yīng)度曲線方面，GWO在第40 次迭代時(shí)能夠達(dá)到最佳與平均適應(yīng)度，最優(yōu)懲罰參數(shù)c為389.1，核函數(shù)參數(shù)g為0. 3，交叉驗(yàn)證的準(zhǔn)確率達(dá)到93.3%。相較于人工隨機(jī)選取參數(shù)而言，GWO進(jìn)行參數(shù)尋優(yōu)有助于節(jié)約時(shí)間、提高準(zhǔn)確率，可以通過利用GWO 優(yōu)化SVM的目標(biāo)函數(shù)，從而使目標(biāo)函數(shù)的準(zhǔn)確率達(dá)到最大。未優(yōu)化與GWO優(yōu)化后的分類結(jié)果如圖8所示，類別標(biāo)簽1、2、3分別表示繞組正常、繞組松動(dòng)和繞組翹曲的實(shí)際狀態(tài)。

圖8 優(yōu)化前后SVM診斷結(jié)果分類

未優(yōu)化SVM的分類結(jié)果如圖8（a）所示，診斷結(jié)果中出現(xiàn)7個(gè)誤診結(jié)果，其中有2 個(gè)“繞組正?！焙? 個(gè)“繞組松動(dòng)”被誤診，未優(yōu)化的SVM分類模型準(zhǔn)確率為76.7%；采用GWO-SVM優(yōu)化的分類結(jié)果如圖8（b）所示，其中僅有1 個(gè)“繞組正常”和2個(gè)“繞組松動(dòng)”被誤診，綜合診斷的準(zhǔn)確率達(dá)到90%?？梢姴捎肎WO-SVM的識(shí)別與診斷模型提高了繞組機(jī)械故障準(zhǔn)確率，因此該模型可應(yīng)用于實(shí)際變壓器長(zhǎng)圓形繞組機(jī)械故障診斷當(dāng)中。

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

研究了長(zhǎng)圓形繞組正常、松動(dòng)和翹曲狀態(tài)下的振動(dòng)位移分布，在不同狀態(tài)下，薄弱位置R點(diǎn)的能量集中程度和分布發(fā)生改變，為實(shí)驗(yàn)研究故障前后變壓器繞組機(jī)械故障的振動(dòng)特性提供分析基礎(chǔ)。

本文采用的GWO-SVM診斷模型在變壓器的機(jī)械故障診斷中具有準(zhǔn)確率高，樣本點(diǎn)少，泛化能力強(qiáng)的優(yōu)勢(shì)，實(shí)驗(yàn)結(jié)果表明，該混合優(yōu)化算法綜合準(zhǔn)確率達(dá)到90%，在變壓器繞組機(jī)械故障診斷中具有一定的應(yīng)用價(jià)值。