Hidrologia para Ingenieros: María Victoria Vélez Otálvaro
Hidrologia para Ingenieros: María Victoria Vélez Otálvaro
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
Todos o casi la mayoría de filósofos de la antigüedad centraron su
atención en los procesos que involucran la producción de agua superficial
y otros fenómenos que tenían que ver con la ocurrencia del líquido en sus
diversas etapas y de su ciclo perpetuo:: del mar a la tierra, y visceversa.
Infortunadamente, las especulaciones eran a menudo falsas. Homero creía
en la existencia de un gran lago subterráneo que alimentaba ríos, mares y
pozos. Es interesante observar que Homero estudió la dependencia del
flujo de la secci¢n transversal y de la velocidad. Este conocimiento se
perdió luego durante el Imperio Romano y las relaciones entre área,
velocidad y caudal permanecieron desconocidas hasta que Leonardo Da
Vinci las descubrió durante el Renacimiento.
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)
8
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 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
del sistema y especificar las entradas y salidas. La figura 2.2 muestra un
esquema de un sistema.
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
10
= I −Q dt
dS = I
Donde la variable I representa las entradas y Q las salidas.
− Q dt
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 (Q s). Considerando
un período suficiente largo, puede afirmarse que el balance global promedio
es:
P = I +ET+Qs (2.1)
11
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.
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
2.2 INTERACCION HOMBRE _CICLO HIDROLOGICO.
Físicos:
13
general evapotranspiran parte del agua que se infiltra en el suelo. Es así
como más de la mitad del agua lluvia sobre la selva amazónica se debe a la
evapotranspiración de arboles y plantas que allí existen.(UNESCO, 1991).
La tala de bosques causa también un aumento de la escorrentía superficial
directa, y produce erosión , lo que a su vez provoca problemas de
sedimentación en corrientes de agua, lagos y embalses. En estos últimos el
problema Es grave: en algunos países como la India, han tenido en 21
embalses, aportes de sedimentos 200% mayores que los inicialmente
calculados al diseñar la obra.(UNESCO, 1995). El aumento de los
sedimentos también puede disminuir la cantidad de oxígeno presente en el
agua, necesario para la vida de peces y fauna acuática. En Colombia esto se
ha visto en la desembocadura de los ríos San Jorge y Sinú, que debido a la
tala incontrolada en sus cabeceras han aumentado enormemente el caudal
sólido que transportan, afectando no 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
La lluvia ácida es un importante problema mundial. La acidez de algunas
precipitaciones en el norte de Estados Unidos y Europa es equivalente a la
del vinagre. La lluvia ácida corroe los metales, desgasta los edificios y
monumentos de piedra, daña y mata la vegetación y acidifica lagos,
corrientes de agua y suelos, sobre todo en ciertas zonas del noreste de
Estados Unidos y el norte de Europa. En estas regiones, la acidificación
lacustre ha hecho morir a poblaciones de peces. Hoy también es un
problema en el sureste de Estados Unidos y en la zona central del norte de
África. (Enciclopedia Encarta, Microsoft, 1998)
15
de agua subterránea constituyen casi el 99% del total de agua dulce del
planeta.
Dulce no 30 4 80
utilizable
16
Dulce utilizable 11 1 20 100
Agua subterránea 0.99 19.7 98.5
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
La atmósfera sirve como almacenamiento temporal del agua evaporada de
corrientes y cuerpos de agua, plantas y seres vivos, y como medio de
transporte en todas direcciones de esta humedad. Se puede dividir en 3
grandes zonas: ionosfera, estratosfera, y troposfera,que se distribuyen en
capas tal como muestra la figura 3.1. La troposfera es la capa más
importante de la estructura del clima, pues contiene alrededor de las 2/3
partes de los gases que conforman la atmósfera y en ella se presentan la
mayoría de los fenómenos climatológicos. Las moléculas de estos gases se
mueven libremente y el campo gravitacional de la tierra evita que ellas
escapen.
ionosfera
80 km
estratosfera
18 km
8 km troposfera
19
El vapor de agua de una muestra de aire, bajo condiciones muy húmedas
ocupa normalmente un 4% del volumen de la muestra y puede llegar a
alcanzar valores de cero o muy cercanos a cero. Además, puede pasar al
estado líquido o sólido y convertirse nuevamente al (estado) gaseoso. En
este último proceso las moléculas toman energía del medio que las contiene.
En el proceso contrario, liberan energía. Este efecto es uno de los más
importantes en el transporte de calor en el equilibrio térmico del globo.
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
La capa de ozono está compuesta de moléculas de oxígeno en la forma O 3
(llamadas ozono). El ozono es producido por la acción de los rayos
ultravioletas en los átomos comunes de oxígeno. Esta capa protege la tierra
de la mayor parte de los rayos ultravioletas, que de otra forma dañarían los
tejidos de animales y seres vivos en general. Recientemente se ha
descubierto un gran agujero (en esta capa) producido por la producción de
freón en la atmósfera. El freón es un gas sintético que contiene átomos de
flúor, carbón y cloro. Estos compuestos son llamados tambien halocarbonos.
Son usados en los aerosoles, en sistemas de refrigeración y aire
acondicionado. Dichos compuestos viajan hacia la troposfera y
eventualmente pueden alcanzar la estratosfera. Como ellos absorben
radiación ultravioleta se descomponen y producen cloro, que ataca las
moléculas de ozono, convirtiéndolas en moléculas de oxígeno,
disminuyendo así el espesor de la capa (de ozono) e incrementando la
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
21
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.
E
λ
= f(λ,T) (3.1) aλ
Donde:
22
La relación Eλ/aλexpresa la capacidad de emisión de un cuerpo. Los cuerpos
con a =1, para cualquier λ se denominan cuerpos negros. Para estos
últimos la ecuación 3.1 toma la forma de:
Eλ = f(λ,T) (3.2)
F =σT (3.3)
Donde:
23
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
En la mayoría de las estaciones meteorológicas la radiación se mide con el
heliógrafo , el cual está constituído por una esfera de cristal pivoteada en un
eje que hace las veces de lente. Los rayos solares al atravesarla chocan
contra una franja de cartulina que se encuentra sobre un semianillo
alrededor de la esfera , graduada convenientemente y situada a la distancia
focal de la lente, donde se marcan las horas. Los rayos solares (cuando hay
sol) queman la cartulina lo que permite evaluar el número de horas de brillo
solar en el día.
24
la atmósfera, una capa semiporosa compuesta de bióxido de carbono y
vapor de agua, ozono, metano y compuestos de flurocarbono.
3.2.2 La temperatura
25
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
e
ρv = 0.622 (3.4)
RT
Donde:
26
ρ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
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.
27
e
q = 0.622 (3.6)
pa
e
r = 0.622 (3.7)
pa e
28
La humedad se mide con un sicrómetro, que consiste en dos termómetros
iguales, uno con el bulbo en contacto directo con el aire y el otro con el
bulbo húmedo. Por algún método, se acelera la evaporación y debido al
enfriamiento producido por el calor latente, necesario para la evaporación,
el termómetro húmedo, que suministra la energía necesaria para el calor
latente de evaporación, marcará una temperatura menor que el bulbo seco.
Con estas temperaturas y tablas sicrométricas se puede encontrar la
humedad relativa. Otro instrumento que también se usa es el hidrógrafo de
cabello o de otras fibras: Con la variación de longitud que experimenta el
cabello al variar la humedad se puede hallar la humedad relativa.
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.
29
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
30
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).
31
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
402 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.
32
En las zonas polares se presentan focos de alta presión, con aire seco y baja
precipitación y evaporación
33
del Niño” para describir una invasión de aguas tibias provenientes del sur
del Perú, que ocurría con intervalos de pocos años, en la época de Navidad,
y que disminuía la pesca en las costas peruanas. Aunque las temperaturas
del mar y el aire empiezan a decrecer típicamente de sus valores máximos
en Abril o Mayo, ellas pueden permanecer anómalamente altas por un año o
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:
34
1902,1904,1911,1913,1918,1923,1925,1930,1932,1939,1951,1953,19571958,
1963,1963,1965,1969,1972,1976-1977,1982-1983,1986-1987,19911992,1997.
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
35
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
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
36
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).
37
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.
38
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
39
1983,1987,1991/1992, años todos que corresponden a eventos Niño, con
excepción de 1961, año que está determinado por el ENSO 1957- 1958-.
(Poveda 1998)
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.
40
más lluviosas del país. La cordillera Oriental recibe el resto de lluvias del
Pacifico y las que llegan del oriente a impulsos de los alisios. La figura 3.7
muestra el mapa de isoyetas de precipitaciones promedias anuales en
Colombia elaborado por el Ideam.
41
FIGURA 3.7 Isoyetas de precipitaciön promedia anual en Colombia
(IDEAM 1999)
42
43
Capítulo 4
GEOMORFOLOGIA DE CUENCAS
44
Las características físicas de una cuenca forman un conjunto que influye
profundamente en el comportamiento hidrológico de dicha zona, tanto a
nivel de las excitaciones como de las respuestas de la cuenca tomada como
un sistema. Así, pues, el estudio sistemático de los parámetros físicos de las
cuencas es de gran utilidad práctica en la ingeniería de la Hidrología, pues
con base en ellos se puede lograr una transferencia de información de un
sitio a otro donde exista poca información: bien sea que fallen datos, bien
que haya carencia total de información de registros hidrológicos, si existe
cierta semejanza geomorfológica y climática de las zonas en cuestión.
45
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.
46
FIGURA 4.1 Divisoria de aguas superficiales y de aguas subterráneas.
Para efectos de balance hídrico si se presenta una situación como la
mostrada en la figura 4.1, el área superficial puede ser mucho menor que el
área total contribuyente al caudal de un río. Si se presentan estructuras
geológicas que favorecen la infiltración de aguas de otras cuencas, es
necesario tener en cuenta estos aportes que pueden ser bastante
significativos.
47
La tabla 4.2 presenta las relaciones Q min/Qmax encontradas para algunas
estaciones limnigráficas localizadas en el departamento de Antioquia
(Colombia).
Q = C An (4.1)
TABLA 4.2 Relaciones entre Qmin /Qmax para algunas cuencas de Antioquia
Q=100.146.A0.716 (4.2)
48
Donde :
A: área de la cuenca en km2
Q: media de los caudales máximos instantáneos en m3/s.
49
cuencas tienen áreas de km2. Este parámetro se simboliza con la letra
mayúscula A.
A
W= (4.4)
L
50
4.1.3 Parámetros de forma de la cuenca
51
FIGURA 4.3 Hidrógrafas según la forma de la cuenca
A
Rf = L2b (4.7)
52
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:
nhh
S h = Lh
nv
h Sv =
Lv
53
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.
Se + Sn
S= x100 % (4.9)
2
54
FIGURA 4.4. Método para hallar la pendiente S, en una cuenca
55
Es la representación de la superficie, en km2 o en porcentaje, comprendida
entre dos niveles, siendo la marca de clase el promedio de las alturas. De
esta forma, con diferentes niveles se puede formar el histograma. El
diagrama de barras puede ser obtenido con los mismos datos de la curva
hipsométrica. Realmente contiene la misma información de ésta, pero con
una representación diferente, y da una idea probabilística de la variación de
la altura en la cuenca, figura 4.6.
0.25
0.20
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.
56
Con base en la forma del perfil altimétrico del río, se pueden inferir rasgos
generales de la respuesta hidrológica de la cuenca en su expresión de la
hidrógrafa, o sea, la variación del caudal con el tiempo. También los perfiles
se usan para estudios de prefactibilidad de proyectos hidroeléctricos,
producción de sedimentos, ubicación de posibles sitios susceptibles de
avalanchas, etc. Generalmente, cuencas con pendientes altas en el cauce
principal tienden a tener hidrógrafas más picudas y más cortas que cuencas
con pendientes menores.
57
manera que las áreas de corte o positivas y de lleno o negativas sean
iguales y mínimas.
El valor del orden del canal principal,ς, en la boca de la cuenca da una idea de
la magnitud del drenaje de la cuenca.
58
FIGURA 4.8 Orden de una cuenca
= Nn (4.10)
Rb
Nn +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.
donde:
59
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
Li = l1rli−1 (4.12)
Ai = A1rai−1 (4.13)
Donde:
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
60
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:
∑l i
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.
61
FIGURA 4.9 Patrones de drenaje
62
Aunque el coeficiente mencionado en último término no se podría
denominar como un parámetro geomorfológico, sí es interesante citarlo por
la importancia que tiene en el manejo de una cuenca.
63
a) Puede haber un decrecimiento de la contribución de las aguas subterráneas a
los ríos con el incremento de Dd, la densidad de drenaje.
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.
64
Las características no estacionarias de los diferentes ríos conforman una
gran gama en sus variaciones, con cambios continuos en el tiempo. Están
determinadas por parámetros tales como:
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
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.
65
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.
S∝ QQbs Dc b (4.15)
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
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
• Varios autores han desarrollado criterios para definir la sección a banca
llena. Ridley (1972), utilizó un índice “bench” para definir el máximo
quiebre en la pendiente de las bancas. Wolman (1955) sugiere usar la
mínima relación de ancho- profundidad.
• Raíces expuestas por debajo de una capa de suelo intacta, que indican
exposición a los procesos erosivos.
69
• El liquen o el moho que crece en las piedras de las orillas se trunca a un
nivel que está por encima de las aguas bajas y presenta un límite inferior
que corresponde al nivel que alcanza el caudal a sección llena.
1 2 ⋅Sf1
3 2
(4.19)
Q = ⋅A⋅RH 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
n : coeficiente de rugosidad de Manning
70
Figura 4.11. Determinación en campo de la sección llena (Dunne y Leopold,
1978).
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
71
El factor de vegetación C, mide los efectos de la cobertura vegetal. Toma
valores de 1- 0.01.Suelos con buena cobertura tienen un C=1, mientras un
suelo con sobrepastoreo puede tomar un valor de 0.1.
El parámetro P tiene en cuenta las prácticas de conservación del suelo. Zonas con
sembrados perpendiculares a la pendiente (mala práctica de
conservación) se consideran con P=1
Capítulo 5
72
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
73
análisis probabilístico se analizan posibles leyes de probabilidad que pueden
describir el comportamiento de las variables de la población. En el análisis
estadístico, se hacen inferencias sobre la variable (la población), usando la
muestra. Por ejemplo, cuando se calcula una media con observaciones
disponibles, se está infiriendo que la media calculada es la media de la
población, lo cual no necesariamente es verdad, pues esto dependerá de la
calidad de la información, del número de observaciones y otros aspectos.
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
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.
74
En resumen, el objetivo principal de la estadística en hidrología es obtener
información de los fenómenos hidrológicos pasados y hacer inferencias
acerca de su comportamiento en el futuro.
m
P(X = x) = (5.1)
n
entre 0 y 1, así:
0 ≤ Pi ≤ 1 (5.2)
.
La probabilidad de un evento cierto es 1:
∑P = 1 i (5.3)
i=1
2) Si X1 y X
B B
2B son eventos independientes y mutuamente excluyentes, B
entonces:
75
Dos eventos son independientes si la probabilidad de ocurrencia de uno no
se ve afectada por la ocurrencia del otro,. y se dice que son mutuamente
excluyentes cuando la ocurrencia de uno imposibilita la ocurrencia del otro.
Los axiomas anteriores permiten la definición de conceptos importantes.
Por ejemplo, si dos eventos X1 y X 2 no son mutuamente excluyentes, la
B B B B
X1 ) = P( X1 ∩ X2 ) (5.7)
P(
X2 P(X2 )
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
76
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í:
77
Por ejemplo, si un caudal de 8098 m 3 /s es excedido en promedio una vez cada
P P
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
78
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)
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:
Se sabe que:
AA ∩ B
P( B ) = 0.8 = P(
) (4)
P(B)
79
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(P b ) B B
P(Pb ∩ F)
P(Pb / F) =
P(F)
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:
80
P(C) = 20/1000= 0.002
d) La probabilidad de contaminación C se puede expresar como:
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 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.
81
diagrama de barras como el mostrado en la Figura 3.1, que representa la
precipitación promedio anual en una estación.
ni (5.10)
f s (xi) =
n
la cual es un estimado de P( x i - ∆x<X<xi ) , la probabilidad de que la
B B B B
82
La suma de los valores de las frecuencias relativas en un punto es la funció n U U U
acumulada de probabilidades.
de probabilidades f X ( x)
U U U U U B B
distribución de probabilidades F X ( x)
U U U U U U B B
dFX(x)
fX(x) = (5.14)
dx
83
Para un valor dado de la variable aleatoria X, Fx ( x) es la probabilidad
B B
FX (xA ) = P (X ≤ xA ) (5.16)
XA = µ + Kσ (5.17)
FX (XT ) = P (X ≤ XT)
= 1 - P (X > XT)
84
1
0.9
0.8
0.7
0.6
F X( 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 (XT ) = 1−
T
O:
1
FX (µ+σK) = 1−
T
Y se obtiene finalmente:
K = 1 FX−1⎛⎜1− 1 ⎞⎟−
µσ⎝T⎠
85
1)
(x)dx = 1 (5.19)
2)
b
3)
b
(5.21)
∫f (x)dx= 0
X
b
1)
(5.22)
∑f(x ) = 1 i
i
2)
xi ≤b
3)
i=j
86
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:
P(AB) = P(A) + P(B) = (1/2)2 P P
(1/2)x-1 P ½ = (1/2) P
x
P P se tiene
entonces que:
x P(X=x) Fx ( x) B B
1 ½ ½
2 ¼ ¾
3 1/8 7/8
probabilidades acumulada.
87
Las propiedades matemáticas de las distribuciones estadísticas pueden ser
definidas en términos de los momentos de la distribución.
Los momentos representan parámetros que tienen significado físico o
geométrico. Se reconocerá fácilmente la analogía entre los momentos
estadísticos y los momentos de área estudiados en mecánica de sólidos.
El r-avo momento con relación al origen se define como:
∞
o en el caso discreto:
n
µr (x)dx (5.27)
n
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.
88
5.4 CARACTERISTICAS ESTADISTICAS BASICAS
E du (5.29)
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:
ˆ 1
µ x= ∑ N
xi (5.31)
N i=1
89
- La varianza σ2 : mide la “variabilidad” de los datos, la dispersión de
P P
∑ )
σˆ x = N i=1 (xi −µx (5.33)
-1
90
E[(X- (x)dx (5.35)
1 3
γ = 3 E[(x-µ) ] σ (5.36)
N ∑(xi -µˆ x )3
γx
∧
i=1 = (5.37)
3
(N-1)(N-2) σˆ x
91
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.
γ<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:
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 yB
la frecuencia relativa es f x ( xi ) .
B B B B B
92
Se elabora entonces la tabla siguiente.
(mm) medio xi B B
(mm) fx ( xi )
B B B B
Σ=100 Σ=138.90
93
es 1-p .Supóngase que se está considerando un período de 3 años. La
probabilidad de excedencia en el año 3 y no en los años 1 y 2 es qqp, dado
que los eventos son independientes año a año. La probabilidad de
excedencia en cualquiera de los 3 años es pqq +qpq + qqp debido a que la
excedencia pudo ocurrir en el 1o., 2o o en el 3o. año. La probabilidad de
excedencia en 3 años está dada como 3q2 p. La probabilidad de dos
P P
n! x n−x
(5.38)
P(X = x) = p (1−p) x!(n − x)!
94
P(X = 1) = 50( 1 )1(1−1/50)49 = 0.37
3 50
50 0 50
año,P, es:
P = (1−p)n−1p (5.40)
95
1
µ= P
2 (1− P) (5.41) σ =
P2
Ejemplo 5.9
F(H)
5 6 7
H(m)
Solución:
96
P(H ≥ HTr=20 ) =1/ 20 = 0.05
HTr=20 = 6.9 m
b) Se puede escribir la siguiente ecuación:
P(HTr=20 sea excedida al menos una vez) =1 - P(H
B B
Tr=20 B no sea excedida) B
5 1 0.4
Ejemplo 5.9
97
El dique III tiene un caudal de diseño con un período de retorno de 25 años.
Asumir que las crecientes en los ríos A y B son estadísticamente independientes y
que las fallas de los diques I y III también lo son.
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:
98
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
99
basan en suposiciones de que la variable aleatoria se distribuye
normalmente. Es por ello que la Normal encuentre tantas aplicaciones en
hidrología: en pruebas de hipótesis, intervalos de confianza, etc.
(x−µx )2
1 − 2
e
fX (x) = σ 2π 2σ x
(5.42)
x
µ = (x−µx )/σx
1 u2
-
e 2
2π
∞
1 -w
2
(u ) = e 2
2π -∞ ∫
fu (u) =
(5.43)
Fudµ
10
0
FIGURA 5.5 Distribución normal.
10
1
5.6.1.1 Estimación de parámetros
1N
(5.44)
µˆ = ∑x i
N i=1
1 N
σˆ = = ∑(xi −µˆ)1/2
5.6.1.2 Factor de N (5.45)
frecuencia
x-µˆ
K= (5.46) σ
1 1
1. Fu (K) = 1− ⇒ K = Fu−1 (1− )
T T
10
2
⎛ 1⎞
2. Usando el valor calculado de P P
⎜1− ⎟ en la tabla 5.1, se lee el valor
⎝ T⎠
de x en la primera columna, que corresponde a K o F-1 µ (1- 1/T) P P B B
XT =µˆ +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 m 3 /s. ¿Si los datos se ajustan
P P P
Solución:
K = Fu -1 (0.99) B PB P
10
3
Supóngase, por ejemplo, que se desea estimar la media de la población, µ.
Asúmase que µ1 y 2 son dos estadísticos (funciones de la muestra aleatoria)
tales que: µ1 < µ2 y P(µ1< µ < µ2) =;. Entonces [µ1 , µ2] es llamado el
intervalo de confianza para la media µ., ; es llamado el nivel de confianza
(nivel de probabilidad) y µ1 y µ2 son llamados los límites de confianza
inferior y superior, respectivamente. Esta definición puede extenderse al
intervalo de estimación de un parámetro cualquiera o a una función del
parámetro.
XT ± u1-αST (5.47)
error estándar.
B
Cada distribución tiene expresiones para hallar el error estándar, por ejemplo,
el de la distribución Normal es:
x
σˆ
ST = N (1 + K / 2)
2
(5.48)
10
4
Ejemplo 5.12
año?. P
Solución:
1−µˆ
P(X ≤ 1) = P(µ ≤ )
σˆ
Ejemplo 5.13
10
5
probabilidad de que la escorrentía anual sea menor que 280 mm en todos los
tres siguientes años.
Solución:
y:
1 ( )
⎡ 1 y -µy 2 ⎤ fX(x) = exp ⎢ - 2
⎥ (5.49) σyx 2π ⎢⎣ 2 σy ⎥⎦
10
6
La forma de la distribución lognormal se muestra en la Figura 5.7.
1 [
⎧⎪ 1 ln (X-xo)-µy ] ⎫⎪
2
10
7
La distribución lognormal con tres parámetros puede aplicarse a eventos
con valores positivos o negativos, siempre que x≥x0 ; mientras que la B
lognormal con dos parámetros solo puede aplicarse a eventos con valores
B
positivos.
1
(5.51)
µˆ Y = ∑log (X )a i
N i=1
12
⎩⎝ 1N⎞⎠Ni=1
⎧⎛ ⎭
⎛ Cvx ⎞ (5.53) x0 = µx ⎜⎜ 1- ⎟⎟
⎝ Cvz ⎠
Donde:
σx Cvx
=
µx
(5.54) σz
Cvz =
µz
Donde:
( 1-w ) 2/3
Cvz = 1/3
w (5.55)
10
8
w= [ - γˆ + ( γˆ x x
2
+4 ) ];
1/2
γx
>0
σX , dadas como:
BB
(5.56)
µY = loga (µX )− 1 σY2
2
⎛ σ ⎞
σY =⎡⎢⎢loga ⎜⎜⎝1+ µXX2 ⎟⎟⎠ ⎤⎥⎥⎦ (5.57)
2
⎣
En este caso, se estiman µX y σX con los datos originales, y con las
B B B B
lognormal.
Ejemplo 5.14
Los caudales medios de un río en una estación hidrométrica han sido modelados
con las siguientes distribuciones:
a) Normal con parámetros µ = 256.7 m3 /s y σ= 191 m3 /s P P P P
Solución:
10
9
P(300[Q[400)= FX (400)- FX (300) B B B B
⎟ ⎜ ⎟
P(300≤Q≤400)= F µ ⎜⎜⎝⎛ B B 400σ-µx ⎟⎠⎞−Fu⎝ ⎛⎜ 300σ−x µx ⎠⎞ ⎟
x
= Fu ( u400
B B B ) - Fu ( u300
B B B B ) donde:B
y
u400 = (400 - 256.7)/191 = 0.75
B B
P(300≤Q≤400)=0.7734 - 0.5871=0.1863
y:
ln(300) = 5.704
ln(400) = 5.99
0.8186 B
11
0
Fµ ( µ5.704 )= (5.704 - 5.228)/0.84 = 0.564
B B B B
ln
(X )= µ + K σ
T y y (5.58)
en donde K F= u−1⎛⎜1− 1⎞⎟
⎝ T⎠
Si se quiere trabajar con la variable no transformada al campo logarítmico se
tiene que:
⎛ 1⎞
Fu−1 ⎜1− ⎟ es el inverso de la función de distribución Normal estandarizada
⎝ T⎠
11
1
acumulada y Cv es el coeficiente de variación
B B
en
donde:
σY (5.62)
ST = δ N
y
1/2
⎛
KT2 ⎞ (5.63) δ = ⎜⎜⎝1+ 2 ⎟⎟⎠
Ejemplo 5.15
Solución:
11
2
1 1
Fu(KT ) = 1- = 1- = 0.99 TI
100
De la tabla 5.1:
KT = Fµ−1(0.99) = 2.33
⎡ ⎡
⎢
⎣
(ln(1+0.33 )) 2 1/2
⎣
( 2 ⎦⎦
)
- ⎢ ln 1+0.332 ⎥⎤⎤⎥ - 1
exp 2.33
K=
0.333
K= 3.028
El valor asociado a un período de retorno de 100 años será:
XT = 15 + 5 x 3.028 = 30.14 m 3 /s
B B P P
ln(XT)±u1-α 2 ST
Se calcula primero δ con la ecuación (5.63) y luego ST con la ecuación (5.60), B B
el resultado es:
1/2
2
0.3246
ST =1.93* 30 = 0.11
Por lo tanto:
11
3
ln (30.28) ± 1.64 * 0.11
= 3.41 ± 0.1875
= [3.2225, 3.5975]
= [e3.2225 , e 3.5975 ] = [25.091, 36.5]
P P P P
6 (5.65)
αˆ = π σˆ
ˆ
β = µ -0.5772αˆ (5.66)
11
4
6
K = - π {0.577 +ln[lnTr -ln(Tr -1)]} (5.67)
σ
ST = δ N (5.69)
δ = [1+1.1396K+1.1K2]1/2 (5.70)
11
5
Esta distribución ha sido una de las más usadas en hidrología. Como la
mayoría de las variables hidrológicas son sesgadas, la función Gamma se
utiliza para ajustar la distribución de frecuencia de variables tales como
crecientes máximas anuales, caudales mínimos, volúmenes de flujo anuales
y estacionales, valores de precipitaciones extremas y volúmenes de lluvia de
corta duración. La función de distribución Gamma tiene dos o tres
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 ⎛ x ⎞β-1 -x
fX (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
dz (5.72)
β-1
11
6
|α |Γ(β) ⎝ α ⎠ ⎝ α⎠
Donde:
xo ≤ x < ; para ; > 0
B B
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)
11
7
Cˆ 2
v
µˆ (5.76)
αˆ =
βˆ
µ σ$ $, y C$ v son
la media, desviación estándar y coeficiente de variación
calculados con la muestra, respectivamente.
Para la distribución Gamma con tres parámetros o Pearson tipo III, los
parámetros, por el método de momentos, pueden estimarse por:
ˆ ⎛ 2 ⎞
β = ⎜⎜⎝ γˆ ⎟⎟⎠ (5.77)
αˆ = σˆ (5.78)
ˆ (5.79)
Xˆ 0 = µˆ −αˆβ
γ$ es el coeficiente de asimetría calculado usando la muestra.
Si se define:
⎛ 1⎞
KT = Fu ⎜⎜1- ⎟⎟ (5.80)
⎝ Tr ⎠
el factor de frecuencia K tiene la siguiente forma:
11
8
Para la distribución Pearson tipo III o Gamma de 3 parámetros, existen
tablas, como la 5.2, que dan el factor de frecuencia en función del
coeficiente de asimetría calculado con la muestra.
Si se tiene que:
XT± u1−α 2 ST
σ
ST = δ (5.82)
N
β-1
localización.
11
9
TABLA 5.2. VALORES DE KT PARA LA DISTRIBUCIÓN PEARSON
B B
12
0
0.0 0.000 0.842 1.282 1.751 2.054 2.326 2.576
ˆ
β = ⎛⎜ 2 ⎞⎟
⎜⎝ γˆ y ⎟⎠
γˆ y
αˆ = σˆ y (5.84)
2
yˆ 0 = µˆ y
ˆ
−αˆβ
12
1
YT = lnXT = µˆ 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
σˆ y (5.86)
ST = δ N
Los límite de confianza se pueden expresar como:
12ST
lnXT ±µ1−α/2 (5.87)
2
5.7 ANÁLISIS DE FRECUENCIA
California:
m
P= (5.88) n
Weibull:
m
P= (5.89) n +1
Hazen:
2m-1
P= (5.90)
2n
12
3
gráficamente en un papel de probabilidad diseñado para la distribución. En
este papel, las ordenadas representan el valor de x en una cierta escala y las
abscisas representan la probabilidad de P(X >x) o P(X< x), el período de
retorno o la variable reducida. Las escalas de las ordenadas y las abcisas son
diseñadas de tal manera que cuando una muestra es de una población con
esa distribución, la gráfica debe ajustarse a una línea recta. El propósito de
este papel es "linealizar" las relaciones de probabilidad para que los datos
puedan ser fácilmente dibujados y usados en extrapolación o propósitos de
comparación. Se puede observar en las páginas siguientes los papeles de
probabilidad correspondientes a las distribuciones Gumbel y Log-Normal.
Ven te Chow propuso que toda muestra se puede ajustar a una expresión como la
siguiente:
X = µˆ +K σˆ (5.91)
12
4
12
5
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.
12
6
de la función de distribución de probabilidades teórica, escogida , Fx (x), tal B
que: B
Dn = Max|FE(x)-Fx(x)| (5.92)
La prueba requiere que el valor D n calculado con la expresión anterior sea menor
B B
12
7
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
:
∑O =∑e i i
12
8
rechaza (para esos niveles de significancia). Si ocurre lo contrario, entonces
B
Ejemplo 5.16
Solución:
1944 77 77 77
1945 72 76 76
1946 76 78 74
12
9
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
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
13
0
77 6 31 0.6458 0.5753 0.0705
Ejemplo 5.17
Año Q m3 /s
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
13
1
1990 2664.9
1991 1392.7
Solución
Distribución Gumbel
σˆ = 1234.58 m3 /s P P
que:
ST =1111.458 m3 /s
B B P P
Para ;=0.05 se obtiene de la tabla 5.1 que T0.95 =1.645 y aplicando la ecuación
B B
Distribución Log-Normal
σˆ y = 0.4504
13
2
Aplicando la ecuación 5.59 para hallar el factor de frecuencia K y utilizando la
tabla 5.1 se halla:
K=2.055
De la ecuación 5.58: ln
QTr=50 B =8.8286 B
y sacando el antilogaritmo
:
QTr=50 =6827 m3 /s
B B P P
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
13
3
Capítulo 6
LA PRECIPITACIÓN
142
6.1 INTRODUCCION
143
- Enfriamiento de una masa por debajo del punto de condensación.
Este enfriamiento debe continuar hasta que (Ta), temperatura del
aire, sea inferior a la del punto de condensación o temperatura del
punto de rocío (Td).
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
pluviógrafo, que es el mismo pluviómetro provisto de un mecanismo de
relojería que le permite 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
6.6.1 Curva de masas de la lluvia.
10
9
8
Precipitación (mm)
7
6
5
4
3
2
1
0
1
13
15
17
19
21
23
11
Tiempo (h)
6.6.2 Hietograma.
148
Es el gráfico que relaciona la intensidad de la lluvia contra el intervalo de
tiempo. Se define la intensidad como la variación de la precipitación con el
tiempo. El intervalo de tiempo depende del tamaño de la cuenca. Para
cuencas pequeñas, se usan intervalos de minutos, y para cuencas grandes,
los intervalos son generalmente de horas.. El hietograma es muy utilizado en
el diseño de 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
13
15
1
11
Tiempo (días)
149
barras, como el mostrado por la figura 6.5, que presenta las precipitaciones
promedios mensuales multianuales de la Estación el Cedral en Pereira
(Risaralda)
150
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.
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 9101112
Meses
M N
N N
151
6.7.2 Consistencia de los registros.
Mc (6.2)
Pcx = Px Ma
donde:
152
Los registros obtenidios de un pluviómetro o de un pluviográfo representan
solamente un dato puntual dentro de un área determinada. Para los análisis
hidrológicos, se requiere conocer la precipitación sobre una región
cualquiera, por lo que se hace necesario convertir los valores puntuales de
varias estaciones en un valor promedio para esta región. Para esto, existen
tres procedimientos, que son los siguientes:
5
10 *
4.5
la estación X en mm
4
1971
4
3.5 1969
3 1967 a
2.5 c
2 19 65
1963
1.5 1961
1
1959
0.5 1957
0
13
15
11
1
153
6.8.2 Polígonos de Thiessen.
∑P i Ai
P = i=1 (6.4)
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
154
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.
155
AT
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.
156
estaciones pluviométricas son casi inexistentes a alturas mayores de 3500
msnm por lo que estimaciones de precipitación, hechas con este mapa para
estas alturas no son muy confiables.
2500
2000
Altura (m)
1500
1000
500
1600 1800 2000 2200 2400 2600 2800 3000 3200 3400
Precipitación (mm/año)
157
Nota: al oriente de la línea punteada la información no es confiable por escasez de estaciones
158
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í:
P = Po e-K A
n
(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.
159
Precipitaciones máximas); el área contenida dentro de cada isoyeta se
determina y se puede dibujar luego, un gráfico de Lámina - Area - Duración.
kTr m (6.7)
i = (c+ n
d)
160
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.
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
15min hora h
m
P=
161
n +1
7) Se construye la curva.
162
FIGURA 6.13 Curva I-D-F para la estación la Rápida (Antioquia)
(Smith,Vélez, 1997)
163
En los modelos lluvia escorrentía también es importante determinar la
distribución de la lluvia de diseño en el tiempo. Tradicionalmente para este
propósito se han utilizado los diagramas de Huff (1967) obtenidos por el
investigador del mismo nombre en Norteamérica,ver figura 6.15
164
de Medellín (Smith, Vélez 1997). La figura 6.16 muestra la distribución de
la lluvia en el tiempo para la estación Boquerón.
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
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
165
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)
K = XMax - P (6.9) σ
donde:
XMax :es la lluvia máxima observada de 24 horas.
166
&P&& :promedio de la precipitaciones máximas observadas. σ
:desviación estándar de las precipitaciones máximas observadas.
167
Capítulo 7
ANÁLISIS DE CAUDALES
7.1 GENERALIDADES
165
El régimen de caudales de una corriente de agua durante un período
determinado es el único término del balance hidrológico de una cuenca que
puede ser medido directamente con una buena precisión. Los otros
elementos de ese balance, como las precipitaciones, la evaporación, etc, no
pueden ser estimados sino a partir de mediciones observadas en distintos
puntos de la cuenca o deducidos de fórmulas hidrológicas, los cuales dan
siempre estimativos muy aproximados. El estudio del régimen de caudales
es indispensable, para los todos los diseños hidráulicos y para muchas obras
civiles en los que ellos son parte importante como las carreteras, puentes,
acueductos, presas, etc. Así, el preámbulo de todo estudio hidráulico de una
cuenca es la instalación de muchas "estaciones de aforo" que permitan
observar, en una serie de años tan larga, como sea posible, los caudales
escurridos en puntos característicos del río principal y, si fuere oportuno, de
sus diversos afluentes. Si embargo, en países como el nuestro, las estaciones
de aforo de caudales son inexistentes en muchos sitios, lo que ha obligado a
recurrir a métodos aproximados para la estimación de los caudales de
diseño, como son los métodos de regionalización. Empero, jamás debe
olvidarse que ningún método por bueno que sea reemplaza la medida directa
de la variable
166
Los métodos para medir caudales pueden clasificares en dos grandes
categorías: métodos directos y métodos indirectos. En estas dos categorías,
los más utilizados son:
Métodos directos:
Método área velocidad
Dilución con trazadores
Métodos indirectos:
Estructuras hidráulicas.
Método área pendiente.
167
Una de los procedimientos más comunes empleados en este método es el
descrito a continuación.
168
y el caudal Qi , correspondiente a la respectiva área de influencia, Ai, es:
Qi = &V&&iAi
QT =∑Qi (7.2)
i=1
169
FIGURA 7.2 Inyección de un volumen conocido de trazador
V1C1
Q= t
(7.3) ∫
t
170
qC o
C2
1
2
despejando el caudal Q :
Es importante anotar que para aplicar este método, se supone que el flujo es
permanente.
171
Los trazadores son de 3 tipos:
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?
Solución:
Q=
Q = 113.6 m3/s
172
7.2.2 Métodos indirectos
Q = CHn (7.5)
173
Y1
Y2
Z1
l
Z2
NR
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/3S1f /2 (7.7) VA = Q = RH
n
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
hf = h1 −h2 Sf =
L
174
2) Si en la fórmula de Manning :
1 2/3A
K = RH 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= K1K2 (7.9)
1
6
• Meyer - Peter & Muller, 1948 n = 0 .038 ⋅ D 90 (7.10)
1
= 0.0411 ⋅ D 65
6
• Raudkivi, 1976 n (7.11)
175
• Garde & Raju, 1978;
1
(7.13)
Subramanya, 1982 n=
0.047D506
• Bray, 1979 n (7.14)
⎜⎝ k ⎟⎠⎥⎦ (7.15)
a = 5.7798k0.139633 (7.15a)
1
R H6
n= 8 ⋅g (7.15b)
f
En estas ecuaciones, D50, D65 y D90 son diámetros característicos del 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.
176
realiza el conteo aleatorio de granos, según procedimiento ideado por
Wolman (1954). Este procedimiento es el siguiente:
3. En cada intercepto, se mide la cara expuesta más larga del grano que allí
se encuentre.
n = 0.0487⋅D50 16 (7.17)
donde:
Ejemplo 7.2.
177
Durante una creciente, las profundidades del agua en un canal rectangular de
10 m de ancho fueron 3 y 2,9 m, en dos secciones separadas por 200 m. La
pendiente del canal es 0,0001. Si n = 0,025, estimar el caudal. Recordar que
el radio hidráulico, RH ,es el área, A, sobre el perímetro mojado, P.
Solución:
Y1=3 m Y2=2.9 m
A1=30 m2 A2=29 m2
P1=16 P2=15.8 m2
RH1=1.875 RH2=1.875
=1824.7 =1738.9
en donde :
So: pendiente del canal.
Se calcula el K, así:
K = K1K2 =1781.3
178
El valor de hf se halla en la última columna con la ecuación 7.6; con ese
valor se empieza la proxima iteracción.
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
179
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.
Alturas de mira en m
29
28
27
26
25
24
23
0 5 10 15 20 25 30
Caudal en m 3 /s
180
14
12
Caudal m 3/s
10
8
6
4
2
0
0 0.2 0.4 0.6 0.8 1
Altura de mira m
V1 + y1 + z1 = V22 + y2 + z2 +hf
2g 2g
Si:
Z=z+y
Y despreciando cabezas de velocidad se obtiene:
Z1 − Z2 −hf = Fn (7.17)
181
yn
M1
y1
y2
1
Estaci ón Mira auxilia r
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
182
Nivel
bajando
Flujo
H
permanente
Nivel subiendo
QM = 1+ 1 dh (7.19)
Qn VwSo 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)
183
La mayoría de los diseños hidrológicos para estructuras hidráulicas,
necesitan considerar los caudales máximos extremos. Por razones obvias, la
medición directa de estos niveles y caudales extremos, rara vez se puede
realizar, por lo que se hace necesario extrapolar la curva de calibración para
hallar los caudales que correspondan a estos niveles. Existen varios métodos
para hacer esta extrapolación. Los dos más utilizados se presentan a
continuación: método logaritmico y método de Manning.
Q = C(H − H0 )n (7.21)
donde:
Q: caudal
H: nivel medido en la mira
H0 : nivel cuando Q es cero C
y n : constantes.
d) Se encuentran C y n
184
d) Se calcula Q para el valor deseado de H
2 /3
A RH
185
Para el diseño de estructuras hidráulicas y en general de obras relacionadas
con el agua, se trabaja con una serie de términos referidos al caudal que es
necesario conocer. Los principales son:
Caudal medio mensual, Qm. Se calcula hallando para cada mes la media
aritmética de los caudales promedios diarios.
186
FIGURA 7.10 Caudal promedio diario
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 de datos, la probabilidad de excedencia , P, de cualquier descarga( o
valor de clase), Q, es:
m
P= ×100 (7.23)
N
187
FIGURA 7.11 Curva de duración
188
% Tiempo
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.
189
El método para hallar esta curva regional es comparar gráficamente las
diferentes curvas de duración existentes en la zona, adimensionalizadas por
el caudal promedio diario correspondiente.
Q
Z=
Q (7.24) medio Donde:
Z: Caudal adimensional
Q: Caudal registrado
Qmedio: Caudal promedio diario multianual
10
190
8
Q/Qmedio
6
0
0 10 20 30 40 50 60 70 80 90 10
p(%) 0
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
191
Total 365 365 366 N=1096
140
12
Q m3/s 0
100
80
60
40
2
0
7.6.CURVA DE MASAS
V = ∫Qdt (7.25)
to
192
cualquier dV punto representa el caudal, Q, para un intervalo de tiempo
determinado.
dt
La diferencia entre dos puntos cualquiera de la curva es el volumen
almacenado, S, para ese período de tiempo, asumiendo que no hay pérdidas
en el embalse, S1 y S2 son los volúmenes de embalse requeridos para un
caudal de diseño determinado durante dos épocas de sequía.
Ejemplo 7.4
193
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:
194
Volumen acumulado en Mm x 10
4500
3
4000
3500
3000 Volumen
de embalse
2500
2000
1500
1000
500
0
105
135
166
196
227
258
288
319
349
15
45
74
Días acumulados
195
Capítulo 8
LA INFILTRACIÓN
8.1 INTRODUCCION.
196
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.
Vv 100 (8.1)
n= ×
Vt
Donde:
Vv: volumen de vacíos
VT: volumen total de la muestra.
197
Mármoles 0.1 - 0.2
Esquistos 1 - 10
Dolomitas 2.9
Granitos 1.5 - 0.02
Basaltos 0.1 - 2.9
⎛ V ⎞
Vv = P2 - P1 ⇒ n = ⎜ 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
198
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
Si se usa el factor F, la geofísica propone la expresión de Archie para
encontrar la porosidad total
C
F= nm (8.3)
199
• La pendiente, que permite saber la uniformidad del material: a mayor
pendiente más uniforme es la granulometría y más infiltración se
presenta. El material correspondiente a la curva 2 (figura 8.1), material
no gradado, favorece más la infiltración que el material correspondiente a
la curva 1, material bien gradado.
• 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.
200
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.
201
La capacidad de campo y el punto de marchitez dependen
obviamente del tipo de suelo y de la vegetación. Por
ejemplo, en la zona del cercano oriente antioqueño
(Rionegro, La Ceja, Marinilla, Guarne, Santuario), se tiene
un perfil de suelo, tal como se muestra en la Figura 8.2
(Universidad Nacional de Colombia, 1997). En este tipo de
suelo, la capacidad de campo es del 40.36 % y el punto de
marchitez el 22.96 %
( porcentajes respecto a una lámina infiltrada unitaria).
8.2.3 La vegetación.
202
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 siguiente relación para determinar la capacidad de infiltración:
f = fc + (f0 −fc )e−Kt (8.4)
203
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.
204
FIGURA 8.4 Infiltración, según Horton
f = (0.5)At−1/2 + B
(8.6)
F = At + Bt
1/2
205
ω∆θ
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.
8.4.1 Infiltrómetros.
206
Estos se usan en pequeñas áreas o cuencas experimentales. Cuando hay gran
variación en los suelos o en la vegetación, el área se divide en pequeñas
zonas uniformes y en cada una de ellas se realizan mediciones. Los
infiltrómetros son de dos tipos: tipo inundación y simuladores de lluvia.
I = P−ESD (8.8)
Ejemplo 8.1
207
FIGURA 8.5 Infiltrómetro tipo inundación
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 =ΠD2 /4 = 962cm2
V
f=
A∆t
∆t min 2 3 5 10 10 30 30 60
208
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
Tiempo min
209
FIGURA 8.7 Método del índice Φ.
N = 0.83 A2 (8.9)
A es el área en Km2.
210
FIGURA 8.9 Separación del flujo base por método b.
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:
211
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)
212
La escorrentía total es la suma de los anteriores valores (5.8 cm), lo que
implica que el índice φ de 0.55 cm/h es el correcto.
213
2
I )
(P - a (8.10)
P e = P - Ia +
S
⎛1000 ⎞
Ia = 0.2 S (8.12)
214
TABLA 8.3 Valores de CN para diferentes condiciones de humedad del
suelo
CN AMC II CN AMC I CN AMC III
100 100 100
95 87 98
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
215
Tabla 8.4. VALORES DE CN PARA AMC II.
---------- >1 68 79 86 89
---------- <1 39 61 74 80
Pastizales Contorneo >1 47 67 81 88
Contorneo <1 6 35 70 79
---------- <1 30 58 71 78
-
Pradera ----------
permanente ----------
----------
Bosques - 56 75 86 91
naturales - 46 68 78 84
Muy ralo - 36 60 70 77
---------- - 26 52 62 69
Ralo ----------
Normal - 15 44 54 61
Espeso
Muy espeso
- 72 82 87 89
Caminos - 74 84 90 92
De terracería
De superficie dura
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 evaporación generalmente se expresa en unidades de milímetro o
centímetro por hora, por día, por mes o por año.
- 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
efecto del viento sobre una superficie de agua depende también de la
profundidad de la masa de agua.
223
medida que la superficie de agua en los suelos esté más profunda. Por ello,
la tasa de evaporación desde la superficie del suelo está limitada tanto por la
disponibilidad de agua, como por la oportunidad de evaporación desde
aquél. Este comportamiento es muy lejano al de la evaporación desde una
superficie libre de agua en contacto con la atmósfera.
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.
224
contener, y a mayor temperatura del agua, mayor facilidad para la
evaporación. Generalmente, la evaporación es mayor en climas tropicales y
es muy baja en las regiones polares (Shaw,1994)
9.2.3 El viento.
9.3 LA TRANSPIRACION
225
hojas, a la atmósfera. La transpiración es gobernada por las mismas
relaciones físicas que rigen la evaporación y afectada por los mismos
factores, además de otros que son característicos del tipo de vegetación.
Desde el punto de vista hidrológico, se puede tomar como un proceso
similar al de la evaporación, excepto que la fuente no es una superficie de
agua.
226
profundidad del nivel freático. Lo que significa que para ellas, el consumo
de agua puede doblarse si la profundidad del nivel freático pasa de 140 cm a
70 cmm. El efecto que tiene la calidad del agua en la transpiración varía
según los tipos de planta. Muchas plantas usan menos agua si la cantidad
de sólidos disueltos aumenta.
227
para transformar parte de la información de los registros de temperatura y
aumentar así, la información sobre la evaporación.
A) Balance hídrico
B) Balance energético
C) Técnica de transferencia de masas
D) Evaporímetro o tanque evaporador
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
228
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.
229
Energía para Radiación
evaporación Radiación solar directa
reflejada
Intercambio de
calor
230
Los métodos de balance energético y transferencia de masas requieren datos
y sobre todo una buena instrumentación, que está disponible la mayoría de
las veces. Por esta razón, se han desarrollado una serie de fórmulas
empíricas, basadas casi siempre en el método de balance energético y
transferencia de masas y en el llamado método aerodinámico, en el que la
tasa de evaporación se puede obtener en función de la velocidad del viento,
que es el factor que transporta el vapor de agua lejos de la superficie del
agua, así:
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+ )(ew −ea ) (9.4)
10
donde:
E: evaporación en pulgadas por día.
Existen tablas para hallar los valores de ew y ea, según la temperatura, como
la tabla 9.1.
TABLA 9.1 Presión de vapor a varias temperaturas
231
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:
232
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
233
Para estimaciones de la evaporación por este método, para intervalos menores
de un año, pueden presentarse errores serios en los resultados. Así, pues, las
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
234
Er=0.8 (0.09083x607031)=44109.3 m3
235
directas con lisímetros y parcelas experimentales, métodos teóricos y
métodos empíricos.
9.6.1 Lisímetros.
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
236
períodos mayores de 10 días, para que el ∆S sea significativo. En general,
se usan períodos de un mes o mayores.
Donde:
237
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
ew = 4.584exp( (9.8)
237.3+ t)
H
n
n = Ha (1+ r)(a + b Nn )−σTa4(0.56− 0.092 ea )(0.10+ 0.90 N ) (9.9)
238
Donde:
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.29cosφ
239
Ta: temperatura media del aire en grados Kelvin,K. ( K = 273 + C)
ea y ew: presiones de vapor del aire y del agua, respectivamente.
Ea: parámetro que tiene en cuenta la velocidad del viento y el déficit de
saturación. Se expresa como:
P
EVPR = 2 para (P/L)∃0.316 (9.11)
P
0.9+ L2
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:
240
TABLA 9.4 Valores promedios de horas de sol, N
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
EVP = 1.6La ( &&)a (9.13)
IT
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
241
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
(9.14)
i = ∑i
1i =
(T&& /
5)1.514
a: constante empírica:
EVP = 2.54KF
(9.16)
F = ∑phTf /100
donde:
242
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
243
HR: humedad relativa media diurna, hallada en función de la humedad a las 8
y a las 14 horas.
244
EVP = 0.473HaC
(9.20)
C = CTCHCuCsCECM
Donde:
TC )+ 0.0476(TC )2 CT
= 0.393+ 0.5592(
20 20
U U 2
(9.21)
Cu = 0.708+ 0.3276( )−0.036( )
96.56 96.56
Hm )+ 0.12(Hm )2 −0.0218(Hm )4
CH = 1.25+ 0.348(
40 40 40
Lat. E F M A M J J A S O N D
N
245
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
H H
CH = 1.25−0.212( R
)−0.038( )
R 5
57.4 57.4
(9.22)
CS = 0.542+ 0.64(S/80)−0.4992(S/ 80) + 0.3174(S/ 80)
2 3
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.
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.
246
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.
247
9.7 LA EVAPORACION EN COLOMBIA
- 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
cordilleras, en niveles superiores a los 3.000 metros sobre el nivel
del mar.
248
entre 300 y 375 mm en el mes. El mes de agosto ocupa el segundo
lugar en cuanto a la magnitud de la evaporación media mensual.
T = 28.3079−0.0056517h (9.23)
donde:
T: temperatura promedio diaria anual en ºC
h:altura sobre el nivel del mar en m.
249
Capítulo 10
RELACIONES PRECIPITACIÓN-
ESCORRENTÍA
250
10.1 INTRODUCCION
251
El hidrograma puede ser mirado como la expresión integral de las
características fisiográficas y climáticas que gobiernan las relaciones entre
precipitación y escorrentía para una cuenca particular. Los factores
climáticos que afectan la forma y el volumen de la hidrógrafa son: la
intensidad y duración de la lluvia y su distribución espacial y temporal sobre
la cuenca. Como factores fisiográficos de más influencia en la hidrógrafa, se
señalan: el área y la forma de la cuenca, la naturaleza de la red de drenaje, la
pendiente de la cuenca y del canal principal.
252
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)
253
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:
254
Otros métodos para obtener la curva de recesión, junto con sus ventajas y
desventajas pueden consultarse en Blandon y Saldarriaga 1997.
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.
255
El hidrograma de escorrentía superficial directa (obtenido por cualquiera de
los métodos explicados en el capítulo 8), presenta la forma que se muestra
en la figura 10.3, en donde se pueden definir los siguientes parámetros:
Donde:
256
Tb :tiempo base
Tc = Tb − D (10.2)
257
Otras definiciones parten del tiempo de rezago, TR. Este parámetro es el
intervalo de tiempo entre el centro de gravedad del pluviograma de lluvia
efectiva y el centro de gravedad de la hidrógrafa de escorrentía superficial
directa. Otros autores reemplazan el centro de gravedad de la hidrógrafa por
el tiempo al pico, TP.
Tc = Tb − Tp (10.4)
258
A
L 0.4 (10.6)
• Williams (1922) Tc = DSo0.2
Tc = 0.066 ⎜⎟ (10.7)
⎝ So ⎠
⎛ ⎞ (10.8)
Tc=5 ⎜ ⎟
⎝ So ⎠
⎝H⎠
260
3
Lc
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⎞
FL
• Bransby - Williams Tc= A0.1So0.2 (10.14)
261
Tc : tiempo de concentración, en horas.
F : 58.5, si el área está en kilómetros cuadrados.
A : área de la cuenca, en kilómetros cuadrados. L
: longitud del cauce principal, en kilómetros.
So : pendiente del canal, en m/km.
L
• Pérez Monteagudo (1985) Tc = (10.15)
VR
VR se expresa como:
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.
262
TABLA 10.1 Tiempos de concentración (h), calculados por diferentes
métodos, en algunas cuencas de Antioquia (Smith, Vélez, 1997)
263
G-8 4.55 4.48 4.56 3.89 5.35 4.80
RG-5 4.83 5.56 4.80 4.48 6.82 5.18
RMS-17 2.23 1.69 2.22 1.93 2.20 1.77
RMS-11 1.05 1.02 1.06 1.15 1.86 0.63
RMS-14 10.04 6.43 9.87 6.43 5.11 13.88
PP-1 3.60 2.12 3.56 2.35 2.48 3.58
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:
TC = 9.88744A−0.093L0.334S−0.544 (10.17)
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 %
264
puede ser una hora, un día, intervalos de minutos, etc. Lo más importante es
que la lluvia sea constante y uniforme en un período unitario. Si la
intensidad de la precipitación efectiva es constante e igual a una unidad en
un tiempo unitario, la lámina de escorrentía directa producida será también
unitaria.
La teoría del hidrograma unitario se empezó a desarrollar con estudios
intensivos, encaminados a producir caudales de crecientes a partir de
registros de lluvias.
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:
265
dos precipitaciones efectivas de la misma duración, pero de distinta
intensidad, cada una producirá una hidrógrafa unitaria proporcional a la otra,
tal como lo muestra la figura 10.4. Por ejemplo si en la figura la lámina de
precipitación P2 es dos veces la lámina P1, entonces las ordenadas de la
hidrógrafa unitaria 2 serán el doble de las de la 1.
266
hidrológicos son seleccionados con cuidado, esta teoría es aceptable para
propósitos prácticos.
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
Q1 = P1Y1
Q2 = P1Y2 + P2Y1
Q3 = P3Y3 + P2Y2 + P3 Y1
. (10.18)
.
.
Qn PmYn−m+1 (10.19)
267
Matricialmente, podría expresarse como:
Q = PY (10.20)
N = L+ M −1 (10.21)
268
FIGURA 10.5 Superposición de HU
269
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.
270
VE: volumen de escorrentía superficial directa.
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)A
Q=(
VT t
pero:
VT = AL, lo que implica que:
VELA VE (10.23)
Q= =
(LA)t t
Vu L A = (1mmA) L A Y
=
VT t (LA) t (10.24)
1mmA
=
t
271
Dividiendo la ecuación 10.24 por la 10.23 y despejando Y, se tiene:
Q
Y= (10.25)
(VE / A)×1mm
272
superposición, el hidrograma de escorrentía directa está dado por la integral
de convolución como:
1t
=
Q(t) ∑U(T,i)P(t-i +1)Dt (10.29)
i=1
273
en donde N representa lo memoria del sistema.
10.4 LA CURVA EN S
274
Las ordenadas de la curva en S serán (a partir de un hidrograma unitario
cualquiera):
S1 = Y1
S2 = Y2
S3= Y3
.
(10.30)
.
.
n
Sn Yi
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
275
t (h) 8 9 10 11 12 13
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.
276
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)
277
ORDENADAS DE LA HIDRÓGRAFA UNITARIA
1 110 110 0
2 98 98 0
3 220 95.4 124.6 83.1
4 512 92.9 419.1 279.4
5 585 90.3 494.7 329.8
6 460 87.8 372.2 248.1
7 330 85.2 244.8 163.2
8 210 82.6 127.4 84.9 9 150 80.1 69.9 46.6
10 105 77.5 27.5 18.3
11 75 75 0
12 60 60 0
13 54 54 0
Qn PmYn−m+1 y
se tiene:
278
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
279
= 18.3 x 2 +46.6 x 1.5
= 106.5 m3/s
Q11 = Y8 P4
= 18.3 x 1.5 = 22.5 m3/s
Ejemplo 10.2
Solución:
280
Q p=2.8 m/s
40 min
100 min
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
281
Qn PmYn−m+1
se tiene:
282
100 0 14 13.5 0.25
14 14 0
Las hidrógrafas unitarias sintéticas pueden ser de tres tipos (Chow, 1994):
283
• Snyder
• Servicio de Conservación de Suelos (SCS)
• Williams y Hann
284
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.
TL = L19000c.8(S*+S1c)0.05.7 (10.31)
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)
TL (10.33) ts =
5.5
285
Donde ts es la duración de la lluvia efectiva en horas.
Caudal pico. El modelo propone calcular el caudal pico por milla cuadrada,
up, como:
640
up = Cp (10.34)
TL
Donde:
up: caudal pico del Hidrograma Unitario por unidad de área, en pie 3/ s.mi2 ,
Cp : coeficiente.
TL :l tiempo de rezago.
640
up = Cp (10.35)
[TL +(T− ts )/ 4]
286
Una vez obtenido el caudal pico por unidad de área de la cuenca, el caudal
pico total se obtiene como:
Up = upA (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 (10.38)
tb = 3+ 3
287
24
Donde:
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 = u1.08 W75 = u1p.08 (10.39)
p
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.
288
El punto C en la Figura 10.8 se obtiene con las ecuaciones (10.36) y (10.37).
El punto F se obtiene con la ecuación (10.38). Los puntos A, B, D y E se
obtienen con las ecuaciones (10.39).
289
Adimensional a partir de una serie de hidrógrafas reales, correspondientes a
cuencas de muy diversos tamaños y ubicadas en distintos sitios de los
Estados Unidos. Las expresiones matemáticas con las que se puede
construir el Hidrograma Unitario de T-horas de duración, a partir del
Hidrograma Unitario adimensional del SCS, fueron deducidas definiendo
sus tiempos característicos y representándolas como un Hidrograma
Unitario Triangular.
8
1 (10.40)
Tb = 0.375 = 2.67 = 3
Tr = 1.67Tp (10.41)
290
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
291
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).
=5
TC 3 TL (10.43)
292
TL (10.44)
En donde LC es la longitud del canal principal en pies, S C 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)
293
Obtención del Hidrograma Unitario a partir del Hidrograma Unitario
Adimensional del SCS. El Hidrograma Unitario de T-horas de duración, a
partir del modelo del SCS se obtiene de la siguiente manera:
• 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.
294
1.0 1.00 2.2 0.207 5.0 0.000
n−1
295
⎛t⎞ ⎡ ⎛t ⎞⎤
)
U = Up ⎜ ⎜ T ⎟⎟exp⎢(1−n ⎜⎜ Tp −1⎟⎠⎟⎥⎥⎦ para t ≤ t0 (10.48)
⎝ p⎠ ⎢⎣ ⎝
⎛ t0 −t ⎞
U = U0 exp ⎜ ⎟ para t0 < t ≤ t1 (10.49)
⎝K⎠
⎛ t1 − t ⎞
U = U1 exp ⎜ ⎟ para t > t1 (10.50)
⎝ 3K ⎠
296
0.124
⎛ L⎞
K ⎜⎜⎝
W ⎟⎟⎠c
c (10.51)
0.133
⎛ L⎞
Tp ⎜⎜⎝
Wcc ⎟⎟⎠ (10.52)
297
FIGURA 10.9. Hidrograma Unitario Instantáneo de Williams y Hann
⎡ 1 ⎤
t0 = Tp ⎢1+ 1 ⎥ (10.53)
⎢⎣ (n−1) 2
⎥⎦
298
t1 = t0 + 2K (10.54)
⎧ ⎡ 1 1 ⎤ 12 ⎫⎪2 ⎪ 1
⎪ p ⎢⎣4K Tp
∗
Up = B Ac (10.56)
Tp
299
Construcción del Hidrograma Unitario Sintético de Williams y Hann.
Para construir el Hidrograma Unitario Sintético de Williams y Hann, se
siguen los siguientes pasos:
100
B
10
0 2 4 6 8 10 12
n
3
• Se calcula el caudal pico, Up, en pie /s, con la ecuación (10.56).
300
• Se discretizan las abscisas en intervalos ∆T. El intervalo debe ser tal que
la duración, T, de la lluvia a la que se le va a aplicar el Hidrograma
Unitario Sintético sea múltiplo entero de ∆T. T= n∆T.
• Se obtienen las ordenadas del Hidrograma Unitario Sintético, Ut,
reemplazando cada abscisa, t, en las ecuaciones (10.48), (10.49) y (10.50.
3
El Hidrograma Unitario Instantáneo, tiene dimensiones de pie /s /pul
para los caudales y de horas para los tiempos.
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
301
• Tb = 600 min
• W50 = 3.0285 h ≈ 180 min
• W75 = 1.7306 h ≈ 100 min
Parámetros Valor
Area de la cuenca, Ac (km ) 2
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
302
unitario de 0.33 horas de duración. Los resultados obtenidos son los
siguientes:
• S = 2.83
• TL = 1.8 h
• Tp = 1.965 h
• Tb = 9.825 h
• Up = 14.2378 m3/s /mm
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:
303
• 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
16
14
Caudal (m3/s)
12
10
8
6
4
2
0
0 200 400 600 800
Tiempo (min)
304
Y HANN
Caudal (m3/s)
18
16
14
12
10
8
6
4
2
0
0 100 200 300 400 500
Tiempo (min)
305
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
En el método racional, también es necesario conocer la tormenta de
diseño, para hallar el caudal asociado a un período de retorno
cualquiera.
298
TABLA 11.1 Tiempos de concentración para la cuenca del río San
Carlos (Smith, Vélez, 1997)
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
Precipitación efectiva.
299
La precipitación efectiva es la parte de la precipitación que produce
escorrentía superficial directa. Los métodos más populares para
determinarla son el del índice φ y el del SCS, explicados en el capítulo
8.
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
el método debe dar buenos resultados, solo en cuencas pequeñas no
mayores de 50 km2. La intensidad se obtiene de las curvas I-D-F, para
una lluvia con un período de retorno igual al tiempo de concentración,
Tc.
301
QTR = 0.278CTRITcTRA (11.2)
en donde:
QTR : Caudal en m3/s, para un período de retorno TR
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.
302
FIGURA 11.1 Isolíneas de coeficientes de escorrentía para un Tr=50
años (Smith, Vélez, 1997)
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
Planas arenosas cultivadas 0.2
Planas arenosas bosques 0.1 Colinas
arcillosas cultivadas 0.5 Colinas arcillosas bosques
0.4
Colinas arenosas cultivadas 0.4
Colinas arenosas bosques 0.3
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
Empresas Públicas de Medellín (Anuario Hidrometeorológico, 1991,
pag 129). Para un período de retorno de 50 años, la intensidad para 20
minutos es de 300 mm/h.
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
305
(mm) 1.05 0.81 0.73 0.15
Período
Método Q m3/s
Snyder 1.56
SCS 1.91
Williams y Hann 0.28
Racional Probabilístico 1.85
Racional 5.29
306
FIGURA 11.3 Hidrógrafa del SCS para la q. El Limón
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
307
FIGURA 11.4 Hidrógrafa de William y Hann para la q. El Limón
308
extrapolaciones para períodos de retorno mayores, aun con series de
registros mas cortas, pero estiman los intervalos de confianza del valor
de interés (Mutreja, 1986)
309
implica que para cada caso se debe hacer un estudio cuidadoso de todas
las variables influenciandola.
310
estimar y usar distribuciones de probabilidad de valores extremos, para
definir caudales máximos asociados a diferentes períodos de retorno.
utilizando la expresión del factor de frecuencia (Chow, 1953):
311
−2 A0.166L1.716S10.068 (11.5) σˆ =
1.507×10
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.
312
El método fue propuesto inicialmente en 1960 por Dalrimple y Benson.
Se basa en el supuesto de que la distribución de las crecientes en
diferentes sitios de una región homogénea es la misma, escalada por un
parámetro llamado Índice de Crecientes, que refleja el tamaño, la lluvia
y la escorrentía características de cada cuenca. Generalmente, se emplea
la media de los caudales máximos instantáneos como índice de
crecientes. Sí se desea estimar el caudal asociado a una determinada
frecuencia, se halla la relación correspondiente del caudal con la media
de los caudales máximos instantáneos usada como índice. Esa relación
se puede regionalizar para una zona hidrológicamente similar. Estas
curvas son llamadas curvas regionales de frecuencia
313
geomorfológicos como el método de la sección llena explicado en el
capítulo 4.
1,80
1,60
1,40
1,2603 F X (x)
1,20 Q/Q2,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
F X (x)
X =ξ+α[−Ln(Fx(x)] (11.6)
314
En donde F(x) es la probabilidad de no excedencia (función acumulada
de probabilidades) del evento x; y α es el gradex (gradiente de valores
extremos).
315
máximos, en vez de un caudal con un período de retorno de 10 años. La
ecuación quedaría entonces de la siguiente forma (si se emplea la
distribución Gumbel):
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.
316
Como ilustración de la diferencia de resultados que pueden obtenerse
con las diferentes metodologías, se presentan los valores hallados para la
cuenca del río San Carlos hasta, la estación Puente Arkansas (Arbelaez,
Vélez, Smith,1997), donde se tiene una serie de registros de caudal de 16
años de longitud.
317
GUMBEL 699.3 790.5 881.1
GRADEX con
701.8 814.1 926.4
Tr=2.33
GRADEX con
724.7 837.0 949.3
TR =10
INDICE DE
596.7 611.9 619.6
CRECIENTES
Ca
pítulo
12
El AGUA SUBTERRÁNEA
318
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.
319
En Urabá, la explotación del banano (segundo renglón agrícola de
exportación de Colombia) depende en gran parte de los recursos de agua
subterránea. En 1993, el volumen diario de explotación de agua
subterránea en la región era de unos 58000 m³, de los cuales el 65% se
usaba en la agricultura y el resto en la industria y abastecimiento público
INGEOMINAS, 1993).
320
A manera de resumen, puede afirmarse que las aguas subterráneas
constituyen un recurso natural de gran importancia para el hombre, en
razón, principalmente, de:
321
superficie imaginaria que coincide con el nivel hidrostático del agua
en el acuífero, Figura 12.1.
322
dejan entre ellos vacíos de muy pequeñas dimensiones, llamados
poros, que pueden ser ocupados por líquidos o gases.
- En las no consolidadas.
- En las compactas fisuradas.
323
• El análisis estructural para las rocas compactas y fisuradas, realizado
tanto en el campo como en el laboratorio.
- Continuidad o discontinuidad
- Isotropía o anisotropía
- Homogeneidad o heterogeneidad
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.
324
almacenar ni transmitir el agua, como ciertas rocas metamórficas no
fracturadas.
- Agua de retención
- Agua gravítica
325
Esta clasificación se basa en el estudio de las distintas fuerzas a que se
pueden ver sometidas las moléculas de agua presentes en los medios
porosos. Desde el punto de vista hidrogeológico, esta distinción es
fundamental, ya que no toda el agua que hay en los vacíos de una
formación determinada es utilizable en forma efectiva, situación que
está ligada al concepto definido anteriormente de porosidad eficaz. La
Tabla 12.1 es una clasificación de los diferentes tipos de aguas
subterráneas que existen, según los medios y la energía necesaria para su
extracción.
326
geológica cualquiera. La parte más importante desde el punto de
vista hidrogeológico es la porosidad eficaz (ne), que mide la
cantidad de agua que el acuífero puede liberar por acción de la
gravedad, es decir, el agua útil o realmente explotable.
T=Kb (12.1)
327
4) Coeficiente de almacenamiento: en los acuíferos libres, el agua
que se extrae a través de un pozo proviene esencialmente del
agua liberada por efecto de la gravedad. En los acuíferos
cautivos, que se encuentran completamente saturados y bajo
otras condiciones de presión, el agua es liberada debido a
fenómenos de elasticidad que operan tanto sobre el material
sólido o armazón del acuífero como sobre el agua misma; esto se
estudiará mas adelante.
328
FIGURA 12.2 Experimento de Darcy
329
Si en lugar de considerar el tubo vertical se le da una cierta inclinación,
figura 12.3 se tiene:
∆h
i= (12.4)
l
∆h (12.5
Q =− KA =−KAi
)
l
330
La 12.5 puede reescribirse:
QV= (12.6
= Ki )
A
- La permeabilidad, K.
- El gradiente hidráulico, i.
- La velocidad del flujo, V.
Ejemplo 12.1
331
cm/s. Si la presión es siempre la misma en la columna, cual es la
conductividad hidráulica ?.
Solución:
−h2 ∆h
Donde el gradiente hidráulico es: i = h1 =
L L
Ahora;
Lsen30o
i= ⇒ i = sen30°
L
0.005
k= o ⇒ k = 0.01 cm/s sen30
332
12.2.1 Estudio de la permeabilidad
Si V = Ki, en donde:
V: velocidad de flujo,
i: gradiente hidráulico,
γ
K =K 0 (12.7) µ
En donde:
333
ν= (12.8) ρ
en donde ρ es la densidad del fluido, y como γ = ρ g, entonces:
γ ρg g
= = µ νρ ν
g
K =K 0 (12.9) ν
K 0 =Cd102 (12.10)
334
Arena media, bien gradada 80 - 120
Grava arenosa, pobremente gradada 80 - 120
Grava arenosa limpia, bien gradada 120 – 150
335
Se puede utilizar ya sea permeámetro de cabeza constante o de
cabeza variable.
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
336
Q =KA ∆h K= Ql l A∆h
h
Q =KA
l
337
FIGURA 12.4 a) Permeamétro cabeza constante. b)Permeamétro
cabeza variable
dh dh AKdt
Q =−a = dt
h a L
h AK
ln =− (t − t0 ) (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.
338
- La muestra puede ser perturbada y no representa, por consiguiente, la
situación en el terreno, por ejemplo, en cuanto hace relación al grado
de compactación.
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
339
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 = Vrne K = Vr ne
i
340
FIGURA 12.5 Esquema típico para un ensayo de bombeo.
h1 = pA +ZA h2 = pB +ZB γ γ
341
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
i =− lim =− (12.16)
∆h→0 l ∂l
342
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.
343
i= −∆h∆h = hA −hB
l
hA = Z1 −h1hB = Z2 −h2
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.
344
isopiezas se puede calcular el gradiente hidráulico para grandes
distancias.
i= 100−90 = 0.004
2500
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)
..
345
denomina también velocidad aparente o velocidad de flujo, para
diferenciarla de la velocidad real, que viene a ser la velocidad de
desplazamiento de las moléculas de agua a través de los poros y que
puede determinarse mediante el uso de trazadores. A esta velocidad real
de flujo se le denomina también velocidad de filtración o velocidad
eficaz o efectiva.
S = Ane (12.17)
346
Ahora, si se considera el flujo a través de una sección A, se tiene que por
la ley de Darcy:
Q = AKi (12.18)
347
REFERENCIAS BIBLIOGRÁFICAS
370
CARDONA, A., LONDOÑO, G.A. 1991. Calibración de los
Parámetros del Modelo de Williams y Hann para Cuencas Antioqueñas.
Tesis (Ingeniero Civil) Universidad Nacional. Medellín. 150 pp
371
GARCON, R.. 1994 The Gradex Method: for a statistically robust
evaluation of extreme value floods based on local characteristics.
Electricite de France- Generation and Transmission Division- Water
Resources Department France.Kluewer Academic
Publisher
Neteherlands,vol 1, pp 285-298.
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
INGEOMINAS. Evaluación del Agua Subterránea en la Región de
Urabá, Departamento de Antioquia. 1993.
373
NOAA. 1999. Página de la Red Global.www.noaa.gov
374
Editores técnicos. 515 pp
375
UNESCO. 1981. Métodos de Cálculo de balance Hídrico. Guía
Internacional de Investigación y Métodos. España: Instituto de
Hidrología de España,. 193 pp
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
Congreso Latinoamericano de Hidraulica. Santiago de Chile. Noviembre
377