CN114114359A - 单星与地基设备联合的再入预报方法、装置和电子设备 - Google Patents
单星与地基设备联合的再入预报方法、装置和电子设备 Download PDFInfo
- Publication number
- CN114114359A CN114114359A CN202210098432.5A CN202210098432A CN114114359A CN 114114359 A CN114114359 A CN 114114359A CN 202210098432 A CN202210098432 A CN 202210098432A CN 114114359 A CN114114359 A CN 114114359A
- Authority
- CN
- China
- Prior art keywords
- reentry
- initial
- satellite
- target
- coefficient
- 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
- 238000013277 forecasting method Methods 0.000 title claims abstract description 20
- 238000012544 monitoring process Methods 0.000 claims abstract description 44
- 238000001514 detection method Methods 0.000 claims abstract description 28
- 230000010354 integration Effects 0.000 claims description 63
- 238000000034 method Methods 0.000 claims description 57
- 238000012937 correction Methods 0.000 claims description 25
- 238000004891 communication Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 3
- 238000013215 result calculation Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 14
- 230000009471 action Effects 0.000 description 4
- 229910052741 iridium Inorganic materials 0.000 description 4
- GKOZUEZYRPOHIO-UHFFFAOYSA-N iridium atom Chemical compound [Ir] GKOZUEZYRPOHIO-UHFFFAOYSA-N 0.000 description 4
- 238000013213 extrapolation Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000002679 ablation Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/26—Acquisition or tracking or demodulation of signals transmitted by the system involving a sensor measurement for aiding acquisition or tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种单星与地基设备联合的再入预报方法,包括:基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;利用初始轨道数据和初始弹道系数,预报再入目标的再入轨迹;根据再入轨迹预报多颗卫星对再入目标的可见性,并根据预报结果调整其中的可见卫星经过再入轨迹的空域时的侦照策略;基于一颗可见卫星对再入目标的监测数据,对初始弹道系数进行修正,得到修正后弹道系数;利用初始轨道数据和修正后弹道系数,获得再入目标的再入预报结果。本发明利用低轨巨型星座重访率高的优势,通过地基探测数据引导天基卫星对再入目标的再入轨迹进行监测来修正弹道系数,从而进行再入目标预报,可提高再入预报精度。
Description
技术领域
本发明属于航天测量与控制领域,具体涉及一种单星与地基设备联合的再入预报方法、装置和电子设备。
背景技术
绕地运行的低轨空间目标在完成任务后,受大气阻力作用,其轨道高度会逐渐降低,除少数航天器会在控制下安全返回底面,其他目标则会无控进入大气层,即发生再入事件。
大质量空间目标在再入过程中并不会被完全烧毁,仍有10~40%的残骸返回地球表面,这会对地表的生命群体、建筑设施、生态环境等产生很大威胁,这种目标也被称为危险再入目标。对危险再入目标保持连续监测,并迭代预报其再入事件及再入地点,是当前有效规避损害、为相关部门决策提供支撑的重要手段。再入预报的主要难点在于轨道确定和大气阻力建模,受限于当前半经验性大气模型的精度,当前短期精密再入预报的相对误差一般为15%左右,即提前10天预报,误差一般为1.5天;提前1天预报,误差为3.6小时。即使是陨落临圈,预报误差也将近15分钟,误差范围达6400余公里。因此,在地基设备布设受限、无法连续观测的情况下,对目标再入信息进行准确判别成为一个难题。
由于监测弧段距离目标最终再入越近,再入时间和落点的预报精度越高。但在不能全球布站的情况下,目前地基设备对再入目标的重访周期较长,部分目标的地基跟踪间隔甚至达16小时以上,因而无法实现持续监测,严重影响再入预报精度。其精度一般为轨道周期量级,落点区域覆盖整个地球,因而对最终再入情况无法进行准确判断。
因此,如何提出一种预报精度较高的判别再入事件方法,是本领域内一个亟待解决的问题。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种单星与地基设备联合的再入预报方法、装置和电子设备。本发明要解决的技术问题通过以下技术方案实现:
第一方面,本发明实施例提供了一种单星与地基设备联合的再入预报方法,包括:
基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;
利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;其中,所述再入轨迹包括多个轨道数据,预报的截止高度为预设高度;
根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略;
基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数;
利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果;其中,所述再入预报结果包括再入时刻、再入窗口和落点经纬度。
在本发明的一个实施例中,所述根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略,包括:
根据所述预报结果中所述多颗卫星对所述再入轨迹的可见情况,确定多个可见卫星;其中,所述可见情况包括卫星对所述再入轨迹的可见时段和可见时长;
根据所述可见情况,以及可见卫星和所述再入轨迹的相对位置关系,调整所述多个可见卫星的侦照角度以监测所述再入目标的再入轨迹。
在本发明的一个实施例中,所述基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数,包括:
利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数。
在本发明的一个实施例中,所述利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数,包括:
S001,解析所述监测数据,并将其转换成点迹序列;其中,所述点迹序列包括多个点数据,每个点数据包括观测时刻、星下点经度和星下点纬度;
S002,设置积分初值为所述初始轨道数据和所述初始弹道系数,积分步长为预设的第一积分步长;
S003,使用数值法进行轨道积分;
S004,判断积分结束时刻与所述一颗可见卫星的首点观测时刻之差是否小于所述第一积分步长;若否,返回S003;若是,执行S005;
S005,将积分步长减小为预设的第二积分步长;
S006,使用数值法进行轨道积分;
S007,将各次积分结束时刻得到的所述再入目标的状态矢量转换为大地坐标;
S008,针对每个大地坐标,判断该大地坐标与所述点迹序列中某点的位置差是否小于或等于预设位置差阈值;若否,返回S006;若是,执行S009;
S009,利用计算所述位置差时的大地坐标和所述点迹序列中某点分别对应的时刻,以及所述首点观测时刻,计算得到经过同一区域元时观测时刻与预报时刻的时间差;
S010,判断积分结束时刻是否大于所述点迹序列中最后的观测时刻;若否,返回S006;若是,执行S011;
S011,利用所述初始弹道系数和得到的多个时间差,计算弹道系数的修正因子;
S012,判断所述修正因子与所述初始弹道系数之商的绝对值是否小于预设的收敛阈值;若否,返回S002;若是,执行S013;
S013,将所述初始弹道系数与所述修正因子求和得到修正后弹道系数。
在本发明的一个实施例中,所述针对每个大地坐标,判断该大地坐标与所述点迹序列中某点的位置差是否小于或等于预设位置差阈值,包括:
其中,所述点迹序列表示为,其中、和分别为所述点迹序列中第个点数据的观测时刻、星下点经度和星下点纬度;大地坐标表示为,其中、和分别为第个大地坐标中的预报时刻、星下点经度和星下点纬度;表示预设位置差阈值。
在本发明的一个实施例中,所述时间差的计算公式,包括:
在本发明的一个实施例中,所述修正因子的计算公式,包括:
在本发明的一个实施例中,所述第一积分步长为10秒,所述第二积分步长为0.1秒。
第二方面,本发明实施例提供了一种单星与地基设备联合的再入预报装置,包括:
轨道数据和弹道系数初始获取模块,用于基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;
再入轨迹预报模块,用于利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;其中,所述再入轨迹包括多个轨道数据,预报的截止高度为预设高度;
卫星侦照策略调整模块,用于根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略;
弹道系数修正模块,用于基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数;
再入预报结果计算模块,用于利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果;其中,所述再入预报结果包括再入时刻、再入窗口和落点经纬度。
第三方面,本发明实施例提供了一种电子设备,包括处理器和存储器,其中,
所述存储器,用于存放计算机程序;
所述处理器,用于执行所述存储器上所存放的程序时,实现本发明实施例所提供的单星与地基设备联合的再入预报方法的步骤。
本发明实施例提供了一种联合利用单颗天基监测卫星与地基探测设备判别再入事件的方案。首先基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;其次利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;接下来根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略;然后基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数;最后利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果。本发明实施例利用低轨巨型星座重访率高的优势,通过地基探测数据引导天基卫星对再入目标的再入轨迹进行监测,进而修正弹道系数,基于修正后弹道系数进行再入目标预报,可提高再入时间和落点的预报精度、实现再入事件的准确判别。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1为本发明实施例提供的一种单星与地基设备联合的再入预报方法的流程示意图;
图2为本发明实施例提供的利用一颗可见卫星对再入目标的监测数据,采用区域元时间差的方法,修正初始弹道系数过程的流程示意图;
图3为本发明实施例所提供的一种单星与地基设备联合的再入预报装置的结构示意图;
图4为本发明实施例所提供的一种电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
现有技术中由于地基探测设备布设受限、监测间隔过长会导致难以准确判别大型空间目标再入事件,而大型再入目标在再入大气层的过程中,会与稠密大气发生剧烈摩擦,导致温度逐渐上升,最终烧蚀解体,其高温和光亮效应可被相关星载敏感器捕获。本发明实施例的发明人研究发现,由于低轨巨型星座重访率高,而大型再入目标的烧蚀过程往往持续近十分钟,星座组网联动则有极大概率发现目标再入过程,进而有望缩短预报时长、提高再入时间和落点的预报精度,实现再入事件的准确判别。
因此,本发明实施例提出了一种单星与地基设备联合的再入预报方法、装置和电子设备,通过联合利用单颗天基监测卫星与地基探测设备判别再入事件,以提高再入事件的预报精度。
第一方面,本发明实施例提供了一种单星与地基设备联合的再入预报方法,可以应用于地基设备或者除地基设备和卫星设备之外的其余设备中。请参见图1,图1为本发明实施例提供的一种单星与地基设备联合的再入预报方法的流程示意图。该方法包括以下步骤:
S1,基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数。
S2,利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹。
其中,所述再入目标的再入轨迹的预报方法属于现有技术,在此不做详细说明。
预报的截止高度为预设高度,比如可以为0km等。可选的一种实施方式中,所述预设高度可以为10km。
S3,根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略。
可选的一种实施方式中,该步骤包括:
S31,根据所述预报结果中所述多颗卫星对所述再入轨迹的可见情况,确定多个可见卫星。
其中,所述可见情况包括卫星对所述再入轨迹的可见时段和可见时长。
S32,根据所述可见情况,以及可见卫星和所述再入轨迹的相对位置关系,调整所述多个可见卫星的侦照角度以监测所述再入目标的再入轨迹。
其中,卫星可见性预报属于现有技术,具体方法在此不做赘述。该步骤的目的是通过卫星可见性预报获取多颗卫星对所述再入轨迹的可见情况,从中确定能够监测到所述再入目标的再入轨迹的部分卫星作为可见卫星,再调整所有可见卫星的光学载荷的侦照角度,使其可监视再入目标的再入轨迹,以此保证至少有一颗可见卫星能够获取到针对所述再入目标的监测数据。
S4,基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数。
由于实际的卫星数量较多,本发明实施例中能够监测到所述再入目标的再入轨迹的可见卫星有多个,仅从中获取一颗可见卫星的监测数据即可。
可选的一种实施方式中,该步骤包括:
利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数。
为了布局清晰,关于所述区域元时间差的方法在后文中予以说明。当然,任何能够利用所述一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正的方法都属于本发明实施例的保护范围。
S5,利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果。
具体的,利用所述初始轨道数据和所述修正后弹道系数,同时考虑地球非球形引力、大气阻力等摄动力因素,使用数值法进行轨道外推,直到轨道高度低于所述预设高度,可以得到所述再入目标的再入预报结果。
其中,所述再入预报结果包括再入时刻、再入窗口和落点经纬度。
该步骤中使用数值法进行轨道外推得到所述再入目标的再入预报结果的过程属于现有技术,在此不做详细说明。
现有技术中,仅地基装备进行监测的情况下,跟踪间隔太长,无法持续监测,导致预报误差较大。仅天基装备进行监测的情况下,难以合理确定监测策略,无法实现再入目标的有效捕获,且单颗天基监测卫星的光学测量精度有限,需要配合专门的再入预报方法,才能达到提高预报精度的目的。
本发明实施例提供了一种联合利用单颗天基监测卫星与地基探测设备判别再入事件的方案。首先基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;其次利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;接下来根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略;然后基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数;最后利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果。本发明实施例利用低轨巨型星座重访率高的优势,通过地基探测数据引导天基卫星对再入目标的再入轨迹进行监测,进而修正弹道系数,基于修正后弹道系数进行再入目标预报,可提高再入时间和落点的预报精度、实现再入事件的准确判别。
请参见图2,可选的一种实施方式中,所述利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数,包括以下步骤:
S001,解析所述监测数据,并将其转换成点迹序列。
该步骤采用现有技术实现,在此不做详细说明。
其中,所述点迹序列包括多个点数据,每个点数据包括观测时刻、星下点经度和星下点纬度。
S002,设置积分初值为所述初始轨道数据和所述初始弹道系数,积分步长为预设的第一积分步长。
S003,使用数值法进行轨道积分。
S004,判断积分结束时刻与所述一颗可见卫星的首点观测时刻之差是否小于所述第一积分步长。
S005,将积分步长减小为预设的第二积分步长。
所述第二积分步长可以选择小于所述第一积分步长的任意数值。可选的一种实施方式中,所述第二积分步长为0.1秒。
S006,使用数值法进行轨道积分。
该步骤是在减小积分步长之前的积分结果基础上,利用所述第二积分步长继续进行轨道积分,具体过程与前文类似。
S007,将每次积分结束时刻得到的所述再入目标的状态矢量转换为大地坐标。
本发明实施例将积分步长调整后,每一次积分结束时刻得到的均转换为大地坐标。具体转换过程请参见相关现有技术,在此不做赘述。得到的大地坐标表示为,其中、和分别为第个大地坐标中的预报时刻、星下点经度和星下点纬度;表示预设位置差阈值。
S008,针对每个大地坐标,判断该大地坐标与所述点迹序列中某点的位置差是否小于或等于预设位置差阈值。
其中,所述点迹序列表示为,其中、和分别为所述点迹序列中第个点数据的观测时刻、星下点经度和星下点纬度;大地坐标表示为,其中、和分别为第个大地坐标中的预报时刻、星下点经度和星下点纬度;表示预设位置差阈值。
根据当前判断结果执行后续操作,若否,返回S006;若是,执行S009。
S009,利用计算所述位置差时的大地坐标和所述点迹序列中某点分别对应的时刻,以及所述首点观测时刻,计算得到经过同一区域元时观测时刻与预报时刻的时间差。
具体的,所述时间差的计算公式,包括:
S010,判断积分结束时刻是否大于所述点迹序列中最后的观测时刻。
若否,返回S006。若是,执行S011。
S011,利用所述初始弹道系数和得到的多个时间差,计算弹道系数的修正因子。
具体的,所述修正因子的计算公式,包括:
S012,判断所述修正因子与所述初始弹道系数之商的绝对值是否小于预设的收敛阈值。
S013,将所述初始弹道系数与所述修正因子求和得到修正后弹道系数。
为了便于理解本发明实施例所提供方案的效果,以具体实验进行说明。
本发明实施例以北美防空联合司令部(NORAD)公开的编号为42821的空间目标作为再入目标进行举例,该目标为一个1U立方体卫星,根据Space-track网站的发布结果,该目标的再入时间为2019年5月21日。本发明实施例基于该再入目标的最后一组TLE(Two-Line Element,两行轨道根数)仿真生成地面跟踪的伪测量数据,作为S1中地基设备获取的探测数据,确定该再入目标的初始轨道数据,并解算出初始弹道系数,结果如表1所示。
执行S3,根据再入轨迹仿真铱星系列卫星对该再入目标的再入轨迹的可见情况,结果如表2所示,表2中仅列出了可见铱星、对应的可见时段和可见时长。根据铱星对再入轨迹的可见时段及相对位置关系,调整这些可见铱星光学载荷的侦照角度,使其可监视该再入目标的再入轨迹。
虽然调整了多颗可见卫星,但是最终哪颗卫星可以监测到再入事件的发光现象并不确定。仿真发现IRIDIUM 921卫星于17时45分发现目标,则利用该可见卫星对所述再入目标的监测数据执行S4。具体过程请参见前文所述,修正后的弹道系数为0.02194m2/kg。
可见,使用本发明实施例方法,能够修正最终的再入时间和落点位置,且使得再入窗口大大缩小。
综上,本发明实施例给出的该种仅单颗卫星监测到目标再入过程的情况下,天地基联合的再入判别方法,能够大大提高再入时间和落点预报精度,缩小残骸散布区域范围。
第二方面,相应于上述方法实施例,本发明实施例还提供了一种单星与地基设备联合的再入预报装置,如图3所示,该装置包括:
轨道数据和弹道系数初始获取模块301,用于基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数。
再入轨迹预报模块302,用于利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;其中,所述再入轨迹包括多个轨道数据,预报的截止高度为预设高度。
卫星侦照策略调整模块303,用于根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略。
弹道系数修正模块304,用于基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数。
再入预报结果计算模块305,用于利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果;其中,所述再入预报结果包括再入时刻、再入窗口和落点经纬度。
进一步的,所述卫星侦照策略调整模块303具体用于:
根据所述预报结果中所述多颗卫星对所述再入轨迹的可见情况,确定多个可见卫星;其中,所述可见情况包括卫星对所述再入轨迹的可见时段和可见时长。
根据所述可见情况,以及可见卫星和所述再入轨迹的相对位置关系,调整所述多个可见卫星的侦照角度以监测所述再入目标的再入轨迹。
进一步的,所述弹道系数修正模块304具体包括:
基于区域元时间差的弹道系数修正子模块,用于利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数。
进一步的,所述基于区域元时间差的弹道系数修正子模块,在利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数时,执行以下步骤:
S001,解析所述监测数据,并将其转换成点迹序列;其中,所述点迹序列包括多个点数据,每个点数据包括观测时刻、星下点经度和星下点纬度。
S002,设置积分初值为所述初始轨道数据和所述初始弹道系数,积分步长为预设的第一积分步长。
S003,使用数值法进行轨道积分。
S004,判断积分结束时刻与所述一颗可见卫星的首点观测时刻之差是否小于所述第一积分步长;若否,返回S003;若是,执行S005。
S005,将积分步长减小为预设的第二积分步长。
S006,使用数值法进行轨道积分。
S007,将各次积分结束时刻得到的所述再入目标的状态矢量转换为大地坐标。
S008,针对每个大地坐标,判断该大地坐标与所述点迹序列中某点的位置差是否小于或等于预设位置差阈值;若否,返回S006;若是,执行S009。
S009,利用计算所述位置差时的大地坐标和所述点迹序列中某点分别对应的时刻,以及所述首点观测时刻,计算得到经过同一区域元时观测时刻与预报时刻的时间差。
S010,判断积分结束时刻是否大于所述点迹序列中最后的观测时刻;若否,返回S006;若是,执行S011。
S011,利用所述初始弹道系数和得到的多个时间差,计算弹道系数的修正因子。
S012,判断所述修正因子与所述初始弹道系数之商的绝对值是否小于预设的收敛阈值;若否,返回S002;若是,执行S013。
S013,将所述初始弹道系数与所述修正因子求和得到修正后弹道系数。
进一步的,S008中,所述针对每个大地坐标,判断该大地坐标与所述点迹序列中某点的位置差是否小于或等于预设位置差阈值,包括:
其中,所述点迹序列表示为,其中、和分别为所述点迹序列中第个点数据的观测时刻、星下点经度和星下点纬度;大地坐标表示为,其中、和分别为第个大地坐标中的预报时刻、星下点经度和星下点纬度;表示预设位置差阈值。
进一步的,S009中,所述时间差的计算公式,包括:
进一步的,S011中,所述修正因子的计算公式,包括:
进一步的,所述第一积分步长为10秒,所述第二积分步长为0.1秒。
关于具体内容请参见第一方面所述的方法部分,在此不做赘述。
本发明实施例所提供的方案中,利用低轨巨型星座重访率高的优势,通过地基探测数据引导天基卫星对再入目标的再入轨迹进行监测,进而修正弹道系数,基于修正后弹道系数进行再入目标预报,能够提高再入时间和落点的预报精度、实现再入事件的准确判别。
第三方面,本发明实施例还提供了一种电子设备,可以为地基设备或者除地基设备和卫星设备之外的其余设备。如图4所示,包括处理器401、通信接口402、存储器403和通信总线404,其中,处理器401,通信接口402,存储器403通过通信总线404完成相互间的通信,
存储器403,用于存放计算机程序;
处理器401,用于执行存储器403上所存放的程序时,实现如第一方面所述的单星与地基设备联合的再入预报方法的步骤。
该电子设备可以为:台式计算机、便携式计算机、智能移动终端、服务器等。在此不作限定,任何可以实现本发明的电子设备,均属于本发明的保护范围。
上述电子设备提到的通信总线可以是外设部件互连标准(Peripheral ComponentInterconnect,PCI)总线或扩展工业标准结构(Extended Industry StandardArchitecture,EISA)总线等。该通信总线可以分为地址总线、数据总线、控制总线等。为便于表示,图中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。
通信接口用于上述电子设备与其他设备之间的通信。
存储器可以包括随机存取存储器(Random Access Memory,RAM),也可以包括非易失性存储器(Non-Volatile Memory,NVM),例如至少一个磁盘存储器。可选的,存储器还可以是至少一个位于远离前述处理器的存储装置。
上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(Network Processor,NP)等;还可以是数字信号处理器(Digital SignalProcessing,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
对于装置/电子设备而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
需要说明的是,本发明实施例的装置、电子设备分别是应用上述单星与地基设备联合的再入预报方法的装置、电子设备,则上述单星与地基设备联合的再入预报方法的所有实施例均适用于该装置、电子设备,且均能达到相同或相似的有益效果。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (10)
1.一种单星与地基设备联合的再入预报方法,其特征在于,包括:
基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;
利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;其中,所述再入轨迹包括多个轨道数据,预报的截止高度为预设高度;
根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略;
基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数;
利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果;其中,所述再入预报结果包括再入时刻、再入窗口和落点经纬度。
2.根据权利要求1所述的单星与地基设备联合的再入预报方法,其特征在于,所述根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略,包括:
根据所述预报结果中所述多颗卫星对所述再入轨迹的可见情况,确定多个可见卫星;其中,所述可见情况包括卫星对所述再入轨迹的可见时段和可见时长;
根据所述可见情况,以及可见卫星和所述再入轨迹的相对位置关系,调整所述多个可见卫星的侦照角度以监测所述再入目标的再入轨迹。
3.根据权利要求1所述的单星与地基设备联合的再入预报方法,其特征在于,所述基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数,包括:
利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数。
4.根据权利要求3所述的单星与地基设备联合的再入预报方法,其特征在于,所述利用一颗可见卫星对所述再入目标的监测数据,采用区域元时间差的方法,对所述初始弹道系数进行修正,得到修正后弹道系数,包括:
S001,解析所述监测数据,并将其转换成点迹序列;其中,所述点迹序列包括多个点数据,每个点数据包括观测时刻、星下点经度和星下点纬度;
S002,设置积分初值为所述初始轨道数据和所述初始弹道系数,积分步长为预设的第一积分步长;
S003,使用数值法进行轨道积分;
S004,判断积分结束时刻与所述一颗可见卫星的首点观测时刻之差是否小于所述第一积分步长;若否,返回S003;若是,执行S005;
S005,将积分步长减小为预设的第二积分步长;
S006,使用数值法进行轨道积分;
S007,将各次积分结束时刻得到的所述再入目标的状态矢量转换为大地坐标;
S008,针对每个大地坐标,判断该大地坐标与所述点迹序列中某点的位置差是否小于或等于预设位置差阈值;若否,返回S006;若是,执行S009;
S009,利用计算所述位置差时的大地坐标和所述点迹序列中某点分别对应的时刻,以及所述首点观测时刻,计算得到经过同一区域元时观测时刻与预报时刻的时间差;
S010,判断积分结束时刻是否大于所述点迹序列中最后的观测时刻;若否,返回S006;若是,执行S011;
S011,利用所述初始弹道系数和得到的多个时间差,计算弹道系数的修正因子;
S012,判断所述修正因子与所述初始弹道系数之商的绝对值是否小于预设的收敛阈值;若否,返回S002;若是,执行S013;
S013,将所述初始弹道系数与所述修正因子求和得到修正后弹道系数。
8.根据权利要求4所述的单星与地基设备联合的再入预报方法,其特征在于,所述第一积分步长为10秒,所述第二积分步长为0.1秒。
9.一种单星与地基设备联合的再入预报装置,其特征在于,包括:
轨道数据和弹道系数初始获取模块,用于基于地基设备获取的探测数据,确定再入目标的初始轨道数据,并解算出初始弹道系数;
再入轨迹预报模块,用于利用所述初始轨道数据和所述初始弹道系数,预报所述再入目标的再入轨迹;其中,所述再入轨迹包括多个轨道数据,预报的截止高度为预设高度;
卫星侦照策略调整模块,用于根据所述再入轨迹预报多颗卫星对所述再入目标的可见性,并根据预报结果调整其中的可见卫星经过所述再入轨迹的空域时的侦照策略;
弹道系数修正模块,用于基于一颗可见卫星对所述再入目标的监测数据,对所述初始弹道系数进行修正,得到修正后弹道系数;
再入预报结果计算模块,用于利用所述初始轨道数据和所述修正后弹道系数,获得所述再入目标的再入预报结果;其中,所述再入预报结果包括再入时刻、再入窗口和落点经纬度。
10.一种电子设备,其特征在于,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现权利要求1-8任一所述的方法步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210098432.5A CN114114359B (zh) | 2022-01-27 | 2022-01-27 | 单星与地基设备联合的再入预报方法、装置和电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210098432.5A CN114114359B (zh) | 2022-01-27 | 2022-01-27 | 单星与地基设备联合的再入预报方法、装置和电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114114359A true CN114114359A (zh) | 2022-03-01 |
CN114114359B CN114114359B (zh) | 2022-04-26 |
Family
ID=80361297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210098432.5A Active CN114114359B (zh) | 2022-01-27 | 2022-01-27 | 单星与地基设备联合的再入预报方法、装置和电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114114359B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002162195A (ja) * | 2000-11-21 | 2002-06-07 | Toshiba Corp | 飛翔体誘導装置 |
WO2016091728A1 (en) * | 2014-12-10 | 2016-06-16 | The Provost, Fellows, Foundation Scholars, And The Other Members Of Board, Of The College Of The Holy And Undivided Trinity Of Queen Elizabeth, Near Dublin | Computer modelling of atmospheric re-entry of objects such as satellites |
CN109323698A (zh) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | 一种空间目标陨落多模型跟踪引导技术 |
CN109992927A (zh) * | 2019-04-27 | 2019-07-09 | 中国人民解放军32035部队 | 稀疏数据情况下小椭圆目标的再入预报方法 |
CN110427727A (zh) * | 2019-08-14 | 2019-11-08 | 中国人民解放军32035部队 | 一种基于单组tle的空间目标轨道寿命快速分析方法 |
CN110471450A (zh) * | 2019-08-29 | 2019-11-19 | 大连理工大学 | 在高度速度剖面内直接规划再入轨迹的方法 |
CN110489879A (zh) * | 2019-08-22 | 2019-11-22 | 中国人民解放军32035部队 | 一种适用于空间环境扰动情况下的空间目标陨落预报方法 |
US20190377770A1 (en) * | 2018-06-01 | 2019-12-12 | Korea Astronomy And Space Science Institute | Method for re-entry prediction of uncontrolled artificial space object |
CN110598270A (zh) * | 2019-08-22 | 2019-12-20 | 中国人民解放军32035部队 | 一种基于编目根数序列的高精度空间目标陨落预报方法 |
CN113788166A (zh) * | 2021-09-16 | 2021-12-14 | 中国科学院国家天文台 | 基于空间物体轨道误差的差速拦截跟踪方法 |
-
2022
- 2022-01-27 CN CN202210098432.5A patent/CN114114359B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002162195A (ja) * | 2000-11-21 | 2002-06-07 | Toshiba Corp | 飛翔体誘導装置 |
WO2016091728A1 (en) * | 2014-12-10 | 2016-06-16 | The Provost, Fellows, Foundation Scholars, And The Other Members Of Board, Of The College Of The Holy And Undivided Trinity Of Queen Elizabeth, Near Dublin | Computer modelling of atmospheric re-entry of objects such as satellites |
US20190377770A1 (en) * | 2018-06-01 | 2019-12-12 | Korea Astronomy And Space Science Institute | Method for re-entry prediction of uncontrolled artificial space object |
CN109323698A (zh) * | 2018-12-03 | 2019-02-12 | 西安四方星途测控技术有限公司 | 一种空间目标陨落多模型跟踪引导技术 |
CN109992927A (zh) * | 2019-04-27 | 2019-07-09 | 中国人民解放军32035部队 | 稀疏数据情况下小椭圆目标的再入预报方法 |
CN110427727A (zh) * | 2019-08-14 | 2019-11-08 | 中国人民解放军32035部队 | 一种基于单组tle的空间目标轨道寿命快速分析方法 |
CN110489879A (zh) * | 2019-08-22 | 2019-11-22 | 中国人民解放军32035部队 | 一种适用于空间环境扰动情况下的空间目标陨落预报方法 |
CN110598270A (zh) * | 2019-08-22 | 2019-12-20 | 中国人民解放军32035部队 | 一种基于编目根数序列的高精度空间目标陨落预报方法 |
CN110471450A (zh) * | 2019-08-29 | 2019-11-19 | 大连理工大学 | 在高度速度剖面内直接规划再入轨迹的方法 |
CN113788166A (zh) * | 2021-09-16 | 2021-12-14 | 中国科学院国家天文台 | 基于空间物体轨道误差的差速拦截跟踪方法 |
Non-Patent Citations (10)
Title |
---|
A. M. DARYA ET AL.: "Analysis of Space Debris Re-Entry over the Arabian Peninsula (2004 to 2018)", 《2018 INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING AND INFORMATION SECURITY (ICSPIS)》 * |
EUN-JUNG CHOI ET AL.: "A Study on Re-entry Predictions of Uncontrolled Space Objects for Space Situational Awareness", 《JOURNAL OF ASTRONOMY AND SPACE SCIENCES》 * |
OKCHUL JUNG ET AL.: "Recurrent neural network model to predict re-entry trajectories of uncontrolled space objects", 《ADVANCES IN SPACE RESEARCH》 * |
SITTIPORN CHANNUMSIN ET AL.: "Analysis of coupled attitude and orbit dynamics for uncontrolled re-entry satellite", 《2019 7TH INTERNATIONAL ELECTRICAL ENGINEERING CONGRESS (IEECON)》 * |
刘劲宏: "基于TLE数据的空间目标大气再入预报", 《测绘学报》 * |
周琳等: "基于可观性分析的高精度空间目标跟踪方法", 《雷达科学与技术》 * |
张炜等: "基于TLE的弹道系数计算方法及应用分析", 《中国空间科学技术》 * |
张炜等: "基于大气模型误差特性的"天宫一号"再入预报", 《航天返回与遥感》 * |
张炜等: "稀疏数据情况下小椭圆轨道目标再入预报方法", 《空间碎片研究》 * |
李君龙等: "临近空间目标跟踪与预报技术研究", 《现代防御技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114114359B (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105224737B (zh) | 一种空间目标轨道改进初值修正方法 | |
CN111551183B (zh) | 基于天基光学观测数据的geo目标多点择优短弧定轨方法 | |
Bennett et al. | Accurate orbit predictions for debris orbit manoeuvre using ground-based lasers | |
Sang et al. | Experimental results of debris orbit predictions using sparse tracking data from Mt. Stromlo | |
US20210356275A1 (en) | Method of satellite precise orbit determination using parallactic refraction scale factor estimation | |
CN112857306B (zh) | 一种视频卫星任意视向点的连续太阳高度角确定方法 | |
Pham et al. | An autonomous star recognition algorithm with optimized database | |
CN103971167B (zh) | 预测光伏发电设备的发电功率的方法和系统 | |
CN112435301A (zh) | 一种基于恒星轨迹的遥感相机在轨几何定标方法 | |
CN111912430A (zh) | 高轨光学卫星的在轨几何定标方法、装置、设备及介质 | |
Song et al. | Autonomous estimation of Allan variance coefficients of onboard fiber optic gyro | |
CN109855652B (zh) | 星载激光测高仪指向角误差为非常数时的在轨标定方法 | |
CN114114359B (zh) | 单星与地基设备联合的再入预报方法、装置和电子设备 | |
Xu et al. | Stellar radiation modeling and image simulation for airborne daytime star sensor | |
CN112540367A (zh) | 空间目标雷达定轨实时识别方法、设备和存储介质 | |
CN115543637B (zh) | 空间目标的关联方法、装置和存储介质 | |
Štěpánek et al. | Impact of orbit modeling on DORIS station position and Earth rotation estimates | |
CN115578650B (zh) | 一种基于地表温度图像预测干旱指数的方法及装置 | |
Yuan et al. | Autonomous estimation of angle random walk of fiber optic gyro in attitude determination system of satellite | |
Yinhu et al. | Rigorous and integrated self-calibration model for a large-field-of-view camera using a star image | |
CN111156990A (zh) | 基于指向自动测定的空间碎片实时天文定位和测光方法 | |
Kuznetsov et al. | Multifunctional astronomical self-organizing system of autonomous navigation and orientation for artificial Earth satellites | |
CN114386282B (zh) | 半分析法的低轨巨型星座轨道动力学分析方法及装置 | |
Paulsen et al. | Calibration and verification of the TERMA star tracker for the NEMO satellite | |
CN114993341A (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 |