一种低频率GNSS电离层闪烁因子的有效性验证方法技术

本发明专利技术公开一种验证低频率GNSS电离层闪烁因子有效性的方法，适用于监测电离层闪烁天文灾害。首先对天文数据和地磁活动数据预处理，然后统计低采样频率GNSS电离层闪烁因子的闪烁探测发生率与ISMR电离层闪烁因子探测的电离层闪烁发生率的统计比较分析，若电离层闪烁发生率相关性高则认为待验证闪烁因子有效性高；再通过太阳辐射和地磁活动指数与待验证闪烁因子和ISMR电离层闪烁因子给出的日、时电离层闪烁发生率的相关性，判断二者相关性，从而再次验证电离层闪烁因子有效性高。其步骤简单，能够有效筛选出较优的闪烁因子，助力其推广应用，进而大幅度降低电离层不规则体监测的成本，有利于对全球电离层闪烁形成完备监测，具有广泛的实用性。具有广泛的实用性。具有广泛的实用性。

【技术实现步骤摘要】
一种低频率GNSS电离层闪烁因子的有效性验证方法


[0001]本专利技术涉及一种验证低频率GNSS电离层闪烁因子有效性的方法，尤其适用于适用于监测电离层闪烁天文灾害的一种低频率GNSS电离层闪烁因子的有效性验证方法。
技术背景
[0002]电离层闪烁是赤道和极区频发的一种天文灾害，是近年来空间科学面临的一大难题，对人类社会可持续发展构成了重大威胁。准确监测电离层闪烁，对于掌握其变化规律，预防电离层闪烁天文灾害的不利影响具有重大现实意义。当前，电离层闪烁监测研究的主流设备通常采用电离层闪烁监测接收机，简称ISMR。ISMR可给出高精度的电离层闪烁因子，它是对电离层闪烁强弱的定量表征。但ISMR价格较贵，且采样频率多设在50Hz，需要较大的存储空间，进而限制了布站数量，尚无法利用ISMR进行全球电离层闪烁研究。
[0003]相比ISMR，测地型接收机价格更为低廉，采样频率低，分布更为广泛，全球布站多大数万个。提取基于测地型接收机的电离层闪烁因子，更有助于对全球电离层闪烁开展完备监测。近来年，研究人员提出了众多的基于低采样频率的GNSS电离层闪烁因子。Pi等利用历元间差分载波相位无几何关系组合估算斜路径总电子含量的方法，消除了低频率误差对提取电离层闪烁信号的影响，构建了ROTI因子。Juan等对无几何关系组合做历元内全部可视卫星加权平均，并引入周跳探测与修复方法，提取了可反映电离层活动强度时空变化信息的AATR因子。Nguyen等采用0.1Hz截止频率的巴特沃斯滤波技术，分离了多路径效应对电离层闪烁因子的不利影响。Ahmed等研究了从高纬度区域垂直总电子含量数据中提取电离层闪烁信号的小波变换方法。赵东升等也基于小波变换方法，构建了适用于北极区域的电离层相位闪烁因子，该因子可以给出每个频率上的相位闪烁幅值。
[0004]以上基于低采样频率电离层闪烁因子的提出丰富了电离层闪烁监测研究的手段，但尚未得到较大规模的应用，其主要原因是这些闪烁因子有效性缺乏长时间的数据验证。尽管研究人员在提出闪烁因子时也进行了一定的有效性验证，但多关注于该方法在几天甚至几个卫星观测弧段对电离层闪烁监测的准确性。这种准确性多将ISMR输出的闪烁因子作为参考值，按弧段将所估计的电离层闪烁因子与参考值进行比较。基于观测时间较短的数据进行的验证结果可能存在特例偏差，难以说明闪烁因子的有效性。目前，尚未见适用于长时间观测数据验证低频率GNSS电离层闪烁因子有效性的统计分析方法。

技术实现思路

