Разработка эффективной конструкции газодинамического реактора на основе математического моделирования

Градов, Дмитрий Владимирович

Бесплатный автореферат и диссертация по наукам о земле на тему
Разработка эффективной конструкции газодинамического реактора на основе математического моделирования
ВАК РФ 25.00.13, Обогащение полезных ископаемых

Автореферат диссертации по теме "Разработка эффективной конструкции газодинамического реактора на основе математического моделирования"

На правах рукописи

ГРАДОВ Дмитрий Владимирович

РАЗРАБОТКА ЭФФЕКТИВНОЙ КОНСТРУКЦИИ ГАЗОДИНАМИЧЕСКОГО

РЕАКТОРА НА ОСНОВЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ (ПРИМЕНИТЕЛЬНО К ПЕРЕРАБОТКЕ УПОРНЫХ ЗОЛОТОСОДЕРЖАЩИХ КОНЦЕНТРАТОВ)

Специальность 25.00.13 - Обогащение полезных ископаемых

Автореферат диссертации иа соискание ученой степени кандидата технических наук

Санкт-Петербург - 2013

005062122

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Национальный минерально-сырьевой университет «Горный»

Научный руководитель:

доктор технических наук, профессор

доктор технических наук, профессор, ФГБОУ ВПО «Национальный минерально-сырьевой университет «Горный», кафедра автоматизации технологических процессов и производств, профессор

кандидат технических наук, ООО "ТОМС инжиниринг", ведущий специалист

Ведущая организация - ЗАО «НПО РИВС»

Защита состоится 3 июля 2013 г. в 16 ч 30 мин на заседании диссертационного совета Д 212.224.03 при Национальном минерально-сырьевом университете «Горный» по адресу: 199106, Санкт-Петербург, 21-я линия, д. 2, ауд. 1303.

С диссертацией можно ознакомиться в библиотеке Национального минерально-сырьевого университета «Горный»

Автореферат разослан 03 июня 2013 года.

Сизяков Виктор Михайлович

Официальные оппоненты:

Шариков Юрий Васильевич

Ленёв Леонид Александрович

УЧЕНЫЙ СЕКРЕТАРЬ

диссертационного совета

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследований. Снижение содержания золота в добываемых рудах приводит к необходимости вовлекать в переработку упорные руды и руды с низким содержанием ценного компонента. Тиосульфатное выщелачивание золота рассматривается в качестве перспективной технологии. Одним из основных процессов, определяющих эффективность технологических схем предприятий обогащения и цветной металлургии при переработке упорных золотосодержащих руд и концентратов, является стадия выщелачивания в каскаде газожидкостных реакторов. Выбор типа и конструктивных особенностей аппарата секции выщелачивания напрямую определяет общую эффективность работы предприятия.

Повышение эффективности извлечения золота из руды напрямую зависит от химии процесса, а также массопереноса. Одним из основных критериев, по которому отдается предпочтение тому или иному способу избирательного растворения является экономическая оценка процесса. Тиосульфатное выщелачивание требует больших операционных затрат, вследствие повышенного, по сравнению с цианированием, расхода реагентов, нестабильности и ограниченности повторного использования тиосульфата, а также ограниченности выбора метода экстракции. С другой стороны тиосуль-фатный процесс работает при атмосферном давлении и температуре до 30°С. После серии техногенных катастроф многие страны начали компанию по отказу от применения цианида путем экономических санкций, что повышает привлекательность альтернативных способов.

В настоящее время ведутся исследования по совершенствованию газожидкостных реакторов. Значительный вклад в решение задач связанных с проектированием реакторов перемешивания внесли: В.М. Семенов, Коровкина Г.П., Л.Н. Брагинский, B.J1. Садовский, В.М. Барабаш, В.Е. Зеленский, A.B. Осипов, В.И. Бегичев. С помощью современных компьютерных систем способных, основываясь на математических моделях, методом конечных объемов просчитывать гидродинамику движения пульпы и газа. Одним из веду-

щих компьютерных пакетов, разработанных для инженеров, является Workbench компании Ansys. С его помощью возможно произвести оптимизацию реактора с целью получения наибольшей площади контакта трех фаз газ-жидкость-твердое вещество.

Цель работы - разработка конструкции газодинамического реактора и определение операционного режима на основе компьютерного моделирования, обеспечивающего повышение показателя извлечения золота в жидкую фазу.

Идея работы - при помощи современного компьютерного пакета разработать математическую модель гидро- и газодинамики газожидкостного реактора, апробированную методом лазерной трассировки и определить наилучшие режимные настройки для повышения эффективности избирательного растворения золота в реакторе.

Основные задачи исследований:

1. Научно-технический анализ реагентов, процессов и аппаратов, применяемых при выщелачивании золота, предварительной обработке золотосодержащих концентратов, а также методов экстракции.

2. Экспериментальные исследования по исследованию кинетики избирательного растворения окисленного золотоносного концентрата в газожидкостном реакторе.

3. Научно-технический анализ математического описания гидро- и газодинамики и разработка модели трехфазной системы в газожидкостном реакторе с помощью компьютерного пакета Workbench компании Ansys.

4. Апробация и корректировка полученной модели с использованием средств лазерной трассировочной визуализации.

5. Оптимизация конструкции газожидкостного реактора посредством скорректированной модели и определение оптимальных условий работы исследуемого реактора.

Методы исследований. В представленной работе были применены методы планирования, постановки и анализа полученных результатов. При выводе основных зависимостей применены положения системного подхода к исследованию технологических

процессов, а также методы математического и физического моделирования.

В экспериментальных исследованиях по избирательному растворению золота из промышленного концентрата, полученного из флотационного передела переработки арсенопиритной золотосодержащей руды, применялись физико-химические методы анализа материала и продуктов. Планирование экспериментов проводилось по методу Бокса-Уилсона. При построении и исследовании модели было применено современное программное обеспечение от ведущей компании, разрабатывающей компьютерные пакеты в области инженерии Апзув. Методика визуализации потоков газа и жидкости, в частности лазерная трассерная визуализация, использовалась для проверки компьютерной модели на адекватность.

Научная новизна:

1. Получены кинетические зависимости и определены технологические показатели процесса аммиачно-тиосульфатного выщелачивания золота в присутствии меди(Н) из окисленного промышленного арсенопиритного концентрата и получена зависимость извлечения золота в раствор от расхода реагентов.

2. Разработана система визуализационной апробации математической модели, описывающей поведение газожидкостных потоков в реакторе перемешивания применительно к тиосульфатному выщелачиванию.

3. Построена независимая от геометрии аппарата компьютерная модель газодинамического реактора.

Основные защищаемые положения:

1. При проектировании конструкции газожидкостных реакторов для избирательного растворения необходимо применять компьютерное моделирование с использованием программных пакетов, базирующихся на методе конечных объемов, результаты которого позволят произвести интенсификацию процесса растворения за счет увеличения площади поверхности контакта фаз, участвующих в реакции.

ремешивания при оптимальных режимных настройках, полученных с помощью компьютерной модели.

Практическая значимость работы:

1. Разработанная методика оптимизации конструкции и операционного режима газожидкостного с помощью компьютерного моделирования, основанного на методе конечных объемов, может найти свое применение при проектировании газодинамических реакторов.

2. Предложенная конструкция аппарата и режим его работы востребованы при переработке упорных золотосодержащих концентратов.

Достоверность научных результатов. Приводимые результаты, выводы и рекомендации подкреплены значительным объемом данных полученных экспериментальным путем и их соответствием теории проектирования аппаратов перемешивания для металлургической и обогатительной промышленности, применением современных методов планирования экспериментов и физико-химического анализа технологических продуктов и материалов, использованием адекватных математических моделей и статистических методов обработки данных с применением компьютерных технологий.

Апробация работы.

Содержание и основные положения диссертации докладывались и обсуждались на XXVII международном конгрессе по обогащению полезных ископаемых (Германия, г. Фрайберг) в 2012 г., на научно-технических советах по работе с аспирантами СПГГИ (ТУ) им. Г. В. Плеханова 2010-2013 г.г., на европейском конгрессе химической инженерии 2013 г.

Личный вклад автора.

Автор самостоятельно выполнил:

1. Постановку задач и разработку общей методики исследований.

2. Анализ современных методов и оборудования, применяемого для перемешивания в жидких средах, а также обзор известных методов и реагентов, используемых при выщелачивании золота.

3. Весь комплекс экспериментальных исследований в лабораторных условиях по избирательному растворению золотоносного концентрата и лазерной визуализации потоков в реакторе.

4. Разработку математической модели поведения потоков в газожидкостном реакторе с помощью компьютерного пакета Workbench и алгоритма апробации полученной модели.

Публикации. По теме диссертации опубликованы 3 научные работы в изданиях, рекомендованных ВАК Минобрнауки России.

Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения и библиографического списка. Содержит 139 страниц машинописного текста, 66 рисунков, 22 таблицы, список литературы из 85 наименований.

Автор выражает благодарность профессору Турунену И. и доценту Арто Л. за помощь и поддержку, оказанную при работе над диссертацией.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

В первой главе представлен описание и аналитический анализ методов и оборудования по избирательному растворению, применяемые при выщелачивании золота из руд и промышленных концентратов. Рассмотрены недостатки и преимущества каждого из них. Описаны механизмы и особенности тиосульфатного выщелачивания.

Приведены характеристика объекта и постановка задачи исследования.

Во второй главе описаны экспериментальные исследования в лабораторном масштабе по тиосульфатному выщелачиванию золота из окисленного промышленного арсенопиритного концентрата в реакторе перемешивания.

Третья глава посвящена описанию основ математического моделирования потоков флюидов. Представлена компьютерная модель на основе пакета Workbench компании Ansys, описывающая движение газожидкостных потоков в реакторе, с целью определения оптимальных режимных настроек аппарата перемешивания.

В четвертой главе описывается основы метода лазерной трассерной визуализации потоков. Показан алгоритм проверки с последующей корректировкой разработанной математической модели реактора, используемого при выщелачивании, на адекватность на основе метода лазерной визуализации газожидкостных потоков.

В пятой главе приводятся результаты исследования полученной модели по определению режимных настроек обеспечивающих оптимальное перемешивание на основе определения максимальной площади соприкосновении трех фаз. Описаны экспериментальные лабораторные исследования с учетом рекомендованных режимных настроек по выщелачиванию арсенопиритного концентрата, предварительно окисленного едким натром. Приведено сравнение эффективности избирательного растворения золота с учетом полученных режимных настроек.

Заключение отражает обобщенные выводы по результатам исследования в соответствии с поставленной целью и решёнными задачами.

ОСНОВНЫЕ ЗАЩИЩАЕМЫЕ ПОЛОЖЕНИЯ 1. При проектировании конструкции газожидкостных реакторов для избирательного растворения необходимо применять компьютерное моделирование с использованием программных пакетов, базирующихся на методе конечных объемов, результаты которого позволят произвести интенсификацию процесса растворения за счет увеличения площади поверхности контакта фаз, участвующих в реакции.

