呂舒藝，李曉華，章李烽，趙文彬，夏能弘，黃堃，楊勇

（1.上海電力大學(xué)電力工程學(xué)院，上海 200090；2.廣州供電局有限公司，廣東廣州 510245）

變壓器作為電網(wǎng)運(yùn)行中重要的設(shè)備之一，其安全穩(wěn)定的運(yùn)行與整個(gè)電力系統(tǒng)的可靠性密切相關(guān)，隨著超高壓直流輸電系統(tǒng)（high voltage direct current，HVDC）在我國的廣泛使用，在變電站附近中性點(diǎn)接地的變壓器中出現(xiàn)了直流偏磁現(xiàn)象[1-2]。直流偏磁會導(dǎo)致變壓器振動異常、噪聲加劇，長期處于偏磁影響下的變壓器發(fā)生繞組變形、松動等故障的概率加大，影響變壓器的穩(wěn)定運(yùn)行[3-4]。因此，有必要對直流偏磁下繞組振動特性進(jìn)行研究和分析。

國內(nèi)外對于直流偏磁下變壓器的偏磁及振動特性研究已有一定基礎(chǔ)。Hassan等人建立了單相變壓器仿真模型，通過諧波特性分析研究了直流偏磁分量對變壓器勵(lì)磁特性的影響[5]。梁貴書等人提出了用兩個(gè)質(zhì)量單元表示繞組的一個(gè)線餅，單元之間由彈簧連接的二維模型。但是該模型將墊塊等效為彈簧時(shí)僅能考慮其彈性特征，未能考慮墊塊的材料特性，同時(shí)無法考慮繞組幾何結(jié)構(gòu)對振動的影響，導(dǎo)致其與繞組實(shí)際振動存在一定誤差[6]。第二類模型是利用有限元進(jìn)行建模分析。羅彤等人采用有限元模型研究了在額定負(fù)載條件下變壓器的電磁特性和振動特性，并通過實(shí)驗(yàn)進(jìn)行了驗(yàn)證[7]。徐永明等人計(jì)算了變壓器在短路時(shí)的短路電動力和在電動力影響下繞組的應(yīng)力特性，得出了繞組不同位置的應(yīng)力分布[8]。但是現(xiàn)有的研究在研究繞組的電動力時(shí)，對于直流偏磁下發(fā)生短路故障時(shí)繞組的振動特性研究較少。變壓器長期處于偏磁影響下時(shí)，即使偏磁電流不大，但長期以往其機(jī)械穩(wěn)定性會受到影響，發(fā)生短路故障時(shí)造成繞組損壞的概率會變大。

文中借助仿真軟件建立了變壓器三維仿真模型，仿真分析了有/無偏磁條件下發(fā)生短路故障時(shí)繞組的短路電流和振動位移特性及其變化規(guī)律，并使用分布反饋光纖激光器（distributed feedback fiber laser，DFB-FL）搭建實(shí)驗(yàn)平臺對模型進(jìn)行了驗(yàn)證。

1 變壓器繞組振動機(jī)理及DFB 激光器測量原理

1.1 繞組振動電磁力分析

當(dāng)變壓器帶負(fù)載運(yùn)行時(shí)，由于繞組中電流的影響，除了在變壓器鐵心中流動的主磁通外，還有漏磁通穿過繞組及結(jié)構(gòu)件。繞組在交變漏磁通的影響下產(chǎn)生交變的電動力，繞組在交變的電動力的作用下會產(chǎn)生振動。振動的大小與漏磁通的大小和繞組中流過電流的大小有關(guān)。

繞組中電流在軸向產(chǎn)生軸向漏磁場Bzt，軸向漏磁場與繞組電流相互作用產(chǎn)生輻向力Fx；同樣，繞組電流在輻向產(chǎn)生輻向漏磁場Bxt，輻向漏磁場與負(fù)載電流相互作用產(chǎn)生軸向力Fz，軸向力向內(nèi)壓縮低壓和高壓繞組，在輻向力的作用下，低壓繞組向內(nèi)壓縮，高壓繞組向外擴(kuò)張。繞組在磁場中的受力如圖1所示。

圖1 繞組在磁場中的受力Fig.1 The force of the winding in the magnetic field

依據(jù)前面對繞組的受力分析可知，作用在繞組上的電磁力可表示為[9]

式中：b為流過繞組的磁密；i為繞組中的電流。

