BDS网络RTK中距离参考站整周模糊度单历元解算方法

祝会忠，徐爱功，高　猛， 杨秋实

辽宁工程技术大学测绘与地理科学学院，辽宁 阜新 123000

The Algorithm of Single-epoch Integer Ambiguity Resolution between Middle-range BDS Network RTK Reference Stations

ZHU Huizhong，XU Aigong，GAO Meng，YANG Qiushi

School of Geomatics，Liaoning Technical University，Fuxin 123000，China



BDS网络RTK中距离参考站整周模糊度单历元解算方法

祝会忠，徐爱功，高猛， 杨秋实

辽宁工程技术大学测绘与地理科学学院，辽宁 阜新 123000

Foundation support： The National Natural Science Foundation of China (Nos.41504010；41474020)；The National High-tech Research and Development Program of China (863 Program) (No.2014AA121501)；The Program for Liaoning Doctoral Foundation (No.20141141)；The Open Research Foundation of Science and Technology on Aerospace Flight Dynamics Laboratory(No. 2014afdl005)

摘要：提出了一种BDS网络RTK中距离(50～100km)参考站间的双频载波相位整周模糊度单历元解算方法。该方法首先利用B1、B2载波相位整周模糊度间的线性关系选取B1、B2载波相位整周模糊度备选值。利用双频载波相位整周模糊度备选值计算双差电离层延迟误差，根据参考站各卫星电离层延迟误差间的空间关系，使用双差电离层延迟误差构建双差电离层延迟误差的线性计算模型。通过双差电离层延迟误差线性计算模型的建立搜索和确定B1、B2载波相位的整周模糊度。经CORS网实测数据试验算例的验证，该方法只需一个历元的观测数据即可确定参考站间双差B1、B2载波相位整周模糊度，且不受周跳影响。

关键词：BDS；网络RTK；单历元；整周模糊度；电离层延迟误差

北斗导航卫星系统(BeiDou navigation satellite system，BDS)是我国具有独立自主产权、自主研发的卫星导航定位系统[1]。目前，北斗导航卫星系统空间星座由5颗地球同步轨道(geo-synchronous orbit，GEO)卫星、5颗倾斜地球同步轨道(inclined geo-synchronous orbits，IGSO)卫星、4颗中高度轨道(medium earth orbit，MEO)卫星组成[2]。2012年12月27日，北斗系统空间信号接口控制文件(ICD)正式版公布，北斗导航卫星系统正式对亚太地区提供无源定位、导航、授时服务。BDS提供的标准服务定位精度为10 m，测速精度为0.2 m/s，授时精度为20 ns[3]。BDS网络RTK方法是提高BDS标准定位精度的重要手段之一。

实现BDS高精度定位都是采用载波相位观测数据[4]。而利用BDS网络RTK方法进行高精度动态定位的首要条件是参考站间载波相位整周模糊度的准确确定，只有正确固定了参考站模糊度才能得到高精度的综合误差或是建立高精度的误差模型。由于网络RTK参考站间距离较长，一般在几十千米以上，采用简单双差组合的方法单历元确定BDS载波相位整周模糊度非常困难，主要是由于电离层延迟和对流层延迟等误差对双差观测值的影响大于模糊度的半个波长。即使在使用双频观测数据和参考站坐标已知的情况下，整周模糊度也难以与误差分离。

目前，国内外学者对GPS网络RTK参考站双差模糊度解算作了大量研究，并取得了很多成果[5-16]。文献[5]研究了GPS长距离静态定位模糊度解算算法；文献[6]使用卡尔曼滤波方法实现了参考站间模糊度的快速解算；文献[7]研究了纳伪概率可控的四舍五入方法进行GPS参考站间的模糊度固定，可提高模糊度整数解的可靠性。由于GPS系统发展比较成熟，研究者们研究和实现了GPS网络RTK参考站整周模糊度的单历元确定[8-9]，文献[8]提出了GPS网络RTK参考站间单历元模糊度搜索方法。但BDS刚刚建立，其信号质量、接收设备等相对于GPS不够完善，同时，BDS星座中包含GEO卫星和IGSO卫星，广播星历的卫星轨道精度相对于GPS较低。所以，利用GPS网络RTK参考站整周模糊度解算方法不能很好地实现BDS参考站的整周模糊度单历元解算。

