Статьи - Статьи

ВЕСТНИК ДАГЕСТАНСКОГО НАУЧНОГО ЦЕНТРА. 2007. № 29. С. 59-68.

УДК 541.183

ИММОБИЛИЗАЦИЯ ТЯЖЕЛЫХ МЕТАЛЛОВ ИЗ ВОДНЫХ РАСТВОРОВ КАРБОНАТЫМИ ПОРОДАМИ ДАГЕСТАНА

Бутаев А.М., Магомедбеков У.Г., Мухучев А.А., Магомедрасулова Х.М.

Оценена возможность использования карбонатных пород Дагестана для извлечения Pb, Cd, Zn, Co, Cu и Ni из водных растворов. Дана термодинамическая и кинетическая интерпретация полученных результатов. Предложен новый механизм иммобилизации тяжелых металлов из водных растворов карбонатными породами.

Известно [1], что в саморегуляции и самоочищении поверхностных и подземных вод важную роль играют природные геохимические барьеры, представляющие собой участки земной коры, где на коротких расстояниях происходит резкое уменьшение интенсивности миграции химических элементов. Эту идею позаимствовали исследователи у природы и с успехом применяют на практике для очистки питьевых и технических вод от тяжелых металлов (ТМ) путем создания искусственных геохимических барьеров c помощью глин, опок, цеолитов, глауконитов, диатомитов, трепелов, торф, гидратированных оксидов, карбонатных пород и других веществ [2-6]. В литературе эти вещества известны под общим названием "природные сорбционные материалы", но по превалирующему механизму иммобилизации металлов их можно разделить на осадители и сорбенты. В первом случае происходит перевод металла из биологически активного состояния в труднорастворимое или комплексное соединение, во втором – его физическое или ионообменное поглощение. Высокие сорбционные свойства, дешевизна и широкая распространенность в природе делают многих природных сорбентов экономически целесообразным сырьем в технологии водоподготовки.

Целью настоящей работы явилось выяснение возможности использования природных карбонатных пород Дагестана для очистки водных объектов от соединений тяжелых металлов. Эти исследования могут в дальнейшем представлять и самостоятельный научный интерес по изучению влияния различных региональных геохимических барьеров на формирование химического состава природных вод Дагестана. Так, горные породы Дагестана представлены широким спектром различных формаций (см., например, [7]), каждая из которых, выступая в качестве геохимического барьера, так или иначе, оказывает влияние на химический состав водной среды, что в одних районах может инициировать эндемические заболевания, обусловленные избытком ТМ в питьевой воде, в других – их недостатком.

 

Материалы и методика исследования

В качестве природного карбонатного материала использовались доломиты и известняки из верхнего мелового отложения Гумбетовского и Унцукульского районов. Их химические составы приведены в табл. 1.

Таблица 1. Химический состав доломитов и известняков

Район отбора проб

Содержание, масс.%

SiO2

Al2O3

Fe2O3

CaO

MgO

Na2O+K2O

п.п.п

Доломиты

Гумбетовский*

1,25

0,30

0,22

31,2

19,2

0,63

47,2

Унцукульский*

0,81

0,38

0,45

31,1

21,3

0,10

45,9

Унцукульский**

0,58

0,24

0,23

30,6

20,8

0,14

47,4

Среднее

0.88

0.31

0.30

31.0

20.4

0.29

46.8

Известняки

Гумбетовсий*

4.52

0.85

0.54

42.6

7.9

0.38

43.3

Гумбетовсий**

3.10

0.47

0.33

43.8

8.2

0.17

44.0

Среднее

3.81

0.66

0.44

43.2

8.0

0.28

43.6

Примечание. п.п.п. – потери при прокаливании (0.5-1.5% Н2О, <0.5% серы и другие летучие вещества, остальное СО2). Анализ выполнен: * – на заводе "Стекловолокно" г. Махачкала, ** – в НПО «Технология» г. Обнинск.

 

При изучении подобных гетерогенных процессов возникает ряд специфических сложностей, влияющих на воспроизводимость экспериментальных данных от опыта к опыту. В основном они связаны с необходимостью обеспечения химической и гранулометрической однородности сорбента в параллельных опытах. В данной работе это достигалось следующим образом. Куски породы подвергались дроблению, измельченный материала промывали в дистиллированной воде для удаления пылевидной составляющей (3-4 разовое взмучивание с последующей декантацией), высушивали при температуре 110±5 °С, а для опытов брали частицы округлой формы диметром 1-2 мм, полученные путем последовательного трехкратного скатывания измельченного и высушенного материала по наклонной плоскости (лист бумаги).

Исследование процесса сорбции ТМ проводили в статических условиях методом отдельных навесок с использованием модельных растворов, приготовленных путем растворения точных навесок солей (в основном нитратов) Cd, Co, Cu, Ni, Zn и Pb классификации «х.ч.» в дистиллированной воде. В конические колбы с 100 мл раствора заданной концентрации вводили разное количество карбонатной породы, таким образом, чтобы весовое отношение твердая фаза:раствор (Т:Ж отношение) составляло 0.01, 0.02, 0.04, 0.07 и 0.1 кг/кг Н2О. Для цинка и свинца снимали изотермы сорбции при Т:Ж = 0.04 кг/кг Н2О. Колбы выдерживали при периодическом перемешивании в лабораторных условиях (при 20±3 оС) в течение определенного времени, после чего раствор фильтровали ("синяя лента") и фильтрат анализировали на атомно-абсорбционном спектрофотометре модели АА-6701 фирмы SHIMADZU (Япония). Относительная погрешность измерений не превышала 2%.

 

