"AÑO DE LA LUCHA CONTRA LA CORRUPCIÓN E IMPUNIDAD”

FACULTAD DE INGENIERÍA

ESCUELA PROFESIONAL DE INGENIERÍA CIVIL

HIDROLOGIA

ALUMNO

 PAUL VELASQUEZ LOZADA

TEMA

ESTUDIO DE TRANSITO DE AVENIDAS DEL RIO PIURA POR EL METODO

DE MUSKINGUM - TRAMO TAMBOGRANDE – PIURA

DOCENTE

HECTOR YAURI QUISPE

PIURA – PERÚ

2019
1. INTRODUCCION
Actualmente, el uso de estos modelos, permiten realizar análisis hidrológicos
para la prevención de inundaciones en áreas aledañas o riveras de los ríos;
además, permiten formular hipótesis que ofrecen cierto grado de confianza
para la toma de decisiones.
El uso de modelos hidrológicos para el estudio de avenidas o crecidas de los
ríos es una herramienta que se ha extendido por todo el mundo, especialmente
en países desarrollados.
El método de Muskingum es un modelo hidrológico, que estudia flujos de
entrada y salida en un punto de control (estación hidrométrica) en un cauce de
un río, y con la ayuda de información hidrológica de caudales y niveles,
previamente registrados, se obtienen los parámetros (K, X) y los coeficientes
(C1, C2 y C3) que son operadores aplicables del método de Muskingum. Este
método, permite obtener una simulación del comportamiento de la crecida en el
tramo de un río, el cual se representa gráficamente mediante una serie
sucesiva de hidrogramas en sentido aguas abajo, en donde la característica
principal, es decir el pico del hidrograma que representa la crecida, sufre un
proceso de descenso por efecto de amortiguamiento producido por la fricción
con las paredes y el fondo del cauce.
La Región Piura, se vio afectada por el fenómeno del Niño Costero, debido a
las inundaciones tanto en el ámbito rural como en el urbano, relacionadas con
la máxima crecida del río Piura, ocurrida el día 27 de marzo del año 2017.
El proyecto tiene como finalidad realizar por medio del método de Muskingum
el estudio del tránsito de avenidas en el rio Piura, puntualmente en las
estaciones hidrométricas de Tambogrande y Puente Sánchez Cerro (Piura), en
el departamento de Piura, como consecuencia de la crecida causada por el
fenómeno del “Niño Costero” (27 Marzo del 2017), tragedia que prácticamente
rebeló dramáticamente el alto grado de vulnerabilidad, que caracteriza a
nuestro país, ante evento naturales que provocan desastres como las
inundaciones

2. ANTECEDENTES
Estudio De Tránsito De Crecientes Con Desarrollo De Una Aplicación
Informática, En El Río

Fonce Afluente Del Río Suárez En El Departamento De Santander.

Tesista: Alonso Homero Contreras

Institución: UNIVERSIDAD DISTRITAL FRANCISCO JOSE DE CALDAS
COLOMIA

Año: 2016

Conclusiones:

 Se desarrolló con éxito una aplicación informática por medio del VBA
que utiliza los datos de niveles horarios anuales, calcula el caudal
correspondiente a dichos niveles y realiza en tránsito de crecientes por
el método de Muskingum para sistemas agregados.
 Esta aplicación informática pretende dar un inicio a la programación de
Tránsito de Creciente, no solamente para tránsitos agregados, sino
también tránsitos distribuidos. Se recomienda realizar un trabajo
multidisciplinar con el área de Ingeniería de Sistemas que ayude a
insertar modelos más complejos que requieran programación más
avanzada.
 Se aplica de manera exitosa el método de tránsito agregado de
crecientes, para el año seleccionado, transitando de manera particular el
caudal generado en un tiempo determinado.
 Después de la correcta lectura e interpretación de los datos
suministrados por el IDEAM, se calculó de manera clara y organizada
los caudales correspondientes al año seleccionado, de los niveles
horarios de cada estación. Determinando el Hidrograma Anual, que
permite realizar el posterior estudio de tránsito de creciente para cada
estación hidrométrica.

Transito hidráulico de avenidas por un cauce natural utilizando el método

de las características

Tesista: Alonso Homero Contreras

