Integrating Multidisciplinary Modeling Tools to Support Upstream Mine Waste Management

Le drainage minier acide (DMA), suivant l'oxydation des sulfures déclenchée par l'oxygène atmosphérique et l'eau libre, occasionne une détérioration de la qualité de l'eau en terme d'acidité et de concentration des métaux et oxyanions. Le DMA constitue un danger environnemental d'ampleur mondiale vu son effet néfaste sur les systèmes aquatiques et les formes de vies fauniques et floristiques. Les contributions scientifiques menées par le gouvernement, l'industrie minière, les universités et les établissements de recherche se concentrent sur l'évaluation, la prévention et le traitement du DMA pour préserver les écosystèmes avoisinant les installations minières. La plupart des contributions scientifiques abordent les aspects avals de la gestion du DMA impliquant les étapes opérationnelles et post-fermeture du cycle minier. Peu de solutions pratiques ont été suggérées pour la phase de développement en raison du manque de rejets solides in-situ et les données expérimentales nécessaires pour envisager lesdites contributions scientifiques. Par conséquent, le concept de gestion en amont a reçu très peu d'attention. De même, les approches de modélisation utilisées pour prévoir le DMA sont largement étudiées. Bien qu'elles offrent de nombreux avantages, la majorité de ces études de modélisation sont réalisées pendant les phases d'exploitation et de fermeture du cycle minier, car elles abordent les stratégies de conception et la performance des scénarios de restauration. En outre, la recherche scientifique basée sur des approches interdisciplinaires pour atténuer le risque environnemental du DMA devrait être davantage mise en évidence et développée. La revue de la littérature souligne trois concepts principaux ; la géométallurgie, le principe de conception pour la fermeture et la réflexion en amont. La géométallurgie a été principalement développée pour dissoudre les limites interdisciplinaires entre le géologue, le métallurgiste et l'ingénieur minier afin d'optimiser les profits économiques et d'atténuer les risques techniques. Plus récemment, divers chercheurs ont inclus les enjeux environnementaux miniers dans la réflexion holistique de la géométallurgie. Le principe de conception pour la fermeture exige que les problèmes environnementaux potentiels soient pris en compte et planifiés avant et pendant les étapes de l'exploitation minière. De même, la réflexion en amont propose d'introduire des pratiques préventives dans les filières de gestion des rejets miniers. Ces pratiques doivent être entreprises dès les premières étapes possibles du cycle minier, appelées étapes en amont. Bien que de nombreuses contributions scientifiques aient abordé les concepts susmentionnés, elles se sont principalement concentrées sur des expériences en laboratoire et/ou sur le terrain. Aucune recherche n'a abordé la réconciliation entre les approches de modélisation pour soutenir l'atténuation des risques environnementaux du DMA. Le lien entre le géologue et le géochimiste environnemental est une préoccupation croissante qui devrait être abordée pour fournir au gestionnaire des rejets miniers de nouvelles options pour aller au-delà des méthodes de gestion conventionnelles. L'objectif principal du présent travail est de dissoudre les limites interdisciplinaires entre les approches de modélisation pertinentes pour améliorer la gestion en amont des rejets miniers solides et la prévention du DMA. Ainsi, trois objectifs principaux ont été définis : (1) lier la modélisation géologique et les attributs environnementaux du DMA pour effectuer une classification spatiale proactive des rejets miniers en fonction de leur risque environnemental inhérent; (2) réconcilier les outils de modélisation cinétique avec les contraintes des stades amonts du cycle minier (ex. le stade de développement) et (3) intégrer la composante spatiale développée selon le premier objectif et la dimension temporelle de la modélisation cinétique pour concevoir une approche de modélisation holistique permettant la classification en amont des rejets miniers et supportant l'atténuation des risques environnementaux. Pour atteindre ces objectifs, l'approche méthodologique a consisté dans un premier temps à lier les informations géologiques, collectées tout au long des campagnes du logging géologique, à la modélisation numérique. Ce lien a été établi par une approche stochastique qui relie les variables discrètes et continues du logging géologique. Le résultat de la simulation stochastique soutient la modélisation géologique 3D en accomplissant la densité spatiale adéquate des données numériques. Cette partie permettait d'établir des modèles numériques 3D décrivant la distribution spatiale d'un contaminant donné contenu dans la roche hôte. Par la suite, l'approche de modélisation cinétique a été réalisée pour simuler le pH résultant des principales réactions de génération et de neutralisation de l'acidité. Le modèle cinétique prend en compte des conditions hautement oxydantes et des réactions contrôlées par la réactivité minérale et la surface disponible. Par conséquent, la diffusion d'oxygène n'a pas été considérée comme l'étape limitant le processus d'oxydation-neutralisation. Le modèle a été calibré et étalonné par rapport à des tests cinétiques expérimentaux dont les conditions opératoires sont conformes à l'hypothèse du modèle. Enfin, les modèles spatiaux et cinétiques susmentionnés ont été intégrés pour permettre une classification dynamique de la roche abritant le minerai. Le modèle spatio-temporel intégré implique le logging géologique, la simulation stochastique, la modélisation géologique 3D, la modélisation cinétique et la modélisation de l'écoulement non saturé. Cette approche holistique décrit la distribution spatiale des principaux minéraux générateurs et neutralisants d'acide et entreprend une modélisation du transport réactif 1D pour chaque constituant volumétrique élémentaire, nommé voxel. Par la suite, une classification en amont des rejets miniers pourrait être effectuée en fonction de la teneur inhérente d'un contaminant donné dans la roche hôte et en fonction du pH qui pourrait être libéré à l'emplacement X, Y, Z du corps minéralisé si l'assemblage minéral correspondant est soumis à des conditions oxydantes. Les résultats de chaque axe sont résumés ci-dessous. L'utilisation de la modélisation géologique 3D pour la gestion des rejets miniers a permis de visualiser la distribution spatiale des contaminants dans un corps minéralisé et sa roche hôte. Par la suite, les responsables en matière de gestion pourraient facilement entreprendre la classification des stériles. À cet égard, le site minier Éléonore a fourni une base de données restreinte des teneurs en arsenic, l'élément le plus délétère dans son environnement minier, pour créer un modèle spatial 3D de la teneur en arsenic. Leapfrog Geo a été utilisé pour effectuer la modélisation géologique 3D et le logiciel de modélisation géostatistique de Stanford (SGeMS) a été utilisé pour entreprendre l'analyse du variogramme spatial. Le résultat de ce travail consiste en un modèle spatial 3D en multi-réalisation de la teneur en arsenic à travers le gisement et la roche encaissante. Chaque réalisation a été évaluée à l'aide des analyses chimiques mesurées pour souligner la fiabilité du modèle. Les résultats ont révélé un vaste halo géochimique d'arsenic qui s'étend jusqu'à 500 m du gisement d'or, avec jusqu'à 94 % des teneurs en arsenic dépassant 50 ppm. Les résultats de la modélisation cinétique à l'aide de PHREEQC ont montré un bon accord avec les données des minicellules d'altération. L'objectif principal de simulation du pH à l'aide de la modélisation cinétique des essais en minicellules d'altération a été atteint. Cependant, le modèle n'inclut pas les processus de rétention géochimique tels que la coprécipitation et la sorption. Étant conscient de ces limites, le modèle cinétique PHREEQC n'est pas conforme aux objectifs de conception liés à la restauration minière. Cependant, il est conforme aux études préliminaires tout au long de la phase de développement, qui a peu bénéficié des outils de modélisation géochimique. Le principal atout du modèle cinétique proposé est la capacité d'entreprendre une analyse paramétrique pour l'identification en amont des risques en se basant sur une base de données restreinte et un raisonnement de modélisation conservateur. À cet égard, les données d'entrée sont constituées de la caractérisation minéralogique habituelle, les tests de mini-cellule d'altération et les taux de réactivité minérale tirés de la littérature, respectant ainsi deux contraintes principales qui orientent l'étape de développement : la disponibilité des matériaux et le coût d'évaluation. Les résultats du modèle intégrant la modélisation spatiale et la modélisation du transport réactif montrent l'évolution spatio-temporelle du pH tout au long du plan central d'exploitation minière. À cet égard, 527 simulations de transport réactif ont été effectuées tout au long du plan minier composé de voxels de 40х40х40 mètres. Le géomodèle spatio-temporel met en évidence l'effet de la réactivité des minéraux neutralisants sur le pH lors de l'oxydation des sulfures. L'intégration de la géologie et la géochimie environnementale est la solution clé pour des opportunités de production plus propres. Le présent projet a fait progresser les connaissances sur la classification en amont des rejets miniers et a mis en place des méthodes prometteuses pour intégrer des approches de modélisation multidisciplinaires en vue d'un meilleur contrôle des rejets solides dans les mines métalliques. Les méthodes appliquées ici ne se limitent pas aux études de cas abordées et pourraient être appliquées à d'autres projets miniers atteignant les stades de développement et/ou d'exploitation.