1.

RESUMEN

La modelación hidrológica en este caso se realiza para el estudio de inundación para un

tramo del Rio Macheta ubicado en el departamento de Cundinamarca Colombia. Se
utilizan datos reales obtenidos de tres estaciones pluviométricas cercanas descargados
de la página del IDEAM para obtener el hidrograma. Ya teniendo caudales se da inicio a
la modelación, mediante el programa HEC-RAS, versión 5.0.3. Primero, se genera un
Modelo Digital de Terreno, (DEM) con ayuda del Software QGIS3 3.16.1. En HEC-RAS se
realiza la modelación hidráulica, es decir se obtienen resultados como el perfil de flujo,
láminas de agua y numero de Froude. Luego de ello, se exportan resultados.

Primero se dedica un capítulo al marco teórico aplicable. Posteriormente se efectúa el

procedimiento completo para luego realizar un análisis del modelo obtenido y así finalizar
con las conclusiones del trabajo realizado.

2. PALABARAS CLAVE

Inundaciones, Cuenca hidrográfica, Caudal, Modelo Digital de Elevación (DEM), HEC-

RAS, QGis3, Raster, Coordenadas planas y geográficas, Magna Sirgas, Hidrograma.
Misión Topográfica Shuttle Radar (SRTM).

3. OBJETIVOS

3.1 OBJETIVO GENERAL

Realizar una simulación de inundación para un tramo del “rio macheta” ubicado en el
departamento de Cundinamarca Colombia.

3.2 OBJETIVOS ESPECÍFICOS

 Aplicar los conocimientos aprendidos en el manejo de las herramientas QGis3 y Hec-

RAS.
 Mostrar el procedimiento ejecutado paso a paso.
 Analizar el modelo y los resultados obtenido
 4. MARCO DE REFERENCIA

El “rio Negro”, pasa alrededor del municipio de Rio Negro del departamento de Santander-
Colombia ubicado a 19 km al norte de Bucaramanga. El municipio fue fundado el 5 de
mayo de 1805 por José Gutiérrez Calderón, Enrique Puyana, Juan Andrés Ortíz y los
hermanos José y Facundo Mutis con el nombre de Santa Bárbara de Rionegro. Tiene una
población de 27.114 habitantes.

Geografía

Descripción Física
El Municipio de Rio negro, está situado estratégicamente en la margen izquierda de la vía
que va de la capital de Santander a la Costa.
Es epicentro turístico por su cercanía con Bucaramanga, rico en fauna y flora, cuenta con
fuentes de agua; ríos y quebradas que son propicios para la recreación, el ecoturismo,
deporte de aventura y la preparación del típico sancocho familiar que se hacen los días de
fiestas y fines de semana por decenas de familias que se desplazan gracias a la facilidad
del transporte intermunicipal.
Rio negro además cuenta con un clima muy agradable, hermosas mujeres, magnifica
hospitalidad de sus gentes residentes en la cabecera Municipal que hace agradable la
estadía del turista que visita la población.

Límites del municipio

Limita por el Norte: Con los Municipios de La Esperanza - Norte de Santander, San
Alberto, San Martín (El Cesar) y El Playón.
Por el Occidente: Con los Municipio de Puerto Wilches y Sabana de Torres.
Por el Sur: Con los Municipios de Lebrija, Girón y Bucaramanga.
Por el Oriente: Con los Municipio de Matanza y Suratá.
Extensión total:1'277.457 km2 Km2
Temperatura media: 28°Cº C
Distancia de referencia: 20 Km de la capital de Santander, Bucaramanga

A continuación, se da una breve descripción de los programas a utilizar en la ejecución

del presente trabajo.

