Математическое моделирование и оптимальное управление природо-промышленными системами

УДК 519.95+330.15+661.56

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ПРИРОДО-ПРОМЫШЛЕННЫМИ СИСТЕМАМИ

© А.П. Зубаков

Ключевые слова: математическое моделирование; производство азотной кислоты; природо-промышленные системы.

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

В рамках математического моделирования была сформулирована идея о необходимости управления природо-промышленными системами (ППС) [1, 2].

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

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

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

- формализация постановки задачи оптимального управления ППС;

- построение кинетических моделей реакций очистки;

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

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

- обеспечение совместимости моделей промышленной и экологической подсистем с учетом вероятностных факторов и ограничений, действующих на технологические и «экологические» переменные;

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

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

до двуокиси азота, охлаждаются в подогревателях воздуха 4, хвостовых газов 5, холодильнике-конденсаторе

6, где происходит конденсация водяных паров, образование неконцентированной азотной кислоты. Затем нитрозные газы поступают в колонну абсорбции 7, где происходит поглощение двуокиси азота паровым конденсатом с образованием продукционной кислоты, которая направляется в продувочную колонну 8 для отдувки окислов азота воздухом. Образующиеся здесь неочищенные хвостовые газы (НХГ) поступают на абсорбцию в подогреватель 5 и далее на установку каталитической очистки, являющейся объектом, определяющим количество окислов азота и углерода, выбрасываемых в окружающую среду.

Рассмотрим более подробно процесс очистки.

На рис. 2 изображена упрощенная схема камеры сгорания и реактора каталитической очистки. Поток природного газа (ПГ), подаваемый на основную горелку камеры сжигания, сгорает в потоке воздуха (Зона П).

Рис. 1. Схема процесса производства азотной кислоты

Рис. 2. Схема установки каталитической очистки

Продукты сгорания смешиваются с выхлопными газами (Зона IV) и природным газом на восстановление. В результате смесь разогревается до температуры зажигания катализатора. На катализаторе АПК-2 происходит реакция восстановления окислов азота (Зона VI). За счет тепловыделения смесь разогревается и поступает на второй слой катализатора (Зона Vin), где происходит реакция окисления СО до СО2. Очищенные хвостовые газы направляются в газотурбинную установку (ГТУ) 11 (рис. 1), в которой используется энергия расширения и охлаждения газов для выработки воздуха высокого давления в компрессоре 12 и нагнетателе 13. Избыток энергии преобразуется в двигателе-генераторе 14 в электроэнергию и выдается в сеть. Остаток энергии ОХГ используется в котле-утилизаторе 15 для выработки пара. Затем ОХГ через выхлопную трубу выбрасываются в атмосферу.

Особенностями технологии производства кислоты являются факты изменения нагрузки на процесс очистки по неочищенным хвостовым газам (НХГ), связанные с режимом отдувки в продувочной колонне; изменение активности катализатора (АПК-2) из-за его обратимой дезактивации сернистыми соединениями, присутствующими в природном газе. Управление процессом очистки НХГ также имеет особенности. В частности, увеличение расхода природного газа (ПГ) в реактор каталитической очистки приводит к образованию восстановительной среды, полной очистке от окислов азота. Однако при этом наблюдается повышенный выброс окиси углерода, проскок не прореагировавшего ПГ, что ведет к загрязнению атмосферы, прямым потерям. Высокие температура и давление, агрессивная среда в установке очистки затрудняют наблюдаемость процесса.

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

I*

Q

Q^ =йНаХГ ,

Y = ¥

C

Qmin - Q - Qmax ; Tmin - T - Tmax),

CE

Здесь I - критерий оптимальности, включающий затраты на очистку НХГ и экологический ущерб, нанесенный природе выбросами данного производства, рассчитанный по методике Госстроя. Уравнение (2) является условием обеспечения очистки для заданного расхода НХГ. Выражение (3) представляет математическую модель процесса каталитической очистки. Соотношение (4) есть модель процесса распространения газовых примесей в атмосфере, позволяющая рассчитать CE (X,, у, ) - значение приведенной концентрации примеси в узлах сетки на заданной области ю(х,,у,) eQ . Выражения (5) являются необходимыми ограничениями на управление технологическим объектом. Условие (6) выражает требование непревы-

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

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

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

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

