Проектирование оптимальной подсистемы реакторного узла процесса изомеризации н-пентана путем решения двухэтапной задачи оптимизации Текст научной статьи по специальности «Математика»

УДК 66.01

И. В. Зайцев, Т. В. Лаптева, Н. Н. Зиятдинов

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

ИЗОМЕРИЗАЦИИ Н-ПЕНТАНА ПУТЕМ РЕШЕНИЯ ДВУХЭТАПНОЙ ЗАДАЧИ ОПТИМИЗАЦИИ

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

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

Key words: optimal process design, optimization under uncertainty, two-stage optimization problem, chance constraints.

Optimal process design problems are concerned with various complexities: taking into account of source data uncertainty and calculation of multidimentional integrals, which underlie the criterion and constraints functions, are necessary. In previous studies we considered approaches and algorithm for solving two-stage optimal design problems with chance constraints. In this paper reactor unit subsystem of n-pentane isomerization process optimal design problem will be solved by these approaches.

Введение

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

Предложенные нами подходы были успешно опробованы на задачах проектирования разных ХТС [1, 2]. В данной работе с помощью разработанных нами подходов и алгоритмов решения задачи двухэтапного оптимального проектирования ХТС [3, 4] будет решена задача оптимального проектирования подсистемы реакторного узла, который входит в ХТС изомеризации н-пентана [5], с учетом имеющейся внутренней и внешней неопределенности в исходной информации. Для решения задачи будет использован разработанный нами программный комплекс решения задачи двухэтапной оптимизации с вероятностными ограничениями, которая имеет вид [6]:

maxg (d,z(0),0) < 0, j = m +1,m + q,

(3)

min EJf(d,z(0),0);T]

d, z(6) в

Pr{g,(d,z(0),0) < 0} >a., j = im,

(1) (2)

где б - вектор конструктивных параметров размерности пб , г(0) - вектор-функция с компонентами

гр (0), р = 1, п2 , 0 - вектор неопределенных параметров размерности п0, f (б, 1(0), 0) - критерий

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

количество проектных ограничений. Ограничения (2) - вероятностные ограничения, которые выполняются с некоторой заданной вероятностью а, (3) - жесткие ограничения, которые выполняются безусловно на всей области неопределенности Т .

Задача (1) является задачей полубесконечного программирования [7, 8]. Нами были рассмотрены подходы и методы, позволяющие решать такие задачи [9]. Для решения задачи (1) мы использовали модифицированный метод внешней аппроксимации [8]. Также нами были рассмотрены разные виды аппроксимации функций г(0) зависимостей управляющих параметров от неопределенных [1].

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

Технологическая схема узла изомеризации н-пентана представлена на рисунке 1. Н-пентановая фракция смешивается с водородосодержащим газом (ВСГ) в смесителе С , нагревается в рекуператоре Т до температуры 306°С за счет теплоты реакционных газов. Затем газосырьевая смесь (ГСС) нагревается в печи П до температуры 360 - 450°С . Процесс изомеризации н-пентана протекает в реакторе в неподвижном слое алюмоплатинового катализатора ИП-62 ВК в среде ВСГ при общем давлении 3,6МПа . Контактный газ (КГ) из реактора Р охлаждается в рекуператоре Т и далее поступает на отде-

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

Рис. 1 - Технологическая схема реакторного узла процесса изомеризации н-пентана

Состав н-пентановой фракции и ВСГ, поступающих на переработку в реакторный узел, приведены в табл. 1. Содержание водорода в ВСГ колеблется, диапазон изменения концентрации водорода в ВСГ приведен в таблице 2.

Таблица 1 - Состав н-пентановой фракции и ВСГ

Водород Азот Метан Этан

Сырье, % вес 0 0 0 0

ВСГ, % вес 93,8 0,3 4,3 0,44

Пропан Изопентан Н-Пентан Гексан

Сырье, % вес 0 3,19 95,96 0,85

ВСГ, % вес 0,17 0,7 0,29 0

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

Математическая модель реактора

Реакция высокотемпературной изомеризации н-пентана на алюмоплатиновом катализаторе ИП-62 ВК [11] протекает через стадии дегидрирования н-пентана в олефин, изомеризации олефина в изоолефин, гидрирования изоолефина в изопентан.

