陳荊洲, 杜志葉, 王 棟, 任君鵬, 阮江軍
(武漢大學電氣工程學院,湖北 武漢 430072)
高壓開關柜作為電力系統(tǒng)中重要的終端執(zhí)行元件,其安全運行決定著電力系統(tǒng)供電的安全性和可靠性[1-3]。在實際運行中,發(fā)熱問題是影響開關柜安全運行的關鍵因素,近年來由于過熱問題引起開關柜故障的案例越來越多,已經引起了生產運營單位和研究機構的廣泛關注[4-8]。高壓開關柜主要采用封閉式結構,散熱性能較差,當工作在高電壓、大電流條件下就會產生溫升,嚴重時將會嚴重影響開關柜的運行狀態(tài)[9]。目前,國內外的工作主要集中在開關柜在線測溫裝置的研發(fā)和改進上,應用較廣泛的測溫方法包括傳統(tǒng)的溫度傳感器測溫法、紅外測溫法和光纖光柵測溫法[10-15]。
為進一步探究開關柜溫度分布特點,合理選擇溫度監(jiān)測傳感器布置位置,提高過熱監(jiān)測的準確性,需要獲得開關柜內溫度的特征。目前在開關柜內溫度分布計算方面,主要有平均溫升簡易算法、熱路方法和溫度場數值模擬算法三種。其中平均溫升簡易算法[16]是一種較粗略的估算方法,該方法形式簡單且易于使用,但是誤差略大于其余兩種方法,在實際應用中有一定的指導價值;熱路方法[17,18]利用熱電類比法建立開關柜某部件的熱路模型然后求解溫度場,但是該方法求解的精度不是很高,需要進一步修正完善;隨著計算機技術的飛速發(fā)展,通過多物理場耦合計算的數值模擬算法逐漸成為主流[19-21]。數值模擬算法雖然一定程度上提高了開關柜溫度場求解的準確度,但目前的研究分析大多對開關柜模型作了極大的簡化或者單獨考慮某個部件的溫度場,沒有對開關柜的整體模型作全面的分析計算,同時也未考慮開關柜的風機作用,會影響開關柜溫度場的求解精確度。
本文以10kV KYN-28A高壓開關柜為研究對象,通過Solidworks建立開關柜的實體模型,利用ANSYS Icepak對開關柜的溫度場、流體場進行數值計算和分析,得到了開關柜整體的溫度分布特征,并且將仿真結果與溫升試驗數據進行對比,驗證了計算方法的準確性。
開關柜內的熱量主要來自于回路載流導體的焦耳產熱,導體產生的焦耳熱主要以熱傳導的方式在固體間傳遞,同時通過熱對流和輻射換熱傳遞到柜體內空氣中,因此數值計算需要同時求解溫度場方程和流體場方程。
開關柜內的熱量傳遞包含三種基本方式:熱傳導、熱對流和熱輻射。三種傳熱方式的方程分別如下所示:
(1)
Q2=hcA(tw-tf)
(2)
Q3=δ0Aεxt(T14-T24)
(3)
式中,Q1、Q2、Q3分別為熱傳導熱量、熱對流換熱量和熱輻射換熱量;A為換熱面積;λ為導熱系數;hc為對流換熱系數;tw為固體表面溫度;tf為周圍空氣溫度;δ0為斯蒂芬—玻爾茲曼常數,其值為5.669×10-8W/(m2·k4);εxt為系統(tǒng)發(fā)射率;T1為物體表面熱力學溫度;T2為環(huán)境溫度的熱力學溫度。
計算流體動力學一般包含三組控制方程,分別是質量守恒方程、動量守恒方程和能量守恒方程。
質量守恒方程:
(4)
動量守恒方程:
(5)
(6)
(7)
能量守恒方程:
(8)
式中,V為流體流速;u、v、w為流速在x,y,z方向上的分量;ρ為流體密度;μ為流體動力粘度;T為流體溫度;p為流體壓力;Cp為定熱容;Su、Sv、Sw為動量守恒方程廣義源項;ST為粘性耗散項。
本文利用Solidworks軟件針對10kV KYN-28A開關柜建立三維模型,實際開關柜的模型極其復雜,直接利用實物結構模型進行仿真計算非常困難,需要根據研究問題的側重點對實際開關柜模型進行簡化,刪除一些對熱分析影響不大的部件,同時對一些復雜的結構建立簡化的等效模型,建立的模型如圖1所示。
圖1 開關柜三維模型Fig.1 Three dimensional model of switchgear
開關柜模型包含四個獨立的隔離室,分別為母線室、架空進線室、斷路器室和儀表室。斷路器和電流互感器部分均采用簡化模型,同時刪除了外殼及隔板上的安裝孔。在Workbeach平臺下,Icepak的標準CAD接口為DesignMolder(DM),需要把模型導入到DM進行修復處理,將所有的部件轉換為Icepak認可的幾何體后,由DM導入到Icepak。開關柜實際模型中有6個風機,其中5個位于頂部外殼,1個位于斷路器底部,Icepak具有獨立的風機模塊,模型導入后需要在Icepak里添加風機模型。
開關柜發(fā)熱的熱源主要包括載流導體和電接觸的焦耳損耗。文獻[22]表明電接觸部分的產熱對于開關柜溫升有重要影響,本文電接觸部分的發(fā)熱量將單獨計算。整個載流回路發(fā)熱功率為:
P=I2R
(9)
式中,I為流過載流導體的電流;R為載流導體電阻。載流導體的電阻由定義式計算,電接觸部分的接觸電阻受接觸表面狀況、接觸壓力等諸多因素的影響,文獻[23]給出了接觸電阻計算的經驗公式:
Rc=Kc/(Fk/9.8)m
(10)
式中,Rc為接觸電阻;Kc為接觸材料系數;Fk為接觸壓力;m為與接觸形式有關的系數。
開關柜內電接觸部分包括母排搭接處、母排與靜觸頭連接處和動觸頭梅花觸指。其中母排搭接處及母排與靜觸頭連接處通過高強螺栓進行連接,均為面接觸,m值取1,同時接觸壓力可通過緊固力矩計算公式獲得;梅花觸指為點接觸,m值取0.5,接觸壓力通過查找開關說明技術參數獲得。經過計算,當開關柜流過額定電流4000A時各部分的熱量計算結果見表1。
表1開關柜各部分發(fā)熱量
Tab.1 Heating power of each part of switchgear
發(fā)熱部位發(fā)熱功率/W母線3309架空進線7911斷路器7223電流互感器117
表1中所示發(fā)熱量均為三相總發(fā)熱量,也包括電接觸的發(fā)熱量。在仿真計算中,母排搭接處的發(fā)熱量換算成體密度加載在母排上,母排與靜觸頭連接處及動觸頭梅花觸指處的電接觸發(fā)熱量加載在接觸面上。
Solidworks建好后的實體模型導入 ANSYS Icepak后進行模型預處理,開始網格剖分,通過精細控制散熱路徑和母排表面的網格尺寸和數量,來準確模擬熱流的傳熱特性和流動特性,剖分后網格數量約為250萬單元。
母線室進線母排、架空進線母排、斷路器靜觸頭、動觸頭、觸臂的材質為銅,絕緣套管、靜觸頭盒材質為環(huán)氧樹脂,柜體外殼為鍍鋅鋼板。銅和鍍鋅鋼板的物理參數參考文獻[24]得到,環(huán)氧樹脂的物理參數參考文獻[25-27]整理得到。各種材質的物理參數見表2。
表2材料物理參數
Tab.2 Physical parameters of materials
材料物理參數熱導率/(W/(M·K))密度/(g/cm3)比熱率/(J/(G·K))電阻率/(10-6Ω·m)銅3928903900179環(huán)氧樹脂027609814-鍍鋅鋼板4678050137
對各個部件定義物理參數后,按表1計算出的發(fā)熱功率進行熱源加載。環(huán)境溫度設為25℃,根據軟件自動計算出的雷諾數和瑞利數選擇湍流模型,外殼設置對流換熱系數,采用經驗值10W/(m2·K)。輸入實際的風機P-Q曲線,計算模型較復雜,網格數量很多,為了加快求解速度采用雙核并行計算。
3.4.1 溫度場分析
額定工況下開關柜的溫度分布如圖2所示,可以看出,斷路器本體是開關柜中溫度較高的部位,最高溫度出現在中間上動觸頭處,可達61℃,是開關柜溫度監(jiān)控的關鍵部位;母排的溫度在45℃左右,A相母排靠近風扇,散熱面積大,溫度較B相和C相母排較低;母排與靜觸頭接頭由于與斷路器直接接觸,同時由于接觸電阻的影響,出現了較高溫升,其余載流回路的溫度在45℃-50℃之間。具體溫度分布見表3。
圖2 開關柜三維溫度分布圖Fig.2 Three dimensional temperature distribution of switchgear
部位名稱平均溫度/℃部位名稱平均溫度/℃上動觸頭565上靜觸頭560下動觸頭559下靜觸頭556母線462架空進線425電流互感器485出線連接處496
3.4.2 流體場分析
開關柜三維流場中三個不同截面的流速矢量圖如圖3所示。其中,圖3(a)為z軸中心的xy平面圖;圖3(b)和3(c)分別為x=0.55和x=1.45時的yz平面圖。從流速矢量圖中可以得出流體流動的方向和速度大小。
圖3 開關柜截面速度矢量圖Fig.3 Velocity vector diagram of cross section
在柜頂風機和斷路器室底部風機的共同作用下,空氣從左部外殼進風口和前后外殼進風口進入柜體內。 從左板進風口進入的空氣初始平均速度約為5.7m/s;而從前后板進風口進入的空氣初始平均速度約為4.6m/s。從前后板進風口進入的空氣大部分流入母線室和架空進線室,然后從母線室頂部風機流出,母線室和架空進線室風路較通暢,熱量散失較快,因而溫升較低。從左板進風口進入的空氣在斷路器底部風機的作用下大部分進入斷路器室, 然后一部分會從斷路器外側繞經儀表室,最后從斷路器室頂部流出,然而斷路器室結構較為封閉,風路狹窄,而且一部分氣流在儀表室下側形成環(huán)流,導致熱流路徑變長,風速降低;另一小部分則通過上靜觸頭盒流至母線室,最后從母線室頂部流出,然而由于隔板的阻流影響,風速也很低??梢园l(fā)現,斷路器室封閉復雜的結構延長了熱流路徑,降低了風速,極大地降低了該處的散熱能力,導致斷路器部分整體溫升較高。
為了驗證仿真計算結果的準確性,針對本開關柜開展了溫升試驗,溫升試驗中,電流為額定電流,即4000A,環(huán)境溫度為24.8℃。由于斷路器及觸頭部分是開關柜溫度較高的部位,因此測點主要布置在斷路器處,溫度通過溫度傳感器測得。共選取了9個溫度測點,分別為后出線連接處、電流互感器處、進線與靜觸頭連接處、下靜觸頭、下動觸頭、上動觸頭、上靜觸頭、母線與靜觸頭連接處、上分支母線。測點布置如圖4所示。取三相平均溫度作為測點溫度,數字1-9依次代表上述9個測點,仿真計算結果與試驗結果對比見表4,對比如圖5所示。
圖4 開關柜溫度測點布置圖Fig.4 Arrangement of temperature measuring points of switchgear
測點編號平均溫度值/℃實測值計算值相對誤差(%)1474496462462485493529554474532556455555559076579565-2475495602853856244945246222
圖5 溫升計算值與實測值對比圖Fig.5 Experiment result compared with simulation result
從表4可以得出,仿真計算與試驗實測的溫度分布趨勢一致,最高溫度均在上動觸頭處。各個測點溫度值與計算溫度值也能較好地吻合,誤差均在5% 以內,工程上可以接受,驗證了本文計算結果的正確性。
(1)建立了開關柜有限元分析三維模型,計算了開關柜各個部位的發(fā)熱功率,并對熱流路徑和流固交界面的網格進行了精細控制,在此基礎上利用Icepak對開關柜的溫度場和流體場進行了數值分析。計算結果表明上動觸頭是開關柜溫度最高的部位,最高溫升達到36℃,是開關柜溫升控制和溫度監(jiān)控的關鍵。
(2)將仿真計算結果與溫升試驗結果進行對比,計算結果與試驗結果吻合的很好,誤差在5%以內,工程上可以接受,驗證了仿真計算的準確性。
(3)對流體場進行了分析,發(fā)現由于斷路器室結構較為封閉,導致熱流路徑延長,風速降低,極大的降低了斷路器室的散熱能力。在接下來的研究中,應根據開關柜實際情況,優(yōu)化斷路器室的風路結構,減小熱流路徑,來控制該處溫升。
[1] 丁健(Ding Jian). 金屬封閉高壓開關柜觸頭發(fā)熱機理分析及預防過熱故障措施探討(Heating mechanism analysis on metal-enclosed high voltage switchgear contact with measures to prevent overheating)[J]. 高壓電器(High Voltage Apparatus),2012,48(8):114-122.
[2] 朱根良(Zhu Genliang). 淺議中壓開關柜事故調查中的故障分析(Simple explanation of faults analysis in incident investigation of medlum voltage switchgear)[J]. 高壓電器(High Voltage Apparatus),2002,38(5):62-63.
[3] 王有元,李寅偉,陸國俊,等(Wang Youyuan,Li Yinwei,Lu Guojun,et al.). 開關柜局部放電暫態(tài)對地電壓傳播特性的仿真分析(Simulation of transient earth voltage aroused by partial discharge in switchgears)[J]. 高電壓技術(High Voltage Engineering),2011,37(7):1683-1688.
[4] 李晶,滿家健,周電波,等(Li Jing,Man Jiajian,Zhou Dianbo,et al.). 金屬封閉高壓開關柜手車觸頭溫度場仿真分析及實驗研究(Simulation analysis and experimental research on temperature distribution of metal-enclosed high voltage switchgear handcart contactor)[J]. 高壓電器(High Voltage Apparatus),2015,51(8):74-79.
[5] 鞏憲鋒,衣紅鋼,王長松,等(Gong Xianfeng,Yi Honggang,Wang Changsong,et al.). 高壓開關柜隔離觸頭溫度監(jiān)測研究(Research on temperature monitoring of isolators in HV switchgear)[J]. 中國電機工程學報(Proceedings of the CSEE),2006,26(1):155-158.
[6] 周巖(Zhou Yan). 高壓開關柜觸頭溫度場的數值仿真分析(Numerical simulation of high voltage switchgear panel contact temperature)[J]. 傳感器與微系統(tǒng)(Transducer and Microsystem Technologies),2008,27(12):102-105.
[7] 宋帆,申春紅,徐建源,等(Song Fan,Shen Chunhong,Xu Jianyuan,et al.). 800kV GIS隔離開關磁場—溫度場計算與分析(Calculation and analysis on magneto- thermal fields of 800 kV GIS disconnector)[J]. 高電壓技術(High Voltage Engineering),2008,34(7):1383-1388.
[8] 熊蘭,徐敏捷,楊子康,等(Xiong Lan,Xu Mingjie,Yang Zikang,et al.). 高壓開關柜電纜室溫度場分析及在線監(jiān)測系統(tǒng)構建(Analysis of temperature field in cable chamber of high voltage switchgear and construction of on line monitoring system)[J]. 電力自動化設備(Electric Power Automation Equipment),2014,34(6):153-157.
[9] 王新彤(Wang Xintong). 基于光纖光柵的開關柜溫度在線監(jiān)測系統(tǒng)的開發(fā)與應用(The development and application of switchgear temperature monitoring system based on fiber bragg grating)[D]. 北京:華北電力大學(Beijing:North China Electric Power University),2010.
[10] 李成榕,馬國明(Li Chengrong,Ma Guoming). 光纖布喇格光柵傳感器應用于電氣設備監(jiān)測的研究進展(Research and development of fiber bragg grating monitoring in electrical engineering)[J]. 中國電機工程學報(Proceedings of the CSEE),2013,33(12):114-122.
[11] 張艷,田競,葉逢春,等(Zhang Yan,Tian Jing,Ye Fengchun,et al.). 基于紅外傳感器的高壓開關柜溫度實時監(jiān)測網絡的研制(Development of temperature real-time monitoring network for high voltage switchgear based on IR sensor)[J]. 高壓電器(High Voltage Apparatus),2005,41(2):91-94.
[12] 鄔鋼(Wu Gang). 光纖在開關柜觸頭溫度監(jiān)測中的應用(Application of optical fiber in temperature monitoring of the contacts of switchgear)[J]. 高電壓技術(High Voltage Engineering),2006,32(2):122.
[13] 趙炳成(Zhao Bingcheng). 高壓開關柜在線監(jiān)測系統(tǒng)的開發(fā)與應用(Development and application of on-line monitoring of high voltage switchgear)[D]. 杭州:浙江大學(Hangzhou:Zhejiang University),2011.
[14] 田毅,黃新波(Tian Yi,Huang Xinbo). 高壓開關柜溫升在線監(jiān)測系統(tǒng)(On-line monitoring system of temperature inside HV switchgear cabinet)[J]. 高壓電器(High Voltage Apparatus),2010,46(3):64-67,71.
[15] 張純璽(Zhang Chunxi). 光纖光柵溫度監(jiān)測系統(tǒng)設計及其在高低壓開關柜中的應用(Design offiber Bragg grating temperature monitoring system and its application in high-voltage & low-voltage switchgear)[D]. 濟南:山東大學(Jinan:Shandong University),2010.
[16] 黃恒禧(Huang Hengxi). 應用高壓開關柜溫升的簡易算法及降低相應發(fā)熱的技術措施(A study of simple arithmetic for applied high voltage switchboard temperature rise and the technical easures to reduce relevant heat)[J]. 福建建筑(Fu jian Architecture & Construction),2008,117(3):91-93.
[17] 趙健康,樊友兵,王曉兵,等(Zhao Jiankang, Fan Youbing, Wang Xiaobing,et al.). 高壓電力電纜金屬護套下的熱阻特性分析(Thermal resistance properties of the part between metal sheath and conductor in high-voltage power cable)[J]. 高電壓技術(High Voltage Engineering),2008,34(11):2483-2487.
[18] 孫凱祺,劉浩,傅晨釗,等(Sun Kaiqi, Liu Hao,Fu Chenzhao,et al.). 220kV 瓷柱式斷路器修正熱路模型研究(Modified thermal circuit model for 220kV porcelain stanchion-type circuit breaker)[J]. 華東電力(East China Electric Power),2013,41(6):1215-1219.
[19] 廖才波,阮江軍,劉超,等(Liao Caibo, Ruan Jiangjun,Liu Chao,et al.). 油浸式變壓器三維電磁-流體-溫度場耦合分析方法(3-D coupled electromagnetic-fluid-thermal analysis of oil-immersed triangular wound core transformer)[J]. 電力自動化設備(Electric Power Automation Equipment),2015,35(9):150-155.
[20] 溫志偉,傅德平,顧國彪(Wen Zhiwei,Fu Deping,Gu Guobiao). 浸潤式混合強迫內冷的蒸發(fā)冷卻汽輪發(fā)電機定子三維溫度場的仿真分析(Simulation of 3D temperature distribution of immersion evaporative cooling turbo-generator combined with forced inner cooling system)[J]. 電工電能新技術(Advanced Technology of Electrical Engineering and Energy),2006,25(3):13-17.
[21] 鄒斌(Zou Bin). 基于有限元法的配電柜溫度場研究(Research on temperature field of switchboard based on finite element method)[D]. 蘇州:蘇州大學(Shuzhou:Shuzhou University),2009.
[22] 王秉政,江健武,趙靈,等(Wang Bingzheng,Jiang Jianwu,Zhao Ling,et al.). 高壓開關柜接觸發(fā)熱溫度場數值計算(Temperature field simulation of contact heating in high voltage switchgear)[J]. 高壓電器(High Voltage Apparatus),2013,49(12):42-48.
[23] 郭鳳儀(Guo Fengyi). 電接觸理論及其應用技術(Electrical contact theory and applied technology)[M]. 北京:北京電力出版社(Beijing:China Power Press),2008.
[24] 徐國政(Xu Guozheng). 高壓斷路器原理與應用(Principle and application of high voltage circuit break-er)[M]. 北京:清華大學出版社(Beijing:Tsinghua University Press),2010.
[25] 程鵬,楊永明,莊曉蕓,等(Cheng Peng,Yang Yongming,Zhuang Xiaoyun,et al.). 立體卷鐵心干式變壓器三維溫度場計算與分析(Calculation and analysis of three-dimensional temperature field of Dry-Type transformer with tridimensional wound core)[J]. 變壓器(Transformer),2013,50(9):12-15.
[26] 蔡定國,姚育成,徐業(yè)彬,等(Cai Dingguo,Yao Yucheng,Xu Yebin,et al.). 基于流固耦合分析的干式變壓器溫度場數值分析(Numerical simulation of dry type transformer temperature field based on the fluid solid coupling analysis)[J]. 華北電力大學學報(Journal of North China Electric Power University),2012,39(2):34-38.
[27] 范春橋(Fan Chunqiao). 10kV 干式變壓器溫度場仿真計算與壽命評估(The simulation of temperature field and life assessment of 10kV Dry-Type transformer)[D]. 北京:華北電力大學(Beijing:North China Electric Power University),2012.