Способ управления процессом полимеризации этиленпропиленовых синтетических каучуков

Изобретение относится к способу управления по показателям качества в производстве этиленпропиленовых каучуков на основе моделей расчета показателей качества, которые адаптируются к текущему технологическому режиму. Способ обеспечивает оперативную адаптацию модели под изменяющиеся характеристики процесса и управление технологическим процессом. Способ включает в себя расчет вязкости по Муни, оценку погрешности расчета вязкости по Муни и расчет константы, обеспечивающей адаптацию модели под меняющиеся характеристики процесса, формирование управляющего воздействия. Технический результат - решение задачи оперативного управления по показателям качества с заданной точностью. 1 ил., 1 пр.

 

Изобретение относится к области нефтехимической промышленности, в частности к способам управления производством этиленпропиленовых каучуков на основе моделей для расчета показателей качества продуктов.

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

Известен способ управления качеством продуктов разделения нефтяных смесей методом ректификации на основе вычисления значений показателей качества по технологическим параметрам [Веревкин А.П., Арсланов Ф.А., Иванов В.И., Махов А.Ф., Муниров Ю.М. Изобретение RU 2065761 C1]. Однако способ управления не рассматривает получение и адаптацию моделей.

Известен способ управления процессом каталитического риформинга, в котором задается средневзвешенная температура на входе в реактор на основе расчета октанового числа по модели; при этом модель для расчета октанового числа корректируется по лабораторным данным [Нефть, газ и нефтехимия, Каталитический риформинг, фирма «Эплайд Отомейшн», 1989, №3, с.104]. При параметрической идентификации модели используется система оперативной адаптации на основе собственной запатентованной технологии фирмы-разработчика.

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

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

В работе [Поплавский В.Ф. Автоматизация технологического процесса полимеризации изопрена в производстве синтетического каучука. - Канд.дисс.- М.: - МХТИ, 1985] оценка вязкости осуществляется с помощью уравнения множественной регрессии:

Mh=g₀+g₁γ+g₂N₂+g₃T

где g₀, g₁, g₂ - коэффициенты, γ - концентрация полимера, N₂ - мощность, потребляемая электродвигателем мешалки реактора; T - температура реакции.

В работе [Солдатов Е.А. Моделирование и алгоритмизация адаптивного управления непрерывными процессами растворной сополимеризации. - Канд.дисс - Воронеж: ВГТУ, - 1996] управление процессом полимеризации предлагается осуществлять по характеристической вязкости, так же как и в предыдущем методе, с использованием математической модели.

В работе [Дорофеев В.И. Управление динамическими режимами процесса получения разветвленного полибутадиена оптимального качества. - Канд.дисс. - Воронеж: ВПИ, - 1988] контроль предлагается осуществлять с применением математической модели, где используются эквивалент момента молекулярно-массового распределения второго порядка, модуль, учитывающий кинетику реакции, среднее время пребывания реакционной массы в аппарате, значение параметрической вязкости, конверсия мономера, температура реакции, среднемассовая молекулярная масса, молекулярный вес мономера и моменты молекулярно-массового распределения первого и второго порядка соответственно.

