李宗耀, 張科科,2, 胡海鷹, 盛 蕾

（1.中國科學(xué)院 上海微小衛(wèi)星工程中心，上海 201203；2.中國科學(xué)院大學(xué)，北京 100039）

SAA域選劃對(duì)天基探測(cè)衛(wèi)星影響簡析

李宗耀1, 張科科1,2, 胡海鷹1, 盛 蕾1

（1.中國科學(xué)院 上海微小衛(wèi)星工程中心，上海 201203；2.中國科學(xué)院大學(xué)，北京 100039）

通過仿真計(jì)算南大西洋異常區(qū)（SAA）高能粒子的通量與能量，分析不同通量閾值與能量閾值對(duì)SAA空間分布的影響，并以科學(xué)實(shí)驗(yàn)衛(wèi)星為研究對(duì)象，對(duì)其歷經(jīng)SAA的時(shí)間百分比進(jìn)行連續(xù)24 h和6個(gè)月的統(tǒng)計(jì)，獲得了SAA對(duì)衛(wèi)星常規(guī)運(yùn)行的工作時(shí)間影響。結(jié)果表明：在50 MeV能量閾值條件下，500 cm-2?s-1較100 cm-2?s-1通量閾值對(duì)SAA空間分布的影響范圍在緯度上相差約23.35%，經(jīng)度上相差約27.06%，歷經(jīng)SAA時(shí)間上縮短約48.0%；在100 cm-2?s-1通量閾值條件下，150 MeV較50 MeV能量閾值對(duì)SAA空間分布的影響范圍在緯度上相差約13.86%，經(jīng)度上相差約17.29%，歷經(jīng)SAA時(shí)間上縮短約31.7%。

南大西洋異常區(qū)；天基觀測(cè)；高能質(zhì)子；時(shí)間百分比

0 引言

南大西洋異常區(qū)（South Atlantic Anomaly, SAA）的中心位置位于西經(jīng)45°、南緯30°，由于地磁弱，其上空出現(xiàn)了范?艾倫輻射帶凹陷，從而導(dǎo)致空間高能粒子在該區(qū)域聚集，是航天器異?；蚬收系母甙l(fā)區(qū)[1-2]。低軌衛(wèi)星特別是攜帶高靈敏度有效載荷的衛(wèi)星飛經(jīng)SAA時(shí)，高能粒子將以表面充電、單粒子翻轉(zhuǎn)等形式對(duì)星上的微電子元件、功能材料、光學(xué)器件等造成輻射損傷[1-5]，甚至導(dǎo)致儀器的損壞。如“哈勃”太空望遠(yuǎn)鏡（HST）、美國“中段空間試驗(yàn)”（MSX）衛(wèi)星和加拿大“藍(lán)寶石”（Sapphire）微小衛(wèi)星等在歷經(jīng)SAA時(shí)盡管采取了關(guān)閉鏡蓋、載荷斷電、軟件算法修正、CCD退火處理等手段，但還是遭受到SAA高能粒子不同程度的輻射損傷[1,6-8]。美國NASA和歐洲ESA等機(jī)構(gòu)針對(duì)空間高能輻射環(huán)境開展了多年的研究，并初步建立了空間輻射環(huán)境模型，開展了空間高能粒子分布計(jì)算的在線服務(wù)。例如 ESA 的SPENVIS[5,9]軟件可以計(jì)算空間高能粒子分布。盡管如此，在衛(wèi)星方案制定和研制初期，如何開展SAA中的粒子輻射對(duì)衛(wèi)星的影響評(píng)估，制定合理有效的規(guī)避策略，事關(guān)衛(wèi)星的在軌工作效率。

我國2013年7月發(fā)射的科學(xué)實(shí)驗(yàn)衛(wèi)星，運(yùn)行于670 km高度的晨昏太陽同步軌道，傾角97°。星上裝載的高靈敏度相機(jī)，在歷經(jīng)SAA時(shí)受到了大量高能粒子輻射的影響，拍攝的圖像中出現(xiàn)大量白色亮點(diǎn)、亮斑和亮線等，導(dǎo)致探測(cè)圖像背景噪聲異常增加，與恒星背景混雜在一起，不僅影響空間碎片目標(biāo)的識(shí)別和定位，還會(huì)降低探測(cè)相機(jī)的性能、壽命和任務(wù)的探測(cè)效率[1]。

