架空導(dǎo)線融冰過程分析

毛天，金梓翔，陳雙濤*

（1-云南電力試驗(yàn)研究院集團(tuán)（有限公司）電力研究院，云南昆明 650217；2-西安交通大學(xué)能源與動力工程學(xué)院，陜西西安 710049）

架空導(dǎo)線融冰過程分析

毛天1，金梓翔2，陳雙濤*2

（1-云南電力試驗(yàn)研究院集團(tuán)（有限公司）電力研究院，云南昆明 650217；2-西安交通大學(xué)能源與動力工程學(xué)院，陜西西安 710049）

覆冰對我國南方電網(wǎng)輸變電設(shè)施造成了相當(dāng)嚴(yán)重的影響，本文綜合考慮了外部環(huán)境與冰層的對流和輻射換熱、重力以及導(dǎo)線表面的水滴流失等因素的影響，建立了短路電流融冰的CFD數(shù)學(xué)物理模型。對不同熱流密度下冰層的融化情況進(jìn)行了模擬。研究結(jié)果表明，熱流密度的線性增大并未造成融冰單位耗時和融冰總耗時的線性減少；在融冰過程最后階段，單位耗時大于其他階段，這種趨勢在小熱流密度融冰過程中尤為明顯。

電網(wǎng)；融冰；CFD模擬；熱流密度

0 引言

近年來，我國南方電網(wǎng)的輸變電設(shè)施時有出現(xiàn)覆冰情況，其危害包括不同期脫冰或不均勻覆冰事故、導(dǎo)線覆冰舞動事故等[1]。架空導(dǎo)線由于覆冰會不同程度的影響輸電線路的可靠性，嚴(yán)重時甚至?xí)绊懙秸麄€電力系統(tǒng)的運(yùn)行。研究融冰方法對減輕冬季輸電線路所受冰雪災(zāi)害意義重大。

常見的除冰方案有機(jī)械除冰法、熱力融冰法、自然脫冰法等。文獻(xiàn)[2]表明，除上述除冰方法外，還有利用電磁脈沖、氣動脈沖、電暈放電、電子凍結(jié)、碰撞前顆粒加熱和凍結(jié)等除冰方法，這些方法大多處在試驗(yàn)階段。而交、直流融冰法是面向大型電網(wǎng)最為可行的兩種融冰方式。劉文濤等[3]基于直流融冰提出了電網(wǎng)大面積冰災(zāi)防御策略。敬華兵等[4]在直流融冰基礎(chǔ)上提出了一種輸電線路柔性直流融冰技術(shù)。短路融冰是利用焦耳效應(yīng)加熱導(dǎo)線使之融冰，用較低電壓提供較大短路電流加熱導(dǎo)線的方法使導(dǎo)線覆冰融化。喬海洋等[5]對短路電流融冰技術(shù)進(jìn)行了概括。短路電流融冰技術(shù)包括不帶負(fù)荷的三相短路融冰技術(shù)、導(dǎo)線對導(dǎo)線兩相短路融冰技術(shù)、導(dǎo)線對地線單項(xiàng)短路融冰技術(shù)、利用短路電磁力除冰技術(shù)等。2008年冰災(zāi)發(fā)展初期，湖南省先后對40條線路采用了交流短路融冰。

由于CFD模型被通常被應(yīng)用在實(shí)驗(yàn)測量難以獲得結(jié)果的場合來揭示物理現(xiàn)象的本質(zhì)，例如潘良高等[6]利用數(shù)值方法研究了微通道內(nèi)氣液兩相流型，而范晨等[7]則利用數(shù)值方法對空氣源熱泵結(jié)霜及除霜過程進(jìn)行了詳細(xì)的分析。本文從融冰的傳熱學(xué)機(jī)理出發(fā)，綜合考慮環(huán)境與冰層的對流輻射換熱影響、重力影響以及導(dǎo)線表面的水滴流失影響的情況，建立了電線短路融冰的冰塊融化過程CFD模型，對不同通電電流下具有不同熱流密度的導(dǎo)線融冰過程進(jìn)行了數(shù)值描述，獲得了不同除冰熱流密度下的融冰時間。

