樹脂澆注干式變壓器三維溫度場仿真計(jì)算

張 牧，高立業(yè)，魏 娟，修春波

（天津工業(yè)大學(xué)電氣工程與自動(dòng)化學(xué)院，天津 300387）

樹脂澆注干式變壓器三維溫度場仿真計(jì)算

張 牧，高立業(yè)，魏 娟，修春波

（天津工業(yè)大學(xué)電氣工程與自動(dòng)化學(xué)院，天津 300387）

為便于診斷干式變壓器的故障，對(duì)干式變壓器進(jìn)行熱點(diǎn)分析.根據(jù)干式變壓器實(shí)際結(jié)構(gòu)建立了三維對(duì)稱熱傳導(dǎo)模型，并基于多物理場有限元軟件平臺(tái)Comsol Multiphysics對(duì)干式變壓器開展溫度場仿真計(jì)算，計(jì)算得出溫度場分布云圖以及鐵芯、繞組溫升曲線.通過比對(duì)實(shí)際數(shù)據(jù)和仿真分析結(jié)果，從而確定溫度分布及最熱點(diǎn)溫度位置在低壓繞組距底部2/3處，可為該類干式變壓器的在線診斷及結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供理論依據(jù).

干式變壓器；流體溫度場；熱點(diǎn)溫度；有限元仿真

樹脂澆注干式變壓器因具備防火、防爆、運(yùn)行可靠、維護(hù)方便、體積小等顯著優(yōu)點(diǎn)，已在城市高層建筑、石油、化工煤礦等場所得到了廣泛的應(yīng)用[1-2].然而，在干式變壓器容量不斷增加的同時(shí)，發(fā)熱和冷卻問題也逐漸的凸顯出來，并成為影響其安全運(yùn)行的主要因素.研究表明，干式變壓器的運(yùn)行不穩(wěn)定或毀壞主要是由于變壓器的繞組溫度超過絕緣材料的耐熱極限[3-4]，其中對(duì)于A級(jí)絕緣的變壓器溫度，每增加8℃，變壓器絕緣老化速度將增加一倍[5].繞組過熱不僅會(huì)造成絕緣的加速老化進(jìn)而影響其壽命，而且會(huì)引發(fā)嚴(yán)重的安全問題，國內(nèi)外每年都會(huì)發(fā)生因?yàn)槔@組過熱造成一系列惡性事故，導(dǎo)致了巨大的經(jīng)濟(jì)損失[6-8].因此對(duì)干式變壓器進(jìn)行熱點(diǎn)分析顯得尤為重要，不僅可以維護(hù)干式變壓器的正常安全運(yùn)行，同時(shí)也為該類干式變壓器的結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供了理論依據(jù).

國內(nèi)外研究者經(jīng)過多年研究，獲得了許多對(duì)干式變壓器生產(chǎn)設(shè)計(jì)具有重要指導(dǎo)意義的研究成果，大部分是關(guān)于干式變壓器穩(wěn)態(tài)運(yùn)行狀態(tài)下的溫度場分布特性的研究[8-11].在數(shù)值仿真計(jì)算方法中，國內(nèi)的許多專家學(xué)者也做了大量的嘗試，在二維溫度場模擬方面取得了一定成果[12-13]，但由于變壓器的熱問題涉及傳熱學(xué)、流體力學(xué)、電磁學(xué)等多學(xué)科知識(shí)，其復(fù)雜性和試驗(yàn)難度限制了對(duì)其全面的研究，尤其是三維溫度場立體建模目前研究還十分有限，尚未有能夠全面反映出干式變壓器的整體溫度場模型.本文針對(duì)以上問題，根據(jù)干式變壓器實(shí)際結(jié)構(gòu)建立了三維面對(duì)稱熱傳導(dǎo)模型，基于Comsol平臺(tái)對(duì)干式變壓器開展溫度場仿真計(jì)算，并結(jié)合傳熱學(xué)理論對(duì)干式變壓器溫度場分布情況進(jìn)行理論分析，找出干式變壓器運(yùn)行中的最熱點(diǎn)溫度的位置，為熱點(diǎn)溫度傳感器及監(jiān)測單元的安裝點(diǎn)提供數(shù)據(jù)支撐.

