基于相控阵雷达波束篱笆的空间碎片数量与分布估计方法

随着载人航天与空间站等航天活动的增多,不能有效防护、也无法定期跟踪和编目的小尺寸(尤其是1～10 cm)碎片的危害越来越受到关注,这些碎片信息的获取依赖于统计采样技术.针对简化的相控阵雷达波束篱笆空间碎片探测模式,提出了一种采用统计技术估计空间碎片总数量以及高度和倾角分布的方法.将碎片穿越波束篱笆的过程用Poisson分布来建模,根据观测时段内穿越波束篱笆目标的平均到达率及测量的轨道高度和倾角数据来估计给定轨道高度范围或倾角范围内碎片的数量,进而得到碎片的总数量以及碎片数量随轨道高度或倾角的分布.在获取雷达散射截面信息时,该方法还可用于估计碎片数量随尺寸的分布.通过仿真实验验证了该方法的有效性.      

Characterizing space debris longitude-dependent distribution based on RAAN perturbation rate

The space debris environment is one of the major threats against payloads. Space debris orbital distribution is of great importance for space debris environment modeling. Due to perturbation factors, the Right Ascension of Ascending Node (RAAN) of space objects changes consistently, causing regular rotation of the orbit plane around Earth’s axis. Based on the investigation of the RAAN perturbation rate of concerned objects, this paper proposes a RAAN discretization method in order to present the space debris longitude-dependent distribution. Combined with two line element (TLE) data provided by the US Space Surveillance Network, the estimated value from RAAN discretization method is compared with the real case. The results suggest that using only the initial orbital data at the beginning of the time interval of interest, the RAAN discretization method is able to provide reliable longitude distribution of concerned targets in the next following period. Furthermore, spacecraft cumulative flux against space debris is calculated in this paper. The results suggest that the relevance between spacecraft RAAN setup and flux output is much smaller for LEO targets than MEO targets, which corresponds with the theory analysis. Since the nonspherical perturbation is the major factor for RAAN variation, the RAAN perturbation rate has little connection with the size of orbital objects. In other words, the RAAN discretization method introduced in this paper also applies to space debris of different size range, proposing a possible suggestion for the improvement of space debris environment engineering models.     

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.     

空间目标预警阈值确定分析

针对目前通常采用计算碰撞概率来进行空间目标碰撞预警的方法,为提高其评判标准的精度以达到提高航天器精准预警的目的,分析了目前空间目标碰撞概率阈值确定过程的总体框架.对阈值确定过程中所涉及的误差协方差矩阵确定方法、碰撞概率计算以及碎片通量计算方法进行了研究.以国际空间站（ISS）为例,通过计算其所在空间的碎片通量,分析其轨道受空间环境的影响情况,并与早期结果进行比较,得出其以往的碰撞概率阈值可能被高估的结论.      

空间碎片自然衰减和太阳活动

分析研究了空间碎片数随太阳辐射流量Ｆ１０．７的变化;给出预报Ｆ１０．７长期变化的计算方法和预测空间碎片数的数学模型。结果显示：①强太阳活动造成空间碎片年增长率下降;②空间碎片数与太阳活动１１年变化密切相关,相关数为０．９;③空间碎片增长率约为发射率的两倍;④若发射率保持不变,则到２０２０年,大于１０ｃｍ的碎片数将达到１４５００;⑤若小碎片的增长为大碎片增长的两倍,则到２０２０年,大于１ｃｍ的碎片数可达１２５０００。     

柱状弹丸撞击防护屏形成碎片云材料状态特性研究

微流星体及空间碎片的高速撞击威胁着长寿命、大尺寸航天器的安全运行，导致其严重的损伤和灾难性的失效。为精确估计微流星体及空间碎片高速撞击防护屏所产生碎片云对舱壁的损伤，必须确定碎片云中三种状态材料的特性，建立了碎片云特性分析模型，分别计算了柱状弹丸撞击防护屏所产生碎片云以及碎片云中弹丸和防护屏材料三种状态物质的质量分布。通过计算分析可见，弹丸以不同速度撞击防护屏所产生碎片云三种状态物质的质量分布是不同的，速度增大，液化和气化增强，对靶件的损伤小。而在速度小于7km/s时，碎片云以固体碎片的形式存在，对靶件的损伤大。     

