Структура контрастной расслоенности "критических" зон базит-гипербазитовых интрузий Кивакка и Бушвельд: Свойства, происхождение тема диссертации и автореферата по ВАК РФ 04.00.08, кандидат геолого-минералогических наук Котов, Сергей Робертович

• Котов, Сергей Робертович
• кандидат геолого-минералогических наук
• 1998

Оглавление диссертации кандидат геолого-минералогических наук Котов, Сергей Робертович

1 Математическое моделирование как метод решения геологических задач

1.1 Модели и их классификация.

1.2 Задачи математической геологии и концептуальные вероятностные модели

2 Модели формирования расслоенности в мафических интрузиях

2.1 Внедрение магмы

2.2 Конвективное перемешивание.

2.3 Механические процессы.

2.4 Механизмы, являющиеся следствием изменения интенсивных параметров.

2.5 Механизмы поздней стадии кристаллизации и охлаждения.

3 Краткая геолого-петрографическая характеристика исследуемых объектов

3.1 Массив Кивакка.

3.1.1 Геологическое положение массива Кивакка в составе Олангской группы интрузий.

3.1.2 Состав и строение интрузии.

3.1.3 Минералогия.

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

3.2 Бушвельдский комплекс.

3.2.1 Геологическое строение.

3.2.2 Минералогия.

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

4 Стохастические и геометрические характеристики ритмической расслоенности

4.1 Контрастная расслоенность в рамках представлений о последовательном слоеобразовании.

4.2 Контрастная расслоенность в рамках представлений об объемном слоеобразовании: геометрический подход

5 Контрастная расслоенность как результат объемного (жидкостного) расслоения магмы: концептуальная вероятностная модель

5.1 Концептуальная схема и математическая формализация

5.2 Тестирование модели.

5.3 Фрактальная размерность как прогнозный критерий.

6 Структурно-вещественные особенности норитовой зоны массива Кивакка, связь с моделью

6.1 Распределение микроэлементов в рамках модели фракционной кристаллизации.

6.2 Структурные взаимоотношения минералов в породах расслоенной серии.

6.3 Элементные и минеральные ассоциации и определяющие их факторы.

Введение диссертации (часть автореферата) на тему «Структура контрастной расслоенности "критических" зон базит-гипербазитовых интрузий Кивакка и Бушвельд: Свойства, происхождение»

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

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

Очевидно, что понимание механизма формирования ритмично - расслоенного строения основных - ультраосновных интрузивов представляет и теоретический, и практический интерес. Неслучайно уже не одно десятилетие в геологическом мире продолжается острая дискуссия по этому поводу. При этом привлекаются самые разнообразные, зачастую взаимоисключающие гипотезы: это представления о повторяющемся послойном внедрении новых порций магмы, либо внедрении уже расслоенного в магматическом очаге расплава; это всевозможные модели направленной (от подошвы к кровле) кристаллизации, осложненной внешними факторами (например, непрерывным или периодическим конвективным перемешиванием); это модели кристаллизационной дифференциации, в которых в качестве основного механизма постулируется разделение во всем объеме интрузива жидкой и новообразующейся твердой фаз в гравитационном поле; это механизмы, связанные с вариациями каких либо интенсивных параметров системы: флуктуации скоростей зародышеобразования, контролируемые диффузией нуклеация и рост кристаллов, рост кристаллов в температурном градиенте, флуктуации фугитивности кислорода или общего давления; это ликвация и жидкостное разделение силикатного расплава; это представления о перераспределении вещества за счет газовой конвекции, о вторичном изменении первично - осадочных слоистых пород или об изменении вдоль системы параллельных трещин и др. [60, 98].

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

Цель работы

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

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

В части предварительного изучения свойств моделируемых объектов и проверки согласия модели работа опирается на результаты экспедиционных работ 1994-1996 г.г. на массиве Кивакка (Сев. Карелия, Россия) и Бушвельдском комплексе (ЮАР), проводившихся при непосредственном участии автора в соответствии с плановой темой лаборатории Математической геологии ИГГД РАН и Института Математической Геологии РАЕН в рамках международного (Россия, Греция, ЮАР) проекта "Хромит" по заказу министерства Науки и Технологий РФ.

