楊蕓,高超,陳楓林,孫瑞杰,李銳海,王黎明

(1.南方電網(wǎng)科學(xué)研究院,廣州 510080；2.清華大學(xué)深圳研究生院,深圳 518055)

棒板長間隙擊穿電壓計算

楊蕓1,高超1,陳楓林2,孫瑞杰1,李銳海1,王黎明2

(1.南方電網(wǎng)科學(xué)研究院,廣州 510080；2.清華大學(xué)深圳研究生院,深圳 518055)

利用模擬電荷法提出了一種仿真模型,模擬了1 m到10 m棒板長間隙放電過程中電暈、流注、先導(dǎo)和末躍的基本過程。確定了先導(dǎo)判據(jù)和流注長度經(jīng)驗公式中的參數(shù),得到了電暈電壓、先導(dǎo)長度、末躍電壓和擊穿電壓的仿真數(shù)據(jù),與經(jīng)驗公式和試驗結(jié)果相符。并在最后歸納出了一個擬合公式,用于判斷1 m到10 m任意間隙長度的擊穿電壓。

棒板間隙；長間隙放電；擊穿電壓；先導(dǎo)發(fā)展模型；模擬電荷法

0 前言

研究電力線路外絕緣的設(shè)計需要有數(shù)量龐大且精確可靠的長空氣間隙放電參數(shù)作為依據(jù),需要開展可靠的長空氣間隙放電試驗[1]。但進(jìn)行實際大小的長間隙放電試驗需要大量資金、考慮復(fù)雜的試驗環(huán)境條件和全面的間隙類型以及試驗儀器的精度和測量,采用實際大小的線路來試驗非常困難。因此,利用先進(jìn)的計算機(jī)技術(shù)和功能強(qiáng)大的仿真軟件,通過數(shù)學(xué)建模,將長空氣間隙的物理過程合理地表達(dá),為研究開辟了一條新的路徑[2]。將理論和仿真結(jié)果相結(jié)合,并配合部分試驗,以此分析并推測各種間隙類型和電壓形式下的空氣放電特性,找到其中的普遍規(guī)律[3]。

棒-板間隙是典型間隙,在輸電線路的絕緣配合中,絕緣子長串與地之間、絕緣子長串與鐵塔之間以及電力變壓器套管與地之間的放電情形都最接近棒-板長間隙放電模型。由于棒-板間隙屬于極不均勻間隙,其擊穿電壓遠(yuǎn)遠(yuǎn)低于其它的間隙類型。綜上,對棒-板間隙的研究具有代表性和實際意義,因此國內(nèi)外專家均把它作為最基本的間隙放電模型進(jìn)行研究[4-7]。

目前國內(nèi)外在間隙擊穿的數(shù)值計算中,對電場采用的數(shù)值分析方法主要有：有限元法、有限差分法、蒙特卡羅法和模擬電荷法等等。建立的典型計算模型主要有：空間電荷模型、能量平衡模型、等離子體模型、電暈云模型和先導(dǎo)發(fā)展模型等等[8]。

1 程序流程與模型設(shè)計

1.1 程序流程

考慮到棒板長間隙放電發(fā)展的物理過程,將其分為四個階段。第一至四階段分別為電暈階段、流注階段、先導(dǎo)階段、末躍階段。每個階段采用各自模型,階段之間由判據(jù)連接。

電壓從初值開始,當(dāng)不滿足判據(jù)條件時,電壓增加,直到能依次通過各個判據(jù),程序結(jié)束,得到最終擊穿電壓。

1.2 模型設(shè)計

棒電極由3組,每組600個電荷來模擬。棒電極直徑0.038 m,長6 m,截面上,電荷在三條互成120度半徑的二分之一處,距圓心0.009 5 m；軸向上,電荷之間距離0.01 m。在圓柱表面取600個參考點(diǎn),每個參考點(diǎn)相距0.01 m,每個參考點(diǎn)的電壓均取外加電壓。

設(shè)流注也成圓柱形,與棒電極有相同的橫截面,用3組電荷模擬,排布與電極模型相似。流注長度用經(jīng)驗公式：

確定S為間隙長度,a為調(diào)整系數(shù)。軸向電荷的數(shù)量與流注長度有關(guān),間隔0.01 m。在流注中心取參考點(diǎn),參考點(diǎn)間隔0.01 m。參考點(diǎn)電壓由近電極端到遠(yuǎn)電極端遞減,考慮到流注內(nèi)場強(qiáng)約4 kV/cm[11],與流注頭部場強(qiáng)約25 kV/cm[12]。最左端參考點(diǎn)電壓取棒電極電壓,向右逐點(diǎn)遞減,先按照4 kV/cm遞減,再按照25 kV/cm遞減,則可得各參考點(diǎn)電壓。

