Геобаротермометрия

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

Проверить нейтральность. На странице обсуждения должны быть подробности.

Стиль этой статьи неэнциклопедичен или нарушает нормы русского языка. Статью следует исправить согласно стилистическим правилам Википедии.

Геобаротермометрия (гео- земля; баро- давление; термо- тепло; метрия- измерение) совокупность методов определения температуры и давления образования минералов и горных пород.

Выделяются два основных направления:

1. Измерение температур (далее Т) и давлений (далее Р) образования современных геологических объектов, например, продуктов извержений вулканов. Главными методами являются физические методы измерения. К последним, например, относятся оптические или радиофизические исследования.

2. Определение древних (палео) Т и Р минералообразования в ранее образованных породах и минералах различного происхождения. При геологических исследованиях является важным методом исследования.

Все методы подразделяются на прямые и косвенные.

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

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

Методы геобаротермометрии

Аналоговые методы

Эта группа включает методы определения Т — Р параметров, основанные на сравнении с ранее известными или определёнными экспериментами данными. Для решения этой задачи используются следующие методы^[1]:

1. Температуры плавления минералов и горных пород;
2. Температуры разложения;
3. Температуры фазовых превращений;
4. Температуры распада твёрдых растворов;
5. Наличие типоморфных минералов;
6. Парагенетические ассоциации минералов;
7. Данные о Т выделения минералов из растворов или расплавов.

Некоторые примеры приведены в таблице:

Формула	Минерал	Tплав°С	Tразл°С	Tфаз.перех°С
Pb	Свинец	327,6
MgO	Периклаз	2800
SiO2	Кварц	1370		573, 867
Fe2SiO4	Фаялит	1205
BaSO4	Барит		01580
FeCO3	Сидерит		0282
ZnSO4	Цинкозит		0740
AsS	Реальгар			000267

Эти методы дают недостаточно точные, порой дают ложные результаты. Например, в природных условиях кварц выделяется в широком интервале Т (от 100 до 1000 °C), тогда как его Т плавления превышает 1000 °C.

Физические методы

Можно выделить две группы методов^[1]:

А.Собственно физические методы включают

1. Методы с применением термопар. Этим методом пользуются при измерении температур жидких лав. Одним из приборов для измерения является платиново-иридиевая термопара.
2. Дистанционные оптические методы на основе тесной связи температуры среды с её цветом. Для этого используются различные оптические пирометры.

Б.Они включают методы анализа газово-жидких включений (далее ГЖВ)^[2],^[3], декрепитации^[4] и расплавных включений (далее РВ)^[5]. В СССР основателем этого метода был профессор МГУ Н. П. Ермаков. Свою первую работу он опубликовал ещё в 1949 г.^[6]. Ранее эти исследования назывались «минералогическая термометрия» [2]. В настоящее время результаты этих исследований объединяются общим понятием — термобарогеохимия.

«Объектом её изучения являются разнообразные по составу и агрегатному состоянию флюидные включения¸ широко распространенные в минералах пневматолитового и гидротермального происхождения и встречающиеся в минералах интрузивных и эффузивных пород. Эти мельчайшие остатки минералообразующей среды несут информацию о процессах¸ происходивших при становлении магматических тел и формировании эндогенных рудных месторождений. Методы термобарогеохимии позволяют определять относительную и истинную температуры образования минералов, количественный и качественный состав растворов и расплавов во включениях¸ давление и агрегатное состояние минералообразующей среды¸ из которой происходила кристаллизация минералов или рудоотложение»^[7].

В результате проведения этих исследований были получены не только достаточно полные характеристики многих природных геологических процессов (магматизма, метаморфизма, гидротермальных образований и пр.), но и данные, не соответствующие общепринятым представлениям о минералообразовании. Например, по А. Н. Заварицкому (1884—1952) [3], [4]^[8] в гранитах кварц (далее Qw) выделяется последним, заполняя промежутки между ранее выделившимися минералами. Но по результатам анализа РВ^[9] и изотопов кислорода температура образования (далее Т_обр) Qw около 700 °C и близка T_обр наиболее раннего минерала — биотита (далее Bio). Таким образом Qw является наиболее ранним минералом.

Термодинамические методы

В эту группу (далее ТДМ) включены изотопные и геохимические геобаротермометры. Они основаны на анализе особенностей распределения элементов и их изотопов между соединениями (минералами) под влиянием внешних термодинамических, как правило температурных и барических, условий.

Изотопные геотермометры

Теория метода описана в многочисленных работах, прежде всего^[10],^[11] и т. д. В советской и русской литературах эти вопросы подробно освещены в работах^[12] и^[13] и др. В его основе лежат реакции изотопного обмена, по смыслу идентичные реакциям ионного обмена. Достоинством их является независимость от давления Р показаний ИГТ в геологически востребованном интервале глубин. Например, для кальцита установлено значение всего 0,03 %о/кбар^[14]. Это обусловлено тем, что размерные параметры кристаллических решёток, содержащих изотопно различные компоненты, отличаются друг от друга очень слабо. Так, расстояние Na — Cl в молекулах Na³⁵Cl и Na³⁷Cl равно соответственно 4,485 и 4,4601Å^[15], между компонентами в молекулах ¹⁰⁷AgBr и ¹⁰⁹AgBr — 2,393068 и 2,393074Å, в соединениях ¹⁰⁷Ag³⁵Cl и ¹⁰⁹Ag³⁷Cl — 2,28074 и 2,28082Å^[16]. В протоэтилене С₂Н₄ расстояние С — Н равно 1,1030Å, а в дейтероэтилене C₂D₄ расстояние С — D равно 1,089Å; для Н₂O параметры решётки равны: a_o = 4,5134Å, с_o = 7,3521Å, для D₂O — соответственно 4,5165Å и 7,3537Å^[17]. Малые изменения межатомных расстояний ведут к слабому изменению объёма кристаллической решётки. Например, объём кристаллической решётки оксалатов NaНСО₂С₄Н₂О и NaDC₂О₄D₂О изменяется от 224,429 до 225,113 (Å)³^[18] и свидетельствуют об относительной независимости изотопных распределений от давления на доступных для геологических наблюдений глубинах, что подтверждено экспериментами с изотопами О в системе СаСО₃— Н₂О^[14] и С в системе СН₄ — Н₂О^[19]. Согласно этим исследованиям распределение изотопов между двумя соединениями (минералами) подчиняется уравнению ${1000\ln\alpha = {B\over {T^2}} + A}\,$, где $\alpha\,$ — показатель (коэффициент) фракционирования, А и В — постоянные; физический смысл их обычно не ясен. Это уравнение строится по теоретическим или экспериментальным данным. В настоящее время температурные исследования проводятся преимущественно с изотопами кислорода, водорода, углерода, серы.

Наиболее успешно ИГТ используются при работе с изотермами, то есть прямыми линиями, построенными по множеству проб, в которых определены величины $\delta\,$ — индексов для двух минералов, образованных при одной и той же Т°С. Тогда угловой коэффициент $S\,$ изотермы определяется выражением

$S= \frac{ln\alpha(M_{Y})}{ln\alpha(M_{X})}\qquad{(1)}\,$,

где $X\,$ и $Y\,$ — оси координат, по которым откладываются значения $\delta\,$ — индексов для каждого изучаемого минерала. В этом случае выделяются несколько классов ИГТ^[20].

Мономинеральные ИГТ

Они основаны на анализе распределений изотопов разных элементов в одном и том же минерале (карбонаты (углерод — кислород), мусковит и биотит (водород — кислород), сульфаты (сера — кислород) и пр.) по уравнению (1)^[20].

Двуминеральные ИГТ

Моноэлементные ИГТ с одновременно образованными минералами^[20].

Этот ИГТ- единственная в настоящее время форма, применяемая при изотопных исследованиях. Их теория хорошо разработана. Примером является ИГТ на паре Qw- Mt (магнетит), которые обмениваются между собой изотопами кислорода по реакции

Qw(¹⁸О) + Mt(¹⁶О) ↔ Qw(¹⁶О) + Mt(¹⁸О). ……………………………………………………………………(2)