1 溫升計(jì)算

干式變壓器運(yùn)行時(shí)鐵芯由于磁滯損耗、渦流損耗及繞組由于電阻和引線損耗等產(chǎn)生熱量.其中一部分熱量儲(chǔ)存在干變的發(fā)熱體中，使發(fā)熱體本身溫度升高，另一部分熱量發(fā)散到周圍介質(zhì)中，使周圍介質(zhì)溫度升高[14-15].熱量向周圍傳遞的方式可分為熱傳導(dǎo)、熱對(duì)流和熱輻射3種基本形式[16].而干式變壓器內(nèi)部的熱量傳遞主要是通過熱傳導(dǎo)來實(shí)現(xiàn)的，同時(shí)在設(shè)備表面及外界環(huán)境間還存在著熱對(duì)流和熱輻射.干式變壓器的溫升主要包括鐵芯溫升和繞組溫升.

1.1 鐵芯溫升計(jì)算

由于鐵芯材料和結(jié)構(gòu)的不同，其經(jīng)驗(yàn)系數(shù)取值也不同.由硅鋼片疊裝的鐵芯，計(jì)算其平均溫升時(shí)，經(jīng)驗(yàn)系數(shù)k、n分別取0.36和0.8[17]，鐵芯溫升按如下經(jīng)驗(yàn)公式（1）計(jì)算：

式中：τc為鐵芯溫升（℃）；qc為鐵芯散熱面單位熱負(fù)荷（W/m2）；Pc為鐵芯損耗（W）；Sc為鐵芯的有效散熱面積（m2）.

鐵芯損耗Pc主要包括磁滯損耗Ph和渦流損耗Pe，由方程組（2）表示為：

式中：η為鐵芯材料的磁滯系數(shù)；f為頻率（Hz）；Bmax為最大磁通密度（T）；d為硅鋼片厚度（cm）；V為鐵芯體積（cm3）.η、Bmax、d、V都是材料的固有屬性，可以通過變壓器出廠手冊(cè)進(jìn)行查詢，屬于已知量，因此，只需測出頻率就可以計(jì)算出鐵芯的損耗.

1.2 繞組溫升計(jì)算

低壓和高壓繞組的發(fā)熱主要是由各自的負(fù)載損耗引起.負(fù)載損耗由基本損耗和附加損耗構(gòu)成，其中基本損耗是指繞組導(dǎo)線的直流電阻損耗，是導(dǎo)致繞組溫度升高的主要原因.因此，本仿真中負(fù)載損耗只考慮直流電阻損耗[18]，以簡化計(jì)算過程，提高運(yùn)算速度.

繞組溫升計(jì)算公式由公式（3）給出：

式中：τw為繞組溫升（℃）；qw為繞組散熱面單位熱負(fù)荷（W/m2）；Pw為繞組損耗（W）；Sw為繞組有效散熱面積（m2）.n，k分別為經(jīng)驗(yàn)系數(shù)，經(jīng)驗(yàn)系數(shù)的取值由其散熱條件確定.

本文所用干式變壓器為樹脂澆注型，低壓繞組表面與氣道直接接觸，k取0.66，高壓繞組采用分段層式，繞組中設(shè)有豎直氣道，有2個(gè)軸向散熱面，k、n分別取0.46與0.8.

繞組損耗計(jì)算公式由公式（4）給出：

式中：Pk為變壓器短路有功功率（kW）；β為負(fù)載系數(shù)；Sjs為變壓器計(jì)算負(fù)荷（kVA）；Se為變壓器額定容量（kVA）.Pk、Se都是變壓器固有參數(shù)，可以通過變壓器出廠手冊(cè)查詢，屬于已知量，因此，只要獲得變壓器的計(jì)算負(fù)荷便可以計(jì)算出繞組損耗.

