Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Рафальский Р.П. Гидротермальные равновесия и процессы минералообразования

.pdf
Скачиваний:
11
Добавлен:
24.10.2023
Размер:
10.76 Mб
Скачать

где Ро — давление чистого насыщенного водяного пара. Разность Р—Р0 условно принимается равной парциаль­ ному давлению газа. Как было показано С. Д. Малининым [70], изучавшим растворимость СОг в воде, в об­

ласти сравнительно

невысоких

значений

Р—Р0

(менее

100 атм) разница между

Kv и k для СОг

возрастает с

понижением температуры

(300—100°С), однако

вплоть

до 100° С остается

достаточно

малой, чтобы ею

можно

было пренебречь в геохимических расчетах. Что касается хорошо растворимых в воде полярных газов, то здесь расхождения между КР и k заметно больше. Так, из экспериментальных данных, характеризующих раствори­ мость в воде сероводорода [71] следует, что расхожде­ ния становятся значительными уже при величинах Р— —Ро, измеряемых первыми десятками атмосфер.

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

ления

фугитивностей газа,

вычисленных

по

соотноше­

нию

(2.2), с парциальными

давлениями.

При

расчете

фугитивностей использовали константы Генри для пре­ дельно разбавленных растворов и экспериментально найденные концентрации газа в жидкой фазе. Парци­ альные давления газа (для тех же экспериментальных условий) приняли равными Р—Р0.

