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

Мясников, Андрей Владимирович

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

Автореферат диссертации по теме "Мониторинг состояния магматических структур вулкана Эльбрус по наблюдениям литосферных деформаций Баксанским лазерным интерферометром"

На правах рукописи

005047337

Мясников Андрей Владимирович

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

25.00.10 - геофизика, геофизические методы поисков полезных ископаемых

автореферат

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

2 7 СЕН 2012

Москва 2012

005047337

Работа выполнена в Московском государственном университете им. М.В.Ломоносова, Государственном Астрономическом институте им. П.К. Штернберга

Научный руководитель: доктор физико-математических наук

Милюков Вадим Константинович,

Официальные оппоненты: доктор физико-математических наук

Собисевич Алексей Леонидович, Федеральное государственное бюджетное учреждение науки Института физики Земли им. О.Ю.Шмидта РАН, заведующий лабораторией Прикладной физики и вулканологии

кандидат физико-математических наук Старовойт Олег Евгеньевич Федеральное государственное бюджетное учреждение науки Геофизическая служба РАН (ГС РАН), Зам. директора по научной работе

Ведущая организация: Федеральное государственное бюджетное

учреждение науки Институт геологии рудных месторождений, петрографии, минералогии и геохимии РАН (ИГЕМ РАН)

Защита состоится _2012 года на заседании диссертационного

совета Д 002.001.01 в Федеральном государственном бюджетном учреждении науки Института физики Земли им. О.Ю.Шмидта РАН по адресу: 123995, ГСП-5, Москва Д-242, Б.Грузинская ул., 10, стр. 1.

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

Автореферат разослан_2012 г.

Ученый с екр етарь дис сер тационного • с об е та ^

к.ф.-м.н Р^,__^ О .В. Пилипенко

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

Вулкан Эльбрус расположен на северном склоне Главного Кавказского хребта центральной части Большого Кавказа и соединен с ним горной перемычкой Хотютау [Авдулов, 1994]. Современные геологические исследования, в частности, радиоуглеродного датирование продуктов его извержений, показывают, что такие извержения происходили многократно и в течение длительного времени (около 1 миллиона лет) [Чернышев, 2001]. Эльбрус относится к классу полигенных стратовулканов центрального типа, которые, как правило, распространены над зонами субдукции или конвергентными границами тектонических плит.

Актуальность темы

Для проведения прогнозной оценки возможной активности вулкана Эльбрус необходимо выяснить состояние и динамику вулканических магматических структур. Особенно важно оценить размер и глубину залегания этих образований и решить вопрос о возможном возобновлении активности Эльбруса. Составленные карты хронологии вулканических циклов и различных этапов практически совпадает у разных групп исследователей [Богатиков O.A., Гурбанов А.Г., Мелекесцев И.В. и др., 1999]. О геометрических параметрах внутренней структуры вулкана до сих пор не существует единого мнения и оценки сильно разнятся. Кроме того, в настоящее время нет ни крупномасштабной исследовательской программы по непрерывному мониторингу внутреннего состояния вулканической структуры, что является самым важным моментом в прогнозе вулканической активности, и, в конечном итоге, извержения, которое для данного типа вулкана будет иметь эксплозивный характер. Отсутствуют разработки, которые учитывают резонансные особенности неоднородных структур вулканической постройки. Исследования, связанные с изучением резонансных особенностей магматического очага и магматической камеры могут быть использованы для определения свойств и состояния внутренних структур вулкана.

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

прилегающим территориям, Эльбрус обладает самой крупной гляциологической системой Кавказа, состоящей из 25 ледников. Суммарная физическая поверхность всех ледников с учетом снежно-ледового покрова вершин Эльбруса составляет 139 кв. км, а суммарный объем льда оценивается в 6 куб. км. Наличие обширного снежно-ледового покрова делает вулкан Эльбрус еще более опасным, так как в случае извержений обязательно возникнет катастрофическая ситуация, связанная с быстрым разрушением и таянием ледников. Это приведет к возникновению лахаров и крупномасштабных наводнений в регионе. Причем в момент извержения в горных ущельях потоки воды до нескольких десятков метров. Таким образом, опасности подвергается не только населенные пункты непосредственно расположенные в Приэльбрусье, а практически весь Северный Кавказ в целом, особенно территории в долинах рек Баксан, Кубань и Малка.

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

Цель работы

На территории Северного Кавказа, в 20 км от вулкана Эльбрус, создана уникальная геофизическая лаборатория ГАИШ МГУ. С помощью длиннобазового лазерного интерферометра-деформографа, расположенного в штольне Баксанской Нейтринной обсерватории, ведется непрерывный мониторинг движения земной коры с 1998 года.

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

• Разработка методики выявления собственных колебаний резонансных структур, возбужденных воздействием сейсмической волны. Определение геометрических и физико-механических параметров магматической камеры Эльбруса на основе разработанных моделей магматических структур вулканов.

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

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

• Анализ влияния метеорологических факторов, таких как атмосферное давление и температура, на измерения деформации литосферы, выработка алгоритма по снижению этих факторов.

Научная новизна

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

2. Проведен спектральный анализ литосферных деформаций, возбужденных удаленными землетрясениями с магнитудами М\у 4.5-9.1, за более чем 7-летний период наблюдений. В результате этого анализа в спектрах литосферных деформаций с помощью резонансного метода выявлены региональные резонансные моды, относящиеся к магматическим структурам вулкана Эльбрус. Получены оценки характерного размера, глубины залегания, состава магматических флюидов близповерхностной магматической камеры вулкана Эльбрус.

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

4. Разработан новый адаптивный многопараметрический алгоритм компенсации термоупругих и барических процессов в измерениях литосферных деформаций Баксанским лазерным интерферометром, входными параметрами которого служат измеренные значения атмосферного давления, а также температуры в трех различных узлах прибора.

5. Исследовано влияние метеорологических возмущений на измерение литосферных деформаций Баксанским лазерным интерферометром в условиях подземного размещения в БНО ИЯИ РАН. Показано, что разработанный алгоритм компенсации обеспечивает достаточно эффективное подавление

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

Защищаемые положения

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

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

3. С помощью авторского пакета программ выявлены резонансные моды магматической камеры вулкана Эльбрус, на основе семилетнего ряда наблюдений. Получены оценки геометрических параметров близповерхностной магматической камеры в структуре Эльбрусского вулканического центра, которые согласуются с полученными ранее результатами [Богатиков и др.,2001; Милюков, 2006].

4. Обнаружена динамика изменения добротности резонансных мод близповерхностной магматической камеры вулкана Эльбрус на 7-ем интервале непрерывных наблюдений.

Личный вклад автора

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

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

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

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

4. Анализ литосферных деформаций вызванных сейсмическими событиями, произошедшими за семилетний период наблюдений и зарегистрированными Баксанским интерферометром-деформографом.

5. Практическое применение резонансного метода и обнаружение в сверхнизком диапазоне собственных мод, которые относятся к внутренней структуре вулкана Эльбрус. Выявлено устойчивое уменьшение добротности магматической камеры как акустического резонатора и доказательство практической возможности мониторинга внутреннего состояния вулкана.

6. В течение более чем 20 лет автор принимал активное участие в работе коллектива по созданию Баксанского лазерного интерферометрического комплекса, обеспечению его работы в режиме долговременных наблюдений.

Практическая значимость

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

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

Основное содержание диссертации изложено в работах

1 Буклерский А.В., Карт A.M., Клячко Б.С., Мележников И.В., Милюков В.К., Мясников А.В., Нестеров В.В., Руденко В.Н. Баксаиский лазерный интерферометр // Измерительная техника. 1995. № 10. С. 5.

