王彥輝， 閔應(yīng)昌， 劉亞櫟， 劉銀萍

(1.南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心， 南京 210044； 2.南京信息工程大學(xué)大氣物理學(xué)院， 南京 210044； 3.中國科學(xué)院西北生態(tài)環(huán)境資源研究院高原大氣物理研究室， 蘭州 730000)

現(xiàn)代信號處理技術(shù)和高速數(shù)據(jù)采集技術(shù)的應(yīng)用提升了雷電探測的能力和水平，閃電甚高頻(very high frequency, VHF)輻射源定位技術(shù)能夠?qū)崿F(xiàn)閃電通道發(fā)展的時空演變過程分析，對閃電物理機(jī)理研究及雷電電磁脈沖防護(hù)具有重要的意義[1-4]。

在整個閃電定位系統(tǒng)中最核心的技術(shù)是閃電定位算法[5]。定位過程主要包括波形預(yù)處理、尋找脈沖、匹配脈沖以及到達(dá)時差方法解算最優(yōu)解[6]。利用閃電輻射源的脈沖信號到達(dá)不同測站的到達(dá)時間差(time difference of arrival, TDOA)進(jìn)行閃電通道定位的關(guān)鍵就是準(zhǔn)確地找到輻射源脈沖，并將不同測站的同一個放電事件產(chǎn)生的脈沖匹配在一起，錯誤的脈沖匹配會得到不準(zhǔn)確甚至錯誤的定位結(jié)果[6]。因此，將閃電產(chǎn)生的相同輻射源脈沖準(zhǔn)確匹配獲得脈沖的到達(dá)時間差對定位結(jié)果至關(guān)重要。Shao等[7]將原始波形經(jīng)過Hilbert變換后獲取歸一化的功率波形，再從歸一化的功率波形中提取脈沖的峰值信息,并進(jìn)一步判斷兩兩脈沖峰的時間差是否小于兩站之間的光程差(距離/光速)，如果兩站之間的時間差小于光程差，那么說明這兩個脈沖的匹配是合理的。Shi等[8]在該方法基礎(chǔ)上做了改進(jìn)，在波形匹配后，首先，判斷兩兩脈沖峰的時間差是否小于兩站之間的光程差，其次，為保證兩個測站的脈沖匹配唯一性，其脈沖到該段波形最大幅值的高度差應(yīng)最接近?；谠疾ㄐ无D(zhuǎn)換的歸一化功率波形進(jìn)行尋峰運算，獲得滿足條件的脈沖峰值幅度和峰值時間；基于兩個子站間脈沖的單脈沖峰值時間差限制條件和幅度相近條件，獲得匹配的脈沖峰值時間，并結(jié)合子站位置信息，組成到達(dá)時間非線性方程組，使用最優(yōu)化算法獲得最優(yōu)解。該方法由于僅僅采用單脈沖峰值時間和幅度進(jìn)行匹配，在脈沖信號比較豐富的情況下，會出現(xiàn)較多的錯誤匹配，導(dǎo)致較低的定位計算效率；同時，采用單脈沖相似幅度的匹配，也會導(dǎo)致部分放電事件不能有效匹配的情況，需要進(jìn)一步改進(jìn)算法。

中外一些定位系統(tǒng)在利用閃電低頻/甚低頻(low frequency/very low frequency, LF/VLF)電場脈沖信號的定位研究中采用了互相關(guān)(cross-correlation)算法來實現(xiàn)閃電信號到達(dá)各個測站的時間估計[9-11]。

劉波等[12]將頻域時間反轉(zhuǎn)(time reversal, TR)聚焦算法應(yīng)用到閃電輻射源定位中，李書磊等[2]將正交傳播算子(orthogonal propagator method, OPM)方法應(yīng)用于閃電寬帶輻射源定位。

為了增加脈沖信號的匹配成功率和提高脈沖峰值時間的精度，從而提高定位精度和精細(xì)化程度，F(xiàn)an等[13]將經(jīng)驗?zāi)B(tài)分解(empirical mode decomposition, EMD)信號處理技術(shù)應(yīng)用到閃電波形處理中，實現(xiàn)多站閃電電場波形的低頻濾波和高頻噪聲降低，并基于處理后的波形，通過逐步縮減窗口進(jìn)行波形相關(guān)匹配后轉(zhuǎn)化形成匹配的功率波形，尋找脈沖獲得匹配的脈沖峰值時間，并采用到達(dá)時差方法解算最優(yōu)解。但由于EMD固有的模式混疊現(xiàn)象，采用EMD方法處理閃電電場信號后，利用波形互相關(guān)進(jìn)行脈沖匹配時仍存在多站波形匹配不正確，導(dǎo)致脈沖信息提取不準(zhǔn)確[14]。

