《实时区域电离层TEC建模、预报及差分码偏差估计》由会员分享，可在线阅读，更多相关《实时区域电离层TEC建模、预报及差分码偏差估计（7页珍藏版）》请在金锄头文库上搜索。
1、实时区域电离层TEC建模、预报及差分码偏差估计畅鑫，张伟武汉大学测绘学院，武汉 430079摘 要：电离层总电子含量（TEC）模型对于导航，精密定位以及其他相关应用有重要意义，能否有 效地消除或减弱电离层延迟误差关系到众多单频GNSS接收机用户导航与定位的精度与可靠性。目前中 国连续地面参考运行（CORS）系统的高速发展给实时精确建立区域电离层模型提供了条件。本文将使 用电离层残差组合观测值和低阶球谐函数模型对区域电离层TEC建模，同时估计差分码偏差（DCBs） 和VTEC。广域定位中，由于区域跨度大，观测站分布较稀疏，平均站间距较大，故选择欧洲均匀分布 的14个IGS观测站将组成一个大型的C
2、ORS网，VTEC模型系数15分钟结算一次，差分码偏差一天结算一 组结果。在与IGS分析中心CODE发布模型的对比中得出，差分码偏差的差值的平均值小于0.35 ns，RMS 小于0.2 ns，VTEC差值基本小于2TECU，作为预报的VTEC模型精度95%在1TECU内，在单频单点伪距静 态定位中，较之CODE模型也有较大改善。关键词：CORS ；电离层；区域模型；预报；总电子含量；差分码硬件偏差1 引言电离层总电子含量（TEC）及其变化不但是电离 层形态学研究的重要资料，也是精密定位、导航和电 波科学中电离层改正的重要参数，在美国取消SA政策 后，电离层延迟成为了影响定位和导航的最大误差源。
3、 在精密定位中，电离层的准确估计将更好的改正 GNSS观测值，同时高精度的电离层估计对空间大气、 地球观测等方面都有重要意义i，2。IGS于1998年采用 Schaer3 等提出的电离层总电子含量数据交换格式文 件IO NEX，同年成立IGS电离层工作组发布了全球电 离层图（GIM）,提供卫星和接收机频率间码延迟偏差 DCB信息。Gao YU等对二维单层模型和三维层析模 型进行了对比分析。萧佐4对电离层模型进行了系统 的分类，将电离层模型分为统计、经验及物理等几种。 Schaerb结合CODE分析中心的全球电离层模型对利 用GPS技术探测电离层理论进行了详细介绍。GPS电 离层探测技术可以反演
4、电离层变化，对电离层物理特 性及其观测进行研究。张小红、李征航等 6人对利用 双频GPS观测数据建立电离层延迟模型进行了深入研 究。袁运斌，欧吉坤利用GPS研究了电离层延迟及 电子浓度变化的规律。章红平着重对利用地基GPS 进行电离层模型建立、数据分析处理，电离层时空变 化的监测进行了研究。GPS差分码偏差（Difference Code Bias，DCB） 包括卫星端和接收机端差分码偏差，是指同时刻同频 率或不同频率不同伪距码观测量之间的时间偏差。 DCB是一个相对量，根据接收机的不同，可分为P1码 /P2码、P1码/C1码及C1码/P2码等。差分码偏差直接影 响C1码、P1码及P2码相对测
5、量精度，该参数的精密 确定对提高GPS精密单点定位精度、GPS时间同步精 度以及GPS电离层监测精度等具有重要作用。近年来，随着我国CORS的迅猛发展，如何利用 区域CORS数据高精度实时电离层模型，尤其是针对 广域定位，观测站分布稀疏，站间距较大的情况下， 准确消去电离层影响对提高单频接收机定位精度具有 重要应用价值。本文将研究基于P4 （电离层残差）组 合，将卫星端和接收机端DCB作为参数参于球谐函数 区域电离层建模的方法，准确估计区域电离层模型及 DCBs，并探讨使用较短时间间隔确定的VTEC模型作 为后一时段的预报模型的精度与可靠性，采用单频单 点伪距静态定位检验其效果。2 GPS电离
6、层探测方法电离层是由电离化的等离子体等组成的距地球表 面50-2000km的大气层区域。根据电子密度可以分为 高度不同的层，一般在350km电子密度达到峰值。对 于二维电离层模型，一般采用薄壳模型，如图1将整个 电离层压缩成一个高度为H没有厚度的薄壳10，以总 电子含量TEC描述其性质，TEC是底面积为1m2的贯穿 整个电离层的柱体中的自由电子数，通常用TECU（1 TECU=1016Ne/m2）表示。Figure 1. Ionospheric Shell双频接收机具有L (1575.42 MHz )和L(1227.6012MHz )两个频率载波相位观测值及其加载的伪距码观Pi= P i +
7、d i + d i + c(T i T ) + d i + d + E ik, j0, jion,k,jtrop, jj k k, jP,k,j(1) L = pi + di + di + c(T i T ) X(bi+ Ni)k, j 0, jion,k,jtrop, jjk, jk,j+ E iL,k,ji为第i颗卫星；j为第j号接收机。P，L，L进行差分，得到P和L 2124组合：P = Pi Pi = (di di ) + DCBi + DCB41, j 2, jion ,1, jion ,2, jjL = Li Li41, j 2, j=(di di)入(bi bi )入(Ni Ni
8、ion,1, jion,2, j1,j2, j1,j(3)2,j)(5)(9)(10)sTEC = 宀 (P cDCBi cDCB )40.28(f2 f2) 4,s12Rmf (z) = cos(arcsin(sin(a z)R+HF = 121,240.28( f 2 f 2)12vTEC = mf (z)QF 乂P cDCBi cDCB )1,24, sj(11)(12)(13)P =A (n, m) P =nmnm测值，其观测方程为：(2)其中：P为GPS伪距观测值；L为GPS载波相位观测 值； P 为接收机与卫星真实几何距离； d 为电离层ion延迟量；d 为对流层延迟量；c为光速；
9、t i为卫星trop钟差； T 为接收机钟差； d 为卫星、接收机码硬件偏 j差； b 为卫星、接收机载波相位硬件偏差； N 为载波 相位模糊度；E为GPS观测值残差；k = 1,2表示L的 两个频率对P，1 离层残差)(4) 其中 DCBi =didi ， DCB =d d 分别为卫星和12j 1,j 2,j接收机差分码偏差，实际为P相对P的偏差量。由于 12伪距精度比较低，故采用L平滑P。平滑之前应使用 44MW ( Melbourne-Wubbena combination )和电离层残差法对 L 和 L 进行周跳探测。下面给出平滑方程11： 12()(、丄( )4, s ) + )4
10、t ( ) + ( )1 t2 t1 t2 t(P ) +5(L )4, s t 14 t , t 11 t a 2(P4)t(W ) = 12 t a2+a(P4)t5 (L ) = (L ) (L )(8)4 t ,t14 t4 t1其中t为历元数，P为载波相位平滑伪距，O为平滑4, s方程的权。由载波相位平滑伪距方程(6)可以看出P4, s 中仍包含待定的DCBi和DCB，它们将与电离层延迟 j参数一起被估计。在忽略二、三阶项的情况下，电离层对GPS伪距观测值所造成的影响可以表示为12：,40.28d = sTECionf 2其中 f 为载波频率， sTEC 为倾斜路径总电子含量将式(9
11、)代入式( 3)，并进行相位平滑伪距得至P = 40.28(丄-丄)sTEC + DCBi + DCB4,sf12 f22j经过相位平滑提取出的DCBs将更为可靠。3实时球谐函数模型建立及DCBs估计 方法易由式(10)得到STEC,进而得至到VTEC：(14)其中z为卫星高度角，R为地球半径,H为电离层薄 壳高度，为准确与CODE结果比较，本文对H，R和 a 均采用CODE的设定，H =506.7km, R =6378km, a =0.9782，该投影函数与JPL的扩展单层模型(Extended Slab Model，ESM)的投影符合的最好13。 球谐函数模型是IGS分析中心之一 CODE
12、生成全 球电离层产品所采用的模型(15阶)。对于区域范围 建模选择低阶球谐函数模型效果更好，以下给出球谐 函数模型，E(P,s)即为VTEC14：E(卩,s) = 工 P (sin 卩)(a cos ms + b sin ms) (15)nmnmnmn = 0 m=02 2n +1 (n - m)!1 + 5(n + m)!0 m、式mf (z )QF =1,2 P cDCBi cDCB4, sjRf 2 f 2cos(arcsin(sin(a z)(itR + H(17)40.28( f 2 f 2)12呢乞 P (sin 0 )(a cos ms + b sin ms)nmnmnmn0 m
13、=0(18)图2.欧洲IGS测站分布Figure 2. Distributions of IGS stations chosen from Europe4.2 DCBs 分析PRN: #l-#32100 101102103104105106107108109DOY4 2 0 2 4 6 8O.C5O.O.O.C5(su)ssaQuID居七J图3卫星DCB与CODE差值Figure 3. Differences in the satellites DCB value from this paperand CODE8 6 4 2 0CJO.C5 6ONSA m POTS S PTBBWTZA其中0
14、是穿刺点的地刺纬度，s =九-九为穿刺点处在 0太阳-地磁参考框架下的地方时角，P为正则化的勒让德多项式，a ，b为球谐函数系数，max_n为球谐nmnm函数最大阶数。5为克罗内克函数。由式(11)(12)和式(15 )可得到：E (0, s)P - cDCBi cDCB4, sj其中 a ,b ,DCBi, DCB 为待估计的未知量，在进行 nm nm j估计时，以2h为一个时段，在一个时段内认为其卫星 及接收机DCB是不变的，一个时段得到一组4阶球谐 函数模型的25个系数，及所有参与解算的GPS卫星和 测站接收机的DCB值。根据最小二乘原理，为了分离 出卫星DCB值还需要加入式(18)中的
15、限制条件：max_ i艺 DCBi 0(19)i =1平差随机模型中权阵为单位阵，即观测值相互独立。球谐电离层模型系数的更新率通常取决于地面基 准站网的密度及分布，一般5-15分钟更新一次13， 但 DCB值变化缓慢，故本文采取15分钟结算一次V TEC 模型，并探讨使用前15分钟结算的V TEC模型作为后 15分钟V ETC预报模型的精度与可靠性。使用累计数 据一天结算一组DCB值。模拟实时过程中，逐历元读 取事后数据，达到15分钟数据时，进行建模估计并发 布，采用观测网中某站作为检测站根据模型结果计算 出V TEC及DCBs，并于CODE结果作对比。4 实验及结果分析4.1 实验数据选取了欧洲德国及其周边地区分布均匀的IG
