Реконструкция процессов гидротермального рудообразования методами равновесно-динамического моделированияНИР

Reconstruction of hydrothermal ore-forming processes using equilibrium-dynamic simulation methods

Источник финансирования НИР

грант РФФИ

Этапы НИР

# Сроки Название
1 1 января 2020 г.-26 декабря 2020 г. Реконструкция процессов гидротермального рудообразования методами равновесно-динамического моделирования (этап 2020 г.)
Результаты этапа: Цель проекта – развитие теории гидротермального рудообразования на основе детальных геохимических исследований, совершенствования программного обеспечения и разработки методов равновесно-динамического моделирования природных процессов на ЭВМ. Задачи исследования направлены на реконструкцию условий рудообразования при формировании гидротермальных объектов различных типов (жильных полиметаллических и эпитермальных золоторудных). Получены новые данные по распределению РЗЭ в пробах по полному разрезу через рудное тело Цагарсар (11 проб, шаг 5 см, рудная зона Цагарсар, штольня 49, месторождение Джими). Установлено, что в образовании жилы приняли участие гидротермальные растворы из разных комбинированных и эволюционирующих источников. Соотношение пород субстрата (метаморфиты PR3-PZ1 и граниты PZ3) в комбинированных источниках можно установить по характерным чертам спектров РЗЭ по интервалам опробования в разрезах. Эти данные лежат в основе новых вариантов равновесно-динамическких моделей жильного рудообразования. В моделях комбинированные источники рудных компонентов сменяются во времени. Построены и исследованы термодинамические модели, по которым можно оценить влияние смешения флюидных потоков из разных источников на процессы жильного полиметаллического рудообразования. Результаты моделирования со сменой во времени областей мобилизации рудных компонентов позволяют получить значительно лучшее согласие с природными данными. Верификация моделей проведена по геохимическим данным для месторождения Джими. С использованием программ серии Optim проведена ревизия термодинамических свойства растворенных форм золота при гидротермальных условиях по новым экспериментальным данным (Trigub et al., 2017; Zotov et al., 2018). Схема комплексообразования золота в сернокислых растворах при гидротермальных условиях проверена путем моделирования эксперимента по данным (Pokrovski et al., 2015; Kokh et al, 2016). Разработана программа HKF_full, позволяющая рассчитывать термодинамические свойства компонентов водного раствора, используемые в модели HKF (Helgeson-Kirkham-Flowers) во всём диапазоне условий, используемых этой моделью. Помимо ионов и комплексов водного раствора программа рассчитывает также термодинамические свойства растворителя (самой воды) используя модели HGK (Haar-Gallagher-Kell) или WP (Wagner-Pruss) в зависимости от выбора пользователя. Программа позволяет также пользователю выбор модели диэлектрической проницаемости воды – UF (Uematsu-Franck) или F&al (Fernandez et al). По умолчанию применяются модели, принятые в HKF.
2 23 марта 2021 г.-28 декабря 2021 г. Реконструкция процессов гидротермального рудообразования методами равновесно-динамического моделирования (этап 2021 г.)
Результаты этапа: Методами равновесно-динамического моделирования оценена возможность переотложения ранее сформированного вещества жильных руд (J2) на карбонатном барьере (J3) при повторной активизации гидротермальной системы. Установлено, что поздние низкотемпературные гидротермальные растворы могут обеспечить пиритизацию, окварцевание и отложение Pb-Zn руд в карбонатных породах. Метод термодинамического моделирования использован для определения роли геологических факторов (кипение/конденсация, охлаждение, смешение, метасоматоз) в образовании золотых месторождений типа «high sulfidation». Показано, что исследованные факторы в эпитермальных условиях действуют в сторону увеличения растворимости золота. Разработана программа AD_full, позволяющая рассчитывать термодинамические свойства растворённых в воде газов. Эти свойства описываются моделью AD (Akinfiev-Diamond) и рассчитываются программой во всём диапазоне условий, используемых этой моделью. Модель Акинфиева-Даймонда представляет растворённые в воде газы в виде нейтральных частиц, имеющих аналогичную газам стехиометрию. Помимо свойств этих частиц программа рассчитывает также и некоторые свойства растворителя (самой воды) используя модель HGK (Haar-Gallagher-Kell) или WP (Wagner-Pruss) в зависимости от выбора пользователя. Программа позволяет рассчитывать свойства газов, растворённых не только в жидкой, но и в газообразной воде.
3 25 марта 2022 г.-28 декабря 2022 г. Реконструкция процессов гидротермального рудообразования методами равновесно-динамического моделирования (этап 2022 г.)
Результаты этапа: Цель – развитие теории гидротермального рудообразования на основе детальных геохимических исследований, совершенствования программного обеспечения и разработки методов равновесно-динамического моделирования природных процессов на ЭВМ. Задачи исследования направлены на реконструкцию условий рудообразования при формировании гидротермальных объектов различных типов (жильных полиметаллических и эпитермальных золоторудных). Получены новые данные по распределению РЗЭ в пробах по разрезу через рудное тело Цагарсар (месторождение Джими). Можно предположить, что комбинированный источник рудных компонентов значительно эволюционирует во времени. В комбинированном источнике, вероятно, преобладают амфиболиты (PR3-PZ1). На начальных этапах формирования рудоносных растворов из амфиболитов предпочтительно выщелачиваются тяжелые РЗЭ. Возможно, что со временем доля кислых пород в источнике (граниты или кристаллические сланцы) увеличивается. Проведено обобщение данных по РЗЭ для месторождения Джими (80 рудных проб и 14 типов вмещающих пород). Показано, что вмещающие породы являлись комбинированными источниками рудных компонентов и РЗЭ. Соотношения масс пород субстрата в источниках могут значительно изменяться на разных временных интервалах формирования рудных тел. В 80% рудных проб суммарное содержание РЗЭ не превышает 35 г/т. Построены и исследованы термодинамические модели, по которым можно оценить влияние смешения флюидных потоков из разных источников на процессы жильного полиметаллического рудообразования. Результаты моделирования со сменой во времени областей мобилизации рудных компонентов показывают хорошее согласие с данным для месторождения Джими. Методами равновесно-динамического моделирования оценена возможность переотложения ранее сформированного вещества жильных руд (J2) на карбонатном барьере (J3) при повторной активизации гидротермальной системы. Установлено, что поздние низкотемпературные гидротермальные растворы могут обеспечить пиритизацию, окварцевание и отложение Pb-Zn руд в карбонатных породах. Проведен анализ возможных факторов отложения золота в условиях вулканогенной гидротермальной системы типа «high sulfidation». С использованием программ серии Optim проведена проверка термодинамических свойства комплексов золота для условий гидротермального процесса. Улучшено качество термодинамических данных для форм нахождения золота в кислых богатых серой растворах при высоких температурах. Показано, что геологические факторы, обычно рассматриваемые как факторы рудоотложения золота (кипение/конденсация, охлаждение, смешение, метасоматоз) в эпитермальных условиях действуют в сторону увеличения растворимости золота. Рассчитана комбинированная модель, описывающая взаимодействие откипевшего гидротермального раствора с андезитами и аргиллизитами и смешение с вадозными водами. В модели получено осаждение кристаллического золота, однако температурный диапазон осаждения и рассчитанные минеральные ассоциации не характерны для месторождений типа «high sulfidation». В некоторых случаях рассмотрение термодинамических свойств веществ, участвующих в природных процессах, является необходимым элементом их изучения. Однако пакет HCh, моделирующий такие процессы методом равновесной термодинамики, не предоставляет полной термодинамической информации об участвующих в них веществах, ограничиваясь лишь расчётом свободной энергии Гиббса, на минимизации которой основано равновесное моделирование. Поэтому в данном проекте пакет HCh был дополнен тремя программами (HKF_full, AD_full и Ind_full), позволяющими пользователям этого пакета получать полную информацию о термодинамических свойствах главных компонентов изучаемых природных процессов: компонентов водного раствора и индивидуальных веществ (твёрдых, жидких и идеальных газов) для всего диапазона поддерживаемых пакетом условий.

Прикрепленные к НИР результаты

Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".