TRÁNSITO DE AVENIDAS EN RÍOS

El tránsito de avenidas en ríos es un procedimiento que sirve para

determinar el hidrograma de salida en el tramo de un río sometido
a un hidrograma de entrada.

Para determinar la variación del volumen en el tramo del río,

también se cumple la ecuación:
𝑑𝑆
𝐼(𝑡) − 𝑄(𝑡) =
𝑑𝑡
Donde:
I: Caudal de entrada
Q: Caudal de salida
𝑑𝑆
: Variación de volumen almacenado en el tiempo.
𝑑𝑡

Para un tiempo dt, el volumen almacenado será el siguiente:

𝑑𝑆
𝐼(𝑡)𝑑𝑡 − 𝑄 (𝑡)𝑑𝑡 = 𝑑𝑡
𝑑𝑡
Para un intervalo de tiempo: Δt comprendido entre tJ+1 y tJ,
tenemos:
𝐼𝐽+1 + 𝐼𝐽 𝑄𝐽+1 + 𝑄𝐽
𝑆𝐽+1 − 𝑆𝐽 = ( ) ∆𝑡 − ( ) ∆𝑡
2 2

Para el modelamiento de una creciente en un río se aplica el

Método de Muskingum, el cual considera que el almacenamiento
resulta de la combinación de un prisma y una cuña, tal como se
observa en la figura:
Durante el avance de la onda creciente, el caudal de entrada es
mayor que el caudal de salida, siendo un almacenamiento de cuña.
Durante la recesión, el caudal de salida es mayor que el caudal de
entrada resultando una cuña, pero negativa. Adicionalmente,
existe un almacenamiento por prisma que está formado por un
volumen de sección transversal constante a lo largo de la longitud
del canal prismático.

Suponiendo que el área de la sección transversal del flujo

creciente es directamente proporcional al caudal en la sección, el
volumen de almacenamiento por prisma sería igual a KQ, donde K
es un coeficiente de proporcionalidad que representa el tiempo de
tránsito de una onda de creciente a través del tramo del río. Por
otro lado, el volumen de almacenamiento por cuña sería igual a:

𝐾 𝑋 (𝐼 − 𝑄)
Donde X es un factor de ponderación que depende de la forma
modelada de almacenamiento por cuña. El valor de X varía desde
0, para un almacenamiento tipo embalse (No tiene remanso), hasta
0,5 para una cuña completamente desarrollada:
0 ≤ 𝑋 ≤ 0,5
En cauces naturales, X se encuentra entre 0 y 0,3 con un valor
medio cercano a 0,2.

El almacenamiento total por lo tanto es la suma de ambos

componentes:
𝑆 = 𝐾𝑄 + 𝐾 𝑋 (𝐼 − 𝑄)
La cual puede reordenarse para dar la función de almacenamiento
para el método de Muskingum:

𝑆 = 𝐾 [𝑋 𝐼 + (1 − 𝑋) 𝑄]
la cual representa un modelo lineal para el tránsito de caudales en
cursos de agua. El almacenamiento en el tiempo J y J+1 puede
expresarse como:
𝑆𝐽 = 𝐾[𝑋𝐼𝐽 + (1 − 𝑋)𝑄𝐽 ]

𝑆𝐽+1 = 𝐾[𝑋𝐼𝐽+1 + (1 − 𝑋)𝑄𝐽+1 ]

Luego, el cambio de almacenamiento durante este intervalo de
tiempo será:

𝑆𝐽+1 − 𝑆𝐽 = 𝐾{[𝑋𝐼𝐽+1 + (1 − 𝑋)𝑄𝐽+1 ] − [𝑋𝐼𝐽 + (1 − 𝑋)𝑄𝐽 ]}

….………. (I)
Sabemos también que podemos expresarlo como:

𝐼𝐽+1 + 𝐼𝐽 𝑄𝐽+1 + 𝑄𝐽
𝑆𝐽+1 − 𝑆𝐽 = ( ) ∆𝑡 − ( ) ∆𝑡
2 2
…………. (II)
Igualando las ecuaciones (I) y (II) y simplificando, obtenemos:

𝑄𝐽+1 = 𝐶1 𝐼𝐽+1 + 𝐶2 𝐼𝐽 + 𝐶3 𝑄𝐽

La cual es la ecuación de tránsito por el método de Muskingum,

