0% encontró este documento útil (0 votos)
32 vistas11 páginas

Intensidad de Lluvias y Recurrencia

Cargado por

agustina
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
32 vistas11 páginas

Intensidad de Lluvias y Recurrencia

Cargado por

agustina
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

1

H0513 – HIDROLOGÍA I Curso 2020

Trabajo Práctico N° 4

RÉGIMEN DE LAS LLUVIAS

Ejercicio 1: Graficar las intensidades máximas I en función de su recurrencia


asociada R, para cada muestra de duración constante (d=cte).

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.

A partir de la recopilación y el análisis de datos de lluvias reales se estableces leyes I-R-d,


que gráficamente se representan mediante familias de curvas y muestran la forma en que se
relacionan dichas variables. En este ejercicio las curvas obtenidas son del tipo I=f(R), las cuales
permiten visualizar cómo varían intensidad respecto a la recurrencia para cada duración.

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]

mayor recurrencia, lo que 200


implica que su probabilidad
de ocurrencia es baja. Por 150
otro lado, se observa que
cuanto más duraderas sean
100
las tormentas, más
recurrentes serán y a la vez,
50
su intensidad será menor.

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'

Agustina Ferrari – 64196/0


2

relación de la Intensidad de la lluvia en función de la Recurrencia para duración


constante y superponer la recta de ajuste en el gráfico obtenido en el numeral 1).

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

Para las componentes de la matriz de ajuste que 0.720 0.519


refieren a las intensidades sí debe hacerse un 0.623 0.388
cálculo para cada duración ya que para una
misma recurrencia R, las intensidades I van 0.544 0.296
disminuyendo a medida que la duración de la 0.477 0.228
tormenta aumenta.
0.419 0.176
0.368 0.135
0.322 0.104

[ ∑ 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”

Agustina Ferrari – 64196/0


3

5 minutos 10 minutos 15 minutos


∑y y*x ∑y*x ∑y y*x ∑y*x ∑y y*x ∑y*x

∑I I * (log R) ∑ I (log R) ∑I I * (log R) ∑ I (log R) ∑I I * (log R) ∑ I (log R)

312.705 275.683 229.141


210.467 175.338 149.196
151.864 131.582 107.919
123.651 105.791 85.627
103.023 87.567 70.614
86.888 73.830 59.031
73.524 61.692 49.477
61.905 52.182 41.787
51.701 43.752 34.921
43.467 36.507 29.548
2827.8 1384.966 2374.8 1179.933 1950.8 970.986
36.536 29.796 24.853
30.258 24.814 20.926
25.264 20.598 17.433
20.515 17.134 14.175
16.790 13.634 11.530
13.003 10.735 8.940
9.810 8.122 6.727
7.009 5.764 4.747
4.473 3.673 2.977
2.113 1.740 1.418

20 minutos 30 minutos 45 minutos


∑y y*x ∑y*x ∑y y*x ∑y*x ∑y y*x ∑y*x

∑I I * (log R) ∑ I (log R) ∑I I * (log R) ∑ I (log R) ∑I I * (log R) ∑ I (log R)


193.176 159.327 128.652
125.504 103.446 85.474
95.750 77.495 61.270
75.833 60.277 47.603
61.639 48.676 39.202
50.326 40.315 32.481
41.462 32.683 26.957
34.411 27.369 22.004
28.739 23.146 18.730
23.876 19.236 16.143
1639.6 823.298 1304.4 666.183 1041.2 536.873
20.135 16.316 13.171
16.843 13.926 11.010
13.934 11.434 8.935
11.358 9.245 7.413
9.308 7.555 5.801
7.428 5.822 4.405
5.699 4.231 3.203
4.050 2.926 2.283
2.595 1.865 1.447
1.231 0.892 0.689

Agustina Ferrari – 64196/0


4

60 minutos 90 minutos 120 minutos


∑y y*x ∑y*x ∑y y*x ∑y*x ∑y y*x ∑y*x

∑I I * (log R) ∑ I (log R) ∑I I * (log R) ∑ I (log R) ∑I I * (log R) ∑ I (log R)

106.439 87.399 71.532


69.135 57.697 44.626
51.128 38.283 32.198
40.689 30.895 25.134
31.723 26.239 19.757
25.952 21.164 16.376
21.232 16.270 13.216
17.981 13.622 10.939
14.609 11.665 8.721
11.245 9.957 6.863
826.8 435.139 650.8 348.420 499.8 274.546
9.801 8.453 5.532
8.482 6.927 4.788
6.936 5.269 3.936
5.582 4.314 3.187
4.632 3.346 2.528
3.555 2.516 1.866
2.615 1.881 1.422
1.801 1.319 0.998
1.100 0.822 0.626
0.502 0.384 0.299

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.

Agustina Ferrari – 64196/0


5

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

Intensidad vs. Recurrencia


250

f(x) = 45.0641608824342 ln(x) + 99.5821026184721 d


=
5'
200 f(x) = 40.9334504183819 ln(x) + 80.7643346141913
L
Intensidad de lluvias [mm/h]

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]

Agustina Ferrari – 64196/0


6

Ejercicio 3: Obtener las expresiones analíticas de cada una de las relaciones


I=f(R); d=cte.

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

Agustina Ferrari – 64196/0


7

Ejercicio 4: Estimar los valores de intensidad máxima de lluvia para 2, 5 y 10


años de Recurrencia.

Estimación de intensidades máximas [mm/h]


