论文部分内容阅读
航空重力测量是一种快速地获取分布均匀、精度良好、大面积的重力场信息的有效方法。其中,捷联式航空重力仪以体积小、重量轻、成本低、可靠性高以及功耗小等优点成为航空重力测量研究的热点。除此之外,捷联式的重力仪还可以进行矢量测量,而重力矢量信息对于地球物理等诸多学科具有重要意义。本文基于国防科学技术大学自主研制的捷联式航空重力仪SGA-WZ,结合捷联式航空重力矢量测量的误差模型,在捷联式航空重力矢量测量误差分析、惯性器件误差对重力矢量测量的影响及补偿、重力扰动对重力矢量测量的影响及消除、飞机运动对重力矢量测量的影响以及重力矢量测量结果的误差分离等五个方面进行了研究。主要研究成果归纳如下:1.对捷联式航空重力矢量测量的原理进行了研究,讨论了基于位置更新的重力矢量测量算法。根据重力矢量测量的误差模型,文中设计了正交试验来分析各影响因素对航空重力矢量测量的影响程度。2.针对惯性器件误差对重力矢量测量的影响问题,提出了IPSE算法。在捷联式航空重力测量的数据处理中通常采用SINS/GNSS组合15状态卡尔曼滤波的方式来提高比力的测量精度,并将惯性器件的误差作为卡尔曼滤波器的滤波状态估计出来再进行补偿。然而,在水平匀速直线飞行条件下,15状态卡尔曼滤波是不完全可观的。本文对15状态卡尔曼滤波的可观性进行了分析,分析了可观的状态及状态的组合,指出了加速度计零偏及姿态误差不能单独地估计出来,只能估计出它们的线性组合。本文在研究SPS算法的基础上提出了IPSE算法,使得器件误差的逐步估计过程更加明晰,文中还对IPSE算法进行了仿真分析和实例验证。实例处理结果表明,重力矢量北向分量的精度从15.57m Gal提高到了4.12m Gal,东向分量的精度从17.36m Gal提高到了14.70m Gal。3.重力扰动的存在会对重力矢量测量带来影响。仿真分析了在导航解算中采用正常重力会对重力矢量的测量带来怎样的误差。针对重力扰动对重力矢量测量造成的误差,本文还分析了重力矢量估计的迭代算法,通过仿真分析指出了迭代算法的使用条件,证明了迭代算法的有效性,并利用迭代算法对航空重力实测数据进行了处理,验证结果表明,通过迭代算法能够使重力矢量的水平分量收敛到与真实值之间存在一个常值性的偏差。去除了直流分量之后,重力矢量北向分量的精度从4.12m Gal提高到了2.57m Gal,东向分量的精度从14.70m Gal提高到了3.06m Gal。4.以飞机为载体的航空重力矢量测量会受到飞机运动的影响。分析指出了长周期运动和慢螺旋运动有可能会对航空重力测量带来影响,其中慢螺旋运动只是引起测量过程中的定位误差,对重力测量真正带来影响的是长周期运动。文中还以两次实际的航空重力测量为例说明了飞机的实际运动情况,进一步确认了长周期运动对飞行轨迹以及比力测量的影响。并给出了飞机在飞行测量过程中动态性的衡量指标。根据重力矢量测量的误差频谱,文中对怎样减小飞机运动带来的重力矢量测量误差进行了一定的讨论,并对时间同步误差和姿态误差进行了估计和补偿,希望能在保证测量分辨率的同时提高测量的精度。最后,文中采用经验模分解的方式对重力扰动结果进行固有模态分解并去除其中的某些成分,最终实现了对飞机运动影响的抑制。5.在得到的重力矢量测量结果中还存在着误差,针对这一问题,本文对重力矢量测量结果的误差分离算法进行了研究。首先研究了波数相关滤波(WCF)方法,该方法针对重复测线或者间隔比较近的平行测线,认为这些测线之间相关性比较低的成分是误差,通过抑制这些相关性较低的成分以实现误差的分离。经过WCF处理之后,重力扰动北向分量的精度从2.49m Gal提高到了1.43m Gal,东向分量的精度从2.81m Gal提高到了1.80m Gal。其次研究了卡尔曼滤波逆向解算的应用,通过综合卡尔曼滤波正向解算和逆向解算的结果来消除测量结果中的误差。与单独采用正向解算结果而得到的重力扰动相比,综合正向解算和逆向解算结果之后,重力扰动北向分量的重复性精度从17.71m Gal提高到了4.87m Gal,东向分量的重复性精度从4.55m Gal提高到了4.02m Gal。最后研究了利用全球重力场模型来校正重力扰动中存在的系统性偏差和低频趋势。通过全球重力场模型计算出的重力矢量中虽然不含有重力场的高频信息,但其包含的长波信息可用于校正重力扰动中的系统性偏差和低频趋势。在山东某地试验中,利用全球重力场模型EGM2008校正之后,重力扰动北向分量的重复性从17.71m Gal改进到了2.74m Gal,东向分量的重复性从4.55m Gal改进到了3.20m Gal。在某海域试验中,利用EGM2008校正之后,重力扰动北向分量的重复性从50.24m Gal改进到了7.46m Gal,东向分量的重复性从48.11m Gal改进到了4.48m Gal。