Учитывая, что лимитирующей стадией реакции является изомеризация олефина [12] и что, кроме основной реакции, процесс сопровождается необратимым гидрокрекингом н-пентана и изопен-тана в легкие углеводороды (метан-бутановая фракция), упрощенный механизм реакции может быть представлен в виде

нС5

иС5

(4)

нСк

-+а1СГ ■

иС

-+а2СГ.

5 ' Г ' " 5

где нС5 , иС5 , Сг обозначают н-пентан, изопентан,

^, М,

3 '

- скорости реакций; ау, а2 - стехиометриче-

ские коэффициенты. Реакции крекинга играют не-

значительную роль вследствие того, что реакция

изомеризации протекает с малым тепловым эффек-

том и высокой селективностью. Этим обусловливается низкий суммарный тепловой эффект процесса, составляющий 84 кДж / кг сырья. Скорости реакций изомеризации и побочных реакций описываются уравнениями [12]

М = РНг , М2 = 1<2Рис5/Рн2 , (5)

М3 = <3РнС5 , М4 = <4рис5 ,

где РнС5, РиС, Р„2 - парциальные давления н-пентана, изопентана и водорода соответственно, МПа ; <0, / = 1,4 - кажущиеся константы скоростей реакций, зависимость которых от температуры выражается уравнением Аррениуса

<, = <0 0ехр{-Е ,/ЯТ), (6)

в котором Я - газовая постоянная, кДж / (моль ■ К); Т - абсолютная температура, К ; Е - энергии активации, кДж / моль ; величины < 10, Е и Я равны [5]:

<20 = 1,98 ■ 1013 кмоль

<30 = 1,56 ■ 10

<40 = 4,69 ■ 10

ч ■ м кмоль ч ■ м3 ■ МПа 17 кмоль

ч ■ м3 ■ МПа = 2

Е. = 217 , Я = 8,3 ■ 10-3.

Е1 = 172 , Е2 = 180 , Е3 = 217 ;

Значения величин <20, <30, <40 вычислены в работе [5], на основе идентификации модели промышленного реактора и включают в себя погрешности измерения, ошибку идентификации и не учитывает изменение активности катализатора со временем. Поэтому будем считать эти значения неточными. Диапазоны изменения этих неопределенных параметров приведены в табл. 2.

Таблица 2 - Диапазоны изменения неопределенных параметров

Параметр $ в" ви

<20 1,9206-1013 1,98-1013 2,0394-1013

<30 1,5132-1017 1,56-1017 1,6068-1017

<40 4,5493-1017 4,69-1017 4,8307-1017

с„2, % вес 90 92,5 95

При математическом описании реактора приняты следующие допущения:

1) аэродинамическая обстановка близка к идеальному вытеснению;

2) тепловой режим адиабатический;

3) радиальный перенос вещества и тепла отсутствует.

Уравнение материального баланса по ключевым компонентам н-пентану и изопентану [5]

СИ

А ( (-<Х„С5 + ХиС5 ) - <3РНг Х„С5), (7)

КГ

dx

= A ^ - xцC5) - ^ xцC5). (8)

Уравнение теплового баланса

dT A

— = — (Q1(k2(kxяc5 - XЦC5 ))

dL Cp 55 (9)

+ Q2k3PH2 x„c5 + Q3к4PH2 xцc5)

При L = 0 : xяc5 = ^ , ^ = ^ , T = Т.

При гидрокрекинге одной молекулы н-пентана или изопентана образуются две молекулы легких углеводородов (а, а2 в уравнении (4) равны 2), и на их образование расходуется соответственно одна молекула водорода.

Уравнения материального баланса по легким углеводородам и водороду [5]

ДХг = 2 • ((^ - х^) - (х^ - Xc)), (10) xг = xlг + Дxг , (11)

Xн2 = xH2 -Даг/2 (12)

В уравнениях (7) - (9) множитель A определяется из выражения

А = Скат¡Т/СР^ , (13)

С = Сс + Свсг -Асг. (14)

Константа равновесия реакции изомеризации к = к1/к2 аппроксимируется уравнением

к = 2,89 + 0,0044(600 - Т). (15)

Конверсия (C, %) и селективность (5 , %) реакции определяются по уравнениям

