周蜜, 左芯睿, 王冬冬, 趙偉翰, 蔡力, 王建國

(1.武漢大學 電氣與自動化學院,湖北 武漢 430072; 2.雷電防護與接地技術(shù)教育部工程研究中心,湖北 武漢 430072; 3.中國艦船研究設(shè)計中心,湖北 武漢 430060)

0 引 言

雷電是雷暴天氣中發(fā)生的一種瞬時放電現(xiàn)象。一次完整的雷電放電過程往往包含十分豐富的雷電流分量信息,如Rakov等[1]報道了一次負地閃放電中電荷向地面轉(zhuǎn)移的3種模式:先導/回擊模式、連續(xù)電流模式和M分量模式。其中,回擊過程會導致大峰值電流、高峰值功率、強電磁輻射以及沖擊波等,這些物理效應(yīng)能對地面設(shè)施產(chǎn)生可觀的破壞作用[2-6]。另外,雷擊放電全過程中的轉(zhuǎn)移電荷量作為評估閃電放電能量或強度的核心參數(shù),是自然雷電引起森林火災和各類金屬接閃器、飛機機翼以及金屬屋頂?shù)炔考l(fā)熱、燒蝕和穿孔等災害的重要原因[6-10]。

雷電電流及其電荷量的測量對雷電防護領(lǐng)域的發(fā)展有重要的現(xiàn)實意義。然而,雷電的發(fā)生在時空上隨機程度很高,目前主要采用人工觸發(fā)閃電[10-13],或在高建構(gòu)物[14-15]上進行定點測量,這種通過直接測量方式收集到的可用數(shù)據(jù)有限,且效率較低。雷云對地放電形成的主通道不被云層遮擋,其輻射光的強度與放電通道中的電流大小密切相關(guān),光學觀測方法是諸多雷電相關(guān)研究中最重要的手段之一[16-17]。

Mitchard等[18-19]在實驗室內(nèi)通過半球形間隙自擊穿產(chǎn)生了閃電電弧光譜,報道了大電流和高電壓條件下產(chǎn)生的雷電弧光譜與自然閃電光譜的相似性,為聯(lián)系實驗室電弧放電與自然界雷電電弧放電提供了潛在的參考依據(jù)。在實驗室開展電弧放電實驗,模擬在波形特征上類似于雷電流的沖擊電流,可以有效把控試驗各階段的工況,相對容易地制造可重復的放電。綜合自然觀測與實驗室模擬試驗結(jié)果,有助于揭示閃電通道的電流空間分布和時間演化特征,提升對于雷電發(fā)展過程的科學認識。

目前針對實驗室放電觀測的各類影響因素已有一些研究[20-21],然而所觀測的電流脈沖仍集中于單一分量。如,利用實驗室放電產(chǎn)生具有較短波前時間和較長衰減時間的電流波形以模擬自然雷電回擊過程,發(fā)現(xiàn)光強僅在后續(xù)緩慢衰減階段與密切跟隨電流波形變化[22-23]。通過石墨間隙放電,產(chǎn)生具有緩慢上升前沿的衰減振蕩電流波形對自然雷電M分量電流進行模擬,發(fā)現(xiàn)光強峰值與電流峰值間顯著的線性關(guān)系[24-25]。有必要針對雷電全過程放電波形,綜合各類分量的波形特性,研究組合電流脈沖與放電通道光強度的關(guān)系。

本文通過在高壓實驗室搭建模擬閃電通道的人工傳輸線模型電路,研究類似自然雷電回擊電流、M分量、連續(xù)電流等不同階段波形疊加的參數(shù)調(diào)控方法。在此基礎(chǔ)上,利用構(gòu)建的實驗室電弧放電電流與通道光輻射強度同步觀測系統(tǒng),使用兩種光譜響應(yīng)不同的光電探測器,開展電弧放電電流與通道光輻射強度的同步觀測,分析不同波形下長間隙放電通道光輻射波形變化特征。最后,對通道電流與光輻射強度波形特征參數(shù)之間的關(guān)系進行探究,討論不同光譜帶寬的光電探測器觀測結(jié)果的差異。

1 試驗布置