2 溫度場的邊界條件

在仿真計(jì)算中，模型合理的物理過程與邊界條件是保證計(jì)算穩(wěn)定收斂、結(jié)果正確的基礎(chǔ).為了使導(dǎo)熱微分方程產(chǎn)生唯一解，則需要對(duì)微分方程設(shè)置符合實(shí)際情況的邊界條件和初始條件.邊界條件指導(dǎo)熱體的表面與周圍介質(zhì)間的熱量交換狀況.在進(jìn)行溫度場分析時(shí)，第三類邊界條件最為常見.本文所研究的干式變壓器三維流體溫度場具有熱對(duì)流和熱輻射過程，是典型的第三類邊界條件，對(duì)于干變繞組間氣道的輻射與對(duì)流相比，可認(rèn)為輻射小到可以忽略[18]，表示如式（5）所示：

式中：h為對(duì)流換熱系數(shù)（W/（m2·K））；Tw為物體表面溫度（K）；Tf為流體溫度（K）.

3 溫度場仿真分析

本仿真對(duì)象為一臺(tái)實(shí)際干式變壓器，型號(hào)為SCB35-1000/35，電壓等級(jí)為35±3×2.5%/10.5 kV，絕緣等級(jí)為F級(jí).此干式變壓器的結(jié)構(gòu)為低壓繞組置內(nèi)，采用銅箔繞制而成，低壓繞組和鐵芯、高壓繞組之間分別設(shè)有氣道進(jìn)行散熱.高壓繞組采用扁銅線繞制成分段多層圓筒式，共分為4段，中間澆注絕緣樹脂，高、低壓繞組兩端各有一定厚度的絕緣端封.干式變壓器三維結(jié)構(gòu)示意圖及三維有限元模型分別如圖1和圖2所示.

圖1 干式變壓器三維結(jié)構(gòu)圖Fig.1 Three-dimensional structure of dry-type transformer

圖2 干式變壓器三維有限元模型Fig.2 Three dimensional finite element model of dry type transformer

為了盡量還原干式變壓器的實(shí)際結(jié)構(gòu)并兼顧仿真計(jì)算的快捷性，在建立仿真模型時(shí)，以干式變壓器中軸平面為對(duì)稱平面建立1/2面對(duì)稱模型，這樣既保證了實(shí)際結(jié)構(gòu)的正確又節(jié)省了仿真計(jì)算時(shí)間，其他結(jié)構(gòu)中的部分簡化為：

（1）變壓器內(nèi)部熱源為鐵芯、低壓繞組、高壓繞組，其截面近似為圓柱體，本仿真將鐵心當(dāng)做圓柱體.

（2）假定變壓器置于一個(gè)氣流方向豎直向上，且流速一定的空間中.

3.1 干式變壓器材料屬性設(shè)定

根據(jù)實(shí)際干式變壓器各部分使用的材料設(shè)置材料的導(dǎo)熱系數(shù)，如表1所示.

表1 干式變壓器材料導(dǎo)熱系數(shù)Tab.1 Material coefficient of thermal conductivity of dry type transformer

3.2 施加載荷

仿真計(jì)算所施加的載荷條件主要包括變壓器運(yùn)行時(shí)各部分產(chǎn)生的損耗以及環(huán)境溫度.仿真數(shù)據(jù)是實(shí)際變壓器廠所提供的1 000 kVA的樹脂澆注式干式變壓器在額定運(yùn)行狀態(tài)下計(jì)算出的變壓器鐵芯、高低壓繞組的損耗.具體載荷數(shù)據(jù)如表2所示.

表2 干式變壓器施加載荷數(shù)據(jù)Tab.2 Dry-type transformer load data

3.3 網(wǎng)格劃分

