Главная
страница 1
скачать файл
УДК 552.3

Новый алгоритм петрохимических пересчетов на основе CIPW нормы.



Д.Г.Егоров

Геологический институт Кольского НЦ РАН. 184200, Мурманская область, Апатиты, ул.Ферсмана, 14, Россия.



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

Петрохимические пересчеты анализов играют заметную роль как для целей петрографических классификаций, так и при изучении состава и генезиса горных пород. Особенно велика их роль при исследовании эффузивных пород, где прямая петрографическая диагностика часто весьма затруднительна. Одним из старейших, но до сих пор не утративших свое значение методов пересчета является разработанный В.Кроссом, Дж.Иддингсом, А.Пирсоном и Х.Вашингтоном метод CIPW, названный так по первым буквам фамилий авторов (Cross et al, 1903). Алгоритм в его классическом виде широко известен и описан в большом количестве справочной петрографической литературы (Ефремова, Стафеев, 1985; Четвериков, 1956), и до сих пор широко используется, особенно за рубежом. Популярность CIPW связана с тем, что положенная в основу метода ступенчатая десиликация действительно позволяет по появлению нормативных минералов отображать важные этапы петрохимической эволюции магматических серий. Тем не менее реально процесс десиликации в природе не существует, схема кристаллизации магмы скорее прямо ему противоположна: отдельные атомы в расплаве связываются друг с другом в элементарные кластеры, те в свою очередь - в более сложные, и т.д., до тех пор, пока для этого процесса силикации хватает кремнезема.

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

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

SiO2  кварц

Al2O3  корунд

Fe2O3  гематит

TiO2  рутил

K2O + кварц  кали-силикат

Na2O + кварц  натр-силикат

кали-силикат + корунд + кварц  калиофилит

натр-силикат + корунд + кварц  нефелин

натр-силикат + гематит + 3 кварц  эгирин

FeO+ рутил  ильменит

FeO+ гематит  магнетит

CaO + корунд + 2 кварц  aнортит

CaO+ рутил  перовскит

2CaO+ кварц  ларнит

2(FeO+MgO)+ кварц  оливин

калиофилит + 2 кварц  лейцит

ларнит +оливин + 2 кварц 2 диопсид

ларнит+ кварц  2 волластонит

лейцит + 2 кварц  ортоклаз

нефелин + 4 кварц  aльбит

перовскит + кварц  сфен

оливин + кварц  2 гиперстен

В сравнении с классическим алгоритмом, описание которого занимает несколько страниц, выигрыш в краткости очевиден (результат расчета идентичен классическому варианту CIPW). Для всех групп пород (насыщенных, недосыщенных и пересыщенных глиноземом) порядок расчета одинаков. На основе предлагаемого алгоритма чрезвычайно просто создавать компьютерные приложения, в том числе макросы электронных таблиц. Главное же достоинство предлагаемого алгоритма - его гибкость и наглядность, в силу чего он становится весьма удобным для внесения изменений (при работе с породами, химизм которых сильно отклоняется от принятого в CIPW стандарта). В частности, весьма просто можно менять последовательность расчета минералов в зависимости от новых экспериментальных и полевых данных. Так, как указывается М.И.Дубровским (1993), в соответствии с приведенными Гиорсо и Кармайклом (1992) данными по реакциям силицификации можно рекомендовать вычисление волластонита непосредственно после оливина, а сфена - после ортоклаза, что дает в большинстве случаев более точное соответствие нормативных расчетов петрологическим закономерностям. В этом случае алгоритм начиная с расчета оливина изменится:

2(FeO+MgO)+ кварц  оливин

ларнит+ кварц  2 волластонит

калиофилит + 2 кварц  лейцит

ларнит +оливин + 2 кварц 2 диопсид

лейцит + 2 кварц  ортоклаз

перовскит + кварц  сфен

нефелин + 4 кварц  aльбит

оливин + кварц  2 гиперстен

гиперстен+волластонит 2диопсид

В таблице приведены результаты расчета на миналы химических анализов недосыщенных кремнеземом горных пород, взятых из справочника (Ефимов, Стафеева, 1985), по классическому CIPW, и в альтернативном варианте. Видно, что в этом случае результаты расчета сильно зависят от принятого алгоритма.

Автор благодарит М.И.Дубровского за плодотворное обсуждение работы, и Российский фонд фундаментальных исследований за финансовую поддержку (проект 96-05-64402).

ЛИТЕРАТУРА

