Informe

Tipo Prueba 1:

Alumno: Manuel Ignacio González

Profesor: Marco Díaz
Asignatura: Diseño de Obras Hidráulicas

Domingo, 17 de mayo de 2020

 Índice

Contenido
Introducción ................................................................................................................................. 3
Marco Teórico .............................................................................................................................. 3
Distribución Gumbel ................................................................................................................. 3
Distribución Normal ................................................................................................................. 4
Distribución Log-Normal .......................................................................................................... 5
Distribución Gamma ................................................................................................................. 6
Prueba de Chi-Cuadrado........................................................................................................... 7
Planteamiento del problema ........................................................................................................ 8
Desarrollo ..................................................................................................................................... 9
Distribución Gumbel ............................................................................................................... 10
Distribución Normal ............................................................................................................... 12
Distribución Log-normal......................................................................................................... 14
Distribución Gamma ............................................................................................................... 17
Resultados y Conclusiones ......................................................................................................... 20

pág. 2
 Introducción

El presente informe tiene como finalidad la resolución de un problema de hidrología, el cual,

en forma resumida consiste en realizar un análisis de frecuencias para una cierta cantidad de
caudales medios anuales considerando las siguientes distribuciones de frecuencias:
 Distribución Gumbel;
 Distribución Normal
 Distribución LogNormal
 Distribución Gamma
Dado lo anterior, se obtendrán caudales para diversos periodos de retorno dados en el
problema.
Todo esto con el objetivo de buscar cual se ajusta mejor utilizando la prueba de chi-cuadrado.

Marco Teórico

Distribución Gumbel: En teoría de probabilidad y estadística la “distribución de Gumbel” es

utilizada para modelar la distribución del máximo (o el mínimo), por lo que se usa para
calcular valores extremos. Por ejemplo, sería muy útil para representar la distribución del
máximo nivel de un río a partir de los datos de niveles máximos durante 10 años.
Su función densidad es:

En la práctica se usa la relación:

Donde y y σn se denominan variables reducidas y se encuentran tabuladas en función del

número de datos de la muestra.

pág. 3
Distribución Normal: En estadística y probabilidad se llama distribución normal,
distribución de Gauss, distribución gaussiana o distribución de Laplace-Gauss, a una de las
distribuciones de probabilidad de variable continua que con más frecuencia aparece en
estadística y en la teoría de probabilidades.
La gráfica de su función de densidad tiene una forma acampanada y es simétrica respecto de
un determinado parámetro estadístico. Esta curva se conoce como campana de Gauss y es el
gráfico de una función gaussiana.
Su función densidad es:

Z: Variable reducida o tipificada

La variable z también es llamada variable estandarizada media cero y desviación estándar
uno.

pág. 4
Distribución Log-Normal: En probabilidades y estadísticas, la distribución normal
logarítmica es una distribución de probabilidad continua de una variable aleatoria cuyo
logaritmo está normalmente distribuido. Es decir, si X es una variable aleatoria con una
distribución normal, entonces exp(X) tiene una distribución log-normal.
La base de una función logarítmica no es importante, ya que log a X está distribuida
normalmente si y solo si logb X está distribuida normalmente, solo se diferencian en un factor
constante.
Se define la variable y como: y = log x

En que:
μy = μ (log x)
σy = σ (log x)
En este caso se debe “reducir” o normalizar la variable normal según la siguiente
transformación:

Finalmente se debe transformar a la variable original x = 10y

pág. 5
Distribución Gamma: En estadística la distribución gamma es una distribución de
probabilidad continua con dos parámetros κ>0 y λ>0 reales cuya función de densidad para
valores x>0 es:

Donde Г() es la función gamma completa, definida

por:

Esta función no es integrable analíticamente, por lo que es necesario aplicar tablas numéricas
ya calculadas. Los parámetros  y  cumplen las siguientes relaciones:

pág. 6
Prueba de Chi-Cuadrado: La prueba chi-cuadrado, también llamada Ji cuadrado (Χ2), se
encuentra dentro de las pruebas pertenecientes a la estadística descriptiva, concretamente la
estadística descriptiva aplicada al estudio de dos variables. La prueba chi-cuadrado es una de
las más conocidas y utilizadas para analizar variables nominales o cualitativas, es decir, para
determinar la existencia o no de independencia entre dos variables. Que dos variables sean
independientes significa que no tienen relación, y que por lo tanto una no depende de la otra,
ni viceversa.
 Permite definir si la función de densidad de frecuencia es apropiada para la muestra.
 Se comparará el valor de la variable obtenida a partir de la muestra X 2v con el valor
para la distribución, con un cierto nivel de confianza X2v,A.
Si X2v  X2v,A se acepta la hipótesis de que los datos tienen la distribución elegida.

Si X2v  X2v,A se rechaza la hipótesis, ya que es muy poco probable.

pág. 7
 Planteamiento del problema

Considerando la información de caudales medios mensuales registrada en la estación