С = (х^ - • 100/х^ , (16)

5 = (xцc5 - <) • 100/(х^ - х^). (17)

В уравнениях: хнС , хиС , хг , хн - мольный поток

н-пентана, изопентана, продуктов гидрокрекинга, водорода соответственно (верхний индекс I относится к начальным условиям); Q1, Q2, Q3 - тепловые эффекты скоростей реакций изомеризации, гидрокрекинга н-пентана и изопентана, кДж / моль ; L - текущая относительная длина реактора (0 < L < 1); Cp

- теплоемкость газосырьевой смеси, кДж /(моль • К); С - суммарный поток газосырьевой смеси, кг / ч ; Скат - весовая загрузка катализатора, кг ; Сс - поток сырья в реактор, кг / ч ; Свсг -поток водородосодержащего газа в реактор, м3 / ч ; рвсг - плотность водородосодержащего газа, кг / м3.

Таким образом, уравнения (6-17) составляют математическую модель реактора.

Математическая модель смесителя

Математическая модель предназначена для учета смешения потоков в емкостях и трубных соединениях. Тепловой эффект смешения не учитывается. Количество С° и состав х° выходного потока рассчитывается по формулам

(18)

в° = 21С ,

где х° 1 - содержание I -го компонента в } -ом выходном потоке, вес. доли, - количество } -го

входного потока, кг /ч .

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

^ = 2/(С°с°), (20)

где Ц - температура } -го входного потока, К, с'р] - теплоемкость } -го входного потока, с° - теплоемкость выходного потока, кДж / (кг • К).

Уравнения (18-20) составляют математическую модель смесителя.

Математическая модель рекуператора

Математическая модель предназначена для вычисления выходных температур потока ГСС (поток 3) и потока реакторных газов (поток 6), которые участвуют в теплообмене в рекуператоре Т. Температура реакторных газов ^ и ГСС tg на выходе из рекуператора рассчитываются по формулам [5]

tR = -Ф(^ - (?), (21)

Г-д = ^ +ФА(^ - ^), (22)

1

Ф =

г + А +1 г

+ -

2

5 = -

взг -1 кР

А =

= 4(1 + А)2 - 4рА

срг

' " Сдс]д

где t - температура теплоносителя, ° C , Ф - функция тепловой эффективности, А г - функции водяных эквивалентов, 5 - число единиц переноса тепла, С - расход теплоносителя, кг/ч, Р - площадь теплопередачи, м2, т] - коэффициент потерь, к -коэффициент теплопередачи, кДж/ м2, р = 1 - индекс противоточности, верхние индексы Я - поток реакторных газов, д - поток ГСС, нижние индексы

I, о - входной и выходной параметры потока.

Математическая модель нагревательной печи

Математическая модель предназначена для расчета расхода топливного газа в печи Стг, кг / ч, обеспечивающего заданный перепад температуры нагреваемой ГСС. Расход топливного газа определяется из уравнения теплового баланса печи

Стг =

Сср - ?)

(23)

х° =2/ С°

(19)

ДН

где ср - средняя теплоемкость смеси в интервале температур - ^), кДж /(кг • К), ДН - теплотворная способность топлива, кДж / кг , т] - КПД печи, С -количество ГСС, поступающей в печь, кг / ч . Постановка задачи

В качестве неопределенных параметров выбраны параметры математической модели, которые послужили настроечными параметрами в процессе идентификации математической модели реактора.

Это предэкспоненциальные коэффициенты в урав-

Л I кмоль

нениях Аррениуса (6), k20, -

k3

ч • м

кмоль

k4

кмоль

•. Также в число неоп-

ч • м3 • МПа ' 4 ч • м3 • МПа

ределенных параметров включили концентрацию водорода в ВСГ, сн , %, поскольку водород расходуется на гидрокрекинг, теряется в отдувках ВСГ, циркулирующего в ХТС изомеризации н-пентана, что необходимо для соблюдения молярного соотношения водород :н-пентан, потери водорода происходят при растворении водорода в изомеризате. Кроме того, ВСГ накапливает в себе побочные продукты реакций. Таким образом вектор неопределенных параметров представляет собой

^ = {k20 , ^ k40 , СН2 }.

Область неопределенности задана в виде T = {в, : в\ < в, < в",i = 1,...,4}. Значения в\ и в" ,

i = 1,4 , приведены в таблице 2.

В качестве поисковых переменных задачи выбраны:

1) масса загрузки реактора катализатором, GKaT, кг , - конструктивный параметр d ;