Зона II

СН4 + 0,5 02 = СО + 2Н2; СО + 0,5 О2 = СО2;

Н2 + 0,5 О2 = Н2О;

Зона IV

СН4+ Н2О = СО + 3Н2; СН4 + 2О2 = СО2+2Н2О;

2Ш +2Н2 = N +2Н2О; 2С0+2Ы0 = 2СО2 + Ы2;

2СО + О2 = 2СО2 +0,5Ы2;

Зона VIII

2СО + О2 = 2СО2; 2С0+2Ы0 = 2СО2 + Ы2;

2С0+Ы02 = 2СО2 + 0,5Ы2;

В [3] подробно рассматривался кинетический механизм каждой брутто-реакции схемы. На основании концепций «квазистационарных» концентраций промежуточных продуктов реакций были получены уравнения для скоростей реакций. Константы рассматриваемых реакций были найдены из литературных источников, подробная библиография которых приведена в [4].

При этом значения некоторых констант уточнялись путем минимизации выражения:

П

Я = 2

і=1

где а1, а 2 - весовые коэффициенты; СЭ, Т^э - экспериментальные значения концентрации и температуры.

Исследование «гидродинамической» структуры технологического объекта методом трассерного эксперимента позволило осуществить выбор структуры математических моделей установки каталитической очистки. Были разработаны математические модели камеры сгорания и реактора каталитической очистки. Также разработана имитационная система прогноза распространения примесей (ИСПРП) для расчета уровней приземных концентраций загрязняющих веществ в исследуемой области.

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

ченные кривые РВП позволили сделать вывод о том, что зоны I, II, III камеры сгорания, зоны V, VII, IX реактора могут быть описаны моделью идеального вытеснения. Слой насыпного гранулированного катализатора хорошо (рис. 3) описывается моделью диффузионного типа. Зона IV конструктивно исполнена, как зона смешения. При построении моделей пренебрегали тепловыми потерями; газы, составляющие смесь, рассматривались как идеальные.

Ниже приведена математическая модель камеры сгорания по зонам (рис. 2):

Зона I

Ci (z) = Ci0 = const, i = 1,3;

T(z) = T0 = const,l1 < z < l1 +l2;

Рис. 3. Кривая РВП трассера в слое катализатора

Зона II

V2 ~Т~ +

dz

dT2 a-4(T2 -T3) А AHlri

-Y-

dz D2P IICPlI 7=1 Ц i P IICPu

= 0;

Зона III

dT3 a-4

= 0;

dz D2pIIICpIII

і1 + і2 + і3 < z < і1 + і2 + і3 + і4; Ci (z) = Ci0 = const, i = 1,3;

здесь Сг- - концентрация компонентов газовой смеси; Т - температура; V - скорость потока газа: Б - диаметр; р- плотность; Ср - теплоемкость; а- коэффициент теплопередачи; АН - энтальпия реакции; ц - молекулярный вес; I - длина; г - текущая длина.

Система уравнений (10) - (11) интегрировалась методом Рунге - Кутта с автоматическим выбором шага. Результаты расчетов, приведенные на рис. 4, показывают хорошее согласование с экспериментальными данными.

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

- тепловые потери незначительны, т. к. реактор конструктивно выполнен как адиабатический;

- каталитические реакции происходят только на поверхности катализатора;

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

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

Зона V

dCf

dz

T

dz

= 0; Cf = Clf (0);

=0;Tf0 = Tf (°);

Сі, мшь/мі*І0-4

T,K* 150

100

50

/ '

p ,

11 5/1 6

v ;

6fJ \ I V b

V2 3 .

10

О

Рис. 4. Распределение концентраций (1 - СН4, 2 - О2, 3 - СО, 4 - Н2, 5 - Н2О, 6 - СО2) и температур Т2 (Зона II), Т3 (Зона III) вдоль оси камеры сгорания

Зона VI

с 2 с# сс# к

