Introducción a la estabilidad de la presa y a la fuga de la lotería

Las presas de colas mineras se encuentran entre las estructuras más grandes de ingeniería de la Tierra, diseñadas para almacenar los subproductos de desechos finos (detalles) generados durante la extracción de minerales. Estos implacables deben permanecer estables durante décadas, a menudo bajo condiciones ambientales duras y cargas sísmicas. El fracaso de una presa de cola puede liberar millones de metros cúbicos de esloria tóxica, causando daños y pérdida de vida.

Las evaluaciones de la estabilidad geotécnica tradicional dependen de métodos de equilibrio y de hipótesis simplificadas de presión poro. Mientras estos enfoques son útiles, a menudo no captan el comportamiento complejo y dependiente del tiempo de la lotería dentro de la presa. Las mezclas son multifase, fluidos no neontonianos que evolucionan con el tiempo a medida que las partículas se resuelven, consolidan y descifran.

Comprender el flujo de la lotería en los muelles

Las manchas son suspensión de partículas sólidas finas (típicamente se realizaron 200 μm) en agua, a menudo con una alta concentración sólida (40%-60% por masa). El comportamiento de flujo de estas suspensiones se rige por la distribución del tamaño de partículas, mineralogía y química del agua. A altas concentraciones, las manchas presentan características no-Newtonian: pueden ser titrópicas (vicosidad dependientes del tiempo)

Modelos retológicos para las rayas

Para modelar el flujo de lodo con precisión en ANSYS Fluent, los ingenieros deben definir un modelo reológico adecuado.

  • ]Bingham plastic] – una relación lineal de estriado vs. de escaneo después de que se exceda el estrés de rendimiento. Adecuado para las esclusas de alto grado que se comportan como un fluido viscoso.
  • Herschel‐Bulkley – un modelo más generalizado que representa el estallido o la vacilación después del rendimiento. Ampliamente utilizado para la cola porque se ajusta a los datos experimentales bien sobre una amplia gama de tasas de esquila.
  • Power law] – útil cuando el estrés de rendimiento es insignificante, pero a menudo inadecuada para las manchas establecidas o altamente concentradas.
  • Casson] – históricamente utilizado para las manchas de barro y minería, aunque menos común hoy.

La caracterización experimental (por ejemplo, mediante un reómetro rotacional) es esencial para obtener los parámetros de estos modelos. Sin datos reológicos precisos, las simulaciones CFD pueden producir resultados engañosos.

Asentamiento y consolidación

Con el tiempo, las partículas se asientan bajo gravedad, formando una cama densa y compactada en el suelo de la presa. Esta sedimentación cambia la reología de la rejilla de la rejilla local y afecta la evolución de la presión poro. En muchas presas de la cola, la lotería se deposita en capas finas (encaje), permitiendo que el agua se decanta y se consolide fases de grano.

El papel de la CDF en el análisis de la estabilidad de las presas

El análisis de estabilidad de una presa de colas implica evaluar el equilibrio entre las fuerzas de conducción (gravidad, presión de visor) y las fuerzas de resistencia (fuerza de ciervos y las colas mismas). La presión de agua poro es un factor crítico porque reduce el estrés efectivo, disminuyendo la fuerza de ciervo. CFD permite a los ingenieros simular la evolución de la presión poro a lo largo de la presa como cálculos esquilibes.

Las capacidades de los medios multifase y porosos de ANSYS Fluent permiten tratar los acoplamientos como una mezcla fluida y sólida que se mueve a través de un medio poroso deformable (las colas consolidadas y la estructura de presa). Coupling CFD con análisis geotécnico finito-elemento (por ejemplo, mediante un intercambio de datos de una sola dirección o de dos vías) proporciona aún más información, pero los campos de cálculo de velocidad de flujo ya pueden

Técnicas de modelado en ANSYS Fluent

