CN109948118B - Correction method, system, medium and terminal of geomagnetic activity index based on MSISE00 atmospheric model - Google Patents

Abstract

The invention provides a method, a system, a medium and a terminal for correcting a geomagnetic activity index based on an MSISE00 atmospheric model, which comprise the following steps: acquiring polar region energy injection data; calculating an average energy injection rate in each three-hour interval based on the polar region energy injection data; calculating a corrected geomagnetic activity index every three hour intervals based on the average energy injection rate; and acquiring the corrected geomagnetic activity index corresponding to a natural day or a specified time interval based on the corrected geomagnetic activity index in every three hour intervals. The method, the system, the medium and the terminal for correcting the geomagnetic activity index based on the MSISE00 atmospheric model correct the geomagnetic activity index input by the MSISE00 model based on the statistical relationship and the change rule between polar region energy injection and the geomagnetic activity index, thereby obviously improving the calculation accuracy of the atmospheric model.

Description

Correction method, system and medium of geomagnetic activity index based on MSISE00 atmospheric model quality, terminal

technical field

The present invention relates to the technical field of geomagnetic activity index, in particular to a method, system, medium and terminal for correcting geomagnetic activity index based on the MSISE00 atmospheric model.

Background technique

The impact of space weather on the upper atmosphere is mainly determined by the magnitude of the energy injection in the Earth's polar regions that accompanies geomagnetic disturbances, and is statistically positively correlated with the strength of geomagnetic activity. the magnitude of change. The polar region energy injection is obtained through the analysis of the results of NOAA series satellites and low-energy particle detectors carried by European Space Agency (ESA) polar-orbiting meteorological satellites. The geomagnetic activity index is obtained comprehensively from the geomagnetic stations on the ground. As far as the changes in the upper atmosphere are concerned, the energy injection in the polar region is undoubtedly a more direct impact index.

Many analyzes and studies have found that when the atmospheric model evaluates the orbital atmospheric density, overestimation and underestimation usually occur. When the international standard Kp index is used to input the atmospheric model, the 1-day density integral error of the calculation result sometimes exceeds 1 times. Therefore, ensuring the accurate operation of atmospheric models has become an important task of space weather. But commonly used atmospheric models are formed on the basis of a large amount of data and research. Due to the early compilation, the geomagnetic activity index is used to characterize the influence of geomagnetic disturbance. However, the input parameters of the atmospheric model represented by MSISE00 are in the form of geomagnetic activity index Kp or Ap.

The latest research results prove that the energy injected into high latitude regions during geomagnetic storms can cause a global increase in atmospheric density. Therefore, using parameters related to energy injection parameters can effectively improve the calculation results of upper atmosphere models. That is to say, it is necessary to make reasonable corrections to the geomagnetic activity index according to the magnitude of energy injection, and then drive the atmospheric model to obtain a more reasonable atmospheric density result.

Contents of the invention

In view of the shortcomings of the prior art described above, the object of the present invention is to provide a correction method, system, medium, and terminal based on the geomagnetic activity index of the MSISE00 atmospheric model, based on the statistical relationship and change between the energy injection in the polar region and the geomagnetic activity index According to the regularity, the geomagnetic activity index input into the MSISE00 model is corrected, thereby significantly improving the calculation accuracy of the atmospheric model.

In order to achieve the above purpose and other related purposes, the present invention provides a method for correcting the geomagnetic activity index based on the MSISE00 atmospheric model, comprising the following steps: obtaining polar region energy injection data; calculating intervals every three hours based on the polar region energy injection data The average energy injection rate within; based on the average energy injection rate, the corrected geomagnetic activity index for every three-hour interval is calculated; based on the corrected geomagnetic activity index for every three-hour interval, a natural day or a specified time interval corresponding to The adjusted geomagnetic activity index of .

In an embodiment of the present invention, calculating the average energy injection rate in every three-hour interval based on the energy injection data in the polar region includes the following steps:

Divide the daily data detection time into 8 three-hour intervals;

计算每三小时区间内的平均能量注入率；其中，N表示该三小时区间内的探测时间点个数，hp_i和t_i分别为第i个探测时间点的能量注入和时间。According to the formula