Для этой реакции выражение ${\delta{^{18}}O}\,$(Qw)-${\delta{^{18}}O}\,$(Mt) = $f(T)\,$ = ${1000\ln\alpha}\,$, используется для определения Т_обр минералов. Это выражение следует из уравнения (1), если S= 1 и C₁= C₂. C₁ и C₂ — соединения, изотопно равновесные исследуемым минералам. Здесь ${\delta{^{18}}O}\,$ ($\delta\,$ — индекс) отражает отношение ${\frac{^{18}O}{^{16}O}}\,$ в минерале, отнесённое к эталону (для кислорода среднее отношение концентраций изотопов в океанической воде -SMOW (Standart Mode of Ocean Water)). Размерность ${\delta{^{18}}O}\,$ — промиле ‰.

В целом, этот вид геотермометра не предназначен для определения Т. Он должен использоваться только как база для построения других ИГТ, то есть имеет феноменологическое значение.

Моно- и биэлементные ИГТ с минералами в произвольных временных соотношениях^[20].

Это разновременно образованные соединения с изотопами как общих, так и различных элементов (например, окисел и сульфид, кислород Qw и водород Bio), причем Т_обр одного из них (эталона) и состав соединения C для него заранее известны. Имеются разновидности.

1.Моноэлементные ИГТ .

Силикаты магматитов и метаморфитов образуют пары минералов с изотопами общего элемента. Например, широко распространённые пары минералов в гранитах Qw — Fs (полевые шпаты).

2. Биэлементные ИГТ (например, ${\delta{^{18}}O}\,$(Qw) и ${\delta{D}}\,$(Bio)).

Наличие нескольких минералов, взятых из одних и тех же проб, позволяет провести контрольную оценку полученных результатов. Например, по данным рядовых определений по парам кварц — мусковит и кварц — калиевый полевой шпат контрольной может быть пара мусковит— калиевый полевой шпат.

Использование этих классов ИГТ позволило определить Т_обр широкого круга конкретных минералов. Некоторые примеры из^[21] приведены в таблице.

Породы	Регионы	Кварц	Биотит	Ильменит	Магнетит	Калишпат	Мусковит	Альбит	Гранат
Пегматит	Калифорния	780	-	-	-	-	490	510	-
Пегматит	США	700	-	-	-	-	-	-	320
Пегматит	США	-	-	680	-	-	-	480	-
Пегматит	Норвегия	760	770	-	660	-	480	440	260
Аплит	Калифорния	710	-	-	-	580	570	-	-
Аплит	Карелия	710	-	-	-	-	-	440	-
Гранит	Кавказ	710	-	-	584	-	-	-	-
Гранодиорит	Япония	700	-	-	-	570	-	-	-
Сланцы	Австрия	700	-	-	-	-	-	-	330
Сланцы	Гренландия	700	-	-	610	-	-	-	-
Метапелит	Альпы	670	-	604	-	-	-	-	-
Ортогнейс	Альпы	650	-	680	-	550	-	-	-
Гнейс	Альпы	700	-	-	-	-	-	-	320

В результате установлено, что Т_обр одноимённых минералов в кислых магматических и метаморфических разных ступеней метаморфизма породах оказались равными. Порядок их выделения отражён в последовательности [5]

(Bio,Qw)(≈750 °C)>Il(≈650)>Mt(600—650)>Kf(≈570)>Mus(450—570)>Al(440—500)>Grn(≈300 °C).

В этой последовательности Il — ильменит, Kf — калиевый полевой шпат; Mus — мусковит; Al — альбит;Grn — гранат. Во всех случаях все силикатные минералы выделяются в изотопном равновесии с водой. Рудные минералы — Mt и Il — выделяются в изотопном равновесии с рутилом (?).

В оливинах, клинопироксенах, плагиоклазах основных и ультраосновных горных пород и метеоритов для Т ≈ 1250 °C выявлено изотопное равновесие по кислороду с углекислым газом (СО₂)[6]

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

Геохимические геобаротермометры

Геохимические геобаротермометры (далее ГГБТ). Термодинамическая база этих методов разработана достаточно полно и учитывает, в частности, неидеальность твёрдых растворов минеральных систем (J. Ganguly, S.K. Saxena^[22] и др.). Её обобщение проведено В. А. Курепиным^[23]. В России эти методы применяются достаточно давно; их родоначальником является профессор МГУ Л. Л. Перчук [7]. Эти методы он называл «геологической термобарометрией». Выделяются классы ГГБТ:

1. геохимические геотермометры.
2. геохимические геобарометры.

На практике различие между ними определяется чувствительностью их к соответствующим параметрам среды- Т и Р. На самом деле различие между ними более существенное: они различаются типом протекающих химических реакций, положенных в основу ГГБТ и ведущих к образованию минералов.

Любая химическая реакция характеризуется определёнными термодинамическими параметрами, в частности, изобарно- изотермическим потенциалом, или свободной энергией Гиббса, $\Delta{G}\,$. Эта величина складывается из трёх составляющих $\Delta{G}= {\Delta{G}_T+\Delta{G}_P+ \Delta{G}_{mix}}\,$, где K — коэффициент распределения элементов между минералами; $\Delta{G}_{T}\,$ —тепловая (основная) составляющая , $\Delta{G}_{P}\,$ -барическая составляющая, то есть связанная с влиянием давления, и $\Delta{G}_{mix}\,$ — энергия смешения компонентов при образовании минерала. Все эти параметры рассчитываются определённым образом: ${\Delta{G}_T = -RT\ln{K}_{react} = AR + BRT}\,$; ${\Delta{G}_P = (P-1)\Delta{V}_{react}}\,$; $\Delta{G}_{mix}\,$ рассчитывается на основе экспериментальных исследований и зависит от типа твёрдого раствора. Здесь R — газовая постоянная; ${\Delta{V}_{react}}\,$ — объёмный эффект реакции. Как и в изотопных геотермометрах физические смыслы параметров ${A}\,$ и ${B}\,$ строго не определёны.

Геохимические геотермометры

Геохимические геотермометры (далее ГГТ) основаны на реакциях ионного обмена элементов в разных средах и поэтому по своей сущности близки ИГТ. Типичным примером является хорошо изученный биотит (Bio) − гранатовый (Gr) геотермометр^[24], основанный на ионом обмене между минералами магния и железа Fe⁺²:

Bio(Fe) + Gr(Mg) ↔ Bio(Mg) + Gr(Fe)………………………………………………………………………..(3)

Для этой реакции определяется величина $ln{K} = \frac{A}{T} + B$ , причём $RTln{K} = -\Delta{G}\,$, $\Delta{G}\,$ — свободная энергия Гиббса реакции ионного обмена. В этом выражении $~K$ — коэффициент распределения атомов Fe⁺² и Mg между биотитом и гранатом

$K = {\frac{\left [{\frac{a^{Fe}_{Gr}}{a^{Mg}_{Bio}}}\right]}{\left[{\frac{a^{Fe}_{Bio}}{a^{Mg}_{Gr}}}\right]}}\,\qquad {(4)}$

${a^{Fe}_{Gr}}$ — коэффициенты активности соответствующих элементов (Fe или Mg) в минерале (Gr или Bio). Практически величины активностей заменяются значениями мольных долей элементов в минерале, то есть

$K = {\frac{\left({\frac{M^{Fe}}{M^{Mg}}}\right)_{Bio}}{\left({\frac{M^{Fe}}{M^{Mg}}}\right)_{Gr}}}\,\qquad{(5)}$

$M^{Fe}\,$ — мольные доли элементов в минерале. В такой записи вместо мольных долей используются значения кристаллохимических коэффициентов. В этих реакциях происходит замещение одного на другой двух изоморфных элементов, для которых размеры атомов достаточно близки друг другу. Из-за этого и размерные параметры минералов, которые образуют эти элементы, различаются незначительно. На практике можно пренебречь и влиянием неидеальности смешения соответствующих твёрдых растворов. Что касается других элементов с неидеальной растворимостью, то в практических исследованиях для изоморфных элементов влиянием неидеальности образуемых ими твёрдыми растворами также можно пренебречь. Вклад параметров, учитывающий неидеальность твёрдых растворов, соизмерим с ошибкой определения термодинамических параметров. Например, по разным источникам значение ${-\Delta{G}_{o}}\,$ (298°К) гроссуляра изменяется в пределах 1488,696- 1509,0 ккал/М (разница (Δ) почти 20 ккал/М, или ≈ 1,34 % от минимального значения); S_o(298К)- 56,7- 76,97 кал/М°К. Для андрадита ${-\Delta{H}_{o}}\,$(298°К)= 1373,5- 1410,3 ккал/М (Δ = 27 ккал/М, или ≈ 1,97 %). Параметры W же Маргулиса имеют значения: в андрадитах W(Mg- Ca)= 3,83; W(Mg- (Fe+ Mn))=3,0; W(Ca-(Fe+ Mn))= 1,00 ккал/М. Так как параметр W используется в сочетании с мольными долями X, которые всегда <1, или X², с ростом Т величина |ΔG| растет, а W- уменьшается, то этот вклад будет ещё меньше.

