Физико-химическая модель образования гипергенных сульфатов на примере Шерловогорского месторождения (Восточное Забайкалье)
Приводятся результаты моделирования процессов гипергенного сульфатообразования в карьере Шерловогорско-го месторождения в диапазоне температур 0-45°С с помощью программного комплекса «Селектор». Установлены разнообразие сульфатов, очередность их появления и поля устойчивости, составлены ряды сульфатов в зависимости от температуры и количества молей воды в моделируемой системе.
Physico-chemical model of formation of Hypergene sulfates on example Sherlovogorsk Deposit (Eastern Transbaikalia).pdf Введение Шерловогорское олово-полиметаллическое месторождение разрабатывалось с 1936 по 1993 г. Итогом деятельности одноименного ГОКа стали оловоруд-ный карьер с техногенным озером, образовавшийся в результате затопления карьера, хвостохранилище, сеть отвалов и складов вскрышных пород и некондиционных руд. Основные рудные минералы месторождения представлены касситеритом, сфалеритом, галенитом, пиритом, пирротином, халькопиритом, арсе-нопиритом; нерудные - кварцем, хлоритом, турмалином, полевыми шпатами, кальцитом, флюоритом и др. [Касаткин, Клопотов, Плашил, 2014]. Гипергенные геохимические процессы в карьере Шерловогорского месторождения развиваются достаточно быстро и охватывают значительную часть территории. Процесс окисления руд и вмещающих пород приводит к появлению большого числа разнообразных минералов групп сульфатов, карбонатов, оксидов, гидроксидов и др. Сульфаты представлены растворимыми и нерастворимыми формами. Растворимые в основном представлены кристаллогидратами с формулой Me2+SO4nH2O. Целью работы являются установление минералов класса сульфатов, образующихся на испарительном геохимическом барьере в температурном диапазоне 0-45°С, условия их образования и стадийность. Исследованию гипергенных процессов уделялось внимание многих ученых на протяжении длительного периода времени [Карасик, 1946; Яхонтова, 1961; Jambor, Boyle, 1962; Nordstrom, 1982; Авдонин, 1984; Двуреченская, 1987; Жданов, 1989; Щербакова, 1995; Юргенсон, 1997; Jambor, Nordstrom, Alpers, 2000; Бортникова, 2001; Белогуб, Щербакова, Ни-кандрова, 2007; Сергутская, 2013; Зверева, Ястрем-ская, Лысенко, 2014; Русаль, 2015 и др.]. Большой интерес в настоящее время представляют процессы гипергенеза, протекающие при отрицательных температурах [Юргенсон, 1997; Еремин, 2004; Павлюкова, © Русаль О.С., 2018 DOI: 10.17223/25421379/8/5 Маркович, 2006; Marion, Kargel, Catling, 2008; Marion, Mironenko, Roberts, 2010; Птицын и др., 2009; Эпова, 2014; Moncur et al., 2015]. Гипергенные сульфаты эфемерные, условия их возникновения и пределы устойчивости достаточно ограничены, поэтому натурные исследования представляются в большей части нереальными. Метод компьютерного моделирования, предложенный И. К. Карповым и реализованный в компьютерную программу [Карпов, 1971; Карпов, 1977; Карпов, 1981; Карпов и др., 1999; Бычинский, Исаев, Тупицин, 2004; Авченко, Чудненко, Александров, 2009; Чудненко, 2010], основан на принципе работы, заключающимся в минимизации энергии Гиббса и построении термодинамической модели, содержащей информацию о качественном и количественном составе раствора и процентном соотношении выпавших из него минералов, pH-Eh-параметрах системы. Объект и методы Объектом исследования выступают гипергенные минералы класса сульфатов, образующиеся на испарительном геохимическом барьере в карьере Шер-ловогорского месторождения. Физико-химическое моделирование процессов образования сульфатов осуществлялось с помощью программного комплекса «Селектор» [Чудненко, 2010]. Расчет проводился для системы 24 независимых компонентов: Mg, Zn, Fe, Mn, Cu, Ni, Co, Ca, S, Si, Al, Pb, C, Cl, F, As, N, P, Sn, K, Na, H, O, e. Учитывались 323 зависимых компонента, в том числе 147 компонентов водного раствора, 19 газов и 157 твердых фаз. Использованы термодинамические базы данных для водных, газовых компонентов и для твердых фаз. Список твердых фаз дополнен согласованными значениями изобарно-изотермического потенциала некоторых возможных минералов сульфатов с кристаллизационной водой, их энтальпией образования из элементов и энтропией в стандартном состоянии (298,15 K; 1 бар). Модель однорезервуарная, открытая по отношению к атмосфере. Расчет химического равновесного состава проводился при температурах от 0 до +45°С и давлении 1 атм. Рассматриваются процессы, происходящие на испарительном барьере в карьере месторождения на контакте порода - поверхностные воды в соотношении 10:1, состав атмосферы C - 0,0104977, N - 54,65, O -14,61099. Минимальное значение выводимого зависимого компонента 10-8. Для расчета использовалось среднее значение состава руд и вмещающих пород, отобранных в карьере месторождения. Результаты исследования и обсуждение Моделируемая система имитирует процесс испарения воды на испарительном геохимическом барьере. В рассматриваемой модели уделяется внимание минералам магния, меди, цинка. При 0°С первыми в системе появляются сульфаты цинка, меди, магния, кобальта, натрия, свинца, кальция, а также кварц, гематит, алунит. Сульфаты меди. Сульфаты меди представлены халькантитом CuSO45H2O, боннатитом CuSO43H2O, пуатвенитом CuSO41H2O и безводным сульфатом меди CuSO4. Зависимости для сульфатов меди построены при концентрации воды 0,5 моль, это концентрация, при которой в системе присутствуют все четыре сульфата. Для халькантита и боннатита схожи кривые изменения концентрации от температуры. Причем при уменьшении 3-водного сульфата одновременно пропорционально идет увеличение 1-водного (пуатвени-та), он образуется непосредственно на боннатите как продукт его дегидратации (рис. 1). Рассмотрение моделируемой системы в диапазоне температур 0-45°С дает возможность построения рядов сульфатов по стадийности их образования в зависимости от концентрации воды от 55 до 0,51 моли. Ряды сульфатов: 0-45°С CuSO4-5H2O + CuSO4-3H2O + CuSO^^O ^ CuSO4-5H2O + CuSO4-3H2O + CuSO41H2O + CuSO4 Рис. 1. Зависимость изменения концентрации сульфатов меди от температуры Температура T указана в °С, Variant - число вариантов решения модели; число молей воды в моделируемой 0,51 моль Fig. 1. The dependence of copper sulfate concentration change on temperature Temperature T is indicated in °С, Variant is the number of model solutions; the number of moles of water in the simulated 0.51 mol Сульфаты магния. Среди сульфатов магния в модели присутствуют меридианит MgSO411H2O, эпсомит MgSO4 7H2O, пентагидрит MgSO4 5H2O, кизерит MgSO41H2O, а также магнийсодержащий минерал группы галотрихита - пиккеренгит MgAl2(SO4)422H2O, причем эпсомит появляется при температуре 10°С, а при 40°С и пентагидрит, кизерит присутствует в ассоциации с пиккеренгитом, что и наблюдалось при натурных исследованиях. Мери-дианит устойчив при низких температурах воздуха. Сульфаты магния представлены на двух графиках, так как поля устойчивости этих минералов разные, при 5 молях воды в системе отсутствует кизерит, а при 0,5 молях - меридианит и эпсомит. Одновременно при уменьшении концентрации меридианита растет концентрация пиккеренита, который, в свою очередь уменьшаясь, при определенном количестве молей воды дает начало появлению и росту эпсоми-ту, а эпсомит - пентагидриту (рис. 2, а). Пента-гидрит в свою очередь переходит в кизерит, при этом пиккеренгит в системе остается, хотя концентрация его уменьшается (рис. 2, b). MgS04(H20)11 MgS04(H20)5 v MgAI2(S04)4(H20)22 -*-MgS04(H20)7 Рис. 2. Зависимость изменения концентрации сульфатов магния от температуры Температура T указана в °С, Variant - число вариантов решения модели. а - число молей воды в моделируемой системе равно 5 моль; b - число молей воды в моделируемой системе равно 0,5 моль b a Fig. 2. The dependence of the change in the concentration of magnesium sulfate temperature Temperature T is indicated in °C, Variant is the number of model solutions; a - the number of moles of water in the simulated system is 5 mol; b - is the number of moles of water in the simulated system is 0.5 mol Рассмотрение моделируемой системы в диапазоне температур 0-45°С дает возможность построения рядов сульфатов по стадийности их образования в зависимости от концентрации воды от 55 до 0,51 моли. Большое разнообразие рядов отмечается для сульфатов магния, поля устойчивости меридианита, эпсоми-та и пентагидрита ограничены температурой. Ряды сульфатов: 0-10°С MgSO4- IIH2O ^ MgSO4- IIH2O + + MgAh(SO4V22H2O ^ MgSO4-1H2O + + MgAh(SO4V22H2O 10°С MgSO4- IIH2O ^ MgSO4- IIH2O + + MgAh(SO4V22H2O + MgSO4-7H2O ^ ^ MgSO4-1H2O + MgAh(SO4V22H2O 15°С MgSO4- 11H2O + MgSO4-7H2O ^ MgSO4- 11H2O + + MgAl2(SO4)4-22H2O + MgSO4-7H2O ^ ^ MgSO4-1H2O + MgAh(SO4V22H2O 20-35°С MgSO4-7H2O ^ MgSO4-7H2O + + MgAl2(SO4)4'22H2O ^ MgAh(SO4V22H2O + + MgSO4-1H2O 40°С MgSO4-7H2O + MgSO4-5H2O ^ MgSO4'7^O + MgAl2(SO4)4-22H2O ^ MgAl2(SO4)4'22H2O + MgSO4lH2O 45°С MgSO4-7H2O + MgSO4-5H2O ^ MgSO4'7^O + MgSO4-5H2O + MgAl2(SO4)4'22H2O ^ MgAl2(SO4)4'22H2O + MgSO4-1H2O Сульфаты цинка. Минералы цинка представлены госларитом ZnSO4 7H2O, бианкитом ZnSO4 6H2O, сандеритом ZnSO42H2O, ганнингитом ZnSO41H2O и безводным сульфатом цинка ZnSO4. График для сульфатов цинка разделен для наглядности изменения концентрации безводного сульфата (рис. 3, а, b). Вообще, присутствие всех 4 сульфатов возможно во всем диапазоне температур и практически при любой концентрации воды от 55 до 0,5 моль, исключение составляет температура 0°С, при которой ZnSO4 появляется при концентраций воды 0,5 моль, что вполне объяснимо. i-ZnS04(H20)7 *ZnS04(H20>6 V ZnS04(H20) ■*ZnSQ4|H2Q)2 Рис. 3. Зависимость изменения концентрации сульфатов цинка от температуры Температура Т указана в °С, Variant - число вариантов решения модели; а - число молей воды в моделируемой системе равно 5 моль; b - число молей воды в моделируемой системе равно 0,5 моль a Ъ Fig. 3. The dependence of zinc sulfate concentration variation on temperature Temperature T is indicated in ° C, Variant is the number of model solutions; a - the number of moles of water in the simulated system is 5 mol; b - the number of moles of water in the simulated system is 0.5 mol Ряды сульфатов: 0°C ZnS04-7H20 + ZnS04-6H20 + ZnS04-2H20 + + ZnS04- 1H20 -> ZnS04-6H20 + ZnS04-2H20 + + ZnS04- 1H20 + ZnS04 5°C ZnS04-7H20 + ZnS04-6H20 + ZnS04-2H20 + + ZnS04- 1H20 -> ZnS04-7H20 + ZnS04-6H20 + + ZnS04-2H20 + ZnS04- 1H20 + ZnS04 -> -» ZnS04-6H20 + ZnS04-2H20 + ZnS04- 1H20 + ZnS04 10-15°C ZnSO4-7H2O + ZnSO4-6H2O + ZnSO4-2H2O + + ZnSO41H2O + ZnSO4 ^ ZnSO4-6H2O + + ZnSO4-2H2O + ZnSO4-1H2O + ZnSO4 20°C ZnS04-6H20 + ZnS04-2H20 + ZnS04- 1H20 + + ZnS04 -> ZnS04- 7H20 + ZnS04-6H20 + + ZnS04-2H20 + ZnS04- 1H20 + ZnS04 -> ZnS04-6H20 + ZnS04-2H20 + ZnS04- 1H20 + ZnS04 25-45°C ZnS04-6H20 + ZnS04-2H20 + ZnS04- 1H20 + ZnS04 Несмотря на большое разнообразие сульфатов, образующихся на геохимическом барьере, в растворе сохраняется повышенная концентрация цинка, меди, магния и других металлов. Это объясняется нахождением природно-техногенной системы в кульминационной стадии процесса окисления, которая характеризуется наличием большего количества водорастворимых сульфатов металлов, сульфиды которых менее устойчивы в окислительных условиях, в первую очередь сфалерит, халькопирит и галенит [Русаль, 2015]. Согласно В.П. Зверевой [Зверева, Ястремская, Лысенко, 2014], окисление сфалерита приводит к образованию высококонцентрированных растворов, но из них не происходит кристаллизации техногенных минералов, а цинк и сера полностью переходят в раствор. В растворе при малом количестве молей воды содержание цинка, меди и магния и серы достаточно высокое в пределах n10-8-n-10-3, рН колеблется в пределах 0,48-1,64, Eh = 1,13-1,22. Заключение Физико-химическое моделирование дает возможность рассмотреть процесс образования гипергенных сульфатов, их количество и стадийность появления. Первыми в системе появляются сульфаты магния, цинка и меди. Сульфаты магния неустойчивые во всем диапазоне температур, среди них выделяются низкотемпературные - меридианит и высокотемпературные - пентагидрит. Кизерит присутствует в системе в ассоциации с пиккеренгитом. Сульфаты цинка и меди стабильны во всем диапазоне температур, лишь безводные их члены при низких температурах появляются в системе при малом количестве молей воды. Образующиеся растворы в моделируемой системе высокоминерализованные, даже при малом количестве молей воды в модели в растворе присутствует большое количество металлов и серы, которые определяют низкие значения рН и высокую реакционную активность.
Ключевые слова
Шерловогорское месторождение,
сульфаты меди,
сульфаты магния,
сульфаты цинка,
программный комплекс «Селектор»,
модель,
ряды сульфатов,
Sherlovogorsk deposit,
copper sulphate,
magnesium sulphate,
zinc sulphate,
program complex "Selektor"Авторы
Русаль Ольга Сергеевна | Институт природных ресурсов, экологии и криологии СО РАН | младший научный сотрудник, лаборатория геохимии и рудогенеза | sergutskaya@mail.ru |
Всего: 1
Ссылки
Авдонин В.Н. Техногенное окисление сульфидов Красногвардейского месторождения на Урале // Материалы по минералогии месторождений Урала. Свердловск : УНЦ АН СССР, 1984. С. 63-69
Авченко О.В., Чудненко К.В., Александров И.А. Основы физико-химического моделирования минеральных систем. М. : Наука, 2009. 232 с
Белогуб Е.В., Щербакова Е.П., Никандрова Н.К. Сульфаты Урала: распространенность, кристаллохимия, генезис. М. : Наука, 2007. С. 30-31
Бортникова С.Б. Геохимия тяжелых металлов в техногенных системах (вопросы формирования, развития и взаимодействия с компонентами экосферы) : автореф. дис. ... д-ра геол.-минерал. наук. Новосибирск, 2001. 48 с
Двуреченская С. С. Минералогия зоны гипергенеза одного из золото-серебряных месторождений Северо-Востока СССР : автореф. дис. ... канд. геол.-минерал. наук. М. : МГУ, 1987. 20 с
Жданов Ю.Я. Минералогия зоны окисления Купольного олово-серебряного месторождения // Геология и полезные ископаемые центральной части Главного металлогенического пояса Северо-Востока СССР. Якутск : Якутский научный центр СО АН СССР, 1989. С. 133-151
Зверева В.П., Ястремская Е.А., Лысенко А.И. Оценка влияния рудничных вод полиметаллических и касситерит-сульфидных месторождений Дальнегорского района на гидросферу методом физико-химического моделирования // Экологическая химия. 2014. № 23 (3). С. 159-166
Ерёмин О.В. Криоминералогенез в зоне окисления Удоканского месторождения : дис. ... канд. геол.-минерал. наук. Чита, 2004. 102 с
Карасик М.А. Поведение кобальта в зоне окисления сульфидных месторождений // Вопросы минералогии, геохимии и петрографии. М. : Изд-во АН СССР, 1946. С. 275-282
Карпов И.К. Оптимальное программирование в физико-химическом моделировании обратимых и необратимых процессов минералообразования в геохимии // Ежегодник-1970. СибГЕОХИ. Иркутск, 1971. С. 372-383
Карпов И.К. Разработка теоретических основ физико-химического моделирования природных процессов минералообразования на ЭВМ // Фундаментальные исследования. Наука о Земле. Новосибирск : Наука, 1977. С. 91-94
Карпов И.К. Физико-химическое моделирование на ЭВМвгеохимии. Новосибирск : Наука, 1981. 247 с
Карпов И.К., Чудненко К.В., Артеменко М.В., Бычинский В.А., Кулик ДА. Термодинамическое моделирование геологических систем методом выпуклого программирования в условиях неопределенности // Геология и геофизика. 1999. Т. 40, № 7. С. 971-988
Касаткин А.В., Клопотов К.И., Плашил Я. Гипергенные минералы Шерловой Горы // Минералогический альманах: Шерловая Гора. 2014. Т. 19, вып. 2. С. 94-137
Павлюкова В. А., Маркович Т.И. Геохимические процессы в криогенных зонах окисления сульфидных месторождений // Химия в интересах устойчивого развития. 2006. № 1. С. 89-96
Птицын А.Б., Абрамова В.А., Маркович Т.И., Эпова Е.С. Геохимия криогенных зон окисления. Новосибирск : Наука, 2009. 88 с
Русаль О. С. Гипергенные сульфаты группы роценита Шерловогорского олово-полиметаллического месторождения Забайкальского края // Вестник Томского государственного университета. 2015. № 400. C. 414-419
Сергутская О. С. Прогноз образования сульфатов магния в карьере Шерловогорского месторождения Забайкальского края // Известия Сибирского отделения секции наук о Земле Российской академии естественных наук. Геология, поиски и разведка рудных месторождений. 2013. № 2 (43). С. 67-71
Чудненко К.В. Термодинамическое моделирование в геохимии: теория, алгоритмы, программное обеспечение, приложения. Новосибирск : Гео, 2010. 287 с
Щербакова Е.П. Сульфатизация земной коры - последствие технической деятельности человека // Уральский минерал. Миасс : ИМин УрО РАН, 1995. № 5. С. 23-27
Эпова Е.С. Геоэкологические аспекты поведения химических элементов в условиях криогенной зоны окисления на примере Удоканского месторождения (Восточное Забайкалье) : дис. ... геол.-минерал. наук. Томск, 2014. 129 с
Юргенсон Г. А. Зона окисления в многолетнемерзлотных породах // Записки Всесоюзного минералогического общества. 1997. Ч. CXXVI, № 5. С. 15-27
Яхонтова Л.К. Поведение кобальта в зоне окисления Дашкесанского железорудного месторождения // Труды минералогического музея АН СССР. 1961. Вып. 11. С. 123-138
Jambor J.L., Boyle R.W. Gunningite, a new zinc sulphate mineral from the Keno Hill-Galena area // Yukon, Canadian Mineralogist. 1962. № 7. P. 209-218
Jambor J.L., Nordstrom D.K., Alpers C.N. Metal-sulfate Salts from Sulfide Mineral Oxidation // Sulfate Minerals. Rev. Miner. Geochem. 2000. V. 40. P. 303-350
Marion G.M. Kargel J.S., Catling D.C. Modeling ferrous-ferric iron chemistry with application to martian surface geochemistry // Geochimica et Cosmochimica Acta. 2008. V. 72. P. 242-266
Marion G.M., Mironenko M.V., Roberts M.W. FrezChem: A geochemical model for cold aqueous solutions // Computers & Ge-osciences. 2010. V. 36. P. 10-15
Moncur M.C. Hydrogeochemical, mineralogical and microbial processes occurring in old sulfide-rich tailing: Thesis doctor of philosophy in earth sciences. Waterloo, Ontario, Canada, 2015. 249 p
Nordstrom D.K. Aqueous pyrite oxidation and the consequent formation of secondary iron minerals // Acid sulfate weathering. 1982. V 3. P. 37-39