當(dāng)有直流偏磁時(shí)，繞組中會有直流分量流過，同時(shí)激磁電流中還包含其他頻率的諧波分量，此時(shí)i可表示為

式中：I0為直流分量；n為偏磁電流中常見的諧波分量，一般為2 次、3 次、4 次、5 次諧波；ω為外部磁場的變化頻率；φ為相位角。將式（2）代入式（1）可得：

單獨(dú)計(jì)算每次諧波電流產(chǎn)生的F并求和可以得出直流偏磁下變壓器繞組所受到得電磁力表達(dá)式[10]如下式：

式中：Im為電流幅值。

通過對式（4）中等號右側(cè)的第2 項(xiàng)、第3 項(xiàng)分析可知：當(dāng)有直流偏磁時(shí)，由于受到電流中暫態(tài)分量的影響，繞組振動中同時(shí)包含偶倍頻和奇倍頻分量。

1.2 繞組振動位移分析

在繞組的振動位移分析中，將繞組按照由絕緣墊塊層層分開的結(jié)構(gòu)進(jìn)行分析，此時(shí)繞組假設(shè)為質(zhì)量為m的質(zhì)量塊，絕緣墊塊假設(shè)為彈簧，用KB，K，KH表示，繞組間的阻尼系數(shù)為CB，C，CH，繞組的位移用z表示，繞組的受力表示為F。此時(shí)建立出繞組的振動位移方程如下式所示：

由方程可知其中未知數(shù)和方程數(shù)相同，即此時(shí)方程有唯一解，但每個(gè)繞組的阻尼系數(shù)和彈性系數(shù)都不相同，若直接使用會使計(jì)算不易收斂。為了簡化計(jì)算，將繞組的振動看作為整體的振動，即z1=z2=…=zn=z，將方程左右兩邊相加可得：

式中：mˉ為數(shù)個(gè)絕緣墊塊質(zhì)量的平均數(shù)；N為絕緣墊塊個(gè)數(shù)。

將式（6）和上文求出的繞組所受電磁力結(jié)合得到：

1.3 DFB激光器測量原理

傳統(tǒng)的壓電式傳感器易受到電磁干擾，影響測量結(jié)果，并且體積大，不宜放置，因此，文中在進(jìn)行試驗(yàn)驗(yàn)證時(shí)，采用分布式反饋光纖激光器（DFB-FL）作為測量設(shè)備。

DFB 光纖激光器上刻有光柵。光纖光柵周期長度為Λ時(shí)，只有滿足布拉格反射條件波長為λB的光波，才能產(chǎn)生激光震蕩，使激光器得到單頻輸出。由于分布式反饋激光器是由光柵選擇單縱模，因而在高速調(diào)制下仍維持單縱模輸出。Λ的計(jì)算式如下式：

式中：Λ為光柵周期；ne為光纖有效折射率；λB為反射波長。

以往的DFB 激光器多用在水聽器、溫度的檢測中[11-12]，作為振動檢測使用較為少見。圖2是DFB 激光器的測振原理，λn表示反射波長。980 nm 泵浦光源發(fā)出的泵浦光經(jīng)980/1 550 nm波分復(fù)用器（wavelength division multiplexing，WDM）的980 nm 端口傳輸?shù)接蒁FB 激光器組成的傳感頭上，利用WDM 能夠?qū)⒍喾N攜帶各種信息的光載波信號整合和分離的能力，可以將由DFB 激光器產(chǎn)生的激光經(jīng)過980/1 550 nm WDM的1 550 nm 端輸出，通過波長解調(diào)儀進(jìn)行振動數(shù)據(jù)的分析。

圖2 DFB激光器測振原理圖Fig.2 DFB laser measurement vibration schematic

DFB 激光器粘貼在懸臂梁上，當(dāng)懸臂梁振動的時(shí)候，刻寫在DFB 激光器上的光纖光柵的中心波長會移動，從而引起DFB 激光器中心波長的移動，振動越大，中心波長的移動越多，從而可以根據(jù)激光器中心波長移動的大小來判斷振動的大小。

2 有限元模型及振動仿真分析

2.1 有限元模型