В приведённом примере использовано распределения двух элементов в двух минералах. Однако имеются разработки, в которых изучается распределение одного элемента между двумя средами, например, распределение Ni между расплавом и кристаллизующимся из него оливином^[25]. Многочисленные примеры показывают лучшую устойчивость ГГТ при использовании системы из двух минералов и двух элементов.

Геохимические геобарометры

Геохимические геобарометры (далее ГГБ) основаны на анализе химических уравнений образования минералов, протекающих в природных условиях. Примеры уравнений геобарометров (Перчук,1977^[26]):

1. Mg₃Al₂Si₃O₁₂ + 3(1-x)MgSiO₃ = 3MgSiO₃xAl₂O₃;
2. MgAl₂O₄+ (1+x)MgSiO₃ = MgSiO₃xAl₂O₃ + xMgSiO₄.

Как видим, они описывают реакции образования новых минералов. Важным свойством этих уравнений является высокое значение объёмного эффекта $\Delta{V}_{react}\,$ реакции. Например, для второй реакции при х=1 $\Delta{V}_{react}\,$ = −18,03 см³/М (или ≈ −67 % относительно оливина (мольный объём $\Delta{V}_{ol}\,$ = 26,89 см³/М — минимальное значение для данного набора минералов). Для биотит — гранатового геотермометра (реакция 1) $\Delta{V}_{react}\,$ = 7,6 см³/М (или ≈6,7 % относительн пиропа (pyr) $\Delta{V}_{pyr}\,$ = 113,3 см³/М). Эти значения даже по абсолютной величине значительно различаются. Если их взять относительно конкретного минерала, то различие будет ещё больше. Это приводит к существенной зависимости направления реакции от величины внешнего давления

Влияние диффузии

Широкое участие диффузии при образовании минералов и влияние её на значения Т_из обусловлено наличием газовых составляющих в качестве компонентов изотопных систем. Как правило это вода, CO₂, реже СН₄.Влияние диффузии на изотопный состав природных соединений оценивается через показатель диффузионного фракционирования ${\alpha_{d}}\,$ по схеме: ${D\rightarrowtail{\alpha_{d}}\rightarrowtail{\alpha_{d}^{n}}}\,$. Здесь ${\alpha_{d} = {\sqrt{\frac {D}{D_T}}} = {\sqrt{\frac {M_T}{M}}}}\,$; здесь D и D_T -коэффициенты диффузии, М и М_T — молекулярные веса легкой и тяжелой компонент; n- число теоретических тарелок, характеризующих глубину диффузионного фракционирования (разделения)^[27],^[28]. При одновременном действии двух механизмов фракционирования с показателями $\alpha_{T}\,$ и $\alpha_{d}^{n}\,$ суммарный эффект разделения определится выражением $\alpha={\alpha_{T}}\cdot{\alpha_{d}^{n}}\,$. Если диффундируют изотопы одного элемента, то приходим к равенству (6),

$S = \frac{1000ln\alpha_{T_{1}} + n1000ln\alpha_{d_{1}}}{1000ln\alpha_{T_{2}}}\qquad{(6)}$

показывающему зависимость Т_из от диффузии. Методика учета влияния диффузии описана в^[29]. Если диффундантом является вода, то

$\alpha_{d} = \sqrt{19/18} = 1{,}027$

При формировании карбонатных систем (скарнов, мраморов, не магматических карбонатитов) величина n изменяется в пределах 0,03- 3,0 для H₂O, CH₄, CO₂. Полностью отсутствует диффузия в магматических карбонатитах. Для компоненты HDO в мусковитах аплитов, сланцев и гранитогнейсов n= 0-1,0; в биотитах аляскитов, гранодиоритов, сланцах, гнейсах- 0-0,65; в роговых обманках гранодиоритов n=1,5.

Влияние диффузии выявлено при анализе изотопного состава кислорода в плагиоклазах, пироксенах и оливинах мафических пород^[28]. Величина глубины фракционирования в молекуле СО₂ составляет 0,19 (то есть n = 0,19).

Результаты определения влияния диффузии в минералах на Т_из приведены в таблице^[21]. В таблице And — андалузит, Sil — силлиманит,Stv — ставролит.

№ п.п.	Минерал	Порода	Регион	Тиз	n	Источник анализов
00001	Роговая обманка	0Гранодиорит	0Невада	553	1,5	Taylor B.E. 1977
00002	000Мусковит	00Пегматит	Калифорния	476	1,0	Taylor B.E et. al. 1968
00003	000Мусковит	000Аплит	Калифорния	470	1,0	Taylor B.E et. al. 1968
00004	000Мусковит	0Сланец And	0Невада	440	1,0	Shich T.N.et.al.1969
00005	000Мусковит	000Сланец	0Греция	436	1,0	Rye R.O. et.at.,1976
00006	0000Биотит	0Гранодиорит	0Невада	760	0,65	Shappard S.M.,1974
00007	0000Биотит	Сланец Sil-Mu	0Намибия	746	0,65	Hoernes S.et.al.1979
00008	0000Биотит	Сланец Stv-And	00Альпы	694	0,65	Hoernes S.et.al.1978
00009	0000Биотит	00Гнейс Bio-Pl	00Альпы	756	0,65	Бутвин В. В., 1984
00010	0000Биотит	00Ортогнейс	00Альпы	667	0,65	Hoernes S.et.al.1978

Основные особенности применения термодинамических геотермометров

Для успешного использования ИГТ и ГГТ необходимо выполнение следующих условий^[20],^[30]:

1. минералы должны быть одновременны;
2. Т_обр минералов должны быть равными;
3. минералы должны контактировать друг с другом;
4. между минералами должен протекать изотопный или ионный обмен; приведённые выше положения способствуют протеканию этого обмена;
5. в результате протекания этой реакции должно установиться термодинамические изотопное или геохимическое равновесия.

Практика показала, что эти критерии даже не проверяются. Она выявила недостатки применения ГГБТ, не подвергнутые, однако, анализу. Многие авторы, например,^[31],^[32],^[33] и др., отмечают частое несовпадение Т_обр минералов, полученных разными авторами и ИГТ (или ГГТ). Якобы для устранения указанных проблем широко распространен ввод в уравнения ГГТ концентраций третьих элементов^[34],^[35]. Это приводит к громоздкости выражений, описывающих ГГТ; так В. И. Ваганов и С. В. Соколов^[32] упоминают выражение ГГТ, содержащее около 24 числовых параметров. В результате с одной стороны для расчетов приходится использовать ЭВМ без серьёзной практической отдачи. С другой стороны точность этих операций оставляет желать лучшего и фактически не изменяется. В связи с этим ряд исследователей, например, B. Wood (1974)^[32], T. Morry (1975)^[32], отказывается от введения таких поправок. Нами уже отмечались многие из этих недостатков^[36],^[37], Эти явления привели также к появлению гипотез о «неравновесном состоянии» при выделении минералов, а также существенном влиянии скорости протекания реакции образования минералов на распределения изотопов или вообще элементов^[38].

Не отрицая влияния указанных факторов, заметим наличие в технологии применения этих методов некоторых особенностей:

1. Изучаются только точечные пробы, то-есть берутся два минерала из одной точки.
2. Все минералы отбираются из валовой выборки, в которой между изучаемыми минералами присутствуют третьи минералы.
3. Одновременность выделения минералов не ясна.
4. Отсутствуют доказательства образования минералов при одних и тех же температурах или давлениях.
5. Наличие термодинамического изотопного и геохимического обмена элементами не выявляется.
6. Существование термодинамических изотопного и геохимического равновесий между минералами под вопросом.

Для широко распространённого Qw — Mt — ИГТ по данным классических представлений Т_обр Qw около 300—400°С, а Mt — 600—700°С, то есть они образованы при разных температурах. По классическим представлениям Qw является наиболее поздним минералом, тогда как акцессорный Mt — один из наиболее ранних минералов. В метаморфических породах температурные и временные соотношения между этими минералами вообще не ясны. Так в ГГТ (2) по изотопным данным Bio образован при ≈700 °C, Gr же по изотопным данным и данным анализа ГЖВ образован при Т ≈ 300—400 °С^[20]. Такая же ситуация с амфибол — гранатовым ГГТ . На время создания ГГТ считалось, что Gr образован при ≈ 700 °C, тогда как Т_обр амфибола намного меньше, то есть эти минералы не одновременны.

7. Выявляются свойства, заставляющие сомневаться в достоверности ГГБ. Например, не ясно, по какому принципу выбираются минералы как с точки зрения их совместного нахождения в породе, так и с точки зрения возможностей применения к ним принципов термодинамики. Ясно одно, что на данный момент, во-первых, доказательств существования этих реакций в реальных природных условиях нет; во-вторых, все эти уравнения составляются произвольно и не учитываются реальные условия образования минералов. Многие ГГТ составляются на основе экспериментальных исследований, в которых не изучен механизм образования минералов.

8. Ошибки в написании формулы геотермометра, имеющей вид

$lnK = {-A\over {T}}+ B$.

8a. В этом уравнении не выявлен физический смысл параметров $~A$ и $~B$. При теоретическом анализе это уравнение обычно путают с равенством

$ln{K}_T = {-{A}_T\over {T}}+ {B}_T$,

в котором параметры $~{A}_T$ и $~{B}_T$ зависят от температуры и для каждого конкретного значения T имеют вполне определённую величину. В уравнении же эти параметры постоянны и не зависят от T.

8б. Не ясен физический смысл параметра $\ln{K}\,$. Это выражение имеет смысл, если величина $K$ является безразмерной величиной. Анализ^[39], показал, что это не так, поэтому выражение $\ln{K}\,$ не имеет физического смысла. Эта ошибка восходит к ошибке в написании химического потенциала ${\mu= \mu_{0} + RT\ln{C}}\,$, где $C\,$ -концентрация некоторого компонента в растворе. Поскольку параметр $C\,$ является размерной величиной, то выражение $\ln{C}\,$ также не имеет физического смысла. На самом деле это уравнение должно иметь вид

$\mu= \mu_{0} + RT\ln{\frac{C}{C_0}}\,$ .

Размерность далее обозначается через квадратные скобки […], то есть размерность $X_i\,$ = $\left[X_i\right]\,$. В термодинамике имеется достаточно много выражений, в которых размерность величин не учитывается. Например, закон Рауля в форме $P_i = {X_i}\cdot{P_0}\,$, $P_0\,$ и $P_i\,$ — давление пара над чистым растворителем и раствором, содержащим $X_i\,$ мольных долей растворенного вещества. Поскольку $\left[X_i\right]\,$ ≠ 0, то равенство физически не корректно. Другой пример — абсолютная энтропия в форме

$S = {{\frac{R}{{\gamma}-1}}\cdot{\ln{P{V}^{\gamma}}} + C}$

где $C\,$ — постоянная интегрирования^[40]; а) поскольку $\left[V\right]\,$≠ 0, то не определены размерности $\left[V^{\gamma}\right]\,$= ? и $\left[{P}\cdot{V^{\gamma}}\right]\,$= ?; б) так как $\left[{P}\cdot{V^{\gamma}}\right]\,$= ?, то остаются не ясными физический смысл и размерность $\ln{\left({{P}\cdot{V^{\gamma}}}\right)}\,$. То же относится и к выражению ^[41],

