Simulación probabilística de ruptura de presas y de diques: modelado dinámico de erosión geotécnica y enrutamiento de ondas hidrodinámicas

Las presas, los embalses y los sistemas de diques de protección contra inundaciones representan algunos de los activos de infraestructura civil más críticos de la sociedad moderna. Aseguran el suministro de agua para el riego agrícola, generan energía hidroeléctrica renovable y protegen de las inundaciones a millones de residentes que viven en llanuras aluviales históricas. Sin embargo, la falla catastrófica de una presa de tierra o un terraplén de dique libera grandes volúmenes de energía potencial en minutos, generando olas de inundación violentas y catastróficas capaces de arrasar con comunidades enteras, puentes y corredores industriales río abajo.

Históricamente, las evaluaciones del riesgo de rotura de presas se basaban en fórmulas de rotura empíricas y estáticas (como las ecuaciones clásicas de Froehlich, MacDonald-Langridge-Monopolis o Von Thun-Gillette) combinadas con rutas hidráulicas unidimensionales (1D). En estos marcos heredados:

• El ingeniero asumió una geometría de brecha fija (forma de brecha trapezoidal, ancho del fondo y pendientes laterales) y un tiempo de desarrollo de brecha fijo (por ejemplo, 1,5 horas) basado en gráficos de dispersión de regresión histórica.
• El flujo de salida de la brecha se dirigió río abajo a través de secciones transversales 1D que ignoraron por completo la extensión del valle de la llanura aluvial 2D, la densidad aparente del transporte de sedimentos y los reflejos de las olas en las paredes del cañón.

Este enfoque empírico estático presenta peligrosas limitaciones de ingeniería. Los terraplenes de tierra no fallan siguiendo esquemas geométricos arbitrarios; fallan a través de complejos, no linealesmecánica de erosión geotécnica: erosión progresiva hacia atrás por desbordamiento, erosión por corte de conductos de tuberías internas y colapsos del círculo de deslizamiento de estabilidad de pendiente rotacional. Además, los supuestos estáticos ignoran la incertidumbre de los parámetros geotécnicos: cohesión de compactación del suelo, distribución del tamaño de grano (D_ {50}), esfuerzo cortante crítico (tau_c) y fluctuaciones del nivel del agua del embalse.

A finales de 2026, las normas regulatorias de seguridad (exigidas por la Comisión Federal Reguladora de Energía, FERC, y el Cuerpo de Ingenieros del Ejército de EE. UU., USACE) requierenSimulaciones probabilísticas de rotura de presas y de diques basadas físicamente. Al combinar la física de erosión de brechas geotécnicas transitorias directamente con solucionadores de ecuaciones de aguas poco profundas bidimensionales (2D) acelerados por GPU y la cuantificación de incertidumbre del conjunto Monte Carlo, los ingenieros de seguridad de presas pueden mapear los peligros de inundación aguas abajo con una precisión matemática sin concesiones.

Este tratado técnico integral explora la física de la erosión geotécnica, la hidrodinámica de las olas cargadas de sedimentos no newtonianos, el muestreo de conjuntos de Monte Carlo y los protocolos de mapeo de peligros de evacuación de emergencia necesarios para ejecutar evaluaciones de rotura de presas de última generación.

1. Mecánica de brechas geotécnicas: desbordamiento versus erosión interna de tuberías

Las presas de terraplén de tierra fallan predominantemente a través de uno de dos mecanismos físicos: desbordamiento hidráulico o tuberías geotécnicas internas.

La ecuación de erosión por esfuerzo cortante excesivo

En los modelos de brecha basados físicamente (como USDA-ARS WinDAM o DL-Breach), el desprendimiento del suelo a lo largo del límite de la brecha se rige por laFormulación de esfuerzo cortante excesivo:

ε_r = k_d (τ_e - τ_c)^a

donde:

• épsilon_rrepresenta la tasa de erosión por desprendimiento del suelo (m/s o cm/hr).
• k_des el coeficiente de erosionabilidad del suelo (cm^3/N·s), determinado mediante pruebas de erosión por chorro (JET) in situ o pruebas de dispersión por orificios sumergidos.
• tau_ees el esfuerzo cortante hidrodinámico efectivo ejercido por el agua que fluye en el límite de la brecha:

τ_e = γ_w R_h S_f

dondegamma_wes el peso unitario del agua,R_hes el radio hidráulico del canal de ruptura, yS_fes la pendiente de la energía.

• tau_ces el esfuerzo cortante crítico del suelo compactado del terraplén (Pa), por debajo del cual se produce cero desprendimiento.
• $a$ es un exponente empírico, normalmente tomado como $a = 1,0$.