文中以一臺匝數(shù)比77/1 785，電壓等級220/5 000 V，額定容量20 kV·A 的單相實(shí)驗(yàn)變壓器為對象進(jìn)行有限元建模分析，主要研究在有/無直流偏磁的影響下變壓器繞組的振動特性。由于真實(shí)的變壓器結(jié)構(gòu)比較復(fù)雜，對模型進(jìn)行了適當(dāng)?shù)暮喕?。忽略了夾件和螺桿等部件，將鐵心和繞組的外部空間設(shè)為空氣，并利用插值函數(shù)的方法定義B—H曲線來定義鐵磁材料，為繞組添加線圈幾何分析并將材料屬性設(shè)定為銅，繞組首末兩端設(shè)置固定約束來模擬實(shí)際情況。絕緣墊塊的彈性模量隨其受到的力的變化而變化，所以在材料設(shè)置中將絕緣墊塊的材料屬性定義為非線性彈性材料。為了提高整體仿真的精確度，對變壓器不同的部位采用不同的網(wǎng)絡(luò)結(jié)構(gòu)，需要重點(diǎn)分析的繞組部位采用較細(xì)化的網(wǎng)格，其他部分采用常規(guī)網(wǎng)格。

為了提高運(yùn)算效率，繞組共分6層，初級和次級繞組每層初始墊塊數(shù)量為4個(gè)。整體模型如圖3所示，共104 515個(gè)域單元。

圖3 變壓器有限元模型Fig.3 Finite element model of transformer

在以往對變壓器繞組振動的研究中，對于在直流偏磁的影響下發(fā)生短路故障時(shí)繞組振動的研究較為少見。文中通過有限元模型和外電路的耦合，建立由低壓、高壓繞組和交流、直流源組成的等效電路，其中直流偏磁通過直流源引入，外電路由交、直流源疊加組成，同時(shí)設(shè)置單相變壓器仿真模型在40 ms 時(shí)發(fā)生短路故障，研究有/無偏磁情況下發(fā)生短路故障時(shí)，繞組振動特性的變化。偏磁仿真電路圖和繞組振動有限元分析流程圖如圖4、圖5所示。

圖4 偏磁仿真電路圖Fig.4 Bias magnetic simulation circuit diagram

圖5 繞組振動有限元分析流程圖Fig.5 Windings vibration finite element analysis flow chart

電場模塊方程為

式中：ε0為自由空間的介電常數(shù)；εr為相對介電常數(shù)；σ為電導(dǎo)率；Je為外部電流密度；V為電勢。

將圖4外電路中計(jì)算出來的電流作為激勵(lì)引入到變壓器模型中的線圈中，實(shí)現(xiàn)電場和磁場的耦合，磁場模塊的方程為

式中：μ0為自由空間的磁導(dǎo)率：μr為相對磁導(dǎo)率：A為矢量磁位。

將電磁場模塊中求得的b，Je等參數(shù)作為初值代入到式（10）中，實(shí)現(xiàn)電磁場和結(jié)構(gòu)場的二次耦合，計(jì)算得到繞組的加速度和應(yīng)力應(yīng)變，結(jié)構(gòu)場的求解域方程為

式中：md2u/dt2為繞組慣性力；ζdu/dt為繞組阻尼力；ku為繞組絕緣彈性力；Ft為繞組受到的隨時(shí)變化的電磁力。

絕緣墊塊的應(yīng)力-應(yīng)變關(guān)系可表示為[17]

式中：σ，ε為絕緣墊塊的應(yīng)力、應(yīng)變；ao，bo為線性系數(shù)和硬化系數(shù)。

由式（12）可推出墊塊彈性模量的計(jì)算公式為

2.2 仿真結(jié)果及分析

變壓器在突發(fā)短路故障時(shí)，短路電流的峰值可達(dá)額定電流的數(shù)十倍之多。為了研究突發(fā)短路故障時(shí)繞組的振動特性，文章在軟件中設(shè)置單相變壓器仿真模型在40 ms 時(shí)發(fā)生短路故障，分析在有/無偏磁情況下發(fā)生短路故障時(shí)繞組的振動特性。在實(shí)際情況中繞組上、下部都有緊固件進(jìn)行約束，文中在仿真軟件中忽略了上、下部的緊固件，通過在繞組上、下底面施加固定約束來模擬緊固件的作用，使仿真和實(shí)際更加貼近。

圖6 是有/無偏磁情況下繞組發(fā)生短路故障時(shí)二次側(cè)的電流波形。從圖6 中可以看出，在無偏磁時(shí)短路電流的峰值約為正常額定電流的17倍，在偏磁情況下短路電流在達(dá)到第1個(gè)峰值時(shí)會發(fā)生偏移，隨著偏磁電流的增大，偏移也隨之增大，經(jīng)過約6個(gè)周期后波形衰減為穩(wěn)態(tài)短路電流。