Institución: UNIVERSIDAD NACIONAL AUTONOMA DE MEXICO. Año: 2010

Conclusiones:

 Para realizar el tránsito hidráulico se utilizó la solución por el método de

las características, empleando el programa NUCARI (ver Anexo A) que
considera el tránsito de avenidas en un cauce de sección rectangular.
 Los métodos hidráulicos para el tránsito de avenidas en cauces al no
contar con una solución analítica para las ecuaciones de Saint-Venant,
se utilizan procedimientos numéricos para su solución como son los
esquemas de diferencias finitas tanto implícitos como explícitos, el
método de las características (el cual se explica en el capítulo 3) o el
método del elemento finito.
 El método de Muskingum se puede aplicar cuando se tiene un
hidrograma de entrada y también si se tienen gastos de ingreso y
egreso.
 Estos métodos describen en forma precisa el proceso de tránsito tanto
en el tiempo como en el espacio, pero requieren de información
detallada de rugosidades, geometría del cauce, condiciones de frontera
y de técnicas numéricas un poco complejas; los principales factores que
determinan la selección de un método para realizar un tránsito son los
datos disponibles y la información de salida que se requiera.

 En general, existen dos técnicas para el tránsito de avenidas en cauces,

los métodos hidrológicos, o de almacenamiento, y los hidráulicos, o de
proceso; los primeros utilizan la ecuación de continuidad y una relación
entre el almacenamiento y el gasto de salida, además de ser menos
complicados, y los métodos hidráulicos utilizan la ecuación de
conservación de la masa y de conservación de la cantidad de
movimiento (ecuación de Saint-Venant) para describir el flujo de agua en
un cauce.
Método De Muskingum Para El Estudio De Tránsito De Crecientes En El
Río Atrato Entre Las
Estaciones “El Siete Y Gindrama” En El Departamento Del Chocó
Tesista: Jessica Lorena García Castiblanco - Mayra Alejandra Rodríguez
Muñoz.
Institución: UNIVERSIDAD DISTRITAL FRANCISCO JOSE DE CALDAS.
Año: 2016
Conclusiones:
 Los valores calculados para la constantes X - K en las estaciones El
Siete y Gindrama son de X=0.1 K=0.5, lo que hace referencia a cauces
naturales ya que el rango de X está entre 0.1 y 0.3.
 En cuanto a los coeficientes del método C0, C1 y C2, éstos varían, no
mostrando uniformidad en los resultados, de lo que se concluye que
cada crecida se comporta diferente a otra. De igual manera para todos
los casos se cumple que la sumatoria de los coeficientes sea igual a uno
(1) como lo exige el método.

 Inicialmente se había planteado el estudio de crecientes mediante los

tiempos de viaje de ondas de crecida, pero dicho método no fue posible
implementarlo para el desarrollo de este proyecto debido a que el
IDEAM no tiene estudios de ondas de crecida, o velocidad de onda para
las estaciones, por tanto, no es posible emplear el método de
Muskingum.
 Se comprueba que, durante el avance de la onda de creciente, el caudal
de entrada es mayor que el caudal de salida, existiendo almacenamiento
en cuña, durante la recesión, el caudal de salida es mayor que el caudal
de entrada, resultando una cuña negativa.
 3. OBJETIVOS
3.1. OBJETIVO GENERAL

Evaluar el exceso de agua que origino que el centro de la

ciudad de Piura quede inundada el día 27 de marzo del 2017,
para lo cual se va usar el método de Muskingum en las siguientes
estaciones: Puente Sánchez Cerro-Piura y Tambogrande.

3.2. OBJETIVOS ESPECÍFICOS

 Construir hidrogramas de las siguientes estaciones: Puente
Sánchez Cerro-Piura y Tambogrande.
 Hallar las variables K y X que corresponden al cauce entre las
estaciones Puente Sánchez Cerro-Piura y Tambogrande.
 Evaluar el tránsito de la creciente máxima en un tramo aguas
debajo de las estaciones Puente Sánchez Cerro-Piura y
Tambogrande.
4. PROCEDIMIENTO
4.1. DEMARCACION DEL AREA DE ESTUDIO

La cuenca del río Piura tiene una superficie de 12,216 Km2.

Las estaciones se encuentran situadas en el rio Piura que es en la vertiente del

