Модель электропроводности земной коры

Рис. 3

 В последние годы на Южной Камчатке выполнены наблюдения МТЗ по двум региональным профилям, которые дают возможность в значительной мере уточнить и дополнить существующие представления о структуре электропроводности земной коры. Методика интерпретации и анализ результатов данных МТЗ подробно изложены в работе [17]. Полученные результаты приведены на рис.3. Они сводятся к следующему.

Геоэлектрическая модель содержит коровый проводящий слой, хорошо развитый на большей части площади, за исключением восточного побережья, где он почти отсутствует. Сопротивление слоя от 5-10 Ом^.м на западе увеличивается до 80 Ом^. м в центральной части профиля, при этом его мощность уменьшается до 15 км. Далее к востоку, в районе Авачинско-Корякской группы вулканов, происходит уменьшение сопротивления коры до 15 Ом^. м на глубине 8-35 км. Здесь из-за сильного эффекта S (узкая проводящая зона в осадочно-вулканогенном чехле) параметры корового проводника определяются неоднозначно. Восточнее по профилю сопротивление корового слоя возрастает до 200 Ом^. м, мощность составляет 25 км.  

Относительно природы корового проводящего слоя можно отметить следующее. Электрическое сопротивление верхних частей консолидированной коры 1000-2000 Ом^. м. При нормальном распределении глубинной электропроводности нижние толщи земной коры имеют удельное сопротивление не менее 500 Ом^. м [5,14], поэтому природу проводящего слоя в земной коре Южной Камчатки сопротивлением в десятки-единицы ом на метр можно объяснить наличием жидких флюидов. Грубые оценки показывают, что достаточно десятых долей процента относительно объема флюидов в связанных каналах, чтобы обеспечить уменьшение сопротивления горной породы до десятков - единиц ом на метр. Можно предполагать, что на большей части площади коровый проводящий слой связан с наличием высокоминерализованных гидротермальных растворов. В зоне V (рис.3), где располагается Авачинско-Корякская группа вулканов, коровый проводник может быть обусловлен наличием магматических расплавов. В восточной части профиля отсутствует низкоомный коровый проводник; здесь же отмечена интенсивная положительная аномалия поля силы тяжести.

По данным электропроводности можно оценить пористость и проницаемость пород земной коры Южной Камчатки. Для этого воспользуемся данными работ [5,6]. Согласно закону Арчи

= _ф /k²

где - удельное сопротивление влажной породы, _ф - удельное электрическое сопротивление флюида; k - пористость. При = 2000 Ом^. м значение _ф = 0,3 Ом^. м, пористость верхних частей земной коры до глубины 10 км равна 0,8%. Эта величина близка к данным по северо-западу США, полученным Л.Л.Ваньяном. О проницаемости верхних частей консолидированной земной коры можно судить по зависимости удельного сопротивления влажных пород от их проницаемости [36]. При сопротивлении коры в 1400 Ом^. м проницаемость ее ~10^-21 - 10^-22 м². При такой проницаемости время проникновения флюида из корового слоя в верхние горизонты земной коры будет ~ 50 млн.лет. Возраст вулканизма Авачинско-Корякской группы вулканов составляет первые миллионы лет. По-видимому, более быстрое проникновение флюидов в верхние горизонты коры в данном районе определяется наличием глубинных разломов. Уже указано, что в зоне V электропроводность коры между осадочно-вулканогенным чехлом, заполняющим грабен, и коровым слоем определяется крайне неоднозначно. Не исключается, что здесь существует вертикальный проводник, обусловленный глубинными разломами, по которым магматические расплавы поступают на поверхность. Возможно, с этими разломами связана современная вулканическая деятельность Авачинско-Корякской группы вулканов.

Объемная скоростная модель

