雷電光電離過(guò)程的仿真研究

龍浩然，朱立波

(南寧供電局，廣西壯族自治區(qū) 南寧 530031)

然而光電離不僅可由外來(lái)入射光的作用而產(chǎn)生，也可由氣體放電內(nèi)部的光輻射而產(chǎn)生。上述方程忽略氣體內(nèi)部的光電離作用，Kline等補(bǔ)充了這一點(diǎn)列出了下列的方程［5－7］:

雷電光電離過(guò)程的仿真研究

龍浩然，朱立波

(南寧供電局，廣西壯族自治區(qū) 南寧 530031)

考慮氣體中的光電離，列出電荷密度方程，建立Matlab數(shù)學(xué)模型，仿真出雷電放電的電荷密度變化，解釋雷電流最大陡度出現(xiàn)時(shí)間為放電一段時(shí)間后的峰值前，而不是在放電最開(kāi)始。

光電離;電荷密度;雷電流最大陡度;峰值

1 引言

雷電流的變化率(di/dt)是人們對(duì)雷電流研究工作中最關(guān)心的一個(gè)參數(shù)，因?yàn)槔纂娏鞫付仍酱螅瑢?duì)電氣設(shè)備造成的危害也越大。目前在各種文獻(xiàn)中最常使用的是雙指數(shù)函數(shù)的雷電波形，它的雷電流變化率是由大變小的，也就是在雷電放電開(kāi)始時(shí)陡度最大。而近年來(lái)觀測(cè)得到大多數(shù)的第一次主放電電流波形在其上升到幅值之前時(shí)比較緩慢，然后再轉(zhuǎn)入陡的部分，雙指數(shù)波形已不符合實(shí)際雷電放電波形，在工程計(jì)算中無(wú)法真實(shí)的體現(xiàn)雷電引起的電磁干擾水平。實(shí)質(zhì)上雷電放電時(shí)，雷電通道相當(dāng)于一條充滿(mǎn)了各種離子的等離子氣體溝，進(jìn)行著各種復(fù)雜的電離反應(yīng)，其中光電離對(duì)于雷電流迅速上升起著很重要的作用［1］。流注理論計(jì)算中考慮了光電離的作用，但因解方程式的困難而簡(jiǎn)化了光電離，使得放電發(fā)展的計(jì)算不理想［2，3］。近年來(lái)計(jì)算機(jī)迅速發(fā)展，出現(xiàn)了很多方便的計(jì)算軟件，計(jì)算光電離在放電中的作用已成為可能。本文從微觀等離子體研究角度，考慮了光電離在雷電放電過(guò)程中的重要作用，仿真出雷電放電波形，與實(shí)際放電波形相符合。

2 光電離計(jì)算模型

Davies1971年開(kāi)始采用特征法計(jì)算了Wagner的實(shí)驗(yàn)結(jié)果，把計(jì)算范圍擴(kuò)展到了更大的電流。電子和離子的密度變化由下列方程表示［4］:

其中:

ne，n+—電子及正離子的密度。

ve，v+—電子及正離子的漂移速度。

ve和電離系數(shù)α—?dú)鈮汉碗妶?chǎng)強(qiáng)度的函數(shù)，有如下關(guān)系:

然而光電離不僅可由外來(lái)入射光的作用而產(chǎn)生，也可由氣體放電內(nèi)部的光輻射而產(chǎn)生。上述方程忽略氣體內(nèi)部的光電離作用，Kline等補(bǔ)充了這一點(diǎn)列出了下列的方程［5－7］:

N是氣體密度，ψ=θ/αNΩ是光電離系數(shù)，由測(cè)量來(lái)決定:

上式中:rp=2rd，它表示0≤r≤2rd范圍內(nèi)光電離產(chǎn)生的離子對(duì)都加入放電的發(fā)展之中，均勻放電時(shí)取rp=rd。

