PRONSTICO DE CRECIDAS EN TIEMPO REAL XIMENA VARGAS

PRONÓSTICO DE CRECIDAS EN TIEMPO REAL XIMENA VARGAS M. PROFESOR ASOCIADO DEPTO. ING. CIVIL FACULTAD DE CS. FÍS. Y MAT. UNIVERSIDAD DE CHILE

AGENDA ANTECEDENTES GENERALES TRABAJOS PREVIOS A LA MODELACIÓN MODELOS DE PRONÓSTICO COMENTARIOS FINALES

Antecedentes Generales Inundación es la ocupación por parte del agua de zonas que habitualmente están libres de ésta. Causas: desbordamiento de ríos, superación de la capacidad de sistemas colectores en zonas urbanas, avalanchas, tsunamis, etc. Efectos: SANTA FÉ, 2003 VALPARAÍSO, 1906

Antecedentes Generales Estaciones de Monitoreo

Antecedentes Generales Estaciones de Monitoreo

Antecedentes Generales Sistema de Transmisión

TRABAJOS PREVIOS A LA MODELACIÓN Mapa de Riesgos

TRABAJOS PREVIOS A LA MODELACIÓN ESTACIONES DE MONITOREO METEOROLÓGICAS FLUVIOMÉTRICAS

TRABAJOS PREVIOS A LA MODELACIÓN ESTACIONES EN TIEMPO REAL

TRABAJOS PREVIOS A LA MODELACIÓN ANÁLISIS DE LA INFORMACIÓN: CANTIDAD CALIDAD PERÍODO CON INFORMACIÓN CONCURRENTE CONSISTENCIA DE LA INFORMACIÓN CARACTERIZACIÓN DE LA CUENCA

MODELACIÓN Modelos Hidrológicos: Distribuidos o Concentrados Caracterización de cuencas Requieren información detallada puede ser mejorada con espacial y temporal de técnicas modernas parámetros y variables

MODELACIÓN Modelos Meteorológicos

MODELACIÓN MODELO MM 5 Distribución temporal y espacial de la Precipitación.

MODELACIÓN MM 5: Características Generales • El Modelo entrega dos dominios de salida: – 45*45 km 2 – 15*15 km 2 • Se determinan 13 variables en cada celda de la grilla PPAC Presión Temperatura HR Td

MODELACIÓN MODELO SE INICIALIZA A LAS 20 HORAS Y LOS RESULTADOS SON CONOCIDOS A LAS 8 AM

MODELACIÓN SITUACIÓN EN TIEMPO REAL

MODELACIÓN SITUACIÓN EN TIEMPO REAL http: //mapas. snet. gob. sv/hidrologia/sat. phtml

MODELOS DE PRONÓSTICO MODELOS CONCEPTUALES MODELOS ARMA MODELOS DE FUNCIÓN DE TRANSFERENCIA MODELOS DE REDES NEURONALES LÓGICA DIFUSA

Modelos Conceptuales

Modelo Sacramento…. Modelo Swat Modelo VFlo

MODELOS DE PRONÓSTICO MODELOS ARMA (p, q) Son poco usados en tiempo real pues es difícil modelar la persistencia y responder, al mismo tiempo, a bruscas variaciones del caudal debido a una lluvia y al decaimiento gradual posterior.

MODELOS DE PRONÓSTICO

MODELOS DE PRONÓSTICO MODELOS de FUNCIÓN de TRANSFERENCIA o ARMAX ENTRADA SISTEMA RESPUESTA SI VARIACIÓN EN LA ENTRADA (X) PROVOCA RESPUESTA RETRASADA (Y) DEBIDO A LA INERCIA DEL SISTEMA Y SE ALCANZA UN NUEVO VALOR DE EQUILIBRIO CAMBIO ES RESPUESTA DINÁMICA MODELO DE FUNCIÓN DE TRANSFERENCIA ES AQUEL QUE DECRIBE ESTA RESPUESTA DINÁMICA