2 Milyukov V.K., Rudenko V.N., Klyachko B.S., Kart A.M., Myasnikov A.V. Wide-Band Laser Inerferometer for Monitoring the Earth Strains. In: Solid State Lasers (Laser Optics' 98). Proceedings of the International Society for Optical Engineering (SPIE), 1998, V. 368, p.l 16-121.

3 Милюков В.К., Руденко В.Н., Клячко Б.С., Карт A.M., Мясников А.В. Широкополосный лазерный интерферометр для мониторинга деформаций Земли. Известия Академии наук, серия Физическая, 1999, т. 63, №6, с.1192-1197.

4 Vadim Milyukov, Alexandra Kozyreva, Boris Klyachko and Andrey Myasnikov. Geophisical Applications of Laser Interferomters: Long-Term Monitoring Crustal Deformations // In: Proceedings of 28th International Cosmic Ray Conference (ICRC2003), Tsukubo, Japan, 2003, pp 3157-3160.

5 B.K. Милюков, А.В.Козырева, Клячко B.C., Мясников А.В. Собственные колебания Земли.// в кн: "Современные методы геолого-геофизического мониторинга природных процессов на территории Кабардино-Балкарии". Нальчик, 2003, с. 189-199. (Russian)

6 В. К. Милюков, Б. С. Клячко, А. В. Мясников, Стриганов П.С., Янин А.Ф., Власов А.Н. Лазерный интерферометр-деформограф для мониторинга движений земной коры. // Приборы и техника эксперимента. 2005. № 6. С. 87-103.

7 Милюков В. К., Мясников А.В. Метрологические характеристики Баксанского лазерного интерферометра // Измерительная техника.

2005. № 12. С. 26-30.

8 Милюков В.К., Мясников А.В. Компенсация термоупругих и барических процессов в измерениях литосферных деформаций Баксанским лазерным интерферометром // Измерительная техника.

2006. № 2. С. 32-35.

9 В. К. Милюков, А.В. Копаев, А.В. Лагуткина, А.П. Миронов, А.В. Мясников. Наблюдения приливных деформаций земной коры в Приэльбрусье. // Физика Земли, 2007, No. 11, С. 21-29.

10 Vadim Milyukov, Alexander Kopaev, Anna Lagutkina, Alexey Mironov, Andrey Myasnikov. Observation of Luni-Solar Tides by the

Long Base Laser Interferometer // Progress of Theoretical Phys. Suppl. No. 172 (2008) pp. 165-168.

11 Vadim Milyukov, Andrey Myasnikov and Alexey Mironov. Monitoring the State of the Magmatic Structures of Elbrus Volcano Based on Observation of Lithosphere Strains // in: Nonlinear Acoustics-Fundamentals and Applications. AIP Conference Proceedings, 1022, New York, 2008, pp. 405-408.

12 Vadim K. Milyukov, Alexander V. Kopaev, Anna V. Lagutkina, Alexey P. Mironov, Andrey V. Myasnikov, Observation of tidal deformations in a tectonically active region of the Northern Caucasus // Marees Terrestres Bulletin ^Informations, no. 145 (2009), 1172911740.

13 V. Milyukov, A. Kopaev, V. Zharov, A. Mironov, A. Myasnikov, M. Kaufman, D. Duev. Monitoring crustal deformations in the Northern Caucasus using a high precision long base laser strainmeter and the GPS/GLONASS network. // J.Geodyn. 49 (2010) 216-223

14 B.K. Милюков, A.B. Мясников. Долговременные наблюдения литосферных деформаций Баксанским лазерным интерферометром в условиях подземного размещения. //Измерительная техника. 2012. № 1.С. 32-35.

15 B.K. Милюков, A.B. Мясников. Влияние термоупругих и барических процессов на измерения литосферных деформаций Баксанским лазерным интерферометром. // Измерительная техника. 2012.(в печати)

Апробация результатов

Результаты исследований, изложенные в диссертации,

докладывались на российских и международных конференциях:

• Международной конференции "Laser Optics" (С-Петербург, 1998);

• Международной конференции по космическим лучам ICRC-2003 (Япония, Тсукубо 2003);

• Генеральных ассамблеях Европейского союза наук о Земле EGU (Австрия, Вена 2005, 2006, 2008,2009, 2011).

• Сагитовских чтениях ГАИШ (Москва 2005, 2007,2009);

• Международной конференции по гравитации, астрофизике и космологи ICGA8 (Япония, Нара, 2007);

• Ломоносовских чтениях МГУ (Москва, 2007)

• Международной конференции по земным приливам ETS-2008 (Германия, Йена, 2008).

• Международном симпозиуме по нелинейной акустике 1БЫА_2008 (Швеция, Стокгольм, 2008).

• Всероссийской научной конференции «Природные процессы, геодинамика, сейсмотектоника и современный вулканизм Северного Кавказа» (Кабардино-Балкария, 2009,2010).

• Международной сейсмологической школе «Современные методы обработки и интерпретации сейсмологических данных» (Владикавказ, 2010)

• V Всероссийский Симпозиум по вулканологии и палеовулканологии. Екатеринбург, 21-25 ноября 2011г. Екатеринбург: Институт геологии и геохимии УрО РАН, 2011.

Структура диссертации

Диссертация состоит из введения, четырех глав, заключения и списка литературы. Общий объем работы составляет 156 страниц, в ней содержится 52 рисунка, 16 таблиц, список цитируемой литературы включает в себя 88 наименований.

Благодарности

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

Автор глубоко благодарен всем сотрудникам Баксанской станции и лаборатории лазерной интерферометрии: П.С.Стриганову и А.Ф.Янину. Особую благодарность автор выражает Н.А.Перелыгину за тяжелый и каждодневный (без выходных) труд по обеспечению режимных непрерывных наблюдений на Баксанском лазерном интерферометре-деформографе.

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

Во введении даны основные понятия о внутренней структуре стратовулканов центрального типа. Дается описание геологических особенностей вулкана Эльбрус и внешнее описание его кратеров. Проводится исторический анализ эруптивной деятельности с оценкой возрастов предыдущих извержений геохронологическими методами [Богатиков и др. , 1999] Имеются все основания полагать, что Эльбрус находится в непосредственной фазе очередного извержения [Уткин и др., 2009]. Вводится понятие резонансной структуры, относящейся к магматическому очагу и магматической камере Эльбруса [Собисевич и др., 2001], на основе которого будет развиваться резонансный метод оценки параметров внутренней структуры.

В первой главе приводится функциональная схема и расположение Баксанского лазерного интерферометра-деформографа (Рис.1), на базе которого были выполнены все наблюдения литосферной деформации за семилетний период (2004-2010 гг.) и дальнейшая апробация резонансного метода в условиях Эльбруского вулканического центра. Дано описание механической части интерферометра, вакуумной системы откачки и условий эксплуатации в подземных условиях штольни.

Рис. 1. Механическая и вакуумные системы Баксанского лазерного интерферометра. 1 - Не-Ме лазер; 2 - телескопическая система; 3, 4, 5 - сильфоны; 6, 7 - вакуумные камеры; 8 - невакуумированный светопровод; 9 - диафрагма. Т1, Т2 и ТЗ - датчики температуры; Р - датчик давления

Приводятся оптическая схема прибора и характеристики частотно-стабилизированного He-Ne лазера ЛГН-303. Описывается система регистрации и принцип программного обеспечения автоматизированного сбора данных.

Сбор данных ведется по 4 каналам, соответствующих разным частотным диапазонам:

1. Низкочастотный канал (геофизический)

2. 1-й высокочастотный канал (сейсмоакустический)

3. 2-й высокочастотный канал (сейсмический)

4. Канал «тремор»

Предусмотрены еще три служебных канала:

1. Температура (в 3 разнесенных точках)

2. Давление внутри системы (степень вакуума)

3. Атмосферное давление

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

Выполнена оценка пороговой чувствительности Баксанского лазерного интерферометра.

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

Выводы к главе 1.

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

2. Обеспечен режим долговременных наблюдений литосферных деформаций Баксанским лазерным интерферометром-деформографом. По техническим характеристикам (неравноплечий интерферометр Майкельсона с измерительным плечом порядка 100 м, разрешающей способностью порядка 10"13 стрейн) и способу размещения (подземный туннель нейтринной обсерватории)

инструмент относится к лучшим лазерным интерферометрам-деформографам, работающим в настоящее время в мире.

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

Запись литосферной деформации (Рис. 2) существенно осложнена двумя факторами: вариациями температуры и атмосферного давления. Отклик интерферометра как механической системы на изменения температуры происходит по двум "каналам" с разной температурной динамикой [Милюков, Мясников, 2006].

Рис. 2. Результаты наблюдений деформации литосферы (НЧ канал Еаксанского лазерного интерферометра) за период с ноября 2004 г. по декабрь 2011 г. Запись осложнена влиянием вариации годовой температуры.

Термоупругим деформациям, с одной стороны, подвергается 75-метровый светопровод с включенными в систему тремя сильфонами и вакуумные камеры (инструментальная деформация), а с другой стороны, - скальная порода и жестко связанные с ней бетонные фундаменты (деформация геофизической среды). Суточные и годовые вариации температуры в разной степени «проникают» к элементам этой сложной системы. Следует отметить, что суточные изменения температуры передаются только по первому каналу.

В первой части 2 главы изучено влияние суточных вариаций температуры и давления. Спектр температуры и давления характеризуются двумя максимумами, приходящимися на 24 и 12

часов. Причем было выяснено, что 12-ти часовой максимум температуры является следствием активной системы вентиляции штольни.

Рис.3. Амплитудные спектры вариаций деформаций (вверху) и температуры (внизу) в суточном диапазоне частот. Основные моды соответствуют суточным волнам лунно-солнечного прилива.

На Рис. 3 показан спектр деформации в суточном диапазоне (2723 Ь). В спектре выделены составляющие солнечно-лунного прилива в этом диапазоне (где СЬ, Оь Мь Рь Бь К, - приливные волны в суточном диапазоне [Мельхиор, 1968]). Оказалась существенно искаженной волна Б], которая совпадает с основной (центральная спектральная линия на Рис. 3) суточной волной температуры (24 Ь). Кроме того, суточная волна температуры имеет сезонную модуляцию (в теплое время года амплитуда суточной волны увеличивается), что привело к появлению боковых линий в спектре (Рис. 3, внизу). Эти боковые линии в точности совпадают с приливными волнами Рь Кь и таким образом вариация суточной температуры оказывает возмущения и на эти составляющие. Подобный анализ был проведен относительно полусуточной составляющей вариации температуры и, было выявлено

ее влияние на приливные волны в этом диапазоне: Т2, 82, Яг, К2 [Мельхиор, 1968].

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

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

Выводы к главе 2

1. Исследовано влияние термоупругих и барических процессов на измерения литосферных деформаций в условиях подземного размещения Лазерного интерферометра деформографа в штольне с принудительной вентиляцией.

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

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

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

Рассматриваются две модели магматической камеры: сферическая модель [Fujita ,1995] и в виде тектонически ослабленной зоны - "трещины", наполненной текучим магматическим веществом (fluid-filled crack model) [Kumagai, 2000]. Для предварительных оценок физико-механических параметров магматической камеры привлекаются геологические данные анализа образцов вулканитов последнего извержения вулкана. Для рассматриваемых давлений 0,5 -4 кбар, и значений объемной компоненты газа 30-70%, добротность мод магматического резонатора [Милюков 2006] находится в диапазоне от 150 до 500 единиц. Оценка скорости звука дает величину 180-270 м/с. Принимая во внимание результаты других методов и геологические данные, было принято решение выбрать частотный диапазон для поиска резонансных мод в интервале 0.01-0.07 Гц (150 -10 сек).

Изложена методика оценивания параметров основных спектральных мод. Реальная запись землетрясения всегда осложнена шумами различного происхождения: сейсмическими,

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

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

Ns(co) = N0 f , (1)

о + (co-L¿)

где N0 = max Ns (of) - максимальное значение спектральной

О)

