空间碎片环境与碰撞预警仿真系统设计

在轨编目航天器数目已超过了9000,为了保证航天任务的顺利完成,必须对可能威胁到任务轨道的空间目标进行分析,文章首先对两行轨道根数和SGP4／SDP4轨道预报模型进行了分析,在此基础上建立了空间目标数据库及空间碎片环境与碰撞预警仿真和数据可视化系统。     

TLE预报精度改进及碰撞预警中的应用

双行根数（Two Line Elements,TLE）是目前唯一公开发布且编目最完备的地球轨道空间目标编目数据,由于其根数的特殊性,必须配合SGP4/SDP4 使用导致其预报精度有限. 通过对TLE+SGP4/SDP4预报误差源的分析,利用历史TLE数据生成准观测数据重新进行轨道拟合,针对不同轨道高度和面积质量比情况,确定相应拟合周期,得出了基于弹道系数先验信息有效确定面积质量比的经验方法,对多个独立时段拟合后轨道根数结合数值方法的预报误差进行统计,结果表明,拟合后的TLE根数预报精度与稳定性有很大程度提升. 利用实例对拟合TLE在碰撞预警中的应用进行分析,认为该方法对航天工程中碰撞预警置信度的提高具有重要意义.      

空间目标碰撞预警的协方差计算与应用

碰撞概率是碰撞预警工程中空间目标危险交会的重要判据之一,其计算精度会受到预报协方差计算精度的影响. 本文统计计算了两种形式的预报协方差. 一是利用精密数值预报模型对卫星精密根数进行预报,将预报根数与精密根数的差作为样本,统计得出1～7天不同预报期的预报协方差;二是采用 SGP4/SDP4预报模型对TLE数据进行预报,将预报根数与TLE根数的差作为样本,统计得出1～7d不同预报期的预报协方差. 分别分析两种方法中包含变轨过程和无变轨情况下的轨道预报精度. 结合2012年某低轨道卫星的危险交会,分析了采用不同协方差时,协方差精度对碰撞预警精度的影响. 协方差计算可为实际碰撞预警工程提供参考.      

基于驻留时间比的近地空间碎片环境建模

为分析近地空间碎片的分布规律,提出了一种以碎片在空间网格内驻留时间为基础的碎片环境统计建模方法.该方法利用多项式拟合和求根方法统计碎片在空间网格内的停留时间,获取模型基础数据,并据此采用多项式预测、插值和时间序列分析等技术,综合分析空间碎片的分布与演化规律.给出了一个基于双行根数(TLE,Two Line Elements)数据的建模实例,该实例通过了ORDEM2000模型的对比验证,并获得了一些更精细的近地空间碎片环境特征.所得建模方法和分析结论可为长期运行的近地航天器轨道设计、碰撞风险评估及防护等提供技术支撑.     

同步轨道通信卫星占频保轨现象分析与实测

为了应对国际电信联盟对卫星频轨资源的相关规定，卫星运营商可能会选择将卫星送至另一个不同的轨位并连续保持 90 天来启用该位置的频率指配.这种变轨行为称作“占频保轨”，是一些同步通信卫星的常规操作.本文以国际通信组织的 INTELSAT 系列卫星为研究对象，基于双行根数 (TLE) 和SGP4/SDP4模型计算了67个INTELSAT (国际通信卫星) 的自发射以来的轨道数据，提出了基于密度和k 最近邻的离群点检测来筛选数据野值的方法，提高了数据准确性；对这些卫星的占频保轨行为进行统计和分析，提出我国卫星运营商应充分利用国际规则，维护自身太空权益的建议.     

基于预报偏差的LEO航天器轨道异常检测

针对由轨道控制、大气环境、碰撞等因素造成的低轨(LEO)航天器轨道突变问题,提出了一种基于预报偏差的轨道异常检测方法。选择LEO轨道的半长轴和倾角作为特征轨道参数,利用SGP4模型长期项对目标的两行轨道要素(TLE)进行预报得到特征轨道参数的预报值,通过对特征轨道参数的编目数据和预报数据进行平滑后求差得到预报偏差序列,基于马氏距离对预报偏差数据的两个分量进行联合异常检测。对Terra卫星2010年的机动事件分析结果同NASA发布的其机动历史相吻合,表明该方法可以有效地检测航天器轨道异常的次数、时间和类型,可应用于空间目标的监视与空间态势的感知。     

基于期望最大算法的空间事件及异常值探测

