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

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

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

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

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

теля

л с ѵ т х

, разность {Рх— Р2) была бы велика, а паде-

 

 

ние температуры Д Т заметно. ,

II. Пусть в трубе с протекающим через нее слева направо

газом (рис. 18 б)

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

диафрагма

с отверстием,

пористая перегородка), снижающее

давление

газа

с р х до р2,

и пусть остальные параметры газа

будут начальными: удельный объем г,, температура Тх, ско­

рость сх, и

конечными:

удельный объем ‘іо, температура Т2>

скорость с2.

 

 

Прилагая

к процессу

первый закон термодинамики

 

 

с2 — с2

 

= и<уих-f- AI А -----------,

 

 

2g

заметим, что в данном случае при течении совершается внеш­ няя работа перемещения по трубе газа. Найдем величину этой работы.

Предположим, что в начальный момент некоторое количе­ ство газа помещалось между сечениями трубы а —а и б—б и

что в этих

сечениях

находились

воображаемые

поршни, и

пусть по прошествии

времени

Ді эти поршни передвинулись

в положение а'а' и б'б', пройдя пути Дх, и Дх2.

давление.на

Так как

отсеченная часть

газа производит

левый поршень против движения

газа, а на правый—по дви­

жению, то левый поршень поглотит работу, правый же произ­ ведет'ее. Если / —сечение трубы, то эта работа выразится как

—f p xb.xx и + //? 2 Дх2 и будет отнесена не к 1 кГ,

а к. неко­

торому количеству G кГ газа. И

работа /,

которую

необхо­

димо ввести в выражение первого закона,

выразится так

^ _ fPx Ахх

 

fp2 Ьх2

"

 

G

'

G

 

 

Найдем величину G. Так как за отрезок времени Д? через сечения а —а и б—б протекает сбъем газа, j c xДт и /с 2Дт и при установившемся течении вес газа, проходящего через эти

» «_... -

50

сечения, должен быть одинаковым, то вес G газа,( производя­ щего работу, найдется так

Q _fcxАх _ / С о А X

Тогда

I/,

 

ѵ2

Ді Af.oi

Дз A >rat>3 _

_

Но

Ci Ax

 

c3 Ax

 

A jr3

 

А Л1!

C[ и

2 »

ТГ

Ax

поэтому

 

 

р24%ш'

l 3= _ plVl

Вновь обращаясь к выражению первого закона, видим, что из-за отсутствия при данном процессе подвода тепла извне мы должны считать, что q — 0, вследствие чего

 

 

 

 

 

AI

2

2

 

 

 

 

'

 

( и 2 — И ,) +

Со

Сі

 

 

 

 

»

 

2

1 _ - о

 

 

 

или

 

 

 

 

2g

 

 

 

 

и2—их—А р хтх-\-А Рг^'і +

с \ - с \

 

 

0

Но

А 22g

1-

 

 

 

А р21\ і2 и

иі + А р хѵ х

=

і х

 

 

 

+

 

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

окончательно получаем

 

 

<

 

«2

 

=

 

 

 

 

 

При малых скоростях сх и с2, как это обычно

имеет место,

можем

принять -

 

9

Л

 

 

 

 

\

 

 

 

 

 

\

 

л * 2" “ -

0,

 

 

вв иду

чего

 

іх— 12 = 0.

 

 

 

 

 

 

 

Как видим, рассматриваемый процесс

непрерывного дрсс-

I селирования протекает без изменения энтальпии

газа.

На основании

(36) можем

написать

 

 

 

Яр (Л - Т2) +

-

Т\) = Арх

 

 

 

\Т \ др (2Ч~Е)-'

 

 

 

 

 

- А р 2

 

 

 

 

К ак и раньше,

из-за малого значения членов правой части

этого

равенства Тх ъ Т 2 должны

быть

близки друг к другу.

4 »

51

Полагая, что

^Т = Т1—Т2 Т —предполагаемое

падение

температуры), представим левую

часть написанного

выраже­

ния

в виде: •

 

 

 

 

ар + —-

(ТI +Тг,) Д7, « ( а р + й Т 1)Д7’ = СР Д7\

где

Ср—теплоемкость идеального

газа при температуре _ Ти

= ар+ ь'т1•

Правую же часть

уравнения,

считая

T2z T u представим как

 

 

 

 

А

'і7 0

'2 +

S )

 

ІРі

Р-і)-

 

 

 

 

 

 

 

RT\1 4-е

 

 

 

 

 

Тогда будем иметь для

всего уравнения

 

 

 

 

 

Ср АГ =

Л

Др Г- -г ;N

 

ІРі - А>).

 

 

 

 

 

 

 

RT\

 

 

 

 

 

откуда

находим

 

 

