RU2379709C1 - Method of determining propagation speed and direction of arrival of ionospheric disturbance - Google Patents

Method of determining propagation speed and direction of arrival of ionospheric disturbance Download PDF

Info

Publication number
RU2379709C1
RU2379709C1 RU2008125941/09A RU2008125941A RU2379709C1 RU 2379709 C1 RU2379709 C1 RU 2379709C1 RU 2008125941/09 A RU2008125941/09 A RU 2008125941/09A RU 2008125941 A RU2008125941 A RU 2008125941A RU 2379709 C1 RU2379709 C1 RU 2379709C1
Authority
RU
Russia
Prior art keywords
arrival
ionospheric
ionospheric disturbance
disturbance
values
Prior art date
Application number
RU2008125941/09A
Other languages
Russian (ru)
Inventor
Эдуард Леонтьевич Афраймович (RU)
Эдуард Леонтьевич Афраймович
Владислав Викторович Кирюшкин (RU)
Владислав Викторович Кирюшкин
Артур Амангалиевич Дисенов (RU)
Артур Амангалиевич Дисенов
Original Assignee
Государственное образовательное учреждение высшего профессионального образования "Иркутское высшее военное авиационное инженерное училище (военный институт)" Министерства обороны Российской Федерации
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 Государственное образовательное учреждение высшего профессионального образования "Иркутское высшее военное авиационное инженерное училище (военный институт)" Министерства обороны Российской Федерации filed Critical Государственное образовательное учреждение высшего профессионального образования "Иркутское высшее военное авиационное инженерное училище (военный институт)" Министерства обороны Российской Федерации
Priority to RU2008125941/09A priority Critical patent/RU2379709C1/en
Application granted granted Critical
Publication of RU2379709C1 publication Critical patent/RU2379709C1/en

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

FIELD: physics.
SUBSTANCE: invention relates to radiophysics and can be used in missile defense systems and airspace monitoring. Time series of complete electronic content obtained from measurements of double frequency receivers of satellite radio navigation systems GLONASS/GPS, which form an extended reception array, are filtered to separate variations corresponding to response of the ionosphere to effect of a source. Hypotheses of direction of arrival and speed of propagation values of a plane front of ionospheric disturbance are verified by forming a directional pattern of the reception array and its scanning in the given "direction of arrival - propagation speed" scan sector through synthesis of an output signal of the reception array with coherent summation of variations of series of complete electronic content of separate elements of the array with time shifting, calculated from verified values of propagation speed and direction of arrival of ionospheric disturbance and distance traversed by the front of ionospheric disturbance between elements of the reception array in the verified direction inside the spherical sphere of the ionosphere of the Earth. Correctness of the verified hypothesis and detection of inospheric disturbance is indicated by the total signal exceeding a given threshold value. Corresponding values of direction of arrival and phase speed of propagation of the ionospheric disturbance are estimated values.
EFFECT: increased sensitivity of detection and accuracy of determining propagation speed and direction of arrival of ionospheric disturbance detected by an extended array of reception stations of satellite radionavigation systems GLONASS/GPS.
3 dwg

Description

Изобретение относится к области радиофизики и может быть использовано в системах дистанционного контроля ядерных и иных взрывов, предупреждения о запусках ракет, наблюдения за сейсмической активностью.The invention relates to the field of radiophysics and can be used in remote control systems for nuclear and other explosions, warnings about missile launches, monitoring of seismic activity.

Существует способ определения направления прихода и скорости перемещения ионосферных возмущений естественного и техногенного характера, основанный на анализе данных о полном электронном содержании (ПЭС) в ионосфере Земли [1]. Он заключается в том, что с помощью решетки, состоящей из пространственно разнесенных приемников спутниковых радионавигационных систем (СРНС), получают временные ряды флуктуации ПЭС в ионосфере Земли для лучей «приемник - навигационный спутник». В результате фильтрации рядов ПЭС для отдельных лучей выделяют характерные возмущения, содержащие отклик ионосферы на воздействие источника. При этом параметры движения фронта волны ионосферного возмущения определяются по значениям взаимных задержек флуктуации ПЭС, зарегистрированных в трех точках. В указанном способе решение об обнаружении ионосферного возмущения принимается по анализу данных каждого приемника в отдельности, без учета информации двух других приемников, а для определения направления прихода волны ионосферного возмущения используется минимально необходимое количество приемников, поэтому данный способ характеризуется низкой чувствительностью и недостаточной точностью измерения направления прихода волны ионосферного возмущения.There is a way to determine the direction of arrival and the speed of movement of ionospheric disturbances of a natural and technogenic nature, based on the analysis of data on the total electronic content (TEC) in the Earth’s ionosphere [1]. It consists in the fact that using a lattice consisting of spatially spaced receivers of satellite radio navigation systems (SRNS), one obtains time series of fluctuations of TEC in the Earth’s ionosphere for “receiver-navigation satellite” beams. As a result of filtering the rows of TECs for individual rays, characteristic perturbations are identified that contain the response of the ionosphere to the action of the source. In this case, the parameters of the motion of the wave front of the ionospheric disturbance are determined from the values of the mutual delays of the TEC fluctuations recorded at three points. In this method, the decision to detect an ionospheric disturbance is made by analyzing the data of each receiver separately, without taking into account the information of the other two receivers, and the minimum number of receivers is used to determine the direction of arrival of the ionospheric disturbance wave, therefore this method is characterized by low sensitivity and insufficient accuracy of measuring the direction of arrival waves of ionospheric disturbance.

