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

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

Estimation of ballistic coefficients of low altitude debris objects from historical two line elements

This paper presents a new method for estimating ballistic coefficients (BCs) of low perigee debris objects from their historical two line elements (TLEs). The method uses the drag perturbation equation of the semi-major axis of the orbit. For an object with perigee altitude below 700 km, the variation in the mean semi-major axis derived from the TLE is mainly caused by the atmospheric drag effect, and therefore is used as the source in the estimation of the ballistic coefficient. The method is tested using the GRACE satellites, and a number of debris objects with external ballistic coefficient values, and agreements of about 10% are achieved.     

Instability of the present LEO satellite populations

Several studies conducted during 1991–2001 demonstrated, with some assumed launch rates, the future unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects. In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new breakup debris due to collisions would exceed the loss of objects due to orbital decay.     

Towards accurate atmospheric mass density determination Using precise positional information of space objects

This paper presents a new method of deriving atmospheric mass densities with a high temporal resolution from precise orbit data of low earth orbiting (LEO) space objects. This method is based on the drag perturbation equation of the semi-major axis of the orbit of LEO space objects which relates the change rate of the semi-major axis to the atmospheric mass density. The effectiveness of the new method is evaluated using the GFZ-ISDC GPS rapid science orbit (RSO) products of the CHAMP satellite over a time period of 3 months. The densities derived using this new method and obtained from accelerometer data are compared and good agreements are achieved. An example of using the derived density to generate good orbit prediction for CHAMP is presented.     

Characterization of the cataloged Fengyun-1C fragments and their long-term effect on the LEO environment

The intentional breakup of Fengyun-1C on 11 January 2007 created the most severe orbital debris cloud in history. The altitude where the event occurred was probably the worst location for a major breakup in the low Earth orbit (LEO) region, since it was already highly populated with operational satellites and debris generated from previous breakups. The addition of so many fragments not only poses a realistic threat to operational satellites in the region, but also increases the instability (i.e., collision cascade effect) of the debris population there.     

远距离轨道接近及绕飞控制技术研究

文章对半长轴、偏心率、倾角、升交点赤经相近的两航天器的远距离轨控接近方案进行了理论分析。首先研究了通过改变轨道半长轴调整航天器升交点赤经变化率,进而控制轨道法向振幅的方法;其次给出控制点的选择时机与相对绕飞椭圆大小的定量关系;然后提出采用径向矢径模之差估计远距离相对运动参数方法;最后通过"神舟七号"伴星与轨道舱的远距离接近及绕飞控制仿真对理论推导进行了验证。     

Accurate orbit predictions for debris orbit manoeuvre using ground-based lasers

Orbit manoeuvre of low Earth orbiting (LEO) debris using ground-based lasers has been proposed as a cost-effective means to avoid debris collisions. This requires the orbit of the debris object to be determined and predicted accurately so that the laser beam can be locked on the debris without the loss of valuable laser operation time. This paper presents the method and results of a short-term accurate LEO (<900 km in altitude) debris orbit prediction study using sparse laser ranging data collected by the EOS Space Debris Tracking System (SDTS). A main development is the estimation of the ballistic coefficients of the LEO objects from their archived long-term two line elements (TLE). When an object is laser tracked for two passes over about 24 h, orbit prediction (OP) accuracy of 10–20 arc seconds for the next 24–48 h can be achieved – the accuracy required for laser debris manoeuvre. The improvements in debris OP accuracy are significant in other applications such as debris conjunction analyses and the realisation of daytime debris laser tracking.     

LEO空间碎片甚短弧角度数据初轨确定方法对比

光学观测是空间目标观测中最常见的一种观测方式。采用扫描模式工作时光学观测得到的观测弧段弧长通常很短,有时甚至不到被观测空间目标运行周期的1%,这样的角度数据被称为甚短弧角度数据。基于近圆LEO空间碎片地基实测场景,研究比较仅利用角度数据进行初始轨道确定常用方法的性能差异,分析观测弧长对不同初轨确定算法的定轨成功率和误差的影响,为初轨确定工作提供参考。对比分析了常用的几种方法,包括Laplace方法、Gauss方法、Gooding方法和近几年提出的距离搜索算法等。大规模实测数据处理结果显示,距离搜索算法的成功率高于90%,初轨半长轴统计误差仅为25 km。初轨结果表明,距离搜索算法定轨成功率高于其他算法。研究成果可为解决空间碎片初轨确定问题提供参考。      

