Introducción

El presente informe se refiere al desarrollo de la prueba 1 de la asignatura en

cuestión, el cual se refiere al diseño de obras hidráulicas, en específico al diseño de
caudales medios anuales en la estación fluviométrica Aconcagua en Chacabuquito.

El diseño de estos caudales son diseños entendiendo el comportamiento de los

caudales a través de diferentes distribuciones, las cuales son; Distribución Normal,
Log Normal, Gumbel y Gamma.

Esto puede realizarse gracias a la ayuda del software Excel, herramienta

matemática indispensable que nos ayudara en todos los cálculos, minimizando de
esta manera el error.

Objetivo General
 realizar el análisis de frecuencia para los caudales medios anuales de una
estación fluviometrica.

Objetivos específicos.

 Entender la primera materia de la asignatura del ramo, correspondiente a

Hidrologia.
 Realizar cálculos y diseños a través de Excel.
 Enunciado
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 6 intervalos de clase para def inir el histograma de la
muestra

Caudales Medios Mensuales Qmm (m3/s) "Registrados en la Estación Aconcagua en Chacabuquito"

AGOS SEPTIE OCTUBR NOVIEM DICIEMB FEBRER
AÑO ABRIL MAYO JUNIO JULIO TO MBRE E BRE RE ENERO O MARZO MEDIO

1950-51 13,90 14,00 9,20 7,40 12,40 13,40 23,40 43,00 76,10 43,40 28,40 20,20 25,4

1951-52 14,40 11,50 12,90 17,70 15,60 16,30 22,20 47,60 57,70 43,40 26,00 18,00 25,3

1952-53 13,80 13,40 13,80 12,20 12,70 19,20 26,20 59,00 77,20 47,40 37,90 23,90 29,7

1953-54 16,70 14,10 15,90 14,10 27,80 42,10 44,50 125,00 179,00 126,00 77,00 37,50 60,0

1954-55 20,10 14,20 14,80 12,00 14,90 19,70 24,00 57,70 61,30 54,10 33,90 21,00 29,0

1955-56 13,50 13,70 16,10 16,30 15,10 17,30 23,60 54,20 54,90 36,80 30,10 20,10 26,0

1956-57 13,30 15,50 9,70 9,30 11,80 14,80 21,20 35,40 30,50 32,50 25,40 19,50 19,9

1957-58 13,70 16,70 17,10 14,20 16,30 16,60 24,20 43,30 60,00 42,30 24,60 21,00 25,8

1958-59 14,50 11,70 16,10 13,30 10,90 14,90 35,40 42,90 40,10 29,70 27,10 20,40 23,1

1959-60 13,70 18,80 14,90 14,50 18,90 22,30 25,80 49,40 69,30 41,40 31,40 19,70 28,3

1960-61 12,20 8,60 13,10 11,30 13,40 16,40 30,60 73,20 83,70 44,60 32,10 28,30 30,6

1961-62 14,60 10,90 13,00 12,90 14,50 18,10 45,50 85,70 96,70 52,80 47,00 26,30 36,5

1962-63 16,10 11,70 12,50 15,20 15,30 16,50 32,90 77,30 67,80 40,90 31,40 20,90 29,9

1963-64 12,70 10,00 10,40 15,20 13,60 27,70 33,20 43,30 149,00 134,40 63,70 32,20 45,5

1964-65 19,00 13,60 10,80 9,30 8,70 11,30 14,20 18,30 20,20 25,60 23,10 19,30 16,1

1965-66 15,00 8,60 7,40 10,10 25,80 27,70 47,80 95,60 83,70 91,50 46,70 26,60 40,5

1966-67 15,80 11,60 10,90 14,00 13,40 20,20 29,20 52,20 53,40 46,30 35,60 20,60 26,9

1967-68 14,00 10,00 8,90 9,10 6,80 9,10 16,20 22,20 33,20 25,80 22,60 13,80 16,0
1968-69 9,10 6,90 6,20 5,90 5,90 6,60 7,40 15,50 16,20 21,90 19,30 9,10 10,8

1969-70 6,40 5,60 7,30 4,90 6,20 8,70 11,20 36,20 71,10 40,80 31,30 19,00 20,7

