Отбензинивающая колонна К-1

  • Вид работы:
    Курсовая работа (т)
  • Предмет:
    Другое
  • Язык:
    Русский
    ,
    Формат файла:
    MS Word
    123,88 Кб
  • Опубликовано:
    2015-06-12
Вы можете узнать стоимость помощи в написании студенческой работы.
Помощь в написании работы, которую точно примут!

Отбензинивающая колонна К-1

1. ТЕХНОЛОГИЧЕСКАЯ СХЕМА УСТАНОВКИ

Отбензинивающая колонна К-1 входит в состав установки АТ с двукратным испарением нефти (рис.1). Эта схема технологически гибкая и работоспособная при любом фракционном составе нефти. Благодаря удалению в колонне К-1 лёгких бензиновых фракций в змеевиках печи, в теплообменниках не создается большого давления и основная колонна К-2 не перегружается по парам.

. ИСХОДНЫЕ ДАННЫЕ ДЛЯ РАСЧЕТА

Проведём технологический расчет отбензинивающей колонны мощностью 6 млн т в год по нефти, разгонка (ИТК) которой представлены в табл.1. В качестве дистиллята предусмотрим отбор фракции легкого бензина Н.К.-85оС. Плотность нефти =0,8393.

. ФИЗИКО-ХИМИЧЕСКИЕ СВОЙСТВА СЫРЬЯ

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

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

Разобьём нефть на 9 фракций (компонентов): 28-58оС, 58-72оС, 72-85оС 85-102оС, 102-140оС, 140-180оС, 180-240оС, 240-350оС и 350-К.К. Три первые фракции 28-58оС, 58-72оС и 72-85оС отбираем в качестве дистиллята и шесть остальных - в качестве остатка (полуотбензиненной нефти).

Принципиальная схема установки АТ перегонки нефти

Рис. 1

Таблица 1

Разгонка (ИТК) нефти

№ фракции            Температуры кипения фракций при 1 ат, °С           Выход на нефть, % масс.               Молекуляр-

ный вес (Мi)

фракций


 



отдельных фракций

суммарный



1

28 - 58

2,09

2,09

0,6510

75

2

58 - 72

2,13

4,22

0,6753

-

3

72 - 88

2,45

6,67

0,6925

-

4

88 - 102

2,28

8,95

0,7049

-

5

102 - 115

2,38

11,33

0,7167

108

6

115 - 128

2,38

13,71

0,7285

-

7

128 - 138

2,41

16,12

0,7372

-

8

138 - 150

2,48

18,60

0,7497

-

9

150 - 162

2,58

21,18

0,7657

134

10

162 - 173

2,44

23,62

0,7748

-

11

173 - 184

2,54

26,16

0,7875

-

12

184 - 192

2,13

28,29

0,7973

-

13

192 - 206

2,55

30,84

0,8085

-

14

206 - 217

2,58

33,42

0,8175

166

15

217 - 228

2,65

36,07

0,8250

-

16

228 - 240

2,62

38,69

0,8325

-

17

240 - 252

2,55

41,24

0,8400

-

18

252 - 264

2,65

43,89

0,8468

-

19

264 - 274

2,69

46,58

0,8523

210

20

274 - 289

2,76

49,34

0,8567

-

21

289 - 302

2,69

52,03

0,8641

-

22

302 - 315

2,69

54,72

0,8705

-

23

315 - 328

2,72

57,44

0,8770

-

24

328 - 342

2,79

60,23

0,8832

270

25

342 - 356

2,86

63,09

0,8891

-

26

356 - 370

3,00

66,09

0,8960

-

27

370 - 386

3,10

69,19

0,9032

-

28

386 - 400

3,27

72,46

0,9108

340

29

400 - 418

3,34

75,80

0,9229

-

30

418 - 434

3,27

79,07

0,9267

-

31

434 - 452

3,27

82,34

0,9368

-

32

452 - 500

3,27

85,61

0,9394

410

33

Остаток

14,39

100,00

-

-

Среднюю температуру кипения компонента tср определяем как среднее арифметическое между начальной и конечной температурой кипения фракции.

Молекулярную массу Мi каждого компонента (фракции) можно определить по данным табл.1 или по формуле Воинова:


где Тср - средняя температура кипения фракции, К.

Относительную плотность компонента  определяем через молекулярную массу по формуле Крэга:


или через относительную плотность :

,

где a - средняя температурная поправка относительной плотности на 1К, определяем по эмпирической формуле Кусакова:


Относительную плотность компонента  определяем по данным табл.1 или по уравнению аддитивности:

,

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

 

Таблица 2

Физико-химические свойства сырья

№ компонента

Пределы выкипания фракции

% масс.

tср, C

Мi

% мольн.

1

28-58

2,09

43,0

65,057

5,62427

0,667814

2

58-72

2,13

65,0

72,845

5,12050

0,687682

3

72-85

1,99

78,5

78,103

4,45738

