EMD與小波分析組合探測電離層VTEC周期方法研究

鮑亞東　劉長建　張學東　劉　宸　馮　緒

1　信息工程大學地理空間信息學院，鄭州市科學大道62號， 450001 2　65014部隊，沈陽市，110000

?

EMD與小波分析組合探測電離層VTEC周期方法研究

鮑亞東1劉長建1張學東2劉宸1馮緒1

1信息工程大學地理空間信息學院，鄭州市科學大道62號， 450001 265014部隊，沈陽市，110000

摘要：將EMD分解與小波分析組合應用于VTEC周期探測，通過對多尺度分解后的IMF分量進行小波分析，探測可能存在的電離層活動周期，從而提高探測成功率。分別對高、中、低3個不同緯度的VTEC序列進行周期探測，驗證了該方法的可行性，并發現可能存在6 d、24 d等周期。

關鍵詞：EMD；小波分析；VTEC；周期探測

電離層垂直總電子含量(vertical total electron content，VTEC)是描述電離層形態和結構的重要參量之一。通過對VTEC的周期性進行研究，不僅有助于深化對等電離層氣候與長期變化趨勢的認識，而且可以推進電離層結構及其變化、建模與預報等研究的發展，進而在一定程度上提高電離層結構分析的準確性和預報結果的精度[1-3]。

國際電離層組織(International GNSS Service，IGS)自1998年以來持續提供可靠的全球電離層格網產品(global ionosphere map，GIMs)，由于其數據覆蓋范圍大、時間分辨率適中，非常適合用于研究全球電離層的變化[1,4]。國內外許多學者利用IGS提供的VTEC數據對電離層的周期變化特征進行研究和分析。余濤等[1]通過對太陽活動高年(2000年)的全球VTEC數據進行傅里葉展開，分析了白天電離層TEC周年和半年變化的全球特征。對VTEC最大值出現季節進行統計顯示，全球大部分地區在春秋月份出現最大值，北半球近極地區在冬季出現最大值；韓吉德等[2]采用譜分析和小波分解的方法對2010年全球電離層VTEC的時空變化特征進行分析表明，高、中、低緯度地區VTEC具有周日和半周年變化現象，在二分點處存在峰值； Hernández-Pajares等[4]計算了1998～2008年的全球電離層總電子含量(global electron content，GEC)，并通過能譜分析發現了344 d、180 d、26 d等多個周期的存在；Oinats等[5]利用JPLG分析中心的VTEC產品對Irkutsk地區2003～2005年的VTEC變化與27 d周期的太陽活動進行了相關性分析，結果表明，其相關系數達到0.89。本文提出將經驗模態分解(empirical mode decomposition，EMD)和小波分析組合應用于VTEC周期探測，嘗試為電離層周期探測提供新方法。

1EMD與小波分析組合探測周期的方法

1.1EMD基本原理

EMD是一種自適應的分析和處理非平穩信號的方法。根據信號自身特性，在時域內從高頻到低頻將復雜的信號進行多尺度分解，得到若干個按頻率高低排列的IMF分量，并能很好地反映信號在任何時間的局部頻率特征[6-9]。EMD分解要求每個IMF分量滿足兩個條件[6,8]：1)零點數目與極值點數目相同或最多相差1；2)函數關于時間軸平均對稱。對信號x(t)進行EMD分解，過程如下：1)對x(t)的極大值和極小值序列分別進行三次樣條插值構造上下包絡曲線，并計算上下包絡線的均值m(t)；2)提取信號細節為d(t)=x(t)-m(t)；3)對信號細節序列d(t)重復操作步驟1)、2)，直到d(t)滿足IMF分量條件，記IMF(t)=d(t)；4)計算分解殘差r(t)=x(t)-IMF(t)；5)對殘差重復操作步驟1)～4)，即可得到多個IMF分量和一個殘差項。

1.2基于小波分析的周期探測理論

小波分析是在小波變換的基礎上通過對信號進行分析進而提取數字信息[10]。對于采樣時間間隔為常數Δt的離散時間序列x(t)，其小波變換形式為:

(1)

式中，a為時間尺度因子，b為平移因子，N為信號數量，ψa,b(t)為小波基函數φ(t)在參數a、b下的小波函數，Wf(a,b)為小波變換系數。對小波系數的平方進行積分，即可得到不同尺度a下的小波方差：

(2)

式(2)計算得到的小波方差能夠反映波動能量隨尺度的變化情況，通過查找小波方差的極大值點對應的尺度大小，即可確定信號序列可能存在的周期。

1.3EMD與小波分析組合探測思路

直接利用小波分析對原始VTEC序列進行周期探測時，由于序列中各種周期信息之間的相互影響，往往無法找到明顯的小波方差極大值點，從而導致周期探測失敗。通過EMD多尺度分解得到的各分量序列既包含原始序列的多重信息，又能將這些信息合理地分配到各分量序列中，從而提高各信息在分量序列中的可探測度。因此，本文提出將EMD與小波分析組合的方式對周期進行探測，以提高周期探測的成功率。其基本思路為：對VTEC序列進行多尺度EMD分解，得到包含序列信息的本征模態函數(intrinsic mode function，IMF)分量。通過對不同尺度下的IMF分量進行小波分析，找到不同分量中的小波方差極大值點，進而實現對VTEC可能存在的周期進行探測。