基于上述原因，采用EMD方法對原始信號進(jìn)行分解后再部分合成，并對合成后的波形信號進(jìn)行歸一化處理，獲得較干凈和平穩(wěn)的閃電電場信號波形。提出了四脈沖序列的高度比例特征和間隔特征作為脈沖匹配條件，對相同閃電輻射源脈沖進(jìn)行匹配，并基于到達(dá)時間差對閃電輻射源進(jìn)行三維定位，描繪出閃電通道的發(fā)展路徑并計算出各個輻射源發(fā)生的時間。

1 閃電輻射源定位方法

利用到達(dá)時間差對閃電輻射源進(jìn)行定位，4個及4個以上測站的同一輻射源的到達(dá)時間差定位出一個輻射源的空間位置及其發(fā)生的時間，多個輻射源點能較好地描繪出閃電通道的基本形狀和放電通道的發(fā)展過程。對于同一個閃電輻射源，(xi,yi,zi)和ti表示第i個測站的三維坐標(biāo)和接收到輻射源信號的時間，(x,y,z)和t表示輻射源空間位置和發(fā)生時間，c為光速，則第i個測站接收到閃電輻射源的時間滿足關(guān)系式為

(1)

每個測站都可得到式(1)所示的一個關(guān)系式，4個及4個以上測站得到一個方程組，解方程組求得輻射源的空間位置(x,y,z)和輻射源發(fā)生的時間t。

2 數(shù)據(jù)預(yù)處理方法

經(jīng)驗?zāi)B(tài)分解(EMD)是一種新的處理非平穩(wěn)信號的方法[15]，能將復(fù)雜信號分解為有限個本征模函數(shù)(intrinsic mode function, IMF)。EMD適合于非線性、非平穩(wěn)信號的分析，也適合于線性、平穩(wěn)信號的分析。經(jīng)驗?zāi)B(tài)分解方法的實質(zhì)是通過特征時間尺度來識別信號中所內(nèi)含的所有振動模態(tài)，為了從原始信號中分解出本征模函數(shù)，經(jīng)驗?zāi)B(tài)分解方法過程如下：找到測站i的原始信號Xi(t)所有的極值點，用3次樣條曲線分別擬合出上下極值點的包絡(luò)線，原始信號減去兩條包絡(luò)線的平均值得到c1，根據(jù)預(yù)設(shè)判據(jù)判斷c1是否為IMF，如果不是，則以c1代替Xi(t)，重復(fù)以上步驟直到c1滿足判據(jù)，則c1就是需要提取的IMF分量，每得到一階IMF分量，就從原信號中扣除該分量，重復(fù)以上步驟，直到信號最后剩余部分rn滿足預(yù)設(shè)判據(jù)，剩余部分rn通常是直流分量、單調(diào)分量和低頻周期分量或3種分量的疊加。

利用EMD方法將原始信號Xi(t)分解成一系列IMF分量以及剩余部分rn的線性疊加，可以表述為

(2)

根據(jù)EMD的信號分解原理，EMD方法將組成原始信號的各尺度分量不斷從高頻到低頻進(jìn)行提取，其中殘差分量rn通常是單調(diào)序列或者常值序列或者低頻周期序列或三者的疊加，也稱為趨勢分量，利用EMD方法對原始信號進(jìn)行分解之后，再由去除趨勢分量之外的其他分量合成得到關(guān)系式

(3)

的信號，獲得較干凈的信號xi(t)，從而獲得非線性和非平穩(wěn)信號中的輻射源脈沖。

閃電輻射源的寬帶電場原始信號經(jīng)過EMD方法分解后再合成，由于不同的測站對信號的放大能力不完全相同，以及測站與輻射源的距離也會影響測站接收到的輻射源信號強度，因此需要對EMD方法處理之后的閃電電場信號進(jìn)行歸一化處理，以使得不同測站的同一輻射源脈沖信號強度更趨于一致，并根據(jù)波形對應(yīng)的GPS(global positioning system)時間實現(xiàn)信號波形的整體對齊。

