Enfoque guia: Varianza y desviacion estandar
¿La planta esta estable antes de perder recuperacion?
Ver concepto general MathPlay →Varianza y desviación estándar
Caso: sobre-espumación en flotación de cobre — detectar inestabilidad antes de perder recuperación
Idea central
Objetivo de esta lección: no construir aún un modelo predictivo sofisticado. Entender si la planta ya perdió estabilidad antes de que caiga la recuperación.
Problema (caso hilo conductor): en una planta concentradora de cobre ocurre sobre-espumación en celdas de flotación. Eso desestabiliza nivel de pulpa, espuma, arrastre, ley de concentrado y recuperación.
Pregunta operativa: ¿Podemos detectar condiciones que anticipan la sobre-espumación y estimar su impacto sobre la recuperación?
Respuesta de este paso (solo con σ): la media de recuperación puede verse “aceptable”, pero si la recuperación salta mucho hora a hora (σ alta), la operación ya está nerviosa — eso suele aparecer antes de un evento fuerte de espuma.
1. Matemática simple
Promedio (media) — μ
Recuperaciones medidas en un turno:
88 % 91 % 87 % 90 % 89 %
Media = suma ÷ cantidad = 89 %
En planta: “¿Cuánto recuperamos en promedio?”
Varianza — σ²
Mide qué tan dispersos están los valores respecto al promedio.
- Restas cada valor menos la media.
- Elevas al cuadrado (para que los negativos no se cancelen).
- Promedias esos cuadrados (en Excel: divides entre n − 1).
Ejemplo con media 89:
| Hora | x | x − μ | (x − μ)² |
|---|---|---|---|
| 1 | 88 | −1 | 1 |
| 2 | 91 | +2 | 4 |
| 3 | 87 | −2 | 4 |
| 4 | 90 | +1 | 1 |
| 5 | 89 | 0 | 0 |
Varianza muestral (como Excel): 10 ÷ 4 = 2,5 (unidades poco intuitivas: “% al cuadrado”)
Desviación estándar — σ
σ = √varianza
√2,5 ≈ 1,58 puntos porcentuales
En planta: “¿Cuánto típicamente se aleja la recuperación del promedio?”
Interpretación rápida (recuperación)
| σ recuperación | Lectura operativa |
|---|---|
| < 1,5 % | Turno bastante estable |
| 1,5 – 3 % | Variabilidad moderada; vigilar |
| > 3 % | Operación nerviosa; revisar aire, nivel, espumante |
Umbrales orientativos para el taller. En planta real se calibran con histórico propio.
2. Indicadores del caso (tu presentación)
Variables del PPT y del archivo datos_flotacion_demo.csv donde tiene sentido calcular σ:
| Prioridad | Indicador | Columna CSV | ¿Por qué σ? |
|---|---|---|---|
| Hoy 1 | Recuperación Cu | recuperacion_cu_pct | KPI central: inestabilidad metalúrgica |
| Hoy 2 | Nivel de celda | nivel_celda_pct | Espuma inestable → celda “nerviosa” |
| Hoy 3 | Flujo de aire | flujo_aire | Régimen de burbujas / espuma |
| Tarea | Ley de relave | ley_relave_cu_pct | Señal temprana: más Cu al relave |
| Tarea | Ley de concentrado | ley_concentrado_cu_pct | Balance ley vs recuperación |
| Tarea | Densidad de pulpa | densidad_pulpa | Estabilidad pulpa–espuma |
| Tarea | P80 | p80_um | Variabilidad de alimentación |
| Tarea | Ley de cabeza | ley_cabeza_cu_pct | Perturbación de feed |
En clase trabajamos solo prioridad 1, 2 y 3. No uses primero σ en observaciones del operador (texto) ni en tonelaje (más planificación que estabilidad de celda).
3. Qué puedes hacer con σ en este caso
- A) Estabilidad por turno — σ de recuperación en turnos A, B y C. Pregunta: ¿misma media pero distinta estabilidad?
- B) Antes vs evento — σ cuando
sobre_espumacion = 0vs= 1(tarea: ventana 3–6 h antes). - C) Dos promedios iguales — dos periodos con media ~85 % pero σ distinta: el operador que solo mira el promedio pierde la señal temprana.
- D) Semáforo de estabilidad — primer “soft sensor” sin IA (tabla abajo).
- E) Puente a la pregunta del PPT — con σ aún no predices el evento, pero detectas: la planta ya no está en régimen estable.
Semáforo (solo recuperación)
| σ recuperación | Estado |
|---|---|
| ≤ 2 | Verde — estable |
| 2 – 3 | Amarillo — vigilar |
| > 3 | Rojo — operación nerviosa |
Datos de referencia (archivo demo)
| Variable | Media | σ |
|---|---|---|
| Recuperación Cu (%) | 85,3 | 3,41 |
| Nivel celda (%) | 69,9 | 8,41 |
| Flujo de aire | 122,3 | 12,42 |
| Turno | Media rec. | σ rec. |
|---|---|---|
| A | 84,6 | 3,88 |
| B | 85,0 | 3,58 |
| C | 86,1 | 2,51 |
Turno C: mejor media y menor σ → más estable. Compara con filas donde obs_operador dice “espuma inestable”.
4. Excel — paso a paso
Abre o descargue CSV demo. Verifica en qué columna está recuperacion_cu_pct (en muchas importaciones es la columna G).
4.1 Toda la planta
| Concepto | Fórmula (Excel en español) |
|---|---|
| Media | =PROMEDIO(G2:G121) |
| Varianza | =VAR.MUESTRA(G2:G121) |
| Desviación estándar | =DESVEST(G2:G121) |
4.2 Por turno
| Turno | Media | σ |
|---|---|---|
| A | =PROMEDIO.SI(C:C,"A",G:G) | Tabla dinámica o filtrar filas |
| B | =PROMEDIO.SI(C:C,"B",G:G) | idem |
| C | =PROMEDIO.SI(C:C,"C",G:G) | idem |
DESVEST.SI no existe en todas las versiones. Usa Tabla dinámica: filas = turno, valores = Promedio y Desv.estánd de recuperación.
4.3 Nivel y aire
Repite con nivel_celda_pct y flujo_aire (revisa letra de columna al importar).
Pregunta en clase: ¿Cuál σ es más alta? ¿Coincide con horas “espuma inestable”?
4.4 Gráfico
- Eje X:
fecha_hora - Eje Y:
recuperacion_cu_pct - Líneas horizontales: media y media ± σ
El gráfico muestra los picos; σ los resume en un solo número.
5. Python (opcional — validar Excel)
Abre o descargue CSV demo.
Si el resultado coincide con Excel (± redondeo), tu pipeline de datos es confiable — igual que revisar SCADA vs laboratorio.
import pandas as pd
df = pd.read_csv("datos_flotacion_demo.csv")
r = df["recuperacion_cu_pct"]
print("Media:", round(r.mean(), 2))
print("Varianza:", round(r.var(), 2)) # ddof=1, como Excel
print("Desv. estándar:", round(r.std(), 2))
print(df.groupby("turno")["recuperacion_cu_pct"].agg(["mean", "std"]))
6. Proyección a Machine Learning
Hoy no entrenas un modelo. Solo entiende para qué servirá σ después:
| Idea | Para el alumno |
|---|---|
| Escala de variables | Aire ~120, pH ~10, recuperación ~85 → en ML se normaliza; σ muestra qué variable “mueve mucho”. |
| Variable derivada | σ móvil en 6 h → feature de inestabilidad reciente. |
| Calidad del dato | σ absurda en sensor → el modelo aprende ruido. |
| Evaluación | Si σ del turno es alta, el error del modelo de recuperación también suele subir. |
7. Investigación Sugerida
Tarea
Calcula σ de ley de relave y entrega medio página: ¿sube cuando la recuperación es inestable?
8. Siguiente paso
Paso 2: covarianza y correlación — ¿qué variables se mueven junto con la inestabilidad?
No avances hasta interpretar σ sin mirar solo el promedio del reporte.
Regístrate al boletín MathPlay
Datos, simulación, visualización, IA y tecnología aplicada a problemas reales de industria y minería.
Sin calendario fijo. Sin spam.
Listo. Te escribo cuando haya algo nuevo.
