0 引言

北斗卫星导航系统(BeiDou navigation satellite system，BDS)是中国正在组建、自主研发、独立运行的卫星导航系统。2012-12-28，我国第二代BDS已具备亚太地区的定位、导航和授时能力，在轨卫星由5颗地球静止轨道(geostationary Earth orbit，GEO)卫星、5颗倾斜地球同步轨道(inclined geosynchronous satellite orbit，IGSO)卫星和4颗中圆地球轨道(medium Earth orbit，MEO)卫星构成。BDS计划于2020年左右建成为覆盖全球的卫星导航系统[1-4]。

相较于单一导航系统，多系统可显著改善观测冗余度，提高定位的可靠性、精度及可用性[5]。虽然兼容多系统观测数据的多系统组合试验观测网(muti-GNSS experiment，MGEX)已基本实现全球分布，但尚处于组网建设中的BDS在全球大部分区域尚不能稳定观测到4颗以上卫星。在BDS及全球定位系统(global positioning system，GPS)联合观测条件下，2个系统的观测信息间具有关联性，利用多系统观测信息，可提高尚处于组网建设中的BDS的轨道解算精度。BDS/GPS融合精密定轨可分为“一步法”和“两步法”，其中“两步法”首先利用GPS精密星历进行精密单点定位，得到测站坐标和对流层参数，将其作为已知信息再计算BDS卫星轨道，而“一步法”在同一平差模型中解算BDS和GPS轨道和地球物理参数。BDS/GPS融合定轨的基本思想是利用双系统的观测数据提高系统间的对流层延迟、接收机钟差、测站坐标等公共参数的求解精度[6]。

精密轨道作为精密定位的重要起算数据，其质量的好坏直接决定了定位结果的精度，精确的卫星轨道信息是评价导航卫星系统可用性的重要指标之一[7]。研究BDS精密定轨，提供高精度的BDS卫星轨道信息，对于推广BDS具有重要意义。文献[8]运用一种基于模糊度固定的多系统融合定轨方法，解算的BDS GEO、IGSO、MEO卫星3维精度分别达到1.236、0.214、0.134 m。文献[9]就BDS的服务性能进行了估计，并对BDS精密定轨技术进行了探讨，并利用区域站解算BDS单天轨道，3维精度达到1～10 m。文献[10]利用BDS卫星观测实验网(BeiDou experimental tracking stations，BETS)的实测数据解算的BDS精密轨道径向精度优于10 cm。文献[11]利用多系统观测数据实现了多模融合定轨。

本文在已有研究的基础上采用“一步法”进行BDS/GPS融合精密定轨，在同一平差模型中同时估计GPS、BDS卫星的轨道参数，并将精密定轨的结果与武汉大学分析中心提供的相应产品进行比较，统计、分析BDS/GPS融合精密定轨的精度。

2.5 两组护理满意度比较 观察组护理满意度为(94.33±10.60)分，对照组(89.46±10.19)分，两组比较差异有统计学意义(t=2.068，P<0.05)。

1 BDS/GPS融合精密定轨数学模型

1.1 观测方程

式中：P为伪距观测值；L为载波相位观测值；ρ为站星间几何距离；dts和dtk分别为卫星和接收机钟差；c为真空中光速；δtg为BDS和GPS系统间信号延迟；dρtron为对流层延迟；dρion为电离层延迟；pr和ps分别为伪距观测值在接收机和卫星端的硬件延迟；lr和ls分别为载波相位观测值在接收机和卫星端的硬件延迟；N为整周模糊度；λ为相应观测值的波长；εP和εL分别为伪距和载波的多路径、观测噪声等未模型化误差[11]。

(1)

不同于单一的导航卫星系统，由于BDS和GPS信号调制方式、频率的差别，2个系统的信号在接收机端会存在系统间信号延迟；故在BDS/GPS融合精密定位和定轨的观测方程中需要顾及BDS和GPS系统间信号延迟量δtg，在实际数据处理过程中以GPS卫星信号解算的接收机钟差为参考，估计BDS相对于GPS的系统间信号延迟。伪距和载波相位观测值的非差观测方程为

进行标贯试验85次，实测锤击数N’=1.0～2.0击，平均为1.1击；经杆长修正后N=0.7～1.7击，平均为0.9击。取样13组，共13件进行统计分析。详见附表4和土工试验报告附表。建议其承载力特征值的经验值为fak=60kPa。

BDS采用由GEO、IGSO和MEO卫星组成的混合星座，其卫星定轨精度与卫星轨道类型有关；因此按照卫星轨道类型对轨道和钟差结果进行统计。图3分别为BDS卫星在卫星轨道坐标系A、C、R方向轨道偏差的RMS统计结果；BDS GEO、IGSO、MEO卫星A、C、R方向轨道差值的平均RMS统计值如表3所示。通过图3和表3可以看出BDS不同轨道类型的卫星轨道精度存在着较明显的差异。

