Реферат
Отчет 89 с., 6 ч., 35 рис., 2 прил.
Ключевые слова – сейсмология, землетрясение, сейсмическая томография, обратные задачи, геоэлектрика, УНЧ-излучение, магнитное пересоединение, магнитосферная суббуря.
Работы, вошедшие в проект, продолжаются в течение длительного времени и выполняются в тесном международном сотрудничестве с партнерами из многих стран (список международных соглашений дан в Приложении). Такое долговременное сотрудничество и возобновление договоров (так же как и заключение новых) является дополнительным подтверждением значимости и успешного выполнения научных исследований.
В соответствии с Техническим заданием и Календарным планом на 2007 год по проекту проводились следующие исследования:
1. Развитие томографических методов восстановления геофизических параметров (руководитель В.Н. Троян) (пункт Технического задания: разработка исследовательских вычислительных программ, позволяющих моделировать волновые поля в слоистых пористых средах и извлекать физические закономерности волновых процессов;)
2. Изучение сейсмоэлектромагнитных процессов, связанных с подготовкой сильных землетрясений в тектонически активных областях (руководитель Н.А. Смирнова) (пункт Технического задания: моделирование сейсмо-электромагнитных процессов)
3. Теоретическое и экспериментальное исследование магнитного пересоединения (руководитель В.С.Семенов) (пункт Технического задания: реализация алгоритма восстановления скорости пересоединения и пересоединившегося магнитного потока по данным о возмущениях магнитного поля на данных системы спутников CLUSTER.)
Кроме того, продолжались работы по традиционным направлениям:
4. Исследование глубинного строения земной коры и верхней мантии по данным сети сейсмостанций (руководитель Т.Б. Яновская )
5. Магнитотеллурические исследования (руководитель А.А. Ковтун);
6. Исследование динамики токового слоя магнитосферы (руководитель Сергеев В.А.)
Объекты исследования проекта в целом это:
внутреннее строение Земли, а именно: анизотропия верхней мантии, электропроводность коры и верхней мантии Фенноскандинавского щита, аномальная электропроводность литосферы на севере Ботнического залива; неоднородные среды с различными упругими свойствами, с горизонтальными и вертикальными градиентами скоростей распространения волн, а также с отражающими границами, на которых скорость терпит разрыв. Кроме того, исследовались естественные ультранизкочастотные (УНЧ) электромагнитные поля Земли и пространственно-временное распределение региональной сейсмичности, способные нести информацию о приближении катастрофических событий (мощных землетрясений);
внешнее строение Земли – структура и эволюция токовых слоев в магнитосфере Земли, взрывная диссипация энергии в космических токовых слоях, параметры, обуславливающие динамическое состояние магнитосферы.
Цель работы - Проект направлен на решение фундаментальных проблем геофизики, состоящих в исследовании свойств естественных электромагнитных полей и упругих колебаний Земли, изучении физической природы их источников, и развитие на этой основе методов исследования глубинного строения земной коры и верхней мантии, способов разведки полезных ископаемых, прогноза землетрясений и диагностики состояния околоземной плазмы.
Методы, применяющиеся для разных задач проекта, очень разнообразны и более подробно описываются в соответствующих разделах проекта. Общими можно назвать методы компьютерного моделирования для решения прямых и обратных задач геофизики, двумерную томографию, контроль численных решений для модельных задач аналитическими методами.
Результаты работы (за 4-ый этап):
1. Выполнено обобщение метода расчета полей волн в слоистых упругих средах на случай упруго-пористых сред, описываемых уравнениями Френкеля-Био, в предположении отличной от нуля вязкости флюида. Учет вязкости жидкой фазы приводит к тому, что обычные волновые уравнения (гиперболического типа), определяющие законы распространения упругих волн, превращаются при наличии значительной вязкости в уравнения, описывающие "диффузию" волн. В вычислительном отношении это приводит к смещению особых точек подынтегральных функций (точек ветвления, полюсов ) в комплексную плоскость.
2. Основываясь на точном решении задачи о волновом поле источника, расположенного в заполненной жидкостью скважине, окруженной пористой средой, были построены диаграммы направленности продольных (быстрой и медленной) и поперечных волн, возбуждаемых в окружающей пористой среде. В разработанных алгоритмах и программах реализованы два механизма взаимодействия жидкости в скважине с пористым флюидом ( "закрытые " и "открытые " поры).
3. Созданный комплекс программ был применен для анализа данных сейсмического эксперимента с месторождения Маллик в Канаде. На основе сопоставления результатов численного моделирования с экспериментальными данными:
- Доказана возможность проведения межскважинного просвечивания при помощи трубных волн в скважине и объемных волн, распространяющихся в пористом резервуаре.
- Предложена методика использования возбуждаемых в модели трубных и объемных волн для обнаружения изменений в резервуаре в процессе добычи.
- Предложена методика контроля координат источников и приемников в процессе проведения межскважинного просвечивания.
- Проведено моделирование сейсмо-электромагнитных процессов для поиска адекватных механизмов генерации УНЧ электромагнитных сигналов в зоне очага землетрясения. Исследованы механизмы генерации и частотно-временные характеристики сейсмо-электромагнитных эмиссий в моделях критической динамики метастабильных тектонических систем. Показано, что включение в развиваемую тектоническую СОК модель параметров электропроводности (например, непроводящие горные породы и проводящие флюиды, заполняющие поры и трещины в залегающих породах, или движение электрических зарядов, образующихся на концах трещин при естественном растрескивании горных пород), приводит к фрактальной динамике электропроводности. Это оправдывает применение нами фрактальных методов при анализе динамики электромагнитных полей Земли в сейсмоактивных регионах и служит свидетельством того, что самоорганизованные критические процессы (СОК процессы) в зонах разломов могут лежать в основе сейсмо-электромагнитных предвестников землетрясений.
- Выполнен фрактальный анализ временных рядов данных по регистрации УНЧ (ультранизкочастотных - f= 0.001-1 Гц) электромагнитных полей Земли в сейсмоактивном (остров Гуам, Тихий Океан, Φm=9° N, Λm=225°) и сейсмоспокойном (дрейфующая станция Северный Полюс-30, Fm=75o N, Lm=172o) регионах на подготовительной фазе мощного землетрясения и большой магнитной бури. Выявлены определенные тенденции в динамике фрактальных характеристик, которые, в случае их подтверждения на независимом материале, могут быть использованы при разработке методов прогноза землетрясений и сильных геомагнитных бурь. Показано, что есть общие закономерности в динамике УНЧ- электромагнитных шумов на подготовительной фазе мощных природных событий – появление фрактальных (коррелированных) флуктуаций со спектром f – 1/b (фликер-шум) в определенном диапазоне частот. Это указывает на усиление пространственно-временных корреляций в критической системе при приближении катастрофических событий. Полученные результаты дают основание для нового направления поиска предвестников землетрясений, основанного на фрактальном анализе вариаций электромагнитных полей Земли.
- В 2007 году проведены детальные исследования импульсного пересоединения на дневной магнитопаузе, в результате которых были впервые получены значения величины скорости пересоединения. Были выбраны два интервала с интенсивной активностью. а) 21 января 2001 г., - 4 случая пересоединения (т.н. FTE события) наблюдались за каспом, восстановленная скорость пересоединения оказалась равной 3.7 +- 0.8 мв/м; б) 14 февраля 2001г. - FTE события наблюдались на высокоширотной магнитопаузе перед каспом, восстановленная скорость пересоединения оказалась равной 2.1 +- 0.4 мв/м для магнитосферных событий и 0.9 +- 0.5 мв/м для FTE событий в переходной области. Во всех случаях длительность импульса пересоединения составляла 1-2 мин. Результаты были сопоставлены с оценками скорости пересоединения по радарным данным в ионосфере. Оказалось, что наши результаты в целом дает более высокие значения скорости пересоединения, что, по-видимому, связано с тем, что радарные измерения дают усредненные по времени значения этого параметра
- Построена модель флэппинга - осцилляций магнитного поля в хвосте магнитосферы, - основанная на неустойчивости двойного градиента: градиента тангенциальной компоненты поля в нормальном направлении и градиента нормальной компоненты вдоль слоя. Оказалось, что если произведение этих градиентов отрицательно, что случается в местах локального утоньшения токового слоя, то слой неустойчив и возмущения нарастают с характерным временем 5-10 минут. Если же произведение этих градиентов положительно, то слой устойчив, по его поверхности бегут волны от центра к флангам со скоростями порядка 50 км/с. Разработанная модель хорошо описывает известные экспериментальные факты флэппинга.
- Разработан метод совместного использования данных о групповых скоростях на трассах эпицентр-станция и о фазовых скоростях на трассах между станциями. Главным моментом в этой методике является определение скоростных разрезов, средних по межстанционным трассам, совместно по данным измерений фазовой скорости и расчетным дисперсионным кривым групповой скорости. Поскольку данных по межстанционным трассам недостаточно для определения латеральных вариаций скорости методом двумерной сейсмической томографии, предлагается использовать полученные средние скоростные разрезы для нахождения начального приближения при инверсии локальных дисперсионных кривых групповых скоростей, получаемых методом двумерной томографии. Разработанный метод применен для уточнения строения верхней мантии в регионе Азии. Данные о фазовых скоростях предоставлены нам институтом земной коры СО РАН. Получены уточненные скоростные разрезы верхней мантии в области, ограниченной 32-56°N, 80-120°E. Оказалось, что обнаруженная ранее 105-градусная граница между областью пониженных и повышенных скоростей проявляется только на глубинах 200-250км. Под Тибетом отсутствует астеносферный слой. Как и в предыдущем исследовании (Yanovskaya, Kozhevnikov, PEPI, 2003), выявилось поднятие к поверхности низкоскоростного вещества в области Байкальского рифта, что подтверждает гипотезу о наличии плюма. Низкие скорости в верхней мантии под Хангайским поднятием указывают на высокие температуры в этой области, что и приводит к поднятию поверхности.
- Модифицированным методом поверхностно-волновой томографии получены новые сведения о особенностях строения верхней мантии Земли в Альпийской зоне по данным фазовых скоростей поверхностных волн. Дисперсионные кривые фазовых скоростей волн Релея и Лява на трассах между станциями получены по разности фазовых спектров от удаленных землетрясений. Скоростные разрезы поперечных волн в верхней мантии до 275 км построены в области, включающей Альпы, Динариды и область Европейской плиты к северу от Альп. Показано, что под Альпами-Динаридами в верхней мантии отсутствует зона пониженной скорости, а на глубинах 100-200 км выявляется область повышенной скорости, совпадающая в плане с северной границей Альпийской зоны. Под Рейнским грабеном на глубинах 100-200 км, наоборот, выделяется зона пониженных скоростей. По характеру распределения скорости в верхней мантии под Альпами делается вывод о несостоятельности недавно предложенной гипотезы о погружении Адриатической плиты под Европу. Данное исследование проводилось в рамках гранта ИНТАС 04-83-381 в кооперации с учеными из Италии.
- Разработан новый численный метод оценки поля поверхностных волн в окрестности каустик, образующихся на больших эпицентральных расстояниях (>90°) за счет латеральной неоднородности и сферичности Земли. Метод имеет то преимущество, что он позволяет рассчитывать волновые поля и в области наложения нескольких каустик. В тех областях, где применим асимптотический метод для анализа поля, результаты, полученные численным методом, хорошо согласуются с результатами асимптотического анализа. Для модели реальной Земли численным методом рассчитаны аномалии амплитуд и фаз на полусфере противоположной эпицентру для периодов волн Релея 40-100 с. Полученные результаты показывают, что при интерпретации реальных волновых полей нельзя пренебрегать аномалиями, обусловленными близостью к каустикам.
- В плане завершения ранее выполненных исследований на Фенноскандинавском щите нами были дополнительно рассмотрены данные магнитотеллурических и магнитовариационных зондирований, полученные по международному проекту ВЕАR, с целью выявления влияния крупной аномалии в электропроводности коры, расположенной вблизи границы между Карельским и Свекофенским геоблоками и простирающейся в северо-западном направлении на сотни километров. Как выяснилось, ее влияние простирается на одну-две сотни километров от оси аномалии и приводит к фиктивному завышению продольной проводимости коры, получаемой при одномерной интерпретации фазовых магнитотеллурических кривых. Возникающая аномалия в магнитовариационных параметрах позволяет более точно оценить истинные параметры аномальной зоны. По результатам исследований подготовлена статья для журнала Вопросы Геофизики и будет сделан доклад на конференции Геокосмос 2008.
- В результате совместной работы компетентной международной группы проведена классификация среднемасштабных транзиентных структур разного типа наблюдающихся в хвосте магнитосферы в активные периоды и сделано заключение, что многообразие структур может быть в основном сведено к 2-3 группам, из которых наиболее важное значение имеет группа явлений создаваемая процессом импульсного локализованного (поперек хвоста) магнитного пересоединения, одним из проявлений которых являются быстрые струйные течения плазмы ориентированные в плоскости силовых линий вдоль хвоста магнитосферы. По материалам совместной работы подготовлен обзор теорий генерации и наблюдений среднемасштабных структур, выполнены исследования топологии магнитных канатов и сравнительный анализ разных методов количественной реконструкции среднемасштабных плазменных структур.
Прогнозные предположения о развитии исследований: