Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЛАЖНОСТИ МАССИВА ПРИ КАПИЛЛЯРНОМ ДВИЖЕНИИ РАСТВОРОВ

Харитонова М.Ю. 1 Михайлов А.Г. 1 Вашлаев И.И. 1
1 Институт химии и химической технологии СО РАН
В статье приведены результаты лабораторного эксперимента и теоретических расчетов по изменению влажности массива пород. В качестве объекта исследований использованы лежалые хвосты обогащения сульфидных медно-никелевых руд Норильского комбината, складированные в долине реки Щучья. Установлена возможность повышения насыщенности породы раствором до требуемой влажности путем чередования горизонтальных слоев с различными капиллярными свойствами. Используя горизонтальные слои с меньшими значениями крупности материалов, не прикладывая дополнительного внешнего давления, можно поднять раствор к поверхности породы и получить необходимую влажность на поверхности. Экспериментально доказана возможность извлечения раствора с поверхности породы с использованием гигроскопического материала, размещенного на поверхности. Представлена математическая модель процесса инфильтрационного увлажнения техногенного массива, сложенного горизонтальными слоями различной крупности. С использованием разработанного алгоритма проведены расчеты влажности породы, сложенной слоями различной толщины и крупности.
хвосты обогащения
влажность
капиллярное поднятие
1. Зосин А.П., Приймак Т.И., Кошкина Л.Б. Экологические аспекты процессов геохимической трансформации минеральных отходов от переработки сульфидных медно-никелевых руд // Экологическая химия. – 2003. – № 12(1). – С. 33–40.
2. Полубаринова-Кочина П.Я. Теория движения грунтовых вод. – М.: Наука, 1977. – 664 с.
3. Чантурия В.А., Макаров Д.В., Макаров В.Н., Васильева Т.Н. Окисление нерудных и сульфидных минералов в модельных экспериментах и в реальных хвостохранилищах // Горный журнал. – 2000. – № 4. – С. 55–58.
4. Чантурия В.А., Макаров Д.В., Макаров В.Н. Изменение нерудных минералов горнопромышленных отходов в процессе хранения под воздействием минеральных кислот // Инженерная экология. – 2000. – № 1. – С. 31–39.
5. Wong H., Morris S., Radke C.J. Three dimensional menisci in polygonal capillaries // J. Colloid Interf. Sci. – 1992. – № 148 (2).– Р. 317–336.
6. Wong H., Morris S., Radke C.J. Two dimensional menisci in nonaxisymmetric capillaries // J. Colloid Interf. Sci. – 1992. – № 148 (1). –Р. 284–287.

Практически во всех хвостах переработки и отвалах некондиционных руд происходят инфильтрационные геологические процессы, вызванные тем, что минеральные отходы подвергаются воздействию атмосферных осадков, нижние слои – воздействию грунтовых вод. В результате почти постоянного притока воды происходит трансформация отходов, растворимые вещества вовлекаются в водную миграцию. Далее может происходить концентрирование растворенных веществ на геохимических барьерах или, наоборот, рассеяние далеко за границами хвостохранилищ. Движение минерализованных растворов в хвостохранилищах осуществляется в условиях двух видов: при полном водонасыщении массива или при неполном водонасыщении в зоне капиллярной каймы. В последнее время вопросы фильтрации водных растворов через отвалы и изменения минералов, происходящие при этом, являются объектом изучения многих исследователей [1, 3, 4 и др.]. Что касается вагопереноса в условиях неполного насыщения в хвостохранилище, то этот вопрос изучен недостаточно, хотя имеет первосте­пенное значение для изучения ряда важных прак­тических проблем. Среди них направленный перенос растворов в техногенных рыхлых массивах; управление уровнем насыщенности массива раствором с использованием инфильтрационных процессов; создание условий для концентрирования компонентов в отвалах; изучение режима влажности массива в связи с решением различных задач обоснования параметров технологии инфильтрационного извлечения полезных компонентов из техногенных объектов (имеется в виду подпитывание процессов направленного движения растворов по капиллярам массива потенциальными ресурсами недр).

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

