地磁感應(yīng)電流作用下大型變壓器的溫升特性計算

朱 濤 王豐華

地磁感應(yīng)電流作用下大型變壓器的溫升特性計算

朱 濤1王豐華2

（1. 上海電力大學(xué)電氣工程學(xué)院 上海 200090 2. 上海交通大學(xué)電氣工程系 上海 200240）

為準(zhǔn)確分析電力變壓器在地磁感應(yīng)電流（GIC）作用下的溫升特性，該文通過建立某500kV三相共體油浸式電力變壓器的電磁-流熱場耦合三維仿真模型，計算變壓器滿負載運行及不同GIC作用下其內(nèi)部構(gòu)件的損耗密度和溫度場分布，并依據(jù)IEEE C57.163標(biāo)準(zhǔn)，計算GIC作用下變壓器內(nèi)部構(gòu)件的動態(tài)溫升曲線。仿真計算結(jié)果與設(shè)計值的良好吻合說明了該文計算結(jié)果的正確性。同時發(fā)現(xiàn)，受GIC影響后，電力變壓器的夾件、拉板、油箱等的渦流損耗分布極不均勻，且隨GIC電流的增加而增大。GIC脈沖作用下的變壓器金屬構(gòu)件局部溫度會急劇上升，在100A直流脈沖下，變壓器主柱拉板的局部溫升達到允許限值。研究結(jié)果可為GIC作用下大型變壓器的耐受性能分析提供重要參考依據(jù)。

大型變壓器 地磁感應(yīng)電流 損耗密度 溫度場 溫升曲線

0 引言

地磁感應(yīng)電流（Geomagnetically Induced Current, GIC）是地磁場發(fā)生強烈擾動（磁暴）時的強感應(yīng)電場在導(dǎo)電體內(nèi)產(chǎn)生的感應(yīng)電流，持續(xù)時間為幾秒到幾小時。由于GIC的頻率小于0.01Hz，對變壓器而言，可將其視為緩慢變化的直流[1]，故與高壓直流輸電引起的變壓器直流偏磁現(xiàn)象類似，GIC在變壓器中的流動會導(dǎo)致鐵心半周飽和及漏磁增大，引起變壓器金屬構(gòu)件的渦流損耗增加，在一定條件下有可能會使變壓器內(nèi)部的局部溫升超過允許的限制[2]。而溫度是影響變壓器壽命的關(guān)鍵因素之一，即變壓器的溫度越大，其預(yù)期壽命就會越短，有必要準(zhǔn)確掌握GIC作用下變壓器的溫升特性，提高其運行可靠性。

GIC通常多發(fā)于高緯度地區(qū)，國內(nèi)外研究者分別從GIC電流水平監(jiān)測及其對電力系統(tǒng)的影響、GIC對電力變壓器的勵磁電流、無功損耗、振動、噪聲、溫升等的影響及抑制策略等方面進行了若干研究[3-6]。其中，對GIC作用下電力變壓器的溫升特性研究主要包括試驗研究和數(shù)值計算等。文獻[7]選取了最易受直流偏磁影響的單相三柱式變壓器進行溫升試驗，發(fā)現(xiàn)其在200AGIC電流作用下的夾件局部溫升約110℃。文獻[8]給出了某110kV單相單柱旁軛變壓器空載和負載試驗中振動、噪聲、繞組熱點溫升和結(jié)構(gòu)件溫升、油色譜等的測試結(jié)果，發(fā)現(xiàn)直流偏磁下變壓器繞組熱點溫升和結(jié)構(gòu)件溫升顯著增大，當(dāng)直流達到36A時，低壓繞組熱點溫升和拉板上部最高溫升均超過90K，且油中出現(xiàn)碳氫氣體。試驗研究雖能夠較貼切實際，但受制于試驗條件如電源容量限制、試驗周期和成本等諸多因素影響，研究者更多的是借助于數(shù)值計算或與試驗研究結(jié)合的方法探究變壓器的溫升特性。根據(jù)變壓器的有限元分析模型，文獻[9]依據(jù)GIC作用下不同鐵心結(jié)構(gòu)的變壓器勵磁電流和磁通分布計算結(jié)果，計算了自耦變壓器鐵心螺栓和油箱夾件的損耗和溫升，發(fā)現(xiàn)鐵心螺栓存在產(chǎn)生局部過熱的風(fēng)險，且變壓器溫升與直流持續(xù)時間線性相關(guān)。文獻[10]計算分析了GIC作用下變壓器的油箱損耗特性及其影響因素（變壓器類型、功率因數(shù)、并聯(lián)磁路等），認為降壓變壓器的油箱損耗大于升壓變壓器，且在功率因數(shù)為0.9時變壓器油箱的損耗和溫升將會最嚴(yán)重。文獻[11-12]計算了GIC作用下單相自耦變壓器的損耗和溫升變化規(guī)律，指出變壓器鐵心以及鄰近夾件在

