[发明专利]一种基于重力异常计算外部扰动重力北向分量的方法

说明书

本发明涉及一种基于重力异常计算外部扰动重力北向分量的方法，其技术特点是：建立基于重力异常计算外部扰动重力北向分量的传统全球积分式：利用移去‑恢复技术，从重力异常中移去参考重力异常，得到残差重力异常；从积分式核函数中移去与参考场对应阶次的核函数球谐表达式，得到截断核函数，使之与残差重力异常的频谱匹配；基于截断核函数和残差重力异常的局域积分，得到残差外部扰动重力北向分量；利用全球位场模型高阶信息进行远区效应补偿；恢复参考外部扰动重力北向分量，得到计算点处高精度的外部扰动重力北向分量。本发明能够精确计算外部扰动重力北向分量，提高了外部扰动重力北向分量的计算精度，可广泛用于物理大地测量领域。

技术领域

本发明属于大地测量学与测量工程技术领域，尤其是一种基于重力异常计算外部扰动重力北向分量的方法。

背景技术

外部扰动重力北向分量是地球重力场逼近建模研究内容的重要组成部分，也是解算大地测量边值问题的主要应用目标之一，在航空航天器飞行轨迹精密计算和空间科学技术研究中具有重要的应用价值。

基于重力异常计算外部扰动重力北向分量的传统积分式要求全球积分，但实际应用中受观测数据覆盖范围限制，无法做到全球覆盖，实际计算过程中外部扰动重力北向分量的全球积分式需要进行适用观测数据保障条件的改化，以保证计算结果的可靠性。目前，迫切需要解决基于重力异常计算外部扰动重力北向分量的传统积分式与实际应用数据保障不匹配的问题，以提高外部扰动重力北向分量的计算精度。

发明内容

本发明的目的在于克服现有技术的不足，提出一种基于重力异常计算外部扰动重力北向分量的方法，解决基于重力异常计算外部扰动重力北向分量的不确定问题，提高外部扰动重力北向分量的计算精度。

本发明解决其技术问题是采取以下技术方案实现的：

一种基于重力异常计算外部扰动重力北向分量的方法，包括以下步骤：

步骤1、建立基于重力异常计算外部扰动重力北向分量的传统全球积分式：

步骤2、顾及实测重力异常数据保障条件，引入全球位场模型，利用移去-恢复技术，从重力异常中移去参考重力异常，得到残差重力异常；从积分式核函数中移去与参考场对应阶次的核函数球谐表达式，得到截断核函数，使之与残差重力异常的频谱匹配；基于截断核函数和残差重力异常的局域积分，得到残差外部扰动重力北向分量；利用全球位场模型高阶信息进行远区效应补偿；恢复参考外部扰动重力北向分量，得到计算点处高精度的外部扰动重力北向分量。

而且，所述步骤1建立的传统全球积分式为：

式中，Δg为球面上流动点q处的已知观测重力异常；σ为单位球面；dσ为单位球面的面积元；R为地球椭球平均半径；r为计算点地心向径；为计算点的纬度和经度；为流动点的纬度和经度；ψ为计算点至流动点之间的球面角距；是计算点至积分流动点之间的空间距离；F_ψ(r,ψ)为积分核函数；α为流动点到计算点的方位角。

而且，所述步骤2得到的计算点处高精度的外部扰动重力北向分量的公式为：

其中δΔg为残差重力异常；为截断核函数；为远区效应计算值；为参考扰动重力北向分量；

所述残差重力异常δΔg的计算式为：

δΔg＝Δg-Δg_ref

其中Δg_ref为由N阶参考场位模型计算得到的参考重力异常，其计算公式为：