HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA.

JANETT FLORES

CAUDAL

1.- INTRODUCCIÓN
Hasta ahora hemos visto las características fisiográficas de la cuenca como volumen de control, las
variables meteorológicas y la precipitación como variable de entrada de agua al sistema.

Entradas I(t) Operador  Salidas O(t)

En este capítulo veremos como medir la variable de salida, el caudal, en el punto de drenaje. Si
podemos medir en forma precisa la precipitación y el caudal, será posible entonces determinar los
procesos internos () que inciden en la relación precipitación – escorrentía, tales como intercepción,
infiltración y otras pérdidas hidrológicas.

Estaciones pluviométricas

Estación fluviométrica

2.- EL CAUDAL
El Caudal o gasto de una corriente se define como el volumen de agua que pasa por una sección dada
de un canal en un tiempo dado y tiene dimensiones de volumen sobre tiempo.

volumen Q  m 
3

Q = Caudal = 
tiempo t  seg 

El caudal puede expresarse como el producto de un área por una velocidad así:
Q=A*V
donde: A : área de la sección (m2) y V: velocidad (m/s).

1
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Una vez que la precipitación alcanza la superficie del suelo, se infiltra hasta que las capas superiores
del mismo se saturan. Posteriormente, se comienzan a llenar las depresiones del terreno, y al mismo
tiempo, el agua comienza a escurrir sobre la superficie hasta llegar a un cauce para finalmente ser
drenada hasta la salida de la cuenca.
Este flujo sobre el terreno se denomina escurrimiento superficial, una vez que alcanza el cauce se
denomina escorrentía.
Si se mide el caudal que pasa por una determinada sección transversal de un río y se grafican los
valores obtenidos contra el tiempo, se obtendría una gráfica que se denomina hidrograma.
Q (m3/s)

Tiempo

El área bajo el hidrograma es el volumen total escurrido y la lámina de escurrimiento será ese
volumen escurrido entre el área de la cuenca.
Un milímetro de escurrimiento es la cantidad de agua necesaria para cubrir el área de drenaje hasta
un milímetro de altura; es una unidad útil para comparar el caudal con la precipitación que lo produjo, se
conoce como "Lámina de Escurrimiento". El mm. es unidad de volumen sólo cuando va asociada con
un área de drenaje específica.

Volumen escurrido = Q * t

donde : Q : caudal medio (m3/s)

t : tiempo (s)

V agua = A drenaje * Lámina de Escurrimiento

V agua
L.Esc.= ---------------- (mm)
A drenaje

2
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

El caudal no puede medirse de manera continua, pues se requiere cierto tiempo de medición para
evaluar un caudal instantáneo Q, entonces para conocer la relación Q vs. t debe recurrirse a otra
variable, en este caso el nivel que si se puede medir en forma continua con un flotador . A partir de la
relación H vs. t se obtiene la gráfica Q vs. t usando una función de transferencia Q vs. H llamada curva
de caudales ó Curva Nivel – Caudal que debe construirse a partir de los datos de campo.

H vs t H vs Q Q vs t

3. MEDICIÓN DEL NIVEL

El término Nivel se refiere a la elevación de la superficie de agua en una sección del río, a partir de un
datum de referencia, que puede ser ó la cota del nivel medio del mar ó el punto para el cual el caudal Q
es igual a cero.

3.1. Miras
La Mira es un instrumento utilizado para medir en forma puntual la elevación de la superficie de
agua. Puede colocarse en forma vertical ó inclinada un ángulo con respecto a la vertical. Se sujeta
firmemente a una estructura permanente, tal como la columna de un puente. Está graduada con marcas
grabadas en la mira, normalmente con precisión de un centímetro. Una parte siempre debe estar
sumergida en el agua, en caso de que sea necesario se coloca una serie de miras para abarcar todo el
rango de alturas, todas referidas al mismo datum.

3
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

3.2. Medidor de niveles máximos.

El medidor de niveles máximos está diseñado para obtener lecturas de los caudales picos de las
crecientes. Consiste de una regla graduada de ancho y largo suficiente para entrar en un tubo
galvanizado de 5 cm de diámetro. El tubo galvanizado es sellado en sus extremos con unas tapas y se
hacen varios orificios de 6 mm. El tubo se monta verticalmente en el canal, referido al datum de la
sección. La mira se inserta por la parte superior y se le añaden virutas de corcho. Al ocurrir la creciente
máxima, el corcho se adhiere a la parte húmeda, y al bajar el agua, el corcho queda como señal del nivel
que alcanzó el agua. El operador remueve la tapa superior, lee el nivel y limpia la regla para su próximo
uso.

AGARRADER
A

ORIFICIOS

AGARRADER
A

Medidor de nivel máximo

3.3. Limnígrafo ó fluviógrafo

Son instrumentos que miden en forma continua los niveles de agua en una sección determinada del
río. En general se usan de dos tipos: (a) de flotador (b) sensor de nivel de agua.
En el caso del uso del flotador, debe construirse a la orilla del río un pozo entre 0.30 a 1,00 m de
diámetro, con una profundidad por lo menos igual al mínimo nivel de agua. Se conecta al río con
tuberías galvanizadas horizontales de 60 cm de diámetro. El lomo de la tubería normalmente se asume
como el punto de caudal mínimo. Si el río es profundo se instalan varios tubos a diferentes
profundidades. El aparato se instala en una caseta ubicada arriba del pozo y que lo protege de la
intemperie. El agua sube por el pozo hasta igualar el nivel de agua en el río (vasos comunicantes) , pero

4
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

reduce las fluctuaciones propias del flujo turbulento del río. El registro puede se hace sobre papel
colocado sobre un tambor que gira con un mecanismo de reloj, o un mecanismo que abre un hueco en
un papel cada cierto intervalo predefinido de tiempo. Así se obtiene la curva H vs. T

4.- MEDICIÓN DEL CAUDAL.

Los caudales no pueden ser medidos en forma directa, por lo que se mide otra variable (nivel de
agua) y se establece la relación funcional entre caudal y nivel para cada sección, evaluada con una curva
denominada curva de caudales o curva nivel - caudal. La construcción de la curva de caudales requiere
de un número representativo de pares de valores nivel - caudal. Como esta curva de caudales puede
cambiar con el tiempo por factores tales como la erosión o la sedimentación del material en el lecho del
río, se hace necesario su chequeo continuo para así poder determinar en forma segura los caudales a
partir de los niveles registradores.
- Aforos.
Los métodos de medición directa del caudal se denominan generalmente, métodos de aforo.
Existen diversos métodos para determinar el caudal de una corriente de agua, cada uno aplicable a
diversas condiciones que varían según el tamaño de la corriente o según la precisión con que se
requieran los valores obtenidos. Entre estos métodos, los más conocidos son:
- Método volumétrico
- Método pendiente - área
- Método de los aforadores de profundidad crítica
- Métodos químicos
- Método Sección - Velocidad

- Método Sección - Velocidad:

El método se basa en la ecuación de continuidad:
Qt =  Qi = ( Ai * Vi )
Donde: Q es el caudal en m3/s, V es la velocidad del flujo en una sección transversal del canal en m/s
y A es el área de la sección en m2 .
La sección del canal se divide en varias subsecciones. El número de espaciamientos dependen del
ancho del canal (Tabla N 1), normalmente el mínimo es de diez. No necesariamente deben ser
igualmente espaciadas. Algunos criterios pueden ser: (a) Los espaciamientos no deben ser mayores de

5
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

1/15 a 1/20 del ancho del cauce. (b) La diferencia de las velocidades entre dos secciones no debe
exceder el 20% (c) El caudal en cada sección debe ser menor del 10% del caudal total.
El área se determina mediante sondeos en diferentes verticales, los cuales permiten calcular las áreas
parciales mediante la construcción de una serie de trapecios; la suma de estas áreas es el área total de la
sección. En cada vertical se mide la velocidad del agua mediante distintos métodos de los cuales el más
práctico y económico es el método con molinete hidráulico ó correntímetro.
El molinete hidráulico es un aparato provisto de una rueda con aspas o copas, la cual gira por el
choque del agua contra ellas permitiendo conocer la velocidad de la misma en el lugar en que está
colocado expresada en revoluciones por segundo N.
V = a + bN
Donde:
N = rev/s
a , b = parámetros de calibración del aparato.
V = velocidad del agua en m/s
Como la velocidad del agua varía de un punto a otro en una misma sección transversal de una
corriente, es indispensable efectuar medidas en diferentes sitios y profundidades para conocer el valor de
la velocidad media en cada trapecio, la cual al multiplicarse por el área parcial correspondiente,
proporciona el caudal que pasa por ese trapecio. Sumando todos los gastos parciales se obtiene el
caudal total que pasa por la sección.
La distribución vertical de la velocidad en un río no es uniforme a lo largo del ancho y de la
profundidad del cauce. La velocidad es mayor en la parte mas profunda del cauce y casi cero en los
límites del flujo. La mayor velocidad ocurre justo debajo de la superficie, en ñla parte mas profunda de
un cauce recto. LA turbulencia o curvaturas pueden ubicar la máxima velocidad en otro lugar del cauce.

3.0

2.5
0
2.0
2.2
.5 1.5
0 2.0
2.2
Distribución de velocidades en un canal.5no uniforme
0