干式變壓器的流體溫度場模型較為復(fù)雜，在進(jìn)行網(wǎng)格劃分時(shí)需要采用多重網(wǎng)格劃分方法.網(wǎng)格尺寸的大小決定著仿真結(jié)果的精確度.網(wǎng)格劃分越細(xì)計(jì)算結(jié)果越準(zhǔn)確，但對(duì)計(jì)算機(jī)的要求也就越高.實(shí)際處理中，往往采用分區(qū)域不同密度剖分方法以實(shí)現(xiàn)計(jì)算精度與速度的兼顧.為了準(zhǔn)確的反映干式變壓器內(nèi)部溫度場分布情況，根據(jù)變壓器的結(jié)構(gòu)，依據(jù)重點(diǎn)部位網(wǎng)格細(xì)劃的原則，分別對(duì)鐵芯、高、低壓繞組及樹脂絕緣設(shè)定不同的網(wǎng)格劃分尺寸.本仿真共有網(wǎng)格單元109 314個(gè).網(wǎng)格剖分如圖3所示.

圖3 干式變壓器網(wǎng)格劃分圖Fig.3 Mesh of dry type transformer

3.4 仿真結(jié)果

本次仿真在室溫28℃以下模擬了干式變壓器在100%負(fù)荷運(yùn)行情況下，干式變壓器三維流體溫度場的分布特性，最后將干式變壓器溫度場仿真結(jié)果與實(shí)測溫度值進(jìn)行比對(duì)，其整體溫度場分布如圖4所示.

由圖4可見，干式變壓器內(nèi)部溫度場分布規(guī)律為從底部到頂部溫度變化基本是由低到高再到低，內(nèi)側(cè)溫度高于外側(cè)，最熱點(diǎn)位于干式變壓器低壓繞組距底端2/3處附近，其值為81.2℃，低于F級(jí)絕緣極限溫升100℃.取變壓器的B相作為研究對(duì)象，分別繪制B相鐵芯、低壓和高壓繞組的溫度分布曲線，如圖5所示.

由圖5可知，低壓繞組相對(duì)于鐵芯和高壓繞組溫升較高，是因?yàn)榈蛪豪@組處于鐵芯和高壓繞組之間，其散熱方式主要通過氣道內(nèi)空氣對(duì)流將熱量帶走，而高壓繞組處于外側(cè)不僅可以通過對(duì)流散熱，還可以以輻射散熱，散熱效果明顯優(yōu)于低壓繞組.由圖5（a）可知，鐵芯底部溫度上升比較緩慢，而后迅速上升，主要原因是外側(cè)的低壓繞組對(duì)其產(chǎn)生的影響.

圖4 溫度場分布圖Fig.4 Temperature field distribution

圖5 干式變壓器B相鐵芯、低壓和高壓繞組溫度分布曲線Fig.5 B iron core dry-type transformers，low voltage and high voltage winding temperature distribution curve

由圖5（c）可知，高低壓繞組上半部的散熱效果比下半部散熱效果差，在分段絕緣的地方溫度稍微下降，但繞組溫度總體沿軸向升高，這是由于絕緣的地方?jīng)]有內(nèi)部生熱，繞組內(nèi)部溫度不能持續(xù)上升，所以溫度略有下降.同時(shí)由于空氣在高低壓繞組之間形成的氣道中流動(dòng)，空氣受熱后溫度升高，體積膨脹，會(huì)沿著氣道向上自由流動(dòng).

最熱點(diǎn)部位絕緣性能的好壞直接影響干式變壓器的安全運(yùn)行情況.仿真結(jié)果表明，干式變壓器的最熱點(diǎn)部位為內(nèi)部絕緣弱的部位，因此，在設(shè)計(jì)干式變壓器時(shí)應(yīng)加強(qiáng)此處的絕緣性能，或加強(qiáng)此處的通風(fēng)效果.同時(shí)，在監(jiān)測干式變壓器時(shí)，可以對(duì)此處溫度進(jìn)行監(jiān)測.