雷达波束驻留模式下空间碎片数量的置信区间估计

碎片数量估计是空间碎片环境统计特征描述的重要内容之一,对于空间碎片环境模型验证、航天器碰撞风险分析以及碎片数量增长趋势预测有重要意义.针对波束指向正东、正南任意仰角的雷达波束驻留(Beam-park)模式(天顶指向是波束指向仰角为90°时的特例),给出了一种估计碎片数量置信区间的方法.对于给定轨道高度范围内一个具有穿越雷达波束可能性(即雷达散射截面足够大,且轨道倾角相对测站纬度足够大)的碎片,将其是否真正穿越波束这一随机事件用(0-1)分布来建模,根据所采集的轨道高度和倾角数据,计算出该轨道高度范围内碎片穿越波束的平均概率,进而采用中心极限定理来估计碎片数量的置信区间.仿真结果表明了方法的有效性.      

Analysis of laser ranges and angular measurements data fusion for space debris orbit determination

In the framework of space debris, the orbit determination process is a fundamental step, both, for researchers and for satellite operators. The accurate knowledge of the orbit of space debris objects is needed to allow space debris characterization studies and to avoid unnecessary collision avoidance maneuvers.The accuracy of the results of an orbit determination process depends on several factors as the number, the accuracy, the kind of processed measurements, their distribution along the orbit, and the object-observer relative geometry. When the observation coverage of the target orbit is not homogeneous, the accuracy of the orbit determination can be improved processing different kind of observables. Recent studies showed that the satellite laser ranging technique can be successfully applied to space debris.In this paper, we will investigate the benefits of using laser ranges and angular measurements for the orbit determination process. We will analyze the influence of the number of used observations, of the covered arc of orbit, of each observable, and of the observation geometry on the estimated parameters. Finally, using data acquired on short observation arcs, we analyze the achievable accuracies for the orbital regimes with the highest space debris density, and to the consequences of the data fusion on catalog maintenance operations. The results shown are obtained using only real data (both angular and laser measurements) provided by sensors of the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald owned by the Astronomical Institute of the University of Bern (AIUB) and for some studies also using ranges provided from other stations of the International Laser Ranging Service (ILRS).     

Impact effects from small size meteoroids and space debris

When the impact risk from meteoroids and orbital debris is assessed the main concern is usually structural damage. With their high impact velocities of typically 10–20 km/s millimeter or centimeter sized objects can puncture pressure vessels and other walls or lead to destruction of complete subsystems or even whole spacecraft. Fortunately chances of collisions with such larger objects are small (at least at present). However, particles in the size range 1–100 μm are far more abundant than larger objects and every orbiting spacecraft will encounter them with certainty. Every solar cell (8 cm² area) of the Hubble Space Telescope encountered on average 12 impacts during its 8.25 years of space exposure. Most were from micron sized particles.     

基于模糊应力—强度模型的空间系绳碰撞可靠性

空间系绳系统由于其特殊的组成结构日益受到关注,空间系绳的碰撞可靠性研究是系绳任务设计的重要一环.本文采用可靠性分析基本原理中的应力—强度模型,根据近地轨道的空间碎片通量和泊松分布方法,进行空间系绳在轨碰撞可靠性的研究分析.根据对空间系绳碰撞可靠性有影响的系绳结构因素,即单股系绳的直径、长度,双股系绳的绳间距、碰撞角等,对空间碎片撞击切割系绳后碰撞点处的系绳剩余截面进行建模,利用模糊应力—强度模型计算空间系绳在撞击后发生毁灭性碰撞的概率,进而根据泊松方法计算空间系绳在轨可靠性随时间的变化,通过仿真分析,对比不同结构空间系绳的有效在轨时间.      

空间系绳碎片碰撞生存能力研究

空间系绳结构和几何尺寸对其在复杂空间环境中抵抗碎片碰撞的能力有显著影响.本文根据系绳-碎片碰撞模型和泊松分布,计算了双股系绳和带状系绳被碎片割断的概率,分析了太空碎片对系绳生存能力可能造成的影响.双股系绳是由两条平行系绳每隔一段距离打结连接构成,带状系绳可视为一种横截面为矩形的特殊单股系绳.首先对单股系绳被碎片割断的概率进行建模,然后构建双股系绳和带状系绳生存能力函数,最后根据空间碎片环境模型,对比分析具有不同结构、不同尺寸的系绳的生存能力.仿真结果表明,相较于单股系绳,双股系绳和带状系绳的生存能力有明显提高.      

