一种电离层f_0F_2和M_3000_F_2长期预测新方法

  • 格式:pdf
  • 大小:326.23 KB
  • 文档页数:4

第23卷 第5期2008年10月 电 波 科 学 学 报CHINESE JOURNAL OF RADIO SCIENCE Vol.23,No.5 October,2008 文章编号 100520388(2008)0520946204一种电离层f0F2和M(3000)F2长期预测新方法3奚迪龙 李建儒☆ 刘玉梅 赵振维(中国电波传播研究所,山东青岛266107)摘 要 将观测站实测数据的月中值按世界时进行排列,分析得到了电离层参数与太阳活动性指数F12之间的相关系数,然后,根据F12的预测数据,可以预测得到观测站的电离层参数。

之后,通过Kriging地理插值法可以预测得到一定区域内观测站以外任一地点的电离层参数。

上述过程形成了一个完整的区域电离层参数长期预测方法。

关键词 回归关系;地理内插;长期预测中图分类号 TN011;P352.3 文献标志码 AA ne w method for long2term prediction of ionospheref0F2and M(3000)F2XI Di2long L I Jian2ru L IU Yu2mei ZHAO Zhen2w ei(China Research I nstitute of Radio w ave Propagation,Qing dao S handong266107,China)Abstract Arranging mont h median values of observed data at every site by U TC,t he correlation coefficient s between iono sp here parameters and t he sun activity in2dex F12are regressed.Then according to t he p redicted F12,t he ionosp here parame2ters at t he observation stations can be p redicted.Afterwards,using t he Kriging ge2ograp hical interpolation met hod,t he iono sp here parameters at any place can be pre2dicted in t his area near China.Then t he complete iono sp here parameters long2termp rediction met hod is formed by all t he above step s.K ey w ords regression relationship;geograp hical interpolation;long2term p redic2tion1 引 言电离层预测,最重要的是电离层F2层临频f0F2和M(3000)F2因子的预测。

对f0F2和M (3000)F2因子的预测,国际上一般采用国际参考电离层(IRI),该模式现在已发展到2007版本。

在我国,上世纪70年代,文献[1]提出了“亚大地区F2电离层预测方法”,对中国地区f0F2和M(3000)F2的预测有了较大的改进,最后发展成为中国参考电离层(CRI),沿用至今。

文献[2]提出使用改进Krig2 ing技术实时重构区域电离层f0F2的分布。

因工作需要,本文对电离层f0F2和M(3000) F2的长期预测方法进行了研究,预测思路主要来之于文献[1],但在电离层拟合参数和地理插值方法两个方面作了变动,与国际参考电离层IRI比较,本文方法(CIM)的预测精度有了较大改进。

3收稿日期:2007212218 ☆E2mail:ljr22s@.9462 电离层参数随太阳活动变化预测与国际参考电离层IRI不同,我们选用了10.7cm太阳射电流量的12个月流动平均(F12)作为太阳活动性指数,来分析F12与电离层参数f0F2和M(3000)F2因子的相关性。

然后使用亚大地区41个电离层探测站多个太阳活动周期的历史观测数据,国外站点从建站到1990年,国内站点从建站到2006年,统一使用世界时,对相关系数C i重新进行了非线性回归,得到了一套适用于亚大地区的相关系数C i。

然后根据F12的预测值,预测得到41个站的电离层f0F2和M(3000)F2的月中值。

电离层探测站位置分布见表1。

表1 电离层探测站位置序号站名纬度经度序号站名纬度经度1TOMSK56.584.922O KINAWA26.3127.8 2NOVOSIBIRSK54.683.223TA IPEI25.0121.5 3PETROPAVLOVSK53.1158.724GUAN GZHOU23.2113.4 4IR KU TSK52.5104.025A HM EDABAD23.072.6 5CHITA52.0113.526HAIKOU20.0110.3 6KARA GANDA49.873.127MAN IL A14.6121.1 7MANZHOUL I49.6117.428TIRUCHIRA PALL I10.878.3 8KHABARROVSK48.5135.129KODAIKANAL10.277.5 9WA KKANA I45.4141.730TRIRANDRUM8.577.0 10CHAN GCHUN43.9125.331SIN GAPORL 1.35103.8 11WUL UMUQ I43.887.732VAN IMO-2.7141.3 12AL MA2A TA43.276.933DRWIN-12.4130.9 13B EI J IN G39.9116.334TOWNSV ILL E-19.7147.0 14A KITA39.7140.135L A REUN ION-21.155.9 15ASH KHABAD37.958.836BRISBAN E-27.5152.9 16L ANZHOU36.1104.037MUNDARIN G-32.0116.2 17KO KUBUNJ I35.7139.438CANB ERRA-35.3149.1 18YAMA GAWA31.2130.539SAL ISBU R Y-34.7136.7 19L ASA29.691.240HOBAR T-42.8147.0 20CHON GQ IN G29.5106.441KER GU EL EN-49.470.3 21DEL HI28.677.2F12的统计方法如下:F12=1126n+5i=n-5f i+12( f n-6+ f n-6)(1)式中,n为时间序列的第n月; f为每月F10.7的月均值。