Результаты исследований

На первом этапе были проведены исследования по оценке эффективности извлечения ТМ из водных растворов доломитом и известняком. Для этого колбы выдерживали в течение 10 суток (время, достаточное для установления равновесия) и по убыли концентрации ТМ в растворе рассчитывали сорбционную емкость поглощения (J), коэффициент распределения металла между раствором и карбонатом (КD) и степень извлечения металла из раствора (W):

J = [(Сор)·V]/m, мг/г; КD = J/Ср, л/г; W = [(Сор)·100]/Со, %,

где Со и Ср – исходная и равновесная концентрация металла, мг/л; V – объем раствора, мл; m – масса образца породы, г. Результаты этой серии экспериментов обобщены в табл. 2 и 3.

Таблица 2. Равновесные концентрации металлов в растворе, мг/л (t=10 суток)

Весовое отношение

Т:Ж, кг/кг Н2О

Доломит

Известняк

Zn

Pb

Zn

Pb

Cu

Ni

Co

Cd

Исходный раствор

0.01

0.02

0.04

0.07

0.10

3.50

-

2.05

1.53

1.07

0.80

3.50

-

0.45

0.27

0.15

0.12

3.50

2.25

1.04

0.55

0.42

0.24

3.50

0.17

0.13

0.06

0.04

0.03

3.50

-

-

0.23

-

-

3.50

-

-

2.06

-

-

0.75

-

-

0.08

-

-

0.75

-

-

0.02

-

-



 

 

 

 

 

 


Таблица 3. Сорбционные свойства карбонатов

Параметр

сорбции

Т:Ж, кг/кг Н2О

0.02

0.04

0.07

0.10

Доломит, Zn

J, мг/г

KD, л/г

W, %

72.5

35.4

44.4

49.2

32.2

56.3

34.7

32.4

69.4

27.0

33.8

77.1

Доломит, Pb

J, мг/г

KD, л/г

W, %

152.5

339

87.1

80.7

299

92.3

47.8

319

95.7

33.8

282

96.6

Известняк, Zn

J, мг/г

KD, л/г

W, %

123.0

118

70.3

73.8

134

84.3

44.0

105

88.0

32.6

136

93.1

Известняк, Pb

J, мг/г

KD, л/г

W, %

168.5

1296

96.3

86.0

1433

98.3

49.4

1254

98.8

34.7

1157

99.1

Известняк, Cu

J, мг/г

KD, л/г

W, %

-

-

-

81.7

355

93.5

-

-

-

-

-

-

Известняк, Co

J, мг/г

KD, л/г

W, %

-

-

-

16.7

209

89.3

-

-

-

-

-

-

Известняк, Cd

J, мг/г

KD, л/г

W, %

-

-

-

18.2

912

97.3

-

-

-

-

-

-

Известняк, Ni

J, мг/г

KD, л/г

W, %

-

-

-

36.0

17

41.1

-

-

-

-

-

-



 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнение абсолютных величин сорбции показывает более высокую поглотительную способность известняка к цинку и свинцу по сравнению с доломитом (табл. 3).Так, при прочих равных условиях, КD известняка в 4 раза выше, чем доломита. Надо, однако, заметить, что с увеличением Т:Ж отношения, т.е. по мере приближения к условиям, неспособных «насытить» сорбент до предела (что, как правило, имеет место в естественных условиях), различия в емкостях поглощения доломита и известняка нивелируются.

На рис. 1 приведены изотермы сорбции цинка и свинца на известняке. Как и следовало ожидать, с ростом концентрации металла в растворе увеличивается сорбционная емкость сорбента. Обратим внимание и на некоторые различия данных табл. 3 и рис. 1 (рисуночные величины J примерно на 10% меньше табличных), что лишний раз указывает на влияние предыстории подготовки сорбента на его сорбционные свойства.

Рис. 1. Изотермы сорбции цинка и свинца известняком (Т:Ж = 0.04, t = 10 суток)

Рис. 1. Изотермы сорбции цинка и свинца известняком (Т:Ж = 0.04, t = 10 суток)

Сорбционная способность карбонатных пород в отношении изученных металлов сильно различается. При этом отметим, что оценка эффективности сорбента для извлечения металлов из водных растворов с помощью удельной сорбционной емкости может привести к ошибочным выводам. Так, оценка эффективности иммобилизации ТМ известняком с помощью J для случая Т:Ж = 0.04 дает следующий ряд сорбции Pb > Cu > Zn > Ni > Cd > Co, а с помощью КD и W – Pb > Cd > Cu > Co > Zn > Ni. Это связано с тем, что параметр J зависит от массы взятого в опыт карбоната.

Таблица 4. Кинетика поглощения цинка и свинца известняком (Т:Ж = 0.04, кг/кг Н2О)

t, час

0.25

0.5

1

2

4

6

8

10

t ® ¥

Zn, Со = 3.5 мг/л

Сt, мг/л

3.23

2.82

2.18

1.80

0.88

0.76

0.51

0.62

0.50

J, мг/г

6.7

17.0

33.0

42.5

65.5

68.5

74.7

72.0

75.0

KD, л/г

