高著秀，孫 健，王玉林，李羅鋼，張耀磊

?

彈道導(dǎo)彈高能粒子輻射環(huán)境適應(yīng)性研究

高著秀，孫 健，王玉林，李羅鋼，張耀磊

（中國運(yùn)載火箭技術(shù)研究院研究與發(fā)展中心，北京，100076）

導(dǎo)彈飛行環(huán)境中包含多種高能帶電粒子，會(huì)對(duì)敏感電子器件產(chǎn)生輻射環(huán)境效應(yīng)影響。針對(duì)3個(gè)典型的彈道獲取高能粒子能譜分布情況，開展輻射環(huán)境及效應(yīng)的分析，探討導(dǎo)彈對(duì)高能粒子輻射環(huán)境的適應(yīng)性，提出防護(hù)設(shè)計(jì)建議，以改善其適應(yīng)性。

彈道導(dǎo)彈；輻射效應(yīng)；環(huán)境適應(yīng)性

0 引 言

彈道導(dǎo)彈飛行任務(wù)環(huán)境中包含輻射帶電子、輻射帶質(zhì)子、太陽宇宙線質(zhì)子等[1~3]多種高能帶電粒子。相對(duì)于長期在軌的衛(wèi)星，導(dǎo)彈飛行時(shí)間較短，一般忽略輻射環(huán)境的影響，因此，也鮮有對(duì)彈道導(dǎo)彈輻射環(huán)境影響的相關(guān)分析和報(bào)道。隨著對(duì)各類飛行器空間輻射環(huán)境效應(yīng)研究的深入，發(fā)現(xiàn)即使飛行任務(wù)時(shí)間短，軌道相對(duì)較低，也會(huì)受到高能粒子的輻射影響。尤其是單粒子效應(yīng)，是一種概率事件，飛行過程中敏感器件受到高能粒子入射時(shí)，達(dá)到其單粒子效應(yīng)閾值即有可能發(fā)生故障，影響任務(wù)安全。另外，導(dǎo)彈結(jié)構(gòu)及材料更加多樣化，電子器件特征尺寸不斷減小、集成度不斷提高，加劇了發(fā)生單粒子事件的風(fēng)險(xiǎn)[4~6]。因此，有必要針對(duì)彈道導(dǎo)彈開展輻射環(huán)境適應(yīng)性的研究，為導(dǎo)彈研制提供輻射環(huán)境效應(yīng)防護(hù)方面的參考。本文針對(duì)3條飛行軌跡經(jīng)過相對(duì)惡劣的輻射環(huán)境的典型彈道導(dǎo)彈任務(wù)開展分析，獲取高能粒子能譜分布情況，計(jì)算在軌吸收劑量及翻轉(zhuǎn)頻次，總結(jié)影響輻射環(huán)境效應(yīng)的因素及相關(guān)規(guī)律，提出防護(hù)設(shè)計(jì)建議。

1 輻射環(huán)境能譜分析

在彈道導(dǎo)彈飛行過程中經(jīng)過的輻射環(huán)境中的帶電粒子主要包括：地球輻射帶質(zhì)子及電子，太陽宇宙射線質(zhì)子及重離子以及銀河宇宙射線質(zhì)子和重離子。其中，對(duì)單粒子效應(yīng)起作用的帶電粒子包括：地球輻射帶質(zhì)子，太陽宇宙線質(zhì)子和重離子，銀河宇宙線重離子；對(duì)總劑量效應(yīng)起作用的帶電粒子包括：地球輻射帶質(zhì)子和電子、太陽宇宙線質(zhì)子。

彈道導(dǎo)彈飛行一般經(jīng)過內(nèi)輻射帶區(qū)域，對(duì)于輻射帶高能粒子，大于10 MeV的高能質(zhì)子基本分布于磁緯度絕對(duì)值小于50°的空間[7]。高度在2 000 km以下時(shí)，輻射帶質(zhì)子主要分布于南大西洋異常區(qū)，輻射帶電子主要分布于南大西洋異常區(qū)和緯度約60°的兩個(gè)帶狀區(qū)域。對(duì)于太陽宇宙線和銀河宇宙線高能粒子主要集中于地磁屏蔽作用相對(duì)較弱的極區(qū)，因此，本文選取3條彈道進(jìn)行分析：a）彈道一，飛行軌跡基本處于南大西洋覆蓋范圍內(nèi)；b）彈道二，飛行軌跡在北緯60o附近；c）彈道三，穿過極區(qū)。3條彈道高度最高點(diǎn)均在2 000 km附近，飛行時(shí)間均約40 min。以上3條彈道均為輻射環(huán)境相對(duì)惡劣的區(qū)域，針對(duì)以上3條典型彈道開展分析，基本可以說明彈道導(dǎo)彈在輻射環(huán)境中的適應(yīng)性。