本文基于SAA的高能粒子分布規(guī)律和試驗(yàn)衛(wèi)星約束條件，通過仿真分析SAA的邊界選劃方法，給出衛(wèi)星歷經(jīng)SAA時(shí)受影響的程度（用時(shí)間百分比表示），為提高衛(wèi)星在軌工作效率及輔助地面預(yù)估研判飛越SAA的影響提供數(shù)據(jù)支撐。

1 SAA高能粒子影響空間大小分析

SAA中的高能粒子主要包括高能質(zhì)子、高能電子和少量的高能量宇宙射線粒子等[10-13]，其中高能質(zhì)子的能量范圍在0.1～400 MeV甚至更大[3]，在空間分布上則主要集中在南大西洋上空；高能電子的能量在1 MeV以下，遍布整個(gè)南北極區(qū)與大西洋上空。在進(jìn)行低地球軌道衛(wèi)星的空間輻射評(píng)估時(shí)，通常以能承受的輻射總劑量（主要受空間高能粒子通量大小的影響）來衡量考慮，而忽略空間高能粒子的短時(shí)轟擊效應(yīng)（主要受空間高能粒子能量大小的影響）。本文使用SPENVIS[5]計(jì)算工具，調(diào)用NASA的輻射帶電子通量模型AE-9[2-3]和質(zhì)子通量模型AP-9[2,14]，對(duì)空間高能粒子輻射環(huán)境進(jìn)行仿真計(jì)算，以得到高能粒子的分布，并根據(jù)試驗(yàn)衛(wèi)星自身約束條件，開展SAA高能粒子輻射通量和能量分布估算。

1.1 SAA高能質(zhì)子輻射通量估算

衛(wèi)星的 CCD類有效載荷在設(shè)計(jì)選用元器件時(shí)，須考慮其輻射總劑量指標(biāo)，它是衡量在軌工作能力的重要因素之一?！肮碧胀h(yuǎn)鏡可利用在軌圖像中有效強(qiáng)宇宙線事件反推評(píng)估高能粒子通量大小[7]，并根據(jù)該方法，估算了試驗(yàn)衛(wèi)星在1600 ms積分時(shí)間下遭受＞50 MeV的高能粒子通量約為500 cm-2? s-1。

SAA的空間粒子分布（以高能粒子輻射通量和能量衡量）在不同的高度、地理位置各不相同。

本文定義了試驗(yàn)衛(wèi)星相機(jī)的高能粒子通量閾值——經(jīng)機(jī)殼防護(hù)后相機(jī)能夠承受高能粒子轟擊而其性能又不會(huì)有明顯降低（單幅圖像的熱像元數(shù)目增幅小于20%）的高能粒子通量（以下簡稱通量閾值）作為評(píng)價(jià)的主要因素之一。考慮試驗(yàn)衛(wèi)星相機(jī)的防護(hù)措施，以50 MeV高能質(zhì)子能量作為輸入條件，仿真計(jì)算得到 100、200、300、500 cm-2?s-1的通量閾值下所對(duì)應(yīng)的SAA空間分布邊界，如圖1所示；各閾值對(duì)應(yīng)的地理分布范圍如表1所示?？梢钥吹剑篠AA的空間分布邊界隨著通量閾值的增加而遞減，若設(shè)定在50 MeV、100 cm-2?s-1條件下計(jì)算得到的 SAA分布邊界大小作為基礎(chǔ)邊界，在50 MeV、500 cm-2?s-1條件下的SAA分布邊界的緯度相對(duì)縮小了23.35%，經(jīng)度相對(duì)縮小了27.06%。

圖1 能量為50 MeV以及不同通量閾值條件下SAA空間分布Fig.1 The SAA spatial distribution for different flux thresholds at 50 MeV

表1 質(zhì)子能量為50 MeV以及不同通量閾值條件下的SAA空間分布邊界Table 1 The SAA spatial distribution boundary for different flux thresholds at 50 MeV

1.2 SAA高能質(zhì)子能量估算

