基于場路耦合低值絕緣子串溫度分布研究

夏 凡

（三峽大學(xué) 電氣與新能源學(xué)院，湖北宜昌 443002）

引言

絕緣子是輸電線路的重要組成部分，其對保障輸電線路電氣性能起關(guān)鍵性作用，其絕緣有效性關(guān)乎輸電線路的安全運行。瓷質(zhì)絕緣子在長期運行過程中容易受到機(jī)電聯(lián)合作用發(fā)生機(jī)械性能和絕緣性能下降的情況，這種現(xiàn)象容易產(chǎn)生低值絕緣子。雷電、雨天時，低值絕緣子容易受到雷擊，其表面將產(chǎn)生強(qiáng)大的電流，電流的熱效應(yīng)容易使低值絕緣子發(fā)生斷裂事故[1-3]，對電力系統(tǒng)的穩(wěn)定運行造成嚴(yán)重威脅。

文獻(xiàn)[4]研究了絕緣子性能劣化時，其溫度將發(fā)生改變，為研究絕緣子熱故障提供了理論支持；文獻(xiàn)[5,6]研究了絕緣子表面電導(dǎo)率對絕緣子串溫度分布的影響；文獻(xiàn)[7]針對干燥帶在絕緣子表面發(fā)熱機(jī)理進(jìn)行了研究，并得出干燥帶長度、寬度以及位置對絕緣子溫度分布影響較大；文獻(xiàn)[8]對絕緣子表面積污進(jìn)行了溫度場研究；文獻(xiàn)[9]研究表明，當(dāng)濕度發(fā)生變化時，絕緣子串溫度也會隨之改變；文獻(xiàn)[10]利用電網(wǎng)絡(luò)法對絕緣子的紅外熱像檢測盲區(qū)進(jìn)行了研究，量化了盲區(qū)的概念；文獻(xiàn)[11]針對零值絕緣子紅外檢測盲區(qū)問題，提出了應(yīng)用盤面特征識別來檢測零值絕緣子。

上述研究絕緣子紅外熱像問題時，用傳統(tǒng)有限元軟件考慮桿塔計算對低值絕緣子串溫度分布的影響較為復(fù)雜，計算成本高，且較少考慮多種因素的影響；而單純利用電網(wǎng)絡(luò)法參考絕緣子雜散電容經(jīng)驗值進(jìn)行仿真存在一定問題，難以適應(yīng)實際現(xiàn)場各種各樣的環(huán)境和絕緣子串模型。故本研究以有限元軟件與電網(wǎng)絡(luò)法相結(jié)合，以桿塔7片絕緣子為例，首先通過Maxwell電容矩計算110 kV絕緣子雜散電容，結(jié)合電網(wǎng)絡(luò)利用MATLAB計算絕緣子阻值對絕緣子串電壓分布影響，進(jìn)而討論低值絕緣子阻值及其位置，風(fēng)速、測量誤差對含低零值絕緣子溫度分布檢測的影響。

1 絕緣子雜算電容理論分析

以常見的瓷絕緣子為例，每片絕緣子由鋼帽、鋼腳、傘裙組成，水泥粘合劑連接其組合部分。其中鋼帽、鋼腳作為絕緣子中導(dǎo)體部分，在計算絕緣子串雜算電容時可以將其視為一個導(dǎo)體。在一個絕緣子片數(shù)為n的塔-絕緣子-導(dǎo)線系統(tǒng)中，每片絕緣子鋼帽和上一片絕緣子鋼腳相連，絕緣串上金具和桿塔相連，絕緣子串下金具和導(dǎo)線相連。為方便計算雜算電容，可以分別將其視為導(dǎo)體，由于桿塔為零電位，可以將桿塔為參考地。在絕緣子片數(shù)為n的塔-絕緣子-導(dǎo)線系統(tǒng)中有n個導(dǎo)體，假設(shè)任意一個導(dǎo)體上的電荷為Qi，電位為Ui，其中i=1、2、…、n。Ci0為導(dǎo)體i對地電容，Cij為導(dǎo)體i和導(dǎo)體j之間相互電容。所以在塔-絕緣子-導(dǎo)線系統(tǒng)中，各導(dǎo)體上電荷、電容、電位之間滿足如下關(guān)系[12]：

將式（1）以矩陣形式表述為：

矩陣C中有：

矩陣C對角線上元素的值Cii為相應(yīng)導(dǎo)體對系統(tǒng)中所有其他導(dǎo)體的部分電容和。故導(dǎo)體i對地電容為：