1970-71 11,40 9,40 8,60 10,70 10,90 12,40 23,40 47,30 46,40 32,20 26,50 15,30 21,2

1971-72 9,60 7,40 5,70 9,90 11,70 15,70 30,30 63,30 45,00 36,40 26,60 14,10 23,0

1972-73 10,00 17,50 25,40 20,70 25,30 36,40 41,40 75,70 189,00 179,00 97,30 57,70 64,6

1973-74 28,40 18,80 16,90 18,70 16,90 14,80 23,80 52,90 57,20 64,80 40,80 27,30 31,8

1974-75 17,30 12,70 15,70 20,20 18,60 18,40 40,10 65,90 72,30 72,80 39,60 24,10 34,8

1975-76 15,40 11,20 10,30 11,40 13,80 17,80 23,70 32,10 53,90 40,00 25,90 18,50 22,8

1976-77 13,60 9,60 9,40 8,10 8,20 11,10 17,30 52,60 55,60 43,80 26,50 20,90 23,1

1977-78 11,70 10,60 10,90 18,30 22,70 42,80 57,50 88,70 123,00 74,10 46,70 27,10 44,5

1978-79 16,60 13,00 11,40 24,70 23,10 22,80 45,00 90,30 168,00 125,00 56,30 32,40 52,4

1979-80 18,60 15,10 11,30 10,50 14,10 18,60 28,80 38,80 59,90 62,90 40,40 31,70 29,2

1980-81 37,40 24,40 23,00 22,40 23,90 25,80 38,80 65,90 111,00 70,50 54,80 33,90 44,3

1981-82 19,10 18,90 15,40 10,90 12,40 12,60 20,10 37,30 33,10 34,40 26,20 16,40 21,4

1982-83 13,80 14,10 27,20 39,90 32,80 45,00 50,60 88,70 182,00 166,00 113,00 54,70 69,0

1983-84 28,90 18,90 15,50 17,80 19,80 20,00 54,40 91,40 116,00 75,00 49,70 26,50 44,5

1984-85 17,30 12,70 10,00 17,30 18,60 24,90 65,00 89,50 118,00 101,00 63,20 47,10 48,7

1985-86 24,60 18,90 15,60 16,30 14,50 12,50 21,00 58,00 57,40 45,10 34,60 23,30 28,5

1986-87 16,30 14,50 47,60 24,90 22,20 25,60 42,00 78,30 148,00 101,00 61,50 37,70 51,6

1987-88 19,20 14,50 18,40 39,80 50,00 35,60 55,00 147,00 176,00 150,00 79,00 44,80 69,1

1988-89 24,70 15,20 11,00 9,50 10,40 10,50 18,50 33,90 33,00 32,30 30,20 17,40 20,6

1989-90 10,70 9,00 6,50 6,70 14,80 23,40 40,00 78,60 59,50 37,80 26,50 16,70 27,5

1990-91 10,50 7,70 7,10 7,00 7,30 11,40 19,00 33,50 29,50 25,50 20,20 15,70 16,2

1991-92 11,90 13,60 16,20 27,40 20,80 33,90 38,30 77,60 85,10 80,20 49,00 32,80 40,6

1992-93 18,90 16,70 21,10 17,70 16,50 22,00 41,40 65,80 74,70 65,00 38,00 26,20 35,3

1993-94 19,80 38,00 31,90 16,30 14,60 17,50 27,70 45,70 54,50 44,10 27,00 21,70 29,9

1994-95 13,80 11,50 10,50 11,70 14,60 18,70 24,60 50,20 56,80 39,50 26,00 18,20 24,7

1995-96 13,40 12,00 10,70 10,50 9,90 16,00 18,20 41,80 39,40 25,50 22,80 19,30 20,0

1996-97 11,90 9,10 8,00 7,40 6,50 6,60 8,60 12,90 16,20 18,20 15,40 13,30 11,2

1997-98 9,50 7,70 25,50 22,00 28,60 37,50 39,20 12,90 128,00 149,10 71,10 46,80 48,2
 Marco Teórico
Para desarrollar el ejercicio, primero debemos entender cómo funcionan las distintas
distribuciones que nos piden en el enunciado y el chi-cuadrado.
Distribución Normal
La distribución normal es una distribución con forma de campana donde las
desviaciones estándar sucesivas con respecto a la media establecen valores de
referencia para estimar el porcentaje de observaciones de los datos.
Su función densidad es