2.1

6.0

15.1

23.6

74.4

90.1

146.6

116.1

150.0

R, %

7.7

19.4

37.7

48.6

74.9

78.3

85.4

82.3

85.7

Pb, Со = 3.5 мг/л

Сt, мг/л

2.24

1.73

1.24

0.55

0.21

-

-

-

0.05

J, мг/г

31.5

44.2

56.5

73.7

82.2

-

-

-

86.2

KD, л/г

14.1

25.6

45.6

134.1

391.7

-

-

-

1725

R, %

36.0

50.6

64.6

84.3

94.0

-

-

-

98.6

Примечание. Знак «t ® ¥» означает, что данные получены путем экстраполяции кривой Ct=f(t).

 

Крайне важным в практическом отношении фактором при сорбционной очистке вод является кинетика процесса, т.е. скорость детоксикации водной среды. Влияние длительности сорбции цинка и свинца известняком на Сt, J, КD и W для случая Т:Ж = 0.04 приведены в табл. 4 и на рис. 2. Они показывают на неплохие кинетические свойства карбоната, однако, при t = const степень извлечения разных металлов различна. Если в течение 4 часов контакта раствора, содержащего 3.5 мг/л металла, с карбонатом степень извлечения свинца из водной фазы составляет 94%, то для цинка эта величина не превышает 75%. Поэтому в случае использования карбонатных пород для детоксикации многокомпонентных водных сред надо ориентироваться на скорость иммобилизации конкретного металла.

Рис. 2. Кинетика изменения концентрации цинка и свинца в растворе, находящемся в контакте с известняком при Т:Ж = 0.04.

Рис. 2. Кинетика изменения концентрации цинка и свинца в растворе, находящемся в контакте с известняком при Т:Ж = 0.04.

Таким образом, полученные данные показывают, что природные карбонатные породы Дагестана с успехом могут быть использованы для очистки вод от ТМ; степень извлечения Pb, Cd, Cu, Co и Zn превышает 90%. При этом известняк заметно эффективнее иммобилизует ТМ из водных растворов, чем доломит. Возможность же использования карбонатных пород для извлечения Ni из водных сред не столь привлекательна; степень его извлечения составляет всего 40%.

Обсуждение экспериментальных данных

Базовая теоретическая предпосылка использования карбонатов для иммобилизации ТМ из водных растворов основана [8] на образование карбоната МеСО3, произведение растворимости (ПРМеСО3) которого меньше, чем ПРСаСО3:

Ме2+ (р) + СаСО3 (тв) ®¬ МеСО3 (тв) + Са2+ (р), (1)

где индекс (р) означает раствор, (тв) – твердое вещество, Ме - металл.

Записав процесс растворения карбоната в виде МеCO3 (тв) ®¬ Ме2+ (р) + CO32- (р), учитывая, что ПРМеСО3 = [Ме2+][СО32-], и, полагая, что активности ионов Са2+ и Ме2+ равны их концентрациям, получим:

КПР = [Са2+]/[Ме2+] = ПРСаСО3/ПРМеСО3, (2)

где КПР – константа равновесия реакции (1), квадратными скобками обозначены моляльные концентрации.

Из (2) следует, что если численное значение КПР > 1, то равновесие реакции (2) сдвинуто в сторону образования МеСО3 (тв), если КПР < 1, то в сторону исходных веществ. Надежные выводы о преобладании в смеси продуктов реакции можно, однако, делать при КПР > 102.

В табл. 5 приведены справочные величины ПР карбонатов. Видно, что они различаются на порядок и более, и это затрудняет анализ экспериментальных данных. Для расчетов мы использовали усредненные величины, при определении которых отбрасывали резко отличающиеся (выделенные курсивом) значения.

Табл. 5. Произведение растворимости карбонатов при 25 оС

 

Карбонат

Образую-

щиеся ионы

ПРМеСО3 по данным

 

KПР

[8]

[9]

[10]

[11]

сред

CaCO3

Ca2+; CO32-

3.3×10-9

4.96×10-9

4.4×10-9

4.7×10-9

4.34×10-9

-

CdCO3

Cd2+; CO32-

1.0×10-12

6.18×10-12

2.5×10-14

5.2×10-12

4.13×10-12

1.05×103

CoCO3

Co2+; CO32-

1.4×10-10

1.40×10-13

1.5×10-10

-

1.45×10-10

2.99×101

NiCO3

Ni2+; CO32-

1.5×10-7

1.42×10-7

1.3×10-7

1.4×10-7

1.40×10-7

3.10×10-2

PbCO3

Pb2+; CO32-

3.6×10-14

1.46×10-13

3.6×10-14

1.5×10-15

3.60×10-14

1.21×105

ZnCO3

Zn2+; CO32-

5.0×10-11

1.19×10-10

5.3×10-11

2.0×10-10

5.15×10-11

8.43×101

CuCO3

Cu2+; CO32-

2.4×10-10

1.40×10-10

2.4×10-10

2.5×10-10

2.18×10-10

1.99×101

Судя по величинам КПР, Pb и Cd должны хорошо извлекаться из водных растворов карбонатами кальция, извлечение Zn, Co и Cu является проблематичным, а Ni вовсе не должен извлекаться. Эти положения на практике выполняются не однозначно.

