基于ATP－EMTP的配電線路耐雷性能計算與仿真分析

龍浩然，朱立波

（南寧供電局，廣西 南寧 530031）

1 引言

在噴射氣流條件下，由于氣流場與電磁場等的耦合作用，卡西（Cassie）電弧模型和麥也爾（Mayr）模型的能量耗散等參數(shù)均發(fā)生了相應的變化。因此需要通過理論建模分析在強氣流噴射條件下，氣體對間隙電弧的強烈消電離作用，高速氣流對電弧溫度的快速冷卻作用，高速氣流對電弧中導電離子的加速擴散作用等，并建立適合噴射氣流條件下的電弧模型。

由于電弧模型只能做相應的定性分析，本文提出的噴射氣流滅弧防雷間隙電弧計算模型并借鑒了K.Hori－nouchi的假定，將電弧假定為一條由若干電弧片斷連接而成的鏈條，電弧片斷的個體運動決定了電弧鏈的整體運動。同時引入弧柱的三層物理結(jié)構(gòu)，沿電弧徑向定量求解溫度場，依據(jù)噴射氣流條件下電弧溫度場的變化求解電弧的特性參數(shù)［1－3］。

2 噴射氣流條件下間隙電弧的流體動力學模型

電弧運動是電流、電磁場、氣流場及熱輻射等相互作用的結(jié)果。初始的溫度和壓力分布決定了等離子體中的電導率，通過施加合理的電位邊界條件，求解電磁場方程，便可確定出電位分布；結(jié)合電導率分布，就決定了電弧中的電流密度分布；同時，電弧作為載流體，一方面，將產(chǎn)生一定的焦耳熱，該值由電導率和電流密度分布決定，另一方面，在外加噴射氣流的作用下會受到噴射力的作用，從而使電弧運動；這樣最終導致等離子體內(nèi)部的能量傳輸，從而引起等離子體參數(shù)的變化。

本文在以往工作的基礎上，對基于歐拉方程組的介質(zhì)特性進行分析計算，并對各種條件進行了研究，建立了以流體動力學為基礎的噴射氣流條件下間隙電弧動態(tài)模型的研究。流體動力學電弧模型以流體動力學方程（Navier－storkes）為基礎，由能量方程、動量力程、安培定律、質(zhì)量連續(xù)性方程、磁感應強度等多方面構(gòu)成。

首先假定間隙中的氣體流動滿足局部熱力學平衡條件，在本文中用柱坐標下軸對稱可壓縮的非定常歐拉方程組來描述流場。則噴射氣流條件下氣體流動控制方程組如下所示［1］:

質(zhì)量守衡方程:

徑向動量守衡方程:

軸向動量守衡方程:

能量守衡方程:

氣體的狀態(tài)方程:

式中:ρ—氣體密度；

v，u—氣體的軸向和徑向速度；

P—氣體壓力；

E—氣體比總能，E=e+0.5（u2+v2）；

Q—體積能源密度項；

T—氣體溫度。

2.1 單個電流元受力分析

如果將單個電流元視為一質(zhì)量為零的直棒導體，根據(jù)經(jīng)典激波理論，當它以速度va運動時，在其前面會產(chǎn)生一個激波波震面。根據(jù)氣體動力學中的質(zhì)量、動量和能量守衡定律，并結(jié)合氣體狀態(tài)方程，可得:

式中:Co為空氣中的音速，331.2m/s；y為激波強度；μ=（γ－1）（γ+1）；γ=1.4。

電流元將受到三個力的作用:氣體噴射力Fp，洛侖磁力Fm和空氣阻力Fs。忽略電流元的質(zhì)量，則電流元的受力情況可以描述為:

式中Po為標準大氣壓0.1013MPa；B為電流元重心處的磁感應強度；r和L分別為電流元的半徑和長度；I為電弧電流；F為氣體發(fā)生器出口噴射力；a為氣體噴射角。

將式（1）與式（2）結(jié)合，可以得到電流元在空氣中的速度，如式（3）所示，速度方向與洛倫磁力方向相同。由式（3）可知，電流元的速度和電流元的長度無關(guān)。

2.2 電流元徑向溫度分析

根據(jù)噴射氣流滅弧防雷間隙電弧在強氣流噴射條件下的物理過程，將空氣中的電流元劃分為3個區(qū)域（見圖1）。

圖1 電流元的區(qū)域劃分

區(qū)域1為導電圓柱體區(qū)域，在區(qū)域內(nèi)電導率δ為溫度的函數(shù)；

區(qū)域2為圍繞區(qū)域1的熱邊界層，在區(qū)域內(nèi)電導率為零；

區(qū)域3為包圍在區(qū)域2外的空氣，溫度恒為環(huán)境溫度300K。

電流元的半徑為區(qū)域1的外邊界半徑。

根據(jù)噴射氣流滅弧防雷間隙電弧在強氣流噴射條件下的運動過程，對所研究的電流元及其周圍的空氣作如下假定:

（1）在強氣流噴射前，所研究的場域內(nèi)氣體處于局部熱力學平衡狀態(tài)，因此可以用統(tǒng)一的熱力學溫度來描述電流元及其周圍空氣；

（2）電流元為軸對稱圓柱體，因此可用二維圓柱坐標描述；

（3）電流元在強氣流噴射條件下運動時，氣體大部分從電流元外面繞流經(jīng)過，穿過電流元內(nèi)部的氣體較少，可以忽略不計；

（4）電流元各參數(shù)沿軸向的變化可忽略不計；

（5）電流元內(nèi)部沿軸向的熱傳遞可忽略不計；

（6）通過電流元軸向的磁感應強度和電場強度保持不變，以重心處數(shù)值表示。

