Principios Y Fundamentos de La Hidrología Superficial - Agustín Breña Puyol & Marco Jacobo Villa (UAM)
Principios Y Fundamentos de La Hidrología Superficial - Agustín Breña Puyol & Marco Jacobo Villa (UAM)
Principios Y Fundamentos de La Hidrología Superficial - Agustín Breña Puyol & Marco Jacobo Villa (UAM)
y
Fundamentos
de la
Hidrología Superficial
http://estudiantesingcivil.blogspot.mx/
CONTENIDO
Presentación ………………………………………………………………………… 3
Introducción ………………………………………………………………………… 6
Prefacio …………………………………………………………………………………… 9
Así mismo, esta publicación aparece dos meses antes del inicio del IV Foro
Mundial del Agua, que se llevará a efecto en la Ciudad de México del 16 al 22
de marzo del 2006, y es congruente con los principios rectores del Foro, así
como con sus objetivos, a saber:
Los autores del libro (Agustín Felipe Breña Puyol y Marco Antonio Jacobo
Villa), la responsable de su adaptación electrónica (Guadalupe Alarcón
Lizardi) y un servidor como editor, se constituyeron como un grupo de
especialistas en su campo de acción, que desean que este trabajo apoye
tanto a los especialistas con todos los recursos disponibles, como a los
estudiantes, profesores e interesados en el tema de los lugares más alejados
del planeta, en donde sea posible consultarlo, ya sea en un disco compacto o
en algún sitio público con acceso a Internet.
4
Extendemos una cordial invitación al lector a consultar el libro y hacer llegar
sus comentarios y sugerencias al Dr. Agustín Breña Puyol o al Dr. Marco
Antonio Jacobo Villa ([email protected]) o a un servidor al correo
[email protected].
5
Introducción
6
Asimismo, para poseer un mayor dominio sobre el proceso de la precipitación
se describen nociones tales como: el role que desempeña una estación
climatológica; los tipos de precipitación que ocurren en una cuenca
hidrológica; los aparatos que se usan para medir la lluvia en forma continua
o discreta; las técnicas que se utilizan para analizar el comportamiento de
registros discretos de lluvia y su consistencia; la descripción pormenorizada
de la distribución geográfica de la precipitación en los diferentes ámbitos
geográficos del país; y algunos otros más. En su fase final, se anexan varios
ejemplos para ilustrar la aplicación de los eventos asociados con el proceso
de la precipitación.
El objetivo del capítulo 4 es describir las interrelaciones que ocurren entre los
procesos del escurrimiento y la lluvia. De acuerdo con el ciclo hidrológico el
escurrimiento se define como la porción de la precipitación que ocurre en una
zona o cuenca hidrológica y que circula sobre o debajo de la superficie
terrestre y que llega a una corriente para ser drenada hasta la salida de una
cuenca o bien alimentar un lago, si se trata de cuencas abiertas o cerradas,
respectivamente.
7
Ahora bien, la infiltración es un proceso complejo ya que esta gobernado por
un conjunto de factores relacionados con: el medio poroso o fracturado del
suelo o bien con una combinación de diferentes medios; los factores que
afectan a la capacidad de infiltración tales como lámina de agua disponible,
conductividad hidráulica y características químicas en la superficie del
terreno, contenido de agua en los estratos de suelo someros y pendiente de
la superficie y rugosidad; y la medición y métodos para calcular la
infiltración.
8
Prefacio
Esta obra tiene como objetivo ofrecer los conceptos teóricos básicos de la
Hidrología Superficial, los cuales serán fundamentales en la adecuada
aplicación de los métodos y técnicas ya aceptados en la Ingeniería
Hidrológica, así como para el desarrollo de nuevas herramientas que ante la
problemática actual se requieren.
En fin, esperamos que este aporte sea de utilidad tanto para los ingenieros
hidrólogos como para aquellos profesionistas interesados en adquirir el
conocimiento elemental de la Hidrología de Superficie, así como para quienes
requieran consultar alguna definición o algún concepto específico relacionado
con el tema.
9
1. Fundamentos básicos
1.1. Definición
1.2. Objetivos
Subsistema atmosférico
Esta precipitación podrá tener lugar tanto en la superficie del océano como
en el continente y, en algunas situaciones, el agua precipitada no se
incorporará a ningún proceso de los subsistemas de agua superficial y
subterránea, regresando a la atmósfera; a esta variable se le da el nombre
de intercepción o pérdidas, y queda ejemplificada por aquella porción de
agua líquida atrapada en las copas de los árboles, techos de las casas, etc.
11
Subsistema del agua superficial
• Con los sistemas geológicos profundos como aquellos que están asociados
a sistemas geotérmicos y que pueden interactuar con los componentes
recarga, flujo subterráneo y flujo subsuperficial.
12
Capas superiores de la atmósfera Sistema adyacente
ES
Subsistema
Ev P It Ev y ET atmosférico
Océano Continente
If FSs Subsistema
de agua
subterránea
R FSb
Figura 1.1. Representación esquemática del sistema global del ciclo hidrológico
13
Tabla 1.1. Volúmenes globales estimados de almacenamiento
Agua Subterránea
- Agua dulce 134.8 10 530 000 0.76 30.10
- Agua salada 134.8 12 870 000 0.93
Lagos
- Agua dulce 1.2 91 000 0.007 0.26
- Agua salada 0.8 85 400 0.006
14
Tabla 1.2. Volúmenes globales estimados de almacenamiento
por componente del ciclo hidrológico.
Precipitación
- Volumen anual km3 /año 458 000 119 000
- Lámina anual mm/año 1 270 800
Evaporación
- Volumen anual km3 /año 505 000 72 000
- Lámina anual mm/año 1 400 484
Escurrimiento total
- Volumen anual km3 /año 47 000
- Lámina anual mm/año 316
15
La ecuación de Reynolds se aplica a las propiedades extensivas (aquellas que
dependen de la masa) del fluido que está bajo estudio. Si N es la propiedad
extensiva total en estudio, entonces:
∂
dN / dt =
∂t ∫ η ρ dV + ∫ η ρ v dA
vc sc
(1.1)
E − S = ± ∆V
(1.2)
Por otra parte, algunos de los procesos de mayor relevancia que ocurren en
el ciclo hidrológico se comentan a continuación:
16
La energía que llega a la superficie terrestre desde el Sol y de los
estratos atmosféricos superiores
Esta energía puede ser absorbida, reflejada o emitida por los diferentes tipos
de superficie terrestre. El porcentaje que se refleja se denomina albedo, y a
la resta de la parte absorbida menos la emitida se le denomina radiación
neta, la cual se expresa como:
R n = R i (1 − α) − R e
(1.3)
4
Re = e σ T
(1.4)
u 1
=
κ ln [z / z 0 ]
*
u
(1.5)
La evaporación
La evapotranspiración
1.5
t
i=
5
(1.6)
18
El escurrimiento superficial
1 2 /3 1/2
v= Rh S
n
(1.7)
Qp = C I A
(1.8)
19
∂θ / ∂t = −∇ [K(ψ)∇ψ ] + ∂K / ∂z
(1.9)
v = − K ∂ h/ ∂ x
(1.10)
20
http://estudiantesingcivil.blogspot.mx/
21
Referencias
22
2. Cuenca hidrológica
23
2.2. Cuencas y regiones hidrológicas en la República Mexicana
24
Para obtener un conocimiento más detallado de los recursos hidráulicos y
hacer más eficaz la administración y la planeación del agua, en el programa
sectorial de mediano plazo denominado Plan Nacional Hidráulico 1995-2000
(CNA, 1996), las 314 regiones hidrológicas con que cuenta el país se
agruparon en 72 subregiones hidrológicas y seis regiones administrativas.
25
Tabla 2.1. Características de las regiones hidrológicas
26
Figura 2.2. Límites geográficos de las regiones hidrológicas
27
Área de la cuenca. Es la proyección del parteaguas a un plano horizontal,
caracterizándose así el tamaño de la cuenca. El valor del área se obtiene de
los mapas topográficos a través del uso del planímetro o de otros métodos.
a) Aspectos generales
Es conveniente indicar que el sistema o red de cauces que drena una cuenca
se clasifican en: dendrítico, rectangular, radial, enrejado y multicuenca. La
figura 2.3 indica la pauta de los diferentes sistemas de cauces.
30
http://estudiantesingcivil.blogspot.mx/
Dendrítico Rectangular
Radial Enrejado
Multicuenca
31
b.2) A partir de la duración de su descarga, los ríos se pueden
clasificar en:
Corriente perenne. Son las que conducen agua durante todo el año.
32
ejemplo para definir el número de orden para cada tramo que conforma la
red de drenaje de una cuenca hidrológica.
Cuenca homogénea; 3 ≤ Rb ≤ 5
33
Figura 2.6. Magnitud de la corriente
Ns
Ds =
A
(2.1)
D s > 500 km / km2 , combinación de un régimen pluvial alto con una litología
fácilmente erosionable.
34
Densidad de drenaje. Es la relación de la suma de las longitudes de las
corrientes (L) entre el valor del área (A) de la cuenca y se calcula con la
expresión:
L
Dd =
A
(2.2)
35
Nx D
Sx =
Lx
(2.3)
Ny D
Sy =
Ly
(2.4)
N D sec θ
Sc = (2.5)
L
∑ (E )
i =1
i
Em =
N
(2.6)
36
donde Em es la elevación media de la cuenca, en msnm; Ei es la elevación i
asociada a la intersección i, en msnm; y N es el número total de
intersecciones.
P
K c = 0.282
A
(2.7)
4 π A
R ci =
P2
(2.8)
37
donde Rci es la relación de circularidad, adimensional; P es el perímetro de la
cuenca, en m; y A es la superficie de la cuenca, en m2.
D A
Re = = 1.128
Lm Lm
(2.9)
38
2.8. Perfil y pendiente del cauce
2
n
S=
1 1 1
+ + .....
S1 S2 Sn
(2.10)
39
2
L
S=
l1 l2 l
+ + ..... n
S1 S2 Sn
(2.11)
2,500
2,400
2,300
Elevaciones, msnm
2,200
2,100
2,000
1,900
0 5 10 15 20 25 30
Distancia, km
40
http://estudiantesingcivil.blogspot.mx/
Solución:
H
S=
L
510
S= = 0.0170 ; esto es
30000
S = 1.7%
41
2
2
L 30,000
S= =
277,163.1714 = 0.0117
l1 l l
+ 2 + ..... n
S1 S2 Sn
S = 0.0117
S = 1.17%
42
Solución:
a) Con el apoyo de los datos de las columnas (1) y (4) de la tabla 2.3, se
elabora la figura 2.8 la cual indica la curva hipsométrica.
1,600
Elevaciones, msnm
1,200
800
400
0
0 10 20 30 40 50 60 70 80 90 100
Porc entaje del área sobre el límite inferior
43
20
15 14.6 14.5
% Área de la cuenca
11.7
11.1
10.6
10 9.3
8.3
7.3
5
3.1 3.0
2.4
2.0 2.0
0
250 350 450 550 650 750 850 950 1,050 1,150 1,250 1,350 1,450
Elevac iones, mnsm
44
Referencias
CNA, Comisión Nacional del Agua (1992). Ley de Aguas Nacionales. Diario
Oficial de la Federación, Comisión Nacional del Agua; México.
CNA, Comisión Nacional del Agua (2004). Estadísticas del Agua en México
2004. Un producto del Sistema Unificado de Información Básica del Agua,
Comisión Nacional del Agua; México.
45
3. Precipitación
46
Posteriormente, al sobrepasar los 90 km, la temperatura aumenta
nuevamente, preservando esta condición hasta los 500 km. A esta región se
le denomina con el nombre de termosfera.
100
T ERMOSFERA
MESOPAUSA
80
MESOSFERA
60
Altitud, en km
40
20
T ROPOPAUSA
T ROPOSFERA
0
-100 -80 -60 -40 -20 0 20
Temperatura, en °C
47
En la troposfera no sólo existen movimientos verticales o convectivos de las
masas de aire; también los hay horizontales o de advección. El movimiento
se puede deber a cambios espaciales y temporales de los elementos
meteorológicos temperatura y presión atmosférica en las proximidades de la
superficie terrestre. La extensión que abarcan estos cambios, así como su
intensidad, da como resultado fenómenos con diferentes escalas, tal como a
continuación se indica.
48
La lluvia, el granizo y la nieve se forman dentro de la misma tormenta. El
que las gotas de lluvia se desplacen hacia la superficie, tiene que ver con el
tamaño; en el caso del granizo, los cristales de hielo son arrastrados por las
fuertes corrientes verticales ascendentes, provocando con ello, un aumento
de tamaño por la colisión con otros núcleos, hasta que el peso es tal que la
atracción gravitacional es mayor que el empuje por la influencia de la
corriente ascendente. En el caso de la nieve, los copos nacen cuando las
gotas de agua atraviesan una capa de aire frío que presenta una
temperatura por debajo de los 0° C y se encuentra cercana a la superficie del
terreno (Llauge, 1976).
Las trombas se forman de una manera similar a los tornados, pero con una
intensidad menor; sin embargo, la gran cantidad de agua precipitada en una
pequeña superficie es lo que hace tan temible a este fenómeno (Simons,
1996).
Los rayos también tienen lugar en el proceso evolutivo de una nube. Son
impronosticables y con periodos de vida de unos cuantos segundos. La
chispa eléctrica que llega a la tierra es quien recibe el nombre de rayo;
mientras que cuando la descarga es entre nubes, se denomina relámpago. Al
sonido que sigue después de la descarga eléctrica se le llama trueno (Llauge,
1976).
Las latitudes medias son aquellas que se encuentran acotadas por los
trópicos y los círculos polares ártico y antártico. En estas zonas, la
interacción entre masas de aire de diferentes condiciones da lugar a
fenómenos que abarcan extensas regiones de la superficie terrestre (de
algunos cientos de miles de kilómetros), denominados frentes.
49
superficie sobre la cual se desplazan, siendo posible clasificarlas según
Dingman (1994) tal como se muestra en la tabla 3.1.
Invierno Verano
Polar c ontinental Temperatura baja Oc éano Ártic o, Canadá, Norte de Canadá, norte de
Humedad baja norte de E.U.A., Eurasia Asia
Tropic al c ontinental Temperatura alta Méxic o, suroeste de E.U.A., Norte de Méxic o, sur
Humedad baja norte de Áfric a, Arabia y suroeste de E.U.A.,
norte de la India Mediterráneo, norte de
Áfric a, Arabia, Asia Central
Polar marítima Temperatura baja Noroeste del Oc éano Norte del Oc éano Atlántic o,
Humedad alta Atlántic o, noroeste del Oc éano Ártic o, norte del
Oc éano Pac ífic o Oc éano Pac ífic o
Tropic al marítima Temperatura alta Parte Central del Atlántic o y Parte Central del Atlántic o y
Humedad alta Pac ífic o, Mar Árabe, Bahía Pac ífic o
de Bengala
Los centros de las masas de aire son áreas de alta presión o anticiclónicas,
con una circulación en el sentido de las manecillas del reloj, en el hemisferio
norte. En el hemisferio sur la circulación es opuesta a las manecillas del reloj.
Por otro lado, los ciclones presentan una circulación en el sentido opuesto a
las manecillas del reloj en el hemisferio norte y en sentido horario en el sur.
50
http://estudiantesingcivil.blogspot.mx/
52
Figura 3.3. Isotermas medias anuales, en °C
53
Unidad orogénic a Región geomórfic a
54
Generalmente se expresa en porcentaje y si la atmósfera presenta un alto
contenido de humedad, la humedad relativa se aproximará al valor de 100%.
Este elemento climatológico también muestra un comportamiento temporal
como la temperatura. En la figura 3.5 se puede observar la variación
geográfica de la humedad relativa en el territorio nacional.
En zonas cercanas a los mares ocurre que, durante el día, el viento se mueve
del continente al cuerpo de agua, y conforme avanza la noche, la dirección
55
se invierte. Caso especial son los vientos asociados a los fenómenos descritos
en la sección 3.1.
56
norte la caída de nieve (caso especial son las montañas más altas, que
permanentemente están cubiertas de nieve), son formas de precipitación
relevantes desde el punto de vista regional.
La figura 3.7 muestra la distribución geográfica de la precipitación pluvial y
para tal efecto se dibujan en un plano de la República Mexicana las isoyetas
(líneas de igual valor de precipitación).
57
Ahora bien, al quedar México comprendido entre los paralelos 14° 31’ y 32°
43’ de latitud norte, toda su extensión se ubica dentro de una de las franjas
de desiertos; se salva de esta fatalidad geográfica gracias a otros factores
que propician la ocurrencia de la precipitación pluvial.
Analizando la figura 3.8 se puede afirmar que los días con precipitación son
mayores en la costra del Golfo de México, Sierra Madre Oriental, Sierra
Madre de Oaxaca, Meseta Central de Chiapas y en la vertiente oriental de la
Cordillera Neovolcánica. En cambio en la costa del Pacífico son pocos los días
con precipitación y lo mismo ocurre con la Sierra Madre Occidental.
58
Heladas. Se denomina helada al proceso de enfriamiento intenso que llega a
producirse a causa de la pérdida nocturna del calor de irradiación terrestre,
situación que se presenta generalmente coincidiendo con la salida del Sol o
pocos minutos después de ella. Tiene mucha importancia sobre la
vegetación, especialmente sobre los cultivos, porque ese descenso de la
temperatura a veces es tan fuerte que provoca la congelación de la savia de
las plantas o de parte de ellas que no tienen eficaz defensa, causándoles
grandes daños e, incluso, la muerte.
59
Figura 3.9. Promedio de días con heladas en el año
61
La humedad relativa se mide por medio de higrómetros; los hay de
absorción, de cabello y electrónicos. También se puede medir a través de un
psicrómetro, el cual está compuesto de dos termómetros: uno de bulbo seco
y otro de bulbo húmedo.
62
Figura 3.11. Estación climatológica automática
La tabla 3.2 (Llauge, 1976) muestra una clasificación de los núcleos por
tamaño representativo y por concentración de partículas en un volumen
dado.
64
La única oportunidad de sobrevivir que tienen las gotas es por la colisión con
otras, situación que favorece al incremento de su volumen; esto prosigue
hasta que el peso es mayor a la influencia de las corrientes ascendentes o a
la evaporación, teniendo lugar la caída como lluvia, granizo o nieve. En su
viaje descendente, el tamaño de las gotas puede aumentar más por
coalescencia, es decir, por la incorporación de gotas más pequeñas que
arrastra a su paso. El tamaño más común de las gotas de lluvia es de
aproximadamente 2.5 mm.
65
• Los pluviómetros pueden estar enterrados o colocados a 1 m de altura
desde la superficie del terreno. El uso de protectores contra viento no
mejora la calidad de la medición.
Dependiendo del objetivo del estudio que se pretenda llevar a cabo, se hará
uso de la información del registro discreto o del continuo. En ambos casos, y
de hecho para cualquier variable hidrológica, se pueden aplicar técnicas para
análisis de valores máximos, mínimos y ordinarios.
66
Figura 3.13. Curva masa de precipitación
67
c) Curvas intensidad-duración-periodo de retorno
n+1
Tr =
m
(3.1)
68
Figura 3.14. Hietogramas de alturas de lluvia e intensidades
Una vez obtenidos los valores máximos por cada duración y por año,
entonces se acomodan de mayor a menor, asignando al valor más grande el
rango 1, y al menor el rango n. Una vez asignado el rango de cada evento,
se procede a estimar el periodo de retorno con el apoyo de la expresión (3.1)
• Método de Chow
hp = a + b log (Tr )
(3.2)
El método más común que permite realizar este proceso esta representado
por la función matemática del tipo siguiente:
m
i = k Tr / dn
(3.3)
o bien
Y = A 0 + A 1 X1 + A 2 X 2
(3.5)
70
http://estudiantesingcivil.blogspot.mx/
A 0 = Ln (k), A 1 = m, A 2 = −n
n n n
∑
i =1
Y = N A 0 + A1 ∑
i =1
X1 + A 2 ∑X
i =1
2
n n n n
∑
i =1
(X1 Y) = A 0 ∑
i =1
X1 + A 1 ∑
i =1
(X12 ) + A 2 ∑ (X
i =1
1) (X 2 )
n n n n
∑
i =1
(X 2 Y) = A 0 ∑
i =1
X 2 + A1 ∑
i =1
(X1 ) (X 2 ) + A 2 ∑ (X
i =1
2
2)
(3.6)
71
3.7. Registros discretos puntuales
n
1
xm =
n
∑x
i=n
i
(3.7)
Para tal proceso, existen tres métodos básicos para realizar el análisis
espacial: precipitación media, curva masa media ajustada y curvas altura de
precipitación-área-duración.
a) Precipitación media
• Media aritmética
72
n
1
hp =
n
∑h
i =n
pi
(3.8)
• Polígonos de Thiessen
n
1
hp =
Ac
∑ (h
i =1
pi Ai )
(3.9)
Con este proceso, cada estación pluviométrica quedara rodeada por líneas
rectas, formando un polígono de forma irregular y el área encerrada por esta
figura será la superficie de influencia de la estación correspondiente Ai.
73
n
1
hp =
Ac
∑ (h
i =1
pi Ai )
(3.10)
74
9. Al graficar los valores afectados por el factor de ajuste, se obtiene la
curva masa media ajustada.
5. Para cada zona encerrada por una isoyeta, empezando por la isoyeta de
mayor valor, se calcula el área encerrada por la isoyeta y la precipitación
media correspondiente;
7. Se calcula una curva masa media para cada área encerrada entre
isoyetas, multiplicando la curva masa de cada estación pluviográfica por el
porcentaje correspondiente obtenido en el inciso anterior;
8. Para cada área encerrada entre isoyetas, se calcula la curva masa media
ajustada (obteniendo el factor de ajuste y multiplicándolo por las curvas
masa);
75
9. Para todas las duraciones de interés, las cuales deben ser múltiplos del
intervalo de tiempo utilizado en el análisis, y para cada área, se calculan
los incrementos máximos de precipitación;
Así, los valores que están fuera de la línea de consistencia, serán afectados
por el factor K.
76
Sharp, 1984). Cabe señalar que existen otros métodos basados en
consideraciones estadísticas y probabilísticas que aquí no se tratan.
1 PA PB PA PC PA PD
PA = + +
3 PB P P
C D
(3.12)
2 P P P P
PI = X T N2 + S2 + E2 + O2
X N X S X E X O
(3.13)
2
donde PI es el valor de precipitación en la estación de interés; X T es la
2 2 2 2
suma cuadrática de las distancias XN , X S , X E y X O las cuales
corresponden a la distancia que va de la intersección de los dos
segmentos de recta perpendiculares, a cada una de las estaciones
ubicadas en uno de los cuatro cuadrantes; PN, PS, PE y PO son los valores
de precipitación para las estaciones ubicadas en los cuadrantes
respectivos.
77
3.9. Distribución geográfica de la precipitación
De la lluvia que recibe nuestro país cada año, apenas una cuarta parte
escurre superficialmente por los cauces de ríos y arroyos. La mayoría se
evapora y una parte menor es absorbida por el suelo o penetra hacia
cavidades subterráneas.
79
Figura 3.15. Distribución geográfica de la precipitación media anual
en la República Mexicana
80
http://estudiantesingcivil.blogspot.mx/
n
∑ (Pxi ) (Pyi ) − n (Px ) (Py )
rxy = i =1
1
n 2 2 n 2 2
2
7. Con los círculos trazados en el mapa, se identifican las zonas que no están
contenidas. Estas son las zonas que requerirán de la instalación de una
estación adicional.
81
b) Coeficiente de variación
2
C
N= v
ε
(3.15)
Este tipo de relación se mide con las curvas altura de lluvia media-área
asociadas a diferentes duraciones, tal como puede observarse en la figura
3.16.
82
Estas curvas describen la variación de la lluvia media a medida que va
aumentando el tamaño del área que abarca una tormenta y existen al
respecto diversos métodos para determinar las curvas altura de precipitación
media-área.
h p(Ai )
FRA =
hp( A = 0)
(3.16)
b) Fórmulas empíricas
0.25
0.25
FRA = 1 − exp − 1.1 d + exp − 1.1 d − 0.01 A
(3.17)
83
c) Método de la curva altura de lluvia media-área
Solución:
84
Tabla 3.3. Discretización de la curva masa de precipitación
0 0.0
1 2.5 2.5
2 6.5 4.0
3 12.5 6.0
4 16.0 3.5
5 20.0 4.0
6 22.0 2.0
7 23.0 1.0
8 24.0 1.0
9 24.0 0.0
10 24.0 0.0
11 24.0 0.0
12 24.0 0.0
85
Ejemplo 3.2. Calcular la altura de precipitación media en una cuenca
hidrológica.
Solución:
a) Método aritmético
86
1
hp = (12 + 9 + 19 + 14 + 23 + 27) = 17.33
6
h p = 17.33 mm
b) Polígonos de Thiessen
La figura 3.20 muestra los polígonos de Thiessen definidos para cada una de
las estaciones, mientras que la tabla 3.4 muestra las áreas parciales y el
proceso para estimar la precipitación media.
87
Tabla 3.4. Áreas parciales de los polígonos de Thiessen
2 2
Estac ión A i , en km hpi , en mm (A i )(hpi ), en km mm
1 4,613 12 55,356
2 1,170 9 10,530
3 2,802 19 53,238
4 4,061 14 56,854
5 3,314 23 76,222
6 1,390 27 37,530
1
hp = (289,730) = 16.7 mm
17,350
h p = 16.7 mm
La figura 3.21 ilustra el trazo de las isoyetas a partir de los datos puntuales
de lluvia registrados en las estaciones, mientras que la tabla 3.5 indica las
áreas parciales y el proceso para estimar la lluvia media.
88
Figura 3.21. Isoyetas de la cuenca hidrológica de estudio
2
Isoyetas A i , en km hpi, en mm (A i)(hpi), en km 2mm
1
hp = (294,635) = 17.0 mm
17,350
h p = 17.0 mm
89
Ejemplo 3.3. Determinar la curva masa media para la cuenca hidrológica
mostrada en la figura 3.22.
Solución:
90
http://estudiantesingcivil.blogspot.mx/
f) Para definir la curva masa media ajustada, los valores de la columna (14)
se multiplican por el factor de ajuste dado por la ecuación:
hp (Isoyetas) 17.0
Fa = = = 1.01796407
hp ( Thiessen)
16.7
91
h) La figura 3.24 muestra la curva masa media ajustada válida para la
cuenca sujeta a este análisis.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14)
Hora hp1 , hp2 , hp3 , hp4 , hp5 , hp6 , A1 hp1 , A 2 hp2 , A3 hp3 , A 4 hp4 , A 5 hp5 , A6 hp6 , hpmedia ,
mm mm mm mm mm mm mmkm
2
mmkm
2
mmkm
2
mmkm
2
mmkm
2
mmkm
2 mm
1 2.5 1.0 0.0 0.0 10.0 5.0 11533 1170 0 0 33140 6950 3.04
2 5.0 2.0 2.5 0.0 15.0 10.0 23065 2340 7005 0 49710 13900 5.53
3 10.0 4.0 7.5 2.5 17.5 15.0 46130 4680 21015 10153 57995 20850 9.27
4 12.0 5.0 12.5 7.5 17.5 20.0 55356 5850 35025 30458 57995 27800 12.25
5 12.0 7.5 13.8 10.0 17.5 22.5 55356 8775 38667.6 40610 57995 31275 13.41
6 12.0 9.0 15.0 12.5 20.0 25.0 55356 10530 42030 50763 66280 34750 14.97
7 12.0 9.0 17.5 14.0 23.0 26.5 55356 10530 49035 56854 76222 36835 16.42
8 12.0 9.0 19.0 14.0 23.0 27.0 55356 10530 53238 56854 76222 37530 16.70
9 12.0 9.0 19.0 14.0 23.0 27.0 55356 10530 53238 56854 76222 37530 16.70
10 12.0 9.0 19.0 14.0 23.0 27.0 55356 10530 53238 56854 76222 37530 16.70
92
Tabla 3.7. Valores de la curva masa media ajustada
Hora hp =( F a ) (hpmedia ),
mm
0 0.00
1 3.10
2 5.63
3 9.44
4 12.47
5 13.65
6 15.24
7 16.71
8 17.00
9 17.00
10 17.00
93
Ejemplo 3.4. Calcular las curvas intensidad-duración-periodo de retorno (i-
d-Tr) para periodos de retorno de 5, 10, 25 y 50 años, en la estación
pluviométrica Departamento del Distrito Federal, D. F.
Solución:
94
Tabla 3.8 Intensidades máximas de lluvia registradas en la estación
pluviométrica Departamento del Distrito Federal, en mm/h
1 25.00 240.0 147.6 129.6 115.5 85.2 59.9 51.0 38.3 30.6 25.5
2 12.50 234.0 140.4 112.8 105.0 76.0 58.5 44.0 33.5 28.2 23.8
3 8.33 184.8 124.8 110.7 105.0 73.5 52.0 43.0 33.0 27.9 23.6
4 6.25 180.0 122.0 104.0 84.0 62.0 52.0 40.6 32.3 26.4 22.0
5 5.00 177.6 120.0 98.8 81.0 60.0 50.1 39.2 32.0 25.8 21.5
6 4.17 175.2 117.0 89.6 76.5 58.5 45.7 39.1 29.6 23.8 20.0
7 3.57 168.0 112.2 80.0 72.0 58.0 45.2 38.2 29.4 23.6 19.7
8 3.13 156.0 102.0 78.0 69.0 56.2 42.7 33.3 25.1 20.4 17.0
9 2.78 144.0 90.0 72.0 69.0 56.0 40.0 32.0 24.0 19.2 16.0
10 2.50 132.0 90.0 71.5 63.5 48.4 37.9 30.7 23.1 18.5 15.4
11 2.27 128.4 87.6 69.6 61.8 48.0 36.0 28.0 21.0 16.8 14.1
12 2.08 120.0 87.0 66.4 60.0 46.0 32.4 26.0 20.0 16.6 14.0
13 1.92 120.0 84.0 66.0 60.0 44.6 30.7 24.5 18.8 16.0 13.9
14 1.79 117.6 82.5 64.0 57.0 41.0 29.7 24.3 18.6 15.0 13.0
15 1.67 108.0 72.0 58.0 47.0 39.0 27.3 22.6 18.0 15.0 12.5
16 1.56 108.0 69.6 52.0 44.4 36.0 26.6 22.3 17.3 14.0 12.5
17 1.47 96.0 66.0 52.0 42.0 34.0 25.3 21.4 16.8 13.8 12.1
18 1.39 96.0 66.0 49.8 42.0 32.0 24.1 20.6 16.7 13.4 12.0
19 1.32 90.0 63.0 48.0 40.5 31.0 24.0 20.0 16.2 13.2 11.5
20 1.25 84.0 60.6 48.0 39.9 31.0 24.0 20.0 15.8 13.1 11.3
21 1.19 84.0 60.0 44.0 37.5 30.0 23.9 19.0 15.0 12.6 11.0
22 1.14 84.0 54.0 44.0 36.0 28.0 22.6 18.0 14.6 12.2 10.8
23 1.09 75.6 48.0 40.0 36.0 28.0 21.3 17.7 13.8 11.5 10.2
24 1.04 39.6 26.1 22.0 18.0 14.0 11.5 9.7 7.8 7.4 7.2
95
0.42
259.823 Tr
i= 0.66
d
(3.19)
96
Tabla 3.9. Procedimiento para determinar las sumas del sistema de ecuaciones
(3.6)
2 2
n Y X1 X2 (X 1 ) (X2 ) X1 X2 X1 Y X2 Y
Ln i Ln T r Ln d
97
Figura 3.25. Curvas intensidad-duración-periodo de retorno de la estación
pluviográfica Departamento de Distrito Federal
98
Figura 3.26. Tormenta convectiva discretizada asociada a una duración de 24 horas
Solución:
90 - 80 10 10
80 - 70 15 25
70 - 60 25 50
60 - 50 30 80
50 - 40 42 122
40 - 30 45 167
30 - 20 180 347
20 - 10 105 452
99
b) Se calcula para cada porción de área asociada entre dos isoyetas
adyacentes, la altura de lluvia media con el procedimiento siguiente:
2
Si A = 0 km ; h p = 90 mm
2 90 + 80
Si A = 10 km ; hp = = 85.0 mm
2
2 15 (75) + 10 (85)
Si A = 25 km ; hp = = 79.0 mm
25
2 25 (65) + 25 (79)
Si A = 50 km ; hp = = 72.0 mm
50
2 30 (55) + 50 (72)
Si A = 80 km ; hp = = 65.6 mm
80
2 42 (45) + 80 (65.6)
Si A = 122 km ; hp = = 58.5 mm
122
100
http://estudiantesingcivil.blogspot.mx/
Tabla 3.11. Procedimiento para estimar el Factor de Reducción por Área (FRA)
0 90.0 1.00
10 85.0 0.94
25 79.0 0.88
50 72.0 0.80
80 65.6 0.73
122 58.5 0.65
167 52.2 0.58
347 38.1 0.42
452 32.7 0.36
101
e) Con los resultados de la tabla 3.11 se construye la curva del Factor de
Reducción por Área (FRA) para la tormenta analizada lo que se indica en
la figura 3.28.
Figura 3.28. Factor de Reducción por Área (FRA) de una tormenta convectiva para
una duración de 24horas en el Valle de México
Con los datos de lluvia registrados en las estaciones se trazaron las isoyetas
y en la figura 3.29 se observa su distribución en la cuenca de estudio.
Asimismo, se trazaron los polígonos de Thiessen con el apoyo de las seis
estaciones pluviográficas. En la figura 3.30 se ilustra los polígonos irregulares
de los pluviógrafos.
102
Por otra parte, la tabla 3.12 presenta las curvas masa de precipitación de los
6 pluviógrafos en intervalos de 6 horas, mientras que la tabla 3.13 indica los
valores de las áreas de los polígonos de Thiessen y de las alturas de
precipitación para una duración de 24 horas.
A 144 887
B 102 1,494
C 54 1,244
D 64 1,888
E 53 837
F 43 995
103
Tabla 3.14. Estimación de la altura de precipitación media-área para el plano
de isoyetas de la figura 3.29 del ejemplo 3.6
104
Figura 3.30. Polígonos irregulares de los pluviógrafos del ejemplo 3.6
Solución:
105
hp A = 144 (1.0) = 144.0 mm
Isoyeta 140 mm ,
hp m = hp A = 144.0 mm
106
Figura 3.31. Estaciones que quedan comprendidas dentro de las isoyetas 60, 100 y
140 mm, respectivamente, superponiendo el plano de polígonos de Thiessen
107
e) Se obtiene la curva masa media (CMM) al sumar las alturas de
precipitación para cada intervalo, multiplicadas por el porcentaje de
influencia de cada estación (columna 8 de la tabla 3.15).
109
Tabla 3.16. Valores de las precipitaciones máximas para las
duraciones de 6, 12, 18 y 24 horas
110
http://estudiantesingcivil.blogspot.mx/
111
Figura 3.32. Curvas altura de precipitación-área-duración
Solución:
112
estación sujeta a este estudio y su lluvia anual acumulada; la columna 10
indica la lluvia anual corregida en la estación.
1 2 3 4 5 6 7 8 9 10
Altura de lluvia anual, en mm Lluvia anual, en mm Estac ion Tepames
Año Estac ión Estac ión Estac ión Suma tres Media Media Lluvia anual, en mm
Buenavista Coquimatlán Ixtlahuac án estac iones ac umulada Registrada Ac umulada Corregida
1 914 857 1,426 3,197 1,066 1,066 1,168 1,168 1,168
2 888 532 741 2,161 720 1,786 755 1,923 755
3 1,082 808 916 2,805 935 2,721 760 2,682 760
4 1,035 931 1,065 3,031 1,010 3,731 1,088 3,770 1,088
5 1,255 984 696 2,935 978 4,710 1,272 5,043 1,272
6 1,178 848 394 2,419 806 5,516 651 5,693 651
7 1,702 948 893 3,543 1,181 6,697 360 6,053 586
8 1,208 890 1,087 3,185 1,062 7,759 1,151 7,204 1,874
9 1,018 848 836 2,702 901 8,659 715 7,919 1,164
10 731 746 598 2,075 692 9,351 509 8,428 828
11 1,057 766 1,337 3,160 1,053 10,404 603 9,031 982
12 952 1,096 1,547 3,596 1,199 11,603 370 9,401 602
113
Figura 3.33. Curva masa doble de precipitación
Solución:
114
hp
i=
d
i d (min)
hp =
60
60 hp
i=
d
(3.21)
0.42 0.34
hp = 4.33 Tr d
(3.22)
115
Figura 3.34. Curvas altura de precipitación-duración-periodo de retorno
(hp-d-Tr) para la estación Departamento del Distrito Federal
116
Referencias
CNA, Comisión Nacional del Agua (1994). El agua y sus senderos. Ríos, lagos
y cascadas de México, México.
CNA, Comisión Nacional del Agua (2004). Estadísticas del Agua en México
2004. Un producto del Sistema Unificado de Información Básica del Agua,
Comisión Nacional del Agua, México.
117
Simons, P. (1996). Weird Weather. Editorial Little, Brown and Company,
Inglaterra.
118
4. Escurrimiento
119
aportarán cantidades suficientes para mantener por algunas semanas más,
después de las últimas lluvias, el gasto en un cauce, formando así las
corrientes intermitentes.
120
http://estudiantesingcivil.blogspot.mx/
4.3. Hidrogramas
121
La figura 4.1 representa un hidrograma anual y si la escala se amplia de tal
manera que se pueda observar el escurrimiento producido por una sola
tormenta, se obtendrá una gráfica como la que se muestra en la figura 4.2.
En este caso el significado de las variables es: t0 es el tiempo de inicio del
escurrimiento directo; tp es el tiempo pico y se define como el tiempo que
transcurre entre el inicio del escurrimiento y el gasto máximo o pico; tb es el
tiempo base y equivale al lapso de tiempo durante el cual ocurre el
escurrimiento directo; y Qb es el gasto base.
122
Los elementos fundamentales del hidrograma son: el gasto antecedente; la
rama ascendente; la cresta o pico; la rama descendente; la curva de
recesión; y el gasto base. A continuación se define cada uno de ellos.
123
El gasto base. Es el valor casi constante al cual tiende la curva de recesión.
Generalmente se asocia con el aporte del agua subterránea; sin embargo,
como se mencionó anteriormente, si la cuenca ha sido alterada en su
funcionamiento hidrológico por alguna obra de control, entonces el flujo base
tendrá su origen en la operación de los sistemas hidráulicos existentes.
Una vez que se define la escala de tiempo para el análisis (gastos horarios,
diarios, mensuales, etc.), se procede a identificar el punto donde se inicia el
gasto antecedente y se traza una línea horizontal.
124
Figura 4.3. Métodos de separación del flujo base y del escurrimiento directo
126
Por su parte, el punto donde se inicia la curva de recesión se localiza al
superponer la curva de recesión del gasto base, dibujada en escala
aritmética, y la rama descendente del hidrograma de análisis. El punto de
inicio de la curva de recesión se localiza donde ambas líneas se separan. En
la fase final, se traza una línea horizontal que comienza en el punto donde
comienza el gasto antecedente y termina en el punto donde se inicia la curva
de recesión. Con este proceso se separa el escurrimiento base del directo y
en la figura 4.3 se indica la aplicación del método.
Cabe hacer mención que los métodos tres presentados también se pueden
analizar sobre una escala semilogarítmica (logaritmo aplicado en el gasto).
−kt
Q t = Q0 e
(4.2)
t
Q t = Q0 k
(4.3)
127
4.5. Aforo de corrientes
Q=VA
(4.4)
128
La medición del nivel se puede hacer a través de una escala, de un
limnímetro, y si se trata de un registro continuo, de un limnígrafo. En la
figura 4.6 se muestra un esquema representativo de una sección de río y de
una sección de canal con sección trapecial.
a) b)
129
(V0.2 + V0.8 )
V =
2
(4.5)
- El valor de velocidad máxima se ubica por arriba del 85% del tirante.
130
http://estudiantesingcivil.blogspot.mx/
Esta velocidad angular se traduce a velocidad del agua (en m/s), utilizando
una fórmula de calibración que previamente se ha determinado para cada
instrumento en particular.
132
En la figura 4.10 se muestra la forma que presenta la curva de calibración
del molinete tipo Price.
V = Cv 2 g ∆h (ρ m − ρ) / ρ
(4.6)
134
líquido del manómetro; ρ es la densidad del agua; y ∆h es la diferencia de
presiones registrada entre las dos ramas del manómetro.
c) Método de flotadores
V = 0.85 Vmedida
(4.7)
135
Si la sustancia se inyecta a un ritmo constante, entonces el gasto se calcula
por medio de la expresión
[
Q T C T − C eq ]
Q=
[C eq − Cb ]
(4.8)
VT [C T − C b ]
Q= ∞
∫ 0 [C d (t) − C b ] dt
(4.9)
a) Método de la velocidad-sección
136
En la figura 4.11 se puede observar un esquema representativo de la
obtención de la sección. Una vez que se conoce la geometría de la sección,
se definen los puntos donde se medirá la velocidad, tanto a lo ancho del
cauce, como a profundidad. Esto se logra por medio del vadeo, en un puente
o en una canastilla.
Con los datos obtenidos se procede a estimar el valor del gasto y para ello se
hace uso de la expresión siguiente (Dingman, 1994):
n
Q= ∑ [Xi +1 − Xi ] [
/ 2 Vi Yi + Vi +1 Yi +1 ]
i =1
(4.10)
137
al 40% del tirante); y Y es la profundidad en cada uno de los puntos de
medición.
2 2 32
[ ] [ h /(h + w)]
2
12
Q = 2 /3 2 g µ 1 + 3 / 2κ (b / B) bh
(4.11)
52
[ ]
12
Q = 8 / 15 2 g tan (θ / 2) µ h
(4.12)
32
[ ]
Q = 2 /3 2 g
12
0.63 b h
(4.13)
138
a). Rectangular b). Triangular
52
Q=µ φD
(4.14)
donde Q es el gasto; y para el resto de las variables ver figura 4.13 (Sotelo,
[ ]
1995). La magnitud de la variable µ = 0.555 + D / 110 h + 0.041 h / D ,
mientras que el valor de φ se obtiene de tablas (Sotelo, 1995).
Los vertedores de pared gruesa son menos precisos que los de pared
delgada en lo que respecta a la medición del flujo; sin embargo, los primeros
son más versátiles, ya que se pueden aplicar en usos múltiples de forma
simultánea y bajo situaciones diversas.
139
La ecuación general del gasto para un vertedor de pared gruesa es (Bos,
1989):
12 32
2
Q= C d C v [2g / 3] b h
3
(4.15)
C d = 0.93 + 0.1 h / L
(4.16)
a). Perspectiva
140
http://estudiantesingcivil.blogspot.mx/
b). Perfil
Para el caso de los vertedores, la función que asocia a las dos variables es la
ecuación general que se muestra a continuación:
32
Q=CBh
(4.18)
Solución:
Figura 4.14. Separación del gasto directo y base con el método de la línea recta
142
b) Método del tiempo fijo. Se traza una línea horizontal desde que inicia el
escurrimiento directo (punto A), hasta el tiempo pico (punto C), tal como
se observa en la figura 4.15.
0.20
N = 0.828 (731) = 3.09 días
N = 3.09 día s
(4.20)
Finalmente, se traza una línea recta que une los puntos C y E. Con este
proceso se separa el gasto directo y base con el método del tiempo fijo.
143
Figura 4.15. Separación del gasto directo y base con el método del tiempo fijo
144
Tabla 4.1. Registro de aforos en la estación hidrométrica Armería, ubicada en
Colima
3
Qi, m /s Elevac ion (h), metros
162.0 0.69
149.9 0.70
197.1 0.96
214.7 1.05
271.4 1.36
283.5 1.50
388.8 2.00
396.9 2.04
716.9 3.31
831.6 3.55
1,265.0 5.13
1,367.6 5.51
1,470.2 5.62
2,106.0 7.53
2,207.3 7.83
2,301.8 8.10
Solución:
Ln h = Ln k + n Ln Q
(4.22)
145
y = a + bx
(4.23)
donde: y = Ln h, a = Ln k, b = n, x = Ln Q
k = 0.0093
n = 0.8817
r = 0.996
0.8817
h = 0.0093 Q
(4.24)
146
Figura 4.16. Curva elevaciones-gastos
Ejemplo 4.3. Gasto que pasa por una sección transversal de una corriente.
Calcular el gasto que pasa por la sección transversal de una corriente que ha
sido aforada por el método del molinete. La figura 4.18 indica la sección
transversal de la corriente, la cual ha sido dividida para su análisis en nueve
secciones denominadas dovelas.
147
1 2 3 4 5 6 7 8 9 10
A B C D E F G H I
5m
1 0.00 0.00
2 0.17 0.35
3 0.33 0.48
4 0.38 0.68
5 0.47 0.91
6 0.41 0.73
7 0.32 0.61
8 0.24 0.33
9 0.18 0.21
10 0.00 0.00
Solución:
n
Q= ∑ [Xi +1 − Xi ] [
/ 2 Vi Yi + Vi +1 Yi +1 ]
i =1
148
Dovela A:
(5 − 0)
[ ]
3
qA = 0 (0) + 0.17 (0.35) = 0.1488 m / s
2
Dovela B:
(10 − 5)
[ ]
3
qB = 0.17 (0.35) + 0.33 (0.48) = 0.5448 m / s
2
Dovela C:
(15 − 10)
[ ]
3
qC = 0.33 (0.48) + 0.38 (0.68) = 1.0420 m / s
2
Dovela D:
(20 − 15)
[ ]
3
qD = 0.38 (0.68) + 0.47 (0.91) = 1.7153 m / s
2
Dovela E:
(25 − 20)
[ ]
3
qE = 0.47 (0.91) + 0.41 (0.73) = 1.8175 m / s
2
Dovela F:
(30 − 25)
[ ]
3
qF = 0.41 (0.73) + 0.32 (0.61) = 1.2363 m / s
2
Dovela G:
(35 − 30)
[ ]
3
qG = 0.32 (0.61) + 0.24 (0.33) = 0.6860 m / s
2
Dovela H:
(40 − 35)
[ ]
3
qH = 0.24 (0.33) + 0.18 (0.21) = 0.2925 m / s
2
Dovela I:
(45 − 40)
[ ]
3
qI = 0.18 (0.21) + 0 (0) = 0.0945 m / s
2
149
Una vez calculado el gasto q para cada dovela, el gasto total que pasa por
la sección transversal de la corriente es:
3
Q = q A + qB + qC + qD + qE + qF + qG + qH + qI = 7.5777 m / s
3
Q = 7.5777 m / s
Solución:
150
http://estudiantesingcivil.blogspot.mx/
VBASE = ∆t QB
donde ∆t = 14 h; QB = 40 m3/s
6 3
VBASE = 14 (3600) (40) = 2.016 x10 m
6 3
VBASE = 2.016 x10 m
151
c). La magnitud del volumen de escurrimiento directo (VED), es igual al área
bajo la curva definida por la figura 4.19. Para el caso específico de un ∆t
constante, se determina con los valores de las figura 4.19 y la expresión
siguiente:
6 3
VED = ∆t Σ Qi = 2 (3600) [40 + 180 + 260 + 160 + 80 + 20] = 5.328 x10 m
6 3
VED = 5.328 x10 m
152
Referencias
153
5. Infiltración
El terreno puede estar formado por un medio poroso (como son los suelos),
por roca consolidada fracturada, o por una combinación de ambos. Las
características particulares de cada uno de los medios que constituyen un
terreno, influye en forma directa sobre el flujo del agua.
a) Medio poroso
154
Cabe mencionar que la condición hidráulica y las propiedades físicas del
medio poroso son fundamentales para su caracterización. A continuación se
presentan algunos de los parámetros físicos de mayor uso:
MS
ρS =
VS
(5.1)
MS
ρm =
VT
(5.2)
Vaire + Vagua ρm
φ= =1−
Vt ρs
(5.3)
155
- Porosidad efectiva. Es el porcentaje de los poros que están conectados
entre sí. Este parámetro se relaciona con el agua que drena
gravitacionalmente o bajo la influencia de una fuerza centrífuga.
Vagua Mt humedo − Mt sec o
θ= =
Vt
(Vt ) ρ agua
5.4)
θ
S=
φ
(5.5)
156
- Tensión mátrica. Cuando el medio es no saturado, entonces la presión del
medio es menor que la atmosférica, por lo que su valor es negativo, es decir,
es una tensión o succión.
157
1000000
100000
Tensión mátrica, en cm
10000
Saturación
1000
Secado
100
10
1
0 10 20 30 40 50
Contenido de humedad, en %
1
Conductividad hidráulica relativa
0.1
0.01
Saturación
0.001
Secado
0.0001
0.00001
0.000001
20 30 40 50
Contenido de humedad, en %
158
de saturación recibe el nombre de punto de entrada del aire, y corresponde
al vaciado de los espacios entre partículas sólidas.
Las curvas que unen las dos condiciones extremas, contenidos de humedad a
saturación y residuales, se llaman curvas principales de secado o de
saturación, según sea el caso. Si en cualquier momento, la condición cambia
en cualquier punto de las curvas principales fuera de los extremos, entonces
se formarán las curvas parciales que tenderán a la curva principal opuesta.
En lo que respecta a la curva de conductividad hidráulica, se tiene un
comportamiento descendente con el aumento del valor absoluto de la tensión
mátrica.
θ − θ
r
Se =
θ − θ
s r
(5.7)
b) Medio fracturado
159
c) Combinación de diferentes medios
160
http://estudiantesingcivil.blogspot.mx/
161
S u p e r fic ie d e l te r r e n o
Z ona vadosa
Fr a n ja c a p ila r N iv e l
fr e á tic o
Z o n a s a tu r a d a
Figura 5.2. Representación esquemática del perfil hídrico del medio poroso
1. Mantener en un área muy reducida una lámina de agua fija definida sobre
la superficie del terreno durante un intervalo de tiempo;
Una parte del cilindro se entierra lo suficiente para evitar que el agua que se
vierta en su interior alcance la superficie del terreno. La parte que queda por
encima de la superficie debe de tener una altura adecuada para evitar que el
agua del interior se derrame.
Por otro lado, si se emplea la variación de una lámina inicial (volumen inicial
conocido) en el tiempo, se debe de contar con un dispositivo adecuado para
ello, como el mostrado en la figura 5.3 b).
163
a). Lámina constante
165
entre 1 y 2 cm. Para emplear los bloques se requiere de la obtención de una
curva de calibración a través del proceso descrito a continuación: en un
recipiente abierto que sea capaz de contener una muestra de suelo de al
menos 2 cm de espesor que contenga al bloque, se satura la muestra y se
mide la resistencia eléctrica en el otro extremo del cable por medio de un
voltímetro; si se permite el cambio del contenido de humedad a través de la
evaporación o el drenado, se podrán realizar diferentes lecturas de
resistencia para cada contenido de humedad.
A − B
θ = 100
B
(5.9)
Aquí se presenta el método del drenaje (referido por Hillel, 1980) para
determinar la relación de conductividad hidráulica y de la difusividad
hidráulica con respecto a la tensión mátrica y el contenido de humedad.
∂θ
∂t
=
∂
∂z
[
K(θ) ∂H ∂z ]
(5.10)
166
L ∂θm ∂H
=K
∂t ∂z L
(5.11)
dθ ∂ψ ∂H
L = −K
dψ ∂t ∂z L
(5.12)
dψ ∂ψ dH
D=K = L /
dθ ∂t dz
(5.13)
∂θ
∂t
[
= −∇ K(Ψ ) ∇{Ψ − z} ]
(5.14)
167
donde ∇ es el operador nabla.
a) Ecuación de Horton
b) Ecuación de Phillip
12
F(t) = S t + K t
(5.16)
dF(t)
f(t) =
dt
(5.17)
168
−1 2
1
f(t) = S t + K
2
(5.18)
12
F(t) = S t
(5.19)
F(t)
F(t) = K t + Ψ ∆θ ln 1 +
Ψ ∆θ
(5.20)
Ψ ∆θ
f(t) = K + 1
F(t)
(5.21)
169
Se presentan dos casos extremos: el primero, cuando la entrada es menor
que la capacidad máxima del suelo de transmitir el agua, es decir, a su
conductividad hidráulica; y el segundo cuando el ritmo de entrada es mayor
que la conductividad hidráulica.
f(t) = 0; para t ≥ t w
C o n te n id o d e h u m e d a d
S u p e r fic ie
d e l te r r e n o R e s id u a l S a tu r a d o
t1
t2
P r o fu n d id a d
t3
Fr e n te h ú m e d o
170
http://estudiantesingcivil.blogspot.mx/
Para estimar el tiempo que transcurre desde el inicio del evento w hasta que
la superficie se satura, se hace uso de las ecuaciones siguientes (Dingman,
1994):
F(t s ) = w t s
(5.22)
F(t s ) = z f (t s ) [φ − θ 0 ]
(5.23)
z f (t s ) [φ − θ 0 ]
ts =
w
(5.24)
K s [Ψf − Ψs ]
w = Ks −
z f (t s )
(5.25)
Ψ∞
Ψf = ∫ 0
K( Ψ ) / K s d Ψ ; con Ψ ∞ representando valores muy grandes;
Ψf = 0.76 Ψf
Ks Ψ
z f (t s ) =
[w − K s ]
(5.26)
K s Ψf φ − θ0 ( )
ts =
[ (
w w − Ks )]
(5.27)
Ψ + y(t)
f(t) = K s 1 + f
z f (t)
(5.28)
172
(
Ψf φ − θ0
)
f(t) = K s 1 +
F(t)
(5.29)
[F(t) − F(t s )] + Ψf (φ − θ0 ) s s (
F(t ) + Ψ φ − θ
0)
t = ∗ Ln
Ks Ks (
F(t) + Ψs φ − θ0)
(5.30)
173
Tabla 5.1. Parámetros representativos de los suelos
174
5.5. Métodos empíricos
Vp = VLL − VED
(5.31)
VED
Pe =
Ac
(5.32)
175
A continuación se presentan algunos de los métodos empíricos más
representativos y de uso más común en la Hidrología.
VED
C=
VLL
(5.33)
176
las características de cada porción del área como un promedio ponderado,
con el apoyo de la expresión:
C1 A 1 + C 2 A 2 + C 3 A 3 ... + C n A n
C=
A 1 + A 2 + A 3 + ... + A n
(5.34)
177
Sin embargo, dado que un gran número de cuencas de nuestro país no están
aforadas y no se dispone de datos, se recurre a métodos estructurados a
partir de la lluvia total y de las características de la cuenca.
508
P − N + 5.08
Pe =
2032
P+ − 20.32
N
(5.36)
N1 A 1 + N2 A 2 + ...Nn A n
N=
Ac
(5.37)
178
Tabla 5.2. Clasificación de suelos
Tipo A Arenas con poco limo y arcilla de tamaño medio (escurrimiento mínimo)
179
Tabla 5.3. Valor del número de escurrimiento N, según uso y tipo de suelo,
condición de la cuenca y pendiente del terreno
Tipo de suelo
Uso de la tierra C ondición de la superficie Pendiente
o cobertura terreno, A B C D
en %
Bosques (sembrados Ralo, baja transpiración ----- 45 66 77 83
y cultivados) Normal, transpiración media ----- 36 60 73 79
Espeso o alta transpiración ----- 25 55 70 77
C aminos De tierra ----- 72 82 87 89
Superficie dura ----- 74 84 90 92
Bosques naturales Muy ralo o baja transpiración ----- 56 75 86 91
Ralo, baja transpiración ----- 46 68 78 84
Normal, transpiración media ----- 36 60 70 76
Espeso o alta transpiración ----- 26 52 62 69
Muy espeso, alta transpiración ----- 15 44 54 61
Descanso (sin cultivo) Surcos rectos ----- 77 86 91 94
C ultivos en surco Surcos rectos >1 72 81 88 91
Surcos rectos <1 67 78 85 89
Surcos en curvas de nivel (contorneo) >1 69 79 84 88
Surcos en curvas de nivel (contorneo) <1 65 75 82 86
Terrazas >1 66 74 80 82
Terrazas <1 62 71 78 81
Leguminosas (sembradas Surcos rectos >1 66 77 85 89
con maquinaria o al voleo) Surcos rectos <1 58 72 81 85
o praderas con rotación Surcos en curvas de nivel >1 64 75 83 85
Surcos en curvas de nivel <1 55 69 78 83
Terrazas >1 63 73 80 83
Terrazas <1 51 67 76 80
Pastizal ------ >1 68 79 86 89
------ <1 39 61 74 80
C ontorneo >1 47 67 81 88
C ontorneo <1 6 65 70 79
Potrero permanente ------ ----- 30 58 71 78
C ereales Surcos rectos >1 65 76 84 88
Surcos rectos <1 63 75 83 87
C ontorneo >1 63 74 82 85
C ontorneo <1 61 73 81 84
Terrazas >1 61 72 79 82
Terrazas <1 59 70 78 81
Superficie impermeable ------ ----- 100 100 100 100
180
http://estudiantesingcivil.blogspot.mx/
N Corrección A Corrección B
10 0.40 2.22
20 0.45 1.85
30 0.50 1.67
40 0.55 1.50
50 0.62 1.40
60 0.67 1.30
70 0.73 1.21
80 0.79 1.14
90 0.87 1.07
100 1.00 1.00
Bosque 55 0.20
Tierra desnuda 10 0.60
Pavimento bituminoso 20 0.85
Campos cultivados 15 0.10
Solución:
C1 A 1 + C 2 A 2 + C 3 A 3 ... + C n A n
C=
A 1 + A 2 + A 3 + ... + A n
C=
[ ]
3.9 (0.20) (0.55) + (0.60) (0.10) + (0.85) (0.20) + (0.10) (0.15)
3.9 (0.55 + 0.10 + 0.20 + 0.15)
C = 0.36
El suelo está formado por arenas finas y limos orgánicos. La longitud del
cauce principal es de 14.4 km con una pendiente media de S=0.0013. Se
tiene una estación climatológica y una de aforo dentro de la cuenca y, con
base en la información recabada en las estaciones, se puede deducir que el
gasto base es nulo. Cuando ocurren tormentas en la zona de captación, la
lluvia antecedente es mayor de 5 cm.
182
Figura 5.7. Curva intensidad-duración-periodo de retorno
Solución:
2
N2 = 60 , A 2 = 0.64 (21.5) = 13.76 km
2
N3 = 100 , A 3 = 0.08 (21.5) = 1.72 km
0.77 0.77
L 14,400
t c = 0.0003245 = 0.0003245 = 6.67 h
S 0.0013
t c = d = 6.67 h
2 2
508 508
P − N + 5.08 40.02 − 78.6 + 5.08
Pe = =
2032 2032
P+ − 20.32 40.02 + − 20.32
N 78.6
Pe = 32.8 cm
184
Figura 5.8. Hietograma e hidrograma
Solución:
3
VED = 126 000 m
185
b) Se calcula la altura de lluvia efectiva total.
Pe = 3.5 mm
Tabla 5.6. Proceso para estimar el valor del índice de infiltración media
ф, P e1 , P e2 , P e3 , P e4 , P e5 , P e6 , ∑P ei ,
mm/h mm mm mm mm mm mm mm
Solución:
2
508
P − N + 5.08
Pe =
2032
P+ − 20.32
N
Pe =
78.6
=
[P − 1.3831]
2032 P + 5.5324
P+ − 20.32
78.6
(5.38)
c) Se procede a graficar la ecuación (5.38) para diferentes valores de Pe y P
y el resultado obtenido se indica en la figura 5.9.
187
Ejemplo 5.5. Calcular el coeficiente de escurrimiento para el ejemplo 5.3.
Solución:
C = 0.19
188
Referencias
189
6. Relaciones lluvia-escurrimiento
190
http://estudiantesingcivil.blogspot.mx/
a) Método de envolventes
191
Figura 6.1. Envolvente mundial de Creager
192
Tabla 6.1. Valores de la envolvente de Creager para gastos máximos
Re g ió n E s t a c ió n C o rrie n t e C
Hid ro ló g ic a
1 C e rc a d e Ne s t o r C a lifo rn ia Río T iju a n a 2 .9 4 7
3 E l O jo d e A g u a Río P u rís ima 1 .1 6 3
8 P it iq u it o I Río La A s u n c ió n 5 .0 0 0
S a n t a T e re s a Río A lt a r 1 2 .0 0 0
9 E l Á g u ila Río Y a q u i 1 6 .0 0 0
E l O ré g a n o Río S o n o ra 6 .0 0 0
Punt o de A gua Río M á t a p e 5 .0 0 0
T re s He rma n a s Río M a y o 3 3 .0 0 0
10 Hu it e s Río F u e rt e 5 8 .0 0 0
Bo d ig u a ra t o Río Bo d ig u a ra t o 5 9 .0 0 0
11 Ba lu a rt e Río Ba lu a rt e 9 9 .0 8 5
A c aponet a Río A c a p o n e t a 1 1 0 .0 0 0
12 Lo s F re s n o s 3 9 Río A n d a má c u a ro 1 4 .0 0 0
P e ñ u e lit o s 2 0 Río d e la E rre 1 4 .0 0 0
Y ago Río S a n t ia g o 1 9 .0 0 0
P a s o d e A n a lc o Río S a n t ia g o 1 3 .0 0 0
La C u ñ a Río V e rd e 6 .0 0 0
La Bo q u illa Río Hu ic ic ila 5 .0 0 0
E l C a imá n Río Bo la ñ o s 5 .0 0 0
Hu a y a n a mo t a II Río Hu a y a n a mo t a 1 1 .0 0 0
T a ra n d a c u a o A rro y o T a ra n d a c u a o 1 0 .3 2 0
13 Ju ma t á n Río In g e n io 3 .0 0 0
E l Re f ilió n Río Hu ic ic ila 5 .0 0 0
P a s o d e A ro c h a Río Hu ic ic ila 2 7 .0 0 0
14 P u e n t e A me c a Río A me c a 1 .0 0 0
La V e g a Río A me c a 1 .0 0 0
P ijin t o Río A me c a 7 .0 0 0
La s Ga v io t a s Río A me c a 9 .0 0 0
Puent es F F CC Río A h u a lu lc o 3 .0 0 0
E l S a lit re Río C o c u la 1 .0 0 0
S a n M a rt ín Hid a lg o Río S a n M a rt ín 5 .0 0 0
C o rrin c h is Río M a s c o t a 7 .0 0 0
La De s e mb o c a d u ra Río M a s c o t a 9 .0 0 0
18 A - 9 E l M o lin o Río Ixt la h u a c a 2 3 .6 0 0
19 T ec pan Río T e c p a n 4 6 .0 0 0
20 E l S a lit re Río O mit lá n 6 4 .0 0 0
22 Ost ut a Río O s t u t a 4 9 .0 0 0
23 S u c h ia t e Río S u c h ia t e 5 3 .0 0 0
24 O rie n t e C e rc a d e l Río A rro y o P in t o 9 1 .8 0 0
S a b in a s Río S a b in a s 1 2 .5 0 0
24 P o n ie n t e E l C u c h illo Río S a n Ju a n 3 7 .8 3 0
25 La E s p e ra n z a Río S a n Ra fa e l o E l C a rriz a l 7 3 .6 0 0
193
Tabla 6.1. Valores de la envolvente de Creager para gastos máximos (continuación)
Re g ió n E s t a c ió n C o rrie n t e C
Hid ro ló g ic a
26 T e mp o a l Río T e mp o a l 4 1 .0 0 0
T e rre rilla s Río C a la b o z o 3 7 .0 0 0
E l P u ja l Río T a mp a ó n 3 2 .0 0 0
T e ma ma t la Río A ma ja c 2 9 .0 0 0
M a g is c a t z in Río Gu a y a le jo 2 3 .0 0 0
Pánuc o Río P á n u c o 2 1 .0 0 0
La s A d ju n t a s Río P á n u c o 2 1 .0 0 0
C . M é xic o - La re d o Río Gu a y a le jo 7 5 .0 0 0
C . T a mp ic o - M a n t e Río Gu a y a le jo 5 4 .0 0 0
Ga lin d o Río Ga lin d o 1 5 .0 0 0
V enados Río M e t z t it lá n 3 7 .0 0 0
T e p e ji Río T e p e ji 7 .0 0 0
T la u t la Río T la u t la 6 .0 0 0
La s Ro s a s Río La s Ro s a s 7 .0 0 0
27 P o z a Ric a Río C a z o n e s 6 1 .0 0 0
A t epat ahua Río A t e p a t a h u a 8 5 .0 0 0
M a rt ín e z d e la T o rre Río Bo b o s 5 3 .0 0 0
Ra u d a l Río M is a n t la 4 7 .0 0 0
28 Ja lc o mu lc o Río La A n t ig u a 3 8 .0 0 0
29 La s P e rla s Río C o a t z a c o a lc o s 3 8 .8 8 0
30 Canasay ab Río C h a mp o t ó n 1 3 .6 0 0
34 La T ra s q u illa Río d e l C a rme n 3 .0 0 0
36 P re s a Lá z a ro C á rd e n a s Río Na z a s 2 8 .0 0 0
37 Lo s P ila re s A rro y o La P a ra d a 0 .9 0 0
E l Grit o A rro y o Gra n d e 6 .2 0 6
E l T u le A rro y o E l T u le 2 .0 0 0
T u la Río T u la 3 .0 0 0
P re s a E l P e a je A rro y o Gra n d e 1 9 .4 0 0
P re s a S a n Jo s é Río S a n t ia g o 1 7 .0 0 0
P re s a A lv a ro O b re g ó n A rro y o A lo q u in e s 4 2 .8 0 0
CL
q= 0.85
(A c + 259)
(6.2)
194
La figura 6.3 muestra para México la envolvente general definida por el
método de Lowry. Para la envolvente mundial se puede tomar un valor de CL
= 3 500. Para las 37 Regiones Hidrológicas de nuestro país, ver tabla 6.2.
195
Tabla 6.2. Valores de la envolvente de Lowry para gastos máximos
Re g ió n E s t a c ió n CL
Hid ro ló g ic a
1 Ba ja C a lif o rn ia No ro e s t e 980
2 Ba ja C a lif o rn ia C e n t ro - O e s t e 530
3 Ba ja C a lif o rn ia S u ro e s t e 2 ,1 9 0
4 Ba ja C a lif o rn ia No re s t e 1 ,0 5 0
5 Ba ja C a lif o rn ia C e n t ro - E s t e 990
6 Ba ja C a lif o rn ia S u re s t e 5 ,1 2 0
7 Río C o lo ra d o 1 ,0 5 0
8 S o n o ra No rt e 760
9 S o n o ra S u r 2 ,1 4 0
10 S in a lo a 3 ,2 9 0
11 P re s id io - S a n P e d ro Z o n a C o s t e ra 4 ,6 3 0
11 P re s id io - S a n P e d ro Z o n a A lt a 470
12 Le rma - S a n t ia g o 1 ,2 9 0
13 Río Hu ic ic ila 760
14 Río A me c a 600
15 C o s t a d e Ja lis c o 5 ,2 7 0
16 A rme ría - C o a h u a y a n a 4 ,9 4 0
17 C o s t a d e M ic h o a c á n 2 ,1 0 0
18 Ba ls a s A lt o 1 ,0 9 0
18 Ba ls a s M e d io y Ba jo 4 ,4 5 0
19 C o s t a Gra n d e d e Gu e rre ro 2 ,1 0 0
196
Tabla 6.2. Valores de la envolvente de Lowry para gastos máximos (continuación)
Re g ió n E s t a c ió n CL
Hid ro ló g ic a
20 C o s t a C h ic a d e Gu e rre ro - Río V e rd e 3 ,1 8 0
20 A lt o Río V e rd e 390
21 C o s t a d e O a xa c a 3 ,0 0 0
22 T e huant e pe c 2 ,1 7 0
23 C o s t a d e C h ia p a s 1 ,1 9 0
24A A lt o Bra v o - C o n c h o s 1 ,0 2 0
24B M e d io Bra v o 5 ,1 7 0
24C Río S a la d o 1 ,4 1 0
24D Ba jo Bra v o 2 ,1 3 0
25 S a n F e rn a n d o - S o t o La M a rin a 2 ,3 3 0
26A A lt o P a n u c o 1 ,3 6 0
26B Ba jo P a n u c o 3 ,0 1 0
26C V a lle d e M é xic o 760
27 T u xp a n - Na u t la 2 ,4 5 0
28 P a p a lo a p a n 1 ,7 5 0
29 C o a t z a c o a lc o s 1 ,8 4 0
30 Grija lv a - Us u ma c in t a 2 ,1 3 0
30 A lt o Grija lv a 610
31 Y uc at án Oest e 370
32 Y u c a t á n No rt e ND
33 Y uc at án Est e ND
34 C u e n c a s C e rra d a s d e l No rt e 230
35 M a p imí ND
36 Na z a s 1 ,5 1 0
36 A guanav al 380
37 E l S a la d o 1 ,3 1 0
197
donde Qp es el gasto pico o máximo, en m3/s; C es el coeficiente de
escurrimiento, adimensional; i es la intensidad de lluvia para una duración
que es igual al tiempo de concentración, en mm/h; A es el área de la
cuenca, en km2.
198
Tabla 6.3. Valores del coeficiente de escurrimiento C
Zonas Comerciales:
Zona comercial 0.70 0.95
Vecindarios 0.50 0.70
Zonas Residenciales:
Unifamiliares 0.30 0.50
Multifamiliares, espaciados 0.40 0.60
Multifamiliares, compactos 0.60 0.75
Semiurbanas 0.25 0.40
Casas habitación 0.50 0.70
Zonas Industriales:
Espaciado 0.50 0.80
Compacto 0.60 0.90
Calles:
Asfaltadas 0.70 0.95
De concreto hidráulico 0.70 0.95
Adoquinadas 0.70 0.85
Praderas:
Suelos arenosos planos (pendientes 0.02 o menos) 0.05 0.10
Suelos arenosos con pendientes medias (0.02-0.07) 0.10 0.15
Suelos arenosos escarpados (0.07 o más) 0.15 0.20
Suelos arcillosos planos (pendientes 0.02 o menos) 0.13 0.17
Suelos arcillosos con pendientes medias (0.02-0.07) 0.18 0.22
Suelos arcillosos escarpados (0.07 o más) 0.25 0.35
199
La intensidad de lluvia (i) se determina con el auxilio de las curvas
intensidad-duración-periodo de retorno (i-d-Tr). La selección de la magnitud
de la intensidad (i) se fundamenta con la estimación de la frecuencia y de la
duración.
L
tc =
3600 v
(6.4)
Las tres hipótesis fundamentales del método del hidrograma unitario son las
siguientes:
200
http://estudiantesingcivil.blogspot.mx/
201
Ahora bien, el hidrograma unitario tradicional tiene dos limitaciones:
202
Figura 6.4. Esquema que indica el cálculo del Hidrograma Unitario Tradicional
203
6.4. Método de la curva o hidrograma S
A
qeq =
3.6 de
(6.6)
204
Figura 6.5. Curva o hidrograma S
3) Se obtienen las ordenadas del hidrograma unitario para una duración d1,
multiplicando los valores del paso anterior por la relación de/d1 (de esta
205
forma se garantiza que el área bajo el hidrograma unitario sea igual al
área de la cuenca por 1 mm).
206
6) Si la duración en exceso de corresponde a la del hidrograma unitario, se
procede a calcular la avenida de diseño, multiplicando las ordenadas del
hidrograma unitario por la lluvia efectiva obtenida en el paso 5. Si no
ocurre esto, hay que usar el método de la curva S para tener el
hidrograma unitario para la duración de interés y luego multiplicar sus
ordenadas por la lluvia efectiva determinada en el punto 5.
(6.7)
b) u (t) = 0 , para t ≤ 0
c) u(t) → 0 , para t → 0
∞
d) ∫ u (t) dt = 1.0
0
∞
e) ∫ u (t) dt = t
0
r
208
donde tr es el tiempo de retraso del hidrograma unitario instantáneo.
209
La hipótesis fundamental del método considera que para un hidrograma
unitario de duración en exceso ∆t, el escurrimiento directo (ver figura 6.9) se
calcula con el sistema de ecuaciones que se describe a continuación:
210
http://estudiantesingcivil.blogspot.mx/
expresión sujeta a Pei = 0 para i > NP; ui = 0 para i > NU, i = 1, 2,..., NQ,
donde:
NU = NQ – NP + 1
(6.9)
La solución del sistema de ecuaciones definido por (6.8) muy pocas veces se
puede obtener, ya que el hidrograma unitario es muy sensible a pequeños
errores de los datos de partida y, en ocasiones, la solución puede dar
resultados negativos.
211
ei = Qm − Qi = Qm − (Pe1 ui + Pe 2 ui −1 + Pe 3 ui − 2 + ... + Pei u1 )
(6.10)
2
S= ∑ ei
(6.11)
NU
φPe Q (τ) = ∑ ui φPePe (τ − i + 1)
i =1
(6.12)
NQ
φPe Q (γ) = ∑ Pel Ql + γ
l =1
(6.13)
NP
φPePe (ρ) = ∑ Pek Pek + P
k =1
(6.14)
212
6.6. Hidrogramas unitarios sintéticos
A
qp = 0.208
tp
(6.15)
213
Figura 6.10. Hidrograma Unitario Triangular
Los tiempos pico (tp) y de recesión (tr), se evalúan con las expresiones:
tp = t c + 0.6 t c
(6.16)
t r = 1.67 t p
(6.17)
A
qp = 0.208
tp
(6.18)
tp = t c + 0.6 t c
(6.19)
215
a). Se escoge un valor de t/tp y con el auxilio de la figura 6.11 se obtiene
q/qp.
c) Método de I-Pai-Wu
n −1
1 t −t /K
Q (t) = e
K(n − 1)! K
(6.20)
A Pe
Qp = 0.278 f (n, t p)
tp
(6.21)
n 1−n
(n − 1) e
f(n, t p ) =
Γ (n)
(6.22)
216
donde n es una constante que depende de las características fisiográficas de
la cuenca en estudio; tp es el tiempo pico; Γ (n) es la función matemática
gamma.
4 tp
n=
K
(6.23)
Por su parte, para estimar los valores de tp y K, I-Pai-Wu propone utilizar las
dos ecuaciones siguientes:
5) Con t/tp del hidrograma sintético de la figura 6.12, se calcula Q/Qp o bien
se determina con la expresión:
217
t
n −1 ( −1)
− (n −1) tp
Q t
= e
Qp tp
(6.26)
218
Figura 6.12. Hidrograma adimensional propuesto por I-Pai-Wu
Pe = Peb
(6.27)
Pe = F Peb
(6.28)
219
Asimismo, es importante recalcar que este método no proporciona un
hidrograma unitario, sino el de escurrimiento directo, ya que el cálculo de Qp,
con la ecuación (6.21), ya considera a la lluvia efectiva.
Qp = qp Pe
(6.30)
A
qe =
3.6 d
(6.31)
qp
Z=
qe
(6.32)
A Z Pe
Qp = 0.278
d
(6.33)
220
http://estudiantesingcivil.blogspot.mx/
Peb Pe
Qp = A Z
d Peb
(6.34)
Peb Pe
X= ; y Y = 0.278
d Peb
(6.35)
Qp = A X Y Z
(6.36)
221
Figura 6.13. Factor de reducción Z
222
Figura 6.14. Definición del tiempo de retraso
0.64
L
t R = 0.0050
S
(6.37)
223
n
L L ca
t R = 0.75174 (C t )
S
(6.38)
0.3
t R = C t (L L ca )
(6.39)
t R = 0.6 t c
(6.40)
224
5) Se calcula la relación d/tR, para determinar el factor Z con ayuda de la
figura 6.13.
8) El valor más grande de los gastos pico calculado con la ecuación (6.36)
corresponde al gasto de diseño.
225
Figura 6.15. Registro simultáneo de lluvia y escurrimiento
Solución:
226
c) Se calcula la altura de lluvia efectiva (Pe)
6 3
VED 5.328 (10 ) m
Pe = = = 6 mm
Ac 6 2
888 (10 ) m
Pe = (7-4) + (9-4) = 3 + 5 = 8 mm ≠ 6 mm
Pe = (7-5) + (9-5) = 2 + 4 = 6 mm
En consecuencia:
φ = 2.5 mm/h
227
Tabla 6.4. Cálculo de las ordenadas del hidrograma unitario
(1 ) (2 ) (3 ) (4 ) (5 )
T ie mp o , E s c u rrimie n t o E s c u rrimie n t o E s c u rrimie n t o O rd e n a d a s d e l HU
h o ra s t o t a l, m3 /s b a s e , m3 /s d ire c t o (2 )- (3 ) (4 )/P e , m3 /s /mm
0 40 40 0 0 .0 0
2 80 40 40 6 .6 7
4 220 40 180 3 0 .0 0
6 300 40 260 4 3 .3 3
8 200 40 160 2 6 .6 7
10 120 40 80 1 3 .3 3
12 60 40 20 3 .3 3
14 40 40 0 0 .0 0
T otal 740 1 2 3 .3 3
f). Con los valores de la tabla 6.4 se grafica el hidrograma unitario tal como
aparece en la figura 6.16
de = 4 horas
228
he = 1 mm
h). Comprobación
VED
Pe =
Ac
3
VED 123.33 (2) (3600) m
Ac = =
he −3
10 m
2 2
A c = 887.9 km ≈ 888 km ; En consecuencia: estimación correcta
Solución:
229
Tabla 6.5. Cálculo del Hidrograma S
(1 ) (2 ) (3 ) (4 ) (5 ) (6 ) (7 ) (8 )
Tie m p o , HU H U d e s p la - H U d e s p la - H U d e s p la - H U d e s p la - H id ro g ra m a S H id ro g ra m a S
h o ra s de = 4 h zado 4 h zado 8 h zado 12 h zado 16 h a ju s ta d o
0 0 .0 0 0 .0 0 0 .0 0
2 6 .6 7 6 .6 7 6 .6 7
4 3 0 .0 0 0 .0 0 3 0 .0 0 3 0 .0 0
6 4 3 .3 3 6 .6 7 5 0 .0 0 5 0 .0 0
8 2 6 .6 7 3 0 .0 0 0 .0 0 5 6 .6 7 5 6 .6 7
10 1 3 .3 3 4 3 .3 3 6 .6 7 6 3 .3 3 6 1 .6 7
12 3 .3 3 2 6 .6 7 3 0 .0 0 0 .0 0 6 0 .0 0 6 1 .6 7
14 0 .0 0 1 3 .3 3 4 3 .3 3 6 .6 7 6 3 .3 3 6 1 .6 7
16 3 .3 3 2 6 .6 7 3 0 .0 0 0 .0 0 6 0 .0 0 6 1 .6 7
18 0 .0 0 1 3 .3 3 4 3 .3 3 6 .6 7 6 3 .3 3 6 1 .6 7
20 3 .3 3 2 6 .6 7 3 0 .0 0 6 0 .0 0 6 1 .6 7
Ac 888
qeq = = = 61.67
3.6 de 3.6 (4)
230
http://estudiantesingcivil.blogspot.mx/
(1 ) (2 ) (3 ) (4 ) (5 ) (6 )
Tie m p o , HU H id ro g ra m a S H id ro g ra m a S (3 ) - (4 ) H U p a ra d e = 2 h
h o ra s de = 4 h a ju s t a d o D e s p la z a d o 2 h (5 )x (4 /2 )
0 0 .0 0 0 .0 0 0 .0 0 0 .0 0
2 6 .6 7 6 .6 7 0 .0 0 6 .6 7 1 3 .3 3
4 3 0 .0 0 3 0 .0 0 6 .6 7 2 3 .3 3 4 6 .6 6
6 4 3 .3 3 5 0 .0 0 3 0 .0 0 2 0 .0 0 4 0 .0 0
8 2 6 .6 7 5 6 .6 7 5 0 .0 0 6 .6 7 1 3 .3 4
10 1 3 .3 3 6 1 .6 7 5 6 .6 7 5 .0 0 1 0 .0 0
12 3 .3 3 6 1 .6 7 6 1 .6 7 0 .0 0 0 .0 0
14 0 .0 0 6 1 .6 7 6 1 .6 7 0 .0 0 0 .0 0
16 6 1 .6 7 6 1 .6 7 0 .0 0 0 .0 0
To ta l 1 2 3 .3 3 1 2 3 .3 3
e). La figura 6.18 presenta el hidrograma unitario válido para una duración
efectiva de 2 horas y una altura de lluvia efectiva de 1 mm.
231
Figura 6.18. Hidrograma unitario válido para una de = 2 h y una he = 1 mm
Solución:
232
Figura 6.19. Hietograma de lluvia en exceso e hidrograma de escurrimiento directo
NU = NQ − NP + 1 ; donde: NQ = 6 y NP = 3
NU = 6 − 3 + 1 = 4
233
NU
φPQ (τ) = ∑u
i =1
i φPP (τ − i + 1) , τ = 0, 1,...., (NU − 1)
τ=0 φPQ (0) = u1 φPP (0) + u2 φPP (−1) + u3 φPP (−2) + u4 φPP (−3)
τ =1 φPQ (1) = u1 φPP (1) + u2 φPP (0) + u3 φPP (−1) + u 4 φPP (−2)
τ=2 φPQ (2) = u1 φPP (2) + u2 φPP (1) + u3 φPP (0) + u 4 φPP (−1)
τ=3 φPQ (3) = u1 φPP (3) + u2 φPP (2) + u3 φPP (1) + u4 φPP (0)
NQ
φPQ (γ) = ∑P
l =1
l Ql + γ , γ = 0, 1,...., (NU − 1)
l > NP, Pl = 0
NP
φPP (ρ) = ∑P
k =1
k Pk + ρ , ρ = 0, 1,...., (NP − 1)
234
14 u1 + 9 u2 + 2 u3 + 0 u 4 = 89
9 u1 + 14 u2 + 9 u3 + 2 u 4 = 127
2 u1 + 9 u2 + 14 u3 + 9 u 4 = 105
0 u1 + 2 u2 + 9 u3 + 14 u 4 = 54
(6.41)
u1 = 2.33
u2 = 5.67
u3 = 2.61
u 4 = 1.33
235
Calcular el hidrograma de escurrimiento directo con el hidrograma unitario
instantáneo (HUI) obtenido en el ejemplo 6.3, para el hietograma de lluvia
en exceso mostrado en la figura 6.21.
Solución:
a). Los datos de la lluvia en exceso (Pe) y ordenadas del hidrograma unitario
instantáneo son:
Pe = 2, 6, 3, 1
ui = 2.33, 5.67, 2.61, 1.33
NQ = NU + NP − 1 , donde NU = 4, NP = 4
NQ = 4 + 4 − 1 = 7
236
c). Se calculan las ordenadas del hidrograma de escurrimiento directo con la
ecuación siguiente:
i
Qi = ∑ uk Pi − k +1
k =1
Q1 = P1 u1 = 2 (2.33) = 4.67
Q2 = P1 u2 + P2 u1 = 2 (5.67) + (6)(2.33) = 25.33
Q3 = P1 u3 + P2 u2 + P3 u1 = 2 (2.61) + 6 (5.67) + 3 (2.33) = 46.33
Q 4 = P1 u4 + P2 u3 + P3 u2 + P4 u1 = 2 (1.33) + 6 (2.61) + 3 (5.67) + 1 (2.33) = 38.00
Q5 = P2 u4 + P3 u3 + P4 u2 = 6 (1.33) + 3 (2.61) + 1 (5.67) = 21.67
Q 6 = P3 u4 + P4 u3 = 3 (1.33) + 1 (2.61) = 6.67
Q7 = P4 u 4 = 1 (1.33) = 1.33
237
Ejemplo 6.5. Hidrograma Unitario Triangular (HUT)
Solución:
A 0.208 (160) 3
qp = 0.208 = = 7.28 m / s / mm
tp 4.57
238
Figura 6.23. Hidrograma unitario triangular
Las características de la cuenca de estudio son: área 215 km2; longitud del
cauce principal 20.92 km; y pendiente media del cauce 0.006.
Solución:
Ac (215)
qp = 0.208 = 0.208 = 8.62 m3 / s / mm
tp 5.19
b). Con ayuda de la figura 6.11 y la relación t/tp, se obtiene q/qp tal como se
indica en la columna (3) de la tabla 6.7.
239
Tabla 6.7. Hidrograma de escurrimiento directo según el método del SCS
240
http://estudiantesingcivil.blogspot.mx/
Solución:
0.356
259.923 Tr
i=
0.558
d
0.356
259.923 (25)
i= = 61.966 mm / h
0.558
(102)
Pb
i=
d
241
2 2
508 508
P − + 5 .08 10.534 − + 5 .08
b N 93
Peb = = = 8.54 cm
2032 2032
Pb + − 20.32 10.534 + − 20.32
N 93
252.0
F= = 1.534
164.3
t p = 0.93 (A c )1.085 (L) −1.233 (S) −0.668 = 0.93 (60)1.085 (15) −1.233 (4.94) −0.668 = 0.963 h
K = 0.73 (A c )0.937 (L) −1.474 (S) −1.473 = 0.73 (60)0.937 (15) −1.474 (4.94) −1.473 = 0.05935
g) Se estima el valor de n:
4 tp 4 (0.963)
n= = = 64.9
K 0.05935
242
n 1−n
(n − 1) e
f (n, t p ) =
Γ (n)
∞
n −1 −x
Γ (n) = ∫ x e dx
0
−n (n − 1 / 2) 1 1 139 571
Γ (n) = e n 2π 1 + + − −
12 n 288 n2 51840 n3 2488320 n 4
donde 1 ≤ m ≤ 2
243
0.278 A c Pe 0.278 (60) (131) 3
Qp = f (n, t p ) = (3.20) = 7260.9 m / s
tp 0.963
n −1 (t / tp −1)
− (n −1)
Q t
= e
Qp tp
Q
= 0.090
Qp
(6.42)
Q p = 7260.9 m3 / s
k). La primera pareja de datos que asocian tiempo y gasto es (0.722 horas,
653.521 m3/s). El proceso se repite tantas veces como sea necesario para
determinar la forma del hidrograma de escurrimiento o avenida de diseño.
244
La tabla 6.8 indica en forma resumida los resultados obtenidos para 10
parejas de datos, mientras que la figura 6.25 muestra la avenida de
diseño obtenida a partir del método de I-Pai-Wu.
3
t/t p t, en h Q/Qp Q, en m /s
245
Ejemplo 6.8. Método de Chow
246
Figura 6.26. Curva i-d-Tr para un periodo de retorno de 10 años
Solución:
a). Con los datos del tipo y uso del suelo, mostrados en las tablas 5.2 y 5.3,
se obtuvieron los valores de N:
c). El número de escurrimiento obtenido con este proceso se debe afectar por
un factor que considere la lluvia antecedente. Es decir, si hubo lluvia cinco
247
días antes se hace una corrección de acuerdo con la cantidad de lluvia que
se haya registrado. Esto se hace con el auxilio de la tabla 5.4, la cual
indica el tipo de corrección A y B, según la magnitud de la lluvia
antecedente.
P = 6 (6.67) = 40.02 cm
P = Pb
Pe = Peb
Peb 328
X= = = 49.2
d 6.67
Pb
Y = 0.278 = 0.278
PEB
248
0.64
10000
tR = 0.0050 = 1.67 h
1.3
t R = 1.67 h
l). La tabla 6.9 señala los resultados obtenidos para otras duraciones de la
tormenta.
Tabla 6.9. Resumen de los cálculos realizados para obtener el gasto de diseño
3
QD = 480.3 m / s
249
Referencias
Kuichling, E. (1989). The relation between the rainfall and the discharge of
sewers in populous districts. Transactions ASCE, Vol. 20.
250
http://estudiantesingcivil.blogspot.mx/
7. Tránsito de avenidas
7.1. Introducción
Ahora bien, para predecir las variaciones temporal y espacial de una onda de
avenida a través de un tramo de río o en un vaso de almacenamiento, o bien
para determinar el hidrograma de salida de una cuenca sobre las que
presentó una determinada lluvia, se usan los procedimientos de tránsito de
avenidas.
252
- Volumen de azolves. Es el que queda abajo del nivel de la toma y se
reserva para recibir el acarreo de sólidos por el río durante la vida útil de la
presa. Es conveniente señalar que el depósito de sedimentos en una presa
no se produce con un nivel horizontal como está mostrado en la figura 7.1,
sino que los sedimentos se reparten a lo largo del embalse, teniéndose los
más gruesos al principio del mismo y los más finos cerca de la cortina. De
hecho, en algunos casos existe movimiento de los sedimentos depositados
dentro del vaso, fenómeno que se conoce como corriente de densidad.
253
En síntesis, se analizarán los métodos para evaluar el volumen útil que debe
tener una presa para satisfacer sus demandas, así como el volumen de
superalmacenamiento requerido para que la presa no presente situaciones
de peligro.
dV
I−O =
dt
(7.1)
Ii + Ii + 1 Oi + Oi + 1 V − Vi
− = i +1
2 2 Δt
(7.2)
254
donde Δt es el intervalo de tiempo, en segundos; i, i+1 son subíndices que
representan los valores de las variables en el tiempo t = i Δt y t = (i +1) Δt,
respectivamente.
⎛ 2 Vi ⎞ 2 Vi + 1
Ii + Ii + 1 + ⎜⎜ − Oi ⎟⎟ = + Oi + 1
⎝ Δt ⎠ Δt
(7.3)
Δt ≤ 0.01 t p
(7.4)
255
40
30
Elevaciones, en m 20
10
0
0 5000 10000 15000 20000 25000
3
Capacidades, en miles de m
Además, es útil establecer una función entre la elevación del agua y el gasto
de salida, lo cual se logra por medio de la curva elevaciones-gastos de salida.
Esta curva se calcula por medio de la ecuación de un vertedor rectangular, la
cual es del tipo siguiente:
3/2
O=CL h
(7.5)
256
La curva elevaciones-gastos de salida puede ser dibujada a partir de la
ecuación (7.5), tomando en cuenta que la elevación del agua es igual a la
suma de la elevación de la cresta del vertedor más la carga de agua. La
figura 7.3 muestra una curva de este tipo.
40
38
Elevaciones, en m
36
34
32
30
0 500 1000 1500 2000
3
Gasto de salida, en m /s
257
Por su parte, para resolver la ecuación (7.3), tomando en cuenta las curvas
elevaciones-capacidades y elevaciones–gastos de salida, se pueden utilizar
tres métodos: semigráfico, numérico y algoritmo de incrementos finitos.
a) Primera etapa
258
6
0
0 20 40 60 80
3 3
2V/∆t + O, en m /s (10 )
b) Segunda etapa
2Vi +1
+ Oi +1 = K
Δt
Δt
Vi +1 = (K − Oi +1 )
2
3) Se considera a i+1 como i, es decir Ii+1, Oi+1 y Vi+1 pasan a ser Ii, Oi y Vi
en el instante siguiente. Se repite el procedimiento hasta terminar todos
los valores de interés del hidrograma de entrada.
259
Tabla 7.1. Secuencia de cálculo para el método semigráfico
Columna 1
Columna 2
Columna 3
Columna 4
260
http://estudiantesingcivil.blogspot.mx/
Columna 5
Columna 7
⎛ 2V ⎞ 2Vi + 1
Ii + Ii + 1 + ⎜⎜ i − Oi ⎟⎟ = + Oi + 1
⎝ Δt ⎠ Δt
1) Se conocen los valores Ii, Ii+1, Vi, Oi y Δt. Por consecuencia se puede
calcular el miembro izquierdo de la ecuación (7.3).
3 /2
Q = C L (h − H)
(7.6)
N
V =Kh
(7.7)
262
donde V es el volumen de almacenamiento; h es la elevación o nivel del
agua; K y N son parámetros que se determinan al ajustar la ecuación
(7.7) a los datos disponibles.
h
H
Figura 7.5. Elevación de la cresta del vertedor (H) y nivel del agua (h)
263
Figura 7.6. Hidrograma de entrada y salida
Δt ⎛ ⎞ Δt ⎛I ⎞
ΔV = ⎜ Ii − Qi ⎟ + ⎜ i + 1 − Qi + 1 ⎟
2 ⎝ ⎠ 2 ⎝ ⎠
(7.8)
ΔV Ii + 1 + Ii Qi + 1 + Qi
= −
Δt 2 2
(7.9)
dV ≅ ΔV
(7.10)
dh ≅ hi +1 − hi
(7.11)
264
La variación de la curva elevaciones-capacidades (volúmenes almacenados)
es:
N
V = K hi
N −1
dV = K N hi dh
(7.12)
dQ ≅ ΔQ ; ΔQ = Q i + 1 − Q i ; ΔQ = Q i + 1 − Q i − Q i + Q i
Qi +1 + Q1 ΔQ
ΔQ + 2Qi = Qi +1 + Qi ; = + Qi
2 2
ΔQ Q i + 1 + Qi
= − Qi
2 2
(7.14)
3/2 1/2
3
Q = C (hi − H) ; dQ = C (hi − H) dh
2
1/2
dQ 3 dQ ΔQ
= C (hi − H) dh ; si =
2 4 2 2
1/2
ΔQ 3
= C (hi − H) dh
2 4
(7.15)
1/2
Qi +1 + Q1 3
− Qi = C (hi − H) (hi +1 − hi)
2 4
265
1/2
Qi +1 + Q1 3
= C (hi − H) (hi +1 − hi ) + Qi
2 4
(7.16)
N −1
1/2
KNhi ⎛h ⎞ Ii +1 + Ii − 3 C (h − H) ⎛h ⎞
⎜ i +1 − hi ⎟ = i ⎜ i +1 − hi ⎟ − Qi
Δt ⎝ ⎠ 2 4 ⎝ ⎠
(7.17)
N 1/2
KN Ii +1 + Ii 3
hi + + C (hi − H) hi − Qi
hi +1 = Δt 2 4
N −1 1/2
KN 3
hi + C (hi − H)
Δt 4
(7.18)
Ii +1 + Ii
a hi + − Qi
hi +1 = 2
a
(7.19)
Las expresiones que permiten definir los valores de la variable auxiliar (a) y
de los gastos de descarga ( Q i ) son:
N −1 1/2
KN 3
a= hi + C (hi − H)
Δt 4
(7.20)
266
3 /2
Qi = C L (hi − H) ; si hi ≥ H
(7.21)
Qi = 0 ; si hi < H
(7.22)
- Curva elevaciones-capacidades
N
V =K h
(7.23)
X = a+b Y
(7.25)
Por otra parte, para estimar los valores óptimos de a y b , se establece que
la suma de los cuadrados de los errores, diferencia entre los valores
observados y estimados, sea mínima. Al realizar este proceso se obtiene un
sistema de ecuaciones normales representado por:
267
ΣX = a N + b ΣY
(7.26)
2
ΣXY = a ΣY + b ΣY
2
(ΣX)(ΣY ) − Σ(Y)(ΣXY )
a=
2 2
NΣY − (ΣY)
(7.27)
NΣXY − Σ(Y)(ΣY)
b=
2 2
N ΣY − (ΣY)
N=b
(7.28)
a
K=e
(7.29)
268
⎡ ⎤
V = K ⎢XI + (1 − X) O⎥
⎣ ⎦
(7.30)
[
ΔV = Vi + 1 − Vi = K X (Ii + 1 − Ii ) + (1 − X) (Oi + 1 − Oi ]
(7.31)
donde t i = i Δt = t1 ; y t i + 1 = (i + 1) Δt = t 2
269
1 ⎛ ⎞ 1 ⎛ ⎞
ΔV = V2 − V1 = ⎜ I1 + I2 ⎟ Δt − ⎜ O1 + O2 ⎟ Δt
2 ⎝ ⎠ 2 ⎝ ⎠
(7.32)
O i = C 0 Ii + C 1 Ii − 1 + C 2 O i − 2
(7.33)
donde:
− KX + 0.5 Δt
C0 =
K − KX + 0.5 Δt
(7.34)
KX + 0.5 Δt
C1 =
K − KX + 0.5 Δt
(7.35)
K − KX − 0.5 Δt
C2 =
K − KX + 0.5 Δt
(7.36)
C 0 + C 1 + C 2 = 1 .0
(7.37)
270
http://estudiantesingcivil.blogspot.mx/
1/6
L
K ≅ 0.014
1/6
H
(7.38)
donde H es el desnivel entre las dos secciones de los extremos del tramo del
río, en km; L es la longitud del tramo del río, en km.
Este criterio se aplica cuando se dispone de una avenida que haya sido
aforada en los extremos del tramo de interés del río (Fuentes y Franco,
1999).
⎡I + I Oi + Oi +1 ⎤
i i +1
ΔV = Δt ⎢ − ⎥
⎢ 2 2 ⎥
⎣ ⎦
(7.39)
4) Se supone un valor de X
Y = X I +(1 − X) O
(7.40)
271
8) Una vez seleccionado el valor adecuado de X, el parámetro K resulta ser
igual a la pendiente de la recta de ajuste (ver figura 7.8).
VT VT
272
Las expresiones que permiten estimar el valor de los parámetros K y X son:
⎛ ⎞
K = 1.41 ⎜ Tp − t p ⎟
⎝ ⎠
(7.41)
t p ⎡⎢ Ip − Op ⎤⎥
X = 0.71
K ⎢ Ip ⎥
⎣⎢ ⎦⎥
(7.42)
2) Se conoce I1 y O1
Oi = C 0 Ii + C1 Ii −1 + C 2 Oi − 2
(7.43)
273
Ejemplo 7.1. Tránsito de una avenida en un vaso de almacenamiento, a
partir del método del algoritmo de incrementos finitos.
3 /2 3/2
Q = CL (hi − H) = 60 (hi − 30)
(7.44)
274
Figura 7.9. Hidrograma de entrada o avenida de diseño
4.5
V = 1.4 h
(7.45)
NK 4.5 (1.4)
= = 0.0035
Δt 1800
275
Solución:
276
Tabla 7.2. Evolución de los gastos de descarga y elevaciones de los niveles del agua
en el vaso de almacenamiento
3 3 3
i I i , m /s I i+ 1 , m /s h i , me t ro s Q i , m /s a h i+ 1 , me t ro s
0 0 100 3 0 .0 0 0 .0 0 5 1 7 .6 0 3 0 .1 0
1 100 400 3 0 .1 0 1 .9 0 5 2 4 .1 4 3 0 .5 7
2 400 500 3 0 .5 7 2 5 .8 2 5 5 3 .9 8 3 1 .3 4
3 500 460 3 1 .3 4 9 3 .0 7 6 0 4 .8 7 3 1 .9 8
4 460 350 3 1 .9 8 1 6 7 .1 7 6 4 9 .4 7 3 2 .3 5
5 350 230 3 2 .3 5 2 1 6 .1 5 6 7 6 .2 5 3 2 .4 6
6 230 140 3 2 .4 6 2 3 1 .5 0 6 8 4 .3 6 3 2 .3 9
7 140 80 3 2 .3 9 2 2 1 .6 9 6 7 9 .1 9 3 2 .2 3
8 80 30 3 2 .2 3 1 9 9 .8 1 6 6 7 .4 8 3 2 .0 1
9 30 0 3 2 .0 1 1 7 0 .9 8 6 5 1 .6 1 3 1 .7 7
10 0 0 3 1 .7 7 1 4 1 .2 9 6 3 4 .5 9 3 1 .5 5
11 0 0 3 1 .5 5 1 1 5 .7 8 6 1 9 .2 7 3 1 .3 6
12 0 0 3 1 .3 6 9 5 .1 6 6 0 6 .2 3 3 1 .2 0
13 0 0 3 1 .2 0 7 8 .8 7 5 9 5 .4 0 3 1 .0 7
14 0 0 3 1 .0 7 6 6 .4 1 5 8 6 .7 0 3 0 .9 6
15 0 0 3 0 .9 6 5 6 .4 4 5 7 9 .3 9 3 0 .8 6
16 0 0 3 0 .8 6 4 7 .8 5 5 7 2 .8 1 3 0 .7 8
17 0 0 3 0 .7 8 4 1 .3 3 5 6 7 .5 7 3 0 .7 1
18 0 0 3 0 .7 1 3 5 .9 0 5 6 3 .0 2 3 0 .6 5
19 0 0 3 0 .6 5 3 1 .4 4 5 5 9 .1 3 3 0 .5 9
20 0 0 3 0 .5 9 2 7 .1 9 5 5 5 .2 6 3 0 .5 4
21 0 0 3 0 .5 4 2 3 .8 1 5 5 2 .0 5 3 0 .5 0
22 0 0 3 0 .5 0 2 1 .2 1 5 4 9 .4 9 3 0 .4 6
23 0 0 3 0 .4 6 1 8 .7 2 5 4 6 .9 3 3 0 .4 3
24 0 0 3 0 .4 3 1 6 .9 2 5 4 5 .0 2 3 0 .4 0
25 0 0 3 0 .4 0 1 5 .1 8 5 4 3 .1 1 3 0 .3 7
26 0 0 3 0 .3 7 1 3 .5 0 5 4 1 .2 0 3 0 .3 5
27 0 0 3 0 .3 5 1 2 .4 2 5 3 9 .9 3 3 0 .3 3
28 0 0 3 0 .3 3 1 1 .3 7 5 3 8 .6 6 3 0 .3 1
29 0 0 3 0 .3 1 1 0 .3 6 5 3 7 .4 0 3 0 .2 9
30 0 0 3 0 .2 9 9 .3 7 5 3 6 .1 3 3 0 .2 7
277
Figura 7.10. Hidrogramas de entrada y de salida obtenidos al realizar el tránsito
de la avenida con el método del algoritmo de incrementos finitos
278
Con el método de Muskingum transitar la avenida entre el tramo A-B y
determinar si el poblado localizado en B se inunda, si en sitio A se presenta
el hidrograma mostrado en la figura 7.12.
279
Solución:
280
http://estudiantesingcivil.blogspot.mx/
281
Figura 7.14 b). Relaciones entre el almacenamiento y los gastos
282
Figura 7.15. Estimación del parámetro K del método de MUskingum
(1 ) (2 ) (3 ) (4 ) (5 ) (6 ) (7 ) (8 ) (9 ) (1 0 ) (1 1 ) (1 2 ) (1 3 ) (1 4 ) (1 5 )
t, I, O, I-O , V a lm a ce n a d o VT X =0 .1 X =0 .3 X =0 .1 3
3 3 3 3 3
h m /s m /s m /s (4 h )m /s (4 h )m /s Y = X I + (1 -X ) O Y = X I + (1 -X ) O Y = X I + (1 -X ) O
0 .1 I 0 .9 O To ta l 0 .3 I 0 .7 O To ta l 0 .1 3 I 0 .8 7 O To ta l
0 20 20 0 0 .0 0 2 .0 1 8 .0 2 0 .0 6 .0 1 4 .0 2 0 .0 2 .6 1 7 .4 2 0 .0
4 106 30 76 3 8 .0 3 8 .0 1 0 .6 2 7 .0 3 7 .6 3 1 .8 2 1 .0 5 2 .8 1 3 .8 2 6 .1 3 9 .9
8 75 43 32 5 4 .0 9 2 .0 7 .5 3 8 .7 4 6 .2 2 2 .5 3 0 .1 5 2 .6 9 .8 3 7 .4 4 7 .2
12 60 56 4 1 8 .0 1 1 0 .0 6 .0 5 0 .4 5 6 .4 1 8 .0 3 9 .2 5 7 .2 7 .8 4 8 .7 5 6 .5
16 55 58 -3 0 .5 1 1 0 .5 5 .5 5 2 .2 5 7 .7 1 6 .5 4 0 .6 5 7 .1 7 .2 5 0 .5 5 7 .6
20 45 58 -1 3 -8 .0 1 0 2 .5 4 .5 5 2 .2 5 6 .7 1 3 .5 4 0 .6 5 4 .1 5 .9 5 0 .5 5 6 .3
24 30 54 -2 4 -1 8 .5 8 4 .0 3 .0 4 8 .6 5 1 .6 9 .0 3 7 .8 4 6 .8 3 .9 4 7 .0 5 0 .9
28 20 47 -2 7 -2 5 .5 5 8 .5 2 .0 4 2 .3 4 4 .3 6 .0 3 2 .9 3 8 .9 2 .6 4 0 .9 4 3 .5
32 20 37 -1 7 -2 2 .0 3 6 .5 2 .0 3 3 .3 3 5 .3 6 .0 2 5 .9 3 1 .9 2 .6 3 2 .2 3 4 .8
36 20 31 -1 1 -1 4 .0 2 2 .5 2 .0 2 7 .9 2 9 .9 6 .0 2 1 .7 2 7 .7 2 .6 2 7 .0 2 9 .6
40 20 27 -7 -9 .0 1 3 .5 2 .0 2 4 .3 2 6 .3 6 .0 1 8 .9 2 4 .9 2 .6 2 3 .5 2 6 .1
44 20 24 -4 -5 .5 8 .0 2 .0 2 1 .6 2 3 .6 6 .0 1 6 .8 2 2 .8 2 .6 2 0 .9 2 3 .5
283
b) El proceso para estimar los volúmenes almacenados que aparecen en la
columna (5) de la tabla 7.3, para los intervalos de tiempo, se indica a
continuación:
0 + 76 76
V0−4 = = = 38
2 2
76 + 32 108
V4−8 = = = 54
2 2
32 + 4 36
V8−12 = = = 18
2 2
4 + (−3) 1
V12−16 = = = 0.5
2 2
− 3 − 13 − 16
V16 − 20 = = = −8.0
2 2
− 13 − 24 − 37
V20 − 24 = = = −18.5
2 2
− 24 − 27 − 51
V24 − 28 = = = −25.5
2 2
− 27 − 17 − 44
V28 − 32 = = = −22.0
2 2
− 17 − 11 − 28
V32 − 36 = = = −14.0
2 2
− 11 − 7 − 18
V36 − 40 = = = −9.0
2 2
− 7 − 4 − 11
V40 − 44 = = = −5.5
2 2
284
− KX + 0.5 Δt − 11.43 (0.13) + 0.5 (4)
C0 = = = 0.04
K − KX + 0.5 Δt 11.43 − 11.43 (0.13) + 0.5 (4)
C 0 = 0.04
Oi = C 0 Ii + C1 Ii−1 + C 2 Oi−2
285
Tabla 7.4. Aplicación del método de Muskingum
In g re s o E g re s o
i Ho ra Ii, 0 .0 4 I i , 0 .2 9 I i-1 , 0 .6 7 O i-2 , O i,
3 3 3 3 3
e n m /s e n m /s e n m /s e n m /s e n m /s
1 0 20 .... .... .... 2 0 .0 0
2 4 120 4 .8 0 5 .8 0 1 3 .4 0 2 4 .0 0
3 8 100 4 .0 0 3 4 .8 0 1 6 .0 8 5 4 .8 8
4 12 80 3 .2 0 2 9 .0 0 3 6 .7 7 6 8 .9 7
5 16 60 2 .4 0 2 3 .2 0 4 6 .2 1 7 1 .8 1
6 20 40 1 .6 0 1 7 .4 0 4 8 .1 1 6 7 .1 1
7 24 30 1 .2 0 1 1 .6 0 4 4 .9 7 5 7 .7 7
8 28 25 1 .0 0 8 .7 0 3 8 .7 0 4 8 .4 0
9 32 20 0 .8 0 7 .2 5 3 2 .4 3 4 0 .4 8
10 36 20 0 .8 0 5 .8 0 2 7 .1 2 3 3 .7 2
11 40 20 0 .8 0 5 .8 0 2 2 .5 9 2 9 .1 9
12 44 20 0 .8 0 5 .8 0 1 9 .5 6 2 6 .1 6
286
Referencias
287
Principios y Fundamentos de la Hidrología Superficial, se
editó en formato electrónico el 10 de enero de 2006,
en la Coordinación General de Vinculación y Desarrollo
Institucional de la Rectoría General de la Universidad
Autónoma Metropolitana, Tlalpan, D.F. México.