GIC作用下溫升較小，但拉板處是最容易快速溫升的變壓器組件。文獻[13]依據(jù)GIC作用下變壓器油箱三維模型的損耗和溫升計算結(jié)果，指出油箱前壁中心處為其損耗最大值和熱點位置，同時依據(jù)變壓器內(nèi)部油的熱對流速度確定了其溫升速率。文獻[14-15]基于某500kV單相變壓器的電路-磁路模型的計算結(jié)果，發(fā)現(xiàn)變壓器內(nèi)部結(jié)構(gòu)件的損耗隨著直流偏磁程度的增加而增大，當(dāng)鐵心飽和時，結(jié)構(gòu)件的渦流損耗增幅將加劇。文獻[16]根據(jù)某殼式變壓器的三維有限元模型分析了變壓器的繞組電流、損耗密度分布及溫升變化，同時指出，如果變壓器采用非磁性材料金屬構(gòu)件，可長時間承受200A直流和短時間承受400A直流電流影響，并通過試驗進行了驗證。綜上，現(xiàn)有研究大都對變壓器空載運行下的損耗進行計算分析，對受GIC影響的變壓器溫升特性的研究尚不多見，且實際GIC為低至中等水平直流下伴隨一些短時間內(nèi)較高的峰值，其對變壓器溫升特性的影響程度尚不清楚，需要進一步研究。

我國雖處于中低緯度地區(qū)，但隨著太陽活動的增強及我國超/特高壓輸電工程的快速發(fā)展，GIC對電力系統(tǒng)安全穩(wěn)定運行的威脅也越來越大。為進一步理解和掌握GIC作用下的變壓器溫升特性，本文以某500kV三相共體變壓器為對象，基于電磁場和流熱場耦合模型，計算其在不同GIC條件下滿負載運行時的損耗和溫升特性，據(jù)此分析評估變壓器的GIC耐受能力，并給出降低變壓器GIC溫升的防范措施建議，為大型主變壓器的安全可靠運行提供參考。

1 變壓器電磁-流熱場計算模型

1.1 變壓器的電磁計算模型

針對GIC作用下的變壓器電磁特性，本文在此采用場路耦合模型進行計算，分別描述如下。

1.1.1 電路方程

對高電壓等級和大容量的有載調(diào)壓自耦變壓器來說，因需要較多電壓檔位進行調(diào)壓，通常采用正反調(diào)壓的調(diào)壓方式。相應(yīng)地，可得GIC或直流電流作用下此類三相自耦變壓器的電路模型如圖1所示。圖中，A、B和C為三相交流電壓源；A1、B1和C1為三相串聯(lián)繞組（高壓）漏感；A2、B2和C2為三相公共繞組（中壓）漏感；A3、B3和C3為三相調(diào)壓繞組漏感；A、B和C為三相負載；DC為直流偏置激勵。此處，通過在變壓器中性點施加直流來進行不同GIC作用下的變壓器電磁特性計算。

圖1 500kV變壓器直流偏磁時電路模型

依據(jù)圖1所示的變壓器電路，可列寫電路方程為

式中，A1、B1和C1為三相串聯(lián)繞組電阻；A2、B2和C2為三相公共繞組電阻；A3、B3和C3為三相調(diào)壓繞組電阻；A1、B1和C1為三相串聯(lián)繞組匝數(shù)；A2、B2和C2為三相公共繞組匝數(shù)；A3、B3和C3為三相調(diào)壓繞組匝數(shù)；A、B和C為三相磁通；A1、B1和C1為三相串聯(lián)繞組電流；A2、B2和C2為三相公共繞組電流；A3、B3和C3為三相負載電流。

1.1.2 磁場方程

在變壓器三維交變電磁場中，由麥克斯韋方程可得到其矢量磁位方程為

依據(jù)式（2）可得到變壓器的磁場分布。此外，變壓器的電路方程計算得出的激勵電流i以外部電流密度的形式耦合到磁場方程中，然后通過硅鋼片的磁化曲線（曲線）得到的鐵心交變磁通反饋到電路方程實現(xiàn)場路耦合，其存在對應(yīng)關(guān)系為

式中，為磁感應(yīng)強度；為磁場強度；為電流密度；為電場強度。