Компьютерное моделирование осуществляется с использованием программного пакета на языке си Fluent компании Ansys, широко используемого инженерами проектировщиками по всему миру, а также ведущими университетами в учебных и исследовательских целях. Данный программный продукт позволяет с высокой точностью прогнозировать поведение газовых и жидкостных потоков в пределах заданной геометрии и их механическое, тепловое, а также химическое взаимодействие с возможностью корректировки основных законов и параметров прогнозируемого взаимодействия.

Программа Fluent способна оперировать только потоками флюидов, т.е. жидкостями, газами и сыпучими средами, но не жест-

кими каркасами и движущимися частями механизмов. Расчет движения частиц основывается на уравнениях Навье - Стокса, описывающих гидродинамику ньютоновских жидкостей в элементарной ячейке (Рисунок 1) сетки разбивки.

""XL

Рисунок 1 - Принцип расчета элементарной ячейки

-1, 3 v I J /

ct ox.

OXl OX ,

On, дії, 2 dut , — + — ——6„

ox..

OX;

3 oXl.

+ PSi + F, (1)

где x, у, z - координаты точки, м; ui jk - мгновенная скорость в направлении х, у, z соответственно, м/с; р - плотность жидкости; р - вязкость жидкости, <5,у - дельта Кронекера; g, - ; Fj - источник.

Пульпа с содержанием твердого 30 % возможно представить как Бингамовский пластик. Таким образом, определив вязкость пульпы лабораторными методами можно, создав квазижидкость с заданной вязкостью, использовать вышеприведенный подход к описанию гидродинамики пульпы при высоких значениях напряжения сдвига в рассматриваемом реакторе перемешивания.

Описание Эйлера поля течения потока в динамики жидкости рассматривает движение жидкости как пространственные домены, через которые жидкость проходит с течением времени. Такая модель позволяет моделировать множество отдельных и к тому же взаимодействующих друг с другом потоков. Которые в свою очередь могут быть жидкостями, газами или твердыми частицами практически в любой комбинации.

Модель Эйлера в программе Fluent не различает многофазные потоки жидкостей и пульп, моделируя их как жидкостные потоки, при этом есть возможность задавать размер и форму частиц вто-

ричной по отношению к жидкости фазе. При расчете динамики потока единое давление потока делиться между всеми фазами, а уравнения скорости и неразрывности решаются для каждой фазы отдельно.

Весь расчет динамики жидкостных потоков с помощью метода конечных объемов строится на уравнениях сохранения массы:

j¡ ip-qPq) + V ■ (aqpqvq) = - ™qp 3 + ^ (2)

где Vq - скорость фазы q, mpq - массообмен из фазы p в фазу q, rhqp - массообмен из фазы q в фазу р, aq - объемная доля фазы q, Pq - плотность фазы q, Sq - источник фазы q, который по умолчанию равен 0, но может быть задан в виде импульса: Jt (aqPqvq) + V ■ (aqpqvqvq) = -aq4p + 4-Tq + aqpj + (3)

ÍRP<1 + ™pqVpq ~ ™<!pVqp) + (Fq + Flift« +

где rq - тензор напряжения фазы q, Rpq - сила взаимодействия между фазами, р - давление, Fq - внешняя сила фазы, FUftiq -сила Архимеда, iwwj - сила воображаемой массы, vpq - межфазная скорость.

Для расчета турбулентности в пакете имеется три способа: с помощью вычисления кинетической энергии турбулентности и энергии рассеивания fe-e, кинетической энергии турбулентности и ее частоты к- w и усредненное уравнение Навьера - Стокса.

Чтобы имитировать процесс перемешивания, построенную геометрическую модель необходимо разбить на две основные зоны (Рисунок 2): зону возмущения (она же зона вращения, непосредственно примыкающая к приводу с импеллером) и зону отклика (все оставшееся пространство в пределах геометрии), на которую передается возмущение. Таким образом, жидкость в первой зоне движется вокруг неподвижных привода и импеллера, передавая остаточные возмущения в зону отклика, и тем самым создается эффект вязкой среды, сопротивляющейся возникшему в ней движению. В описанном способе симуляции процесса перемешивания в реакторе есть недостаток, прежде всего связанный с передачей движения между указанными зонами и зависит от качества разбивки на контрольные

элементы перед расчетом, что приводит к увеличению либо необходимой мощности компьютера, либо времени расчета. Еще одним недостатком является объем зоны возмущения, чем она больше, тем большее искажение имеет место.

Затем геометрия была помещена в среду Fluent для задания условий симуляции: скорость вращения импеллера 200 об/мин по часовой стрелке, содержание твердой фазы (диаметром частиц 74100 мкм и плотностью 2970 кг/м3) 30 % от объема жидкой фазы и расход воздуха - 1 л/мин. После обработки данных программным пакетом были получены результаты в виде контуров и векторов полей скоростей (Рисунок 2). На рисунке 26 при увеличении видно, как смесь увлекаемая лопастями импеллера разбивается об отражательные перегородки и, снижая скорость, огибают препятствие. Последствия недостатков описанного способа можно заметить на рисунке 2а на месте стыка зон возмущения и отклика ввиду расхождений в их разбивки сеткой на ячейки. Проблема устраняется путем более детальной разбивки зон и неизбежно ведет к увеличению времени расчета.

С помощью контуров распределения фаз (Рисунок 3) можно отследить поведение подаваемого через клапан воздуха, уменьшая диапазон содержания газа. Сначала подаваемый газ плотным потоком движется к основанию импеллера. Затем, увлекаемый создаваемым импеллером потоком, газ рассеивается, занимая больший объем с меньшим содержанием, постепенно движется к поверхности и достигает содержания 0,04 %

Использование описанного выше программного продукта позволяет определить следующие оптимальные режимные настройки проектируемого газожидкостного реактора:

- насосная мощность импеллера:

- задержка газа в реакторе: ¿V = KNpn?ds (жидкость) (4)

N = Na + Nk + iV- = + + GraspgH (газ-жидкость) (5)

где р - плотность жидкости, п - частота вращения импеллера, d -диаметр импеллера, g - ускорение свободного падения, w - скорость подачи воздуха, А'Лг - критерий Эйлера, 6газ - расход газа, Ng - мощ-

действия метода трассерной визуализации и последующей обработки с помощью программного обеспечения компании Lavision Апробация компьютерной модели на адекватность осуществляется с помощью метода цифровой трассерной визуализации (Particle Image Velocimetry, PIV). Принцип действия (Рисунок 4) основан на измерении мгновенного поля скорости потока жидкости или газа в заданном сечении. Таким образом, лазер, создающий плоскость в исследуемом потоке, позволяет камере отследить специально подобранные трассирующие частицы, идентичные частицам потока и после обработки изображений софтом построить поля скоростей частиц двигающихся в потоке. Отслеживание движения пузырьков воздуха осуществляется с помощью специального флюоресцирующего вещества.

Процедура апробации требует определить критерии, по которым будет определяться валидность модели. Для определения ва-лидности модели движения флюидных потоков предлагается ис-

ность для перемешивания газожидкостной среды, - мощно-

сти, подводимые кинетической и потенциальной энергиями при подаче газа.

- мертвые зоны реактора, форма и расположение отражателей. _____________________________________

Рисунок 2 - Контуры (а) н векторы (б) полей скоростей

Вс«у Со* ^^ у МадгИоЛг (гаСгг) (пОД <Т«п*=1 ЮОСа*00)

а)