donde:
∆𝑡−2𝐾𝑋 ∆𝑡+2𝐾𝑋 2𝐾 (1−𝑋)−∆𝑡
𝐶1 = 𝐶2 = 𝐶3 =
2𝐾 (1−𝑋)+∆𝑡 2𝐾 (1−𝑋)+∆𝑡 2𝐾 (1−𝑋)+∆𝑡
En la cual se cumple que: 𝐶1 + 𝐶2 + 𝐶3 = 1

Si se tienen los caudales de entrada y salida para el tramo de un

río, se pueden hallar los valores de K y X. De las ecuaciones (I) y
(II) obtenemos:

0,5∆𝑡[(𝐼𝐽+1 + 𝐼𝐽 ) − (𝑄𝐽+1 + 𝑄𝐽 )]
𝐾=
𝑋(𝐼𝐽+1 − 𝐼𝐽 ) + (1 − 𝑋)(𝑄𝐽+1 − 𝑄𝐽 )

Para hallar K y X, se asignan valores de X= 0; 0,1; 0,15; 0,20; etc.

y se grafican los puntos tomando el numerador como ordenada y
el denominador como abscisa. Esto produce un gráfico en forma
de bucle. El valor de X será aquel que produzca un bucle más
parecido a una línea recta, luego, por definición el valor de K será
igual a la pendiente de dicha recta.

Ejemplo 1: Determine el hidrograma de salida para un tramo de

un río sometido al hidrograma de entrada mostrado. Se sabe que el
caudal inicial de salida es 92 m3/s. k=2,3 h y X = 0,15

t I (m3/s) t I (m3/s) t I (m3/s) t I (m3/s)

01:00 92 06:00 560 11:00 634 16:00 247
02:00 117 07:00 640 12:00 571 17:00 184
03:00 180 08:00 680 13:00 477 18:00 134
04:00 320 09:00 700 14:00 390 19:00 108
05:00 450 10:00 675 15:00 329 20:00 90

I (m3/s)
800
700
600
Caudal (m3 /s)

500
400
300
200
100
0
0 2 4 6 8 10 12 14 16 18 20
t (h)
Calculamos los coeficientes:
C1 = 0,0631 C2 = 0,3442 C3 = 0,5927

Reemplazamos estos coeficientes en la expresión de Muskingum:

𝑄𝐽+1 = 0,0631 𝐼𝐽+1 + 0,3442 𝐼𝐽 + 0,5927 𝑄𝐽

Obteniendo los resultados siguientes:

t I (m3/s) J Q (m3/s)
01:00 92 1 92.000
02:00 117 2 93.578
03:00 180 3 107.093
04:00 320 4 145.622
05:00 450 5 224.849
06:00 560 6 323.494
07:00 640 7 424.871
08:00 680 8 515.017
09:00 700 9 583.477
10:00 675 10 629.359
11:00 634 11 645.362
12:00 571 12 636.759
13:00 477 13 604.044
14:00 390 14 546.809
15:00 329 15 479.092
16:00 247 16 412.785
17:00 184 17 341.286
18:00 134 18 274.068
19:00 108 19 215.378
20:00 90 20 170.507

800
700
600
Caudal (m 3 /s)

500
400
300
200
100
0
0 2 4 6 8 10 12 14 16 18 20
t (h)

I (m3/s) Q (m3/s)
Ejemplo 2: A partir del registro de los hidrogramas de entrada y
salida medidos en un tramo de un río se pide determinar el valor
de las constantes K y X.

t (min) I (m3/s) Q (m3/s) t (min) I (m3/s) Q (m3/s)

0 0 0 75 951 968
3 60 0 78 925 965
6 120 13 81 890 956
9 180 42 84 852 938
12 240 81 87 810 919
15 300 127 90 767 884
18 364 178 93 717 851
21 446 231 96 668 812
24 530 293 99 618 769
27 613 363 102 566 725
30 696 437 105 514 677
33 776 514 108 462 629
36 855 593 111 410 579
39 932 672 114 359 528
42 948 757 117 309 478
45 932 822 120 261 427
48 914 861 123 248 373
51 911 879 126 238 332
54 921 888 129 229 302
57 941 897 132 222 278
60 958 910 135 216 260
63 975 924 138 210 246
66 982 940 141 205 235
69 980 954 144 199 225
72 969 964 147 194 217

Aplicando la expresión:

0,5∆𝑡[(𝐼𝐽+1 + 𝐼𝐽 ) − (𝑄𝐽+1 + 𝑄𝐽 )]
𝐾=
𝑋(𝐼𝐽+1 − 𝐼𝐽 ) + (1 − 𝑋)(𝑄𝐽+1 − 𝑄𝐽 )

De los datos vemos que Δt = 3 min = 180 s, luego, obtenemos los

valores correspondientes al numerador y denominador para X =
0,0; 0,1; 0,2; 0,3; 0,4 y 0,5
 t I Q X= 0.00 0.1 0.2 0.3 0.4 0.5
J
(min) (m3/s) (m3/s) Num. Den. Den. Den. Den. Den. Den.
1 0 0 0
2 3 60 0 5400 0.00 6.00 12.00 18.00 24.00 30.00
3 6 120 13 15030 13.00 17.70 22.40 27.10 31.80 36.50
4 9 180 42 22050 29.00 32.10 35.20 38.30 41.40 44.50
5 12 240 81 26730 39.00 41.10 43.20 45.30 47.40 49.50
6 15 300 127 29880 46.00 47.40 48.80 50.20 51.60 53.00
7 18 364 178 32310 51.00 52.30 53.60 54.90 56.20 57.50
8 21 446 231 36090 53.00 55.90 58.80 61.70 64.60 67.50
9 24 530 293 40680 62.00 64.20 66.40 68.60 70.80 73.00
10 27 613 363 43830 70.00 71.30 72.60 73.90 75.20 76.50
11 30 696 437 45810 74.00 74.90 75.80 76.70 77.60 78.50
12 33 776 514 46890 77.00 77.30 77.60 77.90 78.20 78.50
13 36 855 593 47160 79.00 79.00 79.00 79.00 79.00 79.00
14 39 932 672 46980 79.00 78.80 78.60 78.40 78.20 78.00
15 42 948 757 40590 85.00 78.10 71.20 64.30 57.40 50.50
16 45 932 822 27090 65.00 56.90 48.80 40.70 32.60 24.50
17 48 914 861 14670 39.00 33.30 27.60 21.90 16.20 10.50
18 51 911 879 7650 18.00 15.90 13.80 11.70 9.60 7.50
19 54 921 888 5850 9.00 9.10 9.20 9.30 9.40 9.50
20 57 941 897 6930 9.00 10.10 11.20 12.30 13.40 14.50
21 60 958 910 8280 13.00 13.40 13.80 14.20 14.60 15.00
22 63 975 924 8910 14.00 14.30 14.60 14.90 15.20 15.50
23 66 982 940 8370 16.00 15.10 14.20 13.30 12.40 11.50
24 69 980 954 6120 14.00 12.40 10.80 9.20 7.60 6.00
25 72 969 964 2790 10.00 7.90 5.80 3.70 1.60 -0.50
26 75 951 968 -1080 4.00 1.80 -0.40 -2.60 -4.80 -7.00
27 78 925 965 -5130 -3.00 -5.30 -7.60 -9.90 -12.20 -14.50
28 81 890 956 -9540 -9.00 -11.60 -14.20 -16.80 -19.40 -22.00
29 84 852 938 -13680 -18.00 -20.00 -22.00 -24.00 -26.00 -28.00
30 87 810 919 -17550 -19.00 -21.30 -23.60 -25.90 -28.20 -30.50
31 90 767 884 -20340 -35.00 -35.80 -36.60 -37.40 -38.20 -39.00
32 93 717 851 -22590 -33.00 -34.70 -36.40 -38.10 -39.80 -41.50
33 96 668 812 -25020 -39.00 -40.00 -41.00 -42.00 -43.00 -44.00
34 99 618 769 -26550 -43.00 -43.70 -44.40 -45.10 -45.80 -46.50
35 102 566 725 -27900 -44.00 -44.80 -45.60 -46.40 -47.20 -48.00
36 105 514 677 -28980 -48.00 -48.40 -48.80 -49.20 -49.60 -50.00
37 108 462 629 -29700 -48.00 -48.40 -48.80 -49.20 -49.60 -50.00
38 111 410 579 -30240 -50.00 -50.20 -50.40 -50.60 -50.80 -51.00
39 114 359 528 -30420 -51.00 -51.00 -51.00 -51.00 -51.00 -51.00
40 117 309 478 -30420 -50.00 -50.00 -50.00 -50.00 -50.00 -50.00
41 120 261 427 -30150 -51.00 -50.70 -50.40 -50.10 -49.80 -49.50
42 123 248 373 -26190 -54.00 -49.90 -45.80 -41.70 -37.60 -33.50
43 126 238 332 -19710 -41.00 -37.90 -34.80 -31.70 -28.60 -25.50
44 129 229 302 -15030 -30.00 -27.90 -25.80 -23.70 -21.60 -19.50
45 132 222 278 -11610 -24.00 -22.30 -20.60 -18.90 -17.20 -15.50
46 135 216 260 -9000 -18.00 -16.80 -15.60 -14.40 -13.20 -12.00
47 138 210 246 -7200 -14.00 -13.20 -12.40 -11.60 -10.80 -10.00
48 141 205 235 -5940 -11.00 -10.40 -9.80 -9.20 -8.60 -8.00
49 144 199 225 -5040 -10.00 -9.60 -9.20 -8.80 -8.40 -8.00
50 147 194 217 -4410 -8.00 -7.70 -7.40 -7.10 -6.80 -6.50
 X = 0,00