因此，本文针对B1、B2双频载波相位观测数据，提出一种BDS网络RTK参考站间的双频载波相位整周模糊度单历元解算方法。该方法利用B1、B2载波相位整周模糊度间的线性关系选择载波相位整周模糊度备选值，然后通过参考站间双差电离层延迟误差线性计算模型的建立搜索和确定B1、B2载波相位的整周模糊度。

1参考站B1、B2载波相位整周模糊度间的线性关系

BDS系统B1、B2的频率分别为1 561.098 Mhz、1 207.140 Mhz，B1频率与GPS的L1频率接近，B2频率与Galileo E5b频率相同。因为双差组合观测值具有可消除接收机钟差和卫星钟差，大大削弱卫星星历误差、电离层延迟误差、对流层延迟误差等误差影响的优点，双差载波相位模糊度可保持整周特性，所以，本文对BDS载波相位双差模糊度进行单历元解算。通常情况下，BDS网络RTK的参考站都设在比较开阔的地方，卫星的多路径效应较小，可忽略多路径的影响，则参考站的双差B1、B2载波相位观测方程可表示为

(1)

(2)

(3)

2参考站B1、B2载波相位整周模糊度的单历元确定

2.1利用模糊度备选值计算双差电离层延迟误差

根据参考站各BDS卫星电离层延迟误差的空间相关关系，可以建立双差电离层延迟误差的空间计算模型，通过构建双差电离层延迟误差的空间计算模型搜索和确定B1、B2载波相位的整周模糊度。正确的B1、B2载波相位整周模糊度备选值组合可以计算出正确的电离层延迟误差，而错误的B1、B2载波相位整周模糊度组合计算出的电离层延迟误差存在偏差。可利用双频载波相位模糊度备选值组合得到对应的双差电离层延迟误差，然后利用这些双差电离层延迟误差构建多个双差电离层延迟误差模型，其中只有正确的双差电离层延迟误差建立的误差模型才是准确的，即可实现对B1、B2载波相位整周模糊度的搜索。

(4)

(5)

2.2利用双差电离层延迟误差模型搜索模糊度

根据两个参考站电离层延迟误差的空间相关性，可以建立BDS双差电离层延迟误差的空间线性模型如式(6)所示

(6)

如果当前参考站的共视卫星有n颗(n-1>3)，则对于参考站BDS系统B1载波相位的电离层延迟误差，其计算双差电离层延迟误差的线性计算模型有

(7)

以矩阵的形式表示为

li=hm

(8)

利用2.1节中得到的各双差卫星整周模糊度备选值对应的双差电离层延迟误差，计算参考站双差电离层延迟误差线性计算模型的系数，根据最小二乘解算模型系数的计算确定正确的模糊度备选值。将参考站每颗双差卫星模糊度备选值对应的B1、B2载波相位双差电离层延迟误差作为待检测对象。对各颗BDS双差卫星相同频率的待检测电离层延迟误差数值进行组合，组成多组n-1维的双差电离层延迟误差向量，利用各组待检测双差电离层延迟误差向量计算双差电离层延迟误差的线性计算模型系数。根据式(8)，li即为n-1维的双差电离层延迟误差向量，则计算载波相位双差电离层延迟误差模型系数m的误差方程为

Vi=hm-li

(9)

根据最小二乘计算原理，如果进行等权计算，则有

(10)

对于式(10)，也可根据卫星高度角定权，进行加权计算。将式(10)平差计算出的双差电离层延迟误差线性计算模型的系数代入式(9)，可得

(11)

分别对B1、B2载波相位双差电离层延迟误差组合进行上述计算和搜索，确定出满足要求的双差电离层延迟误差，其对应的B1、B2载波相位整周模糊度备选值则为正确的模糊度备选值。