Наиболее близким к решению поставленной задачи, принятым за прототип, является способ обнаружения и измерения фазовой скорости и направления прихода ионосферного возмущения, основанный на анализе данных ПЭС ионосферы от большого количества пространственно разнесенных наземных двухчастотных приемников СРНС ГЛОНАСС/GPS. Сущность его заключается в том, что временные ряды ПЭС, полученные с помощью решетки двухчастотных приемников СРНС ГЛОНАСС/GPS с размером апертуры, соответствующей дальней зоне источника ионосферных возмущений, фильтруют с целью выделения вариаций ПЭС, которые соответствуют отклику ионосферы на воздействие источника. Полученные после фильтрации сигналы сдвигают по времени на величину Δτi и когерентно суммируют, добиваясь максимальной амплитуды суммарного сигнала. Решение об обнаружении возмущения принимают при превышении суммарным сигналом порогового уровня, а направление прихода и величину фазовой скорости волны возмущения находят из решения системы уравнений, описывающих в избранной системе координат семейство плоских волновых фронтов, отстоящих друг от друга на расстояние Δρi, определяемое относительными временными сдвигами Δτi сигналов отдельных приемников [2]. Это позволяет повысить чувствительность обнаружения и точность определения направления прихода ионосферного возмущения. В силу того, что в способе-прототипе расстояния Δρi, пройденные фронтом волны ионосферного возмущения, рассчитываются в топоцентрической системе координат в предположении прямолинейного распространения ионосферного возмущения, использование этого способа в протяженных приемных решетках приведет к ухудшению чувствительности обнаружения ионосферного возмущения и снижению точности определения его пространственно-временных параметров. Это обусловлено влиянием кривизны Земли на точность определения расстояния Δρi в протяженных приемных решетках, что в свою очередь приводит к повышению погрешности определения временных сдвигов Δτi, увеличению степени несинфазности при сложении рядов вариаций ПЭС и, как следствие, к снижению амплитуды суммарного сигнала и точности определения направления прихода и скорости распространения ионосферного возмущения.Closest to the solution of the problem, adopted as a prototype, is a method for detecting and measuring the phase velocity and direction of arrival of the ionospheric disturbance, based on the analysis of ionospheric TEC data from a large number of spatially separated two-frequency ground-based GLONASS / GPS SRNS receivers. Its essence lies in the fact that the time series of TECs obtained using the array of two-frequency receivers of the SRNS GLONASS / GPS with the aperture size corresponding to the far zone of the source of ionospheric disturbances is filtered to distinguish TEC variations that correspond to the response of the ionosphere to the source. The signals obtained after filtering are shifted in time by Δτ i and are coherently summed, achieving the maximum amplitude of the total signal. The decision to detect disturbances is made when the total signal exceeds the threshold level, and the direction of arrival and the magnitude of the phase velocity of the disturbance wave are found from the solution of the system of equations describing in the chosen coordinate system a family of plane wave fronts separated by a distance Δρ i determined by relative time shifts Δτ i signals of individual receivers [2]. This makes it possible to increase the detection sensitivity and the accuracy of determining the direction of arrival of the ionospheric disturbance. Due to the fact that in the prototype method, the distances Δρ i traveled by the front of the ionospheric disturbance wave are calculated in the topocentric coordinate system under the assumption that the ionospheric disturbance is linearly distributed, the use of this method in extended receiving arrays will lead to a decrease in the sensitivity of detection of the ionospheric disturbance and reduce the accuracy of determining it spatio-temporal parameters. This is due to the influence of the Earth's curvature on the accuracy of determining the distance Δρ i in extended receiving gratings, which in turn leads to an increase in the error in determining the time shifts Δτ i , to an increase in the degree of non-phase state when adding series of TEC variations and, as a result, to a decrease in the amplitude of the total signal and accuracy determining the direction of arrival and the propagation velocity of the ionospheric disturbance.

