中短期任務的太陽質子事件通量預報研究

崔延美，師立勤，劉四清

（中國科學院國家空間科學中心，北京 100190）

中短期任務的太陽質子事件通量預報研究

崔延美，師立勤，劉四清

（中國科學院國家空間科學中心，北京 100190）

太陽質子事件通量的預測對航天器的抗輻射加固設計和宇航員的出艙活動具有重要意義。針對1年以下的中短期航天任務，對太陽活躍年和太陽平靜年分別統計了太陽質子事件和大于10 MeV質子事件通量的發生概率，分析得到太陽質子事件通量分布基本符合對數正態分布。在此基礎上，計算出了一定置信度下1年以下不同航天任務期內的質子事件通量分布，為執行中短期航天任務提供了太陽質子事件通量預測的依據。

太陽質子事件；太陽質子通量；航天任務；預報

1 引言

太陽質子事件是最具破壞性的空間天氣現象之一。太陽上發生爆發活動時，會噴射出大量的高能帶電粒子，它們以每秒幾千、幾萬、甚至十幾萬千米的速度，最快十幾分鐘就可以到達地球，使地球周圍的高能帶電粒子數量增加數千倍，甚至上萬倍。由于絕大部分的粒子是質子，因此把這種事件稱之為太陽質子事件。

太陽質子事件威脅著地球磁層和大氣層外航天器的安全，使衛星上的電子元器件失效，還會對宇航員造成傷害，嚴重影響著人類的航天活動。因此，在航天器的設計階段或者在人類從事太空活動任務期間，都需要提前預知太陽質子事件的發生狀況，以采取相應措施規避不必要的風險。然而，太陽質子事件是一種小概率事件，在其產生、加速機理還沒完全確定的條件下，對其進行預報是十分困難的。

對于持續1年以上的航天任務，當前國際上應用比較普遍的質子事件通量模型有King模型［1］、JPL模型［2-3］、ESP［4-9］模型及MSU模型［10］。King模型是建立的第一個太陽高能質子通量的經驗統計模型。模型提出后取得了一定的成效，在很長一段時間被視作執行航天任務期間太陽質子積分通量的標準預報模型。但該模型只采用第20太陽活動周的活躍年份數據，統計數據較少；且模型中對太陽黑子數和年均太陽質子累計通量呈線性相關的假設并不成立。JPL模型解決了King模型存在的問題，它采用了近三個太陽活動周的太陽質子事件數據。在統計分析的基礎上，JPL模型中應用了以下假設：太陽質子事件的發生是隨機的，事件的爆發概率滿足泊松統計分布；質子事件通量服從對數正態分布；太陽質子事件只發生在太陽活動周的連續7個活躍年中，且爆發頻率一致。JPL模型對持續時間超過一年的航天任務能預報不同能量范圍的質子積分通量水平。JPL模型的第二個版本JPL-91是目前公認較好的太陽質子事件長期預報模型，被推薦用于航天任務的設計規劃［11］。ESP模型采用了更充分的觀測數據作為統計基礎，不僅預測太陽質子積分通量，還能夠給出最惡劣情況的太陽質子事件。MSU模型則是由莫斯科州立大學開發的，模型的一個版本《太陽高能粒子通量和峰值流量的概率模型第一部分：質子》已經是俄羅斯國家標準（GOST R 25645.165-2001），并被推薦為國際標準（ISO TS 15391）［12］。在總結吸取國外模型方法的基礎上，李婷婷等人［13］利用人造衛星開始用于探測太陽高能粒子以來至2007年4個太陽活動周的數據，探討了太陽活動性對質子事件發生概率的影響因素，評估了不同太陽活動水平下的質子通量，建立了新模型。

對于持續1年以下的中短期航天任務，國內外還沒有建立起太陽質子事件通量預測的應用模型。在考慮太陽質子事件通量時，通常借鑒歷史上曾經發生的最惡劣事件或者某一次普通質子事件進行分析，缺乏統計依據。本文將利用1967—2014年期間的太陽高能粒子觀測數據，統計1年以下不同時間段中的太陽質子事件的發生，給出1—365天不同時間段的太陽質子事件和質子通量的發生概率。在此基礎上，初步建立起中短期太陽質子事件通量模型，從而為中短期航天任務提供太陽質子事件通量預測的依據。

2 數據處理

我們采用的太陽質子數據主要來自OMNI網站［14］。OMNI網站整合了1967年以來的1小時精度的高能粒子觀測數據，其中2005年之前的數據來自IMP衛星的觀測，2006年之后的數據來自GOES衛星的觀測。本文選用的數據時間段為1967—2014年，能段為＞10 MeV。缺值數據時，采用插值法進行補充。由于2001年—2003年中有大段數據缺值，我們采用GOES衛星觀測的數據［15］進行了補充。

