Hidrologia - David Cedeño

F
TNI\'ERSIDAD TECNOLOGICA DE PANAMA

FACIJLTAD DE IIIGENTERIA CIVIL
Departamento de llidrulica, Sanitaria y

Ciencias Ambientales
Apuntes de
HIDROLOGIA
Preparados por
Ing. Dad Cedeo B.
Verano
L997
CONTENIDO
PRINCIPIOS HIDROLOGICOS
1
Introduccin.
a
1 1
) )
Evolucin de la Hidrologa
| -2
1-
Ciclo Hidrolgico
a
Balance Hdrico
l5
Precipitacin.
2t
.
a) HumedadAtmosfrica . . . b ) Cambios de Fase c) Cantidad de Agua Precipitable d ) Causas y Mecanismos de Formacin de la

Precipitacin
22
26
27
32
) f) g) h)
e 1-
Anlisis de Datos de Precipitacin

Estimacin de Datos Faltantes
34
38
Anlisis de Doble Masa
4l
46 52 54
55
Precipitacin Promedio sobre una Regin
Evaporaciny Transpiracin
) b) c) d) e) f)
a
Mtodo del Balance Hdrico Mtodo de Transferencia de Masa Mtodo de Balance Energtico
Tanque Evapormetro Mtodos Combinados
56 59 60
66
Evapotranspiracin.
1-
Infiltracin
) b)
a
Rata de
Infiltracin Infiltracin . Aforos
69
73
Otros Mtodos para calcular la
1-6
EscorrentaSuperficial
a
78
80
Medicin del Caudal:
2.
ANALISIS DE PRECIPITACION - ESCORRENTIA
2.1RelacinentrePrecipitacinyEscorrenta
87 Racional 90 2 - 2 Anlisis de Hidrogramas 94 a) Componentes del Hidrograma 91 b ) Separacin del Flujo Base y Recesin. 99 Precipitacin Neta y el Hidrograma. c) . 102 d ) Mtodo de Tiempo - Area . . 108 2 - 3 Teoa del Hidrograma Unitario . 109 Derivacin de Hidrogramas Unitarios. a) . 119 b ) Mtodo de la Curva S . Mtodos Matriciales para Desarrollar c) . 128 Hidrogramas Unitarios 2-4 Desarrollo de Hidrogramas Unitarios Sintticos . 137 . l4l Mtodo de Snyder a) b ) Mtodo SCS (Hidrograma Unitario Triangular) . 149 . 158 2-5 Aplicaciones de Hidrogramas Unitarios. . 159 Convolucin de Hidrogramas Unitarios a) a
Mtodo
ANALISIS DE FRECUENCIAS
3-1 Introduccin. .166 a) Variables Aleatorias . 166 b ) Presentacin de los Datos . 168 c) Conceptos de Probabilidad. . 174 3 -2 Yariables Aleatorias y Distribuciones de Probabilidad . 177 a) Momento de una Distribucin . 182 b ) Estimacin de los Momentos a partir de los Datos . 188 c) Ajuste de una Distribucin a los Datos . 195 3 - 3 Perodo de Retorno Intervalo de Recurrencia . I91 a) Clasificacin de los Datos . 200 3 - 4 Modelos Probabilsticos Comunes . 205 a) DistribucinBinomial . .208 b) RiesgoyConfiabilidad .210 Distribucin Exponencial c) . 213 d ) Distribucin Normal . 217 e) Distribucin Log Normal . 223 f ) Distribucin Gamma (2 Parmetros y Pearson Tipo 3) . 231 Distribucin Log Pearson Tipo 3 . . 240 C) h ) Distribucin Gumbel Valor Extremo Tipo I . 244
REFERENCIAS
250
lll
PRINCIPIOS HIDROLOGICOS
INTRODUCCION
Hidrologa es una ciencia multidiciplinaria que estudia la ocurrencia,
circulacin
distribucin del agua sobre la
tierra. El
dominio de la
hidrologa incluye los procesos fsicos, qumicos y las reacciones biolgicas
del agua en los ambientes naturales y aquellos construidos por el hombre.

Debido a la compleja naixalezadel ciclo hidrolgico
y- sus
relaciones con los
patrones climticos, tipos de suelo, topografa y otros factores geolgicos, es
dificil establecer fronteras entre la hidrologa y otras ciencias de la tierra tales

como: Meteorologa, Geologa, Ecologa y Oceanografa entre otras.
La hidrologa es de fundamental importancia para los Ingenieros Civiles
Ambientales, Hidrlogos
y otros cientficos
relacionados con la tierra;
debido a las implicaciones en el medio ambiente de las inundaciones y

sequas, abastecimiento de agua, consideraciones sobre la calidad del agua,
drena-jes
y control de inundaciones.
EVOLUCION DE LA HIDROLOGIA
La
historia inicial de
la
hidrologa incluye las prcticas
de
administracin
de agua en las civilizaciones Egipcias y

1o
Sumerias
(Mesopotamia) establecidas en el Medio Oriente a
largo de los ros Nilo,
Tigris y Efrates y en la China a lo largo del Ro Amarillo. Las excavaciones
Hidrologa / David
Cedeo
arqueolgicas en estos lugares muestran evidencias de estructuras hidrulicas
que fueron construidas para irrigacin y otros proyectos para el control del
agua. Por ejemplo en el ao 4000 A.C. se construy un canal para

transportar agua dulce entre El Cairo y Suez (Alejandra).
Los filsofos griegos fueron los primeros investigadores de
la
hidrologa. Aristteles propuso la teora de que la humedad del aire se

converta en agua dentro de las montaas y que sta era la fuente de todas las
corrientes de agua
Homero sugiri la idea de que exista un mar
subterrneo que era el origen del agua superficial.

Las tcnicas de medicin de caudales se intentaron por primera vez en
los sistemas de acueductos Romanos (97 D.C.) basados en la

transversal del
seccin
flujo.
Estos teoremas continuaron hasta el Renacimiento
cuando Leonardo Da Vinci descubri una relacin adecuada entre rea,

velocidad y caudal (Q
: V.A).
La primera medicin registrada de precipitacin y flujo superficial

fueron hechas en el siglo XVII por Perrault, quien compar la cantidad de lluvia medida y el flujo estimado del ro Sena, mostrando que ambas variables
estaban relacionadas; estas observaciones fueron publicadas en 1694.
El astrnomo ingls Halley (1636 - 1742) utlliz una bandeja pequea
para estimar la evaporacin del Mar Mediterrneo

equivalente
concluy, que era tambin midi la
al flujo de los ros tributarios. Mariote
velocidad del flujo en el ro Sena.

Estos comienzos modestos de la ciencia de la hidrologa establecieron
las bases para los avances de esta ciencia en el siglo XVIII, incluyendo el
Hidrologa / Davd Cedeo

Teorema de Bernoulli, el tubo de Pitot y la frmula de Chezy (1769), que
constituyen la base de la hidrulica y las mediciones de fluidos.
Durante el siglo XIX, ocurrieron avances significativos en hidrologa subterrnea: La ley de Darcy para flujos en medios porosos, la formula de
pozos de
Dupuit-Thiem y tambin
se desarrollo la ecuacin de
flujo capilar
de Hagen-Poiseville.
En hidrologa superficial se desarrollaron muchas frmulas

de las corrientes de agua. Humprey y
instrumentos de medicin que permitieron el inicio de la medicin sistemtica
Abbot (1861) efectuaron mediciones
de caudales en el ro Mississipp, y el U.S.G.S. (United States Geological Survey) estableci un programa para la medicin y registro de caudales del Mississippi en 1888. Adems, la formula de Manning se introdujo en 1889
y el medidor
de corriente fue inventado por Price en 1885. Durante este
perodo el gobierno de los Estados Unidos fund una serie de agencias

hidrolgicas incluyendo el U.S. Army Corps of Engineers (1802), U.S.
Geological Survey (1879) y el Weather Bureau (1891).
El intervalo de tiempo comprendido entre 1900 hasta 1930 se denomin
"Perodo de Empiricismo" (Chow, 1964), debido a la gran cantidad de

frmulas empricas que se desarrollaron; muchas de las cuales resultaron ser
incorrectas. Las agencias del gobierno incrementaron sus esfuerzos en la

investigacin hidrolgica
tcnicas para
se organizaron cierto nmero de sociedades
el avance de esta ciencia. Por ejemplo, el U.S. Bureau of
Reclamation (1902), el Forest Service (1906) V el U.S. Army Engineers
Waterways Experiment Station (1928) fueron organizadas durante este
Hidrologa / David Cedeo
perodo. Adems, The International Association of Scientific Hidrology

(1922) y The Hidrologic Section del American Geophysical Union (A.G.U.) empez antes de 1930.
El
perodo
de 1930 a 1950 se le
denomin "Perodo de
Racionalizacin", el cual produjo un avance significativo en el campo de la hidrologa; ya que las agencias de varios gobierno desarrollaron programas
de investigacin hidrolgica. Entre los avances significativos del perodo tenemos: Sherman (1932) estableci el concepto de hidrograma unitario,
Norton (1933) desarroll la teora
de infiltracin, Theis (1935) formul la
ecuacin de flujo inestable en hidrulica de pozos y
Gumbel (1950) propuso
el uso de la distribucin de valores extremos para el anlisis de frecuencias

de datos hidrolgicos, estabieciendo las bases de la hidrologa estocstica.
En este perodo el U.S. Army Corps of Engineers, el U.S. Weather

Bureau (ahora National Weather Service), el U.S. Department of Agriculture
(U.S.D.A.) a travs del Soil Conservation Service y el U.S. Geological

Survey (U.S.G.S.) aportaron contribuciones significativas
hidrolgica
a la teora
y se continu con el desarrollo de redes de medicin para registrar la precipitacin, evaporacin y escorrenta. Estas agencias
efectuaron estudios vitales
y proporcionaron fondos para la investigacin
privada
universitaria en el rea de hidrologa. Las grandes presas,
embalses, proyectos de control de inundaciones, etc. son el resultado directo de los avances en los campos de Mecnica de Fluidos, Sistemas Hidrolgicos,
Hidrologa Estadstica, Anlisis de Evaporacin, Trnsito de Avenidas

Investigacin de Operaciones.
A partir de 1950,
se conoce como
"Perodo de Teorizacin", ya que
la introduccin de computadoras digitales en hidrologa durante 1960
1970
permiti la simulacin de problemas complejos de sistemas hidrolgicos.
El primer modelo
hidrolgico comprensivo fue desarrollado por

se denomin
Crawford y Linsley (1966) en la Universidad de Stanford y
Stanford Watershed Model (S.W.M.). Este modelo puede simular los

procesos principales del ciclo hidrolgico: precipitacin (P), evaporacin (E),
transpiracin
(T), infiltracin (F),
escorrenta superficial (R)
flujo
subterrneo (G). Otro modelo que ha alterado significativamente el curso de la hidrologa

moderna es el programa
IIEC-I
(1973) desarrollado por el U.S. Army Corps
of Engineers, Hidrologic Engineering Center, Davis, California. Este modelo simula inundaciones a partir de datos de precipitacin utilizando hidrogramas unitarios y funciones elementales de prdidas. Otro modelo que le acompaa es el HEC-2 (1976), desarrollado tambin por el Hidrologic Engineering Center, el cual efecta clculos de los perfiles de la superficie del agua para
una geometra conocida del canal y caudales mximos, los cuales se pueden obtener utilizando
HEC-I.
El Storm Water Management Model (S.W.M.M.) fue desarrollado

por el U.S. Enviromental Protection Agency (E.P.A.) durante 1981 a 1988
y es el modelo ms comprensivo para el anlisis de escorrenta urbana

sistemas de alcantarillados.
en
El modelo ILLUDAS (Illinois Urban Drainage Area Simulator), desarrollado por Terstriep y Stalt (1974) est basado en un modelo del
Hidrologa / David
Cedeo
y
British Road Research Laboratory
utrliza un procedimiento simple de
precipitacin-escorrenta para la simulacin de tormentas

drenajes adecuados.
y el diseo de
Estos modelos se han convertido en

frecuentemente
un
instrumento utilizado
por los investigadores e ingenieros hidrlogos y la lista
representa algunos de los programas de computadora ms poderosos de la hidrologa moderna. El desarrollo de esta herramienta en los ltimos 20 aos
ha ayudado directamente en la coleccin de datos hidrolgicos al permitir la
calibracin del modelo contra los datos observados. Por consiguiente, se ha

avanzado mucho con estos procesos en la comprensin del comportamiento de los sistemas hidrolgicos.
Los modelos hidrolgicos para computadoras desarrollados inicialmente
en 1960
1970 han sido aplicados satisfactoriamente a otras reas que
anteriormente no se estudiaban
que estaban definidas empricamente
solarnente. Por ejemplo, hidrologa urbana, hidrologa de cuencas y planicies

de inundacin, diseo de drenajes, diseo y operacin de embalses, anlisis
de la frecuencia de inundaciones
sequas; adems,
la administracin y
planificacin de la cuenca de los ros han resultado beneficiados por la

aplicacin de los modelos de computadora.
Las limitaciones de los modelos de simulacin incluyen el peligro de
creer que un modelo producir resultados adecuados para todas
las
situaciones. Una confianza excesiva en los programas de computadora en la dcada de 1970 condujo a un tratamiento ms cuidadoso de los modelos hidrolgicos en la dcada de 1980, produciendo un regreso a las aplicaciones
de los modelos que no excedieran la disponibilidad de datos satisfactorios de
entrada. Sin embargo, la simulacin de modelos en hidrologa, al aplicarse

correctamente, producir la aproximacin ms lgica a la comprensin de los
procesos complejos que ocurren en el ciclo hidrolgico y por lo tanto, nos

encontramos en una nueva era en
1a
ciencia de la hidrologa.
EL CICLO HIDROLOGICO
Las componentes bsicas del ciclo hidrolg-ico son las siguientes:
precipitacin (P), evaporacin (E), transpiracin
(T), infiltracin
(F),
escorrenta superficial (R) y flujo subterrneo (G). El ciclo del agua es un

proceso contnuo en el cual el agua se evapora de los ocanos, lagos, ros y
otras fuentes, se mueve en la atmsfera formando masas de aire hmedo y

luego se produce la precipitacin cuando existen condiciones adecuadas. La
lluvia que cae en la superficie de la tierra se dispersa a travs de muchos

medios; pero una porcin es retenida en el suelo cerca del sitio donde cay
y retorna a la atrnsfera por evaporacin, es decir, la conversin del lquido

a vapor de agua
y tambin por transpiracin, la cual consiste en la prdida de
agua a travs de las plantas; esta prdida de agua combinada se denomina:
evapotranspiracin(ET: E + T).
Otra porcin del agua se convierte en flujo superficial o escorrenta
directa (R), la cual abastece las corrientes y ros. Finalmente, la parte restante
del agua entra al suelo como infiltracin (F), la cual puede convertirse en flujo subsuperficial y aparecer posteriormente en canales percolarse hacia las profundidades para recargar el flujo del agua subterrnea.
El
agua
superficial
subterrnea
se mueve hacia elevaciones inferiores y
eventualmente puede descargar en cantidades de agua superficial
el ocano. Sin
embargo, grandes
porciones de agua subterrnea pueden
regresar a la atmsfera a travs del proceso de evapotranspiracin.
aqua
atmosferica
fI
'"fS#,ti,r
f_
I
superlicial {lagos)
aqua
subtrranea
Figura No1: Representacin Esquemtica del Ciclo Hidrolgico.
La estimacin de la cantidad total del agua en la tierra y en los

diferentes procesos del ciclo hidrolgico ha sido un tema de investigacin
cientfica desde hace muchos aos. Sin embargo, la cantidad de datos

escasa, particularmente sobre los ocanos
es
y tambin la cantidad de agua en
los diferentes componentes del ciclo hidrolgico global todava no se conoce
con exactitud.
La Tabla 1: "Cantidades Mundiales Estimadas de Agua" muestra los

volmenes estimados de agua en las diversas fcmas sobre
la
tierra.
Alrededor del 96.5% de toda el agua sobre Ia tierra est en los ocanos. Si
la tierra fuera una esfera uniforme, esta cantidad sera suficiente para cubrir
el planeta hasta una profundidad de aproximadamente
2.6Kn
(1.6 millas).
Del resto del agua, 7.7% es hielo polar, 1.69% es agua subterrnea y
solamente
0.ll%
se encuentra en la superficie
y la atmsfera. El
sistema
atrnosfrico del agua contiene solamente 12,900 Km3 de agua; es decir, menos
de una parte en 100,000 del total de agua de la tierra y es la componente
principal que impulsa la hidrologa del agua superficial.

Del volumen total de agua dulce en la tiera, aproximadamente 2/3 es
hielo polar
y la mayor parte del resto es agua subterrnea hasta una
profundidad de 200 a 600 metros. Debajo de esta profundidad la mayora del

agua subterrnea es salina. Solamente 0.006% del agua dulce est contenida
en los
ros. El agua biolgica,
retenida en los tejidos de las plantas y
animales, contiene aproximadamente 0.003% de toda
el
agua dulce,
equivalente a la mitad del volumen de agua contenida en los ros.
10
TABLA
N'1:
Cantidades Mundiafes Estimadas de Agua