Экспериментальные данные С. Д. Малинина [72J свидетельствуют о том, что при температуре 200° С фугитивности углекислого газа, равной 100 атм, соответ­ ствует парциальное давление 120 атм. Разница между парциальными давлениями и фугитивностями сернисто­

го газа и

сероводорода

при

/=20

атм составляет

соот­

ветственно

2,5 атм (180°С)

и 0,5

атм

(17ГС) .

Это

очень приблизительные

оценки, однако,

вероятно,

они

дают правильное представление о возможных погрешно­

стях, обусловленных

допущением f—p. Из

общих сооб­

ражений следует, что при

температурах,

близких к

200° С, расхождение

между

фугитивностью

и парциаль­

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

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

Таблица 14

 

 

Давление

[2] и

фугитивность

[73]

водяного

пара

 

 

 

на

ПИНИИ насыщения

 

 

t,

°С

100

150

200

250

300

350

р,

апш

1,0

4,7

15,3

39,2

84,8

163,2

f,

атм

1,0

4,5

14,0

33,5

66,1

113,8

давление газа невелико, его коэффициент

фугитивности

может, по-видимому, заметно отличаться

от

единицы,

что является результатом межмолекулярного

взаимо­

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

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

Имеющиеся в литературе сведения относительно ра­ створимости газов в воде при повышенных температурах и давлениях, охватывают все газы, представляющие наи­ больший интерес для геохимии, в том числе СОг, H2S, SO2, Н2 и Ог. Кроме того, опубликованы данные по ра­ створимости ряда других газов [74, 75]. Все это исклю­ чает необходимость далеко идущих экстраполяции кон­ стант равновесия, подобных сделанным в гл. 1.

Температурная зависимость.константы Генри анало­ гична температурной зависимости констант диссоциации слабых электролитов. В большинстве случаев константа

Генри возрастает по мере повышения температуры, про­ ходит через максимум, которому соответствует мини­ мальная растворимость газа в воде, а затем понижается. /Максимальное значение k для газов, растворимость ко­

торых

изучена

при температурах

выше

100° С, нахо­

дится

в интервале

75—200° С. Исключение

составляет

 

 

 

 

 

 

Таблица 15

Коэффициенты в

уравнении

(1.4) температурной

зависимости

 

константы Генри, агл-1000 г/моль

 

Газ

Температур­

А

с

D

Литература

ный интервал,

 

°С

 

 

 

 

 

 

 

 

 

 

 

*

о а

15—343

 

1202

9,69

0,009049

j [12]

н 2

24—343

 

1547

11,63

0,01140

с о ,

0—350

 

1825,5

10,348

0,0093003

H 2 S

10—330

 

1090,7

6,2J3

0,005057

[76]

so .

0—240

 

2065

9,13

0,00764

гелий, отличающийся

самой

низкой

растворимостью,

константа

Генри для которого равномерно

понижается

с ростом

температуры.

 

 

 

 

 

Зависимость

константы

Генри

от

температуры по

крайней мере до 300 °С удовлетворительно

выражается

эмпирическим

уравнением

(1.4).

И. Л. Ходаковский

и др. [12], обработав

критически

отобранные экспери­

ментальные

данные,

опубликованные

в литературе,

определили

числовые

значения коэффициентов в этом

уравнении для ряда газов;

некоторые из них приведены

в табл. 15.

Коэффициенты

в уравнении

температурной

зависимости константы Генри для S0 2 были вычислены путем обработки экспериментальных данных по системе S 0 2 — Н 2 0 [76].

РАСТВОРИМОСТЬ ГАЗОВ В ВОДНО-СОЛЕВЫХ РАСТВОРАХ

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

\ g ^ = Kmm,

(2.4)

где т'

и тТ

— растворимость

газа в воде

и

растворе

электролита

соответственно;

Кт

— коэффициент

Сечено­

ва; т—-моляльная

концентрация электролита.

Подоб­

ное уравнение, в котором вместо моляльности

подстав­

лена ионная сила, имеет вид:

 

 

 

 

 

 

l g —

- / Q i ,

 

(2.5)

где К», —коэффициент высаливания.

 

 

В

первом

приближении

можно принять

(пренебре­

гая взаимодействием между

нейтральными

молекулами

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

стях

в

газовой

фазе

активности газа,

растворенного

в чистой

воде

и растворе электролита,

равны,

имеем

[с учетом выражения

(1.38)]

 

 

 

 

 

 

ЫУ = К^,

 

(2.6)

где

у — коэффициент

активности газа в

растворе

элек­

тролита.

Соотношение (2.6) можно получить из двухпарамет-

рического

уравнения

Дебая — Гюккеля

(уравнение

(1.43)), в

котором первый член равен нулю, так

как

для находящегося в растворе неэлектролита

z + = z _ = 0 .

Результаты многочисленных экспериментальных

ис­

следований, выполненных при низких температурах, сви­ детельствуют о том, что уравнения (2.4)—(2.6) выпол­ няются при концентрациях электролита в растворе по крайней мере до 2,0 М [77], а в ряде случаев и до бо­ лее высоких концентраций. Для повышенных темпера­ тур данных значительно меньше, однако они также по­ зволяют считать, что приведенные уравнения справедли­ вы до указанной концентрации.

Значения коэффициента высаливания меняются в до­ вольно широких пределах в зависимости от составов неэлектролита и находящегося в растворе электролита. По данным, характеризующим растворимость 16 неэлек­

тролитов в растворах 1-1-электролитов

(в том числе кис­

лот

и щелочей)

при 15—30° С,

величина

Кц

меняется

от

0,255 до —0,13

[77]. Средняя

величина

Лц,

выведен­

ная из 151 значения 7<ц> полученного

для

конкретных

растворов, равна

0,074. Коэффициент

высаливания для

водорода, кислорода и углекислого газа в растворах солей щелочных металлов меняется от 0,141 до 0,038 (среднее из 19 значений для указанных растворов 0,086*). Коэффициент высаливания для сероводорода, вычисленный по результатам определения растворимо­ сти H2S в растворах 17 электролитов при 25° С [78, 79], меняется от —0,008 до 0,082 (среднее значение 0,035).

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

/Ср., по данным для растворов SO2 в присутствии

одного

из восьми 1-1-электролитов, меняется от 0,01 до

—0,13

(среднее —0,051).

Имеющиеся данные свидетельствуют о том, что ко­ эффициент высаливания не зависит существенно от дав­ ления при давлениях до нескольких сот атмосфер [70].

Что касается

зависимости от температуры, то, по дан­

ным Эллиса

и Голдинга

[80J, изучавших растворимость

углекислого

газа в растворах

NaCl, величина /Ср. ме­

няется

от 0,094 при 50° С до 0,166 при 300° С и до 0,313

при 350° С проходя через

минимум приблизительно при

150° С

(/Ср. =0,073). К выводу

о наличии минимума на

кривой

температурной

зависимости Кц для растворов

СОг в присутствии СаС!2

приходит также С. Д. Малинин

[64]. Вместе

с тем, Брюну

с соавт. [81], изучавшим

растворимость кислорода и водорода в растворах элек­

тролитов

при

температурах до 250° С, не удалось вы­

явить

зависимость

коэффициента

высаливания

от

температуры.

Последнее указывает

на то, что

если

изменение величины

Ку. и происходит,

то только

в пре­

делах погрешности эксперимента.

 

 

 

В настоящей

работе принято, что коэффициент

выса­

ливания

не зависит

от температуры и солевого состава

раствора и равен для кислорода, водорода и углекисло­ го газа 0,09, для сероводорода 0,07. Эти значения, повидимому, наиболее близко соответствуют средним зна­

чениям Кц при повышенных температурах.

 

На рис. 7 приведены коэффициенты

активности, вы­

численные нами

по уравнению (2.6) на основании

при­

нятых значений

К\х и коэффициенты

активности

угле-

* Получено без учета данных для растворов NaOH и КОН, не представляющих для нас большого интереса, но характеризующихся аномально высокими значениями К^.

кислого

газа [70, 80J,

водорода

и

кислорода

[81J, а

также

сероводорода [8J, основанные

на

эксперимен­

тальных определениях

растворимости

этих

газов

в

вод­

но-солевых растворах.

Как видно

на

графиках, при

низ-

о - о,г о& о,б о,8 р.

Рис. 7. Расчетные и экспериментальные значения коэффициентов активности газов в водных растворах электролитов.

ких ионных силах коэффициенты активности можно при­ нимать равными единице (при ц.=0,5 и t= 100—200° С у=£С1,1). С другой стороны, даже при ионной силе, рав­ ной 1,0, расхождения между экспериментальными и рас­ четными значениями коэффициента активности при 100 и 200°С настолько невелики, что можно использовать в

б Р, П. Рафальскнй

65

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

Данные по растворимости SO2 в растворах электро­ литов при повышенных температурах, по-видимому, от­

сутствуют. Коэффициент

активности

сернистого

газа

при ц.= 1,0, соответствующий среднему

значению

Ку. =

=—0,05 при 25° С, равен

0,89, т. е. довольно близок к

1,0. В дальнейшем принимается, что коэффициент выса­

ливания

для

S02

при

температурах,

превышающих

100° С,

равен

нулю,

а коэффициент

активности в

раство­

ре равен

единице.

 

 

 

 

 

 

РАВНОВЕСИЯ В СИСТЕМАХ С 0 2

Н 2 0

И H2 S — Н 2

0

На

примере системы

С 0 2 — Н 2 0

более подробно рас­

смотрим методы расчета равновесий в двухфазных си­ стемах, состоящих из водного раствора и газовой фазы. Равновесные активности компонентов в растворе и фу­

гитивность

углекислоты в

газовой

фазе

в

этой системе

определяются следующими

реакциями:

 

 

 

С02 (Г ) + Н 3 0 ( ж ) = Н2 СОз(р.р ) ,

йС о2 =

 

;

(2.7)

HiC03(p.p) =

Н^;.Р) -f- НСОз(р.р),

/Сн2соа =

 

а н2 со,

;(2.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

°н+ а с о ? -

 

Н С О І 7 Р . Р ) = Н + Р ) + С О | Г Р - Р ) , Кнгсо3=—

 

', (2.9)

Н А ж )

= Н£.р) + ОН £ . р ) ,

Kw

= a H

+

a 0 H - ,

(2.10)

где ан,со, —сумма

активностей

С 0 2 -р)

 

и Н 2 С 0 3 ( р . р ) .

Допуская,

что

система

открыта

относительно

С 0 2 ,

вычислим активности частиц в растворе при различных заданных фугитивностях С02 ( Г ). При этом условии четы­

ре выражения

констант равновесия реакций (2.7)—•

(2.10) содержат

пять неизвестных. Пятым уравнением,

необходимым для определения последних, служит урав­ нение баланса.

По достижении равновесия ионы в растворе присут­ ствуют в следующих количествах:

 

п н с о - = * — У'

п

н +

^ х

+ У + 2''

 

 

где х — количества

ионов

HCOj- и Н+

образовавшихся

по реакции

(2.8); у — количество ионов

НСО^~, израсхо­

дованных в

реакции

(2.9),

и

равное количествам

ионов

Н+ и СО|~~ , выделившихся

в

результате этой

реакции;

г — количества

ионов

Н+ и

О Н - ,

выделившихся

при

диссоциации

воды.

 

 

 

 

 

 

 

 

 

Путем совместного решения последних четырех урав­

нений получаем уравнение

баланса

 

 

 

 

 

"н+ =

"нсо- +

2 " с о 2 -

