一種電離層周期探測新方法

馬符訊，黃張裕，徐 南

（河海大學 地球科學與工程學院，南京 210098）

1 引言

近年來，電離層時空變化的研究已經成為了研究熱點之一，其影響因素包括太陽活動、地磁活動、日食和地震等［1－3］。本文以電離層時間變化規律中的太陽活動和電離層周期的探測為切入點，對電離層的時間周期性進行研究。掌握電離層活動規律有助于電離層的預報、精確模型的建立；有助于提高全球衛星導航系統（global navigation satellite system，GNSS）的定位精度，進而完善利用GNSS來監測電離層活動的理論和方法；有助于和無線電傳播密切相關的科學技術的研究和應用的發展［4－7］。

己有的太陽活動周期研究方法大多是基于傅立葉變換，將時間序列數據從時域變換到頻域進行分析，在頻域里得到的幅度峰值表現出來，來研究太陽活動在不同時間尺度上的變化，但是這并不能得到完全分離的某個時間尺度上的變化分量和探測出較弱的周期。本文采用的是一種基于小波多尺度分解和快速傅里葉變換相結合的方法來探測太陽活動和電離層周期，小波多尺度分解能將數據按波長的不同分解為不同波長范圍的數據，并根據需求剔除一部分數據從而得到一定范圍波長的數據，具有良好的局域分析能力［8－10］。快速傅里葉變換以其快速運算和信噪比閾值低等特點在頻譜分析領域得到了廣泛應用［11－12］。故本文將兩種方法相結合用于探測太陽活動和電離層周期，通過實驗表明兩者結合探測效果好、速度快、準確率高，可以探測出較弱的周期。

2 基于小波多尺度分解和快速傅里葉變換相結合的方法

小波變換即是用小波序列分解信號，把一個基本的小波經伸縮和平移后得到的小波序列，這種伸縮和平移是在乎母小波函數中同時進行的。母小波在尺度的選取和小波的變換中發揮重要的作用，要充分考慮其適用性，即小波函數具有多樣性，采用不同的母小波分析同一個問題會產生不同的結果。目前大多是通過分析的結果來判定母小波的選取的好與壞，不斷試驗來選擇母小波。莫萊（Morlet）小波定義為

式中，w0是小波數，經過多次試驗這里設為7。對于時間連續的X和xn值，在時間指數n的情況下，小波變換的Wn（s）是小波函數和最初連續時間的卷積。

式中，s表示尺度的大小，δt代表時間間隔，＊表示復數共軛元數，在Morlet小波情況下，小波尺度總等于傅里葉周期［13］。

本文提出的基于小波多尺度分解和快速傅里葉變換相結合的探測周期方法，首先對太陽黑子數原始數據序列進行小波多尺度分解預處理，再根據各層的重構的細節信號，然后對其進行快速傅里葉變換，生成各層細節信號的頻譜圖，再對此提取振幅最大的頻率，最終獲取原始數據序列的周期信息。

其特征在于包括以下步驟：

（1）對太陽黑子數原始數據序列進行小波多尺度分解預處理。

（2）根據各層的重構的細節信號，然后對其進行快速傅里葉變換，生成各層細節信號的頻譜圖像。

圖1 一種基于小波多尺度分解和快速傅里葉變換相結合的探測周期方法的具體算法流程

3）在上述步驟（2）中所得的頻譜圖像中，提取振幅最大的頻率，最終獲取原始數據序列的周期信息。

3 實例分析

本文利用基于小波多尺度分解和快速傅里葉變換相結合的探測周期方法對1988～2012年的總電子含量（total electron content，TEC），太陽黑子數（sunspot number，SSN）及沃爾夫黑子數（F10.7）數據序列進行了分析，探測三者周期，分析了TEC、SSN和F10.7指數三者的相關性。本文利用比利時皇家天文臺的太陽耀斑數據索引中心（sunspot index data center，SIDC）和全球衛星導航系統國際服務協會（international global navigation satellite system service，IGS）提供的數據作為數據源，用于探測太陽活動和電離層的周期，太陽活動參數包括SSN和F10.7指數，電離層參數為TEC。

