-□ □-
Для питомог на одиницю бюмаси швидкостi окиснення шкдливог речо-вини отримано диференцшне рiвняння на макрорiвнi, яке описуе ктетику бю-хiмiчног деструктуризаци. На осно-вi феноменологiчного тдходу до уза-гальнення результатiв стащонарних лабораторних експериментiв запропо-нована утверсальна макроктетична математична модель при субстратному iнгiбiюваннi. Отримана модель пропонуеться як базова основа для кшьтсного опису нестащонарних про-цеЫв у бюреакторах
Ключовi слова: бiохiмiчна деструк-туризащя, питома швидтсть окиснення, феноменологiчний пiдхiд, суб-
стратне iнгiбiювання
□-□
Для удельной на единицу биомассы скорости окисления вредного вещества получено дифференциальное уравнение на макроуровне, описывающее кинетику биохимической деструк-туризации. На основании феноменологического подхода к обобщению результатов стационарных лабораторных экспериментов предложена универсальная макрокинетическая математическая модель при субстратном ингибировании. Полученная модель предлагается в качестве базовой основы для количественного описания нестационарных процессов в установках биоочистки
Ключевые слова: биохимическая деструктуризация, удельная скорость окисления, феноменологический
подход, субстратное ингибирование -□ □-
УДК 504.06+577.1
|DOI: 10.15587/1729-4061.2016.65036
РАЗРАБОТКА УНИВЕРСАЛЬНОЙ МОДЕЛИ КИНЕТИКИ СТАЦИОНАРНОГО ПРОЦЕССА БИООЧИСТКИ С СУБСТРАТНЫМ ИНГИБИРОВАНИЕМ
А. Ю. Бахарева
Кандидат технических наук, доцент* E-mail: baharevaann@gmail.com А. В. Шестопалов Кандидат технических наук, доцент** E-mail: shestopalov.it@khpi.edu.ua О. Н. Филенко Кандидат технических наук, доцент** E-mail: olesyafilenko200@mail.ru Б. Б. Кобылянский Кандидат технических наук, доцент Кафедра охраны труда и экологической безопасности Учебно-научный профессионально-педагогический институт Украинской инженерно-педагогической академии ул. Носакова, 9а, г. Бахмут, Украина, 84500 E-mail: b.kobiliansky@yandex.ua *Кафедра охраны труда и окружающей среды*** **Кафедра химической техники и промышленной экологии*** ***Национальный технический университет «Харьковский политехнический институт» ул. Кирпичева, 21, г. Харьков, Украина, 61002
1. Введение
Проблемой для многих коммунальных предприятий (мусорные свалки, мусороперерабатывающие заводы, канализационные сети, площадки активного ила) является загрязнение окружающей среды дур-нопахнущими продуктами разложения органических веществ, которыми является метан, сероводород и аммиак [1]. Низкая концентрация этих загрязняющих веществ от коммунального сектора делает традиционные физико-химические методы очистки (газофазное окисление, термические методы, абсорбция и т. д.) выбросов экономически не целесообразными. Высокая экологичность, экономичность и универсальность биотехнологических способов очистки промышленных загрязнений определяет их перспективность в сфере защиты окружающей среды. Проектирование процессов биохимической деструкции
вредных веществ требует разработки научно обоснованных методов расчета.
Процессы биохимической детоксикации, реализующиеся в биореакторах и установках периодического действия, являются нестационарными в силу внешнего воздействия. Внешнее воздействие на зону реакции выражается в виде поступления и удаления вредного вещества, а также в виде изменения ее объема. Таким образом, адекватное математическое описание нестационарного процесса будет представлять собой суперпозицию математической модели стационарного процесса и комбинации внешних воздействий, характерных для данной конкретной установки. Очевидно, что научно обоснованная методика расчета и проектирования процесса биоочистки в качестве основной составляющей должна содержать общую, универсальную для данного загрязнителя часть - математическое описание стационарного процесса.
©
Перспективность биотехнологических способов очистки промышленных загрязнений определяет актуальность разработки адекватных научно обоснованных методов проектных расчетов. Данное исследование посвящено разработке универсальной математической модели стационарного процесса, как базовой основы алгоритма расчета и проектирования установок биоочистки.
2. Анализ литературы и постановка проблемы
Сегодня именно биологические методы очистки газов во многих странах считаются наиболее рациональными и экономически обоснованными, так как процессы биохимического окисления и детоксика-ции загрязняющих веществ не требуют реагентов и протекают при температуре и давлении окружающей среды. Поэтому вопросам защиты окружающей среды с помощью биологических методов уделяется большое внимание во всем мире.
Сегодня нашли распространение следующие виды очистки газовых выбросов: биофильтры, биоскрубберы и реакторы с орошаемым слоем [2], которые показали достаточно высокую эффективность. Так, например, в публикации [3] показана возможность очистки от выбросов метана в присутствии неметановых органических веществ в орошаемом биофильтре и исследовано влияние толуола на макрокинетические и микрокинетические параметры процесса окисления метана. В статье [4] показана эффективность очистки выбросов метана от полигона захоронения ТБО, подтверждающая применимость биофильтрации для очистки выбросов от этого загрязнителя. В то же время данные статьи [5] показывают, что применение биологического метода очистки выбросов рентабельно лишь при низких концентрациях и небольшой производительности, так как при высоких концентрациях загрязнителя в воздухе требуются существенные затраты на обеспечение большей площади фильтрования и соответственно увеличение затрат на создание биофильтров.
В работе [6] исследуется возможность применения орошаемого биофильтра для деградации дисульфида углерода и доказана эффективность очистки от СS2 в газовом потоке со скоростью до 16,6 г на м3/час. Исследование [7] научно обосновывает применение орошаемых биофильтров для очистки выбросов подземных автомобильных парков от соединений толуола, эти-ленбензина и оксидов азота. Изучено также влияние накопления биомассы и пыли на гидравлическое сопротивление, влияющее на возрастание эксплуатационных затрат протирочных орошаемых биофильтров [8]. Известны и другие работы, описывающие кинетику и механизмы удаления сульфатов [9] и окислительные процессы в жидкой фазе очистных сооружений [10] с помощью живых организмов.
В литературе также отмечена высокая эффективность биоочистки выбросов от различных органических соединений (ацетон, толуол и бутанол), для которых при начальной концентрацией загрязняющего агента 100 мг/м3, эффективность биофильтра достигает 95 % [11]. Также достаточно хорошо рассмотрены
кинетические характеристики процесса очистки в биопленке сооружений очистки газов [12].
Ряд работ посвящен определению кинетики удаления сероводорода по результатам производственных данных [13] и разработке математических описаний процессов орошаемого биофильтра, основанной на статистическом методе оценки экспериментальных данных [14] или представлении о массопередаче в биофильтре [15].
Научной основой описания процессов биочистки служат закономерности биологической кинетики, являющиеся важнейшей частью исследований в биофизике [16]. Вместе с тем, применение на практике математических моделей биофизики не лишено определенных трудностей и недостатков.
Модели биологической кинетики математически описывают ферментативные реакции различных типов. На сегодняшний день рассмотрены десятки механизмов реакций, наиболее известная и простейшая из которых - модель Михаэлиса-Ментен, применимая для биологический очистки газов в орошаемых быо-фильтрах [17].
Математические модели представляют собой системы дифференциальных уравнений, упрощаемых в результате допущения про стационарность ферментативного акта, до систем обычных линейных уравнений. Решение данных систем всегда представляет собой дробную функцию, для практического использования которой необходимо определить кинетические коэффициенты, сопоставив расчетную формулу с лабораторными экспериментальными данными. Таким образом, математические модели ферментативной кинетики базируются не только на аналитическом представлении механизма реакции, но и на стационарных лабораторных исследованиях. При этом число экспериментально определяемых коэффициентов возрастает по мере усложнения предполагаемого механизма ферментативной реакции, как в случае учета явления субстратного ингибирования [18].
Структура моделей (дробная функция) не позволяет непосредственно графически или статистическим методом получить достоверные значения коэффициентов. Поэтому на практике применяется переход в системы с обратными координатами. При этом в случае субстратного ингибирования вводятся существенные допущения, сводящие исходную зависимость к формуле Моно для низких и к гиперболе для высоких концентраций. Прямой статистический расчет также неизбежно приводит к искажению фактических значений коэффициентов из-за применения системы с обратными координатами. Кроме того, при заранее неизвестном максимальном порядке непродуктивного фермент-субстратного комплекса, влияние уже четвертого комплекса может не улавливаться вследствие погрешности исходных данных. Поэтому учет в исследовательских расчетах сложного механизма субстратного ингибирования может нивелироваться неточностью количественной оценки параметров формул и погрешностью лабораторных экспериментов.
Изложенное свидетельствует о легитимности использования в практических расчетах подхода, в котором намеренное огрубление математического опи-
сания механизма реакции с использованием усредненных параметров дает возможность достоверного определения эмпирических коэффициентов. Такой подход называется феноменологическим и на нем основана химическая кинетика. В химической кинетике используется средняя концентрация, а собственно механизм реакции учитывается только ее порядком. При этом удельная скорость реакции определяется произведением кинетического коэффициента и концентрации, возведенной в степень. В рамках такого подхода возможно аналитически компактно учесть субстратное ингибирование на основе усредненной количественной оценки этого явления, как релаксации скорости реакции с увеличением концентрации удаляемого вещества.
3. Цель и задачи исследования
Цель исследования состоит в разработке адекватной математической модели на макроуровне, адекватно описывающей кинетику биохимической деструкции загрязнений окружающей среды с учетом явления ингибирования. Универсальность модели определяется структурно общими расчетными зависимостями для газообразных и растворенных в воде вредных веществ.
Для достижения поставленной цели необходимо решить следующие задачи:
- на основе феноменологического подхода к анализу и обобщению экспериментальных данных обосновать аналитическую форму макрокинетической модели биохимической деструкции и рассчитать численные параметры модели для различных видов загрязнений;
- доказать адекватность предлагаемой стационарной макрокинетической модели в широком классе биохимической деструкции с субстратным ингибиро-ванием.
4. Материалы и методы исследования кинетических характеристик биохимической деструкции
В колонной установке из оргстекла объемом 5 дм3 отрабатывали технологические параметры биохимического окисления растворенных в воде Н^, NH3 и SO2. Емкость заполняли водой с солями биогенных элементов и соединениями серы и азота, имитирующими продукты растворения в воде Н^ , NH3 и SO2.
Установка включала систему подачи воздуха, обеспечивающую мелкопузырчатую диспергацию. В установку помещали насадку (лавсановые ерши) с плотностью заполнения 150-200 мм/дм3 объема, на которой путем автоселекции формировали специализированный иммобилизированный микробиоценоз. Концентрацию S2-, SO32- и NH4+ в водной среде определяли по методикам, рекомендуемым в специальной литературе: S2- - по реакции образования ме-тиленовой сини, SO32- - титрованием раствора йода в присутствии формальдегида, NH4+ - с реактивом Несслера.
В колонной установке из оргстекла объемом 5 дм3 отрабатывали технологический режим окисления газообразного метана СН4. В установке реакционный объем был герметично изолирован от внешней среды. Загрузочный материал для иммобилизации метаноо-кисляющего микробиоценоза в установке представлял собой лавсановые ерши с плотностью 150-200 мм/дм3 объема. Биомассу метанотрофных бактерий получали на элективной среде. Инокулятором служили донные отложения из лабораторной установки, в которой моделировали лотковую часть канализационного коллектора.
Весовое соотношение распыленной в биоскруббере воды к прошедшему газу варьировали в диапазоне от 1:10 до 10:1.
Концентрацию метана в пробах газовоздушной смеси контролировали на газовом хроматографе «Цвет-250».
В установке объемом 2 дм3 устанавливали показатели биохимического окисления формальдегида в аэробных и анаэробных условиях. Эксперименты по биологической очистке водной среды от формальдегида проводили со свободноплавающей микрофлорой. При исследовании аэробного окисления формальдегида в установку через диспергирующее устройство подавали воздух. Условия диссимиля-ционной денитрификации - анаэробного окисления формальдегида - обеспечивали за счет снижения концентрации кислорода. Для этого при анаэробной очистке воды от формальдегида использовали механическое перемешивание иловой смеси, предотвращающее образование газообразных продуктов, нарушающих процесс биохимического окисления. Концентрацию формальдегида в модельной водной среде контролировали по показателю ХПК (поскольку другие органические соединения в среду не вводили). Проведенные эксперименты моделировали окисление формальдегида на этапе регенерации воды при очистке газообразных выбросов от этого загрязнителя в биоскруббере.
Оборудование и условия проведения экспериментов описаны в [19]. Моделировались стационарные процессы биохимической деструкции без дополнительной подачи загрязнителя в течение эксперимента (биореактор периодического действия). Первичные результаты эксперимента представляли собой точки в координатах «концентрация-время». Каждая точка являлась усреднением по трем-пяти значениям.
5. Результаты исследования кинетических характеристик биохимического деструктурирования
Зависимость «концентрация-время» не является инвариантом по отношению к начальной концентрации в эксперименте. Простота формы кривой допускает многовариантное статистическое описание. Однако, как будет показано ниже, при всей простоте формы, эта зависимость достаточно сложна для аналитического представления. Поэтому для анализа экспериментальных данных принимается удельная скорость окисления, определяемая как количество вредного
вещества, деструктурированного единицей биомассы в единицу времени, для любых пар значений концентрации и времени:
6. Обсуждение результатов исследования и математическое моделирование биохимического деструктурирования
Api w —^ = -ц„ V , At i'
(1)
где ^о - начальная концентрация биомассы; Ар^рм-^ - приращение концентрации; Дti=ti+1 - ^ - приращение времени; Vi - удельная скорость окисления.
Очевидно, что при этом значения Vi следует понимать как средние значения удельной скорости на интервале Арí и соотносить их со средними значениями концентрации на этом же интервале. Тогда аналитической формой описания стационарного процесса биохимической деструкции будет зависимость «удельная скорость - концентрация», инвариантная по отношению к начальной концентрации в эксперименте при постоянной концентрации биомассы.
Статистическая значимость предварительного описания кривых вида «удельная скорость-концентрация» полиномом второй степени доказала существование функциональной зависимости с максимумом [19]. Наличие максимума свидетельствует о субстратном ингибировании во всех рассматриваемых видах биохимической деструкции. Учет механизма субстратного ингибирования вызывает трудности и приводит к ошибкам при определении эмпирических коэффициентов. Феноменологические модели извлекаются из удачной параметризации экспериментальных данных, а не выводятся из теорий с объяснением механизма наблюдаемых явлений.
В нашем случае очевидны следующие явления:
1. Контакт молекулы загрязнителя с микроорганизмом, без которого невозможен биохимический процесс.
2. Угнетающее действие внешней среды на способность микроорганизма к биодеструкции, определяемое концентрацией вредного вещества.
Влияние первого феномена, как и в химической кинетике, может быть оценено степенной зависимостью от концентрации. Содержанием второго феномена является релаксация удельной скорости окисления с увеличением концентрации субстрата. Коэффициент, учитывающий его влияние, будет равен единице при нулевой концентрации загрязнителя, и стремиться к нулю при ее увеличении. В самом общем виде этому соответствует показательная кривая функция, реализующая эффект ингибирования при любой степенной зависимости от концентрации. Таким образом, добавив масштабный коэффициент, получим зависимость вида [20]:
Vp = a-рь
Полученная в предыдущем разделе зависимость (2) является необходимой, но не достаточной для практических расчетов стационарного процесса биохимической деструкции.
В предельном переходе из (1), при А^0, получим уравнение, описывающее процесс биохимической очистки в дифференциальной форме:
dp
dP = -^.
(3)
Уравнение (3) в дифференциальной форме на макроуровне описывает кинетику процесса биоокисления и его решением, при известной удельной скорости, будет недостающая до полного описания процесса зависимость концентрации загрязнения от времени. Таким образом, макрокинетическая математическая модель биохимической очистки представляет собой систему двух функций, количественно отражающих взаимозависимость концентрации удаляемого вещества, времени протекания процесса очистки, удельной скорости деструктурирования вредного вещества и начальной концентрации биомассы, полученных на основе взаимосвязи этих же параметров в дифференциальной форме:
p = fp(t); V = fv(p).
(4)
Результаты статистической обработки, представленные на рис. 1-4 и в табл. 1, свидетельствуют о том, что предложенный вид зависимости для V адекватно описывает экспериментальные данные. Полученные регрессионные формулы для всех исследованных процессов биохимической очистки статистически достоверны по F - критерию и коэффициенту детерминации. Таким образом, предложенная зависимость (2) является универсальной макрокинетической моделью биохимической деструкции.
(2)
45 40
35 30 25 20 15 10 5 0
40 80 120
р - концентрация H^S в воде, mi дм
160
где а, Ь, c - эмпирические коэффициенты.
При этом нахождение эмпирических коэффициентов не вызывает затруднений. Прологарифмировав выражение (2) и произведя замену переменных получим линейное уравнение множественной регрессии, определение коэффициентов которого возможно с помощью LINEST в EXCEL .
Рис. 1. Влияние концентрации Н^ в воде на скорость биосорбционного окисления: сплошная линия — уравнение регрессии вида Vр=aрbе-cр; пунктирная линия — численное
дифференцирование зависимости tp, полученной приближенным интегрированием уравнения регрессии вида Vр=aрbе-cр; точки — экспериментальные данные
45 40 35 30 25 20 15 10
- _
♦-
//
и
р
/
4
/
1
40 в0 120
р - концентрация ЙСЬ в воде, мг/дм3
16С
Необходимо подчеркнуть, что понятие макрокине-тической модели с необходимостью вытекает из определения важнейшего параметра - удельной скорости удаления загрязнения.
Решение уравнения (3) с использованием выражения (2) в замкнутой форме невозможно. Поэтому, для определенного интеграла
= ^--1} ^ 1 ^о ¿V
(5)
применим численное интегрирование, при котором на каждом 1-ом шаге точное значение скорости заменим ее линейным приближением:
Рис. 2. Влияние концентрации SO2 в воде на скорость биосорбционного окисления: сплошная линия — уравнение регрессии вида Vp=apbе-cp; пунктирная линия — численное
дифференцирование зависимости tp, полученной приближенным интегрированием уравнения регрессии вида Vp=apbе-cp; точки ■
■ экспериментальные данные
V = V + Мр, (6)
где к1 - значение производной V в 1-той точке. Проведя замену переменной
АН = г -^ (7)
V Ч
/
/ / /
'У* и
1
" 0 5 10 15 20 25
' р - концентрация ЫНз в воде, мг/дм3
г"*
Рис. 3. Влияние концентрации NHз в воде на скорость биосорбционного окисления: сплошная линия — уравнение регрессии вида Vp=apbе-cp; пунктирная линия — численное
дифференцирование зависимости tp, полученной приближенным интегрированием уравнения регрессии вида Vp=apbе-cp; точки — экспериментальные данные
300
I
и
3 250
а 200
3
о
¡2 - 150
и -Ь о г
и
• 100
50 О
- — - _
-3 / / /
// //
и
1
>
0 5 10 15 20
р - концентрация СН4 в воздухе оО %
Рис. 4. Влияние концентрации СН4 в воздухе на скорость биосорбционного окисления: сплошная линия — уравнение регрессии вида Vp=apbе-cp; пунктирная линия — численное
дифференцирование зависимости полученной приближенным интегрированием уравнения регрессии вида Vp=apbе-cp; точки — экспериментальные данные
и проинтегрировав (5), после преобразований получим формулы для расчета приращений концентрации на каждом ьом шаге численного интегрирования:
1
Ар1 =— (е-ф1 ^Аг -1),
ФР1
_Ь -
Р1
(8)
(9)
Точность и сходимость предложенного алгоритма оценивали по изменению величины концентрации в точке, близкой к концу расчета процесса очистки, так как ошибка может накапливаться по мере увеличения числа шагов интегрирования. На рис. 5 показано изменение величины концентрации в зависимости от шага интегрирования. Полученные данные свидетельствуют о том, что использованный метод является быстро сходящимся. Таким образом, при малых дг мы имеем практически точное в математическом смысле решение, отлично описывающее экспериментальные данные (рис. 6-9).
Изложенное позволяет считать доказанным, что предлагаемую эмпирическую зависимость V на статистически достоверном и достаточном для инженерных расчетов уровне точности можно использовать в качестве составной части макрокинетической математической модели биохимической очистки. В целом, найденные зависимости для удельных скоростей окисления в совокупности с методом численного интегрирования по сути представляют собой макрокинетические математические модели соответствующих биотехнологических процессов. В силу структурной общности зависимостей можно говорить об универсальной модели для различных видов загрязнителей.
Определенным неудобством данных моделей является отсутствие аналитических формул, позволяющих рассчитывать изменение концентрации в зависимости от продолжительности обработки.
Таблица 1
Результаты статистической обработки экспериментальных данных для зависимостей вида «удельная скорость — концентрация»
V - функция р - аргумент Коэффициенты уравнения регрессии вида Vр=aрьe-cр Доверительный интервал отклонений уравнения регрессии от эксперимента, SV Коэффициент детерминации, R2 Отношение расчетного значения F-критерия к табличному
а Ь с
V - удельная скорость окисления З-З2- , мг/г-ч 3,7689 0,6508 0,006172 1,94 0,9833 9,2
р - концентрация ЩЗ в воде, мг/дм3
V - удельная скорость окисления З-ЗОз2-, мг/г-ч 4,3363 0,6685 0,008102 2,45 0,9757 6,3
Р - концентрация 302 в воде, мг/дм3
V - удельная скорость окисления N-NH4+, мг/г-ч 1,1752 0,9318 0,07794 0,49 0,9329 2,2
Р - концентрация NHз в воде, мг/дм3
V - удельная скорость окисления З-содержащих субстратов, мг/г-ч 10,884 0,4899 0,003260 3,29 0,9853 10,5
Р - концентрация 3 содержащих субстратов в воде, мг/ дм3
V - удельная скорость окисления СН4, мл/г-ч 145,79 0,5485 0,06908 5,15 0,9986 74,3
р - концентрация СН4 в воздухе, об. %
V - удельная скорость аэробного окисления формальдегида в воде, мг/г-ч 0,02087 1,3850 0,001993 2,01 0,9457 1,9
Р - концентрация формальдегида в воде, мг/ дм3
V - удельная скорость анаэробного окисления формальдегида в воде, мг/г-ч 2,2563 0,8460 0,001027 0,08 0,99999 >100
р - концентрация формальдегида в воде, мг/ дм3
Рис. 5. Зависимость точности решения от дробности пошагового разбиения продолжительности обработки при численном интегрировании зависимости вида Vр=aрbе-cр
Существует возможность получения приближенного аналитического решения путем замены интегрального неравенства [23] на равенство. В нашем
случае - для начального условия р=р0 можно записать: р р
|р-Ь(1песр ^р |р-ьес|Мр
I р- >
к
ЯЬdP
(10)
Нетрудно увидеть, что для искомого интеграла, находящегося в числителе правой части интегрального неравенства, есть решение. После интегрирования и преобразований получили:
р(р0
t = т(рп -рп)е
1
^а(1 - Ь)
; п = 1 - Ь ; р = с
1 - Ь
2 - Ь .
(11)
(12)
Рис. 6. Влияние продолжительности обработки на концентрацию Н^ в воде: сплошная линия — численное интегрирование уравнения регрессии вида Vр=aрbе-cр; пунктирная линия — зависимость tp, полученная приближенным интегрированием уравнения регрессии вида Vр=aрbе-cр; точки — экспериментальные данные
Рис. 7. Влияние продолжительности обработки на концентрацию SO2 в воде: сплошная линия — численное интегрирование уравнения регрессии вида Vр=aрbе-cр;
пунктирная линия — зависимость tp, полученная приближенным интегрированием уравнения регрессии вида Vр=aрbе-cр; точки — экспериментальные данные
е
?ñ
n
i-,
ri
1» ?n
n
о
tt
CQ 15
5 В 10
cz
S 5
И
■
чч
v
-о 2 4 6 8
I - продолжительность обработки, часов
Рис. 8. Влияние продолжительности обработки на концентрацию NHз в воде: сплошная линия — численное интегрирование уравнения регрессии вида Vp=apbе-cp;
пунктирная линия — зависимость tp, полученная приближенным интегрированием уравнения регрессии вида Vp=apbе-cp; точки — экспериментальные данные
ю о а/
У
fe :
\
\
SO
t
20 40 60
- продолжительность обработки, мин.
Рис. 9. Влияние продолжительности обработки на концентрацию СН4 в воздухе: сплошная линия — численное интегрирование уравнения регрессии вида Vp=apbе-cp; пунктирная линия — зависимость полученная приближенным интегрированием уравнения регрессии вида р;точки — экспериментальные данные
Vp=apbe-cp"
Результаты расчетов по формуле (11) показали (рис. 6-9) высокую статистическую значимость, как по отношению к экспериментальным данным, так и по
соответствию результатам численного интегрирования. Минимальный коэффициент детерминации составлял Я2=0,9918. Методом численного дифференцирования была выполнена не обязательная для оценки практической применимости проверка соответствия формулы (11) дифференциальному уравнению (3). Даже в этом случае можно констатировать статистически достоверное соответствие (рис. 1-4).
Таким образом, показано, что уравнение регрессии вида (2) и взаимосвязанная с ним формула (11) в совокупности являются макрокинетической математической моделью процесса биохимической очистки, так как они адекватно описывают экспериментальные данные и получены на основе дифференциального соотношения (3). Полученная модель является универсальной базовой основой для разработки методик расчета конкретных нестационарных процессов в биореакторах.
7. Выводы
В результате проведенных исследований:
1. Экспериментально доказано наличие субстратного ингибирования для рассмотренных газообразных и растворенных в воде загрязнений окружающей среды. Исследовалась биохимическая деструкция газообразного метана, растворенных в воде сероводорода, аммиака, диоксида серы, а также аэробное и анаэробное биоокисление формальдегида.
2. На основе феноменологического подхода определена структура и рассчитаны численные параметры зависимости удельной (на единицу биомассы) скорости окисления от концентрации вредных веществ. Структура данной зависимости представляет собой произведение масштабного коэффициента, степенной и показательной зависимостей от концентрации.
3. Сформулировано понятие и получены расчетные формулы макрокинетической модели биохимической деструкции с субстратным ингибированием. Доказана адекватность модели для газообразного метана, растворенных в воде сероводорода, аммиака, диоксида серы, а также аэробного и анаэробного биоокисления формальдегида. Для данного класса процессов биоокисления разработанная модель является универсальной.
Литература
Shareefdeen, Z. Hydrogen Sulphide (H2S) Removal in Synthetic Media Biofilters [Text] / Z. Shareefdeen, B. Herner, D. Webb and S. Wilson // Environmental Progress. - 2003. - Vol. 22, Issue 3. - P. 207-213. doi: 10.1002/ep.670220319
Shareefdeen, Z. Biotechnology for Odor and Air Pollution Control [Text] / Z. Shareefdeen, A. Singh. - Springer, Berlin, 2005. doi: 10.1007/b138434
Ménard, C. Kinetics of simultaneous methane and toluene biofiltration in an inert packed bed [Text] / C. Ménard, A. A. Ramirez, M. Heitz // Journal of Chemical Technology and Biotechnology. - 2014. - Vol. 89, Issue 4. - P. 597-602. doi: 10.1002/jctb.4162 Park, S. Y. The Use of Biofilters to Reduce Atmospheric Methane Emissions from Landfills: Part 1, Biofilter Design [Text] / S. Y. Park, K. W. Brown and J. C. Thomas // Water, Soil and Air Pollution. - 2004. - Vol. 155, Issue 1-4. - P. 63-85. doi: 10.1023/b:wate.0000026522.36984.42
Nelson, M. Soil-Based Biofiltration for Air Purification:Potentials for Environmental and Space LifeSupport Application [Text] / M. Nelson, H. Bohn // Journal of Environmental Protection. - 2011. - Vol. 2, Issue 8. - P. 1084-1094. doi: 10.4236/jep.2011.28125 Rojo, N. Carbon disulfide biofiltration: Influence of the accumulation of biodegradation products on biomass development [Text] / N. Rojo, R. Muñoz, G. Gallastegui, A. Barona, L. Gurtubay, F. X. Prenafeta-Boldú, A. Elias // Journal of Chemical Technology and Biotechnology. - 2012. - Vol. 87, Issue 6. - P. 764-771. doi: 10.1002/jctb.3743
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
26^--
Rondeau, A. Biopurification of air containing a low concentration of TEX: comparison of removal efficiency using planted and non-planted biofilters [Text] / A. Rondeau, A. Mandon, L. Malhautier, F. Poly, A. Richaume // Journal of Chemical Technology & Biotechnology. - 2012. - Vol. 87, Issue 6. - P. 746-750. doi: 10.1002/jctb.3730
Andreasen, R. R. Pressure drop in biofilters as related to dust and biomass accumulation [Text] / R. R. Andreasen, R. E. Nicolai, T. G. Poulsen // Journal of Chemical Technology and Biotechnology. - 2012. - Vol. 87, Issue 6. - P. 806-816. doi: 10.1002/ jctb.4049
Papirio, S. Acid Mine Drainage Treatment in Fluidized-Bed Bioreactors by Sulfare-Reducing Bacteria: A Critical Review [Text] / S. Papirio, D. K. Villa-Gomez, G. Esposito, F. Pirozzi, P. N. L. Lens // Critical Reviews in Environmental Science and Technology. - 2013. - Vol. 43, Issue 23. - P. 2545-2580. doi: 10.1080/10643389.2012.694328
Oturan, M. A. Advanced Oxidation Processes in Water/Wastewater Treatment: Principles and Applications. A Review [Text] / M. A. Oturan, J.-J. Aaron // Critical Reviews in Environmental Science and Technology. - 2014. - Vol. 44, Issue 23. - P. 25772641. doi: 10.1080/10643389.2013.829765
Zagorskis, A. An Investigation on the Efficiency of Air Purification Using a Biofilter with Activated Bed of Different Origin [Text] / A. Zagorskis, R. Vaiskünaité // Chemical and Process Engineering. - 2014. - Vol. 35, Issue 4. - P. 435-445. doi: 10.2478/ cpe-2014-0033
González-Sánchez, A. Kinetic Characterization by Respirometry of Volatile Organic Compound-Degrading Biofilms from GasPhase Biological Filters [Text] / A. González-Sánchez, L. Arellano-García, W. Bonilla-Blancas, G. Baquerizo, S. Hernández, D. Gabriel, S. Revah // Industrial & Engineering Chemistry Research. - 2014. - Vol. 53, Issue 50. - P. 19405-19415. doi: 10.1021/ie503327f
Shareefdeen, Z. Hydrogen Sulphide Removal Using a Novel Biofilter Media [Text] / Z. Shareefdeen, A. Aidan, W. Ahmed, M. B. Khatri, M. Islam, R. Lecheheb. F. Shams // World Academy of Science, Engineering and Technology. - 2010. - Vol. 62. -P. 13-16.
Shareefdeen, Z. Kinetics and Modeling of H2S Removal in a Novel Biofilter [Text] / Z. Shareefdeen, A. Aidan, W. Ahmed // Advances in Chemical Engineering and Science. - 2011. - Vol. 1, Issue 2. - P. 72-76. doi: 10.4236/aces.2011.12012 Bonilla-Blancas, W. Application of a novel respirometric methodology to characterize mass transfer and activity of H2S-oxidizing biofilms in biotrickling filter beds [Text] / W. Bonilla-Blancas, M. Mora, S. Revah, J. A. Baeza, J. Lafuente, X. Gamisans et. al. // Biochemical Engineering Journal. - 2015. - Vol. 99. - P. 24-34 doi: 10.1016/j.bej.2015.02.030
Романовский, Ю. М. Математическое моделирование в биофизике [Текст] / Ю.М. Романовский, Н. В.Степанова, Д. С. Чернавский. - Москва Ижевск: Институт компьютерных исследований, 2003. - 402 с.
Veillette, M. Biofiltration of air polluted with methane at concentration levels similar to swine slurry emissions: influence of ammonium concentration [Text] / M. Veillette, A. A. Ramirez, M. Heitz // Journal of Environmental Science and Health, Part A. - 2012. - Vol. 47, Issue 7. - P. 1053-1064. doi: 10.1080/10934529.2012.667327
Surerus, V. Activated sludge inhibition capacity index [Text] / V. Surerus, G. Giordano, L. A. C. Teixeira // Brazilian Journal of Chemical Engineering. - 2014. - Vol. 31, Issue 2. - P. 385-392. doi: 10.1590/0104-6632.20140312s00002516 Кричковська, Л. В. Проектш ршення у розробщ апара^в бюлопчно! очистки газоподiбних викидiв [Текст]: монографiя / Л. В. Кричковська, Л. А. Васьковець, I. В. Гуренко та ш.; за ред. Л. В. Кричковсько!. - Харгав: НТУ «ХП1», 2014. - 208 с. Бахарева, Г. Ю. Розробка макрокшетично! моделi процесу бюлопчно! очистки газопов^ряних сумшей [Текст] / Г. Ю. Бахарева, О. В. Шестопалов, 6. О. Семенов, Н. О. Букатенко // ScienceRise. - 2015. - Т. 2, № 2 (7) - С. 12-15. doi: 10.15587/2313-8416.2015.37057