Analysis of the orbit lifetime of CubeSats in low Earth orbits including periodic variation in drag due to attitude motion

It is estimated that more than 22,300 human-made objects are in orbit around the Earth, with a total mass above 8,400,000 kg. Around 89% of these objects are non-operational and without control, which makes them to be considered orbital debris. These numbers consider only objects with dimensions larger than 10 cm. Besides those numbers, there are also about 2000 operational satellites in orbit nowadays. The space debris represents a hazard to operational satellites and to the space operations. A major concern is that this number is growing, due to new launches and particles generated by collisions. Another important point is that the development of CubeSats has increased exponentially in the last years, increasing the number of objects in space, mainly in the Low Earth Orbits (LEO). Due to the short operational time, CubeSats boost the debris population. One of the requirements for space debris mitigation in LEO is the limitation of the orbital lifetime of the satellites, which needs to be lower than 25 years. However, there are space debris with longer estimated decay time. In LEÓs, the influence of the atmospheric drag is the main orbital perturbation, and is used in maneuvers to increment the losses in the satellite orbital energy, to locate satellites in constellations and to accelerate the decay.The goal of the present research is to study the influence of aerodynamic rotational maneuver in the CubeSat?s orbital lifetime. The rotational axis is orthogonal to the orbital plane of the CubeSat, which generates variations in the ballistic coefficient along the trajectory. The maneuver is proposed to accelerate the decay and to mitigate orbital debris generated by non-operational CubeSats. The panel method is selected to determine the drag coefficient as a function of the flow incident angle and the spinning rate. The pressure distribution is integrated from the satellite faces at hypersonic rarefied flow to calculate the drag coefficient. The mathematical model considers the gravitational potential of the Earth and the deceleration due to drag. To analyze the effects of the rotation during the decay, multiple trajectories were propagated, comparing the results obtained assuming a constant drag coefficient with trajectories where the drag coefficient changes periodically. The initial perigees selected were lower than 400 km of altitude with eccentricities ranging from 0.00 to 0.02. Six values for the angular velocity were applied in the maneuver. The technique of rotating the spacecraft is an interesting solution to increase the orbit decay of a CubeSat without implementing additional de-orbit devices. Significant changes in the decay time are presented due to the increase of the mean drag coefficient calculated by the panel method, when the maneuver is applied, reducing the orbital lifetime, however the results are independent of the angular velocity of the satellite.     

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

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

Apparent rotation properties of space debris extracted from photometric measurements

Knowledge about the rotation properties of space debris objects is essential for the active debris removal missions, accurate re-entry predictions and to investigate the long-term effects of the space environment on the attitude motion change. Different orbital regions and object’s physical properties lead to different attitude states and their change over time.Since 2007 the Astronomical Institute of the University of Bern (AIUB) performs photometric measurements of space debris objects. To June 2016 almost 2000 light curves of more than 400 individual objects have been acquired and processed. These objects are situated in all orbital regions, from low Earth orbit (LEO), via global navigation systems orbits and high eccentricity orbit (HEO), to geosynchronous Earth orbit (GEO). All types of objects were observed including the non-functional spacecraft, rocket bodies, fragmentation debris and uncorrelated objects discovered during dedicated surveys. For data acquisition, we used the 1-meter Zimmerwald Laser and Astrometry Telescope (ZIMLAT) at the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald, Switzerland. We applied our own method of phase-diagram reconstruction to extract the apparent rotation period from the light curve. Presented is the AIUB’s light curve database and the obtained rotation properties of space debris as a function of object type and orbit.     

Feasibility of performing space surveillance tasks with a proposed space-based optical architecture

