AIHR AIIH

XXIII CONGRESO LATINOAMERICANO DE HIDRÁULICA

CARTAGENA DE INDIAS COLOMBIA, SEPTIEMBRE 2008

ESTIMACIÓN DE CURVAS INTENSIDAD-DURACIÓN-FRECUENCIA

A PARTIR DE INFORMACIÓN PLUVIOMÉTRICA

Mario Díaz-Granados y Xiomara Puente*

Centro de Investigación en Ingeniería Ambiental (CIIA), Universidad de los Andes, Bogotá, Colombia; (*) Gradex
Ingeniería S. A., Bucaramanga, Colombia, [email protected]

RESUMEN:

Las curvas intensidad – duración – frecuencia (IDF) son la base para el diseño de estructuras
hidráulicas menores basadas en modelación lluvia – escorrentía. Estas curvas se determinan a partir
de información pluviográfica. Dada la menor cobertura espacial de estaciones pluviográficas en
relación con la densidad de estaciones pluviométricas en Colombia, cobran valor métodos de
estimación de curvas IDF con base en información pluviométrica, i.e. datos diarios de precipitación.
En este trabajo se propone un procedimiento basado en (1) el ajuste de la distribución Gumbel con
estimación de parámetros mediante Momentos Ponderados con Probabilidad a la serie anual de
lluvia máxima diaria, que permite identificar la pendiente de la curva de frecuencia, la cual se
preserva para duraciones menores y además en (2) el uso de ecuaciones de regresión para calcular la
magnitud de las intensidades correspondientes. Los errores de estimación son bajos y la robustez
del método propuesto se verifica mediante la comparación de curvas estimadas con curvas reales en
estaciones pluviográficas no utilizadas en la derivación del método.

ABSTRACT:

Intensity – Duration – Frequency (IDF) curves are the basis for designing minor hydraulic
structures through rainfall – runoff modeling. These curves are determined from recording rainfall
gauges. Due to the lower density of this type of gauges in Colombia in relation to non recording
ones, it is valuable having procedures to estimate IDF curves from daily rainfall data. In this paper a
procedure is proposed based on (1) the Gumbel distribution for fitting maximum daily rainfall
annual series and Probability Weighted Moments for parameter estimation, capturing from it the
slope of the frequency curve which is preserved for lower durations, and (2) the use of regression
equations for calculate their corresponding intensity magnitudes. Estimation errors are low, and the
procedure is verified comparing the estimated curves with the real ones in gauging stations not used
in the derivation of the proposed method.

PALABRAS CLAVES:

Curvas IDF, Análisis de Frecuencia de Precipitación, Intensidad-Duración-Frecuencia.

INTRODUCCIÓN

En el diseño de obras hidráulicas de drenaje se requiere determinar el comportamiento hidrológico

de la cuenca tributaria para la estimación del caudal de diseño. Para tal fin, frecuentemente se
recurre a modelos lluvia – escorrentía, donde un evento extremo de precipitación con un período de
retorno dado se transforma en el correspondiente hidrograma de escorrentía al cual se le asigna el
mismo período de retorno. El evento de precipitación está caracterizado normalmente mediante un
hietograma que describe su distribución temporal, o en modelos simples como el método Racional
mediante la intensidad promedio, cuya determinación en ambos casos parte de la disponibilidad de
las curvas Intensidad – Duración – Frecuencia (IDF). Estas curvas sintetizan el comportamiento
pluvial de la zona, dan la intensidad máxima promedio durante una determinada duración para un
período de retorno dado, y son la base para la modelación lluvia – escorrentía mencionada. En
general, las curvas IDF se determinan a partir de registros pluviográficos identificando las series
anuales de profundidad de precipitación máxima para diferentes duraciones y ajustando curvas de
frecuencia mediante distribuciones de probabilidad como Gumbel, logPearson tipo 3, etc. Sin
embargo, usualmente la red pluviográfica tiene menor cobertura que la pluviométrica, y por lo tanto
es conveniente disponer de métodos aproximados para su estimación en sitios sin información
pluviográfica a partir de datos de estaciones pluviométricas. Este es el objetivo principal del
presente trabajo.

ALGUNOS ANTECEDENTES

Las curvas IDF han sido objeto de muchos estudios. Existen varias ecuaciones paramétricas que
buscan describirlas analíticamente. Por ejemplo, está la ecuación de Bernard: IT,t = aTb/tc (Powell,
1932), donde IT,t es la intensidad máxima promedio para un período de retorno T y una duración t, y
a, b y c son parámetros. Powell sugirió valores genéricos de b y c iguales a 0.25 y 0.75,
respectivamente, con base en sus análisis para varios sitios de EEUU. Froehlich (1995) presenta la
siguiente ecuación: IT,t = dTe/(t+f)g, donde d, e, f y g son parámetros.