與通量閾值的定義類似，本文也定義了 CCD相機(jī)的質(zhì)子能量閾值。應(yīng)用SPENVIS計(jì)算不同能量閾值對(duì)應(yīng)的SAA粒子空間分布。設(shè)定高能質(zhì)子輻射通量閾值為100 cm-2?s-1時(shí)，選擇50、70、100、150 MeV分別作為高能質(zhì)子能量閾值進(jìn)行研究。仿真計(jì)算得到的SAA質(zhì)子空間分布如圖2所示，各閾值對(duì)應(yīng)的地理分布范圍如表2所示。對(duì)比表中數(shù)據(jù)可以發(fā)現(xiàn)：隨著能量閾值的增加，SAA空間分布邊界向中心區(qū)域縮小，且各個(gè)方向的縮小程度不同，南部和東部縮小相比北部和西部較為明顯，在緯度上的縮小程度較經(jīng)度上來得大。若設(shè)定在50 MeV、100 cm-2?s-1條件下計(jì)算得到的SAA分布邊界作為基礎(chǔ)邊界，在150 MeV、100 cm-2?s-1條件下的SAA分布邊界在緯度上相對(duì)縮小了13.86%，經(jīng)度上相對(duì)縮小了17.29%。

圖2 通量為100 cm-2?s-1以及不同能量閾值條件下SAA空間分布Fig.2 The SAA spatial distribution for different energy thresholds at 100 cm-2?s-1

表2 通量為100 cm-2?s-1，不同能量閾值對(duì)應(yīng)的SAA空間分布邊界Table 2 The SAA spatial distribution boundary for different energy thresholds at 100 cm-2?s-1

2 歷經(jīng)SAA時(shí)間統(tǒng)計(jì)

如上文分析，SAA的空間分布與高能粒子通量和能量的大小相關(guān)。低地球軌道航天器各結(jié)構(gòu)表面抗高能粒子輻射的能力不同，將導(dǎo)致其受SAA空間分布邊界的影響也不同。特別是對(duì)高靈敏度的有效載荷，不同的抗輻射能力，決定了其受SAA影響范圍的大小；而確認(rèn)規(guī)避SAA空間分布的大小差異，則直接影響衛(wèi)星在軌工作的時(shí)序與時(shí)間。美國的 MSX衛(wèi)星在進(jìn)行在軌空間目標(biāo)信息收集的地面時(shí)序規(guī)劃時(shí)，將24 h內(nèi)衛(wèi)星經(jīng)過SAA的時(shí)間與持續(xù)時(shí)間加入任務(wù)表中進(jìn)行規(guī)避[10]。可見，SAA空間分布對(duì)衛(wèi)星歷經(jīng)SAA時(shí)間的影響研究可用于衛(wèi)星的在軌安全、任務(wù)時(shí)間規(guī)劃以及工作效率分析。

采用STK軟件仿真計(jì)算衛(wèi)星在6個(gè)月內(nèi)通過表1和表2中確定的SAA空間分布的時(shí)間占衛(wèi)星飛行時(shí)間的百分比。以天為單位考量SAA短期影響，圖3給出了衛(wèi)星在2014年1月1日24 h歷經(jīng)通量閾值為100 cm-2?s-1、能量閾值為50 MeV時(shí)所確定的SAA空間分布邊界的弧段。圓點(diǎn)代表衛(wèi)星過南極點(diǎn)時(shí)間，方框代表衛(wèi)星過SAA持續(xù)時(shí)間。由圖可知衛(wèi)星飛行15圈，其中有9圈歷經(jīng)SAA，從經(jīng)歷弧段的時(shí)間長短可以看出：衛(wèi)星歷經(jīng)SAA的時(shí)間各不相同，即在相同的空間分布邊界條件下，衛(wèi)星在同一天不同的時(shí)段受到SAA的影響不同。

圖3 衛(wèi)星24 h過SAA區(qū)域的弧段Fig.3 The segment of SAA-pass time in 24 hours

以月為單位統(tǒng)計(jì)SAA長期影響，在不同的通量閾值定義的SAA空間分布邊界條件下，衛(wèi)星在2014年1月至6月歷經(jīng)SAA的時(shí)間占衛(wèi)星飛行時(shí)間的百分比如圖4所示。結(jié)果表明：在相同的通量閾值條件下，衛(wèi)星歷經(jīng)SAA的時(shí)間百分比基本處在同一數(shù)值上，最大差異僅為0.13%，說明相同的區(qū)域在時(shí)間軸上對(duì)衛(wèi)星的影響基本保持不變；對(duì)不同的通量閾值條件下的時(shí)間百分比進(jìn)行比較，差異最大的是500 cm-2?s-1下的4.27%和100 cm-2?s-1下的8.21%之間，盡管二者的時(shí)間百分比差值僅約4%，但歷經(jīng)SAA時(shí)間卻縮短了48.0%。