HEC-RAS
El diseño de obras hidráulicas para protección contra inundaciones, protección contra
socavación y/o erosión de obras de captación de agua para cualquier uso y obras de
cruces de cauces naturales o canales requieren un análisis delicado y detallado del
cuerpo de agua donde se realizarán para poder determinar qué tipos de obras civiles se
proponen según el objetivo final del proyecto. Este tipo de estudios se han llevado a cabo
a lo largo de la ingeniería en un principio con modelos matemáticos a mano, unos años
atrás, se han venido desarrollando y mejorando diferentes softwares que simulan el
comportamiento de dichas corrientes para disminuir el tiempo de ejecución de estas
labores y también para mejorar la calidad de los cálculos y así tomar mejores decisiones.
De los modelos más aprobados desde el punto de vista ingenieril, el Software HECRAS
es el más utilizado ya que además de contar con modelos 1D/2D, cuenta con la ventaja
de ser un software de libre acceso, volviéndolo muy accesible por los diferentes usuarios
a través de todo el mundo. Este software de modelación hidráulica fue creado por el US

2
Army Corps of Engineers (El cuerpo de ingenieros de la armada los Estados Unidos) el
cual ha sido integrado con otros tipos de softwares que se utilizarán en este curso; los
cuales son: QGIS v2.18.28 y su PLUG-IN QRAS; adicionalmente se usará AutoCAD como
plataforma de apoyo.

QGIS
QGIS es un Sistema de Información Geográfica (SIG) de Código Abierto licenciado bajo
GNU - General Public License. QGIS es un proyecto oficial de Open Source Geospatial
Foundation (OSGeo). QGIS proporciona una creciente gama de capacidades a través de
sus funciones básicas y complementos. Puede visualizar, gestionar, editar y analizar
datos, y diseñar mapas imprimibles

SRTM
La Misión Topográfica Shuttle Radar (acrónimo en inglés SRTM, de Shuttle Radar
Topography Mission) es un proyecto internacional entre la Agencia Nacional de
Inteligencia-Geoespacial, NGA, y la Administración Nacional de la Aeronáutica y del
Espacio, NASA. Su fin es obtener un modelo digital de elevación de la zona del globo
terráqueo entre 56 °S a 60 °N, de modo que genere una completa base
de mapas topográficos digitales de alta resolución de la Tierra.

EARTHDATA
La NASA estableció la División de Ciencias de la Tierra para comprender y proteger
nuestro planeta mediante el uso de tecnologías avanzadas y nuestra vista desde el
espacio para estudiar la Tierra y mejorar el análisis y la predicción del cambio del sistema
terrestre.

Para cumplir mejor con estos objetivos, el programa Earth Science Data Systems (ESDS)
de la NASA promueve el intercambio completo y abierto de todos los datos, metadatos,
documentación, modelos, imágenes y resultados de investigación, y el código fuente
utilizado para generarlos, manipularlos y analizarlos.

El conocimiento científico es más sólido y procesable cuando se deriva de métodos

transparentes, rastreables y reproducibles. La reproducibilidad se basa en el acceso
abierto a los datos y software utilizados para lograr resultados. Cuanto mayor sea la
disponibilidad de estas herramientas, más rápida y eficazmente las comunidades de
usuarios podrán utilizar la información para responder a preguntas básicas de las ciencias
de la Tierra y proporcionar la base para desarrollar aplicaciones prácticas innovadoras en
beneficio de la sociedad.

3
 5. METODOLOGÍA

A partir de la cartografía obtenido del programa QGis3, se realiza un Modelo Digital de

Terreno o DEM establecido para el “rio Machetá”. Luego se hace un recorte mediante un
raster, herramienta misma del QGis3. Este último se guarda como archivo tipo Tiff para
poder luego abrir el archivo en Hec-Ras. Una vez exportado el archivo, se delimita el flujo
de rio mediante un polígono, se indica la parte baja y alta del rio, se ingresan los datos del
hidrograma o los caudales del rio, los cuales, particularmente para este trabajo fueron
asumidos. También, se ingresan las fechas de inicio y fin del modelo y demás valores de
entrada. Por último, si todo el procedimiento se realizó se forma correcta, se genera el
modelo en 2D del flujo en el rio.

6. DESARROLLO Y PROCEDIMIENTO

6.1 PROCEDIMIENTO EN QGIS

Primero se abre la aplicación y se descargan los complementos SRTM y

Quickmapservices.

Ilustración 1. Descarga de complementos. Fuente: Elaboración propia.

