Hidrologia Estudiantes Ing Civil
Hidrologia Estudiantes Ing Civil
Hidrologia Estudiantes Ing Civil
HIDROLOGIA
para estudiantes de Ingeniería civil
666
Obra auspiciada por
66 CONCYTEC
-
6 •
•
¡ .. - " . - - , , - - _ . , . - - " , '''''''''--
""''i','t :'l";;;;;¡j't
¡. ';",
,;"Ji'¡ltldi' ,¡:¡III:"''''lf
i : ¡ .
t, I ' 1,"
I
,.,
t
! 1 ' t', { )
'
I .... ',...
... .. -
" '" "
.. _-.. __.. _.. _--_._--j
i
H 1 D ROL O G 1 A
para estudiantes de ingeniería civil
LIMA - PERU
Obra auspiciada por el CONCYTEC
Segunda Impresión
A mi esposa Yolanda
A mis hijas Claudia y Mónica
PRESENTAC ION
Como en otras ramas del conocimiento en .Hidro10gía son muy escasas las
publicaciones peruanas. Aparte de las obras del Ing. M.S. Segundo
ga y del Dr. Medardo Malina, ambos de la Universidad Nacional Agraria,
no se conocen de similar envergadura y esto no obstante la sólida
formación y rica de numerosos profesionales peruanos.
Para llenar en parte este vacío es que publica el presente libro.
Sale a luz gracias a los auspicios del CONCYTEC y con él se intenta or-
denar un poco la de la Hidrología en el país y, por qué no,
motivar en algo los trabajos de investigación en este campo.
El libro desarrolla el curso que con el mismo nombre se imparte como
obligatorio a los alumnos eje ingeniería civil de la Pontificia Universl
dad Católica del Perú. Contiene la 'descripción de los principales ele-
mentos del ciclo hidrológico y. los métodos ordinarios de solución a los
problemas hidrológicos 4ue'se:presentan con más fr.ecuencia inge'niero
civil, No obstante que desde el principio se maneja la Estadística co-
mo herrami enta'de' trabajo en el aná1 i s is y sol uéión de los
problemas, el libro contiene un capitulo completo dedicado a 'la Hidro10
gía Estadística donde son con cierto detenimiento 10sm.QOe10s
hidrológicos probabilísticos yen' fonna somera los modelos estocásticos,
También se ha dedicado un capítulo aparte hacer referencia' a la
hidrología peruana. El sobre el particular .es
el r,ea1ízado alrededor de 1980 p.or un equipo de ingenieros peruanos e
italianos, al amparo del Convenio de Cooperación Técnica suscrito por
el Instituto Ita10-Latino Americano (lILA), el Servicio Nacional de Me-
teorología e Hidrología (SENAMHI) y la Universidad Nacional de Ingenie-
ría (UNI). Resultado de dicho trabajo conjunto es la Rub1icación titu-
lada Estudio de la Hidrología del Perú, que contiene una descripción de
las metodologías empleadas y las conclusiones a manera de informaclón
lista a ser utilizada en la elaboración de proyectos hidráulicos en el
país.
Agradezco en primer término al CONCYTEC.y' a sus autoridades', porJ:lue sin
su apoyo no hubiera sido posible la publicación de este libro; Expreso
igualmente mi agradecimiento aL Ing. Manuel Garda Naranjo por su valio
sa contribución bibliográfica; asimismo a la Srta. Martha Calderón y a
la Srta. E1isabeth Ramos, a ambas por su encomiable trabajo en la prepa
ración del original. -
El autor
CONTENIDO
Página
CAPITULO 2 LA PRECIPITACION (4 h)
2.1 Introductión 15
2.2 Medición de la Precipitación 17
2.3 Análisis de los Datos Pluviométricos 19
2.3.1 Estimación de Datos Faltantes 21
2.3.2 Análisis de Consistencia 26
2.3..3 Extens ión del Regi s tro 27
2.4 Estudio de la Cuenca 28
2.5 Precipitación Media en la Cuenca 30
2.6 Masa de la Precipitación Media en la
Cuenca 33
2.7 Curvas Intensidad-Duración-Frecuencia 35
2.8 Problemas 39
CAPITULO 4 LA INFILTRACION (2 h)
4.1 Descripción 61
4.2 Capacidad de Infiltración 63
4.3 Factores de la Infiltración 65
4.4 Medición de la Infiltración 66
iii
Página
4.5 El Ciclo de Escorrentía 66
4.6 Estimativos de la Escorrentía Usando Infiltración 67
4.7 Problemas 71
CAPITULO 6 EL CAUDAL (5 h)
6.1 La Curva de Descarga 99
6.2 Medición de Caudales 99
6.3 Curva de Descarga de Corrientes sin Aforar 104
6.4 Análisis ,de la Información Hidr,ométrica 106
6.5 La 'Curva de Durad ón 107
6.6 La Curva Masa 109
6.7 Problemas 111
1.1 Generalidades
Los proyectos hidráulicos son de dos tipos: los proyectos que se refie-
ren al uso del agua y los que se refieren a la defensa contra los daños
que ocasiona el agua. Los proyectos típicos de uso del agua son los
de abastecimiento de agua potable, los de irrigación y los de
miento hidroeléctrico; comprenden, además, los de navegación, recrea-
ción y otros. Los proyectos típicos de defensa son los de drenaje
no, drenaje vial y drenaje agrícola; comprenden, además, los de encausa
miento de ríos, los de defensa contra las inundaciones y otros.
En el Perú estamos bastante familiarizados con estos dos tipos de pro-
blemas que se presentan con el agua, los de utilización y los de defen-
sa.
El estudio de nuestros recursos hidrológicos corre por cuenta del Esta
do, siendo su objetivo proporcionar a los ingenieros los elementos para
el aprovechamiento y el control del recurso agua.
1
" I/
-0-
/1"
8
III 1
I I I I 5 5
I
¡ I I
¡II 1
he 6
10
I t1b 1
¡ I I I
III
II !
RIO
MAR
1.3 La Atmósfera
El interés de su estudio en Hidrología radica en que en ella tiene lugar
parte del ciclo hidrológico. Se define como aquella capa de aire que ro
dea a la tierra y donde se realiza parte del ciclo hidrológico.
La atmñsfera resulta comportándose como un gran reservorio de vapor de
agua, un sistema amplio de transporte de agua y un gran colector de ca
1or.
1.4 La Temperatura
La temperatura es un factor importante del ciclo hidrológico inter
viene en todas sus etapas. Desde el punto de vista práctico, la
tura interviene como parámetro en las fórmulas para calcular la evapora
ción y en las fórmulas'para calcular las necesidades de agua de riego de
las plantas. Como prácticamente en todas partes hay registros de
ratura, su empleo está plenamente justificado.
Radiómetros
Los instrumentos que miden la intensidad de energía radiante reciben el
nombre genérico de radiómetros, de los cuales hay varias versiones.
En vista de la importancia que tiene la radiación solar se podría pen-
sar que existe una amplia red de radiómetros en el país, pero esto no
es así. Las razones principales son el elevado costo de equipos y la
exigencia de personal especializado para su servicio.
Heliógrafo
El heliógrafo es un instrumento sencillo que mide el número de horas de
insolación en cada día. Consiste de una esfera maciza de cristal y un
papel sensible que va siendo quemado mientras el sol brilla. El número
de horas de sol es un parámetro que interviene en el cálculo de la eva
poración.
5
La humedad atmoférica interesa a la Hidrología por dos motivos: por ser
el origen de las aguas que caen por precipitación y porque determina en
cierto modo la velocidad con que tiene lugar la evaporación.
Tensión de vapor. En toda mezcla de gases cada gas ejerce una presión
parcial independientemente de los otros gases; la atmósfera es una mez
c1a de gases; la presión parcial que ejerce el vapor de agua se llama
tensión de vapor. Se puede escribir:
ea = p - pi
ea tensión de vapor
p presión del aire húmedo
pi presión del aire seco
Ejemplo 1.1
Encontrar a cuántos milibares equivale 1 mm. de Hg.
pa = 760 mm Hg = 1.033 Kg/cm2
1 mm Hg =
1.033
--no =
1.033
--no x 9.81 x 10 5 dinas
cm2 ----cñ12
0.0133 X 10 5 dinas
Cñi2 o barias
= 1. 33 milibares.
6
TABLA 1.1 TENSION DE VAPOR DE SATURACION ( es ) EN
mm. DE MERCUR 10
t O .1 .2 .3 .4 .5 .6 .7 .8 .9
- -a
10 2.15
9 2.-32 2.30 2.29 2.27 2.26 2.24 2.22 2.21 2.19 2.17
8 2.51 2.49 2.47 2.45 2.43 2.41 2.40 2.38 2.36 2.34
7 2.71 2.69 2.67 2.65 2.63 2.61 2.59 2.57 2.55 2.53
6 2.93 2.91 2.89 2.86 2.84 2.82 2.80 2.77 2.75 2.73
5 3.16 3.14 3.11 3.09 3.06 3.04 3.01 2.99 2.97 2.95
4 3.41,. 3.39 3.37 3.34 3.32 3.29 3.27 3.24 3.22 3.18
3 3.67 3.64 3.62 3.59 3.57 3.54 3.52 3.49 3.46 3.44
2 3.97 3.94 3.91 3.88 3.85 3.82 3.79 3.76 3.73 3.70
1 4.26 4.23 4.20 4.17 4.14 4.11 4.08 4.05 4.03 4.00
O 4.5.8 4.55 4.52 4.49 4.46 4.43 4.40 4.36 4.33 4.29
O 4.58 4.62 4.65 4.69 4.71 4.75 4.78 4.82 4.86 4.89
1 4.92 4.96 5.00 5.03 5.07 5.11 5.14 5.18 5.21 5.25
2 5.29 5.33 5.37 5.40 5.44 5.48 5.53 5.57 5.60 5.64
3 5.68 5.72 5.76 5.80 5.84 5.89 5.93 5.97 6.01 6.06
4 6.10 6.14 6.18 6.23 6.27 6.31 6.36 6.40 6.45 6.49
5 6.54 6.58 6.63 6.68 6.72 6.77 6.82 6.86 6.91 6.96
6 7.01 7.06 7.11 7.16 7.20 7.25 7.31 7.36 7.41 7.46
7 7.51 7.56 7.61 7.67 7.72 7.77 7.82 7.88 7.93 7.98
8 8.04 8.10 8.15 8.21 8.26 8.32 8.37 8.43 8.48 8.54
9 8.61 8.67 8.73 8.78 8.84 8.90 8.96 9.02 9.08 9.14
10 9.20 9.26 9.33 9.39 9.46 9.52 9.58 9.65 9.71 9.77
11 9.84 9.90 9.97 10 .03, 10.10 10.17 10.24 10.31 10.38 10.45
12 ·10.52 10.58 10.66 10.72 10.79 10.86 10.93 11.00 11.08 11.15
13 11. 23 11.30 11.38 11.46 11. 53 11.60 11.68 11. 76 11.83 11. 91
14 11.98 12.06 12.14 12.22 12.30 12.38 12.46 12.54 12.62 12.70
15 12.78 12.86 12.95 13.03 13.11 13.20 13.28 13.37 13.45 13.54
16 13.63 13.71 13.80 13.90 13.99 14.08 14.17 14.26 14.35 14.44
17 14.53 14.62 14.71 14.80 14.90 14.99 15.09 15.17 15.27 15.38
18 15.46 15.56 15.66 15.76 15.86 15.96 16.06 16.16 16.26 16.36
19 16.46 16.57 16.68 16.79 16.90 17.00 17.10 17.21 17.32 17.43
20 17.5.3 17.64 17.75 17.86 17.97 18.08 18.20 18.31 18.43 18.54
21 18.65 18.77 18.88 19.00 19.11 19.23 19.35 19.46 19.58 19.70
22 19.82 19.94 20.06 20.19 20.31 20.43 20.58 20.69 20.80 20.93
23 21.05 21.19 21.32 21.45 21.58 21.71 21.84 21.97 22.10 22.23
24 22.27 22.50 22.63 22.76 22.91 23.05 23.19 23.31 23.45 23.60
25 23.75 23.90 24.03 24.20 24.35 24.49 24.64 24.79 24.94 25.08
26 25.31 25.45 25.60 25.74 25.89 26.03 26.18 26.32 26.46 26.60
27 26.74 26.90 27.05 27.21 27 .37 27.53 27.69 27.85 28.00 28.16
28 28.32 28.49 28.66 28.83 29.00 29.17 29.34 29.51 29.68 29.85
29 30.03 30.20 30.38 30.56 30.74 30.92 31.10 31.28 31.46 31.64
30 11.82 32.00 32.19 32.38 32.57 32.76 32.95 33.14 33.33 33.52
7
Humedad absoluta y humedad relativa. Son dos formas de expresar la hume
dad atmósferica. Se denomina humedad absoluta a la masa de vapor de
gua, medida en gramos, contenida en 1 m3 de aire.
ha ::: 216.7 ea ... (1.1)
T
ha humedad absoluta, en gr/m 3
ea tensión de vapor, en milibares
T temperatura absoluta en grados Ke1vin (OC + 273)
Veamos brevemente cómo se obtiene (1.1).
La densidad del aire seco viene dada por la fórmula:
Ps gr/cm 3
ea milibares
R constante del aire seco = 2.87 x 10 3
T °K
La gravedad específica del vapor de agua con respecto al aire seco es
0.622, 1uego :
ea
Pa = 0.622
expresión en 1a- cual ahora ea está en mm. de mercurio. Para una tem
peratura de unos 15°C, T = 288 °K. Luego:
ha::: ea
El valor máximo de ea es la presión de saturación maXlma, es decir unos
56 mm. de mercurio. Quiere decir que a puede alcanzar a 10 más valores
de ::: 56 gr/m 3 , y por 10 general sus valores son menores.
La humedad relativa es la relación entre la tensión de vapor actual y la
tensión de vapor de saturación ala misma temperatura. Se expresa en
porcentaje:
hr ea x 100 % (1. 2)
=
---e-s
Lo más frecuente es que br sea medida. La medición de hr se realiza
mediante 'instrumentos simples llamados' sicrómetros. Entonces la fórmu
la (1.2) sirve más bien para encontrar el valor de ea'
La relativa es el indice que mejor refleja la sensación de hu
medad que experimentan los seres vivos que se encuentran en una
fera húmeda.
Punto de rocío. Es la temperatura a la cual el vapor de agua de aire
que se considera se hace saturante. Para obtenerlo, se usa simplemen
te la TABLA 1.1, buscando en ella la temperatura para la que es igu:!
la a la ea dada. También es susceptible de ser medida mediante
trumentos sencillos llamados higrómetros de evaporación.
Cuando la temperatura desciende durante la noche, a una temperatura que
corresponde a que el vapor de agua de la atmósfera resulte saturante ,
el vapor de agua se condensa en pequeñas gotitas sobre la superficie
de las hojas formando el rocío.
1.8 El Clima
La palabra "clima" deriva de una voz griega que significa inclinación, a
ludiendo seguramente a la inclinación del eje terrestre. Como se sabe
las estaciones tienen lugar debido al movimiento de traslación de la
Ti erra al rededor del sol, con su eje de rotaci'ón i ncl i nado con respecto
al plano de traslación.
Son numerosas las definiciones que existen de clima, pero todas ellas a-
luden al estado medio de la atmósfera. Para la Organización Meteorológ'i
ca Mundial, clima es el "conjunto fluctuante de condiciones atmosféricas
caracterizado por los estados y la evolución del tiempo, en el curso de
un período suficientemente largo en un dominio espacial determinado".
Los elementos que permiten distinguir un-clima de otro son: la temperat,!!
ra, la precipitación, la presión, el viento y la radiación solar. Los
dos primeros son los principales.'
Los factores que condicionan el clima son: la latitud, la altitud, y la
continentalidad. La latitud determina la intensidad de radiación solár,
la altitud determina la temperatura. La continentalidad se refiere a la
10
mayor o menor proximidad de un lugar a los mares. Muchas veces juegan
papel importante en el condicionamiento del clima las corrientes mari
nas. Otros factores de importancia eventual son la orientación, los
vientos dominantes, la naturaleza del terreno y la vegetación.
Clasificación de climas
El objeto de clasificar los climas radica en poder establecer compara-
ciones. Esto es muy importante en Hidrología, porque hace posible a-
plicar las mismas fórmulas en lugares de clima similar. Para el caso
del Perú, es particularmente frecuente que la región del proyecto ca-
rezca de estaciones y que tenga, por eso, que usarse, registro de
otras regiones de clima similar. En climas similares, la temperatura
y la precipitación son similares en magnitud, variación anual y distri
bución.
Puesto que el clima queda definido por una compleja combinación de
mentos, y además viene definido por una no menos combinación de facto
res, es muy difícil intentar una clasificación única de los
mos de clima que se presentan.
Recordemos que, en principio, se distinguen tres zonas en la superficie
terrestre:
- zona tórrida, comprendida entre el Trópico de Cáncer (23°27 I N) y el
Trópico de Capricornio (23°27 S).
I
13
DESARROLLO HISTORICO DE LA HIDROLOGIA
(tomado de la referencia 2)
sigue en la página 42
14
CAPITULO 2 LA PRECIPITACION
.1 Introducción
La principal fuente de humedad para la precipitación la constituye la
desde la superficie de los océanos. Sin embargo, la cerca-
nía a lv} océanos no conlleva una precipitación proporcional, como lo de
muestran.muchas islas desérticas. Son los del clima ya
dos (latitud, altitud, continentalidad, corrientes marinas, vientos domi
nantes) y las barreras orográficas, las que determinan la humedad atmos
férica sobre una región.
Definición
Se define precipitación a toda forma de humedad, que, originándose en
las nubes, llega hasta la superficie terrestre. De acuerdo a esta defl
nición, las lluvias, las pranizadas, las garGas y las nevadas son formas
distintas del mismo fenómeno de la precipitación. En Estados Unidos, la
lluvia se identifica segGn su intensidad, en:
- ligera, para tasas de caída de hasta 2.5 mm/h
- moderada, desde 2.5 hasta 7.6 mm/h
- fuerte, por encima de 7.6 mm/h
Formación
Debido a su calentamiento cerca de la superficie, motivado por diferen-
cias de radia<;ión, las maSé\S de aire ascienden hasta alturas 'de enfria
llegar a esto conlleva
Supomendo que el alre esta saturado,· o caSl saturado,
para que se forme neblina o gotas de agua o cristales de hielQ, se re
quiere la presencia de nGcleos de condensación (en los dos primeros
sos) o de congelamiento (en el tercero). Los nGcleos de condensación -
consisten de productos de combustión, óxidos de nitrógeno y minGsculas
partículas de sal; los nGcleos de congelamiento consisten de minerales
arcillosos, siendo el caolín el más frecuente.
Después de la nucleación se forman finísimas gotitas de diámetro medio
de aproximadamel'lte 0.02 mm. y como las gotas de lluvia tienen ·un diáme-
tro medio de aproximadamente 2 mm., significa que se produce un aumento
del ,orden de un millón de veces en el volumen de las gotitas. Este
enorme aumento de tamaño se produce por la uni6n entre sí de numerosas
gotitas y esta uni6n se explica por:
- la atracción electrostática entre las gotitas que conforman las nubes;
las microturbulencias dentro de la mas¿ de la nube;
- el barrido de las gotitas más finas por las gotas mayores;
- la diferencia de temperaturas: gotas más frías se engrosan a expen
sas de las más calientes.
15
Mantenimiento de la precipitación
Lo que se acaba de exponer explica la formación de las gotas de lluvia
dentro de la masa de la nube, pero esto no quiere decir que las gotas
así formadas llegarán a .la superficie terrestre, o, en otras palabras
que el volumen de agua contenido en la nube es igual al volumen
agua llovida. Mediciones realizadas demuestran que lo normal es
el agua de lluvia que cae a tierra sea mucho mayor que el agua r-: .• teni
da en 1a nube. La úni ca exp 1i caci ón es que 1as nubes se conti
nuamente durante el proceso mismo de la formación de las precipitacio-
nes, lo que significa una alimentación constante a partir del vapor
de agua de los alrededores; esto se produce principalmente:
- cuando existe una turbulencia dentro de la nube que provoca y facilita
la renovación del vapor de agua;
- cuando hay movimiento del aire húmedo desde las partes bajas, es decir
un movimiento vertical ascendente.
La lluvia artificial
De tiempo en tiempo se habla de la lluvia artificial en el Perú, como
una solución al riego de las zonas áridas de la Costa, sin que hasta
ahora se haya logrado concretar algo. Esto se explica por lo compleja
que resulta en realidad la producción de la lluvia artificial. En los
experimentos que se vienen realizando en otros países se usa para el
bombardeo de las nubes tanto el dióxido de carbono sólido (hielo seco)
como el yoduro de plata; ambos agentes actúan como núcleos de congela
miento. El envío se hace por medio de avionetas, globos, cohetes y ge
neradores. Aun cuando el panorama actual no es del todo claro, hay
el optimismo de lograr a corto plazo la producción a costo razonable de
lluvia artificial
Tipús de precipitaciones
Las precipitaciones se clasifican en t,es el factor
sable del levantamiento del aire que favorece el enfriamiento necesario
para que se produzcan cantidades significativas de precipitación.
a) Precipitaciones convectivas. Son causadas por el ascenso de aire cá
lido más liviano que el aire frío de los alrededores. Las diferen
cias de temperatura pueden ser sobre todo el resultado de calenta-
mientos diferenciales en la superficie o en la capa superior de la
capa de aire. La precipitación convectiva es puntual y su intensi-
dad puede variar entre aquella correspondiente a lloviznas ligeras y
aguaceros. '
b) Precipitaciones orográficas. Resultan del ascenso del aire cálido
hacia una cadena de montanas. Las regiones que quedan del otro la
do de las montañas pueden sufrir la ausencia de lluvias," puesto que
todas las nubes son interceptadas y precipitadas en el lado de don
de ellas provienen.
Es el caso de la Selva Alta de nuestro país, la región más'lluviosa,
donde las nubes provienen de la Selva Baja.
c) Precipitaciones ciclónicas. Se producen cuando hay un encuentro de
nubes de diferentes temperaturas: las más calientes son impulsadas a
las partes más altas donde precipitan.
16
En la naturaleza, los efectos de estos tres tipos de enfriamiento están
inter-relacionados y la precipitación resultante no puede identificarse
como de un solo tipo.
e" b
FitjJ.2.1 PLUVIOMETRO
6
5
4
3
2
18
2.1
B. Febrero 16 empezó a las 9 p.m.
17 terminó a las 9.30 a.m.
empezó a las 11 a.m.
terminó a la 1 p.m.
medida a las 6 p.m. = 5.56 pg.
C. Febrero 16 empezó a las 11 p.m.
17 medida a las 6 p.m. = 2.06 pg.
D. Febrero 16 empezó a las 10 p.m.
17 medida a las 8 a.m. = 3.40"
terminó a la 1.30 p.m.
medida a las 6 p.m. = 4.06"
A = pluviómetro registrador
B, C, D = pluviómetros no registradores con medida
diaria a las 6 p.m.
Pulg. 6 ...-------r---.,.-----r-----r----,.-----r---.
e- 2.0 11
E F M A J J A S O N O
1951 156.3 133.8 55.2 6.4 0.0 2.1 3.4 3.1 33.6 1.4 6.8 56.5 458.6
1952 126.3 163.9 152.2 15.0 0.0 0.0 0.0 0.0 0.0 9.0 76.4 76.4 619.2
1953 100.3 255.9 "86.0 31. 9 11.4 0.0 0.0 0.0 0.'0 11. 5 29.9 27.1 554.0
1954 123.6 154.4 221.6 17.8 18.0 6.0 3.0 0.0 5.0 21.3 22.8 160.9 754.4
1955 114.1 81.1 22.6 1.2 1.7 0.0 0.0 0.4 2.2 0.0 11.2 14.0 248.5
1956 122.4 164.7 150.0 3.5 0.0 27.1 0.0 0.0 8.7 5.8 1.6 97.0 580.8
1957 82.3 145.5 123.2 5.3 6.3 0.0 2.5 2.7 13.8 24.2 31.2 24.9 461.9
1958 64.0 131.9 102.6 25.1 7.9 1.0 0.0 0.0 5.0 5.9 9.7213.5 566.6
1959 193.4 66.1 225.6 10.8 0.0 0.6 0.0 6.1 24.1 29.9 68.0 52.9 677 .5
1960 152.4 152.3 91. 1 55.1 28. O 0.8 0.0 24.8 17.7 6.3 103.5 163.7 795.7
Pro- 123.5 145.0 123.0 17.2 7.3 3 . 8 O. 9 . 3. 7 11. O 11. 5 36. 1 88.7 "
571. 7
me-
dio
S
x
=\
V
I ¿ (x - x)
n-1
2 _ "
- V
L X2-
n-1
X ¿x ••• (2. 1 )
20
Hechos los cálculos se obtiene Sx = 159 mm.
Coeficiente de variabilidad.- Se define: y (x) = x 100 (2.2)
x-
Sx 159 x 100
y (x) =
-x x 100 = = 27.8 %
571. 7
Si suponemos que las precipitaciones anuales en Sibayo constituyen una po
blación normal, y que la muestra tiene también una distribución de
cuencias normal como la población, entonces el punto medio de la campana
de Gauss corresponde al valor medio x y los demás valores respecto a la
media se distribuyen así:
- el 50 % de los datos se encuentra entre x- + 2 S
"3 x
- el 68 % de los datos se encuentra entre -x +
Sx
- el 95 % de los datos se encuentra entre -x +
2 Sx
Se interpreta diciendo que, en la localidad de Sibayo:
- es de esperar una precipitación anual comprendida entre 678.0 y 466.0
mm., con un 50 % de probabilidad;
- es de esperar una precipitación anual comprendida entre 731.0 y 413.0
mm., con un 68 % de probabilidad;
- es de esperar una precipitación anual comprendida entre 890.0 y 254.0
mm., con un 95 % de probabilidad.
Luego del trabajo que se acaba de hacer, ya es posible tomar decisiones
para el diseño de alguna estructura hidráulica en particular que se
yecte en la zona de Sibayo, en Puno.
El manejo estadístico de la información pluviométrica, es decir el estu
dio de su comportamiento según un modelo matemático, sólo es posible
realizarlo cuando la información reúne estos tres requisitos: es comple
ta, consistente y de extensión suficiente. Es por eso que, una
ción pluviométrica antes de ser estudiada en su comportamiento debe ser
revisada en estos tres aspectos, en la forma cómo se describe enseguida.
2.3.1 Estimación de Datos Faltantes
Frecuentemente se halla uno con que faltan datos en los registros
de lluvias. Esto se debe a ausentismo del operador o a fallas
instrumentales.
Se llama correlación a la operaclon o procedimiento por medio del
cual se completan los datos faltantes. Para ello se utilizan los
datos de estaciones índices. que sí tienen los datos completos y
que se seleccionan de modo que estén lo más cerca posible y sean
de altitud parecida a la estación en estudio. Distancia y alti-
tud son pues los factores principales para la selección de las
estaciones índice.
21
Métodos de estimación
Método del U.S. Weather Bureau.- Si los datos faltantes son llu
vias diarias, se escogen tres estaciones índice A, B, C.
a) Si la precipitación anual media en cada estación índice ( xA,
xB, xC) está dentro de un 10 % de la correspondiente a la
estación incompleta ( x ), un promedio aritméticQ simple de
las precipitaciones en las estaciones índice da una estimación
adecuada. En los ejemplos que siguen las precipitaciones es
tán en rTDn.
1o 2. 3
Estación x- fl % día j
A 680 10 1.5 15
B 710 40 6.0 20
C 701 31 4.6 25
x 670
X - 15 + 20 + 25 = 20 rrun.
- 3
px = -1
3
(x
-- pA + - x- pB + - x)
-
pC (2.3)
xA xB xC
22
-x =
731 + 690 + 306 + 610 = 584.3
4
/
'
, ,
••• 0.
r =
L ( x-x )( y-y) (2.5)
(n-1) Sx. Sy
23
r = o significa que no existe ningún grado de asociación entre
los valores de x y los valores de y (correlación nula).
r = 1 ,significa que los puntos del diagrama de dispersión se
alinean en una recta de pendiente positiva (correlación
directa óptima).
r = -1 , significa que los puntos del diagrama de dispersión se
alinean en una recta de pendiente negativa (correlación
inversa óptima).
y y y
...
r =O
x le x
yl =a+b(x-x . . .. (2.8)
Ejemplo 2.5
Completar la información para la estación pluviométrica A a par-
tir de las estaciones índice B y C, próximas a A y de altitudes
parecidas, mediante el método de la recta de regresión. Las ci
fras se refieren a precipitaciones anuales en mm.
24
A B C
1967 786 628 765
8 846 708 876
9 1332 1112 1020
1970 918 816 641
1 830 918
2 930 803 781
3 1115 1020 849
4 887 867 807
5 800 1056 875
6 857 847 947
7 930 756 889
8 918 799
9 888 793 871
1980 915 1002 1000
1 817 831 933
2 999 797 849
y' = y- + r Sx x-x )
138.19
. y' = 930 + 0.59 133.67 x - 860 )
25
2.3.2 Análisis de Consistencia
Cualquier cambio en la ubicación como en la exposición de un plu-
viómetro puede conllevar un cambio relativo en la cantidad de llu
via captada por el pluviómetro. El registro completo publicado
representará condiciones inexistentes. Un registro de este tipo
se dice que es inconsistente.
Una forma de detectar las inconsistencias es mediante las curvas
doble másicas.
Una curva doble más i ca' se construye 11 evando en ordenadas los va
lores acumulados de la estación en estudio y en abscisas los valo
res acumulados de un patrón, que consiste en el promedio de va-
rias estaciones índice.
6,000
Itao!
V
79
78/
·76
.7&
-'
74
·75
7
y
v<
1,000
l,A(,
°o 1,000 2,000 3,000 4,000 5,000 6,000
PATRON (mm)
27
20 Relleno de datos faltantes. Se emplea el método de la recta
de regresión.
30 Extensión del registro. Con cualquiera de las dos formas in
cadas.
\
I
J
28
Frecuencia en %
6 10 10 20 30
m.s.n.m
4,000
olígono de fr ecuenc io
da oltitll des
5,000 "'" J
1
............. _ curv P ,
2,000 '< hips pmetri a
1,000
25
J
50 75
, '"1\ 100
Area en k mE
125
29
2.5 Precipitación Media en la Cuenca
A partir de las lluvias medidas ep los pluviómetros es posible calcular
la precipitación media en la cuenca. Singularmente útil resulta la pre
cipitación media anual, o módulo anual, en la cuenca. -
Los pluviómetros deben ubicarse estratégicamente y en número suficiente
para que la información resulte de buena calidad.
El problema entonces se refiere al cálculo de la lámina o altura de agua
que cae en promedio durante 1 año en una cuenca. Existen para ello va-
rios métodos disponibles, de los cuales los más usados son los tres que
se describen a continuación.
Promedio aritmético.- Si pI, p2, ....... , pn son las precipitaciones
anuales observadas en diferentes puntos de la cuenca, entonces la preci-
pitación anual media en la cuenca es:
= p1+p2+.: ..... +pn (212)
P n •. , .
....
", ,
" \
\
I ---/ \
I
I
I
....
..... _- -" /
/
30
Isoyetas.- Se define isoyeta la línea de igual precipitación.
El método consiste en (fig. 2.8):
1° Tra.zar "as isoyetas, interpolando entre las diversas estaciones,
de modo similar a se trazan las curvas de nivel;
2° Hallar las áreas al, a2, ..... , an entre cada 2 isoyetas seguidas;
30 S·1 po,
'\ P1 , ..... , pn son 1 ···
as preclpltaclones anua 1es representa das
por 1as i soyetas respecti vas, entonces':
po +2 pI al + ..... + pn-l 2+ pn a
n
p = ... (2.l4)
al + ....... + an
es la precipitación anual media en la cuenca.
ISOVETAS /200
I
4-- CID
o
o
ID
o
o
31
2° Se hall a la precipitación sobre cada polígono operando- con las 1S0
yetas.
a
h = L hm. al (2.16)
pi =
precip. sobre el polígono x relación de áreas (2.17)
precip. en la estación
Ejemplo 2.6
EST.D
o
4.0.-
32
TABLA 2.2 CALCULO DE LOS PESOS
Para que el método modificado resulte práctico, las isoyetas y los poli
ganas se dibujan una sola vez, para la información de mayor confianza,
a fin de obtener los pesos de las estaciones. De ahí en adelante, para
evaluar la precipitación media en la cuenca sólo se requieren los datos
de precipitación en las estaciones (fórmula 2.15).
El método de los polígonos Thiessen recién descrito, junto el ejem-
plo numérico, están contenidos en la referencia 6, de donde han sido to
mados con ciertos ajustes para su adaptación al presente texto.
33
TABLA 2.3 CURVA MASA DE LA PRECIPITACION MEDIA
(1) (2) (3) (1) (2) (3) (1) (2 ) (3). (1) (2) (3)
O O
1 0.17 0.17 .053 O .053
2 O 0.33 0.16 .050 O 0.15 0.15 .006 .056
3 0.20 0.20 .070 0.52 0.19 .059 0.09 0.09 .026 0.29 0.14. .006 .161
4 0.40 0.20 .070 0.80 0.28 .087 0.17 0.08 .023 0.52 0.23 .009 .189
5 0.73- 0.33 .116 1.20 0.40 .124 0.32 0.15 .. 044 0.84 0.32 .013 .297
6 1.20 0.47 .164 1.41 0.21 .065 .52 0.20 .058 1.01 0.17 .007 .294
7 1.20 O O 1.85 0.44 .136 .52 O O 1. 34 0.33 .013 .149
8 2.05 0.85 .298 2.91 1.06 .329 .89 0.37 .107 2.05 0.71 .028 .762
9 2.80 0.75 .262 3.49 0.58 .180 1.22 0.33 .096 2.47 0.42 .017 .555
10 3.15 0.35 .122 4.19 0.70 .217 '1.37 0.15 .044 3.00 0.53 .021 .404
'11 3.90 0.75 .262 4.79 0.60 .186 1. 70 0.33 .096 3.40 0.40 .016 .560
12 4.20 0.30 .105 5.08 0.29 .090 1.83 0.13 .038 3.63 0.23 .009 .242
13 4.40 0.20 .070 5.18 0.10 .031 1.92 0.09 .026 3.73 0.10 .004 .131
14 4.40 O O 5.18 O O 1. 92 O O 3.83 0.10 .004 .004
15 4.59 0.19 .066 5.49 0.31 .096 2.00 0.08 .023 3.97 0.14 .006 .191
16 4.70 0.11 .038 5.56 0.07 .022 2.04 0.04 .012 4.04 0.07 .003 .075
17 4.73 0.03 .010 5.56 O O 2.06 0.02 .006 4.06 0.02 .001 .017
Los valores de la última columna son los promedios pesados de las preci-
pitaciones horarias en la cuenca, valores con los que se puede dibujar
la curva masa.
34
Pulg.
".14
4.0 Flg. 2.10
CURVA MASA DEL
EJEMPLO 2.7
3.0
2.0
1.0
f) Z 3 . .. . 5 6 7 8 9 10 11 12 13 14 15 16 17 horas
35
Si bien este asunto del análisis de tormentas ha podido posponerse para
ser estudiado en el Capítulo 10, 10 vamos a tratar aquí porque necesita-
mos el resultado del análisis de tormentas para una buena interpretación
de la fórmula racional (Capítulo 7).
El análisis de tormentas se hace a través de siete etapas o pasos.
Paso 1.- Se parte de un p1uviograma, es decir el registro de un p1uvió-
grafo, como el de la fig. 2.2.
Paso 2.- Se hace la siguiente tabulación., a partir del p1uviograma.
11.00
60 0.5 0.5
12.00
50 8.5 10.2
12.50
70 10.0 8.6
14.00
140 4.5 1.9
16.20
36
10
1
• HISTOGRAMA
8
E
.J
16 17
11 lE 13 14 15
t horas
37
Año Período de 'duraci ón (mi n. )
10 30 60 120 240
1950 102 81 64 42 !@'Z\
1951 83, 70 56 33 16
1952 76 61 42 29 19
1953 .W;2IOS M23
1954 61 58 36 28 14
Paso 6.- Procediendo por separado para cada período de duración, se colo
can en orden decreciente, prescindiendo del año, los valores de la tabla
última.
n=30
100
80
60
40
20
2.8 Problemas
Prob 1ema 2. 1
En una cierta cuenca se han instalado pluviómetros en 4 estaciones A, B,
C, D. Las altitudes de las 4 estaciones son parecidas.
La estación A está situada entre las estaciones B, e, D, a las distan -
cias:
A B = 10 Km.
A e = 5 Km.
,A D = 2 Km.
Problema 2.2
En una cierta cuenca se han instalado 4 pluviómetros totalizadores de
lectura mensual. En un cierto mes del año falta una de las lecturas
mientras que las restantes son 37, 43 y 51. Si las precipitaciones me
dias anuales de estos 3 pluviómetros son 726, 752 Y 840 mm., respectiva
mente, y del pluviómetro incompleto 694 mm., estimar la lectura faltan
te de precipitación mensual.
Problema 2.3
La fig. 2.12 representa el registro de un pluviógrafo durante una cierta
tormenta.
Calcular las intensidades de lluvia durante períodos sucesivos de 1 hora
y dibujar el histograma.
39
cm. !5
4 I / I I
3
V / I I
2 V V I I
1/ / I I 1/
o
7 8 9 10 1/ 12
V
13 14 15 16
11
17
V
18
V
19 20 21
Problema 2.4
Las dos figuras de abajo representan los histogramas de dos tormentas.
lmm/h lmm/h
20 20
10 O
O 2 3 t horas O 2 :3 t horas
Problema 2.5
La tabla 2.4 presenta las precipitaciones anuales de la estación X y las
precipitaciones anuales medias de ,una estación patrón.
a) Examinar la consistencia de la información de X.
b) ¿Cuándo ocurrió un cambio en el régimen? Discutir las posibles
sas.
c) Ajustar la información y determinar la diferencia-en la precipitación
anual media de los 36 años en la estación X.
40
TABLA 2.4 DATOS DEL PROBLEMA 2.5
Problema 2.6
En una cuenca se han instalado 4 pluviómetros. En la figura 2.13 se
presentan las precipitaciones medias anuales y las curvas isoyetas, con
sus correspondientes porcentajes de área. Determinar la precipitación
anual media por medio de los polígonos Thiessen y las curvas isoyetas.
''\-
I --""-
2Z.5 % / -
/
/
o I
/
42.2 % El 0°
450 mm /
I
/
I !560 111'" /
/ e /
/
/
/
/
57!5mm
I
sigue en la págini 72
42
CAPITULO 3 EVAPORACION y , '
3.1 Introducción
La evaporación es una etapa permanente del ciclo hidrológico. Hay
ración en todo momento y desde toda superficie húmeda. Considerada como
un fenómeno puramente físico, la evaporación es el pasaje del agua al es
tado de vapor; sin embargo hay otra evaporación, la provocada por la ac
tividad de las plantas y que recibe el nombre de transpiración.
De modo general, la evaporación se puede estudiar por separado, a partir
ae las superficies libres del agua (lagos, embalses, ríos, charcas), a
partir de la nieve, a partir del suelo y a partir de las plantas (trans-
piración). O bien se puede estudiar la evaporación total en una cuenca,
sin tomar en cuenta las formas particulares que adopta; a esta evapora -
ción total se llama evapotranspiración.
Nosotros estudiaremos preferentemente la evaporación en embalses y la
evapotranspiración. La primera, porque el ingeniero tiene interés en
evaluar la cantidad de agua almacenada que se va a perder por evapora-
ción. La segunda, por sus aplicaciones en los proyectos de irrigación.
El fenómeno de la evaporación a partir de los espejos de agua es comple
jo, pero podemos esquematizarlo del modo que sigue. Las moléculas de la
superficie libre adquieren energía cinética por acción de la energía so
lar y vencen la retención de la masa de agua, salen al aire y se acumu-
lan formando una capa encima del agua; para que continúe el proceso es
necesario remover esta capa de vapor de agua y esto lo hace el viento.
El papel de la temperatura es doble: aumenta la energía cinética de ras
moléculas y disminuye la tensión superficial que trata de retenerlas.'
Factores
De todos los factores que intervienen en la evaporación, los principales
son los meteorológicos: radiación solar, temperatura del aire, la pre-
sión de vapor, el viento y en menor grado la presión atmosférica. Estos
factores son los que provocan la evaporación. Debido a que la radiación
solar es el factor más importante, la evaporación varía con la latitud,
época del año, hora del día y condiciones de nubosidad.
La tasa de evaporación desde un suelo saturado es aproximadamente igual
a la evaporación desde una superficie de agua cercana, a la misma tempe-
ratura. Al comenzar a secarse el suelo la evaporación disminuye, y fi-
nalmente cesa porque no existe un mecanismo que transporte el agua desde
una profundidad apreciable.
En cuanto a los efectos de la calidad del agua, puede decirse que la pre
sencia de sales hace disminuir ligeramente la evaporación. En el agua
de mar, por ejemplo, es del orden de 2 % menor que en el agua dulce.
Quiere decir que los efectos de la salinidad pueden despreciarse en la
estimaci6n de la evaporación de un embalse.
43
3.2 Evaporación en Embalses
La medida directa de la evaporación en el campo no es posible, en el sen
tido en que se puede medir la profundidad de un río, la precipitación ,
etc. Debido a esto se han desarrollado una serie de técnicas para esti
mar la evaporación desde la superficie de un embalse.
.5
2
, / 500 20
- /
5 10
t- 8
0.2 "'O 0.0 At.,
15
'-..... 400 15
"j fa
.4
4
L.90
l
0.3
. 3
r -2 1.0
10 ¡:- .2
.5 300
(J1 0.4
·c
0.5 /""
30 20
-.
'1 t·
N
-j(- /"" E
-4
[ 5 .2 .&
1 .111
,.;
.4 200
j i
I l
"j
,J
u -.Jo .02
-5
,j el [ to O
1.0
IJINI'LO Con el nomogromo le encuentran
-&
- .2
'" j j I 0.5 I
1 f- .01
"".99
0.7 . lo. valore. de El' E 2 ' E s
.05
20· . E, : -1.0 .... /day
"D 0.4 E, = 2.5
-s
Ra
u,
550
5
E,
Lu.,o' E.
= 1.8
= 3.1 ",.. /doy I .00J t.1 100 1 I 0.0
RA ... valor de Angot. Es la cantidad de radiación solar, en
calorías por día en un plano horizontal del cm 2., en-
trante en los límites exteriorres de la atmósfera. Es
una función de la posición geográfica y la época del año
(tabla 3.1).
Latitud E F M A M J J A S O N O
Sur
0° 885 915 925 900 850 820 830 870 905 910 890 875
10° 965 960 915 840 755 710 730 795 875 935 955 960
20° 1020 975 885 765 650 590 615 705 820 930 1000 1025
30° 1050 965 830 665 525 460 480 595 750 900 1020 1065
40° 1055 925 740 545 390 315 345 465 650 840 995 1080
50° 1035 865 640 415 250 180 205 325 525 760 975 1075
Ejemplo 3.1
Averiguar el valor de Eo para los siguientes datos:
t = 20 oC
h= 0.7
U2 = 5 m/sg
n 0.4
0=
RA = 550 cal
cm 2 - día
El se lee en la primera parte del nomograma = -1.0 mm/día
E2 se lee en la segunda parte del nomograma = +2.3 mm/día
E3 se lee en la tercera parte del nomograma = +1.8 mm/día
Luego, Eo = El + E2 + E3
= -1.0 + 2.3 + 1.8
= 3.1 rnm/dfa
46
En términos de calor, se expresa Elo :; 60 Eo
Elo .... calor requerido en cal
cm 2 - día
Eo evaporación en mm/día
R .... energra en ca 1
cm 2 - día
C5 •••• constante:; 117.4 x 10- 9 cal
cm 2 - día
T .... temperatura absoluta :; 273 0 + t oC
Rí :; Rc (1 - r) (3.3)
RB :; (0.47 - 0.077
C5 T4 . re¡) (0.20 + 0.80
\j -- On ... (3.4)
a T4 radiación de Stefan-Boltzmann
ea presión de vapor actual en-el aire, en mm de Hg
-1-7
La cantidad neta de energia remanente en la superficie, y disponi
ble para varias pérdidas, es el llamado calor almacenado H: -
H = RI - RB · . .. (3.5)
48
3° Si la temperatura de la masa de agua permanece constante, o
el lago es poco profundo, o se consideran períodos cortos de
10 a 20 días, A S puede despreci arse.
4° El valor de A es negativo cuando un tanque aislado lleno con
agua, en un desierto caliente y seco, en adición al calor di
recto en su superficie recibe también calor en los lados
lar de advección). Se toma comO'cero cuando el embalse es
grande. Estos efectos de borde se pueden pues despreciar.
Reemplazando: H = E'o + K . . . .. (3.9)
t ti t
49
En términos de calor: Ela = 60 Ea
Ela = 21 (es - ea) (0.5 + 0.54 u2 ) ..... (3.11)
En 1a 3.9:
H = Elo + K = Elo + Elo y ti t
e's - ea
Elo = H
1 + y ti - t
e's - ea
Según la 3.10:
ti - t = els - es
A
Reemplazando:
Elo H
=
H
1 + y e s - es
I 1 1 + _:1- els - es
A els - ea A els - ea
Elo H
1 + y (els-ea) - (es-ea)
A e's - ea
Reemplazando:
H
EI o = -------
1+ :(l - ro)
Ela
50
De donde, finalmente:
Elo =
AH + y Ela ..... (3.12)
II + y
Ejemplo 3.2
Averiguar el valor de Eo, por el método del balance energético, para los
mismos datos del ejemplo 3.1
t = 20°C
h = 0.7
u2 = 5 m/seg.
n
O = 0.4
RA = 550 - -cal
-----
cm 2 - día
t = 20 oC - es = 17.53 mm Hg
ea = hes = 0.7 x 17.53 = 12.27 mm Hg
es ea 17.53 - 12.27 = 5.26 mm Hg
T = 20 + 273 293 °K
RA ca1/cm 2
= 550 día (da to)
RC ca1/cm 2
= 216 día (con 3.2)
RI ca1/cm 2
= 203 día (con 3.3)
RB = 91ca1/cm 2 día (con 3.4)
H = 112 ca1/cm 2 día (con 3.5)
Ela = 354 ca1/cm 2 - día (con 3.11 )
Reemplazando en 3.12:
Elo =
1.05 x 112 + 0.49 x 354 189 ca1/cm 2 - día
1.05 + 0.49
Elo
ElO = = 3.15 mm/día
60
[)1
Ejemplo 3.3
Hallar, adicionalmente al ejercicio 3.2, el valor de ti.
La 3.8: K = y c l (tl-t) .f (u)
K = H - Elo de la 3.9
y = 0.49 constante sicrométrica
c = 21
l
de 1a 3.11
f (u) = 0.5 + 0.54 U2 de 1a 3.11
t = 20 oC dato
Resolviendo:
112 189 = 0.49. x 21 t l
-20) (0.5 + 0.54 x 5)
ti = 17.7 oC
Observamos: ti < t , es decir la temperatura en la superficie libre es
menor que en el aire. Este es el conocido efecto de enfriamiento de la
evaporación.
52
2. Fórmula de Meyer (Inglaterra)
53
Las tasas de evaporación así observadas pueden ser consideradas como maXl
mas y dan una buena aproximación del poder evaporante de la atmósfera.
Aplicando a esas tasas máximas diversos coeficientes de reducción, el in
geniero deducirá los valores más probables de las tasas de evaporación
que le interesan (embalses, terrenos desnudos, terrenos cubiertos de vege
tación, etc). -
Tanques de evaporación.-
Se pueden clasificar en tres grupos, según que estén dispuestos en la su-
perficie del suelo, enterrados en éste o flotando.
Tanques colocados en la superficie del suelo.-
Tienen la ventaja de una instalación sencilla y sus resultados no corren
el riesgo de ser falseados por salpicaduras de gotas de lluvia. Son, en
cambio, muy sensibles a las variaciones de la temperatura del aire y a
los efectos de la insolación. .
A este grupo pertenece el tanque llamado Clase A, del U.S. Weather Bureau.
Tiene un diámetro de 121.9 cm. y una profundidad de 25.4 cm. Está cons
truldo de hierro galvanizado no pintado y colocado sobre un bastidor de
madera a unos 15 cm. del suelo. Es el usado entre nosotros. Para ha-
llar la evaporación en el embalse puede emplearse un coeficiente anual
de 0.7. Para perfodos menores los coeficientes son variables.
Tanques enterrados.-
Son menos sensibles a las variaciones de la temperatura del aire y a los
efectos de la insolación, pero en cambio las gotas de lluvia que caen en
su rededor pueden salpicar, y falsear las medidas. Otra dificultad es
que no se podría descubrir a tiempo una pequeña fuga.
A este grupo pertenece el tanque llamado Colorado, grandemente extendido
en el oeste de los Estaqos Unidos. Tiene la forma de un prisma cuya ba
se es un cuadrado de lado 0.914 m. y cuya altura es de 0.462 m. Es ente
rrado en el suelo de manera que sus aristas superiores quedan 0.10 m. so
bre la superficie del suelo.
Tanques flotantes.-
Su instalación y operación pueden resultar algo complicadas además de
costosas. La tendencia es preferir el tanque Colorado instalado en la
orilla del embalse.
Evaporímetro Piche.-
De amplio uso en estaciones evaporimétricas, está constituldo por un tu-
bo cilfndrico de vidrio de 25 cm. de largo y 1.5 cm. de diámetro. El tu
b0 está graduado y encerrado en su parte superior, mientras que su cober
tura inferior está obturada por una hoja circular de papel filtro
lizado de 30 mm. de diámetro y 0.5 mm. de espesor. fijada por capilari -
dad y mantenida por un resorte.
Llenado el aparato de agua destilada, ésta se evapora progresivamente a
través de la hoja de papel filtro, la disminución del nivel del agua en
el tubo permite calcular la tasa de evaporación (en mm. por 24 horas ,
por ejemplo). El aparato se instala bajo cubierta para mantenerlo ale-
jadO de la lluvia.
54
Reducción de la evaporación.-
A fin de tener una idea del orden de magnitud de la evaporación en embal
ses, vamos a considerar los valores que alcanzó la evaporación en Huanca
yo en un par de años (en mm).
E F M A M J J A S O N D TOTAL
1956 159 1.33 136 149 134 151 176 198 186 229 241 249 2,141
1957 197 142 154 143 137 123 207 208 200 189 217 226 2,143
3.4 Evapotranspiración
Del agua que una planta absorbe del suelo sólo una parte muy pequeña se
queda para formar los tejidos de la planta; el resto regresa a la atmós
fera en forma de vapor, configurando la transpiración.' Este fenómeno de
la transpiración constituye una fase muy importante del ciclo hidrológi-
co, porque es el mecanismo principal por medio del cual el agua precipi
tada a tierra regresa a la atmósfera.
Al estudiar el balance hídrico de una cuenca, el interés recae en la de
terminación de las pérdidas totales de agua, es decir por evaporación y
por transpiración. Además, desde el punto de vista práctico es muy difí
cil evaluar por separado cada pérdida. Las pérdidas totales de agua
constituyen la evapotranspiración.
El término "evapotranspiración potencial" fue introducido por Thornth-
waite y se define como la pérdida total de agua que ocurrlrla si en nin
gún momento existiera deficiencia de agua en el suelo para el uso de la
vegetaci ón.
Se define "uso consuntivo" la suma de la evapotranspiración y el agua u-
tilizada directamente para construir los tejidos de las plantas. La dis
55
tinción entre los términos evapotranspiración potencial y uso consuntivo
es más que todo académica, porque las diferencias numéricas caen siempre
dentro de los errores de medición y por 10 común se tratan como términos
sinónimos.
En los proyectos de irrigación interesa hacer un cálculo previo de las
necesidades de agua de los cultivos. Estas necesidades de agua, que van
a ser satisfechas mediante el riego, vienen a constituir la evapotranspi
ración o el uso consuntivo. Para el cálculo de estas cantidades de agua
se han desarrollado métodos basados en datos meteorológicos, de los cua
les los más conocidos son el de Thornthwaite y el de B1aney-Cridd1e. Am
bos se usan entre nosotros.
56
TABLA 3.2 FACTORES DE CORRECCION DE e
Latitud E F M A M J J A S O N D
Sur
5 1.04 0.95 1.04 1.00 1.02 0.99 1.02 1.03 1.00 1.05 1.03 1.06
10 1.08 0.97 1.05 ·0.99 1.01 0.96 1.00 1.01 1.00 1.06 1.05 1.10
15 1.12 0.98 1.05 0.98 0.98 0.94 0.97 1.00 1.0pl 1.07 1.07 1.12
20 1.14 1.00 1.05 0.97 0.96 0.91 0.95 0.99 1.00 1.08 1.09 1.15
25 1.17 1.01 1.05 0.96 0.94 0.88 0.93 0.98 1.00 1.10 1.11 1.18
30 1.20 1.03 1.06 0.95 0.92 0.85 0.90 0.96 1.00 1.12 1.14 1. 21
35 1. 23 1.04 1.06 0.94 0.89 0.82 0.87 0.94 1.00 1.13 1.17 1. 25
40 1.27 1.06 1.07 0.93 0.86 0.78 0.84 0.92 1.00 1.15 1.?O 1.29
45 1. 31 1.10 1.07 0.91 0.81 0.71 0.78 0.90 0.99 1.17 1. 26 1. 36
50 1. 37 1.12 1.08 0.89 0.77 0.67 0.74 0.88 0.99 1.19 1. 29 1.41
u= L U = K. L f (3.18)
57
TABLA 3.3 VALORES DE P
Latitud E F M A M J J A S O N D
Sur
5 8.68 7.76 8.51 8.15 8.34 8.05 8.33 8.38 8.19 8.56 8.17 8.68
10 8.86 7.87 8.53 8.09 8.18 7.86 8.148.27 8.17 8.62 8.53 8 88
15 9.05 7.98 8.55 8.02 8.02 7.65 7.95 8.15 8.15 8.68 8.70 9.10
.
20 9.24 8.09 8.57 7.94 7.85 7.43 7.76 8.03 8.13 8.76 8.87 9.33
25 9.46 8.21 8.60 7.84 7.66 7.20 7.54 7.90 8.11 8.86 9.04 9.58
30 9.70 8.33 8.62 7.73 7.45 6.96 7.31 7.76 8.07 8.97 9.24 9.85
34 9.92 8.45 8.64 7.64 7.27 ·6.74 7.10 7.63 8.05 9.06 9.42 10.08
38 10.15 8.57 8.66 7.54 7.08 6.50 6.87 7.49 8.03 9.16 9.61 10.34
42 10.40 8.70 8.68 7.44 6.85 6.23 6.64 7.33 8.01 9.26 9.82 10.64
--
46 10.69 8.86 8.70 7.32 6.61 5.92 6.37 7.16 7.96 9.37 10.07 10.97
Ejem210 3.4
Mostrar, en forma tabulada, la manera cómo se efectúan los cálcu-
los en el método de Blaney-Criddle.
58
TABLA 3.5 DESARROLLO DEL EJEMPLO 3.4
Culti vo Alfalfa
Lugar Valle de Salinas, California
Mes t p f k(*) u
oC % mm.
Abri 1 14.39 8.85 130.0 0.6 78.0
Mayo 16.95 9.82· 155.7 0.7 109.0
Junio 18.73 ·9.84 164.2 0.8 131.4
Jul i o 20.23 10.00 170.3 0.85 111. 7
Agosto 19.89 9.41 162.0 0.85 137.7
Setiembre 19.23 8.36 141. 3 0.85 120.1
Octubre 16.78 7.84 123.8 0.70 86.6
3.5 Problemas
Prob 1ema 3. 1
Hallar la evapotranspiración potencial, utilizando el nomograma de Pe n-
man, en el siguiente caso. Campo cultivado en la latitud 40 0 S, en Se-
tiembre, temperatura media del aire 20°C, humedad relativa media 70 % ,
insolación relativa 40 %, velocidad media del viento V2 = 2.5 m/sg., va
lor de la relación evapotranspiración potencial a evaporaciór potenciaT
70 %.
Prob 1ema 3. 2
En una cuenca de tamaño medio, las temperaturas medias mensuales en No-
viembre y Diciembre de 1974 fueron 16.1 y 17.9 oC, respectivamente. Da
do que el índice térmico anual fue 66.9 y las duraciones astronómicas
medias mensuales de esos días fueron 15.00 y 16.20 horas/día, respecti-
vamente, hallar la evapotranspiraciónpotencial para cada mes.
59
Problema 3.4
Dos cuencas vecinas, A y B, tienen similar altitud, clima y uso de la tie
rra. No hay agua subterránea ni hacia ni desde las cuencas. -
La cuenca B es, sin embargo, más grande y tiene una precipitación anual -
mayor. Se dispone de la siguiente información:
A B
Area en Km2 200 3,00
Precipitación anual media (mm) 1,200 1,500
Descarga media (m /sg,)
3 <, 4 ?
[iO
CAPITULO 4 LA INFILTRACION
4.1 Desc..ripción
Cuando llueve, parte de la lluvia del comienzo es retenida en la cobertu
ra vegetal como intercepción y en las depresiones del terreno como alma
cenamiento superficial. Conforme continúa la lluvia, el suelo se cubre
de una delgada capa de agua conocida como detención superficial y el
flujo comienza pendiente <abajo hacia los cursos, lo que la es
correntía superficial. Inmediatamente debajo de la superficie tiene lu
gar la escorrentía subsuperficial y las dos escorrentías, la superficial
y la subsuperficjal, constituyen la escorrentía directa. La infiltra-
ción es el paso del agua a través de la superficie del suelo hacia el in
terior de la tierra; la percolación es el movimiento del agua dentro del
suelo y ambos fenómenos, la infiltración y la percolación, están íntima-
mente ligados puesto que la primera no puede continuar sino cuando tiene
lugar la segunda. El agua que se infiltra en exceso de la escorrentía -
subsuperficial puede llegar a formar parte del agua subterránea, la que
eventualmente puede llegar a los cursos de agua.
El agua de un río, en general. puede así estar formada de dos partes.
Una parte de escorrentía (superficial y que recibe 'el
nombre de escorrentía directa y otra parte de agua subterránea que reci
be el nombre de' flujo base.
s
F
\
\
\ p'
/' '"
61
P = 1 + S + E + F + Pn . . . . .. (4. 1 )
P lluvia total
1 intercepción
S almacenamiento superficial
E evaporación desde el suelo
F infiltración
Pn escorrentía directa, también llamada lluvia neta.
Topografía
Naturaleza de Ondulada Inclinada
la superficie S de 5 a 10% S de 10 a 30%
Cultivos generales 0.60 0.72
Cultivos de pastos 0.36 0.42
Cultivos de bosques 0.18 0.21
Areas desnudas 0.80 0.90
62
TABLA 4.2 VALORES DE e (REFERENCIA 1)
mm/h
fo
- - - - - - -=-=-=-,......--
63
Horton encontró que las curvas de infiltración se aproximan a la forma:
fp = fc + (fo - fc) e -kt (4.4)
Escorrent¡'a
directo
Infiltración
64
Esta segunda manera de obtener la escorrentía directa por separación en
el histograma es también sencilla sólo en apariencia. En primer lugar
se requiere una estaci6ncon p1uviógrafo, en segundo lugar el suelo de
la cuenca no es homogéneo de modo que la curva de infiltración no es la
dibujada y en tercer lugar la determinación de la retención tiene tam-
bién sus complicaciones.
1;5
4.4 Medición de la Infiltración
La mayoría de datos sobre tasas de infiltración se obtienen mediante ensa
yos con un infiltrómetro.
En los casos en que existe una amplia variaClon en los suelos o en la ve-
getación dentro del área, la cuenca se divide en subcuencas homogéneas ,
cada una de las cuales está cubierta por un complejo de suelo único. Re
pitiendo varias veces el ensayo se pueden obtener datos fidedignos de ca
da subárea.
Un infiltrómetro es un tubo u otro contorno metálico diseñado para aislar
una sección del suelo. El área efectiva varía desde menos de 1 pie 2 has
ta varias decenas de pie 2 . El procedimiento clásico consiste en inundar
el infiltrómetro, aplicando agua en una lámina de altura constante sobre
el área encerrada y medir el tiempo que tarda en infiltrarse. Esta prác-
tica está siendo reemplazada por los simuladores de lluvia. Como en este
caso no se puede medir directamente la cantidad de agua que penetra el
suelo, la infiltración se calcula como la diferencia entre el agua aplica
da y la escorrentía directa medida. -
66
sobre el cauce
intercepción '
horas
Indices de Infi1tración.-
Las dificultades inherentes al método directo de infiltración, para eva
1uar la escorrentía directa, han conducido al uso de los índices de in
67
filtración. El más simple de estos es el índice como la
tasa de precipitación por encima de la cual el volumen de lluvia es
igual al volumen de escorrentía directa. Ver fig. 4.6.
. :.
horas
Fig.4.6 EL INDICE
. ,- .
1. e. d.
------ . '. .. ..,. . '"
- - '- - -' --=- =- - '1'
68
Ejemplo 4.1
Se trata de determinar el valor del índice para la cuenca de la
ra 2.9 del ejemplo 2.7.
La Columna de la tabla 2.3 del ejemplo 2.7 contiene los valores
de la precipitación horaria sobre la cuenca, con los que se construye
la curva masa en la cuenca. Si además, una estación de aforo ubicada en
el punto de salida de la cuenca ha permitido medir el escurrimiento di-
recto correspondiente a la misma tormenta, encontrándose que es 2.0 pg.,
es posible obtener el índice por tanteos.
La tabla 4.3 muestra el procedimiento seguido en este caso. El valor de
0.17 pg/h supuesto en el tercer tanteo produce la mejor concordancia en
tre la lluvia neta calculada (1.96 pg.) y la escorrentía directa medida
(2.0 Ver fig. 4.9.
Con este procedimiento se pueden obtener los valores del índice para
varias tormentas importantes, siendo el promedio de estos valores el va-
lor medio del índice para la cuenca.
La duración de lluvia neta es el tiempo durante el cual los incrementos
de precipitación exceden al valor del índice Después del tercer
tanteo se encuentra que esta duración es de 9 horas. En la figura 4.10
se ha dibujado el histograma de lluvia neta.
pg./h
."
",,'
. ..
0.4 • ' I •
').2
o 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
horas
69
TABLA 4".3 CALCULO DEL INDICE
O
1 0.05 0.25 0.15 0.17
2 0.06
3 0.16 " 0.01
4 0.19 0.04 0.02
5 0.30 0.05 0.15 0.13
6 0.29 0.04 0.14 0.12
7 0.15 O O O
8 0.76 0.51 0.61 0.59
9 0.56 0.31 0.41 0.39
10 0.40 0.15 0.25 0.23
11 0.31 0.41 0.39
12 0.24 0.09 0.07
13 0.13 O O
14 O O O
15 0.19 0.04 0.02
16 0.08
17 0.02
Total 4.14 1. 37
--
2.15 1.96
pg /h.
0.6-
0.4-
0.2 -
horas
9 horas
-1
70
4.7 Problemas
Problema 4.1
Determinar la ecuación de la curva capaci dad de infiltración para los
siguientes datos observados.
t mino O 5 10 15 20 25 30 35 40 45
fp cm/h 16.0 11.0 7.9 5.7 4.1 2.8 1.9 1.3 1.0 1.0
Prob 1ema 4. 2
La tabla de abajo indica las láminas horarias de tres tormentas que pro
dujeron escorrentías de 14, 23 Y 18.5 mm., respectivamente. Determinar
el índice de la cuenca.
1 2 4 3
2 6 9 8
3 7 15 11
4 10 12 4
5 5 5 12
6 4 3
7 4
8 2
Problema 4.3
Durante una tormenta, las láminas acumuladas de lluvia en sucesivos pe
ríodos de 2 minutos son: 0.2, 0.5, 0.8, 1.4, 1.8 Y 2.1 cm. Dibujar
la curva masa y el histograma.
¿Cuánto vale la lámina de lluvia neta si la pérdida uniforme (índice
es igual a 10.2 cm/h? ¿Cuánto vale el coeficiente de escorrentía C?
71
viene de la página 42
72
CAPITULO 5 EL AGUA SUBTERRANEA
5.1 Descripción
El agua subterránea es de gran importancia, especialmente en aquellos
lugares secos donde el escurrimiento fluvial se reduce mucho en algunas
épocas del Se estima que en Estados Unidos, de toda el agua que
se'usa al año, una sexta parte es agua subterránea, En Lima, por otro
lado, del total de
-
---,
agu'a que
-- -
se- consume un 40%
-
-
Las aguas del subsuelo, como las aguas superficiales, provienen de las
lluvias. No son independientes unas de otras, sino que, por el contra-
rio, están muy ligadas entre sí. Muchas corrientes superficiales reci-
ben agua del subsuelo y, a su vez, el agua del subsuelo se realimenta
de las aguas superficiales.
Veamos un esquema de las condiciones del agua subterránea.
73
Es necesaria la presencia de un estrato impermeable (1').. Las corrien-
tes superficiales pueden ser afluentes (2) o efluentes (3). Debajo de
la superficie, los poros del suelo contienen agua y en .cantidades
variables'; es la zona vados a (4); en ella la presión es menor que la at
mosfé-rica. Después de una lluvia el agua puede moverse hacia abajo a
través de esta zona de aireación; una parte del agua que penetra es re-
tenida por fuerzas de capilaridad y fuerzas moleculares; el resto sigue
bajando hasta la zona de agua subterránea (5); allí la presión es mayor
que la atmosférica y el agua escurre siguiendo las leyes de la hidráu1i
ca. El nivel superior del agua del subsuelo constituye el nivel freátl
ca (6). A ese nivel se presenta un cordón capilar (7), en el cual los
poros del suelo contienen agua que ha ascendido desde el agua subterrá-
nea por la acción capilar.
Acuíferos no confinados
Una formación como la representada en la figura 5.1 constituye un acuí-
fero no confinado. Si se perforan pozos de observación hasta el estra-
to impermeable, el lugar geométrico de los niveles alcanzados es el ni
ve1 freático (figura 5.2).
---
-- --
¡ ,
\1
Zona vadosa
"-
N,F.
----
Acuífero
Acuíferos confinados
Son acuíferos comprendidos entre dos estratos impermeables (figura 5.4).
El flujo es a presión, como en las tuberías.
;--......
Impermeable
7fJ
En vez de un nivel freático se tiene ahora un nivel piezométrico. La
línea de energía, COfilO en el caso de los acuíferos no confinados, se
confunde prácticamente con el nivel piezométrico debido a que la
de velocidad del agua es muy pequeña.
Los acuíferos confinados presentan las ventajas de conducir el agua a
grandes distancias y entregar el agua por encima del nivel del acuífe-
ro, y las desventajas de tener áreas de recarga relativamente pequeñas,
rendi r menos agua y provocar asentami entos de 1 terreno en los 1uga res
de extracción (pozos de bombeo).
_/
de hoyuelo el e \/I,:j.
de coniente
de onticlinol
Arcilla 45 3
Arena 35 25
Grava 25 22
Grava y arena 20 16
Arenisca 15 8
Calizas densas y esquistos 5 2
Cuarcita y granito 1 0.5
Ley de Darcy
Fue Darcy (1856) quien confirmó que,' con excepclon de las grandes cave!'.
nas o fisuras, el agua del subsuelo escurre siempre con movimiento lami
nar. Aceptando las hipótesis del flujo y uniformemente
distribuido en espesor, la ecuación de Darcy se expresa.:
v = Kp. s (5.1 )
1
v' = v;- v = -vp (5.2)
vr
en que p es la porosidad del suelo. Como p es siempre menor que1;v ' es
siempre mayor que v.
Kp depende de las propiedades del líquido y del medio poroso, y se pue-
de expresar como:
Kp = K ]J
(5.3)
(5.4)
Penneabi 1i dad Kp
Material unidades P. intrínseca K
meinzer m/día Darcys
Arcilla 0.01 0.0004 0.0005
Arena 1,000 41 50
Grava 100,000 4,100 5,000
Grava y arena 10,000 410 500
Arenisca 100 4.1 5
Calizas densas y esquistos 1 0.041 0.05
Cuarcita y granito 0.01 0.0004 0.0005
78
Transmisividad.- Llamando Y al espesor del acuífero y B a su ancho, se
puede escrilJi'r:
Q = A . v = A Kp s = B Y Kp s = B T s (5.5)
El producto Kp . Y se reemplaza muchas veces por un único término T qU2
representa la transmisividad del acuífero. Sus dimensiones son L2 /T,
2
por ejemplo m /día. Si Kp se expresa en unidades meinzer, T resulta en
gal/día/pie. Se deduce que 1 m2 /día = 80.5 gal/día/pie.
Deteminación de la permeabilidad.
Hay dos formas de determinar el valor de Kp de un determinado suelo: en
el 1aboratori o y en el campo.
a) En el laboratorio se usan los permeámetros:
A
r-
±L
o= Av
Kp =
= A Kp s
. L
(5.6)
79
300 m
con 10 que:
Q = A v = 10.25 x 0.005 = O.05125_m J /dia por m.1. de canal.
80
suelo es un limo arenoso con Kp = 2,5 m/día; el fondo impermeable queda
a una profund i dad uniforme de 6 m,
A fin de reducir la contaminación de la corriente, el efluente de una
planta de tratamiento no vaciado directamente sino rociado sobre
el terreno a cierta distancia de ella. Después de la infiltración el
hacia abajo como flujo subterráneo y drenará hacia la
corriente.· El flujo subterráneo y la infiltración previa mejoran consi
derablemente la calidad del efluente con lo que la polución de la
rriente disminuye en alto grado. El sistema deberá ser disenado y ope-
rado de modo que se suprima la escorrentía superficial. Si la aplica-
ción de los rociadores es de 2 cm/día, ¿cuál será el mayor ancho W del
area que ser rociada al mismo tiempo?
2 cm/día
Kp
81
del nivel freático en la cima de la colina?
-
1,. _ _
dh vt'¿ry
= Kp . dx . \,1,..f'.
= _ Kp . h dh
. dx
dQ = P . dx Qx = P . ) (5.7)
1- ¡:-:. -1,_
Igualando las dos expresiones de Qx:
P. x dh
= - Kp . h . dx
- Kp . h . dh = P . x . dx
Integrando entre la cima de la colina y el borde de la corriente:
. 1 L
= P -2 (X2)
o
Kp (hf - = P L2
82
= h2 + P L2
(5.8)
2
K"P
.
1
y di
y Kp
83
Q = Av = A. Kp s
. dx
Q x = 2 TI Kp . Y • dy
2 TI Kp Y (d 2 - dI) (5.9)
Q=
L-¡ rI
_ _ r-
t "2
·1,---
J r-----
I - - - - -N.F.
-
Q = A v = A • Kp s =' 2 TI X Y • Kp
Q = 2 TI Kp . y dy
QL
rI
.. = TI Kp (d/ - di)·
(5.10)
Se acostumbra simplificar:
84
I
Q= (5.11)
85
s =
v
Ad
v= dSA
aV _ s A
af- - -ª-.S!.
at
V volumen de agua desplazada por area horizontal A del
acuífero
d altura de la superficie piezométrica sobre el borde infe
rior del acuífero
S constante de almacenamiento del acuífero
A área horizontal del acuífero a la cual se aplica
t tiempo
El signo negativo corresponde al hecho de que d disminuye canfor
me aumenta t.
Para una area elemental anular, a la distancia r del pozo:
aV = 2TIr. dr
af at
Pero El.
at = - -ª.Q.
ar • dr, correspondiendo el signo negativo porque Q
crece con r decreciente.
Reemplazando:
E..Q dr r . dr
ar
=
at S 2
-ª-º-
ar = at
(5.12)
(5.13)
S 2 TI r
atar
= 2 TI T (i9.. + r
ar2
S ad (5.14)
I 3f
86
Acuíferos confinados.- Theis, en 1935, sugirió una solución pa-
ra la ecuación 5.14 basada en la analogía de transmisión del ca-
lor. Su fórmula es:
Zr =
Q . W( u t (5.15)
4nT
00 -u u2 u3
W(u)=f e - du =-O.5772-Lu+u- 2 21+3x3! (5.17)
u u n x.
(5.18)
87
TABLA 5.3 VALORES DE W(u) PARA DIVERSOS VALORES DE u
Ejemplo 5.4
Se desea calcular la caída de la superficie piezométrica a las
distancias 100 m y.200 m de un pozo de bombeo, para un acuífero
confinado con T = 1,000 m2 jdía y S = 0.0001. El pozo es bombea-
do por la días a un ritmo de 1000 m3 jdía.
En la tabla 5.4 se ha dado solución al problema, a partir de un
juego de valores de t y siguiendo el camino recién señalado.
88
TABLA 5.4 SOLUCION DEL EJEMPLO 5.4
t r = 100 r = 200
días u W{u) Zr u W{u) Zr
y la 5.18 como:
r2 4T
(5.20)
log T = 10g S + 10g u
Q 4T
se puede ver que, desde que tr;I y S son cons tantes en un en-
r2
sayo determinado, la relación entre 10g Zr y log T debe ser si
mi1ar a la relación entre W(u) y u.
Así, si se p10tea Zr contra r 2 /t y W(u) contra u en papel doble
logarítmico, las curvas resultantes serán de la misma forma pero
horizontal y verticalmente desfasadas por las constantes
Q 4T .
YS'
89
Ejemplo 5.5
Hallar T Y S en el siguiente ensayo de un acuífero confinado:
Q = 1,000 m3 /día
rl = 100 m
r2 = 200 m
t (días) 0.001 0.005 0.01 0.05 0.1 0.5 1 5 10
Zn (m) 0.083 0.196 0.249 0.376 0.431 0.559 0.614 0.742 0.797
Zr2 (m) 0.017 0.097 0.145 0.267 0:322 0.449 0.504 0.632 0.687
En primer lugar se confecciona la tabla 5.5
90
111
...
O
Lo.
CD
E er=IOOM
o r =200M
Lo.
N
10 3 I(J
,2ft en m 2 f d(o
111
oLo.
I
IÓ
1
c:
Q)
...
N
10"4 I 1 u 10
3 '1 !5
10 10 10
r 2 /t en m2 f dia
91
Reanplazando los valores de Zr, (W(u) en la eCllaClon 5.15 se ob-
tiene T = 1,000 m2 jdía. Reemplazando los valores de u, r 2 jt, T
en la ecuación 5.18 se obtiene S = 0.0001.
Nótese que los dos ejemplos 5.4 y 5 5 se refieren al mismo ensa-
yo, a fin de comprobar resultados. En el ejemplo 5.4 se conocen
T, S Y se determinan los abatimientos. En el ejemplo 5.5 se
usan esos abatimientos para encontrar T, S. En la práctica, los
abatimientos se miden en el terreno, en los pozos de observación.
Acuíferos no confinados
,
abatimiento
resul tante
//",--, ........ -
Para 1 m de agua dulce por encima del nivel del mar, la ecuación
de equilibrio hidrostático se escribe:
Yl h1 = Y2 h2
1.00 (1 + y) = 1.025 Y
1 +y = 1.025 Y
y = 40 m '
93
FIG. 5.16 CONO INVERTIDO
Ecuación de balance
5 l' + (P - Qs - E) + Qg - R = 52 (5.21)
s almacenamiento
P precipitación en el área tributaria del acuífero
Qs escorrentía en la misma área
94
E evapotranspiración para la misma area
Og agua subterránea neta hacia el acuífero
R rendímiento seguro
Todos los términos pueden referirse a valores medios anuales.
5.6.5 Compresibilidad
Los acuíferos confinados presentan alta compresibilidad. El bom
beo provoca un alivio en la presión interior y su resultado pue-
de ser una compresión del acuífero acompañada de un hundimiento
de la superficie del terreno, veces considerable.
95
5.7 Problemas
Problema 5.1
En la estación A, la elevación del nivel de agua es de 642 pies sobre
el nivel del mar. En la estación B, el nivel es de 629 pies. Las esta
ciones están a una distancia de pies. La permeabilidad del
fero es de 300 unidades meinzer y la porosidad es de 14%. ¿Cuál es la
velocidad del flujo en el acuífero?
Problema 5.2
Suponga que hay dos canales, a diferentes niveles, separados por una
franja de terreno de 1,000 ro de ancho, como indica la figura 5.17. La
permeabilidad es de 12 m/día. Un canal está a 2 m por encima del otro
y la profundidad del acuífero es de 20 m debajo del canal inferior has-
ta el estrato impermeable. Encontrar el caudal que entra o sale de.ca-
da canal por metro de longitud. Considerar una precipitación anual de
1.20 m y asumir una infiltración del 60%.
2m.
Kp=12m/dío 20m.
1. L = 1,0 O O m. */
FIG. 5.17 DATOS DEL PROBLEMA 5.2
Problema 5.3
Un pozo de 12 pulgadas de diámetro penetra 80 pies por debajo de la ta-
bla de agua estática. Después de 24 horas de bombeo a 1,100 gal/min,
el nivel freático en un pozo de observación a uné distancia de 320 pies
desciende 1.77 pies, y en otro pozo a 110 pies de distancia desciende
3.65 pies. ¿Cuál es la transmisividad del acuífero?
Problema 5.4
El registro de abatimiento versus tiempo para un pozo de observación a
96
296 pies de un pozo de bombeo (500 gal/min) se tabula abajo. Encontrar
la transmisividad y la constante de almacenamiento del acuífero. Util i
zar el método de Theis.
97
viene de la página 72
98
CAPITULO 6 EL CAUDAL
6.1 La de Descarga
Para llegar a conocer los recursos de una cuenca es necesa-
rio averiguar el caudal, diariamente, a la misma hora, y durante el ma-
yor número posible de años. Así es como se llega a conocer el régimen
de los ríos. Todos los países cuidan de organizar este servicio, esta-
bleciendo estaciones de aforo y publicando los En el Perú
esta labor la realiza principalmente Senamhi.
Los términos caudal, gasto y descarga son sinónimos. Aforar significa
medir caudales. El principal método para aforar corrientes naturales
es el del correntómetro, el cual es descrito en el apartado siguiente.
Después de seleccionar adecuadamente la sección del río, se establece
la sección de aforo y se procede a medir diariamente el caudal; también
se mide el nivel. Luego de un tiempo es posible dibujar la curva de
descarga del río en el lugar de la estación. Es una curva de caudales
versus niveles o alturas de agua. Se usa en proyectos.
Los niveles se miden con limnímetros o limnígrafos instalados a un cos-
tado de la estación de aforo.
Dibujada la curva de descarga pueden suspenderse los aforos directos,
pues entonces con medir el nivel para conocer el caudal. Se re
comienda revisar periódicamente la curva de descarga con mediciones
rectas de caudal.
a
fIG. 6.1 CURVA DE DESCARGA
99
v = a n+b
v velocidad en el punto
R
n nGmerode revoluciones por segundo = t
a,b constantes de calibración.
Para iniciar un aforo es necesario dividir la secci6n transversal (área
majada) en franjas, como indica la figura 6.2, usando verticales.
l' 4
\1/
CAJA DE
CONTEO
BATERIA
CABLE DE
S'USPENSION
CAMARA DE
ALETAS
HELlCE DIRECCIONALES
LASTRE
101
PONTfFJCJA ()NTVRRsln", n
CATOLlf' A UV!. "EHU
BIBLIOTECA
I
I
INGENIERIA
Según la magnitud de la corriente se hace trabajar el correntómetro sus
pendido de un cable o sujeto a una barra que se hinca en el lecho. La
figura 6.5 corresponde a la primera modalidad.
El cable es para mantener el aparato suspendido desde un puente o una
oroya. El lastre es para impedir que sea sacado de posición por la
fuerza de la corriente. En el eje de la hélice hay una serie de finos
engranajes para poder contar el número de revoluciones. La pequeña c!
mara de contacto hace el cambio de 10 revoluciones a una señal luminosa
y otra auditiva. De esta manera lo único que hace el operario es con-
tar el número de señales en un tiempo arbitrario, a fin de obtener n
(número de revoluciones por segundo) en cada puesta en estación del
rato.
Las corrientes moderadas son vadeables. En ese caso se usa la barra,
debiendo el operario hacerse a un lado a fin de no interrumpir la co-
rriente que va a ser registrada.
A = í.!lA = 3.558 m2
Q = í.!lQ = 1.034 m3 /seg
V = Q = 0.29 m/seg
A
102
Unidades
La unidad básica de flujo es el m3 /seg. El volumen de flujo se puede
expresar en m3 , pero como esto lleva a números demasiado grandes se
acostumbra expresar en miles de m (m MC) o en millones de m
3 3
Los caudales pueden expresarse también en m3 /seg/km 2 , para comparar ca-
sos dé flujo en ríos con áreas tributarias diferentes, y son iguales al
caudal en m3 1seg dividido entre el área de drenaje en km 2 .
El mm es la cantidad de agua necesaria para cubrir el área de drenaje
con una profundidad de un milímetro; es una unidad de volumen bastante
útil para comparar caudales con la precipitación que ha sido la causa.
Caudales medios
En época de caudales estables s610-es necesario determinar el caudal
(m 3 /seg) una vez al día, siempre a la misma hora. Ese valor es conside
rado el caudal medio diario. En época de variación de caudales es nece
sario determinar el caudal dos o tres veces al día a fin de obtener el
caudal medio diario. Ahora, el promedio mensual de las desca-rgas me-
dias diarias proporciona la descarga media mensual y el promedio de és-
tas la descarga media anual.
Hi drogramas
Reciben el nombre de hidrogramas los gráficos Q-t, en general. Un hi-
drograma de creciente es el hidrograma que corresponde a una crecida
aislada del río por efecto de una tormenta importante en la cuenca co-
lectora (figura 6.6).
Q
t
FIG. 6.6 HIDROGRAMA DE CRECIENTE
E F M A M J J A s o N o
104
TABLA 6.2 n = 0.030 S = 0.00395
26.2 O O
95.0 50.58
30 95.0 50.58 1.88 450.5
305.0 24.57
35 400.0 75.15 5.33 3,799.6
377.5 12.46
40 777 .5 87.61 8.87 10,371.6
410.0 11.66
45 1,187.5 99.27 11.96 19,333.8
442.5 12.46
50 1,630.0 111. 73 14.59 30,298.5
Lr".. _trol
de l. . . .oiÓ"
d. U.ro.
(A·)
( B)
105
lO
••
- -r-
V
y
- I
-. - L
40
,1 JI
z:
O '5
U
111
54
t7
-: II
JO
2'
ZI
o Z 4 • • 10 12 14 16 " 20 Z2 24 26 Z' 50
106
de la tierra, se recomienda utilizar solamente los registros que.repre-
sentan las condiciones actuales y extenderlos en base a correlaclones.
Q
m3/seg
\
r-- t---
.
r--
o 10 20, 30 40 50 60 70 80 90 100 %
70
40
30
20 ____________ _
B
10 ____________ _
O 20 40 60 80 100
Construcción
El método de construcción de la curva de duración que se va a describir
es el método del año calendario. Se ordenan los caudales medios mensua
les para cada año en forma decreciente y se les asigna un número de
den. Luego se promedian los caudales para un mismo número de orden.
Por Gltimo se grafica: caudales en ordenadas y número de orden o proba
bilidad de excedencia en a b s c i s a s . ' -
N!:: de
orden 1 2 3 10 11 12
Año
1972 4.2 3.9 3.6 0.3 0.1
1973 13.8 13.7 13.3 0.2 0.1
1974 4.5 4.1 3.8 0.4 0.2
1975 12.8 10.6 9.9 0.5 0.3
1987
1988
Promedio 15.7 12.2 11.6 0.4 0.2
% 8-.3 16.7 91. 7 100.0
108
6.6 La Curva Masa
La curva masa, llamada también curva de volúmenes acumulados, es una
curva que se utiliza en el estudio de regularización de los ríos por
dio de emba'lses. Proporciona el volumen acumulado que ha escurrido en
una estación en función del tiempo, a partir de un origen arbitrario.
Es por ello una curva siempre creciente, que contiene a lo más pequeños
tramos horizontal.es o casi horizontales correspondientes a los meses se
coso
Supondremos, para los efectos de explicación, que se ha dibujado la
va masa para los tres años de mayor irregularidad dentro del tiempo de
registros del río (figura 6.12). La idea es estar prevenidos en caso
se presente más adelante un período crítico como éste.
Vo.lumen
acumulado
MMC
A
1977 1978 1979 meses
e
MMC
meses
MMC
_ _-
meses
6.7 Problemas
Problema 6.1
Calcule el caudal con la información dada en la tabla de abajo. Suponga
que la calibración del medidor es de la forma v = a+bn, con a = 0.1 Y
b = 2.2 para v en pie/seg.
Distancia Profundidad
desde la Profundidad del corren- Revo 1u- Tiempo
ori 11 a (pies) tómetro ciones (seg)
(pies) (pies)
2 1 0.6 10 50
4 3.5 2.8 22 55
0.7 35 52
6 5.2 4.2 28 53
1.0 40 58
9 6.3 5.0 32 58
1.3 45 60
11 4.4 3.5 28 45
0.9 33 46
13 2.2 1.3 22 50
15 0.8 0.5 12 49
17 O
111
Problema 6.2
-
A continuación se presentan las descargas medias diarias en metros cúbi
cos por segundo en una estación de medición para un período de 5 días.
¿Cuál es el caudal medio para el períOdo en metros cúbicos por segundo?
¿Cuál es el volumen total durante el período en metros cúbicos? Si el
área tributaria es de 100,000 km 2 , ¿cuál es la lámina de escorrentía
equivalente en mm?
Día 1 2 3 4
Problema 6.3.
Dibujar la curva de duración para los datos Ge la tabla de abajo. Las
cifras son caudales medios mensuales en m3 /seg.
Año 1 Año 2 Año 3
Enero 110 180 193
Febrero 102 118 109
Marzo 97 88 99
Abri 1 84 79 91
70 56 82
Junio 62 52 74
Julio 45 47 68
Agosto 67 35 43
Setiembre 82 60 30
Octubre 134 75 48
Noviembre 205 98 49
Diciembre 142 .127 63
Problema 6.4
La figura representa un hidrograma simplificado.
Calcular y dibujar la curva masa.
Q
m3/seg.
10
1,
I
I
I
o 3 9 18 t horas
112
Prob 1em a 6. 5
Una corriente proporciona los siguientes volúmenés en un período de 80
días en el lugar de un posible reservorio. a) Dibujar la CUY'va masa.
b) Determinar los caudales medio, máximo y mínimo. c) capacidad
de reservorio se necesita para asegurar un caudal regulado igual al
dal medio del período si el reservario arranca el período estando lle-
no? d) ¿Qué cantidad de agua se perdería en este caso por el alivia-
dero de demasías del embalse?
113
viene de la página 98
114
CAPITULO 7 RELACIONES PRECIPITACION-ESCORRENTIA
7.1 Introducción
En el apartado 4.1 se describió cómo el agua de un río, en general,
de estar formada de dos partes: una parte de escorrentía directa y otra
parte de agua subterránea. Si bien ambas provienen de las lluvias, só-
lo la primera obedece a las precipitaciones recientes.
El poder inferir el caudal proveniente de una precipitación tiene múlti
ples aplicaciones. Por ejemplo, permite obtener los caudales en un no
sin estaciones hidrométricas; o extender los registros cortos de cauda-
les a fin de someterlos a análisis
Por éstas y otras razones, un problema clásico en Hidrología está cons-
tituido por la obtención de la escorrentía directa que corresponde a
una determinada lluvia, en un lugar específico. El primer método es a
través del coeficiente de escorrentía C (apartado 4.1). El segundo mé-
todo es mediante la separación en el histograma usando la curva de in-
filtración (apartado 4.2). El tercer método consiste en el empleo de
los índices de infiltración, de los cuales el índice es el más conoci
do (apartado 4.6). Existen todavía otros métodos, como el que usa los
datos de suelos y cubierta vegetal, el método racional y los métodos de
simulación por computadora; de éstos serán descritos en los apartados
que siguen los dos primeros (7.2 y 7.3). El último escapa a los alcan-
ces del texto. Se hace notar que todos los métodos reseñados son para
el cálculo de la escorrentía por tormenta individual; en la práctica se
requiere también el cálculo para períodos largos de tiempo (mensual o
anual), 10 cual es descrito en el apartado 7.4.
115
hidrogramas de descarga de tres cuencas de diferente forma e igual área
sometidas a una lluvia de igual intensidad (figura 7.1). Si cada cuen-
ca se divide en segmentos concéntricos, que se puede asumir tengan to-
dos los puntos a la misma distancia del punto de salida de la cuenca,
se puede ver que la forma A requerirá 10 unidades de tiempo antes que
todos los puntos de la cuenca estén contribuyendo a la descarga. Simi-
larmente B requerirá 5 y C 8 1/2. Los hidrogramas de escorrentía resul
tantes serán similares a los mostrados en la figura 7.1, cada uno
do con la correspondiente letra minúscula. La forma B da una corriente
de ascenso más rápido que las forma C y A, y también de descenso más ra
pi do.
,.,,""
'\5-·--
y'
- -1,.
,4, .. ---. ,
I I
I I \2 --1"
"'3
I
I
I ,.
", r,'
:É':' 1L3-1'?i.!
A B e
?-V
ti
e l/
VVc
...e
C>
I
J/ V
u
'"
cu
O
.,/
O 5 10 15
Tiempo (horas)
Otros.- Aparte de los citados hay otros factores que afectan la tasa y
cantidad de la escorrentía, como el déficit de humedad del suelo, la al
titud (con su efecto sobre la temperatura y la presencia de nieve en in
116
vierno), el uso de la tierra (ya sea área de bosques o tierras de culti
va), la proporción del desarrollo urbano, etc.
Barbecho SR 77 86 91 94
Cultivos en hil eras SR Mala 72 81 88 91
SR Buena 67 78 85 89
C Mala 70 79 84 88
C Buena 65 75 82 86
Cy T rv1a 1a 66 74 80 82
Cy T Buena 62 71 78 . 81
Granos pequeños SR Mala 65 76 84 88
SR Buena 63 75 83 87
C Mala 63 74 82 85
C Buena 61 73 81 84
Cy T Mala 61 72 79 82
Cy T Buena 59 70 78 81
Legumbres tupidas SR Mala 66 77 85 89
o rotación de SR Buena 58 72 81 85
pradera C Mala 64 75 83 85
C Buena 55 69 78 83
Cy T Mala 63 73 80 83
Cy T Buena 51 67 76 80
Pradera o pastizal Mala 68 79 86 89
Regul ar 49 69 79 84
Buena 39 61 74 80
C Mala 47 67 81 88
C Regul ar 25 59 75 83
C Buena 6 35 70 79
Pradera perm. 30 58 71 78
Bosques (lotes de Mala 45 66 77 83
bosque) Regul ar 36 60 73 79
Buena 25 55 70 77
Cascos ranchos 59 74 82 86
Cami nos reves t. 72 82 87 89
Pavimentos 74 84 90 92
SR = hileras rectas
C = por líneas de nivel
T = terrazas
Cy T = terrazas a nivel
119
TABLA 7.2 NUMEROS DE LAS CURVAS DE ESCURRIMIENTO PARA LOS
COMPLEJOS SUELO-VEGETACION
1 La peor 56 75 86 91
II Mala 46 68 78 84
III Media 36 60 70 76
IV Buena 26 52 62 69
V Mejor 15 44 54 61
Grupos de suelos
Vegetación Condición
A B C D
Herbácea Mala 90 94 97
Regular 84 92 95
Buena 77 86 93
Artemisia Mala 81 90
Regul ar 66 83
Buena 55 66
Roble-Tiemblo Mala 80 86
Regular 60 73
Buena 50 60
Junípero Mala 87 93
Regular 73 85
Buena 60 77
120
Condición III Que se présenta cuando ha llovido mucho o poco y han
ocurrido bajas temperaturas durante los cinco días ante-
riores a ]a tormenta, y el suelo está casí saturado.
Los números de la tabla 7.1 y de la tabla 7.2-A son para la condición
media de la cuenca, 11. Los números de la tabla 7.2-B son la con-
dición casi saturada, 111. Los números de las curvas para una condi-
ción precedente pueden convertirse a una condición usando la
tabla 7.4. Por ejemplo, el cálculo dado en la tabla 7.3 da una curva
número 73 en una condición 11.. Por interpolación se pueden obtener los
números de las curvas .correspondientes a la condición I y a la condi-
ción III de las columnas 2 y 3 de la tabla 7.4. Los números de las
vas para la condición I y para la condición 111 son 55 y 89, respectiv,ª-
mente.
Determinación de lá'escorrentía directa
Sea Q escorrentía dir-ecta en lámina de agua, en pulg
P precipitación en lámina de agua, en'pulg
S diferencia potencial máxima entre P y Q, a la que co-
mienza la tormenta.
Por mediciones hechas en cuencas naturales se sabe que Q se aproxima a
P, mientras P aumenta en la tormenta. También que los valores (P-Q) se
aproximan a una constante mientras P continúa aumentando. Las cantida-
des pueden agruparse en la forma:
P-:Q = Q. (7.1 )
'S P
Número representativo =
7,"337 = 73.37 -+ 73
100
1 2 3 4 5
Número de la Números correspondien La curva
curva pa ra 1a tes de la curva para: Va lores comienza
condici6n II S* donde P = *
Condici6n Condición
I 1 II
FIG'. 7.2
SIGNIFICADO DE la
lo
tiempo
Despejando S:
s = 1000 - 10 N (7.7)
N
Ejemplo 7.l
Averiguar la lámina de directa que es de esperar ocurra en
una cuenca natural con número de curva representativa 60, como conse-
cuencia de una tormenta de 20 pulg.
Q caudal de escorrentía
e coeficiente de (tablas 4.1 y 4.2)
i intensidad máxima de la lluvia para un perfodo de duración
igual al tiempo de concefltración, y para la frecuencia desea-
da en el diseño.
A ... área de la cuenca.
Si i está en- m/seg y A en mL , Q resulta en m3 iseg. Si i está.·"en rnm/h y
A en Ha, entonces Q en m3 /seg viene dado pór:
C·i A
Q = 360 (7.9)
TABLA 7.5
125
7.4 Correlaciones Precipitación-Escorrentía
La correlación lluvia-escorrentía más simple es la representación gráfi
ca de dos variables: lluvias promedio contra escurrimientos resultantes
(figura 7.3). La relación típica es una curva ligera que indica un in-
cremento en el porcentaje de escurrimiento con las mayores lluvias. Es
tas relaciones simples no toman en cuenta las condiciones iniciales que
afectan el escurrimiento y, generalmente, hay una dispersión considera-
ble de los puntos con respecto a la curva media.
..
...o
4
8.
:;
Q.
.
...oo 3
o
..
:o
e
!!
."
..
o
Q.
Q.
.2
>
." >
:J ."
::J
o 0·!5 1.0' 1.!5 2.0 o 0.4 0.8 1.2 1.8 2.0 2.4
Puede introducirse una tercera variable para dar una explicación a las
desviaciones que hay en la relación simple. En regiones húmedas, el es
currimiento inicial en la corriente refleja condiciones antecedentes
con bastante claridad y puede servir como un parámetro efectivo (figu-
ra 7.4). .
Otro tipo de correlación lluvia-escorrentía tjene la forma:
M(D) = a. M(A)S
M(D) escorrentí a medi a anual, en mm
M(A) precipitación media anual, en mm
a. , S coefi ci entes propi os de cada cuenca.
Este es el tipo de correlación aplicado en el Estudio de la Hidrología
del Perú (capítulo 11).
Otro tipo de correlación, aplicable en regiones húmedas, es el que se
explica en el ejemplo 7.2 y que tiene la forma:
=
126
Dn es correntia media en el periodo n
An precipitación media en el periodo n
Al"
n· precipitación media en el período n-1
a,b ... coeficientes propios de la cuenca tales que a+b = 1
Ejemplo 7.2
Los datos de la tabla 7.6 se refieren a precipitaciones medias anuales,
en mm, y medios anuales, en m3 /seg, para una cuenca de superfi
2
cie 458 km • Usando el modelo matemático de correlación:
Dn = a An- 1 + b An
a+b = 1
determinar los caudales correspondientes a los años sin datos de cauda 1.
Metodología.-
12 Con a=l, b=O y usando el modelo, hallar para' cada año del periodo
1973-1982 el' valor del escurrimiento O, mm.
22 Repetir varias veces, disminuyendo a y aumentando b cada vez en 0.1
32 El valor correcto de a es aquel que hace minima la suma de las des-
viaciones cuadráticas de Dn y DQ:
o
6
Perfil de la s....,.rficie del 0VUG
( de las puntos de 0vuos mó"imOl
'Elav.!!" 1\ u
i \Ele". 5.47 1°
a .. " o
¡- Morco. de .... ClQuos
4 .d.,lIIII. (nO ro
¡'"- de le47
_:t lI!lIO!
(1)
3
•• ",ti
.Q
11.1
¡
Z
t
I
"Perfil aproa/,. adO
de' fondo
,
I
-
11.1
O
I !! O
Z
4410 4Z0
lo ¡.toneio en
380 340 300 Z60
p io • • la 1000o de la línea central del COUC4I
Z20 110
O 7
IJ
i"l-\-L_ Nivel dal a Ilila 5.69 r I
f-1
/"
el: 6
>
11.1
--- --- Nivel del 00UG 5.47
i--.L - - + - - -
...J 5 \
11.1
04 I \ !/
\ \... l,P-o' \L I
"-
3
-
$acc.on F Secciot. E
O
O 10 ZO 30 040 O 10 ZO 30 40
DISTANCIA - EN PIES
FIG. 7.5 DATOS DEL EJEMPLO 7.3
FIG. 7.6
- G
SIN
G + hV1 = hV 2 + h amp + hf
G + h. 1 = hV 2 + 0.5 (hv 1 - hv 2) + hf
hf = G + 0.5 (hv1 - hv 2)
\
Para proseguir con el ejemplo, se van suponiendo diferentes gastos en
diferentes tanteos hasta que los gastos supuesto y calculado sean igua-
les .. Esto se ilustra en la tabla adjunta. Se usa' el valor medio de K
de 1 primer tanteo (7035). .
Sección
Suponer Q = 460 pcs
E 5.J2 0.51
F 5.62 0.49 -0.02 0.20 0.00408 449
Suponer Q = 450 pcs
E 5.60 0.49
F 5.49 0.47 -0.02 0.20 0.00408 449
130
TABLA 7.7 COEFICIENTES DE RUGOSIDAD (n) PARA CAUCES NATURALES
131
7.6 Problemas
Problema 7.1
Se trata de averiguar la lámina de escorrentía directa que es de espe-
rar ocurra como consecuencia de una tormenta de 25 cm, una cuenca de
las siguientes características:
Porcentaje del área 60
Grupo hidrológico del suelo A B
Cubierta Granos en hileras Pradera
Infiltración Buena Buena
Problema 7.2
Calcular el caudal máximo escorrentía es de esperar ocurra una
vez cada 10 años, en una cuenca de 1.2 km de superficie y 3% de pen-
diente media. La longitud del curso principal es de 2,800 m. Del to-
tal del área el 40% corresponde a cultivos generales y el 60% a pastiza
.les. las curvas de intensidad.duración-frecuencia del
2.7.
CAPITULO 8 HIDROGRAMAS DE CRECIDAS
8.1 Introducci ón
El caudal de una corriente, en general, está constituido de dos partes.
Una de ellas, el flujo base, proviene del agua subterránea y la otra,
la escorrentía directa, proviene de las últimas lluvias. No todas las
corri entes reci ben aporte de agua subterránea, nT-todas 1as preci pita-
tiones provocan escorrentía directa. Sólo las precipitaciones
tes, es decir, intensas y prolongadas, producen un aumento significati-
vo en la escorrentía de las corrientes. La contribución de agua subte-
rránea a las corrientes de agua no puede fluctuar rápidamente debido a
la baja velocidad del flujo.
Las corrientes en cuenca con suelos permeables, y que reciben gran apor
te de agua subterránea, muestran caudales altos sostenidos a lo largo
del año, con una relación baja entre caudales de avenidas (crecidas) . y
caudales medios. Las corrientes en cuencas con suelos de baja permeabi
lidad, y que mas bien aportan agua a los acuíferos,
nes y promedlo, con -caudales muy bajos o nu-
los entre crecientes. El hidrograma A de la figura 8.1 corresponde a
las corrientes del primer tipo, y el hidrograma S a las del segundo ti-
po. Nuestros ríos que desembocan en el Pacífico tienen características
del tipo B.
Q
133
8.2 El Hidrograma Típico
El hidrograma típico de una tormenta aislada (figura 8.2) consta de una
rama ascendente, un segmento de cres ta y una rama descendente o curva
de recesión.
Q romo ramo
ascendente d"c.ndon"
ESCORRENTIA
DIRECTA
AGUA SUBTERRANEA
t
FIG. 8.2 EL HIDROGRAMA TIPICO
La forma de la rama ascendente está influenciada sobre todo por las ca-
racterísticas de la lluvia que causa el ascenso. La forma de la rece-
sión en camhio es bastante independiente de ello y más bien depende de
las características de la cuenca (apartado 7.1). Se asume por 10 gene-
ra l' que e 1 punto de i nfl ex ión de 1a curva de reces i ón coi nci de con el
tiempo al cabo del cual cesa la escorrentía superficial hacia los cur-
sos; de ahí en adelante la curva representa el aporte de agua almacena-
da dentro de la cuenca. El último tramo de la curva de recesión repre-
senta casi completamente el flujo de agua subterránea.
Separación en el hidrograma
En un hidrograma de crecida hay necesidad de separar 10 que es escorren
tía directa y 10 que es flujo base. No existe una forma única
la separación, y puesto que las definiciones de las dos componentes sQP
un tanto arbitrari as los métodos de separación son tambi éil arbitrar; os.
N
Q
ESCORREI\ITlA
DIRECTA
FLUJO BASE
I
I
.'-
" II ,
./
----
- - - --------=- ------=-
--- ... _--.._- ---- ---" ---=- ::-
8
t
FlG. 8.4 METODOS DE SEPARACION
Q
I cm
Vo = Icm
A
- - H. U. de las tI horos
T t
t.
A
/' h,
Hidrograma resultanfe
de lo lluvia h2, h 5.
h cm
1, ti t
L
1
--r---- JI cm.,
t
/H.U de los
139
8.3.3 La Curva S
la curva S de una cuenca se dibuja a partir del H.U. de las t]
horas y sirve para obtener el H.U. de las t 2 horas. Aquí radicé
su enorme importancia: pennite derivar hjdrogramas unitarios é
partir de uno conocido.
Se llama curva S (figura 8.10) el hidrograma de escorrentía di-
recta que es generado por una lluvia continua uniforme de dura-
ción infinita. La lluvia continua puede considerarse formada de
una serie infinita de lluvias de período p tal que cada lluvia
individual tenga una lámina de 1 cm. El efecto de la lluvia con
tinua se halla sumando las ordenadas de una serie infinita de
hidrogramas unitarios de t l horas según el principio de superpo-
sición.
1 1
__ __+-__+-__+-__+-__+-__+-__
Q P
par rodo P
__ __ __+-__+-__+-__ __+-__
______________ __ ________
Q
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _L -_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
hf. x ordenados
de' H.U.
h, x ordenados
de' H.U.
FLUJO BASE
t horas
142
adoptar una duración tipo de lluvia neta (T). Snyder adoptó:
T
T = nL (8.2)
(8.3)
Las ecuaciones 8.1, 8.3, 8.4 definen los tres elementos necesa-
rios pare construir el H.U. para una duración tipo dada por 8.2.
Para cualquier otra duración TD el tiempo de retardo viene dado
por:
TD - T
TLO = TL + 4 (8.5)
143
lOO I
h
te) I \ o t· q*
'I
80
I\k ,o O O
f \ b 3
1\ e 7 1!5
10 I
d 10 35
"\ e 13 56
2., 60 \ f 16 77
o
O'
-3 /JO
e
,
1
m
9
h
i
18 90
20 97
23 100
j 27 96
• 1
lit
, k 30 85
"
.g
40
l\. 1 34 F
e t-- r-d
m 40 50 '
::> 50 \" r--
n 47 33
I r\. o 53 24
20 1 Areq bajo lo cur-
vo::3,300 unida- p 60 16
'\ p 67
10
de volumen.
'"
q
r 84 4 "
f""'-. • 100 O
o
IbJ
1/ '" -o
r
-... t-- • • t Y q en %
O lO 20 30 40 50 60 10 eo 90 100
Unidades de tiempo
Ejemplo 8.1
Elaborar el ridrograma aproximado de una cuenca sin aforar, correspon-
diente a una crecida cuyo caudal pico es 17.64 m3 /seg y su respectivo
volumen 677.000 m3 •
generales:
V
u = 3,300 m3 /unidad (8.6 )
u
t·- 60 q min/unidad (8.8)
10
i ttt
j, o q"
h -
I 1'\ o O Q
15 I9F k b 0.9
c 137 2..6
1/ 1\ d 19S 6.2
\ e 2!5S 9.9
ÓI
..,• I
, f
9
315 13.6
382 IS.9
E
10 I .. I h
i
392 17.1
4!50 17.6
"'O
leS'
\ m j 530 18.9
"o
::J k sea 15.0
u \ 1 66.7 12.7
d \ m 783 8.8
5
r\n n 920 5.8
o 1,040 4,2
1 !'bo p 1,180 2.8
/ A ec balo la curvo
- 1 -1--
q 1,320 1.9
c1 567,000 mI Q r 1,6!50 0.7
h I r
* t en minutos
O lf I , • q en m5Jseg.
O 500 IpOO 1,500
Tiempo en minutos
h cm
t horas
145
precipitación netó h cm.
FIG. 8.16
hidrograma de
escorren t(a direc ta HIDROGRAMA
Qp ( h cm) TRIANGULAR
O
________ ___ t.
= 1,800 (T p Qp + T Qp)
r
1Q4A
10 4 A h
=
1',800 '{T p + Tr} = 5.556
T +
Ah
T
P r
Qp = 2.08T A h (8.9)
P
Para el tiempo de retardo se puede usar la relación empírica:
(8.10)
de modo que:
(8.11)
147
484 A h (8.9a)
Tp
Tb=T p +T r =T p
+aT = (1+a)=2.67T
p . p =18.7horas
18
16
I
/' :--....
\
14
J7Q.62 / 1\
12
( lIu io ne o
I \
10
V
8 / /
/"
1', \
L 1\\
"-
/
v-- / "- ,
6
"'- --
-----
"
4 I v- __ L
j ' .;'
/
'" .....
L ..... " '_ ... .......
"-
....... 1',
,/ / 1"- __ r--- . . . . "-
. ;"
--
,/
/ .........
2 -- r--
.... "- .......
¡"í
/ /
/
,.,/
/
..... ..... " to....
O
O 2 4 6 8 10 12 14 16 18 20 22 24
horas
Q 25
4 hora.
Problema 8.2
Una lluvia constante de 4 horas de duración e intensidad 50 rnm/h produ-
ce un caudal pico de 280 m3 /seg. La tasa de pérdida de la cuenca es
12 mm/h y el flujo base es 20 m3 /seg. Según la teoría del hidrograma
unitario, ¿cuál sería el caudal pico de una lluvia de 4 horas, de ·38
rnm/h, si la tasa de pérdida es 15 mm/hy el flujo base 3 m3 /seg?
Problema 8.3
Una' lluvia uniforme con una intensidad de 50 rnm/h y una duración de 1
hora genera a la salida de cierta cuenca el hidrograma que se muestra
abajo. Calcular el hidrograma causado por una lluvia uniforme de 20
rnm/h y duración 2 horas. El índice es de 10 mm/h.
<1> .
Q
m'/s8g 200
__ ______________ t
2 6 horas
149
Problema 8.4
El hidrograma unitario para una lluvia de 2 horas, de intensidad uni for
me y lámina neta de 10 lTV11 tiene las siguientes ordenadas:
Tiempo (horas) O 1 2 3 4 5 6
Q (m 3 /seg) O 77 155 116 78 38 O
Obtener el hidrograma unitario para una lluvia de 3 horas, de intensi-
dad uniforme y la misma lámina de lluvia neta.
CAPITULO 9 TRANSITO DE AVENIDAS
Las ondas de avenidas se forman por aumento no uniforme del caudal del
río por efecto de una tormenta importante. Para su estudio hay disponi
bles dos métodos: el método hidráulico y el hidrológico. Ambos inten-=-
tan describir los cambios que en el tiempo experimenta la onda de aveni
da. El análisis del paso de estas ondas de avenida constituye materia
de estudio del "tránsito de avenidas" o "flood routing".
El método hidráulico de análisis es bastante complejo por cuanto las
condiciones naturales a que se aplica también lo son: el flujo es no
permanente, la sección transversal es no. uniforme, la rugosidad es va-
riable, etc. El método hidrológico está basado en hipótesis simplifica
torias y consiste básicamente en plantear la ecuación de continuidad
un tramo corto de la corriente.
La ecuación de continuidad, referida a un tramo corto del cauce, puede
escribirse, para un tiempo corto:
Volumen que ingresa - volumen que sale = cambio en el almacenamiento
V.1 - Vo = (9.1)
151
l'.S (9.2)
1 - O = l'.t
1 + 12
1
=
02 S2 01 SI
2 - °1 (2 + M) (2 + M)
11 + 12
2 - °1 = N2 - NI (9.3)
° S
siendo N = 2 + l'.t
o
..
s
o
N
o
F1G. 9.? OBTENC10N DE LA GRAF1CA N-O
152
Para mostrar el procedimiento de cálculo se va a emplear el esquema de
la figura 9.3. Estando el embalse lleno llega el flujo base de la co-
rriente, lo' A este valor inicial del caudal de llegada corresponden
los valores iniciales de nivel de agua, embalsamiento, carga y caudal
de salida. En estas circunstancias se supone llega la avenida
tada por el hidrograma I-t. En este hidrograma se adopta un valor de
y se leen las escorrentías directas 1, con 10 que las cuatro prime-
ras columnas de la tabla 9.1 pueden ser llenadas. El otro valor conocl
do antes· de iniciar los cálculos es la descarga inicial 01. El modus
operandi para llenar la tabla está indicado con flechas en la tabla mis
ma'.
1
Co
Hora I N, O
NI
153
Tránsito en embalses controlados
Para un embalse con compuertas en el vertedero el gasto de salida
de del número de compuertas que estén abiertas (figura 9.4). Entonces
la curva N-O debe ser sustituida por una familia de curvas, en las que
el número de 1 a 3 indica el número de compuertas abiertas.
e c
2 3
s
N
2
o
FIG. 9.4 EMBALSES CONTROLADOS
/1
6
5
j \ I
¡-
en
rn--CD 4 LQC'\ Q
t'
iento ,
'0. " -/
,
,
CII
"O 3 1\ I
i'
\
111
I I
\
,
rdido' ,
E
e 2 IJ I
e olmoc
...
o
CII
I/ ..•.
""
",
"O .,.
"O I .,.
.,.
'--r--- - - -
:::J
o
U
-- -
o
6p.m 6a.m 6p.m. 60.m. 6p.m 60.m. 6p.m. 6 a.m. 6p.m.
4 5 6 7
Horo 1 O
¿s S
6t AT
3/6pm -
560 700 -140
3/12 pm -
660 660 O
4/6om O
2500 1020 1480
4/12 m 1480
5500 1850 3650
4L6pm 5130
5950 2650 3300
41'2 pm 8430'
4200 3350 850
5/6 am 9280
2950 3700 -750
5/12 m 8530
2100 3680 -1580
6950
1470 3100 -1630
5/12pm 5320
1000 2450 - 1450
6/60m 3870
740 2000 - 1260
6/12m 2610
600 1650 -1050
6/6pm 1560
530 1300 - 770
790
155
.1>-
/'
.. 16
,;;-
.
•e 12
I
/
.,/
v
L "
o
(1) V"
4 /
rY'
o
0.2 0.4 0.6 0.8 1.0 I.e 1.4
20
.,
o 16
o
Q
e
ti 12
......
O
)( 8
I
-+
)(
4
L.-J
O
Al macena miento en IO'm 3
156
Si el almacenamiento está dado en m3 /seg x día y los caudales están da-
dos en m3 /seg, K tiene unidades de días.
Si se escribe la ecuación 9.4 para S1 y S2 y se reemplaza en la 9.2a,
se obtiene, resolviéndola para O2 :
O2 = Co 12 + C 11 + C2 01 (9.5)
l
donde:
Co = - K -K xK x- 0.5 lit
+ 0.5
Cl = K -K xK +x 0.5 lit
+ 0.5 lit
K - K x - 0.5 lit
C2 = K .. (9.6)
Se verifica que Co + Cl + C2 = 1
En estas ecuaciones lit es el período de tránsito en las mismas unidades
de K. Con los valores de K, x y t establecidos, se puede calcular los
valores de Co ' Cl y C2' El cálculo del tránsito se reduce a resolver
la ecuación 9.5, con los valores 02 de un período transformándose en
los valores 01 del período siguiente (tabla 9.3).
Hora 1 o
157
9.4 Problemas
Problema 9.1
Dados los hidrogramas tabulados abajo, encontrar el almacenamiento en
el tramo y hacer un gráfico del almacenamiento en cualquier instante co
mo una función del flujo simultáneo de salida. Considere nulo el influ
jo loca l.
1
Fecha Hora pie 3 jseg
O
pie 3 jseg Fecha Hora pie 13 jseg pie O
3 jseg
1 40 40 7 N 245 394
2 N 35 39 M 192 307
M 37 37 8 N 144 235
3 N 125 52 M 118 180
M 340 130 9 N 95 142
4 N 575 287 M 80 114
722 472 10 N 67 93
5 N 740 624 M 56 77
M 673 676 11 N 50 64
6 N 456 638 M 42 55
M 320 574
M medianoche
N mediodía
Problema 9.2
Un embalse pequeño tiene un espejo de agua de 300 acres cuando el agua
está al nivel del aliviadero; las bancas son esencialmente verticales
por encima de este nivel. El aliviadero tiene 15 pies de largo y un
coeficiente de 3.75. Con el hidrograma de entrada del problema ante-
rior, calcular el máximo nivel del embalse y la descarga máxima espera-
da si el embalse está inicialmente a nivel del vertedero a la mediano-
che del día primero (1 acre = 43,560 pies 2 ).
Problema 9.3
Encontrar K y x 'de para el problema 9.1.
Problema 9.4
Utilizando el de salida del problema 9.1 como entrada a un
tramo con K = 27hbras y x = 0.2, encontrar el caudal pico de salida
utilizando el método de Muskingum.
158
CAPITULO 10 ELEMENTOS DE HIDROLOGIA ESTADISTICA
10.1 Introducción
De la referencia 9 extraemos los siguientes comentarios, con ciertos
ajustes para su adaptación al presente texto.
El objetivo básico de la aplicación de la estadística en Hidrología es
el análisis de la información hidrológica en forma de muestras, a fin
de inferir las características con que debe ser esperado en el futuro
el fenómeno que se estudia. El avance en el campo de las computado-
ras y el desarrollo creciente de métodos numéricos han dado una impor-
tancia particular al uso de la estadística en todas las ciencias natu-
rales, especialmente en Hidrología.
Existe en muchos la idea de que la estadística es usada sólo cuando no
es posible dar una solución exacta a un problema hidrológico. En esta
interpretación la solución exacta es una solución determinística' del
problema. Sin embargo, se puede demostrar que la solución determinís-
tica constituye una solución particular de la solución estadística o
probabi lísti ca.
En forma general, la mayoría de los problemas hidrológicos se pueden
agrupar en tres categorías principales de acuerdo al objetivo princi-
pal del proyecto:
a. Diseño de estructuras hidráulicas, siendo necesaria la evaluación
y cuantificación de los valores extremos (máximos y mínimos) del
escurrimiento superficial.
b. Satisfacción de demandas, siendo necesario evaluar y cuantificar
las descargas disponibles en el punto de interés.
c. Diseño y operación de embalses, siendo necesario evaluar y cuanti-
ficar la variación del escurrimiento superficial en todas sus ca-
racterísticas estadísticas, como valores medios, máximos y míni-
mos.
En cada una de las tres categorías menci onadas se presentan diferentes
tipos de problemas, dependiendo la simplicidad o complejidad de la so
lución del tipo, cantidad y calidad de la información disponible, así
como de 1a magni tud del proyecto. Los casos más comunes que se presen
tan en cada una de las tres categorías mencionadas son: --
lo Cuencas con suficiente información hidrológica. Este es el caso
más optimista donde se pueden aplicar todo tipo de metodologías
existentes.
2. Cuencas con escasa información hidrológica. En este caso se pue-
den desarrollar modelos <que relacionen las precipitaciones con las
de?cargas, mediante.el uso de la regresión simple o múltiple, li-
neal o no lineal.
3. Cuencas'sin información hidrológica. Este es el caso más crítico
y el más común, el cual puede resolverse mediante un análisis re-
gional.
159
Conviene recordar aquí que las técnicas probabilísticas se aplican e
la información hidrológica sólo después que ésta ha sido sometida a ur
tratamiento previo en que se atiende: el relleno de la información, el
análisis de consistencia y la extensión del registro. Para este
jo previo resulta fundamental la referencia 8.
1988 97 116 91 50 28 19 12 11 6 5 2 12
1988 116
n = 30 ]8
162
TABLA 10.1 POSICIONES DE TRAZADO
N = 10
m= 1 m=5
Probabilidad de
Nombre Año excedencia (P) P T P T
m 0.10 10 0.500 2
California 1923 N
1930 2m - 1 0.05 20 0.450 2.2
Hazen 2N
m 0.091 0.455 2.2
Weibull 1939 N+ 1
11
Yn
u = x- - a x a
n
a
a =
n
a
x
Reemplazando valores:
y = ax au
a a Yn
= -!!.x n (X'- a
a a X a
X X n
a a
= -a nx - -
a
n-
X + Yn
-
x x
a a
¿x = -nx- + y - y-
a
x a
x n
y - Y
x = x+ n a
a X
(10.2)
n
--
JI. '- X + Ka (10.3)
x
164
Con la ecuación 10.3 es posible hallar los caudales con largos
períodos de recurrencia (avenida centenaria, avenida milena-
ria, avenida diezmilenaria). Esta ecuación es la ecuación de
una línea recta en papel probabilístico de Gumbel (página 168).
Precisamente la manera de comprobar que el modelo de Gumbel es
el apropiado para el problema en estudio consiste en graficar
la recta y plotear los puntos de la muestra; deberá cumplirse
que todos los puntos caen alineados cerca de la recta.
Ejemplo 10.1
Hallar el valor de K para una probabilidad de 1% si se trabaja
con una muestra de 40 años.
P = 0.01
De la ecuación 10.1: y = 4.600
Para n = 40 Yn = 0.54
CI
n
= 1,14
De modo que:
y - yn
K = CI
= 3.55
n
1 100,
Como T = TY = se indica K100 = 3.55
Ejemplo 10.2
Usando el registro de caudales que se acompaña, encontrar las
crecidas de los 10,100 Y 1,000 años, utilizando la distribu-
ción de Gumbel. Los caudales están en pie 3 /seg.
Año Qmáx Año Qmáx Año Qmáx
1969 36,600 1976 53,100 1983 49,600
1970 69,900 1977 58,800 1984 58,600
1971 99,000 1978 64,100 1985 39,700
1972 76,200 1979 77,800 1986 38,200
1973 62,600 1980 71 ,200 1987 103,000
1974 44,200 1981 59,600 1988 4Z,900
1975 49,200 1982 55,100
165
m Q T m Q T
1 103,000 21.0 11 58,600 1.91
2 99,000 10.5 12 55,100 1. 75
3 77 ,800 7.0 13 53,100 1.62
4 76,200 5.25 14 49,600 1.50
5 71 ,200 4.2 15 49,200 1.40
6 69,900 3.5 16 47,900 1.31
7 64,100 3.0 17 44,200 1.24"
8 62,600 2.63 18 39,700 1.17
9 59,600 2.33 19 38,200 1.11
10 58,800 2.1' n = 20 36,600 1.05
Para n = 20 Yn = 0,52
on = 1.06
Para T = 10 T = 100 T = 1,000
P = 0.1 P = 0.01 P = 0.001
Y = 2.25 Y = 4.60 Y ::;; 6.91
10 = 1.63 100 = J.85
K K K = 6.03
1000
Además: -x = 60,720
°x = 18,233.65
De modo que:
50
I
:60
I
70 lO 90 95 .6
P"O.".ILIO"O
97 9.[-.m '9
.IOOJ
t9.5 99,7 99,8 99.51
.J .: 11111111 11 11 1 1111111
-2.0 -,.& -'.0 _0.5 o 0.5 1.0 1.5 1.0 2·5 3.0 3.5 4.0 4.5 5,0 5.5 -6,0 4;.5 7.0
VA.,A.LE l' DA
RE OU
PAPEL DE PRO B A B I L I DA D DV GUMBEL
--
log x = 1og x + K cr 1og x (l0.4 )
169
El método con s i s te en p1otea r 2 puntos (Q, T) y traza r 1a lí-
nea recta que los une, en papel de 'probabilidad log-normal. Se
prolonga esta recta para encontrar el caudal con período de re
torno grande. -
2.33 12 25 40
10 8 18 38 90
25 12 31 75 105
50 15 39 90 110
100 48 100 115
10.3.6 Probabilidad de
Frente a un caso concreto el ingeniero debe decidir el período
de retorno de su creci da de di seña. Para e 11 o debe preci sar
en primer término la vida útil de la obra; luego se preguntará
sobre la probabilidad de ocurrencia de crecidas durante ese pe
ríodo. Escogiendo un porcentaje adecuado de riesgo
rá el período de retorno de la crecida de diseño y procederá a
calcular la magnitud de ésta con la ley de Gumbel o la distri-
171
bución 10g-Pearson tipo 111; o mejor aun, con las dos a fin de
tomar el promedio.
El período de retorno se ha venido manejando como un valor me-
dio. Es decir, el suceso son período de retorno de T años se-
rá excedido una vez cada T años en promedio; pero también será
excedido en más de la mitad de las veces en una serie de T
años.
Aclaremos. Supongamos un evento de T = 10 años; la probabili-
dad de no excedencia en un año cualquiera es (1-0.1) = 0.9 Y
la probabilidad de no excedencia durante un período de 10 años
es 0.9 10 = 0.35. O 10 que es 10 mismo, la probabilidad de ex-
cedencia de un suceso con un período de retorno de 10 años en
un período de 10 años es 65%.
La fórmula a usar es pues:
(1 _ p)n = probabilidad de no excedencia
1 n
(1 - ¡) = probabilidad de no excedencia (10.5)
172
Toma de decisión
El estudio de las avenidas de largo período de retorno tiene
una aplicación importante en el funcionamiento de los vertede-
ros de demasías de las presas de embalse. Por esta razón el
breve comentario que sigue se refiere a estas obras.
La selección de la magnitud de la avenida de proyecto es algo
sumamente complejo y que debe ser encarado por ingenieros
rimentados en estudios de avenidas. Básicamente se toman en
consideración los efectos derivados de un eventual fracaso de
la obra. La determinación de los daños probables no debe limi
tarse a las condiciones que existen en la época de la
ción. Deben, por el contrario, evaluarse los futuros aprove-
chamientos en la llanura de inundación aguas abajo, como vi-
viendas, granjas, caminos, puentes y otros.
Las presas que forman grandes vasos, construidos en los ríos
principales, con un elevado potencial de escurrimiento, sin
da pueden considerarse en la categoría de muy peligrosas. En
estos casos deberá asumirse un criterio conservador, tomando
como base que no se podría tolerar una falla debido a la
ble pérdida de vidas y debido a los daños potenciales que pue-
den alcanzar las proporciones de una catástrofe. En cambio,
las pequeñas presas construidas en corrientes aisladas en
áreas rurales, donde la falla ni amenaza la vida humana ni pro
duciría daños mayores, se pueden considerar en la categoría de
poco riesgo. En este caso, el criterio de diseño se puede es-
en una base mucho menos conservadora.
En ios Estados Unidos, los vertederos de demasías de las gran-
des presas, cuando existe un gran peligro de vidas humanas,
son diseñados para la creciente más desfavorable posible lla-
mada creciente máxima probable (CMP). La experiencia demues-
tra que la OMP es mayor que aquellas que pueden ser obtenidas
utilizando el análisis de frecuencia. A la CMP no se le puede
asignar una probabilidad realista ni es susceptible de un estu
dio económico; el debe estar basado en la decisión de
dar la másima protección. La es muy grande y casl siempre
superior a las posibilidades de control; su consideración en
el estudio sólo sirve para eliminar la posibilidad adicional-
de una falla repentina en la presa; aparte de la protección de
vidas humanas puede ser útil, por ejemplo, colocar las plantas
de agua potable, así como otros servicios públicos, por encima
del nivel de la CMP.
La obtención de la CMP es un proceso largo, laborioso y costo-
so. Se determina a partir de la precipitación máxima probable
(PMP). La PMP se puede obtener a partir del análisis de pre-
cuencia de las precipitaciones máximas. A partir de la PMP
puede obtenerse la CMP mediante técnicas hidrológicas consis-
tentes en el uso de modelos lluvia-escorrentía e hidrogramas
unitarios. Se requiere mucha experiencia. Las condiciones ex
tremas serán condiciones antecedentes muy húmedas, hidrograma
unitario con pico muy pronunciado y embalse lleno al presentar
se la tormenta. -
173
10.3.7 Metodo de Gradex
Los métodos basados en distribuciones teóricas para el estudio
de frecuencia de crecidas (Gumbel y log-Pearson tipo 111), se
basan en la utilizaci6n de registros de caudales pico.
bargo, los datos de caudales son por lo general cortos y rela-
tivamente escasos. Guillot y Duband elaboraron una metodolo-
gía que permite evaluar los gastos máximos de crecidas utili-
zando registros cortos de caudales conjuntamente con registros
largos de lluvias. El método consiste en extrapolar la curva.
de frecuencia de gastos observados, a partir de cierto valor,
mediante una recta paralela a la de precipitaciones máximas
diarias.
Guillot y Duband desarrollaron el método del Grádex apoyándose
en las siguientes hipótesis:
l. Las precipitaciones máximas diarias siguen la ley de Gum-
bel o valores extremos tipo l.
2. La lluvi2 media diaria sobre una cuenca tiene por gridex
el promedio de los grádex de las subcuencas.
Se entiende por gridex a la pendiente de la recta precipi
taci6n-probabilid'ad en un papel probabilístico de Gumbel-::-
Se demues.tra que es igual a la razón 0x/on'
3. La retención, a medida que aumenta la lluvia, tiende hacia
un límite superior, funci6n de las características de la
cuenca y de la humedad del suelo. A medida que se alcanza
un equilibrio entre el volumen de agua que llueve y el es-
currido todo incremento de precipitación tiende a producir
un aumento igual de volumen escurrido.
El método del grádex es aplicable a cuencas de hasta el c;rden
de los 1,000 km 2 de superficie. En su aplicación' se procede
analizando los gastos observados hasta definir la crecida dece
nal y luego se extrapola la recta con el grádex medio de las
precipitaciones sobre la cuenca.
De la ecuación JO.2:
y - y
x = x + n
°n
= y
O¡
A
10. Uti 1izando las curvas detenni nadas en 8 y 9 se puede obte:-
ner una curva de frecuencia para cualquier sección dentro
de la región en estudio.
cuenca húmeda
cuenco topogrófica
Curva de Duración
En el apartado 6.5 fue descrita esta curva y la forma cómo se cohstru"
ye. Ahora vamos a ver de qué manera se obtiene la curva de duración a
partir de registros cortos.
Una curva de duración será más representativa de las condiciones de es
currimi ento en una cuenca mi entras mayor sea el período de tiempo que
abarquen los datos. Si se cuenta con un regi stro corto, se puede dib.H,.
jar una curva de duración representativa de un perfodo más largo me-
diante el siguiente método recomendade por el U.S. Geologica1 Survey.
1. Se dibuja curva de duración para la estación en estudio (A)
con los registros que ella posee. .
2. Se elige otra estación B de régimen fluviométrico similar re
gistro más extenso. Se construye una curva de duración de esta es
tación B para un perfodo de registro concurrente con el de la esti
ción A. -
3. A partir de las curvas dibujadas en 1 y 2 se determinan los cauda"
les para varias probabilidades de excedencia. Estos puntos selle
van a un gráfico y se determina una curva de. relación entre ambas
es tac iones.
4. Se construye luego la curva de duración de B para topo el período
que tiene registros.
5. En base a las curvas determinadas en 3 y 4 se construye la curva
de duración de A representativa del período de larga duración. Pa
ra ello, se entra a la curva determinada en 4 con una probabi1idaa
de excedencia y se obtiene un caudal para B. Con ese caudal se en
179
tra a la curva determinada en 3 y se tiene el caudal para A con la
probabilidad de excedencia dada. Repitiendo este procedimiento
ra varias probabilidades se puede trazar la nueva curva de dura-
ción para la estación A.
p(y/x} =
180
Xb Xm
¡·SUbCOnjUnto y -1
Evento sobre uno magnitud Xm
Subconjunto x ·1
Evento sobre uno magnitud Xb j total de n eventos
..
\
Espacio muestral; total de N eventos
..1
P(y) = P(y/x)
181
Curvas intensidad-duración-frecuencia
En el apartado 2.7 se describió la forma de obtener las curvas intensi
dad-duración-frecuencia. En el manejo de la fórmula racional
do 7.3) se vió la aplicación de estas curvas para obtener el valor de
i que entra en dicha fórmula.
Lamentablemente, los registros pluviográficos son escasos, de modo que
el procedimiento descrito para obtener las curvas pocas veces eS apli-
cable en la práctica. Esta situación ha llevado a algunos
res a buscar procedimientos que permitan superar tal deficiencia, apro
vechando al máximo la información realmente existente. Frederich BelT
(1969) publicó un trabajo en el cual generalizaba las curvas intensi
dad-duración-frecuencia, a partir de datos recogidos principalmente eñ
Estados Unidos. El argumento físico en que se apoyó es el hecho de
que las lluvias extremas de menos de dos horas de duración se deben a
tormentas de tipo convectivo, las cuales poseen características simila
res en totlas las regiones del mundo.
La expresión matemática propuesta por Be11 es la siguiente:
t duración en minutos
T período de retorno en años
T
Pt precipitación caída en t minutos con período de retorno T
años.
La fórmula es aplicable a lluvias de menos de dos horas de,duración y
con períodos de retorno comprendidos entre 2 y 100 años. Estudios he-
chos en diferentes partes del mundo han a valores sensible-
mente iguales para todos los lugares. Se observa que se requiere cono
cer la precipitación de una hora de duración y 10 años de período de
retorno.
En realidad no siempre se cuenta con información de lluvias de una ho-
ra de duración. Pero Espí1dora, obtuvo en Chile que la relación entre
la lluvia máxima diaria y la lluvia de una hora es más 9 menos constan
te e igual a 4.04. Esto hace posible obtener la lluvia PiS que entra
en la fórmula, a partir de las lluvias máximas diarias cuyos registros
son más frecuentes.
El procedimiento completo entonces viene a ser el siguiente. A partir
de los registros de precipitaciones máximas diarias obtener, mediante
un análisis de frecuencia, la magnitud de la lluvia con período de re-
torno 10 aAos. Usando el coeficiente de Espí1dora tener y luego
aplicar la fórmula de Be11. Por último, calcular, a partir de las
nitudes encontradas de lluvia, las intensidades correspondientes a fin
de poder construir las curvas intensidad-duraci6n-frecuencia.
Mapas de 11uvia-duración-frecuencia
En vez de intensidades se grafican otras veces las cantidades de 11u-
via 1 siempre para diferentes duraciones y frecuencias, lográndose con-
fecciona,r mapas regionales e incluso para todo el país. Estados Uni-
dos, por ejemplo, tiene mapas de lluvia para duraciones entre 30 minu-
182
tos y 24 horas, para períodos de retorno entre 1 y 100 años. Estos
pas se elaboran a partir de precipitaciones en parte observadas y en
parte estimadas a través de correlaciones múltiples usando parámetros
como la pendiente, la elevación, las distancias, etc.
Todos los mapas de 11uvia-duración-frecuencia muestran datos puntuales
de 11uvia,pero se supone que estos valores son aplicables para áreas
de hasta 25 km 2 • La precipitación media en una área debe ser menor
que la precipitación máxima puntual en el área.
""
o 1\" 1"- 1"183 diOs
8 t's; 0,2
6
" r-... 12C díos
'\. "-
.. 4 ........... - 0,1
'\.
""
e l'
i'-. c;
GI
.........
.
..
"V
2
\'\ 6Od(os 2
"U
E
o \, '\ "\
f'
1' ......
r----
-{),05 411
Q.
t'--,. "\ " 30 dios
E
o....
'O O,8I
Q.
"V
-
__ o
o 20 dlol - 0,02 .g
o, t::. :l
8
o, 4
'" "\ lIios
- 0,01
l'
.1"-
l i d ío
"",
-
r--- -
o, 2
PE 110 p Rt OA lNo IN 1': I os
O,1 I
1,01 1,1 1,2 1,5 2 3 4 5 6 8 10 20 30 40
X:t
+----------t
Fig. a Fig. b
a 2
t
::z: Var(x t ), t = 1,. 2, ...
Modelos
Un modelo matemático representando un proceso estocástico se
llama modelo estocástico o modelo de serie de tiempo. Tiene
una cierta forma matemática o estructura y un juego de paráme-
tros.
Un modelo de serie de tiempo simple podría estar representado
por una simple función de distribución de probabilidad.
f (x; e)
185
donde es una serie independiente con promedio cero y varian
za (1 €t es la serie dependiente y es el parámetro
modelo.
Si €t de 10.8 fuera representada por el modelo dependiente
10.9 entonces Xt sería también un modelo dependiente. En este
caso los parámetros del modelo de Xt serían u, a y
Desde que los parámetros de los modelos recién considerados son
constantes los modelos son estacionarios, representando series
de tiempo estacionarias o procesos estocásticos estacionarios.
Los modelos no estacionarios resultan cuando dichos parámetros
varían con el tiempo.
Z = lt - u
Entonces el proceso
--
Zt = 1 Zt_l + <P2 Zt_2 + ... + <PpZt_p + a t (l0.11)
(10.12)
(10.13)
5) Modelos no estacionarios
Se trata de una generalización de los modelos previamente
mencionados. Se trata de representar series cuyo comporta
miento no es estacionario y en particular que no varían en
torno a un valor medio, pero que sin embargo tienen un
cierto tipo de comportamiento homogéneo.
Un proceso de esta naturaleza puede ser representado por-
que sus diferencias de orden d conforman un proceso esta-
cionario tipo ARMA, es decir si
Wt = vd Zt
188
donde vd representa las diferencias de orden d del proceso
Zt:
•
W
y = <1>1 W _ + .•. + q¡pWt_p + a
t 1 t '- el at _1 - ... - eqa t _q
(10.14)
Estimación de
Hasta aquí se ha considerado solamente la función teórica de
autocorre1ación que describe un proceso estocástico. En la
práctica se tiene una serie de tiempo finita Zl' ... , Zn, con
N observaciones, de la cual sólo se puede obtener una estima-
ción de la función.
Se han sugerido diversas estimaciones de la función e y cada su
gerencia tiene características y propiedades particulares. Se
ha concluido que la estimación más satisfactoria del valor de
la función de autocorre1ación para un desplazamiento K es:
189
r = (10.17)
1 N
Z =N L Zt
t=l
Ejemplo 10.3
Si la recesión del flujo en el río es de la fonlla
-Kt
Qt = Qo
Qo caudal al inicio del año cronológico
K constante de recesión
Entonces la serie anual dependiente de Zi seguirá el
modelo autorregresivo ,de primer orden, AR(l):
Z.1 - Z= (Z.1 - 1 - Z) + E·
1
Z promedio de Zi
coeficiente de autorregresión o coeficiente auto-
rregres ivo
E·
1
componente estocástica independiente.
Se demuestra que;
= e-K
Sin embargo, qué t.an bueno resulta este modelo AR(l), depende
de qué tan buena es 1a supos i ci ón de recesi ón exponenci al- en- .. ef.
río.
191
Ejemp lo 10.4
Justificación física de los modelos autorregresivos de prome-
dios móviles para simulación de caudales.
precipitación
Xt
evo poroción 1 O",o+b.,,¡;:1
b Xt r
escorrenti'o superficial
(I-a-b) XtcdXt
- ---
-- - -,..---
aguo 8ubterraneo -
c S t-I
o (10.20)
192
CAPITULO 11 ACERCA DE LA HIDROLOGIA EN EL PERU
11.1 Introducción
De acuerdo a los objeti vos previ amente seña 1ados se han descri too en el
texto en forma ordenada los fundamentos y los métodos hidrológicos or-
dinarios. Acerca de estos últimos, se han descrito los métodos
les aplicables"a cuencas con información y los métodos particulares
aplicables a cuencas con escasa o ningena información, que es el caso
de" gran parte de nuestro territorio.
Corresponde tratar en este último capítulo, de modo especifico, los
pectos relevantes de la Hidrología en el Perú. En torno a esto es ne-
cesario decir que el principal trabajo realizado aquí en el país es el
que se hizo al empezar la presente década gracias al aporte del Gobier.
no del Perú y del Ministerio de Asuntos Exteriores de Italia. El tra-
bajo se enmarcó en el Convenio de Cooperación Técnica para el Estudio
de la Hidrología del Perú, suscrito por el Instituto Italo-Latino
ricano (lILA), el Servicio Nacional de Meteorología e Hidrología
(SENAMHI) y la Universidad NaCional de Ingenierfa (UNI), en el . ámbito
de sus respectivas competencias institucionales.
Se asignó como finalidad del Estudio el proporcionar a los ingenieros
elementos necesarios para evaluar:
1. Los recursos hidrológicos disponibles en las diferentes zonas del
país e identificar las posibilidades que existen para su mejor
aprovechamiento.
2. Las máximas avenidas que pueden verificarse a lo largo de los cur-
sos de agua que atraviesan las diferentes zonas del país e identi-
ficar las posibilidades para reducir los daños que ellas pueden
provocar, por medio de adecuadas capacidades de embalses.
Las metodologías adoptadas para elaborar los datos disponibles y las
conclusiones deducidas se encuentran contenidas en la publicación ESTU
DIO DE LA HIDROLOGIA DEL PERU, disponible en Senarnhi desde Diciembre
de 1982.
Es intención del autor concluir el presente texto con una descripción
sucinta de dicha publicación, a fin de facilitar su uso en la solución
de los problemas de interés practico que se presentan en e1 Perú en el
desarrollo de los Proyectos Hidráulicos.
193
indican los pasos a seguir en la solución de problemas especificas qUE
se presentan en la prictica.
El Volumen 1 está dedicado a la Pluviometria y a la evaluación del
aflujo meteóri cóque púede verifi carse con una determi nada probabil i-
dad sobre la cuenca a una sección genérica del cauce.
El Volumen II está dedicado a la Hidrometría y a la evaluación de las
escorrentías que pueden verificarse con determinada probabilidad en
tal sección (Parte II:A) y a la evaluación de la capacidad de embalse
necesaria para regular las escorrentías con finalidades prefijadas de
utilización (Parte II:B).
El Volumen 111 está destinado a la evaluación de los caudales máximos
y de las escorrentías que pueden ocurrir en una sección genérica en
ocación de eventos de máxima intensidad con una probabilidad asignada,
Se cita, por considerarlo de interés, el siguiente párrafo del Prólogo
de la publicación que nos ocupa:
"En base a los datos disponibles ha sido posible lograr los objetivos
previstos para la mayoría de las zonas, mientras que, para las restan-
tes, la información obtenida de·10s datos a disposición resultó total-
mente inadecuada. Al respecto, las conclusiones del Estudio permite
completar el diseño de la red nacional de estaciones hidrometeoro1ógi-
cas ".
El Plano lOA muestra la subdivisión del territorio patrio en zonas p1u
viornétricas. En total son 1, 2, ... 10 zonas pluviométricas; 1, 2, ..
16 subzonas pluviométricas y A, B, .,. L áreas anómalas.
El Plano 28 muestra los 8 grupos de cuencas hidrológicamente similares
en que igualmente ha sido subdividido el territorio.
El Anexo A contiene:
- Cuadro 1, en se indican las estaciones pluviométricas utiliza
das en eles tudio (códi go, cuadrante, nombre, n, cp, Y, Oro, Oc). -
número de años de registro
latitud, en grados
altitud, en m.s.n.m.
distancia a la línea de costa, en km
distancia al borde occidental de la cordillera (línea de de
limitación hacia occidente de las zonas pluviométricas 3 y
4), en km.
Y{h) = m Mn
e ji
con e .. =100 ,el coeficiente de utilización para el uso i en
J1 i el mes j
con Ci' el diagrama de demanda para el uso i
con Pi' la deficiencia porcentual que puede ser aceptada para
el us o i
Pi
con E; mín = (1 - mr) Ei , el valor mínimo al cual puede tolerarse
que el suministro total anual E1 descien-
da en los años de escorrentía mln;ma,
196
Se admite que los ingenieros que laboran en cada sector SUmlnlS-
tran los valores Ei , c i ' Pi' Ei mín' Con esta información y si-
guiendo la metodología señalada en el Anexo se obtienen los diagra
mas de demanda (c) y (c l
). -
Las figuras 11:8-33 a 11:8-40 del Volumen 11, representan los dia-
gramas de disponibilidad deducidos para los diferentes grupos de
cuencas hidrológicamente similares en que se ha dividido el terri-
torio (Plano 2B).
El procedimiento consiste en precisar la ubicación de la cuenca en
el Plano 2B y ver qué diagrama de dispÓnibi1idad le corresponde de
los contenidos en las figuras 11:8-33 a 11:8-40 del Volumen 11.
se asume = 0(0)
in
y se obtiene: Dnl = M(D n) - 1.96 o(D n)
Dn2 = M(D n) + 1.96 o(Dn )
4. Cálculo de T,r D
O, ,
T,r ... valor mlnimo al cual puede descender O una vez cada 'T
años con un riesgo r.
Se calculan los valores T,r D, en función de M(D), utilizando la in
formación contenida en los cuadros Il:8-1 al del Volumen Ir:
200
APENDICE 1 INUNDACIONES, HUAYCOS y SEQUIAS
Se reproduce aquí el excelente trabajo del Ing. Rosenqo Chávez Díaz, pu-
blicado con el mismo titular en la revista El Ingeniero Civil, números 26 (Se-
tiembre-Octubre 1983) y 29 (Marzo-Abril 1984).
El Ingeniero Civil es una revista de información profesional editada en
Lima por el Instituto de Publicaciones de Ingeniería Civil (PUBLICIVIL), insti
tución que alberga en su seno a profesionales peruanos de reconocido prestigio.
El trabajo fue publicado con ocasión de los desastres naturales ocurridos
en ese entonces en el territorio peruano y se reproduce aquí, co'n anuencia de
su autor, porque hay que tener presente, como él mismo señala, que lIeste tipo
de fenómenos han ocurrido muchas veces en tiempos pasados y volverán a ocurrir
muchas veces más en el futuro ll
•
PRIMERA PARTE
Aunque en los años anteriores inmediatos han ocurrido huaycos, algunas lluvias
costeras y sequías parciales, ha sido durante este año 1983, que se ha llegado
al clímax de daños, resultantes de dos fenómenos climáticos contrapuestos.
Inundaciones y huaycos en el Norte y Centro del país y severas sequías en el
Sur. Como secuela de estas violencias climáticas queda, en adelante, una in-
mensa tarea de reconstrucción de ciudades, caminos, canales, etc. y, 10 que es
mu'cho más grave, la recuperación agrícola y pecuaria en todas partes de Costa
y Sierra.
Este tipo de fenómenos han ocurrido muchas veces en tiempos pasados y volverán
a ocurrir muchas veces más en el futuro. Antes de hoy la resonancia de los da
ños fue menguada por la falta de comunicación oportuna, por la lejanía y por=-
"que el menor desarrollo de nuestro país determinó un mucho menor volumen de
daño.
El crecimiento de nuestra población que se hacina hoy en cada vez más grandes
ciudades y en un número mayor de centros poblados, construyendo dentro de la
mayor improvisación y sin planificación razonable, ha significado, finalmente,
exponer sectores importantes de las ciudades a torrenteras e inundaciones.
En cuanto a caminos y vías de comunicación la engañosa seguridad derivada de
varios años consecutivos de años hidrológicos pobres, con es correntías relati-
vamente modestas y el olvido de los esporádicos años con violencias pluviales
y fluviales, dio paso a soluciones donde priman aspectos inmediatos de orden
económico y/o pob1aciona1 dejando a un lado precauciones de seguridad
mente innecesarias. El resultado: caminos deshechos y centenares de millones
de soles de pérdidas en productos, combustibles, vehículos, etc, sin contar
con las trágicas pérdidas de vidas.
En el área agrícola si bien el crecimiento físico ha sido pequeño (en algunos
valles nulo y en otros con retroceso) en cambio las mejores tecnologías
das a la explotación de los suelos ha aumentado la intensidad de ésta, expre=-
sándose en mayores volúmenes de productos de mucho mayor valor económico y so-
cial. De allí que la inundación de áreas de cultivo, las pérdidas de suelos
por erosión y la destrucción de las estructuras de aducción de agua no s610
201
significan fuertes pérdidas de producto bruto vendible sino que exigen inmedia
tas e importantes inversiones, tanto para reponer las estructuras destruidas,
total o parcialmente, como para nuevas obras necesarias. A esto hay que agre-
gar las inversiones, en dinero y tiempo, que significa recuperar los standars
y niveles de producción previos al desastre.
Con ocasión de los fuertes daños ocurridos este año, que ha conmocionado al
país y consecuente con la resonancia que, en estos casos, le dan los medios de
comunicación masivos, han surgido una serie de cuestionamientos y, de las sim-
ples preguntas de años anteriores: ¿Cómo es posible que suceda ésto? ¿Por qué
no se hizo ésto, aquéllo o lo de más allá? ¿Qué ha hecho, hace o hará el Go-
bierno? se ha pasado a la acusación dictada por la desesperación. En estos ar
tículos se trata de explicar, sintetizando, el cómo y por qué de los fenómenos
ocurridos y a partir de su conocimiento tratar de buscar y sugerir soluciones
razonables.
202
1) Restringe fuertemente, al enfriar el agua, las tasas de evaporación
desde el Pacífico, reduciendo proporcionalmente los posibles volúme-
nes de lluvia potencial; y
2) Al enfríar a las capas de aire en contacto con el espejo de agua, re-
duce considerablemente el poder de ascender de esa masa de aire, cuya
pequeña humedad suele condensarse a poca altura formando las nieblas
costeras, que alimentan a las llamadas "lomas" o "pastizales efíme-
ros".
I
!
500
100 A
,
-- " 1\
v
/' ..... A I /" r'\.
o
:=
N
V
o
U V \
1\
! !
Gráfico No. 1
203
1.2 La Corriente del Niño. Su Incidencia
La corriente caliente, llamada del Niño, es un brazo terminal de la de Ca
1ifornia (originada a su vez, por la Kuroshivo, del Japón) que corre de
Norte a Sur, es decir, contraria a la de Humboldt, que se va perdiendo
lentamente al entremezclar sus aguas con las frías de esta última. Cuan-
do estos encuentros ocurren enfrente del Ecuador las condiciones climáti-
cas, de la costa peruana, son "normales", es decir, se producen las llu-
vias orográficas, hay aridez en la costa y las descargas de los ríos ocu-
rren dentro de los límites usuales.
Desde que las corrientes son elementos dinámicos, no existe un área espe-
cífica donde, persistentemente, se efectúe ese encuentro sino que puede
ocurrir en cualquier punto de una amplia extensión de mar. Esas diferen-
tes ubicaciones sucesivas ocasionan "fluctuaciones" climáticas a 10 largo
de los años; en la costa del Pacífico.
Cuando el área de convergencia se desplaza hacia el Sur, más de 10 usual,
se producen precipitaciones de importancia en ciertas áreas de la costa
peruana. Esto se debe a que, permaneciendo las lluvias orográficas más o
menos las mismas y en su área, acusan su presencia y se suman a ellas las
lluvias convectivas procedentes del Océano. En efecto, consecuentes con
la mayor temperatura del agua del mar, las tasas de evaporación son más
altas y las masas de aire, ahora húmedas y calientes, pueden ascender fá-
cilmente. La suma de ambas precipitacio,nes determina mayores y más per-
sistentes masas de lluvia, abarcando además mayor área de cuenca o toda
ella. Cuanto mayores sean los desplazamientos hacia el Sur y la tempera-
tura del agua del Pacífico, mayores en intensidad y frecuencia serán las
lluvias y abarcarán una extensión mayor de costa.
Resumiendo, la acción de la Corriente del Niño, se pOdría definir como li
beradora de las restricciones climáticas actuales de la Costa
con 10 cual ésta recupera su clima típicamente tropical que debería ser
el normal.
1.3 Periodicidad y
Como lo prueflan las· expell'jencias habidas en años pasados, las grandes
cipitaciones y sus secuelas se repiten según algún período, de recurrencia.
Los geólogos y arqueólogos han podido establecer, a partir del estudio de
suelos y sub-suelos, var-ias ocurrencias anteriores fijando, aproximadamen
te, la época en se produjeron. Entre los registrados se encuentran
algunos como el de 1891 que originó el camóio de trayectoria del Río Piu
ra; en 1925 se prOdujeron también desastres importantes (quedó enterrado
como ejemplo, el Muelle y Puerto de Samanco) y hemos llegado a los de
1982 y sobre todo 1983 que no sólo sería el peor de todos sino que abarca
mayor extensión.
La muy sucinta relaci6n hecha se refiere a los sucesos de mayor magnitud,
pues de rangos menores se han produci do muchos: 1931, 1939, 1972, etc.
Esta sucesión de años con problemas ha dado lugar a la búsqueda de una re
lación numérica que permita, eventualmente, predecir la nueva oportunidad
de repetición y se ha hablado de diferentes períodos de recurrencia. In-
fortunadamente el período estadístico es muy corto y no ha sido posible
lograr una relación razonable. Se sabe, con seguridad, que el fenómeno
velverá a ocurrir pero no cuándo ni con qué intensidad. (ver gráfico N2
1) .
204
Cuando se trata de fenómenos naturales que se repiten con intensidad va-
riable, cada vez, y en especial en Hidrología se recurre a una figura con
venciona1 que establece una relación fenómeno/tiempo. Se dice, por ejem:-
p10, que una intensidad 1 corresponde a un período de recurrencia de 20
años, que uno de intensidad 2 a uno de 50 años y que otro de 3 a 100 años
y así sucesivamente. Esta. es solamente una manera de racionalizar esa re
lación fenómeno/tiempo y constituye, en realidad, una escala de magnitu:-
des que sirve, en i"ngeniería, para establecer la relación economía/proyec
to para cubrir una magnitud dada de riesgo. De allí que una obra de cos:-
to n será diseñada para un riesgo de magnitud 1 que, se supone, ocurrirá
una-vez en 25 años. Para cubrir un riesgo 2 es probable que el costo de
la obra alcance 2 (o más) n y podrá servir, eficientemente, para cubrir
un riesgo que se supone será igual a uno de 50 años de recurrencia y así
sucesivamente. El problema grave estriba en que, desde que estos fenóme-
nos están ocurriendo desde hace muchos millones de años, nunca se sabe
cuándo se inició el período correspondiente a un fenómeno del tipo 1,000
años, por ejemplo, y éste puede tener lugar en cualquier momento y arrasa
rá con las obras que, dentro de las ineludibles exigencias económicas to-=-
madas en cuenta, han sido diseñadas para fenómenos de intensidad corres-
pondiente a 50 ó 100 años de recurrencia, de acuerdo a la escala conven-
cional.
El lector se dará cuenta, con mayor claridad, de las dificultades, hoy in
salvables, que no permiten establecer una periodicidad definida o siquie:-
ra medianamente aceptable, al examinar las razones que luego se exponen
relacionadas con un elemento netamente dinámico, dependiente de múltiples
factores: las lluvias.
4
... .... -
(2) _4
.-'-
z 3
.,z .... ....
z
... ...e ce
CL
_
ce
CD
......
...>
CI •::i V ti 1ce
1/)
111
.J
i
., .,ce
z
/ e
ce
CI
i ...
1
.. CL
/'" '"
V ce
ce
ce
Z
.J
ce
. - - 111
4 --;7
.""
ZONJI HUME DA ",
" ""
O
10 ",
"
N",
3
ZONA
. /''''T
I \
Z
..,.. 2 / I
.,.,
I
/
v/
.,z / I /
z 1.!5 .------
I
----- --;;0/
I ,/ '"
....
Q
,['1"
SA NTA EULAL A
.J
ce I .:. /
V__
p.--:q.--} NANA
I 0+ /
d-
. --BI (CA NPO DE NARTE)
La Curva Sutton
Las cuencas peruanas de la vertiente del Pacífico presentan, normalmente,
los tres tipos de cuenca, húmeda, árida y semi-árida. Mientras no reci-
ban la influencia del fenómeno del Niño el complejo hidrogeológico
na dentro de límites "normales", con problemas rutinarios y manejables
si siempre. Cuando la influencia generada por dicho fenómeno es excesi-
va, en uno u otro sentido, ocurren ya sea la tropicalización violenta del
área semi-árida o la sequía extrema. Dentro de los extremos esquematiza-
dos hay, desde luego, diferentes grados de intensidad en función de la
magnitud de las causas y el área comprometida varía, aumentando o dismínu
yendo, constantemente.
207
Hasta aquí se ha tratado, casi exclusivamente, del fenómeno del Niño co-
mo el factor tangible más resaltante pero, como en todo proceso complejo,
y éste lo es en grado sumo, intervienen en grado menos conspicuo los anti
ciclones, el vulcanismo (investigadores acuciosos han establecido ciertas
secuencias de evidente paralelismo entre estos fenómenos y efectuado com-
probaciones interesantes) y otros más.
Las cuencas norteñas solamente presentan zonas húmedas y semi-áridas pues
las lluvias, esporádicamente. llegan a más de 250 mm/año a pocos metros
sobre el mar. Confonne se avanza hacia el Sur (y la incidencia de la co-
rriente fría de Humboldt es cada vez más nítida y única) el límite prome-
dio de la isoyeta 250 se sitúa cada vez más alto llegando, en Moquegua y
Tacna a ocurrir, normalmente, a altitudes superiores a los 4,000 msnm. En
este último caso se define claramente la presencia de los tres tipos y,
mientras más al Sur, las áridas son cada vez de mayor extensión. En la
región de Lima la altura estaría por cerca de los 2,000 msnm (ver gráfico
3).
1oI0ll0COI:H.
208
La cuenca húmeda "normal" del Río Rímac es del orden de 1,245 km 2 , es de-
cir, corresponde al 36.7% del área total y en ella se generan masas anua-
les de escorrentía del orden de 900 millones de m3 /año, o sea con descar-
gas mayores del orden de 85 m3 /s promedio. (Hay más de 55 años de esta-
dísticas). Las zonas semi-áridas con 925 km 2 (27.3% del total) y la arl-
da con 1,220 km 2 (36%) ocupan el resto de los 3,890 km 2 del total de cuen
ca hasta el mar. -
Cuando el fenómeno del Niño hace sentir su acción la isoyeta 250 descien-
de hasta el límite inferior de la semi-árida y se produce es correntía en
un área 175% mayor de lo normal. Puede rendir, entonces,casi 1,500 MMC/
año y producir caudales hasta de 500 m3 /s. (No se dispone, todavía,de ci-
fras estadísticas de 1983).
Si bien estas cifras, promedio, son explícitas y espectaculares, no refle
jan la situación de los caudales generados en quebradas laterales
mente sin agua y que, de un momento a otro, pueden descargar decenas y
ci entos de m3 /s.
Como estos caudales extremos se producen luego de varios años (el 70% der
tiempo ocurren escorrentías normales) arrastran consigo grandes cantida-
des de sedimentos generados y depositados en los años intermedios, por in
temperismo en las áreas yermas y desnudas de las inclinadas laderas cord,
lleranas: son los huaycos, llocllas y yapanas.
Notas:
Como la extensión de este artículo es ya excesiva, en proxlmo nUmero se
tratará de la relación entre estos procesos naturales de la fisiografía,
ecología y geología sugerir planteamientos de soluciones de los per-
tinentes problemas.
CARLOS W. SUTTON, ingeniero peruano (de origen americano) que inició en
el Perú la era de la hidrología aplicada a proyectos de irrigación. Ge-
nial e incansable partió de cero y ha dejado una rica herencia en esta
tecnología específica. Durante casi 40 años estudió la hidrología en
nuestro país y estableció las bases de todos los proyectos peruanos de
irrigación.
SEGUNDA PARTE
En la primera parte de este artículo se trató del fenómeno del Niño, brevemen-
te, sus secuelas: inundaciones, huaycos, yapanas, etc. No es necesario abun-
dar sobre su peligrosidad y cuantía de daños, actuales y su proyección al
ro, cuya gravedad es manifiesta y por ello se hará, en este artículo, sucintos
comentarios relacionados con soluciones para estos problemas.
Parece oportuno, al tratar de pstos asuntos, exponer algunas reflexiones de im
portancia indudable pl s, si bien se refieren al comportamiento humano
tuadamente peruano, en este caso) su incidencia para llevar a buen término las
soluciones pertinentes, es contundente. En efecto, para lograr resultados es
necesario considerar lo siguiente:
1) Decisiones inmediatas y
Los fenómenos naturales que producen estos problemas exigen soluciones
2lJ\J
que requieren muy largos períodos para dar resultados tangibles y satis-
factoriDs, por consiguiente, es imprescindible tomar decisiones inmedia-
tas para lograr resultados en tiempo oportuno. Ojalá que, en esta ,oca-
sión,surta efecto el factor catalizador que significan las calamidades su
fridas en este año y las, al parecer probables, del próximo.
Al respecto conviene resaltar que la magnitud de los daños derivados de
estos fenómenos son, en alguna forma, función del mayor uso del territo-
rio nacional por una población cada vez mayor y cada día con mayores exi-
gencias de todo orden (alimentación, vivienda, comunicación, etc). Los
daños ocurridos en periodos criticos similares anteriores (como ejemplos:
1891 y 1925 en cuanto a huaycos e inundaciones y, en el Sur, la sequia
del periodo 1936/45), no tuvieron ni los alcances ni la resonancia de los
producidos en años recientes (1972 y 1983).
La cifra del Cuadro N2 1 Y el gráfico N2 1 son sumamente reveladoras.
En primer término el promedio nacional ha crecido, en algo más de 100
años, de 2 a 13 Hab/km 2 (6.5 veces) y es probable que el año 2000 haya
llegado a más de 20 Hab/km 2 . Dentro de este contexto el Norte (Tumbes,
Piura, Lambayeque, La Libertad y Cajamarca) ha crecido al mismo ritmo del
promedio nacional (7 veces) llegando a una densidad de 35 Hab/km 2 , deter
minando -con excepción de Lima 140 Hab/km 2 1981 y 20 aumentos- la mayor
ocupación territorial.
Desde que la ocupación del territorio, ya se refiera al crecimiento de
las ciudades, a la mayor exigencia sobre la producción agropecuaria o a
la mayor intensidad dei uso de las vías de comunicación, es proporcional
a la mayor población, resultaría que los daños que en 1940 afligieron a
1 persona, en 1983 han lesionado a 3 personas con intensidad claramente
mayor que la simple relación 3/1. Obsérvese que, de acuerdo a las ten den
cias expresadas geométricamente en el gráfico, al año 2000 es probable
que se haya duplicado la población actual y entonces el daño recaerá so-
bre 6 veces más personas y la relación será mucho mayor que 6/1. Faltan
16 años para el 2000 y las decisiones deben ser tomadas de inmediato para
lograr, a tiempo, realizaciones que eviten los daños futuros o los cir-
cunscriban a limites razonablemente estrechos. No se puede olvidar que
la acción de la naturaleza es continua, persistente y sin pausa y eXlge
del actuar humano 10 mismo para lograr resultados adecuados y consisten-
tes.
30
25
/
1/
U)
20
""z
lo-
15 ,.- c---- ",'
.....- /
--
<t
!: 10
en
<t
% 5 r--------- -- ,
""o 1--
U)
8 --- LIMA
""oz 6 --,.-
- 01 _i co
/
/
/
CENTRO
..J NORTE
..J
i 4 c---.,
- /
/ /
/
SUR
---- /-
f.---- ,- SELVA
.-
2
1
O
- - r- -
......
-- - -- ,
"
ZONA AREA A Ñ O S
(1 ) km 2 1876 1940 1961 19?2 1981 2000 (2)
Norte 113, 5 13 n 28 35 53 6T O
Centro 244, 4 9 12 lb 18 25 6'0
Sur 243, 3 6 7 9 11 16 4'0
Selva 651, 0.15 0.5 1 1.4 2 5 3'0
LIMA 34, 7 24 60 102 140 235 8'0
Aspectos Ecológicos.- Estos, junto con los hidrológicos, son los p'ocesos más
complejos y dificiles. Nuestras cuencas del Pacífico sufren, en este sentido,
limitaciones muy serias:
213
a) Por razón de la presencia de zonas áridas y semiáridas que limitan tremen-
damente la posibilidad de cubierta vegetal, persistente y adecuada, conse-
cuencia de la falta de abastecimiento de agua;
b) Por razón de la altura, pues por encima de los 3,800 a 4,000 msnm, el cli-
ma, de dureza extremada, reduce a escasos tipos de pastos, pequeños y muy
rústicos, la vegetación capaz de vivir en las serranías. A partir de los
4,500 msnm prácticamente no hay vida vegetal.
Estas condiciones naturales dificultan, grandemente, soluciones sencillas
de simple reforestación, por ejemplo. Hay, siempre, exigencias adicionales de
riego artificial, adecuación de las tierras (los andenes, por ejemplo), etc.
Habría que adicionar, en las alturas, los efectos de las heladas (ver más ade-
lante) que afectan muy seriamente contra la supervivencia de las plantas.
Desde este punto de vista la prospección de cochas, bofedales, pequeños
embalses, canales de riego, etc. complementados por andenerías, surcado de la-
deras, etc. tenderán a resolver los problemas para lograr la cubierta vegetal
satisfactoria que, además de la seguridad de acondicionamiento de la cuenca,
representa una fuente de recursos económicos de origen agropecuario.
.
......
Ó
40+-
300
750
..,
.
e
...•
.!
'"e 'E
500
el)
ftl IIJ
IIJ 200 Z
-J IIJ
Q Z
250 ::>
::> -J
100 ..
o
o CO\Jd o.'e 5
L.om,!, >
100
O
5 10 15 20 25
AREAS DE CUENCAS Km 2
ri vaci ón y ri ego junto' con 1as grandes prospecci ones: grandes embalses,
les de fuerza, etc. todas las cuales, aunque podrían ser financiadas y
das individualmente, deben ser planteadas como piezas interdependientes de un
'conjunto en los aspectos concernientes a operación y mantenimiento.
Como conclusión general, deducidas de las.consideraciones anteriores, se
llega a que hay solucicnes reales, efectivas y económicas para los problemas
de esta índole. Que no es razonable ni práctico proponer planteamientos com-
plejos de corto plazo (o cortísimo, como es usual reclamar), llenos de compli-
caciones técnicas, pero que sí es procedente y lógico planteamientos integra-
les, compuestos por muchos proyectos, generalmente poco relevantes y más bien
simples, que actúen como piezas interdependientes e interoperantes de un gran
todo. Que es agente muy importante el edafológico que, a través de seres vi-
vos, activos y persistentes en su acción {la cubierta vegetal) se opone exito-
samente a esa fuerza dinámica, de peligrosidad fluctuante e imprevisible punto
·de aplicación, como es el agua en grandes masas conjugadas con factores fisio-
gráficos y geológicos.
Conviene recalcar que el ingrediente principal, esencial, es la decisión
humana de hacer y, sobre todo, persistir en la acción para lograr las metas
apropiadas con indesmayable voluntad.
DESIERTOS Y SEQUIAS
1.- Consideraciones generales
En términos panorámicos el fenómeno de la sequía, que debería ser llama-
da, con mucha mayor propiedad, deficiencia de ocurre a todo lo lar-
go y ancho del mundo y se presenta en períodos variables en cuanto se re-
fiere a duración y frecuencia. La "sequíall se expresa en forma permanen-
te en los desiertos, es decir, áreas sometidas indefinidamente a una agu-
da escasez de agua. Otra forma corresponde a deficiencias de dicho recu!:,
so durante períodos, los períodos críticos, de mayor o menor duración. En
Perú se presentan ambos casos: de los 150,000 km 2 de Costa hay 80,000 km 2
de los desiertos más áridos del mundo y 63,000 km 2 de cuasi-gesiertos y,
215
apenas, 6,500 km 2 de áreas cultivadas, es decir sólo el 4.5% y, en la Sie
rra, unos 100,000 km 2 (25% de la misma) están sometidas a esporádicas
quías, algunos de los cuales pueden abarcar hasta 6 ó 7 años consecutivos
pobres en lluvias. Esto significa que el 20% del territorio nacional es-
tá afectado por el fenómeno de la sequía.
En el ámbito mundial no menos de 1/3 del área continental (casi 50'000,000
km 2 ) corresponde a desiertos y áreas en proceso de desertificación y, lo
que es más grave, se van extendiendo "inexorable y rápidamente", según 10
establecen las conclusiones de la Primera Conferencia Internacional sobre
Desertificación (Nairobi, 1977). Estos desiertos afectan a no menos de
700'000,000 de seres humanos, de los cuales unos 100 millones son agricu1
tores, quienes han perdido sus tierras de cultivo. Se admite que unos
75 millones adicionales se verán afectados, por este fenómeno, a corto
plazo.
En las áreas desérticas se producen lluvias esporádicas, a veces de gran
intensidad, pero cuyo promedio anual no permite la vida persistente vege-
tal, fluctúa entre 100 mm/año (Sahara, Africa), 150 mm/año (Gobi, Asia) y
alrededor de 200 IlIm/año (Australia). 'Los desiertos americanos son muchí-
simos más áridos: 70 mm/año en América del Norte (Valle de la Muerte,USA)
y unos 20 mm/año en cualquiera de los peruanos y algo menos en Tarapacá y
Atacama.
Como se acaba de indicar, aun en los desiertos hay una cierta intensidad
de precipitación; conviene establecer entonces un criterio a partir del
cual se puede calificar una zona como desértica (sequía grave y permanen-
te) o transitoriamente deficiente de recursos de agua. En principio ese
criterio se basa, esencialmente, en si hayo no suficiente agua, a lo lar
go del tiempo, para que haya vegetación que permita la vida humana y, des
de luego, la vida animal. Ambos beneficiarios, la fauna y la flora, usan
el agua ininterrumpidamente a lo largo de toda su existencia para sus vi-
tales funciones fisiológicas y cuando no la hay suficiente, la planta o
animal, muere. Los animales, incluyendo los humanos, consumen cantidades
reducidas de agua. El hombre necesita alrededor de 1/2 litro/día (toda
la humanidad consumiría apenas 4 millones de m3 /día; el Perú entero unos
36,000 m3 /día) pero necesita alimentos diariamente yesos alimentos (in-
cluso carnes) significan un consumo muchas veces mayor de agua pues las
plantas para producir cosechas requieren entre 4 y 6,000 m3 /Ha/año hasta
12 a 20,000 m3 /Ha/año como consumo neto.
Estas cifras establecen, con claridad, que es el consumo vegetal el deter
minante para establecer la condición de sequía de una área. Si se puede
lograr cosechas se habrá superado esa condición.
De allí que, para el caso de cultivos de secano (riego natural, con llu-
via únicamente) es necesario que las lluvias sean del orden de magnitud
indicado, sobre 400 mm/año y para áreas con riego artificial es necesario
contar con fuentes de agua que permitan conducir a las chacras entre 5 y
1.25 veces los volúmenes (4 a 5 veces en caso riego tradicional: por
surcos, por inundación, etc, y 1.25 veces en el de riegos sofisticados:
aspersión, goteo, etc).
La Sequía en la Costa
En la Costa las 650,000 Has cultivadas están exigiendo probablemente alre
dedor de 10,000 millones de m3 /año, agua que recibe de unos 133,000 km T
de cuencas donde 11 ueve al rededor de 400 'mm/año. Convi ene hacer notar
216
que la masa de lluvia promedio es del orden de 33 millones de m3 /año y al
ser utilizados solamente 10,000 MMC/año se está desperdiciando el 70% de
los recursos promedio anuales.
Como se ve, los apenas 6,500 km 2 cultivados de la Costa se surten de los
ríos que descienden de la Cordillera Occidental que suman unos 33,000
MMC/año, de los cuales solamente se utilizan 10,000 MMC/año y sin embar-
go,'tada cierto lapso de tiempo ocurren sequías o deficiencias de abaste-
cimiento de agua. La razón de estas ocasionales sequías corresponde con
las variaciones de intensidad de la precipitación en la Cordillera y la
correspondiente ubicación de la isoyeta 250 mm/año, como se indicó en la
primera parte de esta exposición. Esto expl ica, por ejemplo, que el Río
Ica descargue algún año 535 MMC (1925) y en otro apenas 77 MMC (1947). En
el primer año hubo excedentes de agua que fueron al mar y en el otro al-
canzó apenas a cubrir el 15% de los requerimientos del valle y determinó,
por consiguiente, pérdidas enormes a la agricultura local.
La solución, en estos casos, estriba en grandes proyectos que teniendo un
alcance plurianual permitan reservar los excedentes de agua de años con
abundancia de ellas para completar los riegos en los años pobres.
Son proyectos más o menos sofisticados que requieren inversiones de impor
tancia variable que pueden incluir algunas (o todas) las estructuras mayo
res como: túneles, presas, trasvases, etc. En forma genérica, se puedeñ
sintetizar indicando que si bien requieren inversiones de mayor magnitud
son, también por la cantidad, calidad y seguridad de la producción, los
de más alto rendimiento económico. Las previsibles deficiencias de abas-
tecimiento de agua quedan, si no totalmente eliminadas reducidas a un mí-
nimo. El relieve nacional asegura, además, el beneficio adicional de am-
plia disponibilidad de energía hidroeléctrica.
3.- Sequías Serranas
Para expl i ca r mejor los conceptos de "sequ ía en su fase de defi ci enci a
11
., ....
'\
\ o"
. .oL
j /__ 'V
,'" ,
I
,
v."
I
r-;:,
\
,,
',... ..,./'\
'" 100 or
I v'''' J" ,,'
''¡
"
,.
//'.j
'v . í
\i CABANILL AS
01
>
\1 (PUNO
fCO
3
.J • 1 --..... )180'7.70 NIVELES LAGO ITITICACA
. / r--. ./
t>-/ ..............
• -1
e -2
./
-; -5 .......... ./"'"
,
'--"'"
E -4
-e
..,
ID o .,.
." o on o 1ft
,..
In o
!
Grófico No. 3
218
ma (5 km 3 /año, promedio). Además de éste tan serio se han presentado
ciclos cortos, incluso de menor intensidad relativa, en períodos no
menores de 3 años, de tal manera que, generalizando, en no menos de
45 años de cada 100 'habrá algún grado de sequía en el Altiplano y de
ellos en no menos de 30 (o sea el 30%, del tiempo) corresponden a gra-
dos muy severos de deficiencias.
b) Las Secuelas.- La incidencia de la sequía, sobre la ecología puneña,
es tremenda. No solamente mueren las plantas de escasa raíz (sobre
todo pastos -sobrevive, miserablemente, sólo el ichu-) sino que, la
falta de nubosidad permite alta insolación la cual, combinada con las
bajas temperaturas de esas (SObre los 3,800 msnm) determinan alta fre
cuencia de severas heladas que completan de arruinar lo poco que
biese quedado.
219
APENDICE 2 PROGRAMA EN PASCAL PARA COMPLETAR INFORMACION PLUVIOMETRICA POR
EL METODO DE .LA·"RECTA DE REGRESION.
program hidro17;
const max=100 ;maxl=l1 ;
type vector=array[I .. maxJ of real;
matrix=array[I .. max,l .. maxlJ of real;
var media,desv:vector; .
a:matriz;
r:array[2. ;maxlJ of real;
m,n,rnl,i,j,jmax:integer;
rmax : rea 1; .
begin
l*leamos la matriz*)
wrltel'orden de la matriz '};readlnÚD,n);
for i:=l to m do
begin
write('fila J,i,';ngrese ',n,'datos: ');
for j:=l to n do read(a[i,jJ};
readln
end;
(*inicialicemos en cero los vectores media,desv,r y la variablem1*)
m1:=0;
for j:=l to n do
begin
medi;- [j J:=O;
desv [j J:=O;
if j>l then r[j J:=O
end;
(*calculemos la media*)
for i :=1 to mdo
if a [i,1 J>=O then
begin
for j:=l to n do
medi a [j J:=medi a [j J+a [i ,j J;
mI :=rn1+1
end;
for j: =1 to n do
medi a [j J: =medi a [j J/oü;
l*calculemos la desviación estándar*}
for i: =1 to -m do
if a [i ,1 ]>=0 then
for j: =1 to n do
desy [j J: =desy [j]+ sqr( a [i ,j]- medi a [j J) ;
for j:=l to n do
begin
desy [j ]: =desy [j JI (ml-1) ;
desv [j J: =sqrt(desv [j J)
end; .
(*calculemos los valores r*)
for i:=l to m do
if a[i,I]>=O then
221
for j:=2 to n do
r[j]:=r[j]+(a[i,1]-media[1])*(a[i,j]- media[j]);
for j:=2 to n do
r[j]:=r[j]/(ml-l)/desy[l]/desy[j];
(*obtengamos el r máximo*)
nnax:"-2;jmax:=0
for j: =2 to n do
if r[j ]>rmax then
begin
nnax :=r[j];
jmax:=j
end;
(*hagamos el relleno de datos*)
for i: =1 to m do
i f a [i ,1 ]<0 then
a[i,l]:=media[l]+r[jmax]*desv[l]/desv[jmax]*(a[i,jmax]
-media[jmax]);
(*escribamos los resultados*)
for i:=1 to m do
wri te 1n Ca [i ,1 ] : 6 :1 ) ;
writel n;
wri te ln(medi a [j]:6:1) ;
write 1n (des y[j ] : 6 :1) ;
writeln;
writeln('estación =',jiuax:2);
wri teln(' r máximo"" ,nnax:3: 2)
end
222
APLICACION DE LOS COMPUTADORES ELECTRONICOS EN HIDROLOGIA
(tomado de la referencia 2)
223
HIDROLOGIA
FE DE ERRATAS
El po I II 0Ie1lbto
mIIMta la'
ó._
la ............; .. de 11» ¡)ji : .... " .......,,_
doI ado hodool6yio» Y11» .. 4IvdoA 0tCIir\IIIw"......,;6n .... pt'OtIM-
cM!.
como 1m.
,.art.nt. hl'rTetnitotlla de lfl\bafO .... al .rtMilII y de .... probIe.
moo.... lulo <IedictOO un CQII"CIItoIo ... 1'odrcIogi. Eatltdl'tiee
do!Ide""" 00I'I áeno _ • •..,¡jO loo n'WOOeIos
pmt>eI_IIO: )1 V el! ron...- ' l O ... 'n ,' ...... j