轉(zhuǎn)子繞組短路故障時發(fā)電機轉(zhuǎn)子不平衡電磁力分析

周國偉 李永剛 萬書亭 張 玉 李和明

（1.華北電力大學(xué)電力工程系 保定 071003 2.華北電力大學(xué)機械工程系 保定 071003）

1 引言

轉(zhuǎn)子繞組短路是發(fā)電機常見的電氣故障之一 ，造成此故障的原因很多。它將導(dǎo)致轉(zhuǎn)子振動，甚至發(fā)展為轉(zhuǎn)子接地、轉(zhuǎn)子繞組燒損、發(fā)電機失磁、發(fā)電機部件磁化等，危及電機和系統(tǒng)的安全。

在轉(zhuǎn)子繞組短路方面，基于短路后徑向磁通線圈上電壓波形的探測方法在很早就由D.R.Albright提出[1]。然后，在此基礎(chǔ)上進一步得到了以探測線圈的方式來監(jiān)測磁通密度的方法[2-4]，同時利用交流電機繞組理論和多回路理論詳細分析了發(fā)電機轉(zhuǎn)子繞組短路時定子繞組并聯(lián)支路內(nèi)感應(yīng)電動勢的諧波特性、定子繞組感應(yīng)電流產(chǎn)生的旋轉(zhuǎn)磁動勢以及轉(zhuǎn)子繞組感應(yīng)電流的諧波特性，給出了適用于一般發(fā)電機結(jié)構(gòu)的轉(zhuǎn)子繞組短路故障特征規(guī)律[5]。文獻[6]通過分析故障時的電磁特性，提出匝間短路引起勵磁電流增大，但無功卻相對減小或不變的征兆，并基于勵磁電流的變化率診斷轉(zhuǎn)子繞組短路故障。文獻[7，8]則基于勵磁電流變化特征，利用人工神經(jīng)網(wǎng)絡(luò)實現(xiàn)了此方法。文獻[9，10]通過分析發(fā)電機在轉(zhuǎn)子繞組匝間短路、定子繞組匝間短路以及氣隙偏心故障時定子繞組并聯(lián)支路環(huán)流特征，提出了一種基于環(huán)流特征的故障類型辨別診斷方法。文獻[11]建立了汽輪發(fā)電機轉(zhuǎn)子繞組短路故障的多回路數(shù)學(xué)模型，并分析了其中電感系數(shù)的計算方法。

上述方法都沒有考慮故障對發(fā)電機徑向振動的影響，但發(fā)電機作為一個整體，機械與電氣相互耦合，電氣故障會引起氣隙磁場畸變，產(chǎn)生不同于正常運行時的電磁力波，從而激起發(fā)電機徑向電磁振動，如文獻[12，13]所述的發(fā)電機組振動超標現(xiàn)象就是由轉(zhuǎn)子繞組匝間短路引起的。因此研究轉(zhuǎn)子繞組匝間短路引起的發(fā)電機徑向振動特征，將可以在現(xiàn)有研究成果的基礎(chǔ)上，更加全面了解故障對發(fā)電機運行特性的影響，進一步為轉(zhuǎn)子繞組短路故障診斷提供更加全面的征兆。文獻[14，15]詳細分析了發(fā)電機轉(zhuǎn)子繞組短路對發(fā)電機定轉(zhuǎn)子徑向振動特性的影響，指出發(fā)電機徑向振動特征與電氣特征一樣，可作為定子繞組匝間短路故障征兆。

為了深入分析轉(zhuǎn)子繞組短路故障引起轉(zhuǎn)子振動的成因，本文以某火力發(fā)電廠QFR—400—2—20型發(fā)電機轉(zhuǎn)子繞組短路故障為工程背景，提出了一種轉(zhuǎn)子繞組短路故障時轉(zhuǎn)子不平衡電磁力的磁動勢疊加計算方法，并對文獻[16]提出的等效磁通計算方法進行改進。然后分別利用等效磁通計算方法和磁動勢疊加計算方法分析計算了不平衡電磁力，指出了兩種計算方法的特點，為發(fā)電機轉(zhuǎn)子繞組短路故障的現(xiàn)場分析和診斷提供了理論依據(jù)。

