Introduction à la stabilité du barrage de Tailings et au débit de boue

Les barrages miniers de résidus miniers sont parmi les plus importants ouvrages d'ingénierie sur Terre, conçus pour stocker les sous-produits de déchets à grains fins (déchets de queue) produits lors de l'extraction minière.Ces réservoirs doivent rester stables pendant des décennies, souvent dans des conditions environnementales difficiles et des charges sismiques.

Les évaluations de stabilité géotechnique traditionnelles reposent sur des méthodes d'équilibre limite et des hypothèses simplifiées de pression interstitielle. Bien que ces approches soient utiles, elles ne permettent pas souvent de saisir le comportement complexe et dépendant du temps du lisier à l'intérieur du barrage. Les boues sont des fluides multiphasés non néotoniens qui évoluent au fil du temps lorsque les particules s'installent, se consolident et se dessèchent. La dynamique des fluides computationnels (CFD) offre un moyen de simuler ces processus avec une grande fidélité.

Comprendre le débit de boue dans les barrages de résidus

Les boues sont des suspensions de particules solides fines (habituellement < 200 μm) dans l'eau, souvent à forte concentration de solides (40%–60% en masse). Le comportement de ces suspensions est régi par la distribution de la taille des particules, la minéralogie et la chimie de l'eau. À de fortes concentrations, les boues présentent des caractéristiques non néotoniennes : elles peuvent être thiotropiques (viscosité dépendante du temps), cisaillement-finage (viscosité diminue avec le taux de cisaillement), ou avoir un stress de rendement qui doit être surmonté avant le début du débit.

Modèles de rhologie pour les tailings

Pour modéliser avec précision le débit de lisier dans l'ANSYS Fluent, les ingénieurs doivent définir un modèle rhéologique approprié.

  • Plastice de bingham – une relation de résistance linéaire au cisaillement par rapport au taux de cisaillement après dépassement d'une contrainte de rendement.
  • Herschel‐Bulkley – un modèle plus généralisé qui explique le cheveu de cisaillement ou d'épaississement après rendement. Largement utilisé pour les résidus parce qu'il s'adapte bien aux données expérimentales sur une large gamme de vitesses de cisaillement.
  • Loi sur la puissance – utile lorsque le stress de rendement est négligeable, mais souvent insuffisant pour les lisiers stables ou fortement concentrés.
  • Casson – utilisé historiquement pour la boue et les boues minières, bien que moins courant aujourd'hui.

La caractérisation expérimentale (p. ex., à l'aide d'un rhéomètre rotationnel) est essentielle pour obtenir les paramètres de ces modèles.

Règlement et consolidation

Au fil du temps, les particules se déposent sous la gravité, formant un lit dense et compacté au sol du barrage. Cette sédimentation modifie localement la rhéologie du lisier et affecte l'évolution de la pression interstitielle. Dans de nombreux barrages de résidus, le lisier se dépose en couches minces (plage), permettant à l'eau de décanter et à des solides de se consolider.

Le rôle du CFD dans l'analyse de stabilité du barrage de Tailings

L'analyse de stabilité d'un barrage de résidus consiste à évaluer l'équilibre entre les forces motrices (gravité, pression de suintement) et les forces résistantes (résistance au cisaillement des matériaux du barrage et des résidus eux-mêmes). La pression de l'eau de pore est un facteur critique car elle réduit la contrainte efficace, abaisse la résistance au cisaillement.

Les capacités multiphases et poreuses de Fluent ANSYS permettent de traiter les résidus comme un mélange solide de fluide se déplaçant dans un milieu poreux déformable (les résidus consolidés et la structure du barrage). L'analyse géotechnique des éléments finis (par exemple par échange de données unidirectionnel ou bidirectionnel) permet de mieux comprendre le problème, mais Fluent seul peut déjà produire des champs de pression, des profils de vitesse et des emplacements libres de surface qui alimentent les calculs de stabilité conventionnels.

Techniques de modélisation dans Anysy Fluent