Сотот Ш Уексе,' м иуаеийс (.иг) ¡ггУ51 (Ттк-

Рисунок 3 Контуры газовой фазы с уменьшением диапазона содержания (|3)

Рисунок 5 ^Результаты лазерной визуализации, а - контур скорости движения гетерогенной смеси (жидкость - твердое), б - контур скорости движения газовой фазы, в - векторное поле скоростей движения гетерогенной смеси (жидкость - твердое), г - векторное поле скоростей движения газовой фазы

пользовать сравнение шаблонов векторов полей скоростей. Для более точного сравнения необходимо выделить матрицу вектор из общего поля и построить график изменения скорости потока по мере движения от импеллера к стенке реактора. Выборку следует производить по вертикальной оси с тем расчетом, чтобы представить каждую завихрений.

Алгоритм апробации: Шаг (1): получение набора полей векторов скоростей газа и жидкости с последующим соединением полученных полей в единое целое. Т.к. рабочее окно установки лазерной визуализации ограничено необходимо выполнить следующую последовательность действий:

- загрузка полей векторов скоростей;

- объединение фрагментов с захлестыванием через соединение радиальной и осевой составляющих векторов скоростей;

- создание симметричного фрагмента половинки реактора относительно оси привода импеллера и соединения фрагментов в единый контур (Рисунок 5);

Шаг (2): определение точек отбора по вертикальной оси с целью отбора репрезентативных матриц векторов (Рисунок 6).

Рисунок 6 - Схема определения точек отбора репрезентативных матриц-векторов Шаг (3): отбор матриц векторов и построение графиков зависимостей осевых и радиальных составляющих векторов скоростей от удаления от оси привода мешалки.

Шаг (4): сравнение общего шаблона движения газожидкостных потоков в реакторе и выделенных графиков с результатами моделирования.

у, мм

Сравнение результаты лазерной трассерной визуализации и моделирования представлены ниже:

■

«feewas^sterse

. ; ;; Щ

! Ж ; jTt^^^M Ш&, \ Щ *

• л

< •

V * Л

\ !.

^ ,, * •< * ' 4 л

' К i.•»> г > - V

t :' £

jgj > * Ч 3

| * 4 л у v S

• Ы-** \ Щ. ш

• 1 1

•1 >

г i 1 i i f / < i i i i i 4 г i f t t i i •. i i i l i

V 1 I i t I M i 1 I ^iUIUHMMII •MUlUISIillfl

•чшммммп ■ ■ • н И I 1 1 » M"

', ' ' U I j ( i H 1 H j I i

j •:; i i U^ ': -Щ11

v -'< i i j s I \ i

ic ; ii......... ц

■ 1 i

, ! I

- I

* iif'Tf':

I®

Рисунок 7 - Поля векторов скоростей жидкой фазы, полученные при помощи: а - программы Fluent компании Ansys, б - техники лазерной трассерной визуализации

PIV

Рисунок 8 — Радиальная скорость жидкой фазы на уровне 110 мм по вертикальной оси

считая от основания реактора

М/С

О 50 1С» ¡50

импеллер стенка

Рисунок 9 - Радиальная скорость жидкой фазы на уровне 145 мм по вертикальной оси

считая от основания реактора

м/с

Рисунок 10 - Радиальная скорость жидкой фазы на уровне 200 мм по вертикальной оси