4 仿真結(jié)果驗(yàn)證

通過利用干式變壓器中預(yù)埋好的溫度傳感器采集干式變壓器在運(yùn)行過程中的關(guān)鍵點(diǎn)溫度信息，并在實(shí)際結(jié)構(gòu)建立的仿真模型中同樣的部位設(shè)定溫度檢測點(diǎn).通過將仿真計(jì)算所得的溫升值與實(shí)際運(yùn)行過程所采集的溫度信息進(jìn)行比較，以此來驗(yàn)證仿真結(jié)果是否能滿足實(shí)際工程的需要.試驗(yàn)熱電偶布置在B相的情況如圖6所示.

圖6 熱電偶布置情況Fig.6 Thermocouple arrangement

熱電偶布置采用對(duì)稱結(jié)構(gòu)設(shè)定，這樣采集可以最大程度的反映出干式變壓器的溫度信息.鐵芯、低壓和高壓繞組之間都留有空氣氣道，分別將熱電偶按干式變壓器縱向高度三等分放置.仿真計(jì)算數(shù)值與實(shí)測數(shù)值比較結(jié)果如表3所示.通過比較結(jié)果可知，仿真計(jì)算所得數(shù)值與溫度傳感器所測數(shù)值基本一致，相對(duì)誤差表示為：

式中：δ為相對(duì)誤差；x*為計(jì)算值；x為實(shí)測值.經(jīng)計(jì)算相對(duì)誤差在3%以內(nèi)，滿足工程應(yīng)用的要求.

表3 計(jì)算值與實(shí)測值的比較Tab.3 Comparison of calculated and measured values

5 結(jié)論

本文針對(duì)二維仿真不能全面反映干式變壓器熱點(diǎn)溫度場的問題，建立了三維軸對(duì)稱熱傳導(dǎo)模型，基于Comsol平臺(tái)對(duì)干式變壓器開展溫度場仿真計(jì)算.仿真結(jié)果表明，變壓器從底部到頂部溫度變化基本是由低到高再到低，內(nèi)側(cè)溫度高于外側(cè)，最熱點(diǎn)位于干式變壓器低壓繞組距底端2/3處附近，最熱點(diǎn)溫度約為81.2℃.結(jié)合傳熱學(xué)與變壓器實(shí)際結(jié)構(gòu)分析，給出了干式變壓器運(yùn)行中的最熱點(diǎn)溫度的位置，為熱點(diǎn)溫度傳感器及監(jiān)測單元的安裝點(diǎn)提供數(shù)據(jù)支撐，最后本文將仿真計(jì)算的結(jié)果與實(shí)測值比較，仿真計(jì)算所得數(shù)值與溫度傳感器所測數(shù)值基本一致，并且誤差均在3%以內(nèi)，滿足工程應(yīng)用的要求.本研究成果可為生產(chǎn)設(shè)計(jì)以及運(yùn)行提供理論基礎(chǔ)與數(shù)值參考.

[1] 姚志松，姚磊.新型配電變壓器結(jié)構(gòu)、原理和應(yīng)用[M].北京：機(jī)械工業(yè)出版社，2006.

[2]袁臣虎，張牧，高圣偉，等.基于DSP+ARM的干式變壓器狀態(tài)監(jiān)測系統(tǒng)硬件設(shè)計(jì)[J].天津工業(yè)大學(xué)學(xué)報(bào)，2011，30（4）：77-80.

[3]李春華，梁作德，蘇寶國，等.干式變壓器溫度場的有限元仿真計(jì)算[J].變壓器，2012，49（10）：7-10.

[4]姚育成，蔡定國，徐業(yè)彬，等.干式變壓器溫度場數(shù)值模擬[J].變壓器，2012，49（8）：2012，49（8）：31-33.