Геометрико-вероятностные свойства рисунков расслоенности исследовались на материале двух детальных пересечений через расслоенную серию массива Кивакка и пяти скважин через "критическую" зону Буш-вельдского комплекса. Общий объем изученных выборок составляет около четырех тысяч замеров мощностей слоев с идентифицированными составами. Также в работе были использованы результаты исследования около 400 шлифов, 30 количественных рентгенофлюоресцентных анализов на Ре, №, Со, Сг, Т1, V, Ъг, У, Бг, Ш). РЬ и 118 количественных рентгеноспек-тральных анализов пород массива Кивакка (Са, Сг, Си, Ре, Оа. Мп, №, Ш>, Бг. Т1, У, Zn), 18 микрозондовых определений составов породообразующих минералов.

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

1. Впервые установлены стохастические и геометрические ограничения на модели формирования контрастной расслоенности. В частности, выявлены фрактальная структура рисунков расслоенности и стохастическое дальнодействие в последовательностях мощностей слоев (ритмов).

2. Впервые разработана концептуальная вероятностная модель объемного (жидкостного) расслоения магматического расплава, показано ее согласие с реальными объектами.

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

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

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

Основные защищаемые положения

1. Ритмическая расслоенность "критических" зон массивов Киваккаи Буш-вельд не может быть результатом простой периодической динамики процесса слоеобразования; в то же время формирование расслоенности определялось действием малого количества (скорее всего не более трех) независимых процессов.

2. Динамика процессов слоеобразования в "критических" зонах массивов Кивакка и Бушвельд обеспечивает геометрические и вероятностные свойства рисунков контрастной расслоенности: 1) фрактальную структуру; 2) стохастическое "дальнодействие", которое характеризуется квазипериодическими автокорреляционными функциями по последовательностям мощностей слоев (ритмов).

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

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

Апробация работы

Основные результаты исследований по теме диссертации докладывались на научных заседаниях лабораторий ИГГД РАН, института математической геологии РАЕН, на международном совещании по проблемам докембрия Европы (Санкт-Петербург, 1995), на международном совещании по проблемам докембрия Северной Евразии (Санкт-Петербург, 1997), на международной конференции по проблемам генезиса магматических и метаморфических пород (Санкт-Петербург, 1998). По теме диссертации опубликовано 6 работ. Внедрение результатов исследований производилось в виде их использования при составлении научных отчетов по тематическим работам института математической геологии (Вистелиус и др., 1995; Деч и др., 1996). Наработки по проблеме использовались автором при чтении лекций по курсам "Математические методы в геологии", "Математическое моделирование природных процессов", "Синергетика" для бакалавров и магистров геологического факультета СПбГУ.

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

Диссертация состоит из Введения, Заключения и 6-ти глав. В связи с тем, что математическое моделирование пока еще не стало привычным методом в геологии, глава 1 посвящена классификации методов моделирования и месту предлагаемой модели в общем их множестве. В главе 2 приводится краткий обзор основных предлагаемых моделей формирования расслоенных интрузий. Краткой геолого-петрографической характеристике исследуемых объектов, а также описанию используемого фактического материала посвящена глава 3. В главе 4 рассматриваются стохастические и геометрические особенности ритмической расслоенности, выявленные для изучаемых объектов. Предполагается, что они должны являться дополнительными ограничениями при создании и математической формализации любой выбранной модели. Концептуальная схема модели, математическая формализация, а также результаты проверки согласия модели тестовым объектам рассматриваются в главе 5. В главе 6 на основе изучения структурно-вещественных особенностей массива Кивакка определяется место и доля участия процесса, отвечающего концептуальной схеме модели.

Заключение диссертации по теме «Петрография, вулканология», Котов, Сергей Робертович

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

А. На ранних этапах — ликвация магматического расплава в объеме расслоенной серии на более основную ("ортопироксенитовую") и менее основную ("базальтоидную") жидкости; их доли составляли 10-15% и 90-85% соответственно. На отдельных уровнях магматической камеры ликваты перераспределялись с образованием относительно меланократовых (40-60% "ортопироксенитовой" жидкости) и лейкократовых (7-0.5% "ортопироксенитовой" жидкости) субгоризонтальных слоев. Последующая кристаллизация обеспечивала на этих уровнях ритмическое чередование норитов, меланоноритов, лейконоритов.

