北斗卫星导航系统的毫米级精度变形监测算法与实现

1、北斗卫星导航系统的毫米级精度变形监测算法与实现 肖玉钢1 ,姜卫平2-,陈华1,袁鹏2,席瑞杰11. 武汉大学测绘学院，湖北武汉430079;2. 武汉大学卫星导航定位技术研究中心，湖北武汉430079收稿日期：2014-12-08;修回日期：2015-07-08基金项目：国家 863计划(2012AA12A209);国家自然科学基金(41374033 )。第一作者简介：肖玉钢(1984 )，男，博士生，研究方向为 GNSS高精度定位定轨算法。通信作者：姜卫平，摘要：研究了北斗卫星导航系统(BeiDou Navigation SatelliteSystem，BDS)毫米级精度变形监测算法。首先

2、改进了 TurboEdit方法，以能够探测到1周的小周跳；针对 BDS星座 结构给出更为高效的独立双差观测值搜索方法；对于模糊度固定，采用决策函数和序贯模糊度固定相结合的方法。在此基础上，研制了 BDS变形监测软件。最后，利用变形监测试验平台的实 测数据，从星座分布、解算精度等方面分析了BDS在变形监测中应用的可行性。结果表明，目前在试验区域内BDS与GPS在卫星几何分布等方面基本相当。BDS的短基线解算精度略低于GPS，但仍可达到平面1 mm以内、高程2 mm以内的精度水 平。关键词：北斗卫星导航系统变形监测软件实现精度分析Research and Realizatio n of Defor

3、mati on Mon itori ng Algorithm with Millimeter Level Precision Based on BeiDou Navigati on Satellite SystemXIAO Yugang 1, JIANG Weiping 一，cheN Hua 1, YUAN Peng2, XI Ruijie 1Abstract : The deformati on mon itori ng algorithm with millimeter level precisi on based on BeiDou Navigati on Satellite Syste

4、m (BDS) was researched. The TurboEdit method was improved to detect small cycle slips, e.g. 1 cycle. Focus ing on BDS con stellati on, a more efficie nt algorithm used to con struct double-differe need observati ons was developed. The Bootstrap+Decisi on function method was utilized to improve the p

5、robability of biases fixing. Based on the improved algorithm above, a deformatio n mon itori ng software based on BDS was achieved. Afterwards, the availability of BDS in the field of deformati on mon itori ng was an alyzed in terms of satellites distribution and precision and accuracy of soluti ons

6、, utilizi ng the observatio ns acquired from the experime ntal platform. The con clusi on was draw n thatcurre ntly BDS is similar to GPS in terms of satellites distribution in the test area. The precision of short baselines derived from BDS is better than 1 mm for the horiz on tal comp onen ts, bet

7、ter tha n 2 mm for the vertical components, which is still a little lower than GPS.Key words : BeiDou Navigati on SatelliteSystemdeformati on mon itori ngsoftwareachieveme nt precisi on an alysis随着各种大型结构体的大量涌现以及滑坡、泥石流等地质灾 害的频繁发生，变形监测研究的重要性日益突出，变形监测理论和技术方法也在迅速发展。全球卫星导航系统（GNSS）具有全天候、高精度等优点，早在 20世纪80年代

8、中后期，就被作为一 种变形监测技术手段。1995年以来，GNSS被用于监测滑坡、 大坝等缓慢变形，如在Pacoima、清江隔河岩等大坝建立的 GPS 变形监测系统，对GPS用于高精度变形监测的可行性进行了分 析，精度达毫米级。结果表明，GPS解算所得变形量与水库蓄水 量具有高度一致性1,2, 3, 4。为降低变形监测系统的建设成本，文献5提出了一机多天线变形监测系统，采用附加的多天线转换开 关实现一台接收机与多个天线的时分单通连接。同时，考虑到监测特点（如监测站坐标已知），文献6 8提出了“无周跳无模糊 度”的变形监测方法。对于桥梁、高大建筑物等结构体的动态变 形，适宜采用实时监测系统进行，文

