Modelado de inundaciones compuestas acopladas: integración de la escorrentía de las cuencas fluviales, los diluvios de tormentas pluviales y las oleadas costeras de ADCIRC

Los estuarios costeros, los deltas y los centros urbanos marítimos representan las regiones ecológicamente más dinámicas y económicamente vitales de la Tierra. Sin embargo, también ocupan el punto de mira de un fenómeno hidrodinámico complejo y devastador:Inundación compuesta. Las inundaciones compuestas ocurren cuando múltiples mecanismos distintos de generación de inundaciones coinciden en el espacio y el tiempo o se suceden en rápida sucesión cronológica, lo que hace que los niveles de agua sean mucho más altos que la suma estadística de los peligros individuales aislados.

Históricamente, las prácticas de ingeniería civil y el mapeo regulatorio de peligros de inundaciones (como los Estudios de Seguros contra Inundaciones del Programa Nacional de Seguro contra Inundaciones de FEMA) evaluaban los peligros de inundaciones de forma completamente aislada:

• Hidrólogos evaluadosinundaciones fluviales(descarga fluvial que sobrepasa las orillas de los cauces principales) utilizando modelos de escorrentía de cuencas como HEC-HMS y motores de enrutamiento de ríos 1D/2D como HEC-RAS, asumiendo límites de marea oceánicos bajos y estáticos aguas abajo.
• Oceanógrafos costeros evaluadosinundaciones oceánicas(marejadas ciclónicas, configuración de olas de viento y mareas reales astronómicas) utilizando modelos de circulación oceánica costera como ADCIRC y SWAN, suponiendo límites de descarga de ríos de tierras altas estáticos o nulos.
• Ingenieros municipales de aguas pluviales modeladosinundaciones pluviales(encharcamientos de precipitación urbanos localizados) suponiendo un drenaje libre y sin obstáculos hacia los estuarios de marea receptores.

En los ciclones tropicales, los diluvios fluviales atmosféricos y las severas tormentas invernales del noreste, estos tres factores atacan simultáneamente. Los vientos huracanados provocan oleadas oceánicas catastróficas hacia la desembocadura de un estuario, elevando la elevación del agua de cola costera de tres a cinco metros. Al mismo tiempo, diluvios torrenciales arrojan pulgadas de precipitación por hora en todo el núcleo urbano, mientras que enormes hidrogramas de ríos de tierras altas corren río abajo hacia la costa. Cuando el oleaje del océano aguas abajo bloquea físicamente la descarga del río para que no escape al mar, se produce un efecto de remanso extremo tierra adentro, produciendo inundaciones catastróficas en las llanuras costeras bajas.

A finales de 2026, la ingeniería hidráulica costera ha logrado un avance monumental a través deModelado de inundaciones compuesto multifísico totalmente acoplado. Al vincular dinámicamente los motores de escorrentía hidrológica de las tierras altas, el enrutamiento hidrodinámico de inundaciones terrestres en 2D y elModelo de acción de olas y marejadas ciclónicas costeras de ADCIRC+SWANA través de una red computacional continua y sin interrupciones, los ingenieros pueden simular dinámicas de inundaciones compuestas con una fidelidad física sin concesiones.

Este completo tratado técnico examina las matemáticas que rigen, las arquitecturas de acoplamiento de límites, la física de las mareas estancadas y las estadísticas de probabilidad conjunta necesarias para modelar los riesgos compuestos de inundaciones costeras.

flowchart TD
  Meteo["Meteorological Forcing\nNOAA GFS / HRRR Atmospheric Pressure P_s(x,y,t)\n10-Meter Wind Velocity Vector W_{10}(x,y,t)"]
  ADCIRC["ADCIRC Ocean Solver\nGeneralized Wave Continuity Equation (GWCE)\nSea surface elevation ζ(x,y,t) & currents (U, V)"]
  SWAN["SWAN Wave Model\nAction Density Spectrum Balance N(σ, θ)\nRadiation stress tensors S_{xx}, S_{xy}, S_{yy} & wave setup ∇·S"]

  Meteo --> ADCIRC
  Meteo --> SWAN
  ADCIRC <-->|Two-Way Dynamic Coupling| SWAN

La ecuación de continuidad de onda generalizada (GWCE)

ADCIRC discretiza la cuenca oceánica a través de una malla triangular no estructurada de elementos finitos que se extiende desde las profundidades del Océano Atlántico o Pacífico hasta las plataformas continentales y las llanuras aluviales costeras. Para garantizar la estabilidad numérica sin amortiguación artificial, ADCIRC reemplaza la continuidad estándar con laEcuación de continuidad de onda generalizada (GWCE):

(∂^2 ζ) / (∂ t^2) + τ_0 (∂ ζ) / (∂ t) + ∇ · J - U H ∇ τ_0 = 0

donde:

• zeta(x,y,t)es la elevación de la superficie libre sobre el datum del geoide (NAVD88).
• H = zeta + hes la profundidad total instantánea de la columna de agua.
• tau_0es el parámetro de amortiguación numérico que suprime el ruido espacial espurio de alta frecuencia.
• jrepresenta vectores de flujo integrados en profundidad que incorporan aceleraciones de Coriolis, gradientes de presión atmosférica en la superficie, potencial de marea astronómico y tensión cortante del viento:

τ_s = ρ_{aire} C_d |W_{10}| W_ {10}

Acoplamiento bidireccional de onda-corriente con SWAN

Las olas de viento generadas por vendavales oceánicos no simplemente se desplazan sobre las marejadas ciclónicas; bombean activamente agua adicional a la costa a través degradientes de tensión de radiación de ondas (nabla · S). Cada 10 a 30 minutos, ADCIRC pasa elevaciones actualizadas de la superficie del mar.zetay los vectores actuales $(U, V)$ a SWAN, mientras que SWAN calcula la divergencia de la tensión de radiación y la transfiere de nuevo a ADCIRC como una fuerza impulsora de impulso adicional, resolviendo el nivel total del agua costera con precisión de centímetros.

3. Topologías de acoplamiento de límites dinámico: acoplamiento unidireccional, flojo y estrecho

La integración de las oleadas oceánicas de ADCIRC con motores 2D de cuencas hidrográficas urbanas y fluviales requiere una selección cuidadosa de la arquitectura de acoplamiento:

La ventaja de la malla continua sin costuras

A finales de 2026, las principales plataformas de modelado (como FloodSim Compound, HEC-RAS 2D con límites costeros y ADCIRC+SWAN no estructurado) favorecen laMalla no estructurada continua sin costurasenfoque.

Empleando elementos triangulares, cuadriláteros y poligonales adaptativos:

• En las profundidades del Océano Atlántico, los elementos son grandes (15 km de diámetro), donde domina la física simple de onda larga.
• A medida que la malla ingresa a las islas de barrera costera y las ensenadas de los estuarios, los elementos se refinan a 50 metros.
• A lo largo de las riberas de los ríos urbanos, las zanjas de drenaje y las calles de las ciudades costeras, los elementos se refinan hasta`1,5 a 3,0 metros`.

Debido a que no existe un límite artificial cero que divide el "océano" del "río", se eliminan las reflexiones numéricas y los errores de conservación. El modelo calcula naturalmente cómo las descargas masivas de los ríos rechazan las marejadas ciclónicas que se aproximan, capturando puntos de estancamiento límite exactos sin ajuste manual de parámetros.

4. Física de remansos estuarinos y parálisis del drenaje

Para entender por qué las inundaciones compuestas causan daños catastróficos al drenaje urbano, es necesario analizar la hidráulica deEmisarios de aguas pluviales alimentados por gravedad.

La física de la asfixia por remansos

En condiciones climáticas favorables, los conductos de aguas pluviales urbanas descargan en estuarios de marea a través de emisarios por gravedad equipados con compuertas abatibles para evitar la entrada de marea durante la marea alta. Cuando una marejada ciclónica extrema eleva el nivel del agua del estuario por encima de la corona del conducto (H_aumento} > z_{invertir} + D), se producen dos condiciones hidráulicas destructivas:

1. Cierre de puerta abatible y parálisis del sistema: La presión hidrostática del océano fuerza a cerrar la compuerta del emisario. La escorrentía pluvial interna recogida en la red de tuberías no tiene adónde ir. En cuestión de minutos, la red de tuberías subterráneas se presuriza, lo que provoca que los depósitos de captación a nivel de la calle se desborden y acumulen aguas pluviales a profundidades superiores a un metro, ¡incluso si el oleaje físico del océano nunca traspasa el malecón costero!
2. Restricción de descarga del orificio sumergido: Si el emisario no tiene compuertas, el nivel del agua del océano impone una enorme cabeza de cola aguas abajo (h_ {tw}), asfixiando la capacidad de descarga del emisario:

Q = C_d A √(2g (h_{aguas arriba) - h_{tw})}

comoh_ {tw}enfoquesh_ {corriente arriba}, la descarga cae hacia cero (Q → 0), paralizando completamente las redes de drenaje urbano.

5. Estadística de probabilidad conjunta y cópula bivariada en riesgo compuesto

La evaluación del riesgo de inundación compuesto requiere ir más allá de los períodos de retorno tradicionales de una sola variable (por ejemplo, "La inundación de los 100 años"). Una inundación compuesta de 100 años rara vez es causada por un evento de lluvia de 100 años que coincide con una marejada ciclónica de 100 años; más bien, es frecuentemente producido por elcoincidencia de una marejada ciclónica moderada de 10 años y un diluvio de lluvia moderado de 10 años, cuya ocurrencia conjunta crea una profundidad de inundación catastrófica de 100 años.

Modelado de dependencia no lineal mediante cópulas de Arquímedes

Debido a que las marejadas ciclónicas y las precipitaciones extremas son impulsadas por los mismos sistemas meteorológicos sinópticos (huracanes y depresiones tropicales), exhiben fuertesdependencia de la cola. Los coeficientes de correlación lineal estándar de Pearson no capturan esta agrupación no lineal.