плотности моды, Q - резонансная частота, 8 - полуширина спектральной линии. Таким образом, на основе метода наименьших

квадратов определяются значения параметров Мп, 5, Л , дающие максимально точное приближение лоренцевского спектрального контура (1) к наблюденной резонансной полосе.

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

Выводы к главе 3.

1. Разработана методика оценки добротности резонансной моды на основе метода максимального правдоподобия.

2. Проведено статистическое испытание метода максимального правдоподобия. Уверенно определяется период (несущая частота) сигнала на любом произвольно выбранном сейсмическом фоне. Ошибка оценки добротности для разных вариантов моделирования находилась в интервале [- 15%: +7%] от заданного значения добротности на входе модели.

В четвертой главе приводятся результаты применения резонансного метода для определения параметров внутренней структуры вулкана Эльбрус по данным записи Баксанского лазерного интерферометра-деформографа.

Раздел 4.1 посвящен геодинамической обстановке сложившейся на Северном Кавказе в целом [Хаин и др., 2010]. С позиции геотектоники, Кавказ входит в Альпийско-Гималайский сейсмический пояс, который простирается от западного Средиземноморья до Восточной Азии, расширяясь в этом направлении, сливаясь с Тихоокеанским поясом, и представляет собой мозаику литосферных микроплит. По современным представлениям между Аравийской и Евразийской плитами имеется по крайне мере две промежуточных плиты; Иранская и Закавказская. Был проведен анализ распределения гипоцентров с 2000 по 2010 годы и выявлена наклонная сейсмофокальная плоскость, что можно интерпретировать как остаточный субдукционный процесс на Северном Кавказе. Что в свою очередь указывает на определенный тип вулканизма в этом регионе.

В разделе 4.2 дается классификация всех сейсмических событий вошедших в обработку. Приведена карта распределения эпицентров

этих землетрясений. В качестве иллюстрации работы цифрового фильтра (ИКО-фильтр) и спектрального анализа, которые будут использоваться для поиска собственных колебаний магматической камеры (программа Global), построен спектр собственных колебаний Земли (СКЗ) после чилийского землетрясения 27.02.2010 г по записи литосферной деформации Баксанским интерферометром. Как видно (см. Рис. 4), все выделенные моды в спектре соответствуют своим теоретическим значениям (серые вертикальные линии [Gilbert, 1975]).

