北京市地面沉降监测高精度数据处理及分析

为了获取高精度的监测点高程变化信息，在借鉴以往观测数据成果、并顾及到点位分布的情况下，布设监测网时，以10个监测点为一个同步环，同步环之间的公共连接点均不少于3个，共组成13个同步环，每年外业观测时间从8月底开始到10月底结束，持续时间近2个月。观测的仪器采用的是天宝的双频接收机和扼流圈天线。

2 沉降监测网数据处理

数据处理主要包括基线解算与网平差工作。本文数据为北京市2020年的观测数据。GPS基线解算采用GAMIT/GLOBK软件包，其作为高精度的数据处理软件，广泛应用于地壳变形的分析研究、高精度GNSS数据处理等。GAMIT/GLOBK软件自动化程度高，通过选取适当数量的国际GNSS服务（International GNSS Service, IGS）站，并使用精密星历解算GPS数据，可求得高精度基线向量和测站坐标，是目前国际上比较流行且解算精度较高的GPS数据处理软件之一[9]。

2.1GPS测量数据的预处理

在高精度GPS数据处理前，需要对GPS原始观测数据进行预处理，包括GPS测量点的变动、GPS天线高的检查、GPS观测数据的质量分析。

GPS测量点的变动检查：主要核对在监测周期内GPS监测点的点位是否发生过移动、GPS接收机和天线是否更换、GPS观测数据是否有效、监测点点名与其实际位置是否一致等。

GPS天线高的检查：高精度地面沉降GPS测量需要有统一的起算基准，针对观测墩的特点，采用在GPS监测墩顶面彼此间隔120°位置上刻一浅槽的量测模式方案，使浅槽底部基本保持在同一水平面上，并将浅槽打磨光滑。GPS天线高的量测是从刻槽的位置垂直向上至GPS扼流圈天线圆盘外沿的最底部（图1）。在每个观测时段测前、测中和测后各量取一次，3次量取的天线高较差在1mm以内时，取其平均值作为最后天线高。

GPS观测数据的质量分析：GPS观测数据的可靠性和完好性是后续进行高精度数据处理的重要的保障，因此在数据处理前应对GPS原始观测数据进行质量分析。采用GPS观测数据质量分析软件TEQC（translation, editing, and quality checking）对本文的GPS观测数据质量进行了全面分析。TEQC软件是通过计算观测数据的多路径情况（MP1和MP2）、数据的完整性（实际观测值的数据量除以理论上的观测值数）和信号的失锁情况等内容来评价GPS监测测站周围的环境对GPS观测值的影响程度。本文对采集的GPS监测网原始观测数据进行了质量检查，并对多路径影响MP1、MP2值做了详细统计，下面主要给出MP1、MP2指标分布图（图2至图3）。

2.2基线解算