$S = {{C_V}\cdot{\ln{T} + R{\ln{V}}} + Const}\,$

поскольку $\left[\ln{T}\right]\,$ = ? и $\left[\ln{V}\right]\,$= ? (V- объём) (правильно $\ln{\frac{T}{T_0}}\,$ или $\ln{\frac{V}{V_0}}\,$).

Реальное выражение для уравнения геотермометра исходит из представлений^[30]

$R{\ln{\frac{K}{K_{pe}}}} = {-}{\Delta{H_{pe}}}{\left[{{\frac{1}{T}} -{\frac{1}{T_{pe}}}}\right]}\qquad{(7)}$,

где $K_{pe}$, ${\Delta{H}_{pe}}\,$ и $T_{pe}\,$ — значения коэффициента распределения, энтальпии реакции и температуры в точке фазового перехода (phaze exchange = pe) (или в любой другой точке с известными названными параметрами). Довольно часто за точку фазового перехода берётся температура плавления минерала. Следовательно,

$A = \frac{\Delta{H}_{pe}}{R}$,

$B = \frac{\Delta{H}_{pe}}{R{T_{pe}}} = \Delta{S}_{pe}$

$\Delta{S}_{pe}\,$ — энтропия реакции обмена в точке фазового перехода.

9. Место геотермометров в ряду методов геологических исследований чётко не определёно, отсюда не чётко определены цели и задачи температурных исследований, а также взаимоотношений между аналоговыми и физическими методами определения температуры с одной стороны и термодинамическими методами — с другой^[20],^[42].

Область применения геобаротермометров

Постановка проблемы.

Основной задачей использования ТДМ является не просто определение Т_обр и Р_обр минералов, а через эту операцию и сопоставление полученных значений Т и P с таковыми, установленными независимыми аналоговыми или физическими методами, определение равновесного минералу соединения, а затем на этой основе — выявление механизма образования минералов, то есть вида химической реакции, при осуществлении которой происходит образование минерала^[43]. Решение этой задачи опирается на аксиомы:

• 1) все минералы выделяются в термодинамически изотопном и геохимическом равновесиях с некоторыми «компонентами» окружающей среды;
• 2) минералы и эти «компоненты» выделяются при разложении некоего протовещества (материнского соединения). Именно в процессе этого разложения достигаются необходимые термодинамические изотопное и геохимическое равновесия.

Намечается последовательность решения задачи:

1) с помощью ФМ определяется Т_обр и Р_обр минерала;

2) с помощью ИГТ и ГГТ определяются соответствующие Т_обр и Р_обр минерала и по сходству их с Т_обр, выявленного ФМ, определяется состав равновесного минералу соединения;

3) по этим данным составляется уравнение выделения минерала и с помощью ГБ определяется кажущееся Р_обр;

4) кажущееся Р_обр сопоставляется с Р_обр, выявленным ФМ, и по близости их друг другу устанавливается достоверность предложенного уравнения, на основе чего это уравнение принимается или отвергается.

Примеры решения задачи о механизме образования минералов

Карбонатная минерализация

Карбонатная минерализация представлена простыми (кальцит и пр.) или двойными (доломит и др.) солями угольной кислоты с наибольшим распространением кальциевых солей, прежде всего кальцит (КЛ), меньше доломит и арагонита. Кальциты и арагониты могут быть органогенными, являясь продуктом жизнедеятельности организмов (кораллы, моллюски и пр.), и неорганической природы (сталактиты, сталагмиты, травертины и пр.). Они существуют при достаточно низких и стабильных температурах (0- 40 °C). Громадная масса кальцитов слагает высокотемпературные образования: скарны, мрамора, гидротермальные жилы, кальцифиры, карбонатиты и пр.

Низкотемпературные кальциты.

В большинстве случаев предполагается выделение кальцитов по уравнению, неоднократно воспроизводимому в экспериментах,

Ca⁺² + 2HCO₃⁻¹ = CaCO₃ + H₂O + СО₂…………………………………………………………………..(8).

Однако в природных условиях механизм выделения не установлен. Имеются и экзотические схемы, например, редукция гипса или серы. Изотопный анализ выявил следующие схемы изотопных равновесий (запись означает: на первом месте стоит соединение, равновесное кальциту по углероду; на втором месте — по кислороду). Далее (*) отмечает обменивающийся изотоп:

1). CaCO₃ — ^*СО₂- Н₂О^*;

соответствует реакции (8). Возможно участие организмов, играющих роль катализаторов.

2). CaCO₃- ^*СН₄- Н₂О^*;

образуется по схеме

Са⁺²+ 2НСО₃⁻¹+ 4Н₂(орг) → СаСО₃+ ^*СН₄+ 3Н₂О^*(?);