Действительно, с одной стороны, наблюдается некоторая корреляция между степенью извлечения металла и произведением растворимости его карбоната (Pb2+ выводится из раствора значительно быстрее, чем, скажем, Zn2+, что можно связать с ПРPbCO3 < ПРZnCO3 < ПРСаСО3), но, с другой стороны, Ni2+ извлекается из раствора (табл. 3), несмотря на то, что ПРNiCO3 > ПРCaCO3. Главное же противоречие между экспериментом и теорией заключается в следующем. Если бы лимитирующей стадией процесса иммобилизации металла была бы только реакция (1), то Т:Ж отношение не должно было влиять на степень его иммобилизации из водного раствора, поскольку концентрация СО32- в растворе не зависит от количества взятого в опыт карбоната; чистая вода, находящаяся в контакте с СаСО3 (тв), содержит [CO32-] = ÖПРСаСО3 = Ö4.34ּ10-9 = 0.066 ммоль/кг Н2О. Экспериментальные данные (табл. 2 и 3) показывают обратное.

С целью более углубленного анализа системы СаСО3 (тв)-Н2О-Ме2+ (р) нами были привлечены положения равновесной химической термодинамики, которая позволяет учитывать участие твердых компонентов в суммарной реакции. Однако и термодинамика имеет ограничения; термодинамические выводы зависят, во-первых, от уровня наших знаний о химическом и фазовом составах исходных веществ и продуктов реакции и, во-вторых, от наличия и достоверности характеризующей их термодинамической информации.

Растворение карбоната кальция в воде можно представить не как прямой процесс диссоциации СаСО3 (тв) ®¬ Са2+ + СО32-, а как гидратацию структурных единиц СаСО3о и их последующую диссоциацию

СаСО3 (тв) ®¬ СаСО3о (р) ®¬ Са2+ (р) + СО32- (р). (3)

Косвенным доказательством реальности процесса (3) является обнаружение спектроскопическими методами в водной фазе ассоциированных частиц, соответствующие составу твердой фазы карбонатов, сульфатов, фторидов и других минералов-солей [8]. Тогда, процесс иммобилизации Ме2+ можно представить в виде:

Ме2+ (р) + СаСО3о (р) ®¬ МеСО3о (р) + Са2+ (р), (4)

{СаСО3} (тв)∙МеСО3 (тв)

Ясно, что в условиях отсутствия в системе СаСО3 (тв) переход МеСО3о (р) ® МеСО3 (тв) должен начинаться с образования собственных кристаллических зародышей, но при наличии в системе СаСО3 (тв) именно СаСО3 (тв) предоставляет центры кристаллизации для такого перехода, и образование МеСО3 (тв) протекает на поверхности СаСО3 (тв). Иными словами, роль СаСО3 (тв) в рассматриваемой системе сводится не только к предоставлению ионов СО32- в водную фазу для образования МеСО3о, но и своей поверхности (центров кристаллизации) для вывода МеСО3о (р) из зоны реакции, что, по-видимому, и определяет влияние Т:Ж отношения на иммобилизацию Ме2+ из водного раствора.

Для термодинамической проверки этих предположений реакцию (4) достаточно представить в виде:

Ме2+ (р) + СаСО3о (р) ®¬ МеСО3 (тв) + Са2+ (р). (5)

Теперь сопоставим термодинамические параметры реакций (1) и (5). Их вычисляли с помощью уравнений:

DHo реакции = åDHo продуктов реакции - åDHo исходных веществ; (6)

DSo реакции = åSo продуктов реакции - åSo исходных веществ; (7)

DGо реакции = åDGo продуктов реакции - åDGo исходных веществ; (8)

lgKo равновесия = - DGo реакции/2.303RT, (9)

где Ко – стандартная (термодинамическая) константа равновесия реакции. С целью проверки точности вычислений энергию Гиббса определяли также по уравнению DGo* реакции = DHo реакции - Т×DSo реакции.

Использованные в расчетах усредненные из табулированных в разных справочниках термодинамические параметры компонентов системы приведены в табл. 6, а рассчитанные значения величин термодинамических параметров реакций (1) и (5), – в табл. 7. Заметив предварительно, что различия между DGo и DGo*, а также между константами равновесия реакции (1), найденными из уравнений (2) и (9), находятся в пределах точности исходных термодинамических данных, проанализируем полученные результаты.

Таблица 6. Энтальпии образования (DHo298), энтропия (So298) и энергия Гиббса образования (DGo298) для ионов и карбонатов в водных растворах (по данным [8, 9, 10, 12, 13]).

Компоненты

раствора

DHo298,

кДж/моль

So298,

Дж/(К×моль)

DGo298,

кДж/моль

Ca2+ (р)

-543.0

-54.5

-553.2

CaCO3 (тв)

CaCO3о (р)

-1207.6

-1220.8

+91.2

-110.1

-1128.4

-1093.6

Cd2+ (р)

-74.2

-66.0

-77.8

CdCO3 (тв)

-754.5

+97.0

-674.0

Cu2+ (р)

+65.1

-98.1

+65.2

CuCO3 (тв)

-594.0

+92.0

-518.0

Co2+ (р)

-63.2

-110.4

-53.8

CoCO3 (тв)

-723.2

+88.4

-642.2

Ni2+ (р)

-59.0

-126.0

-54.7

NiCO3 (тв)

-691.5

+86.1

-612.0

Pb2+ (р)

-1.0

+11.0

-24.0

PbCO3 (тв)