下文給出3條彈道輻射環(huán)境中地球輻射帶電子、輻射帶質(zhì)子、太陽宇宙線及銀河宇宙線質(zhì)子和重離子積分能譜，如圖1~7所示。

圖1為3條彈道輻射環(huán)境中地球輻射帶電子微分能譜分布。輻射帶電子環(huán)境模型采用工程中常用的AE8模型，考慮最惡劣情況，采用太陽活動(dòng)高年模式 AE8MAX，由于構(gòu)造AE8模型時(shí)所用的地磁場模型為IGRF國際參考地磁場模式，因此本文中選用IGRF作為地磁場模式[8]，該模型是一系列描述地球主磁場及其變化規(guī)律的數(shù)學(xué)模型[9,10]。

圖1 地球輻射帶電子微分能譜

從圖1中分析可知，對(duì)于能量小于3 MeV的輻射帶電子，經(jīng)過南大西洋的彈道一中的輻射帶電子通量遠(yuǎn)高于其它兩條彈道。對(duì)于能量大于3 MeV的高能輻射帶電子，穿過高緯區(qū)域的彈道二中的電子通量最高。

圖2為3條彈道的地球輻射帶質(zhì)子微分能譜分布。輻射帶質(zhì)子環(huán)境模型采用AP8MIN，輻射帶質(zhì)子通量在太陽活動(dòng)低年時(shí)較大，因此選太陽活動(dòng)低年時(shí)模型作為輻射惡劣情況，地磁場模型選用IGRF國際參考地磁場模式。

圖2 地球輻射帶質(zhì)子微分能譜

從圖2中分析可知，對(duì)于能量大于1 MeV的輻射帶質(zhì)子，經(jīng)過南大西洋的彈道一中的輻射帶質(zhì)子通量遠(yuǎn)高于其它兩條彈道。對(duì)于能量小于1 MeV的輻射帶質(zhì)子，穿過高緯區(qū)域的彈道二中的質(zhì)子通量最高?？傮w上，穿過極區(qū)的彈道三輻射環(huán)境中輻射帶質(zhì)子和電子通量均較小。

圖3為彈道二和彈道三輻射環(huán)境中太陽宇宙線質(zhì)子微分能譜。圖3a環(huán)境模型選用1989峰通量模式，地磁條件為非地影區(qū)、地磁風(fēng)暴模式；圖3b環(huán)境模型選用國際通用的ESP平均統(tǒng)計(jì)模型，該模型可用于總劑量和最劣事件劑量的預(yù)測(cè)[11]。地磁條件為非地影區(qū)、地磁風(fēng)暴模式。穿過南大西洋的彈道一輻射環(huán)境中太陽宇宙線質(zhì)子通量在平均模式和峰通量模式下均為零。穿過極區(qū)的彈道三輻射環(huán)境中太陽宇宙線質(zhì)子通量略高于穿過60°高緯地區(qū)的彈道二輻射環(huán)境中太陽宇宙線質(zhì)子通量。對(duì)比可知，峰通量模式下太陽宇宙線質(zhì)子積分通量比平均模式下約高3個(gè)量級(jí)，可見不同的太陽活動(dòng)情況下，太陽質(zhì)子的通量變化顯著。

a）峰通量模式

b）平均模式

圖3 太陽宇宙線質(zhì)子微分能譜

