CN116148898A - Ionosphere scintillation index construction method based on GNSS Doppler observation value - Google Patents

Ionosphere scintillation index construction method based on GNSS Doppler observation value Download PDF

Info

Publication number
CN116148898A
CN116148898A CN202211095253.2A CN202211095253A CN116148898A CN 116148898 A CN116148898 A CN 116148898A CN 202211095253 A CN202211095253 A CN 202211095253A CN 116148898 A CN116148898 A CN 116148898A
Authority
CN
China
Prior art keywords
gnss
index
scintillation
satellite
doppler
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211095253.2A
Other languages
Chinese (zh)
Inventor
罗小敏
谢自春
曹光胤
李扬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences
Original Assignee
China University of Geosciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences filed Critical China University of Geosciences
Priority to CN202211095253.2A priority Critical patent/CN116148898A/en
Publication of CN116148898A publication Critical patent/CN116148898A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses an ionosphere scintillation index construction method based on GNSS Doppler observation values, which belongs to the field of space weather monitoring. And analyzing the DI availability by utilizing the ionosphere scintillation receiver products and the scintillation indexes constructed by the geodetic GNSS receivers, and finally giving the DI indexes of the GNSS receivers of different types. The method constructs the scintillation index based on the Doppler observed value of the common geodetic receiver which is widely distributed worldwide and is used for ionosphere scintillation monitoring, can reduce the ionosphere scintillation monitoring cost, and has the characteristics of simplicity, reliability and practicability.

Description

Ionosphere scintillation index construction method based on GNSS Doppler observation value
Technical Field
The invention belongs to the technical field of GNSS space weather monitoring, relates to a method for monitoring ionosphere scintillation by a global navigation satellite system (Global Navigation Satellite System, GNSS), and particularly relates to a method for constructing an ionosphere scintillation index based on GNSS Doppler observation values.
Background
Ionosphere flicker refers to the phenomenon of rapid random fluctuations in amplitude and phase of radio signals passing through ionosphere irregularities, which can affect GNSS high-precision positioning, navigation and time service. With the arrival of the high years of solar activity at the 25 th solar activity week, the probability and intensity of ionosphere scintillation occurrence which is closely related to the solar activity are also increased, so that the ionosphere scintillation monitoring and early warning are of great practical significance.
Ionospheric scintillation monitoring receiver (Ionospheric Scintillation Monitoring Receiver, ISMR) can provide an amplitude scintillation index S by 50Hz sampling 4 And phase scintillation index
Figure BDA0003832602610000011
However, due to the high layout cost, the global number of the mobile phones is only hundreds of. With the development of GNSS technology and the continuous construction of ground monitoring stations, international organizations and regional organizations represented by IGS/MGEX can provide GNSS actual measurement data of approximately 2000 sites worldwide, and can provide abundant observation information for ionospheric scintillation monitoring.
Pi in 1997 proposed the concept of the ionosphere total electron content time rate index (the rate of total electron content index, ROTI). Based on this, many scholars have proposed ionospheric scintillation indices based on GNSS phase observations. Mendillo proposed f in 2000 p An index; forte in 2005 proposed
Figure BDA0003832602610000012
An index; sanz has proposed AATR index in 2014 and Juan has proposed sigma in 2017 IF An index. In addition, luo constructed an amplitude flicker index S in 2020 using carrier-to-noise density observations of a 1Hz geodetic GNSS receiver 4c . The construction of these indices illustrates that it is feasible to base the global widely distributed geodetic receiver data for ionospheric scintillation monitoring.
The ROTI index is widely used in ionospheric scintillation monitoring at present, but in the data preprocessing process, the ROTI index needs to be built based on carrier phase observations of correct detection and repair cycle slip. Incorrect processing of cycle slips can result in the ROTI value being overestimated, affecting the reliability of the corresponding ionospheric scintillation monitoring. In addition, if the cycle slip repair is unsuccessful, the ROTI observation value may be removed, which is not beneficial to the continuity of ionosphere scintillation monitoring. How to effectively deal with cycle slips is still one of the difficulties in the GNSS field.
Disclosure of Invention
Aiming at the problem that the availability of the ROTI index is affected by frequent cycle slip in a strong ionosphere scintillation environment, the invention provides a method for constructing the ionosphere scintillation index based on the GNSS Doppler observation value by utilizing the characteristic that the Doppler observation value is not affected by cycle slip.
In order to achieve the above objective, the present invention provides a method for constructing ionosphere scintillation index based on GNSS doppler observations, comprising:
(1) Acquiring GNSS observation values of different latitude areas;
(2) Performing data preprocessing on the acquired GNSS observation values to obtain the altitude angles of all satellites above the measuring station, and setting a cut-off satellite altitude angle;
(3) Extracting ionosphere scintillation information from GNSS Doppler observation values above a cut-off satellite altitude angle, and constructing an ionosphere scintillation index DI from epoch to epoch by adopting a sliding window method;
(4) The availability of the ionospheric scintillation index DI is analyzed by utilizing the scintillation index constructed by the ionospheric scintillation receiver product and the geodetic GNSS receiver, and the threshold value of the DI index of different types of GNSS receivers is given.
In some alternative embodiments, step (2) comprises:
(2.1) finishing the combination and reading of GNSS high sampling rate observation value data;
(2.2) reading broadcast ephemeris parameters, and calculating the approximate position of the observation epoch satellite through the broadcast ephemeris parameters;
(2.3) obtaining WGS-84 accurate coordinates of the measuring station by a precise single-point positioning technology;
and (2.4) calculating the altitude of all satellites above the measuring station based on the approximate position of the satellites and the WGS-84 precise coordinates, setting the altitude of the cut-off satellite, and intercepting GNSS Doppler observation values above the altitude of the cut-off satellite.
In some alternative embodiments, the mathematical expression of the GNSS doppler observations is:
Figure BDA0003832602610000031
λ 1 and lambda is 2 Respectively the GNSS signal frequency f 1 And f 2 Wavelength of D (D) 1 For GNSS signal frequency f 1 Doppler observations, D 2 For GNSS signal frequency f 2 Doppler observations of>
Figure BDA0003832602610000032
For GNSS signal frequency f 1 Rate of change of geometrical distance between satellite and receiver>
Figure BDA0003832602610000033
For GNSS signal frequency f 2 Rate of change of geometrical distance between satellite and receiver>
Figure BDA0003832602610000034
And->
Figure BDA0003832602610000035
Rate of change of GNSS receiver clock and satellite clock, respectively, +.>
Figure BDA0003832602610000036
For GNSS signal frequency f 1 Rate of change of ionospheric delay error, < +.>
Figure BDA0003832602610000037
For GNSS signal frequency f 2 The rate of change of the ionospheric delay error under,
Figure BDA0003832602610000038
epsilon as the rate of change of tropospheric delay error 1 For GNSS signal frequency f 1 Multipath effects and errors, ε 2 For GNSS signal frequency f 2 The multipath effect and error under c is the propagation speed of light in vacuum.
In some alternative embodiments, step (3) comprises:
(3.1) combining the GNSS Doppler observed values above the cut-off satellite altitude into a double-frequency Doppler observed value combination delta D;
(3.2) selecting the length of the sliding window, wherein the time is 5 minutes;
(3.3) calculating the flicker index DI from the combination of dual-frequency doppler observations Δd epoch by epoch based on the sliding window.
In some alternative embodiments, Δd=λ 1 D 1 -λ 2 D 2 。
In some of the alternative embodiments of the present invention,
Figure BDA0003832602610000039
<>the average value is obtained.
In some alternative embodiments, step (4) comprises:
(4.1) calculating DI and ROTI indices for the geodetic receiver;
(4.2) calculating DI of the flash type receiver to obtain the flash index S 4 And (3) with
Figure BDA00038326026100000310
(4.3) comparing and analyzing each flicker index of different types of receivers, and counting DI of non-flicker days;
(4.4) constructing and obtaining the empirical threshold of DI under different types of receivers.
In general, the above technical solutions conceived by the present invention, compared with the prior art, enable the following beneficial effects to be obtained:
1. the ionosphere scintillation index DI based on the double-frequency Doppler observation value is constructed by using the GNSS Doppler observation value which is not affected by cycle slip, and has important practical value. Compared with the ROTI, the index can avoid the condition that the ROTI index is overestimated or discontinuous due to frequent cycle slip in the strong flicker environment.
2. Compared with a professional scintillation receiver, the ionosphere scintillation index DI constructed based on the Doppler observed value of the common geodetic GNSS receiver can reduce the ionosphere scintillation monitoring cost and has a certain practical value.
Drawings
Fig. 1 is a flow chart of a flicker index DI according to an embodiment of the present invention;
FIG. 2 is a time series of Doppler observations, ΔD and DI provided by an embodiment of the present invention;
FIG. 3 shows a flicker index DI and S according to an embodiment of the present invention 4 、
Figure BDA0003832602610000041
A comparison plot of ROTI;
FIG. 4 shows a flicker index DI and S according to an embodiment of the present invention 4 、
Figure BDA0003832602610000042
Correlation of ROTI.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention. In addition, the technical features of the embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
The specific flow of the embodiment of the invention is shown in fig. 1, and comprises the following processing steps:
step 1: collecting GNSS tracking station data of different latitude areas of the world based on a reference station network such as an IGS network, an MGEX network, a Beidou foundation enhancement network or a land state network;
step 2: the data preprocessing can be realized by the following steps:
(2.1) finishing the combination and reading of GNSS high sampling rate observation value data;
(2.2) reading broadcast ephemeris parameters, and calculating the approximate position of the observation epoch satellite through the broadcast ephemeris parameters;
(2.3) obtaining WGS-84 accurate coordinates of the measuring station by a precise single-point positioning technology;
and (2.4) calculating the altitude of all satellites above the measuring station based on the approximate position of the satellites and the WGS-84 precise coordinates, setting the altitude of the cut-off satellite, and intercepting GNSS Doppler observation values above the altitude of the cut-off satellite.
The mathematical expression of the Doppler observed value is as follows:
Figure BDA0003832602610000051
in the formula (1), lambda 1 And lambda is 2 Respectively the GNSS signal frequency f 1 And f 2 Is a wavelength of (c). D is the doppler observation.
Figure BDA0003832602610000052
Is the rate of change of the geometric distance between the satellite and the receiver. />
Figure BDA0003832602610000053
And->
Figure BDA0003832602610000054
The rates of change of the GNSS receiver clock bias and the satellite clock bias, respectively. />
Figure BDA0003832602610000055
And (3) with
Figure BDA0003832602610000056
The rates of change of ionospheric delay error and tropospheric delay error, respectively. Epsilon is the multipath effect and error.
Step 3: the construction of the flicker index DI can be realized specifically by the following means:
(3.1) combining the GNSS Doppler observed values above the cut-off satellite altitude angle into a double-frequency Doppler observed value combination delta D, wherein the mathematical expression is as follows:
ΔD＝λ 1 D 1 -λ 2 D 2 (2)
(3.2) selecting a sliding window length of 5 minutes;
(3.3) calculating a flicker index DI, the mathematical expression of DI being:
Figure BDA0003832602610000057
wherein, the expression is the average value taking operation.
Step 4: analyzing the availability of DI-index and establishing empirical thresholds based on DI built by different types of receivers can be achieved in particular by:
(4.1) calculating a geodetic receiver DI and a ROTI index, wherein the mathematical expression of ROTI is:
Figure BDA0003832602610000061
Figure BDA0003832602610000062
in the method, in the process of the invention,
Figure BDA0003832602610000063
is a GNSS carrier phase observation. i is the observation epoch. Δt is epoch difference.
(4.2) calculating a scintillation receiver DI to obtain a scintillation index S 4 And (3) with
Figure BDA0003832602610000064
S 4 And->
Figure BDA0003832602610000065
The mathematical expressions of (a) are respectively:
Figure BDA0003832602610000066
Figure BDA0003832602610000067
in SI det 、
Figure BDA0003832602610000068
The signal strength and carrier phase observations after detrending are respectively.
(4.3) performing comparative analysis of the flicker indexes of the different types of receivers (as shown in fig. 2 and 3), and counting DI of non-flicker days;
(4.4) constructing and obtaining the empirical threshold of DI under different types of receivers.
Step 5: analysis of DI index and commonly used scintillation index S 4 、
Figure BDA0003832602610000069
And the correlation of the ROTI (as shown in fig. 4).
According to the technical scheme, the ionosphere scintillation index DI is constructed by using the Doppler observed value which is not influenced by cycle slip, and compared with the common scintillation index ROTI, the DI index avoids the condition that the ROTI index is overestimated or discontinuous due to the cycle slip which frequently occurs in a strong scintillation environment.
It should be noted that each step/component described in the present application may be split into more steps/components, or two or more steps/components or part of the operations of the steps/components may be combined into new steps/components, as needed for implementation, to achieve the object of the present invention.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the invention and is not intended to limit the invention, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the invention are intended to be included within the scope of the invention.

Claims (7)

1. An ionosphere scintillation index construction method based on GNSS Doppler observation values is characterized by comprising the following steps:
(1) Acquiring GNSS observation values of different latitude areas;
(2) Performing data preprocessing on the acquired GNSS observation values to obtain the altitude angles of all satellites above the measuring station, and setting a cut-off satellite altitude angle;
(3) Extracting ionosphere scintillation information from GNSS Doppler observation values above a cut-off satellite altitude angle, and constructing an ionosphere scintillation index DI from epoch to epoch by adopting a sliding window method;
(4) The availability of the ionospheric scintillation index DI is analyzed by utilizing the scintillation index constructed by the ionospheric scintillation receiver product and the geodetic GNSS receiver, and the threshold value of the DI index of different types of GNSS receivers is given.
2. The method of claim 1, wherein step (2) comprises:
(2.1) finishing the combination and reading of GNSS high sampling rate observation value data;
(2.2) reading broadcast ephemeris parameters, and calculating the approximate position of the observation epoch satellite through the broadcast ephemeris parameters;
(2.3) obtaining WGS-84 accurate coordinates of the measuring station by a precise single-point positioning technology;
and (2.4) calculating the altitude of all satellites above the measuring station based on the approximate position of the satellites and the WGS-84 precise coordinates, setting the altitude of the cut-off satellite, and intercepting GNSS Doppler observation values above the altitude of the cut-off satellite.
3. The method of claim 2, wherein the mathematical expression of the GNSS doppler observations is:
Figure FDA0003832602600000011
λ 1 and lambda is 2 Respectively the GNSS signal frequency f 1 And f 2 Wavelength of D (D) 1 For GNSS signal frequency f 1 Doppler observations, D 2 For GNSS signal frequency f 2 Doppler observations of>
Figure FDA0003832602600000012
For GNSS signal frequency f 1 Rate of change of geometrical distance between satellite and receiver>
Figure FDA0003832602600000013
For GNSS signal frequency f 2 Rate of change of geometrical distance between satellite and receiver>
Figure FDA0003832602600000021
And->
Figure FDA0003832602600000022
Rate of change of GNSS receiver clock and satellite clock, respectively, +.>
Figure FDA0003832602600000023
For GNSS signal frequency f 1 Rate of change of ionospheric delay error, < +.>
Figure FDA0003832602600000024
For GNSS signal frequency f 2 Rate of change of ionospheric delay error, < +.>
Figure FDA0003832602600000025
Epsilon as the rate of change of tropospheric delay error 1 For GNSS signal frequency f 1 Multipath effects and errors, ε 2 For GNSS signal frequency f 2 The multipath effect and error under c is the propagation speed of light in vacuum.
4. A method according to claim 3, wherein step (3) comprises:
(3.1) combining the GNSS Doppler observed values above the cut-off satellite altitude into a double-frequency Doppler observed value combination delta D;
(3.2) selecting the length of the sliding window, wherein the time is 5 minutes;
(3.3) calculating the flicker index DI from the combination of the dual-frequency Doppler observations DeltaD epoch by epoch.
5. The method of claim 4, wherein Δd = λ 1 D 1 -λ 2 D 2 。
6. The method of claim 5, wherein the step of determining the position of the probe is performed,
Figure FDA0003832602600000026
<>the average value is obtained.
7. The method of claim 6, wherein step (4) comprises:
(4.1) calculating DI and ROTI indices for the geodetic receiver;
(4.2) calculating DI of the flash type receiver to obtain the flash index S 4 And (3) with
Figure FDA0003832602600000027
(4.3) comparing and analyzing each flicker index of different types of receivers, and counting DI of non-flicker days;
(4.4) constructing and obtaining the empirical threshold of DI under different types of receivers.
CN202211095253.2A 2022-09-05 2022-09-05 Ionosphere scintillation index construction method based on GNSS Doppler observation value Pending CN116148898A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211095253.2A CN116148898A (en) 2022-09-05 2022-09-05 Ionosphere scintillation index construction method based on GNSS Doppler observation value

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211095253.2A CN116148898A (en) 2022-09-05 2022-09-05 Ionosphere scintillation index construction method based on GNSS Doppler observation value

Publications (1)

Publication Number Publication Date
CN116148898A true CN116148898A (en) 2023-05-23

Family

ID=86358837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211095253.2A Pending CN116148898A (en) 2022-09-05 2022-09-05 Ionosphere scintillation index construction method based on GNSS Doppler observation value

Country Status (1)

Country Link
CN (1) CN116148898A (en)

Similar Documents

Publication Publication Date Title
CN101971047B (en) Device and method for the real-time monitoring of the integrity of a satellite navigation system
US9488729B2 (en) Detection of scintillations in signals of global navigation satellite systems caused by ionospheric irregularities
CN107942346B (en) A kind of high-precision GNSS ionized layer TEC observation extracting method
CN107607969B (en) Four-system pseudo range positioning method based on DCB correction
CN102998681A (en) High-frequency clock error estimation method of satellite navigation system
CN116243591B (en) Subnanosecond time service method integrating UTC (k) and Beidou broadcast ephemeris
CN114397425B (en) GNSS-IR soil humidity inversion method based on generalized extension approximation
CN113805208A (en) GNSS-IR height measurement method suitable for navigation receiver
Nie et al. Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS
Ge et al. Mitigation of the multipath effect in BDS-based time transfer using a wave-absorbing shield
Sinha et al. Ionospheric scintillation analysis using ROT and ROTI for slip cycle detection
Tseng et al. Determination of near real-time GNSS satellite clocks for the FORMOSAT-7/COSMIC-2 satellite mission
Yin et al. A novel cycle slips detection model for the high precision positioning
Elgered et al. Assessment of GNSS stations using atmospheric horizontal gradients and microwave radiometry
CN116148898A (en) Ionosphere scintillation index construction method based on GNSS Doppler observation value
CN110554412A (en) GNSS original data quality correction method for slope monitoring
CN114217332A (en) Modeling method and system of GNSS precision positioning random model considering flicker index
CN115902968A (en) PPP terminal positioning method based on Beidou third GEO broadcast enhancement information
CN113156477B (en) High-precision RTK positioning method based on Android smart phone
Hu et al. Towards GNSS Ambiguity Resolution for Smartphones in Realistic Environments: Characterization of Smartphone Ambiguities with RTK, PPP, and PPP-RTK
CN112666575A (en) Ionized layer abnormal disturbance detection method based on Beidou GEO satellite
Gill GNSS Precise Point Positioning using low-cost GNSS receivers
Chen et al. Feasibility and performance evaluation of low-cost GNSS devices for sea level measurement based on GNSS-IR
Pan et al. Machine learning-based multipath modelling in spatial-domain: a demonstration on GNSS short baseline processing
Hankansurijat et al. Atmospheric water estimation using GNSS precise point positioning method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination