Escenarios estocásticos de precipitación y Evaluación
de la amenaza por inundación
Juan Camilo Olaya G.
juanog88@[Link]
27.05.2015
POR QUÉ ESCENARIOS ESTOCÁSTICOS?
Incertidumbre sustancial en cuanto a la caracterización de
la distribución temporal y espacial de la lluvia.
Donde nos
encontramos hoy
Condiciones futuras
plausibles
Evento Disruptivo
Presente Futuro
(Timpe & Scheepers, 2003)
POR QUÉ ESCENARIOS ESTOCÁSTICOS?
• Debido a que la ocurrencia de eventos a lo
largo del tiempo no puede ser pronosticada • Mutuamente excluyentes
con exactitud, y la ventana de tiempo de
análisis es desconocida, un set de eventos • Colectivamente exhaustivos
estocásticos es modelado. Cada uno de ellos
caracterizado con su frecuencia anual de • Admiten representación probabilista
ocurrencia.
Varias medidas de intensidad
Momentos estadísticos
Frecuencia anual de ocurrencia
Cuenca del río Chinchiná – Amenaza por lluvias
intensas
Número total de escenarios: 200
Número de medidas de intensidad: 1
Profundidad de Precipitación [mm]
FORMATO .AME
1. ESCENARIOS
ESTOCÁSTICOS DE
PRECIPITACIÓN
CLASIFICACIÓN DE MODELOS ESTOCÁSTICOS
DE PRECIPITACIÓN
(McMahon et al., 2000)
Categoría del modelo Clasificación del modelo Sub - clasificación
Procesos puntuales simples (e.g Poisson)
Modelos 1D puntuales Procesos agrupados
Procesos híbridos
Cadenas de Markov
Modelos multiestación
No-paramétricos
Modelos estocásticos de precipitación Celdas Agrupados
Banda modificada
Modelos espacio-temporales
Cascada aleatoria no homogénea
Radares
Modelos puntuales
Modelos de desagregación
Redes neuronales
Análisis regional de frecuencia de
precipitación
Profundidad-Duración-Frecuencia Estimación de reducción de área
Modelos sintéticos de tormentas de
diseño Modelación temporal
Transposición estocástica de
tormenta
MODELO ANALÍTICO PROPUESTO (CAPRA)
(ERN-AL., 2010)
MODELO ANALÍTICO PROPUESTO (CAPRA)
1. Análisis puntual (comportamiento puntual en sitios
de estaciones pluviométricas):
Curvas IDF
Curvas PDF
2. Análisis espacial (distribución espacial de la
precipitación):
Patrones espaciales
Curvas PAD
Curvas PADF
MODELO ANALÍTICO PROPUESTO(CAPRA)
20 Estaciones seleccionadas de 41 disponibles
MODELO ANALÍTICO PROPUESTO (CAPRA)
Procedimiento de completitud de información pluviométrica
DATOS INCOMPLETOS DATOS COMPLETADOS
Modelo analítico propuesto (CAPRA)
Procedimiento de completitud de información pluviométrica
Servicio de información ambiental de EEUU
De 4 estaciones cercanas y lo más uniformemente espaciadas
dentro de la región hidrológicamente homogénea
Modelo analítico propuesto (CAPRA)
Curvas PADF (Profundidad – Área – Duración –
Frecuencia)
Análisis espaciales de precipitaciones máximas
Éstos son análisis de precipitación máxima realizados
conjuntamente con todas las estaciones de medición de lluvia
disponibles, para fechas y horas concurrentes. Es decir, son
análisis de eventos históricos de aguaceros fuertes.
Para esto, a partir de las series históricas de profundidad
de precipitación en cada estación:
a) Elaboración de mapas de isoyetas de eventos de lluvia
intensa registrados en la zona de estudio para varias duraciones
Modelo analítico propuesto (CAPRA)
Curvas PADF (Profundidad – Área – Duración –
Frecuencia)
Se seleccionan los eventos extremos a partir de criterios simples de
selección, por ejemplo:
Del total de días con información diaria concurrente en la
mayoría de estaciones, seleccionar los días con lluvia diaria
que cumplen:
… más de 70% de estaciones con lluvia > 0
… precipitación máxima registrada en estas estaciones >
50 mm
… precipitación promedio > 25 mm
El resultado es un conjunto de aguaceros históricos, para
cada año, que pueden considerarse como los más intensos
Modelo analítico propuesto (CAPRA)
Curvas PADF (Profundidad – Área – Duración –
Frecuencia)
Para cada evento… 0
2.3
73,2 104.5
0
33.4
126.3
64.9
102.9
123.7
103.0
105.6
75.1
117.5
0
37.9
41.8
0
72.6
Modelo analítico propuesto (CAPRA)
Curvas PADF (Profundidad – Área – Duración –
Frecuencia)
b) Determinar para cada evento su correspondiente curva PAD
(Profundidad – Área – Duración)
Pmax (mm) Pmax (mm)
Area (km 2 )
22/05/1987 01/05/1989 140
4.6 98.3 131.8
Prec. Max. promedio (mm)
23 94.2 127.3 120
01-05-89
50 88.3 120.1 100
101 81.7 110.8
80
203 77 98.9
300 74.8 91.9 60 22-05-87
500 71.8 82.9
1000 64.4 72.4 40
2000 59 63.7 20
5000 49.6 53.7
10000 41.5 46 0
20000 34.7 40.2 1 10 100 1000 10000 100000
30000 29.8 36.7
Área (km2)
32607 28.4 35.7
Modelo analítico propuesto (CAPRA)
Curvas PADF (Profundidad – Área – Duración –
Frecuencia)
c) Dibujar las curvas PAD de cada año, con una duracuión D
específica. Y dibujar la envolvente superior
Curva PAD para un año y D particulares (Envolvente superior)
140
120
100
80
P (mm)
E1
60
E2
E3
40 E4
E5
20 E6
E7
0
0 200 400 600 800 1000
Área (km2 )
Modelo analítico propuesto (CAPRA)
Curvas PADF (Profundidad – Área – Duración –
Frecuencia)
d) Dibujar las envolventes superiores para todos los años de
registro y realizar el análisis de frecuencia para diferentes
valores de área
P (mm)
100 años
50 años
25 años
5 años
Área (km2)
Modelo analítico propuesto (CAPRA)
Curvas PADF Cuenca Chinchiná
PADF
PADF(3(1Días)
(5 Día)
800
250
1400
700
1200
200
Tr=2
600 Tr=2
1000 Tr=3
Profundidad [mm]
Tr=3
500
150
Tr=5
800 Tr=5
Tr=10
400 Tr=10
Tr=15
600 Tr=15
100 Tr=20
300
Tr=20
Tr=25
400 Tr=25
200 Tr=50
50 Tr=50
200 Tr=100
100 Tr=100
Tr=250
00
00 200
200 400
400 600
600 800
800 1000
1000 1200
1200
Área
Área[Km2]
[Km2]
Modelo analítico propuesto (CAPRA)
Malla de centros de tormenta
Análisis anual de isoyetas para las estaciones disponibles
Modelo analítico propuesto (CAPRA)
Patrones sintéticos espaciales
(ERN-AL., 2010)
Modelo analítico propuesto (CAPRA)
Patrones sintéticos espaciales
Modelo analítico propuesto (CAPRA)
Mapas de amenaza
Precipitación de
1 Día
TR 25 años.
P[mm]
Modelo analítico propuesto (CAPRA)
Mapas de amenaza
Precipitación de
1 Día
TR 50 años.
P[mm]
Modelo analítico propuesto (CAPRA)
Mapas de amenaza
Precipitación de
1 Día
TR 75 años.
P[mm]
Modelo analítico propuesto (CAPRA)
Mapas de amenaza
Precipitación de
1 Día
TR 100 años.
P[mm]
2. EVALUACIÓN DE
LA AMENAZA POR
INUNDACIÓN
CLASIFICACIÓN DE MODELOS DE
INUNDACIÓN
Consideración
No existe un modelo que pueda representar
hidráulicamente todas las cuencas, así como tampoco
existe un modelo que pueda caracterizar de manera
correcta todas las condiciones físicas y de flujo de
todos los cauces. Por lo tanto el modelador debe tener el
criterio necesario para que, con base en las
características particulares de su caso de análisis y las
características y calidad de la información con la que
cuenta, pueda seleccionar el modelo más adecuado.
CLASIFICACIÓN DE MODELOS DE
INUNDACIÓN
El análisis de flujos en ríos y en canales abiertos
generalmente es físicamente basado; sin embargo existen
modelos que no lo son. A continuación se listan los tres tipos
principales de modelos hidrodinámicos:
• Modelos estocásticos, físicos y de sensores remotos
• Modelos conceptuales.
• Modelos hidrodinámicos físicamente basados.
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Existen diversas maneras de clasificar los modelos y aproximaciones
de cálculo. En esta revisión se presenta una clasificación en cuanto a
la dimensión de análisis de cada uno de los modelos.
Fuente: FloodSite, 2009
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Para todos los modelos físicamente basados, se debe tener en cuenta
la ecuación de conservación de cantidad de movimiento de Navier –
Stokes:
Donde:
ρ, es la densidad del fluido
u, es la velocidad de flujo
t, tiempo
p, presión
µ, viscosidad del fluido
F, Factor de gravedad
Al combinar la anterior ecuación (Que describe la conservación del
movimiento en el fluido) con la ecuación de continuidad (Que describe
la conservación de la masa), se obtiene:
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Ecuaciones de flujo bi-dimensional (Saint-Venant)
h (uh) (vh)
0 Conservación de
t x y la masa
1 u u v h h h2 u u
u v 43
g t x x h x x Conservación de
momento
1 v v v h h h2 u u
u v 43
g t x y h y y
Descripción de la totalidad de flujo? Pérdidas?
Turbulencia? [10^-2 mm] – Discretización de malla infinita
Cuenca como sistema complejo (Evapotranspiración, meandros, planicie de
inundación, etc)
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos unidimensionales 1D: HEC-RAS, MIKE 11, ISIS SOBEK,
LISFLOOD-FF
Condición de flujo unidimensional en canales abiertos;
Solucionando las ecuaciones de St. Venant para la conservación
de cantidad de movimiento en una dimensión.
Donde:
Q, caudal de descarga
A, área de sección transversal
g, gravedad
Sf, pendiente de fricción
So, pendiente del cauce
Conservación de la masa
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos unidimensionales 1D: HEC-RAS, MIKE 11, ISIS SOBEK,
LISFLOOD-FF
Ya que las ecuaciones anteriores no tienen soluciones
analíticas, es preciso contar con valores de frontera
(Secciones transversales), que permitan su solución por
medio de métodos numéricos.
Condiciones de frontera: Dependiendo del
régimen de flujo, se establecen condiciones de
frontera “aguas arriba” o “aguas abajo”
Niveles limnigráficos conocidos
Profundidad crítica
Profundidad normal (Pendiente de la línea
de energía) (Propiedades morfométricas)
Curva de flujo vs elevación
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos unidimensionales 1D: HEC-RAS, MIKE 11, ISIS SOBEK,
LISFLOOD-FF
GEOMETRÍA
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos unidimensionales 1D: HEC-RAS, MIKE 11, ISIS SOBEK,
LISFLOOD-FF
RESULTADOS
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Topografía de detalle?
Discretización de la malla (Resolución
espacial)?
Costo computacional?
Su aplicación se recomienda para zonas
que comprenden, como máximo, 3 veces
el ancho del cauce principal (Pender,
2006)
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos bidimensionales 2D: RMA2, TELEMAC 2D, Hydro_as-2D,
LISFLOOD-ff, MIKE 21, TUFLOW, DELFT -HS
Retomando las ecuaciones de St. Venant
h (uh) (vh)
0 Conservación de
t x y la masa
1 u u v h h
h2 u u
u v 43
g t x x h x x Conservación de
momento
1 v v v h h
h2 u u
u v 43
g t x y h y y
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos bidimensionales 2D: RMA2, TELEMAC 2D, Hydro_as-2D,
LISFLOOD-ff, MIKE 21, TUFLOW, DELFT -HS
Q
Flujo bi-dimensional
t (hr)
1 u u v h u u h h
2
u v 43
g t x x h x x
1 v v v h u u h h
2
Dx u v 43
g t x y h y y
Se repite para (i,j+1) (i+1,j+1)
todos los
escenarios Dy M
(i,j) (i+1,j) h (uh) (vh)
0
t x y
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos bidimensionales 2D: RMA2, TELEMAC 2D, Hydro_as-2D,
LISFLOOD-ff, MIKE 21, TUFLOW, DELFT -HS
Proceso de cálculo Resultado
MODELOS HIDRODINÁMICOS
FÍSICAMENTE BASADOS
Modelos bidimensionales 2D: RMA2, TELEMAC 2D, Hydro_as-2D,
LISFLOOD-ff, MIKE 21, TUFLOW, DELFT -HS
Topografía de detalle?
Discretización de la malla (Resolución
espacial)?
Costo computacional?
Su aplicación se presenta en cuencas en
las que predomine la planicie de
inundación sobre la geometría del cauce.
EVALUACIÓN DE INUNDACIÓN
Esquema de funcionamiento de la herramienta ERN-Inundación.
(Ayuda ERN-Inundación, 2011)
EVALUACIÓN DE INUNDACIÓN
Cuenca tributaria generadora
de crecientes
Zona de afectación por
inundaciones
Estimar crecientes en cuencas tributarias y transitar éstas en zona
de inundaciones
EVALUACIÓN DE INUNDACIÓN
precipitación escorrentía
Cuenca
P
aguacero
Q
creciente
t
EVALUACIÓN DE INUNDACIÓN
• FACTORES DE ESCURRIMIENTO (CN)
Estimación del número de curva (CN) para diferentes condiciones
de suelo, a partir de:
1. Mapa georeferenciado de uso de la
tierra y/o cobertura vegetal
EVALUACIÓN DE INUNDACIÓN
• FACTORES DE ESCURRIMIENTO (CN)
Estimación del número de curva (CN) para diferentes condiciones
de suelo, a partir de:
2. Información de condición de la
superficie (IMÁGENES LANDSAT 5TM)
EVALUACIÓN DE INUNDACIÓN
• FACTORES DE ESCURRIMIENTO (CN)
Estimación del número de curva (CN) para diferentes condiciones
de suelo, a partir de:
3. Mapa georeferenciado de tipo de
suelo
EVALUACIÓN DE INUNDACIÓN
• FACTORES DE ESCURRIMIENTO (CN)
Teniendo en cuenta condiciones antecedentes secas:
EVALUACIÓN DE INUNDACIÓN
• FACTORES DE ESCURRIMIENTO (CN)
Teniendo en cuenta condiciones antecedentes intermedias:
EVALUACIÓN DE INUNDACIÓN
• FACTORES DE ESCURRIMIENTO (CN)
Teniendo en cuenta condiciones antecedentes húmedas:
EVALUACIÓN DE INUNDACIÓN
• Proyecto HEC-RAS. Definición de geometría
Herramientas de SIG para la determinación de la geometría del
cauce
Definición de:
• Distancia a la próxima sección transversal.
• N de Manning para las secciones del cauce.
• Límites izquierdo y derecho del cauce.
• Coeficientes de pérdida de energía:
contracción y expansión.
EVALUACIÓN DE INUNDACIÓN
• Proyecto HEC-RAS. Régimen de flujo
Régimen de flujo (Subcrítico, crítico o supercrítico)
Condiciones de frontera: Dependiendo del régimen de flujo, se
establecen condiciones de frontera “aguas arriba” o “aguas
abajo”
Niveles limnigráficos conocidos
Profundidad crítica
Profundidad normal (Pendiente de la línea de energía)
(Propiedades morfométricas)
Curva de flujo vs elevación
Caudales de entrada (Calculados a partir de análisis
Precipitación – Escorrentía)
Suelo Seco
TR 10 años. [m]
Suelo
Intermedio
TR 25 años. [m]
Suelo
Intermedio
TR 40 años. [m]
Suelo
Intermedio
TR 75 años. [m]
Suelo Húmedo
TR 100 años. [m]
Suelo Seco
TR 10 años.
[m/s]
Suelo
Intermedio
TR 25 años.
[m/s]
Suelo
Intermedio
TR 40 años.
[m/s]
Suelo Seco
TR 75 años.
[m/s]
Suelo Húmedo
TR 100 años.
[m/s]