MODELOS DE PRONÓSTICO MODELO DE FUNCIÓN DE TRANSFERENCIA DISCRETO: OBSERVACIONES (X, Y) EN INTERVALOS EQUIESPACIADOS DE TIEMPO APARTE DE X, OTRAS VARIABLES AFECTAN A Y: PERTURBACIONES O RUIDOS MODELACIÓN DE SISTEMA REAL DEBE INCLUIR FUNCIÓN DE TRANSFERENCIA (MODELO DETERMINÍSTICO) Y MODELO DEL RUIDO

MODELOS DE PRONÓSTICO FILTRO DE KALMAN Xt: VARIABLE DE ESTADO Xt = F(Xt-1, t-1) + G(Wt, t) Ec. de estado o del Sistema Función de Transición

MODELOS DE PRONÓSTICO DE CRECIDAS PARA CONSTRUCCIÓN DE AUTOPISTA EN CAUCE RÍO MAPOCHO 4: QA ARRAYÁN EN LA MONTOSA 5: QM MAPOCHO EN LOS ALMENDROS

SISTEMA DE ALERTA DE CRECIDAS EN EL RIO MAPOCHO durante construcción costanera norte. La Montosa Los Almendros Cerro Calán Cerro San Cristóbal San Carlos Channel

QM Peak (m 3/s) Mean QM (m 3/s) t+2 t+6 t+12 t+20 Calibration 18. 7 3. 4 5. 9 19. 8 19. 2 50. 1 Validation 325. 0 21. 3 7. 9 20. 1 18. 1 50. 2 Data Set QA Peak (m 3/s ) Mean QA (m 3/s) Calibration 44. 8 14. 6 Validation 64. 2 9. 1 SE (m 3/s) t+2 t+6 t+12 t+20 3. 5 4. 5 5. 2 6. 0 3. 0 5. 9 8. 0 9. 1

MODELOS DE PRONÓSTICO FT Evento Q OBS [m 3/s] 100 500 1 105 éxito 2 600 éxito 3 95 falla éxito 4 470 éxito 5 267 éxito 6 422 éxito 7 95 falla éxito 8 89 falla éxito 9 188 éxito

SISTEMA DE ALERTA DE CRECIDAS EN EL RIO MAPOCHO durante construcción costanera norte.

MODELOS DE PRONÓSTICO REDES NEURONALES Se enmarcan dentro del área de la Inteligencia Artificial, cuyo principal objetivo es descubrir, describir y Capa oculta 1 Pesos de las Sinapsis capa oculta 2 capa de salida capa de entrada simular la inteligencia humana mediante métodos computacionales. NEURONA Los problemas que interesa estudiar son aquellos que no pueden expresarse mediante un algoritmo, este tipo de problemas tienen como característica común la EXPERIENCIA, la cual es usada por el hombre para resolver problemas.

MODELOS DE PRONÓSTICO Wi 1 Yj-1, 1 La Neurona UMBRAL Ii Wi 2 Yj-1, 2 Xi Yj SALIDA NEURONAL SINAPSIS Win Yj-1, n ESTADO DE EXCITACIÓN INTERNO

MODELOS DE PRONÓSTICO Funcionamiento de la Neurona. • Para iniciar la red, los pesos de las sinapsis toman valores pequeños. • La neurona comienza a operar cuando ocurre una Sinapsis hacia ella. • La neurona calcula su excitación interna Xi= S Wik Yj-1, k + Ii • Dónde el umbral pasa a ser W 0 k. Wi 1 Yj-1, 1 Xi Yj

MODELOS DE PRONÓSTICO • Luego la neurona calcula la salida como: Yi=F(Xi) • La función F puede ser una función no lineal para aprovechar la no-linealidad de la red neuronal.

MODELOS DE PRONÓSTICO • Otras Funciones típicas de transferencia

MODELOS DE PRONÓSTICO ENTRENAMIENTO = Proceso de Aprendizaje Para realizar el entrenamiento de la RN es necesario conocer pares de Entrada (xk) - Salida observados (ok) [Set de Ejemplos] que sean representativos del proceso que se quiere “enseñar” a modelar a la Red Neuronal. Función Objetivo: minimizar Sdk ok = Salida observada. yk = Salida de la Red. El proceso de Entrenamiento se basa fundamentalmente en métodos de Optimización (Min S ei), tales como el del Gradiente o los de Newton.