En qué  y σ son el promedio y la desviación típica, además de la variable reducida

(z) es:

Distribución Log-Normal
La distribución log-normal es una distribución flexible que se relaciona
estrechamente con la distribución normal. Esta distribución puede resultar
particularmente útil para modelar datos que sean aproximadamente simétricos o
asimétricos a la derecha.
Su función densidad.

Donde:

se debe reducir la variable normal

y luego a
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. Su función densidad es.

Se usan las relaciones:

Los valores de promedio y desviación se encuentran tabulados por la siguiente tabla

según la cantidad de muestra.

Distribución Gamma
La función densidad es

Los parámetros cumplen la siguiente relación.

Bondad de Ajuste Chi-Cuadrado.
La Prueba de Bondad de Ajuste Chi Cuadrado es el test de bondad de ajuste más
utilizado. En general un test de bondad de ajuste se utiliza para discriminar si una
colección de datos o muestra se ajusta a una distribución teórica de una
determinada población. En otras palabras, nos dice si la muestra disponible
representa (ajusta) razonablemente los datos que uno esperaría encontrar en la
población.
 Desarrollo
Lo primero que debemos hacer es ordenar los datos por año, considerando desde
marzo a diciembre, para luego obtener el caudal anual sumando los caudales
mensuales, esto con el fin de obtener un caudal anual promedio.
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIO Qma
1950 13,9 14 9,2 7,4 12,4 13,4 23,4 43 76,1 25,4 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 25,3 307,9
1952 43,4 26 18 13,8 13,4 13,8 12,2 12,7 19,2 26,2 59 77,2 29,7 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 60 588,4
1954 126 77 37,5 20,1 14,2 14,8 12 14,9 19,7 24 57,7 61,3 29 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 26 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 19,9 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 25,8 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 23,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 28,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 30,6 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 36,5 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 29,9 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 45,5 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 16,1 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 40,5 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 26,9 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 16 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 10,8 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 20,7 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 21,2 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 23 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 64,6 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 31,8 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 34,8 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 22,8 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 23,1 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 44,5 477,4
1978 74,1 46,7 27,1 16,6 13 11,4 24,7 23,1 22,8 45 90,3 168 52,4 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 29,2 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 44,3 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 21,4 339
1982 34,4 26,2 16,4 13,8 14,1 27,2 39,9 32,8 45 50,6 88,7 182 69 571,1
1983 166 113 54,7 28,9 18,9 15,5 17,8 19,8 20 54,4 91,4 116 44,5 716,4
1984 75 49,7 26,5 17,3 12,7 10 17,3 18,6 24,9 65 89,5 118 48,7 524,5
1985 101 63,2 47,1 24,6 18,9 15,6 16,3 14,5 12,5 21 58 57,4 28,5 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 51,6 522,4
1987 101 61,5 37,7 19,2 14,5 18,4 39,8 50 35,6 55 147 176 69,1 755,7
1988 150 79 44,8 24,7 15,2 11 9,5 10,4 10,5 18,5 33,9 33 20,6 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 27,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 16,2 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 40,6 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 35,3 456,8
1993 65 38 26,2 19,8 38 31,9 16,3 14,6 17,5 27,7 45,7 54,5 29,9 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 24,7 305,2
1995 39,5 26 18,2 13,4 12 10,7 10,5 9,9 16 18,2 41,8 39,4 20 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 11,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 48,2 357,8
1998 149,1 71,1 46,8 267

Con estos datos podemos calcular el promedio, y la desviación estándar.

Tabla 1

