CN110231643B - 高能电子暴事件的预报方法、装置及一种存储介质和设备 - Google Patents

Abstract

本发明提供种高能电子暴事件的预报方法、装置及一种存储介质和设备，包括以下步骤：提取当天的高能电子日通量F_T及中能电子指数数据I_T0；根据高能电子日通量F_T及中能电子指数数据I_T0计算第二天的高能电子通量F_T+1；根据计算得到的第二天的高能电子通量F_T+1，预报高能电子暴事件是否发生，如果F_T+1大于等于高能电子暴发生的高能电子日通量的第一阈值，则预报第二天发生高能电子暴事件，否则预报第二天不发生高能电子暴事件。本发明所记载的技术方案提高对地球同步轨道的高能电子暴事件预报的准确率，准确预报地球同步轨道高能电子通量，特别是高能电子暴事件的发生时间，不仅为地球同步轨道业务卫星的防护提供数据基础，也为运行在该区域的其他卫星的在轨防护提供参考。

Description

高能电子暴事件的预报方法、装置及一种存储介质和设备

技术领域

本发明涉及高能电子环境预报领域，具体涉及高能电子暴事件的预报方法、装置及一种存储介质和设备。

背景技术

空间中，尤其是地球同步轨道中的高能电子能够穿过卫星外层防护材料，沉积到内部的绝缘介质中，如果高能电子通量很高，则会形成局部的高电势状态，又称为“深层充电”，轻者影响卫星的正常运行，严重时甚至可以产生永久损伤。因此人们非常重视高能电子环境预报。地球同步轨道不仅运行着众多业务卫星，而且高电子通量的状态经常发生，因此地球同步轨道中高能电子环境预报显得尤为重要。

能够影响卫星安全的高能电子能量通常很高，卫星防护领域通常用能量大于2MeV的高能电子来表征高能电子环境的状况，国内外也发展了多种模型来预报高能电子的通量。现有预报模型主要基于地磁活动指数的变化、行星际太阳风的状况等相关参数进行统计和数值预报。而基于物理原理的数值仿真方法并不针对地球同步轨道而且不够稳定。上述方法均针对的是形成能量较低电子被加速的可能性和规模，而由于相关数据和物理过程知识的缺乏，基于这些传统空间环境数据的高能电子暴预报存在一定缺陷，特别是无法对高能电子暴的发生时间进行准确的预报。由于高能电子的加速过程与地磁环境扰动的相关性非常复杂，导致现有预报模型在预报高能电子暴事件(日通量大于10⁸cm^-2sr^-1)的开始、发展过程和结束的时间方面存在着一定缺陷。例如预报的准确率过低等。

发明内容

针对现有技术中的高能电子环境预报方面存在的上述不足和缺陷，本发明提供一种高能电子暴事件预报方法，通过该方法提高对地球同步轨道的高能电子暴事件预报的准确率，准确预报地球同步轨道高能电子通量，特别是高能电子暴事件第二天是否发生，以及高能电子日通量大小，不仅为地球同步轨道的业务卫星的防护提供数据基础，同时，也可以为运行在该区域的其他卫星的在轨防护提供参考。

根据本发明的第一方面，本发明提供了一种高能电子暴事件预报方法，包括以下步骤：

提取当天的高能电子日通量F_T及中能电子指数数据I_T0；

根据所述高能电子日通量F_T及所述中能电子指数数据I_T0计算第二天的高能电子通量F_T+1；

根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否发生，如果F_T+1大于等于高能电子暴发生的高能电子日通量的第一阈值，则预报第二天发生高能电子暴事件，否则预报第二天不发生高能电子暴事件。

可选地，根据所述高能电子日通量F_T计算第二天的高能电子通量F_T+1还包括以下步骤：

将提取的当天的所述高能电子日通量F_T与所述第一阈值进行比较；

根据比较结果选择不同的公式计算第二天的所述高能电子通量F_T+1。

可选地，根据比较结果选择不同的公式计算第二天的所述高能电子通量F_T+1还包括以下步骤：