Feynman等人［2］定義：當＞10 MeV的質子通量日均值大于1 PFU時就認為發生了質子事件；當＞10 MeV的質子通量日均值連續兩天低于1 PFU時則認為該質子事件結束。本研究中，由于要考慮的任務期會很短，因此不考慮質子事件的持續性，只對＞10 MeV的質子通量日均值大于1 PFU的時段進行太陽質子通量的計算，不滿足該條件的時段則認為沒有質子事件發生，太陽質子通量計為0。同時，我們將太陽活躍年和太陽平靜年分別進行考慮。太陽活躍年指一個太陽活動周中極大年的前2.5年到后4.5年，活動周中的其余4個年份則為太陽平靜年。在我們的數據統計樣本中，考慮的太陽活躍年包括1967—1972年，1977—1983年，1987—1993年，1998—2004年，2010—2014年，共計32年；其余16年則為太陽平靜年。總體來看，在太陽活躍年，質子事件發生天數多，太陽質子年積分通量高；在太陽平靜年，質子事件發生天數少甚至是不發生，太陽質子年積分通量低（圖1）。

圖1 太陽質子事件發生與太陽活動周的關系Fig.1 The relationship between the production of solar proton events and the solar cycle

3 統計分析

從經驗上來看，任務期越長，遭遇到質子事件的概率越大，可能遭遇到的太陽質子通量越高。而且，在太陽活躍年遭遇到質子事件的概率要高于太陽平靜年。對樣本中所有太陽活躍年和太陽平靜年，我們分別統計分析了太陽質子事件和太陽質子通量的發生概率。

通過時間窗口滑動，可以得到不同任務期的質子事件發生狀態和太陽質子積分通量。具體方法是：假定總樣本天數為N，當任務期為n（1，2，3，4，……365）天時，進行逐日的時間窗口滑動，滑動步長為n-1，可以得到N-n+1組任務期為n天的樣本。在任務期n天內只要有質子事件發生就認定發生了質子事件，對N-n+1組樣本中的太陽質子日積分通量進行累計，可得到n天任務期的太陽質子通量的數值。在此基礎上，統計N-n+1組樣本中發生質子事件的樣本所占地比例，就得到了n天任務期中發生質子事件的概率；以太陽質子通量的數值為樣本，統計超過以一系列通量值為閾值的樣本所占比例，就得到n天任務期中太陽質子通量的發生概率。例如，在23太陽活動周的活躍年1998—2004年中，共有2557天，當任務期n=2天時，通過逐日的時間窗口滑動，得到N-n+1=2557 -2+1=2556組樣本，對這些樣本分別進行太陽質子日通量的累計，就得到2556組任務期為2天的太陽質子通量的數值。其中，有504組樣本有質子事件發生，則發生質子事件的概率為504/2556= 0.197。對太陽質子通量數據，若取107p·cm-2為閾值，則超過該通量值的太陽質子通量數據有226組，得到在2天的任務期中，發生超過107p·cm-2質子通量的概率為：226/2556=0.088。

圖2分別展示不同任務期的質子事件發生率：（a）為太陽活躍年中太陽質子事件發生率隨任務期時間的變化；（b）為太陽平靜年中質子事件發生率隨任務期時間的變化。由圖可見，隨著任務期時間的延長，太陽質子事件發生率逐漸增大。在相同時間長度的任務期中，太陽活躍年中的質子事件發生率要遠高于太陽平靜年。例如，任務期20天時，太陽活躍年中的質子事件發生率約為0.5；太陽平靜年中的質子事件發生率僅約為0.15。在太陽活躍年中，當任務期為100天時，太陽質子事件的發生率已超過了0.9，隨著任務期時間的進一步延長，質子事件的發生率逐步趨近于1，并最終等于1。在太陽活動平靜年，即使對于所考慮的最長任務期365天，太陽質子事件的發生率最大也僅為0.74。

圖2 不同任務期的太陽質子事件發生率Fig.2 The production probability of solar proton events with different space mission periods from 1day to 365 days

質子事件發生率隨任務期時間的散點分布可由指數函數進行擬合，擬合公式如式（1）：

式中y0、A、t分別為指數函數的偏移、幅度與衰減常數。對應太陽活躍年和太陽平靜年，擬合參數的各取值見表1，擬合結果見圖2。

參數太陽活躍年太陽平靜年y00.9920.818 A -0.822-0.750 t -40.588-169.557

表1 太陽質子事件發生率進行指數擬合的參數Table 1 The fitting parameter values of exponential functions for the production probability of solar proton events

圖3給出了在太陽活躍年和太陽平靜年中，任務期時間分別為1天、2天、5天、10天、20天、30天、50天和100天時，太陽質子通量的發生概率。從圖中可以看到，隨著通量值的提高，發生大的太陽質子通量的概率逐漸降低。隨著任務期的時間延長，以相同通量為閾值的太陽質子通量的概率逐漸增大。在相同的任務期中，太陽活躍年中的質子通量發生概率要遠高于太陽平靜年。如，在太陽活躍年，任務期為20天時，發生超過1 ×107p·cm-2質子通量的可能性為0.333，當通量閾值提高到1×109p·cm-2時，發生超過該閾值的太陽質子通量的可能性下降至0.055。當任務期為30天時，發生超過1×107p·cm-2和1×109p·cm-2的太陽質子通量的可能性分別為0.443和0.084，均高于任務期為20天的情況。當在太陽平靜年，任務期為20天時，發生超過1×107p·cm-2和1×109p·cm-2的太陽質子通量的可能性分別為0.102和0.018，均低于太陽活躍年任務期20天的情況。

值得注意的是，并非隨著任務期時間的延長，發生大的太陽質子通量的概率可以一直增大下去。到達一定通量閾值后，太陽質子通量的發生概率會存在一個截止值。如在太陽平靜年，任務期時間從20天到100天，均未出現超過5×109p·cm-2的太陽質子通量發生，見圖3（b）。

4 模型建立

在實際航天任務的評估中，通常要求給出一定置信度下任務期的質子通量。置信水平是以概率分布函數為基礎的統計量。95%的置信度意味著95%的概率不會出現超過給定質子通量的情況，也就是只有5%的可能出現超過給定質子通量的情況，80%的置信度也就意味著只有20%的可能出現超過給定質子通量的情況。

對太陽活躍年和太陽平靜年，在假定不同任務期時間1—365天中，我們分別比較了太陽質子通量分布與對數正態分布的關系。圖4顯示了太陽活躍年期間，假定任務時間為5天和100天時太陽質子通量的分布。圖中縱軸是對數坐標下的質子通量，橫軸是正態分布坐標下樣本中通量小于給定值的比例，也就是質子通量與（i×100）/ （N+1）的關系，i是某一樣本通量值在所有樣本從最小到最大排序中的排序數，N是樣本總數。如果太陽質子的通量分布符合對數正態分布，圖中的數據將分布在一條直線上。從圖中可以看出，多數數據分布在一條直線上，即太陽質子通量分布基本符合對數正態分布，即其概率密度函數為式（2）：

圖3 在不同時間長度的任務期中太陽質子通量的發生概率Fig.3 The production probability of solar proton fluences with different space mission periods from 1 day to 365 days

式中，F為任務期的太陽質子通量，σ是通量對數的標準偏差，μ是通量對數的期望值。σ和μ的數值都可以從太陽質子事件的通量分布圖中擬合獲得。在航天任務中，大的太陽質子通量更為關注，因此我們只對50%概率以上的數據點進行了直線擬合。對于太陽活躍年和太陽平靜年中的其他任務時間段，太陽質子通量分布也基本都符合對數正態分布。對應不同任務期，σ和μ的擬合取值見表2。

圖4 太陽活躍年中任務時間為5天和50天時太陽質子通量的分布Fig.4 The distribution of solar proton fluences in solar active years

表2 太陽質子通量對數正態分布參數Table 2 The parameter values of lognormal distribution for solar proton fluences

利用該方法我們可以計算出在一定置信度下不同航天任務期的質子通量。圖5中給出了太陽活躍年和平靜年期間，任務期時間為50天時的質子通量概率分布。通過查找縱坐標對應曲線點的橫坐標，可以得到某種置信度下的質子通量。如，在太陽活躍年，當任務期時間為50天時，要求給出95%的置信度下的質子通量。通過查找圖5中對應任務期的質子通量概率曲線，可以得到縱坐標1%～95%對應的曲線點橫坐標為3.3×109p·cm-2，該通量值即為滿足要求的質子通量。

圖5 任務期時間為50天時的質子通量概率分布Fig.5 The probability distribution of solar proton fluences with 50 days space mission

5 結論

準確計算任務期內可能遭遇的太陽質子事件通量水平，可以為航天器的抗輻射加固設計和宇航員的出艙防護提供重要的依據。本文利用1967—2014年期間共48年的太陽高能粒子觀測數據，統計分析了太陽質子事件的發生規律。針對持續1年以下的航天任務，我們對太陽活躍年和太陽平靜年分別分析研究了在1—365天不同時間段中的太陽質子事件和大于10 MeV質子通量的發生概率。進而，通過對比太陽質子通量分布與對數正態分布的關系，得到太陽質子通量分布基本符合對數正態分布。以該關系為基礎，我們計算出了一定置信度下，1年以下航天任務期內的質子通量分布，初步建立起了中短期太陽質子事件通量模型。