4
Una vez descargados los complementos, se activan las capas “Google Hybrid” y “Google
Terrain” en la opción “Web” de la barra de inicio y pestaña “QuickMapServices”.

Ilustración 2. Activar capas “Hybrid” y “Terrain”. Fuente: Elaboración propia.

Luego, crear la capa tipo shape.

Ilustración 3. Crear capa tipo shape. Fuente: Elaboración propia.

5
En la creación de capa, guardar el archivo, escoger tipo de geometría polígono y asegurar
la zona de referencia geográfica que aplique y dar en aceptar.

Ilustración 4. Datos creación de capa shape. Fuente: Elaboración propia.

Luego, seleccionar en propiedades del proyecto, el mismo sistema de referencia de

coordenadas utilizado en el paso anterior y el caso de estudio debe estar dentro del área
naranja.

Ilustración 5. SRC. Fuente: Elaboración propia.

6
Luego, se dibuja la cuenca y se guardan los cambios realizados en la capa shape.

Ilustración 6. Dibujar la cuenca. Fuente: Elaboración propia.

Luego, las ubicaciones de las estaciones pluviométricas se crean como una capa de tipo
delimitado por comas, dando clic en capa-añadir capa- añadir capa de texto delimitado.
Se selecciona el archivo donde está la información de las cuencas, en el campo X se
selecciona Longitud en Y va Latitud y en Z la precipitación máxima (la cual se explica en
el siguiente subcapítulo como obtenerla), el sistema de coordenas es de donde se
obtuvieron es

Ilustración 7. Ubicación de las estaciones pluviometricas. Fuente: Elaboración propia.

La ubicación de las estaciones y de la cuenca a analizar atraviesa el municipio de Manta y

Macheta Cundinamarca.

7
 Ilustración 8. Ubicación de estaciones y cuenca. Fuente: Elaboración propia.

Luego, se activa la capa de elevaciones digitales SRTM

Ilustración 9. Descargar elevaciones del earthdata. Fuente: Elaboración propia.

Para descargar los datos es necesario inscribirse en la página web earthdata de la Nasa,
dando clic en el link que aparece en el cuadro emergente. Después de realizar el registro,
colocar usuario y contraseña y dar clic en aceptar.

8
 Ilustración 10. Inscripción en earthdata. Fuente: Elaboración Propia

Luego, se genera una capa SRTM

Ilustración 11. Descarga DEM. Fuente: Elaboración Propia

9
Luego, realizar el corte del raster a partir de una capa por máscara. Seleccionar la capa
de entrada y la capa de máscara y dar clic en ejecutar.

Ilustración 12. Extraer raster. Fuente: Elaboración Propia

Por último, se debe exportar el DEM creado como un archivo de tipo “GeoTiff” y después
de haber seleccionado la carpeta para guardar el archivo, escoger el sistema de
coordenadas EPSG:3116 y dar clic en aceptar.

Ilustración 13. Exportar DEM. Fuente: Elaboración Propi

10
 6.2 PROCEDIMIENTO EN HEC-RAS

Se debe crear una carpeta donde se almacene los archivos creados en QGis como el
DEM y el proyecto como tal que se ejecutará en Hec-Ras. Luego de haber guardado el
proyecto HEC-RAS, se debe abrir el archivo Tiff. Es importante recordar que antes de
empezar con el modelo en Hec-Ras se debe cambiar la configuración de región, dejando
como símbolo decimal el punto, y como símbolo de separación de miles la coma.

Ilustración 14. Apertura del archivo tipo Tiff en Hec-Ras. Fuente: Elaboración Propia

El modelo de la cuenca en Hec-Ras se puede apreciar con elevaciones.

11
 Ilustración 15. Modelo en Hec-Ras. Fuente: Elaboración Propia

Luego, se cierra el crea el modelo, se da clic en “edit” y luego en “Geometric data”. En

este punto se debe realizar un polígono que encierre el rio o quebrada principal, se crea el
inicio y fin del flujo. El polígono creado es el siguiente, se guarda el archivo y se cierra.

Ilustración 16. Geometría del flujo. Fuente: Elaboración Propia.