Por otro lado, en la literatura existen varias propuestas para estimar curvas IDF. Por ejemplo Bell
(1969), propuso la siguiente ecuación para estimar las curvas Profundidad – Duración – Frecuencia
(PDF): PT,t = (alnT+b)/(ctd+e)P2,60, donde la profundidad de precipitación PT,t para un período de
retorno T (entre 2 y 100 años) y duración t (entre 5 y 120 minutos) se estima a partir de la
profundidad de precipitación con un período de retorno de 2 años y una duración de 60 minutos, y
a, b, c, d y e son parámetros. Si esta última profundidad de precipitación no está disponible, Bell
sugirió estimarla a partir del promedio de la lluvia máxima diaria anual M y del promedio anual de
días con tormenta Nt, con resultados menos satisfactorios. Nótese que si PT,t se divide por t se tiene
IT,t. Kothyari y Garde (1992) analizaron 80 curvas IDF en la India encontrando que el parámetro b
de la ecuación de Powell varió entre 0.18 y 0.26 y el parámetro c entre 0.7 y 0.85. Su ecuación
propuesta fue IT,t = CT0.2P2,24 /t0.71, donde P2,24 es la profundidad de precipitación con un período de
retorno de 2 años y una duración de 24 horas y C es una constante. Según Kothyari y Garde, los
errores de estimación con esta fórmula son menores a ±30% en el 95% de los casos. Vargas y Díaz-
Granados (1998) utilizando 165 curvas IDF en Colombia propusieron la ecuación IT,t =
aTbMdNePAf/tc, donde M ya ha sido previamente definido, N es el número promedio anual de días
con lluvia y PA es la precipitación media anual, siendo a, b, c, d, e, f y g parámetros. De acuerdo con
Vargas y Díaz-Granados, los errores de estimación son menores a ±20% en el 46% de los casos y
menores a ±40% en el 86% de los casos. Finalmente Vélez et al. (2002), con base en la teoría de
escalamiento simple e información horaria de 51 estaciones pluviográficas colombianas,
propusieron la ecuación IT,t = aI24exp(bφT)/t0.83, donde I24 es la intensidad de 24 horas, φT es la
inversa de la distribución normal acumulada estándar para una probabilidad de no excedencia igual
a 1/T y a y b son parámetros. Según sus autores, los mayores errores en la validación de esta
fórmula son del orden del 25% y en promedio están entre el 10 y el 15%. Como puede observarse,
la gran mayoría de las fórmulas presentadas pretenden estimar las curvas IDF o PDF a partir de
información de precipitación diaria, buscando de esta manera suplir con información pluviométrica
la no existencia de información pluviográfica. En este trabajo se pretende contribuir a este propósito
con un procedimiento alterno basado en el ajuste de la distribución Gumbel con parámetros
estimados con Momentos Ponderados con Probabilidad (MPP) para dos casos: (1) PDFs e IDFs
multihorarias de 1 a 60 horas de duración t, y (2) PDFs e IDFs para duraciones entre 10 minutos y
2 horas. La información utilizada corresponde a Colombia.

INFORMACIÓN UTILIZADA

Para el desarrollo del procedimiento de estimación propuesto de estimación de curvas IDF

multihorarias se utilizó información de precipitación horaria de 35 estaciones pluviográficas con sus
correspondientes registros diarios concurrentes, cuya distribución según las zonas hidrográficas del
Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM) se muestra en la Tabla 1.
El registro utilizado en cada una de estas estaciones es mayor a 20 años. Adicionalmente, para
propósitos de validación de la metodología propuesta se usaron 3 estaciones adicionales. En el caso
del procedimiento para IDFs entre 10 minutos y 2 horas, se utilizó información de 157 curvas IDF
provenientes de diferentes fuentes, principalmente de estudios de firmas de consultoría y de
entidades prestadoras de servicios públicos. Estas curvas están en una de las 4 siguientes formas: (1)
ecuación tipo Bernard, (2) ecuación tipo Froehlich, (3) datos tabulares y (4) forma gráfica. La
distribución de estas 157 estaciones por zonas hidrográficas del IDEAM se presenta en la Tabla 2.
Para estas estaciones se recopiló adicionalmente información de precipitación diaria
correspondiente a las respectivas series anuales de lluvia diaria máxima, precipitación media anual
y número promedio de días al año con lluvia. Adicionalmente, para propósitos de validación se
consideraron 2 estaciones adicionales.

Tabla 1.- Distribución zonas hidrográficas información multihoraria

Zona hidrográfica Número
estaciones
Alto Magdalena 8
Atrato 3
Catatumbo 2
Cauca 14
Saldaña 1
Magdalena Medio 5
Sogamoso 2

Tabla 2.- Distribución zonas hidrográficas información IDFs de 10 a 120 minutos