A Partially Orthogonal EnKF approach to atmospheric density estimation using orbital debris

A key requirement for accurate trajectory prediction and space situational awareness is knowledge of how non-conservative forces affect space object motion. These forces vary temporally and spatially, and are driven by the underlying behavior of space weather particularly in Low Earth Orbit (LEO). Existing trajectory prediction algorithms adjust space weather models based on calibration satellite observations. However, lack of sufficient data and mismodeling of non-conservative forces cause inaccuracies in space object motion prediction, especially for uncontrolled debris objects. The uncontrolled nature of debris objects makes them particularly sensitive to the variations in space weather. Our research takes advantage of this behavior by utilizing observations of debris objects to infer the space environment parameters influencing their motion.The hypothesis of this research is that it is possible to utilize debris objects as passive, indirect sensors of the space environment. We focus on estimating atmospheric density and its spatial variability to allow for more precise prediction of LEO object motion. The estimated density is parameterized as a grid of values, distributed by latitude and local sidereal time over a spherical shell encompassing Earth at a fixed altitude of 400 km. The position and velocity of each debris object are also estimated. A Partially Orthogonal Ensemble Kalman Filter (POEnKF) is used for assimilation of space object measurements to estimate density.For performance comparison, the scenario characteristics (number of objects, measurement cadence, etc.) are based on a sensor tasking campaign executed for the High Accuracy Satellite Drag Model project. The POEnKF analysis details spatial comparisons between the true and estimated density fields, and quantifies the improved accuracy in debris object motion predictions due to more accurate drag force models from density estimates. It is shown that there is an advantage to utilizing multiple debris objects instead of just one object. Although the work presented here explores the POEnKF performance when using information from only 16 debris objects, the research vision is to utilize information from all routinely observed debris objects. Overall, the filter demonstrates the ability to estimate density to within a threshold of accuracy dependent on measurement/sensor error. In the case of a geomagnetic storm, the filter is able to track the storm and provide more accurate density estimates than would be achieved using a simple exponential atmospheric density model or MSIS Atmospheric Model (when calm conditions are assumed).     

An adaptive threshold method for improving astrometry of space debris CCD images

Optical survey is a main technique for observing space debris, and precisely measuring the positions of space debris is of great importance. Due to several factors, e.g. the angle object normal to the observer, the shape as well as the attitude of the object, the variations of observed characteristics for low earth orbital space debris are distinct. When we look at optical CCD images of observed objects, the size and brightness are varying, hence it’s difficult to decide the threshold during centroid measurement and precise astrometry. Traditionally the threshold is given empirically and constantly in data reduction, and obviously it’s not suitable for data reduction of space debris. Here we offer a solution to provide the threshold. Our method assumes that the PSF (point spread function) is Gaussian and estimates the signal flux by a directly two-dimensional Gaussian fit, then a cubic spline interpolation is performed to divide each initial pixel into several sub-pixels, at last the threshold is determined by the estimation of signal flux and the sub-pixels above threshold are separated to estimate the centroid. A trail observation of the fast spinning satellite Ajisai is made and the CCD frames are obtained to test our algorithm. The calibration precision of various threshold is obtained through the comparison between the observed equatorial position and the reference one, the latter are obtained from the precise ephemeris of the satellite. The results indicate that our method reduces the total errors of measurements, it works effectively in improving the centering precision of space debris images.     

基于解体事件的空间碎片轨道演化算法研究

针对航天器解体事件所生成的空间碎片的演化过程,进行了数学分析,确定了新生成的空间碎片的速度增量,在该增量作用下碎片轨道会发生变更,本文根据该增量得出了空间碎片在轨道变更后的轨道根数,分析了在大气阻力摄动作用下,空间碎片的数目和轨道分布的演化情况,给出了相关结果,结果表明此算法可行。     

Hypervelocity impacts on HST solar arrays and the debris and meteoroids population

