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

Васильева, Евгения Владимировна

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

Автореферат диссертации по теме "Геохимия магматогенно-гидротермальных систем на границе пластично-хрупкого перехода в земной коре"

005011583

ВАСИЛЬЕВА Евгения Владимировна

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

Специальность 25.00.09 - геохимия, геохимические методы поисков полезных ископаемых

Автореферат

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

1 6 0ЕВ 2012

ИРКУТСК 2012

005011583

Работа выполнена в Федеральном государственном бюджетном учреждение науки Геологическом институте СО РАН (ГИН СО РАН)

Научный руководитель:

доктор геолого-минералогических наук

Жатнуев Николай Сергеевич

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

доктор геолого-минералогических наук, профессор

Скворцов Валерий Александрович кандидат геолого-минералогических наук

Бычинский Валерий Алексеевич

Ведущая организация:

Федеральное государственное бюджетное учреждение науки Институт вулканологии и сейсмологии ДВО РАН

Защита состоится " " -гса^т^ 2012 г. в ~/0.££> часов на заседании диссертационного совета Д 603.059.01 при Федеральном государственном бюджетном учреждении науки Институте геохимии им. А.П. Виноградова СО РАН по адресу: 664033, г. Иркутск, ул. Фаворского 1а.

С диссертацией можно ознакомиться в научной библиотеке Федерального государственного бюджетного учреждения науки Института геохимии им. А.П. Виноградова СО РАН, по адресу: 664033 Иркутск, ул. Фаворского 1а.

Автореферат разослан " " (2012 г. e-mail: korol@igc.irc.ru.

Ученый секретарь диссертационного совета к.г.-м.н.

Г.П. Королева

ВВЕДЕНИЕ

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

Целью исследования является создание согласованной модели магматогенно-гидротермальной системы с привносом глубинного и метеорного вещества вблизи зоны пластично-хрупкого перехода (ПХП). Задачи работы:

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

2. Разработка сценариев эволюции МГС в зависимости от положения магматического очага относительно зоны ПХП; расчет физических параметров в каждом узле модельной сетки в соответствии со сценарием.

3. Разработка физико-химической многорезервуарной модели МГС с расчетом исходных равновесных минеральных парагенезисов, а также химического состава пород и флюида в каждом узле модельной сетки:

а) для траектории подъема глубинного флюида в пластической области с прорывом в область хрупких деформаций;

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

4. Разработка методики тектонофизического эксперимента на оптически активных средах в поляризованном свете; проведение серии экспериментов; разработка программного обеспечения для цифрового анализа результатов экспериментов.

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

1. Возникновение и эволюция МГС возможны по трем сценариям:

а) внедрение магмы в область хрупкости пород ведет к смешению ювенильного и метеорного вещества при промывании быстро остывающего очага гидротермальными растворами;

з А

б) при локализации очага в области пластично-хрупкого перехода имее место быстрое остывание апикальной части массива и частично смешение ювенильного и метеорного материала;

в) если очаг локализуется в области пластичности, загрязнение массив; метеорным веществом практически исключается, а в области хрупкост! формируются низкотемпературные гидротермальные системы метеорным составом растворов.

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

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

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

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

2. На основе этой концептуальной модели разработана комплексна компьютерная модель физико-химической эволюции МГС в зоне изменени реологических свойств пород с привносом глубинного и метеорног вещества.

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

Фактический материал. Данные по химическому составу природных сред базы термодинамических данных были любезно предоставлень непосредственно Н.С. Жатнуевым, В.И. Васильевым (Геологический институ СО РАН, Улан-Удэ), К.В. Чудненко, В.А. Бычинским (Институт геохимии им А.П. Виноградова СО РАН, Иркугск) и С.Н. Рычаговым (Институт вулканологт ДВО РАН, Петропавловск-Камчатский).

В работе также использовался фактический материал, самостоятельн полученный автором в результате экспериментальных исследований проводимых при финансовой поддержке гранта НШ-2284.2003.5 и Президиумо СО и ДВО РАН (проект № 117-09-П-СО-08-006). Кроме этого, фактически: материалом работы явились общепризнанные литературные данные по геохими и минеральному составу горных пород и их термодинамическим свойствам. Практическая значимость работы.

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

2. Разработанные модели МГС с учетом изменения плотности и давления

флюида на геохимическом барьере ПХП позволяют по-новому оценивать продуктивность современных геотермальных месторождений, объяснить природу их строения и структуры. 3. Созданная база физико-химических данных для различных сценариев эволюции МГС, каждый из которых включает 80 минеральных парагенезисов с водным раствором и газовой фазой (итого 240 возможных равновесных ассоциаций) позволит поднять будущие модельные работы на более детальный уровень. Личный вклад автора. Автором разработаны концептуальные модели эволюции МГС, в чем оказали помощь д.г.-м.н. Н.С. Жатнуев (ГИН СО РАН) и д.г.-м.н. С.Н. Рычагов (ИВиС ДВО РАН). Автором реализованы геометрический, физический и динамический этапы моделирования. Также автором были выполнены все физико-химические расчеты на ПК «Селектор» с применением сторонних программных продуктов (Васильев, ГИН СО РАН; Арискин, ГЕОХИ РАН; Персиков, ИЭМ РАН). Автором также лично разработана методика тектонофизических экспериментов, изготовлены модельные установки, проведены эксперименты и анализ полученного материала, выполнена интерпретация результатов.

Публикации и апробация работы. Основные положения диссертации опубликованы в 19 печатных работах (в том числе 2 статьи в рецензируемых журналах «Геоинформатика» и «Литосфера» и 1 статья принята в печать в рецензируемый журнал «Вулканология и сейсмология»), а также докладывались на Всероссийской конференции, посвященной 80-летию чл.-корр. РАН Ф.П. Кренделева (Улан-Удэ, 2007), Региональной молодежной конференции, посвященной 50-летию Бурятского научного центра (Улан-Удэ, 2008), Региональной конференции молодых ученых «Современные проблемы геологии, геохимии, геоэкологии Дальнего Востока России» (Владивосток, 2008), Всероссийском научном совещании «Геодинамическая эволюция литосферы Центрально-Азиатского подвижного пояса: от океана к континенту» (Иркутск, 2008, 2009, 2010), Всероссийском совещании «Разломообразование и сейсмичность в литосфере: тектонофизические концепции и следствия» (Иркутск, 2009), Всероссийской молодежной конференции «Структура, вещество и история литосферы» (Сыктывкар, 2009), Международной конференции ЮНЕСКО «Глобальные и региональные проблемы устойчивого развития мира» (Улан-Удэ, 2010), Всероссийской молодежной конференции «Геология Западного Забайкалья» (Улан-Удэ, 2011), Международного симпозиума «Large Igneous Provinces of Asia, Mantle Plumes and Metallogeny» (Иркутск, 2011), Всероссийской молодежной конференции «Современные проблемы геохимии» (Иркутск, 2011), Всероссийской конференции «Граниты и процессы рудообразования» (Москва, 2011), научных сессиях и семинарах Геологического института СО РАН, Института геохимии им. А.П. Виноградова СО РАН, Бурятского научного центра и Бурятского Государственного университета.

