-
在GPS地面沉降监测中,为降低监测成本,可采用单双频混合模式,即大部分监测点使用单频接收机,利用少部分双频接收机数据建立区域电离层改正模型,从而提高单频机的精度[1-2]。欲建立电离层改正模型,须对电离层的影响规律有清楚的认识。电离层活动虽然瞬息万变,但具有明显的年、季、日周期变化规律及一定的空间分布特征[3-4]。目前对电离层变化规律的研究集中于全球或区域电离层总电子含量(total electric content, TEC)的变化[4-5],对高精度相对定位,特别是较短基线应用测量精度影响的相关研究还较少。电离层延迟改正方法主要有3类:(1)经验模型,如KlobHchar模型、Bent模型;(2)区域拟合模型[6-8];(3)参数估计改正方法[9-10]。在单双频混合GPS地面沉降监测中宜采用第2类方法,即利用双频基线提取双差电离层延迟, 通过一定的内插模型内插得到单频基线的电离层延迟[11]。
本文采用广州南沙地面沉降监测网2013年的单天静态解数据,首先分析了电离层对GPS相对定位的影响,以及该影响与基线长度的关系和在时间域的变化规律;然后从观测值域建立电离层延迟改正模型,对单频地面沉降监测点电离层延迟的改正效果进行了分析。
-
单频和双频消电离层组合定位解算方程为[12]:
$$ {L_i} = \rho - {I_i} + {N_i}\;\;\;\;\left( {i = 1,2} \right) $$ (1) $$ {L_c} = \frac{{f^2_1{L_1} - f^2_2{L_2}}}{{f^2_1 - f^2_2}} = \rho + {N_c} $$ (2) 其中,f1和f2分别是GPS信号的2个频率;Li(i=1, 2)为对应频率的相位观测值(m); Lc为L1和L2消电离层一阶项的线性组合观测值(m);ρ是站星距(m);Ni表示Li的相位模糊度(m);Nc为Ni线性组合项,其值为(f12N1-f22N2)/(f12-f22);Ii为对应频率的电离层延迟误差项(m)。
电离层延迟对GPS定位结果影响如图 1所示,S(i)表示接收机可观测到的卫星(i=1, 2…n, n为卫星数),R和R′分别表示单频和双频消电离层组合解算得到的流动站位置,记电离层延迟引起的流动站的点位估值变化为RR′。
-
采用中国广州南沙地面沉降监测网的数据,测区如图 2所示,点位分布如图 3所示。图 3中三角形表示双频监测点,圆点表示单频监测点。数据时间跨度为1 a(其中GD11点2013-01-01-2013-07-12期间的数据缺失),采样间隔20 s。采用事后处理软件HGO进行数据解算。
图 2 南沙GPS地面沉降监测网区域概况图
Figure 2. Regional Overview Diagram of Nansha GPS Ground Deformation Monitoring Network
图 3 南沙GPS地面沉降监测网点位分布图
Figure 3. Point Distribution Plot of Nansha GPS Ground Deformation Monitoring Network
为了对比单双频数据解算质量,选用图 3中所有的双频点构成独立基线,分别进行单频基线定位解算和双频消电离层定位解算,其中单频解算时只采用L1频率。单、双频解算结果减去各自平均值绘于图 4,图 4中蓝色表示单频结果,红色表示双频结果,横轴表示时间,纵轴表示偏差量,单位为m。
从图 4中可以看出,不论是东(east, E)方向还是北(north, N)和高(height,H)方向,双频结果都要优于单频结果,随着基线长度的增加,单频结果残差逐渐增大,而双频结果基本不变。1 a的静态解算结果中单频数据存在部分天数解算不合格的现象,而且随着基线长度的增加,不合格率有上升的趋势,但双频数据全部解算合格,所以可以认定基线解算不合格是由电离层延迟影响造成的。根据表 1中的统计数据, E方向单频结果浮动的幅度较小,基线08-11长度最短,整体残差(用标准偏差(standard deviation, std)表示)最小(std=2.1 mm);而基线05-08长度最长,整体残差也最大(std=6.5 mm)。6条基线随着长度的增加残差也逐渐增大,存在较为显著的线性关系。所有基线双频结果E分量都较为平稳,整体残差很小(std<3 mm),精度较之单频结果平均提高50%。N方向单频结果浮动的幅度比E方向大,但与基线长度的线性相关性较差,双频解算结果整体平稳(std<3 mm),精度较之单频平均提高了58%。H方向不论单频还是双频的解算结果,浮动的幅度都比E、N方向大,且与基线长度有一定的线性相关性,双频结果精度比单频结果平均提高了40%。
表 1 基线坐标分量标准差、解算不合格率统计信息
Table 1. Statistics of Each Baseline's Coordinate Component STD and Solving Reject Ratio
基线 E/mm N/mm H/mm 三维/mm 解算不合格率(%) L1 Lc L1 Lc L1 Lc L1 Lc L1 Lc 08-11(2.5 km) 2.1 1.2 6.1 2.9 4.8 2.9 8.0 4.3 0.5 0.0 08-03(4.4 km) 3.8 2.8 4.2 2.4 5.7 3.7 8.0 5.2 0.7 0.0 05-03(5.9 km) 4.2 2.2 5.9 2.1 8.3 5.1 11.8 5.9 13.3 0.0 05-10(6.8 km) 4.7 2.7 6.4 1.5 9.2 5.9 12.2 6.6 8.4 0.0 05-11(8.6 km) 5.3 1.5 4.2 2.5 9.7 6.1 11.8 6.7 14.1 0.0 05-08(9.7 km) 6.5 1.8 5.0 1.4 11.3 5.2 14.0 5.6 26.0 0.0 -
为进一步研究电离层延迟对单频点定位结果的影响,用双频结果减去单频结果得到电离层延迟引起定位结果的偏差。将6条基线电离层导致E、N、H 3个方向上定位结果的偏差制作统计图, 如图 5。为将3个方向的偏移量绘于同一幅图中,将dN加上0.02 m,dH减去0.03 m。从图 5中可以看出,3个方向中E方向受电离层延迟影响最小,在基线较短(小于5 km)时,上下浮动整体上不超过8 mm,标准差小于3 mm(见表 2);基线长度增加,偏移量的浮动程度也随之增加(相关系数为0.96),标准差最大值达到了6.7 mm。表 2中最后一行为各方向偏差量与基线长度的相关系数。电离层延迟对N、H方向上的影响明显大于E方向,其中N方向的上下浮动的幅度为0~3 cm,标准差为3.9~6.0 mm,但与基线长度相关性很差(相关系数为0.01);H方向上下浮动的幅度为0~5 cm,标准差4.8~9.2 mm,整体随着基线长度增加而增大(相关系数为0.85)。
表 2 定位结果偏差的标准差统计
Table 2. STD Statistics of Positioning Result Bias Caused by Ionosphere
基线 E N H 三维 08-11(2.5 km) 2.2 mm 5.5 mm 4.9 mm 7.7 mm 08-03(4.4 km) 2.8 mm 3.9 mm 4.8 mm 6.8 mm 05-03(5.9 km) 4.1 mm 5.6 mm 6.1 mm 9.2 mm 05-10(6.8 km) 4.1 mm 6.0 mm 6.5 mm 9.7 mm 05-11(8.6 km) 5.0 mm 4.3 mm 6.8 mm 9.4 mm 05-08(9.7 km) 6.7 mm 5.1 mm 9.2 mm 12.4 mm 相关系数 0.96 0.01 0.91 0.85 -
电离层对定位结果的影响除了在空间域有一定的规律外,在时间域也呈现一定的规律。为探究电离层在时间域的变化规律,以基线GD05-GD08 N方向的偏差量为例,取每个月的最大值和平均值分别绘于图 6和图 7,图 6、7中的曲线为6次多项式拟合曲线。可以得出以下结论。
1) 电离层延迟对定位结果影响的最大值的月变化幅度较大,从11月的5 mm到9月的12 mm,增加了一倍多。而月平均值的变化幅度较小,增幅最大为2.4 mm。
2) 电离层延迟对定位结果影响的最大值和平均值的峰值均分布在3月和9月,极值分别达到12.2 mm和4.6 mm。
3) 从6阶多项式拟合的结果来看,电离层延迟的季节性变化趋势比较明显,夏冬季节的影响小于春秋季节。
-
观测值域改正方法是从削弱单频观测值电离层误差的角度出发,首先提取双频测站间的双差电离层延迟信息,然后通过一定的内插模型, 如线性组合模型(linear combination model, LCM)逐历元逐卫星内插得到单频点的电离层延迟,从而降低单频观测值的误差,提高解算精度。本文内插模型采用LCM模型,模型原理如下[13-14]。
LCM模型通过构建虚拟单差线性观测值,能够消除流动站的空间相关误差,如内插系数满足式(3)的3个条件[14],则轨道偏差可以完全消除,电离层、对流层误差也可以大幅削弱。
$$ \sum\limits_{i = 1}^n {{\alpha _i} = 1} ;\sum\limits_{i = 1}^n {{\alpha _i}\left( {{X_u} - {X_i}} \right)} = 0;\sum\limits_{i = 1}^n {\alpha ^2_i = \min} $$ (3) 其中,αi为内插系数; Xu和Xi分别为参考站和流动站的坐标分量。当参考站间的双差模糊度固定以后,可以构建以下公式:
$$ \nabla \mathit{\Delta} {\varphi _{u,n}} - \sum\limits_{i = 1}^{n - 1} {{\alpha _i} \cdot {V_{i,n}}} = \nabla \mathit{\Delta} {\rho _{u,n}} + \lambda \nabla \mathit{\Delta} {N_{u,n}} $$ (4) 其中,▽Δφu, n为双差相位观测值; ▽Δρu, n为双差伪距观测值; ▽ΔNu, n为双差模糊度; Vi, n为基准站生成的双差改正数:
$$ \begin{align} {V_{i,n}} =& \nabla \mathit{\Delta} {\varphi _{i,n}} - \nabla \mathit{\Delta} {\rho _{i,n}} - \lambda \nabla \mathit{\Delta} {N_{i,n}}\\ &\left( {i = 1 \ldots n - 1} \right) \end{align} $$ (5) 上述模型中,获得精确的双差改正数的关键是正确固定模糊度,本文模糊度固定方法采用三步法[15],周跳探测与修复采用三差法[16]。由于该方法需要直接改正电离层延迟误差,所以需要对现有的软件进行二次开发或自编软件进行处理,本文采用自编GPS数据处理软件GPSDMS进行解算,数据解算流程见图 8。
-
为分析电离层改正效果,单频点选取GD07点,双频点选用GD03、GD08、GD10,分别以GD05为基准点进行基线解算。从图 9及表 3可以看出,经过双频点提供的改正信息改正之后,单频点的定位结果有较大改善,改善比例为57%。由于GD07点安装的是单频接收机,无法接收双频信息,所有无法对比改正后的结果与双频结果,但参考表 1可以看出,改正后的效果已经基本接近双频定位结果。
表 3 两种改正方案残差统计
Table 3. Statistics of Two Correction Scenarios Residual
残差 E/mm N/mm H/mm 三维/mm 改善比例(%) 单频未改正 5.2 3.9 8.4 10.6 观测值域改正 2.4 1.8 3.4 4.5 57 -
单频点的电离层改正是单双频GPS混合地面沉降监测的关键,本文从坐标域探讨了电离层对单频GPS点定位结果的影响,以及该影响与基线长度和季节的关系, 发现电离层延迟对高程方向的影响最大,南北方向次之,东西方向最小。随着基线长度的增加,电离层延迟影响增强,其中东西方向和南北方向的点位偏差量与基线长度呈较强的线性相关性;电离层延迟对定位结果的影响呈明显的季节性变化规律,夏冬季节强于春秋季节。采用观测值域改正方案对单频点的电离层延迟进行改正,能改善定位单频点的定位结果。
Influence of Ionosphere Delay on Single Frequency GPS Point and its Correction Method
-
摘要: 电离层延迟是单频GPS地面沉降监测点的最主要误差源,如何削弱该误差的影响是提高单双频混合地面沉降监测系统精度的关键。采用中国广州南沙单双频混合GPS地面沉降监测网数据进行处理分析,结果表明,在低纬度地区即使基线较短,电离层延迟对单频GPS监测精度影响仍然较大,影响程度随基线长度增加而增大,而且在时间域上有明显的季节性变化规律:2、3月份与8、9月份电离层影响较为显著。利用双频点数据从观测值域对单频点电离层延迟误差进行改正,监测精度提高了57%,改善效果明显。Abstract: Ionosphere delay is the dominant error source of single frequency GPS ground deformation monitoring point coordinate precision, how to reduce the error source influence is the key point to enhance the precision of single and dual frequency mixed ground deformation monitoring system. Based on the data of Guangzhou Nansha single and dual frequency mixed GPS ground deformation monitoring network, analysis indicates that ionosphere delay influence is still large, even though the baseline length is short in low latitude region. The impact of ionosphere delay aggravates as baseline length increases, and the rule of seasonal change is quite obvious:significant in February, March, August and September. Dual frequency data is used to correct the ionosphere delay error in observation value domain, as a result, the precision is improved by 57%, indicating that the effect of observation domain correction is obvious.
-
表 1 基线坐标分量标准差、解算不合格率统计信息
Table 1. Statistics of Each Baseline's Coordinate Component STD and Solving Reject Ratio
基线 E/mm N/mm H/mm 三维/mm 解算不合格率(%) L1 Lc L1 Lc L1 Lc L1 Lc L1 Lc 08-11(2.5 km) 2.1 1.2 6.1 2.9 4.8 2.9 8.0 4.3 0.5 0.0 08-03(4.4 km) 3.8 2.8 4.2 2.4 5.7 3.7 8.0 5.2 0.7 0.0 05-03(5.9 km) 4.2 2.2 5.9 2.1 8.3 5.1 11.8 5.9 13.3 0.0 05-10(6.8 km) 4.7 2.7 6.4 1.5 9.2 5.9 12.2 6.6 8.4 0.0 05-11(8.6 km) 5.3 1.5 4.2 2.5 9.7 6.1 11.8 6.7 14.1 0.0 05-08(9.7 km) 6.5 1.8 5.0 1.4 11.3 5.2 14.0 5.6 26.0 0.0 表 2 定位结果偏差的标准差统计
Table 2. STD Statistics of Positioning Result Bias Caused by Ionosphere
基线 E N H 三维 08-11(2.5 km) 2.2 mm 5.5 mm 4.9 mm 7.7 mm 08-03(4.4 km) 2.8 mm 3.9 mm 4.8 mm 6.8 mm 05-03(5.9 km) 4.1 mm 5.6 mm 6.1 mm 9.2 mm 05-10(6.8 km) 4.1 mm 6.0 mm 6.5 mm 9.7 mm 05-11(8.6 km) 5.0 mm 4.3 mm 6.8 mm 9.4 mm 05-08(9.7 km) 6.7 mm 5.1 mm 9.2 mm 12.4 mm 相关系数 0.96 0.01 0.91 0.85 表 3 两种改正方案残差统计
Table 3. Statistics of Two Correction Scenarios Residual
残差 E/mm N/mm H/mm 三维/mm 改善比例(%) 单频未改正 5.2 3.9 8.4 10.6 观测值域改正 2.4 1.8 3.4 4.5 57 -
[1] Janssen V, Roberts C, Rizos C, et al. Experiences with A Mixed-Mode GPS-Based Volcano Monitoring System at Mt. Papandayan, Indonesia[J]. Geomatics Research Australasia, 2001, 74:43-58 https://eprints.utas.edu.au/488/1/janssen_etal_GRA_2001.pdf [2] 谢荣安, 陈玉林, 戴吾蛟, 等. 单双频GPS混合地面沉降自动化监测方案设计[J]. 工程勘察, 2013(12): 48-52 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=gckc201312011&dbname=CJFD&dbcode=CJFQ Xie Rong'an, Chen Yulin. Scheme Design of Ground Subsidence Monitoring Using Mixed Single-and Dual-Frequency GPS Receivers[J]. Geotechnical Investigation & Surveying, 2013(12): 48-52 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=gckc201312011&dbname=CJFD&dbcode=CJFQ [3] 蔡昌盛, 高井祥, 李征航.利用GPS监测电离层总电子含量的季节性变化[J].武汉大学学报·信息科学版, 2006, 31(5):451-453 http://ch.whu.edu.cn/CN/abstract/abstract2462.shtml Cai Changsheng, Gao Jingxiang, Li Zhenghang. Monitoring Seasonal Variations of Ionospheric TEC Using GPS Measurements[J]. Geomatics and Information Science of Wuhan University, 2006, 31(5):451-453 http://ch.whu.edu.cn/CN/abstract/abstract2462.shtml [4] 姚宜斌, 陈家君, 陈鹏, 等. 2003~2006年磁暴期间欧洲区域电离层三维层析及演变分析[J].武汉大学学报·信息科学版, 2014, 39(2):132-136 http://ch.whu.edu.cn/CN/abstract/abstract2944.shtml Yao Yibin, Chen Jiajun, Chen Peng, et al. Analysis of Europe Ionospheric Responses During Magnetic Storms in 2003-2006 Using Ionospheric Tomographic Technology[J]. Geomatics and Information Science of Wuhan University, 2014, 39(2):132-136 http://ch.whu.edu.cn/CN/abstract/abstract2944.shtml [5] 柳景斌, 王泽民, 王海军, 等.利用球冠谐分析方法和GPS数据建立中国区域电离层TEC模型[J].武汉大学学报·信息科学版, 2008, 33(8):792-795 http://ch.whu.edu.cn/CN/abstract/abstract1682.shtml Liu Jingbin, Wang Zemin, Wang Haijun, et al. Modeling Regional Ionosphere Using GPS Measurements over China by Spherical Cap Harmonic Analysis Methodology[J]. Geomatics and Information Science of Wuhan University, 2008, 33(8):792-795 http://ch.whu.edu.cn/CN/abstract/abstract1682.shtml [6] 涂锐, 张勤, 黄观文, 等.利用相位平滑伪距和最小二乘曲面函数建立西安市区域电离层延迟模型[J].武汉大学学报·信息科学版, 2011, 36(2):218-221 http://ch.whu.edu.cn/CN/abstract/abstract471.shtml Tu Rui, Zhang Qin, Huang Guanwen, et al. Establishment of Regional Ionospheric Delay Model in Xi'an Based on Phase Smoothing Pseudo-Range and LS Surface Function[J]. Geomatics and Information Science of Wuhan University, 2011, 36(2):218-221 http://ch.whu.edu.cn/CN/abstract/abstract471.shtml [7] 张小红, 李征航, 蔡昌盛.用双频GPS观测值建立小区域电离层延迟模型研究[J].武汉大学学报·信息科学版, 2011, 26(2):140-143 http://ch.whu.edu.cn/CN/abstract/abstract5157.shtml Zhang Xiaohong, Li Zhenghang, Cai Changsheng. Study on Regional Ionospheric Model Using Dual-Frequency GPS Measurements[J]. Geomatics and Information Science of Wuhan University, 2011, 26(2):140-143 http://ch.whu.edu.cn/CN/abstract/abstract5157.shtml [8] 席超, 蔡成林, 李思敏, 等.一种基于三角分区的广域电离层改正新方法[J].武汉大学学报·信息科学版, 2015, 40(3):390-394 http://ch.whu.edu.cn/CN/abstract/abstract3216.shtml Xi Chao, Cai Chenglin, Li Simin, et al. A New Wide Area Ionospheric Correction Method Based on Triangular Partition[J]. Geomatics and Information Science of Wuhan University, 2015, 40(3):390-394 http://ch.whu.edu.cn/CN/abstract/abstract3216.shtml [9] 涂锐, 黄观文, 张勤, 等. GPS单频机电离层延迟改正新算法[J].武汉大学学报·信息科学版, 2012, 37(6):667-670 http://ch.whu.edu.cn/CN/abstract/abstract230.shtml Tu Rui, Huang Guanwen, Zhang Qin, et al. A New Algorithm on Ionospheric Delay Correction for Single Frequency GPS Receivers[J]. Geomatics and Information Science of Wuhan University, 2012, 37(6):667-670 http://ch.whu.edu.cn/CN/abstract/abstract230.shtml [10] 施闯, 辜声峰, 耿长江, 等.顾及电离层改正精度信息的高精度单频单点定位方法研究[J].武汉大学学报·信息科学版, 2011, 36(7):807-810 http://ch.whu.edu.cn/CN/abstract/abstract600.shtml Shi Chuang, Gu Shengfeng, Geng Changjiang, et al. High-Precision Single-Frequency Point Positioning with Randomness of Ionosphere Delay Correction in Consideration[J]. Geomatics and Information Science of Wuhan University, 2011, 36(7):807-810 http://ch.whu.edu.cn/CN/abstract/abstract600.shtml [11] 唐卫明, 刘经南, 刘晖, 等.一种GNSS网络RTK改进的综合误差内插方法[J].武汉大学学报·信息科学版, 2007, 32(12):1156-1159 http://ch.whu.edu.cn/CN/abstract/abstract2046.shtml Tang Weiming, Liu Jingnan, Liu Hui, et al. A Modified Combined Bias Interpolation Method for GNSS Network RTK[J]. Geomatics and Information Science of Wuhan University, 2007, 32(12):1156-1159 http://ch.whu.edu.cn/CN/abstract/abstract2046.shtml [12] 刘西凤, 袁运斌, 霍星亮, 等.电离层二阶项延迟对GPS定位影响的分析模型与方法[J].科学通报, 2010, 55(12):1162-1167 https://mall.cnki.net/qikan-DKXB201501022.html Liu Xifeng, Yuan Yunbin, Huo Xingliang, et al. Model Analysis Method(MAM) on the Effect of the Second-Order Ionospheric Delay on GPS Positioning Solution[J].Chinese Sci Bull, 2010, 55(12):1162-1167 https://mall.cnki.net/qikan-DKXB201501022.html [13] Han S, C Rizos. GPS Network Design and Error Mitigation for Real-Time Continuous Array Monitoring Systems[C]. The 9th Int. Tech. Meeting of the Satellite Div. of the U. S. ION, Kansas City, Missouri, 1996 https://core.ac.uk/display/21122947 [14] Dai L, Han S. Comparison of Interpolation Algorithms in Network-Based GPS Techniques[J]. Navigation, 2003, 50(4):277-293 doi: 10.1002/navi.2003.50.issue-4 [15] 唐卫明, 刘经南, 施闯, 等.三步法确定网络RTK基准站双差模糊度[J].武汉大学学报·信息科学版, 2007, 32(4):305-308 http://ch.whu.edu.cn/CN/abstract/abstract1875.shtml Tang Weiming, Liu Jingnan, Shi Chuang, et al. Three Steps Method to Determine Double Difference Ambiguities Resolution of Network RTK Reference Station[J]. Geomatics and Information Science of Wuhan University, 2007, 32(4):305-308 http://ch.whu.edu.cn/CN/abstract/abstract1875.shtml [16] 王贵文, 王泽民, 殷海涛.基于三差观测量的实时动态GPS周跳修复方法研究[J].武汉大学学报·信息科学版, 2007, 32(8):711-714 http://ch.whu.edu.cn/CN/abstract/abstract1972.shtml Wang Guiwen, Wang Zemin, Yin Haitao. An Cycle-Slip Correction Method for Real-Time Kinematic GPS Data Based on Triple Differences Observation[J]. Geomatics and Information Science of Wuhan University, 2007, 32(8):711-714 http://ch.whu.edu.cn/CN/abstract/abstract1972.shtml -