圖4~6為3條彈道輻射環(huán)境中太陽宇宙線重離子積分能譜，其中峰通量模式采取的環(huán)境模型為CREME96最壞5 min模式，地磁條件為非地影區(qū)、風(fēng)暴模式。平均模式采取的環(huán)境模型為PSYCHIC模式，該模型具有較寬的能量范圍[12]。地磁條件為非地影區(qū)、地磁風(fēng)暴模式。其中彈道一輻射環(huán)境中太陽宇宙線重離子只有峰通量模式下存在少量通量。3條彈道中，經(jīng)過極區(qū)的太陽宇宙線重離子通量最大。對(duì)比可知，峰通量模式下太陽宇宙線重離子積分通量比平均模式下高約2個(gè)數(shù)量級(jí)。

圖4 彈道一太陽宇宙線重離子積分能譜（峰通量模式）

a）峰通量模式

b）平均模式

圖5 彈道二太陽宇宙線重離子積分能譜

a）峰通量模式

b）平均模式

續(xù)圖6

圖7為3條彈道輻射環(huán)境中銀河宇宙線重離子積分能譜（=2~28）。

a）彈道一

b）彈道二

c）彈道三

由圖7可知，分析中采用的環(huán)境模型為CREME96模式[13]，是目前國內(nèi)外廣泛采用的銀河宇宙線統(tǒng)計(jì)模型，該模型利用解析表達(dá)式描述銀河宇宙線粒子的強(qiáng)度變化。地磁條件為非地影區(qū)、地磁平靜模式（該模式下銀河宇宙線通量較大）。穿過極區(qū)的銀河宇宙線重離子通量最高，穿過南大西洋異常區(qū)的銀河宇宙線重離子通量最低。銀河宇宙線重離子通量比以上各種粒子均低，其對(duì)飛行器的輻射環(huán)境效應(yīng)影響相對(duì)較小。

總結(jié)上述對(duì)3條彈道輻射環(huán)境中高能帶電粒子的分析結(jié)果，彈道一考核輻射帶質(zhì)子和輻射帶電子的影響，彈道二考核能量和通量均相對(duì)較高的各種粒子輻射環(huán)境的影響，彈道三重點(diǎn)考核太陽宇宙線和銀河宇宙線輻射環(huán)境的影響。具體特性見表1。

表1 彈道輻射環(huán)境特性統(tǒng)計(jì)

2 輻射環(huán)境適應(yīng)性分析

2.1 總劑量效應(yīng)適應(yīng)性分析

輻射環(huán)境中吸收的劑量主要來源于地球輻射帶質(zhì)子、輻射帶電子和太陽宇宙線質(zhì)子。根據(jù)第1節(jié)對(duì)3條彈道輻射環(huán)境中獲得的帶電粒子能譜的分析結(jié)果，可以獲得3條彈道的吸收劑量-深度分布曲線。

圖8為3條彈道飛行任務(wù)中的吸收劑量-深度分布曲線，分析時(shí)采用實(shí)心球模型。由圖8可知，穿過南大西洋異常區(qū)的彈道一輻射劑量遠(yuǎn)高于其它兩條彈道，且彈道一輻射劑量主要來源于輻射帶質(zhì)子，彈道二、彈道三的輻射劑量主要來源與輻射帶電子和太陽宇宙線質(zhì)子。假設(shè)只有電子設(shè)備殼體對(duì)電子器件提供防護(hù)，設(shè)備殼體厚度為2 mm，材料為鋁合金。對(duì)照 圖8中的劑量深度曲線，3條彈道中電子設(shè)備內(nèi)部在經(jīng)歷彈道輻射環(huán)境后總的吸收劑量分別為15.1 rad（Si）、2.26 rad（Si）和0.33 rad（Si），常規(guī)電子器件抗總劑量閾值大于1 000 rad（Si），可以滿足任務(wù)要求，電子設(shè)備發(fā)生總劑量效應(yīng)的風(fēng)險(xiǎn)極低，無需采取專門的總劑量效應(yīng)防護(hù)措施。

a）彈道一

b）彈道二

c）彈道三

2.2 單粒子效應(yīng)適應(yīng)性分析

為開展電子設(shè)備內(nèi)部敏感電子器件單粒子效應(yīng)計(jì)算分析，首先應(yīng)基于空間高能帶電粒子能譜及LET譜的計(jì)算。針對(duì)空間高能帶電重離子需計(jì)算其LET譜，而高能質(zhì)子需計(jì)算微分能譜，下面給出3條彈道輻射環(huán)境高能粒子穿過防護(hù)層后的能譜及LET譜。

