不同地磁活动期间地球同步轨道

文章利用地球同步轨道GOES6～12卫星1986-2006年约两个太阳活动周期的高能电子角通量密度数据,分析了在地磁平静期或不同程度磁暴期间高能电子通量的变化,为建立更精确的航天器内部充电动态模型奠定基础。研究结果表明:地球同步轨道高能电子通量随太阳活动周期的变化规律表现为,在太阳活动极大年通量较低,极小年反而较高;高能电子通量随季节的变化规律表现为,在冬至和夏至附近通量达到最低值,春秋分附近达到最高值;Dst指数最小时,对应的高能电子通量一般较低。     

不同地磁活动期间地球同步轨道高能电子通量分析

文章利用地球同步轨道GOES～12卫星1986-2006年约两个太阳活动周期的高能电子角通量密度数据,分析了在地磁平静期或不同程度磁暴期间高能电子通量的变化,为建立更精确的航天器内部充电动态模型奠定基础.研究结果表明:地球同步轨道高能电子通量随太阳活动周期的变化规律表现为,在太阳活动极大年通量较低,极小年反而较高;高能电子通量随季节的变化规律表现为,在冬至和夏至附近通量达到最低值,春秋分附近达到最高值;Dst指数最小时,对应的高能电子通量一般较低.     

地球边缘大气密度时空特性分析

分析了地球边缘大气密度的时变特性及空间分布特性。从物理原因出发并基于NRLMSISE-00和MET-99大气模型计算随不同时空变量变化的相对大气密度,以说明各时空因素对地球边缘大气密度的影响程度。分析结果表明地球边缘大气密度具有较强的时变特性,昼夜和季节均对大气密度产生一定影响。空间特性表现为纬度对大气密度的影响明显,而经度影响微弱。季节不同导致纬度的影响程度不同,纬度因素在夏季和冬季的影响强于春秋。随着纬度的增加,季节对地球边缘大气密度的影响逐渐增强,而当地时间的影响呈减弱趋势。相比较而言,低纬区的昼夜波动强于季节波动,高纬区的季节波动强于昼夜波动。     

中国上空过渡流区大气月平均状态的模拟结果

利用美国海军研究实验室在MSIS系列模型基础上发展的大气模型NRLMSISE-00, 模拟研究了2001年—2013年中国上空过渡流区80～140 km高度的大气状态。对于中国东部和中部过渡流区, 基于模拟数据得到的月平均结果显示大气密度和温度呈现一致的变化趋势, 还表现出与太阳活动的显著关联。谱分析的结果显示, 在90 km高度以上大气密度呈现半年周期变化, 在90～110 km高度范围这种半年变化的幅度随纬度增大。在100 km高度, 上半年出现的半年周期过程中密度的最大、最小值分别出现在3月和6月, 明显超前中、高热层的半年周期过程。     

地球扰动大气模型

本文介绍了地球扰动大气模型和建立该模型的方法。文中详细描述了该模型中,密度变化的各单个分量和风场变化的各单个分量,并给出了它们各自的计算公式。本大气扰动模型中包含了137个大气扰动状态方案,可供返回轨道计算、救生轨道计算和返回控制设计计算分别选用。     

地磁扰动期间低轨目标轨道预报的误差修正

针对地磁扰动期间大气密度变化造成的低轨目标较大的轨道预报误差,提出一种根据POES卫星观测的极光能量注入数据改进短期轨道预报的方法。分析表明CHAMP卫星的沿迹大气密度及轨道衰减与极光能量注入具有较好的相关性。通过线性回归方法,建立轨道半长轴衰减及阻力调制系数的修正公式,并使用修正后的阻力调制系数取代两行元(TLE)中的该系数带入SGP4模型进行位置预报。该方案考虑了外推过程中地磁扰动引起的大气密度响应,能更准确地反映外推过程中大气阻力对轨道的影响。将其应用到2008年CHAMP卫星和国际空间站的轨道预报中,结果表明,半长轴和位置的预报误差可分别降低50%和30%左右。进一步对不同年份、不同轨道高度的目标进行了预报误差修正的分析,验证了该方法的普适性。     

高层大气模型对空间站轨道漂移和寿命的影响分析

本文以轨道摄动分析方法一阶理论为基础，其中大气阻力摄动采用数值积分方法，给出一种可利用各种大气模型进行轨道摄动分析的计算方法，并利用三种高层大气模型（ＣＩＲＡ—７２，ＣＩＲＡ—８６和ＤＴＭ）和三个太阳活动水平（Ｆ１０．７＝１００，１５０和２００）分析比较了大气阻力振动对高度为４００ｋｍ的空间站轨道漂移和寿命的影响，以及估算修正轨道漂移所需的能量。给出的定量分析结果将为空间站或航天飞行器的轨道设计和能量估算提供依据。     

