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

книги / Переработка нефтяных и природных газов

..pdf
Скачиваний:
10
Добавлен:
19.11.2023
Размер:
47.13 Mб
Скачать

Таблица II.1. Коэффициенты для индивидуальных углеводородов в уравнении 11.11

Компонент

Во

Л 0

c„.io-«

ъ

а

с 10-«

а . 10*

7 • 102

Метан

0,04260

1,8550

0,02257

0,00338004

0,4940

0,0025450

0,124359

0,6000

Этан

0,0627724

4,15556

0,179592

0,0111220

0,345160

0,0327670

0,243389

1,180

Пропан

0,0973130

6,87225

0,508256

0,022500

0,94770

0,12900

0,607175

2,200

Изобутан

0,137544

10,23264

0,849943

0,0424352

1,93763

0,286010

1,07408

3,400

к-Бутан

0,124361

10,0847

0,992830

0,0399983

1,88231

0,316400

1,10132

3,400

Изопентан

0,160053

12,7959

1,74632

0,066812

3,75620

0,695000

1,700

4,630

я-Пентан

0,156761

12,1794

2,12122

0,0668120

4,07480

0,82417

1,810

4,750

к-Гексан

0,174813

14,4375

3,31935

0,109131

7,11671

1,51276

2,81086

6,66849

я-Гептан

0,199005

17,5206

4,74574

0,151954

10,36475

2,47000

4,35611

9,00

я-Октан

я-Нонан

2,32095

131,56

3,20417

0,856466

55,1599

0,781821

2,32899

0

н-Декан

6,23189

358,180

0,131900

1,96701

125,122

0,00442954

2,14459

0

СО.,

0,0499101

2,73742

0,00721045

0,136814

0,0149180

0,0847

0,539

Азот

0,0407426

1,053642

0,0023277

0,025102

0,00072841

0,1272

0,53

< 1 »8ркр [3]. При решении уравнения относительно плотности получается несколько значений. Наименьшее из полученных дей­ ствительных значений соответствует плотности паровой фазы, наибольшее — плотности жидкой фазы.

С целью повышения точности расчетов, расширения областей применения уравнения БВР был сделан ряд попыток его усовер­ шенствования: путем корректирования коэффициентов уравне­ ния, представлением их в виде функции от температуры, введе­ нием дополнительных коэффициентов и новых правил смешения при вычислении параметров смеси из параметров чистых компо­ нентов [3—10].

Одной из наиболее точных модификаций уравнения БВР яв­ ляется одиннадцатипараметрическое уравнение состояния, полу­ ченное Старлингом и Ханом [11]. При его разработке одновре­ менно использовали экспериментальные данные по основным тепло­ физическим свойствам (давлению, объему, температуре, энталь­ пии и давлению насыщенных паров) в жидкой и газовой областях с тем, чтобы обеспечить полную согласованность между всеми определяемыми свойствами системы. По уравнению Старлинга— Хана давление Р является функцией температуры Т и мольной

плотности р.

Уравнение имеет вид

 

р = рет + (В'ЯТ -

А, —

+ - |f - - -ff-) р» +

