《万方数据-数字化期刊群》全文上网期刊
CNKI《中国学术期刊(网络版)》全文收录期刊
《中文科技期刊数据库》(维普网)全文收录期刊
超星期刊域出版平台、博看网全文收录期刊
日本科学技术振兴机构数据库收录

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响

呼延奇 蔡震波

呼延奇, 蔡震波. 地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响[J]. 航天器环境工程, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
引用本文: 呼延奇, 蔡震波. 地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响[J]. 航天器环境工程, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
HU Y Q, CAI Z B. Influence of slot region dynamics of Earth radiation belt on the radiation effect of medium-orbit satellites[J]. Spacecraft Environment Engineering, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
Citation: HU Y Q, CAI Z B. Influence of slot region dynamics of Earth radiation belt on the radiation effect of medium-orbit satellites[J]. Spacecraft Environment Engineering, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014

地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响

doi: 10.12126/see.2021.03.014
详细信息
    作者简介:

    呼延奇(1982—),男,博士学位,高级工程师,研究方向为航天器空间环境工程设计。E-mail: zzuhuyan@163.com

  • 中图分类号: P353; V524.3

Influence of slot region dynamics of Earth radiation belt on the radiation effect of medium-orbit satellites

  • 摘要: 中轨卫星运行于地球辐射带槽区,而槽区粒子辐射环境可能存在显著涨落,增加卫星抗辐射设计输入的不确定性。文章利用典型地球辐射带模型,对中轨卫星累积性辐射效应的主要来源进行深入分析,再结合槽区粒子辐射环境动态变化特征,初步量化分析其对中轨卫星遭遇辐射效应的影响。结果表明:槽区粒子辐射环境的动态变化对星表材料及太阳电池辐射损伤的附加影响较小;槽区质子填充事件对8000 km以上高度轨道的电离总剂量有明显影响(但此类事件遭遇概率很低);槽区电子填充事件使10 000 km以上高度轨道的电离总剂量明显增大,这点必须在相应的卫星抗辐射设计要求中予以考虑。
  • 图  1  地球辐射带捕获粒子通量随L值的变化曲线

    Figure  1.  Distribution of fluxes of trapped particles in Earth’s radiation belt vs L value

    图  2  1 mm铝屏蔽下不同能量的捕获电子和捕获质子对总剂量的贡献

    Figure  2.  Contribution of electrons (left) and protons (right) of different energies to the total dose with 1 mm aluminum shield

    图  3  3 mm铝屏蔽下不同能量的捕获电子和捕获质子对总剂量的贡献

    Figure  3.  Contribution of electrons (left) and protons (right) of different energies to the total dose with 3 mm aluminum shield

    图  4  不同轨道上捕获电子和质子对太阳电池辐射年累积等效1 MeV电子通量

    Figure  4.  Annual cumulative radiation on solar cell equivalent to 1 MeV electron flux for electrons (left) and protons (right) on different orbits

    图  5  常规厚度和5倍常规厚度玻璃盖片屏蔽下太阳电池辐射损伤来源分析

    Figure  5.  Particle source of radiation damage on solar cell under glass cover of normal thickness(left) and 5× thickness

    表  1  不同高度轨道卫星总剂量来源分析

    Table  1.   Particle source for total dose on satellite at different altitudes of orbits

    轨道高度/
    km
    不同屏蔽状态下不同能量的捕获电子和质子对总剂量的贡献/%(粒子能量/MeV)粒子类型
    0.04 mm铝屏蔽状态1 mm铝屏蔽状态3 mm铝屏蔽状态10 mm铝屏蔽状态
    3000约3%(2~20)约47%(15~100)约76%(25~300)约98%(50~400)质子
    约97%(0.1~0.6)约53%(0.7~3)约24%(1.7~4)约2%(0.3~4.5)电子
    4000约8%(2~20)约65%(15~100)约84%(25~300)约98%(50~400)质子
    约92%(0.1~0.6)约35%(0.7~3)约16%(1.7~4)约2%(0.3~4.5)电子
    5000约17%(2~20)约74%(15~100)约90%(25~300)约94%(50~400)质子
    约83%(0.1~0.6)约23%(0.7~3)约10%(0.4~4)约6%(0.3~4.5)电子
    6000约30%(2~20)约80%(15~100)约92%(25~300)约87%(50~400)质子
    约70%(0.1~0.6)约20%(0.7~2)约8%(0.15~4)约13%(0.1~1)电子
    7000约44%(2~20)约82%(15~60)约92%(25~300)约77%(50~400)质子
    约56%(0.1~0.6)约18%(0.7~2)约8%(0.15~4)约13%(0.1~1)电子
    8000约56%(2~20)约74%(15~60)约86%(25~200)约58%(50~200)质子
    约44%(0.1~0.6)约26%(0.7~2)约14%(0.15~4)约42%(0.1~2)电子
    9000约66%(2~10)约50%(15~60)约64%(25~100)约36%(50~200)质子
    约34%(0.1~0.6)约50%(0.7~3)约36%(0.1~4)约64%(0.1~2)电子
    10 000约77%(2~10)约23%(15~60)约18%(25~60)约17%(50~100)质子
    约23%(0.1~1)约77%(0.7~4)约82%(0.1~4)约83%(0.1~4)电子
    11 000约86%(2~8)约7%(15~40)约3%(25~40)约6%(50~100)质子
    约14%(0.1~1)约93%(0.7~4)约97%(0.1~4)约94%(0.06~4)电子
    12 000约93%(2~8)约2%(15~40)可忽略可忽略质子
    约7%(0.1~1)约98%(0.7~4)约100%(1.7~4)约100%(0.06~6)电子
    13 000约92%(2~8)可忽略可忽略可忽略质子
    约8%(0.1~1)约100%(0.7~5)约100%(1.7~5)约100%(0.06~6)电子
    14 000约88%(2~8)可忽略可忽略可忽略质子
    约12%(0.1~1)约100%(0.7~5)约100%(1.7~5)约100%(0.06~6)电子
    15 000约81%(2~8)可忽略可忽略可忽略质子
    约19%(0.1~1)约100%(0.7~5)约100%(1.7~5)约100%(0.06~6)电子
    下载: 导出CSV

    表  2  不同高度轨道卫星太阳电池辐射损伤来源分析

    Table  2.   Particle source for solar cell radiation damage on different altitudes of satellite orbits

    0°倾角轨道
    高度/km
    不同屏蔽状态下不同能量的捕获粒子
    对辐射损伤的贡献/%(粒子能量/MeV)
    粒子类型
    常规玻璃盖片5倍常规厚度玻璃盖片
    3000约99%(2~100)约99%(8~100)质子
    约1%(0.2~1)约1%(0.45~2.5)电子
    6000约99.9%(2~100)约99.8%(8~100)质子
    约0.1%(0.2~1)约0.2%(0.45~2.5)电子
    9000约99.97%(2~100)约99.8%(8~250)质子
    约0.03%(0.2~1)约0.2%(0.45~2.5)电子
    12 000约99.9%(2~100)约88%(8~250)质子
    约0.1%(0.2~1)约12%(0.45~5)电子
    15 000约99%(2~6)约10%(8~20)质子
    约1%(0.2~4)约90%(0.45~5)电子
    下载: 导出CSV

    表  3  槽区动态变化对6000 km×0°轨道总剂量影响

    Table  3.   Total dose variations on 6000 km×0° orbit due to slot region dynamics

    屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
    AE8MAXAE9
    99%置信度值
    CRRES扰动值/
    SLOT模型
    99%置信度值
    AP8MAXAP9
    99%置信度值
    CRRES
    扰动值
    0.04 mm
    铝屏蔽
    0.1~0.6 MeV(对总剂量贡献70%)2~20 MeV(对总剂量贡献30%)总剂量
    增加0%
    2.26×1081.58×108无数据3.72×1066.66×1066.37×106
    1 mm
    铝屏蔽
    0.7~2 MeV(对总剂量贡献20%)15~100 MeV(对总剂量贡献80%)总剂量
    增加1%
    1.54×1062.54×106无数据4.88×1045.73×1042.30×104
    3 mm
    铝屏蔽
    0.15~4 MeV(对总剂量贡献8%)25~300 MeV(对总剂量贡献92%)总剂量
    增加23%
    1.64×1089.11×107无数据9.21×1033.32×1045.82×103
    10 mm
    铝屏蔽
    0.1~1 MeV(对总剂量贡献13%)50~400 MeV(对总剂量贡献87%)总剂量
    增加73%
    2.31×1081.62×108无数据2.40×1032.27×1041.31×103
    下载: 导出CSV

    表  4  槽区动态变化对8000 km×0°轨道总剂量影响

    Table  4.   Total dose variations on 8000 km×0° orbit due to slot region dynamic changes

    屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
    AE8MAXAE9
    99%置信度值
    CRRES扰动值/
    SLOT模型
    99%置信度值
    AP8MAXAP9
    99%置信度值
    CRRES
    扰动值
    0.04 mm
    铝屏蔽
    0.1~0.6 MeV(对总剂量贡献44%)2~20 MeV(对总剂量贡献56%)总剂量
    增加25%
    1.63×1086.01×1071.83×1095.65×1061.06×1076.86×106
    1 mm
    铝屏蔽
    0.7~2 MeV(对总剂量贡献26%)15~60 MeV(对总剂量贡献74%)总剂量
    增加78%
    1.68×1061.19×1061.00×1081.68×1042.43×1041.63×104
    3 mm
    铝屏蔽
    0.15~4 MeV(对总剂量贡献14%)25~200 MeV(对总剂量贡献86%)总剂量
    增加99%
    1.16×1083.53×1077.00×1081.04×1039.88×1031.27×104
    10 mm
    铝屏蔽
    0.1~1 MeV(对总剂量贡献42%)50~200 MeV(对总剂量贡献58%)总剂量
    增加118%
    1.68×1086.17×1071.97×1092.33×1024.16×1033.87×103
    下载: 导出CSV

    表  5  槽区动态变化对10 000 km×0°轨道总剂量影响

    Table  5.   Total dose variations on 10 000 km×0° orbit due to slot region dynamics

    屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
    AE8MAXAE9
    99%置信度值
    CRRES扰动值/
    SLOT模型
    99%置信度值
    AP8MAXAP9
    99%置信度值
    CRRES
    扰动值
    0.04 mm
    铝屏蔽
    0.1~1 MeV(对总剂量贡献23%)2~10 MeV(对总剂量贡献77%)总剂量
    增加20%
    9.14×1078.61×1071.06×1095.59×1061.21×1072.91×106
    1 mm
    铝屏蔽
    0.7~4 MeV(对总剂量贡献77%)15~60 MeV(对总剂量贡献23%)总剂量
    增加67%
    1.74×1061.74×1063.31×1073.77×1035.41×1031.73×103
    3 mm
    铝屏蔽
    0.1~4 MeV(对总剂量贡献82%)25~60 MeV(对总剂量贡献18%)总剂量
    增加41%
    9.17×1078.67×1071.07×1091.16×1032.00×1031.24×103
    10 mm
    铝屏蔽
    0.1~4 MeV(对总剂量贡献83%)50~100 MeV(对总剂量贡献17%)总剂量
    增加82%
    9.17×1078.67×1071.07×1092.08×1015.31×1021.03×102
    下载: 导出CSV

    表  6  槽区动态变化对12 000 km×0°轨道总剂量影响

    Table  6.   Total dose variations on 12 000 km×0° orbit due to slot region dynamics

    屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
    AE8MAXAE9
    99%置信度值
    CRRES扰动值/
    SLOT模型
    99%置信度值
    AP8MAXAP9
    99%置信度值
    CRRES
    扰动值
    0.04 mm
    铝屏蔽
    0.1~1 MeV(对总剂量贡献7%)2~8 MeV(对总剂量贡献93%)总剂量
    增加22%
    2.11×1071.72×1088.27×1083.91×1067.74×1061.47×106
    1 mm
    铝屏蔽
    0.7~4 MeV(对总剂量贡献98%)15~40 MeV(对总剂量贡献2%)总剂量
    增加157%
    1.54×1063.67×1065.24×1073.96×1021.00×1031.04×102
    3 mm
    铝屏蔽
    1.7~4 MeV(对总剂量贡献100%)可忽略总剂量
    增加63%
    2.94×1053.85×1054.29×106
    10 mm
    铝屏蔽
    0.06~6 MeV(对总剂量贡献100%)可忽略总剂量
    增加27%
    3.61×1072.69×108无数据
    下载: 导出CSV

    表  7  槽区动态变化对14 000 km×0°轨道总剂量影响

    Table  7.   Total dose variations on 14 000 km×0° orbit due to slot region dynamics

    屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
    AE8MAXAE9
    99%置信度值
    CRRES扰动值/
    SLOT模型
    99%置信度值
    AP8MAXAP9
    99%置信度值
    CRRES
    扰动值
    0.04 mm
    铝屏蔽
    0.1~1 MeV(对总剂量贡献12%)2~8 MeV(对总剂量贡献88%)总剂量
    增加21%
    2.57×1073.02×1085.59×1081.57×1063.19×1065.24×105
    1 mm
    铝屏蔽
    0.7~5 MeV(对总剂量贡献100%)可忽略总剂量
    增加91%
    2.78×1068.46×1065.59×107
    3 mm
    铝屏蔽
    1.7~4 MeV(对总剂量贡献100%)可忽略总剂量
    增加35%
    5.78×1058.05×1055.23×106
    10 mm
    铝屏蔽
    0.06~6 MeV(对总剂量贡献100%)可忽略总剂量
    增加49%
    3.97×1074.67×108无数据
    下载: 导出CSV

    表  8  常规厚度玻璃盖片屏蔽下槽区辐射环境动态变化影响

    Table  8.   Variations of radiations due to slot region dynamics for solar cell glass cover of normal thickness

    0°倾角轨道
    高度/km
    2~100 MeV捕获质子日累积通量/cm-2扰动影响
    AP8MAXAP9 99%置信度值CRRES扰动值
    60003.72×1066.66×1066.37×106等效1 MeV电子通量增加8%
    80005.65×1061.06×1076.88×106等效1 MeV电子通量增加9%
    10 0005.61×1061.21×1072.91×106等效1 MeV电子通量增加12%
    12 0003.92×1067.76×1061.47×106等效1 MeV电子通量增加10%
    14 0001.57×1063.19×1065.24×105等效1 MeV电子通量增加10%
    下载: 导出CSV

    表  9  5倍常规厚度玻璃盖片屏蔽下槽区辐射环境动态变化影响

    Table  9.   Variations of radiations due to slot region dynamics for solar cell glass cover of 5× thickness

    0°倾角轨道
    高度/km
    捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
    AE8MAXAE9
    99%置信度值
    CRRES扰动值/
    SLOT模型
    99%置信度值
    AP8MAXAP9
    99%置信度值
    CRRES
    扰动值
    6000可忽略8~100 MeV等效1 MeV电子
    通量增加0.5%
    3.70×1052.70×1053.90×105
    8000可忽略8~200 MeV等效1 MeV电子
    通量增加5%
    1.57×1052.26×1051.84×105
    10 000可忽略8~250 MeV等效1 MeV电子
    通量增加10%
    4.24×1048.49×1042.00×104
    12 0000.45~5 MeV(对太阳电池损伤贡献12%)8~250 MeV(对太阳电池损伤贡献88%)等效1 MeV电子
    通量增加23%
    3.34×1061.40×1077.85×1078.17×1031.73×1042.28×103
    14 0000.45~5 MeV(对太阳电池损伤贡献12%)8~250 MeV(对太阳电池损伤贡献88%)等效1 MeV电子
    通量增加69%
    4.86×1062.78×1078.09×1078.67×1022.23×1032.28×102
    下载: 导出CSV
  • [1] 陈文磊, 谢伦. >0.3 MeV高能电子注入辐射带槽区事件的分析与预报[J]. 地球物理学报, 2010, 53(12): 2796-2804 doi:  10.3969/j.issn.0001-5733.2010.12.002

    CHEN W L, XIE L. The analysis and forecast of the penetration of >0.3 MeV electrons to the slot region in radiation belt[J]. Chinese Journal of Geophysics, 2010, 53(12): 2796-2804 doi:  10.3969/j.issn.0001-5733.2010.12.002
    [2] BRAUTIGAM D H. CRRES in review: space weather and its effects on technology[J]. Journal of Atmospheric and Solar-Terrestrial Physics, 2002, 64(16): 1709-1721 doi:  10.1016/S1364-6826(02)00121-9
    [3] MILLAN R M, BAKER D N. Acceleration of particles to high energies in Earth’s radiation belts[J]. Space Science Reviews, 2012, 173(1/2/3/4): 103-131
    [4] THORNE R M, SHPRITS Y Y, MEREDITH N P, et al. Refilling of slot region between the inner and outer electron radiation belts during geomagnetic storms[J]. Journal of Geophysical Research: Space Physics, 2007, 112(A6): A06203
    [5] EVANS D, GARRETT H, JUN I, et al. Long-term observations of the trapped high-energy proton population (L < 4) by the NOAA polar orbiting environmental satellites (POES)[J]. Advance in Space Research, 2008, 41(8): 1261-1268 doi:  10.1016/j.asr.2007.11.028
    [6] VETTE J I. The AE-8 trapped electron model environment: WDC-A-RS-91-24[R]. NASA/Goddard Space Flight Center, 1991: 1-138
    [7] SAWYER D M, VETTE J I. AP-8 trapped proton environment for solar maximum and solar minimum: WDC-A-RS, -76-06[R]. NASA/Goddard Space Flight Center, 1971: 1-177
    [8] BRAUTIGAM D H, BELL J T. CRRESELE documentation: PL-TR-95-2128[R]. Phillips Laboratory: Environmental Research Papers 1178, 1995: 1-54
    [9] MEFFERT J D, GUSSENHOVEN M S. CRRESPRO documentation: PL-TR-94-2218[R]. Phillips Laboratory: Environmental Research Papers 1158, 1994: 1-32
    [10] SICARD-PIET A, BOSCHER D, LAZARO D, et al. A new ONERA-CNES slot electron model[J]. IEEE Transactions on Nuclear Science, 2014, 61(4): 1648-1655 doi:  10.1109/TNS.2013.2293346
    [11] GINET G P, O'BRIEN T P, HUSTON S L, et al. AE9, AP9 and SPM: new models for specifying the trapped energetic particle and space plasma environment[J]. Space Science Review, 2013, 179(1): 579-615
  • [1] 吴康, 方美华, 张昭, 周宏涛, 陈建飞, 梁筝, 王彪, 宋定一.  空间辐射对空间站航天员的剂量影响分析 . 航天器环境工程, 2023, 40(5): 475-483. doi: 10.12126/see.2023041
    [2] 刘佳强, 张振龙, 赵班池, 刘洋, 王慧元, 高辉.  地球内辐射带核心区环境特征分析及质子屏蔽的蒙特卡罗模拟 . 航天器环境工程, 2022, 39(4): 355-360. doi: 10.12126/see.2022.04.004
    [3] 施辉, 张海良, 宋思德, 曹利超, 王印权, 刘国柱, 顾祥, 吴建伟, 洪根深, 李明华, 贺琪.  MRAM的辐射效应分析及加固方法简述 . 航天器环境工程, 2021, 38(1): 106-114. doi: 10.12126/see.2021.01.017
    [4] 蔡震波, 曲少杰.  火星探测器全任务期空间环境特征与防护要点 . 航天器环境工程, 2019, 36(6): 542-548. doi: 10.12126/see.2019.06.003
    [5] 高欣, 王俊, 冯展祖, 汪伟, 杨生胜, 尹飞, 黄新宁, 薛玉雄, 把得东.  空间高速InGaAs-PIN光电二极管辐射损伤效应试验研究 . 航天器环境工程, 2019, 36(2): 139-145. doi: 10.12126/see.2019.02.006
    [6] 余永涛, 陈毓彬, 水春生, 王小强, 冯发明, 费武雄.  大容量抗辐射加固SRAM器件单粒子效应试验研究 . 航天器环境工程, 2018, 35(5): 462-467. doi: 10.12126/see.2018.05.010
    [7] 王建昭, 张庆祥, 田岱, 朱安文, 邱家稳.  木星系粒子辐射环境效应及防护关键技术 . 航天器环境工程, 2018, 35(5): 500-510. doi: 10.12126/see.2018.05.017
    [8] 李铮, 于庆奎, 罗磊, 孙毅, 梅博, 唐民.  宇航用双极器件和光电耦合器位移损伤试验研究 . 航天器环境工程, 2017, 34(1): 86-90. doi: 10.12126/see.2017.01.014
    [9] 罗磊, 张洪伟, 董艺, 梅博, 于庆奎.  一款新研制ASIC器件的单粒子效应检测与故障定位 . 航天器环境工程, 2017, 34(2): 202-206. doi: 10.12126/see.2017.02.016
    [10] 杨生胜1, 高欣2, 冯展祖2, 王俊1, 把得东1, 曹洲1空间激光通信系统光电器件辐射效应研究 . 航天器环境工程, 2017, 34(6): 571-576. doi: 10.12126/see.2017.06.001
    [11] 罗磊, 张磊, 祝名, 于庆奎, 唐民.  一种串并转换器件单粒子效应试验研究 . 航天器环境工程, 2013, 30(5): 509-511.
    [12] 杨悦, 王西京, 赵静, 吕铁鑫.  低轨航天器运行环境质子通量预测模型对比研究 . 航天器环境工程, 2013, 30(5): 481-486.
    [13] 高欣1,3, 杨生胜1,3, 冯展祖2,3, 张雷2,3空间CCD图像传感器辐射损伤评估方法 . 航天器环境工程, 2013, 30(6): 596-601. doi: 10.12126/see.2013.06.006
    [14] 周 飞, 李 强, 信太林, 韦锡峰, 张 华.  空间辐射环境引起在轨卫星故障分析与加固对策 . 航天器环境工程, 2012, 29(4): 392-396.
    [15] 姜利祥1,2, 刘宇明2, 刘国青1,2, 李 涛1,2, 丁义刚1,2, 郑慧奇2.  热控涂层GEO带电粒子辐射效应加速试验方法 . 航天器环境工程, 2012, 29(5): 527-531.
    [16] 沈自才.  深空辐射环境及其效应的分析与比较 . 航天器环境工程, 2010, 27(3): 313-320.
    [17] 于庆奎, 张大宇, 张海明, 等.  SRAM FPGA电离辐射效应试验研究 . 航天器环境工程, 2009, 26(3): 229-231.
    [18] 张庆祥, 及 莉, 王 立.  卫星空间辐射效应及防护技术飞行试验初步设想 . 航天器环境工程, 2008, 25(3): 247-250.
    [19] 刘宇明.  空间紫外辐射环境及效应研究 . 航天器环境工程, 2007, 24(6): 359-365.
    [20] 丁义刚.  空间辐射环境单粒子效应研究 . 航天器环境工程, 2007, 24(5): 283-290.
  • 加载中