? totaf
96.5
Descripcin
ocanos
Agua
Area (106 Km? ) 361.3
Volumen
(Km3
de agua
1,338,000,000
?agua dul ce
Subterrnea:
Dul ce
Safina
Humedad del Suel o Capas
134.8 134.8
82.O
16 .0
10,530, 000 L2,870,O00 l-6,500

24 ,023 , 500
0.76
0 _ 93
30.1
0
0.0012
.05
L.7
0.02s
68.6
Pofares
Glaciares y
N
0.30
340,600
1.0
aeve
Lagos:
Dul ce
S
al ina
os
L.2 0.8
9r-,000 85,400
0.007 0.006 0.0008 0.0002 0.0001 0.001

100
.26 .03
Pant ano s
Ra
2.7
148.8 510.0 5l-0.0 510.0 148.8
rr,470
2,120 a,L20
t2 , 900
1, 38s, 984 , 6L0
35 , 029
0.006 0.003
0
Aguas B i of gi cas
Agua
Atmo
fri ca Total de
s
.04
Aguas
Agua Dufce
,2r0
2.5
100
11
TABLA
N"2: Bafance Hdrico Global Anuaf.

LJn1clacl
Km2
Descripcin
Area
Prec
OCEANO
TIERRA
361,300,000
4s8, 000
148, 800, 000

r_19, 000
72 , OO0
ipi
tac in
Km"/ Km' /
ano ano
Evaporac
Es
in
50s,000
corrent.
ia
Mar
Km'/ ano 44 ,7 00
hacia ef
Fluj
o
Subt errneo
m'/ ano
2,200
47 ,
Escorrentia Totaf
Km-
/ ano
OOO
FUENTE: World Water Balance and Iat.er Resources of the Eart.h, UNESCO, 1978.
A pesar de que el contenido de agua en la superficie y en la atrnsfera
es relativarnente pequea en cualquier momento, cantidades inmensas de agua
pasan anuahnente a travs de estos sistemas. La Tabla 2 muestra el balance
hdrico global anual y la Figura 2 (Balance Hdrico Promedio Gtobal Anual)
muestra los componentes principales del ciclo hidrolgico en unidades
relativas a una precipitacin anual sobre la tierra de 100 unidades de

volumen.
Se puede observar que la evaporacin sobre la
tierra consume 61%
de la precipitaciny el resto39%, constituye la escorrentia hacia los ocanos,
principalmente en forma superficial.
La
evaporacin de
los
ocanos
1-2
contribuye con el 90% de la humedad atmosfrica. El antisis del flujo y

almacenamiento del agua en el balance global anual proporciona algunos conocimientos elementales sobre la dinmica del ciclo hidrolgico.
l6t \/
r,\
-;:i+
'* I
lr
,#?
*i
I
(?:
lrY
Figura No2: Balance Hdrico Promedio Global Anual.
El ciclo del agua es muy complejo, pero bajo ciertas condiciones bien
definidas,
la
repuesta de
la
cuenca
a la
precipitacin, infiltracin y
evaporacin se pueden calcular si se establecen suposiciones simples. Por

ejemplo, si la rata de precipitacin sobre una cuenca es menor que la rata de
13
infiltracin y si existe un amplio almacenamiento en la humedad del suelo,

entonces la escorrentia directa en la superficie
canales de drenaje ser
y el flujo resultante en los
cero. Por el contrario, si la precipitacin precedente
ha llenado la capacidad de almacenamiento de la humedad del suelo y la

intensidad de la precipitacin es mucho mayor que la rata de infiltracin y
evaporacin, entonces el volumen de escorrentia superficial ser igual al

volumen de precipitacin. En la mayora de los casos, desafortunadamente,
las condiciones existentes quedaran localizadas eritre estos dos lmites y

debemos medir cuidadosamente o calcular ms de una componente del ciclo
hidrolgico para predecir la respuesta de la cuenca.
El ingeniero hidrlogo debe ser capaz de calcular o estimar las diversas

componentes del ciclo hidrolgico para disear adecuadamente proyectos de
recursos hidrulicos. Muchos de estos proyectos hidrulicos deben ser

diseados para proteccin contra los daos producidos por eventos extremos de inundaciones y sequas y sern operados generalmente tomando en cuenta
estos eventos
crticos. Algunos de los temas tpicos relacionados con la
ingeniera hidrolgica incluyen los siguientes:
l. 2.
Flujos mximos de inundacin esperadas en los vertederos,

alcantarillas y puentes de las carreteras.
las
Capacidad de los embalses requerida para asegurar una cantidad

adecuada de agua para irrigacin ciudades.
y abastecimiento de agua para las
3.
Efectos de embalses, muros de contencin

control de inundaciones en una corriente.
otras estructuras de

4.
L4
Efectos de desarrollo urbano en la capacidad futura de un sistema de

drenaje y los flujos asociados con las inundaciones.
5.
Determinacin de los niveles probables de inundacin para mejorar la proteccin que ofrecen los proyectos construidos por el hombre contra
las inundaciones para promover el establecimiento de zonas con

riesgo de inundaciones.
Ejemplo
1: Estimar el tiempo de residencia de la humedad

global.
atmosfrica
Solucin: El tiempo
de residencia
T,
es la duracin promedio para que una
molcula de agua pase a travs de un subsistema del ciclo hidrolgico.
T=
De la Tabla 1: De la Tabla 2:
s
O
volumen de agua almacenada
rata de flujo
12,900 Km3
s--
a: 458,000 + 119,000 :
12,900 Km3
577
577,000 Km3/ao
T=
= 0.022 aos = 8.2 das
,000 Km3lao
Observacin: Debido al corto perodo de residencia de la humedad en la
atmsfera, es difcil predecir
el
estado del tiempo con varios das de

puede
anticipacin. Este valor del tiempo de residencia es promedio y

mostrar variaciones espaciales considerables.
l_5
BALANCE HIDRICO
Para cualquier sistema hidrolgico, se puede desarrollar un balance
hdrico para tomar en cuenta las trayectorias del flujo y el almacenamiento de
agua. El sistema ms simple es una superficie impermeable

confinada en todos sus bordes
inclinada,
con una sola salida. Un lote pequeo
pavimentado de estacionamiento en un rea urbana satisface este modelo.
Entrada =
Salda = O
Figura N'3: Sistema Hidrolgico Simple

La ecuacin de continuidad (flujo no permanente inestable) para
cualquier sistema hidrolgico es:
r-o-ds
dt
Hdrologa / David Cedeo

donde:
I6
flujo de entrada (vol/tiempo)
o
ds/dt
flujo de salida (vol/tiempo) cambio en
el volumen de almacenamiento por
unidad de tiempo.
Al
aplicar este modelo al lote de estacionamiento (figura 3), la

a
precipitacin se acumula sobre la superficie y eventualmente se descarga
travs de la salida. Si despreciamos la evaporacin durante todo el perodo,
eventualmente toda la precipitacin se convertir en
flujo de salida, pero
estar algo retrasada con respecto al tiempo. La diferencia entre el caudal de entrada acumulado
y el caudal de salida en cualquier momento representa el
cambio en el almacenamiento, el cual se descarga a travs de la salida,

despus que halla finalizado la lluvia.
Figura N'4:
El mismo concepto
Balance hdrico de una Cuenca.
se puede aplicar a una cuenca pequea o grande,
con la diferencia que es ms difcil el anlisis, pues algunos de los trminos
1-'7
correspondientes a los perodos en el balance hdrico pueden ser desconocidos
o difciles de evaluar. Una Cuenca se define como una superficie de tierra que es drenada
a
travs de una salida nica y que esta separada de las otras cuencas por una
divisoria de aguas. Para un perodo de tiempo, el modelo matemtico

conceptual para la cuenca es el balance hdrico (figura 4) el cual se puede
expresar (en unidades de profundidad: cm plg) como:
P-R-G-E-T:AS donde:
P
precipitacin
escorrentia superficial
R G
flujo subterrneo
evaporacin
T
AS
transpiracin
cambio en el almacenamiento
Tambin podemos definir un coeficiente de escorrentia (R / P) como la
razrentre la escorrenta y la precipitacin. Observe que la infiltracin F es

una prdida en la superficie del sistema y una ganancia para el
flujo de agua
general.
subterrnea, por consiguiente se cancela en
el balance hdrico
Adems, las unidades de profundidad (cm plg) representan un volumen de agua cuando se multiplican por el rea de la cuenca.

Ejemplo
18
2:
Balance Hdrico/Unidades de Conversin
Para un mes dado, un lago con una superficie de 300 acres recibe un
caudal de entrada de 15 p3ls y descarga 13 p3ls; adems, el almacenamiento
total se incrementa en 16 acres-pie durante este perodo. Un medidor cerca
del lago registro un total de 1.3 plg de precipitacin durante dicho mes.
Asumiendo que la infiltracin es insignificante, determinar la prdida de agua
por evaporacin sobre el lago (en pulgadas).
I
P = frecipitacion
= evporann
Solucin:
La ecuacin de balance hdrico para la evaporacin del lago

unidades de profundidad) es:
(en
E:I-O+P-F-aS

Donde:
I9
, I:
(ls
pies3lseg) (12 plglpie) (86,400 seglda) (30 das) (300 acres) (43,560 pies2lacre)
35.70 plg
(13 pies3/seg) (12
..\
O
P
plglpie) (36,400 seglda) (30 das)
(300 acres) (43,560 pies2lacre)
: :
30.94 plg 1.3 plg (16 acrespie) (12 plglpie)

300 acres
, -
0.64 ptg
aS
0.64 plg
E:I-O+P-aS E : 35.70 -30.94 + 1.3 - 0.64 :

Ejemplo
5.42pls
3:
Balance Hdrico en una cuenca.
En un ao dado, una cuenca con un rea de 250,000 hectreas recibi
130 cm de precipitacin. El caudal promedio medido en un ro que drena
la cuenca fue de R
30 m3/s. Estimar la cantidad de agua que se pierde
debido a los efectos combinados de evaporacin, transpiracin e infiltracin
hacia el agua subterrnea. Adems, calcular el coeficiente de escorrenta.
20
Asumir que los niveles de agua al inicio (t

iguales; por
1o
0) y al final
(t
0.
tanto el cambio de almacenamiento es AS
1 ao) son
Solucin:
La ecuacin de balance hdrico para la superficie de la cuenca es:
ET+F:P-R-aS
Donde:
* _ (30 m3lsei (86,400 segldla) (365 dlais = 0.3784 m

(2sq000 ha) (1o,ooo m2lha)
37.84 cm
ET+F
ET +
P-R-AS
130
- 37.84
92.16 cm
El coeficiente de escorrenta
es:
R P
37.84 cm
t30
cm
0.29
2a
PRECIPITACION
La precipitacin es la cantidad primaria de entrada del agua en el ciclo
hidrolgico superficial, ya sea en forma de lluvia, nieve
grarizo; y
generalmente se deriva de la humedad atmosfrica. Las masas de aire

hmedo deben estar sometidas a un ascenso, con el enfriamiento resultante,
condensacin y crecimiento de las gotas de agua antes de que la precipitacin
ocurra. La precipitacin
se clasiflca a menudo en tres tipos de acuerdo con
las condiciones que generaron el movimiento verticll de las masas de aire

cargadas de humedad:
A.
Convectiva: Debido al calentamiento intenso del aire al nivel del suelo
el mar, se produce una espansin y el
ascenso vertical del aire
hmedo. Este tipo de precipitacin es caracterstico del trpico.
B.
Ciclnica: Esta asociada con el movimiento de grandes sistemas de

masas de aire, de regiones de alta presin a regiones de baja presin.
(Como en el caso de frentes clidos y fros). Esta diferencia de presin

es creada por el calentamiento desigual de la superficie de la tierra.
C.
Orogrfica: Producida por el ascenso mecnico de las masas de aire

hmedo sobre las montaas.
Hidrologa I David Cedeo
OBOBAFIC,T
[t
YECflfA
FROIITAL
-r-l
I
HIIiTDO
"-.'ff.,'1' l'+
T
TFE F
Figura N"5: Tpos de Precipitacin
HT]MEDAD ATMOSFERICA
La humedad atmosfrica es la fuente requerida para la precipitacin y
se deriva de la evaporacin
y transpiracin.
Mediciones comunes relacionadas con la humedad atmosfrica, o

simplemente humedad, incluyen presin de vapor, humedad especfica,
humedad relativa y temperatura para la formacin de
roco. Bajo condiciones
hmedas, se puede asumir que el vapor de agua satisface las leyes de los
gases ideales, 1o cual permite una derivacin de relaciones sencillas entre
presin, densidad y temperatura.
)a
La Presin Parcial es la presin que actuar sobre la superficie de un

recipiente por un gas particular en una mezcla de gases. La presin parcial
producida por el vapor de agua se denomina Presin de Vapor (e)
se puede
obtener a partir de la ley de Dalton y la ley de los gases ideales de la siguiente manera:
e=
P,RT
0.622
donde:
presin de vapor en mb
densidad de vapor o humedad absoluta en constante de los gases para aire seco
2.87 xl03 (mb. cm3)/(gr. oK)
p'
R R
grlcmr
temperatura absoluta en oK
El factor 0.622 surge de la razt entre el peso molecular del agua
(= 18)
al
peso molecular del aire ("29). Cerca de la superficie de la tierra
la
presin del vapor de agua es
l% al2% delapresin atmosfrica total, donde
la presin atmosfrica promedio es 1,013.2 mb al nivel del mar.
(Nota:
I mb :
100 Pa).
que
La hesin de Saturacin de Vapor (e,) es la presin parcial

ms evaporacin) y es una funcin de la temperatura.
ejerce el vapor de agua cuando el aire esta completamente saturado (no ocurre
24
la
La Humedad Relativa (H) es aproximadamente la raz6n entre

condiciones e igual temperatura. Esta se puede definir como:
presin de vapor de agua a la presin de saturacin de vapor bajo las mismas
H=100e
"
Por consiguiente, 50% de humedad relativa signiflca que la

saturadas a esa temperatura.
atmsfera
contiene 50% de la humedad mxima que podra retener bajo condiciones
La Humedad Especlica (q) es la masa del vapor de agua contenida en

una unidad de masa de aire hmedo y es igual a:
p.,
q=
p-
Observe que la humedad especfica q es adimensional; por 10 tanto, la

densidad del vapor pu
mismas unidades.
Ia densidad del aire hmedo
p. la
deben tener las
Utilizando la ley de Dalton
asumiendo que
atmsfera esta
compuesta de solamente aire y vapor de agua, tenemos:
p^
(P - e) + 0.622
RT
o = P (r-03784) RT\ P)

La ecuacin anterior nos muestra que el aire hmedo es ms liviano que
el aire seco a la misma presin y temperatura por
pu
1o
tanto:
p^
donde:
q e P
Pm
0.622 e - 0.378
humedad especfica (grlgr)
presin de vapor (mb) presin atmosfrica total (rnb)

densidad de la mezcla de aire seco y
vapor de agua (grlcm3) p,

densidad de vapor (grlcm3)
Finalmente, la temperatura para la formacin de roco (TJ es el

valor para el cual una masa de aire llega a estar safurada (e
e.) cuando se
enfra a presin constante y el mismo contenido de humedad. Una relacin

aproximada para la presin de saturacin de vapor de agua e. en funcin de la temperafura para la formacin de roco T es:
" =
2.7459t fO' '
- 4,278.6 "*n I Ta + 242'79

es exacta dentro de un
donde e, est en mb y Tu est en
oC. Esta relacin
rango de ms o menos de 0.57o de los valores observados dentro del intervalo de temperatura de OoC hasta 40oC.
Hidrologa / David
Cedeo
el
26
Las mediciones anteriores de la humedad atmosfrica se utilizan

frecuentemente en
anlisis del estado del tiempo para predecir la
probabilidad de precipitacin sobre un rea en particular.
CAMBIOS DE FASE
Para que el vapor se condense (pase del vapor al estado lquido) y
comience la formacin de la precipitacin, una cantidad de calor conocida como calor latente debe ser removida de la masa de-aire hmedo. El calor
latente de condensacin L" es igual al calor latente de evaporacin
se define como la cantidad de calor requerido para
L", el cual
convertir agua en estado
lquido a vapor a la misma temperatura; por consiguiente:
-|J
597.3
0-57 T
donde los calores latentes L" y L" estan en callgr y la temperatura T est
medida en oC.
Los meteorlogos usan las relaciones de humedad y los conceptos de
calor latente para obtener relaciones de presin-temperatura para el

enfriamiento de las masas ascendentes de aire hmedo. La rata de cambio de
temperatura con la elevacin en la affnsfera se denomina Gradiente Trmico.
Cuando las masas de aire hmedo no saturado, ascienden en la

ansfera, la humedad relativa se incrementa y al alcanzar cierta elevacin,
la masa de aire hmedo est completamente saturada, es decir, que

humedad relativa (H) alcanza
Ia
el I00%. Adicional enfriamiento del aire
resulta en la condensacin de la humedad y el calor latente de condensacin
27
es liberado, calentado el aire y por 1o tanto, disminuye el gradiente trmico
atmosfrico. Este intercambio de calor latente es la principal fuente

energa de los huracanes y ciclones tropicales.
de
Se ha observado que no existe una relacin definida entre la cantidad de vapor de agua
y la precipitacin resultante sobre una regin; es decir que
la
condensacin puede ocurrir formando nubes
sin que se
produzca
precipitacin sobre la superficie de la tierra, por lo tanto es necesario

considerar otros procesos climticos para analizar -los mecanismos de la precipitacin.
CANTIDAD DE AGUA PRECIPTABLE

La estimacin de la cantidad de precipitacin que puede ocurrir sobre
una regin donde existan condiciones favorables en el ambiente es una

informacin
til.
Este valor puede ser obtenido calculando la cantidad de
agua contenida en una columna de la atmsfera que se extiende desde la
superficie de la tierra
y el resultado
se conoce como cantidad de agua
precipitable (D), la cual se expresa generalmente en centmetros o pulgadas;

sin embargo esta cantidad no se puede remover totalnente de la atrnsfera por procesos naturales.
La
de rea
ecuacin para obtener
la
cantidad de agua precipitable en la
atmsfera se puede derivar considerando una columna de aire con una base
A.
La masa total de agua M contenida en sta columna de aire
hmedo entre la elevacin cero y alguna
alfiraZ
se puede expresar como:
28
Y = 1." P'dz
donde pu es la humedad absoluta. La ecuacin fundamental de la hidrosttica
se puede
escribir como:
dp - _ p^gdz
donde p- es la densidad total de lamezcla de aire seco y aire hmedo. Por
consiguiente, podemos despejar dz enla ecuacin hidrosttica, obteniendo:
-dP dz
P^8
Sustituyendo esta expresin en la primera ecuacin y observando que la razn de densidades es la humedad especfica
p, /
p,n ; tenemos:
M=-1
A
g "Po P^
[,o"ar=1[,.ndp g
"'
Recordando que la densidad del agua en estado lquido se puede expresar de la siguiente manera:
Masa
Volumen
M AD
29
donde D es la cantidad de agua precipitable (profrmdidad); finalmente tenemos
la siguiente ecuacin:
D-
pA
I
Pg
JP
[-"00,
Despus de introducir los factores de conversin correspondientes al
sistema mtrico, podemos obtener aproximadamente
la cantidad de
agua
precipitable de la siguiente manera:
D * [r'"uat
donde la profundidad D es en cm, la presin P se mide en mb
y la humedad
especfica q en gm/gm. El siguiente ejemplo ilustra un procedimiento de

integracin numrica de la frmula anterior para el calculo de la profundidad
de agua precipitable.
Ejemplo
4:
Cantidad de Agua Precipitable.
La siguiente tabla muestra los datos de elevacin, temperatura, presin atmosfrica
y presin de vapor. Calcular la profundidad D en cm de la
cantidad de agua precipitable en la atrnsfera contenida en una columna de 5.4
Km de altura.
30
Tabla: Datos Meteorolgicos.
Elevacin (z) Km 0.0 0.6
Temperatura (T)
Presin Atmosfrica (P) mb

1013.0
Presin de Vapor (e) mb
oc
15
11
7.0 5.0
3.8 3.2
942.0 875.0 812.0 753.0 697.0 644.0
t.2
1.8
1 J
-1
2.4
3.0
2.0
1.6
1.1
-5
3.6 4.2 4.8 5.4
-9
-13
-1',|
595.0
550.0 500.0
0.8 0.6 0.4
-20
Procedimiento:
Se divide la columna de la atmsfera en capas de 600 m de alfura, se
utiliza el valor promedio de la humedad especfica q
se calcula el
incremento en la presin atmosfrica AP para integrar numericamente la frmula de la cantidad de agua precipitable de la siguiente manera:
D"tqtp
31
Tabla de Clculos
Elevacin Z
Humedad Especfica
Valor
Promedio
q
Incremento
AP
Producto
9"0.622e1P
x
10-3
q' LP
103mb
L03m
gm/gm
xl0'3 gm/gm
mb
0.0 0.6
1.2 1.8
4.30
33r
2.71
3.81 3.01
71
269.80 201.00
67
63
2.45
1.65 1.43 1.06
2.58
2.Os
t62.54
120.95
2.4 3.0 3.6 4.2 4.8 5.4
59 56
53
t.s4
1.25
86.24
66.25 46.55
0.84
0.68
0.95
49
45
0.76 0.59
34.20
0.50
50
29.s0
Resultados:
D(cm) - E q.tP
1,017.0310-3
0.40
D r
7.02
cm =
plg
Hidrologa / David Cedeo CAUSAS
MECANISMOS DE
FORMACION DE
LA
PRECIPITACION
La condensacin de vapor de agua en gotitas en las nubes ocurre como
resultado del enfriamiento del aire a una temperatura por debajo del punto de
saturacin para el vapor de agua. Esto se logra generalnente a travs del

ascenso vertical a elevaciones donde la temperatura
y la presin son ms
bajas. La mitad de la masa de la atmsfera
esta locplizada hasta una altura
de 18,000 pies (5.48 Km) medida desde la superficie y contiene la mayora

de las nubes y la humedad. La condensacin puede ser producida por:
1.
Enfriamiento dinimico adiabtico (sin perdida de calor hacia los

alrededores).
2.
Mezclas de masas de aire que tienen diferentes temperafuras.
Enfriamiento por contacto.

4.
Enfriamiento por radiacin.
El
enfriamiento dinmico es el mecanismo ms importante en la
produccin de cantidades apreciables de precipitacin. El roco sereno, la

escarcha y la niebla son productores menores de precipitacin
y son causados
por enfriamiento por contacto y radiacin.

Los ncleos de condensacin deben estar presentes para la formacin de gotas en las nubes. El origen de estos ncleos es variado; entre ellos
tenemos: sal de los ocanos, polvo proveniente de suelos arcillosos, producto

de la combustin industrial y automotriz, cenizas volcnicas, etc.
33
y varan en
tamao desde 0.1p a 10r. (Nota:
lp = 1x10-6 metros)
Las gotas en las nubes inicialmente tienen un dimefro promedio de

0.01mm y solamente cuando estas exceden 0.5mm de dimetro es que ocrre
una precipitacin significativa. El proceso para que una gota pequea de
llua
horas.
(lmm) crezca sobre un ncleo de condensacin puede tomar varias
Cuando las masas de aire cargadas de humedad suben, estas se enfran y

expanden, y al ocurrir la saturacir el vapor de agua eomienza a condensarse
en los ncleos activos. El principal mecanismo para el suministro de agua a las
gotas crecientes en las etapas iniciales es la difusin de las molculas de vapor
de agua debido al gradiente de presin hacia las superhcies de las gotas.

Cuando la masa de las gotas se incrementa, estas comienzan a moverse con
respecto a las nubes. Sin embargo, existen otros procesos que afectan el crecimiento de las gotas hasta que alcancen un tamao suficiente (0.5mm -
3.0mm) de manera que superen la resistencia del aire
caigan como
precipitacin en cualquiera de sus formas. Estos mecanismos son el proceso de coalescencia y el proceso de cristales de hielo.
El proceso de coalescencia es considerado el mecanismo dominante de

la precipitacin en forma de
lluvia. Cuando las gotas de agua caen, las ms
pequeas
y lentas son absorbidas por las ms grandes, las cuales tienen una
velocidad de cada mayor, y el tamao de las gotas se incrementan a travs
de la colisin. Esto puede producir una precipitacin

nube).
abundante,
especialmente en cmulos clidos en las regiones tropicales (cierto tipo de
34
El proceso de cristales de hielo provoca la condensacin en los ncleos

congelados debido a las presiones de vapor ms
bajas. Los cristales de hielo
crecen en tamao a travs del contacto con otras partculas
y la colisin
produce
la formacin de nieve en forma de hojuelas.
Estas pueden
transformarse en gotas de lluvia, si al caer entran en contacto con aire en el cual la temperatura se encuentra por encima del punto de congelacin. Los ncleos de condensacin se pueden introducir artificialmente en las
nubes para provocar
la
precipitacin
bajo
- ciertas
condiciones.
Corrientemente se utiliza hielo seco artificiales. En la actualidad esta
yoduro de plata como ncleos

muy activa para
es un rea de investigacin
el control del clima, y todava existen muchos problemas tcnicos y legales por resolver relacionados con la precipitacin inducida o artificial.
ANALISIS DE DATOS DE PRECIPITACION

Los eventos de precipitacin son registrados en localidades especf,cas
utilizando pluviometros
y pluviografos. La interpretacin de los datos

Las variaciones en la
recogidos en las diferentes estaciones de medicin muestra la gran variacin
en el espacio
y el tiempo de la precipitacin.
distribucin y frecuencia de la precipitacin ocurren debido a las estaciones climticas
y alalocalizacin geogrfica, al igual
que las variaciones de un
evento individual de precipitacin se deben al tipo de tormentas, intensidad,
duraciny poca del ao; los vientos prevalecientes y la temperatura relativa
de la tierra con respecto al ocano tambin tienen su efecto sobre

precipitacin.
la

Generalmente se requiere una red de 5 a
35
l0
estaciones de medicin por
cada 100 millas cuadradas (25,000 hectreas) para registrar las variaciones
de la precipitacin. Pero el mantenimiento de estas redes de medicin

costoso y muchas veces ocurren fallas en el equipo; por
veces los registros estn incompletos.
es
lo tanto, algunas
Los datos de precipitacin se pueden utilizar para derivar las curvas de
intensidad-duracin-frecuencia (IDF), las cuales se utilizan generalmente

para obtener las caractersticas de las tormentas de diseo.
Se deben utllizar mtodos estadsticos (tales como la distribucin de
valores extremos) para ulr.lrzar la informacin requerida para la construccin
de curvas
IDF.
Uno de los modelos ms simples para estas curvas fue
propuesto por Steel (1960), el cual tiene la siguiente forma:
t+B
donde:
i
t
intensidad de la lluvia (plg/hr) duracin de la precipitacin (min)

constantes
A,B
La intensidad i representa el valor promedio de la profundidad de la

precipitacin acumulada P dividida por la duracin / registrada, es decir:
36
Los coeficientes A y B varin con la localizacin y el perodo de

retorno T en aos. Estos coeficientes se pueden obtener utilizando el mtodo
de regresin lineal, el cual requiere la transformacin de la funcin a una

lnea recta:
llB = _t iAA
INTENSIDAD tds/h'J
FRECUENCA
'100 aos
50 aos 25 aos
10 os
5 ens
DUFACI0N {minl
Figura N'6: Curvas Tpicas 5:
Intens dad-Duracin- Frecuenc ia

Curvas de Intensidad-Duracin-Frecuencia.
de
Ejemplo
La siguiente tabla muestra la precipitacin acumulada y la duracin

para lluvias con un perodo de recurrencia de una vez cada 10 aos. Ajustar
una curva IDF a los datos utilizando el mtodo de regresin lineal.
31
Tabla: Datos y Calculos
(min)
5 10 15
(plg)
(plg/hr)
'7 1
(hr/plg)
0.60
0.98
1.27 1.90
0.1389 0.1695 0.1961 0.2632
5.9
5.1
30 60
3.8 2.3
1.4
2.30 2.80
o.4348
t20
0.7t43
Solucin:
LacurvalDFtransformada
lnea recta, de la forma tpica:
I li: (l lA)t + (B/A) representauna

b
y -- mx + donde: x : y : m b : : t lli llA: B/A:
variableindependiente variabledependiente
pendiente
ordenadaenelorigen
Efectuando un anlisis de regresin lineal con los datos anteriores, se obtienen los siguientes resultados:
Pendiente:
m
:
0.0050
Ordenada
0.1190

Con wr Coeficiente de Determinacin: Rz de las constantes de 1a curva IDF es: A
m
3E
0.9986; por consiguiente el vaior
0.005
200
23.8
B = b.A
(0.1190)
(200) =
La curva de intensidad-duracin-frecuencia para un perodo de retorno
10 aos resulta ser: 200
t +
donde ia intensidad
23.8
se mide en plg/hora
y el tiempo / en minuros.
ESTIMACION DE DATOS FALTANTES

Muchas estaciones de precipitacin tienen datos taltantes en
sus
registros, causados por la ausencia del observador debido a f'allas en los

instrumentos, por
1o
tanto, a menudo es necesario estimar la precipitacin en
una estacin utilizando los valores registrados en las estaciones localizadas en

1os airededores,
Para la estimacin de los daros faitantes, Faulhus y Kohler (1952) propusieron el uso del promedio aritmetico simple con los datos de tres
estaciones cercanas,
el cual es adecuado cuando ia precipiracin anual

en ms del
rie
cada estacin no
estacin
difiere
i0
% de la precipitacin anual de ia
con registro incompleto; de tal manera que:
39
3"
P"
* Pz * P")
En caso contrario, es necesario ajustar las precipitaciones observadas
utilizando un factor de correccin igual a la razn de las precipitaciones anuales N
/ N
entre la estacin con datos faltantes y las tres estaciones
cercanas. Este procedimiento se denomina mtodo de la razn normal y la

frmula es la siguiente:
+(
la
N p,
NA
Np. N"'u
i/c
/ P")
en
donde
es
precipitacin estimada y ly' es la precipitacin anual
registrada en las estaciones.
Sin embargo, el mtodo ms utilizado en la actualidad para
la
estimacin de datos faltantes de precipitacin fue desarrollado por el National
Weather Service (1972), el cual esta basado en el promedio pesado de los

valores observados en los alrededores. En este caso, el peso
W es el
recproco de la distncia al cuadrado; es decir, la suma de las coordenadas al cuadrado de las estaciones medidas desde el punto de inters:
D2=X2*Y2
Por lo tanto, el factor de peso utilizado en el promedio pesado es:
W
1
D2
40
El valor estimado de la precipitacin utilizando n estaciones cercanas

(aproximadamente una estacin en cada cuadrante) resulta ser:
i=l
E r,w,
t=l
Ew,
Ejemplo
6:
Promedio Pesado
Estimar la precipitacin en la estacin A utilizando los datos de 5

estaciones cercanas.
Localizacin de las estaciones cercanas (sin escala)
41
Datos y Clculos.
Tabla:
P
Coordenadas
Estacin
(millas) (ple)
D2
w.103
P.W.103
X
0
Y
0 0
A
B
C
I
1.6 1.8
1.5
+4
+1
-3
+2 +6
+2
-J
a
20
37
50.0 27.0
80.0 48.6
D
E F
t3
18 8
76.9 55.6
tt5.4
2.0
-J
ttt.2
2I2.5
567.7
t.7
+2
t25.0
SUMA: 334.5
Solucin:
E p,w,
-A
tt
Ew, i=l
fn
ANALISIS DE DOBLE MASA
Latecntca de doble masa se utiliza para verificar la consistencia de los
= l.T|plg
datos de precipitacin. Este mtodo est basado en el hecho de que la precipitacin promedio acumulada para cierto nmero de estaciones no
es
muy sensitiva a los cambios en una estacin individual debido a que los
42
errores se compensan, mientras que los valores acumulados para una estacin
individual son afectados inmediatamente por los cambios que ocurren en la

estacin. Ests variaciones se producen por cambios en la localizacin de la
estacin, tipo de instrumento, mtodo de observacin; los cuales muchas

veces no se indican en los registros publicados. Si al graficar la precipitacin
anual acumulada para la estacin bajo investigacin contra la precipitacin promedio anual acumulado de las otras estaciones se obtiene una lnea recta,
se puede garantizar que los registros completos para esa estacin han sido obtenidos bajo las mismas condiciones; pero si existe un cambio de pendiente, generalmente se puede encontrar una explicacin al fenomeno. (por ejemplo:
la estacin fue movida de sitio). En este caso, los registros anteriores al

cambio de pendiente deben ser ajustados multiplicando por la razn de las
pendientes S2
/ Sl para hacerlos compatibles con los datos ms recientes.
Ejemplo
6:
Anlisis de Doble Masa
La siguiente tabla muestra los datos de precipitacin anual de la

estacin
X y el promedio
de precipitacin anual para 10 estaciones localizadas
en los alrededores (para 25 aos de registro).
a) b)
Determinar la consistencia de los registros de la estacin X; en caso

necesario, indicar el ao donde ocurre el cambio de pendiente. Calcular la precipitacin promedio anual para la estacin X utilizando los datos originales y efectuando los ajustes correspondientes.

Tabla: Datos y Clculos
AO
PRECIPITACION (cm)
43
PRECIPITACION ACUMULADA
ESTACION
ESTACIONX
49
38
PROMEDIO IO BSTAC]ONES
38 25
3'7
PROMED]O
1O
ESTACIONBS
1960 1961
49 87 123
148 181
38 63 100 126 150

181
t962
1963
36
25 35 38
1964 1965 1966

196'7
31
221
34 40 26
33 30 20 25
255 29s
214
244 264
1968 1969
32r
345 393 419
461
289 325
351
1970
19',71
48
36 26 24 21
26 42
l9'72 1973
197
375
49t
523 548
402 434 464 490 514

550 577
32
25
18 12
32 30
197 5
19'16
1977
566
5'78
r978
19'19
T6 18
36
602
618 636
27
25 26
1980
1981
602
628 6s9 691
728
20
6s6
680
1982
1983 19
18
3l
32 37
699
111
1984
= at) -uJ 9 o
F
@
)o OF a tlJ
@
tl
*1
x
ob
I o
= uJ
J
tr
o o
J f
o o l! o
J
I 9
z
o< At
= o 6
t c
tu
=
c\t
8383 !4, ta)

oooo oooa ro
ci
c.t
c\
x NolSvIS3 :VO\rInv{n3v Nol3vildlc3ud

Solucin:
45
a)
El an"lisis de doble masa muestra un cambio de pendiente en el
ao
1974, por lo tanto no existe consistencia en los datos de precipitacin
para la estacin X. Los Valores da las pendientes para cada tramo

son los siguientes:
^ J. ' ., r. '
523
434
1.20
( 7t7 - s23) (728 - 434)
0.66
b)
Precipitacin promedio anual para la Estacin X:
F=
Datos originales:
Li,,
ll
=t
F = 717 =
25
28.68 cmlao
Datos ajustados:
8_1
25
J:9!-+sr + (717-4e1 )]=
te.84cmlao
46
PRECIPITACION PROMDDIO SOBRE TINA REGION
El promedio de la profundidad de la precipitacin
sobre un rea
especfica se requiere a menudo para predecir la respuesta de una cuenca o
para desarrollar la tormenta de diseo. Existe tres mtodos bsicos para obtener los valores promedios sobre el
rea: Promedio Aritmtico, el
Polgono de Thiessen y el Mtodo de las Isoyetas.
El mtodo ms simple
es el promedio aritmtico de los valores P;
observados en las n estaciones de precipitacin loicalizadas dentro de la

cuenca, es decir:
P=
Li,,
h
=t
Este mtodo es satisfactorio si los medidores de precipitacin estn distribuidos uniformemente y las variaciones individuales de las lecturas no difieren mucho de la precipitacin promedio. El Polgono de Thiessen permite la distribucin de la precipitacin de acuerdo a las reas correspondientes a cada estacin. Para construir los polgonos se conectan las estaciones por medio de lneas rectas y se trazan
lneas perpendiculares que bisecten a las lneas conectoras para formar los
polgonos alrededor de cada estacin; a continuacin se miden o calculan las
reas
y laraz6n de las reas Ade
cada polgono dentro de los lmites de Ia
cuenca
y el rea total A, se utiliza para obtener la contribucin
de la
precipitacin P, en cada estacin a la precipitacin promedio sobre la cuenca'

en este caso:
47
E p,A,
P
E,
f)l
A.
l=l
' A,
Este mtodo es nico para cada red de medicin cuando las localizacin
de las
estaciones
es permanente y por consiguiente no permite
la
incorporacin de los efectos orogrficos (tales como cambios en la elevacin
del terreno) en la distribucin de la precipitacin. No obstante, es

probablemente el mtodo ms utilizado de los tres mtodos disponibles para obtener la precipitacin promedio sobre la cuenca.
El mtodo de isoyetas involucra el trazado de contornos de igual

precipitacin, o lneas isoyetas, sobre el reade la cuenca y es el mtodo ms
exacto de los tres; sin embargo, se requiere una gran cantidad de estaciones
de medicin para dibujar las isoyetas con precisin. Los clculos de la

precipitacin promedio sobre la cuenca estn basados en el valor promedio de
profundidad de la precipitacin entre cada par de contornos, luego
se
multiplica por el rea entre las isoyetas para obtener el volumen

total; en otras palabras:
de
precipitacin, finalmente se suman estos productos y se divide por el trea
E i=l E
D/
t=l AI
Volumen
A
4B
El mtodo de la isoyetas puede incluir los efectos orogrficos y la

morfologa de las tormenias; por lo tanto, el ftazado de isoyetas constituye un
mapa adecuado del patrn de la precipitacin.
Ejemplo
8:
Precipitacin Promedio sobre un rea.
Una cuenca de 23.6 millas cuadradas tiene un sistema de cuatro

estaciones de precipitacin, tal como se indica en el mapa (sin escala)
las
profundidades de precipitacin observadas en cada eitacin se muestran en
la tabla de datos. Determinar Ia precipitacin promedio sobre la

utilizando los siguientes mtodos:
cuenca
a) b) c)
Promedio Aritmtico
Polgono de Thiessen
Mtodo de Isoyetas
Tabla de Datos:
ESTACION A
B
C
PRECIPITACION (ple) 2.0

1.8
1.2
1.0

Tabla de Datos:
ESTACION PRECIPITACION (plg) 2.0
1.8
1,2
49
A
B
C
1.0
de la Cuenca (sin escala) mostrando 1a localizacin de fas

Mapa
Estaciones
Solucin:
a)
PromedioAritmtico:
(
D-l J
| ",
+
. P, *
1.O
Po
)
33
1.8
1.2 +
plg
50
b)
Polgono de Thiessen: Tabla de Clculos:
Pi
Ai
Ai/AT
0.064
0.305
(ple) 2.0
1.8 1.2 1.0 Suma:
(mi')
1.5
Pt(At/Ar) (ple)
0.13 0.55
7.2
5.1
0.216
0.415 1.000
o.26 0.42
1.35
9.8
23.6
I.3s plg
Figura N'7: Poligono de Thiessen.
51_
c)
Mtodo de Isoyetas
Tabla de Clculos
ISOYETA (ple)
2.0
5.1 1.8
A.
(mi')
(nls)
Vol=PA lt
(plg-mi')
1.9
9.69
9.8
i.5
t4.'7
r.2
3.1 1.0
1.1
3.41
5.6 0.4
0.'7
3.92
Suma:
N:23.6
P
VOL4
:31.72
Vol,
3172
23.60
Ar
1.34
plg
0,t
+Yalor estmado
Figura N'B: Mtodo de Isoyetas.
Hidrologa / Dad Cedeo
52
EVAPORACION Y TRANSPIRACION
Evaporacin es el proceso por medio del cual el agua en estado lquido
o slido es transformada en vapor de agua, el cual se mezcla con el aire de

la atmslera.
La evapotranspiracin se considera
combinada de vapor de agua
separadamente como la prdida
a travs de la superficie de las plantas
(transpiracin) y la evaporacin de la humedad del suelo.
El conocimiento de los procesos de evaporac-in es importante para

predecir las prdidas de agua debido a la evaporacin que ocurrirn en un
lago o embalse. Aproximadamente
el6l%
de la precipitacin promedio anual
sobre la superficie terrestre regresa a la atmsfera a travs de la evaporacin
evapotranspiracin (tal como se indica en la figura del balance hdrico

a
promedio global anual). Sin embargo, las variaciones en la evaporacin
travs del continente pueden ser muy grandes, ya que existen regiones
desrticas o ridas donde la evaporacin anual puede exceder la precipitacin
promedio anual.
En el caso de evaporacin desde la superficie de un lago, la prdida de
agua es funcin de radiacin solar, temperatura del agua y el aire, diferencia
en la presin de vapor entre el agua y la capa de aire sobre el lago y la

velocidad del viento sobre el
lago. Cuando ocurre evaporacin dentro
de un
sistema cerrado a temperatura constante, la presin dentro del recipiente se incrementa debido al aumento en la presin parcial de
vapor. La evaporacin
contina hasta que la presin de vapor de la capa de aire sea igual a Ia presin
de vapor de la superficie del lquido; en este instante se dice que la masa de
tr2
aire est saturada a esa temperatura y no ocurre ms evaporacin. Este

estado de equilibrio no se alcanzara si el recipiente estuviera abierto a la atrnsfera; en cuyo caso, el lquido se evaporara completamente. Se requiere
energia trmica para incrementar la energa libre de las molculas de agua

para que estas pasen a travs de la interfase gasJquido. La cantidad de calor
requerida para convertir agua en estado lquido a vapor se denomina calor

latente de evaporacin.
Cuando la evaporacin contina sobre una s+tperficie horizontal de

agua, la acumulacin de molculas de vapor de agua produce un incremento
en la presin de vapor en el airejustamente sobre la superficie del agua, hasta
que evenfualmente comience la condensacin. El aire est saturado cuando la rata de condensacin es igual a la rata de evaporacin y adems, la presin
de vapor es igual a la presin de vapor de saturacin. Sin embargo, existen
varios procesos de transporte convectivo que afectan el transporte de vapor
(tales como corrientes de aire vientos) las cuales evitan que ocurra el equilibrio en el ambiente (sistema abierto).
La evaporacin solamente es de gran preocupacin en la planificacin
de grandes proyectos de recursos hidrulicos

abastecimiento de
en los estudios
de
agua.
Durante perodos tpicos de tormentas, con
intensidades de precipitacin de 0.5
plg/hr, la evaporacin
se encuentra en el
orden de 0.01 plg/hr y por lo tanto se puede despreciar en los estudios de

caudales de inundaciones
y en las aplicaciones de diseo de drenaje urbano.
La evaporacin ha sido estudiada extensivamente en los Estados Unidos
54
especiaknente en los proyectos de investigacin de evaporacin efectuados en
el Lago Hefner, Oklahoma, por Marciano y Harbeck (1954).
Existen tres mtodos primarios para estimar la evaporacin desde la superficie de un lago:
a) b) c)
El mtodo de balance hdrico El mtodo de transferencia de masa El mtodo de balance energtico
METODO DEL BALANCE HIDRIC"

EVAPORACION
"A*A-
DETERMINAR LA
El mtodo de balance hdrico para obtener la evaporacin de un lago

esta basado en la ecuacin de continuidad
hidrolgica. Asumiendo que el
cambio en el almacenamiento AS, la escorrentia superficial de entrada I y de
salida O, la infiltracin F hacia el flujo subterrneo
y la precipitacin P
pueden ser medidas; la evaporacin se puede calcular de la siguiente manera:
E=P+I-O-F-A,S
Este procedimiento es simple en teora, pero la evaluacin del trmino
correspondiente a la infiltracin hace que este mtodo sea muy difcil de implementar. Las dificultades con este procedimiento resultan de los errores
en la medicin de la precipitacin
y caudales de entrada y salida, cambios en
el
almacenamiento
rata de infiltracin. Se han obtenido muy buenos
resultados con este mtodo en el Lago Hefrrer con errores del 5% a
l0%. Es
importante sealar que el Lago Hefner fue escogido entre ms de 100 lagos
Hidrologa lDawd
Cedeo
55
embalses, ya que es uno de los tres o cuatro lugares que satisfacen mejor
los requerimientos del balance hdrico.
METODO DE TRANSFERENCIADE MASA

Las tcnicas de transferencia de masa estn basadas principalmente en
el concepto de transferencia turbulenta de vapor de agua desde la superficie
del lquido hacia la atmsfera. Se han desarrollado rumerosas frmulas

empricas para obtener
la rata de evaporacin
co?no una funcin de la
diferencia de presin de vapor
y la velocidad del viento
sobre el lago o
embalse. La mayora de estas ecuaciones se pueden escribir de manera

similar a la ley de Dalton:
E = l"* - eo)la * bu)

donde:
evaporacron.
eea
presin de vapor en la superficie del agua. presin de vapor a cierta altura sobre la superficie velocidad del viento
constantes empricas
1as
a,b
Un obstculo para comparar
diferentes frmulas de evaporacin es
la variabilidad en la medicin de la altura para ea y
z.
Si reducimos todas las
frmulas existentes a efectuar las mismas mediciones a una alfura de 2 metros
(6.5 pies) para la velocidad del viento y la presin de vapor y tomamos en
Hidrologa / D avid Cedeo
56
cuenta la diferencia de alrededor de 30% entre la evaporacin medida en un tanque evapormetro y la evaporacin actual sobre un embalse, la discrepancia
entre las diferentes frmulas se reduce considerablemente. La frmula emprica con la mejor base de datos es para el lago Hefner, la cual tambin funciona para el lago Mead, fue presentada por Harbeck y
Meyers (1970) y tiene la siguiente forma:
E= donde:
E N N N
ew
e2
N url e,
- er)
rata de evaporacin (cm/da)

constante emprica
0.012 para el Lago Heher 0.0118 para el Lago Mead presin de vapor en la superficie del agua (mb)
presin de vapor medida a

superficie (mb)
2 metros sobre
la
u2
velocidad del viento medida a 2 metros sobre Ia

superficie del agua (m/s).
METODO DE BALANCE ENERGETICO
El mtodo ms preciso y complejo para determinar la evaporacin

utiliza el balance energtico de un lago. La ecuacin general para el balance
energtico de un lago en langley/da (1 langley expresar como:
- Lv I callc) se puede
57
o
donde:
Qn Qn
-o.-o=o^-o
radiacin neta absorbida por el cuerpo de agua
transferencia de calor sensible (conduccin y conveccin hacia la atmsfera)
Q"
Qu
energa utilizada para evaporacin
incremento en la energa almacenada en el cuerpo

de agua
Q"
energa transportada por adveccin del caudal de

entrada y salida
Figura N'9: Balance Energtico de un Lago
Por otro lado, la radiacin neta Q" absorbida por el cuerpo de agua es equivalente a:
o =o -o -o-b
58
donde: Q, Q. Qn
radiacin solar de onda corta radiacin reflejada de onda corta

radiacin de onda larga reflejada hacia la atmsfera
Si recordamos que L"