Zona hidrográfica Número Número Número
Zona hidrográfica Zona hidrográfica
estaciones estaciones estaciones
Alto Magdalena 47 Cesar 1 Magdalena Medio 9
Atrato 7 Ciénaga Grande 5 Meta 1
Bajo Magdalena 1 Guajira 6 Micay 1
Catatumbo 4 Guaviare 1 Mira 1
Cauca 40 León 1 Nechí 18
Patía 3 Putumayo 1 Saldaña 1
San Juan 2 Sinú 2 Sogamoso 4
Tomo 1

METODO DE GUMBEL CON MOMENTOS PONDERADOS CON PROBABILIDAD

Greenwood et al. (1979) introdujeron el método de MPP como un procedimiento de estimación

robusto y alternativo a momentos (M) y máxima verosimilitud (MV). La estimación de parámetros
con MPP requiere el cálculo de los momentos Ml,j,k = E[XlFj(1-F)k], donde E[.] es el valor esperado,
l, j y k son números reales y X es una variable aleatoria que sigue una distribución específica F =
F(x) = Pr[X<x] y es expresable en forma inversa. Si j y k son enteros no negativos, el momento
ponderado con probabilidad de orden (l, j, k) es proporcional al l-ésimo momento con respecto al
origen de la estadística de orden (j+1) para una muestra de tamaño n = k + j +1. En el caso en que j
= k = 0 y l un entero no negativo, Ml,0,0 es idéntico al momento convencional de orden l con
respecto al origen. El método es atractivo con la distribución Gumbel pues puede ser resuelto
analíticamente. Landwehr et al. (1979) compararon los métodos de M, MV y MPP para la
distribución Gumbel, concluyendo que MPP es favorable y produce estimativos menos sesgados
que con M y MV, y es estadísticamente robusto especialmente cuando los registros históricos son
cortos, altamente asimétricos o con kurtosis demasiado grandes. La distribución Gumbel está dada
por f(x) = exp[-(x – m)/a –exp(x - m)/a]/a donde a y m son los parámetros de la distribución y x es la
magnitud del evento extremo.

Según Greenwood et al. (1979), la expresión con base en la distribución Gumbel para la magnitud
de la variable con un período de retorno T es:

X T = m − a ln[− ln(1 − 1 T )] [1]

donde

m = M 0 − εa [2]

a = (M 0 − 2M 1 ) / ln(2) [3]

y ε = 0.5772 y M k = M l ,0, k . Mk para k entero no negativo se estima como:

1 n ⎛ n − i ⎞ ⎛ n − 1⎞
Mk = ∑ {xi ⎜⎜ ⎟ /⎜
n i =1 ⎝ k ⎟⎠ ⎜⎝ k ⎟⎠
⎟} [4]

siendo n el número de datos históricos y xi el valor del dato i ordenados de menor a mayor.
Operando la expresión para Mk se obtiene:

1 n 1 n n−i
M0 = ∑ xi
n i =1
M1 = ∑ xi
n i =1 n − 1
[5]

con base en las cuales se puede calcular a y m mediante las ecuaciones 2 y 3.

ANÁLISIS DE LA INFORMACIÓN

La metodología propuesta se basa en la estimación de curvas PDF y no directamente de curvas IDF.

De esta manera se garantiza que la profundidad de precipitación para un determinado período de
retorno T es monotónicamente creciente con la duración t. El análisis de la información conducente
al procedimiento de estimación propuesto de curvas PDF consta de los siguientes pasos para el caso
de la información multihoraria:

Determinación de las series anuales de profundidad máxima de precipitación: con base en los
registros horarios de las 35 estaciones pluviográficas se determinaron las series anuales de
profundidad máxima de precipitación para duraciones t de 1, 2, 3, 4, 5, 6, 7, 8, 9 10, 11, 12, 13, … ,
24, 36, 48 y 60 horas y de profundidad máxima diaria del día hidrológico (7 am a 7am).
Análisis preliminar de frecuencia de las series anuales: para las series anuales de profundidad
máxima de precipitación en 6 estaciones escogidas aleatoriamente entre las 35 disponibles, se
ajustaron las siguientes distribuciones con sus correspondientes métodos de estimación de
parámetros: Gumbel con M, MV y MPP; lognormal de tres parámetros con M y MV, y logPearson
tipo 3 con M, Promedios de Sundry (SAM) y momentos aplicados a los logaritmos de los datos
(WRC). Los análisis de frecuencia se realizaron con el programa computacional Hyfran.
Análisis de bondad de ajuste de distribuciones: con base en comparaciones gráficas de las
distribuciones ajustadas y los datos muestrales, en pruebas de hipótesis sobre los datos muestrales
de independencia (Wald-Wolfowitz) y de estacionaridad (Kendall), en pruebas de hipótesis sobre la
aceptabilidad de las distribuciones ajustadas y en criterios de menor error estándar, se determinó
que la distribución de mejor desempeño para las 6 estaciones y todas las duraciones es la
distribución Gumbel con parámetros estimados con MPP. En consecuencia, se adoptó para todas las
35 estaciones la distribución Gumbel con el método de MPP de estimación de parámetros.
Análisis puntual y global de valores de M0 y M1: Greis y Wood (1981) propusieron un método de
estimación regional de frecuencias de caudales máximos basado en el método de MPP para la
estimación de parámetros y lo ilustraron con Gumbel. Díaz-Granados y Deeb (1988) aplicaron con
algunas modificaciones este método de análisis regional en la cuenca del Alto Magdalena con
resultados satisfactorios comparados con otros métodos. En este trabajo de estimación de curvas
PDF se aplican elementos del método de Greis y Wood. Así, el procedimiento utilizado es el
siguiente:
• Calcular M0, M1 y M1/M0 para cada estación y cada duración t. En la Tabla 3 se muestra un
ejemplo de esta última relación para 8 de las 35 estaciones. Obsérvese la notable
preservación de un valor muy similar de M1/M0 en cada estación para diferentes duraciones
t. Para estimar las curvas PDF esto permite adoptar el valor de (M1/M0)diario estimado con la
serie anual de lluvia diaria máxima (obtenible de registros pluviométricos) como
representativo de otras duraciones t. El rango de variación de M1/M0 en las 35 estaciones
osciló entre 0.409 y 0.459, indicando aparentemente que diferencias en el régimen de lluvia
entre estaciones se manifiestan en valores diferentes de M1/M0.
• Dado lo anterior, es necesario estimar M0 o M1 para diferentes duraciones t. Se decidió en
este trabajo estimar M0 a partir de información pluviométrica. Específicamente se encontró
mediante análisis de regresión múltiple que para cualquier estación y con R2 > 0.91, M0t es
función de M0diario, PA/N y E, donde E es la altitud del lugar en msnm. De esta forma, la
estimación de M0t es:

M 0t = k1 + k 2 M 0 diario + k 3 (PA N ) + k 4 E [6]

donde los coeficientes k1, k2, k3 y k4 se estiman en función de la duración t (en horas) como:

k1 = a1 + b1t + c1t ln t + d1t 0.5 ln t + e1t −1.5 [7]

k 2 = a 2 + b2 t ln t + c2 t 2.5 + d 2 t 0.5 ln t [8]
k 3 = a3 + b3t + c3 ln t + d 3t −0.5 ln t + e3t −1.5 [9]
k 4 = a4 + b4t + c4t 1.5 + d 4t 2 + e4 e − t [10]

En la Tabla 4 se dan los valores de los coeficientes de las ecuaciones 7 a 10.

• Análisis de bondad de ajuste de las curvas PDF: Para cada una de las 35 estaciones se
estimaron sus correspondientes curvas PDF a partir de su correspondiente información
pluviométrica de PA y N, al igual que de su elevación E, y del análisis de frecuencia de sus
series anuales de lluvia diaria máxima (día hidrológico) con Gumbel y MPP para estimar sus
respectivos valores de M0diario y (M1/M0)diario. Mediante la comparación entre las curvas PDF
ajustadas y estimadas se hizo un análisis de bondad de ajuste que indica que el 80% de los
valores estimados de profundidad máxima (para períodos de retorno menores que 50 años)
 tienen errores relativos menores que ±15% para una duración t de 1 hora y menores que
±10% para duraciones mayores. En promedio los errores relativos están alrededor de ±4%
para t = 1 hora y menores a ±2% para duraciones mayores. Adicionalmente, el 80% de los
valores estimados de profundidad máxima tienen errores absolutos menores a ±10 mm para
una duración t de 1 hora y menores a ±12 mm para duraciones mayores. En promedio los
errores absolutos están alrededor de -2 mm y +3 mm para todas las duraciones.

Tabla 3.- Valores de M1/M0 para 8 estaciones