под влиянием организмов, видимо, происходит преобразование аниона НСО₃⁻¹.

3). CaCO₃- ^*СН₂О- СН₂О^*;

характеризуется связью с растительной клеткой

2Ca⁺²+ 4НСО₃⁻¹+ 4Н₂ (орг) → 2CaCO₃+ ^*CH₂O+ CH₂O^*+ 4Н₂O,

хотя механизм возникновения равновесия между КЛ и клеткой не ясен.

4)CaCO₃- ^*СН₄- СО₂^*.

Диагенетические КЛ торфяных зон^[44] — по реакции

2Са⁺² + 4СН₃СОО⁻¹ + 4Н₂ орг) → 2СаСО₃ + 5^*СН₄ + СО₂^*

за счёт микробиологического разложения органических кислот. Практически во всех случаях на формирование изотопного состава кальцитов влияет диффузия преимущественно Н₂O и СО₂.

Высокотемпературные кальциты.

Все гидротермальные кальциты различных месторождений и регионов равновесны изотопной системе

CaCO₃ — ^*СО₂- Н₂О^*

и образованы по реакции (8), осложненной диффузией этих компонентов.

В метаморфитах с известной изначально известняковой природой (мрамора и скарны) по кальциту до 500 °C (скарны Невады и др.) типична изотопная система

СаСО₃-^*СО₂-СО₂^*.

При Т > 500 °C по кальциту(мрамора Невады, Адиронак, Украины, скарны Алдана)- изотопная система

СаСО₃-^*СН₄-СО₂^*.

Кальцифиры Якутии. Им соответствует изотопная система

СаСО₃ — ^*СН₄ — Н₂О^*

В кимберлитах Якутия выявлены изотопные системы:

СаСО₃-^*СО₂-СО₂^* и

СаСО₃-^*СН₄-СО₂^*.

Карбонатиты.

Фации по В. С. Самойлову^[45] хлорит-серицит- анкеритовая (по кальциту Т_из ≈ 250 °C) и амфибол- доломит- кальцитовая (по кальциту Т_из ≈ 288 °C) изотопная система

СаСО₃-^*СО₂- Н₂О^*;

Фация альбит — кальцитовая (по кальциту Т_из = 400—550°С) — в основном изотопная система

СаСО₃-^*СО₂-СО₂^*.

Фация калиево-полевошпат — кальцитовая (по кальциту Т_из = 500—750°С) — изотопная система

СаСО₃-^*СН₄-СО₂^*,

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

Силикатная минерализация

I. В этом разделе рассматриваются минералы кислых магматических и метаморфических пород (кварц, биотит и мусковит, плагиоклазы, калиевые полевые шпаты, гранаты, пироксены). Все эти минералы изучены изотопными методами; они выделяются в равновесии с водой^[21],^[43] за счёт образования при разложении гидратированных комплексов по схеме

(⁺³Si − O − Si⁺³)⁺⁶_распл → 2(H − O − Si⁺³є)_распл^[46].

В окрестности точки Т_крист за счёт полимеризации сиботаксической группы (S⁺³ − O − H)_распл образуется ассоциат H₄SiO₄, разлагаемый по схеме

H₄SiO₄ → SiO₂ + 2H₂O.

Этот механизм подтверждается образованием полевых шпатов (калишпат − альбит)(Reesmon, Keller^[47] ; Рыженко и др.^[48]; Алексеев и др.^[49])

Al(OH)₄⁻ + 3H₄SiO₄ + Na⁺ → [NaAlSi₃O₈ + 8H₂O].

В правой части уравнения присутствуют пары Al − H₂O или Kf − H₂O, соответствуя изотопным данным. Это уравнение нельзя признать точным, поскольку не ясна ситуация с изотопом водорода в воде. Если же привлечь гипотезу И. Г. Ганеева^[50], заменив H₄SiO₄ на Si(OH)₆⁻², то реакция перепишется в приемлемом виде

Al⁺³ + 3Si(OH)₆⁻² + Na⁺¹ + 2Н⁺¹ → [NaAlSi₃O₈ + 8H₂O] + 2OH⁻¹.

Для гранатов имеются и изотопные (кислород) и силикатные анализы. Изученные гранаты находятся в ассоциации преимущественно с биотитом, кордиеритом и плагиоклазом. По изотопным данным гранат выделяется при 300—450°С; результаты анализа ГЖВ дают те же пределы. Биотит и гранат выделяются в равновесии с водой. Установлено, что наилучше согласие наблюдается в равновесии граната с кордиеритом Cor, биотитом и пироксенами. По экспериментам Л. Л. Перчука и др.^[51] при Т= 550—1000°С при совместной кристаллизации ионный обмен между гранатом и кордиеритом отсутствует. По ряду признаков вероятное уравнение образования гранатов имеет вид

… = {Cor + [Grn }+ H₂O]+ ….

Здесь скобки отражают: […] — изотопное; {…} — геохимическое равновесия. Этот вид соответствует полной реакции:

хлорит + кварц → {кордиерит + [гранат} + вода].

II. Данные по анализам изотопов кислорода в минералах из метеоритов и пород (эклогиты, кимберлиты, лерцолиты, толеиты, щелочные базальты, основные и ультросновные породы, гавайяиты, габбро, адамелиты). Из анализа расплавных включений следует, что плагиоклазы (Pl) выделялись при 1100 −1300 °C; клинопироксены (Cpx) (преимущественно диопсиды) — 1100—1300°С; оливины (Ol) −1100- 1300 °C (^[52];^[53] и др.). Эти минералы изотопно равновесны СО₂. На распределение изотопов кислорода в минералах метеоритов влияет диффузия СО₂. Таким образом, поставщиком ¹⁸О в минералы, вероятно, является СО₂, учитывая хорошую растворимость СО₂ в ультраосновных расплавах высоких давлений. Возможно, СО₂ образует с кремнием соединения вида Si(CO₃)₄⁻⁴ и при его переносе играет роль транспорта. Предполагаемый эффект описывается уравнением

2Mg⁺² + Si(CO₃)₄⁻⁴ = Mg₂SiO₄ + 4СО₂.

Именно это обстоятельство объясняет присутствие углекислоты в высокотемпературных расплавных включениях. В этом случае роль СО₂ аналогична роли воды в расплавах.

По геологическим данным пироп Pyr находятся совместно с клино-(CPX, диопсид) и ортопироксенами (OPX, энстатит), редко — кианитом. По геохимическим данным пироп термодинамически равновесен пироксенам^[54][8]. Возможное уравнение образования граната имеет вид

… = … {Grn + Px} + ….

Связь граната с пироксеном генетическая: они образуются при разложения некоего протовещества. Оценка его состава опирается на данные количественных измерений соотношений Pyr/CPX в реальных объектах, чаще всего Pyr/CPX ≈ 1.0. Тогда

Mg₃Al₂Si₃O₁₂ + CaMgSi₂O₆= CaMg₄Al₂Si₅O₁₈

— это кордиерит, изначально образованный, возможно, по реакции

(An)+ (Mg-Ol) + 2(Mg-Opx)

рассмотренной J.M. Mclelland и P.R. Whitney^[55].

Рудная минерализация

В этом разделе пассмотрены рудные минералы, встречающиеся в интрузивных породах: магнетиты, ильмениты и шпинели.

Магнетит, ильменит.

I.Изучен изотопный состав кислорода в акцессорных магнетите (Mt) и ильмените (Il) кислых магматитов и метаморфитов. Равновесие минералов с Н₂О, СО₂ и СО не подтверждается, зато выявлено равновесие с рутилом TiO₂, соответствуя возможному образованию системы Mt(Il)- Ru при разложении ферропсевдобрукита или ильменита (П. Я. Ярош^[56]; П. Р. Бусек, К. Келль,1966; и т. д.) по реакциям

FeT₂iO₅ →Il + Ru;

6Il+ О₂ → 2Mt + 6Ru.

В отдельных случаях возможно образование минерала за счёт разложения ульвошпинели

6Fe₂TiO₄ → [6FeTiO₃ + 2Fe₃O₄].

Тогда по паре Mt-Il изотопная температура образования магнетита равна Т_изот = 853 °C. Для сравнения в нельсонитах Джугджурского и анортозитов Каларского массивов (А. Н. Соляник и др.,^[57]. ) Т_изот(Mt) = 820—1100 °C, а в рудных габбро- 560—900°С.