fluviométrica Aconcagua en Chacabuquito, realizar el análisis de frecuencia para los
caudales medios anuales, considerando las siguientes distribuciones de frecuencia: Normal.
Log-Normal, Gumbel y Gamma.
Según lo anterior, determinar los caudales asociados a los Periodos de Retorno 5, 25, 50, 100,
200 y 250 años.
Determinar la distribución de mejor ajuste utilizando el test de chi-cuadrado, considere los
intervalos de clase necesarios para definir el histograma de la muestra.
Caudales Medios Mensuales Qmm (m3/s) "Registrados en la Estación Aconcagua en Chacabuquito"
AÑO ABRIL MAYO JUNIO JULIO AGOSTOSEPTIEMBRE OCTUBRE NOVIEMBRE DICIEMBRE ENERO FEBRERO MARZO MEDIO
1950-51 13,9 14 9,2 7,4 12,4 13,4 23,4 43 76,1 43,4 28,4 20,2 25,4
1951-52 14,4 11,5 12,9 17,7 15,6 16,3 22,2 47,6 57,7 43,4 26 18 25,3
1952-53 13,8 13,4 13,8 12,2 12,7 19,2 26,2 59 77,2 47,4 37,9 23,9 29,7
1953-54 16,7 14,1 15,9 14,1 27,8 42,1 44,5 125 179 126 77 37,5 60
1954-55 20,1 14,2 14,8 12 14,9 19,7 24 57,7 61,3 54,1 33,9 21 29
1955-56 13,5 13,7 16,1 16,3 15,1 17,3 23,6 54,2 54,9 36,8 30,1 20,1 26
1956-57 13,3 15,5 9,7 9,3 11,8 14,8 21,2 35,4 30,5 32,5 25,4 19,5 19,9
1957-58 13,7 16,7 17,1 14,2 16,3 16,6 24,2 43,3 60 42,3 24,6 21 25,8
1958-59 14,5 11,7 16,1 13,3 10,9 14,9 35,4 42,9 40,1 29,7 27,1 20,4 23,1
1959-60 13,7 18,8 14,9 14,5 18,9 22,3 25,8 49,4 69,3 41,4 31,4 19,7 28,3
1960-61 12,2 8,6 13,1 11,3 13,4 16,4 30,6 73,2 83,7 44,6 32,1 28,3 30,6
1961-62 14,6 10,9 13 12,9 14,5 18,1 45,5 85,7 96,7 52,8 47 26,3 36,5
1962-63 16,1 11,7 12,5 15,2 15,3 16,5 32,9 77,3 67,8 40,9 31,4 20,9 29,9
1963-64 12,7 10 10,4 15,2 13,6 27,7 33,2 43,3 149 134,4 63,7 32,2 45,5
1964-65 19 13,6 10,8 9,3 8,7 11,3 14,2 18,3 20,2 25,6 23,1 19,3 16,1
1965-66 15 8,6 7,4 10,1 25,8 27,7 47,8 95,6 83,7 91,5 46,7 26,6 40,5
1966-67 15,8 11,6 10,9 14 13,4 20,2 29,2 52,2 53,4 46,3 35,6 20,6 26,9
1967-68 14 10 8,9 9,1 6,8 9,1 16,2 22,2 33,2 25,8 22,6 13,8 16
1968-69 9,1 6,9 6,2 5,9 5,9 6,6 7,4 15,5 16,2 21,9 19,3 9,1 10,8
1969-70 6,4 5,6 7,3 4,9 6,2 8,7 11,2 36,2 71,1 40,8 31,3 19 20,7
1970-71 11,4 9,4 8,6 10,7 10,9 12,4 23,4 47,3 46,4 32,2 26,5 15,3 21,2
1971-72 9,6 7,4 5,7 9,9 11,7 15,7 30,3 63,3 45 36,4 26,6 14,1 23
1972-73 10 17,5 25,4 20,7 25,3 36,4 41,4 75,7 189 179 97,3 57,7 64,6
1973-74 28,4 18,8 16,9 18,7 16,9 14,8 23,8 52,9 57,2 64,8 40,8 27,3 31,8
1974-75 17,3 12,7 15,7 20,2 18,6 18,4 40,1 65,9 72,3 72,8 39,6 24,1 34,8
1975-76 15,4 11,2 10,3 11,4 13,8 17,8 23,7 32,1 53,9 40 25,9 18,5 22,8
1976-77 13,6 9,6 9,4 8,1 8,2 11,1 17,3 52,6 55,6 43,8 26,5 20,9 23,1
1977-78 11,7 10,6 10,9 18,3 22,7 42,8 57,5 88,7 123 74,1 46,7 27,1 44,5
1978-79 16,6 13 11,4 24,7 23,1 22,8 45 90,3 168 125 56,3 32,4 52,4
1979-80 18,6 15,1 11,3 10,5 14,1 18,6 28,8 38,8 59,9 62,9 40,4 31,7 29,2
1980-81 37,4 24,4 23 22,4 23,9 25,8 38,8 65,9 111 70,5 54,8 33,9 44,3
1981-82 19,1 18,9 15,4 10,9 12,4 12,6 20,1 37,3 33,1 34,4 26,2 16,4 21,4
1982-83 13,8 14,1 27,2 39,9 32,8 45 50,6 88,7 182 166 113 54,7 69
1983-84 28,9 18,9 15,5 17,8 19,8 20 54,4 91,4 116 75 49,7 26,5 44,5
1984-85 17,3 12,7 10 17,3 18,6 24,9 65 89,5 118 101 63,2 47,1 48,7
1985-86 24,6 18,9 15,6 16,3 14,5 12,5 21 58 57,4 45,1 34,6 23,3 28,5
1986-87 16,3 14,5 47,6 24,9 22,2 25,6 42 78,3 148 101 61,5 37,7 51,6
1987-88 19,2 14,5 18,4 39,8 50 35,6 55 147 176 150 79 44,8 69,1
1988-89 24,7 15,2 11 9,5 10,4 10,5 18,5 33,9 33 32,3 30,2 17,4 20,6
1989-90 10,7 9 6,5 6,7 14,8 23,4 40 78,6 59,5 37,8 26,5 16,7 27,5
1990-91 10,5 7,7 7,1 7 7,3 11,4 19 33,5 29,5 25,5 20,2 15,7 16,2
1991-92 11,9 13,6 16,2 27,4 20,8 33,9 38,3 77,6 85,1 80,2 49 32,8 40,6
1992-93 18,9 16,7 21,1 17,7 16,5 22 41,4 65,8 74,7 65 38 26,2 35,3
1993-94 19,8 38 31,9 16,3 14,6 17,5 27,7 45,7 54,5 44,1 27 21,7 29,9
1994-95 13,8 11,5 10,5 11,7 14,6 18,7 24,6 50,2 56,8 39,5 26 18,2 24,7
1995-96 13,4 12 10,7 10,5 9,9 16 18,2 41,8 39,4 25,5 22,8 19,3 20
1996-97 11,9 9,1 8 7,4 6,5 6,6 8,6 12,9 16,2 18,2 15,4 13,3 11,2
1997-98 9,5 7,7 25,5 22 28,6 37,5 39,2 12,9 128 149,1 71,1 46,8 48,2

pág. 8
 Desarrollo

El primer paso es ordenar los datos de forma correcta para que no exista confusión con los
caudales y años; paso seguido, calcular el caudal medio anual (Qma).
Año ENERO FEBRERO MARZO ABRIL MAYO JUNIO JULIO AGOSTO SEPTIEMBRE OCTUBRE NOVIEMBRE DICIEMBRE Qm a
1950 13,9 14 9,2 7,4 12,4 13,4 23,4 43 76,1 212,8
1951 43,4 28,4 20,2 14,4 11,5 12,9 17,7 15,6 16,3 22,2 47,6 57,7 307,9
1952 43,4 26 18 13,8 13,4 13,8 12,2 12,7 19,2 26,2 59 77,2 334,9
1953 47,4 37,9 23,9 16,7 14,1 15,9 14,1 27,8 42,1 44,5 125 179 588,4
1954 126 77 37,5 20,1 14,2 14,8 12 14,9 19,7 24 57,7 61,3 479,2
1955 54,1 33,9 21 13,5 13,7 16,1 16,3 15,1 17,3 23,6 54,2 54,9 333,7
1956 36,8 30,1 20,1 13,3 15,5 9,7 9,3 11,8 14,8 21,2 35,4 30,5 248,5
1957 32,5 25,4 19,5 13,7 16,7 17,1 14,2 16,3 16,6 24,2 43,3 60 299,5
1958 42,3 24,6 21 14,5 11,7 16,1 13,3 10,9 14,9 35,4 42,9 40,1 287,7
1959 29,7 27,1 20,4 13,7 18,8 14,9 14,5 18,9 22,3 25,8 49,4 69,3 324,8
1960 41,4 31,4 19,7 12,2 8,6 13,1 11,3 13,4 16,4 30,6 73,2 83,7 355
1961 44,6 32,1 28,3 14,6 10,9 13 12,9 14,5 18,1 45,5 85,7 96,7 416,9
1962 52,8 47 26,3 16,1 11,7 12,5 15,2 15,3 16,5 32,9 77,3 67,8 391,4
1963 40,9 31,4 20,9 12,7 10 10,4 15,2 13,6 27,7 33,2 43,3 149 408,3
1964 134,4 63,7 32,2 19 13,6 10,8 9,3 8,7 11,3 14,2 18,3 20,2 355,7
1965 25,6 23,1 19,3 15 8,6 7,4 10,1 25,8 27,7 47,8 95,6 83,7 389,7
1966 91,5 46,7 26,6 15,8 11,6 10,9 14 13,4 20,2 29,2 52,2 53,4 385,5
1967 46,3 35,6 20,6 14 10 8,9 9,1 6,8 9,1 16,2 22,2 33,2 232
1968 25,8 22,6 13,8 9,1 6,9 6,2 5,9 5,9 6,6 7,4 15,5 16,2 141,9
1969 21,9 19,3 9,1 6,4 5,6 7,3 4,9 6,2 8,7 11,2 36,2 71,1 207,9
1970 40,8 31,3 19 11,4 9,4 8,6 10,7 10,9 12,4 23,4 47,3 46,4 271,6
1971 32,2 26,5 15,3 9,6 7,4 5,7 9,9 11,7 15,7 30,3 63,3 45 272,6
1972 36,4 26,6 14,1 10 17,5 25,4 20,7 25,3 36,4 41,4 75,7 189 518,5
1973 179 97,3 57,7 28,4 18,8 16,9 18,7 16,9 14,8 23,8 52,9 57,2 582,4
1974 64,8 40,8 27,3 17,3 12,7 15,7 20,2 18,6 18,4 40,1 65,9 72,3 414,1
1975 72,8 39,6 24,1 15,4 11,2 10,3 11,4 13,8 17,8 23,7 32,1 53,9 326,1
1976 40 25,9 18,5 13,6 9,6 9,4 8,1 8,2 11,1 17,3 52,6 55,6 269,9
1977 43,8 26,5 20,9 11,7 10,6 10,9 18,3 22,7 42,8 57,5 88,7 123 477,4
1978 74,1 46,7 27,1 16,6 13 11,4 24,7 23,1 22,8 45 90,3 168 562,8
1979 125 56,3 32,4 18,6 15,1 11,3 10,5 14,1 18,6 28,8 38,8 59,9 429,4
1980 62,9 40,4 31,7 37,4 24,4 23 22,4 23,9 25,8 38,8 65,9 111 507,6
1981 70,5 54,8 33,9 19,1 18,9 15,4 10,9 12,4 12,6 20,1 37,3 33,1 339
1982 34,4 26,2 16,4 13,8 14,1 27,2 39,9 32,8 45 50,6 88,7 182 571,1
1983 166 113 54,7 28,9 18,9 15,5 17,8 19,8 20 54,4 91,4 116 716,4
1984 75 49,7 26,5 17,3 12,7 10 17,3 18,6 24,9 65 89,5 118 524,5
1985 101 63,2 47,1 24,6 18,9 15,6 16,3 14,5 12,5 21 58 57,4 450,1
1986 45,1 34,6 23,3 16,3 14,5 47,6 24,9 22,2 25,6 42 78,3 148 522,4
1987 101 61,5 37,7 19,2 14,5 18,4 39,8 50 35,6 55 147 176 755,7
1988 150 79 44,8 24,7 15,2 11 9,5 10,4 10,5 18,5 33,9 33 440,5
1989 32,3 30,2 17,4 10,7 9 6,5 6,7 14,8 23,4 40 78,6 59,5 329,1
1990 37,8 26,5 16,7 10,5 7,7 7,1 7 7,3 11,4 19 33,5 29,5 214
1991 25,5 20,2 15,7 11,9 13,6 16,2 27,4 20,8 33,9 38,3 77,6 85,1 386,2
1992 80,2 49 32,8 18,9 16,7 21,1 17,7 16,5 22 41,4 65,8 74,7 456,8
1993 65 38 26,2 19,8 38 31,9 16,3 14,6 17,5 27,7 45,7 54,5 395,2
1994 44,1 27 21,7 13,8 11,5 10,5 11,7 14,6 18,7 24,6 50,2 56,8 305,2
1995 39,5 26 18,2 13,4 12 10,7 10,5 9,9 16 18,2 41,8 39,4 255,6
1996 25,5 22,8 19,3 11,9 9,1 8 7,4 6,5 6,6 8,6 12,9 16,2 154,8
1997 18,2 15,4 13,3 9,5 7,7 25,5 22 28,6 37,5 39,2 12,9 128 357,8
1998 149,1 71,1 46,8 267