3 脈沖匹配方法

本文中使用的數(shù)據(jù)對應(yīng)的采樣率為20 MS/s，單次閃電的數(shù)據(jù)記錄時長為1.2 s，即每個測站的一次閃電信號由2 400萬個數(shù)據(jù)點組成，相鄰兩個數(shù)據(jù)點的時間間隔為50 ns。

信號的局部極大值與其左右極小值的高度差為該脈沖相對高度，以預(yù)設(shè)的脈沖最小相對高度和脈沖最小間隔為條件提取各個測站的脈沖。以明德測站為主站，明德測站的脈沖序列為樣板。由于本文中使用的數(shù)據(jù)其探測站間最遠(yuǎn)距離不超過15 km，即相同輻射源脈沖到達(dá)各個測站的光程差在±0.05 ms范圍內(nèi)，本文以各個測站0.15 ms信號波段的GPS時間相同確定各測站滑動窗口的位置，滑動窗口的數(shù)據(jù)長度為3 000個數(shù)據(jù)點，窗口時間長度為0.15 ms，即在滑動窗口內(nèi)匹配相同輻射源的脈沖信號。各測站完成當(dāng)前窗口的脈沖匹配，滑動窗口移動到主站下一個脈沖所在位置，以主站的該脈沖時間為中心在各測站截取GPS時間相同的0.15 ms信號波段以獲得新的滑動窗口。

對于同一個輻射源的脈沖，歸一化處理之后各測站脈沖的幅度仍然存在差異，對所有輻射源脈沖又不能明確這種差異在某一特定的范圍。本文中采取的方法是：在各測站的滑動窗口內(nèi)，采用四個脈沖序列的高度比例特征和間隔特征作為脈沖匹配的條件，每一次匹配出四組相同輻射源脈沖。

如圖1(a)所示，主站的第j、k、m、n4個脈沖在主站的當(dāng)前滑動窗口內(nèi)；如圖1(b)所示，第i個測站的第ij、ik、im、in4個脈沖在第i個測站的當(dāng)前滑動窗口內(nèi)。主站的第j個脈沖高度為基線，第k個脈沖到基線的高度為hk，第m個脈沖到基線的高度為hm，第n個脈沖到基線的高度為hn；同樣的，第i個測站的第ij個脈沖高度為基線，第ik個脈沖到基線的高度為Hk，第im個脈沖到基線的高度為Hm，第in個脈沖到基線的高度為Hn，hk與Hk的比例為λk，hm與Hm的比例為λm，hn與Hn的比例為λn，當(dāng)λk、λm、λn的兩兩之差在一定的容差范圍內(nèi)，即滿足4個脈沖序列的高度比例特征；同時，主站的第j個脈沖與第k個脈沖的距離為xjk，第k個脈沖與第m個脈沖的距離為xkm，第m個脈沖與第n個脈沖的距離為xmn，第i個測站的第ij個脈沖與第ik個脈沖的距離為Xijik，第ik個脈沖與第im個脈沖的距離為Xikim，第im個脈沖與第in個脈沖的距離為Ximin，當(dāng)xjk與Xijik之差在一定容差范圍內(nèi)，xkm與Xikim之差在一定容差范圍內(nèi)，xmn與Ximin之差在一定容差范圍內(nèi)，即滿足4個脈沖序列的間隔特征。

圖1 脈沖分布示意圖Fig.1 Schematic diagram of pulses distribution

滿足上述的脈沖匹配條件，則第i個測站的第ij、ik、im、in4個脈沖分別與主站的第j、k、m、n4個脈沖匹配成功；下一個測站與主站在當(dāng)前滑動窗口內(nèi)進(jìn)行脈沖匹配，直到最后一個測站完成當(dāng)前窗口內(nèi)的脈沖匹配，滑動窗口向后滑動，進(jìn)行下一組脈沖的匹配，重復(fù)上述步驟，直到滑動窗口移動到信號末尾。

如圖1所示，當(dāng)hk、hm、hn、Hk、Hm、Hn均為非零時，設(shè)

(4)

