СПОСОБ ОБНАРУЖЕНИЯ ИОНОСФЕРНОГО ВОЗМУЩЕНИЯ И ОПРЕДЕЛЕНИЯ МЕСТОПОЛОЖЕНИЯ ЕГО ИСТОЧНИКА Российский патент 2002 года по МПК G01S5/06 G01S13/95 

Изобретение относится к радиотехнике и может быть использовано в системах дистанционного контроля ядерных и иных взрывов, предупреждения о запусках ракет, наблюдения за сейсмической активностью.

Известны способы наблюдения ионосферных возмущений, порождаемых антропогенными и естественными источниками, основанные на регистрации задержек сигналов спутниковых радионавигационных систем ГЛОНАСС/GPS. С помощью нескольких приемников ГЛОНАСС/GPS осуществляют двухдиапазонные измерения задержек навигационных сигналов ГЛОНАСС/GPS, по измеренным задержкам определяют значение полного электронного содержания на трассе "приемник - навигационный искусственный спутник Земли". Производят фильтрацию рядов полного электронного содержания для отдельных приемников в выбранном диапазоне периодов колебаний и при наличии сигнала, превышающего заданный уровень, регистрируют возмущение, вызванное мощными наземными взрывами [1], запусками космического аппарата [2] . Данные способы отличаются низкой чувствительностью и невозможностью определения координат источника возмущений.

Наиболее близким к решению поставленной задачи является способ определения направления прихода и скорости перемещения ионосферных возмущений различной природы, основанный на анализе данных о полном электронном содержании в ионосфере Земли [3].

Сущность его заключается в том, что с помощью решетки, состоящей из трех пространственно разнесенных двухчастотных приемников сигналов спутниковых радионавигационных систем ГЛОНАСС/GPS получают временные ряды флуктуации полного электронного содержания в ионосфере Земли для линий "приемник - навигационный искусственный спутник Земли". В результате фильтрации рядов полного электронного содержания для отдельных приемников выделяют возмущения, содержащие отклик ионосферы на воздействие источника. При этом параметры движения фронта волны возмущений определяются по значениям взаимных задержек флуктуации полного электронного содержания, зарегистрированных в трех точках. В силу того, что при указанном способе решение об обнаружении возмущения принимается по анализу сигнала каждого приемника отдельно, без учета информации двух других приемников, а для определения направления прихода волны возмущения используется минимально необходимое число приемников n=3, способ-прототип характеризуется низкой чувствительностью обнаружения. Кроме того, этот способ обеспечивает определение направления прихода возмущения, но не позволяет определить местоположение источника возмущения. Это обусловлено тем, что расстояние между приемниками в способе выбрано малым, соответствующим дальней зоне источника. При этом фронт волны возмущения для приемной решетки мало отличается от плоского, а информация о местоположении источника возмущения, заложенная в величине кривизны волнового фронта, теряется.

Целью изобретения является повышение чувствительности обнаружения и точности определения местоположения источника возмущения. В сравнении со способом прототипом это достигается за счет того, что для реализации предложенного способа используется решетка с большим числом n(n>3) пространственно разнесенных двухчастотных приемников сигналов ГЛОНАСС/GPS. Апертура решетки выбирается таким образом, чтобы выполнялось условие зоны Френеля

где L - апертура решетки;
R - предполагаемое расстояние до источника возмущения;
λ - длина волны возмущения [3, 4, 5].

Сигналы, полученные в результате фильтрации рядов полного электронного содержания для отдельных приемников решетки, задерживают на время Δτ_i (i - номер приемника) относительно сигнала одного из приемников, выбранного в качестве опорного сигнала, и когерентно суммируют между собой. Решение об обнаружении возмущений принимают при превышении суммарным сигналом заданного порога. Местоположение источника возмущения находят из решения системы уравнений, описывающих в избранной системе координат семейство сферических волновых фронтов, отстоящих друг от друга на расстояние, определяемое относительными временными сдвигами Δτ_i сигналов отдельных приемников, соответствующими абсолютному максимуму суммарного сигнала.

