Mayo-2012
Trnsito de Hidrogramas
Conceptos bsicos
Se trata de conocer cmo evoluciona un hidrograma a medida que discurre a lo largo de un
cauce o a travs de un depsito o embalse.
Tambin se habla de trnsito de avenidas, y se utilizan las expresiones transitar una avenida
o transitar un hidrograma1.(En ingls Hydrograph Routing, Flood Routing o Flow Routing)
Supongamos que en el extremo de un canal seco arrojamos un volumen de agua (Figura 1). El
pequeo hidrograma generado ser inicialmente ms alto y de menor duracin (posicin A del
dibujo) y, a medida que avanza, el mismo volumen pasar por los puntos B y C cada vez con un
hidrograma ms aplanado. Suponemos que no existe prdida de volumen (por infiltracin o
evaporacin), de modo que el rea comprendida bajo los tres hidrogramas ser idntica.
Fig. 1
Calcular el trnsito de un hidrograma es obtener el hidrograma del punto C a partir del
hidrograma del punto A. La utilidad prctica del procedimiento es evidente. Por ejemplo, el
carcter catastrfico de una avenida est relacionado directamente con la altura del pico del
hidrograma (el caudal mximo), de modo que es fundamental calcular cmo ese pico va
disminuyendo a medida que nos movemos aguas abajo.
Fig. 2 (A) (B)
Si la figura 1 evocaba el proceso que se produce en un ro, tambin se estudia el proceso de
trnsito de caudales en embalses o cualquier depsito con una entrada y una salida. Observando la
figura 2 se comprende que un aumento en el caudal de entrada producir tambin un aumento en el
caudal de salida, pero amortiguado por el depsito. Si en el caudal de entrada ( I ) se produjera un
hidrograma similar al de la Figura 1-A, en el caudal de salida (O) se producira un hidrograma
similar a la Figura 1-B 1-C.
Existen diversos procedimientos para efectuar estos clculos, que se agrupan en dos
categoras:
1
Expresin incorrecta en espaol: el verbo transitar es intransitivo
F. Javier Snchez San Romn---- Dpto. Geologa Univ. Salamanca (Espaa) [Link] Pg. 1
Mtodos Hidrolgicos. Se basan en la ecuacin de la continuidad, que para un tramo de un
cauce (o para un embalse) establece que:
Volumen de entrada en un t - Volumen de salida en ese mismo t= almacenamiento
dividiendo por t:
Q entrada - Q salida = almacenamiento/ t (1)
Con las variables indicadas en la figura 2-B, la igualdad (1) podemos expresarla as:
S S i 1
I O i (2)
t
Siendo:
t = intervalo de tiempo entre los tiempos ti-1 y ti
Si-1 = volumen almacenado en el comienzo del t (tiempo ti-1)
Si = volumen almacenado al final del t (tiempo ti)
I = Caudal medio de entrada (durante el intervalo t)
O = Caudal medio de salida (durante el intervalo t)
Para evaluar los caudales medios de entrada ( I ) y salida (O) dentro de cada t hacemos la
media de los caudales de dos tiempos consecutivos. As, la expresin (2) resultara:
I i 1 I i Oi 1 Oi S i S i 1
(3)
2 2 t
Mtodos hidrulicos2. Adems de la ecuacin de la continuidad, utilizan las ecuaciones del
movimiento del fluido, de modo que para cauces o canales en rgimen no permanente hay que
considerar velocidades, seccin, coeficientes de rugosidad del cauce,...
Todos los modelos (programas de ordenador) utilizados en Hidrologa Superficial incluyen el
clculo del trnsito de hidrogramas. No obstante, siempre conviene saber realizar a mano, aunque
sea para casos sencillos, las tareas que despus encomendaremos a las mquinas.
Trnsito en cauces: Mtodo de Muskingum
Entre los mtodos hidrolgicos, posiblemente el ms utilizado en clculos manuales por su
sencillez sea el de Muskingum3 (Chow et al., 1994, p.264; Singh, V.P, 1992, p.680; Wanielista,
1997, p.323; Viessman, 1995, p. 235).
El almacenamiento (S) en un tramo del cauce puede descomponerse en dos partes:
almacenamiento en prisma, que
sera proporcional al caudal de
salida ( O ) y almacenamiento en
cua, que sera funcin de la
diferencia entre el caudal de
entrada y el de salida (I-O), ya
que cuanto mayor sea esa
diferencia, ms pronunciada ser
la cua:
S prisma= K . O (4a)
S cua= K . X . (I-O) (4b)
Sumando las dos expresiones anteriores, se obtiene:
2
Segn Chow et al. (1993), Mtodos hidrolgicos ="Trnsito agregado de crecientes". Mtodos hidrulicos =
"Trnsito distribudo de crecientes"
3
Muskingum no es el nombre de su autor, sino que el mtodo fue desarrollado en los aos 30 por el Servicio de
Conservacin del distrito de Muskingum (Ohio, USA) para prevencin de avenidas.
F. Javier Snchez San Romn---- Dpto. Geologa Univ. Salamanca (Espaa) [Link] Pg. 2
S = K [X I + (1-X) O] (5)
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
Aplicamos (5) a dos tiempos consecutivos ti-1 y ti, separados por un intervalo t:
Si-1 = K [X Ii-1 + (1-X) Oi-1]
Si = K [X Ii + (1-X) Oi] (6b)
Sustitumos las expresiones (6a) y (6b) en la ecuacin (3) y despejando Oi, resulta:
KX 0,5t KX 0,5t K KX 0,5t
Oi I i I i 1 Oi 1 (7)
K KX 0,5t K KX 0,5t K KX 0,5t
Que para el clculo del caudal de salida para el tiempo ti , se esquematiza as:
Oi = C0 Ii + C1 Ii-1 + C2 Oi-1 (8)
donde: Ii-1 , Oi-1 = Caudales de entrada y salida al final del tiempo anterior
Ii , Oi = Caudales de entrada y salida tras este tiempo
C0 = (-KX + 0,5 t) / (K - KX + 0,5 t) (9a)
C1 = ( KX + 0,5 t) / (K - KX + 0,5 t) (9b)
C2 = (K - KX - 0,5 t) / (K - KX + 0,5 t) (9c)
K, X = constantes que dependen de cada tramo de cauce
Ejemplo: el caudal de salida para tiempo = 6 se calcula as:
tiempo (horas) Q entrada ( I ) Q salida ( O )
5 23 19 Oi
6 31 Oi Oi = C0 31 + C1 23 + C2 19
Puede comprobarse fcilmente (sumando 9a+9b+9c) que C0 + C1+ C2 = 1. Esto es til como comprobacin de
los clculos realizados a mano.
K puede asimilarse al tiempo de recorrido de la onda cinemtica de un extremo a otro del
tramo estudiado. K tiene las mismas unidades que para t (horas o das).
X es una constante que en teora puede estar entre 0 y 0,5, pero normalmente vale 0,2 - 0,3. En
primera aproximacin suele tomarse 0,2. Junto con el valor de K, de ella va a depender la mayor o
menor amortiguacin del hidrograma a lo largo del tramo del cauce. Ms adelante explicamos el
clculo de los parmetros K y X.
Si K= t y X = 0,5, el hidrograma de salida es idntico al de entrada pero desplazado a la
derecha un tiempo igual a K
El t elegido debe estar entre K y 2KX (Wanielista, Singh) o entre K y K/3 (Viessman) 4.
Dentro de estos mrgenes, cuanto menor sea el t , mayor es la precisin del mtodo.
Si conocemos estas dos constantes, K y X, podemos calcular los caudales de salida a partir de
los caudales de entrada. Inversamente, si disponemos de los caudales de entrada y salida para el
mismo hidrograma, podremos calcular las constantes K y X para ese tramo de cauce.
4
Inversamente, y de acuerdo con estos lmites, sit =1 da, K puede valer de 1 a 5 das
F. Javier Snchez San Romn---- Dpto. Geologa Univ. Salamanca (Espaa) [Link] Pg. 3
Clculo de caudales de salida, conocidos K y X
Disponemos de los caudales diarios de entrada en un tramo de un cauce, que aparecen en la
primera columna de la tabla adjunta. Deseamos calcular los correspondientes caudales a la salida de
ese tramo sabiendo que K=1,3 das y X=0,3
Calculamos C0, C1 y C2 mediante las expresiones (9):
C0 = 0,0780 ; C1 = 0,6312 ; C2 = 0,2908
Aplicamos la frmula (8) para cada uno de los caudales de entrada, obteniendo los caudales
que aparecen a la derecha.
Q
Representamos grficamente el hidrograma de Entrada Clculo de Q de salida (I)
entrada y el de salida, aprecindose las dos (O)
caractersticas del trnsito: el retardo (desviado 3 3,00
hacia la derecha) y la atenuacin (el caudal mximo 3 C o * 3 + C 1 * 3 + C 2 * 3 = 3,00
o punta del hidrograma ha disminudo): 5 C o * 5 + C 1 * 3 + C 2 * 3 = 3,16
15 C o * 1 5 + C 1 * 5 + C 2 * 3 , 1 6 = 5,24
50 41 C o * 4 1 + C 1 * 1 5 + C 2 * 5 , 2 4 = 14,19
32 C o * 3 2 + C 1 * 4 1 + C 2 * 1 4 , 1 9 = 32,50
40 Q entrada 19 C o * 1 9 + C 1 * 3 2 + C 2 * 3 2 , 5 0 = 31,13
Q salida 6 C o * 6 + C 1 * 1 9 + C 2 * 3 1 , 1 3 = 21,51
30
Caudal
3 C o * 3 + C 1 * 6 + C 2 * 2 1 , 5 1 = 10,28
20 3 C o * 3 + C 1 * 3 + C 2 * 1 0 , 2 8 = 5,12
3 C o * 3 + C 1 * 3 + C 2 * 5 , 1 2 = 3,62
10 3 C o * 3 + C 1 * 3 + C 2 * 3 , 6 2 = 3,18
3 C o * 3 + C 1 * 3 + C 2 * 3 , 1 8 = 3,05
0
3 C o * 3 + C 1 * 3 + C 2 * 3 , 0 5 = 3,02
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
tiempo (das) 3 C o * 3 + C 1 * 3 + C 2 * 3 , 0 2 = 3,00
Clculo de K y X, conocidos los caudales de entrada y salida
Para calcular estos coeficientes para un tramo de cauce necesitamos conocer dos series de
caudales de entrada y salida simultneos de dicho tramo.
Segn HEC (2000), si disponemos de los hidrogramas de entrada y salida, K puede asimilarse
aproximadamente al tiempo observado entre los centroides de ambos hidrogramas, entre sus puntas
o entre los puntos medios de las curvas de crecida de ambos hidrogramas.
Para el clculo de estos parmetros, el proceso es el siguiente: La expresin (5) puede
considerarse como la ecuacin de una recta con pendiente K. Por tanto, si representamos
grficamente el almacenamiento S en funcin de [XI+(1-X)O] debera obtenerse una recta cuya
pendiente sera K .
El procedimiento consiste en elaborar dicho grfico para diversos valores de X (tpicamente:
0,1; 0,2; 0,3; 0,4) y con el que se obtenga lo ms parecido a una recta se tomar como valor de X.
Despus, la pendiente de dicha recta nos proporcionar K.
Previamente, hemos debido calcular el almacenamiento Si . Para ello, en la frmula (3)
despejamos Si :
I I i Oi 1 Oi
S i t i 1 S i 1 (10)
2 2
Con esta expresin calculamos Si para cada incremento de tiempo a partir de los caudales de
entrada y salida.
F. Javier Snchez San Romn---- Dpto. Geologa Univ. Salamanca (Espaa) [Link] Pg. 4
Ejemplo. Disponemos de los caudales de entrada y salida de un tramo de un cauce, datos que
figuran en las tres primeras columnas de la tabla.
La columna Almac(S) se calcula mediante la frmula (10), y las cuatro ltimas columnas
mediante la expresin K [X I + (1-X) O] para los diversos valores de X que se indican en el
cabecero de cada columna.
X I + (1-X) O
tiempo Entrada (I) Salida(O) (Ii-1+Ii)/2 (Oi-1+Oi)/2 Almac(S) X=0,1 X=0,2 X=0,3 X=0,4
0 0 0 0 0 0 0 0 0 0
1 2 0 1,0 0,0 1,0 0,2 0,4 0,6 0,8
2 6 1 4,0 0,5 4,5 1,5 2,0 2,5 3,0
3 12 3 9,0 2,0 11,5 3,9 4,8 5,7 6,6
4 23 7 17,5 5,0 24,0 8,6 10,2 11,8 13,4
5 38 14 30,5 10,5 44,0 16,4 18,8 21,2 23,6
6 29 24 33,5 19,0 58,5 24,5 25,0 25,5 26,0
7 20 26 24,5 25,0 58,0 25,4 24,8 24,2 23,6
8 14 24 17,0 25,0 50,0 23,0 22,0 21,0 20,0
9 8 19 11,0 21,5 39,5 17,9 16,8 15,7 14,6
10 4 14 6,0 16,5 29,0 13,0 12,0 11,0 10,0
11 0 10 2,0 12,0 19,0 9,0 8,0 7,0 6,0
12 0 6 0,0 8,0 11,0 5,4 4,8 4,2 3,6
13 0 3 0,0 4,5 6,5 2,7 2,4 2,1 1,8
14 0 2 0,0 2,5 4,0 1,8 1,6 1,4 1,2
Representamos grficamente S en funcin de cada una de las cuatro ltimas columnas:
X=0,1 X=0,2 X=0,3 X=0,4
S S S S
60 60 60 60
50 50 50 50
D= 20
40 40 40 40
30 30 30 30
20 20 D= 46 20 20
10 10 10 10
0 0 0 0
0 5 10 15 20 25 30 0 5 10 15 20 25 30 0 5 10 15 20 25 0 5 10 15 20 25 30
X*I+(1-X)*O X*I+(1-X)*O X*I+(1-X)*O X*I+(1-X)*O
En este ejemplo, observamos que la que ms se aproxima a una lnea recta es la obtenida
aplicando X= 0,2. Calculamos la pendiente de esa recta:
Pendiente = K = 46/20 = 1,3
Por tanto, para ese tramo los parmetros son: K = 1,3 ; X = 0,2
Las unidades de K son las mismas utilizadas para el tiempo (en la 1 columna), horas o das.
Mediante el modelo HMS
podemos evaluar los parmetros K y
X por ensayo y error: ejecutando
con distintos valores de K y X, hasta
conseguir un hidrograma de salida
similar a las medidas reales. HMS
puede hacer esto automticamente
mediante una optimizacin.
En HMS no es posible asignar
caudales reales (medidos) a la entrada del
cauce, por lo que el hidrograma de entrada
(generado por HMS) tambin deber estar
calibrado con datos reales.
F. Javier Snchez San Romn---- Dpto. Geologa Univ. Salamanca (Espaa) [Link] Pg. 5
Mtodo de Muskingum- Cunge
Cunge combin mtodos hidrulicos con la simplicidad del mtodo de Muskingum
Calcula las dos constantes utlizadas en el mtodo de Muskingum, K y X, mediante parmetros
hidrulicos del cauce.
K = x / c (10)
1 Q
X 1 (11)
2 BS 0 cx
x = longitud del tramo del cauce considerado
c = celeridad = velocidad media . m
m = aproximadamente 5/3 para cauces naturales amplios
S0 = pendiente media del cauce (adimensional)
Q = caudal
B = anchura del cauce
La correcta aplicacin de este mtodo requiere elegir correctamente el t y el x. 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 (US Army Corps of Engineers, 1994).
Bibliografa
Chow, V.T.; D.R. Maidment & L.W. Mays (1993).- Hidrologa Aplicada. McGraw-Hill, 580 pp.
Singh, V.P (1992).- Elementary Hydrology. Prentice Hall, 973 pp.
Viessman, W. & G. L. Lewis (1995).- Introduction to Hydrology. Harper Collins, 4 ed., 760 pp.
Wanielista, M. (1997).- Hydrology and Water Quality Control 2 edicin. Ed. Wiley
En Internet:
HEC (2000).- Hydrologic Modeling [Link] Reference Manual. 138 pp.
[Link]
Mockus, V. & W. Styner (1972).- National Engineering Handbook Part 630, Chapter 17, 100 pp.
National Resources Conservation Service,
[Link]
US Army Corps of Engineers (1994).- Flood Runoff Analisys, Chapter 9, 24 pp.
[Link]
F. Javier Snchez San Romn---- Dpto. Geologa Univ. Salamanca (Espaa) [Link] Pg. 6