MODELOS DE PRONÓSTICO Retropropagación del error Yi Capa de salida di DWi Wji+1=Wji+DWi Oi Matriz de pesos de la capa j modificada

MODELOS DE PRONÓSTICO Inicializar todos los pesos con valores arbitrarios y pequeños. Presentar un estímulo de entrada Determinar la salida de la red Adaptar los pesos según el siguiente algoritmo, partiendo de la capa de salida hasta la primera capa oculta. ganancia

MODELOS DE PRONÓSTICO Si neurona en capa de salida Salida observada Si neurona en capa interna Salida de la neurona j Salida calculada por red Peso entre neurona j, k

MODELOS DE PRONÓSTICO Es importante mencionar que para realizar el Entrenamiento, los datos disponibles, se dividen en: 1. -Set de Entrenamiento [Se usan para ajustar los pesos y bias de la Red] 2. - Set de Validación. [Se usan durante el Entrenamiento para comprobar como se comporta la Red ante datos no “vistos”]. 3. -Set de prueba. [Se usa para ver el comportamiento de la Red ante situaciones nuevas ]. Para asegurar un buen ajuste se utiliza Early Stopping, el cual permite que el Modelo Neuronal se comporte en forma correcta al ser usado con datos “no vistos” en el proceso de Entrenamiento.

RNs han probado (Maier y Dandy, 2000; Dawson y Wilby, 2001) ser un mecanismo versátil tanto para modelar diferentes elementos del proceso lluvia-escorrentía como para el pronóstico de caudales (Campolo et al, 1999; Hu et al, 2001) Wilby et al (2003), entre otros aspectos, examinan el comportamiento de RNs para modelar el proceso lluvia – escorrentía y determinan la arquitectura que puede corresponder a un modelo parsimonioso de balance.

MODELOS DE PRONÓSTICO

1 día de adelanto Datos diarios de ríos de Turquía

MODELOS DE PRONÓSTICO SISTEMA DE ALERTA DE CRECIDAS EN EL RIO BIOBÍO.

MODELOS DE PRONÓSTICO . . . Capa de Salida (solo 1 neurona) Capa de Entrada ARQUITECTURA RED NEURONAL Pesos o conexiones

MODELOS DE PRONÓSTICO Datos disponibles y confección del set de datos requerido. - Biobío en Desembocadura - Biobío en Rucalhue Datos: eventos de crecidas (H Desembocadura 2. 6 [m]) ocurridos en los últimos 25 años. 29 crecidas pluviales: Entrenamiento : 15 crecidas => 5949 ejemplos (50 %) Validación Prueba : 9 crecidas => 2928 ejemplos (25 %) : 5 crecidas => 2973 ejemplos (25 %) Escalamiento de los datos considerando el valor máximo registrado en todas las crecidas.

MODELOS DE PRONÓSTICO Modelo para 1 paso de adelanto (k=1) Variando el número de neuronas en la capa oculta y observando aquella configuración de Red que entregue un menor error de validación. Con lo anterior se obtuvo una Red Neuronal [6 7 1]

MODELOS DE PRONÓSTICO Modelos Directos para 6, 10, 14 y 20 pasos de adelanto. El uso recursivo del modelo Directo a 1 paso , da origen a una alternativa para realizar pronóstico con hasta 20 horas de adelanto

MODELOS DE PRONÓSTICO Resultados obtenidos. Modelos Directos a 1, 6, 10, 14 y 20 pasos de adelanto. Modelo recursivo hasta 20 pasos de adelanto.

MODELOS DE PRONÓSTICO

MODELOS DE PRONÓSTICO Predicción hasta 20 pasos [Crecidas de Prueba]

MODELOS DE PRONÓSTICO Alertas De un total 32 Eventos, 25 Inundaciones: 32 Alertados y solo 7 Falsas Alarmas 29 Evacuaciones y solo 4 Falsas Alarmas 14 horas de anticipación entre Evacuación e inundación.