在一定地点,一定月份、同一时间(世界时),电离层参量随太阳活动性指数变化的回归方程,仍采用CRI模型中的多项式关系: ^X(t)=6K i=0C i(t)F i12(2)式中,^X(t)代表f0F2和M(3000)F2,根据文献[1]的研究,从f0F2、M(3000)F2与太阳活动性指数呈现的相关特性来看,前者K取2,后者K取1就可以。

使用本文预测方法(CIM)与IRI(2007)对国内多个站点、最近13年f0F2和M(3000)F2进行了预测,并与实测数据进行了比较,按照(3)式得到了各站点预测数据的年均方根误差,如图1和图2所示。

图1 各站f0F2两种模式比较749第5期 奚迪龙等:一种电离层f o F2和M(3000)F2长期预测新方法图2 各站M (3000)F 2两种模式比较均方根误差是用来评估电离层参数预测效果的参数之一,其定义为:RM S =6Ni =0(sc i -yb i )2/n (3)式中,n 表示实测数据个数;sc i 表示第i 月实测数据的统计月中值;yb i 表示第i 月的预测值。

对图中六个站点13年预测数据的误差均方根进行了统计,得到本文预测方法与IRI 的平均预测精度比较如表2所示。

表2 电离层预测方法精度比较-预测误差均方根IRICIM f 0F 2(M Hz )0.5400.445M (3000)F 20.1150.080可以看出,新预测方法的预测精度比IRI 在各站、各年都得到了较好改进。

3 电离层参数在地域上的预测电离层参量存在较为明显的地理位置变化。

我们选用Kriging 地理插值方法,利用41个观测站的预测数据来预测观测站以外一定区域内任一地点某月某世界时的电离层参数f 0F 2和M (3000)F 2。

在Kriging 插值方法中,电离层的地域特性用相关距离d ab 来表征:d ab =(S K (L on (a )-L on (b )))2+(S F (L a (a )-L at (b )))2(4)式中,L on (a )、L on (b )和L at (a )、L at (b )分别表示a 、b 站的地理经度和纬度,S K 、S F 是相关距离比例因子,这说明两站的经纬度差异对d ab 的作用可能是不相同的。

用多个电离层观测点的实测值Z i (i =1,2,…,n )来求某个位置的电离层参量Z 0: Z 0=6ni =1WiZ i(5)其中滤波权函数W 由C ・d 来求得。

C 是一个常数。

6VijW j =V j 0 i ,j =1,2,…n 6ni =1W i -1=0 V ij =C 3d ij V j 0=C 3d j 0(6)V ij 是(n +1)(n +1)个单元的矩阵,W i 和V j 0为(n +1)元素的矢量。

由(6)式可得W i ,从而可由(5)式求得预测量Z 0。

参量C 、S F 和K 见表3。

表3 C 、S F 和相关系数KS F 0.20.60.8 1.0 1.6 2.0C 2.0 1.00.30.70.70.2K0.690.840.960.890.730.80根据文献[3]C 和S F 取相关系数K 最大(0.96)所对应的值:S F =0.8,C =0.3。

在实际应用中我们发现经度的距离因子S K 比纬度的距离因子S F 更重要,C 的变化对电离层参数预测精度的影响不明显,我们取C =1。

通过最优化计算S F 取1.0,S K 取0.3比较合适。

见表4。

表4 两组系数K riging 插值计算结果比较站点S K =0.3,S F =1.0S K =1.2,S F =0.8K RMS K RMS 北京0.980.420.950.95长春0.970.510.970.73重庆0.980.620.980.76广州0.980.700.980.89海口0.970.760.980.75兰州0.990.540.980.91拉萨0.970.820.96 1.08满洲里0.960.630.960.64乌鲁木齐0.980.440.970.60 注:在插值某个站的电离层参数时,该站的数据不参与计算。