图(5) / 表ll (9)
计量
  • PDF下载量:  34
  • 文章访问数:  312
  • HTML全文浏览量:  15
文章相关
  • 中图分类号:  P353; V524.3
  • 收稿日期:  2021-03-09
  • 修回日期:  2021-06-06
  • 网络出版日期:  2021-07-08
  • 刊出日期:  2021-06-25

地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响

doi: 10.12126/see.2021.03.014
    作者简介:

    呼延奇(1982—),男,博士学位,高级工程师,研究方向为航天器空间环境工程设计。E-mail: zzuhuyan@163.com

  • 中图分类号: P353; V524.3

摘要: 中轨卫星运行于地球辐射带槽区,而槽区粒子辐射环境可能存在显著涨落,增加卫星抗辐射设计输入的不确定性。文章利用典型地球辐射带模型,对中轨卫星累积性辐射效应的主要来源进行深入分析,再结合槽区粒子辐射环境动态变化特征,初步量化分析其对中轨卫星遭遇辐射效应的影响。结果表明:槽区粒子辐射环境的动态变化对星表材料及太阳电池辐射损伤的附加影响较小;槽区质子填充事件对8000 km以上高度轨道的电离总剂量有明显影响(但此类事件遭遇概率很低);槽区电子填充事件使10 000 km以上高度轨道的电离总剂量明显增大,这点必须在相应的卫星抗辐射设计要求中予以考虑。

