Моделирование потока суспензий в шахтных плотинах для анализа стабильности с использованием Ansys Fluent

Введение в стабильность плотин и поток суспензии

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

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

Понимание потока суспензии в плотинах

Сыпь — это суспензии мелких твердых частиц (обычно < 200 мкм) в воде, часто с высокой концентрацией твердого вещества (40% — 60% по массе). Поведение потока этих суспензий регулируется распределением размеров частиц, минералогией и химией воды. При высоких концентрациях суспензии проявляют неньютоновские характеристики: они могут быть тиксотропными (вязкость, зависящая от времени), тоньше сдвига (вязкость уменьшается со скоростью сдвига) или иметь напряжение выхода, которое должно быть преодолено до начала потока. Эти свойства значительно влияют на то, как суспензия движется через плотину, как она консолидируется и как развиваются давления пор.

Реологические модели для хвостохранилищ

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

Бингемский пластик — линейное отношение «сдвига-стресс» и «сдвига-скорость» после превышения предела урожаев. Подходит для высококонцентрационных суспензий, которые ведут себя как вязкопластическая жидкость.
Herschel-Bulkley — более обобщенная модель, которая учитывает утоление сдвига или утолщение сдвига после выхода. Широко используется для хвостохранилищ, поскольку она хорошо подходит для экспериментальных данных в широком диапазоне скоростей сдвига.
Закон власти (FLT:0) — полезен, когда стресс урожая незначительный, но часто неадекватный для оседлых или высококонцентрированных шламов.
Кассон (FLT:0) — исторически использовался для грязи и добычи суспензий, хотя сегодня он встречается реже.

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

Урегулирование и консолидация

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

Роль CFD в анализе стабильности плотин хвостохранилища

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

Многофазные и пористые возможности ANSYS Fluent позволяют рассматривать хвостохранилища как жидкостно-твердую смесь, перемещающуюся через деформируемую пористую среду (консолидированные хвостохранилища и структура плотины). Связывание CFD с геотехническим анализом конечных элементов (например, посредством одностороннего или двустороннего обмена данными) обеспечивает еще больше понимания, но один только Fluent уже может выводить поля давления, профили скорости и местоположения свободной поверхности, которые подпитывают обычные расчеты стабильности.

Методы моделирования в ANSYS Fluent

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

Многофазные модели потока

Хвосты по своей сути многофазны. ANSYS Fluent предлагает несколько вариантов:

Объем модели жидкости (VOF) — лучше всего при отслеживании четкого интерфейса между суспензией и воздухом (или водой), например, во время пляжного отдыха или в декантном пруду. VOF захватывает динамику свободной поверхности, но не моделирует отдельные фазы частиц; суспензия рассматривается как одна жидкость с определяемой пользователем реологией.
Эвлерийско-Эвлеровская модель — рассматривает каждую фазу (воду, твердые тела) как взаимопроникающий континуум со своей собственной скоростью, температурой и объемной фракцией. Это наиболее строгий подход для плотных суспензий, где взаимодействия частиц и частиц значительны. Он требует моделей закрытия для сопротивления, подъема и гранулированного давления (например, кинетическая теория гранулярного потока). Он вычислительно дорог, но предоставляет подробную информацию о распределении твердых тел.
Модель смеси — упрощенная евлеровская модель, которая решает уравнение одного импульса для смеси и использует алгебраические отношения для скоростей скольжения между фазами.Подходящие при оседании скорости хорошо характеризуются и фракция твердых тел умеренная.
Discrete Phase Model (DPM) – рассматривает частицы как отдельные лагранжевые элементы в непрерывной жидкости. Подходит для загрузки твердых частиц с низким объемом (<10%), но непрактичен для плотных хвостов из-за вычислительной стоимости.

Для типичных высокотвердых хвостохранилищ рекомендуется подход Эйлера-Эулера или модель смеси (с соответствующими затворами перетаскивания и оседания).Модель VOF может использоваться в сочетании с одной из этих моделей для отслеживания свободной поверхности.

Моделирование неньютоновской реологии

ANSYS Fluent позволяет пользователям определять неньютоновскую вязкость с помощью функций, определенных пользователем (UDF), или путем выбора одной из встроенных моделей (Bingham, Herschel-Bulkley и т. Д.).

Выберите соответствующую реологическую модель на основе экспериментальных данных.Сохранить параметры (уровень напряжения, индекс консистенции, индекс потока) в базе данных материалов.
Если фракция твердых веществ изменяется пространственно, используйте UDF, чтобы сделать реологические параметры зависимыми от локальной объемной фракции твердых веществ. Это имеет решающее значение, поскольку суспензия вблизи точки разряда менее концентрирована, чем оседлый слой.
Включите динамические сетчатые или пористые медиа-зоны для регионов, которые консолидировались в твердый пирог.

Распространенной ловушкой является использование ньютоновского приближения вязкости, которое полностью пропускает поведение «доходность-стресс». Небольшой стресс (например, 5 Па) может резко повлиять на рассеивание порового давления и распространение потока.

Пористые СМИ и консолидация

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

Шаги в процессе моделирования

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

Геометрия и сетчатая генерация

Геометрия хвостохранилища включает набережную плотины, дно бассейна, неплотный пруд и трубопровод или спигот. Во многих случаях плотина со временем растет; 3D-модель может включать поэтапную конструкцию с использованием нескольких областей или динамического движения сетки. Сетка должна захватывать критические особенности:

Свободная поверхность (обычно рафинированная вблизи суспензионного воздушного интерфейса).
Зона сброса суспензии (высокоскоростные градиенты).
Внутренняя часть плотины и фундамент (где вычисляется давление пор).

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