+ " о н - •

 

(2-1 1 )

Необходимо

иметь

в

3

что

3

уравнениях

баланса,

виду,

в

подобных уравнению

(2.11), участвуют

не концентрации

компонентов,

а их

количества,

выраженные

в

молях.

В том случае, если все компоненты, участвующие в урав­ нении баланса, находятся в растворе, мольные количе­

ства равны концентрациям*, которые, в свою

очередь,

при

достаточно малых значениях ионной силы

раство­

ра

можно

приравнять к

активностям. Более подробно

эти вопросы рассматриваются

ниже.

 

 

Из выражений констант равновесия реакций

(2.7) —

(2.10) следует:

 

 

 

 

 

 

 

 

 

 

а

_ =

ан+

;

(2.12)

 

 

 

он

 

4

'

 

 

 

 

-

 

 

 

 

 

% 2 сс, =

7к^ ;

(2-13)

 

 

 

 

 

kco2

 

 

 

 

_

3 со3

а н,со,

 

2 со, fco2 .

/о u \

 

 

а н с о Г - — —

 

 

- 1 ^ 7 '

 

( 2 Л 4 )

 

а с о 2 - - ^ С О з а Н С ° Г

=

* H ^ . W C O . .

 

( 2 Л 5 )

 

.

3 .

а н+

 

 

 

*соа .а н+

 

 

* Уравнение баланса при этом сохраняется неизменным при вы­ ражении концентраций как в молярной, так и в моляльной шкале [см. уравнение (1.1)].

Подставляя эти соотношения в уравнение (2.11), после небольших преобразований получаем кубическое урав­ нение

 

H.COafcO.

, J,

\

 

 

 

 

 

(2.16)

решая которое легко найти активность

водородных

ионов. Подставляя

полученные

значения а н +

в выраже­

ния (2.12) — (2.15),

определяем

активность

остальных

компонентов.

 

 

 

 

Как следует из

уравнений (2.8)

и (2.9),

присутствие

углекислоты (так же, как сероводорода) должно приво­

дить к понижению рН водных растворов.

Зависимость

рН раствора от фугитивности углекислого

газа в

систе­

ме СОг—НгО (в отсутствие других компонентов)

пока­

зана на рис. 8, где нанесены также значения рН нейт­ ральной точки при различных температурах: с повыше­ нием температуры влияние углекислоты на рН раствора существенно ослабевает.

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

ев, однако,

подобные

уравнения

можно

значительно

упростить,

пренебрегая

активностями (концентрациями)

второстепенных частиц

[82].

Например,

в

рассмотрен­

ном случае

К^сор> Л'н.соа'

благодаря

чему « н с о - >

^#со2—, а

активность водородных

ионов,

выделившихся

по реакции (2.9), несравненно меньше ан+, обусловлен­ ной реакцией (2.8). Что касается реакции (2.10), то ее вклад в суммарную активность водородных ионов ста­

новится

заметным лишь

при относительно

низких /со,

(около

0,1 атм). Таким образом, при / с о 2 ^ 0,1 атм мож­

но пренебречь образованием ионов

Н+

по

уравнениям

(2.9) и

(2.10). Уравнение

баланса

при

этом предельно

упрощается:

 

 

 

 

п'н+

(2.17)

 

а выражение для вычисления активности водородных ионов приобретает вид:

К н 3 с о 3 ' С 0 3

(2.18)

fecoa

 

Если расчеты относятся к растворам, ионная сила которых не настолько мала, чтобы можно было принять

lg fgjanfM)

Рис. 8. Зависимость рН раствора в системе СОг—Н»0 от фугитивности СОг (вертикальные линии соответствуют рН нейтральной точки).

активности ионов равными концентрациям, как это было сделано в рассмотренном случае, и если вместе с тем используются термодинамические константы равновесия реакций (а не концентрационные константы, относя­ щиеся к данной ионной силе), в уравнениях баланса необходимо учесть коэффициенты активности частиц. Перевод концентраций компонентов в системе СОг—Н2 0 в активности в соответствии с уравнением (1.38) дает

Соседние файлы в папке книги из ГПНТБ