[发明专利]基于滤波原理的卫星重力梯度反演方法

说明书

技术领域

本发明涉及卫星重力梯度学、空间大地测量学、地球物理学、航空航天等交叉技术领域，特别是涉及一种通过滤波原理建立新型卫星重力梯度观测方程，进而精确和快速反演地球重力场的方法。

背景技术

地球重力场及其时变反映地球表层及内部物质的空间分布、运动和变化，同时决定着大地水准面的起伏和变化，因此确定地球重力场的精细结构及其时变不仅是大地测量学、地震学、海洋学、空间科学、国防建设等的需求，同时也将为全人类寻求资源、保护环境和预测灾害提供重要的信息资源。

GOCE（Gravity Field and Steady-State Ocean Circulation Explorer）重力梯度卫星由欧洲空间局（ESA）独立研制，已于2009年3月17日发射升空，主要用于地球中短波重力场的精密测量。GOCE采用近圆（轨道离心率0.001）、近极地（轨道倾角96.5°）和太阳同步轨道，经过20个月的飞行计划，轨道高度由250km降为240km。GOCE采用卫星跟踪卫星高低（SST-HL）和卫星重力梯度（SGG）模式的结合，除基于高轨道的GPS/GLONASS卫星对低轨道的GOCE卫星进行精密跟踪定位（1cm）外，利用定位于质心处的重力梯度仪（3×10^-12/s²）高精度测量卫星轨道高度处引力位的二阶导数，同时基于无阻尼离子微推进器补偿卫星受到的非保守力。早在20世纪80年代，国外便开始制定国际卫星重力梯度计划。由于地球重力场信号随卫星轨道高度的增加而急剧衰减(R_e/r)^l，基于分析卫星轨道运动仅适合于精密确定中长波地球重力场，而卫星重力梯度技术可直接测定地球引力位的二次微分，进而在一定程度上抑制了地球重力场信号的衰减效应，因此卫星重力梯度测量有利于高精度感测地球中短波重力场信号。

在利用卫星重力观测数据反演地球重力场的众多方法中，时域法的优点是直接对卫星观测数据进行处理，不需作任何近似，求解精度较高且能有效处理色噪声；缺点是随着卫星观测数据的增多，观测方程数量剧增，极大地增加了计算量。为了满足下一代卫星重力梯度测量计划中精密和快速解算高阶地球重力场的要求，本发明紧跟国际卫星重力梯度测量的热点和动态，利用新型滤波卫星重力梯度法，精确和快速反演了250阶GOCE地球重力场。

在现有技术中已有利用卫星重力梯度对角张量（V_xx,V_yy,V_zz），通过4个参考球面，在未考虑卫星重力梯度数据滤波情况下，恢复地球重力场的数值模拟方法。由于该方法未考虑卫星重力梯度交叉张量的贡献，而且未对卫星重力梯度观测数据进行有效滤波处理，因此地球重力场反演精度仍未达到预期要求。不同于已有技术，本发明基于卫星重力梯度全张量（V_xx,V_yy,V_zz,V_xy,V_xz,V_yz），通过5个参考球面，采用卫星重力梯度数据滤波处理技术，在不过多增加反演计算量的前提下，进一步提高了GOCE地球重力场的反演精度和计算速度。

发明内容

本发明的目的是：通过滤波原理，建立新型卫星重力梯度观测方程，进而精确和快速反演地球重力场。

为达到上述目的，本发明采用了如下技术方案：

一种基于滤波原理的卫星重力梯度反演方法，包括如下步骤：

步骤1：采集重力梯度卫星观测数据，其中通过重力梯度卫星的星载重力梯度仪采集卫星重力梯度全张量观测数据V_xx,V_yy,V_zz,V_xy,V_xz,V_yz，通过重力梯度卫星的星载GPS/GLONASS复合接收机采集卫星轨道位置观测数据r；

步骤2：建立卫星重力梯度观测模型，具体包括：