9、献 9研制了 GPS变形监测 数据处理软件GPSMON，定位精度为亚厘米级；文献1°集成 GPS、加速度计、伪卫星等技术，建立了桥梁变形实时监测系统。同时，Geo+、Leica、Trimble等测量仪器公司也开发了 GPS 变形监测系统，并得到了广泛的应用。中国自20世纪80年代开始研制北斗系统，并实施了系统建 设的“三步走”规划11。2012年底，北斗卫星导航系统空间信号接口控制文件正式公布，标志着北斗区域导航定位系统的正式 建成。目前BDS在轨卫星16颗，其中工作卫星14颗，包括5 颗GEO卫星，5颗IGSO卫星及4颗MEO卫星。预计到2020 年左右，将建成覆盖全球的北斗卫星导

10、航系统。已有研究表明，北斗二代系统的卫星钟性能、测距精度等与GPS基本处于同一水准12, 13。在变形监测中，尤其对于滑坡、大坝等变形体的监测，受观 测环境所限，测站往往遮挡严重，从而影响卫星的几何构型，降 低变形监测的精度及可靠性。 BDS由于其特殊的星座设计（5颗 GEO卫星以及5颗IGSO卫星的存在），即使在目前尚未完成全 球构网的情况下，在我国大部分区域其卫星可见数也大于GPS等系统。随着BDS卫星的逐步入轨，这一优势会更加明显。虽 然BDS与GPS等类似，但其在信号结构、卫星分布、卫星类型 等方面还存在较大的不同14。同时，利用多种系统进行监测也可 以提高监测结果的精度和可靠性。此外

11、，针对BDS变形监测应用的研究也可以拓展我国北斗系统的应用空间，验证北斗系统在高精度导航定位领域的可用性。因此，研究利用BDS及联合其他系统进行变形监测具有重要的科学和现实意义。本文首先基于双差模型研究了 BDS高精度基线解算方法， 并针对变形监测特点，采用可靠的数据处理策略， 研制开发了一 套基于BDS的变形监测软件平台。同时，本文基于配备了三维 精度测试系统的变形监测试验平台，利用在中国中部某地区采集的实测数据，从卫星分布、解算结果的内、外符合精度等多方面 对北斗系统在变形监测中的可用性进行了探讨分析，并与GPS解算结果进行了比较。1系统算法与实现为探讨利用现有北斗星座进行变形监测的可行性

12、，本文研制了变形监测软件，其数据处理流程如图1所示。本软件目标为实现BDS变形监测数据的准实时（采用时段解模式）、高精度处 理，故软件采用广播星历进行解算。已有研究表明，在短基线条 件下，广播星历误差对基线解算结果的影响较小，可忽略不计15 软件同时具有BDS-GPS联合解算能力，但已有研究表明目前 BDS与GPS高精度解算结果之间存在由于天线相位中心改正等 原因引起的系统偏差16，因此本文仅针对BDS单模结果分析。图 1 数据处理流程 Fig. 1 Data processing flow chart图选项软件包含了从数据处理到结果输出整个流程，主要包括数据 准备、法方程形成与叠加、模糊度固

13、定、时间序列结果分析等4部分。(1) 数据准备主要包括两方面，一是根据广播星历计算卫星位置；二是将观测数据进行周跳探测与标记。周跳探测采用改进的TurboEdit方法17, 18，利用MW、LG组合观测值以及LC观 测值的双差残差分别探测周跳。结果表明，上述周跳探测策略基本不受电离层活动及观测条件的影响，可以探测到1周的小周跳，从而大幅减小了残差编辑的工作量，缩短了数据解算时间。为避免周跳错误修复而引起的灾难性后果，本文对周跳采取只探测、标记而不修复的策略，对所有周跳引入模糊度参数进行估计。(2) 形成双差观测方程，进行法方程叠加，并将模糊度映射为双差模糊度19。双差观测值选取采用参考站-参考

