Paso a paso de la actividad practica de HEC-RES con Qgis
Continuando con el curso Modelación Hidráulica con QGIS y Hec-RAS, el día de hoy
veremos modelación hidráulica en cauces naturales utilizando QGIS, RIVERGIS Y
HEC-RAS
1. Ejecutamos nuestro navegador
PostgreSQL: Windows installers (ver video de instalación)
2. Descargar Qgis en la siguiente Pagina o viendo el siguiente
video
https://qgis.org/es/site/forusers/download.html#
https://www.youtube.com/watch?v=oc3fzeTsDxA
3. Enviar 3 carpetas
En la carpeta HEC se exportarán los resultados, el ráster es la imagen
georreferenciada, y la shapefile en donde tenemos el área de estudio y las
curvas de nivel
4. Abrimos Qgis
5. Cambiamos el sistema de coordenadas
La zona de estudio se encuentra en wgs84 UTM ZONA 17 SUR. Para cambiar
esto se debe dirigir a la parte inferior derecha. Esta en Perú
Escribir el código 3271
Elegimos ese haciendo el cambio de proyección, escogimos ese porque abarca
la zona de estudio que es Perú
6. Cargamos capas
- Hacemos clic derecho a la barra de herramientas flotantes
7. Ponemos la carpeta como favoritos en donde tenemos toda la
información
}
- Arrastramos la imagen georreferenciada de la carpeta raster
- Arrastramos el área de estudio y la capa de curvas de nivel
- Damos un mejor estilo al área de estudio, quitar el relleno y que
solo se muestre el contorno
- Activamos el panel
- Seleccionamos la capa área, relleno sencillo, y en estilo de relleno
es sin relleno
- Color de marca
- Aumentar anchura de marca
- Ahora obtendremos el TIN a partir de las curvas de nivel, esta
interpolación lo tenemos que hacer en función de la elevación
- Activamos en panel
- Atributo de interpolación es en donde están ubicadas las alturas
- La ventaja de Qgis es que Podemos generar una capa temporal, es
decir, que lo podemos ejecutar para un proceso, y una vez que lo
realizamos lo podemos eliminar y no llenarnos de archivos
innecesarios y tenemos un mejor orden en nuestras capas
- Ejecutamos y obtenemos el TIN de las curvas de nivel
- Al superponer la capa área observamos que la interpolación a
sobrepasado al área que yo requiero
- Por lo que vamos a hacer un corte en función del Shape de área
- Como el archivo lo vamos a necesitar, lo guardaremos en la carpeta
raster, lo guardamos como DEM y ejecutamos y .tif
- El resultado es el DEM recortado con el área de estudio, Tenemos el
DEM en función del área de estudio
- Eliminamos el primer ráster
- Ahora empezaremos a trabajar con el complemento RiverGis
- Ejecutamos el complemento RiverGis
- La instalación del complemento riverGis y la conexión de datos
postGis ya se vieron
- Ahora vamos a ejecutar el complemento RiverGis, la base de datos
es la que hemos creado
- Creamos un nuevo esquema y lo llamamos modelo
- Creamos un nuevo esquema que se llame rivergis
- Al poner el esquema Tpology se han activado todos los iconos para
trabajar en Rivergis
- Creamos la base de datos de nuestras capas del rio
- En el caso de que tuviera diques, áreas inefectivas o obstrucción de
bloques los podemos seleccionar
- Observamos que se han creado las 4 capas que hemos seleccionado
anteriormente
- Cerramos riverGis
- Empezamos con la edición de cada capa
- Streamcenterlines
- Crearemos nuestro eje del río, para que sea más fácil y practico
- Hay que tener activado esto
- Configuramos la capa del DEM
- Le damos transparencia del 50%
- Ahora dibujaremos las líneas, el mismo proceso de como se trabaja
en HEC-RAS, las líneas deben ser dibujadas de aguas arriba a aguas
abajo
(la línea central de un curso de agua, como un río o un arroyo. Esta línea
representa el eje principal del flujo de agua)
- Dibujamos el eje del río con base a la imagen georreferenciada
- Observamos que se ha creado el streamcenterline y la etiqueta
- Damos clic en conmutar edición para guardar los cambios
- Ahora crearemos lo banklines (Los banklines son líneas que representan
los límites de los bancos de un río. En el contexto de la extensión RiverGIS
para QGIS, los banklines se utilizan para crear la geometría del modelo de
flujo HEC-RAS a partir de datos espaciales
RiverGIS es un complemento de QGIS que ayuda a crear el modelo de flujo
HEC-RAS. Permite estimar diferentes parámetros hidráulicos y definir un
entorno de simulación. Además, permite visualizar los resultados en RAS
Mapper y en QGIS para editar las capas de información obtenidas, como
mapas de velocidades, profundidad y zonas inundables )
- Desactivamos el DEM
- Conmutamos edición
- Clic en añadir línea
- Trazamos la primera línea de abajo hasta llegar al final, debe estar
dentro del DEM, luego al finalizar el primer tramo debemos dar clic
derecho y guardamos
- Repetimos el proceso para el borde superior de aguas arriba hacia
agua abajo
- Clic en conmutar edición para guardar los cambios
- Ahora editaremos los flowpaths àConmutamos edición à Añadir
línea
(los "flowpaths" podrían referirse a las trayectorias o rutas del flujo de agua
en un sistema fluvial modelado con RiverGIS.)
- En donde le daremos una posible llanura de inundación, el proceso
va a ser el mismo, las líneas deben ser creadas de aguas arriba a
aguas abajo, las trazamos à clic derecho à FpId se deja à Tipo de
línea se debe decir en que margen es la que hemos realizado el
dibujo, si estamos trabajando de aguas arriba a aguas abajo la línea
que hicimos se ubica en la margen derecha
- Para la otra parte hacemos lo mismo
- Se debe escribir en inglés para no tener problemas a la hora de
exportar los resultados
- En HEC-RAS para poder trabajarlo se debe tener tres líneas, margen
derecha, izquierda y eje del río à Si desactivamos el
StreamCenterLineas , observamos que nos falta el flowpath de la
mitad
- Para añadirlo los dirigimos a la extensión de RiverGis
- El flowpath añadido va de acuerdo con el eje de nuestro río
- Si abrimos la tabla de atributos del Flowpaths tenemos margen
derecho e izquiedo y el canal que viene siendo el eje de nuestro río
- Ahora crearemos los XSCutlines (Los "XSCutlines" en RiverGIS o QGIS
se refieren a las líneas de corte de secciones transversales. Estas líneas
representan las secciones transversales de un río o curso de agua y se
utilizan en el modelado hidráulico para definir la geometría del flujo.En
RiverGIS, las secciones transversales están representadas por la tabla
"XSCutlines". Esta tabla contiene atributos como los nombres de los ríos
o tramos, códigos de río, códigos de alcance y atributos relacionados con
las líneas de corte de las secciones transversales. Estas líneas de corte
se utilizan para calcular la posición de los bancos, las distancias a la
siguiente sección transversal aguas abajo a lo largo de un camino de
flujo, y se ejecutan diversas herramientas relacionadas con las
XSCutlines.)
Los bancos se refieren a las posiciones o ubicaciones a lo largo de las líneas de
corte de las secciones transversales de un río.
- A diferencia de HEC-RAS Las líneas de corte se deben realizar de
manera manual, recordad que debe ser de izquierda a derecha, se
empieza de aguas arriba a aguas abajo
- Realizada la edición de todas estas capas nuevamente damos clic
en el complemento Rivergis
- Añadimos la información topológica de mi río y de las secciones
- En este caso empezaremos por nuestro río, y le damos que realice
todo el cálculo topológico y la longitud
- Creo el punto 1 y el punto 2 que es el inicio y el final del cauce
-
- Observamos que se ha realizado el calculo y aparecen dos puntos
que viene a ser el inicio y el final de mi cauce
- En el caso de las secciones, haremos lo mismo
- Nos pedirá el ráster de la interpolación que es nuestro DEM, Le
damos seleccionar y aceptar
- Verificamos en la tabla de atributos del eje de ríos
streamcenterlines, tenemos la longitud total de nuestro cauce
- Ahora revisamos la tabla de atributo de los XSCutlines
- Observamos que nos calculó los Banks de la margen izquierda y
derecha y la longitud que hay en cada uno de ellos en ambas
márgenes
- Los datos que hay en 0, le debemos dar un orden, dando clic en
XsecID
- Viene siendo el punto o corte final que le hemos dado, como ya no
hay otra sección es por eso que lo toma como 0
- Ya luego observaremos nuestros resultados
- Observamos que se ha realizado la exportación
- Muestra el aviso de WARMING
- Es porque no tenemos esa información aparece así
- Nos dirigimos a la carpeta y vemos que se guardó, esto nos
permitirá realizar la simulación hidráulica con HEC- RAS y exportar
los resultados de modelación a nuestro Qgis
AHORA VEREMOS LA MODELACIÓN HIDRAULICA EN HECRAS Y EXPORTAR
LOS RESULTADOS A QGIS
- Ejecutamos el software HEC- RAS
- Creamos un nuevo proyecto
- Nos muestra donde lo guardamos y también el sistema de
unidades, como vamos a trabajarlo en sistema internacional lo
dejamos tal y como esta
- Importamos la geometría del video anterior (Lo que hicimos en
Qgis)
- Nos aparece lo siguiente:
- Que es la ventana de importación y ponemos el sistema
internacional métrico
- Lo dejamos como esta
- Esto también lo dejamos como esta, y le damos finalizar, importar
data
- Como resultado tenemos la línea de nuestro cauce principal, las
secciones trasversales y los Banks (en este caso no nos aparecen)
- Si deseamos visualizar cada una de las secciones, solamente damos
clic en cada sección
Abrimos HEC-RAS
Creamos un nuevo proyecto en File -> new Project… con el nombre
Modelación_Hidraulica
Recuerda guardar todo en una sola carpeta
Para trabajar en unidades del sistema internacional, se debe cambiar la
configuración ya que el software viene con el sistema americano
predeterminado
Ya una vez con la configuración lista, procederemos a cargar la geometría:
Edit -> Geometry Data
lo que nos abrirá la siguiente ventana
Posteriormente se abrirá File -> Import Geometry Data -> GIS Format
Y se buscará el archivo con la geometría que se usará
Recuerde usar el sistema internacional, posteriormente daremos click en next
dejando los parámetros por defecto y daremos click finalmente en finished
import data
Una vez cargado el archivo deberá de aparecer todas las secciones las cuales
son modificables (haciendo click en la sección de interés)
Para una modificación general se puede hacer por medio de la opción de
Tables
Para nuestro ejemplo requerimos modificar el coeficiente de rugosidad de
manning:
Tables -> Manning’n or k values
Al seguir esta ruta se abrirá el siguiente cuadro:
Donde n#1 y n#3 son las márgenes y n#2 representa la parte central del
cause. Para asignar un valor se selecciona la columna deseada y en set values
se asigna un valor (para el caso de nuestro ejemplo se usara 0.025)
Y finalizamos con ok, esto deberá dar como resultado que todos los
coeficientes sean el mismo que se ingresó
Para ambas márgenes en este ejemplo se usara un coeficiente de 0.044
Y finalizamos con OK los datos ingresados. Para guardar se hace por medio de
file -> sabe geometry data as… y finalmente culminamos el apartado de las
condiciones de geometría.
En nuestra consola quedará la geometría por lo que podremos cerrar la
pestaña.
Para este ejemplo trabajaremos con un modelo unidereccional por lo que en
steady Flow data modificaremos en PD 1 (periodo de retorno en 100)
Para ello iremos a options -> Edit profile names y cambiaremos PF 1 por TR
100
Y el valor que le asignaremos al caudal será de 150 mˆ3/ seg
Ya tenemos asignados los caudales. Para las condiciones de contorno debemos
dar click en Reach Boundary Conditions y una vez abierta la pestaña en aguas
arriba como aguas abajo se le asignará el tirante critico
Se finaliza con click en OK y guardamos los datos (de la misma forma que se
guardo la geometría).
En nuestra consola podemos visuakizar que solo falta el Plan , para ello nos
dirigiremos a la pestaña de Run -> steady Flow análisis
Y una vez abierto, procederemos a guardarlo
Y en nombre usaremos TR100
Cambiamos la opción a mixted y activamos la opción Floodplain Mapping
Y por ultimo en options -> Critical Depth output option y activamos la siguiente
opción
Y luego damos click en Compute
