Научная статья на тему 'Информационная система мониторинга ионосферы'

Информационная система мониторинга ионосферы Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
416
99
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СПУТНИКОВАЯ РАДИОНАВИГАЦИОННАЯ СИСТЕМА / ЧАСТОТНО-СЕЛЕКТИВНЫЕ ЗАМИРАНИЯ / НАВИГАЦИОННЫЙ РАДИОСИГНАЛ / МОНИТОРИНГА / ИОНОСФЕРА / SATELLITE RADIO NAVIGATIONAL SYSTEM / FREQUENCY SELECTIVE DYING DOWN / NAVIGATION RADIO SIGNAL / MONITORING / IONOSPHERE

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Катков Константин Александрович, Пашинцев Владимир Петрович, Катков Евгений Константинович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Катков Константин Александрович, Пашинцев Владимир Петрович, Катков Евгений Константинович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INFORMATION SYSTEM FOR MONITORING THE IONOSPHERE

The paper deals with the creation the information system allowing on signals of satellite radio navigational systems to carry out monitoring of an ionosphere and to issue to consumers the correcting information for carrying out navigation definitions and organization of communication are considered. Sizes which measurements allow to eliminate the uncertainty arising at strong ionospheric disturbances are determined.

Текст научной работы на тему «Информационная система мониторинга ионосферы»

УДК 621.391

ИНФОРМАЦИОННАЯ СИСТЕМА МОНИТОРИНГА ИОНОСФЕРЫ

© 2016 К.А. Катков, В.П. Пашинцев, Е.К. Катков

Северо-Кавказский федеральный университет, г. Ставрополь

Статья поступила в редакцию 27.10.2016

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

Ключевые слова: спутниковая радионавигационная система, частотно-селективные замирания, навигационный радиосигнал, мониторинга, ионосфера

Известно [1-4], что при воздействии мощным коротковолновым излучением на ионосферу происходит её радионагрев, что приводит к возникновению искусственных ионосферных образований (ИИО). Их характерной особенностью является появление интенсивных неоднородно-стей электронной концентрации (ЭК), поперечный размер которых варьируется в очень широком диапазоне: от 0,1м до нескольких сотен км. Наибольшее влияние на изменение условий распространения радиоволн через ионосферу оказывают мелкомасштабные (с размерами порядка 100-1000 м) неоднородности ЭК. Последние обуславливают проявление дифракционных свойств ионосферы и возникновение замираний принимаемых сигналов спутниковых радиосистем. В случае передачи простых сигналов замирания носят общий (релеевский или райсовский) характер, а при передаче широкополосных сигналов могут иметь частотно-селективный характер. Это связано с сужением полосы когерентности трансионосферного канала связи. Вследствие этого точность позиционирования спутниковых радионавигационных систем (СРНС), использующих широкополосные радиосигналы при ИИО будет существенно снижаться. При этом потребитель не имеет никакой информации о состоянии трансионосферного канала связи, а навигационная аппаратура потребителей (НАП) не исключает из рабочего созвездия навигационный космический аппарат (НКА), сигнал которого проходит через ИИО, так как сигнал о «нездоровье» НКА в навигационном сообщении отсутствует.

Катков Константин Александрович, кандидат технических наук, заведующий кафедрой информатики. E-mail: [email protected]

Пашинцев Владимир Петрович, доктор технических наук, профессор кафедры информационной безопасности автоматизированных систем. E-mail:

pashintsevp@mail. ru.

Катков Евгений Константинович, аспирант. E-mail: [email protected]

Возникает необходимость иметь достоверную информацию о состоянии ионосферных не-однородностей и их влиянии на характеристики проходящих через нее радиосигналов. В этом могут помочь результаты, полученные в [2]. Опираясь на них можно сказать, что при помощи одного двухчастотного приемника СРНС может проводиться обнаружение ИИО по результатам измерения интенсивности мелкомасштабных неоднородностей ЭК. После обнаружения ИИО возникает задача измерения пространственных координат ИИО, т.е. осуществление пеленгации ИИО.

