Curvas IDF en Hidrología: Análisis y Resultados
Curvas IDF en Hidrología: Análisis y Resultados
CURVAS IDF
1. Introducción........................................................................................................................................... 3
2. Objetivos ................................................................................................................................................ 4
2.1 Objetivo general .................................................................................................................................. 4
2.2 Objetivos específicos ........................................................................................................................... 4
3. Datos de intensidades máximas............................................................................................................. 5
4. Cálculos y resultados ............................................................................................................................. 6
4.1. Funciones de probabilidad para diferentes duraciones ................................................................. 7
4.2. Curvas IDF .................................................................................................................................... 17
4.3. Relaciones intensidad-duración-frecuencia ................................................................................. 18
4.3.1. Método gráfico .......................................................................................................................... 18
4.3.2. Regresión múltiple ..................................................................................................................... 27
4.4. Comparación de los métodos ...................................................................................................... 29
5. Conclusiones ........................................................................................................................................ 30
6. Referencias .......................................................................................................................................... 31
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
1. Introducción
El estudio de precipitaciones en una región es importante debido a las actividades relacionadas con el
aprovechamiento hídrico. Se debe estudiar las precipitaciones y su distribución temporal para fines
meteorológicos e hidrológicos, con lo cual se pueden proporcionar índices para realizar estudios de crecidas
o alimentar los modelos precipitación-escorrentía que mejora la información disponible para el adecuado
diseño y dimensionamiento de obras civiles. Por lo cual es necesario conocer las intensidades de
precipitación para los distintos periodos de retorno.
En el presente informe se pretende mostrar el procedimiento realizado para la construcción de las curvas
IDF para diferentes periodos de retorno y los resultados obtenidos, haciendo uso de los datos de
intensidades máximas de diferentes bandas pluviográficas encontradas y realizando un tratamiento
estadístico de los datos con el fin de realizar un ajuste de los datos mediante la distribución de Gumbel o
de valores extremos.
3
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
2. Objetivos
2.1 Objetivo general
Graficar las curvas IDF para los datos de intensidades máximas dados y obtener las ecuaciones que
describan la relación intensidad-duración- frecuencia
− Determinar las funciones de probabilidades para las intensidades de acuerdo con su duración.
− Calcular los parámetros necesarios para realizar los ajustes de los datos.
− Determinar los valores de intensidad máxima para diferentes periodos de retorno.
4
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Los datos con los cuales se construyen las curvas IDF, corresponden a las intensidades máximas en mm/h
calculados mediante los valores de precipitaciones máximas de 37 bandas pluviográficas para ciertos
intervalos de tiempo (véase Tabla 1).
Intensidades (mm/h)
15 min 30min 1h 2h 3h 6h 12h 24h
32 30 30 27.5 27 16.3 11.67 10.21
90 90 90 45 43.3 23.33 18.3 14.58
63 45 32 29 27.8 26.2 22 18
20.8 14.4 10 5.75 4 2 1 0.5
64 61 60 42 31.3 20.67 11.75 5.92
32 18 16 10.8 11.5 7.8 5.3 2.9
90 90 90 66.5 54.67 34.58 18.42 13.37
54.7 45.3 32.2 20.5 18.1 17.5 12.6 6.3
61.6 60 59 39.5 33 18.7 11.8 5.9
50 50 50 35 34.6 23.3 20.8 14.6
100 100 100 100 86.66 51.66 28.33 23.75
70 70 45 2.87 2.5 2.08 1.71 0.91
50 50 50 35 30 28.33 20 17.5
80 64 33 23.7 20.8 17.5 15.4 8
90.8 77.46 44.7 14.8 9.56 6.36 3.18 1.54
20 20 20 12.5 9.57 5.18 3.17 2.72
30 26.1 20 19.1 15.5 9.8 3.9 2.3
100 100 100 97 65 36.83 18.5 9.38
78.4 65.8 39 25.65 24.93 14.13 8.52 4.37
132 71 36 20 13.3 6.83 3.5 1.75
100 100 100 100 90 57.5 29.17 20.41
92 90 90 65 43.3 26.67 17.92 16.67
120 100 70 48.5 36.67 26.67 14.17 7.08
45.6 36.8 26.25 19.24 13.2 7.73 6.04 5.52
81.2 70 63 41.5 34.3 17.17 11.25 6.75
5 5 5 5 4.33 2.5 1.75 0.87
20 14 8 4.75 3.5 1.95 1.42 0.8
50 50 50 45 40 22.83 11.5 5.57
40 28 27.5 26 23.3 15.8 11.4 10.7
12 9 7 4.3 3.7 2.9 2.4 2
85 85 85 52.5 42.7 24.2 23.7 15.5
32 20 20 9.5 7.2 2.08 0.97 0.26
340 200 132 76 60.66 38.66 32.91 23
60 60 60 42.5 33.3 29.5 14.75 7.38
15.2 10.8 9 6.2 6 3.7 2.07 2.04
160 160 100 100 93.3 90 60 39.16
18 9 4.5 2.405 1.62 0.82 0.413 0.416
Tabla 1 Datos iniciales de intensidades máximas en mm/h para diferentes intervalos de tiempo
5
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
4. Cálculos y resultados
Con los datos de la Tabla 1, se determina una función de probabilidad aproximada a la distribución de
Gumbel, la cual se utiliza para el estudio de los valores extremos, por ejemplo, en este caso corresponde a
las intensidades máximas. La probabilidad que presente un valor inferior a 𝑥𝑖 según la distribución de
Gumbel es:
−𝛼(𝑥𝑖 −µ)
𝐹(𝑥 ≤ 𝑥𝑖 ) = 𝑒 −𝑒 Ecuación (1)
Donde α y µ corresponden a los parámetros del ajuste, y se calculan mediante las siguientes ecuaciones:
𝑆𝑛 Ecuación (2)
µ=
𝑆
𝑌𝑛 Ecuación (3)
𝛼 = 𝑥̅ − 𝑆
𝑆𝑛
𝑆𝑛 y 𝑌𝑛 dependen del número de datos utilizados en cada ajuste y son respectivamente la media y la
desviación estándar de la variable reducida. En este caso, los valores de Sn y Yn se determinan mediante el
cuadro 3.6. del libro “Hidrología en la Ingeniería”, donde se revisan los valores para un tamaño de muestra
de 37 datos (véase Tabla 2).
Con las ecuaciones señaladas anteriormente se calculan los parámetros del ajuste de acuerdo con las
duraciones de las intensidades.
6
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Los datos utilizados para determinar los parámetros corresponden a la columna 1 (15 min) de la Tabla 1, a
partir de estos se calculan los valores mencionados anteriormente para el cálculo de los parámetros de la
función de probabilidad:
Como se mencionó anteriormente, los valores de Yn y Sn se determinan mediante la tabla para un tamaño
de muestra de 37 (véase Tabla 2).
Los datos necesarios para el cálculo de los parámetros se muestran en la siguiente tabla resumen:
n 37
x̄ 69.87
µ2 8182.63
s2 3300.40
s 57.45
Yn 0.5418
Sn 1.1339
Tabla 3 Datos encontrados para duración de 15 min para el cálculo de los parámetros del ajuste.
A continuación, se calculan los parámetros del ajuste haciendo uso de la Ecuación (2) y la Ecuación (3) de la
siguiente manera:
0.5418
𝛼 = 69.87 − 57.45 = 42.42
1.1339
1.1339
µ= = 0.020
57.45
La siguiente tabla, muestra los parámetros obtenidos para el ajuste de la duración de 15 minutos:
7
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Una vez determinados los parámetros, se reemplazan los valores en la Ecuación (1) y se obtiene la función
de probabilidad para el primer intervalo de tiempo:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Al hacer uso de la Ecuación (7) se determinan las intensidades máximas correspondientes a ciertos periodos
de retorno despejando de la siguiente ecuación la probabilidad de que ocurra una intensidad menor o igual
a la que se quiere determinar:
1 1 Ecuación (8)
𝑇𝑟 = =
𝑝(𝑥 ≥ 𝑥𝑖 ) 1 − 𝑝(𝑥 ≤ 𝑥𝑖 )
Donde Tr hace referencia al periodo de retorno. Por ejemplo, para un periodo de retorno de 5 años, el valor
correspondiente a la probabilidad de que ocurra un evento menor o igual a una intensidad máxima Xi es:
1
5=
1 − 𝑝(𝑥 ≤ 𝑥𝑖 )
1
𝑝(𝑥 ≤ 𝑥𝑖 ) = 1 − = 0.8
5
Y la intensidad máxima correspondiente a este periodo de retorno será:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 118.42
10 0.1 0.9 156.44
25 0.04 0.96 204.48
50 0.02 0.98 240.11
100 0.01 0.99 275.49
Tabla 5 Intensidades máximas de diferentes periodos de retorno para una duración de 15 minutos.
De la misma manera se aplica el método para las siguientes duraciones. A continuación, se presentan las
tablas que contienen los resultados obtenidos para cada uno de los pasos mostrados anteriormente.
8
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 2 (30 min) de la Tabla 1 y haciendo uso de las ecuaciones mostradas
anteriormente se calculan los parámetros (véase Tabla 6).
n 37
x̄ 59.34
µ2 5255.14
s2 1733.65
s 41.64
Yn 0.5418
Sn 1.1339
Tabla 6 Datos encontrados para duración de 30 min para el cálculo de los parámetros del ajuste
Función de densidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los
periodos de retorno mostrados en la tabla:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 94.53
10 0.1 0.9 122.08
25 0.04 0.96 156.90
50 0.02 0.98 182.73
100 0.01 0.99 208.37
Tabla 8 Intensidades máximas de diferentes periodos de retorno para una duración de 30 minutos.
9
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 3 (1h) de la Tabla 1 y haciendo uso de las ecuaciones mostradas
anteriormente se llega a los siguientes valores:
n 37
x̄ 49.03
µ2 3517.01
s2 1112.97
s 33.36
Yn 0.5418
Sn 1.1339
Tabla 9 Datos encontrados para duración de 1 hora para el cálculo de los parámetros del ajuste
Mediante los anteriores valores se llega a los siguientes parámetros del ajuste:
Función de densidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los periodos
de retorno:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 77.22
10 0.1 0.9 99.30
25 0.04 0.96 127.20
50 0.02 0.98 147.89
100 0.01 0.99 168.43
Tabla 11 Intensidades máximas de diferentes periodos de retorno para una duración de 1 hora.
10
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 4 (2 h) de la Tabla 1 y haciendo uso de las ecuaciones mostradas
anteriormente se llega a los siguientes valores necesarios para el cálculo de los parámetros del ajuste:
n 37
x̄ 35.69
µ2 2110.14
s2 836.30
s 28.92
Yn 0.5418
Sn 1.1339
Tabla 12 Datos encontrados para duración de 2 horas para el cálculo de los parámetros del ajuste.
Mediante los anteriores valores se llega a los siguientes parámetros del ajuste:
Tabla 13 Parámetros de la función de probabilidad para duración de 2 horas.
Función de probabilidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los
periodos de retorno:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 60.13
10 0.1 0.9 79.27
25 0.04 0.96 103.45
50 0.02 0.98 121.39
100 0.01 0.99 139.19
Tabla 15 Intensidades máximas de diferentes periodos de retorno para una duración de 2 horas.
11
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 5 (3 h) de la Tabla 1 y haciendo uso de las ecuaciones mostradas
anteriormente se llega a los siguientes valores necesarios para el cálculo de los parámetros del ajuste:
n 37
x̄ 29.73
µ2 1477.51
s2 593.38
s 24.36
Yn 0.5418
Sn 1.1339
Tabla 16 Datos encontrados para duración de 3 horas para el cálculo de los parámetros del ajuste.
Mediante los anteriores valores se llega a los siguientes parámetros del ajuste:
Función de probabilidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los
periodos de retorno:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 50.32
10 0.1 0.9 66.44
25 0.04 0.96 86.81
50 0.02 0.98 101.92
100 0.01 0.99 116.92
Tabla 18 Intensidades máximas de diferentes periodos de retorno para una duración de 3 horas.
12
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 6 (6 h) de la Tabla 1 y haciendo uso de las ecuaciones mostradas
anteriormente se llega a los siguientes valores necesarios para el cálculo de los parámetros del ajuste:
n 37
x̄ 19.99
µ2 726.30
s2 326.56
s 18.07
Yn 0.5418
Sn 1.1339
Tabla 19 Datos encontrados para duración de 6 horas para el cálculo de los parámetros del ajuste.
Mediante los anteriores valores se llega a los siguientes parámetros del ajuste:
Función de probabilidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los
periodos de retorno:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 35.26
10 0.1 0.9 47.22
25 0.04 0.96 62.33
50 0.02 0.98 73.54
100 0.01 0.99 84.67
Tabla 21 Intensidades máximas de diferentes periodos de retorno para una duración de 6 horas.
13
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 7 (12 h) de la Tabla 1 y haciendo uso de las ecuaciones se llega a los
siguientes valores necesarios para el cálculo de los parámetros del ajuste:
n 37
x̄ 13.02
µ2 307.40
s2 137.92
s 11.74
Yn 0.5418
Sn 1.1339
Tabla 22 Datos encontrados para duración de 12 horas para el cálculo de los parámetros del ajuste.
Mediante los anteriores valores se llega a los siguientes parámetros del ajuste:
Función de probabilidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los
periodos de retorno:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 22.94
10 0.1 0.9 30.71
25 0.04 0.96 40.53
50 0.02 0.98 47.82
100 0.01 0.99 55.05
Tabla 24 Intensidades máximas de diferentes periodos de retorno para una duración de 12 horas.
14
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Con los datos mostrados en la columna 8 (24 h) de la Tabla 1 y haciendo uso de las ecuaciones se llega a los
siguientes valores necesarios para el cálculo de los parámetros del ajuste:
n 37
x̄ 8.88
µ2 149.85
s2 70.96
s 8.42
Yn 0.5418
Sn 1.1339
Tabla 25 Datos encontrados para duración de 24 horas para el cálculo de los parámetros del ajuste
Mediante los anteriores valores se llega a los siguientes parámetros del ajuste:
Función de probabilidad:
Función de probabilidad
− . (𝑥𝑖− . )
𝐹 𝑥≤𝑥 = 𝑒 −𝑒
Intensidades máximas:
Despejando de la función de la Ecuación (9) se llega a las siguientes intensidades máximas para los
periodos de retorno:
Tr F(x ≥ x i ) F(x ≤ x i ) xi
5 0.2 0.8 16.00
10 0.1 0.9 21.58
25 0.04 0.96 28.62
50 0.02 0.98 33.84
100 0.01 0.99 39.03
Tabla 27 Intensidades máximas de diferentes periodos de retorno para una duración de 24 horas.
15
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Una vez encontrados todos los valores de intensidades máximas para las duraciones y periodos de retorno
correspondientes, se procede a graficar las curvas IDF, para lo cual se realizó las siguientes tablas en las que
se recopila toda la información necesaria.
Tabla 28 Datos de intensidades para ciertos periodos de retorno y duraciones.
Periodo de
Duración, t Intensidad, i
retorno, Tr
(min) (mm/h)
(años)
15 118.42
30 94.53
60 77.22
120 60.13
5
180 50.32
360 35.26
720 22.94
1440 16.00
15 156.44
30 122.08
60 99.30
120 79.27
10
180 66.44
360 47.22
720 30.71
1440 21.58
15 204.48
30 156.90
60 127.20
120 103.45
25
180 86.81
360 62.33
720 40.53
1440 28.62
15 240.11
30 182.73
60 147.89
120 121.39
50
180 101.92
360 73.54
720 47.82
1440 33.84
15 275.49
30 208.37
60 168.43
120 139.19
100
180 116.92
360 84.67
720 55.05
1440 39.03
16
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Se grafica intensidades vs duraciones para cada periodo de retorno, es decir que a cada punto de una misma
curva en la gráfica (véase Gráfica 1) le corresponde un mismo periodo de retorno. Las curvas IDF se
muestran a continuación:
Gráfica 1 Curvas IDF para tiempo en horas.
Curvas IDF
300
250
200
Intensidad, i (mm/h)
Tr = 5 años
150 Tr = 10 años
Tr = 25 años
Tr = 50 años
Tr = 100 años
100
50
0
0 5 10 15 20 25
Duración, t (hora)
Curvas IDF
300
250
200
Intensidad, i (mm/h)
Tr = 5 años
150 Tr = 10 años
Tr = 25 años
Tr = 50 años
Tr = 100 años
100
50
0
0 200 400 600 800 1000 1200 1400
Duración, t (minutos)
17
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Ahora se determina la ecuación que describa la relación entre la intensidad la duración y la frecuencia. Para
esto, se hizo uso de dos métodos, el método gráfico y la regresión múltiple; en ambos se busca encontrar
los parámetros que permitan caracterizar la relación entre las variables mencionadas mediante una
ecuación de tipo:
𝐶 𝑘𝑇𝑟 𝑚 Ecuación (15)
= =
(𝑡 + 𝑡0 )𝑛 (𝑡 + 𝑡0 )𝑛
Donde los parámetros son m, n, k y t0; y las variables son i (intensidad en mm/h), Tr (periodo de retorno en
años) y t (duración en minutos).
Para la determinación de los parámetros por el método gráfico se coloca en coordenadas logarítmicas las
intensidades máximas en función de la duración para cada periodo de retorno, de tal manera que la
ecuación que representa la relación es la siguiente:
El valor de to debe ser un valor escogido por medio del método de ensayo y error, de tal manera que para
el valor de to los puntos graficados se dispongan mejor a lo largo de la línea de tendencia graficada en Excel.
Para lo anterior se presenta en el gráfico el valor de R2, es decir, la correlación de los datos con la regresión
lineal. A medida que la correlación se aproxime más a 1, mayor será la relación de los datos con la
linealización. Debido a que la gráfica no es exactamente una recta se debe agregar un valor t0 para ajustar
los datos de tal manera que se alcance el máximo valor de correlación.
El valor de to se puede determinar variando dicho valor de forma manual o implementando la herramienta
de solver en Excel. Para este caso, se hace uso de la herramienta solver. En caso de que la herramienta de
solver no se encuentra activada, se debe habilitar en el menú de complementos del software (véase Figura
1). Esta herramienta es muy útil, ya que implementa métodos numéricos con los que es posible determinar
el valor de una variable de acuerdo a parámetros o condiciones establecidas.
18
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Figura 1 Opción solver para determinar el máximo valor de to para cada tiempo de retorno.
En el método gráfico se debe determinar el valor de to para el cual el valor de la correlación sea el máximo
posible. Por lo tanto, se implementa la opción solver de Excel para conocer el valor to correspondiente a
cada uno de los periodos de retorno.
Para utilizar la herramienta de solver, es necesario realizar la regresión lineal con los datos disponibles para
después determinar el tiempo de ajuste to para el cual se consigue una mayor cercanía al 1 en la correlación.
Con solver se realiza la iteración de to y de esta manera se conoce el valor exacto del tiempo que ajusta
mejor los datos a la regresión realizada.
Este procedimiento se realiza para cada uno de los tiempos de retorno que se tienen.
19
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Tr 5 años
2.5
2.0
1.5
y = -0.5872x + 3.0633
log i
R² = 0.9984
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
log (t+t0 )
Se compara la ecuación de la anterior gráfica con la ecuación 16 y se observa que el corte con el eje y (b)
corresponde a 𝐿𝑜𝑔10 (𝐶) y la pendiente corresponde a n.
20
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Tr 10 años
2.5
2.0
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
log (t+t0 )
Para un periodo de retorno de 10 años se obtuvieron los siguientes resultados del análisis.
21
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Tr 25 años
2.5
2.0
y = -0.5547x + 3.2165
1.5
R² = 0.996
log i
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
log (t+t0 )
Para un periodo de retorno de 25 años se obtuvieron los siguientes resultados del análisis.
22
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Graficando las coordenadas anteriores se obtiene la siguiente gráfica con su respectiva ecuación:
Gráfica [Link] de coordenadas logarítmicas de intensidad en función de la duración para Tr= 50 años.
Tr 50 años
2.5
2.0
y = -0.5479x + 3.2686
1.5 R² = 0.9953
log i
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
log (t+t0 )
Donde los resultados obtenidos a partir del análisis de la ecuación son para un periodo de retorno de 50
años:
23
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
15 275.49 100 2.44 1.66 30.758 4.05 2.76 2.41 0.0007 0.1467
30 208.37 100 2.32 1.78 30.758 4.14 3.18 2.35 0.0008 0.0685
60 168.43 100 2.23 1.96 30.758 4.36 3.83 2.25 0.0007 0.0287
120 139.19 100 2.14 2.18 30.758 4.67 4.74 2.13 0.0001 0.0075
180 116.92 100 2.07 2.32 30.758 4.81 5.40 2.05 0.0002 0.0001
360 84.67 100 1.93 2.59 30.758 5.00 6.72 1.91 0.0004 0.0167
720 55.05 100 1.74 2.88 30.758 5.01 8.27 1.75 0.0002 0.1001
1440 39.03 100 1.59 3.17 30.758 5.04 10.03 1.60 0.0000 0.2169
Total 16.46 18.54 37.06 44.94 0.0031 0.5851
Graficando las anteriores coordenadas se obtiene la siguiente gráfica con su respectiva ecuación:
Gráfica 7. Grafica de coordenadas logarítmicas de intensidad en función de la duración para Tr= 100 años.
Tr 100 años
3.0
2.5
2.0
y = -0.5429x + 3.3152
log i
1.5 R² = 0.9948
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
log (t+t0 )
A partir de la gráfica anterior se obtiene que para un periodo de retorno de 100 años los resultados de
análisis son:
24
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Resultados
Tr (años) ni n t0 t (min) log C
5 -0.5872 36.447 3.0633
10 -0.5686 34.137 3.1361
25 -0.5547 0.5603 32.342 33.025 3.2165
50 -0.5479 31.439 3.2686
100 -0.5429 30.758 3.3152
Ahora, como se quiere determinar el valor de los parámetros k y m (véase la ecuación 16) se realiza el mismo
análisis para la ecuación 16. Se indica en coordenadas logarítmicas los valores de c en función de los
periodos de retorno. A partir de lo anterior se obtiene:
De los resultados obtenidos anteriormente se realiza la siguiente tabla de los valores de Log C para Log Tr
asociado:
log C log Tr
3.0633 0.6990
3.1361 1.0000
3.2165 1.3979
3.2686 1.6990
3.3152 2.0000
3.3
y = 0.1932x + 2.9374
3.3 R² = 0.9932
log C
3.2
3.2
3.1
3.1
0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
log Tr
25
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
C= 𝑘𝑇𝑟 𝑚
log C= − 𝑜𝑔𝑇𝑟
A partir de la anterior gráfica y de la Ecuación (18) se encuentra que los valores de Log K y de m son el
corte con el eje y y la pendiente, respectivamente (véase Tabla 37).
Resultados
m log k k
0.1932 2.9374 865.695
𝑘 = 102.9374 = 865.695
Finalmente se determinan los parámetros restantes a partir de los resultados anteriores de la siguiente
manera:
∑ 𝐿𝑜𝑔(𝑐 )
𝐿𝑜𝑔(𝐶) =
𝑛
Donde n es el número de periodos de retorno evaluados y Log (Ci) hace referencia a los resultados
mostrados anteriormente (véase la Tabla 35. Valores de Log (c) y Log (Tr). De lo anterior:
3.0633 + 3.1361 + 3.2165 + 3.2686 + 3.3152
𝐿𝑜𝑔(𝐶) = = 3.1999
5
Y a partir de lo anterior se calcula el valor de c de la siguiente manera:
𝑐 = 103.1999 = 1584.654
Una vez hallado el valor de c, se calcula el parámetro n, el cual se calcula como el promedio de los valores
obtenidos a partir del anterior análisis (véase la Tabla 34).
−0.5872 − 0.5686 − 0.5547 − 0.5479 − 0.5429
𝑛= = 0.5603
5
26
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Se determina ahora el valor de t0 mediante el promedio de los datos mostrados en la Tabla 34:
865.70 𝑇𝑟 0.19
=
(𝑡 + 33.02)0.56
4.3.2. Regresión múltiple
El método de la regresión múltiple se realizó mediante una hoja de cálculo en Excel haciendo uso de la
función estimació[Link] la cual determina los valores de los parámetros directamente, sin embargo en este
método no se determina el valor de t0 sino se calculan los parámetros para una ecuación de tipo:
Se elabora una tabla de Excel en la que se contiene la información necesaria para realizar la regresión, es
decir una tabla con los valores de la variable explicada y los valores de las variables explicativas, como se
muestra a continuación:
Periodo de
Intensidad, i Duración, t
retorno, Tr Log (i) Log (Tr) Log (t)
(mm/h) (min)
(años)
118.42 15 5 2.07 0.70 1.18
156.44 15 10 2.19 1.00 1.18
204.48 15 25 2.31 1.40 1.18
240.11 15 50 2.38 1.70 1.18
275.49 15 100 2.44 2.00 1.18
94.53 30 5 1.98 0.70 1.48
122.08 30 10 2.09 1.00 1.48
156.90 30 25 2.20 1.40 1.48
182.73 30 50 2.26 1.70 1.48
208.37 30 100 2.32 2.00 1.48
77.22 60 5 1.89 0.70 1.78
99.30 60 10 2.00 1.00 1.78
127.20 60 25 2.10 1.40 1.78
147.89 60 50 2.17 1.70 1.78
168.43 60 100 2.23 2.00 1.78
60.13 120 5 1.78 0.70 2.08
79.27 120 10 1.90 1.00 2.08
103.45 120 25 2.01 1.40 2.08
121.39 120 50 2.08 1.70 2.08
139.19 120 100 2.14 2.00 2.08
50.32 180 5 1.70 0.70 2.26
27
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Finalmente se escribe la función en una celda de la hoja de cálculo y se llega a los siguientes parámetros:
Log k m n
2,45 0,28 0,43
𝑘 = 102.45 = 280.29
Se reemplazan los valores en la Ecuación (19) y se obtiene la siguiente relación:
280.29 𝑇𝑟 0.28
=
(𝑡)0.43
28
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Una vez encontradas las dos relaciones mediante los dos métodos se realiza la siguiente tabla en la que,
para los valores de duración y periodo de retorno, se encuentran los respectivos valores de intensidad,
con el fin de que se pueda comparar el valor obtenido mediante las ecuaciones con los valores reales.
Periodo de Intensidad
Duración, t Ajuste método Ajuste regresión
retorno, Tr Real
(min) gráfico (mm/h) múltiple (mm/h)
(años) (mm/h)
15 5 118.42 134.98 137.44
15 10 156.44 154.33 166.72
15 25 204.48 184.21 215.21
15 50 240.11 210.61 261.06
15 100 275.49 240.79 316.68
30 5 94.53 115.91 102.10
30 10 122.08 132.53 123.86
30 25 156.90 158.19 159.88
30 50 182.73 180.86 193.95
30 100 208.37 206.77 235.26
60 5 77.22 93.20 75.85
60 10 99.30 106.55 92.01
60 25 127.20 127.19 118.78
60 50 147.89 145.41 144.08
60 100 168.43 166.25 174.78
120 5 60.13 70.52 56.35
120 10 79.27 80.62 68.36
120 25 103.45 96.23 88.24
120 50 121.39 110.02 107.04
120 100 139.19 125.79 129.84
180 5 50.32 58.58 47.36
180 10 66.44 66.98 57.45
180 25 86.81 79.95 74.16
180 50 101.92 91.41 89.96
180 100 116.92 104.51 109.12
360 5 35.26 41.57 35.18
360 10 47.22 47.52 42.68
360 25 62.33 56.73 55.09
360 50 73.54 64.86 66.83
360 100 84.67 74.15 81.07
720 5 22.94 28.88 26.14
720 10 30.71 33.01 31.71
720 25 40.53 39.41 40.93
720 50 47.82 45.05 49.65
720 100 55.05 51.51 60.23
1440 5 16.00 19.83 19.42
1440 10 21.58 22.67 23.56
1440 25 28.62 27.06 30.41
1440 50 33.84 30.94 36.88
1440 100 39.03 35.37 44.74
29
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
5. Conclusiones
Al iterar de forma manual el valor de T0 en cada una de las regresiones no se consigue con exactitud el valor
para el cual se da la mayor correlación posible. Por esta razón resulta indispensable el uso de solver en
Excel. Además, el método manual resulta impreciso por el hecho de que en ocasiones las variaciones en la
correlación son inapreciables a simple vista.
El ajuste realizado mediante la regresión múltiple es más confiable que el ajuste por el método gráfico.
865.70 𝑇𝑟 0.19
=
(𝑡 + 33.02)0.56
El gráfico que representa la anterior relación obtenida es:
Gráfica 9. Curvas IDF representadas por la relación hallada en el método gráfico.
280.29 𝑇𝑟 0.28
=
(𝑡)0.43
Las curvas que representan la anterior relación es:
30
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
HIDROLOGÍA
CURVAS IDF
PRIMER SEMESTRE 2020
Gráfica 10. Curvas IDF representadas por la relación hallada mediante la regresión múltiple.
6. Referencias
Monsalve Sáenz, G. (1995). Precipitación. En G. Monsalve Sáenz, Hidrología en la Ingeniería (págs. 97-
113). Bogotá D.C., Colombia: Grupo Editorial 87.
31