Рис 4. Спектр СКЗ вызванных Чилийским землетрясением 27.02.2010 г. Магнитуда 8.84. Вертикальные линии (oSs\0Ts -ÄjItfTW- теоретические значения сфероидальных и торсионных мод СКЗ в диапазоне 4-12 минут (модель 1066А, [Gilbert, 1975])

Раздел 4.3 посвящен выделению т.н. «семейств» резонансных мод на базе спектрального анализа всех сейсмических событий. Методика основана на том предположении, что резонансные моды, переизлученные региональными резонансными структурами, будут повторяться после каждого землетрясения, в одних и тех же узких интервалах (1 сек) периодов и образовывать устойчивые группы («семейства»). Таким образом, вводится понятие «повторяемость». Резонансные пики, которые не попали в эти группы, являются случайными.

Для этой цели, с помощью специально созданной процедуры (входящей в программный пакет GLOBAL), производилось автоматическое распознавание и «сканирование» всех возможных резонансных мод в полосе 10-150 секунд, что возникали после всех зарегистрированных землетрясений (2004 - 2011гг.).

Ill.......Illllllll.lllllllllllllllllll,.llll..llllll,l,

I... III.....................

п I "i...........i............I............I............I............ ......i..........I............I............i"--г

150 140 130 120 110 100 90 80 70 60 50 40 30 20 10

Период, сек

Рис 5. Результат статистического распределения повторяемости резонансных мод на интервале периодов 10-150 сек

В результате последней операции был получен ряд из 14 тыс. значений частот, предполагаемых в качестве резонансных. Статистическое распределение позволило весь исследуемый интервал периодов 10-150 сек разбить на пять отдельных диапазонов (рис.5), каждый из которых имеет определенный вид повторяемости резонансных мод .

I Диапазон 70- 150 сек. Обнаружены резонансные моды, но не удалось выделить устойчивые «семейства».

II Диапазон 55 - 70 сек. В этом интервален было выделено в общей сложности 10 «семейств» резонансных мод. Низкочастотная мода 1 (см. Рис. 2) была исключена из рассмотрения в силу низкой повторяемости (40%) и низкой оценки средней добротности (179 ед.).

Рис. 6. Гистограмма повторяемости резонансных мод всех рассматриваемых событий в диапазоне периодов 72 -55 сек. Цифрами обозначены номера наиболее часто повторяемых мод, в дальнейшем рассматриваемые как резонансные.

Таким образом, самая интенсивная из низкочастотных мод -мода № 2 (см. рис.6).. Период этой моды 67,9 сек, повторяемость 88%.

Были сделаны оценки периодов и добротности этих мод, и определены их среднеквадратические отклонения (СКО) (см. таблица 1). Так оценка добротности моды №2 составила 237 единиц. Именно в этом диапазоне обнаружены резонансные моды, соответствующие собственным частотам магматической камеры Эльбруса.

Таблица 1. Общий результат статистического анализа сейсмических

событий.

Региональные моды 2 3 4 5 6 7 8 9 10

Период,сек 67.9 66.5 653 64.2 63.1 62.0 60.9 59.6 56.7

СКО(мода),с 0.12 0.11 0.22 0.26 0.19 0.29 0.24 0.19 0.26

Повторяемость 0.88 0.78 0.46 0.57 0.42 0.69 0.29 0.57 0.42

Добротность 237 247 246 264 240 252 269 282 274

СКО(добротн) 47 46 50 43 33 50 45 63 48

III Диапазон 22 -55 сек. Резонансных мод практически не обнаружено

IV Диапазон 14-22 сек. Имеет явно выраженный резонансный характер. Обнаружено множество резонансных мод, но с очень низкой статистикой повторяемости. Тем не менее, высокая добротность этих мод (600-800 ед) позволяет предположить существование газового пузыря в магматическом расплаве. Учитывая , что газ находится в среде с очень низкой скоростью продольной волны (расплав), диаметр этой газовой полости можно оценить 2000-4000м.

V Диапазон 10-14 сек. Отсутствие резонансных мод.

Отдельно был выполнен анализ сейсмических событий из т.н. «ближней зоны» (расстояние до эпицентра землетрясения не превышает 1500 км). В данном случае преследовалась вполне определенная цель: выявить устойчивые семейства резонансных мод, возбужденных землетрясениями умеренной магнитуды (М\у~6), которые заведомо не способны вызвать собственные колебания Земли.

Полученные спектры нормировались по амплитуде максимальной моды в конкретном спектре и затем были просуммированы. В результате, был построен суммарный спектр, изображенный на рис 7. Это позволило выявить три устойчивые моды - 2, 3, 7 (согласно ранее введенной нумерации на рис.6), которые соответствуют только региональным резонансным структурам.

Рис. 7. Наиболее интенсивные резонансные региональные моды, возбужденные всеми рассмотренными умеренными землетрясениями в 1500-км зоне вокруг вулкана Эльбрус. Спектр получен усреднением по ансамблю всех нормированных спектров.

Таким образом, была вновь выделена мода с периодом 67.93 секунды. Рассматривая эту выделенную резонансную моду как самую низкочастотную, была вычислена оценка размера магматической камеры вулкана Эльбрус. Оценка была выполнена для сферической модели (Тчу'йа, 1995; Кита§а1, 2000].

(4.5)

где Т0 - период доминантной моды; а * - скорость звука в магматическом флюиде, а р ^ - плотность флюида; а к и р у -соответствующие параметры для окружающих твердых пород; а , ,ЪХ и С, некоторые числовые коэффициенты [Тиука ,1995]

Для значений Т0= 67.9 с (наиболее интенсивная низкочастотная мода), а{ =200 м/с, р( =2000 кг/м3, =6000 м/с и рх =2700 кг/м3,

получаем Ь ~ 9 км,

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

т Тоа/ Ь » х

2 л

а, + Ьх х tgh

Рг аг

зондирование, гравиметрическая съемка), и подтверждает результат, полученный ранее в работе [Милюков, 2006].

В Разделе 4.3 была изучена динамика изменения добротности магматической камеры. На рис. 8 приведена оценки динамики добротности за период наблюдения с 2004 по 2011 гг. С помощью метода наименьших квадратов (МНК) было установлено, что добротность камеры как магматического резонатора уменьшается, в среднем, на 6 единиц в год.

400

350-

300-

150100-

Рис 8. Динамика добротности региональных мод в период с ноября 2004г. по декабрь 2010г. Значение добротности для каждого землетрясения определялось как среднее значение по всем зарегистрированным модам. Прямая линия динамика добротности магматического очага как акустического резонатора (6+2 единиц/год.).

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

Выводы к главе 4 1. Выделено 10 устойчивых групп резонансных мод (семейств). Выполнена оценка добротности всех резонансных мод входящий в эти группы. На основании этого стало возможным отнести 9 из

2005

2007 2008 Время, годы

2009

2011

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

2. Выполнена оценка размера магматической камеры вулкана Эльбрус с помощью резонансного метода - 9км. Глубина залегания камеры 3-7 км.

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

Заключение

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

Список цитируемой литературы

1. Богатиков O.A., Гурбанов А.Г., Мелекесцев И.В. и др. Проблема активизации вулкана Эльбрус (Северный Кавказ) и возможные ее катастрофические последствия // Глобальные изменения природной среды и климата: Избранные тр. Новосибирск: СО РАН, НИЦ ОИГГМ, 1999. С.145-156.

2. Богатиков O.A., Нечаев Ю.В., Собисевич A.JI. Уточнение структурных особенностей материнского магматического очага и вулканической камерой Эльбруса // Современные математические и геологические модели в задачах прикладной геофизики (Геофизика и математика XXI). М.: ИОФЗ РАН, 2001. С. 249-260.

3. Копаев А. В., А. Г. Гурбанов. Гравиметрические исследования в Геналдонском ущелье: первые результаты. Вестник Владикавказского научного центра, Том 4, №3, 2004 г.

4. Мельхиор П. Земные приливы. М.: "Мир". 1968. 483 с.

5. Милюков В.К., Клячко Б.С., Мясников A.B., Стриганов П.С., Янин А.Ф., Власов А.Н. Лазерный интерферометр-деформограф для мониторинга движений земной коры. // Приборы и техника эксперимента. 2005. № 6. С. 1-17.

6. Милюков В. К., Мясников A.B. Метрологические характеристики Баксанского лазерного интерферометра // Измерительная техника. 2005. № 12. С. 26-30.

7. Милюков В.К. Мониторинг состояния магматических структур вулкана Эльбрус по наблюдениям литосферных деформаций //Вулканология и сейсмология. 2006. № 1.С. 1-13

8. Нечаев Ю.В. Космические технологии в задачах изучения локальных неоднородностей земной коры // Избранные труды ученых ОИФЗ РАН «Геофизика на рубеже веков». М.: ОИФЗ РАН, 1999. С. 276-290.

9. Рогожин Е.А., Собисевич Л.Е., Нечаев Ю.В., Собисевич А.Л., Богатиков O.A., Гурбанов А.Г., Коваленко В.И., Газеев В.М. Поляк Б.Г., Покровский Б.Г., Лаврушин В.Ю., Куликов В.И., Мелекесцев И.В., Кащук Д.Г., Милюков В.К., КопаевА.В. Геодинамика, сейсмотектоника и вулканизм Северного Кавказа. Под ред. акад. Н.П. Лаверова. М.: Изд-во ИФЗ РАН. 2001. 333 с.

10. Собисевич Л.Е., Милюков В.К., Собисевич А.Л. Механико-математический мониторинг магматических структур вулкана Эльбрус // Современные математические и геологические модели в задачах прикладной геофизики. М.: Изд.ОИФЗ РАН, 2001г. С 223249.

11. Собисевич А.Л., Руденко О.В., Милюков В.К., Нечаев Ю.В. Мониторинг наведенных геофизических процессов в гетерогенных структурах геологической среды вулканов центрального типа // Катастрофические процессы и их влияние на природную среду. Том 1. Вулканизм / Под ред. Н.П. Лавёрова. М.: ИЗФ РАН, 2002г. С 365-397.

12. Уткин И.С, Федотов С.А., Уткина Л.И. Оценка тепла, накопленного магматическим очагом вулкана Эльбрус во вмещающих его породах, и возможности его извлечения. // Вулканология и сейсмология, 2009г. № 5, С. 1-21.

13. Шевченко В.И., Гусева Т.В., Лукк A.A., Мишин A.B., Прилепин М.Т., Рейлинджер Р.Э., Хамбургер М.У., Шемпелев А.Г., Юнга С.Л. Современная геодинамика Кавказа (по результатам GPS измерений и сейсмологическим данным) // Физика земли. 1999. № 9. С. 3-18.

14. Хаин В.Е., Ломизе М.Г. Геотектоника с основами геодинамики. М: КДУ, 2010. 560

15. Чернышев И.В., Лебедев В.А., Бубнов С.Н. и др. Этапы магматической активности Эльбрусского вулканического центра (Большой Кавказ): изотопно-геохронологические данные // Доклады АН. 2001. Т. 380. № 3. С. 384-389.

16. Fujita Е., Ida Y., Oikawa J. Eigen oscillation of a fluid sphere and source mechanism of harmonic volcanic tremor // J. Volcanol. Geotherm. Res. 1995. V. 69. P. 365-378.

17. Gilbert F., Dziewonski A.M. An Application of Normal Mode Theory to the Retrieval of Structural Parameters and Source Mechanisms from Seismic Spectra // Philosophical Transactions of the Royal Society of London. A. Math, and Phys. Sci. 1975. V. 278. N 1280. P. 187.

18. Kumagai H., Chouet B.A. Acoustic properties of a crack containing magmatic or hydrothermal fluids // J. Geophys. Res. 2000. V. 105. № В11. P. 25493-25512.

Мясников Андрей Владимирович

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

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

Подписано в печать 12.04. 2012 г. Заказ №159/11А Формат 60 х 90/16. Усл.печ. л. 1.Тираж 120 экз.

Содержание диссертации, кандидата физико-математических наук, Мясников, Андрей Владимирович

ВВЕДЕНИЕ.

Актуальность темы.

Цель работы.

Научная новизна.

Защищаемые положения.

Личный вклад автора.

Апробация результатов.

Практическая значимость.

Структура диссертации.

ГЛАВА 1. БАКСАНСКИЙ ЛАЗЕРНЫЙ ИНТЕРФЕРОМЕТР-ДЕФОРМОГРАФ.

1.1. Лазерные интерферометры-деформографы для геофизических наблюдений.

1.2. Основные узлы Баксанского лазерного интерферометра.

1.2.1. Размещение.

1.2.2. Оптическая схема интерферометра.

1.2.3. Механическая и вакуумная системы.

1.2.4. Источник излучения.

1.2.5. Система регистрации.

1.3. Автоматизированная система сбора данных.

1.3.1. Работа программного обеспечения.

1.3.2. Выходные данные системы.

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

1.5. Пакет программ предварительной обработки и анализа данных.

1.6. Контроль качества регистрируемой информации.

1.7. Предварительная обработка данных.

1.8. Результаты главы 1.

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

2.1. Долговременные наблюдения деформаций и метеорологический фактор.

2.2. Влияние метеорологических факторов в приливном диапазоне.

2.2.1 Суточные вариации.

2.2.2. Полусуточные вариации.

2.3. Влияние метеорологических факторов в годовом диапазоне.

2.4. Адаптивный компенсатор метеорологических возмущений.

2.5. Обработка наблюдательных данных с помощью адаптивного компенсатора.

2.6. Результаты главы 2.

ГЛАВА 3. ОСНОВЫ РЕЗОНАНСНОГО МЕТОДА КОНТРОЛЯ ДИНАМИЧЕСКИХ ПРОЦЕССОВ И СОСТОЯНИЯ МАГМАТИЧЕСКИХ СТРУКТУР ВУЛКАНА ЭЛЬБРУС.

3.1. Основные методы изучения внутренней структуры вулкана Эльбрус.

3.1.1. Метод теплового дистанционного зондирования.

3.1.2. Линеаментный метод.

3.1.3. Метод магнитотеллурического зондирования.

3.1.4. Гравиметрический метод.

3.2. Принципы резонансного метода выявления и контроля состояния магматических структур «спящих» вулканов.

3.3 Предварительные оценки геометрического размера и физико-механических параметров магматической камеры Эльбруса.

3.4. Методика спектрального оценивания основных параметров резонансных мод.

3.5. Результаты главы 3.

ГЛАВА 4. РЕЗОНАНСНЫЕ ХАРАКТЕРИСТИКИ РЕГИОНАЛЬНЫХ НЕОДНОРОДНЫХ СТРУКТУР ПРИЭЛЬБРУСЬЯ.

4.1. Сейсмотектоническая характеристика Северного Кавказа.

4.2. Регистрация вторичных волн от магматической камеры Эльбруса Баксанским лазерным интерферометром.

4.2.1. Характеристика наблюдательных данных 2004-2011гг.

4.2.2. Конструирование цифрового фильтра и предварительная подготовка данных

4.2.3. Оценка резонансных параметров.

4.2.4. Оценка резонансных параметров по землетрясениям ближней зоны».

4.2.5. Оценка размера магматической камеры вулкана Эльбрус.

4.3. Динамика добротности магматической камеры вулкана Эльбрус.

4.4. Результаты главы 4.

Введение Диссертация по наукам о земле, на тему "Мониторинг состояния магматических структур вулкана Эльбрус по наблюдениям литосферных деформаций Баксанским лазерным интерферометром"

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

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

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

Вершины большинства стратовулканов обычно находятся на высоте от 2000 до 2500 м над их основанием; вершины с высотой более 3000 м над основанием встречаются редко. Объем большинства слоистых конусных вулканов над зонами субдукции составляет примерно 200 км3, что справедливо и для вулканической постройки Эльбруса (152 км3) [Уткин, 2009].

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

На рисунке В.1 приведена принципиальная схема развитой системы магматических каналов и очагов для вулканов центрального типа [Федотов и др., 1991].

ЕЗ 1 £3 з

ЕЗ 2 ЁЕЭ 4

Рис. В.1 Схема развитой системы магматических каналов и очагов для вулканов центрального типа. Обозначения: А - астеносфера, Ь - граница литосферы, М -граница земной коры, Н - подкоровый, или нижнекоровый магматический очаг, С -коровый промежуточный очаг, возникающий при плавлении гранитных или осадочных пород основными магмами, В - верхний, периферический очаг, 1 -«базальтовый» слой, 2 - «гранитный» слой, 3 - возможная область плавления в «гранитном» слое вокруг магматического канала, 4 - осадочный слой, I -магматическая астеносферная колонна, II - часть питающего канала, находящаяся в «базальтовом» слое; III - то же, в «гранитном слое», IV - между «периферическим» очагом и кратером [Федотов и др., 1991].

Именно эти структуры, являющиеся своеобразным «сердцем» вулканов и определяют его «жизненные» циклы. Характеризуя последние, укажем, что экспериментальные наблюдения последних лет, проведенные в районе Эльбрусского вулканического центра, дают основания отнести магматический очаг, магматическую камеру, коровые и другие каналы к геологическим структурам резонансного типа, определяющим механизмы подготовки и развития катастрофического извержения [Собисевич, Шумейко и др., 2000; Собисевич, Милюков и др., 2001].

Вулкан Эльбрус приурочен к месту пересечения продольной Сылтранской магмоконтролирующей разломной зоны с поперечным Эльбрусским разломом и «насажен» на древние кристалические породы, образующие горный блок. Диаметр основания вулкана около 14-15 км с относительной высотой 3 км. Западная и восточная вершины, а также общий конус вулкана приурочены к огромной (площадь 230 км2) кальдере обрушения. Общий объем изверженных продуктов за весь обозримый период существования вулкана оценивается [Уткин, 2009] до 190 км3. Как показали расчеты, выполненные по алгоритму С.А. Федотова, объем современной вулканической постройки составляет 152 км3.

Итак, две вершины Эльбруса относятся к разным по типу, размеру и возрасту вулканическим постройкам. Западная вершина - это верхняя часть более старого и крупного стратовулкана, а восточная - является побочным, относительно небольшим, но также многоактным вулканом.

Западная вершина имеет более сложное строение и худшую сохранность, чем восточная. Так, в пределах верхней части постройки западной вершины выделяются два кратера диаметром 1,0 и 0,5 км, а также центральный купол. Весь западный сектор западной вершины разрушен сильными взрывами и гигантскими обвалами, о чем свидетельствуют как минимум три обвальных цирка и открытый на ЗЮЗ обвально-взрывной кратер размером 2,3x1,8 км.

Восточная вершина устроена гораздо проще, там имеются лишь очень свежий по облику вершинный кратер овальной формы размером 0,35x0,3 км, а в ВЮВ секторе -открытый на ВЮВ удлиненный кратер размером 0,6x0,35 км. В ССВ секторе, примерно в 1км от кромки вершинного кратера находится еще один эруптивный центр с небольшим лавовым куполом, из которого излился лавовый поток.

Считается, что вершины Эльбруса - два самостоятельных вулкана, выросших на более древнем вулканическом основании [Богатиков, 1998]. Восточный конус (5621 метров) очень «молод», сохранил правильную вулканическую форму с четко выраженной чашей-кратером. Западный конус (5642 метров) древнее и сильно изменен, почти треть его верхней части разрушена вертикальным разломом.

Начало изученного периода вулканической деятельности Эльбруса относят к концу позднего плиоцена. В истории формирования вулкана выделены докальдерный цикл, кальдерный цикл, начало которого ознаменовалось мощным эксплозивным извержением, приведшим к образованию кальдеры, и посткальдерный цикл, в результате которого сформировался современный стратовулкан Эльбрус [Богатиков, 1998!].

Докальдерные вулканические образования пока достоверно не установлены. о і г і о с +

ШЁ ауюаюмдойв пее^бАтЗма/мп) ^

• шитрм «0

• маошОо < Э4) « іріуушш'мкмл іротар *ТТ* Омжої мшмі Дмшуиаа» 1037>

• ЮЮІИ) гоюси

