Curso de Hidrología.

Docente Eduardo Luis Flores Quispe

ANALISIS DE FRECUENCIA DE TORMENTAS

1. Relaciones intensidad – duración – frecuencia

CHOW et al. (1994), Uno de los primeros pasos que debe seguirse en muchos proyectos de
diseño hidrológico, como el diseño de un drenaje urbano, es la determinación del evento o los
eventos de lluvia que deben usarse. La forma más común de hacerlo es utilizar una tormenta
de diseño o un evento que involucre una relación entre la intensidad de lluvia (o profundidad),
la duración y las frecuencias o períodos de retorno apropiados para la obra y el sitio. En
muchos casos existen curvas estándar de intensidad-duración-frecuencia (IDF) disponibles para
el sitio, luego no hay que llevar a cabo este análisis. Sin embargo es conveniente entender el
procedimiento utilizado para desarrollar estas relaciones. Usualmente los datos se presentan
en forma gráfica, con la duración en el eje horizontal y la intensidad en el eje vertical,
mostrando una serie de curvas, para cada uno de los períodos de retorno de diseño.

2. Ejemplo de calculo

Dados los datos de intensidades máximas en mm/hr, para diferentes duraciones (5, 10, 15, 20,
30, 40, 60, 80, 100 y 120 minutos) para el período 1956-1970, se pide:

- Dibujar las curvas de intensidad – duración – período de retorno (frecuencia)

- Realizar un ajuste utilizando la fórmula utilizada en USA

Dónde:
a, b, K = parámetros de regresión
T = período de retorno en años
Imáx = intensidad máxima en mm/hr
D = duración en minutos

Determinar el valor de la intensidad máxima de una lluvia con una duración de 45 minutos que
tiene un período de retorno de 10 años.

Los datos de intensidades máximas se muestran en el siguiente cuadro.

-1-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

Cuadro 5. Intensidades máximas de lluvia en mm/hr

Minutos
Fecha 5 10 15 20 30 40 60 80 100 120
15-jun-56 120.0 81.5 54.8 41.7
20-ago-56 36.0 34.8 25.7 21.5 18.2 16.8
23-feb-57 51.0 36.0 28.8 20.4 16.1 12.7 12.3 10.9
27-jun-57 84.0
26-jul-57 9.5
16-jun-58 20.7
01-ago-58 11.6 10.1
11-sep-58 23.4 16.5 13.7 11.6
15-sep-58 116.4 40.4 30.3
07-jun-59 48.1 38.0 32.2 27.6
26-sep-59 240.0 162.0 128.0 111.0 84.0 67.5
18-ago-60 122.4
09-sep-60 42.6 31.4 18.0 17.3 15.5
20-sep-60 23.0
28-oct-60 122.4 69.0 52.0 30.0
23-may-61 15.8 13.9
30-jul-61 48.0
31-jul-61 25.6 21.9 18.2 11.8 11.7 10.3
10-sep-61 30.0
11-abr-62 120.0 120.0 84.0 82.8 60.0 48.9 32.7 24.5 19.6 16.4
10-ago-63 10.1
26-nov-63 116.4 58.2 51.6 44.7 36.4 29.0 19.3 14.5 11.6
15-jun-64 36.0 27.0 26.0 24.0 16.6 12.8
22-ago-64 9.5 8.1 7.1 6.3
06-oct-65 78.0 55.2 48.0 48.0 41.8 34.2 29.5 22.9 18.7 15.9
18-oct-66 78.0 60.0 52.0 47.4 46.0 40.5 29.2 22.9 18.9 16.3
17-jun-67 54.0 34.0 23.4 18.8 14.0 11.6 11.6 10.6
22-jun-67 45.0 28.5 18.8
29-jun-68 84.0 48.0 32.0 24.0 16.0 12.0 8.0 6.0
02-jul-68 6.0
12-oct-68 5.3 4.7
07-abr-69 54.0 54.0 45.2 41.7 31.8 25.1 17.8 13.4
19-ago-69 10.0
29-ago-69 12.0 10.0
25-jun-70 120.0 84.0 72.0 57.6 41.0 31.1 22.5 17.0 13.7
04-ago-70 12.0