圖4 衛(wèi)星經(jīng)歷不同通量閾值定義的SAA空間分布邊界時(shí)間百分比Fig.4 The SAA-pass time percentage of the satellite for different flux thresholds

通量閾值的大小由CCD相機(jī)的抗輻射性能決定，而衛(wèi)星表面的抗輻射保護(hù)則決定了衛(wèi)星能夠承受的高能粒子的能量閾值。在由不同能量閾值確立的SAA空間分布邊界條件下，衛(wèi)星歷經(jīng)SAA的時(shí)間百分比如圖5所示。能量閾值越高，則SAA空間分布的邊界越小，而受SAA影響的時(shí)間百分比越低，這與衛(wèi)星抗高能粒子輻射能力增強(qiáng)的趨勢(shì)一致。同一能量閾值條件下，歷經(jīng)SAA時(shí)間百分比趨于一致，最大差異僅為0.13%；不同能量閾值條件下，差異最大的是150 MeV下的5.60%與50 MeV下的8.21%之間，盡管二者的時(shí)間百分比差值僅有2.6%，但歷經(jīng)SAA時(shí)間卻縮短了約31.7%。

圖5 衛(wèi)星經(jīng)歷不同能量閾值定義的SAA空間分布邊界時(shí)間百分比Fig.5 The SAA-pass time percentage of the satellite for different energy thresholds

綜上分析結(jié)果表明，在航天器總體設(shè)計(jì)之初，提高衛(wèi)星有效載荷的通量閾值和能量閾值，對(duì)于降低SAA的影響和提高在軌任務(wù)工作效率具有顯著的效果。

3 結(jié)束語

本文通過對(duì)SAA空間高能粒子的通量與能量的仿真計(jì)算，分析了與通量閾值和能量閾值所對(duì)應(yīng)的SAA空間分布對(duì)衛(wèi)星的影響。通過計(jì)算科學(xué)實(shí)驗(yàn)衛(wèi)星經(jīng)歷這些SAA空間分布的時(shí)間百分比，定量地給出了 SAA空間分布對(duì)衛(wèi)星工作時(shí)間的影響，對(duì)衛(wèi)星安全評(píng)估與任務(wù)規(guī)劃具有指導(dǎo)意義。

本文雖然就SAA空間分布以及對(duì)低軌科學(xué)試驗(yàn)衛(wèi)星的影響給出了定量結(jié)果，但仍有許多工作需進(jìn)一步完善，例如：

1）需利用在軌試驗(yàn)數(shù)據(jù)對(duì)SAA空間分布的計(jì)量結(jié)果進(jìn)行驗(yàn)證；

2）需利用在軌試驗(yàn)數(shù)據(jù)開展SAA空間分布三維模型研究。

（References）

[1]張科科, 李宗耀, 胡海鷹, 等.基于天基光學(xué)探測(cè)圖像初析南大西洋異常區(qū)[J].中國空間科學(xué)技術(shù), 2015, 10(5): 33-39 ZHANG K K, LI Z Y, HU H Y, et al.Analysis of influence of south atlantic anomaly based on spacebased visible detection images[J].Chinese Space Science and Technology, 2015, 10(5): 33-39

[2]LIU C M, LI X J, GENG H B, et al.Incident particle range dependence of radiation damage in a power bipolar junction transistor[J].Chin Phys B, 2012, 21(10): 0104211

[3]謝倫, 濮祖蔭, 焦維新, 等.南大西洋異常區(qū)內(nèi)輻射帶高能質(zhì)子輻射環(huán)境長期變化的研究[J].中國科學(xué) D輯: 地球科學(xué), 2005, 35(7): 658-663

[4]文林, 李豫東, 郭旗, 等.質(zhì)子輻照導(dǎo)致科學(xué)級(jí)電荷耦合器件電離效應(yīng)和位移效應(yīng)分析[J].物理學(xué)報(bào), 2015, 64(2): 024220-1 WEN L, LI Y D, GUO Q, et al.Analysis of ionizing and department damage mechanism in proton-irradiationinduced scientific charge-coupled device[J].Acta Phys Sin, 2015, 64(2): 024220-1

[5]JIN J, WANG X Q, LIN S, et al.Effect of radiationinduced mean wavelength shift in optical fibers on the scale factor of interferometric fiber optic gyroscope at a wavelength of 1300 nm[J].Chin Phys B, 2012, 21(9): 094220