• івгаив /••уіммвмчя б-предлотгкмкмі

• и-іуп ипшги аидп« (0]7)

• й-гимбсжть ¥ тусЫ кДЛС^ООО&ШГКиеГО

• гаоші пктша Этшою —- І »опрлок»тя ожирію гкро

ХТО<КО Ш|Мв . апуиутиви« (іа:шь

•«утре «а«і<разіміив гав

• поспаіі4Сріи) туоинсш іссон плоди V лостоскі»погг€«:о калиіор« <озоаж»и<* иопжг

Рис. В.З. Схематическая карта Эльбрусской кальдеры и вулкана Эльбрус согласно [Богатиков, Гурбанов, 2003]

Кальдера вулкана образовалась по уточненным данным около 800 тыс. лет тому назад. Кальдерный вулканизм связан с образованием крупной (17 х 14 км по бровке ограничивающего уступа) Эльбрусской кальдеры обрушения. Периферический магматический очаг вулкана до образования кальдеры, судя по ее размерам, мог иметь диаметр около 16-18 км. На период, предшествовавший образованию кальдеры, и непосредственно после ее образования, в изверженных продуктах преобладали игнимбриты и ассоциирующие с ним туфы, что говорит об эксплозивном характере извержений того периода.

W, км3/год 0.003 Ч