2 發(fā)電機主要參數(shù)及故障分析

某電廠為350MW×3的天然氣聯(lián)合循環(huán)電廠，具體機組參數(shù)見下表。

發(fā)電機轉(zhuǎn)子勵磁繞組分布如圖1所示，轉(zhuǎn)子有16個繞組，每極靠近大齒的一個線圈有7匝（即1號和16號線圈），其他14個繞組分別有9匝，轉(zhuǎn)子長度L=6.03m，轉(zhuǎn)子半徑R=0.55m，定轉(zhuǎn)子間氣隙長度δ=0.1m。

表 某電廠發(fā)電機主要參數(shù)Tab.Main parameters of generator

圖1 發(fā)電機轉(zhuǎn)子繞組分布Fig.1 Winding distribution in the generator rotor

該機組2007年2月投運，7#瓦（發(fā)電機汽端瓦）振幅較小，過一階臨界轉(zhuǎn)速時振幅為66μm，過二階臨界轉(zhuǎn)速時振幅為93μm，高負荷時振幅45μm。2007年5月～10月，7#瓦過臨界轉(zhuǎn)速時振幅加大，一階增大至100μm附近，二階增大至100～120μm之間，8#瓦（發(fā)電機勵端瓦）增大卻不大。2008年2月10日～3月13日，7#瓦過一階振動達到330μm，8#則90μm左右，1#～6#瓦的振動和以前比較相差不大；在停機降速過程中（發(fā)電機定子、轉(zhuǎn)子均沒有電流電壓），7#瓦不會發(fā)生類似的振動，但在下降過程中轉(zhuǎn)子加888A勵磁電流，7#瓦過一階振動達到287μm，加200A勵磁電流，7#瓦過一階振動不明顯，為84μm；7#瓦過一階振動均是工頻振動，其重復(fù)性非常好。

根據(jù)上述振動數(shù)據(jù)，在停機降速過一階臨界轉(zhuǎn)速時，勵磁電流對7#瓦處的振動有較大影響。勵磁電流對轉(zhuǎn)子所產(chǎn)生的最主要作用是不平衡電磁力，故初步判斷轉(zhuǎn)子繞組匝間短路所導(dǎo)致的不平衡電磁力是產(chǎn)生振動的主要因素之一。然后采用兩極電壓法判定轉(zhuǎn)子勵端8號線圈存在匝間短路，另外采用RSO法也判定存在轉(zhuǎn)子繞組匝間短路。2008年4月中轉(zhuǎn)子返廠，拔線發(fā)現(xiàn)Ⅱ極8號線圈的4、5匝匝間和5、6匝匝間有兩匝靜態(tài)短路（由于匝間的絕緣層跑位導(dǎo)致匝間短路，見圖2、圖3），Ⅱ極的6號線圈有一匝在轉(zhuǎn)速680～1940r/min時有動態(tài)匝間短路（磨損導(dǎo)致轉(zhuǎn)子的動態(tài)匝間短路，匝間R部有燒焦的痕跡，見圖4）。

圖2 8號線圈的4、5匝匝間Fig.2 Interturn between fourth and fifth turn of No.8

圖3 8號線圈的5、6匝匝間Fig.3 Interturn between fifth and sixth turn of No.8

圖4 6號線圈的3、4匝匝間R部Fig.4 Interturn between third and fourth turn of No.6

3 不平衡電磁力的等效磁通計算方法

