УДК 681.54
с1о1:10.15217/1ззп1684-8853.2015.4.59
СИНТЕЗ СИСТЕМЫ УПРАВЛЕНИЯ ТЕХНОЛОГИЧЕСКИМ ПРОЦЕССОМ УДАЛЕНИЯ СОЕДИНЕНИЙ АЗОТА ИЗ СТОЧНЫХ ВОД
Е. К. Грудяева3'1, аспирант
аСанкт-Петербургский государственный электротехнический университет «ЛЭТИ», Санкт-Петербург, РФ
Введение: существующие системы очистки сточных вод активным илом используются для узкого диапазона концентраций загрязнений. Применение регулирования позволяет расширить диапазон, однако при изменении состава воды на входе системы или условий внешней среды требуется постоянная перенастройка параметров регулятора. Цель исследования: синтез системы управления с логико-динамическим регулятором для удаления загрязнителей в широком диапазоне концентраций. Результаты: на основе математического описания модели АЭМ1 построена схема взаимодействия компонентов биохимических реакций. Представлена управляемая математическая модель удаления соединений азота в двухзонном биореакторе с рециклом нитратов, блоком ультрафильтрационных мембран и возвратным иловым потоком в первую зону. Введены показатели качества процессов для системы управления биологической очисткой стоков. На их основе разработана методика параметрической настройки регулятора с логически перестраиваемой структурой. Введены понятия ведущего и ложного ведущего регуляторов. В соответствии с разработанной моделью в среде Ма^аЬ/ЭтиНпк построена компьютерная модель. Результаты моделирования демонстрируют возможность качественной очистки в условиях широкого изменения уровней входных многокомпонентных загрязнений в заданном неизменном объеме биореактора. Практическая значимость: применение логико-динамического регулятора с оптимальными настройками параметров позволяет обеспечивать допустимые концентрации аммония, нитратов и химическое потребление кислорода на выходе на всех рабочих режимах системы управления.
Ключевые слова — нитрификация, денитрификация, активный ил, двухзонный биореактор, рециркуляционный поток, математическое моделирование, АЭМ1, логико-динамический регулятор.
Введение
Потребность в разработке математических моделей (ММ) обусловлена требованиями высокого качества управления технологическим процессом (ТП), ограниченным набором средств измерения, а также временными затратами при использовании лабораторных анализов и сложностью качественного анализа проб иловой массы. Процессы биологической очистки сточных вод активным илом описываются ММ различной степени сложности и адекватности.
Одной из ММ, представляющих основные процессы в активном иле, является модель ASM1 [1, 2]. Эта нелинейная модель характеризует процессы удаления соединений азота в результате жизнедеятельности групп бактерий. Бактерии авто-трофы обеспечивают процесс нитрификации при наличии кислорода, в результате которого аммонийный азот окисляется до нитратного. Процесс денитрификации обусловлен действием бактерий гетеротрофов, которые в бескислородных условиях удаляют нитратный азот, преобразуя его в газообразный. В дальнейшем модель ASM1 берется за основу моделирования.
1 Научный руководитель — доктор технических наук, профессор кафедры автоматики и процессов управления Санкт-Петербургского государственного электротехнического университета «ЛЭТИ» С. Е. Душин.
Математические модели
Рассматриваемый биореактор (рис. 1) состоит из зон денитрификации Д и нитрификации Н, причем, из Н в Д предусмотрен рециркуляционный поток. Возвратный поток из мембранного блока М поступает в зону денитрификации, а фильтрат — на последующее обеззараживание и далее — потребителю (обозначения потоков приводятся ниже в таблице).
Размер пор в мембранах ультрафильтрации блока М составляет 0,03-0,1 мкм. Взвешенные частицы имеют размер более 0,1 мкм, размер бактерий нитрификаторов — от 0,8 мкм, дени-трификаторов — от 1 мкм. Таким образом, принимается, что ультрафильтрационные мембраны фильтруют все взвешенные вещества и бактерии.
Модель должна учитывать химическое потребление кислорода (ХПК) — количество кислорода, потребляемое при химическом окислении содержащихся в воде органических веществ
ев
■ Рис. 1. Схема биореактора
под действием различных сильных окислителен в присутствии катализатора. Показатель ХПК характеризует необходимые затраты кислорода на удаление загрязняющих веществ из сточной воды [2] и рассчитывается как сумма следующих фракций [3]: растворенное легко разлагаемое органическое вещество; растворенное биологически инертное органическое вещество; взвешенное медленно биологически разлагаемое органическое вещество; взвешенное биологически инертное разлагаемое органическое вещество; биомассы бактерий нитрификаторов и денитрификаторов. С учетом принятого допущения о фильтрации мембранами всех взвешенных веществ и бактерий, ХПК после мембранного блока складывается из растворенных веществ — легко разлагаемого органического и биологически инертного.
При разработке ММ принимались следующие ограничения и допущения:
— зоны биореактора характеризуются идеальным перемешиванием; пространственная распределенность потоков не учитывается;
— состав активного ила образован бактериями нитрификаторов и денитрификаторов;
— температура в технологическом пространстве аппарата фиксирована и равна 20 °С;
— на вход биореактора поступают аммонийный азот, нитратный азот, растворенное биологически
разлагаемое органическое вещество, а также взвешенное органически медленно разлагаемое вещество; входные стоки также характеризуются щелочностью;
— в модели мембранного блока процессы нитрификации и денитрификации не учитываются.
Процессы в зонах нитрификации и денитри-фикации характеризуются однотипными уравнениями химической кинетики. Подача кислорода в зону нитрификации позволяет учесть преобладание нитрифицирующих процессов, тогда как в бескислородной зоне эффективнее протекают процессы денитрификации. Схема, дающая представление о взаимодействии компонентов модели АЯМ1, представлена на рис. 2.
Для удобства восприятия составляющие модели, их обозначения и размерности сведены в таблицу.
С учетом обозначений процессы в зонах Д и Н описываются системами нелинейных обыкновенных дифференциальных уравнений.
Модель зоны Д (денитрификация):
ах,
Ьа1
0вх + 0Р + 0в ' V
' V
х,
Ьа2
0
+ ~ ХЬа3 ~ У1
ХЬа1 + гЬа1 (^пк1, ^о1, ХЬа1
■ Рис. 2. Схема взаимодействия компонентов модели АБМ1
■ Таблица
№ Наименование величин и параметров Обозначение Размерность
1 Концентрация бактерий-нитрификаторов ХЬа г ХПК/м3
2 Концентрация бактерий-денитрификаторов ХЬк г ХПК/м3
3 Концентрация аммония и аммонийного азота ^пк г N4+/^
4 Концентрация нитратного азота г Шд/м3
5 Концентрация растворенного биологически разлагаемого органического вещества Я* г ХПК/м3
6 Концентрация взвешенного органического медленно разлагаемого вещества X г ХПК/м3
7 Концентрация растворенного кислорода Яо г 02/м3
8 Концентрация растворенного органического азота г ^м3
9 Концентрация взвешенного органического медленно разлагаемого азота Хпй г ^м3
10 Щелочность Я щел экв/м3
11 Концентрация твердых продуктов распада биомассы г ХПК/м3
12 Концентрация растворимого инертного органического вещества «I г ХПК/м3
13 Концентрация взвешенного инертного органического вещества г ХПК/м3
14 Химическое потребление кислорода ^ХПК г ХПК/м3
15 Удельная объемная скорость роста концентрации вещества ], где ] — составляющие модели ASM1 1 — 11, [...] — соответствующая единица измерения массы г1 г [...]/м3
16 Расход сточных вод входного потока биореактора евх м3/сут
17 Расход рециркуляционного потока нитратов ер м3/сут
18 Расход возвратного потока из мембранного блока ев м3/сут
19 Объем зоны 1 биореактора м3
АХШ1 _ <р ^
& V
ХЬк2 + ~ХЬк3 -
<вх + <р + <в
VI
~ХЬЫ + гЬМ (по1' ^31' ХЬЫ );
& V,
«пй + <~«пк2 + V!
+ 0вх+0р+0вс + «пПЗ ~ ьиП1 +
У1 У1
" гпЫ (пП1' &по1' «31' ®о1' ' ХЬа1' ХЬЫ );
_ <вх явх + <р ч +
Яв явх + Яр + Яв
+—«поЗ ~ «по1 +
+ гпо1 («по1' «31' «о1' «пй1' ХЬа1' ХЬЫ ); 31
¿«31 <вх «вх + <р ч +
-йз +--«з2 +
& V
V
+яв _ явх + яр + яв _ + +11йз3 V йз1 +
+ гзз1 (йпо1' йз1' йо1' ХЬЫ' Хз1 );
¿Хз1 _ <вх ^вх . <р
& V
~Х„ +--Хз2 +
V
+ <в Х <вх + <р + <в Х +
+--Х„я--Х„1 +
v1 33 V 31
+ г.
.31 ((о1' йо1' ХЬа1' ХЬЫ' Х31 );
<Р О" <вх + <р + <в
3о1 __
& V
йо2 + «оЗ -
'оЗ
«о1 +
" го1 (пк1' «31' So1' ХЬа1' ХЬП1 );
р Яв явх + Яр + Яв
1. «па2 +~«паз ~ «па1 +
+ г3
3п41
(по1' «о1' ХЬа1' ХЬЫ' Хпа1' Х31 );
аХпа1 <р „
р Яв явх + Яр + Ов
тт Хпа2 +—Хпаз —-Хпа1 +
У1 У1 У1
+ г.
хпд,1
(«по1' «о1' ХЬа1' ХЬЫ' Хпа1' Х31 );
dSщeл1 _ 0 „вх <
свх _о
т. тт '-'щел т тт к-'щел2
¥1
+<Р1в
я™ + <р +
ЩелЗ '
V,
- Ч + щел1
" гщел1 ^пт' «пЛ' «31' «о1' 8па1' ХЬа1' ХЬН1 );
ахР1 - о? аг V овх + ор + ов
ов
Хр2 + ХР3 -
V
р3
V
хр1 + Гр1 (ХЬа1, ХЬк1 );
СХПК1 - Хз1 + Х1 + Хз1 + х1 + ХхЬа1 + ХхЬк1 ■ (1) Модель зоны Н (нитрификация):
~(ХЬа1 - ХЬа2 ) +
аХЬа2 - овх + ор + ов
а* V
+ гЬа2 (Хпк2, Хо2, ХЬа2 );
АХьн2 - 0вх + 0р + 0в
(ХЬк1 - ХЬк2 ) +
аг
+ гЬк2 (Хпо2, Хз2, Хо2, ХЬк2 );
АБПк2 - 0вх + 0р + 0в
(Хпк1 - Хпк2 ) +
& V2
+ гпк2 (Хпк2, Хпо2, Хз2, Хо2, Хпй2, ХЬа2, ХЬк2 );
(Хпо1 - Хпо2 ) +
аБПо2 - овх + ор + ов а V
+ гпо2 (Хпо2, Хз2, Хо2, Хп^2, ХЬа2, ХЬк2 );
ая32 - овх + ор + ов
а* у2 _(Хз1-Хз2)+
+ гзз2 (Хпо2, Хз2, Хо2, ХЬк2, Хз2 );
ах32 -овх + ор + ов (х х )+ ^Г - V (з1" Хз2)+
+ гхз2 (Хпо2, Хо2, ХЬа2, ХЬк2, Хз2 );
Хо
%2 )
аХо2 - овх + ор + ов аг
+ го2 (Хпк2, Хз2' Хо2, ХЬа2, ХЬк2 );
аХп^2 - овх + ор + ов (Х Х )+
ТТ - 71 - Хпа2) +
аг К2
+ гзпй2 (Хпо2, Хо2, ХЬа2, ХЬк2, Хпй2, Хз2 );
¿Хп^2 0вх + 0р + 0в
±пЛ2
аг
у2
(Хп41 - Хпй2 ) +
+ гхпй2 (Хпо2, Хо2, ХЬа2, ХЬк2, Хпй2, Хз2 ); аХщел2 - 0вх + 0р + 0в
аг
у2
(Хщел1 Хщел2)"
+ гщел2 (Хпк2, Хпо2, Хз2, Хо2, Хпй2, ХЬа2, ХЬк2 );
охр! - ((р1 - хр2)+
аг у2 1 р1 р2)
+ гр2 (ХЬа2, ХЬк2 ); СХПК2 - Хз2 + Х1 + Хз2 + Х1 + ХхЬа2 + ХхЬк2 ■ (2)
В мембранном блоке происходит разделение потока на возвратный и выходной. Модель возвратного потока включает все компоненты потока иловой смеси из биореактора. В потоке фильтрата присутствуют только растворенные компоненты смеси. Таким образом, концентрация нерастворимых веществ в возвратном потоке увеличивается обратно пропорционально расходу.
Модель возвратного потока:
ХЬа3°В - ХЬа2°вх; ХЬк3°В - ХЬк2°В*;
Хпк3 - Хпк2' Хпо3 - Хпо2; Хз3 - Хз2;
Хз3<°В - Хз2°вх; Хо3 - Хо2; Хп43 - Хп42;
Хп43°В - Хп^2°вх; Хщел3 - Хщел2;
ЯВ хр3 - хр2°вх;
СХПК3 - Хз3 + Х1 + хз3 + х1 + ХхЬа3 + ХхЬк3 ■ (3)
Модель выходного потока фильтрата: ХЬа4 - 0; ХЬк4 - 0; Хпк4 - Хпк2;
Хпо4 - Хпо2; Хз4 - Хз2; Хз4 - 0; Хо4 - Хо2;
Хп44 - Хп42; Хп44 - 0;
Хщел4 - Хщел2; Хр4 = 0; СХПК4 - Хз2 + Х1 ■ (4)
В биореакторе зоны нитрификации и дени-трификации принимаются одинаковыми. Соотношение входного и возвратного потоков определяется из технических характеристик ультрафильтрационных мембранных блоков.
Постановка задачи
Для каждой составляющей загрязнения при технологическом расчете системы очистки используются фиксированные средние показатели, по которым вычисляют конструктивные параметры реактора и значение поступающего расхода. Анализ литературных источников [4-6] показал, что для ТП многозонного реактора с рециклом и подачей кислорода принятые за рубежом нормы очистки сточных вод на порядок хуже, чем требуется по нормам, принятым в Российской Федерации [7].
Поскольку состав сточных вод на выходе реактора является многокомпонентным и концентрации компонентов меняются в различных диапазонах, выделяются три составляющие субстрата: аммоний, нитраты и быстро биологически разлагаемое вещество. На основе данных, полученных с 28 очистных сооружений в России и странах ближнего зарубежья, были найдены минималь-
ные, средние и максимальные по уровню значения загрязнений по аммонию, нитратам и легко биологически разлагаемому веществу, определяющие режимы работы системы. Объем всех зон реактора принимался одинаковым и равным 100 м3. Как показали исследования, при подаче на вход реактора даже «средних» загрязнений получаемая на выходе концентрация аммония оказывалась выше установленных норм, что требует обеспечения эффективного его удаления из системы.
Ставится задача разработки ММ управляемых ТП очистки сточной воды, удовлетворяющих требованиям СНиП и показателям качества процессов управления на основе использования логико-динамического регулятора. Кроме того, необходимо, чтобы производительность установки очистки была максимальна при удовлетворении этих условий.
Показатели качества процессов управления
Установившиеся и переходные процессы в динамических системах биоочистки, представленные нелинейными моделями (1)-(4), характеризуются определенными показателями качества управления.
Одними из основных показателей качества, которые непосредственно характеризуют степень очистки сточных вод при использовании ММ АЯМ1, являются установившиеся ошибки по концентрациям аммония в^к, нитратов в^о и ХПК вХПк:
„уст _ ХЗ
сnh
_ Slh - Sl™ > 0; elf _ - SyCT > 0;
VCT _ x3 - хуст > 0
еХПК _ ХХПК хХПК > 0'
(5)
где ХЗк, ХЗо, Х^пк — заданные значения концентраций (уставки) по аммонию, нитратам и ХПК; ХукСт, Х;УСТ, ХХСТк — установившиеся значения концентраций, подлежащие измерениям.
На выходные концентрации накладываются ограничения в виде Хпк < Х^™, Хпо < Хп°п,
ХХПК < ХХПК , где ЯДЛ Х;?°П, ХХПК — допустимые по СНиП установившиеся значения.
При использовании ПИ (ПИД)-регуляторов по ошибкам уставки принимают значения
Хпк = ХпкП, Х1по = ХпоП, ХХПК = ХХХ<ГХК ■ В случае
использования П-регуляторов по ошибкам для уставок можно принять Хпк < ХДОП Х\о < ХДОП
сЗ ^ оДОП
ХХПК <ХХПК■
К важному динамическому показателю очистки относят максимальное положительное отклонение концентрации от установившегося значения. При регулировании процессов это значение
не должно превышать 5 % от допустимого значения:
S„TX < 1,05. Sär; smax < 1,05. sr;
SXnK <1,05 ■ sXIIK • (6)
Другим важным динамическим показателем очистки, характеризующим скорость ее протекания, является время установления Рст переходного процесса, который требуется минимизировать:
tyCT _ min. (7)
Установившееся значение расхода субстрата на входе реактора Q^ целесообразно максимизировать. Это требование обусловлено необходимостью переработки установкой определенного количества сточных вод, т. е. производительность установки должна быть адекватна требованиям нагрузки, что обеспечивает очистку подаваемых стоков в течение продолжительного периода. Теми же соображениями обусловлено требование максимизации интегрального показателя входного расхода, характеризующего суммарный объем стоков, прошедших через систему за период очистки
VBX _ jQBXdi _ max.
(8)
0
В случае неоднозначности выбора настроек регулирования вводятся в рассмотрение также интегральные показатели загрязнений, характеризующие суммарное покомпонентное количество загрязнений на выходе системы за период очистки íu, которые подлежат минимизации в результате регулирования:
*и *и
ХпкЕ _ | Хпк2аг _ шт; ХпоЕ _ { Хпо2°Г _ т*п;
0
0
t„
SXnKI _ j SXnK4dt _ min.
(9)
0
Приведенные показатели качества являются наиболее характерными для рассматриваемых технологических процессов.
Методика параметрической настройки
Методика отражает порядок настройки логико-динамического регулятора системы управления биореактора с выносным мембранным блоком фильтрации. Регулятор входного расхода (рис. 3) состоит из трех ПИД-субрегуляторов Р. Его работа основана на выборе минимального управляющего воздействия.
■ Рис. 3. Схема регулятора
Методика настройки ПИД-субрегуляторов осуществляется в соответствии с показателями качества, приведенными ранее, и сводится к следующему.
1. На вход системы подается субстрат со «средними» концентрациями аммония, нитратов и легко биологически разлагаемого вещества среднего уровня. Коэффициенты интегральной и дифференциальной составляющих субрегуляторов полагаются равными нулю. Отслеживаются значения аммония, нитратов, ХПК на выходах системы, управляющих воздействий каждого субрегулятора, а также результирующего управляющего воздействия регулятора. По графикам процессов изменения концентраций определяется выполнение условий (5).
Если условия (5) нарушаются, то по графикам управляющих воздействий субрегуляторов выявляется регулятор, выход которого совпадает с результирующим воздействием, т. е. воздействие которого является минимальным. Этот регулятор считается ведущим. Воздействие ведущего регулятора определяется формулой
<£х = ш1п(«?х), I = 1, 2, 3,
где Я?х — воздействие 1-го регулятора, причем I = 1 для регулятора концентрации аммония, I = 2 — концентрации нитратов и I = 3 — концентрации ХПК.
2. Для ведущего регулятора производится увеличение пропорционального коэффициента и сравнение показателей качества процессов в порядке приоритетности с показателями, полученными в п. 1.
3. Если в соответствии с (5) значения установившихся ошибок приблизились к нулевым значениям, настройки п. 2 считаются лучше, чем в п. 1.
4. Если значения ошибок не изменились, то берется следующий по приоритетности показатель качества, и в соответствии с (6) определяются лучшие настройки. В случае невозможности определить лучший вариант, сравнение продолжается по каждому следующему показателю качества с (7) по (9).
5. Если увеличение коэффициента приводит к ухудшению показателей качества процессов, то пропорциональный коэффициент следует уменьшать относительно значения, полученного в п. 1. Далее сравнение настроек проводится так же, как в пп. 4 и 5.
6. Изменение пропорционального коэффициента проводится до тех пор, пока не обеспечиваются условия (5) либо по условию (5) лучшие результаты не достигаются, а в соответствии с условиями (6) или (7) изменения коэффициента не принесут лучших результатов.
7. Подключается интегральная составляющая ведущего регулятора, и показатели качества получившихся процессов сравниваются с предыдущими (лучшими) в п. 7.
8. Проводится изменение интегрального коэффициента в целях улучшения качества процессов до тех пор, пока не выполнятся условия (5), а в соответствии с условиями (6) или (7) изменения коэффициента не принесут лучших результатов; в противном случае интегральная составляющая отключается.
9. Подключается дифференциальная составляющая, и настройка проводится аналогично пп. 8 и 9.
10. В процессе изменения коэффициентов (параметров) ведущий регулятор может поменяться. В этом случае продолжается настройка по указанному алгоритму уже для другого ведущего регулятора.
11. При настройке необходимо отслеживать возникновение ложного ведущего регулятора. Управляющий сигнал ложного ведущего регулятора совпадает с результирующим воздействием, однако концентрация загрязнителя, на регулирование которой он направлен, гораздо ниже допустимого значения, т. е. установившаяся ошибка отрицательна, и превышения допустимого значения нет. При этом концентрация другого загрязнителя превосходит допустимое значение.
Управляющее воздействие ложного ведущего регулятора определяется условиями
явх (*) = ш1п(явх (*)), еУст < 0, I * ].
Ложный ведущий регулятор может быть настроен таким образом, что все показатели качества по регулируемой им концентрации удовлетворяются. В данном случае необходимо проводить настройку ложного ведущего регулятора таким образом, чтобы его управляющий сигнал превысил воздействие регулятора концентрации, не удовлетворяющей показателям качества. Как правило, для этого следует в несколько раз увеличить пропорциональную составляющую ложного ведущего регулятора.
После того как регулятор настроен на очистку «средних» загрязнений, аналогично производится настройка для максимальных и минимальных по уровням загрязнений (соответствующих режимов) на входе. Поскольку выделены три компонента загрязнения, то всего возможны 27 динамических ситуаций. Для каждой ситуации находятся свои оптимальные настройки субрегуляторов и свои ведущие регуляторы.
При переходе к настройкам каждой следующей динамической ситуации используются настройки той ситуации, значения концентраций загрязнений на входе которой близки к текущей, а ведущий регулятор совпадает.
Полученная серия настроек регуляторов запоминается и используется в зависимости от того, в каком из режимов (с минимальным, средним или максимальным уровнем) находится каждый из трех субстратов, причем переключение субрегуляторов происходит автоматически.
Параметрический синтез системы управления
Ниже приводятся результаты исследования поведения системы биоочистки при различных уровнях загрязнений на входе.
Процесс изменения концентраций загрязнений на выходе мембранного блока при постоянном расходе на входе биореактора без регулирования отражен на рис. 4. Концентрации аммония и легко биоразлагаемого субстратов на входе максимальны, концентрация нитратов минимальна. Штриховой линией обозначена 5%-я зона допустимых концентраций. Как видно из графиков, концентрации аммония и ХПК во много раз превышают допустимые уровни.
Результаты моделирования управляемых процессов удаления загрязнений из сточных вод с помощью регулятора входного расхода при наилучших настройках параметров приводятся на рис. 5. Даже в этом случае на графике концентрации аммония наблюдается отклонение от установившегося значения, что можно объяснить ограниченной скоростью роста популяции бактерий-нитрификаторов.
Графики концентраций на выходе системы без регулирования при подаче максимальной концентрации нитратов, концентрации аммония и легко биоразлагаемого субстратов на входе представлены на рис. 6. Наблюдается превышение концентрации нитратов в установившемся режиме, а также превышение концентрации по аммонию в переходном режиме.
Реакцию системы управления с использованием логико-динамического регулятора демонстрируют графики на рис. 7. Значения расхода на входе стремятся к нулю при заданных входных концентрациях. Увеличить установившееся значение расхода в этом случае, а также снизить максимальное положительное отклонение концентрации аммония (см. рис. 5) возможно за счет введения дополнительных управляющих координат по расходам потока рецикла и возвратного потока, а также по концентрации растворенного кислорода.
100
50
пк
Qв
100 200 300 400
t, сут
100 200 300 400 500
t, сут
500
0,1
0,05
600 400 200 0
0
СХПК
100 200 300
t, сут
100 200 300
t, сут
500
500
■ Рис. 4. Процессы концентраций загрязнителей и расхода в системе без регулирования
в
по
0
0
0
0
0
400
400
МОДЕЛИРОВАНИЕ СИСТЕМ И ПРОЦЕССОВ
1
0,5 0 1
0,5 0
'пк
Бп
* V--
<Эв
0 100 200 300 400 500
■)вх *, Сут
0 100 200 300 400 500
*, сут
0,2 0,1 0
40 30 20 10
0 100 200 300 400 500
п сут
СХПК
0 100 200 300 400 500
*, сут
■ Рис. 5. Процессы концентраций загрязнителей и расхода в системе управления
5
пк
0,5
0
Qвх
&ГГ"
100 200 300 400
*, сут
100 200 300 400
*, сут
500
500
60 40 20 0
40 30 20 10
0
С
ХПК
0
100 200 300 400 500
*, сут
100 200 300 400 500
*, сут
■ Рис. 6. Процессы концентраций загрязнителей и расхода в системе без регулирования при максимальных входных уровнях
по
1
0
0,2 0,1 0 20 10
пк
4 6
*, сут
0 20 40 60 80 100 120
евх *, СУт
10
10
С
40 30 20 10
ХПК
*, сут
0 20 40 60 80 100 120
0 20 40 60 80 100 120
*, сут
■ Рис. 7. Процессы концентраций загрязнителей и расхода в системе управления при максимальных входных уровнях
по
5
0
8
Для увеличения расхода на входе и для улучшения качества регулируемых процессов следует также применить регулирование расхода рецикла и концентрации кислорода [8].
Существует ограничение по регулированию ХПК концентрацией растворимого вещества, значения которой определяются в результате химического анализа.
Заключение
Логико-динамическое регулирование обеспечивает необходимое качество воды незави-
Литература
1. Henze M., Grady C. P. L., Gujer W., Marais G. v. R., Matsuo T. Activated Sludge Model No. 1: IAWPRC Scientific and Technical Report. — London: IAWPRC, 1987. N 1. — 37 p.
2. Жмур Н. С. Технологические и биохимические процессы очистки сточных вод на сооружениях с аэро-тенками. — М.: АКВАРОС, 2003. — 512 с.
3. Хенце М., Армоэс П., Ля-Кур-Янсен Й., Арван Э. Очистка сточных вод. — М.: Мир, 2009. — 480 с.
4. Гордеева Ю. Л., Гордеев Л. С. Математическая модель непрерывного процесса в биореакторе с рециклом субстрата и биомассы // Вестник АГТУ. Сер. Управление, вычислительная техника и информатика. 2013. № 2. С. 9-18.
5. Weijers S. R. Modelling, Identification and Control of Activated Sludge Plants for Nitrogen Removal. —
симо от концентраций загрязнений исходных стоков в неизменном объеме биореактора без дополнительных подстроек параметров регулирования в процессе работы системы и без внесения конструктивных изменений установки очистки стоков. Рассмотренная методика синтеза позволяет унифицировать разработку систем очистки и сократить время и средства при проектировании биореактора. Возможно снижение производительности установки в случае высоких концентраций загрязнителей, однако потери качества очистки воды на выходе не происходит.
Eindhoven: Technische Universiteit Eindhoven, 2000. — 235 p.
6. Stahl T., Duffy G., Kestel S., Gray M. Dissolved Oxygen Control Based in Real-Time Oxygen Uptake Rate Estimation // Florida Water Resources Journal. 2013. N 4. P. 50-53.
7. СНиП 2.04.03-85. Канализация. Наружные сети и сооружения. — М.: ОАО «ЦПП», 2012. — 97 с.
8. Грудяева Е. К., Душин С. Е. Моделирование управляемых технологических процессов нитрификации и денитрификации в биологической очистке сточных вод// Проблемы автоматизации и управления в технических системах: сб. ст. Междунар. науч.-техн. конф., Пенза, 19-21 мая 2015 г./ под ред. М. А. Щербакова. Пенза: Изд-во ПГУ, 2015. Т. 1. С. 351-355.
UDC 681.54
doi:10.15217/issn1684-8853.2015.4.59
Synthesizing a Control System for Removing Nitrogen Compounds from Wastewater
Grudyaeva E. K.a, Post-Graduate Student, [email protected]
aSaint-Petersburg State Electrotechnical University «LETI», 5, Prof. Popov St., 197376, Saint-Petersburg, Russian Federation
Introduction: The existing systems of sewage treatment by activated sludge are used for a narrow range of pollution concentrations. Control can expand the range; however, when the water composition changes at the system entrance or when the environmental conditions vary, the control parameters have to be continuously readjusted. Purpose: The goal of this work is synthesizing a control system with dynamic logic to remove contaminants in a wide range of concentrations. Results: On the base of ASM1 mathematical model description, a scheme of interaction between biochemical reaction components is developed. A controllable mathematical model is presented for removing nitrogen compounds in a two-zone bioreactor with recycling of nitrates, a block of ultrafiltration membranes and a returnable sludge stream towards the first zone. Process quality indices are introduced for the biological wastewater treatment control system. On their basis, a technique of parametrical setup of the regulator with logically reconstructed structure is developed. For the regulators, the concepts of a leader and a false leader are introduced. According to the developed model, a computer model is constructed in MatLab/Simulink environment. The modeling results demonstrate the possibility of high-quality cleaning when the levels of entrance multicomponent pollutions change in a wide range in a predefined invariable volume of the bioreactor. Practical relevance: Using the logic-dynamic regulator with optimum settings of parameters helps to provide admissible output concentration of ammonium, nitrates and COD at all the operating modes of the control system.
Keywords — Nitrification, Denitrification, Activated Sludge, Two-Zone Bioreactor, Recirculation Stream, Mathematical Modeling, ASM1, Logic-Dynamic Regulator.
References
1. Henze M., Grady C. P. L., Gujer W., Marais G. v. R., Matsuo T. Activated Sludge Model No. 1. IAWPRC Scientific and Technical Report. London: IAWPRC, 1987, no. 1. 37 p.
2. Jmour N. S. Tekhnologicheskie i biokhimicheskie protsessy ochistki stochnykh vod na sooruzheniiakh s aeroten-kami [Technological and Biochemical Processes of Waste Water Treatment on Treatment Plants with Aero-tanks]. Moscow, AKUAROS Publ., 2003. 512 p. (In Russian).
3. Henze M., Harremoes P., La Cour Jansen J., Arvin E. Wastewater Treatment. Biological and Chemical Processes. Berlin, Springer-Verlag, 1997. 383 p.
4. Gordeeva Y. L., Gordeev L. S. Mathematical Model for Continuous Process in Bioreactor with Substrate and Biomass Recycling. Vestnik AGTU, Ser. Upravlenie, vychislitel'naia tekhnika i informatika, 2013, no. 2, pp. 9-18 (In Russian).
5. Weijers S. R. Modelling, Identification and Control of Activated Sludge Plants for Nitrogen Removal. Eindhoven, Technische Universiteit Eindhoven, 2000. 235 p.
6. Stahl T., Duffy G., Kestel S., Gray M. Dissolved Oxygen Control Based in Real-Time Oxygen Uptake Rate Estimation. Florida Water Resources Journal, 2013, no. 4, pp. 50-53.
7. Construction Norms and Regulations 2.04.03-85. Sewerage. External Networks and Constructions. Moskow, OAO "TsPP" Publ., 2012. 97 p. (In Russian).
8. Grudyaeva E. K., Dushin S. E. Modelling of Controlled Technological Processes of Nitrification and Denitrification in a Biological Wastewater Treatment. Sbornik statei Mezhdun-arodnoi nauchno-tekhnicheskoi konferentsii "Problemy avto-matizatsii i upravleniia v tekhnicheskikh sistemakh" [Proc. Int. Conf. "Problems of Automation and Control in Technical Systems"]. Penza, 2015, vol. 1, pp. 351-355 (In Russian).
Научный журнал «ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ» выходит каждые два месяца.
Стоимость годовой подписки (6 номеров) для подписчиков России — 4800 рублей,
для подписчиков стран СНГ — 5400 рублей, включая НДС 18%, таможенные и почтовые расходы.
Подписку на печатную версию журнала можно оформить в любом отделении связи по каталогу:
«Роспечать»: № 48060 — годовой индекс, № 15385 — полугодовой индекс,
а также через посредство подписных агентств:
«Северо-Западное агентство „Прессинформ"»
Санкт-Петербург, тел.: (812) 335-97-51, 337-23-05,
эл. почта: [email protected], [email protected],
сайт: http://www.pinform.spb.ru
«МК-Периодика» (РФ + 90 стран)
Москва, тел.: (495) 681-91-37, 681-87-47,
эл. почта: [email protected], сайт: http://www.periodicals.ru «Информнаука» (РФ + ближнее и дальнее зарубежье) Москва, тел.: (495) 787-38-73, эл. почта: [email protected], сайт: http://www.informnauka.com «Деловая пресса»
Москва, тел.: (495) 962-11-11, эл. почта: [email protected], сайт: http://delpress.ru/contacts.html «Коммерсант-Курьер»
Казань, тел.: (843) 291-09-99, 291-09-47, эл. почта: [email protected],
сайт: http://www.komcur.ru/contacts/kazan/
«Урал-Пресс» (филиалы в 40 городах РФ)
Сайт: http://www.ural-press.ru
«Идея» (Украина)
Сайт: http://idea.com.ua
«ВТЪ» (Узбекистан)
Сайт: http://btl.sk.uz/ru/cat17.html и др.
На электронную версию нашего журнала (все выпуски, годовая подписка, один выпуск, одна статья)
вы можете подписаться на сайтах НЭБ: http://elibrary.ru;
РУКОНТ: http://www.rucont.ru; ИВИС: http://www.ivis.ru/
Полнотекстовые версии журнала за 2002-2014 гг.
в свободном доступе на сайте журнала (http://www.i-us.ru),
НЭБ (http://www.elibrary.ru)
и Киберленинки (http://cyberleninka.ru/
journal/n/informatsionno-upravlyayuschiesistemy).
Печатную версию архивных выпусков журнала за 2003-2014 гг. вы можете заказать в редакции по льготной цене.