pág. 9
Distribución Gumbel
Qma Para el caso de la distribución Gumbel,
Qm a Ordenado
Año Orden P>X F(x) X(Gumbel)
primero debemos hallar las variables
1950 212,8 1 0,02 0,98 771,024857 755,7
reducidas y y σ, las cuales se hallan en la
1951 307,9 2 0,04 0,96 689,6813145 716,4
siguiente tabla.
1952 334,9 3 0,06 0,94 641,58534 588,4

1953 588,4 4 0,08 0,92 607,0885591 582,4

1954 479,2 5 0,1 0,9 580,0327855 571,1
1955 333,7 6 0,12 0,88 557,6744743 562,8
1956 248,5 7 0,14 0,86 538,5496559 524,5
1957 299,5 8 0,16 0,84 521,7840208 522,4
1958 287,7 9 0,18 0,82 506,8131893 518,5
1959 324,8 10 0,2 0,8 493,2513855 507,6
1960 355 11 0,22 0,78 480,8230343 479,2 A través del uso de la interpolación,
1961 416,9 12 0,24 0,76 469,324244 477,4
1962 391,4 13 0,26 0,74 458,5997553 456,8
encontramos que, para los 49 datos, y es
1963 408,3 14 0,28 0,72 448,5284527 450,1 igual a 0,54806 y σ es igual a 1,15892.
1964 355,7 15 0,3 0,7 439,0138818 440,5
1965 389,7 16 0,32 0,68 429,9778273 429,4 Media= 383,173469 Calculamos la
1966 385,5 17 0,34 0,66 421,3558359 416,9
Desv.Est= 134,02057 media, la
1967 232 18 0,36 0,64 413,0940146 414,1
1968 141,9 19 0,38 0,62 405,1466875 408,3 y 0,54806 desviación estándar,
1969 207,9 20 0,4 0,6 397,474647 395,2
σ= 1,15892 alfa () y μ.
1970 271,6 21 0,42 0,58 390,0438225 391,4
1971 272,6 22 0,44 0,56 382,8242475 389,7 α= 0,009
1972 518,5 23 0,46 0,54 375,7892457 386,2 u= 319,794363
1973 582,4 24 0,48 0,52 368,9147744 385,5
Para cada tiempo de
1974 414,1 25 0,5 0,5 362,1788855 357,8
1975 326,1 26 0,52 0,48 355,5612731 355,7
retorno tenemos como resultado los
1976 269,9 27 0,54 0,46 349,0428807 355 siguientes caudales Q m3/s.
1977 477,4 28 0,56 0,44 342,6055517 339
1978 562,8 29 0,58 0,42 336,2317057 334,9 T F(x) y Q(m3/s)=X
1979 429,4 30 0,6 0,4 329,9040247 333,7
5 0,2 1,49994 493,2513855
1980 507,6 31 0,62 0,38 323,6051362 329,1
1981 339 32 0,64 0,36 317,3172761 326,1 25 0,04 3,1985343 689,6813145
1982 571,1 33 0,66 0,34 311,0219127 324,8 50 0,02 3,9019387 771,024857
1983 716,4 34 0,68 0,32 304,6993099 307,9
100 0,01 4,6001492 851,7677716
1984 524,5 35 0,7 0,3 298,3279945 305,2
1985 450,1 36 0,72 0,28 291,8840862 299,5 200 0,005 5,2958121 932,2160688
1986 522,4 37 0,74 0,26 285,3404224 287,7
250 0,004 5,5194576 958,0790176
1987 755,7 38 0,76 0,24 278,6653785 272,6
1988 440,5 39 0,78 0,22 271,8212254 271,6
1989 329,1 40 0,8 0,2 264,7617653 269,9
1990 214 41 0,82 0,18 257,4288025 267
1991 386,2 42 0,84 0,16 249,7466537 255,6
1992 456,8 43 0,86 0,14 241,6131812 248,5
1993 395,2 44 0,88 0,12 232,8842472 232
1994 305,2 45 0,9 0,1 223,3446481 214
1995 255,6 46 0,92 0,08 212,6480336 212,8
1996 154,8 47 0,94 0,06 200,1739009 207,9
1997 357,8 48 0,96 0,04 184,6043943 154,8

