MÉTODOS

ESTADÍSTICOS
Docente:
Ms. Ing. Edgar Gustavo Sparrow Alamo
Alumnos:
 Flores Utrilla Keneth
 Blas Huatoro José
 Condori Villanueva Carlos
 Cabello López JeanMarco
 Colchado Ircañaupa Brayan
 Métodos 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 se requiere tener como datos, el

registro de los caudales máximos anuales, cuanto mayor sea el tamaño del registro, mayor será también

la aproximación del cálculo del 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 quiere inferir un caudal con un periodo de retorno mayor al tamaño de

registro. El problema se origina, en que existen muchos tipos de distribuciones que se apegan a los datos y

que, sin embargo, difieren en los extremos. Esto ha dado lugar a diversos métodos estadísticos,

dependiendo del tipo de distribución que se considere.

 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 para poder aplicarlo 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:
…1

Siendo:
…2

Donde:
Caudal máximo para un periodo de retorno determinado, en m3/s.
N= número de años de registro.
Qi =Caudales máximos anuales registrados, en m3/s.
, Caudal promedio, en m3/s
T= Periodo de retorno.

Constantes función de N, tabla 6.13 (Variables reducidas)

= Desviación estándar de los caudales.
Para calcular el intervalo de confianza, o sea, aquel dentro del cual puede
variar 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:
 

Donde:
N= número de años de registro
=constante en función de ф, tabla 6.14

= Constantes función de N, tabla 6.13

= Desviación estándar de los caudales
(ecuación 2)
 Tabla 6.14 Valores de en función de ф.
Tabla 6.13. Valores de y en función de N.
La zona de ф comprendida entre 0.8 y 0.9
se considera la transición, donde es
proporcional al cálculo con las ecuaciones 3
y 4, dependiendo del valor de ф.
 
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).
 
 Año(1) Caudal m3/s(2) Año(1) Caudal m3/s(2)
EJEMPLO DEL METODO DE GUMBEL. 1970 1660 1985 563
1971 917 1986 520
Se tiene el registro de caudales máximos de 1972 3800 1987 360
30 años para la estación 9-3 Angostura, 1973 1410 1988 367
como se muestra en la tabla 6.15. 1974 2280 1989 658
1975 618 1990 824
En este río se desea construir una presa de
1976 683 1991 850
almacenamiento.
1977 934 1992 1230
Calcular el caudal de diseño para el vertedor 1978 779 1993 522
de demasías, para períodos de retorno 50 y 1979 921 1994 581
100 años respectivamente. 1980 876 1995 557
1981 740 1996 818
1982 1120 1997 1030
1983 610 1998 418
1984 1150 1999 953
 Año(1) Caudal m3/s(2) Q2(M3/S)
1970 1660 2755600
1971 917 840889
1972 3800 14440000
1973 1410 1988100
1974 2280 5198400
1975 618 381924 Paso 1.
1976 683 466489
1977
1978
934
779
872356
606841
Calcular el caudal promedio.
1979 921 848241 𝑁
𝑄
1980 876
𝑄 767376
𝑚 =∑
1981
1982
740
1120
547600
1254400 𝑖= 1 𝑁
SOLUCIÓN.
1983 610 372100
1984 1150 1322500
1985 563 316969 =958.3m3/s
1986 520 270400
1987 360 129600
1988 367 134689
1989 658 432964
1990 824 678976 918338.89
1991 850 722500
1992 1230 1512900
1993 522 272484
1994 581 337561
1995 557 310249
1996 818 669124
1997 1030 1060900
1998 418 174724
1999 953 908209
 Paso 3.
Paso 2. Cálculo de los coeficientes σN, YN

Cálculo de la Desviación estándar de los caudales.

√
𝑁

∑ 2
𝑄 𝑖 − 𝑁 𝑄𝑚
2

𝑖 =1
𝜎𝑄=
𝑁 −1

𝜎 𝑄 = √¿ ¿ ¿

σN 1.11238
YN
0.53622
Paso 4.
Cálculo del Caudal Máximo.

¿
Para los periodos de retorno de 50 y 100 años.

¿
 Para T=50

m3/s

 Para T=100

¿
m3/s
 Paso 6.
Paso 5.
Cálculo del intervalo de confianza.
Cálculo de ф.
Como en ambos casos vemos que
ф=1-1/T ф es mayor que 0.90, Utilizaremos
Para T=50años la ecuación:
ф=1-1/50=0.98 1.14 𝜎 𝑄
∆ 𝑄 =±
Para T=100años 𝜎 𝑁

T=1-1/100=0.99
m3/s
Paso 7.
Cálculo del caudal de diseño.

 Para T=50

 Para T=100
MÉTODO DE
NASH
MÉTODO DE NASH

■ Nash considera que el valor del caudal para un determinado periodo de retorno se puede calcular con la
ecuación:

