平台式重力仪测量数据的卡尔曼滤波处理

平台式重力仪测量数据的卡尔曼滤波处理

蔡体菁;周薇;鞠玲玲

【摘要】根据三轴惯性平台海空重力仪的特点,在重力测量数据处理中,把重力异常作为状态量,建立了扩展卡尔曼滤波状态方程和观测方程,并阐述了三轴惯性平台海空重力仪重力测量数据处理步骤.依据给出的卡尔曼滤波方程,应用平滑卡尔曼滤波方法,对三轴惯性平台海空重力仪的海洋重力测量数据进行了处理.海空重力仪测量结果为:南北重复测线内符合精度为0.2 mGal,东西重复测线内符合精度为0.4 mGal,空间分辨率为650 m.

【期刊名称】《中国惯性技术学报》

【年(卷),期】2015(023)006

【总页数】4页(P718-720,726)

【关键词】海空重力仪;数据处理;卡尔曼滤波;重力异常

【作者】蔡体菁;周薇;鞠玲玲

【作者单位】东南大学仪器科学与工程学院,南京210096;东南大学仪器科学与工程学院,南京210096;东南大学仪器科学与工程学院,南京210096

【正文语种】中文

【中图分类】P631.1

自20世纪90年代以来，海空重力仪发展迅速，出现了一系列商用二轴和三轴稳定平台式海空重力仪[1-2]。二轴稳定平台式海空重力仪有：美国Micro-g LaCoste公司的L&RS、SII海空重力仪、TAGS航空重力仪，俄罗斯中央电气仪

表所的CHEKAN-AM海空重力仪。三轴稳定平台式海空重力仪有：加拿大Sander地球物理公司的AIRGrave航空重力仪，俄罗斯重力技术公司的GT-A航

空重力仪、GT-M海洋重力仪。近十几年来，我国海空重力仪发展也很迅猛[3-4]。中国船舶重工集团公司天津航海仪器研究所和中国科学院测量与地球物理研究所都分别研制出二轴稳定平台式重力仪，北京自动化设备研究所研制出三轴惯性平台海空重力仪。

目前国际上大多数稳定平台式海空重力仪配套的处理软件都是应用IIR滤波器和FIR滤波器来处理重力数据[5]，而GT-1A航空重力仪配套的重力处理软件不仅用

了FIR低通滤波，而且还用了卡尔曼滤波技术[6-7]。国内学者提出用小波方法处

理航空重力测量数据。东南大学在研究航空重力数据处理的FIR滤波、IIR滤波和

小波等方法[8]的基础上，提出用平滑卡尔曼滤波方法处理重力数据。本文针对国

产三轴惯性平台海空重力仪，给出重力测量数据的平滑卡尔曼滤波方法，以及处理重力测量数据的结果。

国产三轴惯性平台海空重力仪是由三轴惯性平台、重力敏感器、二次电源、控制电路、GPS接收机、计算机和显控存储装置等组成的。三轴惯性平台使用的惯性敏

感器是2个两自由度动力调谐陀螺和3个石英挠性加速度计；重力敏感器采用小

量程高精度的石英挠性加速度计；计算机进行导航计算；显控存储装置用来控制海空重力仪完成各个工作流程，存储和显示海空重力仪测量数据、导航参数和状态。海空重力仪工作流程包含准备和测量二个阶段，其流程见图1所示。

准备阶段主要是使惯性器件精度达到要求，完成重力仪的初始对准。具体完成的任务有：① 通电预热：重力仪内部所有分部件均启动工作，并对稳定平台进行温度

控制；② 参数自标定：利用稳定平台框架实现惯性器件的不同空间标定位置，标

定出陀螺、加速度计的关键参数；③ 抗扰动对准：调整稳定平台坐标系精确进入

测量工作状态。

测量阶段包括三个步骤：① 前校测量：完成准备后，在载体不移动的情况下，进

行一段时间的静态测试，自动记录重力仪和卫星数据；② 重力测量：平台系统工