Pacifico sur y su punto de inicio es en la sierra de San Jorge y tiene como
nombre rio Huarmaca, provincia de Morropón y recorre 280 km en el
departamento de Piura, al inicio tiene rumbo Nor

-Este y a mitad de su recorrido cambia a Sur-Este para luego su

desembocadura es en Laguna Ramon-Bajo Piura.

La ubicación geográfica de las estaciones son las siguientes:

- TAMBOGRANDE

ESTACION TAMBOGRANDE

LATITUD 4°55´37" S

LONGITUD 80°20´41" O

DEPARTAMENTO Piura

PROVINCIA Piura

REGION Piura

- PIURA

ESTACION PIURA

LATITUD 5°11´40" S

LONGITUD 80°37´58" O

DEPARTAMENTO Piura

PROVINCIA Piura

REGION Piura
 4.2. MÉTODOS HIDROLÓGICOS
En función del tiempo se calcula y se usa la ecuación de continuidad tomando
un tramo de cauce y por lo tanto se da lo siguiente:

Volumen entrada en un 𝑡𝑖𝑒𝑚𝑝𝑜 ∆𝑡 − 𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑠𝑎𝑙𝑖𝑑𝑜 𝑒𝑛 𝑒𝑠𝑒 ∆𝑡 = ∆

𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖ento

Si la expresión se divide entre ∆𝑡, se tiene:

𝑄𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑄𝑠𝑎𝑙𝑖𝑑𝑎 = ∆𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜⁄∆t (1)

O también:

𝐼 − 𝑂 = ∆𝑆⁄∆𝑡 (2)
𝐼 − 𝑂 = (𝑆2 − 𝑆1) ⁄∆t (3)

Siendo:

𝐼 = 𝑐𝑎𝑢𝑑𝑎𝑙 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑚𝑒𝑑𝑖𝑜 𝑑𝑢𝑟𝑎𝑛𝑡𝑒 ∆𝑡

𝑂 = 𝑐𝑎𝑢𝑑𝑎𝑙 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑚𝑒𝑑𝑖𝑜 𝑑𝑢𝑟𝑎𝑛𝑡𝑒 ∆𝑡

 ∆𝑆 = 𝑆2 − 𝑆1 = 𝑖𝑛𝑐𝑟𝑒𝑚𝑒𝑛𝑡𝑜 𝑑𝑒𝑙 𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜 𝑒𝑛 𝑒𝑙 𝑡𝑖𝑒𝑚𝑝𝑜 ∆t

Para calcular con exactitud los caudales medios de cada ∆ , se debe disponer
de un hidrograma continuo, pero si se conoce solamente un dato de caudal
para cada ∆𝑡 , los caudales medios se pueden evaluar por el promedio de los
caudales de dos ∆𝑡 consecutivos, así:

𝐼1+𝐼2/2 − 𝑂1+𝑂2/2 = 𝑆2−𝑆1/∆t (4)

Una vez conocido el hidrograma de caudales de entrada el problema se reduce

solo en encontrar el hidrograma de caudales de salida.

4.3. MÉTODO DE MUSKINGUM

Método, desarrollado por G. T. McCarthy en 1938, y toma el nombre de

Muskingum del distrito donde fue desarrollado por el Servicio de Conservación
del Distrito de Muskingum, Ohio, USA. Este método se basa en el modelo del
almacenamiento natural en un cauce, y el amortiguamiento de la onda de
crecida por efecto de rozamiento en el fondo y las márgenes del canal natural.
Almacenamiento Prismático: Formado por un volumen de sección transversal
constante a lo largo del cauce. Almacenamiento en Cuña: Formado por la
diferencia entre los caudales de entrada y salida, o bien por la pendiente de la
lámina de agua en el tramo considerado. Durante el avance de la avenida, el
caudal de entrada es mayor que el de salida y se forma lo que se denomina
cuña positiva y durante la recesión, el caudal de entrada es menor al caudal de
salida, formándose una acuña negativa. Para su aplicación es necesario
registrar dos hidrogramas al mismo tiempo, en dos estaciones próximas en el
cauce del río.
El volumen de almacenamiento prismático es proporcional al caudal de salida,
ya que se supone que el caudal de salida es proporcional al área de la sección
del cauce, por lo tanto:

𝑆 𝑝𝑟𝑖𝑠𝑚𝑎 = 𝑆𝑝 = 𝐾𝑂 (5)

Es la diferencia entre las entradas y las salidas:

𝑆𝑐𝑢ñ𝑎 = 𝑆𝑐 = (𝐼 − 𝑂)

𝑆 = 𝑆𝑝 + 𝑆𝑐

𝑆 = 𝐾𝑂 + 𝐾𝑋(𝐼 − 𝑂)

𝑆 = 𝐾[𝑂 + 𝑋(𝐼 − 𝑂)]

𝑆 = 𝐾[𝑂 + 𝑋𝐼 − 𝑋𝑂]

𝑆 = [𝑋𝐼 + (1 − 𝑋)]

Dónde:

S = almacenamiento en el tramo considerado de un cauce

I = caudal de entrada en el tramo

O = Caudal de salida del tramo

K = Representa el tiempo de tránsito de la onda de avenida a través del tramo

X = Depende del almacenamiento en cuña, y varía entre 0 (almacenamiento

tipo embalse) a 0.5 que corresponde a una cuña completamente desarrollada.
Valores prácticos pueden tomarse entre 0.2 – 0.35; 0.2 para tramos situados en
la parte baja de la cuenca y valores mayores para tramos en la parte alta de la
cuenca.

Se tiene que:
𝑂2 = 𝐶1𝐼1 + 𝐶2𝐼2 + 𝐶3𝑂1

4.4. HIDROGRAMAS DE TAMBOGRANDE Y PIURA

5. RESULTADOS:
En este punto se desarrollará lo explicado por el docente Ing. Héctor
Yauri Quispe y tratar de hallar K y X, también hacer un hidrograma a
partir de los datos de Tambogrande hasta Piura.

LIMITES DE LA ESCORRENTIA BASE

TAMBOGRANDE
5500000

5400000

5300000

5200000

5100000
y = 33957x + 5E+06
5000000 R² = 1
4900000

4800000
0 5 10 15 20

Tambogrande
2 4899789
18 5443093

Limites de la Escorrentia Base

PIURA
7000000
6000000
5000000
4000000
3000000 y = 120506x + 3E+06
R² = 1
2000000
1000000
0
0 5 10 15 20 25 30
 Piura
1 3495600
27 6628743

TAMBOGRANDE

12000000.0

10000000.0

8000000.0

6000000.0

4000000.0

2000000.0
Tambogrande
.0 2 4899789
1 3 5 5443093
18 7 9 11 13 15 17

Límites de la
X Hora Escorrentía Base Caudal total (m3/h) Escorrentía directa
(Y)

1 26 4806000.0
2 2.00 5067914.000 4899789 0
3 4.00 5101871.000 5953253 851382
4 6.00 5135828.000 7351260 2215432
5 8.00 5169785.000 8660081 3490296
6 10.00 5203742.000 9895533 4691791
7 12.00 5237699.000 10726530 5488831
8 14.00 5271656.000 11212751 5941095
9 16.00 5305613.000 11370950 6065337
10 18.00 5339570.000 11313739 5974169
11 20.00 5373527.000 11209319 5835792
 12 22.00 5407484.000 10668562 5261078
13 27 5441441.000 9754202 4312761
14 2.00 5475398.000 8786368 3310970
15 4.00 5509355.000 7933896 2424541
16 6.00 5543312.000 7060574 1517262
17 8.00 5577269.000 6190433 613164
18 10.00 5611226.000 5443093 0

PIURA

14000000
12000000
10000000
8000000 y = 33957x + 5E+06
6000000 R² = 1

4000000
2000000
0
1 3 5 7 9 1113151719212325
 Límites de la
X Hora Escorrentía Base Caudal total (m3/h) Escorrentía directa
(Y)