日本學(xué)者渡邊孝在文獻[16]中詳細描述了發(fā)電機轉(zhuǎn)子勵磁繞組匝間短路時引起的不平衡電磁力計算方法，并在國內(nèi)許多發(fā)電廠和電機制造廠故障分析中得到推廣應(yīng)用。如圖5所示，在正常時，N極和S極的磁通分布是對稱的，一旦發(fā)生匝間短路，由于存在安匝差，磁通分布便發(fā)生虛線向?qū)嵕€那樣的變化。根據(jù)N極側(cè)和S極側(cè)的磁通相等的原則，假如氣隙磁導(dǎo)沿圓周均勻分布、勵磁線圈在圓周上均勻分布，則磁通密度為零的點β為

即

式中Nt—每極勵磁線圈的總匝數(shù)；

Ns—引起短路的匝數(shù)。

但由于β的存在，N極對應(yīng)的圓心角由π減小為π-2β，S極對應(yīng)的圓心角由π增加為π+2β，即S極短路Ns時，N極對應(yīng)的勵磁線圈總匝數(shù)已不是Nt，將小于Nt。因此文獻[16]中推導(dǎo)的式（1）需進一步修正。

同樣假如勵磁線圈的分布在圓周上是均勻的，由于每極勵磁線圈的總匝數(shù)為Nt，則單位圓周角對應(yīng)的勵磁線圈匝數(shù)為Nt/π。設(shè)S極短路匝數(shù)為Ns，則N極勵磁線圈的總匝數(shù)N1和S極勵磁線圈的總匝數(shù)N2為

式（1）應(yīng)修正為

即（式（2）得到修正）

在圖2中，由于N極和S極的磁通相等，假如其分布為正弦波，那么式（6）成立

式中BN—N極磁通密度；

BS—S極磁通密度。

N極和S極磁通密度可表示為

式中If—勵磁電流；

μ0—空氣磁導(dǎo)率，μ0=4π×10-7H/m；

δ—氣隙長度。

作用于轉(zhuǎn)子的不平衡電磁力可推導(dǎo)為

式中L—轉(zhuǎn)子長度；

R—轉(zhuǎn)子半徑。

QFR—400—2—20發(fā)電機每極勵磁線圈的總匝數(shù)Nt為70匝，其中S極有3個短路匝，分別為汽端Ⅱ極的8號線圈的4、5匝匝間和5、6匝匝間，該兩處是靜態(tài)固定的匝間短路；汽端Ⅱ極的6號線圈的3、4匝匝間R部有燒焦的痕跡，是磨損導(dǎo)致轉(zhuǎn)子的動態(tài)匝間短路。低速時一般只有靜態(tài)匝間短路，因此Ns1=2，而在轉(zhuǎn)子高速旋轉(zhuǎn)時，才發(fā)生動態(tài)匝間短路，因此Ns2=3。

根據(jù)式（5），求得靜態(tài)匝間短路時的β1；將If=888A、δ=0.1m代入式（7），可分別計算N極和S極磁通密度BN1、BS1；然后將L=6.03m、R=0.55m、BN1、BS1代入式（8），求得由于轉(zhuǎn)子8號線圈4、5匝和5、6匝兩處靜態(tài)匝間短路產(chǎn)生的作用于轉(zhuǎn)子的不平衡電磁力為

同理計算當(dāng)轉(zhuǎn)子高速運轉(zhuǎn)時發(fā)生動態(tài)短路時，轉(zhuǎn)子8號線圈兩處靜態(tài)匝間短路和6號線圈一處動態(tài)匝間短路時產(chǎn)生的作用于轉(zhuǎn)子的總不平衡電磁力為W=8 708.2N。

4 不平衡電磁力的磁動勢疊加計算方法

如圖1所示，發(fā)電機正常運行和轉(zhuǎn)子繞組短路故障時，轉(zhuǎn)子磁動勢都關(guān)于縱軸對稱。根據(jù)轉(zhuǎn)子N極和S極的磁通相等的原則，同時假如氣隙磁導(dǎo)沿圓周均勻分布，以縱軸為轉(zhuǎn)子位置角θr的原點，設(shè)Fk(θr)為轉(zhuǎn)子第k個線圈在轉(zhuǎn)子θr角位置時所產(chǎn)生的磁動勢為

