《三种重力场模型在高精度动态测量中的适用性分析》由会员分享，可在线阅读，更多相关《三种重力场模型在高精度动态测量中的适用性分析（14页珍藏版）》请在金锄头文库上搜索。
1、 三种重力场模型在高精度动态测量中的适用性分析 高士健，高淑照，孙诚彬，李 杨，杨 林三种重力场模型在高精度动态测量中的适用性分析高士健1，高淑照1，孙诚彬2，李 杨1，杨 林1（1. 西南交通大学 地球科学与环境工程学院，成都 611756；2. 临沂市公路事业发展中心 莒南县中心，山东 临沂 276600）为了进一步研究重力场模型在工程测量中的适用性，对动态测点GNSS高程转换的精度进行评估：基于某铁路公司实际测量作业，将地球重力场模型2008（EGM2008）、欧洲新技术改进重力场模型6C4（EIGEN-6C4）、实验重力场模型2019e_2159（XGM2019e_2159）3种重力场
2、模型应用在轨道上动态测点的GNSS高程转换中；并采用精度评定指标和多基准站高程互差法、基准站位差互差法进行对比分析。结果表明：1）XGM2019e_2159重力场模型的精度评定指标略优于EGM2008、EIGEN-6C4模型；2）不同模型的多基准站高程互差存在显著差异，2个测段中XGM2019e_2159重力场模型计算的高程互差都能满足限差要求；3）不同模型的基准站位差互差存在较大差异，XGM2019e_2159重力场模型在各测段的基准站位差互差最小。在实际作业中应结合测区选用合适的重力场模型，以保证动态测点的高程转换精度。地球重力场模型；高程异常；适用性分析；后处理动态；高精度动态测量0 引
3、言轨道线形参数的测量是铁路工务部门进行平顺性检测与维护的基础。目前运营普速铁路的轨道线形测量，多采用常规地面测量方法，存在周期长、现势性差、作业效率低、线上工作受天窗时间限制严重等诸多问题，与铁路运营部门的需求存在较大差距，需要对现有的方法和手段进行改进1-3。随着铁路沿线连续运行参考站（continuous operating reference station, CORS）的建设，部分铁路公司在普速铁路轨道线形测量中开始探讨使用全球卫星导航系统（global navigation satellite system, GNSS）后处理动态（post-processing kinematic,
4、 PPK）测量模式。GNSS PPK测量是一种基于载波相位测量的后处理差分定位技术，能够获得厘米级的平面和高程定位结果4-5。在实际作业中，为保证流动站测量结果的准确性，可以通过多个基准站计算结果的比较进行检核与控制6。多基准站GNSS PPK测量模式下，30 km范围内流动站的平面坐标精度能够满足铁路工程测量规范的要求，一般情况下能够达到厘米级，观测条件良好时可以达到毫米级。但高程精度相对较低，一般仅能够达到分米级，采用高分辨率的地球重力场模型可将精度提高至厘米级。流动站的高程是GNSS大地高转换后的正常高，转换所需的高程异常由地球重力场模型格网值内插得到。地球重力场模型在利用GNSS技术确
5、定正常高过程中发挥着至关重要的作用7，是地学领域研究的热点问题。目前国内外已发布170多种模型。地球重力场模型2008（Earth gravitational model 2008, EGM2008）是美国国家地理空间情报局（National Geospatial-Intelligence Agency, NGA）2008年发布的全球超高阶重力场模型，模型球谐系数的阶次扩展至2190阶，空间分辨率达到5（约9 km）8，它是目前使用最广泛的重力场模型之一。欧洲新技术改进重力场模型6C4（European improved gravity model of the Earth by new te
6、chniques 6C4, EIGEN-6C4）是德国地学研究中心（Deutsches GeoForschungsZentrum, GFZ）2014年发布的全球联合重力场模型，其球谐系数的阶数和次数均为2190阶9，它使用了当时最新的卫星数据。实验重力场模型2019e（the experimental gravity field model 2019e, XGM2019e）是慕尼黑工业大学（Technische Universitt Mnchen, TUM）天文和物理大地测量学研究所2019年发布的全球重力场模型，其球谐系数的最高阶次至5399，空间分辨率为2（约4 km），它包含了地球201
7、4（Earth2014）模型中地形数据导出的重力信息；XGM2019e_2159是球谐系数阶次截断至2190阶的重力场模型10。除了理论研究外，重力场模型的适用性验证也具有重要意义。文献11研究了利用重力场模型直接进行全球定位系统（global positioning system, GPS）高程转换的方法，认为在局部区域基于EGM2008模型的GPS高程转换精度可以达到厘米级。文献12采用地球重力场模型与数学函数模型相结合的“移去-恢复”法来反映高程异常的细节变化，其中多面函数与EIGEN-6C4模型相结合的高程转换方案精度最高。文献13利用GNSS/水准数据与不同重力场模型计算的高程异常值
8、进行对比，发现重力场模型计算的高程异常精度在消除系统偏差后精度最高可达到3.94 cm。文献14通过分析模型的内外符合精度和不同阶次组合的差异，选取可靠的截断阶次确定组合重力场模型，结果表明组合重力场模型的高程异常精度最优可达厘米级。文献15以多类重力场模型为基础，经简单谱组合法或加权谱组合法得到组合重力场模型，并利用GNSS/水准数据进行精化，可得到较高精度的区域似大地水准面。文献16利用水准结合GNSS大地高的时变影响、全球大地水准面测量起算面与水准测量起算之间的区域性基准差，提出了GNSS水准地球重力场模型评价方法。为了研究重力场模型在工程测量中的适用性，本文基于某铁路公司轨道线形测量实
9、际项目，在多基准站GNSS PPK测量模式作业中，分别使用EGM2008、EIGEN-6C4、XGM2019e_2159三种重力场模型进行高程转换，采用精度评定指标、多基准站高程互差法和基准站位差互差法，对不同重力场模型的计算结果进行对比分析。1 原理与方法1.1 高程异常计算（1）1.2 模型精度评价方法1.2.1 精度评定指标也可以采用均方根误差（root mean square error, RMSE）、相关系数（correlation coefficient）、纳什系数（Nash-Sutcliffe efficiency coefficient, NSE）等指标检验重力场模型高程异常的
10、质量，即：1.2.2 多基准站高程互差在多基准站GNSS PPK测量模式作业中，使用重力场模型计算的高分辨率（11）高程异常格网值内插动态测点高程异常进行转换，可以采用多基准站高程互差的方法来保证动态测点高程结果的可靠性。由误差传播定律，高程互差的中误差1.2.3 基准站位差互差2 实验与结果分析2.1 数据来源图1 三基站GNSS PPK测量模式应用于轨道线形测量作业示意图表1 测段信息表2 测段基线向量解算精度 m2.2 数据计算按式（2）式（7）计算重力场模型应用于测区内基准站的精度评定指标。图2为EGM2008、EIGEN-6C4、XGM2019e_2159重力场模型位系数误差的阶方差
11、和累积阶方差。图2 3种重力场模型位系数误差的阶方差和累计阶方差表3 重力场模型解算高程异常质量的评定标准图4 测段2三种重力场模型的多基准站高程互差表4 测段1三种重力场模型的多基准站高程互差 m表5 测段2三种重力场模型的多基准站高程互差 m表6 测段高程异常互差 m2.3 结果分析图5 基准站模型高程异常与实测高程异常差值图6 测区模型间高程异常差值等值线由图6可以看出，在数十公里范围内，不同重力场模型计算的高程异常差异显著，进而导致由不同基准站计算的动态测点高程互差超限。重力场是地球质体密度分布的直接映像，结合大地水准面处处与重力方向垂直的性质，大地水准面起伏的中、短波部分与岩石圈内部
12、负荷及地形有很强的相关性。查阅资料可知，测区位于大别山东段的超高压变质带范围内，岩体质量分布可能存在较大的局部异常变化21。与EGM2008和EIGEN-6C4相比，XGM2019e_2159重力场模型包含了地形数据导出的重力信息，因此能够更好地反映该区域的大地水准面分布。进一步收集该区域的相关资料进行研究，也将具有重要的科学价值。3 结束语为了研究重力场模型在GNSS高程转换中的适用性，选取了某铁路公司轨道动态测量的2个测段，分别计算了EGM2008、EIGEN-6C4、XGM2019e_2159三种重力场模型的精度评定指标和动态测点的多基准站高程互差、基准站位差互差，结果表明：在测区范围内
13、，XGM2019e_2159重力场模型的精度评定指标略优于EGM2008、EIGEN-6C4模型,其中XGM2019e_2159重力场模型的均方根误差相比另2个模型提高了28.6%、19.2%；在达到一定阶数后，XGM2019e_2159和EIGEN-6C4模型的位系数误差阶方差和累积阶方差均明显优于EGM2008模型。测段中不同重力场模型计算的多基准站高程互差存在显著差异，EGM2008重力场模型的高程互差均值仅16.7%满足限差要求，最大值达到40.3 cm；EIGEN-6C4模型的高程互差均值有66.7%满足限差要求，最大值为19.7 cm；XGM2019e_ 2159模型高程互差均值能
14、够满足限差要求，最大值仅为5.3 cm。基准站位差互差反映了模型高程异常与实测高程异常的相关程度，即模型高程异常的区域性符合效果。2个测段中XGM2019e_2159重力场模型计算的高程异常互差均值小于6 cm，区域高程异常差值最接近实测值，模型的区域符合效果最佳。在GNSS PPK测量模式中，动态测点的高程、多基准站高程互差和位差互差都与重力场模型计算的高程异常有关，在不同测段选用合适的重力场模型可以将动态测点的高程精度提高至厘米级。基准站点不同重力场模型间的高程异常存在着明显的区域性差异，如何应用重力场模型进行更高精度的高程转换，还须结合各种误差来源做进一步探讨。1 路计哲. 基于GNSS
15、定位技术的轨道几何参数测量方法J. 铁道建筑, 2018, 58(7): 133-136.2 朱洪涛, 龙辉, 吴维军. 一种既有线偏差约束精测网建网方法J. 铁道标准设计, 2017, 61(2): 20-23.3 徐翔. 关于新时代铁路网发展战略选择的思考J. 铁道经济研究, 2021(4): 32-35.4 周东旭, 唐秋华, 张化疑, 等. 星站差分与PPP技术在深远海调查中的位置服务精度分析J. 海洋通报, 2020, 39(2): 215-222.5 杜锐. GPS PPK技术在远海水下地形测量中的应用J. 全球定位系统, 2017, 42(1): 111-113.6 白立舜, 张宏伟, 聂敏莉. GPS PPK技术和GPS RTK技术在包头市第二次土地调查中的应用与对比分析J. 测绘通报, 2013(7): 49-53.7 李玉平, 闫志学, 裴佳佳, 等. 不同地形条件下超高阶重力场模型精度比较分析J. 工程勘察, 2019, 47(8): 62-66.8 吴寒, 吴燕苹, 邓明镜. 综合EGM2008和移动曲面模型的GNSS高程拟合J. 全球定位系统, 2020, 45(1): 77-81.9 贾雪, 徐炜, 刘超, 等. 顾及地球重力场模型的GNSS高程转