Calculate the average energy injection rate in every three-hour interval; where, N represents the number of detection time points in the three-hour interval, and hp _i and t _i are the energy injection and time of the i-th detection time point, respectively.

In an embodiment of the present invention, the average energy injection rate of the current three-hour interval is input into the HpAp model to obtain the corrected geomagnetic activity index corresponding to the three-hour interval.

In an embodiment of the present invention, the corrected geomagnetic activity index of a natural day is the average value of the corrected geomagnetic activity index of eight three-hour intervals in the natural day; the corrected geomagnetic activity index corresponding to the specified time interval includes the current The corrected geomagnetic activity index for the three-hour period, the corrected geomagnetic activity index for the first three-hour period, the corrected geomagnetic activity index for the first two three-hour periods, and the corrected geomagnetic activity index for the first three three-hour periods Geomagnetic Activity Index, the average of the corrected Geomagnetic Activity Index for the previous 5-12 three-hour intervals, and the average of the corrected Geomagnetic Activity Index for the previous 13-20 three-hour intervals.

Correspondingly, the present invention provides a correction system based on the geomagnetic activity index of the MSISE00 atmospheric model, including an acquisition module, a first calculation module, a second calculation module and a correction module;

The acquiring module is used to acquire polar region energy injection data;

The first calculation module is used to calculate the average energy injection rate in every three-hour interval based on the energy injection data in the polar region;

The second calculation module is used to calculate the corrected geomagnetic activity index for every three-hour interval based on the average energy injection rate;

The correction module is configured to obtain a corrected geomagnetic activity index corresponding to a natural day or a specified time interval based on the corrected geomagnetic activity index in each three-hour interval.

In an embodiment of the present invention, the calculation of the average energy injection rate in every three-hour interval by the first calculation module based on the energy injection data in the polar region includes the following steps:

Divide the daily data detection time into 8 three-hour intervals;

计算每三小时区间内的平均能量注入率；其中，N表示该三小时区间内的探测时间点个数，hp_i和t_i分别为第i个探测时间点的能量注入和时间。According to the formula

Calculate the average energy injection rate in every three-hour interval; where, N represents the number of detection time points in the three-hour interval, and hp _i and t _i are the energy injection and time of the i-th detection time point, respectively.

In an embodiment of the present invention, the second calculation module inputs the average energy injection rate of the current three-hour interval into the HpAp model to obtain the corrected geomagnetic activity index corresponding to the three-hour interval.

In one embodiment of the present invention, in the correction module, the corrected geomagnetic activity index of a natural day is the average value of the corrected geomagnetic activity index of eight three-hour intervals in the natural day; the corrected geomagnetic activity index corresponding to the specified time interval The geomagnetic activity index includes the corrected geomagnetic activity index for the current three-hour interval, the corrected geomagnetic activity index for the previous three-hour interval, the corrected geomagnetic activity index for the first two three-hour intervals, and the corrected geomagnetic activity index for the first three three-hour intervals. The corrected geomagnetic activity index for the hourly interval, the average of the corrected geomagnetic activity index for the previous 5-12 three-hour intervals, and the average of the corrected geomagnetic activity index for the previous 13-20 three-hour intervals.

The present invention provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the above method for correcting the geomagnetic activity index based on the MSISE00 atmospheric model is realized.

The present invention provides a terminal, including: a processor and a memory;

The memory is used to store computer programs;

The processor is configured to execute the computer program stored in the memory, so that the terminal executes the above-mentioned method for correcting the geomagnetic activity index based on the MSISE00 atmospheric model.

As mentioned above, the correction method, system, medium and terminal of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention have the following beneficial effects:

(1) Correct the geomagnetic activity index input by the MSISE00 model based on the statistical relationship and change law between the energy injection in the polar region and the geomagnetic activity index;

(2) Compared with the standard geomagnetic activity index, the corrected geomagnetic activity index significantly improves the calculation accuracy of the atmospheric model and reduces the error;

(3) Contribute to further research on space weather.

Description of drawings

