Introduzione ai tailings Dam Stability e Slurry Flow

Le dighe di coda di mining sono tra le più grandi strutture ingegnerizzate sulla Terra, progettate per immagazzinare i sottoprodotti di scarto finiti (codazioni) generati durante l'estrazione mineraria. Questi imposti devono rimanere stabili per decenni, spesso in condizioni ambientali dure e carichi sismici.

Le valutazioni tradizionali della stabilità geotecnica si basano su metodi di equilibrio limite e supposizioni semplificate di pressione pora. Mentre questi approcci sono utili, spesso non riescono a catturare il complesso, comportamento dipendente dal tempo della slurry all'interno della diga.

Capire il flusso di liquami in tailings dighe

Le slurie sono sospensioni di particelle solide fini (tipicamente < 200 μm) in acqua, spesso con un'alta concentrazione solida (40%–60% da massa). Il comportamento di flusso di queste sospensioni è regolato dalla distribuzione delle dimensioni delle particelle, dalla mineralogia e dalla chimica dell'acqua.

Modelli Rheologici per sartoriali

Per modellare il flusso di liquami con precisione in ANSYS Fluent, gli ingegneri devono definire un modello reologico adatto.

  • Bingham plastic[[] – una lineare traforo e taglio rapporto dopo un rendimento di stress è superato. Adatto per le slurries ad alta concentrazione che si comportano come un fluido viscoplastico.
  • Herschel‐Bulkley[[[]] – un modello più generalizzato che rappresenta un taglio-thinning o un taglio-thickening dopo la resa. Ampiamente usato per pedine perché si adatta ai dati sperimentali su una vasta gamma di velocità di taglio.
  • La legge della potenza[] – utile quando lo stress della resa è trascurabile, ma spesso inadeguato per le slurries stabili o altamente concentrate.
  • Casson] – storicamente usato per fango e miniere di doghe, anche se oggi meno comune.

La caratterizzazione sperimentale (ad esempio, utilizzando un reometro rotazionale) è essenziale per ottenere i parametri per questi modelli. Senza dati reologici accurati, le simulazioni CFD possono produrre risultati ingannevoli.

Impostazione e consolidamento

Nel tempo, le particelle si stabiliscono sotto la gravità, formando un letto denso e compatto al piano di diga. Questa sedimentazione cambia la reologia di liquami localmente e influisce sull'evoluzione della pressione di poro. In molte dighe di coda, la liquami viene depositata in strati sottili (beaching), permettendo all'acqua di decantare e solidi di consolidare.

Il ruolo del CFD in Tailings Dam analisi della stabilità

L'analisi della stabilità di una diga di coda comporta la valutazione dell'equilibrio tra forze di guida (gravità, pressione di infiltrazione) e forze di resistenza (forza di taglio dei materiali di diga e dei pedinamenti stessi). La pressione dell'acqua di poro è un fattore critico perché riduce lo stress efficace, abbassando la forza di taglio. CFD consente agli ingegneri di simulare l'evoluzione della pressione di poro durante la diga come il calcolo di slurry.

Le capacità multimediali multifase e porose di ANSYS Fluent permettono di trattare i codeing come miscela fluido-solid che si muove attraverso un mezzo poroso deformabile (il codeing consolidato e la struttura della diga).

Tecniche di modellazione in ANSYS Fluent

La scelta del modello fisico giusto in Fluent è la pietra angolare di una simulazione di successo.Le sottosezioni seguenti dettagliano gli approcci più rilevanti per il flusso di slurry di coda.

Modelli di flusso multifase

I tailing sono intrinsecamente multifase. ANSYS Fluent offre diverse opzioni:

  • Volume di Fluid (VOF) model[[[] – meglio quando si traccia una interfaccia distinto tra la liqua e l'aria (o l'acqua), ad esempio, durante la spiaggia o nello stagno decantato.
  • Il modello euleriano-eulerico[[[]] – tratta ogni fase (acqua, solidi) come interpenetrating continuo con la propria velocità, temperatura e frazione di volume. Questo è l'approccio più rigoroso per le doghe di calcolo dense dove le interazioni particella-particella sono significative.
  • Modello di miscela[[] – un modello Euleriano semplificato che risolve un'unica equazione di slancio per la miscela e utilizza relazioni algebriche per velocità di scorrimento tra le fasi.
  • Modello di fase discorrevole (DPM) – tratta le particelle come singoli elementi lagrangiani in un fluido continuo. Adatto per il carico di solidi bassi (il volume di 0,0%) ma impraticabile per i code densi a causa di costi computazionali.

Per le tipiche dighe di coda ad alto livello, è consigliato l'approccio euleriano-eulerico o il modello di miscela (con le relative chiusure di trascinamento e di fissaggio) che può essere utilizzato in combinazione con uno di questi modelli per il tracciamento della superficie libera.

Modellazione di Rheologia Non Newtoniana

ANSYS Fluent consente agli utenti di definire viscosità non newtoniana tramite Funzioni definite dall'utente (UDF) o selezionando uno dei modelli incorporati (Bingham, Herschel‐Bulkley, ecc.).

  1. Scegli il modello reologico appropriato basato sui dati sperimentali. Conserva i parametri (sforzo del bambino, indice di consistenza, indice di flusso) nel database dei materiali.
  2. Se la frazione dei solidi varia in senso spaziale, utilizzare un UDF per rendere i parametri reologici dipendono dalla frazione di volume locale dei solidi, ciò è fondamentale perché il liquame vicino al punto di scarico è meno concentrato rispetto al letto fisso.
  3. Attivare zone di rete dinamica o media porose per regioni che si sono consolidate in una torta solida.

Una trappola comune sta usando un'approssimazione della viscosità newtoniana, che manca completamente al comportamento di resistenza alla resa. Una piccola pressione di resa (ad esempio 5 Pa) può influenzare drammaticamente la dissipazione della pressione e la propagazione del fronte di flusso.

Media e consolidamento Poroso

In Fluente, è possibile definire zone porose con permeabilità direzionale e una funzione definita dall'utente per aggiornare la permeabilità basata sulla frazione solida locale. Questo consente la simulazione del processo di consolidamento: inizialmente bassa frazione solida (alta permeabilità) gradualmente transizioni a una torta di bassa permeabilità.

Passi nel processo di modellazione

Un flusso di lavoro robusto assicura che la simulazione sia accurata e computazionalmente fattibile.

Generazione di Geometria e Mesh

La geometria di una diga di coda include l'argine diga, il pavimento del bacino, lo stagno di decant, e la tubazione di scarico di liquami o spigot. In molti casi, la diga cresce nel tempo; un modello 3D può incorporare la costruzione in scena utilizzando più regioni o movimento di mesh dinamico.

  • Superficie libera (tipicamente raffinata vicino all'interfaccia a raggi ultravioletti).
  • Zona di scarico di liquami (gradienti ad alta velocità).
  • Interno e fondazione diga (dove le pressioni pore sono calcolate).

Sono comuni le mesh ibride (tetraedrali in regioni complesse, esadedrale in sezioni rette) e uno studio di indipendenza della griglia dovrebbe essere eseguito per almeno due mesh (coarse/fine) per garantire che le pressioni e i fronti di flusso del poro calcolato non siano sensibili alla mesh-sensibili.

Condizioni e proprietà dei materiali

Le condizioni principali del confine includono:

  • Inlet[] – tasso di flusso di liquami, frazione di volume solido, parametri reologici (o un UDF per variare).
  • Outlet[] – decant laghetto superficie o scarico della punta. In genere impostata come presa di pressione con una testa idrostatica fissa.
  • Walls[] – dam face e lavabo: no-slip è adatto per le fasi fluide, con lo stress della cesoia calcolato dalla reologia del liquame.
  • Simmetria[] – usato se la geometria della diga è simmetrica.

Per la fase dei solidi in un modello euleriano, specificare densità, diametro delle particelle e modelli di viscosità granulari. La fase fluida (acqua) utilizza proprietà standard. Ricordatevi di specificare una dimensione del passo appropriato: schemi espliciti possono richiedere passaggi di tempo molto piccoli (ad esempio, 0.001 s) per mantenere la stabilità nelle regioni ad alto rendimento.

Impostazioni e simulazione del solvente

Per simulazioni transitorie di deposizione e consolidamento dei liquami:

  • Utilizzare il risolutore a pressione con l'algoritmo PISO o SIMPLE per l'accoppiamento a velocità di pressione.
  • Abilitare l'equazione della frazione di volume per il modello Euleriano, o impostare lo schema esplicito VOF se il tracciamento della superficie libera.
  • Selezionare una discorsizzazione spaziale del secondo ordine per frazione di volume e di volume per ridurre la diffusione numerica.
  • Impostare attentamente i fattori di sotto-relaxazione: i flussi di soccorso possono essere sensibili ai fattori di rilassamento inferiori a 0,3.
  • Inizializzare il dominio con un piccolo volume d'aria sopra il letto iniziale di coda.

Eseguire la simulazione per un tempo fisico sufficiente per catturare il ciclo di deposizione completo (ore a giorni di tempo reale). Verificare la convergenza ad ogni punto di tempo: i residui dovrebbero cadere di almeno tre ordini di grandezza.

Post-Processing e convalida

Dopo la simulazione, gli strumenti di post-elaborazione di ANSYS Fluent possono visualizzare:

  • Contours di pressione poro e stress efficace.
  • Volume frazione di solidi per mostrare il profilo spiaggia e consolidare anteriore.
  • Trame vettoriali di velocità di liquami, indicando zone di flusso preferenziale che possono portare a tubazioni.
  • Posizione di superficie libera per il confronto con le foto aeree o immagini di droni.

Confronta le pressioni dei pori calcolati a varie profondità con le letture del piezometro dalla diga reale. Se la simulazione corrisponde ai dati del campo all'interno della tolleranza accettabile (ad esempio, ± 15%), il modello può essere utilizzato per prevedere il comportamento futuro o per testare le modifiche alla progettazione.

Considerazioni pratiche per l'analisi della stabilità

L'obiettivo finale della simulazione CFD è quello di fornire dati per un'analisi quantitativa della stabilità, particolarmente utili per i seguenti output:

  • Campo di pressione della potenza[[[] – direttamente ingerito in software di equilibrio limitato (ad esempio SLIDE, GeoStudio) tramite una mappa di pressione.
  • distribuzione dello stress[[[] – Fluente può calcolare lo stress della cinghia esercitato dalla liquami sul pavimento e sulle pareti della diga.
  • La posizione della superficie libera[ – l'altezza della superficie del faatico all'interno della diga è un classico input per le carte di stabilità.
  • Le vie di scorrimento delle liquami[] – i vettori di velocità possono rivelare l'erosione interna (piping) se le zone ad alta velocità si trovano vicino al nucleo o alla fondazione della diga.

Gli ingegneri spesso combinano i risultati CFD con un'analisi di stress a struttura finita (ad esempio, utilizzando Ansys Mechanical) per calcolare il fattore di sicurezza contro il fallimento scorrevole e rotazionale.

Studi e applicazioni dei casi

Diversi interventi minerari hanno utilizzato ANSYS Fluent per analizzare la loro gestione dei pedinamenti.

  • Una miniera di rame in Cile ha modellato la deposizione sequenziale di pedine in una diga a valle per ottimizzare la posizione della spigotta e ridurre il rischio di sovra-topping. La simulazione ha dimostrato che spostare il punto di scarico 20 m a monte ha ridotto l'altezza della superficie del faratico di 3 m, migliorando il fattore di sicurezza da 1,2 a 1,5.
  • Una miniera d'oro in Africa occidentale ha usato un modello euleriano-eulerico con la reologia di Herschel‐Bulkley per valutare il tasso di dissipazione della pore-pressura in una diga di coda densa. I risultati hanno portato all'installazione di ulteriori strati di drenaggio orizzontale, che hanno accorciato il tempo di consolidamento da 18 mesi a 10 mesi.

Mentre i dettagli sono proprietari, tali esempi dimostrano il valore pratico di CFD in ingegneria di diga di pedine.

Limitazioni e sfide

Nessuna simulazione è perfetta. Le limitazioni chiave dell'approccio descritto includono:

  • Costo computazionale[[] – simulazioni eulerian-eulerian ad alta risoluzione con modelli DEM-coupled possono richiedere giorni su una workstation multi-core.
  • Incertezza reologica[[] – le proprietà di coda cambiano con contenuto minerale, forma delle particelle e chimica dell'acqua. La simulazione è altrettanto buona quanto i dati reologici di input, che potrebbero non essere rappresentativi di tutte le fasi di deposizione.
  • Le leggi costitutive geotecniche semplificate[[] – Fluent non modella in modo nativo il comportamento dello stress-strain della struttura della diga solida.
  • Dati di valutazione[[] – molte miniere non hanno una strumentazione completa (piezometro, inclinometri) da confrontare. Senza convalida, la simulazione rimane uno strumento predittivo con accuratezza sconosciuta.

Nonostante queste sfide, l'uso di CFD per pedine stabilità di diga è cresciuto rapidamente come aumenta la potenza computazionale e come pressione normativa richiede più valutazioni di rischio quantitativo.

Conclusioni

Modellare il flusso di slitte nelle dighe di coda mineraria utilizzando ANSYS Fluent fornisce una base potente per l'analisi della stabilità. Catturando la natura multifase del liquame, la sua reologia non-Newtonian, e l'evoluzione della pressione poro nel tempo, gli ingegneri possono identificare i meccanismi di guasto che i metodi più semplici mancano. Il flusso di lavoro, la creazione di mesh, le condizioni di confine, il di processo di configurazione, i risolutorale, gli ingegneri, i dettagli

Per ulteriori informazioni sulla fisica sottostante e sulle migliori pratiche, consultare ANSYS Documentazione fluida, Consiglio internazionale sulle miniere e metalli (ICMM) Orientamenti di gestione[, e il Geoengineer.org pedine ulteriori risorse dim