Объем работы. Диссертация состоит из введения, пяти глав, заключения и

списка литературы общим объемом 162 страниц печатного текста, 13 таблиц и 106 иллюстраций. Список использованной литературы включает 119 наименований.

Благодарности. Автор считает своим долгом выразить сердечную благодарность в первую очередь своему научному руководителю - д.г.-м.н. Николаю Сергеевичу Жатнуеву за терпеливую плодотворную помощь на протяжении всего времени учебы и работы над диссертацией. Автор душевно благодарит своих учителей: д.г.-м.н. Д.И. Царева, д.г.-м.н. А.Г. Миронова, д.г.-м.н. А.М. Плюснина, к.г.-м.н. В.И. Васильева, к.г.-м.н. Д.А. Орсоева, к.г.-м.н. Г.С. Риппа, кх.-м.н. Т.Т. Врублевскую и других преподавателей. В подготовке и проведении экспериментов непосредственную помощь оказали Г.Д. Санжиев, М.В. Бурцева и к.г.-м.н. Б.Б. Дамдинов, за что им автор очень благодарна. Особую признательность хотелось бы выразить сотрудникам Лаборатории физико-химического моделирования Института геохимии им. А.П. Виноградова СО РАН дх.-м.н. К.В. Чудненко и к.г.-м.н. В.А. Бычинскому за предоставленные версии ПК «Селектор» и специальные базы термодинамических данных. ГЛАВА 1. ОБЪЕКТ И МЕТОДЫ ИССЛЕДОВАНИЯ

1.1. Магматогенно-гидротермалъные системы. Предлагается собственное определение: МГС - это физико-химическая система горных пород с циркулирующим в них флюидом, в которой под воздействием тепла магматического очага, привноса глубинных и метеорных флюидов происходят процессы тепломассопереноса и минералообразования. Кроме стандартных составляющих МГС: магматической системы и системы вмещающих пород, вводится подсистема порового флюида, представляющая собой группу подвижных фаз «раствор + газ» с переменным фазовым и химическим составом.

1.2. Численное физико-химическое моделирование. Доминирование физико-химического направления в геохимических исследованиях закономерно следует из обязательной эволюции вещества в геологических объектах (Жариков, 2005).

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

1. Проблема исходных петрофизических и термодинамических данных.

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

3. Алгоритмическое и методологическое обеспечение задач численного физического и физико-химического моделирования.

4. Создание комплекса рабочих программ, необходимых для численного серийного решения широкого круга геохимических задач.

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

природных объектов.

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

1.3. Методика комплексного компьютерного моделирования геохимических объектов включает четыре взаимосвязанных этапа: геометрический, физический, физико-химический и динамический. На геометрическом этапе обосновывается размерность моделируемой системы, производится разбиение ее пространства на подсистемы и объединение их в зоны, соответствующие природным обособлениям в составе моделируемого объекта. Физический этап предусматривает расчет физических параметров каждой подсистемы, исследование взаимодействия подсистем между собой и установку начальных и граничных условий в системе. В качестве расчетных параметров выступают температура и давление. Физико-химический этап предусматривает расчет исходных равновесных составов подсистем, которые рассматриваются как физико-химические резервуары, изначально не взаимодействующие друг с другом. Расчеты предлагается осуществлять методом минимизации термодинамических потенциалов на базе ПК «Селектор» (Чудненко, 2005, 2007, 2010), описание которого также приводится в данном разделе. Входными данными для таких расчетов будут являться рассчитанные на физическом этапе температура и давление в каждом резервуаре, а также набор и концентрации независимых компонентов. Результатом этапа должны являться расчетные равновесные химические составы каждого резервуара (подсистемы). На динамическом этапе начальными и граничными условиями являются элементы дискретизации временной области. Главное из рассматриваемых событий при исследовании динамики объекта - это временные изменения, поэтому в первую очередь важно предусмотреть методы объекта как функции вектора времени. При этом учитываемые методы можно разделить на две категории: межзонные методы контролируют изменение свойств зон модели, а межрезервуарные методы описывают взаимодействие резервуаров и изменение их свойств во времени.

1.4. Тектонофизические исследования. После обзора состояния тектонофизических исследований и современных взглядов на петрофизическое строение земной коры описываются использованный в работе экспериментальный поляризационно-оптический метод исследования

напряжений в оптически активных средах (Осокина, 1963; Осокина, Цветкова, 1980) и авторский метод анализа результатов экспериментов. Первичными результатами поляризационно-оптического метода являются фото- и видеоизображения. На таких изображениях в строго количественной интерпретации нуждаются напряжения модельного материала, выраженные оттенками определенных цветов. Необходимо учитывать, что фотографии делаются при различных внешних условиях: освещении, влажности, разных коэффициентах преломления стенок модельной емкости, различными фотокамерами. Таким образом, средство анализа должно «уметь» разлагать изображение на цвета, вычислять интенсивность каждого цвета к иметь некий «эталон интенсивности цветов» для независимости от внешних условий. В связи с этим автор работы принял участие в создании программного продукта УЫ! МишОгарЬ. Продукт успешно апробировался при интерпретации экспериментов. Программный продукт соответствует всем вышеприведенным требованиям. Последнее требование достигается контролем соотношения четырех эталонных цветов в конкретных внешних условиях, при которых производится съемка эксперимента. Это белый, красный, синий и зеленый цвета, оцифрованные по КОВ-модели с 256-разрядной интенсивностью цветов (0-255).

ГЛАВА 2. КОНЦЕПТУАЛЬНАЯ МОДЕЛЬ ЭВОЛЮЦИИ МГС 2.1. Реологическая модель. Известно, что в земной коре на глубине 6-20 км хрупкие деформации сменяются пластическими (Ващилов, 1984; Иванов, 1970, 1990). При этом в пределах хрупкой коры флюид находится в открытых трещинах под гидростатическим давлением, а ниже - в изолированных полостях и порах в условиях максимального литостатического давления. На рис. 1 показан возможный интервал ПХП.

Рис. 1. Соотношение литостатического и гидростатического давлений с глубиной в связи с реологией земной коры (по Киссину, 2001 с дополнениями). Р! - I

гидростатическое давление флюида в трещинах в области хрупкого состояния пород; Р2 -давление флюида в порах, равное литостатическому в области пластического состояния; РЗ-зона перехода от Р! к Р2, в зависимости от физических свойств пород. На разных глубинах перехода различен и перепад давлений. I - зона открытой трещиноватости в области хрупких деформаций; II -зона закрытой трещиноватости В верхней, хрупкой и трещиноватой части земной коры, растворы находятся под собственным гидростатическим давлением и при плотности флюида 1.0 г/см' график давления будет соответствовать линии Р,. Реальные давления могут

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