1 導(dǎo)線融冰的數(shù)值模型

1.1 融冰數(shù)學(xué)物理模型的建立

文獻(xiàn)[8]指出，采用短路電流融冰的方法，將冰層從內(nèi)到外進(jìn)行融化，是復(fù)雜的熱力學(xué)過程，必須綜合考慮包括對流、輻射、熱傳導(dǎo)在內(nèi)的諸多因素。本文選用的架空導(dǎo)線型號為LGJ-400/35，導(dǎo)線半徑r = 12 mm，電阻R = 0.06 Ω/km。假定融冰前冰層是與導(dǎo)線同心的均勻分布的圓環(huán)，并取冰層厚度為di= 10 mm。由于冰層在通電流之前的溫度與環(huán)境溫度保持一致，因此取冰層溫度與環(huán)境溫度均為Ta= 268.15 K。導(dǎo)線通電后，產(chǎn)生的焦耳熱先使導(dǎo)線升溫到0 ℃以上，然后融掉導(dǎo)線與冰層接觸面的覆冰，從而在導(dǎo)線表面形成一層水膜，冰層此時在重力作用下下墜一定的距離，而融掉的水會自動下落并透過冰層微小的空隙滲透到大氣環(huán)境中，最終在導(dǎo)線下部形成一個隨融冰時間增加越來越大的水氣層。為了對下墜融冰過程進(jìn)行詳細(xì)的描述，本文針對冰層每次下移1 mm時，不計水平方向的水膜厚度情況下的融冰過程，建立了融冰數(shù)值模型，如圖1所示。

圖1 融冰求解模型

1.2 熱交換系數(shù)h的確定

用h表示外冰層與環(huán)境的總換熱系數(shù)。根據(jù)文獻(xiàn)[9]，熱交換系數(shù)的確定方法如下：

式中：

h1——外冰層對流換熱系數(shù)；

h2——外冰層輻射換熱系數(shù)；

Ri——冰層厚度的兩倍與導(dǎo)線半徑之和，取32 mm；

λa——空氣的熱傳導(dǎo)率，取0.023 8 W/(m·K)；

ε——冰表面的發(fā)射率，取0.9[10]；

σ——斯蒂芬-波爾茲曼常數(shù)，取5.67×10-8W/m2·K4；

Nu1、Nu2——覆冰導(dǎo)線的自然對流和強(qiáng)制對流的Nusselt數(shù)；

B、b——Gr數(shù)決定的系數(shù)。在覆冰條件下，取B為0.48、b為0.25；

Gr、Pr、Re——Grashof數(shù)、Prandlt數(shù)和Reynolds數(shù)；

g——重力加速度，取9.8 m/s2；

ν——空氣動力粘度，取1.328×10-5m2/s；

μ——空氣動粘度系數(shù)，取1.72×10-5kg/(m·s)；

Ca——空氣比熱容，取1,005 J/kg·K；

ρa(bǔ)——空氣密度，取1.293 kg/m3；

νa——風(fēng)速，按三級平均風(fēng)速，取4.4 m/s。

C、n——Reynolds數(shù)(Re)決定的系數(shù)[6]。C、n值如表1所示。

表1 C、n值與Re的關(guān)系

將上述參數(shù)代入式中可最終得出h為38 W/m2·K。

1.3 臨界電流的確定

使得冰層到達(dá)臨界狀態(tài)（冰水混合物273.15 K）的通電電流大小稱之為臨界電流，用Icr表示，而熱流密度Hf與導(dǎo)線電流I有如下關(guān)系：

臨界熱流密度為：

對應(yīng)的臨界電流強(qiáng)度為：

經(jīng)過計算，取臨界熱流密度為190 W/m2，而模擬臨界熱流密度為199 W/m2。圖2為模擬臨界熱流密度下最后1 mm冰層的固液相分界線的局部放大圖。計算采用參數(shù)比理論計算值多的9 W/m2熱流用于保證最后1 mm冰層溫度場穩(wěn)定。而這從理論還是實(shí)際來說，都是合乎傳熱學(xué)原理和現(xiàn)實(shí)情況的。