2.3B1、B2载波相位整周模糊度的检验

将根据B1载波相位双差电离层延迟误差模型搜索出的B1、B2载波相位整周模糊度备选值，与根据B2载波相位双差电离层延迟误差模型搜索出的B1、B2载波相位整周模糊度备选值进行比较，如果二者一致则说明搜索出的B1、B2载波相位整周模糊度为正确的整周模糊度组合。

此外，对于任意两个以上的参考站所组成的闭合基线，双差整周模糊度的代数和在理论上为零。以参考站A、B和C为例，则

(12)

(13)

式中，下标表示参考站。将闭合基线的双差B1、B2载波相位整周模糊度代入式(12)、式(13)进行检验，如果满足该闭合条件则认为被检验的双差B1、B2载波相位整周模糊度是正确的。同时，3条基线同一双差卫星的3个双差电离层延迟误差的代数和也为零，此条件也可用来对解算出的载波相位整周模糊度进行检验。

3算例与分析

本文使用2013年4月29日(BDS时382周，起始时刻周内秒为86 400)华中地区CORS网一天的BDS实测数据进行算法检验，采样间隔为1 s。该试验数据共有3个参考站，测站分布如图1所示，stationA到stationB距离53 km，stationA到stationC距离102 km，stationB到stationC距离81 km。

图1　测站分布Fig.1　Location of stations

表1　stationA—stationB双差模糊度备选值

本文一天的试验数据算法验证中，分别以卫星C03-C07、C09-C07双差卫星为例。卫星C03-C07高度角变化如图2所示，卫星C03-C07电离层穿刺点的经纬度差值如图3所示。C03为GEO卫星，其高度角变化很小，C07为IGSO卫星运动周期约为24 h。在截止高度角选择为15°的情况下，该组一天的试验数据中参考站连续观测卫星的时间达到约14.5 h。由于参考站观测到的GEO卫星电离层穿刺点变化较小，所以，电离层穿刺点的经纬度差值的变化，是C07的穿刺点经纬度变化决定的。IGSO卫星的运行轨道星下点轨迹以子午线为中心，参考站观测到的C07电离层穿刺点经度变化较小，纬度变化较大。所以，卫星C03-C07电离层穿刺点的经纬度差值中，经度方向的差值变化较慢，纬度方向的差值变化稍大。利用本文的方法将B1、B2载波相位整周模糊度准确确定后，卫星C03-C07的B1载波相位观测值双差电离层延迟误差值如图4所示。

图2　C03-C07高度角Fig.2　The elevations of C03-C07

图3　卫星C03-C07穿刺点经纬度差值Fig.3　The longitude and latitude difference of C03-C07 IPP

卫星C09-C07高度角变化如图5所示，卫星C09-C07电离层穿刺点的经纬度差值如图6所示。C09-C07均为IGSO卫星，运动周期约24 h，截止高度角选择为15°，该组一天的参考站试验数据中第1个连续观测到卫星C09的时段长约8 h。从图6中可知， 卫星C09-C07电离层穿刺点的经纬度差值的变化相对于卫星C03-C07较快。卫星C09-C07电离层穿刺点的经纬度差值中，经度方向的差值变化较慢，纬度方向的差值变化稍大。利用本文的方法将B1、B2载波相位整周模糊度准确确定后，卫星C09-C07连续8 h的B1载波相位观测值双差电离层延迟误差值如图7所示。

图4　卫星C03-C07双差电离层延迟误差Fig.4　Double difference ionosphere delay error of C03-C07

图5　卫星C09-C07的高度角Fig.5　The elevations of C09-C07

图6　卫星C09-C07穿刺点经纬度差值　Fig.6　The longitude and latitude difference of C09-C07 IPP

图7　卫星C09-C07双差电离层延迟误差Fig.7　Double difference ionosphere delay error of C09-C07

