GNSS电离层视线方向电子含量(STEC)提取

目录

一、基于GNSS观测值提取STEC

二、伪距法提取

三、相位平滑伪距法提取


一、基于GNSS观测值提取STEC

        GNSS信号主要包括测码伪距和载波相位两种常用观测量,其观测方程可以表示为:

P^{_{r,j}^{s}} = \varrho{_{r,j}}^{s}+c*(dt_{r}-dt^{s})+d_{orb}^{s}+T_{r}^{s}+I_{r,j}^{s}+c*B_{r,j}+c*B_{j}^{s}+rel^{_{r,j}^{s}}+\xi

\phi ^{_{r,j}^{s}} = \varrho{_{r,j}}^{s}+c*(dt_{r}-dt^{s})+d_{orb}^{s}+T_{r}^{s}-I_{r,j}^{s}-c*b_{r,j}-c*b_{j}^{s}+rel^{_{r,j}^{s}}+\lambda^{_{j}}N^{_{r,j}^{s}}+\epsilon​​​​​​​​​​​​​​

        其中:

        P^{_{r,j}^{s}}\phi ^{_{r,j}^{s}}分别表示接收机r观测卫星s在第j个频率上的伪距和相位观测值

       \varrho{_{r,j}}^{s}表示接收机r到卫星s的几何距离

        dt_{r}dt^{s}表示接收机r、卫星s的钟差

        d_{orb}^{s}表示卫星s的轨道误差

        T_{r}^{s}表示接收机r到卫星s观测路径上的对流层延迟

        I_{r,j}^{s}表示接收机r到卫星s观测路径上在第i个频率上的电离层延迟

        B_{r,j}B_{j}^{s}表示接收机端与卫星端在第i个频率上的伪距硬件延迟

        b_{r,j}b_{j}^{s}表示接收机端与卫星端在第i个频率上的相位小数偏差

        rel^{_{r,j}^{s}}表示接收机r与卫星s的相对论效应

        从上式可以看到除了电离层延迟硬件延迟整周模糊度频率有关外,其他误差项均与频率无关。因此可以通过两个或两个以上的无几何距离组合消除频率无关项的误差,进而方便求解电离层观测值

P^{_{r,4}^{s}} =P^{_{r,1}^{s}}-P^{_{r,2}^{s}} =I_{r,1}^{s}-I_{r,2}^{s}+c*B_{r,1}-c*B_{r,2}+c*B_{1}^{s}-c*B_{2}^{s}

\phi ^{_{r,4}^{s}} =\phi ^{_{r,1}^{s}} -\phi ^{_{r,2}^{s}} =-(I_{r,1}^{s}-I_{r,2}^{s})-c*(b_{r,1}-b_{r,2})-c*(b_{1}^{s}-b_{2}^{s})+\lambda _{1}N_{1}-\lambda _{2}N_{2}

        令;

DCB_{_{r,12}}=B_{r,1}-B_{r,2}

DCB_{12}^{s}=B_{1}^{s}-B_{2}^{s}

DPB_{_{r,12}}=b_{r,1}-b_{r,2}

DPB_{12}^{s}=b_{1}^{s}-b_{2}^{s}

        则电离层观测值可简化为:

        P^{_{r,4}^{s}} =P^{_{r,1}^{s}}-P^{_{r,2}^{s}} =I_{r,1}^{s}-I_{r,2}^{s}+c*DCB_{_{r,12}}+c*DCB_{12}^{s}

\phi ^{_{r,4}^{s}} =\phi ^{_{r,1}^{s}} -\phi ^{_{r,2}^{s}} =-(I_{r,1}^{s}-I_{r,2}^{s})-c*DPB_{_{r,12}}-c*DPB_{12}^{s}+\lambda _{1}N_{1}-\lambda _{2}N_{2}

        目前利用上式的电离层观测值进行TEC提取的方法有三种:

  1. 伪距法,即直接利用无几何距离组合获取电离层TEC观测量,方法简单、易于计算,单精度较低。
  2. 相位平滑伪距法,即无几何距离观测值和伪距观测值计算电离层TEC观测值,方法简单,精度较高。
  3. 非组合PPP法,利用非组合PPP模型将电离层延迟作为待估参数估计,相较于相位平滑伪距法,其提取方法复杂,但是精度要比相位平滑伪距法高。

二、伪距法提取

        伪距法提取的关键方程为:

P^{_{r,4}^{s}} =P^{_{r,1}^{s}}-P^{_{r,2}^{s}} =I_{r,1}^{s}-I_{r,2}^{s}+c*DCB_{_{r,12}}+c*DCB_{12}^{s}

        在正确改正了卫星和测站上的DCB后得到STEC观测量,其中星站上的DCB可以通过上海天文台CAS发布的BSX文件获取。

三、相位平滑伪距法提取

        直接将无几何距离伪距值P^{_{r,4}^{s}}和相位值\phi ^{_{r,4}^{s}}求和,即可得到相位模糊度Na,其中Na表示整周模糊度和相位小数偏差的合并项(同样和上面伪距法一样需要估计或者读取相关产品获取DCB):

Na=P^{_{r,4}^{s}}+\phi ^{_{r,4}^{s}}+c*DCB_{_{r,12}}+c*DCB_{12}^{s}

        上式中受伪距观测值噪声的影响,单历元求解的模糊度误差较大。鉴于模糊度在一个连续观测弧段内,如果不发生周跳,其值不会改变。因此可以通过对一个连续弧段取平均求模糊度的方式来削弱伪距观测值噪声的影响。若连续弧段长度为n个历元,k表示历元号则有:

\bar{Na}=\frac{1}{n}\sum (P^{_{r,4}^{s}}+\phi ^{_{r,4}^{s}}+c*DCB_{_{r,12}}+c*DCB_{12}^{s})

        由此得到STEC的提取方程:

\phi ^{_{r,4}^{s}}-Na =-(I_{r,1}^{s}-I_{r,2}^{s})

四、实验结果和分析总结

图1JFNG站上STEC提取结果对比

图2 STEC转换为VTEC结果对比

        由图中可以看出,虽然伪距法得到不同卫星的电离层TEC 观测值的站间单差值貌似完美重合,但只是均值意义上的貌似重合,其波动范围非常大。而相位平滑伪距法得到的电离层TEC观测值每颗卫星弧段内非常平滑,而且波动幅度很小,说明该方法很好的平滑了伪距噪声影响。

评论 20
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值