Análisis de transito de avenidas

2020

ANALISIS DE
TRANSITO DE
AVENIDAS

CATEDRATICO: ING. MANUEL ENRIQUE GARCIANARANJO

BUSTOS

INTEGRANTES:
- CANCHANYA TAIPE BRAYAM (HORARIOPRACTICA 9086)
- RIVEROS QUISPE JULISSA (HORARIOPRACTICA 9085)
- ROSALES CONDORI YANELI (HORARIOPRACTICA 9085)
- VILLAGARAY CHUQUILLANQUI SERGIO (HORARIOPRACTICA 9085)

2020-Huancayo
 Análisis de transito de avenidas

INDICE
1. OBJETIVO: ................................................................................................................. 1
2. DEFINICIÓN: .............................................................................................................. 1
3. MARCO TEÓRICO: ..................................................................................................... 1
3.1 DEFINICIÓN DE UN HIDROGRAMA: ................................................................... 1
3.2 PASO DE UNA ONDA A TRÁNSITO: .................................................................... 1
3.3 MÉTODOS: ......................................................................................................... 2
3.3.1 EL MÉTODO HIDROLÓGICO ........................................................................ 2
3.3.2 EL MÉTODO HIDRÁULICO ........................................................................... 7
3.4 SITUACIONES EN LAS QUE SE PRESENTA: .......................................................... 8
3.4.1 TRÁNSITO DE AVENIDA EN RÍOS................................................................. 8
3.4.2 TRÁNSITO DE AVENIDAS EN EMBALSES ..................................................... 8
3.4.3 TRÁNSITO EN EMBALSES CONTROLADOS ................................................ 11
3.4.4 TRÁNSITO EN CAUCES NATURALES .......................................................... 11
4. MODELAMIENTO EN HEC-HMS .............................................................................. 13
5. EJEMPLO DEL METODO MUSKINGUM EN MATHLAB ............................................. 14
6. VERIFICACION DEL CALCULO DEL CAUDAL DE DISEÑO CON EL METODO DEL
HIDROGRAMA TRIANGULAR .......................................................................................... 15
7. CONCLUSIONES :..................................................................................................... 19
8. RECOMENDACIONES : ............................................................................................ 19
 Análisis de transito de avenidas

ANÁLISIS DE TRÁNSITO DE AVENIDAS

1. OBJETIVO:

Comprender cómo evoluciona un tránsito de avenidas a medida que discurre a lo largo

de un cauce o a través de un depósito o embalse.

2. DEFINICIÓN:

Se conoce también con los nombres de Tránsito de hidrogramas y en inglés como:

Hydrograph routing, Flood routing, Flow routing.
Según Velarde (2008) menciona que el tránsito de avenidas es un procedimiento
matemático para predecir el cambio en magnitud, velocidad y forma de una onda de
flujo en función del tiempo (Hidrograma de Avenida), en uno o más puntos a lo largo de
un curso de agua (Cauce o canal).
Aparicio Mijares (1987) afirma que el sitio donde se miden los escurrimientos o donde
se encuentra una presa para control de inundaciones se localiza varios kilómetros aguas
arriba del punto donde las avenidas pueden causar daos, debido a las condiciones
topográficas y geológicas que deben existir para construir una presa o las que debe
reunir el sitio para instalar una estación hidrométrica.

3. MARCO TEÓRICO:

3.1 DEFINICIÓN DE UN HIDROGRAMA:

Un hidrograma de crecida refleja en realidad el movimiento de una onda al pasar por

una estación. Es necesario tener presente que conforme la onda se mueve hacia aguas
abajo su forma cambia. Estos cambios en la onda se deben a la adición de agua de los
tributarios y a que las velocidades en los diversos puntos a lo largo de la onda no son las
mismas.

3.2 PASO DE UNA ONDA A TRÁNSITO:

Las ondas de avenidas se forman por aumento no uniforme del caudal del río por efecto
de una tormenta importante.
Ilustración 1: Paso de una onda a Tránsito

1
 Análisis de transito de avenidas
3.3 MÉTODOS:

Para su estudio hay disponibles dos métodos: el método hidráulico y el hidrológico.

Ambos intentan describir los cambios que en el tiempo experimenta la onda de avenida.

3.3.1 EL MÉTODO HIDROLÓGICO

Está basado en hipótesis simplificadoras y consiste básicamente en plantear la

ecuación de continuidad en un tramo corto de la corriente.
Se basan en la ecuación de la continuidad, referida a un tramo corto del cauce, puede
escribirse, para un tiempo corto:

Volumen que ingresa - volumen que sale = cambio en el

almacenamiento

𝑉𝑖 −𝑉0 = ∆𝑆
Dividiendo por ∆𝑡:
𝑉𝑖 𝑉0 ∆𝑆
− =
∆𝑡 ∆𝑡 ∆𝑡
∆𝑆
𝐼 ̅ − 𝑂̅ =
∆𝑡
 I= Caudal de entrada medio (durante el tiempo Δt)
 O= Caudal de salida medio (durante el tiempo Δt)
 ΔS= (S2- S1)= Incremento en almacenamiento en el tiempo Δt

a) MÉTODO MUSKINGUM

Entre los métodos hidrológicos, posiblemente el más utilizado en cálculos manuales

por su sencillez es el de Muskingum.

El almacenamiento S en un tramo del cauce puede descomponerse en dos partes:

 Almacenamiento en prisma, que sería proporcional al caudal de salida (O).
𝑆𝑝𝑟𝑖𝑠𝑚𝑎 = 𝐾. 𝑂

 Almacenamiento en cuña, que sería función de la diferencia entre el caudal de

entrada y salida (I-O), ya que cuanto mayor sea esa diferencia, más pronunciada
será la cuña.

𝑆𝑐𝑢ñ𝑎 = 𝐾. 𝑋. (𝐼 −O)

2
 Análisis de transito de avenidas

Sumando las dos expresiones:

𝑆 = 𝐾. [𝑋𝐼 + (𝐼 − 𝑋)O]
𝑆
𝐾=
𝑋𝐼 + (𝐼 − 𝑋)O
Donde:
S= Almacenamiento en el tramo considerado de un cauce
I= Caudal de entrada en ese tramo
O= Caudal de salida de ese tramo
K, X= Constantes para ese tramo de cauce
K= Constante de almacenamiento; tiene dimensiones de tiempo
X=Constante sin dimensiones (para la mayoría de los ríos cae entre 0.1 y 0.3).
Si el almacenamiento está dado en m3/seg x día y los caudales están dados en
m3/seg, K tiene unidades de días.
Si se escribe la ecuación: 𝑆 = 𝐾. [𝑋𝐼 + (𝐼 − 𝑋)O] para S1 y S2
𝑆1 = 𝐾. [𝑋𝐼1 + (𝐼 − 𝑋)𝑂1 ]
𝑆2 = 𝐾. [𝑋𝐼2 + (𝐼 − 𝑋)𝑂2 ]
Reemplazando en:
𝐼1 + 𝐼2 𝑂 1 + 𝑂 2 𝑆1 + 𝑆2
− =
2 2 ∆𝑡
𝑂2 = 𝐶0 𝐼2 + 𝐶1 𝐼1 + 𝐶2 𝑂1
Donde:
𝐾𝑋 − 0.5∆𝑡
𝐶0 =
𝐾 − 𝐾𝑋 + 0.5∆𝑡
𝐾𝑋 + 0.5∆𝑡
𝐶1 =
𝐾 − 𝐾𝑋 + 0.5∆𝑡
𝐾 − 𝐾𝑋 − 0.5∆𝑡
𝐶2 =
𝐾 − 𝐾𝑋 + 0.5∆𝑡
∆t = Es el período de tránsito en las mismas unidades de K

3
 Análisis de transito de avenidas
Por lo tanto:
𝐶0 + 𝐶1 + 𝐶2 = 1

Esta fórmula es útil para la comprobación de los cálculos realizados a mano.

Ilustración 2: Cuadro de cálculo del método de Muskingum

