Научная статья на тему 'Построение кусочно-линейной модели процесса окислительного пиролиза'

Построение кусочно-линейной модели процесса окислительного пиролиза Текст научной статьи по специальности «Химические технологии»

CC BY
105
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АЦЕТИЛЕН / ОКИСЛИТЕЛЬНЫЙ ПИРОЛИЗ / КУСОЧНО-ЛИНЕЙНАЯ МОДЕЛЬ / ЦЕНТР БЕЗОПАСНОСТИ / MATHEMATICAL MODELING / ACETYLENE / REACTOR OF OXIDIZING PYROLYSIS / PIECEWISE AND LINEAR MODEL / AREA OF SAFETY / CENTER OF SAFETY

Аннотация научной статьи по химическим технологиям, автор научной работы — Санаева Наталья Александровна, Санаева Галина Николаевна, Пророков Анатолий Евгеньевич, Богатиков Валерий Николаевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по химическим технологиям , автор научной работы — Санаева Наталья Александровна, Санаева Галина Николаевна, Пророков Анатолий Евгеньевич, Богатиков Валерий Николаевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THE CONSTRUCTION OF THE PIECEWISE-LINEAR MODEL OF THE PROCESS OF OXIDATIVE PYROLYSIS

The article presents a piecewise linear model of the process of obtaining acetylene by oxidative pyrolysis of natural gas. The system of linear constraints, defined area and the center of technological safety solution of a nonlinear programming problem.

Текст научной работы на тему «Построение кусочно-линейной модели процесса окислительного пиролиза»

УДК 681.5

Санаева Н.А., Санаева Г.Н., Пророков А.Е., Богатиков В.Н.

ПОСТРОЕНИЕ КУСОЧНО-ЛИНЕЙНОЙ МОДЕЛИ ПРОЦЕССА ОКИСЛИТЕЛЬНОГО ПИРОЛИЗА

Санаева Наталья Александровна, студент 3 курса бакалавриата факультета кибернетики;

Санаева Галина Николаевна, старший преподаватель кафедры вычислительной техники и информационных

технологий, e-mail: gn_san@mail.ru;

Пророков Анатолий Евгеньевич, к.т.н., доцент, заведующий кафедрой вычислительной техники и информационных технологий;

Новомосковский институт РХТУ им. Д.И. Менделеева, Новомосковск, Россия 301665, Тульская область, г. Новомосковск, ул. Дружбы, д. 8

Богатиков Валерий Николаевич, д.т.н., профессор, профессор кафедры информационных систем, Тверской государственный технический университет, Тверь, Россия

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

Ключевые слова: ацетилен, окислительный пиролиз, кусочно-линейная модель, центр безопасности.

THE CONSTRUCTION OF THE PIECEWISE-LINEAR MODEL OF THE PROCESS OF OXIDATIVE PYROLYSIS

Sanaeva N.A., Sanaeva G.N., Prorokov A.E., Bogatikov V.N.*

D. Mendeleev University of Chemical Technology of Russia, Novomoskovsk, Russia

* Tver State Technical University , Tver, Russia

The article presents a piecewise linear model of the process of obtaining acetylene by oxidative pyrolysis of natural gas. The system of linear constraints, defined area and the center of technological safety solution of a nonlinear programming problem.

Keywords: mathematical modeling, acetylene, reactor of oxidizing pyrolysis, piecewise and linear model, area of safety, center of safety.

Процесс производства ацетилена окислительным пиролизом природного газа протекает в условиях неполноты и неопределенности информации, вызванными изменением характеристик

оборудования в процессе его эксплуатации («закоксовыванием» ряда элементов),

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

окислительного пиролиза.

Технологический процесс, протекающий в реакторе окислительного пиролиза, состоит из трех этапов [1]:

1) перемешивание предварительно подогретых до 650°С исходных компонентов для реакции пиролиза в смесителе реактора; 2) окислительный пиролиз в зоне реакции при температуре около 1500°С; 3) «закалка» продуктов реакции водой с целью охлаждения полученного газа пиролиза до температуры ниже 100°С, поскольку при температуре реакции ацетилен представляет собой

неустойчивое соединение и может разлагаться на водород и сажу.

Количество ацетилена, получаемого

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

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

Для решения задачи повышения уровня технологической безопасности производства ацетилена целесообразно создание специальной системы оценки состояний и диагностики нарушений с использованием кусочно-линейной модели указанного процесса. При этом функционирование процесса окислительного пиролиза рассматривается как последовательная схема смены его состояний на некотором интервале времени 4). Состояние системы в каждый момент времени t из этого интервала характеризуется набором параметров процесса. Область номинальных режимов (область работоспособного состояния) процесса определяется совокупностями

технологических, конструктивных параметров и ограничения, «вырезают» на множестве всех

параметров управления.

г 02_

СН4

на охлаждение и л ГРР 2

на "закалку" НгО газа пир ошпа

Ч>

газораспределительная решетка (ТТР)

НгО_

состояний процесса и-мерную область, в которой процесс не выходит во внештатные ситуации -область всех работоспособных состояний процесса [2].

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

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

зона закалки

вода с сажей

Рис.1. Реактор окислительного пиролиза

ограничения его рабочего функционирования, выход за которые означает переход процесса во внештатную ситуацию. Таким образом, эти

- Смеситель:

= ОСН4РСН4 + Оо2Ро2 " ОсмРс

а?