0.002

0.001

W(t) а ,км ^250

200

F-150

-100

200 175 150 125 100 75 50 25 t, тыс. лет

I I I I I I I

225 200 Т | | м |

175 150 | | | | | | | | 125 100

Г 75

I I I I I I I

50 25 t, тыс. лет

Рис. В.4. Периоды вулканической активности вулкана Эльбрус за последние 225 тысяч лет [Уткин и др., 2009]. а. Кусочно-постоянные зависимости расхода магмы для вулкана Эльбрус за последние 225 тыс. лет (по геологическим данным) W(t) м3/год. Зависимость объема изверженных продуктов от времени: J W{t)dt км3. б. Зависимость размера магматической камеры вулкана Эльбрус от времени с начала современного зарождения (225 тыс. лет тому назад) до настоящего времени.

Образование кальдеры не могло не повлиять на значительное уменьшение объема магматического очага вулкана. Именно это могло привести к периоду длительного затишья в активности вулкана -550 тыс. лет, наступившего после образования кальдеры.

Новый период вулканической активности начался примерно 225 тыс. лет тому назад. По материалам работы [Уткин] была построена зависимость расхода магмы и объема изверженных продуктов за этот период, которая наилучшим образом отражает периоды активности вулкана, (Рис. 1а) а также эволюция радиуса магматического очага в процессе кондуктивного теплообмена (см. рис. 16).

В период от 225 до 180 тыс. лет до н.в. деятельность вулкана была наиболее интенсивной за весь посткальдерный период. Формируется западная вершина вулкана Эльбрус. Именно в этот период вулканической деятельности, средний расход извергнутых продуктов и их количество, примерно в 4 раза превышает те же показатели за период от 12 тыс. лет до н.в. и до нашего времени. Именно такая вулканическая активность требуется для того, чтобы в пассивный последующий период, с 180 до 130 тыс. лет до н.в., магматический очаг не застыл. По данным геологических исследований магматический очаг в этот период времени не застывал. Этот перерыв вулканической деятельности, 50 тыс. лет, самый протяженный за рассматриваемый период. В период 130160 тыс. лет до н.в. вулкан снова активен. Затем опять следует перерыв 25 тыс. лет. И, наконец, последние 35 тыс. лет вулкан Эльбрус находится в очередной активной фазе.

Определение возраста с помощью разных геохронологических методов показало, что отдельные извержения этого периода происходили 35 + 5, 28 + 3, 23 + 2 и около 21 тысячи лет тому назад. За последние 10-12 тыс. лет в голоцене интенсивность деятельности вулкана Эльбрус значительно возрастает, паузы между периодами активности сокращаются, средний расход изверженных продуктов также возрастает. Образована восточная вершина вулкана. Последний, голоценовый период вулканической активности, превосходит по интенсивности предыдущий, длившийся 70 тыс. лет с 130 по 60 тыс. лет до н.в. За период с 8150 ± 40 по 990 ± 60 лет по данным 14 С радиоуглеродного датирования вулканическая активность возобновлялась 8 раз с интервалом 400-1700 лет. Сильное последнее извержение вулкана происходило в I-II вв. нашей эры и менее интенсивное - 900 лет назад. Таким образом, в современную эпоху весьма вероятней новый период активности вулкана Эльбрус.

Актуальность темы

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

Следует помнить, что территория Северного Кавказа это территория с интенсивным развитием хозяйственной деятельности, высокой плотностью населения. Кроме известных факторов катастрофических последствий вероятного извержения, таких как лавовые и пирокластические потоки, угрожающие непосредственно прилегающим территориям, следует сказать, что Эльбрус обладает самой крупной ледниковой системой Кавказа, состоящей из 25 ледников. Суммарная физическая поверхность всех ледников с учетом снежно-ледового покрова вершин Эльбруса составляет 139 км2, а суммарный объем льда оценивается в 6 км . Наличие обширного снежно-ледового покрова делает вулкан Эльбрус еще более опасным, так как в случае будущих любой силы и типа извержений к собственно вулканической опасности обязательно добавится катастрофическая опасность от образования лахаров и крупномасштабных наводнений, которые могут предваряться валом воды до нескольких десятков метров. Таким образом, опасности подвергается не только населенные пункты непосредственно расположенные в Приэльбрусье, а практически весь Северный Кавказ в целом, особенно территории в долинах рек Баксан, Кубань и Малка.

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

Цель работы

Разработка методики выявления собственных колебаний резонансных структур, возбужденных воздействием сейсмической волны. Определение геометрических и физико-механических параметров магматической камеры Эльбруса на основе разработанных моделей магматических структур вулканов.

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

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

Научная новизна

2. Проведен спектральный анализ литосферных деформаций, возбужденных удаленными землетрясениями с магнитудами Мш 4.5-9.1, за более чем 7-летний период наблюдений. В результате этого анализа в спектрах литосферных деформаций с помощью резонансного метода выявлены региональные резонансные моды, относящиеся к магматическим структурам вулкана Эльбрус. Получены оценки характерного размера, глубины залегания, состава магматических флюидов близповерхностной магматической камеры вулкана Эльбрус.

Защищаемые положения

Личный вклад автора

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

В течение более чем 20 лет автор принимал активное участие в работе коллектива по созданию Баксанского лазерного интерферометрического комплекса, обеспечению его работы в режиме долговременных наблюдений.

Апробация результатов

Результаты исследований, изложенные в диссертации, докладывались на российских и международных конференциях:

• Международной конференции "Laser Optics" (С-Петербург, 1998);

• Международной конференции по космическим лучам ICRC-2003 (Япония, Тсукубо 2003);

• Генеральных ассамблеях Европейского союза наук о Земле EGU (Австрия, Вена 2005, 2006, 2008, 2009, 2011).

• Сагитовских чтениях ГАИШ (Москва 2005, 2007,2009);

• Международной конференции по гравитации, астрофизике и космологи ICGA8 (Япония, Нара, 2007);

• Ломоносовских чтениях МГУ (Москва, 2007)

• Международной конференции по земным приливам ETS-2008 (Германия, Йена, 2008).

• Международном симпозиуме по нелинейной акустике ISNA2008 (Швеция, Стокгольм, 2008).

Практическая значимость

Структура диссертации

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Мясников, Андрей Владимирович

Заключение

Диссертационная работа является результатом 20-летней работы автора в лаборатории лазерной интерферометрии ГАИШ МГУ и непосредственно на Баксанской станции в Приэльбрусье, где установлен лазерный интерферометр-деформограф. Этот период времени охватывает непосредственно момент создания прибора до настоящего времени.

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