在計(jì)算電子密度分布時(shí)，將雷電通道看成一長(zhǎng)度100m，(先導(dǎo)的電暈半徑約0.6～6m)通道半徑為1m的導(dǎo)體。x及y分別是圓柱坐標(biāo)系的軸向及徑向坐標(biāo)。x=0，r=0是放電通道軸心，rd是放電通道的半徑。假設(shè) 0≤r≤rd范圍內(nèi)，ne(r，x，t)及n+(r，x，t)是均勻的;r＞rd之外，電子及正離子密度等于零［8］。DL是電子縱向擴(kuò)散系數(shù)，氣體中的電子擴(kuò)散系數(shù)大約在105～106(cm2/sec)，相比較電離和漂移項(xiàng)小很多，可以忽略。

根據(jù)方程(3)，方程(5)可寫(xiě)為:

3 仿真結(jié)果

以上的計(jì)算都是在雷電流到達(dá)峰值前的一段時(shí)間里，也就是對(duì)于雷電波前波形的推導(dǎo)。

接著運(yùn)用matlab根據(jù)式(7)搭建數(shù)學(xué)模型進(jìn)行仿真。在仿真中有些數(shù)據(jù)做了一定的近似取值，但不影響它的大致走向。得出波前電荷分布密度的曲線如圖1所示。

圖1 電荷密度變化曲線

雷電流i的值與先導(dǎo)通道的電荷密度σ及主放電發(fā)展v速度的關(guān)系式為［3］:

由于R一般不超過(guò)30Ω，而雷電通道波阻抗Z0一般在300Ω以上，即R＜＜Z0，這樣得:

從表3還可看出，利用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行日負(fù)荷預(yù)測(cè)，對(duì)一天48個(gè)時(shí)段電動(dòng)汽車(chē)充電負(fù)荷值求平均相對(duì)誤差，值為5.25%；利用BP神經(jīng)網(wǎng)絡(luò)算法和指數(shù)平滑法進(jìn)行電動(dòng)汽車(chē)充電負(fù)荷的滾動(dòng)預(yù)測(cè)平均相對(duì)誤差為5.13%。

所以雷電流的增長(zhǎng)近似于電荷密度的增長(zhǎng)。將圖1波形進(jìn)行坐標(biāo)的轉(zhuǎn)換和一定的幅值處理，于是得到雷電流波前波形如圖2所示。

圖2 到達(dá)峰值前的波形

顯而易見(jiàn)，雷電流的最大增長(zhǎng)率出現(xiàn)在峰值前，當(dāng)先導(dǎo)走完雷電通道，則強(qiáng)烈的中和就出現(xiàn)，雷電流開(kāi)始變小。參照結(jié)合雙指數(shù)波形的后續(xù)波形，就可以得出新的符合實(shí)際放電情況的雷電波形，如圖3所示。

圖3 新雷電波形

圖3的新雷電波形將雷電放電的三個(gè)階段展現(xiàn)得很好，第一階段中先導(dǎo)放電階段的電流非常小，大約幾十安培，對(duì)應(yīng)于圖3波形中0～7μs波形;當(dāng)迎面先導(dǎo)與上行先導(dǎo)接通后進(jìn)入第二階段，也就是主放電階段，電流瞬間達(dá)到幅值20kA，對(duì)應(yīng)于8～10μs的波形;主放電后，電流開(kāi)始持續(xù)減小，進(jìn)入第三階段對(duì)應(yīng)于圖3中10μs以后的波形。

4 仿真效果分析

在得出新雷電波形基礎(chǔ)上，再將傳統(tǒng)使用的雙指數(shù)函數(shù)波形與其進(jìn)行對(duì)比。根據(jù)建筑物防雷設(shè)計(jì)規(guī)范國(guó)家標(biāo)準(zhǔn)CB50057－94本文采用10/35μs波形，即波頭時(shí)間 τ1為10μs，波長(zhǎng)時(shí)間 τ2為 350μs，本文取 α =1/τ2，β =1/τ1，Im取其典型值為 2 × 104A［4］。兩個(gè)波形對(duì)比見(jiàn)圖4。