Fig. 1 shows the flow chart of the correction method of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention in an embodiment;

Fig. 2 shows the comparative schematic diagram of the calculation result of atmospheric density obtained respectively and the measured value based on the MSISE00 atmospheric model according to the standard and corrected geomagnetic activity index in an embodiment;

Fig. 3 shows the comparison schematic diagram of the calculation results of atmospheric density and measured values obtained respectively according to the geomagnetic activity index after standard and correction based on the MSISE00 atmospheric model in another embodiment;

Fig. 4 shows the structural representation of the correction system of the geomagnetic activity index based on the MSISE00 atmospheric model in an embodiment of the present invention;

FIG. 5 is a schematic structural diagram of a terminal of the present invention in an embodiment.

Component designation description

41 Get module

42 First computing module

43 Second computing module

44 Calibration module

51 processors

52 memory

detailed description

Embodiments of the present invention are described below through specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. The present invention can also be implemented or applied through other different specific implementation modes, and various modifications or changes can be made to the details in this specification based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that, in the case of no conflict, the following embodiments and features in the embodiments can be combined with each other.

It should be noted that the diagrams provided in the following embodiments are only schematically illustrating the basic ideas of the present invention, and only the components related to the present invention are shown in the diagrams rather than the number, shape and shape of the components in actual implementation. Dimensional drawing, the type, quantity and proportion of each component can be changed arbitrarily during actual implementation, and the component layout type may also be more complicated.

The correction method, system, medium and terminal of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention correct the geomagnetic activity index input by the MSISE00 model and reduce atmospheric density calculation error, thus significantly improving the calculation accuracy of the atmospheric model.

When dealing with the process of geomagnetic storms affecting the atmospheric density, geomagnetic storms are divided into two types: recurring coronal holes and coronal mass ejections (Coronal Mass Ejection, CME), and statistical analysis of the upper atmosphere density corresponding to geomagnetic storms Variety. The statistical analysis results show that:

(1) For geomagnetic storms caused by coronal holes, the MSISE00 model inputs the measured geomagnetic activity index Ap, and the error between the calculated atmospheric density and the satellite measured value is positive or negative, and the single-day average error is less than 10%;

(2) For the geomagnetic storm caused by CME, if the MSISE00 model inputs the measured Ap value, the calculated value of the model will be lower than the measured value.

Furthermore, the relevant statistical analysis of the high-latitude energy injection corresponding to the geomagnetic storms of the above two generation mechanisms is carried out, and the conclusions obtained are as follows:

(1) For geomagnetic storms caused by recurrent coronal holes, the energy injection is roughly proportional to the geomagnetic activity index Ap;

(2) Under the condition of the same geomagnetic activity index Ap value, the energy injection accompanying the geomagnetic storm caused by CME is significantly higher than that caused by recurrent coronal holes.

Therefore, the present invention uses the energy injection index to correct the geomagnetic activity index Ap, and uses the corrected geomagnetic activity index AP as an input parameter of MSISE00, thereby improving the calculation result of the atmospheric model. Because under the same geomagnetic activity index Ap, the energy injection associated with the geomagnetic storm caused by the CME is higher than that of the geomagnetic storm caused by the recurring coronal hole. The measured geomagnetic activity index Ap is high; since the higher the geomagnetic activity index Ap value calculated under the same conditions of the MSISE00 model, the greater the calculation result of the atmospheric density, so energy injection can effectively solve the problem of underestimation of the atmospheric density during geomagnetic storms caused by CME . Therefore, the key point of the correction method of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention is to obtain the corrected geomagnetic activity index Ap value, namely Ap*.

As shown in Figure 1, in one embodiment, the correction method of the geomagnetic activity index based on MSISE00 atmospheric model of the present invention comprises the following steps:

Step S1. Obtain energy injection data in the polar region.

Specifically, download from address ftp://ftp.swpc.noaa.gov/pub/lists/hpi/swpc_aurora_power_yyyymmdd.txt or address ftp://ftp.swpc.noaa.gov/pub/lists/hpi/power_yyyy.txt Polar region energy injection data. Among them, the first address can obtain the daily energy injection files of the last 30 days, and the second address can obtain the historical data of polar energy injection in units of years.