Duración Recurrencia [años]
Expresión analítica calculada
[minutos] 2 5 10
5 I = 103.76 * log (R) + 99.58 130.8 172.1 203.3
10 I = 94.25 * log (R) + 80.76 109.1 146.6 175.0
15 I = 78.15 * log (R) + 66.05 89.6 120.7 144.2
20 I = 68.73 * log (R) + 54.29 75.0 102.3 123.0
30 I = 59.41 * log (R) + 41.28 59.2 82.8 100.7
45 I = 49.58 * log (R) + 32.08 47.0 66.7 81.7
60 I = 43.10 * log (R) + 23.98 37.0 54.1 67.1
90 I = 36.42 * log (R) + 17.87 28.8 43.3 54.3
120 I = 30.91 * log (R) + 12.53 21.8 34.1 43.4

Agustina Ferrari – 64196/0


8

Ejercicio 5: Graficar los valores, obtenidos en el numeral 4), de intensidad


máxima de lluvia en función de la duración para Recurrencia constante (I=f (d);
R=cte.)

Intensidad vs. Duración


250

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

Intensidad vs. Duración


1000
Intensidad [mm/h]

100

10
1 10 100 1000
Duración [minutos]

R=2 R=5 R=10

Agustina Ferrari – 64196/0


9

Ejercicio 6: Ajustar los puntos graficados en 5) por mínimos cuadrados,


superponer la recta de ajuste en el mismo gráfico.

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

log ( Y +C )=log (a)+b [ log ( x ) ] → Y +C=a X


b

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

de donde se despejan los valores de “A” y “b” a través de la resolución de un sistema de


ecuaciones:

N * A + (∑X) * b = ∑Y

(∑X) * A + (∑(X)^2) * b = ∑ YX

Por último, se definen las equivalencias de las Datos generales de la matriz


variables que serán:
∑X= X^2 = ∑ X^2 =
X = log x
∑ log x (log x)^2 ∑(log x)^2
→ Y = I + C = log (I+C)
log d ∑ log d (log d)^2 ∑(log d)^2

→ A = log a 0.699 0.489


1.000 1.000
1.176 1.383
→ b=b 1.301 1.693
1.477 13.118 2.182 20.783
1.653 2.733
→ X = log (d) 1.778 3.162
1.954 3.819
2.079 4.323

Agustina Ferrari – 64196/0


10

R=2 años R=5 años R=10 años


∑Y=∑ YCoeficientes
= log ∑ Y = ∑ Y = log ∑Ajuste
Y=∑
Y = log y X*Y ∑ X*Y X*Y ∑ X*Y X*Y ∑ X*Y
Recurrencia
log y y log y y log y
(log b ∑ (log A (loga ∑ (log Y=AX^b // log y=log a(log
+ b log x∑ (log
log I ∑2 log I d)*(log-0.577
d)*(log log
2.595I ∑ log I d)*(log 393.83
d)*(log log I ∑ log I d)*(log d)*(log
I) I) I) I)
y=393.830*x^(-0.577) I) I)
2.117 5 1.479-0.524 2.668
2.236 1.563 465.81 y=465.809*x^(-0.524)
2.308 1.613
2.038 10 2.038-0.501 2.721
2.166 2.166 525.96 y=525.956*x^(-0.501)
2.243 2.243
1.952 2.296 2.082 2.448 2.159 2.539
1.875 2.439 2.010 2.615 2.090 2.719
Duración Recurrencia [años] 25 Valor de C
1.772 Expresión2.618
15.793 analítica22.060
calculada
1.918 17.139 2.833 24.110 2.003 17.914 2.959 25.278
[minutos] 2 5 10 2 5 10
5
1.672 2.764
I = 103.76 * log (R) + 99.58
Intensidad
1.824
130.8
vs3.016
Duración
172.1
1.912
203.3 155.8
3.161
197.1 228.3
1.568 2.788 1.733 3.082 1.827 3.248
10 I = 94.25 * log (R) + 80.76 109.1 146.6 175.0 134.1 171.6 200.0
1000.0
1.460 2.853 1.637 3.199 1.735 3.390
15 I = 78.15 * log (R) + 66.05 89.6 120.7 144.2 114.6 145.7 169.2
1.339 2.784 1.533 3.188 1.638 3.405
20 I = 68.73 * log (R) + 54.29 75.0 102.3 123.0 100.0 127.3 148.0
30 I = 59.41 * log (R) + 41.28 59.2 82.8 100.7 84.2 107.8 125.7
45 I = 49.58 * log (R) + 32.08 47.0 66.7 81.7 72.0 91.7 106.7
60 I = 43.10 * log (R) + 23.98 37.0 54.1 67.1 62.0 79.1 92.1
90 I = 36.42 * log (R) + 17.87 28.8 43.3 54.3 53.8 68.3 79.3
120 f(x) = 472.330990425903
I = 30.91 * log (R) + 12.53 21.8 x^-0.395172088612363
34.1 43.4 46.8 59.1 68.4
R² ==0.988618924373446
f(x) 405.996166789129 x^-0.39482025785008
R² = 0.989456503440022 c=0 0.9766 0.9796 0.9806
f(x) = 318.303538336512 x^-0.394133511565913
R² = 0.990900419742075 c=10 0.9853 0.985 0.9848
Intensidad [mm/h]

c=15 0.9877 0.9869 0.9863


100.0
2
R c=20 0.9896 0.9883 0.9876
c=25 0.9909 0.9895 0.9886
c=30 0.9919 0.9904 0.9895
c=35 0.9927 0.9911 0.9902

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).

Agustina Ferrari – 64196/0


11

Ejercicio 7: Obtener las expresiones analíticas de I=f (d); R=cte.

Recurrencia Expresiones analíticas para I = f (d)


[años] [mm/h]

2 I + 25 = (318.3)d^(-0.394)

5 I + 25 = (406)d^(-0.395)

10 I + 25 = (472.33)d^(-0.395)

Agustina Ferrari – 64196/0

También podría gustarte