袁 萍,趙金山,王雪娟,穆亞麗,張 萌

(西北師范大學(xué)物理與電子工程學(xué)院，甘肅省原子分子物理與功能材料重點(diǎn)實(shí)驗(yàn)室，甘肅蘭州 730070)

?

利用Hα譜線Stark加寬研究閃電放電通道的電場和電流

袁 萍,趙金山,王雪娟,穆亞麗,張 萌

(西北師范大學(xué)物理與電子工程學(xué)院，甘肅省原子分子物理與功能材料重點(diǎn)實(shí)驗(yàn)室，甘肅蘭州 730070)

依據(jù)等離子體電場強(qiáng)度的診斷方法，得到了Hα線Stark加寬與閃電放電等離子體電場強(qiáng)度的關(guān)系；同時(shí)利用Borovsky定理對(duì)放電通道電場分布的分析及廣義歐姆定律，得出了閃電放電通道表面徑向電場與通道內(nèi)垂直電場及電流強(qiáng)度的關(guān)系，并分析了青海地區(qū)用無狹縫高速光柵攝譜儀獲得的一次云對(duì)地閃電回?fù)暨^程中光譜Hα線的Stark加寬，對(duì)通道不同位置的電場和電流進(jìn)行了定量計(jì)算，得到了該放電通道表面附近徑向方向電場為5.49×104～6.84×104kV·m-1，通道中電場垂直分量為3.29～4.10 kV·m-1，電流強(qiáng)度為28.98～36.04 kA.結(jié)果表明，通道表面附近的徑向電場遠(yuǎn)大于通道垂直電場，且隨著通道位置的升高，電流及沿通道垂直方向的電場強(qiáng)度均減小.

Stark加寬；電場強(qiáng)度；電流強(qiáng)度

閃電放電過程中，強(qiáng)電流產(chǎn)生的電磁輻射及熱效應(yīng)是導(dǎo)致電力和通訊設(shè)備中斷、工業(yè)控制系統(tǒng)損壞、引起森林火災(zāi)、造成人身傷害等各種雷電災(zāi)害的根源.對(duì)閃電電流和電場的研究是雷電防護(hù)領(lǐng)域所關(guān)注的熱點(diǎn)課題.傳統(tǒng)的方法主要是利用安裝在高塔上的感應(yīng)線圈直接測量電場和電流[1].人工引雷是直接獲得電流和電場觀測資料的主要途徑[2-4],但試驗(yàn)成本較高，并且火箭引雷的成功率受天氣條件等眾多因素的影響.人們也利用光學(xué)觀測得到通道發(fā)光強(qiáng)度強(qiáng)弱來定性地推斷放電電流[5-6].通過建立模型進(jìn)行模擬計(jì)算也是得到電流的一種途徑.其中，最常用的是傳輸線模型，利用閃電輻射電磁場強(qiáng)度估算電流.由回?fù)敉ǖ婪逯惦娏髋c峰值電場之間的關(guān)系，從閃電電場的觀測資料也可以反演得到峰值放電電流[7].現(xiàn)有的關(guān)于閃電電流的資料仍不能滿足相關(guān)研究的需要，這方面的工作對(duì)雷電的防護(hù)和閃電形成過程物理機(jī)制的研究都具有非常重要的意義.

閃電放電通道是典型的等離子體環(huán)境，利用對(duì)閃電光譜的觀測來研究閃電放電過程中的一系列物理特性是一種較為簡捷方便的途徑[8-9].目前，對(duì)地閃回?fù)暨^程的研究大多集中在閃電高速攝像和聲、光、電、磁的綜合觀測以及回?fù)裟Ｊ降陌l(fā)展方面，將光譜與閃電放電特征結(jié)合起來的分析相對(duì)較少[10-11].放電通道內(nèi)輻射原子在外電場作用下發(fā)生能級(jí)分裂，導(dǎo)致光譜線加寬，通過測量譜線寬度可得到放電通道電場大小與電流強(qiáng)度[12-13].