(2)

(3)

式中：PC和LC分别为无电离层组合的伪距和相位观测值；f为观测值的频率[12]。

1.2 动力学原理及数据处理策略

导航卫星的在轨运行是一个受摄运动的过程，动力学定轨的基本思想是依据卫星的受力建立接近卫星真实运动状态的运动方程。卫星的运动方程由一组待求的卫星初始轨道参数表示。几何观测信息用于解算观测时刻的卫星坐标，可通过积分变分方程计算该时刻卫星坐标，利用初始卫星初始轨道参数和变分方程可推导出任意时刻的卫星位置。

卫星运动方程可表达为

(4)

式中：为卫星所受的保守力引起的加速度，包括地球中心引力、地球非球形引力摄动加速度、N体引力摄动加速度、潮汐效应引起的摄动加速度、广义相对论摄动力加速度；为卫星所受非保守力，包括太阳辐射、地球反照辐射和热辐射摄动力加速度、大气阻力摄动加速度；为经验加速度，其不具有实际物理意义，由摄动力模型误差和未模型化的摄动力引起[11]。

BDS/GPS融合精密定轨的动力学模型及数据处理策略与单系统精密定轨模型类似，主要差异在于BDS和GPS观测量精度不同及2系统间信号延迟；同时不同轨道卫星的姿态控制模型不同，GPS卫星姿态控制模型为动偏，BDS GEO卫星姿态控制模型为零偏，当BDS IGSO和MEO卫星轨道角为90o且轨道夹角小于4o时卫星姿态控制模型由动偏转换为零偏[12]。动力学模型如表1所示。数据处理策略如表2所示，其中BDS的IGSO/MEO卫星采用欧洲航天局(European Space Agency，ESA)提供的模型进行卫星天线相位中心偏移(phase center offsets，PCO)和卫星天线相位中心变化(phase center variations，PCV)改正。

表1 BDS/GPS融合精密定轨动力学模型策略

模型 方法模型地球重力场EGM2008模型前8阶潮汐固体潮、海洋潮、极潮N体引力太阳、月球、金星、水星、火星、木星、土星、天王星、海王星、冥王星天体星历jpleph_de405太阳光压BREN5参数地球反照辐射不考虑相对论效应考虑 大气阻力不考虑轨道积分Runge-Kutta单步法,Adams预报校正

表2 BDS/GPS融合精密定轨数据处理策略

参数 方法模型观测量非差无电离层组合LC/PC观测值定权模型高度角定权卫星截止高度角7°卫星天线相位中心改正GPS:igs08.atx改正PCO和PCVBDS:GEO改正PCOIGSO/MEO采用ESA的PCO/PCV改正模型接收机天线相位中心改正GPS:igs08.atx改正BDS:不改正天线相位缠绕改正卫星姿态GPS:动偏BDS:GEO零偏IGSO/MEO动偏、零偏切换对流层延迟Saastamonine模型估计+随机游走电离层延迟双频无电离层组合消除一阶项,高阶忽略潮汐改正固体潮、海洋潮、极潮相对论效应考虑卫星钟差过程估计+白噪声接收机钟差过程估计+白噪声系统间延迟过程估计+白噪声模糊度浮点解初始轨道参数X,Y,Z,VX,VY,VZ,Bren光压参数(D0,Y0,B0,Bc,Bs)地球自转参数估计

2 实验及结果分析

GPS卫星在A、C、R方向偏差的平均RMS统计值如表3所示，通过统计结果可以看出：BDS/GPS融合精密定轨解算出的GPS轨道信息与武汉大学分析中心提供的精密轨道信息在各个方向表现出较好的一致性，卫星各方向的轨道偏差约为1～4 cm；其中R方向的轨道精度最高，偏差约为0～2 cm，R方向卫星的平均RMS值为0.99 cm，C方向卫星轨道的偏差约为1～2 cm，精度为1.33 cm，A方向的轨道偏差相对较大约为1～4 cm。钟差精度以钟差标准差(standard deviation，STD)表示。

图1 跟踪站分布

图2分别给出GPS卫星在卫星轨道坐标系的切向(along-track，A)、法向(cross-track，C)和径向(radial，R)3个方向偏差的均方根(root mean square，RMS)统计结果。

实验采用2017-05-20—2017-05-24(年积日第140天至第144天)共5 d的观测数据，选取全球约110个跟踪站，其中约60个测站可观测到BDS/GPS数据，跟踪站分布如图1所示。采用表2的数据处理策略，取1 d为一个定轨弧长，以武汉大学分析中心提供的事后精密轨道和精密钟差产品为参考，计算BDS/GPS融合定轨解算的轨道和钟差与武汉分析中心产品的偏差，并对结果进行统计、分析。

