一种结合NeQuickG模型的电离层层析方法

本发明属于电离层探测，特别涉及一种结合nequick g模型的电离层层析方法。

背景技术：

1、随着世界各国对于地外空间探测活动的日益频繁，以及卫星通信导航遥感的应用普及，人类对电离层结构特征及其时空变化规律的监测反演需求日益增强，同时电离层监测反演对于理解其各层结构的形成机制也具有重要的科学意义。电离层探测是获取电离层基本信息和结构特征的重要手段，主要包含地面测高仪、雷达、电离层模型、全球导航卫星系统及电离层层析成像(cit)等技术。

2、电离层层析方法展现出了探测电离层电子密度三维分布和变化方面的优势，可以重构电离层空间结构，适合监测电离层电子密度的大尺度空间分布及其变化，具有费用低、操作简单、探测范围广等诸多优势，对于电离层不同尺度结构变化及全球电离层环境监测具有重要意义，但是电离层层析方法的精度方面还有较大的提升空间。

3、电离层层析方法中需要构建映射矩阵，以解决矩阵秩亏的问题。构建映射矩阵的方法有很多，其中一种就是使用经验正交函数(eof)来构建垂直轮廓，并使用基于体素的网格来绑定纬度和经度，利用经验正交函数对电离层垂直剖面形成约束。传统的电离层层析方法使用国际参考电离层(iri)模型生成经验正交函数，所得反演结果的精度较低，所以可使用精度较高的nequick g模型构建经验正交函数并应用到电离层层析当中。与国际参考电离层模型相比，nequick g模型具有以下特点：nequick g模型是三维时变的电离层经验模型，该模型是以epstein公式为基础，以e层、f1层和f2层峰值数据为锚点建立解析函数公式，能够对参数进行近似实时的改正。nequick g模型不仅能计算出电离层中f2层的峰值电子密度(nmf2)和峰值高度(hmf2)，也可以计算出任何位置的电子密度、任意传播信号的倾斜总电子含量(stec)，并且在全球范围内表现出了较高精度，可消除50％以上的电离层误差，所以将nequick g模型的结果生成一个经验正交函数，作为层析技术的映射矩阵，应用到电离层层析反演中，计算得到单点位置处的电子含量ne、倾斜总电子含量是极具意义的。但目前尚未发现相应的方法。

技术实现思路

1、为了解决上述问题，本发明的目的在于提供一种结合nequick g模型的电离层层析方法。

2、为了达到上述目的，本发明提供的结合nequick g模型的电离层层析方法包括按顺序进行的下列步骤：

3、1)确定电离层层析反演区域，并根据gnss观测文件获得电离层层析反演区域内接收站的具体地理位置；

4、2)根据电离层层析反演区域的精密星历文件以及步骤1)获得的接收站的具体地理位置，获取卫星的地理位置，由此确定出卫星的轨迹，并设定电离层层析反演时间；

5、3)根据步骤1)获得的电离层层析反演区域、步骤2)获得的电离层层析反演时间获得月份、世界时，并根据galileo导航电文计算出nequick g模型所需的电离层水平有效参数az；

6、4)根据步骤1)获得的接收站的具体地理位置、步骤2)获得的卫星的地理位置、步骤3)获得的电离层水平有效参数az，利用nequick g模型计算出三维空间区域中单点位置处的电子密度ne；

7、5)利用步骤4)获得的三维空间区域中单点位置处的电子密度ne生成经验正交函数eof作为电离层层析反演的映射矩阵m；

8、6)利用电离层层析算法，基于步骤2)获得的卫星轨迹、步骤5)获得的映射矩阵m，反演三维空间区域中单点位置处的电子密度ne与信号传输路径的倾斜总电子含量stec，之后根据卫星的仰角计算得到垂直总电子含量vtec。

9、在步骤1)中，所述确定电离层层析反演区域，并根据gnss观测文件获得电离层层析反演区域内接收站的具体地理位置的方法是：

10、首先确定电离层层析反演区域，然后找到该区域内接收站的观测文件，将观测文件中多个接收站的经度、纬度、高度数据提取出来并合并到一个文件中，由此获得接收站的具体地理位置。

11、在步骤2)中，所述根据电离层层析反演区域的精密星历文件以及步骤1)获得的接收站的具体地理位置，获取卫星的地理位置，由此确定出卫星的轨迹，并设定电离层层析反演时间的方法是：

12、获取电离层层析反演区域的精密星历文件，然后根据步骤1)中获得的接收站的具体地理位置和精密星历文件，计算出不同卫星在不同时刻所处的包括经度、纬度和高度在内的地理位置，由此确定出卫星的轨迹，同时将电离层层析反演时间间隔设置为30min。

13、在步骤3)中，所述根据步骤1)获得的电离层层析反演区域、步骤2)获得的电离层层析反演时间获得月份、世界时，并根据galileo导航电文计算出nequick g模型所需的电离层水平有效参数az的方法是：

14、根据步骤1)获得的电离层层析反演区域、步骤2)获得的电离层层析反演时间获取所求时刻的月份、世界时，然后在galileo导航电文中找到当日播发的一组电离层水平有效参数az。