本文依據(jù)等離子體環(huán)境中的譜線Stark加寬理論研究放電通道的電場與Hα譜線Stark加寬的關(guān)系，并對(duì)使用無狹縫高速光柵攝譜儀測的1次云對(duì)地閃電回?fù)暨^程的光譜進(jìn)行分析，利用該光譜Hα線的Stark加寬計(jì)算了該閃電過程中通道的電場強(qiáng)度，并由閃電通道中電磁波傳播規(guī)律及歐姆定律計(jì)算得到了閃電通道中的電流強(qiáng)度，分析了電流和沿通道垂直方向電場分量隨通道高度的變化關(guān)系.

1 實(shí)驗(yàn)及理論方法

1.1 實(shí)驗(yàn)設(shè)備

本次拍攝閃電使用無狹縫高速光柵攝譜儀，記錄系統(tǒng)為美國Phantom(M310)公司的高速攝像機(jī)，試驗(yàn)的分光裝置為600條/mm的平面透射光柵，光柵置于攝像機(jī)物鏡前，閃耀波長500 nm，光譜分辨率為1.1 nm.

1.2 譜線Stark加寬與外電場的關(guān)系

閃電通道內(nèi)部是一個(gè)非常復(fù)雜的等離子體環(huán)境[14]，粒子間的相互作用主要是庫侖作用，可以看做粒子處于時(shí)變的電場，即[15]

(1)

(2)

即，原子在外電場中的能級(jí)分裂大小正比于電場強(qiáng)度大小[16-17].因此，氫原子在外電場中的哈密頓量可寫為[16]

(3)

其中后一項(xiàng)，即電場與氫原子電矩的作用可以看做微擾處理.在閃電通道這一等離子體環(huán)境中，由于輻射原子與每一個(gè)與其相互以庫侖力作用的粒子的距離都不相等，即相互作用的大小也不一致，譜線的分裂也不相等，疊加的結(jié)果即譜線的加寬[18].據(jù)準(zhǔn)靜態(tài)近似理論可知，電場隨時(shí)間的變化足夠慢[13]，則某頻率間隔內(nèi)的場強(qiáng)正比于干擾粒子存在的幾率[19]，可得到Stark加寬與電場強(qiáng)度之間的關(guān)系[20]為

(4)

其中，λ0為譜線的中心波長;ε0為真空中的介電常數(shù);m和e分別為單電子質(zhì)量和電荷量;c為真空中的光速;n1和n2分別為躍遷過程中上下能級(jí)的主量子數(shù).閃電通道中光譜輻射源在相應(yīng)時(shí)刻的電場強(qiáng)度為E.對(duì)于Hα線，n1=3,n2=2，上式化為

(5)

1.3 閃電通道的電場分布

根據(jù)Borovsky假定，閃電通道是一個(gè)直的等離子體柱，其中電場是均勻的，同時(shí)閃電通道中產(chǎn)生的電磁波按調(diào)諧函數(shù)eikz-iω t傳播[21-22].

根據(jù)假定，閃電通道為等離子體柱，則根據(jù)等離子體環(huán)境的歐姆定律[21]

(6)

(7)

在柱坐標(biāo)系(r,θ,z,t)中，電場為

(8)

其中，z方向?yàn)檠亻W電通道的軸向；r方向?yàn)殚W電通道的徑向.由假設(shè)閃電通道中的電場可寫為

(9)

代入(7)式得

(10)

由于Eθ與Er，Ez是去耦合的，即Eθ與Ez和Er的波模不同,即通過(5)式得到的電場強(qiáng)度是垂直電場和徑向電場的矢量和的大小.

由(7)式和(9)式可得

(11)

在閃電通道外部當(dāng)r→∞時(shí)，Er→0，且在閃電通道內(nèi)，Ez是有限的，在r=0時(shí)，Ez不為零，即上式的解為I0(γr)和修正的Bessel函數(shù)K0(γr):

(12)

(13)

將上式代入(10)式即可得

(14)

(15)

其中,a和b為單位電場下的復(fù)常數(shù);R為通道半徑;σ為閃電通道的電導(dǎo)率；ω為調(diào)諧函數(shù)的角頻率,是上升時(shí)間的倒數(shù).通常R的取值為0.1～10cm，本文中取R為1cm，σ為2.8×104S·m-1[23]，取回?fù)綦娏鞯牡湫蜕仙龝r(shí)間Δr為2μs，則可得

(16)

k為電磁波波數(shù)，則由色散方程[22]

(17)

