CC BY
21
Халюты, в особенности окрестности целебного источника Аршан-
%алк>тинского, являются традиционным местом отдыха и лечения
жителей Бурятии.
Литература
], Кабата-Пендиас А., Пендиас X. Макроэлементы в почвах и растениях. - М.: Мир, 1989 - 489 с.
2. Колотов Б.А. Гидрогеохимия рудньгх месторождений. - М.: Недра, 1992,- 193 с.
3. Критерии оценки экологической обстановки территорий для выявления зон чрезвычайной экологической ситуации и зон экологического бедствия. - М.: Минприрода РФ, 1992,- 58 с.
4. Сает Ю.А., Алексинская Л.Н., Янин Ь.П. Методические рекомендации по геохимической оценке загрязнения поверхнослных водотоков химическими элементами. - М.: Изд. ИМГРЭ, 1982,73 с.
5. Сает Ю.Е., Онищенко Т.Л., Янин Е.П. Методические рекомендации но геохимическим исследованиям для оценки воздействия на окружающую среду проектируемых горнодобывающих предприятий. - М.: Изд. ИМГРЭ, 1986.- 99 с.
6. Сает Ю.Е,, Ревич Б.А., Янин Е.П. и др. Геохимия окружающей соеды. - М.: Недра, 1990 - 335 с.
7. Ripp G.S., Khodanovich P.Yu., Smirnova О С. Nature of amagmatic carbonate and carbonate-bearing rocks of Transbaikalia // Геодинамика и эволюция Земли. - Новосибирск: Изд-во СО РАН НИЦ ОИГГМ, 1996-С. 162-163.
А.М. Плюсник, В.И. ГУНИН ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ФОРМИРОВАНИЯ ЗОНЫ ОКИСЛЕНИЯ СУЛЬФИДНЫХ МЕСТОРОЖДЕНИЙ
Введение
Окисление сульфидов это один из процессов, который окапает существенное воздействие на миграцию и концентрирование Ряда элементов в экзогенных условиях. Этому вопросу посвящено большое количество работ"(Альбов. 1960; Смирнов, 1963; Крейтер и , 1958). Накоплен богатый фактический материал по структуре, Минеральному составу зон окисления. Меньше внимания уделялось Исследованию динамики трансформации сульфидной минерализации в окислительных условиях. Полученные нами эксперименталь-
ные данные механизма окисления некоторых сульфидных минера лов позволяют восполнить этот пробел. При окислении сульфидои устанавливаются разнообразные окислительно-восстановительные и кислотно-основные условия среды, которые определяют формы на. хождения многих элементов, их миграционную способность (Плкх нин и др., 1988; Плюснин и др., 1990).
Задачей данного исследования являлись изучение методом математического моделирования трансформации различной суль фидной минерализации в экзогенных условиях, вынос и отложение вещества в рудном теле, формирование качественного и количественного состава ореолов и потоков рассеяния вокруг рудных тел.
Исходные данные В основу модели положены данные экспериментального изучения окисления мономинеральных фракций сульфидов. В результате проведения этих работ предложен химический механизм окисления сульфидов через образование ряда метастабильных форм серы, железа, мышьяка и др. элементов. По известной мегодике Панченков, Лебедев, 1985) определены порядки и константы скоростей этих реакций. Используемые уравнения представлены в табл. 1,2.
Таблица 1
Урв-ннмния реакшм окислакия <зульфи до о
f Минерал Уравинияя реакций Коист. 1 скорости 1 1/ЧВС 1
(Плрит О. ЗвРвЯ? ♦ о. Э6ЕО**- о* - 1вн*о - о. 10Рв"' * а гвгыоН) -+ о. оэз4ов * О. 1вН* 1. 3-Ю" 1
О. * О! * О-вНжО - 0. 14Рв** ♦ О. 1^30 4* "*■ О.-ОН"* + О-14Ри(0Н>а 1 Ю~ 5 1
О . 25ГГ«5в ч- 0. 3*?0» ■+■ О 75Н*0 — О 5Г«<ОИ], + 0,43 1.-10'в !
1 пирро-1 тин О. »Рей «■ О О.гГе" -• Й ■+ О.вН^О - 0.4РеГ0Ш> • О. 2Ре(ОИ)а * О . в'* 7. 1* |
Ке» О* ■» О.5Н,С 0.5Ре<0Н>* + О.йРв** О.63«0лв 1-Ю * |
1 Халъко-1 пирит О.АСирвЗ* 0 • 0, -*■ о.вн»о - о.*Гв<он),3 ► о. гстою* ■« О. 5 . 3*Ю * I
а. 5СиРезг +• О« ■+ 0.5Н,0 — о. 5РеЗв •» О.еСи<0И)е 3.0«10~л1
{ Галенит ГХЗ& ■*■ оя -* н.яа о 5?Ь(ОН)в * 0_5Ръв- ♦ 0.5Ь>Ж0ЯЯ* 1
еарьв ♦ о.ггз^Ов* о£ » о вец.о — О вврЬ<ОЙ)* + 0.22РВ" + п о-ю а|
^ Арсвно- 0 44»Аа£ п. ввн* ♦ О« -» Г1.4*ИтО - о 44КвАдП4' 6. О-1 О"а|
О а^ГеАвЯ о-в7Ц«Ав* •* о. азо, -+- нго - о.в7Г«к ♦ о.е7н* ■» 1 . С»10~*1
О, 57реЛаЗ Оя ♦ О 2вН,0 - о. ■* О.57Н* * » 0. 67ГеАяО- я. п-ю'т|
Таблица 2
Уравнения реахиий преопрАэовании продуктов окисления сульфидов
1 1 » Урввнения реакций |Коьготалта. | I гисороСТИ. | 1 \
1РЬ(0Н)в + н* ( ......1 ПО7 |
)Си(ОНЭ щ + 1 >г ибо« - си*" йо**- ■*- гк*о 1 ( 1 10* 1
|Гв(ОН)и ♦ И* * ИЗО« - - Ре** » &Оля- £КвО 1 1 [ 10т (
)Ки<ОП>» ► 0 МГ * О. — (-в 1 ОН) щ* ♦ Г>. " НвО 1 Ют I
! Р е ( ОН > с + (он > " рв{онл Э1 1 1 |1.вЙхЮв |
|РС:<СН}- ИЗО — Рл'* +■ ЗО,*- * КаО 1 1 ^ ю* (
1Гви* ■+ о2 V н рп - о.вбРесон 1я' * о.вен* 1 1 1 Юв 1
1 О. и - 7 5С р Н.о - + НЗП, Г 1 »0.017 |
0. 50 х -*■ О-ЙНаО О . 53 «Ое' ~ ♦ ОН"" 1 1 10.002 1
I О. БРвАйО« Н.О - О.5РеГ* -«- 0.5НгАзОЛ- ► ОН" ' 1 11.0*10"* |
1Н.аяо4- И* ♦ Рв(ОН11а - РвАчО-'2НвО Н*0 1 1 li.o-i.o-4 |
1 о ЗЭГвАмО О. Э '-1С ш « ЙяО - О. ЭЗКаАаСи-2НаО -+• О. ЙбОЯ- I 1. о»ю-4 \
1 |нво ^ н он ~ ¡К„-Ю ** |
1 И* | к:,-! . 2К10"1
Взаимодействие сульфида с окислителем, как правило, описывается несколькими параллельными реакциями, на протекание которых требуется различное количество кислорода. В модели используются, в основном, прямые реакции. Из равновесных применяются только реакции диссоциации воды и серной кислоты.
Перемещение продуктов реакции в растворе осуществляется но законам фильтрации с дисперсией его состава в зависимости от пористости среды.
Описание модели Для описания процесса была взята система дифференциальных уравнений материального баланса массы растворенных веществ в элементе водоносного пласта (Бочевер и др, 1979). Система уравнений переноса и химической кинетики записывается так:
аС/с/1\ V* с}<^с(х О * с?С/<1х- £ аи*(1р1/Ж- 81(1),
J
Где:
С, - Молярная концентрация элементов в поровом водном растворе й моль/л, г = 1 1;
Pj - Молярная плотность минералов сульфидов, гидроокислов и солей в единице обьема в моль/л, j — I -ь 8;
a ij - Сгехиометрический коэффициент для i-ro элемента в /-ой реакции окисления сульфидов или осаждения-растворения гидроокислов, солей.
V - Истинная скорость фильтрации грунтовых вод; D - Коэффициент фильтрационной продольной дисперсии D-V* /, где / характерный размер частиц породы. V-истинная скорость фильтрации;
S, - Прибыль или убыль 1-го элемента в растворе за счет равновесных и окислительно-восстановительных реакций продуктов окисления сульфидов, г - 1 -ь 13.
0
dp,/dt - { S,, X, 8„„ Яя,
+ (Y, кРш * Со2 - „* р„* Y\C, ),
i m I
где, к,, - константа скорости i-й реакции окисления s-ixj сульфида. dij - сгехиометрический коэффициент для s-ro сульфида в i-й реакции. s=l 4-4;
Á ¡ - стехиометрический коэффициент для кислорода 02 в г-й реакции;
кт„- константа скорости ш-й реакции п-го гидроокисла(соли); S„„ - стехиометрический коэффициент для n-го гидроокисла в ш-й реакции. п= i -г 4;
Лт i - стехиометрический коэффициент 1-го элемента в m-ой реакции, 1 < /< 13
0
Si = { Л г! к ,q
±2>*ПС'-ПС. >•
г 1 q
где, кг константа равновесия r-ой равновесной или константа скорости r-ой окислительно-восстановительной реакции. Я Г г стехиометрический коэффициент для 1-го элемента в r-ой реакции. 1 < 13;
, стехиометрический коэффициент для ц-го элемента в г-ой реакции. 1 < й< 13.
Для решения данной системы уравнений задаются начальные и граничные условия.
Начальные условия: С,(х, 0) *» 0, для всех I, /=1,13
р (х). для х в рудном теле
рАх,0)= {
0 для х в пустых породах Молярная плотность минералов сульфидов ръ (*) задается в частях от единицы. Л' - 1 -г 4
Граничные условия:
0, для всех С, кроме 02
С, (0,0= {
Со для О2 ОД, I) = 0,
где С0- концентрация кислорода в грунтовых водах.
Система уравнений решалась численно по явно-неявной схеме с использованием расщепления по конвективной, дисперсионной и кинетической составляющим (Самарский, 1977). В каждый момент времени для всей системы расчитывалась конвективная составляющая по явной схеме, затем дисперсионная по неявной схе-ме(метод прогонки) и в заключение считалась кинетика тоже явно. Расчетная область разбивалась сеткой с равномерным шагом по времени I, и по пространству. Шаг на сетке выбирался так, чтобы скорость фильтрации за выбранный шаг по времени не превышала величину шага по пространственной сетке.
V * t < Л
Количество твердого вещества также как и растворенное рассчитывалось в моль/л. Пористость в точке определялась по ф0р муле:
где: е-пористость в точке, Wt- мольный обьем j-го минерала л/моль, р, - молярная плотность /-го минерала.
Скорость фильтрации по всей расчетной области принималась постоянной и не зависела от изменения пористости. Пористость учитывалась только при образовании твердых осадков реакций, осаждение минералов прекращалось, если пористость достигала заданного минимума.
На основе данной модели разработан пакет программ для ПК типа ШМ в среде MS-DOS на языке FORTRAN, с возможностью представления результатов расчета в текстовом и графическом видах. В программе имеется возможность варьировать составом и соотношением исходных сульфидов, концентрацией кислорода в растворе, интенсивностью водообмена, пористостью вмещающей среды.
Результаты и обсуждение С помощью модели были просчитаны различные варианты возможного окисления сульфидов, результаты и анализ которых приводится ниже.
Проанализируем вначале окисление рудных тел, состоящих из мономинералов сульфидов. На рисунке 1 представлены через одинаковые промежутки времени от начала окисления рудные тела, первоначально состоящие на 80% из пирротина (а,б) и пирита (iy )-оставшиеся 20% занимают вмещающие породы и пористость. При окислении пирротина около 30 объемных процента вещества выносится, оставшиеся 70% заполняет поры вмещающих пород в форме гидроокислов железа. Такая картина наблюдается на всех горизонтах рудного тела. При окислении пирита гидроокислы железа накапливаются только в верхней части рудного тела, где они забивают, практически все пустоты. В нижней части колонки рудного тела их нет, вследствие чего значительно возрастает пористость окисленны* рудных тел. Различие в трансформации рудных тел, первоначально
Рис. ]. Распределение минералов в колонке рудных тел при компьютерном моделировании окисления мономинералов сульфидов. Ус-овные обозначения: 1-пирротин, 2-пирит, 3-гидрооккслы желе-за(-З), 4-халькопирит, 5-гидроокислы меди (+2), 6-вмещающая порода, 7-пористость пород, 8-направление фильтрации кислородона-сыщенных вод. Интенсивность водообмена 1 м/мес, начальное содержание кислорода 0.1 моль/л. Примечание: на рис. а.,с.,е. показаны рудные тела через 200 лет после начала окисления.
На рисунке 2 показаны изменения некоторых параметров Растворов, фильтрующихся через эти рудные тела. При окислении пирита устанавливается кислая среда, но сильная -^рная кислота образуется в растворе через образование слабых по-■¡итионовых кислот, которые получаются непосредственно из сульфидов. Поэтому пик максимума кислотности всегда смещен по набавлению фильтрации от горизонта окисляющихся сульфидов. В ^ультате этого на верхних горизонтах рудных тел накапливается, °сительно, большее количество гидроокислов железа, нижние '^зонты подвержены"более интенсивном}' выщелачиванию. Отме-различия в формах миграции железа, так при кислении пир-
С,т/1 С,т/1
Рис.2. Концентрация растворимых веществ в потоках рассеяния при окислении пирита (1), пирротина (2), халькопирита (3). Интенсивность водообмена 1 м/мес, концентрация кислорода 0.1 моль/л. Условные обозначения: 4-зона окисления, 5-направление фильтрации вод.
ротина в растворе находятся, преимущественно, его гидролизован-ные формы, тогда как при окислении пирита < 1 значительно содержание ионной формы. Вынос серы при окислснии пирита, в основном, происходиг в форме ШО/", БО/ , при окислении пирротина доминируют ее политионатные формы. Окисление рудного тела, состоящего первоначально только из халькопирита, происходит с образованием в твердой фазе гидроокислов железа и меди (рис. 1 Д. е).
Объемное соотношение между ними составляет 1:0.6. В природных условиях в растворах, омывающих рудные тела, всегда пр^ сутствует в той или иной форме диоксид углерода, поэтому гидр» окислы меди преобразуются в малахит, азурит и другие карбонат содержащие минералы.
В целом наблюдается более интенсивный вынос меди "а нижних горизонтах, в результате этого в верхней части рудных те4 образуется обогащенный медью горизонт окисленных руд, а вн№? образуются пустоты выщелачивания. Окисленные рудные тела пр'1'
обретают грибообразную форму. В зоне окисления рудного тела образуется гипергенный пирит. Сера выносится, преимущественно, в форме политионатов. На нижних горизонтах устанавливается кислая реда (рН=3-4).
Лрсенопиритовое рудное тело окисляется с значительным вмноарм вещества. Мышьяк мигрирует в форме РеАьО« , Н^АбОд , ¿ера в форме политионатов, гидросульфатов. Устанавливается кислая сре1а. Из продуктов окисления формируются отложения скородита. Так как его образование проходит через несколько метаста-бильных форм, он рассеивается по зоне окисления.
При окислении галенитового и сфалеритового рудных тел твердых продуктов окисления не отлагается, все вещество выносится. Сера - в форме политионатов, шосульфатов, свинец - в форме Ре , цинк - в форме В природных условиях вынос свинца и цинка контролируется образованием церуссита и смитсонита.
Большое влияние на изменение рудных тел при окислении оказывает интенсивность водообмена. При увеличении интенсивности водообмена, во-первых, к рудному телу поступает большее количество кислорода и, во-вторых, происходит более интенсивное рассеивание растворенных веществ. Рассмотрим окисление рудного тела, состоящего из халькопирита при различной интенсивности водообмена. При водообмене 1 м/сут зона окислительного преобразования за 8 лет захватывает 5 метров (рис. 3 а.). На верхних горизонтах присутствуют и сульфид и продукты его окисления.
При интенсивности водообмена 1 м/мес. зона окисления за 160 лет распространяется всего на 3 метра, граница раздела "исходный сульфид - продукты окисления" - резкая (рис.3 6Л В условиях недостатка кислорода начинают играть заметную роль относительно медленные реакции, для протекания которых требуется меньшее количество кислорода. Фиксируется отложение гипергенного пири-IЧ чего в условиях интенсивного водообмена не наблюдается. В ус I *свиях интенсивного водообмена кислорода достаточно для окисле-чЧч серы до сульфатов, эта форма ее миграции доминирует в ореоле Иссеяния. При замедленном водообмене основной формой мигра-'4п< серы являются политионаты.
в. 4
20 40 60 80%
К
ГШ 2Ш зИ 4и бСК
Рис.3. Распределение минералов в зоне окисления халькопиритового рудного тела при различной интенсивности водообмена. Условные обозначения: 1-халькопирит, 2-гидроокислы железа, 3-гидроокислы меди, 4-пирит, 5-вмещающая порода, 6-пористость пород, 7-направление фильтрации вод. Концентрация кислорода в водах 0.1 моль/л. Интенсивность водообмена: а.-1м/сут, Ь.-1м/мес.
При окислении пирита в условиях замедленного водообмена наблюдается отложение элементарной серы, гидроокислов железа отлагается значительно больше, устанавливается менее кислая среда. Ореол рассеивания продуктов окисления существенно сокращается.
Так как при окислении сульфидов создаются разнообразные условия среды, то при одновременном окислении нескольких суль фидов они влияют друг на друга. Это проявляется в интенсивное растворении одних веществ и высаживании других.
Рассмотрим окисление рудного тела, состоящего из пирита1 галенита. Соотношение между сульфидами 1:1. Проследим прои^ окисления в динамике. Константа скорости окисления у галеШ*1* значительно больше чем у пирита, поэтому первоначально галеон1 задает среду, в которой протекают реакции в растворе. ПроисхоД осаждение гидроокислов железа на всем пространстве зоны оки^1 ния (рис. 4 а. ). В результате реакции сульфат-иона с ионом сви^' образуется англезит. По мере расходования галенита на верхин* ^ ризонтах определяющую роль в формировании среды миграи^ чинает играть окисление пирита (4 б.). В результате чего созД«11
РисАРаспределение минералов в колонке при компьютерном моделировании окисления рудного тела состоящего первоначально из пирита и галенита. Условные обозначения: 1-пирит, 2-галенит, 3-гидроокисль! железа, 4-англезит, 5-вмещающая порода, 6-иористость, 7-напрааление фильтрации вод.
кислая среда, гидроокислы железа начинают интенсивно растворяться, и железо выносится из зоны окисления. Так как на верхних горизонтах рудного тела пирит окислился в какой-то мере одновременно с галенитом, здесь гидроокислы железа сохраняются. В средней части они выщелачиваются, практически полностью. Наблюдается разрыв в сплошности их распространения (4 в ).
В конечном итоге при полном окислении этого рудного тела гидроокислы железа остаются только на верхних горизонтах. Отложение англезита происходит на всех горизонтах рудного тела. В природных условиях происходит его постепенная замена церусси-том. Окисление рудного тела, состоящего из пирита и пирротина, Протекает аналогично. Первоначально условия протекания процесса Сдаются пирротином. Это определяет интенсивное отложение гидРоокислов железа на верхних горизонтах. Затем, по мере расходования пирротина, созд тся более кислая среда, и железо начинает выноситься. В результате этого формируется "шляпообразная" форма Рудных тел.
Иная картина наблюдается при окислении рудного тема. Со стоящего из галенита и халькопирита. Происходит интенсивн окисление галенита, но отложение англезита не столь интенсивно как в рассмотренной ранее системе, так как содержание с)л(), иона в растворе не значительно. Гидроокислы железа и меди отла гаются на всем протяжении рудного тела, на нижних горизонтах от лагается гипергенный пирит. В водном ореоле рассеяния такое руд ное тело проявляется высокими содержаниями свинца и менее вы. сокими меди (при условии отсутствия карбонатов), из форм мигра. ции серы преобладают политионаты.
Сопоставим результаты моделирования с данными изучения зон окисления реальных природных объектов. Для примера рас, смотрим два проявления рудной минерализации, расположенных в Западном Забайкалье в непосредственной близости друг от друга. Время образования оруденения оценивается как палеозойское. Минерализованные зоны кварц-сульфидного прожилкования локагш-!уются в поле гранитизированных эффузивов. Густота прожилкования и количество сульфидов в прожилках изменчивы по просгира нию и падению тел. На первом рудопроявлении рудные тела эндогенной минерализации, в среднем, на 30% состоят из пирита. 5% занимает халькопирит, остальной объем занимает кварц. На втором рудопроявлении эндогенная минерализация на 60% состоит из пирита, 20% занимает пирротин и 15% - халькопирит. На рудопрояв-лениях существует горизонт окисленных руд. На первом рудопроявлении зона окисления местами прослеживается до глубины 100 метров, но железная "шляпа" располагается только в ее верхней части, ее мощность составляет 26 метров. На втором рудопроявлении железная "шляпа" фиксируется до 80 м. В зоне окисления отмечаются гипергенные отложения халькозина, ковеллина, пирита. Проведенные расчеты окисления руд этих рудопроявлений по нашей модели показали, что при окислении рудных тел первого рудопроявления среда миграции элементов задается окислением пирита (рис,5..а) Устанавливается очень кислая среда, происходит интенсивный вынос вещества растворами.
На месте расположения сульфидов формируются пустоты выщелачивания. Отложение гидроокислов железа происходит толь' ко на верхних горизонтах зоны окисления. За 250 лет при интенсив ности водообмена 1 м/сут и содержании кислорода 0.01 моль/л пр'1 изошло отложение гидроокислов железа в пределах двух метров
втором рудопроявлснии устанавливается менее кислая среда, что благоприятствует отложению гидроокислов.
Рис.5. Результаты моделирования рудных гел Замогтинского (а.) и Зангодинского рудопроявлений. Условные обозначения на рис.1. Интенсивность водообмена 1 м/сут, содержание кислорода 0.01 моль/л, продолжительность эксперимента 250 лет.
Пустоты выщелачивания заполняются гипергенными минералами фис.5., б.). За этот же промежуток времени при той же интенсивности водообмена и концентрации кислорода произошло окисление на 10 метров. В верхней части рудного тела фиксируется горизонт обогащенный гипер1енными минералами меди, здесь же обнаруживается гипергенный пирит.
Выводы
Предложена численная математическая модель процесса преобразования рудной минерализации, основанная на кинетике окисления сульфидов и перемещения растворов по законам фильтрации с дисперсией состава.
Путем моделирования установлено следующее.
1. Различия в механизме окисления сульфидов проявляются Не только в химическом составе ореолов рассеяния, но и в особенностях трансформации рудных тел в экзогенных условиях.
2. Прйлодостатке кислорода доминируют реакции, на протекание которых требуется минимальное количество окислителя. В
результате этого в экзогенных условиях могут отлагаться пирит, элементарная сера.
3. На преобразование рудных тел большое влияние оказывает интенсивность водообмена. При замедленном водообмене формируется резкая граница раздела "иервичныс*сульфиды - зона окисления", при интенсивном - размытая. При изменении интенсивности водообмена меняется соотношение между отлагающимися гипергенными минералами.
4. В зоне окисления рудных тел, состоящих из нескольких сульфидов, окислительно - востановительные и кислотно - основные условия среды меняются во времени, происходят то отложение, то вынос гипергенных минералов из зоны окисления, из-за этого окисленные рудные тела продолжают трансформироваться во времени. Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований. Код проекта 96-05-65888.
Литература
1. Альбов М.Н. Вторичная зональность золоторудных месторождений Урала. М.: Госгеолтехиздат.1960.- 215 с.
2. 2. Смирнов С.С. Зона окисления сульфидных месторожде-ний.М.-Л., ОНТИ.- 1963.- 292 с.
3. Крейтер В.М., Аристов В.В., Волынский И.С.,Крестовников А.Н., Кувичинский В.В. Поведение золота в зоне окисления вд-лотосульфидных месторождений.М., 1958,- С.226-236.
4. Плюснин A.M., Миронов А.Г., Константинова И.М., Фомина Г.П. Кинегика окисления сульфидов и растворение золота, ДАН СССР, 1988.- Т.302, 1,-С.192-195.
5. Плюснин A.M., Миронов А.Г., Беломестнова Н.В., Чернигова С.Е. Экспериментальное исследование окисления золотосодержащих сульфидов. Геохимия, 1, 1990-С.51-60.
6. Панченков Г.М., Лебедев В.Г1. Химическая кинетика и катализ М., Химия, 1985,- 592 с.
7. Бочевер Ф.М., Лапшин H.H., Орадовская А.Е. Защита подземных вод от загрязнения. М., Недра, 1979.- 253с.
8. Самарский A.A. Теория разностных схем. М., Наука, 1977.- 653с