Duración
1107008 1107009 1107010 1601506 1602506 2106503 2115502 2116502
[hr]
1 0.438 0.452 0.444 0.405 0.441 0.414 0.438 0.431
2 0.438 0.456 0.445 0.410 0.435 0.425 0.441 0.431
3 0.446 0.457 0.445 0.413 0.442 0.413 0.438 0.425
4 0.446 0.461 0.447 0.416 0.437 0.414 0.435 0.419
5 0.448 0.461 0.451 0.420 0.436 0.412 0.433 0.415
6 0.449 0.461 0.450 0.419 0.434 0.415 0.432 0.410
7 0.448 0.459 0.452 0.421 0.433 0.416 0.434 0.403
8 0.448 0.459 0.454 0.421 0.434 0.417 0.436 0.401
9 0.448 0.460 0.458 0.423 0.432 0.419 0.437 0.401
10 0.449 0.460 0.458 0.424 0.429 0.421 0.438 0.402
11 0.450 0.461 0.459 0.426 0.429 0.420 0.437 0.400
12 0.450 0.462 0.461 0.428 0.429 0.419 0.438 0.400
13 0.453 0.462 0.461 0.429 0.428 0.419 0.438 0.407
14 0.454 0.463 0.461 0.429 0.428 0.420 0.437 0.410
15 0.454 0.462 0.459 0.430 0.428 0.420 0.431 0.409
16 0.454 0.462 0.456 0.430 0.428 0.420 0.431 0.408
17 0.454 0.461 0.455 0.431 0.428 0.418 0.427 0.408
18 0.454 0.460 0.453 0.431 0.429 0.418 0.424 0.410
19 0.454 0.461 0.452 0.431 0.430 0.417 0.423 0.411
20 0.453 0.457 0.448 0.430 0.432 0.417 0.421 0.412
21 0.454 0.456 0.445 0.428 0.433 0.417 0.422 0.413
22 0.454 0.456 0.442 0.428 0.434 0.415 0.420 0.414
23 0.454 0.456 0.438 0.425 0.443 0.413 0.418 0.412
24 0.453 0.453 0.430 0.423 0.445 0.413 0.418 0.412
36 0.446 0.462 0.427 0.432 0.430 0.422 0.416 0.416
48 0.452 0.456 0.416 0.430 0.427 0.415 0.415 0.423
60 0.454 0.446 0.423 0.424 0.424 0.416 0.419 0.420
Diaria 0.457 0.464 0.451 0.416 0.433 0.420 0.426 0.422

Tabla 4.- Valores de coeficientes de las ecuaciones 7 a 10

Subíndice
Coeficiente
1 2 3 4
a -761.087 0.22288 16.809 -0.0129
b 741.947 -0.07566 0.234 0.004169
c -103.788 0.0005667 -5.8561 -0.001365
d -586.892 0.32844 -17.8728 0.0001313
e 38.82 - 2.3229 0.014595

El análisis de la información para el caso de las 157 curvas IDF disponibles es el siguiente:

Determinación de valores de M0 y M1 para cada estación y cada duración t: con base en las
curvas IDF disponibles, éstas se convirtieron en sus respectivas curvas PDF para las cuales, para
cada duración disponible se determinaron los valores de M0 y M1, es decir M0t y M1t. Además, con
estos valores se calculó la correspondiente relación (M1/M0)t. Las duraciones t más frecuentes
definidas en las curvas IDF disponibles corresponden a 10, 20, 30, 40, 50, 60, 120, 240 y 360
minutos.
Determinación de valores de M0 y M1 diarios para cada estación: con cada serie de lluvias
máximas diarias de cada estación con curvas IDF se realizó el ajuste de la distribución Gumbel con
MPP, obteniendo los correspondientes valores de M0diario, M1diario y (M1/M0)diario. De manera similar
al análisis de las 35 estaciones con información multihoraria, se encontró que en términos generales
en cada estación los valores de la relación (M1/M0)diario se preservan en la relación (M1/M0)t.
Estimación de M0t para t de 60 minutos: se adoptó estimar, para cada estación, M0t para t = 60
minutos a partir de M0diario, M1diario y (M1/M0)diario, PA, N y/o E. Esto se hizo mediante análisis de
regresión mutlivariada, tanto para la totalidad de estaciones como para agrupaciones de estaciones
por zonas hidrográficas. En la Tabla 5 se presentan los coeficientes de estas regresiones.
Estimación de M0t para otras duraciones a partir de valores estimados de M0t para t de 60
minutos: para estimar M0t para duraciones t de 10, 20, 30, 40, 50 y 120 minutos, se encontraron
ecuaciones de regresión de la forma M0t = K1M0,60min (M1/M0)diario + K2. Estas regresiones se
establecieron tanto para la totalidad de las estaciones como para agrupaciones de estaciones por
zonas hidrográficas.
Estimación de curvas IDF para las estaciones disponibles: para cada estación, con base en los
tres puntos anteriores se estimaron las correspondientes curvas PDF e IDF, las cuales fueron
comparadas con las curvas PDF e IDF disponibles, encontrando errores promedio desde
aproximadamente 10% hasta 24%, 22%, 20%, 17% y 24% para duraciones de 10, 20, 30, 60 y 120
minutos, respectivamente.

Tabla 5.- Valores de coeficientes lineales para estimación de M0,60min

