WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
PROYECTO:
OPTIMIZACIÓN DEL SISTEMAS DE ACUEDUCTO DEL
CORREGIMIENTO DE SAN ROQUE EN EL MUNICIPIO DE
CURUMANÍ DEPARTAMENTO DEL CESAR
INFORME HIDROLOGÍA
CURUMANÍ JUNIO DE 2021
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
TABLA DE CONTENIDO
INTRODUCCIÓN
1. ASPECTOS GENERALES
Hidrología
Localización General del proyecto
Ubicación especifica del proyecto
2. HIDROLOGÍA
Metodología
Estación pluviométrica seleccionada en la zona de estudio
Histograma de precipitación media mensual
Promedio de días con lluvia anual
Precipitación media anual
Precipitaciones máximas en 24 horas
Curvas de Intensidad – Frecuencia – Duración
Periodo de retorno
Hietograma
Caudales
Datos de cuenca
Intensidad
Tiempo de concentración
Coeficiente de escorrentía
Calculo de caudales
3. ANÁLISIS DE CUERPO DE AGUA FUENTE, MEDIANTE COMPARACIÓN
CON CUENCAS SIMILARES Y CERCANAS
Metodología
Estación limnimétrica
Periodo de retorno
Proyección de caudales
Datos históricos
Histograma
Parámetros estadísticos
Datos generales para cálculo de distribuciones estadísticas
Distribución normal
Distribución Logaritmo Pearson III
Distribución Gumbel
Caudal pico
Caudal mínimo
ANEXO 1: Datos crudos del IDEAM
CONCLUSIONES Y RECOMENDACIONES
BIBLIOGRAFÍA
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
INTRODUCCIÓN
La hidrología es un elemento a considerar en proyectos de abastecimiento de
agua, nos permite, a partir del conocimiento de datos históricos, hacer inferencias
sobre precipitaciones y escorrentías que deben considerarse para la protección de
las obras de captación y demás, que puedan ser afectadas por los cambios de
nivel y caudal de la fuente de abastecimiento y de cualquier cuerpo de agua que
exista en el corredor las obras del sistema.
Dado que la hidrología se basa en procesos estadísticos que permiten realizar
distribución de valores máximos o mínimos, se emplea para estimar las
precipitaciones que tienen probabilidad de ocurrencia en los periodos de retorno,
que a través de normas, se han establecido para el diseño de diversas estructuras
hidráulicas.
Entre mayor sea el periodo de retorno, o tiempo en que se considera o proyecta
que un evento pluvial o de caudal se dé o se repita, mayor será la magnitud del
mismo; por tanto, será de mayor tamaño las obras a ejecutar, por esta razón, y con
base en estudios que han comparado el costo de las posibles reparaciones por la
ocurrencia de dichos eventos y el costo de las obras para mitigarlos, desde la
normatividad se han definido periodos para cada tipo de obra.
Además de los métodos estadísticos, existen métodos simplificados que permiten
modelar las tormentas de diferentes periodos. Un ejemplo son las ecuaciones
regionalizadas de Rodrigo Vargas M. las cuales basan en información histórica y
en parámetros característicos de cada región en Colombia.
En el presente trabajo se empleará los métodos sintéticos para elaborar las curvas
de intensidad vs duración para cada frecuencia, esto para el análisis de la
escorrentía probable que potencialmente puede afectar la estructura de captación.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
1. ASPECTOS GENERALES
El municipio de Curumaní se encuentra ubicado en la zona de los Valles de los ríos
Cesar y Magdalena, donde la precipitación oscila entre 900 y 1500 mm, y su bajo
valor es debido a la acción secante de los vientos Alisios del noreste, que no
encuentran obstáculos orográficos en estos sectores. Presenta una temperatura
promedio anual de 28°C, con máxima de 39°C y mínima de 22°C, dependiendo del
régimen de lluvia anual. Esto permite clasificar la región dentro del piso térmico
cálido ardiente. El régimen de lluvias que se presenta en la zona, de acuerdo con
los datos pluviométricos es bimodal; es decir, existen dos (2) períodos lluviosos los
cuales registran una precipitación media anual de 1700 m; y dos (2) períodos
secos al año. Los períodos lluviosos ocurren regularmente entre los meses de abril
a junio y de agosto a noviembre; el período menos lluvioso tiene lugar en los
meses de diciembre a marzo, siendo el más seco enero. El comportamiento de los
Índices de aridez es deficitario, es decir, mayor a 0.3, lo cual significa que se
presenta déficit de agua en el suelo durante más de seis (6) meses al año.
Límites
Curumaní limita por el norte con el municipio de Chiriguaná, por el oriente con la
república de Venezuela, por el sur con el municipio de Chimichagua y con el
departamento de Norte de Santander y por el occidente con el municipio de
Chimichagua.
Corregimientos
Curumaní cuenta con siete (7) corregimientos:
Champán, El Mamey, Guaimaral, Sabanagrande, San Roque, San Sebastián y
Santa Isabel.
Economía
Su economía se basa en el sector primario, el cual se dedica a las actividades que
trabajan directamente los recursos naturales: agricultura, ganadería, explotación
forestal, pesca y el sector minero. También se incluyen las actividades de
ladrilleras y la actividad comercial.
Hidrografía
De acuerdo al estudio hidrográfico en el Corregimiento de San Roque encontramos
principalmente la Quebrada Simití, fuente abastecedora del acueducto
corregimental, el Río Anime, así como también caños, quebradas y arroyos, donde
sobresalen principalmente los caños: Cerro azul, La victoria, Las mujeres; las
quebradas: Aguas frías y Lamas verdes y los arroyos: Pitalito y El Valle.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Localización general del proyecto
UBICACIÓN GENERAL DEL PROYECTO
Fuente: Propia
Ubicación específica del proyecto
La estructura de la bocatoma se encuentra ubicada a la orilla de la quebrada
Simití, con el siguiente punto de referencia georreferenciado:
COORDENADAS GEOGRÁFICAS
Latitud Norte 9º 16’ 32.66”
Longitud Oeste 73º 27’ 47.12”
COORDENADAS PLANAS
Este 1067509
Norte 1517570
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
2. HIDROLOGÍA
Metodología
En las estructuras hidráulicas es necesario determinar la magnitud de los eventos
hidrológicos que pudieran ocurrir con cierta frecuencia, para valores máximos
durante un período determinado. Para esto, se deben procesar estadísticamente,
los datos observados en las estaciones de medición, con el fin de proyectar los
valores de eventos que podrían suceder en periodos definidos en la normatividad.
Dado que los caudales que interesan al estudio de un sistema como el que nos
ocupa, corresponde a fuentes de agua menores, las cuales no cuentan con datos
estadísticos de sus caudales, de tal manera que permitan inferencias a futuro, se
hace necesario el empleo de metodologías que permitan proyectar caudales que
puedan eventualmente concentrarse en los puntos donde se ubican estructuras
como la captación.
Para una bocatoma, el área o cuenca que concentra su escorrentía en el punto de
captación, se delimita por la divisoria de aguas que presenta la topografía de la
zona.
Para la estimación de los caudales en el caso descrito anteriormente y que es el
que nos ocupa, el manual de drenaje para carreteras del INVIAS incluye el método
simplificado, de las curvas de intensidad, duración, frecuencia, el cual se basa en
una fórmula empírica regionalizada para Colombia, y que tiene la siguiente
síntesis:
No obstante, para el presente trabajo se elaboraron las curvas de intensidad -
duración – frecuencia, a partir de la ecuación propuesta por Vargas y Díaz –
Granados-1, la cual tiene en cuenta en las proyecciones de precipitación, con
relación a la fórmula anteriormente mencionada, los datos del promedio multianual
de días con lluvia al año y el dato de precipitación anual promedio multianual.
1VARGAS M., Rodrigo y DIAZ GRANADOS, Mario. Universidad de los Andes. Curvas sintéticas regionalizadas de Intensidad –
Duración – Frecuencia para Colombia – 1998
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Tb
I = a x x Md x Ne x Ptf
tc
Para el cálculo de los caudales, se empleará el método racional y se hará una
proyección con base en método de transposición, con bese en caudales de
cuencas similares que cuentan con datos estadísticos.
Estación pluviométrica seleccionada en la zona de estudio
Para el análisis se emplean los datos de la estación pluviométrica ubicada en el
corregimiento de Rincón Hondo, Municipio de Chiriguaná, departamento del Cesar.
Esta estación se encuentra funcionando desde 1963 y a la fecha se encuentra
activa, por tanto sus datos son pertinentes para realizar la evaluación de la zona
del proyecto.
A continuación se resumen los datos de la estación seleccionada.
DATOS DE LA ESTACIÓN METEOROLÓGICA
Nombre RINCÓN HONDO
Código 25020260
Municipio CHIRIGUANÁ
Elevación msnm 100
Latitud 9º 23’ 49.3” N
Localización
Longitud 73º 29’ 16.9” W
Período de Registro para
1985 – 2020
estudio
Tecnología Convencional
Fuente: IDEAM
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Histograma precipitación media mensual multianual, estación pluviométrica
de Rincón Hondo
Fuente: Propia
Promedio de días con lluvia anual
Acorde a la información encontrada, el promedio de días con lluvia en la zona es:
AÑO DÍAS CON LLUVIA AÑO DÍAS CON LLUVIA
1985 67 2002 73
1986 68 2002 74
1987 94 2003 78
1988 118 2004 110
1989 106 2005 85
1990 81 2006 112
1991 57 2007 101
1992 82 2008 95
1993 97 2009 158
1994 90 2012 117
1995 105 2013 97
1996 118 2014 79
1997 60 2015 139
1998 93 2017 92
1999 109 2018 91
2000 92 2019 80
2001 65
PROMEDIO DÍAS CON LLUVIA MULTIANUAL 94.1
Fuente: IDEAM
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Precipitación media anual
Con base en los datos de la estación de Rincón Hondo, Ubicada en el municipio de
Chiriguaná, el valor promedio de la precipitación total anual es:
ENERO FEBRERO MARZO ABRIL MAYO JUNIO JULIO AGOSTO SEPTIEMBRE OCTUBRE NOVIEMBRE DICIEMBRE ANUAL
1985 10 0 137 395 313 30 223 168 133 271 155 158 1993
1986 0 41 15 93 394 77 60 228 193 402 22 22 1547
1987 15 30 150 257 223 97 246 247 297 599 156 146 2463
1988 0 25 60 326 367 487 249 597 333 266 344 20 3074
1989 56 86 398 113 300 134 140 362 706 261 283 183 3022
1990 0 0 0 351 233 146 212 263 220 546 226 94 2291
1991 0 0 236 49 149 117 147 80 405 216 99 0 1498
1992 27 0 12 148 290 217 183 183 413 187 389 0 2049
1993 0 33 189 341 305 160 174 255 318 111 397 132 2415
1994 58 120 144 295 240 164 205 396 179 501 170 68 2540
1995 29 10 211 104 125 382 474 635 378 396 256 182 3182
1996 90 62 305 174 595 377 215 217 476 356 185 31 3083
1997 61 0 121 81 220 139 91 115 313 267 253 0 1661
1998 0 107 113 238 437 214 207 279 339 357 382 306 2979
1999 18 211 292 198 275 405 141 249 799 488 558 160 3794
2000 72 77 52 192 161 205 267 117 487 262 321 185 2398
2001 0 0 144 82 316 99 237 81 211 494 152 217 2033
2002 0 37 92 74 304 207 139 42.5 341 177 226 122 1762
2003 0 22 134 100 172 361 133 150 251 360 104 1787
2004 21 126 39 320 275 214 123 110 355 319 273.1 60 2235
2005 138 68 108 143 331 248 112 90 209 484 590 113 2634
2006 102 0 62 210 235 183 93 305 233 400 486 105 2414
2007 0 8 115 378 385 201 302 381 312 644 389 26 3141
2008 0 82 286 280 477 285 294 385 378 350 378 0 3195
2009 161 23 299 277 229 468 244 295 205 288 280 58 2827
2010 0 51 156 273 303 275 515 590 561 641 858 574 4797
2011 75 124 323 365 565 278 264 278 568 1163 1113 339 5455
2012 0 0 52 363 262 344 126 436 440 700 403 0 3126
2013 0 0 200 340 525 193 248 340 509 551 309 85 3300
2014 49 75 123 155 372 37 32 239 383 503 723 176 2867
2015 55 106 225 156 182 125 169 184 128 217 251 79 1877
2016 37 71 144 323
2017 0 37 214 262 387 286 347 390 517 356 293 116 3205
2018 236 0 164 377 284 222 337 174 510.8 466 201 0 2972
2019 32 47 172 100 259 172 262 255 280 410 233 0 2222
2020 0 17 70 74 402 399 590 455 429 2436
PRECIPITACIÓN ANUAL PROMEDIO 2694
Fuente: IDEAM
Precipitaciones máximas en 24 horas
Tal como se mencionó, la estación seleccionada es la denominada Rincón Hondo,
ubicada en el municipio de Chiriguaná, operada por el IDEAM. El período de
registro empleado para el análisis es 1985-2020. Los datos de precipitación
máxima anual en 24 horas y el promedio multianual de la precipitación máxima
anual en 24 horas es:
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
PRECIPITACIÓN PRECIPITACIÓN
AÑO AÑO
MÁXIMA EN 24 HORAS MÁXIMA EN 24 HORAS
1985 110 2002 93
1986 52 2003 93
1987 170 2004 112
1988 105 2005 75
1989 115 2006 137
1990 63 2007 149
1991 130 2008 138
1992 147 2009 107
1993 78 2012 158
1994 84 2013 83
1995 98 2014 143
1996 106 2015 51
1997 74 2017 68
1998 98 2018 115
1999 152 2019 48
2000 97 2020 109
2001 97
PRECIPITACIÓN MÁXIMA EN 24 HORAS, PROMEDIO MULTIANUAL 105.1
Fuente: IDEAM
Datos en mm.
Curvas de Intensidad – Frecuencia – Duración
Para el cálculo de las curvas de intensidad – duración – frecuencia de la estación
Rincón Hondo, ubicada en el municipio de Chiriguaná, se utiliza la metodología
propuesta por Vargas y Díaz – Granados-2, la cual está en línea con la presentada
en el manual de drenaje para carreteras del INVIAS, solo que considera
parámetros adicionales en el cálculo.
De acuerdo con esta referencia, se pueden calcular curvas Intensidad –
Frecuencia – Duración, con la siguiente expresión:
b
I = a x T x Md x N e x Pt f
tc
2VARGAS M., Rodrigo y DIAZ GRANADOS, Mario. Universidad de los Andes. Curvas sintéticas regionalizadas de Intensidad –
Duración – Frecuencia para Colombia – 1998
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
En donde:
I = Intensidad de precipitación, mm/h
T= Período de retorno, años
t = Duración de lluvia, h
M =Promedio del valor máximo anual de precipitación en 24 h, mm
N = Promedio del número de días con lluvia al año
Pt= Precipitación media anual, mm
a,b,c,d,e y f = Parámetros dependientes de la ubicación geográfica de la estación
Para el área de influencia del proyecto, los valores correspondientes son:
a= 8,51
b= 0.21
c= 0.50
d= - 0.01
e= - 0.08
f= 0.28
De esta manera la ecuación queda:
T 0.21
I = 8.51 x x M-0.01 x N-0.08 x Pt0.28
T 0.50
A continuación se presentan los resultados para los diferentes periodos de retorno,
y para intervalos de duraciones de 10 minutos
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
RESULTADOS DE INTENSIDAD – FRECUENCIA – DURACIÓN
PERIODO DE RETORNO
t (MINUTOS) t (horas) Tr: 1 Tr: 5 Tr: 10 Tr: 20 Tr: 50 Tr: 100
INTENSIDAD
10 0.17 127.45 178.70 206.70 239.09 289.82 335.23
20 0.33 91.48 128.26 148.36 171.60 208.01 240.61
30 0.50 74.32 104.20 120.53 139.41 168.99 195.47
40 0.67 64.20 90.01 104.12 120.43 145.99 168.86
50 0.83 57.68 80.87 93.55 108.20 131.16 151.71
60 1.00 52.55 73.68 85.22 98.58 119.49 138.22
70 1.17 48.58 68.12 78.79 91.14 110.47 127.78
80 1.33 45.57 63.89 73.90 85.48 103.62 119.85
90 1.50 42.91 60.16 69.59 80.49 97.57 112.85
100 1.67 40.66 57.02 65.95 76.28 92.47 106.96
110 1.83 38.85 54.47 63.00 72.87 88.33 102.17
120 2.00 37.16 52.10 60.26 69.71 84.50 97.73
130 2.17 35.67 50.02 57.85 66.92 81.12 93.83
140 2.33 34.43 48.27 55.83 64.58 78.28 90.55
150 2.50 33.23 46.60 53.90 62.35 75.58 87.42
160 2.67 32.16 45.09 52.16 60.33 73.13 84.59
170 2.83 31.24 43.80 50.66 58.60 71.03 82.16
180 3.00 30.34 42.54 49.20 56.91 68.99 79.80
190 3.17 29.51 41.38 47.87 55.37 67.11 77.63
200 3.33 28.80 40.38 46.70 54.02 65.48 75.74
210 3.50 28.09 39.38 45.55 52.69 63.87 73.88
220 3.67 27.43 38.46 44.49 51.46 62.38 72.15
230 3.83 26.85 37.65 43.55 50.37 61.06 70.63
240 4.00 26.27 36.84 42.61 49.29 59.75 69.11
250 4.17 25.73 36.08 41.73 48.27 58.52 67.69
260 4.33 25.25 35.41 40.96 47.37 57.43 66.42
270 4.50 24.77 34.73 40.18 46.47 56.33 65.16
280 4.67 24.32 34.09 39.44 45.62 55.30 63.96
290 4.83 23.91 33.53 38.78 44.85 54.37 62.89
300 5.00 23.50 32.95 38.11 44.09 53.44 61.81
310 5.17 23.11 32.40 37.48 43.35 52.55 60.79
320 5.33 22.76 31.91 36.91 42.70 51.76 59.87
330 5.50 22.41 31.42 36.34 42.03 50.95 58.94
340 5.67 22.07 30.94 35.79 41.40 50.18 58.05
350 5.83 21.76 30.52 35.30 40.83 49.49 57.24
360 6.00 21.45 30.08 34.79 40.24 48.78 56.43
Fuente: Propia
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
CURVAS DE INTENSIDAD – FRECUENCIA – DURACIÓN
Fuente: Propia
Periodo de retorno
Acorde al numeral 3 del artículo 53 del RAS 2017, se debe analizar eventos con un
periodo de retorno de 100 años, con el fin de proyectar los muros protectores de la
estructura de captación.
Hietograma
Los hietogramas son herramientas que permiten alimentar el programa Epa Swmm
para un análisis alterno de la cuenca que alimenta la fuente de la cual se capta el
agua para el acueducto de San Roque. Los hietogramas se elaboran acorde con
los periodos de retorno a aplicar en la parte de cálculos hidráulicos.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
HIETOGRAMA 100 AÑOS
Fuente: Propia
Caudales
Para el calculo de caudales, se emplea el metodo racional, el cual se basa en la
formula
Q = 0.00278 x A x i x C
Donde:
Q: Caudal en m3/s
A: Area de drenaje. Corresponde al área de la subcuenca en Ha
i: intensidad de la lluvia en mm/hr
C: Coeficiente de escorrentía
Áreas de drenaje
El área se drenaje hace referencia a la cuencan que alimenta la fuente de
abastecimiento, hasta el punto de captación, la cual se establece con la ayuda de
la imagen satelital de Google earth y con la ayuda de información topográfica del
proyecto.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
CUENCAS RIO SIMITI HASTA LA CAPTACIÓN
Fuente: Google Earth
Datos cuenca
Área: 21.5 Km2
Pendiente media: 7.6
Longitud: 8.74 Km
Ancho Medio: 2.92 Km
Intensidad
La Intensindad se obtiene a partir de las curvas IDF, entrando con el tiempo de
concentración.
Tiempo de concentración
Es el tiempo que se considera necesario para que toda la cuenca contribuya al sitio
de concentración en análisis, en este caso corresponde al sitio de ubicación de la
alcantarilla. Este tiempo se cuenta a partir del inicio de la precipitación.
Para el cálculo del tiempo de concentración existen diversas fórmulas, las cuales
se presentan en el manual de drenaje para carreteras del IVIAS. Para el presente
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
trabajo se utilziará la ecuación de la onda cinemática de Morgali y Linsley, y Aron y
Erborge.
La siguente en la síntesis de la fómula:
Donde:
Tc: tiempo de concentración en minutos
L: Longitud del cauce principal en m.
n: Coeficiente de rugosidad de manning de la superficie.
S: Pendiete media del cauce principal de drenaje de la cuenca en m/m
Para el calculo de los cauales probables en el punto de bocatoma del sistema de
acueducto del corregimiento de San Roque, se obtiene:
En la ecuación de la onda cinética, los perámetros de tiempo de concentración y la
intensisdad pluvial, son interdepoendientes, razon por la que se requiere realizar
iteraciones hasta encontrar el valor adecuado, para este caso: el tiempo de
concentración es:
Para T = 100 años
6.998 x 8740^0.6 x 0.04^0.6
Tc = = 73.9 min
124.5^0.4 x 0.076^0.3
Para T = 10 años
6.998 x 8740^0.6 x 0.04^0.6
Tc = = 94.1 min
68.05^0.4 x 0.076^0.3
Para T = 1 año
6.998 x 8740^0.6 x 0.04^0.6
Tc = = 119.9 min
37.17^0.4 x 0.076^0.3
Corroborado con los datos de las curvas IDF, para un valor de Tc de 73.9, la
intensidad es de 124.5 mm/hr.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Coeficiente de escorrentía
El corefociente de escorrentía se define a partir de la tabla 2.9 del manual de
drenaje para carreteras del INVIAS, en la cual se presenta los coeficientes de
escorrentía para areas rurales.
Fuente: Manual de drenaje para carreteras - INVIAS
Para la modelación, en cada subcuenca se define el coeficiente de escorrentía,
acorde con las caracteristicas de urbanismo de la misma.
Cálculo de caudales
Lo caudales se definen con el metodo racional, el cual tiene en cuenta el área de
drenaje, el coeficiente de escorrentía, el tiempo de concentración, el periodo de
retorno y la consecuente intensidad pluviométrica.
Aplicando la formula precedente, los caudales para el periodo de revisión son:
Q = 0.00278 x A x i x C
PERIODO DE RETORNO CAUDAL
100 años 128
10 año 70
1 año 38
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
3. ANÁLISIS DE CUERPO DE AGUA FUENTE, MEDIANTE COMPARACIÓN
CON CUENCAS SIMILARES Y CERCANAS
Tendiendo en cuenta que el metodo racional no no es muy recomendado para el
análisis de áreas mayores, se procede a la comparacion de la cuenta en
estudio,con una cuenca aledaá y similar en em du topografiaa y forma. Lo anterior
aprovechando la existencia de datos estadisticos de los caudles del cuerpo de
agua que la drena.
Metodología
En las estructuras hidráulicas es necesario determinar la magnitud de los eventos
hidrológicos que pudieran ocurrir con cierta frecuencia, para valores máximos
durante un período determinado. Para esto, se deben procesar estadísticamente,
los datos observados en las estaciones de medición, con el fin de proyectar los
valores de eventos que podrían suceder en periodos definidos en la normatividad.
Con el fin de trascender los datos históricos disponibles en las estaciones de
medición, es necesario ajustar los datos a distribuciones de probabilidad
conocidas, como la Gumbel o Pearson tipo III.
Para el caso de estudio, los caudales históricos se obtienen de la estación más
cercana encontrada, la cual está vigente y tiene un rango de datos adecuado.
Estación limnimétrica seleccionada en la zona de estudio
Para el análisis se emplean los datos de la estación limnimétrica ubicada en el
municipio de Chiriguaná, departamento del Cesar, la cual se encuentra
funcionando desde 1963 y a la fecha se encuentra activa, por tanto, cuenta con
datos suficientes y pertinentes para realizar una proyección.
Datos de la estación seleccionada.
DATOS DE LA ESTACIÓN
Nombre PUENTE CARRETERA
Código 25027890
Municipio CHIRIGUANÁ
Cuenca LA MULA
Elevación msnm 120
Latitud 9º 20’ 24” N
Localización
Longitud 73º 29’ 24” W
Período de Registro usado 1980 – 2018
Fuente: IDEAM
Periodo de retorno
Se asume un periodo de retorno de 100 años.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Proyección de caudales
Para la proyección de los caudales en el punto de la estación limnigráfica
denominada el Reposo, se emplean la distribución de probabilidad de Gumbel, log-
Pearson tipo III y extremos.
Datos históricos
Los datos históricos se toman de
DATOS HISTÓRICOS ESTACIÓN PUENTE CARRETERA CHIRIGUANÁ
AÑO CAUDAL M3 AÑO CAUDAL M3
1980 21.56 1995 17.03
1981 21.78 1996 18.83
1982 16.44 1997 15.74
1984 53.43 1998 24.95
1985 11.78 1999 10.63
1986 13.25 2005 11.22
1987 23.17 2006 11.02
1988 18.90 2007 12.90
1989 14.01 2008 11.21
1990 21.34 2010 11.32
1991 15.84 2011 17.85
1993 15.79 2017 17.30
1994 17.02 2018 27.02
Fuente: Propia
Histograma
Fuente: Propia
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Parámetros estadísticos
Valor medio
El valor medio de la serie histórica corresponde al promedio de la misma
Xm = Xi / N => i de 1 a N
Donde:
Xi: Datos individuales de la serie histórica
N: Número de datos históricos
Desviación estándar
Para el cálculo de la desviación estándar se emplea la siguiente fórmula.
Sm = { 1 / N x ( Xi - Xm ) 2 } 0.5 => i de 1 a N
Donde:
Xi: Datos individuales de la serie histórica
Xm: Valor medio de la serie
N: Número de datos históricos
Coeficiente de variación
Para el cálculo del coeficiente de variación se emplea:
Cv = Sm / Xm
Cs = a / Sm 3
a = { N / (( N – 1 ) x ( N – 2 )) } k
k = ( Xi - Xm ) 3 => i de 1 a N
Donde:
Xi: Datos individuales de la serie histórica
Xm: Valor medio de la serie
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
N: Número de datos históricos
Datos generales para cálculo de distribuciones estadísticas
DATOS DE CAUDALES MÁXIMOS RIO LA MULA
AÑO CAUDAL m POSICIÓN Q ORDEN (Xi) Pem Yi (log(Xi) (Yi - Ym)^2 (Yi - Ym)^3 m / (n+1) 1- m / (n+1) (Xi - Xm)^3
1980 21.56 1 6 53.43 3.70 1.73 0.25 0.13 27.00 0.96 43994.17
1981 21.78 2 5 27.02 7.41 1.43 0.04 0.01 13.50 0.93 703.05
1982 16.44 3 14 24.95 11.11 1.40 0.03 0.00 9.00 0.89 317.48
1984 53.43 4 1 23.17 14.81 1.36 0.02 0.00 6.75 0.85 128.17
1985 11.78 5 21 21.78 18.52 1.34 0.01 0.00 5.40 0.81 48.70
1986 13.25 6 19 21.56 22.22 1.33 0.01 0.00 4.50 0.78 40.42
1987 23.17 7 4 21.34 25.93 1.33 0.01 0.00 3.86 0.74 33.14
1988 18.90 8 8 18.90 29.63 1.28 0.00 0.00 3.38 0.70 0.46
1989 14.01 9 18 18.83 33.33 1.27 0.00 0.00 3.00 0.67 0.35
1990 21.34 10 7 17.85 37.04 1.25 0.00 0.00 2.70 0.63 -0.02
1991 15.84 11 15 17.30 40.74 1.24 0.00 0.00 2.45 0.59 -0.57
1993 15.79 12 16 17.03 44.44 1.23 0.00 0.00 2.25 0.56 -1.32
1994 17.02 13 13 17.02 48.15 1.23 0.00 0.00 2.08 0.52 -1.36
1995 17.03 14 12 16.44 51.85 1.22 0.00 0.00 1.93 0.48 -4.81
1996 18.83 15 9 15.84 55.56 1.20 0.00 0.00 1.80 0.44 -11.98
1997 15.74 16 17 15.79 59.26 1.20 0.00 0.00 1.69 0.41 -12.78
1998 24.95 17 3 15.74 62.96 1.20 0.00 0.00 1.59 0.37 -13.62
1999 10.63 18 26 14.01 66.67 1.15 0.01 0.00 1.50 0.33 -69.84
2005 11.22 19 23 13.25 70.37 1.12 0.01 0.00 1.42 0.30 -116.08
2006 11.02 20 25 12.90 74.07 1.11 0.01 0.00 1.35 0.26 -142.90
2007 12.90 21 20 11.78 77.78 1.07 0.02 0.00 1.29 0.22 -255.82
2008 11.21 22 24 11.32 81.48 1.05 0.03 -0.01 1.23 0.19 -315.55
2010 11.32 23 22 11.22 85.19 1.05 0.03 -0.01 1.17 0.15 -329.66
2011 17.85 24 10 11.21 88.89 1.05 0.03 -0.01 1.13 0.11 -331.10
2017 17.30 25 11 11.02 92.59 1.04 0.03 -0.01 1.08 0.07 -359.13
2018 27.02 26 2 10.6 96.30 1.03 0.04 -0.01 1.04 0.04 -421.55
0.6053 0.1073 42877.8
Fuente: Propia
Distribución normal
N: 26 Ym: 1.2014 valor medio logarítmico
Xm: 17 Sy: 0.1212 Desviación estándar logarítmica
Sm: 5 Cs: 0.1387 Coeficiente asimetría logarítmico
Cv: 0.28 k: 0.0231 (Cs / 6 )
a: 58
Cs: 0.57
Para la distribución normal se emplean las siguientes fórmulas:
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
P: Probabilidad según periodo de retorno
Xt = Caudal para el periodo de retorno
Xt = Xm + Z x Sm
Donde:
Xt: Caudal en m3/s
Z = Variable aleatoria
Xm: Valor medio serie de datos
Comportamiento de los datos históricos con respecto a la distribución
Distribución Logaritmo Pearson III
Para la distribución Logaritmo Pearson III se emplean las siguientes fórmulas:
Yt = Ym + KT x Sy
Donde:
Yt: Logaritmo del caudal
KT = Variable aleatoria
Ym: Valor medio serie de datos
Sy: Desviación estándar logarítmica
Xt = Antilogaritmo de Yt, en m3/s (caudal)
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Distribución Gumbel
Para la distribución Gumbel se emplean las siguientes fórmulas:
Donde:
KT = Variable aleatoria
T = Periodo de retorno
Xt = Xm + KT x Sm
Donde:
Xt = Caudal proyectado para el periodo de retorno en m3/s
Xm = Valor medio de la serie histórica
KT = Variable aleatoria
Sm = Desviación estándar serie histórica
Xt = Antilogaritmo de Yt, en m3/s (caudal)
COMPORTAMIENTO DE LAS PROYECCIONES REALIZADAS
PERIODO
w z Xt Kt Yt Xt KT Xt Yt Xt
RETORNO (T) p
2 0.500 1.18 -0.0000001 18 -0.1969 1.19666 16 -0.1643 17 1.2017 16
3 0.333 1.48 0.4303000 22 0.2354 1.26393 18 0.2538 20 1.2668 18
5 0.200 1.79 0.8415000 25 0.7228 1.33977 22 0.7195 24 1.3393 22
10 0.100 2.15 1.2817000 29 1.3320 1.43456 27 1.3046 29 1.4303 27
20 0.050 2.45 1.6452000 32 1.9081 1.52420 33 1.8658 34 1.5176 33
25 0.040 2.54 1.7511000 33 2.0890 1.55235 36 2.0438 36 1.5453 35
50 0.020 2.80 2.0542000 36 2.6407 1.63819 43 2.5923 40 1.6307 43
100 0.010 3.03 2.3268000 38 3.1814 1.72233 53 3.1367 45 1.7154 52
200 0.005 3.26 2.5762000 40 3.7145 1.80528 64 3.6791 49 1.7998 63
500 0.002 3.53 2.8785000 43 4.4118 1.91378 82 4.3947 56 1.9111 81
1000 0.001 3.72 3.0905000 44 4.9355 1.99526 99 4.9355 60 1.9953 99
Fuente: Propia
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Fuente: Propia
Como se puede ver en la gráfica anterior, la distribución a la que más se acomoda
la serie histórica, es la distribución de Gumbel.
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
Caudal pico
El caudal para el periodo de retorno de 100 años, es de: 45 m3/s.
El caudal para el periodo de retorno de 10 años, es de: 29 m3/s.
El caudal para el periodo de retorno de 1 años, es de: 11 m3/s.
Para el caso de estudio, los caudales históricos se obtienen de la estación más
cercana encontrada, la cual está vigente y tiene un rango de datos adecuado.
Considerando que la estación limnimétrica no está en el punto de obra, para
establecer los caudales del río Cesar más cercanos a esta ubicación, se empleará
el método de transposición de datos, el cual se usa cuando la estación existente,
no se encuentra en el sitio exacto de la obra, pero si pertenece a la misma cuenca.
El cálculo de los caudales se hace empleando la siguiente fórmula:
Qp = Qe x ( Ap / Ae )x
Donde:
Qp: Caudal en el punto del proyecto
Qe: Caudal en la estación limnimétrica
Ap: Área en el punto del proyecto
Ae: Área en la estación limnimétrica proyecto
X: coeficiente estandarizado a falta de datos de investigación
Área de la cuenca de la fuente de abastecimiento del sistema de acueducto de San
Roque.
Ap: 2150 Ha
Área de la cuenca del río La Mula hasta la estación limnimétrica Puente Carretera.
Ae: 34100 Ha
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
CUENCA RIO LA MULA
Fuente Base: Google earth – Fuente análisis: Propia
Entonces
Para T = 100 años Qp = 11.5 m3/s
Para T = 10 años Qp = 7.4 m3/s
Para T = 1 años Qp = 2.8 m3/s
Caudal mínimo
Considerando los datos estadísticos existentes en los últimos 25 años del río La
Mula, es decir desde 1995, con relación a los caudales mínimos anuales
reportados, es de 140 l/s. por tanto y empleando la misma metodología de la
transposición, el caudal mínimo con probabilidad de ocurrir el periodo de retorno,
es de:
140 X ( 22 / 341 ) ^ 0.5 = 35.6 l/s
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
CONCLUSIONES Y RECOMENDACIONES
El analisis de la protección de la estrucutra de captación, serian los obtenidos a
traves del método de transposición para periodos de retorno 100, 10 y 1 año.
Aunque las proyecciones hidrológicas se hacen con base en datos estadisticos,
lo cual implica que no son valores necesariamente ciertos, es imporatnte que se
con abse en los datos resultantes, se analice el sistema, buscando prevver
situaciones para alores extremos.
Los datos de caudales obtenidos mediante el método racional, no se
recomiendan para el análisis de las estruturas, dado que por magnitud de la
cuenca, los valores obtenidos son altos con relación a cuencas similares, tal es
elc aso del río La Mula.
__________________________________
WILSON FERNÁNDEZ RODRÍGUEZ
Ingeniero Civil
WILSON FERNÁNDEZ RODRÍGUEZ – INGENIERO CIVIL UAN
BIBLIOGRAFÍA
Reglamento Técnico para el Sector de Agua Potable y Saneamiento Básico
RAS -2017
Manual de drenaje para carreteras – INVIAS
Curvas sintéticas regionalizadas de Intensidad – Duración – Frecuencia para
Colombia - Rodrigo Vargas M.
Hidrología básica y aplicada – Allen Bateman
Fundamentos de hidrología de superficie – Francisco Javier Aparicio Mijares