ОСН4Рсн4 + Оо2Ро2 - ОсмРсм >АРсм(кX ОСН4РсН4 + Оо2Ро2 _ ОсмРсм <АРсм(к);

- Реакционная зона (для реакции 2СН4 = С2Н2 + 3Н2):

УрР,

тСн ас,

сн4 и^сн.

УрР

рг см

т

тсм ^

= ОсмРс

т

т

СН4 (с о — с

С СН с с

■'СН4 СН.

р см

т

а?

= ОсмРс

т

т

С2Н2 (с0 _ С

С С Н с<

С 2Н 2 С 2Н 2

тН &Сн

V р н2__н2 = О р

ргсм 1 , _ см^с

см

т

р см

тсм а?

т

Н2 (сН — с

) АСН4 (ссн4 , сС2н2, сн2,Т)' ) + ^С2Н2 (сСН4 ' ^^С2Н2 ' ' Т)? )+ АН2 (сСН4 , сС2Н2 , сн2,Т),

ср-

атв

а?

^СН4^СН4+ср(свх _ твых)) —

Ос

Ос

Ос

Рсн

тс

4-сСН4 _ к5Гр (ССн4)2

(

Осм -Рсн^ _ 2к5ГсСн

см 5 Р СНл

тсм 4

сСНП4 >дссн4(к),

рсн4

тс

Рсн,

сСН4 _ к5Ур (Сн4 )2 _ ГОСм ^ - 2к5Грссн4

V тсм

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

сСНХ <Дссн4(к),

тс

4-сСхн2 + къУ, (ССН4 )2

Ос

Рсн,

/ Л

тс

4 + к5УрсС2Н2

ссН4 >Ас_(к),

Осм ^ ссхн2 + к5Ур (сСн4 )2 _ I Осм ^+к, УрсС

т

т

р^С2Н2

сСНХ <Дсс2Н2(к),

Осм ^ снх + к5Ур (сСн4 )2

Ос

т

Рсн

тс

4 снх2 + к5Ур (ССН4)

Осм ^ + к5УрсН2

Ос

т

Рсн

сСННП4 >дс^(к),

тс

^ + к5УрсН2

см

вых

сснх <дсн,(к),

(Осн, Р сн, + Оо2 Ро2 )т'4 + Осн 4Гсн4 _ ОГпсТГПах >АТвых >(Осн4Рсн4 + Оо2Ро2 )сТтах + Осн4гСн4 _ ОГпсТГП'4;

22

- Зона «закалки»:

срУз

1 Т^вых аТГП

а?

: ОГПсР(тГп _ тгпх)_ КГ( _ ТвТ)_ ОГП^

СВЗрВЗУз '

[ВЗ

а?

вх вых ТВЗ _ ТВЗ

)+ КГ( _ ТввЗых),

см

о

GrncpTra - Grnap - KFTB3 - (Grncp + KF) > cpAT^,

Grn cPTi

вых ГП

- KF (T

вых ГП

B3 min

T'B

n

) < AcpATIB

GBm3ncB3PB3 (TBB3X - TBB3IX ) + KF(Trn - TbT ) > cB3PB3F3 ATB

B3

GB3 cB3pB3(TB3 TB3 ) + KF(Trn TB3 ) < CB3pB3V3 ATBB

где р - плотность, кг/м ; G - расход, м3/с;

FCM, Vp, V3 - объемы смесителя, зоны реакции, зоны «закалки», соответственно, м3; mi - молекулярный вес компонента i, кг/моль; С - концентрация компонента, % (об.); с - теплоемкость газа пиролиза, Дж/(кг-К); Т - температура, К;

ki - константа скорости химической реакции, с-1; Ai - скорость изменения количества соответствующих компонентов реакции моль/с; r - удельная теплота сгорания метана, Дж/м3;

г^ 2

F - поверхность теплопередачи, м ; a - коэффициент теплоотдачи, Вт/(м3-К); индексы: см - метано-кислородная смесь; ГП - газ пиролиза; ВЗ - вода на «закалку».

Полученные соотношения позволяют прогнозировать возникновение постепенных отказов, являющихся результатом ухудшения свойств веществ, участвующих в ХТП (например, нестационарность состава природного газа как одного из исходных компонентов для реакции), свойств аппарата (образование сажи на поверхности реактора - «закоксовывание») и т.д., которые приводят к нарушению составленных ограничений и тем самым указывают на появление неполадок в технологической системе [3].

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

Xik(min) < х < Xik(Max), (i = 1,1; k = 1,K), (1)

где i - номер переменной состояния; k - номер технологического режима работы i-й переменной, с какой угодно степенью дифференциации, перечислять на основе чисто комбинаторного перебора возможные не пересекаемые области существования переменных состояния [2].

Определение центра технологической безопасности осуществляется в соответствии с алгоритмом, предложенным в [2], решением задачи нелинейного программирования (2) для определения отклонения координат данной точки от границ области безопасности: i = 1 —

Zd.(z) ^max. (2)

i=1 1

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

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

Работа выполнена при частичной финансовой поддержке гранта РФФИ «Исследование рисков при управлении динамическими процессами в слабоструктурированных и плохо формализуемых средах», проект № 17-07-01368.

Список литературы

1. Антонов В.Н., Лапидус А.С. Производство ацетилена. М.: Химия, 1970. 416 с.

2. Богатиков В.Н., Палюх Б.В., Пророков А.Е., Алексеев В.В. Приложение метода разделения состояний к управлению технологической безопасностью на основе индекса безопасности. Тверь: Тверской государственный технический университет, 2009. 256 с.

3. Санаева Г.Н., Пророков А.Е., Богатиков В.Н. О возможности применения дискретных моделей к диагностике состояния процесса получения ацетилена окислительным пиролизом природного газа // Энергетика, информатика, инновации: тезисы докл. V Междунар. науч.-тех. конф. (Смоленск, 2627 ноября 2015 г.). Смоленск, 2015. Т. 1. С. 225-228.

вых

i Надоели баннеры? Вы всегда можете отключить рекламу.