中国宇航学会深空探测技术专业委员会会刊

高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于GNSS加密网的VLBI电离层时延修正方法

周伟莉 宋淑丽 李培佳 张志斌 黄超 黄勇 王广利

周伟莉, 宋淑丽, 李培佳, 张志斌, 黄超, 黄勇, 王广利. 基于GNSS加密网的VLBI电离层时延修正方法[J]. 深空探测学报(中英文), 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
引用本文: 周伟莉, 宋淑丽, 李培佳, 张志斌, 黄超, 黄勇, 王广利. 基于GNSS加密网的VLBI电离层时延修正方法[J]. 深空探测学报(中英文), 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
ZHOU Weili, SONG Shuli, LI Peijia, ZHANG Zhibin, HUANG Chao, HUANG Yong, WANG Guangli. Ionospheric TEC Correction for VLBI Based on GNSS Density Network[J]. Journal of Deep Space Exploration, 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
Citation: ZHOU Weili, SONG Shuli, LI Peijia, ZHANG Zhibin, HUANG Chao, HUANG Yong, WANG Guangli. Ionospheric TEC Correction for VLBI Based on GNSS Density Network[J]. Journal of Deep Space Exploration, 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056

基于GNSS加密网的VLBI电离层时延修正方法

doi: 10.15982/j.issn.2095-7777.2020.20200056
基金项目: 国家重点研发计划资助项目(2016YFB0501503);国家自然科学基金资助项目(41730109)
详细信息
    作者简介:

    周伟莉(1985– ),女,高级工程师,主要研究方向:GNSS高精度数据处理及应用。通讯地址:上海市南丹路80号中国科学院上海天文台(200030)电话:(021)34775273 E-mail:wlzhou@shao.ac.cn

  • ●The processing strategy of GNSS nearby CVN stations is designed. ●The ionospheric delay correction method based on GNSS inversion is given, which is called SHAO model. and can be used in the VLBI orbit determination. ●When the ionosphere is in the active period, the orbit determination residual based on the SHAO model is 2~3 ns higher than that based on the CODE model in the S-band.
  • 中图分类号: P135,P164