閃電通道可以用圖1所示的分布電路模型來描述RLC傳輸線,其電壓V與電流I的電報方程如下:

圖1 閃電通道的分布電路模型

(1)

(2)

式中:R、L、C分別為單位長度串聯(lián)電阻、串聯(lián)電感和并聯(lián)電容;z′為閃電通道方位縱坐標;t是時間。通道的傳輸線參數(shù)會隨時間和空間發(fā)生變化,即傳輸線是非線性且不均勻的。通道電感隨攜帶z方向電流的通道核半徑的變化而變化,通道電阻與通道電子密度及通道核半徑有關(guān)。

利用集總參數(shù)鏈式網(wǎng)絡(luò),通過構(gòu)建如圖2所示的非均勻參數(shù)分布的人工傳輸線模型電路,將多個高壓脈沖電容(C1～Cn)和電感(L1～Ln)串聯(lián),形成L-C鏈。電容器充電后,經(jīng)觸發(fā)間隙(FS)放電,產(chǎn)生電流脈沖波形,各鏈輸出電流波形都是與前級鏈路波形疊加的結(jié)果。一個完整的雷電流多分量組合脈沖波形包含具有快速上升波頭和緩慢下降波尾的回擊電流,以及后續(xù)疊加有若干個M分量的連續(xù)電流?？傮w思路是,將首鏈(L1-C1鏈)設(shè)計為過阻尼狀態(tài),產(chǎn)生類似回擊電流的波形。通過調(diào)節(jié)中間鏈的電感值,可以改變疊加脈沖的數(shù)量及脈沖的振蕩幅度,制造疊加有模擬M分量電流脈沖的連續(xù)電流波形。

圖2 集總鏈式網(wǎng)絡(luò)電路原理圖

搭建了如圖3所示的模擬雷電多分量電流實驗平臺,輸出多分量組合電流波形,模擬雷電全過程的波形特征及發(fā)展趨勢。該平臺主要由模擬多分量電流發(fā)生器及放電通道同步觀測系統(tǒng)組成。使用電感-電容器組經(jīng)連接構(gòu)成L-C鏈,鏈數(shù)為10鏈。鏈數(shù)越多,形成的電流上升速度越快,同時電流持續(xù)的時間越長,本文實驗中調(diào)試發(fā)現(xiàn)鏈數(shù)為10是較為合適的。將每鏈電容器容量固定為6 μF,調(diào)整電感抽頭可改變每鏈電感值,可調(diào)節(jié)范圍為1～3.6 mH。調(diào)波電阻值設(shè)置為14.6 Ω,作為集中參數(shù)鏈式網(wǎng)絡(luò)電路中的負載電阻串入主回路,可以抑制電流反峰,因為自然界雷電流中除極為罕見的特殊案例外,絕大多數(shù)電流波形是單極性的。調(diào)波電阻越大,對電流反峰的抑制效果越明顯,但電阻過大將導致電流輸出峰值的降低。使用額定容量為15 kVA的充電變壓器,經(jīng)整流器向電容器充電。放電時,通過鎢銅半球電極觸發(fā),由一對安裝于絕緣放電架、直徑為10 mm的石墨棒電極放電,電極間空氣間隙長度為5 mm。

圖3 模擬雷電多分量電流實驗平臺

石墨棒放電間隙的通道光輻射強度由通道光輻射強度觀測系統(tǒng)進行測量,該系統(tǒng)由Thorlabs雪崩光電探測器、延長套管、中性密度濾光片、可調(diào)水平機械狹縫及相應(yīng)連接組件構(gòu)成。因雷電流通道發(fā)光光譜范圍主要集中在紅光敏感帶和近紅外敏感帶,采用了兩種不同型號的光電探測器,其中紅光敏感型光電探測器型號為APD120A2,其帶寬范圍由DC至50 MHz,響應(yīng)峰值為600 nm,可測量200～1 000 nm光譜范圍的光信號;近紅外敏感型光電探測器型號為APD120A,其帶寬范圍由DC至50 MHz,響應(yīng)峰值為800 nm,可測量400～1 000 nm光譜范圍的光信號。使用光學隔振平臺支撐光電探測器,調(diào)整光電探測器的機械狹縫中心對準石墨放電間隙中心,對光信號同步觀測并分析對比兩種型號光電探測器的輸出數(shù)據(jù)。

2 結(jié)果分析與討論

2.1 數(shù)據(jù)概況

通過調(diào)節(jié)每鏈電感值,產(chǎn)生了不同充電電壓下的多分量組合全過程電流波形,并同步觀測其光輻射強度脈沖波形,對波形及波形參數(shù)進行對比分析。試驗輸出的波形數(shù)據(jù)包含3個模擬回擊脈沖和15個模擬M分量脈沖。圖4是充電電壓設(shè)定為20 kV輸出模擬雷電流多分量組合波形,及利用兩種光電探測器同步測量的電弧通道光強度波形。

圖4 模擬雷電流多分量組合波形及同步通道光強度波形

對比Zhou等[26]記錄的火箭觸發(fā)閃電所有放電過程的通道底部電流波形(如圖5所示),試驗中輸出的組合電流幅值難以與真實雷電流分量幅值水平相匹配,但在波形形狀及發(fā)展趨勢上已能夠充分模擬。

圖5 觸發(fā)閃電通道底部電流波形(摘自Zhou等[26])

下文將紅光敏感型光電探測器輸出的光強波形命名為光強度A,近紅外敏感型光電探測器觀測到的波形命名為光強度B。

在模擬回擊脈沖階段,紅光敏感型光電探測器觀測到的光強波形A能更好地跟隨電流變化,而近紅外敏感型光電探測器觀測的波形B在下降沿產(chǎn)生了一定的畸變。對于此類畸變的產(chǎn)生,已有學者對其進行研究。Quick和Krider[27]在處理火箭觸發(fā)閃電的試驗數(shù)據(jù)中發(fā)現(xiàn),某些光輻射強度波形在初始上升階段后衰減比電流的更為迅速,而在后續(xù)衰減階段,光輻射強度波形出現(xiàn)“平臺”或“第二峰值”,電流則不存在類似的“平臺”或“第二峰值”,他們推測,這可能是由于初始高溫發(fā)展階段后,電弧通道后續(xù)冷卻,導致光譜朝波長更長的方向轉(zhuǎn)移。Carvalho等[28]總結(jié)了15次觸發(fā)閃電回擊過程的觀測結(jié)果,發(fā)現(xiàn)光輻射波形在首次峰值之后,緊接著出現(xiàn)第二次峰值,同樣將這一現(xiàn)象歸因于回擊電流到達峰值后通道冷卻過程中相關(guān)譜線強度得到增強。此外,根據(jù)文獻[20]在實驗室產(chǎn)生模擬回擊及連續(xù)電流電弧對此類畸變成因的研究,這一現(xiàn)象可以歸因于重復放電過程中所導致的電極表面不同程度的燒蝕以及電極形狀的差異引起的電弧形態(tài)的扭曲。

在其后的模擬連續(xù)電流階段及疊加M分量階段,兩光電探測器總體與電流波形的變化趨勢基本保持一致,但總體上圖4(c)中近紅外敏感型光電探測器輸出的結(jié)果相對圖4(b)紅光敏感型光電探測器與電流波形吻合得更好。

在圖4中對電流及兩種光電探測器測量的光強度波形脈沖進行了標記,將初始回擊脈沖標記為R1,后續(xù)疊加在連續(xù)電流過程中的M分量脈沖逐個標記為M1、M2、M3及M4,將連續(xù)電流階段基本結(jié)束后產(chǎn)生的M分量脈沖標記為M5。針對上述脈沖波形,對其特征參數(shù)(峰值、脈沖峰值、10%～90%上升時間及轉(zhuǎn)移電荷量)進行測算,研究模擬雷電流多分量組合波形脈沖及同步通道光強度波形脈沖間的關(guān)系。

2.2 幅值關(guān)聯(lián)