MODELOS DE PRONÓSTICO Validación Calibración 6 horas de adelanto Validación 6 horas de adelanto

PRONÓSTICO DE CAUDALES DE DESHIELO

Configuraciones Planteadas. Arrayán en La Montosa. RN (8, 15, 1)

MODELOS DE PRONÓSTICO Lógica Difusa Nace como una técnica perteneciente al área de la llamada inteligencia artificial y ha dado origen a sistemas expertos y sistemas de control automático, área en que se ha utilizado ampliamente. En particular, permite traducir o interpretar un sistema mediante un conjunto de reglas o hipótesis multivaluadas, a diferencia de la lógica binaria, puesto que no impone en las proposiciones o hipótesis valores discretos de las variables de decisión y consecuencia (falso o verdadero).

MODELOS DE PRONÓSTICO Lo difuso puede entenderse como la posibilidad de asignar valores de verdad intermedios entre “falso'' y “verdadero”. Además, este tipo de esquema es “tolerante” debido a que las causas (variables de entrada) que producen cierto efecto (variable de salida) pueden variar dentro de cierto rango sin que el resultado se altere mayormente

MODELOS DE PRONÓSTICO Estructura del FIS MF(k, j) Wi Outi MF 41(V 4) Output MF 31(V 3) ØOut 1 = W 1 * MFout 1 ØMFout 1 = ao 1 + a 11*(V 1) + a 21*(V 2) + a 31 *(V 3) + a 41*(V 4) MF 21(V 2) MF 11(V 1) ØW 1 Ø Ø = MF 11(V 1) * MF 21(V 2) * MF 31(V 3) * MF 41(V 4) Wi = MF 1 i(V 1) * MF 2 i(V 2) * MF 3 i(V 3) * MF 4 i(V 4) Out i = Wi * MFout i (V 1; V 2; V 3; V 4) MFout i = aoi + a 1 i*(V 1) + a 2 i*(V 2) + a 3 i *(V 3) + a 4 i*(V 4) Output = (Out 1+Out 2+Out 3+Out 4) / (W 1+W 2+W 3+W 4) Ø Q / Qfalla = Output / 1000

MODELOS DE PRONÓSTICO

MODELOS DE PRONÓSTICO Tormenta 1 Tormenta 8 Tormenta 3 Tormenta 5 Tormenta 2

MODELOS DE PRONÓSTICO Tabla comparativa ajuste de los diferentes modelos de pronóstico de caudales Modelo Error medio (%) calibración Error medio (%) validación FIS (pronóstico de 3 horas) 5, 34 6, 47 Modelo lineal (pronóstico de 3 horas) 4, 46 5, 55 FIS (aplicación recursiva) (pronóstico de 6 horas) 12, 9 10, 4 Modelo lineal (aplicación recursiva) (pronóstico de 6 horas) 6, 68 9, 3

COMENTARIOS FINALES AREA AFECTADA POR TORMENTA PASO DE TORMENTA SOBRE CUENCA RÍO LLUTA La temperatura del punto más alto de la superficie que registra el satélite, puede derivar de la radiación emitida desde la tierra o de una formación nubosa que cubre la zona observada. Esto puede definirse considerando que las temperaturas bajas son características de las nubes altas y que las temperaturas altas son explicadas por la reflexión de la radiación solar, típicamente alta en el norte del país en la época de verano.

C OME N T A R I OS F IN AL E S Izquierda: 1 Enero 2003 09: 39: 00. Derecha: 11: 39: 00 Izquierda: 1 Enero 2003 20: 39: 00. Derecha: 23: 39: 00

COMENTARIOS FINALES DETECCIÓN DE AREAS NIVALES EN DÍAS PREVIOS A LOS EVENTOS USANDO IMÁGENES PROGRAMA MODIS (MODERATE RESOLUTION IMAGING SPECTRORADIOMETER)

COMENTARIOS FINALES MODELACIÓN DISTRIBUIDA CUENCAS NIVALES SE REQUIEREN PRONÓSTICOS METEOROLÓGICOS