当提取的当天的所述高能电子日通量F_T小于所述第一阈值时，采用公式(1)计算所述F_T+1，当F_T大于等于所述第一阈值时，则采用公式(2)计算所述F_T+1：

其中，F_T小于所述第一阈值且大于高能电子日通量的第二阈值时，Δ＝0；F_T小于所述第二阈值时，Δ＝1.1×(1gF_T-7)；I_T0为零度方向上的中能电子通量指数数据；当I_T0大于零度方向上的中能电子通量指数数据的第三阈值时，I_T0等于所述第三阈值。

可选地，还包括以下步骤：

当F_T大于所述第一阈值时，根据公式(2)计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否结束，否则预报所述高能电子暴事件尚未结束。

可选地，当所述F_T大于所述第一阈值时，提取前一天的中能电子指数数据I_(T-1)0；比较前一天的所述中能电子指数数据I_(T-1)0和当天的所述中能电子指数数据I_T0，判断所述高能电子暴事件是否结束。

可选地，比较前一天的所述中能电子指数数据I_(T-1)0和当天的所述中能电子指数数据I_T0，判断所述高能电子暴事件是否结束还包括以下步骤：

如果所述中能电子指数数据由前一天的I_(T-1)0小于所述中能电子指数数据的第四阈值，上升到当天的I_T0大于中能电子通量指数数据的第五阈值，则判断所述高能电子暴事件结束。

根据本发明的第二方面，本发明提供了一种高能电子暴事件预报装置，包括：

数据提取模块，用于提取当天的高能电子日通量F_T及中能电子指数数据I_T0；

计算模块，与所述数据提取模块通信连接，用于根据所述高能电子日通量F_T及所述中能电子指数数据I_T0计算第二天的高能电子通量F_T+1；以及事件预报模块，与所述计算模块通信连接，用于根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否发生。

根据本发明的第三方面，本发明提供了一种存储有计算机程序的计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行本发明第一方面所述的高能电子暴事件的预报方法。

根据本发明的第四方面，本发明提供了一种设备，该设备包括：

一个或多个处理器；

存储介质，通过总线与一个或多个所述处理器通信连接，用于存储一个或多个程序；

当一个或多个所述程序被一个或多个所述处理器执行时，一个或多个所述处理器执行本发明第一方面所述的高能电子暴事件的预报方法。

如上所述，本发明的高能电子暴事件的预报方法、预报装置及计算机可读存储介质和设备具有如下技术效果：

本发明的方法，包括提取当天的高能电子日通量F_T及中能电子指数数据I_T0，并根据所述高能电子日通量F_T的大小选择不同的公式计算第二天的高能电子通量F_T+1；根据计算得到的第二天的所述高能电子通量F_T+1预报高能电子暴事件是否发生。上述方法根据提取的当天的高能电子日通量的大小选择第二天电子通量的不同的计算公式，提高了计算结果的准确率，例如本发明所述方法对高能电子事件预报的准确率达到96.4％。远远高于现有技术中的方法对高能电子暴事件的预报正确率。

如上所述，本发明的方法通过准确预报地球同步轨道高能电子通量，特别是高能电子暴事件的第二天是否发生，以及高能电子日通量大小，不仅为地球同步轨道业务卫星的防护提供数据基础，同时，由于高能电子通量的增加发生在整个外辐射带，高能电子暴的预报信息也可以为运行在该区域的其他卫星的在轨防护提供参考。

附图说明

通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，在附图中：

图1为本发明实施例一提供的高能电子暴事件的预报方法的流程图。

图2显示为本发明实施例五提供的高能电子暴事件预报装置的示意图。

图3显示为本发明实施例七提供的设备的示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

发明人的研究表明，高能电子暴发生之前，能量较低的电子(又称为“种子电子”)的通量会有显著的增加，而高能电子暴发展过程中，“种子电子”通量的改变也能预示高能电子通量的变化趋势。由此，发明人也发现了高能电子暴结束的先兆特征。将这些先兆特征进行综合，便可形成基于低轨卫星低能电子探测数据的地球同步轨道高能电子通量预报方法。