可知k值約為5.33(1+i)×10-6cm-1.由文獻(xiàn)[23]可知，通道內(nèi)E2沿通道徑向方向變化極小且在通道半徑處連續(xù).由(11)和(14)式可得通道表面附近徑向電場與通道內(nèi)部的垂直電場之間的關(guān)系

(18)

將以上數(shù)據(jù)代入(17)式，取實(shí)部，得Er約為Ez的1.667×104倍，即通道表面附近的徑向電場Er遠(yuǎn)大于沿通道垂直方向電場Ez，即譜線Stark加寬的主要貢獻(xiàn)來自于通道表面附近電場的徑向分量，由(5)式得到的電場強(qiáng)度為通道表面附近電場的徑向分量Er.

1.4 閃電放電電流計(jì)算

由歐姆定律[23]可知，通道內(nèi)沿通道垂直方向的電場和電流關(guān)系為

(19)

由于Ez沿通道方向變化極小，忽略Ez的變化，上式積分得

(20)

通過上式與(5)式、(18)式可得閃電通道中譜線Stark加寬與電流強(qiáng)度的關(guān)系.

2 結(jié)果分析

選取用無狹縫光柵拍攝的一次云對(duì)地閃電回?fù)舫跗诘墓庾V，如圖1所示.圖2為通道某一位置處用相對(duì)強(qiáng)度表示的光譜,圖3為對(duì)沿通道不同位置的光譜中Hα譜線Stark加寬的Lorenz擬合，在此過程中選取譜線較為清晰，通道在一定范圍較直處的點(diǎn)，得到譜線Stark加寬在2.44～3.04 nm，且發(fā)現(xiàn)通道中位置離地面越高的點(diǎn)處譜線Stark加寬越小.

圖1 原始光譜Fig 1Original spectrum

圖2 某一高度處用相對(duì)強(qiáng)度表示的光譜Fig 2Spectrum at given position for discharge channel

圖3 通道不同位置處光譜Hα線的Lorenz擬合曲線Fig 3 Lorenz fitting for Hα line at different position

將圖3得到的不同位置處Hα譜線Stark加寬的數(shù)據(jù)代入(5),(17)與(19)式，得到相應(yīng)電場強(qiáng)度和電流強(qiáng)度數(shù)據(jù)，見表1所示.

表1 不同位置處Hα譜線Stark加寬與 閃電通道的電場和電流強(qiáng)度Tab 1 The width of Hα Stark broadening and electric field strength,current intensity at different position

由表1可以看出,沿通道垂直方向電場強(qiáng)度在3.29～4.10 kV·m-1，通道表面附近徑向方向電場強(qiáng)度在5.49×104～6.84×104kV·m-1，放電通道表面附近徑向電場遠(yuǎn)大于通道垂直電場，說明放電通道的帶電粒子主要集中在通道表面處，與Borovsky報(bào)道的結(jié)果基本吻合[22].由上表亦可以看出，電流強(qiáng)度在28.98～36.04 kA，與本小組以前報(bào)道的云對(duì)地閃電回?fù)敉ǖ缼缀纹骄挡顒e不大[25]；其他文獻(xiàn)[26]報(bào)道的負(fù)地閃中值電流大多約為30 kA，這說明我們的計(jì)算結(jié)果在合理的數(shù)值范圍之內(nèi).

從表1還可以看出，隨著通道中觀測點(diǎn)的增高，閃電放電通道的電場強(qiáng)度和電流強(qiáng)度均減小，Hα線Stark加寬亦隨之減小.其主要原因是回?fù)舻倪^程中，由于先導(dǎo)下行的速度比有效通道的擴(kuò)張速度快，導(dǎo)致有效通道半徑與通道高度呈反比，地面附近通道沒有完全展開，致使閃電通道下部的電子密度較大，粒子間的相互作用亦較大，即通道中的電場和電流都較大，譜線半寬也隨之增大.這一結(jié)論與其他理論模型得出的結(jié)論相一致[27].