Задачей изобретения является повышение чувствительности обнаружения и точности определения скорости распространения и направления прихода ионосферного возмущения, регистрируемого протяженной решеткой приемных станций СРНС ГЛОНАСС/GPS.The objective of the invention is to increase the detection sensitivity and accuracy of determining the propagation velocity and direction of arrival of the ionospheric disturbance recorded by the extended array of receiving stations SRNS GLONASS / GPS.

В сравнении со способом-прототипом решение данной задачи достигается за счет того, что используется протяженная приемная решетка, угловой размер апертуры γ которой в направлении распространения ионосферного возмущения определяется из выраженияIn comparison with the prototype method, the solution to this problem is achieved due to the fact that an extended receiving grating is used, the angular size of the aperture γ of which in the direction of propagation of the ionospheric disturbance is determined from the expression

Figure 00000001
,
Figure 00000001
,

где δ - предельно допустимая относительная погрешность определения расстояния, пройденного фронтом волны ионосферного возмущения. Формируется диаграмма направленности приемной решетки и последовательно ориентируется во всех направлениях заданного сектора обзора в фазовом пространстве [α, V], где α - направление прихода, а V - скорость перемещения ионосферного возмущения. Ориентирование диаграммы направленности приемной решетки в заданном направлении (α, V) осуществляется путем формирования выходного сигнала, который является результатом пространственно-временного суммирования полученных после фильтрации рядов вариаций ПЭС отдельных элементов приемной решетки со сдвигами Δτi (i - номер элемента решетки), зависящими от заданных значений α и V и расстояний Δρi, пройденных плоским фронтом ионосферного возмущения между элементами приемной решетки в заданном направлении α внутри сферического слоя ионосферы Земли.where δ is the maximum permissible relative error in determining the distance traveled by the front of the ionospheric disturbance wave. The radiation pattern of the receiving array is formed and sequentially oriented in all directions of a given viewing sector in the phase space [α, V], where α is the direction of arrival and V is the velocity of movement of the ionospheric disturbance. The directional pattern of the receiving array in a given direction (α, V) is carried out by generating an output signal, which is the result of the spatio-temporal summation of the series of TEC variations of individual elements of the receiving array obtained after filtering with shifts Δτ i (i is the number of the array element), depending on the given values of α and V and the distances Δρ i traveled by the plane front of the ionospheric disturbance between the elements of the receiving array in the given direction α inside the spherical layer of ionospheres s of the Earth.

Решение об обнаружении ионосферного возмущения принимается при превышении суммарным сигналом заданного порогового уровня. Соответствующие значения α и V, для которых суммарный сигнал приемной решетки превысил пороговое значение, считаются оценками направления прихода и фазовой скорости распространения ионосферного возмущения. Предложенный способ обеспечивает обнаружение и определение пространственно-временных параметров ионосферных возмущений, одновременно пришедших на приемную решетку с различных направлений.The decision to detect an ionospheric disturbance is made when the total signal exceeds a predetermined threshold level. The corresponding values of α and V, for which the total signal of the receiving grating exceeded the threshold value, are considered as estimates of the direction of arrival and the phase velocity of propagation of the ionospheric disturbance. The proposed method provides for the detection and determination of spatio-temporal parameters of ionospheric disturbances that simultaneously arrive at the receiving array from various directions.

На фиг.1 представлена геометрия определения координат удаленного точечного источника ионосферного возмущения, на фиг.2 показаны вариации полного электронного содержания ΔIi(t) для отдельных элементов протяженной приемной решетки, а также суммарный сигнал

Figure 00000002
, полученный на выходе решетки. На фиг.3 приведена зависимость амплитуды суммарного сигнала от проверяемых значений направления прихода и скорости распространения ионосферного возмущения, полученная в результате сканирования диаграммы направленности приемной решетки в пределах сектора обзора.Figure 1 shows the geometry of determining the coordinates of a remote point source of ionospheric disturbance, figure 2 shows the variations in the total electronic content ΔI i (t) for individual elements of an extended receiving array, as well as the total signal
Figure 00000002
obtained at the output of the lattice. Figure 3 shows the dependence of the amplitude of the total signal from the verified values of the direction of arrival and the propagation velocity of the ionospheric disturbance obtained by scanning the radiation pattern of the receiving array within the viewing sector.