式中nk—轉(zhuǎn)子第k個線圈的匝數(shù)；

ak—第k個轉(zhuǎn)子槽對應(yīng)的圓心角。則轉(zhuǎn)子磁動勢為

QFR—400—2—20型發(fā)電機轉(zhuǎn)子有16個繞組，編號如圖1所示，對應(yīng)的圓心角分別為α1=0.362π、α2=0.447π、α3=0.532π、α4=0.617π、α5=0.702π、α6=0.788π、α7=0.873π、α8=0.958π。

轉(zhuǎn)子有16個繞組，每極靠近大齒的一個線圈有7匝（即1和16號線圈），其他14個繞組分別有9匝，繞組2～15線圈匝數(shù)用nk表示，繞組1和16線圈匝數(shù)用nk-2表示。求得轉(zhuǎn)子在每個角度下的磁動勢及合成磁動勢為

第8個線圈有2匝短路，可求得轉(zhuǎn)子在[0，π]之間的磁動勢見式（12）。

根據(jù)轉(zhuǎn)子磁動勢關(guān)于縱軸對稱，當(dāng)轉(zhuǎn)子處在[π，2π]時，轉(zhuǎn)子磁動勢與[0， π]范圍內(nèi)的轉(zhuǎn)子磁動勢相對稱，可以據(jù)此得到發(fā)電機完整的磁動勢波形，如圖6所示。

第6個線圈有1匝短路，可求得轉(zhuǎn)子在[0， π]之間的磁動勢見式（13）。該故障下發(fā)電機的完整磁動勢波形如圖7所示。

第8個線圈有2匝短路、第6個線圈有1匝短路，可求得轉(zhuǎn)子在[0， π]之間的磁動勢為見式（14）。該故障下發(fā)電機的完整磁動勢波形如圖8所示。

圖6 第8線圈2匝短路情況下的磁動勢波形Fig.6 Waveform of magnetomotive force with 2 turns short circuit of No.8

圖7 第6線圈1匝短路情況下的磁動勢波形Fig.7 Waveform of magnetomotive force with 1 turn short circuit of No.6

圖8 第8線圈2匝、第6線圈1匝短路情況下的磁動勢波形Fig.8 Waveform of magnetomotive force with 2 turns short circuit of No.8 and 1 turn short circuit of No.6

主磁動勢產(chǎn)生的單位面積徑向電磁力為

由于分布電磁力關(guān)于縱軸對稱，則Qq=0。

取勵磁電流If=888A，將式（12）代入式（17），求解得到8號線圈4、5匝和5、6匝兩處靜態(tài)匝間短路產(chǎn)生的不平衡電磁力為Qd1=3 176.5N。將式（14）代入式（17），求解得到8號線圈兩處靜態(tài)匝間短路和6號線圈一處動態(tài)匝間短路產(chǎn)生的不平衡電磁力為Qd2=1 0767.4N。

根據(jù)計算結(jié)果，當(dāng)靠近橫軸的轉(zhuǎn)子8號線圈兩處靜態(tài)匝間短路時，等效磁通計算方法求得的不平衡電磁力大于磁動勢疊加計算方法求得不平衡電磁力。當(dāng)靠近橫軸的轉(zhuǎn)子8號線圈兩處匝間短路和偏離橫軸的轉(zhuǎn)子6號線圈一處匝間短路時，等效磁通計算方法求得的不平衡電磁力小于磁動勢疊加計算方法求得不平衡電磁力。等效磁通計算方法是近似公式，近似認為線圈在轉(zhuǎn)子上是均勻分布的，并且沒有考慮短路位置對不平衡電磁力的影響。而根據(jù)式（9）～式（11），轉(zhuǎn)子繞組短路位置對不平衡電磁力影響較大，短路點接近磁極大齒，即靠近縱軸對不平衡電磁力影響較大，靠近橫軸則影響較小，而磁動勢疊加計算方法考慮了短路槽位置的影響。