用Hα譜線Stark加寬獲得閃電等離子體通道的電場強(qiáng)度和電流強(qiáng)度是一種有效的方法，但大氣中的氫原子的來源主要是水分子的電離，這與天氣的濕度有一定關(guān)系；且由于本實(shí)驗(yàn)在高原地區(qū)進(jìn)行，云底比較低，而且云層比較薄，經(jīng)過云體折射的雜散光照到光柵表面后，衍射形成無規(guī)則的彩色花斑，這些花斑構(gòu)成的背景會(huì)對(duì)閃電光譜的觀測帶來一定的影響；另外，在選取閃電通道中電磁波頻率即閃電上升時(shí)間和通道半徑時(shí)，取的是典型值，這也會(huì)對(duì)計(jì)算的結(jié)果造成影響，這就需要更加系統(tǒng)和嚴(yán)謹(jǐn)工作對(duì)數(shù)據(jù)和計(jì)算結(jié)果進(jìn)行修正和完善.

3 結(jié)論

利用Griem對(duì)等離子體中電場強(qiáng)度的診斷方法，得到Hα線Stark加寬與閃電放電等離子體電場的關(guān)系，并由Borovsky對(duì)放電通道電場分布的分析及廣義歐姆定律，得出閃電放電通道表面徑向電場與通道內(nèi)垂直電場及電流強(qiáng)度的關(guān)系,對(duì)用無狹縫光柵拍攝的一次云對(duì)地閃電回?fù)暨^程光譜進(jìn)行分析，計(jì)算得到放電通道表面附近徑向方向電場強(qiáng)度在5.49×104～6.84×104kV·m-1,通道內(nèi)垂直電場分量強(qiáng)度在3.29～4.10 kV·m-1，通道內(nèi)電流強(qiáng)度在28.98～36.04 kA，對(duì)通道不同位置處電流強(qiáng)度和附近徑向電場及通道內(nèi)垂直電場求得到的平均值分別為32.53 kA和6.17×104及3.70 kV·cm-1.實(shí)驗(yàn)結(jié)果表明，通道表面附近徑向電場遠(yuǎn)大于通道的垂直電場，且隨著高度的增高通道中沿通道垂直電場和放電電流都在增大，說明放電通道的電子主要集中在通道表面處，且通道下部的等離子體密度較大，粒子間的相互作用更加強(qiáng)烈.與以往只是能觀測到閃電通道外的輻射電場與大氣電場相比，本工作對(duì)閃電的形成機(jī)制的研究具有更重要的意義，對(duì)閃電放電通道電流強(qiáng)度的計(jì)算與觀測不僅在方法上提供了新的方向，也為雷電的防護(hù)與研究提供了依據(jù).

[1] ROMERE C,PAOLONE M,RUBINSTEIN M,et al.A system for the measurements of lightning currents at the Santis Tower[J].ElectrPowSystRes,2012,82:34.

[2] BIAGI C J,UMAN M A,HILL J D,et al.Observations of stepping mechanisms in a rocket-and-wire triggered lightning flash[J].JGeophysRes,2010,115:701.

[3] QIE X S,JIANG R B.Simultaneously measured current,luminosity and electric field pulses in a rocket-triggered lightning flash[J].JGeophysRes,2011,116:587.

[4] SCHOENE J,UMAN M A,RAKOV V A.Return stroke peak current versus charge transfer in rocket-triggered lightning[J].JGeophysRes,2010,115:951.

[5] HILL R D.Channel heating in return stroke lightning[J].JGeophysRes,1971,76:637.

[6] ADACHI T,CUMMER S A,LI J,et al.Estimating lightning current moment waveforms from satellite optical measurements[J].GeophysResLett,2009,36:475.

[7] LU G,CUMMER S A,BLAKESLEE R J,et al.Lightning morphology and impulse charge moment change of high peak current negative strokes[J].JGeophysRes,2012,117:907.

[8] 李忠文,袁萍,喬紅貞,等.利用Hα線的Stark加寬計(jì)算激光誘導(dǎo)Cu等離子體的電子密度[J].原子與分子物理學(xué)報(bào),2010,27(3):420.

[9] 屠昕，陸勝勇，嚴(yán)建華.大氣壓直流氬等離子體光譜診斷研究[J].光譜學(xué)與光譜分析,2006,26(10):1785.

[10] BARNARD J,COOPER J,SMITH E W,et al.The broadening of He I lines including ion dynamic corrections,with application toλ4471 ?[J].SpectroscRadiatTransfer,1974,14:1025.

[11] 林潔麗.高分辨(飽和)分子光譜譜線線型、線寬及其應(yīng)用的研究[D].武漢：中國科學(xué)院武漢物理與數(shù)學(xué)研究所，2000.