Cuando el agua que fluye acelera a través de una muesca de brecha en desarrollo,tau_epuede exceder los 500 Pa, superando ampliamente la resistencia al corte crítica del suelo (tau_c ≈ 5 - 20Pa)—impulsando tasas de erosión volumétrica exponencial.

2. Cálculo del hidrograma de flujo de salida dinámico: hidráulica de orificio y vertedero de cresta ancha

A medida que la muesca de la brecha se ensancha y profundiza, la descarga instantánea del flujo de salidaQ_{incumplimiento}(t)cambia transitoriamente entreflujo de orificio sumergido(durante las primeras tuberías internas) yflujo de vertedero de cresta ancha y superficie libre(una vez que el techo de la brecha se derrumbe).

flowchart TD
    Pool["Reservoir Pool Level: H_res(t)"] --> Geom["Instantaneous Breach Notch Geometry<br/>• B_w(t), Z_s, H_b(t)"]
    Geom -->|Internal Piping Phase| Pipe["Pressurized Orifice Flow<br/>• Q = C_d · A_pipe(t) · √(2g · ΔH)"]
    Geom -->|Open Breach Crest Phase| Weir["Broad-Crested Weir Flow<br/>• Q = C_w · B_w(t) · (H_res - Z_b)^1.5 + C_s · Z_s · (H_res - Z_b)^2.5"]
    Pipe --> Wave["Downstream 2D Hydrodynamic Flood Wave Routing"]
    Weir --> Wave

La formulación combinada de vertedero trapezoidal de cresta ancha

Una vez que se establece una zanja de brecha abierta, el flujo de salida se calcula mediante la ecuación combinada del vertedero de cresta ancha rectangular y triangular:

Q_{incumplimiento}(t) = C_{w1} B_w(t) ≤ft( H_{res}(t) - Z_b(t) )^{3/2} + C_{w2} Z_s ≤ft( H_{res}(t) - Z_b(t) )^{5/2}

donde:

• B_w(t)es el ancho instantáneo del fondo de la muesca de la brecha en expansión (m).
• Z_b(t)es la elevación instantánea del fondo de la muesca de la brecha, erosionándose hacia abajo hasta alcanzar un lecho de roca o estratos de cimentación no erosionables.
• Z_ses la pendiente lateral de la brecha trapezoidal (Z_s = (Delta x) / (Delta z), normalmente entre $ 0,5 $ y $ 1,0 $ para arcillas cohesivas y hasta $ 1,5 $ para arenas no cohesivas).
• C_{w1}, C_{w2}son coeficientes de descarga del vertedero calibrados para los efectos de la sumersión del agua de descarga.

Simultáneamente, la elevación de la piscina del embalseH_res}(t)cae dinámicamente según la ecuación de conservación etapa-almacenamiento del yacimiento:

frac{dH_{res}}{dt} = frac{Q_{entrada}(t) - Q_{brecha}(t) - Q_{aliviadero}(t)}{A_{res}(H_{res})}

dondeA_res}(H_res})es el área de superficie de la piscina del embalse en elevaciónH_res}.

3. Enrutamiento de ondas hidrodinámicas: ondas de rotura de presas no newtonianas en 2D

Una ola que rompe una presa no es agua limpia y clara; es una mezcla violenta y de alta densidad de agua, grava de terraplén pulverizada, sedimentos en suspensión y árboles arrancados de raíz.

En cañones montañosos escarpados o valles fluviales estrechos, la concentración volumétrica de sedimentos (C_v) puede exceder $20\% - 45\%$, transformando la ola de inundación en unasuspensión hiperconcentrada no newtoniana.

Hidrodinámica plástica de Bingham no newtoniana

>

Formulación de tensión de corte aparente:

τ = τ_y + µ_B · (∂u / ∂z)

>

Parámetros:

• τ_y: El límite elástico de la suspensión de sedimentos (Pa) debe superarse para iniciar el movimiento del fluido.

• µ_B: Viscosidad plástica dinámica de Bingham (Pa·s)

• Densidad aparente: ρ_lechada = ρ_w · (1 - C_v) + ρ_s · C_v

El solucionador hidrodinámico de captura de impactos 2D

Para encaminar esta onda violenta y de alta velocidad aguas abajo sin inestabilidad numérica ni difusión no física:

• Los modelos utilizan las ecuaciones de aguas poco profundas 2D conTérminos de origen de resistencia al flujo de plástico Bingham o Herschel-Bulkley (tau_ {bx}, tau_ {por}).
• Los solucionadores aplican esquemas de captura de shock TVD (Disminución de variación total) de segundo orden para resolverPerforaciones de olas de rotura de presasavanzando a velocidades superiores12 a 20 m/spor el valle río abajo.
• Las ecuaciones de momento resuelven las fuerzas centrífugas laterales a medida que la ola recorre curvas cerradas del río, prediciendo con precisión el avance de la onda de superelevación contra las paredes del cañón.