分析發(fā)現(xiàn)雙指數(shù)函數(shù)波形存在著幾點(diǎn)嚴(yán)重不足之處:

(1)幅值在表達(dá)式中無(wú)法體現(xiàn)，且峰值遠(yuǎn)達(dá)不到幅值數(shù)值。

雷電流波形雙指數(shù)函數(shù)表達(dá)式為:

i=I0(e－αt－e－βt)

根據(jù)式子我們?cè)O(shè)其峰值時(shí)間為tm，峰值電流為Im，則有

顯然式中的I0并不是雷電流的幅值，表示的只是某一時(shí)刻的雷電流值，它具體的物理意義不明確。而通常在使用雙指數(shù)波時(shí)卻將I0認(rèn)為是雷電流的幅值Im，這是不對(duì)的。

從圖4雙指數(shù)波形中可明顯看到，它的峰值在17kA左右，遠(yuǎn)達(dá)不到20kA的幅值標(biāo)準(zhǔn)。

圖4 與雙指數(shù)波形的對(duì)比

在雷擊時(shí)室內(nèi)磁場(chǎng)分布的估算中若采用雙指數(shù)雷電流波形會(huì)因峰值不夠、波頭陡度不足而低估室內(nèi)可能出現(xiàn)的磁場(chǎng)水平，且雙指數(shù)函數(shù)在t=0處沒(méi)有連續(xù)的一階導(dǎo)數(shù)，在雷電電磁場(chǎng)計(jì)算中也帶來(lái)諸多不便。

5 結(jié)論

本文仿真出雷電放電時(shí)的電荷密度變化，很好的詮釋了雷電流最大陡度出現(xiàn)時(shí)間為放電一段時(shí)間后的峰值前，而不是在放電最開(kāi)始。新雷電流的表達(dá)式可根據(jù)電荷密度變化曲線結(jié)合實(shí)際觀測(cè)波形進(jìn)行擬合得出，對(duì)解決電力系統(tǒng)中由于陡度而造成的雷害具有重要的參考價(jià)值。

［1］ 谷山強(qiáng)，何金良，陳維江，等．架空輸電線路并聯(lián)間隙防雷裝置電弧磁場(chǎng)力計(jì)算研究［J］．中國(guó)電機(jī)工程學(xué)報(bào)，2006，26(7):140 －145．

［2］ Yos JM．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 TC．Radiation transport in wall-stabilised ni-trogen arcs．JPhys D:Appl Phys，1978，11:1743 －1756．

［4］ 馬桂珍．真空電弧最新控制技術(shù)［J］．真空電子技術(shù)，1998，4:18－20．［5］ 袁海文，陳德桂，辜曉川．微型低壓限流器跑弧區(qū)磁場(chǎng)的計(jì)算及其對(duì)限流性能的影響［J］．西安交通大學(xué)學(xué)報(bào)，1996，30(3):19 －25．

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

［7］ 袁海文，李安云，李榮源．電弧通道磁場(chǎng)分布對(duì)斷弧效果的控制作用研究［J］．北京航空航天大學(xué)學(xué)報(bào)，2002，28(3):280 －282．

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

Artificial Research on Thunderstorm Electricity Photom ization Proctess

LONG Hao-ran，ZHU Li-bo

(Nanning Power Supply Bureau ，Nanning 530031，China)

Considered the photoionization in the gas discharge，Listed charge density equation，established Matlabmathematicalmodel，Simulation of a lightning discharge charge density changes，Current interpretation of the largestminesteep for the discharge of a time after a period of time before the peak in the discharge rather than the beginning．

photoionization;charge density;current steepness of the largestmine;peak

TM86

B

1004－289X(2013)03－0063－04

2013－02－25

龍浩然(1977－)，男，工學(xué)學(xué)士，南寧供電局，長(zhǎng)期從事配網(wǎng)維護(hù)檢修及配網(wǎng)規(guī)劃工作;

朱立波(1983－)，男，工學(xué)碩士，南寧供電局，長(zhǎng)期從事配網(wǎng)自動(dòng)化及檢修維護(hù)工作。