本文算例的观测时间从凌晨开始，随着观测值时间的推移，越接近中午，图4、图7中3条基线双差电离层延迟误差的整体数值呈逐渐增大的趋势。同时，双差的电离层延迟误差与卫星高度角以及两颗求双差卫星的电离层穿刺点的空间位置关系有关，即卫星高度角越大，电离层延迟误差越小，双差卫星穿刺点经纬度的差值绝对值越小(电离层穿刺点的空间位置越接近)，双差电离层延迟误差越小。因此，图4、图7中的B1载波相位观测值双差电离层延迟误差除了随观测时间变化外，还受卫星高度角和双差卫星电离层穿刺点经纬度差值的影响，双差电离层延迟误差有数值增减的抖动，即双差电离层延迟误差会出现小范围内减小或增大的情况。

图4、图7中的双差电离层延迟误差，满足闭合基线电离层延迟误差和等于零，即stationA—stationB、stationB—stationC、stationC—stationA上同一颗双差卫星的电离层延迟误差之和为零，stationC—stationA的电离层延迟误差由stationA—stationC的电离层延迟误差取相反数得到。表2为该算例3个参考站间双差载波相位整周模糊度固定的成功率(卫星截止高度角15°)。方法1为本文的方法，从表2中可以看出该算例利用本文方法的模糊度固定成功率较高。方法2为文献[8]中的单历元模糊度确定方法。由于BDS与GPS的波长不同，且BDS广播星历的轨道精度和伪距观测值的精度较GPS差，所以，利用方法2进行该算例BDS整周模糊度的单历元解算效率和成功率较低。方法3是利用MW组合计算宽巷模糊度，然后计算确定当前历元原始载波相位模糊度。MW组合计算宽巷模糊度的收敛过程需要一段观测时间，表2中的统计信息为宽巷模糊度固定之后单历元确定载波相位模糊度的数目。由于BDS伪距观测值的噪声较大，精度较低，所以，消耗了较多的观测时间用于确定宽巷模糊度，导致很多单历元模糊度不能在当前历元固定。同时，由于BDS系统的广播星历轨道精度较低，即使宽巷模糊度固定之后，原始载波相位模糊度的确定效率和成功率也较低。从表2中可以看出，本文的方法相比方法2、方法3具有较高的BDS载波相位整周模糊度单历元固定成功率。根据算法试验结果可知，本文的方法能够实现BDS中距离参考站载波相位整周模糊度的单历元准确解算，且具有较高成功率和效率。

表2　模糊度固定成功率

4结论

本文方法通过消除电离层延迟误差选取载波相位模糊度的备选组合，并构建双差电离层延迟误差模型进行模糊度搜索，可降低电离层延迟误差对参考站间整周模糊度解算的影响。利用该方法参考站网只需一个历元即可启动，实用性强，可靠性好。同时，本文利用当前历元电离层延迟误差的空间关系建立双差电离层延迟误差改正模型，并且是进行参考站整周模糊度的单历元解算，因此，电离层延迟误差历元间的变化对本文的方法影响较小。由于GEO卫星相对于地球静止、IGSO卫星运动周期长，卫星的电离层穿刺点变化缓慢，历元之间的电离层延迟误差的相关性较强。因此，随着观测时间的延长，在没有电离层剧烈突变的情况下，可以根据电离层延迟误差的变化，对解算出的B1、B2载波相位整周模糊度进行可靠性判断。本文的方法是针对中距离BDS网络RTK参考站整周模糊度的单历元解算，对于长距离BDS网络RTK参考站整周模糊度的单历元解算方法还需进一步研究。

参考文献：

[1]杨元喜. 北斗卫星导航系统的进展、贡献与挑战[J]. 测绘学报, 2010, 39(1): 1-6.