2 結(jié)果與分析

由于其幾何上的對稱性，針對1/2對稱體進(jìn)行模擬。基于CFD模擬軟件、采用融化模型，對不同通電電流下冰層每融1 mm進(jìn)行了模擬，并對其總的融冰時間進(jìn)行分析。從圖2可以看到，冰層最外層的溫度非常接近0 ℃但永遠(yuǎn)達(dá)不到融化狀態(tài)。這也證實(shí)了最后1 mm融冰耗時突增的事實(shí)與溫度場的穩(wěn)定需要9 W/m2熱流密度來維持的合理性。

圖2 臨界熱流密度下最后1 mm冰層的固液分界

為了更好地描述其熱力學(xué)特性，利用熱流密度來描述加入熱量的大小，所對應(yīng)的電流強(qiáng)度可通過式(11)進(jìn)行換算。圖3～圖8為不同熱流密度下不同時間的冰層溫度分布圖。從圖中可以看出，隨著熱流密度的增加，融冰耗時在相應(yīng)地減少。不同熱流密度下冰層僅在融化層周邊從內(nèi)到外（即與導(dǎo)線上半部分接觸的冰層）形成一個較大的溫度梯度分布，而非融化區(qū)域的溫度場幾乎還是與環(huán)境溫度一致，即說明了前一冰層的融化對后面冰層的溫度分布的影響幾乎可以忽略不計，這也說明本文所建立的融冰模型的合理性。

圖3 熱流密度1,000 W/m2冰層溫度分布隨時間變化（單位：K）

圖4 熱流密度2,000 W/m2冰層溫度分布隨時間變化（單位：K）

圖5 熱流密度3,000 W/m2冰層溫度分布隨時間變化（單位：K）

圖6 熱流密度4,000 W/m2冰層溫度分布隨時間變化（單位：K）

圖7 熱流密度5,000 W/m2冰層溫度分布隨時間變化（單位：K）

圖8 熱流密度6,000 W/m2冰層溫度分布隨時間變化（單位：K）

圖9為不同熱流密度下的單位厚度融冰耗時。為了更清楚的顯示熱流密度對融冰過程的影響，對比分析不同熱流密度下的融冰耗時。從圖9可以看出，單位融冰耗時總體上呈現(xiàn)減小趨勢，而隨著熱流密度的增大，單位耗時減小的趨勢逐漸變緩?？傮w上看，前9 mm單位厚度融冰耗時基本一致，而最后1 mm融冰耗時突增。這是因?yàn)榍? mm的融冰過程中，熱流密度都是用于冰層的吸熱（此期間，最外層冰與環(huán)境溫度一致，故不產(chǎn)生對流輻射換熱），而最后1 mm融冰過程的熱流密度中有一部分用于對流輻射耗熱（這段融冰時間內(nèi)最外層冰的溫度由環(huán)境溫度上升到0 ℃），從而用于冰層吸熱的有效熱流密度減小了，故融冰耗時突增；隨著熱流密度的增大，用于對流輻射耗熱的熱流密度占總熱流密度的比例在減小，即用于冰層吸熱的有效熱流密度所占的比例在增加，從而導(dǎo)致了這種融冰耗時突增的程度在不斷減小。圖10顯示了不同熱流密度下總?cè)诒臅r。從中可以看出，隨著熱流密度的線性增加，融冰的總耗時呈非線性減小。且隨著熱流密度持續(xù)增大，總耗時減少的趨勢逐漸變緩，總耗時與融冰厚度的線性度變得越來越好。

圖9 不同熱流密度下的單位厚度融冰耗時對比

圖10 不同熱流密度下總?cè)诒臅r

3 結(jié)論