2實驗及分析

提取75°N、120°E處1998-07-01～2013-04-30每天(共5 418 d，當天 14:00LT)的VTEC值作為觀測序列進行研究。由圖1可以看出，在觀測時段內，VTEC主要存在劇烈變化期和平靜變化期兩個階段，以及較為明顯的上升區和下降區。若以10 TECu為臨界值(圖1虛線所示)將電離層活動分為高年和低年，則電離層活動高年(1999～2003年)序列變化范圍比較大，抖動性比較明顯；電離層活動低年(2004～2011年)電離層VTEC值及其變化范圍都較小，且每年存在明顯的峰值點和谷值點，峰值出現在每年5～7月份(夏季)，谷值出現在每年的11～1月份(冬季)，有一定的周期性變化規律。

根據本文提出的VTEC周期組合探測方法，對VTEC觀測序列值進行多尺度的EMD分解。考慮到周期探測不涉及序列的外延拓，為使上下包絡線能夠包含所有序列值的分解條件，在進行多尺度EMD分解時，將序列的端點值同時作為極大值和極小值插值點，與對應極值點分別對包絡曲線進行插值計算，并以分解效果不明顯為停止準則，對VTEC序列進行多級分解，最終得到1個殘差項和8個IMF分量。由圖2可以看出，EMD分解后得到的殘差曲線基本反映了VTEC的長期變化趨勢和電離層活動高、低年及上升、下降區的時間分布情況。

圖3為不同尺度EMD分解得到的8個IMF分量序列隨時間的變化情況。可以看出，隨著分解尺度的增大，對應IMF序列的變化頻率依次降低，體現出EMD分解可以得到頻率高低排列的本征模態函數的特點。對比IMF1～IMF6分量序列與圖1所示的觀測序列可以發現，電離層活動高年對應的IMF分量變化幅度較大，低年對應的IMF變化幅度較小。該現象表明，EMD分解得到的IMF分量能反映原信號的信息，猜測不同分量中可能隱含著VTEC的周期性信息。為了對該猜測的可能性進行判斷，依次對不同分量進行小波變換。Haar小波形式簡單、計算速度快，是最簡單的正交小波。為了提高計算效率，本文選取式(3)所示的Haar小波作為變換基小波，根據式(1)對各IMF分量進行小波變換：

(3)

根據式(2)計算得到各小波方差隨尺度因子的變化情況，見圖4(a)，各IMF分量的最大小波方差依次為2 d、6 d、12 d、24 d、151.5 d、134 d、179 d、366 d。圖4(b)為小波方差歸一化后的結果，歸一化準則定義為：

(4)

由IMF1～IMF8的歸一化小波方差最大值對應的尺度因子變化情況可以看出，隨著分解尺度的增大，所得到的IMF分量中反映的周期數量不斷增多。其中，IMF1～IMF3主要包含VTEC高頻變化，IMF4、IMF8所反映的24 d、366 d變化分別接近與太陽活動相關的27 d和1 a周期，IMF5、IMF6、IMF7與半年周期接近。由此可見，該方法可以實現對電離層VTEC周期的探測。

根據上述過程，對P1～P5共5個位置處的VTEC序列進行相同方法的周期探測，其坐標分別為(5°N, 120°E)、(40°N, 120°E)、(75°N, 120°E)、 (40°N, 60°W)、(40°N, 60°E)，結果見表1。可以看出，本文提出的組合探測方法在IMF2、IMF4、IMF7分量中的探測結果比較一致，變化周期分別為6 d、24 d和0.5 a；在IMF3分量中主要包含高、中緯度的10 d ～半個月周期和低緯度的22 d周期；IMF5分量主要體現高緯度處的151.5 d周期、中緯度處的1個半月周期和低緯度處的2個月周期變化特征；IMF6分量的周期探測結果與IMF7較接近；IMF8分量主要反映高、低緯度處的1 a變化周期和中緯度處1.5 a變化周期。根據本文方法的探測結果可以認為，全球電離層可能存在6 d、24 d和0.5 a的變化周期。此外，電離層活動在高緯度地區可能存在1/3 a的變化周期，中緯度地區可能存在半個月、1個半月和1.5 a的變化周期，低緯度地區可能存在2個月的變化周期。

3結語

本文提出將EMD和小波分析組合應用于電離層VTEC周期探測，在對VTEC序列進行EMD分解的基礎上，對各尺度下的IMF分量進行小波分析，進而對電離層周期進行探測。通過對高、中、低3個不同緯度處VTEC序列進行實驗，驗證了本文方法應用于VTEC周期探測的可行性，為電離層周期探測提供了新思路。同時發現了全球范圍內6 d、24 d的活動周期，以及高緯度地區1/3 a周期，中緯度地區半個月、1個半月和1.5 a周期，低緯度地區2個月周期存在的可能性。由于本文實驗的時間分辨率為1 d，因此并未對小于1 d的周期進行探測；同時也未能探測到大于1.5 a的周期，這可能與分解停止準則的設置和小波基函數的選擇有關，該問題將在下一步研究中解決。