0,702087

4

85-102

2,77

93,5

84,373

5,74712

0,713809

5

102-140

7,77

121,0

97,037

14,03947

0,736061

6

140-180

8,73

160,0

117,590

13,00545

0,777088

7

180-240

13,64

210,0

148,390

16,11139

0,823778

8

240-350

23,55

295,0

212,225

19,44219

0,869474

9

350-К.К.

37,33

-

397,500

16,45224

0,915580

Итого

-

100

-

-

100

-


Пересчет массовых долей в мольные ведём по формуле:


Результаты расчётов физико-химических свойств сырья отбензинивающей колонны приведены в табл.2.

4. МИНИМАЛЬНОЕ ЧИСЛО ТЕОРЕТИЧЕСКИХ ТАРЕЛОК

Проведём расчет методом температурной границы деления смеси.

Для этого определяем мольный отбор дистиллята Е' по отношению к сырью:

,

где D' и F' - мольный расход дистиллята и сырья в колонне, кмоль/ч.

В нашем случае Е' принимаем равным сумме мольных долей первых трёх фракций, которые должны пойти в дистиллят:

Е' = 0,15202

Определяем самую тяжелую фракцию, которая должна пойти в дистиллят - это третья фракция 72-85оС. Задаемся степенью извлечения этой фракции в дистиллят φD3 = 0,85. Это означает, что 85% этой фракции от потенциального её содержания в нефти пойдёт в дистиллят. В общем случае, чем выше степень извлечения фракции, тем больше требуется теоретических тарелок в колонне.

Степень извлечения этой фракции в остаток φW3:

φW3 =1 - φD3 = 0,15

Содержание данной фракции в дистилляте и в остатке рассчитываем по формулам:


 = 0,85·4,45738/0,15202= 0,24923

=0,15·4,45738/(1-0,15202)= 0,00788

Рассчитываем коэффициент распределения ψi этой фракции:

= 31,60878

Принимаем среднее давление в колонне Pср=4,5 ат = 0,45 МПа.

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

В нашем случае ключевыми фракциями являются третья и четвёртая фракции: 72-85оС и 85-102оС. При среднем давлении в колонне Рср находим температуры кипения этих фракций - Т3 и Т4. Для расчётов используем уравнение Ашворта.

Определяем функцию f(То) всех фракций по формуле:

,

где То - средняя температура кипения фракции при атмосферном давлении (табл.2), К

Например, для первой фракции 28-58оС:

 7,4281630

Результаты расчетов для всех фракций приведены в табл. 3.

 
Таблица 3

Значения параметра f(То) фракций

Параметр

Значение параметра

f(T0)1

7,4281630

f(T0)2

6,6299711

f(T0)3

6,2005121

f(T0)4

5,7691102

f(T0)5

5,0837524

f(T0)6

4,2968571

f(T0)7

f(T0)8

2,5855713

f(T0)9

1,5689269


Вычисляем параметр f(Т) для ключевых фракций по формуле:

,

где Рср - среднее давление в колонне, ат

= 4,68922

=4,36297

Температура кипения фракции при данном давлении:

, К

Получаем Т3 = 412,38 К, Т4 = 429,34 К.

Истинная величина ТЕ находится между Т3 и Т4 и определяется методом подбора такого её значения, которое удовлетворяет следующим условиям:

 

428,216 К(TE)= 4,38361

Рассчитываем при температуре TE коэффициенты относительной летучести ai всех фракций:

;

где Рi - давление насыщенных паров фракции определяем по уравнению Ашворта при температуре TE, ат:

;

Например, для первой фракции:

1==2,78758

Результаты расчетов представлены в таблице 4.

 

Таблица 4

Коэффициенты относительной летучести фракций при температуре ТЕ

Обозначение параметра

Значение параметра

12,78758


21,79813


31,35547


40,97816


50,51985


60,19619


70,04881


80,00304


90,00000345



Определяем минимальное число теоретических тарелок в колонне:

= lg(31,60878)/lg(1,35547)=11,35433

. СОСТАВ ДИСТИЛЛЯТА И ОСТАТКА

Находим коэффициенты распределения всех фракций i :


Например, для первой фракции:

=113568,3604

Рассчитываем составы дистиллята и остатка по формулам:


Например, для первой фракции:

При верном подборе ТЕ выполняются условия:


Результаты расчета составов дистиллята и остатка представлены в табл.5.

Таблица 5

Состав дистиллята и остатка

Пределы выкипания фракции

Рi, ат

ai

Yi


1

28-58

12,54411

2,78758

113568,360

0,36995

3,257.10-6

2

58-72

8,09158

1,79813

782,20603

0,33444

0,00043

3

72-85

6,09963

1,35547

31,60878

0,24923

0,00788

4

85-102

4,40172

0,97816

0,77823

0,04629

0,05948

5

102-140

2,33933

0,51985

0,00059

0,00010

0,16555

6

140-180