1998 267 49 0,98 0,02 162,0514822 141,9

pág. 10
 Gráfico de distribución de Gumbel.

Como último paso, debemos realizar la prueba de bondad de ajuste por el método de chi-
cuadrado usando 6 intervalos.
Prueba de Bondad de ajuste
Rango F.Relativa F.Acum Frecuencia
Intervalos m3/s ni xi m3/s y Ơy Y Fs(xi) Fs(xii) F(xi) Valor Teorico(xi) x²
1 0-142 1 141,9 0,463 0,822 -1,02 0,020408163 0,0204082 0,06246213 0,062462134 1,38737312
2 142 - 284 11 236,97273 0,511 1,013 -0,59 0,224489796 0,244898 0,16464092 0,102178787 7,17408367
3 284-426 21 354,4619 0,526 1,071 0,3 0,428571429 0,6734694 0,47672369 0,31208277 2,13055905
4 426-568 11 488,10909 0,511 0,964 1,27 0,224489796 0,8979592 0,75515548 0,278431786 0,51207219
5 568-710 3 580,63333 0,47 0,851 1,72 0,06122449 0,9591837 0,83605059 0,080895118 0,23437443
6 710-852 2 736,05 0,466 0,836 2,67 0,040816327 1 0,9330913 0,097040707 1,59621534
0,933091301 13,0346778

