-
摘要:
针对中国区域连续运行参考站接收机野值会干扰星基增强系统(SBAS)电离层异常事件提取的问题,提出了一种基于电离层垂直延迟时间梯度的野值检测方法。首先,介绍了电离层延迟数据的提取方法;然后,论述了依据野值和电离层异常的不同时空相关特性进行野值检测的方法,并用单双频定位误差结果验证了野值检测的正确性;最后,对检测结果进行了分析。结果表明:该方法可有效区分野值和电离层异常;在中国大陆构造环境监测网络200多个参考站中,野值检测所剔除的参考站数目在3~10个,对电离层穿透点空间分布的影响在可接受范围内。
-
关键词:
- 中国区域 /
- 星基增强系统(SBAS) /
- 电离层 /
- 延迟值异常 /
- 野值检测
Abstract:To solve the problem that continuously operating reference stations' receiver's outliers in China will interfere with ionospheric anomaly event extraction of satellite-based augmentation system (SBAS), an outlier detection method based on the ionospheric vertical delay time gradient is proposed. First, ionospheric delay data extraction method was introduced. Then, the ionospheric vertical delay time gradient outlier was detected based on the different spatio-temporal correlation characteristics of outlier and ionospheric anomaly. The correctness of outlier detection was verified by single-frequency and dual-frequency positioning error results. Finally, the detection results were analyzed. The results show that the method can effectively distinguish between outliers and ionospheric anomalies; among more than 200 reference stations of crustal movement observation network of China, the number of reference stations rejected by the outlier detection is about 3 to 10, and the effect on the spatial distribution of the ionospheric penetration points is within an acceptable range.
-
对于全球导航卫星系统(Global Navigation Satellite System,GNSS)的单频用户,电离层延迟是最大的定位误差源。为保证航空用户的安全,星基增强系统(Satellite Based Augmentation System,SBAS)为单频用户提供网格点电离层改正数和改正误差[1]。由于观测值不充分、电离层结构复杂、异常事件难以预测等原因,电离层异常是影响SBAS完好性的最大风险项[2]。电离层在不同区域具有不同的结构特点。其中,中国国土范围跨越中低纬地区,与全球平均电离层水平相比,中国电离层变化更为复杂。文献[3]研究了中国区域SBAS电离层网格校正技术,针对中国南部区域的异常特性,提出了南北分区的研究思路。文献[4-7]分析了中国电离层异常现象的产生原因与影响,其中包括赤道异常、电离层暴、电离层闪烁等典型电离层异常现象。为保证SBAS完好性,需要针对中国区域复杂的电离层分布特性,开展中国区域的电离层异常研究工作。
目前,美国利用广域增强系统(Wide Area Augmentation System,WAAS)分布在美国本土、阿拉斯加、夏威夷岛和墨西哥的多个参考站收集的高精度观测数据进行电离层异常事件的提取和分析,以建立适用于SBAS的电离层异常模型和处理方法[8-11]。由于WAAS同一站配备3套接收机,数据质量得到了保障,电离层延迟值提取较为准确。研究国内电离层特性,可使用可靠性较高的连续运行参考站(Continuously Operating Reference Stations,CORS)的实测数据,但其数据质量无法达到WAAS站的水平。袁运斌等[12-13]研究了中国区域电离层延迟精确建模方法,使用CORS数据拟合电离层延迟值时,由于模型的平滑作用,拟合延迟值对接收机野值并不敏感;但进行完好性研究时,由于接收机野值具备和电离层异常相似的“异常”特性,会严重干扰电离层异常事件的提取。因此,在电离层异常研究中,需要对各个参考站的数据进行筛查以排除非电离层因素(例如卫星故障或接收机故障等)导致的数据异常,以保证电离层异常建模准确有效。
针对上述问题,本文定义由于非电离层因素导致的数据异常为“野值”,并提出一种基于电离层垂直延迟时间梯度的野值检测方法。
1. 电离层延迟数据提取
电离层垂直延迟数据是电离层异常研究的基础。本文采用数据质量较为可靠的GPS卫星L1和L2频点的双频码和载波相位观测值进行电离层延迟值的解算。L1和L2双频码伪距和载波相位伪距观测方程可表示为[14]
(1) 式中:r为接收机到卫星的真实距离;c为光速;δtu为接收机钟差;δts为卫星钟差;T为对流层延迟;I为L1频率信号上的电离层延迟。由于电离层延迟与信号频率的平方成反比,因此可通过给I乘上一个比例因子得到L2频率信号上的电离层延迟γI。下标1表示L1频率信号,下标2表示L2频率信号。f1和f2为信号频率;λ1和λ2为波长;d1r和d2r为接收机硬件延迟;d1s和d2s为卫星硬件延迟;ερ1和ερ2为码噪声;N1和N2为载波测量值的整周模糊度;εϕ1和εϕ2为载波相位噪声; ρ1和ρ2为码伪距观测值;ϕ1和ϕ2为载波相位观测值。式(1)中双频码伪距值相减再除以常数项可得到L1的电离层延迟码观测值Iρ1;双频载波伪距值相减再除以常数项可得到L1的电离层延迟载波观测值Iϕ1:
(2) 其中:延迟观测值Iρ1和Iϕ1可作为延迟真实值I的估计,但其中还包含噪声以及码间偏差(Differential Code Biases,DCB)项的干扰,DCB也称为内频偏差,是由于卫星和接收机的硬件电路对不同频率的信号产生的硬件延迟不同导致的[15]。双频码伪距计算得到的电离层延迟值噪声大;载波相位伪距噪声小,却有整周模糊度的问题。因此可利用载波相位平滑码伪距,得到低噪声的电离层延迟观测值,而平滑之前还需要进行周跳探测以确保载波相位观测值的连续性。参考Bernese软件的后处理方法进行周跳探测以及平滑[16],并建立区域垂直总电子含量(Vertical Total Electron Content,VTEC)球谐模型以估计卫星和接收机的DCB值,最后得到准确的电离层垂直延迟数据[17]。
欧洲定轨中心(Center for Orbit Determination in Europe,CODE)是国际上权威的GNSS分析机构,以其发布的全球电离层图(Global Ionosphere Maps,GIM)产品内插得到的延迟值作为参考,与处理得到的电离层垂直延迟值结果进行一致性比对,结果如图 1所示。结果表明,相比于单频定位普遍采用的Klobuchar模型值,本文处理得到的双频观测值更接近CODE电离层产品,证明了该方法可用于中国区域电离层延迟数据的提取。
通过分析结果发现,电离层垂直延迟双频观测值有时会出现跳变、梯度增大、剧烈抖动等异常现象。由式(2)可看出,除电离层因素引起真实延迟值I异常会导致电离层延迟观测值异常外,卫星故障、接收机故障、接收机受到干扰引起噪声增大等非电离层因素也会导致延迟观测值出现异常,这种异常会干扰对电离层真实特性的研究和分析,因此视其为“野值”,需要剔除。
卫星端异常(例如星历故障、星钟故障等)会导致不同参考站的同一颗卫星的观测值发生异常变化;接收机端异常(例如接收机故障、受到电磁干扰等)会导致单个参考站的观测值发生异常变化。由于卫星发生故障的概率非常小,且卫星端具备自主完好性监测能力,目前处理的数据还未出现卫星端异常的事件。因此实际数据中“野值”产生的原因主要是接收机异常,可通过剔除参考站数据的方法进行处理。
2. 野值检测方法
以利用上述方法处理得到的电离层垂直延迟数据作为基础,进行野值检测,剔除由非电离层因素导致的异常数据,并筛选出由电离层因素导致的异常数据。根据电离层活动的物理特性,可对电离层延迟异常值进行初步的识别。电离层异常导致的数据异常一般具有空间相关性和时间相关性,即在一定范围内,多于一个站或一颗星的电离层垂直延迟值在相同时间段内发生趋势相近的异常变化。而非电离层因素引起的数据野值,往往表现为不同卫星在不同时间段的电离层延迟值的跳变。
本文选取每个站每颗卫星的电离层垂直延迟值序列进行差分处理,得到电离层垂直延迟时间梯度(单位:mm/s)。电离层异常或者接收机异常都会导致时间梯度增大。通过对中国大陆构造环境监测网络(简称中国陆态网)200多个站的数据进行处理分析,计算得到每个站所有卫星的电离层垂直延迟时间梯度。以2016年3月26日(磁静期)为例,所有站的所有历元的时间梯度概率分布如图 2所示(样本数为4 896 618)。结果表明90%以上的时间梯度值都集中在(-0.5, 0.5) mm/s的区间内,仅存在少部分由于电离层异常或非电离层异常因素而增大的时间梯度值,最大异常梯度甚至超过20 mm/s。
利用CODE电离层产品计算出的电离层垂直延迟时间梯度最大值为1~2 mm/s,可认为实际电离层延迟正常梯度在其范围之内。由于接收机热噪声以及预处理时弧段划分等原因,CORS数据解算出的双频延迟值时间梯度会产生2~3 mm/s的干扰。因此可认为约5 mm/s以内的时间梯度都为正常值。根据对数据量的不同要求可调整检测门限以满足不同的场景。电离层异常研究中需采用尽可能密集的穿透点延迟数据以捕捉中国区域小时空尺度下的异常事件,因此不能剔除过多的数据。考虑到研究中国区域电离层异常特性的数据需求,以及野值对电离层异常检测的影响,本文选择了10 mm/s作为门限对数据进行筛选。根据电离层异常和非电离层异常的不同特性,以实际经验数据为基础,提出一种野值检测和剔除的方法,其流程如图 3所示。首先,计算单站单星的电离层垂直延迟时间梯度序列并进行遍历,若时间梯度绝对值超过10 mm/s,则判定其为异常值(野值或是电离层异常)。然后,根据该点附近2 h间隔内是否有其他历元的时间梯度绝对值超过5 mm/s来区分电离层异常或是野值。若单个站的野值个数小于2,则剔除该站该星的异常数据,否则剔除该站所有数据。最后,利用剔除站的单双频定位结果进行复核,验证其为野值。
对不同检测结果的参考站数据进行具体分析,可发现由于不同异常因素而导致的不同时间梯度特性。图 4~图 7分别是4个站所有卫星的电离层垂直延迟时间梯度值随时间变化。
图 4是2016年3月26日(磁静期)位于中国北部的HLAR站的电离层垂直延迟时间梯度。该站所有数据在野值检测中都被判定为正常值。由于不受到电离层暴和赤道异常等因素的影响,HLAR站所在区域电离层活动很平静,延迟值变化均匀且缓慢,因此时间梯度大部分都很小,个别不连续的点是由于观测值平滑时弧段划分引起的。
图 5是2016年3月26日QHTR站的电离层垂直延迟时间梯度。该站在野值检测中数据出现多个野值,被判定为剔除站。由图 5可见,该站时间梯度的极端异常值可超过±20 mm/s,而且异常值不连续。通过对该站周边站的电离层延迟值进行一致性比对,证明该异常不具有空间相关性,因此判定其产生异常的原因是该站接收机故障或受到干扰。
图 6是2016年3月26日位于海南的QION站的电离层垂直延迟时间梯度。野值检测中,该站出现了电离层异常。由图 6可以看出,13~16 h有几颗卫星的延迟值发生了异常变化,时间梯度出现了连续异常值,且波动较大。对北纬20°附近的多个站进行分析,都发现了相似的异常现象,其发生时间与区域都与电离层等离子体泡现象相吻合[18]。因此判断该站电离层异常的原因是电离层等离子体泡。
图 7是2015年3月17日(磁暴期)HAQS站的电离层垂直延迟时间梯度。野值检测中,该站也出现了电离层异常。该天的地磁指数Kp最大值为7,全球区域发生了强磁暴,因此中国区域有很大可能性发生了电离层暴。由图 7可看出,在GPS时5~10 h(当地时间13~18 h,日间),电离层最活跃也是延迟值变化最快的时段,多颗卫星的电离层垂直延迟值时间梯度连续增大。对比附近站的数据,发现该异常具有空间相关性。与图 6的电离层异常相比,图 7所示时间梯度变化更为缓慢,延迟值变化具有单调性,判定其电离层异常原因为赤道异常引起的电子密度增强或电离层暴引起的风暴增强密度(Storm Enhanced Density,SED)。
无论是哪种原因引起的电离层异常,其时间梯度异常都具有连续性的特点。而由于非电离层原因导致的数据野值,则不具有连续性的特点。基于这一特点,提出了以上区分野值和电离层异常的检测方法。为证明检测结果的正确性,用单站的定位误差结果进行复核。
3. 定位误差验证
野值检测的难点之一是判别异常值产生原因是否为电离层因素。一般认为,双频观测值可消除电离层异常对定位结果的影响。因此,可利用单双频定位误差结果来验证异常产生的原因。电离层异常会引起单频定位结果异常,而不影响双频定位结果;非电离层因素的“野值”则会引起单双频定位结果都出现异常。
RTKLIB是GNSS标准/精密定位开源程序包,由日本东京海洋大学开发,可支持多种GNSS定位模式[19]。本文利用RTKLIB的后处理单点定位模式对观测数据进行单双频定位解算。单频定位采用广播模型进行电离层改正。
对单双频定位结果进行对比,可验证异常产生的原因。图 8~图 10分别是2016年3月26日YNMZ站GPS时1~1.25 h(以下未作说明均为GPS时),YNMZ站14~14.25 h,QHTR站1~1.25 h 3个方向的单双频定位误差结果。野值检测结果表明,QHTR站由于数据出现较多野值而被剔除;YNMZ站在13~16 h发生了电离层异常。对比图 8未发生异常时段的单双频定位结果,图 9的双频定位结果正常,而单频定位出现较大误差,证明了YNMZ站该时段由于电离层异常引起数据异常的结论。而图 10的单双频定位结果都出现了较大误差,可说明QHTR站由于非电离层因素引发了异常。通过比较单双频定位误差结果,验证了野值检测结果的合理性。
4. 野值检测结果
采用上述方法对2016年不同月份的12天电离层垂直延迟数据进行野值检测,并统计剔除站数,如表 1所示。结果显示,在中国陆态网200多个站中,野值检测所剔除的参考站数目在3~10个,对电离层延迟值数据分布的影响在可接受范围内。
表 1 野值剔除结果Table 1. Outlier elimination results日期 剔除站数 2016-01-08 8 2016-01-17 5 2016-02-22 3 2016-02-28 5 2016-03-06 9 2016-03-13 9 2016-03-26 9 2016-04-09 5 2016-06-18 10 2016-10-13 10 2016-12-16 6 2016-12-25 3 在12天中剔除次数超过2次的参考站如图 11所示。对于某些站,多天的野值检测结果呈现一致性。例如,GSTS站在每天的野值检测中都被判别为剔除站,可能由于其接收机长期故障或者是受到干扰,在后续电离层异常建模研究中,将不采用此站的数据。
5. 结论
1) 野值检测方法可根据电离层异常和野值的不同时空相关特性有效检测出野值,并将其剔除,以保证电离层异常数据的可靠性。
2) 利用电离层垂直延迟时间梯度不仅可以筛选野值,还可以区分不同类型的电离层异常现象(电离层等离子体泡或电子密度增强),后续研究电离层异常时可借鉴该方法。
3) 利用单双频定位误差结果可判别异常值产生的原因是否是电离层因素,也可验证野值检测结果的正确性。
4) 在中国陆态网200多个站中,野值检测所剔除的参考站数目在3~10个,对电离层延迟值数据分布的影响在可接受范围内。
本文野值检测方法将用于中国区域电离层异常研究数据预处理阶段,以保障电离层延迟值数据的真实可靠性。为适应不同接收机以及不同数据要求的场景,仍需适当调整野值检测门限参数。
致谢: 感谢“中国大陆构造环境监测网络”(简称中国陆态网)提供的数据支持。 -
表 1 野值剔除结果
Table 1. Outlier elimination results
日期 剔除站数 2016-01-08 8 2016-01-17 5 2016-02-22 3 2016-02-28 5 2016-03-06 9 2016-03-13 9 2016-03-26 9 2016-04-09 5 2016-06-18 10 2016-10-13 10 2016-12-16 6 2016-12-25 3 -
[1] 张其东.基于Kriging法的电离层网格校正技术研究[D].北京: 北京航空航天大学, 2016: 7-17.ZHANG Q D.Research on the technology of grid ionospheric correction based on Kriging[D].Beijing: Beihang University, 2016: 7-17(in Chinese). [2] ICAO APANPIRG/27.SBAS safety assessment guidance related to anomalous ionospheric conditions, Edition 1.0[EB/OL].(2016-07-27)[2019-04-02]. [3] 李海鹏.基于中国区域的电离层网格校正技术研究[D].北京: 北京航空航天大学, 2018: 55-57.LI H P.Research on ionospheric grid correction technology based on China area[D].Beijing: Beihang University, 2018: 55-57(in Chinese). [4] 孟范伟, 郭英, 刘振.电离层异常原因分析与建模[J].测绘与空间地理信息, 2014, 37(1):215-218. doi: 10.3969/j.issn.1672-5867.2014.01.068MENG F W, GUO Y, LIU Z.The cause analysis and modeling of ionosphere anomalies[J].Geomatics & Spatial Information Technology, 2014, 37(1):215-218(in Chinese). doi: 10.3969/j.issn.1672-5867.2014.01.068 [5] 李春斌.中低纬电离层异常扰动统计分析[D].西安: 西安电子科技大学, 2012: 9-17.LI C B.Statistical study of ionosphere mid/low-latitude anomalous disturbances[D].Xi'an: Xidian University, 2012: 9-17(in Chinese). [6] 孙树计, 陈春, 奚迪龙.中国西部及周边区域电离层暴的统计特性[J].装备环境工程, 2008, 5(1):36-39. doi: 10.3969/j.issn.1672-9242.2008.01.009SUN S J, CHEN C, XI D L.Statistic characteristics of ionospheric storms in west China and the adjacency[J].Equipment Environmental Engineering, 2008, 5(1):36-39(in Chinese). doi: 10.3969/j.issn.1672-9242.2008.01.009 [7] 徐彬, 刘钝.电离层闪烁对GNSS的影响[J].全球定位系统, 2011, 36(3):5-8. doi: 10.3969/j.issn.1008-9268.2011.03.003XU B, LIU D.Ionospheric scintillation effects on GNSS[J].GNSS World of China, 2011, 36(3):5-8(in Chinese). doi: 10.3969/j.issn.1008-9268.2011.03.003 [8] SPARKS L, BLANCH J, PANDYA N.Estimating ionospheric delay using kriging:1.Methodology[J].Radio Science, 2011, 46:RS0D21. [9] SPARKS L, BLANCH J, PANDYA N.Estimating ionospheric delay using Kriging:2.Impact on satellite-based augmentation system availability[J].Radio Science, 2011, 46:RS0D22. [10] SPARKS L, ALTSHULER E.Ionospheric storms of solar cycle 24 and their impact on the WAAS ionospheric threat model[C]//Proceedings of the 29th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS+2016).Washington, D.C.: Inst Navigation, 2016: 682-690. [11] SPARKS L, ALTSHULER E.Improving WAAS availability along the coast of California[C]//Proceedings of the 27th International Technical Meeting of the Satellite Division of the Institute of Navigation(ION GNSS+2014).Washington, D.C.: Inst Navigation, 2014: 3299-3311. [12] 袁运斌, 霍星亮, 张宝成.近年来我国GNSS电离层延迟精确建模及修正研究进展[J].测绘学报, 2017, 46(10):1364-1378. doi: 10.11947/j.AGCS.2017.20170349YUAN Y B, HUO X L, ZHANG B C.Research progress of precise models and correction for GNSS ionospheric delay in China over recent years[J].Acta Geodaetica et Cartographica Sinica, 2017, 46(10):1364-1378(in Chinese). doi: 10.11947/j.AGCS.2017.20170349 [13] 袁运斌, 李子申, 王宁波, 等.基于拟合推估的中国区域电离层延迟精确建模方法[J].导航定位学报, 2015, 3(3):49-55.YUAN Y B, LI Z S, WANG N B, et al.Precise modeling of ionospheric delay over China region based on collocation[J].Journal of Navigation and Positioning, 2015, 3(3):49-55(in Chinese). [14] 刘长建.GNSS电离层建模方法与质量控制研究[D].郑州: 中国人民解放军信息工程大学, 2011: 7-25.LIU C J.Study on modeling method and model quality control of ionosphere based on GNSS[D].Zhengzhou: PLA Information Engineering University, 2011: 7-25(in Chinese). [15] WILSON B D, YINGER C H, FEESS W A, et al.New and improved the broadcast interfrequency biases[J].GPS World, 1999, 10(9):56-66. [16] DACH R, HUGENTOBLER U, FRIDEZ P, et al.Bernese GPS Software Version 5.0[EB/OL].(2011-05-20)[2019-04-02]. [17] JIN R, JIN S G, FENG G P.M_DCB:Matlab code for estimating GNSS satellite and receiver differential code biases[J].GPS Solutions, 2012, 16(4):541-548. doi: 10.1007/s10291-012-0279-3 [18] TAKESHI F, TOSHIAKI T.Investigation and modeling of effects of ionospheric plasma bubbles on GBAS availability[C]//Proceedings of the 27th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS+2014).Washington, D.C.: Inst Navigation, 2014: 835-840. [19] 朱李忠, 吴文会, 王连仲, 等.RTKLIB后处理模块定位精度分析及可用性探讨[J].测绘与空间地理信息, 2018, 41(6):41-43. doi: 10.3969/j.issn.1672-5867.2018.06.012ZHU L Z, WU W H, WANG L Z, et al.RTKLIB post-processing module positioning accuracy analysis and availability discussion[J].Geomatics & Spatial Information Technology, 2018, 41(6):41-43(in Chinese). doi: 10.3969/j.issn.1672-5867.2018.06.012 期刊类型引用(3)
1. 何金花,李伟伟. 北斗三代SBAS在抗干扰天线下的实现及验证. 现代导航. 2024(05): 324-329 . 百度学术
2. 张振军. 基于数据筛选的无人机测绘数据异常检测. 西华大学学报(自然科学版). 2022(04): 66-71 . 百度学术
3. 冯乔. 基于超球面支持向量机的传感器网络数据异常检测分析. 微型电脑应用. 2020(10): 174-176 . 百度学术
其他类型引用(2)
-