聯(lián)合求解式（1）～式（3），可得出直流偏磁時500kV變壓器的勵磁電流、負載電流和磁通分布等，從而為變壓器的損耗及溫升計算提供數(shù)據(jù)支持。

1.2 變壓器的損耗計算

變壓器損耗包括鐵損、銅損和雜散損耗等，其中，鐵損主要由磁滯損耗和渦流損耗組成，有

銅損由直流電阻損耗和渦流損耗組成，有

式中，W為繞組總損耗；R為直流電阻損耗；E為繞組渦流損耗；av為繞組電流密度有效值；為導(dǎo)線截面積；為繞組電阻。

大型變壓器繞組的渦流損耗能達到繞組總損耗的10%以上[17]，特別在受GIC影響后渦流損耗會明顯增大，因此不能忽略。對于漏磁引起的繞組渦流損耗分為徑向和軸向兩部分，其中，繞組不同部位的渦流損耗[18]為

變壓器的雜散損耗主要指油箱、夾件、拉板等金屬部件在漏磁通下產(chǎn)生的渦流損耗。由于變壓器金屬結(jié)構(gòu)不規(guī)則且磁通密度分布不均勻，利用已有解析公式和半經(jīng)驗公式計算的損耗誤差較大，因此，本文在此通過式（2）計算得到的渦流密度可求得各單元的渦流損耗為

式中，e為渦流損耗。

1.3 變壓器的流熱場計算模型

油浸式變壓器運行中的損耗所轉(zhuǎn)化的熱量通過絕緣油的循環(huán)對流進行散熱冷卻，通常將變壓器油視為不可壓縮性流體，認為其流動遵循質(zhì)量守恒定律、動量守恒定律和能量守恒定律，其在直角坐標(biāo)下可分別表示為

變壓器的鐵心、繞組、油箱和夾件等通過熱傳導(dǎo)使熱量從溫度高的地方傳至溫度較低的地方，故熱傳導(dǎo)方程為

式中，為物體密度；為物體比熱容；為單位時間單位體積熱源的發(fā)熱量。

在固液交界面上，固體和液體的熱耦合是靠表面?zhèn)鳠徇M行的，表面?zhèn)鳠岬挠嬎愎綖?/p>

式中，為發(fā)熱量；為表面?zhèn)鳠嵯禂?shù)；w為固體表面溫度；f為流體溫度。

變壓器內(nèi)外部只要有溫度差，就會有熱輻射，物體表面的輻射程度是受表面輻射率影響的，輻射的熱轉(zhuǎn)移斯特凡玻耳茲曼定律為

式中，RA為輻射熱量；為斯特凡玻耳茲曼常數(shù)；R為表面發(fā)射率；R為表面積；S為表面平均溫度；a為周圍環(huán)境溫度。

2 變壓器的電磁-流熱場有限元建模

以一臺型號為OSFSZ-525000/525的三相五柱式自耦變壓器為研究對象，額定容量為525MV·A/ 525MV·A/105MV·A，額定電壓為525kV/345kV/ 34.5kV，空載電流百分數(shù)為0.03%，阻抗電壓百分數(shù)為：高/中6.9%；高/低72.6%；中/低63.5%，繞組匝數(shù)為398/762/132，冷卻方式為自然油循環(huán)。因?qū)嶋H變壓器繞組線餅數(shù)較多且機械結(jié)構(gòu)較為復(fù)雜，且散熱器為內(nèi)部油循環(huán)必不可少的一部分，在建模計算時進行簡化和假設(shè)：

（1）將圓筒繞組進行分層以更為接近繞組線餅和油道的結(jié)構(gòu)，用箱體代替片式散熱器，忽略墊塊、絕緣紙筒等絕緣件。

（2）忽略鐵磁材料的磁滯損耗。

（3）繞組電流密度在導(dǎo)體內(nèi)認為均勻分布。

（4）考慮到非平均損耗密度作為流熱模型熱源的導(dǎo)入，磁路模型和流熱模型的結(jié)構(gòu)應(yīng)保持一致。

圖2為所建立的變壓器三維有限元模型。當(dāng)電路模型計算得到的激勵電流施加到變壓器繞組上時，會產(chǎn)生線圈直流電阻損耗，同時部分漏磁會在變壓器結(jié)構(gòu)件中感應(yīng)出渦流電流，引起結(jié)構(gòu)件發(fā)熱。在電磁暫態(tài)計算模型中邊界條件設(shè)定如下：

（1）磁場求解域邊界為平行邊界條件，即磁動勢的切向分量在邊界處設(shè)置為零。

（2）電磁場計算所用材料屬性為各向同性，見表1，鐵心-曲線如圖3所示。