pág. 11
Distribución Normal
Orden Año Qma Qma ordenado Distribución Normal
1 1950 212,8 141,9 0,000588802
u= 383,173469 Calculamos el
2 1951 307,9 154,8 0,000696968
3 1952 334,9 207,9 0,001265718 σ= 134,02057 promedio y la
4 1953 588,4 212,8 0,001326822 desviación estándar
5 1954 479,2 214 0,001341957 de la muestra.
6 1955 333,7 232 0,001575616
7 1956 248,5 248,5 0,001796679 T F(x) Z Q(x)(m3/s)
8 1957 299,5 255,6 0,00189226 5 0,2000 0,841621234 495,968027
9 1958 287,7 267 0,002044439 25 0,0400 1,750686071 617,801415
10 1959 324,8 269,9 0,002082661 50 0,0200 2,053748911 658,418069
11 1960 355 271,6 0,00210494 100 0,0100 2,326347874 694,951938
12 1961 416,9 272,6 0,002117997 200 0,0050 2,575829304 728,387581
13 1962 391,4 287,7 0,002309616 250 0,0040 2,652069808 738,605377
14 1963 408,3 299,5 0,002449607
15 1964 355,7 305,2 0,002513249 Para cada año de retorno, calculamos su
16 1965 389,7 307,9 0,002542365 respectivo caudal haciendo uso de la
17 1966 385,5 324,8 0,002707346
18 1967 232 326,1 0,002718681 variable z.
19 1968 141,9 329,1 0,002744033
20 1969 207,9 333,7 0,00278066
21 1970 271,6 334,9 0,002789754
22 1971 272,6 339 0,002819345
23 1972 518,5 355 0,002911673
24 1973 582,4 355,7 0,002914832
25 1974 414,1 357,8 0,002923851
26 1975 326,1 385,5 0,002976276
27 1976 269,9 386,2 0,002975965
28 1977 477,4 389,7 0,002973197
29 1978 562,8 391,4 0,002971122
30 1979 429,4 395,2 0,002964763
31 1980 507,6 408,3 0,002924866
32 1981 339 414,1 0,002898515
33 1982 571,1 416,9 0,002883945
34 1983 716,4 429,4 0,002804817
35 1984 524,5 440,5 0,002716491
36 1985 450,1 450,1 0,00262777
37 1986 522,4 456,8 0,00255978
38 1987 755,7 477,4 0,002324875
39 1988 440,5 479,2 0,002302817
40 1989 329,1 507,6 0,001934497
41 1990 214 518,5 0,001787882
42 1991 386,2 522,4 0,001735377
43 1992 456,8 524,5 0,001707148
44 1993 395,2 562,8 0,001212439
45 1994 305,2 571,1 0,001113726
46 1995 255,6 582,4 0,000986026
47 1996 154,8 588,4 0,000921617
48 1997 357,8 716,4 0,000135304
49 1998 267 755,7 6,25164E-05

pág. 12
Finalizando con su respectivo prueba de bondad de ajuste.

Prueba de Bondad de ajuste

Rango F.Relativa F.Acum Fs Variable Valor
Intervalo m3/s ni Anual xi Fs(xi) (xi) normal Z Frec. F(xi) Teorico p(x) x2
1 0-142 1 141,9 0,02040816 0,02040816 -1,80027192 0,03590886 0,03590886 0,32786626
2 142 - 284 11 236,9727273 0,2244898 0,24489796 -1,09088285 0,13766222 0,10175336 7,25428034
3 284-426 21 354,4619048 0,42857143 0,67346939 -0,21423252 0,41518287 0,27752065 4,02853089
4 426-568 11 488,1090909 0,2244898 0,89795918 0,78298146 0,783181 0,36799813 2,74223533
5 568-710 3 580,6333333 0,06122449 0,95918367 1,4733549 0,92967232 0,14649132 2,43189271
6 710-852 2 736,05 0,04081633 1 2,63300276 0,99576832 0,066096 0,4737659
49 1 17,2585714

Gráfico de la distribución normal

hecho a partir de los datos de los
caudales de retorno.

Gráfico de la función densidad,

elaborado a partir de los datos de
cada año.

pág. 13
Distribución Log-normal

