Considerando dos años:

Año Ene. Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
1968 9.92 12.50 13.84 8.05 1.84 1.21 0.80 0.64 0.65 0.81 2.21 4.50
1969 8.00 17.00 25.00 14.45 2.75 1.10 0.82 0.71 0.56 1.86 6.42 10.67

Q V Vacum
t
m3/s MMC MMC
1 9.92 26.57 26.57
2 12.50 30.24 56.81
3 13.84 37.07 93.88
4 8.05 20.87 114.75
5 1.84 4.93 119.68
6 1.21 3.14 122.82
7 0.80 2.14 124.96
8 0.64 1.71 126.67
9 0.65 1.68 128.35
10 0.81 2.17 130.52
11 2.21 5.73 136.25
12 4.50 12.05 148.30
13 8.00 21.43 169.73
14 17.00 41.13 210.86
15 25.00 66.96 277.82
16 14.45 37.45 315.27
17 2.75 7.37 322.64
18 1.10 2.85 325.49
19 0.82 2.2 327.69
20 0.71 1.9 329.59
21 0.56 1.45 331.04
22 1.86 4.98 336.02
23 6.42 16.64 352.66
24 10.67 28.58 381.24

Graficando, tenemos la Curva Masa:

 450
400
Curva Masa
Volumen Acumulado (MMC)
350
300
250
200
150
100
50
0
0 4 8 12 16 20 24

t (meses)

Para calcular el volumen útil de almacenamiento para

regular un caudal QR, se debe representar este caudal
mediante una recta calculando el volumen “V” para un
tiempo conocido “t”, por ejemplo para 8 meses:
V = QR * 8*3600*24*30/106 MMC

QR

12 t 20

También se puede realizar este cálculo de forma analítica.

Ej.: Para el periodo crítico de dos años mostrado, se pide
determinar el volumen útil que se debe almacenar para
regular un caudal de 12 m3/s.

Año Ene. Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
1965 25.30 28.70 15.86 19.84 9.36 4.80 4.50 4.25 3.78 3.55 13.40 16.20
1966 21.62 20.43 33.28 11.15 5.40 3.61 2.91 3.48 3.31 3.27 18.74 25.86
Realizamos el siguiente cuadro:
A B
Q V V.Reg.
T
m3/s MMC MMC A-B
1 25.30 67.76 32.14 35.62
2 28.70 69.43 29.03 40.40
3 15.86 42.48 32.14 10.34
4 19.84 51.43 31.10 20.32 106.68
5 9.36 25.07 32.14 -7.07
6 4.80 12.44 31.10 -18.66
7 4.50 12.05 32.14 -20.09
8 4.25 11.38 32.14 -20.76
9 3.78 9.80 31.10 -21.31
10 3.55 9.51 32.14 -22.63 -110.52
11 13.40 34.73 31.10 3.63
12 16.20 43.39 32.14 11.25
13 21.62 57.91 32.14 25.77
14 20.43 49.42 29.03 20.39
15 33.28 89.14 32.14 57.00 103.16
16 11.15 28.90 31.10 -2.20
17 5.40 14.46 32.14 -17.68
18 3.61 9.36 31.10 -21.75
19 2.91 7.79 32.14 -24.35
20 3.48 9.32 32.14 -22.82
21 3.31 8.58 31.10 -22.52
22 3.27 8.76 32.14 -23.38 -134.70
23 18.74 48.57 31.10 17.47
24 25.86 69.26 32.14 37.12

∑Primer año, diferencias positivas: 121,56 MMC

∑Segundo año, diferencias positivas: 157,75 MMC

El volumen útil de almacenamiento necesario para regular

el caudal de 12 m3/s es: 134,70 MMC.

ESTUDIO DE CAUDALES DE AVENIDAS