0,88285

0,19619

0

0

0,15337

7

180-240

0,21963

0,04881

0

0

0,19000

8

240-350

0,01369

0,00304

0

0

0,22928

9

350-к.к.

1,556.10-5

3,45.10-6

0

0

0,19402

S

-

-

-

-

1,00000

1,00000


. МАТЕРИАЛЬНЫЙ БАЛАНС КОЛОННЫ

Таблица 6

Материальный баланс колонны

ФРАКЦИЯ

С Ы Р Ь Ё

 


% мольн.

кмоль/ч

% масс.

кг/ч

тыс т/г

 

28-58

5,62427

236,04052

2,08843

15356,08824

125,30568

 

58-72

5,12050

214,89827

2,12898

15654,26471

127,73880

 

72-85

4,45738

187,06836

1,98705

14610,64706

119,22288

 

85-102

5,74712

241,19664

2,76767

20350,54412

166,06044

 

102-140

14,03947

589,21173

7,77585

57175,33824

466,55076

 

140-180

13,00545

545,81585

8,72882

64182,48529

523,72908

 

180-240

16,11139

676,16674

13,64575

100336,38235

818,74488

 

240-350

19,44219

815,95470

23,55057

173165,98529

1413,03444

 

350-К.К.

16,45224

690,47140

37,32688

274462,38235

2239,61304

 

Итого

100

4196,82420

100

735294,11765

6000

 

ФРАКЦИЯ

Д И С Т И Л Л Я Т


 


% мольн.

кмоль/ч

% масс.

кг/ч

тыс т/г

 

28-58

36,99471

236,02893

33,51556

15355,33404

125,29953

 

58-72

33,44424

213,37665

33,92609

15543,42239

126,83433

 

72-85

24,92262

159,00811

27,10663

12419,05000

101,33945

 

85-102

4,62866

29,53116

5,43842

2491,64030

20,33178

 

102-140

0,00984

0,06277

0,01329

6,09086

0,04970

 

140-180

0

0

0

0

0

 

180-240

0

0

0

0

0

 

240-350

0

0

0

0

0

 

350-К.К.

0

0

0

0

0

 

Итого

100

638,008

100

45815,538

373,855

 

ФРАКЦИЯ

О С Т А Т О К


% мольн.

кмоль/ч

% масс.

кг/ч

тыс т/г

28-58

0,00033

0,01159

0,00011

0,75419

0,00615

58-72

0,04276

1,52162

0,01608

110,84232

0,90447

72-85

0,78847

28,06025

0,31786

2191,59706

17,88343

85-102

5,94764

211,66547

2,59020

17858,90382

145,72866

102-140

16,55463

589,14896

8,29166

57169,24738

466,50106

140-180

15,33700

545,81584

9,30884

64182,48519

523,72908

180-240

18,99976

676,16674

14,55250

100336,38235

818,74488

240-350

22,92769

815,95470

25,11550

173165,98529

1413,03444

350-К.К.

19,40171

690,47140

39,80724

274462,38235

2239,61304

Итого

100

3558,817

100

689478,580

5626,145


Средняя молекулярная масса дистиллята:

=71,81

Относительная плотность дистиллята:

 =0,6860

Средняя молекулярная масса остатка:

=193,74

Относительная плотность остатка:

=0,8520

7. ТЕМПЕРАТУРНЫЙ РЕЖИМ КОЛОННЫ


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

 
Таблица 7

Расчёт температуры верха колонны

Pi

ki

6,67563

1,66891

0,22167

3,99122

0,99780

0,33518

2,86495

0,34797

1,95387

0,48847

0,09476

0,93070

0,23268

0,00042

0,29669

0,07417

0

0,05800

0,01450

0

0,00223

0,00056

0

Итого

-

1,00000


Расчёт ведётся путём подбора такой температуры, при которой уравнение изотермы паровой фазы превращается в тождество:


где ,

Pi - давление насыщенных паров при Тверха, по уравнению Ашворта,

Рверха - давление вверху колонны, примем равным 4 ат.

Искомая температура Тверха=391,42 К = 118,42ОС.

Температура низа Тниза колонны рассчитывается как температура кипения остатка. Расчёт ведётся путём подбора такой температуры, при которой уравнение изотермы жидкой фазы превращается в тождество:


где ,

Pi - давление насыщенных паров при Тниза ,

Рниза - давление вверху колонны, примем равным 5 ат.

Искомая температура Тниза=526,20 К =253,20ОС.

Таблица 8

Расчёт температуры низа колонны

Pi

Кi

39,71752

7,94350

0,00003

29,43316

5,88663

0,00252

24,26328

4,85266

0,03826

19,4137

3,88274

0,23093

12,60267

2,52053

0,41727

6,474371

1,29487

0,19859

2,501605

0,50032

0,09506

0,375243

0,07505

0,01721

0,003646

0,00073

0,00014

Итого

-

1,00000