На фиг. 1 в декартовой системе координат изображены точки расположения приемников, проекции подионосферных точек лучей "приемник - навигационный искусственный спутник Земли", семейство фазовых фронтов возмущения полного электронного содержания. На фиг.2 изображены вариации полного электронного содержания, полученные для отдельных приемников, соответствующие отклику ионосферы на воздействие источника возмущений, и суммарный сигнал всех приемников.

Исследованиями ряда авторов установлено, что откликом ионосферы на воздействие мощных антропогенных (ядерные и химические взрывы, запуски космических аппаратов) и естественных (землетрясения, извержения вулканов) источников являются распространяющиеся в направлении от источника ударно-акустические волны. Распространение ударно-акустической волны приводит к возмущениям полного электронного содержания в ионосфере. Возмущение полного электронного содержания имеет длину волны λ ≈ 100-200 км, период Т≈200-300 с, амплитуду 10¹⁵-10¹⁶ эл/м² [3, 4, 5].

Каждый из разнесенных в пространстве двухчастотных приемников ГЛОНАСС/GPS осуществляет регистрацию фазовых задержек l₁ и L₂ навигационных сигналов с несущей частотой f₁= 1575,42 МГц, f₂= 1227,6 МГц на трассе "приемник - навигационный искусственный спутник Земли".

На основе измеренных значений задержек l₁ и L₂ находят величину полного электронного содержания для указанной трассы:

где
Z_m - зенитный угол луча "приемник - навигационный искусственный спутник Земли", отсчитываемый от вертикального направления;
λ₁,λ₂ - длина волны навигационных сигналов с частотой f₁ и f₂ соответственно.

Флуктуации значений полного электронного содержания содержат как низкочастотные составляющие, обусловленные движением навигационного искусственного спутника Земли, широтно-долготным и суточным ходом полного электронного содержания, так и высокочастотные колебания ΔI(t), вызванные возмущением полного электронного содержания при распространении ударно-акустической волны от источника. Для выделения этих колебаний ряды полного электронного содержания с выхода каждого приемника подают на вход фильтра с полосой пропускания, соответствующей полосе частот ударно-акустической волны. Полученные отклики ΔI(t) относят к проекциям 2 (фиг.1) на горизонтальную плоскость точек в ионосфере, соответствующих пересечению луча "приемник - навигационный искусственный спутник Земли" с поверхностью на высоте максимума слоя F₂ ионосферы. Координаты х_i; у_i указанной точки находят, зная координаты приемника и навигационного искусственного спутника Земли, с помощью стандартных в системе GPS процедур.

Полученные с выходов фильтров сигналы ΔI_i(t) 3 (фиг.2) задерживают на время Δτ_i относительно сигнала опорного приемника ΔI₀(t) 4 (фиг.2) и затем суммируют. Значения Δτ_i подбирают так, чтобы обеспечить максимальную амплитуду суммарного сигнала ΔI_c(t) 5 (фиг.2). При превышении суммарным сигналом заданного порога принимается решение об обнаружении источника возмущения. Величина порога определяется выбранным критерием обнаружения.

Если размер апертуры решетки L и расстояние до источника возмущений R таковы, что выполняется условие зоны Френеля, то фронт ударно-акустической волны в раскрыве решетки приемников является сферическим. Значения Δτ_i определяют расстояния ρ_i между фазовыми фронтами, проходящими через точки х_i; у_i и имеющими форму окружностей с центром x_u; y_и, где
ρ_i = Δτ_i•V;
x_и; y_и - координаты источника возмущения;
V - фазовая скорость волны возмущения.

Расстояния ρ_i определяются относительно фазового фронта, проходящего через точку х₀; (у₀), соответствующую опорному приемнику. Пусть декартова система координат задана таким образом, что ее начало совпадает с точкой х₀; (у₀), ось Оу системы направлена на север, а ось Ox - на восток (фиг.1). В заданной системе координат уравнение фазового фронта имеет вид
(x-x_u)²+(y-y_u)² = τV,
где τ - время распространения возмущения.