設(shè)α為高度比例差的容差范圍，β為脈沖間隔差的容差范圍，本文中設(shè)定α=50%，β=0.75 μs。若λk、λm、λn滿足關(guān)系式

(5)

則滿足4個脈沖序列的高度比例特征；若xjk、xkm、xmn和Xijik、Xikim、Ximin滿足

(6)

則滿足4個脈沖序列的間隔特征。同時滿足4個脈沖序列的高度比例特征和間隔特征，則完成第i個測站的第ij、ik、im、in4個脈沖分別與主站的第j、k、m、n4個脈沖的匹配。當(dāng)這4個脈沖序列的某一組脈沖相對于第一組脈沖高度為零或趨近于零時，例如當(dāng)hm或Hm為零或者趨近于零時，則不計算該組脈沖的高度比例λm，對于該組脈沖只需滿足四脈沖序列的間隔特征即可認(rèn)為該組輻射源脈沖匹配完成。

本文中將利用一次探測到的閃電寬帶電場信號驗證所提出的多脈沖匹配方法的效果，用于進(jìn)一步補充說明本文的多脈沖匹配方法。

4 實驗分析

4.1 數(shù)據(jù)采集

閃電甚高頻(very high frequency,VHF)輻射源定位系統(tǒng)[3]主要由VHF天線、寬帶電場變化測量天線、帶通濾波器、對數(shù)放大器、高速A/D數(shù)據(jù)采集卡、高精度時鐘(時間精度50 ns)、處理器及數(shù)據(jù)無線傳輸模塊組成。VHF接收系統(tǒng)中心頻率270 MHz，帶寬為6 MHz，高速A/D數(shù)據(jù)采集卡的采樣率為20 MS/s，寬帶電場系統(tǒng)帶寬0～10 M。為把接收到的閃電寬帶電場信號波形不失真的放大，該系統(tǒng)采用了對數(shù)功放電路。閃電輻射的脈沖信號幅度峰值和峰值時間由數(shù)字化模塊記錄緩存，記錄長度根據(jù)基線長度及處理數(shù)據(jù)的分辨率而定，本文采用的單次數(shù)據(jù)記錄時長為1.2 s,系統(tǒng)觸發(fā)及數(shù)據(jù)記錄的時間與GPS高精度時鐘同步。信號的噪聲電平由同步觸發(fā)閾值電路控制，處理器接收、處理和記錄輻射源脈沖峰值絕對時間信息，并實時通過無線數(shù)據(jù)傳輸系統(tǒng)將時間信息發(fā)送到中心站進(jìn)行時差計算及實時處理顯示。

采用和分析的數(shù)據(jù)來自于在青海省大通縣進(jìn)行的雷電綜合觀測實驗，閃電發(fā)生在2011年8月15日0：2：41。大通縣位于青海省東北部，平均海拔約2 600 m，是一個多山地區(qū)，地形復(fù)雜、氣候垂直差異明顯，由于地處內(nèi)陸高原地區(qū)，海拔較高，受高原天氣系統(tǒng)與西風(fēng)帶天氣系統(tǒng)相互作用，加之河谷地形和下墊面性質(zhì)的復(fù)雜性，造成短時強降水較多，并常伴有雷暴和冰雹等強對流天氣。

如圖2所示，雷電綜合觀測實驗由7個觀測站組成，其中以明德站為中心，各測站分布半徑約為8 km。基于時差法的閃電VHF輻射源三維定位系統(tǒng)中,觀測試驗的7個測站的寬帶電場變化測量系統(tǒng)通過無線寬帶通信系統(tǒng)連接組成同步觀測網(wǎng)絡(luò)，采用主站控制方式，并用GPS同步的高精度時鐘(50 ns)記錄觸發(fā)時間，從各測站記錄的寬帶電場變化波形中提取并匹配相同輻射源的脈沖，計算各測站相同輻射源的脈沖信號到達(dá)時間并利用時差法定位技術(shù)實現(xiàn)閃電輻射源的三維定位，精確描繪閃電放電通道，其中，明德、藥草、新莊、苗圃、極樂、斜溝、良教這7個測站的寬帶電場變化測量系統(tǒng)數(shù)據(jù)采樣率為20 MS/s,單次閃電數(shù)據(jù)記錄時長為1.2 s。