60000
50000
40000
30000
Numerador 20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador

X = 0,10
60000
50000
40000
30000
Numerador

20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador

X = 0,20
60000
50000
40000
30000
Numerador

20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador
 X = 0,30
60000
50000
40000
Numerador 30000
20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador

X = 0,40
60000
50000
40000
30000
Numerador

20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador

X = 0,50
60000
50000
40000
30000
Numerador

20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador
En lugar de graficar podemos aplicar el coeficiente de correlación
para determinar qué valor de X hace que la gráfica se asemeje a
una recta. En el caso que estamos analizando, tenemos:

X: 0.0 0.1 0.2 0.3 0.4 0.5

r= 0.98864 0.99578 0.99926 0.99897 0.99494 0.98729

Vemos que X se encuentra entre 0,2 y 0,3. Evaluando otros

valores de X hallamos que para X = 0,24 la gráfica se asemeja
más a una recta:
t I Q X= 0.24
J
(min) (m3/s) (m3/s) Num. Den.
1 0 0 0
2 3 60 0 5400 14.40
3 6 120 13 15030 24.28
4 9 180 42 22050 36.44
5 12 240 81 26730 44.04
6 15 300 127 29880 49.36
7 18 364 178 32310 54.12
8 21 446 231 36090 59.96
9 24 530 293 40680 67.28
10 27 613 363 43830 73.12
11 30 696 437 45810 76.16
12 33 776 514 46890 77.72
13 36 855 593 47160 79.00
14 39 932 672 46980 78.52
15 42 948 757 40590 68.44
16 45 932 822 27090 45.56
17 48 914 861 14670 25.32
18 51 911 879 7650 12.96
19 54 921 888 5850 9.24
20 57 941 897 6930 11.64
21 60 958 910 8280 13.96
22 63 975 924 8910 14.72
23 66 982 940 8370 13.84
24 69 980 954 6120 10.16
25 72 969 964 2790 4.96
26 75 951 968 -1080 -1.28
27 78 925 965 -5130 -8.52
28 81 890 956 -9540 -15.24
29 84 852 938 -13680 -22.80
30 87 810 919 -17550 -24.52
31 90 767 884 -20340 -36.92
32 93 717 851 -22590 -37.08
33 96 668 812 -25020 -41.40
34 99 618 769 -26550 -44.68
 35 102 566 725 -27900 -45.92
36 105 514 677 -28980 -48.96
37 108 462 629 -29700 -48.96
38 111 410 579 -30240 -50.48
39 114 359 528 -30420 -51.00
40 117 309 478 -30420 -50.00
41 120 261 427 -30150 -50.28
42 123 248 373 -26190 -44.16
43 126 238 332 -19710 -33.56
44 129 229 302 -15030 -24.96
45 132 222 278 -11610 -19.92
46 135 216 260 -9000 -15.12
47 138 210 246 -7200 -12.08
48 141 205 235 -5940 -9.56
49 144 199 225 -5040 -9.04
50 147 194 217 -4410 -7.28

Donde, r = 0,99960 Gráficamente:

X = 0,24
60000
50000
40000
30000
Numerador

20000
10000
0
-10000 -55 -35 -15 5 25 45 65 85
-20000
-30000
-40000
Denominador

La pendiente de esta recta será el valor de K. Aplicando regresión

lineal simple hallamos la ecuación de la recta:

Num. = 599,74 Den. – 84,95

Es decir: K = 599,74 s. K ≈ 10 minutos