YANG Yuanxi. Progress, Contribution and Challenges of Compass-BeiDou Satellite Navigation System[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(1): 1-6.

[2]杨元喜, 李金龙, 徐君毅, 等. 中国北斗卫星导航系统对全球PNT用户的贡献[J]. 科学通报, 2011, 56(21): 1734-1740.

YANG Yuanxi, LI Jinlong, XU Junyi, et al. Contribution of the Compass Satellite Navigation System to Global PNT Users[J]. Chinese Science Bulletin, 2011, 56(26): 2813-2819.

[3]SHI Chuang, ZHAO Qile, LI Min, et al. Precise Orbit Determination of BeiDou Satellites with Precise Positioning[J]. Science China Earth Sciences, 2012, 55(7): 1079-1086.

[4]高星伟, 过静珺, 程鹏飞, 等. 基于时空系统统一的北斗与GPS融合定位[J]. 测绘学报, 2012, 41(5): 743-748, 755.

GAO Xingwei, GUO Jingjun, CHENG Pengfei, et al. Fusion Positioning of BeiDou-GPS Based on Spatio Temporal System Unification[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(5): 743-748, 755.

[5]VOLLATH U, LANDAU H, CHEN Xiaoming. Network RTK——Concept and Performance[C]∥Proceedings of the GNSS Symposium. Wuhan: [s.n.], 2002.

[6]周乐韬, 黄丁发, 袁林果, 等. 网络RTK参考站间模糊度动态解算的卡尔曼滤波算法研究[J]. 测绘学报, 2007, 36(1): 37-42.

ZHOU Letao, HUANG Dingfa, YUAN Linguo, et al. A Kalman Filtering Algorithm for Online Integer Ambiguity Resolution in Reference Station Network[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(1): 37-42.

[7]李博峰, 沈云中, 张兴福. 纳伪概率可控的四舍五入法及其在RTK模糊度固定中的应用[J]. 测绘学报, 2012, 41(4): 483-489, 495.

LI Bofeng, SHEN Yunzhong, ZHANG Xingfu. Error Probability Controllable Integer Rounding Method and Its Application to RTK Ambiguity Resolution[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4): 483-489, 495.

[8]高星伟, 刘经南, 葛茂荣. 网络RTK基准站间基线单历元模糊度搜索方法[J]. 测绘学报, 2002, 31(4): 305-309.

GAO Xingwei,LIU Jingnan,GE Maorong.An Ambiguity Searching Method for Network RTK Baselines between Base Stations at Single Epoch[J]. Acta Geodaetica et Cartographica Sinica, 2002, 31(4): 305-309.

[9]祝会忠, 刘经南, 唐卫明, 等. 长距离网络RTK基准站间整周模糊度单历元确定方法[J]. 测绘学报, 2012, 41(3): 359-365.

ZHU Huizhong, LIU Jingnan, TANG Weiming, et al. The Algorithm of Single-epoch Integer Ambiguity Resolution between Long-range Network RTK Base Stations[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 359-365.

[10]高星伟, 陈锐志, 赵春梅. 网络RTK算法研究与实验[J]. 武汉大学学报(信息科学版), 2009, 34(11): 1350-1353.

GAO Xingwei, CHEN Ruizhi, ZHAO Chunmei. A Network RTK Algorithm and Its Test[J]. Geomatics and Information Science of Wuhan University, 2009, 34(11): 1350-1353.

[11]BLEWITT G. Carrier Phase Ambiguity Resolution for the Global Positioning System Applied to Geodetic Baselines up to 2000 km[J]. Journal of Geophysical Research, 1989, 94(B8): 10187-10203.

[12]DONG Danan, BOCK Y. Global Positioning System Network Analysis with Phase Ambiguity Resolution Applied to Crustal Deformation Studies in California[J]. Journal of Geophysical Research, 1989, 94(B4): 3949-3966.

[13]BLEWITT G. An Automatic Editing Algorithm for GPS Data[J]. Geophysical Research Letters, 1990, 17(3): 199-202.

[14]SUN Huangqi, CANNON M E, MELGARD T E. Real-time GPS Reference Network Carrier Phase Ambiguity Resolution[C]∥Proceedings of the National Technical Meeting of the Institute of Navigation. San Diego: Institute of Navigation, 1999: 193-199.

[15]HAN Shaowei. Carrier Phase-based Long-range GPS Kinematic Positioning[D]. Sydney: University of New South Wales, 1997.

[16]DAI Liwen, WANG Jinling, RIZOS C, et al. Predicting Atmospheric Biases for Real-time Ambiguity Resolution in GPS-GLONASS Reference Station Networks[J]. Journal of Geodesy, 2003, 76(11-12): 617-628.

[17]葛茂荣, 刘经南. GPS定位中对流层折射估计研究[J]. 测绘学报, 1996, 25(4): 285-291.

GE Maorong, LIU Jingnan. The Estimation Methods for Tropospheric Delays in Global Positioning System[J]. Acta Geodaetica et Cartographica Sinica, 1996, 25(4): 285-291.

[18]高星伟, 陈锐志, 李夕银. 中性大气对非差伪距定位的影响及其模型改正分析[J]. 测绘学报, 2007, 36(2): 134-140.

GAO Xingwei, CHEN Ruizhi, LI Xiyin. A Study of the Neutral Atmospheric Effects on Pseudorange Positioning with Non-difference Models[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(2): 134-140.

[19]HOPFIELD H S. Two-quartic Tropospheric Refractivity Profile for Correcting Satellite Data[J]. Journal of Geophysical Research, 1969, 74(18): 4487-4499.

[20]SAASTAMOINEN J. Contributions to the Theory of Atmospheric Refraction[J]. Bulletin Géodésique, 1973, 107(1): 13-34.

[21]BLACK H D. An Easily Implemented Algorithm for the Tropospheric Range Correction[J]. Journal of Geophysical Research, 1978, 83(B4): 1825-1828.

(责任编辑：丛树平)

修回日期： 2015-08-30

First author： ZHU Huizhong (1983—), male, PhD, majors in the algorithm and application of GNSS for high-accuracy dynamic positioning.

E-mail： zhuhuizhong@whu.edu.cn

The Algorithm of Single-epoch Integer Ambiguity Resolution between Middle-range BDS Network RTK Reference Stations

ZHU Huizhong，XU Aigong，GAO Meng，YANG Qiushi

School of Geomatics，Liaoning Technical University，Fuxin 123000，China

Abstract：The algorithm of double frequency carrier phase integer ambiguity resolution between middle-range reference stations at single-epoch is proposed. The ambiguity candidates of B1 frequency and B2 frequency carrier phase are selected by linear relationship between carrier phase ambiguities of B1 frequency and B2 frequency carrier phase. The double difference ionosphere delay errors can be calculated by the ambiguity candidates of double frequency carrier phase. Then the linear calculation model of double difference ionosphere delay errors can be established by the double difference ionosphere delay errors according to the spatial relationship between the ionosphere delay errors of each satellites on reference stations. And the ambiguities of B1 frequency and B2 frequency carrier phase are searched and fixed in the establishing of the linear calculation model of double difference ionosphere delay errors. This algorithm is tested by the data of BDS CORS network. The results indicate that the double difference integer ambiguity between reference stations can be fixed at single-epoch, and this algorithm avoids the effect of cycle clips.

Key words：BDS；network RTK； single-epoch；integer ambiguity；ionosphere delay error

作者简介：第一 祝会忠(1983—)，男，博士，研究方向为高精度GNSS定位算法与应用。

收稿日期：2014-10-13

基金项目：国家自然科学基金(41504010；41474020)；国家863计划(2014AA121501)；辽宁省博士启动基金(20141141)；航天飞行动力学技术重点实验室开放基金(2014afdl005)

中图分类号：P228

文献标识码：A

文章编号：1001-1595(2016)01-0050-08

引文格式：祝会忠，徐爱功，高猛，等.BDS网络RTK中距离参考站整周模糊度单历元解算方法[J].测绘学报，2016,45(1)：50-57.DOI:10.11947/j.AGCS.2016.20140525.

ZHU Huizhong，XU Aigong，GAO Meng，et al.The Algorithm of Single-epoch Integer Ambiguity Resolution between Middle-range BDS Network RTK Reference Stations[J]. Acta Geodaetica et Cartographica Sinica,2016,45(1):50-57.DOI:10.11947/j.AGCS.2016.20140525.