14、星以及全局 搜索相结合的方法。若某历元观测存在参考站或参考星，则采用 参考站-参考星法，否则采用全局搜索法，在当前历元所有观测 值中最大限度地搜索函数独立的双差观测值。考虑到BDS星座的特殊性，在选择参考星时首选 GEO卫星，若由于环境遮挡任 何一颗GEO卫星均不为参考星，则依次选择 IGSO、MEO卫星 为参考星。根据基线长度，解算时既可采用单频观测值，也可利用双频观测值形成无电离层组合以消除电离层误差的一阶影响。 形成法方程时需考虑的误差主要有对流层改正、地球自转改正、 天线相位中心改正20、潮汐负载改正等，其他误差认为在双差过程中得到消除。对流层改正采用Saastamoi nen 模型，

15、剩余对 流层影响采用分段线性模型估计。BDS卫星天线相位中心采用MGEX(multi-GNSS experime nt)公布的模型，接收机天线相位中心改正模型未知，在系统中以0代替。求解法方程后得到浮点 解，此时需进行残差编辑，对可能存在的周跳和坏值进行标记， 然后重新形成法方程，直到残差编辑没有周跳或坏值为止。(3) 模糊度固定采用决策函数和序贯模糊度固定相结合的方法进行19。首先计算每个模糊度可以被固定的概率，并固定对应最大概率的模糊度，再更新法方程，重复上述程序，直到模糊度 全部被固定或没有模糊度可以被固定。结果表明，在变形监测应用中，由于基线一般较短，采用此模糊度固定策略基本可将所有

16、的北斗双差模糊度固定。将上述过程得到的整数模糊度回代法方 程解算得到模糊度固定解。(4) 时间序列分析在频域与时域进行。频域分析主要研究原始变形信号中各地球物理效应及GNSS技术类误差的影响，通过设计滤波器从原始结果中提取真正的变形信息。时域分析通过对变形时间序列进行 ARIMA建模以分析变形的内在联系， 预测 变形的发展情况。2数据分析2.1数据采集为了分析利用北斗系统进行变形监测的精度和可靠性，本文搭建了模拟变形平台，利用实测数据从内、外符合精度等多个角度出发对试验结果进行了讨论。试验采用的仪器为TrimbleNetR9型接收机，天线型号为TRM29659.00。此类型接收机可 同时接收G

17、PS以及BDS信号，满足系统间的兼容与互操作要求。本次试验共布设 3个测站（JZ01、JC01、JC02），位于中国 中部某城市。所有测站均为土层观测墩，高出地面3m。其中JZ01 站基座深8m，为钢筋混凝土结构，JC01、JC02站基座深3m， 为钢结构。各测站视野开阔，10。高度角以上基本不存在遮挡物。 测站间基线长度如表 1所示。3个测站均配置有强制对中标志。 另外JC01站装有高程精度测试系统，JC02站装有水平精度测 试系统。两套系统均可以通过旋转螺栓使接收天线在水平或垂直 方向上精确移动。表1测站间基线长度Tab. 1 Len gth of baseli nes基线名称JZ01-

18、JC01JZ01-JC02JC01-JC02长度/m29.2277.1274.5表选项本次试验采集了从2014年8月7日至8月10日（年积日219 222）共4d的数据。采样间隔30s，截止高度角10。数 据采集过程中试验平台的位移量如表 2所示。其中JC01、JC02 站分别在垂直、水平方向（NS方向）上移动。在此过程中JZ01站 保持不动。表2试验平台位移量 Tab. 2 Displacements ofIexperime ntal platformmmDOYNEU2190 0 02201 0 12212 0 2222303表选项2.2数据质量与GPS相比，BDS具有自己的独特性，其包含