Зная координаты ИИО можно определить те НРС, которые проходят через нее, рассчитать ошибку определения времени прихода максимума сигнала , от которой зависит погрешность определения псевдодальности до НКА <зв = с . Многократное повышение погрешности измерения псевдодальности даже в одной радиолинии {ав), вследствие локального характера ИИО, приводит к повышению погрешности позиционирования (□ <ув ), не отвечающей требованиям, предъявляемым к СРНС [4-6]. Для решения задачи снижения погрешностей позиционирования при сильных ионосферных возмущениях предлагается создать информационную систему мониторинга ионосферы по сигналам СРНС, конечной целью которой будет являться определение наличия и характеристик ИИО и повышение точности позиционирования потребителей СРНС, сигналы которых попадают в эту область.

Цель статьи: выработка рекомендаций по созданию информационной системы мониторинга ионосферы по сигналам СРНС.

Фазовые и временные флуктуации сигнала. Известно [5-7], что пространственные (р, г = х, у, г) изменения электронной концентрации в ионосфере (р, г) = N (г) + Ш(р, г) можно описать совокупностью изменения по

высоте (г) ее среднего (фонового) значения и пространственных флуктуаций электронной концентрации Ш(р, г) в мелкомасштабных не-однородностях ионосферы. Последние характеризуются спектральной плотностью мощности флуктуаций электронной концентрации Фт (q), которая хорошо аппроксимируется степенным законом Фш (4) ~ д~у распределения неоднород-ностей пространственного спектра д =2ж/1 с показателем степени V , где I - масштаб ионосферных неоднородностей в диапазоне 1тп < I < То от

минимального (1т1п ~ несколько метров) до максимального (То ~ несколько километров).

При возникновении в слое F ионосферы искусственных возмущений возрастает интенсивность неоднородностей р (с 10-2 до 1). При этом максимальная электронная концентрация N возрастает от 2*10и-2*1012 эл/м3 до 2*1014 эл/м3. В качестве параметра, характеризующего степень ИВИ, обычно принимается среднее квадратичное отклонение (СКО) флуктуаций электронной концентрации (ЭК) в неоднородностях

ионосферы сАN = PNm. При ИВИ эта величина может достигать значений 2*1014 эл/м3 [6, 7].

Погрешность местоопределения потребителя (Ад) прямо пропорциональна погрешности измерения псевдодальности (св = с ■ст) до каждого из НКА, которая, в свою очередь, зависит от погрешности измерения времени прихода НРС

(С). Эта погрешность будет зависеть, как от параметров возмущенной ионосферы, так и от параметров искаженного на ионосферных неоднородностях навигационного радиосигнала [8, 9, 14]. Одним из таких параметров является ошибка слежения за фазой принимаемого НРС.

В случае, когда навигационный приемник не может отследить фазу несущей частоты, радиосигнал теряется. Потеря синхронизации связана со срывами слежения системы фазовой автоподстройки частоты (ФАПЧ) [9]. Согласно [9, 10] дисперсия ошибки слежения за фазой на выходе системы ФАПЧ определяются, как сумма трех слагаемых:

_2 2 2 . 2 С =С? +Ст +Сг,ПР

(1)

/

^ =|80,8су) ЬЛэ 5ес @Г(р/2 -12)