El modelo físico adecuado en Fluent es la piedra angular de una simulación exitosa. Las siguientes subsecciones detallan los enfoques más relevantes para el flujo de rejilla de rejilla de rejilla.

Modelos de flujo multifase

Las barrigas son inherentemente multifase. ANSYS Fluent ofrece varias opciones:

  • ]Volumen del modelo Fluid (VOF) – mejor al rastrear una interfaz distinta entre el aire y la leña (o el agua), por ejemplo, durante el enfriamiento o en el estanque de decantación. VOF captura dinámicas de superficie libre pero no modela fases de partículas distintas; la lotería se trata como un solo fluido con una reología definida por el usuario.
  • Modelo eulerian‐Eulerian – trata cada fase (agua, sólidos) como interpenetrating continua con su propia fracción de velocidad, temperatura y volumen. Este es el enfoque más riguroso para las manchas densas donde las interacciones de partículas son significativas. Requiere modelos de cierre para la arrastre, elevación y presión granular (por ejemplo, la teoría granular).
  • Modelo de la micro] – un modelo eulerio simplificado que resuelve una ecuación de impulso única para la mezcla y utiliza relaciones algebraicas para las velocidades de deslizamiento entre fases. Apto cuando se ajustan las velocidades se caracterizan bien y la fracción de sólidos es moderada.
  • нерентениенитение Modelo de Fase (DPM) seleccionado/fuerteng confianza – trata las partículas como elementos lagrangos individuales en un fluido continuo. Apropiado para la carga de sólidos bajos (volumen de 10%) pero impráctico para colas densas debido al costo computacional.

Para las típicas presas de colas de alta resistencia, se recomienda el enfoque eulerian‐eulerian o el modelo de mezcla (con cierres apropiados de arrastrar y ajustar).El modelo VOF se puede utilizar junto con uno de estos modelos para el seguimiento de la superficie libre.

Modelado de Rheology no Newtonian

ANSYS Fluent permite a los usuarios definir la viscosidad no-Newtoniana a través de funciones definidas por el usuario (UDFs) o seleccionando uno de los modelos incorporados (Bingham, Herschel-Bulkley, etc.) para una simulación de colas:

  1. Elija el modelo reológico adecuado basado en datos experimentales. Almacene los parámetros (estres de rendimiento, índice de consistencia, índice de flujo) en la base de datos de materiales.
  2. Si la fracción de sólidos varía espacialmente, utilice un UDF para hacer que los parámetros reológicos dependan de la fracción de volumen local de sólidos. Esto es crítico porque la lotería cerca del punto de descarga es menos concentrada que la cama fijada.
  3. Permitir zonas de malla dinámica o de medios porosos para regiones que se han consolidado en un pastel de tipo sólido.

Un escollo común es el uso de una aproximación de viscosidad Newtoniana, que falla completamente el comportamiento de la tensión de rendimiento. Un pequeño estrés de rendimiento (por ejemplo, 5 Pa) puede afectar dramáticamente la disipación pore-presión y la propagación frontal del flujo.

Medios de comunicación y consolidación porosos

Una vez que se resuelven los ademas, se comportan como un medio poroso con una disminución de la permeabilidad a medida que se obtiene la consolidación. En Fluent, se pueden definir zonas porosas con permeabilidad direccional y una función definida por el usuario para actualizar la permeabilidad basada en la fracción sólida local. Esto permite simular el proceso de consolidación: inicialmente baja fracción sólida (alta permeabilidad) transiciones gradualmente a un pastel de baja permeabilidad.

Pasos en el proceso de modelado

Un flujo de trabajo robusto garantiza que la simulación sea precisa y computacionalmente factible.

Geometría y generación de malla