Donde:
a, b = constantes en función del registro de caudales máximos anuales.
Qmáx =caudal máximo para un periodo de retorno determinado, en m3/s
T =periodo de retorno, en años.
Los parámetros se estiman utilizando el método de mínimos cuadrados, con la ecuación lineal: Q = a + bx ,
utilizando las siguientes ecuaciones:

Siendo:

Donde:

N= número de años de registro

=caudales máximos anuales registrados, en m3/s

, Caudal medio, en m3/s.

= constante para cada caudal registrado, en función de su periodo de retorno

correspondiente.

, valor medio de las

Para calcular los valores de correspondientes a los se ordenan estos en forma decreciente, asignándole a cada
uno número de orden ; al máximo le corresponderá el valor 1, al inmediato siguiente 2, etc. Entonces, el valor del
periodo de retorno para se calculará utilizando la formula de Weibull con la ecuación:

Finalmente el valor de se obtiene sustituyendo el valor (10) en (9).

El valor dentro del cual puede variar el calculado por la ecuación (6), se obtiene como:

Siendo:
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 Qmáx . Todos los demás términos que intervienen en la ecuación (11) se
obtienen de los datos.
El caudal máximo de diseño correspondiente a un determinado periodo de retorno será igual al caudal máximo
obtenido de la ecuación (6), más el intervalo de confianza calculado según la ecuación (11), es decir:
 EJEMPLO DEL METODO DE NASH

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.

TABLA 6.15. CAUDALES MÁXIMOS

Año(1) Caudal m3/s(2) Año(1) Caudal m3/s(2)

1970 1660 1985 563

1971 917 1986 520
1972 3800 1987 360
1973 1410 1988 367
1974 2280 1989 658
1975 618 1990 824
1976 683 1991 850
1977 934 1992 1230
1978 779 1993 522
1979 921 1994 581
1980 876 1995 557
1981 740 1996 818
1982 1120 1997 1030
1983 610 1998 418
1984 1150 1999 953
 Caudal
m (1) m3/s(2) T (3) T/(T-1) (4) X (5) Q*X (6) Q2 (7)
X2 (8)
Solución: 1 3800 31.0000 1.0333 -1.8465 -7016.61 14440000 3.4095
2 2280 15.5000 1.0690 -1.5381 -3506.97 5198400 2.3659
1. Ordenando en forma descendente, los valores de los caudales de la columna
3 1660 10.3333 1.1071 -1.3545 -2248.54 2755600 1.8348
2, de la tabla 6.15 se obtiene la tabla 6.16 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
2. Cálculos Preliminares
6 1150 5.1667 1.2400 -1.0296 -1183.99 1322500 1.0600
Columna 1: número de orden 7 1120 4.4286 1.2917 -0.9541 -1068.58 1254400 0.9103
Columna 2: caudales máximos ordenados de mayor a menor 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
Columna 3: periodo de retorno obtenido con la fórmula de Weibull
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
12 917 2.5833 1.6316 -0.6724 -616.61 840889 0.4521
13 876 2.3846 1.7222 -0.6269 -549.19 767376 0.3930
Columna 4: cociente 14 850 2.2143 1.8235 -0.5835 -495.98 722500 0.3405
15 824 2.0667 1.9375 -0.5418 -446.40 678976 0.2935
16 818 1.9375 2.0667 -0.5013 -410.08 669124 0.2513
Columna 5 :
17 779 1.8235 2.2143 -0.4619 -359.81 606841 0.2133
18 740 1.7222 2.3846 -0.4232 -313.15 547600 0.1791
19 683 1.6316 2.5833 -0.3849 -262.90 466489 0.1482
Columna 6 : producto de Q*X
20 658 1.5500 2.8182 -0.3468 -228.21 432964 0.1203
 Q  28749 21 618 1.4762 3.1000 -0.3086 -190.71 381924 0.0952

 Q  40595.065
2
22 610 1.4091 3.4444 -0.2699 -164.66 372100 0.0729

De la tabla, se tiene:  X  17.8528 23 581 1.3478 3.8750 -0.2304 -133.88 337561 0.0531
24 563 1.2917 4.4286 -0.1896 -106.74 316969 0.0359
 X  17.6256
2

25 557 1.2400 5.1667 -0.1468 -81.76 310249 0.0215

 QX  25554.28 26 522 1.1923 6.2000 -0.1011 -52.75 272484 0.0102
27 520 1.1481 7.7500 -0.0510 -26.49 270400 0.0026
3. Cálculo del Qm y Xm 25554.28  30 x( 0.5951) x958.3
b
17.6256  30 x(0.5951) 2
b  1206.3152

a  958.3( 1206.3152) x( 0.5951)

a  240.4218
28749 3
Qm   958.3 m
30 s 5. Cálculo del Caudal Máximo.
17.8528
Xm   0.5951
30

4. Cálculo de los Parámetros a y b; Sustituyendo los valores de los parámetros a y b,

en la ecuación (6), se tiene:
T
Qmax  240.4218  1206.3152 log log
T 1
Luego:

Para T = 50 años, Qmax = 2721.5783 m3/s

Para T = 100 años, Qmax = 3087.3680 m3/s

6. Cálculo de las desviaciones estándar y covarianza 8. Cálculo del Caudal de Diseño

Para un T = 50 años  X= -2.0568

Para un T = 100 años  X= -2.3600

S xx  30 x17.6256  (17.8528) 2  210.0455

SQQ  30 x 40595065  (28749) 2  391346949 9. Sustituyendo en la ecuación:
S XQ  30 x(25554.28) 2  28749 x(17.8528)  253378.2528 Q  2 14994.1360  14571.0472( X  0.5951) 2

Para un X= -2.0568  ΔQ = 429.5412

7. Cálculo del Intervalo de Confianza
Para un X= -2.3600  ΔQ = 491.4586

10. Cálculo del caudal de diseño:

391346949 1 1  (253378.2528) 2 
Q  2 2
 ( X  0.5951) x x  391346949   Para T = 50 años, Qd = 2721.5783 + 429.5412 = 3 151.12m3/s
302 x 29 28 210.0455  210.0455 

2
Para T = 100 años, Qd = 3087.3680 + 491.4586 = 3 578.83m3/s
Q  2 14994.1360  14571.0472( X  0.5951)
 Método de Lebediev

Este método está basado en suponer que los caudales máximos anuales son
variables aleatorias. El caudal de diseño se obtiene a partir de la fórmula:

Donde:

caudal máximo probable obtenido para un periodo de retorno determinado, en

= intervalo de confianza, en
Donde:

K = coeficiente que depende de la probabilidad P = 1/T, expresada en porcentaje de que se repita el caudal de
diseño y del coeficiente de asimetría (Tabla 1)

= caudal promedio, en , el cual se obtiene de

= caudales máximos anuales observados, en

N = años de observación
Donde:

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 40 años, se toma el valor de 0.7.

=coeficiente que depende de los valores de (ecuación 3) y de la probabilidad P=1/T, su valor se encuentra de
la figura 6.3

N = años de observación

coeficiente de asimetría, se calcula como:

Por otra parte, Lebediev recomienda tomar los siguientes valores:

Entre estos valores y el que se obtiene de la ecuación 4, se escoge el mayor

Tabla 1

Valores de K
Figura 1: Valores de Er en función de Cv y p
 EJEMPLO MÉTODO DE LEBEDIEV

Para los mismos datos de la tabla 6.15, del ejemplo 6.7, calcular el caudal de diseño
utilizando el método de Lebediev para periodo de retorno de 50 y 100 años.
TABLA 6.15. CAUDALES MÁXIMOS
Año(1) Caudal m3/s (2) Año(1) Caudal m3/s(2)
1970 1660 1985 563
1971 917 1986 520
1972 3800 1987 360
1973 1410 1988 367
1974 2280 1989 658
1975 618 1990 824
1976 683 1991 850
1977 934 1992 1230
1978 779 1993 522
1979 921 1994 581
1980 876 1995 557
1981 740 1996 818
1982 1120 1997 1030
1983 610 1998 418
1984 1150 1999 953
Solución:

1. Obtención del caudal medio mediante la siguiente formula:

2. Cálculos previos
Con los datos de la tabla 6.15 y con el valor de Qm se obtiene, se obtiene la siguiente
tabla:
  
     
Año Caudal m3/s
     
 
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
∑ 28749 14.2049 28.0063
Promedio 958.30
3. Cálculo del coeficiente de Variación Cv

Aplicando la siguiente ecuación:

4. Determinación del coeficiente de asimetría Cs

De la siguiente ecuación se obtiene:

Por otra parte, Levediev recomienda tomar los siguientes valores:

 Cs = 2Cv Para avenidas producidas por deshielo.

 Cs = 3Cv Para avenidas producidas por tormentas.
 Cs = 5Cv Para avenidas producidas por tormentas en
cuencas ciclónicas.
Considerando que la avenida es producida por una tormenta, se tiene:
Cs = 3Cv = 3 x 0.6881 = 2.0643
De estos valores se escoge el mayor, por lo tanto, se tiene:
⸫
5. Obtención del coeficiente k:
 Para el periodo de retorno de 50 años, el valor
de P es:

Con P=2% y Cs=2.8654, de la tabla

mostrada se obtiene k=3.12
 Para el periodo de retorno de 100 años, el valor
de P es:

Con P=1% y Cs=2.8654, de la tabla

mostrada se obtiene k=3.98
6. Cálculo de Er
Para ello se utiliza el siguiente nomograma:

 Para P=2% y Cv=0.6881,

se obtiene Er=0.95

 Para P=1% y Cv=0.6881,

se obtiene Er=1.02
7. Cálculo del caudal máximo
Se utiliza la siguiente ecuación:

 Para T=50 años, se tiene:

 Para T=100 años, se tiene:

8. Cálculo del intervalo de confianza

Para N=30 años se puede tomar A=0.85 y se utiliza la siguiente ecuación:
 Para T=50 años:

 Para T=100 años:

9. Cálculo del caudal de diseño:

 Para T=50 años