Ionospheric TEC Correction for VLBI Based on GNSS Density Network

  • 摘要: 基于全球导航卫星系统(Global Navigation Satellite System,GNSS)反演的电离层模型是甚长基线干涉(Very Long Baseline Interferometry,VLBI)网(Chinese VLBI Network,CVN)电离层时延修正的有效方法。由于国际GNSS服务组织(International GNSS Services,IGS)跟踪网在中国区域较稀疏,针对探月工程高精度需求,提出了加密GNSS跟踪网的处理策略,构建了全球电离层精修正模型(ShangHai Astronomical Observatory,SHAO),并与欧洲定轨中心发布的全球电离层修正模型(Center for Orbit Determination in Europe,CODE)作对比。结果表明:在电离层平静期,采用两种模型计算的CVN站电离层时延相对均方根误差(Root Mean Square Error,RMSE)最大1.23 ns;在电离层活跃期,CVN站的电离层时延RMSE最大达16.80 ns。在电离层处于活跃期时,基于SHAO模型的VLBI定轨残差相比CODE模型在S1波段由6.07 ns提高到3.33 ns,S2波段由9.10 ns提高到7.07 ns,X0波段由0.70 ns 提高到0.60 ns。
    Highlights
    ●The processing strategy of GNSS nearby CVN stations is designed. ●The ionospheric delay correction method based on GNSS inversion is given, which is called SHAO model. and can be used in the VLBI orbit determination. ●When the ionosphere is in the active period, the orbit determination residual based on the SHAO model is 2~3 ns higher than that based on the CODE model in the S-band.
  • 图  1  IGS跟踪站全球分布图

    Fig.  1  Global distribution of IGS stations

    图  2  CMONOC(红色圆点)和CVN站(浅蓝色五角星)分布图

    Fig.  2  Locations of CMONOC(red dot)and CVN (blue pentagram)stations

    图  3  s8531a观测对应两种模型下CVN站电离层时延和时延率的比较,绿色表示两者差值

    Fig.  3  Comparisons of ionosphere delay &rate at CVN stations using SHAO and CODE model on s8531a,the difference of two presents in green

    图  4  s4a24a观测对应两种模型下CVN站电离层时延和时延率的比较,绿色表示两者差值

    Fig.  4  Comparisons of ionosphere delay &rate of CVN stations using SHAO model and the CODE on s4a24a,the difference of two presents in green

    图  5  S1波段两种模型CVN站基线方向定轨残差

    Fig.  5  The difference of orbit determination residuals on S1 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

    图  6  S2波段两种模型CVN站基线方向定轨残差

    Fig.  6  The difference of orbit determination residuals on S2 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

    图  7  X0波段两种模型CVN站基线方向定轨残差

    Fig.  7  The difference of orbit determination residuals on X0 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

    图  8  X1波段两种模型CVN站基线方向定轨残差

    Fig.  8  The difference of orbit determination residuals on X1 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

    表  1  电离层平静期两种模型计算的CVN站电离层时延精度比较

    Table  1  Comparison of ionosphere delay at CVN stations using SHAO or CODE model with the quiet period

    ns  
    时延BiasRMSE
    弧段s8529as8530as8531as8529as8530as8531a
    BJ1.150.700.391.230.840.54
    KM1.030.730.361.110.890.50
    UR0.871.060.911.021.211.10
    TM0.790.510.290.970.710.37
     注:本文电离层时延和时延率均基于2 GHz。
    下载: 导出CSV

    表  2  电离层平静期两种模型计算的CVN站电离层时延率精度比较

    Table  2  Comparison of ionosphere delay rate at CVN stations using SHAO or CODE model with the quiet period

    (ps·s–1)  
    时延BiasRMSE
    弧段s8529as8530as8531as8529as8530as8531a
    BJ0.190.210.230.210.240.27
    KM0.140.140.140.230.190.20
    UR0.110.130.140.220.150.16
    TM0.160.160.170.130.200.20
    下载: 导出CSV

    表  3  电离层活跃期两种模型计算的CVN站电离层时延比较精度

    Table  3  Comparison of ionosphere delay at CVN stations using SHAO or CODE model with the active period

    ns  
    时延BiasRMSE
    弧段s4a24as4a25as4a24as4a25a
    BJ6.733.197.403.56
    KM6.8114.3410.1116.80
    UR2.834.183.084.26
    TM8.627.739.7811.10
    下载: 导出CSV

    表  4  电离层活跃期两种模型计算的CVN站电离层时延率比较精度

    Table  4  Comparison of ionosphere delay rate at CVN stations using SHAO or CODE model with the active period

    (ps·s–1)  
    时延BiasRMSE
    弧段s4a24as4a25as4a24as4a25a
    BJ1.060.691.310.98
    KM2.081.022.851.17
    UR0.500.220.650.30
    TM1.981.232.801.70
    下载: 导出CSV

    表  5  采用两种模型分别计算不同波段的定轨残差

    Table  5  The residuals of orbit determination on different bands using SHAO or CODE

    ns  
    定轨残差s4a24a
    波段频率/HzS12210S22234X08440X18460
    SHAO3.337.070.602.23
    CODE6.079.100.702.23
    下载: 导出CSV
  • [1] 欧阳自远. 我国月球探测的总体科学目标与发展战略[J]. 地球科学进展,2004,19(3):355-357.

    OUYANG Z Y. Scientific objectives of Chinese lunar exploration project and development strategy[J]. Advance in Earth Sciences,2004,19(3):355-357.
    [2] BARBOSA R C. China’s Chang’e -3 and Jade Rabbit Duo Land on the Moon[EB/OL]. (2020-6-15).http://NASAspaceflight.com.
    [3] 董光亮,李海涛,郝万宏,等. 中国深空测控系统建设与技术发展[J]. 深空探测学报,2018,5(2):99-114.

    DONG G L,LI H T,HAO W H,et al. Development and Future of China’s Deep Space TT&C System[J]. Journal of Deep Space Exploration,2018,5(2):99-114.
    [4] 洪晓瑜, 张秀忠, 郑为民, 等. VLBI技术在中国探月工程的应用[J]. 深空探测学报(中英文), 2020, 7(4): 321-331.

    HONG X Y, ZHANG X Z, ZHENG W M, et al. Application of VLBI technology in ChinaLunar exploration project[J]. Journal of Deep Space Exploration, 2020, 7(4): 321-331.
    [5] LANYI G. Total ionospheric electron content calibration using series gps satellite data[J]. Telecommunications and Data Acquisition Report,SEE N86-28275 19,1986,32:1-12.
    [6] 李子申,王宁波,李敏,等. 国际GNSS服务组织全球电离层TEC格网精度评估与分析[J]. 地球物理学报,2017,60(10):3718-33729. doi:  10.6038/cjg20171003

    LI Z S,WANG N B,LI M,et al. Evaluation and analysis of the global ionospheric TEC map in the frame of international GNSS services[J]. Chinese Journal of Geophysics,2017,60(10):3718-33729. doi:  10.6038/cjg20171003
    [7] 欧洲定轨中心全球电离层修正模型[DB/OL].(2020-6-15). http://ftp.aiub.unibe.ch/.
    [8] IGS. IGS data[DB/OL].(2020-6-15).http://www.igs.org/.
    [9] CMONOC. CMONOC data[DB/OL].(2020-6-15). http://www.neiscn.org/.
    [10] ZHOU W L,SONG S L,CHEN Q M,et al. Determination of nighttime VTEC average in the Klobuchar ionospheric delay model[J]. Geodesy and Geodynamics,2018,9:175-182.
    [11] SONG S L,CHENG N,XIE H,et al. Adaptability of GPS/BDS broadcast ionospheric models to solar activities[J]. Advance in Space Research,2019,63:2867-2876.
    [12] CHENG N,SONG S L,XIE H. Investigation of Solar flares impact on GPS/BDS/GALILEO broadcast ionospheric models[J]. Radio Science,2019,54:91-103.
    [13] YANG Z, MORTON J Y, ZAKHARENKOVA I, et al. Global view of ionospheric disturbances impacts on kinematic GPS Positioning solutions during the 2015 St Patrick’s day storm[J]. Journal of Geophysical Research: Space Physics, 2020, 125(7): A02768110.1029/2019JA027681
    [14] 黄勇,李培佳,樊敏,等. 基于地基测量数据的地月L2点探测器轨道确定[J]. 中国科学:物理学 力学 天文学,2018,48(7):079501.

    HUANG Y,LI P J,FAN M,et al. Orbit determination of CE-5T1 in Earth-Moon L2 libration point orbit with ground tracking data[J]. Scientia Sinica Physica,Mechanica & Astronomica,2018,48(7):079501.
  • [1] 黄勇, 李培佳, 胡小工.  VLBI月球定轨和月面定位技术及其应用 . 深空探测学报(中英文), 2020, 7(4): 340-346. doi: 10.15982/j.issn.2095-7777.2020.20200022
    [2] 陈如荣, 张海燕, 金乘进, 高智胜, 朱岩, 岳友岭, 朱凯.  FAST VLBI系统和观测研究 . 深空探测学报(中英文), 2020, 7(2): 136-143. doi: 10.15982/j.issn.2095-7777.2020.20190618001
    [3] 王广利, 洪晓瑜, 刘庆会, 李斌, 马军, 苏彦, 郝龙飞.  “嫦娥4号”高精度VLBI测轨技术 . 深空探测学报(中英文), 2020, 7(4): 332-339. doi: 10.15982/j.issn.2095-7777.2020.20200026
    [4] 郑为民, 张娟, 徐志骏, 刘磊, 童力, 张秀忠.  实时VLBI处理机技术 . 深空探测学报(中英文), 2020, 7(4): 354-361. doi: 10.15982/j.issn.2095-7777.2020.20200027
    [5] 洪晓瑜, 张秀忠, 郑为民, 王广利, 刘庆会, 王娜, 张洪波, 汪敏, 舒逢春.  VLBI技术研究进展及在中国探月工程的应用 . 深空探测学报(中英文), 2020, 7(4): 321-331. doi: 10.15982/j.issn.2095-7777.2020.20200030
    [6] 韩松涛, 谢剑锋, 王美, 李黎.  一种应用于实时深空干涉测量的电离层时延修正方法 . 深空探测学报(中英文), 2019, 6(3): 232-235,300. doi: 10.15982/j.issn.2095-7777.2019.03.005
    [7] 吴伟仁, 刘继忠, 唐玉华, 于登云, 于国斌, 张哲.  中国探月工程 . 深空探测学报(中英文), 2019, 6(5): 405-416. doi: 10.15982/j.issn.2095-7777.2019.05.001
    [8] 段玉虎.  基于波束分离技术的空间时延补偿方法研究 . 深空探测学报(中英文), 2018, 5(2): 124-128,146. doi: 10.15982/j.issn.2095-7777.2018.02.003
    [9] 王美, 陈略, 韩松涛, 路伟涛, 段建锋, 任天鹏, 李黎.  深空测控网干涉测量系统在“鹊桥”任务中的应用分析 . 深空探测学报(中英文), 2018, 5(6): 539-543. doi: 10.15982/j.issn.2095-7777.2018.06.006
    [10] 董光亮, 李海涛, 郝万宏, 王宏, 朱智勇, 石善斌, 樊敏, 周欢, 徐得珍.  中国深空测控系统建设与技术发展 . 深空探测学报(中英文), 2018, 5(2): 99-114. doi: 10.15982/j.issn.2095-7777.2018.02.001
    [11] 刘庆会.  火星探测VLBI测定轨技术 . 深空探测学报(中英文), 2018, 5(5): 435-441. doi: 10.15982/j.issn.2095-7777.2018.05.005
    [12] 任天鹏, 唐歌实, 许柏, 路伟涛, 陈略, 韩松涛, 王美.  面向绕月交会对接的同波束VLBI研究 . 深空探测学报(中英文), 2017, 4(4): 367-372. doi: 10.15982/j.issn.2095-7777.2017.04.009
    [13] 王春锋.  卫星编队自主相对导航与通信一体化系统探讨 . 深空探测学报(中英文), 2017, 4(1): 38-42. doi: 10.15982/j.issn.2095-7777.2017.01.006
    [14] 于登云, 吴学英, 吴伟仁.  我国探月工程技术发展综述 . 深空探测学报(中英文), 2016, 3(4): 307-314. doi: 10.15982/j.issn.2095-7777.2016.04.002
    [15] 刘文俊, 郝万宏, 谌明, 李海涛.  异地天线组阵站间时延差修正技术研究与验证 . 深空探测学报(中英文), 2015, 2(3): 267-271. doi: 10.15982/j.issn.2095-7777.2015.03.013
    [16] 王明远, 平劲松.  基于MARSIS的火星低纬度电离层统计特性研究 . 深空探测学报(中英文), 2014, 1(4): 294-297. doi: 10.15982/j.issn.2095-7777.2014.04.008
    [17] 王震, 王娜, 平劲松.  月球电离层掩星探测研究 . 深空探测学报(中英文), 2014, 1(3): 220-225. doi: 10.15982/j.issn.2095-7777.2014.03.010
    [18] MikhailVasilyevichVASILYEV, EleonoraIvanovnaYAGUDINA.  俄罗斯应用天文研究所月球历表研究现状 . 深空探测学报(中英文), 2014, 1(3): 187-191. doi: 10.15982/j.issn.2095-7777.2014.03.004
    [19] 平劲松, 王明远, 张素君, 简念川, 王震, 鄢建国, 孟桥, 陈从颜, 席禹豪, 唐继斐, 张天翼, 李文潇, 张洪波, 汪敏.  嫦娥系列探月卫星无线电科学实验简介 . 深空探测学报(中英文), 2014, 1(3): 192-199. doi: 10.15982/j.issn.2095-7777.2014.03.005
    [20] 孙靖, 王美, 平劲松.  利用VLBI技术进行深空航天器跟踪的仿真分析 . 深空探测学报(中英文), 2014, 1(3): 226-229. doi: 10.15982/j.issn.2095-7777.2014.03.011
  • 加载中
