CC BY
32
© A.C. Голиков, В.Г. Наумснко, С.О. Фсдосссва, Е.И. Назимко, 2013
УДК 622.794.2:622.765
А.С. Голиков, В.Г. Науменко, С.О. Федосеева, Е.И. Назимко
МОДЕЛИРОВАНИЕ РАБОТЫ ВОДНО-ШЛАМОВОЙ СИСТЕМЫ, ПРОЦЕССОВ ОБЕЗВОЖИВАНИЯ И ФЛОТАЦИИ ПРИ ОБОГАЩЕНИИ УГОЛЬНЫХ ШЛАМОВ
Представлены результаты моделирования процессов накопления шламов в замкнутых технологических циклах углеобогатительных фабрик и взаимодействия фаз при обезвоживании тонкодисперсных осадков в различном механическом состоянии осадка и флотации угольных шламов.
Ключевые слова: моделирование, уголь, накопление шламов, фильтрование, скорость, сдвиг, флотация.
Водно-шламовая система
(ВШС) углеобогатительных фабрик представляет собой сложный комплекс, включающий классификацию частиц по крупности, осветление шламовых вод для их повторного использования, сгущение и обезвоживание шламов, обезвоживание и складирование продуктов обогащения [1]. В таких системах действует множество аппаратов разного назначения, соединенных транспортными магистралями. При этом показатели работы одного аппарата зависят от показателей работы остальных.
В связи с этим особую актуальность приобретает исследование поведения замкнутых технологических циклов, к которым относятся и ВШС. В момент запуска системы установки функционируют в нестационарном режиме, что существенно влияет на изменение сепарационных характеристик сгустителей в сторону увеличения граничного зерна разделения и среднего вероятного отклонения [2— 4]. При циркуляции потоков накапли-
вается тонкий шлам и как следствие повышается плотность и вязкость воды, что отрицательно сказывается на показателях работы гравитационного оборудования, осаждения, флотации, фильтрования.
Изменение равновесной концентрации шламов в системе и продолжительность нестационарного режима ее работы можно проследить, исследуя скорость накопления шлама в оборотной воде с учетом задержки порции шлама в каждом аппарате системы (инерционности аппарата Та) и запаздывания подачи порции шлама по трубопроводам (транспортного запаздывания Тр.), которые в сумме составляют общее запаздывание Т. Инерционность аппарата Та зависит от его конструктивных и технологических особенностей и принята по данным практики. Транспортное запаздывание поступления порций питания разных аппаратов рассчитано по формуле:
Ттр. = Ь Б 3600 п / У, (1)
ей
го
1
КП (мвс }—► гвч ]
Концентрат т
Обезвоживание
Фугат
КП
т
КП
— Т1 > Т2 — Т1 < Т2
„ 1„„„ 2„„„ 3„„„ 4„„„ 5„„„ 6„„„ 7„„„ 8„„„ 9„„„ Номер итерации (время работы системы, с)
б
Рис. 1. Водно-шламовая система в блочном изображении (а) и накопление шлама в системе (б)
Ф
„
а
где L — длина трубопровода, м; 5 — площадь поперечного сечения потока, м ; № — количество пульпы, проходящей по трубопроводу, м3/час; п — количество трубопроводов.
ВШС представляется графом, в котором каждый узел (аппарат) связан с другими узлами дугами (трубопроводами) [3]. Анализ работы систем осветления оборотных вод выполнен путем имитационного численного моделирования с применением компьютерных технологий [4]. Исходной информацией служат коэффициенты распределения шлама в узлах и потоках, определенные на основе опробования и анализа продуктов ВШС, а также запаздывание потоков и топология системы (порядок связи и взаимо действия аппаратов). Количество шлама, поступающего с рядовым углем, принято за единицу, расчеты для потоков проводятся в относительных единицах (долях). ВШС в виде укрупненных блоков и результаты моделирования накопления шлама в питании флотации представлены на рис. 1. На рис. 1, а введены следующие обозначения: 1 - поступление рядового угля, ГО - гравитационное обогащение,
РС - радиальный сгуститель, МВС -мокрая винтовая сепарация, ГВЧ -грохот высокочастотный, Ф - флотация, фильтрование и флокуляция, КП - конечные продукты.
Изменяемыми параметрами приняты транспортные задержки порций шлама Ттр. для коммуникаций, выделенных на рис. 1, а жирными стрелками.
Моделирование проводилось для случаев, когда Ттр1 > Ттр2 и Ттр1 < Ттр2, его результаты показаны в графической форме на рис. 1, б. Анализ полученных данных, показал, что путем изменения сочетания транспортных задержек потоков системы можно ускорить достижение равновесной концентрации шлама, а также достичь эффекта сглаживания приращений (скачков) его количества в потоках.
Количество таких приращений при 7тр1 > Ттр2 (жирная линия на рис. 1, б) больше, чем для Т^ < 7тр2 при близком равновесном содержании шлама 0,68-0,69 долей ед. В этом случае преобладают маршруты вывода шлама с низкими транспортными задержками, стабилизация его накопления происходит быстрее. При исполь-
Х,
Уг
У2
((ыг Яг \ -г-
V /V V / \\ \ / V \ Х^г Г N
--Ч М)
Рис. 2. Взаимодействие двух элементов (исходная схема МДЭ)
зовании маршрутов с высокой задержкой выводящих потоков равновесная концентрация достигается позже. Преобладание быстрых маршрутов вывода шлама предпочтительнее, т.к. в узлах системы происходит множество незначительных изменений его концентрации в потоках. Все сгу-стительно-осветлительные аппараты будут плавно реагировать на эти изменения, и их работа будет более устойчивой.
Полученные результаты помогут разработать технологические рекомендации к построению замкнутых циклов с низкой продолжительностью переходных процессов, что способствует повышению эффективности работы ВШС в целом.
Технологические операции обезвоживания продуктов обогащения играют важную роль не только в работе ВШС, но и всего предприятия. Одним из процессов механического обезвоживания является фильтрование, ко-
X торое исследовалось достаточно широко в различные годы. Исследование физических процессов, сопровождающих формирование определенных состояний осадков и свойств их микроструктуры, позволяет найти резервы снижения влажности осадков [5].
Взаимодействие фаз при обезвоживании тонких осадков представляет значительную трудность для изучения, т.к. является динамичным, находится под влиянием большого количества факторов и происходит в маленьком масштабе. Для численного моделирования сложных процессов, протекающих при фильтровании углесодержащих осадков, разработана компьютерная модель, позволяющая исследовать кинетику взаимодействия фаз. Модель базируется на методе дискретных элементов (МДЭ) [6-9].
В модели приняты сферические частицы радиусами Яг и Я2, которые движутся под действием сил и Г2 и вращаются в прямоугольной системе координат Х-У при действии моментов Мг и М2 (рис. 2). Движение частиц и их взаимодействие рассматривается в дискретные периоды времени. В компьютерной реализации эти периоды моделируются как циклы счета. Координаты центров тяжести частиц Хг, Уг и Х2, У2, скорости и силы рассчитываются на каждом цикле работы программы.
Работа компьютерной программы основана на фундаментальных законах — использованы второй и третий законы Ньютона и закон Гука. Досто-
верность результатов моделирования обеспечивается соблюдением динамического, кинетического и геометрического подобия.
В исходном файле задаются начальные скорости элемента вдоль осей, начальная скорость вращения, ускорения. Кроме того, вводится целый ряд коэффициентов, отражающих свойства среды и элементов: коэффициенты трения качения и скольжения, вязкостного сопротивления, сопротивления кручению, касательной и нормальной жесткости и др. [5-7].
В ходе имитационного исследования циклически выполняется перебор параметров всех элементов. Вычисляются силы, их моменты, поворачивающие элементы, определяются скорости и ускорения. Затем элементы перемещаются на новые позиции в соответствии с действующими на них силами.
Для исследования поведения многих сложных систем, к которым относятся и тонкодисперсные осадки, используют идеализированные среды в виде частиц правильной формы, обычно сферической. Пористость зернистых сред зависит не только от гранулометрического состава, но также от способа укладки частиц. Крайними случаями укладки являются кубическая (ортогональная) и гексагональная. Для исследования изменения свойств осадков при механическом воздействии на частицы рассмотрены различные варианты их упаковки. Наиболее низкую пористость имеют разноразмерные осадки с гексагональной упаковкой. В механике горных пород такие осадки относятся к переконсолидированным. Осадки с ортогональным расположением частиц рассматриваются как недоконсо-
лидированные, т.е. неуплотненные. При механическом воздействии на осадок в виде сдвига происходит переформирование структуры осадка и изменение сети каналов в его объеме. При этом тупиковые каналы становятся проводящими, и влага может удаляться более полно.
Из численного моделирования следует, что при приложении сдвига недоконсолидированные осадки уплотняются. Разноразмерные осадки, т.е. переконсолидированные, имеющие поровую структуру с более низкой пористостью по сравнению с одноразмерными, несколько разуплотняются. Установлено положительное влияние сдвига на переконсолидированные осадки, которое приводит к увеличению каналов между частицами. Необходимо определить, как влияет скорость сдвига на скорость фильтрации.
Для моделирования истечения жидкости (фильтрата), содержащей наиболее тонкие частицы, создана модель, в которой в добавление к элементам, имитирующим частицы осадка, введены блоки очень мелких частиц [7]. Выполнены исследования без наложения сдвига и со сдвигом с различной скоростью, отличающейся на несколько порядков. На основе полученных данных количество фильтрата определено по числу выделившихся из осадка мелких «частиц» жидкости для всех типов осадков в условных единицах без механического воздействия, при медленном и быстром сдвиге, а также средняя скорость фильтрации (рис. 3).
Установлено, что для недоконсо-лидированных (рыхлых) осадков скорость удаления влаги мало зависит от механического воздействия на осадок. Тем не менее, при медленном
а
б
г д е
Рис. 3. Состояние модели: а - исходное; б, в - при развитии медленного сдвига во времени 2 и 4с; г, д, - при развитии быстрого сдвига во времени 2 и 4с; е - средняя скорость фильтрации за весь период моделирования
в
сдвиге рыхлый осадок уплотняется, и скорость фильтрации снижается. Быстрый сдвиг повышает скорость удаления влаги, но она более низкая, чем при отсутствии сдвига [9]. Для плотных осадков скорость фильтрования без приложения сдвига гораздо ниже, чем у рыхлых осадков, ее максимальное значение составляет 79 усл. ед. (рис. 3, е). Приложение к плотному переконсолидированному осадку сдвига с высокой скоростью обеспечивает рост средней скорости фильтрования на 10-15 %.
Отсюда следует, что скорость фильтрации зависит от скорости приложения сдвигающих усилий и от механического состояния осадка. Быстрый сдвиг способствует снижению влажности осадков в любом первона-
чальном состоянии поровой структуры осадка.
Проведенные физические эксперименты для тонкодисперсных угле-содержащих осадков подтвердили качественные результаты, полученные при имитационном исследовании. Это позволило определить рациональный момент механического воздействия на осадок в зависимости от его состояния и типа.
Описанная выше модель (рис. 2) настроена для изучения взаимодействия воздушного пузырька и угольной частицы при различной степени окисленно-сти ее поверхности. Углехимическими исследованиями установлено отрицательное влияние окисленности на адсорбцию различных веществ на поверхности твердой фазы из-за значи-
4000 4500
Абсцисса
0,14 0,12 0,1 0,08 0,06 0,04 0,02 0
Б б
\ \
\ \ \ \ 42,03 = 0 >5х"2,9 9756 268
V V \ \
\ \ \ С
В А
6 8 10 12 14 16 18 20 22 24 Окисленность поверхности, %
Рис. 4. Траектория движения частицы в координатах пузырька (а) и зависимость времени нахождения фаз в контакте от степени окисленности поверхности угольной частицы (б)
тельного количества кислородсодержащих групп, их активного взаимодействия с водой, повышения количества пор и трещин выветривания и в итоге снижения природной гидрофобности угля. При высокой окисленности количество собирателя на поверхности недостаточно для эффективного закрепления частицы на пузырьке, вязкость слоя собирателя низкая. Поэтому именно этот параметр варьировался при имитационном исследовании.
В ходе моделирования установлено, что при высоких значениях окисленно-сти от 28 до 24 % угольная частица скользит по пузырьку, не закрепляясь на нем. На рис. 4 показано изменение траекторий движения частицы относительно пузырька при снижении окис-ленности поверхности с 17 до 13 % и времени нахождения фаз в контакте.
Траектория движения частицы в координатах пузырька начинается в точке 1 при контакте фаз (рис. 4, а). При движении по верхней полусфере пузырька частица разгоняется и при окисленности поверхности 17 % отрывается от поверхности пузырь-
ка в точке 2, отмеченной темным маркером. Снижение окисленности до 15 % способствует торможению и частица успевает подняться по нижней полусфере пузырька до точки 3, что увеличивает время нахождения фаз в контакте и вероятность удержания частицы на пузырьке. Однако снижение скорости в результате повышения вязкости поверхностных слоев оказывается недостаточным и флотокомплекс разрушается. При снижении окислен-ности до 13 % частица поднимается на нижней полусфере пузырька до точки 4, отмеченной крупным белым маркером. Затем скользит, ускоряясь, и поднимается на правой части полусферы до точки 5, отсюда опускается вниз и закрепляется на пузырьке в точке 6. Изменение траектории относительного перемещения частицы свидетельствует о том, что окисленность поверхности, равная 15 %, является критической. При этом значении частица движется по поверхности пузырька и быстро отрывается от него. Снижение
окисленности менее 15 % приводит к образованию флотационного комплекса, о чем свидетельствуют и данные рис. 4, б.
На графике можно выделить три участка: 1) при высокой окисленности поверхности от 15 до 24 % (участок АВ) время нахождения частицы на поверхности пузырька незначительное; 2) снижение окисленности от 15 до 11 % (участок ВС) способствует росту времени контакта элементов и повышению вероятности закрепления частицы; 3) при окисленности менее 11 % (участок СЭ) наиболее вероятно надежное закрепление частицы на пу-
1. Оборотное водоснабжение углеобогатительных фабрик / И.С. Благов, М.А. Борц, Б.И. Вахрамеев, Е.Г. Тресков. М.: Недра, 1980. 215 с.
2. Назимко Е.И. Совершенствование Работы систем осветления оборотных вод углеобогатительных фабрик / Е.И. Назимко, Е.Е. Гарковенко. Днепропетровск: Полиграфист, 2000. 174 с.
3. Nazimko L.I. A Simulation of Slime Circulation and the Effect of Circuit Design / L.I. Nazimko, V.V. Nazimko // Coal preparation. 1996. V 17. № 3-4. Р. 215-232.
4. Назимко Е.И. Исследование работы водно-шламовых систем углеобогатительных фабрик в нестационарном режиме / Е.И. Назимко, А.С. Голиков, К.А. Ёифенко // Вюник Криворiзького техтчного утвер-ситету. 2005. Вып. 6. С. 158-161.
5. Особенности флотации и обезвоживания тонкодисперсных углесодержащих материалов / Е.Е. Гарковенко, Е.И. Назим-ко, А.И. Самойлов, Ю.Л. Папушин. Донецк: Норд-Пресс. 2002. 266 с.
6. Cundall P.A. A discrete numerical Model for granular assemblies / P.A. Cundall,
КОРОТКО ОБ АВТОРАХ -
зырьке и длительное существование комплекса.
Имитационное исследование подчеркнуло необходимость применения наиболее активных реагентов и создания специальных конструкций флотомашин для возможности сепарации трудно флотируемых угольных частиц [10].
Таким образом, современный метод численного моделирования позволяет исследовать взаимодействие фаз при обогащении полезных ископаемых на микроуровне и находить рациональные условия ведения процессов разделения.
- СПИСОК ЛИТЕРАТУРЫ
0.D.L Strack // Geotechnique. 1974. 29, N
1, р. 47-65.
7. Гарковенко Е.Е. Моделирование процесса вакуумного фильтрования угольных шламов // Горный информационно-аналитический бюллетень. Москва. 2004. №1. — С. 317-319.
8. Nazimko L.I. & al. Kinetics of Phases Interaction during Mineral Processing Simulation / L.I. Nazimko, E.E. Garkovenko, A.N. Corchevsky & al // Proceedings of XV International Congress of Coal Preparation. China. 2006. P. 785-798.
9. Назимко Е.И. Моделирование процесса удаления влаги из углесодержащих осадков / Е.И. Назимко, В.Г. Науменко // Збагачення корисних копалин. Днепропетровск. 2011. Вып. 45(86). С. 153—158.
10. Kurchenko I.P. New Coal Flotation Technology Regulated by the Intake Degree of the Air and Reagents / I.P. Kurchenko, Yu.N. Filippenko, G.A. Mavrenko, S.O. Fe-doseeva // Proceedings of XVI International Congress of Coal Preparation. USA. 2010. P. 471-478. ШШ
Голиков Алексей Сергеевич - ассистент, соискатель, Науменко Виктория Георгиевна — ассистент, соискатель,
Назимко Елена Ивановна - доктор технических наук, профессор, зав. кафедрой, Донецкий национальный технический университет, [email protected],
Федосеева Светлана Олеговна - научный сотрудник института УкрНИИуглеобогащение, аспирант ДонНТУ.
