В отечественной геологической науке для решения некоторых задач структурной геологии, реконструкции напряжений и деформаций горных пород, изучения структур рудных полей и месторождений была разработана и применялась методика микроструктурного анализа. Данный анализ впервые возник в 1930 г., после того как Б. Зандер сформулировал его основные положения. Позднее микроструктурные исследования проводились многими, преимущественно советскими, учёными и исследователями. Однако в настоящее время развитие и повсеместное использование при геологических и структурно-геологических исследованиях рассматриваемого анализа, из-за крайней сложности его проведения с технической точки зрения, необходимости сбора и обработки большого объёма информации, а также отсутствия иных независимых методов верификации полученных результатов, практически прекратилось.
В рамках современного этапа развития геологической науки для решения структурных задач всё чаще делается акцент на изучении трещин и микротрещин в ориентированных образцах горных пород, измерении их геометрических параметров и анализе минерального выполнения. Полученные результаты позволяют делать выводы о последовательности образования структур, выявлять этапы тектонических движений, восстанавливать пути миграции флюидов в трещинном пространстве горной породы, реконструировать параметры поля напряжений-деформаций (ПНД), связанные с конкретным тектоническим этапом. На основе накопленного опыта комплексного изучения микроструктур в образцах горных пород авторами создана специальная методика микроструктурного анализа (СММА) и сформулированы теоретические основы её применения.
Цель исследования: разработка и обобщение теоретических основ созданной авторами СММА, а также демонстрация возможностей и результатов её применения на основе использования дополнительных прецизионных методов анализа минерального вещества.
Материалы и методы исследования
Среди микроструктур в горных породах обычно выделяют открытые микротрещины, которые не заполнены минеральным веществом, а также минерализованные микротрещины, которые выполнены либо вторичными минералами, либо вторичными флюидными включениями. Захват вторичных флюидных включений происходит под действием какого-либо позднего (например, гидротермального) процесса после того, как кристаллизация вмещающего их минерала была уже завершена, при наличии в нём доступного проницаемого трещинно-порового пространства [1].
Минерализованные и частично минерализованные трещины являются важными индикаторами протекания процессов миграции флюидов в прошлые геологические эпохи. Вторичные флюидные включения обычно «захватываются» микротрещинами, которые сформировались на фоне ориентированного стресса и потому представляют собой ориентированные закономерным образом генерации. Данные микротрещины при проведении структурных исследований получили название «планарные системы флюидных включений» (ПСФВ) [2].
Пространственная ориентировка каждой генерации ПСФВ также напрямую зависит от параметров действовавшего на момент её формирования ПНД. Поэтому исследования подобных планарных систем в качестве структурных маркеров в сочетании с детальным комплексным изучением других типов систем микроструктур можно использовать для восстановления этапов деформаций и проницаемости пород, реконструкции геометрии путей миграции флюидов, установления динамики изменения термобарических и физико-химических условий, связанных с определёнными этапами деформации геологических тел.
ПСФВ на начальных этапах деформаций формируются как микротрещины отрыва. Плоскость трещины чаще всего ориентирована перпендикулярно оси наименьшего сжатия σ3. При этом вектор максимальной проницаемости лежит в плоскости промежуточной оси σ2 и оси максимального сжатия σ1, где обеспечивается наиболее интенсивная миграция флюидов (рис. 1).
Рис. 1. Ориентировка осей поля напряжений-деформаций в зависимости от ориентировки ПСФВ; а – в трехмерном пространстве; б – в плоскости ориентированного шлифа: 1 – ПСФВ, 2 – открытые или минерализованные микротрещины, 3 – границы минеральных зёрен, Q – зерна кварца. σ1 – ось максимального сжатия, σ1 – ось наименьшего сжатия, σ1σ2 – плоскость промежуточной оси и оси максимального сжатия
Через некоторый промежуток времени смена тектонической обстановки приводит к компрессии, под действием которой активные микроструктуры запечатываются, «захватывая» первые флюидные (газово-жидкие) включения. Следующая тектоническая обстановка сопровождается новым этапом деформаций, следы которого фиксируются в ориентировке ПСФВ второй генерации, а состав и свойства выполняющих их включений характеризуют термобарические, физические и химические условия. Данный процесс повторяется с образованием новых генераций ПСФВ, до окончательного затухания гидротермальной активности.
Как показали в своих работах Смит и Эванс (1984), ПСФВ могут быть обнаружены во многих породообразующих минералах горных пород, вне зависимости от их кристаллографических свойств. При этом лучшую сохранность обеспечивают зерна кварца. Опыты продемонстрировали, что кристаллы полевых шпатов и карбонатов менее устойчивы к гидротермально-метасоматическим преобразованиям, процессам выщелачивания и выветривания, соответственно, сохранность ПСФВ и возможности их изучения в данных минералах малы.
Именно на установление ориентировок генераций ПСФВ в сочетании с методами термобарогеохимии были направлены усилия зарубежных ученых при изучении рудных объектов. При этом другие важные типы микротрещин оставались без внимания. Необходимо отметить, что изучение открытых микроструктур и микротрещин, в которых происходило отложение минерального вещества, не всегда позволяет реконструировать параметры ПНД. Но измерение пространственно-геометрических параметров данных типов микротрещин и изучение характера их минерального выполнения играет решающую роль при конечном расчёте фильтрационных параметров горных пород, так как они, в отличие от ПСФВ, характеризуются большей протяжённостью и апертурой (шириной раскрытия), способствуя увеличению фильтрационной способности горных пород (рис. 1). Кроме того, параметры данных микроструктур необходимо учитывать при реконструкции геометрии трещинно-порового пространства, восстановления направлений и условий движения палеопотоков флюидов. Поэтому авторский подход к реализации СММА подразумевает, помимо установления ориентировок ПСФВ, решение задачи по реконструкции пространственных параметров и характера минерального выполнения минерализованных и открытых микроструктур.
Важным условием для реализации СММА, как и для классической методики, является отбор ориентированных образцов. Далее из образцов обычно изготавливаются ориентированные шлифы. Пространственные параметры различных генераций микроструктур и характер их взаимоотношений могут быть изучены под оптическим микроскопом при помощи столика Федорова. Но при таком подходе измерение геометрических характеристик микроструктур, количество которых может достигать нескольких тысяч в одном ориентированном шлифе, занимает значительное время. Всё это определяет необходимость создания специализированного инструмента, который мог бы ускорить сбор, обработку и анализ пространственных микроструктурных данных. Так авторами был разработан ГИС модуль для реализации СММА, изучения структурных неоднородностей горных пород, установления геометрии трещинно-порового пространства, различных типов и систем микроструктур, их пространственных взаимоотношений, реконструкции последовательности активизации и расчета фильтрационных характеристик [3].
Результаты исследования и их обсуждение
Реконструкция ориентировки осей поля напряжений-деформаций. Благодаря применению СММА, реализованной с помощью авторского программного обеспечения, устанавливаются геометрические параметры микроструктур и их генераций, характеризующих параметры трещинно-порового пространства определённого этапа деформаций.
Результаты выявления ориентировок ПСФВ, с учетом того, что они в большинстве случаев являются микротрещинами отрыва (или «гибридными») и их плоскости в большинстве случаев ориентированы преимущественно перпендикулярно оси наименьшего сжатия, могут быть восстановлены параметры ПНД на конкретном этапе эволюции гидротермальной системы.
В качестве примера можно привести работу авторов по изучению ориентировок ПСФВ с помощью СММА на урановом месторождении Антей (юго-восточное Забайкалье). На 9 и 11 горизонтах месторождения были выявлены две наиболее отчётливо проявленные генерации ПСФВ, имеющие северо-восток-юго-западные и северо-запад–юго-восточные простирания. На основе полученных результатов были построены сводные розы-диаграммы ориентировок ПСФВ (рис. 2) с учетом суммарной длины каждой генерации. Ориентировки сопоставлялись с результатами реконструкции динамики изменения ПНД, выполненной ранее на основе тектонофизических методов изучения макроструктур [4]. Оказалось, что ориентировка выявленных генераций ПСФВ совпадает с восстановленным простиранием осей максимального сжатия σ1, характеризующих второй и третий этапы тектогенеза, имевших место в процессе формирования месторождения. Полученные результаты подтвердили возможность использования ПСФВ в качестве структурных маркеров для реконструкции параметров ПНД, действовавших в палеоусловиях.
Рис. 2. Ориентировка ПСФВ на сводных розах-диаграммах для горизонтов месторождения Антей. n – количество замеров геометрических параметров микроструктур. Красные стрелки – ориентировка оси максимального сжатия σ1 на первом – σ1(I), втором – σ1(II) и третьем σ1(III) этапах тектогенеза, выявленных ранее [4] макроструктурными методами
Восстановление физико-химических параметров гидротермального процесса для определённых стадий деформаций горных пород. Важной особенностью представляемой методики является то, что она позволяет устанавливать связь тектонодинамических процессов (формирование микроструктур) и физико-химических параметров.
Мелкая трещиноватость, рассматриваемая как совокупность небольших по размеру трещин, является не только отражением развития деформаций горных пород, но и формирует необходимые условия для успешного протекания процессов минералообразования.
Для решения задачи установления взаимосвязи структурных и физико-химических параметров недостаточно измерений геометрических параметров, помимо этого требуется использование прецизионных методов анализа минерального вещества, дающих информацию о температурах, давлениях, составе, солёности флюида и т.д. Результаты должны быть сопоставимы со структурными данными, полученными с помощью СММА. Решая данную задачу, авторы смогли сопоставить микроструктурные и тектонофизические данные с результатами определения свойств палеофлюидов (вторичных флюидных включений в микроструктурах) с помощью микротермометрии и состава минерального выполнения микротрещин с использованием сканирующей электронной микроскопии (СЭМ). В дальнейшем результаты могут быть значительно дополнены благодаря применению рамановской спектроскопии [5] вторичных флюидных включений.
На основе синтеза результатов выявления ориентировок ПСФВ с помощью СММА и микротермометрических исследований вторичных флюидных включений, приуроченных к микротрещинам, процесс формирования месторождения Антей был разделён на три структурных эпизода. Каждый эпизод характеризовался определённым набором активных микроструктур, способствовавших ускорению процесса миграции флюидов в рамках гидротермальной активности (таблица), а также уникальными физико-химическими параметрами.
Структурные эпизоды преобразования пород на 9 и 11 горизонтах месторождения Антей по данным СММА и микротермометрических исследований вторичных флюидных включений, формирующих ПСФВ
|
Параметр |
Структурный эпизод I |
Структурный эпизод II |
Структурный эпизод III |
|
Ориентировки ПСФВ |
Субширотные |
Северо-западные (I генерация), северо-восточные |
Субширотные, северозападные (II генерация), субмеридиональные |
|
Температуры гомогенизации флюидных включений, 0C |
380–300 |
300–220 |
220–150 |
|
Тип деформаций |
Упруго-пластичные |
Хрупкие |
|
|
Гидротермально-метасоматические преобразования |
Калишпатиты, альбититы |
Грейзенизация |
Аргилизация, хлоритизация, гематитизация |
Реконструкция внутрирудной тектоники и путей миграции рудоносных флюидов. На основе комплексного изучения ориентировок выполненных рудным веществом микроструктур с помощью СММА в сочетании с различными методами анализа минерального вещества (СЭМ, рентгеноспектральный микроанализ, атомно-абсорбционная спектроскопия, радиография и др.) появляется уникальная возможность реконструкции внутрирудной тектоники и путей миграции исключительно рудоносных флюидов.
Так, применение осколковой радиографии и СЭМ в комплексе с СММА позволило авторам для месторождения Антей воссоздать последовательность активизации микроструктур, выполненных урановорудным веществом и характеризующих внутрирудную тектонику. Рудные скопления локализовались преимущественно в двух системах микротрещин. По тектоническим индикаторам (смещение одной системы относительно другой) было установлено, что первой сформировалась северо-запад-юго-восточная генерация, а затем проявились деформации, приведшие к появлению генерации северо-восток-юго-западного простирания (рис. 3, А). С помощью СЭМ были выявлены отличия в минеральном выполнении данных микроструктур, дополнительно характеризующиеся концентрациями Ti в северо-восточных системах микротрещин (рис. 3, Б). Отсюда следует вывод, что отложение руды месторождения Антей происходило в рамках двух стадий микротектонических смещений, отражающих общий неоднородный тектонический характер функционирования рудоносной гидротермальной системы.
Рис. 3. Реконструкция процессов внутрирудной тектоники с применением осколковой радиографии (А) и изучение минерального выполнения рудоносных микроструктур с помощью СЭМ (Б) в сочетании с результатами СММА
Количественный расчет общих и директивных фильтрационных параметров кристаллических пород в палео- и современных условиях. Помимо прочего, на основе реконструкции пространственных параметров микроструктур с помощью СММА, могут быть количественно рассчитаны фильтрационные характеристики горных пород для каждого выявленного ранее эпизода деформаций. С использованием данного подхода появляется возможность количественно рассчитать значения пористости и проницаемости, формируемые каждой активной системой микроструктур (директивные свойства) и оценить её вклад в общие фильтрационные значения.
Для расчета пористости горных пород, формируемой микротрещинами, авторы предполагают, что микроструктуры, возникшие в рамках одного структурного эпизода, взаимосвязаны. С учётом выявленной геометрии сети микротрещин, количественные расчеты производятся по формуле
где Ф – пористость (эффективная), выраженная в %; L – суммарная длина определённой системы трещин; S – площадь поверхности (в данном случае – площадь ориентированного шлифа); e – апертура (средняя ширина раскрытия) системы трещин.
Расчёты микротрещинной проницаемости осуществляются также исходя из предположения о взаимосвязи систем микротрещин одного эпизода деформаций. Помимо этого, для расчета структурной проницаемости необходимо знать угол падения микротрещин, значения которого пока не представляется возможным точно измерить в плоскости ориентированного шлифа. Поэтому авторами делается допущение о том, что микроструктуры имеют вертикальное падение. Такое допущение абсолютно справедливо и подтверждается сравнением расчётных и полученных в лабораторных условиях значений проницаемости горных пород. Учитывая вышесказанное, микротрещинная проницаемость (Kf), выраженная в м2 или в дарси, может быть рассчитана по формуле
где ε – коэффициент неровности стенок трещин, Nf – объёмная плотность трещин системы f, с – среднее значение половины длины трещин, е – апертура.
Объёмная плотность вычисляется по формуле
где θ – угол падения каждой трещины. Так как трещины перпендикулярны плоскости шлифа – sinθ = 1. Формула для Kf может быть упрощена:
При расчётах трещинной пористости и проницаемости, помимо предложенных формул, могут быть использованы различные подходы и формулы, выведенные на основе экспериментальных данных и результатов математического моделирования [6]. Алгоритмы наиболее известных подходов были заложены в разработанный авторами ГИС модуль реализации СММА.
Значения современной пористости и проницаемости, рассчитанные авторами по приведённым формулам для различных типов горных пород, показали хорошую сопоставимость с результатами лабораторных испытаний на основе петрофизических исследований [7]. Значения проницаемости, рассчитанные для выявленных структурных эпизодов, проявившихся на фоне давлений и температур, хорошо соотносятся с результатами лабораторных испытаний различных типов горных пород при воздействии на них эффективных давлений (до 200 МПа) и высоких температур (до 600 °С) [8].
Выводы
Реализация специальной методики микроструктурного анализа средствами ГИС технологий позволяет осуществлять картирование и маркировку различных типов микроструктур, а также оперативно измерять их пространственно-геометрические параметры непосредственно в ориентированных шлифах образцов горных пород. Комплексный подход к изучению трещинно-порового пространства на основе сочетания структурных и вещественных данных открывает возможности для реконструкции пространственно-временной связи между эпизодами деформаций и физико-химическими условиями протекания процесса эволюции рудоносной гидротермальной системы. Таким образом, при статистически значимом количестве образцов становится возможной реконструкция ориентировки осей ПНД, путей миграции рудоносных флюидов, расчета фильтрационных характеристик, формируемых определённой системой микроструктур или их совокупностью.
Разработанный подход применим для изучения многочисленных геологических объектов и месторождений, относящихся к различным промышленно-генетическим типам.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 18-35-00109.