у[яТ(р/2 -1)

(2)

где 80,8 - коэффициент с размерностью [м3/с2]; /о - частота навигационного радиосигнала; ИЭ -эквивалентная толщина ионосферы (500 км); 0 -зенитный угол НКА; р - фазовый спектральный индекс.

Значение фазового спектрального индекса 1,7<р<2 характеризует спокойную или естественно-возмущенную ионосферу. При ИВИ значение р может возрастать до значений З^р^б. Входящая в выражение (2) величина сш определяет значение СКО флуктуаций ЭК на высоте максимума ионизации ионосферы (к400 км). Для практических вычислений СКО флуктуации фазы сигнала при ионосферных возмущениях (сч,б ) использование выражения (2) затруднительно, так как остается неизвестным максимальный размер ионосферных неоднородностей (!о) и величина фазового спектрального индекса (р). В настоящее время существуют технические средства, позволяющие провести ряд измерений, устраняющих эту неопределенность. Так, аппаратура «№уА1е1 ОРБ-б» позволяет измерить значение СКО инте" Г-.Т7- " / накл\

гральной ЭК в наклонной радиолинии (с^ ) на трассе распространения сигнала. Эта величина определяется выражением [7]

[ТЙ^ Г (Р/2 -1/2)

86С0Г,Г(Р/2 -1) (3)

Зная это и подставив (3) в (2), получим выражение для расчета дисперсии флуктуации фазы принимаемого радиосигнала:

80,8ж

с/0

(4)

где С^ - дисперсия флуктуаций фазы прини-

2

маемого радиосигнала; сч>т - дисперсия ошибки

2

слежения за фазой из-за теплового шума; сч>пр -собственный шум генератора приемника (Спр ~ 0,122 рад).

Величина дисперсии флуктуаций фазы принимаемого навигационного радиосигнала, согласно [11], имеет вид:

Выражение (4) показывает, что измерения

накл

величины Сшт позволит определить значение СКО флуктуации фазы принимаемого НРС

Дисперсия теплового шума, вызванная ионосферными возмущениями будет зависеть от типа используемой системы автоматического регулирования (АРУ) в схеме ФАПЧ [11, 13] и также зависеть от величины С<рз, которая в условиях дальней зоны будет определять индекс мерцаний

1- ехр(-2ст^). Выражение для дисперсии теплового шума будет иметь вид

с2 = В

С <рТ ,„2

И

ехр(С)

[1 + (И2т )-1 ] 2И2Т [1 + (И2т )-1 ]

(5)

где Вп - ширина полосы пропускания приемника; И2 - отношение энергии принимаемого НРС к спектральной плотности мощности шума; Т -

время накопления корреляционного интервала на этапе обнаружения сигнала (20 мс). В выражении (5) также присутствует величина , кото-

накл

рая определяется путем измерения °тТ и выражением (4).

СКО флуктуации фазы сигнала () будет определять ошибку слежения за временем прихода сигнала (аг). Можно предположить, что в условиях сильных ионосферных сцинтилляций ошибка слежения за временем прихода сигнала будет представлять из себя сумму ошибки временных флуктуаций сигнала (аг3) и ошибки определения времени прихода максимума сигнала

СГт=атЗ + °гН

(6)

Проиллюстрировать выражение (6) поможет рис. 1 [6]. На рис. 1а показана аддитивная

смесь гауссовского шума пх(}) и амплитуды напряжения сигнала на выходе согласованного фильтра для случаев прихода сигналов с прямоугольными огибающими и одинаковой длительности Т и разной энергией > £г2. При наличии шумов, когда ап ф 0, момент достижения максимума выходного напряжения согласованного фильтра (/ = г +Т3 + Аг) является случайным и отличается от истинного на величину погрешности аг> 0, которая обратно пропорциональна амплитуде пика указанного напряжения. В данном случае при игюх1 > игвьа2 (ах - произвольная постоянная) погрешности измерения времени запаздывания соотносятся, как: аг1 < аг2 . Такая картина будет наблюдаться в случае, когда временные флуктуации сигнала, обусловленные сильными фазовыми сцинтилляциями, отсутствуют (а « 0).

а) б)

Рис. 1. Влияние изменения параметров выходного сигнала согласованного фильтра на рост погрешности измерения времени его запаздывания а) влияние амплитуды пика выходного сигнала; б) влияние «расплывания» и смещения отклика согласованного фильтра при ЧСЗ сигнала

В случае возникновения ИИО и сопутствующих им сильных фазовых и амплитудных сцинтилляций навигационный радиосигнал будет подвержен частотно-селективным замираниям. В этом случае будет иметь место «расплы-вание» отклика согласованного фильтра (рис. 1б), которое сопровождается уменьшением его амплитуды итъ^1ч < С/геых и «остроты» (т.е. увеличением ширины основания Тхрч > Т). Здесь К -коэффициент уменьшения максимального значения амплитуды сигнала, рч - коэффициент увеличения длительности сигнала на выходе согласованного фильтра вследствие ЧСЗ входного сигнала. Кроме этого, из-за сильных фазовых флуктуаций будет возникать временная флуктуация сигнала (ат3 ф 0), которая, согласно [14], определяется выражением