-701.3

+131.0

-627.4

Zn2+ (р)

-153.5

-109.9

-147.1

ZnCO3 (тв)

-813.4

+82.2

-732.7



Примечание. В тексте и в формулах индекс «298» при DHo, DSo и DGо опущен.

 

Таблица 7. Термодинамические параметры реакций (1) и (5) при Т=298 К

 

Реакция

Параметры реакции

DHo298

DSo298

DGo298

DG*o298

lgKo

Ko

CaCO3 (тв) + Cd2+ (р) = CdCO3 (тв) + Ca2+ (р)

CaCO3о (р) + Cd2+ (р) = CdCO3 (тв) + Ca2+ (р)

-15.5

-2.5

+17.3

+218.6

-21.0

-55.8

-20.7

-62.6

3.68

9.78

4.79×103

6.03×109

CaCO3 (тв) + Co2+ (р) = CoCO3 (тв) + Ca2+ (р)

CaCO3о (р) + Co2+ (р) = CoCO3 (тв) + Ca2+ (р)

+4.6

+17.8

+53.1

+254.4

-13.2

-48.0

-11.2

-58.0

2.31

8.41

2.04×102

2.57×108

CaCO3 (тв) + Cu2+ (р) = CuCO3 (тв) + Ca2+ (р)

CaCO3о (р) + Cu2+ (р) = CuCO3 (тв) + Ca2+ (р)

+5.5

+18.7

+44.4

+245.7

-8.0

-42.8

-7.7

-54.5

1.40

7.50

2.51×101

3.16×107

CaCO3 (тв) + Ni2+ (р) = NiCO3 (тв) + Ca2+ (р)

CaCO3о (р) + Ni2+ (р) = NiCO3 (тв) + Ca2+ (р)

+32.1

+45.3

+66.4

+267.7

+17.9

-16.9

+12.3

-34.5

-2.15

2.96

1.41×10-2

9.12×102

CaCO3 (тв) + Pb2+ (р) = PbCO3 (тв) + Ca2+ (р)

CaCO3о (р) + Pb2+ (р) = PbCO3 (тв) + Ca2+ (р)

-35.7

-22.5

-25.7

+175.6

-28.4

-63.0

-28.1

-74.8

4.98

11.04

9.55×104

1.09×1011

CaCO3 (тв) + Zn2+ (р) = ZnCO3 (тв) + Ca2+ (р)

CaCO3о (р) + Zn2+ (р) = ZnCO3 (тв) + Ca2+ (р)

+4.7

+17.9

+46.4

+247.7

-10.4

-45.2

-9.1

-55.9

1.82

7.92

6.61×101

8.32×107

 

Анализ позволяет сформулировать два важных выводов: а) с термодинамической точки зрения, протекание реакции (5) значительно выгоднее, чем реакции (1); различие Ко достигает 6 порядков; б) извлечение никеля из водных растворов карбонатом кальция термодинамически противоречит реакции (1), но не противоречит реакции (5).

Термодинамика позволяет также предсказать влияние температуры на интенсивность иммобилизации ТМ из водной среды, что важно при прогнозировании влияния сезона года на качество природных вод. Для узкого интервала температуры температурную зависимость Ko можно оценить с помощью уравнения lgKoТ = (-DHo/RT + DSo/R)/2.303.

Рис. 3. Температурная зависимость константы равновесия реакции (5). Уравнения регрессии: lgKoCd = 11.4153 + 0.131/Т; lgKoCo = 13.2873 – 0.9288/Т; lgKoCu = 12.8354 – 0.9778/Т; lgKoPb = 9.1767 + 1.1733/Т; lgKoZn = 12.939 – 0.9357/Т; lgKoNi = 13.9785 – 2.3651/Т.

Рис. 3. Температурная зависимость константы равновесия реакции (5). Уравнения регрессии: lgKoCd = 11.4153 + 0.131/Т; lgKoCo = 13.2873 – 0.9288/Т; lgKoCu = 12.8354 – 0.9778/Т; lgKoPb = 9.1767 + 1.1733/Т; lgKoZn = 12.939 – 0.9357/Т; lgKoNi = 13.9785 – 2.3651/Т.

 

Как видно из рис. 3, повышение температуры способствует увеличению сорбции Ni, Co, Cu, Zn, уменьшению сорбции Pb и практически не влияет на сорбцию Cd карбонатными породами. В практическом отношении это означает, что с увеличением температуры различия между сорбциями исследованных металлов нивелируется. Так, ΔlgKо = lgKоPb - lgKоNi при Т= 0 оС составляет 7.2, при Т = 45 оС – 5.7. Отсюда следует, что влияние сезона на абсолютные величины содержания ТМ в природных водах, контактирующих с карбонатными породами, будет существенным.

Кинетически процесс иммобилизации Ме2+ из водных растворов карбонатами можно анализировать как со стороны раствора, так и со стороны твердой фазы. При этом, однако, нужно учесть, что единственным параметром системы, который может быть измерен точно, является концентрация металла в водном растворе, параметры же, характеризующие сорбент (J и KD), являются расчетными.

Несомненно, начальным актом процесса иммобилизации Ме2+ из водной фазы является гомогенная реакция Ме2+ (р) + СО32- (р) = МеСО3 (тв). С точки зрения формальной химической кинетики ее скорость по стехиометрии является бимолекулярной