本发明中，定量化地总结：1)高能电子暴发生前低轨探测到的大于100KeV的电子(以下简称“中能电子”)的特征；2)高能电子暴过程中，中能电子反映出的地球同步轨道高能电子(以下简称“高能电子”)通量变化趋势；3)高能电子暴的结束的先兆特征。由此形成高能电子暴事件的预报方法。

实施例一

本实施例提供一种高能电子暴事件的预报方法，所述方法包括以下步骤：

提取当天的高能电子日通量F_T及中能电子指数数据I_T0；

根据所述高能电子日通量F_T及所述中能电子指数数据I_T0计算第二天的高能电子通量F_T+1；

根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否发生，如果F_T+1大于高能电子暴发生的高能电子日通量的第一阈值，则预报第二天发生高能电子暴事件，否则预报第二天不发生高能电子暴事件。

现参照图1，对所述方法进行详细描述。如图1所示，首先在开始该方法之后提取当天的高能电子日通量F_T及中能电子指数数据I_T0。国际标准时每日3:00左右可以得到前一日的中能电子和高能电子通量的标准值，也可以通过前一日的实测数据进行计算，本实施例中，参考NOAA(National Oceanic and Atmospheric Administration，国家海洋大气管理局)给出的标准值。在本实施例的优选实施例中，可以通过互联网下载高能电子数据文件和中能电子指数文件。

例如所述高能电子数据文件的下载地址为：

https://services.swpc.noaa.gov/text/daily-particle-indices.txt；下载得到的高能电子数据示例如表1所示。

表1高能电子探测数据文件示例

如表1所示，下载的高能电子数据共分9列，第1-3列为探测日期，第8列为能量大于2MeV的高能电子日通量。在本实施例中取第8列的高能电子日通量的数据。

中能电子指数文件的下载地址为：

http://satdat.ngdc.noaa.gov/sem/poes/data/belt_indices/noaa19/bi_N19_ XXXX.txt；下载得到的中能电子指数文件示例如下表2所示。

表2中能电子指数文件示例

如表2所示，下载的中能电子指数文件中，每13行为一天的NOAA卫星观测数据，数据内容包括8列，在本实施例中，提取中的第8列的外辐射带指数中的第一行，是E>100KeV的0度方向的中能电子指数，即I_T0。

在本实施例的优选实施例中，以2018年2月17日为例，在下载的高能电子探测数据文件和中能电子指数文件中提取2月17日当天的高能电子日通量F_T＝6.2ⅹ10⁷，I_T0＝10.299。比较F_T与第一阈值的大小。在本实施例的优选实施例中，该第一阈值取值为10⁸cm^{- 2}sr^-1。可见2月17日当天的F_T＝6.2ⅹ10⁷<10⁸cm^-2sr^-1，此时，选择如下公式(1)计算第二天2月18日的高能电子日通量：

上述公式中，当F_T小于所述第一阈值且大于高能电子日通量的第二阈值时，Δ＝0；F_T小于所述第二阈值时，Δ＝1.1×(lgF_T-7)；I_T0为零度方向上的中能电子通量指数数据；当I_T0大于零度方向上的中能电子通量指数数据的第三阈值时，I_T0等于所述第三阈值。

在本实施例的优选实施例中规定，上述第二阈值取值为10⁷cm^-2sr^-1，中能电子通量指数数据的第三阈值取值为70。

根据上述规定，将提取的2月17日当天的高能电子日通量F_T＝6.2ⅹ10⁷，I_T0＝10.299，以及Δ＝0带入上述公式(1)，得到第二天2月18日的高能电子日通量F_T+1＝1.02ⅹ10⁸。

第二天2月18日的高能电子日通量F_T+1的计算结果大于第一阈值取值为10⁸cm^-2sr^-1，因此可预报2月18日会发生高能电子暴事件。

