Содержание
Введение 5
Тема исследования и обзор литературы...................................5
Цель диссертационной работы...........................................25
Публикация материалов диссертации и её краткое содержание .. 25 Положения, выносимые на защиту...................................34
1. Распространение волн в приземном волноводе с локальной трехмерной неоднородностью на ионосферной стенке (метод последовательных приближений) 37
1.1. Импедансная модель локально нерегулярного волноводного канала Земля-ионосфера...........................................38
1.2. Задача о поле вертикального электрического диполя в импедансиом волноводе с локальной неоднородностью................46
1.3. Задача о поле вертикального магнитного диполя в импедансном волноводе с локальной неоднородностью............................68
1.4. О применимости скалярного приближения в задаче о распространении радиоволн в локально-нерегулярном волноводе 79
1.5. Задача о поле вертикального электрического диполя в волноводе с локальной неоднородностью, нарушающей геометрическую регулярность ионосферной стенки волноводного канала...................89
1.6. О границах применимости метода последовательных приближений в задаче о распространении радиоволн в локально нерегулярном волноводе...............................................109
1.7. Сравнение расчета на основе сформулированной теории с результатами эксперимента по искусственному воздействию на ионосферную плазму...................................................113
2
2. Асимптотический метод решения задачи распространении
радиоволн в трехмерно неоднородном волноводном канале
(горизонтальные размеры неоднородности произвольны) 119
2.1. Асимптотический метод решения локально нерегулярных волноводных задач.............................................120
2.2. Новый подход к классической задаче о береговой рефракции (прямолинейная линия берега)..................................134
2.3. Задача о береговой рефракции в импедансном волноводе (геометрия береговой линии - произвольна).....................148
3. Распространение волн в приземном волноводе с локальной
трехмерной крупномасштабной неоднородностью (скалярная постановка задачи) 158
3.1. Распространение радиоволн в плоской модели волноводного канала с локальной неоднородностью, лежащей в плоскости одной из стенок волновода...................................................159
3.2. Распространение радиоволн в сферической модели волноводного канала с локальной неоднородностью, не выступающей за пределы стенки волновода...................................................174
3.3. Распространение радиоволн в плоской модели волноводного канала с локальной неоднородностью, в виде усеченного цилиндра 185
3.4. Распространение радиоволн в плоской модели волноводного канала с локальной неоднородностью, имеющей произвольную гладкую поверхность................................................197
3
4. Распространение электромагнитных волн в приземном волноводе с локальной трехмерной крупномасштабной неоднородностью нижней ионосферы (векторная задача) 209
4.1. Возможный подход к решению векторной трехмерной локальнонерегулярной волноводной задачи..............................210
4.2. Деполяризация электромагнитного поля при рассеянии на локальной крупномасштабной неоднородности нижней ионосферы (модельная неоднородность в плоскости стенки волновода)...........229
4.3. Исследование рассеяния электромагнитного поля на ионосферной неоднородности, имеющей вертикальное измерение...............242
Заключение 257
Список литературы 262
4
Введение
Тема исследования и обзор литературы
Диссертация посвящена теоретическим аспектам проблемы распространения радиоволн в околоземном волноводном канале с учетом трехмерных геометрических и электрических особенностей нижней ионосферы и земной поверхности. Основное внимание уделено разработке математических методов моделирования и расчета волновых полей, возбуждаемых гармоническими во времени точечными источниками в слоистых средах с трехмерными локальными неоднородностями. Отметим, что математические аспекты названной физической проблемы включают в себя выбор и построение модели среды распространения, корректную постановку математической задачи, эффективное решение возникающих уравнений и исследование получаемых решений с целыо выяснения зависимостей характеристик электромагнитного поля от величины различных параметров задачи. Под термином "локальная неоднородность" понимается нерегулярность среды распространения, имеющая конечные размеры по всем направлениям, и что особенно важно для задач распространения, в направлении источник-точка наблюдения и в поперечном к нему (внутри рассматриваемого слоя) направлении. При этом мы не ограничиваем себя неоднородностями, размеры которых малы в масштабе какого-либо параметра задачи (например, длины волны исследуемого электромагнитного поля, длины трассы распространения, размера главной зоны Френеля). Учет влияния такой неоднородности приводит к задачам математической физики с не разделяющимися переменными, и исключает возможность сведения задачи к двумерной. Надо также отметить, что речь пойдет о задачах распространения в среде с детерминированными, а не статистическими неоднородностями, и поэтому объектами исследования будут являться
5
непосредственно амплитудно-фазовые характеристики электромагнитного поля, а не его статистические моменты.
К рассматриваемому типу неоднородностей, кроме неоднородностей ионосферы и земной поверхности в задачах околоземного распространения радиоволн, могут быть также отнесены трехмерные возмущения диэлектрических и металлических волноведущих систем, трехмерные включения в слоистых средах, моделирующие нерегулярности в задачах радиолокации и дистанционного зондирования. В различных отраслях современной техники, в частности, в электромагнитной дефектоскопии, разведочной геофизике и других приложениях важное значение имеют исследования возмущений электромагнитного поля, обусловленные наличием рассеивателей различной физической природы, и на этой основе определение информационных параметров этого поля.
Из столь широкого списка возможных приложений развиваемой теории, в качестве основного для данной диссертации, будем иметь в виду задачу околоземного распространения радиоволн низкочастотной части спектра (СДВ и ДВ диапазонов), что будет определять параметры рассматриваемых в работе моделей среды распространения. Важно также заметить, что в указанном диапазоне частот свойства ионосферной плазмы таковы, что оказывается возможным использование импсдансной модели приземного волновода. И последнее замечание к определению круга рассматриваемых в диссертации задач. Внимательное изучение, в том числе и геофизической литературы, показывает действительное наличие потребности понимания поведения радиоволн данного диапазона в волноводном канале Земля-ионосфера, поскольку именно электромагнитные поля в этом диапазоне частот несут огромный объем информации об окружающей среде и, в частности, о естественных и искусственных возмущениях ионосферы, вплоть до возмущений, связанных с лунными полусуточными приливами и еженедельными вариациями в линиях электропередачи [1]. Решение прямых трехмерных задач должно в дальнейшем стать основой для перехода к решению обратной задачи распознавания неоднородностей нижней ионосферы на основе результатов наблюдений амплитуды и фазы как естественных, так и искусственных электромагнитных полей.
6
В настоящее время не существует канонических универсальных методов, которые без сомнения приводили бы к решению трехмерной задачи распространения волн с заданной точностью. И даже бурное развитие вычислительной техники и программной реализации численных алгоритмов не позволяет уповать в данной ситуации на всесилие численных методов. Во-первых, система уравнений Максвелла является системой эллиптического типа [2], для которой известно, что нельзя корректно поставить задачу Коши. Весьма перспективными являются методы, заменяющие эллиптическую систему уравнений на близкую к ней параболическую систему, для которой корректна постановка начальной задачи. Однако, такая замена вносит принципиально не уменьшаемую погрешность, оценить которую достаточно сложно [3]. Во-вторых, дискретизация исходной математической задачи накладывает жесткие требования на объем вычислительных ресурсов (объем оперативной памяти и быстродействие) используемой вычислительной техники, что заставляет ограничиваться в задачах распространения рассеивающими объектами с размерами не более нескольких десятков длин волн. В таких условиях говорить о решении обратной задачи в реальном масштабе времени численными методами - пока не приходится. Это заставляет искать новые подходы на пути синтеза численных и аналитических методов. Последние, в противоположность чисто численным методам, позволяют выявлять основные физические свойства волновых полей в разнообразных ситуациях, указывать простые качественные взаимосвязи между параметрами задачи, что позволяет, используя эту информацию, предсказывать особенности волновых нолей для реальных структур, в ситуациях, когда использование численных методов или нерационально или даже невозможно (например, если требуется указать "разумное" начальное приближение).
В то время как теорию распространения радиоволн в околоземном пространстве (волноводный канал Земля-ионосфера, океанический волновод, ионосфера, магнитосфера Земли) в предположении его однородности и поверхностной регулярности границ можно считать в достаточной степени завершенной [4-14], то уже переход к моделям, учитывающим реальные условия распространения, даже в простейшем виде, когда возможно свести векторные
7
уравнения Максвелла к скалярным уравнениям Гельмгольца, и затем осуществить процедуру разделения переменных, продолжает оставаться достаточно бурно развивающейся областью математической физики, в которой современные исследователи продолжают находить ответы на вопросы, возникающие в различных приложениях. Именно поэтому, описанные в отмеченных монографиях, а также в книгах [15-20], точные аналитические решения задач распространения волн в структурах эталонной геометрической формы, таких как плоскослоистые среды, цилиндрические и сферические области, клин и некоторых других, а также аналитические методы исследования волновых полей, опирающиеся на медленность изменения одних параметров по отношению к другим [21, 22], позволяют получать явные формулы для самих полей либо их характеристик, а также представляют возможность делать качественные заключения о физических свойствах уже нерегулярных структур. Формулы для одиннадцати известных систем координат, допускающих разделение переменных в скалярном уравнении Гельмгольца, к которому сводится большинство задач теории распространения электромагнитных волн в скалярной постановке, приводятся в [23]. При этом отдельные решения для канонических геометрий могут быть использованы для проверки трехмерного решения в предельных случаях, а также в качестве функций Грина в методе интегральных уравнений.
И все же, в большинстве практически значимых случаев, ког да качественных выводов оказывается недостаточно, необходимо искать новые подходы, опираясь на фундамент знаний, заложенных предшественниками. Остановимся более подробно на методах и подходах, которые были непосредственной отправной точкой исследований, представленных в данггой диссертации. В большом количестве работ, затрагивающих проблему оценки влияния на электромагнитное тюле трехмерных неоднородностей среды распространения, и подробный обзор которых будет проведен ниже, основным методом построения решения является метод поверхностных интегральных уравнений или метод теории потенциала. Основное преимущество этого метода состоит в том, что он сводит задачи для уравнений в частных производных в неограниченной трехмерной области к интегральным
8
уравнениям на двумерной поверхности и позволяет автоматически удовлетворять условиям на бесконечности [24, 25]. Большое внимание этому методу уделено в монографиях [10, 20, 26, 27], в работах [28, 29, 30]. Надо признать, что метод интегральных уравнений играет одну из центральных ролей в изучении граничных задач, связанных с дифракцией электромагнитных и акустических волн на ограниченных телах. Это обусловлено главным образом тем, что математическая постановка таких задач приводит к уравнениям в неограниченных областях, и, следовательно, их переформулировка в виде граничных интегральных уравнений не только уменьшает размерность задачи, но и позволяет свести ее к задаче в ограниченной области (по границе рассеивателя). При этом не накладывается ограничений на форму (поверхность рассеивателя должна, конечно, оставаться в классе поверхностей типа Ляпунова) и количество рассеивателей. Именно поэтому методы поверхностных, а также объемных интегральных уравнений являются отправной точкой многих, развиваемых на их основе, численных алгоритмов решения задач дифракции [31, 32, 33, 34, 35, 36]. Этот список ссылок может быть неограниченно продолжен, т.к. литература по численным методам решения задач дифракции чрезвычайно обширна и постоянно пополняется. Укажем также на возможность использования при решении задач дифракции методов, основанных на математической теории сингулярных интегральных уравнений [37, 38]. Эти последние ссылки я привожу только для того, чтобы показать, что выбор пути, по которому мы пошли в данной диссертации, был сделан на основе анализа достаточно широкого спектра имевшихся возможностей.
Характерной особенностью работ [24, 28, 29, 30] является то, что все теоремы касающиеся интегральных уравнений в теории дифракции и рассеяния, проводятся их авторами для линейных операторов, задаваемых на множестве функций из пространства С, а не Ь 2, как это принято делать в чисто математической литературе [39, 40, 41]. В данном случае С(С) - банахово пространство комилекснозначных непрерывных функций, определенных на й и оснащенное нормой ||ф|| = тал:|ф(д:)|.
хеС
А пространство Ь2- пространство комплекснозначных интегрируемых с квадратом
9
функций. Разрабатываемые в отмеченных работах методы применимы, в основном, для задач рассеяния в случае |&г| < 1 в, так называемом, длинноволновом приближении.
Отметим важное, на наш взгляд, отличие двумерных и трехмерных задач дифракции, возникающее при решении их методом поверхностных интегральных уравнений. При решении задач дифракции в рамках этого метода приходится работать с двумя формами математических выражений: это, так называемое, интегральное соотношение и собственно интегральное уравнение. Интегральное соотношение, получаемое с использованием теоремы Грина (ее скалярного варианта [42], либо векторного [24, 43], называемого часто формулой Стрэтгона-Чу), выражает значение искомого поля в любой точке внутри волноводного объема через значения этого же поля и его нормальной производной на поверхности, шраничивающей этот объем. Часто это соотношение называют формулой Кирхгофа (или интегралом Кирхгофа) [10], или интегралом Гельмгольца [20], иногда - интегралом Гюйгенса, поскольку она служит математическим определением принципа Гюйгенса. Чтобы получить интегральное уравнение необходимо осуществить предельный переход точки наблюдения из внутренней точки волноводного объема на его граничную поверхность, которая является, в тоже самое время, и поверхностью интегрирования. Этот предельный переход является ключевым моментом теории потенциала. В результате него выделяется скачок нормальной производной функции искомого поля, и интегральное соотношение, тем самым, преобразуется в интегральное уравнение, в котором области определения функций, стоящих по разные стороны от знака равенства, совпадают. После такого предельного перехода поверхностный интеграл в уравнении уже должен пониматься в смысле главного значения (У.Р.), с "выколотой" точкой, соответствующей точке наблюдения, опущенной на поверхность. В тоже время, чисто формально можно показать [10, 20], что подынтегральное выражение не всегда имеет особенность в этой точке, которая определяется членом с нормальной производной функции Грина. Эга нормальная производная в двумерном случае (дифракция плоской волны на цилиндрическом объекте) определяется радиусом кривизны цилиндрической поверхности в данной предельной точке. В трехмерном
10
случае, также имеется конечное значение этого предела, которое будет определяться уже двумя главными радиусами кривизны отражающей поверхности. Но, если в двумерном случае значение этого конечного предела единственно, то в трехмерном -оно зависит от направления, по которому точка движется к своему предельному значению. А это значит, что функция (нормальная производная функции Грина) в этой точке предела не имеет [44]. Высказанное замечание может оказаться важным при необходимости численной оценки интеграла в формуле Кирхгофа, например, при использовании метода итераций, либо метода Кирхгофа.
Остановимся подробнее на одной из классических задач теории распространения волн, задаче, возникшей в связи с проблемой береговой рефракции радиоволн, задаче, для решения которой был использован метод интегральных уравнений, и аналог которой описан в данной диссертации в качестве примера использования разработанной нами методики. В теории дифракции под задачей о береговой рефракции понимают задачи дифракции на бесконечной плоскости, поверхностный импеданс которой иретерпеваег скачок вдоль некоторой прямой (т.н. задача об имледансной ступеньке). В теории распространения это будет задача о распространении земной волны над двухкусочной трассой. Заметим, что теория дифракции рассматривает эту задачу как частный случай дифракции на клине с различными импедансами сторон и углом раскрыва, равным л. По-видимому, впервые задача о набегании радиоволн на прямолинейный берег моря была сформулирована в [45] с помощью приближенных граничных условий импедансного типа, записанных относительно вертикальной компоненты электрического поля (при этом проводимость моря считалась бесконечно большой). Задача была записана как в виде уравнения Вольтерра, так и в виде неоднородного интегрального уравнения второго рода с ядром, зависящим только от модуля разности аргументов, которое затем было решено приближенно [46]. Как было выяснено, неоднородное интегральное уравнение второго рода с ядром, зависящим только от модуля разности аргументов, принадлежит к классу интегральных уравнений, для которых возможно точное решение. Теория однородного уравнения с разностным ядром была развита в [47], хотя впервые неоднородное уравнение типа Винера-Хопфа было рассмотрено (без
11
вопросов единственности и существования) в [48], где и было представлено его частное решение. Полное математическое исследование этого уравнения было представлено в [49, 50], что позволило получить точное аналитическое решение задачи для случая плоской волны, распространяющейся нормально к прямолинейной линии берега [51, 52]. Отметим, что метод Винера-Хопфа-Фока это один из старейших методов получения аналитического решения в задачах дифракции, он был использован также в [53, 54]. К сожалению, точное решение, получаемое этим методом, записывается в виде некоторых интегралов, форма которых до настоящего времени позволяет получать непосредственные численные результаты только для ограниченного числа специальных случаев, в которых замкнутая форма исходных интегралов допускает возможность преобразований. В частности, в [51, 52] был исследован приближенный вид решения, применимый в случае не слишком косого падения волны на линию берега. Аналогичная задача о дифракции плоской Н-поляризованной электромагнитной волны, набегающей перпендикулярно к береговой линии, в предположении об идеальной проводимости моря была решена независимо в работе [55]. Обобщение на случай конечных значений импедансов обеих полуплоскостей было выполнено в [56, 57, 58]. В акустике проблема расчета звукового давления, возбуждаемого гармоническим источником над неоднородной плоскостью, также является классической и рассматривалась неоднократно [59, 60], а также библиография к [60]. В исследовании [61] путем интегрирования решений для плоских волн была решена задача дифракции цилиндрической волны. Дифракция поверхностных волн при условии их нормального падения на линию скачка поверхностного импеданса стала предметом работ [62, 63]. Решение задачи о возбуждении поверхностных волн нитевидным источником (в общем случае мультипольным), расположенным на линии скачка импеданса, представлено в работе [64]. Работа [65] посвящена дифракции сферической звуковой волны. В случае нормального падения к границе раздела исследование скалярных (акустических) и векторных (электромагнитных) полей проводится практически одинаково. При более общей постановке задачи, когда допускается косое падение волны, такое вырождение снимается, и каждый из перечисленных случаев требует отдельного рассмотрения.
12
Задача о наклонном падении на импедансную ступеньку решена методом Винера-Хопфа-Фока в статье [66] для поверхностной звуковой волны и в статье [67] для поверхностной электромагнитной волны.
В теории распространения радиоволн над неоднородной в электрическом отношении земной поверхностью отметим отдельно работы [68, 69, 70, 71], посвященные задаче о двухкусочной трассе, когда основное внимание уделяется влиянию на поле именно скачка свойств поверхности, над которой распространяется либо земная волна, либо исследуется поле точечного источника. Авторы [68, 69, 70] в двумерном варианте рассматривают распространение над сферой. При этом в [68] наряду с методом интегральных уравнений, когда решение получается на первом шаге метода итераций, в том числе и для случая распространения вдоль береговой линии, используется и метод параболического уравнения. В работах [69, 70] ищется решение для точечного источника. Решение задачи о распространении земной волны, падающей под произвольным углом на линию берега рассмотрено в статье [71] методом интегрального уравнения Вольтерра. Здесь же стоит обратить внимание на работы, в которых рассматривается проблема перевозбуждения мод на границе суша-море в волноводном канале Земля-ионосфера [72, 73]. Примером того, что интерес к теории береговой рефракции не ослабевает и, что задача представляет интерес в связи с разнообразными приложениями, может служить статья [74], которая, в отличие от уже упомянутых работ, посвящена построению приближенной теории для случая, когда берег может быть криволинейным. В [74] в скалярном случае, путем приближенного разделения переменных в уравнении Гельмгольца, решение задачи находится в виде отклонения от решения В.А.Фока для прямолинейного берега. В работах [75, 76] рассмотрен аналог классической задачи о береговой рефракции -трехмерная задача о поле вертикального электрического диполя в нерегулярном импедансном волноводе. Рассмотрены плоская и сферическая модели волноводного каната. Для прямолинейной границы раздела сред получено решение в случае трассы, проходящей вдоль такой границы. Рассмотрен случай произвольной геометрии береговой линии. И последнее, результаты, полученные во всех только что отмеченных работах, имеют смысл только в пределах применимости приближенных
13
импедансных граничных условий [77], которые могут быть определены только в случае, если известно точное решение соответствующей задачи дифракции. В частности, для того случая, когда граница раздела сред и трасса распространения пересекаются под прямым углом, можно воспользоваться результатами работ [78, 79, 80, 81, 82]. В работах [78, 79, 80] представлены результаты строгого исследования двумерных задач дифракции электромагнитных волн в системе из двух прямоугольных клиньев, имеющих общую грань, из которых один предполагается диэлектрическим, а другой - идеально проводящим. Решение было построено в случае |е;п|»1, где е'т - относительная диэлектрическая проницаемость клина. Из
результатов [79, 80] следует, что различия в импедансном и точном решениях задачи проявляются только в непосредственной близости от ребра клиновидной структуры на расстоянии порядка длины волны в диэлектрическом клине. Более общий случай, допускающий конечную проводимость обоих прямоугольных клиньев, был рассмотрен в работах [81, 82]. Задача дифракции была сведена к системе из двух сингулярных интегральных уравнений для преобразованного по Фурье граничного значения магнитного поля, рассеянного в однородное пространство, для которой была установлена однозначная разрешимость и описаны схемы приближенной регуляризации в пространствах Ь2 и С0. Приближенное решение было построено и исследовано численно и аналитически как в дальней зоне, так и вблизи ребра системы, только в случае |е[|»1 и |е'2|»1 - больших но абсолютной величине
относительных комплексных диэлектрических проницаемостях клиньев. Оказалось, что импсдансное решение хорошо аппроксимирует точное решение практически всюду - вплоть до экспоненциально малых расстояний от ребра. Результаты [78, 79, 80, 81, 82] доказывают правомочность применимости импедансной постановки задачи (для случая нормального падения плоской волны на грань клина) всюду за исключением экспоненциально малой окрестности ребра системы. Чтобы закончить с этим вопросом, необходимо остановиться на работах [52, 83], где было проведено уточнение теории скин-эффекта и граничных условий Леонтовича для выпуклых тел, в частности, в условиях распространения радиоволн над земной поверхностью, а
14
также указать на работу [84], обобщающую известные сведения о граничных условиях ЦЦукина-Леонтовича, в том числе и для структур с ребрами.
Теперь будет логично, чтобы закончить с обзором двумерных задач теории распространения радиоволн, остановиться на работах, посвященных анализу электромагнитного поля на многокусочных трассах без учета их поперечной структуры. В англоязычной литературе это приближение получило название mixed-path approximation. Характерная черта всех исследований - это использование приближенных импедансных граничных условий на поверхности, над которой происходит распространение волн. Большинство цитируемых работ рассматривает проблем}' в рамках метода поверхностных интегральных уравнений. Поскольку поперечная структура трассы не учитывается, в исходном интегральном соотношении Кирхгофа проводится приближенное интегрирование по поперечной координате, и задача переформулируется в виде интеграла Вольтерра, для решения которого предлагаются различные схемы. Если отклонения геометрии поверхности от идеальной (рельеф) не учитываются [85, 86, 87, 88, 89], то возможны аналитические представления решения. Однако, такие формулы, удобные для трасс с протяженными участками различных импедансов, ориентированы на использование для решения прикладных задач. В случае большого числа непротяженных участков предпочтительнее будут численные методы решения. В статье [90] рассмотрена ситуация, когда зависимость импеданса от расстояния от излучателя представима в виде ряда по полу целым степеням. Здесь же отметим работы [91, 92], посвященные трассовому приближению в задаче о распространения волн в сферическом анизотропном плавно-нерегулярном волноводе Земля-ионосфера, когда не учитывается влияние на фазу собственной волны рефракционных эффектов, вызванных поперечной нерегулярностью волновода. Приближение построенное в [91] позволяет оценить погрешность трассового приближения и при заданной точности расчета фазы определить расстояния, на которых возможно использование трассового приближения. В работе [92] анализируется влияние поперечной к направлению распространения нерегулярности свойств волноводного канала Земля-ионосфера на электромагнитное поле низкочастотных излучателей, расположенных в нижней
15
ионосфере и на спутниковых высотах. Приводятся численные оценки влияния поперечной нерегулярности на фазу и амплитуду поля при распространении в окрестности экватора.
Остановимся на работах исследующих, одномерные трассы, но с учетом высотного рельефа земной поверхности. Два альтернативных интегральных уравнения типа Вольтерра получают авторы [93] и [94, 95], учитывая высотный профиль трассы в области между источником и приемником. Влияние резко выраженной геометрической неоднородности трассы распространения исследовалось в работах [96, 97] в случае, когда рельеф трассы задавался в виде гауссовой кривой. В случае произвольной функции рельефа чаще применяется одномерное интегральное уравнение [93]. В ряде работ численные расчеты на основе решения интегрального уравнения Вольтерра сопоставляются с результатами экспериментов [98, 99, 100, 101]. При этом в [100] задача усложнена исследованиями распространения импульсного сигнала, [101] переносит исследования в область более высоких частот. Модификация известных уравнений производится в [102], где различаются случаи выпуклой и вогнутой поверхностей, над которыми происходит распространение радиоволн. Аналогичная акустическая задача рассматривается в [103, 104]. Метод
параболического уравнения, родственный методу интегрального уравнения Вольтерра, применяется в статьях [105, 106]. Оба эти метода пренебрегают отраженными на линии скачка импеданса волнами. Обзор известных в литературе примеров использования параболического уравнения в теории распространения волн приводится в [107].
Обзор трехмерных задач теории распространения радиоволн мы ограничим задачами распространения для различных моделей, учитывающих трехмерные локальные неоднородности волноводного канала Земля-ионосфера. Таким образом, за рамками данного обзора будут оставлены исследования задач дифракции и рассеяния на трехмерных включениях в слоистых средах, в диэлектрических и металлических волноведущих системах, обычно проводимые на основе объемных интегральных уравнений и прямых численных методов. Прежде чем описывать, имеющиеся к настоящему времени в литературе примеры решения задач математической физики,
16
относящихся к заявленной проблеме, кратко опишем круг природных явлений, приводящих к нарушению однородности параметров ионосферной плазмы, т.е. ведущих к образованию локальных неоднородностей нижней ионосферы. Здесь необходимо отметить возмущения, обязанные своим происхождением терминатору, землетрясениям, магнитным бурям, образованию спорадического Е-слоя [108], высыпаниям электронов и протонов из радиационных поясов. Размеры таких неоднородностей в горизонтальной плоскости могут быть от десятков и сотен до тысяч километров в диаметре. Обе полярные области ионосферы могут рассматриваться как локальные трехмерные неоднородности на длинных трассах распространения радиоволн. Отметим отдельно, что причиной широко известного явления возмущения амплитуды и фазы регистрируемых СДВ и ДВ сигналов [109, 110, 111], называемого Trimpi эффект, также являются локальные изменения нижней ионосферы. Помимо возмущений естественного характера возможно появление локальных ионосферных неоднородностей, являющихся следствием деятельности человека. Причем такую деятельность можно подразделить на два вида. Во-первых, это события, для которых модификация параметров ионосферы оказывается побочным результатом. К таким можно отнести: старты космических аппаратов [112, 113, 114, 115], промышленные взрывы [116], воздействие мощных наземных источников электромагнитного поля [117], возможное появление в будущем спутниковых энергетических установок [118, 119]. И, во-вторых, это преднамеренное искусственное возмущение плазменной среды в рамках стремительно развивающегося в настоящее время научного направления в исследованиях околоземной среды, так называемыми, активными методами [120]. Это эксперименты, сопровождающиеся впрыскиванием химически активных веществ [121], и эксперименты по направленному облучению ионосферной плазмы мощным высокочастотным электромагнитным полем (ВЧ нагрев ионосферы) [122, 123, 124, 125, 126]. В большой своей части такие исследования носят ярко выраженный экологический характер. Размеры рукотворных локальных неоднородностей конечно уступают размерам естественных неоднородностей, но, в конечном итоге, любой размер важен в сравнении с двумя величинами: либо с длиной волны
17
распространяющегося излучения, либо с длиной трассы распространения. Локальные неоднородности на земной стенке волновода - это неоднородности "островного" и "полуостровного" типа. Их ограниченность в направлении, поперечном к трассе распространения, может существенно влиять как на амплитуду, так и на фазу распространяющегося над такими объектами сигнала, особенно в случае, когда линия трассы проходит вблизи касательного направления к границе неоднородности [127].
Большинство теоретических работ, посвященных анализу трехмерных задач теории распространения, ставило своей целью объяснить наблюдавшиеся на эксперименте вариации амплитуды и фазы, в основном СДВ и СНЧ, радиосигналов рассеянием электромагнитных волн на локальных неоднородностях ионосферы. Предполагая такие неоднородности малыми возмущениями задачи, решение определялось в виде поправок к невозмущенному решению в регулярном волноводе. Из наиболее ранних работ отметим [128], где на основе теории дифракции Френеля производятся оценки влияния на одномодовое СДВ иоле в плоском волноводе малого локального понижения высоты ионосферы. В работах [129, 130, 131, 132, 134] решение получено для неоднородностей малых в масштабе длины волны, распространяющегося поля, а в [135, 136] для неоднородностей малых в масштабе длины трассы распространения. Авторы упомянутых работ используют метод интегральных уравнений. Для описания ионосферной неоднородности привлекаются различные математические модели. Наиболее часто используется идея, высказанная еще в [129], где предлагалось характеризовать неоднородность среды распространения малыми возмущениями постоянных распространения волноводных мод, которые таким образом становились двумерными функциями координат в горизонтальной плоскости распространения волны, аналогично индексу рефракции в пространственно неоднородной среде. Взаимодействием мод пренебрегалось, и в конечном результате рассматривалась только одна основная нормальная волна. Полученное в [129] интегральное уравнение неоднократно использовалось в дальнейшем [130, 131, 136], но только в предположении одномодовости сигнала, что вынуждало авторов ограничиваться диапазоном СНЧ, либо СДВ на очень
18
протяженных трассах. В качестве методов решения двумерного интегрального уравнения использовались либо прямые численные методы [131], либо в качестве решения рассматривался первый шаг метода последовательных приближений (Борновское приближение) [129, 132, 134, 136]. В [132, 134] неоднородность задается возмущением эффективной высоты волноводного промежутка, которое, опять же только для основной моды исследуемого СНЧ поля, через известное дисперсионное соотношение можно пересчитать в возмущение постоянной распространения. Поскольку в области, занятой неоднородностями тина Е$ или вызванными высыпанием электронов, затухание волн СНЧ может возрастать в несколько раз, то изменения такой “постоянной” распространения СНЧ в волноводе Земля-ионосфера несомненно есть. Однако, на наш взгляд, выглядит спорной возможность рассматривать в такой модели возмущение, располагающееся непосредственно над источником (или над приемником), поскольку в непосредственной близости от источника любое, даже самое низкочастотное поле, описывается бесконечной суммой нормальных мод, выделить в которой главную (которая определит постоянную распространения) до некоторых расстояний от источника оказывается невозможным. Это подтверждают результаты работы [132], где в отличие от [133] показана необходимость учета ближних полей (электростатической и индукционной составляющих) в задаче дифракции СНЧ электромагнитного поля на локальной неоднородности полости Земля-ионосфера в случае расположения неоднородности над одним из концов радиотрассы.
Иной подход развивается в задачах о распространении радиоволн над неоднородной в элекгрическом отношении земной поверхностью, когда неоднородность учитывается введением неоднородного поверхностного импеданса, задаваемого как функция координат на граничной поверхности волновода [137, 138, 139, 140, 141, 142]. Получаемые поверхностные интегральные уравнения решаются, в большинстве случаев, либо прямыми численными методами, либо решение определяется в Борцовском приближении. Сравнение Борцовского приближения решения двумерного интегрального уравнения с его численным решением для нерегулярности малой по сравнению с шириной первой зоны Френеля [140] показало
19
их эквивалентность. Анализ известных решений показывает, что наличие в получаемых соотношениях поверхностных интегралов, содержащих быстро осциллирующие функции, приводит к необходимости обращать (если решать задачу прямо, аппроксимируя интегральный операгор в большом числе узлов) получаемые при этом линейные алгебраические системы высокого порядка. Это накладывает слишком жесткие ограничения на размеры области нерегулярности, что в конечном итоге не дает заметных преимуществ прямому численному решению [131] перед решением в виде первого члена ряда теории возмущений [132, 134, 135, 136]. Авторы [141] проводят предварительное асимптотическое интегрирование по поперечной к направлению распространения координате, чтобы упростить дальнейшее решение, оставаясь в дальнейшем в рамках первого шага теории возмущений. Некоторые схемы приближенного интегрирования по поперечной координате рассмотрены в работах [89, 137]. Это позволяет исследовать значительные по протяженности (вдоль трассы распространения) неоднородности, однако, информация об их поперечной структуре в результате такой операции оказывается потерянной. Предложенная в [138] формула для учета влияния на функцию ослабления земной волны нерегулярностей трассы с учетом их поперечных размеров, получена без достаточного математического обоснования и носит полуэмпирический характер. Несколько особняком в предложенном списке стоит работа [139], авторы которой развивают методы параболического уравнения, используя алгоритм задачи Коши. В [139] сравниваются решения, полученные методом сеток в двумерном и трехмерном вариантах, для задачи распространения над неоднородностями земной поверхности в изотропном сферическом волноводе с имиедансными стенками при возбуждении поля вертикальным электрическим диполем. Как известно [3], замена системы уравнений эллиптического типа (уравнения Максвелла) на уравнения параболического типа вносит принципиально не уменьшаемую погрешность, оценить которую достаточно сложно. При этой замене отбрасываются также волны, распространяющиеся во встречном направлении. Численное решение [142] двумерного интегрального уравнения для трехмерной задачи распространения над плоской Землей с неоднородностями поверхностного импеданса сравнивается с экспериментальными
20
данными, полученными в [143]. Результаты находятся в хорошем соответствии везде, за исключением небольшой области расположенной непосредственно за неоднородностью, что авторы [142] объясняют приближенным характером импедансных граничных условий. Метод физического моделирования был использован в [145] для изучения влияния крупно масштабной ионосферной неоднородности, связанной с заметным понижением ионосферы. Отметим, что все перечисленные работы, за исключением [132, 134], используют скалярную постановку задачи, сущность которой состоит в пренебрежении деполяризацией поля при отражении от нерегулярных стенок волновода. Обоснование возможности такого подхода для задач о распространении земной волны содержится в [86]. Что касается задач с ионосферными локальными неоднородностями, то сопоставление результатов скалярной и векторной постановок задачи, проведенное в [144], показало правомочность скалярного подхода в большинстве практически значимых случаев.
Принципиально иной подход к проблеме распространения электромагнитных волн в приземном волноводном канале в сложных условиях горизонтально и вертикально неоднородной анизотропной ионосферы в случае, когда электромагнитное поле и свойства среды зависят от всех трех пространственных координат, развивается автором [146, 147]. Изложенная в отмеченных работах поверхностная теория двумерного телеграфного уравнения всецело учитывает специфику распространения радиоволн СНЧ и КНЧ диапазонов, когда волновод Земля-ионосфера - тонкий, и среди нормальных волн выделяется и остается одна, являющаяся аналогом кабельной волны. Остальные моды - местные и затухают на горизонтальном масштабе расстояния порядка высоты волновода.
Как было отмечено выше, большинство авторов, исследовавших проблему влияния локальных неоднородностей ионосферы на распространение радиоволн СНЧ и СДВ диапазонов, моделировали возмущение волноводного канала с помощью постулируемой зависимости, от координат в горизонтальной плоскости, постоянной распространения основной волноводной моды (введением малого ее возмущения). Учет в рамках той же теории [129, 148] эффекта многомодовости, принципиально изменить ситуацию не может, т.к. перевозбуждение мод по-прежнему не учитывается.
21
В данной диссертации свойства ионосферной стенки волновода описываются неоднородным поверхностным импедансом, определяющим свойства вертикально стратифицированной ионосферы, и задаваемым как функция координат граничной поверхности, в окрестности которой происходит распространение электромагнитной волны. Такой подход не ограничен частотами и расстояниями, допускающими только одномодовое представление поля в регулярном волноводе Земля-ионосфера, он позволяет последовательно рассматривать возмущения, находящиеся непосредственно над источником, когда любое даже самое низкочастотное поле описывается бесконечной суммой нормальных волн, выделить в которой главную (определяющую постоянную распространения) до некоторых расстояний от источника оказывается невозможным. Введением импеданса можно учесть влияния неоднородного ионосферного слоя, ответственного за отражение радиоволн, и тем самым отделить задачу об отражении волн от ионосферы от задачи формирования поля в полости волновода. Импедансная модель волновода Земля-ионосфера является, по-видимому, наиболее простой и часто использовавшейся с начала периода расцвета внимания к проблеме волноводного распространения СДВ [149, 150, 151, 152], в том числе и в задачах, учитывавших влияние геомагнитного поля [153, 154]. Описание отражательных характеристик плазменного ионосферного слоя с помощью импеданса позволяет относить его к эффективной отражающей поверхности, располагающейся на различных высотах внутри области существенной при отражении [155], и производить пересчет импеданса на новую высоту с сохранением характеристик отражающего слоя [156]. При этом, может быть найдена высота, на которой значительно ослабляется зависимость импеданса от угла падения волны на ионизованный слой [157], и, значит, импеданс, отнесенный к этой высоте, будет слабо зависеть от вида поля, падающего снизу на ионосферу. Определение значений импеданса и высоты, на которую его следует относить, можно производить по известным вертикальным профилям электронной концентрации и частоты соударений путем интегрирования системы нелинейных уравнений первого порядка [158, 159,
160]. Из результатов [161] следует, что усложнение импедансной модели, использующее первую производную эквивалентного импеданса по синусу угла
22
падения и определяющее эффективную высоту по минимуму этой производной, обладает наибольшей возможной точностью. Если использование импедансной постановки задачи в СНЧ диапазоне иногда требует дополнительного обоснования [162, 147], то для диапазона СДВ и ДВ (до частот примерно 100 кГц) такая постановка является общепризнанной. Отметим также, что все, только что процитированные результаты, получены на основе исследований либо задачи о падении плоской волны на горизонтально однородное, заполненное плазмой, полупространство, либо задачи о распространении волн в горизонтально однородном волноводе, т.е. в результате решения однородных задач.
В данной работе уже для неоднородной задачи используется та же самая идея о разделении исходной проблемы на две: на задачу построения математической модели (и придания физического смысла ее параметрам), описывающей условия распространения в околоземном пространстве с учетом горизонтальной (продольной и поперечной) неоднородности среды распространения, и задачу дифракции в промежутке между двумя импедансными поверхностями. Обратимся очень коротко к первой задаче - построения модели среды распространения. Не конкретизируя физические причины возникновения ионосферных неоднородностей, можно утверждать, что все они в конечном итоге выражаются соответствующими изменениями, например [163], вертикальных профилей электронной концентрации №е(х) и частоты соударений че(х) электронов (в рассматриваемом диапазоне частот влияние ионов не существенно). Эти профили как функции вертикальной координаты х над каждой точкой земной поверхности (г*, <р) можно использовать для вычисления импеданса ионосферы, отнесенного к некоторому, вполне конкретному [155, 156, 157,
161], уровню над поверхностью Земли и определяющего отражательные характеристики ионосферы. Иными словами, чтобы определить численные значения параметров импедансной модели, адекватно отражающей условия распространения радиоволн в присутствии локальных неоднородностей, в данной диссертации, для каждой точки (г, <р) на земной поверхности рассматривается задача распространения волн в плоском горизонтально-однородном волноводе, образованном Землей и вертикально-неоднородной ионосферой со свойствами (Ne(z) и че(г)), взятыми из
23
реальной ионосферы на вертикали, проходящей через рассматриваемую точку. После этого решается задача распространения радиоволн, возбуждаемых точечным источником, в вакуумном промежутке, ограниченном неоднородными импедансными поверхностями. Работоспособность такой модели была проверена в том числе и путем сопоставления расчетов, сделанных на ее основе в [164], с результатами эксперимента [165]. Условия же применимости модели станут ясны только после решения строгой задачи.
Сформулируем выводы, которые вытекают из обзора литературы.
- Проблема учета влияния локальных трехмерных неоднородностей нижней ионосферы или земной поверхности представляет несомненный интерес, как с теоретической точки зрения, поскольку относится к еще только начинающему развиваться направлению - исследованиям трехмерных волноводных задач, так и с точки зрения приложений к широкому кругу радиофизических, геофизических и экологических проблем.
- Имеющиеся к настоящему времени решения трехмерных волноводных задач охватывают очень узкую область возможных параметров таких неоднородностей, по большей части малых в масштабах либо длины волны, либо длины трассы распространения, либо размеров зоны Френеля. Большинство опубликованных решений - это скалярное приближение для задач дифракции и рассеяния векторных электромагнитных полей.
- Практически не нашли своего отражения в литературе, либо вовсе отсутствуют исследования крупномасштабных трехмерных неоднородностей, исследования векторных эффектов, например, деполяризации поля при рассеянии на трехмерной неоднородности, не учитывается реальная анизотропия приземного волноводного канала, обусловленная свойствами ионосферной плазмы и наличием геомагнитного поля.
24
Основываясь на только что приведенных положениях, сформулируем Цель данной диссертационной работы как развитие теории распространения низкочастотных электромагнитных волн в трехмерном локально нерегулярном околоземном волноводном канале на основе:
1.) разработки нового математического аппарата исследования волновых полей;
2.) построения решений новых задач, либо известных задач, но для областей изменения параметров не исследовавшихся ранее и доведения этих решений до численной реализации;
3.) исследования волновых полей для рассмотренных впервые трехмерных задач теории распространения радиоволн в приземном волноводе с различными модельными локальными неоднородностями.
Публикации материалов диссертации и её краткое содержание.
Все результаты, представленные в диссертации являются новыми и
оригинальными. Основное содержание диссертации опубликовано в статьях [75, 76, 127, 135, 144, 164, 166 - 184, 214] и отражено в докладах на научных семинарах, конференциях и симпозиумах [185 - 200, 215].
Вошедшие в диссертационную работу материалы представлялись на:
- XIV Всесоюзной конференции но распространению радиоволн, Ленинград, 1984, [185];
- XVI Всесоюзной конференции по распространению радиоволн, Харьков, 1990, [188];
- XVII Всероссийской конференции по распространению радиоволн,
Ульяновск, 1993, [191];
- XVIII Всероссийской конференции по распространению радиоволн,
С.Петербург, 1996, [192];
25
- Київ+380960830922