CN-121995543-A - CME-induced geomagnetic storm disaster grade division method
Abstract
The invention discloses a CME-induced geomagnetic storm disaster grade division method, which relates to the technical field of space weather disaster monitoring and comprises the specific steps of data collaborative acquisition, data standardization pretreatment, parameter fusion inversion, trigger source positioning and grade prejudgment correction; according to the invention, multi-source satellite data are synchronously acquired through the solar dynamic astronomical station, the solar earth relation astronomical station and the Pake solar detector, and through a standardized preprocessing flow, the comprehensive parameter inversion formula of the CME is fused by utilizing the multi-source satellite data, so that the comprehensive inversion of multi-dimensional key parameters such as the CME burst time, the source region position, the three-dimensional propagation direction, the initial speed, the acceleration, the quality, the magnetic flux and the like is realized, and the innovation point is directly aimed at the defects of single data source and incomplete inversion parameters in the traditional method, and the precision and the reliability of the CME parameter inversion are improved through the multi-source data fusion and mutual verification mechanism, so that the data support is provided for the positioning and disaster grade division of the follow-up geomagnetic storm triggering source.
Inventors
- BI YI
- DENG LINHUA
- WAN MIAO
- LI YING
- XU TINGTING
Assignees
- 中国科学院云南天文台
Dates
- Publication Date
- 20260508
- Application Date
- 20260106
Claims (9)
- 1. A CME-induced geomagnetic storm disaster grade division method is characterized by comprising the following specific steps: The data collaborative collection comprises the steps of synchronously collecting three types of satellite data through an atmospheric imaging assembly carried by a solar dynamics astronomical platform, a solar spherical layer imaging instrument and a corona instrument carried by a solar earth relation astronomical platform, and a solar wind particle detection instrument and a magnetic field detection instrument carried by a Pake solar detector; Data standardization preprocessing, namely sequentially carrying out preprocessing operations of denoising, unifying formats, regulating time references and adjusting sampling frequencies and aligning observation coordinates on three types of acquired satellite data to obtain standardized data; Performing parameter fusion inversion, namely performing CME multi-dimensional key parameter inversion by combining multi-source satellite data with a CME comprehensive parameter inversion formula based on standardized data, and mutually verifying and calibrating through three types of satellite data; Triggering source positioning, namely calculating an intersection judgment index by utilizing an intersection judgment formula of the geomagnetic storm triggering source by combining the multi-dimensional key parameters of the CME and the space geometrical relationship of the sun and the earth, and judging whether the CME is the geomagnetic storm triggering source according to the intersection judgment index; And (3) level pre-judging and correcting, namely constructing a CME parameter and geomagnetic storm level association data set based on the historical data, and judging the geomagnetic storm level through a geomagnetic storm level judging formula.
- 2. The method for grading geomagnetic storm disaster caused by CME according to claim 1, wherein three types of satellite data comprise multi-band extreme ultraviolet coronal imaging data, stereoscopic observation data and near-sun detection data in the data collaborative collection, the multi-band extreme ultraviolet coronal imaging data are collected through an atmosphere imaging assembly to cover a plurality of extreme ultraviolet wave bands of 131 a, 171 a, 211 a, 304 a and 335 a, stereoscopic observation data are collected through a solar ball layer imager and a coronal instrument to cover a solar center angle range of 0-150 degrees, and near-sun detection data are collected through a solar wind particle detector and a magnetic field detector to comprise parameters of solar wind plasma density, speed, temperature, magnetic field intensity and direction.
- 3. The CME-induced geomagnetic storm disaster grade classification method of claim 1, wherein in the parameter fusion inversion, a multi-source satellite data fusion CME comprehensive parameter inversion formula is: , wherein, For the CME to integrate the characteristic parameters, Based on the preprocessed SDO multiband extreme ultraviolet coronal imaging data, the CME brightness characteristic quantity is obtained by counting the difference value between the gray average value of the CME area and the background gray, Based on the preprocessed stereoscopic observation data, the propagation direction angle of the CME is obtained by analyzing the azimuth angle of the CME in a Japanese coordinate system, Is CME substance characteristic quantity, is obtained by integrating the observed values of solar wind plasma density and magnetic field intensity based on the preprocessed near-sun detection data, 、 、 Is a fusion weight, and And (3) determining through historical data calibration.
- 4. A CME-induced geomagnetic storm disaster grade classification method as set forth in claim 3, wherein said parametric fusion inversion is based on The method comprises the steps of disassembling inversion CME multi-dimensional key parameters from corresponding satellite data sources, wherein the CME multi-dimensional key parameters comprise CME burst time, burst source region longitude and latitude, three-dimensional propagation direction, initial speed, acceleration, quality and magnetic flux, the burst time and burst source region longitude and latitude are obtained based on multi-band extreme ultraviolet coronal imaging data inversion, the three-dimensional propagation direction and initial speed are obtained based on stereoscopic observation data inversion, and the quality and magnetic flux are obtained based on near-sun detection data inversion.
- 5. The CME-induced geomagnetic storm disaster grade classification method of claim 1, wherein in said trigger source positioning, a geomagnetic storm trigger source intersection decision formula is: , wherein, In order to meet the judgment index, For the CME to integrate the characteristic parameters, Is the CME propagation speed, is obtained based on the inversion of the stereo observation data, In order to fix the distance constant of the sun, And obtaining the included angle between the CME propagation track and the earth orbit based on inversion of the stereoscopic observation data.
- 6. The method for grading a geomagnetic storm disaster caused by a CME as set forth in claim 5, wherein said trigger source location is based on Judging whether the CME propagation track can form an intersection with the earth orbit, when When 1.0X10 -4 s -1 is reached, the CME propagation track and earth track have junction, and the CME is determined as triggering source of geomagnetic storm triggering risk, when When 1.0X10 -4 s -1 is reached, the CME propagation track and the earth orbit are not possible to be intersected, and the CME is judged to not cause the magnetotelluric storm, and the grading flow is directly ended.
- 7. The method for grading a geomagnetic storm disaster caused by a CME according to claim 1, wherein in the grading pre-judgment correction, based on the history data of CME events, the CME multidimensional parameter and the geomagnetic storm grade corresponding to each event are sorted, a CME parameter and geomagnetic storm grade association data set is constructed, and geomagnetic storm grade index is calculated by a geomagnetic storm grade judgment formula 。
- 8. The CME-induced geomagnetic storm disaster grade classification method of claim 7, wherein in said grade pre-judgment correction, a geomagnetic storm grade judgment formula is: , wherein, Is the geomagnetic storm level index, 、 For the rank weight, by the CME parameter and geomagnetic storm rank association data set, For the CME to integrate the characteristic parameters, Is the intersection judgment index.
- 9. The CME-induced geomagnetic storm disaster grade division method of claim 8, wherein in said grade pre-judgement correction, based on geomagnetic storm grade index Determining geomagnetic storm level In the case of G1, when In the case of G2 stage, when In the case of G3 stage, when In the case of G4 stage, when In this case, the G5 stage is used.
Description
CME-induced geomagnetic storm disaster grade division method Technical Field The invention relates to the technical field of space weather disaster monitoring, in particular to a geomagnetic storm disaster grade division method caused by CME. Background The geomagnetic storm is taken as a strong space weather phenomenon, influences the electromagnetic environment of the earth, satellite operation, power grid safety and even human health, and along with the continuous deep exploration of the space and the increase of the degree of dependence of human beings, the accurate prediction and evaluation of the disaster level of the geomagnetic storm become particularly important, but the triggering mechanism of the geomagnetic storm is complex and changeable and mainly comes from solar activities, especially coronal material projection (CME) events, so that great challenges are brought to the prediction of geomagnetic storm. However, the traditional magnetic storm prediction technology mainly depends on data of a single satellite or a ground observation station, the methods have obvious defects in the comprehensive performance and the data processing fineness of a data source, specifically, single satellite observation often only can acquire a certain dimension or local characteristic of a CME event, for example, the three-dimensional propagation direction and speed of the CME are difficult to accurately invert only through extreme ultraviolet imaging data, or key parameters such as the quality and magnetic flux of the CME cannot be accurately measured only by depending on stereo observation data, the limitation of the data source causes larger errors when the traditional method inverts the key parameters of the CME in multiple dimensions, the physical characteristics of the CME are difficult to comprehensively and accurately describe, and meanwhile, a collaborative acquisition and fusion inversion mechanism of multi-source satellite data is lacked, so that the traditional method cannot effectively integrate and analyze data from different satellite platforms and observation types when processing complex space weather events, and further the accuracy of geomagnetic storm trigger source positioning and disaster grade division are limited. Therefore, a geomagnetic storm disaster grade division method caused by CME is developed. Disclosure of Invention The invention aims to make up the defects of the prior art and provides a CME-induced geomagnetic storm disaster grade division method, which synchronously collects multi-source satellite data through a solar dynamic astronomical station, a solar-earth relation astronomical station and a Pake solar detector, utilizes the multi-source satellite data to fuse a CME comprehensive parameter inversion formula through a standardized preprocessing flow, the comprehensive inversion of multi-dimensional key parameters such as CME burst time, source region position, three-dimensional propagation direction, initial speed, acceleration, quality, magnetic flux and the like is realized, the innovation point is directly aimed at the defects of single data source and incomplete inversion parameters in the traditional method, the precision and reliability of CME parameter inversion are improved through a multi-source data fusion and mutual verification mechanism, and data support is provided for the positioning and disaster grade division of a follow-up geomagnetic storm triggering source. The invention provides a CME-induced geomagnetic storm disaster grade division method for solving the technical problems, which comprises the following specific steps: The data collaborative collection comprises the steps of synchronously collecting three types of satellite data through an atmospheric imaging assembly carried by a solar dynamics astronomical platform, a solar spherical layer imaging instrument and a corona instrument carried by a solar earth relation astronomical platform, and a solar wind particle detection instrument and a magnetic field detection instrument carried by a Pake solar detector; Data standardization preprocessing, namely sequentially carrying out preprocessing operations of denoising, unifying formats, regulating time references and adjusting sampling frequencies and aligning observation coordinates on three types of acquired satellite data to obtain standardized data; Performing parameter fusion inversion, namely performing CME multi-dimensional key parameter inversion by combining multi-source satellite data with a CME comprehensive parameter inversion formula based on standardized data, and mutually verifying and calibrating through three types of satellite data; Triggering source positioning, namely calculating an intersection judgment index by utilizing an intersection judgment formula of the geomagnetic storm triggering source by combining the multi-dimensional key parameters of the CME and the space geometrical relationship of the sun and the earth, and judging wh