N° 49
Qma total 18775,5
promedio 383,173469
desviación estandar 134,02057
Distribución Normal
Tabla 2
Una vez obtenido el promedio y la desviación
Caudal
Orden Año Medio Dist. Normal
entandar podemos calcular Caudales con los
(m3/s) diferentes periodos de retorno, ya que con estos
19 1968 141,9 0,00059 datos junto a los que la Distribución Normal nos
47 1996 154,8 0,00070
20 1969 207,9 0,00127
entrega podemos lograrlo.
1 1950 212,8 0,00133
41 1990 214 0,00134
18 1967 232 0,00158 N° 49
7 1956 248,5 0,00180
46 1995 255,6 0,00189 Qma total 18775,5
49 1998 267 0,00204 promedio 383,173469
27 1976 269,9 0,00208
21 1970 271,6 0,00210 [Link] 134,02057
22 1971 272,6 0,00212 Tabla 3
9 1958 287,7 0,00231
8 1957 299,5 0,00245 T (años) Prob Z Qma
45 1994 305,2 0,00251 5 0,2 0,842 496,0
2 1951 307,9 0,00254 25 0,04 1,751 617,8
10 1959 324,8 0,00271 50 0,02 2,054 658,4
26 1975 326,1 0,00272 100 0,01 2,326 695,0
40 1989 329,1 0,00274
200 0,005 2,576 728,4
6 1955 333,7 0,00278
250 0,004 2,652 738,6
3 1952 334,9 0,00279
32 1981 339 0,00282
11 1960 355 0,00291
15 1964 355,7 0,00291 En la tabla 3 se aprecian los Qma para cada uno
48 1997 357,8 0,00292 de los periodos de retorno.
17 1966 385,5 0,00298
42 1991 386,2 0,00298 En la tabla 2 se ordenaron los valores de
16 1965 389,7 0,00297
13 1962 391,4 0,00297 distribución normal de menor a mayor con el fin
44 1993 395,2 0,00296 de poder graficarlo como se puede observar en el
14 1963 408,3 0,00292 grafico 1.
25 1974 414,1 0,00290
Grafico 1
12 1961 416,9 0,00288
30 1979 429,4 0,00280
39 1988 440,5 0,00272
36 1985 450,1 0,00263
Dist. Normal
43 1992 456,8 0,00256 0.00350
28 1977 477,4 0,00232
5 1954 479,2 0,00230 0.00300
31 1980 507,6 0,00193 0.00250
23 1972 518,5 0,00179
37 1986 522,4 0,00174 0.00200
35 1984 524,5 0,00171 0.00150
29 1978 562,8 0,00121
33 1982 571,1 0,00111 0.00100
24 1973 582,4 0,00099 0.00050
4 1953 588,4 0,00092
34 1983 716,4 0,00014
0.00000
38 1987 755,7 0,00006 0 200 400 600 800
Prueba de Bondad de ajuste
Los intervalos elegidos son seis y de acuerdo a cada intervalo es que se calcula su
frecuencia relativa, variable normal, entre otros.
Tabla 4

Frec Relativa Frec Variable Frecuencia Valor

Intervalo Rango ni xi X^2
fs (xi) Acum Fs norm z F(xi) teorico
1 0 a 142 1 141,900 0,020 0,020 -1,80 0,036 0,036 0,327866264
2 142 a 284 11 236,973 0,224 0,245 -1,09 0,138 0,102 7,25428034
3 284 a 426 21 354,462 0,429 0,673 -0,21 0,415 0,278 4,028530887
4 426 a 568 11 488,109 0,224 0,898 0,78 0,783 0,368 2,742235328
5 568 a 710 3 580,633 0,061 0,959 1,47 0,930 0,146 2,431892708
6 710 a 852 2 736,050 0,041 1,000 2,63 0,996 0,066 0,473765897
49 17,25857143
Distribución Log-Normal

Tabla 5