El primer paso para la solución es ordenar en orden descendente las series para cada duración,
luego se calcula la probabilidad de excedencia que es inversa al período de retorno. En este
caso se utiliza la fórmula de Weibull aunque existen otras.

-2-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

( )

Dónde:

P(X>x) = probabilidad que la variable aleatoria X (en este caso intensidad máxima) sea mayor a
un valor x.

m = orden de la posición descendente.

n = número total de datos analizados.

El período de retorno (T) es la inversa de la probabilidad de excedencia

( )

En el siguiente cuadro se muestra los valores ordenados descendentemente y los cálculos de

probabilidad y período de retorno.

Cuadro 6. Datos ordenados descendentemente y cálculo del período de retorno

Minutos
m P(X>x) T (años) 5 10 15 20 30 40 60 80 100 120
1 0.06 17.00 240.0 162.0 128.0 111.0 84.0 67.5 48.1 38.0 32.2 27.6
2 0.12 8.50 122.4 120.0 84.0 82.8 60.0 48.9 32.7 24.5 19.6 16.8
3 0.18 5.67 122.4 84.0 72.0 57.6 46.0 40.5 29.5 22.9 18.9 16.4
4 0.24 4.25 120.0 81.5 54.8 48.0 41.8 34.8 29.2 22.9 18.7 16.3
5 0.29 3.40 120.0 69.0 52.0 47.4 41.0 34.2 25.7 21.5 18.2 15.9
6 0.35 2.83 120.0 60.0 52.0 44.7 36.4 31.1 23.0 18.0 17.3 15.5
7 0.41 2.43 116.4 58.2 51.6 42.6 36.0 30.0 22.5 17.0 13.7 12.0
8 0.47 2.13 116.4 55.2 48.0 41.7 31.8 29.0 19.3 14.5 12.0 10.6
9 0.53 1.89 84.0 54.0 45.2 41.7 31.4 25.1 17.8 13.7 11.7 10.3
10 0.59 1.70 84.0 51.0 36.0 28.8 23.4 20.7 16.5 13.4 11.6 10.1
11 0.65 1.55 78.0 48.0 34.0 28.5 23.4 18.8 14.0 12.3 11.6 10.1
12 0.71 1.42 78.0 45.0 32.0 24.0 20.4 18.8 13.9 11.8 11.6 10.0
13 0.76 1.31 54.0 40.4 30.3 24.0 18.2 16.1 12.7 11.6 11.6 10.0
14 0.82 1.21 54.0 30.0 26.0 21.9 16.6 15.8 9.5 8.1 10.9 9.5
15 0.88 1.13 48.0 27.0 25.6 16.0 12.8 8.0 6.0 7.1 6.3
16 0.94 1.06 36.0 12.0 6.0 5.3 4.7

Para determinar los parámetros de ajuste a la fórmula utilizada en USA se debe linealizar esta
aplicando logaritmos.

( ) ( ) ( ) ( )

Así se transforma en un modelo de regresión lineal múltiple

-3-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

Los parámetros a y b se obtiene directamente pero K se obtiene así

( )

Si se aplica el logaritmo en base 10.

Para el proceso de regresión se realiza el ordenamiento en el siguiente cuadro

-4-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

Cuadro 6. Transformación para realizar la regresión múltiple