В предложенном способе определение значения ПЭС ионосферы осуществляется по двухчастотным измерениям псевдодальности между навигационным спутником и наземным приемником СРНС [1, 3]:In the proposed method, the determination of the ionospheric TEC value is carried out by two-frequency pseudorange measurements between the navigation satellite and the ground-based receiver of the SRNS [1, 3]:

Figure 00000003
Figure 00000003

где f1=1575.42 МГц, f2=1227.6 МГц, λ1, λ2 - частоты и длины волн навигационных радиосигналов СРНС; L1λ1 и L2λ2 - фазовый путь трансионосферных навигационных радиосигналов СРНС (L1, L2 - число полных оборотов фазы); z - зенитный угол луча «приемник - навигационный спутник».where f 1 = 1575.42 MHz, f 2 = 1227.6 MHz, λ 1 , λ 2 - frequencies and wavelengths of navigation radio signals SRNS; L 1 λ 1 and L 2 λ 2 - phase path of the transionospheric navigation radio signals SRNS (L 1 , L 2 - the number of full rotations of the phase); z is the zenith angle of the beam "receiver - navigation satellite."

Для устранения пространственной неопределенности, вызванной интегральным характером ПЭС, значение ПЭС соотносится к одной ионосферной точке (ИТ) - точке пересечения луча «приемник - навигационный спутник» с плоскостью на высоте hmax максимума ионизации F-области ионосферы (поскольку эта область вносит основной вклад в формирование ПЭС). На основании известных географических координат приемных станций СРНС, значения азимута и угла места луча «приемник - навигационный спутник» преобразуются в геоцентрические координаты (x, y, z) соответствующей ИТ.To eliminate the spatial uncertainty caused by the integral character of the TEC, the TEC value is related to one ionospheric point (IT) - the point of intersection of the receiver – navigation satellite beam with a plane at a height h max of ionization of the ionosphere F region (since this region makes the main contribution to formation of PES). Based on the known geographical coordinates of the SRNS receiving stations, the azimuth and elevation angle of the beam “receiver - navigation satellite” are converted to the geocentric coordinates (x, y, z) of the corresponding IT.

Совокупность лучей «приемник - навигационный спутник» в заданном регионе образует приемную решетку, каждый i-й элемент которой в момент времени t характеризуется измеренным значением ПЭС Ii(t) и положением соответствующей ионосферной точки xi(t), yi(t) и zi(t). Временные ряды ПЭС отражают как регулярные изменения ПЭС в точке регистрации, так и вариации ПЭС, вызванные ионосферными возмущениями различного характера.The set of receiver-navigation satellite beams in a given region forms a receiving array, each i-th element of which at time t is characterized by a measured TEC value I i (t) and the position of the corresponding ionospheric point x i (t), y i (t) and z i (t). The time series of TECs reflect both regular changes in TECs at the registration point and variations in TECs caused by ionospheric disturbances of a different nature.

Для выделения характерных ионосферных возмущений ряды ПЭС подвергаются процедуре специальной фильтрации в диапазоне периодов, соответствующих масштабу возмущения.To distinguish characteristic ionospheric disturbances, the TEC series undergo a special filtering procedure in the range of periods corresponding to the scale of the disturbance.

Обнаружение и определение пространственно-временных параметров ионосферного возмущения осуществляется путем последовательной проверки гипотез о значениях направления прихода и скорости распространения ионосферного возмущения. Для каждой пары проверяемых значений (α, V) формируется диаграмма направленности приемной решетки и соответствующим образом ориентируется в фазовом пространстве [α, V] за счет синфазного суммирования отдельных рядов ΔIi(t) приемной решетки к некоторому центральному ряду ΔI0(t), выбранному в качестве опорного, с временными сдвигами τi и формирования выходного сигнала приемной решетки:The detection and determination of the spatio-temporal parameters of the ionospheric disturbance is carried out by sequentially testing hypotheses about the values of the direction of arrival and propagation velocity of the ionospheric disturbance. For each pair of checked values (α, V), the radiation pattern of the receiver array is formed and accordingly oriented in the phase space [α, V] due to the in-phase summation of the individual rows ΔI i (t) of the receiver array to a certain central row ΔI 0 (t), selected as a reference, with time shifts τ i and the formation of the output signal of the receiving array:

Figure 00000004
Figure 00000004

где Р - количество элементов приемной решетки.where P is the number of elements of the receiving grid.

