Dam Management in the Era of Climate Change

Resumen Técnico: Gestión de Presas en la Era del Cambio Climático

Formulación del Problema
Este artículo aborda el control óptimo de una presa bajo las presiones duales del cambio climático y las restricciones operativas. El desafío central radica en la gestión de las entradas de agua que exhiben fluctuaciones extremas: periodos de sequía prolongados intercalados con eventos de lluvia intensos y agrupados ("bombas de agua"). Los modelos tradicionales que dependen del movimiento browniano no logran capturar estas dinámicas discontinuas y autoexcitables. En consecuencia, los autores modelan el problema de la gestión de la presa como un problema de control de conmutación estocástica óptima.

El estado del sistema se define mediante dos variables:

1. Nivel de Agua ($H_t$): Debe permanecer dentro de los límites de seguridad (de $h^-$ a $h^+$) para evitar el fallo de la presa (desbordamiento) y asegurar la funcionalidad (requisitos turísticos y ecológicos).
2. Intensidad de la Lluvia ($\lambda_t$): Modelada como la intensidad de un proceso de Hawkes, capturando la naturaleza autoexcitable de la lluvia, donde los eventos pasados aumentan la probabilidad de eventos futuros.

Las variables de control incluyen:

• Estado de la Turbina ($I_t$): Un control de conmutación binaria (0 para cerrada, 1 para abierta) para generar electricidad.
• Flujo del Aliviadero ($\beta_t$): Una variable de control continua que regula la liberación de agua sin generación de energía.

El objetivo es maximizar los ingresos descontados esperados por la producción de energía, minimizando al mismo tiempo las penalizaciones asociadas con niveles de agua por debajo del umbral turístico ($h^-$), por encima del umbral de seguridad ($h^+$), o por no cumplir con los requisitos mínimos de caudal ($\mu$) durante periodos de baja intensidad. El problema termina si el nivel del agua alcanza un punto crítico de fallo ($h_{max}$), incurriendo en una penalización terminal.

Metodología
Los autores proponen un marco matemático que se aparta de los modelos de difusión continua utilizando un proceso de Hawkes marcado para impulsar la dinámica del nivel del agua. La entrada de agua se modela como un proceso de salto compuesto donde la intensidad $\lambda_t$ evoluciona según un proceso de reversión a la media con saltos desencadenados por eventos de lluvia.

• Ecuaciones Diferenciales Estocásticas (SDEs): La dinámica del nivel del agua está gobernada por una SDE discontinua impulsada exclusivamente por el proceso de Hawkes marcado, sin un componente browniano. Esta elección se justifica por la naturaleza de las presas de montaña con cuencas pequeñas, donde las lluvias locales causan cambios abruptos en el nivel.
• Formulación de Conmutación Óptima: El problema se plantea como un problema de control mixto regular y de conmutación. Las funciones de valor $v_0$ y $v_1$ (correspondientes a la turbina inicialmente apagada o encendida) se definen sobre el espacio de estados $D = [0, h_{max}] \times [b, +\infty)$.
• Soluciones de Viscosidad: Debido a la falta de regularidad a priori en las funciones de valor (una consecuencia de la dinámica de saltos y los controles de conmutación), los autores caracterizan la solución utilizando soluciones de viscosidad de un sistema de desigualdades de integrales diferenciales parciales de Hamilton-Jacobi-Bellman (HJB). Establecen un Principio de Programación Dinámica (DPP) débil y demuestran que las funciones de valor son las soluciones de viscosidad únicas para el sistema HJB asociado.
• Aproximación Numérica: Para resolver las ecuaciones HJB, los autores emplean un esquema de diferencias finitas adaptado de problemas de cadenas de Markov controladas (Kushner y Dupuis). Esto implica discretizar el espacio de estados e iterar hasta la convergencia para aproximar las funciones de valor y las políticas óptimas.

Contribuciones Clave