4. Marco probabilístico del conjunto Monte Carlo: cuantificación de la incertidumbre geotécnica

En la ingeniería forense de seguridad de presas, las propiedades geotécnicas exactas nunca se conocen con un 100% de certeza. La calidad de la compactación, el contenido de humedad, la cohesión y la homogeneidad interna de un terraplén de tierra varían a lo largo de su longitud lineal.

En lugar de ejecutar una única simulación determinista del "peor de los casos" que puede ser irrealmente conservadora o peligrosamente optimista, los ingenieros modernos ejecutanConjuntos probabilísticos de Montecarlo.

flowchart TD
  Params["Input Parameter Probability Distributions\n• Reservoir Head: Normal(μ = Full Pool, σ = 1.2m)\n• Soil Erodibility k_d: Log-Normal(μ = 10 cm³/N-s, σ = 0.4)\n• Critical Shear Stress τ_c: Weibull(λ = 8.5 Pa, k = 2.1)\n• Cohesion c': Normal(μ = 25 kPa, σ = 5 kPa)\n• Breach Width B_w: Beta distribution (Froehlich)"]
  LHS["Latin Hypercube Sampling (LHS)\nGenerates 500 Stratified Scenarios"]
  GPU["Automated GPU Solver Cluster\n500 Coupled 2D Simulations (4.2 GPU Hours)"]
  Hazard["Probabilistic Hazard Outputs\nP10, P50, P90 Inundation Maps & Arrival Time Bands"]

  Params --> LHS --> GPU --> Hazard

Métricas probabilísticas de inundaciones por inundaciones

Del conjunto de 500 corridas, el sistema extrae probabilidades de excedencia distribuidas espacialmente:

• Mapa de inundaciones P90: El límite de inundación conservador con un nivel de confianza de $90\%$, utilizado para el diseño de protección de infraestructura civil (plantas de energía nuclear, complejos hospitalarios, almacenamiento de productos químicos).
• Mapa de inundaciones P50: El escenario esperado mediano, utilizado para la zonificación municipal, la evaluación de tarifas de seguros y la ruta de evacuación primaria.
• Incertidumbre en el tiempo de llegada de la onda: Cuantifica la llegada más temprana posible del frente de onda de alta velocidad a centros de población específicos aguas abajo (T_{llegada, P10}), proporcionando a los gestores de emergencias márgenes de seguridad defendibles para las ventanas de autorización de evacuación.

5. Matriz de comparación integral de metodologías de rotura de presas

La siguiente tabla contrasta los métodos empíricos heredados con los marcos modernos de simulación de rotura de presas 2D probabilísticos basados físicamente:

6. Estudio de caso de ingeniería del mundo real: Auditoría de seguridad de presas de tierra de alto riesgo

Para validar este marco probabilístico en un activo de infraestructura civil crítico, FloodSim ejecutó una auditoría de ruptura de presa del Plan de Acción de Emergencia (EAP) para una presa de terraplén de tierra de 85 pies (26 metros) de altura que incauta un depósito de agua potable y control de inundaciones de 32,000 acres-pie (39,5 millones de metros cúbicos) situado a 6 millas río arriba de un valle municipal urbano de 95,000 residentes.

Descubrimientos críticos para la seguridad de la vida:

1. La peligrosa ilusión del "retraso": El modelo 1D heredado predijo que la ola de inundación inicial necesitaría 65 minutos para viajar 6 millas río abajo hasta el centro de la ciudad. Por el contrario, el solucionador de GPU 2D de alta resolución, que tiene en cuenta el estrechamiento del valle y la canalización del flujo supercrítico a lo largo del corredor del río principal, reveló que la perforación hidráulica de alta velocidad golpearía el núcleo municipal ensolo 38 minutos. Si los gestores de emergencias hubieran confiado en el modelo heredado, las órdenes de evacuación habrían llegado casi 30 minutos tarde.
2. Súper socavación del puente interestatal: La simulación 2D mapeó velocidades de flujo localizadas que exceden6,2 m/sen los pilares de un puente de una carretera interestatal, calculando las tensiones cortantes del fondo que socavarían los cimientos del muelle dentro de los 12 minutos posteriores a la llegada de la ola.
3. Recalibración del Plan de Acción de Emergencia: Armada con curvas de llegada probabilísticas exactas, la agencia municipal de gestión de emergencias reubicó cuatro sirenas de evacuación, estableció cierres automatizados de puertas de carreteras y estableció un protocolo certificado de autorización de evacuación de 22 minutos, garantizando la total seguridad de la vida pública.

7. Clasificación de peligros de evacuación aguas abajo: las curvas australianas de vulnerabilidad a peligros

La evaluación del riesgo para la seguridad de la vida pública requiere traducir la profundidad bruta de la inundación ($h$) y la velocidad ($v$) en categorías de peligro para la estabilidad humana. Una profundidad de agua poco profunda de 0,4 metros que fluye a 3,0 m/s es mucho más letal que un charco de agua estancada de 1,0 metro, ya que las fuerzas de arrastre hidrodinámicas arrastran a los peatones y automóviles adultos fuera de las carreteras.

Los planes de acción de emergencia modernos utilizan elCurvas australianas de vulnerabilidad a peligros en la gestión de emergencias (H1 a H6):

<!--ec:bloque {"tipo":"tabla","hasHeaderRow":true,"filas":[{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"CATEGORÍAS DE CLASIFICACIÓN DE RIESGOS (H1 - H6)"}],"isHeader":true}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"Intensidad hidrodinámica total: D \\times V = h \\times \\sqrt{u^2 + v^2} (m²/s)"}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"tipo":"span","text":""}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"[H1: Seguro para todos]: D \\times V \\le 0,4 \m²/s, h \\le 0,3m"}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"Capaz de caminar con seguridad; niños y ancianos seguros."}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"tipo":"span","text":""}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"[H2: Inseguro para pequeños]: D \\times V \\le 0,6 \m²/s, h \\le 0,5m"}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"No es seguro para niños pequeños y ancianos frágiles."}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"tipo":"span","text":""}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"[H3: Inseguro para vehículos]: D \\times V \\le 0,8 \m²/s, h \\le 0,6m"}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"Los sedanes y vehículos ligeros pierden tracción y flotan."}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"tipo":"span","text":""}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"[H4: Inseguro para todas las personas]: D \\times V \\le 1,0 \m²/s, h \\le 1,0m"}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"Adultos arrastrados; todos los vehículos flotando."}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"tipo":"span","text":""}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"[H5: Daño estructural]: D \\times V \\le 4.0 \m²/s, h \\le 4.0m"}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"Los edificios de mampostería no reforzada sufren daños."}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"tipo":"span","text":""}]}]},{"tipo":"tableRow","celdas":[{"tipo":"tableCell","content":[{"type":"span","text":"[H6: Colapso catastrófico]: D \\times V > 4,0 \m²/s"}]}]},{"tipo":"tableRow","celdas":[{"type":"tableCell","content":[{"_type":"span","text":"Falla estructural completa de edificios residenciales"}]}]}]} -->

Exportando capas SIG dinámicas deH_1a través deH_6Con clasificaciones en incrementos de 5 minutos aguas abajo de la presa rota, los administradores de emergencias pueden identificar qué caminos de evacuación siguen siendo transitables de manera segura y en qué momento exacto se convierten en trampas mortales letales.

8. Protocolo de cumplimiento normativo para ingenieros de seguridad de presas

Los ingenieros de seguridad de represas, consultores geotécnicos y directores de servicios municipales que realicen evaluaciones de riesgos de falla de represas en 2026 deben cumplir con esta lista de verificación de cumplimiento de cinco pasos:

1. Realizar pruebas de erosión por chorro in situ (JET): Nunca estime los coeficientes de erosionabilidad (k_d) únicamente a partir de tablas de clasificación de suelos de libros de texto. Realice pruebas estandarizadas ASTM D5852/D6665 para calibrar los parámetros reales de exceso de corte.
2. Modelo de Ampliación de Brecha Geotécnica Dinámica: Descartar los supuestos trapezoidales estáticos. Emplear algoritmos de erosión acoplados con base física que modelen la migración del corte de cabeza hacia atrás, el colapso rotacional de talud lateral y la caída del estanque del embalse.
3. Implemente solucionadores de captura de impactos 2D acelerados por GPU: Ruta de las ondas de salida de ruptura a través de mallas no estructuradas 2D derivadas de LiDAR de alta resolución, lo que garantiza que la topografía crítica del valle, los terraplenes de carreteras y las alcantarillas estén hidroacondicionados.
4. Cuantificar la incertidumbre de los parámetros geotécnicos: Ejecute un mínimo de 250 ejecuciones de conjuntos Monte Carlo a través de distribuciones estratificadas de parámetros de hipercubo latino para generar contornos de peligro defendibles P10, P50 y P90.
5. Clasifique las rutas de evacuación aguas abajo mediante curvas de profundidad-velocidad: Traducir las salidas hidrodinámicas al estándarH_1 - H_6clasificaciones de peligros humanos y vehiculares para brindar inteligencia procesable y que salve vidas a los directores de gestión de emergencias.

A través de la integración disciplinada de la mecánica geotécnica del suelo, la hidrodinámica de aguas poco profundas no newtoniana acelerada por GPU y la cuantificación probabilística de la incertidumbre de Monte Carlo, la simulación moderna de rotura de presas permite a los ingenieros civiles proteger a la sociedad contra los escenarios de falla más catastróficos imaginables.