Función
Funcion distribucion Distribucion
Año Qma Qma ordenados Log Densidad probabilidad teorica Diferencia
1 1950 212,8 141,9 2,1519824 0,09067075 0,004903757 0,02 0,01509624
2 1951 307,9 154,8 2,18977096 0,164154539 0,009604537 0,04 0,03039546
3 1952 334,9 207,9 2,31785449 0,796564768 0,063721378 0,06 0,00372138
4 1953 588,4 212,8 2,32797162 0,87709363 0,072184091 0,08 0,00781591
5 1954 479,2 214 2,33041377 0,89716195 0,074350541 0,1 0,02564946
6 1955 333,7 232 2,36548798 1,208774115 0,111164577 0,12 0,00883542
7 1956 248,5 248,5 2,39532639 1,497515948 0,151508202 0,14 0,0115082
8 1957 299,5 255,6 2,40756085 1,617948252 0,170566384 0,16 0,01056638
9 1958 287,7 267 2,42651126 1,802082842 0,202980006 0,18 0,02298001
10 1959 324,8 269,9 2,43120288 1,84664442 0,211539446 0,2 0,01153945
11 1960 355 271,6 2,43392977 1,87227825 0,216610024 0,22 0,00338998
12 1961 416,9 272,6 2,43552585 1,887181866 0,219610244 0,24 0,02038976
13 1962 391,4 287,7 2,45893986 2,094838753 0,266272184 0,26 0,00627218
14 1963 408,3 299,5 2,47639683 2,231753344 0,304063206 0,28 0,02406321
15 1964 355,7 305,2 2,48458453 2,289213418 0,322574584 0,3 0,02257458
16 1965 389,7 307,9 2,48840969 2,314397032 0,331379705 0,32 0,01137971
17 1966 385,5 324,8 2,51161602 2,441854551 0,386654965 0,34 0,04665496
18 1967 232 326,1 2,5133508 2,449501881 0,390897714 0,36 0,03089771
19 1968 141,9 329,1 2,51732788 2,465983958 0,400672851 0,38 0,02067285
20 1969 207,9 333,7 2,52335621 2,488122334 0,415607072 0,4 0,01560707
21 1970 271,6 334,9 2,52491515 2,493279359 0,419489958 0,42 0,00051004
22 1971 272,6 339 2,5301997 2,508993044 0,432708548 0,44 0,00729145
23 1972 518,5 355 2,55022835 2,543075827 0,483369751 0,46 0,02336975
24 1973 582,4 355,7 2,55108387 2,5436168 0,48554562 0,48 0,00554562
25 1974 414,1 357,8 2,55364034 2,544782225 0,492049937 0,5 0,00795006
26 1975 326,1 385,5 2,58602438 2,501319051 0,574045781 0,52 0,05404578
27 1976 269,9 386,2 2,58681227 2,498941289 0,576015606 0,54 0,03601561
28 1977 477,4 389,7 2,59073041 2,486217077 0,585782358 0,56 0,02578236
29 1978 562,8 391,4 2,59262082 2,479546909 0,590476091 0,58 0,01047609
30 1979 429,4 395,2 2,59681694 2,46352427 0,60084752 0,6 0,00084752
31 1980 507,6 408,3 2,61097938 2,397485728 0,635290721 0,62 0,01529072
32 1981 339 414,1 2,61710523 2,36348654 0,649874832 0,64 0,00987483
33 1982 571,1 416,9 2,6200319 2,346148402 0,656766762 0,66 0,00323324
34 1983 716,4 429,4 2,63286204 2,262301156 0,686343549 0,68 0,00634355
35 1984 524,5 440,5 2,64394591 2,180487716 0,71097268 0,7 0,01097268
36 1985 450,1 450,1 2,65330901 2,105465614 0,731041503 0,72 0,0110415
37 1986 522,4 456,8 2,6597261 2,051312686 0,744379806 0,74 0,00437981
38 1987 755,7 477,4 2,67888241 1,878971532 0,782047415 0,76 0,02204742
39 1988 440,5 479,2 2,68051681 1,863666519 0,7851059 0,78 0,0051059
40 1989 329,1 507,6 2,70552161 1,622325414 0,828711827 0,8 0,02871183
41 1990 214 518,5 2,71474876 1,531511166 0,843262463 0,82 0,02326246
42 1991 386,2 522,4 2,71800317 1,499468557 0,848194478 0,84 0,00819448
43 1992 456,8 524,5 2,71974549 1,482327544 0,850792104 0,86 0,0092079
44 1993 395,2 562,8 2,75035409 1,187060592 0,891607497 0,88 0,0116075
45 1994 305,2 571,1 2,75671216 1,128122313 0,898966965 0,9 0,00103303
46 1995 255,6 582,4 2,76522137 1,051086591 0,908237047 0,92 0,01176295
47 1996 154,8 588,4 2,76967266 1,011718419 0,912827878 0,94 0,02717212
48 1997 357,8 716,4 2,85515558 0,415647689 0,971529614 0,96 0,01152961
49 1998 267 755,7 2,87834942 0,31018603 0,979902726 0,98 9,7274E-05

pág. 14
De lo anterior podemos calcular la desviación estándar y la media, y con eso los caudales para
cada año de retorno (junto a su respectivo gráfico).

u 2,556764
desv.est 0,1567376

T F(x) z P(x) Qma

5 0,2000 0,84162123 2,68867768 488,2898323
25 0,0400 1,75068607 2,83116232 677,8948333
50 0,0200 2,05374891 2,87866367 756,2470031
100 0,0100 2,32634787 2,92139017 834,4305086
200 0,0050 2,5758293 2,9604933 913,0473399
250 0,0040 2,65206981 2,97244305 938,5189571

Y por ultimo la prueba de chi-cuadrado.

Prueba de Bondad de ajuste
F.Relativa F.Acum Variable Frecuencia Valor
Intervalo Rango m3/s ni xi m3/s Log(x) Fs(xi) Fs(xi) Normal Z F(x) teorico p(x) x2
1 0-142 1 141,9 2,2 0,02 0,02 -2,58 0,00 0,00 2,40202403
2 142 - 284 11 236,972727 2,4 0,22 0,24 -1,16 0,12 0,12 4,73522914
3 284-426 21 354,461905 2,5 0,43 0,67 -0,05 0,48 0,36 0,66075122
4 426-568 11 488,109091 2,7 0,22 0,90 0,84 0,80 0,32 1,34781586
5 568-710 3 580,633333 2,8 0,06 0,96 1,32 0,91 0,11 0,96385422
6 710-852 2 736,05 2,9 0,04 1,00 1,98 0,98 0,07 0,57158198
49 1,00 0,98 10,6812565

pág. 15
Por ultimo también es posible graficar la función densidad y la función distribución de
probabilidades.

