CN103091722A - 基于载荷误差分析原理的卫星重力反演方法 - Google Patents
基于载荷误差分析原理的卫星重力反演方法 Download PDFInfo
- Publication number
- CN103091722A CN103091722A CN2013100241732A CN201310024173A CN103091722A CN 103091722 A CN103091722 A CN 103091722A CN 2013100241732 A CN2013100241732 A CN 2013100241732A CN 201310024173 A CN201310024173 A CN 201310024173A CN 103091722 A CN103091722 A CN 103091722A
- Authority
- CN
- China
- Prior art keywords
- delta
- centerdot
- sigma
- error
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明涉及一种精密测量地球重力场的方法,特别是一种基于载荷误差分析原理精确建立K波段测距仪的星间距离误差、GPS接收机的卫星轨道位置误差和轨道速度误差、以及星载加速度计的非保守力误差联合影响累计大地水准面精度的误差模型,进而精确和快速反演地球重力场的方法;该方法地球重力场反演精度高,保证计算精度的前提下有效提高反演速度,卫星重力反演过程简单,计算机性能要求低,卫星观测方程物理含义明确;基于载荷误差分析原理的卫星重力反演方法是计算高精度和高空间分辨率地球重力场的有效方法。
Description
一、技术领域
本发明涉及卫星重力学、空间大地测量学、航空航天等交叉技术领域,特别是涉及一种基于载荷误差分析原理精确和快速反演地球重力场的方法。
二、背景技术
21世纪是人类利用SST-HL/LL(Satellite-to-Satellite Tracking in theHigh-Low/Low-Low Mode)和SGG(Satellite Gravity Gradiometry)提升对数字地球认知能力的新纪元。地球重力场及其时变反映地球表层及内部物质的空间分布、运动和变化,同时决定着大地水准面的起伏和变化。因此,确定地球重力场的精细结构及其时变不仅是大地测量学、地球物理学、地震学、海洋学、空间科学、国防建设等的需求,同时也将为寻求资源、保护环境和预测灾害提供重要的信息资源。
GRACE(Gravity Recovery and Climate Experiment)双星采用近圆和近极地轨道设计,由美国宇航局(NASA)和德国航天局(DLR)共同研制开发。GRACE利用K波段测距仪高精度测量星间距离,利用高轨GPS(Global PositioningSystem)卫星对低轨双星精密跟踪定位,利用高精度SuperSTAR加速度计测量作用于双星的非保守力。GRACE系统既包含两组SST-HL,同时以差分原理测定两个低轨卫星之间的相互运动,因此得到的静态和动态全球重力场精度比CHAMP(Challenging Minisatellite Payload)至少高一个数量级,同时为将来GOCE(Gravity Field and Steady-State Ocean Circulation Explorer)卫星重力梯度测量奠定了坚实基础。
早在20世纪60年代,Baker首次提出了利用SST恢复地球重力场的重要思想。自此以后,国际大地测量学界的许多学者都积极投身于地球重力场恢复的方法与算法的理论研究和数值计算之中。在众多方法中,按照卫星观测方程的建立和求解的不同可分为解析法和数值法。解析法是指通过分析地球重力场和卫星观测数据的关系建立卫星观测方程模型,进而估计地球重力场的精度。解析法的优点是卫星观测方程物理含义明确,易于误差分析且可快速求解高阶地球重力场;缺点是由于在建立卫星观测方程模型时作了不同程度的近似,因此求解精度较低。数值法是指通过分析地球引力位系数和卫星观测数据的关系建立卫星观测方程,并通过最小二乘法拟合出地球引力位系数。数值法的优点是地球重力场求解精度较高;缺点是求解速度较慢且对计算机要求较高。不同于以前的技术,本发明基于载荷误差分析法建立了K波段测距仪的星间距离、GPS接收机的轨道位置和轨道速度、以及加速度计的非保守力误差联合影响累积大地水准面的误差模型,基于关键载荷精度指标的匹配关系论证了误差模型的可靠性,基于美国宇航局喷气推进实验室(NASA-JPL)公布的2009年的GRACE-Level-1B实测误差数据,有效和快速地反演了120阶GRACE地球重力场精度。
三、发明内容
本发明的目的是:基于载荷误差分析法较大程度优化地球重力场反演速度,而且进一步提高地球重力场反演精度。
为达到上述目的,本发明采用了如下技术方案:
基于载荷误差分析原理的卫星重力反演方法包含下列步骤:
步骤一:卫星关键载荷数据采集
1.1)通过星载K波段测距仪获取星间距离误差数据δρ12;
1.2)通过星载GPS接收机获取轨道位置误差数据δr和轨道速度误差数据
1.3)通过星载加速度计获取非保守力误差数据δf;
步骤二:关键载荷误差模型建立
2.1)K波段测距仪的星间距离误差模型
基于能量守恒定律,卫星观测方程可表示为
其中,表示卫星的瞬时速度,表示卫星的平均速度,GM表示地球质量M和万有引力常数G之积,r表示由卫星质心到地心之间的距离,表示由地球扰动位引起的速度变化;V=V0+T表示地球引力位,V0表示中心引力位,T表示扰动位;C表示能量积分常数;公式(1)可变形为
扰动位方差和速度变化方差的关系为
其中,δρ12表示K波段测距仪的星间距离误差,Δt表示采样间隔;
地球扰动位T(r,φ,λ)表示为
地球扰动位的方差表示为
其中,
基于球谐函数的正交性,公式(9)可化简为
其中,δClm,δSlm表示地球引力位系数精度;
大地水准面高的方差为
联合公式(4)、(7)和(12),可得累积大地水准面误差和星间距离误差之间的关系式
2.2)GPS接收机的轨道位置误差模型
在公式(15)两边同时微分可得
将公式(18)代入公式(13)可得累积大地水准面误差和轨道位置误差之间的关系式
2.3)GPS接收机的轨道速度误差模型
将公式(22)代入公式(13),可得累积大地水准面误差和轨道速度误差之间的关系式
2.4)加速度计的非保守力误差模型
由于公式(24)表示如下
将公式(25)代入公式(13)可得累积大地水准面误差和非保守力误差之间的关系式
2.5)关键载荷联合误差模型
联合公式(13)、(19)、(23)和(26),可得星间距离、轨道位置、轨道速度和非保守力误差联合影响累积大地水准面的误差模型
其中,
步骤三:地球重力场反演
第一,首先以0.5°×0.5°为网格分辨率,在地球表面的经度0°~360°和纬度-90°~90°范围内绘制网格;其次,在地球表面的轨迹点位置依次加入δη;最后,将分布于地球表面的δη平均归算于划分的网格点δη(φ,λ);
第二,将δη(φ,λ)按球谐函数展开为
δη在各阶处的方差表示为
本发明是基于载荷误差分析法有利于精确和快速反演地球重力场的特点而设计的,优点为:
1)地球重力场反演精度高;
2)保证计算精度的前提下有效提高反演速度;
3)卫星重力反演过程简单;
4)计算机性能要求低;
5)卫星观测方程物理含义明确。
四、附图说明
图1表示GRACE双星在轨飞行示意图。
图2表示GRACE-A卫星的全球轨道分布图。
图3表示GRACE卫星关键载荷误差在全球地表的分布(单位:μm/s)。
图4表示GRACE关键载荷匹配精度指标论证。
图5表示基于载荷误差分析法反演GRACE累积大地水准面误差。
五、具体实施方式
以下结合附图,以GRACE双星为例,对本发明的具体实施方式作进一步的说明。
基于载荷误差分析原理的卫星重力反演方法:
步骤一:卫星关键载荷数据采集
1.1)通过星载K波段测距仪获取星间距离误差数据δρ12;
1.3)通过星载加速度计获取非保守力误差数据δf。
步骤二:关键载荷误差模型建立
2.1)K波段测距仪的星间距离误差模型
基于能量守恒定律,卫星观测方程可表示为
其中,表示卫星的瞬时速度,表示卫星的平均速度,GM表示地球质量M和万有引力常数G之积,r表示由卫星质心到地心之间的距离,即r取平均值,r=Re(地球平均半径)+H(平均轨道高度),表示由扰动位引起的速度变化;V=V0+T表示地球引力位,V0表示中心引力位,T表示扰动位;C表示能量积分常数。公式(31)可变形为
扰动位方差和速度变化方差的关系为
其中,δρ12表示K波段测距仪的星间距离误差,Δt表示采样间隔。
地球扰动位T(r,φ,λ)表示为
地球扰动位的方差表示为
其中,
基于球谐函数的正交性,公式(39)可化简为
其中,δClm,δSlm表示地球引力位系数精度。
大地水准面高的方差为
联合公式(34)、(37)和(42),可得累积大地水准面误差和星间距离误差之间的关系式
2.2)GPS接收机的轨道位置误差模型
在公式(45)两边同时微分可得
基于公式(47)且星间距离误差δρ12和轨道位置误差δr的关系表示为
将公式(48)代入公式(43)可得累积大地水准面误差和轨道位置误差之间的关系式
2.3)GPS接收机的轨道速度误差模型
将公式(52)代入(43),可得累积大地水准面误差和轨道速度误差之间的关系式
2.4)加速度计的非保守力误差模型
由于双星受到的主要非保守力和星间速度近似同向,且非保守力通常表现为累积误差特性,据平方误差积分准则,星间速度误差和非保守力误差δf的关系表示为
将公式(55)代入公式(43)可得累积大地水准面误差和非保守力误差之间的关系式
2.5)关键载荷联合误差模型
联合公式(43)、(49)、(53)和(56),可得星间距离、轨道位置、轨道速度和非保守力误差联合影响累积大地水准面的误差模型
其中,
步骤三:地球重力场反演
第一,首先以0.5°×0.5°为网格分辨率,在地球表面的经度(0°~360°)和纬度(-90°~90°)范围内绘制网格;其次,按照GRACE卫星轨道(如图2所示)在地球表面的轨迹点位置依次加入δη;最后,如图3所示,将分布于地球表面的δη平均归算于划分的网格点δη(φ,λ)处,其中横坐标和纵坐标分别表示经度和纬度,颜色表示平均归算于网格点处的误差值δη(φ,λ)的大小(μm/s)。
第二,将δη(φ,λ)按球谐函数展开为
δη在各阶处的方差表示为
基于公式(59)计算将公式(60)代入公式(57),可有效和快速反演全球重力场精度。如图4所示,实线、圆圈线、十字线和虚线分别表示单独引入K波段测距仪的星间距离误差1×10-5m、GPS接收机的轨道位置误差3×10-2m和轨道速度误差3×10-5m/s、以及加速度计的非保守力误差3×10-10m/s2反演累积大地水准面误差。基于GRACE关键载荷精度指标的匹配关系,据图中4条曲线在各阶处的符合性可验证本发明基于载荷误差分析法建立的误差模型是可靠的。
如图5所示,虚线表示德国波兹坦地学研究中心(GFZ)公布的120阶EIGEN-GRACE02S全球重力场模型的实测精度,在120阶处反演累计大地水准面精度为18.938cm;实线表示基于关键载荷联合误差模型反演累计大地水准面的精度,在120阶处累计大地水准面精度为18.825cm。通过两条曲线在各阶处的符合性可知,载荷误差分析法是反演高精度和高空间分辨率全球重力场的有效方法之一。
以上具体实施方式仅为本发明的一种实施示例,其描述较为具体和详细,但不能因此而理解为对本发明专利范围的限制。其具体实施步骤顺序和模型参数可根据实际需要进行相应的调整。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。
Claims (1)
1.一种基于载荷误差分析原理的卫星重力反演方法,包含下列步骤:
步骤一:卫星关键载荷数据采集
1.1)通过星载K波段测距仪获取星间距离误差数据δρ12;
1.3)通过星载加速度计获取非保守力误差数据δf;
步骤二:关键载荷误差模型建立
2.1)K波段测距仪的星间距离误差模型
基于能量守恒定律,卫星观测方程可表示为
其中,表示卫星的瞬时速度,表示卫星的平均速度,GM表示地球质量M和万有引力常数G之积,r表示由卫星质心到地心之间的距离,表示由地球扰动位引起的速度变化;V=V0+T表示地球引力位,V0表示中心引力位,T表示扰动位;C表示能量积分常数;公式(1)可变形为
扰动位方差和速度变化方差的关系为
由于因此,公式(6)可表示为
其中,δρ12表示K波段测距仪的星间距离误差,Δt表示采样间隔;
地球扰动位T(r,φ,λ)表示为
地球扰动位的方差表示为
其中,
基于球谐函数的正交性,公式(9)可化简为
其中,δClm,δSlm表示地球引力位系数精度;
大地水准面高的方差为
联合公式(4)、(7)和(12),可得累积大地水准面误差和星间距离误差之间的关系式
2.2)GPS接收机的轨道位置误差模型
在公式(15)两边同时微分可得
将公式(18)代入公式(13)可得累积大地水准面误差和轨道位置误差之间的关系式
2.3)GPS接收机的轨道速度误差模型
将公式(22)代入公式(13),可得累积大地水准面误差和轨道速度误差之间的关系式
2.4)加速度计的非保守力误差模型
星间速度误差和非保守力误差δf的关系表示为
将公式(25)代入公式(13)可得累积大地水准面误差和非保守力误差之间的关系式
2.5)关键载荷联合误差模型
联合公式(13)、(19)、(23)和(26),可得星间距离、轨道位置、轨道速度和非保守力误差联合影响累积大地水准面的误差模型
其中,
步骤三:地球重力场反演
基于载荷误差分析法,利用星间距离误差数据δρ11、轨道位置误差数据δr和轨道速度误差数据以及非保守力误差数据δf反演累积大地水准面误差的过程如下:
第一,首先以0.5°×0.5°为网格分辨率,在地球表面的经度0°~360°和纬度-90°~90°范围内绘制网格;其次,在地球表面的轨迹点位置依次加入δη;最后,将分布于地球表面的δη平均归算于划分的网格点δη(φ,λ);
第二,将δη(φ,λ)按球谐函数展开为
δη在各阶处的方差表示为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310024173.2A CN103091722B (zh) | 2013-01-22 | 2013-01-22 | 基于载荷误差分析原理的卫星重力反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310024173.2A CN103091722B (zh) | 2013-01-22 | 2013-01-22 | 基于载荷误差分析原理的卫星重力反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103091722A true CN103091722A (zh) | 2013-05-08 |
CN103091722B CN103091722B (zh) | 2015-06-17 |
Family
ID=48204522
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310024173.2A Expired - Fee Related CN103091722B (zh) | 2013-01-22 | 2013-01-22 | 基于载荷误差分析原理的卫星重力反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103091722B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108020866A (zh) * | 2017-11-20 | 2018-05-11 | 中国空间技术研究院 | 一种星体重力场反演的方法和系统、以及处理器 |
CN108267792A (zh) * | 2018-04-13 | 2018-07-10 | 武汉大学 | 全球重力场模型反演方法 |
CN109557594A (zh) * | 2018-12-11 | 2019-04-02 | 中国人民解放军火箭军工程大学 | 基于重力异常时变的重力基准图时变修正方法及系统 |
CN111198402A (zh) * | 2020-01-15 | 2020-05-26 | 东华理工大学 | 基于轨道掩膜微分算子的地球重力场模型建模方法 |
CN111308570A (zh) * | 2020-03-04 | 2020-06-19 | 东华理工大学 | 基于载波相位微分速度构建全球重力场的方法 |
CN112729275A (zh) * | 2021-01-08 | 2021-04-30 | 中国船舶重工集团公司第七0七研究所 | 一种利用因子分析的卫星反演图重力适配区选取方法 |
CN112989589A (zh) * | 2021-03-05 | 2021-06-18 | 武汉大学 | 联合grace和gnss的局部地表质量变化反演方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498616A (zh) * | 2009-02-24 | 2009-08-05 | 航天东方红卫星有限公司 | 整星试验中基于应变反馈的载荷输入方法 |
CN102262248A (zh) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | 基于双星空间三维插值原理的卫星重力反演方法 |
CN102305949A (zh) * | 2011-06-30 | 2012-01-04 | 中国科学院测量与地球物理研究所 | 利用星间距离插值建立全球重力场模型的方法 |
CN102313905A (zh) * | 2011-07-18 | 2012-01-11 | 中国科学院测量与地球物理研究所 | 基于星间速度插值原理的卫星重力反演方法 |
CN102393535A (zh) * | 2011-07-20 | 2012-03-28 | 中国科学院测量与地球物理研究所 | 基于双星能量插值原理的卫星重力反演方法 |
-
2013
- 2013-01-22 CN CN201310024173.2A patent/CN103091722B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498616A (zh) * | 2009-02-24 | 2009-08-05 | 航天东方红卫星有限公司 | 整星试验中基于应变反馈的载荷输入方法 |
CN102262248A (zh) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | 基于双星空间三维插值原理的卫星重力反演方法 |
CN102305949A (zh) * | 2011-06-30 | 2012-01-04 | 中国科学院测量与地球物理研究所 | 利用星间距离插值建立全球重力场模型的方法 |
CN102313905A (zh) * | 2011-07-18 | 2012-01-11 | 中国科学院测量与地球物理研究所 | 基于星间速度插值原理的卫星重力反演方法 |
CN102393535A (zh) * | 2011-07-20 | 2012-03-28 | 中国科学院测量与地球物理研究所 | 基于双星能量插值原理的卫星重力反演方法 |
Non-Patent Citations (4)
Title |
---|
WEI ZHENG等: "Improving the accuracy of GRACE Earth’s gravitational field using the combination of different inclinations", 《PROGRESS IN NATURAL SCIENCE》, vol. 18, 31 December 2008 (2008-12-31), pages 555 - 561, XP022614994, DOI: 10.1016/j.pnsc.2007.11.017 * |
周旭华等: "用GRACE 卫星跟踪数据反演地球重力场", 《地球物理学报》, vol. 49, no. 3, 31 May 2006 (2006-05-31), pages 718 - 723 * |
郑伟等: "国际下一代卫星重力测量计划研究进展", 《大地测量与地球动力学》, vol. 32, no. 3, 30 June 2012 (2012-06-30), pages 152 - 159 * |
郑伟等: "基于卫-卫跟踪观测技术利用能量守恒法恢复地球重力场的数值模拟研究", 《地球物理学报》, vol. 49, no. 3, 31 May 2006 (2006-05-31), pages 712 - 717 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108020866A (zh) * | 2017-11-20 | 2018-05-11 | 中国空间技术研究院 | 一种星体重力场反演的方法和系统、以及处理器 |
CN108267792A (zh) * | 2018-04-13 | 2018-07-10 | 武汉大学 | 全球重力场模型反演方法 |
CN109557594A (zh) * | 2018-12-11 | 2019-04-02 | 中国人民解放军火箭军工程大学 | 基于重力异常时变的重力基准图时变修正方法及系统 |
CN111198402A (zh) * | 2020-01-15 | 2020-05-26 | 东华理工大学 | 基于轨道掩膜微分算子的地球重力场模型建模方法 |
CN111308570A (zh) * | 2020-03-04 | 2020-06-19 | 东华理工大学 | 基于载波相位微分速度构建全球重力场的方法 |
CN111308570B (zh) * | 2020-03-04 | 2022-09-06 | 东华理工大学 | 基于载波相位微分速度构建全球重力场的方法 |
CN112729275A (zh) * | 2021-01-08 | 2021-04-30 | 中国船舶重工集团公司第七0七研究所 | 一种利用因子分析的卫星反演图重力适配区选取方法 |
CN112989589A (zh) * | 2021-03-05 | 2021-06-18 | 武汉大学 | 联合grace和gnss的局部地表质量变化反演方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103091722B (zh) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103091722B (zh) | 基于载荷误差分析原理的卫星重力反演方法 | |
CN102262248B (zh) | 基于双星空间三维插值原理的卫星重力反演方法 | |
CN103076640B (zh) | 利用方差-协方差对角张量原理反演地球重力场的方法 | |
CN102313905B (zh) | 基于星间速度插值原理的地球重力反演方法 | |
CN102305949B (zh) | 利用星间距离插值建立全球重力场模型的方法 | |
CN102393535B (zh) | 基于双星能量插值原理的卫星重力反演方法 | |
CN104035138B (zh) | 一种全球及局部海洋扰动重力的精确快速计算方法 | |
CN103018783B (zh) | 重力卫星编队轨道稳定性优化设计和精密反演地球重力场方法 | |
CN102998713B (zh) | 基于功率谱半解析的卫星重力梯度反演方法 | |
CN103513294A (zh) | 一种低低星星跟踪卫星重力场测量性能解析计算方法 | |
CN101881619A (zh) | 基于姿态测量的船用捷联惯导与天文定位方法 | |
CN103076639B (zh) | 基于残余星间速度原理反演地球重力场的方法 | |
CN103163562B (zh) | 基于滤波原理的卫星重力梯度反演方法 | |
CN103093101B (zh) | 基于重力梯度误差模型原理的卫星重力反演方法 | |
CN108020866B (zh) | 一种星体重力场反演的方法和系统、以及处理器 | |
Ratnam et al. | Performance evaluation of selected ionospheric delay models during geomagnetic storm conditions in low-latitude region | |
CN103760537A (zh) | 基于卫星测高数据自身的潮汐校正方法 | |
CN103017787A (zh) | 适用于摇摆晃动基座的初始对准方法 | |
CN107270937A (zh) | 一种离线小波降噪快速初始对准方法 | |
CN103091721B (zh) | 利用不同轨道倾角卫星联合反演地球重力场的方法 | |
CN103575297A (zh) | 基于卫星导航接收机的gnss和mimu组合导航航向角估计方法 | |
CN103091723B (zh) | 降低重力卫星质心调整误差对地球重力场精度影响的方法 | |
Weintrit et al. | A novel approach to loxodrome (rhumb line), orthodrome (great circle) and geodesic line in ECDIS and navigation in general | |
CN103064128B (zh) | 基于星间距离误差模型的地球重力场恢复方法 | |
CN103913169A (zh) | 一种飞行器的捷联惯性/星光折射组合导航方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150617 Termination date: 20160122 |
|
EXPY | Termination of patent right or utility model |