太阳的空间观测始于20世纪50年代,到目前为止,比较著名的太阳探测卫星有OSO, Hinotori, Yohkoh, SOHO, TRACE, Ulysses[1], RHESSI, HINODE, STEREO, SDO和IRIS等。这些卫星上搭载各种设备,通过多波段对太阳进行成像和原位探测,极大地促进了人们对太阳爆发活动过程和物理机制的理解。然而,上述卫星基本上运行在围绕地球转动的轨道上,在太阳的黄道面附近以地球单视线方向观测太阳。唯一一颗脱离地球进入太阳极轨轨道运行的卫星Ulysses,上面没有搭载成像设备,限制了对太阳极区的观测。太阳极区是太阳活动和磁场变化的关键区域,通过对极区的观测和研究,可以更好地了解太阳活动的周期性变化、太阳风的性质、太阳磁场的演化等等。此外,太阳大气中的活动具有非常复杂的三维结构,单视线方向的磁场观测带有与生俱来的缺陷,理论和数值模拟研究已经发展到更为真实的三维。这些都极其需要多视角立体观测来探测其中的三维结构,以便准确把握其中真实的物理过程。文[2]首次提出了利用3颗卫星从不同方向对太阳进行立体观测,是目前精确测量太阳矢量磁场比较可行的方法。通过进一步模拟发现,30 °的夹角是立体磁场观测的最佳角度[3]。因此,为实现立体观测,3颗卫星中至少1颗应部署在倾角大于30 °的太阳极轨轨道上。另一方面,太阳活动的能量来源于磁场自由能的突然释放,磁重联在这一过程中起核心作用。目前对磁重联机制的理解依然比较浅薄,原因之一是缺乏对磁重联过程中产生的电流片结构的原位探测。这需要卫星抵近太阳几倍半径处,穿越电流片结构,直接原位探测电流片中的等离子体的特性。太阳的抵近探测还能提供前所未有的超高分辨率太阳活动图像,将帮助人们解决长期令人迷惑的“日冕加热”问题。因此,太阳的深空探测,包括对太阳极区的极轨探测和对电流片结构以及太阳活动精细结构的抵近探测,成为当前急需发展的空间项目。所幸的是,目前我国的深空探测技术已经克服了日地距离带来的障碍。中国太阳物理学家也提出了自己的“极轨探测”[4]和“抵近探测”计划[5],并且已经在进行关键技术的研究[6]。
将卫星直接发射到深空极其困难,需要的能量远远超过目前任何火箭所具备的能力。国际上已发射的有限几个偏离黄道面的太阳探测卫星无不借助行星的引力。例如首个太阳极轨轨道卫星Ulysses,于1990年10月发射,经过14个月抵达木星。在木星强大的引力作用下,Ulysses卫星轨道发生偏转并抬升到80.2 °,最终进入周期为6.2年的太阳极轨轨道上运行。SolO (Solar Orbiter) [7]是2020年2月发射的抵近探测太阳极区的卫星,借助地球和金星的引力,SolO需要大约3.5年才能使轨道倾角最终抬升到24 °,它的近日点为0.28 AU,远日点为0.91 AU。2018年发射的PSP (Parker Solar Probe) [8]是首颗抵近太阳的探测卫星,它借助金星和地球的引力入轨,最近将抵达距太阳表面9.86倍太阳半径处。2006年发射的STEREO的两颗卫星则是在月球的引力辅助下进入绕日轨道运行。除了借助行星的引力,欧洲空间局(European Space Agency, ESA) 还曾提出利用太阳帆实现太阳极轨探测的概念[9]。太阳帆通常由一个大型薄膜组成,在航天器发射后展开,利用太阳光对其施加的压力推动自身运动。太阳帆虽然推力较小,但是它能够持续提供推力。这一概念早在100多年前就已经提出,首要的困难是太阳帆的制作。太阳帆需要采用轻薄、耐高温、高透明度和高强度的材料,并且制作面积需要足够大。其次,需采用非常精确的控制系统,以确保太阳帆能够准确地朝着目标方向推进。这些要求在技术上带来了极大的挑战。日本航天局曾在2010年成功实施第1次太阳帆飞行任务IKAROS (Interplanetary Kite-craft Accelerated by Radiation of the Sun) [10],使航天器从地球飞到金星。作为技术验证,证明了太阳帆推动卫星确实可行。但是由于太阳光压产生的推力很小,依靠太阳帆的动力提高卫星轨道需要的周期非常长(可达5年)。综上所述,目前太阳极轨和抵近轨道成熟可行的技术方案是借助行星或月球的引力,太阳帆提供了另一种可能。但是这两种方案入轨消耗的时间都很长,入轨后卫星真正工作的时间在整个任务中可能仅约占一半。
小行星是太阳系内围绕太阳运行的一类行星形态的小天体,一般以岩石、金属或者冰体的形态存在,上面没有大气。小行星呈现不规则的形状,它们的尺寸范围从数米到数百千米不等,其中最大的小行星直径超过900 km。根据美国罗威尔天文台提供的小行星数据(网址: http://asteroid.lowell.edu/main/astorb/),目前观测到的小行星数量已超过百万。大部分小行星位于火星和木星之间的小行星带,距离太阳2~4 AU。小行星通常是太阳系内最古老的物体之一,保留了很多关于太阳系形成和演化的重要信息。小行星上含有的稀有矿物质未来可能成为人类太空开采的目标,用来替代地球上稀少的或者已经耗尽的物质。此外,一些近地小行星具有闯入地球大气层的潜在危险,威胁地球上所有的生命。为此,美国航空航天局(National Aeronautics and Space Administration, NASA) 和欧洲空间局分别或者联合制定了小行星防御计划,把具有威胁的小行星拖离地球。2022年9月美国航空航天局实施了名为“双小行星重定向测试”(Double Asteroid Redirection Test, DART) [11]的任务,成功利用卫星正面撞击小行星,使小行星的轨道发生了轻微的改变。因此,研究小行星的形成、结构和轨道特征,对于了解人类所在的太阳系的形成历史以及地球上水和生命的起源问题、小行星资源的利用以及为人类提供防卫措施都具有极其重要的意义。
本文提出一种新颖的太阳极轨和抵近轨道的部署方案,即搭乘某些特殊的小行星直接将卫星带入轨道。该方案能够大大节省轨道部署的时间和经济成本,同时也是首次实现对小行星资源的利用。
1 小行星统计为了实现科学目标,太阳极轨卫星需要部署在倾角至少30 °以上且偏心率较小的轨道。同时考虑到搭乘方便,小行星的轨道与黄道面的交点需要在地球轨道附近。为此,我们选取目前已观测到的所有轨道倾角大于30 °的小行星作为样本,统计它们的轨道倾角、半长轴以及偏心率的分布,样本总数为9 994颗。图 1显示了小行星的轨道倾角和半长轴的分布。由图 1 (a)可见,越靠近黄道面(倾角越小),小行星的数量越多;反之,随着偏离黄道面,小行星的数量越来越少,不过倾角大于60 °的小行星仍有163颗。而小行星轨道的半长轴分布有所不同,主要分布在2.7 AU附近(见图 1(b))。除去这个区域,在1~3.5 AU范围内,小行星轨道的半长轴基本呈现较为均匀的分布。1 AU以内的距离也分布一定数量的小行星。图 2显示了小行星轨道偏心率的分布。由图 2可见,大部分小行星的轨道偏心率集中在0.4以下,甚至有些在0~0.1之间,这说明有不少小行星的轨道非常接近圆形。以上的统计表明,在地球轨道附近存在倾角比较多的偏心率较低的小行星,能够满足太阳极轨卫星的搭乘需求。
|
| 图 1 小行星轨道倾角和半长轴的分布 Fig. 1 Inclination and semi-major axis distribution of asteroid orbits |
|
| 图 2 小行星轨道偏心率的分布 Fig. 2 Eccentricity distribution of asteroid orbits |
为满足太阳抵近探测的轨道需求,我们统计了小行星轨道近日点距离的分布(见图 3)。如图 3,小行星的数量随着近日点距离增加而不断增加。近日点距离在0.55 AU以内的小行星数量总计1 883颗,在0.1 AU距离以内仍存在少量小行星。近日点距离最小的小行星是2005 HC4,距离太阳仅约15倍太阳半径。
|
| 图 3 小行星轨道近日点距离的分布 Fig. 3 Perihelion distance distribution of asteroids orbits |
有些小行星的轨道既具有较高的倾角,也具有较小的近日点距离,可以同时满足太阳极轨探测和抵近探测的轨道需求。图 4显示了轨道倾角大于30 °和近日点距离小于0.5 AU的小行星散布图。这类小行星总计263颗。由图 4可见,在较高倾角和较小近日点距离范围内也存在小行星。
|
| 图 4 小行星倾角和近日点距离的散点图 Fig. 4 Scatter plot of asteroids inclination and perihelion distance |
由于小行星的尺寸微小,对其直径的准确测量较为困难,需要投入大量的观测时间和资源。在小行星数据库中,大部分小行星没有直径的测量数据。在倾角高于30 °的小行星样本中,仅有15颗小行星具有直径的测量数据,它们的直径大约在几十千米左右。随着观测设备能力的提升,未来有望不断发现直径更小的小行星。
根据以上对小行星轨道参数的统计,我们可以发现存在各种各样的小行星轨道,就像一个立体交通系统一样运行在太阳系内。如果能够通过某种方式让卫星搭乘小行星,那么找到一个合适的、满足轨道要求的目标小行星是不存在问题的。特别是对于极轨探测,存在偏心率较小、倾角较高的轨道,卫星进入这种接近圆形的轨道对太阳进行极轨探测极其有利。为此,我们提出了卫星搭乘小行星入轨的方案。
2 太阳深空探测轨道的部署方案具体实施方案是在卫星的前端安装一个套星圆框,框内布置具有高弹性和高韧性的特殊材料制作的大网,类似于太阳帆。这张网用来网住小行星,其直径需要超过小行星的直径。首先,根据任务需求,选择满足轨道参数的大小适当的目标小行星,直径在50~100 m之间最佳。选定小行星后,对其进行精密定轨观测,高精度确定小行星的轨道参数。目前,使用1 m光学望远镜,配合相互逼近技术,对小行星位置的测量精度可以达到6毫角秒[12]。然后,当小行星靠近地球时,通过常规方式将卫星发射到地球轨道附近的绕日轨道上。该轨道在黄道面上,与小行星的轨道相交。卫星通过自主导航与小行星不断接近,在两者轨道交叉处相遇。小行星穿过套星圆框,拉拽弹力网将卫星拖入其轨道。最后,卫星脱离小行星,自主调整姿势,实现预定的探测任务。整个入轨过程仅相当于把卫星发射到地球轨道附近所需的时间(几个月),与小行星交汇后,可在瞬间被拖入其轨道。因此,相比其他方案,该方案可以缩短2~5年的时间。
图 5给出了几个候选小行星的运行轨道。其中粉红色、淡紫色和深蓝色的轨道分别表示水星、金星和地球的轨道。其他3条灰色的轨道则是小行星2003 QQ47, 2005 HC4和Moshup的轨道,分别代表可以满足太阳极轨、抵近探测以及同时满足两者轨道要求的目标小行星。以小行星2003 QQ47为例,它的轨道倾角为62.19 °,半长轴为1.09 AU,偏心率为0.19,轨道运行速度约24.5 km/s。将卫星发射到与地球近似的轨道上。在小行星与卫星交汇时,小行星通过拉拽弹力网将卫星带入其轨道。假设卫星的运行速度与小行星相同,它们之间的夹角近似为小行星的轨道倾角,即大约60 °。那么卫星的动量改变为|I| = m va-m vs =mv,其中va和vs分别为小行星和卫星的速度矢量,v是速度的大小,m是卫星的质量。根据冲量定理,|I| = |F| t,以我国首颗综合性太阳探测卫星“夸父一号”的质量约859 kg作为卫星的质量,如果小行星网作用时间为1 s,以此推算,卫星受到的平均拉力约为2×107 N,换算成地球上的质量大约是2×103 kg。
|
| 图 5 小行星2003 QQ47, 2005 HC4和Moshup的运行轨道 Fig. 5 Orbits of asteroid 2003 QQ47, 2005 HC4 and Moshup |
卫星受到的瞬时作用力非常巨大,这给弹力网的材料特性带来巨大的挑战。我们可以考虑通过增加作用时间,比如增加弹力网的弹性等措施来减小弹力网的拉力。当然,以上仅是简单的估算。考虑到更精细的过程,在小行星运行方向上,弹性张力从0逐渐增加到极大值,然后再逐渐减小至0。在这个过程中,卫星一直加速。因此需要设计某种机制,在适当的时候脱离弹力网。
类似于小行星2003 QQ47这样的轨道倾角在60 °以上,半长轴在1.1 AU以内的高倾角轨道的小行星目前总共有7颗。表 1列出了这7颗小行星的轨道参数。如果轨道倾角扩大到大于30 °,则符合要求的小行星可达200多颗。对于抵近探测来说,交汇时小行星与卫星的运行方向大致相同,实现的难度相比极轨探测要低得多。考虑到弹力网尺寸的限制,小行星的直径最好小于100 m。文[13]统计发现,直径50~140 m、20~50 m、10~20 m以及10 m以下的近地小行星新增发现数逐年增长,已发现直径10 m以上的近地小行星总数仅占理论预测值的不到0.1%[14]。可见,满足条件的目标小行星将越来越多。
| Name | Semi-major axis/AU | Inclination /(°) | Eccentricity |
| Izhdubar | 1.00 | 63.46 | 0.27 |
| 2003 QQ47 | 1.09 | 62.10 | 0.19 |
| 2003 UL9 | 1.01 | 62.94 | 0.23 |
| 2015 XG261 | 0.99 | 65.44 | 0.20 |
| 2016 JP5 | 0.85 | 65.21 | 0.31 |
| 2017 MH | 1.05 | 65.02 | 0.13 |
| 2017 TF4 | 1.09 | 60.52 | 0.24 |
根据调研,这种实现极轨和抵近轨道的概念由我们首次提出。该方案将太阳的深空探测与小行星的探测相结合。与SolO 24 °的轨道倾角和0.53的偏心率相比,该方案可以提供半径约为1 AU的更高倾角的近圆形轨道,这对太阳极轨探测更具优势。有些轨道可以同时进行太阳的极轨探测和抵近探测。该方案入轨时间短,比借助行星引力或者太阳帆的方案可缩短5年以上。此外,该方案还具有灵活性,即便卫星和小行星相遇失败,卫星也可以在其原有轨道继续执行任务,等待下次相遇的机会。一旦掌握了该方案的实施技术和数据,小行星可以重复搭乘。该方案的执行,还能为拖离近地危险小行星积累经验。直径在100~200 m的小行星质量大约是1亿吨左右,比卫星的质量大8个量级,因此卫星对小行星轨道的改变微乎其微,不会对地球引入新的威胁。
该方案主要的难点如下:(1) 弹力网的材料和制作。弹力网需要足够的韧性和轻量化。最近国内刚刚突破的国产化高性能纤维可以用来制作直径12 mm的细绳,能够拉动20多吨的重物,有望满足这一需求。(2) 小行星和卫星轨道的精准测控。精度要高于1 m,否则可能导致卫星和小行星相撞,从而使卫星报废。(3) 卫星及载荷必须具备坚固的结构和极强的抗冲击能力,以承受瞬间巨大的冲击力。同时,需要设计一个将弹力网装入体积有限的容器内并弹出的装置,以及最终脱离小行星的装置。(4) 卫星姿态的自主调整。在卫星受到剧烈冲击后,调整姿态,以进入工作模式。
该方案的核心是搭乘小行星入轨,不排除有比利用弹力网更好的搭乘方法。遵循由简至难的原则,建议整个任务可以分为3步:(1) 发射一颗实验小卫星,对方案进行验证;(2) 发射抵近探测卫星;(3) 发射极轨探测卫星。
3 总结太阳的深空探测是多年来人们一直追求的目标。由于技术上的困难以及经济上的巨大投入,仅有少数项目得以实施。因此,低成本和高效的轨道设计显得极为关键。本文对目前已发现的小行星的轨道进行了统计和分析,发现各种各样的小行星在太阳系内组成了一个立体交通系统。据此,本文开创性地提出一种搭乘特殊小行星实现太阳深空探测的入轨方案。该方案能够大大减少任务执行的时间,节省卫星携带的燃料,从而降低经济成本。方案的实施不仅可以使我国的太阳物理研究跨入国际先进行列,而且能够促进小行星的探测、卫星的测控和材料科学的发展,同时也将首次实现对小行星资源的开发利用。
| [1] | WENZEL K P, MARSDEN R G, PAGE D E, et al. The Ulysses mission[J]. Astronomy and Astrophysics Supplement, 1992, 92(2): 207–219. |
| [2] |
季海生, 汪毓明, 汪景琇. 太阳的立体观测[J]. 中国科学: 物理学力学天文学, 2019, 49(5): 059605 JI H S, WANG Y M, WANG J X. Stereoscopic observations of the Sun[J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2019, 49(5): 059605. |
| [3] |
周茹芸, 汪毓明, 宿英娜, 等. 利用双视角消除太阳矢量磁场观测中的180°不确定性[J]. 天文学报, 2021, 62(4): 104–115 ZHOU R Y, WANG Y M, SU Y N, et al. Using observations of solar vector magnetic field from dual view points to remove the 180° ambiguity[J]. Acta Astronomica Sinica, 2021, 62(4): 104–115. |
| [4] |
邓元勇, 周桂萍, 代树武, 等. 太阳极轨天文台[J]. 科学通报, 2023, 68(4): 298–308 DENG Y Y, ZHOU G P, DAI S W, et al. Solar polar-orbit observatory[J]. Chinese Science Bulletin, 2023, 68(4): 298–308. |
| [5] | LIN J, HUANG S J, LI Y, et al. In situ detection of the solar eruption: lay a finger on the Sun[J]. Chinese Astronomy and Astrophysics, 2021, 45: 301–351. DOI: 10.1016/j.chinastron.2021.08.003 |
| [6] |
黄善杰, 种晓宇, 林隽, 等. 太阳爆发抵近探测器的热防护设计方案研究[J]. 天文研究与技术, 2021, 18(1): 87–100 HUANG S J, ZHONG X Y, LIN J, et al. Study on design scheme of thermal protection of probe for in situ measurements of solar eruption[J]. Astronomical Research & Technology, 2021, 18(1): 87–100. |
| [7] | MÜLLER D, CYR O C, ZOUGANELIS I, et al. The Solar Orbiter mission science overview[J]. Astronomy and Astrophysics, 2020, 642: A1. DOI: 10.1051/0004-6361/202038467 |
| [8] | FOX N J, VELLI M C, BALE S D, et al. The solar probe plus mission: humanity's first visit to our star[J]. Space Science Reviews, 2016, 204(1): 7–48. |
| [9] | APPOURCHAUX T, LIEWER P, WATT M, et al. POLAR investigation of the Sun-POLARIS[J]. Experimental Astronomy, 2009, 23: 1079–1117. DOI: 10.1007/s10686-008-9107-8 |
| [10] | TSUDA Y, MORI O, FUNASE R, et al. Flight status of IKAROS deep space solar sail demonstrator[J]. Acta Astronautica, 2011, 69(9/10): 833–840. |
| [11] | RIVKIN A S, CHABOT N L, STICKLE A M, et al. The double asteroid redirection test (DART): planetary defense investigations and requirements[J]. The Planetary Science Journal, 2021, 2: 173. |
| [12] |
覃永贵, 彭青玉. 相互逼近技术用于小行星702 (Alauda) 的高精度位置测量试验[J]. 天文研究与技术, 2022, 19(2): 118–124 QIAN Y G, PENG Q Y. The experiment for the application of mutual approximation technique to high-precision positional measurement of asteroid 702 (Alauda)[J]. Astronomical Research & Technology, 2022, 19(2): 118–124. |
| [13] |
胡寿村, 赵海斌, 季江徽. 近地小行星发现数目与发现场景的统计分析[J]. 天文学报, 2022, 63(5): 3–25 HU S C, ZHAO H B, JI J H. Statistical analysis of discoveries and discovery scenarios of near-earth asteroids[J]. Acta Astronomica Sinica, 2022, 63(5): 3–25. |
| [14] | HARRIS A W, CHODAS P W. The population of near-earth asteroids revisited and updated[J]. Icarus, 2021, 365: 114452. |