В работе [Дорофеев Д.В., Подвальный С.Л. Применение принципа инвариантности для управления процессом полимеризации бутадиен-стирольного каучука // Системы управления и информационные технологии: Межвуз.сб.науч.тр. - Воронеж, 2001. - С.57-62] предложена математическая модель, состоящая из нескольких модулей, в том числе модуля контроля качества, по которому происходит коррекция общей модели, а в качестве способа управления предложена инвариантная подсистема регулирования процесса полимеризации за счет компенсации основных измеряемых возмущений.

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

Недостатком рассмотренных методов управления по показателям качества (ПК) является то, что адаптация моделей проводится изменением их коэффициентов, а иногда и структуры, что требует наличия большого объема экспериментов и применения процедур аппроксимации.

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

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

Поставленная задача решается тем, что в способе управления процессом полимеризации этиленпропиленовых синтетических каучуков на основе моделей расчета показателей качества, включающем регулирование нагрузки реакторов, соотношения расходов катализатора к сокатализатору, давления в реакторе, уровня заполнения растворителя в реакторе, производится расчет вязкости по Муни Mh, при этом проводят адаптивную подстройку моделей расчета показателей качества продуктов при отклонении наблюдаемого ПК, определяемого по данным лабораторных анализов, от расчетного на величину допустимой погрешности (т.е. при нарушении условия $$\left|M{h}_{лаб}-Mh\right|<{\epsilon }_{доп},$$ где Mh_лаб - значение вязкости, определенное лабораторно, ε_доп - допустимое отклонение); при этом подстройка производится изменением константы модели Mh₀ на величину этого отклонения. Модель вычисления вязкости имеет вид

$$Mh=\Delta Mh+M{h}_0,(1)$$

где ΔMh - суммарное отклонение вязкости от некоторого базового значения, определяемое как

$$\Delta Mh=\Delta M{h}_1+\Delta M{h}_2+\Delta M{h}_3+\Delta M{h}_4,(2)$$

где ΔMh₁, ΔMh₂, ΔMh₃ и ΔMh₄ - отклонения вязкости в зависимости от частных параметров: температуры низа реактора Т_Н, верха реактора Т_В, концентрации этилена С_Э и водорода С_В, соответственно. При этом начальное значение константы модели Mh₀ определяется как величина вязкости, соответствующая некоторому базовому режиму, адаптация модели проводится изменением константы Mh₀=Mh_лаб.-ΔMh.

При управлении по показателям качества продуктов проводится адаптивная подстройка моделей ПК при отклонении наблюдаемого ПК, определяемого по данным лабораторных анализов, от расчетного на величину допустимой погрешности. Подстройка производится изменением константы модели на величину этого отклонения.

Сущность изобретения заключается в том, что для управления по вязкости по Муни используется модель расчета вязкости в режиме реального времени, которая адаптируется путем изменения константы модели Mh₀.

Начальное значение константы модели Mh₀ определяется как величина вязкости, соответствующая некоторому базовому режиму. В базовом режиме: Mh=Mh₀ при $$T_H=T_H^0$$ , $$T_B=T_B^0$$ , $${С}_{Э}={С}_{Э}^0$$ , $${С}_{В}={С}_{В}^0$$ (верхний индекс «0» соответствует базовому режиму).

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

$$\Delta M{h}_1=f\left({Т}_H\right)=a_2\left({Т}_H-{Т}_H^0\right)\left({Т}_H+a_1/a_2-{Т}_H^0\right);\left(3\right)$$

$$\Delta M{h}_2=f\left({Т}_B\right)=b_2\left({Т}_B-{Т}_{{}^B}^0\right)\left({Т}_B+b_1/b_2-{Т}_B^0\right);\left(\text{4}\right)$$

$$\Delta M{h}_3=f\left({С}_{Э}\right)=c_2\left({С}_{Э}-{С}_{Э}^0\right)\left({С}_{Э}+c_1/c_2-{С}_{Э}^0\right);\left(\text{5}\right)$$

$$\Delta M{h}_4=f\left({С}_{В}\right)=d_2\left({С}_{В}-{С}_{В}^0\right)\left({С}_{В}+d_1/d_2-{С}_{В}^0\right),\left(\text{6}\right)$$

где - a₁, a_2, b₁, b_2, c₁, c_2, d₁, d₂ - коэффициенты, определяемые экспериментально. Показано, что эти коэффициенты при изменении технологического режима изменяются незначительно, в то время как Mh₀ может изменяться в значительных пределах.

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

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

Адаптация модели проводится при увеличении погрешности определения по модели Mh вязкости по Муни относительно данных лабораторных анализов (либо измеряемых на потоке) Mh_лаб более чем на заданную величину изменением константы Mh₀

$$M{h}_0=M{h}_{лаб.}-\Delta Mh.\text{(7)}$$

На чертеже приведена принципиальная схема системы управления, реализующая способ управления процессом производства этиленпропиленовых каучуков, включающая в себя реактор полимеризации 1, сепараторы 2 и 3, холодильник-адсорбер 4, испаритель-холодильник 5, насос 6, смеситель 7, блок 8 задачи ограничений, задатчик 9, управляющее устройство 10, блок 11 адаптации моделей, блок 12 измерения, вычислительное устройство 13, а также потоки: растворитель I, этилен II (сырье), пропилен (сырье) III, водород IV, катализатор V, сокатализатор VI, циркулирующий газ от компрессора VII, раствор полимера (продукт) VIII, жидкий аммиак IX, газообразный аммиак X, непрореагировавшие мономеры на факел XI, непрореагировавшие мономеры на узел компримирования XII.

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

Сигналы о значениях расхода катализатора ΔG_кат (FE01), сокатализатора ΔG_сокат (FE02), отдувки ΔF_отд (FE04) и температуры растворителя ΔТ_р (TE03), температуры низа Т_Н (ТЕ05), верха Т_В (ТЕ06), концентрации этилена С_Э (QE07) и водорода С_В (QE08) поступают в вычислительное устройство 13. В вычислительном устройстве 13 рассчитываются по моделям (2)-(5) приращения показателя вязкости по Муни по отклонению измеренных значений относительно базового режима и итоговое значение вязкости по Муни по (1). Рассчитанное значение вязкости по Муни подается на вход управляющего устройства 10. В нем на основе действующих ограничений на параметры режима, задаваемых блоком 8, требуемого качества каучука (Mh_зад), определяемого задатчиком 9, и рассчитанном в вычислительном устройстве 13 значении вязкости по Муни (Mh) по определенному алгоритму производится расчет требуемых расхода катализатора (G_кат), сокатализатора (G_сокат), отдувки (F_отд) и температуры растворителя (T_р) с учетом ограничений. Рассчитанные значения поступают в качестве задания на соответствующие регуляторы сепаратных подсистем, которые воздействуют на клапаны-регуляторы на линиях расхода катализатора G_кат(FC01), сокатализатора G_сокат (FC02), отдувки F_отд (FC04) и на линии откачки газообразного аммиака из холодильника 4 (TC03). Периодически рассчитанные по модели значения вязкости по Муни, сравниваются в блоке 11 адаптации моделей с результатами лабораторных анализов продукта (B) или результатами измерения блоком (B*) 12 и проводится подстройка коэффициентов моделей (C).

Пример реализации предлагаемого способа управления процессом производства этиленпропиленовых каучуков

Расчет значения вязкости по Муни и адаптацию модели ведут в следующей последовательности.

Шаг 1. На стадии формирования модели, конфигурирования системы:

1.1 По каждой измеряемой переменной (фактору) f_i (i=1,2,3,4): f={Т_Н,Т_В,С_Э,С_В} выбирается базовая модель с параметрами $${Т}_{Н}={Т}_{Н}^0$$ ; $${Т}_{В}={Т}_{В}^0$$ ; $${С}_{Э}={С}_{Э}^0$$ ; $${С}_{В}={С}_{В}^0$$ , значение Mh₀ и заданное абсолютное значение допустимых отклонений определения δMh.

1.2 Устанавливается период θ_ПТК определения Mh_лаб, текущее значение таймера θ_Т обнуляется.

Шаг 2. При оперативном управлении в цикле работы контроллера:

2.1 Вводятся измеренные значения факторов f, текущее значение таймера θ_Т.

2.2 По соотношениям (2)-(6) вычисляется суммарное отклонение ΔMh.

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

$$\left|\Delta Mh\right|>\delta Mh$$ ,

то в качестве текущего значения на управляющее устройство 10 отправляется последнее измеренное значение вязкости Mh_лаб и осуществляется переход на п.2.6, иначе выполняется п.2.3.

2.3 По модели (1) вычисляется значение Mh.

2.4 Сравнивается значение θ_ПТК с текущим значением таймера θ_Т.

Если выполняется условие θ_ПТК>θ_Т (это означает отсутствие свежего анализа ПК к текущему моменту времени), то в качестве текущего значения Mh на управляющее устройство 10 отправляется значение Mh и осуществляется переход на п.2.6, иначе выполняется п.2.5 (имеется свежий анализ ПК).

2.5 По соотношению (7) рассчитывается смещение Mh₀, осуществляется сброс таймера и осуществляется переход на п.2.3.

2.6 Конец цикла.

Использование предлагаемого изобретения позволяет решить задачу оперативной адаптации моделей расчета показателей качества получаемых продуктов для целей оперативного управления, что в целом позволит:

1) сократить время на получение моделей для расчета показателей качества получаемых продуктов;