Температура ввода сырья в отбензинивающую колонну составляет обычно 200-220оС. Примем

 = 493 К = 220ОС

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

Расчет доли отгона производим по методу А.М. Трегубова. Для этого путём последовательного приближения подбираем такое значение мольной доли отгона сырья e`, при котором выполняется тождество:

Таблица 9

Расчёт доли отгона сырья на входе в колонну

Компо- ненты

tcр

Mi

ci,

Pi, кПа

-1+1









28-58

43

65

0,0209

0,3210

0,0542

2816,9465

6,2599

5,2599

0,5917

1,5917

0,0341

0,2133

13,88

58-72

65

73

0,0213

0,2923

0,0494

2008,3132

4,4629

3,4629

0,3896

1,3896

0,0355

0,1585

11,55

72-85

78,5

78

0,0199

0,2544

0,0430

1614,9319

3,5887

2,5887

0,2912

1,2912

0,0333

0,1195

9,33

85-102

93,5

84

0,0277

0,3280

0,0555

1255,7645

2,7906

1,7906

0,2014

1,2014

0,0462

0,1289

10,88

102-140

121

97

0,0778

0,8013

0,1361

771,6880

1,7149

0,7149

0,0804

1,0804

0,1259

0,2160

20,96

140-180

160

118

0,0873

0,7423

0,1268

365,0801

0,8113

-0,1887

-0,0212

0,9788

0,1296

0,1051

12,36

180-240

210

148

0,1365

0,9196

0,1585

126,6667

0,2815

-0,7185

-0,0808

0,9192

0,1725

0,0486

7,20

240-350

295

212

0,2355

1,1097

0,1941

17,6184

0,0392

-0,9608

-0,1081

0,8919

0,2177

0,0085

1,81

350-К.К.

437

398

0,3733

0,9390

0,1823

3,2873

0,0073

-0,9927

-0,1117

0,8883

0,2052

0,0015

0,60

Итого

-

-

1,0000

5,7077

1,0000

-

-

-

-

-

1,0000

1,0000

88,57



Расчёты сведены в табл. 9, где

e` - мольная доля отгона;

сi - массовая доля отдельных фракций в нефти;

ci`, xi`, yi` - мольные доли отгона отдельных фракций в сырье, в жидкой и паровой фазах сырья;

Мi - молекулярный вес отдельных фракций;

Рвх - абсолютное давление в зоне питания, примем его равным среднему давлению в колонне 4,5 ат или 450 кПа;

Pi - давление насыщенных паров отдельных фракций при температуре ввода сырья, по уравнению Ашворта;

Т - температура при которой определяется давление паров, 493 К;

Тср - средняя температура кипения фракции, К.

Искомая величина е`=0,1125.

Молекулярные веса компонентов Mi вычисляем по формуле Воинова. По данным таблицы 9 средний молекулярный вес нефти:


Молекулярный вес паровой фазы My==89

Массовая доля отгона:


. МИНИМАЛЬНОЕ ФЛЕГМОВОЕ ЧИСЛО

Минимальное флегмовое число Rmin определяется по уравнениям Андервуда:


где αi - коэффициент относительной летучести по отношению к ключевому компоненту


где Pi - давление насыщенных паров при температуре ввода сырья;

Pk - давление насыщенных паров ключевого компонента (которым задавались в начале расчета);

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

В общем случае при увеличении  левая часть уравнения возрастает.

q - отношение количества тепла Q, которое надо сообщить сырью, чтобы перевести его в парообразное состояние, к скрытой теплоте испарения сырья Qисп:

 или

где JC - энтальпия сырья при температуре ввода;

JП - энтальпия насыщенных паров сырья;

JЖ - энтальпия кипящей жидкости сырья.

При расчёте минимального флегмового числа возможны следующие варианты.

а) Если сырьё вводится при температуре кипения, то e`=0 и q=1.

б) Если сырьё вводится в виде холодной жидкости, не доведенной до температуры кипения, то q>1.

в) Если сырьё вводится в виде насыщенных паров, то e`=1 и q=0.

г) Если сырьё вводится в виде перегретых паров, то q<0.

д) Если сырьё вводится в виде парожидкостной смеси,

то 0<e`<1 и 1-q=e`.

 
Таблица 10

Расчёт минимального флегмового числа

№ компо- нента

Pi при tFai






1

0,0562

2,8169

1,7443

0,1154

0,3699

0,7591

2

0,0512

2,0083

1,2436

0,1823

0,3344

1,1905

3

0,0446

1,6149

1,0000

0,4214

0,2492

2,3561

4

0,0575

1,2558

0,7776

-0,3832

0,0463

-0,3086

5

0,1404

0,7717

0,4778

-0,1611

0,0001

-0,0001

6

0,1301

0,3651

0,2261

-0,0440

0

0

7

0,1611

0,1267

0,0784

-0,0155

0

0

8

0,1944

0,0176

0,0109

-0,0024

0

0