基于精确星光大气折射观测模型的轨道摄动研究

基于星光折射间接敏感地平自主导航方法,改进了现有的大气密度模型和固定高度(25km)的观测模型,建立了自适应连续高度(20km～50km)的星光折射观测模型,并在此观测模型的基础上深入分析研究了影响导航精度的各种轨道摄动力,包括地球形状各阶摄动力、不同高度的大气阻力摄动力、太阳光辐射压力等.提出一种精简的大气模式--静止变标高球面大气,解决了由于高层大气密度的求解复杂,使大气阻力摄动模型不精确的问题,由此取得了较好的导航定位精度.利用UKF和EKF两种不同非线性算法,对考虑各种摄动力后的导航定位精度进行全面的计算机仿真实验,并就仿真结果进行了误差分析,同时研究了各种有关参数对导航精度的影响.     

近地卫星简化轨道预报方法研究

基于航天器轨道动力学的基本原理对轨道动力学模型简化,以保持模型精度,同时降低计算负担.某些简化的天文学模型包括岁差、章动、恒星时模型等.针对轨道坐标系统、轨道动力学模型等,利用引力加速度估计函数(GAAF)代替传统的计算量过大的球谐引力场函数模型.其它简化技术包括日月解析星历模型代替DE405模型、简化的大气密度模型以及用于轨道积分的Bulirsh-stoer·积分器等,可进一步减少计算负担.算例表明,本方法作为一种高精度的轨道积分器,可用于星上计算或者轨道蒙特卡罗仿真.计算结果表明其计算时间远小于全尺度模型的1%,而精度相当.算法由CHAMP卫星的真实星上GPS观测数据验证.     

太阳同步轨道空间粒子辐射剂量探测与研究

利用“实践六号”A星3台辐射剂量仪2004年9月至2007年1月的在轨探测数据,对太阳活动低年,该卫星所处太阳同步轨道不同方向、不同屏蔽厚度下的辐射剂量水平及辐射剂量增长变化特性进行了分析和讨论。结果表明：空间环境扰动会导致日辐射剂量率较之空间环境宁静期有1个数量级的增长,这种增长一般滞后于空间环境扰动事件1～2天,持续时间为十几天至数十天;太阳活动低年,平均日辐射剂量率小于5rad天,与地球辐射带模型计算结果大致相当。这些结果对卫星的在轨安全保障和抗辐射加固设计都具有重要意义。     

分布式卫星轨道构形的大气摄动分析及修正方法

在低轨道运行的分布式卫星受大气摄动的影响,其轨道构形很快遭到破坏,环绕卫星的相对运动轨迹中心不断发生漂移。本研究的目的在于建立一种降低大气摄动对分布式卫星轨道构形影响的补偿方法,以使分布式卫星的轨道构形能够更好地自然维持。研究基于考虑高度变化和太阳周日变化的大气密度模型,得出分布式卫星不同初始相位环绕卫星的长半轴摄动方程。并提出一种补偿大气摄动影响的长半轴修正方法。仿真结果显示,采用此大气摄动补偿方法能够在给定时间内大大降低大气摄动的影响,从而显著提高分布式卫星轨道构形的自然维持能力。     

An evaluation of Jacchia and MSIS 90 atmospheric models with CBERS data

“CBERS” (China Brazil Earth Resources Satellite) was put into orbit on 14^th of October 1999. The atmospheric drag is the major non-gravitational perturbation affecting the control of ground track. As accuracy requirements increase, greater reliance is placed on the empirical techniques. During the initial operational phase of CBERS, the solar activity is almost at its peak. This phenomenon has provided an opportunity to carry out an evaluation of the atmospheric density models. The study puts emphasis on two commonly used atmospheric density models viz. Jacchia and Mass Spectrometer Incoherent Scatter (MSIS). The analysis is based on the decay histories of CBERS to choose the accurate density model. Density models used for orbit propagation are usually derived empirically from actual flight data. Brief synopsis of some of the models is presented along with some of the density tables and orbit solutions of CBERS. Typical plots of density are presented. The study indicated that the drag estimation is relatively precise using MSIS based models. Among them MSIS-90 density model is observed to be a better compromise in terms of accuracy, flexibility and computational aspects. The analysis would be useful in mitigating the impact of solar activity on orbit prediction and maintenance.     

一种地球紫外中心快速提取方法

