CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
6. Comparar el valor del estadstico
D , con el valor crtico D0
P( x)
Donde:
de la Tabla C-4
del Anexo C, con los siguientes criterios de decisin:
m N 1
Si:
D D0
m = posicin;
El ajuste es bueno, al nivel de significacin seleccionado
P( x) = Probabilidad Asignada, o posicin de ploteo
N = Nmero total de datos
D D0
m
3
El ajuste no es bueno, al nivel de significacin seleccionado,
La posicin de ploteo de los caudales para este ejemplo se muestra en la Tabla 10.9.
Tabla 10.9 Posicin de ploteo de caudales medios anuales
P(x)
siendo necesario probar con otra distribucin.
Esta prueba de ajuste no requiere del conocimiento a priori de la funcin de distribucin terica, es aplicable a distribuciones de datos no agrupados y de cualquier distribucin terica.
Comparndola con la prueba Chi-cuadrado, no requiere que la frecuencia absoluta de cada clase sea igual o mayor que 5, esta no es una prueba exacta, sino una prueba aproximada.
Ejemplo 10.1
Sea una serie de caudales medios anuales medidos desde 1983 a 1997 como se observa en la Tabla 10.8.
1988 26 1996 38
45
Tabla 10.8 Caudales Medios Anuales
1989 14 1997 42 28 1990
Ao
1983
1984
1985
1986
1987
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Q (m /s ) 11 13 14 15 22 24 26 28 29 33 33 38 38 42 44 0,0625 0,1250 0,1875 0,2500 0,3125 0,3750 0,4375 0,5000 0,5625 0,6250 0,6875 0,7500 0,8125 0,8750 0,9375
6,25% 12,50% 18,75% 25,00% 31,25% 37,50% 43,75% 50,00% 56,25% 62,50% 68,75% 75,00% 81,25% 87,50% 93,75%
Q ( m3 / s )
13
15
29
11
33
A partir de esto se grafica en los papeles de probabilidad de la ley normal y gumbel, como se observa en la Figura 10.13 y la Figura 10.14 respectivamente.
Papel de Probabilidades de la Ley Normal
Ao
1991
1992
1993
1994
1995
Q ( m3 / s )
33
44
38
22
24
Se pide:
40
a) Determinar la ley a la que mejor se ajusta esta serie utilizando los papales de probabilidad Normal y Gumbel.
Indistintamente a la ley que elegida, aplicar la ley normal y:
35
b) Calcular los caudales de aporte referidos a la ocurrencia de un ao seco, hmedo y un ao medio.
30
25
Considerar como ao seco aquel que tiene una probabilidad de 80% de ser excedido en el transcurso de un ao cualquiera. Un ao hmedo aquel que tiene una probabilidad de 90% de no ser excedido y un ao medio aquel que es excedido en promedio una vez cada dos aos.
20
3 c) Calcular el Periodo de Retorno de los caudales: 48, 25 y 15 m s .
15
d) Determinar el caudal para un periodo de retorno de 10, 20 y 50 aos.
Solucin:
10
0.01 0.05 0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 95.0 98.0 99.0 99.8 99.9
F(x)
Probabilidad de no Excedencia
a) Para determinar a qu ley corresponde, se deber ordenar en forma ascendente y asignar una probabilidad o posicin de ploteo, para este fin se aplicara la formula de Weibull de la Tabla 10.6, de donde se tiene: 255
Figura 10.13. Ajuste Grfico de los Caudales Medios anuales a la Ley Normal y su recta analtica Copyright 2009 by Agustn and Weimar UMSS F.C. y T. - ING. CIVIL
Copyright 2009 by Agustn and Weimar
UMSS F.C. y T. - ING. CIVIL
256
CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
Papel de Probabilidades de la Ley Gumbel
Con la probabilidad de no excedencia de P( X
x) )
0.9 , se determina la
45
variable reducida Z, con ayuda de la tabla de probabilidades acumuladas de donde se tiene:
40
Z 1.285
Z S
1.285
x 27.33 27 33 10.794 10 794
x 41.2
35
Se considera entonces como ao lluvioso la ocurrencia de 41.2 m3 / s ms. Para un ao medio, con un periodo de retorno de:
30
T x)
De los datos se calcula que:
P( P( X
X Q 27,33
x)
P( P( X
25
1 1 P( X
T 1 T
2 2 1 x) 0.5 50% 05 2
Sn 1 10,794 10
x) ) 0.5 , se determina la
20
Con la probabilidad de no excedencia de P( X
15
variable reducida Z, con ayuda de la tabla de probabilidades acumuladas de donde se tiene:
10 98.0 99.0 99.5 99.7 99.8 99.99
Z 0
Z
F(x)
x S
x 27.33
0.1 0.5 1.0
5.0
10.0
20.0 30.0 40.0 50.0 60.0
70.0
80.0
90.0
95.0 96.0 97.0
x 27.33 27.3 10.794 10 794
Probabilidad de no Excedencia
Figura 10.14. Ajuste Grfico de los Caudales Medios anuales a la Ley Gumbel y su recta
analitica
El caudal de aporte que se puede tener para un ao medio es de 27.33 m3 / s segn la ley normal. Ntese que el caudal de ao medio es directamente el valor medio estimado a partir de los 15 valores observados (ley normal)
Observando los papeles de probabilidad Normal y Gumbel, se concluye que los datos se ajustan mejor a una Ley Normal.
b) Para un ao seco, la probabilidad de excedencia P( X
x) =80% )
Aplicando la ecuacin 10.5, P( X
x) P( X
P( X P( x) 0.2
x) 1
Se tiene que la probabilidad de no excedencia es:
c) Para calcular el periodo de retorno de los caudales se determina la variable reducida Z, y con esta la probabilidad de no excedencia con ayuda de la tabla de probabilidades acumuladas, como sigue:
P( X
x) 0.8 1 )
P( X P(
x) 1 0.8
x S
Z 1,91495 0.972252
De los datos se calcula que:
X
x) ) 0.2 , se determina la
Q 27,33
Sn 1 10,794 10
T
Con la probabilidad de no excedencia de P( X
48 27.33 27.3 10 794 10.794 P( X x) P(
variable reducida Z, con ayuda de la tabla de probabilidades acumuladas de donde se tiene:
1 1 P( X
x)
1 0 972 1 0.972252
T 36.0
Considerando muchsimos aos, se espera que el valor de 48 igualado o excedido en promedio una vez cada 36 aos.
Z
x 18,2
0.845
0.845
x 27.33 27 33 10.794 10 794
m3 s , sea
El caudal de aporte define un ao seco es entonces 18.2 m3 / s valores menores ste son considerados ao seco.
x X S
25 27.33 27.3 10 794 10.794
-0,21586
P( X P(
x) 0.416834
Para un ao hmedo, la probabilidad de no excedencia P( X
x) =90% )
De los datos se calcula que:
Q 27,33
Sn 1 10,794 10
257
1 1 P( X
x)
Copyright 2009 by Agustn and Weimar
1 1 0.416834 0 416
T 1.71
UMSS F.C. y T. - ING. CIVIL
Copyright 2009 by Agustn and Weimar
UMSS F.C. y T. - ING. CIVIL
258
CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
Considerando muchsimos aos, se espera que el valor de 25 Contestar las mismas preguntas del ejercicio 10.1 utilizando la ley Gumbel. Solucin: Para aplicar la ley Gumbel se tiene:
x
m s , sea
Ejemplo 10.2
igualado o excedido en promedio una vez cada 1.715 aos.
15 27.33 27.3 10.794 10 794
-1,1425
P( X x) 0.127143 P(
T 1.15
m3 s , sea
Donde:
1 1 P( X
x)
6 S *10 79 10.79 95 *10.795 6
1 1 0 127 0.127143
P( x)
8.4168 8.4168
F ( x) e
Considerando muchsimos aos, se espera que el valor de 15
igualado o excedido en promedio una vez cada 1.15 aos.
X 0.57721
27.333 0.57721*8.4168 3 0.57721*
22.474 22.4747 47
d) Para T=10
De donde se tiene que la probabilidad de no excedencia se expresa como:
T
P( P( X x)
1 1 P( X
10 1 0.9 90% 09 10
x)
P( P( X
x)
T 1 T
P( x)
F ( x) e
x 22.4747 22.47 8 4168 8.4168 8
Con la ayuda de tabla de probabilidades acumuladas se determina la
variable reducida Z:
Z 1.285
1.285
3
x 27.33 27 33 10.794 10 794
x 41.2
a) No es necesario ajustar, porque en el ejemplo anterior ya se determino que corresponde a la ley Normal. b) Para un ao seco, la probabilidad de excedencia P( X Se tiene que la probabilidad de no excedencia es: Despejando x se tiene:
El caudal para un periodo de retorno de 10 aos es de 41.2 m
x) =80% )
P( X P( x) 0.2
Para T=20
*l *ln
ln ln P( x)
x 22.4747 8.4168*ln n x 22.4747 8.4168*ln ln 0.2
ln P( x)
1 1 P( X
x)
x S X
P( P( X
x)
T 1 T
P( P( X
x)
20 1 0.95 95% 0 95 20
x 18.5
Con la ayuda de tabla de probabilidades acumuladas se determina la
variable reducida Z:
Z 1.645
Z
El caudal de aporte que se puede tener para un ao seco es de 18.5 m3 / s segn la ley Gumbel. Para un ao hmedo, la probabilidad de no excedencia P( X
x) =90% )
1.645
3
x 27.33 27 33 10.794 10 794
x 45.1
s
*l *ln
ln P( x) ln
x 22.4747 8.4168*ln
ln P( x)
El caudal para un periodo de retorno de 20 aos es de 45.1 m
x 22.4747 8.4168*ln
m3 / s segn la ley Gumbel.
ln 0.9
x 41.4
El caudal de aporte que se puede tener para un ao hmedo es de 41.4
Para T=50
T
P( P( X x)
1 1 P( X
50 1 0 98 98% 0.98 50
x)
P( X P(
x)
T 1 T
2.055
Para un ao medio, con un periodo de retorno de:
T
1 T 1 P( X
x
2
x)
*l *ln x
P( P( X x) T 1 T P( X P( x) 2 1 0.5 50% 05 2
Con la ayuda de tabla de probabilidades acumuladas se determina la
variable reducida Z:
Z 2.055
m3 s
259
x 27.33 27 33 x 49.5 10.794 10 794
ln P( x) ln
x 22.4747 8.4168*ln
Copyright 2009 by Agustn and Weimar
22.4747 8.4168*ln ln 0.5
ln P( x)
El caudal para un periodo de retorno de 50 aos es de 49.5
25.6
UMSS F.C. y T. - ING. CIVIL
Copyright 2009 by Agustn and Weimar
UMSS F.C. y T. - ING. CIVIL
260
CAPITULO X
3
ESTADISTICA APLICADA A LA HIDROLOGIA
CAPITULO X
ESTADISTICA APLICADA A LA HIDROLOGIA
3
El caudal de aporte que se puede tener para un ao medio es de 25.6 m / s segn la ley Gumbel. El caudal para un periodo de retorno de 10 aos es de 41.4 m
Para T=20
c) Para un caudal de 48
m3 s .
x)
*l *ln ln ln P( x)
e
48 22 4747 22.47 22.4747 8.4168 8 4168 8
T
x
P( X P( x) 0.95296
1 1 P( X
P( P( X
x 22.4747 8.4168*ln
x)
T 1 T P( X P( x)
20 1 0 95 0.95 95% 20
ln P( x)
P( x) F ( x) e
x 22.4747 8.4168*ln
T
m3 s , sea
T x) 1 1 P( X
Para T=50
P( x) F ( x) e
ln 0.95
21.26 21
x 47.474 47.5
m3 s
1 1 P( X
x)
1 0 952 1 0.95296
El caudal para un periodo de retorno de 20 aos es de 47.5
Considerando decenas de aos se espera que el valor de 48
igualado o excedido en promedio una vez cada 21 aos.
3 Para un caudal de 25 m s .
25 22 4747 22.47 22.4747 8.4168 8 4168 8
P( X P(
x)
T 1 T
x
P( X P(
x)
50 1 50
0.98 98% 0 98
P( x) F ( x) e
P( X P( x) 0.476734 0.4767346
P( x) F ( x) e
x
x
T
m3 s , sea
10.6.- CUESTIONARIO
*l *ln
ln P( x) ln
22.4747 8.4168*ln
ln P( x)
22.4747 8.4168*ln
ln 0.98
x 55.317 55.3
m3 s
1 1 P( X
x)
1 0.4767346 1 0 476 0.476734
1.911 2
El caudal para un periodo de retorno de 50 aos es de 55.3
Considerando decenas de aos se espera que el valor de 25
igualado o excedido en promedio una vez cada 2 aos.
Para un caudal de 15
e
15 22 4747 22.47 22.4747 8 4168 8.4168 8
m3 s .
P( X P( x) 0.087999
Describa que entiende por funcin de densidad y funcin de distribucin acumulada Qu diferencia existe entre la funcin de distribucin de probabilidad y la funcin de distribucin de distribucin acumulada? Defina probabilidad de excedencia y probabilidad de no excedencia, utilice en su explicacin la grafica de la funcin de probabilidad acumulada Defina periodo de retorno de dos diferentes maneras Para qu nos sirve y en que se utiliza el periodo de retorno? Qu entiende por riesgo de fallo? Mencione al menos 5 ejemplos en ingeniera civil donde se utilizan el periodo de retorno y el riesgo de fallo Qu entiende por posicin de ploteo, y para que se utiliza? Qu es un papel de probabilidad y cul es su utilidad? Enuncie las funciones de distribucin de probabilidad ms usadas en hidrologa Cules son las aplicaciones de las funciones de probabilidad en hidrologa? Qu diferencia existe entre la funcin de distribucin de probabilidad Normal y Gumbel? 261
Copyright 2009 by Agustn and Weimar UMSS F.C. y T. - ING. CIVIL
P( x)
T
m3 s , sea
F ( x)
1.096 1
P( x) F ( x) e
1 1 P( X
x)
1 0 087 1 0.087999
Considerando decenas de aos se espera que el valor de 15
igualado o excedido en promedio una vez cada ao.
d) Para T=10
T
ln P( x)
1 1 P( X
x)
P( P( X
x)
T 1 T x)
P( P( X
10 1 0 9 90% 0.9 10
x
x 41.416 41.4
*l *ln
ln P( x) ln
22.4747 8.4168*ln
x 22.4747 8.4168*ln
ln 0.9
Copyright 2009 by Agustn and Weimar
UMSS F.C. y T. - ING. CIVIL
262
