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

Богданов, Вадим Васильевич

Бесплатный автореферат и диссертация по наукам о земле на тему
Исследование в геосферных оболочках динамических процессов, инициированных солнечным и литосферным воздействиями
ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Исследование в геосферных оболочках динамических процессов, инициированных солнечным и литосферным воздействиями"

РОССИЙСКАЯ АКАДЕМИЯ НАУК ДАЛЬНЕВОСТОЧНОЕ ОТДЕЛЕНИЕ

ИНСТИТУТ КОСМОФИЗИЧЕСКИХ ИССЛЕДОВАНИЙ И РАСПРОСТРАНЕНИЯ РАДИОВОЛН

На правах рукописи УДК 537 22,550 343 551 510

БОГДАНОВ ВАДИМ ВАСИЛЬЕВИЧ

ИССЛЕДОВАНИЕ В ГЕОСФЕРНЫХ ОБОЛОЧКАХ ДИНАМИЧЕСКИХ ПРОЦЕССОВ, ИНИЦИИРОВАННЫХ СОЛНЕЧНЫМ И ЛИТОСФЕРНЫМ ВОЗДЕЙСТВИЯМИ

25 00 29 - физика атмосферы и гидросферы

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

I 8 СЕН 2008

с Паратунка, Камчатский край - 2008

003446128

Работа выполнена в Институте космофизических исследований I распространения радиоволн Дальневосточного отделения РАН

Официальные оппоненты:

доктор физико-математических наук, член-корр РАН Николаев Алексей Всеволодович

доктор физико-математических наук, профессор Новиков Лев Симонович

доктор физико-математических наук Кичигин Геннадий Николаевич

Ведущая организация: Институт космических исследований РАН

Защита состоится Щ ^^ 2008г в /V часов на заседании диссертационного совета Д 003 03401 при Институте солнечно-земной физики СО РАН по адресу 664033, Иркутск, ул Лермонтова, 126, а/я 291

С диссертацией можно ознакомиться в библиотеке ИСЗФ СО РАН

Автореферат разослан аЛгдсЬ2008г Ученый секретарь диссертационного совета к ф -м н

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Запуски КА с научной аппаратурой в околоземное и межпланетное пространство привели к серьезным открытиям, позволившим иначе взглянуть на геофизические явления и выделить роль Солнца в их формировании Активные процессы на Солнце сопровождаются электромагнитным излучением в широком диапазоне частот, генерацией солнечных космических лучей (CKJ1) и выбросом в межпланетное пространство огромных потоков плазмы с вмороженным в нее магнитным полем названных солнечным ветром Солнечный ветер достигает орбиты Земли и, взаимодействуя с ее магнитным полем, образует магнитосферу, полярные каспы, нейтральный токовый слой в хвосте и тд и является источником многих явлений, которые от внешней границы магнитосферы, определяемой положением ударной волны, через различные геосферные оболочки, трансформируясь и преобразуясь, распространяются к поверхности нашей планеты Эти явления задают направление развитию многих интересных и важных процессов, в той или иной мере влияющих на эко- и биосистемы Земли В связи с этим, в 1988 г была разработана международная программа STEP (Solar-Terrestrial Energy Program), призванная изучить всю цепочку преобразования и трансформации энергии от Солнца к Земле Центральное место в этой программе занимала тема «Общие механизмы в солнечно-земной системе» В результате проведенных исследований пришло понимание, что, несмотря на разницу в масштабах формиру ющихся событий, процессы в фотосфере и хромосфере Солнца во многом аналогичны процессам, которые развиваются в магнитосфере Земли и наблюдаются в экспериментальных лабораторных установках Была осознана важность и необходимость дальнейших исследований по изучению влияния Солнца на околоземное пространство и взаимодействие гео-сферных оболочек Это послужило причиной разработки программы «Солнечная активность и физические процессы в системе Солнце-Земля» (координатор академик Жеребцов Г А ) в рамках Программы фундаментальных исследований Президиума РАН

Многие явления, связанные с макроскопическими перемещениями плазмы, развитием различных плазменных иеустойчивостей, возникающих на магнитопаузе, в нейтральном слое и в районе каспа так или иначе проявляются в различных геосферных оболочках Следовательно, нейтральный слои и высокоширотные каспы в магнитосфере Земли, нейтральные области в районе солнечных пятен, нейтральные токовые слои, образующиеся при 9-разрядах, - те все нестационарные области магнитных полей, в которых реализуются нейтральные линии и нейтральные точки, являются источником многих интересных процессов в астрофизике, геофизике и, по-видимому, в установках лабораторного масштаба Эти яатения представляют собой важнейшую проблему в физике плазмы - проблему пересоединения магнитных силовых линий В этих случаях реализуется универсальный механизм, приводящий к развитию неравновесных процессов и к преобразованию магнитной энергии в энергию плазмы Этот механизм имеет и общефизическое значение, так как открывает путь к исследованию и развитию новых плазменных методов ускорения заряженных частиц Такие теоретические исследования важны и с точки зрения обоснования новых механизмов разогрева в термоядерных исследованиях Кроме того, экспериментально было установлено, что перестройка магнитного поля и эффект появления энергичных электронов (собственно сам вспышечный процесс) наблюдались только при

выполнении некоторых условий, одним из которых является высокая степень вмороженности поля в плазму При макроскопических движениях плазмы эта вмороженность должна проявляться на магнитных поверхностях, охватывающих нулевые точки и линии В аксиально-симметричных полях указанное условие отождествляется с сохранением третьего адиабатического инварианта, который был введен Нортропом и Теллером как магнитный поток, охватываемый дрейфовой поверхностью В таком виде он не совпадает с канонической формой адиабатических инвариантов движения

Наряду с фундаментальными открытиями, осуществленными с помощью КА, появилась серьезная проблема, связанная с надежностью и ресурсом их работы при эксплуатации в натурных условиях Как показал опыт, двигаясь по орбите, КА пересекают различные геосферные оболочки и непрерывно подвергаются воздействию волновой и корпускулярной радиации в виде высокоэнергичных заряженных частиц, солнечных и галактических космических лучей Такое взаимодействие не проходит для аппаратов бесследно и приводит к дифференциальному и объемному заряжению (электризации) и, как следствие, к эрозии диэлектрических и защитных элементов конструкции КА Эксплуатация КА показала, что электризация ответственна за целый ряд сбоев и катастрофических отказов в работе не только научной, но и штатной, отвечающей за нормальную работу всего спутника, аппаратуры Статистическая обработка данных выявила, что ~ 85% сбоев и отказов происходили в период магнитосферных возмущений Однако оставшаяся часть зарегистрированного аномального поведения аппаратуры (~15%) происходила в условиях спокойной магнитосферы Более того, были зарегистрированы не единичные случаи сбоев, возникавших в спокойной магнитосфере накануне сейсмических событий, когда подспутниковая точка ИСЗ проходила через район готовящегося землетрясения Поэтому эксплуатация КА поставила не только перед конструкторами задачу по разработке методов и рекомендаций по их защите от нежелательного действия активной внешней среды, , но и перед учеными была сформулирована серьезная проблема, связанная с необходимостью объяснить причины и механизмы проникновения на ионосферные (сотни км) и в магнитосферные высоты (тысячи км) возмущений, которые формируются в сейсмоопасных областях литосферы

Структура ионосферы, распределение ее параметров по высоте и по составу определяется плотностью атмосферы и ее химическим составом, спектральными характеристиками солнечного излучения Одновременно с регулярными изменениями параметров ионосферы, соответствующих суточному и сезонному ходу, а также определяемых активностью Солнца в сейсмоактивных регионах Земли может наблюдаться также и аномальное поведение распределения таких характеристик ионосферы, как Цч и соответствующих им действующих высот Ь , в эффектах проявления турбулизации слоев ¥2 и Еь (1;-зргеа<1 и Е-^ргеас!) и т д , наблюдаемых накануне землетрясения и задаваемых литосферно-ионосферными взаимодействиями Современными исследованиями установлено, что литосфера активно влияет на физические процессы, протекающие в верхних геосферных оболочках Установлено, что период, предшествующий землетрясению (ЗТ), сопровождается развитием целого ряда явлений, имеющих механическую, гидрохимическую, электромагнитную природу и регистрируемых как на поверхности Земли, так и на ионосферных и магнитосферных высотах Каждое из этих явлений сопровождается формированием ряда определенных признаков, которые часто маскируются процессами, связанными с проявлением активности Солнца

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

Актуальность проблемы Изучение процессов, протекающих в различных геосферах на фоне солнечного и литосферного воздействия, обусловлено практической значимостью формирующихся в них явлений и непосредственно влияющих на физические процессы в околоземном пространстве и на Земле Проблемы, затронутые в данной диссертации, соответствуют программам фундаментальных исследований Президиума РАН на 2005г и находятся на пересечении трех направлений научная программа № 18 «Теплофизика и механика интенсивных энергетических воздействий» в части «Солнечный ветер генерация и взаимодействие с Землей и другими планетами», научная программа № 13 «Изменение окружающей среды и климата природные катастрофы» и программа №30 «Солнечная активность и физические про-

цессы в системе Солнце-Земля» Кроме того, актуальность выбранных исследований подтверждается тем, что основные направления диссертации соответствуют тематике исследований, активно проводящихся в мировых научных центрах, и входят в проблематику крупных современных международных программ (International Living with a Star Program, Climate and Weather of the Sun-Earth System, Lower Thermosphere Coupling Study, Magnetosphere-Ionosphere Coupling Storm/Substorm Effects Mid&Low Latitude), исследования по которым координируются крупнейшими научными организациями (SCOSTEP, COSPAR, URS1, AGU, EGS и др )

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

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

2) построение инженерной модели поверхностной и объемной электризации диэлектриков, используемых в качестве защитных элементов для КА и взаимодействующих в магнитосфере Земли с радиациси волновой и корпускулярной природы,

3) исследование методами вейвлст-преобразования, совмещенными с моделями авторег-рессии-проинтегрированного скользящего среднего (АРПСС), изменений в параметрах ионосферы на примере критической частоты f„F2 с целью обнаружения возможных аномальных эффектов в ссйсмоопасных регионах, которые могли бы быть отождествлены с литосферными процессами,

4) исследование статистическими методами вариаций частоты появления турбулнзацни слоя Е, и увеличения вероятностей турбулизацин слоя F2 (F-spread и E-spread) на фоне сейсмической активности,

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

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

- каноническое представление адиабатического инварианта движения заряженной частицы, дрейфовое приближение, кинетические уравнения Власова (Главы 2,3),

- упрощенная с учетом специфики магнитосферной плазмы кинетическая модель накопления диэлектриком объемного заряда Роуза-Фаулера-Вайсберга для расчета электризации конструктивных элементов КА в условиях магнитосферы (объемная модель электризации) (Глава 3),

- методы дискретного и непрерывного вейвлет-преобразования и класс моделей авторегрессии

- проинтегрированного скользящего среднего - для исследования на примере критической частоты f„F2 динамики ионосферного процесса (Глава 4),

- статистические методы с целью анализа вариаций ионосферных параметров и вероятностей турбулизацин слоев F2 и Es (F-spread и Е- spread) на фоне сейсмической активности (Глава 4)

- элементы теоретико-вероятностного подхода для построения статистической модели на основе обработки конкретного каталога землетрясении и методы вейвлет-анализа, используемые для изучения динамики функции распределения сейсмических событий по глубине и по времени, с целью выявления возможных прогностических признаков готовящегося сейсмического события в конкретной области п-ова на полуострове Камчатка (Глава 5)

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

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

геомагнитной ловушке во время глобальных магнитосферных возмущениях Показано, что нестационарное геомагнитное поле в аксиальном приближении сохраняет близкий к невозмущенному (опорному) полю характер широтной зависимости вдоль магнитной силовой линии в уравнениях соответствующих постоянству двух адиабатических инвариантов Ф=сопз1 и ^сош1 На основе полученных результатов построена динамическую систему координат, описывающая перенос плазмы поперек магнитного поля в условиях нестационарной магнитосферы Данная система координат является дальнейшим развитием принципов построения естественной системы координат Мак-Илвейна,