La sélection du bon modèle physique dans Fluent est la pierre angulaire d'une simulation réussie. Les sous-sections suivantes décrivent les approches les plus pertinentes pour le débit de boue de barrage de résidus.

Modèles de flux multiphasés

Les boues de tailing sont intrinsèquement multiphasées.

  • Modèle de volume de liquide (VOF)[ – meilleur pour suivre une interface distincte entre le lisier et l'air (ou l'eau), p.ex., pendant le plage ou dans l'étang décanté. Le VOF capture la dynamique de surface libre mais ne modélise pas des phases de particules distinctes; le lisier est traité comme un fluide unique avec une rhéologie définie par l'utilisateur.
  • Modèle eulérien – traite chaque phase (eau, solides) comme une continuité interpénétratrice avec sa propre vitesse, température et fraction de volume. Il s'agit de l'approche la plus rigoureuse pour les boues denses où les interactions particule-particules sont importantes. Il faut des modèles de fermeture pour la traînée, le levage et la pression granulaire (p. ex., la théorie cinétique du débit granulaire).
  • Modèle de pixel – un modèle eulérien simplifié qui résout une seule équation de moment pour le mélange et utilise des relations algébriques pour les vitesses de glissement entre les phases.
  • Modèle de phase de dispersion (DPM) – traite les particules comme des éléments lagrangien individuels dans un fluide continu.

Pour les barrages à résidus à haute teneur en soufre, il est recommandé d'utiliser l'approche eulérienne ou le modèle de mélange (avec des fermetures de traînée et de décantation appropriées).

Modélisation de la rhéologie non néotonienne

ANSYS Fluent permet aux utilisateurs de définir la viscosité non néotonienne par l'intermédiaire des fonctions définies par l'utilisateur (FDU) ou en sélectionnant l'un des modèles intégrés (Bingham, Herschel-Bulkley, etc.). Pour une simulation des résidus :

  1. Choisissez le modèle rhéologique approprié basé sur des données expérimentales. Conservez les paramètres (contrainte de rendement, indice de cohérence, indice de débit) dans la base de données des matériaux.
  2. Si la fraction des solides varie spatialement, utilisez un UDF pour rendre les paramètres rhéologiques dépendant de la fraction locale de volume des solides. Ceci est critique parce que le lisier près du point de décharge est moins concentré que le lit stabilisé.
  3. Activer des zones dynamiques de mailles ou de milieux poreux pour les régions qui se sont consolidées en un gâteau de type solide.

Un piège commun est l'utilisation d'une approximation de viscosité Newtonienne, qui manque complètement le comportement de résistance au rendement. Une faible contrainte de rendement (p. ex. 5 Pa) peut affecter de façon dramatique la dissipation de la pression interstitielle et la propagation frontale du débit.

Médias poreux et consolidation

Une fois les résidus déposés, ils se comportent comme un milieu poreux avec une perméabilité décroissante au fur et à mesure que la consolidation progresse. Dans Fluent, vous pouvez définir des zones poreuses avec une perméabilité directionnelle et une fonction définie par l'utilisateur pour mettre à jour la perméabilité basée sur la fraction solide locale. Cela permet de simuler le processus de consolidation : la fraction solide initiale faible (perméabilité élevée) se transforme progressivement en un gâteau à faible perméabilité.

Étapes du processus de modélisation

Un flux de travail robuste garantit que la simulation est à la fois précise et calculable.

Géométrie et génération de mesh

