一种快速筛选地球准卫星小行星探测目标的方法技术

本发明专利技术公开的一种快速筛选地球准卫星小行星探测目标的方法，属于航空航天领域。本发明专利技术实现方法为：已知小行星当前时刻轨道状态的情况下，根据哈密顿平均近似方法，构建给定Kozai参数i

【技术实现步骤摘要】
一种快速筛选地球准卫星小行星探测目标的方法


[0001]本专利技术涉及一种快速筛选地球准卫星目标的方法，尤其涉及一种半解析的准卫星类型探测目标快速判定方法，属于航空航天领域。

技术介绍

[0002]作为研究太阳系起源的“活化石”，小天体保存着太阳系形成演化的原始信息，是目前国际深空探测的热点和行星科学研究的前沿。自上世纪90年代起，世界多个国家和组织进行了多项小天体探测活动，如飞越、着陆、采样返回等，其中近地小天体由于距离地球近、转移消耗少、撞击风险大，成为当前国内外小天体探测的重点。受到地球摄动、太阳热辐射等复杂空间环境影响，存在一类与地球共轨且绕地球运动的近地小天体，如2016 HO3、2004 GU9等。这类小天体可以在地球附近绕行几年到几百年不等，在此期间可以近似看作地球的天然卫星，因此被称为地球准卫星。与其他短暂飞越地球的近地小天体相比，这类小天体能够在相当长一段时间内与地球保持较近距离，因此其探测窗口时间长、转移轨道燃料消耗少，是实施小天体采样返回、小天体捕获等复杂深空任务的理想探测对象，例如2016 HO3小行星被选为我国未来小行星采样返回任务的潜在探测目标。随着近地天体观测计划的实施，近地小天体数目明显增加，从中快速筛选识别出地球准卫星，进而扩充地球准卫星探测目标库，对未来开展地球准卫星探测具有重要意义。
[0003]在先技术[1](参见：De la Fuente Marcos C,De la Fuente Marcos R.Asteroid (469219)2016 HO3,the smallest and closest Earth quasi
‑
satellite[J].Monthly Notices of the Royal Astronomical Society,2016,462(4):3441
‑
3456)在筛选地球准卫星时采用长期数值积分，通过分析其轨道运动参数随时间的变化曲线来判断其是否为地球准卫星，该方法准确性高但计算耗时较大，不适用于大批量小行星分析判定。
[0004]此外，在先技术[2](参见：Di Ruzza S,Pousse A,Alessi E M.On the co
‑
orbital asteroids in the solar system:medium
‑
term timescale analysis of the quasi
‑
coplanar objects[J].Icarus,2023,390:115330)，提出了一种利用平均化方法得到的平面可积模型，可以在不依赖数值积分的情况下判别近平面小行星的共轨类型，但该方法无法适用于轨道倾角较大的地球准卫星。

技术实现思路