9

0,1645

0,0032

0,0020

-0,0004

0

0

Сумма

1,0000

-

-

0,1125

1,0000

3,9970


В нашем случае 1-q=e`=0,1125. Методом подбора находим из первого уравнения Андервуда корень , подставляем его во второе уравнение и определяем Rmin. Результаты расчета приведены в таблице 10.

=0,8943

=3,997-1=2,997

отбензинивающий колонна нефть летучесть

9. ОПТИМАЛЬНОЕ ФЛЕГМОВОЕ ЧИСЛО. ОПТИМАЛЬНОЕ ЧИСЛО ТЕОРЕТИЧЕСКИХ ТАРЕЛОК

Приведём два способа расчёта оптимального флегмового числа.

Графический способ Джиллиленда

а) Задаёмся коэффициентом избытка флегмы i=(1,1…1,8).

б) Рассчитываем флегмовые числа:


Например, 3,2967.

в) Находим параметр Хi :


Например, 0,06975

г) Находим параметр Yi:


Например, =0,58551

д) Находим число теоретических тарелок N из уравнения:


Например, =28,80615

e) Находим величину Ni(Ri+1).

Например, N1(R1+1)= 28,80615·(3,2967+1)=123,7701

Расчёты приведены в таблице 10.

Расчёт параметров Rопт и Nопт

iRixiyiNiNi(Ri+1)






1,1

3,2967

0,0698

0,5855

28,8061

123,7701

1,2

3,5964

0,1304

0,5236

24,9345

114,6078

1,3

3,8960

0,1836

0,4747

22,5177

110,2480

1,4

4,1957

0,2307

0,4348

20,8571

108,3680

1,5

4,4954

0,2727

0,4015

19,6433

107,9485

1,6

4,7951

0,3103

0,3734

18,7153

108,4577

1,7

5,0948

0,3442

0,3491

17,9813

109,5928

1,8

5,3945

0,3749

0,3280

17,3851

111,1697


ж) Строим график Ni(Ri+1)=f(Ri):

График зависимости параметра Ni(Ri+1) от флегмового числа

Рис.2

Минимум на полученной кривой соответствует искомым параметрам: Rопт=4,45; Nопт=19,65; опт=1,5.

Аналитический вариант расчёта (по приближённым уравнениям):


2,9967+0,35=4,3959

опт=1,7

опт=1,7·+0,7=20,0024

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

. МЕСТО ВВОДА СЫРЬЯ В КОЛОННУ. РАБОЧЕЕ ЧИСЛО ТАРЕЛОК

Определяем минимальное число теоретических тарелок в концентрационной части колонны ()

,

где α3 и α4 - коэффициенты относительной летучести компонентов при температуре ввода сырья (см. табл. 10).

=7,7028 ~ 8

Оптимальное число теоретических тарелок в верхней части колонны


Отсюда

=13,3308 ~ 14

Рабочее число тарелок в колонне:


где  - к.п.д. тарелки, примем равным 0,6.

=32,75 ~ 33

Рабочее число тарелок в верхней части колонны


=22,21~ 23

В нижней, исчерпывающей части колонны, таким образом, будет 33-23=10 тарелок. На практике для ввода сырья предусматривают до 5 точек вблизи сечения, определённого по этим уравнениям.

. ВНУТРЕННИЕ МАТЕРИАЛЬНЫЕ ПОТОКИ

а) Верхняя часть колонны.

Количество флегмы, стекающей с тарелок верхней части колонны:

=Rопт·D= 4,45·45815,5=203879 кг/ч

Количество паров, поднимающихся с тарелок верхней части колонны:

+D=203879+45815,5=249695 кг/ч

Объём паров:


 м3/с = 281997,72 м3/ч

Плотность паров:

= 8,85519 кг/м3

Относительная плотность жидкости:


где  - температурная поправка по формуле Кусакова.

Относительная плотность при температуре верха колонны:


Абсолютная плотность жидкости  кг/м3

Объёмный расход жидкости:

 м3/ч

б) Нижняя часть колонны.

Количество флегмы, стекающей с тарелок нижней части колонны:

203879 + 735294(1 - 0,0572) = 897114 кг/ч

Количество паров, поднимающихся с тарелок нижней части колонны:

 = 897114 - 689478,6 = 207635 кг/ч

Объём паров:

 м3/с = 9347 м3/ч

Плотность паров:

=22,2140 кг/м3

Плотность жидкости:


где  - температурная поправка по формуле Кусакова.

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


Абсолютная плотность жидкости  кг/м3

Объёмный расход жидкости:

 м3/ч

. ТЕПЛОВОЙ БАЛАНС КОЛОННЫ

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

Уравнение теплового баланса колонны в этом случае будет иметь вид:

 ,

где QF - тепло, поступающее в колонну с сырьём, кВт;

QB - тепло, подводимое в низ колонны, кВт;

QD - тепло, отводимое из колонны с дистиллятом, кВт;

QW - тепло, отводимое из колонны с остатком, кВт;

QХОЛ - тепло, отводимое в конденсаторе-холодильнике, кВт;

QПОТ - потери тепла в окружающую среду, кВт.


где F, Fж, Fп - массовый расход сырья, жидкой и паровой фаз сырья, кг/ч;

iFж - энтальпия жидкой фазы сырья, кДж/кг;

IFп - энтальпия паровой фазы сырья, кДж/кг;

е = 0,0572 - массовая доля отгона сырья (см. раздел 7);

,

где D - массовый расход дистиллята, кг/ч;

iхол - энтальпия холодного дистиллята при температуре его отвода после конденсатора-холодильника, кДж/кг;

,

где W - массовый расход остатка, кг/ч;

iW - энтальпия остатка при температуре при температуре его отвода из колонны 253,2ОС, кДж/кг;

,

где Lор - количество холодного орошения, подаваемого в колонну, кг/ч;

ID - энтальпия паров дистиллята при температуре верха колонны 118,42оС.

,

где Rопт = 4,45 - оптимальное флегмовое число;

iконд - энтальпия жидкого дистиллята при температуре его конденсации, кДж/кг;

Qконд - теплота конденсации паров дистиллята. Для светлых нефтепродуктов эту величину можно рассчитать по уравнению Крэга:

, кДж/кг.

где Тср.м. - средняя молекулярная температура кипения дистиллята, К;

В общем случае средняя молекулярная температура кипения смеси рассчитывается по формуле:

,

где Тi - среднеарифметическая температура кипения узкой фракции в смеси, К:

xi` - мольная доля узкой фракции в смеси.

