短弧长积分法与平均加速度法比较分析.docVIP

短弧长积分法与平均加速度法比较分析 　　基金项目： 中央高校基本科研业务费专项资金资助项目（SWJTU12BR012） 　　作者简介： 黄强（1986-），男，博士研究生，研究方向为卫星重力测量，电话E-mail：406599573@ 　　文章编号： 0258-2724（2013）03-0455-06DOI： 10.3969/j.issn.0258-2724.2013.03.010 　　摘要： 　　为比较短弧长积分法和平均加速度法在卫星重力反演中的效果，用Fortran程序语言分别基于短弧长积分法和平均加速度法编写了2套计算程序，对GOCE（重力场和静态海洋环流探测）卫星精密科学轨道进行了反演计算，得到了2组90阶次的重力场模型Model_Arc（短弧长积分法）和Model_Acc（平均加速度法）.结果表明：Model_Acc大地水准面阶误差最大在80阶仅4 cm，而Model_Arc在90阶处达到10 cm；2个模型总体上都比用CHAMP和GRACE卫星的轨道数据反演的模型精度高，Model_Acc在80阶次以前比GO_CONS_GCF_2_TIM_R3精度高，但平均加速度法耗时比短弧长积分法长. 　　关键词： 　　短弧长积分法；平均加速度法；重力反演；GOCE（重力场和静态海洋环流探测） 　　中图分类号： P223文献标志码： A 　　利用轨道数据反演重力场模型的方法中，较常用的方法有能量守恒法、点加速度法、平均加速度法、短弧长积分法和动力学法等[1-4].其中，短弧长积分法和平均加速度法在反演过程中不需要速度向量，避免了数值微分计算后得到的低精度速度向量代入到反演计算中.另外，这2种方法得到的观测方程本身就是线性的，不需要先验重力场模型，也不需要迭代计算.因此，在轨道反演方法中，平均加速度法和短弧长积分法比其他方法具有一定的优势. 　　短弧长积分法和平均加速度法都基于牛顿运动方程，并且平均加速度法可以在短弧长积分法的基础上推导出来，因此平均加速度法又称为超短弧长积分法（每段弧长仅3个历元），但这2种方法的数值积分计算较复杂.游为等对短弧长积分法进行了研究，并提出了改进的短弧长积分法[5]；赵齐乐、柳响林等对平均加速度法进行了研究，取得了不错的效果[6]. 　　GOCE（gravity field and steady-state ocean circulation explorer）[6-14]卫星已于2009年3月17日发射升空，采集的数据包括卫星轨道数据和重???梯度数据，结合2种类型的观测数据，可以反演出兼顾中长波和中短波信息的高精度重力场模型，在数据释放以前武汉大学钟波等就开始了对GOCE卫星的模拟研究[15]. 　　由于短弧长积分法与平均加速度法存在一定程度的联系，因此分别基于短弧长积分法和平均加速度法编写了2套Fortran计算程序，计算了62 d的GOCE卫星轨道数据，反演了2组90阶次的重力场模型，对2种方法的计算结果进行了比较，并将计算模型同基于CHAMP（challenging mini-satellite payload）、GRACE（gravity recovery and climate experiment）、GOCE的一些高精度模型进行了比较分析，可供进一步联合GOCE卫星轨道和梯度数据反演重力场模型参考. 　　1 　　计算方法 　　1.1 　　短弧长积分法 　　短弧长积分法是基于牛顿运动方程，将卫星轨道表示成Fredholm积分方程形式的边界值问题，并采用积分方程的方法求解地球重力场模型.短弧长积分法的思想最初应用于轨道的确定，1969年Schneider 和 Reigber将该方法应用于地球重力场模型的确定，2005年Mayer首次将短弧长积分法应用于CHAMP卫星真实数据的处理，随后被应用到GRACE卫星数据确定重力场模型.这里略去繁琐的公式推导过程，直接给出短弧长积分法的观测方程： 　　值得注意的是，非保守力只在局部指北坐标系下具有一定的规律性，而方程又是在惯性系下联立的，因此，必须进行相关坐标系间的转换. 　　1.2 　　平均加速度法 　　加速度法分为点加速度法和平均加速度法2种，前者需要通过轨道内插和微分计算，后者则不需要，因此平均加速度法比点加速度法更严密. 　　平均加速度法可以在短弧长积分法的基础上推导出来.假设短弧长积分法每个弧段仅有3个历元，将其代入到短弧长积分法的方程中，经过对积分时间的转换，可以得到平均加速度法的观测方程： 　　平均加速度法的加速度由相邻3个历元差分计算，从而避免了点加速度法中数值微分计算和因数据缺失导致的大量历元信息不能使用的情况.另外，平均加速度法的积分计算仍采用移动窗口多项