圖2 變壓器三維有限元模型

表1 電磁場材料屬性

Tab.1 Electromagnetic field material properties

圖3 鐵心B-H曲線

因自然油循環(huán)變壓器是一內(nèi)部油流為自然對流的封閉容器，受熱膨脹上升后的油從油箱出口流出經(jīng)過外部冷卻裝置散熱后從油箱入口流入。在流熱物理模型中邊界條件設(shè)定如下：

（1）油箱外壁表面?zhèn)鳠嵯禂?shù)設(shè)為15W/(m2·K)，環(huán)境溫度恒定為20℃，散熱器表面?zhèn)鳠嵯禂?shù)為143W/(m2·K)。

（2）鐵心熱源以總損耗形式施加，其余熱源以非平均損耗密度形式施加，流場采用層流模型，流熱場模型采用瞬態(tài)求解。

（3）變壓器固體材料屬性為常數(shù)，絕緣油材料屬性為隨溫度變化的擬合函數(shù)，分別見表2和表3。

表2 流熱場固體材料屬性

Tab.2 Solid material properties

表3 流熱場絕緣油材料屬性

Tab.3 Insulating oil material properties

變壓器多物理場有限元模型包括磁、溫度場單元實體類型和流場單元類型，各單元離散化程度選擇線性單元類型，計算網(wǎng)格要同時適用于電磁和熱流場物理屬性特征，因此網(wǎng)格要盡量規(guī)則且單元大小和增長率不宜過大。變壓器網(wǎng)格剖分如圖4所示，對繞組、鐵心柱、散熱器、夾件等采用掃掠網(wǎng)格剖分，其中繞組網(wǎng)格為六面體單元類型，絕緣油網(wǎng)格主要為四面體單元類型。在計算GIC作用下變壓器的溫升特性時，首先通過場路耦合計算出變壓器滿負載運行不同GIC影響時的損耗，然后將變壓器各

圖4 變壓器網(wǎng)格剖分

部件的損耗導(dǎo)入熱流場模型作為傳熱接口的熱源計算其溫升，流熱場對電磁場無影響，該過程為順序耦合。據(jù)此可得本文的多物理場耦合計算過程如圖5所示。此處，場路模型采用的是直接耦合法，熱場中的熱源對流體加熱使流體受熱膨脹產(chǎn)生流動，流場流速的大小影響熱場表面?zhèn)鳠岬拇笮?，因此流熱場模型亦采用直接耦合法?/p>

圖5 電磁-流熱場耦合計算過程

3 結(jié)果分析

3.1 變壓器損耗

分別對變壓器額定負載下和GIC作用下的變壓器損耗計算結(jié)果進行分析。

3.1.1 額定負載下變壓器電磁及損耗結(jié)果

由于在高中運行調(diào)壓最小分接的工況下變壓器損耗最大，所引起的溫升最為嚴(yán)重[19]，因此選擇此運行方式進行損耗和溫升的仿真計算。

圖6為變壓器高中最小分接額定負載運行時的A相電流和磁感應(yīng)強度計算結(jié)果。其中，圖中所給磁感應(yīng)強度云圖對應(yīng)于B相繞組所在的鐵心柱的磁感應(yīng)強度達最大值時刻。此處，根據(jù)變壓器銘牌參數(shù)，其對應(yīng)的高壓繞組、中壓繞組和負載的電流有效值分別為577.4A、456.2A和1 033.6A，且變壓器鐵心的磁感應(yīng)強度約在1.7T左右。顯然，計算結(jié)果與變壓器設(shè)計參數(shù)吻合良好，說明了計算結(jié)果的準(zhǔn)確性。

表4為變壓器高中最小分接額定負載運行時各部件的總損耗計算結(jié)果與廠家的相關(guān)設(shè)計值。由表4可見，本文計算結(jié)果與設(shè)計值吻合良好，說明了本文計算模型及計算結(jié)果的合理性和準(zhǔn)確性。其中，漏磁分布的計算誤差是繞組損耗結(jié)果產(chǎn)生偏差的主要原因。同時，由于在有限元建模時未考慮油箱加強筋以及箱蓋和側(cè)壁等連接處為非直角連接等因素，致使油箱損耗計算結(jié)果偏高。此外，因為未考慮連接螺栓等產(chǎn)生的附加損耗、拉板和夾件的磁滯損耗等原因?qū)?dǎo)致鐵心、夾件和拉板損耗導(dǎo)致的計算結(jié)果偏低。

圖6 變壓器的電流和磁感應(yīng)強度

表4 變壓器各部件的損耗計算結(jié)果