在噴射氣流滅弧防雷間隙裝置中，氣體發(fā)生器產(chǎn)生的氣流將橫吹電弧，如果直接將電弧假設為圓柱坐標、軸對稱，則將會引起較大的誤差。采用鏈式電弧模型，將電弧分成若干非常短的電流元，由氣體橫吹引起的電弧彎曲被反映為電弧鏈形狀的改變和電流元不斷的分裂為多個電流元，于是由此假設所引起的誤差就隨著電流元長度的減小而減小。將同樣長度的電弧分成的電流元數(shù)目越多，假設所引起的誤差越小。

3 數(shù)學模型控制方程

描述電流元物理特性的基本方程是能量守恒方程。對應電流元的不同區(qū)域，能量守恒方程的表達形式略有不同。在區(qū)域1中，一方面外電路輸入焦耳熱，另一方面由于噴射氣流產(chǎn)生傳導、對流和輻射而散失能量。對于區(qū)域2，沒有焦耳熱輸入，但由于強氣流迫使電弧運動，從而強制電弧對流、熱傳導和熱輻射。這一區(qū)域從區(qū)域1帶走大部分熱量，并輸送到區(qū)域3中。區(qū)域3是一個開域，溫度恒為環(huán)境溫度。

相應于區(qū)域1和區(qū)域2的能量守恒方程見式（4）和式（5）。

區(qū)域1:

區(qū)域2:

式中:T—電流元弧柱溫度；

t—時間變量；

r—徑向變量；

ρ—空氣密度；

Cp—定壓比熱；

δ—電導率；

E—電場強度；

k—熱導率；

u—輻射系數(shù)；

α—對流換熱系數(shù)［4－5］。

針對模型控制方程進行求解計算。電弧特性參ρ、Cp、δ、k和u的變化是受到溫度、壓力等諸多因素的綜合影響的，但本文在對控制方程進行求解的過程中，均假定它們是溫度T的單值函數(shù)，這種假定較之電弧特性參數(shù)為常數(shù)的假定更符合電弧物理現(xiàn)象的實際情況，對控制方程的求解又進行了必要的簡化。在具體的計算中，參數(shù)ρ、Cp、δ、k隨溫度變化的曲線從參考文獻［5］獲得，輻射系數(shù)u根據(jù)參考文獻［3］提供的曲線確定。對流換熱系數(shù)α根據(jù)圓柱體外橫向繞流的原理，在已知電流元運動速度的情況下通過雷諾準則、努塞爾準則等關(guān)系式求得［6－8］。

4 基于Mayer電弧模型的ATP－EMTP仿真

4.1 仿真模型建立

本節(jié)針對雷擊過電壓的不同保護方式，立足于分析開關(guān)電弧的電壓和電流進行了研究，根據(jù)噴射氣流條件下改進的Mayer電弧模型方程建立了基于ATPEMTP的仿真模型，以不同電壓等級電流接地系統(tǒng)，對防雷間隙進行了仿真。ATP－EMTP的仿真模型如圖2。

圖2 基于ATP－EMTP的防雷保護間隙電弧仿真模型

4.2 仿真結(jié)果

圖3 間隙工頻電弧電流、電壓波形

圖4 噴射氣流速度V=850m/s，P=15MPa時的電流、電壓波形

圖5 噴射氣流速度V=2800m/s，P=34MPa時的電流、電壓波形

圖6 噴射氣流速度V=4600m/s，P=46MPa時的電流、電壓波形

5 結(jié)論

通過ATP－EMTP仿真結(jié)果表明:當單相接地故障電弧的電容電流經(jīng)消弧線圈補償，故障點殘流為20～30A時，噴射氣流速度能夠在20～40ms內(nèi)完全將間隙電弧吹滅。仿真結(jié)果驗證了在高速氣體的沖擊下，在開放式大氣環(huán)境中，噴射氣流滅弧過程中電弧弧柱溫度、電流、電場強度和平均電導率等參數(shù)的變化情況，強氣流對電弧強烈的消電離作用，對電弧溫度的快速冷卻作用以及對電弧中導電離子產(chǎn)生加速擴散作用。使得電弧的溫度持續(xù)降低直至電弧熄滅。從而能夠較好的保護絕緣子和導線免受電弧灼燒，起到降低雷擊跳閘率的作用。

［1］谷山強，何金良，陳維江，等.架空輸電線路并聯(lián)間隙防雷裝置電弧磁場力計算研究［J］.中國電機工程學報，2006，26（7）:140 －145.

［2］Yos J M.Transport properties of nitrogen，hydrogen，oxygen and air to 30000K.Wilmington Massachusetts:Research and Advanced Development Division AvcoCcorp，1963:39 －40.

［3］Shayler P J，F(xiàn)ang M T C.Radiation transport in wall－ stabilised nitrogen arcs.J Phys D:Appl Phys，1978，11:1743 －1756.

［4］馬桂珍.真空電弧最新控制技術(shù)［J］.真空電子技術(shù)，1998，4:18－20.

［5］袁海文，陳德桂，辜曉川.微型低壓限流器跑弧區(qū)磁場的計算及其對限流性能的影響［J］.西安交通大學學報，1996，30（3）:19－25.

［6］張晉，陳德桂，付軍.三相短路情況下塑殼斷路器開斷特性的計算機模擬［J］.低壓電器，1999，1:10 －14.

［7］袁海文，李安云，李榮源.電弧通道磁場分布對斷弧效果的控制作用研究［J］.北京航空航天大學學報，2002，28（3）:280 －282.

［8］陳旭，陳德桂.新型滅弧系統(tǒng)對背后擊穿現(xiàn)象的抑制作用［J］.低壓電器，2000，3:7 －9.