2) управляющие параметры г(в) :

- температура на входе в реактор, Т1, °C ;

- расход электролитического водорода, GH , кг / ч .

В качестве ограничений задачи проектирования оптимального реактора узла изомеризации н-пентана использованы следующие мягкие ограничения:

1) вероятностное ограничение на конверсию

Pr{0,55 - conv(d, 2{в)в) < 0} > 0,95, (24)

2) вероятностное ограничение на селективность

Pr {0,99 - sel(d, z(e),e) < 0} > 0,95, (25)

3) вероятностное ограничение на соотношение водород: н- пентан

Pr{2 - CM(d,z(e),e) < 0}> 0,95 . (26)

А также жесткое ограничение на соотношение толщины стенки, диктуемое [13]

0,001DB < s(d,z(e),e) < 0,1DB + ^ • ^экспл Рассмотрим теперь ограничения на управляющие переменные. Для расхода электролитического водорода границы изменения составят

0 < GH2 (в) < 100 , (27)

для температуры [11]

350 < Т1 (в) < 450. (28)

Ограничения на управляющие переменные являются жесткими ограничениями. Хотя в силу ограничений (24-26) заданная вероятность выполнения соответствующих проектных требований будет соблюдена, при этом теоретически возможны точки области неопределенности, в которых значения ограничений могут не удовлетворять возможным требованиям на качество продукта. В связи с этим добавим следующие жесткие ограничения на конверсию и селективность:

0,97 - sel(d,z(e),e) < 0 . (30)

Ограничения (29), (30) гарантируют выполнение нижней границы проектных требований на конверсию и селективность на всей области неопределенности.

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

C(d, г{в),в) = (ЦаДат + ЦстАуст (d, z(6),6))/ 80000 +ЦG(d,z(ß),в) + ЦН2GH¿ (d,z(0),0), (31) где Цкат - цена катализатора, руб / кг, GKaT - масса катализатора, кг , Цст - цена стали, руб / кг, тст -необходимая для постройки реактора масса стали, кг, ЦН2 - цена электролитического водорода,

руб / кг, wВС2г - весовое количество водорода, содержащееся в ВСГ, кг / ч, GH2 - расход водорода, кг / ч, G^ - расход топливного газа в печи, кг / ч, Цтг - цена топливного газа, руб / кг.

Исходя из характеристик катализатора ИП-62 ВК [11], определяем срок эксплуатации системы 4кспл = 10 лет . С учетом проведения ежегодных остановов на планово-предупредительный ремонт среднее время работы ХТС 8000 ч / год или 80000 ч за все 10 лет работы системы.

Итак, задача оптимального проектирования реактора узла производства изопентана формулируется следующим образом:

Найти такие значения поисковых переменных GKan , Т1 (в), GH2 (в), при которых значение математического ожидания критерия (31) на области неопределенности T будет минимальным, а также будут выполняться ограничения (24-26) с заданной вероятностью и ограничения (27), (28) на всей области неопределенности T .

Мы будем рассматривать постановку задачи проектирования оптимального реакторного узла в виде двухэтапной задачи оптимизации с жесткими и мягкими ограничениями (1). Формализованная постановка задачи примет вид

min EJC(d,z(e),e);T] (32)

G„ ,z(в)

Pr{0,55 - conv(d,z(e),e) < 0} > 0,95, Pr {0,99 - sel(d,z(e),e) < 0} > 0,95 , Pr{2 - CM(d,z(e),e) < 0} > 0,95,

0 < GH2(e) < 100 , Уве T,

