地球物理:大地测量学综述002

球谐与椭球谐函数在物理中的应用解析
本文深入探讨了球谐与椭球谐函数的基本概念及其在物理学中的广泛应用,包括量子力学中的原子轨道描述、地球物理学中的地球磁场模型以及电磁场理论中的场分布计算等。通过实例解析,阐述了这两种谐函数如何帮助科学家理解和解决复杂问题。

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
《 球 谐 与 椭 球 谐 》 《球谐与椭球谐》
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

高原峡谷区“移去-恢复法”高程拟合可靠性分析 摘要:高原峡谷地区地形起伏剧烈,重力场变化复杂,给GNSS高程拟合带来了巨大挑战。本文以云南省某个典型高原峡谷矿山为实验区,基于XGM2019e超高阶全球重力场模型,系统性地比较了2种高程异常拟合方案:(1)直接使用拟合点高程异常值进行拟合;(2)利用“移去-恢复法”结合XGM2019e模型和拟合点进行拟合。通过各项精度指标对比分析,表明“移去-恢复法”能有效削弱系统误差,其拟合精度显著优于纯几何拟合的方法,验证了该方法在高山区进行高精度高程拟合的可靠性与实用性。 关键词:移去-恢复法;XGM2019e;高程异常;高程拟合;精度分析;高原峡谷 0引言 高程系统是现代大地测量学地球物理学及各类工程建设的基础。GNSS技术虽能精确测定点的三维坐标,但其提供的大地高(H)是以参考椭球面为基准的,与工程建设中通常使用的以似大地水准面为基准的正常高(Hr)存在高程异常(ζ)差异(ζ = H - Hr)。如何精确求定高程异常,实现大地高向正常高的转换,是GNSS技术替代传统繁重水准测量、全面提升测绘效率的关键【1】。 在平坦地区,高程异常变化平缓,可采用平面或低次曲面拟合达到厘米级精度。然而,在高原峡谷地区,剧烈的构造运动、巨大的高差起伏和复杂的地质结构导致了重力场的强烈变化,使得高程异常场呈现出高频、非线性特征【2】。“移去-恢复法”是解决这一问题的有效策略。其核心思想是:利用全球重力场模型精确计算并移除高程异常的中长波分量(趋势项),再利用局部的GNSS/水准数据,通过数学拟合方法对残余的短波分量(细节信号)进行建模恢复【4】。该方法既利用了GGM的全局稳定性,又结合了局部数据的细节补充能力,被公认是区域高程异常求精的最佳方法之一。 XGM2019e是由欧洲空间局等相关机构于2019年联合发布的最新一代超高阶全球重力场模型,其阶次高达2190【5】。它融合了GOCE、GRACE、GRACE-FO卫星重力数据,以及最新的地面、船测和航空重力数据【6】。然而,XGM2019e在特定复杂地形区的具体表现,尤其是与“移去-恢复法”结合后的实际效能,仍需通过详实的实证研究进行验证。 本文选取一个具有代表性的云南某高原峡谷矿区作为实验区,开展基于XGM2019e模型的“移去-恢复法”高程拟合,并从多个维度评估该技术路线在极端地形条件下的可靠性。 1. 实验区与数据源 1.1 实验区概况 实验区位于云南省东北部金沙江流域,是典型的高原深切峡谷地貌,区域相对高差达1400米以上(见图1)。区域内山体陡峭,坡度多在30°以上,V型谷特征明显。如此剧烈的地形起伏,导致地球重力场发生显著畸变,高程异常场变化复杂,对高程拟合的精度构成了严峻挑战,是检验“移去-恢复法”可靠性的理想实验场。 图1 实验区地形 Fig.1 1.2 数据来源与处理 实验所需数据主要包括以下几类: (1)GNSS和水准点数据:实验区内矿山控制测量已布设了31个高等级GNSS控制点,所有点均联测了国家四等及以上水准,已知大地高和正常高成果,点位较均匀的分布在山谷、坡地、山脊等不同地形单元。实验将31个已知点分为两组,其中拟合点16个,检核点15个(见图2)。 (2)全球重力场模型:采用XGM2019e全球重力场模型,计算实验区内任意点的高程异常模型值。 图2 点位分布图 Fig.2 2. 实验方法 为全面评估“移去-恢复”法的效果,设计以下2个实验方案。 方案一(纯几何法):直接使用16个拟合点的实测高程异常ζ,采用克里金插值法构建区域高程异常模型,内插出拟合点和检核点的高程异常值ζ_direct。 方案二(“移去-恢复法”),结合XGM2019e全球重力场模型和拟合点高程异常进行拟合。其主要过程为: (1)利用XGM2019e模型计算拟合点的模型值ζ_GGM。 (2)利用拟合点的真实值ζ_meas和模型值ζ_GGM计算拟合点的残差 Δζ = ζ_meas - ζ_GGM。 (3)利用拟合点的坐标和Δζ,采用克里金插值法构建区域残差场模型。 (4)对于任意待求点,计算其模型ζ_GGM并内插出其Δζ,最终高程异常 ζ_RR = ζ_GGM + Δ_ζ。 3.实验结果分析 3.1模型在该区域的整体表现 使用XGM2019e模型计算所有点的高程异常值ζ_GGM,与实测值ζ进行比较,通过统计,其残差最大值、最小值、平均值、标准差和RMSE结果如表1所示。 表1 模型高程异常残差统计 Tab.1 最大值 最小值 平均值 标准差 RMSE 0.629 0.221 0.444 0.097 0.454 这表明即使是超高阶重力场模型,在如此复杂的地形下,也无法独立达到厘米级精度要求。在高原峡谷这类复杂地形区,受模型误差、地形效应及密度异常影响,其高程异常计算值存在明显的系统偏差和局部偏差,难以直接满足高精度高程转换的需求。 3.2不同方案精度对比 2个方案检核点拟合成果残差图如图3所示。 图3 不同方案检核点拟合成果残差图 Fig.3 分别将2个方案拟合点和检核点上的拟合值与实测值比较,统计得到内符合精度和外符合精度(见表2)。 表2 不同方案内、外符合精度统计(单位:米) Tab.2 方案选择 内符合精度 外符合精度 方案一(纯几何法) 0.019 0.431 方案二(“移去-恢复”法) 0.001 0.022 采用“移去-恢复法”的方案,其精度相比纯几何法有数量级上的提升(从分米级到厘米级),充分证明了该方法能显著提升高程拟合精度,拟合结果完全满足高山区似大地水准面精度的要求。 3.3 残差分布与地形相关性分析 为进一步分析误差来源,绘制了方案二检核点上的残差分布图,并将其与地形图叠加。发现残差绝对值较大的点(>2cm)主要分布在两个区域:一是坡度很陡的坡地,二是山谷附近。这暗示了误差可能与地形复杂度(坡度、坡向等)密切相关。全球模型和数学拟合均难以完全模拟极端地形引起的超短波重力效应,这可能是限制精度进一步提升的主要原因,也指明了未来引入高精度地形改正的必要性。 4. 结论 本文通过对一个典型的高原峡谷实验区进行实证分析,得出以下主要结论: (1)可行性验证:研究证实,基于XGM2019e全球重力场模型的“移去-恢复”法,完全适用于高原峡谷这类极端复杂的地形区域,并能获得高精度的高程拟合结果。 (2)精度可靠性:该方法将高程异常拟合的外符合精度从纯几何法的±0.431m大幅提升至±0.022m,其成果可靠性高,能够满足高等级工程测量和科学研究的需要。 (3)效能优势:“移去-恢复”法的核心优势在于物理与几何的融合。XGM2019e模型提供了较为准确的物理背景场,而局部几何拟合则灵活地补偿了模型的局部偏差和短波信号。二者相辅相成,克服了各自单独使用的局限性。 综合来看,“移去-恢复法”结合XGM2019e模型与克里金插值,是高原峡谷区GNSS高程拟合的有效且可靠的技术途径,在实践中应优先采用。未来的研究可聚焦于更多不同规模、不同地质背景的高原峡谷区,融合高分辨率的局部地形和重力数据,以进一步优化残差场的拟合效果。 参考文献 【1】 宁津生, 罗志才, 李建成. 现代大地测量基准及其进展[J]. 测绘学报, 2017, 46(10): 1555-1564. 【2】 李建成, 褚永海, 邢乐林, 等. 高程异常计算的误差分析与精度评估[J]. 武汉大学学报(信息科学版), 2012, 37(9): 1021-1025. 【3】 Pavlis N K, Holmes S A, Kenyon S C, et al. The development and evaluation of the Earth Gravitational Model 2008 (EGM2008)[J]. Journal of Geophysical Research: Solid Earth, 2012, 117(B4). 【4】 翟振和, 孙中苗. 移去-恢复技术在多波束测深数据融合中的应用[J]. 测绘科学, 2016, 41(3): 36-40. 【5】 Zingerle P, Pail R, Gruber T, et al. The combined global gravity field model XGM2019e[J]. Journal of Geodesy, 2020, 94(7): 1-12. 【6】 徐新禹, 詹总谦, 李星星, 等. XGM2019e全球重力场模型在中国地区的精度评价[J]. 测绘地理信息, 2022, 47(1): 1-6. 请将以上论文扩充到5000字以上。
最新发布
09-11
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值