Tab.4 Total loss of several components of transformer

圖7為電力變壓器的雜散損耗密度分布?？梢?，變壓器油箱渦流損耗主要集中分布在繞組上下方和兩側(cè)，頂部損耗密度最高達66kW/m3，最大損耗密度出現(xiàn)在靠近繞組的夾件局部位置，達700kW/m3。

圖7 額定下變壓器雜散損耗密度分布

3.1.2 GIC作用下時的變壓器損耗

以50A直流的情形為例進行分析，對應(yīng)的變壓器雜散損耗密度分布如圖8所示。由圖可見，損耗密度最大值出現(xiàn)在繞組端部對應(yīng)的夾件區(qū)域，約為3.5MW/m3。由于主柱拉板位于繞組和鐵心主柱之間，直流偏磁時，其損耗密度大幅度增加。變壓器油箱的損耗密度最大值出現(xiàn)在繞組正上方的油箱頂部區(qū)域，最大值約為610kW/m3。

圖8 50A直流下變壓器雜散損耗密度分布

表5列出了在變壓器滿負載時不同GIC條件下各部件總損耗的變化，可知變壓器各部件的總損耗隨GIC的增加而增大。由于自耦變壓器一、二次繞組都有直流的通入，因此其損耗都會有所增加，且在大直流下?lián)p耗的增幅變大，但調(diào)壓繞組相比于高、中壓繞組在直流影響下的損耗增加較小。夾件、油箱等雜散損耗在直流影響下的增長倍數(shù)遠大于繞組損耗，考慮其損耗分布的不均勻性，需要關(guān)注所引起的局部溫升過高的問題。

表5 GIC作用下變壓器各部件總損耗結(jié)果

Tab.5 Total loss of transformer components under GIC

3.2 溫升結(jié)果

分別對變壓器額定負載下和GIC作用下的變壓器熱流場計算結(jié)果進行分析。

3.2.1 額定負載下變壓器溫升結(jié)果

圖9為將損耗結(jié)果導(dǎo)入流熱模型計算得到的變壓器溫度基本穩(wěn)定后的溫度場、流場分布和繞組熱點溫升曲線結(jié)果。由圖可見，由于絕緣油的流動將熱量帶到油箱頂部通過出口流入散熱器，熱油經(jīng)過散熱器冷卻后流入油箱底部，各繞組溫度從底部到頂部逐漸增高，各繞組的熱點位置位于繞組頂部第一層處，由流場分布可知，在散熱器出入口附近流速較高，約為0.06m/s。由溫升曲線可知，各繞組熱點溫度隨時間呈指數(shù)形式上升，在10h后基本達到穩(wěn)定。變壓器繞組熱點溫度計算結(jié)果與設(shè)計值對比見表6，兩者吻合良好，驗證了流熱場模型的合理性。

圖9 變壓器熱流場計算結(jié)果

表6 額定下繞組熱點溫升計算結(jié)果

Tab.6 Hot spot temperature rise of windings under rated conditions （單位: ℃）

3.2.2 GIC作用下變壓器溫升結(jié)果

實際GIC通常表現(xiàn)為數(shù)小時內(nèi)連續(xù)大量低至中等水平的窄脈沖，期間并夾雜著小于數(shù)分鐘的高峰值。根據(jù)IEEE C57.163中的建議，在計算變壓器受GIC影響下的各部件的溫升曲線時，GIC應(yīng)為半連續(xù)直流電流，對于每種直流情況，可采用如圖10所示的周期性動態(tài)直流[20]，以此通過施加動態(tài)損耗計算變壓器各部件溫度分布和溫升曲線。圖中，peak為每種GIC的大小，base=30%peak，其中時間b1= 60min，b2=20min，p=2min，通過此周期性動態(tài)電流可近似地模擬數(shù)小時的GIC事件。

在變壓器高中最小分接滿負載運行時溫度達到穩(wěn)態(tài)時，通入每種模擬的GIC電流，重復(fù)兩個周期后來表示整個GIC事件。以表5中最為嚴(yán)重的peak= 100A時情況為例，計算得出溫度最大時的溫度分布如圖11所示。由圖可知，較大的直流電流通入繞組，繞組整體溫度也升高明顯，由于三相磁路不對稱等因素所引起三相繞組損耗的差異，A、C相繞組熱點溫度最高比B相大5℃左右，繞組熱點位于繞組上端部，其中中壓繞組溫升最高，其熱點溫度為123.6℃，因繞組渦流損耗增加較大，繞組最低溫度位置并不位于最下端。由于主柱拉板損耗密度比較集中，且拉板所處位置散熱能力弱，在100A直流脈沖下，變壓器最高溫度出現(xiàn)在主柱拉板處，高達133℃。非平均熱源的加載使油箱熱點位置也對應(yīng)于其渦流損耗集中的位置，因油箱有較好的散熱條件，其熱點溫升遠低于夾件和繞組。