在本实施例的另一优选实施例中，以2018年2月18日为例，在下载的高能电子探测数据文件和中能电子指数文件中提取2月18日当天的高能电子日通量F_T＝3.2ⅹ10⁸cm^-2sr^-1，I_T0＝42.38。比较F_T与第一阈值的大小。在本优选实施例中，该第一阈值取值同样为10⁸cm^-2sr^-1。可见2月18日当天的F_T＝＝3.2ⅹ10⁸cm^-2sr^-1>10⁸cm^-2sr^-1，由该提取结果可知2月18日已经在发生高能电子暴事件。并且此时，选择如下公式(2)计算第二天2月19日的高能电子日通量：

上述公式中，I_T0为零度方向上的中能电子通量指数数据；当I_T0大于零度方向上的中能电子通量指数数据的第三阈值时，I_T0等于所述第三阈值。

在本实施例的优选实施例中规定，中能电子通量指数数据的第三阈值取值为70。

根据上述规定，将提取的2月18日当天的高能电子日通量F_T＝3.2ⅹ10⁸cm^-2sr^-1，I_T0＝42.38，以及Δ＝0带入上述公式(2)，得到第二天2月19日的高能电子日通量F_T+1＝1.51ⅹ10⁹cm^-2sr^-1。

第二天2月19日的高能电子日通量F_T+1的计算结果大于第一阈值10⁸cm^-2sr^-1，因此可预报2月18日发生的高能电子暴事件将会在2月19日持续并发展。

在本实施例的优选实施例中，由于所述方法的目的在于预报高能电子暴事件是否发生，当高能电子日通量大于等于所述第一阈值时，均发生高能电子暴事件，因此，当计算结果大于等于所述第一阈值时，可直接将所述计算结果记为所述第一阈值，据此便可预报高能电子暴事件将会发生或者将会持续发展。

提取的上述高能电子数据文件中显示，2月18日的高能电子日通量为3.2ⅹ10⁸cm^{- 2}sr^-1；2月19日的高能电子日通量为：4.7ⅹ10⁸cm^-2sr^-1，即2月18日发生高能电子暴事件，并且在2月19日持续发生高能电子暴事件。可见，上述方法的预报结果与实际观测的结果相同，准确率高。

如果由上述公式(1)得到的第二天的高能电子日通量的计算结果小于上述第一阈值，则可以预报第二天不发生高能电子暴事件。

实施例二

本实施例同样提供一种高能电子暴事件的预报方法。在本实施例中，例如选择2018年5月27日为例。在下载的高能电子探测数据文件和中能电子指数文件中提取5月27日当天的高能电子日通量F_T＝1.2ⅹ10⁸cm^-2sr^-1，中能电子通量指数I_T0＝0.2404。5月27日的F_T大于所述第一阈值的10⁸cm^-2sr^-1，因此可知在5月27日当天发生高能电子暴事件。由此选择上述公式(2)计算第二天(即5月28日)的高能电子通量F_T+1。将提取的5月27日的上述F_T＝1.2ⅹ10⁸cm^-2sr^-1及I_T0＝0.2404带入公式(2)，求得第二天的F_T+1＝7.92ⅹ10⁷cm^-2sr^-1。该F_T+1小于所述第一阈值，因此，预报高能电子暴事件在5月28日结束。

如上所述，如果由上述公式(2)得到的第二天的高能电子日通量的计算结果小于上述第一阈值，则可预报当天发生的高能电子暴事件将在第二天结束。

上述实施例的方法根据提取的当天的高能电子日通量的大小选择第二天高能电子日通量的不同的计算公式，提高了预报结果的准确率。

实施例三

本实施例同样提供一种高能电子暴事件的预报方法，该方法与实施例一提供的高能电子暴事件的预报方法的不同之处在于，本实施例提供了另外一种预报高能电子暴事件结束的方法，所述方法包括以下步骤：

当提取的当天的高能电子日通量F_T大于上述第一阈值时，提取前一天的中能电子指数数据I_(T-1)0；

比较前一天的所述中能电子指数数据I_(T-1)0和当天的所述中能电子指数数据I_T0，判断所述高能电子暴事件是否结束。如果所述中能电子指数数据由前一天的I_(T-1)0小于所述中能电子指数数据的第四阈值，上升到当天的I_T0大于中能电子通量指数数据的第五阈值，则判断所述高能电子暴事件结束。

