Динамические модели управляемых биохимических процессов очистки сточных вод Текст научной статьи по специальности «Математика»

УДК 681.54

DOI: 10.17586/0021-3454-2015-58-9-732-737

ДИНАМИЧЕСКИЕ МОДЕЛИ УПРАВЛЯЕМЫХ БИОХИМИЧЕСКИХ ПРОЦЕССОВ ОЧИСТКИ СТОЧНЫХ ВОД

Е. К. Грудяева, С. Е. Душин, Н. Н. Кузьмин

Санкт-Петербургский государственный электротехнический университет „ЛЭТИ" им. В. И. Ульянова (Ленина), 197022, Санкт-Петербург, Россия E-mail: [email protected]

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

Ключевые слова: очистка сточных вод, трехзонный биореактор, активный ил, нитрификация, денитрификация, регулирование.

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

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

Типовая технологическая схема очистных сооружений с предвключенной денитрифи-кацией приведена в работе [3]. В ней предусмотрен также рецикл из зоны нитрификации (Н) в зону денитрификации (Д). Недостаток схемы состоит в том, что кислород из зоны Н поступает в зону Д, препятствуя эффективному удалению нитратов. В схеме, представленной на рис. 1, учитывается этот недостаток, поскольку рецикл направлен из последующей зоны (3), в которой кислород отсутствует.

X

^,вх ^,вх ^,вх Snh , ^ , Sno , Q

вх

Д Зона 1

Н Зона 2

Д

Зона 3

o

Рис. 1

На рисунке приняты следующие обозначения*: 0вх — расход сточных вод на входе биореактора; 0р — расход рециркуляционного потока; — концентрация растворенного кислорода на входе в реактор; , Бпо/ — концентрации соответственно аммонийного и нитратного азота в /-й зоне реактора, /=1, 2, 3; — концентрация легко биологически разлагаемого вещества.

Анализ моделей процессов нитрификации и денитрификации. Как показали исследования [4], математические модели процессов нитрификации и денитрификации характеризуются тремя состояниями равновесия, одно из которых физически неосуществимо. В зависимости от значений входных концентраций субстратов и кислорода, а также времени нахождения в зонах реактора модели переходят в состояние равновесия: точку вымывания ила из реактора или точку роста популяции (у каждой модели — собственная особая точка). При этом точка роста может быть как устойчивым узлом, так и устойчивым фокусом. Также в работе [4] для модели денитрификации определены наиболее эффективные, относительно очистки, интервалы времени реакции и значения входных концентраций субстратов.

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

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

Постановка задачи и моделирование управляемых процессов. На сегодняшний день для каждой субстанции загрязнения при технологическом расчете системы очистки используются фиксированные средние концентрации загрязнений. На основе этих данных рассчитываются конструктивные параметры реактора и значение входного расхода. Проведенные исследования [5—7] для управляемых моделей многозонного реактора с контуром рециркуляции и подачей кислорода показали результаты на порядок хуже, чем требуется по российским нормам [8].

Моделирование технологического процесса проводилось при входных значениях концентраций загрязнений, усредненных по 28 объектам очистки. Были определены средние и максимальные значения концентраций загрязнений по аммонию, нитратам и легко биологически разлагаемому веществу. Объем всех зон принимался одинаковым и равным 100 м3. Расходы жидкости на входе реактора и в контуре рециркуляции фиксированы. Поскольку при подаче на вход реактора „средних" загрязнений получаемая на выходе концентрация аммония выше установленных норм, требуется обеспечить эффективное его удаление из системы.

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

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

В настоящей статье используются обозначения, принятые в работах [1, 2].

управляет входной концентрацией кислорода с учетом ошибки концентрации аммония. Зависимость концентрации аммония от расхода кислорода не учитывалась.

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

— зоны биореактора характеризуются идеальным перемешиванием иловой массы;

— состав активного ила образован нитрифицирующими и денитрифицирующими бактериями;

— температура в аппарате фиксирована и равна 20 °С;

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

Модель управляемого трехзонного биореактора с подачей кислорода во вторую зону представлена системой из трех групп уравнений — по 11 уравнений в каждой группе (всего 33 уравнения) [1, 2]:

СХЪаЦ) _ ир х 0вх + ир

^ =Р~Хъа3 -

(ХЪа(г) - ХЪа(г-1) ) + ГЪа (Бпк(г), Бо(г), ХЪа(г) ) ;

СХ,

Ък(г) ир —— _р—

С V

СБ,

Сг

ХЪк3 - 0 т+ (ХЪк(г) - ХЪк(г-1)) + гЪк (Бпо(г), Б5(г), Бо(г), ХЪк(г)) ;

евх+ир (Б Б )+

т/ ^пк(г) Бпк(г-1) ^+

г 01 ч v

41

V

Бвх + ^_б

Бпк + тг Бпк3