同樣針對(duì)2 mm等效鋁防護(hù)層開展分析，由于只有穿過南大西洋異常區(qū)的彈道一輻射環(huán)境中存在能量較高的輻射帶質(zhì)子，因此，只有彈道一輻射帶質(zhì)子穿過2 mm防護(hù)層，其它兩條彈道輻射環(huán)境中輻射帶質(zhì)子能量較小，其未能穿過2 mm防護(hù)層。圖9為彈道一輻射帶質(zhì)子穿過防護(hù)層后能譜。

圖9 彈道一輻射帶質(zhì)子穿過防護(hù)層后能譜

圖10為彈道二和彈道三輻射環(huán)境中太陽宇宙線質(zhì)子穿過2 mm厚防護(hù)層后微分能譜圖（彈道一無太陽宇宙線質(zhì)子通量）。圖10a環(huán)境模型采用平均統(tǒng)計(jì)模型ESP模式；圖10b環(huán)境模型采用1989峰通量模式。

a）平均模式

b）峰通量模式

圖10 太陽宇宙線質(zhì)子穿過防護(hù)層后能譜

由圖10分析結(jié)果可知，太陽宇宙線質(zhì)子穿過2 mm等效鋁防護(hù)層后，在通量峰通量模式下比平均統(tǒng)計(jì)模式約高3個(gè)數(shù)量級(jí)。

圖11為彈道二和彈道三輻射環(huán)境中太陽宇宙線重離子穿過2 mm厚防護(hù)層后LET譜。其中，圖11a離子模型采用PSYCHIC模式；圖11b離子模型采用CREME96 最壞5 min模式。由圖11的分析結(jié)果可知，太陽宇宙線重離子穿過2 mm等效鋁防護(hù)層后，在峰通量模式下比平均統(tǒng)計(jì)模式高約3個(gè)數(shù)量級(jí)。彈道一輻射環(huán)境中只有在峰通量模式下有少量太陽宇宙線重離子穿過防護(hù)層。

a）平均模式

b）峰通量模式

圖11 太陽宇宙線重離子穿過防護(hù)層后LET譜

圖12為3條彈道輻射環(huán)境中銀河宇宙線重離子穿過2 mm厚防護(hù)層后LET譜圖。其中，離子環(huán)境模型采用CREME96模式。彈道二和彈道三輻射環(huán)境中銀河宇宙線重離子穿過防護(hù)層后LET譜基本一致，且遠(yuǎn)高于穿過南大西洋異常區(qū)的彈道一輻射環(huán)境中銀河宇宙線重離子穿過防護(hù)層后LET譜。

圖12 銀河宇宙線重離子穿過防護(hù)層后LET譜

為了對(duì)電子器件抗輻射加固選型提供參考，本文以兩款SRAM型FPGA器件為代表，對(duì)其任務(wù)期間發(fā)生單粒子效應(yīng)的頻次進(jìn)行分析。在軌單粒子翻轉(zhuǎn)概率的計(jì)算結(jié)果見表2，分析時(shí)采取太陽活動(dòng)峰年時(shí)模式，考慮發(fā)生太陽耀斑時(shí)的粒子情況，太陽宇宙線質(zhì)子與重離子都采取平均統(tǒng)計(jì)模型和“5 min峰通量”模式進(jìn)行比較。

表2 單粒子翻轉(zhuǎn)概率分析

依據(jù)工程設(shè)計(jì)經(jīng)驗(yàn)，為確保任務(wù)安全，在單粒子效應(yīng)防護(hù)中應(yīng)進(jìn)行3倍冗余設(shè)計(jì)，所分析的3個(gè)彈道的飛行任務(wù)中器件發(fā)生單粒子效應(yīng)的頻次應(yīng)小于12次/（器件·天）。器件1在3個(gè)彈道的飛行任務(wù)中，發(fā)生單粒子翻轉(zhuǎn)的次數(shù)均大于12次/（器件·天）。器件2在彈道二、彈道三飛行任務(wù)中，平均統(tǒng)計(jì)模式下在軌翻轉(zhuǎn)次數(shù)小于12次/（器件·天），峰通量模式下在軌翻轉(zhuǎn)次數(shù)均大于12次/（器件·天）。彈道一飛行任務(wù)中，無論平均統(tǒng)計(jì)模式還是峰通量模式在軌發(fā)生單粒子翻轉(zhuǎn)次數(shù)均大于12次/（器件·天）。可見器件2的抗輻射性能優(yōu)于器件1，這兩款器件單粒子翻轉(zhuǎn)概率均過高，不滿足抗輻射要求。