Los caudales de avenidas son caudales excepcionales que
se presentan en un río en un determinado periodo de tiempo.
PERIODO DE RETORNO: T
Es el número de años que debe transcurrir en promedio para
que se presente un caudal igual o mayor a una magnitud
dada. Decir: “El caudal de avenida en un río para un periodo
de retorno de 50 años es 180 m3/s”, es equivalente a decir
que en promedio un caudal de esa magnitud o mayor será
transportado por el río una vez cada 50 años.

𝑇 = 1⁄𝑃
Donde: T: Periodo de retorno, en años
P: Probabilidad de que un caudal transportado
por el río sea mayor o igual a QT

RIESGO: R
Se define como riesgo a la probabilidad que un caudal
determinado ocurra una vez en “n” años sucesivos.

Sea “p” la probabilidad de que ocurra un caudal

determinado en un año cualquiera, entonces:
(1-p): Representa la probabilidad de que dicho caudal no
ocurra ese año.
n: Vida útil de la obra.

La probabilidad de que dicho caudal no ocurra en “n” años

sucesivos, será: (1-p)n
Por lo tanto: 𝑅 = 1 − (1 − 𝑝)𝑛

1 𝑛
También: 𝑅 = 1 − (1 − 𝑇)
Ejemplo: Determinar el riesgo o probabilidad de falla de
una estructura hidráulica, que tendrá una vida útil de 50
años, si se diseña considerando un caudal de avenida para
los siguientes periodos de retorno: 25 años, 50 años, 100
años, 1000 años y 10 000 años.

Aplicando la ecuación del Riesgo para cada periodo de

retorno hallamos:

T (años) 25 50 100 1000 10000

R (%) 87,01 63,58 39,5 4,88 0,50

Presentamos a continuación algunos periodos de retorno

recomendados para el cálculo del caudal de avenidas en el
diseño de estructuras hidráulicas:
 FACTOR DE CORRECCIÓN DE FULLER: Ff
Es un factor que permite estimar el caudal máximo
instantáneo a partir del caudal máximo anual.
2,66
𝐹𝑓 = 1 + 𝑄𝑀𝐼 = 𝐹𝑓 . 𝑄𝑀𝐴 Donde:
𝐴 0,30

A: Área de la cuenca, km2

QMI: Caudal máximo instantáneo
QMA: Caudal máximo anual.

Para determinar los caudales de avenidas en una estación se

deberá establecer el modelo probabilístico que mejor
represente al comportamiento aleatorio de los caudales del
río, para ello se deberán realizar pruebas denominadas de
Bondad de Ajuste, considerando un conjunto de funciones
probabilísticas. La función que mejor se ajuste al
comportamiento aleatorio de los caudales será la que nos
permita calcular la ocurrencia de un caudal máximo para un
determinado periodo de retorno.

PRUEBAS DE BONDAD DE AJUSTE

Estudiaremos las pruebas de bondad de ajuste siguientes:
* Chi Cuadrado * Smirnov - Kolmogorov
Las funciones de probabilidad estudiaremos para el cálculo
de los caudales máximos son:

* Distribución Normal * Distribución Log Normal

* Distribución Pearson III *Distribución LogPearson III
* Distribución Gumbel
 PRUEBA DE BONDAD DE AJUSTE
CHI CUADRADO: χ2

La expresión general de Chi Cuadrado viene dada por:

𝑘
(𝜃𝑖 − 𝑒𝑖 )2
𝜒𝐶2 = ∑
𝑒𝑖
𝑖=1

Donde:
χ2C: Valor calculado de Chi cuadrado a partir del registro de
caudales máximos instantáneos.
θi: Número de valores observados en intervalo de clase “i”
ei: Número de valores esperados en intervalo de clase “i”
k: Número de intervalos de clase

El valor teórico de Chi Cuadrado (χ2t) se obtiene de la tabla

de distribución del mismo nombre considerando:
- Nivel de significación: α, por lo general α = 5%
- Grados de libertad: k - h – 1. Donde:

h = 2 Distribución Normal, Log Normal y Gumbel

h = 3 Distribución Pearson III, Log Pearson III