圖6 短路時(shí)二次側(cè)電流波形Fig.6 Secondary side current waveforms during short circuit

由于不同偏磁電流下短路電流峰值的偏移較整體來說不易從圖中看出變化趨勢，將不同偏磁情況下短路電流的峰值單獨(dú)取出作圖，如圖7所示。從圖7中可知，隨著偏磁電流的增加，短路電流峰值呈現(xiàn)出逐步穩(wěn)定增加的趨勢，當(dāng)偏磁電流接近4 A 時(shí)，增加逐漸平緩，這是因?yàn)槲闹醒芯康膶?shí)驗(yàn)變壓器容量較小，所能承受的偏磁電流幅值較低，在偏磁電流達(dá)到4 A 左右時(shí)，就已接近所允許通入的偏磁電流幅值。當(dāng)偏磁達(dá)到2 A 時(shí)，短路電流峰值較無偏磁情況下增加了約5.43%。

圖7 不同偏磁電流下短路電流峰值Fig.7 Short-circuit current peak under different bias currents

圖8是不同偏磁電流下短路電流各次諧波的含量及變化趨勢，電流基頻為50 Hz。由圖8 可知，短路電流諧波頻譜中2 次、3 次、4 次諧波含量較多，同時(shí)含有少量其他倍頻諧波分量。

圖8 不同偏磁電流下短路電流各次諧波Fig.8 Short-circuit current harmonics under different bias currents

由于在繞組端部漏磁較大，繞組會受到較大的短路電動力，所以以第1餅為例，分析此處在有/無偏磁情況下發(fā)生短路故障時(shí)所受到的應(yīng)力，結(jié)果如圖9 所示。由圖9 可知，在發(fā)生短路故障時(shí)，繞組受到巨大的短路電動力，其中包含暫態(tài)分量和穩(wěn)態(tài)分量，隨著暫態(tài)分量的逐漸衰減，短路電動力逐漸趨于穩(wěn)定。在不同的偏磁電流下，繞組受到的短路電動力最大值變化如圖10 所示。由圖10可知，隨著偏磁電流的增加，繞組受到的最大短路電動力逐步增加，當(dāng)偏磁2 A 時(shí)，較無偏磁情況下的最大短路電動力增加了約6.68%。

圖9 短路沖擊下低壓繞組1餅處的應(yīng)力分布Fig.9 Stress distribution of first low voltage winding under short-circuit impulse

圖10 不同偏磁電流下短路電動力最大值Fig.10 Short-circuit electric power maximum under different bias currents

圖11 是不同偏磁電流下繞組短路電動力的諧波分布。由圖11 可知，在不同的偏磁電流下，短路電動力各次諧波中150 Hz，200 Hz，250 Hz 等諧波分量含量較多，同時(shí)增長趨勢較其他諧波分量也更明顯。

圖11 不同偏磁電流下短路電動力各次諧波Fig.11 Short-circuit electric dynamic harmonic under different bias currents

圖12 是有/無偏磁情況下發(fā)生短路故障時(shí)低壓繞組1 餅處的振動位移波形。由圖12 可知，當(dāng)發(fā)生短路故障時(shí)繞組振動出現(xiàn)更加復(fù)雜的變化特性，位移特性和電流及應(yīng)力一樣，也存在暫態(tài)分量和穩(wěn)態(tài)分量，在短路后約10 ms 時(shí)達(dá)到最大值，隨著暫態(tài)分量的衰減逐漸趨于穩(wěn)定。

圖12 短路沖擊下低壓繞組1餅處的位移波形Fig.12 Displacement waveforms of the low-voltage winding 1 under the short-circuit impact

在不同偏磁電流下突發(fā)短路故障時(shí)，位移最大值的變化趨勢不易從波形圖中直接看出，故將不同偏磁電流下突發(fā)短路故障時(shí)繞組位移的最大值取出分析，如圖13 所示。由圖13 可知，繞組最大位移在偏磁電流的影響下呈現(xiàn)出和應(yīng)力相似的變化趨勢。當(dāng)偏磁2 A 時(shí)，繞組位移最大值增加了約6.75%。

圖13 不同偏磁電流下繞組位移最大值Fig.13 Maximum winding displacement under different bias currents