19、5颗地球静 止轨道卫星（GEO）和5颗倾斜地球同步轨道卫星（IGSO）。地球静 止轨道卫星与测站相对静止，而倾斜地球同步轨道卫星与测站的 相对位置变化也与GPS卫星等不同，因此有必要对各测站能观 测到的北斗卫星数及其 GDOP值进行分析。试验中3个测站相距较近，能观测到的卫星数也十分相似， 故此处以JZ01站为例分析能观测到的北斗卫星数。分析中选用JZ01站219日GPS时424时的观测数据。结果如图2所示， 横坐标是以历元形式表示的观测时间， 共2400历元的观测数据。 由图2可知，JZ01站在分析时段内最少可观测到 7颗北斗卫星， 最多12颗，其数量均多于 GPS（最少6颗，最多11颗）。

20、这主 要是由北斗卫星导航系统特殊的星座设计引起的。随着 BDS星座的逐步完善，其在亚太地区卫星可见性方面的优势会更加明显。图 2 JZ01 站卫星可见性 Fig. 2 Satellite visibility for site JZ01图选项几何精度因子是衡量定位精度的重要系数，代表GNSS测距误差造成的接收机与卫星间的距离矢量放大因子。图3为JZ01站BDS与GPS在分析时段内的GDOP值序列。结合图 2可看 出，与GPS相比，虽然BDS在卫星可见数方面略占优势，但其 GDOP值的平均值在观测时段内大于 GPS。与BDS相比，GPS 的GDOP值变化更频繁。GDOP值序列的断续表示系统可见星

21、 数的变化，由此引入更多的模糊度参数需要解算，增加数据处理负担。Oduuc图 3 JZ01 站 BDS 与 GPS GDOP 值对比 Fig. 3 Comparison ofGDOP at site JZ01 for GPS and BDS respectively图选项 2.3定位精度本文主要从基线解算中误差、基线各分量重复性以及变形监测试验系统调整前后基线分量较差等方面来分析和讨论BDS基线解算精度。为了对比分析，本文将GPS解算的结果作为参考值。GPS数据的解算利用DDMS 4软件进行。解算时段长度为 4h。由于每天9时（北京时间）左右调整变形监测试验系统的位移 量，因此舍弃第一个时段（

22、812时）的数据，每天只统计5个时 段的结果。下述对基线解算中误差、 重复性的讨论均基于年积日 219日，基线（JZ01-JC02）的结果进行。GPS与BDS每个测段的基线分量中误差如图4所示，其中横坐标表示不同系统（BDS与GPS）、不同时段（15时段）所对应 的解算结果。由图 4可看出，在本试验的基线长度和观测环境 条件下，BDS的基线水平分量精度0.7mm左右，高程分量精度 1.5mm左右，略低于GPS的结果精度。与GPS类似，利用BDS解算得到的高程分量精度约是水平分量的一半。 同时需要注意的 是，BDS的N方向精度明显低于E方向，这与文献21的结论一 致，这可能主要是由BDS的星座结

23、构引起的。2 10口n 口匸ulS3 H.B图 4 基线分量中误差 Fig. 4 Standard deviations of baseline comp onents图选项表3为利用5个测段的结果统计得到的各系统对应基线分量 的重复性。由表 3可得到与图4类似的结果，即BDS的E方 向精度最高，N方向次之，高程分量精度最低，约为水平分量的 一半。表3各系统基线分量重复性 Tab. 3 Repeatability ofbaseli ne comp onents derived from GPS and BDS系统N分量E分量U分量L分量GPS0.20.21.00.1BDS0.90.41.40.

24、9表选项表4为各系统所对应的平台实际调整量与解算结果的较差统 计。由表4中较差平均值可看出，解算结果存在明显的系统误 差。通过对不同系统、不同调整量所对应的较差平均值取平均后 可得此系统误差约为水平方向0.7mm ，高程方向0.5mm。后用测微器（可精确至0.01mm）的量测结果证实了上述系统误差的存 在，且两种手段所得系统误差值具有较好的一致性。推断此系统误差由试验平台自带的刻度误差引入，是由于试验平台的制造工艺产生的。因此下述对基线分量较差的分析均是根据去除系统误 差后的结果进行的。表4平台实际调整量与解算结果较差统计Tab. 4 Differeneebetwee n the platfo