Гиорсо М.С., Кармайкл И.С. Моделирование магматических систем: петрологические приложения. // Термодинамическое моделирование в геологии. М.:Мир, 1992. Стр. 485-518.

Дубровский М.И. Физико-химические модели кристаллизации магматических оливин-нормативных пород нормальной щелочности. С-Пб.:Наука, 1993. 215 с.

Ефремова С.В., Стафеев К.Г. Петрохимические методы исследования горных пород. М: Недра, 1985. 511 с.

Четвериков С.Д. Руководство к петрохимическим пересчетам. М.,Госгеолтехиздат,1956. 246 c.

Cross W., Iddings J.P., Pirrson L.V. and Washington H.S. Quantitative classification of igneous rocks. University of Chicago Press. 1903.

Таблица. Результаты расчета на миналы химических анализов горных пород по классическому CIPW, и в альтернативном варианте. Все данные приведены в весовых процентах.






мельтейгит

мелилитолит




SiO2

41.9




36.1




TiO2

2.21




3.35




Al2O3

12.2




3.5




Fe2O3

6.41




7.07




FeO

4.32




6.41




MnO

0.22




0.09




MgO

5.45




9.98




CaO

16.6




31.36




Na2O

5.1




1.86




K2O

2.66




0.11




сумма:

97.07




99.83







CIPW

альтерн.

CIPW

альтерн.

ильменит

4.2

4.2

6.37

6.37

магнетит

8.24

8.24

10.26

10.26

гематит

0.73

0.73







ортоклаз




13.18







анортит

2.55

2.55

0.88

0.88

лейцит

12.32

1.99

0.51




нефелин

23.37

23.37

8.52

8.52

калиофиллит










0.37

волластонит

12.92

33.31




29.92

диопсид

29.27




26.92




геденбергит







0.53




ларнит

3.48




36.93

25.63

форстерит




9.51

8.66

17.41

фаялит







0.22

0.44



















УДК 552.3

Новый алгоритм петрохимических пересчетов на основе CIPW нормы.

Д.Г.Егоров (Геологический институт Кольского НЦ РАН. 184200, Мурманская область, Апатиты, ул.Ферсмана, 14, Россия.)



Предлагается новый алгоритм вычисления CIPW нормы, краткий, наглядный, удобный для внесения изменений при работе с щелочными породами, и любыми другими, схема кристаллизации которых отличается от принятого в CIPW стандарта. Библиогр. 3 наим., 1 табл.

The new algoritm of the CIPW norms calculation.
D.G. Yegorov (Geological Institute, Kola Scientific Center, Russian Academy of Sciences, 14, Fersman Street, Apatity, 184200, Russia).
The proposed original algorithm is calculated CIPW norms. The algorithm consist of certain steps directed to combine normative minerals starting from smaller molecules to bigger ones in order to "consume" step by step all silica presented in an analysis. It is assumed that proposed scheme more close follows natural sequence of a magma crystallization, than the conventional algorithm which prescripts the successive steps of a desilication of minerals presented by complex molecules. The algorithm is easy changed for the corresponding with the different schemes of the magma cristallisation.

“Петрология”: ссылка - Гиорсо, Кармайкл: M.S.Ghiorso, I.S.E.Carmichael. Thermodynamic modeling of geological materials: mineral, fluids and melts. // Reviews in Mineralogy, 1987. Vol. 17.
скачать файл



Смотрите также:
Новый алгоритм петрохимических пересчетов на основе cipw нормы
67.77kb.
Рассмотрена задача оценивания параметров регрессионных моделей. Для ее решения предлагается новый алгоритм, основанный на совместном использовании метода максимального правдоподобия и техники построения кривых Пирсона
19.5kb.
Строительные нормы
1301.64kb.
Тема: Ученые-лингвисты, занимающиеся проблемами культуры речи
29.43kb.
В статье разработан оригинальный алгоритм построения равновесного по Нэшу набора стратегий в игре нескольких лиц «в перемещения». Алгоритм реализован в виде действующей программы
74.89kb.
Лабораторная работа №3 на тему «Хэш функции: алгоритм n-хэш»
66.65kb.
Религиозная мораль
18.19kb.
Алгоритм определения астрономического азимута оптико-электронным астровизиром
108.87kb.
1 Настоящие нормы следует соблюдать при проектировании стальных строительных конструкций зданий и сооружений различного назначения
5254.67kb.
Лексические нормы русского языка
70.71kb.
«новый год в прямом эфире»
22.18kb.
Наркотик и организм человека
64.65kb.