UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
MODELACION RIO SEK’E
1. ANTECEDENTES
La Cuenca Katari, es una de las más habitadas de Bolivia. La misma sirve como principal fuente de abastecimiento de
agua para el Municipio cíe El Alto, tanto para consumo como para industrias, alimentación, curtiembre, servicios,
comercio, agricultura y ganadería. El Municipio de El Alto se encuentra ubicado dentro las subcuencas Achicala,
Sallani, Seque y Seco, Tuni Condoriri y Boopy, las primeras cuatro pertenecen a la cuenca Katari y la última a la
cuenca Beni.
Para el estudio del Modelaje del río Seke y Revisión de Diseño de Elementos de los Proyectos: Canalización río Seke
Tramo II y Canalización río Seke Tramo Norte del Drenaje Pluvial de la Ciudad del Alto, se tiene contemplado la
necesidad del estudio de modelación hidráulica sobre un tramo del río Seke para fines de conocimiento de
parámetros y perfiles hidráulicos y en este caso hacemos uso del modelado del rio en el programa HEC RAS, los
mismos pudiendo constituirse en elementos base para el diseño de obras de canalización de dos tramos importantes
de éste río e identificados como: Tramo II Tramo Norte.
El primer tramo tiene una longitud aproximada de 1400 m, el segundo 1800 m y el tercero 565 m; entre los tres
tramos existe un tramo ya canalizado mediante muros longitudinales. Entre tanto el criterio asumido para la
modelización es la consideración de los tres tramos en conjunto, para el cual debe definirse parámetros hidráulicos
de flujo y las características de los perfiles hidráulicos de flujo para poder definir las características de este. Para ello,
una de las informaciones de entrada a considerarse es la crecida, la misma es definida sobre información física
derivada de la imagen satelital, información meteorológica en estaciones cercanas. La otra información es referida a
la geometría del cauce del río Seke, particular del tramo indicado y las propiedades del material del lecho.
2. INTRODUCCION
La modelación del río Seke y Revisión de Diseño de Elementos de los Proyectos: Canalización río Seke Tramo II y
Canalización río Seke Tramo Norte del Drenaje Pluvial de la Ciudad del Alto, se tiene contemplado la necesidad del estudio
de modelación hidráulica sobre un tramo del río Seke para fines de conocimiento de parámetros y perfiles hidráulicos, los
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
mismos pudiendo constituirse en elementos base para el diseño de obras de canalización de dos tramos importantes de
éste río e identificados como: Tramo II Tramo Norte.
2.1. OBJETIVOS
El principal objetivo de este proyecto es la modelación del rio seke mediante la aplicación de la herramienta HEC RAS que
realiza modelos digitales, en este caso se realizara la modelación del rio seke con un caudal determinado dado como dato
académico. Además de poder definir los parámetros hidráulicos de flujo
Además de poder realizar la modelación se puede hacer la comparación de métodos de modelación, referente a los
programas ARC GIS y AUTOCAD CIVIL 3D.
2.1.1. OBJETIVOS ESPECIFICOS
Definir la modelación digital del rio seke para modelar sus perfiles hidráulicos, mediante la morfología y condiciones
hidráulicas que posee este de este y también debido a la obra de canalización.
Cambiar la sección o incorporar diques en las secciones donde se presenta un flujo de régimen supercrítico
(Número de Froude > 1), para reducir el número de Froude y así tener flujos subcríticos.
Hallar las secciones del lecho del río cada 100m en CIVIL 3d y exportar estas secciones al modelo HEC-RAS para un
posterior análisis en este programa.
2.3. CUENCA DE ESTUDIO
2.3.1. MORFOLOGIA DE CUENCA DE RIO SEKE
El cauce principal nace en las faldas del nevado Huayna Potosí, donde se presentan fuertes pendientes, incluso mayores al
50%; en cambio en la parte baja la pendiente no superan el 10%. Parte del cauce principal atraviesa la ciudad de El Alto,
donde se encuentran los dos tramos de interés para la canalización.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
2.3.2. MAPA DE USO DE SUELOS
Mediante datos proporcionados por la página de la alcaldía podemos extraer el mapa de uso de suelos ya que este posee
imágenes satelitales del uso del suelo.
3. MARCO TEORICO
3.1. METODO DE APLICACIÓN
3.1.1. FLUJO GRADUALMENTE VARIADO
Para la aplicación al modelo a determinar, nos basamos en la ecuación del flujo
establecido por Bernoulli analizando el flujo en cuanto a sus energías causadas por altura piezometrica y altura de
velocidad.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
El flujo gradualmente variado, es un flujo permanente cuya profundidad varía de manera gradual a lo largo del canal. Se
tendrán en cuenta las siguientes hipótesis:
La pérdida de altura en una sección es igual que la de un flujo uniforme con las mismas características de velocidad
y radio hidráulico.
La pendiente del canal es pequeña (<10%). Esto quiere decir que la profundidad del flujo puede medirse
verticalmente o perpendicularmente al fondo del canal y no se requiere hacer corrección por presión ni por arrastre
del aire.
El canal es prismático.
Los coeficientes de distribución de la velocidad y el de rugosidad son constantes en el tramo considerado.
El método de aplicación de Bernoulli es basado en energías y se resuelve por método estándar por pasos y método directo.
Para poder analizar el flujo que transcurre en el rio es necesario establecer un método de cálculo para aplicarlo
considerando las condiciones del flujo necesarias como las geométricas. En este caso consideramos la modelación numérica
mediante el programa HEC RAS.
3.2. TIPOS DE PERFILES DE FLUJO
Cualquiera que sea el valor del caudal, la pendiente y la forma de la sección transversal, éstas desigualdades dividen el canal
en tres zonas en las cuales se puede desarrollar el perfil del flujo. Estas zonas son:
A. Zona 1. Cuando el perfil se encuentra por arriba de la línea superior
B. Zona 2. Cuando el perfil se encuentra entre las dos líneas
C. Zona 3. Cuando el perfil se encuentra entre la línea inferior y la plantilla del canal.
Clasificación por zonas para una pendiente subcrítica.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
La clasificación de los perfiles de flujo variado está basada en la pendiente del canal y en la zona en que se localiza el perfil,
como lo muestra la figura 1. En el caso de pendiente positiva (el fondo del canal desciende en la dirección del flujo), se
puede establecer un flujo uniforme de tirante Yn, por lo cual dicha pendiente podría ser:
Suave yn > yc perfil tipo M
Crítica yn = yc perfil tipo C
Pronunciada yn < yc perfil S
Horizontal yn tiende a ∞ , perfil H
Adversa yc < la horizontal, perfil A
3.3. NUMERO DE FROUDE
El número de Froude (Fr) es un número adimensional que relaciona el efecto de las fuerzas de inercia y las fuerzas
de gravedad que actúan sobre un fluido. Debe su nombre al ingeniero hidrodinámico y arquitecto naval inglés
William Froude (1810 - 1879). De esta forma, el número de Froude se puede escribir como:
L - parámetro de longitud [m]
v - parámetro de velocidad [m/s]
g - aceleración de la gravedad [m/s²]
El flujo se clasifica como:
Fr<1, Flujo subcrítico o tranquilo, tiene una velocidad relativa baja y la profundidad es relativamente grande,
prevalece la energía potencial. Corresponde a un régimen de llanura.
Fr=1, Flujo crítico, es un estado teórico en corrientes naturales y representa el punto de transición entre los
regímenes subcrítico y supercrítico.
Fr>1, Flujo supercrítico o rápido, tiene una velocidad relativamente alta y poca profundidad prevalece la energía
cinética. Propios de cauces de gran pendiente o ríos de montaña
3.4. METODOLOGIA DE APLICACION
3.1.1. CARACTERISTICAS DEL MODELO HEC RAS
Una de los objetivos del programa HEC RAS es la modelación en base a parámetros topográficos del rio, como se hizo en
este caso. Y uno de los motivos fue determinar las características en cuanto a crecidas del rio con un caudal dado en este
caso 3.5 m3/s.
A partir de la entrada de datos: geometría del cauce, el caudal de crecida, granulometría; en calidad de datos de salida del
cómputo es posible la generación: del perfil longitudinal de profundidad de (o nivel de agua), el perfil de pendiente de energía,
perfil de condición de flujo. El cálculo de tales perfiles es realizado mediante una solución iterativa, por medio del método del
paso estándar. Por otro, el programa al calcular el perfil de profundidad y al compararse con la cota de las márgenes, es
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
posible visibilizar el desborde. Por otro lado, es posible introducir secciones mejoradas para evitar el fenómeno anterior
igualmente; además poder avaluar el resto de los perfiles, de los cuales depende la estabilidad del curso canalizado.
3.1.2. INFORMACIÓN DE CALIBRACION
Tomaremos como datos de entrada
LONGITUD DE RIO= 3765.04 [M]
COEFICIENTE DE MANNING n=0.021 determinado de experimentación en informe de docente
UBICACIÓN DE BORDES DE LA SECCION:
Para este punto se consideró un ancho efectivo de sección de 50 m a cada lado especialmente para poder visualizar la
sección topográfica, ya que se utilizó el meto de exportación de DIBUJO CAD CIVIL 3D A HEC RAS, este proceso no permite
que la sección no corte al borde del canal marcado, en caso contrario no lo reconocería como sección.
INFORMACION DE LAS ESTRUCTURAS DE PASO. Para los tramos de pendientes del perfil de la sección
Los tramos de pendiente se tomaron en cuenta como si fuesen en tres tramos de pendientes uniformes. Estos valores de
pendiente es recomendable generarlo de plano de pendientes de la topografía para tener los valores exactos de valores.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
cota superior − cota inferior
So =
longitud de tramo
Para los tres tramos considerados: 0+000-1+400 So=0.02
1+400-3+200 So=0.021
3+200-3+765.04 So=0.024
La división de tramos de estudio del rio seke son:
TRAMO UBICACION LONGITUDES SECCIONES CONTENIDAS
1 TRAMO SUP 1400 100 A 1400
2 TRAMO MEDIO 1800 1400 A 3200
3 TRAMO INFERIOR 565 3200 A 3700
3.1.3. GEOMETRÍA DE CAUCE
Determinada en ARC GIS nos proporciona la topografía y está dividida cada 100 metros teniendo 37 secciones del rio
analizado. Este procedimiento permitió modelar el perfil en la modelación del HEC RAS de las 37 secciones modeladas con
más detalle de la topografía logrando la captura de más puntos tomados en cuenta
RIOSEKE Plan: Plan 02 09/05/2018
.021 .021 .021
4097.0 Legend
EG PF 1
4096.5
WS PF 1
4096.0 C rit PF 1
Elevation (m)
Ground
4095.5
Bank Sta
4095.0
4094.5
4094.0
0 10 20 30 40 50 60
Station (m)
SECCION DE TOPOGRAFIA
La geometría del cauce fue determinada mediante ayuda del programa ARC GIS, con la herramienta GEO RAS.
La herramienta GEO RAS es una muy útil herramienta para la determinación del cauce del rio SEKE ya que simplifica el
trabajo de estar manipulando las secciones de manera individual, a una manera más automatizada; ayudando a determinar
la sección con más precisión.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
En caso de que se hubiese tomado las secciones con puntos transcritos la aproximación de la sección sería más
esquematizada.
Realizamos el siguiente proceso de datos
ESTOS RESULTADOS ESTAN EN BASE A LA VERSION DE ARC GIS [Link]
LA VERSION DE GEO RAS ES 10.3.1.
IMPORTAMOS LA IMAGEN DEM PROPORCIONADA POR EL DOCENTE EN AUTOCAD EXPORTANDOLA A FORMATO
DEM.
GENERAMOS CURVAS DE NIVEL EN ARC GIS CON EL COMANDO ARCTOLLBOX.
IMPORTAMOS LA UBICACIÓN ESQUEMATIZADA DEL RIO SEK’E.
UBICACMOS EN ANCHO DEL CANAL DEL RIO.
INICIAMOS LA APLICACIÓN GEO RAS.
DEFINIMOS EL RIO.
DEMARCAMOS EL BORDE DERECHO LUEGO EL IZQUIERDO; ESTOS PASOS DEFINEN DONDE INICIA EL RIO Y DONDE
ES AGUAS ABAJO.
TRAZAMOS LAS SECCIONES CADA 100 METROS.
NOS ASEGURAMOS QUE LAS SECCIONES NO SE INTERSECTEN NI ESTEN MENOS QUE EL ANCHO DEL CANAL DEL RIO
NI EL CANAL DE POSIBLES EXCEDENCIAS.
EXPORTAMOS EL FORMATO A HEC RAS.
Proceso de ejecución en HEC RAS:
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
Para proceder al modelado de la sección analizada debemos proseguir de la siguiente manera:
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
También es posible realizar el mismo procedimiento de civil 3d a HEC RAS
Con las herramientas ya mencionadas y el cálculo ya establecido procedemos al modelado de la sección del canal del RIO
SEK’E.
3.1.4. CONDICIONES HIDRÁULICAS CONSIDERADAS
Las condiciones hidráulicas que se consideraron fueron las de pendiente
En condiciones normales de flujo se tomó pendientes aguas arriba como debajo de valor 0.024 considerada en tramos las
cuales se introdujeron de la siguiente manera:
0+000-1+400 So=0.02
1+400-3+200 So=0.021
3+200-3+765.04 So=0.024
En cuanto a los coeficientes de contracción y expansión se tomó en cuenta:
Que son determinador por defecto en el programa.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
4. RESULTADOS DE LA MODELACION
Después de haber realizado el cálculo con todos los parámetros necesarios para el cálculo procedemos a analizar los
resultados obtenidos por el programa.
4.1. Perfiles de flujo
Los perfiles de flujo son parámetros por el cual se observan los niveles de la superficie del flujo en todas las secciones. En este
perfil se observa la modelación realizada por el programa que considero las condiciones morfológicas y condiciones
hidráulicas ya mencionadas.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
4.2. PARAMETROS HIDRAULICOS
A continuación, en la figura. Se presentan los parámetros hidráulicos en una tabla, con sus respectivas unidades en el
Sistema Internacional (SI).
5. ANALISIS DE RESULTADOS
Para el uso de los programas es necesario tener cuidado con las unidades ya sea en HEC RAS y el ARC GIS ya que podría
causar error en las secciones de cálculo y modificar el caudal de diseño.
Para las secciones se observa que para un caudal de 3.5 m3/s las secciones no sobrepasan el nivel de aguas critico en la
mayoría de las secciones.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
Distancia Sección Profundidad de flujo Posibilidad de desborde
(m) No (m)
300-400 3-4 0.57 POSIBLE DESBORDE
En esta sección el caudal sobrepasa el nivel crítico de aguas por lo que es necesario la aplicación de una obra auxiliar.
Para poder hacer cumplir la condición de flujo suscritico se procedió a considerar que el flujo tendrá gradas de caída en la
sección del rio los cuales ayudan en la caída del flujo y este disminuye la energía de caída produciendo un menor flujo y este
será suscritico analizando el número de froude.
Se modificaron las secciones debido a que las velocidades de flujo en los tramos eran muy elevadas, estas velocidades
eran erosivas, como se observa en la figura anterior que pasaba del nivel de la sección produciendo erosion, para la
progresiva 3+600 se tiene una velocidad de 2.03 m/s, lo cual el flujo en esa sección pertenece al régimen supercrítico,
realizada ya la modificación en dicha sección se puede ver en la siguiente figura que la velocidad disminuyó al igual que el
número de Froude, teniendo como valores 1.06 m/s y para el numero de froude 0.98 respectivamente. Estos valores nos
indican que ahora en la sección se tiene un flujo de régimen subcrítico.
Para la sección en la progresiva 3+400 también se realizó el mismo criterio, modificando la sección, en la Fig.5 se puede
ver que inicialmente se tní un Núm de Froude de 1,01 (supercrítico), al cambiar la sección el Núm de Froude disminuyó a
un valor de 0,99 (subcrítico); al igual que las velocidades estas también disminuyeron de 1.74 (m/s) a 1,03 (m/s).
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
La modificación de la sección de la que se describió anteriormente se produjo una alteración en su base, pero no se hizo
este análisis en todas las secciones, y este método sirvió para que las condiciones hidráulicas que posee como el número
de froude cambien. Esta modificación llevo a cabo que las secciones analizadas se incorporara gradas de disipación de
energía que ayuden a disipar la energía y se procedió de la manera mencionada anteriormente incorporándolas cada 50
m entre secciones donde sea necesario. Y se llegó a la siguiente tabla después de haber realizado el procedimiento
mecionado:
Resultando la incorporación de 18 gradas de disipación de energía.
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
Procedimiento de colocación degradas de disipación de energía
PERFIL LONGITUDINAL
6. CONCLUSIONES
Se determinaron los parámetros y perfiles hidráulicos del río y como se puede observar en la Figura la mayor
parte de las secciones dieron como resultados números de Froude por encima de la unidad, por lo que se requirió
de un análisis adicional ya mencionado en análisis de resultados.
Como se mencionó anteriormente respecto a los números de Froude superior a la unidad, se modificó la sección
o se incorporaron gradas de disipación de energia para cambiar el régimen del flujo y poder trabajar con flujos
subcríticos en lugar de flujos supercríticos. A continuación se mostrará los parámetros hidráulicos finales:
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
PAGINA | DE
UNIVERSIDAD MAYOR DE SAN ANDRÉS PROYECTO 1 MODELACION RIO SEKE
FACULTAD DE INGENIERÍA DOC ING. RAMIRO PILLCO ZOLA
INGENIERÍA CIVIL
ESTUDIANTE: ALVAREZ SEGALES KEVIN DOUGLAS
LABORATORIO DE HIDRAULICA
REYNALDO QUISPE LECONA
7. LITERATURA
HIDRAULICA I SUYO LARUTA, HERNAN; ENERGIA ESPECIFICA
MECANICA DE FLUIDOS MC GRAW HILL
APUNTES DE CATEDRA HERBAS CAMACHO CARLOS
HIDRAULICA SOTELO AVILA
PAGINA | DE