Si existen datos de otras avenidas

 K y X pueden ser estimados haciendo un gráfico de S versus [𝑋𝐼 + (𝐼 − 𝑋)O]
para varios valores de X.
 X es aquel que hace tomar a la curva la forma más cercana a una línea recta.
 K es el valor de la pendiente de la recta
Ilustración 3: Método de Muskingum

b) MÉTODO DEL HIDROGRAMA UNITARIO TRIANGULAR

Mockus (1) desarrolló un hidrograma unitario sintético de forma triangular. De la

geometría del hidrograma unitario, se escribe el gasto pico como:

1
Mockus, Victor. Use of storm and watershed characteristics in syntetic unit hidrograph análisis and application. US.
Soil Conservation Service, 1957.
4
 Análisis de transito de avenidas

0.555 A
qp 
tb
Donde:
A = área de la cuenca en km2
tp = tiempo pico en horas
qp= descarga pico en m3/s/mm.
Del análisis de varios hidrogramas, Mockus concluye que el tiempo base y el tiempo
de pico tp se relacionan mediante la expresión:
tb= 2.67 tp
A su vez, el tiempo de pico se expresa como:

de
tp  t r
2
Donde:
de = duración en exceso
tr = tiempo de retraso, el cual se estima mediante el tiempo de concentración tc.

t r  0.6 t c
o bien con la ecuación:
0.64
 L 
t r  0.005 
 S

Donde L es la longitud del cauce principal en metros, S su pendiente en % y tr el

tempo de retraso en horas.
La duración en exceso con la que se tiene mayor caudal pico, a falta de mejores
datos, se puede calcular aproximadamente como:

Para cuencas grandes o como de=tc para cuencas pequeñas.

El caudal máximo se determina tomando en cuenta la precipitación efectiva Pe.

Qmax = qp x Pe
5
 Análisis de transito de avenidas
Pe puede ser calculada tomando en cuenta los números de escurrimiento
propuesto por el U.S. Soil Conservation Service:

2
 508 
 P  N  5.08
Pe 
2032
P  20.32
N
Donde:
Pe : Precipitación efectiva (cm)
N : Número de curva
P : Altura de lluvia (cm)
La determinación del Número de Curva (N), se efectuó tomando en consideración
la información recabada de la visita de campo en aspectos referentes a las
condiciones del suelo y el uso de estos, relacionándolos con los factores
meteorológicos locales.
La subcuenca implicada en la evaluación, en términos generales, se pueden
clasificar como de suelo tipo C, esto se deduce según el Cuadro siguiente:

Para obtener el valor de N de la subcuenca de captación, se debe tener en cuenta

la descripción y tipo de cobertura, el siguiente cuadro muestra los diferentes valores
de número de curva de escorrentía para las diferentes coberturas:
Numero de curva para grupos
Condición
Descripción y tipo de cobertura de suelos hidrológicos
Hidrológica
A B C D
Pastos, forraje para pastoreo Mala 68 79 86 89
Regular 49 69 79 84
Buena 39 61 74 80
Prados continuos, protegidos de
30 58 71 78
pastoreo, y generalmente segado ----
para heno
Mala 48 67 77 83
6
 Análisis de transito de avenidas
Maleza mezclada con pasto de Regular 35 56 70 77
semilla, con la maleza como Buena 30 48 65 73
principal elemento
Combinación de bosques y pastos Mala 57 73 82 86
(huertas o granjas con árboles) Regular 43 65 76 82
Buena 32 58 72 79
Bosques Mala 45 66 77 83
Regular 36 60 73 79
Buena 30 55 70 77
Predios de granjas, construcciones, 59 74 82 86
---
veredas, caminos y lotes
circundantes
Fuente: U.S. Soil Conservation Service

3.3.2 EL MÉTODO HIDRÁULICO

Además de la ecuación de la continuidad, utilizan las ecuaciones del movimiento del

fluido, de modo que para causes o canales en régimen no permanente se utilizan
ecuaciones diferenciales.
De análisis es bastante complejo por cuanto las condiciones naturales a que se
aplica, también lo son: el flujo es no permanente, la sección transversal es no
uniforme, la rugosidad es variable, etc.