美国ELSET数据库提供的TLE数据是目前使用最广泛的数据,在热层大气密度反演、弹道系数估计、碰撞预警等领域具有重要作用。受空间环境扰动、空间事件以及TLE产生过程等共同影响,ELSET数据库包含大量亟待清理的异常值和识别的空间事件,例如发布错误的TLE、轨道根数异常和Bstar异常。现有方法在清理异常轨道根数时缺乏统一性,需要使用不同的技术,清理流程较为繁杂,并且仅适用于特定轨道区域的少数目标。为克服现有方法的弊端,提出了一种基于期望最大算法的滑动窗口–多项式拟合预报方法,对含有轨道机动的碎片以及受空间环境影响的碎片进行异常值与空间事件探测。研究表明,该方法能够灵活处理不同空间环境下的异常值与空间事件探测,具有普适性,适用于所有轨道碎片。      

基于多项式近似的空间碎片群体轨道预报算法

快速准确地分析空间碎片群轨道演化行为对于其他在轨航天器碰撞规避至关重要。在各摄动力的作用下,空间碎片群演化运动呈现出复杂的非线性特征。空间碎片群体个体数量巨大,如果通过对空间碎片群中每个空间碎片进行轨道积分来分析群体预报的方法会导致计算量过大。针对该问题,提出一种基于多项式近似的轨道快速预报分析方法。该方法将空间碎片群分为少量的标称碎片和其他大量关联碎片。针对标称碎片的轨道预报采用数值积分求解保证预报精度;而针对其他大量的关联碎片轨道预报问题,采用多项式泰勒展开半解析方法求解,从而在保证预报精度的前提下有效减少空间碎片群轨道预报的计算量。为了验证方法的有效性,对不同空间碎片群进行了轨道预报仿真。仿真结果表明,当轨道预报精度设定在1m范围内时,多项式近似算法的计算量较蒙特卡洛方法计算效率提高了2.2~17.2倍,验证了所提出方法的有效性。     

利用TLE数据分析LEO卫星轨道异常的新方法——————综合判据法

及时准确地发现在轨卫星的轨道异常意义重大. 通过有效的异常算法, 能够找出发生轨道异常的碎片或航天器, 为空间碎片碰撞预警系统分析和验证碰撞事件提供数据支持. 通过对利用TLE (Two Line Elements)数据分析LEO在轨卫星轨道异常的方法研究, 提出了一个利用单个卫星相邻根数时间差控制加综合判据的判别方法. 分析表明, 相对于取单一因素阈值的判别方法, 综合判据法能够最大限度地减少漏判, 并且保持相对较高的判断准确率.      

一种基于TLE数据的轨道异常分析方法

空间在轨物体的轨道异常是航天工程及预警领域普遍关注的问题,及时发现轨道异常意义重大,通过分析空间物体的轨道异常,可以及时发现和识别规避事件或碰撞事件,还可以了解监测网的能力.本文提出一种基于TLE数据的简单的轨道异常分析方法——长半轴变化法.该方法快速有效,应用到低轨在用卫星和美俄解体碎片的异常分析中,异常物体正确识别率可达到100%;对美俄解体碎片进行轨道异常分析后得出,美国空间监视网可以稳定探测90%以上的解体碎片.      

地球卫星运动中坐标系附加摄动与参考系选择问题

针对地球赤道面的摆动,曾提出过两种空间地心赤道坐标系,即历元地心平赤道坐标系和混合形式的轨道坐标系.在后一种坐标系中虽然可避免计算因地球引力位变化所引起的坐标系附加摄动的麻烦,但这又带来了两种坐标系之间的不统一和转换的麻烦.针对实际现状,直接在历元地心平赤道坐标系中给出相应的坐标系附加摄动解,所有与此有关的转换问题,均可在历元地心平赤道坐标系中进行,而无需为了避免坐标系附加摄动的计算去通过轨道坐标系来实现.从而建议,不必再引进轨道坐标系,在卫星定轨和预报中,所有的转换问题均可采用同一坐标系,即历元(目前是J2000.0)地心平赤道坐标系.      

Propagation errors analysis of TLE data

Orbit position uncertainty is an important factor for collision avoidance issues. For a single object with high frequency historical data, we can attain its position uncertainty easily. But sometimes data is not enough for errors analysis, orbits need to be classified. In this paper error analysis is made from two-line element sets data (TLEs). The Simplified General Perturbations-4 (SGP4) propagator was used. Statistical errors of debris and R/B are given for lower-altitude orbits which are classified by perigee altitude and eccentricity. The errors results and analysis for SSO (the Sun synchronous orbit) typical orbits are obtained. At last atmospheric drag as a main cause of downrange errors in lower-altitude orbit is analyzed. BSTAR in TLEs is modified to improve prediction precision.     

满州里上空电离层之平均特性及暴时特性

本文对1983年满州里(磁纬38.07°)的测高仪资料和子午仪卫星资料进行了处理.前者主要是对F₂层临界频率f₀F₂,而后者是对微分多普勒的总含量N_T进行分析.得到满州里上空电离层之平均特性,以及磁暴期间电离层的一些异常变化.结果发现,满州里电离层的平均特性具有典型的季节异常,且f₀F₂和N_T的变化呈现较好的一致性.满州里上空电离层的扰动大都发生在两分点季和冬季,且正相大多发生在冬季,双相大多发生在两分点季,负相大多发生在夏季和两分点季.得出1983年满州里暴时N_T和f₀F₂相对于月中值偏移量的均方根值δN_T和δf₀F₂都与磁Ap指数均方根值δAp之间存在一元线性相关.      

Variation of hmF2 and NmF2 deduced from DPS-4 over Multan (Pakistan) and their comparisons with IRI-2012 & IRI-2016 during the deep solar minimum between cycles 23 & 24

We report the results of ionospheric measurements from DPS-4 installed at Multan (Geog coord. 30.18°N, 71.48°E, dip 47.4°). The variations in F2-layer maximum electron density NmF2 and its peak height hmF2 are studied during the deep solar minimum between cycles 23 & 24 i.e 2008–2009 with comparisons conducted with the International Reference Ionosphere (IRI) versions 2012 & 2016. We find that the hmF2 observations peak around the pre-sunrise and sunrise hours depending on the month. Seasonally, the daytime variation of NmF2 is higher in the Equinox and Summer, while daytime hmF2 are slightly higher in the Equinox and Winter. High values of hmF2 around midnight are caused by an increase of upward drifts produced by meridional winds. The ionosphere over Multan, which lies at the verge of low and mid latitude, is affected by both $E\times B$ drifts and thermospheric winds as evident from mid-night peaks and near-sunrise dips in hmF2. The results of the comparison of the observed NmF2 and hmF2 for the year 2008–2009 with the IRI-2012 (both NmF2 and hmF2) and IRI-2016 (only hmF2) estimates indicate that for NmF2, IRI-2012 with Consultative Committee International Radio (CCIR) option produces values in better agreement with observed data. Whereas, for hmF2, IRI-2016 with both International Union of Radio Science (URSI) and CCIR SHU-2015 options, predicts well for nighttime hours throughout the year. However, the IRI-2012 with CCIR option produces better agreement with data during daytime hours. Furthermore, IRI-2012 with CCIR option gives better results during Equinox months, whereas, IRI-2016 with both URSI and CCIR SHU-2015 options predict well for Winter and Summer.     

Comparison of diurnal,seasonal and latitudinal effect of MUF VR and NmF2 VR during some solar cycle epochs

The diurnal and seasonal changes of the variability (VR) of Maximum Useable Frequency (MUF) are compared with those of peak electron density (NmF2) at Ibadan (7.4°N, 3.9°E, 6°S dip) in the African sector. Also compared is the latitudinal effect on both characteristics by combining data from Singapore (1.3°N, 103.8°E, 17.6°S dip) in the East Asian sector and Slough (51.5°N, 359.4°E, 66.5°S dip) in the European sector. MUF VR is found to be about half of NmF2 VR at all the hours and seasons and during the solar cycle epochs considered for the three stations. While nighttime MUF VR is greater in June Solstice and September Equinox during both low and moderate solar activities and in September Equinox and December Solstice during high solar activity, nighttime NmF2 VR is greater in June Solstice and September Equinox during high solar activity and greater at the equinoxes during low and moderate solar activities. This signifies a shift in nighttime MUF peak VR from the middle six months during low and moderate solar activities to the last half of the year during high solar activity. Daytime VR of both characteristics are not observed to show any seasonal variation. MUF VR and that of NmF2 are found to increase and decrease alternately with the Zurich sunspot number (Rz) for Ibadan and Singapore. For Slough, the VR of both characteristics increases with Rz during the first half of the day. It then increases and decreases alternately with Rz during the remaining hours of the day. While nighttime MUF VR decreases with latitude, just like nighttime NmF2 VR, no latitudinal effect is found for daytime VR of both characteristics.     

Improved orbit predictions using two-line elements

The density of orbital space debris constitutes an increasing environmental challenge. There are two ways to alleviate the problem: debris mitigation and debris removal. This paper addresses collision avoidance, a key aspect of debris mitigation. We describe a method that contributes to achieving a requisite increase in orbit prediction accuracy for objects in the publicly available two-line element (TLE) catalog. Batch least-squares differential correction is applied to the TLEs. Using a high-precision numerical propagator, we fit an orbit to state vectors derived from successive TLEs. We then propagate the fitted orbit further forward in time. These predictions are validated against precision ephemeris data derived from the international laser ranging service (ILRS) for several satellites, including objects in the congested sun-synchronous orbital region. The method leads to a predicted range error that increases at a typical rate of 100 m per day, approximately a 10-fold improvement over individual TLE’s propagated with their associated analytic propagator (SGP4). Corresponding improvements for debris trajectories could potentially provide conjunction analysis sufficiently accurate for an operationally viable collision avoidance system based on TLEs only.