Б. На этапе кристаллизации формирование лейкократовых и мелано-кратовых слоев обеспечивалось также фракционированием кристаллизующихся из "базальтоидной" жидкости фаз (OpxII, Pill, CpxII). Доля OpxII для нефракционированных составов варьирует от 0.33 до 0.22 (в среднем — 0.25, табл.6). Фракционирование приводит к изменению этого соотношения в закристаллизовавшихся породах (0.14-0.17 в лейкократовых и 0.60-0.86 в меланократовых). Фракционирование проявляется в пределах локальных уровней магматической камеры, не испытавших, как правило, предварительного ликвационного расслоения, и приводит к формированию норитов, меланоноритов, лейконоритов (анортозитов). В участках сопряжения обоих процессов — расслоения двух ликвировавших жидкостей и фракционирования, — формируется наиболее контрастная ритмика (пироксенит, меланонорит, лейконорит) [26].

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

По всей видимости лучше всего известные особенности строения "критической" зоны объясняются в рамках идеи о формировании отдельных циклов путем внедрения более примитивных порций магмы в магму, кристаллизующую плагиоклаз и пироксен [87]. Известно, что при внедрении, менее вязких жидкостей в более вязкие обычным явлением бывает образование так называемых "вязких пальцев" одной жидкости в другой, имеющих фрактальную структуру ([59], глава 4). Можно предположить, что сходная ситуация возникает при внедрении новой порции магматического расплава в относительно более холодную и вязкую жидко-кристаллическую магматическую кашу. Это представление таким образом согласуется с выявленными в работе геометрическими и вероятностными особенностями рисунка расслоенности. Определенное участие могли также принимать и процессы жидкостной дифференциации. Предложенная вероятностная модель может использоваться для Бушвельдского массива как феноменологический описательный инструмент.

7 Основные результаты и выводы (Заключение)

На основе изучения вероятностных и геометрических свойств контрастной расслоенности "критических" зон массивов Кивакка и Бушвельд установлено, что она не может быть результатом простой периодической динамики слоеобразования; в то же время формирование расслоенности определялось действием малого количества (скорее всего не более трех) независимых процессов.

Разрезы массивов Кивакка и Бушвельд характеризуются квазипериодическими автокорреляционными зависимостями по последовательностям мощностей слоев и ритмов и фрактальной структурой рисунка контрастной расслоенности. Для массива Кивакка локальные максимумы автокорреляционных функций соответствуют лагам 75, 175 (профиль 101), 24, 45, 70 (профиль107), для Бушвельда — лагам 24, 45, 70 (скважина Ш6). Оценка фрактальной размерности для "критической" зоны массива Кивакка и Бушвельд составляет соответственно Б = 0.44 ± 0.03 и Б = 0.26 ± 0.04. Предложенная математическая модель объемного расслоения расплава хорошо согласуется с такой геометрией контрастной расслоенности "критических" зон.

Список литературы диссертационного исследования кандидат геолого-минералогических наук Котов, Сергей Робертович, 1998 год

1. Андреев Н.С , Мазурин О.В., Порай-Кошиц Е.А., Роскова Г.П., Фи-липович В.Н. Явления ликвации в стеклах. J1. Наука, 1974. 220 с.

2. Барков А.Ю., Леднев А.И., Клюнин С.Ф., Захаров A.A. Платиноме-тальная минерализация в массиве Кивакка Северной Карелии //в кн. Новые данные по минералогии редких элементов Кольского п-ва. Апатиты: КНЦ АН СССР, 1991. с.87-89

3. Батороев К.Б. Аналогии и модели в познании. Новосибирск: Наука, 1981. 319с.

4. Безмен Н.И., Сук Н.И. Базит-гипербазитовое расслоение (эксперимент) // Геология и разведка, 1983, N6., с.48-52

5. Безмен Н.И., Персиков Э.С., Сук Н.И. Кинетические особенности расслоения основных магматических расплавов при их плавлении под давлением водорода, //в сб: Динамические модели физической геохимии. Новосибирск: Наука, 1982, с. 129-134

6. Боуэн Н.Л. Эволюция изверженных пород. ОНТИ, М. Л., 1934

7. Брусиловский Б.Я. Теория систем и система теорий. Киев: Вища школа, 1977. 191 с.

8. Вентцель Е.С., Овчаров Л.А. Теория случайных процессов и ее инженерные приложения. М: Наука, 1991. 383 с.

9. Воробьева O.A. О первичной расслоенности Ловозерского щелочного массива // в сб. Производительные силы Кольского полуострова, т.1. М: Изд-во АН СССР, 1940

10. Воронин Ю.А., Еганов Э.А. Методологические вопросы применения математических методов в геологии. Новосибирск: Наука, 1974. 86 с.

11. Высоцкий Ю.А., Котов С.Р. Структура расслоенности базитового массива Кивакка (Северная Карелия) // Докембрий Северной Евразии. Тез. докл. Санкт-Петербург, 1997. С.24.