Тогда для точек x_i можем записать



Неизвестные координаты источника возмущений (х_и; у_и) можно найти из численного решения данной системы уравнений одним из известных методов [6].

ЛИТЕРАТУРА
[l] Calais E. , Minster B.J., Hofton M.A., Hedlin M.A.H. Ionospheric signature of surface mine blasts from Global Positioning System measurements.//Geophys. J. Int. 1998. V. 132. P. 191-202.

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

[3] Э. Л.Афраймович, Е.А.Косогоров, А.В.Плотников. Детектирование с помощью GPS-решеток ударно-акустических волн, генерируемых при запуске ракет. Труды VI международной научно-технической конференции "Радиолокация, навигация, связь". - Воронеж, 25-27 апреля 2000 г. Том 1, с. 462-474. Прототип.

[4] Li Y.Q., Jacobson A.R., Carlos R.C., Massey R.S., Taranenko Y.N., Wu G. The blast wave of the Shuttle plume at ionospheric heights. // Geophys. Res. Lett. 1994. V. 21. P. 2737-2740.

[5] Blanc E., Jacobson A.R. Observation of ionospheric disturbances follwing a 5-kt chemical explosion. 2. Prolonged anomalies and stratifications in the lower thermosphere after shock passage. //Radio Science. 1989. V. 24. P. 739-746.

[6] Н. С. Бахвалов, Н. П.Жидков, Г.М.Кобельков. Численные методы. - М. Наука". 1987.

Использование: в системах дистанционного контроля ядерных и иных взрывов, предупреждения о запусках ракет, наблюдения за сейсмической активностью. Сущность изобретения: временные ряды полного электронного содержания, полученные с помощью пространственной решетки двухчастотных приемников сигналов спутниковых радионавигационных систем ГЛОНАСС/GPS с размером апертуры, соответствующей зоне Френеля источника возмущений, фильтруют с целью выделения вариаций полного электронного содержания, соответствующих отклику ионосферы на воздействие источника, полученные после фильтрации сигналы сдвигают по времени на величину Δτ_i и когерентно суммируют, добиваясь максимальной амплитуды суммарного сигнала; решение об обнаружении возмущения принимают при превышении суммарным сигналом порогового уровня, а местоположение источника возмущения находят из решения системы уравнений, описывающих в избранной системе координат семейство сферических волновых фронтов, отстоящих друг от друга на расстояние, определяемое относительными временными сдвигами Δτ_i сигналов отдельных приемников. Достигаемым техническим результатом является повышение чувствительности обнаружения и точности определения местоположения источника возмущения. 2 ил.

Способ обнаружения ионосферного возмущения и определения местоположения его источника, основанный на анализе данных о полном электронном содержании в ионосфере Земли, которые получают в результате обработки сигналов, принятых двухчастотными приемниками спутниковой радионавигационной системы ГЛОНАСС/GPS, расположенными в узлах решетки, с последующим формированием временных рядов полного электронного содержания и их фильтрацией в диапазоне периодов колебаний, соответствующих отклику ионосферы на воздействие источника ионосферного возмущения, отличающийся тем, что используют решетку с числом n (n>3) пространственно разнесенных приемников и с апертурой, удовлетворяющей условию зоны Френеля, а сигналы, полученные после фильтрации временных рядов полного электронного содержания, когерентно суммируют с временными сдвигами Δτ_i которые обеспечивают максимальную амплитуду суммарного сигнала, решение об обнаружении ионосферного возмущения принимают при превышении суммарным сигналом порогового уровня, а координаты источника ионосферного возмущения (х_и; у_и) находят из решения системы уравнений, описывающих в избранной системе координат семейство сферических волновых фронтов, отстоящих друг от друга на расстояние, определяемое измеренными относительными временными сдвигами Δτ_i сигналов отдельных приемников.