Математическая модель, алгоритм

Исследуемая порода представляет собой материал с совокупностью контактирующих между собой пор различных размеров. Это материал со сложной структурой, и детальное описание его затруднено. Для характеристики реального материала, который мог бы быть описан математически, была применена несколько упрощенная модель, не меняющая основные свойства массива. Структура исследуемого массива представляет собой прямые цилиндрические капилляры, размер которых характеризуется некоторым средним значением RК, а промежутки между ними заполнены твердой фазой из сферических частиц с эффективным диаметром DЭ. При изучении породы, сложенной двумя горизонтальными слоями, рассматривали модель, представляющую собой систему контактирующих между собой капилляров различных радиусов RК1 и RК2.

Для такой модели правильной упаковки сферических частиц одинакового радиуса получены значения эффективного диаметра частицы Dэ = 0,0634 мм, радиуса капилляров Rк = 0,006 мм. Эти величины приближенно отражают реальную структуру исследуемого материала и использованы в дальнейших расчетах.

Рассмотрим массив, сложенный горизонтальными слоями с разным размером частиц и капилляров. Нижний слой массива контактирует с горизонтом раствора, соответствующим горизонту подземных вод. До уровня дневной поверхности слой массива имеет разные степени влажности. Эта зона массива между горизонтом подземных вод и дневной поверхностью представляет интерес в оценке возможности движения растворов. Такая оценка может быть представлена характеристикой степени влажности. Для разных значений диаметров капилляров и крупности частиц значения влажности будут разными. Кроме того, материал разной крупности может быть уложен слоями, как горизонтальными, так и наклонными. Задачу оценки влажности породы, сложенной горизонтальными слоями с разными капиллярными свойствами, представили следующим образом. Предположим, для слоя массива 1, характеризующегося набором основных параметров (размер пор, диаметр зерен, пористость, высота капиллярного поднятия – RК1, DЭ1, m1, hК1), известно распределение влажности W1(X) в зоне капиллярной каймы на высоте Х от уровня гравитационной зоны. Для слоя 2, характеризующегося набором параметров (RК2, DЭ2, m2, hК2), известно распределение влажности W2(X) в зоне капиллярной каймы на высоте Х от уровня гравитационной зоны. Слои массива уложены горизонтально, толщина слоев А1 и А2. Требуется определить влажность слоя 2 на высоте X = (А1 + А2) при различных значениях А1 и А2. При такой постановке задачи необходимо оценить возможность поднятия раствора в слое 1 до границы раздела двух слоев, возможность его капиллярного поднятия во втором слое, возможность поднятия раствора до свободной поверхности и значение влажности на поверхности породы.

Построение решения основывалось на следующих положениях, сформированных известными учеными. Распределение влажности в капиллярной кайме на высоте X от гравитационной зоны может быть определено по эмпирической формуле С.Ф. Аверьянова [2]:

Eqn157.wmf

где m – пористость; w0 – наименьшая влагоемкость; hК – максимальная высота капиллярного поднятия, ось Х направлена по вертикали вверх.

Вода в вертикальном капилляре поднимается под действием капиллярных сил, вызванных силами поверхностного натяжения на границе контакта породы и воды (уравнение Жюрена). Отсюда величина капиллярного давления (ΔpK):

Eqn158.wmf (1)

где σ – коэффициент поверхностного натяжения.

Жидкость будет всасываться в капилляр до тех пор, пока капиллярное давление не уравновесится обратным воздействием гидростатического давления (рг) столба жидкости высотой h (уравнение Лапласа):

Eqn159.wmf

где ρw – плотность воды; q – ускорение свободного падения.

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

Eqn160.wmf (2)

