孫 洋,林文娟,楊作鵬,李興剛
(黑龍江省電力科學(xué)研究院,哈爾濱 150030)
水輪發(fā)電機阻尼繞組能夠穩(wěn)定發(fā)電機轉(zhuǎn)速,削弱負序磁場、抑制轉(zhuǎn)子震蕩[1-2],對水輪發(fā)電機的安全可靠運行起著十分重要的作用。但在水輪發(fā)電機實際運行中,不可避免的會出現(xiàn)不對稱運行狀態(tài)從而產(chǎn)生負序磁場[3-5]。負序磁場會在阻尼繞組中感生很大的渦流,導(dǎo)致阻尼繞組發(fā)熱甚至燒損[6]。因此,對不對稱運行工況下,水輪發(fā)電機阻尼繞組渦流損耗進行深入研究非常必要[7]。
有限元法作為一種高效能的數(shù)值計算方法可以將連續(xù)的求解域離散為一組單元[8-9],它克服了傳統(tǒng)解析法難以全面考慮定子開槽、磁極轉(zhuǎn)動等因素對阻尼繞組渦流損耗的影響[10-11]?;诖?,本文建立不同裝機容量水輪發(fā)電機仿真模型,采用場路耦合有限元法分析了氣隙、節(jié)距、阻尼方式等因素對不同裝機容量水輪發(fā)電機阻尼繞組渦流損耗的影響,以給出有效的應(yīng)對措施,保障機組安全可靠運行。
用于分析計算的三臺水輪發(fā)電機裝機容量分別為225 MW、500 MW和1000 MW,其基本參數(shù)如表1所示。
在不影響計算結(jié)果準確度的情況下,對水輪發(fā)電機數(shù)學(xué)模型進行簡化,取一對磁極作為求解區(qū)域。其建立的有限元模型示例如圖1所示。
采用場路耦合有限元模型求解時,求解場域所滿足的邊值條件為
式中:μ為介質(zhì)磁導(dǎo)率;Az為矢量磁位Z軸分量,它在圓弧邊界AB和CD上分別滿足第一類齊次邊界條件,在直線邊界AD和BC上滿足周期性邊界條件;Jz為源電流密度的Z軸分量;σ為介質(zhì)的電導(dǎo)率;vx為速度的X軸分量。
表1 三臺水輪發(fā)電機基本參數(shù)Tab.1 Basic data of generator
圖1 有限元模型示例Fig.1 Finite element model
導(dǎo)體中的渦流密度地計算表達式為
式中,γ為電導(dǎo)率,A為磁矢位。
一根阻尼繞組渦流為
式中,k為一個磁極阻尼繞組的剖分單元總數(shù),IC為剖分單元的電流,Ie為剖分單元的面積。
每個剖分單元表達式為
中國標準規(guī)定額定容量大于125 MVA的空氣冷卻水輪發(fā)電機能承受的最惡劣不對稱運行工況為9%負序。此時定子繞組中的三相電流為:
9%負序運行時,225 MW水輪發(fā)電機阻尼繞組渦流損耗如表2所示。其中,第4根阻尼繞組渦流損耗最大,為83.2 W;那么經(jīng)計算可知,一個磁極阻尼繞組損耗為193.6 W;整臺水輪發(fā)電機的阻尼繞組總損耗為3.1 kW。
表2 225 MW水輪發(fā)電機阻尼繞組渦流損耗Tab.2 Eddy current loss of 225 MW hydro-generator
9%負序運行時,500 MW水輪發(fā)電機阻尼繞組渦流損耗如表3所示。其中,第6根阻尼繞組渦流損耗最大,為198.5 W;那么經(jīng)計算可知,一個磁極阻尼繞組損耗為525.2 W;整臺水輪發(fā)電機的阻尼繞組總損耗為42.0 kW。
表3 500 MW水輪發(fā)電機阻尼繞組渦流損耗Tab.3 Eddy current loss of 500 MW hydro-generator
9%負序運行時,1000 MW水輪發(fā)電機阻尼繞組渦流損耗如表4所示。其中,第9根阻尼繞組渦流損耗最大,為1310.9 W;那么經(jīng)計算可知,一個磁極阻尼繞組損耗為3880.6 W;整臺水輪發(fā)電機的阻尼繞組總損耗為217.3 kW。
不同裝機容量水輪發(fā)電機阻尼繞組渦流損耗最大值比較分析如圖2所示。
由圖2可知,隨著水輪發(fā)電機額定容量的增加,阻尼繞組渦流損耗不是成倍增長,而是水輪發(fā)電機裝機容量越大,阻尼繞組渦流損耗最大值的增長越多。
表4 1000 MW水輪發(fā)電機阻尼繞組渦流損耗Tab.4 Eddy current loss of 1000 MW hydro-generator
圖2 不同容量水輪發(fā)電機阻尼繞組渦流損耗比較Fig.2 Eddy current loss of different installed capacity
不同容量水輪發(fā)電機單個磁極阻尼繞組平均渦流損耗比較分析如圖3所示。
圖3 不同容量水輪發(fā)電機單個磁極阻尼繞組平均渦流損耗比較Fig.3 Average eddy current loss of a single pole on different installed capacity
從圖3可知,225 MW水輪發(fā)電機單個磁極阻尼繞組平均渦流損耗為48.4 W,500 MW水輪發(fā)電機單個磁極阻尼繞組平均渦流損耗為87.5 W,1000 MW水輪發(fā)電機單個磁極阻尼繞組平均渦流損耗為431.2 W。
由于阻尼繞組內(nèi)渦流損耗與定子開槽產(chǎn)生的齒諧波有關(guān),因此調(diào)整定子、轉(zhuǎn)子氣隙大小可以控制阻尼繞組內(nèi)渦流損耗。下面分別對三臺額定容量不同的水輪發(fā)電機進行分析計算。
225 MW水輪發(fā)電機氣隙大小設(shè)計值為28 mm,調(diào)整氣隙大小±10%,保持其他參數(shù)不變,分別計算氣隙為25.2和30.8 mm時阻尼繞組渦流損耗,計算結(jié)果如表5所示。
500 MW水輪發(fā)電機的氣隙大小設(shè)計值為31 mm,么調(diào)整氣隙大小±10%,分別計算氣隙為27.9 mm和34.1 mm時阻尼繞組渦流損耗,計算結(jié)果如表6所示。
表5 氣隙變化時阻尼繞組渦流損耗Tab.5 Eddy current loss of 225 MW hydro-generator air gap variation
表6 氣隙變化時阻尼繞組渦流損耗Tab.6 Eddy current loss of 500 MW hydro-generator air gap variation
1000 MW水輪發(fā)電機的氣隙大小設(shè)計值為36 mm,調(diào)整氣隙大小±10%,計算氣隙為32.4 mm和39.6 mm時阻尼繞組渦流損耗,計算結(jié)果如表7所示。
表7 氣隙變化時阻尼繞組渦流損耗Tab.7 Eddy current loss of 1000 MW hydro-generator air gap variation
根據(jù)三臺水輪發(fā)電機的計算結(jié)果可知,水輪發(fā)電機氣隙越小,阻尼繞組渦流損耗越大。適當增加氣隙,可以減小阻尼繞組渦流損耗。
為研究節(jié)距對水輪發(fā)電機阻尼繞組渦流損耗的影響,分別建立三臺不同容量機組的水輪發(fā)電機模型。改變其節(jié)距分別為0.5t、0.8t、1.2t。分別計算其阻尼繞著渦流損耗,所得結(jié)果如圖4所示。
圖4 各個容量水輪發(fā)電機阻尼繞組渦流損耗Fig.4 Eddy current loss of different installed capacity
橫向比較圖4中阻尼繞組渦流損耗可知,阻尼繞組節(jié)距越大,那么渦流損耗也越大??v向比較圖4中阻尼繞組渦流損耗可知,當發(fā)電機容量成倍增大時,其渦流損耗成倍增大。而且裝機容量越大,其阻尼繞組渦流損耗受節(jié)距影響也越明顯。
一般水輪發(fā)電機阻尼系統(tǒng)連接方式有兩種,分別為全阻尼和半阻尼。全阻尼連接時,所有磁極的阻尼繞組連接在一起,半阻尼連接時,以磁極為單位的阻尼繞組分別連接在一起。
建立三臺裝機容量的水輪發(fā)電機全阻尼模型,并與半阻尼連接時水輪發(fā)電機阻尼繞組渦流損耗相比較,得出各水輪發(fā)電機組阻尼繞組渦流損耗計算結(jié)果,如圖5所示。
圖5 不同阻尼方式阻尼繞組渦流損耗Fig.5 Eddy current loss of different damping mode
由圖5可知,三臺不同容量水輪發(fā)電機阻尼繞組渦流損耗規(guī)律一致,全阻尼連接方式阻尼繞組渦流損耗大于半阻尼方式阻尼繞組渦流損耗。裝機容量越大,阻尼繞組渦流損耗受阻尼連接方式影響也越明顯。
本文采用場路耦合有限元法,分別對225 MW、500 MW、1000 MW水輪發(fā)電機阻尼繞組渦流損耗進行計算,并通過建立不同仿真模型,分析了氣隙、節(jié)距、阻尼方式對水輪發(fā)電機阻尼繞組渦流損耗的影響,以及采取的應(yīng)對措施。
1)在同一磁極上,沿順時針最后一根阻尼繞組渦流損耗最大。那么這根阻尼繞組也最容易發(fā)熱損壞。
2)水輪發(fā)電機阻尼繞組渦流損耗十分可觀,當裝機容量增大時,阻尼繞組渦流損耗會成倍的增加。
3)氣隙越小,阻尼繞組渦流損耗越大,設(shè)計時應(yīng)適當增大氣隙,減小阻尼繞組渦流損耗。
4)節(jié)距越小,阻尼繞組渦流損耗越小,設(shè)計時應(yīng)適當減小節(jié)距,減小阻尼繞組渦流損耗。
5)半阻尼方式比全阻尼方式阻尼繞組渦流損耗要小,水輪發(fā)電機應(yīng)采用半阻尼方式。
[1] 王國海.三峽右岸全空冷水輪發(fā)電機關(guān)鍵技術(shù)研究[J].中國電機工程學(xué)報,2009,29(15):74 -79.WANG Guohai.Research on key technologies of installed in three gorges right bank station fully air cooled hydro-generators[J].Proceeding of the Csee,2009,29(15):74 -79.
[2] 戈寶軍,李明哲,孫玉田.同步發(fā)電機不對稱運行工況阻尼繞組電流的計算[J].中國電機工程學(xué)報,2013,33(27):154 -160.GE Baojun,LI Mingzhe,Sun Yuntian.Calculation of Damper Winding current for synchronous generators under asymmertric conditions[J].Proceeding of the Csee,2013,33(27):154 -160.
[3] A.I.MEGAHED,O.P.MALIK.Simulation of internal faults in synchronous generators[J].IEEE Trans.on EC,1999,14(4):1306-1311.
[4] V.A.KINISTY.Calculation of internal fault currents in synchronous machines[J].IEEE Trans.on PAS,1965,84(5):381 -389.
[5] MATTIAS WALLIN,OHAN BLADH,URBAN LUNDIN.Damper winding influence onunbalanced magnetic pull in salient pole generators with rotor eccentricity[J].IEEE Transactions on Magnetics,2013,49(9):5158 -5165.
[6] 梁艷萍,孫洋,陳晶,等.1000 MW水輪發(fā)電機不對稱運行時阻尼條渦流損耗計算分析[J].沈陽工業(yè)大學(xué)學(xué)報,2012,34(2):138-143.LIANG Yanping,SUN Yang,CHEN Jing,et al.Calculation and analysis on eddy current loss of damping bar under asymmetrical operation of 1000 MW hydri- generator[J].Journal of Shenyang University of Technology,2012,34(2):138 -143.
[7]王成立,高慧,盧建剛.水輪發(fā)電機定子故障下的阻尼條損耗計算[J].浙江大學(xué)學(xué)報,2012,48(4):770 -776.WANG Chengli,GAO Hui,LU Jiangang.Loss of damper bars in hydro - generator under different stator fault[J].Journal of Zhejiang University,2012,48(4):770 -776.
[8] 李偉力,楊雪峰,顧德寶.空冷汽輪發(fā)電機冷卻氣流風量對定子內(nèi)流體的影響[J].中國電機工程學(xué)報,2009,29(21):53 -61.LI Weili,YANG Xuefeng,GU Debao.Influence of air current flow change on fluid flow and heat transfer of air-cooled turbo-generator with multipath ventilation[J].Proceeding of the Csee,2009,29(21):53-61.
[9] 湯蘊璆.電機內(nèi)電磁場[M].北京:科學(xué)出版社,1981.TANG Yunqiu.The electromagnetic field within electrical machinery[M].Beijing:Science Press,1981.
[10] 李明哲,孫玉田,范吉松,等.大型水輪發(fā)電機穩(wěn)態(tài)負序能力的研究[J].大電機技術(shù),2011(7):59 -62.LI Mingzhe,SUN Yutian,F(xiàn)AN Jisong,et al.The calculation and analysis of the transient negative-seuence capacity of the large hydrogenerator[J].Large Electric Machine and Hydraulic Turbine,2011(7):59 -62.
[11] 孫宇光,王祥珩,桂林,等.場路耦合法計算同步發(fā)電機定子繞組內(nèi)部故障的暫態(tài)過程[J].中國電機工程學(xué)報,2004,24(1):136-141.SUN Yuguang,WANG Xiangheng,GUI Lin,et al.Transient calculation of stator’s internal faults in synchronous generator using FEM coupled with multi- loop method[J].Proceedings of the CSEE,2004,24(1):136 -141.