2 На основе упрощенной кинетической модели Роуза-Фаулера-Вайсберга разработана самосогласованная модель объемного заряжения диэлектриков, которая позволяет рассчитывать в них распределения электрического поля и соответствующего потенциала, наводимых реальными потоками магнитосферной плазмы, и выработать на стадии конструирования соответствующие защитные меры для минимизации пагубного воздействия электризации на КА в натурных условиях Показано, что наиболее опасными зонами для высокоэллиптических КА являются участки орбиты, пересекающие силовые линии, основания которых проектируются в авроральный овал Использование элементов построенной динамической системы координат позволило на конкретном примере геомагнитной бури (по данным ИСЗ «ОГО-3») проследить на восстановительной фазе возмущения динамику расширения магнитосферы и рассчитать на основе экспериментально полученных потоков заряженных частиц объемное заряжение кварцевых защитных покрытии солнечных батарей КА, находящегося на геостационарной орбите Предложен возможный механизм появления в плазмосфере холодных и горячих протонов соответственно на главной и восстановительных фазах магнитной бури На основе модели объемной электризации диэлектриков предложен возможный механизм формирования сбоев и отказов в работе систем КА, пролетающего над сейсмоактивным районом накануне землетрясения

3 На основе статистического анализа параметров ионосферы спорадических слоев Еб и ¥2, подверженных диффузионному расширению, так называемые Ь^ргеаб, и Р-5ргеас1, обнаружено, что за 1-2 дня до землетрясения (М>4) частота появления [й-зргсас! увеличивается в предполуночные часы, т е турбулизация плазмы спорадических Еэ-слоев усиливается В то же время в сейсмоионосферных явлениях характерно уменьшение значения турбулизации слоя ¥2 (Р-эргеа!!) перед землетрясениями (М>6 0) Статистическая достоверность составляет 0 95 для всех обнаруженных эффектов

4 Впервые для анализа динамики ионосферных параметров на примере критической частоты №2 применены современные математические методы вейвлет-преобразования совмещенные с классом моделей авторегрессии - проинтегрированного скользящего среднего с целью обнаружения влияния на их поведение процессов литосферного происхождения

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

Виедренне и востребованность результатов работы. Представленные в диссертации результаты по динамике магнитосферной плазмы в дрейфовом приближении с учетом сохранения трех адиабатических инвариантов движения используются при прогнозе сбоев и отказов на высокоэллиптических КА производства НПО им С А Лавочкина Модель объемной электризации, руководство для конструкторов по минимизации влияния заряжения изделий статическим электричеством в условиях натурной эксплуатации, практические рекомендации и тд применяются в НПО на стадии проектирования (Акт внедрения от 26 10 2005 г ) На разработку устройства для зарядов электризации диэлектриков диссертантом получено авторское свидетельство № 1589789 от 1990г Методика прогноза сейсмических событий на п-ове Камчатка на основе ионосферных и литосферных прогностических признаков используется при составлении еженедельного прогноза Камчатским отделением Федерального центра прогнозирования землетрясений (г Петропавловск-Камчатский) Работы по диссертации выполнялись в рамках про-

ектов РФФИ «Разработка концепции применения космических средств с целью экологического мониторинга, предупреждения стихийных бедствий и оценки ущерба от их воздействия на Камчатке» № 97-05-96576-п (исполнитель), «Исследование высокочастотных мехаиоэлектричс-ских явлений возникающих в областях напряженного состояния земной коры» № 03-05-65302-а (исполнитель), а также при поддержке интеграционных проектов СО РАН и ДВО РАН «Экспериментальные исследования влияния космофизических и геофизических возмущении на формирование аэрозоля в средней и верхней атмосферах» (исполнитель), «Исследование лнто-сферно-атмосферных связей и механизмов взаимодействия геосфериых оболочек при различных уровнях солнечной и геомагнитной активности» (руководитель)

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

Апробация работы Материалы диссертации докладывались на объединенных семинарах Института прикладных проблем механики и математики АН УССР в 1978г и 1980г , конференциях «Электризация КА» (Москва 1987 1988, 1990) и рабочих совещаниях в НИИЯФ (МГУ), российско-японском семинаре «Проблемы геодинамики и прогноза землетрясении», ИТиГ, Хабаровск, 2000, па международных конференциях «Сотнечно-земные связи и электромагнитные предвестники землетрясений» в 1998 2001, 2004, 2007 (с Паратунка, Камчатская область). International Workshop on Seismo-Electromagnetics of N AS DA (IWSE 2000), Tokyo, Japan, 2000, «Солнечно-земная физика», Иркутск, 2004, XI Joint International Symposium "Atmospheric and Ocean Optics Atmospheric Physics" Томск, 2004, European Geosciences Union (EGU) General Assembly, Vienna, Austria, 2005, International Conference on Soft Computing and Measurements Санкт-Петербург 2005, 2006 2007, IAMAS Scientific Assembly, Beipng, China, 2005 European Gcosciences Union (EGU), General Assembly, Vienna, Austria, 2006, International Union for Geophysics and Geodesy (IUGG) XXIV General Assembly, Perugia, July 2-13, 2007, Italy, Физика плазмы в межпланетном пространстве, ИКИ, 5-8 февраля, 2008г

IIa защиту выносятся

1 Доказательство эквивалентности канонической и потоковой форм представления третьего адиабатического инварианта движения заряженной частицы в аксиальных магнитных полях в дрейфовом приближеиии и в приближении среднего дрейфа, полученные допотнитсль-ные инвариантные формы, необходимые для построения системы координат (кинетичеснои модели переноса), справедливой в переменных магнитных полях (см п 3)

2 Обоснованный механизм нарушения третьего адиабатического ипвариа]па в дрейфовом приближении и приближении среднего дрейфа для аксиальных и слабо отличных от аксиальных магнитных полей обусловленный ускорением радиального дрейфа R Произведена оценка точности сохранения третьего инварианта, связанная с возможностью пренебречь кумулятивными эффектами этого ускорения в переменном поле

3 Проведенный анализ характера деформации инвариантиых оболочек в геомагнитной ловушке во время глобальных магнитосферных возмущений, и полученный на основе этого анализа с привлечением дополнительных инвариантных форм (см п 1) возможный механизм появления в плазмосфере холодных и горячих протонов соответствегаго на главной и восстановительных фазах магнитной бури

4 Обобщенная система координат, которая при условии сохранения трех адиабатических инвариантов переводит ф>нкцию распределения частиц в реальном деформированном поле в соответствующую функцию невозмущенного (опорного) поля путем адиабатической вариации магнитного поля, что позволяет судить не только о степени вмороженности потоков в магнитное поле, но и о проникновении новых частиц в ловушку

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

6 Впервые выполненные исследования динамики изменений критической частоты f„F2 методами дискретного и непрерывного вейвлет-преобразования совмещенными с моделями авторегрессии - проинтегрированного скользящего среднего, дающие возможность определить аномальное поведение ионосферных параметров, обусловленное литосферными процессами перед сильными сейсмическими событиям

7 Впервые проведенный статистический анализ параметров спорадических слоев Es и F2 ионосферы, подверженных диффузионному расширению, так называемые Es-spread и F-spread

8 Впервые предложенная статистическая модель сейсмичности п-ова Камчатка, выделенный на ее основе методами вейвлет-анализа прогностический признак Кроноцкого землетрясения, который сформировался в последний год перед событием, инвариантная форма закона повторяемости сейсмических событий в вероятностной интерпретации, ее связь с законом Гутенберга-Рихтера

Публикации Результаты выполненной работы опубликованы в российских, зарубежных изданиях и тематических сборниках Основная часть результатов диссертации представлена в работах [1-51] и включает две монографии, публикации в журналах «Космические исследования», «Метрология», «Доклады АН», «Natural Hazards and Earth System Sciences», «Известия ВУЗов», в сборниках докладов конференций, и выполненных в виде отчетов с НПО им С А. Лавочкина. Работы автора, выполненные по теме диссертации, но не вошедшие в нее в полном объеме и не включенные в основные выводы, опубликованы в [52-69] Автор имеет авторское свидетельство на изобретение и Акт внедрения своих разработок в НПО им С А Лавочкина [7,13]

Личный вклад автора Работы, представленные в диссертации, выполнены автором как самостоятельно, так и в сотрудничестве с коллегами из Института космофизических исследований и распространения радиоволн ДВО РАН, НПО им С А Лавочкина по программе «Электризация», ЦНИИмаш, Новосибирского государственного университета, Института космофизических исследований и аэрономии СО РАН (г Якутск), Института физики Земли РАН и др Теоретические исследования по вспышечным процессам автор выполнял согласно пункту 18 2 1 координационного плана Академии наук СССР по тематике работ лаборатории Молекулярной акустики Московской государственной академии приборостроения и связи (бывший Всесоюзный заочный машиностроительный институт) По этому пункту теоретически и экспериментально (на примере квадруполыюго магнитного поля с нулевой точкой) исследовался процесс преобразования магнитной энергии в энергию мощных магнитогидродинамических движений плазмы, которые сопровождают процесс пересоединения магнитных силовых линий Автором лично выполнено теоретическое исследование правомерности замены канонической формы третьего адиабатического инварианта на потоковую форму в полях слабо отличных от аксиально-симметричных магнитных полей, построена модель объемного заряжения диэлектриков, на основе теоретико-множественного подхода и соответствующей обработки каталога землетрясений построена статистическая модель сейсмичности п-ова Камчатка и получена вероятностная интерпретация закона повторяемости Гутенберга-Рихтера, справедливой в любом сейсмоактивном регионе Земли, для анализа аномального поведения геофизических параметров в сейсмоактивных регионах применены методы вейвлет-разложения. Работы велись в рамках программы фундаментальных исследований Президиума РАН на 2005г № 30 Применение научных разработок автора по проблеме электризация КА в НПО им С А. Лавочкина оформлено соответствующим Актом внедрения

Объем работы и структура диссертации- Диссертация состоит из введения, пяти глав, заключения, приложений, списка литературы, включающего 289 наименований Общий объем диссертации 265 страниц машинописного текста Изних-237стр текста, список литературы -20 стр, приложения - 8 стр Диссертация содержит 54 рисунка, 7 таблиц

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