25、rm adjustme nts and the estimatedresultsmm实际调最大值最小值平均值标准差整量水平垂直水平垂直水平垂直水平垂直1GPS0.91.40.1-0.40.50.40.30.7BDS1.51.70.5-0.60.80.50.40.82GPS0.91.00.6-0.50.70.30.20.6BDS0.80.40.5-1.10.60.20.10.63GPS0.91.10.6-0.20.70.60.10.6BDS1.21.70.60.20.90.70.20.6表选项图5为试验平台调整前后各测段基线分量较差。由图5可以看出，对于3mm的变形，无论此变形发生在水平方向或高

26、程方 向，基于目前星座的BDS均可轻易识别。当此变形为2mm时， 水平方向仍可轻易识别，但高程方向的较差已不十分明显。进一步，当变形量为1mm时，水平方向仍可以分离出此变形，但高 程方向的基线分量较差表现出较大的随机性，已不足以提供明确的变形信息。因此，结合本文上述对基线分量估值中误差以及重 复性的讨论，认为基于目前的 BDS以及本文所实现的变形监测 软件平台可达到水平1mm，高程2mm左右的监测精度。同时 由图5也可以看出，目前的BDS解算精度略低于GPS。结合本 文前述对BDS数据质量的分析，认为其由试验区域内BDS较大 的GDOP值引起。虽然试验中 BDS的可见星数等与GPS基本 相当，

27、但不完善的星座结构导致 BDS较大的GDOP值，从而造 成BDS的解算精度略低于 GPS。但可以预测，随着北斗系统星 座的逐步构建，这一差距会逐渐缩小， 最终接近甚至在某些区域 超过GPS精度。M n m W m m聘餌234512345GPSBDS2345123432JO1 LUUVAj鏗BDSGPS234GPS234BDS注：上、中、下图分别对应年积日220、221、222日的结果，虚线表示平台实际位移值。图5试验平台调整前后各测段基线分量较 差 Fig. 5 Difference of baseline components before and afteradjusti ng the

28、experime ntal platform图选项3结论本文主要分析了基于目前的北斗星座进行高精度变形监测的可行性。通过实测数据分析与讨论，得到如下结论：在试验区域仲国中部某城市），目前BDS卫星可见数略大于 GPS。但由于BDS星座的特殊设计，其对应的GDOP值略差于 GPS；通过分析变形监测试验平台的实测数据，认为目前 BDS在短基线 条件下能够达到平面1mm以内，高程2mm以内的精度水平， 可满足大部分变形监测工程的需要。 总之，目前的北斗卫星导航 系统能够满足高精度变形监测工程的需求，可以在生产实践中推 广应用。参考文献1 HUDNUT K W, BEHR JA. Contin uou

29、s GPS Mon itori ngof Structural Deformati on at Pacoima Dam,CaliforniaJ. Seismological Research Letters, 1998, 69：299-308.2 BEHR J A, HUDNUT K W, KING N E. Mon itori ngStructural Deformati on at Pacoima Dam, Califor niaUsi ng Contin uous GPSC/Proceedi ngs of the 11thInternational Technical Meeting o

30、f the SatelliteDivisi on of the In stitute ofNavigatio n.Nashville,Te nn essee: s. n., 1998: 59-68.3 姜卫平，文U经南.GPS技术在隔河岩大坝监测中的应用研究 J.武汉测绘科技大学学报,1998, 23(S1): 48-49. JIANG Weiping, LIU Jingnan. Study on Application of GPS inthe Geheya n Dam Deformatio n Mon itori ngJ.Jour nal of Wuha n Tech ni cal Uni

31、 versity of Survey ingand Mappi ng, 1998, 23(S1): 48-49.4 姜卫平，刘鸿飞，刘万科，等.西龙池上水库GPS变形监测 系统研究及实现J.武汉大学学报(信息科学版),2012, 37(8): 949-952. JIANG Weipi ng, LIU Hon gfei, LIUWan ke, et al. CORS Developme nt for Xilon gchi DamDeformati on Mon itori ngJ. Geomatics andIn formation Scie nee of Wuhan Un iversity,