считая от основания реактора Современный программный пакет Fluent от компании Ansys обладает достаточным математическим аппаратом для описания, прогнозирования и вывода необходимых параметров системы перемешивания в газожидкостном реакторе. Но, как и любая модель требует апробации на практике через разработанные для каждого отдельного случая критерии. Средства визуализации могут служить для получения контрольных параметров при проведении апробации компьютерной модели перемешивания. При сравнении результатов моделирования и PIV удобно использовать компоненты векторов полей скоростей на плоскости сечения реактора лазером. Критериями проверки модели на адекватность могут служить осевая и радиальная составляющие компоненты векторов полей скоростей, отобранных на контрольных участках изучаемого сечения реактора. Результаты апробации показали высокую сходимость, что говорит о возможности применения модели для оптимизации условий перемешивания в газожидкостном реакторе для повышения уровня мас-сообмена.

2. С целью повышения эффективности извлечения золота и снижения экологической нагрузки следует применять ам-миачно-тиосульфатное выщелачивание в реакторе перемешивания при оптимальных режимных настройках, полученных с помощью компьютерной модели.

Тиосульфат - это химическое вещество, широко применяемое в фотографии и фармацевтической промышленности. И в отличие от цианида не обладает отравляющими все живое свойствами. Главный компонент процесса является по своей сути обычным удобрением, что открывает перспективу использования хвостов в сельском хозяйстве. В реакции образуются два комплекса Аи(8203)~ и Аи(8203)23\ последний из которых более стабилен. Щелочная среда препятствует разложению тиосульфата, что имеет дополнительное преимущество в уменьшении растворимости примесей и в особенности железа. Однако механизм выщелачивания усложняется гомогенными реакциями между Си(Н) и тиосульфатом, что в свою очередь приводит к снижению концентрацию обоих компонентов. К тому же происходит разложение некоторого количества тиосульфата до тетратионата по следующей реакции:

2Си(Ш3)2+ + вЭгОз2- 2Си(8203)35" + + 8Ш3 (6) Тиосульфатное выщелачивание является трехфазным процессом газ, жидкость и твердое (окислитель, растворитель и материал).

Исследования аммиачно-тиосульфатного выщелачивания ар-сенопиритного концентрата в присутствии двухвалентной меди проводились по схеме на рисунке 11. Полученные результаты позволили описать кинетику реакции (Рисунок 12). На основе метода планирования Бокса-Уилсона на кубе позволили определить оптимальные параметры реакции.

После процедуры апробации математическая модель, описывающая перемешивание трехфазной системы в реакторе была применена к реактору, в котором проводились испытания по выщелачиванию. Время задержки воздуха в реакторе и площадь соприкосновения фаз являлись контрольными параметрами при определении скорректированных аппаратных настроек. Стартовыми условиями были условия проведения опытов по выщелачиванию промышлен-

ного арсенопиритного окисленного концентрата. В результате серии симуляций, были определены положение (клиренс) и скорость вращения импеллера.

Кек

Анализ

Флэьтрат

Атомно-адсорбцио нный анализ

АиДл? ДсиД

Приготовление проб

во

.-А

взОз-'

_ ... Ионная хромато-

Реакгор Фильтрация графш

Рисунок 11 - Схема проведения опытов по выщелачиванию

60 80 100 ПО Врсмнотйора пробы,мин

180

Рисунок 12- Кинетика извлечения золота при выщелачивании тиосульфатно-аммиачным раствором арсенопиритного концентрата После проведения серии симуляций на модели детальной разбивки были определены скорректированные условия перемешивания (7 см от поверхности дна реактора и 735 об./мин.), обеспечивающие более равномерное распределение твердой фазы внутри реактора при максимально возможном времени задержки газа.

Выходные данные симуляций были опробованы на поверочных экспериментах в автоклавном выщелачивании и показали рост эффективности извлечения золота в жидкую фазу на 8 % при доверительном интервале 2-3 %. Сравнение результатов исследований по выщелачиванию окисленного арсенопиритного концентрата с учетом корректировки по средствам компьютерного моделирования и без него представлены на рисунке 13. "¿00

s 80

2 60

1 І j ;

і ! ........" _________ц ..........1 *——-' 1—— —...... І

! Í т^':............. г ir--"

¡.......... ....... j І X 1 Í XV щ

........л/ / 0 ♦ 1 ■»■ 2

і // > 1 // ! _____J

Г...... L.............і...................... 1 1

40 60 80 ¡00 120

Время отбора пробы, мин

140

160

180

Рисунок 13 - Эффективность извлечения золота при выщелачивании тиосульфатно-аммиачным раствором арсенопиритного концентрата до и после применения компьютерного моделирования. 1 - оптимальные режимные настройки аппарата, найденные с помощью компьютерной модели, 2 - параметры найденные с помощью

метода Бокса-Уилсона

ЗАКЛЮЧЕНИЕ

Основные результаты выполненных исследований:

1. Получена кинетика и определены технологические показатели избирательного растворения золота из окисленного арсенопиритного промышленного концентрата.

2. Построена математическая модель гидродинамики гетерогенной смеси в реакторе перемешивания на базе компьютерного па кета Workbench компании Ansys и произведена апробация получен

ной модели с помощью метода лазерной визуализации путем сравнения полей векторов скоростей.

3. Разработан алгоритм процесса апробации средствами лазерной визуализации и установлены критерии оценки валидности математической модели перемешивания газожидкостного реактора.

4. Определены оптимальные условия процесса аммиачно-тиосульфатного выщелачивания золота в газодинамическом реакторе в присутствии двухвалентной меди из окисленного арсенопирит-ного промышленного концентрата.

