CC BY
73
УДК 531/534: [57+61]
Российский Журнал
www.biomech.ru
ПРИМЕНЕНИЕ ПАКЕТА ПРИКЛАДНЫХ ПРОГРАММ ДЛЯ РЕШЕНИЯ ЗАДАЧ, СВЯЗАННЫХ С ОЧИСТКОЙ ПОЧВЫ ОТ НЕФТИ И ТЯЖЕЛЫХ МЕТАЛЛОВ ПОСРЕДСТВОМ БИОСУРФАКТАНТА
О.А. Коростина1, М.С. Куюкина2,3, И.Б. Ившина2,3, Л^. Костина2,
М.А. Осипенко1, Ю.И. Няшин1
1 Кафедра теоретической механики Пермского национального исследовательского политехнического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: oma@theormech.pstu.ac.ru Институт экологии и генетики микроорганизмов Уральского отделения Российской академии наук, Россия, 614081, Пермь, ул. Голева, 13, e-mail: kuyukina@iegm.ru
3 Кафедра микробиологии и иммунологии Пермского государственного национального исследовательского университета, Россия, 614900, Пермь, ул. Букирева, 15, e-mail: ivshina@iegm.ru
Аннотация. Разработан пакет прикладных программ, позволяющий решать следующие задачи биомеханики: вычисление времени отмывания биогенным поверхностно-активным веществом (биосурфактантом) загрязненной нефтью почвы, находящейся в накопительном резервуаре, вычисление коэффициента сопротивления при фильтрации биосурфактанта в почве, загрязненной солями тяжелых металлов, прогнозирование эффективности биологической очистки (биоремедиации) почвы. Описываются основные блоки, составляющие пакет прикладных программ, а также кратко формулируются биомеханические и математические модели, на основе которых производятся расчеты в данном пакете.
Ключевые слова: биомеханика, пакет прикладных программ, загрязнение почвы, биоремедиация, нефтезагрязненная почва, тяжелые металлы, сорбция, десорбция, ЯЛсСососсиз-биосурфактант, математическое моделирование.
Введение
В последние годы с целью ликвидации последствий углеводородных загрязнений почвы все чаще применяют эффективные и экологически безопасные технологии биологической очистки (биоремедиации), основанные на использовании микроорганизмов и продуктов их жизнедеятельности, в частности биогенных поверхностно-активных веществ (биосурфактантов).
В лаборатории алканотрофных микроорганизмов Института экологии и генетики микроорганизмов УрО РАН получены высокоактивные биосурфактанты, синтезируемые бактериями рода Rhodococcus. Полученные биосурфактанты с высокой нефтеэмульгирующей и нефтеотмывающей способностью успешно применяются в процессах биоремедиации загрязненных почвогрунтов [1]. Лабораторные
© Коростина О.А, Куюкина М.С., Ившина И.Б., Костина Л.В., Осипенко М.А., Няшин Ю.И., 2011 Коростина Ольга Андреевна, соискатель кафедры теоретической механики, Пермь
Куюкина Мария Станиславовна, д.б.н., в.н.с. лаборатории алканотрофных микроорганизмов, профессор кафедры микробиологии и иммунологии, Пермь
Ившина Ирина Борисовна, д.б.н., профессор, чл.-корр. РАН, зав. лабораторией алканотрофных микроорганизмов, профессор кафедры микробиологии и иммунологии, Пермь
Костина Людмила Викторовна, к.б.н., инженер лаборатории алканотрофных микроорганизмов, Пермь Осипенко Михаил Анатольевич, к.ф.-м.н., доцент кафедры теоретической механики, Пермь Няшин Юрий Иванович, д.т.н., профессор, зав. кафедрой теоретической механики, Пермь
09806267
исследования показали, что Rhodococcus-биосурфактант способствует эффективной десорбции нефтяных углеводородов [6, 7] и солей тяжелых металлов [2, 3, 8, 9] из загрязненной почвы.
Для подбора эффективных приемов и создания оптимальной схемы биоремедиации необходимо построение биомеханических и математических моделей, адекватно описывающих процесс очистки загрязненного почвогрунта от нефти и тяжелых металлов с помощью биосурфактанта. Ранее авторами построены согласующиеся с экспериментальными данными биомеханические модели: модель проникновения нефти в почве, модель проникновения биосурфактанта в нефтезагрязненной почве, модель проникновения биосурфактанта в почве, загрязненной солями тяжелых металлов [2, 4].
Прогнозное математическое моделирование - эффективный инструмент современной экологической биотехнологии, направленный на адекватное описание и прогнозирование природных процессов в условиях техногенного воздействия. В частности, построение математических моделей позволяет находить оптимальное решение, адекватное описанию процесса загрязнения почв нефтью и нефтепродуктами, прогнозирование последствий нарушения почвенных процессов и выбор оптимальной стратегии биоремедиации [4]. Следует отметить, что очищение почвы от загрязняющих веществ (углеводородов нефти) - процесс комплексный, в котором участвуют физические, химические и биологические факторы. При этом действие всех факторов взаимообусловлено и синергично. Исследование сложных механизмов биоремедиации требует создания пакета прикладных программ на основе биомеханических и математических моделей, описывающих процессы очистки почвы от нефти и солей тяжелых металлов под действием биосурфактантов.
В настоящей работе представлен пакет прикладных программ, разработанный на основе биомеханических и математических моделей процесса фильтрации Rhodococcus-биосурфактанта в нефтезагрязненной почве и почве, загрязненной солями тяжелых металлов.
Описание пакета прикладных программ
Пакет прикладных программ «Очистка загрязненной нефтью и тяжелыми металлами почвы с помощью биосурфактанта» (рис. 1) позволяет решать следующие биомеханические задачи: вычисление времени отмывания биосурфактантом
загрязненной нефтью почвы, находящейся в накопительном резервуаре; вычисление коэффициента сопротивления к фильтрации биосурфактанта в почве, загрязненной солями тяжелых металлов, прогнозирование эффективности биоремедиации. Пакет состоит из трех блоков, функционирование которых описано ниже.
Блок «Расчет времени отмывания биосурфактантом загрязненной нефтью почвы, находящейся в накопительном резервуаре» предназначен для определения времени прохождения заданного количества биосурфактанта в загрязненной нефтью почве, находящейся в накопительном резервуаре, а также для определения средней скорости вытекания биосурфактанта из накопительного резервуара.
Загрязненный нефтепродуктами почвогрунт находится в резервуаре, имеющем форму прямоугольного параллелепипеда (рис. 2). На верхней границе (открытая поверхность) давление атмосферное; на дне находится вакуумный насос, который откачивает биосурфактант под давлением, равным атмосферному; боковые стенки и часть дна (10 < х < 11) - непроницаемые поверхности.
Открытая поверхность орошается биосурфактантом. В нефтезагрязненном почвогрунте происходит фильтрация биосурфактанта, осуществляющаяся в два этапа.
Рис. 1. Пакет прикладных программ «Очистка загрязненной нефтью и тяжелыми металлами почвы
с помощью биосурфактанта»
На первом этапе имеется фронт проникновения биосурфактанта. На втором этапе, после достижения биосурфактантом дна резервуара, происходит его откачивание через трубу с помощью вакуумного насоса.
Процесс фильтрации описывается уравнениями безынерционного течения идеальной несжимаемой жидкости под действием силы тяжести, объемной силы сопротивления, пропорциональной скорости, и силы сопротивления, сосредоточенной на фронте проникновения и зависящей от скорости. Для решения использованы математические модели: задача Коши для обыкновенного дифференциального
уравнения и краевая задача для уравнения Лапласа в прямоугольнике. А именно, на первом этапе фильтрации
gh0/а + ( g - kdl(ї)/dt) І(ї) - / (dl (ї)/dt) = 0,
где g - ускорение свободного падения, ^ - толщина слоя биосурфактанта на открытой поверхности почвогрунта, а - пористость почвогрунта, k - коэффициент сопротивления при фильтрации, І - глубина проникновения биосурфактанта, /(•) -зависящая от скорости сила сопротивления, сосредоточенная на фронте проникновения биосурфактанта. На втором этапе фильтрации
Ар(х, у, z) = 0,
где А - оператор Лапласа, р - давление биосурфактанта. Краевые условия для давления на границе резервуара
Ф/йх|х=0 = ф/&Ц = ф/ф|у=0 = ф/Ф|у=4 = ^ ф/&|2=0 = ^ при у > l0,
где І1 - ширина резервуара, І2 - длина резервуара, І0 - ширина трубы, р -эффективная плотность биосурфактанта. Так как на открытой поверхности резервуара
Рис. 2. Модель резервуара с загрязненным почвогрунтом
Рис. 3. Схематическое изображение почвенной колонки
имеется слой биосурфактанта толщиной , а откачивание биосурфактанта производится вакуумным насосом, то
р 2=/з = Р0 + Р 0 ^0, Р\2= 0 = 0 При У < 10 ’ где 13 - высота резервуара, р0 - атмосферное давление, р0 - плотность чистого
биосурфактанта. Последнее условие заменяется на приближенное
Ф/дАг=0 = Р(Ш _ ё) При У < 10 ,
где и следует найти из соотношения
1(0' С р ^ у , 0) йхйУ = °.
Величина и есть средняя (по области 0 < х < I', 0 < у < 10) скорость откачивания
биосурфактанта. Решение задачи для р имеет вид [5]
Р( У,2) = Р о +Р 0 gh0 +Р . ъп10
/
(/3 - 2) -
2 У
2рШ2
sm
я
п=1
/2
где
и =
2 и ЯПІз
п ch-----------
з, япМзлА соз ЯПУ
£+-1
k Ы 3
а
У
/о + 212
12 Я 31 о/з
2 яп/0 , яп/3
, ^---------------------3
‘п=1~3 / 2 / 2
Ею '
,^11
12 Я~ 1013 П П3 12 12
Блок «Определение параметра k фильтрации биосурфактанта в почве, загрязненной солями тяжелых металлов» предназначен для определения к - коэффициента сопротивления со стороны среды при фильтрации биосурфактанта в почве на основе экспериментальных данных по фильтрации биосурфактанта в почве, загрязненной солями тяжелых металлов.
/
/
2
2
Эксперимент проводится в модельных почвенных колонках. Почва, загрязненная солями тяжелых металлов (никеля и молибдена), находится в стеклянной колонке. Одна колонка схематично показана на рис. 3. Здесь h0 - высота столба налитой сверху жидкости (водного раствора RhodOcoccus-биосурфактанта), l - глубина проникновения жидкости в почву, t - время; поверхность почвы отвечает координате x = 0. Параметры колонок - диаметр d и высота L слоя почвы (рис. 3).
Процесс фильтрации биосурфактанта в модельной почвенной колонке описывается уравнениями безынерционного течения идеальной несжимаемой жидкости под действием силы тяжести и объемной силы сопротивления, пропорциональной скорости. Высота h столба жидкости как функция времени оказывается согласно теории различной для солей никеля и молибдена, а именно для солей никеля
h(t) = al - a21,
а для солей молибдена
h(t) = a3 exp(-a4t) + a5, где один из параметров at (в каждом случае) выражается через k. Обработка экспериментальных данных методом наименьших квадратов позволяет найти этот параметр и, следовательно, коэффициент k.
Данный блок состоит из полей для ввода параметров (для каждого металла предназначено отдельное поле), полей для вывода результата и диаграммы, на которой отображаются экспериментальные точки и теоретическая кривая, построенная по результатам расчета.
Блок «Расчет динамики сорбции и десорбции солей тяжелых металлов в почве при фильтрации в ней раствора биосурфактанта» предназначен для прогнозирования эффективности биоремедиации. Данный блок состоит из поля для ввода параметров, поля для вывода результата и диаграммы, на которой отображаются экспериментальные точки и теоретическая кривая, построенная по результатам расчета.
Процессы сорбции и десорбции солей тяжелых металлов изучались в модельных почвенных колонках. Наблюдалась динамика извлечения солей никеля и молибдена из почвы под действием водной эмульсии RhodOcoccus-биосурфактанта (4 г/л, штамм Rhodococcus ruber ИЭГМ 231), синтетического биосурфактанта Tween 80 и дистиллированной воды [2]. Загрязненная солями тяжелых металлов модельная почва находилась в стеклянных колонках (см. рис. 3). Выше нулевой отметки находилась жидкость, ниже нулевой отметки - загрязненная почва и проникшая в почву жидкость; поверхность почвы отвечает координате x 0 .
Количественное определение ионов никеля и молибдена, десорбированных из почвы в водную фазу, проводилось с помощью спектрофотометра диметилглиоксимным и роданидным методом, соответственно. Эффективность процесса десорбции с течением времени t выражалась в массовой доле ^(t) десорбированной соли от общего начального количества сорбированной соли.
Для решения использована математическая модель: начально-краевая задача для уравнения переноса. Математическая модель очистки почвы от солей тяжелых металлов основана на уравнениях материального баланса, описывающих десорбцию тяжелых металлов от почвы и повторную сорбцию их к почве:
dm/dt + V dm/dx = an - Pm, dn/dt = -an + Pm, x > 0, t > 0, n| =n0, ml A = 0, m| A = 0,
11=0 0 ’ 11=0 ’ I x=0 ’
где х - координата, отсчитываемая вдоль колонки, t - время, т - концентрация десорбированной соли, п - концентрация сорбированной соли, а - коэффициент сорбции, Р - коэффициент повторной сорбции, V - скорость течения жидкости в загрязненной почве, п 0 - начальная концентрация соли в почве. Решение записанной задачи имеет вид [5]
апо / \ апо
т ( х, t) =---(1 - ехр( - (а + РУ) )----ехр(-Р х / V ) х
а + Р а + Р
0, 0 < t < x/V,
t - x/V
| exp(-ax)I0 (2^/aPxx/V)(a + Pexp(-(a + P)(t- x/V-x)))dx, t > xjV,
a а
n(x, t) = n о 1----- +------ exp( - (a + p)t)
a+P a+P
aP n 0 a+P
exp(-P x / V) х
0, 0 < t < x/V,
t - x/V
| exp(-ax) I0 ( 2tJ aP xx/V )(l - exp (-(a + P) (t - x/V-x))) d x, t > xjV,
где 10 - модифицированная функция Бесселя. Экспериментально измеряемая
зависимость ) выражается следующим соотношением:
V |д* '
ц(t) =------I т(L, х^х,
Ln •'
0 0
где L - длина колонки, д * = lim д^) . Тогда
t
|i(t) =
aV д,
(a + P) L
0, 0 < t < L/V
1 - exp(-(a + P)t) t • a + P t-xjV
exp (-PL/V) | Д(t, x) dx, t > L/V 0 _
где
Д(t, x) = exp(-ax)I0 (2V0PLVV) a (t - L/V - x) +—^ [l - exp (-(a + P) (t - L/V - x))]
Заключение
В настоящей работе представлен пакет прикладных программ «Очистка загрязненной нефтью и тяжелыми металлами почвы с помощью биосурфактанта». Описаны основные блоки, составляющие представленный пакет прикладных программ, биомеханические и математические модели, на основе которых осуществляется решение задач в каждом из блоков пакета программ.
Благодарности
Исследования поддержаны грантами Министерства образования и науки РФ (ГК 16.518.11.7069, 2011-1.3-513-082) и Программы Президиума РАН “Молекулярная и
клеточная биология”.
Список литературы
1. Ившина И.Б., Куюкина М.С., Костарев С.М. Применение экологически безопасной экспресс-технологии очистки нефтезагрязненных почв и грунтов // Нефтяное хозяйство. - 2003. - № 9. -С. 116-118.
2. Коростина О.А., Куюкина М.С., Ившина И.Б., Костина Л.В., Осипенко М.А., Няшин Ю.И. Моделирование процессов фильтрации растворов солей тяжелых металлов и биосурфактанта в почве // Российский журнал биомеханики. - 2009. - Т. 13, № 3 (45). - С. 32-43.
3. Костина Л.В., Куюкина М.С., Ившина И.Б. Методы очистки загрязненных тяжелыми металлами почв с использованием (био)сурфактантов // Вестник Пермского университета. Биология. - 2009. -Вып. 10 (36). - С. 95-110.
4. Куюкина М.С., Ившина И.Б., Осипенко М.А., Няшин Ю.И., Коростина О.А. Модель нефтеотмывания загрязненного почвогрунта под действием Rhodococcus-биосурфактанта // Российский журнал биомеханики. - 2006. - Т. 10, № 1 (31). - С. 59-67.
5. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Изд-во МГУ, 1999. - 798 с.
6. Kuyukina M.S., Ivshina I.B., Makarov S.O., Litvinenko L.V., Cunningham C.J., Philp J.C. Effect of biosurfactants on crude oil desorption and mobilization in a soil system // Environment International. -2005. - Vol. 31. - P. 155-161.
7. Kuyukina M.S., Ivshina I.B., Osipenko M.A., Korostina O.A., Nyashin Yu.I. Mathematical modelling of crude oil-contaminated soil washing process using Rhodococcus-biosurfactant // Proc. Third IASTED Int. Conference “Environmental Modelling and Simulation”. Honolulu, USA. - 2007. - P. 77-81.
8. Miller R.M. Biosurfactant-facilitated remediation of metal-contaminated soils // Environmental Health Perspectives. - 1995. - Vol. 103, No. 1. - P. 59-62.
9. Mulligan C.N., Wang S. Remediation of a heavy metal-contaminated soil by a rhamnolipid foam // Engineering Geology. - 2006. - Vol. 85. - P. 75-81.
THE USE OF THE APPLICATION PACKAGE FOR SOLVING THE PROBLEMS CONCERNING OIL- AND HEAVY METAL-CONTAMINATED SOIL WASHING BY MEANS OF BIOSURFACTANT
O.A. Korostina, М^. Kuyukina, I.B. Ivshina, L.V. Kostina,
M.A. Osipenko, Y.I. Nyashin (Perm, Russia)
The application package has been designed for solving the following biomechanical problems: the calculation of the duration of oil-contaminated soil washing in the reservoir by means of the biogenic surface-active compound (biosurfactant), the calculation of the resistance coefficient for the biosurfactant filtration in heavy metal-contaminated soil, the prediction of the soil biological refinement (bioremediation) effectiveness. The basic package blocks are described. The biomechanical and mathematical models that govern the calculation are outlined.
Key words: biomechanics, application package, soil contamination, bioremediation, oil-contaminated soil, heavy metals, sorption, desorption, Rhodococcus-biosurfactant, mathematical modelling.
Получено 07 сентября 2011