Возможен и другой механизм: образованием Mt — за счёт разложения ильменита по реакции

3FeTiO₃ + 2O⁻² → [Fe₃O₄ + 3TiO₂].

Тогда Mt находится в изотопном равновесии с рутилом (Ru). В этом случае Mt образуется при Т_изот ≈ 450 °C. Такие Т_изот(Mt) вполне возможны. Так на рудопроявлении р. Кюэричи жилообразные магнетит-гемоильменитовые руды образованы при Т=430-570°С^[57].

В метаморфических породах Il и Mt формируются в равновесии с Ru при Т_из =400 — 500 °C. Если же рассматривать Il как продукт разложения ульвошпинели, то в ассоциации с Mt их Т_из = 458 °C. Магнетит не может быть образован за счёт разложения Il, поскольку в противном случае температуры образования (Т_из = 1100 −2000 °C) геологически не реальны.

В железорудной формации Biwabik в Сев. Миннесоте^[58] скарнового типа изучена пара магнетит- кварц. Полученные данные дают температуру образования магнетита в 500—550 °C при условии равновесия его с СО₂. По материалам Синякова В. И. (1978), А. М. Дымкина и др.^[59] по результатам декрепитации Т_обр(Mt) в скарнах колеблется в пределах 420—530 °C. В скарнах^[58] при Т_изот = 525 °C магнетиты образуются при распаде сидерита по реакции

(а) 3FeCO₃ + O⁻² → [Fe₃O₄ + 3CO₂].

Здесь не ясна ситуация с углеродом: нет соединения, в равновесии с которым (по углероду) находился бы CO₂. В. Н. Загнитко и др.^[60], И. П. Луговая^[61], ссылаясь на эксперименты, приводят более приемлемые реакции

(б) 3FeCO₃ = [Fe₃O₄ + 2CO₂] + CO (безводные среды с удалением газа);

(в) 6FeCO₃= [2Fe₃O₄ + 5CO₂] + C (медленное удаление газа).

Таким образом, по изотопным данным намечаются схемы образования Mt:

• 1) разложение ильменита;
• 2) разложение ульвошпинели;
• 3) разложение сидерита (или железосодержащих карбонатов).

II. Геохимические (петрохимические) данные. Изучены преимущественно магнетиты Украинского щита. Mt преимущественно скарнов равновесен пироксенам относительно элементов (Fe, Mg). Следовательно, исходное уравнение должно иметь вид

.… = … + {Px + [Mt} + CO₂] + ….

Н. А. Елисеев^[62] при описании контактовых роговиков упоминает реакцию

CaMg(CO₃)₂ + 2SiO₂ = CaMg(SiO₃)₂ + 2CO₂.

Если вместо доломита взять анкерит Ca(Mg, Fe)(CO₃)₂, брейнерит (Mg, Fe)CO₃ или сидероплезит (Fe,Mg)CO₃, то при метаморфизме карбонатных пород можем получить искомую реакцию, например,

3Ca₂MgFe(CO₃)₄ + 6SiO₂= 3CaCO₃ +{3CaMg(SiO₃)₂ + [Fe₃O₄} + 8CO₂] + CO

(Ca₂MgFe(CO₃)₄ — анкерит), хотя здесь не ясна ситуация с углеродом кальцита.

Шпинелиды.

Хотя магнетит также относится к шпинелидам, однако эдесь будут рассмотрены другие шпинелиды. Изучены хромиты и хромшпинелиды из трубок Удачная и Мир (В. В. Ковальский^[54][9]). Т_обр ≈ 1240 °C (М. Н. Таран и др.,^[63]), Р не имеет четких значений. Распределение элементов для наиболее вероятной схемы образования минералов описывается уравнением

ln[(Cr+Fe⁺³)/Al] = −1.27 ln(Mg/Fe⁺²) + 2.65.

При интерпретации учитывались разные минералы, а также свободные элементы в твёрдом, расплавленном и газообразном состояниях и их одновалентные газообразные ионы на интервале Т=800-1600°С и Р= 1-45 кбар. В результате параметры уравнения соответствуют Т= 1200 °C и Р ≈ 25-30 кбар для реакций обмена:

FeCr₂O₄ + Mg = MgCr₂O₄ + Fe;

0,5MgAl₂O₄ + Cr = 0,5MgCr₂O₄ + Al.

В этих выборках не определена изотопная природа кислорода из-за отсутствия соответствующих данных. Присутствие элементов в свободном состоянии не совсем понятно. По данным О. Л. Кускова и Н. И. Хитарова^[64] это состояние элементов наблюдается при сверхвысоких давлениях (более 100 кбар). По экспериментам^[65] при ударном воздействии на Bi и флогопит образуются свободное Fe и шпинелиды.

Сульфидная минерализация

Изучались преимущественно сульфиды (чаще всего пириты, галениты, сфалериты, пирротин) и реже сульфаты (обычно бариты, гипсы и ангидриты). Распространенность их различная: сульфиды более характерны для высокотемпературных образований, хотя пириты (и марказиты) весьма часто образуются в экзогенных условиях. Как правило, высокотемпературные (гидротермы, магматические) сульфиды имеют облегченное (около 0 %о) значение δ³⁴S и узкие пределы её колебаний. В отдельных случаях величина δ³⁴S может опускаться до более низких значений (до -(20 — 25)%о). В целом, считается, что основная масса сульфидов образована при участии глубинной (метеоритной) серы. Сульфаты- преимущественно зоны окисления и низкотемпературные, хотя встречаются и высокотемпературные (ангидриты) образования, например, в норильских Cu — Ni рудах. Для сульфатов более характерны высокие значения величины δ³⁴S, достигающих +(25 ÷ 30) %о.

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

Выявлено статистическое преобладание распределений: с уменьшением Т и удалением от рудного тела в жильных и вкрапленных месторождениях происходит обеднение, а в массивных рудах — обогащение изотопом S³⁴. Эксперименты позволяет предположить, что в рудах первого типа сульфиды изотопно равновесны соединениям окисных, а второго — восстановленных форм серы. Наблюдается слабая зависимость состава форм S от морфологического типа руд. Сульфиды в большинстве случаев изотопно равновесны ионам S⁻², SO₄⁻², HSO₄⁻¹, H₂S. Формирование сульфидов происходит по гипотетической реакции вида

… → … + МеS (или МеS₂) + S⁻² + … .

В печати обычно обсуждаются реакции типа

Me(HS)₃⁻¹ → MeS + HS⁻¹ + H₂S [К.Грейнджер, 1976] или

Me^+x(HS)_n^x-n → MeS + HS⁻¹ + H₂S [И. Л. Ходаковский,^[66]].

Однако здесь присутствует ион HS⁻¹, не выявляемый изотопными данными. Комплексы Me(HS)₃⁻¹ имеют вид [Zn(HS)₃]⁻¹; [Pb(HS)₃]⁻¹; Zn(HS)₂ и пр. Эти комплексы больше соответствуют существованию вместо иона HS⁻¹ полисульфидов вида [MeS₂]⁻² , [MeS₄]⁻² и др., связываемых с миграцией металлов Hg, As, Sb, Mo и др. [В. В. Щербина,^[67] 1962] и растворимых в хлоридных и гидрокарбонатных растворах [Н. И. Говоров и др.^[68]]. Но для этих условий необходим высокий потенциал серы, что не всегда имеет место. Более приемлема гипотеза И. Г. Ганеева^[69], по которой сера выделяется из бессульфидных растворов по реакции

Pb(OH)₄ + 2H₂S → PbS + S⁻² + 4H₂O.

Но здесь не ясны дальнейшая судьба иона S⁻² и изотопная природа воды.

Надо заметить, учитывая опыт интерпретации данных по кальцитам, возникновение в правой части уравнений ионов HS⁻¹ и HSO₄⁻¹ сомнительно; кроме того, ион HSO₄⁻¹ должен вести к образованию гипса, что практически не наблюдается. Появления этих ионов обусловлено, видимо, недоучетом некоторых факторов, возможно, влияние SO₃ или диффузии H₂S.