設(shè)先導(dǎo)也成圓柱形,與電極有相同的橫截面,3組電荷模擬,排布與棒電極模型相似。先導(dǎo)的長度會變化,每當(dāng)軸向長度增加,先導(dǎo)的模擬電荷數(shù)量也增加,整個模型的電荷量重新分布。在先導(dǎo)中心取參考點(diǎn),參考點(diǎn)間隔0.01 m。參考點(diǎn)電壓由近電極端到遠(yuǎn)電極端遞減,考慮到先導(dǎo)內(nèi)場強(qiáng)約1.5 kV/cm[14],最左端參考點(diǎn)電壓取棒電極電壓,按照1.5 kV/cm向右逐點(diǎn)遞減,則可得各參考點(diǎn)電壓。

2 程序設(shè)計

2.1 模擬電荷量計算

圖1 棒電極模擬電荷分布

設(shè)圓柱截面半徑為r,q11距p1長度p=r/2,設(shè)棒電極右端為x=0,間隙長度為l。則軸線上電荷q(i)的坐標(biāo)為d(i)=-0.01?i+0.005.

由對稱性知q11=q12=q13,設(shè)其電量均為q1。

由疊加原理,對600個未知電荷 [q1q2…q600]’與600個參考點(diǎn)電勢 [u1u2…u600]’有方程組：

對參考點(diǎn)pi計算電勢：

pi受qj1影響的電勢為：

pi受qj1影響的電勢為：

pi受qj2和qj3影響的電勢均為：

pi受qj2和qj3的鏡像電荷影響的電勢均為：

系數(shù)矩陣中,系數(shù)ai,j為：

2.2 流注計算

如果電壓值足夠使程序進(jìn)入流注階段,則根據(jù)流注長度經(jīng)驗公式與流注內(nèi)部電壓梯度的特征,一次性計算流注存在后的電場情況和流注最前端電勢。并以此判斷完整的流注是否能存在即能否開始發(fā)展先導(dǎo)。

2.3 先導(dǎo)計算

如果電壓值足夠使程序進(jìn)入先導(dǎo)階段,則先導(dǎo)的起始長度即為原流注的長度。先導(dǎo)的計算是分步進(jìn)行的,在一個電壓值下,若電場與帶電離子濃度滿足判據(jù)條件,則先導(dǎo)增加單位長度,重新進(jìn)行模擬電荷的分布,重新計算電場與帶電離子濃度。若滿足判據(jù)條件,則重復(fù)上述過程；若不滿足判據(jù)條件,則增加外加電壓值,重復(fù)上述過程。直到先導(dǎo)最前端電勢滿足發(fā)生末躍的條件,跳出循環(huán)。

2.4 帶電粒子濃度計算

電暈電流與該處電場E有關(guān)系式[16],其中h是常數(shù),取值10-11(Am2/V-2),Ec是先導(dǎo)起始的電場判據(jù),取值31 kV/cm。電子崩集中發(fā)生在1 cm的距離里,而流注的發(fā)展速度約為108cm/ s[17],即可認(rèn)為頭部區(qū)域儲存電荷的時間約是t= 1/108=10-8s。再由帶電離子數(shù)目的計算公式n= It/q0得到。

3 參數(shù)和判據(jù)的選擇

3.1 流注長度經(jīng)驗公式

流注長度經(jīng)驗公式為

其中l(wèi)為間隙長度,a為間隙系數(shù),與間隙類型與間隙長度有關(guān)。在3 m到10 m的間隙長度下,對a從0.5到1取值,觀察計算結(jié)果。

表1 a=0.6仿真結(jié)果

表2 a=0.8仿真結(jié)果

表3 a=1仿真結(jié)果

比較發(fā)現(xiàn)當(dāng)間隙長度在5 m到10 m范圍內(nèi)時,取a=0.8是比較合理的。a=0.8時,在整個長度范圍內(nèi)誤差均較小,且比較平均。因此仿真取a=0.8.

同理,當(dāng)間隙長度小于5 m時,取a=0.5比較合理。

3.2 帶電離子濃度

先導(dǎo)放電維持條件對帶電粒子濃度的要求是1 013/cm3到1 015/cm3,在這之間都是合理的。因此在程序中,我引入了可變的與流注頭部電勢有關(guān)的調(diào)整系數(shù),使仿真過程中擊穿電壓更符合實際。

表4 帶電粒子濃度系數(shù)的選擇對擊穿電壓的影響