5. По результатам лабораторных работ установлено повышение показателя извлечения золота в тиосульфатный раствор после оптимизации режимных настроек процесса и положения импеллера с помощью разработанной компьютерной модели.

ОСНОВНЫЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИИ

ОПУБЛИКОВАНЫ В СЛЕДУЮЩИХ РАБОТАХ:

1. Градов Д.В. Компьютерное моделирование физических процессов в газожидкостном реакторе/ Д.В. Градов, В.М. Сизяков, И. Турунен, А. Лаари // Записки Горного института, 2013.Т.202.С.284-287

2. Градов Д.В. Совершенствование процесса тиосульфатного выщелачивания золота по средствам оптимизации массообменных процессов в газожидкостных реакторах/ Д.В. Градов, В.М. Сизяков, И. Турунен, А. Лаари // Записки Горного института, 2013.Т.202.С.69-71

3. Градов Д.В. Исследование влияния процесса окисления концентрата на тиосульфатное выщелачивание золота/ Д.В. Градов, В.М. Сизяков, И. Турунен, А. Лаари // Записки Горного института, 2013 .Т.202.С.72-74

РЩ Горного университета. 31.05.2013. 3.329. Т.100 экз. 199106 Санкт-Петербург, 21-я линия, д.2

Текст научной работыДиссертация по наукам о земле, кандидата технических наук, Градов, Дмитрий Владимирович, Санкт-Петербург

Федеральное государственное бюджетное образовательное учреждение высшего

«НАЦИОНАЛЬНЫЙ МИНЕРАЛЬНО-СЫРЬЕВОЙ УНИВЕРСИТЕТ «ГОРНЫЙ»

РАЗРАБОТКА ЭФФЕКТИВНОЙ КОНСТРУКЦИИ ГАЗОДИНАМИЧЕСКОГО РЕАКТОРА НА ОСНОВЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ (ПРИМЕНИТЕЛЬНО К ПЕРЕРАБОТКЕ УПОРНЫХ ЗОЛОТОСОДЕРЖАЩИХ КОНЦЕНТРАТОВ)

Специальность 25.00.13 - Обогащение полезных ископаемых Диссертация на соискание ученой степени кандидата технических наук

профессионального образования

На правах рукописи

П4?ГН ТАГЮ П7

Градов Дмитрий Владимирович

Научный руководитель: доктор технических наук, профессор В.М. Сизяков

Санкт-Петербург - 2013

Оглавление

Введение.............................................................................................................................4

ГЛАВА 1 АНАЛИТИЧЕСКИЙ ОБЗОР, ОПИСАНИЕ И ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЙ.............................................................................................................8

1.1 Основы процесса выщелачивания...............................................................................8

1.1.1 Основные агенты выщелачивания..........................................................................................9

1.1.2 Обзор методов и оборудования по избирательному растворению....................................10

1.1.3 Методы предварительной обработки руд............................ ................................................14

1.1.4 Способы извлечения и концентрирования ценного компонента.......................................18

1.2 Механизм цианирования............................................................................................22

1.3 Механизм тиосульфата ого выщелачивания.............................................................29

1.3.1 Химический механизм аммиачно-тиосульфатного выщелачивания................................29

1.3.2 Роль веществ в химии процесса............................................................................................32

1.3.3 Оптимальные условия............................................................................................................34

1.3.4 Недостатки аммиачно-тиосульфатного выщелачивания...................................................37

1.3.5 Пути совершенствования процесса......................................................................................37

1.3.6 Расход реагентов.....................................................................................................................41

1.3.7 Нейтрализация тиосульфата..................................................................................................43

1.4 Характеристика объекта исследований.....................................................................44

1.5 Постановка задачи исследований..............................................................................46

ГЛАВА 2 ИССЕДОВАНИЯ ПРОМЫШЛЕННОГО ФЛОТАЦИОННОГО КОНЦЕНТРАТА АРСЕНОПИРИТНОЙ РУДЫ НА ИЗБИРАТЕЛЬНОЕ РАССТВОРЕНИЕ АММИАЧНО-ТИОСУЛЬФАТНЫМ РАССТВОРОМ В ПРИСУТСТВИИ МЕДИ..................................................................................................47

2.1 Планирование эксперимента......................................................................................47

2.2 Проведение эксперимента и результаты...................................................................62

Выводы по главе...............................................................................................................67

ГЛАВА 3 МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ДВИЖЕНИЯ ПОТОКОВ В ГАЗОЖИДКОСТНОМ РЕАКТОРЕ................................................................................68

3.1 Теоретическое описание основ применяемого компьютерного моделирования на основе метода конечных элементов................................................................................69

3.2 Разработка модели гидрогазодинамики заданной геометрии..................................84

3.2.1 Создание геометрической модели исследуемого реактора................................................84

3.2.2 Разбивка на контрольные объемы........................................................................................87

3.2.3 Задание параметров модели перемешивания и граничных условий.................................90

3.2.4 Результаты симуляции реактора перемешивания...............................................................93

Выводы по главе...............................................................................................................95

ГЛАВА 4 АПРОБАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ С ПОМОЩЬЮ ЛАЗЕРНОЙ ТРАССЕРНОЙ ВИЗУАЛИЗАЦИИ............................................................96

4.1 Основы техники лазерной визуализации гидрогазодинамических потоков..........97

4.2 Описание установки лазерной трассировочной визуализации и результаты измерений.......................................................................................................................109

4.3 Апробация компьютерной модели реактора...........................................................115

