Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Проектирование, монтаж и эксплуатация теплотехнологического оборудования.pdf
Скачиваний:
12
Добавлен:
29.11.2025
Размер:
12.59 Mб
Скачать

200

Анализ оптимального числа корпусов многоступенчатой выпарной установки целесообразно осуществить с помощью ПЭВМ.

2. Методические указания по расчету схемы установки для разделения бинарной смеси взаимно растворимых компонентов

Задание Рассчитать ректификационную установку непрерывного действия для раз-

деления бинарной смеси при следующих данных:

вид раствора;

производительность установки по исходной смеси G f ;

массовое содержание летучего компонента в исходной смеси a f ;

то же в дистилляте ad ;

то же в кубовом остатке aw ;

давление греющего пара р;

температура готового продукта td;

температура охлаждающей воды на входе в установку t1;

то же на выходе t2.

2.1 Расчет схемы ректификационной установки

Принципиальная схема установки представлена на рисунке 2.1.

Смесь поступает в колонну, предварительно подогретой до температуры кипения на данной тарелке.

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

G f

= Gd + Gw

 

,

(2.1)

 

 

 

G f a f = Gd ad + Gw aw

 

 

где G f , Gd , Gw – расход, соответственно, исходной смеси, готового про-

дукта и кубового остатка; a f ,

ad , aw - содержание летучего компонента в ис-

ходной смеси, готовом продукте, кубовом остатке, %.

201

 

 

 

 

 

II

 

 

 

 

 

охлаждающая

 

 

 

 

 

вода

 

 

 

I

 

III

 

 

VII

IV

охлаждающая

 

 

 

вода

 

 

 

 

Gd, ad дистиллят

конденсат

 

греющий пар

 

V

 

 

 

 

 

 

VI

конденсат

Gw, aw

 

 

 

 

 

 

 

кубовый остаток

 

VIII

 

 

Gf, af

 

 

 

 

исходная

 

 

 

 

 

смесь

 

 

Рисунок 2.1 – Принципиальная схема ректификационной установки

I – ректификационная колонна, II – дефлегматор, III – сепаратор, IV – конденсатор - холодильник, V – сборный бак готового продукта (дистиллята), VI – подогреватель исходной смеси кубовым остатком, VII – подогреватель исходной смеси паром, VIII – сборный бак кубового остатка

Находятся молярные доли летучего компонента в исходной смеси, дистилляте и кубовом остатке по формулам:

 

 

 

 

 

 

 

 

a f

 

 

x f

=

 

 

 

 

 

 

μa

 

 

a f

+

100 a f

 

 

 

 

μa

μb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ad

 

 

xd

=

 

 

 

 

 

 

μa

 

 

ad

+

100 ad

 

 

 

 

 

 

 

 

 

 

 

 

μa

 

 

μb

 

 

 

 

 

 

 

 

 

aw

 

 

xw =

 

 

 

 

 

μa

 

 

aw

+

100 aw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μa

 

 

 

 

μb

 

;

(2.2)

;

(2.3)

(2.4)

где µa и µb – молярная масса летучего и нелетучего компонентов смеси. Определение числа тарелок n в ректификационной колонне осуществляет-

ся графоаналитическим методом с помощью фазовой диаграммы заданной бинарной смеси. Для этого в квадрате Y, X-координат (рисунок 2.2) проводится диагональ. Из литературы (табл. XLIII [7]) выписываются равновесные составы

R =1,3 Rмин + 0,3

202

жидкости и пара для бинарной смеси заданных веществ и их температура насыщения. Перечисленные характеристики оформляются в виде таблицы 2.1.

Таблица 2.1 Равновесный состав смеси

X, моль-%

0

5

10

20

30

40

50

60

70

80

90

100

 

 

 

 

 

 

 

 

 

 

 

 

 

Y, моль-%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t, оС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Строится диаграмма равновесия для заданной бинарной смеси (рисунок 2.2). Откладываются на оси 0Х точки, соответствующие составам кубового остатка хw, исходной смеси хf и дистиллята хd. Проводятся через эти точки вертикали. Находится молярная доля летучего в парах, равновесных жидкости питания y f . Обозначаются точки М и W. Определяется минимальное флегмовое число по формуле

R

=

xd

y f

,

(2.5)

y f

 

мин

 

x f

 

где xd – молярная доля летучего компонента в дистилляте (готовом продукте); xf – то же в исходной жидкости (питании) колонны; y f - то же в паре, равновесном с жидкостью питания.

Рабочее флегмовое число

(2.6)

Проводится рабочая линия концентраций, для чего на оси ординат откладывается отрезок ON, величина которого определяется по формуле

ON =

xd

 

.

(2.7)

R + 1

 

 

 

Соединяются точки N и M, а также точки W и F отрезками (точка F получается при пересечении двух рабочих линий и соответствует составу исходной смеси). Отрезки FM и WF – рабочие линии, соответственно, верхней (укрепляющей) и нижней (исчерпывающей) частей колонны.

 

 

 

 

 

 

 

 

 

203

 

 

 

 

 

 

 

 

 

 

 

Диаграмма равновесия

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

90

 

 

 

 

 

 

 

 

 

 

 

 

 

80

 

 

 

 

 

 

 

 

 

М

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

70

 

 

 

 

 

 

 

 

 

 

 

 

%

60

f

 

 

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

 

 

 

 

 

кривая

моль

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y,

40

 

 

 

F

 

 

 

 

 

 

равновесия

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

W

 

x

 

 

 

 

 

xd

 

 

 

0

x

20

f

40

50

60

70

90

100

 

 

w10

30

80

 

 

 

 

 

 

 

X, моль-%

 

 

 

 

Рисунок 2.2 - Графическое определение числа теоретических тарелок

Проводятся последовательно из точки М горизонтальные до линии равновесия и вертикальные до рабочей линии отрезки между кривой равновесия и рабочими линиями MF и FW. Последнюю горизонталь провести так, чтобы она пересекла вертикальную прямую x = xW. Число полученных при построении ступеней соответствует числу теоретически необходимых тарелок ректификационной колонны nт.

Действительное число тарелок определяется по формуле

nд =

nт

,

(2.8)

 

 

ηт

 

где ηт – КПД тарелки, определяемый по графику (рисунок 3-7) [2] по вязкости разгоняемой жидкости µж и относительной летучести α или по уравнению (2.9), описывающему этот график.

 

 

ηт = 0,49 (μжα)0,245,

(2.9)

μ

 

=μ

 

aср

 

+ μ

 

100 aср

,

(2.10)

ж

a 100

b

100

 

 

 

 

 

 

 

 

g α = 9

 

Тb Тa

,

(2.11)

 

 

 

 

 

 

 

 

 

 

 

Тb +Тa

 

204

где µa и µb – вязкость летучего и нелетучего компонентов при температуре кипения исходной смеси; Тa и Тb – абсолютная температура кипения чис-

тых компонентов; аср – среднее массовое содержание летучего компонента в колонне

а

ср

=

аd +aw

(2.12)

2

 

 

 

Тепловой баланс колонны позволяет определить расход греющего пара на процесс ректификации.

Приходные статьи теплового баланса: с исходной смесью

Q1 = G f c f t f ;

(2.13)

где c f - теплоемкость исходной смеси (жидкости питания),

t f – темпера-

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

a f ;

c

 

= c

 

a f

+

c

100 a f

,

(2.14)

f

a 100

100

 

 

 

b

 

 

здесь ca и cb – соответственно теплоемкость летучего и нелетучего компо-

нента;

с греющим паром

Q2 = D (iп iк),

(2.15)

где iп и iк – энтальпия греющего пара и конденсата, определяются из таб-

лиц [9] по давлению пара; с флегмой

Q3 = Gd R cd td ;

(2.16)

где td – температура дистиллята, находится как температура насыщения по

аd .

205

c = c

 

ad

+ c

100 ad

d

a 100

b

100

Расходные статьи теплового баланса:

с парами, выходящими из колонны

Q4 = Gd (R +1)(rd +cd td );

r

= r

ad

+ r

100 ad

,

 

100

d

a 100

b

 

где ra и rb – соответственно теплота парообразования летучего и нелетучего компонента [7];

с кубовым остатком

Q5 = Gw cw tw ;

(2.17)

(2.18)

(2.19)

(2.20)

где tw – температура кубового остатка, находится как температура насы-

щения по aw ;

c

 

= c

 

aw

+ c

100 aw

,

(2.21)

 

w

 

a 100

b

100

 

 

потери теплоты в окружающую среду принимаются 5 % от общего количества тепла.

Из уравнения теплового баланса определяется расход греющего пара на

процесс ректификации

 

1,05 (Q4

 

+Q5 Q1 Q3).

 

D =

 

(2.22)

 

 

 

 

 

 

 

 

 

 

(iп iк)

 

 

 

Расход воды, поступающей в дефлегматор,

 

 

 

 

 

 

 

 

G

R r

 

 

 

 

 

 

Gд =

с

d(t

2

dt ).

 

 

(2.23)

 

 

 

 

 

 

 

в

1

 

 

 

Расход воды, поступающей в конденсатор-холодильник,

 

 

 

 

 

 

 

 

 

 

 

 

G

к- х

=

Gd (rd +сd (td td ))

(2.24)

 

 

 

 

 

 

 

св

(t2 t1)

 

 

 

 

 

 

 

 

 

 

 

 

 

Общий расход охлаждающей воды в установке

206

 

G = Gд +Gк- х

(2.25)

Диаметр ректификационной колонны

D =

 

4 V

 

,

(2.26)

π ωр

где V – секундный расход паров, движущихся по колонне; ωр – рабочая скорость паров.

Объем паров находится из уравнений объединенного закона газового состояния

 

рV

=

рн Vн

,

 

 

T

T

 

 

 

 

 

 

 

 

 

 

н

 

 

откуда

 

 

 

 

 

 

V =

 

рн Vн T

 

,

(2.27)

 

 

 

 

 

 

рT

 

 

 

 

 

 

н

 

 

где р, V, Т - давление, объем, абсолютная температура пара при заданных условиях; рн, Vн, Тн - то же при нормальных условиях; Vн = 22,4 Z, здесь Z – молярный расход паров по колонне,

 

 

Z =

Gd (R +1)

,

 

(2.28)

 

 

 

 

 

 

 

 

 

 

μсм

 

 

 

 

μ

 

=μ

 

aср

+ μ

 

100 aср

,

(2.29)

см

a 100

b

100

 

 

 

 

 

 

 

 

Т = tср +273 ,

(2.30)

здесь µa и µb – соответственно мольная масса летучего и нелетучего компонентов.

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

t

ср в

=

td +t f

;

t

ср н

=

t f +tw

;

 

 

 

2

 

 

2

 

207

 

 

 

 

 

 

t

ср

=

 

tср в +tср н

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

а

ср в

=

аd + а f

;

а

ср н

=

 

а f + аw

;

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

ср

=

аср в + аср н

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Предельная скорость паров в колонне

 

 

ωпр

определяется по графику (рисунок 3-9 [2])

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

h. На графике ρ

 

и ρ - плотность паров и

 

 

 

 

 

1

 

 

 

 

 

 

3

 

п

ж

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

жидкости. Плотность паров летучего ком-

 

 

 

 

 

 

ρп/ρж

понента

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

п а

=

μа Тн

 

,

 

 

 

 

 

 

 

 

 

 

 

22,4 Т

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

здесь Тd = td + 273 , µа

мольная

масса

летучего

Тн - температура при нормальных условиях. То же для нелетучего компонента

ρп b = μb ТТн ,

22,4 w

(2.31)

(2.32)

1h = 0,6 м

2h = 0,4 м

3h = 0,3 м

4h = 0,2 м

(2.33)

компонента;

(2.34)

здесь Тw = tw + 273 ; µb – мольная масса нелетучего компонента.

 

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

 

 

 

 

 

ρ

ж

=

ρа +ρb

.

(2.35)

 

 

2

 

 

Рабочая скорость паров

 

 

 

 

 

ωр = (0,8...0,9)ωпр

(2.36)

208

2.2 Выбор оптимального варианта ректификационной установки для перегонки бинарной смеси

Рабочее (действительное) флегмовое число R > Rмин (формула 2.6). От значения R зависят капитальные затраты и эксплуатационные расходы на ректификацию (рисунок 2.3). Эксплуатационные расходы (линия 1) прямо пропор-

руб

 

 

 

 

циональны R и определяются расходом теплоноси-

 

 

 

 

 

 

 

 

теля (греющего пара) на испарение жидкости в кубе-

 

 

 

3

 

испарителе. Капитальные затраты (кривая 2) в зави-

 

 

 

 

 

 

 

 

2

 

симости от R имеют минимум,

соответствующий

 

 

 

 

минимальному рабочему объему колонны. Суммар-

 

 

 

 

 

 

 

 

1

 

ные (приведенные) затраты в зависимости от значе-

 

 

 

 

 

 

 

 

Rопт

R

ния флегмового числа (кривая 3) также будут иметь

Rмин

 

 

 

 

 

 

Рисунок 2.3 - Определение опти-

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

мального флегмового числа по ми-

 

 

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

 

минимумом капитальных затрат.

 

 

 

 

 

 

 

 

 

В качестве критерия оптимизации рационально принять минимум приве-

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

 

 

 

 

 

 

П = ЕК + Э, руб/год,

(2.37)

где Е - нормативный коэффициент эффективности капиталовложений, который можно принять равным 0,15 год-1; К - капитальные затраты, руб; Э - эксплуатационные затраты, руб/год.

Для исследования влияния флегмового числа на геометрические размеры ректификационной колонны необходимо по методике описанной выше определить действительное число тарелок при различных значениях флегмового числа R. Для принятых значений флегмового числа определить площадь сечения колонны по формуле

S =

Gd (R +1)

,

(2.38)

 

 

ρ

п

w

 

 

 

 

п

 

 

209

где wп – скорость паров, отнесенная к полному поперечному сечению колонны; ρп – плотность паров.

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

 

wп = С

 

ρж

 

,

 

 

 

 

(2.39)

 

 

ρп

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

С – коэффициент, зависящий от конст-

С

 

 

 

 

 

 

 

 

 

 

 

 

 

0,10

 

 

 

 

 

 

 

рукции

тарелок (колпачковые, ситчатые)

 

и рас-

 

 

 

 

 

В

 

 

0,08

 

 

 

 

 

 

 

стояния между ними h (принять h = 0,2…0,6 м и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,06

 

 

 

 

 

Б

 

оставлять постоянным при изменении R); опреде-

 

 

 

 

 

 

 

0,04

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ляется

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

на

 

рисун-

0,02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

ке 2.3;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100 300

500 700

 

 

 

 

 

 

 

 

 

 

 

 

h, мм

ρж и ρп – плотность жидкости и паров; определить при средней концентрации жидкости питания аf по формуле:

Рисунок 2.3 – Значение коэффициента С в формуле (2.39)

А, Б – колпачковые тарелки, В - ситчатые тарелки

ρж = ρf = ρa

 

a f

+ ρb

100 a f

,

(2.40)

100

100

 

 

 

 

здесь ρa и ρb – плотность жидких летучего и нелетучего компонентов смеси при температуре tf (табл. III [7]); tf – температура кипения смеси при концентрации xf (таблица 2.1);

 

μср 273

 

ρп =

22,4 (273 + t f ),

(2.41)

здесь μср – средняя молярная масса при концентрации в парах летучего y f

μср =μа y f + μb (1 y f ),

(2.42)

где μа и μb – молярная масса летучего и нелетучего компонентов

(табл. XL [7]).

Для полученных сечений колонны объем ее активной части

 

 

 

210

 

 

 

 

 

 

 

V = S H = S h (n 1),

 

 

(2.43)

 

крышка

 

где Н – высота активной части колонны;

h – рас-

 

 

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

 

 

 

 

 

 

 

 

 

 

тарелка

 

 

 

 

 

 

 

 

 

Геометрические

размеры

определяют

металлоем-

h

цилиндриче-

кость колонны, а значит ее стоимость. Капитальные за-

 

 