+ (ыгг — а —

р"+ а (а +

р* +

(1 + VP2) exp (—W!) (11.20)

где В0, Л0, CQ, Dо, Е0, b, а, с, d, а , у — параметры уравнения состояния.

Были получены [11] параметры для углеводородов от Q до С3, а также для С02, H2S и азота.

При расчете термодинамических свойств смеси параметры В0,

А 0, С0 и др. определяют по

следующим

правилам смешения:

Во

%iBoi

( 11.21)

i

( 11.22)

(11.23)

(11.24)

(11.25)

(11.26)

(11.27)

(11.28)

32

A .= E

/

(11.29)

*

(П.ЗО)

 

 

£« =

/

(11.31)

i

' ‘ /

где индексы t и / — номера компонентов от I до N, кц — параметры бинарного взаимодействия между t-м и /-м компонентом. Значения этих параметров приве­ дены в табл. II.2. Величина £,-;* приближается к нулю для компонентных пар, образующих раствор, близкий к идеальному, она равна нулю, когда i = /. В рас­ творах, сильно отличающихся от идеального, например, когда один из компо­ нентов легкий углеводород или неуглеводород, величина А// значительно отли­ чается от нуля.

Для расчета термодинамических свойств смеси предложена [11 ] также обобщенная корреляция, в которой уравнение состояния имеет тот же вид, что и уравнение (11.20), применяют правила смешения (II.21"—11.31), но при этом параметры уравнений В0, Ап и т. д. для чистых компонентов представлены как функции от фактора ацентричности со, критической температуры Т кр и кри­ тической плотности рйр.

Ркр^о/ = А± ~Ь 5jC01

(II .32)

9KVlAail(RTкр{) =

Л2 -1- 5а©{*

(11.33)

pKPiC0i/(RTnp() ~

Аз +

 

(11,34)

pfcp.Ve ~ А 4 +

 

(II-35)

P24 >(» l = ^ +

V i

 

<П-36)

p V ' / ( RT«Pi) =

'1« +

B«“ ‘

(П.37)

h p f A * 1**?,) -

Ь +

 

<»-38)

PSp“ i/(« r Kp() =

4> +

 

(” -39)

h S,lD„ J ( R T i t i ) = A s +

B ^ l

(11.40)

 

 

P V K , , ) = A , + V i

 

dl-41)

 

РкрД

/(S lip ,) = -4» + « Л

 

(4 -42)

, Значения

параметров

Л ;- и

В { в уравнениях (11.32)—(11.42)

(j — 1, 2, ...

11) для

нормальных углеводородов от метана до

-октана приведены

ниже:

 

 

 

Индекс

 

 

 

Индекс

 

 

параметра

Ai

 

Bi

параметра

Ai

Bi

(/)

 

(/)

1

0,443690

0,115449

7

0,0705233

— 0,044448

2

1,284380

-0 ,9 2 0 0 7 3

8

0,504087

1,32245

3

0,356306

1,70871

9

0,030745

0,179433

4

0,544979

— 0,270896

10

0,073283

0,463492

5

0,528629

0,349261

11

0,006430

— 0,022143

6

0,484011

0,754130

 

 

 

2 Берлин М. А. и др.

 

 

 

 

зз

Таблица II.2. Значения

параметров бинарного

взаимодействия кц

в уравнениях (11.22) (11.31)

 

 

 

Метан

Этилен

Этан

Пропи­

Пропан

Изо­

«-Бутан

Изо­

н-Пен-

 

лен

бутан

пентан

тан

Метан

0,0

1,0

1,0

2,1

2,3

2,75

3,1

3,6

4,1

Этилен

1,0

0,0

0,0

0,3

0,31

0,41

0,45

0,5

0,6

Этан

1,0

0,0

0,0

0,3

0,31

0,41

0,45

0,5

0,6

Пропилен

2,1

0,3

0,3

0,0

0,0

0,31

0,35

0,4

0,45

Пропан

2,3

0,31

0,31

0,0

0,0

0,31

0,35

0,4

0,45

Изобутан

2,75

0,41

0,41

0,31

0,31

0,0

0,0

0,08

0,1

я-Бутан

3,1

0,45

0,45

0,35

0,35

0,0

0,0

0,08

0,1

Изопентан

3,6

0,5

0,5

0,4

0,4

0,08

0,08

0,0

0,0

я*Пентан

4,1

0,6

0,6

0,45

0,45

0,1

0,1

0,0

0,0

Гексан

5,0

0,7

0,7

0,5

0,5

0,15

0,15

0,0

0,0

Гептан

6,0

0,85

0,85

0,65

0,65

0,18

0,18

0,0

0,0

Октан

7,0

1,0

1,0

0,8

0,8

0,2

0,2

0,0

0,0

Нонан

8,1

1,2

1,2

1,0

1,0

0,25

0,25

0,0

0,0

Декан

9,2

1,3

1,3

1,1

U

0,3

0,3

0,0

0,0

Ундекан

10,1

1,5

1,5

1,3

1,3

1,3

0,3

0,0

0,0

Азот

2,5

7,0

7,0

10,0

10,0

11,0

12,0

13,4

14,8

со 2

5,0

4,8

4,8

4,5

4,5

5,0

5,0

5,0

5,0

H2S

5,0

4,5

4,5

4,0

4,0

3,6

3,4

2,8 |

2,0

 

Гексан

Гептан

Октан

Нонан

Декан

Ундекан

Азот

с о 2

H2S

Метан

5,0

6,0

7,0

8,1

9,2

10,1

2,5

5,0

5,0

Этилен

0,7

0,85

1,0

1,2

1,3

1,5

7,0

4,8

4,5

Этан

0,7

0,85

1,0

1,2

1,3

1,5

7,0

4,8

4,5

Пропилен

0,5

0,65

0,8

1,0

1,1

1,3

10,0

4,5

4,0

Пропан

0,5

0,65

0,8

V

 

 

 

 

 

1,0

1,1

■ 1,3

10,0

4,5

4,0

Изобутан

0,15

0,18

0,2

0,25

0,3

1,3

11,0

5,0

3,6

«-Бутан

0,15

0,18

0,2

0,25

0,3

0,3

12,0

5,0

3,4

Изопентан

0,0

0,0

0,0

0,0

0,0

0,0

13,4

5,0

2,8

н-Пентан

0,0

0,0

0,0

0,0

0,0

0,0

14,8

5,0

2,0

Гексан

0,0

0,0

0,0

0,0

0,0

0,0

17,2

5,0

0,0

Гептан

0,0

0,0

0,0

0,0

0,0

0,0

20,0

5,0

0,0

Октан

0,0

0,0

0,0

0,0

0,0

0,0

22,8

5,0

0,0

Нонан

0,0

0,0

0,0

0,0

0,0

0,0

26,4

5,0

0,0

Декан

0,0

0,0

0,0

0,0

0,0

0,0

29,4

5,0

0,0

Ундекан

0,0

0,0

0,0

0,0

0,0

0,0

32,2

5,0

0,0

Азот

17,2

20,0

22,8

26,4

29,4

32,2

0,0

0,0

0,0

о о 2

5,0

5,0

5,0

5,0

5,0

5,0

0,0

0,0

3,5

H2S

2,0

0,0

0,0

0,0

0,0

0,0

0,0

3,5

0,0

 

Таблица 11.3.

Физические

свойства чистых веществ,

используемых

в обобщенном

уравнении состояния (11.20)

 

 

 

 

 

 

Критическая

 

 

 

Вещество

Критическая

плотность

Молекуляр­

Фактор

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

Ркр

ная

масса,

ацентрично-

 

 

°F

фунт-моль

 

M i

стн О);

 

 

 

куб. фут *

 

 

 

Метан

 

—116,43

0,6274

16,042

0,0130

Этан

 

90,03

0,4218

30,068

0,1018

Пропан

 

206,13

0,3121

44,094

0,1570

Изобутан

 

274,69

0,2373

58,120

0,1830

«-Бутан

 

305,67

0,2448

58,120

0,1970

Изопентан

 

369,00

0,2027

72,146

0,2260

«-Пентан

 

385,42

0,2007

72,147

0,2520

«-Гексан

 

453,45

0,1696

86,172

0,3020

«-Гептан

 

512,85

0,1465

100,198

0,3530

«-Октан

 

563,79

0,1284

144,224

0,4120

«-Нонан

 

610,50

0,1150

128,250

0,4750

«-Декан

 

651,90

0,1037

142,276

0,5400

к-Ундекан

 

692,31

0,0946

156,300

0,6000

Этилен

 

49,82

0,5035

28,050

0,1010

Пропилен

 

197,40

0,3449

42,080

0,1500

Азот

 

—232,60

0,6929

28,016

0,0350

Диоксид углерода

87,80

0,6641

44,010

0,210

Сероводород

 

212,70

0,6571

34,076

0,105

Циклогексан

 

535,60

0,2027

84,156

0,210

Бензол

 

552,00

0,2401

78,108

0,215

Толуол

 

605,50

0,1924

92,134

0,260

* 1 фунт-моль/куб. фут = 16,016 кмоль/м*.

При расчете свойств по уравнениям Старлинга—Хана должны быть использованы значения T KV , ркр. и со,., приведенные

в табл. II.3.

Уравнение (11.20) может быть использовано для определения всех термодинамических свойств. Поскольку оно является моди­ фикацией БВР, для определения плотности можно использовать тот же метод, что и при применении уравнения БВР.

По рекомендации авторов [11] уравнение (11.20) может быть использовано для расчетов удельного объема, энтальпии, энтро­ пии и летучести каждого компонента в широких пределах темпе­ ратур от минус 162 до 260 °С и давлений от 0,01 до 56,0 МПа. При этом получена хорошая сходимость с экспериментальными дан­ ными в расчетах для смесей легких углеводородов в области криогенных температур, где разработанные ранее методы (в том числе БВР) дают большую погрешность.

Обобщенная корреляция была использована для расчета про­ мышленных процессов переработки углеводородов, в том числе процессов выделения азота и гелия из природного газа, криоген­ ной переработки природных газов с использованием в схеме тур­ бодетандеров в области температур от —157 до 143 °С, процессов

36

низкотемпературной сепарации легких углеводородов, низко­ температурной абсорбции в области температур от —73 до 315 °С, а также для расчета холодильных циклов со смешанным хлада­ гентом. Полученные результаты показали, что применение обоб­ щенной корреляции обеспечивает достаточную для инженерных расчетов точность [11].

Наряду с уравнениями БВР и Старлинга—Хана уравнение состояния Редлиха—Квонга (РК) является одним из наиболее надежных для расчета термодинамических функций углеводород­ ных газов. Это довольно простое двухпараметрическое уравне­

ние

имеет вид

[12]

 

 

 

 

 

 

Р =

RTl(v - Ь ) -

у = - - M\v (v - ь)]

(11.43)

где v — мольный

объем,

а,

Ь — константы уравнения.

 

Длячистых

компонентов они имеют значения

 

 

 

 

 

= З Д ^ / Р к р ,

(Ч « )

 

 

 

 

bi = Q^RTup jPup _

(11.45)

где

и Qf, соответственно

равны 0,4278 и 0,0807

 

При применении уравнения к смесям коэффициенты а и b

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

аддитивности

 

 

 

 

 

ЯСМ=

2J 0/1/;

(11.46)

 

 

 

 

^см =

iS ^ibl

(II’47)

где

ai, b[ — коэффициенты

г-го компонента.

 

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

Чу—Праузнитц [13, 14] для увеличения точности расчета термодинамических свойств предложили определять Qa и для каждого компонента, они рекомендовали новое правило сме­ шения для определения постоянной а смеси и ввели константу бинарного взаимодействия klv

| Константы Йа и Йй, одинаковые для всех компонентов в пер­ воначальном уравнении, были заменены константами й а и Й&.,

37

определяемыми для каждого компонента по волюметрическим данным насыщенных паров чистых компонентов в широком ин­ тервале температур — от температуры кипения до критической температуры. Значения безразмерных констант Qa и (для насыщенных паров в уравнениях 11.49, 11.51 и 11.52) приведены ниже:

Компонент

аа

0,0867

Метан . . . .

0,4278

А зо т ................

0,4290

0,0870

Э т а н ................

0,4340

0,0880

Пропан . . .

0,4380

0,0889

Изобутан . . .

0,4420

0,0898

«-Бутан . . .

0,4450

0,0906

Циклогексан

0,4440

0,0903

Бензол . . . .

0,4450

0,0904

Изопентан . .

0,4450

0,0906

Компонент

«а

 

а ь

 

 

 

Диоксид

угле-

0,4470

 

0,0911

рода .

. . .

 

«-Пентан

. . .

0,4510

 

0,0919

«-Гексан

. . .

0,4590

 

0,0935

н-Гептаи

. . .

0,4680

 

0,0952

к-Октан

 

0,4760

 

0,0968

«-Нонан

 

0,4874 *

0,0985 *

н-Декан

 

0,5000

*

0,1000 *

* Константы экстраполированы на основании их зависимости от молекулярной массы соответствующего компонента.

Модифицированная форма уравнения была применена для опи­ сания жидкой фазы. С этой целью константы Qa. и были опре­

делены для каждого компонента в жидкой фазе и применено немного отличное правило смешения параметра а1}-. Значения констант Qa и насыщенных жидкостей в уравнениях 11.49, 11.51, 11.53 приведены ниже:

Компонент

 

Qb

.Метан . . . .

0,4546

0,0872

А зо т ................

0,4540

0,0875

Э т а н ................

0,4347

0,0827

Пропан . . .

0,4138

0,0802

Изобутан . . .

0,4100

0,0790

«-Бутан . .

0,4184

0,0794

Бензол . . . .

0,4100

0,0787

Циклогексан

0,4060

0,0787

Изопентан . .

0,3970

0,0758

Компонент

aа

 

%

 

Диоксид угле­

0,4184

 

0,0794

 

рода . . . .

 

 

«-Пентан . . .

0,3928

 

0,0767

 

«-Гексан . . .

0,3910

 

0,0752

 

н-Гептан . . .

0,3900

*

0,0740

*

к-Октан . . .

0,3900

0,0736

н-Нонан . . .

0,3910

*

0,0738

*

«-Декан . . .

0,3930

0,0755

* Константы экстраполированы на основании их зависимости от молекулярной массы соответствующего компонента.

В модификации Чу—Праузнитца используют следующие пра­ вила смешения:

^см — 2

y ^ i

(П .48)

 

i

 

 

bt =

R T

j P ктр^

(11.49)

° CM = 2

2

у м р и

(11.50)

«/

 

“ « = 0 . / R2J^ PV P«PI

(11.51)

а

 

t) R 2r5p(J/ 2PKP(/ ( » я паровой фазы)

(11.52)

38

aij — |д--(^at- + ®aj) RTl\ PiJ (vKpi -f

j [0,291 —0,04 (%-|-©у)]

(для

жидкой фазы)

 

 

(11.53)

РкР1]- =

zKPijRTKp.j/vKpil

 

(11.54)

(укр3;. +

укр3;.)] 3 (Для

паровой

фазы)

(11.55)

= 0,291 — 0,08 ((0|*-J- (Oj)/2

 

(11.56)

крij — у

TKp.TKPj (I — kij)

 

(11.57)

Константа ki}- характеризуется взаимодействием молекул t-го

компонента с молекулами

/-го компонента,

она

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

не зависит от температуры, плотности и состава и определяется по экспериментальным данным бинарных систем. Значения кон­ стант kVl для некоторых бинарных систем приведены в табл. II.4.

Для повышения надежности и точности уравнения РК в ус­ ловиях низких температур и в присутствии азота в смесях угле­ водородов предложена [15] модифицированная форма уравне­ ния РК, в которой один из постоянных параметров уравнения а заменен на параметр d, зависящий от температуры, с индивиду­ альными коэффициентами для различных веществ.

Уравнение представлено в

следующем

виде:

Р = #7У(о — Ь) — RTbd/lv (v -f d)]

(11.58

где

 

(11.59)

ь = J

bat

i

 

 

bi 0,0867R T Kp ^ j P K^ ^

(11.60)

d — 2

dai

(II.6I)

i

 

 

dt = A t + Bi (TKP(/ T — 1) + Ct (Т КР{/Т -

l)2

 

(11.62)

где Ai, Bi, Сi — коэффициенты индивидуальных веществ.

 

 

 

Зависимость dt от приведенной температуры

была

получена

в результате обработки экспериментальных данных

по

равнове­

сию жидкость—пар и энтальпиям.

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

метра,

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

Ли—Эдми-

стера

[16]:

 

 

Р = RT/(v — Ь) — a![v (v — Ь)] + bc/[v (v — b)(v+ b)]

(II .63)

39

Таблица 11.4. Характеристические константы k{j для бинарных систем в уравнении (11.57)

Метан

Этан

Пропан

Бутан-я илиизо­ бутан

Пентан-« илиизо­ пентан

Гексан-«

Гептан-«

i

Нонан-«

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ё

 

 

 

 

 

 

 

 

9

 

са

X

О)

ч

i

Цикло­ гексан

Бензол

Азот

Диоксид углерода

Метан

 

0

0,01

0,02

0,04

0,06

0,08

0,10

0,12

0,11

0,12

0,08

0,08

0,03

0,05

Этан

 

0,01

0

0

0,01

0,02

0,03

0,04

0,05

0,07

0,08

0,03

0,03

0,05

0,08

Пропан

1

0,02

0

0

0

0,01

0,01

0,02

0,03

0,04

0,05

0,01

0,02

0,09

0,11

«-Бутан

0,04

0,01

0

0

 

0

0

0,01

0,02

 

 

 

0,12

 

Изобутан

|

0

0,03

0

0,01

0,16

«-Пентан

1

0,06

0,02

0,01

0

0

0

0

0

0,01

0,02

0

 

0,12

 

Изопентан

J

0,01

0,18

 

 

«-Гексан

 

0,08

0,03

0,01

0

0

0

0

0

0,01

0,01

0

0,01

0,16

0,20

«-Гептан

 

0,10

0,04

0,02

0

0

0

0

0

0

0

0

0,01

0,18

0,23

н-Октан

 

0,12

0,05

0,03

0,01

0

0

0

0

0

0

0

0,01

0,21

0,28

«-Нонан

 

0,11

0,07

0,04

0,02

0,01

0,01

0

0

0

0

0

0

0,23

0,27

«-Декан

 

0,12

0,08

0,05

0,03

0,02

0,01

0

0

0

0

0

0

0,25

0,29

Циклогексан

 

0,08

0,03

0,01

0

0

0

0

0

0

0

0

0

Бензол

 

0,08

0,03

0,02

0,01

0,01

0,01

0,01

0,01

0

0

0

0

Азот

 

0,03

0,05

0,09

0,12

0,12

0,16

0,18

0,21

0,23

0,25

_

0

Диоксид углерода

 

0,05

0,08

0,11

0,16

0,18

0,20

0,23

0,28

0,27

0,29

0