Бесплатный автореферат и диссертация по геологии на тему
Термодинамика минералов и геологическая термобарометрия
ВАК РФ 04.00.02, Геохимия

Автореферат диссертации по теме "Термодинамика минералов и геологическая термобарометрия"

л п

•гч — л 1« М

У4 .3

Национальная Академия наук Украины Институт геохимии, минералогии и рудообразования

На права! рукописи УДК 636.7:549

КУРЕПИН Виктор Александрович

ТЕРМОДИНАМИКА МИНЕРАЛОВ И ГЕОЛОГИЧЕСКАЯ ТЕРМОБАРОМЕТРИЯ

Специальность 04. 00. 02 - геохимия

Автореферат диссертации на соискание ученой степени доктора геологических наук

Киев 1995

Диссертацией является рукопись. Работа выполнена в Инс

титуте геохимии, минералогии и рудообразования НАН Украины

Официальные оппоненты: доктор химических наук, член-корреспондент РАН В.С.Урусов

доктор геолого-минералогических наук, член-корреспондент НАН Украины Р.Я.Белевцев

доктор геолого-минералогических наук А.А.Вальтер

Ведущая организация : Институт экспериментальной

заседании специализированного совета Д. 50.08.01 при Институте геохимии, минералогии и рудообразования HAH Украины. Адрес: 252680 , Украина, Киев-142, пр. Палладина, 34 для телеграмм Киев-142, Минерал С диссертацией можно ознакомиться в библиотеке ИГМР HAH Украины

Автореферат разослан " " _ 1995 г.

Ученый секретарь специализированного совета

минералогии РАН (Черноголовка)

Защита состоится

доктор геолого-минералогических наук

С. Г. Кривдик

ч

ВВЕДЕНИЕ

Геологическая термобарометрия, начало которой положили Д.С.Коржинский, Г.Рамбырг и Т.Барт, благодаря работам Л. Л. Перчука, С.К.Саксены. Б.Дж. Вуда и многих других исследователей превратилось в интенсивно развивающееся научное направление, применение методов которого составляет непременную часть современных петрологических исследований. Для дальнейшего прогресса этого направления необходимо повышение надежности термодинамической основы геологической термобаромет-рии, определение и повышение точности сведений о РТ-зависи-мости констант равновесия для реакций между компонентами и соотношений активность-состав для минералов.

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

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

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

Таким образом, актуальность темы определяется 1) необходимостью расширения и укрепления термодинамической базы данных для геологической термобарометрии и других областей применения термодинамики минералов, 2) необходимостью теоретического обоснования оценки соотношений активность-состав в минералах с гетеровапентным изоморфизмом.

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

Задачи исследования: 1) изучение связи между пространственными и энергетическими параметрами минералов, между изменением этих параметров при изоморфных замещениях; 2) теоретическое Исследование термодинамики минералов с гетерова-лентным изоморфизмом; 3) разработка термодинамических моде-

лей минералов на кристалпохимической основ«; 4) термодинамический анализ распределения компонентов между минералами и способности минералов концентрировать определенные элементы; 5) исследование распределения компонентов между минералами и определение термодинамических условий образования высокоглиноземистых бердичевских гранитов как полигона для проверки термодинамических моделей минералов и калибровки геотермометров.

Новизна и научная значимость результатов.

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

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

Проведено теоретическое исследование термодинамики минералов с учетом ближнего порядка, обусловленного различием зарядов изоморфных атомов при гетеровалентном изоморфизме, и разработана термодинамическая модель локально электронейтрального (ЛЭН) твердого раствора, а также модель гетерова-лентного твердого раствора с равновесным нарушением локальной электронейтральности. Выведены уравнения вклада конфигурационной энтропии в активность компонентов минералов при отсутствии и наличии ближ-

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

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

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

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

3. Более электроположительные элементы концентрируются в минералах с меньшим энергетическим эффектом йС^д замещения атомов А атомами В и перераспределяются в минералы с большим энтальпийным аффектом ДНР* при повышении температуры. Показателями минералов - концентраторов определенных элементов могут служить объемные эффекты ВА-замещения атомов. Концентраторы более электроположителных элементов имеют

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

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

5. Бердичевские граниты возникли в интервале РТ-условий от ~500 МПа и 630"С до 670 МПа и 790"С при андапузит-силли-манитовом типе метаморфизма с температурным градиентом около 30"/км, который обусловил вертикальную зональность этих пород. Такое же соотношение между температурой и давлением метаморфизма характерно для кордиеритсодержащих гнейсов других районов Украинского щита.

Все научные выводы диссертации сделаны и сформулированы автором диссертационной работы.

Практическое значение работы определяется важной ролью достоверности термодинамической базы данных и оценок активность-состав для термодинамических расчетов в геологической

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

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

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

Автпм благодарен за помощь в исследовании минералов спи-

циальными методами сотрудникам ИГМР АН Украины С. В. Геворкьян (ИК-спектроскопия), A.M. Калиниченко (метод ПМР), Г. А.Малюк (маис-спектрометрия ), В.С.Мельникову (рентгеноструктурный анализ), Э.В.Польшину (метод ЯГР ), М. Н.Тарану (оптическая спектроскопия), Д.В.Уточкину (оптические исследования), сотрудникам лаборатории химического анализа и лаборатории рентгено-флгоресцентного анализа ИГМР, С. Б. Степченко за помощь в математических расчетах, сотрудникам ИЗМ РАН (Черноголовка) Л.Я.Арановичу и К.К.Подлесскому за помощь в проведении микрозондового анализа минералов. Автор признателен А. Я.'Парфеновой за помощь в оформлении диссертации.

Апробация работы. Основные результаты опубликованы в монографии, 40 статьях, 23 тезисах докладов на III, IV, V и VI Всесоюзных симпозиумах по изоморфизму (Киев,1974; Казань, 1977; Черноголовка, 1981; Звенигород,1988), X и XI Всесоюзных совещаниях по экспериментальной минералогии (Киев, 1978; Черноголовка, 1986). I и II Всесоюзных симпозиумах "Термодинамика в геологии" (Суздаль, 1985; Миасс,1988 ), VI и VII Всесоюзном петрографическом совещаниях (Ленинград, 1S81; Новосибирск, 1986), Симпозиуме Уппсальского университета "Систематика термодинамических данных" (Уппсала, Швеция, 1990), II Междунар. симпозиуме "Термодинамика природных процессов" (Новосибирск, 1992), XVI съезде ММА (Пиза,1994), Ежегодных семинарах экспериментаторов (ГЕОХИ, Москва), Выездных сессиях комиссии по изоморфизму БМО АН СССР.

Структура и объем работы. Диссертация состоит из введения, 6 глав и заключения, содержит 205 страниц машинописного текста, включает 67 таблиц, 47 рисунков и список литературы из 302 наименований.

Глава 1. СВЯЗЬ МЕЖДУ ПАРАМЕТРАМИ МИНЕРАЛОВ РАЗНЫХ ИЗОСТРУКТУРНЫХ РЯДОВ Многие пространственные и энергетические параметры минералов. такие как межатомное расстояние, мольный объем, теплоемкость, энтропия, удобно в первом приближении считать аддитивными величинами, равными сумме значений параметров компонентов. На принцип аддитивности опираются системы ионных и атомных радиусов, системы инкрементов объема (Blitz,1934). теплоемкости и энтропии (Robinson and Haas, 1987), энергии Гиббса и энтальпии образования (Chermak and Rimstidt,1989), которые были предложены для приближенного расчета этих параметров минералов.

Обнаружение закономерностей в отклонениях от аддитивности дает возможность значительно более точно предсказывать пространственные и энергетические параметры минералов. Неаддитивность проявляется при сравнении какого-либо параметра изоструктурных однотипных минералов элемента А, образующих ряд АХ (X - Xj, Хг, Х3,...), и аналогичных минералов ряда ВХ.

Исследованы соотношения межатомных расстояний R, мольного объема V, модуля объемного сжатия Кр, теплоемкости Ср , энтропии S и энтальпии образования ÄHf у наиболее изученных соединений - щелочных галогенидов. а также соотношения некоторых из этих параметров у оксидов и халькогенидов щелочных и щелочноземельных металлов и актиноидов, у соединений переходных металлов и других соединений. Все эти соотношения указывают на закономерность отклонений энергетических и пространственных параметров минералов от аддитивности.

Условием корреляции между свойствами F минералов рядов АХ