ская часть

 

 

 

 

 

 

 

 

 

траты складываются из стоимости колонны Цк, стоимо-

 

D

 

сти трубопроводов, арматуры, КИП, фундаментов, затрат

 

 

 

 

 

 

на доставку и монтаж установки, которую можно оце-

 

днище

 

нить в 60…80 % от стоимости колонны, и стоимости

 

 

 

вспомогательного оборудования (испарителя Ци, дефлегматора Цд, подогрева-

теля исходной смеси Цпи, холодильников дистиллята и кубового остатка Цхд и

Цхк, насосов Цн1,

Цн2 )

 

 

 

 

 

 

 

К =1,7 Цк + Цвсп.об

 

 

(2.44)

В

стоимость

 

вспомогательного оборудования,

Цвсп.об ,

достаточно

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

флегмового числа или конструкции колонны другое оборудование остается

практически неизменным. В этом случае

 

 

 

 

К =1,7 Цк + Ци + Цд

(2.45)

Стоимость колонны Цк определяется как произведение массы колонны

Мк на цену за единицу массы цк , т.е.

Цк = цк Мк

(2.46)

Масса колонны равна сумме массы корпуса и всех тарелок

Мк = Мкорп + n Мтар

(2.47)

211

Масса корпуса складывается из масс цилиндрической части, крышки, дни-

ща

 

 

 

Мкорп = Мцил + Мкр + Мдн

 

(2.48)

Масса цилиндрической части колонны определяется ее высотой, зависящей

от числа тарелок n и межтарельчатым расстоянием h , а также диаметром ко-

лонны D

и толщиной стенки обечайки δ

 

 

 

 

 

 

М

 

 

 

+ H

 

δ ρ,

(2.49)

 

цил

= π D (n 1)h + H

сеп

 

 

 

 

 

куб

 

 

где

ρ - плотность

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

для

стали

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

ρ = 7850 кг/м3 ; H

сеп

- высота сепарационного пространства (расстояние от верх-

 

 

 

 

 

 

 

ней тарелки до крышки), принять Hсеп = D , м; Hкуб - высота куба-испарителя

(расстояние от нижней тарелки до днища), принять Hкуб = 2 D , м; δ

- толщина

 

 

 

 

 

 

стенки, принять δ = 0,008...0,012 м; D =

4 S

.

 

 

 

 

 

 

 

π

 

Массу крышки и днища можно приближенно рассчитать по формуле

 

 

Мкр + Мдн 2 D δ ρ

(2.50)

Масса тарелки Мтар зависит от ее конструкции и диаметра и может быть определена из [5].

Стоимость испарителя Ци и дефлегматора Цд зависит от поверхности их

теплообмена F

и F , соответственно, и цены за 1 м2

кожухотрубного тепло-

и

д

 

обменного аппарата цт.а

 

 

Ци = цт.а Fи; Цд = цт.а Fд

(2.51)-(2.52)

212

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

Э = Цп D + Цв G τ, руб/год

(2.53)

где D и G - расход греющего пара в испаритель и охлаждающей воды в дефлегматор (находятся из уравнений теплового баланса соответствующего те-

плообменного аппарата), кг/ч; Цп

и Цв - цена пара и охлаждающей воды,

руб./кг, (постоянно корректируются); τ

- число часов использования установки,

ч/год.

 

 

 

 

 

Тепловая нагрузка на испаритель, Вт

 

 

Q

и

=

D

r ,

(2.54)

3,6

 

 

гр

 

где rгр - теплота конденсации греющего пара.

Поверхность теплообмена испарителя

F

=

 

Qи

 

,

(2.55)

k

 

 

и

 

и

t

и

 

 

 

 

 

 

где kи - коэффициент теплопередачи от конденсирующегося греющего па-

ра к кипящему раствору, Вт/(м2·К), [5]; tи - температурный напор в испарителе, оС,

tи = tгр tw ,

(2.56)

здесь tгр - температура греющего пара, оС, определяется по давлению па-

ра [9].

Тепловая нагрузка на дефлегматор, Вт

Q

д

=

Gd

R r ,

(2.57)

3,6

 

 

d

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]