Mairh 08 Transito Hidrologico de Avenidas
Mairh 08 Transito Hidrologico de Avenidas
Mairh 08 Transito Hidrologico de Avenidas
dS
I (t ) Q (t )
dt
Conocido el hidrograma de entrada l(t), el hidrograma de
salida Q(t) no puede obtenerse directamente de la ecuación
dado que tanto Q como S son incógnitas. Se necesita una
segunda relación que es la función de almacenamiento, para
relacionar S, l y Q; el acoplamiento de la función de
almacenamiento y la ecuación de continuidad proporciona
una combinación de dos ecuaciones y dos incógnitas que
permite resolver el problema.
Función de Almacenamiento
De manera general, la descarga de salida de un sistema
hidrológico depende de la cantidad de agua almacenada en el
propio sistema; la relación entre esas dos magnitudes es
conocida como la función de almacenamiento. Una forma de
determinar la función de almacenamiento es tomando en
consideración que la cantidad de almacenamiento, en la
ecuación de continuidad, aumenta y disminuye con el tiempo en
respuesta a I y Q y a sus tasas de cambio con respecto al
tiempo. De manera general la función de almacenamiento
puede expresarse por la función arbitraria, siguiente:
dI d 2 I dQ d 2 Q
S f I , , 2 ,...Q, , 2 ,...
dt dt dt dt
El sistema se resuelve reemplazando la forma linealizada de
la función de almacenamiento dS/dt en la ecuación de
continuidad e integrar la ecuación diferencial resultante para
obtener Q(t) como una función de l(t).
Otro método de solución es por diferencias finitas, para lo cual
el tiempo se divide en intervalos finitos y la ecuación de
continuidad se resuelve recursivamente desde un tiempo inicial
incluyendo en cada paso la función de almacenamiento.
La forma de la función de almacenamiento depende de la
naturaleza del sistema que está siendo analizando. Por ejemplo
para el tránsito en embalses el almacenamiento puede
asumirse solo como una función no lineal de Q: S = f(Q), donde
la función f(Q) se determina relacionando el almacenamiento y
la salida del embalse con el nivel del agua en éste. Para el
tránsito en ríos el almacenamiento se relaciona linealmente
con I y Q en el método de Muskingum: S = f(I,Q).
La relación entre el caudal de salida y el almacenamiento en un
sistema hidrológico tiene bastante influencia en el tránsito de
avenidas. Esta relación puede ser invariable o variable,
dependiendo si el canal o embalse es ancha y profunda
comparada con su longitud en la dirección de flujo.
En el caso de embalses llenos, donde la velocidad del flujo es
baja y el caudal de salida es fijo para una elevación de nivel de
agua (sin control de compuertas o control fijo de compuertas) la
relación entre el caudal de salida y el almacenamiento puede
considerarse invariable. Si la posición de las compuertas de
control cambia, el caudal y la elevación de la superficie de agua
en la presa cambian, y el efecto se propaga hacia aguas arriba
en el embalse para crear una superficie de agua temporalmente
variable hasta que se establece una nueva elevación de
equilibrio de la superficie de agua a través del embalse.
0 0
0 t 0 t
P
Q R, P Q
R
0 0
0 S 0 S
a) Relación invariable b) Relación variable
Q Q entrada Q Q entrada Q
Q salida Q salida
t t t
Tiempo total Tiempo de redistribución Tiempo de traslación
R
I
Q Q
I
I Q
Q
S j 1 Q j 1 I j I j 1 S j Q j S j 1 Q j 1 I j I j 1 S j Q j
Q j
t 2 2 t 2 t 2 2 t 2
Definiendo la función de S Q I j I j 1
almacenamiento, como: N N j 1 N j Qj
t 2 2
Para calcular Qj+1 a partir de la ecuación anterior, se necesita
una función que relacione el almacenamiento con el caudal de
salida o sea: N con Q. Para desarrollar esta función se usan las
curvas o ecuaciones de cota-volumen y cota-caudal de salida
tal como se muestra en la siguiente figura. La relación entre la
cota de la superficie de agua y el almacenamiento en el
embalse puede determinarse mediante estudios topográficos de
campo. La relación cota-caudal se deduce de las ecuaciones
hidráulicas que relacionan altura de agua y caudal, como las
ecuaciones de descarga para varios tipos de vertederos y de
estructuras de salida descritos anteriormente. El valor de t se
toma como el intervalo de tiempo del hidrograma de caudal de
entrada. Para un valor dado de la cota de la superficie de agua,
se determinan los valores de almacenamiento S y del caudal de
salida Q, luego se calcula el valor de N y se dibuja con el valor
del caudal de salida Q en el eje vertical.
Q de salida Q de Salida
Q Q
Función de
H Altura de Agua S Q Almacenamiento
N
t 2
Almacenamiento
H Altura de Agua
t (h) 0 2 4 6 8 10 12 14 16 18 20
I (m3/s) 5 8 15 30 85 160 140 95 45 15 10
Solución:
Como el caudal de ingreso antes de la AMP es de 5 m3/s,
asumimos que a través del vertedero pasa este caudal inicial;
por lo tanto la altura de agua correspondiente a este caudal
es: Q = 64h3/2, por lo tanto h = (5/64)2/3 = 0.1828 m.
A = 1.0 km2
Conocido el área de la h=3m
superficie de agua a nivel de
la cresta del vertedero y el
A = 0.8+0.2h/3
área correspondiente a h = 3 h
m; se puede determinar, por
interpolación lineal, el área
A = 0.8122 km2
correspondiente a cualquier h = 0.1828 m
valor de h, con la siguiente h = 0.0 m
A = 0.8 km2
Q = f(N)
Del mismo modo 400.00
establecemos la relación 350.00 Q = 0.0004N2 + 0.4397N - 5.9083
Q = f(N) ajustando los 300.00 R² = 0.9994
Cabe mencionar
que la fila
correspondiente a
t = 0, corresponde
a las condiciones
iniciales del
problema, donde:
Q = 5 y N = f(Q)
TRANSITO HIDROLOGICO EN RIOS
Cuña de Cuña de
almacenamiento almacenamiento
Prisma de
Prisma de almacenamiento Prisma de
almacenamiento almacenamiento
Método de Muskingum
Fue desarrollado por McCarthy en 1938 y utiliza la ecuación
de continuidad y la relación de almacenamiento que
depende del caudal de entrada y de salida. El
almacenamiento dentro del tramo en un tiempo dado puede
ser expresado por (Chow, 1959):
S
b xI m / n 1 x Q m / n
am/n
Donde I y Q son relacionadas a a y n de la ecuación de
Manning con a y n constantes. El parámetro x define la
ponderación relativa de I y Q en la determinación del
volumen de almacenamiento en el tramo considerado.
El método de Muskingum asume que m/n = 1 y b/a = K,
resultando en una relación lineal de la forma:
S K xI 1 x Q
Donde K es una constante de proporcionalidad referida al
tiempo de tránsito en el tramo y x es un factor de ponderación
dentro del rango 0 x 0.5. Para el caso de tránsito lineal en
embalses se considera x = 0; en canales uniformes y lisos x =
0.5, lo que implica igual influencia del caudal de entrada y
salida y resulta en una onda de traslación pura. Un valor
típico para la mayoría de canales naturales es x = 0.2. No se
requiere gran precisión en la determinación de x debido a que
los resultados del método son poco sensibles a este valor.
Colocando la ecuación en diferencias finitas, se tiene:
S j 1 S j K xI j 1 I j 1 x Q j 1 Q j
Combinando las ecuaciones se obtiene finalmente:
Kx 0.5t
Q j 1 C1 I j 1 C2 I j C3Q j C1
K Kx 0.5t
Kx 0.5t K Kx 0.5t
C2 C3 C1 + C2 + C 3 = 1
K Kx 0.5t K Kx 0.5t
Ejemplo
Transitar la onda de avenida generada por el hidrograma mostrada en las
columnas (1) y (2) de la siguiente tabla, mediante el método de Muskingum.
Considerar K = 3 h; x = 0.1; t = 1; el caudal de salida inicial es 35.20 m3/s
Solución
Se determinan los coeficientes C1, C2 y C3 : C1 = 0.0625; C2 = 0.2500; C3 =
0.6875; donde: C1 + C2 + C3 = 1. El proceso de tránsito se describe en la
tabla siguiente, la misma que puede ser verificado manualmente o mediante
planilla de cálculo:
Tiempo Entrada C1*Ij+1 C2*Ij C3*Qj Salida
tj (h) Ij (m3/s) Q (m3/s)
(1) (2) (3) (4) (5) (6)
0 35.20 35.20
1 58.70 3.67 8.80 24.20 36.67
2 135.30 8.46 14.68 25.21 48.34
3 272.50 17.03 33.83 33.23 84.09
4 440.90 27.56 68.13 57.81 153.49
5 598.70 37.42 110.23 105.53 253.17
6 670.40 41.90 149.68 174.05 365.63
7 695.10 43.44 167.60 251.37 462.41
8 683.90 42.74 173.78 317.91 534.43
9 620.70 38.79 170.98 367.42 577.19
10 534.60 33.41 155.18 396.82 585.40
11 456.00 28.50 133.65 402.47 564.62
12 386.20 24.14 114.00 388.17 526.31
Tiempo Entrada C1*Ij+1 C2*Ij C3*Qj Salida
tj (h) Ij (m3/s) Q (m3/s)
(1) (2) (3) (4) (5) (6)
13 300.70 18.79 96.55 361.84 477.18
14 235.80 14.74 75.18 328.06 417.98
15 177.90 11.12 58.95 287.36 357.43
16 140.50 8.78 44.48 245.73 298.99
17 112.30 7.02 35.13 205.55 247.70
18 95.30 5.96 28.08 170.29 204.32
19 73.00 4.56 23.83 140.47 168.86
20 60.50 3.78 18.25 116.09 138.12
21 51.40 3.21 15.13 94.96 113.30
22 42.20 2.64 12.85 77.89 93.38
23 35.20 2.20 10.55 64.20 76.95
24 35.20 2.20 8.80 52.90 63.90
25 35.20 2.20 8.80 43.93 54.93
Determinación de los Coeficientes de Muskingum
0.5t I j 1 I j Q j 1 Q j
ecuaciones anteriores:
K
x( I j 1 I j ) 1 x Q j 1 Q
j
Solución
Las tres primeras columnas corresponden a datos medidos en
campo. La columna 4 es el valor acumulado del
almacenamiento obtenido de la diferencia entre I y Q. Las
columnas restantes corresponden a los valores de caudales
ponderados, obtenidos para 4 valores diferentes de x.
El valor del coeficiente x se determina como aquella que
genera el lazo o bucle más concentrada posible alrededor de
una línea recta, en este caso x = 3, y K es la inversa de la
pendiente de la recta de mejor ajuste del bucle
correspondiente a x = 3.
t I Q S [xI + (I - x)Q]
(h) (m3/s) (m3/s) (m3/s-h) para diferentes valores de x
x = 0 x = 0.1 x = 0.2 x = 0.3
1 0.0 0.0 0.0 0.0 0.0 0.0 0.0
2 35.0 10.0 25.0 10.0 12.5 15.0 17.5
3 95.0 50.0 70.0 50.0 54.5 59.0 63.5
4 60.0 60.0 70.0 60.0 60.0 60.0 60.0
5 35.0 50.0 55.0 50.0 48.5 47.0 45.5
6 20.0 35.0 40.0 35.0 33.5 32.0 30.5
7 15.0 25.0 30.0 25.0 24.0 23.0 22.0
8 5.0 15.0 20.0 15.0 14.0 13.0 12.0
9 0.0 10.0 10.0 10.0 9.0 8.0 7.0
10 0.0 5.0 5.0 5.0 4.5 4.0 3.5
11 0.0 0.0 0.0 0.0 0.0 0.0 0.0
Relación entre el almacenamiento y el caudal ponderado
a) x = 0 b) x = 0.1
70.0 70.0
60.0 60.0
50.0 50.0
xI-(1-x)Q
xI-(1-x)Q
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 20.0 40.0 60.0 0.0 20.0 40.0 60.0
S S
c) x = 0.2 d) x = 0.3
70.0
70.0 60.0
Pendiente = 0.8827
60.0 K = 1/0.8827 = 1.13
50.0
50.0
xI-(1-x)Q
xI-(1-x)Q
40.0 40.0
30.0 30.0
20.0 20.0
10.0
10.0
0.0
0.0 20.0 40.0 60.0 0.0
-10.0 0.0 20.0 40.0 60.0
S
S
El método de Muskingum-Cunge;
Básicamente relaciona matemáticamente los valores de K y x en función de las
características físicas del río. El tránsito agregado de crecientes se usa comúnmente
para manejar relaciones caudal-almacenamiento variable. Este método modela el
almacenamiento volumétrico de creciente en un canal de un río mediante la combinación
del almacenamiento de cuña y prisma.
PLANO X
–T
PROCESO DE
CALCULO
RESUTADO
S
RESULTADOS