СИСТЕМНЫЙ АНАЛИЗ X
УДК 551.46.08
МОДЕЛИРОВАНИЕ ЗАМКНУТОЙ СИСТЕМЫ УПРАВЛЕНИЯ «ПРИРОДА-ТЕХНОГЕНИКА»
Р. И. Сольницев,
доктор техн. наук, профессор Г. И. Коршунов, доктор техн. наук, профессор А. А. Шабалов,
ассистент
Санкт-Петербургский государственный университет аэрокосмического приборостроения
Рассмотрены вопросы создания замкнутой системы управления «Природа-техногеника», предназначенной для эффективного снижения загрязняющих веществ, выбрасываемых промышленными предприятиями в атмосферу. Представлена математическая модель замкнутой системы управления и модели составляющих ее звеньев. Приведены результаты анализа системы на основе моделирования.
Концепция «Природа-техногеника», предложенная и развитая в работах [1-6], предусматривает создание замкнутой системы управления «Природа-техногеника» (ЗСУПТ) для широкого класса объектов, где требуется минимизация концентрации загрязняющих веществ (ЗВ) техногенного характера.
Отличительные особенности концепции представлены тремя критериями:
— минимизация или полное исключение «человеческого фактора»;
— управление концентрациями ЗВ на основе определения их максимальных значений;
— сохранение технологии основного производства в условиях минимизации концентраций ЗВ.
В работе рассмотрена задача моделирования ЗСУПТ «Природа-техногеника», решение которой необходимо для создания таких систем в различных областях человеческой деятельности.
Объектом управления ЗСУПТ является масса (концентрация) ЗВ в газовоздушной смеси, перемещающаяся от устья «трубы» промышленного предприятия до точки, где расположен измерительный преобразователь. В работе рассматривается ЗСУПТ на примере только одной составляющей ЗВ — диоксида серы Я02.
Ущерб, наносимый окружающей среде диоксидом серы, его характеристики и составляющие хорошо известны. Дополнительно можно отметить, что при производительности котла ТЭЦ 100 т пара/ч выброс Я02 составит 2-3 т/сут, и выпадение серы будет иметь место в круге с радиусом до 100 км от ¿-го источника, что крайне отрицательно сказывается на здоровье людей, растениях, зданиях и других составляющих окружающей среды этой ТЭЦ.
Задачей математического моделирования является воспроизведение на ЭВМ динамики и процессов функционирования ЗСУПТ (рис. 1) и дальнейший
■ Рис. 1. Схема замкнутой системы управления «Природа-техногеника»
анализ на основе моделирования. Для решения этой задачи необходимо составить математические модели всей системы и ее отдельных звеньев.
Для качественной оценки динамики ЗСУПТ рассмотрим уравнение материального баланса, выражающее неразрывность изменения массы или концентрации ЗВ:
Лу
Лі
= КіХ - К2 у - К у,
(1)
где х — среднее значение количества сжигаемого в ТЭЦ топлива при производстве основного продукта — пара; у — среднее значение количества массы ЗВ, выбрасываемого ТЭЦ в атмосферу; К1 — коэффициент пропорциональности, определяющий содержание в топливе ЗВ (в частности, содержание серы в сланцах составляет 3,5 %, мазуте — 2,5 %, угле — от 0,5 до 4 %); К2 — коэффициент, зависящий от синоптических, температурных, химических и других внешних параметров и определяющий поглощение Я02 внешней средой; К3 — коэффициент, выражающий величину позиционной обратной связи системы управления в ЗСУПТ. Уравнение (1) имеет аналитическое решение
(
У =
Уо
К
К2 + Кз
Х
е-( К+Кз)і +_
К
К2+К
-х, (2)
где у0 — масса ЗВ в момент времени t0 = 0.
Из этого следует очевидный результат: выбором регулируемого коэффициента К3 можно принципиально снизить массу (концентрацию) ЗВ до предельно допустимой величины [5].
В случае заданных как функции времени х(^ и К2(^ уравнение (1) примет вид
Лу
Лі
+ [К(і) + Кз]у = КіХ(і).
(3)
Решение этого уравнения также находится аналитически двумя квадратурами:
-} [ К2(і)+Кз]Лі
У = Уое
+
-\ [К2(і)+Кз]ЛіТ / [ К2(і)+Кз]Лі
+ е 0 | К1х(і)е0 Лі,
о
(4)
где Т — время наблюдения.
Выражение (4) позволяет контролировать результаты моделирования и получать качественные оценки ЗСУПТ. Дальнейшее уточнение предложенной модели достигается содержательным представлением коэффициентов К1 и К2, которые определяются через связи выброса с расходом топлива (К1) и с синоптическими, диффузионными, химическими и конвекционными процессами (К2). Приближенное значение К2 можно определить двумя аддитивными составляющими:
К2 = К2ф + К2£; К2ф = 0,6...0,9;
Кгъ=-
н2
2а?
где К2ф — составляющая от фильтрации ПГУ; — время выброса; Н — высота выброса (высота трубы и высота факела выброса); V — скорость ветра; ау, аг — среднеквадратические отклонения факела от его оси.
Переменные х(і) и К2(і) обычно известны по результатам мониторинга.
Для оценки характеристик ЗСУПТ по устойчивости, качеству, динамике процессов регулирования предлагается математическая модель ЗСУПТ на основе моделей отдельных звеньев, представляющих собой передаточные функции звеньев ЗСУПТ, показанной на рис. 1.
Структурная схема ЗСУПТ согласно рис. 1 изображена на рис. 2.
Приведем математические модели звеньев ЗСУПТ в виде передаточных функций отдельных звеньев.
Передаточная функция переноса ЗВ с выхода топливных газов до устья трубы (W1) представлена выражением
У1(і) = ^(р)о(і), Wí =
К
(5)
где р =— при нулевых начальных условиях; Лі
а(і) = х(і) - и*(і),
(6)
где х — возмущающее воздействие; и — управляющее воздействие.
Передаточная функция преобразования массы ЗВ в составе факела при переносе от устья трубы до точки измерения (W2) в простейшем случае может быть представлена выражением
■ Рис. 2. Структурная схема ЗСУПТ ИНФОРМАЦИОННО-УПРАВЛЯЮШИЕ СИСТЕМЫ 37
№ 2, 2008
У2(t) = W2(p)y1(t), W2 = -К—е-Р\ (7)
1 + Т2 р
где К4 — коэффициент пропорциональности, зависящий от превращений параметров атмосферы, таких как давление, влажность, начальная концентрация ЗВ, температура окружающей среды, коэффициенты диффузии и поглощения ЗВ подстилающей поверхностью; Т2 — постоянная времени инерционных процессов переноса ЗВ; т1 — время чистого запаздывания при переносе количества вещества от устья трубы до точки измерения ЗВ датчиком. В общем случае величина т1 определяется в результате решения уравнения турбулентной диффузии при заданных краевых и начальных условиях. В простейшем случае при расчете т1 значениями горизонтальных коэффициентов диффузии можно пренебречь, а движение воздушных потоков считать однородными в рассматриваемой области пространства. В этом случае т1 рассчитывается по формуле
Т1 = S/v, (8)
где S — расстояние от устья трубы до точки измерения максимума концентрации ЗВ датчиком; V — средняя скорость ветра на высоте факела (Н).
Передаточная функция измерительного преобразования (W3) представлена выражением
У* = те3(Р)У2, Wз =-К^-е-рТ2, (9)
1+Т3 р
где т2 — время измерительного преобразования, зависящее от принципа действия и конструкции используемого измерительного устройства. Для предполагаемых к использованию полупроводниковых сенсоров время реакции газочувствительной характеристики не превышает 5 с и время восстановления также не более 5 с. Измерительный преобразователь располагается в точке максимума концентрации ЗВ и выполняет полный цикл измерения за время, определяемое параметрами сенсора.
Передаточная функция процесса формирования, накопления, обработки, преобразования и передачи данных измерения (W4) представлена выражением
У*>) = W^(рШ, W4 =-К^~е~РТ3, (10)
1 + Т4 р
где т3 — суммарное время, затрачиваемое на процессы формирования, накопления, обработки, преобразования и передачи данных измерения ЗВ; 2 — ошибка, вычисляемая по формуле
^ = У2 - У 20, (11)
где у20 — допустимая величина концентрации, в частности ЗВ, на расстоянии S от устья трубы
при определенных погодных условиях и параметрах функционирования предприятия (ТЭЦ), которая всегда меньше ПДК.
Параметры передаточной функции W4 определяются быстродействием аналого-цифрового преобразования и накопления в памяти информационно-вычислительного устройства, а также временным запаздыванием, зависящим от скорости передачи, времени обработки и накопления данных, которое может составлять единицы или десятки секунд. Выбор стандартных протоколов GSM и GPS обеспечивает стабильность передачи данных, однако возможно использование специальных протоколов.
Управляющее воздействие на агрегат очистки от загрязняющих веществ (W5) представлено выражением
u = (K + Kzp + K9/p)y**. (12)
Значения коэффициентов передачи по пропорциональной составляющей, по первой производной и изодромной составляющей выбираются из условий устойчивости и качества регулирования концентрации ЗВ. Синтез управлений в ЗСУПТ приведен в работах [3, 5]. Здесь использован ПИД-регулятор. Применение других управлений — инвариантных, импульсных, экстремальных — предусмотрено при дальнейшей разработке ЗСУПТ.
Передаточная функция очистного агрегата (W6), включающего усилительно-преобразующее устройство, представлена выражением
u*(t) = W6(p)u(t), W6(p) = -^4 (13)
1 + T5 p
где K10, T5 зависят от типа и конструкции управляющих ПГУ.
Время реакции для различных типов ПГУ (рукавных, электрофильтров, импульсных, каталитических) составляет от единиц до десятков секунд.
На основе математических моделей (5)—(13) строится полная математическая модель (передаточная функция) ЗСУПТ, оценивается ее устойчивость, динамические и статические ошибки при вариации Ki, Tj, ts, i = 1 ^10, j = 1^3, s = 1^3.
Передаточная функция ЗСУПТ по ошибке регулирования z/y20 имеет вид
г = p(1 + T p)(1 + T2 p)(1 + T3 p)(1 + T4 p)(1 + T5 p) (14)
У20 p(1 + T1 p)(1 + T> p)(1 + T3 p)(1 + T4 p)(1 + T5 p) + a ’
где a = K1K4K5КбКю(К9 + K7p + K8p2)e~p(Ti+T2+4
Передаточная функция по отношению к возмущению x:
z K1K4 K5e~ p(T1+T2) p(1+T4 p)(1+T5 p) (15)
x p(1+T1 p)(1+T2p)(1+T3p)(1+T4p)(1+T5p)+a'У )
Наименование X К1 Т1 *1 К4 Т 2 *2 К5 Т3 *3 к6 Т4 * У 20 К7 К8 К9 К10 Т5
Размерность т - с с - с с - с с - с т - - - - с
Значение 1000 0,01 1,6 60 1 10 6 1 1 0 1 1 1 2 0 0,35 1 10
Из (14) и (15) следует, что ЗСУПТ является астатической по отношению к ошибке регулирования У*20 и по отношению к возмущению х.
В таблице приведены номинальные значения параметров, полученные из различных источников [1-6]. На рис. 3 показано изменение ошибки 2 при ступенчатом входном воздействии функции X для разомкнутой (кривая 1) и замкнутой (кривая 2)систем.
На основе данных таблицы можно показать, что процессы измерения выполняются в режиме реального времени и параметры интегратора обеспечивают устойчивое регулирование при периодическом получении данных измерения по каналу связи. В дальнейшем осуществлялось моделирование поведения системы, заданной передаточной функцией, при изменении значений отдельных параметров. По результатам моделирования определены параметры управления, обеспечивающие устойчивость, требуемые величины затухания, динамической и статической ошибок.
Рисунки демонстрируют изменение процесса регулирования от изменения:
— величины чистого запаздывания при переносе количества вещества от устья трубы до точки измерения ЗВ датчиком (рис. 4, а);
— вида топлива (рис. 4, б);
■ Рис. 3. Изменение ошибки 2 при ступенчатом изменении воздействия функции Х, с замкнутой и разомкнутой обратной связью
— параметров интегратора (рис. 4, в);
— выбранных комбинаций величины чистого запаздывания и параметров интегратора (рис. 4, г);
— коэффициента усиления регулятора (рис. 4, д).
Анализ зависимостей позволяет ввести допуски на значения параметров. Результаты моделирования показывают возможность создания эффективной системы управления, способной обеспечивать достаточное качество управления концентрацией ЗВ в газовоздушной смеси.
В дальнейшем учет диффузионных и конвекционных процессов в ЗСУПТ осуществляется на основе уравнения классической турбулентной диффузии:
дУд + „ дУд + „ 3У + „ дУд =
-----+ Vг---+ V----+ Vр---—
дt 1 зс п дп 5 35
—+±кдУ±+±кдУ±-„У, (16)
3 с ^ 3 с дп п 3п 35 5 35 у
где Уд — концентрация ЗВ; к,, кп, ^ — горизонтальные и вертикальная составляющие коэффициента турбулентного обмена; V,,, Vп, v^ — горизонтальные и вертикальная составляющие скорости перемещения ЗВ; а — коэффициент, определяющий изменение концентрации за счет процессов химического превращения примеси, а также ее осаждения (пропорционален К2). Дополняя это уравнение возмущающим воздействием х, получим
дУд + ^ ^ + vпдУ + „М —
дt 1 д{ п дп 5 35
——кМ + А К3»!. + А%дУ< - к* х - К* У. (17)
д С 5 дС дп п дп 35 5 35 1 2 * ;
Переход от модели (17) к модели вход-выход на основе функции Грина предложен в работе [3].
Предложенная в статье структура, математические модели и параметры ЗСУПТ и ее звеньев позволяют получить в результате моделирования качественные оценки устойчивости, регулирования, часть из которых приведена в работе. Также приведена оценка толерантности ЗСУПТ при изменении ее коэффициентов и параметров. Авторы использовали упрощенную модель преобразования массы ЗВ в составе факела в виде звена чистого запаздывания при принятых допущениях, один датчик, устанавливаемый в точке максимума концентрации ЗВ, и ПИД-регулятор. Это позволило по-
t, С
С
t, с
г)
t, с
д)
t, с
■ Рис. 4. Изменение ошибки 2 при ступенчатом изменении воздействия функции Х: а - при различных значения т1; б - при различных значениях К1; в - при различных значениях К9; г - при различных значениях К9 и т1; д - при различных значениях К7
лучить качественные характеристики процесса регулирования и перейти к экспериментам, уточнению значений параметров, а при необходимости — функций звеньев. Для реализации конкретных проектов ЗСУПТ необходима дальнейшая детализация математических моделей и моделирование с целью выбора и обоснования технических решений, в том числе выбора значений допусков параметров, обеспечения расположения измерителя в точке максимума концентрации, уточнения типа и параметров канала связи. Полученные результаты подтверждают концепцию построения ЗСУПТ для минимизации ЗВ и дают основания для опытно-конструкторской разработки этой системы.
Литература 4. Solnitsev R. I. Human factor minimization in the
1. Solnitsev R. I. The instrumentation in ecology and human safety // IEHS’96. ISA — SPb. Russian sect. / SPb. SUAI, 1996. P. 16-18. 2. Solnitsev R. I. The simulation of “Nature-technogenic” system // IEHS’98. ISA — SPb. Russian sect. / SPb. SUAI, 1998. P. 8-10. 3. Solnitsev R. I. Creation of “Nature-technogenic” control systems on the base of information technologies // IEHS’02. ISA — SPb. Russian sect. / SPb. SUAI, 2002. “Nature-technogenic” system // IEHS’04. ISA — SPb. Russian sect. / SPb. SUAI, 2004. P. 15-17. 5. Сольницев Р. И. Построение замкнутых систем «При-рода-техногеника»: Тр. ХХХШ МНТК IT+S&E’06, Украина, Крым, Ялта—Гурзуф, 20-30 мая 2006 г. // Открытое образование. С. 404-408. 6. Solnitsev R. I., Korshunov G. I., Klotchkov I. B. The “Nature-technogenic” closed system — innovational project // IEHS’07. ISA — SPb. Russian sect. / SPb. SUAI, 2007. P. 15-20.
P 1 9 1 7
/
ПАМЯТКА ДЛЯ АВТОРОВ
Поступающие в редакцию статьи проходят обязательное рецензирование.
При наличии положительной рецензии статья рассматривается редакционной коллегией. Принятая в печать статья направляется автору для согласования редакторских правок. После согласования автор представляет в редакцию окончательный вариант текста статьи.
Процедуры согласования текста статьи могут осуществляться как непосредственно в редакции, так и по е-таП (80x@mail.ru).
При отклонении статьи редакция представляет автору мотивированное заключение и рецензию, при необходимости доработать статью — рецензию. Рукописи не возвращаются.
Редакция журнала напоминает, что ответственность за достоверность и точность рекламных материалов несут рекламодатели.