CN111198402A - 基于轨道掩膜微分算子的地球重力场模型建模方法 - Google Patents
基于轨道掩膜微分算子的地球重力场模型建模方法 Download PDFInfo
- Publication number
- CN111198402A CN111198402A CN202010043914.1A CN202010043914A CN111198402A CN 111198402 A CN111198402 A CN 111198402A CN 202010043914 A CN202010043914 A CN 202010043914A CN 111198402 A CN111198402 A CN 111198402A
- Authority
- CN
- China
- Prior art keywords
- satellite
- representing
- orbit
- earth
- field model
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Navigation (AREA)
Abstract
本发明公开了基于轨道掩膜微分算子的地球重力场模型建模方法,涉及卫星重力探测技术领域,包括步骤:对轨道状态信息、星载加速度计和卫星姿态数据进行数据预处理;利用轨道掩膜微分算子计算卫星运动速度继而获得卫星动能和旋转位;通过卫星动能、旋转位、保守力位和非保守力位建立满足能量守恒定律的观测方程,在最小二乘准则下估算位系数改正值,再加上参考重力场模型位系数即可确定最终的地球重力场模型。相较于传统利用轨道牛顿多项式微分速度确定地球重力场模型的方法,本发明提出的方法有效抑制高频误差放大实现降噪,同时灵活性和计算效率高,获得高精度卫星速度观测值,进一步改进地球重力场模型精度水平。
Description
技术领域
本发明涉及卫星重力探测技术领域,特别是涉及基于轨道掩膜微分算子的地球重力场模型建模方法。
背景技术
地球重力场是地球最基本的一个物理场,反映了地球邻近空间及其内部物质迁移运动状态,因此确定高精度地球重力场模型是研究地球相关学科的热点问题。利用卫星重力测量技术改善重力场中长波信号,以及基于卫星跟踪卫星技术确定地球重力场模型的能量守恒法简洁高效被广泛采用。
能量守恒法是根据能量守恒定律建立卫星运动速度与地球重力场模型位系数间的线性函数关系,其关键问题是获取高精度的卫星运动速度值。计算卫星运动速度常用方法有经典的牛顿微分法,即通过牛顿插值的一次微分公式计算卫星速度,导致高频误差放大,降低速度观测数据的信噪比,进而影响重力场模型的建模精度。
发明内容
本发明实施例提供了基于轨道掩膜微分算子的地球重力场模型建模方法,可以解决现有技术中存在的问题。
本发明提供了基于轨道掩膜微分算子的地球重力场模型建模方法,包括以下步骤:
对重力专用卫星的轨道状态信息、星载加速度计数据和原始姿态数据进行数据预处理,获得连续精密轨道和非保守力加速度;
构建以计算时刻为中心的窗口函数,利用掩膜微分算子计算所述窗口函数中心时刻的卫星运动速度,并根据卫星运动速度计算卫星动能,同时根据卫星的位置和地球自转平均速度计算旋转位;
将所述连续精密轨道结合参考重力场模型计算正常重力位,通过背景模型确定潮汐位,根据所述正常重力位和潮汐位获得保守力位,对所述非保守力加速度通过沿轨积分得到非保守力位;在所述卫星动能中扣除旋转位、保守力位和非保守力位,根据能量守恒定律建立观测方程;针对所述观测方程利用最小二乘准则估算球谐位系数改正值,将所述球谐位系数改正值加上参考重力场模型的位系数确定最终的地球重力场模型。
本发明中的基于轨道掩膜微分算子的地球重力场模型建模方法,基于轨道掩膜微分算子确定卫星运动速度,继而计算卫星动能,再在卫星动能中扣除旋转位、保守力位和非保守力位,最后根据能量守恒定律确定地球重力场位系数。相较于传统利用轨道牛顿多项式微分速度确定地球重力场模型的方法,本发明提出的方法能够有效抑制高频误差放大实现降噪,同时灵活性和计算效率高,获得高精度卫星速度观测值,进一步改进地球重力场模型精度水平。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中方法的简化流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1,本发明提供了基于轨道掩膜微分算子的地球重力场模型建模方法,该方法包括以下步骤:
步骤1,数据预处理:对重力专用卫星的轨道状态信息、星载加速度计数据和原始姿态数据进行数据预处理,获得连续精密轨道和非保守力加速度。
具体地,对轨道状态信息进行内插和粗差探测处理后,得到连续精密轨道;
对连续精密轨道进行时空基准转换,将连续精密轨道转换至惯性系;
对星载加速度计数据和连续精密轨道进行时间同步,同时对原始姿态数据进行线性内插,计算惯性系下的非保守力加速度。
在进行时空基准转换时,采用经典的春分点法,即通过岁差、章动、地球自转和地球极旋四个矩阵实现。
步骤2,利用轨道掩膜微分算子计算卫星动能和旋转位:将连续精密轨道转换至惯性系后,利用掩膜微分算子计算窗口中心时刻的卫星运动速度观测值,进一步得到对应的卫星动能,同时结合卫星的位置和地球自转平均速度计算旋转位。
基于轨道掩膜微分算子确定卫星运动速度观测值,其原因在于有效抑制高频误差放大实现降噪,同时灵活性和计算效率高,获得高精度卫星速度观测值,进一步改进地球重力场模型精度水平。由步骤1获得的连续精密轨道采用轨道掩膜微分算子计算卫星运动速度,其推导过程如下。
重力专用卫星的轨道时间序列采用多项式逼近:
式中,x(t0+δm)N表示观测时刻卫星的轨道位置矢量;t0表示计算历元;δm表示计算历元与实际观测历元的差值;N表示多项式阶数;ck表示待定的多项式系数;M表示观测个数;Δt表示观测数据采样间隔;[·]表示向下取整符号。
为降低多项式逼近有限阶次引起的截断误差,观测历元往往关于中心时刻对称分布,要求多项式阶数偶数以及观测个数奇数。此时,对应的逐历元方程组式表述成:
上式的矩阵形式记作:
xM×1=BM×(N+1)c(N+1)×1
其中xM×1表示轨道位置向量;BM×(N+1)表示范德蒙特系数阵;c(N+1)×1表示待估的多项式系数。
需要说明的是,当观测个数恰好等于多项式阶数,称作插值过程;当观测个数远多于多项式阶数,称作拟合逼近过程。可见,多项式阶数和观测个数的不同组合形式使得算法灵活性高。
利用经典最小二乘平差法解算上式得到:
对轨道时间序列多项式微分可得:
仅考虑中心时刻,此时δ=0,则滤波核函数W中的元素γ=[0 1 0 … 0]T。
由此可得:
式中c1表示多项式系数的第二项;υ表示掩膜微分算子。由此可知,卫星运动速度观测值实质上仅与多项式系数c1有关,换句话说,卫星运动速度值仅是滤波核函数W第2行各元素与观测值的线性组合。
需要说明的是,当多项式阶数N、观测个数M、观测数据采样间隔Δt给定时,此时滤波核函数W为常数阵,使得掩膜微分算子υ各元素亦为常数,因此微分滤波系数仅需计算一次即可,有利于提高计算效率。
本步骤中,改进的解算策略采用轨道掩膜微分算子计算卫星运动速度:首先是在计算时刻左右两边各选取相同个数的观测值,构成以计算时刻为中心的窗口函数,然后对连续精密轨道进行最小二乘多项式拟合,得到轨道时间序列的多项式形式,通过线性组合滤波核函数与观测值计算窗口中心时刻的卫星运动速度;随后移动窗口,计算下一时刻,以此类推,组成平滑波谱,从而计算卫星运动速度的整个数据集。
步骤3,基于能量守恒定律估计地球重力场位系数:将连续精密轨道结合参考重力场模型计算正常重力位,同时通过背景模型获得潮汐位,将正常重力位和潮汐位求和得打保守力位,对非保守力加速度通过沿轨积分得到非保守力位。在卫星动能中扣除旋转位、保守力位和非保守力位后,根据能量守恒定律建立观测方程,最后,针对观测方程利用最小二乘准则估算球谐位系数改正值,将该球谐位系数改正值加上参考重力场模型的位系数即可确定最终的地球重力场模型。推导过程如下:
引力位及其梯度分量:
其中T表示扰动位,E0表示能量积分常数,U0表示正常重力位,表示卫星动能,表示卫星运动速度矢量,表示由地球自转引起的旋转位,ω表示地球自转角速度,Vt表示三体问题全部的直接和间接潮汐位,包括日月引力位、固体/海洋/大气、海洋/大气负荷等潮汐影响因素,Vc表示非保守力位,(r,θ,λ)分别表示地固球坐标下的地心向径、地心余纬和地心经度,GM和R分别表示地心引力常数和地球平均半径,l和m分别表示球谐级数展开式的阶和次,L表示确定地球重力场模型的最大阶数,和分别表示完全规格化的l阶m次球谐位系数改正值,表示完全规格化的l阶m次缔合勒让德函数。
需要说明的是,基于轨道掩膜微分算子确定地球重力场模型的方法是先由轨道掩膜微分算子计算卫星运动速度,继而获得卫星动能,再在卫星动能中扣除旋转位、保守力位和非保守力位后在惯性系下根据能量守恒定律建立观测方程,对该观测方程使用最小二乘准则估算球谐位系数改正值,最后在球谐位系数改正值的基础上加上参考重力场模型的位系数即可确定最终的地球重力场模型。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (5)
1.基于轨道掩膜微分算子的地球重力场模型建模方法,其特征在于,包括以下步骤:
对重力专用卫星的轨道状态信息、星载加速度计数据和原始姿态数据进行数据预处理,获得连续精密轨道和非保守力加速度;
构建以计算时刻为中心的窗口函数,利用掩膜微分算子计算所述窗口函数中心时刻的卫星运动速度,并根据卫星运动速度计算卫星动能,同时根据卫星的位置和地球自转平均速度计算旋转位;
将所述连续精密轨道结合参考重力场模型计算正常重力位,通过背景模型确定潮汐位,根据所述正常重力位和潮汐位获得保守力位,对所述非保守力加速度通过沿轨积分得到非保守力位;在所述卫星动能中扣除旋转位、保守力位和非保守力位,根据能量守恒定律建立观测方程;针对所述观测方程利用最小二乘准则估算球谐位系数改正值,将所述球谐位系数改正值加上参考重力场模型的位系数确定最终的地球重力场模型。
2.如权利要求1所述的基于轨道掩膜微分算子的地球重力场模型建模方法,其特征在于,对所述轨道状态信息、星载加速度计数据和原始姿态数据进行数据预处理的具体步骤为:
对所述轨道状态信息进行内插和粗差探测处理,得到所述连续精密轨道;
对所述连续精密轨道进行时空基准转换,将所述连续精密轨道转换至惯性系;
对所述星载加速度计数据和连续精密轨道进行时间同步,同时对所述原始姿态数据进行线性内插,计算惯性系下的非保守力加速度。
3.如权利要求1所述的基于轨道掩膜微分算子的地球重力场模型建模方法,其特征在于,计算窗口函数中心时刻的卫星运动速度的具体方法为:
在计算时刻左右两边各选取相同个数的观测值,构成以计算时刻为中心的窗口函数,然后对连续精密轨道进行最小二乘多项式拟合,得到轨道时间序列的多项式形式,引入掩膜微分算子后计算窗口函数中心时刻的卫星运动速度;随后移动窗口,计算下一时刻的卫星运动速度,以此类推,计算卫星运动速度的整个数据集。
4.如权利要求3所述的基于轨道掩膜微分算子的地球重力场模型建模方法,其特征在于,利用掩膜微分算子计算卫星速度的方法为:
重力专用卫星的轨道时间序列多项式表示为:
式中,x(t0+δm)N表示观测时刻卫星的轨道位置矢量,t0表示计算历元,δm表示计算历元与实际观测历元的差值,N表示多项式阶数,ck表示待定的多项式系数,M表示观测个数,Δt表示观测数据采样间隔,[·]表示向下取整符号;
与上式对应的逐历元方程组式为:
上式的矩阵形式记作:
xM×1=BM×(N+1)c(N+1)×1
其中xM×1表示轨道位置向量;BM×(N+1)表示范德蒙特系数阵;c(N+1)×1表示待估的多项式系数;
利用经典最小二乘平差法解算上式得到:
对所述轨道时间序列多项式微分可得:
仅考虑中心时刻,此时δ=0,则滤波核函数W中的元素γ=[0 1 0 … 0]T,由此可得:
式中c1表示多项式系数的第二项,υ表示掩膜微分算子;
5.如权利要求1所述的基于轨道掩膜微分算子的地球重力场模型建模方法,其特征在于,根据能量守恒定律建立的观测方程为:
利用下面引力位及其梯度分量表达式计算所述球谐位系数改正值:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010043914.1A CN111198402B (zh) | 2020-01-15 | 2020-01-15 | 基于轨道掩膜微分算子的地球重力场模型建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010043914.1A CN111198402B (zh) | 2020-01-15 | 2020-01-15 | 基于轨道掩膜微分算子的地球重力场模型建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111198402A true CN111198402A (zh) | 2020-05-26 |
CN111198402B CN111198402B (zh) | 2021-12-07 |
Family
ID=70744764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010043914.1A Active CN111198402B (zh) | 2020-01-15 | 2020-01-15 | 基于轨道掩膜微分算子的地球重力场模型建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111198402B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965124A (zh) * | 2021-02-08 | 2021-06-15 | 中国人民解放军92859部队 | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 |
CN117826271A (zh) * | 2024-01-03 | 2024-04-05 | 西安交通大学 | 一种地震短临前兆信息组网监测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011109175A2 (en) * | 2010-03-01 | 2011-09-09 | Geco Technology B.V. | Gravity measurements using seismic streamers |
CN102608669A (zh) * | 2012-02-22 | 2012-07-25 | 北京航空航天大学 | 一种具有移动与转动自由度的重力梯度柔性敏感结构 |
CN103091722A (zh) * | 2013-01-22 | 2013-05-08 | 中国科学院测量与地球物理研究所 | 基于载荷误差分析原理的卫星重力反演方法 |
CN107315200A (zh) * | 2017-05-03 | 2017-11-03 | 浙江大学 | 一种光力驱动的高精度绝对相对重力测量仪 |
CN108267792A (zh) * | 2018-04-13 | 2018-07-10 | 武汉大学 | 全球重力场模型反演方法 |
-
2020
- 2020-01-15 CN CN202010043914.1A patent/CN111198402B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011109175A2 (en) * | 2010-03-01 | 2011-09-09 | Geco Technology B.V. | Gravity measurements using seismic streamers |
CN102608669A (zh) * | 2012-02-22 | 2012-07-25 | 北京航空航天大学 | 一种具有移动与转动自由度的重力梯度柔性敏感结构 |
CN103091722A (zh) * | 2013-01-22 | 2013-05-08 | 中国科学院测量与地球物理研究所 | 基于载荷误差分析原理的卫星重力反演方法 |
CN107315200A (zh) * | 2017-05-03 | 2017-11-03 | 浙江大学 | 一种光力驱动的高精度绝对相对重力测量仪 |
CN108267792A (zh) * | 2018-04-13 | 2018-07-10 | 武汉大学 | 全球重力场模型反演方法 |
Non-Patent Citations (1)
Title |
---|
邹贤才: "GOCE卫星重力梯度校准与无阻尼控制效果分析", 《测绘学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965124A (zh) * | 2021-02-08 | 2021-06-15 | 中国人民解放军92859部队 | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 |
CN112965124B (zh) * | 2021-02-08 | 2022-10-11 | 中国人民解放军92859部队 | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 |
CN117826271A (zh) * | 2024-01-03 | 2024-04-05 | 西安交通大学 | 一种地震短临前兆信息组网监测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111198402B (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113819906B (zh) | 一种基于统计相似度量的组合导航鲁棒滤波方法 | |
CN111156987B (zh) | 基于残差补偿多速率ckf的惯性/天文组合导航方法 | |
CN108827310B (zh) | 一种船用星敏感器辅助陀螺仪在线标定方法 | |
CN110018507B (zh) | 一种基于星座间作差的组合精密单点定位方法及系统 | |
CN111198402B (zh) | 基于轨道掩膜微分算子的地球重力场模型建模方法 | |
CN111366984B (zh) | 一种基于重力卫星星间激光测距系统确定引力场模型的方法 | |
WO2017215026A1 (zh) | 一种基于高度约束的扩展卡尔曼滤波定位方法 | |
CN110567454A (zh) | 一种复杂环境下sins/dvl紧组合导航方法 | |
CN110567455B (zh) | 一种求积更新容积卡尔曼滤波的紧组合导航方法 | |
CN107607977B (zh) | 一种基于最小偏度单形采样的自适应ukf组合导航方法 | |
CN108731700B (zh) | 一种视觉惯性里程计中的加权欧拉预积分方法 | |
US5774831A (en) | System for improving average accuracy of signals from global positioning system by using a neural network to obtain signal correction values | |
CN112129300A (zh) | 位置间动力学约束的低轨卫星星载gnss精密定轨方法及系统 | |
CN111722295B (zh) | 一种水下捷联式重力测量数据处理方法 | |
CN108508463B (zh) | 基于Fourier-Hermite正交多项式扩展椭球集员滤波方法 | |
CN110988941A (zh) | 一种高精度实时绝对定轨方法 | |
CN115451952A (zh) | 一种故障检测及抗差自适应滤波的多系统组合导航方法和装置 | |
CN110554443A (zh) | 基于载波相位观测值和点加速度法确定地球重力场的方法 | |
CN115657097A (zh) | 一种基于轨道约束的leo几何法定轨模糊度快速重收敛方法 | |
CN113587926A (zh) | 一种航天器空间自主交会对接相对导航方法 | |
CN110703355B (zh) | 一种星载加速度计的校准方法及装置 | |
CN109813316B (zh) | 一种基于地形辅助的水下载体紧组合导航方法 | |
CN114167472A (zh) | Ins辅助gnss ppp精密动态导航定位方法及系统 | |
CN113008229A (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN111308570B (zh) | 基于载波相位微分速度构建全球重力场的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |