激光测距对静态重力场恢复精度的影响分析
激光测距对静态重力场恢复精度的影响分析
作者:朱永超;万晓云;于锦海;梁磊;张国庆;
高精度的地球重力场模型对经济、军事以及地球科学研究等具有十分重要的作用。为此,国外发射了多颗具有重力场探测效能的地球观测卫星,如CHAMP、GRACE、GOCE以及海洋测高卫星等。其中,CHAMP、GRACE和GOCE统称为重力卫星,GRACE卫星相对于CHAMP和GOCE卫星有更加广泛的应用前景,如海洋、南极冰盖、陆地水等地球科学方面的应用。 GRACE卫星已于2017年10月坠落,其后续卫星GRACE Follow-on已于2018年5月发射,近期Earth system mass transport mission(e.motion)和Next Generation satellite Gravimetry Mission(NGGM)卫星计划正在设计研究阶段。该类卫星的基本原理是通过高低跟踪和低低跟踪距离测量来探测地球重力场所引起的卫星星间距离的变化,并以此来反演地球重力场,其中高低跟踪测量主要通过星载GPS接收机来完成,低低跟踪则主要依靠微波测距完成。对于GRACE卫星而言,卫星低低跟踪数据的精度对重力场模型的精度起着至关重要的作用,因此在新的卫星计划中,除了微波测距外,还将搭载高精度的激光测距来进行更高精度的测距实验。由此,便产生这样的问题:激光测距对重力场模型恢复精度的提升会带来多大帮助,特别是当其他载荷的精度没有提升时。
综合以上很多结论是各类新技术的综合效应,同时不同文献讨论的对象不同,有的关注的是时变重力场,而有的结论则主要针对的是静态重力场。时变重力场强调实时性,一般需要以较短的时间为单位来解算重力场,如以月为单位,这就对各类时变背景场模型的精度提出了更高的要求。例如:当前所有的GRACE时变重力场模型均存在南北条带误差问题,而静态重力场则不存在,这表明静态重力场由于能够使用更多的数据,因此能够降低对各类背景场模型的精度要求。另外一点也需要注意的是定轨定速数据的权重问题,本文认为二者与星间距变率数据和加速度计数据的权重应该有所区别,对于不同的解算方法也会有差别。针对以上问题,本文限制了相应的前提,即:(1)主要讨论地球静态重力场模型的精度,而不考虑重力场解算过程中的时变误差;(2)主要评估采用激光测距后,相对已有的GRACE卫星而言,其余载荷和轨道均不变的情况,重力场模型的精度提升有多大。在以上约束条件下,本文主要通过高精度的数值反演方法讨论本文开头所提出的问题。
2 基本原理
2.1 动力学法原理
GRACE卫星反演重力场的方法主要有Kaula线性摄动法、动力学法、基线法、能量守恒法、加速度法和短弧长法等。从理论上讲动力学方法来自于对卫星轨道方程组的线性化,理论上相对是严密的,精度相对其他方法更高,因此本研究采用动力学方法进行模拟计算。
动力学积分法的基本原理是基于牛顿运动方程,将卫星轨道数据线性化表达成状态转移矩阵(即任意时刻的状态向量(卫星或速度)对初始时刻状态向量的偏导数矩阵)和参数敏感矩阵(即任意时刻的状态向量对未知参数的偏导数,如引力场位系数、加速度校正参数和力模型未知参数)的形式,来建立卫星跟踪数据与重力场位系数间的关系式,然后基于最小二乘准则解算引力场位系数。
2.2 星间变率的线性化推导
GRACE卫星主要测量两颗卫星之间的速度在两颗卫星的水平线上的分速度,即星间变率。星间变率测量的精度高于GPS测量位置差分得到的速度。相比于卫星精密轨道数据,星间变率的优势是对地球重力场中高阶次较为敏感,并且能够减少SST测量系统振荡器的不稳定性。但星间变率与重力场位系数等参数之间是非线性的,本文在于锦海等人(2017)推导的轨道扰动方程组的形式上进行星间变率公式线性化的推导。
2.3 卫星精密轨道数据与星间变率数据融合
由于两颗卫星在同一轨道上一前一后分别运行,若单独使用星间变率数据,则势必会造成局部法方程的欠定问题。因此必须要将卫星精密轨道数据与星间变率数据一起参与计算才可解算正确模型。一方面有卫星精密轨道数据提供卫星位置信息来降低局部法方程欠定的问题。另一方面轨道观测值对卫星重力场的低阶项较为敏感,能够提供一定精度的重力场信息。由于轨道的相对误差比较大,因此会将误差传递到法方程中。因此必须加相应的权比降低其对方程的影响,综上所述,必须融合卫星精密轨道数据及星间变率观测值来反演GRACE卫星的重力场,而实质是通过一定的比例关系将精密轨道方程与星间变率方程相加后进行求解,因此可以直接建立起两个法方程之间的关系。
3 数值计算结果分析
通过数值实验来讨论论文开头提出的问题。
(1)首先来验证有关程序的正确性;
(2)讨论不同情况下反演精度。
轨道数据,模拟GRACE的轨道参数与真实轨道参数相同,然后并进一步得到星间距变率数据,再利用上节所述的方法进行反演计算。为了减小计算量,模拟轨道使用EGM08模型前60阶完整阶次来当作真实场,数值积分方法我们使用Adam预测校正方法,所使用的软件为Intel FORTRAN 64位软件。在解算时不加入任何误差,得到了一组重力场参数,将真实的重力场与反演计算得到的重力场相减。
“EGM08”代表真实重力场EGM08模型的KUALA曲线,“EGM08-EGM96”代表真实重力场EGM08模型减去参考模型EGM96以后得到的曲线,即不同重力场模型之间的误差曲线,“Error-free”代表利用动力学方法反演得到的重力场与模型EGM08相减后的误差曲线,该曲线表明了本实验计算程序的正确性,能够满足本研究的计算精度要求。(2)解算流程与计算方法与(1)一致。区别在于首先将重力场的阶数增加到120阶。然后对观测数据加入几组误差,即GRACE卫星的设计精度指标,以及GRACE Follow On 的测量精度指标,主要是提高了星间变率的测量精度。
“EGM08”代表重力场模型EGM08的KUALA曲线,“GRACE”代表GRACE卫星的设计精度指标,“Follow-1、Follow-2、Follow-3、 Follow-4”分别代表星间变率计算精度提高到5.0×10–7 m/s、1.0×10–7 m/s、5.0×10–8 m/s、1.0×10–8 m/s。该计算结果表明当星间距变率提高到1.0×108 m/s时,相对于GRACE现有的精度指标,重力场恢复精度能够提高1个量级。保持GRACE现有的轨道和其他载荷精度,当星间变率精度依次从1.0×10–6 m/s提高到5.0×10–7 m/s、1.0×10–7 m/s、5.0×10–8 m/s、1.0× 10–8 m/s时,累积大地水准面误差则在120阶时从85.14 cm依次降低为33.09 cm、7.33 cm、3.70 cm、 3.59 cm。这表明当采用高精度的激光测距后,采用低低卫卫跟踪模式,地球静态重力场模型的探测精度有望比GRACE提高1个量级。
4 结 论
本文主要讨论了基于轨道扰动方程对星间变率进行了线性化的推导,并通过一定的权值将线性化的星间变率方程组与精密轨道方程组融合,通过实例验证了公式及程序的正确性。通过仿真计算验证了激光测距对静态重力场恢复精度提升的贡献。结果表明:当激光测距的精度达到50 nm/s时,轨道和其余指标与GRACE卫星相同时,静态重力场的恢复精度可提高约1个量级。但当激光测距精度提高至10 nm/s时,计算结果与精度为50 nm/s的计算结果无明显差别,这表明过高的星间变率测量精度相对于其他指标而言精度有冗余。因此,建议在其他测量技术精度未有提高的前提下,只需将激光测距的精度提高至50 nm/s即可。
值得说明的是:以上工作所讨论的对象为地球静态重力场,而非时变重力场。静态重力场反演可利用多年的数据,如1年甚至5年,因此部分时变误差可被大大削弱和平滑。而对于时变重力场,由于其解算以月甚至更短时间为单位,对时变环境误差的扣除提出了更高的要求(例如:固体潮、海潮、大气潮等),因此仅仅提高星间测距,对于时变重力场的解算不如静态重力场的帮助明显。时变重力场的解算除了提高测距精度外,还需发射多颗不同轨道倾角的卫星来提供观测数据。而无论是对于静态重力场,还是时变重力场的解算,通过激光测距来提高星间距离测量精度,是提高下一代重力卫星反演精度的关键因素之一。因此,及时开展基于激光测距技术的下一代重力卫星关键技术研究具有重要意义。
本文未能对时变重力场的精度提升问题进行系统的分析,下一步可对重力卫星测量的其他精度指标以及背景场模型(如海潮、大气潮等)误差对时变重力场恢复精度的影响进行研究,从而为大幅提升时变重力场的精度提供可行的卫星设计方案。
本文章转自爱学术(aixueshu.com),如有侵权,请联系删除