- dC/dt = kCMe2+ CCO32-. (10)

Но поскольку концентрация СО32- в растворе, находящейся в равновесии с твердой фазой СаСО3, является величиной постоянной, то можно написать:

- dC/dt = k1CМе2+. (11)

где k1 = kCCO32-. Это уравнение предполагает, что СМе2+ ® 0 при t ® ¥. Однако эксперимент (см. табл. 4) показывает, что при t ® ¥, СМе2+ стремится к некоторой равновесной величине С¥. Поэтому перепишем (11) в виде

- dC/dt = k2(Ct-C¥), (12)

интегрирование, которого дает:

2.303lg[(Ct-Cр)/(Cо-C¥)] = k2t, (13)

где Со, С¥ и Ct – концентрация металла в растворе начальная, конечная (равновесная) и к моменту времени t.

Рис. 4. Выполнимость уравнения (13) для описания процесса иммобилизации цинка и свинца из водного раствора известняком. Уравнение регрессии для Zn: у=0.0045-0.2192×t (С¥ = 0.5 мг/л); для Pb: у=-0.0665+0.3948×t (С¥ = 0.05 мг/л).

Рис. 4. Выполнимость уравнения (13) для описания процесса иммобилизации цинка и свинца из водного раствора известняком. Уравнение регрессии для Zn: у=0.0045-0.2192×t (С¥ = 0.5 мг/л); для Pb: у=-0.0665+0.3948×t (С¥ = 0.05 мг/л).

Обработка экспериментальных данных в координатах уравнения (13), показала (рис. 4) его хорошую выполнимость; процесс иммобилизации ТМ из водных растворов известняком подчиняется уравнению реакции первого порядка относительно концентрации металла в растворе. Константа скорости иммобилизации цинка составляет 0.096 час-1, свинца = 0.174 час-1, т.е. свинец удаляется из раствора почти в 2 раза быстрее, чем цинк.

Хотя к сорбционным процессам не применим закон действующих масс, нетрудно показать, что и кинетика процесса адсорбции также протекает по первому порядку. (Ясно, что здесь речь идет о формально-кинетическом порядке). При справедливости допущения, что скорость адсорбции пропорциональна свободной площади поверхности сорбента (вернее, концентрации центров кристаллизации) можно записать dq/dt = kадс(1 - q), или в интегральной форме q = 1 – eхр(kадсt), где q – поверхность сорбента, kадс – константа скорости адсорбции. Проведя аналогию с другими силикатными материалами [14], логично было рассмотреть реакцию (4) как процесс самоторможения за счет образования пленки (корки) МеСО3 на поверхности СаСО3 (тв). Тогда скорость изменения массы взятого в опыт карбонатной породы в интегральной форме запишется в виде:

{2.303lg[100/(100 - J)]}/t = b(J/t), (14)

где b – коэффициент торможения (экранирования) процесса.

Экспериментальные данные в координатах уравнения (14) приведены на рис. 5. Отсечение некоторого участка на оси ординат свидетельствует, что процесс самоторможения начинается после осаждения на сорбенте 5-10 мг МеСО3 на 1 г СаСО3 (тв), а коэффициент торможения PbCO3 (тв) примерно в 1.3 раза выше, чем ZnCO3 (тв). Физический смысл коэффициента b заключается в том, что по мере протекания процесса зона реакции передвигается от СаСО3 (тв), и в конечном счете его скорость должна лимитироваться либо диффузией воды к поверхности СаСО3 (тв) ПН2О = b1оt), либо диффузией СО32- ионов от поверхности СаСО3 (тв) в объем раствора ПСО32- = b2оt), где b1 и b2 - коэффициенты массопередачи, равные отношению коэффициентов диффузии Н2О и СО32- в МеСО3 (тв) к толщине пленки МеСО3 (тв). Очевидно, что в природных условиях, когда карбонатные породы контактируются с водной средой веками, описание процесса (4) перейдет от кинетической задачи к задаче динамики.

Рис. 5. К проверке выполнимости уравнения (14) для описания кинетики изменения массы известняка за счет адсорбции МеСО3о. Уравнение регрессии для Zn: у = 0.045 + 0.0035×J/t; для Pb: у = 0.114 + 0.0044×J/t.Рис. 5. К проверке выполнимости уравнения (14) для описания кинетики изменения массы известняка за счет адсорбции МеСО3о. Уравнение регрессии для Zn: у = 0.045 + 0.0035×J/t; для Pb: у = 0.114 + 0.0044×J/t.

Рис. 5. К проверке выполнимости уравнения (14) для описания кинетики изменения массы известняка за счет адсорбции МеСО3о. Уравнение регрессии для Zn: у = 0.045 + 0.0035×J/t; для Pb: у = 0.114 + 0.0044×J/t.

 

Одним из основных критериев оценки адсорбционных свойств сорбента является изотерма адсорбции. Для описания изотерм сорбции привлекаются уравнения Лэнгмюра {J = zbC/(1+bC)}, Темкина {J = z/g·ln[(1+bo¢C)/(1+bo¢e-gC)]}, Фрумкина {J = zbo¢¢egJC/(1+bo¢¢egJC)}, Фрумкина-Шлыгина {J = A+z/g·lnC}, Фрейндлиха {J = K·C1/n} и другие. В этих уравнениях: g = a/RT, A = z/g∙ln bo¢, b = kадс/kдес=boeQ/RT, bo¢ = boeQm/RT, bo¢¢ = boeQo/RT, где Q, Qm и Qo – теплота адсорбции, максимальная теплота адсорбции и теплота взаимодействия адсорбированных молекул; z – число мест (центров кристаллизации) на адсорбенте, С – концентрация металла в растворе, a, bо, К, n –постоянные.