根據(jù)圖2所示各測站的經(jīng)緯度和海拔高度，以明德站為中心站，確定各測站相對于中心站的坐標(biāo)，進(jìn)行閃電輻射源定位的輻射源空間位置也是以明德站為三維坐標(biāo)原點，各個測站的相對坐標(biāo)分別為明德(0，0，0)，藥草(3 178.44，2 921.02，131.538)，新莊(-3 355.01，5 572.89，64.680 7)，苗圃(-6 013.87，2 442.51，44.694 9)，極樂(-8 325.07，-2 473.84，196.083)，斜溝(-4 165.98，-4 407.31，242.115)，良教(883.798，-4 003.41，-9.318 57)，其中，各測站的相對坐標(biāo)以米為單位。

圖2 青海省大通縣測站地理位置Fig.2 Location of stations in Datong County, Qinghai Province

4.2 數(shù)據(jù)預(yù)處理結(jié)果

如圖3(a)所示，7個測站單次閃電的數(shù)據(jù)記錄時長為1.2 s，數(shù)據(jù)采樣率為20 MS/s，其中，明德站的觸發(fā)前數(shù)據(jù)占比為28%，藥草、新莊、苗圃、良教4個測站的觸發(fā)前數(shù)據(jù)占比為30%，極樂、斜溝兩個測站的觸發(fā)前數(shù)據(jù)占比為35%。

使用EMD方法對原始閃電數(shù)據(jù)分解再部分合成，從原始信號中提取出圖3(b)所示的信號分量，去除了圖3(c)所示的信號分量，去除的信號分量包括直流分量、單調(diào)分量和低頻周期分量。每個測站原始信號經(jīng)過EMD方法分解再部分合成，再對合成的信號進(jìn)行歸一化處理，得到的信號最大值為1,最小值為-1，并根據(jù)閃電信號波形對應(yīng)的GPS時間完成信號波形的整體對齊，如圖3(d)所示。

圖3 數(shù)據(jù)預(yù)處理結(jié)果Fig.3 Results of data preprocessing

4.3 輻射源脈沖匹配結(jié)果

如圖4所示為該次閃電的其中一個四脈沖序列匹配過程，其中，圖4(a)所示為7個閃電探測站在相同GPS時間采集到的閃電電場波形；圖4(b)是圖4(a)中標(biāo)出信號的幅度放大波形，明德(主站)的四脈沖序列與藥草(子站)的四脈沖序列進(jìn)行匹配，分別用紅虛線和藍(lán)虛線標(biāo)出明德測站與藥草測站的四脈沖序列；圖4(c)為該四脈沖序列的最終匹配結(jié)果，可以看到，7個測站中有6個測站實現(xiàn)了脈沖匹配，實現(xiàn)6站同步，可以計算出相應(yīng)的閃電輻射源位置。

圖4 閃電四脈沖序列匹配過程Fig.4 Four pulses sequence matching process of lightning

如圖5(a)所示為稀疏的大脈沖的匹配結(jié)果，如圖5(b)所示為密集的小脈沖的匹配結(jié)果，對于稀疏的大脈沖和密集的小脈沖，使用本文方法能獲得較好的脈沖匹配效果。以明德站的脈沖序列為樣板，根據(jù)滑動窗口內(nèi)4個脈沖序列的高度比例特征和間隔特征，一次匹配出四組相同輻射源的脈沖，滑動窗口再次向后移動到主站的下一個脈沖所在位置，完成其他輻射源脈沖的匹配，并計算出各個測站相同輻射源脈沖的到達(dá)時間。如圖5(c)所示為第884組輻射源脈沖匹配時滑動窗口的信號波形，可以看到，各個測站的脈沖特征差異較大，多個輻射源脈沖信號在時間上交錯，可能是由于閃電的多個分支通道同時發(fā)展，使得輻射源天線在較短時間內(nèi)同時接收來自多個閃電通道的輻射源脈沖信號。如圖5(d)所示為圖5(c)所示滑動窗口內(nèi)的脈沖匹配結(jié)果，其中極樂測站未匹配出輻射源脈沖的原因是該處接收到的輻射源信號較弱，干擾噪聲大于輻射源脈沖信號，而良教測站未匹配出輻射源脈沖的原因是四脈沖序列中存在可能其他分支通道的輻射源脈沖，且不同分支通道相距較遠(yuǎn)。對于2011年8月15日0:2:41的該次閃電共匹配出1 897組相同輻射源脈沖。