350 < Т (в) < 450, Уве T, 0,52 - conv (d, z^^) < 0 , Уве T, 0,97 - sel (d, z^)^ < 0 , Уве T, 0,001DB < s(d, z^) < 0,1DB + ^ • 4кспл, Ув е T, где z(^ = {Gh2(в); Т'(в)}.

0,52 - conv(d,z^)^ < 0:

(29)

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

Результаты

Используя подходы и алгоритмы описанные в настоящей работе была решена задача оптимального проектирования (32), поставленная в настоящей работе. Результат решения приведен в таблице 3.

Таблица 3 - Результаты решения задачи (32)

С - С

Д = ^-д

Подход Проектирование на основе учета коэффициента запаса Проектирование на основе ДЭЗО

Загрузка кат., Скат , кг 15884,7 13507

Высота реактора, Л, м 11,285 10,691

Диаметр реактора, d, м 1,881 1,782

Толщина стенок, 5, мм 28,04 26,6

Приведенные затраты, руб / год 17 070 863 14 637 134

Относит. экономия, % 14,26

Оценка экономии за период экспл., руб - 24 337 290

Полученные результаты сравнивались с подходом проектирования, который традиционно применяется при проектировании аппаратов химической технологии. При данном подходе сначала решается задача номинальной оптимизации без учета имеющейся неопределенности при значениях неопределенных параметров в, = в^, затем полученные значения конструктивных параметров умножают на соответствующий эмпирический коэффициент запаса узап (для химической промышленности обычно принимают узап = 1,25 [14]).

Из таб. 3 видно, что учет неопределенности с помощью двухэтапной задачи оптимального проектирования позволил сократить затраты (31) на 14,26% . Относительная экономия вычисляется по формуле

Сз

где Сз - затраты (31) при использовании коэффициента запаса узап, Сд - затраты (31) при предлагаемом нами подходе проектирования на основе ДЭЗО.

Таким образом, оценка экономии за весь период эксплуатации реактора составляет порядка 24,34 млн руб . Решение задачи заняло около трех минут.

Литература

1. Н.Н. Зиятдинов, И.В. Зайцев, Т.В. Лаптева, Вестник Казан. технол. ун-та, 15, 16, 247-250, (2012).

2. Н.Н. Зиятдинов, Г.М. Островский, И.В. Зайцев, Т.В. Лаптева, Л.Р. Хисамутдинова. Вестник Казан. технол. ун-та, 14, 10, 223-231 (2011).

3. И.В. Зайцев, Т.В. Лаптева, Н.Н. Зиятдинов, Вестник Казан. технол. ун-та, 16, 1, 251-256, (2013).

4. G.M. Ostrovsky, N.N. Ziyatdinov, T.V. Lapteva, Chem. Eng. Science, 65, 2373-2381, (2010).

5. Г.М. Островский, Ю.М. Волин, Н.Н. Зиятдинов. Методы оптимизации химико-технологических процессов. КДУ, Москва, 2008. 424 с.

6. G.M. Ostrovsky, N.N. Ziyatdinov, T.V. Lapteva, I.V. Zayt-sev, Chem. Eng. Science, 66, 3815-3828, (2011).

7. M. Lopez, G. Still, European J. Oper. Res, 2, 180, 491-518, (2007).

8. F. Guerra-Vazquez, J.-J. Ruckmann. Semi-Infinite Programming. - Wiley Encyclopedia of Operations Research and Management Science, (2011).

9. Т.В. Лаптева, Н.Н. Зиятдинов, И.В. Зайцев, Вестник Казан. технол. ун-та, 15, 24, 139-146, (2012).

10. Г.М. Островский, Н.Н. Зиятдинов, Т.В. Лаптева. Оптимизация технических систем. КНОРУС, Москва, 2012. 432 с.

11. Катализатор изомеризации н-пентана ИП-62 ВК [Электронный ресурс] / Олкат, 2013. - Режим доступа: http://www.olkat.ru/ip-62vk

12. Н.Р. Бурсиан, Технология изомеризации парафиновых углеводородов. Л., Химия, 1985. 192 с.

13. ГОСТ Р 52857.2-2007. Сосуды и аппараты. Нормы и методы расчета на прочность. Расчет цилиндрических и конических обечаек, выпуклых и плоских днищ и крышек. - М.: Стандартинформ, 2009. - 41 с.

14. С.И. Дворецкий, А.Ф. Егоров, Д.С. Дворецкий, Компьютерное моделирование и оптимизация технологических процессов и оборудования. Учебное пособие. Изд-во ТГТУ, Тамбов, 2003. 224 с.

© И. В. Зайцев - асс. каф. системотехники КНИТУ, [email protected]; Т. В. Лаптева - к.т.н., доцент той же кафедры; Н. Н. Зиятдинов - д-р техн. ноаук, проф., зав. . каф. системотехники КНИТУ, [email protected].