a) MÉTODO MUSKINGUM-CUNGE
Cumge combino métodos hidráulicos con simplicidad del método de Muskingum.
Calcula las dos constantes utilizadas en el método de Muskingum, K y X, mediante
parámetros hidráulicos del cauces.
𝐾 = ∆𝑥/𝑐
1 𝑄
𝑋 = (1 − )
2 𝐵𝑆0 𝑐∆𝑥
Donde:
∆𝑥 =Longitud del tramo del cauce considerado
𝑐 = Celeridad=velocidad media*m
𝑚 = Aproximadamente 5/3 para cauces naturales amplios
𝑆0 = Pendiente media del cauce (adimensional)
𝑄 = Caudal
𝐵 = Anchura del cauce
Tener en cuenta:

7
 Análisis de transito de avenidas
La correcta aplicación de este método requiere elegir correctamente el ∆𝑡 y el ∆𝑥.
Para ello se dividirá el tramo estudiado en subtramos, de modo que el caudal de
salida de uno de ellos será el caudal de entrada del siguiente.

3.4 SITUACIONES EN LAS QUE SE PRESENTA:

3.4.1 TRÁNSITO DE AVENIDA EN RÍOS

Supongamos que en el extremo de un canal seco arrojamos un volumen de agua,

entonces el hidrógrama se presentara de la siguiente forma:

Arrojamos un Hidrograma El mismo volumen pasa por

caudal de pico alto y B y C con el hidrograma
agua corta duración cada vez más aplanado

Canal seco

Suponiendo que no hay pérdida de volumen por infiltración o evaporación, o sea

que el área bajo los tres hidrogramas es igual. Entonces se podrá obtener el
hidrograma en el punto C a partir del hidrograma en el punto A.

3.4.2 TRÁNSITO DE AVENIDAS EN EMBALSES

Un aumento en el caudal de entrada producirá también un aumento en el caudal de

salida, pero amortiguado por el deposito.

8
 En el caudal de Análisis de transito de avenidas
entrada se producirá
un hidrograma
elevado

En la salida se
producirá un
hidrograma más
aplanado

Se supone aquí, por simplicidad, que el embalse es no controlado, es decir sin

compuertas de descarga, de manera que la descarga se efectúa por un vertedero
de desborde o aliviadero. En general, el almacenamiento en el embalse hace que se
modifique más marcadamente la forma de onda que en un cauce natural de
longitud equivalente.
∆𝑆
De 𝐼 ̅ − 𝑂̅ = se puede escribir:
∆𝑡

𝐼1 + 𝐼2 𝑂 1 + 𝑂 2 𝑆1 + 𝑆2
− =
2 2 ∆𝑡
𝐼1 + 𝐼2 𝑂 1 + 𝑂 2 𝑆1 + 𝑆2
− 𝑂1 = + − 𝑂1
2 2 ∆𝑡
𝐼1 + 𝐼2 𝑂 2 𝑆2 𝑂 1 𝑆1
− 𝑂1 = ( + ) − ( + )
2 2 ∆𝑡 2 ∆𝑡
𝐼1 + 𝐼2
− 𝑂1 = 𝑁2 − 𝑁1
2
𝑂 𝑆
𝑆𝑖𝑒𝑛𝑑𝑜: 𝑁 = +
2 ∆𝑡
El valor ∆t se llama período de tránsito, y su valor se fija de antemano en la forma
que se indica después. Como la descarga tiene lugar por un aliviadero, 0 es función
de la carga H sobre la cresta del aliviadero; pero una carga H define una cota de la
superficie libre, por lo que la curva de descarga C-O es conocida.
Por otro lado, un nivel de agua define un valor del almacenamiento, es decir la curva
C-S es también conocida. De donde se concluya que la curva N-O se puede siempre
construir para un embalse (basta ve que cada par de valores 0, S están relacionados;
en efecto a un valor de 0 corresponde un valor de C, a ese valor de C corresponde
un S, con lo que el par de valores 0, S definen un único valor de N).
𝑂 𝑆
𝑃𝑜𝑟 𝑙𝑜 𝑡𝑎𝑛𝑡𝑜 𝑑𝑒: 𝑁 = +
2 ∆𝑡
Ilustración 4: Obtención de la gráfica N-O