2.2. Механизм движения полостей в пластичной среде. Давление флюида, находящегося в замкнутом поровом пространстве ниже переходной зоны, соответствует литостатическому давлению, но оно должно существенно превышать таковое в голове протяженных по вертикали полостей. Чем протяженнее полость, тем значительнее должна быть разница между литостатическим и флюидным давлением вверху полости; эту разницу будем называть избыточным давлением флюида (ИДФ) (рис. 2).

Давление, МП«___

1 ь, 2 -—^-Н \ р.'.. р,'.. Рыг.1 ч

1, 1V 1> ! ^ч /У] V\ \ р; \ а/ 4 К4 V

П \ 1 .'НИИ \ \ [У/___________________

р{\

О 10 20 30 Гидростатическое давление столба флюида в трещине (Ял).МПа

Рис. 2. Схематическая диаграмма (по Жатнуеву, 2005), поясняющая образование ИДФ. И - расстояние от поверхности земли до головной части полостей 1 и 2; - до хвостовой (нижней) части полости 1 и - до хвостовой части полости 2; /| и /2 - вертикальная протяженность полостей 1 и 2; Рк\ и Л2- гидростатическое давление столба флюида в треполостях 1 и 2 (при плотностях флюида 0.5 и 1.0 г/см3); Р/ - величина литостэтического давления в головной части полостей 1 и 2; Р' и Л - величины литостатических давлений в хвостовой части полостей 1 и 2 (при плотности породы 2.7 г/см3); Р'тг.1,о, Р'„„гЛ5 - величина избыточного давления флюида в головной части полости 1 при плотностях флюида 1.0 и 0.5 г/см1; Р2тг.1л Р2 mr.ii.s- величина избыточного давления флюида в головной части полости 2 при тех же плотностях флюида; Р', и Г/абсолютные давления флюида на стенки полостей 1 и 2 в головной части при плотностях флюида 1.0 и 0.5 г/см3.

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

наблюдаются зоны (участки, блоки, горизонты) с контрастными петрофизическими свойствами пород. Это связано с наличием на различных глубинах, вплоть до 10-12 км (Кольская сверхглубокая..., 1984), зон разуплотнения пород, насыщенных газово-жидкими флюидами. В зависимости от свойств пород МГС, возбуждаемые магматическим очагом, могут развиваться согласно трем сценариям (рис. 3).

Рис. 3. Внедрение магматических очагов в разные реологические зоны земной коры и схемы развития конвективных потоков гидротерм в зависимости от положения источника тепла. А - локализация очага магмы в зоне хрупких деформаций; В - локализация магмы на границе между зонами упругих и пластических деформаций; С - локализация магматического очага в зоне пластичных деформаций. Большими стрелками показана конвекция растворов в зоне хрупких деформаций, маленькими - кондуктивный тепловой поток от магматического резервуара, расположенного в зоне пластичности пород.

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

ГЛАВА 3. МОДЕЛЬ ПОДЪЕМА ФЛЮИДА ЧЕРЕЗ ПХП 3.1. Описание модели. Очевидно, что граница пластичной и хрупкой сред, которая имеет резко различные давления флюида и петрофизические свойства пород в своей верхней и нижней части, является принципиальным геохимическим барьером. Наша модель позволяет оценить эволюцию флюида и вмещающих пород на этом геохимическом барьере. При переходе к реальным значениям глубины, давления и плотностей и при учете температурных условий были разработаны различные варианты эволюции заполняющего полости флюида. На их основе построены шесть модельных сценариев массопереноса

для многорезервуарной термодинамической модели подъема глубинного флюида к поверхности через ПХП. Глубины ПХП приняты равными 5 км (сценарии 1-3), 6.6, 6.9 и 8 км (сценарии 4-6) при локализации магматического очага на глубине 8 км. Сценарии характеризовались различной удаленностью тренда вертикального подъема флюида от оси магматического очага и, соответственно, различным перепадом плотности и давления флюида. В зоне пластических деформаций флюид, находящийся под литостатическим давлением, последовательно поднимался до границы ПХП. В зоне хрупких деформаций последовательный проточный реактор усложнялся привносом метеорных вод из открытых трещин под гидростатическим давлением. Набор независимых компонентов модели отвечал системе 8^А1-Ре-М§-Са-Маг-К-Н-О-С-СЬ-в. Составы резервуаров рассчитывались как андезитобазальтовые '(Богатиков и др., 1987). Состав исходного флюида был принят аналогичным составу глубинных термальных вод о. Кунашир по (Справочник по геохимии, 1990). Состав метеорных вод был усреднен из данных по составам морской воды со средней соленостью 35%о. Термодинамическое моделирование проводилось на базе ПК «Селектор». При расчетах использовался метод последовательного проточного реактора для группы подвижных фаз «раствор + газ». Равновесия в резервуарах рассчитывались методом минимизации энергии Гиббса. 3.2. Результаты моделирования. Главным результатом моделирования явилось подтверждение геохимического барьера ПХП и определение форм осаждения вещества флюида. На рис. 4, 5 и 6 показаны соответственно расчетные изменения масс свободного флюида, фазы водного раствора и газовой фазы в зависимости от глубины. Показательно изменение концентраций компонентов фазы водного раствора.__________

а свободного флюида при глубине ПХП 5 км

ГГ;

Масса свободного флюида при глубине ПХП 6.6 - 8 км

—9—Сечение 1 ■ ■■ Сечение 2 «-■Сечение 3

«ООП !ЪОО 2000 ?5СЮ 2000 Масса, г

—Ч— Селение 1 (ПХП 8 КЫ) —«— Сечение 2 (ПХП 6.9 Ш) —*—Сечение 3 (ПХП 6 6 км)

12СБ 1Ь00 '¿ООО

Мзсса, г

Рис. 4. Расчетное изменение массы свободного флюида по глубине.

Рис- 5. Расчетное изменение массы водного раствора по глубине.

Масса водного раствора при глубине ПХП 5 км

Сечение 1 —•—Сечение 2 —*—Сучение ^

Масса газовой фазы при глубине ПХП 5 км

Сечение ( -■•—Сечение 2 —»—Сечение А

1000 2000 3000

Масса, г

Масса водного раствора при глубине ПХП 8.6 - 8 км

Сечение 1 (ПХП 8 км) -■—Сечение 2 (ПХП 6.9 км) Сечение а (ПХП 6.6 им)

Масса газовой фазы при глубине ПХП 6 6 - 8 км

1000 2000 Масса, г

Рис. 6. Расчетное изменение массы газовой фазы по глубине. Наибольшее падение концентраций характерно для растворенных форм НС1, вЮ* 802, Н^, С02, ИаОН, СО, КС1 и №Н8Ю3. Приблизительно на одном уровне остаются концентрации 02, 8032~, MgC\+ и НС03~. Концентрации в растворе Ре02~, РеОН2+, РеС12 и на геохимическом барьере ПХП, напротив, возрастают. Содержание минералов во вмещающей породе при прохождении флюида через геохимический барьер ПХП также изменяется. Во всех сценариях отлагаются а-кварц, кристобалит, халцедон или тридимит. Первый сценарий характеризуется интенсивным гидротермальным изменением пород в зоне ПХП. Значительно повышается содержание ставролита, появляются Ре-хлорит,

■ектолит и гейландит; исчезает альбит и гематит.Во втором и третьем сценариях интенсивно отлагаются цеолиты, например, K-Na-мерлиноит; ставролит и пектолит растворяются, магнетит сменяется гематитом. Четвертый сценарий на геохимическом барьере качественно повторяет первый, отличаясь интенсивностью и исчезновением диопсида, а в пятом и шестом сценариях, при качественном сходстве со вторым и третьим, содержание диопсида изменяется мало (на 2-3%).

3.3. Основные выводы: 1) ПХП является принципиальным геохимическим барьером, расположенным геометрически ниже конвективной ячеи МГС и, наряду с расположением магматического очага, определяет как исходный состав всей системы, так и режим поступления в нее глубинного вещества; 2) Переход от литостатического к гидростатическому давлению, при прорыве из зоны пластичности в зону хрупкости, приводит к адиабатическому расширению флюида со значительным охлаждением и осаждением растворенного вещества.

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