До (- •("

 

 

 

 

 

 

 

Д Т =

 

- Ь

(Лі — Рі)-

 

 

 

 

 

 

7?rJ+£

 

 

 

 

 

 

 

 

 

 

 

 

-При

рассмотрении этой формулы видим, что при непрерыв­

ном дросселировании газа

могут представиться три

случая:

,

если

д, С + s'

 

,

 

имеет

место при низких

темпера-

з)

___ ■--- >о, что

 

 

 

RT\ "

 

 

 

 

 

 

 

 

 

турах

Ти

то АТ получается

положительным

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

газа будет

падать;

 

 

 

 

 

 

 

 

 

' 2)

если

gf> r2,t --- < ’6,

что имеет

место при

высоких темпе-

 

 

 

RT^*

 

 

 

 

 

 

 

 

 

, ратурах

 

то ДТ

получается, отрицательным

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

^газа будет

повышаться;

 

 

,

 

 

 

 

3) наконец, при - „--у— = Ь температура газа не будет

RT\'

изменяться. Эту температуру, в нашем случае определяемую как

1

а0(2 + в). 1+Е

Т

bR

называют „температурой инверсии .

7. ГАЗОВЫЕ СМЕСИ

Теория смесей идеальных газов базируется на одном из двух законов—Дальтона или Амага. Первый рассматривает давление Я, образуемое газовой смесью объема V и темпера­ туры Т, состоящим, из суммы парциальных давлений Р\ ком­

понентов ^Я = £ Я і)- Второй рассматривает объем смеси V

при давлении Я и температуре Т, состо'ящ ім' из, сумм парци­

дев закона моігут быть

альных объемов Ѵі компонентов ( V — ^ Ѵ і). Оба закона при­

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

R = Ri, где gi и Ri—весовые доли и газовые постоянные

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

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

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

В уравнении Ван-дер-Ваальса

раскроем скобки и, прене­

брегая малым членом — , получим:

V -

 

 

р _ R T __а_

ЬР_

V

V 2

V

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

R.T

 

 

Р через — •, будем иметь:

 

 

V

 

 

p _ R T

а

bRT

V

V -

V -

Уравнение это отнесено к 1 кГ газа. Если в него вместо удельного объема ѵ ввести абсолютный V Gv, м3,* где G— вес газа, кГ, и представить в виде

P = I - G R — — G2 а + — G 4 R ,

I/

Ѵ‘і

Кз

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

Ru ал, Ь\, и вес G\

Р> = ^ öi Ri - - L G-щ +

Gib-Л.

Но так как по Дальтону Р = ^ Р і, то:

£ [ - Г - С ' й - — О'га + Д Of *,/?,]

или

Р = ~ I оА- -у, I au, + X £ а\ь-л.

Сравнивая это выражение с написанным для смеси, видим, что должны иметь место следующие равенства:

(?« = ! № ОЯ = I Gf щ; G2bR = X GibiRi.

53

Введя, в них вместо G и G, весовые доли компонентов

g = — , будем иметь:

G

 

 

 

b = Z ^bi

 

 

R = ^ ,g \R ü

а =

 

 

 

(39)

Эти выражения и должны связать параметры

R, а , b

смеси и Ri, au b\ компонентов.

 

(25)

1

 

 

Далее возьмем уравнение вида

 

 

 

 

RT

а ' , ,

 

 

 

 

V = --------------\-Ь

 

 

 

 

 

. . Р

RT

 

 

 

 

и введем в него вместо удельного

ѵ абсолютный

объем

 

 

V" =

G-V.

 

 

 

 

К, = — GR-----

-

+ G b /

 

 

 

Р

T R

 

 

 

Считая, что уравнение составлено для газовой

смеси» рас­

сматриваемой по Амага

при Р и Т смеси, напишем его для

і-го компонента [смеси

имеющего

параметры Ѵ\, Ru

аи b{ и

вес G{

 

 

 

 

 

 

 

Vi =

j r

GlRl ~

± ^

L

+ Gibu

 

 

и так как по Амага

V = £ VI,

 

 

 

 

 

то'

 

 

 

 

 

 

 

G |Д[ + G\b\

Ri

или

V ■

Сравнивая это выражение с написанным для смеси, видим

[G/? = 5 ]G i/?1; ^ = £ £ £ 1 ; G b = ^ G K

Вводя в эти равенства вместо G и G, весовые доли g, =

%г > будем иметь

и

R=2£іЯі; а= Z

b= Zё'ь'-

(40)

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

между

па­

раметрами смеси и компонентов по Амага.

отношении

Сравнивая формулы (39) и (40),

видим, что в

R они дают одинаковые результаты, выражения же

для а и b

получаются

различными.

 

Покажем, что полученная из закона Дальтона зависимость

Z

2

приводит к ошибочному результату.

0

 

g"i аI

5 4

 

 

 

і

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

а смеси

мы должны были бы написать a =

Но по

•формуле

вместо равенства а = а,-, мы всегда .будем

получать

■а аі, например,

при g \ — gi = 0,5 а = 0,25 а-,.

 

Этот результат,

объясняемый вхождением под знак суммы

квадрата концентраций, показывает, что определение величины ■а с помощью закона Дальтона может привести к значительным юшибкам.

Физическая причина—почему формальное приспособление уравнения Ван-дер-Ваальса к смесям дает погрешности—заклю­ чается в том, что уравнение, будучи написанным для какого-

нибудь индивидуального компонента смеси с членом a\g\

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

Выражения (40), основанные на законе Амага, с математи­ ческой точки зрения не приводят к грубым ошибкам, но сам

закон, основанный на равенстве V = £ V) и именуемый часто

„законом аддитивности объемов“, для .реальных газов экспе­ риментально оправдывается лишь приближенно и потому определение величин а и b по формулам (40) также будет неточным.

Для нахождения величины а для смесей пользуются пра­ вилом Ван Лаара, установленным им па опытным данным и названным „правилом аддитивности квадратных корней из а “.

Согласно этому правилу

Ѵ а =

Ѵ щ .

Что касается выражения для постоянной b смеси, то оно получается при следующих рассуждениях. Так как для инди­ видуального газа Ь—есть объем, занимаемый молекулами 1 кГ газа, то при весовых концентрациях компонентов gi и объемах Ь\ занимаемых молекулами 1 кГ последних, постоянная b смеси

должна равняться b = £ g\b\, что совпадает с результатом,

полученным из закона Амага.

,

1

Итак,

получаем

 

#

=

af = (Z g , KöT)2;

b ^ ^ g i b i .

(41)

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

55

а =

то

для

нахождения

параметров

а0

и е для смесей

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

 

 

 

 

 

 

 

 

'

^ =

( L e

< V ß r

 

.

 

(42)

Если смесь состоит из газов, удовлетворяющих закон соот­

ветственных

состояний,

т. е. е = пост.,

то

для

них

получаем:.

 

 

е = е,-

и а0 =

( £ # V~Q\f ■

 

 

Для смеси, у компонентов которой

параметры

еі получа­

ются

различными,

величины а0 и е могут быть найдены только

при задании интервала температур Ти

. ., 7\, между которыми

протекает

процесс, рассматриваемый

с применением к нему

уравнения

(42).

 

Тогда

 

 

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

равенства эти обращаются

в д е э уравнения с двумя неизвест­

ными а0 и в, которые отсюда и можно определить.

Коэффициент сжимаемости при известных а0 и е может быть,

найден по уравнению

 

Z = 1 — IP — “о

'

Н1 Г2+е

а уравнениями состояния будут служить уравнения (30) и (31)_

8. КАЛОРИЧЕСКИЕ ВЕЛИЧИНЫ СМЕСЕҢ

Поскольку найденные для смесей параметры а0, е, Ь, под ставлремые в уравнения (30) и (31), фактически обращают смесь в воображаемый индивидуальный газ, выведенные для таких газов выражения и, і и s сохраняют силу и для смесей, с той оговоркой, что входящие в них коэффициенты аѵ, ар и b’ также должны быть взяты для смесей. Последние, как известно, выражаются так:

ßv=

Г'Ді ^ѵі;

ap =

£1 g iapc,

'

где avi, api и

bi—коэффициенты теплоемкостей

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

 

g і—их весовые доли.

должны отсчитываться от

Так как калорические

величины

некоторого стандартного

состояния

Р0 и Г,,, то нужно

обеспе­

чить условие,

при котором

для всех компонентов смеси это

состояние было бы одинаковым и реальным: Для этого необ­ ходимо выбирать Р0 н Т 0 такими, чтобы ни для одного из ком­

5 6

понентов смеси температура Т0 не оказалась ниже темпера­ туры насыщения при данном давлении.

■Ввиду того, что коэффициенты аѵ и Ь' не зависят от дав­ ления, следует, задавшись температурой Т0, равной, например, 273° К, подобрать давление Р0 так, чтобы оно было не выше давления насыщения при этой температуре для наиболее высококипящего компонента. .

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

Очевидно, Р0 = 1 • ІО4 кГ/'м*Т0—273 °К возможно только для смесей таких газов, которые имеют температуру кипения при атмосферном давлении ниже 0° С.

ГАЗЫ, СИЛЬНО ОТЛИЧАЮЩИЕСЯ ОТ ИДЕАЛЬНЫХ

Мы видели,

что все уравнения состояния, описывающие

свойства газов

в широком интервале изменяемости давлений

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

сложны и тем они сложнее, чем больше совер­

шенны с точки зрения правильности определения по ним пара­ метров Р, г, Т. Это затрудняет математические операции с такими уравнениями, в частности определение по ним кало­ рических величин. Последние, как правило, получают много­

членные выражения, неудобны?

для технических

расчетов.

В последнее время большую

популярность

приобрели

графоаналитические методы рассмотрения газов, сильно отли­

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

известных диаграмм

сжимаемости

и специальных

графических

построений.1

1. НАХОЖДЕНИЕ КАЛОРИЧЕСКИХ ФУНКЦИЙ і и S

А. Энтальпия і

г

 

 

 

 

Пусть дано состояние газа, определяемое параметрами Р и Т.

Возьмем

состояние

этого

газа

при исходных условиях, а.

за параметры примем давление Р0,

столь низкое, что газ можно

считать идеальным,

и

температуру Т0, которая должна быть

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

насыщения при этом давлении (точка а

на рис. 19). Тогда если г—есть энтальпия в искомом состоянии, изображаемом на рис. 19 точкой С пересечения изобары Р и: изотермы Т, то величина і может быть найдена при переходе

5 7

"из точки а

в точку с

по пути abc. (точка b отвечает состоянию

при давлении Я0 и температуре Т).

Тогда

 

 

 

 

 

 

і

+ h

 

 

 

 

 

(43)

где /р—изобарические

изменения энтальпии

при давлении

Р0;

it—изотермическое изменение энтальпии при температуре Т .

 

 

 

 

Величину ір легко опреде­

 

 

 

 

лить, так как при низком дав­

 

 

 

 

ление,

как для идеального га­

 

 

 

 

за,

будем

иметь

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

*'р = \CfdT

(44)

 

 

 

 

или

 

 

То

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

if = af ( T - T Q) + ^ - ( T * - r l ) .

Рис.' 19. Метод определения

I и 5

Для

нахождения

величины

газа при

высоких давлениях

it

воспользуемся

известным

 

 

 

 

выражением

(22):

 

 

 

' др }т

ZRT-, имеем:

 

 

 

 

 

Заменяя

в нем ѵ

на

 

 

 

 

 

 

=

zART - А Т

L

О

 

Jp

 

 

мли

\ др)Т

 

 

дТ

 

 

 

 

 

 

д.гТ)

 

 

 

 

 

(ІІ Л

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\ д р

 

 

дТ

ü '

 

 

 

Окончательно

 

 

AART- 1/£г

 

 

 

 

 

 

/ д і _ \

 

 

 

 

 

 

 

Vдр

Р

\ дТ

 

 

 

 

откуда изотермическое

изменение

энтальпии:

 

 

' ^~ARrh W dP-

Po p

Заменяя в подынтегральном выражении размерные пара­ метры Р — гс Рк и Т = хТк, а также Р0= гс0Як на безразмерные ■гс, т и гс0, можно представить последнее выражение в, виде

( 45)

-58

или

it = -

ARTf и

(46)

 

где

 

 

 

А

h

dZ

(46')

* . ) d * •

 

О*. jT.

 

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

При наличии такого графика А энтальпия і какого-нибудь газа'в ^ состояния (Р , Т) легко может быть 2 = 1

определена путем снятия с этого

графика при %—— тг0 = — И . X —

Рк Рк

=— значения f\ вычислением по

Тк ,

(44) и (46) величин ір и it и по

(43)—энтальпии і. Коэффициенты а р и b', входящие в выражение для !ір, предполагаются известными.

Построение функций f\ произ- рис. гоОпределение

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

Пусть вертикаль тх=пост. (рис. 20)—есть одна из изотерм рисунка 20, пересекаемая рядом изобар «=пост. Подводя к последним в точках их пересечения с изотермой ряд каса­ тельных и снимая в масштабах г и т отрезки а и Ь, найдем

величины (— ) как отношение — . Умножая последние на

Vдт/г. ь

— , получим ряд частных значений подынтегральной функции

ТС

« = — (— ] для данной изотермы тх при различных іс, после

 

/ я

 

 

 

 

эту

*

 

в координатах,

нетрудно будет построить

изотерму

чего тс \ От

 

 

 

 

 

 

 

 

 

<Р, я (рис. 21).

 

 

 

 

 

 

 

Произведя

далее

по

рис.'21

графическое

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

построенной

кривой

= пост, и определяя значёние

 

 

 

 

ТС

 

ТС

 

 

 

 

 

 

 

о

 

о

 

 

 

как взятую

в

масштабах

и іс площадь

oaxb, можем по.ст-5

роить изотерму тх в координатах

(рис,

22).

 

 

 

 

 

 

)

 

 

 

59

 

 

 

 

 

 

 

 

 

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