12
El siguiente paso es crear el hidrograma, para esto se da clic en “edit” y luego en
“Unsteady Flow data”. La línea de aguas arriba tendrá una condición de flujo hidrográfico,
y la línea de aguas abajo tendrá una condición de profundidad normal.

Ilustración 17. Parametrización del hidrograma. Fuente: Elaboración Propia.

Los datos utilizados para el hidrograma propuesto son:

Tabla 2. Datos para hidrograma. Fuente: Elaboración Propia.

Tiempo Caudal Tiempo Caudal Tiempo Caudal

(hr) (m3/s) (hr) (m3/s) (hr) (m3/s)
1 120.0 41 171.6 81 147.0
2 121.0 42 173.5 82 144.0
3 122.0 43 175.5 83 141.0
4 123.0 44 177.4 84 138.0
5 124.0 45 179.4 85 135.0
6 125.0 46 181.4 86 132.0
7 126.0 47 183.3 87 129.0
8 127.0 48 185.3 88 126.0

13
 Tiempo Caudal Tiempo Caudal Tiempo Caudal
(hr) (m3/s) (hr) (m3/s) (hr) (m3/s)
9 128.0 49 187.2 89 123.0
10 129.0 50 189.2 90 120.0
11 130.0 51 191.2 91 119.0
12 131.0 52 193.1 92 118.0
13 132.0 53 195.1 93 117.0
14 133.0 54 197.0 94 116.0
15 134.0 55 199.0 95 115.0
16 135.0 56 201.0 96 114.0
17 136.0 57 202.9 97 113.0
18 137.0 58 204.9 98 112.0
19 138.0 59 206.8 99 111.0
20 139.0 60 210.0 100 110.0
21 140.0 61 207.0 101 109.0
22 141.0 62 204.0 102 108.0
23 142.0 63 201.0 103 107.0
24 143.0 64 198.0 104 106.0
25 144.0 65 195.0 105 105.0
26 145.0 66 192.0 106 104.0
27 146.0 67 189.0 107 103.0
28 147.0 68 186.0 108 102.0
29 148.0 69 183.0 109 101.0
30 150.0 70 180.0 110 100.0
31 152.0 71 177.0 111 99.0
32 153.9 72 174.0 112 98.0
33 155.9 73 171.0 113 97.0
34 157.8 74 168.0 114 96.0
35 159.8 75 165.0 115 95.0
36 161.8 76 162.0 116 94.0
37 163.7 77 159.0 117 93.0
38 165.7 78 156.0 118 92.0
39 167.6 79 153.0 119 91.0
40 169.6 80 150.0 120 90.0

Teniendo en cuenta que es un ejercicio académico los valores del hidrograma no son reales, sin
embargo, se toma como referencia el caudal máximo en el rio macheta tomado del “ESTUDIO DE
IMPACTO AMBIENTAL PROYECTO UPME-03-2010, SUBESTACIÓN CHIVOR II – Y NORTE 230 kV Y
LÍNEAS DE TRANSMISION ASOCIADAS "

14
La pendiente a utilizar es 0.035 y se calcula dividiendo la longitud del cauce principal
sobre la diferencia de altura entre el punto localizado aguas arriba y el punto aguas abajo.

Ilustración 18. Cauce Principal Fuente: Elaboración Propia.

Para determinar está pendiente de forma sencilla se crea una capa de tipo shape y
longitud, se dibuja el cauce principal y se exporta en formato KML, este archivo se abre en
Google Earth y se da la opción de ver perfil longitudinal.

15
 Ilustración 19. Verificación del hidrograma. Fuente: Elaboración Propia.
Luego, se procede a ingresar los datos para correr el programa. Para esto, se da clic en
“Run” y luego en “Unsteady Flow Analysis”.

Ilustración 20. Ingreso de datos para correr el programa. Fuente: Elaboración Propia.

16
Por último, se verifica el modelo, nuevamente dando clic en “GIS Tools”. El modelo
obtenido fue:

Ilustración 21. Modelo Final. Fuente: Elaboración Propia.

