CN110764121A - 一种gnss实时全球电离层tec建模方法 - Google Patents
一种gnss实时全球电离层tec建模方法 Download PDFInfo
- Publication number
- CN110764121A CN110764121A CN201911040108.2A CN201911040108A CN110764121A CN 110764121 A CN110764121 A CN 110764121A CN 201911040108 A CN201911040108 A CN 201911040108A CN 110764121 A CN110764121 A CN 110764121A
- Authority
- CN
- China
- Prior art keywords
- time
- real
- model
- global
- tec
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种电离层建模方法,尤其是涉及一种GNSS实时全球电离层TEC建模方法;该方法主要采用了两种数据源建立实时全球电离层TEC模型,其中预报全球电离层TEC模型作为背景场描述全球电离层TEC的整体趋势,GNSS实时观测数据能够反映各个观测站附近电离层TEC的实时变化和局部变化;结合两种数据源既能够保证全球范围反演电离层TEC的可靠性,又能提高GNSS实时观测数据覆盖地区的模型精度。
Description
技术领域
本发明涉及一种电离层建模方法,尤其是涉及一种GNSS实时全球电离层TEC建模方法。
背景技术
电离层时延误差是实时单频单点定位的主要误差之一,其在天顶方向上对GPS L1频率的影响可达29.2m,在卫星高度角较低时甚至超过87.6m,目前用于实时电离层时延修正的模型主要为广播电离层模型,然而该类模型通常为经验模型,例如Klobuchar模型、NeQuick模型、IRI模型等,这些模型的电离层时延修正率一般只有50%到70%。
基于全球分布的几百个GNSS观测站,能够反演得到高精度的全球电离层TEC模型,其中最终电离层产品通常滞后6到9天,快速电离层产品通常滞后1到2天,部分电离层分析中心也发布了全球电离层TEC预报产品,但是这些产品一般只能描述全球电离层TEC的整体变化趋势,不能准确反映电离层TEC的实时变化和局部变化。
以上的广播电离层模型和GNSS事后或预报全球电离层模型并不满足日益增长的高精度的定位需求,因此建立实时的、高精度的、稳定可靠的全球电离层时延修正模型是当前的研究热点,IGS于2002年成立实时工作组,2007年6月开展IGS实时试验计划,2013年初步形成实时数据与产品服务,国际大地测量协会于2016年成立实时电离层监测工作组,目前,IGS各个分析中心正在进行实时全球电离层TEC模型算法研究和精度验证,取得了一些初步进展。
随着IAG/IGS实时试验计划的启动,越来越多的IGS观测站能够提供实时数据服务,基于全球分布的GNSS实时观测站,建立GNSS实时全球电离层TEC模型成为可能。然而,由于目前IGS实时观测站数量不足且分布不均匀,仅仅依靠GNSS实时观测数据难以实现高精度的全球电离层TEC建模,在电离层观测信息缺乏的地区会产生明显的模型误差;目前已知的实时全球电离层TEC模型算法并不成熟,仍然需要进行大量的研究工作。
发明内容
本发明主要是针对当前GNSS实时观测站数量有限且分布不均匀导致仅仅依靠GNSS实时观测数据无法解算出实时的高精度的全球电离层TEC模型而提出,该发明提供了一种基于GNSS观测数据实时建立高精度全球电离层TEC模型的方法。
本发明的上述技术问题主要是通过下述技术方案得以解决:
一种GNSS实时全球电离层TEC建模方法,其特征在于:利用全球分布的GNSS事后观测数据建立事后全球电离层TEC模型,根据事后全球电离层TEC模型数据库建立预报全球电离层TEC模型;利用全球有限分布的GNSS实时数据流提取电离层穿刺点处的电离层观测信息,并结合已经生成的预报全球电离层TEC模型生成当前历元的实时全球电离层TEC模型;整个计算步骤如下:
步骤1,根据GNSS事后双频载波和伪距观测数据提取平滑的电离层观测信息,利用球谐函数建立事后GNSS全球电离层TEC模型,同时解算出球谐函数系数和卫星与接收机端的差分码偏差;步骤2,根据当前时间以前的球谐函数模型系数建立模型系数数据库,并采用最小二乘配置法预报全球电离层TEC模型;步骤3,根据预报球谐函数模型系数时间内插出当前历元的球谐函数模型系数;步骤4,根据GNSS实时数据流提取出电离层观测信息,其中卫星和接收机端的差分码偏差直接采用2天前生成的差分码偏差产品;步骤5,结合预报球谐系数和剔除了差分码偏差的实时观测数据采用最小二乘参数估计法建立实时全球电离层TEC模型。
所述的步骤2更详细的包括以下信息:采用最小二乘配置法预报未来2天的球谐函数模型系数;球谐系数数据库的滑动窗口设置为30天;周期项选取显著的天与亚天周期项;预报的球谐系数包括趋势项和随机项。
所述的步骤3更详细的包括以下信息:直接采用预报球谐系数参与实时全球电离层TEC模型参数解算,避免了根据预报球谐系数计算格网点电离层VTEC引入的计算误差,而且只需要确定预报球谐系数与实时观测数据流的权重分配。
所述的步骤4更详细的包括以下信息:GNSS实时数据流的卫星和接收机端的差分码偏差直接采用2天前估计的差分码偏差产品;对于新增加的GNSS实时观测站,事后差分码偏差产品并不包括相应的卫星和接收机端差分码偏差产品,在本次实时全球电离层TEC建模中剔除相关数据,并将新增观测站加入事后全球电离层TEC建模中,其生成的差分码偏差产品用于下一天实时全球电离层TEC建模。
所述的步骤5更详细的包括以下信息:采用最小二乘方法估计实时全球电离层TEC模型系数,模型的时间分辨率为5分钟;每5分钟存储一次实时观测数据,参与建模的实时观测数据采用15分钟的滑动窗口。
在上述的一种GNSS实时全球电离层TEC建模方法中,主要采用了两种数据源建立实时全球电离层TEC模型,其中预报全球电离层TEC模型作为背景场描述全球电离层TEC的整体趋势,GNSS实时观测数据能够反映各个观测站附近电离层TEC的实时变化和局部变化;结合两种数据源既能够保证全球范围反演电离层TEC的可靠性,又能提高GNSS实时观测数据覆盖地区的模型精度。
因此,本发明具有如下优点:1.采用预报全球电离层TEC模型作为背景场,克服了当前GNSS实时观测站数量有限且分布不均匀导致无法实时解算高精度全球电离层TEC模型的缺点;2.在实时全球电离层TEC建模中,引入了事先估计的卫星和接收机端差分码偏差,在实时数据流预处理阶段剔除了差分码偏差的影响,避免了同时解算球谐函数模型系数与差分码偏差参数估计的相互影响;3.采用加权最小二乘估计实时全球电离层TEC模型系数,能够基于方差分量估计方法最优化确定预报电离层背景场与实时观测数据的权重,确保电离层TEC模型的可靠性和时变性。
附图说明
图1是本发明数据处理流程图;
图2是球谐函数模型系数预报流程图;
图3是实时观测数据电离层TEC提取流程图;
图4是实时全球电离层模型参数估计流程图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例1:
一种GNSS实时全球电离层TEC建模方法,利用全球分布的GNSS事后观测数据建立事后全球电离层TEC模型,根据事后全球电离层TEC模型数据库建立预报全球电离层TEC模型;利用全球有限分布的GNSS实时数据流提取电离层穿刺点处的电离层观测信息,并结合已经生成的预报全球电离层TEC模型生成当前历元的实时全球电离层TEC模型;本发明数据处理流程如图1所示;以下根据算法步骤详细介绍本算法:
步骤1,根据GNSS事后双频载波和伪距观测数据提取平滑的电离层观测信息,利用球谐函数建立事后GNSS全球电离层TEC模型,同时解算出球谐函数系数和卫星与接收机端的差分码偏差;采用球谐函数建立全球电离层TEC模型的基本公式表示为:
步骤3,根据预报球谐函数模型系数时间内插出当前历元的球谐函数模型系数;采用分段线性方法内插当前历元的模型系数,公式表示为:
。
步骤5,结合预报球谐系数和剔除了差分码偏差的实时观测数据采用最小二乘参数估计法建立实时全球电离层TEC模型;采用最小二乘估计法建立实时全球电离层TEC模型的函数模型和随机模型分别表示为;函数模型:式中,表示观测向量,表示预报的球谐系数,表示设计矩阵,表示待估计的球谐系数; 随机模型:式中,表示先验单位权中误差,表示实时观测值的权重,表示虚拟观测值的权重。
实施例2:
在步骤2中,电离层TEC具有典型的周期特性,具体表现为球谐函数模型系数具有显著的周期性。利用球谐函数模型系数的周期性变化规律,基于傅里叶级数可以预报未来多天的球谐函数模型系数,进而生成预报全球电离层TEC模型;本发明中周期项选取显著的天与亚天周期项,采用当前时间前30天的历史数据作为滑动窗口进行外推2天预报;为了刻画电离层TEC的瞬时变化,预报的球谐系数即包括趋势项又包括随机项,因此本发明采用最小二乘配置法进行球谐函数系数参数估计;球谐函数模型系数预报流程如图2所示。
实施例3:
在步骤4中,根据直接获取的GNSS实时观测数据提取的电离层信息即包括了电离层TEC信息又包括了卫星和接收机的差分码偏差信息;差分码偏差是影响实时电离层建模的主要误差之一,且差分码偏差具有短时间稳定特性,因此通常直接采用已有的差分码偏差产品;本发明中获取实时电离层TEC信息时直接采用2天前生成的卫星和接收机差分码偏差产品;对于新增加的GNSS实时观测站,事后差分码偏差产品并不包括相应的卫星和接收机差分码偏差产品,在本次实时全球电离层TEC建模中剔除相关数据,并将新增观测站加入事后全球电离层TEC建模中,其生成的差分码偏差产品用于下一天实时全球电离层TEC建模;实时观测数据电离层TEC提取流程如图3所示。
实施例4:
在步骤5中,本发明中采用最小二乘方法估计实时全球电离层TEC模型系数;为了增加实时电离层模型的时间分辨率,将模型时间分辨率设置为5分钟;因此,每5分钟存储一次剔除了差分码偏差的实时观测数据,并保存为文本格式;由于5分钟观测数据弧段较短,为了增加实时观测数据对实时电离层模型的影响,合并当前历元前三个时段的观测文件,即采用15分钟实时观测数据的滑动窗口参与建模;实时全球电离层模型参数估计流程图如图4所示。
Claims (4)
1.一种GNSS实时全球电离层TEC建模方法,其特征在于:利用全球分布的GNSS事后观测数据建立事后全球电离层TEC模型,根据事后全球电离层TEC模型数据库建立预报全球电离层TEC模型;利用全球有限分布的GNSS实时数据流提取电离层穿刺点处的电离层观测信息,并结合已经生成的预报全球电离层TEC模型生成当前历元的实时全球电离层TEC模型,整个计算步骤如下:
步骤1,根据GNSS事后双频载波和伪距观测数据提取平滑的电离层观测信息,利用球谐函数建立事后GNSS全球电离层TEC模型,同时解算出球谐函数系数和卫星与接收机端的差分码偏差;步骤2,根据当前时间以前的球谐函数模型系数建立模型系数数据库,并采用最小二乘配置法预报全球电离层TEC模型;步骤3,根据预报球谐函数模型系数时间内插出当前历元的球谐函数模型系数;步骤4,根据GNSS实时数据流提取出电离层观测信息,其中卫星和接收机端的差分码偏差直接采用两天前生成的差分码偏差产品;步骤5,结合预报球谐系数和剔除了差分码偏差的实时观测数据采用最小二乘参数估计法建立实时全球电离层TEC模型。
2.根据权利要求1所述的一种GNSS实时全球电离层TEC建模方法,其特征在于,在步骤2中,利用电离层TEC的周期性变化规律,采用最小二乘配置法预报未来2天的球谐函数模型系数;其中球谐系数数据库的滑动窗口设置为30天,周期项选取显著的天与亚天周期项,预报的球谐系数即包括趋势项又包括随机项。
3.根据权利要求1所述的一种GNSS实时全球电离层TEC建模方法,其特征在于,在步骤4中,GNSS实时观测数据的卫星和接收机端的差分码偏差直接采用2天前估计的差分码偏差产品。
4.根据权利要求1所述的一种GNSS实时全球电离层TEC建模方法,其特征在于,在步骤5中,模型的时间分辨率设置为5分钟,其中每5分钟存储一次实时观测数据,参与建模的实时观测数据采用15分钟的滑动窗口 。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911040108.2A CN110764121A (zh) | 2019-10-29 | 2019-10-29 | 一种gnss实时全球电离层tec建模方法 |
PCT/CN2019/123842 WO2021082191A1 (zh) | 2019-10-29 | 2019-12-07 | 一种全球导航卫星系统实时全球电离层电子浓度总含量建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911040108.2A CN110764121A (zh) | 2019-10-29 | 2019-10-29 | 一种gnss实时全球电离层tec建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110764121A true CN110764121A (zh) | 2020-02-07 |
Family
ID=69334754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911040108.2A Pending CN110764121A (zh) | 2019-10-29 | 2019-10-29 | 一种gnss实时全球电离层tec建模方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110764121A (zh) |
WO (1) | WO2021082191A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11686850B2 (en) | 2021-07-15 | 2023-06-27 | Qualcomm Incorporated | Ionosphere grid history and compression for GNSS positioning |
CN113985454A (zh) * | 2021-10-23 | 2022-01-28 | 闽江学院 | 一种顾及方位角的电离层投影函数模型的建模方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130044026A1 (en) * | 2011-02-14 | 2013-02-21 | Trimble Navigation Limited | GNSS Signal Processing with Ionosphere Model for Synthetic Reference Data |
CN103792546A (zh) * | 2012-10-31 | 2014-05-14 | 中国科学院光电研究院 | 一种增量式电离层折射误差修正方法 |
CN104101888A (zh) * | 2014-05-28 | 2014-10-15 | 中国科学院光电研究院 | 一种适合于中国区域的gnss电离层延迟精确建模方法 |
CN105301601A (zh) * | 2015-10-09 | 2016-02-03 | 中国科学院光电研究院 | 一种适用于全球区域的gnss电离层延迟三维建模方法 |
CN106202617A (zh) * | 2016-06-24 | 2016-12-07 | 南京信息工程大学 | 一种gnss、hy‑2和cosmic数据融合建立全球电离层格网vtec模型的方法 |
US20160377730A1 (en) * | 2014-12-16 | 2016-12-29 | Trimble Navigation Limited | Navigation satellite system positioning involving the generation of correction information |
CN106405576A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | 一种基于gnss载波相位观测值的电离层监测与建模方法 |
CN108919316A (zh) * | 2018-08-28 | 2018-11-30 | 中国科学院地质与地球物理研究所 | 一种基于局部球对称假设的单站多系统硬件延迟估计方法 |
CN110058268A (zh) * | 2019-04-30 | 2019-07-26 | 桂林电子科技大学 | 一种基于gnss的电离层闪烁闭环测试系统及方法 |
CN110275185A (zh) * | 2019-07-11 | 2019-09-24 | 武汉大学 | 基于gnss和geo卫星的电离层投影函数建模方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893714B (zh) * | 2010-07-09 | 2013-01-23 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统广播电离层时延修正方法 |
KR101387659B1 (ko) * | 2012-06-01 | 2014-04-21 | 한국표준과학연구원 | 글로벌 위성 항법 시스템을 이용한 전리층 내의 총전자량과 차분코드바이어스 측정시스템, 측정방법 및 그 기록매체 |
CN103197340B (zh) * | 2013-04-01 | 2015-05-13 | 东南大学 | 一种格网化的电离层总电子含量实时监测方法 |
CN110275186B (zh) * | 2019-07-11 | 2020-04-03 | 武汉大学 | Leo卫星增强的gnss电离层归一化与融合建模方法 |
-
2019
- 2019-10-29 CN CN201911040108.2A patent/CN110764121A/zh active Pending
- 2019-12-07 WO PCT/CN2019/123842 patent/WO2021082191A1/zh active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130044026A1 (en) * | 2011-02-14 | 2013-02-21 | Trimble Navigation Limited | GNSS Signal Processing with Ionosphere Model for Synthetic Reference Data |
CN103792546A (zh) * | 2012-10-31 | 2014-05-14 | 中国科学院光电研究院 | 一种增量式电离层折射误差修正方法 |
CN104101888A (zh) * | 2014-05-28 | 2014-10-15 | 中国科学院光电研究院 | 一种适合于中国区域的gnss电离层延迟精确建模方法 |
US20160377730A1 (en) * | 2014-12-16 | 2016-12-29 | Trimble Navigation Limited | Navigation satellite system positioning involving the generation of correction information |
CN105301601A (zh) * | 2015-10-09 | 2016-02-03 | 中国科学院光电研究院 | 一种适用于全球区域的gnss电离层延迟三维建模方法 |
CN106202617A (zh) * | 2016-06-24 | 2016-12-07 | 南京信息工程大学 | 一种gnss、hy‑2和cosmic数据融合建立全球电离层格网vtec模型的方法 |
CN106405576A (zh) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | 一种基于gnss载波相位观测值的电离层监测与建模方法 |
CN108919316A (zh) * | 2018-08-28 | 2018-11-30 | 中国科学院地质与地球物理研究所 | 一种基于局部球对称假设的单站多系统硬件延迟估计方法 |
CN110058268A (zh) * | 2019-04-30 | 2019-07-26 | 桂林电子科技大学 | 一种基于gnss的电离层闪烁闭环测试系统及方法 |
CN110275185A (zh) * | 2019-07-11 | 2019-09-24 | 武汉大学 | 基于gnss和geo卫星的电离层投影函数建模方法 |
Non-Patent Citations (1)
Title |
---|
邓哲民等: "融合IGS和CMONOC的GPS数据反演实测电离层研究", 《中国科学院上海天文台年刊》 * |
Also Published As
Publication number | Publication date |
---|---|
WO2021082191A1 (zh) | 2021-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108761574B (zh) | 基于多源信息融合的降雨量估算方法 | |
CN112070286B (zh) | 复杂地形流域的降水预报预警系统 | |
CN108491682B (zh) | 降雨径流预报系统 | |
CN106468774A (zh) | 一种应用于星基增强系统的星历星钟改正参数及空间信号完好性参数方法 | |
Habarulema et al. | Application of neural networks to South African GPS TEC modelling | |
CN111653088A (zh) | 一种车辆出行量预测模型构建方法及预测方法和系统 | |
CN101893714B (zh) | 全球卫星导航系统广播电离层时延修正方法 | |
Pacione et al. | Combination methods of tropospheric time series | |
Mo et al. | Estimating the extreme wind speed for regions in China using surface wind observations and reanalysis data | |
Tu et al. | A real-time ionospheric model based on GNSS Precise Point Positioning | |
CN108830419B (zh) | 一种基于ecc后处理的梯级水库群入库流量联合预报方法 | |
CN110399634B (zh) | 一种基于天气系统影响的预报区域确定方法及系统 | |
CN110764121A (zh) | 一种gnss实时全球电离层tec建模方法 | |
CN110727923B (zh) | 一种短期潮位预报方法 | |
CN107748372A (zh) | 基于gnss高精度系统的卫星钟差实时估计系统及方法 | |
CN112540390A (zh) | 一种航天器的跟踪预报计算方法及装置 | |
CN108549955B (zh) | 一种充电桩异常率确定方法及装置 | |
Otache et al. | ARMA modelling of Benue River flow dynamics: comparative study of PAR model | |
CN110909447A (zh) | 一种高精度电离层区域短期预报方法 | |
JP2006092058A (ja) | 流量予測装置 | |
Fabbro et al. | GNSS positioning error forecasting in the Arctic: ROTI and Precise Point Positioning error forecasting from solar wind measurements | |
Ge et al. | Reducing satellite orbit error effects in near real‐time GPS zenith tropospheric delay estimation for meteorology | |
CN110648007B (zh) | 一种降雨预报指数校正方法和系统 | |
Carmichael et al. | Model intercomparison study of long range transport and sulfur deposition in East Asia (MICS-Asia) | |
CN108510098B (zh) | 基于卫星遥感的输电线路走廊雪深估算与预警方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20231017 |
|
AD01 | Patent right deemed abandoned |