ГЛАВА 4. КОМПЛЕКСНАЯ КОМПЬЮТЕРНАЯ МОДЕЛЬ МГС Согласно методике, на геометрическом этапе моделирования была определена размерность модели и размер системы - двумерный вертикальный разрез МГС, включающий зоны: нижнюю - пластических деформаций, верхнюю - хрупких деформаций и магматический очаг, находящийся на разной глубине согласно концептуальным сценариям. Система считалась симметричной относительно оси магматического очага, поэтому мы ограничились полуплоскостью с четырьмя вертикальными сечениями на глубину 20 км. В каждом из четырех сечений было определено 20 подсистем размером 1><1 км. Глубина ПХП принята равной 8 км. На физическом этапе были рассчитаны распределение температурных полей (рис. 7), гидростатическое давление для хрупкой зоны и литостатическое - для пластичной (рис. 8). Плотность флюида рассчитывалась по эмпирическому уравнению, полученному экспоненциальной аппроксимацией усредненных данных для воды и растворов угле- и кремнекислоты (Burncham etc., 1969; Рид и др., 1982): р>ш = 1.0616е °т7Г ■

На физико-химическом этапе был определен набор независимых компонентов модели (система Si-Al-Fe-Mg-Ca-Na-K-H-O-C-S-Cl-F-Zn-Cu-Pb-Sn). Термодинамическое моделирование проводилось на базе ПК «Селектор». Исходные равновесные парагенезисы с водным раствором и газовой фазой рассчитывались методом минимизации энергии Гиббса. В целом по первому сценарию расчет показал возможность равновесного существования 50 минералов, 4 компонентов газовой фазы и 97 компонентов водного раствора (для хрупкой зоны); и 24 минералов, 2 компонентов газовой фазы и 96 компонентов водного раствора (для пластичной зоны). По второму сценарию

способны равновесно существовать 46 минералов, 5 компонентов газовой фазы и 91 компонента водного раствора (для хрупкой зоны); и 35 минералов, 9 компонентов газовой фазы и 106 компонентов водного раствора (для пластичной зоны). По третьему сценарию определились 29 минералов и 73 компонента водного раствора (для хрупкой зоны, газовая фаза отсутствует); и 39 минералов, 11 компонентов газовой фазы и 106 компонентов водного раствора (для пластичной зоны).

Сценами I

Сценашй II

Сценарий 111

Об" X 5:

Ю >»

х. %

>ч

э-1 ; \ 10]/ \ 111 'V

13 ^

^явИШ

16! .

171—| ЁрЩ