17
 7. DISCUSIÓN DE RESULTADOS Y CONCLUSIONES

Ilustración 22. Profundidad Modelo Final. Fuente: Elaboración Propia.

Como se evidencia en la ilustración 20 la simulación empieza a las 00:00 horas del 01 de

agosto de 2020 hasta las 23:00 horas del 05 de agosto de 2020, es decir, nueve (05) días
aproximadamente.

En la herramienta RAS Mapper, el cual muestra la simulación del modelo generado, se

observa que la profundidad del rio se encuentra en 0 y 15 metros, con una velocidad de
flujo inferior a 15 m/s según las escalas de colores como se evidencia en las ilustraciones
22 y 23.

18
 Ilustración 23. Velocidad Modelo Final. Fuente: Elaboración Propia.

Por otro lado, se determina la fisiografía y los parámetros morfométricos de la cuenca a partir de
datos de entrada obtenidos como se muestra a continuación.

Primero en QGIS, una vez creada la capa shape de la cuenca, se da clic derecho sobre la capa en
la parte izquierda de la pantalla y luego en “Abrir tabla de atributos”. En el cuadro que aparece se
da clic en el lápiz y luego en “Abrir calculadora de campos”. Aquí es posible obtener el área y el
perímetro de la cuenca.

Par obtener la longitud del cauce principal, la cota de entrada y salida y la longitud de los cauces
secundarios, se crea una capa tipo shape y se dibujan los cauces en el mapa, esta capa se exporta
en formato KML. El archivo resultante es compatible y se puede abrir en Google Earth. En este
programa en perfil de elevación se pueden ver toda la información necesaria.

19
 Ilustración 24. Área y Perímetro de la cuenca

Ilustración 25. Longitudes y pendiente de cauces en Google Earth. Fuente: Elaboración Propia.

Los datos de entrada para calcular los parámetros morfometricos de la cuenca se muestran en la
siguiente tabla.
Parámetro Valor
Área (Ha): 105.61
Perímetro (km): 45.494
Longitud Cuenca 23.95
Cota Inicial Cauce Principal
2426
(m):
Cota final Cauce Principal (m): 1414
Longitud Cauce Principal (km): 24.4
Rama 1 (km): 3.49
Rama 2 (km): 3.58
Rama 3 (km): 4.47
Rama 4 (km): 1.87
Rama 5 (km): 1.52
Rama 6 (km): 0.73
Longitud Total Cauces (km): 40.06

Tabla 3. Parámetros de entrada cuenca Rio Macheta.

20
Los datos se ingresan en la página web [Link] en la sección hidrología y
se obtienen los siguientes resultados.

Ilustración 26. Ingreso de datos en página Web.

Resultados
Ancho de la cuenca (km): 4.3283
Factor de Horton (kf): 0.1774
Relación de elongación (Re) 0.4751
Coef. compacidad (Kc): 1.2395
Relación de circularidad (Rc): 0.6412
Densidad de drenaje (Dd): 0.3793
Const. estabilidad de cauce (Cc): 2.6363
Densidad de corrientes (Dc): 0.2632

Tabla 4. Resultados parámetros morfométricos cuenca Rio Machetá

De forma similar se calcula en la página web el tiempo de concentración y se obtiene los siguientes
resultados.

21
 Tiempo de Concentración
Temez (hr): 2.5951
California C-P (hr): 2.6442
Giandotti (hr): 3.0531
Pilgrim McDermott (hr): 4.465
Williams (hr): 7.0318
Ven Te Chow (hr): 2.0577
Kirpich (hr): 2.6413
Scs-Ranser (hr): 1.7122
Clark (hr): 13.6435
Passini (hr): 2.5167
Johnstone-Cross (hr): 43.1413
Pérez (hr): 2.2876

Tabla 5. Tiempos de concentración cuenca Rio Machetá

22
REFERENCIAS

Edgar O. Ladino Moreno. Slides diseño y construcción de canales.

[Link]

[Link]

[Link]
[Link]/informacion_general.shtml

Taller de simulación dinámica de inundaciones con Hec-RAS 5.0

[Link]

23