Временной сдвиг τi определяется как разность времени tj j-го отсчета i-го суммируемого ряда ПЭС и времени t0 регистрации ионосферного возмущения центральным элементом приемной решетки τi=tj-t0 и выбирается исходя из минимизации выражения, описывающего динамику распространения возмущения:The time shift τ i is defined as the difference between the time t j of the jth count of the i-th summed TEC series and the time t 0 of recording the ionospheric disturbance by the central element of the receiving array τ i = t j -t 0 and is selected based on minimizing the expression describing the dynamics of the propagation of the disturbance :

Figure 00000005
Figure 00000005

где Δρi - расстояние, пройденное фронтом волны между i-м и центральным элементом приемной решетки.where Δρ i is the distance traveled by the wave front between the ith and the central element of the receiving grating.

Для протяженных приемных решеток расстояние Δρi рассчитывается с учетом кривизны Земли. С этой целью в заданном направлении α прихода волны ионосферного возмущения на высоте hmax задается удаленный точечный источник (обозначен точкой Е на фиг.1), который будет являться полюсом ортодромической системы координат, экватор которой (сплошная жирная линия на фиг.1) проходит через центральный элемент приемной решетки (точка А на фиг.1). Тогда фронт волны, распространяющийся от удаленного точечного источника и проходящий через i-й элемент приемной решетки (точка В на фиг.1), будет представлять собой широтный круг (жирная прерывистая линия), параллельный экватору полученной ортодромической системы. Такая модель соответствует плоской волне ионосферного возмущения, распространяющейся на сфере Земли.For extended receiving gratings, the distance Δρ i is calculated taking into account the curvature of the Earth. To this end, in a given direction α of the arrival of the ionospheric disturbance wave at a height h max, a remote point source (indicated by point E in Fig. 1) is set, which will be the pole of the orthodromic coordinate system, the equator of which (solid bold line in Fig. 1) passes through the central element of the receiving grate (point A in figure 1). Then the wave front propagating from a remote point source and passing through the ith element of the receiving array (point B in Fig. 1) will be a latitudinal circle (a thick broken line) parallel to the equator of the resulting orthodromic system. Such a model corresponds to a plane wave of ionospheric disturbance propagating over the Earth's sphere.

Геоцентрические координаты (xe, ye, ze) удаленного источника ионосферного возмущения определяются с использованием правил сферической тригонометрии. При этом рассматривается сферический треугольник, вершиной А которого является центральный элемент приемной решетки с известными координатами (x0, y0, z0). Вершиной С этого треугольника является северный полюс геоцентрической системы координат (0, 0, R+hmax), где R - радиус Земли. Необходимо определить координаты третьей вершины Е, которая и будет являться удаленным источником. Чтобы удаленный источник Е являлся полюсом ортодромической системы координат, угловой размер стороны АЕ сферического треугольника задается равным π/2. В полученном сферическом треугольнике известны две стороны АС и АЕ, а также угол между ними ∠CAE=α, что является типовой задачей решения сферического треугольника. С использованием теоремы косинусов сторон сферического треугольника определяется третья сторона и координаты (xe, ye, ze) удаленного источника Е.The geocentric coordinates (x e , y e , z e ) of a remote source of ionospheric disturbance are determined using the rules of spherical trigonometry. In this case, we consider a spherical triangle, the vertex A of which is the central element of the receiving grating with known coordinates (x 0 , y 0 , z 0 ). The vertex C of this triangle is the north pole of the geocentric coordinate system (0, 0, R + h max ), where R is the radius of the Earth. It is necessary to determine the coordinates of the third vertex E, which will be the remote source. So that the remote source E is the pole of the orthodromic coordinate system, the angular size of the side AE of the spherical triangle is set to π / 2. In the obtained spherical triangle, two sides AC and AE are known, as well as the angle between them ∠CAE = α, which is a typical problem for solving a spherical triangle. Using the cosine theorem of the sides of a spherical triangle, the third side and the coordinates (x e , y e , z e ) of the remote source E.

Расстояние, пройденное фронтом волны между i-м и центральным элементом приемной решетки определяется как разность расстояний АЕ и BE (фиг.1) и записывается в виде:The distance traveled by the wave front between the ith and the central element of the receiving array is determined as the difference between the distances AE and BE (Fig. 1) and is written as:

Figure 00000006
Figure 00000006

где (xi, yi, zi) координаты i-го элемента приемной решетки в момент времени tj.where (x i , y i , z i ) are the coordinates of the ith element of the receiver array at time t j .