2) повысить точность поддержания показателей качества на заданных уровнях;

3) обеспечить требования минимальной жесткости режима и минимального запаса на качество получаемых продуктов.

Способ управления процессом полимеризации этиленпропиленовых синтетических каучуков на основе моделей расчета показателей качества, включающий регулирование нагрузки реакторов, соотношения расходов катализатора к сокатализатору, давления в реакторе, уровня заполнения растворителя в реакторе, расчет вязкости по Муни Mh, отличающийся тем, что проводят адаптивную подстройку моделей расчета показателей качества продуктов при нарушении условия
$$\left|M{h}_{лаб}-Mh\right|<{\epsilon }_{доп},$$
где Mh_лаб- значение вязкости, определенное лабораторно, ε_доп - допустимое отклонение, изменением константы Mh₀ в модели вычисления вязкости
Mh=ΔMh+Mh₀,
где ΔMh - суммарное отклонение вязкости от некоторого базового значения
ΔMh=ΔMh₁,+ΔMh₂+ΔMh₃+ΔMh₄,
ΔMh₁, ΔMh₂, ΔMh₃ и ΔMh₄ - отклонения вязкости в зависимости от частных параметров: температуры низа реактора Т_Н, верха реактора Т_В, концентрации этилена С_Э и водорода С_В, соответственно, при этом начальное значение константы модели Mh₀ определяется как величина вязкости, соответствующая некоторому базовому режиму, адаптация модели проводится изменением константы Mh₀=Mh_лаб- ΔMh.