6
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Desde el punto de vista práctico se han definido algunas reglas para ubicar los puntos de medición de
las velocidades en la vertical con el fin de calcular la velocidad media en la vertical que representará la
velocidad media en la sección considerada.

Distribución vertical de velocidades en un cauce.

Tabla N 1. Espaciamientos deseables para diferentes tipos de cauces.

Características del canal Número mínimo de Ancho de espaciamientos
Observaciones - vertical (m)
Canal menor de 15 m de 15 1 a 1,5
ancho, o con cambios
bruscos de sección
Canal con ancho de cauce 15 3a6
entre 15 y 90 m
Canal con ancho de cauce 15 6 a 15
entre 90 y 180 m
Canal con ancho de cauce 25 -
mayor de 180 m
a) La velocidad promedio en la vertical v vertical para cauces poco profundos que no excedan los 3
metros, se mide a 0,6 veces la profundidad del flujo desde la superficie, v 0,6. Se mide un solo punto.
b) La velocidad promedio en la vertical v vertical para cauces moderadamente profundos se calcula
como
v 0,20 + v 0,80
v vertical = ---------------------
2

7
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

donde v 0,20 es la velocidad medida a 0,2 de la profundidad del flujo desde la superficie y v 0,80 es la
velocidad medida a 0,8 de la profundidad del flujo desde la superficie.

El caudal en cada sección se calcula como _

Qi = Ai * v vertical

ei  ei  1
Ai  * Pi
2

5.- CURVA DE GASTO ó CURVA NIVEL -CAUDAL

Cada vez que se realiza un aforo se calcula el valor del caudal instantáneo que pasa por la sección del
río considerada. A su vez se mide el valor del nivel relacionado con ese caudal. Si se realizan varios
aforos al año se obtendrán varios pares de valores Q, H. Los pares de valores nivel y caudal se grafican
en papel aritmético o logarítmico produciendo las denominadas curvas de gasto o relación nivel-caudal.
Estas curvas deben estar definidas al menos por 10 puntos y es recomendable que la desviación de
cualquiera de ellos de la curva, no sea en ningún momento mayor del 5% para obtener una definición
satisfactoria. Para la porción más baja de la curva de caudales, se recomienda graficar esta porción en
escala mayor en la parte izquierda superior del papel.

0,19 0,22
2,25 0,63
FIG. 4,2
5,74 CURVA0,99
DE CAUDALES
9,53 1,25
3,5
14,5 1,52
3
26,03 1,98
NIVEL (mts)

2,5
40,07 2,41
2
53,84 2,65
1,5
64,53 3,03
1
0,5
0
0 20 40 60 80
3
CAUDAL (m /seg)

8
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

- Extensión de la Curva de Gasto

Frecuentemente sucede que para las crecidas más altas, no se conoce la relación nivel-caudal, y se
hace necesario extrapolar las curvas existentes. Dado que la curva de nivel-caudal posee una forma
exponencial, para proceder a la extrapolación debe transformarse en una línea recta. Para ello existen
varios métodos, entre ellos:
- Método Logarítmico.
- Método de Stevens o A D
- Método Pendiente - Nivel - Caudal.
El método más usado es el Logarítmico.

- Método logarítmico
Este método consiste en transformar a coordenadas logarítmicas los valores de Caudal y Nivel, lo
cual debe dar por resultado una línea recta, haciendo la aproximación de que la curva se comporta como
una potencial. Sin embargo esta condición no se cumple por la variabilidad de la pendiente del lecho, lo
que da por resultado que la curva sea una combinación de muchas curvas potenciales.
En general y como resultado de la investigación de una gran cantidad de curvas de gastos, la relación
entre el caudal y el nivel puede ser representada aproximadamente por la ecuación:

Q = b * (h ± z)a

Tomando logaritmos a la ecuación Q = b * (h - z)a se tendrá :

Log Q = a * Log(h ± z) + Log b
Esta es la ecuación de una recta referida a los ejes coordenados logQ y log (h ± z) siendo "a" la
pendiente y Log(b) el intercepto. Normalmente se plotea log(h± z) vs. logQ, por lo tanto el valor de "a"
sería el inverso de la pendiente geométrica y el valor de "b" es igual al caudal cuando Log(h± z) = 1,
obteniéndose entonces la ecuación de la curva, la cual puede expresarse en la forma primitiva:

Q = b * (h± z)a.

Siendo esta ecuación de forma potencial, donde Q es el caudal en m3/s, h es la altura de mira en m,
b, z y a, son constantes. Si se plotea Q y (h-z) en un papel bilogarítmico los puntos deberían formar una
recta. Si la ecuación se cumpliera estrictamente para todas las alturas del cauce, se tendría que el valor
de z representaría la altura de la escala para caudal igual a cero, entonces el valor de z seria tal que Q =
0 y (h-z) = 0 sin embargo, para caudales mínimos la curva se aparta de la forma de una línea recta,
existiendo algún caudal muy pequeño cuya altura no es posible medir en los aforos y por ello se fija el