Criterio de Decisión:

Si: χ2C ≤ χ2t: El registro de caudales se ajusta a la

distribución de probabilidades para el nivel
de significación considerado.
Si: χ2C > χ2t: El registro de caudales no se ajusta a la
distribución de probabilidades para el nivel
de significación considerado. Se debe probar
con otra distribución.

Este método requiere agrupar los datos en intervalos de

clase, su aplicación es más adecuada para ajustes a la
distribución Normal puesto que ha sido desarrollado a base
a datos normales e independientes. En la práctica se utiliza
para verificar el ajuste a cualquier distribución de
probabilidades.

Ejemplo: Para el siguiente registro de caudales máximos

anuales, en m3/s, registrado en una estación cuya cuenca es
250 km2, se pide aplicar la prueba Chi Cuadrado y
determinar si la información se ajusta a la Distribución
Normal para un nivel de significación del 5%.
Año QMA Año QMA Año QMA Año QMA Año QMA
1951 156.49 1961 71.59 1971 145.96 1981 115.63 1991 182.65
1952 146.89 1962 100.33 1972 218.08 1982 144.57 1992 93.31
1953 106.82 1963 182.38 1973 249.93 1983 173.91 1993 156.29
1954 121.46 1964 80.73 1974 137.95 1984 76.23 1994 108.61
1955 136.62 1965 154.64 1975 188.34 1985 228.01 1995 149.4
1956 216.82 1966 137.22 1976 96.49 1986 157.55 1996 136.56
1957 93.97 1967 236.62 1977 218.08 1987 121.85 1997 92.12
1958 170.86 1968 107.88 1978 174.04 1988 143.71 1998 125.03
1959 73.84 1969 74.3 1979 154.83 1989 104.64 1999 115.63
1960 133.25 1970 139.8 1980 130.6 1990 146.16 2000 150.13

Aplicando el Factor de Fuller al registro de caudales

máximos anuales, determinamos los caudales máximos
instantáneo, para un área de 250 km2:
2,66
𝐹𝑓 = 1 + = 1,51 𝑄𝑀𝐼 = 1,51 . 𝑄𝑀𝐴
2500,30
Luego los caudales máximos instantáneos serán los
siguientes:

Año QMA Año QMA Año QMA Año QMA Año QMA
1951 236.30 1961 108.10 1971 220.40 1981 174.60 1991 275.80
1952 221.80 1962 151.50 1972 329.30 1982 218.30 1992 140.90
1953 161.30 1963 275.39 1973 377.39 1983 262.60 1993 236.00
1954 183.40 1964 121.90 1974 208.30 1984 115.11 1994 164.00
1955 206.30 1965 233.51 1975 284.39 1985 344.30 1995 225.59
1956 327.40 1966 207.20 1976 145.70 1986 237.90 1996 206.21
1957 141.89 1967 357.30 1977 329.30 1987 183.99 1997 139.10
1958 258.00 1968 162.90 1978 262.80 1988 217.00 1998 188.80
1959 111.50 1969 112.19 1979 233.79 1989 158.01 1999 174.60
1960 201.21 1970 211.10 1980 197.21 1990 220.70 2000 226.70

Nuestro objetivo es calcular χ2C y compararlo con χ2t

𝑘
( )2
𝜃𝑖 − 𝑒𝑖
𝜒𝐶2 = ∑
𝑒𝑖
𝑖=1

 Cálculo del número de intervalos de clase: k

Según Yevjevich: 𝑘 = 1 + 1,33 𝐿𝑛 𝑁

N: Número de datos de la muestra, en este caso N = 50

𝑘 = 1 + 1,33 𝐿𝑛 50 = 6,20 ≈ 6

 Cálculo de la amplitud del intervalo de clase: ΔQ

𝑄𝑚𝑎𝑥 − 𝑄𝑚𝑖𝑛
Δ𝑄 =
𝑘−1
En este caso:
Qmax = 377,39 m3/s Qmin = 108,10 m3/s