Orden Año Caudal Medio y = ln (x) Distr. Log Normal Para la distribución Log normal se calcula
19 1968 141,9 2,152 0,091 y=ln(x) para cada año, así como también la
47 1996 154,8 2,190 0,164 distribución log normal, ordenándose de
20 1969 207,9 2,318 0,797
1 1950 212,8 2,328 0,877 menor a mayor. Esto se observa en la tabla
41 1990 214 2,330 0,897 5.
18 1967 232 2,365 1,209
7 1956 248,5 2,395 1,498 De la tabla 5, se obtiene el promedio y la
46 1995 255,6 2,408 1,618
49 1998 267 2,427 1,802 desviación estándar total de y=ln(x).
27 1976 269,9 2,431 1,847
21 1970 271,6 2,434 1,872 Promedio 2,56
22 1971 272,6 2,436 1,887 0,157
Desv Estandar
9 1958 287,7 2,459 2,095
8 1957 299,5 2,476 2,232
45 1994 305,2 2,485 2,289
Con estos valores se obtiene el caudal
2 1951 307,9 2,488 2,314 medio anual para cada periodo de retorno
10 1959 324,8 2,512 2,442
26 1975 326,1 2,513 2,450
(tabla 6)
40 1989 329,1 2,517 2,466 Tabla 6
6 1955 333,7 2,523 2,488
3 1952 334,9 2,525 2,493
32 1981 339 2,530 2,509 T (años) Prob Z P (x) Qma (m3/s)
11 1960 355 2,550 2,543 5 0,2 0,842 2,689 488,290
15 1964 355,7 2,551 2,544
48 1997 357,8 2,554 2,545
25 0,04 1,751 2,831 677,895
17 1966 385,5 2,586 2,501 50 0,02 2,054 2,879 756,247
42 1991 386,2 2,587 2,499
100 0,01 2,326 2,921 834,431
16 1965 389,7 2,591 2,486
13 1962 391,4 2,593 2,480
200 0,005 2,576 2,960 913,047
44 1993 395,2 2,597 2,464 250 0,004 2,652 2,972 938,519
14 1963 408,3 2,611 2,397
25 1974 414,1 2,617 2,363
De la distribución normal de la tabla 5 se
12 1961 416,9 2,620 2,346
30 1979 429,4 2,633 2,262 puede obtener el grafico 2.
39 1988 440,5 2,644 2,180
36 1985 450,1 2,653 2,105 Grafico 2
43 1992 456,8 2,660 2,051
28 1977 477,4 2,679 1,879
5 1954 479,2 2,681 1,864 Distr. Log Normal
31 1980 507,6 2,706 1,622
23 1972 518,5 2,715 1,532 3.000
37 1986 522,4 2,718 1,499 2.500
35 1984 524,5 2,720 1,482 2.000
29 1978 562,8 2,750 1,187 1.500
33 1982 571,1 2,757 1,128
1.000
24 1973 582,4 2,765 1,051
4 1953 588,4 2,770 1,012 0.500
34 1983 716,4 2,855 0,416 0.000
38 1987 755,7 2,878 0,310 0 200 400 600 800

Distr. Log Normal

Prueba de Bondad de ajuste
Los intervalos elegidos son seis y de acuerdo a cada intervalo es que se calcula su
frecuencia relativa, variable normal, entre otros.
Tabla 7

Frec Relativa fs Frec Acom Fs Valor teorico

Intervalo Rango ni xi m3/s log x Variable norm z Frecuencia F(xi) X^2
(xi) (xi) p(xi)

1 0 a 142 1 141,9 2,2 0,02 0,02 -2,58 0,00 0,00 2,40202403

2 142 a 28411 236,9727273 2,4 0,22 0,24 -1,16 0,12 0,12 4,73522914
3 284 a 42621 354,4619048 2,5 0,43 0,67 -0,05 0,48 0,36 0,66075122
4 426 a 56811 488,1090909 2,7 0,22 0,90 0,64 0,80 0,32 1,34781586
5 568 a 7103 580,6333333 2,8 0,06 0,96 1,32 0,91 0,11 0,96385422
6 710 a 8522 736,05 2,9 0,04 1,00 1,90 0,98 0,07 0,57158198
49 10,6812565
Distribución de Gumbel
Así como en las anteriores distribuciones, también se ordenan de menor a mayor
los datos de caudales medios. Como se muestra en la tabla 8.
Tabla 8