В нашем случае средняя молекулярная температура кипения дистиллята:

Тср.м. = 316 . 0,3699 + 338 . 0,3344 + 351,5 . 0,2492 + 366,5 . 0,0463 +

+ 394 . 0,0001 = 334,5 К.

= 332,45 кДж/кг.

Примем температуру дистиллята после конденсатора-холодильника и, следовательно, температуру подачи орошения tхол = 35°C.

Энтальпии жидких нефтепродуктов при соответствующих температурах рассчитываются по уравнению Крэга:

, кДж/кг.

Энтальпии паров нефтепродуктов рассчитываются по уравнению Уэйра и Итона:

, кДж/кг.

Например, энтальпия жидкой фазы сырья, поступающего в колонну при 220ОС:

 кДж/кг.

Результаты расчёта энтальпий потоков:

iF = 496,88 кДж/кг


при tF = 220°C

IF = 814,38 кДж/кг


при tF = 220°C

ID = 593,81 кДж/кг


при tD = 118,42°C

iхол = 74,51 кДж/кг


при tхол = 35 °C

iW = 582,25 кДж/кг


при tW = 253,2°C


Количество холодного орошения:

 кг/ч

Рассчитываем тепловые потоки:

QF = 735294,118 . 0,0572 . 814,38 + 735294,118 . (1-0,0572) . 496,88 =

кДж/ч = 105196,27 кВт

QD = 45815,538 . 74,51 = 3413744,7 кДж/ч = 948,26 кВт

QW = 689478,58 . 582,25 = 401448506 кДж/ч = 111513,58 кВт

QХОЛ = (45815,538 + 130521,12) . (593,81 - 74,51) = 91571622 кДж/ч =

,56 кВт

Примем потери тепла в колоне 5%:

Qпот = (948,26 + 111513,58 + 25436,56).5/95 = 7257,81кВт

Тепло, необходимое подвести в низ колонны:

QB = 145156,21 - 9514,41 - 95681,86 = 39959,94 кВт

Таблица 11

Тепловой баланс колоны

Поток

t, °С

Энтальпия, кДж/кг

Расход, кг/ч

Количество тепла, кВт


ПРИХОД:

С сырьём:





паровая фаза

220,0

814,38

42058,8

9514,41

жидкая фаза

220,0

496,88

693235,3

95681,86

В низ колонны




39959,94

Итого




145156,21


РАСХОД:

С дистиллятом

35,0

74,51

45815,54

948,26

С остатком

253,2

582,25

689478,58

111513,58

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




25436,56

Потери




7257,81

Итого




145156,21


13. ДИАМЕТР КОЛОННЫ

Диаметр колонны рассчитывается по наиболее нагруженному сечению по парам. В нашем случае в верхней части колонны расход паровой фазы больше в 7,8327/2,5964 = 3,02 раза, чем в нижней (см. раздел 11).

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

Таблица 12

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

Диаметр колонны, м

Расстояние между тарелками, мм

до 1,0

200-300

1,0-1,6

300-450

1,8-2,0

450-500

2,2-2,6

500-600

2,8-5,0

600

5,5-6,4

800

более 6,4

800-900


Расстояние между тарелками принимается в зависимости от диаметра колонны (см. табл.12). На практике указанные рекомендации не всегда выполняются. Для большинства колонн расстояния между тарелками принимаются таким образом, чтобы облегчить чистку, ремонт и инспекцию тарелок: в колоннах диаметром до 2 м - не менее 450 мм, в колоннах большего диаметра - не менее 600 мм, в местах установки люков - не менее 600 мм. Кроме этого, в колоннах с большим числом тарелок для снижения высоты колонны, её металлоёмкости и стоимости расстояние между тарелками уменьшают.

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

Диаметр рассчитывается из уравнения расхода:

, м

где VП - объёмный расход паров, м3/с;

Wmax - максимальная допустимая скорость паров, м/с

, м/с

где Сmax - коэффициент, зависящей от типа тарелки, расстояния между тарелками, нагрузки по жидкости;

rж и rп - плотность жидкой и паровой фазы, кг/м3.

Сmax = K1 . K2 . C1 - К3(l - 35)

Значение коэффициента С1 определяем по графику в зависимости от принятого расстояния между тарелками (см. приложение). С1 = 1050.

Коэффициент К3 = 5,0 для струйных тарелок, для остальных тарелок К3 = 4,0.

Коэффициент l находится по уравнению:

,

где LЖ - массовый расход жидкой фазы в верхней части, кг/ч;

Коэффициент К1 принимается в зависимости от конструкции тарелок:

Колпачковая тарелка.................................................................... 1,0

Тарелка из S-образных элементов............................................... 1,0

Клапанная тарелка....................................................................... 1,15

Ситчатая и струйная тарелка....................................................... 1,2

Струйная тарелка с отбойниками................................................ 1,4

Коэффициент К2 зависит от типа колонны:

Атмосферные колонны................................................................. 1,0

Ваккумные колонны с промывным сепаратором в зоне питания        1,0

Вакуумные колонны без промывного сепаратора..................... 0,9

Вакуумные колонны для перегонки

пенящихся и высоковязких жидкостей........................................ 0,6

Абсорберы.................................................................................... 1,0

Десорберы..................................................................................... 1,13


Сmax = 1,15 . 1,0 . 1050 - 4(132,75 - 35) = 816,5

 = 0,562 м/с

Диаметр колонны:

 м

Полученный диаметр округляется в большую сторону до ближайшего стандартного значения. Для стальных колонн рекомендованы значения диаметров от 0,4 до 1,0 м через каждые 0,1 м, от 1,2 до 4,0 м через 0,2 м, далее 2,5 м, 4,5 м, 5,0 м, 5,6 м, 6,3 м, от 7,0 до 10 м через 0,5 м, от 11,0 до 14,0 м через 1,0 м, от 16,0 до 20,0 м через 2,0 м.

Итак, примем диаметр колонны DK = 4,5 м.

Проверяем скорость паров при принятом диаметре колонны:

 м/с

Она находится в допустимых пределах (0,4-0,7 м/с) [5] для колонн под давлением и расстоянии между тарелками 600 мм.

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

м3/(м . ч),

где LV - объёмный расход жидкости, м3/ч;

n - число потоков на тарелке;

W - относительная длина слива, обычно находится в пределах 0,65-0,75.

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

. ВЫСОТА КОЛОННЫ

Высота колонны рассчитывается по уравнению:

НК = H1 + Hк + Ни + Нп + Н2 + Нн + Но, м

где Н1 - высота от верхнего днища до верхней тарелки, м;

Нк - высота концентрационной тарельчатой части колонны, м;

Ни - высота исчерпывающей, отгонной тарельчатой части колонны, м;

Нп - высота секции питания, м;

Н2 - высота от уровня жидкости в кубе колонны до нижней тарелки,м;

Нн - высота низа колонны, от уровня жидкости до нижнего днища, м;

Но - высота опоры, м.

Высота Н1 (сепарационное пространство) принимается равной половине диаметра колонны, если днище полукруглое, и четверти диаметра, если днище эллиптическое. Полушаровые днища применяют для колонн диаметром более 4 метров. Поэтому Н1 = 0,5 . 4,5 = 2,25 м.

Высоты Hк и Ни зависят от числа тарелок в соответствующих частях колонны и расстояния между ними:

Нк = (Nконц - 1)h = (23 - 1)0,6 = 13,2 м

Ни = (Nотг - 1)h = (10 - 1)0,6 = 5,4 м

где h = 0,6 м - расстояние между тарелками.

Высота секции питания Нп берётся из расчёта расстояния между тремя-четырьмя тарелками:

Нп = (4 - 1)h = (4 - 1)0,6 =1,8 м

Высота Н2 принимается равной от 1 до 2 м, чтобы разместить глухую тарелку и иметь равномерное распределение по сечению колонны паров, поступающих из печи. Примем Н2 = 1,5 м.

Высота низа (куба) колонны Нн рассчитывается, исходя из 5-10 минутного запаса остатка, необходимого для нормальной работы насоса в случае прекращения подачи сырья в колонну:

 м