9
 Análisis de transito de avenidas

Procedimiento de cálculo para el embalse:

 Emplear el siguiente esquema
Ilustración 5: Tránsito de avenida

Estando el embalse lleno llega el flujo base de la corriente, 𝐼0 este valor inicial del
caudal de llegada corresponden los valores iniciales de nivel de agua,
embalsamiento, carga y caudal de salida. En estas circunstancias se supone llega la
avenida representada por el hidrograma I-t.
 Pasos de cálculo

 En este hidrograma se adopta un valor de ∆t

 Se leen las escorrentías directas I
 El otro valor conocido es la descarga inicial
 Para llenar la tabla está indicado con flechas en la tabla misma.
Ilustración 6: Cuadro de cálculo para el embalse

10
 Análisis de transito de avenidas

3.4.3 TRÁNSITO EN EMBALSES CONTROLADOS

Para un embalse con compuertas en el vertedero el gasto de salida depende del

número de compuertas que estén abiertas. Entonces la curva N-O debe ser
sustituida por una familia de curvas, en las que el número de 1 a 3 indica el número
de compuertas abiertas.
Ilustración 7: Embalses controlados

Procedimiento de cálculo para embalses controlados:

La operación de análisis es similar al que se realiza en el tránsito de embalse, salvo
que el número de compuertas debe tabularse.
Si no hay cambios en la abertura de las compuertas durante el tiempo de estudio el
procedimiento es idéntico a la tabla usada en el tránsito de embalse, ya que todos
los valores se leen en la curva que representa la abertura constante de las
compuertas.

3.4.4 TRÁNSITO EN CAUCES NATURALES

Ahora el almacenamiento no es función única del caudal de salida. El

almacenamiento S viene a ser el volumen de agua en el cauce en cualquier instante.

11
 Análisis de transito de avenidas
∆𝑆
Como la ecuación 𝐼 ̅ − 𝑂̅ = ∆𝑡 involucra únicamente a ∆S, no se necesita conocer los
valores absolutos de almacenamiento.
Al pasar una onda de avenida por un tramo, hay cierto aumento en el
almacenamiento antes de que haya cualquier incremento en el gasto de salida.
Después que la cresta de la onda ha entrado al tramo, el almacenamiento puede
empezar a disminuir aunque el gasto de salida esté todavía aumentado.
Ilustración 8: Tránsito en cauces naturales

Procedimiento de cálculo para el embalse:

 Los hidrogramas de entrada y de salida para el tramo en estudio se dividen en

intervalos de tiempo pequeños
 Se determina para cada período valores promedios de 1 y de O
∆𝑆
 Se calculan los correspondientes valores de ∆S con la ecuación 𝐼 ̅ − 𝑂̅ =
∆𝑡

Ilustración 9: Cálculos de almacenamiento

12
 Análisis de transito de avenidas

 Con los valores de S calculados, se representan en una gráfica versus los gastos
simultáneos de salida.
Ilustración 10: Curva de almacenamiento

 Usualmente aparece que el almacenamiento es ligeramente mayor durante el

tiempo de ascenso del nivel, que durante el tiempo de descenso.

4. MODELAMIENTO EN HEC-HMS

El HEC-HMS es un programa de simulación hidrológica tipo evento, lineal y

semidistribuido, desarrollado para estimar los hidrogramas de salida en una
cuenca o varias subcuencas (caudales máximos y tiempos pico) a partir de
condiciones extremas de lluvia, aplicando para ello algunos de los métodos de
cálculo de histogramas de diseño, pérdidas por infiltración, flujo base y
conversión en escorrentía directa que han alcanzado cierta popularidad en los
EE. UU y por extensión en nuestro país.

El modelo HEC-HMS requiere la siguiente información:

13
 Análisis de transito de avenidas
 Información acerca de la precipitación histórica o de diseño.
 Información acerca de las características del suelo.
 Información morfométrica de las subcuencas.
 Información hidrológica del proceso de transformación de lluvia en escorrentía.
 Información hidráulica de los tramos del cauce principal.

5. EJEMPLO DEL METODO MUSKINGUM EN MATHLAB

Resultado del método Muskingum con el programa MATHLAB:

Ilustración 11 Datos iniciales para el programa

Ilustración 12 Datos de Q de entrada y resultados de Q de salida

14
 Análisis de transito de avenidas
Ilustración 13 Gráfico final de hidrogramas

Resultado con el programa HEC-HMS:

Ilustración 14 Modelamiento en HEC-HMS

Interpretación:

 El uso del mathlab nos da resultados más exactos que los valores calculados
en el HEC-HMS.

6. VERIFICACION DEL CALCULO DEL CAUDAL DE DISEÑO CON EL METODO DEL

HIDROGRAMA TRIANGULAR

Datos usados del expediente técnico: “RECONSTRUCCIÓN DE LA REPRESA EN LA LAGUNA DE

YANACOCHA DE LA LOCALIDAD DE VILLA DE PASCO, DISTRITO DE FUNDICION DE TINYAHUARCO,
PROVINCIA DE PASCO – PASCO”

15
 Análisis de transito de avenidas

Cálculos realizados en el Expediente Técnico por dos métodos:

1° Método Hidrograma Triangular

Ilustración 15 Cuadro N°28 del Expediente Técnico utilizando el Método Hidrograma Triangular

16
 Análisis de transito de avenidas

2° Modelamiento HEC-HMS:

Ilustración 16 Cuadro N°29 del Expediente Técnico utilizando HEC-HMS

Conclusión por el proyectista en el Expediente Técnico:

Finalmente, el caudal de diseño de máximas avenidas para un periodo de retorno de 500 años será el mayor valor estimado por los dos métodos
mencionados anteriormente, esto para dar un mayor margen de seguridad a la estructura a proyectar.
El valor de 12.22 m3/s (Método de Hidrograma Triangular) será usado para el diseño y dimensionamiento del vertedero de demasías.

17
 Análisis de transito de avenidas

Cálculo propio realizado en Excel:

18
 Análisis de transito de avenidas

Interpretación:

Los calculados realizados en Excel nos da un valor de 12.11 m3/s para el

caudal de diseño de máximas avenidas para un periodo de retorno de 500
años.

Comparación de los resultados del expediente por el método de Hidrograma

Triangular, HEC-HMS y el cálculo realizado en Excel:

 Al igual que la comparación del método Muskingum en el mathlab y el modelado

en HEC-HMS; en esta comparación del expediente nos muestra que los
resultados obtenidos por el método del Hidrograma Triangular nos da resultados
más exactos que los valores calculados en HEC-HMS y Excel.

7. CONCLUSIONES :

 El tránsito de avenidas es útil para el diseño y/o verificación del vertedero de

excedencia y para el cálculo de la altura de la presa.

 Los métodos tradicionales para el tránsito de avenidas producen errores de

aproximación que son introducidos durante la consulta de curvas que
depende de la interpretación personal.

 La ecuación principal de Muskingum tiene una solución y puede resolverse

mediante métodos numéricos o hasta algebraicos.

 El método de Hidrograma Triangular sirve para cuencas de superficie mayor

a 10 Km2, calculando la descarga máxima estimada.

8. RECOMENDACIONES :

 Los métodos directos (Muskingum y Hidrograma Triangular) son una mejor

alternativa porque permiten resultados precisos, total automatización,
mejora los procedimientos y requiere mínima carga de trabajo.

 Cuando se utilice programas ya creados en Mathlab es recomendable usar la

misma versión cuando se creó dicho programa ya que no cuenta con
actualizaciones para las futuras versiones de Mathlab.

 Comprobar los resultados del HEC-HMS con algún otro método directo para
llegar al resultado máximo ya que no es un programa totalmente
recomendado por dar resultados mínimos.

19