УДК 519.8 (075)
Профессор В.Ф. Лебедев, доцент С.В. Подкопаева, доцент Е.А. Хромых
(Воронеж. гос. ун-т. инж. технол.) кафедра информационных и управляющих систем.
тел. (473) 255-387-75
E-mail: helen_hrom@mail.ru
Professor V.F. Lebedev, associate Professor S.V. Podkopaeva, associate Professor E.A. Khromykh
(Voronezh state university of engineering technologies) Department of information and control systems products. phone (473) 255-38-75 E-mail: helen_hrom@mail.ru
Распределение шихты по реакторам каскада в производстве низкомолекулярных сополимеров бутадиена со стиролом
Distribution of the mixture through the reactors cascade in the production of low molecular copolymers of butadiene with styrene
Реферат. При полимеризации бутадиена со стиролом основным фактором, ограничивающим производительность каскада реакторов, является теплосъем. При этом время пребывания мономеров в реакторе значительно превышает время, которое на основании кинетических закономерностей необходимо для завершения процесса. С целью увеличения производительности осуществляется распределение шихты по реакторам каскада. Необходимо распределить поток шихты по реакторам каскада таким образом, чтобы образующийся полимер имел одинаковую вязкость на выходе из каждого реактора. Разработан алгоритм распределения шихты по реакторам каскада с учетом параметров синтеза (температуры в реакторе, приведенной скорости подачи шихты, соотношения модификатора и инициатора в комплексе), количества реакторов и заданной динамической вязкости. В соответствии с разработанным алгоритмом проведен расчет скоростей подачи шихты в каждый реактор каскада. Показано, что расход шихты в каждый полимеризатор зависит от общей производительности установки и от количества реакторов в каскаде. Представленный алгоритм распределения исходной шихты по реакторам каскада позволяет обеспечить максимальную производительность установки и заданное качество получаемого полимера. С целью определения степени превращения мономеров и температурного режима протекания процесса полимеризации при рассчитанной скорости подачи шихты в первый полимеризатор каскада осуществлен расчет основных технологических параметров процесса полимеризации в режиме реального времени с использованием математической модели. Анализ результатов моделирования показывает, что в течение первого часа концентрация мономеров не превышает 0,085 моль/л, что соответствует степени превращения мономера 99 %, при этом температура в реакторе соответствует оптимальному режиму - от 65 до 85 0C.
Summary. In the polymerization of butadiene with styrene heat removal is the main factor limiting the output of the cascade reactor. Thus the residence time of the monomers in the reactor exceeds significantly the time necessary to complete the process on the basis of kinetic regularities. To increase the output of the mixture the distribution in the reactor cascade is made. It is necessary to distribute the flow of mixture through the reactor of the cascade to have the resulting polymer of the same viscosity at the outlet of each reactor. The algorithm of distribution of the mixture in the reactor cascade with regard to the synthesis parameters (temperature in the reactor and the feed rate of the mixture, the ratio of the modifier and the initiator in the complex), the number of reactors and a determined dynamic viscosity is developed. In accordance with the developed algorithm the calculation of the velocity of the mixture feed in each reactor of the cascade is made. It is shown that the flow of mixture in each polymerization unit depends on the overall output of the installation and the number of reactors in the cascade. The algorithm for the distribution of the initial mixture in the reactor cascade is developed to provide maximum output of the installation and set the quality of the obtained polymer. To determine the degree of conversion of monomers and temperature conditions of the process of polymerization under the calculated speed of the feed mixture in the first polymerization cascade the basic technological parameters of the polymerization process in real time mode, the calculation using a mathematical model is made. The analysis of the simulation results shows that during the first hour, the concentration of monomers does not exceed of 0.085 mol/l, which corresponds to the degree of conversion of monomer to 99 %, while the temperature in the reactor corresponds to the optimal mode - from 65 to 85 0C.
Ключевые слова: полимеризация, математическое моделирование, стирол, бутадиен, технологический режим, шихта
Keywords: polymerization, mathematical modeling, styrene, butadiene, technological mode, mixture
© Лебедев В.Ф., Подкопаева С.В., Хромых Е.А, 2014
Процесс полимеризации бутадиена со стиролом протекает в каскаде полимеризаторов. Основным фактором, ограничивающим производительность каскада, является теплосъем, который осуществляется хладоагентом, подаваемым в рубашки охлаждения полимеризаторов. Как правило, при этом время пребывания мономеров в реакторе в 3-4 раза превышает время, которое на основании кинетических закономерностей [3] необходимо для завершения процесса.
С целью увеличения производительности установки осуществляется «дробная» подача шихты во все полимеризаторы каскада. Высокая скорость процесса обеспечивает практически 100%-ную конверсию мономеров на выходе из каждого полимеризатора.
С использованием методов математической статистики найдена зависимость между показателем качества (динамическая вязкость) и параметрами синтеза для заданных интервалов изменения исследуемых технологических параметров [2] вида:
Л
дин
= А - В[ ЮК ] - С ■ Т + Д ■ V
\LiBu ]
пр '
(1)
где А, В, С, Д - коэффициенты; Т - температура в реакторе, 0К; Упр - приведенная скорость подачи шихты; [ЯОК ]/[ПВи ] - соотношение модификатора и инициатора в комплексе.
Эта зависимость позволяет осуществить выбор параметров процесса для получения полимера с заданной динамической вязкостью. Достоверность установленной зависимости подтверждена при выпуске опытных партий в условиях опытного завода [1].
Для получения полимера заданного качества в условиях «дробной» подачи шихты необходимо распределить поток шихты по реакторам каскада таким образом, чтобы образующийся полимер имел одинаковую вязкость на выходе из каждого реактора.
Из зависимости (1) получим:
Л
дин
- А + В
У =■
пР1
[ ЯОК ] [ LiBu ]
+ С ■ Т
Д
(2)
пользованием следующих технологических критериев:
— производительность установки;
— расход дорогостоящего катализатора;
— качество получаемого каучука;
Для решения поставленной задачи выполним вывод расчетных соотношений в общем виде.
Пусть Р - производительность установки полимеризации, кг/час; G1 - скорость подачи шихты в первый полимеризатор, л/мин; G2 -скорость подачи шихты в второй полимеризатор, л/мин; О3 - скорость подачи шихты в третий полимеризатор, л/мин; Gn - скорость подачи шихты в п-ый полимеризатор, л/мин; N - скорость подачи инициатора в первый полимеризатор, моль/мин.
В связи с дополнительной подачей шихты в полимеризаторы каскада, начиная со второго, концентрация инициатора в полимеризаторах непостоянна и определяется из следующих соотношений:
N = Д 1 О
N2 =
N
КЗ + 0г)
(3)
N.
N
Р(О + 02 +...+Оп)
где N12, , п - концентрация инициатора в ^том полимеризаторе, моль/л; в - коэффициент, учитывающий изменение плотности полимеризата.
Приведенной скоростью полимеризации называется отношение скорости подачи мономера, поступающего в полимеризатор (моль/л), к количеству инициатора в полимеризаторе (моль).
По определению приведенной скорости:
УпР1 = ^ О1
УпР2 =
N1 Уг а О
пол 2
N2 Упол
(4)
где Т - температура в ьом полимеризаторе, Упр I - приведенная скорость подачи шихты в >тый полимеризатор.
Необходимо определить режим «дробной» подачи шихты в реакторы каскада с ис-
У =
прп N У '
п пол
где а - коэффициент для перевода л/мин в моль/мин; Упол - объем полимеризатора, л.
По1ставляя в зависимости (4) вместо N1, ..., N их значения из формул (3), получим:
пр1
V.
пр2
00.
NV '
пол
,аР02(01 + 02)
" NVпол '
_ор0И.(01 + 02 +...+0п)
прп NV
пол
Из формул (5) выразим:
NV =
пол
ора12
пр1
= ар 02(° + а2) = =
_ар 0„ .р + 02 +... + а) V
прп
Обозначим:
т1 =02/01, т2=03Ю1; ... ,
тп-1=0у]0и
тогда из (6) получим:
Vпp2 _ 02 .(01 + 02) _
V
пр1
01
+
012 02
V 01 У
= т1 + т1
V
пр3
0з.(0 + 02 + 0з)
V
пр1
012
02
0з 02 0з — + —.— +
01 01 01 101 у
т2 +т1 т2 + т2
V.
0п .(01 + 02 +... + 0п)
V
пр1
02
01 01 01
^!
V 01 у
= т п-1 +т1 тп-1 + ... + тп-1
(5)
(6)
(7)
Из системы уравнений (7) выразим:
V
т1 =.
т2 =
V
пр2 1 1 +---
4 2
пр1
V
прз (1 + т1)
1+т
пр1
(8)
Упрп (1 + т1 + т2 +... + тп_г)
прп
тп, = I —!— +
п-1 11 V
' пр1
1 + т1 + т2 +... + тп_ 2
2
Общая производительность определяется из формулы:
Р=у(01+02+^+0п), (9)
где у - коэффициент для перевода кг/час в л/мин.
Разделим левую и правую часть уравнения (9) последовательно на 01, 02, ..., 0п, после преобразований получим:
01= Р
02 =
0
у(1 + т1 + т2 + .. •. + тп-1)
Рт1
у(1 + т1 + т2 + .. -. + тп-1) (10)
Ртп-1
у(1 + т1 + т2 + ... + тп_ 1)
Рассчитав приведенные скорости подачи шихты в каждый полимеризатор в соответствии с параметрами процесса, количеством полимеризаторов и заданной динамической вязкостью по формуле (7), определяем т1, т2,..., тп по формулам (8) и по формулам (10) рассчитываем
0Ь 02, 0п.
Как следует из уравнений (10), расход шихты в каждый полимеризатор зависит от общей производительности установки и от количества реакторов в каскаде.
Представленный алгоритм распределения исходной шихты по реакторам каскада позволяет обеспечить максимальную производительность установки и заданное качество получаемого полимера.
Приведен расчет скоростей подачи шихты в три полимеризатора каскада согласно уравнениям (1-10).
4
2
2
4
<
Производительность установки равна 800 кг/ч. Исходные данные: расход модификатора 28 кг/ч, расход инициатора 19,2 кг/ч, концентрация модификатора в растворе 0,046 моль/л, концентрация инициатора в растворе 0,17 моль/л, плотность раствора модификатора 0,8 кг/л, плотность раствора инициатора 0,8 кг/л, динамическая вязкость 8 спз, коэффициенты зависимости (1) - А=0,44, В=0,69, С=0,046, Д=0,0088. Результаты расчета приведены в таблице 1.
Т а б л и ц а 1 Распределение шихты по реакторам каскада
№ реактора Температура в реакторе, °С Расход шихты в реактор, кг/ч Расход шихты в реактор, л/мин
1 70 377,556 5,034
2 75 236,842 3,158
3 80 185,602 2,475
С целью определения степени превращения мономеров и температурного режима протекания процесса полимеризации при рассчитанной скорости подачи шихты в 1-ый полимеризатор каскада осуществлен расчет основных технологических параметров процесса полимеризации в режиме реального времени с использованием математической модели (11).
Математическая модель [3] имеет следующий вид:
<М=МП — К ■ М — О, ■ м,
^ У р У
— = К„ ■ М - ■ Р, dг р У
(
dT _ dг
- А1 ■
- В1 ■
dУ,„
12 ■ Т ■ У
1 Э у т
Кр ■ м -
О*
У
Р
У ■ d ■ с
тшш
О ■ d ■ с -(Т - Т )
Ш Ш Ш \ Ш /
(11)
У ■ d ■ с
т
Кр ■ Рст \Т - Тл )
^ ст
У ■ d ■ с
т
dг
КР =
= А ■ О,.
с С ^
кат
V Ут У
12
■ К0 ■
р =
ст
2У_
+ Я
дна'
МП = й„ ■ О,,, -1000 ■
ш ш
ГМДСт+МД6утЛ
М
М
где М - концентрация мономера, моль/л; Р - концентрация полимера в системе, моль/л; Т - температура реакционной массы, К; Ут - текущий объем реактора, л; МП - скорость подачи мономера, моль/мин; Кр - константа скорости роста цепи, 1/мин; К0 - предэкспо-ненциальный множитель, л1/2/(моль1/2мин); Е - энергии активации стадии роста цепей, кДж/моль; R - универсальная газовая постоянная, кДж/(моль^К); Ош - объемная подача мономера, л/мин; dш - плотность шихты (смеси мономеров), кг/л; МДст, МДбут - массовые доли мономеров (стирола и бутадиена) в шихте; Мст, Мбут - молекулярные массы стирола и бутадиена; Скат - текущее количество катализатора (инициатора) на стадии роста полимерных цепей, моль; Тш - температура подаваемой шихты, К; Тхл - средняя температура хладагента, К; сш - удельная теплоемкость шихты, кДжДКкг); Кр - коэффициент теплопередачи, кДж/(м2■мин■К); рст - переменная площадь поверхности теплообмена, м2; г - радиус цилиндрической части реактора, дм; Sдна - площадь поверхности дна реактора, м2; А - коэффициент, учитывающий изменение плотности реакционной массы при образовании полимера; А1 - коэффициент, учитывающий наличие охлаждения реакционной массы вследствие подачи в реактор холодной шихты; В1 - коэффициент, учитывающий наличие охлаждения реакционной массы вследствие подачи хладагента в рубашку; г- время, мин.
Начальные условия:
М(0)=то;
Р(0)=Ро; Т(0)=То;
Ут(0)=Уо.
Результаты расчета для первого реактора каскада представлены на рисунках 1-4.
М,
0.085 •••• 0.05
Рисунок 1. Зависимость концентрации мономера М (моль/л) от времени т (мин) ( - текущий такт квантования по времени)
бут у
<
Е
Г
Рисунок 2. Зависимость температуры реакционной смеси Т (°С) от времени т (мин) ( - текущий такт квантования по времени)
Рисунок 4. Зависимость объема реакционной смеси V (л) от времени т (мин) ( - текущий такт квантования по времени)
Анализ результатов показывает, что в течение первого часа концентрация мономеров не превышает 0,085 моль/л, что соответствует степени превращения мономера 99 %, при этом температура в реакторе соответствует оптимальному режиму - от 65 °С до 85 °С.
Рисунок 3. Зависимость концентрации полимера Р (моль/л) от времени т (мин) ( - текущий такт квантования по времени)
ЛИТЕРАТУРА
REFERENCES
1 Лебедев В.Ф., Подкопаева. С.В., Гайво-ронская И.П. Модернизация АСУ процессом полимеризации в производстве низкомолекулярных каучуков в Воронежском филиале «НИИСК» // Вестник ВГТА. 2009. № 2. С. 36-39.
2 Подкопаева. С.В., Гайворонская И.П. Автоматизированная система управления процессом полимеризации в производстве низкомолекулярных каучуков // Производство и использование эластомеров. 2005. № 3. С. 3-4.
3 Битюков В.К. Математическое моделирование объектов управления в химической промышленности (теория и практика). [и др.]. Воронеж: ВГУИТ, 2011. 196 с.
1 Lebedev V.F., Podkopaeva S.V., Gaivor-onskaia I.P. Modernization of automated control systems for the polymerization process in the production of low molecular weight rubber in the Voronezh branch "NIISK". Vestnik VGTA. [Bulletin of VSTA], 2009, no. 2, pp. 36-39. (In Russ.).
2 Podkopaeva S.V., Gaivoronskaia I.P. Automated control system for the polymerization process in the production of low molecular weight rubbers. Proizvodstvo i ispol'zovanie elasto-merov. [The production and use of elastomers], 2005, no. 3, pp. 3-4. (In Russ.).
3 Bitiukov V.K. Matematicheskoe mod-elirovanie ob"ektov upravleniia v khimicheskoi promyshlennosti [Mathematical modeling of the control objects in the chemical industry]. Voronezh, 2011. - 196 p. (In Russ.).