圖11 Ipeak=100A時變壓器部件溫度分布

繞組熱點和夾具、油箱熱點在peak分別為50A、100A時溫升曲線如圖12所示，由圖可知，在base直流階段，繞組溫度隨時間逐漸上升，而夾具和油箱熱點溫度短時間急劇上升，這是因為在額定情況下，夾具和油箱損耗值較低，溫度不會出現(xiàn)明顯升高，但由于其對直流的存在很敏感，在直流通入后渦流損耗值會成倍地增加。

圖12 繞組、夾具和油箱熱點的動態(tài)溫升曲線

在每一個peak直流脈沖階段，由于階躍損耗的作用，熱點溫度急劇升高，當(dāng)階躍損耗開始下降時，熱點溫度就會急劇降低，在熱點溫度降到一定數(shù)值后其下降速度會有所減緩。相對于環(huán)境溫度20℃，在第四個50A直流脈沖階段，繞組熱點溫升為82.8℃，夾具熱點溫升為73.3℃，油箱熱點溫升為69.9℃，其中100A直流脈沖階段繞組熱點溫升達103.6℃，夾具熱點溫升高達113℃。

由圖12b可知，對于變壓器初始溫度較高時，如在環(huán)境溫度40℃下，此時變壓器繞組熱點溫度將超過IEEE/IEC中的長時GIC事件絕緣材料耐熱水平限值140℃，若直流脈沖高達100A時的GIC持續(xù)時間過久時，將會導(dǎo)致絕緣材料惡化，而夾具熱點短時溫度將達到153℃，接近短時GIC絕緣耐受限值160℃。因此，在GIC脈沖為100A時該變壓器仍有足夠承受能力，但在某些條件下會出現(xiàn)變壓器內(nèi)部局部過熱而接近絕緣耐受限值。相應(yīng)地，在更高GIC水平下該變壓器的安全運行將會受到嚴(yán)重威脅。

4 變壓器GIC溫升防范措施建議

由溫升結(jié)果可見，GIC對變壓器的影響具有短時和溫升劇烈的特點，即使變壓器具有足夠的抗GIC干擾的能力，但變壓器的累積效應(yīng)會加速其老化，甚至發(fā)生故障和損壞，為降低GIC引起變壓器局部過熱，優(yōu)化變壓器設(shè)計，可采用以下防范措施：

（1）因GIC引起變壓器渦流損耗的嚴(yán)重增加是變壓器溫升的主要原因，故可以通過改善變壓器構(gòu)件的結(jié)構(gòu)或材料進行防范，繞組線圈導(dǎo)線采用逐層換位式連接；拉板、夾件等結(jié)構(gòu)件適當(dāng)開槽或采用非導(dǎo)磁性材料；對于高漏磁區(qū)域（如夾件和拉板的連接處）的局部過熱問題，可通過在局部加設(shè)磁屏蔽來降低局部渦流損耗值。

（2）GIC引起變壓器漏磁增加的原因是變壓器鐵心磁通半波飽和，因此對于直流偏磁或GIC高發(fā)地區(qū)，在變壓器制造過程中建議采用高導(dǎo)磁材料鐵心，通過增加鐵心飽和裕度來提高變壓器GIC電流的耐受閾值。若條件允許，也可在變壓器內(nèi)部安裝補償繞組，當(dāng)檢測到GIC流入到變壓器繞組時，使補償繞組產(chǎn)生反向直流磁動勢來抵消GIC產(chǎn)生的磁動勢，從而降低鐵心的飽和程度。

（3）對運行中的變壓器來說，針對GIC主要通過中性點侵入變壓器進而對其造成威脅，可在變壓器中性點采用相關(guān)直流偏磁抑制措施，如在中性點串接小電阻、電容器等，或者在中性點安裝反向補償裝置，根據(jù)檢測到的GIC脈沖電流，自動向變壓器中性點注入與GIC大小相同的反向直流電流。

5 結(jié)論

本文基于電磁-熱流場計算模型，以一臺三相五柱自耦變壓器為研究對象，計算了在不同GIC影響下變壓器的損耗和溫升，主要結(jié)論如下：

1）在滿負載運行工況下計算變壓器各部件的損耗和繞組熱點溫升，與設(shè)計值吻合較好，誤差在±3.5%以內(nèi)，驗證了本分析方法的有效性。