N° T D Imáx log(T) log(D) log(Imáx) N° T D Imáx log(T) log(D) log(Imáx) N° T D Imáx log(T) log(D) log(Imáx)
1 17.00 5 240 1.23 0.70 2.38 55 2.83 30 36.4 0.45 1.48 1.56 108 1.55 80 12.3 0.19 1.90 1.09
2 17.00 10 162 1.23 1.00 2.21 56 2.83 40 31.1 0.45 1.60 1.49 109 1.55 100 11.6 0.19 2.00 1.06
3 17.00 15 128 1.23 1.18 2.11 57 2.83 60 23 0.45 1.78 1.36 110 1.55 120 10.1 0.19 2.08 1.00
4 17.00 20 111 1.23 1.30 2.05 58 2.83 80 18 0.45 1.90 1.26 111 1.42 5 78 0.15 0.70 1.89
5 17.00 30 84 1.23 1.48 1.92 59 2.83 100 17.3 0.45 2.00 1.24 112 1.42 10 45 0.15 1.00 1.65
6 17.00 40 67.5 1.23 1.60 1.83 60 2.83 120 15.5 0.45 2.08 1.19 113 1.42 15 32 0.15 1.18 1.51
7 17.00 60 48.1 1.23 1.78 1.68 61 2.43 5 116.4 0.39 0.70 2.07 114 1.42 20 24 0.15 1.30 1.38
8 17.00 80 38 1.23 1.90 1.58 62 2.43 10 58.2 0.39 1.00 1.76 115 1.42 30 20.4 0.15 1.48 1.31
9 17.00 100 32.2 1.23 2.00 1.51 63 2.43 15 51.6 0.39 1.18 1.71 116 1.42 40 18.8 0.15 1.60 1.27
10 17.00 120 27.6 1.23 2.08 1.44 64 2.43 20 42.6 0.39 1.30 1.63 117 1.42 60 13.9 0.15 1.78 1.14
11 8.50 5 122.4 0.93 0.70 2.09 65 2.43 30 36 0.39 1.48 1.56 118 1.42 80 11.8 0.15 1.90 1.07
12 8.50 10 120 0.93 1.00 2.08 66 2.43 40 30 0.39 1.60 1.48 119 1.42 100 11.6 0.15 2.00 1.06
13 8.50 15 84 0.93 1.18 1.92 67 2.43 60 22.5 0.39 1.78 1.35 120 1.42 120 10 0.15 2.08 1.00
14 8.50 20 82.8 0.93 1.30 1.92 68 2.43 80 17 0.39 1.90 1.23 121 1.31 5 54 0.12 0.70 1.73
15 8.50 30 60 0.93 1.48 1.78 69 2.43 100 13.7 0.39 2.00 1.14 122 1.31 10 40.4 0.12 1.00 1.61
16 8.50 40 48.9 0.93 1.60 1.69 70 2.43 120 12 0.39 2.08 1.08 123 1.31 15 30.3 0.12 1.18 1.48
17 8.50 60 32.7 0.93 1.78 1.51 71 2.13 5 116.4 0.33 0.70 2.07 124 1.31 20 24 0.12 1.30 1.38
18 8.50 80 24.5 0.93 1.90 1.39 72 2.13 10 55.2 0.33 1.00 1.74 125 1.31 30 18.2 0.12 1.48 1.26
19 8.50 100 19.6 0.93 2.00 1.29 73 2.13 15 48 0.33 1.18 1.68 126 1.31 40 16.1 0.12 1.60 1.21
20 8.50 120 16.8 0.93 2.08 1.23 74 2.13 20 41.7 0.33 1.30 1.62 127 1.31 60 12.7 0.12 1.78 1.10
21 5.67 5 122.4 0.75 0.70 2.09 75 2.13 30 31.8 0.33 1.48 1.50 128 1.31 80 11.6 0.12 1.90 1.06
22 5.67 10 84 0.75 1.00 1.92 76 2.13 40 29 0.33 1.60 1.46 129 1.31 100 11.6 0.12 2.00 1.06
23 5.67 15 72 0.75 1.18 1.86 77 2.13 60 19.3 0.33 1.78 1.29 130 1.31 120 10 0.12 2.08 1.00
24 5.67 20 57.6 0.75 1.30 1.76 78 2.13 80 14.5 0.33 1.90 1.16 131 1.21 5 54 0.08 0.70 1.73
25 5.67 30 46 0.75 1.48 1.66 79 2.13 100 12 0.33 2.00 1.08 132 1.21 10 30 0.08 1.00 1.48
26 5.67 40 40.5 0.75 1.60 1.61 80 2.13 120 10.6 0.33 2.08 1.03 133 1.21 15 26 0.08 1.18 1.41
27 5.67 60 29.5 0.75 1.78 1.47 81 1.89 5 84 0.28 0.70 1.92 134 1.21 20 21.9 0.08 1.30 1.34

-5-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