表3 BDS/GPS融合定轨轨道RMS和钟差STD统计结果

系统卫星A方向RMS/cmC方向RMS/cmR方向RMS/cm钟差STD/nsGPS1.901.330.99 0.07BDSGEO171.9219.5112.630.41IGSO24.9617.7612.520.30MEO19.3311.027.480.18

粒煤是当前安顺煤矿最小的块煤品种，安顺煤矿于2011年安装振动筛筛分粒煤，筛板为长条筛，条间距1cm，所筛粒煤所占商品煤比例在15%左右，每年有10万吨左右的产量。然而，振动筛有处理能力大的优点，也存在筛分质量差的缺陷，长条筛筛条易变形或破损，导致部分1-2cm颗粒进入粉煤中。通过粉煤筛分试验，粉煤中1cm以上颗粒约占15%，0.6cm-1cm粒度约占30%，有很大的挖潜空间。

图2 GPS轨道精度

患者术后咳嗽、喷嚏、用力排便以及坐起站立的时机过早可使腹压骤然升高，硬膜囊内压力瞬间增大，变薄的硬脊膜撕裂或硬脊膜裂重新裂开，造成脑脊液漏［7］。

图3 BDS轨道精度

导航卫星观测值是卫星与测站之间的相对时间延迟，数据解算时不能同时确定所有卫星和接收机钟差，在数据结解算时首先需要固定某一基准钟的钟差，然后计算其余导航卫星和地面接收机的相对钟差；虽然计算出的卫星钟差结果不是各卫星的绝对钟差，但只要基准钟的精度优于1×10-6 s，相对卫星钟差对定位结果是等价的，在定位模型中相对钟差的系统性偏差的平均项会被接收机钟差完全接收，而不影响位置参数的解算精度。为更加客观地评价解算的卫星钟差的精度，采用二次差比较的方法对卫星钟差进行精度评价，即在一次差的基础上选择一颗参考卫星，将所有的卫星钟差与该参考星做差，消除基准钟引起的卫星钟系统误差[13-14]。

基于双频观测数据建立无电离层组合，消除一阶电离层影响，在式(1)的基础上建立BDS/GPS融合精密定轨的无电离层伪距(PC)和无电离层相位(LC)观测方程，可分别表示为

通过图3可以看出，BDS GEO卫星(C01～C05)精密轨道的精度明显次于IGSO卫星(C06～C10)和MEO卫星(C11～C14)，这主要是由于GEO卫星的静地几何特性导致的。GEO卫星轨道在A方向的精度达到了米级，偏差约为1～3 m；在R、C方向的精度为分米级，其中在C方向各卫星的平均RMS值为19.5 cm，R方向精度相对较高，偏差约为0～20 cm。IGSO卫星精密轨道在A、C、R方向的RMS统计结果全部为分米级，在A方向轨道偏差约为10～40 cm；C方向的平均RMS统计值为17.76 cm；同样R方向的精度相对较高，平均RMS统计值为12.56 cm。BDS MEO卫星定轨结果与武汉大学分析中心提供的精密轨道信息表现出较好的一致性，MEO卫星轨道精度明显高于IGSO卫星和GEO卫星，其A和C方向的精度为分米级，在R方向的偏差约为5～10 cm，A、C和R方向的RMS统计值分别为19.33、11.01和7.48 cm。

图4和图5分别为融合精密定轨解算的GPS和BDS卫星钟差选用GO1号卫星作为参考星，与武汉大学分析中心提供的精密钟差做二次差的STD统计结果；GPS卫星和BDS各星座卫星的二次差的平均STD值如表3所示。通过统计结果可以看出，融合定轨解算出的GPS卫星钟差表现出较好的稳定性，各颗卫星钟差的统计精度均优于0.15 ns，卫星的平均STD统计为0.07 ns。

图4 GPS卫星钟差精度

图5 BDS卫星钟差精度

BDS卫星钟差二次差的STD统计结果同样和卫星轨道类型有关。MEO卫星钟差的稳定性最好，各卫星的钟差均优于0.2 ns；IGSO卫星的钟差约为0.2～0.4 ns，各卫星的平均STD为0.30 ns；GEO卫星钟差误差相对较大，约为0.3～0.5 ns，其平均钟差精度达到0.41 ns。造成GEO卫星钟和轨道存在较大误差的原因是，GEO卫星与测站之间的相对位置几乎保持不变，使得各种误差不能有效地分离出来。