Несмотря на значительный объем сейсмической информации, в настоящее время практически отсутствует скоростная модель Камчатки, включающая области омывающих ее акваторий. Сейсмические наблюдения ГСЗ-КМПВ выполнены в наземной части полуострова с использованием методики кусочно-непрерывного профилирования с большими пропусками в наблюдениях [2,8,9 и др.]. Они были сосредоточены, в основном, в районах Авачинско-Корякской и Ключевской групп вулканов, т.е. территория Камчатки в целом слабо изучена этими методами. Кроме того, длина сейсмических профилей не позволила выполнить корректные оценки скоростных параметров по поверхности Мохоровичича (М). По этой причине остается открытым вопрос о мощности коры под Камчаткой и сейсмических характеристиках зоны перехода кора-мантия. Более длинные ветви годографов получены по наземно-морскому профилю, пересекающему Авачинский залив [8,27 и др.], однако полученные результаты по ряду причин (в основном, редкой сети наблюдений) также крайне противоречивы [1]. Кроме того, следует отметить большие сложности в выполнении сейсмических наблюдений за сигналами от взрывов и необходимость огромных финансовых вложений, что практически делает невозможным их осуществление в настоящее время. Вместе с тем на Камчатке существует детальная сеть сейсмологических станций, которые регистрируют тысячи землетрясений, просвечивающих земную кору и верхнюю мантию под полуостровом и прилегающими акваториями. Эти данные также мало использованы для оценок скоростных параметров. Это, безусловно, было связано со значительными ошибками в определении координат гипоцентров землетрясений и отсутствием вычислительного аппарата для решения обратной трехмерной задачи. Видимо, по этой причине существуют лишь единичные работы, благодаря которым получены первые схемы распределения скорости Vp под Камчаткой и акваториями, а также в области пересечения Курило-Камчатской и Алеутской островных дуг [3,22,28]. Эти работы продемонстрировали четкую корреляцию скоростных аномалий с глубинной тектоникой региона и подтвердили перспективность развития данного направления на Камчатке.  

Здесь приведены результаты томографического восстановления объемной скоростной модели земной коры и верхней мантии Камчатки в ее южной части. Для этого использовались времена вступлений P-волн от близких землетрясений, которые, в основном, сосредоточены в сейсмофокальном слое. Они взяты из инструментального каталога землетрясений Камчатки (в окончательной обработке) за период 1992-1998 гг., составленный Камчатской опытно-методической сейсмологической партией (КОМСП) ГС РАН. Количество зарегистрированных событий в каталоге ~ 13000, сейсмических лучей, просвечивающих кору и мантию от гипоцентров до станций, ~ 100000. Однако объем использованной информации был значительно ограничен точностью входных данных, которая в общем достаточно низкая. Причина и области ошибок детально проанализированы в [30]. Для расчетов были выбраны землетрясения, точность определения времени в очаге для которых не достигала 0,4 сек, точность определения координат эпицентров не превышала по горизонтали +20 км и по вертикали - +10 км. Также накладывались ограничения по энергетическому классу (К> 11), четкости вступления P и S волн, числу станций, участвующих в определении гипоцентров, Т₀ и др. В обработку принимались события с глубиной, не превышающей 45 км. В связи с этим глубина проникновения сейсмических лучей для выбранных событий не превышала 90 км. Эта цифра ограничивает скоростную модель по глубине. Исходя из перечисленных выше ограничений, из каталога было выбрано только 1400 землетрясений. На основе анализа осредненных наблюденных годографов для выбранных событий выполнена оценка средней мощности земной коры под Камчаткой и весь объем литосферы разбит на три слоя:

 0-20 км (верхняя кора), средняя скорость в слое 5,5 км/c;

 20-35 км (нижняя кора), средняя скорость в слое 7,0 км/с;

 35-90 км (слой верхней мантии), средняя скорость в слое 8,0 км/с.  

Исходя из точности входных данных, горизонтальный размер блоков, в которых определяется значение невязки скорости ( Vp), составляет 50х50 км. Методика решения томографической задачи изложена в [18]. Многократным перемещением рассчетной схемы значения Vp получены в ячейках, составляющих по горизонтали ~ 17х17 км. Расчеты выполнены совместно с сотрудниками ОИФЗ РАН И.А.Саниной и М.Ю.Степановой.

Рис. 4

 Фрагменты восстановленной трехмерной скоростной модели (для района Южной Камчатки) представлены на рис.4. Они иллюстрируют существенно новую информацию о распределении скоростных неоднородностей изученного объема среды. Поля невязок свидетельствуют о резкой скоростной дифференциации литосферы под Камчаткой, анизотропии упругих параметров, различии структурных планов скоростных аномалий земной коры и слоя верхней мантии. В целом поля невязок в слоях можно разбить на отдельные зоны (блоки), которые характеризуются различными скоростными свойствами. На рис. 4а нанесены примерные границы этих блоков, которые в пределах Южной Камчатки мы будем именовать согласно [33].  