上表中的誤差指的是與Gallet計算公式：

計算電壓之間的誤差。

3.3 流注先導(dǎo)起始場強(qiáng)判據(jù)

流注的起始場強(qiáng)選用臨界場強(qiáng)23 kV/cm,先導(dǎo)的起始場強(qiáng)選用臨界場強(qiáng)31 kV/cm。

3.4 氣壓修正系數(shù)

氣壓越高,擊穿電壓越低。擊穿電壓受氣壓的影響有公式：

其中U50為氣壓為p時的擊穿電壓,P0取101.3 kPa,n為氣壓特征指數(shù),與間隙距離有關(guān),且正極性受影響的程度大于負(fù)極性。

氣壓與海拔高度有關(guān)系式：

正極性下,高度每升高1 km,擊穿電壓約降低8.7%,負(fù)極性下約為8%。因此,取氣壓修正指數(shù)為：b=1-0.087h,其中h為海拔高度,單位km。

4 正極性仿真結(jié)果與分析

4.1 仿真結(jié)果

表5 正極性仿真結(jié)果

4.2 仿真結(jié)果分析

4.2.1 擊穿電壓

表6 正極性擊穿電壓分析

未在2.5 m、4 m、6 m和8 m間隙長度做試驗。

觀察上表,仿真值與經(jīng)驗公式計算值,以及氣壓修正后與試驗值基本相符,誤差在合理范圍內(nèi)。

從曲線圖上看,仿真值在各理論值的上下限范圍內(nèi),仿真結(jié)果比較精確,具有可信度。

圖2 擊穿電壓仿真值擬合曲線

二次函數(shù)擬合公式為：

4.2.2 先導(dǎo)長度

試驗表明,在能出現(xiàn)先導(dǎo)的長間隙放電中,先導(dǎo)長度通常發(fā)展到間隙的三分之二左右[9],便發(fā)生末躍,間隙擊穿。

表7 正極性先導(dǎo)長度分析

根據(jù)上表統(tǒng)計,先導(dǎo)長度占間隙的60%-80%,在三分之二左右,與試驗現(xiàn)象相符。

4.2.3 末躍電壓

試驗表明,在能出現(xiàn)先導(dǎo)的長間隙放電中,先導(dǎo)上的電壓降通常達(dá)到擊穿電壓的百分之五十到六十,即末躍電壓大概是擊穿電壓的百分之四十到五十,便發(fā)生末躍,間隙擊穿。

表8 正極性先導(dǎo)長度分析

根據(jù)上表統(tǒng)計,末躍電壓占間隙擊穿電壓的40%～60%居多,與試驗相符。

5 負(fù)極性仿真結(jié)果與分析

5.1 仿真結(jié)果

表9 負(fù)極性仿真結(jié)果

5.2 仿真結(jié)果分析

5.2.1 擊穿電壓

表10 負(fù)極性擊穿電壓分析

仿真值與理論值相符,誤差在合理范圍內(nèi)。

5.2.2 極性效應(yīng)

表11 極性效應(yīng)

由上表可知,正極性流注起始電壓高于負(fù)極性流注起始電壓,正極性擊穿電壓遠(yuǎn)遠(yuǎn)低于負(fù)極性擊穿電壓。與極不均勻間隙的極性效應(yīng)的理論相符,程序具有可信度。

6 結(jié)束語

本文建立的模型合理且比較準(zhǔn)確,其優(yōu)勢在于模擬了長間隙放電過程的各個階段,且以先導(dǎo)發(fā)展為重心,突出了先導(dǎo)在長間隙放電中的重要性。得到的數(shù)據(jù)包括流注起始電壓、流注長度、先導(dǎo)長度、末躍電壓和擊穿電壓,非常全面。但此模型只能模擬最基本的情況,且主要用于1 m到10 m的間隙長度。在模型中沒有考慮溫度與濕度或者雨天、酸雨、山火等更為復(fù)雜的環(huán)境條件,也沒有考慮先導(dǎo)發(fā)展路徑的曲折。通過改變電極模型,應(yīng)該能推廣模擬其他典型間隙的放電發(fā)展；改變一些參數(shù),應(yīng)該也能推廣模擬其他電壓形式下的放電發(fā)展。這些都是可以改進(jìn)和完善的地方。

[1] Garrara,G＆Thione,L.Switching surge strength of large air gaps：A physical approach[J].IEEE Trans,1976, Vol.95,2(1)：512-524

[2] 唐耀宗,陳亞珠.對Carrara-Thione的長空氣間隙正極性操作波放電模型的評論 [J].高壓電器,1984(1)：41 -48.