本文在綜合考慮外部環(huán)境與冰層的對流輻射換熱影響、重力影響以及導(dǎo)線表面的水滴流失影響的情況下建立了融冰的物理模型。并基于FLUENT軟件，對不同熱流密度/通電電流強(qiáng)度下冰層每融化1 mm的過程進(jìn)行了模擬，得到了不同熱流密度下的不同時刻融冰的溫度分布云圖，同時對單位厚度融冰耗時和總?cè)诒臅r進(jìn)行了對比分析。結(jié)果表明：

1）融冰的最后，由于與外部環(huán)境的對流換熱增強(qiáng)，導(dǎo)線產(chǎn)生的熱量無法全部用于融冰，一定電流密度下的融冰單位耗時顯示出增加的趨勢，這種趨勢在低熱流密度融冰過程中更加明顯；

2）不同熱流密度下的單位融冰耗時總體上呈減小趨勢，但是隨著熱流密度的增大，單位耗時減小的趨勢逐漸變緩；

3）隨著熱流密度線性增加，融冰的總耗時呈非線性減小，且隨著熱流密度持續(xù)增大，總耗時減少的趨勢逐漸變緩。

[1] 常浩, 石巖, 殷威揚(yáng), 等. 交直流線路融冰技術(shù)研究[J].電網(wǎng)技術(shù), 2008, 32(5): 1-6.

[2] 楊小弟. 架空輸電線覆冰和融冰計算模型研究[D]. 北京: 華北電力大學(xué), 2012.

[3] 劉文濤, 和識之, 陳亦平, 等. 基于直流融冰的電網(wǎng)大面積冰災(zāi)防御策略[J]. 電力系統(tǒng)自動化, 2012, 36(11): 102-106.

[4] 敬華兵, 年曉紅. 輸電線路柔性直流融冰技術(shù)[J]. 高電壓技術(shù), 2012, 38(11): 3060-3066.

[5] 喬海洋, 歐陽明, 何友光, 等. 采用電阻發(fā)熱的帶電融冰技術(shù)探討[J]. 電氣時空, 2009(12): 17-18.

[6] 潘良高, 徐琛, 柏祥華, 等. 微通道內(nèi)氣液兩相流型的數(shù)值模擬[J]. 制冷技術(shù), 2014, 34(4): 8-12.

[7] 范晨, 梁彩華, 江楚遙, 等. 空氣源熱泵結(jié)霜/除霜特性的數(shù)值模擬[J]. 制冷技術(shù), 2014, 34(1): 18-25.

[8] 楊世銘, 陶文銓. 傳熱學(xué)[M]. 3版. 北京: 高等教育出版社, 1998.

[9] 范松海. 輸電線路短路電流融冰過程與模型研究[D].重慶: 重慶大學(xué), 2010.

[10] 劉和云. 架空導(dǎo)線覆冰與脫冰機(jī)理研究[D]. 湖北: 華中科技大學(xué), 2001.

Analysis of Ice-melting Process for Aerial Conductor

MAO Tian1, JIN Zi-xiang2, CHEN Shuang-tao*2
(1-Electric Power Research Institute, Yunnan Electric Power Test & Research Institute Group Co. Ltd., Kunming, Yunnan 650217, China; 2-School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China)

Ice-freezing has the negative influence on the electric transmission and transformation equipments of southern power grid in China. Considering the effects of heat convection, radiation heat transfer, gravity and detachment of water droplet from the wire surface, a short circuit ice-melting CFD model is set up in this paper. The ice-melting processes under different current intensity were simulated. The results show that, the unit time consuming and whole process time consuming did not decrease linearly with the linear grow of current intensity. The unit time consuming at the end of the ice melting process is the largest, and the same trend is more obvious in the lower current intensity ice-melting process.

Power grid; Ice-melting; CFD simulation; Heat flux

10.3969/j.issn.2095-4468.2015.06.102

*陳雙濤（1983-），男，副教授。研究方向：低溫工程。聯(lián)系地址：西安市咸寧西路28號西安交通大學(xué)能源與動力工程學(xué)院，郵編：710049。聯(lián)系電話：029-82664921。E-mail：stchen.xjtu@mail.xjtu.edu.cn。