Однако, как показали Wong и др. [5,6], для капилляров различных форм поперечного сечения численное значение величины а зависит от угла контакта жидкости с капиллярной стенкой и ограничено уравнениями (1) и (2).

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

Результаты исследования и их обсуждение

С использованием предложенного алгоритма (рис. 1) проведены расчеты влажности породы, сложенной слоями различной толщины и крупности. Расчеты проводились для реального материала – хвостов обогащения сульфидных медно-никелевых руд, складированных в долине р. Щучья, которые являются отходами горно-металлургического комбината Норильского промрайона. Установлена возможность повышения влажности породы, если она искусственно отсортирована на слои с различным гранулометрическим составом и слои уложены от самого крупного класса крупности к самому мелкому классу. При такой упорядоченной укладке слоев при одинаковой толщине фильтрующего слоя 84 см влажность на поверхности возрастает с 32,5 % до 33,8 % (рис. 2). Увеличение на первый взгляд незначительное, но при увеличении толщины фильтрующего слоя можно добиться более существенного возрастания влажности. Так, например для исследуемого материала с крупностью частиц Dэ = 0,063 мм и толщине слоя 3 м выщелачивающий раствор под действием капиллярных сил не поднимется до поверхности. Влажность породы порядка 7 % будет находиться на высоте примерно 2,5 м от уровня грунтовых вод. Если расположить 2 м исследуемого материала и на него сверху 1 м породы, состоящей из более мелких частиц (Dэ = 0,059 мм), то при общей толщине слоя 3 м влажность на поверхности составит 7 %.

pic_75.tif

Рис. 1. Алгоритм определения влажности на поверхности породы, сложенной двумя горизонтальными слоями с разными свойствами капиллярности

pic_76.wmf

Рис. 2. Влажность в зависимости от высоты относительно гравитационной зоны для пород различной крупности: 1 – слой с эффективным диаметром частиц dэ = 0,063 мм; 2 – упорядоченное расположение слоев от dэ = 1,5 мм до dэ = 0,03 мм; 3 – слой с эффективным диаметром частиц dэ = 0,03 мм

С использованием только капиллярных свойств породы в приведенном примере раствор поднимается на 0,5 м. При увеличении толщины мелкого слоя до 2 м и снижении толщины крупного слоя до 1 м влажность на поверхности составит уже 18 %. Таким образом, установлено, что, используя горизонтальные слои с меньшими значениями крупности материалов, не прикладывая дополнительного внешнего давления, можно поднять раствор к поверхности породы и получить необходимую влажность на поверхности.

Результаты теоретических расчетов были проверены на лабораторном эксперименте. Отходы обогащения руд объемом 5500 см3 помещались в цилиндрическую колонку диаметром 10 см и высотой 100 см. Снизу в колонку подавался раствор. На поверхности породы разместили слой гигроскопического материала толщиной 3, со средним радиусом капилляров 0,003 мм, образно говоря, «насос». Применяя такой «насос», в лабораторном эксперименте удалось получить более высокую влажность исследуемой породы при одной и той же высоте всасывания. В условиях проводимого эксперимента гигроскопический верхний слой удерживал до 180 г раствора.

Выводы

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

Рецензенты:

Гилев А.В., д.т.н., профессор кафедры горных машин и комплексов, ФГАОУ ВПО «Сибирский федеральный университет», г. Красноярск;

Брагин В.И., д.т.н., профессор кафедры обогащения полезных ископаемых, ФГАОУ ВПО «Сибирский федеральный университет», г. Красноярск.

Работа поступила в редакцию 11.04.2013.


Библиографическая ссылка

Харитонова М.Ю., Михайлов А.Г., Вашлаев И.И. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЛАЖНОСТИ МАССИВА ПРИ КАПИЛЛЯРНОМ ДВИЖЕНИИ РАСТВОРОВ // Фундаментальные исследования. – 2013. – № 6-3. – С. 634-637;
URL: https://fundamental-research.ru/ru/article/view?id=31562 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674