Анализ изменений электропроводности геоэлектрической среды на Юго-Восточной Камчатке

Размещено на http://www.allbest.ru/

УДК 550.837

Институт вулканологии и сейсмологии ДВО РАН

ИЗМЕНЕНИЯ ЭЛЕКТРОПРОВОДНОСТИ ГЕОЭЛЕКТРИЧЕСКОЙ СРЕДЫ НА ЮГО-ВОСТОЧНОЙ КАМЧАТКЕ

Мороз Ю.Ф.

Мороз Т.А.

Большой научный и прикладной интерес представляют результаты изучения динамики глубинной электропроводности на Тихоокеанском побережье Камчатки, где по долгосрочному прогнозу ожидается катастрофическое землетрясение с М? 8.5 [9]. Здесь с 2001 года осуществляется мониторинг глубинной электропроводности в пунктах Верхняя Паратунка, Тундровый, Шипунский на побережье Авачинского залива (рис.1). Пункт Шипунский приближен к зоне субдукции, где в большей мере проявились изменения глубинной электропроводности в течение рассматриваемого 13 летнего периода. Выявленные закономерности в поведении электропроводности геологической среды на п-ове Шипунском находят качественное согласие с изменением сейсмичности в районе Камчатки. Полученные результаты и возможная природа изменения глубинной электропроводности рассмотрены в настоящем сообщении.

Рис.1. Схема расположения пунктов наблюдений электромагнитного поля Земли.

На врезке изображено положение района исследований.

1 - район исследований; 2 - глубоководный жёлоб; 3 - вулканы; 4 - пункты наблюдений; 5 и 6 - продольный и поперечные разломы на п-ове Шипунский [3 ]. электропроводность литосфера сейсмичность геологический

Методика наблюдений и обработки данных электромагнитного мониторинга

На всех пунктах приемные каналы регистрации электромагнитного поля ориентированы по направлениям С-Ю и В-З, а также по направлениям осей геоэлектрической неоднородности среды. Длина электрических диполей принята различной для регистрации сигналов от внутриземных и ионосферных источников электромагнитного поля.

Эта длина определяется входным импедансом среды на различных частотах и интенсивностью магнитного поля [6].] Эксплуатация пунктов наблюдений и предварительная обработка данных выполняется Камчатским филиалом Геофизической службы РАН. [7]. Более глубокая обработка полученной информации осуществляется с использованием специальной программы, создающей входные данные для расчётов передаточных магнитотеллурических функций в широком диапазоне периодов. Эти входные данные, созданные автоматически, представляют собой синхронные временные ряды электрического и геомагнитного полей заданной продолжительности. По ним для различных временных интервалов определяется магнитотеллурический импеданс [1]. Для определения тензора импеданса и его главных значений использована программа, алгоритм которой основан на известной робастной методике, разработанной Дж.Ларсеном. [10].

Интерпретация временных рядов электропроводности литосферы

Обратимся к среднегодовым значениям электрического сопротивления, полученным в результате осреднения 5 суточных значений кажущегося сопротивления и фазы импеданса. Графики данных параметров изображены на рис. 2.

Рис.2. Временные ряды среднегодовых значений кажущегося электрического сопротивления и фазы импеданса на периоде 900 с в пп.Верхняя Паратунка (а), Тундровый (б) и Шипунский (в).

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

В п. Верхняя Паратунка среднегодовые значения кажущегося сопротивления и фазы импеданса практически не меняются. Отдельные отклонения импедансов и фаз от графиков находятся в пределах ошибок их определения. В пп. Шипунский и Тундровый проявились изменения электропроводности геологической среды по главным направлениям на протяжении 13 лет. В большей мере они выражены в п. Шипунский. Так, по направлению вкрест зоны субдукции кажущееся электрическое сопротивление меняется от первых сотен до тысячи Омм. При этом вдоль зоны субдукции оно закономерно уменьшается от 5 тысяч до 3 тысяч Омм. Значения фазы импеданса меняются почти на 20 градусов, что в несколько раз превышает ошибку её определения.

Анализ показывает, что изменение электропроводности геологической среды не связано с метеорологическими факторами. Также нельзя объяснить его влиянием внешнего ионосферного, магнитосферного источника. Важно отметить, что для определения магнитотеллурического импеданса использованы вариации электрического и геомагнитного полей с когерентностью не менее 0.8. В пользу отсутствия влияния внешнего источника свидетельствуют данные мониторинга электропроводности в пп. Верхняя Паратунка, где кажущееся сопротивление и фаза импеданса практически не меняются на протяжении 13 лет. Можно полагать, что изменение электропроводности геологической среды в п. Шипунский связано с глубинными процессами. Об этом свидетельствует тот факт, что наряду с изменением кажущегося электрического сопротивления также происходит изменение фазы импеданса. Как мы уже отмечали, п-ов Шипунский расположен в зоне субдукции, где ярко выражены процессы, связанные с сейсмичностью. Землетрясения в зоне субдукции вызваны напряжениями, возникающими при погружении океанической плиты под Камчатку. Земная кора в результате такого воздействия может испытывать сжатие или растяжение, что может повлиять на поведение электропроводности земной коры. О возможности такого влияния свидетельствуют результаты сопоставления гистограммы числа землетрясений с К 11 в районе Камчатки и графика кажущегося электрического сопротивления на периоде 900 с по направлению вкрест простирания зоны субдукции (рис.3).

Рис.3. Сопоставление временного ряда поперечного кажущегося сопротивления (а) и гистограммы землетрясений с К ? 11 (б).