Coeficientes αi:
Zona hidrográfica o M0diario-
región M0diario M1diario PA N PA/N E Constante
M1diario
Atrato 0,808 5,936
Sabana Bogotá 0,571 -0,00411 13,544
Catatumbo 0,00723 43,974
Guajira 0,412 0,384 -7,93
Patía, Mira Micay 1,435 -0,01224 45,114
Meta, Guaviare Tomo -1,118 183,954
Sinú, Ciénaga Grande,
0,292 4,443 -14,635
Cesar, Bajo Magdalena
Cuenca Nechí 1,018 -0,272 66,851
Cuenca La Vieja 0,775 -0,01906 12,04
Alto Magdalena
2,326 -0,0142 41,176
(Cundinamarca, Tolima)
Cauca (Valle) -3,049 0,528 0,06263 -2,902 54,805
Magdalena Medio
(Caldas, Antioquia, 3,666 -2,344 3,186 -9,384
Cund., Santander)
Sogamoso 0,208 -0,00228 -7,432
Cauca (Risaralda, Cauca,
2,083 -0,218 2,142 51,984 -539,419
Antioquia)

METODOLOGÍA DE ESTIMACIÓN DE CURVAS IDF

Con base en los análisis para las 35 estaciones con información multihoraria, la metodología
propuesta en este trabajo para la estimación sintética de curvas IDF es la siguiente:

1. Obtener de información diaria de precipitación en una o más estaciones en la zona de

interés, procurando que éstas tengan un registro suficientemente largo (más de 15 años).
Para cada estación i estimar PAi (mm/año) y Ni, y determinar Ei (msnm).
2. Para cada estación, determinar de la correspondiente serie anual de precipitación máxima
diaria, expresando ésta en mm/día.
3. Para cada serie anual de lluvia máxima diaria, ajustar la distribución Gumbel con MPP,
calcular los respectivos valores de M0diario y M1diario con la ecuación 5, y calcular el valor de
(M1/M0)diario.
4. Establecer si es adecuado definir un valor regional de (M1/M0)diario si se ha considerado más
de una estación.
5. Definir el conjunto de valores de t de interés, para t igual o mayor que 1 hora.
6. Para cada valor de t mayor o igual a 1 hora calcular los correspondientes valores de los
coeficientes k1, k2, k3 y k4, mediante las ecuaciones 7 a 10 y la Tabla 4.
 7. Para cada estación i y cada valor de t mayor o igual a 1 hora calcular los valores de M0t
usando la ecuación 6.
8. Para cada estación i y cada valor de t calcular M1t = (M1/M0)diario M0t, usando el valor
regional de (M1/M0)diario o en su defecto los valores puntuales diarios.
9. Con base en las ecuaciones 1 a 3 estimar la curva de frecuencia de profundidad de
precipitación máxima para una duración t y valores seleccionados del período de retorno T.
10. Conformar las curvas PDF para cada estación i.
11. Dividir las curvas PDF por la duración t para finalmente obtener las curvas IDF.

Para el caso de las 157 estaciones con curvas IDF disponibles, el procedimiento de estimación de
éstas es el siguiente (hacer paralelo con procedimiento anterior):

1. Igual.
2. Igual.
3. Igual.
4. Igual.
5. Definir el conjunto de valores de t de interés dentro de los siguientes valores posibles: 10,
20, 30, 40, 50, 60 y 120 minutos.
6. Para t igual a 60 minutos estimar M0t con base en los coeficientes de la Tabla 5, es decir
M0,60min = α1M0diario + α2M1diario + α3(M0diario - M1diario) + α4PA + α5N + α6PA/N + α7E+ α8 .
7. Para los otros valores de t de interés, definidos en el paso 5, estimar los correspondientes
valores de M0t con base en las ecuaciones de la Tabla 6.
8. Igual.
9. Igual.
10. Igual.
11. Igual.

VERIFICACIÓN DE LA METODOLOGÍA PROPUESTA

Para ilustrar la verificación de la metodología propuesta de estimación sintética de curvas PDF

derivada de las 35 estaciones con información multihoraria se seleccionaron 2 estaciones no
utilizadas en los análisis previos pero con el mismo tipo de información multihoraria. Para estas
estaciones se aplicó la metodología propuesta para 1 y 6 horas (estación 1) y 3 y 9 horas (estación
2). En la Tabla 7 se comparan las curvas PDF estimadas con el procedimiento propuesto a partir de
información exclusivamente pluviométrica con aquellas ajustadas a las series anuales horarias con
Gumbel y MPP. Se observa una muy buena representación de las curvas PDF. En el caso de
derivación a partir de curvas IDF para duraciones bajas, en la Tabla 8 se muestra un ejemplo de
buena concordancia para duraciones de 10, 20, 30, 60 y 120 minutos. En general, para este caso la
representación es menor que para las estaciones multihorarias.

CONCLUSIONES Y RECOMENDACIONES

Se ha presentado un procedimiento para estimar curvas IDF y PDF a partir de información