Among them, the file header and data content of the polar region energy injection input file are as follows:

Wherein, in the data content, each row is data of a time point, including 5 columns of data in total. The first column indicates the satellite and detection position, and the N/S in brackets indicate the data of the South/Arctic monitoring respectively; the second column indicates the time, expressed as DDDHHMM, DDD is the serial number of the date in the year, and HHMM indicates the hour and minute at that time; The third column indicates the energy injection rate in the polar region, in GW; the fourth and fifth columns indicate the activity level and normalization coefficient, respectively.

Step S2, calculating the average energy injection rate in every three-hour interval based on the energy injection data in the polar region.

Specifically, the daily data detection time is divided into eight three-hour intervals of 0-2, 3-5, 6-8, 9-11, 12-14, 15-17, 18-20, and 20-23. Calculate the average energy injection rate for each three-hour interval according to the following formula:

Among them, N represents the number of detection time points within the three-hour interval, and hp _i and t _i are the energy injection and time of the i-th detection time point, respectively.

Step S3 , calculating the corrected geomagnetic activity index ap* for each three-hour interval based on the average energy injection rate.

Specifically, the average energy injection rate of the current three-hour interval is input into the HpAp model to obtain the corrected geomagnetic activity index ap* corresponding to the three-hour interval. Among them, the HpAp model is an empirical formula obtained by counting the corresponding relationship between the 3-hour geomagnetic activity index ap and the energy injection hp of this period between 1996 and 2016:

hp=9.6+0.95ap- ^0.00053ap2

Through the above formula, the corrected geomagnetic activity index ap* value can be calculated from the measured hp value, that is, the average energy injection rate.

Step S4. Obtain the corrected geomagnetic activity index Ap* or ap* corresponding to a natural day or a specified time interval based on the corrected geomagnetic activity index ap* in every three-hour interval.

In the MSISE00 atmospheric model, there are two input methods for the geomagnetic activity index:

1) When the intensity of the local magnetic storm is small, if the average value of the geomagnetic activity index every 3 hours in a natural day is less than 80, the geomagnetic activity index of a natural day is used as input;

2) When a geomagnetic storm other than situation 1) occurs, that is, a strong geomagnetic storm, the geomagnetic activity index of each three-hour interval needs to be input sequentially.

Therefore, in the present invention, it is necessary to obtain a corrected geomagnetic activity index corresponding to a natural day or a specified time interval.

Specifically, the corrected geomagnetic activity index Ap* of a natural day (UT) is the average value of the corrected geomagnetic activity index ap* of eight three-hour intervals within the natural day. For the current three-hour interval, ap*(1) is the corrected geomagnetic activity index for this interval, ap*(2) is the corrected geomagnetic activity index for the previous three-hour interval, and ap*(3) is the previous two The corrected geomagnetic activity index for the three-hour interval, ap*(4) is the corrected geomagnetic activity index for the first three three-hour intervals, and ap*(5) is the corrected geomagnetic activity index for the first 5-12 three-hour intervals Average Geomagnetic Activity Index, ap*(6) is the average of the corrected Geomagnetic Activity Index for the previous 13-20 three-hour intervals.

The method for correcting the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention will be further explained below through specific examples.

Taking December 15, 2006 as an example, perform the following steps:

Step 1: First go to the download address to download the energy injection file in 2006, then open the header of the file, and intercept the energy injection data of the day; then process the file to keep only the time of day (namely HHMM item) and the third column of polar region energy Inject the item, and rewrite the time item, that is, t=hh+mm/60.

Step 2: Calculate the average energy injection rate of the eight 3-hour intervals of the day.

得到该区间的平均能量注入率为145.9GW。Taking the interval 0-2 as an example, there are altogether 18 detection time points in this interval, that is, N=18. Among them, t ₁ =0.033, t _N =2.947, according to the formula

The average energy injection rate in this interval is 145.9GW.