La geometría de una presa de colas incluye el terraplén de presa, el suelo de la cuenca, el estanque de decantación, y el gasoducto de descarga de lodo o el spigot. En muchos casos, la presa crece con el tiempo; un modelo 3D puede incorporar la construcción escenificada utilizando múltiples regiones o movimiento de malla dinámica.

  • Superficie libre (estéticamente refinada cerca de la interfaz de slurry-air).
  • Zona de descarga de lodos (gradientes de alta velocidad).
  • Interior y fundación de presas (donde se computan las presiones de poro).

Las mallas híbridas (tetraedral en regiones complejas, hexahedral en secciones rectas) son comunes. Un estudio de la rejilla-independencia debe realizarse por lo menos dos mallas (en bruto/fina) para asegurar que las presiones de poro computadas y los frentes de flujo no sean sensibles a la malla.

Condiciones de la frontera y propiedades materiales

Las principales condiciones de los límites son:

  • Inlet] – tasa de flujo de lodo, fracción de volumen sólido, parámetros reológicos (o un UDF para variar).
  • Expeto] – decantar superficie de estanque o drenaje de dedo. Típicamente se establece como salida de presión con una cabeza hidrostática fija.
  • Aguas – suelo de la cara y la cuenca de la presa: no es adecuado para las fases del fluido, con el estrés de la grieta computado de la reología de la llaga.
  • Simetría ] – utilizada si la geometría de la presa es simétrica.

Para la fase de sólidos en un modelo eulerio, especificar la densidad, el diámetro de las partículas y los modelos de viscosidad granular. La fase de fluido (agua) utiliza propiedades estándar. Recuerde especificar un tamaño apropiado del tiempo: los esquemas explícitos pueden requerir pasos de tiempo muy pequeños (por ejemplo, 0.001 s) para mantener la estabilidad en regiones de alta resistencia.

Ajustes y simulación de Solver

Para simulaciones transitorias de la deposición y consolidación de la lotería:

  • Utilice el solucionador basado en presión con el algoritmo PISO o SIMPLE para el acoplamiento de velocidad por presión.
  • Permite la ecuación de fracción de volumen para el modelo eulerio, o establece el esquema explícito de VOF si rastrea la superficie libre.
  • Seleccione una segunda orden de discretización espacial en viento para el impulso y la fracción del volumen para reducir la difusión numérica.
  • Establecer los factores de sub-relaxación cuidadosamente - flujos de resistencia- síield puede ser sensible a factores de relajación por debajo de 0.3.
  • Inicializar el dominio con un pequeño volumen de aire sobre la cama de colas iniciales. Usar regiones de parche para definir la altura inicial de la lotería.

Ejecute la simulación para tiempo físico suficiente para capturar el ciclo de deposición completo (horas a días de tiempo real). Compruebe la convergencia en cada paso del tiempo: los residuos deben caer por al menos tres órdenes de magnitud.

Procesos y validación post-

Después de la simulación, las herramientas de procesamiento postprocesamiento de ANSYS Fluent pueden visualizar:

  • Contornos de presión poro y estrés efectivo.
  • Volumen de sólidos para mostrar el perfil de playa y el frente de consolidación.
  • Parcelas vectoriales de velocidad de lodo, indicando zonas de flujo preferencial que pueden conducir a tuberías.
  • Ubicación de la superficie libre para comparación con las fotos aéreas o imágenes de drones.

Es esencial la validación. Compare las presiones poro computadas a varias profundidades con lecturas de piezometros de la presa real. Si la simulación coincide con los datos de campo dentro de la tolerancia aceptable (por ejemplo, ±15%), el modelo puede ser utilizado para predecir comportamiento futuro o para probar modificaciones de diseño.

Consideraciones prácticas para el análisis de la estabilidad

El objetivo final de la simulación de CFD es proporcionar datos para un análisis de estabilidad cuantitativa. Los siguientes productos son particularmente útiles:

  • Campo de presión poro] – directamente ingerido en el software de límite-equilibrio (por ejemplo, SLIDE, GeoStudio) a través de un mapa de presión. Las altas presiones poro en la masa de colas reducen la fuerza efectiva y aumentan el riesgo de licuación estática.
  • Distribución del estrés de la manta – El fluido puede calcular el estrés de la manguera ejercido por la lotería en el suelo de la presa y las paredes. Comparando estos elementos con la fuerza de la base de los materiales indica dónde es probable que el fracaso.
  • Ubicación de la superficie libre – la altura de la superficie farmacéutica dentro de la presa es una entrada clásica para los gráficos de estabilidad. Una superficie excesivamente alta farmacéutica sugiere drenaje insuficiente.
  • Senderos de flujo de lodo – vectores de velocidad pueden revelar la erosión interna (piping) si se encuentran zonas de alta velocidad cerca del núcleo de presa o la fundación.

Los ingenieros a menudo combinan los resultados de la CFD con un análisis de estrés finito de elemento (por ejemplo, usando Ansys Mechanical) para calcular el factor de seguridad contra el fallo de deslizamiento y rotación. Este enfoque combinado proporciona una imagen más completa que cualquiera de los métodos solos.

Casos de estudios y aplicaciones

Varias operaciones mineras han utilizado ANSYS Fluent para analizar su gestión de la cola. Por ejemplo:

  • Una mina de cobre en Chile modeló la deposición secuencial de las colas en una presa de estanqueidad para optimizar la ubicación del empuje y reducir el riesgo de sobre-topping. La simulación mostró que mover el punto de descarga 20 m río arriba redujo la altura de superficie férática en 3 m, mejorando el factor de seguridad de 1.2 a 1.5.
  • Una mina de oro en África Occidental utilizó un modelo eulerian‐eulerian con la retología Herschel-Bulkley para evaluar la tasa de disipación pore‐pressure en una presa de colas de paredes gruesas. Los resultados llevaron a la instalación de capas de drenaje horizontal adicional, que acortaron el tiempo de consolidación de 18 meses a 10 meses.

Aunque los detalles son propietarios, estos ejemplos demuestran el valor práctico de la CFD en la ingeniería de relevos.

Limitaciones y desafíos

No es perfecta la simulación. Las limitaciones clave del enfoque descrito incluyen:

  • Costo computacional – simulaciones Eulerian‐Eulerian de alta resolución con modelos de coupled DEM pueden tomar días en una estación de trabajo multi-core. Escalar a modelos transientes 3D a gran escala puede requerir computación de racimo.
  • Incertidumbre retológica – las propiedades de la cola cambian con contenido mineral, forma de partículas y química de agua. La simulación es tan buena como los datos de la reología de entrada, que pueden no ser representativos de todas las fases de deposición.
  • leyes constitutivas geotécnicas simplificadas] – Fluent no modela nativamente el comportamiento de estrés-estrés de la estructura de presas sólidas. Coupling CFD con FEA geotécnica es todavía una tarea especializada que requiere una cuidadosa asignación de datos.
  • Datos de validación] – muchas minas carecen de instrumentación integral (piezometers, inclinometers) para comparar. Sin validación, la simulación sigue siendo una herramienta predictiva con precisión desconocida.

A pesar de estos desafíos, el uso de CFD para la estabilidad de las presas se ha incrementado rápidamente a medida que aumenta la potencia computacional y a medida que la presión regulatoria exige más evaluaciones de riesgo cuantitativos.

Conclusión

Modelo de la corriente de las rayas en las presas de las colas mineras utilizando ANSYS Fluent proporciona una base poderosa para el análisis de la estabilidad. Al capturar la naturaleza multifase de la lotería, su reología no neonquista, y la evolución de la presión poro durante el tiempo, los ingenieros pueden identificar mecanismos de falla que los métodos más simples faltan.

Para más información sobre la física subyacente y las mejores prácticas, consulte ANSYS Fluent documentation, the International Council on Mining and Metals (ICMM) tailings management guidelines, and the Geoengineer.org tailings dam resources[LTolic]