pluviométrica, útil en sitios en donde no existe información pluviográfica. Mediante el ajuste de la
distribución Gumbel con MPP a la serie anual de lluvia máxima diaria (día hidrológico) y el cálculo
de M0 y M1 se captura el comportamiento (pendiente) de las curvas de frecuencia que se mantiene
aproximadamente constante para todas las duraciones. Esto constituye un aspecto novedoso entre
los procedimientos existentes para estimar curvas IDF y PDF de manera sintética. Las relaciones
multivariadas encontradas para estimar M0t a partir de valores derivados de información de lluvia
diaria completa la metodología.
 Tabla 6.- Coeficientes ecuaciones lineales para estimar M0,t
Coeficientes
Zona hidrográfica o región M0,t
(M0diario/ M1diario) * M0,60min Constante
M0,5min 0,041 6,48
M0,10min 0,109 5,26
M0,15min 0,147 4,03
M0,20min 0,211 4,14
Todas M0,30min 0,272 4,01
M0,40min 0,339 2,44
M0,50min 0,388 1,36
M0,90min 0,453 2,571
M0,120min 0,528 -0,271
M0,10min 0,133 3,458
Atrato M0,20min 0,232 1,889
M0,30min 0,305 0,624
M0,10min 0,095 5,717
Sabana Bogotá M0,20min 0,21 3,629
M0,30min 0,277 3,167
M0,10min -0,39 67,78
Catatumbo M0,20min -0,28 64,1
M0,30min -0,099 48,38
M0,10min 0,086 5,826
Guajira M0,20min 0,092 10,06
M0,30min 0,202 10,5
M0,10min 0,119 2,919
Patía, Mira y Micay M0,20min 0,197 3,681
M0,30min 0,278 2,102
M0,10min 0,062 12,51
Meta, Guaviare y Tomo
M0,20min 0,134 14,05
M0,5min 0,168 -13,9
M0,10min 0,093 7,148
Sinú, Ciénaga Grande, Cesar, Bajo Magdalena
M0,20min 0,198 5,25
M0,30min 0,216 11,11
M0,10min 0,173 1,41
Nechí M0,20min 0,307 -1,75
M0,30min 0,418 -6,105
M0,10min -0,026 20,93
La Vieja M0,20min 0,065 20,92
M0,30min 0,154 17,96
M0,10min 0,107 4,221
Alto Magdalena (Cund, Tolima) M0,20min 0,189 5,448
M0,30min 0,277 3,134
M0,10min 0,143 0,823
Cauca (Valle) M0,20min 0,238 0,711
M0,30min 0,306 0,349
M0,10min 0,025 12,74
Sogamoso M0,20min 0,199 5,419
M0,30min 0,301 2,041
M0,10min 0,169 1,677
Magdalena Medio (Caldas, Antioquia, Cund., Santander) M0,20min 0,284 0,12
M0,30min 0,359 -1,353
M0,10min 0,16 0,826
Cauca (Risaralda, Cauca, Antioquia) M0,20min 0,287 -2,317
M0,30min 0,363 -3,57

Los errores de estimación de las curvas PDF son relativamente bajos, lo cual hace que el método
propuesto sea robusto como se ilustra con la estimación de curvas PDF en estaciones diferentes a
las utilizadas en los análisis de derivación. En razón de la sectorización en zonas hidrográficas o
regiones, es evidente la dependencia climatológica de las curvas IDF, en particular para cortas
duraciones. Se considera que la metodología tiene un valor agregado importante, aunque es
recomendable involucrar en los análisis de derivación más estaciones, ojalá de otras regiones del
país.

Tabla 7.- Comparación de curvas PDF reales y estimadas

Estación 2613506 Estación 2620502
M0diario [mm] 80.909 M0diario [mm] 74.367
PA/N 10.343 PA/N 10.4562
E [msnm] 1600 E [msnm] 1600
(M1/M0)diario [mm] 0.43626 (M1/M0)diario [mm] 0.4425
1 hora 6 horas 3 horas 9 horas
Parámetro Parámetro
estimada "real" estimada "real" estimada "real" estimada "real"
M0t [mm] 45.86 46.39 75.48 74.35 M0t [mm] 62.25 66.62 72.72 75.95
M1t [mm] 20.01 20.35 32.93 31.40 M1t [mm] 27.54 28.84 32.18 33.61
(M1/M0)t 0.436 0.439 0.436 0.422 (M1/M0)t 0.443 0.433 0.443 0.443
a 8.43 8.21 13.88 16.66 a 10.33 12.89 12.07 12.60
m [mm] 40.99 41.65 67.47 64.73 m [mm] 56.29 59.18 65.76 68.68
T PDF [mm] PDF [mm] T PDF [mm] PDF [mm]
[años] estimada "real" estimada "real" [años] estimada "real" estimada "real"
1.5 40.2 40.9 66.2 63.2 1.5 55.3 58.0 64.6 67.5
2 44.1 44.7 72.6 70.8 2 60.1 63.9 70.2 73.3
2.33 45.9 46.4 75.5 74.4 2.33 62.3 66.6 72.7 76.0
3 48.6 49.1 80.0 79.8 3 65.6 70.8 76.6 80.1
5 53.6 54.0 88.3 89.7 5 71.8 78.5 83.9 87.6
10 60.0 60.1 98.7 102.2 10 79.5 88.2 92.9 97.0
20 66.0 66.0 108.7 114.2 20 87.0 97.5 101.6 106.1
25 68.0 67.9 111.9 118.0 25 89.3 100.4 104.3 109.0
30 69.5 69.4 114.5 121.1 30 91.2 102.8 106.6 111.3
40 72.0 71.8 118.5 126.0 40 94.3 106.6 110.1 115.0
50 73.9 73.7 121.6 129.7 50 96.6 109.5 112.8 117.8
60 75.5 75.2 124.2 132.8 60 98.5 111.8 115.1 120.1
70 76.8 76.5 126.3 135.4 70 100.1 113.8 116.9 122.1
80 77.9 77.6 128.2 137.6 80 101.5 115.6 118.6 123.8
90 78.9 78.6 129.9 139.6 90 102.7 117.1 120.0 125.3
100 79.8 79.4 131.3 141.4 100 103.8 118.5 121.3 126.6