Решение о правильности проверяемой гипотезы принимается при превышении суммарным сигналом заданного порогового уровня. При этом считается, что обнаружено ионосферное возмущение, а соответствующие значения α и V, для которых суммарный сигнал приемной решетки превысил пороговое значение, считаются оценками направления прихода и фазовой скорости распространения обнаруженного ионосферного возмущения.The decision on the correctness of the tested hypothesis is made when the total signal exceeds a predetermined threshold level. In this case, it is believed that an ionospheric disturbance is detected, and the corresponding values of α and V, for which the total signal of the receiving array exceeded the threshold value, are considered as estimates of the direction of arrival and the phase propagation velocity of the detected ionospheric disturbance.

Таким образом, предложенный способ обеспечивает более точное, с учетом кривизны Земли, определение расстояния Δρi в протяженных приемных решетках, что в свою очередь приводит к повышению точности определения временных сдвигов Δτi, повышению амплитуды суммарного сигнала за счет увеличения степени синфазности сложения рядов вариаций ПЭС и, как следствие, к повышению чувствительности обнаружения ионосферного возмущения и повышению точности определения направления прихода и скорости распространения ионосферного возмущения.Thus, the proposed method provides a more accurate, taking into account the curvature of the Earth, determination of the distance Δρ i in extended receiving gratings, which in turn leads to an increase in the accuracy of determining the time shifts Δτ i , to an increase in the amplitude of the total signal due to an increase in the degree of in-phase addition of series of TEC variations and, as a result, to an increase in the sensitivity of detection of an ionospheric disturbance and an increase in the accuracy of determining the direction of arrival and propagation velocity of an ionospheric disturbance.

Кроме того, проверка всех возможных гипотез о значениях скорости распространения и направления прихода волны в заданном секторе обзора обеспечивает обнаружение и определение пространственно-временных параметров ионосферных возмущений, одновременно пришедших на приемную решетку с различных направлений.In addition, verification of all possible hypotheses about the values of the velocity of propagation and the direction of arrival of the wave in a given viewing sector ensures the detection and determination of the spatio-temporal parameters of ionospheric disturbances that simultaneously arrive at the receiving array from different directions.

ЛИТЕРАТУРАLITERATURE

1. Afraimovich, E.L., E.A.Kosogorov, N.P.Perevalova, and A.V.Plotnikov, The use of GPS arrays in detecting shock-acoustic waves generated during rokcet launchings J. Atmos. Solar-Terr. Phys., V.63, 1941-1957, 2001.1. Afraimovich, E.L., E. A. Kosogorov, N. P. Perevalova, and A. V. Plotnikov, The use of GPS arrays in detecting shock-acoustic waves generated during rokcet launchings J. Atmos. Solar-Terr. Phys., V.63, 1941-1957, 2001.

2. Афраймович Э.Л., Чернухов В.В., Кирюшкин В.В. Способ обнаружения, измерения фазовой скорости и направление прихода ионосферного возмущения. Патент Российской Федерации №2189052, зарегистрирован в Государственном реестре изобретений Российской Федерации, г.Москва, 10 сентября 2002 г., (прототип).2. Afraimovich E.L., Chernukhov V.V., Kiryushkin V.V. A method for detecting, measuring the phase velocity and the direction of arrival of the ionospheric disturbance. Patent of the Russian Federation No. 2189052, registered in the State register of inventions of the Russian Federation, Moscow, September 10, 2002, (prototype).

3. Calais E., Minster J.B. GPS detection of ionospheric perturbations following a Space Shuttle ascent // Geophys. Res. Lett. V.23. P.1897-1900, 1996.3. Calais E., Minster J.B. GPS detection of ionospheric perturbations following a Space Shuttle ascent // Geophys. Res. Lett. V.23. P.1897-1900, 1996.

Claims (1)

