Metodos Estadisticos
Metodos Estadisticos
Metodos Estadisticos
METODOS ESTADÍSTICOS.
Los métodos estadísticos, se basan en considerar que el caudal máximo anual, es una variable
aleatoria que tiene una cierta distribución. Para utilizarlos se requiere tener como datos, el registro
de caudales máximos anuales, cuanto mayor sea el tamaño del registro, mayor será también la
aproximación del cálculo de caudal de diseño, el cual se calcula para un determinado periodo de
retorno.
Por lo general, en los proyectos donde se desea determinar el caudal de diseño, se cuenta con pocos
años de registro por lo que la curva de distribución de probabilidades de los caudales máximos, se
tiene que prolongar en su extremo si se requiere inferir un caudal con un periodo de retorno mayor
al tamaño del registro. El problema se origina en que existen muchos tipos de distribuciones que se
apegan a los datos, y sin embargo, difieren en los extremos. Esto ha dado lugar a diversos métodos
estadísticos, dependiendo del tipo de distribución que considere.
Gumbel
Nash
Levediev
Log-Pearson III
Gumbel y Nash consideran una distribución de valores extremos, con la única diferencia, que el
criterio de Nash es menos rígido que el de Gumbel, pues permite ajustar la distribución por mínimos
cuadrados. Por otra parte, Levediev considera una distribución Pearson tipo III. En forma práctica, se
recomienda escoger varias distribuciones y ver cual se ajusta mejor; esto requiere que se tengan los
datos necesarios para poder aplicar alguna prueba estadística, como la prueba de bondad de ajuste.
MÉTODO DE GUMBEL.
Para calcular el caudal máximo para un periodo de retorno determinado se usa la ecuación:
σQ
Q max =Q m− (Y −lnT )…1
σN N
Siendo:
σQ=
Donde:
√ ∑ Q2i −N Q 2m …2
i=1
N−1
T= Periodo de retorno.
Para calcular el intervalo de confianza, o sea, aquel dentro del cual puede variar Q max dependiendo
del registro disponible se hace lo siguiente:
1. Si ф=1-1/T varía entre 0.20 y 0.80, el intervalo de confianza se calcula con la fórmula:
σQ
∆ Q=± √ Nα σ m …3
σN √N
Donde:
El caudal máximo de diseño para un cierto periodo de retorno, será igual al caudal máximo con la
ecuación (1), más el intervalo de confianza, calculado con (3) ó (4).
Q d =Q max + ∆ Q …5
Se tiene el registro de caudales máximos de 30 años para la estación 9-3 Angostura, como se muestra
en la tabla 6.15.
Calcular el caudal de diseño para el vertedor de demasías, para períodos de retorno 50 y 100 años
respectivamente.
Caudal Caudal
Año(1) m3/s(2) Año(1) m3/s(2)
1970 1660 1985 563
1971 917 1986 520
1972 3800 1987 360
SOLUCIÓN.
Caudal
Año(1) m3/s(2) Q2(M3/S)
1970 1660 2755600
1971 917 840889
1972 3800 14440000
1973 1410 1988100
1974 2280 5198400
1975 618 381924
1976 683 466489
1977 934 872356
1978 779 606841
1979 921 848241
1980 876 767376
1981 740 547600
1982 1120 1254400
1983 610 372100
1984 1150 1322500
1985 563 316969
1986 520 270400
1987 360 129600
1988 367 134689
1989 658 432964
1990 824 678976
1991 850 722500
1992 1230 1512900
1993 522 272484
Paso 1.
28749
Qm = =958.3m3/s Q m2=958.32=¿918338.89
30
Paso 2.
σQ=
√ ∑ Q2i −N Q 2m
i=1
N−1
40595065−30 (918338.89)
σQ=
√ 30−1
=670.6893
Paso 3.
σN 1.11238
YN 0.53622
Paso 4.
σQ
Qmax =Qm− (Y −lnT )
σN N
Para T=50
670.6893
Q max =958.3− (0.53622−ln 50)
1.11238
Q max =2993.68m3/s
Para T=100
670.6893
Q max =958.3− (0.53622−ln 100)
1.11238
Paso 5.
Cálculo de ф.
ф=1-1/T
Para T=50años
ф=1-1/50=0.98
Para T=100años
T=1-1/100=0.99
Paso 6.
Cálculo del intervalo de confianza. Como en ambos casos vemos que ф es mayor que 0.90,
Utilizaremos la ecuación:
1.14 σ Q
∆ Q=±
σN
1.14∗670.6893
∆ Q=± =687.34 m3/s
1.11238
Paso 7.
Q d =Q max + ∆ Q
Para T=50
Q d =2993.68+687.34
Qd =3681.02 m3/ s
Para T=100
Qd =3411.60+ 687.34
Q d =4098.94 m3/ s
MÉTODO DE NASH
Nash considera que el valor del caudal para un determinado periodo de retorno se puede calcular
con la ecuación:
T
Q max =a+ bloglog …6
T +1
Donde:
Los parámetros a y b se estiman utilizando el método de mínimos cuadrados, con la ecuación lineal:
Q=a+bX , utilizando las siguientes ecuaciones:
a=Qm−b X m … 7
∑ X i Q i−N X m Qm
b= i=1 N …8
2 2
∑ X i −N X m
i=1
Siendo:
T
X i =loglog …9
T+1
Donde:
Para calcular los valores de X i correspondientes a los Q i se ordenan estos en forma decreciente,
asignándole a cada uno número de orden m i; al Q i máximo le corresponderá el valor 1, al inmediato
siguiente 2, etc. Entonces, el valor del periodo de retorno para Q i se calculará utilizando la formula
de Weibull con la ecuación:
N +1
T= … 10
mi
El valor dentro del cual puede variar el Q max calculado por la ecuación (6), se obtiene como:
1
∆ Q=±2
Siendo:
√
N 2 ( N −1)
+(
S qq
X− X m )2 N−2
S xx
∗1
S qq (
−
S xq2
S xx
… 11)
2
S xx =N ∑ x i2−( ∑ x i )
2
Sqq=N ∑ Qi2−( ∑ Qi )
S xq=N ∑ Q i x i−¿ ( ∑ Q i )( ∑ x i ) ¿
De la ecuación (11), se ve que ∆ Q solo varía con X, la cual se calcula de la ecuación (9), sustituyendo
el valor del periodo de retorno para el cual se calculó el Q max . Todos los demás términos que
intervienen en la ecuación (11) se obtienen de los datos.
Q d =Q max + ∆ Q
Para los mismos datos de la tabla 6.15, del ejemplo 6.7, calcular el caudal de diseño utilizando el
método de Nash para periodo de retorno de 50 y 100 años.
Solución:
1.- Ordenando en forma descendente, los valores de los caudales de la columna 2, de la tabla 6.15 se
obtiene la tabla 6.16
T
Q max =a+ bloglog T =(1+ N )/m X= log log
T +1 T
T −1
m (1) Caudal m3/s(2) T (3) T/(T-1) (4) X (5) Q*X (6) Q2 (7)
X2 (8)
31.000
1 3800 0 1.0333 -1.8465 -7016.61 14440000 3.4095
15.500
2 2280 0 1.0690 -1.5381 -3506.97 5198400 2.3659
10.333
3 1660 3 1.1071 -1.3545 -2248.54 2755600 1.8348
4 1410 7.7500 1.1481 -1.2219 -1722.83 1988100 1.4930
5 1230 6.2000 1.1923 -1.1170 -1373.88 1512900 1.2476
6 1150 5.1667 1.2400 -1.0296 -1183.99 1322500 1.0600
7 1120 4.4286 1.2917 -0.9541 -1068.58 1254400 0.9103
8 1030 3.8750 1.3478 -0.8873 -913.90 1060900 0.7873
9 953 3.4444 1.4091 -0.8270 -788.12 908209 0.6839
10 934 3.1000 1.4762 -0.7717 -720.81 872356 0.5956
11 921 2.8182 1.5500 -0.7205 -663.57 848241 0.5191
N
a=Qm−b X m
∑ X i Q i−N X m Qm
b= i=1 N
∑ X i2−N X m2
i=1
b= -1206.30
a= 240.44
T
Q max =a+ bloglog
T +1
2 2
S xx =N ∑ x i2−( ∑ x i ) Sqq=N ∑ Qi2−( ∑ Qi )
1 X= log log
∆ Q=±2
√ S qq
N 2 ( N −1)
+( X− X m )2 N−2
S xx
∗1
S qq −
S xq2
(
S xx
El valor de X se calcula para cada periodo de retorno
)
T
T −1
T= 50 años X= -2.0568
T= 100 años X= -2.3600
MÉTODO DE LEVEDIEV.
Este método está basado en suponer que los caudales máximos anuales son variables aleatorias
Pearson tipo III. El caudal de diseño se obtiene a partir de la fórmula:
Qd =Qmax + ∆ Q …12
Donde:
Qmax =Qm ( K C v +1 ) … 13
A Er Qmax
∆ Q=± … 14
√N
Los términos que aparecen en las ecuaciones anteriores tienen el siguiente significado.
A=Coeficiente que varía de 0.7 a 1.5, dependiendo del número de años del registro. Cuantos más
años de registro haya, menor será el valor del coeficiente. Si N es mayor de 40años, se toma el valor
de 0.7.
N 3
Qi
C s=
∑
i=1
( Qm
−1 ) … 15
N C v3
N 2
Qi
C v=
√ ∑
i =1
( Qm
N
−1 ) …16
N= Años de observación.
∑ Qi
Qm = i=1 …17
N
Q max =Caudal máximo probable obtenido para un periodo de retorno determinado, en m3/s.
Para los mismos datos de la tabla 6.15, del ejemplo 6.7, calcular el caudal de diseño utilizando el
método de Levediev para periodo de retorno de 50 y 100 años.
SOLUCION
1.- Obtención del Caudal Medio
N
Qi
Qm =∑
i=1 N
2 3
Q Q Q
Qm
-1 ( Qm
−1 ) ( Qm
−1 )
Q (m3/seg)
AÑO
1970 1660 0.7322 0.5362 0.3926
1971 917 -0.0431 0.0019 -0.0001
1972 3800 2.9654 8.7933 26.0754
1973 1410 0.4714 0.2222 0.1047
1974 2280 1.3792 1.9022 2.6236
1975 618 -0.3551 0.1261 -0.0448
1976 683 -0.2873 0.0825 -0.0237
1977 934 -0.0254 0.0006 0.0000
1978 779 -0.1871 0.0350 -0.0065
1979 921 -0.0389 0.0015 -0.0001
1980 876 -0.0859 0.0074 -0.0006
1981 740 -0.2278 0.0519 -0.0118
1982 1120 0.1687 0.0285 0.0048
1983 610 -0.3635 0.1321 -0.0480
1984 1150 0.2000 0.0400 0.0080
1985 563 -0.4125 0.1702 -0.0702
1986 520 -0.4574 0.2092 -0.0957
1987 360 -0.6243 0.3898 -0.2434
1988 367 -0.6170 0.3807 -0.2349
1989 658 -0.3134 0.0982 -0.0308
1990 824 -0.1401 0.0196 -0.0028
1991 850 -0.1130 0.0128 -0.0014
1992 1230 0.2835 0.0804 0.0228
1993 522 -0.4553 0.2073 -0.0944
1994 581 -0.3937 0.1550 -0.0610
1995 557 -0.4188 0.1754 -0.0734
1996 818 -0.1464 0.0214 -0.0031
1997 1030 0.0748 0.0056 0.0004
1998 418 -0.5638 0.3179 -0.1792
1999 953 -0.0055 0.0000 0.0000
SUMATORI
A 28749 14.2049 28.0063
PROMEDIO 958.30
N= 30
N 2
Qi
C v=
√ ∑
i =1
( Qm
N
−1 )
Cv= 0.6881
C s=
∑
i=1
( Qm
−1 )
N C v3
Cs=
2.8652
Cs=
2.0643
Cs=
2.8652
A=Coeficiente que varía de 0.7 a 1.5, dependiendo del número de años del registro. Cuantos más
años de registro haya, menor será el valor del coeficiente. Si N es mayor de 40años, se toma el valor
de 0.7.
A= 0.85
A Er Q max
∆ Q=±
√N
Distribución estándar para análisis de frecuencia de caudales máximos anuales en los Estados Unidos
(Benson 1968).
La transformación Qd = Log QT se usa para reducir la asimetría; en caso de que la asimetría para esta
situación valga cero la distribución log Pearson III se reduce a una log normal.
Qd=log Q T
Siendo:
´
log Q T = LogQ+ K σ log Q
LogQ=
´ ∑ logQi / N
Donde:
σ Log Q = desviación estándar de los logaritmos de la serie Qi, cuya fórmula es:
´ 2 /(N −1) ]1 /2
σ logQ= [∑ (log Q i − LogQ)
Para los mismos datos de la tabla 6.15, del ejemplo 6.7, calcular el caudal de diseño utilizando el
método de Log - Pearson III, para periodo de retorno de 50 y 100 años.
SOLUCION
N= 30
´
LogQ = 2.9174 m3/seg
3.- Calculo de Desviación estándar de los logaritmos de la serie Qi, cuya fórmula es (𝝈 𝑳𝒐𝒈 𝑸)
1 /2
´ 2 /(N −1) ]
σ log Q= [∑ (log Q i − LogQ)
σ Log Q= 0.2216
Cs logQ 0.9366
´
log Q T = LogQ+ K σ LogQ Qd=log Q T
log QT
Periodo de Retorno Qd unidad
T= 50 años 3.4744 2980.93 m3/seg
T= 100 años 3.5777 3782.21 m3/seg
4000.00
3500.00
3000.00
GUMBEL
CAUDAL (M3/SEG)
2500.00
Linear (GUMBEL)
NASH
2000.00
Linear (NASH)
LEVEDIEV
1500.00
Linear (LEVEDIEV)
1000.00 LOG - PEARSON III
Linear (LOG - PEARSON III)
500.00 REGISTRO
0.00
0 10 20 30 40 50 60 70 80 90 100
TIEMPO (AÑOS)
En el gráfico T vs. Q, se observa que la distribución que más se acerca a la distribución registrada, es
la distribución por el Método de Levediev, por lo cual asumiremos esta distribución para calcular el
Qd.
CAUDAL DE DISEÑO
T (años) Qd (m3/s)
50 3460.28
100 4149.90
BIBLIOGRAFÍA
Villón Bejar, Máximo: hidrología. Segunda Edición: editorial Villón, Febrero del 2002. Lima-
Perú
http://docs.google.com/viewer?
a=v&q=cache:QFuPMyK8k50J:intranet.catie.ac.cr/intranet/posgrado/Hidro2006/Presentacio
nes/Capitulo%25206b.ppt
http://ocw.upm.es/ingenieria-agroforestal/climatologia-aplicada-a-la-ingenieria-y-
medioambiente/contenidos/tema-7/METODO-DE-GUMBEL.pdf