
摘要:计算位于不同纬度的3个MGEX测站在不同太阳活动及地磁场活动情况下GPS、GLONASS、Beidou和Galileo卫星导航定位系统的电离层高阶项(HOI)延迟,分析高阶项延迟对多系统PPP结果的影响,总结不同太阳活动和地磁活动对PPP的影响。结果表明,电离层延迟随纬度的增加而减小;太阳活动是影响电离层高阶项延迟的主要因素,在太阳活动指数较低时,二阶项延迟不超过2cm,三阶项延迟约为1mm;当太阳指数较高时,二阶项延迟超过2.5cm,三阶项延迟可达到5mm;当磁暴发生时,高阶项延迟会存在一定程度的增长。将电离层高阶项延迟改正后,太阳活动指数较低时点位精度提升效果明显。
电离层延迟是卫星导航定位中的一个重要的误差来源,采用组合双频观测值可消除电离层一阶项延迟的影响,但剩余的高阶项对定位的影响仍然可达2~4cm。国内外学者对GPS观测值中电离层高阶项延迟进行了大量研究并取得丰富成果[1,2,3,4,5]。
在多系统精密单点定位中,各系统的观测值都需进行电离层高阶项改正。本文通过计算2014年与2017年不同纬度的3个MGEX测站中GPS、GLONASS、Galileo以及Beidou(BDS)卫星导航定位系统观测值的电离层高阶项延迟,分析电离层高阶项延迟对多系统精密单点定位结果的影响。
1、电离层延迟计算原理
Li载波频率上顾及电离层延迟的伪距和载波相位观测方程可以表示为:
式中,PRi和Phi分别为伪距和载波相位观测值,其二、三阶电离层延迟可表示为[2,3]:
式中,TEC为信号传播路径上的总电子含量,本文使用全球电离层图(GIM)进行估计;A≅80.6m3/s2,e=1.60218×10-19C,me=9.10939×10-31kg;B可采用国际地磁参考场IGRF12计算得到[6];θ为穿刺点处磁场强度与信号传播方向的夹角;比例因子δ近似取值0.66[7];Ne,max为电离层单层模型高度处最大电子密度[7]。
2、实验分析
2.1实验数据与方法
实验选取分布于不同纬度的3个MGEX测站(KIRU、VILL、MAL2),在不同的太阳活动和地磁活动下计算HOI,其中太阳活动的相对高低通过太阳射电通量(solarradioflux)指数F10.7进行判断,磁暴等级通过Dst指数以及Ap指数进行判断[8]。表1为实验数据选取的时间段与相应的太阳活动指数和地磁指数的分布区间,其中Dst取当天最小值,F10.7取当天20∶00(UT)数值。
为了探讨电离层高阶项延迟对PPP的影响,对高阶项改正前后的观测值使用PPPH软件按相同的解算策略进行对比处理[9],最终求得高阶项改正前后的定位结果在N、E、U方向上的绝对值互差,解算过程中所需数据从IGS中获取。
表1实验数据信息
2.2高阶项延迟计算结果与分析
下文将分别展示测站KIRU、VILL、MAL2观测到的数据中GPS的L1载波、GLONASS的G1载波、Galileo的E1载波和BDS的B1载波上电离层高阶项延迟结果。
2.2.1低太阳活动、无磁暴发生情况下HOI
图1(a)为2017-05-09从GIM中获取的单日平均VTEC(GIM各格网点的单日平均VTEC,以下简称单日平均VTEC),从图中可以看出,该天单日平均VTEC不超过25TECu。图2为由2017-05-09KIRU、VILL、MAL2站4个系统的观测数据计算的电离层二阶项延迟结果,从图中可以看出,在4个系统中,二阶项延迟最大值不超过15mm。
图1单日平均VTEC全球分布
图22017-05-09电离层二阶项延迟
2.2.2低太阳活动、有磁暴发生情况下HOI
图1(b)为2017-09-13单日平均VTEC分布图,从图中可以看出,该天单日平均VTEC不超过30TECu。图3为电离层二阶项延迟结果,在太阳活动较低时,VTEC的变化使各个系统二阶项延迟最大值在磁暴发生时略大于无磁暴情况。由于在低太阳活动情况下,三阶项延迟最大值约为1mm,数值较小,故不展示三阶项。
图32017-09-13电离层二阶项延迟
2.2.3高太阳活动、无磁暴发生情况下HOI
图1(c)为2014-09-11单日平均VTEC分布图,从图中可以看出,该天单日平均VTEC低于60TECu,VTEC随纬度增加而减小。相较于图1(a),由于太阳活动增强,单日平均VTEC最大值增长超过1倍。图4为该天的二阶项延迟结果,各个系统的电离层二阶项延迟最大值均超过25mm。从图5可以看出,三阶项延迟的最大值达到5mm。
图42014-09-11电离层二阶项延迟
图52014-09-11电离层三阶项延迟
2.2.4高太阳活动、有磁暴发生情况下HOI
图1(d)为2014-02-19单日平均VTEC分布图,相较于图1(b),由于太阳活动增强,该天单日平均VTEC增长也超过1倍,最大值为70.9TECu。图6为该天二阶项延迟结果,除Galileo系统因卫星数量太少无法判别外,其他3个系统的二阶项延迟的最大值均超过30mm。从图7可以看出,三阶项延迟中最大值超过5mm。相较于2014-09-11的结果,二、三阶项延迟均有一定增长。
图62014-02-19电离层二阶项延迟
图72014-02-19电离层三阶项延迟
2.3PPP结果与分析
表2为不同时间段的统计结果,ΔN、ΔE、ΔU分别为高阶项延迟改正前后的定位结果在N、E、U方向上的绝对值互差;ΔP为高阶项延迟改正前后的3D位置绝对值差值;“天数”为各测站总共统计的天数,“≥0天数”为N、E、U坐标值在同一方向上的绝对值差值大于等于0的天数;“提升(%)”为“≥0天数”除以“天数”的百分比;“Max(mm)”为在相应统计期间内以mm为单位的最大值;“Min(mm)”为在相应统计期间内以mm为单位的最小值。图8为各时间段内ΔN、ΔE、ΔU、ΔP(图中dN、dE、dU、dP),其中(a)~(d)为KIRU站的结果,(e)~(h)为VILL站的结果,(i)~(l)为MAL2站的结果,位于同一列表示时间段相同。
2.3.1在低太阳活动、无磁暴发生情况下的坐标结果
表2中2017年第129~135天的统计结果表明,各测站在高阶项延迟改正后的定位精度均得到提升;在不同方向上,N方向定位精度提升效果最好,U、E方向波动较大。从图8可以看出,各测站在各个方向上由高阶项延迟引起的坐标绝对值差值基本都在1mm以内,MAL2站在该时间段的最大值为1.1mm。
2.3.2在低太阳活动、有磁暴发生情况下的坐标结果
表2中2017第250~256天的统计结果表明,各测站在高阶项延迟改正后的定位精度均有所提升,E方向定位精度提升效果最好。图8中该时间段的结果显示,在太阳活动较低时期,相比于无磁暴发生的情况,磁暴发生时的坐标绝对值差值存在0~1mm的变化,其中低纬度的MAL2站变化更为明显。
2.3.3在高太阳活动、无磁暴发生情况下的坐标结果
表2点位结果统计
图8电离层高阶项改正前后的坐标绝对值差值
表2中2014年第248~254天的统计结果显示,较太阳活动较低时,KIRU与VILL测站的定位精度更差,其原因可能为太阳活动增强时,电离层TEC变化较大,从GIM产品计算的TEC与这2个测站上空电离层中实际电子含量相差较大,从而导致误差变大[10]。从图8可以看出,由于太阳活动增强,各方向的差值也随之变大,最大值超过3mm,这是由于太阳活动增强使高阶项延迟增大。
2.3.4在高太阳活动、有磁暴发生情况下的坐标结果
表2和图8中2014年第49~55天的统计结果显示,在太阳活动较高时期,磁暴发生时各测站精度提升效果存在较大变化,同时坐标绝对值差值也出现比太阳活动较低时更大的变化。由高阶项延迟计算公式可以推断,造成这种现象的原因可能为在太阳活动较高时,电离层中总电子含量较高,同等级的磁场变化会引起比太阳活动较低时更大的高阶项延迟变化。在3个方向中N方向受影响更大,这可能是由于南北方向电子含量变化更加剧烈[1]。
3、结语
本文分别在不同太阳活动和地磁活动下计算GPS的L1载波、GLONASS的G1载波、Galileo的E1载波和BDS的B1载波电离层高阶项延迟,通过对不同纬度的MGEX测站在不同太阳活动和地磁活动条件下的数据进行PPP处理,分析电离层高阶项延迟改正后的观测值对定位精度的影响。
在不同太阳活动和地磁活动下的电离层高阶项延迟数据显示,太阳活动是影响电离层高阶项延迟的主要因素,在太阳活动指数较低时,二阶项延迟不超过2cm,三阶项延迟约为1mm;当太阳活动指数较高时,二阶项延迟超过2.5cm,三阶项延迟可达到5mm;当磁暴发生时,高阶项延迟都存在一定程度的增长,定位结果会产生0~1mm的变化;在不同纬度上,高阶项延迟随纬度增加而减小。坐标绝对值差值和高阶项延迟的变化存在相同的趋势,将电离层高阶项延迟改正后,坐标绝对值差值最大超过3mm;在太阳活动指数较低时,对电离层高阶项延迟进行改正后的点位精度提升效果更好;在不同方向上,N方向的定位结果受磁暴和太阳活动的影响更大。
综上所述,电离层高阶项延迟是卫星导航定位中的一个重要的误差来源,在进行高精度定位时,应充分考虑电离层高阶项延迟对定位的影响。
参考文献:
[1]朱春春,李征航,屈小川,等.高阶电离层延迟对GPS双差观测值和基线向量的影响[J].大地测量与地球动力学,2015,35(1):81-86.
[5]周要宗,匡翠林,窦邵华.高阶电离层延迟对中国区域双差定位影响研究[J].大地测量与地球动力学,2019,39(5):506-510.
[8]徐文耀.地磁活动指数的过去、现在和未来[J].地球物理学进展,2009,24(3):830-841.
谢名强,徐佳佳,周厚香,李淑慧,徐泽远,王子钦.电离层高阶项延迟对多系统PPP的影响[J].大地测量与地球动力学,2020,40(07):751-755.
基金:国家自然科学基金(41574011).
分享:
土壤磁学是将现代岩石磁学理论、技术方法应用到土壤学而产生的分支学科,通过研究土壤中磁性矿物的类型、含量及颗粒磁畴状态等特征,探索其在环境系统中的溶解、迁移和转化机制,利用物质在磁性特征上的联系及其反映的环境内涵,进而研究不同时空尺度的环境问题、环境过程和作用机制.随着中国黄土高原黄土/古土壤的磁学研究日趋成熟,南方土壤的磁性研究逐渐受到学者重视.
2020-09-05由于电离层层析反演存在着不适定问题,因此许多国内外学者提出了各种方法来解决。文献[1,2]针对GPS电离层层析反演中的不适定问题,先后提出了多种层析算法,取得了较好的反演结果。虽然算法的优化改进可在一定程度上解决不适定的问题,但并不能从根本上消除该问题所带来的负面影响。目前,越来越多的学者通过增加数据源的思路来解决该问题。
2020-07-09在太阳风-磁层-电离层耦合系统中,随着太阳风和行星际磁场的变化,近地空间的电流系统也会发生改变.在高纬地区,电离层电流的快速变化会引起地球表面的磁场发生扰动.地磁扰动是由于日地系统中太阳风能量、动量和物质输运至地球并与近地空间环境相互作用产生的一种地磁活动现象,地磁扰动主要包括磁暴、地磁亚暴和地磁脉动等变化磁场表征现象。
2020-07-09通过实际干扰数据分析可以看出,核自适应滤波算法作为一种非线性滤波方法,能通过实时调节滤波参数从而抑制干扰,为一种完全自适应的无需人工干预的滤波算法。相比于其他算法,该算法对于地磁数据的处理更具有适用性。本文首次将核自适应滤波算法应用于地磁干扰抑制研究,以期为地磁干扰抑制提供新方法。
2020-07-09在多系统精密单点定位中,各系统的观测值都需进行电离层高阶项改正。本文通过计算2014年与2017年不同纬度的3个MGEX测站中GPS、GLONASS、Galileo以及Beidou(BDS)卫星导航定位系统观测值的电离层高阶项延迟,分析电离层高阶项延迟对多系统精密单点定位结果的影响。
2020-07-08中国幅员辽阔,北方多地地磁台站位于高纬度、高海拔地区,冬季温度低,可能达不到地磁等观测仪器正常工作所需环境温度,影响观测数据质量,台站人员应实时监测并及时处理观测数据异常现象,应具备相应的职业素养,明确地磁台站记录数据干扰类型和表现形态,能快速准确地进行数据预处理,确保地磁观测仪器正常工作。
2020-07-06干扰使得观测数据精度降低,对观测资料造成很多大的影响,故在日常观测中对各种干扰的识别尤为重要。本文选取2015年1月~2018年12月的观测数据,分别统计观测环境因素、观测系统影响及人为干扰等各类干扰情况,分析各类干扰对地磁观测数据的影响特征,总结得出都兰地震台地磁FHD观测数据的干扰因素及干扰特征,并提出各类干扰的应对措施,以提高数据的观测质量。
2020-07-06固原4.6级地震是1998年7月29日海原4.9级地震后西海固地区发生的最大一次地震,也是海原4.9级地震以来西海固地区发生的唯一4级以上地震。本文主要对固原4.6级地震发生前的异常和预测指标应用等进行回顾与总结,以期为宁夏地区中等有感及破坏性地震的预测和震情监视跟踪工作提供借鉴。
2020-07-06通过相关方法将地磁数据中的地磁异常提取出来,再对提取出的异常数据进行相关分析。目前,地震预测以观测事实归纳总结为主,辅以物理模式演绎法,将一些地磁场影响因素通过数学物理建模进行定量的扣除。本文中通过相关系数分析、回归分析,确定温度对温泉地磁场子夜均值影响的数学表达式,从而剔除温度对地磁场数据分析中的影响。
2020-07-06在勘探过程中,钻探所得岩心是研究储层裂缝最基础也是最佳的工具,其最能直观反映储层的裂缝形态[8]。如果获得了定向煤心,便可通过CT扫描等方法来对煤样内部空间裂缝进行观测和定量描述,借此来确定煤层裂缝的空间展布,从而进一步服务于煤层气田开发。本文结合了油气勘探中岩心定向方法,提出了将古地磁岩心定向技术用于煤层气勘探开发中的新思路。
2020-07-06我要评论
期刊名称:现代物理知识
期刊人气:548
主管单位:中国科学院
主办单位:中国科学院高能物理研究所
出版地方:北京
专业分类:科学
国际刊号:1001-0610
国内刊号:11-2441/O3
邮发代号:2-824
创刊时间:1976年
发行周期:双月刊
期刊开本:大16开
见刊时间:1-3个月
影响因子:0.735
影响因子:0.645
影响因子:1.369
影响因子:0.874
影响因子:0.385
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!