用时间序列分析理论预报电离层.pdfVIP

用时间序列分析理论预报电离层 李伟超1,3程宗颐2冯初刚2李志刚 (1．中科院国家授时中心·陕硬西安·710600；2．上海天文台·上海·200030) (3．中国科学院研究生院·北京·100039) 摘要提出用时间序列分析理论中的ARMA(p，I，q)模型对电离层进fS-预报，即本文所称的 直接法，或者用球谐函数对电离层电子含量进行拟舍所得的拟合系数所构成的时问序列，再 用时间序列分析理论中的ARMA(p，I，q)模型对拟合系数进行预报，然后用所预报的球谐函数 系数算出电子含量，这称为间接法．利用2004年1月1日至2005年1月31IGS所给出的泌 电离层资料对本文提出的两种方法进行检验和比较．结果表明，在大约12天以内，上述两 种方法的预报结果基本一致，长于12天的预报时，间接法的精度高于直接法．在20天的时 间内，预报值与已知值之差小于3个TECU的电离层点数占总点数的百分比基本上都在80％ 左右；随着时间的推移，这个比值的下降对直接法并不明显，而对间接法则有所下降．直接 法适用于区域洼的10天左右的电子含量的预报；而间接法适用于全球性的短，中期(20天 以内)的全球电子含量的预报． 关键词 电离层；电子含量；时间序列；电离层预报 1 引言 地球上空的电离层对卫星通讯，航空，航天，天气预报等方面都产生重要的影响。因此测量电离 层中的电子含量，研究其变化规律并对它加以预报一直是研究的热门课题。目前世界上有几个电离层 电子含量测量中心，他们几乎每天公布其所得地球表面一定高度上空的电子含量。电离层的预报模型 75％，且其在夜间的预报具有明显的不合理性。因此有必要进行电离层屯子含量的精确预报的研究。 2用ARMA(p，q)模型预报电离层的电子浓度 从我们所用时间序列的物理背景看，我们选用更为普遍适用的ARMA(p，q)模型[5]。对模型参数的 确定我们基本采用模型参数的逆函数估计法。对于预报我们采用建立在平稳线性最小方差预报理论基 础上的对ARMA(p。CO模型进行预报。在预报之前，我们对每一组时间序列进行预处理。首先娃零均值 化，然后是用参数检验法进行稳定性检验。对非平稳序列做差分处理，使之平稳化。做差分处理时． 我们还考虑季节性变化囡素的作用。因此在实际计算时，有时是用Ag．MA(p，q模型预报，而更多的采 用的是ARIMA(p，d，q)模型预报。 3实际的预报 资料是每两小时、在经度方向每5度、纬度方向每2。5度给出一个电子浓度以及它的精度(RMS)。 法来预报。一是直接法；二是间接法。 3．1直接法预报电子含量 直接法就是对每个网格点上的电子含量所形成的时间序列进行直接预报，也就是要预报5184个 预报值，绿色为IGS所给已知值。横坐标表示时间，以2小时为单位。纵坐标表示电子含量，以TECU 为单位。我们对5184个时间序列都作了类似的预报。从而得到以2小时为间隔的20天的全球电子含 量的预报。 三广■『_■■■。”『] ：挝∑丛蚣盟瞧丛：△丛△△．丛堂-_』。， 图1 电离层球壳上纬度15。，经度2500的空间点处的电子含量的预报20天的示意图 3．2间接法预报电离层电子含量 间接法就是对某时刻的距地面某一高度的球面上的电子含量用如下球面函数进行拟合： max^H (1) n 得到了a。，及k。的时间序列。反过来，如果己知某时刻的n。。及“。，根据(1)式可计算该时刻同 样高度上纬度为a，经度为卢处的电子浓度。因此如果我们能根据已有的a。。及k。的时间序列预报 未来时刻的口。，及k。，就等于预报了未来时刻的全球同样高度上的电子浓度。(1)式中的a，p可以 在地理框架，日固框架，太阳一地磁框架等中量度，都可利用(1)式进行进行电离层的间接预报。 3．3在地理框架，日固框架，太阳．地磁框架中的预报结果及与直接预报结果的比较 取前面直接法预报所用的相同资料利用(1)式分别在地理框架，日l司框架，太阳一地磁框