2）變壓器在GIC作用下，夾件、拉板等雜散損耗增加明顯且隨直流偏磁程度增加而增大，局部損耗密度遠大于整體損耗密度，大直流的注入不僅會引起繞組渦流損耗的增加還會導(dǎo)致直流電阻損耗明顯增大，進而引起繞組整體特別是端部熱點溫度升高。

3）在不同GIC事件作用下的變壓器溫升仿真結(jié)果發(fā)現(xiàn)，該變壓器在GIC水平為100A時依舊有足夠的耐受能力但已接近允許限值，其中在主柱拉板的局部溫升最為嚴(yán)重。

[1] 王澤忠, 黃天超. 變壓器GIC-Q動態(tài)關(guān)系分析[J]. 電工技術(shù)學(xué)報, 2021, 36(9): 1948-1955.

Wang Zezhong, Huang Tianchao. Analysis of GIC-Q dynamic relationship of transformer[J]. Transactions of China Electrotechnical Society, 2021, 36(9): 1948- 1955.

[2] 王豐華, 周翔, 高沛, 等. 基于繞組熱分布的改進油浸式變壓器繞組熱點溫度計算模型[J]. 高電壓技術(shù), 2015, 41(3): 895-901.

Wang Fenghua, Zhou Xiang, Gao Pei, et al. Improved thermal circuit model of hot spot temperature in oil-immersed transformers based on heat distribution of winding[J]. High Voltage Engineering, 2015, 41(3): 895-901.

[3] 王澤忠, 李明洋, 宣夢真, 等. 單相四柱式變壓器直流偏磁下的溫升試驗及仿真分析[J]. 電工技術(shù)學(xué)報, 2021, 36(5): 1006-1013.

Wang Zezhong, Li Mingyang, Xuan Mengzhen, et al. Temperature rise test and simulation of single-phase four-column transformer under DC-bias[J]. Transa- ctions of China Electrotechnical Society, 2021, 36(5): 1006-1013.

[4] 李冰, 王澤忠, 劉海波, 等. 直流偏磁下500kV單相變壓器振動噪聲的試驗研究[J]. 電工技術(shù)學(xué)報, 2021, 36(13): 2801-2811.

Li Bing, Wang Zezhong, Liu Haibo, et al. Experiment on vibro-acoustic characteristic of 500kV single- phase transformer under DC-bias[J]. Transactions of China Electrotechnical Society, 2021, 36(13): 2801- 2811.

[5] 王澤忠, 黃天超. 變壓器地磁感應(yīng)電流-無功功率動態(tài)關(guān)系分析[J]. 電工技術(shù)學(xué)報, 2021, 36(9): 1948- 1955.

Wang Zezhong, Huang Tianchao. Analysis of geo- magnetically induction current-reactive power dynamic relationship of transformer[J]. Transactions of China Electrotechnical Society, 2021, 36(9): 1948-1955.

[6] 王澤忠, 司遠, 劉連光. 地磁暴對電力系統(tǒng)穩(wěn)定性的影響[J]. 電工技術(shù)學(xué)報, 2021, 37(7): 1780-1788.

Wang Zezhong, Si Yuan, Liu Lianguang. Influence of geomagnetic storms on the stability of power system[J]. Transactions of China Electrotechnical Society, 2021, 37(7): 1780-1788.

[7] Takasu N, Oshi T. An experimental analysis of DC excitation of transformers by geomagnetically indu- ced currents[J]. IEEE Transactions on Power Delivery, 1994, 9(2): 1173-1182.

[8] 張書琦, 汪可, 李金忠, 等. 單相單柱旁軛變壓器空負載工況下的直流偏磁關(guān)鍵性能測試[J]. 中國電機工程學(xué)報, 2019, 39(14): 4334-4344, 4281.

Zhang Shuqi, Wang Ke, Li Jinzhong, et al. Tests of DC bias key performances of power transformer with single-phase three-limb core under no-load and rated- load conditions[J]. Proceedings of the CSEE, 2019, 39(14): 4334-4344, 4281.

[9] Price P R. Geomagnetically induced current effects on transformers[J]. IEEE Power Engineering Review, 2002, 22(6): 62-62.

[10] Zhang Bing, Liu Lianguang, Liu Yilu, et al. Effect of geomagnetically induced current on the loss of transformer tank[J]. IET Electric Power Applications, 2010, 4(5): 373-379.

[11] Picher P, Bolduc L. Study of the acceptable DC current limit in core-form power transformers[J]. IEEE Transactions on Power Delivery, 1997, 12(1): 257-265.

