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

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

.pdf
Скачиваний:
5
Добавлен:
19.10.2023
Размер:
5.43 Mб
Скачать

Обозначив через /га, содержание компонента к системе, по­ лучим из уравнений (69), (70) и (68), (71) выражения для определения состава паровой и жидкой фаз в системе через величину е\

 

 

 

 

Ъ =

1

 

е

 

 

 

(72)

 

 

 

 

 

к.

 

 

 

 

 

 

 

 

 

 

т\

«I

 

 

 

(73)

или через s =

1 е

 

1-г(кі— 1) ф

 

 

 

 

 

 

Щ

 

 

 

 

 

 

 

 

 

Уі=

 

 

 

 

(72')

 

 

 

 

 

1

1

- 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Хі — ■

гп,

 

 

 

 

(73')

 

 

 

 

+ (1 —

ЛГі) в

 

 

 

 

 

 

 

 

 

К,

 

 

 

 

Нетрудно

заметить,

что уравнения (72) и (73)

обращаются

в (72') и (73'), если в первых заменить е

на

1 — е.

 

Картина

испарения

и

конденсации, выявленная нами для

■бинарной смеси,

в основном

остается такой же и для много­

компонентных

смесей

при

 

других, конечно,

количественных

показателях

и при условии, что процессы

рассматриваются при

назких

давлениях, чтобы

 

можно было

для

них

применить

закон Рауля—Дальтона.

 

 

жидких и паровых фаз про­

Благодаря

различию состава

цессы испарения и конденсации можно проводить

так,

чтобы

добиться

почти

полного

разделения смеси на слагающие ее

компоненты. Так, например,

из пропанбутановой смеси

полу­

чают технически чистые пропан и бутан,

из водного раствора

спирта—спирт^крепостью 98° Си т. д. Способы

разделения сме­

сей рассматриваются, в специальном разделе („Ректификация“) химической термодйнамики.

8. ИСПАРЕНИЕ И КОНДЕНСАЦИЯ МНОГОКОЛ1ПОНЕНТНЫХ СМЕСЕЙ ПРИ ВЫСОКИХ ДАВЛЕНИЯХ

Как было указано выше, закон Рауля—Дальтона, выражаю­

щийся' формулой —

= к — , можно применять

т,олько при

х

^

Р

 

низких давлениях. Уже

при средних давлениях, по сравнению

•с критическими, этот

закон заметно нарушается

и требует

поправок. Отклонения эти сопровождаются отклонением урав­ нений состояния реальных газов от уравнений идеальных газов. Поэтому возникает мысль—нельзя ли внести такие изменения в

уравнения состояния реальных газов,

при которых формально

ру

при условии замены

•сохранился бы вид выражения к —

упругости

Ру и общего давления Р другими зависящими от

них параметрами. Эта

задача была разрешена Льюисом в раз­

работанной

им теории

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

9.ТЕОРИЯ ФУГИТИВНОСТИ

Известно, что характеристическая функция—изобарный по­ тенциал—выражается формулой

0 = i — TS,

полученной для изобаро-изотермического процесса, которому подчиняются фазовые превран^ения индивидуальных веществ. Мы видели, что в процессах испарения и конденсации смесей при Я=пост. температура t становится переменной и, очевидно,' при Т^пост. давление Р будет также переменным.

Для последующего можем написать

 

(^Ф \

=

\

_ г |AS\

 

Так как

\ дР Iт

 

\д Р

\дР / т

 

 

 

 

 

 

ді ) = AV - A T Z - и ( Щ

' дТ р

дР т

 

\

дТ р

дР

то после подстановки и сокращений получаем

 

 

 

 

!

= Аѵ.

(74)

Для идеального

газа

 

т

 

 

 

 

 

.0 = і -

TS = av Т -f-y- Tz—iQ—T (ap I n T + b ' T -

поэтому

-

AR ln P - S 0),

 

 

 

 

 

 

Из выражений (74) и (75)

получаем

 

 

 

/д\пР\

_ у

(76)

 

 

[

дР h

RT

 

 

 

Лыоис предложил сохранить-" вид этой формулы и для реальных газов при условии замены давления Р величиной /, которую он назвал „фугитцвностью“ или „летучестью“, написал

ln /

\ _ ѵ_

'(77)

[

дР

/т— RT

 

Пользуясь выражением

(77), величину фугитивности можно

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

нений состояния для реальных газов.

1

9 1

Однако решение их сложно и ввиду несовершенства Урав­

нений недостаточно надежно при высоком

давлении, прибли­

жающемся к критическому. Летучесть

таким

путем .вычисля­

ется для отдельных

углеводородов и в

литературе приводится

в виде диаграмм. Более популярным

является метод нахож­

дения летучести

при помощи

коэффициента сжимаемости в

обобщенной форме, пригодной

как для

отдельных углеводо­

родов, так и для смесей.

 

 

 

через г заменим в фор-

Для нахождения выражения для /

муле

ZR.T

Тогда получим

 

 

 

 

(77) V на —— •

 

 

 

 

 

 

 

(д \ n f \

_

г

,

 

 

 

откуда

 

V дР )т

 

Р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d ln / =

dP

 

 

 

 

и

 

 

 

 

Р

 

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

предполагается,

что

значения / 0 и Р0 взяты при низких:

давлениях. Так как

при

этом

реальный

газ стремится перей­

ти в идеальный, то величина / 0 не должна

отличаться от дав­

ления Р0. Считая

поэтому, что / 0= Р0,

можем

написать:

 

 

 

 

 

р

 

 

 

 

 

Ро

Произведя замену Z на (1 — Z), будем иметь.

рр

І п Х

_ С dp

V ~ * d f i

 

In Р

Ро

1 р

J р

 

Ро

откуда

Ро

Ро

 

 

 

 

Д•

II >|с-ч

 

 

 

 

1

 

где

 

 

 

 

 

 

 

Л =

f

 

d«.

 

 

' . А*

 

.Тогда величина / выразится через

/= P-e~fs.

'Для определения величины / можно применить ту же диа­ грамму, что и для определения f s. В этой диаграмме по оси

ординат кроме шкалы для f s имеется и шкала величины т,.

92' \

определяющая отношение

Поэтому, сняв с диаграммы для

заданных параметров

те ,и т

значения

7,

величину / следует

определять из соотношения >

 

 

 

 

/ = ТЯ.

 

 

Для газов, мало отличающихся от

идеальных, величину /

можно найти более просто. Для них

 

 

Z — 1 — 0 Я, где 0

•—

=/(Г).

Поэтому напишем

 

 

я г

 

 

 

 

 

р

 

 

 

 

l n ^ = f 1 ЬР d P ln Р

- Ѳ ( Я - Я о ) ,

Po J

р

Ро

 

 

Р о

 

 

 

 

откуда

 

 

 

 

- 1ПХ = _ 0 ( Я - Я 0) = Д - Д 0,

Потенцируя, получаем

_/ = ег-г„ _

Р

Так как в нашем случае степень Z — Z0 мало отличается от 0 ;

то, разлагая величину ег ~ г0 в ряд и пренебрегая членами раз­ ложения, начиная с третьего, получаем

 

ем ° = 1 + ( Z - Z 0),

j

Следовательно,

для газов, мало отличающихся от идеал ь

ных,

(

Г = { \ + ( г - г ^ ) ] р

 

 

 

 

или ввиду малого отличия величины Z0, взятой при давлении

Я0, от

1,

f = ZP.

.

 

 

 

 

Итак, по Льюису, замена. давления фугитивностью дает возможность внести поправку и в закон Рауля—Дальтона.

Поправка эта заключается в том, что в равенстве к = ^

значение константы к определяется формулой

где / ж—фугитивность жидкости;

 

/ п—фугитирность

пара,

по известным правилам,

Обе эти величины

определяются

из которых опишем два наиболее простых.

1. Для нахождения величины./с, любого компонента систе­

мы, находящейся при 'давлении Я и

температуре Т, .опреде-

93

 

 

Р

 

 

т

 

 

Рк1 и Г,іі—его

критические па-

ляют -я:, = —— и

*! = — , где

 

 

 

Ркі

 

0.-І

 

 

 

(см. приложение

7)

нахо-

рамет.ры. Затем

по

диаграмме

дят вёличину fi, и из равенства f

 

/?

вычисляют

фугцтив-

= —

ность

пара /" — ^

Р

при давлении

Р.

 

насыщения Р\, отве-

Найдя по кривой упругости давление

чающее

температуре

Т

и

 

 

 

 

 

 

 

рУ

 

 

соответствующее uf = — , вновь

обращаются к диаграмме

(приложение № 7)

Ркі

т, на­

и по іг[

и

ходят

fi , а

затем

из

равенства

 

 

/,

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

жидк

сти / г

 

 

 

 

 

 

 

 

 

 

 

Р\

 

 

 

 

 

=

rt-P 3!

при давлении

PJ. Тогда

константу

найдут

из соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

•л-

-

 

 

 

 

 

 

(а.

Этот способ считается

приемлемым

во всех случаях, если

давление Р\

насыщения выше

давления

Р системы.

Если Ру\

ниже

Р,

то

в расчет фугитивности жидкости

вносят поправку,

а именно: сначала находят фугитивность

жидкости f t описан­

ным способом при давлении Р\, а затем из выражения

 

 

 

 

 

 

In (/а

 

 

 

 

і

 

 

 

 

(б>

 

 

 

 

 

 

( П

 

 

 

RT

 

1 dP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

РУ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определить

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

 

при давлении Р.

 

 

В

этом

выражении ѵ-,—удельнуй

объем

жидкости, мало*

изменяющийся при изменении давления и при Т = пост.

Предположив поэтому,

что ^| = пост.,

получаем:

 

 

 

 

 

 

ІП

(->«

 

 

— (Я-*Я?)

 

 

 

(в)

 

 

 

 

 

(/Г).!

 

 

RT 4

 

 

1

 

 

 

 

 

 

 

 

 

 

Рі

 

 

 

 

 

 

 

 

 

 

 

или при

замене і\

на г, RT (г'і—средний

коэффициент сжимае­

мости

жидкости компонента

в пределах

 

изменения

давления,

от P t

до

Р)

будем

иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(/Г)р

_

/ Р _ у

 

 

 

 

 

 

 

 

 

 

 

 

(/Г)р

 

' Uy)

 

 

 

 

(р)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Практика

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

выражающихся десятками атмосфер,

правая часть выражения

94

(в) будет мала,- а выражение (г)—незначительно отличается от

единицы. Тогда

фугитивности (/Г)Р и ( / f )р[ получаются

блгз-

кими друг к другу и указанной поправки

часто не делают.

2. Для нахождения фугитивности жидкости для

углеводо­

родов Ашворт построил

диаграмму (приложение № 12), где в

зависимости

от

температуры

t

системы и

числа углеродных,

атомов в молекуле данного

углеводорода непосредственно

снимается (на абсциссе)

величина фугитивности жидкости ft-

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

пара / “ по-прежнему

находят,

как

было

описано, *с помощью диаграммы

(приложение № 6) по

пара-

метрам к, =

Р

т

(Рк1 и Тк\ критические

параметры

-—- и т, = —

 

к!

Т'к,

 

 

 

 

 

 

данного компонента, а А—давление системы). Наконец, /г, оп­ ределяется как к, = ff

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

10. ДИАГРАММЫ КОНСТАНТ РАСПРЕДЕЛЕНИЯ ДЛЯ КОМПОНЕН­ ТОВ СИСТЕМЫ, ОБЛАДАЮЩИХ ВЗАИМНОЙ РАСТВОРИМОСТЬЮ

Для нахождения

констант распределения

широко исполь­

зуют диаграммы, с

которых непосредственно

можно снимать

при любом давлении Р

и температуре t системы величины к.

_

 

РУ

и напишем ее

в виде

Возьмем зависимость к = —

 

 

 

р

 

 

 

 

 

*

lg к =

lg Ру — lg Р.

 

(д)-

Если по оси абсцисс

откладывать lg Р,

а

по о$и ординат

lg к, то формула представит собой уравнение

прямой или, точ­

нее сказать,—семейства прямых,

параметром у которых служит-

величина упругости,

зависящая

от рода углеводорода и его

температуры t. Угловой коэффициент прямых равен единице.

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

Для какой-нибудь из этих прямых ав, как это следует из уравнения, отрезки ас и cb должны представлять собой вели­ чины lg Ру.

 

Рис. 37. Зависимость

констант углеводородов

 

 

 

 

 

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

и давления:

 

 

 

 

 

 

реальные ^величины;------ — величины, вычи­

 

 

 

 

сленные по закону Рауля—Дальтона; і----------- вычисленные по

 

 

 

 

фугнтнвностям

 

 

 

 

 

 

 

 

 

 

 

Такими были бы

диаграммы,

 

 

 

если системы, в которые входит

 

 

 

данный компонент,

подчинялись

 

 

 

бы закону Рауля—Дальтона.

 

 

 

я-,,

Опытное уточнение величины

 

 

 

основанное

на

определении

 

 

 

для компонентов

величин у, и л:,

 

 

 

и

нахождении

их

 

отношения

 

 

 

к, = ■

показывает,

что

вели­

 

 

 

 

Х\

 

 

 

 

 

кривыми,

 

 

 

чины lg к, выражаются

 

 

 

которые лишь при

низких

дав­

 

 

 

лениях почти совпадают с наши­

 

 

 

ми

прямыми

(сплошн'ые

линии

Рис.

3S. Зависимость

констант

на рис. 37 и 38).

При высоких

давлениях

кривые

значительно

равновесия углеводородов от

 

давления:

 

изгибаются и-при гіекотором дав­

■■■■-

опытные величины;

— — — вели­

лении Рс сходятся в одной точке

чины, вычисленные по закону Р а у л я -

(/6=1). Это

давление

Рс носит

 

Дальтона

 

 

 

 

название

давления

схождения.

По существу оно является критическим давлением смеси тако­ го состава, при котором данная температура t является кри­ тической.

I •

96

Наши кривые не случайно сходятся в точке, имеющей зна­ чение lg к = 0 и к = 1. Дело в том, что в критическом состоя­ нии системы для любого ее компонента не существует разли­ чия между паровой и жидкой фазами. Следовательно, в этом состоянии у = X или к = 1.

По опытным данным обычно строятся диаграммы, как на pile. 38. Но в практике удобно пользоваться диаграммами такого типа, как на рис. 37, полученными в результате обра­ ботки предыдущих данных для разных систем. Очевидно, каж-• дая кривая на рис. 37 будет отвечать критической темпера­ туре некоторой системы, включающей данный компонент.

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

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

Подобные диаграммы констант равновесия для нормальных парафинов от метана до декана опубликованы американской фирмой NGAA для ряда давлений схождения.

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

давления Рс на величины к для

углеводородов, кроме метана,

можно пренебречь.

!

При высоких давлениях, как видно из диаграмм, для нахож­ дения констант равновесия необходимо знать давление схож­ дения данной системы. Если система бинарная, то давление Рс может быть определено по огибающим фазовых диаграмм (см. приложение 13), соответственно температуре смеси. На этом же графике—зависимость tKуглеводородов от молекулярного веса JA.

Если же смесь многокомпонентная, то"для определения

величины Рс

необходимо знать параметры

хотя бы одного ее

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

состояния при высоком давлении, например, дав­

ление

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

начала конденсации

(точки росы) или

кипения (точки кипения) смеси.

 

 

 

При этом условий давление Яс

находят методом

подбора

для удовлетворения

равенств V — = 1 (если известны P u t

точки

росы)

или

 

КI

 

 

S Кі хі 1 (если

известны Р и t точки кипе­

ния).

Константы агі,

входящие в последние выражения,

опреде-

Ф

а - 7

9 7

ляются для каждого компонента смеси по диаграмме констант, соответствующей принятому давлению схождения при Р и і смеси, а уі или Х\— принимаются 'равными мольной концентра­ ции данного компонента в смеси.

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

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

Мы видели, что теория фугитивности дает средство для корректировки закона Рауля—Дальтона. Поэтому интересно

сопоставить результаты вычислений констант к через завися- /Ж

мостъ к — ~ с полученными опытными кривыми. На рис. 37

кривой ав'е со штрихпунктирным окончанием показаны резуль­ таты вычислений констант через фугитивность. Как видим, этот

метод действительно

корректирует

прямые а в е

при

низких

в

средь их

давлениях, при

е ы с о к и х

.же Еблизи

критической

 

 

 

 

 

в

кри­

точки он дает неправильные результты и не приводит

тической точке значения к к единице.

 

 

 

Из сказанного следует, что при

состоянии систем —вблизи

критической точки—единственным пока способом

нахождения

констант к

является

экспериментальный.

 

 

 

 

11. ФАЗОВЫЕ ДИАГРАММЫ ГОМОГЕННЫХ СИСТЕМ

 

 

 

Возьмем

бинарную систему любого состава и го вышеизло

женной методике, задавшись давлением, определяем точку начала кипения а и точку росы в, которые расположатся в координатах РТ на одной горизонтали (рис. 39).

При изменении давления и после необходимых расчетов видим, что эти точки будут перемещаться по кривым а —а', в—в', если расчет производить по закону Рауля—Дальтона.

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

Линии точек кипения и росы (сплошные линии) сходятся в критической точке к системы с координатами Рк и Тк, образуя петлеобразную кривую, называемую ортобарной кривой. Эту кривую также часто называют фазовой диаграммой или Р —Т- диаграммой. По левую сторонуѵэтой кривой ак располагается область жидкой фазы системы, по правую вк—область газовой фазы, а в серединеі кривой будем иметь двухфазное паро­ жидкостное состояние системы.

98

Если при данном давлении Я будем стклэдыгать на рис. 39 равновесные темпера-туры t, отвечающие различным степеням конденсации е, то на прямой ав получим ряд точек, обозна­ чающих степень_влажности системы е—0,1; 0,2; 0,3...

Рис. 39. Фазовая диаграмма в.

Рис. 39а. Изохоры смеси пропан-н-

коордниатах Р, Т

гептана (вес. состав gj—-0,75)

Так как изменяемость t с величинойе, кск зто видно из рис. 35 и 36, происходит почти по прямой линии, то степень влажности б при данном давлении Я достаточно точно можно

определить по рис. 39—отношение отрезков е == по край­

ней мере в том интервале давлений, для которого можно при­ нять справёдливым заюн Рауля—Дальтона.

Соединяя одноименные точки влажностей, найденные для различных давлений Р, можно построить на рис. 39 семейство кривых постоянной влажности (изоплеры), которые сойдутся в критической точке системы (см. внутреннее построение между пограничными кривыми внутри „петли“).

'На Р—Я-диаграмме также можно показать изменение удель­

ного объема системы построением

ряда

нзохор

(вид

изохор

см. 39а).

1

 

 

 

-)§

На рис. 39а, п редставлякщем

Я— Я-диаграмму смеси про-

пан-н-гептана с

весовым содержанием

пропана 75%

(пропан

98,5%-ный), видно характерное

различие в формах изохор

жидкой и газовой фаз, переходящих в область

двухфазного

состояния.

 

 

 

 

 

Изохоры в области жидкого состояния (на рис. 39а изохоры /) более крутые, чем в области газообразного состояния (на

7 *

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