[6]VETTE J I.The AE-8 trapped electron model environment: NSSDCWDC-A-R&S 91-24[R].Greenbelt Maryland: NASA Goddard Space Flight Center, 1991

[7]SAWYER D M, VETTE J L.AP-8 trapped radiation environment for solar maximum and solar minimum: NSSDC/SDC-A-R&S 76-06[R].Greenbelt Maryland: NASA Space Flight Center, 1976

[8]趙慧杰, 何世禹, 孫彥錚, 等.100 keV質(zhì)子輻照對(duì)空間 GaAsPGe太陽電池光電效應(yīng)的影響[J].物理學(xué)報(bào), 2009, 58(1): 404-409ZHAO H J, HE S Y, SUN Y Z, et al.Effect of 100 keV proton irradiation on photo- emission of GaAs/Ge space solar cells[J].ACTA Physica Sinica, 2009, 58(1): 404-409

[9]SPENVIS[DB/OL].[2016-05-06].http://www.spenvis.oma.be/spenvis

[10]電離效應(yīng)[DB/OL].[2016-08-05].http://baike.baidu.eoln/view/3074227.htm

[11]程秀圍.質(zhì)子輻照對(duì)空間太陽望遠(yuǎn)鏡光學(xué)性能的影響[D].南京: 江蘇大學(xué), 2010: 11-25

[12]SUZUKI N, RUBIN D, LIDMAN C, et al.The Hubble Space Telescope cluster supernova survey: improving the dark-energy constraints abovez>1 and building an early-type-hosted supernova sample[J].The Astrophysical Journal, 2012, 746(10): 85

[13]SHARMA J, WISEMAN A, ZOLLINGER G.Improving space surveillance with space-based visible sensor[R].Massachusetts Inst of Tech Lexington Lincoin Lab, 2001

[14]路偉.Geant4的開發(fā)及其在空間輻射效應(yīng)分析中的應(yīng)用[D].長沙: 國防科學(xué)技術(shù)大學(xué), 2007

（編輯：肖福根）

The influence of south atlantic anomaly boundary on space-based observation satellite

LI Zongyao1, ZHANG Keke1,2, HU Haiying1, SHENG Lei1
(1.Shanghai Engineering Center for Microsatellites, Chinese Academy of Sciences, Shanghai 201203, China; 2.University of Chinese Academy of Sciences, Beijing 100039, China)

By simulating the South Atlantic Anomaly (SAA), the ranges of the SAA with different flux thresholds and energy thresholds are obtained.As an example, the time percentage of the space-based observation satellite that passes the SAA region is calculated in 24 consecutive hours and 6 months, and the influence of the SAA field on the operation time is evaluated.The results show that, for the 50 MeV energy threshold, with the maximum influence of the flux threshold on the spatial domain, the latitude difference is about 23.35%, the longitude difference is about 27.06%, and the SAA pass time difference is about 48% between 500 cm-2?s-1and 100 cm-2?s-1flux thresholds.For the 100 cm-2?s-1flux threshold, the latitude difference is about 13.86%, the longitude difference is about 17.29%, and the SAA -pass time difference is about 31.7% between 150 MeV and 50 MeV energy thresholds.

SAA; space-based observation; energetic protons; time percentage

V520.6

：A

：1673-1379(2016)05-0479-05

10.3969/j.issn.1673-1379.2016.05.004

李宗耀（1988—），男，碩士學(xué)位，研究方向?yàn)樾l(wèi)星總體載荷的研究與設(shè)計(jì)；E-mail: zyli_microsat@163.com。通信作者：張科科（1980—），男，碩士學(xué)位，副研究員，研究方向?yàn)槲⑿⌒l(wèi)星總體設(shè)計(jì)；E-mail: 13816686945@139.com。

2016-03-17；

：2016-09-05

中國科學(xué)院國防創(chuàng)新基金項(xiàng)目“基于衛(wèi)星網(wǎng)絡(luò)的分布式星上自主任務(wù)規(guī)劃研究”（編號(hào)：CXJJ-155085）；上海市科學(xué)技術(shù)委員會(huì)揚(yáng)帆計(jì)劃“雙星編隊(duì)望遠(yuǎn)鏡衛(wèi)星的高精度協(xié)同智能控制研究” 項(xiàng)目（編號(hào)：YF141413600）