一是项目法人质量管理意识薄弱。工程建设前期，个别项目法人因各种原因未能及时办理质量监督手续，存在补办质监手续的现象，影响工程的顺利实施。二是部分施工单位对质监工作不够重视，不按程序签字，资料整理不规范，甚至存在事后补充资料的现象。部分工程一线施工人员多数是没有经过培训的农民工，工程质量得不到充分保证。三是监理人员水平参差不齐，部分监理人员业务知识能力欠缺，作用发挥不充分，资料审核把关不严，不能发现问题或发现问题处置不当。

3 结束语

本文根据非差无电离层组合值和动力学定轨原理，基于全球分布的跟踪网观测数据，研究实现BDS/GPS融合精密定轨，并以武汉大学分析中心提供的BDS和GPS精密轨道和钟差产品作为参考，统计、分析BDS/GPS融合精密定轨解算的精密轨道和钟差的精度，实验结果显示：

1)BDS卫星精密轨道精度和卫星轨道类型有关，MEO卫星的轨道精度相对较高，其轨道R方向的精度为厘米级。IGSO卫星在各个方向的精度均为分米级。GEO卫星由于卫星的静地性导致其轨道误差较大，A方向的轨道精度为米级，R和C方向精度为分米级。

他们俩趴在窗台上，隔着玻璃看窗外奔走的行人和烤鸭店拉下卷闸门的老板，台风已经逼近这里，头顶上黑压压的乌云越来越低，大雨轰然落下。

2)BDS的MEO卫星钟差稳定性优于IGSO和GEO卫星，统计结果优于0.18 ns。由于GEO卫星的静地性使其卫星钟差不能有效地从轨道误差中分离出来，导致GEO卫星钟差存在较大误差，各GEO卫星钟差STD的平均值为0.41 ns。IGSO卫星钟差STD的平均值为0.30 ns，精度介于MEO卫星和GEO卫星之间。

3)GPS卫星在轨道坐标系A、C、R方向的精度都较高，3维精度优于4 cm，其中R方向的精度能达到1 cm左右。卫星钟稳定性好，钟差STD达到0.07 ns。

持刀人看见已经高声向他警告了几遍不要靠近自己，可他还是继续向自己走来，心里就不断嘀咕老子耍横他狗日的耍不要命，看来是遇到比我还狠的人了。想弃刀逃跑吧，又觉得自己要挟的钱还没有到手，心有不甘。瞬间权衡了一下处境后，持刀人决定先下手为强，挥刀就朝牛黄丸扑来。

参考文献

[1] 杨元喜，李金龙，王爱兵，等.北斗区域卫星导航系统基本导航定位性能初步评估[J].中国科学：地球科学，2014，44(1)：72-81.

[2] 朱永兴，冯来平，贾小林，等.北斗区域导航系统的PPP精度分析[J].测绘学报，2015，44(4)：377-383.

[3] 张小红，左翔，李盼，等.BDS/GPS精密单点定位收敛时间与定位精度比较[J].测绘学报，2015，44(3)：250-255.

[4] ZHAO Q，WANG C，GUO J，et al.Assessment of the contribution of BeiDou GEO，IGSO，and MEO satellites to PPP in Asia-Pacific region[J].Sensors，2015，15(12)：29970-29983.

[5] 戴小蕾，施闯，楼益栋.多GNSS融合精密轨道确定与精度分析[J].测绘通报，2016(2)：12-16.

[6] 张睿，杨元喜，张勤，等.BDS/GPS联合定轨的贡献分析[J].武汉大学学报(信息科学版)，2017，42(5)：600-608.

[7] 郭斐.GPS精密单点定位质量控制与分析的相关理论和方法研究[M].武汉：武汉大学出版社，2016：46-59.

[8] 刘伟平，郝金明，李建文，等.多GNSS融合的北斗卫星精密定轨[J].测绘学报，2014，43(11)：1132-1138.

[9] MONTENBRUCK O，STEIGENBERGER P，HUGENTOBLER U，et al.Initial assessment of the COMPASS/BeiDou-2 regional navigation satellite system[J].GPS Solutions，2013，17(2)：211-222.

[10] 施闯，赵齐乐，李敏，等.北斗卫星导航系统的精密定轨与定位研究[J].中国科学：地球科学，2012，42(6)：854-861.

[11] 李敏.多模GNSS融合精密定轨理论及其应用研究[D].武汉：武汉大学，2011：17-25.

[12] 郭靖.姿态、光压和函数模型对导航卫星精密定轨影响的研究[D].武汉：武汉大学，2014：81-82.

[13] 楼益栋，施闯，周小青，等.GPS精密卫星钟差估计与分析[J].武汉大学学报(信息科学版)，2009，34(1)：88-91.

[14] 刘伟平，郝金明，于合理，等.导航卫星精密轨道与钟差确定方法研究及精度分析[J].测绘通报，2014(5)：5-7.