Step 3: Call the HpAp model, input the average energy injection rate obtained in the second step, and obtain the 3-hour ap* of 8 intervals, as shown in the following table.

Table 1. Polar region energy injection and ap* on December 5, 2006

period of time 0-2 3-5 6-8 9-11 12-14 15-17 18-20 21-23 Hpi(GW) 146 133.8 132.5 130.6 94.4 65.9 47.9 42.3 ap* 338.2 299.7 295.4 289.1 183.3 102.9 61.2 49.3

Step 4: Obtain the Ap* required by the MSISE00 atmospheric model of the day, that is, calculate the average value of the 8 ap* of the day, and get the Ap* of 120. In addition, if it is necessary to obtain the parameters of the specified time period separately, it can be obtained by processing the ap* of the corresponding time period. For example, the ap*(I) in the 9:00-11:00 period is: ap*(1), the ap* of the current three-hour interval; ap*(2), the ap* of the previous three-hour interval; ap* (3), the ap* of the first two three-hour intervals; ap* (4), the ap* of the first three three-hour intervals; ap* (5), the average of the ap* of the first 5-12 three-hour intervals ; ap*(6), the average value of ap* for the first 13-20 three-hour intervals.

The correction method of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention will be further verified through specific experiments below.

Specifically, based on the MSISE00 atmospheric model, using the standard geomagnetic activity index Ap and the geomagnetic activity index Ap* retrieved from the energy injection in the polar region, the 21 natural days (UT) with Kp≥7 in 2003-2006 The orbital atmospheric density was calculated and compared with the CHAMP satellite detections at the same time. The final result is: the calculation error using Ap* in the 21 sample days is improved to varying degrees compared with the standard error using the standard Ap value, and the effective rate is 100%. However, due to the different physical processes of geomagnetic storms, the maximum improvement can be Up to 15%, the minimum is only 4.2%.

Among them, in the geomagnetic storm on April 14, 2006, the Kp value reached 7 (ap=150) in two time intervals, and the corresponding ap* was calculated according to the energy injection in the polar region on that day, and the results are shown in the following table:

Table 2. Detection values of geomagnetic storm process on April 14, 2006

period of time 0-2 3-5 6-8 9-11 12-14 15-17 18-20 21-23 ap 27 80 150 150 48 48 27 27 ap* 68.2 176.6 216.3 276.8 130.1 106.3 87.2 40.3

Substituting the above results into the MSISE00 atmospheric model, the calculation results of atmospheric density obtained by two different methods are compared with the measured values. The results are shown in Figure 2, where the unit of atmospheric density is g/ ^cm3 . It can be seen that the average relative error of the calculation results using ap* is 21.03%, while the average relative error of the calculation results using ap is 35.82%.

In the geomagnetic storm on December 15, 2006, the Kp value reached 8 (ap=250) in two time intervals. The corresponding ap* was calculated according to the energy injection in the polar region of the day, and the results are as follows:

Table 3. Detection values of geomagnetic storm process on December 15, 2006

Substituting the above results into the MSISE00 atmospheric model, the calculation results of atmospheric density obtained by two different methods are compared with the measured values. The results are shown in Figure 3, where the unit of atmospheric density is g/ ^cm3 . It can be seen that the average relative error of the calculation results using ap* is 12.38%, while the average relative error of the calculation results using ap is 24.24%.

As shown in FIG. 4 , in one embodiment, the correction system of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention includes an acquisition module 41 , a first calculation module 42 , a second calculation module 43 and a correction module 44 connected in sequence.

The acquisition module 41 is used to acquire polar region energy injection data;

The first calculation module 42 is used to calculate the average energy injection rate in every three-hour interval based on the energy injection data in the polar region;

The second calculation module 43 is used to calculate the corrected geomagnetic activity index for every three hour interval based on the average energy injection rate;

The correction module 44 is configured to obtain a corrected geomagnetic activity index corresponding to a natural day or a specified time interval based on the corrected geomagnetic activity index in each three-hour interval.

Wherein, the structure and principle of the acquisition module 41, the first calculation module 42, the second calculation module 43 and the correction module 44 correspond to the steps in the correction method of the geomagnetic activity index based on the above-mentioned MSISE00 atmospheric model, so no longer repeat.