1 Dia 26 3495600
2 2.00 3241012 3643579 402567
3 4.00 3361518 4173706 812188
4 6.00 3482024 4666100 1184076
5 8.00 3602530 5101553 1499023
6 10.00 3723036 5409076 1686040
7 12.00 3843542 5558870 1715328
8 14.00 3964048 5948502 1984454
9 16.00 4084554 6338395 2253841
10 18.00 4205060 6610421 2405361
11 20.00 4325566 6926736 2601170
12 22.00 4446072 7412193 2966121
13 Dia 27 4566578 7966432 3399854
14 2.00 4687084 8445430 3758346
15 4.00 4807590 9006072 4198482
16 6.00 4928096 9947049 5018953
17 8.00 5048602 10862134 5813532
18 10.00 5169108 11508961 6339853
19 12.00 5289614 12188286 6898672
20 14.00 5410120 11975185 6565065
21 16.00 5530626 11383685 5853059
22 18.00 5651132 11069027 5417895
23 20.00 5771638 10685041 4913403
24 22.00 5892144 10180787 4288643
25 Dia 28 6012650 9355074 3342424
26 2.00 6133156 8038680 1905524
27 4.00 6253662 6628743 0
TRANSITO y = 120506x + 3E+06
R² = 1
8000000

7000000

6000000

5000000

4000000

3000000

2000000

1000000

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

Series1 Series2
 C1 0.57894737 C1 0.57894737
C2 0.47368421 C2 0.47368421
- -
C3 C3
0.05263158 0.05263158

SALIDA PIURA SIMULADO K= 1

x= 0.1
∆t = 2
SALIDA
INGRESO INGRESO
Hora Piura
Piura Tambogrande
simulado
1 Dia 26 0 0
C1
2 2.00 0.00 402567 0 0.57894737
3 4.00 403286.21 812188 851382
4 6.00 1521094.94 1184076 2215432 C2 0.47368421
5 8.00 2855859.00 1499023 3490296 -
C3
6 10.00 4092816.63 1686040 4691791 0.05263158
7 12.00 5100861.23 1715328 5488831
8 14.00 5723480.78 1984454 5941095
9 16.00 6011399.85 2253841 6065337
10 18.00 6024990.96 2405361 5974169
11 20.00 5905947.16 2601170 5835792
12 22.00 5559866.68 2966121 5261078
13 Dia 27 4796149.49 3399854 4312761
14 2.00 3812786.92 3758346 3310970
15 4.00 2864671.16 4198482 2424541
16 6.00 1971612.52 5018953 1517262
17 8.00 1065091.87 5813532 613164
18 10.00 298932.22 6339853 0
19 12.00 -15733.27 6898672
20 14.00 828.07 6565065
21 16.00 -43.58 5853059
22 18.00 2.29 5417895
23 20.00 -0.12 4913403
24 22.00 0.01 4288643
25 Dia 28 0.00 3342424
26 2.00 0.00 1905524
27 4.00 0.00 0
 SALIDA Piura simulado
7000000

6000000

5000000

4000000

3000000

2000000

1000000

0
0 5 10 15 20 25 30
-1000000

6. CONCLUSIONES

 La aplicación del método Muskingum, permite conocer de manera

particular las variables y registros de una estación hidrométrica con el fin
de obtener información histórica suficiente para modelaciones
hidrológicas, cuyo análisis y resultados permiten tomar decisiones en
problemas de ingeniera hidráulica de corrientes naturales de flujo no
permanente.

 El método de Muskingum tiene iteraciones que son complejas, lo cual se

tiene que hacer uso de programas especializados y para este caso se
usó Excel.

 El método de Muskingum permite hallar la crecida de ondas en un punto

del rio, mientras se tenga registro hidrométrico.
 Se verifica que el caudal de entrada es mayor que el caudal de salida
durante el avance de la onda de creciente y durante la recesión el
caudal de salida es mayor al caudal de entrada que resulta una cuña
negativa.

 Los valores de las constantes “X” y “K” encontrados para cada una de
las estaciones Tambogrande y Sánchez Cerro – Piura, son de X=0.1
K=1 horas, lo que indica que se trata de cuencas con baja pendiente y
gran caudal, puesto que valores de “X” entre 0.3 a 0.5, son ríos de zonas
de alta pendiente y bajo caudal.

 La diferencia entre los valores de C1, C2 y C3, muestran cierta

uniformidad en los resultados, por lo que se puede concluir que cada
crecida es particularmente diferente. Por otro lado, como lo dicta el
método de Muskingum, en todos los casos se cumple, que la sumatoria
de los coeficientes C1, C2 y C3 es igual a la unidad.