[12] GRIEM H R.Plasma Spectroscopy[M].NewYork:McGraw Hill,1964.

[13] SHEFER R E,BEKEFI G.Stark broadening induced by the intense electric field of a pulsed relativistic magnetron[J].PhysFluids,1979,22:1584.

[14] 馬騰才,胡希偉,陳銀華,等.離子體物理原理[M].合肥:中國科學(xué)技術(shù)大學(xué)出版社,1988.

[15] 孔策.離子體光譜學(xué)導(dǎo)論[M].北京：北京大學(xué)出版社,2013.

[16] 邱德仁.原子光譜分析[M].上海:復(fù)旦大學(xué)出版社,2001.

[17] 張昌莘.在均勻電場中氫原子光譜的分裂規(guī)律[J].安徽師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,28(4):407.

[18] GRIEM H R．SpectralLineBroadeningbyPlasmas[M].New York：Academic Press，1974.

[19] CZARNETZKI U,LUGGEH?LSCHER D,D?BELE H F.Sensitive electric field measurement by fluorescence-dip spectroscopy of rydberg states of atomic hydrogen[J].PhysRevLett,1998,81:4592.[20] GRIEM H R.Stark broadening of isolated spectral lines from heavy elements in a plasma[J].PhysRev,1962,128:515.

[21] 陳渭民.雷電學(xué)原理[M].北京：氣象出版社，2006.

[22] JOSEPH E,BOROVSKY.An electrodynamics description of lightning return strokes and dart leaders:Guided wave propagation along conducting cylindrical channels[J].JGeophysRes,1995,100(D2):2697.

[23] 王雪娟，袁萍，岑建勇，等.依據(jù)光譜研究閃電放電通道的半徑及能量傳輸特性[J].物理學(xué)報(bào),2013,62(10):4066.

[24] 劉文芳，李桂林.廣義歐姆定律[J].吉林師范大學(xué)學(xué)報(bào)(自然科學(xué)版),1986(2):48.

[25] 郭逸瀟,袁萍.NI 4935 ?譜線Stark加寬估算閃電放電等離子體電子密度[J].高原氣象,2009，26(5)：864.

[26] 孫景群.大氣電學(xué)基礎(chǔ)[M].北京：科學(xué)出版社,1987.

[27] BABA Y,ISHII M.Lightning return-stroke model incorporating current distortion[J].IEEETElectromagnC,2002,44:476.

(責(zé)任編輯 孫對(duì)兄)

Study on the electric field and current of lightning channel by Stark-broadening

YUAN Ping,ZHAO Jin-shan,WANG Xue-juan,MU Ya-li,ZHANG Meng

(Key Laboratory of Atomic and Molecular Physics & Functional Materials of Gansu Province,College of Physics and Electronic Engineering,Northwest Normal University,Lanzhou 730070,Gansu,China)

According to Griem’s method of diagnosis electric field in plasma,the relationships between the Stark-broadening of line and the electric field of lightning plasma are obtained.Meanwhile,using the analysis of electric field distribution in lightning discharge channel performed by Borovsky and the generalized Ohm’s law,the relationships between the radial electric field near the channel and the vertical electric field in the channel,and between the radial electric field near the channel and the current are derived.On the basis of above relationships and using the spectrum of lightning plasma,the electric field and the current at lightning different channel height are calculated.The results show that the radialelectric fields near the channel are in the range of 5.49×104～6.48×104kV·m-1;the vertical electric fields are in the range of 3.29～4.10 kV·m-1;the currents are in the range of 28.98～36.04 kA;the radial electric field near the channel is more stronger than that the vertical electric field;the current and the vertical electric field decrease with increasing channel height.

Stark-broadening;electric field strength;current intensity

10.16783/j.cnki.nwnuz.2016.01.008

2015-10-20;修改稿收到日期:2015-12-03

國家自然科學(xué)基金資助項(xiàng)目(11475139,11365019)，甘肅省科技計(jì)劃資助項(xiàng)目(1506RJZA119)

袁萍(1963—)，女，甘肅古浪人，教授，博士，博士研究生導(dǎo)師.主要研究方向?yàn)樵咏Y(jié)構(gòu)與雷電物理.

E-mail:yuanp@nwnu.edu.cn

P 427.3

A

1001-988Ⅹ(2016)01-0035-05