(cal/gm)
representa
el calor latente de vaporizacin
y p larazn
entre la perdida de calor por conduccin y la prdida de
calor por evaporacin, tenemos que:
Q,*Q,-Qs oZ,(1 +P;

y p es la densidad
del agua
donde E es la rata de evaporacin (cm/da)
(gm/cm3). La razn de Bowen p se utiliza como una medida del calor

sensible transferido y puede ser calculada de esta manera:
='.,
(r-r\ro\
l.;=j l,r*,l ='li=)

presin atmosfrica (mb) temperatura del aire (oC) temperatura de la superficie del agua (oC) presin de saturacin de vapor a la temperatura de
(r-r\
donde: P : Tu T, es ea y :
la superficie del agua (mb)

presin de vapor del aire (mb)
constante psicomtrica (mb/oC)

Nota:
59
0.66 P / 1,000
La aplicacin del mtodo de balance energtico requiere la medicin de

la radiacin total de entrada neta. La razn de Bowen fue propuesta debido
a que la transferencia de calor sensible no puede ser calculada
fcilmente. El
mtodo fue aplicado al Lago Hefner y al Lago Mead y fue utilizado para
evaluar los coeficientes empricos para el mtodo de transferencia de masa y
para interpretar los datos de evaporacin para un tanque evapormetro

colocado en
el Lago Hefner. El
mtodo de balance energtico
es
tericamente el ms preciso, pero requiere la coleccin de grandes cantidades

de datos atmosfricos. Para evitar este problema, se han desarrollado otros
procedimientos, tales como
el
tanque evapormetro para estimar la
evaporacin de un lago poco profundo y los mtodos combinados.
TANQUE EVAPORTMETRO La evaporacin puede ser medida utilizando un tanque estandarizado

tipo A, el cual es un tanque cilndrico abierto de hierro galvanizado de 4 pies
de dimetro
10 pulgadas de profundidad, colocado a 12 pulgadas sobre el una
suelo. Para estimar la evaporacin, el tanque se llena de agua hasta pulgadas. El nivel de la superficie de agua se mide diariamente
altura de 8 pulgadas y se debe rellenar cuando la profundidad desciende a 7
la
evaporacin se calcula como la diferencia entre los niveles observados,

ajustados para tomar en cuenta la precipitacin medida en un pluviometro
cercano. La evaporacin en un tanque evapormetro es mayor que la
Hidrologa / David Cedeo evaporacin actual en el lago
60
y debe ser ajustada para tomar en cuenta la

de
radiacin
y los efectos del intercambio
calor. El factor de ajuste se
denomina coeficiente del tanque, el cual varia de 0.64 hasta 0.81 con un valor
promedio de 0.70. Observe que este coeficiente varia con la exposicin a la

radiacin y las condiciones climticas y debe ser utilizado solamente para una estimacin aproximada de la evaporacin de un lago por medio de la siguiente
frmula:
EL = Cr'E,
donde:
EL
CT ET
evaporacin estimada en el lago coeficiente del tanque (Cr = 0.7) evaporacin medida en el tanque evapormetro
METODOS COMBINADOS
Penman (1948) fue el primero en utilizar las mejores caractersticas de
los mtodos de transporte de masa y balance energtico para derivar una

relacin para la evaporacin de la superficie del agua de un lago que fuera
relativamente sencilla de calcular. La ecuacin de Penman [(en unidades de
eneryial (nea'tiempo)
es :
F.= -h
A +Y
O+
a +y
vE
6L
donde: Eh A :
flujo de calor latente debido a la evaporacin

pendiente de la grfica de presin de saturacin de
vapor es en funcin de la temperatura T (mb/oC)

constante psicomtrica (mb/oC)
Qn Eu
radiacin neta absorbida poder de secado del aire
De puede uttl:zar la siguiente frmula para el clculo del flujo de calor
latente debido a la evaporacin
Eh = p L"E
donde: E,
flujo de calor latente debido a la

lener ga I (r
ea'
evaporacin
tiemp o)l
densidad del agua (masa/volumen)
L" E
calor latente de vaporizacin,
generalnente
evaluado a la temperafura del aire (energa/masa). rata de evaporacin (profundidad/tiempo)
En la prctica, es comn medir el parmetro A a la temperatura del aire
y no a la temperatura de la superficie del agua. Este parmetro A (en mb/oC)

se puede obtener diferenciando
la expresin para la presin de saturacin de
vapor en funcin de la tempetratura T (en oC); es decir:
A =
d""
dr
(2.7489xrc8)'G,27g.6)
(r.rorrn\2
\/
'\r*z+zts) "*o
q,zts.a
52
Segn Brutsaert (1982), el poder de secado del aire [en unidades de

energa
/ (rea.tiempo)l
se puede evaluar de la siguiente manera:
Eo = pL"(".0,)(,,,-"")
donde:
E^
p
poder de secado del aire

densidad del agua
L"
calor latente de evaporacin

constantes empricas de transferencia
a,b
u
er"
velocidad del viento

presin de vapor de saturacin a la temperatura del
aire
e2
presin actual del vapor en el aire
"
H
too
sa
humedad relativa en porcentaje
%)
La ecuacin de Penman tiene la ventaja de que la temperafura del agua

o del suelo no se requiere en los clculos. Se ha encontrado que esta ecuacin
es muy
til para estudios de evapotranspiracin, en los cuales
es muy
difcil
la
determinar
la
temperatura superflcial
de la vegetacin. Cuando
temperatura de la superficie del agua se puede medir, el procedimiento del

balance energtico - razn de Bowen es probablemente mejor porque evita la
necesidad de utilizar los coeficientes empricos de transferencia (a
bu).
Hidrologa
lDaidCedeo
Evaporacin utilizando la ecuacin de Penman
63
Ejemplo
9:
La ecuacin emprica de transferencia de masa para cierto lago es:
E=
0.0106f, \
I +0.r,\(, /\'
-e ') l
donde E se mide en pulgadas/da, u en millas/horas y las presiones de vapor

en
mb. Estimar la evaporacin
de ese lago utilizando la ecuacin de Penman,
para una temperatura del aire de 90oF, velocidad hel viento de 20 MPH,
humedad relativa de30%, un flujo de radiacin neta de 400 langley/da y una
presin atmosfrica de l000mb (nota: 1 Ly
: I callcm2).
Datos: T u : H Qn P :
Solucin:
90oF
20MPH
30% 400
Lylda
1,000 mb
a)
Coeficiente emprico de transferencia (a
b u):
E = ('.u,)("--,,)
00roe
(r.0,,)(,--".)
Por lo tanto, el valor del coeficiente emprico es:
/\/\/\+
b/t = 0.0106 I l+0.lrl \/\/\/mb'dta

1a
= o.oloe
f t * ol.zol= ornts
ptg
64
b)
Temperatura (T):
T("C) -
rr".F)-32
lI
s lno-rr) = s\ i
i2.2.c
c)
Constante Psicomtrica (y):
0.66P _ 1,000
0.66(1,000)
1,000
= 066
mb
0c
d) ^
Pendiente (A):
du,
dT
2.7489
to. .
/ \z \T+2a2.7e)
"'.n
)'
(4,278.6) f
av
'l
- 4.27s.6
t
--
T*zztg
1.1761 x 1 012
^=
(rr, .,o 2 .79
- 4,278.6 32.2 + 242.79
) = 2.72
mb
).c
e)
Presin de Saturacin de Vapor a la temperatura del aire (e,"):
"o = 2.7489x
108 exp
( --o''''r u ) -^'lr*242.7s)
,o
2.748er ros exp (
!''!|-u --) 32.2 +242.79 )
48.7 mb
Hidrologa i Dad
Cedeo
H 1oo "so 30 (48.1 ) =
65
Presin actual del Vapor del aire (e):
e = '
loo
l4.4mb
g)
Calor Latente de Vaporizacin (L"):
L, = h)
ss7.3
- 0.57 T =
se7.3
0.s7
(32.2)
sls
Poder de Secado del aire (E"):
pL"( o.r,\(,""-,,\ Eo = --\"

i \
sct
r)
".=
[' #\,,
#)(00,,* #h)(or,*u
=
144nb)
cat'ptg.l( ,ro:y\ =(uroon "' \ cm3.da)\ pts )
1.576 cat = t.si6 cm2'da
Lv
dta
i)
Flujo de Calor Latente de Vaporizacin (ecuacin de Penman):
E.=Lo*\E /' +y
^
'tl
A +T
E, '
2.72 loool) + dla ) 2.72 + 0.66 \
( ,.rru tr_) dial 2.72 + 0.66 \
0.6
66
E.
uro L!
da
630
cm2.da
Rata de Evaporacin (E):
Eh = p L"E
Por consiguiente, tenemos que:
E.
630
cal
PL"
('#) ("' cm)

cal _t
cm2 . da
1.09
cm
dla
Conversin de unidades (sistema ingls):
E = ( ,0, ,,'\( \ dral\
rprs ) = 254cm)
043
plg
dla
EVAPOTRANSPIRACION Evapotranspiracin @T), algunas veces llamado uso consuntivo
evaporacin total, es la combinacin de evaporacin sobre la superficie del
suelo y la transpiracin a travs de los poros (estomas) de las hojas de las plantas. Los mismos factores que afectan la evaporacin de una superficie de
agua
(tal como un lago)
tambin gobiernan
la
evapotranspiracin,
principalmente el abastecimiento de energa y el transporte de vapor. En
adicin, un tercer factor afecta el mecanismo de evapotranspiracin: el

abastecimiento de humedad en la supercie de evaporacin. Cuando el suelo
se seca, la rata de evapotranspiracin se reduce a un nivel inferior al que
existira en un suelo bien irrigado.
La capacidad de campo del suelo es el contenido de humedad por

encima del cual el agua drena por gravedad y el punto de marchitez es el
contenido de humedad por debajo del cual las plantas no pueden extraer agua
del suelo.
Para la mayora de las plantas, la transpiracin ocurre solamente

durante las horas de luz solar mediante el proceso de fotosntesis, la cual
produce variaciones diurnas en el nivel fretico poco profundo en zonas con
vegetacin densa. La evapotranspiracin alcanza un valor mximo si el suministro de agua hacia las plantas y superficie del suelo es ilimitado. La prdida mxima posible esta limitada por condiciones meteorolgicas y se
denomina eyapotranspiracin potencial (Thornthwaite, 1948)
es
aproximadamente igual a la evaporacin que ocurrira en una superficie

grande de agua, tal como un lago. Por consiguiente, los mtodos utilizados
para arnlizar la evaporacin discutidos anteriormente, tambin se pueden usar para predecir la evapotranspiracin potencial.
68
INFILTRACION
El proceso de infiltracin ha sido ampliamente estudiado y representa
un mecanismo importante para el movimiento del agua hacia el suelo bajo la
accin de la gravedad y fuerzas capilares. Horton (1933) demostr que

cuando la rata de precipitacin i excede la rata de infiltracin/, el agua se
infiltra en las capas superficiales del suelo en una proporcin que

generalmente disminuye con el tiempo. Para cualquier suelo, existe una curva limtrofe que define la rata de infiltracin mx-ima posible en funcin
del tiempo. La rata de infiltracin depende de manera muy complicada con
la intensidad de precipitacin, tipo de suelo, condicin de la superficie y

cobertura de la vegetacin.
Cuando existe una precipitacin excedente, es decir,
la rata de
precipitacin es mayor que la rata de infiltracin, la infiltracin seguir la
curva limtrofe mostrada en la Figura 10, la cual se denomina curva de

capacidad de infiltracin del
suelo. En esta grfrca
se puede observar que la
capacidad disminuye con el tiempo hasta que alcarlza un valor constante. Esta
disminucin se produce por el llenado de los poros del suelo con agua,
reduciendo la succin capilar. Por ejemplo, en pruebas controladas se ha
demostrado que esta disminucin es ms rpida y el valor constante es menor
para suelos arcillosos que para suelos arenosos.
bv
i,f
ALh{ACENAMIENTO INICIAI {INTEHCEPN Y AWACNAMIENf EN DEPfiESIONES]
O
INTENSIOAD DE PREEHTACIBN VOLUMEN DE
ESCCFFENiIA
lNRLIFACIIJN
VOLUMEN DE
II}
iNFIL;RAI:IN
Figura N'10: Mdelo Conceptual de -Lnl].Icracron de -Horton.

RATA DE INFILTRACION
El concepto hidrolgico de capacidad de infiltracin es emprico y est
basado en observaciones efectuadas en la superficie del suelo. Cuando la intensidad de precipitacin
es mayor que la rata de

1a
infiltracin f, Horton
(1940) sugiri la siguiente forma para
ecuacin de infiltracin:

donde:
70
f
f"
f" k
capacidad de infiltracin (plg/hr) capacidad inicial de infiltracin (plg/hr) capacidad final de infiltracin (plg/hr) constante emprica (hr-1)
El volumen total de infiltracin F

ecuacin de Horton y esta dado por:
se puede obtener integrando la
F(t\ = r"t
. ("*"1 ,, -,-*,,
de
Una limitante de la ecuacin de Horton es que la capacidad
infiltracin disminuye como una funcin del tiempo, sin tomar en cuenta la
cantidad de agua disponible para infiltracin. Es decir, que la ecuacin asume
la formacin de lagunas charcos de agua en la superficie que del suelo y una
reduccin en la capacidad de infiltracin, independientemente de que la

intensidad de precipitacin
exceda no el valor calculado para la capacidad
inicial de infiltracin f". Por ejemplo, es muy comn que la capacidad de

infiltracin de suelos arenosos sea mayor que la intensidad de precipitacin,
con valores de la capacidad final
que se encontraran en el rango de 10 a
20 plglhr. En muchas ocasiones, las lluvias muy fuertes no alcanzan estos

valores; en consecuencia, toda la lluvia se infiltrara en el suelo, es decir que
f : . Por consiguiente, la capacidad de infiltracin

duracin de la infiltracin.
debe reducirse en
proporcin al volumen acumulado de infiltracin, no en proporcin a la
'7L
Rubin y otros (1963, 1964) demostraron que las curvas observadas por
Horton se pueden predecir tericamente si se conoce Ia intensidad de la lluvia,
las condiciones iniciales de la humedad del suelo y las curvas caractersticas
para el suelo no saturado. Ellos indicaron que la rata de infiltracin final es
numricamente equivalente
a la
conductividad hidrulica para suelos
saturados. Adicionalmente, Rubin mostr que la formacin de charcos de

agua
lagunas en las superficie ocurrir solamente si
la duracin de
la
precipitacin es mayor que el tiempo requerido para que el suelo se sature en
la superficie.
Eiemolo 10:
Ecuacin de infiltracin de Horton.
Se estima que la capacidad inicial de infiltracin
de una cuenca tiene
unvalor de 1.5 plg/hry la constante emprica
se asume que es 0.35
hr
r;
adems, se ha observado que la capacidad de equilibrio
f"
es 0.2 plg/hr.
Utilizar la ecuacin de Horton para encontrar:
a) b)
Los valores de la capacidad de infiltracin

instantes:
para los siguientes
10 min, 30
min,l hr, 2hr y 6 hroras.
El volumen total de infiltracin durante el perodo de 6 horas.
Observacin: Durante el intervalo de tiempo 0 <
t < 6 hr, la intensidad de la
precipitacin
es mayor que la rata de
infiltracin
/;
es
decir:
i > /.

Solucin:
1)
f [plsrhr]
K = 0-35 hr
t [horasJ
a)
Capacidad de infiltracin/ (Ecuacin de Horton):
r = r, - (r,-r"),
f = 0.2 + (
1.5
n'
- 0.2)e-035t -
0.2 + 1.3 exp(-0.35t)
(horas)
(ple/hr)
u6
U2
1
t.43
1.29
l.r2
0.85
2
6
0.36
'73
b)
Volumen de infiltracin
-F
(integracin de la ecuacin de Horton en el
intervalo 0 < t < 6 horas):
F = {"' f o, = [,' t 0.2 + t.3e
o'35t
) d
F=
F =
1.2
or,
-(-J.3-') \ o.ls /
"-0.,,,
0
- 337 exp(-2.10) + 3.71 =
4.46
plg
OTROS METODOS PARA CALCULAR LA INFILTRACION
Se han desarrollado otras frmulas para calcular la infiltracin

utilizando soluciones analticas para la ecuacin de flujo no saturado. Por
ejemplo, Philip (1957) desarrollo una ecuacin de la siguiente forma:
| f = 2 r-rrz +
F = A ttt2 + Bt
donde: f F
: : A, B :
capacidad de infiltracin (plg/hr)
volumen acumulado de infiltracin (plg) constantes relacionadas con movimiento del agua
el tipo de suelo y
'74
i.
[plgJhrl
PFECIPITACIf]I Tf]TAL
rHErcEf
YOLUItEH DE IHFILTRACIH
t [horasl
Figura N"11: Mtodo del Indice
Por otro lado, cuando no existen mediciones detalladas de las prdidas
de agua
en e1 caso de cuencas urbanas, las cuales son altamente
impermeables; el uso de procedimientos empricos producen resultados

satisfactorios en la mayora de estas situaciones. Se ha observado que la infiltracin representa un porcentaje variable de la precipitacin total que cae enuna cuenca. En la mayora de los esrudios de drenaje ubano y control de
inundaciones se utiliza la ecuacin de Horton, o en su reemplazo, mtodos
ms simples para predecir los volmenes de infiltracin. El ndice iD es el
75
mtodo ms elemental
se calcula encontrando la diferencia entre la
precipitacin total y la escorrientia superhcial registrada en un hidrograma de descarga. El mtodo del ndice
asume que la prdida de agua se distribuye
uniformemente durante el evento de precipitacin. El mtodo del ndice o

para infiltracin se ilustra por medio del siguiente ejemplo.
Ejemplo 11:
Mtodo del ndice para calcular la infiltracin.
Utilice los datos de precipitacin mostrados en la tabla para determinar

el ndice
de una cuenca que tiene un rea de 0.875 millas cuadradas, si el
volumen de escorrentia medido fue de 228.7 acres-pie; adems, calcular la profundidad de precipitacin total y la infiltracin (ambas en plg).
Tabla de Datos:
Intervalo de Tiempo (hr)
Intensidad de Precipitacin
(ple/hr)
1.4
o-2 2-5
5-7 7-10
10 12
2.3
1.1
0.7
0.3
Hidrologa / Dad Cedeo Solucin:
76
El primer paso requerido para la solucin del problema involucra ia

construccin de una grfica con los datos de precipitacin.
a)
Profrutdidad de Escorrentia Superficial R
R=
Volumen
Area
(228.7 acres.pies) (43,560 pies2/acres) (12 ptg/pie)
(0.875 mi2
( 5,280 pieslmilla)2
4.9 ptg
R =
77
b)
Indice o:
La rcta de infiltracin ndice o

lnea para la cual
se puede encontrar por ensayo y de agua por encima de la
error. La escorrentia superficial R es el volumen
i : : observe que en algunos perodos el ndice es

o
< 1.1 plg/hr; tenemos:
mayor que la intensidad de precipitacin i. Asumiendo que el rango para la

rata de infiltracin es 0.7 plg/hr <
R=
(1.4
- o)(2 -0) +(2.3 - o) (s-2)+(1.1 -o)(7-s)

= 11.90 - 7O
plglhr
4.9 = (2.8 + .9 + 2.2) - (2 + 3 +2)

11.90
4.9
1.0
Observe que este valor esta dentro del rango asumido. Estos clculos
nos indican que debajo de la lnea atrazos para la
cual o
1.0 plg/hr, la
a
precipitacin se infiltra en el suelo y que el volumen encima de la lnea

trazos corresponde a la escorrentia superficial.
c)
Precipitacin Total P:
P= 1.4(2-0) + 2.3(5 -2) *
1.1(z
-s)
+ 0.7(10-7) + 0.3(12-10)
P =
2.80 + .90 + 2.20 + 2.10 + 0.60
14.60
plg
78
d)
Infiltracin F:
F = P - R = 14.6 - 4.90 = 97 plg

Mtodo alterno para el clculo de la infiltracin:
' = o(7 - 0)+0.7(10 - 7) * 0.3(12 F = 7 +2.1 + 0.6 =

9.7 ple
10)
ESCORRENTIA ST]PERFICIAL Cuando la precipitacin cae sobre
la
superficie de
la tierra,
se
distribuye de diferentes maneras; inicialmente la precipitacin comienza a rellenar las depresiones del suelo, infiltrarse para recargar la humedad del
suelo y agua subterrnea, o viajar como flujo subsuperficial hasta alcanzar
una corriente de agua. El almacenamiento en las depresiones se satisface en los perodos iniciales de la tormenta, seguido a continuacin por la capacidad
de humedad del suelo. Eventualmente comienza el flujo superficial

escorrenta, precipitacin
el cual ocurre
solamente despus que
la
intensidad de la
sobrepase la capacidad de infiltracin
f (i > f).
al
flujo
El
concepto clsico de generacin de corrientes debido
superficial sobre la tierra fue propuesto por Horton (1933), quien indic que
el flujo superficial estaba distribudo de manera generalizada sobre el terreno.
Posteriormente otros investigadores analizaron la gran heterogeneidad que

existe en cuencas naturales e introdujeron el concepto de contribucin parcial
del rea superficial (Betson, 1964). Este concepto reconoce que solamente

algunas porciones de la cuenca contribuyen regularmente al
'79
flujo superficial
hacia las corrientes y que no ms del I0% aproximadamente del rea de la

cuenca en estado natural contribuye al flujo superficial. En ambientes urbanos
con grandes zonas impermeables, el porcentaje de contribucin al flujo

superficial puede ser mucho mayor.
Un segundo concepto importante en la generacin de escorrenta

superficial es el movimiento de agua debajo de la superflcie del terreno en las
capas superiores del suelo sin alcawar la zona de saturacin, el cual se
denomina
flujo subsuperficial. Freeze (1972) concluy que el
flujo
subsuperficial era una componente significativa solamente en terrenos con

pendientes convexas que abastecen canales profundos
y solamente en el caso
de suelos muy permeables. En pendientes cncavas, los valles saturados

creados por el ascenso del nivel fretico contribuyen con el
flujo superficial,
el cual generalmente excede el flujo subsuperficial.
El flujo superficial producido por la precipitacin excedente se mueve

hacia abajo en direccin de la pendiente de la superficie del terreno hasta
alcanzar los pequeos canales de drenaje quebradas, los cuales fluyen hacia
corrientes mas grandes, transformandose generalmente en ros. Cuando el
flujo alcanza la corriente principal, las velocidades y
las profundidades del
flujo se pueden medir en una seccin transversal particular a travs del

tiempo,
1o
cual nos permite obtener el hidrograma, es decir, una grfica de
la descarga caudal en funcin del tiempo. La forma actual
y los tiempos del
hidrograma estn determinados en gran parte por el tamao, forma, pendiente
almacenamiento en
la cuenca; y por la intensidad y duracin de la
80
precipitacin. Estos factores se analizan con mayor detalle en la siguiente

seccin donde se estudian las relaciones entre precipitacin y escorrenta.
Despus que termina la precipitacin, el volumen almacenado en la cuenca se
libera hacia las corrientes, completando el ciclo de la tormenta.
Los canales pueden contener cierta cantidad de flujo base, el cual

proviene del flujo subterrneo
ausencia de precipitacin.
las contribuciones del suelo, an en la
La
descarga producida
por la precpitacin
excedente, es decir, la precipitacin total menos todas-las prdidas, constituye
el hidrograma de escorrenta directa. Por lo tanto, se considera que el

hidrograma total est formado por la escorrenta directa mas el flujo base. La
duracin de la precipitacin determina la porcin del rea de la cuenca que
contribuye al flujo mximo, mientras que la intensidad de la precipitacin determina la magnitud del caudal mximo resultante. Si la precipitacin
mantiene una intensidad constante por un perodo muy largo de tiempo, se
produce un almacenamiento mximo y se alcanza una condicin de equilibrio para la descarga. Esta condicin de equilibrio se logra en muy raras ocasiones
en la naturaleza debido a la variacin de la intensidad y distribucin en el tiempo y el espacio de la precipitacin sobre la cuenca.
MEDICION DEL CAT]DAL: AFOROS

Para determinar el caudal en un ro se utiliza una tcnica denominada
aforo, la cual consiste en dividir el ancho total de la corriente en un nmero

conveniente de secciones
y la velocidad
media en cada seccin se mide
81
utilizando un molinete. Estas mediciones se pueden efectrar por vadeo cuando los ros son poco profundos, desde un bote, puente cableva.
Se ha observado que la velocidad media V en una seccin ocurre

aproximadamente
0.6
D, medida
desde la superficie del agua, por lo tanto
la velocidad se mide a ese nivel con el molinete. Sin embargo, cuando la

profundidad
D enla
seccin es mayor de cierto valor, se recomienda utilizar
el promedio de las velocidades medidas
0.2
Dy
0.8 D.
Tambin se puede estimar el caudal en una corriente de manera

rudimentaria utilizando un flotador; en este daso el caudal
ser igual a la velocidad
Q en cada seccin
del flotador corregida utilizando un coeficiente
C, el cual tiene un valor aproximado de 0.85, multiplicada por el rea,4 de
la seccin transversal del ro sobre la cual

flotador; es decir:
se midi
la velocidad con
el
Q, = C'V,'A,
Por lo tanto, el caudal total es:
o^ =
Io
Cuando escogemos un sitio para el establecimiento de una estacin de
medicin, podemos obtener informacin sobre el caudal para diferentes

niveles del agua en la seccin de aforo, lo cual nos permitir desarrollar una
curva de calibracin, es decir, una relacin entre el nivel y la descarga.
82
F_t
-l!.'t'o*
D
'i
,r
IF.
a.|
"-l
'/
Seccin Transversal Tpic
it \-
ca uda
bl Curva de Calibracin
Figura Nql2: Determinacin de1 Caudal.
Ejemplo
12:
el
Determinacin del Caudal
Obtener
caudal total
y la velocidad
promedio en
ia
seccin
transversal de un ro, utilizando la informacin que se muestra en la siguiente
tabla, obtenida por medio de un aforo por vadeo.
B3
Tabla: Datos del Aforo

Estcin (m)
Profundidad: D (m)
o_28
Lectura *
(m)
Velocidad: V
(m/seg)
3.00
3.75
0.000 0.13 0.37 0.09
0.32 0.46
0.329
0.448 o.397 0.430
4.50
5.00
0.56
0.45
0.
11
o.310
0.468
6.50
0.69
0.55 0.14
0.443
0.527
7.00
0.75
0.60 0.15
0.428
0.458 0.357
7.75
0.73
0.58 0.15
8.75
o.62
0.50
0.428
0.3s2 0.458
0.r2
i0.00
o.s2 0.42 0.10
10.50 o.41
0.329 0.448 0.302

0.410 0.000
0.33 0.08
12.75 14.00
0.33 0.00
0. 13
* Nota de Aclaracin: l lectura en la varilla
de vadeo para la medicin de la velocidad
est tomada con respecto al fondo. Cuando D < 0.40 m se efectu una sola medicin a 0.40 D; pero si D > 0.40 m, se efectuaron dos mediciones a 0.80 D 0.20 D; obsewe que en este caso la suma de las dos lecturas es igual a la profundidad.
a4
Tabia: Clculo del Aforo

Esucin (m)
3.00 3.7s 4.50 5.00 6.50 7.00
'7.75
Profundidad D (m) 0.28 0.32 0.46 0.s6 0.69 0.15

0.13
Velocidad V, (m/s) 0.1097 0.3290 0.422s 0.4000

0.4555 0.4775 0.4075
bi
Ancho (m) 0.315 0.750 0.62s

1.000 1.000 o.625 0.875 1.125 0.875
Ai
Area (m')
Caudal Q (m3/seg) 0.0115
0.1050 0.2400 0.2875 0.5600 0.6900 0.4688 0.6388
0.0790
0.t2r5
0.2240 0.3143
0.2239 0.2603
8.15 10.00 10.50
0.62 0.52
0.41 0.33
0.3900
0.3935
0.697s
0.45s0 0.5638 0.57'/5 0.0516
o.2720 0.1790 0.2114

0.2368
0.3750 0.4100
0.1367
t.375
t.'750 0.625
t2.15
14.00
0.00
0.0070
Suma:
Observaciones:
Ib; :11.00 IAi:5.3355 IQi :2.1407
a)
Columna 3: La velocidad media es el valor nico rnedido el valor promedio de las dos velocidades observadas en cada estacin. En ambas orillas del ro (estaciones de los extremos) se utiliza 1/3 V de la velocidad registrada en las
estaciones adyacentes.
b)
Colurna 4: El ancho correspondiente a cada estacin es la diferencia entre las

distancias a los puntos medios entre cada estacin, con excepcin de los extremos
donde se utiliza el punto rredio entre la estacin adyacente y la estacin de la orilla. Colno regla general tenemos que:
b = r/2 (ESTti total de la seccin transversal del ro).
EST:
Adernas, observe que la suma de todos los anchos debe ser igual al espejo (ancho
85
c)
Columna 5: El rea de cada seccin es igual a la profundidad multiplicada por el ancho; es decir:
Ai = b.' Di
Excepto en las orillas sin profundidad donde se utiliza un cuarto (1/4) de la
d)
profundidad de la estacin adyacente. Columna 6: El caudal en la secciu es el producto de la velocidad media y el rea; por consiguieffe:
o = v..a
Solucin:
a)
Ancho total B (espeio) de la seccin transversal del ro:
_ EST B = E.S/ tnt.ral Jlnat

B = 14.00 - 3.00 =
Ib
17.0O m
b)
Area total
Ar:
Ar = LA,- = c)
Caudal total
5.3355m2
Qr:
Q' = L Q' = d)
21407
m3
/seg
Velocidad Promedio en la Seccin Transversal:
v = 8' Ar
2 7407 m3
lseg =
o.4or2 mrseg
5.3355 m2
ANALISIS DE PRECIPITACION - ESCORRENTIA

RELACION ENTRE PRECIPITACION Y ESCORRENTIA
Cuando la precipitacin sobrepasa la rata de infiltracin en la superficie
del terreno, el exceso de agua comienza a acumularse como almacenamiento superficial enpequeas depresiones del terreno originadas por la topografa.;
eventualmente el flujo escurre sobre la superficie del terreno en algunas porciones de la cuenca
y el flujo
se concentra rpidamente en arroyos
canales pequeos, los cuales fluyen a su vez hacia corrientes ms grandes
ros. Tal como
se mencion anteriormente, el
flujo subsuperficial y el flujo
base tambin contribuyen en cierta proporcin con el hidrograma total de

descarga durante un evento de lluvia.
Los hidrlogos e ingenieros civiles especializados en

cuenca
recursos
hidrulicos, estn interesados en la cantidad de escorrentia generada en una
por un patrn dado de precipitacin. Se han efectuado muchos
intentos para analizar estadsticamente los datos histricos de precipitacin,
evaporacin
datos de
flujo en corrientes con la finalidad de desarrollar
relaciones de prediccin entre estos procesos. Factores tales como

precipitacinprecedente, humedad del suelo, inltracin variable y respuestas
diferentes de la escorrentia superficial con las estaciones del ao, hacen que
el desarrollo de estas relaciones sea muy difcil.

Una gran cantidad de investigadores ha intentado desarrollar relaciones
de precipitacin
- escorrenta
que se puedan aplicar a cualquier regin
cuenca bajo cualquier serie de condiciones. Sin embargo, estos mtodos
87
deben ser utilizados con extrema precaucin debido a la variabilidad de los
factores que afectan la evaluacin de la escorrenta a partir de un volumen

conocido de precipitacin.
Las relaciones simples de precipitacin
escorrenta deben ser
utilizadas solamente en estudios de planificacin de recursos hidrulicos

cuando se requiere una estimacin rstica de la respuesta de la cuenca. Es
importante sealar que se requiere un conocimiento detallado de la magnitud
y distribucin en el espacio (rea) y tiempo de la preeipitacin y escorrenta
para el anlisis completo de los proyectos de control de inundaciones y

estudios de planicies de inundacin, especialmente en las regiones afectadas
por el drenaje urbano. Una de las frmulas mas simples de precipitacin
escorrenta se
denomina Mtodo Racional, el cual permite la prediccin del caudal mximo
de la siguiente manera:
o = c'i.A
donde:
C
coeficiente de escorrenta, el cual vara con el uso de la tierra
intensidad de la precipitacin para la frecuencia
perodo de retorno seleccionado

igual al tiempo de concentracin /"
tc
una duracin
tiempo que demora la lluvia que cae en
el
punto
mas remoto de la cuenca en viajar hasta la salida

A
rea de la cuenca
8B
Existen muchos mtodos empricos para estimar la magnitud del tiempo

de concentracin; una de las primeras frmulas fue propuesta por Johnston
y Cross
(9a\
y es la siguiente:
t'
principal
= orl
/"
''l' lrl-r)
est en horas, la longitud del canal
Donde el tiempo de concentracin
en millas y la pendiente promedio del canal ,S en pies/milla.
El Mtodo Racional se le atribuye generalmente a Mulvaney (185i),

quien describi el procedimiento en una publicacin tcnica en Irlanda. El
mtodo est basado en
la
suposicin de que una rata de precipitacin
consfante y uniforme producir
la esconentia mxxna cuando fodas las partes
de la cuenca estan contribuyendo con el caudal; observe que esta condicin se satisface cuando la duracin de la lluvia es mayor igual al tiempo de
concentracin. El Mtodo Racional fue el precursor del concepto
de
hidrogramas de tormentas y el resto de este captulo se dedicar al desarrollo de la teora de hidrogramas y a la aplicacin de estos mtodos al anlisis de
lluvias complejas en cuencas grandes.
Ejemplo 13:
Mtodo Racional
Se desea determinar el caudal
mximo
tormenta con un perodo de retorno
("n m3/s) producido por una
10 aos. El rea de la cuenca es
A:
40 hectreas y el tiempo de concentracin
es
t" :
25 minutos.
B9
Suponga que la curva de intensidad-duracin-frecuencia para este perodo de
retorno es la siguiente:
I 5
t" *
donde
36
esta en plg/hr
y /" en minutos. Los datos sobre el uso de la tierra

se muestran en la siguiente tabla.
y coeficientes de escorrenta
Tabla de Datos
Uso
Area
(ha)
Coeficiente de escorrenta 0.40 0.60 0.15
Residencial
30
3
Comercial
Parques
Solucin:
a)
Coeficiente de escorrenta:
E c,' A,
E
,t,
0.40. 30 + 0.60' 3 + O.15 . 30 + 3 + 7
c -
1485 =
40
0.37
90
b)
Intensidad:
323
tc+36
323
25+36
5.3
plglhr
c)
Caudal mximo:
o = c-i.A
Qo = 0.37
(r' ,"'r,
0.0254 m
lp lg
.
40
ha.
I o q 00m
I ha
)(
hr
Qo
= (,,r,nrrru m3 lhr
5 .5 J m3
J, 600 s
ls c
ANALISIS DE HIDROGRAMAS
El proceso de escorrenta superficial es el resultado de una combinacin de condiciones fisiogrficas y meteorolgicas de la cuenca y representa los
efectos combinados de la precipitacin, prdidas hidrolgicas,
flujo sobre la
superficie del terreno, flujo subsuperficial
y flujo subterrneo. Segn
Sherman (1932), los factores climticos que influyen en la forma del

hidrograma y volumen de escorrenta son:
1. 2. 3.
Patrn e intensidad de la lluvia
Distribucin de la precipitacin sobre la cuenca

Duracin de la tormenta
7-
9't
y los factores fisiogrficos de mayor importancia son:

1.
Tamao y forma del rea de drenaje Naturaleza del sistema de drenaje

Pendiente del terreno y del canal principal
2. J. 4.
El almacenamiento por retencin en la cuenca
Durante un evento de precipitacin, las prdidas hidrolgicas tales
como infiltracin, almacenamiento en depresiones ] almacenamiento por

retencin deben ser satisfechas primero antes que comience la escorrenta
superficial. Al incrementarse la profundidad del agua retenida en la

superficie, el flujo superficial sobre el terreno comenzar en algunas
porciones de la cuenca. El agua eventualmente se mover hacia los pequeos
arroyos, quebradas y ros que constituyen el sistema de drenaje de la cuenca.

Parte del agua que se infiltra en el suelo se mueve lateralmente a travs de las
capas superiores del suelo hasta que alcarua una corriente
superficial; a sta
porcin de la escorrenta se le denomina flujo subsuperficial. Otra porcin
de la precipitacin que se filtra en el suelo alcarua el nivel fretico,

generalmente varios metros debajo de la superficie del suelo y contribuir
a
la escorrenta como flujo base si el nivel fretico intersecta la superficie de la
corriente de agua en el canal.
92
FFECIPITACION ESCOHFENTIA PRECIPITACIEN DE INTENSIDAD UNIFOFME INFILTRACII]N ALMACENAI"iIENT0 DE DEPFESINES ALMACENAMIENTO

POR RETENCION
ESCOHRENTIA
PF!DUEIDA PIF
ESCOBFENTIA DIFECTA
EL ALMACENAMIENTO
FETENIDO TIEMPo
a) Distribucion de una Precipi'ta.in Unilorrn".

LLUVIA PUNTO DE
DES|}EGA
INFLEXION
b) Hidrograma de Equilibrio.
TIEMPO
Figura N"l-3:
Fenmeno
de Escorrenta Superficial
se
En la parte (a) de la figura del fenmeno de escorrenta directa
ilustra la distribucin de la precipitacin uniforme para una duracin finita.

Pero si la precipitacin contina a intensidad constante por un tiempo infinito.
entonces se produce una descarga de equilibrio, es decir que el flujo de

entrada es igual al de salida y esto se logra cuando el rea total de la cuenca

es contribuyendo al
93
Esta condicin de equilibrio se observa en muy
flujo.
raras ocasiones en la nafural eza, excepto en cuencas urbanas muy pequeas,
debido a las variaciones en la intensidad y duracin de la lluvia.
El flujo base en un canal natural se debe a las contribuciones del agua

subterrnea poco profunda y es otro de los componentes del
hidrograma. En
cuencas naturales grandes, el flujo base representa una fraccin significativa
de la descarga total, mientras que se puede despreciar en cuencas pequeas
urbanizadas donde el flujo superficial sobre el terreno es predominante.
El flujo base es separado y sustrado del hidrograma total para obtener

el hidrograma de escorrenta directa. El volumen de agua del hidrograma de escorrenta directa debe corresponder a la precipitacin excedente neta, la cual se obtiene sustrayendo de la precipitacin total la infiltracin y las otras
prdidas de agua (almacenamiento en depresiones, intercepcin, etc.). El hidrograma de escorrenta directa representa la respuesta de la cuenca a la
lluvia excedente, mientras que la forma y los tiempos del hidrograma

escorrenta directa estn relacionados con la intensidad
de
y duracin de la
precipitacin, as como tambin a los factores fisiogrficos que afectan el

almacenamiento.
El hidrograma en la mayora de los casos esta constitudo por un

miembro ascendente, segmento de la cresta y una porcin descendente
recesin. La pendiente del miembro ascendente esta determinada en gran

parte por la intensidad de la tormenta y el punto de inflexin en el segmento
descendente generalmente indica el instante en que termina la escorrenta
94
superficial directa. En la recesin, la descarga se produce por el movimiento

de agua subterrnea almacenada en el subsuelo que fluye hacia la superficie.
II{TEIISIDAD
DEIIUJIA
?EECI}ITAEIOT{ ]IETA
VI]LIJI.fEIf NX ESI]OEEEI{TIA DISETA
CIESIA
sxGl.{EttTo
SCEHDE TT
SXG!4EI{TO DESCI'DEIITE
ESI]OEEXIITIA
DIXXCIA
DESCAEGA
TII{ IlE XSCOENEMIA DINXMA
,,'
FLUJO BASX
xEcEsIolf
II{ICIO DE ESCOXEEIITIA
DIXECTA
Figura N"14: Hidrograma de Descarga Total
COMPONENTES DEL HIDROGRAMA

Un hidrograma esta formado por varios componentes: flujo superficial
sobre el terreno, flujo subterrneo
flujo base y flujo subsuperficial,
producido por el agua infiltrada que se almacena temporalnente en las capas superiores del suelo
y entra posteriormente de manera lateral al canal. La
contribucin relativa de cada componente al hidrograma depende de la

intensidad de la precipitacin i medida con respecto a la rata de infiltracin
95
/,
as como tambin el almacenamiento de humedad en el suelo en relacin
con la capacidad de campo; la cual definimos anteriormente como la cantidad

de agua retenida en el suelo despus que la gravedad drena el exceso de agua.
Observe que no ocurrir escorrenta superficial en el caso de que la

intensidad de la precipitacin sea menor que la rata de infiltracin
(i < f );
al igual que el flujo subsuperficial y el flujo subterrneo sern nulos si la

humedad almacenada en el suelo es menor que la capacidad de campo, ya que todava existe almacenamiento adicional de humedad en el suelo, a menos que
el flujo subterrneo se produzca por almacenamientos de agua a largo plazo. En la situacin de que tanto la precipitacin sea mayor que la infiltracin, as
como la humedad del suelo sea mayor que la capacidad de campo (lo cual es
tpico en casos de eventos de tormentas grandes), la escorrenta directa
superficial, el flujo subsuperficial
y el flujo base flujo subterrneo

La precipitacin
en el canal
contribuirn con el hidrograma de descarga total.
es un componente tambin, pero generalmente es una
fraccin muy pequea
de la descarga total.
La escorrenta superficial es un factor muy importante cada vez que

llueve considerablemente. Las precipitaciones de gran intensidad, as como
tambin
la
urbanizacin
deforestacin
de las cuencas, simplemente
magnif,ca el caudal mximo y disminuye el tiempo para alcanzar la descarga
mxima. El flujo subsuperficial puede ser un factor grande en tormentas de

moderada intensidad sobre cuencas con capas delgadas de suelo sobre roca
impermeable. Si el agua subterrnea fluye hacia la corriente durante perodos de lluvia fuerte, la corriente se denomina efluente emanacin; pero si el
96
flujo de la corriente es hacia el sistema de agua subterrnea, como en el caso

de condiciones de sequa, la corriente es llamada influjo o afluencia.
En la prctica, es muy comn considerar que ia descarga total se divida
en dos partes solamente: escorrenta directa
y flujo base. La escorrenta
directa puede incluir una porcin considerable del flujo subsuperf,rcial,

mientras que el flujo base esta constitudo principalrnente por el flujo de agua subterrnea.
DESCAFGA O
I ,(
I
nr ot L.l LLuvn
DIRECTA PRECIPIT E EL SUPERfICIAL

I
ILUJO SURSPEBFICIAL
I I
lltlfPtl
Figura N"15: Componentes del Hidrograme
9l
SEPARACION DEL FLUJO BASE Y RECESION

Existen muchas tcnicas para separar la escorrenta directa del flujo
base, las cuales estn basadas en el anlisis de las curvas de recesin del agua
subterrnea. En algunos casos, la curva de recesin se puede describir por

medio de una ecuacin exponencial de la siguiente forma general:
llt = qo' " '

donde: descarga inicial especificada descarga en un tiempo posterior t constante de recesin
base de los logartmos naturales
Qo
Qt
k
e
Las ecuaciones de esta forma son utilizadas a menudo en ingeniera
para describir un agotamiento de primer orden

transformar a una lnea recta: log,o
y la misma se puede
(q,) =
logro
(4,) - | k.log,o(e)l /
Por
consiguiente,
al
graficar
la
ecuacin
de
recesin
en
papel
semilogartrnico, obtendremos una lnea recta; la diferencia entre esta lnea
y la curva correspondiente al hidrograma de descarga total graficado en el

mismo papel, representra el hidrograma de escorrenta directa. En la prctica
se
utilizan otros mtodos para separar el flujo base. El mtodo mas simple
consiste en dibujar una lnea horizontal desde el punto donde comienza la
9B
escorrenta superficial hasta intersectar el hidrograma de recesin. Otro

mtodo sugerido por Linsley, Kohler y Paulus (1949), el cual es utilizado
frecuentemente, consiste en extender la curva de recesin anterior al inicio de
la lluvia
hasta un punto situado debajo de Ia descarga mxima
luego
conectar este punto con el hidrograma en el lugar donde termina la escorrenta
superficial, el cual se puede estimar utilizando la siguiente frmula:
i/ = e Ao2
donde:
tiempo en das despus de la descarga mxima

rea de drenaje de la cuenca
factor de conversin que depende de las unidades

de A (millas cuadradas kilmetros cuadrados)
1.0 para A en mi
0.8 para A en Km2
Todos estos mtodos tienen la desventaja de que son arbitrarios y en

cierta medida inexactos. La separacin del flujo base es un arte ms que una
ciencia; en muchos casos de inters prctico tales como drenaje urbano, el

flujo base se desprecia a menudo porque representa una fraccin muy pequea
de la descarga total. Generalmente, el flujo base es mas importante
en
corrientes naturales y ros grandes debido a la contribucin a la descarga total
a lo largo de las orillas del canal o riberas del ro proveniente del flujo de
agua subterrnea. Independientemente del mtodo seleccionado para la

separacin del flujo base, el hidrlogo debe ser consistente con el mtodo
99
utilizado, de manera que los hidrogramas de una tormenta se puedan

comparar con los otros, al igual que los hidrogramas de escorrenta directa
de diferentes cuencas.
OESCARGA Q
FIN DE
IA ESCOBFENTIA
DIBECTA
LINEA FECTA
HOFUO NTAL
TIEMP,O
a
Figmra N"16: Separacin del Flujo

Base.
PRECIPITACION NETA Y EL HIDROGRAMA

Anteriormente analizamos el fenmeno de escorrenta superficial y
distribuimos
la precipitacin total en varias
componentes: infiltracin,
almacenamiento en depresiones, almacenamiento por retencin y escorrenta
directa. Podemos escribir una ecuacin hidrolgica de continuidad de

siguiente forma:
la
100
Precipitacin
Total
Almacenamiento en
Depresiones
+ Eyaporacin
Infiltracin
Escorrenta
SuperFrcial
Donde hemos asumido que el almacenamiento por retencin se convertir en

escorrenta superficial despus que termira la
lluvia. A menudo
es importante
determirnr la distribucin con respecto al tiempo del exceso de precipitacin

o precipitacin neta, la cual es igual volumen total de escorrenta directa; es
decir el volumen de escorrenta directa superficial ms volumen almacenado retenido, el cual escurrir de la cuenca durante un int-ervalo de tiempo mayor que la duracin de la tofinenta. En otras palabras:
Precipitacin Excdente Neta
Precipitacin Total - Prdidas Hidrolgicas
Precipitacin Excedente Neta
Escorrenta Directa Superficial
Generalmente, los mtodos empleados para determinar la precipitacin excedente incluyen el mtodo de Horton, incluyendo una prdida inicial por almacenamiento en depresiones, y el mtodo del ndice Q, el cual consiste en
una prdida constante durante el perodo de tormenta. En la prctica, los

coeficientes de prdida por infiltracin son difciles de estimar, por lo tanto,
el
procedimiento simple
del ndice S es el mtodo utilizado
ms
frecuentemente debido a la escasez de datos para la distribucin de infiltracin
en el tiempo. Observe que el ndice tiende a subestimar las prdidas al
inicio de la tormenta y sobreestimar las prdidas al final.

Tan pronto como la precipitacin neta ha sido determinada para una
cuenca, entonces el siguiente paso para convertir la precipitacin neta en
101
escorrenta directa superficial es un problema fundamental de la ciencia de la
hidrologa.
El
hidrograma de descarga resultante se forma
por
las
contribuciones del flujo superficial sobre el terreno y del flujo en el canal
llegando en diferentes tiempos desde todos los puntos de la cuenca. Los

tiempos relativos de viaje para el flujo superficial y el flujo en el canal estn relacionados con el tamao de la cuenca; el tiempo de flujo superhcial es ms
significativo para cuencas pequeas, mientras que el tiempo de viaje en el

canal predomina en cuencas grandes.
PRClPlTAC!{]!' I l' FttT R.t G t
OJRVA DE
IN FILf F.Al N DE HOFION
TEtrlPO (
Figura N'17: Distribucin de la Precj-pitacin Total
L02
METODO DE TIEMPO - AREA

Una manera interesante de comprender el proceso de conversin de la
precipitacin excedente en un hidrograma consiste en utilizar el concepto de
Histograma de Tiempo - Area desarrollado por Clark (1945). Este mtodo

asume que el hidrograma de descarga se produce por la traslacin pura de la
escorrenta directa hacia la salida, ignorando los efectos de almacenamiento

en la cuenca. Si urn precipitacin de intensidad uniforme se distribuye sobre
la cuenca, el agua comienza a fluir primero de lai reas inmediatamente

adyacentes a
la salida y el porcentaje del rea total contribuyendo a la
descarga se incrementa progresivamente en el tiempo.
S:roNA
PFECIPfTACIOi!
SAUA
PI
P{
P5
TE'PO
a).
CUENCA
b} . HDROGFAMA
DE
PRECIPTTACIN i{ETA
c)
. HISTOGFAMA 0E TIEMPO-AREA
d)
. HIDROGR,{MA DE DESCARGA DIRECTA
Figura N"18:
Mtodo de Tiempo-Area
103
Por ejemplo, en la figura anterior, la escorrenta superficial proveniente
del rea
A,
alcanza la salida primero, seguida de las contribuciones de Ar,
luego A, y finalmente
d.
Este proceso para la generacin de los caudales se
puede expresar de la siguiente manera:
Q, - Pi Ar+P.-r'Ar+
donde:
+Pr-l+ .'A: L l
Vi<i
i
Qi
Pi
ordenada del hidrograma en el tiempo
precipitacin excedente en el tiempo i

ordenada del histograma tiempo-rea en el instante j
Observe que el nmero de ordenadas del hietograma no tiene que ser necesariamente igual al nmero de ordenadas del histograma. Por ejemplo,
esta frmula nos indica que la precipitacin neta del perodo
Pr, sobre el rea

para
4., P, sobre A, y P, sobre A, arriban a la salida simultneamente

producir la descarga
Q.. El hidrograma
de descarga directa se desarrolla
evaluando todos los caudales Qr, Qz, Qs,...,Q" donde Qn
0 indica el fin de
la escorrenta directa.
El
concepto de tiempo
- rea proporciona una idea bsica para
comprender el fenmeno de escorrenta, pero sus aplicaciones son limitadas debido a que el hidrograma debe ser ajustado posteriormente para incluir los
efectos de almacenamiento en la cuenca;
adems, las lneas isocronas
(lneas de igual tiempo de viaje hasta la salida) son muy difciles de construir.
104
Un concepto ms general en la prctica es la teora del hidrograma
unitario,
la cual es reconocida como una de las contribuciones mas
importntes relacioradas con la prediccin de Ia escorrenta superficial. Esta
teora, combinada con los mtodos de infiltracin y el anlisis del trnsito de avenidas crecidas en canales y embalses, es suficiente para investigar los
efectos de la precipitacin y almacenamiento en cuencas pequeas y grandes.
Observe que el mtodo de tiempo - rea es un caso especial del hidrograma
unitario.
Ejemplo 14:
Mtodo de Tiempo - Area
Una cuenca se divide en cuatro secciones A, B, C y D. Los datos de

rea y el tiempo de viaje desde cada regin hasta la salida se muestran en la
siguiente tabla;
Seccin
A
100
1
D
100
rea (acres) Tiempo hasta la estacin G (hr)
200
2
300
3
La
escorrenta desde cada seccin contribuir con la descarga directa
en la estacin localzada en la salida (punto G). Considere que una

precipitacin neta de 0.5 plg/hr de intensidad est cayendo uniformemente sobre toda el rea de la cuenca durante un intervalo de tiempo de 5 hr.
Encuentre el hidrograma de descarga directa en la estacin G en respuesta a
la lluvia y asuma que no hay prdidas por almacenamiento.

Solucin:
105
En ste caso, la frmula para obtener el caudal directo usando el

Mtodo de Tiempo - Area es la siguiente:
Q, = P,'At * P,-r'A, * P,-r'A, * P,-r'Ao

La frmula anterior, nos permite obtener los caudales (escorrenta directa)
para diferentes instantes:
eoo
Qt Qz Q. Qn Q, Qo Q, Qr
P,'A,
Pr.A, + Pr.A2 Pr.A, + pz.A2 + p1 .A3 Po.A, * Pr. A2 + P2.A' * P,.Ao Pr.A, * Po.A2 + P3.A' * Pr. Ao Pr.A, * Po.A3 + P3.A4 Pr.A, + P4.A4
Pr'Ao
eno
Observacin: La conversin de unidades en la multiplicacin
(e : p . A) es:
l4ll
,\
"rg i
=[ \ n,
ete.o",nl
\ I
/ \
1.0083
.l
hr
r,eoo
r"g t aue
*
.qz,seopes2. t
pie
\
I
t2 pts
pie3/seg
acre.plglhora
1 acre. plg/hr

Tabla de Clculos
Tiempo (horas)
0
1
106
Pn t"
Area
(acres)
P''{
Pr4 Pr.4 Po'4 P''{
(plg/hr)
Q'
(ptls)
0
0.5
0_5
100
50
100 150
50 50
100 150 150
2 3
200
300
100
0.5 0.5 0.5
50
100 150
300
50
100
4
5
50
350
50
100 150
50
350 300
6
7 8
50
150 50
200
50 0
50
o
(rictlrcgt
350 300 250 200
t50
100
50
Grfica de1 Hidrograma de Escorenta Drecta

Alternativa:
L07
El sistema de ecuaciones lineales simultneas descrito anteriormente

puede resolver utilizando lgebra lineal de la siguiente manera:
se
0 = tpt
o o Y.2 o -3
Q1 Q5
P1
000
P2 P3 P4
Pro
P2
o 0
Pr
ar
a2
A3 A4
P3 P2 Pr P+ P3
P2 P3 P4
P0 0 0
o,
o-
Ps P4
P5
8,
0 0P5
Por consiguiente, las ordenadas del hidrograma de escorrenta directa

(solucin del sistema) son:
g1
Q2
%0 00 Yz Yz 00
Y20 %%
100
50 150
ogn
o_
300 350 350

3
200
300
100
oo_
g8
0% 00 ll r/ 00 0%
00
200
50
108
TEORIA DEL HIDROGRAMA TINITARIO

Sherman (1932) introdujo la teora del hidrograma unitario (HU), el
cual se define como el caudal de descarga de una cuenca que se produce por
una pulgada de escorrenta directa generada uniformemente sobre el rea de
drenaje por una precipitacin de intensidad constante durante el perodo de
duracin de la lluvia (en el sistema mtrico generalmente se adopta una

profundidad de un centmetro).
Muchas de las suposiciones inherentes del hidr-grama unitario tienden
a limitar sus aplicaciones en una cuenca de acuerdo a Johnstone y Cross
094\;
entre ellas tenemos:
1.
Se asume que la precipitacin excedente de igual duracin producir hidrogramas con tiempos bases equivalentes, independientemente de la
intensidad de la lluvia.
2.
Se asume que las ordenadas de escorrenta directa para una tormenta

de duracin dada son directamente proporcionales a los volmenes de
la precipitacin neta.
3. 4.
Se asume que el tiempo de distribucin de escorrenta directa

independiente de la precipitacin antecedente.
es
Adems tambin se asume que la distribucin de la precipitacin sobre
la cuenca es la misma para todas las tormentas de igual duracin.
La proposicin clsica de la teora del hidrograma unitario
se puede
resumir brevemente: los sistemas hidrolgicos son lineales e invariables con el tiempo (Dooge, 1973). La propiedad de proporcionalidad y el principio de
109
superposicin se aplican al hidrograma unitario, pero estas suposiciones no
fueron cuestionadas seriamente hasta la mitad de la dcada de 1950.
Las superposiciones de linearidad e invariancia con el tiempo no son

estrictamente correctas para una cuenca, pero las adoptamos frecuentemente
cuando nos son tiles; puesto que existen ejemplos no lineales de flujo en
canales abiertos, modelos de laboratorio de escorrenta y demostraciones en
el campo. Sin embargo, las suposiciones lineales todava se utilizan porque
son relativamente simples
son la base de los mtodos disponibles ms
desarrollados; adems, los resultados obtenidos con mtodos lineales son

aceptables para la mayora de los propsitos requeridos en ingeniera.
DERIVACION DE HIDROGRAMAS UNITARIOS Las cuencas con estaciones de medicin nos permiten registrar el
hietograma de precipitacin y el hidrograma de descarga producido por la tormenta sobre la regin. El hietograma de precipitacin es una grfica de
la intensidad de la lluvia en funcin del tiempo y el hidrograma es una grfica
del caudal de descarga en funcin del tiempo. En la siguiente figura,

muestra el desarrollo del hidrograma unitario parmetros relacionados con
se
tambin se ilustran los
el tiempo, tales como la duracin de la
precipitacin neta o excedente y el tiempo hasta la cspide valor mximo de

la descarga. El hidrograma de descarga total esta constitudo por un miembro
ascendente, segmento de la cresia, segmento descendente y curva de recesin.
110
Los aspectos de tiempo del hidrograma estn caractefizados por los

siguientes parmetros:
1.
Tiempo de retraso (tr): es el tiempo desde el centro de masa de la precipitacin excedente o neta hasta la cspide del hidrograma de
descarga.
2. 3.
Tiempo de ascenso
(t):
es el tiempo desde el
inicio de la precipitacin
excedente hasta la cspide del hidrograma.
Tiempo de concentracin (t"): es el tiempo de eqtrilibrio para la cuenca,

es decir cuando la descarga es igual a la entrada. Tambin representa
el tiempo de viaje de una onda para propagarse desde el punto

distante de la cuenca hasta la salida.
ms
4.
Tiempo base (t"): duracin del hidrograma de escorrenta directa.
Las siguientes reglas generales se deben observar en el desarrollo de

hidrogramas unitarios para cuencas con estaciones de medicin:
1
Las tormentas seleccionadas deben tener una estructura simple con

distribuciones relativamente uniformes en el espacio y tiempo.
2. 3. 4.
El tamao de la cuenca generalmente debe estar dentro de 1,000 acres a 1,000 millas cuadradas.
El volumen de escorrenta directa debe encontrarse en el rango de 0.5

pulgadas a 2 pulgadas.
La duracin de la precipitacin excedente debe ser aproximadamente
del25% al30% del tiempo de retraso.

5.
Se debe analizar cierto nmero de tormentas de duracin
111_
similar para
obtener un hidrograma unitario promedio para esa duracin de lluvia

neta efectiva. 6.
El paso anterior se debe repetir para precipitaciones

diferentes duraciones.
excedentes de
Los siguientes pasos son fundamentales para el desarrollo de un

hidrograma unitario a partir de un evento simple de precipitacin:
l. 2. 3.
Analizar el hidrograma y separar el flujo base.
Medir el volumen total de escorrenta directa y convertirlo a pulgadas

centmetros de profundidad sobre la cuenca.
Convertir la precipitacin total a precipitacin neta
excedente
utilizando algn mtodo tpico de infiltracin y obtener la duracin D

de la lluvia para identificar el tiempo para el hidrograma unitario
y el
hidrograma de escorrenta directa.
4.
Dividir las ordenadas del hidrograma de escorrenta directa por el

volumen en pulgadas centmetros de la lluvia neta y graficar los
resultados para obtener el hidrograma unitario de la cuenca. Se asume
que el tiempo base
t"
es constante para tormentas de igual duracin.
5.
Verificar el volumen del hidrograma unitario para estar seguro de que

corresponde a una pulgada un centmetro
hid rograma en caso necesario.
y ajustar los valores del
1-72
LLt'\TA fplg/i'rl
Rf10s
t0
cr[Ar ETa
u o
[p'lsul I so
129
80
{0
a.
Hielogrm dc Llua I Hidrograma de Dcscarga Total
vd.hn de h^/b nd (i -0)
D=R=
49
100 p
ls
R=2 PLGDE
ESCOhRE}{NA
DI*CTA
I
I t ; ; hr i
tB-- 11
b. Precipitacin
ioir,o'l
Excedente I Hidrograma dc Ocscrrga 0irecta
Hidrologa
/ David Cedeo
113
vofunen unrl*io de
lwia
plg
-asnp3h
R= 1 PLG DE
EScOhRENTIA
DIRECTA
I
Hidrograma Unilario para una duracin D = 2 horas
Figura No19: Desarrollo del Hidrograma Unitario
Debido a las suposiciones de linearidad esenciales en el desarrollo del hidrograma unitario, es necesario resaltar que debemos tener mucho cuidado
al aplicar el mtodo del hidrograma unitario cuando las condiciones existentes
violan el principio de superposicin. Grandes variaciones en la duracin e

intensidad de
la precipitacin se reflejan en la forma de hidrograma
resultante. A menudo se suman hidrogramas unitarios de corta duracin,

retrasados el uno con respecto al otro, para generar un hidrograma de una duracin mayor. Si la variacin de la intensidad de la precipitacin es muy
grande durante la duracin de la tormenta, la suposicin de linearidad puede que no siga siendo vlida en sta situacin.
1-1-4
El volumen de escorrenta directa y la distribucin sobre el rea de la

cuenca de la escorrenta superficial pueden causar variaciones en la forma del
hidrograma. Para cuencas grandes, generalmente es mejor desarrollar

hidrogramas unitarios para sub-cuencas, luego se suman y tambin se atrasan
de acuerdo al tiempo de viaje de la onda para generar el hidrograma de la

cuenca a
total. Los caudales mximos
en los hidrogramas unitarios obtenidos
partir de eventos de tormentas muy pequeas, a menudo son de magnihrd
inferior que aquellos obtenidos utilizando tormentas nxs grandes debido a las
diferencias en el flujo subsuperficial y a los tiempos de flujo en el canal. Las
suposiciones de linearidad se pueden verificar generalmente comparando

hidrogramas para tormentas de diferentes magnifudes. Si no existe linearidad,
entonces los hidrogramas unitarios derivados deben ser utilizados nicamente
para generar eventos de magnitudes similares y se debe tener mucho cuidado
al ttllizar los hidrogramas unitarios para extrapolar eventos extremos.
pesar de estas limitaciones, cuando se utilizan los hidrogramas
unitarios en combinacin con los mtodos de trnsito de avenidas crecidas
para la prediccin de los niveles de inundacin en cuencas pequeas y

grandes, la exactitud de los resultados se ha incrementado en la ltima dcada
con el uso de programas de computadora.
Ejemplo 15:
Hidrograma Unitario
Obtener el hidrograma unitario para una cuenca con un rea 1,135 .5 acres, donde se ha medido el siguiente hidrograma de descarga total y el
hietograma de precipitacin que origin la respuesta de la cuenca:
115
/F[H#]
t5
2 6
a.
Hietograma de lluvia
10
b. Hidrograma de Descarga
Asuma que el ndice de infiltracin 0 para sta tormenta es de 0.5

pulgadas/hora y que el flujo base es constante y tiene una magnifud de 100
pies3/segundo. Adems, indicar el tiempo base tu, el tiempo de atraso to y la
duracin D de la lluvia neta.
Solucin:
a)
Precipitacin neta:
R = (t-0)A/
neta
= (1.s-0.s)(3-1) =
2ptg
El resultado anterior representa una precipitacin con una intensidad
(i - 0)
1 plg/hr durante dos horas; por 10 tanto la duracin de la
precipitacin neta es:

D
r-
2 plg
I plglhr
= 2 horas
Hidrologa
/ David Cedeo
116
Adems, observe que ra interuidad de la precipitacin neta que generar
el hidrograma unitario debe transformarse de la siguienrc manera:
I pls
D
/2 plg/ hr
para que el volumen de escorrenta directa sea equivalenre a una unidad.
fpls/t'rl
Prc cip ita
cin Nela
= 0.5 plslhr
Prdidas por lnfillracin

r th4
b)
Separacin del flujo base
La ordenada del hidrograma de escorrenta directa eo se obtiene

substrayendo el flujo base eo de la escorrenta total
e,,
es decir:
Qo = Q,-Q,
Lt7
c)
Hidrograma Unitario
Las ordenadas del hidrograma de escorrenta directa se dividen por
1a
precipitacin neta para obtener el hidrograma unitario de la cuenca; por lo

tanto, podemos establecer la siguiente proporcin:
o
rpls
o.
R
Tabla de Datos y Clculos

Tiempo
Caudal Totl
Q, (p3lseg)
t
(hr)
Flujo Base Q,
(p3lseg)
Escorrenta Direct
Hidrograma
Qo:Q-Qo
(p3lseg)
Unitario Q,
(p3lseg)
0
1
100 100
100 100 100
0
0
0
100
2
3
300 700 1,000 800 600 400 300
200
600 900 700
100 100 100 100

100 100 100 100 100
300 450 350
4
5
s00
300
2s0
150 100
7
8
200
100
9
10
11
200
100 100
50
0
0 0
L Q,
1,750 P'ls
118
Finalmente se dibuja el hidrograma unitario para la cuenca

determinan los tiempos requeridos:
y se
Tiempo Base: Tiempo de Ascenso: Tiempo de Retraso: Duracin:
tB
tP
thr
3hr 2hr
tL
2hr
0.2
0.5 0.0
hr+l
--T-
I = 1lD = 112 plglhr
450
te= 3 h
ll3f.l
,oo
200
l5[
I0[
t [hrl
lg =9 hr
Hidrograma Unitario para una duracin de la lluvia neta de 2 horas
LL9
d)
Verificacin:
Cuando el intervalo de tiempo At es constante, podemos calcular el
volumen utilizando la regla de Simpson para verificar los resultados de

siguiente manera:
1a
Volumen
= R.A = (LQtt
Por consiguiente, la profundidad de escorrenta en el hidrograma unitario es:
It
O )ar
A
(1,7s0f /s) (1hr) (3600se/hr)

(1,7
3 5.5 a c re
$ (43,5
60 pz / a c r e)
R =
(0.0832
pie) x (t2plgll pie) =
plg
METODO DE LA CT'RVA
Un hidrograma unitario para una cuenca particular est definido para

una duracin especificada de la precipitacin neta. La propiedad lineal del hidrograma ruritario se puede utllizar para generar un hidrograma unitario de
una duracin diferente. Por ejemplo, si se conoce el hidrograma unitario de
una cuenca para una duracin de la precipitacin efectiva de una hora, el

hidrograma unitario que resultara de una tormenta de dos horas de duracin
se puede generar utilizando dos hidrogramas
unitarios de una hora cada uno,

desfasados una hora uno con respecto al otro, sumando las ordenadas
L20 y luego
dividiendo el resultado entre dos. De esta manera, una pulgada de escorrenta directa en una hora, se ha distribuido uniformemente en dos horas de duracin
se ha derivado el hidrograma unitario para dos
horas. Este procedimiento
de desfasamiento esta limitado a mltiples enteros de la duracin original del
hidrograma unitario.
Figura N"20: Hidrograma Unitario por Desfasamiento

El mtodo de la Curva S permite la construccin de hidrogramas
unitarios de cualquier duracin. En este caso se supone que se conoce un hidrograma unitario de duracin D y que se desea generar un hidrograma
unitario para la misma cuenca con una duracin distinta
D'.
Elprimer paso
para obtener las ordenadas de
la
curva
,S
es el de sumar una serie de
Hidrotoga / David Cedeo

hidrogramas unitarios de duracin
L2L
D,
desfasados uno con respecto a otro un

a
tiempo D. El hidrograma de escorrenta directa resultante corresponder
una precipitacin de intensidad:

hidrograma de equi librio.
1/D: el cual es equivalente a un
Si desplazamos la curva
,S
en el tiempo una magnitud
D'
y restamos
las ordenadas entre las dos curvas S, el hidrograma resuitante se producir
por una precipitacin de intensidad
l/D,
que ocurre con una duracin
D';
por lo lanto, debemos ajustar las ordenadas multipiicaldo por la razn D/D'
para obtener el hidrograma unitario para una duracin D'.
FtrcFFAc{tM|A p,i1,4 Et MEro
t{B{ffi.'
Figura N'21 : Hidrograma Unitario por el Mtodo de la Curva S
Hidrologa / Davicl Cecleo
122
J
En muchas ocasiorres se ha observado que despus que la curva
alcanz.a un vlloL mximo, el valor de la suma oscila alrededor de este valor
mximo. Esto
se clebe a errores pequeos en la duracin del hidrograma
unitario original, debido al comportamiento no
lineal. Cuando
se presen[an las
las oscilaciones, el hiclrtogo debe efectuat' ajustes para eliminar

fluctuaciones en la derivacin cle curva hacia clelante y hacia atrs fue sugericlo por
s. Un mtodo modil'icado de curva,s 'lauxe (1978), colllo una fflanera
de emparejar y eliminar el probtema de las fluctuaciolnes.
tequilibrio = lg
{p
llnl
Hid rogram
lJnitrio P6ra
lidrograms Uniltrios
defasados 2 hor
s
Figura N"22: Curva
Hidrologa / David Cedeo Ejemplo 16:

Mtodo de la Curva
S
L23
Uttlizar los siguientes datos de las ordenadas del hidrograma unitario

de dos horas de duracin de la precipitacin neta para obtener un hidrograma
unitario de tres horas, aplicando el mtodo de la curva S.
Ordenadas del Hidrograma Unitario

Tiempo (horas)
10
D:2hr
(p3lseg)
H.U.
75
2so
300
275
200
100
75
50
25
Solucin: Generacin de la curva S: Los datos originales del hidrograma unitario para D
2 hr se utilizan
para generar la Cwva,f, desfasando varios hidrogramas dos horas y sumando las ordenadas. La Curva,S representa un nmero infinito de esas adiciones.
Sin embargo, generalmente es necesario repetir el proceso hasta que se

alcance un valor constante; el cual ocurre generalmente en el instante:
t = t"-D
donde t" es el tiempo base del hidrograma unitario original y de la precipitacin neta.
es la duracin
Hidrologa
/ David
Cedeo
1-24
Tabla de Clculos de la Curva

Tiempo (hr)
(D :2
0
75
H.U.
hr)
Ordenadas Desfasadas 2 horas
Suma
(p'ls)
0
75
0
1
2
3
250
300 275
0
75
250
375 0 75
4
5
2s0
300 275
525 575
200
100
2s0
300 275
0
75
625
7
8
't5 50
25 0
200
100
650
0
'75
2s0
300
275
615 675 0 75 675 675
9
10
11
75
200
100 75
50
25
250
300 275
200
100 75
12 13
50
25
250
300
n
75
675 675
200
Alternativa:
La curva
,S
tambin se puede obtener utilizando un procedimiento
condensado del mtodo anterior, que consiste en desfasar las ordenadas

acumuladas (suma) un intervalo de tiempo igual a la duracin D.
Hidrologa / David Cedeo Mtodo Alterno para el Clculo de la Curva

Tiempo (hr)
0
1
].25
S
H.U.
Ordenadas Acumuladas desfasadas 2 horas
Suma
(D:2hr)
0 75
(p'ls)
0
75
2
3
2so
300
275
0
75
zso
375
525 575
625
4 5 6 7
8
2so
375
200
100 75
525 575
625
6s0
675
50
25
6s0
675
675 67s
675 675
10
1i
12
61s
675
b.
Hidrograma Unitario para una duracin diferente: El hidrograma unitario para
D'
3 horas.
D'
3 hr, se puede obtener a partir de la
curva S, desplazando la curva S tres horas y restando las ordenadas de la curva original. Estos valores deben ser multiplicados por la raz6n entre la
duracin D del hidrograma unitario original y la duracin
D'
del hidrograma
que se desea derivar; esto nos dar el valor de las ordenadas del hidrograma
unitario derivado. En este caso, el ajuste de las ordenadas es:
L26
D/D' :213
Los clculos de las ordenadas del H.U. para
siguiente tabla.
D'
3 hr se presentan en la
Ordenadas del Hidrograma Unitario obtenidas a partir de la Curva S
Tiempo (hr) t,
1
Curva S Original
Curva
Desplazada
Diferencia 0
75
D':3hr
(p'ls)
0.00 50.00
166.6'7
H.U.
(p'ls)
0
75
3hr
2
3
250 375 52s 57s

625 0 75
250
375
250.00 300.00 216.67

166.67 83.33 66.67 JJ,JJ 16.67
4
5
450 32s 2so

125
2s0
375
525
5'75
7
8
650
675 675
i00
50
25 0
62s 650
675 675
l0
11
67s
675 675
0.00 0.00
12
Hidrologa / David Cecleo
r2't
o tpshl
1000 800 600 4(l0 200
curva
original
curya S desplozada
3 hr
I fh4
a. Deslasmiento
de las curvas S
o[p /sl
300
200
100
t [h4
10
1l
b, Hidrograma Unilario para I horas de duracin
Derivacin de Hidrogramas Unitarios para cualquier duracin a partir de fa Curva S
L2B
METODOS MATRICIALES PARA DESARROLLAR HIDROGRAMAS
T]NITARIOS
El proceso para obtener un hidrograma unitario producido por
una
precipitacin excedente de perodos mltiples utilizando mtodos matriciales

se ilustra para el caso mostrado en la siguiente figura.
t I hrl a.
Hidrogrma de Precpitcn Excedcnte
123t5 b.
Hdtogram Unilario [pra una duracin de
t lhrl
h4
nrolq
1234561
c, Hdrograma
dr Escorrcnt Drccl
Figura N"23: Hidrograma Unitario por Mtodos Matriciafes .
L29
Observe que la precipitacin de la lluvia neta o efectiva es constante en

cada perodo; por lo tanto, la cantidad de escorrenta directa en cada intervalo se puede obtener de la siguiente manera:
P = (r-)Ar
En este caso, el sistema de ecuaciones resultantes para la escorrenta
directa producida por la precipitacin en los tres perodos es:
Q' Q,
Q,
Qu
Pr'Ut
Q,
Qu
Pr.U, + P1 .U2 P, .U, + P2.U2 + Pr.U, + Pr.U, + Pr.U, + P2.U4

P.'Uo
P1
.U3
Pl .U4
donde Q representa las ordenadas del hidrograma de escorrenta directa, U

las ordenadas del hidrograma unitario y P la precipitacin excedente durante
el perodo de tiempo. La relacin entre el hidrograma unitario original, la

precipitacin excedente y el hidrograma de escorrenta directa es la siguiente:
n = j+i-l
donde n es el nmero de ordenadas del hidrograma de descarga,
es el
nmero de ordenadas del hidrograma unitario y i es el nmero de perodos de
precipitacin excedente.
Hidrologa / David
Cedeo
13
El
sistema de ecuaciones se puede escribir de manera compacta

:
utilizando notacin matricial
tJ = tPl t71
donde:
a)a
vector ( 6 x 1 ) del hidrograma de descarga directa:

Q,
Q,
o^
tt
8o
o-
o-
b)P:
Matriz ( 6 x 4 ) de precipitacin neta:
Pt 0 0
0 o
PrPr o
tPl
P3P2Pto o P3P2Pl o 0 PrP, 0 0 0P,

c)
131
1 ) de las ordenadas del
vector ( 4 x
hidrograma unitario:
ur
u2
tul
u3
u4
Los mtodos clsicos para la solucin del sistema de
ecuaciones
requieren el clculo del inverso de la matriz de precipitacin: [P]-1, la cual existe cuando Ia maftiz es cuadrada con determinante distinto de cero. En
este caso, la matriz de precipitacin [P] no satisface sta condicin, pero
multiplicando por la matriz traspuesta [Pr], se puede generar una matriz simtrica cuadrada [PrP], que si tiene inversa y de esta manera podemos
obtener la solucin del sistema:
Vr Pt ' lul = fp, el

La cual se puede resolver explcitamente para obtener las ordenadas del
hidrograma unitario de la siguiente forma:
lul = lP' Pl-t v' gl

El procedimiento descrito anteriormente fue desarrollado por Newton
y Vineyard (1967),los cuales utilizaron el mtodo de Gauss para la solucin del sistema de ecuaciones simultneas. Sin embargo, hay que tener cuidado
1-32 la
con los procedimientos numricos, ya que debido al truncamiento,
propagacin de errores puede producir resultados pocos satisfactorios,

especialmente en las ordenadas finales del hidrograma unitario. El siguiente
ejemplo demuestra el procedimiento para obtener el hidrograma unitario

utilizando matrices.
Ejemolo 17: Obtener
Mtodos Matriciales
un
hidrograma unitario para una duracin de una hora
producido por una precipitacin excedente de tres horas de duracin, con tres
perodos de intensidad constante, la cual gener el hidrograma de escorrenta
directa cuyos valores se muestran en la tabla de datos y en la grfica.
Tabla: Precipitacin Neta

Intervalo de Tiempo (horas)
Intensidad Efectiva
0-1
2
r-2
3
2-3
1
(ple/hr)
Tabla: Descarga Directa

t
(horas)
(p'ls)
400
1,200
1,400
900
425
150
25
llidrologa / David Cedeo
133
i[ptgfhr]
a. llieloqrama de Precipitacin Neta
b. Hidrograrna de Descarga Direcla
Diagramas de Precpitacin
y Escorrentia
Solucirr:
a)
Cantidad de Datos:
Valores de Caudal Directo:

Perodos de Precipitacin: Ordenadas Unitarias:
j j j
)+1
+1
1-34
b)
Vector de descarga
Q (n x 1):
Q,
400
1,200 1,400 900 425
150 25
Q,
o^
t8l
Qo
o_
o -6,
o_
c)
Matriz de precipitacin P
(nxj
00 00 00 Pto
):
Pt 0
0 o
P, P,.
tPl
P3 P2 Pl
o P3P2 0 0P3 P" P, 000 P, P, 000 oP1

U
UL
20000 32000 t3200 01320 00132 00013 00001
d)
Vector del Hidrograma Unitario
(j x I ):
u2
tul =
U^
u4
U. f

e)
135
Procedimiento de Algebra Lineal:
tPl . I71
lPr pl . l
Matriz transpuesta
fP' 0l
(j x n ):
231 o23 002 000 000 0000 1000 310h 2310 o237
vrl
c)
Multiplicacindematrices
(j x j):
149 2 0 0 9149 2 0 2 9149 2 0 2 914 9 0 0 2 914
(
IP, PI
h)
Multiplicacin de una matriz por un vector

5,800
j x I ):
7,500
< 7
ot<
))<
1,325
Hidrologa
David Ccdco
136
i)
Sisten-ra de ecuaciones lineales simultneas:
l4 9 2 914 9 2 914 029 002

j)
Vector solucin
00 20 92 149 9t4
ul
u2 U. u4 U-
5,800
'7,500
5,925
3,225 1,325
(i x I )
U,
u2
200
300 l s0
75 25
tul
u)
un
U.
[h4
0
on y's"q
{;
o tpls es]
200
2 3 4
300
150 75
2S
0
5
6
Ordenados y Grafica rlel tlidr0gfti ta Unilario pEra
de
hr.
ull du.acn
L31
DESARROLLO DE HIDROGRAMAS TINITARIOS SINTETICOS

Los procedimientos presentados en la seccin anterior para derivar los
hidrogramas unitarios en cuencas con estaciones de medicin han sido

aplicados satisfactoriamente cuando existen datos adecuados de precipitacin
y escorrenta. Por otro lado, Ios mtodos para derivar hidrogramas unitarios
para cuencas sin estaciones de medicin estn basados en frmulas tericas
empricas que relacionan el caudal mximo y las caractersticas de tiempo
de Ia cuenca. El resultado de estos mtodos se ilenomina Hidrograma
Unitario Sinttico y los mismos ofrecen a los hidrlogos e ingenieros cierta

cantidad de procedimientos para desarrollar hidrogramas unitarios para
cualquier cuenca. Sin embargo, todas estas frmulas tienen suposiciones que
limitan su uso, por consiguiente deben ser aplicadas con precaucin. En

algunos casos, tambin se puede intentar una calibracin de la cuenca cuando
hay cuencas adyacentes con estaciones de medicin. Los Hidrogramas Unitarios Sintticos se han desarrollado a lo largo de
dos tendencias: una asume que cada cuenca tiene un hidrograma unitario
nico y la otra tendencia supone que todos los hidrogramas unitarios pueden ser representados por una familia nica de curvas una ecuacin nica.
La primera categora de desarrollo esta basada en el mtodo racional modificado para incluir las curvas de tiempo-rea para una cuenca en particular. Clark (1945) supuso que la respuesta de la cuenca sera producida
por el trnsito de las curvas de tiempo-rea a travs de un elemento de

almacenamiento lineal, el cual tiende a atenuar el caudal mximo y a retrasar
el hidrograma. Cada hidrograma unitario sera nico para la cuenca; por
Hidrologa
/ David
Cedeo
138
consiguiente, esta tcnica representa una mejora sobre
el
mtodo del
hidrograma de tiempo-rea. En general, se requieren expresiones empricas
para transformar los mtodos de Clark en tcnicas tiles para desarrollar

hidrogramas unitarios para cuencas reales.
La segunda metodologa para el desarrollo de hidrogramas unitarios

asume representaciones matemticas para la forma del hidrograma unitario.
El mtodo ms popular fue desarrollado por el Soil Conservation Service

(SCS) en 1964; el cual define un Hidrograma Urtario Adimensional de
descarga en funcin del
tiempo. Debido a que el volumen es fijo, solamente
se requiere un parmetro para determinar el hidrograma unitario en su

totalidad; el cual sera cualquiera de los siguientes valores: el tiempo al pico
t"
el caudal mximo Qo.
La mayora de los investigadores han sealado que existen varios

parmetros importantes para determinar
la forma y los tiempos del

de ascenso to y el tiempo
hidrograma unitario de una cuenca. Entre los parmetros utilizados ms a

menudo tenemos el tiempo de retraso
t, el tiempo
de concentracin
t".
El tiempo base t" tambin se incluye para definir la
duracin de la escorrentia directa. Estos parimetros de tiempo deben estar
relacionados con las caractersticas
de Ia cuenca para el
desarrollo
satisfactorio de hidrogramas unitarios sintticos. Otros parmetros de

importancia incluyen el caudal mximo Qp y el rea de la cuenca A.
La mayora de los mtodos para derivar
hidrogramas unitarios
sintticos asumen que el hidrograma unitario de una cuenca representan los
efectos combinados del tamao, pendiente, forma
caractersticas de
139
almacenamiento. Bajo estas condiciones, si los factores entre dos cuencas son
iguales y permanecen constantes con el tiempo, la respuesta ser idntica para ambas cuencas. Para dos cuencas del mismo tamao, si la pendiente de una
es mayor que la otra si la forma de la cuenca es ms concentrada (la
razn entre la longitud
cambiar
y el ancho es menor), la forma del hidrograma del tipo a (cuenca natural) al tpo (cuenca parcialmente
desarrollada); tal como se indica en la figura de los factores que modifican ei
hidrograma unitario. Las caractersticas de almacenamiento en la cuenca

esn
relacionadas con la pendiente, tipo de suelo, topografa, resistencia en
el canal y forma de la seccin transversal del canal. Por otro lado, si existe almacenamiento en un embalse en combinacin con Ia canalizacin aguas
abajo del embalse, el hidrograma tambin cambiar del tipo a al tipo b,

debido a un tiempo de ascenso menor por la reduccin en la resistencia en el
canal y menores prdidas por infiltracin debido al desarrollo de la cuenca.
Firnlmente, si la forma de la cuenca fuera alargada y estuviera completamente

urbanizada sin ninguna capacidad de almacenamiento (embalse); la respuesta
de la cuenca sera un hidrograma tipo c (cuenca completamente desarrollada); observndose un tiempo de ascenso mucho menor y un caudal mximo mayor que el hidrograma unitario del tipo
c.
Los primeros mtodos de hidrogramas
unitarios sintticos generalmente no consideraban los efectos del desarrollo de la cuenca; pero frmulas empricas ms modernas permiten incluir los cambios producidos por
la urbanizaciln de la cuenca, tales como
el
de
incremento en el porcentaje de reas impermeables

alcantarillas para el drenaje pluvial.
y la construccin
140
(a) cuenca
(b) Par:ialJtrn! desarrollada
() c pretanehte dsi!o.ltda
Figura N'24:
Factores que Modifican e1 Hidrograma Unitario de una Cuenca .

sintticos
La mayora de los mtodos para hidrogramas unitarios

relacionan el tiempo de retraso
t,
el tiempo de ascenso to del hidrograma
con las medidas de longitud y forma del canal principal de drenaje de la

cuenca. Algunos mtodos tambin relacionan estos parmetros de tiempo con el inverso de la pendiente del canal. Por lo tanto, para una longitud mayor
del canal y una pendiente menor, se debe observar un tiempo de ascenso

mayor en el hidrograma.
144
Tambin se puede presentar generalmente una segunda relacin entre

el caudal mximo Qp y el rea A de la cuenca entre Qp y el inverso de t'-
t, del hidrograma, lo cual nos indica que la magnitud
de Qo es proporcional
al rea tributaria; es decir que reas ms grandes producen un caudal mximo
mayor. Por otro lado debemos satisfacer continuidad, lo que implica que un
caudal mayor Qn esta asociado a un tiempo de retraso tL menor para mantener
el volumen de escorrentia directa del hidrograma unitario.
METODO DE SNYDER
Snyder (1938) fue el primero en desarrollar un hidrograma unitario
sinttico basado en estudios de cuencas localizadas en los Montes Apalaches
(U.S.A.) que tenan un rea de drenaje que variaba de 10 a 10,000 millas

cuadradas. Las relaciones empricas obtenidas por Snyder son:
tL = C,(L
L"\03
Donde:
tL
tiempo de retraso de la cuenca (hr) longitud del canal principal desde la salida hasta la
L
L"
divisoria de aguas (millas)
Iongifud a lo largo del canal principal desde la

salida hasta el punto ms cercano al centroide de la
cuenca (mi)
cr
coeficiente que variaba de 1.8 a 2.2 enlos Montes

Apalaches

La descarga mxima es:
L42
/ o =640co " ,,.
donde: Qo A Cp
caudal mximo del hidrograma unitario (p3lseg) rea de la cuenca (millas cuadradas) coeflciente de almacenamiento que variaba de 0.4 a
0.8. Observe que los valores mayores de Co estn

asociados con valores menores de Cr.
Y finalmente para el tiempo
base (en das):

f_
toB = 31---a
donde
t,
esta en
horas. La expresin anterior para el tiempo
base produce
resultados adecuados para cuencas grandes, pero valores excesivos para

cuencas pequeas. En este ltimo caso, es conveniente asumir de 3 a 5 veces
el tiempo de retraso tL para definir el tiempo base t" en el trazado

hidrograma unitario; es decir:
del
3<a<5
tL
r43
Q,L
s_
D ___!
Qp
ats
oso
Figura No25: Hidrograma Unitario Sinttico de Snyder
El procedimiento anterior define un hidrograma unitario sinttico para

una duracin de precipitacin excedente
D
t-
(horas) igual a:
D=L
5.5
Hidrologa / David
Cedeo
(hr):
a44
Para una duracin de precipitacin neta diferente
D' (hr), se debe ajustar la
frmula para el tiempo de retraso
t,_
ti = tr * o.2s (D/ - D)
donde
r/ (hr) representa el tiempo de retraso ajustado para la duracinD'. Al aplicar este mtodo emprico a cualquier cuenca, se debe tener
otra. Por consiguiente
se puede
cuidado con la seleccin de los coeficientes C, y Cp,ya que se ha observado

que varan considerablemente de urn regin a
efectuar una calibracin utilizando cuencas localizadas en los alrededores con estaciones de medicin para obtener el valor de estos coeficientes antes de
aplicar el mtodo a otras cuencas en la regin.

Despus que se han calculado los tres parmetros:
tr,
Qp
y t",
se puede
dibujar un hidrograma unitario de tal manera que el rea bajo la curva

represente 1 plg de escorrenta directa sobre la cuenca. Para construir el
hidrograma unitario por el mtodo de Snyder se pueden tttlizar puntos adicionales que corresponden al 50Vo y 75% del caudal mximo
Q. El ancho W (en
horas) para estos caudales, se puede obtener de manera emprica por medio
de las siguientes frmulas:
470
(+)"
1-45
w_. 5u -
830
lQ,
\^
\r.l
I
donde el caudal
esta en p3/s
y el rea A en millas cuadradas. Los tiempos
correspondientes al ancho W se deben distribuir de la siguiente manera: 1/3

antes de QrV 213 despus de Qo (en proporcin de
I a2).
Ejemplo
18:
Mtodo de Snyder
Utilice el mtodo de Snyder con la finalidad de obtener el hidrograma unitario sinttico, para una duracin de 1.0 hora de la precipitacin neta, en
una cuenca con las siguientes caractersticas:
100 millas cuadradas 1.8
cr ce
L L"
Solucin:
0.6
18 millas 10 millas
a)
Tiempo de Retraso
tL = C, (L L")03 =
1.8 (18
" 10)03 =
8.55 hr
a46
b)
Duracin de la Precipitacin Excedente

f_
D =
5.5
1.55
hr
c)
Aiustes para
D'
hora
ti = ttno.2s(Dt-D) = d)
Caudal Mximo
8.55
+02s(t -1.5s) =
B.4hr
Qo
640 C^
tL
640
x 0.6 x
100
4'57o Piel lseg
e)
Asumir Tiempo base (cuenca pequea):
tB " 44. = 4x8.4 =
33.6hr
Anchos
0,". =
t_470_470
0.75
Q- =
o.75
x 4,570 =
3,427.5
pief /seg
[2.33 hr : 4.67 hr]
7s
(QotAr,
(4,s70lroo)r
7
1
hr
Hidrologa
David Cedeo
r4'7
= 050Qe = 0.50x 4,570 =

830
2,285
piel /sett
wro
830 (Q,/Ar:
(4,s70l100)t
12.4
hr
14.13
hr:8.27 hrl
D',=1h
0p = 45IB
/l
075=3{27
L
I
0 2285 -50 =
f
I
4.77 6.57
8.9 13.57 1t 17
a = 30'9
hr
Hidrograma
de t hr por el- Mtodo de Snyder.
Unitario para una duracin
1-48
C)
Verificacin:
Rt
volumen
A
rea
prg
t2 ( 228s )'(4 ( 3427 .s
.11)
228s ).( 6.s7 - 4.77
s ,449
.72
r 12
,t4t .2s
) u2 (3427.s + 4570 ).( 13.57 - 8.9 )

U2 ( 4570
3427 .5
).( 8.9 - 6.57
9,317.08
: 1t2(228sX33.6-r7.t7) :
u2(228s + 3427.s).(r7.r7 -13.s7)
- t8,674.t6
10,282.s0
t8,77r.27
67
Volumen
,636 (pies3lseg) . hr
R=
too mi2
( s,zeo prt\t
1.05
plg > t plg
\ t nitt"
h)
Clculo del Volumen Unitario:
vot
= RA = (ro,r. tee)
\
lzptg
)
23\32O,OOO
roo mi2.
I
\
s.zso
p",
tmiua )]
)'l
I/ol =
Or
.(
lhr
3,600 seg
64,s33.3i eties3 /seg) hr

Ajuste del Tiempo Base:
L49
Vol =
(67,636
18,771.27) + % (2,28s) {tB

48,864.73
17.17)
64,533.33
+ 1,142.50 (tB - 17.17)
Por lo tanto, tenemos que:
tB
I 5,668.61
l7
1,142.50
.17 =
13
.71 +
77
.17
tB =
30.88
horas " 30.9 hr
METODO SCS (Hidrograma Unitario Triangular)

Este mtodo fu desarrollado por el Soil Conservation Service (1957),
una divisin del U.S. Departament of Agriculture,
est basado en un
hidrograma adimensional, el cual fu confeccionado utilizando una gran

cantidad de hidrogramas de diferentes cuencas que variaban en tamao y Iocalizacin geogrfica. El hidrograma est representado por un tringulo simple, con una duracin de la precipitacin neta D (hr), tiempo de ascenso
/" (hr), tiempo de descenso
ro
ftr) y caudal
mximo Q" (pies3/seg).
150
Figura N"26:
Hidrograma Unitarro
triangular del
SCS
Para el hidrograma de forma triangular mostrado en la figura anterior,
el volumen de escorrenta directa es:
vor
Por consiguiente:
= Q't' * 22
Q'to
2.Vol
to * t*

Utilizando una gran cantidad de hidrogramas, se encontr que:
151
: 'R3'P
_t5
Por lo tanto, el caudal mximo resulta ser:

t,
075
Vol
tP
Recordando que Vol
R. A, para una precipitacin excedente neta de I plg
(Hidrograma Unitario) tenemos :

484 A
P
tP
Donde
Qp A tp
caudal mximo (pies3/seg) rea de la cuenca (mi2)
tiempo de ascenso (horas)
En la figura del hidrograma triangular se puede observar que:
t- =
'2
-+tr
Donde: D tL
duracin de la precipitacin neta (hr)

tiempo de retraso (hr) medido desde el centroide de
la lluvia hasta Q"
L52
El tiempo de retraso t, se puede estimar utilizando cualquiera de las

frmulas empricas desarrolladas por el SCS, tales como:
tL
ro8 ,s
reoo
1;07
/t
Donde: tL L Y S
tiempo de retraso (hr) distancia hasta la divisoria de aguas (pies)

pendiente promedio de la cuenca en porcentaje (%)
retencin potencial mxima (plg)
.9
1'ooo
CN
-to
Donde Cly'es el nmero de curva que depende del tipo de suelo y uso de la
tierru. La siguiente tabla muestra los nmeros de curva (CA) para diferentes
usos de la tierra, donde los grupos hidrolgicos de suelos son:
A B C D
humedad en el suelo.
suelo arenoso, bien drenado suelo arcilloso-arenoso suelo arcilloso suelo arcilloso plstico (se expande al mojarse)
Los valores de Cly' mostrados en la tabla suponen condiciones normales de

Tabla CN
Nmeros de Curva para escorrenla en terrenos utiJizados para agricultura, zonas rurales y urbanas (condicin de humedad artecedente tipo tr) Grupo Hidrolgico de Suelo
153
Descripcin del Uso de la
Tiea
B
Terrenos cuhivados : Sin mtodos de conservacin de suelos Con mtodos de conservacin de suelos
Pastos campos de pastoreo: Bn mal estado En buen estado
72 62
81 88 7L 78
9r
81
68 79 86 39 - 61, 74
30 58 71
89 80 78
Praderas en buen estdo

Bosques forestas: Vegetacin poco densa con cubiefta pobre, sin maleza
Buena condicin (vegetacin protegida del pastoreo) Espacios abiertos, prados, parques, campos de golf, cementerios, etc. Buena condicin: el csped cubre el 75% 6 t6s del fuea Condicin regular: el csped cubre del 50 al75% del fuea reas comerciales o de negocio (85 % impermeables)
45 25 39 49
89
81
66 55 61 69
92
88
77 70 74 79
94
91
83
'7'7
80 84 95
93
Distritos o zonas industriales (72% imJ'enneables)

Zons residenciales : Tamao promedio del 1/8 acre o
lote menos l/4 acre 113 tcre l/2 acre I acre
porcentaje impermeable
65 % 38 % 30 % 25 % 20 %
77 61 5'7 54 51
98
85 75 72 '70 68
98
90 83 81 80 79
98
92
8',7
86 85 84
Patios de estcionamiento pavimentados, areas techadas, calzadas o vlas de acceso, etc. Calles y caminos: Pavimentadas con cunetas y alcantarilla pluviat Grava
98
98
76
'72
Tierra
98 85 82
98 89 87
98
9l
89
Fuente: SCS Tecbnical Report, "Urball Hydroloy for Sflall Watersheds " (1986)
Hidrologa
/ David Cedeo
r54
El SCS tambin desarroll un hidrograma unitario adimensional para

transformar el hidrograma triangular en una curva continua utilizando los
parmetros tn y
Qp. Este hidrograma unitario adimensional
se muestra en la
siguiente figura, junto con la tbla de los valores'. t
I tp y Q / Qo.
Coordenadas del Hidrograma Unitaio Adimensional del SCS
tltp
0.000 0.100 0.200
0.300 0.400 0.500 0.600
Q/Q"
0.000
0.015 0.075
QP
0.700
0.800
0.160 0.280 0.430 0.600 0.770

0.890 0.970
1.000
0.900
1.000 1. 100
t.200
1.300 1.400 1.500 1.600
1.800
0.980 0.920 0.840
2.000 2.200 2.400 2.600

2.800 3.000
0.750 0.660 0.560 0.420 0.320 0.240 0.180

0.130 0.098 0.075 0.036 0.018 0.009
3.500 4.000 4.500

5.000
0.004
Fiq. N" 27: Hidrog'rama UniLaric Adi mens ional
Hidrologa / David Cedeo Ejemplo 19:

Mtodo SCS
155
Derivar un hidorgrama unitario triangular utilizando el mtodo SCS

para una duracin de la precipitacin neta de 1.5
es de 100 mi2,
horas. El rea de la cuenca
constitudapor 50% de suelo tipo arenoso-arcilloso y el resto
es suelo arcilloso; el
ffi%
de la tierra tiene uso residencial (lotes de 1/4 acres)
el
40% restante se utiliza para pastos (en buena condicin); la pendiente
promedio de la cuenca es 100 pies/milla y la longitud de la corriente principal

es de 18 millas.
Solucin:
a) Datos:
l.-5 hr
100 mi2 18 millas
L
Y
' 5,280 pies/milla
95,040 pies
lOOpies/milla'
1 mi11a/5,280pies
. 100
1.9%
b)
Nmero de Curva:
Uso de la tierra
Residencial (lotes de 1/4 acres)
Pastoreo
Tipo de suelo
B C
Porcentaje del rea
CN
75
83 61
0.6'0.s =
0.6' 0.5 = 0.4. 0.s :
0.3 0.3 0.2 0.2

1.0
B
C
0.4'0.s =
suma
74
Hidrologa / David
Cedeo
74.4
15 6
CN =
0.3
x75+0.3 x83+0.2x61 +0.2x74 =

( nota: redondear al valor entero
CN = 74
ms prximo )
c)
Retencin Potencial Mxima:
^s
l'ooo - to = 3.s ptg to = t{9,0 = CN 74
d)
tL
Tiempo de Retraso:
o'a 15
+ 1oz
t,soo
tfl
(95,040)08
(3.5
y' 1.9
1oz
10.5 hr
r,900
e)
Tiempo de Ascenso:
tp =
*.,,
l-5 *lo.s = tt.25hr
Caudal Mximo:
^ Q, '
484
484
tP
"
100
11.25
4,300
Piel lseg
g)
Tiempo Base:
D-
Volumen
Area
%t"Q,
15'7
tB
2AR 8,
70o mi
5,280p
lmi
tB =
h) Alternativa:
30 hr
tR=
tB
55 zt,
=
tR =
(
3
1.25)
18.75
hr
11.25
+ 18.75 =
3O
hr
D=f.5hr
i=
UD
op= 43oo P3s
tL= 10.5 hr
tp =
11
.25 hr
te= 30 hr
Hidrograma Unitario lrianglr Fara 1.5 hr de precipitacin neta
158
APLICACIONES DE HIDROGRAMAS T]NITARIOS

Cuando se ha calculado el hidrograma unitario para una duracin dada
de la precipitacin excedente se han utilizado mtodos sintticos para

desarrollar el hidrograma de respuesta de una cuenca bajo ciertas condiciones
psiogrficas, estos hidrogramas se pueden utilizar para diferentes clculos hidrolgicos; entre ellos tenemos los siguientes:
1)
Obtener hidrogramas de descarga para precipitaciones excedentes de
mltiples perodos.
2)
Disear hidrogramas para tormentas con un perodo de recurrencia

seleccionado
(T
10, 25, 50 100 aos) utilizando procedimientos de
convolucin y desfasamiento.
3)
Determinar las variaciones en el hidrograma de descarga debido a los
cambios en el uso de la tierra, modificaciones en los canales de drenaje, almacenamiento en la cuenca y otros factores que afectan el
hidrograma unitario.
4)
Los hidrogramas para cuencas grandes que consisten en

subcuencas se pueden simular
varias
por procedimientos de
adicin,
desfasamiento y trnsito de avenidas para los caudales producidos por
hidrogramas unitarios en cada subcuenca. Los efectos en la variacin
en el patrn de intensidad de la precipitacin y uso de la tierra
se
pueden comprobar en la respuesa hidrolgica para la cuenca completa.
159
5)
Los mtodos de trnsito de avenidas se pueden utilizar para analizar la
atenuacin del caudal mximo
retraso de Ia onda que sufrir el
hidrograma de entrada al pasar por un embalse

obtener el hidrograma de salida.
y de esta manera
CONVOLUCION DE HIDROGRAMAS TNITARIOS El procedimiento para derivar un hidrograma de descarga a partir de

un perodo mltiple de precipitacin excedente repres&a un caso general del mtodo de tiempo-rea. Las ordenadas del hidrograma unitario se multiplican
por la precipitacin neta, se desfasan y suman en secuencia para producir las

ordenadas del hidrograma de descarga directa. Cada hidrograma unitario se incorpora en el instante correspondiente al inicio del perodo de
lluvia. El flujo
base se puede aadir para producir un
hidrograma ms realista de descarga total cuando los valores tpicos de flujo

base esln disponibles para la cuenca bajo
estudio. Se debe tener cuidado en
los clculos, especialmente se debe verificar que los incrementos de tiempo
de la precipitacin excedente correspondan a Ia duracin del hidrograma
unitario. La ecuacin en forma discreta que gobierna el hidrograma

tormentas mltiples se denomina ecuacin de convolucin:
para
o = Iru =l
n-t+l
160
Donde: Q, Pj Uj
ordenada del hidrograma de descarga
precipitacin excedente
ordenada del hidrograma unitario
(nota:j:n-i+1)
Observe que cuando la precipitacin tiene perodos mltiples con
intensidad constante, pueden ocurrir intervalos sin precipitacin.
Ejemplo
20:
Derivacin de hidrograma para pero?os mltiples de lluvia.
Derivar el hidrograma de descarga total para una cuenca que se

produce como respuesta a una precipitacin con perodos multiples de
intensidad constante, utilizando los siguientes procedimientos:
a) b)
Mtodo de Convolucin
MtodosMatriciales.
El hietograma de la precipitacin excedente (despus de considerar las
prdidas por infiltracin y evapotranspiracin) y el hidrograma unitario para
una duracin D
t hr se muestran en la siguiente figura; adems,
asuma
que el flujo base es de 25 piesr/seg (constante).
Valores de la Precipitacin Neta

Intervalo de tiempo (horas)
0-1
1.0
1-2
1.5
2-3
0.0
3-4
0.5
P""," (ple)

Ordenadas del Hidrograma Unitario
16L
(hr)
D:1hr
(p'ls)
H.U.
150
300
450
360
270
180
90
Hietograma de Precipitacin Neta
Hidrog'rama Untario
paraD=1hora

Solucin:
1-62
a)
Mtodo de Convolucin:
= )-
=1
Tabla de Clculos
t (hr)
0
1
Pi
uj
(p3lseg) 0
150
(plg)
1.0
1.5
P,U
P,U
P.Uj
PoU
a,
(p3lseg)
0
150
0 0
22s
0 0 0 0 0 150
2
3
0.0
0.5
300
300 450 360 270

180
525
450
360 270
180
4s0
675
0
75 150 225 180
900 1,110 960

810
4
5
540
40s 270
135
7
8
90
90
0
0 0 0 0
s40
270
90 45 0
135 90
45 0
9
10
11
1-63
b)
Mtodo Matriciales:
t1J = rPr.rl
b-1)
Matriz de Precipitacin:
tPl
1000000 312 l 0 0 0 0 0 03/2 I 0 0 0 0 1/2 0 3/2 I 0 0 -0 0t/203/2 I 0 0 0 01120312 I 0 0 0 0t1203/2 0 0 0 0 t/20312 0 0 0 0 0ll2 0 000000112
1
b-2)
Vector del Hidrograma Unitario:

150
300 450
tul
360 270
180 90
L64
b-3)
Vector de Escorrenta Directa:

150
525
900 1,110
Pt - tul
960 810
540_
270
90 45
Hidrograma de Descarga Total
-19
11
1-65
c)
Hidrograma de Escorrenta Total:
O = O.+O.
Tiempo
Caudal Directo
Q
t
(hr) 0
1
Flujo Base Qr
(pies3/seg) 25 25 25 25 25 25 25 25 25 25 25 25 25
Caudal Total Q,
(pies3/seg) 25 175
(pies3/seg) 0
150
2
3
525 900
1, 110
550
925 1,13s
985 835 565
4
5
960 810
6 7
8
540
2'70
295
115
9
10
11
90
45
70
25 25
12
Observacin:
El caudal total mximo producido por la lluvia

sobre la cuenca es Qma*
de perodos mltiples
1,135 pies3/seg
ANALISIS DE FRECI]ENCIAS
INTRODUCCION
Muchos procesos en hidrologa deben ser analizados y explicados de
manera probabilstica debido a la
taixaleza aleatoria (casual) de los mismos.

en
Por ejemplo, es muy difcil predecir la escorrentia y la precipitacin
sentido puramente determinstico (en el pasado, presente fufuro), debido a que es imposible conocer cuantitativamente todos los parmetros que afectan
estos procesos. Afortunadamente, existen mtodos estadsticos disponibles para orgarnzar, presentar y reducir los datos observados de manera de que se
facilite su interpretacin y evaluacin. En esta seccin estudiaremos algunos mtodos estadsticos que nos permitan cuantificar la incertidumbre de los
datos hidrolgicos en un marco de referencia estocstico.
VARIABLES ALEATORIAS
Una variable aleatoria es un parmetro (por ejemplo: precipitacin, escorrentia, etc.) que no se puede predecir con ceteza, debido a que las
variables aleatorias son el resultado de un proceso casual o incierto. Estas
variables deben ser tratadas estadsticamente como discretas
contnuas.
La mayora de los datos hidrolgicos son contnuos y deben ser analizados probabilsticamente utilizando distribuciones contnuas de frecuencias; por
ejemplo, los valores de descarga en un hidrograma pueden ser cualquier

nmero positivo real, el cual depende del aparato de medicin, es decir, que
pueden tomar cualquier valor dentro de cierto rango. Sin embargo estos datos

se presentan a menudo de manera discreta debido al proceso de
L67
medicin (por
ejemplo: valores mensuales, diarios u horarios de escorrenta).
Las variables aleatorias discretas solamente pueden tomar valores

enteros dentro de un conjunto de valores especificados; como ilustracin, podemos observar que al lanzar una moneda se producir como resultado:
cara sello, mientras que tirar un dado de seis caras producir cualquier
nmero entero del conjunto
1, 2, 3, 4, 5, 6
l.
Los datos observados en un
aparato de medicin deben ser tratados como valores discretos, ya que han
sido redondeados
cuantificados durante el proceso de coleccin de la
informacin. Sin embargo, esto es muy inconveniente para el anlisis de

frecuencias debido a la gran cantidad de intervalos de valores que se deben
considerar. Por ejemplo, si el rango de la magnitud del caudal de descarga

vara de 0 a 5,000 pie3/seg y los datos se han redondeado a valores enteros, existirn 5,0) intervalos discretos que se deben considerar. En este caso, es
mucho mas fcil suponer que los registros son contnuos. No obstante,
algunas veces se aplican distribuciones de frecuencia discretas a variables
contnuas (por ejemplo: profundidad de Ia precipitacin); pero la mayora de
las aplicaciones de distribuciones discretas en hidrologa ocurren en el caso de variables aleatorias que representan el nmero de eventos que satisfacen algn criterio; por ejemplo, el nmero de inundaciones que pasarn de cierta
magnitud durante cierto perodo de aos.
168
PRESENTACION DE LOS DATOS
Los datos observados a menudo se presentan en la forma de un

diagrama de barras histograma. La altura y forma general del histograma
es
til para caracterzar los datos y permite la seleccin de una distribucin
de frecuencia para ser aplicada a los datos; por ejemplo, la distribucin puede
ser sesgada (oblicua) o simtrica. Utilizando como ejemplo los datos de

escorrenta, el rango de caudales se divide en intervalos de clase y el nmero
de observaciones (frecuencia) correspondiente a cada intervalo de clase es
tabulado. El ancho de cada intervalo de clase debe ser suficientemente

pequeo para que se pueda observar el comportamiento de los datos
y a la vez
suficientemente grande para que el patrn no sea confuso. El valor utilizado
para los intervalos de clase puede alterar la impresin de los datos en el observador; por consiguiente, este valor deber poder ser alterado fcilmente
en los programas estadsticos de computadora de manera que los ingenieros
puedan comparar varias opciones. Panofsky y Brien (1968) sugirieron la

siguiente frmula como una ayuda para la seleccin:
k = 5log,o(z)
donde k es el nmero de intervalos de clase
es el nmero de datos
observados. Los intervalos de clase no tienen que ser del mismo ancho, ya
que algunas veces es conveniente agrupar los datos en intervalos mas grandes
de ancho variable.
1-69
Una manera de discretizar los datos contnuos de flujo sera la de

asignar una variable discreta aleatoria a cada intervalo de clase. Cualquier valor de caudal dentro de un intervalo de clase sera asignado al valor discreto
correspondiente
al intervalo,
generalmente
el punto medio, el cual
se
denomina marca de clase. En este caso la marca seal de clase se graficar en el eje de las abcisas del histograma.
Si las ordenadas del histograma se dividen por el nmero total
de
observaciones, el resultado se denomina frecuencirrelativa (probabilidad)

de cada intervalo de clase, de tal forma que la suma de las ordenadas sea igual
a 1, produciendo un mtodo alternativo para graficar los datos.
Para
comparar el ajuste de una distribucin terica a la distribucin emprica, hay
que normalizar los datos; es decir, dividir la frecuencia relativa entre la

amplitud del intervalo de clase y el resultado recibe el nombre de frecuencia
relativa normalizada simplemente densidad. En el histograma
de
densidad, el rea total de la grfica es igual a 1, lo cual corresponde al caso
terico. Otra manera de presentar los datos es la distribucin de la frecuencia
acumulada, la cual representa la suma de las frecuencias relativas del

histograma hasta un intervalo dado
y es la probabilidad de que un valor
cualquiera a 1o largo de las abcisas sea menor o igual a la magnitud de ese

dato. Tanto las frecuencias relativas como las frecuencias acumuladas se usan extensivamente en hidrologa, lo cual ser ilustrado en el siguiente ejemplo.

Ejemplo
L't
21:
Histogramas de Frecuencia.
La siguiente tabla muestra los caudales mximos anuales registrados en

una estacin de aforo en un ro durante un perodo de 31 aos. Construir los
histogramas de frecuencia relativa y acumulada.
Tabla: Serie Anual Mxima

Caudal
Caudal
Ao
Mrximo
fn3lseq)
9
Ao
19
Mximo
ln3 /seo)
1945 L946
L947
,840
6r
6,260 1,360 1,000 2,770

1, 3
5, l-70
L9 62
L,620
23s 15, 600
1963
L9 64
L948
t9 49 1950
19 51 1"952
]-965 :-966
L9 67
400
4,740
427
,2L0
1-l-0
1,
3,310
1968
L9 69
5,230
4,300
1953
4,400 7,760
2,
1954
1955
r-956
t97
L97
0
L
2,820
1,900
3, 980
520
340
t97
1957 1958 1959 1960
5,440
3,000
3
L97 3 L97 4 L97 5
6,560
4 , 71,0
,690
3,460
1-0, 300

Solucin:
L1L
a)
Ordenamiento o clasificacin de datos (de mayor a menor caudal):
Tabla: Datos Ordenados

Posicin
Caudal
Posicin
Caudal
15,600 10,300
17
18
3,310
2
1
3,210
3,000
9,840
19
4
5
7,760
6,560
20
21
2,820
2,770
6 1
8 9
6,260 5,440 s,230
22
23 24 25
')
\x\
1,900
t,620
1,400 1,300 1,110 1,000 427
5,t70
4,'140
10
11
26 27
28
4,710 4,400
4,300 3,980 3,960 3,460
12 13
29
30
31
t4
15 16
340 235
L'72
b)
Intervalos de Clase:
k = 5log,o (n) = 5log,o (31) =
7.45
Se utilizaran 8 intervalos de clase y vamos a considerar que el rango del
caudal es:
0<Q<
16,000 pie3/seg; por 10 tanto, la amplitud (constante) de
cada intervalos de clase es:
LQ c)
16'000-o = 8-
2,oooPies3/seg
Clculo de la Frecuencia, Frecuencia Relativa, Densidad y Frecuencia Acumulada:
Tabla: Distribucin de Frecuencias

Intervalo de
Clase (p3lseg)
Marca de
Clase (p3lseg)
Frecuencia
Frecuencia Relativa
Densidad
Frecuencia Relativa Acumulada
0 - 2,000 2,000 - 4,000

4,000 - 6,000 6,000 - 8,000 8,000 10,000
12,000 14,000
1,000
9 9 7
3
1
0.29 0.29 0.23 0.10 0.03 0.03 0.00 0.03
0.000145 0.000145
0.000115
0.29
3,000 5,000 7,000 9,000

11,000 13,000
0.s8
0.81 0.91
0.000050 0.000015 0.00001s 0.000000 0.000015
10,000 12,000 14,000 16,000
0.94 0.97 0.97

1.00
0
1
1s,000
Suma: n :31
L73
0.30 0.25
0 0.1 0.1
0.05 0.00
246
'10 12 14
CAUDAL
16
(10"p"/s)
HISTOGRAI|IA DE FRECUEHCIA RELATIVA

FRECUENCIA
RELATIVA ACUMULADA 1.0 0.8
0.6
0.4
o2
2
4 6 810 12 14
16
HISTOGRANA DE FRECUENCIA RELATIVA ACUIIIULADA
Diagramas de Barra
r'74
CONCEPTOS DE PROBABILIDAD
Considere un experimento aleatorio discreto con n resultados posibles:
{ *,,, *",...., x" }. Los resultados son mutuamente excluyentes si dos

cualquiera de ellos no pueden ocurrir simultneamente. Tambin se dice que
son exhaustivos colectivamente si tienen en cuenta todos los resultados

posibles del experimento. La probabidad de un evento X, se pueden definir
como el numero de ocurrencias del evento despus de un gran nmero de

ensayos. Esta probabilidad puede ser estimada de la siguiente manera:
(Xt)
n. I
n
donde
es el nmero de ocurrencias (frecuencia) del evento
Xi en un total de
n experimentos. Por consiguiente n,/ n es la frecuencia relativa probabilidad

de ocurrencia de X,. Una probabilidad discreta es simplemente la probabilidad de un evento
discreto. Si definimos P(X) como la probabilidad de un evento aleatorio X,,
las siguientes condiciones se satisfacen para las probabilidades discretas de

esos eventos cuando se consideran todos los N resultados posibles:
0 < P(x) <
= EPrx) t=l
1-l5
La probabilidad de la unin (simbolizada por u) de dos eventos

mutuamente exclusivos Xr y Xz es la suma de las probabilidades de cada uno,
es decir:
P(\ x) = P(x) + P(x2)

Los eventos X, y Yr son independientes si la ocurrencia de uno no
afecta la ocurrencia del otro. La probabilidad de
ocurrencia de ambos
eventos
interseccin (simbolizada
por n) cuando los
eventos son
independientes es el producto de sus probabilidades individuales:
P(xtn
Y) = P(x)'
P(r,)
Para eventos que no son ni independientes, ni tampoco mutuamente

exclusivos, se cumple la siguiente relacin:
P(4u Y) = P(x,) * P(rr) - P(\n
Y)
La probabidad condicional de un evento X, despus que el evento Y,

ha ocurrido es:
P,1Y,)
P(X,
=
yr
P(Y)
Hidrologa / David
Si los eventos
Cedeo
son independientes, entonces:
1-1 6
X, y Y,
Ejemplo
22:
Probabilidad Condicional
Y,,
Definimos X, como la condicin de que llover en un da dado y

representa la condicin de que habr relmpagos "n-rr.r
du cualquiera. Las
probabilidades de esos eventos son:
P(X,) P(Y,) P(Y1 lX1)
0.3 0.1
0.25
(probabilidad de lluvia
30%)
(probabilidad de truenos
l0%)
(si hay lluvia,la probabilidad condicional de truenos es de25%)
Calcular la probabilidad de que llueva y truene simultneamente ese da:
P(\n Y) = P(\
I Xr) ,.
P(Xr)
0.25
x 0.3
o.075
Adems, la probabilidad de que llueva truene durante ese da es:
P(xtu Y) = P(x,) + P(Y) - P(4a Y)

P(Xl
uY) = 0,3+0.1 -0.075 =
0.325
1-'7'7
VARIABLES ALEATORIAS Y
PROBABILIDAI)
DISTRIBUCIONES
DE
El comportamiento de una variable aleatoria discreta contnua se

puede describir por medio de su distribucin de probabilidad. Podemos asignar un valor numrico a cada resultado posible de un experimento de
acuerdo a una funcin de masa de probabidad (variables discretas) a una
funcin de densidad de probabilidad (variables contnuas).
En
hidrologa,
las variables
aleatorias
-discretas
se
utilizan
corrientemente para describir el nmero de ocurrencias que satisfacen cierto
criterio; por ejemplo, el nmero de caudales de descarga que pasarn de cierta magnitud (los cuales producirn inundaciones)
tormentas que ocurren por ao en una localidad.
Ocasionalmente, es conveniente tratar las variables contnuas de manera
el nmero de
discreta. Por ejemplo, en el caso de los caudales mximos (histogramas de

frecuencia) podemos utilizar las marcas de clase para describir la magnitud
del evento correspondiente a cada intervalo de clase. En este caso, si la
variable
representa
los
caudales discretos, podemos definir las
probabilidades de cada evento utilizando las frecuencias relativas:
(Xt)
n. I
De acuerdo a los resultados del ejemplo de histogramas (tabla
de
distribucin de frecuencias), tenemos las siguientes probabilidades en el caso discreto:

P(1,000) P(3,000) P(5,000) P(7,000) P(e,000) 0.29 0.29 0.23 0.10 0.03 0.03
L7B
P(l1,000)
P(13,000) P(15,000)
0.00
0.03
Observe que estos valores deben satisfacer los axiomas de probabilidad

mencionados anteriormente
:
< P(X) <

k
P(Xt)
donde k es el nmero de intervalos de clase. Adems, para probabilidades discretas tenemos las siguientes relaciones:
P(a<X<b) = t a<Xb
P(X)
1-7 9
Por lo tanto, podemos definir la funcin de distribucin acumulada como:
F(x) = P(X<x) = E Pfx,l X,,,

Usando los valores de probabilidad anteriores, tenemos que:
F(s,000) = P(X <s,000) = P(1,000)+P(3,000)+P(s,000)

F(5,000)
= 0.29+o.29+0.23 =
0.81
Las variables aleatorias contnuas se utilizan generalmente en hidrologa
para representar los eventos de caudal, lluvia, volumen, profundidad y
tiempo. Los valores utilizados no
estn restringidos a nmeros enteros, sin
embargo, muchas veces las variables contnuas se redondean a cantidades enteras. Para una variable aleatoria contnua, el rea bajo la funcin de densidad de la probabilidad
/(l)
representa la probabilidad; es decir:
f P(Xt<X<X) = Jf?lx xL
La funcin de densidad de probabilidad no es una probabilidad por s misma y tiene unidades correspondientes al inverso de las unidades de la
variable
X.
Sin embargo, a diferencia de otros clculos en ingeniera, estas
unidades generalmente
se ignoran.
Por el contrario, la funcin de densidad acumulada es de mucho

inters ya que es una probabilidad. La funcin de densidad acumulada para
r_80
variables contnuas se define de manera similar al caso de una variable

discreta:
f F(Xt) = P(-- <X<Xt\ = Jf?x

Algunas propiedades de esta funcin incluyen:
x,
o<F(x)<1.0
P(Xr<X<X) = F(X)-F(X)
Utilizando los datos de caudales mximos del ejemplo de histogramas
de frecuencia, podemos calcular algunas probabilidades considerando que el
caudal es una variable aleatoria contnua; por ejemplo:
P(Q<4,ooo) P(4,000<0<6,000) =
=
0.81
o.s8
-0.58 =
0.23
P(Q>6,000) = r-P(8 <6,000) = 1-0.81 = o.le
continuacin, podemos observar que al dividir las frecuencias
relativas por la amplitud del intervalo de clase, obtenemos la densidad o

frecuencia relativa normalizada y el rea total bajo el histograma de densidad
es 1 (una unidad), 1o cual satisface la siguiente definicin para variables

aleatorias contnuas:
1Bi
Area
_J r f(x) dx =
por consiguiente, la densidad es equivalente a una funcin emprica

densidad de probabilidad.
de
La seleccin del tipo de funcin de densidad de ia probabilidad para

representar los datos hidrolgicos es difcil, ya que existen varias opciones
que se ajustan a la fbrma del histograma de densidad. En las siguientes

secciones describiremos aigunas de las funciones tericas de densidad de
probabilidad ms utilizadas para analizar variables hidrolgicas.
Fig. 28: FUNCIOH DE 0EHSIDAD DE PROBAB|LIOAD. {para variables a lBato rias continuas}
L82
MOMENTO DE T]NA DISTRIBUCION

Una funcin de probabilidad de masa una funcin de densidad de
probabilidad son una forma funcional cuyo momento est relacionado con sus parirnetros; por
lo tanto, si
se pueden encontrar los momentos, entonces
tambin se pueden obtener los parmetros de la distribucin. A su vez, los momentos tambin nos indican la forma de la distribucin.
Para una distribucin discreta, el ensimo momento alrededor del

origen puede ser definido como:
Muo
= \x,rtqx i=l
y para una distribucin contnua como:
Moo
= I_: ,r f(x) dx
(,r)
El primer momento alrededor del origen es la media

para una distribucin discreta tenemos:
tv
valor
promedio, el cual es el valor esperado y se denota como E(X). Por lo tanto
E(n = |
x,r1x,
=1
Hidrologa / David Cedeo y para una distribucin contnua:
183
E(x) = I:xf@)dx
La media es una medida de tendencia central y tambin se le llama
parmetro de localizacin, debido a que indica dnde est situado en el eje
el centroide de la distribucin.
La esperanza
es un operador lineal, de tal rodo que si
ay
son
constantes, se satisfacen las siguientes propiedades:
E(a) =
E(bx) =
b E(x)
E(a+bx) = ct+bE(x)
Los momentos de orden mayor generalmente no se necesitan. Sin
embargo, los momentos centrales alrededor de la media se pueden definir
para una funcin discreta de probabilidad de masa como:
Mk =
E 6, - r,loptx,l
=l
Hidrologa / David Cedeo y para una funcin contnua de densidad de probabilidad como:
144
Mk = [: (,-r,)'f(x)dx
Estos momentos centrales representan el valor esperado de la diferencia
entre .r y el valor promedio pr, elevado a una potencia

el primer momento central es cero, debido a:
k.
Es evidente que
[-"U-t)f(x)dx
El
= [--xf(x)dx-u[:-f@)dx y
segundo momento central se denomina varianza
es muy
importante en estadstica, pues es una medida de la variabilidad. Para

variables aleatorias discretas tenemos:
I/AR(X)
= o2 = El(x - p)'l = E, tx, j=l

:
F)2 p(x,)
Para variables aleatorias contnuas
VAR(x)
= o2 =
E(x
- p)l'?l = I :" (x -
r,)2
f(x)
dx
Lavarianza es el valor esperado de la desviacin al cuadrado alrededor del

promedio y represena la escala dispersin de la distribucin. Una medida
185
equivalente es la desviacin estndar (o), la cual simplemente es la raz

cuadrada de la varianza. Utilizando las propiedades de la esperanza, tenemos
el siguiente resultado:
Var(x) = EIQ-tt)21 = E(x2-2tx+1t2

Var(x) = E(r')-2pE(x)+t2
= E(x2)-p2
Var(x) = E(x')-lEQ)
Observe que la varianza no es un operador lineal. Algunas relaciones
tiles incluyendo las constantes a y son:

Var
(a\ =
Var(b
x) = b2 Var(x)
Var(a + bx) = b2 Vdr(x)

En estadstica se pueden obtener momentos de mayor orden cuando se
necesiten, pero en hidrologa el parmetro ms utilizado corrientemente es la
asimetra (1), la cual se define como la raznentre el tercer momento central
y la desviacin estndar elevada al cubo:

M3
o3

La asimetra
186
y es el parmetro que define Ia forma de la distribucin.
Observe que
si la distribucin es simtrica alrededor del promedio, la
asimetra es cero.
Algunas veces es conveniente tener una medida normalizada de la

dispersin de la distribucin. El coeficiente de variacin (CV) se define
como la razrentre la desviacin estndar y el valor promedio:
CV
F
Una medida adicional de tendencia central es la mediana
(Xrr, la cual
no es un momento, pero representa el valor de X para el cual la funcin de distribucin acumulada es 1/2, es decir:
(X^) =
0.5
Otro parmetro de inters ocasional, que tampoco es un momento, es la moda de la distribucin. Este es el valor de X para el cual la funcin de
densidad de probabilidad (o la funcin de probabilidad de masa) alcanza su
valor mximo. La relacin entre el promedio, la mediana y la moda se

ilustra en la siguiente
figura. La mayora
de las distribuciones de inters de presentan
hidrologa son unimodales, pero en algunas ocasiones se

distribuciones mixtas que son bimodales.
Los momentos y parmetros analizados en esta seccin se refieren a las
distribuciones
de probabilidad fundamentales y se pueden derivar
1-8'l
analticamente. Existen frmulas generales para las frrnciones de probabilidad
de masa (discreta)
densidad de probabilidad (contnua) que se pueden
sustituir en las sumatorias integrales para evaluar los momentos en funcin
de los parmetros de la distribucin.
Estas relaciones proporcionan un
mtodo simple para obtener valores estimados de los momentos de la

distribucin.
MEDIA= MEDIANA = MODA
S I Ii4
ETR ICA
MODA MEDIANA
[NEDIA
SESGO POSTIVO
AS IMETRICAS
SESGO NEOATIVO
FlG. 29: EFECTOS DE LA ASlfitETRlA EH LA FUHCIOH DE DEHSDAD DE trROEAEILIT]AD.
1BB
ESTIMACION DE LOS MOMENTOS A PARTIR DE LOS DATOS
Si se conocen los valores numricos para los parmetros de una

distribucin, es posible generar una serie de variables aleatorias
Xl, X2, ...Xo
que pertenecen a la funcin de probabilidad de masa (discreta) la funcin
de densidad de probabilidad (contnua). Dicha serie de longitud infinita

constituira la poblacin de todas las variables aleatorias que pertenecen a las funciones mencionadas cuyos parmetros deflnen los momentos.
Los datos hidrolgicos observados a menudo -son el resultado de una

mezcla de procesos fsicos (por ejemplo, la escorrenta es el resultado de la
precipitacin) y por
1o
tanto incorporan una combinacin de distribuciones de

de
1o
probabilidad. En adicin, los datos observados estn sujetos a errores

medicin y no se ajustarn perfectamente a cualquier distribucin. Por
tanto, los valores de los momentos para la poblacin total calculados a partir
de los datos permanecern desconocidos. Sin embargo, se pueden obtener

estimaciones de dichos valores a partir de los datos observados para los tres
momentos de mayor importancia en hidrologa.
Si el
nmero
de
observaciones independientes de una variable aleatoria representa el tamao

de la muestra n, entonces una estimacin de la media valor promedio, ser:
i=x
Ex, r=l
Los momentos de mayor orden pueden estar sujetos a desviaciones

parciales en sus valores estimados. Una estimacin imparcial es aquella para
189
la cual el valor estimado de la muestra es igual al valor esperado de la

poblacin. Por ejemplo,
es una estimacin se puede demostrar que
16l
= tt ; por lo tanto
imparcialparc p.
i E(n = El llI yr,

t
=l
E i-l
nr
n
_
(Xt)
E (X)
I
11
g Lt= =t
62 =
52
n-1
6, - i)' E, i=l
'2
x,' - "f
n-l
-I
(en vez del valor intuitivo de n) elimina la
Donde el divisor n
desviacin parcial; es
decir: E(,s') = o2. Observe que para muestras ms
grandes un divisor gaal
a n producira
casi la misma respuesta; pero en la
prctica se prefiere el estimador imparcial.
Debido a que los momentos estimados son una funcin de variables aleatorias, estos resultan tambin variables aleatorias. La varianza de la
estimacin del promedio se puede calcular de la siguiente forma:
'2
Var
(X)
190
Por lo tanto, si interpretamos la varianza del promedio como una

medida del error en la estimacin del promedio; sta se reduce al incrementar el tamao de la muestra. Esta proposicin generalmente es vlida para todas las estimaciones de momentos.
La asimetra presenta problemas especiales debido a que involucra la
sumatoria de las desviaciones al cubo de la media y por consiguiente est

suj
eta a errores ms grandes en su evaluacin. Una estimacin imparciai
aproximada es el coeficiente de asimetra:
cs,=f
Desafortunadamente,
(n-r)(n-2)
la
(X,-E
'3
correccin apropiada para
la
parcialidad
depende de la distribucin respectiva. Una correccin
para
C^S, que se
utiliza en la distribucin Pearson Tipo
es:
CS, El valor corregido tambin
u).r, ,* n)
tttlizar con reserva en las otras
se puede
distribuciones cuando se requiera este parmetro. La estimacin de la

asimetra calculada con la ltima frmula se denomina estimacin en sitio, lo cual significa que la estimacin incorpora solamente datos de un solo registro de mediciones.
Las cuencas en desarrollo, con niveles variables de urbanizacin y

deforestacin, presentan un problema muy difcil en la estimacin de los
191
momentos debido a que la escorrenta tiende a aumentar con el tiempo. Por
lo tanto, las distribuciones de frecuencia pueden ser no estacionarias (varan

con el tiempo). En esta situacin, se puede utllizat la simulacin por medio
de modelos utilizando programas de

condicin dada de desarrollo.
computadoras para desarrollar
distribuciones de frecuencias para caudales en cuencas urbanas para una
Ejemplo 23:
Momentos de una Serie Anual Mxima
Utllizar la serie de 31 aos de caudales mximos registrados en la

estacin de aforo, usada anteriormente en el ejemplo de los histogramas de
frecuencia, para evaluar la media
y la desviacin
estandar de los datos
originales y transformados; adems, estimar la asimetra corregida.
Observacin:
En la estimacin de los parmetros se utilizaron
1os
datos originales:
x--o
y transformados de la siguiente manera:
Y =
LoCro
(Q)
Estos valores se muestran en la siguiente tabla, al igual que los

resultados de las diferentes sumatorias obtenidos por medio de una hoja
electrnica de clculo.
Hidrologa / David Cedeo Tabla: Datos Originales y Transformados (Tamao de la Muestra: n 31)
192
Ao
1945 1946 1947 1948 1949 1950
1951
Caudal:
X=
Y:LogQ
3.993 3.713 3.210 2.371
9,840 5,170 1,620

235
15,600
4.t93
3.616 2.630 3.520 3.643 3.890 3.401
4,740
427 3,310
t9s2
1953
t954
1955
1956 1957 1958 1959 1960

1961
4,400 7,160 2,520 340 s,440 3,000 3,690

10,300
2.53r
3.736 3.471 3.567
4.0t3
3.797
6,260
1,360 1,000 2,770 1,400
t962
1963
3.t34
3.000 3.442 3.146 3.507 3.045
t964
1965
t966
1967 1968
3,210
1, 110
t969 t970
s,230 4,300 2,820

1,900 3,980 6,560
3.1t9
3.633 3.450
3.2'79
t97t
t972 t973 t974 t975
Sumatorias:
3.600
3.8t7
3.673 3.539
4,710 3,460
Ex : 128,462 Ex' :861,184,754

E1x-l':
1.690188x10"
Ey : EY' : E(v-i)' :
107.3454 377.0935 - 1.9956

Solucin:
193
a)
Media:
x = Li n
i4
*, = +31
e2q462\
4,t43.e4p3/s
I = Lf n
b)
Yariarza:
--t
t, = + 3l
U07:4s4)
3.462s
s; = ;+ln,t-,,']
E'= +ls6t,ts4,7s4 -3rx
S"' =
q+t+t.o+2
1
10,961,550 .46 p6 /s2
s'2 v
,s,,t
I li",-,/f - n-r l7-'

-3rx
'
I
30
p77.og3s
12.+azt21
0.1794
10/1'
c)
DesviacinEstndar:
,s=
v
10,96
1,5 5
046
= =
3,310.82 p3 /s
s=
d)
Coefi ciente de Asimetra:
0.1794
0.4236
cs.
- x)' Ii.;', =l
.'3 J
(n-I) (n-2)
31
c.s-
7.6902
7012
?Ox29
(3,3 10.s2)3
1.659
cs.ty
t=l
E g, - f)'
,s3 v
(n-r)(n-2)
31
c.s.ty
2e)
1.9956
(30
"
(o.4?3q3
- 0.936
195
e)
Asimetra Corregida:
cs, =
('.*).',
/
(r.65e)
cs", = ( r . 6 )' rl
1.e8r
cs,., Y =(r.6'l
\
t /
"(-0.e36)= -r.n7
AJUSTE DE T]NA DISTRIBUCION A LOS DATOS

Un uso intuitivo para la estimacin de los momentos es el ajuste de la
distribucin de probabilidad a los datos por medio de la igualacin de los parmetros estimados obtenidos de los datos a la forma funcional de la
distribucin. Por ejemplo, la distribucin normal tiene dos parmetros:

Media r
Yarianza o2; el ajuste por el mtodo de momentos para la
distribucin normal simplemente estima los valores:
0=O
6=so
utilizando los resultados obtenidos para el valor promedio y la desviacin

estndar en el anlisis de
la
muestra.
1-96
Los mtodos grficos y el mtodo de probabilidad mxima son los dos

procedimientos altemos al mtodo de momentos para ajustar las distribuciones
a los
datos. A pesar de que estadsticamente se ha observado que el mtodo
de probabilidad mxima es superior al mtodo de momentos, generalmente

es mucho ms complicado de calcular que el mtodo de momentos
y por Io
tanto se encuentra ms all de esta introduccin a la hidrologa estadstica.
El objetivo principal de determinar los parmetros de una distribucin

es el de evaluar la funcin de densidad acumulada.- En algunos casos, se
puede lograr el mismo objetivo sin calcular los parmetros actuales de la
distribucin. En lugar de eso, la distribucin es evaluada utilizando el factor

de frecuencia
K,
el cual se define como:
X_X
s,
El valor de K
es una funcin del valor deseado para
la funcin
de
densidad acumulada y tambin dependen de la asimetra. Por lo tanto, si se
conoce el factor de frecuencia
K para el valor
deseado de
la funcin
de
densidad acumulada y la asimetra, se puede obtener el valor correspondiente
de X, de la siguiente manera:
X = X+KS,
Posteriormente durante
el
anlisis
de diferentes
distribuciones de
probabilidad, se ilustrar el uso de los factores de frecuencia K.
Hidrologa / David
Cedeo
L97
PERIODO DE RETORNO O INTERVALO DE RECT]'RRENCIA

La manera mas comn en hidrologa para indicar la probabilidad de un
evento es la asignacin de un perodo de retorno intrvalo de recurrencia para ese evento. El perodo de retorno se define de la siguiente manera:
Un evento anual m.ximo tiene un perodo de retorno ( intnalo de

recurrencia) de T aos s su magntud es igualada excedda, en promedio,
cada T aos. El recproco de T es la probabilidad de que el evento ser igualado excedido en un ao cualquera.
Por consiguiente, la inundacin con un perodo de retorno de 50 aos,
tiene una probabilidad de 0.02 2% d,e ser igualada excedida un ao

cualquiera. Es importante sealar que el perodo de retorno no implica nada acerca de las secuencias en el tiempo de los eventos hidrolgicos. Esta
inundacin de 50 aos no ocurre una vez cada 50 aos; por el contrario, se debe esperar que ocurran, en promedio, aproximadamente 20 inundaciones
con un perodo de retorno de 50 aos durante los prximos mil aos. Por
ejemplo, suponiendo que los eventos sean independientes, podran ocurrir dos
inundaciones con un perodo de retorno de 50 aos en dos aos consecutivos
sta secuencia de eventos tendra Ia siguiente probabilidad de ocurrencia:
P(Qrol Qr) = P(?r) " P(Qs) =
0.02
x 0.02 =
0.0004
El concepto de perodo de retorno implica eventos independientes y

generalmente se obtiene analizando los datos de una serie de aaudales
198
mximos anuales ( precipitacin mxima,
etc). El evento de mayor
magnitud en un ao dado, se asume que es independiente del mayor evento en otro ao cualquiera. Tambin es posible aplicar dicho anlisis a los n
eventos independientes mayores en un intervalo de z aos, sin considerar el
ao en que ocurrieron los mismos. En este caso, si el segundo evento en
magritud en un ao es ms grande que el evento de mayor magnitud en otro
ao, este se debe incluir en el anlisis de frecuencia. Esta serie de los n

valores mayores (independientes) se denomina Seriede Excedencia Anual,
en oposicin a la Serie
Anual Mxima. Ambos tipos de series
se
utilizan en
hidrologa, con muy pequeas diferencias en los eventos con perodos de

retornos grandes (eventos raros). Existen ms problemas asegurando la
independencia de los eventos cuando se utiliza la serie de excedencia anual,
pero para los eventos con perodos de retorno bajo, esta serie produce
resultados ms realistas para los perodos de retorno de eventos con la misma
magnitud que cuando se utiliza la serie de valores mximos anuales. Segn
Chow (1964), la relacin entre los perodos de retorno basados en Ia

excedencia anual
I" y los valores mximos f

T
es:
Ln(T^)-Ln(T.-l)
ttllizar todos los datos registrados
Otra opcin para el anlisis es la de
en la serie histrica sin considerar la independencia de los eventos. Esta serie
se denomina Serie Completa
y un ejemplo tpico es el registro diario de
caudales, la cual se llizar para determinar los caudales mximos anuales
L99
y los caudales de excedencia anuales. La informacin de frecuencia obtenida

de las series completas generalmente se muestran en una curva de duracin
de caudales (Searcy, 1959), la cual es una grfica de la magnitud del caudal
en funcin del porcentaje del tiempo que este valor del caudal es igualado o
excedido. La informacin obtenida del anlisis de la serie completa no
se
puede relacionar directamente con los perodos de retorno debido a que los
valores registrado en la serie completa no son independientes necesariamente.
CAUDAL
7o del
tiempo en que el caudal es superado
Fig.3O: Cun/a de Duracin

Los perodos de retorno pueden ser asignados a eventos mnimos (por
ejemplo: sequas) de una manera completamente anloga, pero
la
interpretacin "igualado excedido" significa "igual ms severo que esta
magnitud". Por consiguiente, un caudal bajo de 20 aos es aquel que tiene

una probabilidad e 5% de que la magnitud del caudal en un ao cualquiera
ser menor igual a este valor.
200
CLASIFICACION DE LOS DATOS

Para graficar los caudales cualquier otro dato hidrolgico, primero
debemos ordenar los valores desde
hasta n (numero de aos registrados) en
orden decreciente de magnitud (eventos extremos de inundaciones). Es decir que el dato de mayor magnitud tendr una posicin magnitnd tendr una posicin
: I
y el de menor
: n. El rango ru y el nmero de aos n se i.
ttlllzanpara calcular la posicin de los datos y para efectuar una estimacin

emprica de la frecuencia
Fy
el perodo de retorno
En hidrologa la clasificacin de posicin ms comn para obtener el perodo de retorno (
I ) es por medio de la frmula de Weibull:

T=
n+1
por consiguiente, la frecuencia ( F ) es: D-rfr

I T
n+7
Observeque:
1-
F(Q) = P(Q<Q) = I - P(Q>Q)i porloranro se satisface la siguiente relacin: 1/1 = P(Q , Q).El procedimiento
anterior fue analizado en detalle por Gumbel (1958); sin embargo, la

venerable frmula de Weibull ha sido ampliamente criticada debido a que no
proporciona una estimacin adecuada para la funcin de densidad acumulada

de tal manera que el valor esperado
E(')
sea igual al
valor terico para los

,?
datos observados de mayor magnitud, escogido entre el total de
muestras
20r
para cualquier distribucin; exceptuando la uniforme. Es decir, que las

distribuciones corrientemente usadas en los anlisis hidrolgicos para estudios
de frecuencia producen un valor desviado parcializado para la frecuencia
al utilizar la formula de Weibull. Gringorten (1963) propuso una forma generalizada para tratar de corregir estas deficiencias; sugiriendo las
siguientes frmulas para el perodo de retorno y la frecuencia:
n+1-2a
m-a m-a n+7-2a
F = 1-
donde el valor del parmetro a depende de la distribucin: es igual a0.375

para la distribucin normal log normal, 0.44paru la distribucin de Gumbel
se sugiere el valor de 0.40 para ste parmetro, como un compromiso
adecuado para las situaciones en que se desconoce la distribucin exacta de

la muestra.
A pesar de que las ltimas frmulas son mejores desde el punto
de vista terico, la formula de Weibull todava goza de gran aceptacin en hidrologa.
Eiemolo
24:
la
Perodo de Retorno y Frecuencia
Usar la frmula de Weibull con los datos de caudales mximos

registrados en
estacin de aforo (utilizados en diferentes ejemplos
desarrolados en esta seccin) para encontrar el perodo de retorno y la

frecuencia.

Solucin:
En la siguiente tabla se muestran los caudales ordenados de mayor
menor magnitud, mostrando su perodo de retorno utilizando la formula de
Weibull:
1-
n+l
m
y la frecuencia evaluada de la siguiente manera:
F = l-
Tabla: Clculos
Posicin
m
1
Caudal Q (p3ls) 15,600 10.300
Perodo de Retorno T (aos)
Frecuencia F
32.00
16.00 14.67
0.969 0.938 0.906 0.875 0.844 0.813

0.781
2
3
9,840
4 5 6
7,760 6,s60 6,260 5,440

5,230 5,170
8.00 6.40 s.33

A<1
7
8
4.00 3.56 3.20
0.750
9
10
0.119
0.688
4,740

Tabla: Continuacin
m
11
203
a
4,710 4,400
4,300 3,690
3,5 80
2.9t
2.67 2.46 2.29
0.6s6
o.625
12 13 14 15
o.s94
0.562
0.531
2.r3
2.00
1
16
3,460
3,310
0.500
l7
18
,88
0.469
0.438 0.406 0.315
3,2r0
3,000
1.78 1.68 1.60
19
20
2,820
2,770
2l
22 23
t.52
r.46
1.39
1.33
0.344
0.313
2,520
1,900
0.281 0.250 o.219

0. 188
24
25
t,620
1,400 1,360 1,110 1,000 427
t.28 t.23
1. 19
26 27
28
0.156
t.t4
1.10 1.07
1.03
0.12s
0.094 0.063 0.031
29 30
31
340 23s
244
1.mlD
,"""-'
0.cm
0.&m
0.7(m
I
I
ct 0.dm l!
r'
/ ./
I
-t
{ o 0.sm z 5 (,
IL 0.fm
I I
I
0.cm
it I
o2lp
0.1(m
J
/
ry
o^{rm
6{m gno rm
c
tDlL(cD
1ffi
4m
FlG.31: FUNCION DE DENSIDAD ACUMUIDA (OJIVA)
245
MODELOS PROBALISTICOS COMT]NES

En hidrologa se utilizan muchas funciones de probabilidad de masa
para variables discretas
funciones de densidad de probabilidad para
variables contnuas, pero solamente se van ha describir en esta seccin las
ms comunes. Es muy difcil reducir el nmero de distribuciones para

variables contnuas, ya que se pueden seleccionar diversas funciones para el
anlisis
de
caudales mximos;
sin embargo, las ms comunes
son:
Distribucin Normal, Log Normal, Gamma (de 2-parmetros), Pearson
Tipo 3, Log Pearson Tipo 3 y Gumbel Valor Extremo Tipo
I.
En
adicin a las funciones anteriores, se incluir Ia definicin de Distribucin
Binomial para variables aleatorias discretas y la Distribucin Exponencial
para variables contnuas debido a su simplicidad para ilustrar algunos

conceptos de probabilidad.
A menudo, el objetivo de un anlisis discreto es asignar probabilidades

para el nmero de ocurrencias de un evento, mientras que el objetivo de un
anlisis contnuo es el de determinar la probabilidad de la magnitud de un
evento. En el caso discreto, puede existir inters en la funcin
de
probabilidad de masa y en la funcin de distribucin acumulada; mientras que

en el caso contnuo, el valor de la funcin de densidad de probabilidad solo
nos permite comparar el ajuste de la distribucin terica a la distribucin emprica, mientras que la funcin de densidad acumulada para las variables
aleatorias contnuas nos define la probabilidad de ocurrencia de los eventos.
.n
lil
ql
t\
q)
lll \(dl
rrl
6l
,- ll *
ll
l:r
Ill
x*
la
r,lA '*la*
llrt 1n b
. lc
\o
uo
l
8
VI
v)
d
o a
AA
ol
cl
b0
o
A
R
VI
o
t :{
8
I
\d
c)
d
to
EI (!l
dt
N
-ol (t -l t til o.t
r
I
A
I
a.l
q)
h
I
()t
t cl EI
A X o
X !)
't
k'
c-
l
0)
ol
(Jl
a)
a OJ
c.)
cl .Ol
ol
,)l
c-l l
.<
ll
-lr loll
le tct
ll
...
ce-
t-r
tr.,,l
ut
El
c)
ll
h
o
d
\ot
ol
rl .rl
o
C
Cd
li
< d)
o rd
c)
q-.r
-ot
ot
X E]
z
J
9!
.
-ts t.s
co. d
|-\
e..
()
t-l
sl
EI .Fl il dl
Orl
Bl'.
t^ lx .t
I
(-
Io -.>
cat'
I
lx
l
lr
t
tsl*"
.<
C.r
l lx t.:
l\r
I .o-
cl
,Fr I \O F
arl*
I
a\
I
Fr/r
I
X
lt
qJ
C-
13
8
tu
vl R
VI
.n
.5FI '5lx
c)
Ol
bJN
A
uJ
)s
^,
8
I
*
d
d
-
'=l
:l
Ll il :l ^Fl lrt rl ) -l '=l
^l
q x o
I I
El 5 Ct ol
-a r
E - =
!r Hl ^\l
ca 0) X
!l ,-rl dl
q.,
-, l,-. lo l'-, Gl'-. rl

el
q)
cah;
u, ;.r
a
I
h
J
I
xl
>'
o.l 4. !l '<l
'gt
I I I
A X o
a
-t .Ol
o * -/
cn
tJ
<
o ti
()
l \C
ca)
F'; l H F;: 'El al

.ttl xl
l-{l
.xt
;jl
.!+
cd
ti l
.r,
P.
F
l-r
C)
i
a
()
F
()
E]
lJ.r
^oo
\-/
(!
bo
o?u
204
DISTRIBUCION BINOMIAL
Es muy comn examinar ura secuencia de eventos independientes para
la cual el resultado de cada evento puede ser xito o fracaso; por ejemplo, la inundacin con un perodo de retorno
de
aos ocurre o
no. Esta secuencia
consiste
en ensayos independientes de Bernoulli, para los cuales la
probabilidad de xito de cada ensayo es una constante igual a la fraccin p. Utilizaremos la distribucin Binomial para averiaguar cul es Ia probabilidad
de que ocurran exactamente
xitos
en n ensayoi
de Bernoulli. Observe
que sta ser la nica distribucin de tipo discreto analizada en esta seccin
y la misma
se utilizada en algunas ocasiones en Hidrologa, especialmente
para definir el riesgo y la confiabilidad de un proyecto.
La probabilidad de que ocurran

fracasos:
xitos seguidos
por ( n - x )
fracasos es el producto de las probabilidades de ocurrencia de xitos
S:
(1-p)
en
los
eventos independientes; es decir:
p' (l - p)'-'
Pero los ensayos de Bernoulli representan una secuencia posible para la
ocurrenciade
.r
xitos
y (n-x )
fracasos; por lo tanto debemos considerar
todas las secuencias posibles, incluyendo aquellas en las cuales los xitos no
ocurren consecutivrmente. El nmero de combinaciones maneras posibles
de escoger
ir
eventos entre
los
posibles resultados sta dado por el
coeficiente binomial:
nl
xl (n - x)l
209
La probabilidad deseada es el producto de la probabilidad de cada una
de las secuencias y el nmero de maneras en que pueden ocurrir

Distribucin Binomial:
las
secuencias; por 10 tanto, el resultado (funcin de masa de probabilidad) es la
(;)
p* (l - p)n' ; V = 0,1,2,3,...,n
La notacin B (n,p) indica una distribucin Binomial-con parmetros n y p. Para ste tipo de distribucin, la media valor esperado es:
t = E(x) =
la varianza es:
nP
02 = Var(x) = np(1 -P)

y la asimetra
es:
| - 2p Lnp (r - p)l%
En la expresin anterior se puede observar claramente que la distribucin es simtrica si la probabilidad de xito es
0.5. En este caso en particular,
la funcin de distribucin acumulada es:
F(x) = P(x<x) =
l(",)o,r-ob-t)
2L0
La evaluacin de la f,rncin anterior es tediosa para valores grandes de
y valores intermedios de x, pero afortunadamente existen tablas en los
libros de estadstica con los valores calculados para F(x).
RIESGOS Y CONFIABILIDAD Enuna secuencia de n aos, cul es la probabilidad de que el evento

con un perodo de retorno Z ocurra solamente una vez? La probabilidad
deocurrenciadel evento enun ao cualquiera es h f?accin

nmero de ocurrencias en la secuencia de n aos es B
de que el evento ocurra al menos una vez
: I I T y el
(n,p). La probabilidad
en n
ensayos se denomina Riesgo.
Este valor se puede calcular rpidamente considerando que riesgo es igual a

1 menos la probabilidad de que el evento no ocurra
en
aos; es decir:
Riesgo
= 1-(1 -p)'=
(' +)' (' +)'

diseo
Por otro lado, la Confiabilidad se define como 1 menos el riesgo; por lo

tanto tenemos la siguiente expresin:
Confiabilidad=
(l-pf
Los conceptos de riesgo y confiabilidad son importantes en el

hidrulicas.
hidrolgico, especialmente durante las etapas de construccin de obras

Ejemplo
21-].
25:
Diseo de Ataguas
Se va a construir una atagua para proteger una urbanizacin en una
planicie de inundacin mientras se completa un proyecto de canalizacin. La

afagtra se diseara para una inundacin que tiene un perodo de retorno de 20
y el proyecto de canalizacin requiere 3 aos para ser terminado. Describir el proceso estocstico y calcular las siguientes probabilidades
aos
durante la construccin del proyecto
a) b) c) d) e)
El nivel del agua sobrepasar la atagua en un ao cualquiera.

La atagua contendr la inundacin de diseo (Confiabilidad).
El nivel del agua sobrepasar la atagua exactamente una vez.

La atagua ser sobrepasada por lo menos vnavez (Riesgo).
La atagua ser sobrepasada solamente durante el ltimo ao.
Solucin:
En este anlisis se va ha investigar el xito fracaso del atagua para

contener los niveles de inundacin; por 1o tanto se trata de un caso discreto
y el proceso
se puede
definir por medio de una distribucin binomial:
B(n,p) = B(3,00s)
donde donde
es la duracin de la construccin en aos
y la proporcin
es el inverso del perodo de retorno Z.

Clculo de las Probabilidades:
21-2
a)
La
atagua falla en un ao cualquiera:
P(o>o..,\ Dry drteno. = . b)

Confiabilidad de la atagua:
I T
I
20
o.o5
Confiabilidad=
Confiabitidad
q' = (l -P)'
0.86
= (l -0.0s)3 = (0.95)3 =
c)
La atagua falla exactamente una vez en tres aos:
F(x=r) =
f''l \ *i
p'(r-p)(n-x' = f ''l p0-p)2 \ I/

(0.05)
F(x=r) =
+;
(0.e5),
= 0.135
d)
Riesgo de la medida de proteccin:
Riesgo= l-Confiabilidad= l-0.86 =
0.14
e)
La
atagua falla durante el ltimo ao (regla de multiplicacin):
Prob
= qxqxp = (r-p)"p =
(0.9s)2(0.0s)
0.045
21-3
DISTRIBUCION EXPONENCIAL
Considere una secuencia de sucesos aleatorios, de tal manera que los
eventos son independientes,
el
proceso es estacionario (es decir, los
parmetros del proceso no cambian con el tiempo) y no es posible tener ms
de una ocurrencia del evento en cualquier instante. Estas condiciones

describen un proceso de Poisson, el cual es representativo los eventos de tormentas. Si la variable aleatoria
representa el intervalo de tiempo entre
los eventos, se pueden decir que la variable
tiene una distribucin
exponencial con la siguiente funcin de densidad de probabilidad:
f(t) = ),eL';
La media de la distribucin es:
l>
t.r = E(t) _1
y la varianza
es:
a2 =
Var(t)
Esta distribucin tiene la propiedad caracterstica de que su promedio

es igual a su desviacin estndar;
por
10
tanto, el coeficiente de variacin
es
: l. La distribucin es obviamente sesgada hacia la derecha, con una asimetra constante y : 2.

CV
,)
Hidrologa / David
Cedeo
2L4
La funcin de densidad acumulada se puede evaluar analticamente,

integrando la funcin de densidad de probabilidad; es decir:
F(f) =
Jo'
exp
(- )" t) dt
= I - exp(-1. /)
t : * ,la funcin de
las
En la expresin anterior, es evidente que cuando
distribucin acumulada es:
F(-) = 1;
1o cual satisface una de
propiedades de este tipo de funcin, ya que el rea
total bajo la funcin
de
densidad de probabilidad debe ser igual

probabilidad de todo el espacio muestral.
a la unidad y
representa la
La distribucin exponencial puede ser manipulada fcilmente en forma
analtica
en consecuencia se utiliza algunas veces para aproximar
distribuciones asimtricas ms complejas. Tambin se emplea ocasionalmente

para encontrar la correspondencia entre la precipitacin total en la cuenca y
el caudal en el ro; pero se utiliza ms a menudo para obtener los intervalos

de tiempo entre los eventos. Algunas funciones de distribucin estrechamente
relacionadas con Ia distribucin exponencial son la distribucin de Poisson,
la cual define el nmero de ocurrencias en el tiempo

Gamma, la cual define el intrvalo de tiempo entre
y la distribucin
eventos.
Ejemplo 26:
Intrvalo de tiempo entre tormentas
Durante el transcurso de un ao, ocurren aproximadamente 110 eventos

independientes de tormentas en cierta localidad, con una duracin promedio
de 5.3 horas cada
una. Describir el proceso
ignorando las variaciones
2]-5
durante las estaciones del ao y calcular algunas probabilidades, las cuales se
van a definir posteriormente.
Solucin:
El proceso se puede describir utilizando una distribucin exponencial,
la cual tiene las siguientes funciones de densidad de probabilidad f

densidad acumulada
(t)
F(t)
f(t) = )'e-L';
F(t) = I -
t>o
exp(-.1.r)
El intrvalo de tiempo promedio (en horas) entre eventos de tormentas

en un ao es (1 ao
8,760 horas):
8,760 - 110 " 5.3 _ t=ff=743hr
El parmetro
de la distribucin exponencial est relacionado con el
valor promedio anterior, ya que:
E(t) -
I
)L
Por lo tanto, utilizando el mtodo de momentos, tenemos:
.r.
= 1 r
1
74.3
o.or35hrl

Probabilidades:
216
a)
Cual es la probabilidad de que existan al menos 4 das (96 horas) entre
tormentas?
P(t>96hr) = r -F(96) = exp(-0.0135x9) =
0.27
b)
Cul es la probabilidad de que el intervalo entre dos tormentas sea exactamente igual a 12 horas
?
P(t=t2hr) = F(12)-F(12) =
Observe que la probabilidad de que una variable aleatoria contnua sea exactamente igual a un valor especfico siempre es cero.
c)
Cul es la probabilidad de que la separacin entre dos tormentas sea
menor igual a 12 hr?
P(t<12hr) = F(12\ = 1-exp(-0.0135x12) P(t<l2hr) = l-0.85 =

o.l5
211
DISTRIBUCION NORMAL
La Distribucin Normal
tambin se conoce con
el
nombre de
Distribucin de Gauss y es de fundamental importancia en probabilidad y
estadstica. Una razn de esto, es que el Teorema del Lmite Central

establece que bajo condiciones muy generales, la distribucin de la suma de
una gran cantidad de variables aleatorias se puede aproximar a
una
distribucin normal, independientemente de de la distribucin verdadera de
la muestra. Muchos procesos fsicos se pueden ide-alizar como la suma de

procesos individuales; por lo tanto, la distribucin normal encuentra muchas aplicaciones en hidrologa y otras areas de estadstica, tales como prueba de
hiptesis, intrvalos de conflanza y control de calidad.
La funcin de densidad de probabilidad para la distribucin normal

(curva en forma de campana) es la siguiente:
f (x'l
, lt;
.,.0
|
(x - tt)'
2a2
-oo < f <
oo
Los prametros de la distribucin son el promedio
y la varianza o2; las
variables aleatorias que tienen esta distribucin se denotan como
N(r,
o2).
Debido a la naturaleza simtrica de la distribucin normal, esta no es

adecuada para analizar eventos extremos, tales como caudales mximos; pero es satisfactoria para describir totales anuales.
Hidrologa / David
Cedeo
2L8
La evaluacin de la funcin de la densidad acumulada requiere un

cambio de variable:
z = ( o
t)
Donde z es la variable normal estandarizada y tiene una distribucin N (0, 1);

entonces la funcin de densidad acumulada para la variable transformada es:
F(z) = |
,lt; '\ -l 2 )
/ "\ I -u'l exr.tl
du
Desaforhrnadamente, esta integracin no se puede efecfuar fcilmente,
pero existen tablas de F(z) en funcin de la variable z enla mayora de los libros de estadstica. Obsrvese que la variable estandarizada e corresponde
al factor de frecuencia ft utilizado en otras distribuciones. La magnitud de -r para un perodo de retorno dado, se puede encontrar fcilmente utilizando el
siguiente procedimiento
11
Calcular la Funcin Normal Acumulada:
F(Z) =
lT
Obtener el valor de Z en una Tabla de Distribucin Normal
J-
Entonces:
X = X+2.5"
Hidrologa
/ David
Cedeo
279
Funcin de l)istribucin Normal Acumulada
F(z)
I _,, ""nl , )aw

)
0.0t 0.0 0.t

0.1 0.1 0.1
0.02
001 0.i lt0 0.55t7 0.59t0
004
0.5
0.0i
0.J t99
0.06
0.5000
0.5198 0.5791
0.6
l]9
0.6554
0.69 r5
0..i
0.6 0.1 0.8 0.9
1.0
0.1257 0.7580
0.7631
0.8ri9
0.84r 0.8641
Ll
t.2
0.8849
0.90J 2
LJ
1.4
0.9r92
0.9132
0.9452
t.5
t.6
1.7
0.9i54
0.9641
t.8 t.9
2.0
2.1
0.97r1
0.9112 0.9821
0 986r
2.2
2.1 2.4 2.5
0.9891 0.9918
0.9918
0.99.r3
2.6
2.1
2.8
0.99s
0.99 t4 0.9981
O.998 t'
?.9
1.0
0.7910 0.7919 0.1961 0.8186 0.8212 0.82 !8 0.3418 0.8461 0.8485 0.865 0.8686 0.8708 0.8869 0.8838 0.8907 0.9049 0.9066 0 9082 o.9zo7 0.9222 0.9:36 0.9145 0.915? 0.9170 0.9463 0.9414 0.9484 0.9564 0.9571 0.9582 0.949 0.96i6 t\.9664 0.9719 0.9126 0.9132 0.9't't8 0.978J 0.9788 0.9826 0.9810 0.9814 0.9864 0.9868 0.9811 0.9896 0.9898 0.990t 0.9920 0.9922 0.9925 0.9940 0.99 I 0.9941 0.9955 0.9956 0.9951 0.9966 0.9961 0.9s68 0.99?5 0.9976 0.99"t1 0.9982 0.9982 0.9e8.1 0.9987 0.9987 0.9988
0.76
0.504c 0.5080 0.5138 0.5.;'8 0.J812 0.5871 0.6211 0.6255 0.6591 0 6613 0.6950 0.6985 0.7291 0.71t4
0.764.1
0.5557 0.594S 0.6700 0.7054 0.7389
r60
0.6i93
0.6664
0.70
0.611!
0.j596 0.i987
0.6168
0.61i6
0.1088
0.521e 0.5279 0.5.16 0.56?5 0.6026 0.6064 0.6406 0.6441 0.617: 0.6808
0.7
0.51t9
0.57 0.6
0.5i59
0.5751
t4
t0_1
0.6rJl
0.6480
0.6$44
0.6ir7
0.6879
'9 0.7357
0.1613
o.iii; 0.799i 09264

0.3508 0.8729 r).8q25 0.9099
0.925
0.i422
0.7131
0.8021 0.3289
0.351I
0.8749
0.3944 0.9 t r5 0.916s 0.9394

0.9505
0.9.182
0.9611 0.9678 0.9718 09744 0.9193 0.9198 0.9818 0.9842 0.9875 0.9878 0.9904 0.9906 0.9927 0.9929 0.9945 0.9946 0.9959 0.9960 0.9969 0.9970
rr.99.17
0.9978 0.9984 0.9989
0.959r
0.9495
0.9599
0.9984 0.9988
0.7164 0.1194 t).7821 0.78_i: 0.8051 0.8078 0.810 0.81i l 0.81l5 0.3J0 0.8i65 0.8189 0.85i 0.8577 0.8599 0.36: I 0.8770 0.8790 0.88 r0 0.33i0 0.8962 0.8980 0.8997 0.901i 0.9tJ I 0.9147 0.9162 0.9177 0.9219 0.9292 0.910 0.9,119 0.9406 0.9418 0.9429 0.94,1t 0.q5ti 0.9525 0.9515 095J: 0.4608 0.96 t6 0 q25 n.glj 0.9686 0.969t 0.9699 0.910 0.9750 0.9156 0.9?61 t\.9161 0 qRoJ 0.s80R 0.s8 t 2 o 98 I 7 0.9846 0.9850 0.9854 0.9857 0 9RR I 0.98R4 0.'rr87 0 o!9t) 0.9909 0.991 I 0.99tJ 0.991 0.991r 0.9912 0.9914 0.99,\6 0.9948 0.9945 0.9951 0.9ci2 0.996 r 0.9962 0.996J 0.996.r 0.997t 0.9971 0.9973 0.9914 0.99?9 0.9979 0.9980 0.9q81 0.9985 0.9935 0.9986 0.9986 0.9989 0.9989 0.99q0 0.9990
D1 0.7t57 0.7190 0.ltt4 0.7454 0.7486 0.7i I7 0.7:tS

Eiemplo
220
27:
Distribucin Normal
Emplear la distribucin normal en combinacin con los datos de

caudales mximos registrados en la estacin de aforo (utilizados en los
ejemplos anteriores desarrollados en esta seccin) para calcular:
a) b)
El caudal con perodo de retorno de 100 aos (Q,oo).

La probabilidad de que el caudal ser menor igual a 10,000 pies3/seg
y su respectivo perodo de retorno.
Solucin: Los parmetros de la distribucin normal estimados en trminos de los momentos de la muestra son:
= O = = sa =
4,744p3/s
3,311 p3
ls
y la tuncin de densidad de probabilidad
es:
f(e\ =
a)
---="*nl-t0-9r' ,sn" sotfzr"

L
T:
-.1
]
-*<e<*
100 aos:
Caudal para un perodo de retorno
F(Q,oo)
= 1- L = t- ,# = o.ee
)41
Interpolando en las tablas de distribucin normal acumulada, tenemos:
Z =
2.326
Por lo tanto la magnitud del caudal es:
Qroo
= Q+Z'Sn =
4,144+2.326x3,3
11 =
11,850p3/s
b)
Variable Normal Estandarizada:
z = O-O so
1o'ooo-4'144
3,3
1l
1.769
Utllizar las tablas de la distribucin normal acumulada para obtener:
F(Z) = F(1.76e) = 0 9615; por consiguiente,

P(8<1O,OOOp3/s)
la probabilidad es:
0.915
= F(Z) =
y el perodo de retorno para un caudal Q
10,000 pies3/seg es:
_l I_
Observacin:
r-F(z)
1 1-o.els
_ ^Zb AnOS
La Distribucin Normal no se ajusta satisfactoriamente a los datos de

caudales mximos registrados (histograma de densidad)
debido a la asimetra
de la distribucin emprica.
J
H4
{!
Y. az
9A ur ti< 44
a
i
&.. (,g
E <t: d
<;J
fJ
Zi! o4. co
tri =d
E?
f-l Z , ij
?,
Vz
UX )2.
o
.
!E
Q
++{qq4sB HHgHgHfi EREEggHE
223
DISTRIBUCION LOG NORMAL

Considere un clculo hipottico de escorrenta, en el cual la escorrenta
es igual al producto de funciones de diferentes factores aleatorios, tales como
precipitacin, rea tributaria, coeficiente de prdidas, evaporacin, etc. En

general, si una variable aleatoria
es el resultado del producto de una gran
cantidad de otras variables aleatorias (mecansmos de multiplicacin),

entonces la distribucin de los logarnos de x ser aproximadamente normal,
ya que el logaruno de
estar cnstirudo por
h sna de los logartmo de
los factores contribuyentes. En hidrologa, es muy fcil concebir muchos

factores contribuyentes a la escorrenta que son aleatorios y de los cuales
existe muy poca informacin determinstica; por lo tanto, un mecansmo multiplicativo para la escorrenta puede ser una suposicin razonable. En
general, se puede considerar que una variable aleatoria tiene una distribucin
log normal si el logartmo de la variable aleatoria est distribudo

normalmente; es decir que
si I :
log
0)
tiene una distribucin normal,
entonces
tiene una distribucin log normal. La funcin de densidad de
probabilidad para la distribucin log normal es:

log,o (e)
f(x)
-0 "-rl
Fr)'
, o ,lin
2 o2 v
>0
donde
log,o
6) y e es la base de los logartmos nafurales; observe

y
se
que la funcin de densidad de probabilidad para la distribucin log normal

est limitada a la izquierda por el valor cero, tiene una asimetra positiva
.'t a
puede utilizar fcilmente debido a su relacin con la distribucin normal, ya
que la variable transformada tiene la siguiente funcin de densidad de

probabilidad normal:
f(y)
exo
o" rfz
'l
- r,..1' l f-tv | '

2
o,'
para la cual existen tablas que nos permiten obtener directamente la funcin
de densidad acumulada, ya que la distribucin normalse utiliza ampliamente
en estadstica.
La asimetra
coeficiente de
de la distribucin log normal es una funcin del
variacin CV
p": CV3
3 'CV *
pero la asimetra de la variable transformada Y

se espera, debido a
log (X) es cero, tal como
que
est distribuida simtricamente.
A pesar de que la funcin de densidad de probabilidad de X (variable

original) se puede derivar fcilmente, esta funcin se requiere en muy pocas ocasiones. En lugar de eso se encuentran los momentos
de
I:
log
(E
se utiliza la distribucin normal para la variable transformada
I.
Esta
transformacin se puede hacer de dos maneras, ya que se pueden usar

logartmos naturales de base 10. El mtodo de momentos requiere que el
momento de los datos originales sea calculado y relacionado con el momento
de los datos transformados. Cuando se utilizan logartmos naturales en la
transformacin:
ln (X) , stas relaciones son:
22s
CV2
o2
exp
(o,2
u2
1)
t^=
encontrar Fy
/
exp
F..
"'\ * --2-
2)
Las ecuaciones anteriores se pueden resolver simultaneamente para
arz. Si se utilizan logartmos de base 10, las frmulas
anteriores se deben modificar reemplazando la bise de los logartmos

raturales (e
2.718) con logartmos de base 10. Observe que el logartmo
del promedio no es el promedio de los logartmos; por el contrario,
ly = log(X,)
donde X, es la mediana de
X.
El mtodo de momentos para la estimacin de

es
los parmetros conserva el momento de los datos originales, 1o cual

necesario para la simulacin de caudales en ros.
Otro mtodo es el de convertir el promedio y Ia variarza a logaritmos
naturales. Este segundo mtodo es ms comn
est relacionado con
estimaciones de probabilidad mxima para los parmetros
de
I.
En este
caso, primero se calculan los logaritmos de todos los caudales y luego se

obtiene el promedio y la variarva de los datos transformados; si los momentos
de los datos transformados son calculados por ambos mtodos, la semejanza
de los resultados se incrementar al aumentar el tamao de la muestra.

Ejemplo
z6
28:
Distribucin Log Normal
Utilizar los datos de caudales mximos registrados en la estacin de

aforo
aplicar la distribucin log normal (usando logartmos de base l0)
para calcular:
a) b)
El caudal con un perodo de retorno de 100 aos.

La probabilidad de que el caudal sea menor o igual a 10,000 pies3/seg
y su respectivo perodo de retorno.

Adems, comparar los parmetros estimados utilizando los momentos de los datos originales y los obtenidos con los datos transformados.
Solucin: Los parmetros de la distribucin log normal estimados utilizando los momentos de los datos transformados
Y=
log,o
(p)
son:
f =
'Sr =
3463
0 424
y la funcin de densidad de probabilidad
es:
fQ\
2.3 Q Sy
\ft;
.-li#]
T
a)
Caudal para un perodo de retorno
100 aos
F(r,oo) = 1-
L = 1---l
0.99
22'7
Interpolando en la tabla de distribucin normal acumulada, tenemos:
z =
Iroo = i*Z'5,
2.326
Por lo tanto, la variable transformada
l,oo =
logro (0,00 )es:
= 3.463+2.326x0.424 =
4.449
y la magnitud del caudal se puede obtener de la siguiente forma:

Q
,oo = 16(rroq) =
16(a
aas)
= 28:80
p3 I s
b)
Transformacin del caudal
0:
10,000 p3/s:
Y = log,o(10,000) = 40
Normalizacin de la variable transformada:
- = y-y Z ,"
4-3.463
1.267
Utllizar la tabla de la distribucin normal acumulada para obtener:
F(Z) = F(l.267) = 0.8975; por consiguiente,

P
la probabilidad es:
s) = F(z) = 0.8e75 y el perodo de retorno para el caudal Q : 10,000 p3ls es:

(Q <
ro,ooo p3
I
T =
I 1-F
1-0.8975
= 9.8 aos ' lO aos
228
Comparacin del Valor de los Parmetros: Coehciente de Variacin:
CVt
a_ Fx
':
,.9x
3,3 10.82
0.799
4,t43.94
Coeficiente de Variacin al cuadrado:
CVrt =
Yarianza:
exp
(of - t' "
exp(S]'?- l)
S: = tn(Cl/xl)*r = ln[(0.799)2]+r =
Desviacin Estndar:
0.ss1
sf
Media Valor Promedio:
O.742
Px =
exp
I
\
Fv
+) " "*nlr * +)
h(4,143.e4) -
=X
Por lo tanto, tenemos:

,s2 Y = tn(X) - -1!= 2
o'551 =
2
8.054
Cambio de Base de los Logartmos:

logro
(m) =
logro
@) ^
ln(m) =
0.4343
ln(m)
ln (m)
2.3026
aa
Por consiguiente, el Valor de los Parmetros por el Mtodo de Momentos

(utilizando logartmos de base 10) es:
y=8'054=3.49s
2.3026
s. ' Resumen:
2.3026
0742 =
0.322
Valor del
Parimetro
Y
Mtodo de Momentos 3.498

0.322
Usando Datos Transformados
3.463 0.424
,s/
Observe que los valores obtenidos para los parmetros por ambos mtodos
son parecidos; pero el mtodo de datos
transformados
Y=
logro
(p)
es
el ms comun. Sin embargo, la diferencia es notable al calcular la magnitud
del caudal con un perodo de retorno de 100 aos:
f,oo = l*Z'Sr
Q
3.498
+2.326x0.322 =
17,687 p3 ls
4.247
,oo = 1(v'oo) = 16(+zn7) =

^EJ Q
< z E9 a: Fz
J(J
OF
aA
!,
F,l
O o c)
CL
8<
a o l o
a o
c; rrJ
^ z. oFf
En
Fr t,r gEl o A
.l
FZ zr
EEl
|i z. zo
(J
U
e Z a\6 f-l
fe
l-r
rU
3 o C] o
3
tJ
o (tr
3 o o o q
t.J
3 Lr o q
o tr.J o (f c) o
,i
o <l
trl
o o o c) o
(oI
23a
DISTRIBUCION GAMMA
Esta distribucin es utilizada ampliamente en hidrologa simplemente
debido a su forma y propiedades matemticas bien conocidas; entre ellas tenemos las caractersticas de estar limitada por la izquierda
y un sesgo
positivo, 1o cual es conveniente para eventos extremos de inundaciones

(adems, la distribucin garma tambin se puede utilizar con un sesgo
negativo para eventos de sequas). La distribucin Gamma puede ser de dos parmetros (Gamma - 2) de tres prametros (Gammi - 3), la cual recibe el
nombre de Pearson Tipo

el caso
de Pearsontipo
3. A pesar de que los coeficientes: ,1, y p en 3: 1, p y e requeridos para definir la funcin de

media
.r,
densidad de probabilidad son funciones simples que utilizan la
desviacin estndar
y la asimetra T; es muy comn en aplicaciones
hidrolgicas evaluar la funcin de densidad acumulada con factores de

frecuencia, evitando de sta manera la integracin de la funcin de densidad.
Existen tablas para los factores de frecuencia
K, los cuales son funcin de la
asimetra corregida
CJ y del perodo de retorno Z; por 1o tanto, para
estimar la magnitud del caudal con un perodo de retorno deseado, se calculan
los momentos de los datos: media g y desviacin estndar sn y luego se
utiliza la siguiente frmula:
Q, = Q+K(CS.T)'Sn
Este procedimiento es satisfactorio para los perodos de retorno
mostrados en la tabla, pero no es apropiado interpolar en la tabla para otros
z 1l
intrvalos de recurrencia.
El procedimiento correcto sera efectuar
una
interpolacin grfica en el diagrama de caudales en funcin de perodos de
retorno. Adems,
este es el procedimiento para el problema inverso de
encontrar la probabilidad (y por consiguiente, el perodo de retorno)

correspondiente a un caudal de una magnitud dada. Alternativamente, los dos tres parmetros de la distribucin Gamma
se pueden calcular
y a continuacin
se utilizan las tablas de
la funcin de
densidad acumulada que aparecen en algunos librcs de estadstica para determinar la probabilidad.
La distribucin Gamma de dos parmetros (Gamma - 2) corresponde
al establecimiento del lmite a la izquierda de la funcin de densidad en el valor cero
(. :
0); en ste caso no se requiere calcular la asimetra de la
muestra, ya que se puede utilizar la siguiente relacin para estimar ese valor:
cs = t = 2.cv2
Por lo tanto, podemos obtener el coeficiente de asimetra requerido para
encontrar el factor de frecuencia
de la siguiente manera:
cs
,(+
Facto de Frecuenca K para las Distribucnes: camma, Pearson Tpo 3 y Log Pearson Tpo 3.
lhtrvelo d Rcurrnc cn
Ao!
3.0 2,9 2,4 2.7 2.6 2.5 2,3
2.2
2.1 2.O
1.9 1.8 1.6 1.5 1.4 1.3 4.2

'1.1
't.o
0.9 0.8 o.7 0.6 0.5
0.4
0.3 o.2
0,1
o.o
-0.1
4,2
-0.3
4,4
-0.5 -0.6
4.7
J).8
.9
-1.0
-1.1
-1.2 -1.3 -1.4
.1.5
-1.6 -1.7 -1.8 -1.9
-2.O -2.1
-2.2
-2.4
-2.5 -2.6 -2.7 -2.9
3.0
4.667 4.690 4.714 4.740 4-769 4.799 4.832 4.867 4.905 {.9,15 {.990 '-1.037 -1.087 -1.14 -1.197 t_256 f .318 .3a3 -,1.4d9 -1.518 -1.588 -1.660 .733 -1.806 -l.aao -1.995 -2.029 -2.104 -2.17a -2.252 -2.326 -2.400 -2.472 -2.544 -2.615 -2.686 -2.755 -2.424 -2.491 -2.957 3.022 3.087 3,149 ,211 3.271 .330 3.388 -3.444 3.499 3.553 3.605 3.556 3.705 3.753 .800 3.845 3.489 3,932 -3.973 -4.013 -4.051
-0.396
.390 .384
).376
o.48
0.t50
1.195
2.277
1.210
1.224
2,275
2.272 2.267
0.479
0.,1!t9
3.134 3.114 3.093

3.071
4.368 4.360
.0.351
-0.341
0.518
0.537
0.555 0.574 0.592 0.509 o.627
0.513
1.234 1.250 1.262

1.274
2.262 2.256
3.044
3.023
2.28
"0.330
t2a4
1.294 1.302
1.310
2.26
2.230 2.219
2.207
.319
-0.307
2.997 2.970 2.942 2.912

2.a41
4.294
4.2A2
l.3la
1.324
1.329
.268
4.254
4.2&
4.225 4.210
-0.195
0.560 0.675 0.690 0.705 0.719
2.193 2.179
2.163
2.48
2.415 2.7aO
2.16
1.339 1.340 1.341
2.124 2.104
2.O47
{.1a0 {.164
o.732 o.745
0.754 0.769 0.780 0.790 0.800 0.808 0.816
2.743 2.706 2.666 2.426

2.585 2.542 2.498 2.453
4.r
4,132
4.116
1.3&
1.339 1.336
og9 -0.0a3 -0.066 -0.050 -0.033
2.056 2.O43 2.018 1.993 1.967 1.939

1.9'10
2.67
2.359
2.31,1
o.424
0.830 0.836
1.317 1.309
1.301
1.440
L849
1.818 1.785 1.751
2.261 2.211
2.159
2.107 2.O54 2.000 1.945
4.017
0.000 0.017 0.033 0.050 0.066 0.083 0.099 0.116 0.132 0.148
1.292
1.242 1.270
o.442
0.446 0.850 0.853 0.855 0.856
1.716
1.880 1.543 1.506 1.567
l.25A
1.26
1.231
L890 La34
1.777
1.216
1.200
0.857 0.857
0.856 0.854 0.852 0.8rl8 0.844 0.83a 0.832 0.825 0.817 0.808 0.799 0.788 o.777 0.765
1.52a
1.720
1.563 1.606 1.549
t.183
1.166 1.147
.448 ,1.48
1.&7
1.366
0.164
0.1a0 0.195
1.124
1.107 1.086 1.064
1.O41
1.492
1.379
1.324
0.210 0.225
1.2&
1.194 1.157
1.324 1.270
1.217 1.155 1.116 1.069 1.O23 0.940 0.939 0.900
0.2&
0.254
0.268
l.ot8
0.994 0.970 0.945 0.920 0.495 0.469 0.444 0.a19 0.795
o.771
l.ll5
't.075 'f.035
0.996 0.959 0.923 0.488 0.455 0.423 0.793
0.2A2 0 294
0.307 0.319 0.330
0.341
o.752
0.739
0.864
0.830 0.794 0.764
0.351 0.360 0.36a 0.376 0.384 0.390 0.396
o.725
o.711 0.696
0.681
o.747 o,724
0.764
0.738 o.712 0.683 0.656
0.656
0.651
o.702
0.681
o.740 o.714
0.649 0.666
4,051 4.970 4.013 4.904 3.973 4.A47 3.932 4.783 3.449 4.714 3.74{t 4.652 3.aoo 4.584 3.753 4.515 3.705 4.444 3.656 4.372 3.605 4.294 3,553 4.223 3.99 4.147 3.444 4.069 3.344 3.990 3.330 3.910 3.271 3.a28 3.211 3.745 3.'149 3.661 3.087 3.575 3.022 3.89 2.957 3.&1 2.a91 3.312 2.424 3.223 2.755 3.132 2.686 3.041 2,615 2.949 2,544 2 456 2.472 2.763 2.&O 2.670 2.326 2.574 2.252 2.82 2.174 2.33a 2.,104 2.294 2.029 2.201 1.955 2.104 l.aao 2.016 1.806 1.926 1.733 1.837 1.660 1.749 l.saa I 664 1.514 1.58.1 ,1.449 1.501 1.343 1.424 1.31A 1.351 '1.256 1.2a2 1.197 1.216 1,1& 1.155 1.087 1.097 1.037 1.044 0.990 0.995 0.9,16 0.949 0.905 0.907 0.867 0.869 oa32 0 933 0.799 0.400 0.769 0.769 0.740 0.741 0.714 0.714 0.690 0.690
0.636
0.660
Fuentef Haan, C.T., 1377, Slaai.aical Methods in Hlldrology,lowt g.ate ljnvrsity Press, Anes.

Ejemplo
234
26:
Distribucin Gamma
Utllizar los datos de caudales mximos registrados en la estacin de

aforo para calcular la magnitud del caudal con un perodo de retorno de 100 aos, utilizando las siguientes distribuciones de 2 y 3 parmetros:
a) b)
Gamma - 2
Pearson Tipo 3
Adems, describir la funcin de densidad de probabilidad y sus parmetros.
Advertencia: El problema inverso de encontrar et periodo de retorno para un
caudal de 10,000 p3/s requiere una grfica de caudales en funcin de Ia

probabilidad cuando se utiliza la distribucin Gamma; por lo tanto, no lo
vamos a calcular en esta ocasin.
Solucin:
a)
Distribucin Gamma - 2:
Media y Desviacin Estndar de la Muestra
.
O =
4,144 p3 /s
So = 3,3llP3/s
Coeficiente de Asimetra:
cs = ,lLl
\2
loJ
= zf ,rr'1' =
\+'t+t)
2(0.7ee)2
1.277
I I I I I
Hidrologa
David Cedeo
a?tr
Factor de Frecuencia (obtenido interpolando en la tabla):
K = K(CS=r.227,r=loo) =
3.1e7
Magnitud del Caudal para un Perodo de Retorno T
100 aos:
Qtoo
-- Q+KxSo
Q
= 4'144+3'197
14,730 P3 ls
x3,111
l
I I I I I I
Funcin
ro, =
de
Densidad
de
Probabilidad para
la
Distribucin
Gamma - 2:
e'xQ \9 O9-t z ; f@) = '" f(B)

Parmetros de la Distribucin Gamma - 2:
Q>o
L =
/
O s.'
12
(3.31
=
1 )2
3.280,. lo-a
I t I I I I
p = l_9_l - [ \ \ sn /
Observacin:
\2
4,t44
\r
3,31r /
l'=
1566
La funcin de densidad de probabilidad para la distribucin
Gamma - 2 se ajusta de manera adecuada a la distribucin emprica de los caudales mximos registrados (histograma de densidad), tal como se puede
apreciar en la siguiente grfica.
'q
riE
a_
'
z f
to
ac
)d
<9
=u1 !r= F r-t
I LJ
'f
8oa
d,lr
Ov
e t n
a o
f
2.8
l?
)6Fl
o o e
ao
F- r.l vl a
l.1 (4
2A tl r.t 'a, A o7. <U 4?. <5 Ff{ :

Q
dz. t- f
ti FFSIEF
(ol
gEEE
rr uj r.
vr,.rd
q88
.t.
231
b)
Distribucin Pearson Tipo 3: Media, Desviacin Estndar y Coeficiente de Asimetra Corregido,

utilizando los momentos de la muestra:
0 =
So =
4,144 p3 ts 3,311 p3
ls
cs =
1.981
Factor de Frecuencia (obtenido interpolando en Ia tabla):
K = rK(CS=1.981,2=100) =
3.595
Magnitud del caudal para un perodo de retorno
Z:
100 aos:
Qroo
= 8"K'Sn
Q too
= =
4,144 +3.595x3,3
ll
.,
16'050 P3 ls
Funcin de Densidad de Probabilidad para la Distribucin

Pearson Tipo 3:
f (Q)
LF (O
- eF-r ,-L(Q-.) ;
f(p)
Q>e
I I
I
I
I
I
Hidrologa i David Cedeo

Parmetros de la Distribucin pearson Tipo 3:
??a
B=
(+)'= (t+)'=,0,e3
_@
sn
L=
I
3.049
a " l0-
- t, t[l =
4,t44
-:,:u / r-orx =
801.35
p3ls
Observacin:
En la siguiente grfica se muestra el ajuste de la funcin de densidad de probabilidad para ra distribucin pearson tipo 3 a la distribucin emprica
de los caudales mximos registrados (histograma de densidad)
en la misma
.
podemos apreciar el lmite a la izquierda
: e>e"
B0O
p3 ls
J
-l
..
a-
e x
<o .u
7,
rpr AZ.
vl2.I
<J U 2 tr a,o l.l
I.)J FF] d9 (t
H \-
o* o
a
J
7.4
eo E< ()
Uc ;r? q
ir
2.
Z
<^ &tt r- v)
-4
l-t
<U lZ P
U
?7
3 o o
ul tl
3
o a o o
o q
o ll o o o
ci
r.1
3 o
3
lrl
f
u.t
CJ
c)
t].t
o <J
tJ
(f
o c)
".1
o C]
o o
c)
r-t
ci
o
a
(oli
l
I I
240
3
DISTRIBUCION LOG PEARSON TIPO
Cuando se aplica la distribucin Gamma de tres parmetros a los
I I
I
logartmos de los datos de variables aleatorias (datos transformados),

generalmente se
le
denomina distibucin Log Pearson tipo
3.
Esta
distribucin tiene un papel muy importante en hidrologa debido a que ha sido recomendada por el U.S. Interagency Advisory Comittee on Water Data
(1982) como la distribucin de probabilidad que debe ser aplicada a los

caudaies de inundacin.
I
I
La forma de la distribucin Log Pearson tipo 3 es bien flexible debido
a que utiliza tres
parmetros
su uso es completamente similar a la
I I
distribucin Log Normal discutida previamente; sin embargo los momentos

de los datos originales y transformados no esln relacionados en esta ocasin.
Por el contrario, los datos son transformados tomando logartmos (de base
base 10) de la siguiente manera:
y = loe(Q) y la distribucin
Gamma - 3
I I
se aplica exactamente como se hizo en la seccin
anterior. Esto significa que
las magnitudes de los caudales se pueden calcular directamente para los

perodos de retorno mostrados en la tabla, pero el problema inverso de
determinar el perodo de retorno ( la funcin de densidad acumulada)
l
I
correspondiente
grficamente.
una magnitud dada solamente se puede resolver
I
I
t I
I I I
Hidrologa / David
Ejemplo
Cedeo
Distribucin Log pearson Tipo 3
24L
27:
Utilizar los datos de caudales mximos registrados en la estacin de

aforo para calcular la magnirud del caudal con un perodo de retorno de 100
aos aplicando la distribucin Log Pearson tipo 3; dems, describir la funcin
de densidad de probabilidad
y los parmetros de la distribucin.
I
I
)o luc
ron I
-r*nsformacin
de los Datos (logartmos de base r0):
I
I
los datos transformados:
y =
log,o
(0)
Media, Desviacin Estndar y Coeficiente de Asimetra Corregido para
I
I
I I
I I
I
Y =
3.463
0.424
sv =
c,s =
-1.117
Factor de Frecuencia (obtenido interpolando en la tabla):
K = K(CS=
-1.117,
f = 100) =
1.504
Variable Transformadai I,oo = logro (0roo )
tI
I,oo = f *,KrSy =
3463 +1.504x0.424
4.101

Por lo tanto, la magnitud del caudal paraT
242
100 aos es:
Qro, = 16rtt =
16(a
lot) =
12,6.10 P3
/s
Funcin de Densidad de Probabilidad para la Distribucin
Log Pearson Tipo 3:
f (o)
Donde
0 f (p)
; Q > e',
se
= log,o(Q)
laconstantenumrica ln(10) = 2.3
debe al cambio de base de los loeatmos.
Parmetros de la Distribucin Log Pearson Tipo 3:
p (+)'
l=
r.rtt
. '
\z
3206
,lT
2.3
sY
2.3 (0.424)
1.836
= z:ql -s, 1r- I =
23(3.463-0.424r|;.zoe
>=
6.227
Lmite de la Funcin a la Izquierda:
O > e'=
exp(6.227)
506.22 p3 ls
-,
6f FH
<o 6a
AY .tz ?= U
Fd g -.
hla
:f
A
Fl
<
f. 2 A
rJ
2Fl
r. a
uf
A
<A 4z
HHg q38
xlJo
!tt<.(i
+sBs
HHHH Hqsge
(.rF_-jo
++se8
(o)
244
DISTRIBUCION GLIN{BEL (Valor Extremo Tipo I)

Est funcin de densidad de probabilidad surge de la teora de valores
extremos desarrollada por Gumbel (1958). La distribucin tiene una forma

adecuada; pero al igual que la distribucin normal, carece de lmites en ambos
extremos; sin embargo, la posibilidad de un valor menor de cero (por

ejemplo, caudales negativos) es de muy poca preocupacin en su aplicacin.
Unavez ms, no se requiere la funcin de densidad de probabilidad, ya que

la funcin de densidad acumulada tiene la siguiente forma particular de una doble exponenciacin:
F(x) = erp{-exp[-a(x - z)]] ;
-c<x<6
tablas
La funcin anterio puede ser evaluada fcilmente sin necesidad de y en la misma, el parmetro d es un factor de escala y el parmetro
promedio r y la
es la moda. Estos parmetros estn relacionados con el
desviacin estndar
por medio de las siguientes expresiones:
u=P
"G
Las cuales proporcionan una forma de estimacin de los parmetros por el
mtodo de momentos. Se ha demostrado tambin que no es necesario efectua
l
I
I I
I
I
I

correcciones para tomar en cuenta el nmero de datos
245
en Ia evaluacin de
c y u.
La asimetra
de la distribucin de Gumbel tiene un valor
constante igual a 1.1396.
El argumento del exponencial interno:

frecuencia
y = d (x - u) se denomina
de retorno
la variable reducida y debido a que la misma es una funcin nica de la
F=
|-
1/T (laeual est relacionada con el perodo
I),
estr variable se puede :utilizar para obtener una grfica de la distribucin
en forma lineal utilizando factores de escala adecuados.
Si la funcin de densidad acumulada se utiliza en combinacin con las
relacones de momento para obtener
la variable
en trminos de
Ia
frecuencia
del perodo de retorno 7, se puede derivar facilmente un
factor de frecuencia K:
K = -0.7797Y 0.5772- o77e.t,nl,"(+)

]
K=
-0.77e7
05772
- o77s7"l'"(,=)
p)/
o
El factor de frecuencia: K = (X
se puede
utllizar para obtener

ax
.
directamente el valor de la variable X; es
decir: X = v * K
Si la funcin de densidad acumulada se evala para el valor promedio:
X=p
se encuentra que su magnitud
perodo de retorno
es: F(f.) = 0.57;por 1o tanto, el para la media ,r es: T, = Il(1 - F) = 2.33 aos. St
existe una grfica de caudales en funcin del perodo de retorno, algunas
246
veces se usa la informacin anterior para estimar de manera rustica el valor
promedio de los caudales, suponiendo que los mismos obedecen a una

distribucin de Gumbel.
9=t1
Utilizar
Distribucin Gumbel
nuevamente los datos de caudales mximos registrados en la
estacin de aforo (serie anual mxima) en combinacitTcon la Distribicin de
Gumbel para encontrar la magnitud del caudal con un perodo de retorno de
100 aos; adems, encontrar la probabilidad de que la magnirud del caudal

sea menor igual a 10,000
p'ls y su respectivo perodo de retorno'
Solucin:
Media y Desviacin Estndar de la Muestra:
O =
So =
4,744 p3 /s 3,377 p3/s
Funcin de Densidad de Probabilidad:
fQ) = d exp\-*tO - u) -exp[-c tO - r)l\

Parmetros de la Distribucin Gumbel:
<Q<
tnG
3JnG
3.874
^ lo-a slp3
I
I
I'
I
u = e- 0'5772 = 4,144Factor de Frecuencia
0.s772
3.874x IO-a
2,653.9 p3 /s
K para Z -
100 aos:
K = -0.450 -o77s7rl,,,lrqg)l = 3.137 ssl]
[ \
Magnirud del Caudal:
Qroo
= Q*KSB = 4,144+3.137 x3,3 1l =
14,5j0p3/s
Funcin de Distribucin Acumulada:
F(Q) =
Probabilidad:
exp{
-"*p[- "@ -
u)Jl
F(10,000) = exp {-"rp[ -3.874x l0-4(10,000 _ 2,653.g)]l

p (Q < 10,000 p3ls)
Perodo de Retorno:
= '(10,000) =
0.943
{l
-ll
I
| -F
I -
0943
18 aos
(JO Z F.(
t
o ()
a
u)X ZY
.
2t( dbr
I-
r-t
)D
oz
9i'
O it 4 ri
Zj
o ?
Xoi 'th J
ce 4.
'\
ca
r'J
zc Y 9A v
EH
()
1i
F+qqqqqqs ur W ooooooc';;r.u
eooooo9 9oooo^ roaxo
rlj
uJ
ru
rjJ
.-+.
;'-,j:;;;x-c
(ol
Hidrologa / David Cedeo Comparacin de Resultados

Caudal Pronosticado para 100 aos
249
Distribucin Normal Log Normal:

Datos Transformados Mtodo de Momentos Gamma - 2
Pearson Tipo 3
T:
Asimetra Calculada
1.98
11,850 p3ls
28,130 p3ls 17,690 p3ls 14,730 p3ls 16,050 p3ls 12,610 p3/s 14,530 p3/s
- 1.12
1.28 1.98
Log Pearson Tipo 3

Gumbel
t.l2
1.98
Observacin:
Los caudales pronosticados para un perodo de retorno de 100 aos

varan al utilizar cualquiera de las 6 distribuciones de probabilidad estudiadas
en sta seccin, las cuales fueron aplicadas a la misma serie anual mxima (caudales registrados en la estacin de aforo). Es importante sealar que el mejor ajuste se debe determinar grficamente comparando la funcin terica
de densidad de probabilidad con la funcin emprica (histograma
de
densidad). Adems, podemos observar que la asimetra tambin afecta el

ajuste de los datos a una funcin de densidad de probabilidad y en el caso de
la distribucin normal, la funcin terica es simtrica, pero el mtodo de

momentos estima una asimetra corregida diferente para los datos.
REFERENCIAS
and Floodplain Anafiss' Bedient, P. B. y N. C. Hubert ,1gg2' Hydrology Reading' Edi.tion, Addison - Wesley Publisghing Co''
Second Massachusetts.
Chow, V. T., D. R. Maidmenty
L'W'
Mays' 1988' AppliedHydrology'
McGraw - Hill Publishing Co', New York'
Linsley,R.K.,M.A.Koh]eryJ'L'H'Paulus'1982'Hydrologyfor
New York' Engineers,Thid Edition, McGraw - Hill Book Co" Handbook of Hydralogy '
Maidment, D. R. (edito)
1'993,
McGraw - Hill Inc., New York'

England' 1979, Hydrology, Pergamon Pess' Oxford'
Raudkivi,
A'.
J.,
Shaw, E.
M.,lgg4, Hydrology in Practice' Third Edrtion'
Chapman & Hai1, London. 1989' Introduccon to
Viessman, W., G.
L' Lewis y J' W' Kapp'
New York' Hydrology, Third Edition, Harper Collins Publishers'

Hidrologia - David Cedeño

Cargado por

Copyright:

Formatos disponibles

Hidrologia - David Cedeño

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Hidrologia - David Cedeño

Cargado por

Copyright:

Formatos disponibles

F

TNI\'ERSIDAD TECNOLOGICA DE PANAMA

Departamento de llidrulica, Sanitaria y

Ing. Dad Cedeo B.

a) HumedadAtmosfrica . . . b ) Cambios de Fase c) Cantidad de Agua Precipitable d ) Causas y Mecanismos de Formacin de la

Anlisis de Datos de Precipitacin

Anlisis de Doble Masa

Precipitacin Promedio sobre una Regin

Infiltracin Infiltracin . Aforos

Otros Mtodos para calcular la

Medicin del Caudal:

ANALISIS DE PRECIPITACION - ESCORRENTIA

distribucin del agua sobre la

hidrologa incluye los procesos fsicos, qumicos y las reacciones biolgicas

del agua en los ambientes naturales y aquellos construidos por el hombre.

relaciones con los

patrones climticos, tipos de suelo, topografa y otros factores geolgicos, es

dificil establecer fronteras entre la hidrologa y otras ciencias de la tierra tales

relacionados con la tierra;

debido a las implicaciones en el medio ambiente de las inundaciones y

hidrologa incluye las prcticas

de agua en las civilizaciones Egipcias y

(Mesopotamia) establecidas en el Medio Oriente a

largo de los ros Nilo,

Tigris y Efrates y en la China a lo largo del Ro Amarillo. Las excavaciones

arqueolgicas en estos lugares muestran evidencias de estructuras hidrulicas

agua. Por ejemplo en el ao 4000 A.C. se construy un canal para

Los filsofos griegos fueron los primeros investigadores de

hidrologa. Aristteles propuso la teora de que la humedad del aire se

Homero sugiri la idea de que exista un mar

subterrneo que era el origen del agua superficial.

los sistemas de acueductos Romanos (97 D.C.) basados en la

Estos teoremas continuaron hasta el Renacimiento

cuando Leonardo Da Vinci descubri una relacin adecuada entre rea,

La primera medicin registrada de precipitacin y flujo superficial

El astrnomo ingls Halley (1636 - 1742) utlliz una bandeja pequea

para estimar la evaporacin del Mar Mediterrneo

concluy, que era tambin midi la

al flujo de los ros tributarios. Mariote

velocidad del flujo en el ro Sena.

Hidrologa / Davd Cedeo

En hidrologa superficial se desarrollaron muchas frmulas

instrumentos de medicin que permitieron el inicio de la medicin sistemtica

Abbot (1861) efectuaron mediciones

de corriente fue inventado por Price en 1885. Durante este

perodo el gobierno de los Estados Unidos fund una serie de agencias

El intervalo de tiempo comprendido entre 1900 hasta 1930 se denomin

"Perodo de Empiricismo" (Chow, 1964), debido a la gran cantidad de

incorrectas. Las agencias del gobierno incrementaron sus esfuerzos en la

se organizaron cierto nmero de sociedades

el avance de esta ciencia. Por ejemplo, el U.S. Bureau of

Reclamation (1902), el Forest Service (1906) V el U.S. Army Engineers

Waterways Experiment Station (1928) fueron organizadas durante este

Hidrologa / David Cedeo

perodo. Adems, The International Association of Scientific Hidrology

Norton (1933) desarroll la teora

de infiltracin, Theis (1935) formul la

ecuacin de flujo inestable en hidrulica de pozos y

Gumbel (1950) propuso

el uso de la distribucin de valores extremos para el anlisis de frecuencias