CC BY
8
ÂecmAua геаАлуа, октябрь, 2020 г., № 10
УДК 552.122 DOI: 10.19П0^еоу.2020Л0.5
МОДЕЛИРОВАНИЕ ПЕТРОГРАФИЧЕСКИХ СТРУКТУР
Памяти д. г.-м. н. Ю. А. Ткачева
Ю. Л. Войтеховский1' 2, А. А. Захарова1
1Санкт-Петербургский горный университет, Санкт-Петербург 2Геологический институт ФИЦ КНЦ РАН, Апатиты Уоу(вккоузк1у_ УиЬ @рет.&рш I. ги
В статье изложена методика 2Э-моделирования структур биминеральных горных пород в терминах вероятностей бинарных межзерновых контактов. Показано, что типовые петрографические структуры (вкрапленная, цепочечная, пойкилито-вая, порфировая) и текстуры (массивная, полосчатая, шлировая) закономерно заполняют барицентрическую диаграмму (р,,, Рц, ру) вероятностей контактов. Указанные структуры и текстуры соотнесены с равновесием Харди — Вайнберга. Обсуждены проблемы моделирования полиминеральных петрографических структур и текстур. Предложен новый тип структурных диаграмм. Статья посвящена памяти д. г.-м. н. Ю. А. Ткачева — известного российского специалиста в области моделирования структур осадочных горных пород.
Ключевые слова: кристаллическая горная порода, петрографические структуры: вкрапленная, цепочечная, пойкили-товая, порфировая; петрографические текстуры: массивная, полосчатая, шлировая.
MODELING PETROGRAPHIC STRUCTURES
Dedicated to the memory of Dr. Sci. Yu. A. Tkachev
Yu. L. Voytekhovsky1' 2, A. A. Zakharova1
!Saint-Petersburg Mining University, Saint Petersburg 2Geological Institute FRC KSC RAS, Apatity
The article describes a technique for 2D modeling of structures and textures of bimineral rocks in terms of the probabilities of binary intergranular contacts. It is shown that typical petrographic structures (disseminated, chained, poikilitic, porphyric) and textures (massive, banded, schlieren) regularly fill the barycentric diagram (p¡¡, j pj of contact probabilities. These structures and textures are correlated with the Hardy-Weinberg equilibrium. The problems of modeling polymineral petrographic structures and textures are discussed. A new type of structural diagrams is suggested. The article is dedicated to the memory of Dr. Sci. Yu. A. Tkachev — a well-known Russian specialist in the field of modeling the structures of sedimentary rocks.
Keywords: crystalline rock, petrographic structures: disseminated, chained, poikilitic, porphyric;petrographic textures: massive, banded, schlieren.
Введение
Структура и текстура — атрибуты любой горной породы, содержащие разнородную информацию о ее генезисе. Петрографические словари определяют их как совокупные характеристики форм (панидио-, гипи-дио- и ксеноморфно-зернистые структуры), абсолютных и относительных размеров (тонко-, мелко-, средне- и крупно-, равномерно- и неравномерно-зернистые, порфировидные и порфировые структуры) слагающих горную породу минеральных зерен, а также пространственных отношений их самих (пойкилитовая структура) и агрегатов (массивная, шлировая, полосчатая текстуры) [5]. Отдавая должное поколениям петрографов, описавших сотни структур и текстур, заметим, что они не укладываются в математическую систему не столько из-за разнообразия, сколько из-за намерения вместить в одну характеристику несколько разнородных измеримых параметров и даже генетические (!) представления. Типичный пример — лепидогранобласто-вая структура, указывающая на чешуйки слюды, изо-
метричные зерна кварца и бластез горной породы. Мы видим в этом недостаток сегодняшней петрографической номенклатуры. (Можно ли представить номенклатуру точечных или пространственных групп симметрии, в которую входили бы утверждения о генезисе измеряемых кристаллов?) Тем не менее подчеркнем, что дальнейшие подходы предлагаются в дополнение, а не вместо существующих.
Принципы моделирования
Всякая модель представляет объект лишь в тех аспектах, которые полагаются всеобщими и принципиальными. Все более полное описание достигается пошаговым добавлением в модель более частных аспектов. В этом исследовании мы руководствовались следующими принципами.
Номенклатура петрографических структур и текстур сложилась при их наблюдении в 20-поверхностях: обнажениях, штуфах, шлифах под микроскопом. Моделирование в 3Б привело бы к сложной пробле-
Для цитирования: Войтеховский Ю. Л., Захарова А. А. Моделирование петрографических структур // Вестник геонаук. 2020. № 10(310). C. 38—42. DOI: 10.19110/geov.2020.10.5.
For citation: Voytekhovsky Yu. L., Zakharova A. A. Modeling petrographic structures. Vestnik of Geosciences, 2020, 10(310), pp. 38—42. DOI: 10.19110/geov.2020.10.5.
VesTnK oF GeoscEnces, October, 2020, No. 10
ме стереологической реконструкции, то есть распознаванию в сечениях. Поэтому моделирование выполнялось сразу в 2Б.
Горная порода — системный объект с определенными элементами (минеральными зернами различных видов) и пространственными отношениями (контактирования/неконтактирования). Сложность системы более зависит от специфики отношений, чем от разнообразия элементов (хотя ясно, что отношений без элементов не бывает). В этом исследовании мы избегаем описания элементов (форм и размеров минеральных зерен) и делаем акцент на описании их отношений контактирования.
Диагностика минеральных видов и характеристика форм и размеров минеральных индивидов — минералогический, а не петрографический уровень исследования. Последний состоит в узнавании горной породы по ее модальному составу (из которого следует химический состав, а если не следует, то потому, что составы минералов варьируют в широких пределах) и типу организации (структуре и текстуре), что в обоих аспектах подразумевает статистический подход.
Статистический подход содержится и в традиционном описании горной породы, ведь мы всегда отбираем «представительный» штуф, в котором достигается ее узнаваемый облик. В этом исследовании для характеристики пространственных отношений минеральных зерен мы используюем вероятности бинарных контактов. Моделирование выполнено для бими-неральной горной породы. В этом простейшем случае вероятности контактов и расчетные тренды можно изобразить в барицентрической диаграмме {ри, р^, Рц|. Возможности моделирования полиминеральных структур обсуждены отдельно.
Принципиально важная черта организации горных пород (а также металлов и сплавов) состоит в том, что 4 зерна контактируют в точке, 3 — по линейному ребру, 2 — по некоторой поверхности [1; 6, с. 34-35]. В случайное сечение первые не попадают, вторые выглядят как тройные контакты в точке, последние — как бинарные контакты вдоль линии (рис. 1). Но для любого планарного графа выполняется соотношение Эйлера: п — е + V = 1, где п — число клеток, е — ребер, V — вершин. В нашем случае (представляя граф без края — так мыслится и кристаллическая решетка) 3v = 2е. Тогда е = 3п — 3 » 3п при п ® ¥. Средняя координация клетки: 2е / п = 6 - 6/п ® 6 при п ® ¥.
Поэтому статистическое моделирование петрографических структур возможно на гексагональной сетке (рис. 1). На первый взгляд, это противоречит сказан-
ному выше об игнорировании форм минеральных зерен. Но для дальнейшего важны не гексагоны (их можно заменить точками в центрах), а усредненные характеристики сетки: общее число бинарных контактов 3п и средняя координация клетки 6. Историческая аналогия — переход О. Браве от форм «полиэдрических молекул» Р.-Ж. Гаюи к размещению в пространстве их центров, что привело к открытию 14 типов решеток. Удивительно, что правильное кристаллографическое разбиение плоскости на гексагоны проявилось «в среднем» в неправильном рисунке межзерновых границ в сечении горной породы.
Вкрапленные структуры, массивные текстуры
Пусть п — общее число минеральных зерен, тогда 3п — общее число бинарных контактов в структуре. Пусть кп — число черных зерен, рассеянных в белой матрице, здесь и далее 0 < к < 1 (оценка точной верхней границы для к всякий раз составляет задачу). Тогда числа различных бинарных контактов равны: = 0, = 6кп, = 3п - 6кп. Делением на 3п найдем соответствующие вероятности: рчч = 0, рчб = 2к, рбб = 1 - 2к. Предположительно, максимальная концентрация черных зерен достигается для к = 1А (рис. 2). Серия вкрапленных структур для 0 < к < 1А образует ограниченный непрерывный тренд на ребре (рбб, рчб) барицентрической диаграммы (рис. 3). Равномерное рассеяние черных зерен среди белых создает массивную текстуру.
Пойкилитовые структуры
Пусть п — общее число ойкокристаллов, т — хада-кристаллов (рис. 2). Общее число контактов: = 0, = т, = 3п. Соответствующие вероятности: рчч = 0, рчб = т/(3п+т) = 1/(1+3п/т), рбб = 3п/(3п+т) = 1/(1+т/3п). Если хадакристаллов относительно мало, т. е. т/3п ® 0, то рчб ® 0, рбб ® 1. Если хадакристаллов относительно много, т. е. 3п/т ® 0, то рчб ® 1, рбб ® 0. Тренд пойкилитовых структур занимает все ребро (рбб, рчб) барицентрической диаграммы (рис. 3).
Цепочечные структуры, полосчатые текстуры
Пусть п — общее число минеральных зерен, тогда 3п — общее число бинарных контактов в структуре. Пусть кп — число черных зерен в цепочках. Тогда = кп, = 4кп, = 3п - 5кп. Соответствующие
Рис. 1. Межзерновые границы в петрографическом шлифе и гексагональная сетка Fig. 1. Intergranular boundaries in a petrographic section and a hexagonal net
ÂecfoUwtc teaHàtyc, октябрь, 2020 г., № 10
Рис. 2. Модели структур и текстур: 1 — вкрапленная, 2 — пойкилитовая, 3 — цепочечная, 4 — с ветвящимися кластерами, 5 — шлировая, 6 — порфировая
Fig. 2. The models of the structures and textures: 1 — disseminated, 2 — poikilitic, 3 — chained, 4 — with brunching clusters, 5 — schlieren, 6 — porphyric
вероятности: рчч = к/3, рчб = 4к/3, рбб = 1 — 5к/3. Предположительно, максимальная концентрация черных зерен достигается для к = У2 (рис. 2). Серия цепочечных структур для диапазона 0 < к < У2 образует непрерывный тренд (рис. 3). Превращением цепочек в черные полосы и варьированием (через к) ширины белых полос между ними моделируются полосчатые текстуры. Их тренды покрывают угол барицентрической диаграммы ниже тренда цепочечной структуры.
Структуры с ветвящимися кластерами
В работе [2, с. 10—13] описаны структуры с ветвящимися кластерами зерен кварца в полевошпатовой матрице. Выполним моделирование для кластеров из 7 зерен (рис. 2). Пусть п — общее число минеральных зерен в структуре, тогда 3п — общее число бинарных контактов. Пусть кп — общее число кварцевых зерен в кластерах. Тогда число кластеров равно кп/7. Число контактов: = кп, = 4кп, = 3п - 5кп. Соответствующие вероятности: рчч = к/3, рчб = 4к/3, Рбб = 1 - 5к/3. Результат совпал с предыдущим, так как цепочки и ветвящиеся кластеры отличаются лишь одним зерном — центром ветвления, что не сказалось на статистике контактов.
Шлировые текстуры
Пусть п — общее число минеральных зерен, тогда 3п — общее число бинарных контактов в структуре. Пусть кп — число зерен, собранных в шлиры по 7
Рис. 3. Барицентрическая диаграмма вероятностей (ри, pjj, pij). Структуры и текстуры: 1 — вкрапленная, 2 — цепочечная, 3 — порфировые (d — диаметр фенокриста), 4 — шлировые (количества зерен в шлирах указаны числами), 5 — с ветвящимися кластерами, 6 — пойкилитовые; 7 — линия равновесия Харди — Вайнберга
Fig. 3. Barycentric diagram of (ри, pjj, pjj) probabilities. Structures and textures: 1 — disseminated, 2 — chained, 3 — porphyric (d — diameter of phenocryst), 4 — schlieren (quantities of grains in the schlierens are indicated by the numbers), 5 — with brunching clusters, 6 — poikilitic; 7 — the Hardy—Weinberg equilibrium line
VsNK F GeoSciEnceS, October, 2020, No. 10
(рис. 2). Тогда число шлиров в структуре кп/7, число контактов: Кчч = 12кп/7, Кчб = 18кп/7, Кбб = 3п — 30кп/7, соответствующие вероятности: рчч = 4к/7, рчб = 6к/7, рбб = 1 — 10к/7. Тренды шлировых текстур для разныж размеров шлиров и значений к покрывают барицентрическую диаграмму ниже тренда цепочечной структуры (рис. 3).
Порфировидные и порфировые
структуры
Эти структуры удобно моделировать как предыдущие, у которых в шлирах стерты межзерновые границы. Характеристика модельных шлиров дана в таблице. Здесь d — диаметр шлира (диаметр белого зерна принят за 1), п^) — число зерен в шлире (аппроксимируется величиной я^2/4), у^) = Кчч — число границ в шлире, р^) = Кчб = 6d — число внешних границ шлира, Е = + Кчб /2 = 3п^).
Характеристика модельных шлиров Characteristics of model shlierens
d 3 5 7 9 11
n(d) 7 19 37 61 91
v(d) 12 42 90 156 240
P(d) 18 30 42 54 66
S 21 57 111 183 273
Пусть п — общее число минеральных зерен, тогда 3п — общее число бинарных контактов в структуре. Пусть кп — число зерен в шлирах, тогда 3кп — числа контактов в кп/п^)-шлирах. Стирая внутренние границы, превратим шлиры в порфировые вкрапленники. Число оставшихся контактов: Кчч = 0, Кбб = 3п — 3кп — кпр^) / 2п^) = 3п(1-к) — 12кп / pd, Кчб = 2knp(d) / 3n(d) = 16кп / лd. Деля на их сумму, получим вероятности: рчч = 0, рчб = 16к / [3я^(1-к) + 4к], рбб = 1 — рчб.
Тренды порфировидных и порфировых структур в совокупности покрывают ребро (рбб, рчб) барицентрической диаграммы (рис. 3). Каждый начинается вблизи вершины рбб и простирается к вершине рчб тем далее, чем меньше диаметр кластера. При фиксированном d предельное значение к вычисляется из условия рчб > 0 (рис. 4).
Обсуждение и выводы
Выше предложена барицентрическая диаграмма нового типа, показывающая вероятности бинарных межзерновых контактов (в простейшем случае бими-неральной горной породы), существующих во всех кристаллических горных породах и тем самым выражающих фундаментальные пространственные отношения между слагающими их минеральными индивидами. При заданном виде кластера моделируются любые петрографические структуры и текстуры. При этом их канонические типы закономерно заполняют поле диаграммы трендами, определяемыми видом кластера и его концентрацией в структуре. (Левая часть диаграммы на рис. 3 симметрично заполняется структурами «белые кластеры в черной матрице».) Фигуративные
Рис. 4. Соотношение d и k для порфировых структур Fig. 4. Relationship between d and k for porphyric structures
точки реальных петрографических структур и текстур должны занимать промежуточные положения между модельными трендами.
Около вершины рчб внутри диаграммы есть поле, в которое модельные тренды не проникают. Это означает, что невозможна петрографическая структура без мономинеральных контактов (типа шахматной доски), образуемых хотя бы одним минералом (таковы модельные вкрапленные и пойкилитовые структуры). В реальных петрографических структурах обычны мономинеральные контакты обоих видов. Для расчетных трендов могут быть оценены предельные точки, вместе ограничивающие запрещенное поле.
Тренды различных структур (например, цепочечной и с ветвящимися кластерами) и поля текстур (полосчатых и шлировых) в поле диаграммы перекрываются. Совпадение статистик бинарных контактов и, как следствие, термодинамических характеристик [3] сети межзерновых границ обращает внимание на фундаментальный вопрос: почему горная порода следует той или иной стратегии кристаллизации, приводящей к той или иной структуре?
Интересен вопрос о соотношении различных структур и текстур с рассмотренным ранее [4] равновесием Харди - Вайнберга, выражающим гомогенное перемешивание минеральных индивидов в горной породе (массивная текстура). На диаграмме видно, что тренды вкрапленных, пойкилитовых и цепочечных структур, начинаясь вблизи вершины рбб диаграммы, на значительном протяжении следуют вблизи линии равновесия. Явно отклоняются от нее ярко выраженные порфировые, пойкилитовые и цепочечные структуры, а также полосчатые и шлировые текстуры.
Возможно ли расширение метода на полиминеральные структуры? Принципиальных препятствий к этому нет. Проблема состоит в визуализации барицентрических диаграмм и поверхностей равновесия Харди — Вайнберга. Даже для трехминеральной горной породы это 2-мерная поверхность в 6-вершинном (5-мерном) симплексе. В этом случае более наглядно последовательное рассмотрение полиминеральной структуры в виде «минерал — матрица» для каждого из слагающих горную породу минералов и сравнение результатов на барицентрической диаграмме.
ÂecmAute геаАлук, октябрь, 2020 г., № 10
Расчет вероятностей бинарных межзерновых контактов — рутинная процедура даже для биминеральной горной породы. Мы изучаем возможность применения для этих целей линейного метода Розиваля в компьютерном варианте и с применением отечественного анализатора структуры МИУ-5.
Литература
1. Беленький В. 3. Геометрико-вероятностные модели кристаллизации. Феноменологический подход. М.: Наука, 1980. 88 с.
2. Бескин С. М, Ларин В. Н., Марин Ю. Б. Редкометал-льные гранитовые формации. Л.: Недра, 1979. 280 с.
3. Бродская Р. Л. Термодинамические (кинетические) критерии формирования и эволюции структуры минеральных агрегатов // Зап. ВМО. 1988. № 5. С. 623-633.
4. Войтеховский Ю. Л. Типизация неправильных пространственных разбиений // Вестник ИГ Коми НЦ УрО РАН. 2018. № 8. С. 46-50. <М: 10.19110/2221-1381-2018-8-46-50.
5. Левинсон-Лессинг Ф. Ю, Струве Э. А. Петрографический словарь. М.: Госгеолиздат, 1963. 448 с.
6. Салтыков С. А. Стереометрическая металлография. М.: Гос. науч.-тех. изд. лит. по черной и цветной металлургии, 1958. 447 с.
References
1. Belenkiy V. Z. Geometriko-veroyatnostnye modeli kristallizatsii. Fenomenologicheskypodhod (Geometric-probabilistic models of crystallization. Phenomenological approach). Moscow: Nauka, 1980, 88 pp.
2. Beskin S. M., Larin V. N., Marin Yu. B. Redkometalnye granitovye formatsii (Rare-metal granitic formations). Leningrad: Nedra, 1979, 280 pp.
3. Brodskaya R. L. Termodinamicheskiye (kineticheskiye) kriterii formirovaniya i evolutsii struktury mineralnykh agregatov (Thermodynamic (kinematic) criteria of formation and evolution of mineral aggregates structures). Proc. Rus. Mineral. Soc, 1988, No. 5, pp. 623-633.
4. Voytekhovsky Yu. L. Tipizaciya nepravilnykh prostranstven-nykh razbieniy (Typization of irregular spatial tilings). Vestnik IG Komi SC UB RAS, 2018, No. 8, pp. 46-50. doi: 10.19110/22211381-2018-8-46-50.
5. Levinson-Lessing F. Yu., Struve E. A. Petrograficheskiy slovar (Petrographic dictionary). Moscow: Gosgeolizdat, 1963, 448 pp.
6. Saltykov S. A. Stereometricheskaya metallografya (Stereometric metallography). Moscow: State Publ. House on Ferrous and Non-Ferrous Metallurgy, 1958, 447 pp.
Поступила в редакцию / Received 19.10.2020