本文初步建立的太陽質子事件通量模型已經能夠為航天器執行中短期任務提供太陽質子事件通量預測。不過，該模型只提供了＞10 MeV的太陽質子通量的分布，對于更高能量的質子通量分布并未涉及。在以后的工作中，將繼續完善模型。值得注意的是：模型適用的任務期時間為1年以下，對1年以上航天任務期的太陽質子事件通量預測則需參考太陽質子通量的長期預測模型。

（

）

［1］ King，J.H.Solar proton fluences for 1977-1983 space missions［J］.Journal of Spacecraft and Rockets，1974，11（6）：401-409.

［2］ Feynman J，Armstrong T P，Dao-Gibner L，et al.New interplanetary proton fluence model［J］.Journal of Spacecraft and Rockets，1990，27（4）：403-410.

［3］ Feynman J，Spitale G，Wang J，et al.Interplanetary proton fluence model：JPL 1991［J］.Journal of Geophysical Research-All Series，1993，98（A8）：13-281.

［4］ Xapsos M A，Summers G P，Shapiro P，et al.New techniques for predicting solar proton fluences for radiation effects applications［J］.Nuclear Science，IEEE Transactions on，1996，43（6）：2772-2777.

［5］ Xapsos M A，Summers G P，Burke E A.Extreme value analysis of solar energetic proton peak fluxes［J］.Solar Physics，1998，183（1）：157-164.

［6］ Xapsos M A，Summers G P，Burke E A.Probability model for peak fluxes of solar proton events［J］.Nuclear Science，IEEE Transactions on，1998，45（6）：2948-2953.

［7］ Xapsos M A，Summers G P，Barth J L，et al.Probability model for worst case solar proton event fluences［J］.Nuclear Science，IEEE Transactions on，1999，46（6）：1481-1485.

［8］ Xapsos M A，Summers G P，Barth J L，et al.Probability model for cumulative solar proton event fluences［J］.Nuclear Science，IEEE Transactions on，2000，47（3）：486-490.

［9］ Xapsos M A，Barth J L，Stassinopoulos E G，et al.Characterizing solar proton energy spectra for radiation effects applications［J］.Nuclear Science，IEEE Transactions on，2000，47（6）：2218-2223.

［10］ Nymmik R A.Averaged energy spectra of peak flux and fluence values in solar cosmic ray events［C］//International Cosmic Ray Conference.1993，3（7）：29.

［11］ Shea M A，SmartD F，AdamsJ H，et al.Toward a descriptive model of solar particles in the heliosphere［C］//Interplanetary Particle Environment.JPL Publication，1988，88-28.

［12］ ISO 15391.Space Environment（Natural and Artifical）Probabilistic model for fluences and peak fluxes of solar energetic particles：Part I Protons［S］.2004.

［13］ 李婷婷，師立勤，劉四清.太陽質子通量模型研究［J］.空間科學學報，2010，30（3）：205-210.

Li tingting，Shi liqin，Liu siqing.Research on the Solar Proton Fluence Model［J］，Chin.J.Space Sci.，2010，30（3）：205-210.（in Chinese）

［14］ NASA/Goddard Space Flight Center.Paths to Magnetic field，Plasma，Energetic particle data relevant to heliospheric studies and resident at Goddard's Space Physics Data Facility ［EB/OL］.2015.http：//omniweb.gsfc.nasa.gov/

［15］ National Centers for Environmental Information.GOES Space Environment Monitor-Data Access［EB/OL］.2015.http：// www.ngdc.noaa.gov/stp/satellite/goes/dataaccess.html

Study on Forecast of Solar Proton Event Fluence during Short and Medium Term Missions

CUI Yanmei，SHI Liqin，LIU Siqing
（National Space Science Center，Chinese Academy of Sciences，Beijing 100190，China）

Forecasting the solar proton event fluence is very important to guarantee the safety of spacecrafts and astronauts.Within various periods less than 1 year，the occurrence probabilities of the solar proton events and their fluences were analyzed for solar active years and solar quiet years，respectively.It is found that the probability distributions of solar proton event fluence basically satisfy the log-normal distributions.Based on this，solar proton event fluence can be forecasted with some given confidence levels for different duration space missions ranging from 1 day to 365 days.

solar proton events；solar proton fluences；space mission；forecast

P35

A

1674-5825（2015）06-0597-06

2015-02-13；

2015-09-25

崔延美（1980-），女，博士，副研究員，研究方向為空間環境預報。E-mail：ymcui@nssc.ac.cn