It should be noted that it should be understood that the division of each module of the above device is only a division of logical functions, and may be fully or partially integrated into one physical entity or physically separated during actual implementation. Moreover, these modules can be implemented in the form of calling software through processing elements, or can be implemented in the form of hardware, or some modules can be implemented in the form of calling software through processing elements, and some modules can be implemented in the form of hardware. For example, the x module can be a separate processing element, and can also be integrated in a chip of the above-mentioned device. In addition, the x module can also be stored in the memory of the above-mentioned device in the form of program code, and can be invoked by a certain processing element of the above-mentioned device to execute the function of the above-mentioned x module. The implementation of other modules is similar. All or part of these modules can be integrated together, and can also be implemented independently. The processing element mentioned here may be an integrated circuit with signal processing capabilities. In the implementation process, each step of the above method or each module above can be completed by an integrated logic circuit of hardware in the processor element or an instruction in the form of software. The above modules may be one or more integrated circuits configured to implement the above method, for example: one or more specific integrated circuits (Application Specific Integrated Circuit, referred to as ASIC), one or more microprocessors (Digital Singnal Processor, DSP for short), one or more Field Programmable Gate Arrays (Field Programmable Gate Array, FPGA for short), and the like. When one of the above modules is implemented in the form of a processing element scheduling program code, the processing element may be a general-purpose processor, such as a central processing unit (Central Processing Unit, CPU for short) or other processors that can call program codes. These modules can be integrated together and implemented in the form of a System-on-a-chip (SOC for short).

The computer-readable storage medium of the present invention stores a computer program, and when the program is executed by a processor, the above method for correcting the geomagnetic activity index based on the MSISE00 atmospheric model is realized. Wherein, the storage medium includes: ROM, RAM, magnetic disk, U disk, memory card or optical disk, etc., which can store program codes.

As shown in FIG. 5 , in an embodiment, the terminal of the present invention includes: a processor 51 and a memory 52 .

The memory 52 is used to store computer programs.

The memory 52 includes various media capable of storing program codes such as ROM, RAM, magnetic disk, U disk, memory card or optical disk.

The processor 51 is connected to the memory 52, and is used to execute the computer program stored in the memory 52, so that the terminal executes the above-mentioned correction method of the geomagnetic activity index based on the MSISE00 atmospheric model.

Preferably, the processor 51 can be a general-purpose processor, including a central processing unit (Central Processing Unit, referred to as CPU), a network processor (Network Processor, referred to as NP), etc.; it can also be a digital signal processor (Digital Signal Processor, referred to as DSP), Application Specific Integrated Circuit (ASIC for short), Field Programmable Gate Array (Field-Programmable Gate Array, FPGA for short) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

In summary, the correction method, system, medium, and terminal of the geomagnetic activity index based on the MSISE00 atmospheric model of the present invention are based on the statistical relationship and change law between the energy injection in the polar region and the geomagnetic activity index, and the geomagnetic activity index input by the MSISE00 model. Correction; Compared with the standard geomagnetic activity index, the corrected geomagnetic activity index significantly improves the calculation accuracy of the atmospheric model and reduces the error; it is helpful for further research on space weather. The invention effectively overcomes various shortcomings in the prior art and has high industrial application value.

The above-mentioned embodiments only illustrate the principles and effects of the present invention, but are not intended to limit the present invention. Anyone skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or changes made by those skilled in the art without departing from the spirit and technical ideas disclosed in the present invention should still be covered by the claims of the present invention.

Claims (10)

1. a correction method based on the geomagnetic activity index of MSISE00 atmospheric model, it is characterized in that, may further comprise the steps: Obtain polar region energy injection data; calculating the average energy injection rate in every three-hour interval based on the energy injection data in the polar region; calculating a corrected geomagnetic activity index for each three-hour interval based on the average energy injection rate; The corrected geomagnetic activity index corresponding to a natural day or a specified time interval is obtained based on the corrected geomagnetic activity index in each three-hour interval. 2. the correction method of the geomagnetic activity index based on MSISE00 atmospheric model according to claim 1, is characterized in that, Calculating the average energy injection rate in every three-hour interval based on the energy injection data in the polar region includes the following steps: Divide the daily data detection time into 8 three-hour intervals; According to the formula