表1 數據源及其相關信息

利用該方法分析164a的太陽黑子數序列，探測出太陽活動周期如表2所示。

表2 探測出的太陽周期

由表2可知，振幅出現2個峰值，分別對應27d和10.95a的周期。說明這2個周期為太陽活動顯著周期。下面對此法進行推廣，利用該方法對1988～2012年的TECSSN和F10.7數據序列進行了分析。

表3 1988～2012年的TEC、SSN和F10.7數據序列分析結果

從表3中我們可以看出TEC、SSN和F10.7三者的周期具有較好的一致性，下面進一步分析TEC、SSN和F10.7指數三者的相關性，根據下面的相關系數公式求出三者中兩兩的相關系數。

表4 TEC、SSN和F10.7指數三者的相關性

從表4中，我們可以看出，3個相關系數全部大于0.7，故三者高度線性相關；作為影響電離層變化的因素，F10.7和太陽黑子數可以相互替代；在電離層模型中引入參數F10.7更為合理。

4 結論

電離層時空變化存在一定的周期性，本文提出一種基于小波多尺度分解和快速傅里葉變換相結合的探測周期方法，對太陽活動和電離層周期進行探測，由此總結電離層時間變化規律，可以得到以下三點結論：

1）太陽黑子數出現的周期有：41a、11a、8.2a、5.5a、3.8a、2a、1a、216d、97d、38d、27d、12d、5.8d和2.5d。

2）F10.7、電離層TEC和太陽黑子具有相似的周期，且都具有11a和27d的顯著周期。

3）F10.7、電離層TEC和太陽黑子三者具有較強的相關性。

綜上，本文采用一種基于小波多尺度分解和快速傅里葉變換相結合的探測周期方法，用于探測太陽活動和電離層周期，通過實驗表明兩者結合探測效果好、速度快、準確率高，可以探測出較弱的周期，有利于電離層時空變化更進一步的研究。

［1］ ZHU Fu－ying，WU Yun，ZHOU Yi－yan，et al.A Statistical Investigation of Pre－earthquake Ionospheric TEC Anomalies［J］.Geodesy and Geodynamics，2011，2（1）：61－65.

［2］ ZHU Fu－ying，WU Yun.Anomalous Variations in Ionospheric TEC Prior to the 2011Japan Ms9.0Earthquake［J］.Geodesy and Geodynamics，2011，2（3）：8－11.

［3］ 姚宜斌，陳鵬，吳寒，等.2011年3月11日日本地震震前電離層異常變化分析［J］.科學通報，2012，57（5）：355－365.

［4］ 林劍.GPS電離層觀測技術及其在地震中的應用研究［D］.武漢：武漢大學，2011.

［5］ 柳景斌，王澤民，章紅平，等.幾種地基GPS區域電離層TEC建模方法的比較及其一致性研究［J］.武漢大學學報：信息科學版，2008，33（5）：479－483.

［6］ 李強，吳健，許正文，等.利用EMD方法提取太陽活動周期成分［J］.空間科學學報，2007，27（1）：1－6.

［7］ 黃磊，邵雪梅.青海德令哈地區近400年來的降水量變化與太陽活動［J］.第四紀研究，2005，25（2）：184－192.

［8］ 劉少明，申重陽，孫少安，等.小波多尺度分解特征分析［J］.大地測量與地球動力學，2004，24（2）：34－41.

［9］ 徐愛功，李娜，張濤.時間序列分析在地鐵沉降觀測中的應用［J］.測繪科學，2013，38（3）：57－60.

［10］ 王建波，欒元重，許君一，等.小波分析橋梁變形監測數據處理［J］.測繪科學，2012，38（3）：79－81.

［11］ 郭濱濱.快速傅里葉變換信號識別能力分析［J］.現代電子技術，2007，30（7）：58－60.

［12］ 陳振恒，高蕾，張穎燕.用快速傅里葉變換實現譜分析的相位噪聲測量［J］.現代電信科技，1994，7：011.

［13］ 占臘生，何娟美，葉藝林，等.太陽活動周期的小波分析［J］.天文學報，2006，47（2）：166－174.