图(8) / 表 (5)
计量
  • 文章访问数:  19
  • HTML全文浏览量:  4
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-15
  • 修回日期:  2020-08-10
  • 网络出版日期:  2020-11-09
  • 刊出日期:  2020-11-09

基于GNSS加密网的VLBI电离层时延修正方法

doi: 10.15982/j.issn.2095-7777.2020.20200056
    基金项目:  国家重点研发计划资助项目(2016YFB0501503);国家自然科学基金资助项目(41730109)
    作者简介:

    周伟莉(1985– ),女,高级工程师,主要研究方向:GNSS高精度数据处理及应用。通讯地址:上海市南丹路80号中国科学院上海天文台(200030)电话:(021)34775273 E-mail:wlzhou@shao.ac.cn

  • ●The processing strategy of GNSS nearby CVN stations is designed. ●The ionospheric delay correction method based on GNSS inversion is given, which is called SHAO model. and can be used in the VLBI orbit determination. ●When the ionosphere is in the active period, the orbit determination residual based on the SHAO model is 2~3 ns higher than that based on the CODE model in the S-band.
  • 中图分类号: P135,P164

摘要: 基于全球导航卫星系统(Global Navigation Satellite System,GNSS)反演的电离层模型是甚长基线干涉(Very Long Baseline Interferometry,VLBI)网(Chinese VLBI Network,CVN)电离层时延修正的有效方法。由于国际GNSS服务组织(International GNSS Services,IGS)跟踪网在中国区域较稀疏,针对探月工程高精度需求,提出了加密GNSS跟踪网的处理策略,构建了全球电离层精修正模型(ShangHai Astronomical Observatory,SHAO),并与欧洲定轨中心发布的全球电离层修正模型(Center for Orbit Determination in Europe,CODE)作对比。结果表明:在电离层平静期,采用两种模型计算的CVN站电离层时延相对均方根误差(Root Mean Square Error,RMSE)最大1.23 ns;在电离层活跃期,CVN站的电离层时延RMSE最大达16.80 ns。在电离层处于活跃期时,基于SHAO模型的VLBI定轨残差相比CODE模型在S1波段由6.07 ns提高到3.33 ns,S2波段由9.10 ns提高到7.07 ns,X0波段由0.70 ns 提高到0.60 ns。

