劉 坤,崔運(yùn)光,王 牣
基于傳輸線法的高速鐵路貫通地線暫態(tài)研究
劉 坤,崔運(yùn)光,王 牣
提出了基于傳輸線法建立高速鐵路貫通地線的等效電路模型,根據(jù)電磁場理論確定等效電路模型中各分布參數(shù)具體數(shù)值的計算方法。在ATP-EMTP中計算分段導(dǎo)體長度不同時,2根平行貫通地線上軸向電壓峰值分布情況。通過與 CDEGS中仿真結(jié)果的對比,得到雷電流作用下基于傳輸線法的貫通地線等效電路模型中分段導(dǎo)體長度選取的最佳值。研究結(jié)果表明當(dāng)分段導(dǎo)體長度為4 m時,計算誤差最小,能夠滿足工程應(yīng)用的要求,也驗(yàn)證了本文研究方法和研究結(jié)果的有效性。
高鐵防雷;傳輸線法;貫通地線;均勻分段導(dǎo)體;軸向電壓
雷擊鐵路沿線接觸網(wǎng)、建筑物等裝置時會在接地系統(tǒng)中產(chǎn)生高幅值的過電壓,威脅鐵路沿線設(shè)備和人員的安全。貫通地線可以為雷電流提供有效的泄流通道,所以,研究貫通地線在雷電流作用下暫態(tài)特性對解決高速鐵路防雷與接地相關(guān)問題具有重大意義[1]。迄今為止,通過建立貫通地線等效電路模型分析其雷擊暫態(tài)特性的研究還很少,本文借鑒電力系統(tǒng)中變電站接地網(wǎng)的研究成果對貫通地線進(jìn)行暫態(tài)模型的建立與研究。
目前接地網(wǎng)暫態(tài)特性分析方法主要有3種:電路法[2]、電磁場法[3]和傳輸線法[4]。其中電路法基于準(zhǔn)靜態(tài)假設(shè),將接地網(wǎng)簡化為一組集中參數(shù)和分布參數(shù)電路模型,簡單易懂,但是該方法只適用于 暫態(tài)電流頻率不高、接地網(wǎng)規(guī)模較小的情況。電磁場法基于麥克斯韋方程組對暫態(tài)過電壓進(jìn)行精確求解,忽略的因素最少,但其建模復(fù)雜,計算量非常大。傳輸線法是基于以上方法的改進(jìn),建立分布參數(shù)模型,該方法建模靈活,計算效率高,可以有效分析規(guī)模較大的接地網(wǎng)雷擊暫態(tài)響應(yīng)特性,但基于傳輸線法研究貫通地線雷擊暫態(tài)特性時存在分布參數(shù)如何計算和分段導(dǎo)體長度難以確定的問題。
本文首先介紹傳輸線理論以及國內(nèi)高速鐵路復(fù)線線路路基段貫通地線的實(shí)際埋設(shè)方式。提出基于傳輸線法建立符合真實(shí)情況的貫通地線等效電路模型,利用電磁場理論的相關(guān)知識計算等效模型中的各個分布參數(shù)。針對傳輸線法中分段導(dǎo)體長度難以確定的問題,在電磁暫態(tài)計算軟件ATP-EMTP中計算得到雷電流作用下傳輸線模型中分段導(dǎo)體長度取不同值時,2根平行貫通地線上軸向電壓峰值的分布情況。通過將 CDEGS的仿真結(jié)果與EMTP中的計算結(jié)果進(jìn)行對比驗(yàn)證,得到基于傳輸線法的等效模型中分段導(dǎo)體選取的最佳長度,同時也證明了基于傳輸線法對貫通地線進(jìn)行暫態(tài)模型建立與研究的有效性。
1.1 傳輸線法
傳輸線理論用“電路”中的電壓和電流等效傳輸線中的電場和磁場[5]。埋設(shè)于土壤中的水平接地體在低頻電流作用下,可將其等效為對地電導(dǎo)。而在高頻電流作用下,接地體上會出現(xiàn)分布參數(shù)效應(yīng)。圖1為均勻傳輸線上任意一段分段導(dǎo)體dz的等效電路,可以在線上任一點(diǎn)z處取dz進(jìn)行暫態(tài)電壓、電流的討論。
圖1 分段導(dǎo)體dz的等效電路模型圖
圖1中,R1,L1分別為均勻傳輸線分段導(dǎo)體的電阻和電感,G0,C0分別為均勻傳輸線分段導(dǎo)體的電導(dǎo)和電感,I為分段導(dǎo)體上流過的電流,U為分段導(dǎo)體在z處的電壓。
1.2 基于傳輸線法的貫通地線模型
由貫通地線的實(shí)際埋設(shè)方式可知,貫通地線可以視為無限長水平接地體,在激勵源為雷電流的情況下,必須考慮接地體中電流的不均勻分布和波傳播的過程。將貫通地線等效為均勻有損長線,考慮平行貫通地線間的等效電容和等效電導(dǎo),將貫通地線進(jìn)行均勻分段處理,則貫通地線可視為由無數(shù)個均勻分段導(dǎo)體串聯(lián)而成的分布參數(shù)電路,基于上述傳輸線法建立平行雙根貫通地線的等效電路模型如圖2所示??紤]到銅的導(dǎo)磁率非常?。ń咏谡婵眨?,且雷電流的等值頻率不超過1 MHz[6],同時2根貫通地線間的距離較遠(yuǎn),因此本研究忽略了平行分段導(dǎo)體間的互感作用。
圖2中,R1,L1和R2,L2分別為2根貫通地線分段導(dǎo)體自身的電阻和電感,G0,C0分別為貫通地線分段導(dǎo)體的對地電導(dǎo)和對地電感,G12,C12分別為平行貫通地線分段導(dǎo)體間的等效電導(dǎo)與等效電容,I1,I2分別為兩平行分段導(dǎo)體上流過的電流,U1,U2分別為兩平行分段導(dǎo)體在 z處對地的電壓。分段導(dǎo)體參數(shù)值可根據(jù)電磁場理論的相關(guān)知識求得。
圖2 平行雙根貫通地線等效電路模型圖
2.1 貫通地線相關(guān)參數(shù)
國內(nèi)高速鐵路綜合接地系統(tǒng)中,貫通地線采用銅絞線制成,忽略集膚效應(yīng)對導(dǎo)線自身電阻的影響,則R1,R2可用下式求出:
式中,ρ為銅的電阻率,l為分段導(dǎo)體長度,S為貫通地線的橫截面積,一般情況下,S取70 mm2。
分段導(dǎo)體自身電感L1,L2可用下式求出:
式中,μ0為空氣導(dǎo)磁率,a為貫通地線的半徑。
在計算分段導(dǎo)體對地電導(dǎo)G0時,用接地體的接地電阻來近似表示其對地電導(dǎo),因?yàn)橥ǔKf的接地電阻包括接地引線的電阻、接地引線與接地裝置的接觸電阻、接地體本身的電阻,接地體和土壤間的接觸電阻以及土壤的散流電阻,其中散流電阻遠(yuǎn)遠(yuǎn)大于其他4種電阻的總和。根據(jù)平均電位法可推導(dǎo)出埋深為h的水平接地體的接地電導(dǎo):
對地電容根據(jù)靜電場中已知的電容公式得到:
式中,ρs,εs分別為土壤的電阻率和介電常數(shù)。
采用鏡像原理移去土壤與空氣的分界面,則分段導(dǎo)體間的等效電容為
等效電導(dǎo)依舊可以通過式(4)比擬得到。
實(shí)際中貫通地線可以視為無限長水平接地導(dǎo)體,而這在 EMTP的計算中是無法實(shí)現(xiàn)的,且雷電流在貫通地線中的有效散流長度較長,理論上貫通地線上軸向電壓可以無限衰減下去直至無窮小,所以在不影響計算精度的前提下,只選取有限長度的貫通地線并在其末端串接上一個適當(dāng)大小的等效阻抗來近似無限長貫通地線。文獻(xiàn)[7]研究并得到了高速鐵路綜合接地系統(tǒng)中貫通地線的接地阻抗的計算表達(dá)式:
由式(6)可知,當(dāng)貫通地線的半徑、埋設(shè)深度和土壤電阻率保持不變時,接地電阻會隨著貫通地線的增長而減小,接地電抗則反之,且當(dāng)貫通地線超過5 km時,二者的值都會趨于穩(wěn)定值。實(shí)際中貫通地線沿鐵路埋設(shè),其長度遠(yuǎn)遠(yuǎn)大于5 km,所以在傳輸線模型串接的末端阻抗就可以取上述研究中得到的穩(wěn)定值。
在基于傳輸線法的貫通地線等效電路模型中,分段導(dǎo)體長度的選取對貫通地線上軸向電壓分布計算精度的影響比較大。文獻(xiàn)[8]將接地體視為數(shù)個導(dǎo)體元件串聯(lián)而成的均勻分布參數(shù)電路,研究了對不同埋深、不同土壤電阻率和不同雷電流參數(shù)下接地體的有效散流長度,但并未給出研究中分段導(dǎo)體長度選取的具體依據(jù)。原則上,分段導(dǎo)體尺寸越小,即分段數(shù)越大,計算結(jié)果越接近于真實(shí)情況。但同時,除了計算量會隨著分段數(shù)的增加而急劇增加外,由于參數(shù)計算基于諸多假設(shè),計算精度也有可能會隨著累積誤差的增加而減小。另一方面,如果分段導(dǎo)體尺寸很大,即分段數(shù)很小,分布參數(shù)模型就會很接近集中參數(shù)模型,從而限制其在沖擊電流及高頻電流條件下的應(yīng)用。文獻(xiàn)[9]認(rèn)為導(dǎo)體尺寸的選取取決于注入雷電流的等值頻率中的最大頻率f,土壤電阻率ρs,土壤介電常數(shù)εs和空氣導(dǎo)磁率μ0。當(dāng)分段導(dǎo)體長度l滿足關(guān)系(2a / 10)≤l≤(λ / 10)時,可以忽略分段導(dǎo)體上的推遲效應(yīng),即均勻分段導(dǎo)體上各點(diǎn)電流與電位都相同。其中a是分段導(dǎo)體的半徑,λ是頻率為 f的正弦電流作用于埋設(shè)在土壤中的導(dǎo)體上時,其周圍產(chǎn)生的電磁波波長。
基于上述公式,取分段導(dǎo)體長度分別為2,3,4,5,6,8 m,計算得出不同分段導(dǎo)體長度下傳輸線模型的相關(guān)參數(shù)。
2.2 雷電流參數(shù)及波形
取土壤電阻率ρs= 100 Ω·m,土壤介電常數(shù)εs= 10×8.854×10-12F/m,空氣導(dǎo)磁率μ0= 4π×10-7H/m。貫通地線即銅的電阻率ρ = 0.017 5 Ω·mm2/m,貫通地線半徑 a = 4.72 mm。雷電流采用幅值為10 kA,波頭時間2.6 μs,半峰值時間50 μs的標(biāo)準(zhǔn)雷電波。雷電流等效波形采用雙指數(shù)函數(shù)波形,其表達(dá)式為
式中,I0為雷電流峰值,α,β為常數(shù),其值與雷電流相關(guān)參數(shù)有關(guān)。
目前,相關(guān)研究通常使用雷電流等值頻率這一概念,實(shí)際上雷電流是由無數(shù)個不同大小的頻率分量組成,其中含有少量的高頻分量。雷電流通過貫通地線時,正是由于其中含有不同大小的頻率分量,會對均勻分段導(dǎo)體長度的選取產(chǎn)生較大影響。
電磁暫態(tài)計算軟件ATP-EMTP以節(jié)點(diǎn)法及差分法為計算基礎(chǔ),完全適用于基于傳輸線法的2根埋地平行貫通地線的等效電路模型。雷電流從線路中點(diǎn)注入,由于電路模型具有對稱性,所以只需要計算注入點(diǎn)一側(cè)的電位分布情況。由于均勻分段導(dǎo)體長度不同,在 EMTP中電壓分布情況的計算點(diǎn)也不同。
通過在CDEGS中搭建實(shí)際模型對各種電力系統(tǒng)接地問題進(jìn)行仿真分析,其通用性和準(zhǔn)確性已經(jīng)在工程實(shí)踐中得到了驗(yàn)證[10]。本文利用CDEGS仿真得到的結(jié)果作為真實(shí)情況與基于傳輸線法得到的計算結(jié)果進(jìn)行對比分析與驗(yàn)證。
當(dāng) 10 kA標(biāo)準(zhǔn)雷電流通過接地端子侵入高速鐵路綜合接地系統(tǒng)中的貫通地線后,產(chǎn)生最高約為48 kV的過電壓,雷電流在近似無限長的貫通地線中傳播,傳播過程中產(chǎn)生的過電壓不斷衰減,衰減速率隨著傳播距離的增加而減小,同時會在與之平行的另一根貫通地線上感應(yīng)出相當(dāng)幅值的過電壓,最高約為 9 kV。相比于直接注入雷電流的貫通地線上過電壓的衰減速率,感應(yīng)過電壓的衰減速率明顯減小。隨著與雷電流注入點(diǎn)的距離變大,電壓波形上升陡度略有減小,下降陡度急劇減小。將EMTP中計算得到的直接注入雷電流的貫通地線和產(chǎn)生感應(yīng)電動勢的貫通地線上軸向電壓峰值分布情況與CDEGS中的仿真結(jié)果進(jìn)行對比,見圖3。
圖3中,曲線組1表示直接注入雷電流的貫通地線上軸向電壓峰值分布情況的結(jié)果對比,曲線組2表示產(chǎn)生感應(yīng)電動勢的貫通地線上軸向電壓峰值分布情況的結(jié)果對比。通過對比分析可知,均勻分段導(dǎo)體長度的選取對產(chǎn)生感應(yīng)電動勢的貫通地線上軸向電壓峰值分布影響較小,而對直接注入雷電流的貫通地線上軸向電壓峰值分布影響較大,本文認(rèn)為這與雷電流的頻率組成有極大的關(guān)系,在傳播過程中,雷電流中含有的高頻分量相比于低頻分量更容易衰減,所以在距離雷電流注入點(diǎn)較近的區(qū)域,高頻分量衰減是雷電流衰減的主要因素,而在距離雷電流注入點(diǎn)較遠(yuǎn)的區(qū)域,低頻分量衰減是雷電流衰減的主要因素。由圖3可知當(dāng)均勻分段導(dǎo)體長度取4 m時,計算誤差最小,直接注入雷電流的貫通地線上軸向電壓峰值誤差不超過10%,感應(yīng)電壓峰值誤差不超過16%,且偏差值穩(wěn)定。
圖3 不同長度分段導(dǎo)體電壓峰值分布情況對比圖
表1為均勻分段導(dǎo)體長度取4 m時直接注入雷電流貫通地線上軸向電壓峰值結(jié)果對比,表2為感應(yīng)電壓峰值結(jié)果對比。
表1 直擊過電壓峰值對比表
表2 感應(yīng)電動勢電壓峰值對比表
綜上分析,得出如下結(jié)論:
(1)基于傳輸線法可以有效建立符合真實(shí)情況的高速鐵路貫通地線等效電路模型,明確了等效電路模型中各分布參數(shù)的計算方法。
(2)當(dāng)雷電流通過接地端子侵入高速鐵路綜合接地系統(tǒng)中后,會在相互平行的貫通地線上產(chǎn)生高幅值的過電壓,雷電流在近似無限長的貫通地線中傳播,傳播過程中過電壓不斷衰減,衰減速率隨著傳播距離的增加而減小,相比于直接注入雷電流的貫通地線上過電壓的衰減速率,感應(yīng)電壓的衰減速率明顯減小。
(3)當(dāng)均勻分段導(dǎo)體長度取4 m時,計算結(jié)果誤差最小,直接注入雷電流的貫通地線上軸向電壓峰值的計算誤差不超過10%,感應(yīng)電壓峰值誤差不超過16%,且偏差值維持相對穩(wěn)定。
[1] 楊崗. 客運(yùn)專線綜合接地系統(tǒng)方案研究[J]. 鐵道工程學(xué)報,2006,(7):76-80.
[2] 張麗萍,袁建生,李中新.變電站接地網(wǎng)不等電位模型數(shù)值計算[J].中國電機(jī)工程學(xué)報,2000,20,(1):1-3.
[3] Arnautovski-Toseva V,Grcev L.High frequency current distribution in horizontal grounding systems in two-layer soil[C].IEEE International Symposium on Electromagnetic Compatibility,Boston,2003.
[4] 齊磊,崔翔,趙志斌.復(fù)雜接地系統(tǒng)沖擊接地特性的數(shù)值計算及試驗(yàn)[J].電網(wǎng)技術(shù),2006,30,(13):66-69.
[5] Tapan K M.Impulse impedance of grounding systems and its effects on tower cross-arm voltage[D].Tennessee:Tennessee Technological University,2006.
[6] 張波,崔翔,趙志斌,等. 大型變電站接地網(wǎng)的頻域分析方法[J]. 中國電機(jī)工程學(xué)報,2002,22,(9):59-63.
[7] 高國強(qiáng),董安平,張雪原,等. 高速鐵路綜合地線的接地效果[J]. 西南交通大學(xué)學(xué)報,2011,46,(1):103-108.
[8] He Jinliang,Gao Yanqing,Zeng Rong,et al.Effective length of counterpoise wire under lightning current[J].IEEE Trans. on Power Delivery,2005,20,(2):1585-1591.
[9] 楊琳,吳廣寧,曹曉斌. 變電站接地網(wǎng)暫態(tài)建模方法[J].電網(wǎng)技術(shù),2012,36,(5):161-165.
[10] 周建飛,周年光,孔淑琴,等. 接地金屬網(wǎng)對高壓輸電線路下方工頻電場強(qiáng)度及感應(yīng)電壓改善研究[J].電網(wǎng)技術(shù),2011,35,(7):199-202.
The equivalent current model of high-speed running-through earth wire based on transmission line method is put forward, and the calculation method of different distributed parameters in equivalent current model is defined by electromagnetic field theory. The axial peak voltage distributions of the two parallel running-through earth wires are calculated in ATP-EMTP when segmented conductors are different in lengths. The optimum values for selecting lengths of segmented conductors in equivalent current model based on transmission line method can be obtained by comparison between calculation results and CDEGS simulation results. The research results show that the calculation error is smallest when the length of segmented conductor is 4m and is able to meet the requirements of engineering applications. And these also verify that the research methods and results in this paper are effective.
Lightning protection of high-speed railway; transmission line method; running-through earth wire; uniform segmented conductor; axial voltage
U226.8+3
:B
:1007-936X(2015)03-0031-05
2014-09-28
劉 坤.西南交通大學(xué)電氣工程學(xué)院,碩士研究生,電話:18208160585;
崔運(yùn)光.西南交通大學(xué)電氣工程學(xué)院,碩士研究生;
王 牣.西南交通大學(xué)電氣工程學(xué)院,教授。