Изотермы сорбции цинка и свинца на известняке были линеаризованы в координатах уравнений Лэнгмюра, Фрумкина-Шлыгина и Фрейдлиха (рис. 6).

Рис. 6. Изотермы сорбции цинка и меди в координатах уравнений Лэнгмюра (C/J – C), Фрумкина-Шлыгина (J – lgC) и Фрейдлиха (lgJ – lgC). Уравнения регрессии: С/JZn = 0.0043+0.0141·C, R = 0.9956; C/JPb = 0.0026+0.0114·C, R = 0.9968; JZn = 47.11 + 30.3164·lgC, R = 0.9998; JPb = 63.83 + 30.4533·lgC, R = 0.9950; lgJZn = 1.6621+0.2527·lgC, R=0.9975; lgJPb = 1.7999+0.1895·lgC, R=0.9981.Рис. 6. Изотермы сорбции цинка и меди в координатах уравнений Лэнгмюра (C/J – C), Фрумкина-Шлыгина (J – lgC) и Фрейдлиха (lgJ – lgC). Уравнения регрессии: С/JZn = 0.0043+0.0141·C, R = 0.9956; C/JPb = 0.0026+0.0114·C, R = 0.9968; JZn = 47.11 + 30.3164·lgC, R = 0.9998; JPb = 63.83 + 30.4533·lgC, R = 0.9950; lgJZn = 1.6621+0.2527·lgC, R=0.9975; lgJPb = 1.7999+0.1895·lgC, R=0.9981.Рис. 6. Изотермы сорбции цинка и меди в координатах уравнений Лэнгмюра (C/J – C), Фрумкина-Шлыгина (J – lgC) и Фрейдлиха (lgJ – lgC). Уравнения регрессии: С/JZn = 0.0043+0.0141·C, R = 0.9956; C/JPb = 0.0026+0.0114·C, R = 0.9968; JZn = 47.11 + 30.3164·lgC, R = 0.9998; JPb = 63.83 + 30.4533·lgC, R = 0.9950; lgJZn = 1.6621+0.2527·lgC, R=0.9975; lgJPb = 1.7999+0.1895·lgC, R=0.9981.

Рис. 6. Изотермы сорбции цинка и меди в координатах уравнений Лэнгмюра (C/J – C), Фрумкина-Шлыгина (J – lgC) и Фрейдлиха (lgJ – lgC). Уравнения регрессии: С/JZn = 0.0043+0.0141·C, R = 0.9956; C/JPb = 0.0026+0.0114·C, R = 0.9968; JZn = 47.11 + 30.3164·lgC, R = 0.9998; JPb = 63.83 + 30.4533·lgC, R = 0.9950; lgJZn = 1.6621+0.2527·lgC, R=0.9975; lgJPb = 1.7999+0.1895·lgC, R=0.9981.

 

Как видно, эти модели, несмотря на разную теоретическую базу (адсорбция на однородной-неоднородной поверхности, линейное-нелинейное распределение активных центров поверхности по теплотам адсорбции, отсутствие-наличие взаимодействия между адсорбированными молекулами) не позволяют предпочесть какую-либо из них; критерий R практически одинаков и превышает 0.995. Скорее всего, это обусловлено небольшим (всего четыре точки на кинетическую кривую) объемом экспериментального материала. Но, тем не менее, некоторые выводы можно сделать. Например, из уравнения Лэнгмюра следует, что константы адсорбции цинка и свинца на известняке в (3-6)ּ105 раз больше констант их десорбции, а практически одинаковый угол наклона анаморфозы уравнения Фрумкина-Шлыгина указывает на независимость числа активных центров кристаллизации сорбента от природы сорбируемого металла.

Влияние природы иона Ме2+ на адсорбцию, по-видимому, проявляется через теплоту адсорбции, но экспериментально проверить это невозможно. Однако такое влияние косвенно можно оценить с помощью физико-химических свойств металла (табл. 8).

Таблица 8. Тип оболочки, радиусы (r), поляризуемость (α), ионный потенциал (IP), параметр Мисоно (Y) металлов

 

Ион

 

Тип оболочки

r·10-8 см

(к.ч. = 4)

r·10-8 см

(к.ч. = 6)

α·10-24 см3

(к.ч. = 6)

IP, нм-1

(к.ч. = 6)

 

Y, нм

Ni2+

d оболочки

0.055

0.69

0.30

29.0

0.252

Cu2+

d оболочки

0.057

0.73

0.15

27.4

0.284

Co2+

d оболочки

0.058

0.74

0.51

27.0

0.270

Zn2+

d10 оболочки

0.058

0.75

0.24

26.7

0.240

Cd2+

d10 оболочки

0.078

0.95

0.92

21.1

0.303

Pb2+

d и f оболочки

0.098

1.19

3.48

16.8

0.393



Ионный потенциал, т.е. поверхностный заряд иона (IP = Z/r), можно использовать для оценки степени «поверхностной диссоциации», и по этому признаку исследованные металлы располагаются в ряд: Pb > Cd > Zn » Co » Cu > Ni.