在本实施例的优选实施例中规定，上述第四阈值取值为1.5，第五阈值的取值为6.1。

实施例四

本实施例同样提供一种高能电子暴事件的预报方法，在本实施例中，对所述方法的预报结果的误差及上述方法的正确率进行计算。

根据实施例一和实施例二所述的高能电子暴事件的预报方法，对2018年1月1日至6月30日的181天的高能电子暴事件进行了预报检验，统计结果如下表3所示：

表3 2018年1月1日至6月30日高能电子日通量统计结果

由上表3的统计可知，2018年上半年共有57个高能电子暴日，对这些事件进行统计预报精度统计，预报和实测的对数通量的平均相对误差用下式表示：

其中，F_P为计算得到的高能电子日通量，F_O为观测到的高能电子日通量。

经上述公式(3)计算，2018年1月1日-6月30日的对数通量的平均误差为2.88％。

表3的高能电子暴发生预报结果如下表4所示：

表4 2018年1月1日至6月30日高能电子暴统计结果

	预报有	预报无
			实测有	54	2
实测无	1	124

由上述表4的统计结果可计算所述方法的预报准确率R＝54/(54+2)＝96.4％。

由此可见，所述方法的预报准确率较高，通过准确预报地球同步轨道高能电子通量，特别是高能电子暴事件的发生时间和大小，不仅为地球同步轨道业务卫星的防护提供数据基础，同时，由于高能电子通量的增加发生在整个外辐射带，高能电子暴的预报信息也可以为运行在该区域的其他卫星的在轨防护提供参考。

实施例五

本实施例提供了一种高能电子暴事件预报装置，如图2所示，该装置包括：

数据提取模块，用于提取当天的高能电子日通量F_T及中能电子指数数据I_T0；

计算模块，与所述数据提取模块通信连接，用于根据所述高能电子日通量F_T及所述中能电子指数数据I_T0计算第二天的高能电子通量F_T+1；以及

事件预报模块，与所述计算模块通信连接，用于根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否发生。

在本实施例的一优选实施例中，所述事件预报模块中存储有高能电子暴发生的高能电子日通量的第一阈值，所述事件预报模块首先比较计算得到的第二天的所述高能电子通量F_T+1，如果F_T+1大于等于高能电子暴发生的高能电子日通量的第一阈值，则预报第二天发生高能电子暴事件，否则预报第二天不发生高能电子暴事件。

实施例六

本实施例提供一种存储有计算机程序的计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行上述实施例一至实施例三所述的高能电子暴事件的预报方法。

实施例七

本实施例提供一种设备，如图3所示，该设备包括处理器，及通过总线与处理器通信连接的存储介质。图中仅示出一个处理器，但是应该理解的是，该设备可以包括多个处理器。所述存储介质用于存储一个或者多个程序。当一个或多个程序被处理器执行时，处理器执行上述实施例一到实施例三所述的高能电子暴事件的预报方法。

综上，本发明的高能电子暴事件的预报方法、预报装置及计算机可读存储介质和设备具有如下技术效果：

本发明的方法，包括提取当天的高能电子日通量F_T及中能电子指数数据I_T0，并根据所述高能电子日通量F_T的大小选择不同的公式计算第二天的高能电子通量F_T+1；根据计算得到的第二天的所述高能电子通量F_T+1预报高能电子暴事件是否发生。上述方法根据提取的当天的高能电子日通量的大小选择第二天电子通量的不同的计算公式，提高了计算结果的准确率，例如本发明所述方法对高能电子事件预报的准确率达到96.4％。远远高于现有技术中的方法对高能电子暴事件的预报正确率。

如上所述，并发明的方法通过准确预报地球同步轨道高能电子通量，特别是高能电子暴事件第二天是否发生，以及高能电子日通量的大小，不仅为地球同步轨道业务卫星的防护提供数据基础，同时，由于高能电子通量的增加发生在整个外辐射带，高能电子暴的预报信息也可以为运行在该区域的其他卫星的在轨防护提供参考。