针对地球紫外中心指向高精度提取问题,提出一种基于梯度统计的快速、低存储需求的紫外地平圆盘中心提取算法。首先,考虑到星载计算机与高内存消耗的滤波算法、实时导航需求的冲突,采用结合Sobel边缘算子与局部二值模式（LBP）算子的改进边缘快速提取方法,有效剔除背景噪声并准确提取地球紫外临边特征。然后对非完整临边边缘采用最小二乘拟合得到中心的精确位置。实验结果表明,该方法抗噪声,可实现亚像素级地球中心提取,并显著降低存储需求和计算时间开销。对于1046×746的图像,该算法需要的存储空间仅为1100kByte,运算时间在20ms内,满足自主导航的需求。     

预估太阳电池阵在轨最高温度的经验公式

针对距地高度为300～1000 km的太阳同步轨道,分析了一维转动展开式太阳电池阵的温度影响因素,给出了计算太阳电池阵最高温度的经验公式。此公式同样适用于低地球轨道的卫星,且计算结果可包络卫星热控状态变化带来的影响,在太阳电池阵设计阶段可为电源分系统提供太阳电池阵最高温度的初步数据。     

Cosmic-ray propagation through the disturbed heliosphere

The stationary, source-free, one-dimensional solution to the cosmic-ray heliospheric transport equation was solved using Runge–Kutta methods. These calculations have been applied to the treatment of Forbush decreases and the evolution of the solar cycle through the heliosphere. Forbush decreases calculated by this method produce spectra that are in good agreement with both the theoretical heliocentric potential and the potential as determined from neutron monitor data. During solar modulation, the effect of the travelling disturbance is to flatten the gradient of the cosmic ray flux in the inner heliosphere at solar minimum and steepen the gradient at solar maximum, as the incoming cosmic rays see solar activity up to a year earlier. The resulting spectra at the earth's orbit are in good agreement with the heliocentric approximation. The heliocentric approximation is an accurate and convenient method of accounting for solar modulation.     

The heliospheric magnetic field and its relation to the temperature,density, and velocity of solar plasma: Experimental evidence

Using daily and hourly data on solar plasma parameters at the Ulysses spacecraft orbit and at 1 AU it is demonstrated that there is a simple relationship between plasma temperature and density with the heliospheric magnetic field (HMF). A mathematical expression connecting HMF with plasma temperature and density is suggested. Correlation coefficients and regression equations for measured and calculated magnetic fields are presented for the 1990–2009 period according to Ulysses spacecraft data and for 2003–2010 at 1 AU (OMNI database). The roles played by density, temperature, and high-speed solar wind streams in forming the magnetic-field peaks are demonstrated using hourly data of OMNI2 and Ulysses.     

The Earth's Magnetosphere Response to Solar Wind Events according to the INTERBALL Project Data

The results of studying the interaction of two types of the solar wind (magnetic clouds and solar wind of extremely low density) with the Earth's magnetosphere are discussed. This study is based of the INTERBALL space project measurements and on the other ground-based and space observations. For moderate variations of the solar wind and interplanetary magnetic field (IMF) parameters, the response of the magnetosphere is similar to its response to similar changes in the absence of magnetic clouds and depends on a previous history of IMF variations. Extremely large density variations on the interplanetary shocks, and on leading and trailing edges of the clouds result in a strong deformation of the magnetosphere, in large-scale motion of the geomagnetic tail, and in the development of magnetic substorms and storms. The important consequences of these processes are: (1) the observation of regions of the magnetosphere and its boundaries at great distances from the average location; (2) density and temperature variations in the outer regions of the magnetosphere; (3) multiple crossings of geomagnetic tail boundaries by a satellite; and (4) bursty fluxes of electrons and ions in the magnetotail, auroral region, and the polar cap. Several polar activations and substorms can develop during a single magnetic cloud arrival; a greater number of these events are accompanied, as a rule, by the development of a stronger magnetic storm. A gradual, but very strong, decrease of the solar wind density on May 10–12, 1999, did not cause noticeable change of geomagnetic indices, though it resulted in considerable expansion of the magnetosphere.     

Vertical structure of the magnetic field in active regions of the sun at coronal heights

High-altitude measurements and magnetic field distribution in active regions are always an important problem to verify existing models at heights from units to hundreds of Mm. Optical methods of analysis of the magnetic structure work well only at the photosphere level. At the same time, the progress of radio astronomy methods of analyzing the solar radio emission [1] and of the theory of solar plasma radiation [2] facilitates introduction of new methods for measuring the magnetic field strength at various altitudes in the solar atmosphere. In this paper we use multi-wave observations of polarization of the radio emission of active regions in the microwave range together with precise magnetosphere data in order to develop a method of projection for measuring the height distribution of the magnetic field.