Основные итоги выполненной работы:

• Детально изучен механизм влияния температуры и атмосферного давления на различные узлы механической части лазерного интерферометра-деформографа. Дано решения этой задачи в разных частотных диапазонах, начиная с полусуточного приливного диапазона до годового диапазона.

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

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

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

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

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

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

Что касается резонансного метода, то его развитие видится в направлении изучения влияния лунно-солнечного прилива на акустические свойства магматической камеры как резонатора. Согласно теории [Мельхиор, 1968], вследствие радиальных приливных смещений происходит последовательное объемное увеличение и уменьшение полостей в земной коре. Отсюда следует, что происходит последовательное изменение давления внутри камеры, что, в свою очередь, приводит к обратному изменению добротности камеры как механического резонатора и, вероятно, к смещению собственных частот резонатора.

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

Так автор видит дальнейшее совершенствование инструментальной техники, развитие банка данных и резонансного метода в частности.

Библиография Диссертация по наукам о земле, кандидата физико-математических наук, Мясников, Андрей Владимирович, Москва

1. Авдулов М.В. О геологической природе гравитационной аномалии Эльбруса // Известия АН СССР. Сер. Геол. 1962. №9. С. 67-74.

2. Авдулов М.В. Строение земной коры по данным гравиметрии на Центральном Кавказе // Советская геология. 1963. № 9. С. 73-89

3. Авдулов М. В. Методика и результаты количественной интерпретации Эльбрусского гравитационного минимума, Физика Земли №5, 1994 г.

4. Авдулов М. В., Короновский Н. В., О геологической природе Эльбрусского гравитационного минимума, Вестн. Моск. ун-та. серия 4. геология. №3, 1993 г.

5. Авсюк Ю.Н. Приливные силы и природные процессы. М.: ОИФЗ РАН, 1996. 188с.

6. Алешин В.А., Дубов H.H., Яковлев А.П. Геофизический лазерный деформограф штольневого типа // Известия РАН. Сер. Физика Земли. 1993. № 4. С. 62.

7. Арбузкин В.Н., Компаниец М.А., Швец А.И., Греков И.Т., Литовко Г.В. и др. Отчет о комплексных геолого-геофизических исследованиях по Приэльбрусскому профилю. ФГУП «Кавказгеолсъемка». Ессентуки, 2002. 120с.

8. Атлас ледников Эльбруса. Часть 1. Фотоснимки ледников. М.: Изд-во МГУ, 1965. 20с.

9. Багаев С.Н., Орлов В.А., Фомин Ю.Н., Чеботаев В.П. Лазерные деформографы для прецизионных геофизических измерений // Известия АН СССР. Сер. Физика Земли. 1992. № 1. С. 85.

10. Бат Маркус. Спектральный анализ в геофизике. М.: "Недра". 1980. 535 с.

11. Белоусов В.В. Большой Кавказ как тектоническая лаборатория // Проблемы геодинамики Кавказа. М.: Наука, 1982. С. 9-13.

12. Белоусов В.В. Тектоносфера Земли: взаимодействие верхней мантии и коры // МГК АН СССР. М., 1991. 71с.

13. Бердичевский М. Н. Электрическая разведка методом магнито-теллурического профилирования, М., 1968.

14. Богатиков O.A., Мелекесцев И.В., Гурбанов А.Г. и др.//Катастрофические палеолахары вулкана Эльбрус (Северный Кавказ) ДАН. 1998. Т. 363. № 4. С. 515-517.

15. Богатиков O.A., Мелекесцев И.В., Гурбанов А.Г. и др. Радиоуглеродное датирование голоценовых извержений вулкана Эльбрус (Северный Кавказ) ДАН. 1998. Т. 363. № 2. С. 219-221.

16. Богатиков O.A., Мелекесцев И.В., Гурбанов А.Г. и др. Эльбрусская кальдера (Северный Кавказ) ДАН. 1998. Т. 363. № 4. С. 515-517

17. Богатиков О.В. , Мелекесцев И.В., Гурбанов А.Г. и др. Катастрофическая плейстоценовая и голоценовая активность вулканического центра Эльбруса (Северный Кавказ). \\ Докл. РАН 1998г. Т. 362. № 4. С. 518-521.

18. Богатиков O.A., Нечаев Ю.В., Собисевич A.JI. Использование космических технологий для мониторинга геологических структур вулкана Эльбрус. ДАН. 2002. Т. 387. № 3. С. 1-6.

19. Богатиков О. А., Гурбанов А. Г. Комплексные исследования Эльбруского и Казбекского вулканических центров: мониторинг и прогноз. Вестник Владикавказского научного центра, Том 3, №2, 2003 г.

20. Богатиков O.A., Залиханов М.Ч., Карамурзов Б.С. и др. Природные процессы на территории Кабардино-Балкарии // Под ред. акад. Н.П. Лаверова. М.: Изд-во ИГЕМ РАН. 2004. 438 с.

21. Бондур В.Г., Зверев А.Т. Метод прогнозирования землетрясений по результатам линеаментного анализа космических изображений. Изв. ВУЗов. Геодезия и аэрофотосъемка. 2005,№1, с.76-83

22. Буклерский A.B., Карт A.M., Клячко Б.С., Милюков В.К., Мясников A.B. и др. Баксанский лазерный интерферометр // Измерительная техника. 1995. № 10. С. 5.

23. Васильева О.В., Карабутов A.A., Лапшин Е.А., Руденко О.В. Взаимодействие одномерных волн в средах без дисперсии. М.: Изд. МГУ, 1983. 82с.

24. Власов А.Н., Перебякин В.А., Поляков С.Ю., Привалов В.Е. Долговременная стабильность и воспроизводимость частоты He-Ne лазера с внутренними зеркалами // Квантовая электроника. 1986. Т. 13. № 2. С. 320.

25. Гаджиев Т.Г., Нечаев Ю.В., Потапова Е.И., Саттарова В.М. Карта глубинных структур Кавказа по космическим данным, масштаб 1 : 1 000 000. Баку: Изд-во ПО Азербайджангеодезия ГУГК СССР, 1989. 4 листа.

26. Газеев В.М., Носова A.A., Сазонова Л.В.*, Гурбанов А.Г., Докучаев А.Я. Петрогенетическая интерпретация ассоциаций минералов-вкрапленников плейстоценовых- голоценовых вулканитов Эльбруса (Северный Кавказ) // Вулканология и сейсмология (в печати)

27. Гурбанов А. Г., Богатиков О. А., Карамурзов Б. С., Цуканова Л. Е., Лексин А. Б., Газеев В. М., Мохов А. В., Горностаева Т. А., Жариков А. В., Шмонов В. М.,

28. Докучаев А.Я., Горбачева С. А., Шевченко А. В. Необычные виды дегазации из расплавов периферических магматических камер «спящего» вулкана Эльбрус(Россия). Геохимические и минералогические особенности. // Вулканология и сейсмология. 2011, № 4, с. 3-20

29. Звягинцев Л.И., Гурбанов А.Г. Петрофизические типы арисцийских гранитоидов Главного Кавказского хребта и их рудоносность // Изв. АН СССР. Сер. геол. 1986. № 9. С. 13-26.

30. Кадик A.A., Луканин O.A., Лапин И.В. Физико-химические условия эволюции базальтовых магм в приповерхностных очагах,- М.: Наука, 1989.- 346 с.

31. Катастрофические процессы и их влияние на природную среду. / Под редакцией академика Лаверова Н.П. Том 2. Сейсмичность. М.: Региональная общественная организация,ученных по проблемам прикладной геофизики, 2002. 506с.

32. Козырева A.B., Милюков В.К. Оценка резонансных характеристик магматического очага вулкана Эльбрус по деформационным наблюдениям // Вестн. Моск. ун-та. Физика. Астрономия. 2003. № 3. С. 62-67.

33. Копаев А. В., А. Г. Гурбанов. Гравиметрические исследования в Геналдонском ущелье: первые результаты. Вестник Владикавказского научного центра, Том 4, №3, 2004 г.

34. Короновский Н.В. Геологическое строение и история развития вулкана Эльбрус // Оледенение Эльбруса,- М.: МГУ, 1968,- С. 15-72.