Orden Año Caudal Medio De la tabla 9, según nuestra cantidad de muestra es que
19 1968 141,9
47 1996 154,8
podemos calcular ỹ y σn, iterando con 49 muestras entre
20 1969 207,9 45 y 50.
1 1950 212,8
41 1990 214
Tabla 9
18 1967 232
7 1956 248,5
46 1995 255,6 N ỹ σn
49 1998 267
27 1976 269,9 0 0 0,00
21 1970 271,6
22 1971 272,6 10 0,4950 0,9500
9 1958 287,7
15 0,5130 1,0210
8 1957 299,5
45 1994 305,2 20 0,5240 1,0630
2 1951 307,9
10 1959 324,8
25 0,5310 1,0910
26 1975 326,1 30 0,5360 1,1120
40 1989 329,1
6 1955 333,7 35 0,5400 1,1280
3 1952 334,9
32 1981 339
40 0,5440 1,1410
11 1960 355 45 0,5460 1,1520
15 1964 355,7
48 1997 357,8 49 0,5484 1,1592
17 1966 385,5
42 1991 386,2
50 0,5490 1,1610
16 1965 389,7
13 1962 391,4 Con ỹ y σn, se calcula α y μ. El promedio y la desviación
44 1993 395,2
14 1963 408,3
estándar es la misma con la que hemos trabajo en las
25 1974 414,1 distribuciones anteriores.
12 1961 416,9
30 1979 429,4 Tabla 10
39 1988 440,5
36 1985 450,1
43 1992 456,8 Promedio 383,173469
28 1977 477,4
5 1954 479,2
Desv. Estand 134,02057
31 1980 507,6 ỹ 0,5484
23 1972 518,5
37 1986 522,4 σn 1,1592
35 1984 524,5
29 1978 562,8 α 0,00864942
33 1982 571,1
24 1973 582,4
μ 319,770363
4 1953 588,4
34 1983 716,4
38 1987 755,7
Una vez obtenido los datos que se muestran en la tabla 10, se pueden calcular los
caudales medios anuales, con los diferentes periodos de retorno.
Con Qma =(Y- ỹ)/ σn * devesta + promedio
Tabla 11

T (años) Prob Y Qma (m3/s)

5 0,2 1,50 493,19
25 0,04 3,20 689,57
50 0,02 3,90 770,89
100 0,01 4,6 851,62
200 0,005 5,30 932,04
250 0,004 5,52 957,90

De la tabla 11 se obtiene el grafico de distribución de Gumbel (grafico 3).

Grafico 3

Distribución Gumbel
1200.00
1000.00
800.00
600.00
400.00
200.00
0.00
0 50 100 150 200 250 300

Prueba de Bondad de ajuste

Los intervalos elegidos son seis y de acuerdo a cada intervalo es que se calcula su
frecuencia relativa, variable normal, entre otros.
Tabla 12

Frec Relativa Frec Acom Frecuencia F Valor teorico

Intervalo Rango ni m3/s xi landa landa*sigma landa X^2
Fs (xi) Fs(xi) (xi) p (xi)
1 0 a 142 1 141,9 0,463 0,822 -1,02 0,02040816 0,02040816 0,06246213 0,06246213 1,38737312
2 142 a 284 11 236,972727 0,511 1,013 -0,59 0,2244898 0,24489796 0,16464092 0,10217879 7,17408367
3 284 a 426 21 354,461905 0,526 1,071 0,3 0,42857143 0,67346939 0,47672369 0,31208277 2,13055905
4 426 a 568 11 488,109091 0,511 0,964 1,27 0,2244898 0,89795918 0,75515548 0,27843179 0,51207219
5 568 a 710 3 580,633333 0,47 0,851 1,72 0,06122449 0,95918367 0,83605059 0,08089512 0,23437443
6 710 a 852 2 736,05 0,466 0,836 2,67 0,04081633 1 0,9330913 0,09704071 1,59621534
49 13,0346778
Distribución Gamma
Tabla 13
Nuevamente se ordena de menor a mayor
Orden Año
Caudal Medio Distribución con respecto al caudal medio y se calcula su
(m3/s) Gamma
distribución, en este caso gamma.
19 1968 141,9 0,00041
47 1996 154,8 0,00058 Promedio 383,173
20 1969 207,9 0,00154
1 1950 212,8 0,00164 Desv Estan 134,021
41 1990 214 0,00167
18 1967 232 0,00203 Con el promedio y desviación estándar,
7 1956 248,5 0,00233 calculamos α y β.
46 1995 255,6 0,00245
49 1998 267 0,00263 α 8,17425
27 1976 269,9 0,00267
21 1970 271,6 0,00270 β 46,8757
22 1971 272,6 0,00271
9 1958 287,7 0,00289 Teniendo esos datos se encuentran los Qma
8 1957 299,5 0,00300 para cada periodo de retorno. Se muestra en
45 1994 305,2 0,00304
2 1951 307,9 0,00306 la tabla
10 1959 324,8 0,00313
Tabla 14
26 1975 326,1 0,00313
40 1989 329,1 0,00314
6 1955 333,7 0,00314
3 1952 334,9 0,00314
T (años) Prob Qma (m3/s)
32 1981 339 0,00314
11 1960 355 0,00311 5 0,2 489,06
15 1964 355,7 0,00310 25 0,04 646,72
48 1997 357,8 0,00310
17 1966 385,5 0,00293
50 0,02 705,69
42 1991 386,2 0,00292 100 0,01 761,55
16 1965 389,7 0,00289
200 0,005 815,04
13 1962 391,4 0,00288
44 1993 395,2 0,00285 250 0,004 831,84
14 1963 408,3 0,00272
25 1974 414,1 0,00266 De la tabla 12 obtenemos el grafico 4.
12 1961 416,9 0,00263
30 1979 429,4 0,00249 Grafico 4
39 1988 440,5 0,00236
36 1985 450,1 0,00224
43 1992 456,8 0,00216 Distribución Gamma
28 1977 477,4 0,00191 0.00350
5 1954 479,2 0,00189
31 1980 507,6 0,00156
0.00300
23 1972 518,5 0,00144 0.00250
37 1986 522,4 0,00140 0.00200
35 1984 524,5 0,00137
0.00150
29 1978 562,8 0,00101
33 1982 571,1 0,00094 0.00100
24 1973 582,4 0,00085 0.00050
4 1953 588,4 0,00080 0.00000
34 1983 716,4 0,00021
0 200 400 600 800
38 1987 755,7 0,00014
Prueba de Bondad de ajuste
Los intervalos elegidos son seis y de acuerdo a cada intervalo es que se calcula su
frecuencia relativa, variable normal, entre otros.
Tabla 15

Intervalo Rango ni Xi anual Frec Relativa Fs Frec Acum Frecuencia F Valor teorico p X^2
1 0 a 142 1 141,9 0,020408163 0,020408163 0,010338586 0,010338586 0,480570837
2 142 a 284 11 236,9727273 0,224489796 0,244897959 0,124853293 0,114514707 5,175156209
3 284 a 426 21 354,4619048 0,428571429 0,673469388 0,459305166 0,334451874 1,297843043
4 426 a 568 11 488,1090909 0,224489796 0,897959184 0,798311496 0,339006329 1,895503802
5 568 a 710 3 580,6333333 0,06122449 0,959183673 0,917625764 0,119314268 1,38580824
6 710 a 852 2 736,05 0,040816327 1 0,986221221 0,068595457 0,551236572
49 10,7861187
 Resultados y conclusiones
Para rechazar si las distribuciones se adaptan a la muestra tenemos que ver el valor
que nos dio la prueba de bondad de ajuste X^2 para cada una de las distribuciones
y compararla con Xv^2.
Para calcular Xv^2, lo hacemos con el comando de Excel [Link](A;B)
Con A= (1-nivel de confianza) y B= grados de libertad.
Vamos a suponer un nivel de confianza de 95% por lo tanto A=(1-0.95)= 0.05
B=(Columnas-1)(filas-1) = (2-1)(6-1)=(1)(5)=5.
Los valores de Xva^2 son los calculados en las tablas
Tabla 16

Xv^2 Xva^2 Xv^2 < Xva^2

Normal 17,2585714 11,0704977 Se rechaza
Log-normal 10,6812565 11,0704977 No se rechaza
Gumbel 13,0346778 11,0704977 Se rechaza
Gamma 10,7861187 11,0704977 No se rechaza

 Con esto podemos decir que se rechaza la distribución normal y la de

Gumbel.

 La distribución que más se acerca es la de log-normal, ya que, si dos o más

distribuciones son aceptadas, se elige como más probable la que tenga
menor valor de Xva^2. Por consiguiente, los valores caudales que más se
acercan para los periodos de retorno 5, 25, 50, 100, 200 y 250.

T (años) Prob Z P (x) Qma (m3/s)

5 0,2 0,842 2,689 488,290
25 0,04 1,751 2,831 677,895
50 0,02 2,054 2,879 756,247
100 0,01 2,326 2,921 834,431
200 0,005 2,576 2,960 913,047
250 0,004 2,652 2,972 938,519