[3] E.O.Selim,R.T.Waters,et al.Space Charge Modeling in Impulse Corona[J].Proceeding of the 4th ISH.1983,41：0-2.

[4] Vereshchagin I.P,Beloglovsky A.A,Vinokurov V.N,et al. Model of Impulse Streamer Corona Formation[C].Eleventh International Symposium on High Voltage Engineering,1991, 3：248-251.

[5] R.T.Waters.A Thermodynamic Model of the Leader in LongAir Gaps[J].Research on Long Air Gaps Discharge at Renardieres 1973 Result,Electra,1974,35：124-129

[6] 萬啟發(fā),霍鋒,謝梁,等.長空氣間隙放電特性研究綜述[J].高電壓技術(shù),2012,38(10)：2499-2505.

[7] 丁玉劍,李慶峰,廖蔚明,等.高海拔地區(qū)典型長空氣間隙的操作沖擊放電特性和海拔校正 [J].高電壓技術(shù), 2013(06)：1141-1146.

[8] 文習(xí)山,陳慈萱,解廣潤.長間隙放電過程的物理模型[J].高電壓技術(shù),1990(2)：1-6.

[9] Hutzler,B.,Hutzler-Barre,D.Leader Propagation Model for Predetermination of Switching Surge Flashover Voltage of Large Air Gaps[J].IEEE Trans.,1978,Vol.PAS-97,4：1087-1096.

[10] 馬乃祥.長間隙放電 [M].北京：中國電力出版社,1998.

[11] Les RENARDIERES GROUP.Research on long air gap discharges at Les RENARDIERES[J].1973 Results,ELECTRA N.35 July 1974：49-156.

[12] 王文端,張海燕.長間隙放電流注發(fā)展的新的計算模型[J].現(xiàn)代電力,1995(3)：8-14.

[13] 洪川.沖擊電壓下棒-板長空氣間隙放電中空間電場的計算研究 [D].2007,重慶大學(xué).

[14] 胡建林,等.基于有限元弱解式的棒-板長空氣間隙先導(dǎo)放電空間電場仿真研究 [J].中國電機(jī)工程學(xué)報, 2008(10)：148-154.

[15] H.N.GARCIA,B.HUTZLER.Electrical breakdown in long air gaps：The final jump[C].3rd International Conference on gas discharges LONDON Sept.1974：206-210.

[16] J.S.TOWNSEND.Electricity in gases[M].Oxford University Press.NEW-YORK,1914.

[17] 梁曦東.高電壓工程 [M].清華大學(xué)出版社.

[18] 張海燕,王文端,先導(dǎo)放電發(fā)展模型的研究 [J].現(xiàn)代電力,1996(1)：37-41.

[19] 洪川,胡建林,孫才新,等.棒-板長空氣間隙擊穿電場穩(wěn)態(tài)條件的仿真 [J].高電壓技術(shù),2007,33(11)：75 -79.

[20] 蔣興良,等.棒-板長空氣間隙在低氣壓下雷電沖擊特性及電壓修正 [J].中國電機(jī)工程學(xué)報,2005(25)：152 -156.

陳楓林 (1989)，男，博士研究生，清華大學(xué)深圳研究生院，從事高電壓技術(shù)方面研究工作。

Research on Calculation of Rod-plane Long Air Gaps Breakdown Voltage

YANG Yun1,GAO Chao1,CHEN Fenglin2,SUN Ruijie1,LI Ruihai1,WANG Liming2
(1.Electric Power Research Institute,CSG,Guangzhou,510080,China；2.Graduate School at Shenzhen,Tsinghua University,Shenzhen,Guangdong 518055,China)

In this paper,charge simulation method is used to propose an emulation model that simulates the basic course of corona, stream,leader and last step in discharge process of rod-plate long air gap from 1 m to 10 m.Parameters of leader criterion and stream length empirical formula are fixed,and simulation results of corona voltage,leader length,last step voltage and breakdown voltage are obtained,conforming to the empirical formula and experiment results.At last,a fitting formula is generalized to judge the breakdown voltage of all the air gaps from 1 m to 10 m.

rod-plane gap；long air gap discharge；breakdown voltage；leader develop model；charge simulation method

TM73

B

1006-7345(2015)01-0145-06

2014-10-11

楊蕓 (1981)，女，碩士，工程師，南方電網(wǎng)科學(xué)研究院，從事高電壓技術(shù)方面研究工作 (e-mail)yangyun＠csg.com。

高超 (1983)，男，碩士，高級工程師，南方電網(wǎng)科學(xué)研究院，從事高電壓試驗技術(shù)方面研究工作 (e-mail)gaochao＠csg.com。