Calculate the average energy injection rate in every three-hour interval; where, N represents the number of detection time points in the three-hour interval, and hp _i and t _i are the energy injection and time of the i-th detection time point, respectively. 3. the correcting method of the geomagnetic activity index based on MSISE00 atmospheric model according to claim 1, is characterized in that, the average energy injection rate of current three-hour interval is input into HpAp model to obtain the corrected geomagnetic field corresponding to this three-hour interval activity index. 4. the correction method of the geomagnetic activity index based on MSISE00 atmospheric model according to claim 1, is characterized in that, the corrected geomagnetic activity index of a natural day is the corrected geomagnetic activity of 8 three-hour intervals in this natural day The average value of the index; the corrected geomagnetic activity index corresponding to the specified time interval includes the corrected geomagnetic activity index for the current three-hour interval, the corrected geomagnetic activity index for the previous three-hour interval, and the corrected geomagnetic activity index for the first two three-hour intervals. Corrected Geomagnetic Activity Index, corrected Geomagnetic Activity Index for the first 3 three-hour intervals, average of the corrected Geomagnetic Activity Index for the first 5-12 three-hour intervals, and corrections for the first 13-20 three-hour intervals The average value of the subsequent geomagnetic activity index. 5. A correction system based on the geomagnetic activity index of the MSISE00 atmospheric model, comprising an acquisition module, a first calculation module, a second calculation module and a correction module; The acquiring module is used to acquire polar region energy injection data; The first calculation module is used to calculate the average energy injection rate in every three-hour interval based on the energy injection data in the polar region; The second calculation module is used to calculate the corrected geomagnetic activity index for every three-hour interval based on the average energy injection rate; The correction module is configured to obtain a corrected geomagnetic activity index corresponding to a natural day or a specified time interval based on the corrected geomagnetic activity index in each three-hour interval. 6. the correction system of the geomagnetic activity index based on MSISE00 atmospheric model according to claim 5, it is characterized in that, described first calculation module calculates the average energy injection rate in every three hours interval based on described polar region energy injection data Include the following steps: Divide the daily data detection time into 8 three-hour intervals; According to the formula

Calculate the average energy injection rate in every three-hour interval; where, N represents the number of detection time points in the three-hour interval, and hp _i and t _i are the energy injection and time of the i-th detection time point, respectively. 7. the correcting system of the geomagnetic activity index based on MSISE00 atmospheric model according to claim 5, is characterized in that, described second calculation module inputs the average energy injection rate of current three-hour interval into HpAp model to obtain this three-hour interval The corresponding corrected geomagnetic activity index. 8. the correction system based on the geomagnetic activity index of MSISE00 atmospheric model according to claim 5, is characterized in that, in described correction module, the geomagnetic activity index after the correction of a natural day is 8 three-hour intervals in this natural day The average value of the corrected geomagnetic activity index; the corrected geomagnetic activity index corresponding to the specified time interval includes the corrected geomagnetic activity index of the current three-hour interval, the corrected geomagnetic activity index of the previous three-hour interval, the previous two The corrected geomagnetic activity index for each three-hour interval, the corrected geomagnetic activity index for the first three three-hour intervals, the average of the corrected geomagnetic activity index for the first 5-12 three-hour intervals, and the first 13-20 The average of the corrected geomagnetic activity index for the three-hour interval. 9. A computer-readable storage medium, on which a computer program is stored, it is characterized in that, when the program is executed by a processor, it realizes the geomagnetic activity index based on the MSISE00 atmospheric model described in any one of claims 1 to 4 Calibration method. 10. A terminal, comprising: a processor and a memory; The memory is used to store computer programs; The processor is configured to execute the computer program stored in the memory, so that the terminal executes the method for correcting the geomagnetic activity index based on the MSISE00 atmospheric model according to any one of claims 1 to 4.

Images