假設(shè)高壓導(dǎo)線為導(dǎo)體n，則導(dǎo)體i(i≠n)對導(dǎo)線n的電容為Cin，即各片絕緣子對導(dǎo)線的雜算電容為矩陣Cin。

2 110 kV直線貓頭塔瓷絕緣子串模型

2.1 模型建立

110 kV直線貓頭塔帶中相瓷絕緣子串的桿塔選擇B7型直線桿塔，以圓柱體分別模擬金具和導(dǎo)線，絕緣子采用 xp-160瓷質(zhì)絕緣子，xp-160 瓷質(zhì)絕緣子具體技術(shù)參數(shù)見表1，各部分材料介質(zhì)的相對介電常數(shù)見表2[13]。

表1 XP-160絕緣子技術(shù)參數(shù) mm

表2 各介質(zhì)材料相對介電常數(shù)

2.2 雜散電容計算

采用comsol軟件計算麥克斯韋電容矩陣。為計算麥克斯韋電容矩陣，需要對各導(dǎo)體電壓進(jìn)行設(shè)置，由于各導(dǎo)體電容大小與自身結(jié)構(gòu)和材料屬性有關(guān)，與加在導(dǎo)體上電壓高低無關(guān)，可對桿塔設(shè)置零電位，其他導(dǎo)體設(shè)置2 V，結(jié)合公式（1）～（4）計算麥克斯韋電容矩陣。其中以距離導(dǎo)線最近絕緣子為1號絕緣子，其余絕緣子編號以此類推，通過仿真計算，各片絕緣子對地電容和對導(dǎo)線電容分別如圖1、2所示。

圖1 對地電容分布

圖2 對高壓導(dǎo)線電容分布

從圖1、2可以看出，絕緣子對地雜散電容和對導(dǎo)線雜散電容值總體滿足電容值與距離呈非線性關(guān)系。絕緣子距離桿塔橫擔(dān)、導(dǎo)線距離越遠(yuǎn)，雜散電容越小，其中絕緣子對地雜散電容最大值與絕緣子上方金具長短有關(guān)，金具越長代表絕緣子最上片與橫擔(dān)距離越大，相應(yīng)絕緣子對桿塔雜散電容越??；同理絕緣子對導(dǎo)線雜散電容最大值與絕緣子下方金具長短有關(guān)，絕緣子下方金具越長，相應(yīng)絕緣子對導(dǎo)線雜散電容越小。

2.3 絕緣子串路模型及其電壓分布驗證

單片絕緣子可以等效為電容和電阻的并聯(lián)，故110 kV桿塔7片絕緣串可建立如圖3所示的等效模型。

圖3 絕緣子等效電路圖

依據(jù)圖3絕緣子串等效電路圖，在MATLAB建立絕緣子等效電路模型，參考文獻(xiàn)[14]中測得絕緣子自身電容為50 pF，雜散電容由comsol電場計算獲得。絕緣子在潔凈情況下，電阻很大，取絕緣子自身電阻為10 GΩ，絕緣子等效電路模型中U設(shè)置為63.5 kV，利用MATLAB計算110 kV桿塔7片絕緣子串電壓分布，并與實際測量值、有限元仿真軟件仿真值進(jìn)行對比，如圖4所示。從圖4可以看出，MATLAB計算值、實際測量值和有限元計算絕緣子電壓分布值均滿足U型分布，并且這3種計算結(jié)果均相差無幾，從而驗證場路耦合在計算絕緣子分布電壓具有可行性和正確性。

圖4 與實測電壓分布比較

3 低值絕緣子對絕緣子溫度分布影響

在低于500 kV輸電線路中，當(dāng)瓷懸式絕緣子被檢測電壓低于50%標(biāo)準(zhǔn)規(guī)定值或者所測量絕緣子阻值低于300 MΩ時，將被判定為劣化絕緣子[15]。取靠近導(dǎo)線側(cè)絕緣子為1號絕緣子，分別取1號絕緣子為不同阻值，計算不同阻值對絕緣子串電壓分布的影響，如圖5所示。

圖5 1號絕緣子不同絕緣子阻值對電壓分布的影響

從圖5可以看出，當(dāng)1號絕緣子阻值為300 MΩ時，絕緣子串電壓分布和正常絕緣子串電壓分布并無明顯差別，即300 MΩ并不能作為判定絕緣子發(fā)生劣化的依據(jù)，當(dāng)1號絕緣子阻值為25 MΩ時，1號絕緣子分布電壓相對正常絕緣子的分布電壓下降約為50%，即該值可以作為絕緣子發(fā)生劣化判據(jù)。

3.1 絕緣子發(fā)熱理論分析