基线解算是将基准网点与监测网点分别通过与3个北京地面沉降基岩基准点(BJFS、BJSH、JIXN）联网计算，进而将监测网的基线分别纳入国际地球参考框架2005（international terrestrial reference frame2005, ITRF2005）中。利用全球参考框架的稳定性与连续观测数据的完好性，以及基岩基准点的高精度三维地心坐标来进一步提高GPS基线解的精度。

GAMIT基线解文件按单时段天解的形式给出，其基线解为双差固定解。基线解文件的主要内容包括该时段所包含的测站名称、观测历元总数、采样间隔、工作开始时间、各GPS基线在空间直角坐标系与站心坐标系中的解算结果以及基线解的方差-协方差矩阵等[6-8]。

2020年度GPS测量的基线解算采用解算策略：将102个监测网点、3个基岩基准点（BJFS、BJSH、JIXN）综合在一起进行基线解算，获取监测网的基线向量结果。GAMIT主要参数设置如表1所示[9-14]。

2.3网平差方案

监测网的基线向量网平差主要采用拟稳平差。网平差中均以3个北京地面沉降基岩基准点（BJFS、BJSH、JIXN）作为起算基准。

平差方案采用拟稳平差主要是为了避免因对基准点作强约束致使平差观测量受到歪曲而影响精度，同时也为了不会因完全放弃固定基准，造成不能提供所求的变形量。拟稳平差法即对网中相对稳定3个或多个基准点，在最小范数条件下作拟稳秩亏自由网平差，使稳定未知量拟合于其稳定值，而对整网在最小二乘条件下作经典平差。选用拟稳平差方案，网中既保存了相对稳定的变形监测参考基准，可以较好地达到监测变形的目的，又不会使平差坐标值因基准点的少量变形和误差影响而产生扭曲，保持了网中点坐标的精度。

2.4精度分析

GPS测量的精度评定从GPS基线重复边、环闭合差以及大地高精度3个方面进行。

1）基线重复边。以GPS监测网基线重复数来评定基线解算的精度。基线重复数能够反映时段解(单天解)之间的内符合精度，是评定GPS定位结果的重要指标之一，计算公式为

由表2统计结果可以看到，GPS监测网基线重复数精度全部满足限差要求。

2）同步环和异步环闭合差。GAMIT软件采用的是网解（即全组合解），其同步环闭合差在基线解算时已经进行了分配。对于GAMIT软件基线解的同步环检核，可以把基线解的标准化均方根（normalized root mean square, NRMS）值作为同步环质量好坏的一个指标，一般要求NRMS值小于0.3mm，表3给出了监测网的NRMS值的统计信息。由表3中统计结果可以看出，所有基线NRMS值均小于0.2mm，这说明GPS同步基线向量网的内符合精度较高。

同步环闭合差和异步环闭合差的值能够反映GPS外业观测质量以及基线解算质量的可靠性。同步环闭合差可以反映出一个同步环数据质量的好坏，而异步环闭合差则可以反映出整个GPS网的外业观测质量和基线解算质量的可靠性，相比较于同步环闭合差，异步环闭合差对GPS成果质量更为重要。

表4给出了异步环闭合差的统计结果，异步环闭合差主要反映GPS向量网自身的内符合精度，特别是重复基线是否存在异常或超限。依据GPS规范来检测异步环闭合差是否超限。由表4结果可知，GPS测量网异步环闭合差全部合格。

3）GPS站点大地高精度分析。GPS测量成果精度主要通过平差后GPS点大地高精度进行分析。由于北京市地面沉降GPS监测站点形变信息主要基于各站点大地高H获取，因此GPS站点大地高H精度与北京市地面沉降密切相关，其精度高低直接影响到北京市地面沉降GPS监测的可靠性。图4为监测点大地高精度分布图，通过对监测点大地高精度进一步统计得到：2020年监测网大地高误差最小值为0.5mm，最大值16.9mm，平均值3.6mm。

表5为GPS监测网点大地高精度区间分布统计表，从表中可知大地高精度小于10mm的占比为95.1%。

3 沉降监测结果分析

根据上述的沉降监测网数据处理方案，在同一参考基准框架下，由2019、2020年两期GPS测量平差结果，得到各GPS监测点相应于监测周期的大地高变化量。根据监测点的大地高变化量得到了2019—2020年监测点高程变化量分布图，见图5。

从图5中可以看到，2019—2020年高程变化量较大的区域：昌平、顺义及海淀区。在这些区域中，多数的监测点年度高程变化量在0～50mm之间，少量的监测点年度高程变化量大于50mm。同时也可以看到，海淀区的监测点数量较少，尤其是在和昌平交界的海淀北部区域。市区东部的朝阳、通州地区。大兴的礼贤、榆垡区域。综上可知，基于GNSS对北京市地面沉降的监测结果较好地反映出北京年度的地面沉降变化信息。

4 结束语

空旷的通视环境是保证高精度GPS监测的首要条件，但是北京市的城市发展和环境变化却严重影响了部分GPS站点的观测环境，导致少量点的精度比较差。为了确保北京市GPS地面沉降监测工作的高效性、稳定性和连续性，可以从以下2个方面改进：尽可能避开植被茂盛的季节进行GPS监测，减弱GPS多路径误差，提高测量精度；对地面沉降发育区及时进行GPS站点的修建。

全球卫星导航系统在北京市地面沉降监测中已应用多年，取得了很好的监测成果，对全方面了解北京市地面沉降区域发展动态、发展规律，研究北京市平原区内部平面运动、发现地裂缝发展区域提供了基础。同时，我们也应该不断推动自主建设、独立运行的北斗卫星导航系统（BeiDou navigation satellite system, BDS）在地面沉降、地裂缝监测等地质灾害监测领域的应用。

作者简介：

曹炳强1，刘智强1，简程航1，鲁泽宇1，张双成2，罗 强1

（1. 北京市地质工程勘察院，北京 100048；2. 长安大学，西安 710054）

曹炳强（1990—），男，山东潍坊人，硕士，工程师，研究方向为GNSS数据处理及工程监测。

引文格式：曹炳强，刘智强，简程航，等. 北京市地面沉降监测高精度数据处理及分析[J]. 导航定位学报, 2021, 9(6): 125-129.（CAO Bingqiang, LIU Zhiqiang, JIANChenghang, et al. High precision data processing and analysis of landsubsidence monitoring in Beijing[J]. Journal of Navigation and Positioning, 2021, 9(6): 125-129.）DOI:10.16547/j.cnki.10-1096.20210619.

- END -返回搜狐，查看更多