Параметр Мисоно (Y = 10IZr/(IZ+1Z1/2), где IZ - потенциал ионизации, Z – степень окисления металла, r – ионный радиус в нанометрах) представляет собой меру склонности катиона к поляризации и используется для оценки «жесткости» кислот Льюиса [15]. Жесткие кислоты характеризуются высокой плотностью заряда, высокой электроотрицательностью и низкой способностью к поляризации; они стремятся образовывать стабильные комплексы с большим увеличением энтропии и с электростатическим взаимодействием. Мягкие кислоты, напротив, характеризуются низкой плотностью заряда, малой электроотрицательностью и высокой способностью к поляризации; они образуют стабильные комплексы с уменьшением энтальпии и ковалентным взаимодействием. На самом деле (см. табл. 7), образование CdCO3 (тв) сопровождается уменьшением энтальпии, а образование NiCO3 (тв) и CoCO3 (тв) – большим увеличением энтропии. Если расположить исследование металлы по «степени жесткости», то получим ряд: Pb > Cd > Zn > Co > Cu > Ni. Можно предположить, что взаимодействие кислот Льюиса Ме2+ с основанием СО32- приведет к образованию комплексов МеСО3о, стабильность которых возрастает в том же ряду.

Теперь сравним экспериментальный ряд адсорбции Ме2+ известняком с рядами, прогнозируемыми термодинамическими параметрами и физико-химическими свойствами металлов. Такое сравнение для случая Т:Ж = 0.04 дано в табл. 9.

Таблица 9. Экспериментальные и прогнозируемые ряды сорбции Ме2+ известняком.

Параметр

оценки

 

Ряды адсорбции

Опыт

ПРМеСО3

DHo

DSo

DGо

IP

Y

Pb > Cd > Cu > Co > Zn > Ni

Pb > Cd > Zn > Co » Cu > Ni

Pb > Cd > Co » Zn > Cu > Ni

Pb > Cd > Cu » Zn > Co > Ni

Pb > Cd > Co > Zn » Cu > Ni

Pb > Cd > Zn » Co » Cu > Ni

Pb > Cd > Cu > Co > Ni > Zn



 

 

 

 

 

 

 

 

В целом имеет место довольно хорошее согласие прогнозируемых рядов адсорбции с опытным; степень согласия составляет. Конечно, можно достичь и лучшего согласия, например, путем подбора «нужных величин» ПР, термодинамических параметров, чисел координации. Но, строго говоря, вряд ли какой-то один параметр сможет одновременно учесть все многообразие явлений, протекающих в системе СаСО3 (тв) – Н2О – Ме2+. Теория процесса адсорбции Ме2+ карбонатными породами далеко нетривиальна.

 

Литература

1. Перельман А.И. Геохимия. М.: Высш. шк., 1989. 528 с.

2. Смирнов А.Д. Сорбционная очистка воды. Л.: Химия, 1982. 168с.

3. Лупейко Т.Г., Баян Е.М., Горбунова М.О. Использование техногенного карбонатсодержащего отхода для очистки водных растворов от ионов никеля (II) // Журнал прикладной химии, 2004. Т. 77, Вып. 1. С. 87-92.

4. Тарасевич Ю.А. Природные сорбенты в процессах очистки воды. Киев, 1981. 231 с.

5. Никифоров А.Ю., Ильина Л.А., Сударушкин А.Т. Использование природного минерала доломита и его термомодифицированных форм для очистки сточных вод от катионов тяжелых металлов // Изв. вузов. Химия и химическая технология. 1999. Том 42, № 4. С. 138-141.

6. Годымчук А.Ю., Решетова А.А. Исследование процессов извлечения тяжелых металлов на природных минералах // Электронный научно-информационный журнал "Вестник Отделения наук о Земле РАН". 2003. № 1.

7. Бутаев А.М., Гуруев М.А., Магомедбеков У.Г., Осипова Н.Ф., Магомедрасулова Х.М., Магомедова А.Д., Мухучев А.А. Тяжелые металлы в речных водах Дагестана // Вестник ДНЦ РАН, 2006, № 26. С. 43-50/

8. Крайнов С.Р., Рыженко Б.Н., Швец А.М. Геохимия подземных вод. Теоретические, прикладные и экологические аспекты. М.: Наука, 2004. 677 с.

9. Гороновский И.Т., Назаренко Ю.П., Некряч Е.Ф. Краткий справочник по химии. Киев: Наукова думка, 1987. 695 с.

10. Лидин Р.А., Андреева Л.Л., Молочко В.А. Константы неорганических веществ: справочник. М.: Дрофа, 2006. 685 с.

11. Дикерсон Р., Грей Г., Дж. Хейт. Основные законы химии. Т.1. М.: Мир, 1982. 652 с.

12. Новиков Г.И. Основы общей химии. М.: Высш. шк., 1988. 431 с.

13. http://www.chem.msu.ru

14. Бутаев А.М. Прочность стекла. Ионообменное упрочнение. Махачкала, 1997. 253 с.

15. Спозито Г. Распределение потенциально опасных следовых металлов // Некоторые вопросы токсичности ионов металлов. Под ред. Х.Зигель, А.Зигель. М.: Мир, 1993. С. 9-24.

 
Яндекс.Метрика
Яндекс цитирования Valid XHTML 1.0 Transitional