注释:
1)  ●The processing strategy of GNSS nearby CVN stations is designed. ●The ionospheric delay correction method based on GNSS inversion is given, which is called SHAO model. and can be used in the VLBI orbit determination. ●When the ionosphere is in the active period, the orbit determination residual based on the SHAO model is 2~3 ns higher than that based on the CODE model in the S-band.

English Abstract

周伟莉, 宋淑丽, 李培佳, 张志斌, 黄超, 黄勇, 王广利. 基于GNSS加密网的VLBI电离层时延修正方法[J]. 深空探测学报(中英文), 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
引用本文: 周伟莉, 宋淑丽, 李培佳, 张志斌, 黄超, 黄勇, 王广利. 基于GNSS加密网的VLBI电离层时延修正方法[J]. 深空探测学报(中英文), 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
ZHOU Weili, SONG Shuli, LI Peijia, ZHANG Zhibin, HUANG Chao, HUANG Yong, WANG Guangli. Ionospheric TEC Correction for VLBI Based on GNSS Density Network[J]. Journal of Deep Space Exploration, 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
Citation: ZHOU Weili, SONG Shuli, LI Peijia, ZHANG Zhibin, HUANG Chao, HUANG Yong, WANG Guangli. Ionospheric TEC Correction for VLBI Based on GNSS Density Network[J]. Journal of Deep Space Exploration, 2020, 7(4): 362-370. doi: 10.15982/j.issn.2095-7777.2020.20200056
    • 测轨分系统是中国探月工程和深空探测的重要组成部分,其任务是采用甚长基线干涉测量(Very Long Baseline Interferometry,VLBI)测量技术进行探测器轨道确定[1-4]。我国的VLBI测控网包括北京(BJ)、天马(TM)、昆明(KM)和乌鲁木齐(UR)4个射电望远镜台站和位于上海的1个VLBI数据处理中心。VLBI技术基于最小二乘法拟合相位信息,获取理论时延和时延率,在对探测器进行跟踪测量的过程中,无线电信号穿过中性大气和电离层时会产生时延,从而引起时延误差。中性大气时延修正可以通过模型和全球定位系统(Global Positioning System,GPS)观测数据计算天顶时延,然后采用映射函数投影到探测器观测方向,获取斜路径方向的中性大气时延量,或者采用微波辐射计直接观测斜路径时延。

      全球导航卫星系统(Global Navigation Satellite System,GNSS)和电离层有着紧密的联系,卫星信号在穿过电离层时会受到其折射的影响带来一定的时延误差;利用该折射的影响可通过GNSS反演获得电离层时延量。由于GNSS具有全天候、高精度和易布站优点,国际上主要采用GNSS反演获取的全球电离层修正模型消除电离层时延的影响。早在1986年,Lanyi就利用GPS卫星数据,采用多项式模型构建了区域电离层改正模型[5],从而开启了利用GNSS数据反演电离层模型的新方法。

      国际GNSS服务(International GNSS Services,IGS)下设电离层工作组,包括欧洲定轨中心(Center for Orbit Determination in Europe,CODE)、美国喷气推进实验室(Jet Propulsion Laboratory,JPL)、欧洲航天局(European Space Agency,ESA)等多家机构。他们对全球发布高精度的电离层模型产品,包括:基于球谐函数建立的模型、基于三次样条函数建立的模型和基于三维像素建立的电离层模型等。李子申等对IGS各家电离层模型精度进行了评估分析[6]。从评估结果可以看出,欧洲定轨中心提供的基于球谐函数建立的电离层模型产品精度相对最优。近年来,国际GNSS监测评估系统(international GNSS Monitoring & Assessment System,iGMAS)也在提供电离层模型产品[7]

      从当前电离层模型研究的总体状态来看,要想提高电离层时延修正精度,可以从观测网处理策略和模型算法等方面开展深入研究。针对探月工程VLBI测站的高精度应用需求,主要从观测网处理策略方面开展研究,以提高电离层时延修正精度。目前国际上基于GNSS反演的全球电离层修正模型(CODEZ)采用中国区域的测站作为全球GNSS服务站,但是目前在中国区域的IGS GNSS跟踪网较稀疏,且有时数据中断。为了提高中国区域电离层时延修正精度,即提高VLBI技术定轨精度,采用中国区域GNSS网加密的处理策略,构建了全球电离层时延精修正模型(ShangHai Astronomical Observatory ,SHAO)。由于目前欧洲定轨中心提供的电离层模型精度代表了目前国际先进水平,本文基于探月工程实际任务,对比了SHAO全球电离层时延精修正模型与欧洲定轨中心全球电离层修正模型的精度。

    • 本文采用GNSS观测数据反演构建了全球电离层精修正模型,该模型采用15 × 15阶球谐函数表达式,并且2 hr提供一组球谐系数。

    • 为了便于电离层模型化,通常将整个电离层看作是围绕地球一定高度的薄层,薄层高度300~500 km之间。基于GNSS反演的电离层模型其伪距$ {P}_{i} $和相位$ {\varnothing }_{i} $的观测方程分别为

      $$ {P}_{i}=\rho -\frac{40.3}{{f}_{i}^{2}}TEC+{b}_{i}^{s}-{b}_{ri}+\Delta $$ (1)
      $$ {\varnothing }_{i}=\rho +\frac{40.3}{{f}_{i}^{2}}TEC-{N}_{i}{\lambda }_{i}+{b}_{i}^{s}-{b}_{ri}+\Delta $$ (2)

      其中:${\rm{\rho}} $为接收机至卫星的几何距离,包含卫星到地面几何距离、对流层时延、卫星钟差和接收机钟差及其与频率无关的项等;$ {b}_{i}^{s} $$ {b}_{ri} $分别为频率$ i $的卫星硬件时延和接收机硬件时延;$ {\lambda }_{i} $为对应频率$ i $的波长;$ {N}_{i} $为整周模糊度;$ \Delta $为测量噪声和多路径效应等。

      对GNSS观测数据进行反演,模型时间域使用分段线性函数(Piece-Wise Linear,PWL)进行约束,处理间隔为2 hr,电离层模型一天分为13组,实现分段建模;采用最小二乘法参数估计获取电离层球谐函数表达(SHAO)。全球电离层精修正模型球谐函数表达式为

      $$\begin{split} {{VTEC}} = & \sum\nolimits_{n = 0}^{{n_{\rm max}}} \sum\nolimits_{m = 0}^n {\tilde P_{nm}}\times \\ & {\sin \beta } \left( {{{\tilde C}_{nm}}\cos \left( {m\alpha } \right) + {{\tilde S}_{nm}}\sin \left( {m\alpha } \right)} \right) \end{split}$$ (3)
      $${M}\left( {{z}} \right) = \left( {1 - {\rm{si{n^2}z}/}}\left( {1 + {H_{ion}}/{ R_{\rm E}}} \right) \right)^{ - 1/2}$$ (4)

      其中:VTEC为电离层垂直总电子含量(Vertical Total Electron Content);nm为球谐函数的阶次,${n}_{\rm max}$为球谐函数的最高阶数(本文取15);$ \alpha $为穿刺点经度;$ \beta $为穿刺点纬度;$ {\tilde{P}}_{nm} $为正规化勒让德函数;$ {\tilde{C}}_{nm} $$ {\tilde{S}}_{nm} $为待估球谐函数的系数,即全球电离层精修正模型系数;Mz)表示与高度角z相关的投影函数;$ {H}_{ion} $表示电离层单层高度(本文取450 km);${R}_ { \rm E}$表示地球平均半径;z为卫星高度角。

      基于探月工程探测器星历和CVN测站坐标,获取探测器的高度角和方位角,计算CVN站与探测器方向上的电离层时延和时延率,为VLBI技术定轨消除电离层时延影响。

      本文利用偏差Bias和均方根误差RMSE比较两种模型性能,定义如下:

      $${{Bias}} = \left\langle {Io{n_s} - Io{n_c}} \right\rangle $$ (5)
      $${RMS\!E} = \sqrt {\left\langle {{{\left( {Io{n_s} - Io{n_c}} \right)}^2}} \right\rangle } $$ (6)

      其中:$\left\langle {} \right\rangle $表示给定周期的平均值;$ {Ion}_{s} $表示用SHAO模型计算的电离层时延值;$ {Ion}_{c} $为采用欧洲定轨中心发布的事后全球电离层修正模型[8]计算的电离层时延值。

    • 本文构建的全球电离层时延精修正模型采用的GNSS跟踪站数据包括IGS跟踪站网[9]和中国大陆构造环境监测网络(Crustal Movement Observation Network Of China ,CMONOC)[10]。IGS跟踪站网是一个大规模、高精度的GNSS全球网络,本文选取了全球均匀分布的112个IGS测站的观测资料,其位置信息如图1所示。为了提高CVN站的电离层时延修正精度,采用CMONOC观测数据进行加密,选取了CMONOC的123个测站参与建模,这些GNSS站几乎覆盖整个中国大陆,CMONOC站和4个CVN站位置分布如图2所示。

      图  1  IGS跟踪站全球分布图

      Figure 1.  Global distribution of IGS stations

      图  2  CMONOC(红色圆点)和CVN站(浅蓝色五角星)分布图

      Figure 2.  Locations of CMONOC(red dot)and CVN (blue pentagram)stations

    • 目前,常用的电离层时延修正方法是双频消除法和GNSS反演模型法。对于双频消除法,利用S双频计算电离层改正公式为

      $$ {\tau }_{0}=\frac{{\tau }_{2}{f}_{2}^{2}-{\tau }_{1}{f}_{1}^{2}}{{f}_{2}^{2}-{f}_{1}^{2}} $$ (7)

      其中:$ {f}_{1} $$ {f}_{2} $表示频段的两个频率;$ {\tau }_{1} $$ {\tau }_{2} $是对应两个频率的观测时延值。

      由误差传播定律,电离层时延形式误差为

      $$ {\sigma }_{{\tau }_{0}}=\sqrt{\frac{{f}_{1}^{4}{\sigma }_{{\tau }_{1}}^{2}+{f}_{2}^{4}{\sigma }_{{\tau }_{2}}^{2}}{{\left({f}_{2}^{2}-{f}_{1}^{2}\right)}^{2}}} $$ (8)

      其中:$ {\sigma }_{{\tau }_{1}} $$ {\sigma }_{{\tau }_{2}} $分别为两个频段下的观测时延精度。

      电离层时延的影响与观测频率的平方成反比,如在22 GHz波段,1 TECU(Total Electron Content Unit)对应于约0.8 mm的误差,而在2.2 GHz波段就会达到约8 cm。因此,电离层对低频观测时延的影响较大,高频(10 GHz以上)影响为小量。“嫦娥4号”(CE-4)中继星VLBI信标采用S波段(2 GHz),电离层时延是观测时延的主要误差源。CE-4和“嫦娥5号”再入返回飞行试验(CE-5T1)定轨时延精度要求分别优于5 ns和4 ns;对于S波段,电离层时延修正精度要优于2 ns。 对于S频段,S1和S2信号对应的频率分别为2 210 MHz和2234 MHz。假设两个信标时延的测定精度均为1 ns,则由此可计算得出$ {\sigma }_{{\tau }_{0}} $的值为66 ns,等效于20 m。在实际的探月工程中,由于对定轨具有实时性需求,深空VLBI观测所用单通道带宽以及整个频带跨度无法像测地VLBI一样做宽,且因实时任务无法保证探测器观测信标全部开启,并不适合采用双频法实时消除电离层影响。

      电离层时延与太阳活动、纬度和季节变化紧密相关[11-13],4个CVN测站纬度差异较大,需对4站电离层时延分别进行比对分析。太阳活动存在一个11年周期,根据电离层活跃情况,分电离层平静期和活跃期两种方案进行模型比对。

    • 选取了CE-4中继星任务中2018年5月3次观测s8529a、s8530a和s8531a进行精度分析。根据射电源和探测器星历,计算射电源和探测器的高度角方位角信息,然后结合SHAO和CODE模型的电离层球谐系数文件计算得到对应2 GHz的各CVN站电离层时延和时延率。并按观测弧段分别进行统计,分别给出观测代码s8529a、s8530a和s8531a 3次观测的TM、BJ、KM和UR测站时延和时延率结果,如表1和表2所示。

      表 1  电离层平静期两种模型计算的CVN站电离层时延精度比较

      Table 1.  Comparison of ionosphere delay at CVN stations using SHAO or CODE model with the quiet period

      ns  
      时延BiasRMSE
      弧段s8529as8530as8531as8529as8530as8531a
      BJ1.150.700.391.230.840.54
      KM1.030.730.361.110.890.50
      UR0.871.060.911.021.211.10
      TM0.790.510.290.970.710.37
       注:本文电离层时延和时延率均基于2 GHz。

      表 2  电离层平静期两种模型计算的CVN站电离层时延率精度比较

      Table 2.  Comparison of ionosphere delay rate at CVN stations using SHAO or CODE model with the quiet period

      (ps·s–1)  
      时延BiasRMSE
      弧段s8529as8530as8531as8529as8530as8531a
      BJ0.190.210.230.210.240.27
      KM0.140.140.140.230.190.20
      UR0.110.130.140.220.150.16
      TM0.160.160.170.130.200.20

      图3给出了利用SHAO和CODE模型在电离层平静期频率2 GHz的CVN站电离层时延值。两种模型计算的测站时延均在10.0 ns以内,时延率在0.9 ps/s以内。表1表2分别是电离层平静期采用两种模型分别计算的CVN站电离层时延和时延率比较。从结果可以看出,两种模型计算的时延值RMSE在1.3 ns以内,时延率RMSE在0.5 ps/s以内,不同纬度带各CVN站差异较少。

      图  3  s8531a观测对应两种模型下CVN站电离层时延和时延率的比较,绿色表示两者差值

      Figure 3.  Comparisons of ionosphere delay &rate at CVN stations using SHAO and CODE model on s8531a,the difference of two presents in green

    • 2014年电离层处于活跃期,结合再入返回飞行试验(CE-5T1)任务,选取了2014年10月s4a24a和s4a25a两次观测进行精度分析。图4给出了利用SHAO和CODE模型分别计算的CVN站电离层时延和时延率值。计算方式同上,将SHAO和CODE模型的电离层球谐系数文件分别基于CE-5T1任务时探测器星历得到对应2 GHz的各CVN站电离层时延和时延率。按观测弧段进行统计,得到了s4a24a和s4a25a弧段的测站(TM、BJ、KM和UR)时延和时延率结果。

      图4可知在电离层活跃期,TM和KM测站电离层时延影响尤其明显,电离层时延达到100 ns,BJ和UR测站时延相对较少,但也接近60 ns。这是由于KM和TM测站纬度较低,电离层时延受太阳活动影响较大,BJ和UR站纬度相对较高,受太阳活动影响相对较小导致。表3和表4给出了两种模型计算的电离层时延和时延率偏差和均方根误差。两种模型计算的电离层时延RMSE最大达到16 ns,时延率RMSE达到2.8 ps/s。尤其是相对离赤道更近的KM站,差异比较明显。

      图  4  s4a24a观测对应两种模型下CVN站电离层时延和时延率的比较,绿色表示两者差值

      Figure 4.  Comparisons of ionosphere delay &rate of CVN stations using SHAO model and the CODE on s4a24a,the difference of two presents in green

      表 3  电离层活跃期两种模型计算的CVN站电离层时延比较精度

      Table 3.  Comparison of ionosphere delay at CVN stations using SHAO or CODE model with the active period

      ns  
      时延BiasRMSE
      弧段s4a24as4a25as4a24as4a25a
      BJ6.733.197.403.56
      KM6.8114.3410.1116.80
      UR2.834.183.084.26
      TM8.627.739.7811.10

      表 4  电离层活跃期两种模型计算的CVN站电离层时延率比较精度

      Table 4.  Comparison of ionosphere delay rate at CVN stations using SHAO or CODE model with the active period

      (ps·s–1)  
      时延BiasRMSE
      弧段s4a24as4a25as4a24as4a25a
      BJ1.060.691.310.98
      KM2.081.022.851.17
      UR0.500.220.650.30
      TM1.981.232.801.70
    • 由于电离层平静期两种模型计算的电离层时延相差在1 ns 左右,VLBI技术定轨残差的差异不大,本文只比对电离层活跃期基于两种模型计算的电离层时延对VLBI技术定轨的影响。

      CE5T1探测器于协调世界时(Coordinated Universal Time,UTC)2014年10月24日 08:29:15和UTC 2014年10月25日08:24:43分别进行了第一次和第二次中途修正,在两次中途修正期间多次进行调姿喷气,故选取第一次中途修正前的弧段(2014年10月24日3点20分到7点)进行轨道确定。VLBI技术的具体轨道确定策略见参考文献[14]

      图58给出了S1、S2、X0和X1波段分别采用两种电离层模型计算的探测器4个波段上CVN站基线方向上定轨残差。可见,S波段的定轨残差要明显大于X波段,这是由于电离层时延与频率的平方成反比引起的。表5给出了不同波段探测器定轨残差精度差异,采用SHAO模型比CODE模型定轨残差精度在S1波段从6.1 ns提高到3.3 ns,S2波段从9.1 ns提高到7.1 ns,分别提升了约90%和20%,保证了VLBI测量精度满足指标要求。采用SHAO模型比CODE模型定轨残差精度在X0波段提高0.1 ns,X1波段精度相当。所以,在S波段观测时,基于加密GNSS网构建的全球电离层精修正模型可以有效提高VLBI技术测定轨精度。

      图  5  S1波段两种模型CVN站基线方向定轨残差

      Figure 5.  The difference of orbit determination residuals on S1 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

      图  6  S2波段两种模型CVN站基线方向定轨残差

      Figure 6.  The difference of orbit determination residuals on S2 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

      图  7  X0波段两种模型CVN站基线方向定轨残差

      Figure 7.  The difference of orbit determination residuals on X0 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

      图  8  X1波段两种模型CVN站基线方向定轨残差

      Figure 8.  The difference of orbit determination residuals on X1 band at CVN stations base-line using ionosphere delay calculated by SHAO or CODE

      表 5  采用两种模型分别计算不同波段的定轨残差

      Table 5.  The residuals of orbit determination on different bands using SHAO or CODE

      ns  
      定轨残差s4a24a
      波段频率/HzS12210S22234X08440X18460
      SHAO3.337.070.602.23
      CODE6.079.100.702.23
    • 本文首先介绍了探月工程VLBI测轨分系统中介质误差修正的几种方法,针对探月工程高精度实际需求,基于中国区域GNSS加密网构建了全球电离层精修正模型SHAO。该模型在iGMAS上海天文台分析中心平台开发,结合探月工程CE4中继星任务和CE5T1任务,比较了SHAO模型和CODE模型在不同时期电离层时延和定轨残差方面的差异,得出以下几个结论。

      1)电离层处于活跃期时,两种模型在CVN站上的电离层时延差异RMSE最大可达16.80 ns。相对于CODE模型,利用CMONOC数据加密获取的电离层时延进行VLBI定轨的残差在S1波段由6.07 ns提高到3.33 ns,提升了约90%;在S2波段由9.10 ns提高到7.07 ns,提升了约20%;X0波段由0.70 ns提高到0.60 ns,提升了约14%;X1波段精度相当。所以,在S波段观测时,基于加密GNSS网构建的全球电离层精修正模型可以有效提高VLBI技术测定轨精度。

      2)电离层处于活跃期时,不同纬度各CVN站差异较大。TM和KM测站电离层时延影响尤其明显,电离层时延达到100 ns,BJ和UR测站时延相对较少,接近60 ns。

      3)电离层平静期,采用全球电离层精修正模型SHAO与欧洲定轨中心发布的CODE模型计算的CVN站电离层时延RMSE最大1.23 ns,时延率RMSE优于0.5 ps/s,精度相当。不同纬度各CVN站差异较少。

      致 谢

      感谢CE-4中继星任务和CE-5T1任务VLBI测轨分系统VLBI中心提供的分析数据,感谢iGMAS提供的数据和技术支持。

参考文献 (14)

目录

    /

    返回文章
    返回