Tabla 8.- Comparación de curvas IDF reales y estimadas, Estación 1601002

Intensidad [mm/hr]
T
[años] 10 min 20 min 30 min 60 min 120 min
“real” est. “real” est. “real” est. “real” est. “real” est.
10 205,9 194,5 136,3 140,1 104,4 111,1 69,6 75,1 46,4 49,6
25 246,5 232,5 162,4 167,5 124,7 132,8 84,1 89,8 55,1 59,3
50 275,5 260,7 182,7 187,8 142,1 148,9 95,7 100,7 60,9 66,4

AGRADECIMIENTOS

Se agradece a los investigadores del Posgrado en Recursos Hidráulicos de la Universidad Nacional

de Colombia, sede Medellín, por su autorización para utilizar la información horaria de 35
estaciones de medición de precipitación, a la Empresa de Acueducto y Alcantarillado de Bogotá,
EAAB-ESP, por información de estaciones de su jurisdicción, a Cenicafé por permitir tener acceso
a reportes y artículos a través de su página web, y a las firmas consultoras e instituciones que
suministraron directa o indirectamente curvas IDF.

REFERENCIAS

Díaz-Granados, M. y A. Deeb (1988), “Análisis regional de frecuencias usando Momentos Ponderados con
Probabilidad”, VIII Seminario nacional de Hidrálica e Hidrología, Barranquilla, 1988.
Díaz-Granados, M. y X. Puente (2006), “Estimación de curvas Intensidad-Duración-Frecuencia a partir de
información pluviométrica”, XVII Seminario Nacional de Hidráulica e Hidrología, Popayán, 2006.
Froehlich, D. C. (1995), “Intermediate Duration – Rainfall – Intensity Equations”, Journal of Hydraulic
Engineering, 121, 1995. pp. 751-756.
Greenwood, J. A., J. M. Landwehr, N. C. Matalas y J. R. Wallis (1979), “Probability Weighted
Moments: Definition and relation to parameters of several distributions expressable in inverse form”, Water
Resources Research, Vol. 15, No. 5, Octubre 1979, pp. 1049 – 1054.
Greis, N. P. y E. F. Wood (1981), “Regional Flood Frequency Estimation and Network Design”, Water
Resources Research, Vol. 17, No. 4, Agosto 1981, pp. 1167 – 1177.
Kothyari, U. y R. Garde (1992), “Rainfall Intensity_Duration-Frequency Formula for India”, Journal of
Hydraulic Engineering, Febrero 1992, Vol. 118, No. 2.
Landwehr, J. M., N. C. Matalas y J. R. Wallis (1979), “Probability Weighted Moments compared with
some tradicional techniques in estimating Gumbel parameters and quantiles”, Water Resources Research,
Vol. 15, No. 5, Octubre 1979, pp. 1055 – 1064.
Powell, R. W. (1932), Discussion of “Formulas for Rainfall Intensities of Long Durations” de M. M.
Bernard, Trans., ASCE, 96, 1932.
Vargas, R. y M. Díaz-Granados (1998), “Curvas Sintéticas Regionalizadas de Intensidad – Duración –
Frecuencia para Colombia”, XIII Seminario Nacional de Hidráulica e Hidrología, Cali, Agosto de 1998.
Vélez, J., G. Poveda, O. Mesa, L. Salazar, S. Vieira, J. Mejía, C. Hoyos y D. Quevedo (2002),
“Estimación de Curvas Intensidad – Frecuencia – Duración en Colombia Aplicando Diferentes
Metodologías”, XV Seminario Nacional de Hidráulica e Hidrología, Medellín, Agosto de 2002.