12. Галенко П.К. К феноменологичекой теории локально-неравновесной кристаллизации сплавов // ДАН, 1994, т.334, N.6, с. 707-709

13. Генкин А.Д., Дистлер В.В., Лапутина И.П. Хромитовая минерализация дифференцированных трапповых интрузий и условия ее образования / / Условия образования магматических рудных месторождений. М.: Наука, 1979. С.105-126

14. Горстко А.Б. Познакомьтесь с математическим моделированием. М: Знание, 1991. 160 с.

15. Горяинов П.М., Егоров Д.Г., Иванюк Г.Ю. О структурно-вещественной самоорганизации в архейских железорудных ансамблях (Кольский полуостров) // ДАН. 1992. Т.322, N6.

16. Груза В.В. Методологические проблемы геологии. Ленинград: "Недра", 1977. 181 с.

17. Делицын Л.М., Мелентьев Б.Н. Механизм формирования полосчатых горных пород // Бюл. МОИП. Отд. геол. 1979, Т.54, Вып.4, С.102-113.

18. Иванцов Г.П. "Диффузионное" переохлаждение при кристаллизации бинарного сплава // ДАН, 1951, t.LXXXI, N.2, с. 179-182

19. Карри X. Основания математической логики. М: Мир, 1969. 568 с.

20. Клюнин С.Ф., Гроховская T.JL, Захаров A.A., Соловьева Т.В. Геология и перспективы платиноносности Олангской группы массивов (Сев. Карелия) // в сб. Геология и генезис месторождений платиновых металлов, М: Наука, 1994. 301 с.

21. Кноринг Л.Д., Деч В.Н. Геологу о математике. Ленинград:"Недра", 1989. 208 с.

22. Котов С.Р. Стохастические и геометрические ограничения на модели формирования расслоенности в базитовых интрузиях. СПб., 1998, Деп.в ВИНИТИ N659-B98 от 06.03.98.-13с.

23. Крамбейн У., Кауфмен М., Мак-Кеммон Р. Модели геологических процессов (введение в математическую геологию). М: Мир, 1973. 150 с.

24. Крамбейн У., Грейбилл Ф. Статистические модели в геологии. М: Мир, 1969. 398 с.

25. Красногцеков П.С., Петров A.A. Принципы построения моделей. М: Изд-во МГУ, 1983. 264 с.

26. Кульбак С. Теория информации и статистика. М: Наука, 1967, 408с.

27. Лавров М.М. Гипербазиты и расслоенные перидотит габбро - нори-товые интрузии докембрия Северной Карелии. Л.: Наука, 1979. 136 с.

28. Лебедев А.П. О типах дифференциации в траппах Сибирской платформы // Изв. АН СССР. Сер.геол., N2, 1957, с.55-74

29. Лукашевич В.К. Модели и метод моделирования в человеческой деятельности. Минск: Наука и техника, 1983. 119 с.

30. Магматизм и металлогения рифтогенных систем восточной части Балтийского щита (ред. акад. А.Д.Щеглова). СПб: Недра, 1993. 244с.

31. Маракушев A.A., Иванов И.П., Римкевич B.C. Экспериментальное воспроизведение ритмической магматической расслоенности // Докл.АН. 1981. Т.258. N.1. С.183-185.

32. Маракушев A.A., Безмен Н.И. Эволюция метеоритного вещества, планет и магматических серий. М: Наука, 1983. 185 с.

33. Минерал ы (справочник), ред. акад. Ф.В.Чухрова. М: Наука, 1981.

34. Моделирование и познание. (Под ред. В.А.Штоффа). Минск: Наука и техника, 1974, 211 с.

35. Новик И.Б. О моделировании сложных систем. М: Мысль, 1965. 335 с.

36. Осипов М.А. Формирование расслоенных плутонов с позиций термоусадки. М: Наука, 1982, 100 с.

37. Очерки физико-химической петрологии (экспериментальное исследование проблем магматизма). Вып. XIII. ред. Жариков В.А., Федькин В.В., М: Наука, 1985, 208 с.

38. Прокопцев Н.Г. Рудообразование в связи с ликвацией базальтовой магмы // Изв. АН СССР. 1990. Серия геол. N 11. С. 66-78.

39. Прусаков Г.М. Математические модели и методы в расчетах на ЭВМ. М: Наука, 1993. 141 с.

40. Седов JI.И. Теоретические модели. СКОПЛЕ, 1975. 27 с.