[0005]针对现有技术的不足之处，提供一种基于统计规律验证低采样频率GNSS电离层闪烁因子有效性的方法，利用低采样频率GNSS电离层闪烁因子在多种条件下与参考值的规律统计分析，解决了难以对电离层相位闪烁因子进行长时间验证的技术难题。
[0006]为实现上述技术目的，本专利技术的一种低采样频率GNSS电离层闪烁因子的有效性验证方法：
[0007]S1、获取电离层闪烁监测接收机ISMR的电离层闪烁因子数据、日出日落时间数据
以及表征太阳和地磁活动物理参数，然后对以上数据及参数进行预处理，剔除多路径效应，以及f10.7指数和Ap指数数据中存在断点的不利影响；
[0008]S2、选择具有电离层闪烁监测接收机的观测站获取的1Hz采样频率数据，通过文献获取待验证的低采样频率电离层闪烁因子，利用相关性分析方法，统计待验证的低采样频率电离层闪烁因子的闪烁探测发生率，与通过ISMR电离层闪烁因子探测的电离层闪烁发生率，之后将闪烁探测发生率与电离层闪烁发生率进行比较分析，若待验证闪烁因子探测出的电离层闪烁发生率与ISMR给出的电离层闪烁发生率相关性高且差值不大于0.5％，则认为待验证的低采样频率电离层闪烁因子有效性高；
[0009]S3、最后对太阳辐射和地磁活动相关性统计分析，通过分别研究太阳辐射和地磁活动指数与日、时电离层闪烁发生率的相关性，若活动指数与日、时电离层闪烁发生率的相关性较高，太阳辐射与日、时电离层闪烁发生率的相关性较高，说明由待验证闪烁因子探测出的电离层闪烁因子的发生的驱动因素与ISMR电离层闪烁因子给出的驱动因素相同，由此可以判断待验证的低采样频率电离层闪烁因子有效性高，反之则判断待验证的低采样频率电离层闪烁因子有效性低。
[0010]进一步，ISMR电离层闪烁因子数据、日出日落时间数据以及表征太阳和地磁活动物理参数的预处理，包括如下三步：
[0011]首先剔除多路径效应对ISMR电离层闪烁因子的影响：即利用IGS提供的精密星历文件和测站的坐标计算选择出的具有电离层闪烁监测接收机的观测站上空每个卫星的高度角，删除高度角小于30
°
的ISMR电离层闪烁因子，从而提取出受多路径影响较小的ISMR电离层闪烁因子，并将ISMR电离层闪烁因子作为评估低采样频率GNSS电离层闪烁因子有效性的参考值；
[0012]然后利用Jean Meeus所提出的天文算法计算用来验证低采样频率GNSS电离层闪烁因子的测站太阳时角，并根据选择的一处配备有ISMR接收机的GNSS测站位置所处时区，计算太阳日出和日落的地方时；
[0013]最后用加拿大自然资源局提供的f10.7指数表征太阳活动强弱，用德国地学中心提供的Ap指数和并址磁测站提供的当地地磁数据表征全球和当地的地磁活动强弱，由于f10.7和Ap指数数据中存在断点问题，因此找出断点所在，并将缺失的数据赋值为空，从而避免其不良影响。
[0014]进一步，待验证的低采样频率GNSS电离层闪烁因子和ISMR的闪烁发生率为：
[0015]以天为单位获取待验证闪烁因子的日电离层闪烁发生率P，发生率的计算方法如下，
[0016][0017]其中，P表示闪烁发生率，发生闪烁的历元即该历元的待验证闪烁因子幅值大于ISMR电离层闪烁因子的阈值0.2rad，待验证的低采样频率GNSS电离层闪烁因子的阈值通过现有文献获取，该阈值用以判定某一观测历元是否有闪烁发生。
[0018]进一步，通过计算待验证低采样频率GNSS电离层闪烁因子给出的日电离层闪烁发生率P与ISMR给出的日电离层闪烁发生率的差值及相关性：认为90％差值的分布在
±
0.5％以内，且相关性在0.9以上则为该待验证低采样频率GNSS电离层闪烁因子的探测每日的电
离层闪烁；
[0019]进一步，利用选择出具有电离层闪烁监测接收机的观测站的日出和日落时间为基准，统计待验证闪烁因子和ISMR电离层闪烁因子在1年时间内在日落后和日出后的每小时电离层闪烁发生率，由于电离层闪烁在日出和日落后的发生具有一定的规律性，即峰值发生的时间固定；若测站位于北纬66
°
34
’
以北和南纬66
°
34
’
以南的测站，在进行该项验证时需要考虑到极昼和极夜现象，此时可以仅以当地时为基准，比较峰值的发生时刻；当待验证闪烁因子探测出的发生峰值时刻与ISMR电离层闪烁因子给出的峰值时刻差值在半小时内时，可认为待验证闪烁因子能探本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种低采样频率GNSS电离层闪烁因子的有效性验证方法，其特征在于：S1、获取电离层闪烁监测接收机ISMR的电离层闪烁因子数据、日出日落时间数据以及表征太阳和地磁活动物理参数，然后对以上数据及参数进行预处理，剔除多路径效应，以及f10.7指数和Ap指数数据中存在断点的不利影响；S2、选择具有电离层闪烁监测接收机的观测站获取的1Hz采样频率数据，通过文献获取待验证的低采样频率电离层闪烁因子，利用相关性分析方法，统计待验证的低采样频率电离层闪烁因子的闪烁探测发生率，与通过ISMR电离层闪烁因子探测的电离层闪烁发生率，之后将闪烁探测发生率与电离层闪烁发生率进行比较分析，若待验证闪烁因子探测出的电离层闪烁发生率与ISMR给出的电离层闪烁发生率相关性高且差值不大于0.5％，则认为待验证的低采样频率电离层闪烁因子有效性高；S3、最后对太阳辐射和地磁活动相关性统计分析，通过分别研究太阳辐射和地磁活动指数与日、时电离层闪烁发生率的相关性，若活动指数与日、时电离层闪烁发生率的相关性较高，太阳辐射与日、时电离层闪烁发生率的相关性较高，说明由待验证闪烁因子探测出的电离层闪烁因子的发生的驱动因素与ISMR电离层闪烁因子给出的驱动因素相同，由此可以判断待验证的低采样频率电离层闪烁因子有效性高，反之则判断待验证的低采样频率电离层闪烁因子有效性低。2.根据权利要求1所述的一种低采样频率GNSS电离层闪烁因子的有效性验证方法，其特征在于ISMR电离层闪烁因子数据、日出日落时间数据以及表征太阳和地磁活动物理参数的预处理，包括如下三步：首先剔除多路径效应对ISMR电离层闪烁因子的影响：即利用IGS提供的精密星历文件和测站的坐标计算选择出的具有电离层闪烁监测接收机的观测站上空每个卫星的高度角，删除高度角小于30
°
的ISMR电离层闪烁因子，从而提取出受多路径影响较小的ISMR电离层闪烁因子，并将ISMR电离层闪烁因子作为评估低采样频率GNSS电离层闪烁因子有效性的参考值；然后利用JeanMeeus所提出的天文算法计算用来验证低采样频率GNSS电离层闪烁因子的测站太阳时角，并根据选择的一处配备有ISMR接收机的GNSS测站位置所处时区，计算太阳日出和日落的地方时；最后用加拿大自然资源局提供的f10.7指数表征太阳活动强弱，用德国地学中心提供的Ap指数和并址磁测站提供的当地地磁数据表征全球和当地的地磁活动强弱，由于f10.7和Ap指数数据中存在断点问题，因此找出断点所在，并将缺失的数据赋值为空，从而避免其不良影响。3.根据权利要求1所述的一种低采样频率GNSS电离层闪烁因子的有效性验证方法，其特征在于待验证的低采样频率GNSS电离层闪烁因子和ISMR的闪烁发生率为：以天为单位获取待验证闪烁因子的日电离层闪烁发生率P，发生率的计算方法如下，其中，P表示闪烁发生率，发生闪烁的历元即该历元的待验证闪烁因子幅值大于ISMR电离层闪烁因子的阈值0.2rad，待验证的低采样频率GNSS电离层闪烁因子的阈值通过现有文
献获取，该阈值用以判定某...

【专利技术属性】
技术研发人员：赵东升，李旺，张秋昭，唐旭，王潜心，张克非，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：