Accurate debris and meteoroid flux models are crucial for the design of manned and unmanned space missions. For the most abundant particle sizes smaller than a few millimetres, knowledge of the populations can only be gained from in situ detectors or the analysis of retrieved space hardware. The measurement of impact flux from exposed surfaces improves with increased surface area and exposure time.A post-flight impact investigation was initiated by the European Space Agency to record and analyse the impact fluxes and any potential resulting damage on the two flexible solar arrays of the Hubble Space Telescope. The arrays were deployed during the first Hubble Space Telescope servicing mission in December 1993 and retrieved in March 2002. They have a total exposed surface area of roughly 120 m², including 42 m² covered with solar cells. This new Hubble post-flight impact study follows a similar activity undertaken after the retrieval of one of the first solar arrays, in 1993. The earlier study provided the first opportunity for a numerical survey of damage to exposed surfaces from more than 600 km altitude, and of impacts from particles larger than 1 mm. The results have proven very valuable in validation of important flux model regimes. The second set of Hubble solar arrays has again provided an unrivalled opportunity to measure the meteoroid and debris environment, now sampled during a long interval in low Earth orbit, and to identify changes in the space debris environment since the previous survey. The retrieved solar array wings exhibit thousands of craters, many of which are visible to the naked eye. A few hundred impacts have completely penetrated the 0.7 mm thick array. The largest impact features are about 7–8 mm in diameter. The cover glass of the solar cells is particularly well suited to the recognition of small impact features by optical and electron microscopy. In this paper, we present the first results of the impact survey. Data upon the abundance of craters of specific measured size ranges are plotted as cumulative flux curves, and compared to the results of model predictions. The most significant change to the particle flux since 1993 is a decrease in the small debris population.     

弹丸撞击防护屏形成碎片云速度特性研究

微流星及空间碎片的高速撞击威胁着长寿命，大尺寸航天器的安全运行，导致其严重的损伤和灾难性的失效，为精确估计微流星及空间碎片主速撞击防护屏产生的碎片对舱壁的损伤，必须确定碎片云速度特性。文章在冲量和能量守恒的基础上，建立了碎片速度性分析模型，研究了碎片云的速度特性，得到了碎片云材料传播及碎片云喷射角随弹丸撞击速度的变化规律。     

Predicting the vulnerability of spacecraft components: Modelling debris impact effects through vulnerable-zones

The space environment around the Earth is populated by more than 130 million objects of 1 mm in size and larger, and future predictions shows that this amount is destined to increase, even if mitigation measures are implemented at a far better rate than today. These objects can hit and damage a spacecraft or its components. It is thus necessary to assess the risk level for a satellite during its mission lifetime. Few software packages perform this analysis, and most of them employ time-consuming ray-tracing methodology, where particles are randomly sampled from relevant distributions. In addition, they tend not to consider the risk associated with the secondary debris clouds. The paper presents the development of a vulnerability assessment model, which relies on a fully statistical procedure: the debris fluxes are directly used combining them with the concept of vulnerable zone, avoiding the random sampling the debris fluxes. A novel methodology is presented to predict damage on internal components. It models the interaction between the components and the secondary debris cloud through basic geometrical operations, considering mutual shielding and shadowing between internal components. The methodologies are tested against state-of-the-art software for relevant test cases, comparing results on external structures and internal components.     

卫星发射对空间碎片环境影响分析

未来航天发射情况直接影响空间碎片环境，必须对其进行合理规划，以维护外空长期可持续发展.利用中国自主建立的空间碎片长期演化模型（SOLEM），结合蒙特卡洛方法，量化分析了空间物体发射数量、发射质量、发射面积等因子对未来空间碎片环境的影响，进一步研究了大型星座造成的未来空间物体碰撞次数和碎片数量的增加.仿真结果可为合理规划未来的航天发射规模提供理论依据.      

空间站轨道高度选择

有许多因素影响空间站轨道高度的选择,包括:任务要求,辐射环境,微流星和空间垃圾,空间站构型,运载器能力,发射窗口以及空间站的补给要求等。从补给运输系统性能考虑,要求空间站的轨道高度低一些,但是从轨道维持出发,希望空间站轨道高度高一些。该文讨论了补给、轨道维持和太阳活动对轨道高度的影响,给出了空间站的最佳运行高度。文中还以载人飞船和航天飞机为运输系统进行了数值计算,得到的主要结论是空间站的最佳运行高度不是常数,而是随着太阳活动情况和补给频繁次数变化。在太阳活动高年期间,空间站轨道高度应提高。当空间站的补给次数增加时,空间站的最佳运行高度应该下降。