Proyecto de Exposicion para Ingenieria Civil
Proyecto de Exposicion para Ingenieria Civil
Proyecto de Exposicion para Ingenieria Civil
FACULTAD DE CIENCIAS Y TECNOLOGIA CARRERA DE INGENIERIA CIVIL DEPARTAMENTO DE OBRAS HIDRAULICAS Y SANITARIAS
ESTIMACION Y ACTUALIZACION DE LAS ECUACIONES INTENSIDAD - DURACION - FRECUENCIA PARA LAS PRINCIPALES CIUDADES DEL DEPARTAMENTO DE TARIJA
INTRODUCCION AL PROYECTO ESTIMACIN Y ACTUALIZACION DE ECUACIONES IDF PARA LAS PRINCIPALES CIUDADES Y POBLACIONES DEL DEPARTAMENTO DE TARIJA
INTRODUCCION
MOTIVO DE LA REALIZACIN DEL PRESENTE PROYECTO ?
El motivo principal radica que en nuestro medio no se cuenta con datos de caudales aforados en cantidad y calidad deseable para la realizacin de proyectos, para esto se debe recurrir a mtodos indirectos; estos mtodos requieren de intensidades de lluvia asociados a una probabilidad de ocurrencia (tormenta de diseo) que se pueden determinar a partir de mediciones de precipitaciones. Por tal motivo se busca determinar un parmetro importante que es la intensidad de la lluvia por medio de modelos matemticos que relacionan la Intensidad, la Duracin y la Frecuencia. Y con las ecuaciones estimadas, se dejara al buen criterio del diseador el empleo de las mismas en la bsqueda de parmetros de diseo, como es el caudal, que es bsico para dimensionar las estructuras encargadas de transportar y retener el recurso agua, de tal forma prevenir efectos indeseables.
INTRODUCCION
MOTIVO DE PROYECTO ? LA REALIZACIN DEL PRESENTE
Y como sabemos que es mas complicado determinar las intensidades de duracin corta, es por eso que en el presente proyecto por medio de la bibliografa consultada se encontr un modelo que relaciona la intensidad con la duracin llamado Talbot, en el cual se ha demostrado que da buenos resultados para lluvias de 5 min a 20 min de duracin. Por tal motivo se determinaron los coeficientes del modelo para las once ciudades del departamento en estudio, Y mediante estos coeficientes se determinara la intensidad que es un parmetro importante en el diseo urbano, como por ejemplo el dimensionamiento de un sistema de alcantarillado pluvial.
INTRODUCCION
QUE ES UNA TORMENTA DE DISEO?
Se entiende por tormenta al conjunto de lluvias que obedecen a una misma perturbacin meteorolgica y de caractersticas bien definidas. De acuerdo a esta definicin una tormenta puede durar desde unos pocos minutos hasta varias horas y an das; pueden abarcar extensiones de terrenos muy variables, desde pequeas zonas hasta extensas regiones. Tambin una tormenta de diseo puede definirse mediante un valor de profundidad de precipitacin en un punto, mediante un hietograma de diseo que especifique la distribucin temporal de la precipitacin durante una tormenta, o mediante un mapa de isoyetas que especifique el patrn espacial de la precipitacin, entonces la tormenta de diseo es la entrada al sistema, y los caudales resultantes se calculan utilizando procedimientos de lluvia-escorrenta y trnsito de caudales. Las tormentas de diseo pueden basarse en informacin histrica de precipitacin en un sitio o pueden construirse utilizando las caractersticas generales de la precipitacin en regiones adyacentes.
INTRODUCCION
OBJETIVO GENERAL DEL PROYECTO
Determinar los coeficientes de las ecuaciones Intensidad - Duracin Frecuencia para las principales ciudades y poblaciones del departamento de Tarija con el fin de estimar un parmetro fundamental en el diseo de las obras hidrulicas por ejemplo; en el sistema alcantarillado pluvial, sumideros y cunetas. OBJETIVOS ESPECIFICOS DEL PROYECTO Construccin de las curvas Intensidad Duracin Frecuencia, para las principales ciudades del departamento estudiadas en cuestin. Validar las metodologas propuestas mediante la interpretacin de las bandas pluviograficas en la ciudad de Tarija. Comparar los resultados de las ecuaciones Intensidad Duracin Frecuencia, mediante las ecuaciones como las del alcantarillado pluvial (En ciudades que hubieran) y con la formula de ENDE para as dar la validez a los resultados. Crear un software que permita calcular los coeficientes de las metodologas propuestas en el presente proyecto Estimacin y Actualizacin de ecuaciones Intensidad Duracin Frecuencia para las principales ciudades del departamento de Tarija. Generar tablas de ecuaciones de Intensidad Duracin Frecuencia, en funcin del periodo de retorno en la metodologa de Talbot. Proponer ecuaciones de modelo de intensidad en funcin de la duracin de la lluvia y el periodo de retorno.
PRECIPITACIN.- La precipitacin, es toda forma de humedad que originndose en las nubes, llega hasta la superficie terrestre. La precipitacin incluye la lluvia, la nieve y otros procesos mediante los cuales el agua cae a la superficie terrestre, tales como el granizo y nevisca. INTENSIDAD.- La intensidad, es la cantidad de agua cada por unidad de tiempo. De las tormentas interesa la intensidad mxima que se haya presentado. Es decir, la altura mxima de agua cada por unidad de tiempo PERIODO DE RETORNO.- El periodo de retorno se define como el intervalo promedio de tiempo dentro del cual un evento de magnitud dada puede ser igualada o excedida por lo menos una vez en promedio DURACIN DE LA LLUVIA.- La duracin, es el tiempo que transcurre entre el comienzo y el fin de la tormenta, que vara segn la distribucin espacial.
PLUVIOGRAFO
Dnde: Xi = valor a asumir por la variable aleatoria. e = base de los logaritmos neperianos.
Dnde: Dc = supremo de las diferencias entre la frecuencia observada y terica acumulada. Fn (X)i = frecuencia observada acumulada. F(X)i = frecuencia terica acumulada por Gumbell. Una vez obtenido el valor de Dc, se recurre a la tabla de valores crticos de D en la prueba de bondad de ajuste de Kolmogorov Smirnov, con un tamao de muestra N y asumiendo un 95% de confianza. Luego se decide si el ajuste es o no adecuado, considerando las siguientes afirmaciones: Si el valor Dcalculado < Dtabla, se acepta y se aplica el modelo (el ajuste es adecuado). Si el valor Dcalculado > Dtabla, se rechaza el modelo y se busca otro modelo(el ajuste no es adecuado).
Dnde: n = Nmero de registros S = ndice de desviacin calculado Si= Nmero de valores de xj>xi para i< j <n Ti= Nmero de valores de xj<xi para i< j <n Luego se elige un nivel de significancia o valor de confiabilidad en funcin al cual se definir la condicin de homogeneidad de la serie. Este ndice se relaciona con un valor de Vcrit a travs de la funcin de distribucin normal, que se muestra a continuacin en la siguiente Tabla: f
Caso A: La serie de puntos encaja perfectamente en una lnea recta, lo que indica proporcionalidad, y por lo tanto, la estacin que se analiza es consistente. Caso B: Series de rectas paralelas. Lo cual nos indica proporcionalidad, aunque existan aos que estn medidos por exceso o defecto. Caso C: Cuando se forman dos rectas de diferentes pendientes, se tiene un caso tpico de error sistemtico. La correccin se realiza por la relacin de pendientes del tramo ms antiguo ya que la experiencia demuestra en un 80% el periodo ms moderno es el correcto. Caso D: La estacin presenta un tramo central de mayor o menor pendiente; en el 95 % de los casos, dicho tramo se midi incorrectamente, por lo que habr que corregirlo para homogeneizar la serie.
El Puente Zona Andina Yunchara San Lorenzo Cercado Zona del Valle Central San Lorenzo Uriondo Padcaya Bermejo Zona Sub-Andina Padcaya Entre Ros Yacuiba Zona Chaquea Carapar Villamontes
2da Seccin Prov. Mndez 3ra Seccin Prov. Avils 1ra Seccin Prov. Mndez Prov. Cercado 1ra Seccin Prov. Mndez 1ra Seccin Prov. Avils 1ra Seccin Prov. Arce 2da Seccin Prov. Arce 1ra Seccin Prov. Arce Prov. O'Connor 1ra Seccin Prov. Gran Chaco 2da Seccin Prov. Gran Chaco 3ra Seccin Prov. Gran Chaco
Precipitacin Anual :
Ao ms lluvioso : 1945 con Panual : 960 mm Ao menos lluvioso : 1983 con Panual : 315 mm Temperaturas Ambiente (media aritmtica) Verano : Invierno : 17.7 C 18.6 C 39.3 C en octubre de Temperatura Mxima Extrema : 1977 Temperatura Mnima Extrema : d-7.8 C en Julio de 1975 Los vientos predominantes son los Suroeste con 4 nudos 20.5 C 13.7 C Otoo : Primavera :
Siendo la varianza que caracteriza la zona pluviomtrica, donde se agrupan los k pluvimetros. Si las k estaciones pertenecen a una misma zona pluviomtrica, los k puntos definidos por las coordenadas (ni; S^2(h)) debern recaer en el interior de la franja de control. Verificando dicha prueba se podr admitir, desde el punto de vista estadstico que, se cumple que s'^2(h) = s^2(h) y que las variaciones de S^2(h) respecto de que s'^2(h), se deben a errores de muestreo, aceptndose, entonces que s(h) es constante en toda la zona.
hdt = Altura de lluvia Mxima Diaria Ed = Moda Ponderada Kd = Caracterstica Ponderada T = Periodo de Retorno htT = Altura de lluvia Mxima Horaria t= Tiempo de Duracin de la Lluvia = Constante que en nuestro medio se adopta generalmente 0.2 = Equivalente de lluvia Diaria que depende de la Magnitud de la cuenca
Para
Ac < 20 km^2
= 2 hrs
El coeficiente en si es la duracin mxima que tubo lluvia en un da, por lo tanto este valor se lo puede medir en funcin a un registro de datos pluviogrficos que se realiza diariamente. Segn lo averiguado en el SENAMHI se constat que en la
ciudad de Tarija se tubo un evento extraordinario que ocurrio el 2 de febrero del 2000 con una lluvia concentrada que duro
aproximadamente de 8 a 9 hrs donde debido a evento ocurri el colapso el puente de Santa Ana. Bien pero tambin nos informaron que la lluvia puede alcanzar una duracin de hasta 18 hrs pero considerando las lluvias bajas como las lloviznas y no a una verdadera tormenta. Pero generalmente una tormenta aqu en la ciudad de Tarija dura de 2 a 4 hrs. El equivalente de lluvia diaria depende tambin de las caractersticas de la precipitacin en la zona, que puede influir en la duracin mxima de una lluvia en la zona. Por lo tanto en funcin a la bibliografa consultada y por lo averiguado por el SENAMHI se considero el valor de 12 hrs para como el critico. Ahora bien respetando el criterio de otros ingenieros que optan por adoptar como valor de =18 hrs, =22 hrs y inclusive =24 hrs. Donde este criterio no es compartido por mi
persona porque una lluvia nunca duro 22 hrs de forma continua. Pero algo importante que cabe notar, es que mayor valor de la altura de lluvia horaria es menor y viceversa, esto se debe a que es inversamente proporcional por lo tanto se respeta el valor de adoptado en el proyecto.
Dnde: Hp= es la altura o lmina de lluvia registrada, en cm o mm. d= es la duracin de la tormenta, en horas. c= es un coeficiente; y n es un exponente menor que la unidad.
Dnde: =c*n; y =1n. Este modelo de ecuacin intensidad se lo calcula para determinado periodo de retorno y su procedimiento para calcular dicha ecuacin se muestra a continuacin:
Se realiza un cambio de variable parecido para los coeficientes a calibrar para lo cual se tiene la siguiente expresin:
Dnde: A = log a; B = b; Ahora podemos observar que la expresin inicial se volvi lineal por lo cual para determinar los coeficientes A y B se resuelve el siguiente sistema de mnimos cuadrados:
Una vez resuelto el sistema de ecuaciones se debe reemplazar el cambio de variable hecho anteriormente para as hallar los verdaderos coeficientes es decir: Y de esta forma se hallan los coeficientes a y b que pertenecen a una familia de curva. MODELO DE TALBOT.- Otro modelo que interrelaciona la intensidad de la lluvia y la duracin ha sido desarrollado por Talbot (1989) y su estructura general est representada por: Dnde: i= es la intensidad de la lluvia; d= es la duracin de la tormenta; Y son parmetros que se determinan por anlisis de regresin.
I = Intensidad mxima en mm/hora a y b = Los parmetros a y b se obtienen aplicando una regresin por aproximaciones. D = Duracin de la precipitacin en minutos.
Y volviendo al cambio de variable obtendremos de esta forma los coeficientes para la ecuacin Intensidad Duracin Frecuencia para la ciudad en estudio en funcin del periodo de retorno determinado. g
simultneamente la intensidad, la duracin y el perodo de retorno con una familia de curvas representada por la ecuacin siguiente:
MODELO
PROPUESTO
POR
BERNARD.-
Relaciona
Dnde: i= Es la intensidad de la lluvia, en mm/h; T = Periodo de retorno en aos; d = Es la duracin de la tormenta, en min; , y son parmetros que se determinan al ajustar los datos de la ecuacin.
MODELO DE SHERMAN: Este modelo calcula un nuevo coeficiente llamado que busca ajustarse de mejor manera a la familia de curvas.
Donde: n es el numero de los sumandos. Resolviendo y reemplazando el cambio de variable tenemos: de la Grafica
Dnde: I = Intensidad en (mm/hrs). T = Periodo de retorno en (aos). t = Duracin de la lluvia en (min). La ecuacin mostrada fue elaborada por un grupo de ingenieros hidrlogos en la ciudad de La Paz en un estudio hecho a nivel nacional denominado Estimacin de ecuaciones matemticas de Intensidad Duracin Frecuencia para las principales ciudades de Bolivia entre otras ciudades para las que se estimo fueron tambin la ciudad de Santa Cruz, Trinidad, La Paz, El Alto, Sucre, Oruro, Potos, Cobija y Cochabamba.
Dnde: A, B y C son parmetros de ajuste de la funcin adoptada. t = duracin del evento en minutos. I= Intensidad en (mm/hrs). El ajuste para las recurrencias que van desde 2 y hasta 50 aos ha tomado como punto de referencia el trabajo realizado por el Ing. Csar Litwin (Consultora Berger, Proyecto de desages pluviales de las Avenidas Chaco y Mc Lean para la ciudad de Cochabamba, Programa Municipios, 2002), mientras que la funcin para TR = 100 aos ha sido calibrada con el Programa Solver de Excel (Microsoft Corp.), y los resultados de ese estudio se presentan en la tabla siguiente:
2 5 10 20 50 100
tabla de Coeficientes de ENDE Otra forma que se utilizara para validar los coeficientes calculados mediante los modelos Bernard, Sherman para la ciudad de Tarija ser por medio de bandas pluviogrficas. En nuestro caso existen las estaciones pluviogrficas de El Tejar y AASANA en el aeropuerto de la ciudad aunque actualmente ambas estn cerradas.
Luego de calcular las curvas Intensidad Duracin Frecuencia se procedi a aplicar metodologas propuestas para el presente proyecto mediante aproximaciones de mtodos numricos como el ajuste de curvas o aproximacin de curvas por regresin mltiple para lo cual se obtuvieron los siguientes resultados en funcin a las curvas IDF:
Ciudad en Estudio Tarija Yacuiba Bermejo Padcaya Iscayachi Villamontes Carapari Valle de la Concepcin Entre Ros San Lorenzo El Puente 0.162751 0.178261 0.159804 0.200696 0.130236 0.155675 0.167209 0.154788 0.150840 0.125369 0.186157 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 896.1416 1622.3812 1496.2761 838.9521 593.4699 1646.4446 1412.0461 746.4157 1242.9759 782.8759 495.3705
El modelo de Sherman parte del modelo de Bernard donde aparte de calcular los coeficientes , y calcula adems se estima un coeficiente denominado que trata de ajustar de mejor manera la familia de curvas. Por lo tanto luego de calibrar este coeficiente para cada ciudad se obtuvo los siguientes resultados: Ciudad en Estudio
Tarija Yacuiba Bermejo Padcaya Iscayachi Villamontes Carapari Valle de la Concepcin Entre Ros San Lorenzo El Puente 3 2 2 2 3 2 2 3 2 3 4 0.162751 0.173453 0.159804 0.200696 0.130236 0.155675 0.164920 0.158729 0.150840 0.125369 0.170281 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 0.80 896.1416 1713.1456 1496.2761 838.9521 593.4699 1646.4446 1450.0212 732.4949 1242.9759 782.8759 465.5254
Este modelo fue creado para duraciones inferiores a los 20 minutos y superior a los 5 minutos. El modelo parte previo la estimacin de las curvas de Intensidad Duracin Frecuencia. Es un modelo se lo calibra mediante mnimos cuadrados, donde el coeficiente a, b depende del periodo de retorno. Es as que se llegaron a los siguientes resultados:
Modelo de Talbot aplicado a la Ciudad de Tarija: Para lluvias de corta Duracin (5 min - 20 min) Periodo de Retorno 2 5 10 25 50 75 100 Coeficientes a 1936.61970 2397.19596 2665.60575 3206.18471 3554.59720 3758.40544 3903.00969 b 2.6941 2.6941 2.5619 2.6941 2.6941 2.6941 2.6941
Duracin de lluvia en (min) 5 251.70 311.56 352.51 416.71 461.99 488.48 507.28 7 199.77 247.28 278.77 330.74 366.68 387.70 402.62 9 165.61 204.99 230.55 274.17 303.97 321.39 333.76 10 152.56 188.84 212.20 252.57 280.02 296.08 307.47 11 141.42 175.05 196.55 234.13 259.57 274.45 285.01 12 131.80 163.14 183.05 218.20 241.91 255.78 265.62 13 123.40 152.75 171.29 204.29 226.49 239.48 248.69 14 116.01 143.60 160.95 192.06 212.93 225.13 233.80 15 109.45 135.48 151.78 181.20 200.89 212.41 220.58 17.5 95.90 118.71 132.87 158.77 176.02 186.11 193.28 20 85.34 105.63 118.15 141.28 156.63 165.61 171.98
2 5 10 25 50 75 100
Coeficientes a 1280.93674 1500.86261 1667.23011 1887.15598 2053.52348 2150.84223 2219.89098 b 2.6941 2.6941 2.6941 2.6941 2.6941 2.6941 2.6941
Modelo de Talbot aplicado a la Ciudad de Carapari: Para lluvias de corta Duracin (5 min - 20 min) Periodo de Retorno 2 5 10 25 50 75 100 Coeficientes a 3050.67370 3807.00599 4379.14928 5135.48158 5707.62487 6042.30724 6279.76816 b 2.6941 2.6941 2.6941 2.6941 2.6941 2.6941 2.6941
Modelo de Talbot aplicado a la Ciudad de El Puente: Para lluvias de corta Duracin (5 min - 20 min) Periodo de Retorno 2 5 10 25 50 75 100 Coeficientes a 1005.49246 1261.92994 1455.91740 1712.35488 1906.34234 2019.81774 2100.32981 b 2.6941 2.6941 2.6941 2.6941 2.6941 2.6941 2.6941
T 2 5 10 20 50 100
Duracin de lluvia en (Hrs) 0.5 66.02 76.64 85.79 96.03 111.48 119.08 124.79 0.75 47.73 55.41 62.02 69.43 80.60 86.09 90.22 1 37.92 44.02 49.27 55.16 64.03 68.39 71.67 1.5 27.41 31.82 35.62 39.88 46.29 49.45 51.82 5 10.46 12.15 13.60 15.22 17.67 18.87 19.78 8 7.18 8.34 9.34 10.45 12.13 12.96 13.58 12 5.19 6.03 6.75 7.56 8.77 9.37 9.82
2 5 10 20 50 75 100
Duracin de lluvia en (Hrs) 0.5 61.17 71.01 79.49 88.98 103.29 110.34 115.63 0.75 45.33 52.62 58.90 65.94 76.54 81.76 85.68 1 36.47 42.33 47.39 53.05 61.58 65.78 68.93 1.5 26.70 31.00 34.70 38.84 45.09 48.17 50.48 5 10.38 12.05 13.49 15.10 17.53 18.72 19.62 8 7.15 8.30 9.29 10.40 12.07 12.89 13.51 12 5.18 6.01 6.73 7.53 8.74 9.34 9.79
2 5 10 20 50 75 100
Podemos observar que los modelos de Bernard y Sherman se ajustan en gran manera a las interpretaciones de las bandas pluviogrficas. Siendo que el modelo de Sherman tiene mayor coeficiente de correlacin en comparacin al modelo de Bernard, pero ambos tienen un coeficiente de correlacin arriba de noventa y ocho por ciento por lo cual ambas ecuaciones son vlidas.
CURVAS IDF
Periodo de Retorno (aos)
Duracin de lluvia en (Hrs) 0.5 46.35 52.97 58.60 64.83 74.09 78.60 81.96 0.75 39.26 44.87 49.64 54.91 62.76 66.58 69.43 1 34.81 39.78 44.00 48.68 55.63 59.02 61.55 1.5 29.29 33.47 37.03 40.97 46.82 49.67 51.79 5 17.36 19.84 21.94 24.28 27.74 29.43 30.69 8 14.12 16.13 17.85 19.75 22.57 23.94 24.96 12 11.81 13.49 14.93 16.51 18.87 20.02 20.88
2 5 10 20 50 75 100
CURVAS IDF
Periodo de Retorno (aos)
Duracin de lluvia en (Hrs) 0.5 24.48 35.01 44.07 55.17 72.17 86.90 0.75 18.31 25.53 31.26 38.30 48.61 57.25 1 14.71 20.26 24.47 29.64 37.02 43.14 1.5 10.64 14.51 17.32 20.71 25.42 29.31 5 3.83 5.23 6.18 7.22 8.56 9.70 8 2.54 3.49 4.13 4.79 5.63 6.35 12 1.77 2.46 2.92 3.37 3.93 4.41
2 5 10 20 50 100
Podemos deducir en funcin a todas las grficas de correlacin mostradas anteriormente que la ecuacin IDF estimada tanto por Bernard como por Sherman arrojan buenos resultados, con respecto a las intensidades del pluviogrfo eso quiere decir que la variacin que hay no en gran magnitud, debido a su alta correlacin. Pero comparando con otras ecuaciones IDF (Alcantarillado Pluvial y ENDE) con la interpretacin de las bandas pluviogrficas podemos ver que la diferencia es un poco mas grande en funcin a su coeficiente de correlacin que tienen estas ecuaciones.
Nunca consideres al estudio como una obligacin, sino como una oportunidad que este maravilloso mundo del saber y que Dios nos dio