Изменение сейсмичности в общем виде выражено в поведении электропроводности геологической среды в пункте Шипунский. Характерно, что с 2009 по 2013 год ежегодно возрастает количество событий в районе Камчатки. Этому возрастанию отвечает увеличение кажущегося электрического сопротивления в несколько раз. Фаза импеданса также изменяется примерно на 20. Такое поведение электропроводности геологической среды можно объяснить ростом напряжённости литосферы, выраженное усилением сейсмичности. Возникает вопрос, каким образом изменение напряжённости литосферы влияет на её электропроводность. Известно, что существенное влияние на электропроводность горных пород оказывает наличие высокоминерализованных растворов в порах, образующих связанные каналы. Такими каналами являются многочисленные трещины и разломы. На полуострове Шипунском выявлены разломы северо-восточного и северо-западного простираний, т.е. по направлениям вдоль и вкрест простирания зоны субдукции [3]. Разломы изображены на рис. 1..

Попытаемся оценить возможные изменения электропроводности геологической среды в связи с различной электропроводностью разломов на п-ове Шипунский. С этой целью выполнено численное моделирование магнитотеллурического поля с использованием трёхмерной пробной модели южной Камчатки, включающей разломы на п-ове Шипунский. Геоэлектрические параметры верхнего слоя на суше приняты из работы [4,5]. Мощность осадочно-вулканогенного чехла составляет 3 км, а удельное электрическое сопротивление 20 Ом • м. В модель включены разломы, которые выявлены по геолого-геофизическим данным [3]. Один из этих разломов проходит вдоль п-ова Шипунский (рис.1).Удельное электрическое сопротивление морской воды принято 0.25 Ом • м. Глубинные геоэлектрические разрезы континента и дна океана приняты из работы [2] для теплового потока изучаемой территории в 60 мВт/м2 и возраста дна омывающих морей и океана в 40 млн.лет [8]. В результате моделирования по программе Макки [11] получены кривые МТЗ по главным направлениям для моделей с проводящими и непроводящими разломами и локально-нормальные кривые. В случае с проводящим разломом вдоль п-ова Шипунский электрическое сопротивление по направлению вкрест сейсмофокальной зоны уменьшится примерно в 2 раза за счёт индукции электрических токов, концентрирующихся в разломе. В то же время продольное сопротивление в направлении по простиранию сейсмофокальной зоны возрастёт в 1.8 раза за счёт гальванического эффекта S, вызванного контрастом электропроводности между блоком, где расположен пункт Шипунский, и примыкающим проводящим разломом. В случае с проводящими разломами поперёк п-ова Шипунский электрическое сопротивление по направлению вкрест простирания сейсмофокальной зоны увеличилось в 1.5 раза за счёт эффекта S, а продольное сопротивление уменьшилось в 1.3 раза за счёт индукционного эффекта, связанного с проводящими поперечными разломами. Таким образом, данные моделирования подтверждают изменения продольного и поперечного электрического сопротивления высокоомных частей литосферы из-за влияния проводящих глубинных разломов. Проводимость разломов определяется степенью их насыщенности минерализованными растворами, которая зависит, по-видимому, от напряжённости литосферы в зоне субдукции. Можно полагать, что изменение электрического сопротивления верхних частей литосферы связано с увеличением или уменьшением степени насыщенности разломов минерализованными флюидами при растяжении или сжатии земной коры при геодинамических процессах в зоне субдукции.

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

1. Бердичевский М.Н., Дмитриев В.И. Модели и методы магнитотеллурики. М.: Научный мир. 2009. 680 с.

2. Ваньян Л.Л. Электромагнитные зондирования. М.: Научный мир, 1997. 219 с.

3. Карта полезных ископаемых Камчатской области. Масштаб 1 : 500000/ Под ред. Литвинова А.Ф., Патоки М.Г., Марковского Б.А. СПб.: ВСЕГЕИ, 1999. Л. XIII.

4. Мороз Ю.Ф. Электропроводность земной коры и верхней мантии Камчатки. М.:Наука, 1991. 181 с.

5. Мороз Ю.Ф., Нурмухамедов А.Г. Магнитотеллурическое зондирование Петропавловского геодинамического полигона на Камчатке // Вулканология и сейсмология. 1998. № 2. С.77-84.

6. Мороз Ю.Ф. О методике поисков краткосрочных предвестников сильных землетрясений в низкочастотном теллурическом поле Камчатки//Физика Земли. 1994.N9. С.88-90.

7. Мороз Ю.Ф., Бахтиаров, В.Ф.,Гаврилов, В.А., Левин, В.Е.,Попруженко С.В. О мониторинге электротеллурического поля с целью прогноза сильных землетрясений на Камчатке //Вулканология и сейсмология.1995. N4-5.С.139-149.

8.Смирнов Я.Б., Сугробов В.М. Земной тепловой поток в Курило-Камчатской и

9. Федотов С. А,, Соломатин А, В., Чернышов С. Д. Долгосрочный сецсмический прогноз для Курило-Камчатской дуги на IX 2011-VIII 2016 гг.; верноятное место, время и развитие следующего сильнейшего землетрясения Камчатки с М ? // Вулканология и сейсмология.2012. № 2 С.3-26.

10. Larsen I.C., Mackie R.L., Manzella A.,et al. Robust smooth magnetotelluric transfer functions // Geophysical Journal International. 1996. V. 124. P. 801-819

11. Mackie K.L. ,Smith.J.T., Madden T.R. Three-dimensional electromagnetic modeling using finite difference equations: the magnetotelluric example // Radio Science. 1994. № 29. P. 923-935.

Размещено на Allbest.ru