Первая глава, являясь обзорной, посвящена описанию состояния вопросов в контексте с решаемыми в диссертации задачами Кратко представлена морфология гслиогеофизических оболочек начиная с фотосферы на Солнце и заканчивая литосферой Земли Описываются активные процессы на Солнце (хромосферпые попытки эруптивные протуберанцы, корональпые выбросы масс) оказывающие через корпускулярное и волновое излучение, непосредственное влияние на протекание физических процессов в межпланетном и околоземном космическом пространстве в различных слоя\ атмосферы и литосфере Кроме того, активные процессы на Солнце тесно связаны с возникновением -85% различных нештатных ситуации на борту КА. влияя на надежность работы штатной и научной аппаратуры и ресурс спутника Рассматривается структура магнитосферы, которая формируется за счет взаимотеиствня солнечного ветра и межпланетного магнитного поля с магнитным полем Земли При этом особо отмечается общность в механизмах формирования активных процессов на Сошце в районах полярных каспов и нейтральною токового слоя в хвосте, те в таких магнитных конфигурациях, которые обта-дают нейтральными линиями и точками Такие магнитные конфигурации реализмотся н в наземных экспериментальных установках при изучении динамики вспышечных процессов в остроугольных магнитных лов\шках Рассматриваются приближенные методы решения уравнения движения заряженной частицы в скрещенных электрических и магнитных полях, приводящие к понятию трех адиабатических инвариантах движения [I, .1 Ф Отмечается что если первые два

инварианта ц = и Л = | 111V ц (11 совпадают с канонической формой адиабати-

ческих инвариантов I то третий Ф = ^ Н <1 Б сразу был введен как поток магнитного

поля охватываемого оболочкой по которой дрейфуют ведущие центры частиц, и его сохранение скорее определялось наличием аксиальной симметрии в области захвата магнитосферы Земли (дипотьное приближение), чем каноническим условием сохранения Доказательство его постоянства проводилось в пренебрежении самим дрейфом, без которого этот инвариант теряет смысл (Глава 2)

Наличие трех инвариантов объясняет механизм захвата и удержания частиц магнитным полем На основе постоянства первых двух адиабатических инвариантов |[ и 1 для описания распределения корпускулярной радиации в области захвата магнитосферы Земли в дипольном приближении используется так называемая система координат Мак-Илвейна Данная система координат справедлива только в постоянных магнитных полях Учет сохранения третьего адиабатического инварианта Ф дает возможность построил, динамическую систему координат, справедливую в переменных полях, удовлетворяющих условию постоянства Ф (Глава 3)

При своем движении по околоземной орбите КА проходит через области магнитосферы, которые в иевозмущенных условиях характеризуются различными плотностями и энергиями заряженных частиц Взаимодействие внешних защитных диэлектрических поверхностей с активной внешней средой приводит к объемному и дифференциальному заряжению до высоких отрицательных потенциалов Изменение внешних условий при пересечении пограничных слоев приводит к перераспределешпо потенциалов и как следствие, к формированию электрических разрядов, что ответственно за отдельные сбои или отказы в работе штатной и научной аппаратуры а также к деградации внешних поверхностей КА В случае повышения активности Солнца. связанной с формированием мощных хромосферных вспышек, образованием эруптивных протуберанцев и корональпыми выбросами солнечной массы, условия эксплуатации значительно ухудшаются и частота сбоев и отказов, а также степень их влияния на нормальное функционирование КА значительно усиливается Для правильного выбора защитных материалов, выбора оптимальных орбит с точки зрения решаемых задач и уменьшения влияния активных внешних факторов космического пространства, необходимо построите достаточно простой математической модели объемной электризации диэлектрических покрытии (Глава 3)

Физические процессы, формирующиеся в геосферной цепочке «литосфера-нейтральная атмосфера-ионосфера-магнитосфера», имеют планетарный масштаб, и если в проблеме взаимодействия солнечного ветра с магнитным полем Земли достигнуто определенное понимание, то

в изучении процесса взаимодействия между литосферой и располагаемыми выше геосферными оболочками делаются только первые шаги (Полагают, что инжекция радона и аэрозолей в приземную атмосферу в сейсмоактивных районах приводит к перераспределению атмосферного электрического поля Эти изменения передаются на ионосферные высоты, что приводит к аномальному поведению параметров ионосферы Кроме того, инжекция радона в приземный слой может быть ответственной за свечение атмосферы в инфракрасной части спектра (см, например, [61, 69])) В то же время из общих соображений ясно, что если открытая физическая система находится вдали от равновесия, то она чувствительна к любым изменениям внешних полей Следовательно, резкие колебания в космической погоде, обусловленные повышением активности Солнца, могут привести к срабатыванию в геосферных оболочках неких природных механизмов, которые могут стать причиной возникновения катастрофических явлений таких, как тайфуны, наводнения, оползни, цунами, землетрясения и тд Если такого механизма не существует, то задача выполнения краткосрочного и тем более оперативного прогноза в рамках механического подхода становится проблематичной В этом случае необходимо точно знать размеры сейсмогенерирующей области, неоднородность взаимно смещающихся поверхностей разломов, прочность пород, развивающиеся в них напряжения, а, следовательно, и силу, ответственную за рост напряжений, те множество характеристик, экспериментально определить которые известными методами невозможно В такой ситуации перспективным представляется следующий путь поиск аномального поведения геофизических параметров в явлениях, которые могут регистрироваться на земной поверхности, ионосферно-магнитосферных высотах и их систематизация применительно к данному сейсмоопасному региону К таким параметрам относятся критические частоты £^2, и fc.Es, действующие высоты Ь , Р-зргеасЗ, Ез-5ргеа(1 и тд (Глава4)

Для изучения внешних причин, инициирующих подготовку сейсмического явления и ответственных за активизацию подготовительных процессов в литосфере, отклики которых проявляются в магнитосфере и на ионосферных высотах, необходимо сейсмический режим того или иного региона представить в удобной для исследований математической форме Такой формой является интегральная функция распределения или плотность распределения и построенная на их основе статистическая модель сейсмичности п-ова Камчатка (Глава 5)

Вторая глава посвящена теоретическому исследованию третьего адиабатического инварианта движения заряженной частицы в произвольном аксиально-симметричном и слабо отличном от аксиального магнитных полях Этот инвариант вводится как магнитный поток Ф, охватывающий поверхность, по которой дрейфуют частицы (Нортроп, Теллер, 1960) Сохранение этого инварианта обычно отождествляется с условием вмороженности, согласно которому плазма перемещается в пространстве на магнитных поверхностях постоянного магнитного потока Введенный в потоковой форме этот инвариант не совпадает с канонической формой адиабатических инвариантов движения I Более того, доказательство сохранения Ф на основе постоянства обобщенного импульса в аксиальном поле проводилось в пренебрежении самим дрейфом частицы, что не корректно, поскольку в этом случае теряется смысл этого инварианта как в потоковой, так и канонической формах (например, Хесс, Редерер, 1972) Следует отметить, что на основе введенного таким образом третьего инварианта и условия его нарушения строилась теория формирования радиационных поясов Земли Однако динамика заряженных частиц в радиационных поясах и во внешних областях каспа хорошо описывается обратимыми кинетическими уравнениями Власова, и процессы переноса в них определяются точностью сохранения всех трех инвариантов движения Так, например, на основе экспериментов, выполненных на установке ТС-3 (ФИАН), было доказано, что перенос плазмы поперек магнитного поля с последующим формированием токового слоя происходит с высокой степенью вмороженности инвариантных магнитных оболочек в плазму Следовательно, помимо полу феноменологического введения третьего адиабатического инварианта в потоковой форме необходимо выполнить строгий математический анализ и с точки зрения замены канонической формы I на потоковую Ф, и с точки зрения точности сохранения инварианта Исследование проводилось на основе лагранжевого формализма раздельно для частиц, двигающихся в плоскости экватора аксиально-симметричного поля ( Уи =0), и для частиц, имеющих составляющие скорости вдоль силовой линии (Для полей, обладающих слабой неаксильностью, рассматривались движения частиц только с Уц=0) Тк для аксиально-симметрических полей в цилиндрической системе координат азимутальный угол ф является циклической координатой, то соответствующий обобщенный импульсРф =соп51 В дрейфовом приближении (р,«К) обобщенный импульс имеет вид

P =mRuT + eRA» (1)

9 С

В (1) и,,- азимутальная дрейфовая скорость Дрейфовому приближению соответствует случай mRu<p «(е c)RAcp Но пренебрежение членом mRu9 в (1) недопустимо, поскольку теряется смысл третьего инварианта в формах I и Ф Проведенный анализ третьего адиабатического инварианта показал, что 1 = const и Ф = const в стационарных полях и в нестационарных при отсутствии радиального ускорения В то же время I и Ф в переменных полях являются абсолютными инвариантами, если R = const, что не следует из дрейфовых уравнений uR = = R =с Е/Н ^ const в общем случае Однако условие адиабатической инвариантности

ч>

«1, (2)

ограничивающего скорость изменения поля за период дрейфа Т.ф, дает неравенство Цф » Ur, что позволяет в рамках дрейфового приближения пренебречь ускорением инвариантной поверхности Ur = R , как величиной более высокого порядка малости по сравнению с «дрожанием» этого приближения Поэтому точность сохранения Ф удовлетворяет адиабатическому приближению для небольших интервалов времени, когда можно пренебречь кумулятивными эффектами, связанными с малым, конечным ускорением CIr На основе представленных в явном виде инвариантов Ф и I для дрейфового движения в плоскости экватора в диссертации получены дополнительные инвариантные формы Rv = const и R:H = const

Рассмотрение третьего инварианта в форме Ф и I в дрейфовом приближении с Vu Ф О показывает, что в этом случае необходимо перейти к приближению среднего дрейфа В результате получаем

РФ =m RUi > + - iRA,/ = const (3)

'с 1

(Ru^^^+^-J-J^onst (4)

1 т"

с \<а,/шти

где приняты следующие обозначения ( } = —— J dt, J=m |v[,dt=const - второй адиабати-

ческий инвариант, к - некоторая геометрическая характеристика поля, шл - ларморовская частота Возможность считать (3) и (4) постоянными величинами связана с усредняемостью к и к/со , вдоль силовой линии Поскольку положение «зеркальных» точек в постоянном магнитном поле не меняется, то I~ ,'Ru?, = const и Ф~i RA/ -const В случае нестационарных полей, но при выполнении условия (2), что соответствует неравенству Иф » UR, то, как доказано в Приложении А, к и k/Mt являются усредняемыми величинами Поэтому и в случае переменных полей I- Яцф = const и Ф~;ЯАФ/ = const, как это следует из (3) и (4) с точностью

J=const для небольших интервалов времени Анализ экспериментальных данных о переносе плазмы во внешних областях остроугольного магнитного поля, удовлетворяющего дрейфовому приближению, должен показать возможность описания процессов переноса на основе предположения об адиабатической инвариантности потока Ф с учетом поправок, связанных с обязательным его нарушением вследствие кумулятивных эффектов ускорения радиального дрейфа, т е соответствует вмороженности плазмы в дрейфовые оболочки Наличие вмороженности экспериментально подтверждено на установках ТС-3 (ФИАН) и, как доказано, является определяющим условий формирования вспышечного процесса Именно эти поправки, дающие возможность оценить нарушение вмороженности плазмы и, как следствие, возникающую дополнительную диссипацию при распространении магнитогидродинамических волн можно найти на основе развитой в этой главе теории

Для анализа движения в дрейфовом приближении в несимметричных полях была выбрана подвижная система координат К', которая вращается относительно цилиндрической системы К с угловой скоростью Эф/St, локально связанной с угловой скоростью ведущего центра частицы Центр вращения К' совпадает с центром цилиндрической системы координат К Переход к описанию движения заряженной частицы в подвижной системе координат К' относительно кинематических соотношений эквивалентен переходу к аксиально-симметричному переменному полю, приводящему к возникновению электрического поля и, как следствие, к радиальному дрейфу Тем самым задача свелась к уже решенной задаче, и возможность пренебречь ускорением R приводит к требованию наличия слабой аксиальной несимметрии (1/Н)(ЭН/Эф)«1

Третья глава посвящена построению кинетической модели переноса плазмы, справедливой в возмущенном магнитном поле, и модели объемной электризации диэлектрических материалов под действием факторов космического пространства (ФКП) в виде корпускулярного и волнового (ультрафиолетового и рентгеновского) излучения

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

Для интерпретации результатов измерения заряженных частиц во внутренней области магнитосферы используется система координат Мак-Илвейна С учетом сохранения первого адиабатического инварианта Ц и второго J и в предположении о подобии распределения частиц разной энергии на одной инвариантной поверхности дает возможность свести функцию распределения с шестью фазовыми переменными к функции от двух и осуществить следующий переход f(r,q>,X,vr, v<p,va.) ->f(r,X) При этом оказалось наиболее удобным перейти к системе двух переменных Н и L, Н - магнитное поле, I^Ro/R,,, параметр, характеризующий дрейфовую оболочку (R0 - расстояние от центра магнитного диполя до вершины силовой линии в плоскости экватора, R31( - средний радиус Земли) Однако в нестационарных полях система координат (H,L) теряет свой смысл В то же время на основе сохранения всех трех адиабатических инвариантов можно построить кинетическую модель переноса (динамическую систему координат), которая справедлива в переменных полях Это позволяет представить функцию

распределения в виде f=f(L,V, Ct°), где V и а°- соответственно скорость и питч-угол в плоскости экватора В основе построения кинетической модели лежит представление об адиабатической деформации магнитного поля, и модель приобретает смысл, когда возможно однозначное отождествление возмущенной инвариантной поверхности с опорной (невозмущенной), например, дипольной Как показывает анализ, для однозначной идентификации поведения магнитных оболочек и дрейфующих по поверхности этих оболочек заряженных частиц необходимо применение теоремы, сущность которой сводится к следующему если напряженность поля Н на поверхности магнитной оболочки, характеризуемой вмороженной в эту оболочку плазмой, уменьшается, то в области уменьшения Н радиус кривизны магнитной силовой линии растет, и, наоборот, в случае возрастания Н радиус кривизны уменьшается (К Лонгмайр, 1966) Аналогом этой теоремы для возмущенного аксиального поля в приближении Ф~ const является, полученная в Главе 2, инвариантность величин Rv и R2H На основе сформулированной теоремы выбираются уравнения, описывающие возмущенные силовые линии опроного (например, ди-польного) магнитного поля r-R0cosk X, (R„ -экваториальное расстояние до вершины силовой линии, Х- широта, k=2+l 5[1-( t/т)2], т- характерное время восстановления поля после максимальной деформации t = 0, а к = 3 5) Выполненные на примере дипольпого поля оценки, основанные на дополнительных инвариантных формах, показывают, что в возмущенных уравнениях J=const и <D=const можно сохранить невозмущенный (дипольный) вид широтной зави-

симости На основе этих оценок и с учетом сохранения всех трех инвариантов построена кинетическая модель переноса плазмы поперек диполыюго поля Физический смысл этих моделей заключается в том, что они дают возможность на основе проведенных измерений направленного J и всенаправленного потоков I плазмы получить соответствующие невозмущенные потоки путем адиабатической вариации магнитного поля от возмущенного к невозмущенному полю и осуществить переходы вида

(5)

А[гД,Н,а,у) — ^

Я г, А.,Н,а, V)/* (5)

Пересчет реальных потоков в невозмущенные согласно (5) или (5') позволяет сделать вывод не только о вмороженности потоков плазмы в поле, но и о возможном проникновении частиц как с других инвариантных оболочек, так и инжекции частиц от внешних источников

Вторая часть третьей главы посвящена построению моделей расчета влияния факторов космического пространства (ФКП) в виде корпускулярного и волнового (ультрафиолетового и рентгеновского) излучения на надежность и ресурс эксплуатации высокоэллиптических КА

Комплексное действие ФКП приводит как к поверхностной дифференциальной электризации защитных покрытий, так и к объемному заряжению диэлектрических элементов конструкций КА Поскольку стоимость эксплуатации космических аппаратов значительно меньше стоимости их разработки и вывода на заданную орбиту, то возникает естественная задача о значительном увеличении срока их активного существования При этом важную роль в обеспечении длительной и, одновременно с этим, безотказной работы КА играет конструкционная стойкость используемых материалов и элементов, а также разработка надежных методов защиты штатной и научной аппаратуры от действия отмеченных факторов Одним из направлений в этих работах является построение моделей поверхностной и объемной электризации, которые позволяют оценить возможные уровни заряжения и рассчитать распределения электрических потенциалов и полей в диэлектриках при их взаимодействии с потоками космической плазмы Эта информация необходима конструкторам и разработчикам космической техники для того, чтобы уже на стадии проектирования в случаях, когда расчетные уровни заряжения превышают некоторые пороговые уровни заряжения (обычно > 100 В), предусмотреть применение различных способов и методов защиты штатной аппаратуры изделия от пагубного проявления различных эффектов, сопутствующих процессу электризации

Для построения строгой теории взаимодействия потока заряженных частиц с диэлектриком необходимо воспользоваться теорией накопления объемного заряда Роуза-Фаулера-Вайсберга (РФВ) На основе системы уравнений модели РФВ строились различные теории, которые описывают такие явления, возникающие при облучении высокоомных диэлектриков потоком заряженных частиц, как электропроводность, диэлектрические свойства, ползучесть, долговечность и тд Однако эта система уравнений сложна и мало пригодна для проведения оценочных расчетов объемной электризации диэлектриков Более того, предположения, лежащие в основе построения модели РФВ и проводимых на их основе расчетов, не соответствуют условиям, с которыми сталкиваются аппараты, взаимодействующие с реальной магнитосфер-ной плазмой (Например, предположения о монохроматичности потока заряженных частиц и о длине пробега электронов А., для которой полагается, что она много больше толщины облучаемого материала Ь и т д ) Однако реальная магнитосферная плазма достаточно хорошо описывается максвелловским распределением (более адекватное описание дается двухмаксвеллов-ским распределением частиц по скоростям) При этом концентрация электронов высокоэнергичной магнитосферной плазмы колеблется от 0 1 до 10 частиц на см3 при температуре эВ Наличие малой концентрации потоков первичных электронов позволяет в рамках оценочных моделей пренебречь уровнем динамического равновесия, обусловленного скоростью образования в диэлектриках электрон-дырочных пар и их рекомбинацией Следовательно, величиной установившегося равновесия, а вместе с этим и самим процессом генерации электрон-дырочных пар в нашем случае можно пренебречь Поэтому наведенная проводимость диэлектрика будет обусловлена термализованными электронами, что в приближении максвелловской плазмы значительно упрощает расчеты распределений профилей электрического поля и соответствующего потенциала При этом легко учитываются изменения в максвелловской функции распределения, которые возникают за счет кулоновского взаимодействия потока первичных магнитосферных электронов с отрицательно заряженным до потенциала ф0 диэлектриком Это приводит, с одной стороны, к уменьшению начальной концентрации электронов за счет оттал-

кивания низкоэнерюгическои части первичного потока, с друюй - к торможению высокоэнергетической части и к уменьшению энергии конкретного электрона В натурных условиях наибольшую опасность представляет собой отрицательное заряжение диэлекфических материалов, которое при отсутствии освещенности возникает в областях магнитосферы, содержащих электроны большой энергии, но малой плотностью, (т е в приближении так называемой «толстой» плазменной оболочки, когда радиус Дебая \т> мною больше линейною размера сну шика К) В этом случае для расчетов объемной электризации вдоль координаты «х» для плоских диэлектрических матсриатов толщиной Ь, используемых в космических технотогиях, в диссертации получены следующие соотношения

С Р р

ехр(-^)

Е(х) = -

ехр(-то)

Фт

ЛсХР("10хо>

Фт

1-ехр

Фт

Фо =-

ехр(---°) _Фт

т 10 шх°

ехр(-

■ и

ехр

Р^Р.

^(х

ХоР)

Фг

(6)

(7)

Jeo ехр(——) -ток электронов (А/см2), достигающих поверхность заряженного до Ф

потенциала ф0 диэлектрика, ефт - тепловая энергия (э!3),;ео - невозмущенный ток,

^-Шх'

® (В)

а(х) = от 10 " - радиационная проводимость

диэлектрика, обусловленная термализованпыми электронами, СТТ - темновая (собственная) электропроводность

(Ом/см), х0 =(1/р)[-(фтД) 1п(1-Т12;)],",-глубина, на которой в диэлектрике токи проводимости и электрическое поле обращаются в нуль (см), р-плотность диэлектрика (г/см3), т^— суммарный коэффициент вторичной эмиссия, I;, р- численные коэффициенты, получаемые при аппроксимации экспериментальной зависимости длины пробега электронов от энергии в полимерных материалах, используемых в космической технологии Действие кулоновского отталкивания в выражениях (6) и (7) учитываются множитетем ехр(-ф0/фт), действие токов вторичной эмиссии -ехр(-^рРх^/фт)

На рис 1 приведены рассчитанные профили распределений Е(х) и ф(х) в кварцевом стекле и пленке ПЭТФ (лавсан) Кварцевые стекла используются в качестве защитных покрытии на солнечных батареях, а материал ПЭТФ при изготовлении ЭВТИ (электровакуумная теплоизоляция) Электроны магниюсферной плазмы имеют температуру в фт = 50(кэВ), а шк .)ое=109(А см-2) Причем величина гока, посгунающего на диэлектрик, лежит в области максимальных значений, которые соответствуют сильным магнитосферным возмущениям Параметры диэлектриков кварцевое стекло - плотность р=2 7(г см'3), теневая проводимость ст,=4 10"'6 (Ом см"1), лавсан - р =14(1 см-3), от=4 10~'6 (Ом см-1) Показатель степени поправочного коэффициента для радиационной проводимости Д = 2, суммарный коэффи-

Рис 1 Профили распределения потенциала ф(х) и электрического поля Е(х) в объеме диэлектрика

циснт вторичной эмиссии 41=0.3. Толщина кварцевого стекла и материала из лавсана Ь =

]20цт = =1.2-10~2см.

В заключении третьей главы на основе данных, полученных спутником ОГО-3 па восстановительной фазе (9-13 июля 1966г.) конкретной магнитной бури, на примере геостационарного спутника рассмотрено практическое применение построенных в данной главе моделей (рис.2). Орбита геостационарного КА располагается вне области максимума кольцевого тока

Перед магйитной бурей 7 июля

(хЮ7) Ясмг*с*ср)"'

После бури / \ Перёд магнитной

13 июля \ бу!^ей 7 июля

Рис.2. Динамика расширения магнитосферы на восстановительной фазе магнитной бури и соответствующие изменения в распределении потоков протонов, а) в начале фазы восстановления (9 июля); Ь) в конце фазы восстановления (13 июля). Тремя звездочками на рис.Ь) отмечены максимальные значения потоков для L=6 в соответствующие дни восстановительной фазы: верхняя звездочка соответствует 9 июля; нижняя - 11 июля; средняя - 13 июля.

(L ~ 3-4), поэтому здесь напряженность магнитного поля на главной фазе развития магнитной бури будет возрастать, а на восстановительной фазе уменьшаться до исходных значений опорного (диполыюго) поля. При этом из условия постоянства дополнительных инвариантных форм RfH и Rv следует, что на главной фазе заряженные частицы будут двигаться к Земле (к токовому кольцу) и при этом нагреваться, а на восстановительной фазе удаляются от Земли (от кольца) и остывать. Во внутренней области кольцевого тока всё происходит наоборот; на главной фазе магнитной бури частицы начинают двигаться от Земли к формирующемуся токовому слою и остывают, а на восстановительной фазе - от слоя, т.е к Земле, и нагреваются. Следовательно, на главной фазе магнитной бури в области по обе стороны от расположения максимума протонного пояса (L —3.5) складываются естественные условия, при которых за счет встречного к нему движения заряженных частиц начинается формироваться и усиливаться кольцевой ток. В то же время на восстановительной фазе заряженные частицы двигаются от кольцевого тока, вызывая тем самым его ослабление. Одновременно с этим становится ясен механизм появления в плазмосфере на главной фазе магнитной бури холодных протонов и, наоборот, их нагрев на восстановительной фазе. Следует отметить, что для низкоэнергичных частиц плазмосферы ий~(сЕ/1Г)»иф и на временах развития магнитной бури о третьем инварианте, а тем более его сохранении, говорить не приходится. Поэтому R2II и Rv уже не равны соответствующим константам, а имеют более сложную функциональную зависимость. Однако общий характер зависимости сохраняется; с уменьшением поля Н радиус R растет, а скорость v (энергия) уменьшается; с увеличением Н радиус R уменьшается и скорость (энергия) возрастает.

Четвертая глава посвящена исследованию изменений в параметрах ионосферы в сейсмоактивных регионах. Одновременно с ре1улярными изменениями параметров ионосферы, соответствующих их суточному и сезонному ходу и задаваемых активностью Солнца, в сейсмоактивных регионах Земли наблюдается также и аномальное повеление в распределении характеристик ионосферы, таких как критические частоты f„F2, f„E и f„Es, действующие высоты h', в эффектах проявления турбулизации слоев F2 и Es (F-spread и Es -spread) и т.д. Ионосфера, в общей последовательности геосферных оболочек, представляет собой достаточно чувстви-

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

Рис.3. Результаты вейвлет-рздложения критической частоты Г„Р2 для Камчатки и Магадана за 2002г

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

Па основе новой математической конструкции вейвлет-преобразования была исследована динамика ионосферного процесса, выделены основные характерные для нее компоненты. Сопоставление обрабатываемых рядов с данными сейсмического каталога и сравнение их с асейсмической зоной за данный период позволило выявить разномасштабные аномальные эффекты в ионосфере, которые предшествуют наиболее сильным землетрясениям в районе П,-Камчатского в радиусе Я~200км и к >12 (М>5). Процесс обработки данных базировался на основе непрерывного и дискретного вейвлет-прёобразования.

Непрерывной обработке были подвергнуты часовые данные критической частоты Г0Р2 за 2000-2002гг. На рис. 3 приведены результаты преобразования за 2002г. для Магадана. Стрелками показаны моменты возникновения землетрясений на Камчатке рассматриваемого класса. В 2002г. произошло несколько сейсмических событий, как в весенний, так и в зимний периоды. Из результатов обработки видно (рис.3), что накануне всех событий наблюдаются многомасштабные аномалии, обнаруживающие себя, начиная с 80-го масштабного уровня. Ввиду проявления данных аномалий на более низких частотах в отличие от предыдущих лет, из сигнала была выделена аппроксимирующая компонента 6-го масштаб-

ною уровня дискретного ьейвлеч-разложения и дополнительно исследована. Результаты разложения данной компоненты показаны на рис.4. Ярко проявившие себя в этой компоненте аномальные эффекты наблюдаются в весенний период, начиная с 60-го масштабного уровня, а в осенний - на 15-43 масштабных уровнях, выделены па рис.3 пунктирной линией.

Рис.4. Результаты непрерывного вейвлет-разложения аппроксимирующей компоненты 6-го масштабного уровня дискретного вейвлет-разложения критической частоты СР2 для Камчатки за 2002г

Дискретным разложением показано, что наилучшая аппроксимация временного ряда Г,1*2 на основе авторегрсссиопной модели спокойного периода относительно сейсмической активности режима (отсутствие событий с к>12,5) достигается на четвертом уровне разложения. Сравнение с аналогичным разложением, но для сейсмоактивного периода, показало невозможность аппроксимации и, как следствие, модельного представления на этом уровне разложения. Это говорит о значительном усложнении временного ряда 1"„Р2 и влияния процессов литосферного происхождения на параметры ионосферы.

В заключении четвертой главы статистическими методами изучены явления в ионосфере, связанные с землетрясениями, имеющими глубину не более 50км и магнитуду М>4. На основе исследования 15-минутных данных АИС с. Паратунка (<р = 52,97°N. /1 = 158,25°Е; п-ов Камчатка) анализировались ночные Е^-вргеаё, связанные с турбулизацией слоя Е,. Выбор довольно слабых землетрясений определялся тем, что спорадические слои Ея располагаются на высоте Н~ 100км, и сейсмо-иносферные эффекты должны проявляться в большей мере от слабых событий на этих высотах, чем на высотах слоя Р2 (~300км). Показано, что даже перед слабыми ЗТ за 1-2 дня турбулизация слоя Е, усиливается, и число эффектов Е^-эргеас! увеличивается. Максимальный рост эффекта происходит за 2 часа до местной полуночи. Далее па основе данных японской ионосферной станции Акита (Токио) были проанализированы вариации вероятности появления турбулизации слоя Р (Р-эргеас!) перед ЗТ с М>5.5. Н<80км. Обнаружено, что перед ЗТ вероятность появления Р-зргеас1 уменьшается. Статистическая достоверность увеличения числа Е,-$ргеас1 перед полночью и уменьшения вероятности возникновения Р-хргеас! накануне ЗТ составляет более 0.95.

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

вероятностных меюдов для адекватного описания сейсмическою процесса и сильных землетрясении (Гуфельд И Л , Гусев ГА) Следовательно, необходимо отразить степень случайности, иначе говоря, меру возможности появления юю или иною результата В этой связи необходим математический аппарат, позволяющий в цифре отобразить динамику сейсмического режима При теоретико-множественном представлении сейсмическое явление рассматривается как элементарное событие со, в пространстве элементарных событий А В свою очередь, каждое событие О), характеризуется системой случайных величин энергетическим классом к, широтои ф, долготой^., глубиной Ь и временем I При этом делается естественное предположение о независимости случайных сейсмических исходов от системы случайных величин Время, являющееся неограниченной величиной и не обладающее математическим ожиданием, из системы случайных величин исключено, что соответствует переходу к сейсмосгагическому режиму на рассматриваемом временном интервале Сейсмичность всего региона (или его выбранной части) рассматривается как полная группа событии и описывается в виде распределений условных и безусловных вероятностей Р в частотном представлении Случайные события определяются как различные комбинации системы случайных величин к,фД,И в множестве

Р Такой подход позволяет представить модель сейсмических событий п-ова Камчатка за период инструментальных наблюдений как вероятностное пространство трех объектов {Й,Р,Р} На основе этой модели в качестве примера приводятся расчеты распределений вероятностей для различных случайных событий

Если закон распределения системы случайных величин задан в аналитической форме посредством функции распределения Р(фД,11,к) или ее плотности к), то с вероятно-

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

Для непрерывных величин, определяющих сейсмическое событие, плотность распределения можно представить и как мноюмерную производную от Р(фД,И,к), и как произведение условных и безусловных функций Г следующим соотношением

Г(фД,Ь,к) = Э1Р(ф,Х,Ь,к)/ЭфЭХ.ЭЬЭк = Г(ф) Г(%) Г(Ь]фД)Г(к|А.,ф,Ь) (8)

В (8) приняты обозначения Дф) - безусловная плотность распределения сейсмических событии в зависимости от ф, Г(Х,1<р) - плотность распределения сейсмических событии по X при условии, что их широта равна ф, ¡"(Ь'фД) - плотность распределения ссйсмичсских событий по Ь при условии, что их широта и долгота соответственно равны ф и X, Г(к|\,ф 1)) - плотность распределения сейсмических событий по к при условии, что их долгота, широта и глубина соответственно равны X, ф и Ь Зная аналитический вид плотности распределения (9), можно вычислить вероятность того, что сейсмическое событие попадет в заданные интервалы по широте Д<р, = ф, — ф|_1, долготе АХ, =Х1 - X ^, глубине ЛЬ „, = И т - Ь т_[ и по энергетическому классу Дк„ = кп -к„_,

ф2 Ь2 к,

Р(Дф,,ДХрЛ11п,Лкп)= |с]ф|г(к,ф,X,Ь)с1к-

Ф, ь, к,

=Р(ф, Д,, Ь„, к„) - Р(фм Дн, Ьт_|, к„_,) = ^

=Р(Дф, )Р(Д>.,|Аф, )Р(ДЬ т\АХ,, Дф, )Р( Дк „ | ДЬ т, АХ,, Дф,),

где д.тип-индексы, соответствующих интервалов случайных величин Второе равенство в (9) написано согласно определению функции Г(фД,11,к) Выражения (8-9), являясь функциями четырех переменных, не могут быть представлены графически, хотя численные значения Р(Дк,,ДфрДА,т,Д11п) на основе каталога данных сейсмических событий в часюгном пред-

Рис 5 I иаснрамма и стлажснпый 1рафик плотности распределения ^к) для сейсмических событий гипоцентры, которых распределены на площади 5](а) Плотность распределения ((к) в двойном логарифмическом масштабе (Ь) По оси абсцисс(к) цифрами 01 мечены середины шпервалов (к,± 0 5)

С[аалении вычисляются просю Подобным образом вычисляются и безусловные законы распределения для всех случайных величин к,ф,Х,Ь, а также различные комбинации для условных законов распрелеления от этих переменных Статистическая обработка каталога по формуле (9) даст возмол

вопрос о вероятности того, что сейсмическое сооытие произошло в том или ином заданном ишерваде геофафичсскич координат, глубины и энертешческою класса и в заданный ишервал времени, но и получить численные значения ступенчатой функции распределения НДф, АХ, ДЬ, Дк) Положительный ответ па этот вопрос важен потому, что, как установлено эксперимешально, с ростом чис ш собьний п и уменьшением интервала Д ошосшельная час-юта V стремится к своему математическому аналогу Р, а Р(Дф, ДХ.ДИ.Дк) - кусюйчивому непрерывному распределению 1 (ф,X, Ь, к) Иначе говоря, функция Б определяет потенциальные возможности сейсмического режима произвольной области в терминах вероятности Этот вывод естественно распространить и на численные значения вероятностей, определяющих возможности возникновения будущих собьний, для коюрых координаты гипоцентров и энергетические классы попадают в заданные интервалы, те выполнить прогноз При таком подходе каталог землетрясений, представленный в виде математического объекта трех элементов {П. Р,Р), является опорной моделью, на фоне которой можно отслеживать изменения в сейсмической активности, т е регистрировать вариации в распределении средних вероятностей в локальных областях сейсмического региона за различные интервалы времени Следовательно, в рамках предлагаемого подхода прогноз, подразумевающий ответ на вопросы "где" и "какой ишенсивносш" ожидаегся сейсмическое событие СО,, представляет собой предсказание вероятности его возникновения Р в пространстве элементарных событий £2 Причем 12 и средиие вероятности Р постоянно уточняются за счет естественного увеличения числа случайных событий

На примере закона повторяемости представленного в виде плотности распределения Г(к) для представительного класса к^, >9, расчетами показано, что очищенный от афтершо-ков каталог камчатских землетрясений для заданного интервала времени, но для различных интервалов координат однороден по пространству и тождественен плотности распределения, соответствующей неочищенному каталогу сейсмических событий Отсюда следует справедливость закона повторяемое ги в вероятностном предстаалении и для распределения по энергетическим классам афтершоков Статистическая интерпретация закона повторяемости позволяет записать его в единой форме для произвольного сейсмоактивного объема, минимальность которою задаегся минима1ьным числом сейсмических событий, определяющих устойчивость средних частот(вероятноегей)(рис 5а)

Г(к1-) = Г(к;,2> 1(ГТ|:1к'2 к'!>, (10)

пли

^(к12)=1ЕГ(к-)-712(к'2-С) (10)

Выражения (10) и (10 ) записаны для следующих диапазонов изменения к 95<к,<125 с у1 (к), =9 5) и 12 5 < к, < 16 су, (к2 = 12 5) При явном подобии с известным уравнением графика повторяемости огличие выражения (10) заключается в гам, чю в нем отсутствует зависимость от сейсмической активности А0 Эта величина характеризует число событии единицы поверхности в интервате энертетического класса к,± 05 и изменяется в зависимости от рассматриваемой площади В то же время исходным в (10) является плотность распределения 1'(кД

95 105115115 155145155

| —51 N.1-21845

9! 105 11 5 125 13 5 14 5 15 5

кттость не только ответить на

которая, как показано в диссертации, не зависит от рассматриваемой площади. В двойном логарифмическом масштабе выражение (10") представлено на рис.5Ь. Излом графика приходится на к=12.5.

Нормирование интеграла от функции распределения А^к) на единицу, дает возможность связать углы наклона у, (к'<12.5)и у2 (к2> 12.5) в законе повторяемости в единое выражение:

f(kö> ю^

у, In 10

1-10

-YiCkJ-ki+^j)

= 1.

(11)

у 2 In 10

На рис.6a) для эпицентров, расположенных на площадях Si, S2, S3 произошедших за временной интервал Т (с 01.01.62 г. по 31.12.99 г.), изображены для начальных значений к0= 9.5 сглаженные гистограммы плотностей распределения f(k), соответствующие неочищенному (штриховые линии) и очищенному от афтершоков каталогов. (Данные по землетрясениям Камчатского региона были предоставлены КФ ГС РАН (г. Петропавловск-Камчатский). Этот каталог за период Т с эпицентрами, расположенными на площади S<, был очищен от афтершоков сильных камчатских событий (алгоритм Молчан Г.М., Дмитриевой O.E., программа В.Б.Смирнова; данные предоставлены В.А. Салтыковым)). При Äk = 1 вероятности Р(к) численно совпадают с f (к). Несмотря на бессистемное, имеющее хаотический характер, поведение индивидуальных результатов временного ряда, связанного с попаданием случайных значений энергетического класса в конкретный интервал kj± 0.5, выявляется явная устойчивость соответствующих частот. Для заданного минимального значения ко можно отметить статистическое совпадение графиков независимо от рассматриваемой площади и условия однородности (или нет) каталога.

Sl:No=l6735 SI;N о=1)571 S3;N 0-1121 S1;N н=21845 S2;N и- 14837 S3:N H=1270

No - 1 6735 N o=1 1571 N o=1 121 N и - 2 1 8 4 5 N н 14 8 5 7 N н = I 2 7 0

.1 53 - 168' п л.) .. 1 5б°- 1 63" и

9.5 10.5 11.5 12.5 13.5 14.5 15.5

Рис.6. Плотность распределения ^к) закона повторяемости землетрясений очищенного и неочищенного от афтершоков каталогов, гипоцентры которых расположены на площадях Б], Я, и (а). Интегральная функция распределения Р(к) закона повторяемости для эпицентров сейсмических событий очищенного и неочищенного каталогов на площадях й|, Эг и Бз (Ь). N0 — число событий очищенного и Н, - число событий неочищенного каталогов, ко - начальное значение энергетического класса.

Незначительные отличия возника- " >« "5 12.5 13.5 ы.з 15.5

ют для графика повторяемости с наименьшей площадью в-,, для которой число сейсмических событий более чем на порядок меньше, чем у площадей в, и . На рис.бЬ) представлены соответ- „

ствующие функции распределения Р(к). Из рис. 6а,Ь) следует, что независимо от площади все графи- „ ки практически совпадают. Причем, чем больше число рассматриваемых событий N. тем лучше это

0.01..

совпадение, что находится в полном соответствии с предельными

теоремами теории вероятностей. ^.Усредненные но площадям, 82,83 распределения 1вР(к)

очищенного и неочищенного от афтершоков каталогов с ко = о.5, 9.3, Следовательно, МОЖНО говорить О ю.5 и 11.5 (пунктирные линии). Сплошные линии пример усреднения

статистическом описании закона ск0=9.5и 10.5 по площадям ,Б2,8з отдельно для очищенного и

неочищенного каталогов. 20

повторяемое™ сейсмических событий в вероятностном представлении на основе инвариантных уравнений (10) или (10') безотносительно от рассматриваемой площади Это является важной отличительной особенностью в сравнении с графиками повюряемости, в основе коюрых лежит понятие сейсмической активности единицы поверхности А0, которая изменяется в зависимости от рассматриваемой площади При таком подходе закон повторяемости получает вероятностную интерпретацию, не зависит от рассматриваемой площади и характеризует свойства среды при ее разрушении На рис 7 представлены в логарифмическом масштабе усредненные графики распределения Г(к) для очищенного и неочищенного каталогов для ко = 8 5, 9 5, 10 5 и 115 Обращает на себя внимание практическое совпадение значений 1(к) для разных к0 Расчет показывает, что средине значения вероятностей для различных начальных к0, полученные на основе обработки очищенного и не очищенного каталогов землетрясений дтя 8,»52>5з»54(Дк=1), практически одинаковые Г^3,з4 (1^=8 5)=0 648,33(К> = 9 5) = =0 662, БзОч, = Ю 5) = 0 674, з2 ЭдСко = 11 5) = 0 678 (выбор и 84(52° - 52 5° с ш , 159° - 159 5° в д) определяется тем, что они соответствуют площадям надежной регнетра-ции событий с к Й8 5) Из (11) становится очевидным факт подобия графиков на рис 7 с практически одинаковыми, полученными расчетным путем начальными значениями 1§Р(к'0) Действительно, так как при интегрировании плотности распределения 1"(к) от фиксированного значения к<, до бесконечности и нормировании интеграла на единицу, смещение нижнего преде та на конечное ±Ак не приведет к изменению начальных значений АСк'оУ-РОЛ) В диссертации показано, что вероятностная форма закона повторяемости (10) легко сводится к эмпирическому закону Гугенбсрга-Рнхтера ^ N =А - ЬМ (М - магннтуда события) Причем коэффициенты Л и Ь легко рассчитываются Л - логарифм числа сейсмических событий И, попадающих в начальный интервал 0 335 (ко±1) для заданной площади Б и интервала времени Т, Ь - произведение коэффициента наклона у в (11) на множитель «1 5», характеризующего связь для камчатских землетрясений класса к е магнитудой М (к = 4 6+1 5М (Феоотов СА )) Расчет коэффициентов А и Ь для участков наклона о различными у позволяет представить закон повторяемости двумя соотношениями, находящимися в хорошем согласии с известными выражениями Гутсн-берга-Рихтера,

^Ы(М) = 6 397 - 0 684М (М<5 27,к„<12 5), (12)

1ё1ч[(М) = 7 68 - 0 925М (М>5 27, к0>12 5) (13)

Теоретико-множественный подход дает возможность вычислить плотность распределения сейсмических событий по глубине Г(Ь) для заданной площади в заданный интервал времени и на основе веивлет-анализа исследовать вариации в распределении вероятностей возникновения сейсмических событий по глубине и по времени, например, перед Кроноцким землетрясением Следуя общим положениям работ Соболева Г А, было выбрано временное окно ДТ =5 лет и с шагом Д8=1 год осуществлялось скольжение по каталогу однородных по представительному энергетическому классу (К > 9) событии с 05 12 85 г по 05 12 97 г до момента главного толчка В качест-ве исследуемой области брался параллелепипед с основанием 100x100 км и высотой Н=100 км Центр основания имел координаты [54 2° - 162 2°] и совпадал со средней частью афтершоковой области На основе каталога сейсмических событии п-ова Камчатки было рассчитано распределение вероятностей Р(ДЬ), характеризующих попадания сейсмических событий в заданные интервалы глубин ДЬ=1км до Н=100км за период с 18 07 77 по 05 12 97 г Это распределение выбиралось в качестве опорного распределения Р(ДЬ), для которого рассчитывались соответствующие вейвлет-коэффициелты и с которыми сравнивались вейвлет^коэффициенты соответствующих распределений вероятностей Р(АЬ), рассчитанных для соответствующих временных окон с ДЬ=1км Число событий, попадаемых в каждое временное окно, изменялось от 115 (1988-1993 гг) до 165 (1992-1997 гг) Общее число событий N,7 97=501

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

(Гу/ОМ = И —гдеЧ'Ь - базисный вейвлет,

Г(I) - исследуемый сигнал, а, Ь Е Я, а ^ 0 Известно, что одним из основных преимуществ этого метода является его способность к очень точному анализу локальных свойств в поведении сигнала Кроме того, учитывая зиачи-

Рис.8. (Продолжение рисунка и пояснения к нему на следующей странице)

тельную внешнюю схожесть используемых в данной конструкции базисных функций с локальными областями полученного сигнала (наличие особенностей пикообразной формы, характеризующееся широким спектром частот), предлагаемая методика является наиболее подходящей для решения поставленной задачи. Разложение производилось до 32-го масштабного уровня. Поскольку энергия коэффициентов вейвлет-преобразования дает распределение энергии процесса по масштабам, был выполнен расчет суммы вейалет-коэффициентов по всем масштабным уровням, характеризующей распределение энергии исследуемого сигнала по глубине (рис.8 по центру):

1=1

Результаты обработки общей статистики (1977-1997г.г.), показывают, что наибольшая плотность распределения вероятностей возникновения анализируемых сейсмических событий наблюдается на глубинах 20-40км (рис.8). Анализируя результаты преобразований, вычисляемые по временным окнам, видно, что на малых глубинах (от 1 до 15км), начиная с 1993г., происходит некоторое перераспределение вероятностей возникновения сейсмических событий. На это же указывает результат расчета суммы вейвлет-коэффициентов соответствующих окрестностей. Сопоставление полученной аномалии с данными каталога сейсмических событий за рассматриваемый временной период дает основание предполагать, что полученный аномальный эффект вызван событиями 1992-1993гг., произошедшими южнее Кроноцкого (выделено пунктирными линиями и отмечено стрелками). Расчеты показали, что наблюдаются существенные вариации в распределениях вероятностей возникновения сейсмических событий по глубине в анализируемых временных интервалах на глубинах 4-6км и 18-20км. Причем за один год непосредственно перед Кроноцким событием (окно 1992-1997г.г.) постепенное нарастание интенсивности событий на глубинах 18-20км (1991-1996г.г.) сменяется резким увеличением интенси-1987-1992гг.

1988-1993гг.

1989-1994гг.

1991-1996гг.

Рис.8. Результаты вейвлет-преобразоваиия рядов распределения вероятностей возникновения сейсмических событий по глубине (слева); суммарные значения коэффициентов разложения полученных масштабных уровней по глубине (по центру); исходные данные (справа).

вности событий на глубинах 4-6км и достигает своего максимального значения в последний год накануне события.

Завершается глава описанием качественной модели формирования сейсмического события в очаговой зоне, которая рассматривается как открытая физическая система, непрерыв-

1992-1997 гг.

но обменивающаяся с окружающей средой энергией, импульсом, энтропией и т д Произвольная открытая физическая система постоянно находится под действием малых возмущений Исходя из этого, можно предположить, что ЗТ есть результат проявления некоторой совокупности процессов в литосфере являющейся нелинейной неустойчивой системой и которая находится под действием поля фона внешних возмущений Данный подход позволяет рассматривать сейсмическое событие с точки зрения теории катастроф и применить к нему методы нелинейной термодинамики необратимых процессов, развитые И Пригожиным

В заключении сформулированы основные результаты диссертации

1 Доказана эквивалентность канонической и потоковой форм представления третьего адиабатического инварианта и предложен возможный механизм его нарушения, обусловленный кумулятивными эффектами ускорения радиального дрейфа Показано, что физический смысл третьего адиабатического инварианта проявляется в переменных полях Доказано, что точность сохранения третьего адиабатического инварианта в аксиальных полях различна для случаев экваториального и неэкваториалыюго движения и не совпадаете точностью дрейфового приближения, а соответствует самостоятельному приближению

2 Найдены дополнительные адиабатические формы R2H~ const и Rv~ const в произвольных аксиальных магнитных полях, соответствующие случаю Ф~ const На основе сохранения трех адиабатических инвариантов построена и обоснована динамическая система координат (L, a°,V) которая переводит функцию распределения частиц в реальном (деформированном) поле путем адиабатической вариации магнитного поля в соответствующую функцию для невозмущенного (опорного) поля Это дает возможность судить не только о степени вмо-роженностн плазмы в переменное магнитное поле, но и об инжекции новых частиц в геомагнитную ловушку Для расчета динамики восстановительной фазы конкретной магнитной бури в области орбиты геостационарного КА применены на практике элементы системы координат

(L,a°,V)

3 Решена самосогласованная задача, в которой по распределению тока термализован-ных электронов в диэлектриках находится распределение электрического поля, поддерживающего заданный ток На этой основе построена модель электризации позволяющая рассчитать максимальные уровни отрицательного заряжения, которое диэлектрики в реальных условиях эксплуатации могут приобрести, и выработать соответствующие методы защиты на стадии проектирования Выполнены практические расчеты электризации конкретных диэлектриков, используемых в космической технике для высокоэллиптических и геостационарных КА

4 Определены наиболее опасные участки траектории высокоэллиптических КА, пересекающие силовые линии магнитного поля Земли, проектирующиеся в авроральный овал В период геомагнитных возмущений вдоль этих линий движутся потоки высокотемпературной плазмы, которые приводят к объемному и дифференциальному заряжению КА Эта информация позволяет прогнозировать сбои и отказы на конкретных КА, которые в периоды геомагнитных возмущений работают на этих участках орбиты

5 Впервые показано, что появление в плазмосфсре на главной фазе магнитной бури холодных протонов, поднимающихся вверх (от Земли) обусловлено электрическим дрейфом в

скрещенных вихревом электрическом и у бывающем магнитном полях (rotE = —) При

с at

этом протоны, теряя энергию, остывают На восстановительной фазе магнитной бури магнитное поле возрастает, вихревое электрическое поле меняет направление, и протоны со скоростью электрического дрейфа двигаются уже вниз (к Земле) и, приобретая энергию, разогреваются

6 На основе современной математической конструкции вейвлет-преобразования совмещенной с известной моделью авторстрсссии - проинтегрированного скользящего среднего была исследована динамика ионосферного процесса, выделены основные характерные для нее компоненты Сопоставление обрабатываемых рядов с данными сейсмического каталога и сравнение их с асейсмической зоной (г Магадан) за данный период позволило выявить разномасштабные аномальные эффекты в ионосфере которые предшествуют наиболее сильным землетрясениям в районе Петропавловска-Камчатского в радиусе R~200km и к> 12 (М>5) Дискретным разложением показано, что наилучшая аппроксимация временного ряда f0F2 на основе авторегрессионной модели в сейсмически спокойный период (отсутствие событий с к>12,5) достигается на четвертом уровне разложения Сравнение с аналогичным разложением, но для сейсмоактивного периода, показала невозможность аппроксимации и, как следствие модсльно-

го представления на этом уровне разложения Это говорит о значительном изменении временного ряда f„F2 и влиянии процессов литосферного происхождения на параметры ионосферы

7 На основе статистическою анализа нарамефов ионосферы спорадических слоев Es и F2, подверженных диффузионному расширению, так называемые Es-spread и F-spread, обнаружено, что за 1-2 дня до землетрясения (М>4) частота появления Es-spread увеличивается в преднолуночные часы, т е турбулизация плазмы спорадических Es-елоев усиливается, а турбу-лизацня F-слоя за неделю уменьшается сеиемического события Статистическая достоверность сое la и л мет Ü 95 для всех обнаруженных эффектов

8 Пред южен теоретико-множественный подход к описанию ЗТ произвольною сейсмоактивного региона, позволяющий представить модель сейсмических событии п-ова Камчатка за период инструментальных наблюдений как вероятностное пространство трех обьек-тов{П,Т, Р}, и на основе обработки соответствующего каталога землетрясении вычистять распределения событий и в терминах случайности исследовать различные стороны сейсмического режима

9 Получено вероятное гное представление закона повторяемости есйемических событий, которое позволяет представить его в инвариантной форме (б 5 1а) или (6 5 lb) независимо от рассматриваемой площади и показана его связь с законом Гутенбсрга-Рнхтсра Получены рекуррентные соотношения для определения по каталогу событий Камчатскою региона периоды возникновения катастрофических землетрясений

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

Основные результаты диссертации представлены в етедующих работах

1 Богданов В В, Плетнев В Д К вопросу о точности сохранения третьего адиабатическою инварианта движения заряженной частицы в аксиально-симметричных полях I//Космические исследования 1972 Т 10 Вып 3 С 358-367

2 Богданов В В, Плетнев В Д К вопросу о точности сохранения третьего адиабатического инварианта двилсения заряженной частицы в аксиально-симметричных полях II //Космические исследования 1972 Т 10 Вып 4 С 528-531

3 Богданов В В, Плетнев В Д К вопросу о динамической системе координат в геомагнитной лову шке // Космические исследования 1974 Г 12 Вып 3 С 380-386

4 Богданов ВВ, Киселев МП, Абгарян АЛ О качественной теории измерительных сис-тем//Метрология 1975, Л"» 10

5 Богданов В В, Киселев МИ, Абгарян АЛ О качественной теории измеритетьных систем с адиабатическими инвариантами // Метрология, 1975, № 11

6 Богданов ВВ К вопросу о точности сохранения третьего адиабатического инварианта движения заряженной частицы в магнитном поле // Оптико-физические измерения - М Изд-во стандартов, 1977 ВНИИОФИ, С 8-10

7 Резников СБ, Богданов В В и другие, указанные в описании Устройство для измерения зарядов электризации диэлектриков Авторское свидетельств» JV« 1589789 1990

8 Богданов В В Вероятностная интерпретация закона повюряемости землетрясений на примере Камчатского региона // Докл АН 2006 Т 408 № 3 с 393-397

9 Liperovskaya Е V, ParrotМ, Bogdanov V V, Meister С-У, Rodhn М У, Liperovsky VA On variations of fj;2 and F-spread before strong earthquakes in Japan // Natural Hazards and Earth System Sciences 2006,6 P 735-739

10 Liperovskaya E V .Meister С-V, Pokhotelov OF, Parrot M, Bogilaitov VV, Vasü'eva NE On E,-spread effects in the ionosphere connected to earthquakes // Natural Hazards and Earth System Sciences 2006,6 P 741-744

11 Богданов ВВ, Геппенер ВВ, Манорикова OB Моделирование нестационарных временных рядов геофизических параметров со сложной структурой - С -Петербург ЛЭТИ 2006 - 108с

12 Богданов В В Динамика магнитосферной плазмы в дрейфовом приближении -Владивосток Дальнаука, 2006 -138 с

13 АКТ реализации в НПО им С А Лавочкина научных результатов докторской диссертации Богданова В В от 26 10 2005т

14 Богданов В В, Павлов АВ «Программа для расчета вероятностей сейсмических событий» М ВНТИЦ, 2007, № гос рег 50200702092 Свидетельство об отраслевой регистрации разработки № 9093 от 10 10 2007

15 Мандрнкова ОВ, Богданов В В Исследование сигнала критической частоты на основе конструкции вейвлет-разложения // Известия ВУЗов Северо-Кавказский регион Естественные науки Науки о Земле Спецвыпуск 2007 С 90-93

16 МандриковаО В, Богданов В В, Полозов Ю Л Моделирование сигналов со сложной структурой на основе конструкции вейвлет-преобразования//Информационные технологии 2008 №4 С 12-19

17 Мапдрикова ОВ, Богданов ВВ Аппроксимация статистического распределения случайной величины на основе конструкции вейвлет-преобразования// Известия ВУЗов Северо-Кавказский регион Технические науки 2008 №3, С 53-57

18 Богданов В В Ориентировочные расчеты параметров зарядки кварцевого обтекателя электричеством (качественное представление) Отчет о НИР/ НПО им С А Лавочкина -М, 1985- 19 с

19 Богданов В В Оценка параметров электризации в натурных условиях системы 141323 в составе изделия 14112 и выбор способов и требований по защите элементов конструкции от воздействия статического электричества Отчет о ПИР/ НПО им С А Лавочкина -N1,1986 - 17 с

20 Богданов В В Экспериментальное определение удетьного и переходного сопротивления конструктивных элементов системы 141323 Отчет о НИР/НПО им С А Лавочкина -М ,1988 - 12 с

21 БассолоВС, Богданов В В, Инжечевская И В , Ильин В Г, ОсиповНК, Трищенко ЛД Моделирование крупномасштабных магнитосферных процессов 1989 Препринт ИКИ АН СССР, Пр-1480 40 с

22 Богданов ВВ, Волк ИП Разработка инженерной методики расчета поверхностной электризации высокоомных диэлектриков Этап I 1989 Отчет о НИР Паратунка, ИКИР ДВО РАН - 51 с

23 Богданов В В, Волк ИП Разработка инженерной методики расчета объемной электризации высокоомных диэлектриков Этап II 1989 ОпетоНИР Паратунка, ИКИР ДВО РАН -42 с

24 Богданов В В, Вож IIП Разработка инженерной методики расчета электризации диэлектрика, учитывающей кулоиовское взаимодействие потоков ионизованных частиц с заряженным диэлектриком Этап III 1989 Паратунка, ИКИР ДВО РАН -29 с

25 Атюнтов ИЛ, Богданов В В, Вершинин Е Ф Исследование и научное обоснование возможности создания информационной системы оперативного глобального прогноза землетрясений с использованием данных от космических и наземных средств измерений 1989 Научно-технических отчет, с Паратунка, ИКИР ДВО РАН - 145 с

26 Амиантов IIН, Балебанов ВМ, Богданов В В, Покровскии ВН. Чесноков АГ Контроль электромагнитных предвестников наземными средствами с целью осуществления прогноза землетрясений 1990 Препринт ИКИ АН СССР, Пр-1709 -24 с

27 Богданов В В Сейсмическое явление как процесс разрушения неустойчивой структуры на фоне внешних возмущений // Солнечно-земные связи и предвестники землетрясений Тез докл II международное совещание 14-19 августа 2004г с Паратунка - II-Камчатский 2004 С 12-13

28 Богданов В В Применение теории катастроф и методов нелинейной термодинамики Пригожина для описания сейсмических явлений // Солнечно-земные связи и предвестники землетрясений Тез докл II международное совещание 14-19 августа 2004г с Паратунка - П-Камчатский 2004 С 14-16

29 Богданов В В, Бузевич А В, ВиницкиП А В, Дружин Г И Купцов А В , Подбельский ПН, Смирнов ЭС, Шевцов БМ О влиянии солнечной активности на атмосферные и сейсмические процессы Камчатки // Комплексные сейсмологические геофизические исследования Камчатки Сборник к 25-летию Камчатской опытно-методической партии ГС РАН - П -Камчатский Камчатский печатный двор, 2004 С 259-278

30 E F Vershtnm, У V Bogdanov, G1 Drichin, К Yumolo Electromagnetic precursors of Kamchatka earthquakes according to the data ot ground complex observations data//lnternational Workshop on Seismo-EIectromagnetics of NASDA (IWSE 2000) Tokyo September 19-22, 2000 P 32

31 Вершинин E Ф, Дружин Г И Богданов В В, Юмошо К Сенсмоэлектроматнитные предвестники землетрясении по комплексным наземным наблюдениям и связанные с ними явления в ионосферно-магннтосферной плазме // Проблемы геодинамики и прогноза землетрясений Первый Российско-японский семинар ИТиГДВОРАН Хабаровск 26-29 сентября 2001 С 170-178

32 Богданов В В Вероятностная интерпретация закона повторяемости сейсмических событий на примере Камчатского региона // Солнечно-земные связи и предвестники землетрясений III международное совещание 16-21 августа 2004г с Паратунка [сб докл ] - П-Камчатский 2004 С 25-30

33 Богданов В В, Мандрикова О В Вариации в распределении вероятностей потока событий по глубине и времени перед сильными землетрясениями (на примере Кроноцкого землетрясения) // Солнечно-земные связи и предвестники землетрясений III международное совещание 16-21 августа 2004г с Паратунка [сб докл ] - П-Камчатский 2004 С 31-35

34 Богданов В В, Мандрикова О В Непараметрическое моделирование ионосферных данных в базисе вейвлетов на примере критической частоты f()F2 // Солнечно-земные связи и предвестники землетрясений III международное совещание 16-21 августа 2004г с Паратунка [сб докл ] - П -Камчатский 2004 С 35-42

35 Богданов В В, Мандрикова О В Исследование временных рядов распределений вероятностей сейсмических событий методом вейвлет-разложения (Кроноцкое землетрясение, 05 12 97, Камчатка) // Междунар конф по мягким вычислениям и измерениям, Санкт- Петербург 2005 Т2 С 162-164

36 Мандрикова О В, Богданов В В, Павлов А В Изменения в распределении вероятностей возникновения землетрясений по глубине накануне Кроноцкого землетрясения (05 12 97 г) // Меж- дунар конф по мягким вычислениям и измерениям, 27-29 нюня 2006г, С-Петербург [сб докл] -С-Петербург,2006 Т 2 С 153-157)

37 Богданов В В, Мандрикова О В, корнев ЕВ, Шумилова АЛ, Конюхов А В Об аномалиях в ионосферных процессах Камчатского региона накануне сейсмических событий // Вестник КГТУ Выпуск 4 -П Камчатский, КГТУ 2005 С 52-56

38 Богданов В В , Мандрикова О В, Поюзов Ю А Выявление аномальных изменений в ионосферных параметрах на основе модельного представления сигналов критической частоты f„F2 // Солнечно-земные связи и предвестники землетрясений IV межд конф, 14-17 авг 2007г, с Паратунка Камч обл [сб докл ] - П-Камчатский ИКИР ДВО РАН, 2007 Ч1 С 256259

39 Богданов В В, Мандрикова О В Анализ изменений в распределении вероятностей сейсмических событий по глубине на основе вейвлет-разложений // Солнечно-земные связи и предвестники землетрясений IV межд конф , 14-17 авг 2007г, с Паратунка Камч обл [сб докл] -П-Камчатский ИКИР ДВО РАН, 2007 41 С 260-263

40 Богданов В В, Павлов А В Изучение изменении доверительных интервалов для вероятностей попадания сейсмических событий в заданные интервалы энергетического класса накануне крупных землетрясений // Солнечно-земные связи и предвестники землетрясений IV межд конф, 14-17 авг 2007г, с Паратунка Камч обл [сб докл] - П-Камчатский ИКИР ДВО РАН, 2007 Ч I С 264-268

41 Липеровская Е В, Парро М, Богданов В В, Мейстер К В, Родкин М В, Липеро&ск ий В А О возмущениях t„F2 в среднеширотной ионосфере перед сильными землетрясениями // Солнечно-земные связи и предвестники землетрясений IV межд конф, 14-17 авг 2007г, с Паратунка Камч обл [сб докл ] - П-Камчатский ИКИР ДВО РАН, 2007 Ч I С 367-372

42 Лтеровская ЕВ, Богданов В В, Липеровский В А , Похотелов OA , С шина А С, Олифиров А В Статистический анализ возмущений околополуденной критической частоты foE ионосферы в связи с землетрясениями по материалам станций Ташкент и Петропавловск-Камчатский // Солнечно-земные связи и предвестники землетрясений IV межд конф, 14-17 авг 2007г , с Паратунка Камч обл [сб докл ] - П-Камчатский ИКИР ДВО РАН, 2007 Ч I С 372-377

43 Липеровская Е В, Богданов В В, Родкин МВ, Мейстер К В, Силина А С, Мандрикова О В Статистический анализ возмущений критической частоты foF2 ионосферы за несколько суток до и после землетрясениями - сопоставление результатов, полученных по мате-

риалам японской станции Токио и станции Петропавловск-Камчатский // Солнечно-земные связи и предвестники землетрясений IV межд конф 14-17 авг 2007г, с Паратунка Камч обл [сб докл ] - П-Камчатс-кий ИКИР ДВО РАН, 2007 ЧI С 378 -384

44 Липероеская Е В, Богданов В В, Родкин М В, Мейстер К В Васильева Н Э, Опи-фиров А В Суточная зависимость возмущений в спорадическом слое Es ионосферы в связи с землетрясениями по материалам станций вертикального зондирования Петропавловск-Камчатский и Токио // Солнечно-земные связи и предвестники землетрясений IV межд конф, 14-17 авг 2007г, с Паратунка Камч обл [сбдокл] - П-Камчатский ИКИР ДВО РАН, 2007 Ч I С 385-390

45 Богданов В В, Мандрикова О В, Пав те А В, Шумилова А Л Вероятностная интерпретация закона Гуттенберга-Рихтера и оценка времен повторения сильных камчатских землетрясений с магнитудой равной 9 и 9 5 баллов // Междунар конф по мягким вычислениям и измерениям, 25-27 июня 2007г, С-Петербург [сб докл] - С-Петербург, 2007 Т2 С 177-181

46 Мандрикова ОВ, Богданов В В Исследование вариаций распределения сейсмических событий по глубине на основе применения конструкции вей влет-преобразования // Междунар конф по мягким вычислениям и измерениям 25-27 июня 2007г, С-Петербург [сб докл] -С-Петербург, 2007 Т2 С 209-212

47 Богданов В В, Мандрикова О В Вейвлет-анализ распределений вероятностей сейсмических событий Камчатского региона накануне Кроноцкого землетрясения // Вестник КГТУ Выпуск4 П-Камчатский, КГТУ 2005 С 56-59

48 Мандрикова О В., Богданов В В, Поюзов ЮА Анализ изменений, происходящих в ионосферных параметрах накануне землетрясении на Камчатке, на основе вейвлет-преобразования сигналов критической частоты // Междунар конф по мягким вычислениям и измерениям, 25-27 июня 2007г , С-Петербург [сб докл] -С-Петербург, 2007 Т2 С 213-216

49 Мандрикова О В Богданов В В Многоструктурные модели геофизических сигналов // Распознавание образов и анализ изображений Новые информационные технологии, VIII межд конф 7-13 окт 2007 г, Йошкар-Ола [сб докл] - Map ГТУ Йошкар-Ола 2007 С 2021

50 Liperovska)>a Е V, Bogdanov V V, Meister С-V, Liperovsky VА , Rodkm М V Day-time variations of foF2 connected to strong earthquakes // XXIV General Assembly of the Int Union of Geodesy and Geophysics Perugia, Italy, July 2-13,2007

51 Богданов В В Динамика заряженных частиц в нестационарной магнитосфере с учетом сохранения трех адиабатических инвариантов//Физика плазмы в солнечной системе, конф 5-8 февраля, 2008, Москва [сб тезисов] -ИКИ РАН, 2008 http //solarwind cosmos ru/txt/tes 2008 pdf С 45

52 Амиантов И Я, Богданов В В, Вершинин Е Ф идр Анализ возможностей создания комплексной системы прогноза землетрясений по данным наземных и космических измерений параметров солнечной активности, магнитосферы и электромагнитного поля в широком диапазоне частот 1990 Отчет о НИР Паратунка, ИКИР ДВО РАН -120 с

53 Амиантов И Н, Богданов В В., Вершинин Е Ф и др Разработка алгоритмов краткосрочных прогнозов аномальных явлений на основе электромагнитных предвестников 1991 Отчет о НИР Паратунка, ИКИР ДВО РАН - 24 с

54 Богданов В В Обнаружение и выделение оперативных электромагнитных предвестников для прогноза сильных землетрясений на Камчатке 1995 Научно-технический договор № 1-95 ИКИР ДВО РАН - 14 с

55 Маричев В Н, Богданов В В, Живетьев ИВ Шевцов Б М Образование аэрозольных слоев в стратосфере под воздействием солнечной активности // Геомагнетизм и аэрономия 2004 Т 44, №6 С 841-848

56 Шевцов Б М, Маричев В Н, Богданов В В, Живетьев И В Образование аэрозольных слоев в стратосфере под воздействием магнитной активности//Солнечно-земная физика межд конф, 20-25 сент 2004 г, Иркутск [тезисы] -Иркутск ИСЗФ СО РАН, 2004 С 66

57 VN Marichev, V V Bogdanov, I V Jivetyev, В M Shevtsov_ Influence of solar activity on formation of aerosol layers in the stratosphere //Proe XI Joint International Symposium "Atmospenc and Ocean Optics Atmospheric Physics" Tomsk, 2004 - P 169

58 VN Marichev, V V Bogdanov, IV Jivetyev, В M Shevtsov "Formation of aerosol layers in the stratosphere under the influence of solar activity according to the Kamchatka observation data" //Proc XI Joint International Symposium "Atmospenc and Ocean Optics Atmospheric Physics" Tomsk, 2004 P 172

59 В M Shevtsov, VN Marichev, VV Bogdanov, I V Jivetyex, V VBychkov About mechanisms of the aerosol layers formation in the stratosphere under the influence of magnetic activ-ity//Proc International Association of Meteorology and Atmospheric Sciences Beijing, China August 2005 A-43

60 Богданов BB, Калинин BH Феноменологическое уравнение бесстолкновитслыюй плазмы в дрейфовом приближении / Применение ультраакустики к исследованию вещества [сб научных тр ]/ред дф-мн ВФ Ноздрев-М ВЗМИ 1981 Вып 31 С 95-106

61 Liperovsky VA , Meister С-V, Liperoxskaya Е V, Davdov VF, Bogdanov V V On the possible influence of radon and aerosol injection on the atmosphere and ionosphere belore earth-quakes//Natural Hazards and Earth System Sciences 2005 V 5 N6 P 783 - 789

62 Jlunepoeckuü BA, Muxatauii BB, Шевцов Б M, УмархоОжаев P M, Богданов В В, Мейстер KB Электрические поля и инфракрасное излучение в аэрозольных облаках перед землетрясениями // Солнечно-земные связи и предвестники землетрясений IV медд конф, 14-17 авг 2007г, с Паратунка Камч обл [сб докл ] П-Камчатский ИКИР ДВО РАН, 2007 Ч I С 391-397

63 Мандрикова О В, Богданов В В, Конюхов А В Оптимизация обучения нейронной сети на основе вейвлет-преобразовання на примере модельного представления ионосферного сигнала // Междунар конф по мягким вычислениям и измерениям, 27-29 июня 2006г, С-Петербург [сб докл] - С-Петербург, 2006 Т1 С 223-226

64 Мандрикова О В Богданов В В, Полозов Ю А Моделирование сигнала критической частоты на основе нейронных сетей с целью поиска аномального поведения // Междунар конф по мягким вычислениям и измерениям, 27-29 июня 2006г , С-Петербург [сб докл ] - С-Петербург, 2006 Т2 С 158-160

65 Богданов В В Применение вариационных методов неравновесной термодинамики для получения диссипативной системы МГД-уравнений в приближении ведущего центра // Солнечно-земные связи и предвестники землетрясении IV межд конф , 14-17 авг 2007, с Паратунка Камч обл [сб докл ] П-Камчатский ИКИР ДВО РАН, 2007 Ч I С 37-42

66 Богданов В В, Шумилова А Л Применение методов вейвлет-разложения для анализа аномального поведения параметров электрического поля перед землетрясением по данным ИСЗ «Demeter» // Достижения в спутниковой океанографии изучение и мониторинг окраинных морей Азии 3-6 окт2007 г, Владивосток [сб материалов] Владивосток Дальна)ка, 2007 С 19-20

67 Богданов В В, Шумилова A J1 Предварительные результаты обработки измерений электромагнитного поля в диапазоне )-20 кГц, зарегистрированных на ИСЗ «Demeter» (Франция) накануне землетрясения в Коряки 20 04 06 // Проблемы комплексного геофизического мониторинга Дальнего Востока России I регион научно-техн конференция 11-17 ноября 2007, Петропавловск-Камчатский [тезисы] -Петропавловск-Камчатский, ГС РАН, 2007 С 111

68 Богданов В В, Павлов А В Поиск влияния литосферных процессов на ионосферный слой методами непрерывного и дискретного вейвлет-преобразования // Проблемы комплексного геофизического мониторинга Дальнего Востока России I регион научно-техн конференция 11-17 ноября 2007г, Петропавловск-Камчатский [тезисы] - Петропавловск-Камчатский, ГС РАН, 2007 С 108

69 Липеровский В А, Мейстер KB, Леперовская ЕВ, Богданов В В, Давыдов ВФ Электрогравитационные эффекты при инжекции радона и наличии аэрозолей в надразломных зонах перед землетрясениями // Актуальные проблемы мониторинга рисков чрезвычайных ситуаций Научно-методическая конференция 11 октября 2006 [сб материалов] Химки, академия гражданской защиты, 2007 С 69-83

Отпечатано ИП Романснко М И Заказ № 128/256-512 Тираж 150 экз Лиц ПД № 14-2А 002305 г Петропавловск-Камчатский, ул Ленина, 46

Текст научной работыДиссертация по наукам о земле, доктора физико-математических наук, Богданов, Вадим Васильевич, с. Паратунка, Камчатский край

поскольку в аксиально-симметричном поле г/рл = const вдоль траектории заряженной частицы (см. выражение (2.3.8)).

В работах [Хесс В., 1972; Hamlin D.A., 1961] показано, что в дипольном магнитном поле ф(а0) является медленно меняющейся функцией экваториального угла а0 между векторами V и Я. Кроме того, условие применимости дрейфового приближения требует, чтобы изменение сс0, связанное с нестационарностью магнитного поля, было незначительным в течение Ткш, то есть 5ТК0Л ~ е за время Ткол. Отсюда, пренебрегая изменением а0 в течение одного полного колебания, получим

1 t /V

= const • (к

Т / \ (АЛ5)

ТКОЛ К)

с точностью приближения среднего дрейфа.

3. Для осреднения к введем угол 0, связанный с периодом выражением

2п

Q = Y~t+<p = Cim» -t+Ф

(А. 16)

Выбирая начальную точку траектории на экваторе, найдем, что ср=0. Очевидно, что в стационарном аксиально-симметричном поле Т^я^сог^ и сШ

= Q .

j, кол

В нестационарном магнитном поле

de = Q +do^t dt кол dt

(А. 17)

(А. 18)

Уравнение дрейфового приближения для dvII/dt имеет вид (в правой части (А.19) дифференцирование для юл производится по дуге силовой линии 8) dv,т

_ yj аи>л , ку„с(Е-§ф)

dt 2юл 5S г Н Дня стационарного магнитного поля

(А.19)

dv„ dt

Пусть

v„ = Vn cos

у; ^л

2юл as

2п

t = v„ cos9,

(А.197)

(А.20)

Информация о работе

Богданов, Вадим Васильевич
доктора физико-математических наук
с. Паратунка, Камчатский край, 2008
ВАК 25.00.29

Диссертация

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

Автореферат

Похожие работы