+ Гпк (Бпк(г), Бпо(г), (г), Бо(г), БпС(г), ХЪа(г), ХЪк(г) ) ;

сб,

по(г)

Сг

_р

0

вх

и р

овх _ о

-гт по ' -г г ^по3 V У

евх+ир (Б -Б )+

V ^по(г) Бпо(г-1) ^+

+ Гпо (Бпо(г), Б5(г), Бо(г), БпС(г), ХЪа(г), ХЪк(г) ) ;

сб

5 (г)

Сг

_р

0

вх

и р

-бвх +—б

V 5 ' V "53

V г г У

0вх + и р

V (Б5(г) - Б5(г-1) ) + Г55 (Бпо(г), Б5(г), Бо(г), ХЪк(г), Х5(г) ) ;

сбо(/1 _р0рБ 0вх + ир(Б -Б )г,+ ири

С "Р р. Бо3 ^^ ) Бо(г+ V у и

СХ 5(г) _ ир Х 0вх + ир

-_ Р-Х 5 3--

сСг V 5

+Го (Бпк(г), Б5(г), Бо(г), ХЪа(г), ХЪк(г) );

V ((5(г) - Х5(г -1) ) + ГХ5 (Бпо(г), Бо(г), ХЪа(г), ХЪк(г), Х5(г) ) ;

СБ

пС (г) и р

-— _Р-

Сг V

СХ,

пС (г) и р

-— _Р-

Сг V

0 вх + и р

БпС 3 V ((пС (г) - БпС (г-1) ) + Г5пС (Бпо(г), Бо(г), БпС (г), ХЪк(г), ХпС (г), Х5 (г) ) ;

'г

0 вх + и р

ХпС 3 V (бпС (г) - ХпС (г -1) ) + ГхпС (Бпо(г), Бо(г), ХЪа (г), ХЪк(г), ХпС (г), Х5(г) ) ;

^ат(,-) _ ир _ евх + иР

~1Г= 3 V""

(Ба1к(г) - Ба1к(г-1) ) +

ЖХ Ж

+Га1к (Бпк(1), Бпо(г), Б.?(г), Бо(г), (г), Х6а(г), ХЪН(1) ) ;

_ Р ир ХР3 - ^ т+-(Хр(г) - Хр(г-1) ) + Глр (ХЬа(1), ХЬЙ(г) ) , г _ 1, 2, 3 :

где

Р_

[1, г _ 1; 0, г _ 2,3,

V _

[1, г _ 2;

0, г _ 1,3,

ир _

ср, б > ср; бр, 0 < бр < ср; ио Ч 0, бр < 0,

Со (0Х > Со;

С, 0 < С < Со;

0, С < 0,

при этом если индекс какой-либо переменной концентрации Хг или принимает нулевое значение (г _ 0 ), то соответствующая переменная равна нулю ( Х0 _ 0, $0 _ 0 ); р — коэффициент, учитывающий концентрации элементов рециркуляционного потока и их концентрации во входном потоке; V — коэффициент, учитывающий концентрации растворенного кислорода, подаваемого в зону 2 (нитрификации); г — скорости изменения концентраций в результате химических реакций, представляющие собой нелинейные функции концентраций; ир — сигнал управления расходом при рецикле; ио — сигнал управления концентрацией кислорода; Ср, Со — ограничения сверху для управляющих сигналов ир и ио; 0 — температура иловой смеси в реакторе; V — объем г -й зоны; концентрации аммонийного и нитратного Бпо азота являются регулируемыми на выходе переменными.

На рис. 2 приведены графики переходных процессов концентраций аммонийного и нитратного азота на выходах двух- и трехзонного реакторов (кривые 1 и 2 соответственно) равного объема при одинаковых постоянных воздействиях на входе. Как видно из графиков, трехзонный реактор при одном и том же расходе входного потока обеспечивает лучшую очистку от соединений азота. По нитратному азоту отклонение переходного процесса от установившегося значения меньше и время установления также меньше.

£„/,. г/м 10,95

10,90

10,85 10,80

1 //

1

Бпо, г/м 1,5

0,5

1

Г ~~ И 1 ч 2

:---- — - -

п 1

0

400

600

800 Т, сутки

20

40

60

Т, сутки

Рис. 2

На рис. 3 представлены графики концентраций аммонийного и нитратного азота на выходе трехзонного реактора с регулируемым рециклом из зоны нитрификации и регулируемой концентрацией растворенного кислорода в этой зоне. Переходные процессы получены при постоянном уровне подачи „средних" загрязнений (кривая 1) и залповой подаче в течение 5 дней максимальных загрязнений в состоянии равновесия на 200-е сутки (кривая 2), а также при залпе максимальных загрязнений в начале процесса (кривая 3).

1

0

Рис. 3

Заключение. Разработанные динамические математические модели управляемых процессов нитрификации и денитрификации для трехзонного биореактора могут быть положены в основу построения установок биологической очистки сточных вод. Проведенные исследования с использованием предложенных моделей показали, что регулирование расхода рециркуляционного потока и концентрации растворенного кислорода в зоне нитрификации при фиксированном входном расходе позволяет обеспечить допустимые концентрации загрязнений на выходе при „средних" концентрациях на входе. Система регулирования способна справляться с залповыми нагрузками, подаваемыми в различные моменты времени. Трехзон-ный реактор обеспечивает лучшую очистку, чем двухзонный за счет снижения негативного влияния кислорода при рецикле на денитрификацию. Таким образом, за счет предлагаемой схемы регулирования система может обеспечивать необходимое качество очистки при сравнительно широком диапазоне входных концентраций загрязнений. Это позволяет проектировать реакторы меньшего объема и, следовательно, снизить стоимость очистной установки.

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

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

1. Henze M., Grady C. P. L., Gujer W., Marais G. v. R., Matsuo T. Activated Sludge Model No. 1: Scientific and Technical Report. London: IAWPRC, 1987. 37 p.

2. Хенце М., Армоэс П., Ля-Кур-Янсен Й., Арван Э. Очистка сточных вод. М.: Мир, 2009. 480 с.

3. Большаков Н. Ю. Биотехнологии на страже экосистемы // СтройПРОФИ. 2013. № 2(11). С. 54—56.

4. Грудяева Е. К., Душин С. Е. Анализ динамики процесса нитрификации в модели биологической очистки сточных вод ASM1 // Изв. СПбГЭТУ „ЛЭТИ". 2015. Вып. 1. С. 12—17.

5. Гордеева Ю. Л., Гордеев Л. С. Математическая модель непрерывного процесса в биореакторе с рециклом субстрата и биомассы // Вестник АГТУ. Сер.: Управление, вычислительная техника и информатика. 2013. № 2. С. 9—18.

6. Weijers S. R. Modelling, Identification and Control of Activated Sludge Plants for Nitrogen Removal. Eindhoven: Techn. Univ. Eindhoven, 2000. 235 p.

7. 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.

8. СНиП 2.04.03-85. Канализация. Наружные сети и сооружения. М.: ОАО „ЦПП", 2012. 97 с.

Сведения об авторах

Елизавета Камаловна Грудяева — аспирант; СПбГЭТУ „ЛЭТИ", кафедра автоматики и процессов управления; E-mail: [email protected] Сергей Евгеньевич Душин — д-р техн. наук, профессор; СПбГЭТУ „ЛЭТИ", кафедра автоматики и

процессов управления; E-mail: [email protected]

Николай Николаевич Кузьмин — канд. техн. наук, профессор; СПбГЭТУ „ЛЭТИ", кафедра автоматики и

процессов управления; E-mail: [email protected]

Рекомендована кафедрой Поступила в редакцию

автоматики и процессов управления 22.04.15 г.

Ссылка для цитирования: Грудяева Е. К., Душин С. Е., Кузьмин Н. Н. Динамические модели управляемых биохимических процессов очистки сточных вод // Изв. вузов. Приборостроение. 2015. Т. 58, № 9. С. 732—737.

DYNAMIC MODELS OF THE CONTROLLED BIOCHEMICAL PROCESSES OF WASTEWATER TREATMENT

E. K. Grudyaeva, S. E. Dushin, N. N. Kuzmin

Saint Petersburg Electrotechnical University "LETI", 197022, Saint Petersburg, Russia

E-mail: [email protected]

Nonlinear dynamic models of nitrification and denitrification processes in the bioreactor with the corresponding zones and the controlled recirculation contour are presented. The models are based on ASM1 model proposed by Henze to describe the microbiological processes in active sludge that involve removal of nitrogen compounds from wastewater, change of alkalinity, concentration of soluble substratum and inert compounds, as well as bacteria populations. Graphs of transient processes obtained with computer modeling are demonstrated.

Keywords: wastewater treatment, three-zone bioreactor, active sludge, nitrification, denitrification,

control.

Data on authors

Elizaveta K. Grudyaeva — Post-Graduate Student; ETU "LETI", Department of Automation and

Control Processes; E-mail: [email protected] Sergey E. Dushin — Dr. Sci., Professor; ETU "LETI", Department of Automation and Control

Processes; E-mail: [email protected] Nikolay N. Kuzmin — PhD, Professor; ETU "LETI", Department of Automation and Control

Processes; E-mail: [email protected]

For citation: Grudyaeva E. K., Dushin S. E., Kuzmin N. N. Dynamic models of the controlled biochemical processes of wastewater treatment // Izvestiya Vysshikh Uchebnykh Zavedeniy. Priborostroenie. 2015. Vol. 58, N 9. P. 732—737 (in Russian).

DOI: 10.17586/0021-3454-2015-58-9-732-737