[0005]本专利技术的主要目的是提供一种快速筛选地球准卫星小行星探测目标的方法，用于已知小行星当前时刻轨道状态的情况下，根据哈密顿平均近似方法，构建给定Kozai参数i
m
的地球共轨运动值的数据库，利用等值曲面计算得到临界共振角临界准卫星共轨状态小行星的哈密顿函数和临界碰撞地球小行星的哈密顿函数建立地球准卫星相空间分布条件。从小行星数据库中初步筛选符合条件地球共轨小行星候选对象，计算对应的小行星的哈密顿函数和共振角用准卫星的相空间分布条件判断小行星的动力学类型是否为地球准卫星，快速筛选地球准卫星小行星探测目标，提高地球准卫星探测目标的筛
选的效率。由于不依赖长期数值积分，本专利技术具有筛选速度快，效率高的特点。
[0006]本专利技术的目的是通过下述技术方案实现的。
[0007]本专利技术公开的一种快速筛选地球准卫星小行星探测目标的方法，基于日地圆型限制性三体问题，根据哈密顿平均近似方法，通过构建数据库的方式，计算地球准卫星相空间分布条件；从小行星数据库获取小行星数据，筛选满足判断条件的地球潜在共轨小行星；根据小行星当前轨道数据计算小行星的哈密顿函数和共振角判断在相空间的分布，利用准卫星的分布条件判断小行星的动力学类型是否为地球准卫星，快速筛选地球准卫星小行星探测目标。
[0008]本专利技术公开的一种快速筛选地球准卫星小行星探测目标的方法，包括如下步骤：
[0009]步骤一：将Kozai参数i
m
进行离散，对于每一个给定的i
m
，在可行范围内遍历小行星偏心率e、近地点俯角ω和根据哈密顿平均近似方法，计算不同e、ω和下对应的值，构建给定Kozai参数i
m
的地球共轨运动哈密顿函数值的数据库。
[0010]地球当前时刻日心黄道J2000惯性系下六个轨道根数：a
p
为轨道半长轴，e
p
为轨道偏心率，i
p
为轨道倾角，Ω
p
为黄道升交点经度，ω
p
为近日点幅角和平近点角M
p
时间匀速变化(周期为1年)。小行星当前时刻的六个轨道根数：a为半长轴，Ω为黄道升交点经度，e为轨道偏心率，i为轨道倾角，ω为近日点幅角，M为平近点角。小行星与地球间的共振角其中λ＝M+Ω+ω表示小行星的黄道经度，λ
p
＝M
p
+Ω
p
+ω
p
表示行星的黄道经度，小行星的Kozai参数i
m
表示为
[0011]在圆型限制性三体问题下，小行星的哈密顿函数的公式为其中可积项表达式为
[0012][0013]通过平均化近似方法，地球引力摄动项计算为
[0014][0015]其中μ为太阳和行星的质量系数，μ＝m
earth
/(m
sun
+m
earth
)，m
sun
为太阳质量系数，m
earth
为地球质量。r为在日心黄道J2000惯性系下小行星位置矢量，r
p
为在日心黄道J2000惯性系下地球位置矢量。
[0016]已知地球当前时刻日心黄道J2000惯性系下六个轨道根数a
p
＝1AU(AU为天文单位)，轨道偏心率e
p
＝0，轨道倾角i
p
＝0，升交点经度Ω
p
＝0，近日点幅角ω
p
＝0和平近点角M
p
＝λ
p
，小行星当前时刻的六个轨道根数为，半长轴a＝1AU，黄道升交点经度Ω＝0，共振角轨道偏心率e和近日点幅角ω为已知量，平近点角轨道倾角根据二体运动理论，日心黄道J2000惯性系下小行星位置矢量r由已知的小行星轨道根数转换得到；日心黄道J2000惯性系下地球位置矢量r
p
由已知的地球轨道根数转换得到。计算式和式，得到对于给定的Kozai参数i
m
、偏心率e、近地点俯角ω和共振角对应的的值。
[0017]将i
m
∈[2
°
,8本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种快速筛选地球准卫星小行星探测目标的方法，其特征在于：包括如下步骤，步骤一：将Kozai参数i
m
进行离散，对于每一个给定的i
m
，在可行范围内遍历小行星偏心率e、近地点俯角ω和根据哈密顿平均近似方法，计算不同e、ω和下对应的值，构建给定Kozai参数i
m
的地球共轨运动哈密顿函数值的数据库；步骤二：根据步骤一得到的给定Kozai参数i
m
的地球共轨运动哈密顿函数值的数据库，利用该数据库，获取给定对应的在相空间内的等值曲面，记录对应的等值曲面在范围内的个数为N；对于给定Kozai参数i
m
的数据库，选取等值曲面个数N为2的最小为临界值，此时在相空间内的两个等值曲面间的两个距离最近点记为和则计算临界即利用等值曲面计算得到临界共振角步骤三：对于每一个给定的Kozai参数i
m
重复步骤二，计算得到不同i
m
对应的值，进而通过样条插值方法得到，随i
m
的变化曲线；步骤四：对于步骤一离散的每一个给定的i
m
，利用临界准卫星共轨状态小行星的轨道根数和临界碰撞地球小行星的轨道根数，计算临界准卫星共轨状态小行星的哈密顿函数和临界碰撞地球小行星的哈密顿函数得到准卫星的相空间分布条件判断；步骤五：从小行星数据库中初步筛选地球共轨小行星候选对象，提取小行星日心黄道J2000惯性系下的六个轨道根数、轨道根数对应的时刻以及不确定度指标U；计算每一个小行星的Kozai参数i
m
，判断小行星的i
m
是否满足i
m
∈[2
°
,85
°
]，排除不满足条件的小行星数据；步骤六：利用步骤四提取的小行星日心黄道J2000惯性系下六个轨道根数数据，通过星历得到对应时刻的地球日心黄道J2000惯性系下六个轨道根数，利用相对轨道根数计算方法对小行星轨道根数进行处理，根据步骤一哈密顿平均近似方法计算对应的小行星的哈密顿函数和共振角步骤七：通过步骤三得到的随i
m
的变化曲线和通过步骤四得到的和随i
m
的变化曲线，利用样条插值的方法计算小行星对应i
m
的和结合步骤五得到的小行星的哈密顿函数和共振角结合相空间分布条件判断小行星是否为准卫星型小行星，即实现快速筛选地球准卫星小行星探测目标。2.如权利要求1所述的一种快速筛选地球准卫星小行星探测目标的方法，其特征在于：步骤一实现方法为，地球当前时刻日心黄道J2000惯性系下六个轨道根数：a
p
为轨道半长轴，e
p
为轨道偏心率，i
p
为轨道倾角，Ω
p
为黄道升交点经度，ω
p
为近日点幅角和平近点角M
p
时间匀速变化(周期为1年)；小行星当前时刻的六个轨道根数：a为半长轴，Ω为黄道升交点经度，e为轨道偏心率，i为轨道倾角，ω为近日点幅角，M为平近点角；小行星与地球间的共振角其中λ＝M+Ω+ω表示小行星的黄道经度，λ
p
＝M
p
+Ω
p
+ω
p
表示行星的黄道经度，小行星的Kozai参数i
m
表示为在圆型限制性三体问题下，小行星的哈密顿函数的公式为其中可积项
表达式为通过平均化近似方法，地球引力摄动项计算为其中μ为太阳和行星的质量系数，μ＝m
earth
/(m
sun
+m
earth
)，m
sun
为太阳质量系数，m
earth
为地球质量；r为在日心黄道J2000惯性系下小行星位置矢量，r
p
为在日心黄道J2000惯性系下地球位置矢量；已知地球当前时刻日心黄道J2000惯性系下六个轨道根数a
p
＝1AU(AU为天文单位)，轨道偏心率e
p
＝0，轨道倾角i
p
＝0，升交点经度Ω
p
＝0，近日点幅角ω
p
＝0和平近点角M
p
＝λ
p
，小行星当前时刻的六个轨道根数为，半长轴a＝1AU，黄道升交点经度Ω＝0，共振角轨道偏心率e和近日点幅角ω为已知量，平近点角轨道倾角根据二体运动理论，日心黄道J2000惯性系下小行星位置矢量r由已知的小行星轨道根数转换得到；日心黄道J2000惯性系下地球位置矢量r
p
由已知的地球轨道根数转换得到；计算式和式，得到对于给定的Kozai参数i
m
、偏心率e、近地点俯角ω和共振角对应的的值；将i
m
∈[2
°
,85
°
]进行离散，等间距取X个离散点；对于每一个给定的Kozai参数i
m
，对小行星轨道根数e、ω和在其可行范围内(e∈[0,sini
m
]、ω∈[0
°
,360
°
]和)等间距取Z个离散点；遍历离散的偏心率e、近地点俯角ω和共振角得到给定的i
m
的哈密顿函数数据库。3.如权利要求2所述的一种快速筛选地球准卫星小行星探测目标的方法，...

【专利技术属性】
技术研发人员：丁蓥，齐毅，乔栋，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：