没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
电离层大气分子是受太阳高能辐射及宇宙射线激发电离产生的高层大气,距地高度约
60~1 000 km,电离层是日地空间的重要组成部分和研究对象,电离层总电子含量(total
electron content,TEC)是描述电离层性质的重要物理参数
[1-2]
。电磁波信号穿过电离层
时,受电离层影响,传播路径和速度发生改变,成为卫星导航中重要的误差源之一
[3-4]
。为
了提高导航定位精度,需要对观测量进行电离层延迟改正,以国际全球导航卫星系统服务
组织(International Global Navigation Satellite System Service,IGS)为代表的全球电离
层云图(global ionosphere maps,GIM)是用于电离层延迟改正的常用方法之一
[5-6]
,GIM
以格网的形式给出了全球范围内电离层总电子含量 TEC 及相应的均方根(root mean
square,RMS)值,适用于全球范围电离层延迟改正。对于区域用户而言,如以 1 000
km 半径作为活动范围的用户,GIM 中区域以外的 TEC 信息没有用到而且也没有必要,相
比于 GIM,利用尽可能少的测站建立区域内电离层 TEC 模型具有更高的效率和更适用的
范围,因此,建立单站区域电离层 TEC 模型,分析其适用范围和精度对区域内用户提高导
航定位精度尤为重要。研究人员分析了 IGS 及欧洲定轨中心(Center for Orbit
Determination in Europe,CODE)、中国科学院、美国喷气推进实验室(Jet Propulsion
Laboratory,JPL)和西班牙加泰罗尼亚理工大学等分析中心的 GIM 精度,以及近年来中
国 GNSS 电离层延迟建模的进展,为电离层 TEC 精度分析提供了一定的参考
[7-8]
;基于球
谐函数和球冠谐函数建立了中国区域电离层 TEC 模型,提供了一定的电离层区域建模参考
[9-10]
;实现了中国电离层 TEC 同化现报系统
[11]
;基于多项式曲面拟合建立并分析了区域电
离层的模型精度
[12]
。以上研究的区域范围较广,所用数据模型也是高阶次模型,而基于单
站低阶次球谐函数建立区域电离层 TEC 模型的相关研究较少。
本文基于 2~15 不同阶次球谐函数分别建立了欧洲区域 16 个 GNSS 跟踪站单站区域
电离层 TEC 模型,并与 CODE、iGMAS 和 IGS 的全球电离层 TEC 格网产品进行了对
比,分析了单站区域电离层 TEC 模型的适用范围和精度。
1. 单站区域电离层 TEC 建模方法
在 GNSS 信号传播过程中,假设沿信号传播路径上的所有电子都集中在无厚度的理
想薄层上,薄层高度取 450 km,并基于薄层假设利用单站全球定位系统(global
positioning system,GPS)双频观测数据,采用载波相位平滑伪距的方法计算传播斜路径
上的倾斜电离层总电子含量(slant total electronic content,STEC)值,根据投影映射函
数将 STEC 转化为测站上方的垂直电离层总电子含量(vertical total electron content,
VTEC),并与球谐函数模型联立,求解出相应的球谐函数系数,最终将球谐函数系数转
化为区域电离层 TEC 格网文件
[13]
。
双频接收机的测距码和载波相位观测方程
[9,14]
为:
Pk,ij=ρij+c×(dti-dtj)+Ik,ij+Tij+ c×rk,i+c×skj+εP,k,ijLk,ij=ρij+c×(dti-dtj)-Ik,ij+Tij+ εp,k,ij-λkNk,ij ]]>
式中,kk 为频率编号;jj 和 ii 分别为卫星编号和接收机编号;PP 和 LL 分别为测距
码和载波相位观测值;ρρ 为接收机和卫星间的几何距离;c 为真空中的光速;dtidti 和
dtjdtj 分别为接收机钟差和卫星钟差;TT 和 II 分别为信号传播路径上的对流层延迟和电离
层延迟;rr 和 ss 分别为接收机延迟偏差和卫星硬件延迟偏差;εε 为测距码的观测噪声;
λλ 为波长;NN 为整周模糊度。分别对 P
1
和 P
2
、L
1
和 L
2
进行差分,可得 P
4
和 L
4
组合:
P4=P1,ij-P2,ij= (I1,ij-I2,ij)+cDj+cDiL4=L1,ij-L2,ij=-(I1,ij-I2,ij)-λ(N1,ij-N2,ij) ]]>
式中,Dj=sj1−sj2Dj=s1j-s2j 和 Di=r1,i−r2,iDi=r1,i-r2,i 分别为卫星码偏差和接收机码
偏差。采用相位观测量 L
4
对伪码观测量 P
4
进行平滑,伪距测量中电离层延迟量为:
I=40.28f2δSTEC ]]>
式中,I 为电离层延迟量;f 为载波相位频率;δSTECδSTEC 为观测路径上倾斜电离层
TEC。将式(3)代入式(2)可得:
P4=40.28(1f12-1f22)δSTEC+cDj+cDiL4=-40.28(1f12-1f22)δSTEC-λ(N1,ij-N2,ij) ]]>
为了提高伪距观测值的精度,采用 Hatch 滤波进行载波相位平滑伪距,对于连续观测
弧段,平滑方式为:
P4,st=k-1k[P4,st-1+L4t- L4t-1]+1kP4t,t>1P4,s1=P41,t=1 ]]>
式中,P4,sP4,s 为平滑后的伪距观测值;t 为历元;k 为从平滑开始的总历元数。
将 P4,sP4,s 代入式(4)中替代 P4P4 可得:
P4,s=40.28(1f12-1f22)δSTEC+cDj+cDi ]]>
式中,f1f1 和 f2f2 分别为载波相位 L
1
和 L
2
的频率。卫星和接收机硬件延迟偏差在一
天内可视为无变化,基于减少待估参数和提高建模精度考虑,本文采用 CODE 解算的卫星
硬件延迟偏差参与建模。
将 STEC 通过单层模型投影函数最终转化为测站上空垂直电离层 VTEC
[15]
,转化过程
为:
F(z)=δSTECδVTEC=11-(ReRe+Hi·sinz)2δVTEC=f12f2240.28(ReRe+Hi·sinz)1-(ReRe+Hi·sinz)2× (ρ2'-ρ1') ]]>
式中,F(z)F(z)为穿刺点处的投影函数;δVTECδVTEC 为测站上空垂直电离层总电子
含量;zz 为接收机天顶距;ReRe 为地球半径;HiHi 为薄层高度。
球谐函数是研究全球范围场力的重要数学模型工具
[16]
,基于电离层薄层假设利用球谐
函数进行单站区域电离层 TEC 建模,即用较少的参数很好地拟合区域电离层 TEC 的整体
情况
[4,17]
。球谐函数为:
δVTEC(φ,λ)=∑n=0nmax∑m=0nP˜nmsinφ×A˜nmcos(mλ)+ B˜nmsin(mλ) ]]>
式中,φφ 与 λλ 分别为穿刺点处的纬度和经度;P˜nmsinφP˜nmsinφ 为正规化后的
n 阶 m 次勒让德函数;A˜mnA˜mn 与 B˜mnB˜mn 分别为球谐函数模型系数。球谐函数阶次
越低,所求系数越少,解算效率越高,但建模范围也有限,即低阶球谐函数在进行大范围
建模时拟合结果较差,系数个数 NpNp 为:
Np=nmax+12-nmax-mmaxnmax-mmax+1 ]]>
将式(7)中第 2 式和式(8)联立,利用最小二乘拟合求解球谐函数模型系数
A˜mnA˜mn 与 B˜mnB˜mn,通过函数转换,将球谐系数转换为通用电离层 TEC 格网格式
[18]
。
2. 电离层 TEC 精度分析方法
CODE 发布的电离层 TEC 最终产品精度为 2~8 TECU,时间滞后约 7 d
[7]
,是利用全
球 300 多个 IGS 跟踪站,基于 15×15 阶球谐函数建立全球电离层 TEC 模型解算得到的;
IGS 发布的电离层 TEC 最终产品精度为 2~8 TECU,时间滞后为 12 d
[19]
,是综合了
CODE 和 JPL 两个分析中心的电离层产品得到的;国内 iGMAS 发布的电离层 TEC 最终产
品精度为 2~8 TECU,时间滞后约 12 d
[20]
,是由 8 个 iGMAS 分析中心的电离层产品加权
得到的。因为这 3 个机构发布的电离层产品比较权威,所以用其作为电离层 TEC 精度比
较分析的基准。
利用单个测站 GNSS 数据建立电离层 TEC 模型,并生成电离层 TEC 格网转换
(ionosphere map exchange format,IONEX)格式文件,与 CODE、IGS 和 iGMAS 的
GIM 格式相同,即将测区按纬度和经度分辨率为 2.5°×5°格网化,给出了每个格网点的电
离层 TEC 值和相应的 RMS 值,每天分为采样间隔为 2 h 的 13 幅电离层 TEC 格网图。取
以观测站为中心,格网经纬度范围分别为 0°×0°、10°×5°、10°×10°、10°×15°、20°×15°
和 20°×20°内所有格网点上电离层 TEC 的平均值和 RMS 的平均值,并分别与 CODE、
iGMAS 和 IGS 的电离层 GIM 产品中相应的 TEC 的平均值和 RMS 的平均值进行对比分
析,格网经纬度范围如图 1 所示。
图 1 格网经纬度范围
Figure 1. Latitude and Longitude Range of Grid
剩余14页未读,继续阅读
资源评论
罗伯特之技术屋
- 粉丝: 4427
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功