- Уп-*- - вк (С# - с1к) = 0,

z dz2

i = 1,3;

dz

є

f

d 2Tf dTf

f - V f

* ТЛ

(PviCPvi)f dz2 V dz

Fk -a k

є- (PvicPvi ) f

(Tk - Tf) = 0;

F ___

Pk —

є

f

d 2Tk

Fk -a k

(pviCpvi)f dz (1 -є) • (pviCpVI)k

n

-Y (-AHr)-Ф-r (Cik, Tk) = 0;

i=1

Ф = f (Cs, T );

(Tk - Tf) -

ас#

= У (С г0 - С# (0)) = 0;

ат# # а,

-Ууі (руіСруІ)/ (Т0 - т (0)) = 0;| г = 0

1 =0; аг

ас#/ ат{

о*-Г = -х/~г = 0;|г = ІУІ; аг аг

Здесь В2 - эффективный коэффициент диффузии; X -эффективный коэффициент теплопроводности; р -коэффициент массопередачи; е - доля свободного объема катализатора; ф - функция влияния каталитического яда; / - индекс, показывающий принадлежность к газовой фазе; к - индекс, показывающий принадлежность к фазе катализатора.

Математическое описание VII, IX зон (рис. 2), представленное моделью идеального вытеснения, подобно описанию V зоны реактора. Структура модели для зоны VIII реактора, описываемой одномерной диффузионной моделью, аналогична модели зоны VI.

Решение уравнений (16) - (17) осуществлялось численно-аналитическим методом с «замороженными» правыми частями. Система уравнений (18) решалась методом Ньютона - Рафсона. Уравнение (19) было преобразовано к задаче Коши с последующим интегрированием методом Рунге - Кутта.

На рис. 5 приведены характерные профили температур и концентраций компонентов газовой смеси вдоль оси реактора. Максимальная относительная погрешность расчетов по модели в сравнении с данными экспериментов не превышает 6,4 % для реального диапазона изменения режимных параметров процесса, что позволяет использовать ее для целей управления.

Рис. 5. Распределение концентраций с#: (1 - СН4'0,1; 2 -

020,2; 3 - N00,1; 4 - N0^0,1; 5 - СО'0,1; 6 - 00^0,25; 7 -Н2'0,1; 8 - N210; 9 - Н20) и температур (Т/Т ) вдоль оси реактора

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

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

V = Бир1п/ К

Бц Бг Бг Бц

здесь Бц = [Р(цДР(ц2),...] - смешанная стратегия выбора модели из имеющейся в ИСПРП совокупности; Бр = [д(г1),^(г2),...,#(Гп)] - дискретные стратегии

смены классов стабильности атмосферы;

г=1 j=1

где К(Бц,Бг) - средняя эффективность прогноза; Еу - значение оценки точности прогноза по модели ц - при условии ее применения в классе гу Оценка

л

прогноза С , найденная с гарантированной точностью

V в условиях неопределенности смены класса стабильности атмосферы:

ои си

с = | с/£|/(с )ас /1/S|^(с )ас,

здесь /Бц (С) - неоднородная плотность вероятности, полученная в результате выбора оптимальной стратегии Бц . Данный подход особенно уместен в

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

г=1

где ру - значения дискретной вероятности установления определенного класса стабильности. Относитель-

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

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

Согласно постановке задачи оптимального управления ППС (1) - (6) выражение для расчета критерия оптимальности I:

І = (ЦПГбПГ + Уэк) ^ щіп,

где Ц ПГ - стоимость природного газа, Уэк - ущерб от выхлопных газов производства азотной кислоты. Ущерб рассчитывался по методике Госстроя СССР.

г=1

здесь Уэ - укрупненная величина удельного ущерба от 1 т условно-эквивалентного загрязнителя; Ву, Су - коэффициенты, учитывающие природно-климатические и социально-демографические условия.

а, =•

с, пдкс,

ПДК, ПДК г

управления ППС для нормальных условий рассеивания примесей в атмосфере.

В ходе выполнения рассматриваемой работы были достигнуты следующие результаты.

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

Рис. 6. Зависимость изменения экологического ущерба от количества израсходованного природного газа

где А, - коэффициент эквивалентности ,-го загрязнителя.

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

С 2 (Хх, у х ) рассматриваемых примесей значения

сПдк с заданной вероятностью, т. е.:

Р(с2 (ХХ, у *) < с 1дК ) > Р, X = 1, N.

Из условия малого загрязнения атмосферы, согласно действующим нормам (показатель загрязнения Р =

= 0,1), можно оценить Рзад (Рзад = 0,9).

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

На рис. 7 приведены оптимальные значения расходов природного газа на восстановление и сум-

глПГ

марного расхода природного газа б , расхода воз-

ґ~\В „НХГ ттл/'т-'

духа б при изменении нагрузки б по НХ1, полученные в результате решения задачи оптимального

£}в,

м3/ч

8000

7000

6000

5000

- /

I г ч. -

-

30000

40000

50000

1600

1200

800

400

Рис. 7. Оптимальные расходы природного газа на восстанов-

пПГ глПГ

ление бв , суммарного расхода б и расхода воздуха

гъ глИХГ

бв при изменении нагрузки по б

2. Разработаны кинетические модели, адекватно описывающие совокупность реакций очистки газов рассматриваемого производства.

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

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

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

6. Решена задача оптимального управления установкой каталитической очистки при различной нагрузке на нее, переменной экологической обстановке.

7. Разработано алгоритмическое, информационное и техническое обеспечение системы оптимального управления ППС.

8. Получены авторские свидетельства [6-8] на локальные системы автоматической стабилизации найденных оптимальных режимов управления процессом очистки.

ЛИТЕРАТУРА

1. Зубаков А.П., Толстых С.С., Скляревский А.М., Блох Б.М. Имитационная система прогноза распространения примеси в атмосфере // Проблемы охраны окружающей среды в районах с интенсивно развивающейся промышленностью: докл. II Республиканской конф. Кемерово, 1982. С. 17-18.

2. Зубаков А.П., Попов Н.С., Громов Ю.Ю. Система имитационного моделирования и прогноза распространения примесей в атмосфере // Проблемы разработки автоматизированных систем наблюдения, контроля и оценки состояния окружающей среды: докл. II Всесо-юз. науч.-техн. конф. Казань, 1983. С. 12-13.

3. Зубаков А.П., Попов Н.С. Математическое моделирование процесса каталитической очистки нитрозных газов в задаче оптимального проектирования // Автоматизированное проектирование в химическом машиностроении: сб. науч. тр. / МИХМ. М., 1984. С. 7781.

4. Зубаков А. П. Оптимальное управление очисткой нитрозных газов производства азотной кислоты с учетом экологических ограничений: автореф. дис. ... канд. техн. наук. М.: Моск. ордена Трудового Красного Знамени ин-т химич. машиностроения, 1987. 16 с.

5. Зубаков А. П. Моделирование процессов распространения примесей в атмосфере при информационной неопределенности природных факторов // Актуальные проблемы информатизации и информ. технол.: материалы 2 Тамб. межвуз. науч. конф. Тамбов, сент. 1998 г. Тамбов, 1998.

6. А. с. № 1269429 (СССР) Способ регулирования процесса очистки нитрозных газов от окислов азота в производстве слабой азотной кислоты / Зубаков А.П., Бодров В. И., Попов Н.С., Кравченко С.П., Иванова Н.Н., Скляревский А.М. М., 1985.

7. А. с. № 1234356 Способ регулирования процесса очистки хвостовых газов от окислов азота в производстве азотной кислоты / Зуба-ков А.П., Бодров В.И., Попов Н.С., Шарапов А.Н., Иванова Н.Н. М., 1986.

8. А. с. № 1680619 Способ управления процессом каталитической очистки нитрозных газов в производстве неконцентрированной азотной кислоты / Зубаков А.П., Бодров В.И., Кравченко С.П., Скляревский А.М. М., 1991.

Поступила в редакцию 26 марта 2009 г.

Zubakov A.P. Mathematical modeling and optimum control of natural-industrial systems. In the work the approach to optimum control of the natural-industrial systems based on a use of methods of mathematical modeling is considered. As an example usage of the given approach to nitric acid manufacture is considered.

Key words: mathematical modeling; nitric acid manufacture; natural-industrial systems.