где rж - абсолютная плотность остатка при температуре низа колонны (см. раздел 11);

Fк =  - площадь поперечного сечения колонны, м2.

Штуцер отбора нижнего продукта должен находится на отметке не ниже 4-5 м от земли, для того, чтобы обеспечить нормальную работу горячего насоса. Поэтому высота опоры Но конструируется с учётом обеспечения необходимого подпора жидкости и принимается высотой не менее 4-5 м. Примем Но = 4 м.

Полная высота колонны:

НК = 2,25+13,2+5,4+1,8+1,5+5,25+4 = 33,4 м

. ДИАМЕТРЫ ШТУЦЕРОВ

Диаметры штуцеров определяют из уравнения расхода по допустимой скорости потока:

 , м

где V - объёмный расход потока через штуцер, м3/с;

Величина допустимой скорости Wдоп принимается в зависимости от назначения штуцера и фазового состояния потока (м/с):

Скорость жидкостного потока:

на приёме насоса и в самотечных трубопроводах.................. 0,2-0,6

на выкиде насоса....................................................................... 1-2

Скорость парового потока:

в шлемовых трубах и из кипятильника в колонну.................. 10-30

в трубопроводах из отпарных секций..................................... 10-40

в шлемовых трубах вакуумных колонн.................................. 20-60

при подаче сырья в колонну.................................................... 30-50

Скорость парожидкостного потока при подаче сырья в колонну

(условно даётся по однофазному жидкостному потоку)........ 0,5-1,0

Рассчитанный диаметр штуцера далее округляется в большую сторону до ближайшего стандартного значения:

Таблица 13

Стандартные значения диаметров штуцеров

Dу, мм

Dу, мм

Dу, мм

Dу, мм

Dу, мм

Dу, мм

10

50

200

600

1400

2600

15

65

250

800

2800

20

80

300

900

1800

3000

25

100

350

1000

2000


32

125

400

1100

2200


40

150

500

1200

2400



Приложение

График зависимости коэффициента С1 от расстояния между тарелками Нт

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

- кривая нормальных нагрузок для колпачковых тарелок;

- кривая для вакуумных колонн без ввода водяного пара и для стриппинг-секций атмосферных колонн;

- кривая для вакуумных колонн с вводом водяного пара и для десорберов;

- кривая для абсорберов;

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

Список рекомендуемой литературы

1. Ахметов С.А. Технология глубокой переработки нефти и газа: Учебное пособие для вузов. - Уфа: Изд-во «Гилем», 2002. - 672 с.

. Мановян А.К. Технология первичной переработки нефти и природного газа: Учебное пособие для вузов. - М.: Химия, 2001. - 568 с.

. Танатаров М.А., Ахметшина М.Н., Фасхутдинов Р.А. Технологические расчеты установок переработки нефти. - М.: Химия, 1987. - 352 с.

. Сарданашвили А.Г., Львова А.И. Примеры и задачи по технологии переработки нефти и газа. - М.: Химия, 1973. - 272 с.

. Эмирджанов Р.Т., Лемберанский Р.А. Основы технологических расчетов в нефтепереработке и нефтехимии. - М.: Химия, 1989. - 192 с.

. Александров И.А. Ректификационные и абсорбционные аппараты. - М.: Химия, 1979. - 280 с.

. Александров И.А. Перегонка и ректификация в нефтепереработке. - М.: Химия, 1981. 352 с.

. Багатуров С.А. Основы теории и расчета перегонки и ректификации. - М.: Химия, 1974. - 440 с.

. Кузнецов А.А., Кагерманов С.М., Судаков Е.Н. Расчеты процессов и аппаратов нефтеперерабатывающей промышленности. - Л.: Химия, 1974. - 344 с.

. Расчеты основных процессов и аппаратов нефтепереработки: Справочник / Под ред. Е.Н.Судакова. - М.: Химия, 1979. - 569 с.

. Основные процессы и аппараты химической технологии: Пособие по проектированию / Под ред. Ю.И.Дытнерского. М.: Химия, 1983. - 272 с.

. Справочник нефтепереработчика: Справочник / Под ред. Г.А. Ластовкина, Е.Д.Радченко и М.Г.Рудина. - Л.: Химия, 1986. - 648 с.

. Рудин М.Г. Карманный справочник нефтепереработчика. - Л.: Химия, 1989. - 464 с.

. Скобло А.И., Молоканов Ю.К., Владимиров А.И., Щелкунов В.А. Процессы и аппараты нефтегазопереработки и нефтехимии. - М.: ООО "Недра-Бизнесцентр", 2000. - 677 с.

. Колонные аппараты. Каталог ВНИИнефтемаш. - М.: Изд. ЦИНТИхимнефтемаш, 1992. - 26 с.

Похожие работы на - Отбензинивающая колонна К-1

 

Не нашли материал для своей работы?
Поможем написать уникальную работу
Без плагиата!