15、在步骤4)中，所述根据步骤1)获得的接收站的具体地理位置、步骤2)获得的卫星的地理位置、步骤3)获得的电离层水平有效参数az，利用nequick g模型计算出三维空间区域中单点位置处的电子密度ne的方法是：

16、在利用nequick g模型计算三维空间区域中单点位置处的电子密度ne过程中，epstein函数被用作构建电离层的基本分析函数，其中包含了电离层剖面中的e层峰值高度hme、f1层峰值高度hmf1、f2层峰值高度hmf2这三个定位点以及相应的临界频率foe、fof1和fof2和传播因子m(3000)f2在内的参数；

17、epstein函数的解析式表示为：

18、

19、式中，就电离层剖面而言，x1表示某层的峰值振幅，y1表示该层的峰值高度，z1表示该层的厚度，w1表示该层的高度；

20、nequick g模型的输入参数为接收站的包括经度、纬度、高度在内的具体地理位置、卫星的包括经度、纬度、高度在内的地理位置、世界时和电离层有效水平参数az，以及ccir－xx.asc文件提供的f2层临界频率fof2和传播因子m(3000)f2，modip.asc文件提供的地磁纬度，输出参数为空间区域中单点位置处的电子密度ne。

21、在步骤5)中，所述利用步骤4)获得的三维空间区域中单点位置处的电子密度ne生成经验正交函数eof作为电离层层析反演的映射矩阵m的方法是：

22、将步骤4)中获得的一系列三维空间区域中单点位置处的电子密度ne按经度、纬度、高度整合为一个矩阵，然后对该矩阵使用svds函数进行奇异值分解，获得奇异值和奇异向量，再对奇异向量进行矩阵标准化，获得经验正交函数eof作为电离层层析反演的映射矩阵m。

23、在步骤6)中，所述利用电离层层析算法，基于步骤2)获得的卫星轨迹、步骤5)获得的映射矩阵m，反演三维空间区域中单点位置处的电子密度ne与信号传输路径的倾斜总电子含量stec，之后根据卫星的仰角计算得到垂直总电子含量vtec的方法是：

24、gnss接收机的rinex文件包含的gps可观测数据有:

25、

26、

27、

28、

29、其中，p1和p2为精码中的码伪距，p0为电离层自由伪距，l1和l2为载波相位值，f1与f2分别为gps的第一载波信号频率和第二载波信号频率，λ1与λ2分别为gps的第一载波信号频率f1与第二载波信号频率f2对应的载波波长，n1与n2分别为第一载波信号频率f1与第二载波信号频率f2观测的整周模糊度，ε1与ε2分别为卫星和接收机的硬件误差分量；tec为电子总含量；

30、由于相位和码距的观测数据不同，因此由式(2)至式(5)可推导出如下两式：

31、

32、

33、可通过式(6)中的噪声项计算得到倾斜总电子含量stec，通过式(7)中与其相关联的整周模糊度的偏置项也可得到倾斜总电子含量stec；

34、基于上述倾斜总电子含量stec，利用多分辨率电离层层析算法反演获得电离层的电子密度x；方法是首先设gnss接收机与卫星之间的传输路径为h，再将电离层格网化，将无线电波穿过每个格网时在垂直方向上的截距记为h1，h2，…，hi，……，假设每个格网中的电子密度为确定的值x1，x2，…，xi，…，那么单位格网内的总电子含量为hixi，无线电波传输路径上的倾斜总电子含量stec为射线穿过格网的截距与该网格内的电子密度的乘积之和，即stec＝h1x1+h2x2+…+hixi+…；则所有无线电波射线上的倾斜总电子含量stec的方程为hx＝z，其中x即为所求的电离层的电子密度，h为所有射线穿过格网在垂直方向上的截距矩阵；由于射线穿过格网的路径并不完备，出现了矩阵秩亏的问题，为了解决该问题，利用步骤5)中由nequick g模型得到的映射矩阵m来构建垂直轮廓，这样经验正交函数eof就对垂直剖面形成了约束；然后对矩阵映射m进行正则化，这时方程就由hx＝z变为(hm)x＝z，然后引入正则化矩阵r，构建正规方程hth+krtr＝htz，其中k为用户定义的正则化常数，最后利用最小残差法计算上述正规方程的解，即可得到作为电离层层析反演结果的电离层的电子密度x；

35、将电离层层析反演得到的三维空间区域中单点位置处的电子密度ne沿信号传播路径进行积分，即可得到电离层层析算法的倾斜总电子含量stec，之后根据已知的卫星仰角即可得到垂直总电子含量vtec：

36、

37、

38、式中，re＝6371.2km是地球的平均半径，h1＝450km是电离层穿透点(ipp)的高度，e是卫星仰角，单位为弧度。

39、本发明提供的结合nequick g模型的电离层层析方法具有如下有益效果：由于nequick g模型精度较高，本发明先是利用nequick g模型构建经验正交函数来构建垂直轮廓，eof对垂直剖面形成了约束，将经验正交函数作为映射矩阵解决电离层层析中的矩阵秩亏问题，使层析反演的精度达到一定的水平。