а

р5

(7)

Подставляя (4) в (7), получаем выражение, позволяющее по измерениям СКО интегральной ЭК концентрации (аАмЛ) оценить СКО временных флуктуаций НРС при ионосферных возмущениях:

40 4

' накл

7Г а кт

с/0

(8)

Для оценки ошибки определения времени прихода максимума сигнала (агН) при сильных ионосферных возмущениях и возникающих при этом ЧСЗ принимаемого НРС воспользуемся выражением, полученным в [6]

=

^ =

(1 + )

2яЛК? И2

(9)

где ЛКс - ширина спектра НРС; ЛКк - полоса когерентности трансионосферного канала связи определяемая как

Л^ =

^ 2 (1 + ¿72) '

(10)

< =

(3г2 - 3гИэ + И2) с2 8ее2 ©

192/4

(11)

^ =

1

(

2^/ ^2жИ2 ЛК

1 +

:(1 + ^72)

\3/4

(12)

Анализ выражения (12) показывает, что ошибка слежения за временем прихода сигнала в схеме выбора максимума (&г), также, как и фазовая (а<р), зависит от СКО флуктуации фазы принимаемого НРС ). В свою очередь, согласно (4) величина ^ может быть получена с помощью измерений СКО интегральной ЭК концентрации в наклонной радиолинии (&ЛаКЛ )• Анализ выражения (12) показывает, что для определения ошибки слежения за временем прихода НРС необходимо измерять ширину полосы когерентности трансионосферного канала связи (ЛКк) и СКО

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ПЭС в каждой радиолинии (0"ЛлГ ). Это, очевидно, потребует модификации существующей НАП. Такие модифицированные навигационные приемники (МНП) могут стать частью информационной системы мониторинга ионосферы.

Информационная система мониторинга ионосферы. Для точного местоопределения в условиях сильных ионосферных возмущений аппаратура потребителя должна решать следующие задачи: определять величину ПЭС ионосферы N; находить значение СКО флуктуации ЭК в

- накл

неоднородностях ионосферы вдоль трассы каждого НРС; определять величину полосы когерентности трансионосферного канала ЛКк; вычислять значение СКО измерения псевдодальности (&п) до каждого НКА. Решить все эти задачи модификацией только навигационных приемников представляется экономически невыгодным. Более перспективным будет путь создания информационной системы, способной выполнять

функции мониторинга ионосферы и выдавать управляющие сигналы на наземный сегмент СРНС, находящийся в зоне прямой видимости данной системы.

Схематично подобная система изображена на рис. 2. Система состоит из базовой станции (БС) и нескольких модифицированных стационарных навигационных приемников (МНП1-МНП4). Все МНП имеют связь с базовой станцией (БС), которая, в свою очередь, выдает управляющие сигналы и корректирующие поправки наземным потребителям СРНС. На схеме эти потребители изображены в виде отдельных автомобилей.

где г - расстояние от верхней границы ионосферы до точки приема при вертикальном распространении радиоволн (600 км); и - характерный размер ионосферных неоднородностей (400 м). С учетом (7), (9), (10) выражение (6) примет вид

Рис. 2. Система мониторинга ионосферы

Схема информационной системы мониторинга ионосферы приведена на рис. 3. Она должна включать в себя следующие блоки и элементы:

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

накл и г

мерения атТ с частотой не ниже 50 Гц.

2. Базовая станция в составе:

- 2.1 Блок вычисления фазовых ошибок.

- 2.2 Блок детектирования области повышенной ионизации.

- 2.3. Блок выдачи информации потребителям.

В качестве потребителей могут выступать другие приемники СРНС. Работа подобной системы будет заключаться в следующем.

1. Каждый МНП, интегрированный в информационную систему, принимает НРС видимых НКА. В каждой радиолинии оценивается ширина полосы когерентности трансионосферного канала (ЛКк). Это можно делать как внутри МНП, так и на БС. Далее делается вывод о наличии или отсутствии области повышенной ионизации (ОПИ) на трассе распространения сигнала.

2. При наличии ОПИ навигационная задача решается с помощью адаптивного алгоритма [5]. При необходимости возможен ввод псевдоспутника в центре Земли [15]. Полученные ошибки местоопределения (Ад) выдаются потребителям.

3. В блоке вычисления фазовых ошибок, рассчитывается ошибка слежения за фазой сигнала (ар), которая отправляется в блок выдачи информации.

4. Вычисленная в блоке вычисления фазовых ошибок величина СКО флуктуации фазы сигнала

(аре) поступает в блок вычисления временной

ошибки, где вычисляется СКО определения времени прихода сигнала (аТ).

5. В блоке детектирования ОПИ по данным, полученным из МНП, с помощью алгоритма аппроксимации определяется приблизительная форма и размеры ОПИ, через которые проходят навигационные сигналы и величина полосы когерентности на трассах распространения этих сигналов.

Рис. 3. Информационная система мониторинга ионосферы

6. Вся необходимая информация выдается потребителям СРНС, которые находятся в зоне радиовидимости БС. Пакет выходной информации включает в себя следующее: поправки к вектору пространственно временных координат Ад (могут быть использованы близко расположенными (не далее 1 км) потребителями для корректировки своего местоположения); ширина полосы когерентности трансионосферного канала А^ (позволяет сделать вывод о наличии ЧСЗ в канале связи и возможностей смены режима связи); координаты ОПИ в зоне радиовидимости потребителя Q (данная информация требует пересчета в топоцентрическую систему координат каждого потребителя). Необходима для выбора рабочего созвездия космических аппаратов для организации сеанса навигации или связи; величина ошибки слежения за фазой сигнала ар (позволит потребителю скорректировать работу системы ФАПЧ); величина ошибки слежения за временем

прихода сигнала аг. Следует отметить, что набор выходных параметров может быть расширен в зависимости от требований потребителей или по специальному запросу. Кроме этого целесообразно, чтобы потребители сообщали базовой станции свои приблизительные координаты. Это позволит самой информационной системе выбирать для каждого потребителя рекомендации на организацию сеанса связи или навигации.

Выводы: информационная система мониторинга ионосферы позволит обеспечить работоспособность наземной аппаратуры СРНС при возникновении ИИО в слое F ионосферы. Элементами такой системы будут модифицированные навигационные приемники, позволяющие определять полосу когерентности канала связи (А^) и СКО флуктуаций ЭК в ионосфере (оАаГ). Далее эти МНП интегрируются в создаваемую информационную систему мониторинга ионосферы, которая сможет выполнять следующие функции:

1. Определение наличия ИИО в ионосфере.

2. Проведение навигационных измерений в условиях искусственных ионосферных возмущений с выдачей поправок близко расположенным потребителям.

3. Проведение анализа состояния ионосферы

Предлагаемая информационная система обеспечит снижение погрешностей позиционирования потребителей СРНС (Ад) при возникновении сильных ионосферных возмущений локального характера.

СПИСОК ЛИТЕРАТУРЫ:

1. Лобанов, Б.С. Исследование возможности создания в ионосфере объемных образований, эффективно взаимодействующих с электромагнитным излучением в сверхшироком диапазоне частот // Теория и техника радиосвязи. 2009. №3. С. 16-24.

2. Пашинцев, В.П. Обнаружение искусственных ионосферных образований с помощью спутниковых радионавигационных систем / В.П. Пашинцев, С.А Коваль, В.И. Стрекозов, М.Ю. Бессмертный // Теория и техника радиосвязи. 2013. №1. С. 112-116.

3. Пашинцев, В.П. Влияние ионосферы на измерение времени запаздывания сигнала в спутниковых радионавигационных системах / В.П. Пашинцев, М.В. Гамов // Радиоэлектроника. 2002. Т. 45, №12. С. 3-13.

4. Пашинцев, В.П. Спутниковая навигация при ионосферных возмущениях / В.П. Пашинцев, К.А. Катков, Р.П. Гахов и др. - Ставрополь: СевКавГТУ, 2012. 259 с.

5. Пашинцев, В.П. Адаптивный алгоритм определения вектора пространственно-временных координат / В.П. Пашинцев, К.А Катков, Р.П. Гахов // Вопросы радиоэлектроники. 2013. Т. 4, №1. С. 138-150.

6. Пашинцев, В.П. Оценка погрешности измерения псевдодальности в спутниковых радионавигационных системах при возмущениях ионосферы в слое F / В.П. Пашинцев, М.А. Солчатов, А.М. Спирин, К.А. Катков // Физика волновых процессов и

радиотехнические системы. 2007. Т. 10, №6. С. 813.

7. Маслов, О.Н. Модели трансионосферных радиоканалов и помехоустойчивость систем космической связи / О.Н. Маслов, В.П. Пашинцев // Приложение к журналу «Инфокоммуникационные технологии». Вып. 4. - Самара: ПГАТИ, 2006. 357 с.

8. Rino, C.L. A power law phase screen model for ionospheric scintillation // 1, Weak scatter, Radio Sci. 1979. V. 14(6). P. 1135-1145.

9. Beniguel, Y. Scintillations effects on satellite to Earth links for telecommunication and navigation purposes / Y. Beniguel, B. Forte, S.M. Radicella et al. // Annals of geophysics. 2004. Supplement to vol. 47. № 2/3. P. 1179-1199.

10. Conker, R.S. Modeling the effects of ionospheric scintillation on GPS/Satellite-Based Augmentation System availability / R.S. Conker, M.B. El-Arini, C.J. Hegarty, T. Hsiao // Radio Science. 2003. Vol. 38, No. 1. P. 1-23.

11. Катков, К.А. Влияние параметров возмущенной ионосферы на ошибку слежения за фазой навигационного радиосигнала / К.А. Катков, В.П. Пашинцев, Е.К. Катков // Современная наука и инновации. Вып. №2(14). - Ставрополь - Пятигорск, 2016. С. 52-64.

12. Рытов, С.М. Введение в статистическую радиофизику. Часть 2. Случайные поля / С.М. Рытов, Ю.Н. Кравцов, В.И. Татарский. - М.: Наука, 1978. 468 с.

13. Knight, M.F. A Comparison of Predicted and Measured GPS Performance in an Ionospheric Scintillation Environment / M.F. Knight, M. Cervera, A. Finn // Proceedings of ION GPS-99, The Twelfth International Technical Meeting of the Satellite Division of the Institute of Navigation. 1999. P. 1437-1450.

14. Тепляков, И.М. Ионосферные искажения цифровых сигналов с широкополосной модуляцией // Радиотехника. 1984. № 8. С. 8-13.

15. Катков, К.А. Использование псевдоспутника в центре Земли в спутниковых радионавигационных системах / К.А. Катков, С.А. Скорынина, М.С. Окулова // Двойные технологии. 2009. №4. С. 5663.

INFORMATION SYSTEM FOR MONITORING THE IONOSPHERE

© 2016 K.A. Katkov, V.P. Pashintsev, E.K. Katkov

North Caucasus Federal University, Stavropol

The paper deals with the creation the information system allowing on signals of satellite radio navigational systems to carry out monitoring of an ionosphere and to issue to consumers the correcting information for carrying out navigation definitions and organization of communication are considered. Sizes which measurements allow to eliminate the uncertainty arising at strong ionospheric disturbances are determined.

Key words: satellite radio navigational system, frequency selective dying down, navigation radio signal, monitoring, ionosphere

Konstantin Katkov, Candidate of Technical Sciences, Head of the Informatics Department. E-mail: [email protected]; Vladimir Pashintsev, Doctor of Technical Sciences, Professor at the Automated Information Security Systems Department. Email: [email protected]; Evgeniy Katkov, Post-graduate Student. E-mail: [email protected]

i Надоели баннеры? Вы всегда можете отключить рекламу.