由于沒有該型發(fā)電機（日立公司設(shè)計）的詳細結(jié)構(gòu)尺寸和電磁參數(shù)，沒能完成有限元分析。有待于條件成熟后對該型發(fā)電機轉(zhuǎn)子繞組匝間短路引起的轉(zhuǎn)子不平衡電磁力進行有限元分析，并與本文的等效磁通法和磁動勢疊加方進行比較分析。

5 結(jié)論

本文通過分析某火力發(fā)電廠發(fā)電機轉(zhuǎn)子繞組短路故障，詳細分析了作用在轉(zhuǎn)子的不平衡電磁力，得出：

（1）在仔細分析文獻[16]提出的轉(zhuǎn)子繞組短路故障時轉(zhuǎn)子不平衡電磁力計算方法的基礎(chǔ)上，提出了一種改進的等效磁通計算方法。

（2）提出了一種轉(zhuǎn)子繞組短路故障時轉(zhuǎn)子不平衡電磁力的磁動勢疊加計算方法。

（3）以某發(fā)電廠QFR—400—2—20型發(fā)電機轉(zhuǎn)子繞組短路故障為背景，分別利用等效磁通計算方法和磁動勢疊加計算方法分析計算了不平衡電磁力。結(jié)果表明等效磁通計算方法沒有考慮短路位置對不平衡電磁力的影響，因此當(dāng)靠近轉(zhuǎn)子橫軸的槽繞組短路時，以等效磁通計算方法計算得到的不平衡電磁力大于以磁動勢疊加計算方法計算得到的不平衡電磁力。當(dāng)偏離轉(zhuǎn)子橫軸的槽繞組短路時，則相反。

[1] Albright D R.Interturn short-circuit detector for turbine-generator rotor windings[J].IEEE Transactions on Power Apparatus and Systems,1971,PAS-90(2):478-483.

[2] Ramirez Nino J,Pascacio A.Detecting interturn short circuits in rotor windings[J].IEEE Computer Applications in Power,2001,14(4): 39-42.

[3] 劉慶河,蔡維錚,徐殿國.汽輪發(fā)電機轉(zhuǎn)子繞組匝間短路在線檢測方法的研究[J].中國電機工程學(xué)報,2004,24(9): 234-237.Liu Qinghe,Cai Weizheng,Xu Dinaguo.Study of on-line detection of interturn short circuit in turbogenerator rotor windings[J].Prodeedings of the CSEE,2004,24(9): 234-237.

[4] 阮羚,周世平,周理兵.大型汽輪發(fā)電機轉(zhuǎn)子匝間短路在線監(jiān)測方法的研究及應(yīng)用[J],中國電機工程學(xué)報,2001,21(supp): 60- 63.Ruan Ling,Zhou Shiping,Zhou Libing.The research and application of the on-line monitor method of large steam generator rotor winding inter turn shortcircuit[J].Proceedings of the CSEE,2001,21(supp):60-63.

[5] 張超,夏立,吳正國,等.同步發(fā)電機轉(zhuǎn)子繞組匝間短路故障特征規(guī)律分析[J].高電壓技術(shù),2010,36(6): 1506-1512.Zhang Chao,Xia Li,Wu Zhengguo,et al.Analysis on fault characteristics law of interturn short circuit in synchronous generator rotor winding[J].High Voltage Engineering,2010,36(6): 1506-1512.

[6] 李永剛,李和明,趙華.汽輪發(fā)電機轉(zhuǎn)子繞組匝間短路故障診斷新判據(jù)[J].中國電機工程學(xué)報,2003,23(6): 112-116,169.Li Yonggang,Li Heming,Zhao Hua.The new criterion on inter turn short-circuit fault diagnose of steam turbine generator rotor winding[J].Proceedings of the CSEE,2003,23(6): 112-116,169.