4.4 定位結(jié)果對比與閃電過程分析

如圖6所示，圖6(a)和圖6(b)分別是本文中多脈沖匹配方法與傳統(tǒng)的單脈沖匹配方法的定位結(jié)果。圖6中閃電VHF輻射源3D定位系統(tǒng)的顯示界面分為5個視圖，從上到下、從左到右依次為：輻射源高度隨時間變化圖、東西方向的垂直投影圖、輻射源數(shù)量隨高度變化圖、水平投影圖、南北方向的垂直投影圖。其中，輻射源數(shù)量隨高度變化圖的縱坐標(biāo)表示輻射源高度，橫坐標(biāo)表示相應(yīng)高度的輻射源數(shù)量，從圖6可以看出閃電輻射源的高度分布特征；垂直投影和水平投影的坐標(biāo)是以明德測站為原點，表示輻射源相對于明德測站的三維位置。本文中的閃電信號預(yù)處理方法及脈沖匹配方法實現(xiàn)了較多相同輻射源的脈沖匹配，利用輻射源脈沖的到達(dá)時間差定位方法計算輻射源三維坐標(biāo)及其發(fā)生時間，獲得了較清晰的閃電發(fā)展通道，定位出了一次閃電的多個放電發(fā)展過程，其中一些放電過程是單脈沖匹配方法未定位出的，定位結(jié)果如圖6(a)所示，根據(jù)顏色區(qū)分不同放電過程發(fā)生的時間及通道長度。

單脈沖匹配方法未將原始信號中的直流分量、單調(diào)分量和低頻周期分量去除，并且是基于多個測站的單脈沖匹配，依據(jù)脈沖到達(dá)時的GPS時間最接近匹配相同輻射源的脈沖，其定位結(jié)果如圖6(b)所示，部分通道定位效果不明顯，甚至不能定位出閃電的部分放電通道。從整體定位效果來看，本文提出的輻射源脈沖匹配定位方法效果更好。

如圖6(a)所示，閃電起始區(qū)域的東南角的閃電通道的輻射源數(shù)目較少，該時段的閃電信號波形如圖5(c)所示，部分測站接收到的輻射源脈沖信號較弱，小于該處的噪聲干擾，或是四脈沖序列中存在可能其他分支通道的輻射源脈沖，不同放電通道距離較遠(yuǎn)，且輻射源信號天線同時接收到來自多個放電通道的輻射源脈沖信號，使得多個放電通道的多個輻射源脈沖信號在時間上交錯，在這種情況下四脈沖序列同時匹配的匹配效果不夠理想，而傳統(tǒng)的單脈沖匹配方法則具有一些優(yōu)勢。

圖5 脈沖匹配結(jié)果Fig.5 Results of pulses matching

閃電發(fā)生于2011年8月15日0：2：41，利用相同輻射源到達(dá)時間差對閃電輻射源進(jìn)行定位，顏色顯示先后順序為藍(lán)、青、綠、黃、紅、洋紅，顏色代表閃電輻射源發(fā)生的時間先后，如圖6(a)所示。閃電起始位置在明德站的西北方向，相對坐標(biāo)大概(-2 km,8 km)處，第一個閃電輻射源發(fā)生在41.715 000 s。從41.840 s開始，閃電放電通道從閃電起始區(qū)域向西南方向發(fā)展，該通道在41.939 s結(jié)束，通道發(fā)展時間約100 ms，通道長度約7.8 km，通道發(fā)展速度約7.8×104m/s。該通道從42.169 35 s開始出現(xiàn)第一次再放電過程，通道長度約6 km，該次再放電過程發(fā)展較快，持續(xù)時間約7 ms，速度約8.6×105m/s，與曹冬杰等[16]的研究結(jié)果在同一個量級。

圖6 閃電輻射源定位結(jié)果對比Fig.6 Results of lightning radiation localization

