CC BY
35
УДК 538.566
ОСНОВНЫЕ ОСОБЕННОСТИ СТРУКТУРЫ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ, СОЗДАВАЕМЫХ В НИЖНЕЙ ИОНОСФЕРЕ НИЗКОЧАСТОТНЫМИ ПРИЗЕМНЫМИ АНТЕННАМИ. РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ РАСЧЕТОВ
С.Т. Рыбачек, М.А. Белянский
Настоящая работа является продолжением статьи [1] и содержит результаты численных расчетов, а также анализ напряженностей электрического и магнитного полей в нижней области ионосферы. Исследуется влияние на структуру полей высоты расположения приемника, условий распространения, расстояния между корреспондирующими пунктами и частоты.
Ключевые слова: электромагнитные поля, анизотропная ионосфера, низкочастотные антенны.
Используемые модели и параметры задачи
Расчеты были выполнены в диапазоне частот /=0,5-15 кГц для высот расположения приемника в ионосфере Н = 50 - 500 км и ряда трасс, отличающихся протяженностью вдоль поверхности Земли К = а6, удельной проводимостью подстилающей поверхности ст, а также параметрами ионосферы. Компоненты геомагнитного поля задавались равными Н^щ = 0,148 Э, Н0ф = 0, Н^)Г = -0,473 Э. Для
описания свойств ионосферы были использованы летние модели электронной концентрации N [2] для дневных (/ = 13 ч) и ночных (/ = 1 ч) условий распространения, а также модели эффективных частот соударений Vе из работ [3] на высотах до 120 км и [4] - от 120 до 500 км. Результаты приводятся в виде модулей компонент напряженностей электрического и магнитного полей Е^ и (в мкВ/м), создаваемых наземными диполями - вертикальным (радиальным) относительно границ раздела электрическим диполем (ВЭД) и ориентированным по орту ед магнитным диполем (ГМД) с токовыми моментами
106 Ам2 и 1010 Ам2 соответственно (% = г, 6, ф). Поля диполей других типов и ориентаций, как было показано в [1], отличаются постоянными множителями от соответствующих полей вертикального электрического и горизонтального магнитного диполей.
Результаты расчетов и их интерпретация
Результаты расчетов, приведенные на рис. 1-3, выполнены для частоты 10 кГц и трассы, характеризующейся протяженностью 3300 км и проводимостью 10-1 См/м.
На рис. 1 представлены зависимости модулей компонент Е^ и вертикального электрического
диполя от высоты Н расположения точки наблюдения в нижней области дневной ионосферы. Отметим прежде всего, что зависимости модулей компонент Еф (Н) (рис. 1, а) и Нг (Н) (рис. 1, б), как и следовало ожидать, качественно совпадают, поскольку из уравнений Максвелла для ионосферы вытекает, что для излучателя любого типа на всех рассматриваемых высотах справедливы выражения для компонент полей отдельной нормальной волны
=vsEфS/ (кг), (1)
где 5 - номер нормальной волны. В сверхнизкочастотном диапазоне, где справедливо одномодовое представление поля, для 5 = 1 соотношение (1) описывает полное поле и является строгим. В диапазоне сверхдлинных волн для полного поля в рассматриваемой ситуации, полагая г и а и учитывая, что vs и ка , можно написать приближенное выражение
Нг и Еф . (2)
Расчеты по алгоритмам, описанным в [1], показали, что относительное отличие приводимых ниже значений | Еф | и | Нг | не превышает двух процентов. Из уравнений Максвелла также следует, что на
малых высотах, где свойства ионосферы близки к свойствам свободного пространства, для полного поля можно написать аналогичное (2) приближенное равенство
Ег и-Нф, (3)
которое в рассматриваемом примере приводит к качественному совпадению модулей Ег и Нф на высотах от 50 до 65 км.
Из графиков на рис. 1 также вытекает, что при расположении диполя на высотах ~ 50-70 км компоненты первичной поляризации Ег, Еб, Нф по модулю превышают значения компонент вторичной поляризации Нг, Нб, Еф, появление которых вызвано деполяризующим действием геомагнитного поля.
Это превышение обусловлено близостью свойств ионосферы к свойствам свободного пространства на указанных высотах.
102
101
80 М, км
80 Н,кп
80 Н, км
г
60 80 И, км а б в
Рис. 1. Высотные зависимости | Е^ | (а, в) и \ | (б, г) в интервале Н = 52 -100 км:
а, б - ВЭД, 10 кГц, день; в, г - ГМД, 10 кГц, день
В окрестности Н « 75 км модули касательных компонент электрического поля (рис. 1, а), как и модули касательных компонент магнитного поля (рис. 1, б), начинают сближаться. С увеличением высоты эти компоненты совпадают, а зависимости | Ег (Н) | и | EQ д (Н) | отличаются постоянным множителем. Этим особенностям поведения полей можно дать достаточно удовлетворительную интерпретацию, если предположить, что в области высот Н > 75 км поле описывается одной собственной ионосферной волной [1]. Поля такой волны подчиняются уравнению, по форме совпадающему с уравнением, описывающим распространение плоских волн в однородной магнитоактивной среде. Это позволяет ожидать, что на больших удалениях от приземного излучателя, создающего сферическую волну, свойства ионосферных полей в нашей задаче и задаче распространения плоских волн будут близки.
Действительно, в области применимости квазипродольной аппроксимации [5], определяемой не-2
равенством | Ут | << 12Уь (и - X) |, связь между компонентами полей плоских волн в принятых в магни-тоионной теории обозначениях имеет следующий вид:
V Ев = - Не = +-, Ее
п2 = 1 Х , X =
= +-
УТХ
е2 N
(и+Ть )(и - X)
V
Нд
Е
(4)
2
и = 1 +.
ю
и + УЬ Ш&оЮ'
Здесь е, ш - заряд и масса электрона соответственно; п - индексы рефракции собственных ионосферных волн; Уь и Ут - продольная (радиальная) и поперечная составляющие вектора У, направленного противоположно вектору геомагнитного поля Но. Так как в северном полушарии для распространяющейся вверх собственной волны компонента Ноь = Ног отрицательна, то УЬ является положительное числом,
Уь = 1 еНо1 |. шю
В рассматриваемом диапазоне частот в области высот Н « 9о-5оо км днем и «19о-5оо км ночью свойства ионосферы таковы, что имеют место неравенства | и | << | << X, которые приводят к следующим отношениям компонент полей:
Ев/Ег = -Уь/Ут = -Ноь/Нот , Нд/Ее = п, п2 = + X¡Гь . (5)
Как известно [5], в задаче распространения плоских волн в однородной анизотропной среде продольная компонента магнитного поля равна нулю. В рассматриваемой авторами задаче распространения сферических волн, создаваемых приземными диполями в ионосфере, соответствующая радиальная ком-
= п
понента магнитного поля отлична от нуля. Найдем поэтому, используя (2), (4) и (5), отношения компонент магнитных полей, соответствующие сверхдлинноволновому и сверхнизкочастотному диапазонам
Я0/ Яг = -п , Я0/ Яг = - кгп/V . (6)
2
Из формулы для п в (5) следует, что значения индексов рефракции п могут быть либо вещественными, либо мнимыми. Вследствие этого, как уже отмечалось в работе [1], из двух распространяющихся вверх собственных волн существенно распространяющийся характер имеет лишь волна с вещественным (и положительным для нашей зависимости от времени) индексом рефракции, равным
|п| =у/ Х/\¥ь | = еу1\в\М/ (ю | И ог |) , \Н0/ Яг\* 5,3^ М/ (/\И 0г |) . (7)
где с - скорость света в свободном пространстве. В формуле (7), удобной для расчетов, электронная концентрация задается количеством электронов в кубическом сантиметре, а продольная составляющая геомагнитного поля Н0г - в эрстедах.
Возвращаясь к рис. 1, отметим, что, как показывают расчеты, при И = 87 км относительное отличие касательных компонент равно « 5% . На высотах И > Н1 = 87 км это отличие с увеличением высоты уменьшается, а модуль отношения касательной и радиальной составляющих равен примерно 3,2 и совпадает в соответствии с выражением (5) с модулем отношения продольной и поперечной составляющих геомагнитного поля. Однако отношение | Яд ф/Яг | в этой области не остается постоянным и определяется, как следует из отношений (6), модулем индекса рефракции (7).
На рис. 2 приводятся зависимости | Е^ (Н) | (а) и | (И) | (б) ГМД, для которого в отличие от ВЭД первичными компонентами являются Еф, Яг и Яд. При этом для ГМД, как и для ВЭД, в рассматриваемом диапазоне высот выполняется соотношение (2), а в начале ионосферы на высотах от 50 до « 70 км - соотношение (3). На фиксированной высоте, большей Н1 = 87 км, модуль каждой компоненты поля ГМД отличается постоянным множителем от аналогичной компоненты полей ВЭД.
80 90 100 150 И, км
80 90 100
150 Н, км б
Рис. 3. Высотные зависимости | Е^ | (а) и | | (б) в интервале Н=80-220 км (ВЭД, 10 кГц, ночь)
В качестве иллюстрации поведения полей в ночное время на рис. 3 приведены зависимости | Е^ (И) | (а) и | (И) | (б) ВЭД. (В областях 80-100 км и 100-220 км, которые разделены вертикальной
штриховой линией, масштаб по оси абсцисс разный). В начале ионосферного слоя, как и в дневное время, модули компонент полей ведут себя в соответствии с поляризацией падающей на ионосферу волны. В отличие от дня, ночью область сближения модулей касательных компонент полей обладает большей по сравнению с дневной протяженностью - их относительное отличие достигает величины, не превышающей 5%, лишь на высотах И > ^ = 185 км.
На рис. 4-6 приводятся модули компонент полей Е0 и . Напомним, что на высотах И, больших 87 км днем и 185 км ночью, компонента Е0 описывает одновременно компоненту Еф, а компонента Яф описывает составляющую Яд.
Высотные зависимости модулей касательных составляющих полей Е0 (И) и Яф(И) радиального электрического диполя для частоты 10 кГц и разных удалений К = а0 приводятся на рис. 4, а, и рис. 4, б,
а
соответственно. В областях 80-100 км и 100-500 км (они разделены вертикальной штриховой линией) масштаб по оси абсцисс разный. Непрерывными линиями изображены ночные зависимости, штриховыми - дневные. Дневные зависимости модулей обеих компонент на высотах Н > 100 км, в отличие от области 80-100 км, слабо зависят от Н. Ночные зависимости на высоте Н и120 км имеют экстремумы, обусловленные максимумом профиля электронной концентрации N (Н) на этой высоте. Увеличение удаления Я на заданной высоте при Н > 100 км приводит, как правило, к уменьшению амплитуд полей. При этом увеличение Я от 1 до 2 Мм дает в рассматриваемом случае уменьшение амплитуд полей приблизительно в 1,7 раз в дневной и в 3,5 раза в ночной ионосфере.
\Е„
80 90 100 200 300 400 И, км 80 90 100 200 300 400Н, км а б
Рис. 4. Высотные зависимости | Eq | (а) и | Нф | (б) для удалений R = 1; 2 и 3,3 Мм в интервале H=80-500 км. Непрерывные линии - ночь, штриховые - день (ВЭД, 10 кГц)
Иф 102.
Н = 500 км
101-
/, кГц
. \: .......
■ ■ • * \ *—>—i—i—i—i—i - J—*—>
0,6 0,8 1¡ 5 ю f, кГц
б
Рис. 5. Частотные зависимости | Ед | (а) и | Нф | (б) для высот Н = 80, 100, 200, 500 км в диапазонах частот / = 0,5-1 кГц и 5-15 кГц (ВЭД, день) Частотные зависимости модулей касательных компонент Ед (/) и Нф (/) радиального электрического диполя для трассы протяженностью 3300 км изображены на рис. 5 для дневной и на рис. 6 для ночной ионосферы. Зависимости приводятся в диапазонах частот 0,5-1 кГц и 5-15 кГц (они разделены вертикальной штриховой линией), масштаб по оси абсцисс в этих диапазонах разный. В области частот 1-5 кГц результаты расчетов в которой не приводятся, зависимости | Ед (/) | и | Нф (/) | могут меняться
на несколько порядков, что обусловлено наличием так называемых волноводных минимумов. Дневные зависимости | Ед (/)| на всех рассматриваемых высотах имеют максимум в окрестности / и 10 кГц, а зависимости | Нф(/) | - в окрестности / и 9 кГц. В ночных условиях (рис. 6) на частотах, больших примерно 9 кГц, зависимости | Ед (/)| и | Нф (/)| начинают осциллировать вследствие включения последующих нормальных волн. На частоте, примерно равной 12 кГц, наблюдается минимум, обусловленный сменой ведущих мод.
I''el
ю3-
f, кГц
0,6 0,8
/, кГц
б
Рис. 6. Частотные зависимости | Eq | (а) и \ Нф | (б) для высот H = 80, 100, 200, 500 км в диапазонах частот f = 0,5-1 кГц; 5-15 кГц (ВЭД, ночь) Заключение
Результаты численных расчетов дают представление о структуре электромагнитных полей в ионосфере и позволяют оценить параметры наземных антенн, эффективно возбуждающих ионосферные поля.
Одной из существенных особенностей структуры полей является то, что в области высот H > И\ модуль отношения касательной и радиальной компонент магнитного поля определяется выражением (7). Расчеты показывают, что в значительной области частот и широт отношение | Hq | / | Hr | >> 1. Так, на
частоте 0,5 кГц это отношение на широтах от 80 до 10 градусов и высотах 250-300 км, соответствующих максимуму электронной концентрации, может достигать 100-300.
Это обстоятельство выходит за рамки рассматриваемой задачи и дает возможность оценить эффективность возбуждения приземного волновода горизонтальными электрическими и магнитными диполями по сравнению с соответствующими вертикальными. Что касается магнитных диполей, для которых такая эффективность заметно превышает значения, соответствующие электрическим диполям, из обобщенной теоремы взаимности следует, что отношение модулей радиальных компонент электрического поля, создаваемого на поверхности Земли горизонтальным и вертикальным магнитными диполями (hm, vm), расположенными в ионосфере, определяется такой же правой частью, что и в (7):
\Ehm \l\Elm\* 5,37 V N/ (f | H 0r |) .
С учетом приведенных выше расчетов можно сделать вывод, что, по крайней мере, при расположении диполя на высотах 250-300 км на частоте 0,5 кГц горизонтально ориентированный магнитный диполь излучает намного более эффективно, чем вертикально ориентированный магнитный диполь.
Литература
1. Рыбачек С.Т., Белянский М.А. Основные особенности структуры электромагнитных полей, создаваемых в нижней ионосфере низкочастотными приземными антеннами // Научно-технический вестник СПбГУ ИТМО. - 2011. - № 5 (75). - С. 48-53.
2. Rawer K., Bilitza D., Ramakrishnan S. Goals and status of the International Reference Ionosphere // Rev. Geophys. Space Phys. - 1978. - V. 16. - P. 177-181.
3. COSPAR International Reference Atmosphere. - Berlin, 1972. - 450 p.
4. Фаткуллин М.Н., Зеленова Т.И., Козлов В.К. и др. Эмпирические модели среднеширотной ионосферы. - М.: Наука, 1981. - 256 с.
5. Budden K.G. Radiowaves in the ionosphere. - Cambridge Univ. Press, 1961. - 542 p.
Рыбачек Светлана Тимофеевна Белянский Максим Анатольевич
Санкт-Петербургский государственный университет кандидат физ.-мат. наук, доцент, [email protected]
Санкт-Петербургский государственный университет аспирант, [email protected]
а