На Кызылсайском редкометальном месторождении (Ю.Казахстан) изучен пирит кварц- кальцитовых и кварцевых пострудных жил (размер индивидов пирита достигает 5 — 6 мм). Для серы пиритов характерны низкие величины δ³⁴S (до (-30 %о)), уменьшение с δ³⁴S падением Т (по данным гомогенизации ГЖВ и изотопии кальцитов). Эти признаки говорят о возможной роли в выделении Py окисленных форм S, образованных из-за разложения гипотетических тиосульфатов или сульфитов по возможным схемам:

Fe(S₄O₆) → FeS₂ + 2SO₃, или

Fe(S₂O₃) + H₂O → FeS₂ + SO₄⁻² + 2H⁻¹ (Б.Грейнджер, К. Г. Уоррен, 1976).

Однако в этих случаях не объясним изотопный состав кислорода SO₃ и SO₄⁻² и отсутствие парагенных пиритам гипсов в кварц- кальцитовых прожилках. Возможность существования тиосульфатов и сульфитов в гидротермальных условиях показана в работах С. В. Кушнир (1989 г.) и Т. М. Сульжиевой^[70] и др.

Заключение

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

Примеры:

1.При анализе влияния СО₂ на плавление силикатов рассматривались^[71] реакции карбонатизации пироксенов (Ens) и оливинов (Ol) под давлением:

(1) Ol + Dio + СО₂ → Ens + Dol;

(2) Ol + Dol + СО₂ → Ens + Mgt (магнезит);

(3) Ens + Dol + СО₂ → Mgt + Qw.

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

2. В. И. Лучицкий^[72], описывая замещение роговой обманки (далее Amp), приводит реакцию

5Amp+ 7W → 2Ep+ Chl+ Act+ Qw+ 2H₄Mg₂Al₂SiO₉

(Act — актинолит, W — вода) и пишет, что

"Обыкновенно одновременно развивается эпидот Ep (более высокотемпературный (здесь и далее выделено нами)) и хлорит Chl (более низкотемпературный)… "" (вид заключений, весьма частый в петрологии).

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

Другие примеры приведены в таблице.

№ п.п.	Формула реакции	Источники
0001	Px + Spin → Grn + Ol	Leeman W.P., 1970
0002	Cor + 4Hyp → 2Grn + 3Qw	Hensen B.J., 1971
0003	8Stav + 12Qw → 4Grn + Chl + 30Kya	Федькин В. В., 1975
0004	2Bio + 7Cor → 10Grn + 6Mu + 3Qw	Thompson A.B., 1976
0005	4Bio + 3Cor → 6Grn + 4Kf + 4W	Currie K.L., 1971
0006	3Bio + 4Mu + 27Qw → 8Grn + 14Kf + 14W	Thompson A.B., 1976

В этой таблице практически все минералы в правых частях уравнений в реальных условиях образуются при разных термальных условиях, а ряд пар минералов, например, Grn + Chl, — и при разных барических условиях. В таблице: Px — пироксен, Sp — шпинель, Kya — кианит, Stav -ставролит.

Примечания