Выводы по главе.............................................................................................................122

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

МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ............................................................123

Выводы по главе.............................................................................................................129

Заключение.....................................................................................................................130

Список литературы.........................................................................................................132

Введение

Снижение содержания золота в добываемых рудах приводит к необходимости вовлекать в переработку упорные руды и руды с низким содержанием ценного компонента. Тиосульфатное выщелачивание золота рассматривается в качестве перспективной технологии. Одним из основных процессов, определяющих эффективность технологических схем предприятий обогащения и цветной металлургии при переработке упорных золотосодержащих руд и концентратов, является стадия выщелачивания в каскаде газожидкостных реакторов. Выбор типа и конструктивных особенностей аппарата секции выщелачивания напрямую определяет общую эффективность работы предприятия.

Повышение эффективности извлечения золота из руды напрямую зависит от химии процесса, а также массопереноса. Одним из основных критериев, по которому отдается предпочтение тому или иному способу избирательного растворения является экономическая оценка процесса. Тиосульфатное выщелачивание требует больших операционных затрат, вследствие повышенного, по сравнению с цианированием, расхода реагентов, нестабильности и ограниченности повторного использования тиосульфата, а также ограниченности выбора метода экстракции. С другой стороны тиосульфатный процесс работает при атмосферном давлении и в температурном диапазоне 25 30°С. После серии техногенных катастроф многие страны начали компанию по отказу от применения цианида путем экономических санкций, что повышает привлекательность альтернативных способов.

В настоящее время ведутся исследования по совершенствованию газожидкостных реакторов. Значительный вклад в решение задач связанных с проектированием реакторов перемешивания внесли: В.М. Семенов, Коровкина Г.П., JI.H. Брагинский, B.JI. Садовский, В.М. Барабаш, В.Е. Зеленский, A.B. Осипов, В.И. Бегичев. С помощью современных компьютерных систем способных,

основываясь на математических моделях, методом конечных объемов просчитывать гидродинамику движения пульпы и газа. Одним из ведущих компьютерных пакетов, разработанных для инженеров, является Workbench компании Ansys. С его помощью возможно произвести оптимизацию реактора с целью получения наибольшей площади контакта трех фаз газ-жидкость-твердое.

Научная новизна результатов диссертационного исследования заключается в следующем:

- Получены кинетические зависимости и определены технологические показатели процесса аммиачно-тиосульфатного выщелачивания золота в присутствии меди(П) из окисленного промышленного арсенопиритного концентрата и получена зависимость извлечения золота в раствор от расхода реагентов.

- Разработана система визуализационной апробации математической модели, описывающей поведение газожидкостных потоков в реакторе перемешивания применительно к тиосульфатному выщелачиванию.

- Построена независимая от геометрии аппарата компьютерная модель газодинамического реактора.

Практическая значимость работы:

- Разработанная методика оптимизации конструкции и операционного режима газожидкостного с помощью компьютерного моделирования, основанного на методе конечных объемов, может найти свое применение при проектировании газодинамических реакторов.

- Предложенная конструкция аппарата и режим его работы востребованы при переработке упорных золотосодержащих концентратов.

Достоверность научных результатов. Приводимые результаты, выводы и рекомендации подкреплены значительным объемом данных полученных

экспериментальным путем и их соответствием теории проектирования аппаратов перемешивания для металлургической и обогатительной промышленности, применением современных методов планирования экспериментов и физико-химического анализа технологических продуктов и материалов, использованием адекватных математических моделей и статистических методов обработки данных с применением компьютерных технологий.

Автор самостоятельно выполнил:

1. Постановку задач и разработку общей методики исследований.

4. Разработку математической модели и алгоритма ее апробации трехфазной системы в газожидкостном реакторе с помощью компьютерного пакета Workbench.

Апробация работы

Содержание и основные положения диссертации докладывались и обсуждались на XXVII международном конгрессе по обогащению полезных ископаемых (Германия, г. Фрайберг) в 2012 г., на научно-технических советах по работе с аспирантами СПГГИ(ТУ) им. Г. В. Плеханова 2010-2013 г.г., на Европейском Конгрессе Химической Инженерии 2013 г.

Основные защищаемые положения:

1. При проектировании конструкции газожидкостных реакторов для избирательного растворения необходимо применять компьютерное моделирование при помощи современных программных пакетов, базирующихся на методе конечных

объемов, что обеспечивает интенсификацию процесса растворения за счет увеличения площади поверхности контакта фаз, участвующих в реакции.

2. С целью повышения эффективности извлечения золота и снижения экологической нагрузки следует применять аммиачно-тиосульфатное выщелачивание в реакторе установленного типа перемешивания при оптимальных режимных настройках, полученных с помощью компьютерной модели.

ГЛАВА 1 АНАЛИТИЧЕСКИЙ ОБЗОР, ОПИСАНИЕ И ПОСТАНОВКА

ЗАДАЧИ ИССЛЕДОВАНИЙ

1.1 Основы процесса выщелачивания

Выщелачивание представляет собой процесс избирательного растворения компонентов и перевод их из твердой фазы в жидкую посредствам растворителя. Процесс называется прямым, если осуществляется растворение полезного компонента из руд, концентратов, металлургических продуктов или обратным, если происходит растворение пустой породы или ненужных в дальнейшем переделе компонентов с целью повышения концентрации полезного компонента.

Руда

Концентрат Металлы

Рисунок 1.1- Общая схема гидрометаллургического процесса[46]

1.1.1 Основные агенты выщелачивания

Растворитель обычно выбирается по следующим критериям:

- Растворимость. Полнота и скорость растворения материала

- Избирательность. Агент должен растворять только желаемый компонент

- Экономичность. Потери реагента неизбежны, а значит, экономическая целесообразность должна соблюдаться

- Способность к восстановлению позволяет использовать реагент несколько раз и тратить меньше усилий на его нейтрализацию в последующих стадиях

- Корозийность. Его способность разрушать металлы снижает время эксплуатации оборудования, а также повышает его стоимость посредством повышения стандартов металлов при заказе аппаратов.

Наиболее идеальным растворителем является водам из-за её экономических плюсов: нейтральная рН и безопасность в эксплуатации. Но, к сожалению, её активность по отношению к минералам ограничена. Поэтому применяются другие растворители, которые можно разделить на следующие группы:

Кислоты - наиболее используемые реагенты при выщелачивании и в свою очередь могут выступать в качестве:

- Окислителей

- Восстановителей

- Нейтральных веществ, т.е., только атом водорода участвует в реакции

Ярким примером является "Царская водка", состоящая на одну четвертую из

азотной и на три четверти из соляной кислот, способна растворять золото.

Щелочи, такие как КаОН или 1ЧН4ОН, стоят дороже серной кислоты, но благодаря их незначительной корозийности применяются для выщелачивания руд с большим содержанием пустой породы, которая в свою очередь не реагирует со щелочами.

Водные растворы солей - наиболее часто применяемые при получении золота на обогатительных и металлургических предприятиях. Самым распространенным в промышленности и изученным является цианистый раствор натрия или калия, последний используется реже. Его приготовление происходит при адсорбции цианистого газа водным раствором гидроксида натрия:

ИСК + ЫаОН №СЫ + Н20 (1.1)

Окислители и восстановители. При избирательном растворении золота иногда необходимо растворить некоторые минералы, для чего требуется введение дополнительных окислителей, например кислород, озон, пероксид водорода, трех валентное железо, диоксид серы и другие.[19]

1.1.2 Обзор методов и оборудования по избирательному растворению

На практике выбор метода зависит от производительности и типа перерабатываемой руды, а также концентрации извлекаемого компонента. При выщелачивании золота обычно применяют следующие методы:

Кучное выщелачивание. Для осуществления данного метода необходимо подготовить и выровнять площадь земли, которая затем покрывается специальным слоем, обладающим водонепроницаемыми и антикоррозийными свойствами. На подготовленный участок складируется руда, и монтируются трубопроводные системы для подачи и орошения выщелачивающего раствора и окислителей, например кислорода. После определенного периода выщелачивания полученный насыщенный раствор направляется на последующие стадии разделения и концентрации. Восстановленный раствор снова поступает на стадию кучного выщелачивания. Типичная схема кучного выщелачивания представлена на рисунке 1.2.

Рисунок 1.2 - Общая схема кучного выщелачивания[79]

Просачивание или чановое выщелачивание чаще всего применяется для пористых руд, в которых ценный компонент не связан с другими компонентами в минеральных агрегатах и легкодоступен для раствора. Конструктивно чан для описываемого способа представляет собой кузов с перфорированным или сеточным покрытием дна для отделения насыщенного раствора от руды (Рисунок 1.3). Руду загружают в емкость, при этом более мелкую фракцию кладут ближе ко дну, а более крупную сверху. Подаваемый сверху раствор, образуя каналы в крупной фракции, попадет на так называемый фильтр из более мелкой фракции. Насыщенный раствор, также как и в кучном выщелачивании, собирается и подается на дальнейшие этапы разделения и восстановления раствора для повторной его подачи в голову процесса

выщелачивания. Процесс является периодическим и после выщелачивания руда выгружается из чана, после чего операции повторяются. Процесс длится от двух до четырех дней. Из преимуществ можно назвать минимальное потребление рабочего раствора, высокое содержание ценного компонента в насыщенном растворе, а также отсутствие относительно дорогих сгустителей или фильтров.

Рисунок 1.3 - Чановое выщелачивание или просачивание[80]

Выщелачивание в реакторах с перемешиванием используется для руд со

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

мм

WnM

□

ft*

Overflow Inlet

Oxygen

Rubber lining

Acid

Manhole

Outlet

Drain

тм

Oullel

Oxygon

Drau

□

Steam

Г I t

Air arid H2SC>4

Level indicator —

Pulp^nput

...........i

Pulp

discharge

Рисунок 1.4 - Типичное оборудование для выщелачивания, а - реактор с механическим перемешиванием, б - реактор с пневматическим перемешиванием[18], в - вертикальный автоклав с механическим перемешиванием, г - вертикальный автоклав с пневматическим перемешиванием[46]

Непрерывное перемешивание обеспечивает доступ свежего еще не вступившего в реакцию раствора к обнажаемой поверхности растворяемого минерала и отвод насыщенного раствора. В настоящее время существует множество типов сравнительно недорогих реакторов. В зависимости от условий реакции существует два класса реакторов для перемешивания, которые делятся по типу осуществляемого перемешивания: механические и пневматические. [20] В механических используется агитатор различной конструкции, приводимый в движение двигателем (Рисунок 1.4 слева). По сравнению с пневматическими (Рисунок 1.4 справа), механические требуют больших капитальных и операционных затрат. В пневматических реакторах переме�

Информация о работе

Градов, Дмитрий Владимирович
кандидата технических наук
Санкт-Петербург, 2013
ВАК 25.00.13

Диссертация

Разработка эффективной конструкции газодинамического реактора на основе математического моделирования - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы