一种应用于船载重力仪的格值校正方法技术

本发明专利技术公开一种应用于船载重力仪的格值校正方法，以船载重力仪往返启航港且中途停靠补给港的闭合重力数据资料为基础数学模型，建立观测重力值的误差补偿校正模型，以对观测重力值进行误差补偿；并基于补偿后的重力观测值与绝对重力值的对应关系求格值，并对船载重力仪数据进行补偿；能够消除由格值误差引起的系统误差，并获得高精度的航次调查成果数据，对海洋地球物理调查，尤其是深远海航次的海洋重力测量具有重要指导价值和实践意义。

A method of grid value correction applied to shipborne gravimeter

【技术实现步骤摘要】
一种应用于船载重力仪的格值校正方法
本专利技术属于船载重力仪观测数据系统差校准及数据处理的应用领域，具体说是一种应用于船载重力仪的格值校正方法。
技术介绍
作为地球重力测量学研究的重要组成部分，海洋重力测量为研究地球形状和地球内部构造、保障航天和远程武器发射、勘探海洋矿产资源等提供重力资料。船载重力仪作为海洋地球物理调查的关键核心设备，除受外界坏境及现场操作人员水平等因素影响会引起各种系统和偶然误差外，往往在设备使用过程中，船载重力仪本身的格值精度也会发生变化。此前受造船技术及船舶续航等因素影响，海洋重力调查主要集中在近海区域，往返统一港口或是就近的港口，重力仪格值在购买之初也会由重力仪生产商标定给出。在实际调查过程中，通过往返母港的方式可以轻易判别船载重力仪的零点漂移情况，但是格值是否正确及其精度如何往往容易被主观忽略。随着全球深远海航次的逐年增多，越来越多的航次在调查过程中需停靠多个码头，客观上具备长基线观测的条件及要求。目前国内海洋重力调查仍大范围使用国外的海洋重力仪，传统方法是将重力仪运回至生产地进行校准或选取陆地的两个长基线点进行校准，然而船载重力仪作为一种精细的仪器设备极为容易在搬运途中受碰撞等因素影响，导致格值的二次变化。随着科学技术的不断发展，深远海重力调查越来越普遍，中间往往会停靠一个或多个码头，基于这种情况，传统方法并不能有效解决实际问题，故而本方案提出一种基于多港口停靠的重力数据模型的船载重力仪格值校正方法。
技术实现思路
本专利技术提出一种应用于船载重力仪的格值校正方法，基于船载重力仪往返启航/返航港且中途停靠补给港的重力数据资料为数学模型，校正计算单个航次内船载重力仪的格值，有效消除由格值误差引起的误差，并获得高精度的航次调查成果数据。本专利技术是采用以下的技术方案实现的：一种应用于船载重力仪的格值校正方法，包括以下步骤：步骤S1、以船载重力仪往返启航/返航港且中途停靠补给港的重力数据为基础数学模型，建立观测重力值的误差补偿校正模型，以对观测重力值进行误差补偿；步骤S2、基于补偿后的观测重力值与绝对重力值的对应关系求取格值，并对船载重力仪数据进行补偿。进一步的，所述步骤S1中对观测重力值进行误差补偿时具体包括以下步骤：步骤S11、分析并分别获得启航码头、补给码头和返航码头的重力观测点与对应的各码头基点由于纬度变化导致的正常重力场值变化；步骤S12、将启航码头重力观测数据归算至重力基点高程面，得到归算后的启航码头的重力基点处的重力仪数据。进一步的，所述步骤S11具体通过以下方式获得由于纬度变化导致的正常重力场值变化：gL＝0.000814*sin2L*dL其中：其中：dL为船载重力仪与基点的南北向距离，单位为米(m)；a为码头(启航码头、补给码头和返航码头)与正北方向夹角，rad；gL为正常重力场变化，10-5m/s2；L为重力观测点纬度，rad。进一步的，所述步骤S12中，采用以下公式进行计算：G基启＝G启+0.3086*h/K0+gL/K0其中：h为船载重力仪重心到基点的高度，单位为米(m)；A为船左舷到水面距离，单位为米(m)；B为船右舷到水面距离，单位为米(m)；C为船载重力仪重心到甲板的距离，单位为米(m)；D为码头重力基点到水面距离，单位为米(m)；G启为启航码头的重力观测值，格；G基启为归算到重力基点处的重力仪数据，格；K0船载重力仪原始格值，10-5m/s2/格。进一步的，所述步骤S2具体包括以下步骤：步骤S21、基于步骤S1分别分析获得补给码头和返航码头归算至重力基点处的重力仪数据G基补和G基返，采用启航码头和返航码头的重力观测数据计算获得零点漂移率，并补偿补给码头的观测重力值G；步骤S22、计算获得船载重力仪格值K：其中：K为执行本次计算之后的船载重力仪格值，10-5m/s2/格；G绝启、G绝补分别为启航码头和补给码头基点处的绝对重力值，采用绝对重力仪直接测得。进一步的，所述步骤S21具体采用以下方式分析：其中：t1、t2、t3分别为重力观测值G启、G补、G返的观测时间，G补为补给码头的重力观测值，G返为返航码头的重力观测值。进一步的，所述步骤S2之后还包括：步骤S3：基于步骤S2获得的格值，采用多次迭代步骤S1至步骤S2的方式求取更新后的格值，直至相邻两次迭代所得到的格值K的差值绝对值小于允许误差为止，进而得到校正后的格值。与现有技术相比，本专利技术的优点和积极效果在于：本方案以船载重力仪往返启航/返航港且中途停靠补给港的重力数据资料为基础数学模型，建立观测重力值的误差补偿校正模型，以对观测重力值进行误差补偿；并基于补偿后的重力观测值与绝对重力值的对应关系求格值，并对船载重力仪数据进行补偿；能够消除由格值误差引起的系统误差，并获得高精度的航次调查成果数据，对海洋地球物理调查，尤其是深远海航次的海洋重力测量具有重要指导价值和实践意义。附图说明图1为本专利技术实施例中的船舶往返及基点比对示意图；图2为本专利技术实施例中格值校正前船测数据与卫星数据对比图；图3为本专利技术实施例中格值校正后船测数据与卫星数据对比图。具体实施方式为了能够更清楚的理解本专利技术的上述目的和优点，下面结合附图对本专利技术的具体实施方式做详细地描述，需要说明的是，本实施例中，启航港和返航港可以是同一港口，也可以是绝对重力值差值较小的两个不同港口；启航码头为启航港内的某一码头，返航码头和补给码头同理。本实施例提出一种应用于船载重力仪的格值校正方法，包括以下步骤：步骤一、以船载重力仪往返启航港且中途停靠补给港的重力数据资料为基础数学模型，建立观测重力值的误差补偿校正模型，以对观测重力值进行误差补偿，具体包括以下步骤：(1)分析并分别获得启航码头、补给码头和返航码头的重力观测点与对应的各码头基点由于纬度变化导致的正常场值变化：gL＝0.000814*sin2L*dL其中：dL为船载重力仪与基点的南北向距离，单位为米(m)；a为码头(启航码头、补给码头和返航码头)与正北方向夹角，rad；gL为正常重力场变化，10-5m/s2；L为重力观测点纬度，rad。(2)将启航码头(重力观测值G启)重力观测数据归算至重力基点高程面，采用如下公式计算：G基启＝G启+0.3086*h/K0+gL/K0其中：h为重力仪重心到基点的高度，米(m)；如图1所示，A为船左舷到水面距离，米(m)；B为船右舷到水面距离，米(m)；C为重力仪重心到甲板的距离，米(m)；D为重力基点到水面距离，米(m)；G基启为归算到重力基点处的重力仪数据，格；K0重力仪原始格值，10-5m/s2/格。步骤二、基于补偿后的观测重力值与绝对重力值的对应关系求取格值，并对船载重力仪数据进行补偿，具体包括：(3)按照步骤(1)和步骤(本文档来自技高网...

【技术保护点】
1.一种应用于船载重力仪的格值校正方法，其特征在于，包括以下步骤：/n步骤S1、以船载重力仪往返启航/返航港且中途停靠补给港的重力数据为基础数学模型，建立观测重力值的误差补偿校正模型，以对观测重力值进行误差补偿；/n步骤S2、基于补偿后的观测重力值与绝对重力值的对应关系求取格值，对船载重力仪数据进行补偿。/n

【技术特征摘要】
1.一种应用于船载重力仪的格值校正方法，其特征在于，包括以下步骤：
步骤S1、以船载重力仪往返启航/返航港且中途停靠补给港的重力数据为基础数学模型，建立观测重力值的误差补偿校正模型，以对观测重力值进行误差补偿；
步骤S2、基于补偿后的观测重力值与绝对重力值的对应关系求取格值，对船载重力仪数据进行补偿。


2.根据权利要求1所述的应用于船载重力仪的格值校正方法，其特征在于：所述步骤S1中对观测重力值进行误差补偿时具体包括以下步骤：
步骤S11、分析并分别获得启航码头、补给码头和返航码头的重力观测点与对应的各码头基点由于纬度变化导致的正常重力场值变化；
步骤S12、将各码头重力观测数据归算至重力基点高程面，得到归算后的各码头的重力基点处的重力仪数据。


3.根据权利要求2所述的应用于船载重力仪的格值校正方法，其特征在于：所述步骤S11具体通过以下方式获得正常重力场值变化：



gL＝0.000814*sin2L*dL
其中：dL为船载重力仪与基点的南北向距离；a为码头与正北方向夹角，所述码头包括启航码头、补给码头和返航码头；gL为正常重力场变化；L为重力观测点纬度。


4.根据权利要求3所述的应用于船载重力仪的格值校正方法，其特征在于：所述步骤S12中，归算后的启航码头的重力基点处的重力仪数据采用以下公式进行计算：



G基启＝G启+0.30...

【专利技术属性】
技术研发人员：马龙，郑彦鹏，
申请(专利权)人：自然资源部第一海洋研究所，
类型：发明
国别省市：山东;37