Intensidad de Lluvias y Recurrencia
Intensidad de Lluvias y Recurrencia
Trabajo Práctico N° 4
El régimen de lluvias consiste en el estudio de las relaciones existentes entre las tres variables
que rigen el fenómeno que son:
Intensidad “I” [mm/h]: la tasa de la caída (cantidad) de lluvia que se mide en un intervalo
de tiempo.
Recurrencia “R” [años]: el intervalo de tiempo medio entre dos lluvias que alcanzan o
superan un valor extremo de intensidad considerado.
Duración “d” [minutos]: el intervalo de tiempo durante el que se produce la máxima
intensidad media de lluvias.
En el gráfico obtenido
Intensidad vs. Recurrencia
se puede ver que las lluvias 250
más intensas son aquellas de
Intensidad de lluvias [mm/h]
0
1 10 100
Ejercicio 2: Ajustar por Recurrencia [años]
mínimos cuadrados la d=5' d=10' d=15' d=20' d=30' d=45' d=60' d=90' d=120'
Para obtener una recta de ajuste para la ley I=f(R) se ordenan los valores de intensidad
máxima de forma decreciente para cada duración. En cada caso se utilizará el método de ajuste por
mínimos cuadrados, y resolviendo la matriz se obtendrán los valores de pendiente y ordenada de
una recta que se ajustará a la variación de los datos del registro pluviográfico.
El ajuste se realizará por medio de ecuaciones del tipo y=a+b [ log ( x ) ] . Con el método de
mínimos cuadrados las incógnitas a despejar serán "a" y "b". Éstas representarán el comportamiento
general de la variación de la intensidad de lluvias en función de la recurrencia para cada intervalo
de duración particular. La equivalencia de variables será entonces:
y=I
log (x) = log (R)
x ∑x x^2 ∑ x^2
Los valores referidos únicamente a la
recurrencia R se emplean en el ajuste de todas ∑ (log
log R ∑ log R (log R)^2
las rectas porque son los mismos para cada R)^2
duración. 1.322 1.748
1.021 1.043
0.845 0.714
[ ∑ log R a
]
8.058 5.614
N 0.281 0.079
=
∑ log R ∑ log R 2 b
Matriz para resolver el ajuste lineal
0.243 0.059
0.208 0.043
∑I 0.176 0.031
∑ log R∗I 0.146 0.021
0.118 0.014
0.092 0.008
0.067 0.004
N * a + (∑ log R) * b = ∑ I 0.043 0.002
(∑ log R) * ade
+ (∑ log R)^2a *resolver
b = ∑ Idel
(logcual
R) 0.021 0.000
Sistema ecuaciones
se despejarán los coeficientes “a” y “b”
Coeficientes Función obtenida por Funciones logarítmicas Los valores de “a” obtenidos
Duración obtenidas mediante la línea de por mínimos cuadrados son
b a mínimos cuadrados tendencia de Excel
iguales a los que arroja la
5 y = 45.064 ln(x) + 99.582 función logarítmica de Excel,
minutos
103.76 99.58 I=103.76logR+99.58
pero no sucede lo mismo
10 y = 40.933 ln(x) + 80.764 con los valores del
minutos
94.25 80.76 I=94.25logR+80.76 coeficiente “b”.
Sin embargo, existe una
15 y = 33.941 ln(x) + 66.052 relación constante entre
minutos
78.15 66.05 I=78.15logR+66.05
ellos para cada duración y
20 y = 29.849 ln(x) + 54.288 tras analizar el
minutos
68.73 54.29 I=68.73logR+54.29 comportamiento entre las
funciones log10(x) y ln(x), se
30 y = 25.802 ln(x) + 41.283 concluye que las pendientes
minutos
59.41 41.28 I=59.41logR+41.28
son equivalentes y ambas
45 y = 21.533 ln(x) + 32.083 dan el mismo resultado “y=I”
minutos
49.58 32.08 I=49.58logR+32.08 para cada valor de “R=x”.
Por esta razón puede
60 y = 18.717 ln(x) + 23.976 decirse que las rectas
minutos
43.10 23.98 I=43.10logR+23.98
obtenidas por mínimos
90 y = 15.817 ln(x) + 17.866 cuadrados son correctas.
minutos
36.42 17.87 I=36.42logR+17.87
La razón constante entre los
120 y = 13.425 ln(x) + 12.535 coeficientes “b” en cada
minutos
30.91 12.53 I=30.91logR+12.53
caso es 2.3.
FUNCIONES I = f ( R )
Recurr
5' 10' 15' 20' 30' 45' 60' 90' 120'
encia T
[años] I=103.76log I=94.25logR I=78.15logR I=68.73logR I=59.41logR I=49.58logR I=43.10logR I=36.42logR I=30.91logR
R+99.58 +80.76 +66.05 +54.29 +41.28 +32.08 +23.98 +17.87 +12.53
21.00 236.8 205.4 169.4 145.2 119.8 97.6 81.0 66.0 53.4
10.50 205.5 177.0 145.9 124.5 102.0 82.7 68.0 55.1 44.1
7.00 187.3 160.4 132.1 112.4 91.5 74.0 60.4 48.7 38.7
5.25 174.3 148.6 122.3 103.8 84.1 67.8 55.0 44.1 34.8
4.20 164.3 139.5 114.8 97.1 78.3 63.0 50.9 40.6 31.8
3.50 156.0 132.0 108.6 91.7 73.6 59.1 47.4 37.7 29.4
3.00 149.1 125.7 103.3 87.1 69.6 55.7 44.6 35.3 27.3
2.63 143.1 120.3 98.8 83.1 66.2 52.9 42.1 33.1 25.5
2.33 137.8 115.5 94.8 79.6 63.2 50.3 39.8 31.3 23.9
2.10 133.0 111.1 91.2 76.4 60.4 48.1 37.9 29.6 22.5
1.91 128.7 107.2 88.0 73.6 58.0 46.0 36.1 28.1 21.2
1.75 124.8 103.7 85.1 71.0 55.7 44.1 34.5 26.7 20.1
1.62 121.2 100.4 82.3 68.6 53.7 42.4 33.0 25.5 19.0
1.50 117.9 97.4 79.8 66.4 51.8 40.8 31.6 24.3 18.0
1.40 114.8 94.5 77.5 64.3 50.0 39.3 30.3 23.2 17.1
1.31 111.8 91.9 75.3 62.4 48.3 37.9 29.1 22.2 16.2
1.24 109.1 89.4 73.2 60.6 46.7 36.6 27.9 21.2 15.4
1.17 106.5 87.1 71.3 58.9 45.3 35.4 26.9 20.3 14.6
1.11 104.1 84.9 69.5 57.3 43.9 34.2 25.9 19.5 13.9
1.05 101.8 82.8 67.7 55.8 42.5 33.1 24.9 18.7 13.2
o
f(x) = 33.9405959669019 ln(x) + 66.0518966942698 ga
150 ri
f(x) = 29.848718835971 ln(x) + 54.2881013908153 th
m
ic
f(x) = 25.8015748011556 ln(x) + 41.2828053325408 (d
100 =
f(x) = 21.5331034907185 ln(x) + 32.0828437983254 5'
)
f(x) = 18.7167590417276 ln(x) + 23.9756837449435
f(x) = 15.8166752887782 ln(x) + 17.8662121286287 d
50
f(x) = 13.4251332082738 ln(x) + 12.5349452399848 =
1
0'
0
1.00 10.00 100.00
Recurrencia [años]
Duración
Expresión analítica calculada
[minutos]
5 I = 103.76 * log (R) + 99.58
10 I = 94.25 * log (R) + 80.76
15 I = 78.15 * log (R) + 66.05
20 I = 68.73 * log (R) + 54.29
30 I = 59.41 * log (R) + 41.28
45 I = 49.58 * log (R) + 32.08
60 I = 43.10 * log (R) + 23.98
90 I = 36.42 * log (R) + 17.87
120 I = 30.91 * log (R) + 12.53
200
Intensidad [mm/h]
150
100
50
0
0 20 40 60 80 100 120 140
Duración [minutos]
R=2 R=5 R=10
100
10
1 10 100 1000
Duración [minutos]
En este caso, al tratarse de un ajuste que se adaptará a una ley de recurrencia constante tal que
I= f (d), la ecuación que se deberá utilizar será una del tipo
donde "C" es una constante que sirve para alinear los puntos que corresponden a los mayores
valores de Y. La resolución por ajuste de curvas por mínimos cuadrados se fundamenta en que para
una ecuación de forma potencial puede aplicar una linealización de manera que:
b
y=a ∙ x (1)
ln ( y ) =ln ( a ) +b ∙ ln (x) (2)
Y = A +b X (3)
Luego la ecuación (3) es la que se utiliza para formar la matriz del procedimiento que será:
(∑NX ∑X A
2
∑X b )( ) ( )
= ∑Y
∑ XY
N * A + (∑X) * b = ∑Y
(∑X) * A + (∑(X)^2) * b = ∑ YX
10.0
1 10 100 1000
Duración [minutos]
R=2 Power (R=2) R=5 Power (R=5) R=10 Power (R=10)
Las expresiones encontradas se corresponden con el caso de C=0. Para encontrar qué valor
de C ajusta mejor a los datos, se procederá a utilizar la línea de tendencia proporcionada por Excel,
la cual a través de R2 nos indicará los nuevos coeficientes de la función potencial para cuando se
acerque lo más posible a 1, sin perder de vista en el gráfico el acercamiento entre los datos (esto es
para no perder precisión).
2 I + 25 = (318.3)d^(-0.394)
5 I + 25 = (406)d^(-0.395)
10 I + 25 = (472.33)d^(-0.395)