28 5.67 80 22.9 0.75 1.90 1.36 82 1.89 10 54 0.28 1.00 1.73 135 1.21 30 16.6 0.08 1.48 1.22
29 5.67 100 18.9 0.75 2.00 1.28 83 1.89 15 45.2 0.28 1.18 1.66 136 1.21 40 15.8 0.08 1.60 1.20
30 5.67 120 16.4 0.75 2.08 1.21 84 1.89 20 41.7 0.28 1.30 1.62 137 1.21 60 9.5 0.08 1.78 0.98
31 4.25 5 120 0.63 0.70 2.08 85 1.89 30 31.4 0.28 1.48 1.50 138 1.21 80 8.1 0.08 1.90 0.91
32 4.25 10 81.5 0.63 1.00 1.91 86 1.89 40 25.1 0.28 1.60 1.40 139 1.21 100 10.9 0.08 2.00 1.04
33 4.25 15 54.8 0.63 1.18 1.74 87 1.89 60 17.8 0.28 1.78 1.25 140 1.21 120 9.5 0.08 2.08 0.98
34 4.25 20 48 0.63 1.30 1.68 88 1.89 80 13.7 0.28 1.90 1.14 141 1.13 5 48 0.05 0.70 1.68
35 4.25 30 41.8 0.63 1.48 1.62 89 1.89 100 11.7 0.28 2.00 1.07 142 1.13 10 27 0.05 1.00 1.43
36 4.25 40 34.8 0.63 1.60 1.54 90 1.89 120 10.3 0.28 2.08 1.01 143 1.13 15 25.6 0.05 1.18 1.41
37 4.25 60 29.2 0.63 1.78 1.47 91 1.70 5 84 0.23 0.70 1.92 144 1.13 20 0.05 1.30 #¡NUM!
38 4.25 80 22.9 0.63 1.90 1.36 92 1.70 10 51 0.23 1.00 1.71 145 1.13 30 16 0.05 1.48 1.20
39 4.25 100 18.7 0.63 2.00 1.27 93 1.70 15 36 0.23 1.18 1.56 146 1.13 40 12.8 0.05 1.60 1.11
40 4.25 120 16.3 0.63 2.08 1.21 94 1.70 20 28.8 0.23 1.30 1.46 147 1.13 60 8 0.05 1.78 0.90
41 3.40 5 120 0.53 0.70 2.08 95 1.70 30 23.4 0.23 1.48 1.37 148 1.13 80 6 0.05 1.90 0.78
42 3.40 10 69 0.53 1.00 1.84 96 1.70 40 20.7 0.23 1.60 1.32 149 1.13 100 7.1 0.05 2.00 0.85
43 3.40 15 52 0.53 1.18 1.72 97 1.70 60 16.5 0.23 1.78 1.22 150 1.13 120 6.3 0.05 2.08 0.80
44 3.40 20 47.4 0.53 1.30 1.68 98 1.70 80 13.4 0.23 1.90 1.13 151 1.06 5 36 0.03 0.70 1.56
45 3.40 30 41 0.53 1.48 1.61 99 1.70 100 11.6 0.23 2.00 1.06 152 1.06 10 0.03 1.00 #¡NUM!
46 3.40 40 34.2 0.53 1.60 1.53 100 1.70 120 10.1 0.23 2.08 1.00 153 1.06 15 0.03 1.18 #¡NUM!
47 3.40 60 25.7 0.53 1.78 1.41 101 1.55 5 78 0.19 0.70 1.89 154 1.06 20 0.03 1.30 #¡NUM!
48 3.40 80 21.5 0.53 1.90 1.33 102 1.55 10 48 0.19 1.00 1.68 155 1.06 30 0.03 1.48 #¡NUM!
49 3.40 100 18.2 0.53 2.00 1.26 103 1.55 15 34 0.19 1.18 1.53 156 1.06 40 12 0.03 1.60 1.08
50 3.40 120 15.9 0.53 2.08 1.20 104 1.55 20 28.5 0.19 1.30 1.45 157 1.06 60 0.03 1.78 #¡NUM!
51 2.83 5 120 0.45 0.70 2.08 105 1.55 30 23.4 0.19 1.48 1.37 158 1.06 80 6 0.03 1.90 0.78
52 2.83 10 60 0.45 1.00 1.78 106 1.55 40 18.8 0.19 1.60 1.27 159 1.06 100 5.3 0.03 2.00 0.72
53 2.83 15 52 0.45 1.18 1.72 107 1.55 60 14 0.19 1.78 1.15 160 1.06 120 4.7 0.03 2.08 0.67
54 2.83 20 44.7 0.45 1.30 1.65