и ВХ является близкий характер химической связи соединений, объединяемых в один ряд. Отклонение от единицы углового коэффициента Ьр в линейных соотношениях между значениями какого-либо параметра Р минералов АХ и ВХ

Р(ВХ) - ар + ЬР Р(АХ) (1.1)

отражает степень отклонения от аддитивности. Обнаружено, что у таких соединений как щелочные галогениды или щелочноземельные оксиды и халькогениды, у которых и катионы, и анионы являются электронными аналогами, существует корреляция между коэффициентами % и Ьр уравнений для группы рядов (рис.1), которую можно выразить линейным уравнением

ар - а + Рр Ьр . (1.2)

Если АХ и ВХ - однотипные изоструктурные соединения с близким характером химической связи, то а - -р.

Существует общая закономерность для всех рассмотренных параметров , кроме V и для однотипных изоструктурных

соединений, у которых как катионы, так и анионы относятся к одной подгруппе или семейству элементов. Чем крупнее общий ион в ряду соединений ВХ, тем меньше различие параметров у членов одного ряда и поэтому меньше коэффициенты Ьр линейных уравнений (1.1). Различие Я, Кр. ср и Б в рядах щелочных га-логенидов с общим анионом уменьпается от фторидов МР к йоди-дам М1 и соответственно уменьшаются коэффициенты Ьр линейных соотношений Р(МР)-Р(МХ). То же характерно для рядов галоге-нидов с общим катионом и для рядов других рассмотренных соединений. Это отражает отрицательная величина коэффициентов Рр уравнения (1.2) для Я оксидов и халькогенидов лантанндов (Рв—607.6 пм ), для ср щелочных галогенидов (рс—53.87 Дж/К), Б щелочноземельных оксидов и халькогенидов (р3—90.0

II (ИХ).пи

/ СвС1

у ^

у '

' ' Л«| у

^ о х

• ' МлВг т' Ы1

/ ' А«вг

ЫС1 1.1Г

К«С1 0*>С1

г

Я (МЛ.ПЫ »«Г Л|гГ кг Ю* Саг

_I_и_:_I_

-V-

зоо

гзо

Б

№Х

. МС1

ч-*

КХ М1

ий

СгХ

Ся)

/

/

/ / КВг '

щ

/ №£1

250

Рис.1. Соотношение между межатомными расстояниями I? в щелочных фторидах МГ и других щелочных галогенидах МХ о общим анионом X - С1,Вг,I (А) и соотношение между коэффициентами а« и Ья линейных уравнений (1.1) для Н(МР)-ЩМХ) и 1?(1ЛХ)-Г?(МХ) (Б). Различаются соотношения между Я у галогенидов щелочных металлов и галогенидов серебра.

- И -

