ГИАБ. Горный информационно-аналитический бюллетень / MIAB. Mining Informational and Analytical Bulletin, 2019;(12):191-204
УДК 579.66:51-7 DOI: 10.25018/0236-1493-2019-12-0-191-204
КИНЕТИЧЕСКИЕ МОДЕЛИ ОКИСЛЕНИЯ ДВУХВАЛЕНТНОГО ЖЕЛЕЗА АЦИДОФИЛЬНЫМИ ХЕМОЛИТОТРОФНЫМИ МИКРООРГАНИЗМАМИ (ОБЗОР)
Т.С. Хайнасова
Научно-исследовательский геотехнологический центр Дальневосточного отделения РАН, Петропавловск-Камчатский, Россия, e-mail: khainasova@yandex.ru
Аннотация: Исследования в области биогеотехнологии металлов ведутся уже в течение нескольких десятков лет. Несмотря на большой объем полученной информации, внимание исследователей уделяется интенсификации химико-технологических процессов, созданию высокопроизводительных аппаратных конструкций и повышению качества существующих производств. Для этого используется, в том числе метод математического моделирования. Многообразие химических, биологических, электрохимических, эксплуатационных различий и нелинейный процесс извлечения металлов в биовыщелачивающей среде осложняет моделирование подобных систем. Окисление двухвалентного железа ацидофильными хе-молитотрофными микроорганизмами — важная реакция в бактериально-химическом выщелачивании сульфидных руд. В статье представлены краткие результаты обзора существующих математических моделей, описывающих кинетику данной реакции и удельную скорость роста биомассы. Уравнения учитывают влияние таких факторов, как двухвалентное железо, клетки бактерий, температура, рН, газовый состав, тип биореактора, а также ингибирующее действие трехвалентного железа и прочих металлов (цинк, никель). Показано, что большинство из существующих моделей основывается на кинетике ферментативных реакций и представлено в виде модификаций уравнений Михаэлиса-Ментен и Моно. Несмотря на обилие исследований, касающихся данной тематики, до сих пор не существует комплексной модели, учитывающей все нюансы окисления.
Ключевые слова: биовыщелачивание, биоокисление, двухвалентное железо, скорость окисления, скорость роста микроорганизмов, кинетические модели, математическое моделирование, ацидофильные хемолитрофные микроорганизмы
Для цитирования: Хайнасова Т. С. Кинетические модели окисления двухвалентного железа ацидофильными хемолитотрофными микроорганизмами (обзор) // Горный информационно-аналитический бюллетень. - 2019. - № 12. - С. 191-204. DOI: 10.25018/0236-1493-2019-12-0191-204.
Kinetic model of ferrous iron oxidation with acidophil chemo-trophylitic microorganisms (review)
T.S. Khainasova
Geotechnological Scientific Research Center, Far Eastern Brunch of Russian Academy of Sciences, Petropavlovsk-Kamchatski Russia, e-mail: khainasova@yandex.ru
© Т.С. Хайнасова. 2019.
Abstract: Investigations into the bio-geo-technology of metals have been carried out for a few decades. With large bulk of data collected, researchers focus on intensification of chemical-engineering processes, design of high-production equipment and improvement of quality of existing practices. To this effect, the method of mathematical modeling is used among other things. The diversity of chemical, biological, electrochemical and operational features, as well as nonlinearity of metal extraction in the bio-leaching environment complicates modeling of such systems. Ferrous iron oxidation with acidophil chemo-trophylitic microorganisms—a critical reaction in the bacterial-chemical leaching of sulfide ore. This article gives a brief review of the existing mathematical models describing the kinetics of the said reaction and the specific velocity of the biomass growth. The equations take into account influence of such factors as ferrous iron, cells of bacteria, temperature, pH, gas content, type of bioreactor, as well as the inhibitory action of ferric iron and other metals (zinc, nickel). It s shown that the majority of the existing models are based on the kinetics of enzymatic reaction, and are represented by modifications of the Michaelis-Menten and Monod equations. Despite the abundant research in the given scientific area, no integrated model including all nuance of oxidation yet exists.
Key words: bioleaching, biooxidation, ferrous iron, oxidation rate, rate of growth of microorganisms, kinetic models, mathematical modeling, acidophil chemo-trophylitic microorganisms.
For citation: Khainasova T. S. Kinetic model of ferrous iron oxidation with acidophil chemo-tro-phylitic microorganisms (review). MIAB. Mining Inf. Anal. Bull. 2019;(12):191-204. [In Russ]. DOI: 10.25018/0236-1493-2019-12-0-191-204.
Введение
Исследования в области биогеотех-нологии металлов ведутся уже в течение нескольких десятков лет. Накоплен большой объем информации, касающейся биологических, физико-химических и технологических особенностей бактериально-химического выщелачивания и окисления ценных компонентов [1—5]. Существующие технологии требуют развития с целью повышения эффективности промышленного производства. Для этого используется, в том числе метод математического моделирования, который позволяет глубже изучить закономерности технологических процессов [6—9].
Множество химических, биологических, электрохимических, эксплуатационных параметров и нелинейный процесс извлечения металлов в биовыщелачива-ющей среде осложняет моделирование подобных систем [10]. Растворение сульфидных минералов при участии микроорганизмов — сложный химический многоступенчатый процесс. Его характер зави-
сит от ряда различных факторов (способ переработки, размер частиц руды или флотоконцентрата, плотность пульпы, скорость перемешивания, температура, адсорбция-десорбция, рН, газовый состав, состав питательных веществ и ингибиторов процесса, качественный состав, концентрация, биологическая активность, адаптационная способность и скорость роста клеток микроорганизмов и т.д.).
На сегодняшний день механизм бактериально-химического выщелачивания сульфидных руд подробно изучен. Он неодинаков для различных минералов, но связан, главным образом, с деятельностью железо- и сероокисляющих микроорганизмов, которые в результате своего метаболизма обеспечивают ускоренный перевод ценных металлов из нерастворимых форм (сульфидов) в растворимые (сульфаты). Исключение составляет биоокисление — процесс, в котором целевой металл, будучи в окисленной форме, в раствор не переходит, а остается в руде. При этом растворению подверга-
ются побочные химические элементы: мышьяк, железо, сера и прочие.
Роль микроорганизмов сводится к биологической регенерации трехвалентного железа ^е3+) — главного окисляющего химического агента в растворении минералов, в результате окисления двухвалентного железа ^е2+) (1), а также к образованию серной кислоты и поддержанию низких значений рН среды в ходе микробного окисления элементной серы (2) и ее восстановленных соединений [11]:
4Fe2+ + 02 + 4Н+ бактерии > 2Fe3+ + Н20,
(1)
0,125S8 +1,5 02 +
+ Н20
бактерии
> SO42- + 2Н+.
(2)
V = v„
К + Э
(3)
где V — скорость реакции, vm — максимальная скорость реакции, Э — концентрация субстрата, Кт — концентрация субстрата, при которой скорость реакции составляет половину максимальной скорости (константа Михаэлиса-Ментен).
Чаще всего экспериментальные данные выражаются через уравнение Моно (4). Оно описывает скорость роста бактерий в зависимости от концентрации лимитирующего субстрата:
бХа
л
= Кг
к+s
(4)
Описать выщелачивание сульфидных руд, учитывая все особенности процесса взаимодействия жидкой, твердой, газовой фаз и клеток микроорганизмов, как единую или обобщенную модель достаточно сложно.
Цель настоящей работы — провести анализ математических моделей, описывающих кинетику окисления двухвалентного железа ацидофильными железо-окисляющими микроорганизмами.
Анализ существующих
математических моделей
В основе существующих кинетических моделей окисления двухвалентного железа в периодическом и непрерывном режимах с использованием, как правило, бактерий А^№юЬасИи ferro-oxidans, Leptospirillum ferrooxidans и L. ferriphilum заложены уравнения Ми-хаэлиса-Ментена и Моно.
Уравнение Михаэлиса-Ментен (3), описывающее кривую субстратного насыщения в ходе моносубстратной ферментативной реакции, имеет вид:
где ц — удельная скорость роста биомассы; Ха — концентрация активной биомассы; t — время; Э — концентрация субстрата; цтах — максимальная удельная скорость роста биомассы; К, — концентрация субстрата, при которой удельная скорость роста биомассы принимает значение половины максимальной (константа Моно) [12—14].
Необходимо отметить, что уравнение Моно применимо только к начальным скоростям и не учитывает ингибирующие эффекты, влияющие на процесс. В то время как наличие токсического воздействия естественным образом может снижать кинетику.
Согласно ферментативной кинетике, в работах [12, 15, 16] представлены модели процессов окисления Fe2+, учитывающие явление ингибирования. В работах [15, 16] приведены соответствующие моделям уравнения, в основе которых заложены математические выражения скоростей роста бактерий (5) и окисления субстрата (6) в следующих видах:
^ = У(-^ I-ьх, dt { dt
Л
к + Б
(5)
(6)
где dX/dt — скорость роста бактерий (г VS/л ■ ч); dS/dt — скорость окисления
субстрата (г Fe2+/л ■ ч); Э — концентрация субстрата (г Fe2+/л); X — концентрация биомассы клеток (г VS/л); Y — коэффициент выхода (экономический коэффициент) (г VS/г Fe2+); qm — максимальная удельная скорость окисления субстрата (г Fe2+/г VS ■ ч) и Кз — константа полунасыщения (г Fe2+/л).
Распространенным типом ингибиру-ющего воздействия является так называемое самоингибирование, в котором скорость разложения субстрата, катализируемое ферментами, замедляется по причине высокой его концентрации. Модель субстратного ингибирования (модель Холдейна) описывается следующим образом (7) [16]:
dS
qmsx
(7)
dt Кв + Б + Б2 / К,
где К — константа самоингибирования (г/л). 1
Ингибирование, согласно теории кинетики ферментативных реакций, может быть конкурентным и неконкурентным. В первом случае конкурентный ингибитор связывается с активным центром фермента. Как правило, конкурентные ингибиторы по структуре похожи на субстрат и могут вытесняться из фермент-ингибиторного комплекса избытком субстрата. Такое взаимодействие не приводит к денатурации или инактивации фермента [13]. Степень ингибирования пропорциональна степени связывания конкурентного ингибитора с ферментом. Таким образом, уравнение, учитывающее это явление, можно представить так
qmsx
, (8)
[16]:
^в __
61 "(1 + 1С / Кис ж + в + в2/ К,
где 1с — начальная концентрация конкурентного ингибитора ^е3+) (г/л); КИс — константа конкурентного ингибирова-ния (г/л).
В случае неконкурентного ингибиро-вания ингибитор связывается не с актив-
ным центром фермента, а с другим его участком, формируя различные комплексы и приводя к изменению конформа-ции энзима. В данном случае скорость окисления железа описывается согласно формуле [16]:
£ = ,-^-г .(9)
dt (к, + S + S2 / K)(1 + Inc / K,„ )
где /nc — начальная концентрация неконкурентного ингибитора (например. Zn2+ или Ni2+) (г/л); Kiinc — константа неконкурентного ингибирования (г/л).
Третьим типом ингибирования является неконкурентоспособное (бесконкурентное) ингибирование. при котором происходит взаимодействие ингибитора только с комплексом, образованным ферментом и субстратом. Описывается согласно уравнению (10) [16]:
dS =_q^SX_
dt Ks +(S + S2/ K.) (1 + Iuc / Kiiuc) '
(10)
где /uc — начальная концентрация неконкурентоспособного (бесконкурентного) ингибитора (г/л); Kiiuc — константа неконкурентоспособного (бесконкурентного) ингибирования (г/л).
Ниже кратко представлены различные модели биоокисления двухвалентного железа и скорость роста бактерий на субстрате. которые основаны на вышеперечисленных кинетических уравнениях. Модель Kumar, Gandhi Кумар и Ганди предложили модель биоокисления двухвалентного железа бактериями T. ferrooxidans (сейчас Aci-dithiobacillus ferrooxidans) при температуре 30 °С, оптимальных значениях рН и явлении осаждения трехвалентного железа, принимая во внимание гибель микроорганизмов по причине воздействия Fe3+. Скоростями окисления Fe2+, осуществляющегося за счет мембран мертвых клеток. при рН 2,2—2,5 авторы пренебрегали.
Было показано увеличение рН в ходе окисления Fe2+ согласно стехиометри-ческому уравнению его окисления (1), а также понижение рН в результате осаждения железа и формирования ярозита. Авторы принимали во внимание вовлечение протонов водорода в реакции образования некоторых продуктов. Коэффициенты, полученные эмпирическим путем, были использованы для установления связи изменения рН с окислением субстрата и осаждением продуктов реакции. Кумар и Ганди также предположили, что скорости осаждения железа были пропорциональны концентрациям Fe3+ и обратно пропорциональны концентрациям протонов водорода. Максимальная удельная скорость роста зависела от рН; при этом оптимальное его значение составляло 2,2. Таким образом, предложенное кинетическое уравнение выглядело следующим образом:
Ито exP
K {pH - pHo )2 |Fe2+]
Ks + \_Fe2+ ]
(11)
На основании вышесказанных предположений были выведены следующие дифференциальные уравнения (12—15) для периодического процесса:
dX.
d- = -,d (Fe3+ J Xr,,
dFe2
dFe3 dt dpH
dt
= "ß:
dt
dFe2 dt dFe2 dt
= ~VXT
YXT /Fe2
(12)
(13)
Kf [Fe3+ ](10f ,(14) -ß2Kf [Fe3+] (10)
pH
(15)
Кинетические параметры были определены экспериментально. Данная модель с высокой точностью предсказывала концентрации двухвалентного железа, но слабо согласовывалась с данными по Fe3+ и рН [17, 18].
Модель Lacey, Lawson
Лейси и Лоусон исследовали действие температуры в диапазоне 20—31 °С на окисление двухвалентного железа T. ferrooxidans в периодическом режиме. В работе рассматривали неограниченное влияние кислорода, углекислого газа и питательных веществ на скорость роста бактерий. Исследователями было предложено уравнение (16), описывающее процесс следующим образом [19]: dS = ymSX dt Y (K + S )'
Авторами предложен способ, который позволял оценивать соответствующие константы скорости цт и K без непосредственного учета бактериальной популяции. Было обнаружено, что он дает значения цт и K, которые согласуются со значениями, точно определенными путем статистического регрессионного анализа экспериментальных данных.
Модель Kelly, Jones
В модели Келли и Джонса учитывалось влияние трехвалентного железа на рост микроорганизмов. Исследование проводили в периодическом режиме в колбах при температуре 30 °С и рН 1,6. Удельную скорость роста (ц), согласно предложенной модели, рассчитывали в соответствии со следующим уравнением [20]:
__[Fe2+]_
й Ks (1 + Fe3+ /Kp) + [Fe2+ ](l + Fe3+ /Kc)
(17)
Модель Braddock и др.
Модель Брэддока и соавторов представляет собой расширенную версию уравнения Келли и Джонса (17) за счет введения предельных (пороговых) значений. В модели учитывается тот факт, что при низких концентрациях железа бактериальный рост не происходит. Процесс окисления проводили при температуре 22 °С в непрерывном режиме.
1
Уравнение, вычисляющее удельную скорость роста бактерий (ц), имело следующий вид [20]:
11 к.
Модель Liu и др.
В работе Лю и соавторов представлены данные о влиянии инертных суспендированных частиц, концентраций трехвалентного железа и растворенного кислорода на кинетику окисления двухвалентного железа T. ferrooxidans в непрерывном режиме при 35 °С и рН 1,8. Процесс описывается уравнением [20, 21]:
ц = йта" (19)
* [Fe2+ ] + Ks (1 + Kj [Fe3+])
Модель Lizama, Suzuki
Авторы модели изучали окисление двухвалентного железа штаммом T. ferrooxidans в ходе конкурентного ингибирования концентрациями трехвалентного железа и клетками. Анализ кинетики выявил возможность синергетического ингибирования двумя компонентами. Скорость реакции окисления (v), согласно предложенной модели, выражается следующим образом [20, 22]:
v =
ks'[C ]\ Fe2
\_Fe2+ ] + Km {l + [C] / К,.' + \Fe3+ ] / Kif + ([C]\Fe3+ ]) / )аК/К,)
(20)
Модель Pesic и др.
В работе исследовали влияние температуры, рН, начальной концентрации бактерий и двухвалентного железа в присутствии T. ferrooxidans на скорость окисления. В ходе процесса непрерывно измеряли окислительно-восстановительный потенциал раствора. Было установлено, что процесс окисления Fe2+ можно выразить следующим образом [23, 24]:
d [Fe21 г пг
—= 1,62 • 1011 Cbact [н+1 [Fe21 pO2e-58J7/RT. (21)
Уравнение справедливо при атмосферных условиях, [Fe2+] < 0,001М, T < 25 °С и значениях рН выше 2,2. Если рН ниже 2,2, то скорость окисления не зависит от концентрации водорода.
Модель Karamanev, Nikolov (1991)
Авторами была предложена модель (22) для сравнения скоростей биоокисления Fe2+ в биопленочном биореакторе и в биореакторе со свободноплавающими клетками при различных концентрациях субстрата [23, 25].
f) Xfx 1
rs(sus) YS0 x
(22)
где X = 1 — 5/Э0.
Основное предположение состояло в том, что скорость процесса в обоих режимах культивирования описывается одним и тем же кинетическим уравнением, на-
пример, моделью Моно. Когда субстрат полностью расходуется (х = 1), соотношение между скоростью реакции в случае биопленки и скоростью суспендированной культуры становится равным ХЛх /YS0. Экспериментально обнаружено, что ХЛх /У50 = 4,3.
Модель Karamanev, (1992)
Модель (23) описывает непрерывный режим окисления железа при температуре 29 °С, рН 1,8-2,0 [20].
[Fe2+] + Ks +[Fe3+]
Ks K
(23)
K
Модель Huberts
Модель Хьюбертса (24) описывает непрерывный режим окисления железа при температуре 30 °С и рН 2,0 [20]:
с \
~rFe2 - а1
Po
л
V Гв + po2 ,
\_Fe2+ ]
^Fe2+ ] + KR
\Fe3+]^
1 +
K
(24)
;)
Модель Gomez и др.
Описывает процесс окисления двухвалентного железа T. ferrooxidans в периодическом режиме при температуре 30 °С и начальном рН 2,0. Модель (25) основывается на уравнении, предложенном Лю и соавторами (19), и позволяет вычислять удельную скорость роста как функцию концентраций субстрата (Fe2+) и продукта (Fe3+). При этом предполагается, что бактерии конкурентно ингибируются Fe3+. Уравнение представлено в виде [26]:
VmaxFe2 +
Ио =-
KS + Fez
+ KFe
(25)
Модель Crundwell
Крандвелл сформулировал теоретическую модель, основанную на фундаментальной хемоосмотической теории электрон/протонного транспорта у А. ferrooxidans. Признавая тот факт, что механизм служит генератором энергии для метаболической активности клетки, процесс окисления железа был выражен следующим образом [20, 24]:
-r
= k
[Fe-] / [н + ]
v
-[Fe2+]/[H +] + K [Fe3+]
Ke 2
\ Fe
Po
у
kB + Po.
(26)
2
Модель Das и др.
Авторы исследовали процесс окисления двухвалентного железа с учетом влияния рН, концентраций Fe2+, Fe3+ и микроорганизмов. Было установлено, что скорость окисления возрастала с увеличение концентрации биомассы клеток, но уменьшалась с ростом Fe3+ в растворе. С увеличением концентрации Fe2+ до 0,6 г и рН до 2,0 скорость реакции возрастала. После преодоления порога данных значений параметров отмечали тенденцию к снижению кинетики. Осаждение железа происходило
2
2
при увеличении рН, концентрации бактерий и Fe2+. Константы реакции для каждого случая были подсчитаны, используя уравнение скорости первого порядка. Скорость окисления, согласно модели, была описана так:
к' (рН)( концентрация биомассь!)
Скорость окисления железа = --———-- , (27)
(Fe ДРе J
где к' — константа скорости реакции псевдо первого порядка [27]. Модель Nemati, Webb
Описывается кинетика биоокисления двухвалентного железа Т. ferrooxidans с учетом влияния температуры, концентраций Fe2+ и бактерий. В модели (28) авторы учитывали ингибирование последних двух параметров:
c/[Fe2+] К0е^[х][Ре2+]
dt
(28)
Fe
1 + K
где К0 (6438 кг ■ м"3 ■ ч'Укл ■ мл"1); Fa (68,4 кДж ■ моль"1); К'т (0,0672 кг ■ м"3); K'¡ (2,68 ■ 107 кл ■ мл-1); а (26,1 кг ■ м-3); (3 (7,8 ■ 108 кл ■ мл-1) — коэффициенты кинетической модели в отсутствии трехвалентного железа [28].
Модель Hansford
Упрощенная версия модели ингибированного трехвалентным железом роста бактерий Leptospirillum (уравнение Келли и Джонса (17)) была предложена Хэнсфордом. Удельная скорость потребления двухвалентного железа рассчитывалась следующим образом:
1 max
g 2+ =-%-=- ■ [20] (29)
[Fe3
'^[Fe2
Модель Boon и др.
Бун и соавторы предполагали, что потребление субстрата (Fe2+) и кислорода — ключевые процессы. Показана строгая стехиометрическая связь между потреблением кислорода и бактериальным окислением двухвалентного железа, позволяющая вычислять скорость последнего.
Эксперименты проводили в периодическом режиме. Установлено, что возникает пороговый уровень Fe2+, при котором потребление кислорода прекращается. В связи с этим для расчетов было использовано уравнение Моно для конкурентного ингиби-рования трехвалентным железом и пороговой концентрации двухвалентного железа в следующем виде [29]:
о =_q°2'max__(30)
2 , Ks Ks[Fe3+]
Известно, что энергия, образуемая из окисления двухвалентного железа, используется бактериями для роста. Этот процесс описывает уравнение Пирта:
l'y
Fe утаи S X ,
'SX
где -г ,+ — скорость окисления двухвалентного железа (моль Ре2+ ■ л-1 ■ с-1); гх — скорость бактериального роста (С-моль ■ л-1 ■ с-1); Узхтах — максимальный коэффициент выхода биомассы в присутствии двухвалентного железа (С-моль ■ моль Ре2+-1); т3 — коэффициент сохранения биомассы в присутствии двухвалентного железа (моль Ре2+ ■ С-моль ■ л-1 ■ с-1); Сх — концентрация биомассы (С-моль ■ л-1).
Уравнение (32) связывает между собой удельную скорость потребления кислорода биомассой ^0 , моль О2 ■ С-моль-1с-1) и удельную скорость роста клеток (ц, с-1):
Уп = + т0. (32)
°2 утах 0 /
'пх
Значения У0хтах (максимальный коэффициент выхода биомассы по кислороду, С-моль ■ моль 02-1) и гпо (коэффициент сохранения в присутствии кислорода, моль 02 ■ С-моль-1 ■ с-1) получают путем измерений и подставляют в уравнение:
11 m,
+ . (33)
о
Y Ymax ц
'ox 'ox г1
В данном уравнении YOX — это выход биомассы по кислороду (С-моль ■ моль O2-1), который вычисляют из скоростей потребления кислорода и углекислого газа:
YOX = rCO2 /rO2.
Если брать в расчет степень редукционного баланса (-rFe2+ = -4rO — 4,2rCO , в котором -rFe2+ — скорость окисления двухвалентного железа (моль Fe2+ ■ л-1 ■ с-1);-rO — скорость потребления кислорода (моль 02 ■ лр-ра-1 ■ с-1); -rCO — скорость потребления углекислого газа (моль С02 ■ лр-ра-1 ■ с-1)), коэффициенты максимального выхода для железа (YSXmax) и кислорода (YOXmax) взаимосвязаны и достаточно определить только один из них:
4 у max
у max _ 4 'SX (34)
0X 1 - 4,2YSXax '
То же самое верно для mo и mS:
m0 = ms. (35)
0 4
Если применяется уравнение Пирта, то удельная скорость роста напрямую зависит от удельной скорости потребления кислорода. При условии, что скорости максимального удельного роста и потребления кислорода равны, справедливо следующее уравнение:
Vmax = ^ + m0. (36)
'ox
Учитывая уравнения (30), (32) и (36) Бун и соавторы вывели следующую кинетическую модель:
И + m Ymax
. __"ma^ 0 OX__m Ymax /Q7\
' =, + ks + (37)
[Fe2*]"[Fe2+], K, [Fe2+]-[Fe=']i
Это уравнение учитывает тот факт, что при концентрациях двухвалентного железа, близких к пороговому значению ([Fe2+]t), будет иметь место отрицательная удельная скорость роста (в результате разрушения бактерий) и продолжаться потребление кислорода и Fe2+.
Модель Мегаипе и др.
Авторы исследовали влияние окислительно-восстановительного потенциала на скорость окисления двухвалентного железа A. ferrooxidans при температуре 30 °С и рН 1,8. Результаты позволили сформулировать кинетическую модель, которая учитывает влияние Eh, концентрации двухвалентного железа и ингибирующий эффект трехвалентного железа. Было установлено, что при значениях Eh ниже 650 мВ окислительная активность бактерий в основном зависит от концентрации двухвалентного железа. При значениях Eh между 650 и 820 мВ скорость окисления регулируется, главным образом, ингибирующим воздействием трехвалентного железа. При Eh более 820 мВ окислительная активность сильно усугубляется и полностью ингибирует-ся при 840 мВ. Уравнение, описывающее данный процесс [20, 30]:
К* ехр
2ИТ
{ -Е0)]|
1 - ехр
пГ ИТ
{ - Е)
К *
1Ге ]
пГ ИТ
{ - Е0)
(38)
Модель Penev, Kaгamanev
Исследование кинетики проводили в биореакторе с циркулирующем слоем в периодическом режиме при умеренных и высоких концентрациях железа, оптимальном значении рН (ниже 1) и с использованием бактерии L. ferriphilum. По результатам экспериментов была разработана кинетическая модель, в основе которой заложено уравнение Моно с учетом конкурентного ингибирования.
Путем различных преобразований конечный вид уравнения был представлен в следующем виде [31]:
Яз =■
Ятех 3
(39)
Б + Кз + КЗ Модель Шгт '1 и др.
Кинетику окисления железа изучали в биореакторе со взвешенным слоем при использовании культуры с доминирующей бактерией L. ferriphilum. Комбинированная модель с учетом конкурентного и неконкурентного ингибирования была разработана для процесса при одновременном присутствии в растворе Fe3+, Zn2+ и М2+. Уравнение скорости окисления Fe2+ с подсчитанными кинетическими константами имело следующий вид [32]:
_ ^ -dt ~
(2,15 ± 0,25) [Ре2+"
0,413 (1 + Fe3+ / 0,828) + ^е2+ ] + ^е2+ ]2 / 8,65 (1 + 7п2+ /49,1 + 1 + №2+ / 62,6)
Заключение
Биоокисление двухвалентного железа — ключевая реакция в бактериально-химическом выщелачивании сульфидных руд. Разработка его кинетических моделей важна для интенсификации химико-технологических процессов, создания высокопроизводительных аппаратных конструкций и повышения качества существующих
производств. Из-за нелинейности процессов в биовыщелачивающих системах и влияния различных физико-химических и биологических факторов на кинетику окисления математические модели претерпевают изменения уже в течение 100 лет.
Различными авторами исследовалась кинетика окисления с участием небольшой группы ацидофильных железоокис-ляющих бактерий, как правило, представленных монокультурами А. ferrooxidans, L. ferrooxidans и L. ferriphilum. В ходе процессов определяли скорости окисления двухвалентного железа и роста биомассы микроорганизмов с учетом влияния таких факторов, как двухвалентное железо, бактериальные клетки, темпера-Номенклатура
Символы и обозначение Единицы измерения
а — эмпирическая константа
Р1, Р2 — эмпирические коэффициенты л ■ г-1
[С], СЬас(, X — концентрация бактерий г ■ л-1, мг
[Ре2+], Б — концентрация двухвалентного железа г ■ л-1, ммоль
[Ре2+] — пороговая концентрация двухвалентного железа моль ■ л-1
[Ре3+] — концентрация трехвалентного железа г ■ л-1, ммоль
[Н+] — концентрация протонов водорода кг ■ м-3
К1 — безразмерная константа модели
К — константа скорости реакции г ■ л —1 ■ ч-1
К' — константа ингибирования (К), деленная на количество ингибитора на клетку мг/мл
К , К , К, К,, К,, К , К — константы ингибирования р' С 1' Г if' г
К, Кт, КРб2+, К — константа полунасыщения, константа сродства к субстрату (Ре2+) г ■ л-1, моль ■ л-1
к, qFe2+max — максимальная удельная скорость потребления двухвалентного железа кг ■ м-3 ■ ч-1, моль Ре2+ ■ (С-моль)-1 ■ ч-1
к3' — к3-кратное количество фермента на клетку
ц, цс — удельная скорость роста биомассы ч-1
цт — максимальная оптимальная удельная скорость роста биомассы ч-1
Цт, Цтах — максимальная удельная скорость роста биомассы ч-1
ц — константа скорости гибели биомассы клеток л ■ г-1 ■ ч-1
рН0 — оптимальное для роста бактерий значение рН (2,2)
тура, рН, газы, тип биореактора, а также ингибирующего действия трехвалентного железа и прочих металлов (цинк, никель). Большинство из существующих моделей основывается на кинетике ферментативных реакций и представлено в виде модификаций уравнений Михаэли-са-Ментен и Моно. В основе некоторых уравнений лежит хемоосмотическая теория и учитывается механизм переноса электронов и протонов в ходе окисления Fe2+. Также кинетика выражается через уравнение Пирта. Несмотря на обилие исследований, касающихся данной тематики, до сих пор не существует комплексной модели, учитывающей все нюансы окисления.
pO — парциальное давление кислорода атм
qO — удельная скорость потребления О2 моль О2 ■ (С-моль)-1 ■ ч-1
qO ,тах — максимальная удельная скорость потребления О2 моль О2 ■ (С-моль)-1 ■ ч-1
qS — удельная скорость окисления двухвалентного железа г ■ л-1 ■ ч-1
qmax — максимальная удельная скорость окисления двухвалентного железа г ■ л-1 ■ ч-1
R — универсальная газовая постоянная
-Гре2., qFe2. — удельная скорость потребления двухвалентного железа ммоль ■ ч-1, кг ■ м-3 ■ ч-1, моль Fe2+ ■ (С-моль)-1 ■ ч-1
гз(Лх) — скорость реакции в биопленочном биореакторе г ■ л-1 ■ ч-1
г ^ — скорость реакции в биореакторе со свободноплавающими клетками г ■ л-1 ■ ч-1
Б0 — концентрация субстрата на входе г ■ л-1
Б — концентрация субстрата на выходе г ■ л-1
— концентрация общего железа г ■ л-1
T — абсолютная температура К
t — время ч
Хт — общая концентрация биомассы кл ■ л-1
Хцх — концентрация клеток в биопленке г ■ л-1
Y — биомасса, образующаяся на грамм окисленного железа г/г
YXT/Fe2+ — рост свободноплавающих клеток в присутствии двухвалентного железа кл ■ г Fe2+-1
СПИСОК ЛИТЕРАТУРЫ REFERENCES
1. Watling H. R. Review of biohydrometallurgical metals extraction from polymetallic mineral resources. Minerals. 2015;(5)1:1-60. DOI: 10.3390/min5010001.
2. Rawlings D. E., Johnson D. B. The microbiology of biomining: development and optimization of mineral-oxidizing microbial consortia. Microbiology. 2007;153:315-324. DOI: 10.1099/ mic.0.2006/001206-0.
3. Olson G. J., Brierley J. A., Brierley C. L. Bioleaching review part B: progress in bioleaching: applications of microbial processes by the minerals industries. Applied Microbiology and Biotechnology. 2003;63:249-257. DOI 10.1007/s00253-003-1404-6.
4. Johnson D. B. Minireview. Biodiversity and ecology of acidophilic microorganisms. FEMS Microbiology and Ecology. 1998;27:307-317. DOI: 10.1111/j.1574-6941.1998.tb00547.x.
5. Watling H. R. The bioleaching of nickel-copper sulfides. Hydrometallurgy. 2008;91:70-88. DOI: 10.1016/j.hydromet.2007.11.012.
6. Casas J. M., Martinez J., Moreno L., Vargas T. Bioleaching model of a copper-sulfide ore bed in heap and dump configurations. Metallurgical and Materials Transactions B. 1998;29B: 899-909.
7. Botane P., Brochot S., D'Hugues P., Spolaore P. Material size distribution in concurrent bioleaching and precipitation: Experimental procedure and modeling. Hydrometallurgy. 2013;133:7-14. DOI: 10.1016/j.hydromet.2012.11.008.
8. Ahmadi A., Ranjbar M., Schaffie M., Petersen J. Kinetic modeling of bioleaching of copper sulfide concentrates in conventional and electrochemically controlled systems. Hydrometallurgy. 2012;127-128:16-23. DOI: 10.1016/j.hydromet.2012.06.010.
9. Petersen J. Modelling of bioleach processes: Connection between science and engineering. Hydrometallurgy. 2010;104:404-409. DOI: 10.1016/j.hydromet.2010.02.023.
10. Ahmadi A., Hosseini M. R. A fuzzy logic model to predict the bioleaching efficiency of copper concentrates in stirred tank reactors. International Journal of Nonferrous Metallurgy. 2015;4:1-8. DOI: 10.4236/ijnm.2015.41001.
11. Schippers A., Sand W. Bacterial leaching of metal sulfides proceeds by two indirect mechanisms via thiosulfate or via polysulfides and sulfur. Applied and Environmental Microbiology. 1999;65(1):319-321.
12. Rittman B. E., McCarty P. L. Environmental Biotechnology: Principles and Applications. McGraw-Hill, 2001. 754 p.
13. Комов В. П., Шведова В. Н. Биохимия. — М.: Дрофа, 2008. — 638 с.
13. Komov V. P., Shvedova V. N. Biokhimiya [Biochemistry], Moscow, Drofa, 2008, 638 p.
14. Моделирование микробных популяций. Лекция 11. [Электронный ресурс]. URL: http:// www.library.biophys.msu.ru/LectMB/lect11.htm (дата обращения: 16.09.2019).
14. Modelirovanie mikrobnykh populyatsiy. Lektsiya 11 [Microbial population modeling. Lecture 11], available at: http://www.library.biophys.msu.ru/LectMB/lect11.htm (accessed 16.09.2019).
15. Nurmi P., Ozkaya B., Kaksonen A. H., Tuovinen O. H., Puhakka J. A. Inhibition kinetics of iron oxidation by Leptospirillum ferriphilum in the presence of ferric, nickel and zinc ions. Hydro-metallurgy. 2009;97:137-145. DOI: 10.1016/j.hydromet.2009.02.003.
16. Nurmi P. Oxidation and control of iron in bioleaching solutions. Thesis for the degree of Doctor of Technology, 2009. p. 83.
17. Kumar S. R., Gandhi K. S. Modelling of Fe2+ oxidation by Thiobacillus ferrooxidans. Applied Microbiology and Biotechnology. 1990;33:524-528.
18. Haddadin J., Dagot C., Fick M. Models of bacterial leaching. Reviews. Enzyme and Microbial Technology. 1995;17:290-305.
19. Lacey D. T., Lawson F. Kinetics of the liquid-phase oxidation of acid ferrous sulfate by the bacterium Thiobacillus ferrooxidans. Biotechnology and Bioengineering. 1970;12:29-50. DOI: 10.1002/bit.260120104.
20. Ojumu T. V., Petersen J., Searby G. E., Hansford G. S. A review of rate equations proposed for microbial ferrous-iron oxidation with a view to application to heap bioleaching. Hydrometallurgy. 2006;83(1):21-28. DOI: 10.1016/j.hydromet.2006.03.033.
21. Meruane G., Salhe C., Wiertz J., Vargas T. Novel electrochemical-enzymatic model which quantifies the effect of the solution Eh on the kinetics of ferrous iron oxidation with Acidithio-bacillus ferrooxidans. Biotechnology and Bioengineering. 2002;80(3):280-288. DOI: 10.1002/ bit.10371.
22. Lizama H. M., Suzuki I. Synergistic competitive inhibition of ferrous iron oxidation by Thiobacillus ferrooxidans by increasing concentrations of ferric iron and cells. Applied and Environmental Microbiology. 1989;55(10):2588-2591.
23. Jensen A. B., Webb C. Ferrous sulphate oxidation using Thiobacillus ferrooxidans: a review. Process Biochemistry. 1995;30(3):225-236.
24. Nemati M., Harrison S. T. L., Hansford G. S., Webb C. Review. Biological oxidation of ferrous sulphate by Thiobacillus ferrooxidans: a review on the kinetic aspects. Biochemical Engineering Journal. 1998;1:171-190.
25. Karamanev D., Nikolov L. A comparison between the reaction rates and free suspended cells bioreactors. Bioprocess Engineering. 1991;6:127-130.
26. Gomez J. M., Caro I., Cantero D. Kinetic equation for growth of Thiobacillus ferrooxi-dans in submerged culture over aqueous ferrous sulphate solutions. Journal of Biotechnology. 1996;48:147-152.
27. Das T., Chaudhury G. R., Ayyappan S. Use of Thiobacillus ferrooxidans for iron oxidation and precipitation. BioMetals. 1998;11:125-129.
28. Nemati M., Webb C. Inhibition effect of ferric iron on the kinetics of ferrous iron. Biotechnology Letters. 1998;20(9):873-877.
29. Boon M., Ras С., Heijnen J. J. The ferrous iron oxidation kinetics of Thiobacillus ferrooxidans in batch cultures. Applied Microbiology and Biotechnology. 1999;51:813-819.
30. Meruane G., Salhe C., Wiertz J., Vargas T. Novel electrochemical-enzymatic model which quantifies the effect of the solution Eh on the kinetics of ferrous iron oxidation with Acidithio-
bacillus ferrooxidans. Biotechnology and Bioengineering, 2002;80:280-288. DOI: 10.1002/ bit.10371.
31. Penev K., Karamanev D. Batch kinetics of ferrous iron oxidation by Leptospirillum ferriphi-lum at moderate to high total iron concentration. Biochemical Engineering Journal, 2010;50:54— 62. doi:10.1016/j.bej.2010.03.004.
32. Nurmi P., Özkaya B., Kaksonen A. H., Tuovinen O. H., Puhakka J. A. Inhibition kinetics of iron oxidation by Leptospirillum ferriphilum in the presence of ferric, nickel and zinc ions. Hydro-metallurgy, 2009;97:137-145. doi:10.1016/j.hydromet.2009.02.003. [¡223
ИНФОРМАЦИЯ ОБ АВТОРЕ
Хайнасова Татьяна Сергеевна — канд. биол. наук, старший научный сотрудник, e-mail: khainasova@yandex.ru, Научно-исследовательский геотехнологический центр ДВО РАН.
INFORMATION ABOUT THE AUTHOR
T.S. Khainasova, Cand. Sci. (Biol.), Senior Researcher, e-mail: khainasova@yandex.ru, Geotechnological Scientific Research Center, Far Eastern Brunch of Russian Academy of Sciences, 683002, Petropavlovsk-Kamchatsky, Russia.
A
ОТДЕЛЬНЫЕ СТАТЬИ ГОРНОГО ИНФОРМАЦИОННО-АНАЛИТИЧЕСКОГО БЮЛЛЕТЕНЯ
(СПЕЦИАЛЬНЫЙ ВЫПУСК)
ПРОБЛЕМЫ ОСВОЕНИЯ ГЕОРЕСУРСОВ ДАЛЬНЕГО ВОСТОКА
(2019, СВ 30, 160 с.)
Выпуск содержит статьи, отражающие современное состояние исследований по проблемам освоения георесурсов Дальнего Востока. Рассмотрен широкий круг актуальных вопросов, связанных с добычей и переработкой минерального сырья, ведением работ в сложных горно-геологических и горнотехнических условиях, и намечены направления эффективного и безопасного освоения недр Дальнего Востока. Представлены новые результаты фундаментальных и прикладных исследований в области геотехнологии, глубокой переработки минерального сырья, геомеханики, горной теплофизики, геоэкологии и других разделов горных наук, играющих важную роль для эффективного и безопасного освоения недр. Для научных и инженерно-технических работников горно-добывающей отрасли, аспирантов и студентов горных специальностей вузов.
PROBLEMS OF DEVELOPMENT OF GEORESOURCES OF THE FAR EAST
The issue contains articles reflecting the current state of research on the development of geo-re-sources of the Far East. A wide range of topical issues related to the extraction and processing of mineral raw materials, the conduct of work in complex geological and mining conditions, and outlined directions for effective and safe development of the subsoil of the Far East. New results of fundamental and applied research in the field of Geotechnology, deep processing of mineral raw materials, geomechanics, mining Thermophysics, Geoecology and other branches of mining Sciences, which play an important role for effective and safe development of subsoil, are presented. For scientific and engineering workers of the mining industry, graduate students and students of mining specialties of universities.