CN-121978714-A - High-credibility ionosphere scintillation monitoring method and system combining multiple indexes

CN121978714ACN 121978714 ACN121978714 ACN 121978714ACN-121978714-A

Abstract

The invention discloses a combined multi-index high-credibility ionosphere scintillation monitoring system, which belongs to the field of ionosphere scintillation monitoring and comprises the steps of calculating ionosphere scintillation intensity indexes based on historical observation data of a reference station, calculating scintillation intensity threshold values in each satellite altitude interval according to the relative distribution condition of the ionosphere scintillation intensity indexes and satellite altitude angles and combining phase scintillation intensity indexes and amplitude scintillation intensity indexes of all satellites, and calculating ROTI in the ionosphere scintillation intensity indexes, Training DI and AATR to obtain random forest classifier, obtaining ionosphere scintillation intensity and scintillation intensity probability, calculating residual variation rate of observation values of all satellites based on historical observation data of reference station, and evaluating reliability of ionosphere scintillation intensity and probability. The invention judges the current ionosphere scintillation intensity by utilizing a random forest method, gives intensity probability along with the judging result and provides credibility information for users.

Inventors

Chu Ruitao
Xu Chaohan
YAO YIBIN
KONG JIAN
PENG WENJIE
ZHANG QI

Assignees

武汉大学
湖北珞珈实验室

Dates

Publication Date: 20260505
Application Date: 20251229

Claims (10)

1. A combined multi-index high-confidence ionosphere scintillation monitoring method, comprising: based on the obtained historical observation data of the reference station, calculating to obtain ionosphere scintillation intensity indexes of all satellites, and simultaneously obtaining phase scintillation intensity indexes and amplitude scintillation intensity indexes of the satellites obtained by all special ionosphere scintillation receivers; calculating a scintillation intensity threshold value in each satellite altitude angle section according to the relative distribution condition of ionosphere scintillation intensity indexes of all satellites and satellite altitude angles and combining the phase scintillation intensity indexes and the amplitude scintillation intensity indexes of all satellites; based on the scintillation intensity threshold value in each satellite altitude interval, for the ROTI in the ionosphere scintillation intensity index, Training DI and AATR respectively to obtain respective random forest classifiers, and finally summarizing output results of all random forest classifiers to obtain ionosphere scintillation intensity and scintillation intensity probability; Based on the obtained historical observation data of the reference station, calculating the residual change rate of observation values of all satellites, carrying out reliability assessment on the output ionosphere scintillation intensity and scintillation intensity probability, and taking the ionosphere scintillation intensity and the reliability assessment result together as a scintillation monitoring result.
2. A combined multi-index high reliability ionospheric scintillation monitoring process as recited in claim 1, characterized in that the method further comprises: the method comprises the steps of calculating the residual change rate of an observation value of a satellite in real time by adopting a satellite orbit, a satellite clock error and a pseudo-range residual correction updated in real time, performing cross comparison with the output ionosphere scintillation intensity, adjusting the division standard of a scintillation intensity threshold according to a cross comparison result, and dynamically updating a random forest classifier.
3. A combined multi-index high reliability ionospheric scintillation monitoring process as recited in claim 1, the method is characterized by calculating a flicker intensity threshold in each satellite altitude interval, and comprising the following steps: based on the ROTI in the ionospheric scintillation intensity index, DI, classifying according to satellite altitude angles, and calculating the median and absolute median deviation in each satellite altitude angle interval corresponding to each scintillation intensity index; by selecting different absolute median deviation multiplication coefficients as ionospheric scintillation thresholds of different intensities.
4. The method of claim 1, further comprising calculating a scintillation intensity threshold for each satellite altitude interval, and further comprising: Based on AATR of the ionospheric scintillation intensity indices, it is mapped to an amplitude scintillation intensity index, specifically AATR corresponding to when weak scintillation occurs, denoted as a first dataset, and AATR corresponding to when strong scintillation occurs, denoted as a second dataset; the median of the first dataset was counted as a threshold for AATR to generate weak flicker and the median of the second dataset was counted as a threshold for AATR to generate strong flicker.
5. The method for combined multi-index high-reliability ionospheric scintillation monitoring as recited in claim 4, characterized in that the method comprises the following steps of, Training, DI and AATR separately, to obtain respective random forest classifiers, comprising: the ROTI in the ionospheric scintillation intensity index, DI, respectively using local time, satellite altitude angle and puncture point geomagnetic latitude as input parameters, using scintillation intensity corresponding to scintillation intensity threshold value in each satellite altitude angle interval as output parameters, obtaining ROTI by training, Random forest classifiers corresponding to DI respectively; And taking AATR of the ionosphere scintillation intensity indexes and local time and standard station geomagnetic latitude as input parameters, taking scintillation intensities corresponding to scintillation intensity thresholds in the first data set and the second data set as output parameters, and obtaining the AATR random forest classifier through training.
6. The method for combined multi-index high-confidence ionosphere scintillation monitoring of claim 1, wherein calculating the residual change rate of observations for all satellites comprises: Calculating observation value residuals of all satellites based on the acquired historical observation data of the reference station; And calculating the change rate of the observation value residual error by using the adjacent epoch residual error after correcting the influence of the satellite altitude angle based on the observation value residual error of all satellites.
7. The method for combined multi-index high-reliability ionosphere scintillation monitoring as in claim 1, wherein the reliability evaluation of the output ionosphere scintillation intensity and scintillation intensity probability comprises: Comparing the median and standard deviation of the residual error change rate of the observed value when different ionosphere scintillation intensities occur in the statistical data with the historical statistical value; and calculating the probability of flicker occurrence according to a cumulative distribution function of normal distribution, taking the probability as a reliability evaluation result, and taking the probability and the flicker intensity as a flicker monitoring result.
8. A combined multi-index high-confidence ionosphere scintillation monitoring system, comprising: The system comprises a scintillation intensity index acquisition module, a reference station acquisition module and a satellite acquisition module, wherein the scintillation intensity index acquisition module is used for calculating ionosphere scintillation intensity indexes of all satellites based on the acquired historical observation data of the reference station, and simultaneously acquiring phase scintillation intensity indexes and amplitude scintillation intensity indexes of all satellites acquired by a special ionosphere scintillation receiver; the scintillation intensity threshold calculation module is used for calculating the scintillation intensity threshold in each satellite height angle interval according to the relative distribution situation of ionosphere scintillation intensity indexes of all satellites and satellite height angles and combining the phase scintillation intensity indexes and the amplitude scintillation intensity indexes of all satellites; a training and outputting module for outputting the ROTI in the ionosphere scintillation intensity index based on the scintillation intensity threshold value in each satellite altitude interval, Training DI and AATR respectively to obtain respective random forest classifiers, and finally summarizing output results of all random forest classifiers to obtain ionosphere scintillation intensity and scintillation intensity probability; The scintillation monitoring module is used for calculating the residual change rate of the observation values of all satellites based on the acquired historical observation data of the reference station, carrying out reliability assessment on the output ionosphere scintillation intensity and scintillation intensity probability, and taking the ionosphere scintillation intensity and reliability assessment result together as a scintillation monitoring result.
9. An electronic device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the combined multi-index high-reliability ionosphere scintillation monitoring method of any one of claims 1-7.
10. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor implements the steps of the combined multi-index high-confidence ionosphere scintillation monitoring method of any one of claims 1-7.

Description

High-credibility ionosphere scintillation monitoring method and system combining multiple indexes Technical Field The invention belongs to the field of ionosphere scintillation monitoring, and particularly relates to a combined multi-index high-credibility ionosphere scintillation monitoring method and system. Background The phenomenon in which the ionosphere is affected by multiple factors to cause rapid and significant changes in physical characteristics such as internal electron density or ion composition is called ionosphere disturbance. The ionosphere disturbance of small scale is mainly ionosphere scintillation caused by plasma bubbles, the structure size is between tens of meters and tens of kilometers, the duration is short, but the ionosphere electron density in the region in the disturbance time period is changed by tens of times or even hundreds of times, the amplitude and phase stability of signals are disturbed, the tracking of a receiver to Beidou satellite signals is seriously influenced, cycle slip and even lock losing are caused, and meter-level errors are caused in positioning. The ionosphere scintillation monitoring method applied to the current stage scale mainly comprises the step of using one or more devices such as a special ionosphere scintillation monitor (Ionospheric Scintillation Monitor Receiver, ISMR), a general GNSS receiver, an ionosphere vertical measuring instrument, a satellite-based incoherent scattering radar, a ground-based radar and the like to form a monitoring network to detect an ionosphere and search an ionosphere mutation area. The phase scintillation intensity index can be calculated using ISMRAnd amplitude flicker intensity indexThe effect of ionospheric scintillation on the phase and amplitude of the receiver signal is monitored directly. The ionosphere monitoring network is used for monitoring and researching the ionosphere scintillation prone area by utilizing ISMR by national or regional institutions such as China space environment foundation comprehensive monitoring network, china academy of sciences national space science center, australia space environment prediction center, brazil ionosphere anomaly monitoring network, european ionosphere monitoring network, italian national geophysical and volcanic institute and the like so as to reduce the influence on production and life and the like. The ionospheric scintillation intensity monitoring index available to the earliest proposed universal receiver is the ionospheric total electron content change rate index (the Rate Of TEC Index, ROTI), and several subsequent scholars have proposed different ionospheric scintillation monitoring indices, including the electron content change rate index along the arc segment (the Along-ARC TEC RATE, AATR), the carrier-to-noise ratio S4 index [ ]) And ionospheric scintillation Doppler Index (DI), etc., which are all characterized by separating ionospheric delays in GNSS signals and calculating their intensity of variation over a window. In the scintillation event discrimination process, it is currently generally considered that 0.2<When <0.4, weak flicker occurs,Strong flicker occurs at >0.6, but whenWhen the dynamic precision single-point positioning error exceeds 0.8m, the user positioning application is seriously affected, and meanwhile, the current ionosphere scintillation state cannot be accurately judged through a fixed threshold value due to the reasons of satellite altitude angle, solar activity, plasma bubble motion and the like. Since the ROTI index cannot be clearly quantified to determine whether ionospheric scintillation occurs, multiple scholars choose to classify the scintillation by machine learning, but focus on using ISMR receiver observations, and use general receiver observations for scintillation discrimination has been less studied. Disclosure of Invention Aiming at the problems of high equipment cost, low layout density and the like in ionosphere scintillation monitoring by using ISMR, a vertical measuring instrument, a radar and other equipment, the invention uses a relatively low-cost, large-quantity and widely-distributed general GNSS receiver to perform ionosphere scintillation monitoring, but uses the ROTI obtained by monitoring by using the general GNSS receiver,Parameters such as DI and AATR have the problems of non-uniform scintillation monitoring intensity standard, inconsistent monitoring result, no credibility reference of the monitoring result and the like, and a combined multi-index high-credibility ionosphere scintillation monitoring method is provided, which is realized by combining ROTI,And indexes such as DI, AATR and the like, judging the current ionosphere scintillation intensity by using a random forest method, and giving intensity probability along with a judging result to provide credibility information for a user. According to an aspect of the present disclosure, there is provided a method for combined multi-index high-confidence