Reconstrucción numérica del aluvión de San Carlos de Apoquindo (Las Condes, 8-oct-2026)

Documento vivo de metodología: describe datos, modelo, calibración, revisión externa e incertidumbre de forma reproducible.

Estado: en curso (ver «Registro de cambios» al final).

1. Objetivo

Estimar el recorrido del flujo de barro desde la salida de la quebrada hasta la zona urbana afectada, el caudal y volumen de mezcla

compatibles con el alcance observado, y los indicadores hidrodinámicos (profundidad, velocidad, presión dinámica, tiempo de llegada)

sobre calles y edificios, con un nivel de confianza explícito.

2. Observaciones disponibles (ver fuentes/FUENTES.md)

Solo hay observaciones cualitativas de prensa y del municipio (calles cerradas, polígono de afectación, origen informado,

dirección surponiente). No hay aforos, marcas de agua, ni cartografía oficial del área inundada a la fecha del estudio.

Esto limita la calibración a la extensión y obliga a reportar los parámetros como rangos compatibles, no como valores únicos.

3. Datos geográficos y preparación (extraccion/, prep.py)

| Capa | Fuente | Resolución | Uso |

|---|---|---|---|

| Terreno sin edificios (DTM) | FABDEM V1-2 | 1" (~30 m) | elevación base, remuestreo bilineal a 4 m |

| Superficie (DSM) | Copernicus GLO-30 | 1" | altura de bloques = DSM − DTM |

| Calles, cauces, edificios | OpenStreetMap (Geofabrik, ver extraccion/osm_pbf_version.txt) | vectorial | canales urbanos, cauce entubado, plantas de bloques |

• Grilla de trabajo UTM 19S (EPSG:32719), 4 m.
• Bloques: altura = media de (DSM − DTM) en la planta si la planta ≥ 900 m² y la diferencia > 2,5 m; si no, building:levels × 3 m;

• Calles rebajadas 0,15 m (cordón) con ancho por categoría (18/14/11/8/6 m).
• El cauce de la quebrada continúa entubado bajo la urbanización (OSM tunnel=culvert); hipótesis de trabajo H1: el ducto se obstruyó

4. Modelo (sim.py)

Aguas someras 2D (ANUGA 4.0.1, volúmenes finitos, algoritmo DE1) con mezcla homogénea y resistencia de flujo de barro de

O'Brien & Julien (1988) por paso fraccionado semi-implícito. Ecuaciones y coeficientes en el encabezado de sim.py.

Malla no estructurada: 60–120 m² por triángulo en el corredor del abanico a Padre Hurtado, 600 m² fuera.

4b. Verificación numérica (verif_parada.py, verif_conv2.py)

• Parada (Bingham en plano inclinado, S = 0,05, τ_y = 100 Pa, ρ = 1.700 kg/m³): espesor teórico de reposo h_c = τ_y/(ρ g S) = 0,1199 m;

• Reinicio: al agregar mezcla aguas arriba (t = 900 s) el depósito se reactiva y avanza 70 m; volumen conservado (746,67 + 120 m³).
• Convergencia temporal (DE1): alcance del frente 455 / 387 / 352 / 331 / 324 m para CFL 1 / 0,5 / 0,25 / 0,1 / 0,05

Consecuencia de método: la calibración se hace a CFL 1 (barrido) y las corridas finales se repiten a CFL 0,25 y 0,1 con malla

refinada para estimar el error numérico de alcance (índice de convergencia de malla) y reportarlo como banda.

• DE0 se descarta: no responde al CFL y deja depósitos más delgados que la teoría (0,098 vs 0,120 m).

4c. Volumen previo independiente (cuenca.py, volumen_previo.py)

• Cuenca aportante a H1 (FABDEM, D8): A = 2,18 km², relieve R = 890 m (1.869 → 979 m s.n.m.).
• Lluvia: DMC EMA 330019 Tobalaba — 8-oct total 34,8 mm; 15–21 h local 13,8 mm (el evento ~17:00); 49,2 mm el 6–7 oct (suelo saturado);

• Balance lluvia–escorrentía (Monte Carlo, P 12–35 mm, C 0,4–0,9, C_v 0,30–0,50): V_mezcla mediana 53 mil m³ (p5–p95: 26–96 mil).
• Modelo USGS para cuencas quemadas (Gartner et al. 2014; incendio dic-2025/ene-2026 ~1.000 ha): mediana 48 mil m³ (28–77 mil).
• Rango previo para calibrar: 30–78 mil m³. Caudal pico empírico de flujo de barro (Mizuyama et al. 1992) ≈ 100 m³/s para V = 53 mil m³.

5. Calibración (calibrar.py, barrido.sh)

Barrido de caudal pico Qp, volumen V, concentración Cv y tipo de suelo; métrica F = media armónica de sensibilidad sobre calles

observadas y precisión respecto del polígono municipal, penalizada si el frente no alcanza Padre Hurtado. Cada corrida queda en

5b. Hipótesis de origen

• H1: obstrucción del ducto bajo el abanico (Portillo × La Plaza / General Blanche × San Carlos de Apoquindo) — Chilevisión, La Tercera.
• H2: desborde/colapso del canal en Francisco Bulnes Correa (cruce del cauce entubado con esa calle, inicio de calle Quebrada Honda) —

6. Revisión externa

Codex (modelo gpt-6-astra), solo física y matemática, nivel publicación: dos rondas sobre hidrodinámica y reología (A1, A2) y dos rondas

sobre indicadores de salida (B1, B2). Enunciados e informes completos en revision/.

7. Incertidumbre y confianza

(por completar con los resultados de la calibración y la revisión)

8. Registro probatorio

captura en fuentes/capturas/ y hash sha256 (fuentes/registrar.py). Permite contrastar después la predicción del modelo con lo

observado y reproducir exactamente qué información se tenía en cada momento.

Registro de cambios

• 2026-10-09: extracción de datos, terreno, modelo ANUGA + reología, métricas de calibración; ronda A1 de revisión lanzada.
• 2026-10-09: ronda A1 recibida (dictamen: correcciones mayores). Aplicado: coeficientes FLO-2D corregidos (Glenwood 1, Aspen 1 y 2;

en el eje del cauce, hidrograma validado y t_fin ≥ T_B, balance de volumen, ½ρU², llegada al final del paso, desfase de 1 m de las

máscaras, Monte Carlo de terreno disponible (--sigma, --ell). Verificación de parada/reinicio y convergencia temporal. Volumen

previo independiente (lluvia DMC + USGS). Registro probatorio de consultas. Hipótesis H2 (Fco. Bulnes Correa) agregada.

Revisión externa (Codex gpt-6-astra)