Key Technology on Establishment of an Independent CoordinateSystem for Water Diversion Project in Central Yunnan
-
摘要: 目的 云南滇中引水工程具有测区跨度长、平均海拔高、高差变化大等特点,独立坐标系统的建立难度大。本文根据投影变形理论,结合滇中引水工程的特点,提出了适用性的建立独立坐标系统的关键技术解决方案。用1°带的分带方式替代常规的国家3°带分带,有效地控制工程的高斯投影变形;推导了新的适用于高原长距离工程的抵偿投影面求取方法,并结合滇中引水工程实际给出了优化的抵偿投影面设置方案。使用本文建议的独立坐标系统,工程综合变形的最大值能得到最优控制。Abstract: Objective A water diversion project in Central Yunnan exhibited several typical characteristics of theoverall survey region,such as the extraordinarily long span,the high average altitude,and large ele-vation changes.Therefore,difficult to build independent coordinate system for this project.Based onthe theory of projection distortion,we propose a key technology solutions composed of two parts.First,using the 1°projection zone to replace the regular 3°projection zone,Gauss projection distor-tion was effectively controlled.Second,a method for selecting a compensation projection plane wasdeducted to satisfy the needs of this water diversion project.An optimized solution for this project wasfound using the proposed independent coordinate system and efficiently controlled,the maximum com-prehensive deformation for this water diversion project in Central Yunnan.
-
全球定位系统（global positioning system，GPS）是目前最成熟、应用最广泛的全球卫星导航系统，北斗卫星导航系统（BeiDou satellite navigation system，BDS）也建成了覆盖亚太地区的区域导航定位系统，并于2018-12-27初步建成三代基本系统，开始提供全球服务。GPS/BDS融合定位与单系统相比，参与定位解算的卫星更多，卫星几何构型更强。多频观测数据可以提升载波相位整周模糊度的解算性能[1]，从而提高定位的稳定性。GPS/BDS双系统、单系统网络实时动态定位（real-time kinematic，RTK）技术已经成为实时动态定位领域的热点问题。
国内外学者对GPS网络RTK定位和全球卫星导航系统（global navigation satellite systems，GNSS）模糊度解算进行了大量研究[2-13]。文献[2]针对单系统GPS网络RTK载波相位整周模糊度解算、参考站网误差内插、流动站载波相位整周模糊度解算模型进行了详细阐述；文献[14]采用站间距离为数十公里的GPS/BDS等多系统参考站网对比分析了多系统RTK的定位精度；文献[15]考虑了中长距离基线误差的影响，提出了自适应抗差滤波算法；文献[16]提出了BDS三频网络RTK参考站载波相位整周模糊度的解算方法；文献[17]实现了BDS/GPS双系统参考站载波相位整周模糊度的快速解算。相对于参考站间距为几十公里的网络RTK方法，参考站间距离大于100 km的长距离网络RTK可以扩大作业面积，节省成本。电离层延迟、对流层延迟等误差的空间相关性受长距离网络RTK参考站间距离的影响较大，使得整周模糊度与双差残余误差的分离难度增加，参考站网误差改正模型的建立相对困难。目前，有关GPS/BDS双系统长距离网络RTK定位算法的研究较少。
本文利用长距离参考站网GPS/BDS多频伪距和载波相位观测值，研究了一种长距离GPS/BDS网络RTK方法。首先对参考站网的宽巷模糊度和多频载波相位整周模糊度进行解算；然后利用固定的参考站网载波相位整周模糊度计算参考站网观测值的双差观测误差，根据长距离参考站网观测误差特性的不同，利用分类误差内插法分别计算流动站色散性误差（电离层延迟误差）和非色散性误差（对流层延迟误差为主），以改正流动站的观测误差；最后解算流动站GPS/BDS单系统、双系统的整周模糊度和位置参数，并采用实测CORS（continuous operating reference station）网数据进行算法验证和实验分析。
1 长距离参考站网载波相位整周模糊度解算
长距离GPS/BDS网络RTK方法首先要准确解算参考站网GPS/BDS载波相位的整周模糊度。参考站的坐标精确已知，则参考站网双差载波相位和伪距观测方程分别为[17]：
$$ \begin{array}{*{20}{l}} {{\lambda _i}^S \cdot {\rm{\Delta }}\nabla {\varPhi _i}^S = {\rm{\Delta }}\nabla {\rho ^S} - {\lambda _i}^S \cdot {\rm{\Delta }}\nabla {N_i}^S - }\\ \ \ \ \ {{\rm{\Delta }}\nabla {I_i}^S + {\rm{\Delta }}\nabla {T^S} + {\rm{\Delta }}\nabla {\varepsilon _i}^S} \end{array} $$ (1) $$ {\rm{\Delta }}\nabla {P_i}^S = {\rm{\Delta }}\nabla {\rho ^S} + {\rm{\Delta }}\nabla {I_i}^S + {\rm{\Delta }}\nabla {T^S} + {\rm{\Delta }}\nabla {\varepsilon _i}_P^S $$ (2) 式中，${\rm{\Delta }}\nabla $表示双差；$\varPhi $表示载波相位观测值；P表示伪距观测值；上标S表示卫星，GPS、BDS分别简写为G、B；下标i代表频率；ρ表示测站到卫星的几何距离；λ表示波长；N表示载波相位模糊度；I表示电离层延迟误差；T表示以对流层延迟为主的非色散性误差；ε代表载波相位观测噪声；εP代表伪距观测噪声。
首先解算GPS/BDS第一、二频率的宽巷整周模糊度，然后建立引入大气误差参数的载波相位整周模糊度解算模型。采用单系统分别估计的方式进行宽巷整周模糊度的确定。利用Melbourne-Wubbena（M-W）组合观测值不受大气误差、卫星轨道误差影响的特性，解算长距离参考站间的宽巷整周模糊度。各系统M-W组合观测值和宽巷模糊度可由式（1）、（2）推导得到，其分别表示为：
$$ \begin{array}{l} {\rm{\Delta }}\nabla {\rm{MW}}_{12}^G = \frac{{c \cdot {\rm{\Delta }}\nabla \varPhi _1^G - c \cdot {\rm{\Delta }}\nabla \varPhi _2^G}}{{f_1^G - f_2^G}} - \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \frac{{f_1^G \cdot {\rm{\Delta }}\nabla P_1^G + f_2^G \cdot {\rm{\Delta }}\nabla P_2^G}}{{f_1^G + f_2^G}} \end{array} $$ (3) $$ {\rm{\Delta }}\nabla N_{12}^G = \frac{{\left( {f_1^G - f_2^G} \right) \cdot {\rm{\Delta }}\nabla {\rm{MW}}_{12}^G}}{c} $$ (4) $$ \begin{array}{l} {\rm{\Delta }}\nabla {\rm{MW}}_{i, i + 1}^B = \frac{{c \cdot {\rm{\Delta }}\nabla \varPhi _i^B - c \cdot {\rm{\Delta }}\nabla \varPhi _{i + 1}^B}}{{f_i^B - f_{i + 1}^B}} - \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \frac{{f_i^B \cdot {\rm{\Delta }}\nabla P_i^B + f_{i + 1}^B \cdot {\rm{\Delta }}\nabla P_{i + 1}^B}}{{f_i^B + f_{i + 1}^B}} \end{array} $$ (5) $$ {\rm{\Delta }}\nabla N_{i, i + 1}^B = \frac{{\left( {f_i^B - f_{i + 1}^B} \right) \cdot {\rm{\Delta }}\nabla {\rm{MW}}_{i, i + 1}^B}}{c} $$ (6) 式中，${\rm{\Delta }}\nabla {\rm{MW}}_{12}^G$表示M-W组合观测值；${\rm{\Delta }}\nabla N_{12}^G$表示GPS宽巷模糊度；${f_1^G}$、${f_2^G}$表示GPS频率；${\rm{\Delta }}\nabla {\rm{MW}}_{i, i + 1}^B$代表BDS的第i、i+1频率形成的M-W组合观测值；${\rm{\Delta }}\nabla N_{i, i + 1}^B$表示BDS宽巷模糊度；${f_i^B}$、${f_{i + 1}^B}$表示BDS第i、i+1两个频率，其中i=1, 2；c为光速。其他符号含义与式（1）、（2）相同。
采用式（4）、式（6）计算宽巷整周模糊度。由于伪距观测噪声会影响宽巷整周模糊度的计算精度，因而采用对各历元模糊度的计算值取平均的方式提高宽巷模糊度的精度，并利用M-W组合观测值对GPS、BDS多频观测数据进行预处理，确保式（3）、式（5）观测值的连续性。各系统双差宽巷模糊度计算值和中误差分别为[18]：
$$ {\rm{\Delta }}\nabla \hat N_{12}^G = \left\langle {{\rm{\Delta }}\nabla N_{12}^G} \right\rangle $$ (7) $$ {\sigma _{{\rm{\Delta }}\nabla \hat N_{12}^G}} = \sqrt {\frac{{\mathop \sum \nolimits^ {{({\rm{\Delta }}\nabla N_{12}^G - {\rm{\Delta }}\nabla \hat N_{12}^G)}^2}}}{n}} $$ (8) $$ {\rm{\Delta }}\nabla \hat N_{i, i + 1}^B = \left\langle {{\rm{\Delta }}\nabla N_{i, i + 1}^B} \right\rangle $$ (9) $$ {\sigma _{{\rm{\Delta }}\nabla \hat N_{i, i + 1}^B}} = \sqrt {\frac{{\mathop \sum \nolimits^ {{({\rm{\Delta }}\nabla N_{i, i + 1}^B - {\rm{\Delta }}\nabla \hat N_{i, i + 1}^B)}^2}}}{n}} $$ (10) 式中， < > 表示求平均值；${\rm{\Delta }}\nabla \hat N_{12}^G$表示${{\rm{\Delta }}\nabla N_{12}^G}$的多历元平均值；${\rm{\Delta }}\nabla \hat N_{i, i + 1}^B$表示${{\rm{\Delta }}\nabla N_{i, i + 1}^B}$的多历元平均值；n是观测数据的历元数。双差宽巷整周模糊度根据双差宽巷模糊度计算值及其中误差确定，采用的宽巷整周模糊度确定准则为：双差宽巷模糊度计算值与最接近的整数差小于整周模糊度波长的0.25周；双差宽巷模糊度取整成功率超过99.9%；双差宽巷整周模糊度的中误差小于0.01周。
${{\rm{\Delta }}\nabla N_{12}^G}$、${{\rm{\Delta }}\nabla N_{i, i + 1}^B}$准确确定后，不采用宽巷组合或无电离层组合计算窄巷整周模糊度或载波相位模糊度，而是根据宽巷整周模糊度数值构建载波相位整周模糊度约束方程，添加到法方程中解算载波相位整周模糊度。GPS/BDS双差宽巷整周模糊度与载波相位整周模糊度的关系如下：
$$ \left\{ {\begin{array}{*{20}{l}} {{\rm{\Delta }}\nabla N_{12}^G = {\rm{\Delta }}\nabla N_1^G - {\rm{\Delta }}\nabla N_2^G}\\ {{\rm{\Delta }}\nabla N_{i, i + 1}^B = {\rm{\Delta }}\nabla N_i^B - {\rm{\Delta }}\nabla N_{i + 1}^B} \end{array}} \right. $$ (11) 参考站观测值的对流层延迟误差使用对流层投影函数和天顶延迟误差表示[19]。两个距离较近的参考站间同颗卫星的投影函数较接近，可将其天顶对流层延迟误差归化为相对天顶对流层延迟误差参数。因此，使用相对天顶对流层延迟误差作为待估参数，采用分段常数进行估计，双差对流层延迟误差参数表示为：
$$ \left\{ {\begin{array}{*{20}{l}} {{\rm{\Delta }}\nabla {T^G} = {\rm{\Delta Ma}}{{\rm{p}}^G} \cdot {\rm{RZTD}}}\\ {{\rm{\Delta }}\nabla {T^B} = {\rm{\Delta Ma}}{{\rm{p}}^B} \cdot {\rm{RZTD}}} \end{array}} \right. $$ (12) 式中，${{\rm{\Delta Map}}}$是星间投影函数之差，投影函数采用全球投影函数（global mapping function，GMF）；RZTD（relative zenith tropospheric delay）是参考站间相对天顶对流层延迟误差参数。
由GPS、BDS参考站双差载波相位观测方程可得：
$$ \left\{ \begin{array}{l} \lambda _1^G \cdot {\rm{\Delta }}\nabla \varPhi _1^G = {\rm{\Delta }}\nabla {\rm{Ma}}{{\rm{p}}^G} \cdot {\rm{RZTD}} - {\rm{\Delta }}\nabla I_1^G - \lambda _1^G \cdot {\rm{\Delta }}\nabla N_1^G + {\rm{\Delta }}\nabla {\rho ^G} + {\rm{\Delta }}\nabla \varepsilon _1^G\\ \lambda _2^G \cdot {\rm{\Delta }}\nabla \varPhi _2^G = {\rm{\Delta }}\nabla {\rm{Ma}}{{\rm{p}}^G} \cdot {\rm{RZTD}} - {\rm{\Delta }}\nabla I_2^G - \lambda _2^G \cdot {\rm{\Delta }}\nabla N_2^G + {\rm{\Delta }}\nabla {\rho ^G} + {\rm{\Delta }}\nabla \varepsilon _2^G\\ \lambda _i^B \cdot {\rm{\Delta }}\nabla \varPhi _i^B = {\rm{\Delta }}\nabla {\rm{Ma}}{{\rm{p}}^B} \cdot {\rm{RZTD}} - {\rm{\Delta }}\nabla I_i^B - \lambda _i^B \cdot {\rm{\Delta }}\nabla N_i^B + {\rm{\Delta }}\nabla {\rho ^B} + {\rm{\Delta }}\nabla \varepsilon _i^B\\ \lambda _{i + 1}^B \cdot {\rm{\Delta }}\nabla \varPhi _{i + 1}^B = {\rm{\Delta }}\nabla {\rm{Ma}}{{\rm{p}}^B} \cdot {\rm{RZTD}} - {\rm{\Delta }}\nabla I_{i + 1}^B - \lambda _{i + 1}^B \cdot {\rm{\Delta }}\nabla N_{i + 1}^B + {\rm{\Delta }}\nabla {\rho ^B} + {\rm{\Delta }}\nabla \varepsilon _{i + 1}^B \end{array} \right. $$ (13) 第k历元参考站GPS、BDS求解载波相位模糊度和其他未知参数的方程组为：
$$ L = \mathit{\boldsymbol{{A}}} \cdot Y $$ (14) 式中，A表示待估参数系数矩阵；L表示常数项；Y表示待估参数。
利用多历元观测数据对式（14）中的未知参数进行估计。电离层延迟误差随着卫星穿刺点的变化和时间的推移而变化，因此采用参数消去法消去电离层延迟误差参数，主要对整周模糊度参数和对流层延迟误差参数进行估计。首先使用最小二乘法估计双差载波相位整周模糊度浮点解，然后利用模糊度浮点解和方差-协方差矩阵解算双差载波相位整周模糊度。当参考站数量大于两个时，参考站间双差整周模糊度的数值和为零，可利用该闭合条件检验搜索出的双差载波相位整周模糊度结果，以提高载波相位整周模糊度确定的准确性。
2 长距离参考站双差误差分类内插计算
GPS/BDS长距离网络RTK方法实现厘米级定位的关键是通过参考站网计算出的误差改正数削弱流动站观测误差残差的影响，使残差不影响流动站载波相位整周模糊度的准确确定和高精度定位的实现。通常网络RTK采用各类观测误差综合处理的方式进行参考站网误差的改正。在长距离的情况下，使用综合误差处理的方式很难精确求出流动站的误差改正数，需要根据误差特性的不同，对双差观测误差进行分类计算。区域观测误差分类内插方法的实现过程如图 1所示。
图 1中，A、B、C为3个参考站，u为流动站用户，平面1、2分别是参考站和流动站所在的平面和中心电离层平面，在各自平面上内插计算流动站双差非色散性误差和双差色散性误差。误差内插面2上，点Ip、Iq为卫星电离层穿刺点，下标p、q为卫星号。图 2为双差非色散性误差在误差内插面1上的测站分布，其中X、Y表示参考站和流动站的平面坐标，A、B、C、u表示测站。
在确定参考站网双差载波相位整周模糊度之后，可以计算出参考站AB、AC的双差非色散性误差${{\rm{\Delta }}\nabla m_{AB}^S}$、${{\rm{\Delta }}\nabla m_{AC}^S}$，则流动站的双差非色散性误差${\rm{\Delta }}\nabla m_{Au}^S$计算值为：
$$ \begin{array}{l} {\rm{}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{\Delta }}\nabla m_{Au}^S = \left( {\begin{array}{*{20}{c}} {{X_u} - {X_A}}&{{Y_u} - {Y_A}} \end{array}} \right) \cdot \\ {\left( {\begin{array}{*{20}{c}} {{X_B} - {X_A}}&{{Y_B} - {Y_A}}\\ {{X_C} - {X_A}}&{{Y_C} - {Y_A}} \end{array}} \right)^{ - 1}} \cdot \left( {\begin{array}{*{20}{c}} {{\rm{\Delta }}\nabla m_{AB}^S}\\ {{\rm{\Delta }}\nabla m_{AC}^S} \end{array}} \right) \end{array} $$ (15) 双差色散性误差的内插计算模型与式（15）一致，不同的是色散性误差内插公式的系数是由参考站、流动站的共视卫星电离层穿刺点处的平面坐标计算得到的。当参考站间载波相位整周模糊度确定之后，计算参考站网的双差电离层延迟误差，并在中心电离层上进行流动站双差电离层延迟误差的内插计算。
3 流动站整周模糊度解算及位置确定
流动站GPS/BDS载波相位和伪距双差观测方程分别为：
$$ \begin{array}{*{20}{l}} {{\lambda _i}^S \cdot {\rm{\Delta }}\nabla {\varPhi _i}^S = {\rm{\Delta }}{\mathit{\boldsymbol{H}}^S} \cdot {\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}} + {\rm{\Delta }}\nabla {\rho ^S} - {\lambda _i}^S \cdot {\rm{\Delta }}\nabla {N_i}^S - }\\ \ \ \ \ \ \ \ \ \ \ {{\rm{}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{\Delta }}\nabla {I_i}^S + {\rm{\Delta }}\nabla {T^S} + {\rm{\Delta }}\nabla {\varepsilon _i}^S} \end{array} $$ (16) $$ \begin{array}{*{20}{l}} {{\rm{\Delta }}\nabla {P_i}^S = {\rm{\Delta }}{\mathit{\boldsymbol{H}}^S} \cdot {\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}} + {\rm{\Delta }}\nabla {\rho ^S} + {\rm{\Delta }}\nabla {I_i}^S + }\\ \ \ \ \ \ \ \ \ \ \ {{\rm{}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{\Delta }}\nabla {T^S} + {\rm{\Delta }}\nabla {\varepsilon _i}_P^S} \end{array} $$ (17) 式中，${{\rm{\Delta }}{\mathit{\boldsymbol{H}}^S}}$表示流动站坐标参数${{\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}}}$的系数矩阵。利用长距离参考站网计算流动站GPS/BDS观测值的误差改正数，可以恢复模糊度参数的整数特性，有利于实现流动站的厘米级定位，但需准确确定GPS/BDS多频载波相位的整周模糊度。各系统宽巷和L1、B1载波相位整周模糊度的解算同时进行，流动站坐标参数和整周模糊度参数解算的观测方程矩阵形式为：
$$ \left\{ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{l}} {{\rm{\Delta }}\nabla \mathit{\boldsymbol{L}}_1^G = {\rm{\Delta }}{\mathit{\boldsymbol{H}}^G} \cdot {\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}} + \mathit{\boldsymbol{H}}_N^G \cdot {\rm{\Delta }}\nabla \mathit{\boldsymbol{N}}_1^G}\\ {{\rm{\Delta }}\nabla \mathit{\boldsymbol{L}}_1^B = {\rm{\Delta }}{\mathit{\boldsymbol{H}}^B} \cdot {\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}} + \mathit{\boldsymbol{H}}_N^B \cdot {\rm{\Delta }}\nabla \mathit{\boldsymbol{N}}_1^B} \end{array}}\\ {\begin{array}{*{20}{l}} {{\rm{\Delta }}\nabla \mathit{\boldsymbol{L}}_W^G = {\rm{\Delta }}{\mathit{\boldsymbol{H}}^G} \cdot {\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}} + \mathit{\boldsymbol{H}}_W^G \cdot {\rm{\Delta }}\nabla \mathit{\boldsymbol{N}}_W^G}\\ {{\rm{\Delta }}\nabla \mathit{\boldsymbol{L}}_W^B = {\rm{\Delta }}{\mathit{\boldsymbol{H}}^B} \cdot {\rm{ \mathsf{ δ} }}\mathit{\boldsymbol{X}} + \mathit{\boldsymbol{H}}_W^B \cdot {\rm{\Delta }}\nabla \mathit{\boldsymbol{N}}_W^B} \end{array}} \end{array}} \right. $$ (18) 式中，${{\rm{\Delta }}\nabla \mathit{\boldsymbol{L}}}$表示经过参考站网观测误差改正数改正后的包括几何距离计算值的已知向量；${{\rm{\Delta }}\mathit{\boldsymbol{H}}}$表示坐标参数系数矩阵；${{\rm{\Delta }}\nabla \mathit{\boldsymbol{N}}}$表示整周模糊度估计参数向量；${{\mathit{\boldsymbol{H}}_N}}$、${{\mathit{\boldsymbol{H}}_W}}$表示相应${{\rm{\Delta }}\nabla \mathit{\boldsymbol{N}}}$的系数矩阵。下标W表示与宽巷模糊度有关，1表示与各系统第一频率有关。如果进行流动站单系统定位，只需将另一个系统的位置参数和模糊度参数的系数矩阵设为0。采用最小二乘法解算式（18），宽巷整周模糊度按照系统不同，用最小二乘模糊度降相关平差法（least-square ambiguity decorrelation adjustment，LAMBDA）分别搜索确定[20]。利用固定的宽巷整周模糊度约束载波相位整周模糊度参数的解算，添加到载波相位整周模糊度解算的法方程中，加速载波相位整周模糊度的搜索，进而计算流动站的位置坐标参数。
长距离GPS/BDS网络RTK的实现流程如图 3所示。首先利用参考站网的GPS/BDS多频观测数据进行宽巷整周模糊度的计算，然后通过顾及大气误差的载波相位整周模糊度解算模型解算双差载波相位整周模糊度。宽巷整周模糊度确定之后，将载波相位整周模糊度间的整数关系添加到参数估计的法方程中。将观测误差进行分类处理，分别计算流动站的电离层延迟误差和非色散性误差，并改正流动站的GPS/BDS观测误差，使用流动站多频载波相位整周模糊度解算方法确定GPS/BDS单系统、双系统的载波相位整周模糊度和位置参数。
![]()
图 3 长距离GPS/BDS网络RTK流程图Figure 3. Flowchart of Network RTK Between Long Range with GPS/BDSGPS/BDS双系统融合定位需要统一的时空基准[17]，BDS和GPS的时间系统都是原子时，无闰秒，且秒长定义相同，不同的是系统时间起算点，除相差1 356周外，还保持14 s的系统差，可根据这两项差异进行时间变换。国家大地坐标系（China geodetic coordinate system 2000，CGCS2000）和WGS-84的坐标系统构成相同，只是4个参考椭球参数的扁率有微小差异。为了实现两个系统之间的坐标基准统一，利用参考站网中各参考站CGCS和WGS-84坐标系统下的已知坐标解算BDS和GPS的坐标转换参数。如果采用统一定轨后的双系统预报星历，那么参考站和流动站都不需要进行坐标基准的统一。
4 实验与分析
采用华中地区3个参考站A、B、C和1个流动站U的GPS/BDS实测数据，采集时间为2018-03-31 24 h的观测数据，采样率为1 Hz。参考站间基线AB、BC、CA的长度分别是132 km、154 km、127 km，测站的地理分布如图 4所示。
首先计算参考站网GPS和BDS的宽巷整周模糊度，以BDS的C14卫星和GPS的PRN10卫星（以G10表示）第一、二频率宽巷模糊度为例，其计算结果分别如图 5、图 6所示。从图 5、图 6的计算结果可知，利用宽巷整周模糊度计算方法获得的模糊度计算值能够快速收敛，并能准确确定宽巷整周模糊度。以参考站间基线的双差模糊度数值代数和为0作为闭合条件，对宽巷整周模糊度固定结果进行检验，进一步提高模糊度固定的可靠性。M-W组合观测值能够有效消除长距离参考站间观测误差（观测噪声除外）的影响，对于长距离参考站间GPS/BDS宽巷整周模糊度的解算具有很强的适用性。
在宽巷整周模糊度准确确定的情况下，使用引入大气误差参数的载波相位模糊度解算模型可以确定载波相位的整周模糊度，并进行参考站网观测误差的计算，图 7~10分别为参考站卫星C14、G10的双差电离层延迟误差和非色散性误差的计算结果。由于C14、G10为中圆轨道卫星，在24 h的观测时长中仅有部分时段可见，只有可见时段可以观测到数据，其余不可见时段未观测到数据，因此图 7~10中的数据空白部分表示此阶段卫星不可见。由图 7~10可见，各误差相邻历元间的变化较小（除基准卫星变换引起的数值抖动外），历元之间以对流层延迟为主的非色散性误差变化量符合正常情况下的对流层变化规律，说明载波相位整周模糊度的解算正确；双差误差的数值随着双差卫星对的改变而变化，基准卫星变换时不影响流动站的误差改正精度。在参考站网载波相位整周模糊度准确固定的前提下，流动站的误差改正精度不受到测站基准卫星变换的影响，测站间距离和区域观测误差改正模型是影响观测误差改正效果的主要因素。分别选取各系统当前历元中的高度角卫星作为基准卫星，当相邻历元的基准卫星发生变换时，对之前历元的宽巷模糊度计算信息或已固定载波相位整周模糊度的GPS、BDS双差卫星进行基准卫星的变换，并利用当前弧段的观测值继续进行双差整周模糊度的解算。
由于多数历元的大气误差影响大于载波相位半个波长，因此电离层延迟误差和非色散性误差既是影响长距离参考站网载波相位整周模糊度准确解算的主要因素，也是参考站网利用误差改正模型计算流动站观测误差改正信息和解算流动站模糊度时必须考虑的因素。利用§2的误差分类内插方法，计算流动站的电离层延迟误差和非色散性误差改正数，并对流动站观测值进行改正。使用导航定位数据分析软件（position and navigation data analysist，PANDA）和测站的已知坐标，准确计算出流动站GPS/BDS两系统L1、B1载波相位的观测误差，对内插出的电离层延迟、非色散性误差改正数总和进行对比。利用综合误差内插法和误差分类内插法得到的卫星C14、G10的误差改正数残差分别如图 11、图 12所示。
![]()
图 11 两种方法的C14卫星误差改正数残差Figure 11. Residuals of Error Correction of C14 Satellite by Two Methods![]()
图 12 两种方法的G10卫星误差改正数残差Figure 12. Residuals of Error Correction of G10 Satellite by Two Methods对于长距离误差计算和改正，使用分类误差内插计算方法单独对中心电离层高度处的电离层延迟误差进行内插计算更符合实际情况。从图 11、12的误差改正数残差图可知，虽然两种方法均可获得流动站厘米级精度的误差改正信息，并能恢复流动站模糊度的整周特性，有利于载波相位整周模糊度的准确解算，但误差分类内插法的精度高于综合误差内插法。对流动站载波相位整周模糊度解算的几种策略进行了实验对比，初始化时间对比结果如表 1所示。
表 1 不同解算策略的整周模糊度初始化时间/sTable 1. Initialization Time of Integer Ambiguity with Different Solution Strategies/s解算策略 定位系统 GPS/BDS GPS BDS 先解算宽巷整周模糊度，再解算载波相位整周模糊度 16 27 68 直接估计L1载波相位整周模糊度 — 116 — 直接估计B1载波相位整周模糊度 — — 201 直接估计L1/B1载波相位整周模糊度 118 — — 从表 1可知，与直接估计流动站L1、B1双差整周模糊度相比，先解算宽巷整周模糊度，再解算载波相位整周模糊度的策略在解算速度上具有明显优势，这是因为多频组合可形成更长波长的虚拟观测值，利用该特性可加快整周模糊度的解算速度，所以建议选择先估计GPS/BDS宽巷整周模糊度、再解算L1、B1双差载波相位整周模糊度的解算策略。
在载波相位整周模糊度解算正确的情况下，流动站观测误差经参考站网误差改正数改正之后，其位置参数解算的精度直观反映了参考站网观测误差的改正数是否可靠。使用GPS/BDS两系统L1/B1载波相位观测值解算坐标参数，计算出流动站BDS、GPS单系统、GPS/BDS双系统的定位结果与已知坐标在东（east，E）、北（north，N）、高程（up，U）方向的差值，结果分别如图 13~15所示。
从图 13~15的定位结果差值可以看出，长距离GPS/BDS网络RTK方法解算出的流动站坐标与流动站已知的准确坐标在E、N、U 3个坐标分量上的差值均为厘米级。BDS、GPS、GPS/BDS定位结果差值的均方根（root mean square，RMS）如表 2所示。
表 2 定位结果差值的RMS/mTable 2. RMS of Positioning Errors/m定位系统 E分量 N分量 U分量 BDS 0.010 0.014 0.035 GPS 0.013 0.009 0.028 GPS/BDS 0.011 0.010 0.030 由上述实验结果可以看出，长距离参考站网GPS/BDS双系统网络RTK方法可以有效消除流动站GPS、BDS的观测误差。3种模式下流动站的定位精度都达到了厘米级，GPS的定位精度整体优于BDS的定位精度。且从表 1中可以看出，双系统融合的模糊度解算时间相对于单系统有所缩短，可见GPS/BDS双系统融合有利于加快流动站模糊度的解算速度。
5 结语
本文研究了一种基于长距离参考站网的GPS/BDS双系统网络RTK方法，并使用长距离CORS网观测数据对该方法进行了实验验证和结果分析。首先采用长距离参考站网GPS/BDS多频伪距和载波相位观测数据计算和确定双差宽巷整周模糊度，然后利用包含大气延迟误差参数的模糊度解算模型准确解算出GPS/BDS双系统的多频载波相位整周模糊度。使用参考站网观测误差分类内插法，按照观测误差特性不同分别计算出流动站的色散性误差和以对流层延迟误差为主的非色散性误差。分类误差内插法更符合长距离参考站网误差改正的特点，能够恢复流动站GPS/BDS载波相位模糊度的整数特性，进而完成GPS/BDS双系统载波相位整周模糊度的准确解算，最终实现流动站位置的高精度解算。实验结果表明，长距离GPS/BDS双系统网络RTK定位方法可以实现GPS/BDS单系统、双系统的厘米级定位。本文在GPS/BDS双系统融合定位时，对GPS、BDS卫星都是按照高度角进行权值的确定，未考虑系统性能和各卫星观测数据质量等方面的差异，下一步将着重研究双系统网络RTK定位中各类卫星权值的优化分配算法。
-
[1] Featherstone W E,Sproule D M.Fitting AUS-Geio98to the Australian Height Datum Using GPSleveling and Least Squares Collocation:Applicationof a Cross-Validation Technique[J].Surveying Re-view,2006,38(301):573-582[2] Shi Yimin.Modern Geodetic and Control Survey[M].Beijing:Publishing House of Surveying andMapping,2008(施一民.现代大地控制测量[M].北京:测绘出版社,2008)[3] Wei Ziqing.Re-approaching the Geodetic Coordi-nate System[J].Geomatics and Information Sci-ence of Wuhan University,2011,36(8):883-886(魏子卿.大地坐标系新探[J].武汉大学学报·信息科学版,2011,36(8):883-886)[4] Yang Aiming,Yan Jianguo,Yang Chenghong,etal.Study on Construction Surveying Control Net-work of South-to-North Water Diversion MiddleRoute Project[J].Yangtze River,2010,41(19):30-33(杨爱明,严建国,杨成宏,等.南水北调中线工程施 工 测 量 控 制 网 系 统 研 究 [J].人 民 长 江,2010,41(19):30-33)[5] Li Shi’an,Liu Jingnan,Shi Chuang.Study of theEllipsoid Transformation on the Establishment of1050 第39卷第9期 张 辛等:滇中引水工程独立坐标系统建立的关键技术研究Local Independence Coordinate System Using GPSTechnique[J].Geomatics and Information Scienceof Wuhan University,2005,30(10):888-891(李世安,刘经南,施闯.应用GPS建立区域独立坐标系中椭球变换的研究[J].武汉大学学报·信息科学版,2005,30(10):888-891)[6] Wang Lei,Guo Jiming,Shen Lili,et al.EllipsoidExpansion Method Considering Ellipsoid ParallelMisalignment and Its Application to ProjectionPlane Transformation[J].Geomatics and Informa-tion Science of Wuhan University,2013,38(6):725-733(王磊,郭际明,申丽丽,等.顾及椭球面不平行的椭球膨胀法在高程投影面变换中的应用[J].武汉 大学学报 · 信 息 科 学 版,2013,38(6):725-733)[7] Fan Yizhong,Wang Jigang,Zhao Lihua.On theOptimal Selection of the Compensating ProjectingPlane[J].Bulletin of Surveying,2000,(2):20-21(范一中,王继刚,赵丽华.抵偿投影面的最佳选取问题[J].测绘通报,2000,(2):20-21)[8] Wang Huainian.Theoritical Derivation of OptimalCompensation Projecting Plane[J].Bulletin ofSurveying,2004(10):18-19(王怀念.最佳抵偿投影面的理论推导[J].测绘通报,2004(10):18-19)[9] Sanchez L,Brunini C.Achievements and Challen-ges of SIRGAS,Geodetic Reference Frames[C].International Association of Geodesy Symposia,Munich,Germany,2006[10]Xu Shaoquan,Cheng Wenming,Huang Xuebin,etal.The Investigation of the Landslides Monitoringin the Three Gorges Reservoir Region by ApplyingGPS[J].Journal of Hydraulic Engineering,2003,(1):114-118(徐 绍 铨,程 温 鸣,黄 学 斌,等.GPS用于三峡库区滑坡监测的研究[J].水利学报,2003,(1):114-118)[11]Yang Yuanxi.Chinese Geodetic Coordinate System2000[J].Chinese Science Bulletin,2009,54:2714-2721 -
期刊类型引用(12)
1. 谢发荣，邓钰祺，胡健民，孙伟强，易堃. 基于CRPA技术的北斗/GPS对时系统误差校正. 电子设计工程. 2025(04): 146-150 . 
百度学术
2. 罗春林，廖伟，沙二，蒋军海，杨理云. 建设项目使用林地国家级保护植物调查工作思考. 林业建设. 2025(01): 22-25 . 百度学术
3. 祝会忠，郑建行，路阳阳. 一种多频多系统部分模糊度固定算法. 测绘科学. 2024(01): 1-14 . 百度学术
4. 李兆哲. 四系统组合的RTK定位对不同长度基线性能的评估研究. 测绘标准化. 2024(03): 65-70 . 百度学术
5. 胡主宽，范斌斌，周颖熙，张勇，柯维. BDS-3/BDS-2/GPS短基线RTK性能对比分析. 测绘地理信息. 2023(03): 14-19 . 百度学术
6. 祝会忠，张杰，赵齐乐，徐爱功. BDS-3/BDS-2多频数据的长距离非差网络RTK算法. 武汉大学学报(信息科学版). 2023(07): 1146-1159 . 百度学术
7. 彭长虹，张君兰，胡红，田富强，熊友，王莲冀. 植保无人驾驶航空器购置补贴分类分档与补贴额研究. 中国农机化学报. 2023(06): 251-256 . 百度学术
8. 张雨露，范磊，刘江梅，李桢. 商业低轨通信星座纳入国家综合PNT体系的可行性分析. 导航定位学报. 2022(02): 26-36 . 百度学术
9. 杨俊山，项鑫，熊晓峰，王艳利，李小强. 矿区地表形变PPP监测自适应滤波算法及其应用. 测绘通报. 2022(09): 129-133 . 百度学术
10. 陈星宇，袁运斌. 一种直接估计大气延迟的网络RTK算法研究. 数据与计算发展前沿. 2022(04): 45-54 . 百度学术
11. 杨俊山，刘素洁，杜春苗. 基于提升小波变换的矿区地表形变监测数据处理研究. 河南科学. 2022(09): 1472-1477 . 百度学术
12. 舒宝，何元浩，王利，周星，张勤，黄观文. 一种适用于大尺度卫星导航定位基准站的网络RTK方法. 武汉大学学报(信息科学版). 2021(11): 1609-1619 . 百度学术
其他类型引用(5)
下载:
计量
- 文章访问数: 1888
- HTML全文浏览量: 39
- PDF下载量: 584
- 被引次数: 17
