一种太阳质子事件峰值通量预测方法与流程

1.本发明涉及太阳质子事件技术领域，具体为一种太阳质子事件峰值通量预测方法。


背景技术：

2.太阳质子事件是一种伴随太阳爆发出现的太阳发出的高能质子的急剧增加的现象，国际上通常以一次事件中能量大于10mev的高能质子的通量的最大值来确定事件的级别，而且，通过对卫星异常的统计，质子事件的影响与质子事件的峰值通量密切相关。
3.太阳质子事件的峰值通量不仅与太阳爆发事件的特征密切相关，与爆发点在日面上的位置也有联系，因此，通过观测太阳爆发活动的特征能够预测质子事件的有无，但无法准确预测质子事件的峰值通量。


技术实现要素：

4.本发明的目的在于克服现有技术的不足，提供一种太阳质子事件峰值通量预测方法。
5.本发明的目的是通过以下技术方案来实现的：
6.一种太阳质子事件峰值通量预测方法，包括以下步骤：
7.s1，实时获取太阳高能质子数据；
8.s2，从实时太阳高能质子数据中分别选取最近多组起始能量大于10mev以及起始能量大于50mev的高能质子数据，所述起始能量大于10mev的高能质子数据即为f10，所述起始能量大于50mev的高能质子数据即为f50，根据能谱因子计算公式得到能谱因子ri；
9.s3，根据能谱因子ri与f10的关系式lg(ri)＝c d*f10，利用最小二乘法得到常数c、d，根据多组f10的平均值以及lg(ri)的平均值利用相关系数计算公式得到相关系数rr；
10.s4，根据相关系数rr的计算结果，如果rr大于预设的阈值，则根据质子事件最大值的预测值公式f10
max
＝(lg(r
min
)-c)/d，式中，r
min
的值取决于当年太阳活动水平，为固定值，以此得到质子事件最大值的预测值f10
max
，如果rr小于等于预设的阈值，则返回步骤s1重新获取太阳高能质子数据；
11.s5，重复步骤s1-s4，获得多个质子事件最大值的预测值，当连续x个质子事件最大值的预测值之间的相对误差小于p时，计算这x个质子事件最大值的预测值的平均值，作为此次质子事件的峰值通量的预测值。
12.进一步的，所述步骤s2中，获取最近10组f10以及f50，分别记为f10i以及f50i，下角标i为提取通量的序数，分别为t到t-9，根据能谱因子计算公式ri＝f50i/f10i得到10个能谱因子ri。
13.进一步的，所述步骤s3中，相关系数rr的计算公式为：式中，x和分别表示f10i及其平均值，y和则分别表示lg(ri)及其平均值。
14.进一步的，所述步骤s4中，相关系数rr的预设阈值为0.9，如果rr大于0.9，则根据质子事件最大值的预测值公式计算质子事件最大值的预测值f10
max
，如果rr小于等于0.9，则返回步骤s1重新获取太阳高能质子数据。
15.进一步的，所述步骤s4中，将太阳活动水平分为太阳活动极大期、太阳活动中期、太阳活动极小期，将太阳活动极大年所在年份设为e，所述太阳活动极大期所在年份包括e、e 1、e-1，所述太阳活动中期所在年份包括e 3、e 2、e-3、e-2，其余年份为太阳活动极小期。
16.进一步的，所述太阳活动极大期，r
min
＝-1.8，所述太阳活动中期，r
min
＝-1.45，所述太阳活动极小期，r
min
＝-1.0。
17.进一步的，所述步骤s5中，当连续5个质子事件最大值的预测值之间的相对误差小于10％时，计算这5个质子事件最大值的预测值的平均值，将这个平均值作为此次质子事件的峰值通量的预测值。
18..一种数据处理装置，所述装置包括：
19.数据获取模块：用于获取太阳高能质子数据；
20.数据计算模块；用于根据太阳高能质子数据计算质子事件的峰值通量的预测值。
21.一种计算机存储设备，包括一个或多个存储器，其上存储有可执行程序，一个或多个处理器，用于执行所述存储器中的所述可执行程序，以实现所述一种太阳质子事件峰值通量预测方法的步骤。
22.本发明的有益效果是：
23.1、本发明专利中，通过一系列可实施步骤能有效确定太阳质子事件的峰值通量的预测值。
24.2、本发明专利的峰值通量计算流程启动的前提是质子事件的发生，因此，有助于避免误报现象发生。
附图说明
25.图1为本发明流程图；
26.图2为本发明专利的实施例流程判断图；
27.图3为本发明专利的能谱指数与f10的相关特征示意图；
28.图4为本发明专利实施例中预测的质子事件峰值(对数)与实测值(对数)的相关性示意图。
具体实施方式
29.下面结合本发明的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施。
30.在本发明的描述中，需要理解的是，术语“逆时针”、“顺时针”“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
31.一种太阳质子事件峰值通量预测方法，包括以下步骤：
32.s1，实时获取太阳高能质子数据，太阳高能质子数据是通过美国空间天气预报中心下载获取的；
33.s2，从实时太阳高能质子数据中分别选取最近多组起始能量大于10mev以及起始能量大于50mev的高能质子数据，所述起始能量大于10mev的高能质子数据即为f10，所述起始能量大于50mev的高能质子数据即为f50，根据能谱因子计算公式得到能谱因子ri；
34.s3，根据能谱因子ri与f10的关系式lg(ri)＝c d*f10，利用最小二乘法得到常数c、d，根据多组f10的平均值以及lg(ri)的平均值利用相关系数计算公式得到相关系数rr；
35.s4，根据相关系数rr的计算结果，如果rr大于预设的阈值，则根据质子事件最大值的预测值公式f10
max
＝(lg(r
min
)-c)/d，式中，r
min
的值取决于当年太阳活动水平，为固定值，以此得到质子事件最大值的预测值f10
max
，如果rr小于等于预设的阈值，则返回步骤s1重新获取太阳高能质子数据；
36.s5，重复步骤s1-s4，获得多个质子事件最大值的预测值，当连续x个质子事件最大值的预测值之间的相对误差小于p时，计算这x个质子事件最大值的预测值的平均值，作为此次质子事件的峰值通量的预测值。
37.进一步的，所述步骤s2中，获取最近10组f10以及f50，分别记为f10i以及f50i，下角标i为提取通量的序数，分别为t到t-9，根据能谱因子计算公式ri＝f50i/f10i得到10个能谱因子ri。
38.进一步的，所述步骤s3中，相关系数rr的计算公式为：式中，x和分别表示f10i及其平均值，y和则分别表示lg(ri)及其平均值。
39.进一步的，所述步骤s4中，相关系数rr的预设阈值为0.9，如果rr大于0.9，则根据质子事件最大值的预测值公式计算质子事件最大值的预测值f10
max
，如果rr小于等于0.9，则返回步骤s1重新获取太阳高能质子数据。
40.进一步的，所述步骤s4中，将太阳活动水平分为太阳活动极大期、太阳活动中期、太阳活动极小期，将太阳活动极大年所在年份设为e，所述太阳活动极大期所在年份包括e、e 1、e-1，所述太阳活动中期所在年份包括e 3、e 2、e-3、e-2，其余年份为太阳活动极小期。
41.进一步的，所述太阳活动极大期，r
min
＝-1.8，所述太阳活动中期，r
min
＝-1.45，所述太阳活动极小期，r
min
＝-1.0。
42.进一步的，所述步骤s5中，当连续5个质子事件最大值的预测值之间的相对误差小于10％时，计算这5个质子事件最大值的预测值的平均值，将这个平均值作为此次质子事件的峰值通量的预测值。
43..一种数据处理装置，所述装置包括：
44.数据获取模块：用于获取太阳高能质子数据；
45.数据计算模块；用于根据太阳高能质子数据计算质子事件的峰值通量的预测值
46.一种计算机存储设备，包括一个或多个存储器，其上存储有可执行程序，一个或多个处理器，用于执行所述存储器中的所述可执行程序，以实现所述一种太阳质子事件峰值通量预测方法的步骤。
47.具体的，通过对第23、24太阳活动周的质子事件的通量和能谱的变化过程进行了
统计分析，发现了太阳质子事件的能谱在事件进行过程中呈现出规律性的“软化”特征，即能量高的质子占比逐渐减小，直到达到质子事件峰值，此时，高能质子能谱达到“最软”。而且能谱的变化也有规律可循，能谱指数多数呈指数下降。因此，通过分析质子事件初期高能质子能谱及其发展趋势，能够实现对质子事件峰值的预测，本实施例计算了发生在2003年10月28日的一次质子事件的预报过程，我们以50分钟的滑动窗口(10个数据点)连续计算了参数c、d和rr，根据设定的rr的阈值，计算了本次质子事件的峰值通量，如表1所示，表中峰值通量的单位是cm-2
s-1
sr-1
。
48.表1 2003年10月28日部分时间段参数计算结果
49.时间utcdrr峰值通量12:35-9.57e-02-5.88e-03-9.50e-01230.483212:40-1.23e-01-4.78e-03-9.52e-01277.812912:45-1.72e-01-3.38e-03-9.23e-01378.469616:40-3.04e-01-1.13e-04-9.17e-0110162.116:45-3.02e-01-1.13e-04-9.26e-0110140.1516:50-2.94e-01-1.15e-04-9.26e-0110043.6616:55-2.95e-01-1.15e-04-9.18e-0110065.1917:25-4.98e-01-4.90e-05-9.62e-0119416.5717:30-4.84e-01-5.25e-05-9.83e-0118396.5317:35-4.64e-01-5.79e-05-9.57e-0117011.1317:40-4.77e-01-5.51e-05-9.26e-0117662.5117:45-5.07e-01-5.34e-05-9.03e-0117641.0717:50-4.99e-01-5.64e-05-9.49e-0116878.6117:55-4.96e-01-5.74e-05-9.44e-0116608.3218:00-5.09e-01-5.52e-05-9.62e-0117046.5818:05-5.22e-01-5.26e-05-9.44e-0117655.05
50.本次质子事件起始时间为12:10ut，12:35ut开始的3个数据点满足峰值通量计算条件，但其后相关系数rr不满足条件，因此不能作为峰值通量的预测结果；16:40-15:55ut尽管计算的峰值通量稳定，但没有达到稳定条件要求的5组数据；17:25-17:45ut，开始，预测的5个峰值通量平均值为：18925.56cm-2
s-1
sr-1
，每个预测值的相对误差小于15％，因此，对此次质子事件的峰值预测为：18925cm-2
s-1
sr-1
。该次质子事件的实测峰值通量达到时间为10月29日6:15ut，数值为29500cm-2
s-1
sr-1
。数据显示，这是因为激波到加速造成的，持续时间很短，很难预测。本例中，预测到了质子事件峰值通量将超过10000cm-2
s-1
sr-1
，提前量超过10小时，预报时的e》10mev质子通量为2000cm-2
s-1
sr-1
左右；
51.针对第24太阳活动周(2009-2020)的全部质子事件，按照所形成太阳质子事件峰值通量的计算流程，逐个进行计算和统计，得到的结果如表2所示。
52.表2 第24太阳活动周质子事件预警方案应用效果
[0053][0054][0055]
可以看到，多数质子事件(23/34)的峰值通量都可以通过本专利提供的方法预测太阳质子事件的峰值通量，计算失败的事例中(11个)，只有1个峰值通量仅为136cm-2
s-1
sr-1
，预测的太阳质子事件峰值与实测峰值正相关，如图3所示，二者取对数后的相关系数为0.88，6次峰值通量超过500cm-2
s-1
sr-1
的质子事件的峰值通量都实现了较好的预测。
[0056]
以上所述仅是本发明的优选实施方式，应当理解所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明并非局限于
本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。