35. Ле Пишон К., Франшто Ж., Боннин Ж. Тектоника плит // Изд. «Мир» М. 1977г.

36. Лучицкий И.В. Основы палеовулканологии. М.: Наука, 1971. Т.1. 480с.

37. Маков Ю.Н., Руденко О.В., Собисевич Л.Е., Собисевич А.Л. О резонансных явлениях в геофизической среде // Развитие методов и средств экспериментальной геофизики. М.: ИОФЗ РАН, 1996. В. 2. С. 194-200.

38. Малахов А. Гравитационные аномалии вулканических районов. // Измерительная техника. 1990. № 11. С. 9.

39. Мельхиор П. Земные приливы. М.: "Мир". 1968. 483 с.

40. Милюков В.К., Руденко В.Н., Клячко Б.С., Карт A.M., Мясников A.B. Широкополосный лазерный интерферометр для мониторинга деформаций Земли // Известия Академии наук. Сер. Физическая. 1999. Т. 63. № 6. С. 11921197.

41. Милюков В. К., Мясников A.B. Метрологические характеристики Баксанского лазерного интерферометра // Измерительная техника. 2005. № 12. С. 26-30.

42. Милюков В.К. Мониторинг состояния магматических структур вулкана Эльбрус по наблюдениям литосферных деформаций //Вулканология и сейсмология. 2006. № 1.С. 1-13

43. Нечаев Ю.В. Космические технологии в задачах изучения локальных неоднородностей земной коры // Избранные труды ученых ОИФЗ РАН «Геофизика на рубеже веков». М.: ОИФЗ РАН, 1999. С. 276-290.

44. Расцветаев JI.M. О геологической природе линеаментов, выявленных на космических изображениях Кавказа // Известия вузов. Геол. и разведка, 1974, №12, с. 58-66.

45. Розен О.М., Федоровский B.C. Коллизионные гранитоиды и расслоение земной коры,- М.: Научный мир, 2001.- 188 с. (Тр. ГИН РАН, Вып. 545).

46. Саваренский Е.Ф., Нерсесов И.Л., Кармалеева P.M., Латыниа Л.А.// Физика Земли. 1966. № 5. С. 33.

47. Собисевич А.Л. Избранные задачи математической геофизики и вулканологии. М.: ИФЗ РАН, 2010.-464 с.

48. Собисевич Л.Е.,Милюков„В.К., СобисевинА.Л.Механико-математинеский. мониторинг магматических структур вулкана Эльбрус // Современные математические и геологические модели в задачах прикладной геофизики. М.: Изд.ОИФЗ РАН, 2001г. С 223-249.

49. Теплофизические свойства технически важных газов при высоких температурах и давлениях. Справочник / Под ред. Кузнецова В.М., Зубарева В.Н., Козлова А.Д. и др. М.: Энергоатомиздат, 1989. 232 с.

50. Толстых М.Л., Наумов В.Б., Гурбанов А.Г. и др. Состав магматических расплавов вулкана Эльбрус и Казбек (Кавказ) по данным изучения включений в минералах // Геохимия. 2001. № 4. С. 441-448.

51. Трифонов В.Г. Неотектоника Евразии. М.: Научный мир, 1999. 252с.

52. Уткин И.С, Федотов С.А., Уткина Л.И. Оценка тепла, накопленного магматическим очагом вулкана Эльбрус во вмещающих его породах, и возможности его извлечения. // Вулканология и сейсмология, 2009г. № 5, С. 121.

53. Федотов С.А. О входных температурах магм, образовании, размерах и эволюции магматических очагов вулканов // Вулканология и сейсмология. 1980. №4.С. 3-29.

54. Федотов С.А. Расчет питающих каналов и магматических очагов вулканов, имеющих устойчивые размеры и температуру // Вулканология и сейсмология. 1988. № 3. С. 3-7

55. Федотов С.А., Горицкий Ю.А. Тепловой расчет цилиндрических питающих каналов и расхода магм для вулканов центрального типа. 4.1 // Вулканология и сейсмология. 1979. №6. С.78-93.

56. Федотов С.А., Гусев А.А., Чернышова Г.В., Шумилина JI.C. Сейсмофокальная зона Камчатки: (Геометрия, размещение очагов землетрясений и связь с вулканизмом) //Вулканология и сейсмология. 1985. №4. С. 91-107.

57. Флоринский И.В. Визуализация линеаментов и структур центрального типа: количественные топографические подходы. Пущино, ПНЦ РАН, 1992, 48с.

58. Шевченко В.И., Гусева Т.В., Лукк А.А., Мишин А.В., Прилепин М.Т., Рейлинджер Р.Э., Хамбургер М.У., Шемпелев А.Г., Юнга С.Л. Современная геодинамика Кавказа (по результатам GPS измерений и сейсмологическим данным) // Физика земли. 1999. № 9. С. 3-18.

59. Шенгелая Г. Ш. Гравитационная модель земной коры Кавказа. М., «Наука», 1984 г.

60. Хаин В.Е. Тектоника линеаментов и неомобилизм / Geologica Balcanica, 8.3, Sofia, Sept.,1978, р.З 8.

61. Хаин В.Е., Лобковский Л.И. Условия проявления остаточной мантийной сейсмичности альпийского пояса Евразии // Геотектоника. 1994. - № 3. С. 1220.

62. Хаин В.Е. Тектоника континентов и океанов (год 2000). М.: Научный мир, 2001.-606 с.

63. Хаин В.Е., Ломизе М.Г. Геотектоника с основами геодинамики. М: КДУ, 2010. 560

64. Большой Кавказ в альпийскую эпоху. // Под ред. Ю.Г. Леонова. М.: ГЕОС, 2007. 368

65. Хитаров В.И. К оценке активности вулкана Эльбрус // Докл. АН СССР, 1984, т.275, №4, с. 952-984.

66. Якушева А.Ф., Хаин В.Е., Общая геология. М.: Мир, 1988.

67. Agnew D. С. Strainmeters and Tiltmeters // Rev. Geophys. 1986. V. 24 . № 3. P. 579624

68. Chouet B.A. A seismic model for the source of long-period events and harmonic tremor. // In: Volcanic Seismology. Ed. by Gasparini P., Scarpa R. and Aki K. Springer Verlag, New York. 1992. P. 133-156.

69. Crescentini L., Amoruso A., Fiocco G., Visconti G. Installation of a high-sensitive laser strainmeter in a tunnel in central Italy // Rev. Sci. Instrum. 1997. V. 68. № 8. P. 3206-3210.

70. Fujita E., Ida Y., Oikawa J. Eigen oscillation of a fluid sphere and source mechanism of harmonic volcanic tremor // J. Volcanol. Geotherm. Res. 1995. V. 69. P. 365-378.

71. Kopaev A.V., Milyukov V.K. Geodynamical Investigation Program in Baksan Canyon Area // In: Geodesy and Physics of the Earth. Proc. VII Internat. Sympos. Potsdam. 1992. P. 147.

72. Kopaev A., Milyukov V. Environmental effects in tide strain observations near the Mt. Elbrus, Central Caucasus // Marees Terrestrees. Bull. d'Inform. 2002. N 137. P. 10909-10916

73. Kumagai H., Chouet B.A. Acoustic properties of a crack containing magmatic or hydrothermal fluids // J. Geophys. Res. 2000. V. 105. № B11. P. 25493-25512.

74. Milyukov V., Azarova V., Goliaev Yu., Solovieva T. Sensing of Earthquake Precursors with Laser Interferometer // Proceedings of SPIE. 2000. V. 4101, P 117121.

75. Milyukov V. Baksan Laser Interferometer Observations: Gravitational and Astrophysical Aspects // In: Gravitation and Astrophysics. Ed. by Liao Liu. Word Scientific. 2000. P. 89-103.

76. Takemoto S., Araya A., Akamatsu J. et al. A 100 m laser strainmeter system installed in a 1 km deep tunnel at Kamioka, Gifu, Japan // Jour, of Geodynamics. 2004. V. 38 P. 477.

77. Wu Shu-Chao, Chen Fan, Fan Shu-Hua, Luo Iun. Phase leading of temperature variation in cavity caused by heat conduction between air and rock. // Chit.phis.lett. 2003. № 12, P. 20.

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

Мясников, Андрей Владимирович
кандидата физико-математических наук
Москва, 2012
ВАК 25.00.10

Диссертация

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

Автореферат

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