41. Седов Л.И. Механика сплошной среды. М: Наука, 1976. 2т. 528, 560 с.

42. Синицин A.B. Региональная тектоника и металлогения раннего докембрия. Л: Недра, 1990, 493 с.

43. Скиба В.И., Пушкарев Ю.Д. Физико-химическая природа расслоен-ности интрузивных тел // ДАН, 1985, т.280, N.2

44. Скрипниченко В.А. О природе ритмичной расслоенности базит-гипербазитового массива Кивакка (Северная Карелия) // Бюл. МО-ИП, отд. геол., 1984, т.59, вып.5, с.46-54

45. Смолькин В.Ф. Коматиитовый и пикритовый магматизм раннего докембрия Балтийского щита. СПб: Наука, 1992. 273 с.

46. Смолькин В.Ф., Борисов А.Е., МаракушевА.А. Признаки дифференциации и несмешивания расплавов в пикрито-базальтах Печенги // ДАН, 1987, т.294, N.3, с. 669-673

47. Турбин А.Ф., Працевитый Н.В. Фрактальные множества, функции, распределения. Киев: Наукова думка, 1992. 207 с.

48. Уэйджер JI., Браун Г. Расслоенные изверженные породы. М.: Мир, 1970. 529 с.

49. Фракталы в физике. Под ред. Л.Пьетронеро и Э.Тозатти). М: Мир, 1988. 672 с.

50. Федер Е. Фракталы. М: Мир, 1991. 260 с.

51. Френкель М.Я., Ярошевский A.A., Арискин A.A., Бармина Г.С., Коптев-Дворников Е.В., Киреев B.C. Динамика внутрикамерной дифференциации базитовых магм. М: Наука, 1988. 216 с.

52. Френкель М.Я., Ярошевский A.A., Коптев-Дворников Е.В., Арискин A.A., Киреев B.C., Бармина Г.С. Кристаллизационный механизм формирования расслоенности пластовых интрузивов // Записки ВМО, вып.З, 1985, часть 114, с.257-274

53. Френкель М.Я. Формулировка и алгоритм решения на ЭВМ системы уравнений переноса тепла и вещества в процессе становления пластового интрузива // Геохимия, N.4, 1978, с.547-559

54. Хорафас Д.Н. Системы и моделирование. М: Мир, 1967. 420 с.

55. Чехмир A.C., Симакин А.Г., Эпельбаум М.Б. Динамические явления во флюидно- магматических системах. М: Наука, 1991. 142 с.

56. Шапкин А.И. Автоколебательный режим общей конвекции в магматической камере и процесс ритмического расслоения магматического интрузива// Геохимия, 1993, N.9, с. 1297-1307

57. Шапкин А.И. Ритмическая седиментация кристаллической фазы в магматической камере // Геохимия, 1992, N.7, с. 939-949

58. Шарапов В.Н., Черепанов А.Н. Динамика дифференциации магм. Новосибирск: Наука, 1986. 188 с.

59. Шарков Е.В. Происхождение критических зон крупных расслоенных интрузивов //В кн. Геология и генезис месторождений платиновых металлов. М: Наука, 1994. с.35-45

60. Шарков Е.В. Петрология расслоенных интрузий. JI: Наука, 1980. 120 с.

61. Штофф В.А. О роли моделей в познании. JI: Изд-во ЛГУ, 1963. 128 с.

62. Шустер Г. Детерминированный хаос. М: Мир, 1988. 240 с.

63. Эволюция изверженных пород, ред. Х.Йодера. М: Мир, 1983. 528 с.

64. Юдин Б.А. Окисные железо-титановые и железные руды магматических формаций Карелии и Кольского полуострова. Петрозаводск: Карелия, 1987. 209 с.

65. Ярошевский А.А. О происхождении ритмических структур изверженных горных пород // Геохимия, 1970, N.5, с.562-574

66. Amelin Yu.V., Heaman L.M., Semenov V.S. U-Pb geochronology of layered mafic intrusions in the eastern Baltic shield: implications for the timing and duration of paleoproterozoic continental rifting // Precamb. Res., 1995, 75, p.31-46

67. Balashov Yu.A., Bayanova T.B., Mitrofanov F.P. Isotope data on the age and genesis of layered basic-ultrabasic intrusions in the Kola Peninsula and northern Karelia, northeastern Baltic Shield // Precambrian Research, 64, 1993, p.197-205