La géométrie d'un barrage de résidus comprend le remblai du barrage, le plancher du bassin, l'étang décanté et le pipeline de déversement ou le spigot. Dans de nombreux cas, le barrage pousse au fil du temps; un modèle 3D peut intégrer la construction par étapes en utilisant plusieurs régions ou un mouvement dynamique du maillage.

  • Surface libre (généralement raffinée près de l'interface lisier-air).
  • Zone de décharge de boue ( gradients de vitesse élevée).
  • Intérieur et fondation du barrage (où les pressions interstitielles sont calculées).

Les mailles hybrides (tétraédrales dans des régions complexes, hexaédrales en sections droites) sont fréquentes. Une étude d'indépendance de grille devrait être effectuée pour au moins deux mailles (couronne/fine) afin de s'assurer que les pressions interstitielles et les fronts de débit calculés ne sont pas sensibles aux mailles.

Conditions de frontière et propriétés du matériau

Les principales conditions de délimitation sont les suivantes :

  • Inlet – débit d'écoulement de lisier, fraction de volume solide, paramètres rhéologiques (ou un UDF pour les varier).
  • Outlet – surface de l'étang décanté ou drain des orteils. Généralement réglé comme sortie de pression avec une tête hydrostatique fixe.
  • Tailles – face de barrage et plancher de bassin : aucun glissement n'est approprié pour les phases de fluide, avec une contrainte de cisaillement calculée à partir de la rhéologie du lisier.
  • Symmétrie – utilisée si la géométrie du barrage est symétrique.

Pour la phase des solides dans un modèle eulérien, préciser les modèles de densité, de diamètre des particules et de viscosité granulaire. La phase du fluide (eau) utilise des propriétés standard. N'oubliez pas de préciser une taille appropriée des étapes de temps : des schémas explicites peuvent nécessiter de très petites étapes de temps (p. ex. 0,001 s) pour maintenir la stabilité dans les régions à haut rendement de stress.

Paramètres du solvant et simulation

Pour les simulations transitoires de dépôts et de consolidation de lisier:

  • Utilisez le résolveur à pression avec l'algorithme PISO ou SIMPLE pour le couplage pression-vitesse.
  • Activer l'équation de fraction de volume pour le modèle eulérien, ou définir le schéma explicite de VOF si vous suivez la surface libre.
  • Sélectionnez une discrétisation spatiale du deuxième ordre vers le haut pour réduire la diffusion numérique.
  • Réglez soigneusement les facteurs de sous-relaxation – les flux de stress de rendement peuvent être sensibles aux facteurs de relaxation inférieurs à 0,3.
  • Initialiser le domaine avec un petit volume d'air au-dessus du lit de résidus initial. Utilisez des régions de patch pour définir la hauteur initiale du lisier.

Effectuer la simulation pendant suffisamment de temps physique pour saisir le cycle complet de dépôt (heures à jours de temps réel). Vérifier la convergence à chaque étape de temps : les résidus devraient diminuer d'au moins trois ordres de grandeur.

Après traitement et validation

Après la simulation, les outils post-traitement d'ANSYS Fluent , peuvent visualiser :

  • Contours de pression interstitielle et de stress efficace.
  • La fraction de volume des solides pour montrer le profil de plage et le front de consolidation.
  • Les traces vectorielles de vitesse de lisier, indiquant les zones d'écoulement préférentiel qui peuvent conduire à la tuyauterie.
  • Emplacement de la surface libre pour comparaison avec des photos aériennes ou des images de drones.

La validation est essentielle. Comparez les pressions interstitielles calculées à différentes profondeurs avec les valeurs du piézomètre obtenues à partir du barrage réel. Si la simulation correspond aux données de terrain dans des limites acceptables (p. ex. ±15%), le modèle peut être utilisé pour prédire le comportement futur ou pour tester les modifications de conception.

Considérations pratiques concernant l'analyse de stabilité

L'objectif ultime de la simulation du CFD est de fournir des données pour une analyse quantitative de la stabilité.

  • champ de pression de pore – directement ingéré dans un logiciel d'équilibre de limite (p. ex. SLIDE, GeoStudio) via une carte de pression.
  • Distribution de la contrainte de cisaillement – Fluent peut calculer la contrainte de cisaillement exercée par le lisier sur le plancher et les murs du barrage.
  • Situation de la surface libre – La hauteur de la surface phréatique à l'intérieur du barrage est une entrée classique pour les cartes de stabilité.
  • Les voies de débit lisier[ – les vecteurs de vitesse peuvent révéler une érosion interne (pipage) si des zones à grande vitesse sont trouvées près du noyau ou de la fondation du barrage.

Les ingénieurs combinent souvent les résultats du CFD avec une analyse de contrainte par éléments finis (p. ex., en utilisant Ansys Mechanical) pour calculer le facteur de sécurité contre les défaillances de glissement et de rotation.

Études de cas et applications

Plusieurs exploitations minières ont utilisé le système ANSYS Fluent pour analyser leur gestion des résidus.

  • Une mine de cuivre au Chili a modélisé le dépôt séquentiel de résidus dans un barrage de remplissage de vallée pour optimiser l'emplacement du spigot et réduire le risque de suremblayage. La simulation a montré que le déplacement du point de rejet 20 m en amont a réduit la hauteur de surface phréatique de 3 m, améliorant le facteur de sécurité de 1,2 à 1,5.
  • Une mine d'or en Afrique de l'Ouest a utilisé un modèle eulérien avec rhéologie Herschel-Bulkley pour évaluer le taux de dissipation de la pression interstitielle dans un barrage à résidus à parois épaisses. Les résultats ont conduit à l'installation de couches de drainage horizontales supplémentaires, ce qui a réduit le temps de consolidation de 18 mois à 10 mois.

Bien que les détails soient exclusifs, de tels exemples démontrent la valeur pratique du CFD dans le domaine de l'ingénierie des barrages à résidus.

Limites et défis

Aucune simulation n'est parfaite. Les principales limites de l'approche décrite sont les suivantes :

  • Coûts informatiques[ – simulations eulériennes à haute résolution avec modèles couplés DEM peuvent prendre des jours sur un poste de travail multi-cœur. L'échafaudage vers des modèles transitoires 3D à grande échelle peut nécessiter un calcul en grappes.
  • Incertitudes rhumatologiques – Les propriétés des résidus changent avec la teneur en minéraux, la forme des particules et la chimie de l'eau. La simulation n'est que aussi bonne que les données de rhéologie d'entrée, qui peuvent ne pas être représentatives de toutes les phases de dépôt.
  • Lois constitutives géotechniques simplifiées – Fluent ne modélise pas nativement le comportement de contrainte-souche de la structure de barrage solide. L'accouplement CFD avec FEA géotechnique est encore une tâche spécialisée nécessitant une cartographie minutieuse des données.
  • – de nombreuses mines ne disposent pas d'instruments complets (piézomètres, incluinomètres) pour se comparer. Sans validation, la simulation reste un outil prédictif avec une précision inconnue.

Malgré ces difficultés, l'utilisation de CFD pour la stabilité des barrages de résidus s'est accrue rapidement à mesure que la puissance de calcul augmente et que la pression réglementaire exige des évaluations plus quantitatives des risques.

Conclusion

La modélisation du débit de boues dans les barrages miniers à l'aide d'ANSYS Fluent constitue une base puissante pour l'analyse de la stabilité. En saisissant la nature multiphase du lisier, sa rhéologie non néotonienne et l'évolution de la pression interstitielle au fil du temps, les ingénieurs peuvent identifier les mécanismes de défaillance qui manquent à des méthodes plus simples. Le flux de travail — création de géométrie, maillage, conditions de limite, installation du solveur, post-traitement et validation — exige une attention particulière aux détails, mais le bénéfice est une compréhension plus approfondie du comportement des barrages et de la capacité à concevoir des réservoirs plus sûrs et plus fiables.

Pour plus de renseignements sur la physique sous-jacente et les meilleures pratiques, consulter ANSYS Documentation sur les fluides, les lignes directrices du Conseil international des mines et des métaux (ICMM) sur la gestion des résidus[ et les ressources du barrage de résidus Geoengineer.org. On peut trouver des renseignements supplémentaires sur la rhéologie des boues par l'intermédiaire RheoSense Mining slurry applications.