Under ESA contract an industrial consortium including Aboa Space Research Oy (ASRO), the Astronomical Institute of the University of Bern (AIUB), and the Dutch National Aerospace Laboratory (NLR), proposed the observation concept, developed a suitable sensor architecture, and assessed the performance of a space-based optical (SBO) telescope in 2005. The goal of the SBO study was to analyse how the existing knowledge gap in the space debris population in the millimetre and centimetre regime may be closed by means of a passive optical instrument. The SBO instrument was requested to provide statistical information on the space debris population in terms of number of objects and size distribution. The SBO instrument was considered to be a cost-efficient with 20 cm aperture and 6° field-of-view and having flexible integration requirements. It should be possible to integrate the SBO instrument easily as a secondary payload on satellites launched into low-Earth orbits (LEO), or into geostationary orbit (GEO). Thus the selected mission concept only allowed for fix-mounted telescopes, and the pointing direction could be requested freely. Since 2007 ESA focuses space surveillance and tracking activities in the Space Situational Awareness (SSA) preparatory program. Ground-based radars and optical telescopes are studied for the build-up and maintenance of a catalogue of objects. In this paper we analyse how the proposed SBO architecture could contribute to the space surveillance tasks survey and tracking. We assume that the SBO instrumentation is placed into a circular sun-synchronous orbit at 800 km altitude. We discuss the observation conditions of objects at higher altitude, and select an orbit close to the terminator plane. A pointing of the sensor orthogonal to the orbital plane with optimal elevation slightly in positive direction (0° and +5°) is found optimal for accessing the entire GEO regime within one day, implying a very good coverage of controlled objects in GEO, too. Simulations using ESA’s Program for Radar and Optical Observation Forecasting (PROOF) in the version 2005 and a GEO reference population extracted from DISCOS revealed that the proposed pointing scenario provides low phase angles together with low angular velocities of the objects crossing the field-of-view. Radiometric simulations show that the optimal exposure time is 1–2 s, and that spherical objects in GEO with a diameter of below 1 m can be detected. The GEO population can be covered under proper illumination nearly completely, but seasonal drops of the coverage are possible. Subsequent observations of objects are on average at least every 1.5 days, not exceeding 3 days at maximum. A single observation arc spans 3° to 5° on average. Using a simulation environment that connects PROOF to AIUB’s program system CelMech we verify the consistency of the initial orbit determination for five selected test objects on subsequent days as a function of realistic astrometric noise levels. The initial orbit determination is possible. We define requirements for a correlator process essential for catalogue build-up and maintenance. Each single observation should provide an astrometric accuracy of at least 1”–1.5” so that the initially determined orbits are consistent within a few hundred kilometres for the semi-major axis, 0.01 for the eccentricity, and 0.1° for the inclination.     

低轨导航增强卫星星座设计

针对目前全球低轨卫星快速发展的现状,对低轨导航增强卫星星座设计方法进行了详细的研究。首先推导了轨道高度与可视球冠的关系,结合太空垃圾分布,从覆盖范围、经济性及碰撞风险几方面联合确定了轨道高度。然后推导了用户仰角与轨道倾角的关系,分析了实现南北极点覆盖的轨道倾角。接着结合铱星星座,推导出单一星座构型无法实现全球范围内均匀的可见星和精度衰减因子（Dilution of Precision,DOP）值分布。最后提出了一种组合低轨卫星星座设计方法。结果表明,该方法设计的组合星座在实现全球覆盖的同时,能够实现可见星数量与DOP值在全球范围内的均匀分布。     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

基于天基光学探测图像初析南大西洋异常区影响

文章针对目前在轨运行的高灵敏度探测类可见光遥感器经过南大西洋异常区（SouthAtlanticAnomaly,SAA）时出现的探测能力降低、目标无法识别等现象,通过卫星在轨获取的SAA天基探测图像研究图像中出现的异常现象机理;提出了基于CCD探测器固有暗像元区的图像特征变化开展SAA影响评估的方法,根据在轨图像暗像元区的均值和标准差、图像区的热像素及活像元等特征参数的统计分析,判断经过SAA过程中连续获取的图像中出现显著影响时的拐点,通过曲线拟合估算出该轨道高度的SAA空间分布特征和边界,依此进一步提出利用在轨天基探测图像给出测量影响该轨道高度卫星SAA区域边界的简单方法;最后,基于SAA对天基探测器件的影响分析结果和认识,给出了高灵敏度探测类卫星避免或消除SAA影响的改进措施和建议。文章采用的数据分析和评估方法为提升SAA对天基探测影响的认识,以及下一步有效开展空间碎片探测工作提供参考。     

Yarkovsky-Schach effect on space debris motion

The Yarkovsky-Schach effect is a small perturbation affecting Earth satellites and space debris illuminated by the Sun. It was first applied to the orbit of LAGEOS satellites as an explanation of the residuals in orbital elements. In this work, we carry out several numerical integration tests taking into consideration various orbit and rotation parameters, in order to analyse this effect in a broader context. The semi-major axis variations remain small and depend on the spin axis attitude with respect to the Sun. We show that the force amplitude is maximised for orbits inclined with $i\approx$?20–30°. We also observe the influence on other orbital elements, notably on the orbit inclination. However, these effects are clearly observed only on long timescales; in our simulations, we propagated the orbits for 200?y. The Yarkovsky-Schach effect is thus confirmed to have a minuscule magnitude. It should be taken into account in studies requiring high-precision orbit determination, or on expanded timescales.     

Orbital debris–debris collision avoidance

We focus on preventing collisions between debris and debris, for which there is no current, effective mitigation strategy. We investigate the feasibility of using a medium-powered (5 kW) ground-based laser combined with a ground-based telescope to prevent collisions between debris objects in low-Earth orbit (LEO). The scheme utilizes photon pressure alone as a means to perturb the orbit of a debris object. Applied over multiple engagements, this alters the debris orbit sufficiently to reduce the risk of an upcoming conjunction. We employ standard assumptions for atmospheric conditions and the resulting beam propagation. Using case studies designed to represent the properties (e.g. area and mass) of the current debris population, we show that one could significantly reduce the risk of nearly half of all catastrophic collisions involving debris using only one such laser/telescope facility. We speculate on whether this could mitigate the debris fragmentation rate such that it falls below the natural debris re-entry rate due to atmospheric drag, and thus whether continuous long-term operation could entirely mitigate the Kessler syndrome in LEO, without need for relatively expensive active debris removal.     

Analysis of tape tether survival in LEO against orbital debris

The low earth orbit (LEO) environment contains a large number of artificial debris, of which a significant portion is due to dead satellites and fragments of satellites resulted from explosions and in-orbit collisions. Deorbiting defunct satellites at the end of their life can be achieved by a successful operation of an Electrodynamic Tether (EDT) system. The effectiveness of an EDT greatly depends on the survivability of the tether, which can become debris itself if cut by debris particles; a tether can be completely cut by debris having some minimal diameter. The objective of this paper is to develop an accurate model using power laws for debris-size ranges, in both ORDEM2000 and MASTER2009 debris flux models, to calculate tape tether survivability. The analytical model, which depends on tape dimensions (width, thickness) and orbital parameters (inclinations, altitudes) is then verified with fully numerical results to compare for different orbit inclinations, altitudes and tape width for both ORDEM2000 and MASTER2009 flux data.     

全方向长期伴飞卫星集群轨道设计与仿真

主要研究在J2摄动的影响下,在中心卫星周围全方向上部署长期相对有界伴飞卫星的轨道设计方法。首先,利用J2不变的条件,由已知中心卫星轨道根数建立伴飞卫星的半长轴、偏心率及轨道倾角之间关系。然后,通过对伴飞卫星和中心卫星之间初始相对位置及轨道根数之间关系的分析,得到伴飞卫星的全部初始轨道根数,使其能够部署到中心卫星周围的任意方向上,从而能实现卫星集群在中心卫星的周围全方向长期有界伴飞。最后,通过轨道仿真验证了该方法的有效性。     

大气阻力参数修正对低轨空间目标轨道预报精度的改进

对于低轨空间目标, 大气阻力是影响轨道预报精度的主要摄动力. 本文提出了一种 基于空间环境数据和神经网络模型的空间目标大气阻力参数修正方法, 基于目 标的历史两行元根数, 通过模拟得到外推一天轨道预报中预报结果与观测数据 符合最好的阻力调制系数, 分析表明其与太阳F_10.7指数和地磁Ap指数具有很好的相关性. 根据已有数据, 构建神经网络模型, 实现对阻力调制系数 的补偿计算, 从而改进低轨目标外推一天的轨道预报. 结果表明, 神经网络模 型相比两行元根数能够更及时地对空间环境变化进行响应. 将该方案应用于天 宫一号和国际空间站的外推一天轨道预报, 验证了方案的正确性和普适性, 对 地磁扰动引起的较大预报误差改进效果更好, 误差能够降低50%~60%; 平均而言, 预报精度可以提高约30%, 改进成功率达到80%左右.