32、2012, 37(8): 949-952.5 DING X L, CHEN Y Q, HUANG D F, et al. SlopeMon itori ng Using GPS: A Multi-a ntenna ApproachJ.GPS World, 2000, 11(3): 52-55.李征航，张小红，朱智勤.利用GPS进行高精度变形监测的新模型J.测绘学报,2002, 31(3): 206-210. LIZhe ngha ng, ZHANG Xiaoho ng, ZHU Zhiqi n. A NewModel of High Accuracy Deformatio n Mon itor

33、i ngwith GPSJ.Acta Geodaetica et Cartographica Sin ica, 2002, 31(3): 206-210.7 王新洲，花向红，邱蕾.GPS变形监测中整周模糊度解算的新方法J.武汉大学学报(信息科学版),2007, 32(1):24-26. WANG Xin zhou, HUA Xia ngho ng, QIU Lei. ANew Method for In teger Ambiguity Resolution in GPS Deformati on Mon itori ngJ. Geomatics andIn formation Scie nee

34、 of Wuhan Uni versity, 2007, 32(1):24-26.8 张小红，李征航，徐绍铨.高精度GPS形变监测的新方法及 模型研究J.武汉大学学报(信息科学版),2001, 26(5): 451-454. ZHANG Xiaoh on g, LI Zhe ngha ng, XUShaoqua n. A New Model for High AccuracyDeformati on Mon itor with GPSJ.Geomatics andIn formation Scie nee of Wuhan Un iversity, 2001,26(5): 451-454.9

35、陈永奇，JAMES L.单历元GPS变形监测数据处理方法的研 究J.武汉测绘科技大学学报,1998, 23(4): 324-328,363. CHEN Yongqi, JAMES L. Development of the Methodology for Si ngle Epoch GPS Deformatio nMon itori ngJ.Jour nal of Wuha n Tech ni cal Uni versityof Surveyi ng and Mappi ng, 1998, 23(4): 324-328, 363.10 MENG Xiaolin, ROBERTS G W, CO

36、SSER E, et al.Real-time Bridge Deflecti on and Vibrati onMonitoring Using an IntegratedGPS/Accelerometer/PseudoliteSystemC/Proceedi ngs of the 11th In ter nati onalSymposium on Deformati on Measureme nts,Intern ati onal Federati on of Surveyors (FIG),Commissi on 6-E ngin eeri ng Surveys, Worki ng Gr

37、oup11 杨元喜.北斗卫星导航系统的进展、贡献与挑战J.测绘学报,2010, 39(1): 1 -6. YANG Yuanxi. Progress,Con tributio n and Challe nges of COMPASS/BeiDouSatellite Navigatio n SystemJ.Acta Geodaetica etCartographica Sin ica, 2010, 39(1): 1-6.12 SHI Chua ng, ZHAO Qile, HU Zhiga ng, et al. PreciseRelative Positi oning Usi ng Real T

38、rack ing Data from COMPASS GEO and IGSO SatellitesJ. GPS Solutio ns, 2013, 17(1): 103-119.13 STEIGENBERGER P, HUGENTOBLER U, HAUSCHILD A, et al. Orbit and Clock An alysis of COMPASS GEO and IGSO SatellitesJ. Journal of Geodesy, 2013, 87(6):515-525.14 中国卫星导航系统管理办公室.北斗卫星导航系统空间信号接口控制文件2.0EB/OL. (2013-1

39、2-27)2014-11-15. n.China Satellite Navigati on Office. BeiDou Navigati onSatellite System Signal in Space In terface Con trol Docume nt2.0EB/OL.(2013-12-27)2014-11-15. 15 姜卫平，刘经南，叶世榕.GPS形变监测网基线处理中系统 误差的分析J.武汉大学学报(信息科学版),2001,26(3): 196-199, 238. JIANG Weipi ng, LIU Jingnan, YEShir ong. The Systematical Error An alysis of Baseli neProc