作于动态导航状态，自动记录重力仪和卫星的数据，是重力测量的主体工作；③

后校测量：在完成测量后，载体回到出发点，重新进行一段时间的静态测试，自动记录重力仪和卫星数据。最后将所有数据存贮到硬盘设备中保留。

国产三轴惯性平台海空重力仪的重力传感器是被垂直安装在一个带舒拉调谐的三轴稳定平台上的。在航行中，校正回路可以将GPS的位置和速度信息提供给平台作

阻尼，平台导航坐标系为游动方位坐标系。它提供的信息为：三个加速度计的输出，载体的姿态、温度、平台参数，GPS移动站的星历、伪距、多普勒频移、载波相

位等原始数据。利用以上信息以及GPS基站的原始信息，能够计算得到测线上的

自由空间重力异常。自由空间重力异常Δg可以通过对惯性导航比力方程的变形得到，其方程为

式中：fup为加速度计的垂向分量，g0为正常重力场，为载体的垂向加速度，ωe 为地球旋转角速度，VE和VN为载体的东向和北向速度，φ为载体的纬度，Rλ和Rφ分别为参考椭球子卯酉圈和午圈曲率半径。

在工作时，海空重力仪的三轴惯性平台作舒拉或阻尼振荡，会引起平台倾斜，即平台存在一个失准角，方程(1)可改写为

式中：fz是重力传感器输出，fx、fy是水平加速度计输出，β、γ是平台失准角。在地球物理学中，重力异常可以用随机统计模型来描述。本文采用如下近似模型来描述重力异常：

式中：ρ为噪声强度，q为白噪声。

令方程(2)和方程(3)可以改写为卡尔曼滤波状态方程形式：

状态向量为 X=[h, W ,Δ g, U]T。

用差分 GPS的高度信息作为卡尔曼滤波的观测信息，观测方程为

式中：h为高度，H为测量矩阵，ν为测量噪声。式(4)和式(5)可以写成随机线性

离散系统形式：

式中：Xi是 ti时刻的状态向量；Φi+1,i是 ti时刻到 ti+1时刻的状态转移矩阵；

Bi是系统输入矩阵；ui是系统输入向量；Gi是系统噪声输入矩阵；wi是系统噪声向量；Zi是ti时刻的观测向量；Hi是观测矩阵；vi是观测噪声向量；wi和vi均

是零均值高斯白噪声且相互独立，其噪声强度分别为Qi和Ri。

应用平滑卡尔曼滤波方法求解重力异常。方程(6)的卡尔曼滤波为

式中：初始条件为是先验状态估计量，i是后验状态估计量，Ki是增益矩阵，Pi/i-

1是先验状态估计量的误差协方差阵， Pi是后验状态估计量的误差协方差阵。

在进行上述滤波的同时，需要存储i/i-1、Pi/i-1、i和 Pi，然后根据以下方程进行

滤波平滑：

式中：初始条件为是平滑后的状态估计量，是平滑后状态估计量的误差协方差阵。有关平滑卡尔曼滤波的理论与方法见参考文献[9]。

2015年6月国产三轴惯性平台海空重力仪在青岛进行了近海重力测量试验。海空重力仪放在近海科考01号船上，GPS基站设在码头岸上，船以13 n mile/h的平均速度在近海上航行，沿南北和东西两条航线做重复线测量。试验结束后，对GPS基站数据和海空重力仪测量数据进行事后重力数据处理。

三轴惯性平台海空重力仪重力测量数据后处理步骤如下：首先利用GPS基站和GPS移动站的原始数据，进行差分GPS的伪距、多普勒频移和载波相位计算方法，得到船航行的位置和速度；其次根据GPS提供的位置、速度信息和三轴惯性平台

海空重力仪的输出信息，用卡尔曼滤波技术估计出惯性稳定平台失准角、重力传感器误差和加速度计误差；然后对各种误差进行修正，再利用上面给出的平滑卡尔曼滤波，计算出测线上的自由空间重力异常。图2给出了三轴惯性平台海空重力仪

重力数据处理流程。