English Abstract

呼延奇, 蔡震波. 地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响[J]. 航天器环境工程, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
引用本文: 呼延奇, 蔡震波. 地球辐射带槽区粒子环境动态变化对中轨卫星辐射效应的影响[J]. 航天器环境工程, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
HU Y Q, CAI Z B. Influence of slot region dynamics of Earth radiation belt on the radiation effect of medium-orbit satellites[J]. Spacecraft Environment Engineering, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
Citation: HU Y Q, CAI Z B. Influence of slot region dynamics of Earth radiation belt on the radiation effect of medium-orbit satellites[J]. Spacecraft Environment Engineering, 2021, 38(3): 333-343 doi:  10.12126/see.2021.03.014
    • 近年来,轨道高度5000~15 000 km的中地球轨道区域卫星的应用越来越受到关注。中地球轨道一般指高度在2000 km以上、36 000 km(GEO)以下的轨道。这类轨道的特点是大气阻力较小、地面覆盖范围广、可视时间长,具备很多其他轨道所不具备的优越条件,是用于全球导航和移动通信的主要轨道高度。

      中地球轨道卫星位于地球外辐射带和内辐射带之间的槽区。因为等离子层顶产生的哨声波在槽区引起很强的粒子投掷角散射,使得粒子不容易被捕获,形成了一个低辐射通量的区域,所以长期以来人们一般认为槽区是近地恶劣辐射环境中的一个安全区域,是适合航天器运行的“安全岛”[1]。然而,CRRES卫星在轨监测的数据显示,在强烈的太阳风暴期间,地球辐射带粒子环境存在显著的强度变化及空间位置重新分布,尤其是在槽区,甚至出现新的质子带和电子带,这深刻改变了人们的原有认识[2-5]

      如果槽区空间粒子环境的动态变化是短时存在的,则从卫星抗辐射设计角度不需太多关注;而如果这种动态变化引起的辐射强度变化达到或超过1个数量级,且持续时间长达数月(或者虽持续时间较短但动态变化频繁),则在卫星的抗辐射设计中必须采取相应的措施,以确保卫星能够耐受在轨可能遭遇的最恶劣辐射环境。

      本文利用国际上典型的地球辐射带模型,在详细分析槽区粒子辐射环境动态变化特征的基础上,系统研究这种动态变化对卫星遭受的累积性辐射效应的影响,从而有助于为中轨卫星的抗辐射设计确定合理的环境条件及约束。

    • 地球辐射带由被地球磁场捕获的带电粒子(主要为质子和电子)组成,其强度明显集中在内辐射带和外辐射带。内辐射带是靠地球最近的捕获带电粒子区域,其空间范围大致为L=1.2~2.0(L为磁壳参数,是地球空间某点的磁力线与赤道面交点处的地心距离与地球平均半径RE之比,RE=6371 km);外辐射带是离地球较远的捕获带电粒子区域,其空间范围大致为L=3.0~8.0;内、外辐射带之间存在一个粒子辐射通量较低的槽区,其空间范围为L=2.0~3.0,中心位于L≈2.5处。

      图1为当前航天工程设计中采用的标准地球辐射带模型AE8[6]和AP8[7]给出的捕获电子和捕获质子通量在磁赤道面上随磁壳参数L的变化曲线。AE8和AP8模型是长期平均静态模型,反映的是辐射带粒子通量的长期平均结果。从图1可以看出:辐射带槽区主要表现在能量>500 keV的电子通量分布上,在2.0<L<3.0的槽区内,电子通量相比内、外辐射带电子通量存在一定的下降,而且随着电子能量的升高,其通量下降的幅度越来越大;高能捕获质子主要分布在内辐射带,能量>50 MeV的质子强度中心位于L=1.5处(对应的高度约为3000 km),能量>50 MeV的质子在槽区内的通量较内辐射带中心区域存在显著的下降,然而能量<10 MeV的低能质子在槽区内的通量却处于较高水平。

      图  1  地球辐射带捕获粒子通量随L值的变化曲线

      Figure 1.  Distribution of fluxes of trapped particles in Earth’s radiation belt vs L value

      初步分析可知,高度在5000~15 000 km的中轨道所在的辐射带槽区中的捕获粒子通量并非在全能段上均处于较低的水平:能量<500 keV的电子和能量<10 MeV的质子的通量水平并没有明显下降,尤其是能量<10 MeV的质子通量实际上非常高;而对于能量>500 keV的高能电子和能量>50 MeV的高能质子,其在槽区内的通量相比内辐射带确实存在明显的下降。因此,人们普遍认为的槽区粒子辐射水平低主要体现为高能质子和高能电子的通量相对较低。

    • 槽区辐射环境动态变化与太阳活动周期、具体每次太阳风暴的状态以及空间位置等均有明显相关性,因此要获得工程可用的分析结果,需借助于通过对大量卫星实测数据进行同化及处理后建立的地球辐射带粒子环境动态模型。

      AE8MAX(AE8MIN)和AP8MAX(AP8MIN)模型分别给出太阳活动高年(低年)不同空间位置处的电子和质子能谱,然而无法描述槽区的辐射环境动态变化;其可以作为一个参照基准,通过其他动态模型的分析结果来反映槽区辐射环境的动态变化情况。

      目前,可以用于分析辐射带动态变化特性的地球辐射带粒子环境模型包括CRRES模型、SLOT模型以及AE9/AP9模型。

      1)CRRES模型

      基于CRRES卫星在太阳活动高年期间14个月的探测数据,美国空军研究实验室(AFRL)发展了能部分反映地磁扰动影响的准动态地球辐射带模型,包括质子带模型CRRESPRO[8]和电子带模型CRRESELE[9]。其中,CRRESPRO模型可以给出能量范围1~100 MeV、空间范围L=1.15~5.5的捕获质子通量,并分为CRRESAVG和CRRESACTIVE,分别描述辐射带的平静和扰动状态;而CRRESELE模型则可以给出在6种不同的地磁扰动水平下,能量范围0.7~5 MeV、空间范围L=2.5~6.5(赤道面高度9500~35 000 km)的捕获电子通量。CRRES模型还给出了CRRES卫星任务期内遭遇的平均和最恶劣捕获电子通量。

      2)SLOT槽区电子环境模型[10]

      法国国家空间研究中心ONERA基于NOAA(美国国家海洋和大气管理局)的POES卫星,以及HEO-1、HEO-3、ICO和CRRES卫星的电子通量监测数据,专门建立了一个槽区电子环境模型——ONERA-CNES,可以实现对不同置信度条件下,空间范围2.0<L<4.0、能量范围0.1~3 MeV的捕获电子通量的分析。

      3)AE9/AP9模型[11]

      NASA对经典辐射带模型AE8/AP8进行更新,推出最新一代地球辐射带模型——AE9/AP9/SPM。相比于以往的地球辐射带模型给出的均是卫星长期探测数据的平均值,对于粒子通量的涨落范围以及随时间的变化不能给予描述,使得工程应用中对粒子环境不确定性的分析缺少必要数据支持的不足,AE9/AP9/SPM模型首次给出了对粒子通量不确定性的定量描述。

      AE9/AP9/SPM模型给出3种不同粒子通量计算选项,分别是平均值或指定百分位数的量值、静态扰动条件下的量值以及考虑地球辐射带动态变化的蒙特卡罗模拟计算。利用AE9/AP9模型不同置信度条件下的槽区捕获粒子通量分析结果,可以对槽区辐射环境的动态变化情况进行描述。

    • 为了分析槽区辐射环境动态变化对卫星抗辐射设计的附加影响,须首先对槽区高度上运行的卫星遭受的辐射效应的来源进行分析。对于运行于3000~15 000 km高度范围的卫星,其辐射效应主要来源于捕获电子和捕获质子。本文首先采用地球辐射带模型AE8/AP8对不同高度轨道卫星遭遇的辐射效应的来源开展分析识别。

      进行中轨卫星抗辐射设计时,对于单粒子效应、充放电效应等瞬时辐射效应通常考虑最恶劣的环境条件,故槽区辐射环境的动态变化对此类效应的防护设计不产生影响;而对于电离总剂量效应、位移损伤效应等累积效应通常考虑长期平均的辐射环境条件,则可能受到槽区辐射环境动态变化的影响。因此,本文重点关注累积辐射效应。

    • 中轨卫星上元器件及材料遭受的电离总剂量损伤主要来自地球辐射带捕获电子和捕获质子。为了分析中轨卫星上应用的元器件或材料遭受的电离总剂量,选择等效铝屏蔽厚度为0.04、1、3、10 mm的4种屏蔽状态,分别模拟星表热控材料、卫星舱板内部电缆以及舱内设备内部等不同应用条件下的屏蔽状态。空间中的带电粒子是连续能谱的,对于不同屏蔽厚度下的元器件或材料,其电离总剂量的主要贡献来自于不同能段的电子和质子:低能粒子射程短,无法穿透卫星舱板屏蔽,只能在星表材料中沉积电离剂量;对于屏蔽厚度较大的元器件或材料,其电离总剂量主要来自于高能粒子产生的次级辐射(如轫致辐射)。

      本文利用AE8/AP8模型分析不同高度、0°倾角的中地球轨道上的日累积捕获电子和捕获质子通量,然后利用Shieldose-2软件分析了不同铝屏蔽厚度下各轨道上不同能量的捕获电子和质子对总剂量的贡献。图2图3分别给出了1 mm和3 mm铝屏蔽状态下,不同能量的捕获电子和质子对总剂量贡献的累积分数。将特定轨道上的总剂量归一化到1,然后将能量小于某值的电子及质子产生的总剂量除以轨道累积总剂量,即获得相应能量粒子对总剂量的贡献分数;卫星轨道上的总剂量是整个连续能谱的电子和质子产生的电离总剂量的叠加。

      图  2  1 mm铝屏蔽下不同能量的捕获电子和捕获质子对总剂量的贡献

      Figure 2.  Contribution of electrons (left) and protons (right) of different energies to the total dose with 1 mm aluminum shield

      图  3  3 mm铝屏蔽下不同能量的捕获电子和捕获质子对总剂量的贡献

      Figure 3.  Contribution of electrons (left) and protons (right) of different energies to the total dose with 3 mm aluminum shield

      表1对不同轨道高度上来自不同能量捕获电子和质子对总剂量的贡献进行了详细分析和归纳。表中数据显示,星表材料的总剂量主要来自于1~20 MeV的质子和100~600 keV的电子,且随着轨道高度的变化,质子和电子对总剂量的贡献权重不断变化。而经过舱板1 mm铝屏蔽作用,对总剂量有贡献的质子能量范围主要在15~100 MeV,电子能量范围在700 keV~几MeV。随着屏蔽厚度的增加,低能质子无法穿透相应的屏蔽,故对总剂量有贡献的质子能量下限不断提高,而电子在屏蔽材料中减速产生的轫致辐射使其对总剂量的贡献逐渐增加;总体来看,经过3~10 mm铝屏蔽后,对总剂量有贡献的质子能量范围一般在几十~几百MeV,电子能量范围在几百keV~几MeV。

      表 1  不同高度轨道卫星总剂量来源分析

      Table 1.  Particle source for total dose on satellite at different altitudes of orbits

      轨道高度/
      km
      不同屏蔽状态下不同能量的捕获电子和质子对总剂量的贡献/%(粒子能量/MeV)粒子类型
      0.04 mm铝屏蔽状态1 mm铝屏蔽状态3 mm铝屏蔽状态10 mm铝屏蔽状态
      3000约3%(2~20)约47%(15~100)约76%(25~300)约98%(50~400)质子
      约97%(0.1~0.6)约53%(0.7~3)约24%(1.7~4)约2%(0.3~4.5)电子
      4000约8%(2~20)约65%(15~100)约84%(25~300)约98%(50~400)质子
      约92%(0.1~0.6)约35%(0.7~3)约16%(1.7~4)约2%(0.3~4.5)电子
      5000约17%(2~20)约74%(15~100)约90%(25~300)约94%(50~400)质子
      约83%(0.1~0.6)约23%(0.7~3)约10%(0.4~4)约6%(0.3~4.5)电子
      6000约30%(2~20)约80%(15~100)约92%(25~300)约87%(50~400)质子
      约70%(0.1~0.6)约20%(0.7~2)约8%(0.15~4)约13%(0.1~1)电子
      7000约44%(2~20)约82%(15~60)约92%(25~300)约77%(50~400)质子
      约56%(0.1~0.6)约18%(0.7~2)约8%(0.15~4)约13%(0.1~1)电子
      8000约56%(2~20)约74%(15~60)约86%(25~200)约58%(50~200)质子
      约44%(0.1~0.6)约26%(0.7~2)约14%(0.15~4)约42%(0.1~2)电子
      9000约66%(2~10)约50%(15~60)约64%(25~100)约36%(50~200)质子
      约34%(0.1~0.6)约50%(0.7~3)约36%(0.1~4)约64%(0.1~2)电子
      10 000约77%(2~10)约23%(15~60)约18%(25~60)约17%(50~100)质子
      约23%(0.1~1)约77%(0.7~4)约82%(0.1~4)约83%(0.1~4)电子
      11 000约86%(2~8)约7%(15~40)约3%(25~40)约6%(50~100)质子
      约14%(0.1~1)约93%(0.7~4)约97%(0.1~4)约94%(0.06~4)电子
      12 000约93%(2~8)约2%(15~40)可忽略可忽略质子
      约7%(0.1~1)约98%(0.7~4)约100%(1.7~4)约100%(0.06~6)电子
      13 000约92%(2~8)可忽略可忽略可忽略质子
      约8%(0.1~1)约100%(0.7~5)约100%(1.7~5)约100%(0.06~6)电子
      14 000约88%(2~8)可忽略可忽略可忽略质子
      约12%(0.1~1)约100%(0.7~5)约100%(1.7~5)约100%(0.06~6)电子
      15 000约81%(2~8)可忽略可忽略可忽略质子
      约19%(0.1~1)约100%(0.7~5)约100%(1.7~5)约100%(0.06~6)电子

      综上可知:

      1)星表材料的总剂量来源与舱板和设备内部的存在显著差异——5000 km及以下高度轨道星表材料的总剂量主要来自捕获电子;9000 km以上高度轨道的总剂量主要来自于捕获质子;5000~9000 km高度范围,质子和电子对星表材料的总剂量均有贡献,其权重此消彼长。

      2)屏蔽厚度为1 mm状态下,8000 km以上高度轨道的总剂量主要来自于捕获电子;8000 km以下高度轨道的总剂量主要来自于捕获质子,但捕获电子也占有一定份额。

      3)对于屏蔽厚度3~10 mm的状态,电子和质子对总剂量的贡献基本上是各主导一段——9000 km以下高度轨道的总剂量主要来自于捕获质子;9000 km以上高度轨道的总剂量主要来自于捕获电子。

    • 空间高能粒子与材料间除了通过电离相互作用交换能量外,还可以通过非电离相互作用交换能量,产生位移效应和非电离能量损失。受位移效应影响最大的主要是利用少数载流子工作的器件,例如太阳电池等。

      20世纪80年代,美国喷气推进实验室(JPL)发展了等效注量法,通过采用不同能量的质子和电子进行辐照试验,然后利用获得的试验数据,将不同能量的空间带电粒子对太阳电池造成的位移损伤与一定通量的1 MeV电子对太阳电池的损伤进行比对,以获得不同能量和种类的粒子与1 MeV电子产生同样损伤的等效系数;最终把太阳电池在轨期间由连续能谱的空间带电粒子辐射导致的位移损伤等效为1 MeV电子损伤通量,以定量描述太阳电池遭受的位移损伤。

      目前,航天器上广泛采用的是三结GaAs太阳电池,选用的玻璃盖片厚度通常为90 μm和120 μm,按照常用的掺铈抗辐照玻璃的密度2.56 g/cm3计算,太阳电池玻璃盖片的屏蔽面密度相当于0.023~0.031 g/cm2。本文取玻璃盖片屏蔽面密度为0.023 g/cm2,分析不同高度轨道上的捕获电子和捕获质子对太阳电池辐射损伤的贡献。

      图4给出了3000~15 000 km、0°倾角轨道上不同能量的捕获电子和捕获质子产生的年均太阳电池辐射损伤的等效1 MeV电子通量。从图中可以看出,与总剂量的来源截然不同,在3000~15 000 km高度轨道上,捕获质子产生的太阳电池辐射损伤比捕获电子产生的高2个数量级,且随着轨道高度的变化这种差异并没有显著改变。需要特别说明的是,在6000~15000 km高度轨道区域及常规厚度玻璃盖片屏蔽下,捕获质子产生的平均太阳电池辐射损伤等效1 MeV电子通量可达1016~1017/ (cm2·a)。若按照5年在轨寿命考虑,太阳电池遭受的辐射损伤将非常严重。

      图  4  不同轨道上捕获电子和质子对太阳电池辐射年累积等效1 MeV电子通量

      Figure 4.  Annual cumulative radiation on solar cell equivalent to 1 MeV electron flux for electrons (left) and protons (right) on different orbits

      为了进一步明确对太阳电池产生辐射损伤的主要贡献来自哪个能段的粒子,针对不能通量质子和电子在不同轨道上对太阳电池产生的辐射损伤进行分析,这里考虑常规厚度玻璃盖片(相当于屏蔽面密度0.023 g/cm2)和5倍于常规厚度的玻璃盖片(相当于屏蔽面密度0.115 g/cm2)2种情况,具体结果见图5

      图  5  常规厚度和5倍常规厚度玻璃盖片屏蔽下太阳电池辐射损伤来源分析

      Figure 5.  Particle source of radiation damage on solar cell under glass cover of normal thickness(left) and 5× thickness

      表2详细归纳了3000~15 000 km、0°倾角轨道上太阳电池辐射损伤的来源。从中可以看出:在常规厚度玻璃盖片屏蔽下,捕获质子是太阳电池辐射损伤的主要来源;与质子相比,捕获电子产生的损伤可以忽略。采用5倍常规厚度的盖片玻璃后,对于12 000 km高度以下的卫星,其太阳电池辐射损伤仍旧主要来自8.0~几百MeV的质子;而在12 000~15 000 km高度范围,0.45~几MeV的捕获电子对太阳电池辐射损伤的贡献逐渐占据主导地位。

      表 2  不同高度轨道卫星太阳电池辐射损伤来源分析

      Table 2.  Particle source for solar cell radiation damage on different altitudes of satellite orbits

      0°倾角轨道
      高度/km
      不同屏蔽状态下不同能量的捕获粒子
      对辐射损伤的贡献/%(粒子能量/MeV)
      粒子类型
      常规玻璃盖片5倍常规厚度玻璃盖片
      3000约99%(2~100)约99%(8~100)质子
      约1%(0.2~1)约1%(0.45~2.5)电子
      6000约99.9%(2~100)约99.8%(8~100)质子
      约0.1%(0.2~1)约0.2%(0.45~2.5)电子
      9000约99.97%(2~100)约99.8%(8~250)质子
      约0.03%(0.2~1)约0.2%(0.45~2.5)电子
      12 000约99.9%(2~100)约88%(8~250)质子
      约0.1%(0.2~1)约12%(0.45~5)电子
      15 000约99%(2~6)约10%(8~20)质子
      约1%(0.2~4)约90%(0.45~5)电子
    • 星上元器件及材料周围的屏蔽状态不同,则对其产生总剂量效应的主要辐射源的能段不同。因此,本文首先分析槽区电子或质子填充事件后相应能段的粒子通量变化情况。以AP8MAX和AE8MAX模型分析的结果代表卫星当前设计的平均值,捕获质子通量扰动值采用AP9模型99%置信度值和CRRESPRO模型扰动值之间的最大值,捕获电子通量扰动值采用AE9模型99%置信度值、CRRESELE模型扰动值以及SLOT模型99%置信度值之间的最大值。同时,考虑槽区扰动状态在卫星寿命期内可能出现的概率,对于5年寿命的中轨卫星,考虑遭遇1次类似CRRES卫星经历的质子填充事件(持续时间半年,约占寿命期的10%);对于电子填充事件,按照在轨时间的5%作保守估计。最终,槽区辐射环境动态变化对卫星辐射效应的影响可按

      $$ \begin{split} {{\rm{factor}}}=&\dfrac{{\text{捕获电子通量扰动值}}}{{\text{捕获电子通量设计值}}}\times {W_{\rm{e}}}\times {P_{\rm{e}}}+\\ &\dfrac{{\text{捕获质子通量扰动值}}}{{\text{捕获质子通量设计值}}}\times {W_{\rm{p}}}\times {P_{\rm{p}}} \end{split} $$

      进行估计。其中:WeWp分别代表捕获电子效应和捕获质子效应的贡献权重;PePp分别代表遭遇电子扰动和质子扰动的时间概率。

    • 表3给出槽区辐射环境动态变化对6000 km×0°轨道卫星总剂量影响的分析过程及结果。可以看出,以AE8/AP8模型获得的捕获电子和质子通量作为基准,槽区粒子填充事件对该轨道卫星总剂量的影响主要表现在屏蔽厚度≥3 mm的情况下,这是由于此种情况下高能质子通量的增加会引起相应的总剂量增加约70%;而对于其他屏蔽状态的器件及材料,影响可以忽略不计。

      表 3  槽区动态变化对6000 km×0°轨道总剂量影响

      Table 3.  Total dose variations on 6000 km×0° orbit due to slot region dynamics

      屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
      AE8MAXAE9
      99%置信度值
      CRRES扰动值/
      SLOT模型
      99%置信度值
      AP8MAXAP9
      99%置信度值
      CRRES
      扰动值
      0.04 mm
      铝屏蔽
      0.1~0.6 MeV(对总剂量贡献70%)2~20 MeV(对总剂量贡献30%)总剂量
      增加0%
      2.26×1081.58×108无数据3.72×1066.66×1066.37×106
      1 mm
      铝屏蔽
      0.7~2 MeV(对总剂量贡献20%)15~100 MeV(对总剂量贡献80%)总剂量
      增加1%
      1.54×1062.54×106无数据4.88×1045.73×1042.30×104
      3 mm
      铝屏蔽
      0.15~4 MeV(对总剂量贡献8%)25~300 MeV(对总剂量贡献92%)总剂量
      增加23%
      1.64×1089.11×107无数据9.21×1033.32×1045.82×103
      10 mm
      铝屏蔽
      0.1~1 MeV(对总剂量贡献13%)50~400 MeV(对总剂量贡献87%)总剂量
      增加73%
      2.31×1081.62×108无数据2.40×1032.27×1041.31×103
    • 表4给出槽区辐射环境动态变化对8000 km×0°轨道卫星总剂量影响的分析过程及结果。可以看出,在该轨道上,捕获电子通量和20 MeV以上的高能质子通量均可能增加1个数量级,故对等效铝屏蔽厚度≥1 mm的器件及材料,其遭受的总剂量可能增加78%~118%,受辐射带动态变化的影响比较显著。

      表 4  槽区动态变化对8000 km×0°轨道总剂量影响

      Table 4.  Total dose variations on 8000 km×0° orbit due to slot region dynamic changes

      屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
      AE8MAXAE9
      99%置信度值
      CRRES扰动值/
      SLOT模型
      99%置信度值
      AP8MAXAP9
      99%置信度值
      CRRES
      扰动值
      0.04 mm
      铝屏蔽
      0.1~0.6 MeV(对总剂量贡献44%)2~20 MeV(对总剂量贡献56%)总剂量
      增加25%
      1.63×1086.01×1071.83×1095.65×1061.06×1076.86×106
      1 mm
      铝屏蔽
      0.7~2 MeV(对总剂量贡献26%)15~60 MeV(对总剂量贡献74%)总剂量
      增加78%
      1.68×1061.19×1061.00×1081.68×1042.43×1041.63×104
      3 mm
      铝屏蔽
      0.15~4 MeV(对总剂量贡献14%)25~200 MeV(对总剂量贡献86%)总剂量
      增加99%
      1.16×1083.53×1077.00×1081.04×1039.88×1031.27×104
      10 mm
      铝屏蔽
      0.1~1 MeV(对总剂量贡献42%)50~200 MeV(对总剂量贡献58%)总剂量
      增加118%
      1.68×1086.17×1071.97×1092.33×1024.16×1033.87×103
    • 表5给出槽区辐射环境动态变化对10 000 km×0°轨道卫星总剂量影响的分析过程及结果。可以看出,除星表材料外,该轨道上的总剂量主要来源于捕获电子,而捕获电子通量可能增大1~2个数量级,导致器件及材料遭受的总剂量可能增加20%~82%。

      表 5  槽区动态变化对10 000 km×0°轨道总剂量影响

      Table 5.  Total dose variations on 10 000 km×0° orbit due to slot region dynamics

      屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
      AE8MAXAE9
      99%置信度值
      CRRES扰动值/
      SLOT模型
      99%置信度值
      AP8MAXAP9
      99%置信度值
      CRRES
      扰动值
      0.04 mm
      铝屏蔽
      0.1~1 MeV(对总剂量贡献23%)2~10 MeV(对总剂量贡献77%)总剂量
      增加20%
      9.14×1078.61×1071.06×1095.59×1061.21×1072.91×106
      1 mm
      铝屏蔽
      0.7~4 MeV(对总剂量贡献77%)15~60 MeV(对总剂量贡献23%)总剂量
      增加67%
      1.74×1061.74×1063.31×1073.77×1035.41×1031.73×103
      3 mm
      铝屏蔽
      0.1~4 MeV(对总剂量贡献82%)25~60 MeV(对总剂量贡献18%)总剂量
      增加41%
      9.17×1078.67×1071.07×1091.16×1032.00×1031.24×103
      10 mm
      铝屏蔽
      0.1~4 MeV(对总剂量贡献83%)50~100 MeV(对总剂量贡献17%)总剂量
      增加82%
      9.17×1078.67×1071.07×1092.08×1015.31×1021.03×102
    • 表6给出槽区辐射环境动态变化对12 000 km×0°轨道卫星总剂量影响的分析过程及结果。可以看出:该轨道上星表材料的总剂量主要来自2~8 MeV的质子,此能量段的质子受槽区填充事件的影响较小;舱内元器件及材料的总剂量主要来自捕获电子,而捕获电子通量可能增大1个数量级,导致器件及材料遭受的总剂量可能增加27%~157%。

      表 6  槽区动态变化对12 000 km×0°轨道总剂量影响

      Table 6.  Total dose variations on 12 000 km×0° orbit due to slot region dynamics

      屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
      AE8MAXAE9
      99%置信度值
      CRRES扰动值/
      SLOT模型
      99%置信度值
      AP8MAXAP9
      99%置信度值
      CRRES
      扰动值
      0.04 mm
      铝屏蔽
      0.1~1 MeV(对总剂量贡献7%)2~8 MeV(对总剂量贡献93%)总剂量
      增加22%
      2.11×1071.72×1088.27×1083.91×1067.74×1061.47×106
      1 mm
      铝屏蔽
      0.7~4 MeV(对总剂量贡献98%)15~40 MeV(对总剂量贡献2%)总剂量
      增加157%
      1.54×1063.67×1065.24×1073.96×1021.00×1031.04×102
      3 mm
      铝屏蔽
      1.7~4 MeV(对总剂量贡献100%)可忽略总剂量
      增加63%
      2.94×1053.85×1054.29×106
      10 mm
      铝屏蔽
      0.06~6 MeV(对总剂量贡献100%)可忽略总剂量
      增加27%
      3.61×1072.69×108无数据
    • 表7给出槽区辐射环境动态变化对14 000 km×0°轨道卫星总剂量的影响分析过程及结果。可以看出,与12 000 km×0°轨道类似,该轨道上的槽区质子填充事件对星表材料的总剂量影响较小;而捕获电子通量可能增大1个数量级,导致该轨道上的舱内器件及材料遭受的总剂量可能增加35%~91%。

      表 7  槽区动态变化对14 000 km×0°轨道总剂量影响

      Table 7.  Total dose variations on 14 000 km×0° orbit due to slot region dynamics

      屏蔽状态捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
      AE8MAXAE9
      99%置信度值
      CRRES扰动值/
      SLOT模型
      99%置信度值
      AP8MAXAP9
      99%置信度值
      CRRES
      扰动值
      0.04 mm
      铝屏蔽
      0.1~1 MeV(对总剂量贡献12%)2~8 MeV(对总剂量贡献88%)总剂量
      增加21%
      2.57×1073.02×1085.59×1081.57×1063.19×1065.24×105
      1 mm
      铝屏蔽
      0.7~5 MeV(对总剂量贡献100%)可忽略总剂量
      增加91%
      2.78×1068.46×1065.59×107
      3 mm
      铝屏蔽
      1.7~4 MeV(对总剂量贡献100%)可忽略总剂量
      增加35%
      5.78×1058.05×1055.23×106
      10 mm
      铝屏蔽
      0.06~6 MeV(对总剂量贡献100%)可忽略总剂量
      增加49%
      3.97×1074.67×108无数据
    • 与总剂量效应分析类似,为了分析槽区辐射环境动态变化对太阳电池辐射损伤效应的影响。我们同样采用AP8MAX和AE8MAX模型获得的相应能段的粒子通量作为基准,再利用AP9模型、SLOT模型及CRRES模型的99%置信度值或扰动值来描述槽区扰动情况的粒子通量,进而分析其对太阳电池辐射损伤效应的影响。

      表8给出6000~14 000 km高度轨道上,常规厚度玻璃盖片屏蔽下,槽区质子填充事件引起的太阳电池辐射损伤等效1 MeV电子通量的变化。可以看到,由于该轨道区域卫星的太阳电池辐射损伤99%以上来自2~100 MeV的质子,故考虑遭遇类似CRRES卫星所遭受的质子填充事件且持续半年时间,其对太阳电池辐射损伤等效1 MeV电子通量的影响只增加10%左右。而太阳电池电参数的衰退与等效1 MeV电子通量的对数相关,因此,可以认为槽区辐射环境动态变化对太阳电池辐射损伤(常规厚度玻璃盖片屏蔽情况)的影响非常小。

      表9给出6000~14 000 km轨道上,采用5倍常规厚度玻璃盖片屏蔽下,槽区辐射环境动态变化引起的太阳电池辐射损伤等效1 MeV电子通量的变化。可以看到,在12 000 km以下高度,太阳电池辐射损伤主要来自8 MeV以上的质子,相应的槽区质子填充事件引起的太阳电池辐射损伤等效1 MeV电子通量增大不超过23%。同样地,考虑到太阳电池电参数的衰退与等效1 MeV电子通量的对数相关,可以认为槽区辐射环境动态变化对太阳电池辐射损伤(5倍常规厚度玻璃盖片屏蔽情况)的影响非常小。当轨道高度提升至14 000 km时,太阳电池辐射损伤主要来自0.5 MeV以上的电子,而槽区捕获电子通量可能发生数量级的涨落,因此太阳电池辐射损伤等效1 MeV电子通量可能增大近70%。但整体来看,对于6000~14 000 km高度中地球轨道,槽区辐射环境动态变化对太阳电池辐射损伤效应的影响较小。

      表 8  常规厚度玻璃盖片屏蔽下槽区辐射环境动态变化影响

      Table 8.  Variations of radiations due to slot region dynamics for solar cell glass cover of normal thickness

      0°倾角轨道
      高度/km
      2~100 MeV捕获质子日累积通量/cm-2扰动影响
      AP8MAXAP9 99%置信度值CRRES扰动值
      60003.72×1066.66×1066.37×106等效1 MeV电子通量增加8%
      80005.65×1061.06×1076.88×106等效1 MeV电子通量增加9%
      10 0005.61×1061.21×1072.91×106等效1 MeV电子通量增加12%
      12 0003.92×1067.76×1061.47×106等效1 MeV电子通量增加10%
      14 0001.57×1063.19×1065.24×105等效1 MeV电子通量增加10%

      表 9  5倍常规厚度玻璃盖片屏蔽下槽区辐射环境动态变化影响

      Table 9.  Variations of radiations due to slot region dynamics for solar cell glass cover of 5× thickness

      0°倾角轨道
      高度/km
      捕获电子日累积通量/cm-2捕获质子日累积通量/cm-2扰动影响
      AE8MAXAE9
      99%置信度值
      CRRES扰动值/
      SLOT模型
      99%置信度值
      AP8MAXAP9
      99%置信度值
      CRRES
      扰动值
      6000可忽略8~100 MeV等效1 MeV电子
      通量增加0.5%
      3.70×1052.70×1053.90×105
      8000可忽略8~200 MeV等效1 MeV电子
      通量增加5%
      1.57×1052.26×1051.84×105
      10 000可忽略8~250 MeV等效1 MeV电子
      通量增加10%
      4.24×1048.49×1042.00×104
      12 0000.45~5 MeV(对太阳电池损伤贡献12%)8~250 MeV(对太阳电池损伤贡献88%)等效1 MeV电子
      通量增加23%
      3.34×1061.40×1077.85×1078.17×1031.73×1042.28×103
      14 0000.45~5 MeV(对太阳电池损伤贡献12%)8~250 MeV(对太阳电池损伤贡献88%)等效1 MeV电子
      通量增加69%
      4.86×1062.78×1078.09×1078.67×1022.23×1032.28×102
    • 本文在分析槽区扰动事件引起的捕获粒子环境条件变化的基础上,进一步分析了不同屏蔽状态下的中轨卫星遭遇的总剂量效应和太阳电池辐射损伤效应的主要贡献粒子类别和能段,从而对槽区扰动变化对卫星辐射效应的影响开展了初步评估;且分析时所采用的扰动状态下的环境条件均是最恶劣情况,故分析结果应该能够覆盖卫星在轨实际遭遇的情况。结果表明:

      1)对于星表或屏蔽厚度较小的位置,如星表热控材料或太阳电池片,槽区扰动变化对其整个寿命期内相应的辐射效应影响较小。

      2)6000~8000 km高度轨道,卫星舱内元器件及材料的总剂量主要来源于捕获质子,故如果卫星遭遇强烈的槽区质子填充事件(此类事件发生的条件比较苛刻,遭遇的概率较低),可能引起相应的总剂量变化最大增加约1倍。相对而言,6000 km高度轨道受槽区扰动事件影响较小,而8000 km高度轨道则会受到明显影响,但此类事件在卫星寿命期内不一定出现。

      3)10 000 km以上高度轨道,卫星舱内元器件及材料的总剂量主要来源于捕获电子,而槽区电子填充事件在卫星寿命期内发生较频繁,但持续时间一般较短,且不可能一直维持较高的扰动水平。按照卫星1年遭遇18天(占比5%)槽区电子填充事件的概率估计,此类事件可能使10 000~14 000 km高度轨道卫星所遭受的总剂量最大增加约1.5倍。

参考文献 (11)

目录

    /

    返回文章
    返回