GPS变形监测中提取大变形量的单历元解算方法研究.docxVIP

GPS变形监测中提取大变形量的单历元解算方法研究 广西聚海工程勘察设计有限公司 摘要：提出一种综合利用宽巷组合观测值与L1观测值提取测站大变形 量的GPS单历元变形监测方法。充分利用变形监测点已知以及模糊度之间的线性 相关性等特点，建立了按历元解算监测点变形量的基木数学模型，并分析了影响 变形量精度的主要误差来源。通过对试验数据的计算分析，对数学模型的正确性 进行了论证。 关键词：GPS；变形监测；整周模糊度；单历元；宽巷观测值 GPS由于其高精度和连续自动监测能力，已广泛应用于结构物变形监测。 根据被监测对象的特点，GPS变形监测分为周期性重复测量、固定连续GPS测站 阵列和实时动态监测3种监测模式。对于大坝、水库、滑坡体、构造运动等变形 连续、缓慢、量级小的结构物，数据处理时宜采用静态解模式，通过多个历元的 观测资料解算。而大桥、高层建筑物等结构物，在外部荷载作用下摆动变形量级 较大且剧烈，应采用实时动态监测数据处理模式。目前，实时动态监测数据处理 主要采用。 OTF （或ARCTF）方法，常用的有最小二乘搜索法（LS）、快速模糊度搜 索法（FARA）、最小二乘模糊度降相关平差法（LAMBDA）、Choirsky分解搜索计 算法等。这些方法的模糊度解算均需多个历元的观测数据，且需进行周跳的探测 与修复。当卫星失锁时OTF方法必须重新初始化，使监测过程变得十分复杂。部 分学者提出了实时动态监测的单历元解算方法，充分利用GPS变形监测系统中各 监测点精确己知的信息，从GPS数据处理的观测方程入手，避免了常规处理中的 周跳探测与模糊度分离等步骤。但该类算法大多仅适用于瞬间变形值较小的情况, 若变形监测点产牛不规则的大变形，如悬索桥在风力荷载作用下的摆动，这些方 法则失效。木文提出一种适合于提取大变形量的GPS变形监测单历元解算算法。 该算法先利用宽巷观测值对大变形量进行初步提取，将监测点初始坐标进行粗略 改正，然后利用L1相位观测值采用“无周跳无模糊度”算法进一步提取剩余变 形量。分析表明，在监测点产生大变形吋该算法能够取得很好的效果。 1单历元变形监测的数学模型 1.1利用宽巷组合获取变形量模型 由于变形监测中基线一般较短，经过双差后的电离层、对流层残差影响 能够得到很好的消除，多路径效应影响减弱。另外，当设定一定的截止高度角以 及采用扼流圈技术时，多路径效应影响可以认为为0,于是GPS双差相位观测方 程可写为： 提出一种“无周跳无模糊度”的监测方法，认为当变形量小于lambda;/ 2 （lambda;为载波波长）吋，就能根据计算出的双差值中的整周数来确定观测 的双差值中的整周数，从而避开“周跳探测与修复”、“整周模糊度确定”等棘手 问题。如图1 , A为基准站，B为监测站，Bprime;为变形后的B点位置，p和 q为同步观测卫星，以高度角较大的卫星q 为参考星，则A、B两站的相位双差观测值为： 在变形之前，监测站B的坐标已知，可以利用A、B两点的坐标及卫星 星历计算出的卫星位 置求得双差距离： 由式（3）可计算出双差相位： 在短基线数据处理中，经过双差后电离层、对流层、多路径效应等误差 得到很好的消除，则式（2）、（4）代入式（1）组成双差观测方程: 由此可见，当B点坐标精确已知时，双差模糊度Delta;N可以通过 Delta;phi; c — Delta;phi; 0 直接取整来获得： 然而，若监测站B发生变形，产生了位移delta; u至出点，仍然利用B 点坐标计算双差观测值就不能满足式（5）。它们的差值可以用 表示： 则中包含了变形信息。根据“无周跳无模糊度”算法的基本思想，当 变形量 小于载波半波长lambda;/ 2吋，式（7）等号右端 直接取整确定的双 差模糊度Delta;N仍为正确的双差模糊度，即变形量并没有对相位双差观测值 的整数部分造成影响，仅影响了观测值的小数部分。因此，中仅包含了不足一 波段的小数部分，确切地说， [-0.5, 0.5]o据式（1）列出误差方程如下： 而如果变形值大于lambd可/2、小于lambda;,对某些高度角较低的 卫星来说，变形量对 的影响可能会大于0.5周，那么 应为绝对值大于0.5周的 小数值。如果仍然利用“无周跳无模糊度”的解法，对 直接取整（即四舍五入）， 求得的双差模糊度会与正确的Delta;N存在1周的偏差，此吋通过式（7）计算 得到的绝对值仍然在[-0.5, 0.5]之间，而符号也发生了变化。那么，该颗卫星 所列出的方程对该历元整个方程组造成污染，经过最小二乘平差就得不到正确的 变形量。为了能够更进一步地提取较大变形量，本文提出当变形量大于lambda;/ 2、小于lambda;时，对利用监测点近似坐标计算岀的双差模糊度进行plusmn;l 周的搜索，确定正确