18ч ( I

18 Я •

20 4-1--&

123« 1^34 "I г л ч

Рис. 7. Расчетное распределение температуры (°С) в модельной плоскости.

Для всех фаз и зависимых компонентов были построены пространственные и Р-Т-диаграммы содержаний. Практически все диаграммы показывают резкую смену минеральных парагенезисов, составов раствора и газов на геохимическом барьере ПХП. Области наибольшего содержания свободного флюида (рис. 9) локализуются преимущественно в нижней, пластичной зоне и в зоне магматического очага в двух первых сценариях - при его внедрении в хрупкую среду, что свидетельствует о присутствии некоторого количества расплава. Максимальные содержания газовой фазы флюида и области концентраций зависимых компонентов водного раствора флюида приурочены к зонам

магматического очага.

Сценарий I

Сценарий II

Сценарий III

300

\ \

N V

у 1 (с

"8. С

79 22С0 1011

131

154

3400

4000

17-Е

4800-

18^:.

3400

4000

4600

13-1

19-В

ЧЛИ1 -

20?

5260 -

с;

300-

600

2?0Й

Рис. Я. Расчетное распределение полей давления (бар) в модельной плоскости.

В хрупкой зоне изобары проведены через 100 бар, в пластичной зоне - через 200 бар. Распределение компонентов твердых фаз контролируется геохимическим барьером ПХП и Р-Т-условиями подсистем. Минералы кремнезема в системе представлены а-кварцем (рис. 10), тридимитом, кристобалитом, халцедоном и аморфным кремнеземом, причем в нижней, пластичной зоне равновесен только аморфный кремнезем в небольших количествах. Остальные минералы кремнезема равновесны выше геохимического барьера ПХП. Полевые шпаты представлены, в основном, плагиоклазами и микроклином. По равновесным содержаниям альбита и анортита можно судить о номере плагиоклазов в той или иной подсистеме. Альбит характерен для хрупкой зоны, анортит же приурочен к подсистемам магматического очага. Максимальные содержания альбита - до 44 масс. %, анортита-до 28 масс. %. Микроклин существует на глубинах 13-18 км (до 12 масс. %), и точечно - в хрупкой зоне (до 8 масс. %). В периферийных частях магматического очага равновесен гематит (до 8 масс. %), в меньших количествах он присутствует в зоне хрупких деформаций (до 2.5 масс. %).

Сценарий II

Сиенаоий III

1 2 3

Рис. 9. Расчетные диаграммы распределения свободного флюида, масс. %. Сиенавдй I Сиенаоий II С цен алий III

Сценарий I

Рис. 10. Расчетные диаграммы распределения а-кварца, масс. %

Пироксены представлены эгирином, геденбергитом, диопсидом. Эгирин равновесен в пластичной зоне (до 14 масс. %), геденбергит - на средних глубинах (до 5.4 масс. %), причем их изоконцентраты огибают магматический очаг снаружи, не пересекая его границы. Среди амфиболов отметим антофиллит, о составе которого можно судить по содержаниям крайних членов изоморфного ряда: Ре-антофиллит - до 4 масс. % и Г^-антофиллит - до 15 масс. %. Флюорит разных модификаций характерен практически для всей системы; его распределение сложно интерпретировать. В первом сценарии он частично отсутствует непосредственно под геохимическим барьером ПХП и на глубинах 16-17 км, во втором - на глубинах 9-17 км, в третьем - в области между магматическим очагом и линией ПХП. Максимальные содержания флюорита -до 0.12 масс. %. В средне- и низкотемпературных областях зоны хрупких деформаций равновесны некоторые цеолиты, например, пумпеллиит. Из рудных минералов отметим халькозин, существование которого строго контролируется положением магматического очага (рис. 11).

Сценарий I

Сценарий II

Сценарий Ш

оою

0.003 0.008 0.007 О 006 0 005 0.004 0.003 0.002 0.001 0 000

Рис. 11. Расчетные диаграммы распределения халькозина, масс. %.

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

На динамическом этапе использовался метод последовательного проточного реактора для группы подвижных фаз (раствор + газ) с 12 резервуарами в каждом сценарии, расположенными на траектории движения флюида в конвективной ячее МГС зоны хрупких деформаций. Модель усложнялась привносом в первый резервуар глубинного флюида, подъем которого через пластичную среду по

предложенному механизму был приповерхностные резервуары поверхность (рис. 12).

Сценарий 1

рассчитан ранее (глава 3), а также привносом в метеорного раствора и выносом флюида на

1 _3_ 41

_L 0

2 Я

т ш Я

9 31 ш

5 я

6 ¡1 I Я

7 ш 2; SS

8 Т

Ш _

11 .1

14 1

м И

L "L

20 і

1 2 3 4

1 е «:

2 ? 6

8 5

4 9 4

5 .10 3

6 11 г

7 12 S

8 і й

9 1

13 14

1 - 3 4

1 0

2 Р

81 6 "J

5 1(1 4

6 11 3

7 12 2

8 1

12 13

Рис. 12. Схемы сценариев динамической модели. Цифрами в сетке даны номера резервуаров. Нулевые значения -внешняя среда. Глубинный флюид (О внизу) попадает в первый резервуар (1), дальше развивается сценарий проточного реактора, то есть перемещение группы подвижных фаз из резервуара в резервуар (1-12). Массообмен с поверхностью происходит в первом и третьем сценариях через седьмой резервуар (внешняя среда 0), во втором сценарии - через шестой и седьмой резервуары. Конвективные ячеи сценариев расположены в хрупкой зоне, в первом и втором сценариях они захватывают трещиноватые части магматических очагов.

Равновесия рассчитывались методом минимизации энергии Гиббса. Расчет производился до установления стационарного распределения вещества в конвективной ячее МГС (прекращения изменения составов резервуаров) при фиксированных Р-Т-условиях. После установления стационарного распределения вещества в конвективной ячее МГС расчет показал возможность равновесного существования в резервуарах первого сценария 48 минералов, 10 компонентов газовой фазы и 68 компонентов водного раствора. В резервуарах второго сценария расчет выявил 36 минералов, 5 компонентов газовой фазы и 71 компонент водного раствора. В резервуарах третьего сценария определились 27 минералов и 63 компонента водного раствора (газовая фаза отсутствует). Для компонентов были построены диаграммы изменения содержаний в зависимости от глубины и Р-Т-условий. Поведение флюида в резервуарах опишем в виде диаграмм изменения содержания в нем независимых компонентов для первого и двенадцатого резервуаров каждого сценария относительно исходного флюида. Для первого и второго сценариев - это диаграммы изменения составов фазы водного раствора и газовой фазы, для третьего - только фазы водного раствора, так как газовая фаза в его резервуарах отсутствует. Последний результат можно объяснить относительно малыми температурами в зоне хрупких деформаций третьего сценария (66-291°С) при давлениях 195-724 бар. Заметим, что фазовый анализ (раствор/газ) в нашем случае очень сложен, кривая кипения флюида переменного в зависимости от Р-Т—условий состава может быть проверена только экспериментально, что выходит за рамки нашей работы. На рис. 13, 14 и

15 показаны соответственно диаграммы изменения составов водного раствора для 1-3 сценариев.

Н С О Р №М А1 Э С1 КСаРеСигпБпРЬ

1.0Е+18

| 1.0Е+14

| 1,0Е+10

х 1.0Е+06

ф 1.0Е+02

х 1.0Е-02

(D X 03 S

1.0Е-06 1.0Е-10

Глубинный флюид ■ Вэзервуар 1 Резервуар 12

Рис. 13. Изменение концентрации независимых компонентов фазы водного раствора в первом и последнем резервуарах конвективной ячеи МГС по первому сценарию относительно глубинного флюида, мг/кг Н20.

Рис. 14. Изменение концентрации независимых компонентов фазы водного раствора в последнем резервуарах конвективной ячеи МГС по второму сценарию относительно глубинного флюида, мг/кг Н20.

Н СО F Na М Al Si S CI Са Fe Си Zn Sn Pb

1 1.0E+16

g_ 1.0ЕИ2

| 1.0E+08

x 1,0E+04 о

ф 1.0E+00

§ 1.0Б04

ф 1.0E-08

3 1.0E-12

—«— Глубинным флюид —Й— Рэзервуар 1 —Рэзервуар 12

Рис. 15. Изменение концентрации независимых компонентов фазы водного раствора в первом и последнем резервуарах конвективной ячеи МГС по третьему сценарию относительно глубинного флюида, мг/кг НгО.

Н С О F Na М Al Si S CI К CaFe CuZnSnPb

I 1.0E+18

g. 1.0E+14 S 1.0E+10

1,0E+O6 1.0E+02 1.0E-02 1.0E-06 1.0E-10

Глубинный флюид —в— Резервуар 2 —л— Рэзервуар 12

первом и

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

Рис. 16. Диаграмма распределения а-кварца, масс. %

3,00

2,50

0,50

0,00

1 2 3 4 5 6 7 8 9 10 11 12 Резервуары

—Ф— Сценарий 3 А Сценарий 2 —■—Сценарий 1

•Сценарий 3 -Сценарий 2 -Сценарий 1

1 2 3 4 5 6 7 8 9 10 11 12 Резервуары

Рис. 17. Диаграмма распределения 1У^-хлорита М&81205(0Н)4, масс. %.

Разработанная комплексная компьютерная модель МГС вблизи зоны (ПХП) с учетом привноса глубинного флюида и взаимодействия с метеорным раствором позволила выявить основные закономерности существования МГС в данной реологической обстановке, определить роль ПХП и магматического очага в эволюции вещественного состава системы, локализовать области концентрации

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

1. Моделирование показало, что активно участвуют в самоизоляции МГС коллоидные соединения кремнекислоты, высаживающийся на ПХП силикагель и образующиеся из него минералы кремнезема - опал, тридимит, кристобалит, халцедон и кварц. Такой механизм взаимодействия «вода-порода» реализуется в геотермальных процессах в недрах Северо-Парамуширской магматогенно-гидротермальной системы (Рычагов и др., 2002). Как показало бурение скважины ГП-3 до глубины 2.5 км, над интрузивным телом габбро-диоритов образуется зона кипения металлоносных гидротерм, которое приводит к осаждению в пористых вулканогенно-осадочных породах силикатного геля с последующей раскристаллизацией его в минералы кремнезема.

2. Натурные и теоретические исследования (Ващилов, 1984; Corbett, Leach, 1998; Рычагов, Белоусов, 2002; Рычагов, 2003 и многие другие) показывают хорошее соответствие природной фактуры с расчетными данными. Несмотря на то, что внедрение полифазной интрузии происходит в пределах верхней части земной коры, здесь за счет длительной эволюции сложного магматического тела реализуются механизмы сценариев 1 и 2 и, возможно, частично сценария 3.

3. Расчетной зоне ПХП от остывающего магматического тела к вмещающим горным породам соответствует и так называемая «брекчиевая мантия» (Рычагов, 1989) - комплекс различного типа брекчиевых образований (от интрузивных и автомагматических до гидротермальных и метасоматических брекчий) в эндо- и экзоконтактовых частях интрузий.

ГЛАВА 5. ТЕКТОНОФИЗИЧЕСКИЕ ИССЛЕДОВАНИЯ ЭВОЛЮЦИИ ФЛЮИДА Цель главы - физическая проверка состоятельности гипотезы движения флюидозаполненных полостей в пластичной среде.

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

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

5.3. Основные выводы. Эксперименты, проведенные в 2006-2011 гг. позволили уверенно обосновать гипотезу миграции флюидозаполненных полостей. В

гравитационном поле в пластической вмещающей среде происходит миграция флюидозаполненных полостей по механизму гидроразрыва, скорость которой зависит от величины избыточного давления. Для начала миграции полостей необходима критическая высота полости, обеспечивающая критическое ИДФ, равное прочности породы. Траектория движения зависит от напряженного состояния среды, в которой процесс развивается. Полость при движении стремится в зону растягивающих напряжений. При наличии барьеров, прочность на разрыв которых превышает величину ИДФ в полости, последняя некоторое время распространяется по латерали до момента понижения избыточного давления до уровня прочности вмещающей среды. На основе анализа более семидесяти экспериментов было рассчитано среднее критическое ИДФ, необходимое для гидроразрыва модельной среды: 375.39 Па. Средняя критическая длина полости в модели составляет 6.01 см, что при природной прочности пород 130 кПа соответствует -42.84 м, а при прочности 3 МПа соответствует -988.59 м.

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

ЗАКЛЮЧЕНИЕ По существу проведенных исследований можно сделать выводы:

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

2. . Предложено три сценария (модели) формирования циркуляционных МГС, в

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

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

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

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

5. На основании сравнения полученных расчетных данных с данными комплексных геолого-геофизических, гидрогеологических, минералого-геохимических и других исследований МГС Южной Камчатки и острова Парамушир показано высокое соответствие модели внедрения магматических расплавов в хрупкую земную кору строению современных высокотемпературных МГС (прогрессивного этапа развития) на глубинах более 1.0-1.5 км, и рудоносных вулканоплутонических комплексов миоцен-плиоценового возраста, эродированных на различную глубину в разных геологических блоках. В структуре вулкано-плутонических комплексов зоной хрупко-пластического перехода служит брекчиевая мантия (Рычагов, 1989) интрузивных габбро-диоритовых тел, контролирующая потоки высокотемпературных металлоносных газо-водных флюидов.

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в рецензируемых журналах:

1. Васильев В.И., Чудненко К.В., Жатнуев Н.С., Васильева КВ. Комплексное компьютерное моделирование геологических объектов на примере разреза зоны субдукции // Геоинформатика, №3,2009. - С. 15-30.

2. Васильев В.И., Жатнуев Н.С., Рычагов С.Н., Васильева Е.В., Санжиев Г.Д. Массоперенос и минералообразование в магмагогенно-гидрогермальных системах по результатам численного физико-химического моделирования // Литосфера, №3,2010. - С. 145-152.

3. Жатнуев Н.С., Рычагов С.Н., Васильев В.И., Васильева КВ. Влияние реологических свойств земной коры на локализацию рудообразуюших гидротермально-магматических систем // Вулканология и сейсмология (принята в печать 19.09.2011).

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

4. Васильева КВ., Васильев В.И., Жатнуев Н.С. Тектонофизическое моделирование динамики флюидосодержащих трещин в литосфере // Материалы Всероссийской конференции с иностранным участием «Геохимия и рудообразование радиоактивных, благородных и редких металлов в эндогенных и экзогенных процессах», посвященной 50-летию СО РАН и 80-легаю чл.-корр. РАН Ф.П. Кренделева, г. Улан-Удэ: Изд-во БНЦ СО РАН, 2007. - С.122-125.

5. Васильева Е.В., Васильев В.И. Эволюция флюидозаполненных трещин в литосфере по результатам тектонофизического моделирования // Современные проблемы геологии, геохимии и геоэкологии Дальнего Востока России. Материалы II региональной конференции. -Владивосток: Дальнаука, 2008. - С. 111-113.

6. Васильева ЕВ., Васильев В.И., Жатиуев НС. Миграция флюидозаполненных трещим и литосфере по результатам тектонофизического моделирования // Геодинамическая эволюция литосферы Центрально-Азиатского подвижного пояса (от океана к континенту): Материалы всероссийского совещания. Вып. 6. - Иркутск: ИЗК СО РАН, 2008.Т. 1. - С. 54-55.

I. Жатиуев Н.С., Тубанов Ц.А., Васильева Е.В. Упруго-пластичный переход как концентратор напряжений в земной коре // Разломообразование и сейсмичность в лиггосферс: тектонофизические концепции и следствия: Материалы всероссийского совещания. - Иркутск: ИЗК СО РАН, 2009. - Т. 2. - С. 22-23.

8. Васильев В.И, Жатиуев Н.С., Васильева Е.В. Численная термодинамическая модель подъема флюидозаполненных трещин через зону пластично-хрупкого перехода // Геодинамическая эволюция литосферы Центрально-Азиатского подвижного пояса (от океана к континенту): Материалы всероссийского совещания. Вып. 7. -Иркутск: ИЗК СО РАН, 2009. - Т.1. -С. 48-50.

9. Жатнуев Н.С., Васильева ЕВ. Трещинно-флюидный механизм мантийных плюмов // Геодинамическая эволюция литосферы Центрально-Азиатского подвижного пояса (от океана к континенту): Материалы всероссийского совещания. Вып. 7. - Иркутск: ИЗК СО РАН, 2009. -Т.1.-С. 104-106.

10. Васильева ЕВ., Васильев В.И., Санжиев Г.Д. Процесс миграции флюидозаполненных трещин в литосфере по результатам тектонофизического моделирования // Геодинамическая эволюция литосферы Центрально-Азиатского подвижного пояса (от океана к континенту): Материалы всероссийского совещания. Вып. 8. - Иркутск: ИЗК СО РАН, 2010. - Т. 1. - С. 46-47.

II. Васильев В.И., Васильева ЕВ. Использование комплексной методики для компьютерного моделирования распределения тепла и вещества в зонах субдукции // Геохимия, петрология и рудоносность базит-ультрабазитовых комплексов: Материалы II всероссийской молодежно!! школы-семинара. - Иркутск: ИГХ СО РАН, 2010. - С. 19-21.

12. Жатнуев Н.С., Васильева Е.В. Трещинно-флюидный механизм мантийных плюмов // Геохимия, петрология и рудоносность базит-ультрабазитовых комплексов: Материалы H всероссийской молодежной школы-семинара,-Иркутск: ИГХ СО РАН, 2010. - С. 14-16.

13. Васильева Е.В., Жатнуев Н.С., Васильев В.И., Санжиев Г.Д. Возможность участия мантийного флюида в формировании гидротерм Байкальской рифтовой зоны по результатам тектонофизического и численного физико-химического моделирования И Геология Западного Забайкалья: Материалы всероссийской молодежной конференции. - Улан-Удэ: Изд-во БГУ, 2011.-С. 36-38.

14. Васильев В.И., Борхонова Е.В., Чернявский М.К., Васильева Е.В. Численная физико-химическая модель формирования термальных вод Горячинского месторождения // Геология Западного Забайкалья: Материалы всероссийской молодежной научной конференции. - Улан-Удэ: Изд-во БГУ, 2011.-С. 26-32.

15. Васильев В.И., Васильева Е.В. Прогноз антропогенного изменения состава термальных вод Горячинского месторождения II Геология Западного Забайкалья: Материалы всероссийской молодежной научной конференции. - Улан-Уда: Изд-во БГУ, 2011. - С. 32-35.

16. Vasiliev V.l., Zhatnuev N.S., Vasilieva Е. V., Sanzhiev G.D. The physical and chemical computation of the conditions of oceanic plate dehydration under the subduction at rates of 1.0-5.0 cm/year // Large Igneous Provinces of Asia, Mantle Plumes and Metallogeny: Proceedings of International Symposium. - Irkutsk: IOC SB RAS, 2011. - P. 284-287.

17. Vasiliev V.l., Zhatnuev N.S., Vasilieva E. V., Sanzhiev G.D. The viscosity of melt in magma pocket at the brittle-plastic border in the lithosphere // Large Igneous Provinces of Asia, Mantle Plumes and Metallogeny: Proceedings of International Symposium. - Irkutsk: ЮС SB RAS, 2011. - P. 288-291.

18. Васильева E.B., Санжиев Г.Д., Васильев В.И. Поведение корового флюида на геохимическом барьере пластично-хрупкого перехода литосферы по результатам численного физико-химического моделирования Н Современные проблемы геохимии: Материалы Всероссийской молодежной конференции. - Иркутск: ИГХ СО РАН, 2011. -С. 217-221.

19. Васильев В.И., Жатнуев Н.С., Васильева Е.В., Санжиев Г.Д. Численная физико-химическая модель эволюции рудоносного флюида и вмещающих граиитоидов на геохимическом барьере пластично-хрупкого перехода земной коры II Граниты и процессы рудообразования: Материалы Всероссийской конференции, посвященной памяти академика В.И. Коваленко. - Москва: ИГЕМ РАН,2011.-С. 39-41.

Подписано в печать 31.01.2012 г. Формат 60x84 1/16. Бумага офсетная. Объем 1,4 печ. л. Тираж 100. Заказ № 22.

Отпечатано в типографии Изд-ва БНЦ СО РАН. 670047 г. Улан-Удэ ул. Сахьяновой, 6.

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

61 12-4/78

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ НАУКИ ГЕОЛОГИЧЕСКИЙ ИНСТИТУТ

СИБИРСКОГО ОТДЕЛЕНИЯ РАН

ГЕОХИМИЯ МАГМАТОГЕННО-ГИДРОТЕРМАЛЬНЫХ СИСТЕМ НА ГРАНИЦЕ ПЛАСТИЧНО-ХРУПКОГО ПЕРЕХОДА

В ЗЕМНОЙ КОРЕ (ФИЗИКО-ХИМИЧЕСКИЕ МОДЕЛИ)

Специальность: 25.00.09 - «Геохимия, геохимические методы поиска полезных ископаемых»

Диссертация на соискание ученой степени кандидата геолого-минералогических наук

УДК 550.34.013.4:550.424.6+550.347.097.35

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

Васильева Евгения Владимировна

Научный руководитель: д.г.-м.н. Н.С. Жатнуев

Улан-Удэ

-2012

ОГЛАВЛЕНИЕ

Введение 3

Глава 1 .Объект и методы исследования 11

1.1. Магматогенно-гидротермальные системы (МГС) -

1.2. Численное физико-химическое моделирование 18

1.3. Методика комплексного компьютерного моделирования 23

1.3.1. Геометрический этап 24

1.3.2. Физический этап 27

1.3.3. Физико-химический этап 33

1.3.4. Динамический этап 40

1.4. Тектонофизические исследования 41

1.4.1. Поляризационно-оптический метод исследования напряжений 44

1.4.2. Метод анализа результатов исследований 48 Глава 2. Концептуальная модель эволюции магматогенно-гидротермальных систем 50

2.1. Реологическая модель -

2.2. Механизм движения полостей в пластичной среде 56

2.3. Сценарии эволюции магматогенно-гидротермальной системы 62 Глава 3. Модель подъема флюида через пластично-хрупкий переход 71

3.1. Описание модели -

3.2. Результаты моделирования 82

3.3. Основные выводы 86 Глава 4. Комплексная компьютерная модель магматогенно-гидротермальной системы 88

4.1. Геометрический этап -

4.2. Физический этап 89

4.3. Физико-химический этап 94

4.4. Результаты статического моделирования 106

4.5. Динамический этап 117

4.6. Результаты динамического этапа 120

4.7. Основные выводы 128 Глава 5. Тектонофизические исследования эволюции флюида 130

5.1. Подготовка экспериментов -

5.2. Описание экспериментов 140

5.3. Основные выводы 147 Заключение 148 Литература 151

ВВЕДЕНИЕ

Актуальность работы

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

В настоящее время все еще преобладают взгляды на формирование МГС в виде чисто конвективного тепломассопереноса в трещиновато-пористой среде, контролируемого теплом магматического очага. Такой подход не учитывает реологические свойства вмещающей среды и связанные с ними геохимические барьеры, а также оставляет «за кадром» химическую термодинамику взаимодействия флюид/порода как в конвективной ячее, так и в процессе привноса глубинного и метеорного вещества. Поэтому на современном этапе геохимических исследований чрезвычайно важно рассмотреть возникновение и эволюцию МГС именно с реологических позиций, позволяющих задать обоснованные входные данные для расчета химико-термодинамического взаимодействия сред в различных реологических об-становках.

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

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

Цель работы

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

Задачи работы

были поставлены в соответствии с выбранной целью:

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

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

а) для траектории подъема глубинного флюида в пластической области с прорывом в область хрупких деформаций;

Фактический материал

Фактические данные по химическому составу природных сред и базы термодинамических данных были любезно предоставлены непосредственно Н.С. Жатнуевым, В.И. Васильевым (Геологический институт СО РАН, Улан-Удэ), К.В. Чудненко, В.А. Бычинским (Институт геохимии им. А.П. Виноградова СО РАН, Иркутск), С.Н. Рычаговым (Институт вулканологии ДВО РАН, Петропавловск-Камчатский) и другими учеными.

В работе также использовался фактический материал, самостоятельно полученный автором в результате тектонофизических экспериментов на оптически активных средах в поляризованном свете, проводимых при финансовой поддержке гранта НШ-2284.2003.5 и Президиумов СО и ДВО РАН (проект № 117-09-П-СО-08-006).

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

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

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

2. На основе этой концептуальной модели разработана комплексная компьютерная модель физико-химической эволюции МГС в зоне изменения реологических свойств пород с привносом глубинного и метеорного вещества.

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

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

1. Разработан 20-вариант многорезервуарного физико-химического моделирования на программном комплексе «Селектор», что позволяет получать более корректные модели геологических систем с участием твердых и флюидных фаз.

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

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

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

Автором разработаны концептуальные модели эволюции МГС, в чем оказали помощь д.г.-м.н. Н.С. Жатнуев (ГИН СО РАН) и д.г.-м.н. С.Н. Рычагов (ИВиС ДВО РАН). Автором реализованы геометрический, физический и динамический этапы моделирования. Также автором были выполнены все физико-химические расчеты на ПК «Селектор» с применением сторонних программных продуктов (Васильев, ГИН СО РАН; Арискин, ГЕОХИ РАН; Персиков, ИЭМ РАН). Автором также лично разработана методика тектоно-физических экспериментов, изготовлены модельные установки, проведены эксперименты и анализ полученного материала, выполнена интерпретация результатов.

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

1. Возникновение и эволюция магматогенно-гидротермальных систем возможны по трем сценариям:

б) при локализации очага в области пластично-хрупкого перехода имеет место быстрое остывание апикальной части массива и частичное смешение ювенильного и метеорного материала;

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

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

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

Публикации и апробация

Основные положения диссертации опубликованы в 22 печатных работах (в том числе 2 статьи в рецензируемых журналах «Геоинформатика» и «Литосфера» и 1 статья принята в печать в рецензируемый журнал «Вулканология и сейсмология»), а также докладывались на Всероссийской конференции, посвященной 80-летию чл.-корр. РАН Ф.П. Кренделева (Улан-Удэ, 2007), Региональной молодежной конференции, посвященной 50-летию Бурятского научного центра (Улан-Удэ, 2008), Региональной конференции молодых ученых «Современные проблемы геологии, геохимии, геоэкологии Дальнего Востока России» (Владивосток, 2008), Всероссийском научном совещании «Геодинамическая эволюция литосферы Центрально-Азиатского подвижного пояса: от океана к континенту» (Иркутск, 2008, 2009, 2010), Все-

российском совещании «Разломообразование и сейсмичность в литосфере: тектонофизические концепции и следствия» (Иркутск, 2009), Всероссийской молодежной конференции «Структура, вещество и история литосферы» (Сыктывкар, 2009), Международной конференции ЮНЕСКО «Глобальные и региональные проблемы устойчивого развития мира» (Улан-Удэ, 2010), Всероссийской молодежной конференции «Геология Западного Забайкалья» (Улан-Удэ, 2011), Международного симпозиума «Large Igneous Provinces of Asia, Mantle Plumes and Metallogeny» (Иркутск, 2011), Всероссийской молодежной конференции «Современные проблемы геохимии» (Иркутск, 2011), Всероссийской конференции «Граниты и процессы рудообразования» (Москва, 2011), научных сессиях и семинарах Геологического института СО РАН, Института геохимии им. А.П. Виноградова СО РАН, Бурятского научного центра и Бурятского Государственного университета.

Объем и структура работы

Диссертация состоит из введения, пяти глав, заключения и списка литературы общим объемом 162 страниц печатного текста, 13 таблиц и 106 иллюстраций.

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

Во второй главе предлагается концептуальная модель магматогенно-

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

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

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

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

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

Список литературы содержит 119 наименований.

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

Автор считает своим долгом выразить сердечную благодарность в пер-

вую очередь своему научному руководителю - д.г.-м.н. Николаю Сергеевичу Жатнуеву за терпеливую плодотворную помощь на протяжении всего времени учебы и работы над диссертацией. Автор душевно благодарит своих учителей: д.г.-м.н. Д.И. Царева, д.г.-м.н. А.Г. Миронова, д.г.-м.н. A.M. Плюсни-на, к.г.-м.н. В.И. Васильева, к.г.-м.н. Д.А. Орсоева, к.г.-м.н. Г.С. Риппа, к.г.-м.н. Т.Т. Врублевскую и других преподавателей. В подготовке и проведении экспериментов непосредственную помощь оказали Г.Д. Санжиев, М.В. Бурцева и к.г.-м.н. Б.Б. Дамдинов, за что им автор очень благодарна. Особую признательность хотелось бы выразить сотрудникам Лаборатории физико-химического моделирования Института геохимии им. А.П. Виноградова СО РАН д.г.-м.н. К.В. Чудненко и к.г.-м.н. В.А. Бычинскому за предоставленные версии ПК «Селектор» и специальные базы термодинамических данных.

ГЛАВА 1. ОБЪЕКТ И МЕТОДЫ ИССЛЕДОВАНИЯ

1.1. Магматогенно-гидротермальные системы

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

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

Достаточно полный обзор литературы по гидротермальным системам приводится в работах В.Д. Пампуры, который, следуя за Дж. Элдером, рассматривает их как «гидрогеологические резервуары напорных газоводных растворов, возникшие в областях аномального теплового потока при смешении метеорных вод с высокоэнтальпийными флюидами и потоками магматических газов над периферическими магматическими очагами или зонами локального плавления корового субстрата» (Пампура, 1981, 1985).

Г.Л. Поспелов понимал под гидротермальной системой «фильтрующуюся термогидроколонну», при этом имелась ввиду «газогидродинамическая система», включающая только жидкое и газообразное вещество, циркулирующее во вмещающих породах (Поспелов, 1962).

В.И. Белоусов, В.М. Сугробов, В.В. Аверьев и С.И. Набоко считают, что «в областях тектоно-магматической активности длительное время существует аномальный тепломассопоток, реализующийся в благоприятных гидрогеологических условиях в виде гидротермальных систем с областями разг-

рузки термальных вод на поверхности Земли» (Набоко, 1980).

Н.С. Жатнуев, С.Н. Рычагов и И.Ф. Делемень рассматривают гидротермальную систему в виде закономерно органи�

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

Васильева, Евгения Владимировна
кандидата геолого-минералогических наук
Улан-Уде, 2012
ВАК 25.00.09

Диссертация

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

Автореферат

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