UNIDAD VIII: Propagación de Crecidas
PROPAGACIÓN DE CRECIDAS EN RÍOS
• Es un procedimiento para determinar el tiempo y la
magnitud del caudal en una sección de un curso de agua
utilizando hidrogramas conocidos en uno o más puntos
aguas arriba;
• Tránsito agregado (hidrológico): el flujo se calcula como
una función del tiempo únicamente, en un lugar
particular;
• Tránsito distribuido (hidráulico): el flujo se calcula como
una función del espacio y del tiempo a través del
sistema;
• Un onda de crecida viajando a lo largo de un canal
aumenta su tiempo de base (disipación) y, si no hay
aumento de volumen, disminuye su cresta (atenuación).
RELACIONES ENTRE CAUDAL
Y ALMACENAMIENTO
ONDA MONOCLINAL ASCENDENTE
Esquema de definición para el análisis de una onda monoclinal ascendente
ONDA MONOCLINAL ASCENDENTE
Onda estacionaria y flujo de exceso (q’).
q = A1 ⋅ (u − v 1 ) = A 2 ⋅ (u − v 2 )
´
Siendo u: celeridad de la onda
A: área de la sección transversal del canal
q=
( A1 ⋅ v1 − A2 ⋅ v 2 ) (q1 − q 2 )
=
A1 − A2 ( A1 − A2 )
La velocidad de una onda monoclinal es una
función de la relación que hay entre el área y el
caudal del río. Dado que la velocidad aumenta
generalmente con el nivel del agua, las curvas
de área caudal son cóncavas hacia arriba.
ONDA MONOCLINAL ASCENDENTE
Relación área – descarga típica para un cauce
y su influencia en la celeridad de la onda
ONDA MONOCLINAL ASCENDENTE
• La velocidad de la onda es mayor que la
velocidad del agua en la mayoría de los
canales;
• A mayor caudal inicial, para un caudal
final dado la onda viajará más rápido;
• Para una onda de pequeña amplitud con
respecto a la profundidad, es válida la
Ley de SEDDON, verificada en el Río
Mississipi. dq 1 dq
u = = ⋅
dA B dy
Siendo B: ancho del canal
ONDA MONOCLINAL ASCENDENTE
• En un Canal rectangular muy ancho a
partir de Chezy y suponiendo la
profundidad igual al radio hidráulico
1 1
v = C ⋅ y 2
⋅S 2
3 1
q = A⋅v = v⋅B ⋅ y = C ⋅B ⋅ y 2
⋅S 2
Donde S: pendiente de la superficie del agua.
dq 3 1 1 3
= ⋅C ⋅ B ⋅ y ⋅ S = ⋅v ⋅ B
2 2
dy 2 2
y sustituyendo en la ecuación de Seddon:
3
u = ⋅v
2
DISTINTOS TIPOS DE ONDAS EN CANALES
• ONDA ABRUPTA DE TRASLACIÓN: la
celeridad tiene componente principal en
las fuerzas de inercia, u = v1 ± g ⋅ D
• ONDA CINEMÁTICA: no se atenúa ni se
dispersa pero sí cambia de forma debido a
la dependencia de u sobre q , y también
implica que q = f ( y ) (Predominio del
almacenamiento, u se deduce de la
ecuación de continuidad)
• ONDA DINÁMICA: depende también de la
influencia de la inercia (Traslación pura)
ONDAS DE TRASLACIÓN
Ejemplo del movimiento de una onda de traslación, río North Plate cerca de
Bridgeport y Lisco, Nebraska, Estados Unidos
ONDAS DE ALMACENAMIENTO
Reducción del caudal debido al almacenamiento en el embalse amortiguador
de Englewood, río Ohio, Estados Unidos
Análisis de Crecidas del Río Salado (Pcia. Sta. Fe)
4500
4000
3500
(m3/seg.)
3000
2500
2000
1500
1000
500
0
E F M A M J J A
Meses
Año 1973 Año 1998 Año 2003
Característica de la crecida del año 2003:
rápido traslado hasta la ciudad de Santa Fe:
de 700 m3/s a 4.000 m3/s en solo 7 días
ECUACIONES DE SAINT VENANT
Movimiento gradualmente variado
ECUACIONES DE SAINT VENANT
para flujo gradualmente variado no permanente
HIPÓTESIS
– Aceleraciones transversales despreciables;
– Coeficiente de CHEZY válido para una sección y
tirante dados.
MOMENTUM o Cantidad de Movimiento:
∂v ∂v
• Fuerzas de Inercia: A ⋅ dx ⋅ ρ ⋅ +v⋅
∂t ∂x
∂h
• Fuerzas Gravitatorias: A ⋅ dx ⋅ γ ⋅
∂x
ρ
• Fuerzas Viscosas: Cf ⋅ χ ⋅ ⋅ v 2 ⋅ dx
2
ECUACIONES DE SAINT VENANT
para flujo gradualmente variado no permanente
Como el sistema está en equilibrio la suma de sus
fuerzas debe ser = 0:
ρ ∂h ∂v ∂v
Cf ⋅ χ ⋅ ⋅ v ⋅ dx+ A⋅ dx⋅γ ⋅ + A⋅ dx⋅ ρ ⋅ + v⋅ = 0
2
2 ∂x ∂t ∂x
Y dividiendo todo por A ⋅ dx ⋅ γ resulta:
1 χ ρ 2 ∂h 1 ∂v ∂v
⋅ Cf ⋅ ⋅ ⋅ v + + +v⋅ =0
2 A γ ∂x g ∂t ∂x
Ecuación de CANTIDAD DE
MOVIMIENTO DE SAINT VENANT
V 2
∂h 1 ∂v ∂v
+ + ⋅ + v⋅ = 0
C 2
⋅R ∂x g ∂t ∂x
∂h 1 ∂v ∂v
V = C ⋅ R ⋅ − − ⋅ + v⋅
∂x g ∂t ∂x
∂h ∂z ∂y ∂y
Siendo: = + = −I 0 +
∂x ∂x ∂x ∂x
∂y 1 ∂v ∂v
V = C ⋅ R ⋅ I 0 − − ⋅ +v⋅
∂x g ∂t ∂x
• Cuando Io = If se produce la forma de ONDA
CINEMÁTICA.
ECUACIÓN DE CONTINUIDAD de SAINT VENANT
• “La variación del flujo en el tramo de
análisis es IGUAL a la variación temporal
de la sección”.
• En el balance de volúmenes que ingresan
y egresan del tramo el agua es
considerada incompresible.
∂Q ∂A
+ =0
∂x ∂t
ECUACIÓN DE CONTINUIDAD de SAINT VENANT
• La Ecuación de Continuidad discretizada permite realizar el
balance de volúmenes ingresados y egresados, detectando el
aumento del nivel de agua y volumen almacenado en el tramo
a medida que aumenta el caudal. Luego durante la etapa de
recesión de la crecida el tramo de río debe producir una
cantidad de agua equivalente a ese volumen almacenado:
∆Q ∆A
⋅ ∆x + ⋅ ∆x = 0
∆x ∆t
∆S
= − ∆Q = I − O
∆t
(I 1 + I 2 ) (O1 + O2 )
S 2 − S1 = ⋅ ∆t − ⋅ ∆t
2 2
MÉTODO DE MUSKINGUM
• Utiliza la CONSTANTE DE ALMACENAMIENTO (K –
en horas -, también referida como constante de
PROPAGACIÓN) y el FACTOR DE PONDERACIÓN
DE LA CUÑA DE ALMACENAMIENTO (x,
adimensional).
S Total = S Pr isma + S cuña
S Total = K ⋅ O + K ⋅ X ⋅ (I − O )
S Total = K ⋅ [X ⋅ I + (1 − X ) ⋅ O ]
Estando 0 ≤ X ≤ 0 , 5
ALMACENAMIENTO POR PRISMA Y POR CUÑA
PERFILES POSIBLES
ANTE EL PASO DE UNA CRECIDA
Método de Muskingum
• La Ecuación de Propagación de MUSKINGUM:
O2 = C 0 ⋅ I 2 + C1 ⋅ I 1 + C 2 ⋅ O1
Donde los Coeficientes resultan:
C0 =
(− K ⋅ x + 0,5 ⋅ ∆t )
[K ⋅ (1 − x ) + 0,5 ⋅ ∆t ]
C1 =
( K ⋅ x + 0,5 ⋅ ∆t )
[K ⋅ (1 − x ) + 0,5 ⋅ ∆t ]
C2 =
[K ⋅ (1 − x ) − 0,5 ⋅ ∆t ]
[K ⋅ (1 − x ) + 0,5 ⋅ ∆t ]
y debe verificarse que: C 0 + C1 + C 2 = 1
CÁLCULO DE LOS PARÁMETROS DE
ALMACENAMIENTO
MÉTODO DE MUSKINGUM
• Los numeradores de los coeficientes Co y C2
deben ser POSITIVOS para que haya
ESTABILIDAD NUMÉRICA en los cálculos, como
condición necesaria, conocida como
CONDICIÓN DE COURANT, lo que da lugar al
análisis de la ZONA DE VALIDACIÓN de dichos
coeficientes:
∆xi
∆t ≤
Ck
Esto que significa que la Celeridad de la Onda
Cinemática ( Ck ) y la longitud del tramo
( ∆x ) de análisis condicionan la elección del
intervalo de tiempo.
EFECTOS DEL CAMBIO EN LA CONSTANTE X
PROPAGACIÓN DE CRECIDAS EN
EMBALSES: Método Directo o de PULS
• También identificado como de PISCINA NIVELADA en
“Hidrología Aplicada” (Ven Te Chow, Maidment y Mays,
1994).
• Información topográfica del vaso del embalse
• Construcción de relaciones H-
H-S (volumen de
almacenamiento) y H-H-Q (caudal erogado por vertedero de
embalse), permiten diseñar relaciones entre
Almacenamiento (S) y Caudal de Salida (O).
• La ecuación de continuidad permite resolver las
incógnitas que plantea la propagación, a partir de conocer:
• HIDROGRAMA DE ENTRADA AL EMBALSE
• CONDICIÓN INICIAL DEL EMBALSE (H S , y O1)
• CONSIDERAR EMBALSE SIMPLE (Vertedero sin regulación,
o compuertas abiertas en una posición fija).
CONSTRUCCIÓN DE LA RELACIÓN
ALMACENAMIENTO (S) - CAUDAL DE
SALIDA (Q)
PROPAGACIÓN DE CRECIDAS EN
EMBALSES: Método Directo o de PULS
S 2 − S1 =
( I1 + I 2 )
⋅ ∆t −
(O1 + O2 )
⋅ ∆t
2 2
S2 S1
2⋅ + O2 = (I 1 + I 2 ) − O1 + 2 ⋅
∆t ∆t
Donde Si representa el almacenamiento de agua en el
embalse, Ii los caudales del hidrograma de la crecida de
entrada al embalse, Oi los caudales de salida por el
vertedero del embalse y ∆t el intervalo de tiempo.
Cuando se trata de un embalse con REGULACIÓN:
S
+ O 2 = (I 1 + I 2 ) − 2 ⋅ O R + 2 ⋅ 1 − O1
S2
2⋅
∆t ∆t
OR representa la salida regulada por accionar de
compuertas y los valores de Oi siguen siendo la salida no
controlada.
VERIFICACIÓN DEL MÉTODO DE PULS
HIDROGRAMA UNITARIO DE CLARK
• Aplicable en cuencas sin información de caudales o con muy escasa información
hidrométrica disponible;
• La propagación en la cuenca produce ATENUACIÓN del pico del hidrograma y
DISIPACIÓN del tiempo de base;
• Un HIETOGRAMA DE PRECIPITACIÓN NETA reemplaza al Hidrograma de
Entrada que planteaba el Método de Muskingum;
• El Tiempo de Retardo (entre la finalización de la lluvia neta y la finalización del
flujo por tierra) se asimila al Tiempo de Concentración de la Cuenca;
• Trazado de CURVAS ISOCRONAS: líneas que unen puntos de igual tiempo de
traslado hasta la desembocadura, definidas en cantidad según relación entre
Tiempo de Concentración y Tiempo Unitario de la Cuenca;
• El Diagrama TIEMPO-
TIEMPO-ÁREA permite definir el HIDROGRAMA UNITARIO
INSTANTÁNEO (para un t = 0);
• La Ecuación de Continuidad y la Ecuación Simplificada de Muskingum (X = 0)
permiten estimar la Constante K, cuando hay información disponible de algunos
hidrogramas de crecidas.
HIDROGRAMA UNITARIO DE CLARK
HIDROGRAMA UNITARIO DE CLARK
Analizar el hidrograma unitario instantáneo
permite deducir una expresión para la constante de
almacenamiento (K) ,
a partir de la ecuación de continuidad:
S = I ⋅ ∆t − Q ⋅ ∆t
Como a partir del punto A del hidrograma cesa el flujo por tierra y no
hay precipitación, entonces S = −Q ⋅ ∆t
o con diferenciales dS = −Q⋅ dt
dS
También expresable como = −Q
dt
Y como el almacenamiento que se produce en la cuenca es lineal, se
desprecia la cuña de almacenamiento (X=0)
Simplificando la expresión de Muskingum resulta S = K ⋅ Q
o también dS = K ⋅ dQ Q
=
Y combinando ambas expresiones K ⋅ dQ = −Q ⋅ dt
K
dQ
−
Por lo tanto dt
HIDROGRAMA UNITARIO DE CLARK
La Ecuación Simplificada de Propagación, resulta:
Q2 = 2 ⋅ C1 ⋅ I 2 + C 2 ⋅ Q1
• El K puede estimarse como 0,5 ⋅ Tc siendo éste
último calculado por: 0 , 593
A
Tc [horas ]= 0 , 335 ⋅
i
• O lo sugerido por LINSLEY:
b⋅L⋅ A
K [horas] =
S
b: coeficiente que varía entre 0,01 y 0,03
L: longitud del Canal Principal (Km)
A: área de drenaje (Km²)
S: pendiente media del canal principal.
HIDROGRAMA UNITARIO DE CLARK
Conversión de un hidrograma instantáneo en un hidrograma de duración finita
HIDROGRAMA UNITARIO DE CLARK
Período C
Diagrama de tiempo Vs. escorrentía para una tormenta larga
TABLA DE CÁLCULO
DEL HIDROGRAMA UNITARIO DE CLARK
T Ai I = Ai Pe/T Ii . C1 Qi . C2 Qi+1
(horas) (km²) (m³/s) (m³/s) (m³/s) (m³/s)
1 A1 I1 I1 . C1 0 (4+5)
2 A2 I2 I2 . C1 (4+5).C2 Q3
3 A3 I3 I3 . C1 Q3 . C2 Q4
4 A4 I4
5 A5 I5