-6-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

Los parámetros de regresión se pueden obtener aplicando el método de mínimos cuadrados, que
consiste en minimizar el error al cuadrado, entre los datos y el modelo lineal.

3. Análisis de Regresión Múltiple

Este tipo se presenta cuando dos o más variables independientes influyen sobre una variable
dependiente. Ejemplo: Y = f(x, w, z).
Dispone de una ecuación con dos variables independientes adicionales:

Se puede ampliar para cualquier número “m” de variables independientes:

Para poder resolver y obtener a, b1 y b2 en una ecuación de regresión múltiple el cálculo se

presenta muy tedioso porque se tiene atender 3 ecuaciones que se generan por el método de
mínimo de cuadrados:

Para poder resolver se puede utilizar programas informáticos como AD+, SPSS y Minitab y Excel.

4. El coeficiente de determinación múltiple (R2)

Mide la tasa porcentual de los cambios de Y que pueden ser explicados por x1, x2 y x3
simultáneamente.

SCregresión = suma de cuadrados de la regresión

SCtotal = suma de cuadrados totales

∑ (̂ ̅)
∑ ( ̅)

En el siguiente cuadro se muestran los resultados de análisis de regresión múltiple.

-7-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

Cuadro 7. Resultados de regresión múltiple con el software Minitab 15

The regression equation is
log(Imáx) = 2.19 + 0.557 log(T) - 0.643 log(D)
154 cases used, 6 cases contain missing values

Predictor Coef SE Coef T P

Constant 2.19035 0.02227 98.36 0.000
log(T) 0.55730 0.01766 31.56 0.000
log(D) -0.64322 0.01333 -48.25 0.000

S = 0.0723159 R-Sq = 95.7% R-Sq(adj) = 95.7%

Analysis of Variance
Source DF SS MS F P
Regression 2 17.6561 8.8281 1688.10 0.000
Residual Error 151 0.7897 0.0052
Total 153 18.4458

El valor de es K = 102.19 = 154.88

El coeficiente de determinación hallado es R2 = 0.957 muestra que el 95.7% de la variación de la

intensidad es explicado por el período de retorno y la duración. Se considera como adecuado el
modelo.

Los valores a y b se determinaron directamente, por tanto, la fórmula de Intensidad máxima es

La intensidad máxima para una duración de 45 min y 10 años de período de retorno es

Con la fórmula de Intensidad máxima obtenida, se puede obtener curvas Intensidad-Duración-

Período de retorno.

Los siguientes datos de Intensidad máxima en mm/hr se obtuvieron aplicando la ecuación hallada

Cuadro 7. Intensidad máxima obtenida con la formula

Duración (min.)
T
(años) 5 10 15 20 30 40 60 80 100 120
0.25 25.42 16.28 12.54 10.43 8.03 6.68 5.14 4.28 3.70 3.29
0.5 37.40 23.95 18.45 15.34 11.82 9.82 7.57 6.29 5.45 4.85
2 80.95 51.84 39.94 33.20 25.58 21.26 16.38 13.61 11.79 10.49
5 134.86 86.36 66.54 55.30 42.61 35.42 27.29 22.68 19.65 17.47
10 198.41 127.06 97.90 81.36 62.69 52.10 40.15 33.37 28.91 25.71
20 291.90 186.93 144.03 119.70 92.23 76.66 59.06 49.09 42.53 37.82

El siguiente gráfico muestra las curvas intensidad – duración – período de retorno.

-8-
Curso de Hidrología. Docente Eduardo Luis Flores Quispe

Curvas Intensidad-Duración-Período de retorno

350.00

300.00

250.00
Intensidad máxima (mm/hr)

T = 0.25 años
200.00
T = 0.5 años
T = 2 años
150.00
T = 5 años
T = 10 años
100.00
T = 20 años

50.00

0.00
0 20 40 60 80 100 120 140
Duración (min.)

-9-