對(duì)比分析結(jié)果可知：

a）對(duì)于輻射帶質(zhì)子，應(yīng)主要考慮彈道是否穿過南大西洋異常區(qū)，穿過該區(qū)域時(shí)，如果防護(hù)層厚度較小，有發(fā)生單粒子效應(yīng)的風(fēng)險(xiǎn),其它情況下輻射帶質(zhì)子能量較小，穿過防護(hù)層后能量及通量更小，基本不會(huì)造成單粒子翻轉(zhuǎn)。

b）對(duì)于太陽宇宙線及銀河宇宙線，在極區(qū)的通量相對(duì)較高。另外，發(fā)生單粒子效應(yīng)的概率取決于是否爆發(fā)太陽活動(dòng)及爆發(fā)太陽活動(dòng)的強(qiáng)度，發(fā)生太陽耀斑時(shí)比平均模型下的計(jì)算結(jié)果約高3個(gè)數(shù)量級(jí)。

c）由表1中對(duì)兩款器件的分析結(jié)果比較可知，器件的抗輻射水平直接影響著單粒子效應(yīng)的風(fēng)險(xiǎn)，應(yīng)選擇滿足任務(wù)抗輻射指標(biāo)要求的電子器件。

3 輻射防護(hù)設(shè)計(jì)建議

綜合以上分析結(jié)果，結(jié)合衛(wèi)星常規(guī)輻射防護(hù)工程經(jīng)驗(yàn)，針對(duì)彈道導(dǎo)彈總劑量效應(yīng)和單粒子效應(yīng)防護(hù)提出以下設(shè)計(jì)建議：

a）總劑量效應(yīng)防護(hù)設(shè)計(jì)建議。

對(duì)于彈道導(dǎo)彈，即使單機(jī)外部沒有其他防護(hù)結(jié)構(gòu)，僅有單機(jī)自身約2 mm殼體的屏蔽情況下，單機(jī)在軌吸收的輻射劑量即達(dá)到非常小的量級(jí)，目前的結(jié)構(gòu)設(shè)計(jì)可以滿足防護(hù)要求，不需要進(jìn)行額外的總劑量效應(yīng)防護(hù)設(shè)計(jì)。

b）單粒子效應(yīng)防護(hù)設(shè)計(jì)建議。

對(duì)于彈道導(dǎo)彈單粒子效應(yīng)的防護(hù)可參照衛(wèi)星常用單粒子效應(yīng)防護(hù)設(shè)計(jì)要求。結(jié)合器件的抗輻射指標(biāo)及環(huán)境參數(shù)，對(duì)彈道導(dǎo)彈電子器件抗單粒子效應(yīng)提出防護(hù)設(shè)計(jì)要求如下：

1）器件抗SEU（單粒子翻轉(zhuǎn)）的LET閾值大于15 （MeV·cm2）/mg時(shí)可直接使用，器件抗SEL（單粒子鎖定）閾值大于75 （MeV·cm2）/mg時(shí)可直接使用；

2）器件抗SEU的LET閾值小于15 （MeV·cm2）/mg或無LET閾值數(shù)據(jù)的，應(yīng)進(jìn)行抗SEU防護(hù)設(shè)計(jì)，并對(duì)其防護(hù)效果進(jìn)行評(píng)估，確認(rèn)滿足要求方可使用；器件抗SEL閾值在37~75 （MeV·cm2）/mg時(shí)，需進(jìn)行防護(hù)設(shè)計(jì)，對(duì)其防護(hù)效果評(píng)估后，方可使用；

3）器件抗SEL閾值小于37 （MeV·cm2）/mg時(shí)，不得使用。

4 結(jié)束語