從閃電起始區(qū)域向西北方向發(fā)展的閃電通道中輻射源脈沖數(shù)量較少。從42.326 6 s開始該通道再次放電，從輻射源定位結(jié)果來看，該通道再放電過程是沿著原來的通道快速向閃電起始區(qū)域發(fā)展，其發(fā)展方向與該通道第一次放電過程的發(fā)展方向相反，其通道長度約3.6 km，持續(xù)時間約3.5 ms，速度約為1×106m/s，曹冬杰等[16]對云閃的反沖流光速度估計為106～107m/s，與這里的結(jié)果基本一致。

在閃電起始區(qū)域西北方向的通道再次放電結(jié)束后間隔約1.58 ms，從閃電起始區(qū)域向西南方向發(fā)展的閃電通道出現(xiàn)第二次再放電過程，沿著第一次放電的通道向西南方向發(fā)展，該通道第二次再放電過程的放電通道長度約4 km，持續(xù)時間約16 ms，通道發(fā)展二維速度約2.5×105m/s，與Liu等[17]和張泉等[18]的研究結(jié)果基本一致。最后一個閃電輻射源發(fā)生在42.347 675 s，整個閃電過程持續(xù)時間約632.675 ms。

5 結(jié)論

采用經(jīng)驗?zāi)B(tài)分解(EMD)方法對原始信號進(jìn)行分解后再部分合成，并對合成后的波形信號進(jìn)行歸一化處理，獲得較干凈和平穩(wěn)的閃電電場信號波形。在采用EMD方法處理閃電寬帶電場信號之后，不采用波形互相關(guān)算法匹配脈沖，而是在單脈沖峰值時間差限制條件和幅度相近條件匹配閃電輻射脈沖的基礎(chǔ)上，提出了四脈沖序列的高度比例特征和間隔特征作為脈沖匹配條件，對相同閃電輻射源脈沖進(jìn)行匹配。并且，利用一次探測到的閃電寬帶電場信號，以本文方法與單脈沖峰值時間差限制條件和幅度相近條件匹配閃電輻射脈沖的方法進(jìn)行了脈沖匹配效果的比較，所提出的多脈沖匹配方法具有較好的自適應(yīng)性和準(zhǔn)確性，實現(xiàn)了更多組相同輻射源的脈沖匹配。對于本文中使用的這一次閃電數(shù)據(jù)，本文中所提出的多脈沖匹配方法最終定位出了1 897個閃電輻射源，而單脈沖匹配方法最終定位出了727個閃電輻射源。

通過上述閃電輻射源多脈沖匹配定位方法的提出和驗證，得到以下結(jié)論。

(1)通過EMD方法對閃電輻射信號進(jìn)行分解后再部分合成，去除了原始信號中的單調(diào)分量、直流分量和低頻周期分量，可以使信號波形更干凈和平穩(wěn)，便于后續(xù)的脈沖提取和匹配。

(2)對分解再合成的信號進(jìn)行歸一化處理，可以使不同測站的同一輻射源信號強度更趨于一致。

(3)以滑動窗口內(nèi)4個脈沖序列的高度比例特征和間隔特征作為各測站間相同輻射源的脈沖匹配條件，每一組脈沖匹配時四脈沖序列的高度比例特征和間隔特征不同，該方法具有較好的自適應(yīng)性和準(zhǔn)確性，能實現(xiàn)較多組相同輻射源的脈沖匹配。

(4)從整體定位效果來看，本文中提出的閃電輻射源多脈沖匹配定位方法效果較好，獲得了較清晰的閃電發(fā)展通道，定位出了單脈沖匹配定位方法未定位出的一次閃電的多個放電發(fā)展過程及閃電的多個放電通道形狀。

閃電輻射信號的預(yù)處理及同一輻射源的脈沖匹配不僅對閃電三維定位具有重要意義，也能夠幫助研究同一個輻射源發(fā)出的電磁脈沖到達(dá)不同測站的幅度和寬度差異，對雷電電磁脈沖防護(hù)具有一定的意義。在閃電發(fā)展過程中，多個雷電分支通道同時發(fā)展使閃電輻射源天線接收信號發(fā)生錯亂，滑動窗口內(nèi)不同分支通道的輻射源脈沖雜糅在一起，單個輻射源脈沖幅度較弱或者多個輻射源脈沖信號在時間上交錯，這種情況下的相同輻射源脈沖匹配效果不夠理想，該問題值得進(jìn)一步研究，這將是下一步重點研究的內(nèi)容。