68. Bezmen N.I., Elevich V.Ya. Petrochemical types of the PGE-bearing massifs and their pedogenesis, //in book: International platinum (ed. by Laverov N.P., Distler V.V.), S.-Petersburg: Theophrastus publ., 1998, p. 8-17

69. Bezmen N.I. Hydrogen in magmatic systems // Experiment in Geosciences. 1992. Vol.1. N.2 P. 1-33.

70. Cameron E.N. The Lower zone of the Eastern Bushveld Complex in the Olifants River trough // J.Petrology, 1978, p.437-462

71. Campbell I.H. A study of macro-rhythmic layering and cumulate processes in the Jimberlana intrusion, western Australia. Parti: The Upper Layered Series // J.Petrology 18, 1977, p. 185-215

72. Cheng Q. Multifractal modeling and lacunarity analisis // Mathematical Geology. 1997. V.29. N.7. p. 919-932

73. Drever H.I. Immiscibility in the pictitic intrusion at Igdlorssuit West Greenland // Intern. Geol. Congr. 21 Sess. Proc.pt XIII. Copenhagen, 1960. P.47-58

74. Eales H.V., Cawthorn R.G. The Bushveld Complex// in book: Layered intrusions. (Ed. by Cawthorn R.). Amsterdam: Elsevier, 1996. p.181-229

75. Gorring M.L., Naslund H.R. Geochemical reversals within the lower 100 m of the Palisades sill, New Jersey // Contr. Miner. Petrol. 119, 1995, p.263-276

76. Hort M., Marsh B.D., Spohn T. Igneous layering through oscillatoty nucleation and crystal settling in well-mixed magmas.// Contrib. Mineral. Petrol. 1993. Vol. 114. P.425-440.

77. Hughes C.J. Igneous petrology. Amsterdam: Elsevier, 1982. 551 p.

78. Irvin T.N. Crystallization sequences in the Muskox intrusion and other layered intrusions: II Origin of chromitite layers and similar deposits of other magmatic ores. Geochim. Cosmochim. Acta., 39, 1975, p.991-1020

79. Irvine T.N., Keith D.W., Todd S.G. The J-M platinum-palladium reef of the Stillwater complex, Montana: II. Origin by double-diffusive convective magma mixing and implications for the Bushveld complex // Economic geology, v.78, N.7, 1983, p.1287-1334

80. Irvin T.N. Terminology for layered intrusions // Journal of Petrology, v.23,N.2, 1982, p.127-162

81. Ivanov D.N.,Kotov S.R. Metablastesis as the texture-forming process in some tonalites // Theophrastus Contributions to Advansed Studies in Geology, Athens: Theophrastus Publication, 1996, v.l, p.113-124

82. Layered intrusions. (Ed. by Cawthorn R.). Amsterdam: Elsevier, 1996. 531 p.

83. Mandelbrot B.B. The fractal geometry of nature. New-York: Freeman, 1983. 540 p.

84. Marsh B.D. Crystal capture, sorting, and retention in convecting magma // Geol. Soc. Am. Bull. 100, 1988, p. 1720-1737

85. Marsh B.D. On convective style and vigor in sheet-like magma chambers // J.Petrology 30, 1989, p. 479-530

86. McBirney A.R., Nakamura Y. Immiscibility in late-stage magmas of the Skaergaard intrusion // Yrbk. Carnegie Inst. Wash., 1974, p.348-352

87. McBirney A.R., Noyes R.M. Crystallization and layering of the Skaergaard intrusion // J.Petrology 20, 1979, p. 487-554

88. McBirney A.R. Further considerations of double-diffusive stratification and layering in the Skaergaard intrusion // J.Petrology 26, 1985, p. 9931001

89. Naslund H.R., McBirney A.R. Mechanisms of formation of igneous layering // in book: Layered intrusions. (Ed. by Cawthorn R.). Amsterdam: Elsevier, 1996. p. 1-43

90. Rollinson H. Using geochemical data: evaluation, presentation, interpretation. Ney-York: Longman, 1993. 350 p.

91. Sparks R.S., Huppert H.E., Koyaguchi T., Hallworth M.A. Origin of modal and rhythmic igneous layering by sedimentation in a convecting magma chamber // Nature, 361, 1993, p.246-249

92. Turchenko S.I. Precambrian metallogeny related to tectonics in the eastern part of the Baltic Shield // Precambrian Research, 58, 1992, p.121-141

93. Watson E.B. Two-liquid partition coefficients: experimental data and geochemical implications // Contrib. Mineral. Petrol., 56, 1976, p.119-134