[7] 李之昆,馬宏忠.基于ANN的發(fā)電機轉(zhuǎn)子繞組匝間短路診斷方法[J].高電壓技術(shù),2004,30(1):31-32,38.Li Zhikun，Ma Hongzhong.Diagnosis for turn-to-turn short circuit of rotor windings in turbo generator[J].High Voltage Engineering,2004,30(1): 31-32,38.

[8] 黨曉強,劉俊勇,雷霞.水輪發(fā)電機轉(zhuǎn)子繞組故障的智能在線識別[J].高電壓技術(shù),2007,33(8): 160-164.Dang Xiaoqiang,Liu Junyong,Lei Xia.Intelligent method of on-line identifying fault of hydro-generator rotor windings[J].High Voltage Engineering,2007,33(8): 160-164.

[9] 黃順禮,黃春陽.定子相繞組橫聯(lián)電勢分析及勵磁繞組匝間短路故障監(jiān)測[J].東方電氣評論,1998,12(3): 142-145.Huang Shunli,Huang Chunyang.Analysis of transverse electromotance in stator phase-windings and monitoring fault of shorted turns in field winding[J].Dongfang Electric Review,1998,12(3):142-145.

[10] 李和明,萬書亭,李永剛.基于定子繞組并聯(lián)支路環(huán)流特性的發(fā)電機故障識別方法[J].電力系統(tǒng)自動化,2005,29(6): 75-78,107.Li Heming,Wan Shuting,Li Yonggang.The generator fault diagnosis method based on stator winding parallel branches circulating current[J].Automation of Electric Power Systems,2005,29(6): 75-78,107.

[11] 李和明,侯紀勇,李俊卿,等.汽輪發(fā)電機勵磁繞組短路多回路數(shù)學(xué)模型及電感系數(shù)[J] .高電壓技術(shù),2010,36(10): 2576-2581.Li Heming.Hou Jiyong.Li Junqing,et al.Multiloop mathematical model and lnductance of lurbo-generator with field winding lnter-turn short circuits[J].High Voltage Engineering,2010,36(10): 2576-2581.

[12] 朱玉壁,洪水盛.平圩電廠1號發(fā)電機10號軸振動分析[J].中國電力,2000,33(10): 45-47.Zhu Yubi,Hong Shuisheng.Vibration analysis of shaft No.10 for generator No.1 in Pingwei power plant[J].Electric Power,2000,33(10): 45-47.

[13] 李鵬,張秀閣,代國超.轉(zhuǎn)子匝間短路引起發(fā)電機組振動的分析及處理[J].華中電力,2008,21(2):8-10,13.Li Peng,Zhang Xiuge,Dai Guochao.Analysis and treatment of vibration of turbo generator unit induced by turn-to-turn short circuit of rotor windings[J].Central China Electric Power,2008,21(2): 8-10,13.

[14] 萬書亭,李和明,李永剛.轉(zhuǎn)子匝間短路對發(fā)電機定轉(zhuǎn)子振動特性的影響[J].中國電機工程學(xué)報,2005,25(10): 122-126.Wan Shuting,Li Heming,Li Yonggang.Analysis of generator vibration characteristic on rotor winding inter-turn short circuit fault[J].Proceedings of the CSEE,2005,25(10): 122-126.

[15] 武玉才,李和明,李永剛.在線檢測發(fā)電機轉(zhuǎn)子繞組匝間短路的新方法[J].高電壓技術(shù),2009,36(12):2698-2704.Wu Yucai,Li Heming,Li Yonggang.New on-line diagnosis method for generator rotor winding inter-turn short circuit[J].High Voltage Engineering,2009,36(12): 2698-2704.

[16] 渡邊孝.汽輪發(fā)電機轉(zhuǎn)軸的熱振動平衡法（熱振動的解析方法）[J].國外大電機,1997(6): 14-21.Watanabe,Takashi.Thermal vibration balancing method for turbine generator rotor (rotor vibration analysis method due to thermal unbalance[J].External Big Motor,1997(6): 14-21.