INGENIERÍA · AGUA · MODELOS DIGITALES DE ELEVACIÓN

Del terreno al vaso: cómo simular un embalse con un DEM

Un ejercicio técnico y auditable a partir del Embalse Puclaro: cota de agua, conectividad, barrera de presa, superficie, volumen y contraste con antecedentes conocidos.

La pregunta

Si conocemos la topografía y fijamos una cota de agua, ¿podemos reconstruir razonablemente la superficie que quedaría hidráulicamente conectada aguas arriba de una presa y estimar cuánto volumen almacenaría?

La respuesta es sí, con una condición importante: no basta con seleccionar “todo lo que está bajo una cota”. Hay que imponer la conectividad hidráulica y evitar que el agua atraviese el eje de presa.

Abrir simulador interactivo →

Caso de validación · Puclaro

Altura presa de referencia83 m
Capacidad de referencia200 hm³
Superficie de referencia760 ha
DEM de trabajo≈30 m

Estos valores son los parámetros de referencia incorporados en el código original de validación y no sustituyen antecedentes oficiales de diseño.

El problema parece sencillo: escoger una elevación y colorear los píxeles inferiores a ella. Pero ese procedimiento suele inundar depresiones desconectadas, valles vecinos e incluso sectores aguas abajo. La clave está en transformar el ejercicio topográfico en un problema de conectividad espacial.

1. La lógica del modelo

1 · DEMRepresenta la elevación del terreno mediante una grilla raster.
2 · CotaSe identifican celdas cuya elevación es menor o igual al nivel de agua elegido.
3 · MuroEl eje de presa se convierte en una barrera impermeable dentro de la grilla.
4 · SemillaDesde un punto aguas arriba se recorre solamente el terreno bajo cota conectado.
5 · MétricasSe integran área, profundidad y volumen de las celdas conectadas.

En el prototipo original desarrollado en Google Earth Engine, la conectividad se resolvió mediante una superficie de costo y cumulativeCost. En la versión web de DEITARIA el mismo principio se reproduce localmente en JavaScript mediante un recorrido de celdas vecinas sobre teselas de elevación.

2. Por qué la barrera del muro importa

Un DEM no “sabe” que existe una presa. Si el algoritmo ve celdas bajas a ambos lados del cierre, puede conectarlas. Por eso el modelo rasteriza una línea de muro con un ancho definido y la vuelve impermeable para el recorrido.

Interpretación: el muro virtual no intenta reproducir cada detalle estructural de la presa. Su función es topológica: separar el dominio aguas arriba del aguas abajo para que el vaso sea hidráulicamente coherente.

3. Qué se calcula para cada cota

IndicadorCálculo conceptual
Superficie inundadaNúmero de celdas conectadas × área de cada celda.
VolumenSuma de (cota de agua − cota de terreno) × área de celda.
Profundidad mediaVolumen / superficie.
PerímetroAproximación desde los bordes exteriores de la máscara inundada.
Altura hidráulicaCota del agua − mínima elevación DEM en el corredor del muro.

4. Validar antes de extrapolar

Puclaro sirve aquí como caso de control porque el modelo puede contrastarse con tres magnitudes de referencia incorporadas en el ejercicio: 83 m de altura de presa, 200 hm³ de capacidad y 760 ha de superficie máxima. El objetivo no es forzar matemáticamente el DEM para que entregue esos números, sino observar cuánto se aproxima con una metodología reproducible.

Una buena coincidencia simultánea en superficie y volumen es mucho más informativa que ajustar solamente una de las dos variables. Si la superficie coincide y el volumen no, el problema puede estar en las elevaciones del DEM, la definición del muro o la geometría de fondo. Si ocurre lo contrario, puede existir una envolvente espacial equivocada compensada por profundidades excesivas.

5. Qué significa —y qué no significa— la altura calculada

La “altura hidráulica calculada” del prototipo es una diferencia altimétrica obtenida del DEM: nivel de agua menos cota mínima detectada en el eje virtual. Es útil para diagnóstico y comparación, pero no equivale necesariamente a la altura estructural real de una presa, que depende de fundaciones, coronamiento, geometría de obra y criterios de diseño.

6. De un embalse conocido a alternativas nuevas

Una vez validado el procedimiento en un caso real, la misma estructura permite analizar emplazamientos potenciales: mover el punto de cierre y el eje de muro, escoger una semilla aguas arriba y construir curvas cota–superficie–volumen. Esa información permite comparar alternativas antes de entrar a etapas de ingeniería más costosas.

En DEITARIA este tipo de herramienta se entiende como un filtro de prefactibilidad y exploración técnica: permite visualizar rápidamente qué controla la capacidad de un valle, dónde aparecen saltos de volumen y qué tan sensible es una alternativa a la cota de operación.

7. Limitaciones que deben permanecer visibles

La precisión está condicionada por la resolución y calidad vertical del DEM, la representación simplificada del muro, la conectividad de la grilla, la ausencia de batimetría bajo el agua actual y las diferencias entre terreno observado y condiciones previas a la construcción. La herramienta no reemplaza topografía de detalle, estudios geotécnicos, hidrología de diseño, sedimentación, seguridad de presas ni evaluación ambiental.

8. Implementación web

La versión publicada utiliza Leaflet para el mapa interactivo, OpenStreetMap como base cartográfica y Terrain Tiles distribuidas como datos abiertos en AWS para disponer de elevaciones sin depender de una sesión autenticada de Google Earth Engine. El cálculo se ejecuta en el navegador, por lo que el prototipo puede alojarse como HTML + JavaScript convencional.

Resultado práctico: una idea que comenzó como script técnico en Earth Engine puede convertirse en una herramienta web pública, visual y auditable, manteniendo visible la relación entre topografía, cota, vaso y capacidad.

Fuentes y método

Metodología base: código de simulación/validación de Embalse Puclaro desarrollado para DEITARIA en Google Earth Engine. Cartografía: © OpenStreetMap contributors. Librería web: Leaflet. Elevación de la versión web: Terrain Tiles, Registry of Open Data on AWS / Tilezen. Los valores de 83 m, 200 hm³ y 760 ha corresponden a los parámetros de referencia usados en el ejercicio y deben contrastarse con documentación oficial si se usan fuera de este contexto de validación.