圖14 是不同偏磁電流下低壓繞組1 餅處的振動位移頻譜分布。由圖14可知，其諧波分量分布和應(yīng)力的諧波分量相似，包含150 Hz，200 Hz，250 Hz 等諧波分量，并且隨著偏磁電流的增加，有著較為明顯的增長趨勢。

圖14 不同偏磁電流下低壓繞組1餅處的振動頻譜Fig.14 Vibration spectrum of cake at low voltage winding 1 under different bias currents

通過對直流偏磁下變壓器繞組振動的仿真結(jié)果分析可知，短路電流、短路電動力、振動位移都包含暫態(tài)分量和穩(wěn)態(tài)分量，當(dāng)偏磁電流達(dá)到2 A 時(shí)，最大短路電流增加了5.43%，最大短路電動力增加了6.68%，隨著偏磁電流的增加，短路電流頻譜中除了占據(jù)主要分量的基波外，還包含部分的100 Hz，150 Hz 及少量其他倍頻諧波分量，短路電動力除了基頻分量外還包含部分150 Hz，200 Hz，250 Hz 等諧波分量。雖然直流偏磁造成的短路電流和振動位移偏移量較整體而言沒有很大，但是變壓器在實(shí)際運(yùn)行中長期受到偏磁的影響時(shí)，因其包含各種諧波分量，會造成變壓器振動更加復(fù)雜，長期以往會對機(jī)械穩(wěn)定性和使用壽命造成影響。

3 試驗(yàn)研究

為了驗(yàn)證的方便，選擇無偏磁情況下繞組的加速度仿真結(jié)果和實(shí)驗(yàn)結(jié)果進(jìn)行對比。以一臺容量20 kV·A、電壓等級5 000/220 V 的實(shí)驗(yàn)變壓器為研究對象建立了變壓器振動測試平臺。由于傳統(tǒng)的壓電式傳感器容易受到電磁場等外部干擾的影響，本次試驗(yàn)采用DFB 激光器進(jìn)行測試。實(shí)驗(yàn)平臺由高靈敏度DFB 光纖振動傳感器、傳輸光纖、解調(diào)儀、筆記本組成。實(shí)驗(yàn)裝置如圖15所示。

圖15 振動測試裝置現(xiàn)場圖Fig.15 Vibration test device site map

以仿真模型中繞組1/5 處的振動為例，仿真此位置振動加速度及頻譜分析，如圖16所示。在試驗(yàn)時(shí)將繞組低壓側(cè)短路，測得的加速度及頻譜分析如圖17所示。

圖16 繞組加速度仿真值Fig.16 Winding acceleration simulation value

圖17 繞組加速度實(shí)測值Fig.17 Winding acceleration measured value

對比加速度的仿真和實(shí)測值可知，頻譜結(jié)果都是以100 Hz 為基準(zhǔn)值，并且占主要分量，這與已有結(jié)論一致。仿真頻譜中含有少量的低次諧波信號，主要都分布在500 Hz 以下，這是由于仿真模型中考慮了材料的非線性，使得加速度曲線中含有毛刺所致。仿真和實(shí)驗(yàn)加速度曲線的幅值和振動波形基本一致，總體上仿真結(jié)果和實(shí)驗(yàn)結(jié)果吻合較好，證明了所建的變壓器繞組仿真模型的準(zhǔn)確性。

4 結(jié)論

文中借助于有限元軟件對不同偏磁情況下變壓器發(fā)生短路故障時(shí)的繞組振動進(jìn)行了仿真分析，結(jié)果表明：

1）所搭建的變壓器三維仿真模型基于多物理場的耦合，能夠有效實(shí)現(xiàn)不同偏磁電流下發(fā)生短路故障時(shí)變壓器繞組振動的仿真。

2）在不同偏磁電流下發(fā)生短路故障時(shí)，短路電流及短路電動力瞬間增大，且隨著暫態(tài)分量的衰減逐漸平穩(wěn)。繞組振動特征呈現(xiàn)出更加復(fù)雜的變化，其頻譜分量中包含150 Hz，200 Hz，250 Hz等低次諧波分量。

3）通過試驗(yàn)和仿真結(jié)果的對比，驗(yàn)證了仿真模型及分析方法的正確性，文中的模型對分析直流偏磁下變壓器繞組的振動和變壓器的在線監(jiān)測有一定的借鑒意義。