參考文獻

[1]余濤, 萬衛星, 劉立波,等. 利用IGS數據分析全球TEC的周年和半年變化特性[J]. 地球物理學報, 2006, 49(4): 943-949(Yu Tao, Wan Weixing, Liu Libo, et al. Using IGS Data to Analysis the Global TEC Annual and Semiannual Variation[J]. Chinese J Geophys,2006, 49(4): 943-949)

[2]韓吉德, 王祖順, 王春青. 全球電離層時空變化特性分析[J]. 測繪地理信息, 2012, 37(6): 26-29(Han Jide, Wang Zushun, Wang Chunqing. Analysis of Temporal and Spatial Change in Global Ionosphere[J]. Journal of Geomatics, 2012, 37(6): 26-29)

[3]馬符訊, 黃張裕, 徐南. 一種電離層周期探測新方法[J]. 導航定位學報, 2014, 2(3): 30-32(Ma Fuxun, Huang Zhangyu, Xu Nan. A Method Based on the Combination of Wavelet Multi-Scale Decomposition and Fast Fourier Transform to Detect the Ionosphere Cycle[J]. Journal of Navigation and Positioning, 2014, 2(3): 30-32)

[4]Hernández-Pajares M, Juan J M, Sanz J, et al. The IGS VTEC Maps: A Reliable Source of Ionospheric Information since 1998[J]. J Geod, 2009, 83(3-4): 263-275

[5]Oinats A V, Ratovsky K G, Kotovich G V. Influence of the 27-Day Solar Flux Variations on the Ionosphere Parameters Measured at Irkutsk in 2003-2005[J]. Advances in Space Research, 2008, 42(4): 639-644

[6]Huang N E, Shen Z, Long S R, et al. The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Non-Stationary Time Series Analysis[J]. Mathematical Physical & Engineering Sciences, 1998, 454(1 971): 903-995

[7]李強, 吳健, 許正文, 等. 利用EMD方法提取太陽活動周期成分[J]. 空間科學學報, 2007, 27(1): 1-6(Li Qiang, Wu Jian, Xu Zhengwen, et al. Extracting of Periodic Components in the Solar Activity Using EMD Method[J].Chin J Space Sci, 2007, 27(1): 1-6)

[8]王堅, 高井祥, 王金嶺. 基于經驗模態分解的GPS基線解算模型[J]. 測繪學報, 2008, 37(1): 10-14(Wang Jian, Gao Jingxiang, Wang Jinling. GPS Baseline Solution Based on Empirical Mode Decomposition[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(1): 10-14)

[9]甘雨, 隋立芬, 王冰. 經驗模態分解閾值消噪方法及其在慣性導航系統數據處理中的應用[J].測繪學報, 2012,41(4): 504-509(Gan Yu, Sui Lifen, Wang Bing. EMD Threshold De-Noising and Its Applications in INS Data Processing[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4): 504-509)

[10]徐晨, 趙瑞珍, 甘小冰. 小波分析·應用計算[M]. 北京: 科學出版社, 2004(Xu Chen, Zhao Ruizhen, Gan Xiaobing. Wavelet Analysis·Applications Algorithm[M]. Beijing: Science Press, 2004)

Foundation support:National Natural Science Foundation of China,No. 41374041.

About the first author:BAO Yadong, postgraduate, majors in GNSS data processing and ionospheric modeling, E-mail：919668518@qq.com.

A Study on Method to Detect the Period of VTEC by Combining EMD with Wavelet Analysis

BAOYadong1LIUChangjian1ZHANGXuedong2LIUChen1FENGXu1

1School of Surveying and Mapping, Information Engineering University, 62 Kexue Road, Zhengzhou 450001,China 265014 Troops, Shenyang 110000,China

Abstract:The combination of EMD and wavelet analysis is used to detect the period of VTEC. The possible period of ionospheric activity is detected by analyzing the IMFs under multi-scale with wavelet analysis. The probability can be increased with this new method. The VTEC series at different latitudes released by IGS is used as an example. The feasibility and validity of the method promoted in this paper is proven. The possibility of 6 d, 24 d and others is discovered and a new approach is proposed for period detection of VTEC.

Key words:EMD; wavelet analysis; VTEC; period detection

收稿日期：2015-06-30

第一作者簡介：鮑亞東，碩士生，主要研究方向為GNSS數據處理與電離層建模，E-mail：919668518@qq.com。

DOI：10.14075/j.jgg.2016.06.003

文章編號：1671-5942(2016)06-0481-04

中圖分類號：P228

文獻標識碼：A

項目來源：國家自然科學基金(41374041)。