1. ↑ ^{1 2} Войткевич Г. В., Мирошников А. Е., Поваренных А. С., Прохоров В. Г. Краткий справочник по геохимии. М.: Недра, 1979
2. ↑ Ермаков Н. П. Исследование минералообразующих растворов. Харьков: изд-во Харьковского университета.1950
3. ↑ Ермаков Н. П. Геохимические системы включений в минералах. М.: Недра, 1972.
4. ↑ Карский Б. Е., Зорин Б. И., Фортунатов С. П. Исследование структурных разновидностей Мамских слюдоносных пегматитов методом декрепитации. //Изв. ВУЗ-ов. Геология и разведка, 1974. № 3.
5. ↑ Ермаков Н. П., Долгов Ю. А. Термобарогеохимия. — М.: Недра, 1979.
6. ↑ Ермаков Н. П., Критерии познания генезиса минералов и среда рудообразования, Львов, 1949 (Минералогический сборник, № 3. Приложение первое)
7. ↑ Мельников Ф. Н., Прокофьев В. Ю., Шатагин Н. Н. Термобарогеохимия. М.: изд-во Академический проект, 2008, 220 с ISBN 978-5-8291-0958-5/
8. ↑ Заварицкий А. Н. Изверженные горные породы. М.: изд-во АН СССР, 1961.
9. ↑ Термодинамический режим метаморфизма. Л.: Наука, 1976.
10. ↑ Bigeleisen Jacob. Quantum mechanical foundations of isotopechemistry. In book: Isotopes and chemical principles. Red. Peter A. Rock. Ser. 11. Washington, 1975.
11. ↑ Bottinga Y., Javoy M. Comments on oxygen geothermometry.//Earth and Planet. Sci. Let., 1973, 20, P. 250—265.
12. ↑ Гриненко В. А., Гриненко Л. Н. Геохимия изотопов серы. М.: Наука,1974. 271 с.
13. ↑ Носик Л. П. Изотопные методы при изучении минералообразования. М.: Наука, 1986. 222 с.
14. ↑ ^{1 2} Clayton R.N. and et al. Limits on the effect of pressure on isotopic fractionation.//Geochim. Cosmohim. Acta,1975,39,8,1197-1201
15. ↑ Matcha R.L. Theoretical analysis on the electronic structure and molecular propeties of the alkali halides.III Sodium chloride.//J/Chem/Phys.,1968, V.48, 1. P.335 — 343.
16. ↑ Krisher L.C., Norris W.G. Microwave spectrum of silver bromide.//J. Chem. Phys.,1966, 44, 1, 391—394.
17. ↑ Lauric V.W. Note on the determination of molecular structure from spectroscope data.//J.Chem. Phys., 1958, 28, 4, 704—716
18. ↑ Tellgrem R., Olovsson J. Hydrogen bond studies. XXXXVI, The crystal structures of normal and deuterated solid hydrogeh oxalate monohydrate NaHCO₂O₄ •H₂O and NaDCO2O4•D2O.// J. Chem. Phys., 1971,54, 1, 127—134
19. ↑ Harting P. Der thermodynamische kohlenstoffisotopiceffekt im system CH₄ — H₂O //Isotopenprexis. 1981. Vol. 17. № 4. P. 164—167.
20. ↑ ^{1 2 3 4 5 6 7} Макаров В. П. Изотопные геотермометры./XIII научный семинар «Система „Планета Земля“». — М.: МГУ, РОО «Гармония строения Земли и планет», 2005, стр. 228—253.
21. ↑ ^{1 2 3} Макаров В. П. Изотопные геотермометры./XIII научный семинар «Система „Планета Земля“». — М.: МГУ, РОО «Гармония строения Земли и планет», 2005, стр. 228—253.
22. ↑ Саксена С. Термодинамика твёрдых растворов породообразующих минералов. М.: Мир, 1975
23. ↑ Курепин В. А. Термодинамика распределения компонентов между породообразующими силикатами. //Геохимия, 1982. 1. С. 71 — 98
24. ↑ Перчук Л. Л. Биотит — гранатовый геотермометр.//Докл. АН СССР, 1967, Т.177, № 2. С.411-414.
25. ↑ Вуд Б., Фрейзер Д. Основы термодинамики для геологов. М.: Мир, 1981. 184 с.
26. ↑ Перчук П. П. Пироксеновый барометр и «пироксеновые геотермы»//Доклады АН СССР, 1977, 233, 6. С.1196-1199
27. ↑ Шемля М., Перье Ж. Разделение изотопов. — М.: Атомиздат, 1980. С.169.
28. ↑ ^{1 2} Макаров В. П. О механизме выделения минералов. /Материалы XVI научного семинара «Систама планета Земля» М.:РОО «Гармония строения Земли и планет», 2008, С.265 — 300. ISBN 978-5-397-00196-0
29. ↑ Кучер М. И., Макаров В. П. и др.//Геоэкология. 1996. № 2. С.21.
30. ↑ ^{1 2} Макаров В. П. Некоторые свойства геохимических геотермометров./XIV научный семинар «Система „Планета Земля“». М.: РОО «Гармония строения Земли и планет», 2007 °C.142 — 163.
31. ↑ Перчук Л. Л., Кротов А. В. Петрография слюдистых сланцев пояса Тана в южном тектоническом обрамлении Лапландского гранулитового комплекса.//Петрология, 1998, 6,2. С.165-196
32. ↑ ^{1 2 3 4} Ваганов В. И., Соколов С. В. Термобарометрия ультраосновных парагенезисов. М.: Недра, 1988. 149 с.
33. ↑ Венедиктов В. М. и др. Породообразующие пироксены Украинского щита. Киев: Наукова думка, 1979
34. ↑ Перчук Л. Л. Коррекция биотит-гранатового термометра для случая изоморфизма Mn= Mg + Fe.//Докл. АН СССР, 1981, Т.256.№ 2. С.441- 442.
35. ↑ Перчук Л. Л., Аранович Л. Я. Усовершенствование Bi- Gr- геотермометра: коррекция на содержание фтора в биотите.//Докл. АН СССР, 1984, Т.277, № 2. С.471- 475.
36. ↑ Макаров В. П. Некоторые свойства геохимических геотермометров./ XIV научный семинар «Система „Планета Земля“». М.: РОО «Гармония строения Земли и планет», 2007. C.142 — 163
37. ↑ Макаров В. П. Методические аспекты геохимической геобарометрии./ XIII Геологический Съезд республики Коми. «Геология и минеральные ресурсы Европейского Северо- Востока России: новые результаты и перспективы». — Сыктывкар, 1999, стр. 60- 63..
38. ↑ Устинов В. И., Гриненко В. А. Изотопные равновесия в минеральных ассоциациях.// Геохимия. 1990, 3, С. 307—315.
39. ↑ Макаров В. П. К определению «размерности» физических величин (метрологический аспект)./Мат-лы международной научно-практической конференции «Перспективные инновации в науке, образовании, производстве и транспорте». Одесса: издание «Черноморье», 2007, Т. 15. С.5 — 12.
40. ↑ Румер Ю. В. и др. Термодинамика. Статистическая физика и кинетика. М.: Наука, 1972
41. ↑ Сивухин Д. В. Общий курс физики. Термодинамика и молекулярная физика. М.: Наука, 1973
42. ↑ Макаров В. П. Некоторые свойства геохимических геотермометров./ XIV научный семинар «Система „Планета Земля“». М.: РОО «Гармония строения Земли и планет», 2007 °C.142 — 163.
43. ↑ ^{1 2} Макаров В. П. О механизме выделения минералов. /Материалы XVI научного семинара «Систама планета Земля» М.:РОО «Гармония строения Земли и планет», 2008, С.265 — 300. ISBN 978-5-397-00196-0
44. ↑ Anderson T.F. et al. A stable isotope study on the origin of permineralized peat zones in the Herrin Coal.//J. Geology, 1980, V.88,№ 6, P.713- 722.
45. ↑ Самойлов В. С. Карбонатиты (фации и условия образования). М.: Наука, 1977
46. ↑ Анфилогов В. Н., Бобылев И. Б. Термодинамический анализ растворимости воды в силикатных расплавах.//Геохимия, 1985, 9. С.1277 — 1289
47. ↑ Reesmon A.J., Keller W.D. Calculationof apparent stsndart5 free energies of formation of six rock — forming silicate minerals from solubility data///Amer.Miner., 1965, V.50, 10. P.1729 — 1739
48. ↑ Рыженко Б. Н., имельникова Г. Л.. Шваров Ю. В. Моделирование на ЭВМ формирования химического состава природных растворов при взаимодействии в системе вода -порода.//Геохимия, 1981, 4,. С. 481—495.
49. ↑ Алексеев В. А. и др. Кинетика альбитизации калиевого полевого шпата. Условия и механизм реакции.//Геохимия,1989,4.С.559-570
50. ↑ Ганеев И. Г. Строение и свойства гидротермальных растворов. Форма миграции минерального вещества.//Известия АН СССР, сер. Геологическая, 1977, 3, С.22-35.
51. ↑ Перчук Л. Л., др. — Биотит-гранат- кордиеритовое равновесие и эволюция метаморфизма. М.: Наука, 1983
52. ↑ Наумов В. Б. Определение концентрации и давления летучих компонентов в магматических расплавах по включениям в минералах.//Геохимия, 1979, 7. С.997- 1007.
53. ↑ Барсуков В. Л., Когарко Л. Н., Романчев Б. П. Физико-химические параметры дифференциациищёлочно-базальтовых магм.//Геохимия, 1981, 12. С.1816 — 1844
54. ↑ ^{1 2} Макаров В. П. О механизме образования минералов в кимберлитах.2. Природа гранатов и шпинелидов [1]
55. ↑ Mclelland J.M., Whitney P.R. A generalized garnet-forming reaction for metaigneous rocks in the Adirondacks. //Contr.Miner.Petrol.,1980, V.72, 2. P.111 — 122
56. ↑ Ярош П. Я.О выделении рутила при метаморфических изменений ильменита.// Записки Всесоюзного минералогического Общества, 1955,Т.84,№ 4.С.434 — 442
57. ↑ ^{1 2} Соляник А. Н., Коханова В. Г.. Сапин В. И., Соляник В. А.//Геология рудных месторождений, 1984, 4. С.27- 33.
58. ↑ ^{1 2} Perry E.C., Bonnichsen B. Quartz and magnetite: oxygen-18 — oxygen- 16 fractionation in metamorpposed Biwabik Iron formation.//.Science. 1966. V.153, 3735. P. 528—529.
59. ↑ Дымкина А. М.,Мазуров М. П., Николаев С. М. Петрология и особенности формирования магнетитовых месторождений Ирбинского рудного поля (В. Саян).//Труды Института геологии и геохимии СО АН СССР, вып. 268. Новосибирск: Наука, 1975.
60. ↑ Загнитко В. М., Луговая И. П. Изотопная геохимия карбонатных и железо-кремнистых пород Украинского щита. Киев: Наукова Думка, 1989. 315 с
61. ↑ Луговая И. П. Изотопный состав кислорода карбонатов железа и продуктов их термической диссоциации./Геол. ж., 1973, Т.33,3. C.42-48
62. ↑ Елисеев Н. А. Метаморфизм. М.: Недра, 1963.
63. ↑ Таран М. Н., Бондарь Ю. В.. Мацюк С. С.. Польшин Э.Вэ Спектроскопические исследования оливина и хромита из метеорита Марьялахти. //Минер. журн. 1989, Т.11, 5. С. 80- 87.
64. ↑ Кусков О. Л., Хитаров Н. И. Термодинамика и геохимия ядра и мантии Земли. М.: Наука, 1982. 278 с.
65. ↑ Мороз Т. Н. и др. О некоторых аспектах использования методов колебательной спектроскопии и порошковой дифрактометрии в геологических исследованиях./Сб. Геология и минеральные ресурсы европейского Северо-Востока России. Новые результаты и новые перспективы. Т.IV.Сыктывкар: КолНЦ Ин-та Геологии,1999.С.33-35.
66. ↑ Ходаковский И. Л. О гидросульфидной форме переноса тяжёлых металлов в гидротермальных растворах.//Геохимия,1966, 8. С.990 — 971
67. ↑ Щербина В. В. Способы выяснения форм переноса химических элементов в геохимических процессах.//Геохимия, 1962, 11. С.945 — 952
68. ↑ Говоров Н. И., Благодарова Н. С., Мукосеева З. Л. Роль карбонатов и галогенидов щелочей в гидротермальном переносе Pb и Zn./Геохимия и минераллгия магматических образований. Владивосток: изд-во ДВГИ, 1966. С. 40 — 47.
69. ↑ Ганеев И. Г. Строение и свойства гидротермальных растворов. Форма миграции минерального вещества.//Известия АН СССР, сер. Геологическая, 1977, № 3, С.22-35.
70. ↑ Сульжиева Т. М., Волков И. И. Тиосульфаты и сульфиты в термальных и гидротермальных водах.//Геохимия, 1982, 7. С.1043-1047
71. ↑ Кадик А. А., Френкель М. Я. Декомпрессия пород коры и верхней мантии как механизм образования магм. М.: Наука. 1982
72. ↑ Лучицкий В. И. Петрография. Т.2. Горные породы. Л.: Госгеолиздат, 1949. С. 366

Ссылки

• http://alkaline2002.narod.ru
• http://www.093/help-rus-student.ru/text/131/htm
• http://www.mineral.spmi.ru/01brief/zavaritsky.html
• http://www.rbtl.ru/zzzz/zavaric.htm
• http://alkaline2002.narod.ru
• http://alkaline2005.narod.ru
• http://www/msu.ru/deps/petro/perchuk.htm