內(nèi)部穿透性電流發(fā)熱、介質(zhì)損耗電流發(fā)熱和表面泄露電流發(fā)熱是絕緣子發(fā)熱的3個主要部分[16]，當(dāng)絕緣子潔凈時，絕緣子溫升主要由介質(zhì)損耗電流發(fā)熱功率決定；當(dāng)絕緣子表面污穢層濕潤時，此時絕緣子發(fā)熱主要由絕緣子表面泄露電流決定，但介質(zhì)損耗電流發(fā)熱功率仍然不可忽略。根據(jù)牛頓冷卻定律，絕緣子相對環(huán)境溫升可以表達(dá)為：

式中：A為散熱面積，在該模型里指絕緣子鋼帽面積，對于XP-160型絕緣子鋼帽表面積A＝0.03 m2；h為散熱系數(shù)，當(dāng)溫度在5℃～35℃、風(fēng)速在0～0.2 m時，h為15 W/m2·℃。

3.2 低值絕緣子阻值對絕緣子溫度分布的影響

分別對1號絕緣子取不同阻值，如圖6所示。從該圖可以看出，絕緣子整體溫度分布呈現(xiàn)非線性分布，當(dāng)1號絕緣子發(fā)生劣化時，對應(yīng)溫度高于整串絕緣子中其他位置絕緣子。其中，當(dāng)1號絕緣子阻值為50 MΩ對應(yīng)溫度最高，參考值曲線中正常絕緣子溫度最低。

圖6 1號絕緣子不同阻值溫度分布

3.3 低值絕緣子位置對絕緣子溫度分布影響

按照距離高壓導(dǎo)線遠(yuǎn)近，取三個特殊位置絕緣子，分別為1號、4號和7號絕緣子，并令其電阻值為50 MΩ，得到同一阻值不同位置劣化絕緣子對絕緣子串溫度分布影響，如圖7所示。從圖7可以看出當(dāng)絕緣子發(fā)生劣化時，絕緣子串總體溫升會抬升。這是因為絕緣子劣化時，絕緣子串總體電阻減少，導(dǎo)致絕緣子串溫度抬升。絕緣子串發(fā)生劣化時，相對環(huán)境而言，首端絕緣子溫度最高，最高可達(dá)3.4℃，尾端絕緣子溫度次之，為2.5℃，中部絕緣子溫度最低，最低為1.5℃。

圖7 劣化絕緣子位置對絕緣子串溫度影響

3.4 風(fēng)速對低值絕緣子溫度影響

散熱系數(shù)h與風(fēng)速v有關(guān)，為研究風(fēng)速對低值絕緣子溫度影響，通過改變散熱系數(shù)來模擬風(fēng)速的變化，如圖8所示。從圖8可以看出，風(fēng)速越大，對絕緣子溫度分布影響越大，當(dāng)風(fēng)速v=4 m/s時，4號劣化絕緣子溫升為0.3℃，正常絕緣子溫升為0.2℃，當(dāng)考慮測量誤差時，容易將正常絕緣子誤判為零值絕緣子。

圖8 風(fēng)速對絕緣子串溫度影響

3.5 測量誤差對低值絕緣子溫度影響

在絕緣子紅外成像檢測過程中存在多種因素影響，導(dǎo)致絕緣子實際溫度與檢測溫度存在測量誤差，以風(fēng)速v=2 m/s，4號劣化絕緣子阻值為50 MΩ為例，通過改變測量誤差研究測量誤差對劣化絕緣子串溫度的影響，如圖9所示。

圖9 測量對絕緣子串溫度影響

通過圖9可知，測量誤差對絕緣子溫度分布影響較大，當(dāng)測量誤差為ΔT=0.3℃時，4號劣化絕緣子溫升為0.45℃，正常絕緣子溫升為0.1℃。

4 結(jié)論

（1）本研究提出的場路耦合法能夠根據(jù)實際絕緣子桿塔模型計算出雜散電容，綜合考慮低值絕緣子阻值、位置以及風(fēng)速和測量誤差多種因素影響，進(jìn)而計算低值絕緣子串溫度分布，相比傳統(tǒng)有限元仿真，更加準(zhǔn)確高效。

（2）當(dāng)?shù)椭到^緣子于導(dǎo)線側(cè)且阻值為50 MΩ時，低值絕緣子溫度最高，易于檢測。

（3）當(dāng)考慮風(fēng)速對低值絕緣子串溫度影響時，風(fēng)速超過2 m/s，絕緣子串溫度分布變化明顯，低值絕緣子檢測難度大。

（4）當(dāng)考慮測量誤差時，測量誤差越大，計算的低值絕緣子溫度較低，容易造成誤檢。