上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明，本领域技术人员可以在不脱离本发明的精神和范围的情况下作出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (8)

1.一种高能电子暴事件的预报方法，包括以下步骤：

提取当天的高能电子日通量F_T及中能电子指数数据I_T；

将提取的当天的所述高能电子日通量F_T与高能电子暴发生的高能电子日通量的第一阈值进行比较；

根据比较结果选择不同的公式计算第二天的所述高能电子通量F_T+1：当提取的当天的所述高能电子日通量F_T小于所述第一阈值时，采用公式(1)计算所述F_T+1，当F_T大于等于所述第一阈值时，则采用公式(2)计算所述F_T+1：

根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否发生，如果F_T+1大于等于所述第一阈值，则预报第二天发生高能电子暴事件，否则预报第二天不发生高能电子暴事件；

其中，F_T小于所述第一阈值且大于高能电子日通量的第二阈值时，Δ＝0；F_T小于所述第二阈值时，Δ＝1.1×(lgF_T-7)；I_T0为当天的零度方向上的中能电子通量指数数据；当I_T0大于零度方向上的中能电子通量指数数据的第三阈值时，I_T0等于所述第三阈值。

2.根据权利要求1所述的预报方法，其特征在于，还包括以下步骤：

当F_T大于所述第一阈值时，根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否结束。

3.根据权利要求2所述的预报方法，其特征在于，如果第二天的所述高能电子通量F_T+1小于所述第一阈值，则预报所述高能电子暴事件结束，否则预报所述高能电子暴事件尚未结束。

4.根据权利要求1所述的预报方法，其特征在于，还包括以下步骤：

当所述F_T大于所述第一阈值时，提取前一天的零度方向上的中能电子指数数据I_(T-1)0；

比较前一天的所述零度方向上的中能电子指数数据I_(T-1)0和当天的所述零度方向上的中能电子指数数据I_T0，判断所述高能电子暴事件是否结束。

5.根据权利要求4所述的预报方法，其特征在于，比较前一天的所述零度方向上的中能电子指数数据I_(T-1)0和当天的所述零度方向上的中能电子指数数据I_T0，判断所述高能电子暴事件是否结束还包括以下步骤：

如果所述零度方向上的中能电子指数数据由前一天的I_(T-1)0小于所述中能电子指数数据的第四阈值，上升到当天的I_T0大于中能电子通量指数数据的第五阈值，则判断所述高能电子暴事件结束。

6.一种高能电子暴事件预报装置，其特征在于，包括：

数据提取模块，用于提取当天的高能电子日通量F_T及中能电子指数数据I_T；

计算模块，与所述数据提取模块通信连接，用于将提取的当天的所述高能电子日通量F_T与高能电子暴发生的高能电子日通量的第一阈值进行比较，并根据比较结果选择不同的公式计算第二天的所述高能电子通量F_T+1；以及

事件预报模块，与所述计算模块通信连接，用于根据计算得到的第二天的所述高能电子通量F_T+1，预报高能电子暴事件是否发生；

其中，当提取的当天的所述高能电子日通量F_T小于所述第一阈值时，采用公式(1)计算所述F_T+1，当F_T大于等于所述第一阈值时，则采用公式(2)计算所述F_T+1：

F_T小于所述第一阈值且大于高能电子日通量的第二阈值时，Δ＝0；F_T小于所述第二阈值时，Δ＝1.1×(lgF_T-7)；I_T0为当天的零度方向上的中能电子通量指数数据；当I_T0大于零度方向上的中能电子通量指数数据的第三阈值时，I_T0等于所述第三阈值。

7.一种存储有计算机程序的计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行权利要求1-5中任一项所述的高能电子暴事件的预报方法。

8.一种用于高能电子暴事件的预报的设备，其特征在于，包括：

一个或多个处理器；

存储介质，通过总线与一个或多个所述处理器通信连接，用于存储一个或多个程序；

当一个或多个所述程序被一个或多个所述处理器执行时，一个或多个所述处理器执行权利要求1-5中任一项所述的高能电子暴事件的预报方法。

Images