pág. 16
Distribución Gamma
A partir de los caudales medios anuales, calculamos la distribución gamma como se aprecia
en la tabla. También su media, desviación estándar y sus variables alfa y beta.
Año Qma Qma ordenado Distr.Gamma u= 383,173469
1 1950 212,8 141,9 0,000407008
2 1951 307,9 154,8 0,00057702 Desv.Est= 134,02057
3 1952 334,9 207,9 0,001542171 α= 8,17425045
4 1953 588,4 212,8 0,001641791
β= 46,8756704
5 1954 479,2 214 0,001666175
6 1955 333,7 232 0,002025747
7 1956 248,5 248,5 0,002332297
8 1957 299,5 255,6 0,002453448
9 1958 287,7 267 0,002630952
Y calculamos los caudales para cada periodo de
10 1959 324,8 269,9 0,002672414 retorno.
11 1960 355 271,6 0,002695982
12 1961 416,9 272,6 0,002709585 T P(x) Qma
13 1962 391,4 287,7 0,002890561
5 0,2 489,0588204
14 1963 408,3 299,5 0,002998457
15 1964 355,7 305,2 0,003039688 10 0,1 561,8125368
16 1965 389,7 307,9 0,003056727 25 0,04 646,7211119
17 1966 385,5 324,8 0,003127299
18 1967 232 326,1 0,003130191
50 0,02 705,6915071
19 1968 141,9 329,1 0,003135516 100 0,01 761,549544
20 1969 207,9 333,7 0,003140071 200 0,005 815,0400379
21 1970 271,6 334,9 0,003140552
22 1971 272,6 339 0,003140024 250 0,004 831,842585
23 1972 518,5 355 0,003107323
24 1973 582,4 355,7 0,003104836
25 1974 414,1 357,8 0,003096873
26 1975 326,1 385,5 0,002928325
27 1976 269,9 386,2 0,002922714
28 1977 477,4 389,7 0,002893804
29 1978 562,8 391,4 0,002879263
30 1979 429,4 395,2 0,002845634
31 1980 507,6 408,3 0,002719035
32 1981 339 414,1 0,002658442
33 1982 571,1 416,9 0,002628345
34 1983 716,4 429,4 0,002488392
35 1984 524,5 440,5 0,002358294
36 1985 450,1 450,1 0,002242998
37 1986 522,4 456,8 0,002161685
38 1987 755,7 477,4 0,001911478
39 1988 440,5 479,2 0,00188981
40 1989 329,1 507,6 0,001558433
41 1990 214 518,5 0,001438466
42 1991 386,2 522,4 0,001396737
43 1992 456,8 524,5 0,001374543
44 1993 395,2 562,8 0,001006722
45 1994 305,2 571,1 0,000936755
46 1995 255,6 582,4 0,00084719
47 1996 154,8 588,4 0,000802281
48 1997 357,8 716,4 0,000214638
49 1998 267 755,7 0,000136144

pág. 17
Por último, la prueba de chi-cuadrado y los respectivos gráficos.

Prueba de Bondad de ajuste

Frecuencia Frecuencia Valor Teorico
Intervalos Rango m3/s ni xi m3/s Fs(xi) F.Acum Fs(xii) F(xii) p(xii) x²
1 0-142 1 141,9 0,02040816 0,020408163 0,01033859 0,010338586 0,48057084
2 142 - 284 11 236,972727 0,2244898 0,244897959 0,12485329 0,114514707 5,17515621
3 284-426 21 354,461905 0,42857143 0,673469388 0,45930517 0,334451874 1,29784304
4 426-568 11 488,109091 0,2244898 0,897959184 0,7983115 0,339006329 1,8955038
5 568-710 3 580,633333 0,06122449 0,959183673 0,91762576 0,119314268 1,38580824
6 710-852 2 736,05 0,04081633 1 0,98622122 0,068595457 0,55123657
49 1 0,986221221 10,7861187

Distribución gamma vs Qma.

Distribución para los distintos

periodos de tiempo.

pág. 18
Con cada distribución lista y como dato adicional, agregamos un gráfico comparando todas
las distribuciones.

pág. 19
 Resultados y Conclusiones

Como paso final debemos verificar cuál de todos los resultados de cada distribución se adapta más a
la prueba de chi-cuadrado, como se ve en la siguiente tabla:
2
DistribucionesX v X2v,A Resultado
Normal 17,2585714 11,0705 No cumple
Log-Normal 10,6812565 11,0705 Si cumple
Gumbel 13,0346778 11,0705 No cumple
Gamma 10,7861187 11,0705 Si cumple

Tomando un valor de 95% de confianza y haciendo uso de la tabla de distribución de chi-

cuadrado tenemos el siguiente valor:

X2v,A = 11,0705
Valor también obtenible por medio de una función de Excel.
Observando la tabla podemos darnos cuenta que hay dos resultados que, si cumplen con la
prueba, pero solo podemos escoger uno. Eso se decide escogiendo el resultado menor,
quedándonos, así como resultado que la distribución de Log-normal es la que más se adapta
a nuestra prueba, con X2v = 10,68125646.
Como dato adicional de nuestro resultado, agregamos su respectivo histograma:

pág. 20