[发明专利]一种基于重复线校正海空重力仪格值的新算法

摘要

本发明涉及一种基于重复线校正海空重力仪格值的新算法，其主要技术特点是：依据航空重力测量的基本数学模型，建立重复测点处格值校正模型；由最小二乘原理求得格值偏差的最或然估值，并对海空重力仪格值进行补偿。本发明通过分析海空重力测量系统误差的形成机理，发现海空重力仪格值标定误差是引起系统性测量偏差的主要因素之一，利用重复测线检测校正海空重力仪格值的计算模型和补偿方法，利用航空重力实际观测网数据对该方法的合理性和有效性进行了数值验证，试验表明本发明对消除海空重力测量系统性偏差具有显著作用，该方法简单，能够满足海空重力测量用户的需求。

主权项

1.一种基于重复线校正海空重力仪格值的新算法，其特征在于包括以下步骤：步骤1、依据航空重力测量的基本数学模型，建立重复测点处格值校正模型；步骤2、由最小二乘原理求得格值偏差的最或然估值，并对海空重力仪格值进行补偿；所述步骤1的具体实现方法为：设在东西方向上布设两条正反向的第1重复测线和第2重复测线，在第1重复测线和第2重复测线的重复测点P处的绝对重力值分别表示为：g_p1＝g_b+(CS_p1‑CS₀)‑δa_V1+δa_E1+δa_H1+δa_F1+δa_A1g_p2＝g_b+(CS_p2‑CS₀)‑δa_V2+δa_E2+δa_H2+δa_F2+δa_A2此时，重复测点P处的重力互差为：δg_p12＝CS_p1‑CS_p2‑δa_V1+δa_V2+δa_E1‑δa_E2+δa_H1‑δa_H2+δa_F1‑δa_F2+δa_A1‑δa_A2其中，g_b代表比对基点处的重力值，C为格值；所述步骤2的具体实现方法为：假设不符值的趋势性部分主要由仪器格值偏差ΔC所致，其随机性部分则归结为动态环境和各项改正建模误差的综合影响Vg_p12，得到如下观测方程：δg_p12+ΔC(S_p1‑S_p2)+Vg_p12＝0对应于重复测线上的每一个重复测点均建立观测方程，当存在n个重复测点时，分别以ΔC作为待定参量，以重复测点i处的重力互差δg_i12作为观测量，以重复测点i处的综合影响Vg_i12作为观测误差的改正数，则由最小二乘原理可求得格值偏差ΔC的最或然估值如下：估值的精度估算式为：其中：求得格值偏差估值后，按下式计算海空重力仪格值的校正值：使用校正值代替原值C重新计算各个测点的绝对重力值，从而消除格值偏差给海空重力测量成果带来的系统性影响。