CN116736353B - 一种全球-区域-局域多尺度的电离层精细化建模方法 - Google Patents

一种全球-区域-局域多尺度的电离层精细化建模方法 Download PDF

Info

Publication number
CN116736353B
CN116736353B CN202310969531.0A CN202310969531A CN116736353B CN 116736353 B CN116736353 B CN 116736353B CN 202310969531 A CN202310969531 A CN 202310969531A CN 116736353 B CN116736353 B CN 116736353B
Authority
CN
China
Prior art keywords
ionosphere
regional
global
model
satellite
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.)
Active
Application number
CN202310969531.0A
Other languages
English (en)
Other versions
CN116736353A (zh
Inventor
李子申
邱立杰
王宁波
汪亮
王亮亮
王晨旭
孟庆波
侯福荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qilu Aerospace Information Research Institute
Aerospace Information Research Institute of CAS
Original Assignee
Qilu Aerospace Information Research Institute
Aerospace Information Research Institute of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qilu Aerospace Information Research Institute, Aerospace Information Research Institute of CAS filed Critical Qilu Aerospace Information Research Institute
Priority to CN202310969531.0A priority Critical patent/CN116736353B/zh
Publication of CN116736353A publication Critical patent/CN116736353A/zh
Application granted granted Critical
Publication of CN116736353B publication Critical patent/CN116736353B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

本发明公开了一种全球‑区域‑局域多尺度的电离层精细化建模方法,属于卫星定位大气领域,根据北斗监测站骨干网构建全球电离层模型时,结合区域监测站加密网,通过全球电离层改正和精度先验信息,建立区域尺度电离层模型,进一步进行电离层局域分区处理,在保障电离层参数一致性和连续性的前提下,构造单星非差局域电离层改正模型,从而满足用户不同尺度下电离层延迟改正需求。本发明保证在不同范围作业下对多尺度电离层延迟的精细改正,扩展了传统电离层模型的适用范围。

Description

一种全球-区域-局域多尺度的电离层精细化建模方法
技术领域
本发明属于卫星定位大气领域,具体涉及一种全球-区域-局域多尺度的电离层精细化建模方法。
背景技术
电离层是空间大气环境的重要组成部分,其变化不断影响着人类的生产和生活活动。全球卫星导航系统(Global Navigation Satellite System,GNSS)发射的信号穿越电离层时,电离层异常效应将引起卫星通信与导航信号的深度衰落,产生几米到几十米的误差,会严重制约北斗精准定位的载波相位模糊度固定并影响其收敛时间,进而影响终端定位用户服务的连续性和可用性。
电离层精细化监测处理是实现超高精度、超快固定北斗精准定位服务的核心。构建电离层精细化模型可为导航定位用户提供电离层延迟改正,加快导航定位收敛速度,实现对空间天气的精准监测。通过顾及全局-区域-局域等不同尺度空间状态域参数的时变特征与稳定特点,对于进一步提高电离层模型的准确性,具有巨大的现实意义。因此,有必要开展电离层梯度与闪烁效应的精细化建模。
发明内容
为解决上述技术问题,本发明提供一种全球-区域-局域多尺度的电离层精细化建模方法,将全域北斗基准站进行分类,筛选出分布合理且稳定性好的基准站组成骨干网,其目的是用于精确估计全球电离层参数,然后按照定位区域和范围,充分利用区域北斗基站等组成分区加密网,分级分区计算多尺度的高精度大气参数。本发明基于全球“骨干网”及部分区域“加密网”北斗基准站的固定解电离层观测,辅以背景场电离层提供的先验信息,利用球谐函数实现全球电离层改正信息的实时精准建模,综合基准站覆盖区域及无观测区域建模精度情况,精化处理保障全球电离层改正及精度信息的精确一致,实现空间电离层建模由全球范围向区域尺度的不断趋近。
本发明实现其目的所采取的技术方案是:
一种全球-区域-局域多尺度的电离层精细化建模方法,包括如下步骤:
步骤(1)全球电离层模型采用GNSS双频伪距和载波相位观测值计算电离层总电子含量;
步骤(2)仅考虑电离层延迟的一阶项,计算电离层延迟误差;
步骤(3)将倾斜总电子含量STEC转换为垂直方向上的总电子含量VTEC;
步骤(4)选取低阶球谐函数模型,对垂直方向上的总电子含量VTEC进行建模;
步骤(5)引入全球实时电离层改正及精度先验信息,通过虚拟参考星进行约束,结合背景场电离层模型信息,构建电离层残差序列,采用滑动窗口的方法初始化基准估计,进行电离层观测基准统一;
步骤(6)以调整球谐函数为核心,对区域快速电离层进行精化,建立区域尺度下的电离层趋势项模型;
步骤(7)基于步骤(6)中的电离层趋势项模型,构建附加先验电离层约束的区域改正模型并进行迭代处理,进一步准确建立区域快速电离层精度模型;
步骤(8)利用步骤(7)中的区域快速电离层精度模型对区域尺度电离层进行处理,为局域尺度电离层建模提供实时改正及精度信息,结合区域尺度的改正信息,精细处理单星电离层趋势项与随机项,构建残差序列,对残差距离间隔进行划分,计算残差项半方差;
步骤(9)筛选变异函数进行拟合,实现单星精度指标模型的准确构建。
进一步地,上述步骤(1)中,采用GNSS定位的观测方程表达式如下:
式中,P为GNSS伪距观测值;L为GNSS相位观测值;表示卫星,/>表示接收机;/>为卫星和接收机之间的真实几何距离;/>为电离层延迟误差;/>为对流层延迟误差;为接收机钟差;/>为卫星钟差;d为卫星和接收机的差分码偏差;k为频率;b为卫星和接收机的相位超前;N为整周模糊度;ε为随机误差;/>为穿刺点的日固经度,/>为穿刺点的地理经度,/>为太阳的地磁经度。
进一步地,所述步骤(2)中的电离层延迟误差表示为:
式中,f为卫星载波频率;STEC为沿卫星信号路径上的倾斜总电子含量。
进一步地,所述步骤(3)中,垂直方向上的总电子含量VTEC为:
式中,z为卫星的高度角;R为地球的半径;H为电离层薄层高度,即电子密度峰值所在的高度,α为0.9782;为经过相位平滑处理得到的伪距差分观测量,/>和/>为卫星和接收机差分码偏差;c表示光速,/>,/>表示卫星载波频率。
进一步地,所述步骤(4)中的低阶球谐函数模型为:
式中,为球谐函数展开式的最高阶数,/>为n阶m次归一化缔合勒让德函数,Φ为电离层穿刺点的地理纬度或者地磁纬度,/>为穿刺点的日固经度,/>为穿刺点的地理经度,/>为太阳的地磁经度,/>和/>为待估的球谐系数。
进一步地,所述步骤(6)中,利用球冠半角确定覆盖范围,利用球冠极点确定投影变换,利用球冠阶次确定分辨率。
有益效果:
本发明针对电离层建模方法中已有的基准不统一、处理不便捷、精度不自洽等突出问题,提出了一种全球-区域-局域多尺度的电离层精细化建模方法,该方法可以保证在不同范围作业下对多尺度电离层延迟的精细改正,该技术扩展了传统电离层模型的适用范围。
附图说明
图1为本发明的一种全球-区域-局域多多尺度的电离层精细化建模方法流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明在根据北斗监测站骨干网基准站构建全球电离层模型时,结合区域监测站加密网基准站,通过全球电离层改正和精度先验信息,建立区域尺度的电离层模型,进一步进行电离层局域分区处理,构造单星非差局域电离层改正模型,从而满足用户多尺度的电离层延迟改正需求。
如图1所示,本发明的一种全球-区域-局域多尺度的电离层精细化建模方法具体包括如下步骤:
步骤(1)全球电离层模型采用GNSS双频伪距和载波相位观测值计算电离层总电子含量(total electron content,TEC),其采用GNSS定位的观测方程表达式如下:
式中,P为GNSS伪距观测值;L为GNSS相位观测值;表示卫星,/>表示接收机;/>为卫星和接收机之间的真实几何距离;/>为电离层延迟误差;/>为对流层延迟误差;/>为接收机钟差;/>为卫星钟差;d为卫星和接收机的差分码偏差;k为频率;b为卫星和接收机的相位超前;N为整周模糊度;ε为随机误差,c表示光速;/>为穿刺点的日固经度,/>为穿刺点的地理经度,/>为太阳的地磁经度。上述参数包括星历数据和站坐标。
步骤(2)计算一阶项电离层延迟误差:
仅考虑电离层延迟的一阶项,电离层延迟误差可以表示为:
式中,f为卫星载波频率;STEC为沿卫星信号路径上的倾斜总电子含量。考虑到伪距观测值的噪声较大,通常采用载波相位平滑伪距的方法提高实测STEC的估计精度。
步骤(3)将倾斜总电子含量STEC转换为垂直方向上的总电子含量VTEC:
式中,z为卫星的高度角;R为地球的半径;H为电离层薄层高度,即电子密度峰值所在的高度,一般选取350~450km,α为0.9782;为经过相位平滑处理得到的伪距差分观测量,/>和/>为卫星和接收机差分码偏差,一天之中一般可认为其是一个常数;c表示光速,/>,/>表示卫星载波频率。
步骤(4)对垂直方向上的总电子含量VTEC进行建模,选取低阶球谐函数模型:
式中,为球谐函数展开式的最高阶数,/>为n阶m次归一化缔合勒让德函数,Φ为电离层穿刺点的地理纬度或者地磁纬度,/>为穿刺点的日固经度,/>为穿刺点的地理经度,/>为太阳的地磁经度,/>和/>为待估的球谐系数。
所述步骤(4)实现全球电离层谐函数精确建模。
步骤(5)引入全球实时电离层改正及精度先验信息(即先验精度驱动),通过虚拟参考星进行约束,辅以背景场电离层模型信息(即背景场信息约束),构建电离层残差序列,采用滑动窗口的方法初始化基准估计,进行电离层观测基准统一。
步骤(6)采用调整球谐函数为核心,对区域快速电离层进行精化,建立区域尺度的电离层趋势项模型,其中球冠半角确定覆盖范围,球冠极点确定投影变换,球冠阶次确定分辨率,即实现区域电离层观测基准统一与精化建模。
步骤(7)基于步骤(6)中的趋势项模型,构建附加先验电离层约束的区域改正模型并进行迭代处理,进一步准确建立区域快速电离层精度模型(即建立区域尺度电离层模型),满足大站间距条件下区域电离层的精准改正需求。
步骤(8)充分利用步骤(7)中精准改正的区域快速电离层模型对区域尺度的电离层进行处理,为局域尺度电离层建模提供实时改正及精度信息,结合区域尺度的改正信息,精细处理单星电离层趋势项与随机项(即局域单行电离层趋势项与随机项精细处理),构建残差序列,对残差距离间隔进行划分,计算残差项半方差。
步骤(9)筛选变异函数进行拟合,实现单星精度指标模型的准确构建,电离层模型精度由区域范围向局部尺度进一步趋近。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种全球-区域-局域多尺度的电离层精细化建模方法,其特征在于,包括如下步骤:
步骤(1)全球电离层模型采用全球卫星导航系统GNSS双频伪距和载波相位观测值计算电离层总电子含量;
步骤(2)仅考虑电离层延迟的一阶项,计算电离层延迟误差;
步骤(3)将沿卫星信号路径上的倾斜总电子含量STEC转换为垂直方向上的总电子含量VTEC;
步骤(4)选取低阶球谐函数模型,对垂直方向上的总电子含量VTEC进行建模;
步骤(5)引入全球实时电离层改正及精度先验信息,通过虚拟参考星进行约束,结合背景场电离层模型信息,构建电离层残差序列,采用滑动窗口的方法初始化基准估计,统一电离层观测基准;
步骤(6)调整球谐函数,对区域快速电离层进行精化,建立区域尺度下的电离层趋势项模型;
步骤(7)基于步骤(6)中的电离层趋势项模型,构建附加先验电离层约束的区域改正模型并进行迭代处理,进一步准确建立区域快速电离层精度模型;
步骤(8)利用步骤(7)中的区域快速电离层精度模型对区域尺度电离层进行处理,为区域尺度电离层的建模提供实时改正及精度信息,结合区域尺度的改正信息,精细处理单星的电离层趋势项与随机项,构建残差序列,对残差距离间隔进行划分,计算残差项半方差;
步骤(9)筛选变异函数进行拟合,实现单星精度指标模型的准确构建。
2.根据权利要求1所述的一种全球-区域-局域多尺度的电离层精细化建模方法,其特征在于,上述步骤(1)中,采用全球卫星导航系统GNSS定位的观测方程表达式如下:
式中,P为GNSS伪距观测值;L为GNSS相位观测值;上标表示卫星,下标/>表示接收机;为卫星和接收机之间的真实几何距离;/>为电离层延迟误差;/>为对流层延迟误差;/>为接收机钟差;/>为卫星钟差;/>为卫星的差分码偏差,/>为接收机的差分码偏差;k为频率;b为卫星和接收机的相位超前;N为整周模糊度;/>为伪距随机误差,/>为相位随机误差;/>为穿刺点的日固经度,/>为穿刺点的地理经度,/>为太阳的地磁经度。
3.根据权利要求2所述的一种全球-区域-局域多尺度的电离层精细化建模方法,其特征在于,所述步骤(2)中的电离层延迟误差表示为:
式中,f为卫星载波频率;STEC为沿卫星信号路径上的倾斜总电子含量。
4.根据权利要求3所述的一种全球-区域-局域多尺度的电离层精细化建模方法,其特征在于,所述步骤(3)中,垂直方向上的总电子含量VTEC为:
式中,z为卫星的高度角;R为地球的半径;H为电离层薄层高度,即电子密度峰值所在的高度,α为0.9782;为经过相位平滑处理得到的伪距差分观测量,/>和/>为卫星差分码偏差和接收机差分码偏差;c表示光速,/>,/>表示卫星载波频率。
5.根据权利要求4所述的一种全球-区域-局域多尺度的电离层精细化建模方法,其特征在于,所述步骤(4)中的低阶球谐函数模型为:
式中,为球谐函数展开式的最高阶数,/>为n阶m次归一化缔合勒让德函数,Φ为电离层穿刺点的地理纬度或者地磁纬度,/>为穿刺点的日固经度,/>为穿刺点的地理经度,/>为太阳的地磁经度,/>和/>为待估的球谐系数。
6.根据权利要求5所述的一种全球-区域-局域多尺度的电离层精细化建模方法,其特征在于,所述步骤(6)中,利用球冠半角确定覆盖范围,利用球冠极点确定投影变换,利用球冠阶次确定分辨率。
CN202310969531.0A 2023-08-03 2023-08-03 一种全球-区域-局域多尺度的电离层精细化建模方法 Active CN116736353B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310969531.0A CN116736353B (zh) 2023-08-03 2023-08-03 一种全球-区域-局域多尺度的电离层精细化建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310969531.0A CN116736353B (zh) 2023-08-03 2023-08-03 一种全球-区域-局域多尺度的电离层精细化建模方法

Publications (2)

Publication Number Publication Date
CN116736353A CN116736353A (zh) 2023-09-12
CN116736353B true CN116736353B (zh) 2023-11-07

Family

ID=87917100

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310969531.0A Active CN116736353B (zh) 2023-08-03 2023-08-03 一种全球-区域-局域多尺度的电离层精细化建模方法

Country Status (1)

Country Link
CN (1) CN116736353B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104101888A (zh) * 2014-05-28 2014-10-15 中国科学院光电研究院 一种适合于中国区域的gnss电离层延迟精确建模方法
CN110568459A (zh) * 2019-08-28 2019-12-13 桂林电子科技大学 基于igs和cors站的区域电离层tec实时监测方法
CN111123300A (zh) * 2020-01-13 2020-05-08 武汉大学 近实时大范围高精度电离层电子密度三维监测方法及装置
US10852439B1 (en) * 2020-04-30 2020-12-01 Beihang University Global ionospheric total electron content prediction system
CN112034489A (zh) * 2020-07-20 2020-12-04 中国科学院空天信息创新研究院 一种基于多源数据融合的全球电离层格网生成方法
CN114545458A (zh) * 2021-10-23 2022-05-27 闽江学院 一种高精度电离层实时建模方法
CN115857058A (zh) * 2022-11-25 2023-03-28 国网思极位置服务有限公司 一种电离层数据分析模型建设方法及其终端
CN115982564A (zh) * 2022-12-05 2023-04-18 广东电网有限责任公司 一种区域电离层电子密度计算方法、装置及计算机设备

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104101888A (zh) * 2014-05-28 2014-10-15 中国科学院光电研究院 一种适合于中国区域的gnss电离层延迟精确建模方法
CN110568459A (zh) * 2019-08-28 2019-12-13 桂林电子科技大学 基于igs和cors站的区域电离层tec实时监测方法
CN111123300A (zh) * 2020-01-13 2020-05-08 武汉大学 近实时大范围高精度电离层电子密度三维监测方法及装置
US10852439B1 (en) * 2020-04-30 2020-12-01 Beihang University Global ionospheric total electron content prediction system
CN112034489A (zh) * 2020-07-20 2020-12-04 中国科学院空天信息创新研究院 一种基于多源数据融合的全球电离层格网生成方法
CN114545458A (zh) * 2021-10-23 2022-05-27 闽江学院 一种高精度电离层实时建模方法
CN115857058A (zh) * 2022-11-25 2023-03-28 国网思极位置服务有限公司 一种电离层数据分析模型建设方法及其终端
CN115982564A (zh) * 2022-12-05 2023-04-18 广东电网有限责任公司 一种区域电离层电子密度计算方法、装置及计算机设备

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
An end-to-end approach for near real time ionosphere monitoring over mid-latitudes from GPS data using kriging interpolation and IGS products;Elhadi Takka;2018 IEEE/ION Position, Location and Navigation Symposium (PLANS);全文 *
GNSS区域电离层建模与精度分析;袁浩鸣;中国硕士学位论文全文数据库 基础科学辑;全文 *
利用球冠谐分析方法和GPS数据建立中国区域电离层TEC模型;柳景斌;武汉大学学报· 信息科学版;全文 *
基于拟合推估的中国区域电离层延迟精确建模方法;袁运斌;李子申;王宁波;张宝成;霍星亮;;导航定位学报(03);全文 *
多模GNSS实时电离层精化建模及其应用研究;张瑞;中国博士学位论文全文数据库 基础科学辑;全文 *
高精度的连续电离层延迟建模研究;范国清;王威;郗晓宁;魏立栋;;空间科学学报(02);全文 *

Also Published As

Publication number Publication date
CN116736353A (zh) 2023-09-12

Similar Documents

Publication Publication Date Title
US6756938B2 (en) Satellite positioning system receivers and methods therefor
CN109709591B (zh) 一种面向智能终端的gnss高精度定位方法
EP2356403B1 (en) Using magnetometer with a positioning system
CN110275185B (zh) 基于gnss和geo卫星的电离层投影函数建模方法
US9557422B1 (en) Systems, methods, devices and subassemblies for creating and delivering a GNSS augmentation service
CN1756965B (zh) 用于无线辅助混合定位系统中位置确定的方法和装置
US20110063167A1 (en) Using magnetometer with a positioning system
WO2010000156A1 (zh) 运用广播调频信号定位地理位置的方法及系统
CN110441795B (zh) 一种基于时空结构信息的中国区域电离层vtec精确建模方法
CN104101888A (zh) 一种适合于中国区域的gnss电离层延迟精确建模方法
CN103018758A (zh) 基于gps/ins/agps的移动差分基站方法
US9612340B1 (en) Systems, methods, devices and subassemblies for creating and delivering crowd-sourced GNSS models
Ren et al. Global ionospheric modeling using multi-GNSS and upcoming LEO constellations: Two methods and comparison
US20170026797A1 (en) Systems and methods for using doppler measurements to estimate a position of a receiver
US20080249713A1 (en) Gps position accuracy using feedback from a map database
Ma et al. Indoor and outdoor positioning system based on navigation signal simulator and pseudolites
CN109633718B (zh) 一种归一化加权最小二乘导航定位方法
Ren et al. Multi-GNSS contributions to differential code biases determination and regional ionospheric modeling in China
CN113985454A (zh) 一种顾及方位角的电离层投影函数模型的建模方法
Rose et al. The use of ionospheric tomography and elevation masks to reduce the overall error in single-frequency GPS timing applications
CN110146904B (zh) 一种适用于区域电离层tec的精确建模方法
Guo et al. Instantaneous velocity determination and positioning using Doppler shift from a LEO constellation
Pinell et al. Receiver architectures for positioning with low earth orbit satellite signals: a survey
Bahadur et al. Real-time single-frequency multi-GNSS positioning with ultra-rapid products
CN111913199B (zh) 基于变分模态分解的手机gnss数据噪声提取方法

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
GR01 Patent grant
GR01 Patent grant