Estudio hidrológico del embalse de Velilla del Río Carrión

Trabajo sobre SIG
Título	Estudio hidrológico del embalse de Velilla, Río Carrión
Autores	Guillermo Alfaya Franklin Christian Balic Stefanovic Gonzalo Alonso Felices Beceiro Manuel Huerta Duque
Asignatura	Sistemas de Información Geográfica Aplicados a la Ingeniería Civil
Curso	Curso 18/19
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

2 Metodología

Las hojas utilizadas son la 0105, 0106, 0131, 0132 y 0081 de la colección de mapas del IGN 1:50000 ráster.

Los datos de las estaciones meteorológicas se han obtenido de la página de AEMET.

El procedimiento seguido ha sido el siguiente:

Combinación de las hojas ráster en una única capa utilizando el complemento de Clipper que proporciona QGIS 2.18.16. Las hojas utilizadas son la 105, 106, 131, 132 y 81. Crear directorio Grass para poder trabajar con el programa dentro de QGIS. Exportar a Grass la capa creada y obtener las subcuencas con un número mínimo de celdas para cada cuenca de 35000. Convertir la capa Cuencas en vectorial y exportarla a QGIS. De todas las subcuencas existentes, se escoge la correspondiente creando una nueva capa vectorial, usando la opción selección de objeto y el editor de capas vectoriales. Para extraer las subcuencas de la cuenca seleccionada se repite el proceso anterior, quedando 5 subcuencas en una única capa vectorial. Para el estudio hidrológico se necesitarán algunas propiedades geométricas de la cuenca. Se añade el área a la tabla de atributos de la cuenca. Las coordenadas de los pluviómetros y su elevación se extraen del Excel descargado de AEMET. Tras pasarlas a formato .txt se pueden introducir como una capa vectorial de puntos en el proyecto. Se realiza un buffer de 13000 m desde el centroide de la cuenca para seleccionar los pluviómetros más cercanos. Se crea una nueva capa vectorial con la herramienta de intersección entre las capas del buffer y de los pluviómetros. Con la herramienta vectorial polígonos de Voronoi se obtienen los polígonos de Thiessen. Se crean capas con cada una de las 5 subcuencas, las cuales intersectan con los polígonos, obteniendo así la parte de cada uno correspondiente a cada subcuenca. Se extraen las áreas a una hoja de Excel para poder calcular, dentro de cada subcuenca, qué porcentaje de precipitación aporta cada pluviómetro. Trabajando otra vez en Grass, se extraen los ríos de un tamaño mínimo de celda de 1000 unidades. Se convierte en líneas vectoriales dentro de Grass y se exporta a QGIS como capa vectorial. Se intersecta la capa de ríos con la cuenca y las subcuencas. Para poder obtener las propiedades del río en cada subcuenca se combinan los elementos para poder asociar cada tramo de río con una subcuenca. Se obtienen las longitudes de los ríos para exportarlas a Excel. Se hacen perfiles sobre la capa de MDT combinada y a partir de las elevaciones y de la longitud del río, se calculan las pendientes.

4 Conclusiones

A partir de los datos obtenidos se podría deducir la precipitación media y máxima diaria en cada subcuenca, a lo que habría que aplicar un coeficiente de reducción por área para poder asumir que la lluvia vaya a ser homogénea en cada subcuenca. Todo ello debe realizarse a partir de los datos de lluvia para cada pluviómetro usando la distribución estadística que mejor se adapte para conocer el valor en periodos de retorno altos. Por lo tanto, la obtención de estos datos resulta fundamental para comenzar con el estudio hidrológico y el cálculo de los hietogramas. Para el estudio de esta cuenca se han usado los pluviómetros más cercanos que se han obtenido al realizar un buffer de 13000 m, quedando 7 en el área de influencia. Sin embargo, sólo uno de ellos está dentro de la cuenca, por lo que el estudio habría sido más preciso si hubiese más pluviómetros dentro.

5 Anejos

Los mapas empleados durante el estudio han sido los siguientes: Se pueden adjuntar archivos usando el enlace Subir archivo que aparece a la izquierda.