關(guān)于通道光強度幅值與雷電流幅值間的關(guān)系,國內(nèi)外研究者通過對自然閃電和人工觸發(fā)閃電進行觀測,如Diendorfer等[29]在對奧地利Gaisberg塔上行閃電初始階段電流的觀測和研究中發(fā)現(xiàn)電流幅值與相應(yīng)的光強度幅值呈線性相關(guān)。Gomes和Cooray[30]測量了250 mm和500 mm長火花間隙的輸出光信號和通道電流,發(fā)現(xiàn)在上升部分光學信號緊跟電流波形,且光強度上升時間與電流上升時間、光強度幅值與電流幅值均呈線性關(guān)系。利用國際空間站上599～900 nm波長的光度計,Kikuchi等[31]提取了11個閃電事件的可見光閃電數(shù)據(jù),分析表明負極性回擊絕對光強度與峰值電流呈明顯的正比關(guān)系。

借鑒上述研究進展,以及結(jié)合我們試驗得到的數(shù)據(jù)趨勢,針對產(chǎn)生的所有模擬回擊電流及模擬M分量電流及其對應(yīng)光強度的峰值與峰值脈沖進行關(guān)聯(lián)分析,分別繪制強制過零點線性擬合回歸線于圖6(a)和(b)。圖6(a)中可以看出兩種光電探測器觀測到的光輻射強度峰值與電流峰值具有一定線性相關(guān)性,對于圖6(b)所示的電流與光強峰值脈沖,同樣呈現(xiàn)較為強烈的線性相關(guān)性。其中,光強度B的峰值與電流峰值以及光強脈沖峰值與電流脈沖峰值的相關(guān)系數(shù)分別為0.979和0.994,而光強度A得到的數(shù)據(jù)進行線性擬合的結(jié)果相關(guān)性較弱,分別為0.842和0.982,進一步展現(xiàn)了兩類光電探測器觀測結(jié)果的差異。

圖6 模擬雷電流多分量組合波形及同步通道光強度波形峰值關(guān)系

Wang等[32]發(fā)現(xiàn)在雷電流發(fā)展的不同階段,由于波形特征與放電通道環(huán)境的差異,光強與電流之間的關(guān)系也存在差異,在圖6中同樣可以觀察到,相對于模擬M分量電流峰值和峰值脈沖,模擬回擊電流的峰值及峰值脈沖存在不規(guī)律偏移的現(xiàn)象,表明回擊電流和M分量電流與各自對應(yīng)的光強度應(yīng)遵循不一樣的規(guī)律。除去回擊數(shù)據(jù),僅對模擬M分量電流與光強峰值及峰值脈沖數(shù)據(jù)進行線性擬合分析,同樣強制回歸線過零點,如圖7(a)及圖7(b)所示。與圖6(a)相比,圖7(a)中光強度峰值與電流峰值的線性相關(guān)性更強,相關(guān)系數(shù)增大到更接近于1。但是,擬合回歸線斜率發(fā)生了明顯減小:光強度A從0.294變?yōu)?.175,減少了40%;光強度B則從0.342變?yōu)?.318,減少了7%。上述變化表明,即使針對同一電流水平,對于上升迅速、變化率更高的回擊電流,其通道輻射的光強度要比變化更慢的M分量的高,其中光強度A的這一效應(yīng)比光強度B的更加顯著。

圖7 模擬雷電流M分量波形及同步通道光強度波形峰值關(guān)系

計算M分量脈沖峰值時,電流需要減去背景連續(xù)電流水平,光強度也需要對應(yīng)減去背景光強水平。此時,光強度A的擬合線斜率從圖7(a)中的0.175增大為圖7(b)中的0.204,增幅為16%,同時線性相關(guān)性變差;光強度B從圖7(a)中的0.318增大為圖7(b)中的0.370,增幅為17%,也發(fā)生線性相關(guān)性變差的現(xiàn)象。上述結(jié)果表明,同樣針對同一電流水平,對于存在一定變化率的M分量電流,其通道輻射的光強度要比變化更慢甚至不變化連續(xù)電流的高,但是光強度A和光強度B的這一效應(yīng)相當。

上述結(jié)果進一步說明,對于不同分量的波形,其輻射產(chǎn)生的光強度與電流幅值的相關(guān)關(guān)系存在差異,這種差異與光強度探測儀器的光譜范圍有關(guān)。比較線性擬合相關(guān)系數(shù)數(shù)值大小,綜合考慮不同電流分量的線性斜率變化,可以推斷近紅外敏感型光電探測器觀測到的光強峰值(光強度B)與電流峰值的擬合結(jié)果更好。