9
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

valor más bajo (z), a fin de poder obtener un punto tan bajo que sirva de guía para trazar la curva, en su
parte inferior. Por consiguiente es necesario considerar a "z" como el valor que debe usarse para que al
plotear (h-z) vs Q (papel bilogarítmico) la relación sea una línea recta.

Ejemplo.
Se observa una creciente excepcional en una corriente. El caudal pico ocurrió la noche anterior. No
fue posible efectuar una medición de caudal ( aforo ). El máximo nivel observado fue de 4,42 mts. Las
mediciones de nivel - caudal realizadas en mediciones anteriores están dadas en la tabla 4.1. Asuma que
la relación nivel - caudal es de la forma:

Q = b * ( h ± z )a
Se pide obtener un estimado del caudal de la creciente.

Tabla N2. TABLA DE DATOS

3
NIVEL ( m ) Q (m /s)
0,22 0,19
0,63 2,25
0,99 5,74
1,25 9,53
1,52 14,50
1,98 26,03
2,41 40,07
2,76 53,84
3,03 64,53
3,92 116,39

Solución:
Se grafica la curva nivel - caudal usando los datos de la tabla, donde se observa que no hay el
valor de Q = 0

10
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

0,19 0,22
2,25 0,63 FIG. 4.4
5,74 0,99
CURVA NIVEL - CAUDAL
9,53 1,25
414,5 1,52
26,03 1,98
40,07 2,41
3
53,84 2,76
NIVEL (mts)

64,53 3,03
2

0
0 10 20 30 40 50 60 70

CAUDAL (m/s)

Se estima un valor de z = 0 y la graficamos en escala log.log ( los datos de las columnas (2) y (5) de
la tabla 4.2 ).
0,19 0,22
2,25 0,63 F I G. 4.5
5,74 0,99
9,53 1,25
10
14,5 1,52
26,03 1,98
40,07 2,41
53,84 2,76
NIVEL

64,53 3,03
1
11,39 3,92

0,1
0,1 1 10 100

CA UDA L

Se observa en la figura 4.5 que la línea tiene una ligera curvatura hacia arriba en el límite inferior.
Esto se debe a que el valor de la constante z no es cero.
Examinando la figura 4.4, sugiere que el nivel correspondiente a un Q = 0 puede encontrarse entre z
= 0,01 y z = 0,10

11
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Se estima un valor de z = 0,090 y se obtiene la columna (3) restando de la (1) 0,090 (valor fijado
para z ) a cada valor de nivel (tabla 4.2). Estos valores se vuelven a graficar (columna (3) y (5)) en la
escala log.log.(fig.4.6).
0,19 0,18
2,25 0,54 F I G. 4 .6
5,74
10 0,9
9,53 1,16
14,5 1,43
26,03 1,89
40,07 2,32
53,84 2,67
64,53 2,94
NIVEL

1
116,39 3,83

0,1
0,1 1 10 100 1000
CA UDA L

Se observa que esta vez la curvatura de la línea es invertida.

La gráfica de (h± z ) Vs. Q en escala log.log debe ser una línea recta; y en ambos casos tienen
curvaturas en sentidos opuestos. Esto significa que entre estas dos líneas se encuentra la que define una
recta, y el valor de z está entre 0 y 0,090.
Luego se estimó z = 0,04 y se obtiene la columna 4, se grafica (fig.4.7) y se obtiene una línea recta.

12
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

0,19 0,19
2,25 0,59 F I G. 4. 7
5,74 0,95
9,53
10 1,21
14,5 1,48
26,03 1,94
40,07 2,37
53,84 2,72
NIVEL

64,53 2,99
1
116,39 3,88

0,1
0,1 1 10 100 1000

CA UDA L

Se determinan los valores de b y a. El valor de a es el inverso de la pendiente geométrica.

log Q  log Q
2 1
a=
log(h2  z)  log(h1  z)

log 116,39  log 0,19

a=  2,09
log 3,88  log 0,18
a = 2,09
El valor de " b " es igual al caudal cuando el valor de Log(h -z) = 1
Por sustitución:

Q = b * ( h ± z )a

Q = b * (h ± 0,04)2,09
Seleccionando un par de valores de la curva de caudales, y sustituyendo en la ecuación, obtenemos el
valor de b.
Par de valores seleccionado: (2,25; 0,63)

2,25 = b * (0,63 - 0,04)2,09 => b= 6,77

a ) Nivel = 4,42 m

Q = 148,34 m3/s

13
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