本文針對(duì)3個(gè)典型的彈道任務(wù)開展了輻射環(huán)境及效應(yīng)的分析，這3條彈道分別穿過輻射環(huán)境中各種高能粒子相對(duì)惡劣的環(huán)境，分析結(jié)果基本可以說明彈道導(dǎo)彈在輻射環(huán)境中的適應(yīng)性。依據(jù)分析結(jié)果，在對(duì)彈道導(dǎo)彈輻射防護(hù)設(shè)計(jì)中，基本可以不考慮總劑量效應(yīng)對(duì)電子器件及材料的影響；對(duì)于單粒子效應(yīng)，當(dāng)防護(hù)層厚度有限時(shí)，應(yīng)選用抗輻射性能滿足要求的電子器件，以確保在爆發(fā)大的太陽活動(dòng)事件條件下，也不會(huì)因產(chǎn)生單粒子效應(yīng)而發(fā)生故障。

[1] 焦維新, 田天. 空間環(huán)境研究的現(xiàn)狀與展望[J]. 航天器環(huán)境工程, 2007, 24(6): 337-340.

[2] 劉必鎏, 楊平會(huì), 蔣孟龍, 等. 航天器單粒子效應(yīng)的防護(hù)研究[J]. 航天器環(huán)境工程, 2010, 27(6): 693-697.

[3] 黃本城, 童靖宇. 空間環(huán)境工程學(xué)[M]. 北京: 中國科學(xué)技術(shù)出版社, 2010.

[4] 蔡震波. 新型航天器抗輻射加固技術(shù)的研究重點(diǎn)[J]. 航天器環(huán)境工程, 2010, 27(2): 173-176.

[5] 丁義剛. 空間輻射環(huán)境單粒子效應(yīng)研究[J]. 航天器環(huán)境工程, 2007, 24(5): 283-290.

[6] 馮穎, 劉忠健. 單粒子效應(yīng)對(duì)飛行器的影響分析及防護(hù)技術(shù)[J]. 強(qiáng)度與環(huán)境, 2011, 38(1): 26-30.

[7] 陶宏任, 林瑞淋, 師立勤, 劉四清, 張效信, 龔建村. 質(zhì)子輻射帶輻射中心區(qū)域模型[J]. 空間科學(xué)學(xué)報(bào), 2015, 35(3): 293-305.

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

[9] 郭才發(fā). 空間地磁場應(yīng)用的若干關(guān)鍵技術(shù)研究[D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 2014.

[10] 白春華, 徐文耀, 康國發(fā). 地球主磁場模型[J]. 地球物理學(xué)進(jìn)展, 2008, 23(4): 1045-1057.

[11] 沈自才, 閆德葵. 空間輻射環(huán)境工程的現(xiàn)狀及發(fā)展趨勢(shì)[J]. 航天器環(huán)境工程, 2014, 31(3): 229-240.

[12] 余迎慶, 王詠梅, 付利平, 孔春陽, 李小銀, 王月. 太陽同步軌道真空紫外光電倍增管的高能粒子屏蔽方法研究[J]. 空間科學(xué)學(xué)報(bào), 2012, 32(4): 532-536.

[13] 馮彥君, 華更新, 劉淑芬. 航天電子抗輻射研究綜述[J]. 宇航學(xué)報(bào), 2007, 28(5): 1071-1080.

Research of Energetic Particles Environment Adaptability of Ballistic Missile

Gao Zhu-xiu, Sun Jian, Wang Yu-lin, Li Luo-gang, Zhang Yao-lei

(R&D Center, China Academy of Launch Vehicle Technology, Beijing, 100076)

There are various energetic particles in flight environment of missile. They will affect the sensitive electronic apparatus. Based on the energy spectrum of energetic particles along three trajectories, the radiation environment and effects are analyzed in this article. Also, the energetic particles environment adaptability of ballistic missiles is studied, and correlative protection design suggestions are proposed.

Ballistic missile; Radiation effects; Environmental adaptability

1004-7182(2017)06-0106-07

10.7654/j.issn.1004-7182.20170623

V412.1

A

2016-12-20；

2017-06-08

高著秀（1977-），女，博士，工程師，主要研究方向?yàn)榭臻g環(huán)境總體設(shè)計(jì)