Способ определения скорости распространения и направления прихода ионосферного возмущения, основанный на анализе данных о полном электронном содержании в ионосфере Земли, которые получают в результате обработки сигналов, принятых двухчастотными приемниками спутниковой радионавигационной системой ГЛОНАСС/GPS, с последующим формированием временных рядов полного электронного содержания и их фильтрацией в диапазоне периодов колебаний, соответствующих отклику ионосферы на воздействие источника ионосферного возмущения, отличающийся тем, что используют протяженную приемную решетку и последовательно проверяют гипотезу о значениях направления прихода и скорости распространения плоского фронта ионосферного возмущения путем формирования диаграммы направленности приемной решетки и ее сканирования в заданном секторе обзора пространства «направление прихода - скорость распространения ионосферного возмущения» за счет синтеза выходного сигнала приемной решетки при синфазном суммировании рядов вариаций полного электронного содержания отдельных элементов решетки с временными сдвигами, рассчитанными исходя из проверяемых значений направления прихода и скорости распространения ионосферного возмущения и расстояний, пройденных фронтом ионосферного возмущения между элементами приемной решетки в проверяемом направлении внутри сферического слоя ионосферы Земли, решение о правильности проверяемой гипотезы и обнаружении ионосферного возмущения принимается при превышении суммарным сигналом заданного порогового уровня, соответствующие значения направления прихода и фазовой скорости распространения ионосферного возмущения считаются оценочными значениями. A method for determining the propagation velocity and direction of arrival of the ionospheric disturbance, based on the analysis of data on the total electronic content in the Earth's ionosphere, which are obtained as a result of processing signals received by two-frequency receivers by the GLONASS / GPS satellite radio navigation system, followed by the formation of time series of the full electronic content and their filtering in the range of oscillation periods corresponding to the response of the ionosphere to the action of an ionospheric disturbance source, different t m, that they use an extended receiving grating and sequentially test the hypothesis about the values of the direction of arrival and the propagation velocity of the plane front of the ionospheric disturbance by forming the radiation pattern of the receiving lattice and scanning it in a given sector of the “direction of arrival - propagation velocity of the ionospheric disturbance” space due to the synthesis of the output signal of the receiving array during in-phase summation of the series of variations in the total electronic content of individual elements of the array with time by the measured shifts calculated on the basis of the verified values of the arrival direction and propagation velocity of the ionospheric disturbance and the distances traveled by the front of the ionospheric disturbance between the elements of the receiving lattice in the checked direction inside the spherical layer of the Earth’s ionosphere, the decision on the correctness of the hypothesis being tested and the detection of the ionospheric disturbance is made when the total signal exceeds the specified signal threshold level, the corresponding values of the direction of arrival and the phase velocity of propagation ionospheric disturbances are considered estimated values.
RU2008125941/09A 2008-06-25 2008-06-25 Method of determining propagation speed and direction of arrival of ionospheric disturbance RU2379709C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2008125941/09A RU2379709C1 (en) 2008-06-25 2008-06-25 Method of determining propagation speed and direction of arrival of ionospheric disturbance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2008125941/09A RU2379709C1 (en) 2008-06-25 2008-06-25 Method of determining propagation speed and direction of arrival of ionospheric disturbance

Publications (1)

Publication Number Publication Date
RU2379709C1 true RU2379709C1 (en) 2010-01-20

Family

ID=42120946

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2008125941/09A RU2379709C1 (en) 2008-06-25 2008-06-25 Method of determining propagation speed and direction of arrival of ionospheric disturbance

Country Status (1)