COMENTARIOS FINALES Cobertura Nival A) Observada, B) Simulada, C) Error. 18 de Mayo 2002

COMENTARIOS FINALES Espesor Simulado – Julio 2002

COMENTARIOS FINALES An Integrated Neuro-Fuzzy-Statistical Approach to Hydrological Modelling L. See, R. J. Abrahart y S. Openshaw de Universidad de Leeds Presentan 4 metodologías para generar en forma continua la escorrentía y efectuar pronósticos durante crecidas. Se aplican individualmente predictores estadísticos, modelos conceptuales, un modelo lingüístico difuso simple y redes neuronales con lógica difusa usando información de los ríos Ouse y Wye, UK. Luego, se integran usando promedio simple, aproximación Bayesiana y lógica difusa, obteniendo que la aproximación Bayesiana con lógica difusa resulta superior a las otras metodologías.

COMENTARIOS FINALES % veces en que el modelo se comportó mejor Cuenca Ouse en Skelton Modelo lógica difusa simple Entrenamiento Validación HNN 36. 2 33. 6 ARMA 21. 1 22. 2 SLFM 23. 5 26. 4 Naive 19. 2 17. 8 Cuenca Alto Wye Modelo 1984 1985 1986 TOPMODEL 1. 2 3. 0 2. 7 NN 1 4. 4 9. 7 8. 1 NN 2 6. 8 6. 6 6. 9 NN 3 8. 1 5. 9 5. 3 ARMA 60. 3 55. 3 60. 1 Naive 19. 2 19. 5 16. 9

COMENTARIOS FINALES Funciones de membresía para la entrada en Skelton Se recomienda el modelo HNN la mayor parte del tiempo, especialmente en la rama ascendente del hidrograma y algunas veces en la recesión. El modelo de Lógica Difusa Simple (SLFM) sólo se recomienda en la rama ascendente cuando los niveles son ALTOS o MUY ALTOS. Predicciones Naive se recomiendan durante el máximo o cuando el nivel cambia poco. El modelo ARMA se recomienda en la curva de recesión y cuando hay ascensos rápidos del nivel de agua en Dt de 6 horas.

COMENTARIOS FINALES %Veces en que Modelo Alerta [hrs. Promedio ] Modelo Tipo Alarma Temp (E) Temp (V) Correc (E) Correc (V) Tarde (E) Tarde (V) HNN Ind local 21 [1. 0] 30 [1. 7] 31 50 48 [1. 0] 20 [1. 5] área 14 [1. 3] 14 [1. 0] 48 72 38 [1. 9] 14 [1. 0] local 0 0 100 [3. 4] 90* [3. 8] área 0 0 100 [3. 6] 100 [3. 9] local 7 [2. 5] 10 [3. 0] 17 10 76 [3. 5] 80 [2. 8] área 0 0 10 29 90 [2. 7] 71 [3. 2] local 0 0 100 [6. 0] área 0 0 100 [6. 0] local 0 0 0 10 100 [3. 2] 90 [3. 1] área 0 0 100 [3. 4] 100 [3. 6] local 17 [1. 0] 30 [1. 3] 31 50 52 [1. 1] 20 [1. 5] área 9 [1. 0] 14 [1. 0] 48 57 43 [1. 6] 29 [1. 5] ARMA SLFM Naive Promedio FBM Ind Ind Multip Alarma local h= 3 m Alarma área h=3, 5 m

EN CADA SITUACIÓN DEBE ANALIZARSE CUAL ES EL MODELO MÁS ADECUADO SEGÚN INFORMACIÓN DISPONIBLE Y CARACTERÍSTICAS DEL PROCESO P-Q “…desatando una lluvia interminable sobre la ciudad, el valle y los cajones cordilleranos ahítos, provocando el 3 de junio (1783) una riada que alarmó a los vecinos transformados, desde ese momento, en temerosos vigías del nivel de las aguas. En los días siguientes la lluvia no mermó, sino que arreció de tal modo que el 16 de junio el Mapocho era un mar bravío, con color y olor de tierra, que bajaba de los cerros con ruido apocalíptico, como si los Andes estuvieran desplomándose”