Hid Ro Log A para Ingeniero S
Hid Ro Log A para Ingeniero S
Hid Ro Log A para Ingeniero S
INGENIEROS
Facultad de Minas
Sede Medellín
2000
Capítulo 1
GENERALIDADES
1
La palabra hidrología proviene de las raíces griegas “hidros” agua y
“logos” tratado o sea que la hidrología es la ciencia del agua. Según Ven
Te Chow (1994), la hidrología estudia el agua existente en la tierra, su
distribución, sus propiedades físicas y químicas y su influencia sobre el
medio ambiente.
2
1.1 BREVE HISTORIA DE LA HIDROLOGÍA
En el siglo I DC, Marcos Vitruvio sentó las bases para la teoría del ciclo
hidrológico. El lanzó la hipótesis de que la lluvia y la nieve caían en las
áreas montañosas, se infiltraban, y más tarde esa agua aparecía en las
tierras bajas y en las corrientes.
3
La llegada de lo que puede llamarse hidrología moderna empieza con los
estudios de los pioneros como Perrault, Mariotte y Halley en el siglo
XVII. Perrault obtiene medidas de la precipitación en la cuenca del Sena
para un período de 3 años. Midió también evaporación y capilaridad.
Mariotte midió la velocidad en el Sena. Estas fueron luego empleadas
para calcular el caudal.
4
Capítulo 2
EL CICLO HIDROLÓGICO
5
El ciclo hidrológico es un proceso continuo en el que el agua de los océanos
es transportada a la atmósfera, después a la tierra, para regresar luego al
océano. Hay muchos subciclos. La evaporación del agua tierra adentro, su
subsecuente precipitación sobre la tierra y el retorno de esta agua al océano
son un ejemplo. La energía necesaria para este proceso es suministrada por el
sol. El ciclo hidrológico puede representarse esquemáticamente según la figura
2.1.
6
FIGURA 2.1 El ciclo hidrológico
7
superficie y la otra cae al suelo por la acción de la gravedad. De esta última,
una parte es interceptada por la vegetación o por otras superficies y luego se
evapora igualmente, la otra entra en contacto directo con el suelo o con las
extensiones de agua libre (lagos y mares principalmente). Aquella que cae al
suelo puede alimentar la red hidrográfica directamente (escorrentía
superficial), almacenarse en algunas partes de la superficie (almacenamiento
superficial) o penetrar en el suelo y subsuelo (infiltración) para de esta
manera alimentar las reservas de las capas acuíferas.
La energía calorífica proveniente del sol evapora parte del agua que está en
las superficies libres (lagos, mares, etc.), en la red hidrográfica y en la
superficie del suelo. De igual manera una parte del agua que está en la zona
no-saturada del suelo se evapora directamente o se absorbe por las plantas, las
cuales a su vez la devuelven a la atmósfera en forma de vapor, mediante el
mecanismo de transpiración. Estos fenómenos son generalmente designados
con el término de evapotranspiración. En cuanto al agua que se ha infiltrado
en las zonas más profundas de la corteza terrestre y que alimenta las capas
acuíferas, puede pasar a alimentar el flujo de los ríos, o retornar directamente
a los mares a causa de la percolación, para reiniciar de nuevo el ciclo.
• Las fuentes principales de energía que mueven el ciclo del agua a través de
sus diferentes fases son la energía calorífica proveniente del sol y la
energía gravitacional.
Precipitación (P)
Evapotranspiración (ET)
Escorrentía superficial (ESD)
Infiltración (I)
En el desarrollo de estos procesos el agua se encuentra, como ya se anotó, en
diferentes estados y en distintos almacenamientos. Así una partícula de agua
que dé el ciclo completo y que se encuentra en el océano pasa, en razón de la
8
evaporación, a la atmósfera, después de un cierto tiempo se condensa y cae a
la superficie bajo la forma de lluvia y luego, como escorrentía superficial o
subterránea, puede volver al mar. Sin embargo no necesariamente en una
cuenca, ante una lluvia, se presentan simultáneamente todas las fases del
ciclo: por ejemplo, en un terreno muy pendiente e impermeable , la
escorrentía superficial puede ser la única fase que se presente durante el
evento de lluvia, y la evapotranspiración e infiltración pueden ser
insignificantes frente a ésta.
9
SISTEMA
ENTRADAS
SALIDAS
En hidrología, como en muchas otras áreas, las leyes físicas que rigen un
determinado proceso son difíciles de determinar o muy complejas de aplicar, o
la geometría del sistema es muy complicada o la falta de homogeneidad es
muy grande. Los hidrólogos generalmente tratan con un sistema abierto, tal
como el mostrado por la figura 2.1. La cuantificación del agua presente en este
sistema puede hacerse con una simple ecuación de balance de masas así:
ds
= I−Q
dt
10
prácticamente nula y que por consiguiente las entradas son iguales a las
salidas.
Ahora bien, si el objeto de estudio no es ya el ciclo global del agua sino una
zona más reducida como, por ejemplo, una cuenca hidrográfica, cabe (aquí)
también la aplicación del concepto de sistema y se puede considerar dicha
cuenca como un subsistema que hace parte del sistema más complejo
constituido por el ciclo global del agua. Es decir, que a nivel de la cuenca,
dada una entrada (precipitaciones por ejemplo), el sistema produce una salida
(caudal en un punto dado). Esquemáticamente:
Entradas Salidas
Precipitación SISTEMA Caudal en
la estación
Así, las precipitaciones (P) que caen sobre una cuenca durante un tiempo
dado se reparten de la manera siguiente: la infiltración (I), la
evapotranspiración real (ET) y la escorrentía superficial (Qs). Considerando
un período suficiente largo, puede afirmarse que el balance global promedio
es:
P = I + ET + Qs (2.1)
La infiltración total (I) puede tenerse como formada por dos fracciones: una
que sirve para renovar el almacenamiento de agua en la zona de aireación (Ia)
y otra denominada infiltración eficaz (Ie), que alimenta las aguas
subterráneas; de esta última, una parte alimenta el caudal de los cursos de
agua como caudal de base.
11
geográficos de la cuenca hidrográfica son los límites de los acuíferos, ver
figura 2.2.
P = I + ET + Qs + Q w + Qex + q (2.2)
en donde:
Qw : Escorrentía subterránea
Qex : Aguas extraídas o explotadas (bien sean aguas subterráneas o
superficiales)
q : Pérdidas y entradas hacia o desde otras cuencas vecinas.
12
Físicos:
13
solo las colonias de peces sino también los arrecifes coralinos de zonas como
Isla Fuerte e Islas del Rosario. El problema de sedimentos afecta la navegación
en los ríos. En Colombia, el río Magdalena que era navegable por buques de
gran calado hace 50 años desde Puerto Berrío, es hoy prácticamente
descartado para la navegación en casi la totalidad de su recorrido.
Químicos:
14
recorrido urbano han perdido toda señal de flora y fauna naturales. La
contaminacción está ligada íntimamente a las condiciones biológicas de
corrientes y cuerpos de agua, tal como se ha anotado anteriormente.
15
Tabla 2.1 Distribución del agua en la tierra.
16
Capítulo 3
ELEMENTOS DE METEOROLOGÍA
17
Los principales factores que afectan el almacenamiento (y movimiento)del
agua en la atmósfera son la temperatura, presión y radiación solar. Las
diferentes condiciones de la atmósfera (humedad, temperatura, presión y
movimiento del viento), en cualquier lugar y durante períodos
determinados, es lo que se conoce como clima. Se estudiarán en este
capítulo los elementos que influyen en el clima del planeta y se precisaran
algunos aspectos del clima colombiano.
3.1 LA ATMOSFERA
18
Desde el punto de vista hidrológico, los gases más importantes de la
atmósfera son el vapor de agua, el anhídrido carbónico y el ozono.
ionosfera
80 km
estratosfera
18 km
8 km
troposfera
19
TABLA 3.1 Gases presentes en la atmósfera.
Gases % en volumen
Nitrógeno 78.
Oxígeno 20.95
Argón 0.93
Anhídrido carbónico* 0.03
Neón 18.18*10-4
Helio 5.24*10-4
Criptón 1.14*10-4
Hidrógeno 0.5*10-4
Oxido Nitroso* 0.5*10-4
Metano 1.5*10-4
Xenón 0.07*10-4
Gas sulfuroso* ----
Gas sulfídrico1*
1
* valores variables
20
intensidad de los rayos ultravioleta que alcanzan la tierra. Esto se ha
traducido en un aumento de los casos de cáncer de piel en algunas regiones
del mundo afectadas por el fenómeno.
El sol tiene una temperatura promedia de 6000 grados centígrados; los gases
calientes, incandescentes, que comprimen su superficie, emiten energía en
forma de radiaciones electromagnéticas constituidas por un espectro de
ondas de un amplio rango de longitudes, que viajan a una velocidad de
300000 km/s. La energía que es emitida radialmente desde el sol requiere de
aproximadamente 8 1/3 minutos para viajar los 150 millones de kilómetros
que lo separan de la tierra.
Aunque la energía solar viaja a través del espacio sin pérdidas significativas
de energía, los rayos divergen y se mueven en todas direcciones, lejos del
sol. La intensidad de la radiación decrece inversamente con el cuadrado de
la distancia. La tierra solo recibe de una a dos billonésimas partes de la
energía que el sol emite. Esta cantidad es aproximadamente1.94 cal/cm2*
min. La energía se expresa también en Langleys/minuto. Un langley (Ly) es
igual a una caloría por cm2 y por minuto.
21
En el universo todos los cuerpos irradian energía y además pueden
absorberla. La relación entre la energía emitida y la absorbida se puede
calcular por la ley de Kirehhof así:
Eλ
= f (λ , T) (3.1)
aλ
Donde:
E λ = f ( λ , T) (3.2)
4
F = σT (3.3)
Donde:
22
Albedo= Rr/RI
Donde:
Rr : radiación reflejada.
R I : radiación incidente
Nubes 5-84
Hormigón 17-27
Carreteras pavimentadas 5-10
Bosque cerrado 5-10
Arena blanca 35
Vegetales y pastos 20-25
Suelo sin vegetación 10-25
Nieve limpia 70-90
Agua 5-35
23
3.2.1 El efecto invernadero
3.2.2 La temperatura
24
temperaturas mayores que las inferiores. Este fenómeno se llama inversión
de temperatura. Se favorece cuando el cielo está nublado, con humedades
bajas y vientos en calma, que evitan que se produzca la turbulencia,
necesaria, para mezclar las masa de aire de diferente temperatura. Un
período de tan sólo tres días de poca mezcla atmosférica puede llevar a
concentraciones elevadas de productos peligrosos, en áreas de alta
contaminación y, en casos extremos, producir enfermedades e incluso la
muerte. En 1948 una inversión térmica sobre Donora, Pennsylvania, produjo
enfermedades respiratorias en más de 6.000 personas ocasionando la muerte
de veinte de ellas. En Londres, la contaminación segó entre 3.500 y 4.000
vidas en 1952, y otras 700 en 1962. La liberación de isocianato de metilo a
la atmósfera durante una inversión térmica fue la causa del desastre de
Bhopal, India, en diciembre de 1984, que produjo al menos 3.300 muertes y
más de 20.000 afectados (Encarta, 1998).
3.2.3 La humedad
25
completamente saturada. La presión en este momento es llamada presión de
vapor, para esa temperatura. Aplicando la ley de los gases perfectos, la
presión de vapor, e, se puede expresar como:
e
ρ v = 0.622 (3.4)
RT
Donde:
ρv: densidad de la masa de vapor de agua, en g/m3.
e: presión de vapor, en milibares.
T: temperatura absoluta, en grados Kelvin
R: constante de los gases ideales.
La densidad del vapor de agua, ρv, es llamada también humedad absoluta de
la atmósfera
ρv e
HR = 100 = 100 (3.5)
ρs es
Donde:
HR: humedad relativa.
es: presión de saturación
ρv: densidad de la masa de aire o humedad absoluta.
ρs: densidad de la masa de aire saturado.
e: presión de vapor.
es: presión de vapor saturada
26
La presión de saturación es función directa de la temperatura y es bastante
sensible a ésta; se presentan grandes variaciones a nivel diario, debidas a
esta dependencia, ver figura 3.2.
e
q = 0.622 (3.6)
pa
27
Relación de mezcla de masa: se define como:
e
r = 0.622 (3.7)
p ae
Si la tierra fuera una esfera que no rotara sobre su propio eje, la circulación
atmosférica sería tal como la mostrada por la figura 3.3. El aire del Ecuador,
más caliente y menos denso se elevaría y viajaría hacia los polos, donde una
vez enfriado descendería hacia la baja atmósfera y retornaría al Ecuador.
Este modelo de circulación se conoce con el nombre de circulación
Hadley.
28
FIGURA 3.3 Patrón ideal de circulación atmosférica
En las latitudes 308 N y 308S, parte del aire elevado en el Ecuador desciende
de nuevo hacia la superficie, ya que es mas frío que el que está en las capas
inferiores. Parte del aire se devuelve para el Ecuador y parte continúa hacia
el norte
29
FIGURA 3.4 Esquema de la circulación atmosférica
La latitud de los 308 es llamada “de los caballos”(El origen de este nombre
se debe a que los españoles llevaron caballos al Nuevo Mundo, y en esta
zona encontraron vientos más suaves, lo que los obligaba a viajar más
despacio. Por falta de alimento los animales tenían que ser arrojados al
mar). En la latitud de los caballos el aire generalmente es descendente, hay
poca nubosidad, la atmósfera es seca y la precipitación y los vientos son
escasos. En esta franja se encuentran los más grandes desiertos del mundo.
Los movimientos descendentes de masas de aire caliente alrededor de la
Tierra han producido dos cinturones de desiertos, uno a lo largo del trópico
de Cáncer, en el hemisferio norte, y otro a lo largo del trópico de
Capricornio, en el hemisferio sur. Entre los desiertos del hemisferio norte
cabe destacar el de Gobi, entre China y Mongolia, los desiertos del suroeste
de Norteamérica, el Sahara en el norte de África y los desiertos Arábigo y
de Irán en el Oriente Próximo. A lo largo del cinturón sur se encuentra la
Patagonia en Argentina, el desierto de Kalahari en el sur de África, y el
Gran DesiertoVictoria y el Gran Desierto de Arena en Australia, ver figura
3.5 (Enciclopedia Encarta, 1999).
30
FIGURA 3.5 Grandes desiertos del mundo.( Enciclopedia Encarta
1999)
La parte de aire que se dirige de los caballos hacia el polo, entre latitudes
408 y 608 , se encuentra con las corrientes frías del polo en dirección
contraria y se forma un frente polar; el aire caliente más liviano que viene
desde la latitud de los caballos es obligado a levantarse sobre la masa fría.
En las latitudes 60 norte y sur, se desarrolla un franja de baja presión con
climas variables, de precipitación moderada y en algunos casos alta; allí se
presentan grandes tormentas.
En las zonas polares se presentan focos de alta presión, con aire seco y baja
precipitación y evaporación
31
ascenso y posterior condensación de estas masas de aire húmedo. Los
patrones de circulación de los alisios sugieren que ellos convergen en algún
lugar cerca del Ecuador. El encuentro de los alisios del norte y los alisios
del sur se da en una franja estrecha, denominada frente intertropical de
convergencia, FIC. Esta zona es uno de los más importantes sistemas
climáticos que actúan en los trópicos y es decisiva en la caracterización del
clima en diversas áreas de la región tropical. El FIC se encuentra en una
zona donde hay interacción de fenómenos atmosféricos y climáticos, así: Es
la zona de confluencia de los alisios, la zona donde se presentan las
máximas temperaturas de la superficie del mar, la zona de máxima
convergencia de masas de aire y la zona donde se presenta la mayor
cantidad de nubes convectivas. Sin embargo estos fenómenos no
necesariamente ocurren simultáneamente en la misma latitud. Las
características climáticas asociadas al FIC, viajan en dirección Norte sur a lo
largo del año. El FIC alcanza su posición más al norte (88N) durante el
verano del hemisferio norte y su posición más al sur (18 N) durante el mes
de Abril. Además de esta oscilación anual el FIC, presenta oscilaciones de
alta frecuencia con períodos que varían de días a semanas.
32
más. Los científicos reservan el término del El Niño, para significar eventos
mensuales asociados con las anomalías de la temperatura de la superficie del
mar, que pueden durar desde meses hasta más de un año.
Durante los últimos cuarenta años, nueve "Niños" han afectado la costa de
América del Sur. En la mayoría de ellos la temperatura del agua no sólo se
elevó en la costa sino también en las islas Galápagos y a lo largo de una
franja de 5000 millas sobre el Pacífico Ecuatorial. En los eventos más
débiles, las temperaturas se elevaron sólo 1 a 2 grados centígrados Sin
embargo, los eventos fuertes como " El Niño" de 1982 - 83 dejaron una
profunda huella en las condiciones climáticas a lo largo del mundo entero.
Según Kiladis y Díaz(1989), los años Niño y Niña de este siglo son los
siguientes:
AÑOS NIÑO:
1902,1904,1911,1913,1918,1923,1925,1930,1932,1939,1951,1953,1957-
1958,1963,1963,1965,1969,1972,1976-1977,1982-1983,1986-1987,1991-
1992,1997.
33
AÑOS NIÑA:
1903,1906,1908,1916,1920,1924,1928,1931,1938,1942,1949,1954,1964,19
70-1971,1973,1975,1988
Para entender como "El Niño" afecta los océanos, es necesario comprender
primero cómo los vientos superficiales movilizan el agua durante los años
normales y cómo dichos movimientos afectan la temperatura de las aguas y
las cantidades de nutrientes químicos disponibles en la cadena alimenticia.
Se considera, para mejor comprensión, dos regiones separadas: el Pacífico
Ecuatorial, que se extiende en dirección oeste desde las Islas Galápagos
hasta pasada la línea de cambio de fecha y las aguas costeras del Perú y sur
del Ecuador. Los vientos del este que soplan a lo largo del Ecuador y los
vientos surestes que soplan a lo largo de las costas del Perú y Ecuador
tienden a arrastrar con ellos el agua superficial. La rotación de la tierra a su
vez desvía las corrientes superficiales a la derecha en el Hemisferio Norte y
hacia la izquierda en el Hemisferio Sur. Las aguas superficiales se desvían
del Ecuador en ambas direcciones y se alejan de la costa. Conforme se
movilizan y se aleja el agua superficial, agua más fría y rica en nutrientes de
mayores profundidades asciende a la superficie para reemplazarla,
generando el fenómeno conocido como afloramiento. Ambos
afloramientos, el ecuatorial y el costero, se encuentran concentrados en fajas
relativamente estrechas de menos de 100 millas de ancho
Durante los años de "El Niño", cuando los vientos del este se repliegan al
Pacífico Oriental, el océano responde de la siguiente manera:
• La termoclina (línea que separa las capas frías y calientes del
océano) a lo largo del Ecuador se aplana, levantándose en el oeste y
hundiéndose en el este más de 100 metros bajo la superficie a
suficiente profundidad, evitando así que el afloramiento costero sea
34
capaz de traer a la superficie aguas frías, más ricas en nutrientes, de
las capas inferiores;
• El afloramiento ecuatorial se reduce, disminuyendo aún más el
abastecimiento de nutrientes a la cadena alimenticia.
• La lengua fría de la temperatura superficial del mar se debilita o
desaparece
• El nivel del mar se aplana, bajando en el oeste y subiendo en el este,
produciéndose ondas superficiales hacia el este a lo largo del
Ecuador.
Cuando los vientos del este soplan con toda intensidad, el afloramiento de
agua fría a lo largo del pacífico ecuatorial enfría a su vez el aire que se
encuentra por encima, haciéndolo demasiado denso para elevarse lo
suficiente y producir vapor de agua que se condense para formar nubes y
gotas de lluvia. Como resultado, durante los años normales, esta franja del
océano se queda visiblemente sin nubes y la lluvia, en el cinturón ecuatorial,
se limita mayormente al extremo oeste del Pacífico, cerca de Indonesia.
Pero cuando los vientos del este se debilitan en las etapas iniciales de un
evento de "El Niño", el afloramiento disminuye y el océano se calienta, al
igual que el aire húmedo que se encuentra por encima (de éste).
35
FIGURA 3.6 Condiciones normales y condiciones Niño (NOA 1999)
36
El aire se vuelve suficientemente boyante para formar nubes espesas que
producen fuertes lluvias a lo largo del Ecuador. El cambio en las
temperaturas del océano hace que la zona de mayores lluvias, característica
del pacífico oeste, se traslade en dirección este. Los giros y cambios en el
diálogo océano-atmósfera sostenido en el Pacífico pueden traer
repercusiones en las condiciones climáticas en regiones muy distantes
alrededor del mundo (este fenómeno se conoce como teleconexiones). Este
mensaje de cobertura mundial es transmitido a través de cambios en la
precipitación pluvial en los trópicos, afectando los patrones del viento sobre
gran parte del mundo.
37
ESTACIÓN CARTAGO (RÍO LA VIEJA)
600
500
Caudal máximo instantáneo (m³/s)
400
300
200
100
0
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Mes
Niño Niña Total
38
3.5 EL CLIMA EN COLOMBIA
Colombia está situada en los trópicos, al norte del Ecuador y al sur de los
128N. Además de los factores generales de circulación global y del cambio
en la posición aparente del sol durante el año, la topografía, la convección
profunda, la cercanía a las costas y la vegetación son los factores principales
que gobiernan el clima local en los trópicos, dado que la circulación general
es débil en gradientes de presión, temperatura, humedad y velocidad del
viento.
En Colombia se dan todos los extremos climáticos: existe la zona más seca
en el extremo norte de la Guajira, y las zonas más lluviosas del planeta en la
costa Pacífica. Se presentan desde los lugares más cálidos hasta los glaciares
más fríos.
39
FIGURA 3.7 Isoyetas de precipitaciön promedia anual en Colombia
(IDEAM 1999)
40
Capítulo 4
GEOMORFOLOGIA DE CUENCAS
41
En las ciencias de la tierra ha sido reconocida la influencia de la geomorfología
en la interacción de la geología, el clima y el movimiento del agua sobre la
tierra. Esta interacción es de gran complejidad y prácticamente imposible de
ser concretada en modelos determinísticos, y se debe tomar como un proceso
de comportamiento mixto con una fuerte componente estocástica.
42
preferiblemente hayan sido elaborados bajo los mismos criterios cartográficos.
De esta forma se contaría con resultados homogéneos, que podrían ser
comparados en estudios posteriores al estudio mismo de las cuencas.
Toda cuenca en estudio debe estar delimitada en cuanto a su río principal, tanto
aguas abajo como aguas arriba. Aguas abajo, idealmente por la estación de
aforo más cercana a los límites de la cuenca en que se está interesado siendo el
punto de la estación el (punto) más bajo en el perfil del río. Aguas arriba, por
otra estación que sea el punto más alto en el perfil del río en que se incluya el
área en estudio, o por las cabeceras del río si es el caso del estudio de la cuenca
desde el nacimiento.
43
contorno o de divorcio de la cuenca. Realmente la definición de dicha línea no
es clara ni única, pues pueden existir dos líneas de divorcio: una para las aguas
superficiales, que sería la topográfica, y otra para las aguas subsuperficiales,
que sería determinada en función de los perfiles de la estructura geológica,
fundamentalmente por los pisos impermeables (Fig 4.1).
44
producen hidrógrafas con variaciones en el tiempo más suaves y más llanas.
Sin embargo, en cuencas grandes, se pueden dar hidrógrafas picudas cuando la
precipitación fué intensa, en las cercanías, aguas arriba, de la estación de aforo.
Q = C An (4.1)
45
Antioquia
Donde :
A: área de la cuenca en km2
Q: media de los caudales máximos instantáneos en m3/s.
46
3
Q1 ⎛ A1 ⎞ 4
= ⎜⎜ ⎟⎟ (4.3)
Q2 ⎝ A 2 ⎠
47
4.1.2 Longitud, perímetro y ancho.
A
W= (4.4)
L
48
4.1.3 Parámetros de forma de la cuenca
-0.136
A
2
=A (4.5)
L 2
49
cual se han hecho numerosos esfuerzos para tratar de cuantificar este efecto
por medio de un valor numérico.
A
Rf = (4.7)
L2b
50
círculo que contenga la misma área A de la cuenca hidrográfica:
P
K = 0.282 (4.8)
A
La razón para usar la relación del área equivalente a la ocupada por un círculo
es porque una cuenca circular tiene mayores posibilidades de producir
avenidas superiores dada su simetría. Sin embargo, este índice de forma ha
sido criticado, pues las cuencas en general tienden a tener forma de pera.
Son muy importantes, ya que el relieve de una cuenca puede tener más
influencia sobre la respuesta hidrológica que la forma misma de la cuenca.
Los parámetros relativos al relieve son:
51
bastante utilizado es el siguiente: se monta sobre la cuenca una cuadrícula de
tamaño conveniente, se cuentan los cortes de las curvas de nivel con los ejes
horizontal y vertical de la cuadrícula, respectivamente, y se tiene:
nhh
Sh =
Lh
nvh
Sv =
Lv
donde:
h es la diferencia de cotas entre curvas de nivel.
nh es el número de cruces de las curvas de nivel con líneas de igual
coordenada este.
nv es el número de cruces de las curvas de nivel con líneas de igual
coordenada norte.
Sh y Sv son la pendiente horizontal y vertical de la cuenca, respectivamente.
+
S = Se Sn x100 % (4.9)
2
52
FIGURA 4.4. Método para hallar la pendiente S, en una cuenca
53
FIGURA 4.5 Curvas hipsométricas características
54
está referida al nivel
0.30
0.15
0.10
0.05
0.00
1000- 1500- 1600- 1700- 1800- 1900-
1500 1600 1700 1800 1900 2000
Intervalo de Alturas
del mar. Este valor puede ser encontrado usando la curva hipsométrica o el
histograma de frecuencias altimétricas. Esta estimación se realiza por una
media aritmética ponderada en el caso del histograma, o de la curva
hipsométrica calculando el área bajo la curva y dividiéndola por el área total.
55
Figura 4.7 Hidrógrafas según el perfil altimétrico del cauce principal.
56
d) Usando cualquiera de los métodos anteriores, pero sin tener en cuenta
toda la trayectoria del cauce principal, ignorando por lo tanto de un 10%
a un 15% de los tramos extremos (nacimiento y desembocadura).
El valor del orden del canal principal,ς, en la boca de la cuenca da una idea de
la magnitud del drenaje de la cuenca.
57
superior:
Nn
Rb = (4.10)
N n +1
Donde:
Nn es el número de canales de orden n y Nn+1 es el número de canales de
orden n+1.
El valor "medio" de bifurcación, Rb, de una cuenca se determina mediante la
pendiente de la recta que resulta de graficar el logaritmo decimal del número
de corrientes de cada orden en el eje de las ordenadas y el orden de las
corrientes en el eje de las abscisas, por medio de un ajuste de mínimos
cuadrados. El valor "medio" se toma como el antilogaritmo de la pendiente
de la recta ajustada a las parejas de valores.
58
Con base en estudios estadísticos se le ajustó la siguiente relación:
k -n
Nn = R b _ log N n = (k - n) log R b (4.11)
donde:
K: orden mayor de los canales de la cuenca en estudio
n: orden del canal en el cual estamos interesados
Nn: número de canales para el orden n
L i = l 1rli −1 (4.12)
A i = A 1rai −1 (4.13)
Donde:
59
Ai: área de las cuencas de orden i
A1: área de las cuencas de orden 1
ra: es la relación de área de corrientes definida como la relación del área
promedio de las corrientes de un orden i, sobre el área promedio de las
corrientes de orden inmediatamente inferior
Está definida como la relación, Dd, entre la longitud total a lo largo de todos
los canales de agua de la cuenca, en proyección horizontal, y la superficie
total de la hoya:
∑ li
Dd = (4.14)
A
donde:
Σli : longitud total de todos los canales de agua en km
A: área, en km2
li: longitud de cada cauce
Para las unidades citadas, se han encontrado valores mínimos de Dd del orden
de 7, valores promedios en el rango de 20 a 40 y valores máximos del orden de
400.
Valores bajos de Dd, por lo común, están asociados con regiones de alta
resistencia a la erosión, muy permeables, y de bajo relieve. Valores altos son
encontrados en regiones de suelos impermeables, con poca vegetación, de
relieve montañoso.
60
producen en la naturaleza.
61
régimen hidrológico, pues está relacionada con la erosión, temperatura y
evaporación de la región.
62
Se puede añadir que las propiedades geomorfológicas del subsuelo (como las
de los acuíferos), normalmente son parámetros que varían en las escalas de
tiempo geológico, y para el caso de la hidrología, pueden ser tomadas como
invariantes.
Además, algunas de las variables citadas con anterioridad son encontradas por
observaciones hechas sobre la cuenca y estimadas por medios estadísticos, y
desde tal punto de vista deben ser miradas. Aún más, en la definición de los
parámetros geomorfológicos, no se intenta dar la idea de relaciones biunívocas.
Como conclusión del análisis aquí considerado se puede decir que no existe
una relación única entre los parámetros físicos de la cuenca y las variables
hidrológicas, no obstante ellos pueden dar una orientación cualitativa en forma
y magnitud de las diferentes variables hidrológicas en el tiempo. Es claro que
en gran parte las características físicas de una cuenca son debidas a la acción
del agua y que por este hecho es factible pensar en la existencia de una relación
fuerte entre ellas a nivel determinístico, sin embargo,: la carencia de una
relación fuerte se debe fundamentalmente a la diferencia entre las escalas de
tiempo de los procesos dinámicos de la hidrología y a la geología. Además de
la fuerte componente estocástica de varios de los fenómenos hidrológicos.
63
variables y se puede aun llegar a rechazar un determinado ajuste. El ajuste de
la función se puede lograr mediante la técnica de regresión y correlación lineal
multivariada. Se puede obtener por ejemplo una función entre el caudal
máximo anual y algunos parámetros morfológicos como área, densidad de
drenaje, coeficiente de forma, etc. Además podría tenerse en cuenta la
precipitación entre las variables independientes así:
Qp = f(A,Dd, k, P). Con una función de este tipo y si se tiene una cuenca sin
datos de caudal se podría estimar el caudal y su intervalo de confianza.
Como la descarga tanto sólida como líquida de un río es una variable con un
marcado carácter aleatorio y como, además, la cuenca hidrográfica que lo
64
conforma presenta variaciones en los tipos y tamaños de suelos, vegetación,
etc., la predicción de los cambios futuros de una corriente no se puede definir
con patrones determinísticos (aunque existen leyes físicas que describen los
fenómenos locales). Así, el proceso morfológico de los ríos debe ser tomado al
menos con una gran componente de carácter aleatorio.
Cuando una corriente tiene los sedimentos que es capaz de transportar se dice
que el cauce está en equilibrio. Si se produce una sobrecarga de sedimentos
generada por cualquier causa, empieza un proceso de agradación o
sedimentación del lecho. En este caso, el río no tiene la suficiente energía para
transportar el material sólido que lleva y éste entonces se deposita en su cauce.
Si por el contrario hay una deficiencia, el fenómeno que se presenta es el de
degradación o erosión del lecho. En este caso el río tiene energía suficiente
para transportar el material sólido y para socavar el cauce.
65
hombre, se puede utilizar la conocida ecuación de Lane (1955):
Q bs Db
S∝ (4.15)
Qc
Donde:
S: pendiente del río
Qs :caudal sólido
D: diámetro del material del lecho
Q: caudal líquido.
a,b,c son exponentes que dependen del tipo de corriente.
66
FIGURA 4.10. Efectos de la construcción de una presa
4. Alineamiento artificial.
67
través de las relaciones entre caudal, carga de sedimentos, ancho, profundidad
y velocidad media. Leopold y Madodock (1953) desarrollaron las ecuaciones
principales de la geometría hidráulica; estas ecuaciones permiten hacer
inferencias sobre el comportamiento de un río, si se conocen algunas de las
variables; son las siguientes:
B = AQb (4.16)
V = KQm (4.17)
D = CQf (4.18)
Se define el caudal a sección llena como aquel caudal que fluye llenando el
cauce, sin derramar sobre las llanuras de inundación. La sección a banca
llena se ha definido de varias maneras:
68
• Si hay una llanura de inundación bien desarrollada, la altura de su
superficie puede ser considerada como la que determine el nivel de la
sección llena. Sin embargo, la definición de sección llena es mucho más
difícil si el cauce no está bien definido, por ejemplo, cuando las bancas
no tienen la misma elevación, en ríos trenzados, donde la diferencia entre
el cauce del río y la planicie de inundación no es tan obvia, y en
secciones complejas donde se presenten diferentes niveles de terrazas.
69
• En una zona de depósito a una altura incipiente de inundación, la altura
asociada con la parte más alta de los depósitos recientes (barras
puntuales y barras intermedias, pero no terrazas).
• Raíces expuestas por debajo de una capa de suelo intacta, que indican
exposición a los procesos erosivos.
1 2 1
Q= ⋅ A ⋅ R H3 ⋅ S f 2 (4.19)
n
donde:
Q : caudal en m3/s
A : área de la sección en m2
RH : radio Hidráulico en m
Sf : pendiente de la línea de energía
70
n : coeficiente de rugosidad de Manning
71
E = R × K × LS × C × P (4.20)
Donde:
E: pérdida de suelo calculada en ton/acre-año
R: índice de erosión. Depende de la energía cinética de las tormentas
máximas de 30 minutos de duración. En E.U hay mapas con valores de este
índice
72
Capítulo 5
REVISIÓN DE CONCEPTOS DE ESTADÍSTICA Y
PROBABILIDAD
73
La planeación y el diseño de proyectos relacionados con el agua necesitan
información de diferentes eventos hidrológicos que no son gobernados por
leyes físicas y químicas conocidas, sino por las leyes de azar. Por ejemplo,
el caudal de un río varía día a día y año tras año, y no puede predecirse
exactamente cual será su valor en un período de tiempo cualquiera. En el
caso del diseño de un puente, el estudio hidrológico determinaría la
creciente asociada con una probabilidad crítica(se busca determinar el caso
crítico), la cual se supone representa el riesgo para el puente. Esto solo
puede determinarse a través del análisis probabilístico y estadístico basado
en los registros hidrológicos del pasado.
Se dice que una variable aleatoria es discreta si ella sólo puede tomar
U U
74
puede asumir cualquier valor y es entonces una variable aleatoria continua
En este caso la ley de probabilidades asigna medidas de probabilidad a
rangos de ocurrencia de la variable aleatoria.
Por ejemplo supóngase que se tienen registros del caudal del río Magdalena
durante un período de 50 años. Son factibles dos tipos de análisis:
descriptivo y de inferencia. El primero se realiza sin ninguna referencia a su
población, de la cual se tiene una muestra de 50 años. Consiste,
básicamente, en calcular propiedades estadísticas, como media, varianza y
otras. En el segundo, la muestra se analiza para inferir las propiedades de su
población, lo cual ayudará a derivar las características probabilísticas del
caudal. El primero es una aplicación de los métodos estadísticos que
75
requieren poca decisión y poco riesgo. El segundo involucra riesgos y
requiere una total comprensión de los métodos empleados y el peligro
involucrado en la predicción y estimación de las variables.
m
P( X = x ) = (5.1)
n
76
1) La probabilidad de ocurrencia de un evento, Pi, siempre tiene un B B
0 ≤ Pi ≤ 1 (5.2)
.
La probabilidad de un evento cierto es 1:
∑P
i =1
i =1 (5.3)
entonces:
P ( X 1 ∪ X 2 ) = P( X 1 ) + P( X 2 ) (5.4)
P ( X1 ∪ X 2 ) = P ( X1 ) + P ( X 2 ) − P ( X1 ∩ X 2 ) (5.5)
77
La probabilidad de que dos eventos independientes ocurran de manera
simultánea es el producto de las probabilidades individuales así:
P ( X1 ∩ X 2 ) = P ( X1 ) × P ( X 2 ) (5.6)
X1 X1 ∩ X 2
P( ) = P( ) (5.7)
X2 P( X 2 )
Ejemplo 5.1
Supóngase que el río Cauca alcanza cada invierno un nivel de creciente con
una frecuencia relativa de 0.2. En el Cauca hay un puente cuya
probabilidad de falla en los estribos es 0,3 y la experiencia muestra que
cuando hay creciente, las probabilidades de esta falla suben a 0,5. Las
probabilidades son:
P(creciente) = P(C) = 0,2
P(no creciente) = P(C) = 0,8
P(falla) = P(F) = 0,3
P(no falla) = P(F) = 0,7
P (falla dada creciente) = P(F/C)= 0,5
Se desea conocer la probabilidad de falla del puente.
Solución:
El puente falla (queda inutilizado) cuando falla en los estribos o cuando hay
creciente; esto se puede denotar así:
78
P( C ∪ F ) = P( C ) + P( F ) − P( C ∩ F )
P ( C ∩ F ) = P( C ) × P ( F )
C
Reemplazando valores, se obtiene:
P ( C ∩ F ) = 0. 2 × . 0 . 5 = 0. 1
Al reemplazar este valor en la expresión de unión de probabilidades, se
concluye finalmente que P(C∪F)=0.4
79
una situación de riesgo, pues la obra se diseña para soportar cierta avenida
máxima , y crecientes mayores le podrían hacer daño o incluso destruirla. El
riego R puede entonces escribirse como:
1 n
R = 1 - (1 - ) (5.9)
Tr
La confiabilidad se define como el complemento del riesgo (Confiabilidad =
1-R). Se quiere que la obra tenga un riesgo pequeño de dañarse o, lo que es
lo mismo, una alta confiabilidad.
Ejemplo 5.2
1 25
R = 0.1 = 1 - (1 - ) T
Tr
Reemplazando los valores de Tr y n se obtiene: B B
TR = 238 años
B B
Ejemplo 5.3
Una presa por gravedad puede fallar por deslizamiento (A), por crecientes
(B), o por ambas. Asumir que :
1) La probabilidad de falla por deslizamiento es dos veces la
probabilidad de falla por creciente: P(A)=2 P(B)
80
2) La probabilidad de falla por deslizamiento, dado que ha habido creciente,
es 0.8
3) La probabilidad de falla de la presa es de 1*10-3 P
Solución:
La presa queda inutilizada cuando se presenta una falla por deslizamiento o
cuando hay una creciente, lo que puede expresarse como:
P( A ∪ B ) = 0.001 = P( A ) + P( B ) − P( A ∩ B ) (1)
Se sabe que:
A A∩B
P( ) = 0.8 = P( ) (4)
B P( B )
81
Ejemplo 5.4
fecales. Se pregunta:
a) Cuál es la probabilidad de que un sistema seleccionado al azar
resulte con contaminación fecal?
b) Suponiendo que un sistema se encuentre contaminado con materias
fecales, cuál es la probabilidad de que también esté contaminado con
plomo?
c) Cuál es la probabilidad de que un sistema seleccionado al azar esté
contaminado?
d) Suponiendo que la probabilidad de contaminación hallada en el
numeral anterior no es satisfactoria, y que se desea que no exceda de
0.01, ¿cuál es el valor permisible para la probabilidad de
contaminación por materias fecales, asumiendo que el valor de la
probabilidad condicional hallada en el numeral b aún se puede
aplicar?
Solución:
Llamemos P(F) a la probabilidad de contaminación por materia fecal, P(Pb) B B
P(Pb ∩ F)
P( Pb / F ) =
P(F)
82
y P(Pb) = 5/1000. Reemplazando, se obtiene que:
B B
P(PBI/F) = 2/17
c)Se pregunta en este numeral el valor de P(C); este valor establece la
probabilidad de que un circuito esté contaminado con plomo o con materias
fecales. Como hay 15 circuitos contaminados con materias fecales y 5
contaminados con plomo, se tiene entonces que:
P(C) = 20/1000= 0.002
d) La probabilidad de contaminación C se puede expresar como:
P( C) = P( F ∪ Pb ) − P( F ) + P( B ) − P( F ∩ Pb ) (1)
y se conoce el valor de la probabilidad condicional:
P( Pb / F ) = 2 / 17 = P( Pb ∩ F ) (2)
P( F )
Resolviendo la (1) y la (2) simultáneamente se halla que:
P(F) = 0.00567
se representa por una letra mayúscula la variable aleatoria, y por una letra
minúscula, un valor específico, una relación o una muestra de la variable.
P(X = a) indica la probabilidad de que la variable aleatoria X tenga un valor
de a; similarmente, P(a<X<b) indica la probabilidad que la variable
83
aleatoria X esté en el intervalo [a, b] .Si se conoce la probabilidad P(a<X<b)
para todos los posibles valores de a y b, se dice que se conoce la
distribución de probabilidades de la variable X.
84
ni
f s ( xi ) = (5.10)
n
la cual es un estimado de P( xi -∆x<X<xi), la probabilidad de que la variable
B B B B
F S ( x i ) = ∑ f S ( x j) (5.11)
j =1
acumulada de probabilidades.
de probabilidades fX(x)
U U U U B B
(x)
f X (x) = lim f S (5.12)
n→∞
∆x → 0
∆x
La función de frecuencia acumulada se convierte en la función acumulada U U U U U U
85
cuya derivada es la función de densidad de probabilidad:
dFX (x)
f X (x) = (5.14)
dx
P(X ≤ x) = FX (x) =
−∞
∫f X (u)du (5.15)
FX (xA ) = P (X ≤ x A ) (5.16)
86
Donde µ representa la media y σ es la desviación típica de la variable
hidrológica. XT es el valor de la variable aleatoria asociada a un ‘período de
B B
1
0.9
0.8
0.7
0.6
FX(x)
0.5
0.4
0.3
0.2
0.1
0
1 2 3 4 5 6 7 8 9 10
x
De donde:
1
FX ( X T ) = 1 −
T
87
O:
1
FX ( µ + σK ) = 1 −
T
Y se obtiene finalmente:
1 −1 ⎛ 1⎞
K= FX ⎜ 1 − ⎟ − µ
σ ⎝ T⎠
∫f
-∞
X (x)dx = 1 (5.19)
2)
b
P(a ≤ X ≤ b) = ∫ f X (x)dx (5.20)
a
3)
b (5.21)
∫ f X (x)dx = 0
b
88
1)
(5.22)
∑ f (x ) = 1
i
i
2)
xi ≤b
P( a ≤ X ≤ b ) = ∑ f (x )
xi ≥a
i (5.23)
3)
i= j
P( X ≤ x j ) = ∑ f ( x i ) (5.24)
i =1
Lo que implica que las probabilidades se definen solo como áreas bajo la
función de densidad de probabilidades, FDP, entre límites finitos.
Ejemplo 5.5
Hallar la función de distribución acumulada para una variable aleatoria que
se define como el número de veces que se lanza una moneda, hasta que
aparece cara.
Solución:
La probabilidad de que caiga cara en cualquier ensayo es ½ y es
independiente de la probabilidad de que caiga sello.
Si A es el evento de que caiga sello en el primer ensayo y B (es el evento)
de que caiga sello en el segundo ensayo, la probabilidad que suceda A y B
es:
89
P(AB) = P(A) + P(B) = (1/2)2 P P
(1/2)x-1 ½ = (1/2)x
P P P P
1 ½ ½
2 ¼ ¾
3 1/8 7/8
probabilidades acumulada.
µr′ = ∫x
-∞
r
f X (x)dx (5.25)
90
o en el caso discreto:
n
µr′ = ∑ xri f X ( xi ) (5.26)
i =1
µr = ∫ (x - µ ) f X (x)dx
r
(5.27)
-∞
n
µ r = ∑ f X ( x i )( x − µ ) r (5.28)
i =1
Rara vez se necesita calcular más de tres momentos. Estos son usados para
estimar los parámetros y describir las características de la distribución.
91
En general, las características estadísticas básicas se calculan como el valorU U U U
∞
E[g( X )] = ∫ g( u )f X ( u )du (5.29)
−∞
(FDP) de la variable X
Las principales características son:
- La media µ: representa el valor esperado de la variable misma. Para
una variable aleatoria X, la media E(X) es el primer momento con
respecto al origen; es una medida de la tendencia central de la
distribución:
∞
E(X) = µ = ∫ xf
-∞
X (x)dx (5.30)
1 N
µˆ x = ∑ xi (5.31)
N i =1
92
∞
E[(X - µ ) ] = σ 2 = ∫ (x - µ ) f X (x)dx
2 2
(5.32)
-∞
∑
2 2
σˆ x = ( x i − µ x ) (5.33)
N - 1 i =1
∞
E[(X - µ ) ] = ∫ (x - µ ) f X (x)dx
3 3
(5.35)
-∞
93
FIGURA 5.3 Distribución de probabilidades con diferente desviación
estándar.
La asimetría se hace adimensional dividiendo la anterior ecuación
por σ3 y se obtiene así, el coeficiente de asimetría γ:
P P
1 3
γ = 3
E[(x - µ ) ] (5.36)
σ
∧ i =1 (5.37)
γx = 3
(N - 1)(N - 2) σ̂ x
Como se muestra en la Figura 5..4, para γ>0, asimetría positiva, los datos
se concentran a la derecha y para γ<0, asimetría negativa, los datos se
concentran a la izquierda.
94
γ <0 γ >0
fX (x)
x
µ
Ejemplo 5.6
Frecuencia
Intervalo en mm Absoluta
100-110 10
110-112 16
120-130 9
130-140 10
140-150 20
150-160 15
160-170 20
Solución:
95
En total se tiene 100 valores, para cada intervalo se halla el valor medio o
marca de clase y se le asigna una frecuencia relativa, la cual es la frecuencia
absoluta sobre el número total de valores (100). El valor medio de cada
intervalo es xi y la frecuencia relativa es fx(xi).
B B B B B B
(mm) fx(xi)
B B B B
Σ=100 Σ=138.90
96
un fracaso, esto es, solo hay dos eventos mutuamente excluyentes para un
experimento. Además, los experimentos sucesivos son independientes y la
probabilidad permanece constante de ensayo a ensayo. Un ejemplo en
hidrología sería la probabilidad de que un día sea lluvioso o seco. La
distribuciones de este tipo más usadas en hidrología son la distribución
binomial y la geométrica.
n!
P( X = x ) = p x (1 − p )n − x (5.38)
x! (n − x)!
97
µ = np
σ 2 = np(1 − p )
(5.39)
(q − p )
γ=
npq
Ejemplo 5.7
98
P[una o más crecientes en 50 años] = 1 - P[no crecientes en 50 años] o lo
que es lo mismo:
50
P[una o más crecientes en 50 años]= 1 − (1 / 50) 0 (1 − 1 / 50) 50 = 0.64
0
5.5.2 Distribución Geométrica.
año,P, es:
P = (1 − p )n −1 p (5.40)
1
µ=
P
(1 − P ) (5.41)
σ2 =
P2
Ejemplo 5.9
99
b) Cuál es la probabilidad de que durante los próximos 20 años la altura
hallada en el numeral anterior sea excedida al menos una vez?.
c) Cuál es la probabilidad de que durante los próximos 5 años este valor
sea excedido exactamente una vez?
F(H)
5 6 7
H(m)
Solución:
P( H ≥ H Tr=20 ) = 1 / 20 = 0.05
100
(7 − H Tr=20 )(0.95)
0.05 =
2
Despejando el valor de H, se obtiene finalmente:
H Tr=20 = 6.9 m
5
P( H Tr = 20 = 1) = (0.05)1 (0.95) 0.4 = 0.024
1
Ejemplo 5.9
101
a) Cuál es la probabilidad de inundación en un año cualquiera producida
solamente por el río A.
b) Cuál es la probabilidad de inundación de la planicie en un año?
c) Cuál es la probabilidad de que no haya inundación en los próximos 4
años?
Solución:
102
P( A ∪ B ) = 0.145 + 0.04 − 0.145 × 0.04 = 0.179
c) La probabilidad de inundación, P, en cualquier año, es 0.179, como se
explicó en el numeral anterior, y la probabilidad ,q, de no inundación
será entonces:
q =1 -P =1 - 0.179 =0.821
y la probabilidad de no inundación en 4 años será entonces:
P(no inundación en 4 años) =(0.821)4 =0.454
P P
Ejemplo 5.10
103
La mayoría de las variables hidrológicas son variables aleatorias continuas.
Enseguida se describen brevemente las distribuciones de probabilidades más
usadas en análisis de frecuencia de estas variables.
( x −µ x )2
1 −
f X (x) = e 2σ 2 x
(5.42)
σ x 2π
104
µ = (x − µ x ) / σ x
se obtiene como FDP y como función acumulada de la variable µ:
1 - u2
f u (u) = e 2
2π
∞ (5.43)
1 2
-w
Fu (u) =
2π
∫e
-∞
2 dµ
105
- Aproximar la distribución de probabilidades de errores aleatorios .
106
5.6.1.1 Estimación de parámetros
1 N
µˆ = ∑ xi
N i =1
(5.44)
N
1
σˆ = = ∑ ( x i − µˆ )1 / 2 (5.45)
N
x - µ̂
K= (5.46)
σ
107
⎛ 1⎞
2. Usando el valor calculado de ⎜1 − ⎟ en la tabla 5.1, se lee el valor
⎝ T⎠
P P
X T = µˆ + Kσˆ
Ejemplo 5.11
Se tiene una estación con 30 años de datos de caudales medios anuales con
media de 117 m3/s y desviación estándar de 94 m3/s. ¿Si los datos se ajustan a
P P P P
Solución:
K = Fu-1 (0.99) B PB P
108
valor estimado de la media, y si es pequeña, habrá, por el contrario, mucha
confianza en ese valor estimado. Con ese fin se utilizan los llamados intervalos
de confianza.
XT ± u1-α ST (5.47)
109
en donde ; es el nivel de probabilidad, u1-α es el cuantil de la distribución
B B
estándar.
Cada distribución tiene expresiones para hallar el error estándar, por ejemplo,
el de la distribución Normal es:
σˆ x
(1 + K 2 / 2)2
1
ST = (5.48)
N
Ejemplo 5.12
Los caudales medios anuales de un río con media 1.5 m3/s y desviación P P
Solución:
1 − µˆ
P( X ≤ 1) = P(µ ≤ )
σˆ
1 − 1 .5
P (µ ≤ ) = P(µ ≤ −0.83)
0 .6
110
FIGURA 5.6 Simetría de la distribución normal.
Ejemplo 5.13
Solución:
280 − 356
P ≤ 280) = P(µ ≤ ) = P(µ ≤ −0.997)
76.2
y:
111
contribuyente, pérdidas, coeficiente de evaporación, etc. En general, cuando
la variable aleatoria X es el producto de un gran número de otras variables
aleatorias, la distribución de los logaritmos de X puede aproximarse a la
Normal, ya que los logaritmos de X son la suma de los logaritmos de los
factores contribuyentes. Si se tiene una variable aleatoria X y ln X = Y se
ajusta a una distribución Normal, se dice que la variable aleatoria X es
lognormalmente distribuida.
f X (x) =
1
exp ⎢ -
(
⎡ 1 y - µy )
2
⎤
⎥
2 (5.49)
σ y x 2π ⎢⎣ 2 σy ⎥⎦
112
Se ha demostrado que la distribución lognormal puede aplicarse en un amplio
número de eventos hidrológicos, especialmente a aquellos casos en los cuales
la variable tiene un límite inferior, la distribución empírica no es simétrica y
los factores que causan los eventos son independientes y multiplicativos.
f X (x) =
1 ⎧⎪ 1 ln (X - xo ) - µ y
exp ⎨ -
[ ]
2
⎫⎪
⎬ (5.50)
2π (X - xo ) σy ⎪⎩ 2 σy ⎪⎭
donde los parámetros µy, σy y xo son llamados los parámetros de escala, forma
B B B B B B
y localización respectivamente.
con dos parámetros solo puede aplicarse a eventos con valores positivos.
1 N
µˆ Y = ∑ log a ( X i )
N i =1
(5.51)
12
⎧⎛ 1 ⎞ N 2⎫
σˆ Y = ⎨⎜ ⎟ ∑ [log a ( X i ) − µˆ Y ] ⎬ (5.52)
⎩⎝ N ⎠ i =1 ⎭
113
positivo. En este método, el segundo momento de Z = X - xo no depende de B B
x0, esto es, σ²z = σ²x y µz = µx - x0, entonces el límite inferior xo se puede
B B B B B B B B B B B B B B
expresar como:
⎛ Cv x ⎞
x0 = µ x ⎜⎜ 1 - ⎟⎟ (5.53)
⎝ Cv z ⎠
Donde:
σx
Cv x =
µx
(5.54)
σ
Cv z = z
µz
Donde:
(1- w ) 2/3
Cvz = 1/3
w
(5.55)
1
[
w = - γˆ x + γˆ 2x + 4
2
( )
1/2
]; γx > 0
1 2
µ Y = log a (µ X ) − σ Y (5.56)
2
1
⎡ ⎛ σ ⎞⎤
2 2
σ Y = ⎢log a ⎜⎜ 1 + X 2 ⎟⎟ ⎥ (5.57)
⎢⎣ ⎝ µ X ⎠ ⎥⎦
114
En este caso, se estiman µX y σX con los datos originales, y con las
B B B B
lognormal.
Ejemplo 5.14
Calcular la probabilidad de que el caudal medio esté entre 300 y 400 m3/s P P
Solución:
P(300[Q[400)= FX(400)-FX(300)
B B B B
⎛ 400 - µ x ⎞ ⎛ 300 − µ x ⎞
P(300≤Q≤400)= Fµ ⎜⎜ ⎟⎟ − Fu ⎜⎜ ⎟⎟
⎝ σx ⎠ ⎝ σx ⎠
B B
= Fu (u400) - Fu (u300)
B B B B B B B B
donde:
y
u400 = (400 - 256.7)/191 = 0.75
B B
115
lo que implica que:
P(300≤Q≤400)=0.7734 - 0.5871=0.1863
P(300≤Q≤400)=FY(ln(400))-FY(ln(300))
B B B B
⎛ ln(400) − µ Y ⎞ ⎛ ln(300) − µ Y ⎞
= Fu ⎜⎜ ⎟⎟ − Fu ⎜⎜ ⎟⎟
⎝ σY ⎠ ⎝ σY ⎠
y:
ln(300) = 5.704
ln(400) = 5.99
ln (X T ) = µ y + K σ y (5.58)
116
en donde K = Fu −1 ⎛⎜ 1 − ⎞⎟
1
⎝ T⎠
Si se quiere trabajar con la variable no transformada al campo logarítmico se
tiene que:
⎡ 1/2 ⎡ ln (1 + Cv )⎤ ⎤
2
exp⎢K T (ln(1 + Cv 2 )) - ⎢ ⎥⎥ - 1
⎣ ⎣ 2 ⎦⎦ (5.59)
K=
Cv
donde:
-1 ⎛ 1⎞
K T = Fu ⎜⎜ 1 - ⎟⎟ (5.60)
⎝ Tr ⎠
−1 ⎛ 1⎞
Fu ⎜ 1 − ⎟ es el inverso de la función de distribución Normal estandarizada
⎝ T⎠
acumulada y Cv es el coeficiente de variación
B B
en donde:
σY
ST = δ (5.62)
N
117
1/2
⎛ 2
⎞
δ = ⎜⎜ 1 + K T ⎟⎟ (5.63)
⎝ 2 ⎠
Ejemplo 5.15
Solución:
1 1
Fu (K T ) = 1 - = 1- = 0.99
TI 100
De la tabla 5.1:
K T = Fµ−1 (0.99) = 2.33
⎡ 1/2 ⎡ ln (1 + 0. 33 )⎤ ⎤
2
exp⎢ 2.33 (ln(1 + 0. 332 )) - ⎢ ⎥⎥ - 1
⎣ ⎣ 2 ⎦⎦
K=
0.333
118
K= 3.028
El valor asociado a un período de retorno de 100 años será:
el resultado es:
1/2
⎛ 2
⎞
δ = ⎜⎜ 1 + 2.33 ⎟⎟ = 1.93
⎝ 2 ⎠
0.3246
ST = 1.93 * = 0.11
30
Por lo tanto:
ln (30.28) ± 1.64 * 0.11
= 3.41 ± 0.1875
= [3.2225, 3.5975]
= [e3.2225, e3.5975] = [25.091, 36.5]
P P P P
119
Ellas difieren entre sí por el valor del parámetro de forma. La expresión
general de la función de densidad de probabilidades para la distribución
extrema tipo I o Gumbel es:
1 ⎡ x -β ⎛ x - β ⎞⎤
f X (x) = exp⎢- - exp⎜ - ⎟⎥ (5.64)
α ⎣ α ⎝ α ⎠⎦
6
αˆ = σˆ (5.65)
π
βˆ = µ - 0.5772αˆ (5.66)
6
K=- {0.577 + ln[lnTr - ln(Tr - 1)]} (5.67)
π
120
Los límites de confianza por el método de momentos para un nivel de
probabilidad ; son:
X T ± u 1- α 2 S T (5.68)
σ
ST = δ (5.69)
N
δ = [1 + 1.1396K + 1.1 K 2]
1/2
(5.70)
121
parámetros. La última función es llamada también Distribución Pearson tipo
III. La distribución Gamma está relacionada con otras distribuciones muy
conocidas como las distribuciones Chi-cuadrado y la exponencial negativa,
que son casos particulares de la distribución Gamma.
La distribución Gamma de dos parámetros tiene una función de densidad de
probabilidades de la forma:
β -1
1 ⎛ x ⎞ - αx
f X (x) = ⎜ ⎟ e (5.71)
| α | Γ(β ) ⎝ α ⎠
Donde:
El parámetro β siempre es mayor que cero, mientras que ; puede ser positivo o
negativo. La función Gamma completa está dada por:
∞
Γ(β) = ∫ zβ-1 e-z dz (5.72)
0
β -1
1 ⎛ x - xo ⎞ ⎛ x - xo ⎞
f X (x) = ⎜ ⎟ exp⎜ - ⎟ (5.73)
| α | Γ (β ) ⎝ α ⎠ ⎝ α ⎠
Donde:
xo ≤ x < ; para ; > 0
B B
122
-: < x ≤ xo para ; < 0
B B
parámetro de localización.
La Figura 5.9 muestra formas de la función de densidad de probabilidades
Gamma para ; > 0.
µ = αβ (5.74)
σ 2 = α 2β (5.75)
123
1
βˆ = 2
Ĉ v
µˆ (5.76)
αˆ =
βˆ
2
⎛ 2⎞
βˆ = ⎜⎜ ⎟⎟ (5.77)
⎝ γˆ ⎠
γˆ
αˆ = σˆ (5.78)
2
ˆ
X̂0 = µˆ − αˆ β (5.79)
γ$ es el coeficiente de asimetría calculado usando la muestra.
Si se define:
⎛ 1⎞
K T = Fu ⎜⎜ 1 - ⎟⎟ (5.80)
⎝ Tr ⎠
el factor de frecuencia K tiene la siguiente forma:
124
2 3 4
2 γˆ 1 3 ⎛ γˆ ⎞ 2 ⎛ γˆ ⎞ ⎛ γˆ ⎞
K ≈ K T + (K t − 1) + (K T − 6K T ) ⎜ ⎟ − (K T − 1) ⎜ ⎟ + K T ⎜ ⎟
6 3 ⎝6⎠ ⎝6⎠ ⎝ 6 ⎠ (5.81)
Si se tiene que:
XT ± u 1− α 2 S T
σ
ST = δ (5.82)
N
β -1
1 ⎡ ln(x) - y o ⎤ ⎡ ln (x)- y o ⎤
-⎢
f x (x) = e ⎥ (5.83)
x α Γ(β) ⎢⎣ ⎥ ⎣ α ⎦
α ⎦
parámetro de localización.
125
126
TABLA 5.2. VALORES DE KT PARA LA DISTRIBUCIÓN PEARSON
B B
127
FIGURA 5.10 Distribución Log-Pearson Tipo III. (Salas, 1992).
128
YT = ln XT = µˆ y + K σˆ y (5.85)
γ Tr=2
B B Tr=5
B B Tr=10
B B Tr=20
B B Tr=50
B B Tr=100
B B
129
σˆ y
ST = δ (5.86)
N
Los límite de confianza se pueden expresar como:
ln X T ± µ 1− α / 2S T (5.87)
California:
m
P= (5.88)
n
130
Weibull:
m
P= (5.89)
n +1
Hazen:
2m -1
P= (5.90)
2n
Ven te Chow propuso que toda muestra se puede ajustar a una expresión
como la siguiente:
X = µˆ + K σˆ (5.91)
131
132
133
5.8 BONDAD DE AJUSTE DE UNA DISTRIBUCION DE
PROBABILIDADES
Por ejemplo, las distribuciones Pearson tipo III y Log-Pearson tipo III
requieren la estimación del coeficiente de asimetría de datos muestrales. Esto
puede ser una razón suficiente para preferir cualquier otra distribución, ya que
este parámetro tiene un comportamiento muy sesgado, por lo cual se
necesitaría una gran cantidad de registros para tener un estimado más o menos
confiable, y dichos registros no se consiguen fácilmente en nuestro medio. Por
otra parte, las distribuciones de dos parámetros tienen un valor fijo o ignoran
la asimetría de la población, lo cual tampoco es conveniente.
134
El estadístico Smirnov - Kolmogorov, D, considera la máxima desviación de
la función de distribución de probabilidades empírica de la muestra, FE(x), de
la función de distribución de probabilidades teórica, escogida , Fx (x), tal que: B B
Supongase que en una muestra se tengan una serie de posibles eventos E1, E2, B B B B
....Ek que ocurren con frecuencias observadas de O1, O2, .....Ok. Si se tiene
B B B B B B B B
135
N ;=0.20 ;=0.10 ;=0.05 ;=0.01
5 0.45 0.51 0.56 0.67
10 0.32 0.37 0.41 0.49
15 0.27 0.30 0.34 0.40
20 0.23 0.26 0.29 0.36
25 0.21 0.24 0.27 0.32
30 0.19 0.22 0.24 0.29
35 0.18 0.20 0.23 0.27
40 0.17 0.19 0.21 0.25
45 0.16 0.18 0.20 0.24
50 0.15 0.17 0.19 0.23
N∃50 1.07 1.22 1.36 1.63
N N N N
k
⎡ (O − e ) 2 ⎤
χ2 = ∑ ⎢ i i ⎥ (5.93)
i =1 ⎣ ei ⎦
donde:
∑ Oi = ∑ ei
Si ξ2 =0, significa que las distribucion teórica y empírica ajustan exactamente,
P P
significancia ;de 0.05 o 0.01 ( el nivel de confianza se define como 1-;, siendo
136
frecuentemente utilizados niveles de confianza del 95%), se puede decir que
las frecuencias observadas difieren significativamente de las frecuencias
esperadas y entonces la hipótesis Ho se rechaza (para esos niveles de
B B
Ejemplo 5.16
Solución:
137
Año Junio Julio Agosto
1944 77 77 77
1945 72 76 76
1946 76 78 74
1947 74 74 83
1948 78 80 76
1949 75 79 74
1950 75 73 70
1951 73 78 78
1952 82 81 77
1953 79 80 78
1954 78 83 80
1955 69 80 79
1956 74 77 77
1957 75 76 74
1958 72 76 74
1959 72 75 76
138
T fT(t)
B B FE(t) FE(t)N FT(t)
B B FE(t)N -FT(t)
B B
68 0 0 0 0.0035 0.0045
69 1 1 0.0208 0.0084 0.0124
70 1 2 0.0417 0.0197 0.022
71 0 2 0.0417 0.0409 0.0008
72 3 5 0.1042 0.0778 0.0264
73 2 7 0.1458 0.1357 0.0101
74 7 14 0.2917 0.2206 0.0711
75 4 18 0.3750 0.3264 0.0486
76 7 25 0.5208 0.488 0.0328
77 6 31 0.6458 0.5753 0.0705
78 6 37 0.7708 0.6950 0.0758
79 3 40 0.83333 0.7995 0.0338
80 4 44 0.9167 0.8770 0.0397
Ejemplo 5.17
139
Se tienen los caudales máximos instantáneos de la estación RP-3 en el Río
Murrí, en el departamento de Antioquia. Se desea encontrar el caudal de un
período de retorno de 50 años hallado con las distribuciones Gumbel,
Lognormal de dos parámetros y Log Pearson tipo III.
Año Q m3/s
P P
1978 3239.0
1979 3431.7
1980 4577.9
1981 3612.0
1982 4151.8
1983 1949.0
1984 2342.9
1985 1345.0
1986 1862.2
1987 1652.8
1988 4220.0
1989 4958.4
1990 2664.9
1991 1392.7
Solución
Distribución Gumbel
140
De la ecuación 5.67 se halla el factor de frecuencia K=2.5924
que:
ST=1111.458 m3/s
B B P P
(4329.37 ≤QTr=50=6158≤7986.07) B B
Distribución Log-Normal
ln QTr=50=8.8286 B B
y sacando el antilogaritmo :
QTr=50=6827 m3/s
B B P P
141
Para un ;=0.05 se obtiene de la tabla 5.1 T0.95=1.64. Finalmente :
B B
(4814.4≤QTr=50=6827≤9679.84)
B B
Se tiene que:
µˆ = 2957.2
σˆ = 1234.6
γˆ = 0.1702
K=2.144
QTr=50=5604 m3/s
B B P P
Con la ecuación 5.82 y con la tabla 5.3 se obtiene un error estandar ST=809.05
B B
(4273≤QTr=50=5604≤6934.9)
B B
142
Capítulo 6
LA PRECIPITACIÓN
142
6.1 INTRODUCCION
143
condiciones:
Precipitación orográfica.
144
La masa de aire se encuentra con una barrera y es obligada a ascender,
siguiendo los accidentes naturales del terreno, tales como las montañas. Por lo
general, el lado de la montaña contra el que choca el viento es la zona lluviosa,
mientras el otro lado es más seco..
Por ejemplo, cuando una masa de aire frío procedente de los polos se encuentra
con una masa de aire caliente, estas dos masas no se mezclan y forman una
discontinuidad: la masa de aire frío, más densa, se sitúa debajo de la de aire
caliente. Cuando una masa de aire se empieza a mover, su posición anterior es
ocupada por un frente. Un frente cálido se forma cuando aire caliente
reemplaza el aire frío y un frente frío se forma cuando el aire frío desplaza las
masa de aire caliente.
Los frentes se extienden grandes áreas, a veces hasta de más de 3000 km2
(Mutreja,1986)
145
La precipitación se mide generalmente con pluviómetros, que son recipientes
estandarizados en los cuales puede medirse la lámina precipitada. La figura 6.1
muestra un esquema de este tipo de mecanismos.
146
marcar en un tipo especial de papel la variación de la precipitación con el
tiempo, como se muestra en la figura 6.2.
La precipitación también se estima por medio de fotos de satélite; el color y la
forma de las nubes permiten a los expertos estimar la cantidad de agua
precipitada que aquellas podrían producir. Los radares también permiten hacer
estimaciones sobre la cantidad de lluvia que produciría una masa de nubes: el
radar emite ondas electromagnéticas y la velocidad con que las nubes reflejan
estas ondas depende del tamaño y densidad de las gotas de agua de la nube.
147
Es un gráfico de la precipitación acumulada contra el tiempo, en orden
cronológico. Es la curva que se obtiene directamente del pluviógrafo, ver
figura 6.3. Las curvas de masa se usan para extraer información sobre la
magnitud, duración e intensidad de una tormenta.
10
9
8
Precipitación (mm)
7
6
5
4
3
2
1
0
1
3
5
7
9
11
13
15
17
19
21
Tiempo (h) 23
6.6.2 Hietograma.
148
tormentas, para el estudio de caudales máximos, y se deriva de la curva de
masa. El área bajo el hietograma representa la precipitación total recibida en
ese período
Precipitación acumulada (mm)
180
160
140
120
100
80
60
40
20
0
1
11
13
15
Tiempo (días)
149
Antes de iniciar los estudios de los registros de lluvia de una estación
determinada, es necesario verificar la continuidad y la consistencia de éstos. La
continuidad puede romperse porque faltan datos, debido a daños en el
pluviómetro; falta de recolección de los registros, etc.
Los registros faltantes se pueden estimar usando los de las estaciones vecinas,
utilizando la precipitación normal como estandar de comparación. La
precipitación normal es el promedio de la precipitación anual, mensual o diaria
cuando se tienen una longitud de registros de al menos 30 años.
150
Nx ⎡ P1 P2 P ⎤
Px = ⎢ + + ........ m ⎥ (6.1)
M ⎢⎣ N 1 N N m ⎥⎦
450,00
Precipitación promedio
400,00
350,00
mensual (mm)
300,00
250,00
200,00
150,00
100,00
50,00
0,00
1 2 3 4 5 6 7 8 9 10 11 12
Meses
151
Se escogen N estaciones vecinas a la estación X que se analiza. Los registros
de la precipitación media acumulada anual de la estación X se comparan con
los registros de las precipitaciones medias anuales acumuladas de las
estaciones vecinas, figura 6.6 . Si se observa un cambio de pendiente, esto
indica un cambio en el régimen de la estación x. Los valores de la estación X,
a partir del cambio, (año 65 en la figura) se corrigen usando la relación:
Mc
Pcx = Px (6.2)
Ma
donde:
152
+ + + ... Pn
P = P1 P 2 P 3 (6.3)
n
5
4 .5
4
la estación X en mm *10
1971
4
3 .5 1969
3 a
1967
2 .5
c
2 1965
1963
1 .5
1961
1
0 .5 1959
1957
0
1
11
13
15
P r e c ip it a c ió n a n u a l a c u m u la d a e n e s t a c io n e s
4
v e c in a s e n m m * 1 0
153
n
∑P
i =1
i Ai
(6.4)
P=
AT
donde:
n: Número de estaciones usadas en el análisis.
AT: área total de la cuenca
Ai: área de influencia de la estación i
Este método determina las áreas de influencia usando únicamente un criterio
geométrico, sin tener en cuenta influencias climáticas o topográficas.
6.8.3 Isoyetas
Las isoyetas son las líneas que unen los puntos de igual precipitación. Para la
aplicación de este método, se dibuja la cuenca a escala y se ubican las
estaciones de precipitación con sus valores respectivos.
154
isoyetas, el valor promedio de la precipitación para un área A será:
P1 + P2 P + P3 P + Pn
a1 ( ) + a2 ( 2 ) + ......a n −1 ( n −1 )
2 2 2 (6.5)
P=
A
Suponiendo que se tengan las isoyetas, tal como se muestra en la Figura 6.8,la
precipitación promedio será:
Este método, permite, si la persona que lo está usando conoce el área tener en
cuenta variaciones locales de la precipitación, topografía, etc. Sin embargo, en
regiones montañosas tropicales, como son las de los Andes colombianos, para
aplicar con éxito esta metodología es necesario contar con
un buen número de estaciones, pues la precipitación varía con la altura, en
distancias muy cortas.
155
La figura 6.9 muestra para una zona del Departamento de Risaralda la
variación de la precipitación con la altura En ésta se tienen variaciones de
altura desde los 1100 msnm a orillas del río Cauca hasta 5500 msnm en las
cumbres del Parque de Los Nevados. Se observa que la precipitación
promedia anual aumenta con la altura hasta aproximadamente los 2000 msnm
y a partir de esta cota empiezan a disminuir los valores promedios anuales.
2500
2000
Altura (m)
1500
1000
500
1600 1800 2000 2200 2400 2600 2800 3000 3200 3400
Precipitación (mm/año)
156
Nacional 1997)
114000
112000
0
110000
SUPÍA
PENSILVANI
108000
MISTRAT
Ó
106000
MANIZALE
S
104000
MARSELL
STA.
DOSQUEBRADA
PEREIR
A
102000
ESCAL
A
0 10 20 30 km
100000
78000 80000 82000 84000 86000 88000 90000 92000 94000
Nota: al oriente de la línea punteada la información no es confiable por escasez de estaciones
157
ingeniero necesita conocer el área que drena al punto de control. Las curvas
profundidad- área- duración relacionan láminas de lluvia y áreas, para
tormentas de diferentes duraciones. Para una lluvia de duración dada, la
lámina promedio decrece con el área en forma exponencial, así:
n
P = Po e-K A (6.6)
donde:
P::lámina, en cm
A:área en, km2
Po: registro más alto de lluvia en el centro de la tormenta
K y n:constantes para una región dada.
158
6.10 CURVAS INTENSIDAD - FRECUENCIA - DURACION
m
kTr
i= (6.7)
(c + d ) n
donde:
i: intensidad, en mm/h
Tr: período de retorno, en años.
d: duración, en minutos
k, m, c y.n: son parámetros que dependen de la zona donde esté ubicado el
pluviógrafo.
159
Para obtener las curvas son indispensables los registros pluviográficos. El
procedimiento para construir las curvas es el siguiente:
Lluvia 15 min
Marzo 15 de 1967 26 mm
Octubre 22 de 1967 29 mm
Febrero 10 de 1968 12 mm
3) A cada precipitación se le calcula la intensidad en mm por hora. Por
ejemplo para las de 15 minutos, el 1de Marzo de 1967:
26 60 mm mm
i= = 104
15 min hora h
m
P=
n+1
160
permitirá asociar intensidades a períodos de retorno.
7) Se construye la curva.
440
420
Periodo
400
de retorno
100 años
380
50 años
360
25 años
340
320 10 años
300 5 años
Intensidad [ mm / h ]
260
240
220
200
180
160
140
120
100
80
60
40
20
0
0 10 20 30 40 50 60 70 80 90 100 110 120
Duración [ minutos]
161
de los métodos que existen para ello, es tratar de hallar coeficientes regionales
para la ecuación 6.7. En el Departamento de Antioquia se aplicaron este y otro
procedimiento para hallar las curvas I-D-F en sitios sin información
pluviográfica. Una completa descripción de la metodología puede hallarse en
Smith,Vélez (1997)
162
ajustan a comportamientos como los hallados por Huff, razón por la cual,
éstos deben usarse con cuidado en zonas tropicales como Colombia.
100
90
80
PROBABILIDAD
70 DE
% Precipitacion
EXCEDENCIA
10%
60 30%
50%
50 70%
90%
40
30
20
10
0
0 20 40 60 80 100 120
% Tiempo
Los valores extremos de eventos como lluvias y caudales han sido objeto de
estudio de los hidrólogos durante muchos años debido a sus efectos
catastróficos sobre el entorno humano Por esta razón, el diseño de algunas
163
estructuras hidráulicas, tales como vertederos en presas, se debe realizar con
caudales cuyas probabilidades de falla sean de casi cero. Para tales
estructuras, la máxima precipitación que puede esperarse en esa zona se usa
para calcular caudales máximos. Esto significa, que hay un límite máximo de
precipitación que puede físicamente caer sobre una región en un tiempo dado.
Esta es la llamada precipitación máxima probable, PMP.
Desde este punto de vista, la PMP puede definirse como la precipitación que
produce sobre una cuenca, con riesgo cero de ser excedida. Se usan dos tipos
de metodologías para hallar la PMP: La primera usa métodos estadísticos
aplicados a los registros de lluvias extremas y el segundo método estudia los
mecanismos físicos que producen las máximas tormentas. Este último método
es dominio de la meteorología. Detalles de esta metodología se pueden
consultar en Weisner, C.J. (1970).
PMP = P + K σ (6.8)
donde:
XMax :es la lluvia máxima observada de 24 horas.
164
&P&& :promedio de la precipitaciones máximas observadas.
σ :desviación estándar de las precipitaciones máximas observadas.
165
Capítulo 7
ANÁLISIS DE CAUDALES
165
7.1 GENERALIDADES
166
alarmante la falta casi total de estaciones de medida en las áreas urbanas y
semirrurales de la mayoría de ciudades colombianas, ocasionando un completo
desconocimiento del comportamiento hidráulico de pequeñas corrientes,
responsables muchas veces de inundaciones y tragedias en las épocas
invernales.
Métodos directos:
Método área velocidad
Dilución con trazadores
Métodos indirectos:
Estructuras hidráulicas.
Método área pendiente.
Con muy pocas excepciones, las medidas continuas de caudal en el tiempo son
muy costosas , por lo que se relaciona el caudal con el nivel del agua, que se
puede medir mucho más fácil que el caudal. Las curvas que relacionan estos
niveles con el caudal son las llamadas curvas de calibración, cuya obtención se
discutirá más adelante.
167
puede obtener luego el caudal. El lugar elegido para hacer el aforo o medición
debe cumplir los siguientes requisitos:
- La sección transversal debe estar bien definida; en lo posible, no debe
presentarse agradación o degradación del lecho.
168
FIGURA 7.1 Sección transversal para el método área velocidad
169
aplicación del método área velocidad o medidas con estructuras hidráulicas,
como en corrientes muy anchas o ríos torrenciales. Se puede implementar de
dos maneras así:
170
V1C1
Q= t2
(7.3)
∫ (C
t1
2 − CO )
Se inyecta un trazador en una sección dada, a un caudal constante qo, con una
concentración de trazador Co, así:
qCo
C2
1
2
Q
QC1 + qCo = (Q + q )C 2
despejando el caudal Q :
q(C 2 − C o )
Q= (7.4)
( C1 − C 2 )
171
Es importante anotar que para aplicar este método, se supone que el flujo es
permanente.
3) Materiales radioactivos: los más usados son el yodo 132, bromo 82, sodio.
Ejemplo 7.1
Una solución de sal común con una concentración de 200g/l fue descargada en
un río con un caudal constante de 25 l/s. El río tenía inicialmente una
concentración de sal de 10 ppm. Aguas abajo, se midió una concentración de
45 ppm. Cuál es el caudal en el río?
172
Solución:
Q = 113.6 m3/s
Q = CHn (7.5)
173
7.2.2.2 Método área-pendiente.
Y1
Y2
Z1
l
Z2
NR
V12 V2
h1 + = h 2 + 2 + hf (7.6)
2g 2g
donde: h= Y+Z y hf son las pérdidas de energía que se pueden hallar usando la
fórmula de Manning:
1 2 / 3 1/ 2
VA = Q = R H Sf (7.7)
n
174
donde:
V: velocidad en m/s
RH: radio hidráulico en m
Sf: pendiente de la línea de energía
A: área de la sección transversal en m2
n: coeficiente de rugosidad de Manning
hf
h f = h 1 − h 2⇒ S f =
L
2) Si en la fórmula de Manning :
1 2/ 3
K= RH A
n
el caudal puede expresarse como:
Q = KS1f / 2 (7.8)
Se encuentra un valor promedio de K para las dos secciones, el cual puede
hallarse con la media geométrica, así:
K = K 1K 2 (7.9)
175
continúa hasta que dos valores sucesivos de las pérdidas hidráulicas difieran
en muy poco.
• Cano, 1988 ⎡ ⎛ 4R ⎞ ⎤
= a⎢ − log(1.352a ) + log⎜
1
⎟⎥ (7.15)
f ⎣ ⎝ k ⎠⎦
a = 5.7798 k 0.139633 (7.15a)
1
R H6
n= (7.15b)
8⋅g
f
176
material del lecho, hallados a partir de su curva granulométrica; R es el
radio hidráulico y f es el factor de fricción de la ecuación de Darcy -
Weisbach.
- Reposo, k = 0.54 D50, para cascajos, piedras y rocas con diámetro medio
mayor de 0.03 m; para tamaños menores, el coeficiente aumenta de 0.54
a 1.0.
- Movimiento, k = 0.56 D50, para tamaños medios del sedimento mayores
de 0.03 m; el coeficiente aumenta de 0.56 a 0.78 para tamaños menores
de 0.03 m.
3. En cada intercepto, se mide la cara expuesta más larga del grano que allí
se encuentre.
177
Adicionalmente, se debe tomar una muestra de finos del fondo del cauce
para realizar la curva granulométrica completa.
n = 0.0487 ⋅ D 50
1
6 (7.17)
donde:
Ejemplo 7.2.
Solución:
Y1=3 m Y2=2.9 m
A1=30 m2 A2=29 m2
178
P1=16 P2=15.8 m2
RH1=1.875 RH2=1.875
1 1
K1 = × 30 × (1.875)1 / 2 K1 = × 29 × (1.835)1 / 2
0.025 0.025
= 1824.7 = 1738.9
en donde :
So: pendiente del canal.
Se calcula el K, así:
K = K 1K 2 = 1781.3
179
7.3 RELACIONES NIVEL-CAUDAL
El objetivo de aforar una corriente, durante varias épocas del año en una
sección determinada, es determinar lo que se conoce como curva de
calibración de la sección. Esta permite transformar niveles de agua, leídos
con una mira, en caudales. Las curvas se construyen a partir de los aforos
hechos durante un período largo de tiempo, de tal manera que se tengan
niveles bajos y altos del río. La curva tiene la forma mostrada en la figura 7.4.
1)Curva de remanso. Las curvas de remanso son perfiles del tipo M1, que se
presentan debido a la existencia de una sección de control, por ejemplo, una
presa o un vertedero. Si hay curvas de remanso, la misma altura de mira, H,
puede corresponder a dos caudales diferentes. Hay estaciones con muchos
años de registro, que son influidas por la "cola" de los embalses formando
remansos que afectan los registros de la estación de aforo. Para no perder la
serie, este problema se puede resolver instalando otra estación auxiliar, aguas
abajo y se sigue el siguiente procedimiento. Se toman lecturas de los niveles
en las dos miras, y F es la diferencia entre niveles, ver figura 7.7.
180
29
Alturas de mira en m
28
27
26
25
24
23
0 5 10 15 20 25 30
Caudal en m3/s
181
14
12
Caudal m3/s
10
8
6
4
2
0
0 0.2 0.4 0.6 0.8 1
Altura de mira m
V1 V2
+ y1 + z1 = 2 + y 2 + z 2 + hf
2g 2g
V: velocidad
y: profundidad del flujo.
z: cabeza de posición
hf: pérdidas de energía, que pueden estimarse aplicando la ecuación de
Manning, como ya se explicó.
Si:
Z=z+y
Y despreciando cabezas de velocidad se obtiene:
182
Z1 − Z 2 − h f = Fn (7.17)
yn
M1
y1
y2
1
Estación Mira auxiliar
2
Q F
= ( )m (7.18)
Qn Fn
Donde Qn es el caudal normal para una altura de la mira H dada y m es un
exponente con un valor cercano a 0,5
183
Nivel
bajando
Flujo
H
permanente
Nivel subiendo
QM 1 dh
= 1+ (7.19)
Qn Vw S o dt
Donde:
So : pendiente del canal
dh/dt: tasa de cambio del nivel del agua con el tiempo
Vw: velocidad de la onda de creciente; se asume que:
Vw = 1.4V (7.20)
184
V= velocidad halada con la ecuación de Manning
Q = C( H − H 0 )n (7.21)
donde:
Q: caudal
H: nivel medido en la mira
H0 : nivel cuando Q es cero
C y n : constantes.
185
b) Se asumen diferentes valores de H0 y se grafican log Q vs log(H-H0)
d) Se encuentran C y n
2/3
A RH
186
b) De la gráfica anterior, para un nivel máximo observado, H, se obtiene
A R H2 / 3
Caudal medio mensual, Qm. Se calcula hallando para cada mes la media
aritmética de los caudales promedios diarios.
187
aforo tenga limnígrafo. Si no es así, se habla de caudal maximo promedio
anual el cual es menor que el máximo instantáneo anual..
Los datos de caudal medio anual, mensual o diario se pueden usar para
construir la curva. Los caudales se disponen en orden descendente, usando
intervalos de clase, si el número de valores es muy grande. Si N es el número
188
de datos, la probabilidad de excedencia , P, de cualquier descarga( o valor de
clase), Q, es:
m
P= × 100 (7.23)
N
% Tiempo
FIGURA 7.11 Curva de duración
189
3) Cuando se dibuja en papel logarítmico la curva de duración, se
obtiene una línea recta, al menos en la región central. De esta
propiedad, se hallan varios coeficientes que expresan la variabilidad
del flujo en el río y pueden usarse para describir y comparar varias
corrientes.
Flujo Natural
Con Embalse
% Tiempo
190
Por medio de esta curva, se definen los siguientes caudales característicos:
- Caudal característico máximo: Caudal rebasado 10 días al año.
- Caudal característico de sequía: Caudal rebasado 355 días al año.
- Caudal de aguas bajas: caudal excedido 275 días al año o el 75 % del
tiempo.
- Caudal medio anual: es la altura de un rectángulo de área equivalente al
área bajo la curva de duración.
Existen muchos ríos del país que no tienen registros de caudal; es imposible,
entonces, obtener la curva de duración. Sin embargo, si se construye una
curva de duración regional que represente el comportamiento de una zona
hidrológicamente homogénea, se pueden hallar caudales de diseño en
regiones donde se tenga poca o ninguna información.
Q
Z= (7.24)
Q
medio
Donde:
Z: Caudal adimensional
Q: Caudal registrado
Qmedio: Caudal promedio diario multianual
191
estaciones que presentaron un comportamiento más uniforme, figura 7.13
10
8
Q/Qmedio
0
0 10 20 30 40 50 60 70 80 90 100
p(%)
192
Solución:
140-120.1 0 1 5 6 6 0.55
120-100.1 2 7 10 19 25 2.28
100-80.1 12 18 15 45 70 6.38
193
140
120
Q m3/s 100
80
60
40
20
0
0.55 6.38 21.51 60.62 87.78 99.45
P%
7.6.CURVA DE MASAS
t
V = ∫ Qdt (7.25)
to
194
FIGURA 6.15 Curva de masas
Ejemplo 7.4
195
Volumen acumulado
Mes Días acumulados Q m3/s
Mm3 x 103
15 110 142.5
Enero
31 95 273.5
45 85 376.4
Febrero
59 71 462.2
74 63 543
Marzo
90 52 615
105 41 668.8
Abril
120 31 709.0
135 20 734.9
Mayo
151 18 759.7
166 20 785.68
Junio
181 42 840.08
196 125 1002.08
Julio
212 270 1375.08
227 410 1907.08
Agosto
243 460 2543.08
258 405 3068.08
Septiembre
273 250 3092.08
288 140 3573.58
Octubre
304 96 3707.38
319 63 3788.88
Noviembre
334 55 3860.08
349 56 3932.68
Diciembre
365 100 4070.68
Solución:
196
4500
Volumen acumulado en Mm x 10
3
4000
3500
Volumen
3000
de embalse
2500
2000
1500
1000
500
0 105
135
166
196
227
258
288
319
349
15
45
74
Días acumulados
197
Capítulo 8
LA INFILTRACIÓN
198
8.1 INTRODUCCION.
Los principales factores que afectan la infiltración en una cuenca son: textura,
estructura y condición del suelo, vegetación y características de la lluvia.
199
Vv
n= × 100 (8.1)
Vt
Donde:
Vv: volumen de vacíos
VT: volumen total de la muestra.
200
V
⎛ ⎞
V v = P2 - P1 ⇒ n = ⎜ V v ⎟%
⎝V⎠
b) Medidas in situ. Aunque hay varios métodos para realizar estas medidas:
diagrama neutrón, bombardeo del terreno con rayos gamma, velocidad
del sonido, el más común es inferir porosidades a través de medidas de
resistividad del terreno, como se explica a continuación. Con excepción
de las arcillas, los minerales usualmente presentes en el terreno son
aislantes y la electricidad circula en el terreno por medio de la fase
líquida. La resistividad es entonces, una función de la porosidad. Se
define una factor de formación F así:
Resistividad de la roca
F= (8.2)
Resistividad del aguacontenidaenlaroca
201
Si se usa el factor F, la geofísica propone la expresión de Archie para
encontrar la porosidad total
C
F= (8.3)
nm
202
• Parámetros granulométricos. Hay una serie de parámetros que además de
permitir comparar unas curvas con otras permiten igualmente
caracterizar la posición de la curva y el grado de homogeneidad de la
granulometría. Dichos parámetros son los siguientes:
• El diámetro eficaz d10, que es aquél para el cual el 10% del peso de la
muestra tiene un diámetro inferior.
203
favorecen más la infiltración que suelos con granos angulosos.
Por esta razón, para la construcción de filtros es norma utilizar
material de río, en vez del material de cantera, mucho más
anguloso y por tanto menos permeable.
204
a una presión de 15 atmósferas y hallando después su grado
de humedad. (SAI, 1987; UNESCO, 1981).
8.2.3 La vegetación.
205
FIGURA 8.3 Perfil del suelo típico del Oriente antioqueño
Este proceso fué estudiado por primera vez por Horton, en 1950. El propuso la
206
siguiente relación para determinar la capacidad de infiltración:
f = f c + (f 0 − f c )e − Kt (8.4)
donde:
f: capacidad de infiltración en un tiempo t en mm/h.
fc: capacidad de equilibrio, que se da cuando el suelo está completamente
saturado en mm/h.
f0: capacidad inicial en mm/h.
t: tiempo en horas
K: constante que representa la tasa de decrecimiento de esa capacidad.
∂θ ∂ ⎡ ∂ω(θ ) ⎤ ∂K (θ)
= − ⎢K (θ) −
∂z ⎥⎦
(8.5)
∂t ∂z ⎣ ∂z
207
FIGURA 8.4 Infiltración, según Horton
donde:
θ: contenido de humedad en porcentaje.
Ψ(θ): presión de succión capilar en cm.
K(θ): conductividad hidráulica no saturada en cm/s.
z: distancia bajo la superficie del terreno.
Resolver la anterior ecuación es difícil, por las dificultades que se presentan
para la determinación de la conductividad hidráulica no saturada y la presión
de succión capilar. Philip (1957) resolvió la ecuación de Richards para el caso
donde la intensidad de la precipitación es mayor que que la capacidad de
infiltración del terreno. La ecuación de Philip tiene la forma de:
f = (0.5)At −1 / 2 + B
(8.6)
F = At 1 / 2 + Bt
donde:
f: capacidad de infiltración en mm/h.
F: lámina infiltrada en mm.
208
A, B: constantes que dependen del tipo de terreno.
ω ∆θ
f (t ) = K ( + 1) (8.7)
F( t )
donde:
∆υ=n-υi.
υi: humedad en el tiempo i.
K: conductividad hidráulica que varía con el tiempo.
F: Lámina infiltrada en un tiempo t.
209
metodologías generales para estimarla en una cuenca: 1) haciendo medidas
directas en el campo por medio de infiltrómetros, 2) utilizando índices
obtenidos de los hidrogramas y usando métodos empíricos, tal como el
propuesto por el Soil Conservation Service.
8.4.1 Infiltrómetros.
I = P − ESD (8.8)
210
Ejemplo 8.1
t min 0 2 5 10 20 30 60 90 150
Vol
0 300 650 1190 1950 2500 3350 3900 4600
cm3
Solución:
A = Π D 2 / 4 = 962cm 2
211
y la capacidad de infiltración se halla como:
V
f=
A∆ t
∆t min 2 3 5 10 10 30 30 60
f cm/min 0.16 0.12 0.11 0.08 0.06 0.03 0.02 0.02
0.2
0.15
f cm/min
0.1
0.05
0
2 3 5 10 10 30 30 60
Tiem po m in
212
el índice φ. El índice φ se define como la intensidad media por encima de la
cual todo lo precipitado se transforma en escorrentía superficial directa. El
principio en que se basa el método es la separación en el hidrograma de la
parte que corresponde a la escorrentía superficial directa del llamado flujo
base. Esta parte, expresada en lámina, debe ser igual al área por encima del
índice en el hietograma. Gráficamente puede verse en la Figura 8.6.
Para separar el flujo base de la escorrentía superficial directa existen varios
métodos:
N = 0.83 A2 (8.9)
213
A es el área en Km2.
214
FIGURA 8.10 Separación del flujo base por método c.
Ejemplo 8.2
T(h) 1 2 3 4 5 6 7 8
i cm/h 0.4 0.9 1.5 2.3 1.8 1.6 1 0.5
Solución:
215
FIGURA 8.11 Hietograma de la tormenta, ejemplo 8.2.
I = 10 - 5.8 = 4.2 cm
y el índice φ será:
y el índice φ será:
T (h) 1 2 3 4 5 6 7 8
ESD 0 0.35 0.95 1.75 1.25 1.05 0.45 0
(cm)
216
para determinar las tormentas de diseño, cuando se estudian caudales
máximos. Sin embargo, sus resultados deben mirarse con mucho cuidado, pues
los coeficientes de este método fueron desarrollados para zonas de Estados
Unidos muy diferentes a las tropicales de Colombia.
2
(P - Ia )
Pe = (8.10)
P - Ia + S
217
⎛ 1000 ⎞
S=⎜ - 10 ⎟ (8.11)
⎝ CN ⎠
Ia = 0.2 S (8.12)
218
90 78 96
85 70 94
80 63 91
75 57 88
70 51 85
65 45 82
60 40 78
55 35 74
50 31 70
45 26 65
40 22 60
35 18 55
30 15 50
219
Tabla 8.4. VALORES DE CN PARA AMC II.
Caminos
De terracería - 72 82 87 89
De superficie dura - 74 84 90 92
219
94
Capítulo 9
LA EVAPORACIÓN
220
9.1 GENERALIDADES
Uno de los componentes más importantes del ciclo hidrológico, y tal vez el
más díficil de cuantificar es la evapotranspiración, que es el operador natural
mediante el cual parte de la precipitación que cae sobre la tierra regresa de
nuevo al sistema atmosférico, en forma de vapor de agua. Aproximadamente,
el 75% de la precipitación promedio anual retorna a la atmósfera por medio de
la evapotranspiración (Mutreja 1986). En climas tropicales, como el
colombiano, las pérdidas por evaporación en embalses, ríos y canales son
significativas respecto al agua total disponible. El agua que regresa a la
atmósfera en forma de vapor puede provenir de varios lugares como:
- Suelo húmedo.
221
La evapotranspiración se define, pues, como la suma de la evaporación desde
la tierra (superficies de agua, suelo o vegetal) y la transpiración de los
vegetales y seres vivientes.
- Radiación solar
- Temperatura del aire y del agua
- Presión
- Humedad relativa
- Viento
- Calidad de agua de la fuente de evaporación
- Geometría de la superficie del agua
222
suministro de agua potable e irrigación, etc. Mucho mayor será su importancia,
cuando se hace el análisis en zonas áridas y/o de alta evaporación.
La evaporación desde los suelos tenderá a ser más significativa mientras más
frecuentes sean las aplicaciones del agua, ya sea por medio de la irrigación o
de la precipitación.
223
9.2 FACTORES QUE AFECTAN LA EVAPORACION
E = C(e w − e a ) (9.1)
donde :
E: evaporación en mm/d
ew y ea: presiones de vapor en mm de mercurio
C. constante.
Esta ecuación fue deducida por Dalton (1820). Según ella, la evaporación
continúa hasta que ew=ea, cuando ew∃ea, se produce la condensación del vapor
de agua.
9.2.2 La temperatura.
9.2.3 El viento.
224
agua. Para grandes cuerpos de agua, se necesitan velocidades del viento muy
altas para crear tasas máximas de evaporación (Subramaya,1984).
9.3 LA TRANSPIRACION
225
de las superficies de las raíces, tipo de hojas, densidad de las plantas y
conductividad hidráulica del suelo, etc.
226
9.4 MEDIDA DE LA EVAPORACION SOBRE UNA SUPERFICIE DE
AGUA
A) Balance hídrico
B) Balance energético
C) Técnica de transferencia de masas
D) Evaporímetro o tanque evaporador
227
E = S + (Q1 + Q 2 ) + P − ESD − I (9.2)
donde:
S: cambio en el agua almacenada en el embalse en el intervalo de tiempo t
Desde el punto de vista teórico, este método es muy bueno, pero desde el
punto de vista práctico, es bastante inexacto y tendría una aplicación nula
para intervalos de tiempo pequeños: por ejemplo, del orden de una semana o
menos para el caso de un embalse, pues la magnitud del error sería grande.
Los errores en el cálculo de la evaporación se deben en esencia a los errores
en los cálculos de las variables independientes, principalmente en el de la
infiltración. Cuando este método se usa y el valor de la infiltración es del
mismo orden o mayor que el de la evaporación, se introducen grandes
errores provenientes de la dificultad de estimar la infiltración; da buenos
resultados para cuencas con características muy definidas y homogéneas.
Durante épocas de sequía y para infiltraciones pequeñas respecto a la
evaporación, da excelentes resultados, pues la precipitación es nula y los
caudales de entrada y salida son muy bajos o nulos. Hay que tener en cuenta
que la evaporación y su estimación en tiempo de sequías, no necesariamente
es una condición representativa de todas las épocas.Bajo condiciones
óptimas, la exactitud puede ser del orden del 5 al 10 por ciento.
228
En pocas palabras, es la evaluación de la energía almacenada en el sistema, en
función de la energía que entra y sale de él. La ecuación se puede plantear:
Intercambio de
calor
229
9.4.3 Técnicas de transferencia de masa
E = kf (u )(e w − e a ) (9.3)
donde:
f(u): función de la velocidad, u, del viento
k: constante
Existen en la literatura muchas expresiones empíricas de este tipo, como las de
Meyer, Penman y otros. La fórmula de Meyer tiene la forma:
u
E = k (1 + )(e w − e a ) (9.4)
10
donde:
E: evaporación en pulgadas por día.
u: velocidad del viento en mph, a 25 pies por encima de la superficie del lago.
k: coeficiente que vale 0.36 para lagos ordinarios y 0.50 para lagos poco
profundos.
Existen tablas para hallar los valores de ew y ea, según la temperatura, como la
tabla 9.1.
230
TABLA 9.1 Presión de vapor a varias temperaturas
Presión de vapor
Presión de vapor
TºF Pulagadas de
milibares
mercurio
32 6.11 0.18
40 8.36 0.25
50 12.19 0.36
60 17.51 0.56
70 24.79 0.74
80 43.61 1.03
90 47.68 1.42
100 64.88 1.94
Ejemplo 9.1
Solución:
ew = 0.58 pulgadas de Hg
ea = 1.29 x 0.20 = 0.26 pulgadas de Hg
Aplicando la expresión de Meyer, se obtiene:
E = 0.36 (0.58 - 0.26)(1+ 10/10) = 0.23 pul/día
231
resultados aceptables para intervalos de tiempo no muy pequeños, con
coeficientes relativamente estables. El tipo de tanque más común es el llamado
Clase A, construído de lámina galvanizada, sin pintar de unos 121 cms de
diámetro y de aproximadamente 25,4 cms de profundidad, apoyado sobre unas
vigas de madera a una altura sobre la superficie del suelo de unos 25,5 cms. La
regla de operación del evaporímetro clase A es simplemente mantenerlo lleno
de agua hasta una profundidad máxima de 20 cms y una mínima de 17,5 cms,
llevando a cabo al menos una lectura diaria, corregida por la precipitación
acumulada durante el período de tiempo entre observaciones.
Er = KET (9.5)
donde:
Er: evaporación real
ET: evaporación en el tanque
K: constante de proporcionalidad, que generalmente varía entre 0,60 y 0,85,
con un valor promedio general de 0,70
232
informacion adquirida para intervalos estacionales deben ser analizada con
cuidado.
Ejemplo 9.2
Semana 1 2 3 4
Precipitación 0 26.42 46.74 10.07
mm
Agua añadida 22.86 1.02 -17.8 23.6
mm
Solución:
Semana 1 2 3 4
ET mm 18.29 21.95 23.15 27.44
ET=90.83 mm
Er=0.8 (0.09083x607031)=44109.3 m3
233
9.5 EVAPOTRANSPIRACION POTENCIAL
9.6.1 Lisímetros.
234
Tiene un sistema para medir el agua drenada, el agua escurrida; se debe tener
además, un pluviométro para medir la precipitación. Se puede entonces
plantear la ecuación de balance hídrico:
donde:
P: precipitación, medida con un pluviómetro vecino
I: agua infiltrada medida con el lisímetro
ESD: escorrentía superficial directa, medida en la superficie del terreno
∆S: cambio en el almacenamiento de agua al interior del terreno, medido
usualmente con una sonda de neutrones. Se deben tomar medidas en períodos
mayores de 10 días, para que el ∆S sea significativo. En general, se usan
períodos de un mes o mayores.
235
lisímetros, presenta la ventaja de que se mide la evaporación en las
condiciones naturales del terreno. Sin embargo, los resultados pueden verse
alterados por los movimientos del agua subterránea hacia y desde la parcela,
que muchas veces son imposibles de cuantificar con certeza.
AH n + Ea γ
EVP = (9.7)
A+γ
Donde:
EVP: evapotranspiración diaria en mm
γ: constante sicrométrica=0.49 mm Hg/ºC
A: pendiente de la curva: presión de vapor vs temperatura, a la temperatura del
aire, según la siguiente ecuación:
12.27t
e w = 4.584 exp( (9.8)
237.3 + t )
t: temperatura en ºC
ew:presión de vapor en mm de mercurio.. Se da en mm Hg/ C.
236
La tabla 9.2 da los valores de A para algunas temperaturas usuales.
n
H n = H a (1 + r )(a + b Nn ) − σTa ( 0.56 − 0.092 e a )(0.10 + 0.90
4
) (9.9)
N
Donde:
Arbustos 0.15-0.25
Tierra 0.05-0.45
Nieve 0.45-0.90
Agua 0.05
237
TABLA 9.3 Valores promedios mensuales de Ha, en mm de agua
evaporable por día
Latitud E F M A M J J A S O N D
N
0 14.5 15 15.2 14.7 13.9 13.4 13.5 14.2 14.9 15.0 14.6 14.3
10º 12.8 13.9 14.8 15.2 15 14.8 14.8 15.0 14.9 14.1 13.1 12.4
20º 10.8 12.3 13.9 15.2 15.7 15.8 15.7 15.3 14.4 12.9 11.2 10.3
30º 8.5 10.5 12.7 14.8 16 16.5 16.2 15.3 13.5 11.3 9.1 7.9
40º 6.0 8.3 11 13.9 15.9 16.7 16.3 14.8 12.2 9.3 6.7 5.4
50º 3.6 5.9 9.1 12.7 15.4 16.7 16.1 13.9 10.5 7.1 4.3 3.0
a = 0.29 cos φ
u2 (9.10)
E a = 0.35(1 + )(e w − e a )
160
238
P
EVPR = para (P/L)∃0.316
P2 (9.11)
0.9 + 2
L
Si (P/L)′0.316 , EVP=P
donde:
EVPR: evaporación real en mm/año
P: precipitación en mm, del año considerado.
L: se expresa como:
Latitud E F M A M J J A S O N D
N
0º 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1 12.1
10º 11.6 11.8 12.1 12.4 12.6 12.7 12.6 12.4 12.9 11.9 11.7 11.5
20º 11.1 11-.5 12.0 12.6 13.1 13.3 13.2 12.8 12.3 11.7 11.2 10.9
30º 10.4 11.1 12.0 12.9 13.7 14.1 13.9 13.2 12.4 11.5 10.6 10.2
40º 9.6 10.7 11.9 13.2 14.4 15.0 14.7 13.8 12.5 11.2 10.0 9.4
50º 8.6 10.1 11.8 13.8 15.4 16.4 16.0 14.5 12.7 10.8 9.1 8.1
10T&& a
EVP = 1.6L a ( ) (9.13)
IT
239
Donde:
EVP: evapotranspiración potencial mensual en cm.
La: ajuste para el número de horas de brillo solar y días en el mes, relacionado
con la latitud. Se da en la siguiente tabla.
Lat. E F M A M J J A S O N D
N
0º 1.04 0.94 1.04 1.01 1.04 1.01 1.04 1.04 1.01 1.04 1.01 1.04
10º 1.0 0.91 1.03 1.03 1.08 1.06 1.08 1.07 1.02 1.02 0.98 0.99
15º 0.97 0.91 1.03 1.04 1.11 1.08 1.12 1.08 1.02 1.01 0.95 0.97
20º 0.95 0.90 1.03 1.05 1.13 1.11 1.14 1.11 1.02 1.0 0.93 0.94
25º 0.93 0.89 1.03 1.06 1.15 1.14 1.17 1.12 1.02 0.99 0.91 0.91
30º 0.90 0.87 1.03 1.08 1.18 1.17 1.20 1.14 1.03 0.98 0.89 0.88
40º 0.84 0.83 1.03 1.11 1.24 1.25 1.27 1.18 1.04 0.96 0.83 0.81
12
i = ∑i
1 (9.14)
&& / 5)1.514
i = (T
a: constante empírica:
240
EVP = 2.54KF
(9.16)
F = ∑ p h Tf / 100
donde:
EVP: evapotranspiración potencial e en cm, en la estación de crecimiento.
K: coeficiente que depende del tipo de cultivo,tabla 9.6.
ph:porcentaje de horas de sol al día. Es función de la latitud, tabla 9.7.
T&&f : temperatura promedio mensual en ºF
Donde:
EVPr: evapotranspiración real en mm/d
h: cota sobre el nivel del mar, en m
241
HR: humedad relativa media diurna, hallada en función de la humedad a las 8 y
a las 14 horas.
EVP = 0.473H aC
(9.20)
C = C T C H Cu Cs C E C M
242
Donde:
TC T
CT = 0.393 + 0.5592( ) + 0.0476( C ) 2
20 20
U U 2
Cu = 0.708 + 0.3276( ) − 0.036( ) (9.21)
96.56 96.56
H H H
C H = 1.25 + 0.348( m ) + 0.12( m ) 2 − 0.0218( m )4
40 40 40
Lat. E F M A M J J A S O N D
N
0º 8.5 7.66 8.49 8.21 8.50 8.22 8.50 8.49 8.21 8.50 8.22 8.50
10º 8.13 7.47 8.45 8.37 8.81 8.60 8.86 8.71 8.25 8.34 7.91 8.10
15º 7.94 7.36 8.43 8.44 8.98 8.8 9.05 8.83 8.28 8.26 7.75 7.88
20º 7.74 7.25 8.41 8.52 9.15 9.0 9.25 8.96 8.30 8.18 7.58 7.66
25º 7.53 7.14 8.39 8.61 9.33 9.23 9.45 9.09 8.32 8.09 7.40 7.42
30º 7.30 7.03 8.38 8.72 9.53 9.49 9.67 9.22 8.33 7.99 7.19 7.15
35º 7.05 6.88 8.35 8.83 9.76 9.77 9.93 9.37 8.36 7.87 6.97 6.86
40º 6.76 6.72 8.33 8.95 10.02 10.08 10.22 9.54 8.39 7.75 6.72 6.52
HR H
C H = 1.25 − 0.212( ) − 0.038( R )5
57.4 57.4
(9.22)
CS = 0.542 + 0.64(S / 80) − 0.4992(S / 80) 2 + 0.3174(S / 80) 3
CE = 0.97 + 0.03(h / 304.88)
Donde:
TC: temperatura media diaria en ºC, la "opción por defecto" es 20ºC.
U: velocidad promedio diaria del viento en Km/d, la "opción por defecto" es
96.56 Km/d.
243
Hm: humedad relativa promedio medida a mediodía, la "opción por defecto" es
40%.
HR: humedad relativa promedio diaria, obtenida con los registros de las 11 y
las 17 horas; la "opción por defecto" es 57.4%.
S: porcentaje de brillo solar, S=n/N. N se obtiene con la tabla 9.4. La "opción
por defecto" es 80%.
h: elevación de la estación. La opción por defecto es 304.88
CM: coeficiente de uso del agua por la vegetación, función de la latitud y de la
época del año. En climas tropicales, puede tener valores promedios cercanos a
0.9.
Los resultados obtenidos con todas estas expresiones difieren mucho, pues,
generalmente, en las ecuaciones empíricas, los coeficientes son calibrados para
un determinado lugar. Los resultados de la evapotranspiración calculada con
datos climáticos de la Estación la Selva, en Rionegro (Antioquia), tabla 9.8,
muestran esta situación.
244
Método 4 Hargreaves.
Método 5: Thornthwaite
Método 6: Evaporación medida con tanque.
Método 7: evaporación de tanque utilizando un coeficiente de 0.8.
- Las regiones con los valores más bajos, por debajo de los 900
milímetros al año, se encuentran localizadas en las partes altas de las
245
cordilleras, en niveles superiores a los 3.000 metros sobre el nivel del
mar.
donde:
T: temperatura promedio diaria anual en ºC
h:altura sobre el nivel del mar en m.
246
Capítulo 10
RELACIONES PRECIPITACIÓN- ESCORRENTÍA
247
10.1 INTRODUCCION
248
10.2 ANALISIS DE HIDROGRAMAS
249
La forma de la rama ascendente del hidrograma (segmento BD, en la figura
10.1) depende de la duración, de la intensidad y distribución espacial de la
lluvia y de las condiciones antecedentes del suelo.
Qt = Qo e− K t (10.1)
Esta curva señala la disminución gradual del caudal de una corriente de agua,
cuando ha cesado la escorrentía total. En muchas hoyas hidrográficas, esta
curva representa una característica bien definida del aporte de aguas
subterráneas al caudal total, pero existen algunos casos en los cuales no se
puede encontrar una ley general que permita expresar la curva de agotamiento
en forma matemática y es necesario hacer una serie de aproximaciones que
simplifiquen el problema, sacrificando en gran parte la precisión de los
resultados.
250
Hay varios procedimientos para estimar la curva de recesión de una cuenca.
Básicamente, se busca hallar el coeficiente k. Uno de esos procedimientos es
la llamada curva maestra.. El procedimiento para hallarla es el siguiente:
La figura 10.2, tomada de Mccuen 1989, muestra los resultados obtenidos con
este procedimiento.
Otros métodos para obtener la curva de recesión, junto con sus ventajas y
desventajas pueden consultarse en Blandon y Saldarriaga 1997.
251
La curva de recesión se usa para estimar caudales mínimos asociados a un
período de retorno determinado. Si se conoce el coeficiente de recesión k en
función de parámetros del acuífero y parámetros morfométricos de la cuenca
y se conocen las duraciones máximas de períodos sin lluvia para un período
de retorno determinado, es posible estimar, con la ecuación 10.1, caudales
mínimos asociados a diferentes períodos de retorno. En Antioquia y
Risaralda, se hicieron estudios para determinar caudales mínimos con la curva
de recesión. (Vélez, Blandon, Saldarriaga, 1998; Vélez , Smith, 1998,
Universidad Nacional, 1998). Se obtuvieron relaciones del tipo: k=f(A, d, e,
L) para el coeficiente de recesión k., donde A es el área de la cuenca, d, es la
densidad de drenaje, e es el espesor del acuífero y L es la longitud de la
corriente principal.
252
FIGURA 10.3 Parametros de la hidrógrafa de escorrentía superficial
Donde:
253
interés. El tiempo de concentración (Tc) es uno de los parámetros más
importantes en los modelos precipitación - escorrentía, pues la duración de
la tormenta de diseño se define con base en él. La duración crítica de la
lluvia debe asumirse como igual al tiempo de concentración, pues para
duraciones menores que Tc, no toda el área de la cuenca contribuye; y para
duraciones más grandes que Tc, no hay incremento en el área contribuyente;
en cambio la intensidad de la lluvia de una frecuencia dada disminuye. Se
asume que para duraciones menores que el tiempo de concentración, el
efecto de la reducción en el área contribuyente es mayor que el del
incremento en la intensidad de la lluvia. (Smith, Vélez,1997). Es clara,
entonces, la dependencia de este parámetro de variables morfométricas,
tales como el área, pendiente de la cuenca, longitud del cauce principal, etc.
Tc = Tb − D (10.2)
254
Se ha definido el tiempo de concentración, con base en el tiempo de rezago
TR , así:
Tc = Tb − Tr (10.3)
Tc = Tb − Tp (10.4)
• Témez (1978) ⎛ L ⎞
0.75
L A0.4
• Williams (1922) Tc = (10.6)
DSo 0.2
255
D : diámetro de una cuenca circular, con área A, en millas.
• Kirpich (1990) ⎛ L ⎞
0.77
Tc = 0.066 ⎜ ⎟ (10.7)
⎝ So ⎠
Tc=5⎜ ⎟ (10.8)
⎝ So ⎠
256
A : área de la cuenca, en kilómetros cuadrados.
L : longitud del cauce principal, en kilómetros.
So : diferencia de cotas entre puntos extremos de la corriente sobre L,
en m/m.
Lc 3
K =
H
Lc : distancia desde el sitio de interés al punto en el cual la corriente
principal corta la divisoria, en kilómetros.
H : diferencia de cotas entre puntos extremos de la corriente, en pies.
0.35
⎛ LL ⎞
• Linsley Tc = Ct ⎜⎜ 0.5 ⎟⎟ (10.12)
⎝ So ⎠
• Snyder Tc = C t (L L )
0.3 (10.13)
257
FL
• Bransby - Williams Tc= (10.14)
A S o 0 .2
0.1
L
• Pérez Monteagudo (1985) Tc = (10.15)
VR
VR se expresa como:
0 .6
⎛ AE ⎞
VR = 72⎜ ⎟ (10.16)
⎝ L ⎠
VR : Velocidad, en km/h.
AE : diferencia de cotas, en kilómetros, en el cauce principal, desde el
punto más alto al más bajo.
La tabla 10.1 (Smith, Vélez, 1997) muestra para algunas cuencas del
departamento de Antioquia. Los resultados tan diferentes que producen
estas expresiones.
258
TABLA 10.1 Tiempos de concentración (h), calculados por diferentes
métodos, en algunas cuencas de Antioquia (Smith, Vélez, 1997)
259
Campo y Múnera (1997) hicieron, también en Antioquia, un análisis del
tiempo de concentración, a partir de más de 2100 registros simultáneos de
pluviógrafos y limnígrafos, y encontraron que las expresiones que más se
ajustaban a los tiempos de concentración reales (hallados con los registros)
fueron las de Témez, Kirpich, Giandotti y Pérez. En el mismo trabajo, se
calibraron para diferentes zonas, ecuaciones para el tiempo de
concentración, en función de parámetros morfométricos de las cuencas,
como la siguiente, aplicable para la zona de los embalses del oriente
antioqueño:
donde:
TC: tiempo de concentración, en horas.
A: área de la cuenca, en Km2.
L: longitud del cauce principal, en km.
S: pendiente promedio de la cuenca, en %
260
En 1932, Sherman desarrolló un gráfico unitario, asumiendo que para una
duración dada de lluvia efectiva, el tiempo base de la hidrógrafa debería
permanecer constante. Las hipótesis en que está basada la teoría de la
hidrógrafa unitaria son las siguientes:
261
FIGURA 10.4 Hidrógrafas unitarias proporcionales
En el caso de que no se tengan lluvias aisladas, sino una lluvia larga con
intensidades diferentes, se divide esa lluvia en duraciones iguales, con una
intensidad constante para cada duración. Se supone que la precipitación en
intervalos iguales, produce hidrogramas proporcionales de escorrentía
superficial directa. La suma de esos hidrogramas dará el hidrograma total,
figura10.5
262
Q1 = P1 Y1
Q 2 = P1 Y2 + P2 Y1
Q 3 = P3 Y3 + P2 Y2 + P3 Y1
. (10.18)
.
.
Q n −1 = Pm Yn − m + Pm −1Yn −m +1
m
Q n = ∑ Pm Yn −m +1 (10.19)
1
Q = PY (10.20)
en donde Q es una matriz Nx1; P es una matriz NxM, y P es una matriz Mx1.
N = L + M −1 (10.21)
263
FIGURA 10.5 Superposición de HU
264
10.3.1 Obtención de la hidrógrafa unitaria para lluvias simples.
El procedimiento es el siguiente:
Y = ( P , P ) −1 Q (10.22)
Si se denomina:
L: lámina total precipitada.
LESD: lámina de escorrentía superficial directa.
A. área de la cuenca.
VE: volumen de escorrentía superficial directa.
265
Vu: volumen correspondiente a un HU, producido por una lluvia efectiva de 1
mm
VT: volumen total precipitado
Q: ordenadas del hidrograma de escorrentía superficial directa.
Y. ordenadas del HU.
t: tiempo
Se puede plantear la siguiente ecuación:
VE L
Q=( )( )A
VT t
pero:
VT = AL , lo que implica que:
VE LA VE
Q= = (10.23)
(LA )t t
Vu L (1mmA ) L
Y= A= A
VT t (LA ) t
(10.24)
1mmA
=
t
Q
Y= (10.25)
( VE / A ) × 1mm
266
hidrograma unitario de T horas de duración, obtenido usando solo una
tormenta, puede dar resultados poco confiables por lo que se aconseja
promediar los hidrogramas unitarios de varias tormentas de la misma
duración, para obtener el hidrograma unitario único para esa duración y esa
cuenca. Este promedio no debe hacerse tomando la media aritmética de las
coordenadas concurrentes, puesto que si los máximos no ocurren al mismo
tiempo, el máximo medio puede ser inferior a algunos de los máximos
individuales. El procedimiento adecuado es calcular el promedio de los
caudales máximos y el de sus tiempos. A continuación se traza el hidrograma
unitario medio, adecuandolo a la forma de los hidrogramas unitarios
obtenidos, de tal manera, que pase por el caudal medio máximo calculado y
que tenga un volumen (en lámina de agua) igual a la unidad.
t
Q(t) = ∫ U(w)P(t - w)dw (10.26)
0
267
t
Q(t) = ∫ U(t - w)P(w)dw (10.27)
0
t
1
T t∫-T
U(T, t) = U(w)dw (10.28)
N
Q(t) = ∑ U(T, i)P(t - i + 1)Dt (10.29)
i =1
10.4 LA CURVA EN S
268
FIGURA 10.6 Curva en S
S1 = Y1
S 2 = Y2
S3 = Y3
.
(10.30)
.
.
n
S n = ∑ Yi
1
269
Ejemplo 10.1
Dada una lluvia de duración unitaria con una profundidad efectiva de 1.5
pulgadas, que produce el hidrograma mostrado en la tabla 10.2 en una cuenca
de 1.98 millas2 , se desea:
a)Determinar la hidrógrafa unitaria.
t(h) 1 2 3 4 5 6 7
Q 110 98 220 512 585 460 330
(p3/s)
t (h) 8 9 10 11 12 13
Q(p3/s) 210 150 105 75 60 54
t(h) 1 2 3 4
Solución:
El primer paso sería separar el flujo base del hidrograma. La figura 10.7
muestra el hidrograma total. El método más sencillo para separar el flujo
base, tal como se vió en capítulos anteriores, es trazar una línea recta desde el
punto donde empieza la rama descendente hasta donde termina el hidrograma
de escorrentía superficial directa. Si se escoge como punto inicial el de
coordenadas (2,98) y como punto final (11, 75), la recta punteada en la figura
será la separación entre el flujo base y la escorrentía superficial directa.
270
700
600
500
Q (p /s)
400
3
300
200
100
0
1 2 3 4 5 6 7 8 9 10 11 12 13
t (h)
271
ORDENADAS DE LA HIDRÓGRAFA UNITARIA
272
Q1 = Y1 P1
= 83.1 x 0.4 = 33.24 m3/s
Q2 = Y2 P1 + Y1 P2
= 279.4 x 0.4 + 83.1 x 1.1 = 203.2 m3/s
Q3 = Y3 P1 + Y2 P2 + Y1 P3
= 329.8 x 0.4 + 279.4 x 1.1 + 83.1 x 2 = 605.5 m3/s
Q4 = Y4P1 + Y3 p2 + Y2 P3 + Y1 P4
= 248.1 x 0.4 + 329.8 x 1.1 + 279.4 x 2 + 83.1 x 1.5
= 1145.5 m3/s
Q5 = Y5 P1 + Y4 P2 + Y3 P3 + Y2 P4
= 163.2 x 0.4 + 248.1 x 1.1 + 329.8 x 2 + 279.4 x 1.5 = 1416.9 m3/s
Q6 = Y6 P1 + Y5 P2 + Y4 P3 + Y3 P4
= 84.9 x 0.4 + 163.2 x 1.1 + 248.1 x 2 + 329.8 x 1.5
= 1204.4 m3/s
Q7 = Y7 P1 + Y6 P2 + Y5 P3 + Y4 P4
= 46.6 x 0.4 + 84.9 x 1.1 + 163.2 x 2 + 248.1 x 1.5
= 810.6 m3/s
Q8 = Y8 P1 + Y7 P2 + Y6 P3 + Y5 P4
18.3 x 0.4 + 46.6 x 1.1 + 84.9 x 2 + 163.2 x 1.5
= 473.2 m3/s
Q9 = Y8 P2 + Y7 P3 + Y6 P4
= 18.3 x 1.1 + 46.6 x 2 + 84.9 x 1.5
= 240.7 m3/s
Q10 = Y8 P3 + Y7 P4
= 18.3 x 2 +46.6 x 1.5
= 106.5 m3/s
273
Q11 = Y8 P4
= 18.3 x 1.5 = 22.5 m3/s
Ejemplo 10.2
Solución:
Qp=2.8 m/s
40 min
100 min
274
2.8 40
=
Yi ti
Yi m3/s-mm 0 0.7 1.4 2.1 2.8 2.3 1.9 1.4 0.9 0.47 0
T (min) 0 10 20 30 40 50 60 70 80 90 100
P1=47.5 mm y P2=22.5 mm
m
Q n = ∑ Pm Yn −m +1
1
se tiene:
275
entre las dos , multiplicando el resultado por 10/20. Se halla la siguiente
tabla:
Las hidrógrafas unitarias sintéticas pueden ser de tres tipos (Chow, 1994):
276
2) Aquéllas basadas en hidrógrafas unitarias adimensionales: la del Soil
Conservation Service y la de William y Hann.
• Snyder
• Servicio de Conservación de Suelos (SCS)
• Williams y Hann
277
El modelo pretende derivar un Hidrograma Unitario, a partir de algunas
características físicas de la cuenca, para ser aplicado en las cuencas donde
no se posea registro de caudal.
Snyder asume que el rezago es constante para una cuenca, ya que depende
de algunas de sus características físicas y no está determinado por el tipo de
lluvia o sus variaciones.
L0c.8 (S + 1)0.7
TL = (10.31)
1900 * S c0.5
En donde:
TL : tiempo de rezago en horas.
Lc: longitud del canal principal en pies.
S : factor de retención o almacenamiento calculado en términos del número
de curva, CN
Sc: pendiente de la cuenca en porcentaje.
TL = 0.6TC (10.32)
278
Duración de la lluvia seleccionada por Snyder. Snyder consideró lluvias
que estuvieran de acuerdo con el tamaño de la cuenca, definiendo para ellas
una duración de 1/5.5 veces el tiempo al pico de la cuenca, es decir:
TL
ts = (10.33)
5 .5
Caudal pico. El modelo propone calcular el caudal pico por milla cuadrada,
up, como:
640
u p = Cp (10.34)
TL
Donde:
up: caudal pico del Hidrograma Unitario por unidad de área, en pie3/ s.mi2 ,
Cp : coeficiente.
TL :l tiempo de rezago.
640
u p = Cp
[TL + (T − t s ) / 4]
(10.35)
279
En donde T es la duración (en horas) de la precipitación efectiva a la cual se
le va a calcular el Hidrograma Unitario Sintético.
Una vez obtenido el caudal pico por unidad de área de la cuenca, el caudal
pico total se obtiene como:
Up = up A (10.36)
En donde:
A: es el área de la cuenca en mi2.
UP: es el caudal pico del Hidrograma Unitario Sintético, en pie3/s/pul.
T
Tp = + TL (10.37)
2
En donde:
TL
tb = 3 + 3 (10.38)
24
Donde:
280
tb: es el tiempo base en días, con TL, en horas.
Esta ecuación da una estimación razonable del tiempo base para cuencas
grandes, pero produce valores excesivamente altos para cuencas pequeñas.
En éstas, el tiempo base puede calcularse, en forma aproximada, como 3 a 5
veces el tiempo al pico.
770 440
W50 = W75 = (10.39)
u1p.08 u1p.08
En donde:
up :es el caudal pico por unidad de área en pie3/s mi2.
W50 :es el intervalo de tiempo en horas, correspondiente al 50% del caudal
pico.
W75: es el intervalo de tiempo en horas, correspondiente al 75% del caudal
pico.
281
FIGURA 10.7 Hidrograma Unitario Sintético de Snyder.
.
El autor es muy claro al advertir que “el modelo tal como se dedujo, solo es
válido para las cuencas que fueron objeto de su estudio y que la aplicación
del mismo a cuencas de otros sitios debe ser precedida por la deducción de
los coeficientes Ct y Cp, y quizás hasta de las mismas ecuaciones” (Snyder,
1938).
282
sus tiempos característicos y representándolas como un Hidrograma
Unitario Triangular.
1
Tb = = 2.67 = 8 (10.40)
0.375 3
Tr = 1.67 Tp (10.41)
283
FIGURA 10.8. Hidrograma Unitario Adimensional del S.C.S.
Caudal pico. El caudal pico en pies cúbicos por segundo por pulgada se
puede calcular como:
484 ∗ A
Up = (10.42)
Tp
284
Cualquier modificación en el Hidrograma Unitario Adimensional que
conlleve cambios en el porcentaje del volumen de escorrentía bajo su rama
creciente produce variaciones en el factor de forma asociado al Hidrograma
Unitario Triangular, y por tanto la constante también cambia.
Para las cuencas consideradas por el SCS, el factor del caudal pico varió
desde 300, en terrenos llanos, hasta 600, en zonas de pendiente empinadas.
De lo anterior, se deduce que si se utiliza un Hidrograma Unitario
Adimensional diferente al derivado por el SCS, el factor de caudal pico
cambia de valor, y, por consiguiente, dicho caudal será distinto del que se
obtiene con la ecuación (10.42).
TC = 5 TL (10.43)
3
L0.8 ∗ (S + 1)
0 .7
TL = c (10.44)
1900 ∗ S c0.5
285
En donde LC es la longitud del canal principal en pies, SC es la pendiente
promedio de la cuenca en porcentaje y S es el factor de retención o
almacenamiento en pulgadas, definido por la ecuación 8.9
T
Tp = + TL (10.45)
2
TL = 0.6TC (10.46)
Los valores del tiempo de rezago según esta expresión son los mismos
obtenidos para el Hidrograma Unitario de Snyder.
T = 0.133TC (10.47)
286
• Se calcula la duración, T, con la ecuación (10.47).
• Se multiplica este tiempo al pico por cada una de las abscisas del
Hidrograma Unitario Adimensional del SCS (columna 1 de la Tabla
10.3), consiguiendo así las abscisas del Hidrograma Unitario en horas.
287
Nótese que el Hidrograma Unitario así derivado corresponde a una duración
de T-horas, la cual es función de TC.
n −1
⎛ t ⎞ ⎡ ⎛ t ⎞⎤
U = Up ⎜ ⎟ exp⎢(1 − n )⎜ − 1 ⎟ ⎥ para t ≤ t 0 (10.48)
⎜T ⎟ ⎢⎣ ⎜T ⎟⎥
⎝ p⎠ ⎝ p ⎠⎦
⎛t −t⎞
U = U 0 exp⎜ 0 ⎟ para t 0 < t ≤ t 1 (10.49)
⎝ K ⎠
⎛t −t⎞
U = U 1 exp⎜ 1 ⎟ para t > t 1 (10.50)
⎝ 3K ⎠
288
Según las tres ecuaciones anteriores, el Hidrograma Unitario Sintético de
Williams y Hann queda bien determinado de la siguiente manera:
0.124
⎛L ⎞
K = 27.0 ∗ A c ∗S − 0.777
∗ ⎜⎜ c ⎟⎟
0.231
cp (10.51)
⎝ Wc ⎠
0.133
⎛L ⎞
Tp = 4.63 ∗ A 0.422
c ∗S − 0.46
cp ∗ ⎜⎜ c ⎟⎟ (10.52)
⎝ Wc ⎠
289
FIGURA 10.9. Hidrograma Unitario Instantáneo de Williams y Hann
⎡ 1 ⎤
t 0 = Tp ⎢1 + 1 ⎥ (10.53)
⎣⎢ (n − 1) 2 ⎦⎥
290
t 1 = t 0 + 2K (10.54)
2
⎧ ⎡ ⎤
1
2⎫
⎪ 1 1 1 ⎪
n = 1+ ⎨ +⎢ + ⎥ ⎬
⎪ 2(K Tp ) ⎣⎢ 4(K Tp ) (K Tp )⎦⎥ ⎪
2
(10.55)
⎩ ⎭
B ∗ Ac
Up = (10.56)
Tp
291
RELACION ENTRE B Y n
1000
100
B
10
0 2 4 6 8 10 12
n
292
Ejemplo 10.3 (Tomado de Smith, Vélez, 1997)
Solución:
a) HU de Snyder.
• TL = 1.8730 h
• Tr = 0.3405 h
• TLC = 1.8704 h
• Up = 9.7511 m3/s /mm
• Tp = 2.0407 h ≈ 120 min
• Tb = 600 min
• W50 = 3.0285 h ≈ 180 min
• W75 = 1.7306 h ≈ 100 min
293
TABLA 10.4 Parámetros morfométricos de la cuenca del Río San
Carlos, hasta la estación Puente Arkansas
Parámetros Valor
2
Area de la cuenca, Ac (km ) 134.3
Longitud del río principal, Lcp (Km) 27
Longitud de la cuenca, Lc (Km) 16.5
Longitud al centroide, Lca (Km) 14
Altura promedia, ∆h (m.s.n.m) 1449
Pendiente del canal, Scp (%) 3.14
Pendiente del canal, Scp (pie/milla) 256.1
Relación Longitud-ancho de la cuenca, Lc/Wc 2.03
Número de curva, CN 77.95
Tiempo de concentración, TC (min) 180
b) HU del S.C.S
• S = 2.83
• TL = 1.8 h
• Tp = 1.965 h
• Tb = 9.825 h
• Up = 14.2378 m3/s /mm
294
HIDROGRAMA UNITARIO DE
SNYDER (AEI)
12
Caudal (m3/s)
Tiempo (min)
c) HU de Williams y Hann
Con base en las ecuaciones presentadas en la sección 10.4.3 y siguiendo el
procedimiento de cálculo presentado al final de esa sección, se obtienen los
siguientes resultados:
• K = 0.987
• Tp = 2.1 h ≈ 120 min
• K / Tp = 0.47
• n = 9.23
• B = 650
• Up = 17.9054 m3/s /mm
• t0 = 2.83 h ≈ 160 min
• t1 = 4.804 h ≈ 280 min
• U0 = 11.9266 m3/s /mm
• U1 = 1.6141 m3/s /mm
295
HIDROGRAMA UNITARIO DEL S.C.S
16
14
Caudal (m3/s)
12
10
8
6
4
2
0
0 200 400 600 800
Tiempo (min)
296
HIDROGRAMA UNITARIO DE WILLIAMS
Y HANN
18
Caudal (m3/s)
16
14
12
10
8
6
4
2
0
0 100 200 300 400 500
Tiempo (min)
297
Capítulo 11
ESTUDIO DE CRECIENTES
296
Una creciente es un evento que produce en un río o canal niveles muy
altos, en los cuales el agua sobrepasa la banca e inunda las zonas
aledañas. Las crecientes causan daños económicos, pérdidas de vidas
humanas y transtornan toda la actividad social y económica de una
región.
297
tiene el mismo período de retorno de la tormenta, lo que no siempre
corresponde a la realidad, pues tormentas de períodos de retorno bajos
pueden causar crecientes de períodos de retorno mucho mayores, si se
dan en el suelo condiciones extremas de humedad antecedente.
298
muestran en la tabla 11.1 los tiempos de concentración calculados para
la cuenca del río San Carlos, ejemplo 10.4
Expresiones Tc en minutos
Témez 172.02
Williams 292.07
Kirpich 147.60
Johnstone-Cross 147.71
C.C.P 154.81
Giandotti 171.13
S.C.S Ranser 155.89
Linsley 283.67
Snyder 276.90
Ventura Heras 172.02
Brausby-Williams 344.01
299
Precipitación efectiva.
CIA
Q= (11.1)
3.6
en donde:
Q : Caudal pico en m3/s.
C : Coeficiente de escorrentía.
I : Intensidad en mm/h.
A : Area de la cuenca en km2.
300
una lluvia con un período de retorno igual al tiempo de concentración,
Tc.
en donde:
QTR : Caudal en m3/s, para un período de retorno TR
301
CTR : Coeficiente de escorrentía, para un TR
ITcTR : Intensidad en mm/h, para un período de retorno, TR, y una
duración, Tc.
A: Área en Km2.
0
8 30'N
0
8 00'N
CAUCASIA
APARTADO
7 030'N
0.24
0.
24
7 000'N
DABEIBA
GUADALUPE
VEGACHI
0.20
6
0.1
6 030'N
0.2
0
URRAO BELLO
0.12
0.
12
0.16
MEDELLIN
0.0
8
6 000'N
FREDONIA
JARDIN
76 30'W
76 00'W
74 00'W
75 00'W
74 30'W
W
75 30'
0
0
0
0
0
0
5 30'N
302
Ejemplo 11.1
Area=0.1587 Km2
Longitud de la cuenca=1081.1 m
Pendiente del cauce= 21.76%
Longitud al centroide=552.44 m
Cota máxima =1912 m
Cota mínima =1777 m
AREAS URBANAS
AREAS AGRICOLAS
303
Solución:
Se puede observar que hay una banda de valores entre 0.142 y 0.179
horas, mientras los otros son demasiado altos o demasiado bajos, y no
corresponden a las características de la cuenca. Se promedian los valores
de la banda y se obtiene un tiempo de concentración de 11 minutos.
Como la duración de la lluvia de diseño debe ser igual o mayor al tiempo
de concentración, se propone una duración de 20 minutos.
304
b) Hidrógrafas unitarias sintéticas. Se calculan las hidrógrafas unitarias
de Snyder, SCS y William y Hann de acuerdo con las ecuaciones y la
metodología propuesta en el capítulo anterior. Las figuras 11.3, 11.4
y 11.5 muestran las HU de Snyder, SCS y William y Hann,
respectivamente.
Tiempo
(min) 5 10 15 20
acumulado
Tiempo
(%) 25 50 75 100
acumulado
P acumulada (%) 83.75 91.25 98.75 100
P acumulada (mm) 46.06 50.80 54.31 55
P efectiva
(mm) 1.05 1.86 2.59 2.74
acumulada
P efectiva
(mm) 1.05 0.81 0.73 0.15
Período
305
Después de aplicar la ecuación de convolución a cada una de las
hidrógrafas, se obtienen los siguientes valores:
Método Q m3/s
Snyder 1.56
SCS 1.91
Williams y Hann 0.28
Racional Probabilístico 1.85
Racional 5.29
Puede verse entonces la gran diferencia entre los resultados hallados con
cada método, lo que obliga al diseñador a una evaluación cuidadosa
antes de tomar una decisión final sobre el caudal de diseño. En este caso,
puede parecer más confiable el resultado hallado con el método racional
probabilístico, pues el coeficiente de escorrentía se encontró con datos
propios de la región, mientras en el caso de las hidrógrafas unitarias
sintéticas y el método racional tradicional, se trabaja con coeficientes
que no corresponden a la zona y hallados para climas y morfometrías
totalmente diferentes a la colombiana
306
FIGURA 11.4 Hidrógrafa de William y Hann para la q. El Limón
307
registros mas cortas, pero estiman los intervalos de confianza del valor
de interés (Mutreja, 1986)
308
11.3 DISEÑO DE CRECIENTES CON INFORMACIÓN ESCASA
309
Donde:
Donde:
µ̂ : media de los caudales máximos instantáneos, en m3/s.
σ̂ : desviación estandar de los caudales máximos instantáneos, en m3/s.
A: área de la cuenca, en Km2.
L: longitud del cauce principal, en Km.
S0: pendiente del cauce principal, en %.
Pr: precipitación promedia anual, en mm.
310
0
8 30'N
18
17
0
8 00'N 34
28
8
27
7
30 35
0
7 30'N
29
ZONA 1 32
78
68
26
3 70
33
83 82 79
84 62
5
0
4 42 80 ZONA 6
7 00'N
59
31 61 36
69
2 81 63
41 66
65
76 16
77 60
73 72 ZONA 7
15 35
45 50
5 46 33 34
9
4 46
24 58
6 ZONA 5 48
3 53 52
0
6 30'N 67 43
25
44 11
2 74 57
40 12
29
ZONA 2 7 71 64
11 56 13
14
47 31 30 17
15
8 75 27 24
12 32 18
11 16
26
54 10 20 25
ZONA 8
51 10
55 13
19 ZONA 4 28
21 23
0
6 00'N 23 22
37 14
22
21
39 ZONA 3
38 49 19
20
76 30'W
76 00'W
75 30'W
75 00'W
74 30'W
74 00'W
0
0
0
5 30'N
311
Para la aplicación del método se requieren una serie de pasos como:
312
1,80
1,60
1,40
1,2603 FX (x)
1,20 Q/Q 2,33 = 0,5307e
2
R = 0,9926
1,00
Q/Q2,33
0,80
0,60
0,40
0,20
0,00
0,010 0,100 1,000
FX(x)
313
Q TR = Q 10 + α q Ln( TR / T10 (11.7)
La media de los caudales máximos puede hallarse por medio del método
de la banca llena, como se explicó anteriormente, o con ecuaciones de
regionalización.
314
ARBOLETES
0
8 30' N
0
8 00' N
CAUCASIA
APARTADO
0
7 30' N
0 N
7 00'
DABEIBA
15
GUADALUPE 15
10
VEGACHI
20
20
15 25
0 N
6 30'
15
20
15
BELLO
URRAO
15 MEDELLIN 20
ENVIGADO SAN
0 N
6 00'
10
FREDONIA
JARDIN
76 76 75 75 74 74
030' W 000' W 030' W 000' W 030' W 000' W
0 N
5 30'
315
Los métodos Gradex y Racional producen resultados muy
similares a los hallados con el análisis de frecuencia, utilizando la
distribución Gumbel y Lognormal de dos parámetros.
316
Ca
pítulo 12
El AGUA SUBTERRÁNEA
317
El agua subterránea se utiliza en zonas rurales y aglomeraciones urbanas
para el abastecimiento de agua potable, y en proyectos agropecuarios para
riego y para uso animal; igualmente, muchas industrias consumidoras de
grandes cantidades de agua hacen uso de este recurso. Como ya se vió en
el capítulo 2, los recursos de agua subterránea son la mayor reserva de
agua dulce que tiene la tierra.
318
INGEOMINAS, 1993).
319
• Su protección contra la polución superficial, a la cual están
constantemente sometidas las aguas de ríos y lagos.
• Su estrecha relación con las corrientes de superficie, ya que éstas
obtienen de las capas acuíferas la mayor parte de su caudal.
320
piezométricos existente entre los mismos.
321
ocupados por líquidos, principalmente agua, y gases, en particular aire y
vapor de agua. La existencia de dichos vacíos constituye lo que se
denomina la porosidad. Además de la existencia de los poros, es
necesario considerar la interconexión que existe entre ellos y que
permite la circulación del agua, fenómeno que se denomina
permeabilidad.
- En las no consolidadas.
- En las compactas fisuradas.
- Continuidad o discontinuidad
- Isotropía o anisotropía
- Homogeneidad o heterogeneidad
322
Un medio permeable con vacíos interconectados, en el sentido del flujo,
es continuo. Las rocas porosas y las rocas compactadas, con sistemas de
fisuras o fracturas interconectados, constituyen medios continuos. Las
rocas con macrofisuras y las rocas kársticas son medios discontinuos.
323
12.1.1 Tipos de agua en el suelo.
- Agua de retención
- Agua gravítica
324
TABLA 12.1 Tipos de Aguas Subterráneas.
325
3) La transmisividad, T. Si se considera un acuífero cuya zona
saturada tiene una altura b, la transmisividad se define como el
volumen de agua por unidad de tiempo (o caudal) que pasa a
través de una sección vertical de ancho unitario y de altura b,
bajo el efecto de una unidad de gradiente hidráulico y a una
temperatura de 20 °C.
T=K b (12.1)
326
12.2 LA LEY DE DARCY
327
Darcy demostró que el caudal, Q, es proporcional a la pérdida de carga,
inversamente proporcional a la longitud del lecho de arena y proporcional
al área de la sección y a un coeficiente que depende de las características
del material.
Q =K A(h1 − h 2 ) / l (12.2)
328
ahora, ∆h/l no es otra cosa que la pendiente de la línea de carga, o sea la
pérdida de carga por unidad de longitud del lecho filtrante, la cual se llama
gradiente hidráulico. Este gradiente se puede simbolizar mediante la letra
i, entonces:
∆h
i= (12.4)
l
∆h
Q =− K A = − K Ai (12.5)
l
Q
V= = Ki (12.6)
A
- La permeabilidad, K.
- El gradiente hidráulico, i.
- La velocidad del flujo, V.
329
Estos parámetros serán estudiados detenidamente a lo largo de este
capítulo.
Ejemplo 12.1
Solución:
h 1 − h 2 ∆h
Donde el gradiente hidráulico es: i= =
L L
Ahora;
L sen 30o
i= ⇒ i = sen30°
L
330
0.005
k= ⇒ k = 0.01 cm/s
sen 30 o
Si V = Ki, en donde:
V: velocidad de flujo,
i: gradiente hidráulico,
γ
K =K 0 (12.7)
µ
En donde:
331
Si se considera la viscosidad cinemática (ν) del fluido, se tiene que:
µ
ν= (12.8)
ρ
en donde ρ es la densidad del fluido, y como γ = ρ g, entonces:
γ ρg g
= =
µ νρ ν
g
K =K 0 (12.9)
ν
K 0 =Cd 10
2
(12.10)
332
Arena fina 40 - 80
Arena media, bien gradada 80 - 120
Grava arenosa, pobremente gradada 80 - 120
Grava arenosa limpia, bien gradada 120 – 150
333
TABLA 12.2 Valores de la permeabilidad (Rodríguez, 1981).
Rocas K (m/día)
Grava limpia 1000
Arena gruesa limpia 10 - 1000
Mezcla de arena 5 - 10
Arena fina 1-5
Arena limosa 0.1 - 2
Limo 0.001 - 0.5
Arcilla > 0.001
V
Q=
t
∆h Ql
Q =K A K=
l A∆h
334
variable, en el que la muestra está sometida a cargas más altas, por medio
de un tubo largo de sección delgada, Figura 12.4b. En el permeámetro de
cabeza variable, se mide el descenso del volumen de agua durante un
tiempo t, y luego por integración se calcula el coeficiente K.
h
Q =K A
l
335
dh dh A Kdt
Q = −a =
dt h a L
ln
h
=−
AK
(t − t 0 ) (12.11.)
h0 aL
Si se dibuja en papel semilogarítmico log h vs t, se obtiene una recta en la
cual la pendiente es proporcional a K.
336
Los ensayos de permeabilidad por inyección se realizan principalmente
en terrenos no saturados. Si se tiene régimen permanente y H es la
diferencia de carga que corresponde a un caudal Q, la permeabilidad
está dada por:
Q
K =m (12.12)
H
Donde m es un coeficiente que depende de la forma y el diámetro de la
rejilla filtrante. Para que la ecuación 12.11 sea dimensionalmente
homogénea, es necesario que m tenga dimensiones L-1. Si la zona
filtrante es cilíndrica y está en el fondo del pozo, se tiene que:
α
m= (12.13)
D
l
α=
l 1 (12.14)
2π +
D 4
V = Ki V = Vr ne ne
K = Vr
i
337
FIGURA 12.5 Esquema típico para un ensayo de bombeo.
pA pB
h1 = + ZA h2 = + ZB
γ γ
338
Al nivel de h en cada punto se le denomina nivel piezométrico. Al tubo
que permite determinar la altura de la columna de fluido se le denomina
piezómetro y al lugar geométrico de todos los niveles piezométricos se le
llama superficie piezométrica.
− ∆h h 1 − h 2
i= = (12.15)
l l
∆h ∂h
i = − lim =− (12.16)
∆h→0 l ∂l
339
Para el caso de un acuífero libre, la superficie piezométrica coincide
prácticamente con el límite de la zona saturada; se le denomina superficie
freática. El agua en estos puntos está sometida a la presión atmosférica.
entonces:
− ∆h
i= ∆h = h A − h B
l
h A = Z1 − h1 h B = Z 2 − h 2
340
Lo anterior implica que:
∆h = (Z1 − Z 2 ) + (h1 − h 2 )
Como su nombre lo indica, las curvas isopiezas son líneas de igual altura
piezométrica; se les denomina también líneas equipotenciales. Estas
curvas son las que resultarían si se cortara la superficie piezométrica por
planos horizontales igualmente espaciados.
100 − 90
i= = 0.004
2500
341
Una gran cantidad de información útil para estudios hidrogeológicos, se
puede extraer de los mapas de isopiezas. Generalmente, los planos de
isopiezas se referencian a los planos topográficos, lo que permite
determinar la profundidad del agua subterránea, la dirección del flujo,
hacer inferencias sobre la permeabilidad de las formaciones, etc. La figura
12.9 muestra el mapa de isopiezas de una región del oriente
antioqueño.(Universidad Nacional, 1997)
..
342
sección a través de la cual circula realmente el fluido. Esa sección eficaz
será entonces:
S = An e (12.17)
Q = AK i (12.18)
343
344
REFERENCIAS BIBLIOGRÁFICAS
370
CAMPO, J.M., MUNERA,J.C. 1997. Determinación de tiempos de
concentración y coeficientes de escorrentía para algunas cuencas de
Antioquia. Tesis (Ingeniería Civil). Universidad Nacional. Medellín.
371
GEOLOGICAL SURVEY . 1983. Computation of Discharge. Vol 2.
U.S Geological Survey Water Suply Paper 2175.
GREEN, W.H., AMPT, G.A. 1911. Studies of soil physics, 1: The flow
of air and water through soils. J. Of Agriculture Science, vol 4, no 1, pp
1-24.
372
LANE, E.W. 1955. The importance of fluvial geomorphology in
Hydraulic Engineering,. Proc. ASCE, 81, paper 745, pp 1-17.
373
PHILIP, J.R. 1957. The theory of infiltration: I. The infiltration
equations and its solution. Soil Sci, vol. 83, pp 345-357.
374
SHAW, M.E. 1994. Hydrology in Practice. London. Chapman &
Hall.569 pp.
375
UNIVERSIDAD NACIONAL DE COLOMBIA.1997. Evaluación de
oferta y demanda hídrica para la subregión No. 1 del Departamento de
Risaralda. Medellín, 1997
VÉLEZ,M.V.,.SMITH.,R.A.,.PÉREZ,C.A.,FRANCO,C.J.,BOLAÑOS.,
H.D. 1995. Aplicación de algunas técnicas al diseño hidrológico con
información escasa / // En: Avances en Recursos Hidráulicos, No. 3.
p.35-54.
376
WORLD METEOROLOGICAL ORGANIZATION. 1970. Guide to
hydrometeorological practices. 2d.ed.,WMO no. 168, tech. Pap 82, pp
III-8 III-11. Geneve.
377