2.3 波形時間關(guān)聯(lián)

Gomes和Cooray[30]通過在實驗室產(chǎn)生長間隙閃電,并同步觀測電流與光強波形,發(fā)現(xiàn)光強波形在上升階段能較好地跟隨電流波形,認為可以通過光學手段推算電流上升時間。通過繪制圖8所示的電流10%～90%上升時間與對應(yīng)光強度10%～90%上升時間線性擬合關(guān)系圖,線性擬合回歸線被強制過零點,兩種光強探測器輸出的擬合波形都呈現(xiàn)出很強的線性相關(guān)性,光強度上升時間均隨著電流上升時間增加呈正比例增長。光強度A和光強度B波形上升時間關(guān)系的線性擬合表達式分別為y=0.926x和y=0.938x,表明光強度上升時間總體上比對應(yīng)電流上升時間約短7%。由于模擬回擊脈沖波前時間較模擬M分量脈沖波前時間小2～3個數(shù)量級,近紅外敏感型光電探測器輸出波形中的畸變對擬合結(jié)果的影響可以忽略不計,所得兩回歸線的數(shù)據(jù)差異較小,回歸線斜率幾乎一致,均略小于1,表明光強波形上升時間較電流上升時間略短。

圖8 模擬雷電流多分量組合波形10%～90%上升時間與同步通道光強度10%～90%上升時間關(guān)系

2.4 轉(zhuǎn)移電荷量關(guān)聯(lián)

圖9為模擬雷電流多分量電荷轉(zhuǎn)移量與通道光強度-時間積分的線性擬合關(guān)系圖,其中轉(zhuǎn)移的電荷量在0.51～2.93 C之間。由于回擊脈沖持續(xù)時間較短,此處同樣忽略近紅外敏感型光電探測器觀測到的回擊波形畸變的影響。根據(jù)擬合結(jié)果可知,兩種光電探測器對應(yīng)光輻射強度-時間積分的相關(guān)系數(shù)較高,說明其對應(yīng)光輻射強度-時間積分與模擬連續(xù)電流轉(zhuǎn)移電荷量具有較強的線性關(guān)系。其中光強度B積分得到的擬合結(jié)果更好。試驗結(jié)果與文獻[33]報導的觸發(fā)閃電回擊在1 ms內(nèi)傳輸?shù)降孛娴碾娏鬓D(zhuǎn)移電荷量與相對光強度對時間積分之間存在較強的線性相關(guān)性的結(jié)論較為一致,為通過光學觀測反演雷電流轉(zhuǎn)移電荷量提供了依據(jù)。值得注意的是,文獻[33]中的電弧是自然界回擊電流電弧,而本文是在實驗室中產(chǎn)生的模擬雷電流多分量組合波形。

圖9 模擬雷電流多分量組合波形轉(zhuǎn)移電荷量與同步通道光強度-時間積分關(guān)系

3 結(jié) 論

通過搭建集中參數(shù)鏈式網(wǎng)絡(luò)電路,產(chǎn)生與波形特征與自然界雷電流全過程非常相似的模擬雷電流多分量組合脈沖,提供了產(chǎn)生模擬雷電流多分量波形的新思路。采用電流傳感器和兩種不同的光電探測器,同步測量了3組電流,共3個模擬回擊脈沖和15個模擬M分量脈沖與對應(yīng)通道光強度,開展了多分量組合脈沖電流與通道光強度特征參數(shù)關(guān)系研究,所得結(jié)論如下:

1)電流和兩種光電探測器輸出的光強度的峰值以及脈沖峰值間均存在較為明顯的線性關(guān)系;

2)光強度10%～90%上升時間與電流10%～90%上升時間之間也大致線性相關(guān);

3)電流轉(zhuǎn)移電荷量與光強度-時間積分間呈現(xiàn)較為明顯的線性相關(guān)性;

4)近紅外敏感型光電探測器觀測到的光強波形能更好地跟隨實測電流波形的變化,有望更準確地對放電轉(zhuǎn)移電荷量進行反演。