Дж/К) Знание коэффициента ftF для группы соединений позволяет рассчитать значения параметра F для всех членов любого ряда соединений, если известно значение F хотя бы одного соединения это ряда. Расчет этим способом коэффициентов линейных уравнений и параметров элементарной ячейки а - 2R и мольного объема V оксидов и халькогенидов лантаноидов МХ (X -0. S, Se, Те и М - U, Ра, Pa. Np, Ara, ТЮ типа NaCl позволил обнаружить инверсию в соотношениях а(МХ) - а(НО), при которой положительные значения коэффициента b сменяются отрицательными. Если у оксидов R (и V) возрастают от U0 к ThO в последовательности, соответствующей приведенному выше ряду актиноидов, у сульфидов и особенно у теллуридов различие R (и V) у членов одного ряда становятся значительно меньшими, то у теллуридов R (и V) уменьшаются от UTe к ThTe. Это показывает приближенность . системы ионных радиусов, неаддитивность межатомных расстояний и их зависимость от особенностей химической связи.

С увеличением размера общего иона в ряду однотипных изоструктурных соединений МХ различие V у членов ряда и, соответственно, коэффициент bv может как возрастать, как у щелочных галогенидов, так и уменьшаться, как у оксидов и халькогенидов лантаноидов. Соотношения между V некоторых рядов, например щелочноземельных оксидов МО и гидрогалогенидов МНХ, а также между ДНГ щелочных галогенидов сильно отклоняются от линейных и значительно точнее описываются квадратичным уравнением (1.3). К более точной оценке неизвестных значений параметров приводит сравнение рядов соединений с более близкой химической связью.

- 12 -

Глава 2. СВЯЗЬ МЕЖДУ ЭНЕРГЕТИЧЕСКИМИ И ПРОСТРАНСТВЕННЫМИ ПАРАМЕТРАМИ МИНЕРАЛОВ.

Проведено исследование соотношений между изменением энергетических и пространственных характеристик соединений АХ (X = Х[.Хо,Х3,..) при полной замене в них атомов А атомами В. Различие соответствующих характеристик структурных аналогов АХ и ВХ можно назвать эффектом полного изоморфного ВА-замещения вне зависимости от реальных пределов и возможности ВА-изоморфизма. Таковы знтальпийный АНрд - ДНГ(ВХ) -АНГ(АХ), »энтропийный ДБвл - 5(ВХ) - Б(АХ), объемный АУрд -У(ВХ) - У(АХ) и другие эффекты ВА-замещения.

Изменение энергетических и пространственных параметров при ВА-замещении, которое было бы постоянным при аддитивности, в действительности непостоянно. При замене ионов Р ионами Вг в щелочных фторидах МР (М - Ы, N8, К, №) при переходе от 1ЛХ к ЯЬХ знтальпийный эффект ДНРА постепенно уменьшается от 263.7 до 161.4 кДж/моль. При этом изменение межатомного расстояния АЯрд, которое при постоянстве ионных радиусов г должно быть постоянным и равным 2,(гВг-Гр), изменяется от 73.7 до 61.4 пм, а объемный эффект АУрд возрастает от 15.22 до 21.76 смя/'моль.

Связь между изменением энергии и межатомных расстояний. Анализ соотношений между энергетическими и пространственными параметрами в щелочных галогенидах, щелочных и щелочноземельных оксидах и сульфидах, а также в ряде других соединений показал, что в бинарных изоструктурных соединениях ряда АХ при замене атомов А атомами В изменения энергетических параметров Аи (энергии связи й, энергии атомизации Е, энтальпии образования АНГ) и изменения межатомных расстояний

чоо

m

гоо

шо

к •3

Ы -100 <t

-200

-300

-го

t

BrH

16

"BU& v»\ Y

Y

IjVV |

4g Na

r"

, V • ЪХат

AgLi^N

RrF

4

VIO

ч 50" î»

V

20

50

ta Ак,лм

Рис. 2. Связь между изменением энергии атомиаации ЛЕ и изменением межатомного расстояния ДО при полной изоморфном замещении 1-8 - соединения щелочных металлов MX структурного типа NaCl (1 - MF, 2 - MCI. 3 - МВг, 4 - MI, Б - LiX, б - NaX. 7 - КХ, 8 - RbX, 9 - CsX); 10-14 - соединения щелочноземельных металлов MX структурного типа NaCl (10 - MgX, 11 - СаХ, 12 -SrX, 13 - BaX)i 14-1? - соединения щелочных металлов МгХ структурного типа CaFg (14 - LI2X, 16 - NaeX. 16 - КгХ, 17 -Rb^X). Указаны типы замещения

- 14 -

Д{? связаны линейной зависимостью

Ди - ац + Ьи Д)? , (2.1)

которая для ДЕвл - ДКВА показана на рис. 1. Условием такой корреляции является близкий характер химической связи между атомами А и их окружением у всех соединений, объединяемых в ряд АХ. В частном случае ионы X являются электронными аналогами.

Связь между изменением энергии и мольного объема.

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

ДНВЛ - (ан + Ьн ( ДУВЛ)")П , (2.2)

с эмпирическими коэффициентами ал. Ьн. Коэффициенты шип постоянны для группы замещений и равны 1 или -1. Такое же соотношение существует между ДС^д и АУВЛ'.

Соотношения между энергетическими и объемными эффектами определенного ВА-замещения сильно зависят от особенностей химической связи между изоморфными атомами и их окружением. При одной и той же координации атомов А(и В) соотношения ДНВА - ДУвд, а также (^д- ДУрд различны у соединений элементов главных подгрупп Периодической системы и у соединений, содержащих атомы переходных металлов или элементов подгрупп 16 и 26 (Си, 7п и др.). Если не участвующие в замещении атомы относятся к элементам главных подгрупп и при ВА-замещении не происходит существенного изменения химической связи, например при взаимных замещениях ионов главных подгрупп или пе-

реходных металлов , то чем больше ДУСЛ, тем меньше АНВд и ДС^д Напротив, при замещении ионов элементов главных подгрупп ионами переходных металлов чем больше ДУВА , тем больше ДНВА и Дйвд.

Связь между изменением энтропии, теплоемкости.

межатомного расстояния и мольного объема £

При замене атомов А атомами В бинарных изоструктурных соединениях ряда АХ изменения межатомного расстояния АйВд и изменения теплоемкости АСрВд и энтропии ЛБрд связаны линейными соотношениями типа (2.1). Связь между изменениями мольного объема ДУвд, теплоемкости ДСрВА и энтропии Д^д также можно приближенно выразить этими же соотношениями.

Если при ВА-замещении не происходит значительного изменения характера химической связи, например, при замещениях электронных аналогов, то чем больше Д\ГВА. тем меньше ДСрВА и ДБвд. Напротив, при значительном изменении ионности связи, например при замещении элементов главных подгрупп элементами побочных подгрупп, чем больше Д\ГВд, тем больше ДСрдд и ДБрд.

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

Расчет энтальпийного эффекта ДНвд по объемному эффекту ДУВд ВА-замещения с помощью уравнения (2.2) дает возможность оценить энтальпию образования ДНГ соединения ВХ. используя сведения о ДНГ (АХ),

ДНГ (ВХ) - ДНГ (АХ) + ДНВА . (2.3) В большинстве случаев такие оценки мало отличаются от экспериментальных определений (табл.1). Случаи значительного расхождения этих величин(например у №2) могут быть связаны с

погрешностью эксперимента или с особенностями минералов и • заслуживают специального изучения. Установленные соотношения относятся только к решеточной энергии кристаллов, и поэтому из рассмотрения должна быть исключена энергия магнитного ра-

Таблица 1

Сравнение эмпирических значений дНг(ГеХ) соединений железа с результатами расчета по уравнению (2.2)

1 1 лн^ц, -ДНГ, кДж/моль

Общая I Простр. см3 1 кДж

формула 1 группа 1 мех 1 ГеХ

* 1 1 1 эмп. расч.

М1г 1 Р3ш1 -4.04 248 364 1 116 115.8

МВгг | Р3т1 -2.38 272,9 524.3 1 251.4 252.0

МС)? | 1?ЗЛ1 -0,62 299,6 641.3 1 341.7 338.0

04 Г РсЗЗд . -0,05 312,6 1428.4 1 1115.8 1113.1

МТШд ! яз 0,83 339,5 1572.8 1 1233.3 1237.6

М(0Н)2 1 Р3и1 0,87 356,5 925.5 1 56 9 589.6

МгБШ^ 1 РЬпт 1,32 346,6 2171.9 1 1478,8 1476. 9

МС0? I ЯЗс 1,35 351,3 1112,5 1 761,2 764.5

ДО^О» 1 ; РЬса 1.55 350,6 3089, 4 1 2388,2 2383.0

СаМ81г06 1 С2/с 1.76 356,7 3260,2 1 2843,5 2900.9

КМ8эА181З01 0(ОН)2 1.73 358.1 6215.3 1 5140.9

МоА1451$01в Ссст 1.92 363.4 9161.5 1 8434.7

Щ 1 Р4/юпш 2.34 375.7 1124.2 1 661 748.5

Параметры уравнения (2.2): а„ - 3,163*10"3, Ьн - -2.146*10"4.

в - 1, п « -1.

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

Глава 3. ТЕРМОДИНАМИЧЕСКИЕ МОДЕЛИ ГЕТЕРОВАЛЕНТНЫХ ТВЕРДЫХ РАСТВОРОВ Разработана модель гетеровалентного твердого раствора с локально электронейтральным размещением атомов, в котором отсутствуют заряженные дефекты, а также модель гетеровалент-ного твердого раствора с равновесным количеством таких дефектов.

Модель локально электронейтрального (ЛЭН) твердого раствора (АтВ^у)" (ХуY<-т)ь с гетеровалентным изоморфизмом AX-BY

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

Sconf Spand + SlEN . (3.1)

из которых первый (Srand) отражает хаотичное размещение атомов в позициях а, а второй CS^ е N ^ отражает возможность различной ориентировки ассоциатов АХ и BY при каждом из вариантов размещении атомов А и В в позициях а.

Комбинаторная задача оценки S^en решена на примере клино-пироксенового твердого раствора (Са,^..!)"2 (Mg,Alj.x)U1 SljOg ряда диопсид-жадеит, для которого найдено

- 18 -

Slen « x(l-x)Rln3 (3.2)

Поскольку в клинопироксене каждая позиция Mi окружена тремя позициями М2, т.е. координационное число К(М1-М2) - 3 ,то для Sconf ЛЭН твердого раствора (АхBt ¡t)a(XI[Y1 t)b с любым k принято уравнение

Sconf " -Rtxlnx + (l-x)ln(l-x)] + x(l-x)Rlnk . (3.3) которое согласуется с Sconf в твердых растворах при ничтожном содержании одного из компонентов. Scctlf тем больше, чем больше к . причем при к - 1 она такая же, как у оппозиционного раствора (АцВ^у ) , а при большом к близка к Sconf регулярного твердого раствора (АхВ1.х)а(ХхУ1.ж)ь.

Заселенность ближайшей координационной сферы. В ЛЭН твердом растворе (AzBj .х)а (X,Y1.„)b вероятность того, что ближайшими соседями атома А или В являются k - j атомов X и j атомов У ( 1 < j < к ). равна соответственно

CkV x^'Ml-x)1 и (W х^^О-х)1 . где Ск! - число сочетаний из к по j , равное kl/(J!(k-j)!). Нарушающие баланс зарядов группировки A(Yk) и В(Хк) становятся невозможными. Исследование методом ЯГР природных омфа-цитов (Aldridge et al., 1978) показало отсутствие группировок Fe?+(NaNaNa). нарушающих локальный баланс зарядов, что согласуется с моделью ЛЭН твердого раствора .

Энтальпия смешения. Анализ энергии регулярного и ЛЭН твердых растворов (АцВ!-,)®как суммы энергии взаимодействия между атомами, занимающими ближайшие позиции а, ближайшие позиции b и ближайшие позиции а и b , привел к следующим выражениям теплоты смешения:

АНMr,e - х(1 - xHU, + Uj> - U3) (3.4)

- х(1 - xUU, + U2 - (1 - 1/k) U3) (3.5)

где иг, Иг и из - энергия квазихимических реакций ДА + Во - 2АВ (1), XX + УУ - 2ХУ (2) и АУ + ВХ - АХ + ВУ (3). Поскольку из < 0, то чем больше к , тем меньше различие теплоты смешения регулярного и ЛЭН твердых растворов при прочих равных условиях и тем меньше выигрыш энергии в результате размещения атомов о соблюдением локального баланса зарядов.

Анализ 5СОпГ и ДНМ показывает, что чем меньше координационное число 1с , тем больше различие термодинамических свойств регулярного и ЛЭН твердых растворов {АхВ1_х)а(ХхУ1-х)ь.

Твердые растворы о равновесным нарушением локальной адектрояейтрадьностя. Повышение температуры делает термодинамически устойчивым состояние минерала о нарушенным локальным балансом зарядов атомов и более высокой энтальпией и энтропией, чем при соблюдении локальной электронейтрадьнос-ти. Анализ термодинамики твердого раствора

(АхВ1-х)а(Хх-аУ1-х-вХ*вУ*з)Ь. У которого группировки АУ* и НХ* имеют нескомпенсированные электростатические заряды, привел к следующим выводам.

ЛЭН состояние твердого раствора является термодинамически устойчивым при температурах Т < Н*/21? , где Н* - энтальпия квааихимической реакции образования заряженных дефектов АХ + ВУ - АУ* + ВX* .

Цри Т > Н*/2Я количество заряженных дефектов возрастает о ростом температуры и твердый раствор находится в состоянии, промежуточном между ЛЭН и неупорядоченным состояниями. В отличие от разбавленных твердых растворов при высокой концентрации обоих компонентов влияние ближнего порядка на термодинамические овойотва проявляется даже при высоких температурах (Т > н'/клг).

Многокомпонентные многопоаиционные ЛЭН твердые растворы. • Если в эквивалентных позициях твердого раствора нет изовалентных ионов разного сорта, то его Sconf можно рассматривать как сумму энтропийных вкладов Si и Sj, соответственно определяемых 1) числом перестановок ЭН-группировок, соответствующих миналам или ЭН-ассоциатам и 2) числом вариантов взаимной ориентировки ассоциатов при фиксированном положении одного из атомов каждого ассоциата. Первый вклад Sj в простейшем случае можно принять равным энтропии идеального молекулярного раствора

5} - Е N,lnN, , (3.6)

где Nj - мольная доля минерала 1. Второй вклад Sg определяется особенностями кристаллической структуры; в координационных соединениях он может быть значителен, тогда как в низкосимметричных некоординационных соединениях мал и им можно пренебречь. Поэтому в зависимости от структурных особенностей ЛЭН твердого раствора (А,В....(X,Y....)„ величина его Sconf может быть близкой к ее оценке как по модели регулярного , так и модели молекулярного раствора.

Если изовалентные атомы разного сорта занимают эквивалентные позиции, например атомы А", А",... в позициях а или атомы X', X".... в позициях Ь, то возможность их перестановок друг на место друга увеличивает конфигурационную энтропию на

5з - А'1п(А'/(А'+А"+...)) + А"1п(А"/(А'+А"+...)) (3.7) для атомов каждой валентности в каждом типе эквивалентных позиций. При сильной корреляции между размещением атомов разных позиций, например А' и X', вклад S3 определяется по вариантам размещения только атомов А' при зависимом размеще-

нии атомов X' или наоборот.

Вклад в активность компонентов.

Регулярный многопозиционный твердый раствор. В выражениях вклада БсопГ в активность компонентов (Темкин,1946; Томпсон, 1967) не учтен часто встречающийся случай, когда эквивалентные позиции в миналах заняты атомами разного сорта, что обычно при гетеровалентном изоморфизме. При этом как многокомпонентный многопозиционный твердый раствор (АаВь...)т(ХхУу...)п... , так и его миналы (Аа1ВЬ1. ..)„ (Хх1Уу1 ...)„..., (Аа2Вь;>. ..)„ (Хх2Уу2. ..)„... обладают конфигурационной энтропией. С учетом этого выведено общее уравнение вклада 5С0Пг в активность компонентов регулярного многокомпонентного многопозиционного твердого раствора

а, - П П (г/г!)21 (3.8)

X г

где г и 21 - числа атомов 2 , занимающих позиции X , в формульной единице твердого раствора и компонента 1 соответственно.

ЛЭН твердый раствор типа (АуВ< )" (ХуУ, _»)*' с гетерова-лентным изоморфизмом АХ ° ВУ . При хаотичном размещении ас-социатов АХ и ВУ в ЛЭН твердом растворе вклад Б! в активность компонента АХ определяет уравнение

аАХ - х к"(1"х)г (3.9)

Его величина близка к оценке по модели регулярного раствора, г .

равной х , если к велико, и к оценке по модели молекулярного раствора при малом к. Последнее характерно для многих низкосимметричных минералов, таких как биотит. Многокомпонентные многопозиционные ЛЗН твердые растворы

(А*. А", В',В",... )я (Х',Х",У,У",.,. )„ с гетеровалентным • изоморфизмом АХ » ВУ. Если пренебречь Бг, вклад 5сопг в активность компонента А'„Х'П выражается уравнением а(А'рХ'п) - НЛг>Хп(А'/(А'+А,,4...))"(Х'/(Х,+Х"+...))п (3.10) в которое входит суммарная концентрация НА|,хп в твердом растворе миналов А'ПХ'П и однотипных с ними миналов, а также отношения формульного числа ионов каждого сорта (А',Х') к сумме чисел изовалентных ионов, зачинающих те же структурные позиции.

Глава 4. ТЕРМОДИНАМИЧЕСКИЕ МОДЕЛИ НЕКОТОРЫХ МИНЕРАЛОВ Возможность оценить конфигурационную энтропию минерала и ее вклад в активность компонентов, опираясь на сведения или допущения о закономерностях размещения атомов в кристаллической структуре, позволяет создать простейшую модель твердого раствора, которая может быть уточнена, если имеющиеся данные позволяют учесть особенности других термодинамических функций и прежде всего нелинейную зависимость энтальпии от состава или катионного распределения. Это показано для ряда минералов.

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

Одновременный расчет по равновесным данным констант равновесия реакций с участием РеО и функций смешения несте-хиометрического вюстита не позволяет определить реальные термодинамические характеристики при использовании модели молекулярного раствора. Он приводит к отрицательной величине параметра смешения V при рассмотрении вюстита как твердого

раствора FeO-FeOj (Goel et al.. 1980; Fei and Saxena, 1986). Подобный расчет при ином выборе мольных количеств указывает на отрицательную (для Fe0-Fe02/3) или знакопеременную (для Fe0-Fe03/4) избыточную энергию Гиббса и теплоту смешения.

Только предварительный расчет функций AGf. AHf и ÄSf нестехиометрического вюстита позволяет определить реальные значения его функций смешения. Бюститовый твердый раствор Fe0-Fe3/i0 имеет положительную теплоту смешения (V/!2 »3.0 и W21 - 3.9 кДж/моль при 1000°С), а его энтропия смешения почти в два раза меньше, чем при хаотичном размещения катионов в эквивалентных позициях.

Простые шпинели. Анализ термодинамики простых шпинелей (At - J Bjj ) [ A, _х ] 04 на основе регулярной модели с учетом энергии взаимодействия между ближайшими катионами приводит к уравнению катионного распределения

-RTlnx2/(i-x)(2-х) - ÄH* + (l-2x)W . (4.1) где АН* - разница энтальпии обращенной В[АВ]04 и нормальной А СБа 104 шпинелей. Параметр W отражает нелинейную зависимость энтальпии от степени обращенности, причем

* " *!дв) + 1'2WUB] " W(abi . (4.2)

где W(Ав), ab]. W(abi - анергия квазихимических реакций 1/2АА + 1/2ВВ - AB. 1/2[АА) + 1/2[ВВ] - [AB] и АСА] + В[В] -В[А] + А[В] соответственно. Это уравнение катионного распределения (Курепин, 1975,1988), которое согласуется с моделью 0'Нила и Навроцки (O'Neill and Navrotsky, 1983), описывает температурную зависимость катионного распределения, его стремление к х - 2/3 при повышении температуры и предсказывает фазовые переходы I рода обращенных шпинелей в нормальные при пониженных температурах. Регулярная модель пригодна

для описания термодинамики гетеровалентного изоморфизма в атом минерале с координационной структурой.

Геленит. Термодинамическое описание внутрикристалличес-кого распределения катионов в гелените и в простых шпинелях одинаково, если принять структурную формулу геленита Са2 (А1х811-хГ(51хА1г.х[Ь8. По сведениям о катионном распределении в гелените (Сигаловская, 1988) находим АН* - -24.1 кДж/моль и \</=64.1 кДж/моль, что отражает значительную неидеальность смешения катионов А1 и и нелинейную зависимость энтальпии геленита от х.

Кордиерит. Присутствие молекул Н20 , С02 и других летучих компонентов в расширенных участках каналов, которые имеются в структуре кордиерита, увеличивает Бсоп! в связи с возможностью различного заселения этих полостей и превращает кордиерит в твердый раствор. При нахождении в такой полости одной молекулы Но0 или СО* структурную формулу Ми-кордиерита можно представить в виде

Мё2А14516018 [<Н20)п(С0г)„!!!-„-„ ] . где П - вакантная полость. Конечными членами твердого раствора являются кордиериты Сог(Н20), Сог(С02) и СогШ) , у которых все полости заняты молекулами Н20, С02 или вакантны. Вклад БсопГ в активность этих компонентов при хаотичном распределении молекул Н20 и С02 по занимаемым ими позициям равен п, ш и 1-п-ш соответственно. Из экспериментальных данных по растворимости Н20 (М1пка1й е1 а!.. 1979) и С02 (Ара-нович и др., 1981) в кордиерите следуют оценки констант равновесия реакций гидратации, и карбонатизации Согм? + НоО - СогМс(Н20) , 1ПКН20 - 5763/Т - 13.651, СогМг + СОг -СогМг<СОг> , 1пК,;0г - 1930/Т-12.73.

Близость рассчитанной энтальпии АН0 реакции гидратации (47.9 кДж/моль при 1000°С) к результатам калориметрического изучения (41.8*1.6 кДж/моль в стандартных условиях), которые указывают на линейную зависимость энтальпии от n (Carey, Navrotsky, 1992), подтверждает достоверность модели.

Рассчитана растворимость Н20 и С02 в кордиерите как функция температуры, давления и состава флюида, а РТ-зависи-мость доли вакантных полостей (1-п-т) в кордиерите выражена уравнениями

RT ln(l-n-m) - А + ВТ + СТ2 +. (D + ЕТ +FT2)Р (4.3) с эмпирическими коэффициентами А,В,С,..., рассчитанными для интервалов давления и разного состава водно-углекислотного флюида.

Эта модель кордиерита позволила объяснить и количественно оценить влияние парциального давления летучих компонентов на кордиеритсодержащие равновесия. Поскольку вхождение летучих компонентов в кордиерит понижает активность его компонента Сог(П) , то следствием этого является :вменение РТ-условий кордиеритсодержащих равновесий и изменение состава минералов, равновесных с кордиеритом при постоянных РТ-условиях.

Биотит. Сделан вывод, что из возможных схем изоморфного вхождения титана в высокотемпературном биотите реализуется схема М2+(ОН)2 « ТЮ2, поскольку 1) в биотите высокотемпературных пород (мантийные вкрапленники в базальтах, мигматиты и граниты гранулитовой фации) отношение (TI + M2+)/Si близко к 1, что характерно только для этой схемы, и 2) существует отрицательная корреляция между содержанием TiOj и Н20 в биотите (Ушакова, 1971). При такой схеме изоморфизма соотноше-

ние между суммой октаэдрических и тетраэдрических катионов V и числом атомов Т1* при расчете формулы на заряд 22 должно быть близко к Е* - 7 - 0,63611*, что обнаруживается в биотитах из бердичевских гранитов.

На допущение, что основными компонентами биотита являются флогопит, истонит, мусковит и безводный компонент ЩМ^Т!) (А1Б1з)01г, который можно назвать титаноксифлогопи-том (ПрЫ), и различия соотношений между атомами в этих ми-нап _ опирается предложенный метод расчета кристаллохими-ческой формулы биотита и суммарной концентрации его однотипных компонентов. Формульные коэффициенты кристаллохимической формулы биотита, полученные любым методом 1, должны быть умножены на коэффициент

а - 21/(( А1[ + 2( Еу,! + Б!, * 4 )) . где Ту >, - сумма октаэдрических катионов. Суммарная концентрация однотипных миналов определяется по числам атомов определенного сорта в кристаллохимической формуле и равна (3 -А1 - 2ТП/2, 3 - Б!. А1 + - 7 и Т1 соответственно для однотипных миналов на основе флогопита, истонита, мусковита и титаноксифлогопита. Этот метод не требует точного определения содержания Н20. как кислородный, не занижает суммарный заряд катионов, как расчет на заряд +22, не игнорирует октаэдрических вакансий, как катионный метод.

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

На разные допущения о закономерностях внутрикристалличес-кого размещения атомов опираются следующие рассмотренные модели биотитового твердого раствора:

AI. Многопозиционный регулярный твердый раствор, в котором атомы, занимающие эквивалентные позиции размещены беспорядочно. Вклад Sconf в активность компонентов определяется уравнением (3.8).

А2. Такой же твердый раствор, но без энергетически невыгодного соседства тетраэдрических катионов AI. что может быть реализовано при образовании доменов с подрешеткой (AI. Si)г (Sij>) Также используется уравнение (3.8).

A3. ЛЭН твердый раствор. Вклад Sconf в активность компонентов определяется уравнением (3.10).

Аналогичные модели Bl, В2 и ВЗ учитывают распределение катионов между транс- и циспозициями, а аналогичные модели CI, С2 и СЗ учитывают присутствие фтора в ОН-позициях и особенности взаимного твердого раствора.

Оценка активности компонентов в биотитах бердичевских гранитов показала, что модели, предполагающие разные закономерности размещения атомов и значительно различающиеся оценки Sconf, могут давать близкие оценки активности компонентов, что связано с малыми различиями энтропии смешения ÄSM. Примером может служить применение моделей AI. А2 и A3 для оценки активности аннитового и флогопитозого компонентов. Примером другого рода, когда оценки активности компонентов сильно зависят от принятой модели, является оценка активности оксианнитового и титаноксианнйтового компонентов с помощью моделей AI и A3.

Применение различных моделей биотита для термодинамического анализа экспериментальных данных Уонса и Ойгстера (Wones and Eugster.1965) по равновесию Ei + San + Mt показало, что лишь модель A3 пригодна для выражения температурной

зависимости кажущейся константы равновесия реакции Oan » San + Mt по данным для разного давления и разных кислородных бу-

константь

феров. Прочие модели пригодны для выражения этой •'лишь при определенном давлении и кислородном буфере.

Глава 5. ТЕРМОДИНАМИКА РАСПРЕДЕЛЕНИЙ КОМПОНЕНТОВ МЕВДУ МИНЕРАЛАМИ Соотношение между химическими потенциалами и константа равновесия. При равновесии минералов, как и при всяком ином гетерогенном равновесии, компоненты равновесных фаз связаны соотношениями (реакциями) типа

IL + mM + nN + ... = 0 , (5.1)

и подобные соотношения существуют между химическими потенциалами этих компонентов

1HL+ тци ♦ + ...- О . (5.2) Подстановка выражений химических потенциалов vi) - G°, +RTlnat приводит к уравнению константы равновесия

К(Р.Т) - exp(-ÜG°/RT) - (at,)"1 (а«)""(а«)"".... (5.3) где Дй° - энергия Гиббса реакции, a aj - активность компонента i.

Если РТ-зависимость констант равновесия известна, то, оценив их величину nb составу сосуществующих минералов, можно определить термодинамические условия минерального равно-

1

весия. Вначале был применен анализ обменных реакций AY + ВХ • АХ + BY между минералами переменного состава (А.В)Х и (A,B)Y для оценки температуры (Ramberg and De Vore, 1951; Kretz, 1961; Перчук. 1970 и др.), а затем была показана (Wood. 1975; Курепин, 1976) возможность оценки давлений на основе анализа реакций переноса, т.е. реакций между компонентами с изменением количественных соотношений между равно-

весными фазами. Реакции с участием летучих компонентов несут информацию о флюидном режиме.

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

К'(Р.Т.Х) = (Х^'ЧХм)""^)"0.... (5.4)

которая зависит не только от РТ-условий, но и от состава минералов.

Геотермометры и геобарометры и их чувствительность. Из РТ-зависимости констант равновесия от температуры и давления (й1пК/ <ТГ)р - АН°/Шг и (сПпК/ &?)■[ - /ИТ следует, что реакции с большим энтальпийным эффектом ДН° наиболее пригодны для определения температуры, а с большим объемным эффектом Д\/° для определения давления. В работе приведены термодинамические функции реакций между компонентам;! минералов метапелитовой системы,

Чувствительность геотермометров и геобарометров можно определить как способность к изменению состава минералов парагенезиса при изменении температуры или давления. Анализ константы равновесия обменной реакции между минералами (А1 _хВх)X и (А^к-ВяОУ показал, что при х • 1 - х' изменение состава при изменении температуры эх/эТ максимально при (1-2х)[1пх/(1-х)-ЛБ0/2И]-х(1-х)(Ча+Иь)/СТ+1 - 0 , (5.5) где и Уь - параметры смешения твердых растворов. При близости твердых растворов к идеальности наиболее чувствительным является геотермометр с энтальпией обменной реакции ДН°

25-5 кДж/моль (при ДБ"/ИМ). Чем больше неидеальность твердых растворов, тем меньшую энтальпию обменной реакции имеет наиболее чувствительный геотермометр вплоть до АН0 ~ 0 при О^И^/Ш" - 4.

В метапелитовой системе наиболее чувствительными являются

Ре и

вгСог геотермометр, у которого при близости смешения к идеальности ДН"~25 кДж/моль), и геобарометр СогСгБПЦ . который отличается наибольшим объемным эффектом минальных реакций (А1/4 >50 смя/моль). Напротив, такой геотермометр как 01-0рх или такой геобарометр как ОЮрхЦ (ДУ°~3 сма/моль) являются малочувствительными и оттого ненадежными показателями термодинамических условий.

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

Термодинамика Ре-кордиерита и геобарометр СогСгБПО.

Между применяемыми вариантами калибровки геобарометра СогСгБПО существует значительное противоречие, связанное с разной оценкой наклона йР/сГГ равновесия Сог+А1и+811+0, что ведет к различной оценке РТ-соотношений при петрогенезе. Проведен расчет термодинамических функций Ре-кордиерита по сведениям о термодинамике Мё-кордиерита, энтальпийном АНГеМг

или энтропийном ASpeMg эффекте FeMg-замнщения и сведениям об условиях равновесия Cor+Alm+Sil+Q, общим для всех экспериментальных работ (РН20~ЗШ МПа при 750°С). Последующий расчет термодинамики реакции

2/ЗА1ш + 4/3SÍ1 + 5/6Q - CorFe . (5.6) указал на положительный наклон tíP/dT в согласии с экспериментальными данными Л. Я. Арановича и К. К. Подлесского (1S82). Результаты расчета для реакции (5.6)

АС'.Дж/моль - 14497 - 25.62Т + 53.25(Р-0.1) (5.7) близки к ее термодинамической характеристике

AG0,Дж/моль - 16231 - 25.03Т + 53.25(Р-0.1) . (5.8) полученной нами по последним экспериментальным данным для равновесия Cor+Alffl+Sil+Q (Mukhopadyay and Holdaway, 1994).

Полученное уравнение геобарометра

Р. Ша - ( А - АН0 4 Т(В+ AS°) + 2R ln(NfeGr/ HF,Cor)) +

СТ2) / (AV° - D + ЕТ + FT2) (5.9)

учитывает зависимость давления минерального равновесия от состава водно-углекиилотного флюида, включает термодинамические функции AH0, AS0 и AV° реакции (5.6) и коэффициенты А, В, С,.. уравнения (4.3), определяющего зависимость доли вакантных полостей í-n-m в структуре кордиерита от РТ-усло-вий при заданном составе водно-углекислотного флюида. Для интервала 800-1100 К и 400-800 МПа при а„го-0.2 А—14237, В-21.338. С—9.67 10'3. D-64.141. Е—7.885 10"г, F-30.795 10"6.

Закономерности распределения компонентов между минералами. Закономерности распределения компонентов между минералами, которые рассматривались В.С.Соболевым (1948), Г.Рам-бергом (Ramberg, 1952), Л.Л.Перчуком (1970). В. С. Урусовым

(1975) и другими исследователей, рассмотрены с новых позиций в учетом связи между эффектами,изоморфного замещения. Энергия Гиббса, энтальпия и объемный эффект обменной реакции

АУ + ВХ = АХ + ВУ . (5.10)

представляют собой разность соответствующих эффектов полного изоморфного ВА-замещения в минералах МХ и МУ (М - А,В), в частности

Лй0' - 0И) - ДСВА (МХ). (5.11)

Поскольку обменная реакция смещается в правую сторону (ДС<0) при Лйрд(МУ) < ДСВЛ(МХ), то для определения минералов-концентраторов элемента достаточно сравнить величины Д^д минералов. Корреляция между энергетическими и объемными эффектами изоморфного замещения позволяет использовать для этой цели величины ДУвд- С учетом этого приходим к следующим правилам концентрации компонентов в минералах с одинаковой координацией изоморфных атомов.

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

Установленный Г.Рамбергом ряд распределения Ре и между оксидами и силикатами, отражающий увеличение отношения

Mg/Fe у сосуществующих минералов в последовательности оксиды - ортосиликаты - пироксены - амфиболы - слоистые силикаты и увеличение ионности связи между изоморфными атомами и их окружением, является частью более широкого ряда (табл. 1, рис.3), вдоль которого постепенно изменяют величину эффекты Ре1^-замещения Анвл, ДС^д и ДУВЛ и увеличивается способность концентрировать Более сильными концентраторами железа, чем оксиды и силикаты, являются соединения с меньшими эффектами ДСГвИ1 и ДУрвМе, такие как хлорид МС^, бромид МВг2. иодид М12, а наиболее сильным концентратором магния является фторид М?2 с наибольшими эффектами АСрвНе и ДУремг и наи~ большей ионностью связи.

Расположение минералов в ряд в соответствии с направлением перераспределения изоморфных ионов А и В при понижении температуры дает ряд перераспределения (Перчук. 1970). Поскольку такое перераспред еление определяется энтальпией обменных реакций, согласно соотношению <11пК/с1Т -ДН° /ИТ2 . а ДН°-ДНвА(МУ)-ДН°8А(МХ) . то в ряду перераспределения компонентов минералы расположены в порядке увеличения энтальпий-ного эффекта ДНВА. Чем больше удалены друг от друга члены этого ряда, тем больше ДН° и больше изменение константы равновесия при изменении температуры. Поскольку между ДСвд, ДНВЛ и ДУВА существует корреляция, то ряд перераспределения подобен ряду распределения, а объемный эффект ДУВд дает возможность определить положение минерала в этих рядах при условии одинаковой координации изоморфных атомов.

ДН, дв, кДк

Ш2

Мс

^Срх

а

■ МЬ

□о

" 11т

о

■ МС12

а

" МВГ2

J_|_

см3

-2

а

Рис. 3. Соотношение между объемным АУрвм? и энергетическими эффектами РеМв^замещения в октаэдрических позициях минералов в стандартных условиях (АНр-вМе ~ светлые и ДОреМ? ~ черные 1свадраты).

Глава 6. РАСПРЕДЕЛЕНИЕ КОМПОНЕНТОВ МЕЖДУ МИНЕРАЛАМИ И ТЕРМОДИНАМИЧЕСКИЕ УСЛОВИЯ МИНЕРАЛЬНЫХ РАВНОВЕСИЙ В БЕРДИЧЕВСКИХ ГРАНИТАХ И ДРУГИХ КОРДИЕРИТСОДЕРЖАЩИХ ПОРОДАХ УКРАИНСКОГО ЩИТА Исследование высокоглиноземистых нижнепротерозойских бер-дичевских гранитов, для которых характерна ассоциация Сг+ Сог+В1, которые занимают значительные площади в западной части УЩ и образовались в интервале от низов амфиболитовой фации до верхов гранулитовой фаций, дает возможность установить термодинамические условия петрогенеза и использовать полученные данные для проверки пригодности термодинамических моделей минералов, а также сопоставления и согласования геотермометров и калибровки новых.

При переходе от низкотемпературных гранитов к высокотемпературным уменьшается железистость граната, кордиерита и биотита. Для определения содержания летучих компонентов в кордиерите использован оптический метод, метод ПМР и масс-спектрометрический анализ. Характерно более высокое содержание СОг и более низкое Н20 в кордиеритах бердичевских гранитов (Н20 0.4-0.7 мае. %, . С02 0.8-1.0 мае. 50, чем в кордиеритах пегматитовых жил (Н20 0.5-0.8 мае. 56. С02 0.2-0.3 мае. 55).

С увеличением магнезиальности биотита увеличивается суммарная концентрация титаноксифлогопита и однотипных миналов, уменьшается суммарная концентрация истонита и однотипных с ним миналов при довольно постоянной суммарной концентрации флогопита и аннита (~50%). По данным ЯГР, отношение Рея + /(Ре) в биотите увеличивается от 3 - 49£ в наиболее высокотемпературных гранитах до 12% в низкотемпературных.

- 36 -

Анализ распределения Fe и Mg между минералами.

Термодинамические условия минеральных равновесий.

Температура. Применение калибровки GrCor геотермометра И.В.Лаврентьевой и Л.Л.Перчука (1981) указывает.на температурный интервал от 630 до 790"С минерального равновесия в бердичевских гранитах. Поскольку между калибровками GrBi геотермометра существуют значительные расхождения, проведена его калибровка по показаниям Gi-Cor геотермометра:

lnK'o (GrBi) - 1784/Т - 0,454 (г - 0.916) (6.1)

3 +

или с учетом присутствия Fe в биотите

InKo (GrBi) - 2177/Т - 0,772 (г - 0,965) . (6.2) где К' - (Fe/Mg)Gr/(Fe/Mg)Bl иК- (Fe/Mg^MFe^/Mg),,, . Полученные результаты близки к калибровке (Feriy and Spear, 1978).

Обнаружена температурная зависимость железистости граната. кордиерита и биотита в бердичевских гранитах:

t. °С - 1219 - 708 Fe/(Fe ♦ Mg)6r. (6.3)

t, °C - 1026 - 1092 Fe/(Fe + Mg)Cor. (6.4) t. °C - 1028 - 768 Fe/(Fe + Mg)Bi. (6.5) Использование нескольких согласованных показателей температуры делает более надежными результаты термометрии.

Давление определено с помощью предложенной калибровка геобарометра Сог+Gr+Si1+Q, которая учитывает зависимость давления равновесия от состава водно-углекислотного флюида.

Образование бердичевских гранитов происходило в значительном интервале давлений от ~500 МПа при 630°С до 670 МПа при 790°С (рис.4). РТ-соотношение в интервале 630-790°С выражается уравнениями

Р.МПа - -367 + 1.31 t°G ' (6.6)

Рис.4. РТ-условия равновесия Сог+(Зг±ЗШС1 в берди-чевских гранитах (черные квадраты и горизонтальные черточки для анго 0.2 и 1.0 соответственно) и гнейсах УЩ (светлые квадраты для анго 0.2).

- 38 -

Р.МПа - -149 + 1.06 1°С (6.7)

при анго 0.2 и 1.0 соответственно, что указывает на температурный градиент около 30°С/км. Различия в оценках давления, связанные с неопределенностью состава флюида, незначительны при температуре около 800°С и меньше 100 МПа при температуре около 600°С. Более высокотемпературные граниты образовались в более глубинных частях Земной коры, изотермы на площади выходов гранитов характеризуют степень поднятия отдельных участков и глубину эрозионного среза. На более тщательно изученном участке к югу от г.Бердичев они оконтуривают куполовидное поднятие. Флюидный режим. Геогигрометр СгВ1Рбр51Щ указывают на низкую активность Н20 (0.1-0.3 или 0.2-0.5 при использовании моделей биотита А1 и АЗ соответственно). Окислительно-восстановительные условия соответствовали интервалу от РЩ до

рма+2.

Содержание других компонентов в минералах как отражение термодинамических условий Кордиерит. Кордиериты с таким содержанием летучих компонентов, которое установлено в кордиеритах бердичевских гранитов, устойчивы при значительно меньшием давлении, чем то, на которое указывает СогСгёПО геобарометр. Это объясняется частичной потерей летучих компонентов кордиеритом на его пути к Земной поверхности. Вместе с тем состав летучих компонентов отражает меньшую долю СОр во флюиде при образовании пегматитов по сравнению с гранитами. РТ-условия устойчивости кордиерита с установленным содержанием НоО и С0р образуют границу, которую бердичевские граниты могли пересечь лишь при сравнительно низких температурах, предотвращающих даль-

- 39 -

нейшую дегазацию кордиерита по кинетическим причинам.

Титан в биотите. Между числом атомов TI в формульной единице биотита бердичевских гранитов и температурой GrCor равновесия существует зависимость (п-41)

ln[Ti] ± 0.03 - 3.360 - 4635/Т . (6.8) которую можно использовать для целей геотермометрии этих пород.

Содержание титана в биотите регулируется реакцией между компонентами минералов парагенезиса Gr+Bi+Fsp+Q+Ilm с участием титаноксианнита

TiAnn + l/3Sil + 2/3Q - i/3Alm + San + Ilm . (6.9) Выражение температурной зависимости кажущейся константы равновесия этой реакции сильно зависит от выбранной модели биотита. Применение регулярной модели AI дает большой разброс точек на диаграмме 1пК'-1/Т и приводит к выводу об отсутствии зависимости 1пК' от температуры. Напротив, применение ЛЭН модели биотита A3 позволяет получить уравнение геотермометра

In (aÄl„1/3asan/aTlAnn) - -9379/Т ♦ 11.817 ±0.38. (6.10) Трехвалентное железо в биотите. В характерном для гранитов парагенезисе Bi+Fsp+Mt содержание трехвалентного железа в биотите регулируется реакцией

Oan - San + Mt , (6.11)

в которой участвует оксианнитовый компонент. Модель A3 ЛЭН твердого раствора биотита использована для расчета температурной зависимости кажущейся константы равновесия по сведениям о составе минералов и температурах GrCor-равновесия (п-15)

InK'i 0.09 = -0.9223/Т + 13.46 . (6.12)

- 40 -

Применение регулярной модели биотита А1 дает зависимость 1пК'* 0.32 - -701В/Т + 16.558 , (6.11)

с большим разбросом точек на диаграмме 1пК - 1/Т. Таким образом, предпочтительность применения модели АЗ ЛЭН твердого раствора биотита подтверждается на природном материале.

Термометрия гр'анитоидов Украинского щита по составу биотита. Найденные соотношения между температурой и содержанием Ре3+ и Т1 в биотите бердичевских гранитов были применены для термометрии гранитоидов Украинского щита. Использовались усредненные химические характеристики биотита различных грани-тоидных .комплексов УЩ, определенные по опубликованным данным. Несмотря на значительную неточность такого расчета его результаты указывают на более высокие температуры минерального равновесия (>650°С) в гиперстенсодержащих гранитах гай-воронского, токмакского, букинского и восточно-приазовского комплексов, а также в гранитах коростенского и синюхинского комплексов по сравнению с гранитоидами других комплексов, у которых особенности состава биотита указывают на меньшую температуру.

Титансодержащий кварц как геотермометр. Использование такого сквозного минерала как кварц в качестве геотермометра значительно увеличивает возможности геологической термометрии. Установлена зависимость мольной доли Т102 в кварце от температуры Сг-Сог равновесия в бердичевских гранитах

1п Ят,02 ЧМ9 - - 6846/7 -1,528 . (6.14) Для реакции (6.15) РеТЮ3 + 1/3 А128105 + 2/35Ю2 - 1/3 Ре3А12813012 + (Т102)д. которая регулирует содержание титана в кварце в ассоциации Сг+БП+й+Пш , определена температурная зависимость кажущей-

ся константы равновесия при РТ-условиях образования берди-чевских гранитов

In (NI102Q NFe6r) - - 5248/Т - 3,484 , (6.16) а затем рассчитана РТ-зависимость кажущихся констант равнр-весия реакций с участием кварца (табл.2), которые можно ис-

Таблица 2

Кажущиеся константы равновесия реакций с участием компонентов титансодержащего кварца In К' - а/Т + ь + С(Р - 0,1)/Т

Реакция I а b с

Ilm + 1/3 Sil + 2/3 Q - I

1/3 Alm + (TiO2)0 I -5456 -3,060 -0,34 (TiCfe )Ru - (TI02)Q I -6029,7 -2,119 -1,12

Ilm. + Q - Fs + (TiOjJg I -6022.1 -2,934 0,763 пользовать для термометрии ильменит- и рутилсодержащих пород, в частности, чарнокитов и эндербитов.

Сравнение условий минеральных равновесий в металели-тах и бердичевских гранитах. Термобарометрия кордиеритсодержащих гнейсов Днестровс-ко-Бугского, а также других районов УЩ (Ингуло-Ингулецкий, Приазовский) указывает на то же соотношение между температурой и давлением (рис. ), какое характерно для минеральных равновесий в бердичевских гранитах и которое можно считать характеристикой древней геотермы. Поэтому все зависимости между составом минералов и температурой, установленные для бирдичевских гранитов, пригодны для термометрии кордиеритсодержащих гнейсов отмеченных районов УЩ.

- 42 -ЗАКЛЮЧЕНИЕ

1. Обнаружены закономерности в отклонениях энергетических и пространственных параметров минералов от аддитивности. Существует линейная или нелинейная зависимость между параметрами членов изоструктурного ряда АХ (Х-Хг,Хг,Х3....) с близким характером химической связи и аналогичных рядов ВХ. Существует корреляция между коэффициентами линейной зависимости. если у всех рассматриваемых соединений как катионы, так и анионы являются атомами элементов одной подгруппы или одного семейства. Зта связь между энергетическими и пространственными параметрами минералов позволяет с довольно высокой точностью определять их неизвестные значения , ликвидировать многие пробелы, провести ревизию имеющихся термодинамических констант.

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

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

4. Опираясь на связь между закономерностями размещения

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

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

6. Определены термодинамические условия образования бердичевских гранитов, выбранных в качестве полигона для испытаний термодинамических моделей минералов и калибровки геотермометров. Интервал температур и давлений (630-790"С и ~500-670 МПа) соответствует андалузит-силлиманитовому типу метаморфизма с температурным градиентом около 30"/км. обусловившим вертикальную зональность. Такие же РТ-соотношения и. следовательно, вертикальная зональность установлены для метаморфических толщ, содержащих гранат-кордиеритовые гнейсы, в разных районах Украинского щита.

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

вых показателей температуры, в частности, на основе растворимости Т1 в кварце и биотите, Реа+ в биотите.

ОСНОВНЫЕ ОПУБЛИКОВАННЫЕ РАБОТЫ ПО ТЕМЕ ДИССЕРТАЦИИ Монография

1. Курепин В.А. Термодинамика минералов переменного состава

и геологическая термобарометрия. Киев: Наук, думка, 1981,' 160 с. Статьи:

2. Курепин В.А. Распределение Ре и Mg между сосуществующими

силикатами в горных породах Украинского щита. Геохимия и рудообразование, 1973, N2 , 48-57.

3. Курепин В.А. Свободные энергии реакций в системах М10 -

М*0 - БЮ^ и связь между ними. Геохимия, 1973, N 1, 125-130.

4. Курепин В.А. Термодинамическая модель неупорядоченного

твердого раствора. Докл. АН СССР, 1975, т. 220, N 2. с. 462-464.

5. Курепин В.А. Мольные объемы соединений Мп, Ре,Со, № и

Са и связь между ними. Конституция и свойства минералов. 1975, N 9, с. 130-133.

6. Курепин В.А. Термодинамический анализ внутрикристалли-

ческого распределения катионов в простых шпинелях. Геохимия, 1975, N 6, 844 с.

7. Курепин В. А. Активность компонентов, термодинамические

характеристики реакций и фазовые равновесия в системе Ре-0 при высоких температурах и давлениях. Геохимия, 1975, N 10, 1475-1483.

8. Курепин В.А. К методике термодинамического анализа мине-

- 45 -

ральных равновесий с участием фаз переменного состава. Геохимия, 1976, N 2, 289-29В. 9. Курепин В.А. Термодинамическая модель гетеровалентного твердого раствора при условии локальной электронейтральности. Геохимия, 1978. N 1, 16-33.

10. Курепин В. А. Гранат-кордиеритовый парагенезис в горных

породах Украинского щита как показатель термодинамического режима метаморфизма. В сб.: Метаморфизм Украинского щита и его обрамления. Киев: Наук, думка, 1978, 44-45.

11. Курепин В.А. Термодинамический анализ минеральных равно-

весий в пироксенсодержащих горных породах. В кн.: Породообразующие пироксены Украинского щита. Киев, Наук, думка, 1979, 115-127.

12. Курепин В. А. Термодинамика водосодержащего кордиерита и

минеральных равновесий с его участием. Геохимия, 1979, N 1. 49-60.

13. Курепин В. А. К термодинамике окислов и силикатов железа

-при высоких температурах и давлениях. Геохимия, 1980. N3. 381-387.

14. Курепин В. А. Ортопироксен в железистых породах как гео-

барометр. Минерал, журн.. 1980, N 3, 63-69.

15. Курепин В. А., Полыпин З.В., Алибеков Г. И. Внутрикристал-

лическое распределение катионов ?ел* и А1 в клино-пироксене СаГея+А15Юл. Минерал.журнал, 1981, N 1, 83-88.

16. Курепин В. А. Термодинамика распределения компонентов

между породообразующими силикатами. Геохимия, 1982, N 1, 71-98

17. Курепин В.А. Термодинамическая модель гетеровалентного

твердого раствора с равновесным нарушением локальной электронейтральности. Геохимия, 1983, N 12, 1688-1698.

18. Курепин В. А. Содержание Нг0 и С0г в кордиерите как пока-

затель термодинамических условий минералообразова-ния. Геохимия. 1984, N 8. 1125-1134.

19. Курепин В. / , Польшин Э.В. Трехвалентное железо в биоти-

те и термодинамика гранат - биотитового равновесия. В сб.: "Термодинамика в геологии. I Всесоюзный симпозиум. Суздаль 1985". Черноголовка, 1985, т. 2, 51

20. Курепин В.А., Малюк Г.А., Калиниченко A.M.. Уточкин Д.В.

Летучие компоненты в кордиерите из бердичевских гранитов (Украинский щит). Минерал, «урн., 1986, т. 8. N 2, 70-82.

21. Курепин В. А. Термодинамические модели гетеровалентных

твердых растворов. Геохимия, 1986, N 10, 1402-1411.

22. Курепин В.А. К термодинамике гетеровалентных твердых

растворов в системе MgSi03 - AI0O3. Геохимия, 1986, N 11. 1606-1613

23. Курепин В. А. Термодинамика твердых растворов и устойчи-

вость биотита в системе Fe-Mg-KAlSi308-0-H. Геохимия. 1987, N 9, 1337-1346.

24. Курепин В.А. Термодинамика внутрикристаллического расп-

ределения катионов в простых шпинелях. Геохимия, 1988, N 5, 688-697.

25. Курепин В. А. Термодинамика распределения катионов и тер-

модинамические функции бинарных шпинелей. Геохимия, 1988, N 9. 1302-1312.

26. Курепин В.А. Термодинамика минералов и реакций между их

компонентами в системе РеО-А140з-5102-Н20-0г. Геохимия. 1990, N 7. 988-999.

27. Курепин В.А. 0 кристаплохимической формуле биотита и

концентрации его компонентов. Минерал, журн.. 1991. т. 13. N 2. с. 3-10.

28. Курепин В.А. Калибровка геобарометра кордиерит + гранат

+ силлиманит +■ кварц с учетом вхождения Ня0 и COg в кордиерит. Геохимия, 1991, N 2. 250-257.

29. Курепин В.А. Термодинамические условия образования гра-

нат - кордиерит - биотитовой ассоциации в берди-чевских гранитах (Украинский щит). Минер, журн.,

1991. т. 13, N 1. 76-87.

I -•

30. Курепин В.А. Титансодержащий кварц как высокотемператур-

ный геотермометр. Минер, журн.. 1992, N 3. 31-39.

31. Курепин В.А. Связь между энергетическими и пространс-

твенными эффектами полного изоморфного замещения. Геохимия, 1992, N 6, 787-799.

32. Kurepin V. А.• Correlation between energy and space

parameters of minerals. II Internat.Symp. "Thermodynamics of Nat. Processes". Novosibirsk.

1992. 33

. 33. Kurepin V.A. A correlation between energy and spatial parameters of minerals as a means for enlargement of thermodynamic database. IMA 16th General Meeting. Abstracts. Pisa. 1994. 47 .

-'48 -

КУРЕПИН В.А. Термодинамика минералов и геологическая термобарометрия.

Диссертация на соискание ученой степени доктора геологических наук по специальности 04.00.02 - геохимия. Институт геохимии, минералогии и рудообразования HAH Украины. Киев,1995.

Защищается 65 * печатных работ по исследованию связи между термодинамическими и <кристаллохимическими характеристиками минералов и оценке их термодинамических констант, термодинамике твердых растворов , термодинамике распределения компонентов между минералами и ее использованию для определения термодинамических условий минеральных равновесий.

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

KUREPIN V.A. Thermodynamics of minerals and geological thermobarometry.

Doctoral dissertation (Geology) on speciality 04.00.02 -Geochemistry. Institute of Geochemistry, Mineralogy and Ore Formation,Ukrainian National Academy of Sciences, Kiev,1995.

65 scientific works on investigations concerning relationships between thermodynamic and crystal-chemical characteristics of minerals and evaluation of their thermodynamic constants, thermodynamics of mineral solid solutions, thermodynamics of distribution of components

between minerals and Its application for determination of thermodynamic condition of natural mineral equilibria are defended. New relationships between energetic and spatial parameters of minerals of different isostructural series, between energetic and spatial effects of atomic substitutions in minerals have been revealed. Thermodynamic models of minerals including models of heterovalent solid solutions have been developed and applied for aims of geological thermobarométry.

Ключов1 слова: термодинамика м!нерал1в, геологi4Ha тер-мобарометр1я, термодинам1чт константи, тверд1 роэчжи, м1-неральн1 р!вноваги

Keywords: mineral thermodynamic constants.

thermodynamics, solid solutions.

geothermobarometry, mineral equilibria