[5] 陳洪波，李旻.基于傳感器的輸變電設(shè)備狀態(tài)監(jiān)測技術(shù)[M].西安：西安交通大學(xué)出版社，2014.

[6]熊蘭，趙艷龍，楊子康，等.樹脂澆注干式變壓器溫升分析與計(jì)算[J].高電壓技術(shù)，2013，39（2）：265-271.

[7]ESLAMIAN M，VAHIDI B，ESLAMIAN A.Thermal analysis of cast-resindry-typetransformers[J].Energy Conversion and Management，2011，52（3）：2479-2488.

[8]EBRAHIM R，DAVOOD A.Analysis of temperature distribution in cast-resin dry-type transformers[J].Electrical Engineering，2007，89（5）：311-319.

[9] 陳世省，陳秀菊.樹脂澆注干式變壓器繞組最熱點(diǎn)位置的分析[J].變壓器，2002，39（9）：8-11.

[10]范春橋，張重遠(yuǎn).10 kV干式變壓器溫度場仿真計(jì)算與壽命評(píng)估[D].保定：華北電力大學(xué)，2012.

[11]劉小波，姚芳.干式變壓器漏磁場及溫度場的分析與研究[D].天津：河北工業(yè)大學(xué).2012.

[12]趙艷龍，熊蘭.樹脂澆注干式變壓器溫度場數(shù)值計(jì)算及在線監(jiān)測系統(tǒng)研究[D].重慶：重慶大學(xué).2013.

[13]范春橋，張重遠(yuǎn)，李文博，等.干式變壓器溫度場的仿真計(jì)算與分析[J].技術(shù)應(yīng)用，2012，25（3）：76-79.

[14]李巖，劉永志，井永騰.隔爆干式變壓器溫度場有限元分析[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)，2012，24（4）：59-63.

[15]劉永志，李巖，井永騰.基于Fluent軟件的電力變壓器溫度場計(jì)算與分析[J].高壓電器，2012，48（10）：6-11.

[16]曹玉璋.傳熱學(xué)[M].北京：北京航空航天大學(xué)出版社，2001：2-56.

[17]NING Wang，DING Xueqian.Three-dimensional finite element analysis on fluid thermal field of dry-type transformer[C]// 2012 Second International Conference on Instrumentation& Measurement，Computer，Communication and Control，2012.

[18]路長柏，郭振巖，劉文里，等.干式電力變壓器理論與計(jì)算[M].沈陽：遼寧科學(xué)技術(shù)出版社，2003.

Simulation of three-dimensional temperature field of cast-resin dry-type transformers

ZHANG Mu，GAO Li-ye，WEI Juan，XIU Chun-bo
（School of Electrical Engineering and Automation，Tianjin Polytechnic University，TianJin 300387，China）

To conveniently diagnose the faults of dry-type transformer mostly due to abnormal temperature rise directly or indirectly，and the simulation of the hot-spot temperature has important theoretical and engineering application value.According to the actual structure of dry-type transformer，a three-dimensional axisymmetric heat conduction model is established，and the simulation of temperature field calculation of dry-type transformer is carried out based on Comsol Multiphysics to calculate the temperature field distribution and iron core，windings temperature rise curve.Through the comparison between the actual data and the results of the simulation analysis，the temperature distribution and the hottest spot temperature at position 2/3 from the bottom in the low voltage winding is determined.It provides a theoretical reference for dry-type transformer diagnosis and structural optimization of online design.

dry-type transformers；fluid temperature field；hot-spot temperature；finite element simulation

TM412

A

1671－024X（2015）03－0062－05

10.3969/j.issn.1671-024x.2015.03.014

2014-11-27

國家自然科學(xué)基金資助項(xiàng)目（61203302）；天津市基礎(chǔ)與前沿技術(shù)研究計(jì)劃（14JCYBJC18900）

張 牧（1962—），男，教授.研究方向?yàn)闄z測技術(shù)與自動(dòng)化裝置.E-mail：zhangmu@tjpu.edu.cn