王藝霖,吳 凡,李 智,王 泰
(1.東北電力大學(xué) 電氣工程學(xué)院,吉林 吉林 132012;2.吉林市源源熱電有限責(zé)任公司,吉林 吉林 132012)
電能在當(dāng)代是不可或缺的能源,隨著社會的不斷發(fā)展,對供電的可靠性要求也越來越高,然而大自然現(xiàn)象雷擊造成了電力系統(tǒng)經(jīng)常發(fā)生故障。其因是由于雷擊輸電線路產(chǎn)生很高的過電壓,使絕緣子串閃絡(luò),引起線路跳閘停電等事故[1-2]。因此當(dāng)今求解桿塔橫擔(dān)及絕緣子串上雷電過電壓是研究輸電線路耐雷性能的首要任務(wù)。
雷電放電時(shí),主放電發(fā)生過程之前,要先進(jìn)行先導(dǎo)放電,這個(gè)過程中存在大量的負(fù)電荷。此時(shí)會在地面感應(yīng)出大量的正電荷,兩種電荷快速運(yùn)動中會形成等離子區(qū),而且密度很高,并沿著先導(dǎo)通道快速向上傳播。因此可以把這個(gè)通道看成一個(gè)由電感和分布電容等值的導(dǎo)電通道,稱為雷電通道。把雷電放電過程用數(shù)學(xué)模型等值成一個(gè)電流源,其等值電路如圖1所示,根據(jù)電容和電感值可以求出它的波阻抗Z0,一般研究的雷電流幅值在100~300 kA。
圖1 雷電流源等值電路Fig.1 Equivalent circuit of lightning current source
雷擊輸電線路的主要方式包括反擊和繞擊[3]。反擊指雷擊避雷線或線路桿塔塔頂,繞擊指雷繞過避雷線直擊線路導(dǎo)線。反擊是造成線路故障跳閘的主要原因,其危害更大。研究輸電線路反擊耐雷性能的方法主要有區(qū)間組合法、蒙特卡洛法、行波法等。影響線路反擊耐雷性能的因素主要有絕緣子串伏秒特性、桿塔模型、輸電線路感應(yīng)過電壓、工頻電壓和接地電阻等[4]。
絕緣子串伏秒特性通過U50%放電電壓擬合得出,用以擬合的數(shù)據(jù)來源于實(shí)驗(yàn)得出的離散的數(shù)據(jù),通過拉格朗日插值多項(xiàng)式對其進(jìn)行曲線擬合,得到相對應(yīng)的伏秒特性曲線[5]。假設(shè)已知絕緣子串伏秒特性曲線上數(shù)據(jù)點(diǎn)(t0,u0)、(t1,u1)、(t2,u2)、…(tn,un),其中 un為 tn時(shí)刻所對應(yīng)電壓值。則它的n次插值基函數(shù)為:
根據(jù)式(1)可知,這個(gè)絕緣子串的伏秒特性曲線的n次插值多項(xiàng)式為
由上述插值式(1)、式(2)可知,所得到數(shù)據(jù)點(diǎn)越多,擬合出的伏秒特性曲線越準(zhǔn)確,精度越高。
根據(jù)中國高壓現(xiàn)行標(biāo)準(zhǔn),感應(yīng)過電壓Ui計(jì)算公式為[6]
式中:a為感應(yīng)過電壓的系數(shù),kV/m;hc為導(dǎo)線的平均高度,m;hg為避雷線的平均高度,m;k0為避雷線和導(dǎo)線間的耦合系數(shù)。
另外,線路桿塔模型一般被看成各段為分布參數(shù)的桿塔,即為單相無損線,因此通過貝杰龍算法等值得出計(jì)算模型。
當(dāng)發(fā)生反擊時(shí),雷電流會在桿塔各處產(chǎn)生很高的過電壓,當(dāng)這個(gè)過電壓等于絕緣子U50%沖擊閃絡(luò)電壓時(shí),即可求出對應(yīng)的耐雷水平I。
貝杰龍(Bergeron)算法是行波法的推廣延伸,該方法是將分布參數(shù)線路波過程的求解方法和集中參數(shù)電路暫態(tài)過程求解方法相結(jié)合所得到的數(shù)學(xué)模型化簡計(jì)算方法(即特性線法和梯形法相結(jié)合)[7]。利用線路的波過程將線路化簡成只含有電阻和電流源的集中參數(shù)網(wǎng)絡(luò),進(jìn)而求解網(wǎng)絡(luò)的暫態(tài)過程。在求解網(wǎng)絡(luò)暫態(tài)過程時(shí),首先將連續(xù)的時(shí)間轉(zhuǎn)換成離散的時(shí)間間隔,從發(fā)生擾動的初始時(shí)刻t0開始計(jì)算,一個(gè)時(shí)間一個(gè)時(shí)間進(jìn)行迭代,每一次計(jì)算時(shí),都利用t時(shí)刻之前的狀態(tài)作為下一時(shí)刻的初始狀態(tài),計(jì)算出這一時(shí)刻的電壓和電流狀態(tài)。周而復(fù)始,得出所需的電壓值和電流值。
在求解雷電暫態(tài)時(shí),將輸電線路桿塔各段都看成分布參數(shù)線路,即為單相無損線。其對應(yīng)的貝杰龍等值電路如圖2所示。設(shè)線路兩端的節(jié)點(diǎn)分別為k和 m,其對應(yīng)的電壓、電流分別為 uk(t)、ikm(t)、um(t)、imk(t),如圖 2(a)所示。通過線路波過程將其化簡成集中參數(shù)電路,如圖2(b)所示。
圖2 單相無損線路貝杰龍等值電路Fig.2 Equivalent circuit of single-phase lossless line of Bergeron
其中受控電流源值為:
在求解雷電反擊在桿塔上產(chǎn)生的雷電過電壓時(shí),首先利用以上等值化簡方法,將整個(gè)輸電線路桿塔等值成只包括集中參數(shù)電阻(其阻值等于線路波阻抗Z)和等值受控電流源的網(wǎng)絡(luò),確定網(wǎng)絡(luò)節(jié)點(diǎn),進(jìn)而生成網(wǎng)絡(luò)所對應(yīng)的節(jié)點(diǎn)導(dǎo)納矩陣Y,在雷擊的初始狀態(tài),電流源列向量只有雷電流值I,其余電流源初值均為0,進(jìn)行節(jié)點(diǎn)電壓方程求解,并隨著時(shí)間的推進(jìn)不斷更新電流源列向量,以進(jìn)一步求解網(wǎng)絡(luò)每個(gè)時(shí)間所對應(yīng)的各節(jié)點(diǎn)電壓。
節(jié)點(diǎn)電壓方程為
當(dāng)雷擊輸電線路桿塔或附近避雷線時(shí),其對應(yīng)桿塔的絕緣子串會承受很高的過電壓,這個(gè)過電壓包括以下四個(gè)分量即
式中:uins(t)為絕緣子串上過電壓;ucr(t)為橫擔(dān)上電壓;ui(t)為感應(yīng)過電壓;upf(t)為導(dǎo)線上工頻電壓;uco(t)為的耦合過電壓。
當(dāng)絕緣子串端電壓大于絕緣子串的U50%放電電壓時(shí),絕緣子串發(fā)生閃絡(luò)。絕緣子串發(fā)生閃絡(luò)具體時(shí)間的判定方法如圖3所示。
圖3 絕緣子串發(fā)生閃絡(luò)的時(shí)間Fig.3 Insulator string flashover time
當(dāng)絕緣子串上的電壓波形高于絕緣子串的伏秒特性波形時(shí),認(rèn)為絕緣子串此時(shí)發(fā)生閃絡(luò),這樣就基本能夠得出絕緣子串閃絡(luò)的時(shí)間。
直流輸電線路電壓等級高,容量大,傳輸能力強(qiáng),近年來在中國得到了比較廣泛的應(yīng)用[8-9]。桿塔模型及其等值模型如圖4所示。
圖4 桿塔模型及其等值模型Fig.4 Tower model and equivalent model
根據(jù)貝杰龍等值模型進(jìn)行VB編程計(jì)算,將得到絕緣子串兩端電壓,與絕緣子串的U50%放電電壓對比,從而判斷絕緣子串是否閃絡(luò)。具體流程圖如圖5所示。
圖5 輸電線路貝杰龍算法反擊計(jì)算程序框圖Fig.5 Counter calculation program diagram of Ber geron algorithm on transmission lines
為驗(yàn)證VB程序計(jì)算反擊耐雷性能的準(zhǔn)確性,以內(nèi)蒙古東部地區(qū)±500 kV伊穆直流線路152號桿塔型號為J51-33為例,驗(yàn)證模型的準(zhǔn)確性,其桿塔塔高為43.5 m,呼高為33 m,接地電阻9.7 Ω選取波形是2.6/50 μs、幅值為100~240 kA的負(fù)極性雷電流來進(jìn)行計(jì)算。VB程序所得出187 kA時(shí)絕緣子串兩端電壓值如圖6所示。EMTP軟件仿真所得出187 kA時(shí)絕緣子串兩端電壓值如圖7所示。
圖6 VB程序所得絕緣子串兩端電壓值Fig.6 Insulator string voltage value obtained at both ends by using VB program
圖7 EMTP所得絕緣子串兩端電壓值Fig.7 Insulator string voltage value obtained at both ends by using EMTP
從圖6、圖7可以看出,VB程序所求的絕緣子串兩端電壓波形和EMTP仿真所求波形在形狀趨勢幅值上大體相同,且二者所求耐雷水平均為188 kA,可見VB程序所求結(jié)果基本準(zhǔn)確,誤差較小。
桿塔接地電阻的大小直接影響桿塔的耐雷水平,原因是由于接地電阻越大塔頂電位越高,絕緣子串越容易閃絡(luò),現(xiàn)將152號桿塔J51-33分別取其接地電阻為5 Ω、10 Ω、15 Ω、20 Ω、25 Ω,求解其所對應(yīng)的耐雷水平,如表1所示。
表1 不同接地電阻下桿塔耐雷水平Tab.1 Tower lightning withstand level under different grounding resistance
由表1可知,隨著接地電阻的增加,桿塔的耐雷水平逐漸降低,且接地電阻越大時(shí),接地電阻的變化幅度對桿塔耐雷水平的影響也越大。
1)本文使用貝杰龍算法利用VB編程求解了500 kV直流輸電線路的反擊耐雷水平,誤差較小。
2)隨著接地電阻的增加,桿塔的耐雷水平逐漸減弱,故為保證線路的耐雷性能,其接地電阻不宜大于 10 Ω。
[1] 魯鐵成.電力系統(tǒng)過電壓[M].北京:中國水利水電出版社,2009.LU Tiecheng.Power system overvoltage[M].Beijing:China Water& Power Press,2009.
[2] 谷定燮.500 kV輸變電工程設(shè)計(jì)中雷電過電壓問題[J].高電壓技術(shù),2000,26(6):60 -62.GU Dingxie.The research on lightning overvoltage of 500 kV transmission engineering design [J].High Voltage Engineering,2000,26(6):60 -62.
[3] 吳健兒.淺談輸電線路雷電繞擊與反擊[J].浙江電力,2007,26(4):43-46.WU Jianer.Discussion on lightning shielding failure and counterattack of transmission lines[J].Zhejiang Electric Power,2007,26(4):43-46.
[4] 魯志偉,馬文婧.新立變電站接地系統(tǒng)安全性能研究[J].東北電力大學(xué)學(xué)報(bào),2012,32(6):28 -32.LU Zhiwei,MA Wenjing.Research on safety performance of substation grounding system[J].Journal of Northeast Dianli University,2012,32(6):28 -32.
[5] 崔濤.輸電線路防雷計(jì)算中絕緣子串閃絡(luò)判據(jù)研究[D].武漢:華中科技大學(xué),2009.CUI Tao.Research on insulation flashover criterion in lightning protection calculation for transmission lines[D].Wuhan:Huazhong University of Science and Technology,2009.
[6] 姜玉宏,常俊.配電線路感應(yīng)雷過電壓計(jì)算[J].黑龍江電力,2014,36(5):443 -446.JIANG Yuhong,CHANG Jun.Calculation of lightning overvoltage for distribution lines[J].Heilongjiang Electric Power,2014,36(5):443-446.
[7] 胡勁松.貝杰龍算法在墨江500 kV變電站雷擊過電壓計(jì)算中的應(yīng)用[D].成都:四川大學(xué),2005:2-10.HU Jinsong.Application of Bergeron method in lightning overvoltage Calculation for Mojiang 500 kV substation[D].Chengdu:Sichuan University,2005:2 -10.
[8] 張翠霞,葛棟,殷禹.直流輸電系統(tǒng)的防雷保護(hù)[J].高電壓技術(shù),2008,34(10):2070 -2074.ZHANG Cuixia,GE Dong,YIN Yu.Lightning protection of HVDC transmission system[J].High Voltage Technology,2008,34(10):2070-2074.
[9] 曾南超.高壓直流輸電在我國電網(wǎng)發(fā)展中的作用[J].高電壓技術(shù),2004,30(11):11 -12.ZENG Nanchao.Role of HVDC transmission in the power system development in China[J].High Voltage Technology,2004,30(11):11-12.