377,39 − 108,10
Δ𝑄 = = 53,86 𝑚3 /𝑠
6−1
  Cálculo del primer del intervalo de clase:
Límite Inferior: Qmin – ΔQ/2; Límite Superior: Qmin + ΔQ/2
Reemplazando:
Límite Inferior: 108,10 – 53,86/2 = 81,17
Límite Superior: 108,10 + 53,86/2 = 135,03
Podemos ahora calcular los 6 intervalos de clase, las marcas
de clase, como el promedio de los límites superior e inferior
de cada intervalo y el número de valores observados en cada
intervalo de clase: θ
k Linf Lsup Mclase θ
1 81.17 135.03 108.10 5
2 135.03 188.89 161.96 14
3 188.89 242.75 215.82 19
4 242.75 296.61 269.68 6
5 296.61 350.47 323.54 4
6 350.47 404.33 377.40 2

Para hallar el número de valores esperados en cada intervalo

de clase (e) debemos aplicar la siguiente relación:
𝑒𝑖 = 𝑁 [𝐹(𝑆𝑖 ) − 𝐹(𝐼𝑖 ) ] Donde:
F(Si): Función de distribución de probabilidad en el límite
superior del intervalo de clase “i”.
F(Ii): Función de distribución de probabilidad en el límite
inferior del intervalo de clase “i”.
Para determinar la función de distribución Normal de
probabilidades debemos primero determinar la variable
estandarizada (Z):
𝑄𝑖 − 𝑄𝑚
𝑍𝑖 =
𝑆
 Zi: Variable estandarizada para el caudal Qi
Qm: Promedio de los valores agrupados
S: Desviación estándar de los valores agrupados.

∑𝑘
𝑖=1 𝑄𝑖 . 𝜃𝑖 ∑𝑘 2
𝑖=1(𝑄𝑖 − 𝑄𝑚 ) 𝜃𝑖
𝑄𝑚 = 𝑆= √
𝑁 𝑁−1

Donde: Qi: Marca de clase de cada intervalo

Reemplazando valores obtenemos:

Qm = 211,51 S = 66,05

Los cálculos se muestran en la siguiente tabla:

Linf Lsup Mclase θ Zinf Zsup
81.17 135.03 108.10 5 -1.9735 -1.1580
135.03 188.89 161.96 14 -1.1580 -0.3425
188.89 242.75 215.82 19 -0.3425 0.4730
242.75 296.61 269.68 6 0.4730 1.2885
296.61 350.47 323.54 4 1.2885 2.1040
350.47 404.33 377.40 2 2.1040 2.9195

De la tabla de distribución normal hallamos los valores de

F(Z) para los límites inferior y superior de cada intervalo de
clase y completamos los cálculos:

Linf Lsup Mclase θ Zinf Zsup F(Zinf) F(Zsup) e (θi - ei)2/ei

81.17 135.03 108.10 5 -1.9735 -1.1580 0.0242 0.1234 4.96 0.00
135.03 188.89 161.96 14 -1.1580 -0.3425 0.1234 0.366 12.13 0.29
188.89 242.75 215.82 19 -0.3425 0.4730 0.366 0.6819 15.80 0.65
242.75 296.61 269.68 6 0.4730 1.2885 0.6819 0.9012 10.97 2.25
296.61 350.47 323.54 4 1.2885 2.1040 0.9012 0.9823 4.06 0.00
350.47 404.33 377.40 2 2.1040 2.9195 0.9823 0.9982 0.80 1.80

Sumando la última columna hallamos: χ2C = 4,99

De la tabla de distribución Chi-Cuadrado con:
 G. Libertad: k – h – 1: 6 – 2 – 1 = 3
 Nivel de significación: α = 5%
Obtenemos: χ2t = 7,81
Comparando vemos que: χ2C < χ2t es decir:

“El registro de caudales se ajusta a la distribución

Normal de probabilidades con un nivel de
significación del 5%”