TABLA 4.2.
NIVEL (m) NIVEL-0,0 NIVEL-0,090 NIVEL-0,04 Q
(1) (2) (3) (4) (5)
0,22 0,22 0,13 0,18 0,19
0,63 0,63 0,54 0,59 2,25
0,99 0,99 0,90 0,95 5,74
1,25 1,25 1,16 1,21 9,53
1,52 1,52 1,43 1,48 14,50
1,98 1,98 1,89 1,94 26,03
2,41 2,41 2,32 2,37 40,07
2,76 2,76 2,67 2,72 53,84
3,03 3, 03 2,94 2,99 64,53
3,92 3,92 3,83 3,88 116,39

6.- Ejemplo de obtención de un hidrograma a partir de la banda del fluviógrafo

Sobre la cuenca del río Palmar cae una tormenta que produce una crecida. Esta es registrada por el
fluviógrafo ubicado en la estación Las Múcuras. Con el fluviograma resultante (fig.4.8) y la curva nivel
vs. caudal de la sección (fig.4.9), obtenga la lámina de escurrimiento directo, si el área de la cuenca es
de 390 Km2.

0 0,4
3 0,6 FIG.4.8
8 0,8 CURVA NIVEL vs. CAUDAL
15,5
2 1
28 1,2
48 1,4
1,6
70 1,6
)t 110 1,8
1,2
(m 135 2
L
E
I 0,8
V
N
0,4

0
0 20 40 60 80 100 120 140
CAUDAL (m3/s)

14
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Solución:
Se obtienen los datos del fluviograma para construir la curva H vs. t. La escala horizontal representa
el tiempo en horas, dividido cada 12 horas, y la escala vertical representa el nivel en metros, divididos
cada 0,05 metros.
Si elegimos un intervalo de lectura igual a seis horas (tomando en cuenta el pico) tenemos:

12 18 24 6 9 12 18 24 6 12 18 24
t (hr)
0,71 0,70 0,70 1,20 2,01 1,61 1,31 1,20 1,12 1,06 1,01 0,97
H (m)

12 0,71
18 0,7 FIG. 4.9
24 0,7 CURVA NIVEL vs: TIEMPO
30 1,2
2,433 2,01
36 1,61
2
42 1,31
NIVEL (mt)

1,648 1,2
54 1,12
1,2
60 1,06
0,866 1,01
72 0,97
0,4

0
6 12 18 24 30 36 42 48 54 60 66 72 78

TIEMPO (hr)

Entrando con valores de t en la curva nivel vs. tiempo obtenemos diferentes valores de nivel, con los
cuales entramos en la curva nivel vs. caudal y obtenemos valores de caudales, con estos datos se forman
pares de valores caudal - tiempo con los cuales se construye la curva (Q vs. t).

15
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Area = V.E.D. = 4228200 m3

V.E.D.
L.E.D. = --------------
A cuenca

L.E.D.= 10,84 mm  11 mm

16
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

7.- PERIODO DE RETORNO (Tr).

Período de retorno es el intervalo medio de años que ocurre entre eventos que igualan o exceden
una magnitud dada.

n 1
Tr =
m

Donde:

m : Número de orden en una lista de mayor a menor de los datos.

n : Número de datos.

PROBABILIDAD DE OCURRENCIA (P) .

Probabilidad de ocurrencia es la probabilidad de que un evento sea igualado o excedido en un

intervalo de tiempo dado y analíticamente se expresa como P = 1/Tr. La probabilidad de no ocurrencia
para ese mismo período, vendrá dada por P' = 1 - P ; sustituyendo tenemos :

1
P' = 1 -
Tr

Probabilidad de ocurrencia de un evento dentro de su período de retorno.

a.- Probabilidad de que un evento X ocurra cualquier año

1
P(x) =
Tr

b.-Probabilidad de que un evento X no ocurra cualquier año.

P` (x) = 1 - P(x)

1
P’(x) = 1 -
Tr

c.-Probabilidad de que un evento ocurra en el año r-nésimo.

r 1 1 r 1 1
( P') *P=(1- ) *
Tr Tr

17
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

d.-Probabilidad de que un evento no ocurra en "n" años consecutivos.

n 1 n
( P )n = ( 1 - P ) = ( 1 - )
Tr

RIESGO..

Es la probabilidad de que un evento de una magnitud dado sea igualado o excedido dentro del
período de diseño.

1 n n : período de diseño
R=1-(1- )
Tr

Ejemplo :

Determinar el riesgo de falla de una obra que tiene una vida útil de 25 años, si se diseña para un
período de retorno de 50 años.

Solución :

En este caso Tr = 50 años y n = 25 años.

Sustituyendo en la ecuación

1 n
R=1-(1- )
Tr

1 25
R=1-(1- )
50

R = 0,40 R = 40 %

El riesgo es de 40 %, es decir, se tiene una probabilidad del 40 % de que la obra falle durante su
vida útil.

18
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

8.- DISTRIBUCION GUMBEL O DE VALORES EXTREMOS TIPO I .

(Función de Distribución Exponencial Doble).

Gumbel basado en que un evento máximo puede definirse como el máximo valor de una serie de
365, si la distribución de los eventos máximos es ilimitada, propuso que la probabilidad de no ocurrencia
de un valor igual o mayor que X se ajusta a una distribución doblemente exponencial.

 ( x   )
e
P'(X<x) = e
 ( x   )
e
P (Xx) = 1 - e
Donde los parámetros  y ß se estiman como:

1,2825
=
S

ß = X - 0,45 * S

Donde:

X : es la media.

S : desviación típica.

Digamos entonces que:

( x-ß ) = Y

1
P (X>x) =
Tr

19
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Por lo cual tenemos que:

1  e y
 1 e
Tr
1 y
1-  e e
Tr

Linealizando:

1
Ln(1  )   e y
Tr

1
Ln( Ln(1  ))   y
Tr

1
 Ln( Ln(1  ))  y
Tr

1
 Ln( Ln(1  ))
Tr  x  


1
 Ln(  Ln(1  ))
x Tr  


donde :

Tr : Período de retorno

x : Caudal para un período de retorno

Ejemplo:

Los caudales máximos anuales registrados en la estación hidrométrica La Alcabala en el río san
Pedro se muestran en la tabla 7.2. Se plantea construir cerca de este sitio un bordo para protección
contra inundaciones. ¿Cuál debe ser el caudal de diseño si se desea que el período de retorno sea de 60
años?

20
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Valores de caudales instantáneos máximos anuales

AÑO X (m3/s)
1954 2230
55 3220
56 2246
57 1804
58 2737
59 2070
60 3682
61 4240
62 2367
63 7061
64 2489
65 2350
66 3706
67 2675
68 6267
69 5971
70 4744
71 6000
72 4060
73 6900
74 5565
75 3130
76 2414
77 1796
78 7430

Solución:

Para 25 años de registro, de la tabla 7.2 se tiene:

S = 1825,91 m3/s

X = 3886,0 m3/s

Por lo tanto las ecuaciones

1,2825
 y ß = X - 0,45*S
S

 = 0,000702 s/ m3

ß = 3064,34 m3/s

21
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Para Tr igual a 60 años tenemos:

1
y =  Ln(  Ln(1  )) = 4,0860
Tr
y
x = ----- + x = 8884,79

Q máx.= 8884,79 m3/seg.

9.- DISTRIBUCION LOG - PEARSON TIPO III .

La distribución de Log-Pearson tipo III plantea que el valor de una variable para su período de

retorno dado se expresa como LogxTr = Logx - KSLogx.

Es una distribución de tres parámetros : Logx , Slogx y Cs

donde :

Logx = media de los logx.

Slogx = desviación típica de los log x.

K = factor que se determina en función de Cs y Tr.

Cs = coeficiente de asimetría

n * (log x  log x ) 3
Cs =
(n  1) * (n  2) * ( Sloxx ) 3

1
Tr =
P

Pasos a seguir en el método Log - Pearson.

1.- Se construye la serie anual de la variable.

2.- Se determinan sus logaritmos.

3.- Se determina la media de los logaritmos.

4.- Se determina la desviación típica de los logaritmos.

22
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

5.- Se determina el coeficiente de asimetría de los logaritmos .

6.- Con Cs y Tr se entra a la tabla 7.4 y se obtiene K.

7.- Se determina el Log(X) para un Tr determinado.

8.- Se obtiene el valor de la variable para Tr determinado.

Ejemplo:

Resolver el ejemplo anterior utilizando la distribución Log Pearson tipo III, para un periodo de
retorno de 50 anos.

Solución:

En la tabla se presentan los pasos 1 y 2

X
AÑO (m3/s) Log X
LogX - LogX
1954 2230 3,3483 -0,1935
55 3220 3,5079 -0,0339
56 2246 3,3514 -0,1904
57 1804 3,0017 -0,5401
58 2737 3,4373 -0,1045
59 2070 3,3160 -0,2258
60 3682 3,5661 0,0243
61 4240 3,6274 0,0856
62 2367 3,3742 -0,1676
63 7061 3,8489 0,3071
64 2489 3,3960 -0,1458
65 2350 3,3711 -0,1707
66 3706 3,5689 0,0271
67 2675 3,4272 -0,1145
68 6267 3,7971 0,2553
69 5971 3,7760 0,2342
70 4744 3,6761 0,1343
71 6000 3,7782 0,2364
72 4060 3,6085 0,0667
73 6900 3,8384 0,2966
74 5565 3,7455 0,2037
75 3130 3,4955 -0,0463
76 2414 3,3827 -0,1591
77 1796 3,2543 -0,2875
78 7430 3,8710 0,3292

 88,5458  -0,1792

23
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

LogX = 3,5418

Slog X = 0,2173

Cs = - 0,0254

Entrando en la tabla 7.4 con el valor de Cs y Tr obtenemos el valor de K.

K = 2,0404.

Log X = LogX + K * Slog X

Log X = 3,5418 + 2,0404 * 0,2173

Log X = 3,5418

X 50 anos = AntLog (3,5418) = 9664,96 m3/s

Debe comprobarse si la distribución teórica se ajusta a los valores observados, para esto se usan
métodos gráficos y pruebas estadísticas.

24
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

FACTOR K PARA LA DISTRIBUCIÓN LOG-PEARSON TIPO III.

Cs Periodo de retorno en años
2 5 10 25 50 100 200 1000
3,0 -0,396 0,420 1,180 2,278 3,152 4,051 4,970 7,250
2,5 -0,360 0,518 1,250 2,262 3,048 3,845 4,652 6,600
2,2 -0,330 0,574 1,284 2,240 2,970 3,705 4,444 6,200
2,0 -0,307 0,609 1,302 2,219 2,912 3,605 4,298 5,910
1,8 -0,282 0,643 1,318 2,193 2,848 3,499 4,147 5,660
1,6 -0,254 0,675 1,329 2,163 2,780 3,388 3,990 5,390
1,4 -0,225 0,705 1,337 2,128 2,706 3,271 3,828 5,110
1,2 -0,195 0,732 1,340 2,087 2,626 3,149 3,661 4,820
1,0 -0,164 0,758 1,340 2,043 2,542 3,022 3,489 4,540
0,9 -0,148 0,769 1,339 2,018 2,498 2,957 3,401 4,395
0,8 -0,132 0,780 1,336 1,998 2,453 2,891 3,312 4,250
0,7 -0,116 0,790 1,333 1,967 2,407 2,824 3,223 4,105
0,6 -0,099 0,800 1,328 1,939 2,359 2,755 3,132 3,960
0,5 -0,083 0,808 1,323 1,910 2,311 2,686 3,041 3,815
0,4 -0,066 0,816 1,317 1,880 2,261 2,615 2,949 3,670
0,3 -0,050 0,824 1,309 1,849 2,211 2,544 2,856 3,525
0,2 -0,033 0,830 1,301 1,818 2,159 2,472 2,763 3,380
0,1 -0,017 0,836 1,292 1,785 2,107 2,400 2,670 3,235
0,0 0,000 0,842 1,282 1,751 2,054 2,326 2,576 3,090
-0,1 0,017 0,836 1,270 1,716 2,000 2,252 2,482 2,950
-0,2 0,033 0,850 1,258 1,680 1,945 2,178 2,388 2,810
-0,3 0,050 0,853 1,245 1,643 1,890 2,104 2,294 2,675
-0,4 0,066 0,855 1,231 1,606 1,834 2,029 2,201 2,540
-0,5 0,083 0,856 1,216 1,567 1,777 1,955 2,108 2,400
-0,6 0,099 0,857 1,200 1,528 1,720 1,880 2,016 2,275
-0,7 0,116 0,857 1,183 1,488 1,663 1,806 1,926 2,150
-0,8 0,132 0,856 1,166 1,448 1,606 1,533 1,837 2,035
-0,9 0,148 0,854 1,147 1,407 1,549 1,660 1,749 1,910
-1,0 0,164 0,852 1,128 1,366 1,492 1,588 1,664 1,800
-1,2 0,195 0,844 1,086 1,282 1,379 1,449 1,501 1,625
-1,4 0,225 0,832 1,041 1,198 1,270 1,318 1,351 1,465
-1,6 0,254 0,817 0,994 1,116 1,166 1,197 1,216 1,280
-1,8 0,282 0,779 0,945 1,035 1,069 1,087 1,097 1,130
-2,0 0,307 0,777 0,895 0,959 0,980 0,990 0,995 1,000
-2,2 0,330 0,752 0,844 0,888 0,900 0,905 0,907 0,910
-2,5 0,360 0,711 0,771 0,793 0,798 0,799 0,800 0,802
-3,0 0,396 0,636 0,660 0,666 0,666 0,667 0,667 0,668

10.- PRUEBAS DE DE AJUSTE.

Las pruebas de ajuste consisten en verificar si los valores observados X y sus respectivos
Períodos de Retorno Tr = ( n+1)/m se ajustan a los valores calculados según la distribución teórica
escogida, Gumbel ó Log-Pearson.

25
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

Para la Distribución de Gumbel se diseñó un papel especial , donde al plotear en el eje y la

variable X y en el eje x la Probabilidad de Ocurrencia o Período de Retorno los valores observados, ellos
se alinearán en caso de ajustarse a la distribución doblemente exponencial. En el caso de Log-Pearson es
mas difícil por ser una distribución de tres parámetros por lo tanto se pueden plotear los valores en una
hoja de probabilidad log-normal sólo en el caso de que el coeficiente de asimetría sea muy parecido o
igual a cero.Por ejemplo, ploteando en papel Gumbel los valores de caudal vs. período de retorno (fig.
7.1) presentados en la tabla 7.2, se observa que existe un ajuste.

Existen además Pruebas estadísticas de ajuste de datos , tales como la Prueba de Chi-Cuadrado y
la Prueba no paramétrica de Smirnov Kolgomorov.

.- Prueba Kolmogorov - Smirnov.

Esta prueba consiste en comparar el máximo valor absoluto de la diferencia D entre la

probabilidad observada Po (xm) y la estimada P(xm).

D = máx | Po(xm) - P(xm) |

Con un valor crítico d que depende del número de datos y el nivel de significación seleccionado
(tabla 7.5), los valores de  que se usan normalmente son de 10; 5 y 1 %. El valor de , en la teoría
estadística, es la probabilidad de rechazar la hipótesis nula . Si D < d, se acepta la hipótesis nula, Ho.

H0 = La función de distribución de probabilidad es P (,,..).

VALOR CRÍTICO DEL ESTADISTICO DE SMIRNOV - KOLGOMOROV, PARA VARIAR

VALORES DE N, Y VALORES CORRIENTES USADOS EN HIDROLOGIA.


N 0,20 0,10 0,05 0,01
5 0,45 0,51 0,56 0,67
10 0,32 0,37 0,41 0,49
15 0,27 0,30 0,34 0,40
20 0,23 0,26 0,29 0,36
25 0,21 0,24 0,27 0,32
30 0,19 0,22 0,24 0,29
35 0,18 0,20 0,23 0,27
40 0,17 0,19 0,21 0,25
45 0,16 0,18 0,20 0,24
50 0,15 0,17 0,19 0,23

26
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

La función de la probabilidad observada se calcula como:

n 1
Tr =
m

Donde:

m : Número de orden del dato xm en una lista de mayor a menor.

n : Número total de datos.

xm : Variable.

y la función teórica en el caso de Gumbel es:

 ( x   )
e
P (Xx) = 1 - e
Para el caso del ejemplo anterior se ha elaborado una tabla.

En la columna 1 se han escrito los caudales máximos anuales registrados, ordenados de mayor a
menor, en la columna 3 se calculan los valores de la probabilidad observada, en la columna 4 se calculan
por el método de Gumbel los valores de la probabilidad estimada teórica y en la columna 5 se muestran
los valores absolutos de las diferencias entre Po(xm) y P(xm).

Se ha encerrado en un rectángulo el valor Dmáx. Según esta prueba se aceptaría la distribución

para un nivel de significancia de  = 0,05, para el cual el valor crítico d es 0,27 con n = 25 (de la tabla).

D = 0,1390

d = 0,27

D < d , es decir, 0,1390 < 0,27 , por lo tanto se acepta la hipótesis, la distribución se ajusta.

27
HIDROLOGÍA PROFA. SUSANA DE BAUTISTA – PROFA. JANETT FLORES

TABLA DE AJUSTE
Xm m P0 P( Xm) P0 -- P( Xm)
7430 1 0,0385 0,0455 0,0070
7061 2 0,0769 0,0586 0,0183
6900 3 0,1154 0,0654 0,0500
6267 4 0,1538 0,1001 0,0537
6000 5 0,1923 0,1194 0,0729
5971 6 0,2308 0,1217 0,1091
5565 7 0,2692 0,1586 0,1106
4744 8 0,3077 0,2646 0,0431
4240 9 0,3462 0,3546 0,0084
4060 10 0,3846 0,3916 0,0070
3706 11 0,4231 0,4712 0,0481
3682 12 0,4615 0,4769 0,0154
3220 13 0,5000 0,5920 0,0920
3130 14 0,5385 0,6152 0,0767
2737 15 0,5769 0,7159 0,1390
2675 16 0,6164 0,7314 0,1160
2489 17 0,6538 0,7764 0,1226
2414 18 0,6923 0,7938 0,1015
2367 19 0,7308 0,8045 0,0737
2350 20 0,7692 0,8083 0,0391
2246 21 0,8077 0,8308 0,0231
2230 22 0,8462 0,8342 0,0120
2070 23 0,8846 0,8661 0,0185
1804 24 0,9231 0,9114 0,0117
1796 25 0,9615 0,9126 0,0489

28