В верхней коре (0-20км) значения скорости Vp изменяются от 5,2 до 6,0 км/с. Скоростные аномалии, имеющие изометрическую форму, вытянуты в цепочки, ориентированные вкрест простирания Камчатки (рис.4А). Они имеют продолжение в акватории. Северная из них, которая характеризуется наиболее высокой скоростью, располагается в районе мыса Шипунский. В южной части выделяется цепочка аномалий пониженной скорости. Она пересекает сушу и продолжается в шельфовой зоне в районе Мутновско-Вилючинского блока. В акватории (на продолжении Асачинского блока) она граничит с аномалией повышенных значений скорости.  

В нижней коре (20-35км) выделяется высокоскоростная аномалия в районе Шипунского блока (рис.4В). Здесь отмечены наиболее высокие значения скорости Vp (6,7км/с). Южнее Авачинского блока слой нижней коры является низкоскоростным, при этом низкие значения скорости характерны как для района Восточно- Камчатского вулканического пояса, так и Центральной депрессии и шельфовой зоны. Следует отметить корреляцию высокоградиентной скоростной зоны в южной части Авачинского залива по всей мощности коры. Она располагается примерно на продолжении Асачинского блока в океан. Последний, согласно [33], расположен на границе зоны перехода Южной Камчатки к структурам Центральной части полуострова. Согласно полученным данным она продолжается в акватории.

В слое верхней мантии (35-90км) выделяются две области (блоки), различные по скоростным свойствам (рис.4С). Граница между ними проходит ~ на уровне Вилючинско-Мутновской группы вулканов. Северный блок характеризуется пониженными значениями скорости. Здесь располагаются действующие вулканы (в частности, Авачинко-Корякской группы). Слой верхней мантии в южном блоке в целом более высокоскоростной. При этом самые высокие значения скорости отмечены в южной части акватории Авачинского залива и самой южной оконечности полуострова.

Рис. 5

 Анализируя поля невязок в целом по всей изученной мощности литосферы, можно сделать вывод, что высокоградиентную зону, разделяющую Южную Камчатку на суше и в акватории на две различные по упругим свойствам области, можно характеризовать как сквозькоровый глубинный разлом, имеющий продолжение в мантии. Обращают на себя внимание некоторые закономерности в распределении вулканических групп на скоростных схемах. Так в верхней коре они, как правило, приурочены к границам раздела контрастных скоростных неоднородностей и не обязательно связаны с низкоскоростными аномалиями. В связи с этим можно согласиться с выводами, сделанными в [34], что в соотношении вулканизма с тектоникой примат принадлежит последней. В поле невязок скорости верхней мантии тоже не наблюдается связи вулканов с определенными аномалиями. В северной части они связаны с низкоскоростной мантией (~7,3-7,5км/с), а в южной-нормальными значениями Vp (~8,0км/с).  

Обратимся к сопоставлению полученных результатов с данными других геолого-геофизических исследований. Северо-западное простирание цепочек скоростных аномалий в верхней коре находит отражение на суше в виде Петропавловск-Малкинской зоны поперечных дислокаций [7]. На поперечное простирание третичных структур в свое время обратил внимание А.Н.Тихонов [7]. Этот вопрос широко дискутировался в научной литературе, и предполагается, что северо-западное простирание меловых и более древних структур не является доминирующим. Однако, из приведенных выше скоростных моделей видно, что в изученном блоке литосферы Южной Камчатки и прилегающей акватории развиты структуры северо-западного простирания, которые, возможно, имеют продолжение в Тихий океан. Следует отметить, что аномалиям пониженных скоростей в земной коре на юге Авачинского залива отвечают аномалии пониженных значений поля силы тяжести, обусловленных отрицательной структурой в рельефе океанического дна (рис.5) [21].

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

В верхней мантии область пониженных скоростей в районе вулканов Авачинско-Корякской группы выражена также повышенной электропроводностью, что связывается с подъемом астеносферного слоя до глубины 50 км (см. рис.1). Аномальные свойства верхней мантии на глубинах 35-90 км могут быть обусловлены наличием магматических расплавов.