Los ingenieros de riesgos hidráulicos utilizan elCópula de Gumbel-Hougaard, que exhibe dependencia de la cola superior (lambda_U > 0):

C_θ(u, v) = exp ≤ft( - ≤ft[ (-ln u)^θ + (-ln v)^θ ]^{1/θ} )

donde:

• u = F_X(x), v = F_Y(y)representan las transformaciones de probabilidad marginal uniforme.
• $\theta \ge 1$ representa el parámetro de dependencia de la cópula (calibrado mediante la correlación de rangos de Kendalltau = 1 - theta^{-1}).

Muestreando pares de eventos sintéticos(Q_i, zeta_i)Desde la cópula ajustada, los ingenieros ejecutan conjuntos de Monte Carlo a través de solucionadores hidrodinámicos 2D acoplados, construyendoCurvas "Y" del período de retorno bivariado (T_ {Y}) que cuantifican el verdadero riesgo compuesto con rigor estadístico.

6. Matriz comparativa integral de arquitectura de inundaciones compuestas

La siguiente tabla contrasta los enfoques de modelado aislados con el marco de simulación compuesto totalmente acoplado diseñado por FloodSim:

7. Estudio de caso de ingeniería del mundo real: inundación por huracán en el delta de un río de marea

Para demostrar la necesidad crítica del modelado compuesto acoplado, FloodSim ejecutó un estudio de validación hidrodinámica compuesta posterior al evento para un delta de un río costero de 120 kilómetros cuadrados que alberga un puerto comercial activo de aguas profundas, terminales químicas industriales y 185.000 residentes.

La dinámica compuesta revelada

La simulación combinada reveló que el desastre máximo no fue causado únicamente por el oleaje del océano:

1. El punto de estancamiento hidrodinámico: En la hora 11 de la tormenta, la marejada oceánica que avanzaba hacia adentro se encontró con la ola de descarga del río de 1.800 m³/s que avanzaba hacia afuera en River Mile 8.4. Los vectores de impulso opuestos crearon un límite de estancamiento estacionario donde la velocidad del flujo cayó a cero (tu ≈ 0).
2. Aumento de elevación del remanso: Al no poder escapar al mar, la descarga del río se acumuló detrás de la presa de marea, lo que obligó a que los niveles de agua en un parque industrial del interior alcanzaran`+5,27 metros NAVD88`—un completo1,42 metros más altoque la marejada ciclónica medida en la costa!
3. Remediación de defensa de infraestructura: Utilizando el modelo acoplado, los ingenieros portuarios diseñaron un sistema activo2,2 metrosdique de retroceso y un dique de alta capacidad80m^3/sestación de bombeo pluvial de flujo axial en la confluencia del estuario, salvaguardando exitosamente más1.200 millones de dólaresen activos de procesamiento petroquímico contra futuros huracanes de categoría 3+.

8. Protocolo de ingeniería de cinco fases para la ejecución de auditorías de inundaciones compuestas

Las empresas de ingeniería civil y costera que auditen el riesgo de inundaciones de estuarios en 2026 deben cumplir con este protocolo de ejecución de cinco fases:

1. Establecer datos de elevación geodésicos sin interrupciones: Armonizar todos los sondeos batimétricos (estudios hidrográficos de la NOAA) y la topografía LiDAR de tierras altas en un único datum vertical continuo (por ejemplo, NAVD88 con modelos geoides actualizados de la época 2026). Nunca mezcle la bajamar media inferior (MLLW) local con puntos de referencia terrestres.
2. Construir una malla multiescala no estructurada: Diseñar una malla continua que se extienda desde la plataforma continental profunda hasta los centros urbanos, refinando la densidad de nodos a lo largo de islas barrera, canales de navegación y crestas de diques.
3. Ejecutar un acoplamiento bidireccional entre olas y oleaje costero: Asegúrese de que ADCIRC pase niveles de agua instantáneos a SWAN cada 15 minutos, alimentando las tensiones de radiación de las olas nuevamente al solucionador de impulso hidrodinámico.
4. Implementar precipitación hidrológica directa sobre malla: En lugar de aplicar hidrogramas fluviales agrupados en los límites aguas arriba, aplique precipitación de radar NEXRAD distribuida espacialmente directamente en las celdas computacionales 2D, teniendo en cuenta las pérdidas dinámicas por infiltración del suelo mediante ecuaciones de Green-Ampt o Horton.
5. Calibrar la probabilidad bivariada conjunta: Ajustar las distribuciones de valores extremos marginales a los registros históricos de mareógrafos y caudalímetros, utilizando cópulas de Arquímedes para establecer límites rigurosos del período de retorno de múltiples riesgos.

A través de la integración disciplinada de modelos de circulación oceánica, hidrodinámica de aguas poco profundas y estadísticas de cópula multivariadas, los modelos de inundaciones compuestos acoplados permiten a las comunidades costeras anticipar, resistir y adaptarse a las tormentas multirriesgo más extremas del siglo XXI.