[12] Torriano F, Campelo H, Picher P, et al. Experimental and numerical thermofluid study of a disc-type transformer winding scale model[C]//XXIIth Inter- national Conference on Electrical Machines, Lausanne, 2016: 2833-2839.

[13] Vakhnina V V, Shapovalov V A, Kuznetsov V N, et al. The influence of geomagnetic storms on thermal processes in the tank of a power transformer[J]. IEEE Transactions on Power Delivery, 2015, 30(4): 1702- 1707.

[14] 韓金華, 楊曉輝, 王豐華, 等. 500kV電力變壓器直流偏磁損耗特性的仿真研究[J]. 高壓電器, 2015, 51(5): 40-46.

Han Jinhua, Yang Xiaohui, Wang Fenghua, et al. Simulation study of DC bias loss characteristics of 500kV power transformer[J]. High Voltage Apparatus, 2015, 51(5): 40-46.

[15] 楊曉輝, 王豐華, 段若晨, 等. 500kV電力變壓器直流偏磁耐受性能的仿真研究[J]. 電力系統(tǒng)保護與控制, 2014, 42(14): 60-66.

Yang Xiaohui, Wang Fenghua, Duan Ruochen, et al. Simulation study of DC bias enduring ability of 500kV power transformer[J]. Power System Pro- tection and Control, 2014, 42(14): 60-66.

[16] Nishiura R, Yamashita S, Kano S. Simulation analysis of geomagnetically-induced currents (GIC) effects on shell-form transformers[C]//IEEE Power & Energy Society General Meeting, Vancouver, BC, 2013: 1-5.

[17] Xia Yunyan, Han Yongsen, Liu Jinhui. Calculation of winding loss and temperature rise of large power transformer under DC magnetic biasing[C]//Inter- national Conference on the Properties and Appli- cations of Dielectric Materials, Xi’an, 2018: 590-593.

[18] 尹忠東, 魏文思. 基于有限元算法的變壓器繞組渦流效應(yīng)分析及損耗計算[J]. 智慧電力, 2018, 46(5): 71-77.

Yin Zhongdong, Wei Wensi. Analysis of eddy current effect and loss calculation of transformer winding based on finite element algorithm[J]. Smart Power, 2018, 46(5): 71-77.

[19] 劉弘景, 張丹丹, 李偉, 等. 大型油浸式變壓器繞組溫度場的有限元分析[J]. 高壓電器, 2019, 55(12): 83-89.

Liu Hongjing, Zhang Dandan, Li Wei, et al. Finite element analysis for temperature field of winding in large oil-immersed transformer[J]. High Voltage Apparatus, 2019, 55(12): 83-89.

[20] IEEE Std. C57.163-2015 IEEE guide for establishing power transformer capability while under geomag- netic disturbances[S]. 2015.

Calculation of Temperature Rise of Large Transformer under Geomagnetically Induced Current

12

（1. College of Electrical Engineering Shanghai University of Electric Power Shanghai 200090 China 2. Department of Electrical Engineering Shanghai JiaoTong University Shanghai 200240 China）

To accurately analyze the feature of temperature rise of power transformer under the action of geomagnetically induced current (GIC), an electromagnetic-flow and heat field coupling 3D simulation model for a 500kV three-phase integrated oil-immersed power transformer is proposed in this paper. The loss density and temperature field distribution of internal components of power transformer under different GICs are calculated when the transformer is operated at rated load. Meanwhile, the dynamic temperature rise curve of internal components of power transformer under different GICs is obtained according to the IEEE C57.163 standard. The simulation results are in good agreement with the designed results. In addition, the eddy current loss distribution of the clamp, pulling plate and oil tank of the power transformer is uneven after being affected by GIC, and increases with the increase of the GIC. With the action of GIS pulse, the local temperature of metal components of the transformer rises rapidly and approaches the permissible limit for the pulling plate of main column at the 100A DC pulse. The results can provide a reference for the tolerance performance analysis of large transformer under the action of GIC.

Power transformer, geomagnetically induced current (GIC), loss density, temperature field, temperature rise

10.19595/j.cnki.1000-6753.tces.201312

TM411

朱 濤 男，1995年生，碩士研究生，研究方向為電力變壓器直流偏磁、變壓器溫度場。E-mail: 1344098230@qq.com

王豐華 女，1973年生，博士，副教授，研究方向為電力變壓器振聲監(jiān)測及裝備智能化、電力變壓器直流偏磁、電力系統(tǒng)接地技術(shù)等。E-mail: fhwang7723@sjtu.edu.cn（通信作者）

2020-09-29

2021-01-30

（編輯 崔文靜）