1. Modelado Novedoso de la Entrada de Agua: El artículo introduce el uso de procesos de Hawkes marcados para modelar las entradas de agua de las presas, capturando explícitamente la agrupación de eventos de lluvia extrema y la duración de los periodos de sequía. Esto aborda un vacío en la literatura existente donde se suelen utilizar el movimiento browniano o procesos de Poisson estándar, los cuales no logran reproducir la naturaleza "autoexcitable" de los patrones meteorológicos impulsados por el clima.
2. Dinámicas Discontinuas: A diferencia de trabajos previos (por ejemplo, Chevalier et al. [16]) que modelaban los niveles de agua como procesos de difusión, este estudio utiliza un proceso de salto puro impulsado por la intensidad de Hawkes. Esto refleja mejor la dinámica de las pequeñas presas alpinas, donde la entrada de agua es altamente discontinua.
3. Caracterización de Soluciones de Viscosidad: Los autores proporcionan una caracterización rigurosa de las funciones de valor óptimas como soluciones de viscosidad de desigualdades integrales diferenciales, manejando las complejidades introducidas por el proceso de salto y el control de conmutación.
4. Análisis de Sensibilidad: El estudio investiga el impacto del parámetro de autoexcitación ($c$), que cuantifica el grado de agrupación en los eventos de lluvia. Este parámetro sirve como un indicador de la intensidad del cambio climático en el modelo.

Resultos Numéricos
Se realizaron experimentos numéricos utilizando parámetros representativos de presas de montaña (por ejemplo, altura < 100 m). Los hallazgos clave incluyen:

• Estructura de la Política Óptima: El control óptimo del aliviadero es de tipo "bang-bang", alternando entre tasas de flujo mínimas y máximas.
• Comportamiento de Conmutación: Para las configuraciones de parámetros específicas probadas (notablemente costos de conmutación altos), nunca es óptimo detener la producción de electricidad una vez que la turbina está funcionando. Sin embargo, la decisión de iniciar la producción depende del nivel del agua y la intensidad.
• Sensibilidad a la Autoexcitación: Un hallazgo crítico es que, a medida que el parámetro de autoexcitación $c$ aumenta (simulando una mayor variabilidad climática extrema), el nivel de agua óptimo dentro de la presa disminuye.
  • Interpretación: Aunque un $c$ más alto aumenta tanto el riesgo de sequía (que requiere niveles de agua más altos) como el riesgo de desbordamiento (que requiere niveles más bajos), los resultados numéricos indican que el riesgo de desbordamiento predomina. Los gestores responden manteniendo los niveles del embalse más bajos para crear capacidad de amortiguación para las repentinas "bombas de agua", incluso a costa de un mayor riesgo de sequía.
• Cambio de Rol: Los resultados sugieren un cambio de paradigma en la gestión de presas: bajo una alta variabilidad climática, las presas pueden funcionar cada vez más como mecanismos de protección contra inundaciones en lugar de ser la principal reserva de agua.

Significado y Reivindicaciones
El artículo afirma ser el primer estudio que integra explícitamente las restricciones físicas y los impactos del cambio climático (modelados mediante procesos de Hawkes) en un marco de optimización de la gestión de presas. Los autores enfatizan que los modelos de difusión estándar ya no son válidos para las zonas templadas que experimentan fluctuaciones meteorológicas extremas.

La importancia del trabajo reside en:

• Proporcionar una herramienta matemática para analizar las compensaciones entre la protección contra inundaciones y el suministro de agua en un clima cambiante.
• Demostrar que la estrategia de gestión óptima es altamente sensible a la naturaleza de agrupación de las lluvias, lo que sugiere que las "curvas de regla" tradicionales pueden necesitar revisión.
• Destacar que la naturaleza autoexcitable de los eventos extremos obliga a una postura de gestión conservadora (niveles de embalse más bajos) para mitigar los riesgos de fallos catastróficos, alterando potencialmente el rol económico primario de las presas de almacenamiento de energía/agua hacia la mitigación de inundaciones.

Los autores reconocen limitaciones, señalando que la intensidad de Hawkes no es directamente observable y requiere estimación a partir de datos hidrológicos, lo cual sigue siendo un desafío en un clima no estacionario. También señalan que el análisis de sensibilidad se centró principalmente en el parámetro de autoexcitación, dejando los estudios de sensibilidad económica y regulatoria más amplios para investigaciones futuras.