Граничные условия и материальные свойства

Ключевые граничные условия включают:

Впуск — скорость поступления суспензии, фракция твердого объема, реологические параметры (или UDF для их изменения).
Пустошь — декантная поверхность пруда или сток пальца ноги.Обычно устанавливается как выходное отверстие давления с неподвижной гидростатической головкой.
Стены — поверхность плотины и дно бассейна: нет-скальзывание подходит для жидких фаз, при этом напряжение сдвига вычисляется из суспензионной реологии.
Симметрия — используется, если геометрия плотины симметрична.

Свойства материала должны быть определены для каждой фазы. Для фазы твердых тел в модели Эйлера, укажите плотность, диаметр частиц и гранулированную вязкость модели. Жидкая фаза (вода) использует стандартные свойства. Не забудьте указать соответствующий размер шага времени: явные схемы могут потребовать очень небольших шагов времени (например, 0,001 с) для поддержания стабильности в высокодоходных напряженных областях.

Настройки Solver и моделирование

Для временного моделирования осаждения и консолидации суспензии:

Используйте растворитель на основе давления с алгоритмом PISO или SIMPLE для связи скоростей давления.
Включите уравнение объемной фракции для модели Эйлера или установите явную схему VOF при отслеживании свободной поверхности.
Выберите пространственную дискретизацию второго порядка для импульса и объемной фракции, чтобы уменьшить численную диффузию.
Установите факторы пониженной релаксации тщательно — потоки стресса могут быть чувствительны к факторам релаксации ниже 0,3.
Инициировать домен с небольшим объемом воздуха над исходным хвостовым ложем. Используйте участки пластыря для определения начальной высоты суспензии.

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

Последующее производство и валидация

После моделирования инструменты постобработки ANSYS Fluent могут визуализировать:

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

Валидация необходима. Сравните вычисленные поровые давления на различных глубинах с показаниями пьезометра с фактической плотины. Если моделирование соответствует полевым данным в пределах приемлемой допуска (например, ±15%), модель может быть использована для прогнозирования будущего поведения или для тестирования модификаций конструкции.

Практические соображения для анализа стабильности

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

Поле порового давления — непосредственно поглощается в программном обеспечении с предельным равновесием (например, SLIDE, GeoStudio) через карту давления. Высокие поровые давления в массе хвостов снижают эффективную прочность и повышают риск статического сжижения.
Распределение напряжения сдвига — Fluent может вычислить напряжение сдвига, создаваемое суспензией на полу плотины и стенах. Сравнение их с прочностью сдвига материалов фундамента указывает, где вероятна неисправность.
Свободное расположение поверхности — высота поверхности фреатики в пределах плотины является классическим входом для диаграмм устойчивости. Чрезмерно высокая поверхность фреатики предполагает недостаточный дренаж.
Пути потока суспензии — векторы скорости могут выявить внутреннюю эрозию (трубопровод), если зоны высокой скорости находятся вблизи ядра плотины или основания.

Инженеры часто комбинируют результаты CFD с анализом напряжений конечных элементов (например, с использованием Ansys Mechanical) для вычисления коэффициента безопасности от скольжения и отказа вращения. Этот комбинированный подход обеспечивает более полную картину, чем любой из методов в одиночку.

Тематические исследования и приложения

Несколько горнодобывающих предприятий использовали ANSYS Fluent для анализа управления хвостохранилищами.

Медный рудник в Чили смоделировал последовательное осаждение хвостохранилищ в плотине для оптимизации местоположения шиповника и снижения риска перенасыщения. Моделирование показало, что перемещение точки разряда на 20 м вверх по течению уменьшило высоту поверхности фреатики на 3 м, улучшив коэффициент безопасности с 1,2 до 1,5.
Золотая шахта в Западной Африке использовала модель Эйлера-Эулера с реологией Гершеля-Булкли для оценки скорости рассеивания порового давления в плотностной хвостохранилище.Результаты привели к установке дополнительных горизонтальных дренажных слоев, что сократило время консолидации с 18 месяцев до 10 месяцев.

Хотя детали являются собственностью, такие примеры демонстрируют практическую ценность CFD в проектировании хвостохранилищ.

Ограничения и вызовы

Не существует идеального моделирования. Основные ограничения описанного подхода включают:

Расчетная стоимость — моделирование с использованием моделей с высоким разрешением Эйлера-Эулера с DEM-связью может занять несколько дней на многоядерной рабочей станции. Масштабирование до полномасштабных 3D переходных моделей может потребовать кластерных вычислений.
Реологическая неопределенность — свойства хвостов изменяются с содержанием минералов, формой частиц и химией воды. Моделирование так же хорошо, как и входные данные реологии, которые могут быть не репрезентативными для всех фаз осаждения.
Упрощенные геотехнические конститутивные законы — Fluent изначально не моделирует стресс-деформацию структуры плотины.Связь CFD с геотехническим FEA по-прежнему является специализированной задачей, требующей тщательного картирования данных.
Данные проверки — многим шахтам не хватает комплексного приборостроения (пьезометров, инклинометров) для сравнения. Без валидации моделирование остаётся прогностическим инструментом с неизвестной точностью.

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

Заключение

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

Для дальнейшего чтения о базовой физике и лучших практиках, проконсультируйтесь с документацией Ansys Fluent , Международным советом по горному делу и металлам (ICMM) руководящими принципами управления хвостохранилищами и Ресурсы хвостохранилища Geoengineer.org . Дополнительную информацию о реологии шламов можно найти через Приложения для добычи шлама RheoSense .