Country Link
RU (1) RU2379709C1 (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2560094C2 (en) * 2013-12-03 2015-08-20 Федеральное государственное бюджетное учреждение "Институт прикладной геофизики имени академика Е.К.Федорова" (ФГБУ "ИПГ") Method of determining propagation speed and direction of arrival of ionospheric perturbation
RU2560525C1 (en) * 2014-06-25 2015-08-20 Александр Васильевич Тертышников Method of determining position of epicentral area of source and propagation speed of travelling ionospheric disturbances
RU2624634C1 (en) * 2016-03-29 2017-07-05 Федеральное государственное казенное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф. Можайского" Министерства обороны Российской Федерации Method of determining speed of distribution and direction of ionospheric perturbation
RU2624911C1 (en) * 2016-06-28 2017-07-10 Федеральное государственное бюджетное учреждение науки Институт солнечно-земной физики Сибирского отделения Российской академии наук (ИСЗФ СО РАН) Method of detecting ionospheric pertrubations caused by launching the space ships
RU2648237C2 (en) * 2015-04-27 2018-03-23 МИНИСТЕРСТВО ОБОРОНЫ РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ КАЗЁННОЕ ВОЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ Военная академия ракетных войск стратегического назначения имени Петра Великого Method of wave front position measurement
RU2655164C2 (en) * 2016-07-04 2018-05-24 Федеральное государственное бюджетное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф.Можайского" Министерства обороны Российской Федерации System for determining propagation speed and direction of arrival of ionospheric disturbance
RU2683113C1 (en) * 2018-03-26 2019-03-26 Александр Васильевич Тертышников Method of determining characteristics of auroral ovals and state of magnetic field of earth
CN111158021A (en) * 2018-11-08 2020-05-15 千寻位置网络有限公司 Ionosphere interference estimation method and system and early warning terminal
CN112114334A (en) * 2020-09-01 2020-12-22 华东交通大学 GNSS ionosphere abnormal disturbance parameter extraction method
RU2787952C1 (en) * 2021-10-18 2023-01-13 Федеральное государственное унитарное предприятие "18 Центральный научно-исследовательский институт" Министерства обороны Российской Федерации Method for determining radio signal arrival direction

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2560094C2 (en) * 2013-12-03 2015-08-20 Федеральное государственное бюджетное учреждение "Институт прикладной геофизики имени академика Е.К.Федорова" (ФГБУ "ИПГ") Method of determining propagation speed and direction of arrival of ionospheric perturbation
RU2560525C1 (en) * 2014-06-25 2015-08-20 Александр Васильевич Тертышников Method of determining position of epicentral area of source and propagation speed of travelling ionospheric disturbances
RU2648237C2 (en) * 2015-04-27 2018-03-23 МИНИСТЕРСТВО ОБОРОНЫ РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ КАЗЁННОЕ ВОЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ Военная академия ракетных войск стратегического назначения имени Петра Великого Method of wave front position measurement
RU2624634C1 (en) * 2016-03-29 2017-07-05 Федеральное государственное казенное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф. Можайского" Министерства обороны Российской Федерации Method of determining speed of distribution and direction of ionospheric perturbation
RU2624911C1 (en) * 2016-06-28 2017-07-10 Федеральное государственное бюджетное учреждение науки Институт солнечно-земной физики Сибирского отделения Российской академии наук (ИСЗФ СО РАН) Method of detecting ionospheric pertrubations caused by launching the space ships
RU2655164C2 (en) * 2016-07-04 2018-05-24 Федеральное государственное бюджетное военное образовательное учреждение высшего образования "Военно-космическая академия имени А.Ф.Можайского" Министерства обороны Российской Федерации System for determining propagation speed and direction of arrival of ionospheric disturbance
RU2683113C1 (en) * 2018-03-26 2019-03-26 Александр Васильевич Тертышников Method of determining characteristics of auroral ovals and state of magnetic field of earth
CN111158021A (en) * 2018-11-08 2020-05-15 千寻位置网络有限公司 Ionosphere interference estimation method and system and early warning terminal
CN112114334A (en) * 2020-09-01 2020-12-22 华东交通大学 GNSS ionosphere abnormal disturbance parameter extraction method
RU2787952C1 (en) * 2021-10-18 2023-01-13 Федеральное государственное унитарное предприятие "18 Центральный научно-исследовательский институт" Министерства обороны Российской Федерации Method for determining radio signal arrival direction

Similar Documents

Publication Publication Date Title
RU2379709C1 (en) Method of determining propagation speed and direction of arrival of ionospheric disturbance
Groves et al. Intelligent urban positioning using multi-constellation GNSS with 3D mapping and NLOS signal detection
US9588246B2 (en) Data collection system, marine seismic survey system and method of estimating a corrected sound speed
Braasch Multipath
US20150116146A1 (en) Anomaly detection using an antenna baseline constraint
CN105334510A (en) GNSS-R ground surface detection apparatus and method
Bürgmann et al. Space geodesy: A revolution in crustal deformation measurements of tectonic processes
RU2615984C2 (en) Assessment method of arrival navigation signals direction to receiver after reflection from walls in satellite positioning system
Zhao et al. Investigation on underwater positioning stochastic model based on acoustic ray incidence angle
US7446705B1 (en) Method and apparatus for determining parameters for a parametric expression characterizing the phase of an acquired signal
KR102292187B1 (en) Data signal treating method for GNSS receiver, computer readable medium and GNSS receiver system
KR101784178B1 (en) Scatterometer system for ocean parameters monitoring
RU2560094C2 (en) Method of determining propagation speed and direction of arrival of ionospheric perturbation
CN110286354A (en) The method, apparatus and computer readable storage medium of multi-target detection and differentiation
CN113945955A (en) Method and system for improving sea surface measurement high precision based on atmospheric delay error correction
Cario et al. Analysis of error sources in underwater localization systems
Kbayer et al. Robust GNSS navigation in urban environments by bounding NLOS bias of GNSS pseudoranges using a 3D city model
Su et al. Distributed sensing of ionospheric irregularities with a GNSS receiver array
RU2436134C1 (en) Method for rapid investigation of atmosphere, earth's surface and ocean
Groves et al. Enhancing micro air vehicle navigation in dense urban areas using 3D mapping aided GNSS
RU2568937C2 (en) Space navigation system and method
Kunitsyn et al. Earthquake prediction research using radio tomography of the ionosphere
RU2536320C1 (en) Method of navigation of aircrafts
CN115047503A (en) Method and system for calibrating mirror reflection point of land-based satellite-borne GNSS (Global navigation satellite System) reflection signal
RU2189051C2 (en) Technique detecting ionospheric disturbance and locating position of its source

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20100626