i

UNIVERSIDAD FRANCISCO DE PAULA SANTANDER OCAÑA

Documento Código Fecha Revisión
FORMATO HOJA DE RESUMEN PARA TRABAJO DE F-AC-DBL-007 10-04-2012 A
GRADO
Dependencia Aprobado Pág.
DIVISIÓN DE BIBLIOTECA SUBDIRECTOR ACADEMICO i(161)

RESUMEN – TRABAJO DE GRADO

STPHANNIE CAROLINA CAVIEDES NUMA

AUTORES
JULIAN PÉREZ RINCÓN
FACULTAD FACULTAD DE CIENCIAS AGRARIAS Y DEL AMBIENTE
PLAN DE ESTUDIOS PLAN DE ESTUDIOS DE INGENIERÍA AMBIENTAL
DIRECTOR WILSON ANGARITA CASTILLA
CAPACIDAD MÁXIMA DE ASIMILACIÓN DEL RÍO
ALGODONAL POR LA DESCARGA DE VERTIMIENTOS DE
TÍTULO DE LA TESIS
AGUAS RESIDUALES ENTRE LOS MUNICIPIOS DE
ÁBREGO Y OCAÑA, NORTE DE SANTANDER
RESUMEN
(70 palabras aproximadamente)

EN ESTE TRABAJO SE PRESENTA LOS RESULTADOS DE LA CAPACIDAD DE

AUTODEPURACIÓN DEL RÍO ALGODONAL FRENTE A LA DESCARGA DE

VERTIMIENTOS DE AGUAS RESIDUALES DOMÉSTICAS, A TRAVÉS DEL MODELO

MATEMÁTICO DE CALIDAD DEL AGUA QUAL2K, CON EL FIN DE EVALUAR EL

COMPORTAMIENTO DE LOS PARAMETROS DE CALIDAD Y QUE SIRVA COMO UNA

HERRAMIENTA DE APOYO A LA PLANIFICACIÓN DE LA GESTIÓN INTEGRAL DEL

RECURSO HÍDRICO.

CARACTERÍSTICAS
PÁGINAS: 162 PLANOS: ILUSTRACIONES: 61 CD-ROM: 1
 ii

CAPACIDAD MÁXIMA DE ASIMILACIÓN DEL RÍO ALGODONAL POR LA DESCARGA

DE VERTIMIENTOS DE AGUAS RESIDUALES ENTRE LOS MUNICIPIOS DE ÁBREGO

Y OCAÑA, NORTE DE SANTANDER

AUTORES:

STPHANNIE CAROLINA CAVIEDES NUMA

JULIAN PÉREZ RINCÓN

Proyecto de investigación para optar el título de ingeniero ambiental

Director:

WILSON ANGARITA CASTILLA

UNIVERSIDAD FRANCISCO DE PAULA SANTANDER

FACULTAD DE CIENCIAS AGRARIAS Y DEL AMBIENTE

INGENIERÍA AMBIENTAL

Ocaña, Colombia Julio, 2018

iii
 iv

Agradecimientos

A nuestras familias por su apoyo incondicional, por ser el motor de lucha y la voz de aliento

ante las adversidades. Gracias por creer en nosotros.

A nuestro director Wilson Angarita Castilla por orientarnos y apoyarnos en el proceso de

investigación y ejecución del proyecto, por impulsar y reforzar nuestros conocimientos y así

cumplir con cada uno de nuestros objetivos.

Al grupo de investigación MINDALA por sembrar en nosotros el espirito de investigación y por

acercarnos a la realidad ambiental que se vive en nuestra región.

A nuestros compañeros y amigos Heiner C. Vergel Cañizares, Diego A. Díaz Jiménez, Arley

Arenas Sánchez, María F. Mora Sanjuan, Héctor Cabrera Rodríguez, Camila F. Sánchez

Hamilton, y Tania L. Jaime Blanco por brindarnos su apoyo y por ser partícipes en la

realización del proyecto.

 v

Tabla de contenido

Capítulo 1 Capacidad máxima de asimilación del río Algodonal por la descarga de

vertimientos de aguas residuales entre los municipios de Ábrego y Ocaña, Norte de

Santander .............................................................................................................................. 1

1.1. Problema de investigación ................................................................................................... 1

1.2. Planteamiento del problema .............................................................................................. 1

1.3. Formulación del problema ................................................................................................... 4

1.4. Objetivos de investigación ................................................................................................... 5

1.4.1. Objetivo general. ......................................................................................................... 5

1.4.2. Objetivos específicos. ................................................................................................... 5

1.5. Justificación ........................................................................................................................ 6

1.6. Formulación de la Hipótesis ................................................................................................ 7

1.6.1. Problema. ...................................................................................................................... 7

1.6.2. Pregunta de investigación. ............................................................................................ 8

1.7. Delimitaciones .................................................................................................................... 8

1.7.1. Delimitación operativa. .............................................................................................. 8

1.7.2. Delimitación conceptual. ............................................................................................ 9

1.7.3. Delimitación geográfica. ............................................................................................ 9

1.7.4. Delimitación temporal. ............................................................................................... 13

Capítulo 2 Marco referencial .............................................................................................. 14

 vi

2.1. Marco histórico ................................................................................................................ 14

2.2. Marco conceptual ............................................................................................................. 17

2.2.1. El agua. ....................................................................................................................... 17

2.2.2. Manejo del recurso hídrico. ........................................................................................ 19

2.2.3. El ciclo hidrológico. .................................................................................................... 19

2.2.4. Autodepuración. .......................................................................................................... 20

2.3. Marco teórico ................................................................................................................... 24

2.3.1. Autodepuración de cauces naturales ........................................................................... 24

2.3.2. Tipos de modelos. ....................................................................................................... 26

2.3.3. Modelos de calidad de agua. ....................................................................................... 27

2.3.4. Dimensiones espaciales de los modelos de sistemas. ................................................. 28

2.3.5. Parámetros empleados en la modelación. ................................................................... 29

2.3.6. Constantes universales. ............................................................................................... 30

2.3.7. Construcción de modelos dinámicos. ......................................................................... 31

2.4. Marco legal ....................................................................................................................... 32

Capítulo 3 Metodología de la investigación ........................................................................ 35

3.1. Diseño de la metodología de investigación .................................................................... 35

3.2. Población y muestra ......................................................................................................... 36

3.2.1. Población. .................................................................................................................. 36

3.2.2. Muestra. ..................................................................................................................... 36

 vii

3.3. Diagrama de la metodología ........................................................................................... 36

3.4. Técnicas de recolección de la información .................................................................... 37

3.5. Recolección de la información climatológica .................................................................... 41

3.6 Morfometría de la cuenca ................................................................................................... 41

3.7. Selección del modelo de simulación .................................................................................. 42

3.8. Metodología de análisis de calidad del agua...................................................................... 43

3.8.1. Ubicación de los vertimientos ..................................................................................... 43

3.8.2. Cargas contaminantes ................................................................................................. 43

3.8.3. Selección de las estaciones de monitoreo ................................................................... 45

3.8.4. Tiempos de viaje ......................................................................................................... 45

3.8.5. Aforo ........................................................................................................................... 46

3.8.6. Toma de muestras ....................................................................................................... 47

3.9. Metodología de la modelación ........................................................................................... 48

3.10. Instrumentos utilizados para la recolección de la información. ....................................... 48

3.11. Procesamiento y análisis de la información ................................................................. 50

3.11.1. Procesamiento de la información climatológica ................................................... 50

3.11.2. Morfometría del río ................................................................................................... 50

3.11.3. Verificación del modelo de simulación. ............................................................... 50

3.11.4. Estaciones de monitoreo......................................................................................... 50

3.11.5. Calibración del modelo de calidad del agua. ........................................................ 51

 viii

3.11.6. Características hidráulicas ........................................................................................ 54

3.11.7 Escenarios prospectivos de calidad del cauce principal del río Algodonal. .............. 55

Capítulo 4 Resultados y discusión ...................................................................................... 55

4.1. Características naturales del río Algodonal. ...................................................................... 55

4.1.1. Climatología. ............................................................................................................... 55

4.1.2. Características morfométricas de la cuenca ................................................................ 60

4.2. Selección del modelo de simulación .................................................................................. 73

4.3. Análisis de la calidad del agua ........................................................................................... 79

4.3.1. Vertimientos significativos que afectan la corriente principal ................................... 79

4.3.2. Captaciones ................................................................................................................. 80

4.3.3. Cargas contaminantes ................................................................................................. 80

4.3.4. Características de las estaciones de monitoreo ........................................................... 84

4.3.5. Tiempos de viaje ......................................................................................................... 90

4.3.6. Aforo ........................................................................................................................... 92

4.3.7. Resultados de los análisis de calidad del agua ............................................................ 92

4.3.8. Perfiles de calidad ....................................................................................................... 94

4.4. Modelación de la capacidad de autodepuración .............................................................. 104

4.4.1. Topología del río Algodonal. .................................................................................... 104

4.4.2. Hojas de entrada al modelo. ...................................................................................... 104

4.4.3. Escenario actual de calidad del cauce principal del río Algodonal. ......................... 110
 ix

4.4.4. Calidad del agua en el río Algodonal frente a los escenarios prospectivos. ............. 116

Capítulo 5 Conclusiones ................................................................................................... 118

Capítulo 6 Recomendaciones ............................................................................................ 121

Capítulo 7 Referencias ...................................................................................................... 122

Capítulo 8 Apéndice ......................................................................................................... 131

Apéndice A. Protocolo de monitoreo del río algodonal ......................................................... 131

Apéndice B. Etapas de la campaña de aforo y muestreo. ....................................................... 134

Apéndice C. Variables fisicoquímicas y microbiológicas y sus métodos de evaluación ....... 136

Apéndice D. Plan de muestreo para el análisis físico-químico y microbiológico. ................. 136

Apéndice E. Recursos humanos. ............................................................................................. 137

Apéndice F. Recursos institucionales. .................................................................................... 137

Apéndice G. Recursos financieros. ......................................................................................... 138

Apéndice H. Trabajo en campo. ............................................................................................. 139

Apéndice I. Ficha técnica y etiqueta. ...................................................................................... 140

Apéndice J. Laboratorio. ......................................................................................................... 142

 x

Lista de tablas

Tabla 1 Coordenadas de la cuenca del río Algodonal ................................................................. 10

Tabla 2 Distribución porcentual del río Algodonal por municipio .............................................. 11

Tabla 3 Tipología y clasificación de modelos............................................................................... 26

Tabla 4 Estaciones climatológicas del IDEAM. ........................................................................... 41

Tabla 5 Información de los vertimientos. ..................................................................................... 44

Tabla 6 Histórico de tiempo de los municipios de Ábrego y Ocaña. ............................................ 47

Tabla 7 Descripción de equipos y herramientas utilizadas. ......................................................... 48

Tabla 8 Coeficiente de rugosidad Manning para canales de corriente natural (Chow et al. 1988).

....................................................................................................................................................... 54

Tabla 9 Valores típicos de las aguas residuales domésticas. ....................................................... 55

Tabla 10 Parámetros morfométricos de la cuenca del río Algodonal. ......................................... 61

Tabla 11 Pendiente media de la cuenca. ...................................................................................... 66

Tabla 12 Parámetros relacionados con el relieve de la cuenca. .................................................. 66

Tabla 13 Áreas entre curvas de nivel. ........................................................................................... 68

Tabla 14 Datos de entrada al modelo. .......................................................................................... 75

Tabla 15 Características climatológicas de las estaciones en las horas del muestreo en la

campaña 1. .................................................................................................................................... 76

Tabla 16 Características climatológicas de las estaciones en las horas del muestreo en la

campaña 2. .................................................................................................................................... 77

Tabla 17 Concesiones del río Algodonal ...................................................................................... 80

Tabla 18 Cálculo del caudal vertido en m3/s en el municipio de Ábrego. .................................... 82

Tabla 19 Cálculo del caudal vertido en m3/s en el corregimiento de La Ermita.......................... 83

 xi

Tabla 20 Cálculo de carga contaminante del vertimiento de la laguna de oxidación en kg/día. 83

Tabla 21 Tabla E.2.6 (RAS 2000) Aportes per cápita para aguas residuales domésticas. ......... 84

Tabla 22 Información geográfica de las estaciones de muestreo. ................................................ 84

Tabla 23 Caracterización hidráulica de las estaciones. .............................................................. 86

Tabla 24 Características organolepticas de las estaciones de monitoreo. .................................. 87

Tabla 25 Obtención de los tiempos de viaje en tiempo seco......................................................... 91

Tabla 26 Obtención de los tiempos de viaje en tiempo de lluvia. ................................................. 91

Tabla 27 Datos hidráulicos en las estaciones de monitoreo. ...................................................... 92

Tabla 28 Hora de la toma de muestras en tiempo seco. ............................................................... 92

Tabla 29 Datos de calidad del agua en las estaciones de monitoreo en tiempo seco. ................ 93

Tabla 30 Hora de la toma de muestras en tiempo de lluvia. ........................................................ 93

Tabla 31 Datos de calidad del agua en las estaciones de monitoreo en tiempo de lluvia. .......... 94

Tabla 32 Efectos de la temperatura en la vida acuática. ............................................................. 96

 xii

Lista de figuras

Figura 1. Mapa de localización del río Algodonal. ...................................................................... 10

Figura 2. La confluencia del río Oroque y río Frío da origen al río Algodonal. .......................... 12

Figura 3. Diagrama de la metodología. ........................................................................................ 37

Figura 4. Esquema de segmentación del Q2k para un río con afluentes. ..................................... 52

Figura 5. División del tramo en segmentos. ................................................................................. 52

Figura 6.Flujo en equilibrio en el elemento. ................................................................................ 54

Figura 7. Localización de las estaciones climatológicas. ............................................................. 56

Figura 8. Relación de la precipitación y evaporación total anual (mm) de la estación Univ Fco P

Santand.......................................................................................................................................... 57

Figura 9. Relación de temperatura y humedad de la estación Univ Fco P Santand. .................... 57

Figura 10. Relación de la precipitación y evaporación de la estación Abrego Centro Admo. .... 58

Figura 11. Relación de temperatura y humedad de la estación Abrego Centro Admo. ............... 58

Figura 12. Relación de la precipitación y evaporación de la estación Playa La. ......................... 59

Figura 13. Relación de temperatura y humedad de la estación Playa La. .................................... 59

Figura 14. Parámetros morfométricos de la cuenca. .................................................................... 62

Figura 15. Ancho máximo de la cuenca. ...................................................................................... 63

Figura 16. El mapa representa la orientación de la cuenca a través de colores. .......................... 64

Figura 17. Mapa de pendientes (%) de la cuenca. ........................................................................ 67

Figura 18. Análisis hipsométrico de una cuenca. Fuente: (Campos Aranda, 1998) .................... 68

Figura 19. Mapa de elevaciones cada 100 metros. ....................................................................... 69

Figura 20. Curva hipsométrica. .................................................................................................... 70

Figura 21. Histograma de frecuencias altimétricas. ..................................................................... 70

 xiii

Figura 22. Orden de los cauces (Horton, 1945). .......................................................................... 72

Figura 23. Localización de las estaciones de monitoreo. ............................................................. 85

Figura 24. Distancia entre las estaciones en el primer muestreo. ................................................ 90

Figura 25. Distancias entre estaciones en el segundo muestreo. .................................................. 90

Figura 26. Caudal de las estaciones en los dos periodos de muestreo. ........................................ 95

Figura 27. Temperatura de las estaciones en los dos periodos de muestreo. ............................... 96

Figura 28. pH de las estaciones en los dos periodos de muestreo. ............................................... 97

Figura 29. Oxígeno disuelto en las estaciones en los dos periodos de muestreo. ........................ 98

Figura 30. La DBO rápida de las estaciones en los dos periodos de muestreo. ........................... 98

Figura 31. DBO lenta de las estaciones en los dos periodos de muestreo. .................................. 99

Figura 32. Solidos suspendidos inorgánicos de las estaciones en los dos periodos de muestreo. 99

Figura 33. Coliformes totales de las estaciones en los dos periodos de muestreo. .................... 100

Figura 34. Amonio de las estaciones en los dos periodos de muestreo. .................................... 101

Figura 35. Nitratos de las estaciones en los dos periodos de muestreo. ..................................... 101

Figura 36. Conductividad de las estaciones en los dos periodos de muestreo. .......................... 102

Figura 37. Detritos (MOP) de las estaciones en los dos periodos de muestreo. ........................ 102

Figura 38. P inorgánico de las estaciones en los dos periodos de muestreo. ............................. 103

Figura 39. Alcalinidad de las estaciones en los dos periodos de muestreo. ............................... 103

Figura 40. Topología del río Algodonal. .................................................................................... 104

Figura 41. Ejemplo de la hoja Qual2k........................................................................................ 105

Figura 42. Ejemplo de la hoja Headwater. ................................................................................. 105

Figura 43. Ejemplo de la hoja Reach. ........................................................................................ 106

Figura 44. Ejemplo de la hoja Air Temperature......................................................................... 106

 xiv

Figura 45. Ejemplo de la hoja Cloud clover............................................................................... 107

Figura 46. Ejemplo de la hoja Shade.......................................................................................... 107

Figura 47. Contantes de calibración proporcinadas por el modelo. ........................................... 107

Figura 48. Ejemplo de la hoja Point Sources. ............................................................................ 108

Figura 49. Ejemplo de la hoja Hidraulis data. ............................................................................ 109

Figura 50. Ejemplo de la hoja Temperatura data. ...................................................................... 109

Figura 51. Ejemplo de la hoja WQ data. .................................................................................... 110

Figura 52. Celdas de calibración de la hoja WQ data. ............................................................... 110

Figura 53. Tiempo de viaje río Algodonal (m3/s) en la campaña 1. .......................................... 111

Figura 54. Oxígeno disuelto OD río Algodonal (mgO2/L) en la campaña 1. ............................ 112

Figura 55. Temperatura del río Algodonal (°C) en la campaña 1. ............................................. 112

Figura 56. Potencial de hidrógeno pH río Algodonal en la campaña 1. .................................... 113

Figura 57. Tiempo de viaje río Algodonal (m3/s) en la campaña 2. .......................................... 113

Figura 58. Oxígeno disuelto OD río Algodonal (mgO2/L) en la campaña 2. ............................ 114

Figura 59. Temperatura del río Algodonal (°C) en la campaña 2. ............................................. 114

Figura 60. Potencial de hidrógeno pH río Algodonal en la campaña 2. .................................... 115

Figura 61. Comportamiento del OD frente a los escenarios prospectivos. ................................ 117
 xv

Introducción

La presente investigación hace referencia al tema de autodepuración de cauces naturales,

que se puede definir como un fenómeno natural complejo que se da en el medio acuático, y que

se puede atribuir a la población microbiana (tanto de origen animal como vegetal) que utilizan y

mineralizan las sustancias orgánicas aportadas por los vertimientos (Molero, Saez, & Soler,

1981).

La principal característica de este fenómeno es la acción del oxígeno (Demanda Biológica

de Oxígeno y Oxígeno Disuelto). Por lo que cuando hay suficiente cantidad de oxígeno, las

sustancias orgánicas se oxidan convirtiéndose en materia mineral o son destruidas por

organismos aerobios. Asimismo, las características del vertimiento, del caudal y de la

oxigenación (turbulencia) del agua influyen en la capacidad de autodepuración del río.

En la actualidad podemos evidenciar los problemas ambientales derivados de la

inadecuada recolección, tratamiento y disposición de los vertimientos generados por las

diferentes actividades antrópicas. Actividades que originan de forma excesiva, problemas tanto

de salud pública como de calidad del recurso hídrico (Terry, 2007) y (Duarte, Pico, & Estupiñan,

2014). Como consecuencia de esto, la situación ambiental de los cuerpos hídricos del país y a

nivel mundial es insostenible, puesto que a diario son vertidos más del 80% de las aguas

residuales con altas cargas contaminantes a nuestros cuerpos de agua (WWAP, 2017),

disminuyendo o eliminando totalmente la capacidad de asimilación y autodepuración de los ríos,

arroyos y quebradas.

La investigación de esta problemática se basa en la capacidad de autodepuración del río

Algodonal, ya que éste es una de las principales fuentes de agua para el municipio de Ocaña y
 xvi

sin embargo, la falta de planificación permite que el río Algodonal sea receptor del agua residual

del municipio de Ábrego, del corregimiento de la Ermita (Ocaña) y de algunos otros afluentes

que presentan un tratamiento ineficiente, generando así una gran descarga de agua residual

doméstica directamente al río (Alcaldía Ocaña, 2016-2019).

Profundizar sobre este tema es de interés académico, con el fin de fortalecer el ámbito

profesional en el área de gestión integral del agua y gestión de residuos líquidos.

En función de la autodepuración, el comportamiento del oxígeno se representa a través de

modelos que representan la transformación física, química y biológica en los cuerpos de agua.

De esta manera, la investigación inició con la recolección de información primaria y secundaria

sobre las características del río Algodonal y con la identificación de los parámetros

morfométricos de la cuenca a través de un Modelo de Elevación Digital (DEM). Luego se

realizaron aforos a las estaciones de monitoreo localizadas a lo largo del río Algodonal y se

analizaron las características fisicoquímicas y microbiológicas, según los tiempo de viaje de la

masa de agua. Finalmente se determinó la capacidad de autodepuración con el programa de

modelación seleccionado.

Durante la investigación en campo, uno de los obstáculos fue acceder a las estaciones de

monitoreo en las horas correspondientes del tiempo de viaje, y realizar el aforo del caudal en las

horas de la noche.

Teniendo en cuenta lo planteado anteriormente, el objetivo de este trabajo es conocer la

capacidad máxima de asimilación del río Algodonal frente a la descarga de vertimientos de aguas

residuales, para promueve la planificación del recurso hídrico en función de la dinámica real del

río Algodonal, como por ejemplo, la capacidad de remoción que debe tener las Plantas de
 xvii

Tratamiento de Aguas Residuales y los Planes de Ordenamiento de Cuencas Hidrográficas, entre

otros que permitan ejecutar acciones claras sobre la recuperación y la conservación del recurso

hídrico, y que se identifiquen y controlen los factores que afecten el ecosistema.

 1

Capítulo 1 Capacidad máxima de asimilación del río Algodonal por la descarga de

vertimientos de aguas residuales entre los municipios de Ábrego y Ocaña, Norte de

Santander

1.1. Problema de investigación

En la actualidad, existe la necesidad de proteger y garantizar la sostenibilidad y el

aprovechamiento de la gestión integrada del recurso hídrico, debido a la inminente crisis del

agua que ocasiona el crecimiento poblacional y la contaminación (UNDP-SIWI, 2009).

La planificación de este recurso se está viendo afectada por las presiones sobre la oferta

hídrica proveniente del uso, aprovechamiento, administración y demandas de diferentes actores,

y a aún ante los esfuerzos por instaurar políticas de agua en el país, Colombia evidencia una

crisis de gobernabilidad que no favorece la gestión integral del recurso hídrico (Zamudio, 2012),

convirtiendo las fuentes hídricas en receptores de cargas contaminantes.

Así mismo, el río Algodonal es receptor de vertimientos y a su vez, es fuente de

abastecimiento de agua para el consumo humano de la población del municipio de Ocaña, lo cual

nos lleva a preguntarnos cuál es la carga contaminante que puede asimilar el río Algodonal entre

los municipios de Ábrego y Ocaña.

1.2. Planteamiento del problema

Tener acceso al agua de calidad es esencial para el desarrollo de una población y para

garantizar protección al medio ambiente (Vilanova, Santìn, & Pedret, 2017). Organismos

internacionales reconocen que el manejo adecuado de las aguas residuales, mejora las

condiciones de salubridad, promueve el desarrollo de la población, preserva al recurso hídrico y

a los ecosistemas (UN Water, 2014). Del mismo modo, el control de la contaminación del agua y
 2

el tratamiento de aguas residuales se ha convertido en un requisito indispensable para varios

países en desarrollo y los países industrializados que buscan mantener la calidad del agua en

condiciones óptimas (Enderlein, Enderlein, & Williams, s.f.). En el caso de Colombia, el

problema de las aguas residuales se debe a:

La inadecuada o inexistente recolección, tratamiento y disposición de los vertimientos

generados por actividades como la agricultura y la industria, y por las aguas residuales de

origen doméstico, las cuales se están generando de forma sucesiva y elevada, trayendo

como consecuencia problemas de salud pública y calidad de agua. (Ministerio de

Ambiente, Vivienda y Desarrollo Territorial, 2004, pág. 5)

Esta problemática comienza a comportarse de manera insostenible, debido a que las

fuentes hídricas receptoras de vertimientos están alcanzando su capacidad de asimilación.

El río Algodonal es una de las Subzonas hidrogeográficas del Río Catatumbo (HIMAT,

1978) e (IDEAM, 2013); posee una extensión de 74.639,8 hectáreas (CORPONOR, 2011) y es

una de las principales fuentes hídricas del departamento Norte de Santander. Tanto el municipio

de Ocaña, como La Playa y Ábrego, tienen su área urbana incluida en esta cuenca, permitiendo

así dimensionar la importancia de todos los procesos económicos, sociales y ambientales que se

generan en estos centros poblados y que causan impactos ambientales a este sistema hídrico

(POMCHRA, 2010).

De acuerdo a una síntesis ambiental que realizó la Corporación Autónoma Regional de la

Frontera Nororiental en el año 2011 se estimó la oferta y la demanda del Río Algodonal para los

tres municipios de su jurisdicción (Ábrego, La Playa y Ocaña). En su informe se determinó que

el mayor aportante de agua a la cuenca del Río Algodonal es el municipio de Ábrego, debido a
 3

sus altos niveles de precipitación sobre su superficie: por otro lado, el municipio de La Playa

posee menos área de escurrimiento y precipitaciones; en cuanto al municipio de Ocaña, la oferta

de agua depende de los niveles de precipitación en su área. En cuanto a la demanda del recurso,

el municipio de Ábrego tiene el consumo más alto debido a la utilización del agua para las

actividades agrícolas para las cuales se emplean excesivas cantidades. Para el municipio de La

Playa, la demanda es poco significativa dentro del uso y/o consumo del recurso. Ocaña por su

parte, al poseer mayor cantidad de habitantes en relación a los anteriores municipios, su demanda

hídrica equivale el sector doméstico y el agrícola.

Uno de los mayores impactos negativos al Río Algodonal se da por la no prestación de un

servicio de alcantarillado en la población rural, y en cuanto a la zona urbana, se da por la

necesidad de la reposición e instalación de una red de alcantarillado en lugares determinados y

por la ineficiencia del sistema actual de tratamiento de aguas residuales, conformado por la

laguna de oxidación y la PTAR. A esto se le suma las descargas de aguas residuales generados

por la planta de beneficio animal municipal (Alcaldía Ábrego, 2016-2019). Quizás, el aspecto de

mayor relevancia sea la falta de caracterización de los pesticidas y los residuos domésticos del

área rural.

Aguas abajo se encuentra la captación principal de agua para el abastecimiento de agua

potable del municipio de Ocaña por la empresa de servicios públicos ESPO (Alcaldía Ocaña,

2016-2019). Ésta se ve amenazada por el aumento de las descargas por vertimientos aguas arriba,

como los vertimientos generados por el corregimiento de La Ermita.

Con el fin conocer la calidad del agua del río Algodonal y simular su dinámica natural, se

pretende analizar las características actuales del recurso hídrico, las cargas contaminantes de los

vertimientos más significativos al río Algodonal y para realizar una modelación que permita
 4

predecir o simular el oxígeno disuelto y algunos parámetros fisicoquímicos y microbiológicos

que representen el comportamiento del sistema del río Algodonal y permita tomar decisiones

sobre la planificación y control de vertimientos en el recurso hídrico.

1.3. Formulación del problema

¿Cuál es la capacidad máxima de asimilación del río Algodonal por la descarga de

vertimientos de aguas residuales entre los municipios de Ábrego y Ocaña?

 5

1.4. Objetivos de investigación

1.4.1. Objetivo general.

Determinar la capacidad máxima de asimilación del río Algodonal por la descarga de

vertimientos de aguas residuales entre los municipios de Ábrego y Ocaña, Norte de

Santander.

1.4.2. Objetivos específicos.

 Seleccionar un modelo de calidad del agua en el río Algodonal mediante información

primaria y secundaria sobre las características naturales de río y a partir de la

información disponible.

 Analizar la calidad del cuerpo de agua y las cargas contaminantes de los vertimientos al

río Algodonal con el fin de compararlos con los límites permisibles a nivel nacional.

 Determinar la capacidad de autodepuración de la fuente hídrica modela da desarrollando

el programa de calidad del agua.

 6

1.5. Justificación

Desde el año 1984 en Colombia se contempla la Gestión del Recurso Hídrico (Decreto

N° 1594), pero no se tiene en cuenta la capacidad máxima de asimilación que tiene los diferentes

cuerpos de agua, ni se evalúa previamente la afectación que se puede presentar por los

vertimientos y demás usos potenciales (Alean Vera, 2009). Este problema de planificación,

puede ocasionar que exista una presión insostenible sobre una fuente hídrica que puede conducir

a su desaparición.

Los aspectos mínimos del Ordenamiento del Recurso Hídrico se reglamentan en el

artículo 6 del Decreto N° 3930 de 2010, donde la autoridad ambiental deberé, entre otros “aplicar

y calibrar modelos de simulación de la calidad del agua, que permitan determinar la capacidad

asimilativa de sustancias biodegradables o acumulativas y la capacidad de dilución de sustancias

no biodegradables y/o utilización de índices de calidad del agua, de acuerdo con la información

disponible”. Estos modelos de simulación de los impactos que cause cualquier vertimiento en el

cuerpo de agua y/o al suelo, deben ser en función de la capacidad de asimilación y dilución del

cuerpo de agua receptor, y de los usos y criterios de calidad establecidos en el Plan de

Ordenamiento del Recurso Hídrico (PORH), y serán responsabilidad de los generadores de

vertimientos.

La falta del conocimiento científico sobre el recurso hídrico, es el panorama que limita la

creación de acciones en torno al ordenamiento, que permitan apoyar el manejo del recurso y que

integre las tres dimensiones del desarrollo sostenible: el aspecto social, económico y ambiental,

con el fin de conservar los ecosistemas (Alean Vera, 2009). Por esto, es necesario entender que

la fuente hídrica además de ofrecer el agua para el consumo humano y para las actividades
 7

productivas, se hace importante mantener un remanente de agua para la preservación de la

biodiversidad, productividad y estabilidad (Domínguez, Gonzalo, Vanegas, & Moreno, 2008).

El aporte más significativo de este trabajo es determinar el comportamiento del río

Algodonal ante la descarga de los vertimientos de aguas residuales con altas cargas de materia

orgánica (Grupo de Investigación en Asustos Ambientales y Urbanos [MINDALA], 2017),

donde se identifique la cantidad de oxígeno disuelto presente en el cuerpo receptor. Esto permite

conocer la capacidad máxima de asimilación del río Algodonal, de tal modo que cumpla con los

estándares de calidad definidos por la autoridad ambiental competente.

El presente proyecto es una aproximación del estado actual del Río Algodonal, teniendo

en cuenta la implementación de nuevas tecnologías y normatividad actual vigente que favorezcan

la preservación del recurso, y que permita entender la dinámica natural del río y que sirva de

herramienta para proponer acciones en pro de los cambios en la calidad que este pueda presentar.

1.6. Formulación de la Hipótesis

1.6.1. Problema.

La disponibilidad de agua se enfrenta a fuertes presiones por la descarga de

contaminantes (WWDR-UNESCO, 2015).

El río Algodonal se enfrenta a presiones por la descarga directa de aguas residuales,

principalmente del municipio de Ábrego y del corregimiento de la Ermita, y algunas indirectas

de sus tributarios como la descarga de aguas residuales en la quebrada El Playón del municipio

de La Playa de Belén. Sumado a esto, existen sistemas de capacitación legal e ilegal a lo largo

del río, debido al aumento de la demanda del recurso hídrico para consumo humano y las
 8

actividades agropecuarias de la región, lo que ocasiona que el caudal disminuya a lo largo del

gradiente longitudinal del río de manera significativa.

1.6.2. Pregunta de investigación.

¿Cuál es la capacidad máxima de asimilación del río Algodonal frente a la descarga de

vertimientos de aguas residuales entre los municipios de Ábrego y Ocaña, Norte de Santander?

1.7. Delimitaciones

1.7.1. Delimitación operativa.

Durante la ejecución de este proyecto se realizarán técnicas y procedimientos que se

ajustan a cada una de las actividades que han sido programadas para lograr los objetivos

propuestos, para la cual se hace necesario:

 Reconocimiento de la zona de estudio.

 Georreferenciación de las estaciones.

 Realizar la cartografía de la zona.

 Selección del modelo.

 Socialización con las partes interesadas.

 Recolección de la información primaria.

 Aforos y muestreo.

 Análisis de las muestras de agua en el laboratorio.

 Calibración y validación del modelo.

 Modelación.

 Análisis e interpretación.
 9

1.7.2. Delimitación conceptual.

En la elaboración de este proyecto se tendrán en cuenta los conceptos más pertinentes

para la contextualización del proyecto y que a su vez permitan justificar la investigación.

Además, serán el soporte técnico y facilitara la comprensión de la metodología utilizada.

 Agua.

 Administración del agua.

 Autodepuración de cauces naturales.

 Demanda Bioquímica de Oxigeno.

 Oxígeno Disuelto.

 Calidad del agua.

 Contaminación del agua.

 Vertimientos.

 Gestión Integral del Recurso Hídrico.

 Conceptos de modelación.

 Verificación del modelo de calidad del agua.

 Calibración del modelo de calidad del agua.

 Validación del modelo de calidad del agua.

 Modelación de la calidad del agua.

1.7.3. Delimitación geográfica.

El Área de Influencia Directa (AID) abarca la parte media del Río Algodonal, entre los

municipios de Ábrego y Ocaña.

 10

La Cuenca Hidrográfica del río Algodonal se encuentra ubicada en el departamento Norte

de Santander, en los municipios de Ábrego, La Playa, Teorama, Ocaña y dos municipios del

departamento del Cesar: Gonzáles y Río de Oro (POMCHRA, 2010). El río Algodonal está

situado en la cordillera oriental, y al occidente del departamento Norte de Santander. En la zona

norte limita con las cuencas de los ríos Catatumbo y río de Oro (Cesar), en el oriente con la

cuenca del río Tarra, en la zona sur y occidente limita con las cuencas del río Lebrija Regidor, la

cuenca de Río de Oro (Cesar) y la del Río San Alberto (Cesar) (POMCHRA, 2010, p. 2).

Tabla 1
Coordenadas de la cuenca del río Algodonal

Planas Geográficas
Delimitación
X Y Latitud Longitud

Norte 1083090 1421230 08°24’26,10563” 73°19’35,41165”

Sur 1094355 1359565 07°50’58,35814” 73°13’31,48436”

Este 1073670 1402000 08°14’00,76747” 73°24’44,36086”

Oeste 1102965 1372465 07°57’57,60893” 73°08’49,52800”

Fuente: Grupo Diagnóstico POMCHRA, Río Algodonal 2006.

Figura 1. Mapa de localización del río Algodonal.

Fuente: autores del proyecto.

 11

1.7.3.1. Municipios a los que pertenece.

El área de influencia de la cuenca corresponde a cuatro municipios de la Provincia de

Ocaña, Ábrego, Teorama y la Playa de Belén.

Tabla 2
Distribución porcentual del río Algodonal por municipio

Porcentaje de área municipal

Municipio Área total municipal
dentro de la cuenca
Ábrego 137.819,3 30%

La Playa de Belén 24.109,1 6.1%

Teorama 93.010,1 0.2%

Ocaña 64.223,3 42,7%

Fuente: Grupo Diagnóstico POMCHRA, Río Algodonal 2006

El río Algodonal nace en la confluencia de río Frío y río Oroque, y estos a su vez nacen

en el Páramo de Jurisdicciones. El área total del río Algodonal es de 74.639,8 hectáreas (has),

siendo su longitud de 62,7 kilómetros. Su altitud varía desde 3.680 m.s.n.m. hasta los 1.200

aproximadamente.
 12

Figura 2. La confluencia del río Oroque y río Frío da origen al río Algodonal.
 13

1.7.4. Delimitación espacial.

La unidad de análisis corresponde a la selección de 5 estaciones de muestreo,

seleccionadas de acuerdo a la facilidad de acceso y según las zonas de vertimientos más

significativas.

Cabecera: corresponde a aguas arriba de la descarga de la laguna de estabilización del

municipio de Ábrego.

Estación 1: corresponde a 100 metros aguas abajo de la descarga de agua residual de la

laguna de estabilización del municipio de Ábrego y de un afluente.

Estación 2: cerca de la estación de servicio de Terpel.

Estación 3: Trituradora Guayabal E.A.T., aguas abajo del efluente del municipio de la

Playa de Belén.

Estación 4: Algodonal Club, aguas debajo de los vertimientos del corregimiento de la

Ermita.

Estación 5: en la captación de E.S.P.O. S.A.

1.7.4. Delimitación temporal.

El tiempo de ejecución del proyecto es de 4 meses, contados a partir del 28 de febrero del

año 2018, fecha de la aprobación de la propuesta de investigación y concluye en el mes de Julio

con la publicación y sustentación del documento final. Los parámetros serán analizados en dos

campañas de monitoreo, durante el periodo más crítico del rio (época seca) y durante el periodo

de lluvias, en cada uno de las estaciones que fueron estratégicamente definidos a lo largo del rio

Algodonal.
 14

Capítulo 2 Marco referencial

2.1. Marco histórico

El agua es el recurso natural más utilizado por el hombre a lo largo de su historia. El

desarrollo de los asentamientos humanos se ha realizado históricamente en las zonas cercanas a

los ríos, ya que eran fuente de abastecimiento de agua para el consumo doméstico y sus

actividades económicas. Pero al mismo tiempo son zonas de vertimientos de desechos sólidos y

aguas residuales, lo que genera importantes problemas de contaminación que inciden sobre el

medio ambiente y en la calidad de vida de las personas.

El evidente deterioro de las fuentes hídricas y al tratarse de un bien escaso e

indispensable para la vida de los seres vivos, ha permitido considerar el agua como un

patrimonio de la humanidad (Romo, 2003). En la actualidad, la preservación, la conservación y

la defensa de los recursos naturales son una necesidad del hombre y una obligación jurídica, por

lo que es necesario realizar estudios sobre las causas, los efectos y posibles programas de

prevención a los problemas ambientales causados por las actividades antrópicas. En el caso de

las fuentes hídricas cualquier proyecto para la protección de los ríos, quebradas o arroyos, o para

definir cómo debe tratarse el agua residual y su nivel de eficiencia, se debe realizar en función de

la capacidad de autodepuración natural, es decir que se identifique la capacidad que tenga la

fuente de asimilar los vertimientos de aguas y restituir si calidad.

En el intento por describir este fenómeno, Streeter & Phelps (1925) desarrollaron el

primer modelo matemático que toma en consideración los mecanismos principales de transporte

de contaminante y la depuración natural mediante el estudio del oxígeno disuelto y la demanda

bioquímica de oxígeno (Castillo, 2004). Su discusión teórica era demostrar que la capacidad de

una corriente para recibir y oxidar las aguas residuales depende sobre sus recursos de oxígeno
 15

(Streeter & Phelps, 1958). Del mismo modo, Michael Jacobs en su libro sobre la economía

verde, explica que la capacidad que tenga un tramo de río de asimilar residuos orgánicos depende

de la cantidad y del tipo de microorganismos que contenga, conjuntamente con otras

características tales como el caudal y la turbulencia (Jacobs, 1996). Tal como sucede con un

recurso renovable, la pureza del medio receptor permanecerá constante siempre que la tasa de

descargas no exceda a la tasa de absorción del medio receptor.

Actualmente los sistemas naturales están al límite de su capacidad de autoregeneración

como consecuencia de una población creciente que cada vez produce y consume más (Contreras

& Molero, 2012). Es por esto que en el siglo XVIII el interés por la depuración y el tratamiento

de las aguas para uso público alcanza un grado importante que continúa hasta nuestros días,

sobre todo en los países desarrollados, donde es posible aplicar el mecanismo de depuración y

tratamiento más eficaces, pero a su vez más costosos.

Los modelos hidrodinámicos y de calidad del agua integran insumos de fuentes puntuales

y no puntuales para determinar los impactos sobre la calidad del agua en los cuerpos de agua

receptores. Y permiten a los científicos e ingenieros determinar las capacidades de asimilación

del cuerpo de agua, determinar el nivel de mejores prácticas de gestión o predecir el tiempo

requerido para que un sistema se recupere después de haber sido alterado (United States

Environmental Protection Agency (EPA), 2015). Esta herramienta ha sido utilizada en muchos

países para la determinación de la capacidad de autodepuración de los naturales y artificiales.

Entre las investigaciones que se realizaron en Colombia se pueden mencionar a:

Castro Huertas (2015) quien aplicó el modelo QUAL2Kw en la modelación de la calidad

del agua del río Guacaica-Manizales, con el fin de mejorar la planificación del recurso hídrico y

evitar los conflictos por su uso. Se seleccionaron 8 estaciones de monitoreo sobre el río y 4
 16

estaciones adicionales sobre la quebrada Olivares (tributario de Guacaica) con el fin de conocer

la influencia y presión que se ejerce sobre este cuerpo de agua, a causa de vertimientos y

actividades económicas que afectan la corriente durante su paso por el municipio de Manizales.

Se alimentó el modelo de calidad de agua QUAL2Kw con parámetros fisicoquímicos y

microbiológicos tomados en campo en las campañas de monitoreo, descripción de los tramos del

río, constantes hidráulicas, selección de las constantes cinéticas de calibración, condiciones

meteorológicas, y la información fisicoquímica correspondiente a fuentes puntuales y difusas.

Los resultados permitieron evidenciar un tramo crítico sobre la corriente, el cual tiene problemas

de contaminación desde su nacimiento debido a las actividades antrópicas que se realizan en la

zona, sin embargo, el impacto más representativo que se observa es en el punto en el cual el río

recibe a la quebrada Olivares con una alta carga contaminante, lo cual ocasiona que los niveles

de todos los contaminantes se incrementen considerablemente y deterioren la calidad del agua

del río Guacaica en ese tramo (Castro Huertas, 2015).

Del mismo modo en la ciudad de Bogotá, se aplicó el modelo QUAL2Kw de la calidad

del agua del río Tunjuelo en dos escenarios de implementación del plan de saneamiento, a partir

de información secundaria recopilada en diferentes entidades. En la metodología utilizaron el

protocolo de modelación propuesto por el proyecto HarmoniQua que consta de 5 fases: Plan de

modelación, datos y conceptualización, preparación del modelo, calibración y validación y, por

último, simulación y evaluación. La metodología se desarrolló como estaba prevista, en la

calibración buscaron un número de corridas que permitieran la convergencia del coeficiente

utilizado (Nash-Sutcliffe), y se dejan los valores recomendados en el modelo QUAL2kw,

excepto, el número de poblaciones np y las generaciones ngen (Castaño Ossa, 2015).

 17

En el año 2016 (Jiménez & Cortes, 2016) realizaron un trabajo de investigación sobre la

calidad del agua del río Chicamocha a través del programa de modelación Qual2kw. Usaron los

monitoreos al río Chicamocha por parte de CORPOBOYACA en el año 2015, y desarrollaron el

proyecto partiendo del análisis de información y aplicación de índices de calidad de agua;

segundo aplicaron dos modelos, el estático y dinámico; y tercero, simularon escenarios de

remoción de la planta de tratamiento de agua residual. El resultado más relevante del proceso de

modelación fue la identificación del foco de contaminación en el tramo que concerniente al

vertimiento proveniente de Tunja y la generación de los respectivos escenarios de remoción de la

PTAR; además de esto otro resultado fue la modelación estática que permite observar el

comportamiento de los parámetros monitoreados en el río. Sin embargo, la modelación dinámica

no arrojo los resultados esperados por falta de datos.

2.2. Marco conceptual

2.2.1. El agua.

El agua es el elemento más importante para la vida en la Tierra, ya que el 70 % del

planeta y de los seres vivos están constituidos por agua, esencial para la producción agrícola y la

seguridad alimentaria. Sin embargo, el agua está disminuyendo a un ritmo alarmante, siendo este

un reto para lograr el desarrollo sostenible.

El agua es un recurso natural renovable que ha fluido por toda la tierra en su proceso “el

gran ciclo hidrológico”, pero se ha venido deteriorando por la contaminación y el mal uso del

agua. En la Agenda 21 (Conferencia de las Naciones Unidas sobre Ambiente y Desarrollo Río de

Janeiro, 1992) señala entre otros que: “la escasez generalizada de recursos de agua dulce, su

destrucción gradual y su creciente contaminación, así como la implantación progresiva de

 18

actividades incompatibles en muchas regiones del mundo, exigen una planificación y una

ordenación integradas de los recursos hídricos”.

En el año 2014, el informe Global Risks se reportó la “crisis del agua” como la tercera

preocupación global más importante, debido a que más del tercio de la población a nivel mundial

viven en países que presentan insuficientes fuentes hídricas y se prevé que esta cifra se puede

duplicar en el año 2025. (World Economic Forum, 2014)

2.2.1.1. Contaminación de las aguas.

Se le llama contaminación del agua a aquella que altera los parámetros fisicoquímicos

naturales del agua. Tiene lugar en espacios físicos que alojan agua ya sean lagos, mares,

acuíferos o aguas subterráneas. Las mayores causas de contaminación del agua son producto de

acciones antrópicas, cuyos contaminantes y productos químicos deterioran directamente la

calidad del agua. La constante modificación de las propiedades bioquímicas del agua repercutirá

negativamente a los seres vivos, ya que tres cuartas partes del planeta están cubiertas por agua, lo

cual hace que ésta sea el recurso natural más importante del medio ambiente.

2.2.1.2. Calidad del agua.

La calidad del agua representa las características propias del agua, que son requisitos o

condiciones de las especies bióticas que habitan ese ecosistema. Sin embargo, los cuerpos de

agua reciben y transportan cargas contaminantes producto de diferentes actividades

socioeconómicas y que en la mayoría de ocasiones son vertidas sin previo tratamiento. Estas

acciones se están incrementando diariamente lo que afecta la calidad del agua.

La calidad del agua es una situación que preocupa a todos los países, por su repercusión

en la salud de la población, considerando como factores de riesgo a los agentes infecciosos, los
 19

productos químicos tóxicos y la contaminación radiológica. Sin embargo, la definición de la

calidad del agua es muy engorrosa dada la cantidad de variables que se utilizan para describir el

estado de los cuerpos de agua en términos cuantitativos.

2.2.2. Manejo del recurso hídrico.

El manejo racional del agua es fundamental para el desarrollo económico y social. Sin

embargo, este desarrollo entra en conflictos con la preservación del medio ambiente. Lo ideal es

considerar el desarrollo sostenible del recurso hídrico, donde el desarrollo económico se lleve a

cabo considerando la preservación del medio ambiente.

El Ministerio de Ambiente y Desarrollo Sostenible define la Administración del Recurso

Hídrico como parte de la gestión integral del recurso hídrico como “la aplicación de diversos

instrumentos técnico y normativos a través de los cuales se realiza la gestión sostenible que parte

del conocimiento sobre la disponibilidad del agua en términos de calidad y cantidad, así como su

descarga a los cuerpos receptores luego de ser utilizada en diversas actividades” (Min Ambiente,

s.f.). Debido a esto, se hace necesario estudiar los componentes básicos del recurso hídrico

como lo es la hidráulica, la hidrología y la ecología de forma sistémica, buscando un desarrollo

autosostenido (Menéndez, 2010).

2.2.3. El ciclo hidrológico.

El ciclo hidrológico es la recirculación continua del agua comprendida entre la atmósfera,

la superficie terrestre, el subsuelo y los océanos. El agua puede estar en cualquiera de sus tres

estados: sólido, líquido o gaseoso. El ciclo hidrológico comienza con la precipitación, ésta se da

en forma de lluvia, nieve o hielo. La precipitación puede caer sobre vegetación o directamente

sobre el terreno, en éste último, se reúne toda el agua precipitada causando el fenómeno de
 20

escorrentía, el cual permite que parte del agua se infiltre al suelo y que otra parte se evapore y

vuelva a la atmósfera. Los flujos superficiales causados por la escorrentía, se concentran en uno

conformando así las quebradas, arroyos y ríos que eventualmente desembocan en lagos y

océanos. Y luego se produce la evaporación en éstos cuerpos hídricos (Menéndez, 2010).

2.2.4. Autodepuración.

De acuerdo con Rueda (2005), el ciclo integral del uso del agua integra el conjunto de

procesos físicos, químicos y biológicos que sufre el agua, desde que el curso de un río se

interrumpe en una captación de agua para consumo humano, hasta que ésta es devuelta a un

sistema hídrico natural y se acomoda a la perturbación introducida por el vertimiento.

Los ríos se consideran reactores biogeoquímicos, ya que presentan cambios en la

composición o calidad del agua debido a los procesos naturales que se presentan a lo largo de su

cauce (Rueda, 2005). Estos procesos tienen la capacidad de eliminar sólidos en suspensión,

materia orgánica, cierto tipo de contaminantes, entre otros, esta capacidad se conoce como

autodepuración. De esta manera, la autodepuración se define como la capacidad que tiene una

fuente hídrica de absorber o asimilar las diferentes cargas contaminantes que recibe en su

recorrido natural (Duarte, Pico, & Estupiñan, 2014).

La autodepuración es un proceso propio de las aguas naturales, consiste en recuperar las

propiedades naturales del agua a través de mecanismos de sedimentación de partículas

suspendidas, procesos físicos y biológicos que actúan para la degradación de la materia orgánica

presente, la cual sirve de alimento a las algas, haciendo que éstas aumenten su actividad

fotosintética reoxigenando el agua y retornando el equilibrio natural.

2.2.5. Comportamiento de las sustancias en el medio acuático.

 21

Cualquier tipo de sustancia que esté en contacto con los cuerpos de agua genera

alteraciones en las características naturales de la fuente hídrica, provocando algunos fenómenos.

Entre estos tenemos la adveccion, difusión y dispersión.

La adveccion es el transporte de alguna sustancia en un fluido producto de las corrientes

o el caudal del cuerpo hídrico; la difusión es cuando una partícula en un fluido se puede mover

de un lugar a otro con respecto al movimiento de las moléculas del agua en el tiempo y la

dispersión es provocada por el movimiento del fluido a través del medio en diferentes

direcciones. En los ríos, canales y estuarios predomina este movimiento debido al cambio

significativo de velocidades. En estos casos la concentración de la sustancia puede variar, pero su

composición química se mantiene.

Y, por otro lado, las reacciones son producto de cambios en las características químicas,

físicas y biológicas que sufre el contaminante al reaccionar en el agua, donde se transforma en

otro compuesto.

En los procesos naturales de autodepuración de los ríos se distinguen cuatro zonas según

su contaminación y fase de depuración (Revista Ambientum, 2002) y (Universidade da Coruña,

2010-2011):

Zona de degradación próxima al vertido, tienden a desaparecer los organismos más

sensibles (algunos peces y algas) y se presentan otros organismos más resistentes. El agua

muestra un aspecto de suciedad, aumenta la DQO y disminuye el oxígeno disuelto. En ésta zona

se da inicio a la degradación realizada por la flora microbiana.

 22

En la zona de descomposición activa, surgen las aguas sucias, ennegrecidas, con

espumas, y con olores ofensivos. Se presenta una descomposición anaerobia que genera una

emanación de gases al aire.

La tercera, es la zona de recuperación. Retornan la vegetación y el agua se limpia. Esto

ocurre debido al aumento del oxígeno disuelto producto de la actividad fotosintética de la

vegetación, la cual ayuda a degradar los compuestos contaminantes.

Como último existe la zona de aguas limpias, en ésta se dan las características físico-

químicas naturales del agua y la presencia de animales y vegetación propia de la naturaleza del

cauce.

2.2.5. Conceptos de modelación.

Los contaminantes en el agua tienen un comportamiento dinámico, es decir, que está en

constante movimiento. El transporte de estos contaminantes ha sido estudiado con profundidad y

esto, ha permitido establecer numerosos modelos matemáticos que describen su comportamiento.

Estos modelos han mostrado ser de una gran utilidad en la planificación de recursos y en el

adecuado manejo de la mayoría de los problemas de contaminación del agua (Castillo, 2004). La

modelación de un proceso permite proponer descripciones matemáticas simples para mejorar el

análisis, la comprensión y la predicción del comportamiento de un proceso complejo

dependiendo de los diferentes parámetros que lo componen (Pérez & Torres, 2008).

Según Bedoya (2007) (Bedoya Contreras, 2007), la aplicación de un modelo matemático

que permita crear una herramienta de simulación del comportamiento, se debe realizar en tres

fases. Primero la representación conceptual, la cual abarca una abstracción gráfica del prototipo,

considera la particularidad de las propiedades geométricas que van a ser modeladas y se diseña la
 23

manera por la cual serán integrados tanto física como funcionalmente. La segunda fase hace

referencia a la representación funcional, esta vincula todas las características físicas, procesos, y

condiciones de borde en cada parte de las ecuaciones matemáticas. Esto implica la definición

precisa de cada variable y sus relaciones con todos los otros parámetros que caracterizan el

modelo o su relación entrada-salida. La última fase es la representación computacional, proceso

por el cual, el modelo funcional es traducido a fórmulas matemáticas y procedimientos

computacionales requeridos para la solución del problema.

2.2.4.1. Modelación del Oxígeno Disuelto.

En la modelación de Oxígeno Disuelto se debe tener en cuenta las fuentes y las pérdidas

de oxigeno producidas por los diferentes procesos presentes en el transporte del agua residual a

través de los alcantarillados. Entre las fuentes de oxígeno se consideran: el re aireación desde la

atmosfera y el oxígeno disuelto en vertimientos. Y entre las pérdidas de oxigeno están la

oxidación de residuos carbonáceos y los que contienen nitrógeno, y la demanda de oxigeno de

los sedimentos (Diaz, 2004).

La solución completa del modelo de OD está dada por las soluciones parciales de los

procesos descritos anteriormente.

2.2.4.2. Modelación de la Demanda Bioquímica de Oxigeno.

Según Díaz (2004), la modelación de la Demanda Bioquímica de Oxígeno es una

herramienta vital en la predicción de la calidad del agua, ya que la DBO es un parámetro que

mide la cantidad de oxígeno consumido al degradar la materia orgánica. Esta pérdida progresiva

de la DBO se modela con la ecuación de decaimiento de primer orden:

 24

𝑑𝑋
= −𝐾𝑑 𝑋
𝑑𝑡

donde 𝐾𝑑 es la tasa de descomposición de la materia orgánica y la X la concentración de

DBO en mgL-1.

2.3. Marco teórico

2.3.1. Autodepuración de cauces naturales

La capacidad de asimilación de los vertimientos de aguas residuales o autodepuración es

un tema que se ha tratado de diversas maneras y al cual se la ha buscado una solución de análisis

por medio de modelos que representan el comportamiento de los contaminantes dentro de un

cuerpo de agua, con el objeto de mejorar la planificación del recurso hídrico.

Los cuerpos de agua poseen la capacidad de recuperarse después de sufrir un suceso de

contaminación de cargas orgánicas e inorgánicas, a través de la absorción y eliminación de la

contaminación de los vertimientos que estos reciben, pero esta capacidad es limitada. La

autodepuración de los ríos se debe a la mineralización de la materia orgánica del vertimiento por

las bacterias presentes en el cauce del rio, y a su vez las bacterias son fagocitadas por los

protozoos y éstos por organismos superiores (Torres, 2013). El factor limitante del río se da

cuando el vertimiento sobrepasa la capacidad de autodepuración y comienza a presentarse

turbidez, malos olores, cambios en el pH, la disminución en las concentraciones de oxígeno

disuelto, el aumento de la demanda bioquímica de oxígeno, etc.

Uno de los parámetros más usados para seguir el proceso de autodepuración del río es la

concentración del oxígeno disuelto, ya que éste es usado por el consumo de oxigeno usado por

los microorganismos para bioxidar la materia orgánica (desoxigenación), o por el aporte de

oxigeno del aire hacia el interior del río (reoxigenación). Streeter y Phelps (como se citó en Gil,
 25

2006), estudiaron la evolución de la concentración y la simplificaron con un balance en las

etapas de desoxigenación y reoxigenación.

En el proceso de desoxigenación se disminuye la concentración de oxígeno disuelto, es

decir, existe un déficit:

𝐷 = 𝑂𝑠 − 𝑂

Dónde:

D Déficit de concentración de oxígeno disuelto.

Os Concentración de oxígeno disuelto a saturación.

O Concentración de oxígeno disuelto en el agua.

La velocidad de aumento del déficit de oxígeno disuelto es proporcional a la DBO

presente:

𝑑𝐷
( ) = 𝑘1 . 𝐷𝐵𝑂
𝑑𝑡 𝑑

La evolución de la DBO es una reacción de primer orden, por lo que su evolución sigue

su ecuación integrada:

𝐷𝐵𝑂(𝑡) = 𝐷𝐵𝑂𝑢 . 𝑒 −𝑘1.𝑡

Entonces la variación de la concentración del oxígeno disuelto en esta etapa de desoxigenación se

expresa por:

𝑑𝐷
( ) = 𝑘1 . 𝐷𝐵𝑂𝑢 . 𝑒 −𝑘1.𝑡
𝑑𝑡 𝑑

El movimiento del agua en su recorrido por el cauce provoca una transferencia gas-

liquido, la reoxigenación. De modo que la reoxigenación es contraria al déficit de oxígeno:

 26

𝑑𝐷
−( ) = 𝑘2 . 𝐷
𝑑𝑡 𝑟

El primer modelo biogeoquímico que se construyó fue el modelo desarrollado por W. B.

Streeter y E. B. Phelps en 1925, el cual se utiliza para estudiar las relaciones de Oxígeno

Disuelto y DBO en ríos, y fue obtenido a través de mediciones en campo. Este modelo utilizó

una ecuación de balance de oxígeno para una descarga puntual a una tasa constante con

condiciones hidrodinámicas estacionarias y para un tramo con geometría uniforme (Jørgensen &

Bendoricchio, 2001).

Con el fin de resolver los problemas relacionados con los contaminantes de agua, en los

años 70 se desarrollaron modelos para la calidad del agua, usando variables como oxígeno

disuelto, DBO, bacterias termoresistentes, entre otras (Trento, 2001).

2.3.2. Tipos de modelos.

Tabla 3
Tipología y clasificación de modelos

Tipo de modelo Caracterización

Modelos de investigación Usados como una herramienta de investigación.
Modelos de manejo Usados como herramienta de manejo y gestión ambiental.
Modelos de determinísticos La predicción de valores se obtiene exactamente.
La predicción de valores depende de una distribución de
Modelos de estocásticos
probabilidad.
Modelos de Las variables definidas del sistema se cuantifican por la medida
compartimientos de ecuaciones diferenciales de tiempo dependientes.
Modelos de matrices Uso de matrices en la formulación matemática.
Modelos de reduccionistas Incluye todos los detalles relevantes como sea posibles.
Modelos de holísticos Usa principios generales.
 27

Las variables que definen el sistema no son tiempo-

Modelos de estadísticos
dependientes.
Las variables que definen el sistema están en función del
Modelos de dinámicos
tiempo o del espacio.
Los parámetros se consideran en función del tiempo y del
Modelos de distribución
espacio.
Los parámetros se están dentro de una ubicación espacial y
Modelos de agrupación
temporal, se consideran como constantes.
Modelos lineales Se usan ecuaciones de primer orden.
Modelos no lineales Una o más ecuaciones de segundo orden.
Las fuentes, las variables de estado y las salidas se
Modelos causales
interrelacionan mediante relaciones causales.
Los disturbios de las entradas afectan solamente las respuestas
Modelos de caja negra en la salida.
La causalidad no se requiere.
Fuente: Jørgensen (como se citó en Mancera Pineda, Peña Salamanca, Giraldo Henao, &

Santos Martínez (2003).

2.3.3. Modelos de calidad de agua.

El agua presenta un comportamiento dinámico en sus contaminantes, debido a que estos

están en movilidad constante, ya sea por las corrientes del agua, o por la oscilación de las

moléculas y partículas coloidales. Este fenómeno de translación de contaminantes en los cuerpos

de agua ha sido analizado, dando como resultado múltiples ecuaciones matemáticas que

describen este comportamiento (Castillo, 2004). Actualmente, existen varios modelos

matemáticos, que fueron desarrollados para estudiar, analizar, y determinar el comportamiento

de uno o varios contaminantes del agua a lo largo su cauce (Anguiano Valencia, 2012). Muchos

proyectos desarrollados sobre la identificación y estimación de los efectos de la calidad y

 28

cantidad de los cuerpos de agua, permite actividades de evaluación y estudios de impacto

ambiental (Casas C. & Vargas, 2008).

La ejecución de los modelos de simulación de la calidad del agua, permiten realizar una

gestión integral del recurso hídrico porque permiten su planificación de manera efectiva, pues

con éste puede observarse fenómenos como la capacidad de auto-depuración de las corrientes y

efectos del cambio de caudales (Lozano, Zapata, & Peña, 2003). Del mismo modo, tienen

diferentes utilidades como: la utilidad predictiva, que puede ser usada en la evaluación del

impacto de las decisiones de manejo en la cuenca o en la predicción de la evolución de la calidad

del agua bajo diferentes escenarios de desarrollo y control; y la utilidad como herramienta de

planificación, en el establecimiento de objetivos de calidad ambiental y de calidad de emisiones

para alcanzar objetivos de utilidad ambiental, o en la distribución de costos ambientales

asociados al efecto de diferentes descargas.

2.3.4. Dimensiones espaciales de los modelos de sistemas.

En cuanto al ámbito espacial de la operación de los modelos matemáticos, estos pueden

ser adimensional, unidimensional, bidimensional o tridimensionales, dependiendo en gran

medida de sus atributos físicos del medio que se va a simular.

Modelos adimensionales: corresponden a aquellos que proporcionan la evolución del

sistema con el tiempo en un sistema que carece de dimensiones espaciales. Estos modelos

suponen que el medio es homogéneo en todas las direcciones.

Modelos unidimensionales: se usan para representar flujos de aguas en ríos, siendo la

dirección considerada el sentido del escurrimiento. Estos proporcionan la evolución del sistema
 29

con el tiempo a lo largo de una dirección en el espacio, considerando que el medio es homogéneo

en las otras dos direcciones.

Modelos bidimensionales: se utilizan en cauces de ríos con gran anchura, en los cuales las

concentraciones de contaminante son variables a lo largo de la ribera del río. Para este se una un

sistema cartesiano que represente el sentido del flujo y a la dimensión lateral.

Modelos tridimensionales: su aplicación corresponde principalmente a estudios de aguas

subterráneas y en sistemas más complejos de aguas superficiales. Necesitan de mayor

información que los otros modelos, y necesitan mayor tiempo computacional, por lo que su uso

se hace necesario en problemas de gran magnitud y cuando se dispone de recursos suficientes

para su aplicación.

2.3.5. Parámetros empleados en la modelación.

Los parámetros empleados en la modelación dependen de los contaminantes que afecten

el equilibrio natural de un cuerpo de agua, estos pueden ser conservativos o no biodegradables y

no conservativos o biodegradables. Los parámetros más empleados para la determinación del

grado de contaminación son: Oxígeno Disuelto (OD) y Demanda Bioquímica de Oxígeno

(DBO).

El OD se refiere a la cantidad de oxigeno presente en el agua, el cual disminuye en la

medida que aumenta el grado de contaminación. La DBO es la cantidad de oxígeno consumida

por los microorganismos presentes en el agua para estabilizar la materia orgánica, en un cierto

tiempo y a una determinada temperatura.

Shanahan P. (como se citó en García, 2008) explica que en el caso de la modelación del

oxígeno disuelto, se debe tener en cuenta el coeficiente de reaireación, en función de la

 30

temperatura y las características hidráulicas de la corriente (profundidad y velocidad). Este

parámetro puede verse afectado por otro tipo de descargas, lo que implica que estos valores

obtenidos de la calibración no sean transferibles o representativos a otras condiciones.

2.3.6. Constantes universales.

Según Jørgensen & Bendoricchio (2001) los modelos de calidad de agua necesitan tres

pasos para su correcta utilización, la verificación, la calibración y la validación.

Verificación: la verificación es una prueba de la lógica interna del modelo. Las preguntas

típicas en la fase de verificación son: ¿El modelo reacciona como se espera? ¿El modelo es

estable a largo plazo? ¿El modelo sigue la ley de conservación masiva? ¿El uso de unidades es

consistente? La verificación es en cierta medida una evaluación subjetiva del comportamiento

del modelo.

Calibración: la calibración es un intento de encontrar la mejor coincidencia entre los

datos calculados y observados por la variación de algunos parámetros seleccionados. Puede

llevarse a cabo por prueba y error o mediante el uso de un software desarrollado para encontrar

los parámetros que mejor se ajustan entre los valores observados y calculados. En algunos

modelos estáticos y en algunos modelos simples, que contienen solo unos pocos parámetros bien

definidos o directamente medidos, la calibración puede no ser necesaria.

Validación: la validación se refiere en la demostración de la exactitud (sin error

sistemático) del modelo construido, dependiendo del tipo de aplicación (Organización de las

naciones unidas para la agricultura y la alimentación (FAO), 2004) y (Universidad del Cauca,

s.f.). Sin embargo, no todos los modelos se pueden evaluar completamente, Dee (como se citó en

FAO, 2004) dice que debemos tener en cuenta cuatro aspectos para la validación de un modelo:
 31

Validación conceptual, si el modelo es representativo con el sistema estudiado.

Validación de algoritmo, si este traduce bien los conceptos del modelo a través de las

fórmulas matemáticas.

Validación de códigos informáticos, es la aplicación de fórmulas matemáticas en el

lenguaje informático

Validación funcional, es la verificación del modelo frente a las diferentes observaciones

de manera independiente.

2.3.7. Construcción de modelos dinámicos.

Carstensen et al. 1997; Jørgensen, 1994; Kowal, 1971 (como se citó en Llavador, 2002)

mencionan que Para la construcción de modelos matemáticos dinámicos, se puede llevar a cabo

las siguientes etapas:

1. Adquisición de información y conocimientos previos sobre el sistema.

2. Especificación del objeto del modelo. Formulación del problema e identificación de las

variables internas que definen el estado del sistema a representar.

3. Identificación de las variables externas que afectan al sistema (p. ej. condiciones

ambientales) y delimitación del intervalo de valores de estas dentro del cual puede operar el

modelo.

4. Establecimiento de relaciones entre las diferentes variables del modelo; tanto entre las

variables internas entre sí como entre estas y las variables externas.

5. Selección de la estructura matemática del modelo y estimación de parámetros.

 32

6. Verificación del modelo y ajuste a las mediciones preliminares de que se dispone.

7. Validación del modelo a partir de nuevas observaciones del sistema y estimación de la

incertidumbre de las predicciones.

2.4. Marco legal

El agua es el recurso natural más importante para la existencia de la vida. Es esencial para

el fortalecimiento de los ecosistemas, para el sostenimiento de los bienes y servicios ambientales

y para el desarrollo de las actividades económicas. Este recurso se está viendo amenazado por las

prácticas insostenibles de manejo y uso del agua, es por esto que la normatividad de colombiana

ha establecido una serie de políticas públicas para la reglamentación del recurso hídrico, a

continuación, se encuentras la normatividad aplicable:

Con el fin de atender el ordenamiento del recurso hídrico en el año 1984, el artículo 24

del Decreto N° 1594 reglamenta el establecimiento de modelos de simulación de la calidad,

donde deberá tener por lo menos los parámetros de DBO5, DQO, SS, pH, T, OD, Q, datos

hidrobiológicos y Coliformes.

La Constitución Política de Colombia, también llamada Constitución Verde es la ley

máxima del país, vela por los derechos de los ciudadanos y los recursos naturales del país en sus

artículos 8, 79, 80 y 95, donde establece que es obligación del estado colombiano proteger la

diversidad del ambiente, prevenir y controlar los factores de deterioro ambiental, y a los

miembros de la comunidad nacional les reconoce como deberes la protección los recursos

culturales y naturales del país, y velar por la conservación de un ambiente sano, entre otros

(Asamble Nacional Constituyente, 1991).

 33

En 1993 se establece la Ley General Ambiental de Colombia (Ley N° 99, 1993) que en su

artículo 31 sobre las funciones de las Corporaciones Autónomas Regionales ejercerá entre otras,

las funciones de evaluación, control y seguimiento ambiental de los usos del agua y demás

recursos naturales renovables, en este caso lo relacionado con el vertimiento de sustancias o

cualquier tipo de residuos a las aguas, que puedan causar daño o poner en peligro el normal

desarrollo sostenible de los recursos naturales renovables. Y en su artículo 42 y 43 de las rentas

de las corporaciones, las tasas retributivas y compensatorias por la utilización directa o indirecta

del agua, para introducir aguas negras o servidas de cualquier origen, por las consecuencias

nocivas de las actividades. Por otro lado, el gobierno nacional crea la Ley 373 de 1997 que

establece el programa para el uso eficiente y ahorro del agua.

Con la expedición del Decreto N° 3100 (2003) y un año más tarde el Decreto N° 3440

(2004) se reglamenta las tasas retributivas por la utilización directa del agua como receptor de

vertimientos puntuales en las cuencas con prioridad por sus características de calidad. Este

recaudo será invertido exclusivamente a proyectos de inversión de descontaminación hídrica y

monitoreo de calidad de agua.

La Resolución N° 1433 (2004) reglamenta el artículo 12 del Decreto 3100 sobre los

Planes de Saneamiento y Manejo de Vertimientos PSMV. El Plan deberá formularse teniendo en

cuenta la información disponible sobre calidad y uso de las corrientes, tramos o cuerpos de agua

receptores. Así mismo, en el artículo 4 de la presentación de la información se incluirá el

diagnostico, la caracterización de las descargas, las proyecciones de la carga contaminante

generada, recolectada, transportada y tratada, por vertimiento y por corriente, los objetivos de

reducción, entre otras.

 34

Con el fin de contribuir con la sostenibilidad del recurso hídrico, mediante una gestión y

un uso eficiente y eficaz, articulados al ordenamiento y uso del territorio y a la conservación de

los ecosistemas que regulan la oferta hídrica este proyecto se rige en el marco de la Política

Nacional para la Gestión Integral del Recurso Hídrico del año 2010.

En cuanto a usos del agua y residuos líquidos se expide el Decreto N° 3930 (2010)

(Decreto N° 3930, 2010). El cuál en su capítulo III establece los lineamientos relacionados con el

Ordenamiento del Recurso Hídrico, y en el artículo 6, decreta que la autoridad ambiental deberá

tener en cuenta “la aplicación y calibración de modelos de simulación de la calidad del agua, que

permitan determinar la capacidad asimilativa de sustancias biodegradables o acumulativas y la

capacidad de dilución de sustancias no biodegradables y/o utilización de índices de calidad del

agua, de acuerdo con la información disponible”, entre otros aspectos del ordenamiento del

recurso. De igual manera, en el artículo 7 el Ministerio de Ambiente, Vivienda y Desarrollo

Territorial (MAVDT) deberá expedir dentro de los próximos ocho meses, la Guía Nacional de

Modelación del Recurso Hídrico, con base en los insumos que aporte el Instituto de Hidrología,

Meteorología y Estudios Ambientales (IDEAM). Sin embargo, hasta la fecha no se ha publicado

la guía, entonces en el parágrafo del mismo artículo menciona que la autoridad ambiental podrá

seguir aplicando los modelos de simulación existentes, utilizando, por lo menos los siguientes

parámetros: Demanda bioquímica de oxígeno a cinco (5) días DBO5, Demanda química de

oxígeno DQO, Sólidos suspendidos SS, Potencial del ion hidronio pH, Temperatura T, Oxígeno

disuelto OD, Caudal Q, datos hidrobiológicos, coliformes totales y coliformes fecales. En el

artículo 8 se reglamenta la aplicación de los modelos de simulación de la calidad del agua para la

proyección de escenarios con las mejores condiciones ambientales.

 35

El artículo 12 del Decreto N° 2667 (2012) durante la consulta para el establecimiento de

la meta global de carga contaminante, “la autoridad ambiental presentará los escenarios de

metas, de acuerdo al análisis de las condiciones que más se ajusten al objetivo de calidad vigente

al final del quinquenio y la capacidad de carga del tramo o cuerpo de agua definidos a partir de

evaluaciones y/o modelaciones de calidad del agua”.

En el área de influencia del alcance de la Resolución N° 1514 (2012), sobre los Términos

de Referencia para la Elaboración del Plan de Gestión del Riesgo para el Manejo de

Vertimientos (PGGRMV), se definirá con base a los resultados de la modelación del vertimiento

en el agua, teniendo en consideración el escenario más crítico que implique el vertimiento de

aguas sin condiciones de tratamiento y una zona de seguridad que se deberá definir de acuerdo

con la capacidad de asimilación del cauce.

Finalmente en el año 2015 se expide la Resolución 631 por la cual se establecen los

nuevos parámetros fisicoquímicos y sus valores límites máximos permisibles en los vertimientos

puntuales de aguas residuales domésticas (ARD) a cuerpos de aguas superficiales y a los

sistemas de alcantarillado (Resolución N° 631, 2015).

Capítulo 3 Metodología de la investigación

3.1. Diseño de la metodología de investigación

La metodología de la investigación que se pretende desarrollar es básica, porque

incrementa el conocimiento científico a partir la revisión de teorías durante la construcción de los

marcos referenciales en temas de cómo funciona la asimilación de los vertimientos en los cauces

naturales a través de la modelación, con la finalidad de formular o corregir las teorías existentes.
 36

El enfoque de la investigación es cuantitativo porque utilizamos el método científico

donde, a través de un proceso sistémico, nos permite tener resultados medibles como los

parámetros de calidad del agua y las mediciones de las características de la fuente hídrica.

Además, es un tipo de investigación descriptivo, que observa, describe, analiza y predice el

comportamiento de los parámetros fisicoquímicos del río Algodonal.

Por lo tanto, el diseño de la investigación es no experimental, ya que no se va a controlar

las variables externas puesto que estas hacen parte del análisis de la investigación.

3.2. Población y muestra

3.2.1. Población.

La población que se va a tratar corresponde al cuerpo de agua natural del río Algodonal,

teniendo en cuenta los parámetros fisicoquímicos e hidrológicos del cuerpo de agua del río

Algodonal.

3.2.2. Muestra.

La elección de los puntos de monitoreo en el río Algodonal se realizará principalmente

sobre la zona de influencia de los vertimientos de aguas residuales más significativos que estén

en función de la topografía de la zona, la facilidad del acceso y la relación de la distancia entre

punto y punto.

3.3. Diagrama de la metodología

A continuación, se representa gráficamente la metodología:

 37

• Datos de
entrada al Selección del modelo de
modelo
• Verificación simulación

• Elección de los
puntos de monitoreo
• Obtención de
tiempos de viaje
Análisis de la calidad del agua • Campañas de aforo
y muestreo
• Interpretación de los
datos

• Calibración
• Validación Determinación de la capacidad
• Análisis de de autodepuración
resultados

CONCLUSIONES O RECOMENDACIONES

Figura 3. Diagrama de la metodología.

Fuente: autores del proyecto.

3.4. Técnicas de recolección de la información

Este proyecto de investigación se inicia con la recolección de la información preliminar,

iniciando con la búsqueda de las entidades que poseen responsabilidades sobre el río Algodonal,

y la información sobre vertimientos, concesiones y usos del agua.

Información preliminar.

Entidades que tienen responsabilidades sobre el río Algodonal e información

suministrada.
 38

IDEAM. Fortalece la capacidad tecnológica, científica, administrativa y financiera para

producir la información hidrológica, meteorológica y ambiental de manera oportuna y con la

calidad que requieran la ciudadanía, los sectores públicos y privados del país (Instituto de

Hidrología, Meteorología y Estudios Ambientales, s.f.). El IDEAM publicó la Guía para el

Monitoreo de Vertimientos, Aguas Superficiales y Subterráneas con base a la cual se realizaron

las campañas de monitoreo en el río Algodonal.

CORPONOR. Ejecuta políticas, planes, programas y proyectos sobre medio ambiente y

recursos naturales renovables (Corporación Autónoma Regional de la Frontera Nororiental, s.f.).

La corporación publicó El Plan De Ordenación y Manejo de la Cuenca Hidrográfica del río

Algodonal en el año 2010 y nos proporcionó información sobre las concesiones y vertimientos

en el río Algodonal.

UFPSO. La universidad debe, entre otros, desarrollar las capacidades institucionales

promoviendo impactos positivos a la región, el medio ambiente y la comunidad (Universidad

Francisco de Paula Santander (Ocaña), s.f.). La información relacionada sobre estudios de la

calidad del agua del río Algodonal se encuentra en el repositorio institucional.

Alcaldía de Ábrego. La alcaldía debe velar por el buen funcionamiento del sistema de

tratamiento de las aguas residuales de este municipio. Se utilizó el Plan de Desarrollo Municipal

de Ábrego 2016-2019 como fuente de información.

Alcaldía de Ocaña. La alcaldía de Ocaña brinda información sobre el saneamiento

básico del corregimiento de la Ermita. Se utilizó el Plan Básico de Desarrollo Territorial 2002-

2011.
 39

ESPO S.A. La empresa de servicios públicos de Ocaña posee el registro histórico de los

análisis fisicoquímicos y microbiológicos, que se han realizado al caudal del río Algodonal en la

captación de la empresa.

Fuentes de información utilizadas:

Para la revisión documental se tomó como referencia los trabajos de investigación, libros

y artículos científicos que se encontraban en la base de datos de Google Schoolar y Scopus,

también en algunos repositorios institucionales como el de la Universidad Nacional de

Colombia, Universidad de La Salle, Unisangil, Pontificia Uniersidad Javeriana, Universidad de

Granada, Universidad Francisco de Paula Santander y el repositorio de la Universidad Santo

Tomás. En los informes de las Corporaciones Autónomas Regionales, en páginas de la FAO,

UNESCO, IDEAM, Min. Ambiente, Naciones Unidas, y en artículos de investigación de las

revistas de como RCCP de Colombia y la revista electrónica de la Agencia de Medio Ambiente

de Cuba, entre otras.

El libro “Control de la Calidad del Agua”. Procesos fisicoquímicos., de Wlater J.

Weber Jr., donde explica la dinámica de los procesos fisicoquímicos de una masa natural de

agua o en un reactor de tratamiento de aguas residuales (Walter J., 1979).

En el artículo “Changes in the self-purification capacity of the Magdalena River due

to climate change” explica que uno de los factores que afecta la autodepuración es el

cambio climático, ya que el aumento de las temperaturas y las impredecibles precipitaciones

están modificando considerablemente la capacidad de autodepuración de los ríos (Montes,

Navarro, Domínguez, & Jiménez, 2013).

 40

Una de las referencias para la comprensión de los parámetros de la calidad del agua

es el libro “Parámetros Físicos, Carga Orgánica y Nutrientes” de Quím. María Luisa

Castro de Esparza en el año 1995, en donde se especifica la definición, la aplicación en el

agua potable, en las aguas residuales, en la vida acuática y en la determinación en terreno y

en el laboratorio de los parámetros fisicoquímicos de calidad del agua, además se exponen

los rangos de los parámetros que pueden indicar contaminación de la fuente o el valor

adecuado para garantizar la vida acuática, entre otros (Castro de Esparza, 1995).

En el proyecto de investigación “Modelación de la Calidad del Agua en Corrientes

Naturales Mediante Sistemas No Lineales Acoplados” se presenta una primera

aproximación para el desarrollo de una herramienta para la modelación acoplada de la

calidad del agua, mediante ecuaciones no lineales. Además, se especifican algunos

conceptos de la modelación que facilitan el análisis de la información (García, 2008).

Análisis simplificado de oxígeno disuelto en el río Ubaté por el modelo qual2k

(Santamaría, 2013).

Informes de la síntesis ambiental de Norte de Santander, archivos de concesiones de

la Corporación autónoma regional de la Frontera Nororiental.

La caracterización morfológica de la cuenca del río Algodonal obtenida a través de

un Modelo de Elevación Digital DEM procesada con el software ArcGis 10.5.

Estudios fisicoquímicos y microbiológicos del río Algodonal realizados por la

empresa Servianalítica Profesionales S.A.S.

Datos obtenidos en campo (aforo de caudal).

 41

Observaciones in situ (propiedades organolépticas del río Algodonal en las

estaciones de muestreo).

Protocolo para el monitoreo y seguimiento del agua (IDEAM, 2007)

La información recolectada en esta etapa es la base para la selección de los parámetros

necesarios para medir la capacidad de autodepuración del río Algodonal. Además, se

desarrollarán las actividades programadas en el cronograma siguiendo con el orden de los

objetivos así:

3.5. Recolección de la información climatológica

Se consultó la información del catálogo nacional de estaciones meteorológicas del

IDEAM correspondiente al periodo comprendido del año 1997 al año 2017, ya que de los años

anteriores se presentaba datos incompletos. También se verificó que las estaciones estuvieran

activas. En la tabla 4 se discriminan las estaciones identificadas y usadas en el presente estudio.

Tabla 4
Estaciones climatológicas del IDEAM.

COORDENADAS
CÓD. DE
NOMBRE CLASE CATG MPIO CORRIENTE
ESTACIÓN
LATITUD LONGITUD
ABREGO
16055040 Meteorológica CP 8.087222 -73.223056 Ábrego Algodonal
CENTRO ADMO
La
16055060 PLAYA LA Meteorológica CO 8.2175 -73.235 Algodonal
Playa
UNIV FCO P
16055100 Meteorológica CO 8.237778 -73.322222 Ocaña Algodonal
SANTAND
Nota. La información fue tomada del catálogo de estaciones del IDEAM, 2018. Fuente: Autores del proyecto.

3.6 Morfometría de la cuenca

Para la realización de los parámetros morfométricos se utilizó como insumo Modelos de

Elevación Digital DEM e imágenes satelitales; la red hídrica del río Algodonal corregida por

(Arenas Sánchez & Díaz Jimenez, 2018) e información de la cartografía base del SIGOT. Los
 42

DEM corresponden a AP_24872_FBD_F0140_RT1 y AP_24872_FBD_F0150_RT1 (Fecha:

2010-09-25) del satélite ALOS PALSAR FBD en la plataforma www.asf.alaska.edu que abarcaba

el área de la cuenca, para realizar un mosaico (Dataset: © JAXA / METI, 2010). La imagen satelital

de la cuenca fue descargada de la plataforma SASPlanet del satélite de BING con un zoom Z20.

La cartografía del SIGOT fue la base para la localización espacial de la cuenca a nivel nacional.

Las curvas hipsométricas, el mapa de pendientes y el TIN cada 10 metros se generaron a partir

del DEM.

3.7. Selección del modelo de simulación

Para la revisión de la literatura, se recopilará información primaria y secundaria

relacionada con la autodepuración de cauces naturales y modelos de simulación, junto con

información de las características climatológicas presente en el catálogo nacional de estaciones

del Instituto de Hidrología, Meteorología y Estudios Ambientales IDEAM e información

morfométrica de la cuenca del río Algodonal a través de un Modelo de Elevación Digital y

procesada en ArcGIS 10.5.

Se revisó los tipos de modelos de simulación de la calidad de agua que permitan

determinar la capacidad asimilativa de sustancias biodegradables 3.9y que estén certificados,

para cuantificar el movimiento y la concentración de contaminantes en lagos, arroyos, estuarios y

ambientes marinos.

La información recolectada en esta etapa es la base para la selección de los parámetros

necesarios para medir la capacidad de autodepuración del río Algodonal.

 43

3.8. Metodología de análisis de calidad del agua

3.8.1. Ubicación de los vertimientos

Se realizó un reconocimiento en campo para buscar vertimiento de aguas residuales

domesticas alrededor de las estaciones de monitoreo, por la facilidad de acceso al río Algodonal.

Adicional a esto, se realizó un sobrevuelo de la zona de la Laguna de estabilización en Ábrego

con un drone DJI Phantom para identificar el origen del agua residual de un vertimiento aguas

debajo de la Laguna.

Los vertimientos significativos corresponden a los de un caudal Q superior a 1 litro por

segundo, según el criterio de los autores del proyecto.

3.8.2. Cargas contaminantes

El (Decreto N° 2667, 2012) en el capítulo II se define Carga Contaminante diaria (Cc)

como “el resultado de multiplicar el caudal promedio por la concentración de una sustancia,

elemento o parámetro contaminante por el factor de conversión de unidades (de mg/s a kg/h) y

por el tiempo diario de vertimiento del usuario, medido en horas por día”

𝐶𝑐 = 𝑄 𝑥 𝐶 𝑥 0,0036 𝑥 𝑡

𝑘𝑔 𝐿 𝑚𝑔
= 𝑥 𝑥 0,0036 𝑥 ℎ
𝑑í𝑎 𝑠 𝐿

Con el fin de conocer el caudal de agua residual promedio vertido al río Algodonal, se

recoleto la información descrita en la Tabla 5. En el caso de Ábrego, el ingeniero Dairo Torrado

Ascanio director de la Unidad de Servicios Públicos de Ábrego brindó la información del

número de suscriptores, el caudal de entrada de las dos captaciones, la continuidad del servicio y

confirmó que el barrio Alto Viejo posee sistema de alcantarillado pero sin conexión a la PTAR
 44

(Torrado, 2018). Para conocer la cantidad de suscriptores del Barrio Alto Viejo, se contaron las

viviendas a través de un sobrevuelo de la zona con un drone DJI Phantom y se restó el dato a los

suscriptores totales para obtener el dato de las viviendas conectadas a la Planta de Tratamiento

de Aguas Residuales PTAR. Adicionalmente, el día 22 de septiembre de 2017 se realizó, junto

con el Grupo de Investigación en Asuntos Ambientales y Urbanos MINDALA, un aforo al

caudal de entrada y salida de la Laguna, y una caracterización fisicoquímica al vertimiento

puntual de la laguna de oxidación (PTAR), en la que se determinó que el caudal de entrada y

salida era 35,3 L/s.

En el caso de la Ermita, el señor Donaciano Galvis encargado de la operación de la planta

de tratamiento de agua potable permitió medir el caudal de entrada a la planta de tratamiento de

agua e indicó las horas de continuidad del servicio y el número de suscriptores.

Tabla 5
Información de los vertimientos.

Municipio de Ábrego
Vertimiento
Vertimiento Corregimiento
Datos Casco PTAR
del B. Alto La Ermita
urbano Laguna de
Viejo
oxidación
Suscriptores 3860 57 104
Caudal de entrada a la
56 1,32
PTAP L/s
Continuidad h/día 16 12
Caudal vertido m3/s 0.0353
 45

3.8.3. Selección de las estaciones de monitoreo

La metodología para el reconocimiento de las estaciones de monitoreo se tendrá en

cuenta las consideraciones del Protocolo para el Monitoreo y Seguimiento del Agua del

(IDEAM, 2007) con el fin de obtener mediciones confiables y que permitan facilitar la logística

para el desplazamiento del personal y los equipos de hidrometría.

a) La estación debe estar situada en un tramo recto de la corriente. Cuya longitud del

tramo tendrá como mínimo 5 veces el ancho de la sección.

b) La corriente debe tener flujo uniforme y paralelo a las márgenes de la corriente.

c) La sección debe ser profunda.

d) La pendiente longitudinal del cauce debe ser uniforme, para evitar desequilibrios

de la velocidad del flujo (Manning)

e) El lecho del río debe tener geometría regular, cauce estable y no tener obstáculos.

f) Se debe evitar los lechos fangosos.

Adicional a esto, se realizó una descripción de las características de las estaciones en

función de los sentidos del observador.

3.8.4. Tiempos de viaje

El tiempo de viaje 𝑡̅ es el tiempo que tarda una sustancia disuelta en ser conducida por el

movimiento del medio y ser dispersada. Se puede estimar a partir de ensayos con trazadores o

aproximadamente como la razón entre la longitud del tramo 𝐿 y la velocidad media del agua 𝑣̅

(Autoridad Nacional de Licencias Ambientales ANLA, 2013) (p.5).

 46

𝐿
𝑡̅ =
𝑣̅

3.8.5. Aforo

La técnica de aforo fue descrita en el Protocolo de Monitoreo del río Algodonal

adoptada del Protocolo para el Monitoreo y Seguimiento del Agua (IDEAM, 2007) por los

autores del proyecto. Y en virtud de la topografía del terreno, se medirán diferentes

profundidades a lo largo de la sección transversal y su correspondiente distancia. Luego se

graficarán los datos con ayuda de un Excel para obtener el perfil transversal de cada estación de

monitoreo.

a) La verificación de que las estaciones de monitoreo cumplan con las consideraciones del

IDEAM.

b) Se inicia delimitando un transepto a lado y lado de la sección con una cuerda aforada

cada metro o a consideración del técnico (la precisión de las mediciones del caudal

depende en gran parte del número de verticales que se tomen para la ejecución de las

mediciones para el aforo, profundidad, velocidad, toma de muestras de agua, etc). Para la

realización del perfil transversal es necesario conocer las alturas en cada sub división de

la sección. Finalmente se toman las velocidades con el correntómetro en medio de cada

división, sumergido a la mitad de su altura. Nota: en caso de tener más de un metro de

profundidad se toman varias mediciones.

c) La velocidad promedio resultante de las mediciones, será tomada como la velocidad

media para la estimación de los tiempos de viaje.

d) Para la toma de muestras será necesario ubicarse en medio del río. Se llenarán los 3

recipientes, el primero de plástico con tapa rosca de capacidad de 2 litros para los
 47

parámetros fisicoquímicos, el recipiente de vidrio con tapa rosca para los análisis

microbiológicos y el Winkler para la determinación del oxígeno disuelto.

e) Para la preservación de las muestras se refrigeraron y se llevaron a una cava sellada. En

el caso de la prueba de Oxígeno disuelto es necesario preservarla con sulfato maganoso,

alcaliyoduro acida y ácido sulfúrico concentrado (un mililitro de cada reactivo).

f) Luego se enviaron al laboratorio de la empresa Servianalítica Profesionales S.A.S. en el

municipio de Ocaña, para su respectivo análisis.

3.8.6. Toma de muestras

Se realizó dos campañas de aforo. Los aforos realizados en los dos periodos de muestreo,

estuvieron en el contexto de variabilidad climática. Se tuvo en cuenta los datos históricos

mensuales de tiempo en Ábrego y Ocaña, en los que se describe que el periodo seco corresponde

a los meses de diciembre hasta marzo y los de lluvia de abril hasta el mes de noviembre.

Tabla 6
Histórico de tiempo de los municipios de Ábrego y Ocaña.

Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic

T. media (°C) 19,3 19,5 20,4 20,8 21,1 20,9 20,8 20,8 20,8 20,7 20,2 19,6
Ábrego

Pp (mm) 16 33 48 153 231 110 89 146 196 248 95 22

T. media (°C) 20,4 20,5 21,2 21,7 22 21,8 21,7 21,6 21,7 21,5 21 20,4
Ocaña

Pp (mm) 10 32 33 101 191 89 101 121 193 200 64 27

Nota: la tabla muestra los datos históricos de temperatura media y de precipitación en los
municipios de Ábrego y Ocaña. El color amarillo muestra la época de menos lluvias y el color
azul la época de más lluvias. Fuente: 1 (CLIMATE-DATA.ORG, 2018)
 48

3.9. Metodología de la modelación

La calibración y validación del modelo de calidad de agua, dependeré de los datos

adicionales que requiere el modelo como la posición geográfica de las estaciones, donde se

indiquen los parámetros de cada estación teniendo como base los puntos de monitoreo inicial y

final (Pelletier & Chapra, 2008).

 Ejecución del modelo de calidad del agua seleccionado.

 Análisis de los resultados obtenidos

3.10. Instrumentos utilizados para la recolección de la información.

En las Tabla 7 se discriminan los instrumentos utilizados según las actividades

programadas en cada objetivo específico.

Tabla 7
Descripción de equipos y herramientas utilizadas.

Actividad Equipos y herramientas

 Revisión documental.
 Revisión de la información de los permisos de vertimientos.
Información preliminar
 Reconocimiento in situ de los puntos de vertimientos más
significativos y las captaciones.
 Revisión de bases de datos.
 Revisión del estado del arte.
Selección del modelo de
 Información disponible.
simulación
 Características geomorfológicas del río algodonal.
 Características fisicoquímicas de río algodonal.
 GPS.
Georreferenciación  Imágenes satelitales.
 Google Earth Pro.
Selección de las estaciones de  Protocolo para el Monitoreo y Seguimiento del Agua.
monitoreo  Facilidad de acceso.
 49

“Tabla 7 continuación”

Actividad Equipos y herramientas

 Cartera de aforo.
Aforo  Planillero.
 Personal: dos técnicos (inspector y aforador).
 Correntómetro con molinete.
 Para medir el nivel del agua se utilizará una estaca afora cada
10 cm o jalones de colores vivos (rojo, blanco).
 Para la medición de las velocidades se empleará un
correntómetro.
 Cuerda.
 Cinta métrica.
 Linterna.
Obtención de los tiempos de  La razón entre la longitud del tramo 𝐿 y la velocidad media
viaje del agua 𝑣̅ .
 Datos del aforo.
Muestreo de los parámetros
 Marcador permanente.
fisicoquímicos
 Etiquetas.
 Cinta adhesiva.
 Cinta de papel.
 Termómetro ambiental y para medir la temperatura del agua.
 6 recipientes de plástico, tapa rosca de capacidad de 2 litros.
 6 recipientes Winkler.
 6 recipientes de vidrio para análisis microbiológico.
 Papel kraft o papel aluminio.
 Elementos de protección personal.
 GPS
 Cámara digital.
 Decámetro.
 Software ArcGIS.

Fuente: autores del proyecto.

 50

3.11. Procesamiento y análisis de la información

3.11.1. Procesamiento de la información climatológica

Se graficaron los datos de precipitación vs evaporación total media, y la temperatura

media anual vs la humedad relativa de cada una de las estaciones seleccionadas. Sin embargo,

los registros meteorológicos proporcionados por el IDEAM están incompletos.

3.11.2. Morfometría del río

Todo el procesamiento de la información se hizo con el Software ArcGIS 10.5, de la

plataforma de ESRI utilizando la licencia de la Universidad.

Las características físicas de la cuenca tienen gran importancia sobre la respuesta

hidrológica; sin embargo, la escasez de información cartográfica y sus distintas escalas hace que

el estudio sea limitado y aumente la incertidumbre sobre la confiabilidad de los parámetros.

3.11.3. Verificación del modelo de simulación.

Se procede a la verificación del modelo a través de una lista de chequeo y por último, se

definen las variables meteorológicas a tener en consideración para el correcto funcionamiento del

modelo.

La información recolectada en esta etapa es la base para la selección de los parámetros

necesarios para medir la capacidad de autodepuración del río Algodonal.

3.11.4. Estaciones de monitoreo.

Las estaciones seleccionadas cumplirán con las condiciones mencionadas anteriormente,

del Protocolo para el Monitoreo y Seguimiento del Agua. Se elegirán 5 estaciones por criterio de

los autores entre la laguna de estabilización de Ábrego y la planta de agua potable de la Empresa
 51

de Servicios Públicos de Ocaña –ESPO-. Se tratará en lo posible que las distancias sean similares

o equidistantes.

Luego de seleccionadas las estaciones, se describirán sus características hidrológicas y

topográficas (latitud, longitud, velocidad, elevación, distancia, distancia acumulada, caudal). Y

se tomará los datos de temperatura ambiental, humedad relativa y temperatura del agua.

Para medir la longitud del tramo de cada estación se utilizó la herramienta Ruta de

Google Earth, y la velocidad media del agua se determinará, junto con el aforo del río, a través

de un correntómetro (molinete digital medidor de caudal flowatch, JDC).

3.11.5. Calibración del modelo de calidad del agua.

En el Q2k el color se utiliza para indicar si la información se va a introducir por el

usuario o la salida por el programa:

 Azul pálido designa valores de las variables y parámetros que han de ser

introducida por el usuario.

 De color amarillo pálido designa datos que el usuario entra. Estos datos

son mostrados en los gráficos generados por Q2K.

 Verde pálido designa valores de salida generados por Q2K.

 Colores sólidos oscuros se utilizan para las etiquetas y no debe ser

cambiado.

Segmentación e hidráulica: para la calibración del modelo se tiene en cuenta que las

principales ramas del sistema (es decir, el tallo principal y cada uno de los tributarios) se

denominan segmentos. Esta distinción tiene importancia práctica debido a que el software
 52

proporciona gráficos de los resultados del modelo en base a segmentos. Es decir, el software

genera parcelas individuales para el tallo principal, así como cada uno de los afluentes.

HW#1

HW#2 1
6 2
7 Tr
ib
8 1 3
2 12
ib 4
Tr 11 13
14 5
10
9 15 16
HW#3

Main stem
17
18
19
Tr
HW#4 ib
22 3 20
23
24 21
25 26

27
28
29

(a) A river with tributaries (b) Q2K reach representation

Figura 4. Esquema de segmentación del Q2k para un río con afluentes.

Cualquier modelo de alcance puede ser dividida en una serie de elementos igualmente

espaciados. Como en Figura 5, Esto se hace simplemente especificando el número de elementos

que se deseen.

n=4

Reach Elements

Figura 5. División del tramo en segmentos.

Si se desea, cualquier segmento puede dividirse en una serie de n elementos de igual longitud.

En resumen, la nomenclatura utilizada para describir la forma en que organiza Q2K

topología río es el siguiente:

 53

 Reach (alcance). Una longitud de río con características hidráulicas

constantes.

 Element (elemento). unidad de cálculo fundamental del modelo que consta

de una longitud subdivisión igual de un alcance.

 Segment (segmento). Una colección de tramos que representan una rama

del sistema. Estos consisten en el tallo principal, así como cada tributario.

 Headwater (cabecera). El límite superior de un segmento de modelo.

Balance de flujo: como se describe en la sección anterior, la unidad más fundamental de

Q2K es el elemento. Un balance de flujo en estado estacionario se implementa para cada

elemento del modelo como:

𝑄𝑖 = 𝑄𝑖−1 + 𝑄𝑖𝑛,𝑖 − 𝑄𝑜𝑢𝑡,𝑖

Donde:

Qi = flujo de salida desde el elemento i en el elemento aguas abajo i + 1 [m3 / d]

Qi-1 = flujo de entrada desde el elemento corriente arriba i - 1 [m3 / d]

Qin, i es el flujo de entrada total en el elemento de punto y fuentes no puntuales [m3 / d]

Qout, i es el flujo de salida total desde el elemento debido a punto y los retiros no

puntuales [m3 / d]

Por lo tanto, la salida de aguas abajo es simplemente la diferencia entre la entrada y la

fuente ganancias menos las pérdidas de abstinencia.

 54

Qin,i Qout,i

Qi1 Qi
i1 i i+1

Figura 6.Flujo en equilibrio en el elemento.

3.11.6. Características hidráulicas

Una vez que se calcula el flujo de salida para cada elemento, la profundidad y la

velocidad se calculan en una de tres maneras: vertederos, curvas de calificación, y ecuaciones

Manning. El programa decide entre estas opciones de la siguiente manera:

 Si se introducen altura y el ancho del vertedero, se implementa la opción

vertedero.

 Si la altura del vertedero y la anchura son cero y se introducen coeficientes

de la curva de clasificación. Se implementa la opción curva de gasto.

 Si no se cumple ninguna de las condiciones anteriores, Q2K utiliza la

ecuación de Manning.

Tabla 8
Coeficiente de rugosidad Manning para canales de corriente natural (Chow et al. 1988).

MATERIAL n
Canales de corriente Natural
0.03-0.05
Limpia, bobinado y algunas malas hierbas
0.05
Las malas hierbas y piscinas, sinuosas
Arroyos de montaña con los cantos 0.04-0.10
rodados
 55

3.11.7 Escenarios prospectivos de calidad del cauce principal del río Algodonal.

Una de las utilidades de los modelos de simulación es la capacidad predictiva, ante los

escenarios que establezca el autor. Sin embargo, es necesario que el modelo esté calibrado con

las condiciones naturales de la corriente principal donde influyen los vertimientos.

Los valores de calidad a modelar (Tabla 9) en los escenarios son valores de aguas

residuales domesticas sin tratamiento establecidos por la literatura en función del

comportamiento del río Algodonal en tiempo seco (campaña 1).

Tabla 9
Valores típicos de las aguas residuales domésticas.

Concentración
Compuesto
Mínima Promedio Máxima
SST mg/L 100 220 350
DBO5 mgO2/L 110 220 400
Nitrógeno orgánico
8 15 35
mg/L
Nitrógeno
12 25 50
amoniacal mg/L
Alcalinidad mg
510 100 200
CaCO3/L
Nota: Adaptado de Metcalf y Eddy 1991 (como se citó en (Jiménez, 2001).

En el modelo se ingresó los valores de SSI como el 30% de los SST. Por otro lado, la DBO5

se tomó como DBOf. El caudal vertido utilizado en la simulación corresponde a 0,1 m3/s, el

oxígeno disuelto 4 mgO2/L y un valor de pH de 6.

Capítulo 4 Resultados y discusión

4.1. Características naturales del río Algodonal.

4.1.1. Climatología.
 56

Se solicitó información del catálogo nacional de estaciones del IDEAM sobre

precipitación, temperatura, humedad relativa y evaporación de los últimos 20 años, recolectadas

en las estaciones activas más cercanas al área de estudio. En la figura 7 se muestra el mapa de las

estaciones climatológicas utilizadas para la caracterización climática en el límite de la Cuenca

Hidrológica del Río Algodonal.

Figura 7. Localización de las estaciones climatológicas.

 57

Los datos proporcionados por el IDEAM presentan escasez de datos en algunos meses

por lo que el promedio anual se ve afectado, provocando picos inconsistentes en las gráficas.

Valores totales anuales de Precipitación y Evaporación

Estación: UNIV FCO P SANTAND

Precipitación Total Anual (mm) Evaporación

1600
1400
1200
1000
800
600
400
200
0
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017
Años

Figura 8. Relación de la precipitación y evaporación total anual (mm) de la estación Univ Fco P
Santand.

Temperatura media anual y Humedad Relativa anual (IDEAM)

Estación: UNIV FCO P SANTAND

Temperatura Media Anual Humedad Relativa (%)

23 100
Temperatura media (°C)

Humedad Relativa (%)

22.5 80

22 60

21.5 40

21 20

20.5 0
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017
2018

Años

Figura 9. Relación de temperatura y humedad de la estación Univ Fco P Santand.

 58

En este periodo (1997-2017) se presentaron lluvias de más de 600 mm (600 litros de agua

en un m2) y una precipitación máxima de 1462 mm en el año 2010. La temperatura se mantuvo

en un rango de 21 a 22,5°C.

Valores totales anuales de Precipitación y Evaporación

Estación: ABREGO CENTRO ADMO

Precipitación Total Anual (mm) Evaporación Total Anual (mm)

1600.0
1400.0
1200.0
1000.0
800.0
600.0
400.0
200.0
0.0
2001

2012
1997
1998
1999
2000

2002
2003
2004
2005
2006
2007
2008
2009
2010
2011

2013
2014
2015
2016
2017
Años

Figura 10. Relación de la precipitación y evaporación de la estación Abrego Centro Admo.

Temperatura media anual y Humedad Relativa anual (IDEAM)

Estación: ABREGO CENTRO ADMO

Temperatura Media Anual Humedad Relativa (%)

23 100
Temperatura media (°C)

Humedad Relativa (%)

22.5 80

22 60

21.5 40

21 20

20.5 0
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017

Años

Figura 11. Relación de temperatura y humedad de la estación Abrego Centro Admo.

 59

En esta estación se presentaron lluvias de más de 800 mm (800 litros de agua en un m2) y

una precipitación máxima de 1571 mm y 1592 mm en los años 2010-2011 respectivamente. La

temperatura durante 15 años (1999-2014) se mantuvo en un rango de 21,5 a 22,5°C.

Valores totales anuales de Precipitación y Evaporación

Estación: PLAYA LA

Precipitación Total Anual (mm) Evaporación Total Anual (mm)

1600
1400
1200
1000
800
600
400
200
0
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017
Años

Figura 12. Relación de la precipitación y evaporación de la estación Playa La.

Temperatura media anual y Humedad Relativa anual (IDEAM)

Estación: PLAYA LA

Temperatura Media Anual Humedad Relativa (%)

21.5 100
Temperatura media (°C)

Humedad Relativa (%)

21 80

20.5 60

20 40

19.5 20

19 0
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017

Años

Figura 13. Relación de temperatura y humedad de la estación Playa La.

Nota: Los datos representados en las figuras 5 hasta la 10, que corresponden a los años en los que el
IDEAM no contaban con registros completos fueron bordeados con color negro por los autores del proyecto. De
 60

igual manera el color degradado de la figura de temperatura, relaciona de color azul los datos de temperatura media
cercana o menores a los 21°C y los de color rojo a los datos cercanos a los 22°C.

Dentro del horizonte temporal de esta estación se presentaron lluvias de más de 400 mm

(400 litros de agua en un m2) y una precipitación máxima de 1425 mm en el año 2010. La

temperatura se mantuvo bajo los 20,5°C hasta el año 2010, en 2011 bajó de nuevo, pero siguió

subiendo hasta el año 2017 en 21,3°C.

Las características climatológicas de las estaciones correspondientes a Bosque Seco

Premontano bs-PM según (Holdridge 1967).

4.1.2. Características morfométricas de la cuenca

Como bien hemos expuesto, la autodepuración se está viendo amenazada por la presión

ecosistémica, la contaminación, cambio climático… pero también la disponibilidad espacial y

temporal de los recursos hídricos están siendo afectados significativamente. Por esto necesitamos

una descripción de las características de nuestras cuencas hidrológicas, que permita

caracterizarlas y que sirva de base para planificar el recurso en cuanto a su uso y disponibilidad

en proyectos futuros. (Ver figura 10)

4.1.2.1. Parámetros relacionados a la forma de la cuenca.

La forma de la cuenca afecta la velocidad del escurrimiento superficial de una corriente.

La metodología para determinar la forma de la cuenca permite comparar la forma de la cuenca

con figuras geométricas conocidas y permite comparar los resultados adimensionales con otras

cuencas en las que se pueda tener mayor información histórica de su comportamiento

hidrológico (Lux Cardona, 2013).

 61

Tabla 10
Parámetros morfométricos de la cuenca del río Algodonal.

Parámetros Valor Interpretación

Área de la cuenca (A) 499,99 km2 Cuenca
Perímetro de la cuenca (P) 134,17 km
Longitud axial o longitud máxima de la
42,46 km
cuenca (Lm)
Ancho máximo de la cuenca 17,05 km
Ancho de la cuenca (B) 11,77 km
Ancho medio de la cuenca 16,11 km
Orientación
Factor de forma de Horton (Kf) 0,28 Ligeramente achatada
Coeficiente de compacidad o índice de
1,68 Cuenca a rectangular oblonga
Gravelius (Kc)
Índice de alargamiento (Ia) 2,49 Moderadamente alargada
Cuenca con proporciones
Relación de elongación (Re) 0,59
accidentales
Fuente: autores del proyecto.

El tamaño relativo del área de la cuenca define, aunque no de manera rígida, los nombres

de sector, microcuenca, subcuenca o cuenca. El perímetro de la cuenca corresponde a la longitud

del parte aguas y la longitud axial es la distancia que existe entre la desembocadura de la cuenca y el

punto más alejado de la cuenca en una recta paralela al cauce principal. El ancho máximo de la

cuenca fue medido de extremo a extremo en una línea recta horizontal, el ancho de la cuenca es la

relación entre el área y la longitud de la cuenca y el ancho medio de la cuenca es la relación entre el

área y la longitud del cauce principal.

 62

Figura 14. Parámetros morfométricos de la cuenca.

 63

Figura 15. Ancho máximo de la cuenca.

 64

Figura 16. El mapa representa la orientación de la cuenca a través de colores.

Fuente: autores del proyecto.

 65

Para su generación se utilizó como base las curvas de nivel cada 10 metros. Fuente: autores del
proyecto.

El factor forma de Horton intenta medir cuan cuadrada (alargada) puede ser la cuenca. En

este caso, el factor forma es bajo lo que significa que tiene menos probabilidad de crecientes ya que

su forma ligeramente achatada no tiende a concentrar el escurrimiento de las lluvias intensas. E l

coeficiente de compacidad compara la forma de la cuenca con una circunferencia, por lo que un

valor de 1 corresponde a ser redonda. Su forma oblonga aleja la tendencia de concentrar

volúmenes de agua de escurrimiento (mediana susceptibilidad a la torrencialidad). Así mismo, el

índice de alargamiento muestra el comportamiento de la forma de la cuenca respecto a su

tendencia a ser alargada; este define que la cuenca es moderadamente alargada.

Finalmente la Relación de elongación, según estudios realizados (Summerfield, 1991) si

𝑅𝑒 ≈ 1.0 la cuenca es plana, y si ≥ 0,6 𝑅𝑒 ≤ 0,8 la cuenca es de relieve pronunciado.

4.1.2.2. Parámetros relacionados al relieve de la cuenca.

Permite relacionar la oferta hídrica y el movimiento del agua a lo largo de la cuenca. De

la altura y la elevación depende en gran medida la cobertura vegetal, el clima, la biota y otras

características del relieve en el territorio.

La pendiente media (Alvord y Horton) corresponde al valor medio del declive del terreno

y la inclinación, respecto a la horizontal, de la vertiente sobre la cual se ubica la cuenca.

∑ 𝑃𝑟𝑜𝑚𝑒𝑑𝑖𝑜 𝑝𝑜𝑟 𝑜𝑐𝑢𝑟𝑟𝑒𝑛𝑐𝑖𝑎

𝑆𝑚 =
∑ 𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑜𝑐𝑢𝑟𝑟𝑒𝑛𝑐𝑖𝑎
 66

Tabla 11
Pendiente media de la cuenca.

CÁLCULO DE LA PENDIENTE MEDIA DE LA CUENCA

INTERVALO DE PENDIENTE (%)
Pendiente NÚMERO DE
N x Sm %
Nº media del OCURRENCIA
INFERIOR SUPERIOR (%) Pendientes
intervalo (N)
(Sm)
1 0 5 2,5 1627 4067,5 9,21604169
2 5 12 8,5 3642 30957 20,6298856
3 12 18 15,0 2623 39345 14,8578226
4 18 24 21,0 2113 44373 11,9689589
5 24 32 28,0 2303 64484 13,0452022
6 32 44 38,0 2417 91846 13,6909482
7 44 100 72,0 2929 210888 16,5911408
17654 485960,5
Fuente: autores del proyecto.

Tabla 12
Parámetros relacionados con el relieve de la cuenca.

Parámetros Valor Interpretación

3477,66
Cota mayor de la cuenca (CM)
m.s.n.m.
1206,24
Cota menor de la cuenca (Cm)
m.s.n.m.
Altitud o elevación media de la cuenca 2341,95
(Hm) m.s.n.m.
Fuertemente accidentado
Pendiente media de la cuenca (S) 27,53%
(P5)
Coeficiente de masividad (CM) 4,68 Moderadamente montañosa
Coeficiente orográfico (CO) 10969,65 Relieve accidentado
Fuente: autores del proyecto.

La cuenca tiene un relieve fuertemente accidentado donde se pueden presentar problemas

de drenaje y sedimentación.
 67

Figura 17. Mapa de pendientes (%) de la cuenca.

 68

El coeficiente de masividad representa valores bajos (>35) en cuencas planas o

moderadamente montañosa. El coeficiente orográfico expresa el potencial de degradación de la

cuenca, el valor es > 6 representa un relieve accidentado.

Tabla 13
Áreas entre curvas de nivel.

Porcentaje
Nº Intervalo entre Elevación Porcentaje
Área (km2) de área
ORDEN curvas de nivel media acumulado
entre C.N.
8 2750 - 3477 2963 12,64 2,53% 2,53%
7 2500 - 2750 2625 11,71 2,34% 4,87%
6 2250 - 2500 2375 19,08 3,82% 8,69%
5 2000 - 2250 2125 34,72 6,94% 15,63%
4 1750 - 2000 1875 82,30 16,46% 32,09%
3 1500 - 1750 1625 164,83 32,97% 65,06%
2 1250 - 1500 1375 171,38 34,28% 99,33%
1 < 1250 1125 1,61 0,32% 99,66%
498,27 99,66%
Fuente: autores del proyecto.

La curva hipsométrica representa de manera gráfica la variación altitudinal de una cuenca, en

la que a cada altura le corresponde un respectivo porcentaje del área ubicada por encima de esa

altura. Las curvas superior refleja una cuenca con un gran potencial erosivo; la curva

intermedia es característica de una cuenca en equilibrio (Gaspari, y otros, 2013).

Figura 18. Análisis hipsométrico de una cuenca. Fuente: (Campos Aranda, 1998)
 69

Figura 19. Mapa de elevaciones cada 100 metros.

 70

Los datos de elevación son significativos, ya que la altitud afecta el comportamiento de la

temperatura y la precipitación.

Curva Hipsométrica

3050

2800

2550
Elevación media

2300

2050

1800

1550

1300

1050
0.00% 10.00% 20.00% 30.00% 40.00% 50.00% 60.00% 70.00% 80.00% 90.00% 100.00%
Porcentaje de área sobre la elevación media

Figura 20. Curva hipsométrica.

Describe a relación entre la propiedad altimétrica de la cuenca en un plano y su elevación.

Fuente: autores del proyecto. Nota: La cuenca está en fase de madurez (cuenca en equilibrio).

Histograma de frecuencias altimétricas

40.00%

30.00%
% de área

20.00%

10.00%

0.00%
0 - 1250 1250 - 1500 1500 - 1750 1750 - 2000 2000 - 2250 2250 - 2500 2500 - 2750 2750 - 3477
Intervalo de Alturas (m.s.n.m.)

Figura 21. Histograma de frecuencias altimétricas.

 71

Es la representación de la superficie, comprendida entre dos cotas, siendo la marca de clase el

promedio de las alturas. Realmente la curva hipsométrica y el histograma contienen la misma
información pero con una representación diferente (Ibáñez Asensio, Moreno Ramón, & Gisbert
Blanquer, 2008). Fuente: Autores.

4.1.2.3. Parámetros relacionados al perfil de la cuenca.

Parámetros Valor Interpretación

Longitud del cauce principal (Lc) 31,05 km Largo
Cota mayor del cauce principal (CMc) 1430,00 m.s.n.m.
Cota menor del cauce principal (Cmc) 1210,00 m.s.n.m.
Longitud total de cauces (Lt) 322,76 km

4.1.2.4. Parámetros relacionados al drenaje de la cuenca.

Parámetros Valor Interpretación

Orden de los cauces (Horton, 1945) 4,00
Longitud de los cauces de orden uno (L1) 154,93 km
Densidad del drenaje Horton (Dd) 0,65 Baja (cuenca mal drenada)
Constante de estabilidad del río (C): 1,55 Alto
Coeficiente de torrencialidad (CT) 0,15 Bajo
Densidad de corriente (Dc) o Densidad
0,21 corrientes/Km2
hidrográfica (Dh):

La constante de estabilidad del río propuesta por Schumm (1956) es el valor inverso a la

densidad de drenaje y representa una medida de erodabilidad de la cuenca. Por lo que, una menor

densidad de drenaje y altos valores de la constante de estabilidad es propia de materiales del

terreno resistentes a la erosión, suelos muy permeables, y pendientes suaves. Por lo que implica

una elevada capacidad de infiltración.

El Coeficiente de Torrencialidad resulta de la relación entre el número de cauces de

agua de orden uno y el área de la cuenca. Si el coeficiente es bajo significa que son áreas que se
 72

activan como respuesta a lluvias de larga duración cubriendo zonas extensas, sin ocasionar

procesos erosivos ni crecientes.

La densidad de corriente o hidrográfica corresponde al número de corrientes sobre el

área de la cuenca.

Figura 22. Orden de los cauces (Horton, 1945).

 73

4.2. Selección del modelo de simulación

Como se mencionó anteriormente, los ríos funcionan como reactores biogeoquímicos que

absorben o asimilan las distintas cargas contaminantes que este puede recibir en su recorrido

natural y, al mismo tiempo recupera las propiedades naturales del agua. En general, “cuanto más

caudalosos, rápidos y más aislados estén los cursos de agua y no hayan sido utilizados, más

capaces son de tolerar una cantidad mayor de aguas residuales”; sin embargo, “una cantidad

excesiva de cualquier tipo de contaminante produce daños a todo el sistema” (Suárez, 2008, pág.

2).

Estos contaminantes presentan un comportamiento dinámico, ya que sufren diferentes

fenómenos. Unos fenómenos físicos que se refieren a los elementos que se sedimentan en el

fondo del río por la diferencia de densidades, otros que son arrastrados por la corriente que

permanecen en suspensión y otros de baja densidad que se mantienen en la superficie como las

grasas y los detergentes. También podemos encontrar fenómenos de degradación biológica, que

son los que realizan los microorganismos y fenómenos químicos.

El trabajo de la modelación consiste en describir matemáticamente el comportamiento o

las interacciones que se producen en los sistemas acuáticos, en la que se relaciona las entradas de

contaminación y la calidad del agua.

Usualmente, la elección del modelo a implementar depende de la cantidad de información

con que se cuenta, la intención de la modelación y los procesos a modelar, esto implica que

existe una restricción, que puede ser ajustada en la calibración de los modelos. Los modelos

hidrodinámicos y de calidad de agua integran las entradas de fuentes puntuales y no puntuales

para determinar los impactos en la calidad del agua en la recepción de cuerpos de agua. Estos

modelos permiten a los científicos e ingenieros determinar las capacidades de asimilación del
 74

cuerpo de agua, determinar el nivel de mejores prácticas de gestión o predecir el tiempo

requerido para que un sistema se recupere después de ser alterado (United States Environmental

Protection Agency (EPA), 2015), ejemplo de ello son:

 Programa de simulación de análisis de calidad del agua “Water Quality Analysis

Simulation Program” (WASP)

 Código ambiental de dinámica de fluidos “Environmental Fluid Dynamics Code”

(EFDC)

 Modelo de calidad de agua corriente “Water Quality Model” (QUAL2K)

 Modelo hidrodinámico y de calidad de agua dinámico unidimensional “A Dynamic One

Dimensional Hydrodynamic & Water Quality Model” (EPD-RIV1)

De esta manera, el Modelo de Calidad de Agua Corriente “Water Quality Model”

(QUAL2K), fue diseñado por Greg Pelletier, Steve Chapra, y Hua Tao del Departamento de

ecología, Olympia, Washington. Tal como está configurado actualmente, un libro de Excel sirve

como interfaz para QUAL2K, es decir, todas las entradas y salida, así como la ejecución del

modelo se implementan dentro de Excel. Todas las funciones de interfaz están programados en

lenguaje de macros de Excel: Visual Basic para Aplicaciones (VBA). Todos los cálculos

numéricos se implementan en Fortran 90 para la velocidad de ejecución.

La selección del modelo se efectuó considerando aspectos como:

 Definición del problema a estudiar

 Características generales del modelo: tipos de procesos, tipos de métodos de

solución, tipo de cuerpo de agua, dimensión.

 75

 Capacidad de simulación de los parámetros básicos de calidad: oxígeno disuelto,

demanda biológica de oxígeno. Además, tiene en cuenta los contaminantes de los

vertimientos.

 La calibración, permite adaptarse a las condiciones de diferentes corrientes.

 Antecedentes de aplicabilidad: se ha demostrado el éxito en otras regiones.

 Viabilidad económica: es un modelo de calidad de agua gratuito.

Resumen de los datos de entrada al modelo y chequeo realizado

Tabla 14
Datos de entrada al modelo.

Información
Hojas de trabajo Descripción existente o asumida
por la literatura
Información general:
-Nombre del río
-Zona horaria
Qual2 k -Duración de la simulación Si
-Método de solución
Etc.

Información de cabecera y
condiciones de contorno aguas
Hoja de cabecera (Headwater) Si
abajo.

Información relacionada con la

cabecera del río y sus alcances:
-Distancias
Alcance (Reach) Si
-Elevaciones
Etc.

Condiciones iniciales en la
Condiciones iniciales (Initial
columna de agua en cada tramo. Si
Conditions)
Air Temperature, Dew Point
Variables meteorológicas. Si
Temperature, Wind Speed, Shade
Parámetros de velocidad del
Rates modelo y opciones para la Si
calibración automática.
 76

“Tabla 14 continuación”

Información
Hojas de trabajo Descripción existente o asumida
por la literatura
Información de ubicación, caudal
Hydraulics Data, Temperature y calidad, medida en las
Si
Data, WQ Data estaciones de monitoreo sobre el
río.

Definición de las variables meteorológicas

Las variables meteorológicas son precipitación (mm), temperatura del aire (°C), humedad

relativa (%), velocidad del viento (m/s), temperatura de rocío (°C), sombra (%), nubosidad (%).

Se tomó como referencia los datos tomados en campo y con información disponible en la

web.

Tabla 15
Características climatológicas de las estaciones en las horas del muestreo en la campaña 1.

Velocidad
Temperatura Punto de
Hora Humedad % del viento Nubosidad
aire °C rocío
m/s
12:00 a. m. 22.8 78 18.9 4.02 40%
01:00 a. m. 22.2 83 18.9 2.68 10%
02:00 a. m. 21.1 88 18.9 3.13 10%
03:00 a. m. 20.0 88 17.8 2.24 10%
04:00 a. m. 21.1 88 18.9 0.89 10%
05:00 a. m. 20.0 94 18.9 0.00 40%
06:00 a. m. 20.0 94 18.9 1.34 65%
07:00 a. m. 20.0 94 18.9 0.89 65%
08:00 a. m. 22.8 83 20.0 2.24 40%
 77

“Tabla 15 continuación”

Velocidad
Temperatura Punto de
Hora Humedad % del viento Nubosidad
aire °C rocío
m/s
09:00 a. m. 26.1 69 20.0 2.68 65%
10:00 a. m. 27.2 61 18.9 2.24 40%
11:00 a. m. 28.9 58 20.0 2.24 40%
12:00 p. m. 31.1 49 18.9 2.24 40%
01:00 p. m. 31.1 49 18.9 3.58 40%
02:00 p. m. 32.8 41 17.8 6.26 40%
03:00 p. m. 32.2 40 17.2 6.26 65%
04:00 p. m. 32.2 38 16.1 4.02 65%
05:00 p. m. 30.0 48 17.8 8.05 40%
06:00 p. m. 27.8 54 17.8 8.05 10%
07:00 p. m. 26.1 61 17.8 5.36 40%
08:00 p. m. 25.0 65 17.8 3.58 10%
09:00 p. m. 22.8 73 17.8 3.13 10%
10:00 p. m. 22.8 73 17.8 3.58 10%
11:00 p. m. 22.8 73 17.8 3.13 10%

Tabla 16
Características climatológicas de las estaciones en las horas del muestreo en la campaña 2.

Velocidad
Temperatura Punto de
Hora Humedad % del viento Nubosidad
aire °C rocío °C
m/s
12:00 a. m. 28.9 58 20.0 2.68 40%
01:00 a. m. 21.1 94 20.0 2.68 10%
02:00 a. m. 21.1 94 20.0 2.24 10%
03:00 a. m. 22.2 88 20.0 1.34 40%
 78

“Tabla 16 continuación”

Velocidad
Temperatura Punto de
Hora Humedad % del viento Nubosidad
aire °C rocío °C
m/s
04:00 a. m. 21.1 94 20.0 1.34 65%
05:00 a. m. 20.0 94 18.9 1.34 40%
06:00 a. m. 20.0 94 18.9 1.34 65%
07:00 a. m. 21.1 49 10.0 2.24 10%
08:00 a. m. 22.8 88 21.1 0.89 10%
09:00 a. m. 25.0 83 22.2 4.02 40%
10:00 a. m. 27.8 62 20.0 1.34 40%
11:00 a. m. 30.0 58 21.1 2.24 40%
12:00 p. m. 21.1 94 20.0 2.68 10%
01:00 p. m. 30.0 55 20.0 2.24 40%
02:00 p. m. 31.1 70 25.0 2.24 65%
03:00 p. m. 30.0 55 20.0 2.24 65%
04:00 p. m. 28.9 58 20.0 6.26 65%
05:00 p. m. 27.2 65 20.0 4.47 65%
06:00 p. m. 26.1 69 25.6 4.02 10%
07:00 p. m. 25.0 69 18.9 3.13 65%
08:00 p. m. 22.8 94 22.2 2.24 10%
09:00 p. m. 22.8 94 22.2 2.68 10%
10:00 p. m. 22.2 100 22.2 2.68 65%
11:00 p. m. 22.2 100 22.2 1.34 65%

Nota: todos los datos fueron tomados de la página web (Weather Underground, 2018).
 79

4.3. Análisis de la calidad del agua

4.3.1. Vertimientos significativos que afectan la corriente principal

Durante el proceso de reconocimiento en campo de los vertimientos emitidos al río

Algodonal, se pudo determinar tres vertimientos más significativos, en función de su caudal

vertido. El primero es generado por la laguna de oxidación del municipio de Ábrego, la cual

recoge gran parte de las aguas servidas del municipio, ésta se encuentra en un pésimo estado

debido al poco o nulo mantenimiento por parte de la administración municipal, lo cual genera

que la carga contaminante de las aguas residuales domésticas no alcancen un porcentaje óptimo

de remoción y se generan olores ofensivos al medio (Alcaldía Ábrego, 2016-2019).

Así mismo, se identificó a cien metros aguas debajo de la laguna de oxidación, un

vertimiento indirecto generado por el barrio Alto Viejo, éste posee un sistema de alcantarillado

que no se encuentra conectado al canal de alimentación de la laguna de oxidación, sino que

vierte directamente a una quebrada, la cual posteriormente desemboca a unos 300 metros en el

río Algodonal. Cabe destacar que las aguas residuales generadas por este barrio no están siendo

tratadas, por lo que justo en la desembocadura del río se puede evidenciar una mancha oscura

que afecta la calidad paisajística.

Por último, se obtuvo información del sistema de alcantarillado de La Ermita identificada

en el Plan Básico de Ordenamiento Territorial (2002-2011) del municipio de Ocaña (numeral

5.18.5 de los servicios públicos básicos del suelo suburbano de La Ermita) donde el 97,2 % del

total de las viviendas no posee sistema de alcantarillado, y el 2,8 % de las viviendas poseen

pozos sépticos (Alcaldía Ocaña, 2002-2011). Luego, realizamos una visita al corregimiento de la

Ermita donde se verificó que La Ermita no posee sistema de alcantarillado y por tanto ningún

sistema de tratamiento, algunas casas tienen pozos sépticos, pero éstos ya perdieron su
 80

efectividad. Y las demás viviendas vierten directamente al río Algodonal. En conclusión, se

generan vertimientos difusos en toda la ribera del río. Es importante mencionar, que el

corregimiento de la Ermita se encuentra a menos de cinco kilómetros por río de la planta de

tratamiento de agua potable (PTAP) de ESPO S.A.S. y que ésta problemática podría afectar la

calidad del agua.

4.3.2. Captaciones

En la tabla 17 se describen las concesiones otorgadas por la corporación, y de color rojo

se señalan las que no están vigentes.

Tabla 17
Concesiones del río Algodonal

COORDE
SOLICITANTE

ECONOMICA

(QUEBRADA)

EXPEDICION

NADAS
CORRIENTE

CONCESION
OTORGADO
ACTIVIDAD
MUNICIPIO

DURACION
CAUDAL
VEREDA

FECHA

AÑOS
LPS

X Y

Diana
Rio
Carolina Quebrad 23/06/201 108357 139873
Ocaña Domestico Algodona 0,0255 5
Correa a Seca 7 3 0
l
Ríos
Consumo
Arisolina Humano,
Rio
Pacheco abrevadero El 0,4257 30/05/201 108366 139838
Ocaña Algodona 5
De Animal, Rincón 5 7 5 5
l
Arevalo riegocultiv
o
Abrevadero
Alexande de Rio
Abreg 23/06/201 109005 139120
r Rosado Animales, Oropoma Algodona 0.4103 5
o 7 2 8
Gomez Riego De l
Cultivo
Fuente: adoptado de los archivos de Corponor 2018, por los autores del proyecto.

4.3.3. Cargas contaminantes

 81

Para estimar la carga contaminante es necesario conocer el caudal de agua residual

vertido al río Algodonal. Para esto fue necesario completar la información de la tabla 5.

Suscriptores: se contaron las viviendas del barrio Alto Viejo a través de un sobrevuelo de

la zona con un drone DJI Phantom y se restó el dato a los suscriptores totales para obtener el dato

de las viviendas conectadas a la Planta de Tratamiento de Aguas Residuales PTAR.

Calculo población: para el municipio de Ábrego, se usó el dato de la población urbana

estimada por el DANE para el año 2015 (Alcaldía Ábrego, 2016-2019, pág. 37). Para la Ermita,

se multiplicó por cuatro la cantidad de suscriptores, para tener un dato aproximado.

Caudal de entrada: En la Ermita fue medido directamente en la PTAP de la Ermita

(Atores del proyecto).

Para calcular el caudal diario se utilizó la siguiente fórmula:

𝐿 𝐿 3600 𝑠 𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑑𝑎𝑑 (ℎ)

𝑄 = 𝑄𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ∗ ( )∗( )
𝑑í𝑎 𝑠 1 ℎ𝑜𝑟𝑎 1 𝑑í𝑎

𝐿 𝐿 3600 𝑠 16 ℎ)
Ábrego: 𝑄 = 56 ∗ ( ) ∗ (1 𝑑í𝑎 )
𝑑í𝑎 𝑠 1 ℎ𝑜𝑟𝑎

𝐿 𝐿 3600 𝑠 12 ℎ)
La Ermita: 𝑄 = 1,32 ∗ ( ) ∗ ( 1 𝑑í𝑎 )
𝑑í𝑎 𝑠 1 ℎ𝑜𝑟𝑎

Dotación por habitante o Caudal per cápita

𝐿
𝐿 𝑄 ⁄
𝑄 ∗ ℎ𝑎𝑏 = 𝑑í𝑎
𝑑í𝑎 𝑃𝑜𝑏𝑙𝑎𝑐𝑖ó𝑛

La dotación neta máxima por habitante según la altura sobre el nivel del mar de la zona

atendida es de 130 𝐿⁄ℎ𝑎𝑏 ∗ 𝑑í𝑎 para la Ábrego y la Ermita según la (Resolución N°330, 2017)
 82

(Por la cual se adopta el Reglamento Técnico para el Sector de Agua Potable y Saneamiento

Básico [RAS])

Caudal vertido. El vertimiento aproximado que genera el barrio Alto Viejo fue el

resultado de la multiplicación del Caudal L/día por la población del barrio. Por último, se

sumaron los dos caudales para un Caudal total en Ábrego de 48,53 L/s.

Por otro lado, se asume que el caudal de entrada de la Ermita en 12 horas, multiplicado

por el factor de retorno es el caudal total de salida, pero en diferentes puntos de vertimiento.

𝑚 3
𝐿 1 𝑑í𝑎 1𝑚3
𝑄𝑣 ⁄𝑠 = 𝑄 ⁄𝑑í𝑎 ∗ 𝐹𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑅𝑒𝑡𝑜𝑟𝑛𝑜 ∗ (86400 𝑠) (1000 𝐿)

Tabla 18
Cálculo del caudal vertido en m3/s en el municipio de Ábrego.

Vertimiento
PTAR Vertimiento del
Datos Casco urbano
Laguna de B. Alto Viejo
oxidación
Suscriptores 3860 3803 57
Población 17009 15212 228
Caudal de entrada a la PTAP L/s 56
Continuidad h/día 16
Caudal L/día 3225600,00
Dotación por habitante L/hab* día 189,64
Caudal vertido m3/s 0.04853 0.0353 0.01201
 83

Tabla 19
Cálculo del caudal vertido en m3/s en el corregimiento de La Ermita

Datos
Suscriptores 104
Población 416
Caudal de entrada a la PTAP L/s 1,32
Continuidad h/día 12
Caudal L/día 56842,11
Dotación por habitante L/hab* día 137,08
Factor retorno 0.8
Caudal vertido m3/s 0.000528

El valor de la carga contaminante kg/día para el agua de salida en la laguna de oxidación

se presenta en la siguiente tabla:

Tabla 20
Cálculo de carga contaminante del vertimiento de la laguna de oxidación en kg/día.

Datos obtenidos Carga contaminante diaria Cc kg/día

DBO5 (mg O2/L) SST (mg SST/L) DBO5 (mg O2/L) SST (mg SST/L)
14,8 21,2 45,14 64,66
Caudal promedio de AR (L/s)
35,3

Nota: cálculo realizado según la definición de carga contaminante del Decreto 2667 de 2012.

En el caso de la carga unitaria de origen doméstico para los vertimientos de la Ermita,

donde no existe sistema de alcantarillado, debe utilizarse los valores de la tabla E.2.6., del RAS

2000 título E.
 84

Tabla 21
Tabla E.2.6 (RAS 2000) Aportes per cápita para aguas residuales domésticas.

Valor sugerido
Parámetro Intervalo Valor sugerido
kg/hab*día
DBO5, 20°C, g/hab/día 25-80 50 0,05
Sólidos en suspensión, g/hab/día 30-100 50 0,05

4.3.4. Características de las estaciones de monitoreo

En la tabla 22 se describen la información geográfica de las estaciones de monitoreo

definidas a lo largo del río Algodonal.

Tabla 22
Información geográfica de las estaciones de muestreo.

Coordenadas WGS84
Distancia Elevación
Nombre Latitud Longitud
km m.s.n.m.
Grados Minutos Segundos Grados Minutos Segundos

Cabecera 0 1356 8 5 44.287 73 13 18.315

Estación 1 0.204 1355 8 5 44.092 73 13 23.824

Estación 2 6.549 1336 8 6 34.143 73 14 59.763

Estación 3 8.776 1306 8 8 40.640 73 16 55.890

Estación 4 8.776 1253 8 11 22.490 73 18 58.970

Estación 5 4.934 1205 8 13 11.150 73 19 14.000

Fuente: autores del proyecto.

 85

Figura 23. Localización de las estaciones de monitoreo.

 86

Tabla 23
Caracterización hidráulica de las estaciones.

Cota
Ancho de Cota final Longitud Pendiente Manning Descripción de la
Nombre inicial
fondo (m) (m.s.n.m.) (m) (m/m) (n) corriente
(m.s.n.m.)
Meandros con
Cabecera 12.35 1356 1355 204 0.00490 0.05
muchas piedras
Vertimiento Pequeños árboles y
0.5 1356 1355 203 0.00493 0.05
1 arbustos sin follaje
Vertimiento Escasos arbustos y
2.8 1367 1354 294 0.04422 0.05
2 pasto abundante
Estación Meandros con
23.5 1355 1336 6549 0.00290 0.05
1 muchas piedras
Estación Meandros con
12.25 1336 1306 8776 0.00342 0.05
2 muchas piedras
Estación Meandros con
11.7 1306 1253 8776 0.00604 0.05
3 muchas piedras
Vertimiento Meandros con
12 1249 1244 720 0.00694 0.05
3 muchas piedras
Estación Meandros con
13.6 1253 1205 4934 0.00973 0.05
4 muchas piedras
Estación Meandros con
33.25 1205 1198 1266 0.00553 0.05
5 muchas piedras

Nota: se tuvo en cuenta los coeficientes de Manning para cauces naturales. Fuente: autores del
proyecto.

En función de los sentidos, lo siguiente:

 87

Tabla 24
Características organolepticas de las estaciones de monitoreo.

Estación Acceso Estado del recurso Registro fotográfico

El agua en este punto estaba presentaba

Por la entrada de
condiciones de turbiedad porque se evidenciaban
la laguna de
actividades de dragado de material de arrastre
Cabecera oxidación del
aguas arriba de la estación; no se presentaron
municipio de
olores y el sustrato de la estación era de cantos y
Ábrego.
grava.

La estación uno E1 recibe los vertimientos de la

laguna de estabilización y un vertimiento del
municipio de Ábrego. En el primer periodo de
muestreo se evidenció un tercer vertimiento que
correspondía a una obstrucción en el canal de
traslado del agua residual hasta la laguna por
colmatación de sólidos, por lo que el agua residual
se rebosó y llegó al río Algodonal sin ningún
1
tratamiento. En este punto el agua presentaba mal
olor, con una velocidad de corriente lenta, se creó
una mancha en el lateral del río de color oscuro y
se estaba turbia por las actividades de dragado, se
encontró con presencia de peces, el sustrato
presentaba gran cantidad de sedimentos sobre los
cantos y la grava, por lo que aguas abajo seguían
dragando el material de río.
 88

“Tabla 24 continuación”

Estación Acceso Estado del recurso Registro fotográfico

En la estación dos E2 el agua seguía presentando

Por la estación de turbiedad, el sustrato era de cantos y grava, sin
2
servicio de Terpel. presencia de olores ni color, la velocidad del río
aumentó.

Por la trituradora En la estación tres E3 se encontró el agua incolora

3
Guayabal. e inodora, con un sustrato de gravas y cantos.
 89

“Tabla 24 continuación”

Estación Acceso Estado del recurso Registro fotográfico

Por la vía que de

En la cuarta estación E4 tampoco se presentó
la Ermita llega al
malos olores, pero con turbiedad. El sustrato del
4 conjunto
río en ese punto está compuesto de una capa
residencial
gruesa de sedimentos.
Algodonal Club.

Por la Planta de En la estación cinco E5, no presenta olores ni

Tratamiento de turbiedad, el sustrato es un canal artificial de
5 Agua Potable de concreto pues está diseñado para captar el agua en
ESPO en Ocaña. la planta de tratamiento de agua potable de la
empresa ESPO S.A.S.
 90

4.3.5. Tiempos de viaje

Las distancias entre las estaciones del primer periodo de muestreo (tiempo seco) y las del

segundo periodo (tiempo de lluvias) fueron medidas a través de Google Earth Pro.

Figura 24. Distancia entre las estaciones en el primer muestreo.

Figura 25. Distancias entre estaciones en el segundo muestreo.

Se modificó las distancias entre las estaciones del segundo muestreo en función de la facilidad en
la toma de muestras debido a la crecida del río. Fuente: Google Earth.
 91

El tiempo de viaje de la masa de agua, indica el momento aproximado en que la masa de

agua ha recorrido la distancia entre una estación y otra. En el caso del aforo y la toma de

muestras en la primera campaña de monitoreo, se realizó primero en la estación de Cabecera, por

lo que la distancia entre las estaciones y la velocidad para determinar el tiempo de viaje, fueron:

distancia entre la cabecera y estación 2, y la velocidad media de cabecera.

Tabla 25
Obtención de los tiempos de viaje en tiempo seco.

14 de
Fecha: 13 de marzo 13 de marzo 13 de marzo 13 de marzo 14 de marzo
marzo
Estación
Estación: Estación 1 Cabecera Estación 2 Estación 3 Estación 4
5
Distancia: Entre C y E1: Entre C y E2: Entre E2 y E3: Entre E3 y E4: Entre E4 y E5: No
(m) 204 6587 8940 8804 4909 aplica
Velocidad
media: 0,2587 0,4647 0,5000 0,3000 0,2667 0,3912
(m/s)
Tiempo de
No
viaje 0:13:09 3:56:15 4:58:00 8:09:07 5:06:49
aplica
(hh:mm:ss)

Por otro lado, el aforo y la toma de muestras de la segunda campaña inicio con la

Estación 1, luego cabecera y así sucesivamente las demás estaciones río abajo, dependiendo de la

hora del tiempo de viaje.

Tabla 26
Obtención de los tiempos de viaje en tiempo de lluvia.

20 de
Fecha: 19 de abril 19 de abril 19 de abril 19 de abril 19 de abril
abril
Estación
Estación: Cabecera Estación 1 Estación 2 Estación 3 Estación 4
5
Distancia: Entre C y E1: Entre E1 y E2: Entre E2 y E3: Entre E3 y E4: Entre E4 y E5:
No aplica
(m) 204 6469 9251 8519 4603
Velocidad
media: 0,74737 0,41818 0,52000 0,34000 0,41818 0,6212
(m/s)
Tiempo de
viaje 0:04:33 4:17:49 4:56:30 6:57:36 3:03:27 No aplica
(hh:mm:ss)
 92

4.3.6. Aforo

La primera campaña de aforo corresponde a periodo seco (mes de marzo) y la segunda

corresponde al periodo de lluvias (mes de abril). En la tabla 27 se muestran los resultados de la

campaña de aforo.

Tabla 27
Datos hidráulicos en las estaciones de monitoreo.

13-14 de marzo 19-20 abril

Incremento

Incremento
Distancia

%
Nombre acumulada Q H U Q H U
(km) (m3/s) (m) (m/s) (m3/s) (m) (m/s)

Cabecera 0 1.35 0.21 0.46 5.04 0.32 0.75 3.69 273.33%

Estación 1 0.204 1.24 0.2 0.26 4.47 0.46 0.42 3.23 260.48%
Estación 2 6.753 1.39 0.26 0.5 4.85 0.43 0.52 3.46 248.92%
Estación 3 15.529 1.19 0.3 0.3 3.62 0.66 0.34 2.43 204.20%
Estación 4 24.305 2.41 0.62 0.27 5.46 0.80 0.42 3.05 126.39%
Estación 5 29.239 1.55 1.16 0.39 4.79 0.27 0.62 3.24 209.03%
Nota: las variables "H" y "U" hacen referencia a la profundidad y a la velocidad del cauce,
respectivamente. En la tabla se presenta la velocidad media en cada estación monitoreada y el
valor del caudal. Además, se presenta el incremento y la media aritmética (promedio) del caudal.
Fuente: 2. Autores del proyecto.

4.3.7. Resultados de los análisis de calidad del agua

Tabla 28
Hora de la toma de muestras en tiempo seco.

13 de 13 de 13 de 14 de 14 de
Fecha: 13 de marzo
marzo marzo marzo marzo marzo
Estación: Estación 1 Cabecera Estación 2 Estación 3 Estación 4 Estación 5
Hora de toma 9:20:00 a. 1:16:15 p. 6:14:15 p. 2:23:21 a. 7:30:10 a.
9:40:00 a. m.
de muestra: m. m. m. m. m.
 93

Tabla 29
Datos de calidad del agua en las estaciones de monitoreo en tiempo seco.

Parámetro Unidades Cabecera E1 E2 E3 E4 E5

Temperatura (T) °C 21,3 21,7 22,9 24,7 23,1 22,1

Conductividad μohms 62,5 78,2 70,9 75,4 78 83,9

Sólidos
mgD/L 12 17 11 21 14 13
Inorgánicos (SSI)
Oxígeno disuelto mgO2/L 8,1 8 8,8 7,4 8,2 7,5
DBO lenta mgO2/L 1,3 1,1 2,3 2,4 2,9 3,1
DBO rápida mgO2/L 0,9 3,9 1,7 1,5 1,5 1,4
NH4- nitrógeno μgN/L 60 960 0 90 500 100
NO3- nitrógeno μgN/L 30 50 90 120 210 120
Fósforo inorgánico
μgP/L 70 320 180 150 30 190
(SRP)
Detritos o Material
Orgánico mgD/L 0,1 0,3 0,4 0,3 0,1 0,2
Particulado MOP
Patógenos Ufc/100mL 4000 124000 10000 8000 92000 38000

Alcalinidad mgCaCO3/L 35 38 40 40 38 38
pH s.u. 7,11 6,81 7,11 7,13 7,16 7,1
Caudal m3/s 1,35 1,24 1,39 1,19 2,41 1,55
Nota: datos obtenidos en la campaña de muestreo en tiempo seco. Fuente: autores del proyecto.

Tabla 30
Hora de la toma de muestras en tiempo de lluvia.

19 de abril 19 de abril 19 de abril 19 de abril 19 de abril 20 de abril

Fecha:
de 2018 de 2018 de 2018 de 2018 de 2018 de 2018
Estación: Cabecera Estación 1 Estación 2 Estación 3 Estación 4 Estación 5
Hora de
8:40:00 a. 7:37:00 a. 11:54:49 a. 4:51:20 p. 11:48:56 p. 2:52:23 a.
toma de
m. m. m. m. m. m.
muestra:
 94

Tabla 31
Datos de calidad del agua en las estaciones de monitoreo en tiempo de lluvia.

Parámetro Unidades Cabecera E1 E2 E3 E4 E5

Temperatura °C 19,2 19,1 21,4 23,3 22,1 21,8

Conductividad μohms 40,10 42,60 48,80 55,50 58,30 58,80

Sólidos Suspendidos
mgD/L 18 21 20 24 19 20
inorgánicos
Oxígeno disuelto mgO2/L 7,80 7,40 7,20 6,60 7,60 7,50
DBO lenta mgO2/L 3,00 3,70 2,00 2,60 3,00 1,30

DBO rápida mgO2/L 0,80 0,10 2,70 1,20 2,30 2,90

NH4- nitrógeno μgN/L 50 850 10 60 420 70

NO3- nitrógeno μgN/L 30 60 110 130 260 150

Fósforo inorgánico
μgP/L 40 270 140 110 30 160
(SRP)
Detritos o Material
Orgánico Particulado mgD/L 0,10 1,70 1,60 0,10 1,30 1,40
MOP
Patógenos Ufc/100mL 3000 115000 7600 6000 74000 17000

Alcalinidad mgCaCO3/L 20 23 25 25 28 28
pH s.u. 6,35 6,92 6,37 6,89 6,91 6,89
Caudal m3/s 5,04 4,47 4,85 3,62 3,18 4,79

Nota: datos obtenidos en la campaña de muestreo en tiempo de lluvia. Fuente: autores del
proyecto.

4.3.8. Perfiles de calidad

La elaboración de los perfiles de calidad de agua permite evaluar gráficamente el

comportamiento de la calidad y cantidad del agua del cauce principal del rio algodonal

(CORPONARIÑO, 2014). De esta manera, las gráficas representan la calidad del agua en cada

estación.
 95

El punto de partida es la estación Cabecera y la desembocadura es la Estación 5. A

continuación, se muestra el análisis de los parámetros.

En la figura 26 se puede evidenciar el aumento del caudal en la época de lluvia, y un

particular comportamiento de la cuenca, donde la Cabecera posee un mayor caudal que la E1 a

204 metros de distancia. Así mismo, el caudal en Cabecera aumentó 3.69 m3/s, en E1 aumentó

3.23 m3/s, E2 3.46 m3/s, E3 2.4 m3/s, E4 3.05 m3/s (126 %) y E5 aumentó 3.24 m3/s. Es decir, el

caudal aumentó más del 200% a excepción de E4.

Q m3/s
6

4 E4
E1 E2 E5
E3 Campaña 1 Q m3/s
2 C
Campaña 2 Q m3/s
0
-1 4 9 14 19 24 29
Distancia entre estaciones Km

Figura 26. Caudal de las estaciones en los dos periodos de muestreo.

Los valores de temperatura del agua están influenciados por la temperatura ambiente.

Esta puede afectar el desarrollo de la vida acuática como se muestra en la tabla 32 Además, las

temperaturas elevadas incrementan la velocidad de las reacciones químicas y disminuye la

concentración de oxígeno disuelto (OD) (López del Pino & Martín Calderón, 2017).
 96

Tabla 32
Efectos de la temperatura en la vida acuática.

Temperatura
Efectos en la vida acuática
(°C)
5 Cesan su actividad las bacterias nitrificantes autótrofas
15 Las bacterias productoras de metano detienen su actividad
25-35 Desarrollo óptimo de la actividad bacteriana
Se detiene los procesos de digestión aerobia y de
50
nitrificación

Temperatura °C
30

25 Campaña 1
E1 Temperatura del
E3 agua °C
20 E4 E5
C E2 Campaña 2
15 Temperatura del
-1 4 9 14 19 24 29 agua °C

Distancia entre estaciones Km

Figura 27. Temperatura de las estaciones en los dos periodos de muestreo.

Generalmente, el agua posee una temperatura óptima para las reacciones de depuración.

En la segunda campaña la temperatura del agua disminuyó.

El pH es una medida de acidez del agua. Los valores límites establecidos por la

legislación se encuentran en el rango comprendido entre 6,5 y 8,5. Algunos factores que pueden

afectar el pH son la polución, el material orgánico, la eutrofización, los componentes químicos

del agua (agua residual, industrial, etc), la composición del sustrato como, por ejemplo: la

dolomita, piedra caliza y la calcita, tienen agua con un nivel de pH básico; y los minerales de
 97

sulfuro, cuarzo y granito otorgará un pH ácido o significativamente menor al agua (Cusiritati,

2018).

pH
8
7.5 C E2 E3 E4 E5
E1
7
Campaña 1 pH
6.5
Campaña 2 pH
6
-1 4 9 14 19 24 29
Distancia entre estaciones Km

Figura 28. pH de las estaciones en los dos periodos de muestreo.

En la campaña 1 el pH tuvo un comportamiento optimo, mientras que en la campaña 2 el

pH bajó, excepto en la E1.

El oxígeno disuelto OD, es una medida de la cantidad de moléculas de oxígeno libre en el

agua y de su concentración depende si se produce la autodepuración. La turbulencia del agua, la

difusión atmosférica, la fotosíntesis de las plantas acuáticas favorecen la oxigenación.

Hay muchos factores que afectan la concentración de oxígeno disuelto en el agua, como

la temperatura, el flujo de la corriente, las plantas acuáticas, la materia orgánica en

descomposición. También se pierde debido a la nitrificación, la respiración fitoplancton, y la

oxidación de la demanda biológica de oxígeno. Los niveles de oxígeno disuelto oscilan entre 7 y

12 mg/l O2. Los niveles bajos de OD pueden encontrarse en áreas donde el material orgánico

está en descomposición (Pérez Carreño, 2015).

 98

OD mgO2/L
E2
9
E4
E1
8 E3 E5
C Campaña 1 OD
7 mgO2/L
Campaña 2 OD
6 mgO2/L
-1 4 9 14 19 24 29
Distancia entre estaciones Km

Figura 29. Oxígeno disuelto en las estaciones en los dos periodos de muestreo.

Como se observa en la figura 29, en la E1 baja el valor de la concentración de oxígeno en

ambas campañas, esto producto de la descomposición de la materia orgánica de los dos

vertimientos en el municipio de Ábrego: vertimiento de la laguna V1 y Vertimiento del barrio

Alto Viejo V2. De igual forma en la E3 y E5 baja la concentración de oxígeno disuelto.

La DBO corresponde a la medida de la cantidad de oxigeno utilizado por los

microorganismos en la mineralización de la materia orgánica biodegradable, en condiciones

aeróbicas, generalmente en un periodo de cinco días a 20 °C (Navarro R., 2007).

DBOf mgO2/L
5 E1
4
3 Campaña 1 DBOf
E2 E3 E4 E5
2C mgO2/L
1
Campaña 2 DBOf
0 mgO2/L
-1 9 19 29
Distancia entre entaciones Km

Figura 30. La DBO rápida de las estaciones en los dos periodos de muestreo.
 99

La DBOf es la DBO rápidamente biodegradable, compuesta principalmente por ácidos

volátiles y azucares, donde la velocidad de degradación es muy rápida (Ronzano & Dapena,

2002).

DBOs mgO2/L
E1
4 E4 E5
3 E2 E3
2 C Campaña 1 DBOs
mgO2/L
1
0 Campaña 2 DBOs
-1 9 19 29 mgO2/L

Distancia entre estaciones Km

Figura 31. DBO lenta de las estaciones en los dos periodos de muestreo.

La DBOs es la DBO lentamente degradable, que está compuesta por la materia orgánica

degradable (partículas y coloides). Donde su degradabilidad tiene una cinética del orden de diez

veces menor que la anterior (Ronzano & Dapena, 2002).

Los sólidos suspendidos inorgánicos son substancias difícilmente degradables, llamados

también minerales. Dentro de estos se incluyen arenas, aceites y sales minerales disueltas en el

agua potable (Trapote Jaume, 2013).

SSI mgD/L
25
E1
20
E3 E5 Campaña 1 SSI
15 mgD/L
E2 E4 Campaña 2 SSI
10 C
-1 4 9 14 19 24 29 mgD/L

Distancia entre estaciones Km

Figura 32. Solidos suspendidos inorgánicos de las estaciones en los dos periodos de muestreo.
 100

Las bacterias coliformes se encuentran en mayor abundancia en la capa superficial del

agua o en los sedimentos del fondo. Por su amplia diversidad el grupo coliformes ha sido divido

en dos grupos: coliformes totales y coliformes fecales, estas últimas se encuentran en los

intestinos y en el excremento de los humanos y animales (Análisis microbiologico Coliformes

Totales.). Entre las fecales tenemos Escherichia coli.

Patógenos Ufc/100mL
150000 E1
120000 E4
Campaña 1
90000 Patógenos
60000 E5 Ufc/100mL
30000 C E2 E3
Campaña 2
0 Patógenos
-1 4 9 14 19 24 29 Ufc/100mL
Distancia entre estaciones Km

Figura 33. Coliformes totales de las estaciones en los dos periodos de muestreo.

Los tres principales problemas ambientales generados por la presencia de compuestos de

Nitrógeno en el sistema acuático son: el aumento de la acidez, el desarrollo de eutrofización y el

aumento de las concentraciones hasta niveles tóxicos que limitan su uso principalmente como

fuentes de agua para consumo humano (Cárdenas Calvachi & Sánchez Ortiz, 2013).

Aznar dice que: “La oxidación aeróbica de los compuestos amoniacales y órgano

nitrogenados, conduce a la formación de nitritos y posteriormente de estos en nitratos, por lo que

un elevado contenido en nitratos y simultáneamente bajo en amonio, indica que se trata de un

agua contaminada hace tiempo”. (Aznar Jiménez, 2000)

El amonio es un compuesto inorgánico de nitrógeno que se forma durante la degradación

biológica de compuestos orgánicos de nitrógeno. En general, la presencia de amoníaco libre o

 101

ion amonio es considerado como una prueba química de contaminación reciente y peligrosa, es

por esto que la medición del nitrógeno amoniacal es utilizado para indicar la contaminación del

agua y la presencia de aguas residuales no tratadas.

NH4-N μgN/L
1000 E1
750

500 E4 Campaña 1 NH4-N

μgN/L
250
Campaña 2 NH4-N
C E3 E5
0 E2 μgN/L
-1 4 9 14 19 24 29
Distancias entre estaciones Km

Figura 34. Amonio de las estaciones en los dos periodos de muestreo.

En presencia del oxígeno el amonio se puede convertir en nitrato (nitrificación). Los

nitratos provienen de la contaminación orgánica (deshechos de origen humano o industrial) o de

la contaminación por abonos químicos (fertilizantes, estiércol…). El valor máximo permitido por

la legislación para un agua de consumo de una red de abastecimiento es de 50 mg/L. (PanaChlor,

2018)

NO3-N μgN/L
300
250 E4
200
E3 E5 Campaña 1 NO3-N
150 E2
100 E1 μgN/L
50 Campaña 2 NO3-N
C
0 μgN/L
-1 4 9 14 19 24 29
Distancia entre estaciones Km

Figura 35. Nitratos de las estaciones en los dos periodos de muestreo.

 102

La conductividad es una medida de la propiedad que tiene el agua de dejar pasar la

corriente eléctrica. Permite conocer la mineralización de un agua, debido a la concentración de

sales disueltas y presencia de iones. El valor máximo permitido por la legislación para un agua

de consumo de una red de abastecimiento es de 2500µS/cm. (PanaChlor, 2018)

Conductividad μohms
90 E5
E1 E3 E4
80 E2
Campaña 1
70
Conductividad
60 μohms
50 C
Campaña 2
40
Conductividad
-1 4 9 14 19 24 29 μohms
Distancia entre estaciones Km

Figura 36. Conductividad de las estaciones en los dos periodos de muestreo.

El material orgánico particulado (detritos), son residuos sólidos que provienen de la

descomposición de fuentes orgánicas. Aunque es considerada materia muerta es aún materia

orgánica en descomposición, ya que los seres vivos se alimentan de ella degradándola aún más

(Ceccon, 2013, pág. 212).

MOP mgD/L
2
1.5
E1
1 Campaña 1 MOP
E2
E3 E5 mgD/L
0.5 C E4
Campaña 2 MOP
0 mgD/L
-1 4 9 14 19 24 29
Distancia entre estaciones Km

Figura 37. Detritos (MOP) de las estaciones en los dos periodos de muestreo.
 103

El fosforo inorgánico, está representado por el fosforo presente en minerales primarios.

Las bacterias satisfacen normalmente sus necesidades de fósforo extrayéndolo del detritus

orgánico. Sin embargo, algunas bacterias lo extraen directamente desde el agua.

Pinor μgP/L
400
300 E1

200 Campaña 1 Pinor

E2 E5
E3 μgP/L
100
C Campaña 2 Pinor
E4
0 μgP/L
-1 4 9 14 19 24 29
Distancia entre estaciones Km

Figura 38. P inorgánico de las estaciones en los dos periodos de muestreo.

La alcalinidad en el agua es su capacidad de neutralizar ácidos. La alcalinidad del agua se

debe entonces principalmente a sales de ácidos débiles y bases fuertes y, estas sustancias actúan

como amortiguadoras para resistir la caída de pH resultante a la adición de ácidos. Por lo general,

se debe fundamentalmente a su contenido de carbonatos, bicarbonatos e hidróxidos aunque otras

sales o bases también contribuyen a la alcalinidad (Vanessa Callejas & Lopez, 2015).

Alcalinidad mgCaCO3/L
50 E2 E3
E1 E4 E5
40
Campaña 1
30 C
Alcalinidad
20 mgCaCO3/L
10 Campaña 2
0 Alcalinidad
-1 4 9 14 19 24 29 mgCaCO3/L
Distancia entre estaciones Km

Figura 39. Alcalinidad de las estaciones en los dos periodos de muestreo.

 104

4.4. Modelación de la capacidad de autodepuración

4.4.1. Topología del río Algodonal.

El perfil de modelación incluye únicamente los aportes o entradas al río Algodonal que se

consideraron como significativos para la modelación.

Figura 40. Topología del río Algodonal.

4.4.2. Hojas de entrada al modelo.

 105

Qual2k.

Figura 41. Ejemplo de la hoja Qual2k.

Headwater.

Se asume que el comportamiento del río es uniforme, por lo que se agrega los datos de los

parámetros fisicoquímicos y microbiológicos en todas las horas del modelo.

Figura 42. Ejemplo de la hoja Headwater.

 106

Reach.

Información relacionada con la cabecera del río y sus alcances (segmentos). Y por

razones de funcionamiento del modelo Q2k, la distancia se representa en orden decreciente de

E5 a Cabecera.

Figura 43. Ejemplo de la hoja Reach.

Air Temperature.

Se introduce la temperatura del aire horaria.

Figura 44. Ejemplo de la hoja Air Temperature.

Cloud clover.

Porcentaje de nubosidad horario.

 107

Figura 45. Ejemplo de la hoja Cloud clover.

Shade.

Porcentaje de sombra se tomó con un valor de cero porque no existía arboles grandes y

con follaje alrededor de las estaciones.

Figura 46. Ejemplo de la hoja Shade.

Rates.

Constantes.

Figura 47. Contantes de calibración proporcinadas por el modelo.

 108

Point sources.

Se ingresan los valores medios de calidad del agua de los vertimientos puntuales y las

captaciones.

Figura 48. Ejemplo de la hoja Point Sources.

Hydraulics Data.

Información de caudal, ancho de la lámina de agua, velocidad media y tiempo de viaje.

Tiene unas celdas adicionales de color verde y rosa para calibración del Q2k.
 109

Figura 49. Ejemplo de la hoja Hidraulis data.

Temperature Data

Temperatura media del agua en las estaciones de monitoreo.

Figura 50. Ejemplo de la hoja Temperatura data.

WQ Data.

Información relacionada los parámetros de calidad del agua de la cabecera del río y sus

alcances (segmentos).
 110

Figura 51. Ejemplo de la hoja WQ data.

Celdas de calibración de la hija WQ Data.

Figura 52. Celdas de calibración de la hoja WQ data.

4.4.3. Escenario actual de calidad del cauce principal del río Algodonal.

Una vez aplicado y calibrado el modelo de simulación Qual2k con los datos de campo e

información disponible en la literatura, se obtuvieron los resultados del estado actual de calidad

del cauce del río principal desde aguas arriba de la laguna de oxidación del municipio de Ábrego,

hasta la captación de agua de la empresa ESPO S.A.S. del municipio de Ocaña, para la época

seca (campaña 1) y la época de lluvias (campaña 2).

 111

Por otro lado, solo se analizó el vertimiento puntual de la laguna de oxidación de Ábrego

en cuanto a DBO5, caudal vertido y temperatura. En este caso, se hallaron los demás parámetros

por literatura. También se tuvo en cuenta el caudal de las captaciones que afectaron el caudal de

la corriente principal que poseían concesión vigente.

A continuación, se presentan las gráficas obtenidas para tiempo de viaje, oxígeno

disuelto, temperatura, pH, en la campaña 1.

Río Algodonal (Ábrego - Ocaña, Colombia)

(3/13/2018) Mainstem
1.20

1.00

0.80

0.60

0.40

0.20

0.00
30 25 20 15 10 5 0
trav time, d Travel time data (d)

Figura 53. Tiempo de viaje río Algodonal (m3/s) en la campaña 1.

 112

Río Algodonal (Ábrego - Ocaña, Colombia)

9 (3/13/2018) Mainstem

6
30.0 25.0 20.0 15.0 10.0 5.0 0.0
DO(mgO2/L) DO (mgO2/L) data
DO(mgO2/L) Min DO(mgO2/L) Max
Minimum DO-data Maximum DO-data

Figura 54. Oxígeno disuelto OD río Algodonal (mgO2/L) en la campaña 1.

Río Algodonal (Ábrego - Ocaña, Colombia)

35 (3/13/2018) Mainstem

30

25

20

15
30.00 25.00 20.00 15.00 10.00 5.00 0.00
Temp(C) Average Mean Temp-data Temp(C) Minimum
Temp(C) Maximum Minimum Temp-data Maximum Temp-data

Figura 55. Temperatura del río Algodonal (°C) en la campaña 1.

 113

Río Algodonal (Ábrego - Ocaña, Colombia)

9.0 (3/13/2018) Mainstem

8.0

7.0

6.0

5.0
35 30 25 20 15 10 5 0
pH pH data
pH Min pH Max
Minimum pH-data Maximum pH-data
pHsat

Figura 56. Potencial de hidrógeno pH río Algodonal en la campaña 1.

A continuación, se presentan las gráficas obtenidas para para tiempo de viaje, oxígeno

disuelto, temperatura, pH en la campaña 2.

Río Algodonal (Ábrego - Ocaña, Colombia)

(4/19/2018) Mainstem
1.00

0.80

0.60

0.40

0.20

0.00
30 25 20 15 10 5 0
trav time, d Travel time data (d)

Figura 57. Tiempo de viaje río Algodonal (m3/s) en la campaña 2.

 114

Río Algodonal (Ábrego - Ocaña, Colombia)

9 (4/19/2018) Mainstem

6
30.0 25.0 20.0 15.0 10.0 5.0 0.0
DO(mgO2/L) DO (mgO2/L) data DO(mgO2/L) Min
DO(mgO2/L) Max Minimum DO-data Maximum DO-data
DO sat

Figura 58. Oxígeno disuelto OD río Algodonal (mgO2/L) en la campaña 2.

Río Algodonal (Ábrego - Ocaña, Colombia)

30 (4/19/2018) Mainstem

25

20

15
30.00 25.00 20.00 15.00 10.00 5.00 0.00
Temp(C) Average Mean Temp-data Temp(C) Minimum
Temp(C) Maximum Minimum Temp-data Maximum Temp-data

Figura 59. Temperatura del río Algodonal (°C) en la campaña 2.

 115

Río Algodonal (Ábrego - Ocaña, Colombia)

9.0 (4/19/2018) Mainstem

8.0

7.0

6.0

5.0
30 25 20 15 10 5 0
pH pH data
pH Min pH Max
Minimum pH-data Maximum pH-data
pHsat

Figura 60. Potencial de hidrógeno pH río Algodonal en la campaña 2.

Las gráficas permiten relacionar los datos de entrada (estaciones de monitoreo)

representados con puntos negros con la tendencia de calibración del modelo representado en

líneas.

La línea tendencial del tiempo de viaje está directamente influenciada por el coeficiente

de Manning en cada estación de monitoreo. En las dos campañas presenta un comportamiento

lineal.

El comportamiento del oxígeno disuelto simulado concuerda con los datos de las

estaciones y representa decrecimiento en los vertimientos de aguas residuales domesticas al

cauce principal del río Algodonal.

La temperatura en el río Algodonal en la campaña 1 aumentaron en función de los

vertimientos. Sin embargo, en la campaña 2 se mantuvo la temperatura desde cabecera hasta la

 116

última estación debido a que el vertimiento no era significativo comparándolo con el nivel de

caudal del río.

En general, el potencial de hidrogeno se mantiene en un rango óptimo para los procesos

de depuración.

De esta manera, el escenario actual de la quebrada permite la depuración de los

vertimientos de aguas residuales, gracias a las características fisicoquímicas e hidráulicas del río.

Los demás parámetros ingresados al modelo, no fue posible calibrarlos debido a la

insuficiente información sobre las características del agua residual doméstica, los afluentes y las

extracciones (captaciones).

4.4.4. Calidad del agua en el río Algodonal frente a los escenarios prospectivos.

Bajo las condiciones actuales en tiempo seco se simulan 3 escenarios de comportamiento y

calidad del agua:

 Escenario 1: río Algodonal frente a la descarga de vertimientos de agua residual

domestica sin tratamiento previo, con concentraciones mínimas.

 Escenario 2: río Algodonal frente a la descarga de vertimientos de agua residual

domestica sin tratamiento previo, con concentraciones medias.

 Escenario 3: río Algodonal frente a la descarga de vertimientos de agua residual

domestica sin tratamiento previo, con concentraciones máximas.

Enseguida, se realiza la descripción de los resultados obtenidos en los escenarios en función

del Oxígeno disuelto:

 117

Comportamiento del oxigeno disuelto frente a los escenarios

Época seca
ESTANDAR DO(mgO2/L) MIN DO (mgO2/L) MED DO (mgO2/L) MAX DO (mgO2/L)

8.00
mg O2/L

7.00

6.00
0.00 5.00 10.00 15.00 20.00 25.00 30.00
Distancia entre estaciones km

Figura 61. Comportamiento del OD frente a los escenarios prospectivos.

Los escenarios simulados afectan la calidad del agua del tramo de estudio a medida del

aumento de las concentraciones de los vertimientos. Aunque, el río representa un

comportamiento de depuración aceptable.

 118

Capítulo 5 Conclusiones

La selección del modelo se basó en las características de la cuenca a modelar, las

características generales del modelo, la capacidad de simulación de los parámetros básicos de

calidad, la facilidad de calibración, sus antecedentes y la viabilidad económica.

Con respecto a la caracterización del rio Algodonal, se identificaron 3 vertimientos

significativos de aguas residuales: el de la laguna de estabilización de ábrego y el del barrio Alto

viejo no vierten directamente al río Algodonal. Por otro lado, el vertimiento 3 que corresponde al

corregimiento de la ermita si vierte directamente pero no existe un sistema de alcantarillado que

las agrupe, dificultando su localización y aforo.

En la campaña de aforo se caracterizó la hidráulica de las estaciones en función de los

parámetros que exigía el modelo para calcular la profundidad y la velocidad según la ecuación de

Manning. Asimismo, se evidencio un aumento mayor al 200% en el caudal y problemas de

sedimentación en el río.

En cuanto a la calidad del agua representada en el perfil de calidad, se evidencia que el

comportamiento del río se mantiene estable en la distancia en cuanto a la temperatura. El

oxígeno disuelto se ve afectado en los vertimientos, pero luego se recupera. El pH en tiempo de

lluvias presenta un decrecimiento, pero se estabiliza en la estación 3. En cuanto a la DBOf o

rápidamente biodegradable y la DBOs o lentamente degradable mantuvo un comportamiento

similar en las dos campañas de monitoreo.

Los sólidos suspendidos inorgánicos y los detritos aumentaron en la campaña 2 producto

de los problemas de sedimentación de la cuenca.

 119

Las grafica de las coliformes (bacterias) permitió determinar que los patógenos no se ven

afectadas por el aumento del caudal. Se resaltan que las unidades formadoras de colonia de las

coliformes fecales no se identificaron, sin embargo, la gráfica muestra un aumento significativo

en las estaciones aguas debajo de los vertimientos identificados. Asimismo, se resalta la cantidad

de coliformes presentes en la bocatoma de la planta de tratamiento de agua potable de la empresa

ESPO S.A.S. donde por normatividad las coliformes totales y fecales deben ser 0.

Las coliformes, el amonio y los altos niveles de nitratos indican contaminación del agua y

presencia de aguas residuales no tratadas.

En cuanto a la conductividad, no sobrepasa el límite máximo permitido para el uso del

agua de consumo humano.

El contenido de alcalinidad está por encima de la concentración mínima para mantener la

vida acuática (20 mg de CaCO3/L).

Las dos campañas de monitoreo realizadas al río Algodonal y la información utilizada por

la literatura, permitieron calibrar y confirmar el modelo Q2k en cuanto a tiempo de viaje,

oxígeno disuelto, temperatura, pH.

Los resultados permitieron evidenciar que Cabecera- E1 es un tramo crítico. Con

problemas de contaminación generados por vertimientos de aguas residuales de origen domestico

con y sin tratamiento.

Como se ha mostrado, los resultados de la modelación representan una línea base en

términos de calidad del agua en el río Algodonal, puesto que no existen experiencias previas de

modelación en la implementación del modelo. De ello resulta necesario decir que lo resultados

del modelo estuvieron influenciados por la información adicional de los vertimientos y

 120

captaciones, siendo estas aún escazas para la modelar el comportamiento del flujo del agua y los

demás parámetros hidráulicos.

A partir de los resultados del comportamiento de los escenarios de simulación propuestos,

se puede determinar que la capacidad de depuración del río Algodonal es buena. Además, cabe

resaltar que no existen vertimientos que pongan en riesgo los parámetros fisicoquímicos, pero sí,

los microbiológicos (coliformes). Por otro lado, se debe tener en cuenta la aplicación de modelos

de calidad del agua como lo menciona el Decreto 3930 de 2010 con el fin de planificar el

ordenamiento del recurso hídrico y proponer medidas de prevención de la contaminación en la

fuente hídrica. También se concluye que, a pesar del aumento de los vertimientos, no se afecta la

capacidad máxima de asimilación del río Algodonal.

El modelo de calidad del agua Qual2k, puede ser ajustado a otras campañas de monitoreo

por sus antecedentes de éxito. En especial, por la recomendación de la Corporación Autónoma

Regional de Nariño como herramienta de planificación del recurso hídrico.

Finalmente se concluye que la campaña en tiempo seco permite la calibración del modelo

y la campaña en tiempo de lluvia permite la confirmación del mismo.

 121

Capítulo 6 Recomendaciones

Se recomienda que se revise detalladamente los datos de entrada al modelo de calidad,

con el fin de obtener los datos de calidad del agua con análisis de laboratorio y no incurrir a

supuestos de literatura.

Se aconseja caracterizar de manera detallada todos los vertimientos, quebradas o

captaciones, que influyan en el caudal y los parámetros de calidad del río a modelar.

Aunque el modelo exige los parámetros de forma horaria se puede realizar la

caracterización del río en varios días, pero en la misma época (régimen de lluvias o tiempo seco),

tratando de que se caracterice en el menor tiempo posible y se considera que los parámetros de

calidad son constantes durante todo el día. Ya que tomar parámetros continuos en 24 horas y en

diferentes estaciones requiere capital humano, económico y tecnológico.

 122

Capítulo 7 Referencias

Alcaldía Ábrego. (2016-2019). Plan de Desarrollo Municipal "Conmigo la agente gobierna".

Ábrego.
Alcaldía Ocaña. (2002-2011). Plan Básico de Ordenamiento Territorial. Ocaña.
Alcaldía Ocaña. (2016-2019). Plan de desarrollo municipal "Es la hora de Ocaña".
Alean Vera, L. E. (2009). Modelación de la capacidad máxima de asimilación de vertimientos de
carga orgánica en la ciénaga de Mesolandia en el departamento del Atlántico. Tesis de
maestría, Bogotá. Recuperado el 23 de Agosto de 2017, de
http://www.bdigital.unal.edu.co/8792/
Análisis microbiologico Coliformes Totales. (s.f.). Obtenido de Calidad microbiológica:
http://www.calidadmicrobiologica.com.co/microbiologia/coliformes-totales
Anguiano Valencia, J. R. (2012). Determinación del coeficiente de reaireación en un tramo
corto de un río caudaloso con flujo turbulento, mediante el uso de los trazadores
hexafluoruro de azufre (SF) y Rodamina WT. Tesis de maestría, Universidad Michoacana
de San Nicolás de Hidalgo, Morelia. Recuperado el 31 de Octubre de 2017, de
https://posgrado-mcia.umich.mx/tesis/2012/2012-0006.pdf
Arenas Sánchez, A., & Díaz Jimenez, D. A. (2018). Estructuración de una propuesta de
esquema de pagos por servicios ambientales para la protección y recuperación del
páramo de Jurisdicciones, Norte de Santander. Trabajo de grado, Universidad Francisco
de Paula Santander, Ocaña.
Asamble Nacional Constituyente. (13 de Junio de 1991). Constitución Política de Colombia.
Colombia.
Autodepuración y tratamiento de las aguas residuales. (6 de Febrero de 2018). Obtenido de Sites
Google: https://sites.google.com/site/hidrosfera11/home/autodepuracion-y-tratamiento-
de-las-aguas-residuales
Autoridad Nacional de Licencias Ambientales ANLA. (2013). Metodología para la definición de
la longitud de influencia de vertimientos sobre corrientes de agua superficial.
Recuperado el 5 de Marzo de 2018, de
http://www.anla.gov.co/documentos/institucional/Metodologia_-
_Longitud_de_Influencia_de_Vertimientos.pdf
Aznar Jiménez, A. (2000). Determinación de los parámetros fisico-químicos de calidad de las
aguas. Revista interdisciplinar de gestión ambiental, 12-19.
Bedoya Contreras, J. F. (2007). Modelos de simulacion de transporte de metales pesados en la
cuenca baja del río Tunjuelo. Tesis de grado, Universidad de La Salle, Facultad de
ingeniería ambiental y sanitaria, Bogotá. Recuperado el 15 de Octubre de 2017, de
 123

http://repository.lasalle.edu.co/bitstream/handle/10185/14332/T41.07%20B391m.pdf?seq
uence=1
Campos Aranda, D. F. (1998). Procesos del ciclo hidrológico. Mexico: Editorial Universitaria
Potosina.
Cárdenas Calvachi, G. L., & Sánchez Ortiz, I. A. (2013). Nitrógeno en aguas residuales:
orígenes, efectos y mecanismos de remoción para preservar el ambiente y la salud
pública. Scielo, XV(1).
Casas C., W., & Vargas Bejarano, C. (24 de Diciembre de 2008). Modelos matemáticos de
simulación de la calidad del agua en Colombia. Prinicipios y aplicaciones. Obtenido de
Biblovirtual MinAMbiente:
http://biblovirtual.minambiente.gov.co:3000/DOCS/MEMORIA/MMA-0013/MMA-
0013-CAPITULO6.pdf
Castaño Ossa, J. M. (2015). Modelación de la calidad del agua del agua del río Tunjuelo en dios
escenarios de implementación del plan de saneamiento de Bogotá. Tesis de grado,
Universidad Nacional de Colombia, Departamento de ingeniería civil y agrícola, Bogotá.
Recuperado el 15 de Octubre de 2017, de http://www.bdigital.unal.edu.co/52643/
Castillo G., J. (2004). Modelos de calidad de aguas. CI51D-1 Contaminación de Recursos
Hidráulicos 2004, Primavera. Recuperado el 15 de Octubre de 2017, de https://www.u-
cursos.cl/ingenieria/2004/2/CI51D/1/material_docente/
Castro de Esparza, M. L. (1995). Parámetros Físicos, Carga Orgánica y Nutrientes. Lima, Perú.
Recuperado el 3 de Marzo de 2018, de
http://www.ib.usp.br/limnologia/textos/CEPIS%20parame.pdf
Castro Huertas, M. A. (2015). Aplicación del QUAL2Kw en la modelación de la calidad del
agua del río Guacaica, departamento de Caldas, Colombia. Trabajo de grado,
Universidad Nacional de Colombia, Departamento de Ingeniería Química, Manizales.
Recuperado el 31 de Agosto de 2017, de http://www.bdigital.unal.edu.co/51032/
Ceccon, E. (2013). Restauración en Bosques tropicales: fundamentos ecológicos, prácticos y
sociales.
CLIMATE-DATA.ORG. (28 de Febrero de 2018). Clima: Ocaña. Obtenido de Climate-data.org:
https://es.climate-data.org/location/50064/
Conferencia de las Naciones Unidas sobre Ambiente y Desarrollo Río de Janeiro. (1992).
Programa XXI. Recuperado el 2018 de 01 de 25, de
http://www.un.org/spanish/esa/sustdev/agenda21/index.htm
Contaminación-Agua.org. (s.f.). Contaminación-Agua.org. Obtenido de Contaminación-
Agua.org: https://contaminacion-agua.org/definicion/
Contreras López, A., & Molero Meneses, M. (2012). Ciencia y Tecnología del Medioambiente.
Madrid: Univercidad Nacional de Educación a Distancia. Obtenido de
 124

https://books.google.com.co/books?id=xtvNPP8X6KAC&pg=PA21&dq=La+capacidad+
de+autoregeneraci%C3%B3n+de+un+r%C3%ADo&hl=es&sa=X&ved=0ahUKEwiljrve
15TZAhUxq1kKHXKQDu4Q6AEIJzAA#v=onepage&q&f=false
CORPONARIÑO. (2014). Plan de ordenamiento territorial del cauce principal de la quebrada
LA Fragua. San Juan de Pasto. Recuperado el 1 de Julio de 2018, de
http://corponarino.gov.co/expedientes/descontaminacion/porhqlafragua.pdf
CORPONOR. (2011). Sintesis ambiental del Norte de Santander.
Corporación Autónoma Regional de la Frontera Nororiental. (s.f.). Objetivos y Funciones.
Ocaña. Obtenido de http://corponor.gov.co/es/index.php/nuestra-entidad/la-
entidad/objetivo-y-funciones
Cusiritati. (23 de Mayo de 2018). Los factores naturales que afectan el pH en agua. Obtenido de
Cusiritati: http://www.cusiritati.com/pBoMaQ9mN/
Dataset: © JAXA / METI. (2010). ALOS PALSAR [L1.0]. Recuperado el 15 de Abril de 2018,
de https://www.asf.alaska.edu], se accede a través de ASF DAAC.
Decreto N° 1076. (26 de Mayo de 2015). Diario Oficial de Colombia. Bogotá, Cundinamarca,
Colombia.
Decreto N° 1594. (26 de Julio de 1984). Diario Oficial 36700 de la República de Colombia.
Bogotá, Colombia.
Decreto N° 2667. (21 de Diciembre de 2012). Diario Oficial 48651 de la República de Colombia.
Bogotá, Colombia.
Decreto N° 3930. (25 de Octubre de 2010). Diario Oficial 47837 de la República de Colombia.
Bogotá, Colombia.
Diaz Merchán, B. E. (2004). Modelación de la calidad del agua en el interceptor río Bogotá en
los tramos Fucha-Tunjuelo-Canoas. Tesis de maestría, Universidad de los Andes,
Departamento de ingeniería civil y ambiental, Bogotá. Recuperado el 15 de Octubre de
2017, de http://oab2.ambientebogota.gov.co/es/documentacion-e-
investigaciones/resultado-busqueda/modelacion-de-la-calidad-del-agua-en-el-interceptor-
rio-bogota-en-los-tramos-fucha-tunjuelo-canoas
Domínguez Calle, E. A., Gonzalo Rivera, H., Vanegas Sarmiento, R., & Moreno, P. (Junio de
2008). Demanda-oferta de agua y el índice de escasez de agua como herramientas de
evaluación del recurso hídrico colombiano. Revista de la Academia Colombiana de
Ciencias Exactas, Físicas y Naturales, XXXII(123), 195-212. Recuperado el 20 de
Octubre de 2017, de https://www.researchgate.net/publication/228463075
Duarte Silva, C. C., Pico Pinto, G. S., & Estupiñan Pinto, R. A. (2014). Determinación de la
capacidad de autodepuración en un tramo del río Fonce en San Gil con bse en el balance
de oxígeno disuelto. Revista Matices Tecnológicos, VI, 30-35. Recuperado el 15 de
 125

Octubre de 2017, de http://publicaciones.unisangil.edu.co/index.php/revista-matices-

tecnologicos/article/view/257/220
Enderlein, U. S., Enderlein, R. E., & Williams, W. P. (s.f.). Capitulo 2: Requisitos de la calidad
del agua.
García Quintero , H. A. (2008). Modelación de la calidad del agua en corrientes naturales
mediante sistemas no lineales acoplados. Tesis de maestría, Pontificia Universidad
Javeriana. Recuperado el 31 de Octubre de 2017, de
https://repository.javeriana.edu.co/bitstream/handle/10554/3809/GarciaQuinteroHermesA
lfonso2008.pdf?sequence=1
García Quintero, H. A. (2008). Modelación de la calidad del agua en corrientes naturales
mediante sistemas no lineales acopladas, empleando dos algoritmos genéticos en el
procesos de calibración y solución del sistema. Tesis de Maestría, Pontificia Universidad
Javeriana, Bogotá. Recuperado el 3 de Marzo de 2018, de
https://repository.javeriana.edu.co/bitstream/handle/10554/3809/GarciaQuinteroHermesA
lfonso2008.pdf?sequence=1&isAllowed=y
Gaspari, F. J., Rodríguez Vagaríaa, A. M., Senisterra, G. E., Denegri, G., Delgado, M., &
Besteiro, S. (2013). Caracterización morfométrica de la cuenca alta del río Sauce Grande,
Buenos aires, Argentina. AUGMDOMUS, 143-158.
Gil Rodríguez, M. (2006). Depuración de aguas residuales: modelización de procesos de lodos
activos (Vol. XXIV). Madrid: CSIC. Recuperado el 29 de Enero de 2018, de
https://books.google.com.co/
Grupo de Investigación en Asustos Ambientales y Urbanos [MINDALA]. (2017). Análisis de
calidad del agua al vertimiento de la Laguna de Oxidación. Ábrego.
HIMAT. (1978). Resolución 00337 "Por la cuál se adopta un sistema de codificación para las
estaciones Hidrometeorológicas". Obtenido de
http://www.ideam.gov.co/documents/14691/16498/Resolucion+337_1978.pdf/f17e6dfa-
794c-40ca-90d3-e1c9ec1edf42
Ibáñez Asensio, S., Moreno Ramón, H., & Gisbert Blanquer, J. M. (2008). Morfología de las
cuencas hidrográficas. Universidad Politecnica de Valencia, Valencia.
IDEAM. (2007). Protocolo para el monitoreo y seguimiento del agua. Recuperado el 6 de Marzo
de 2018, de
http://documentacion.ideam.gov.co/openbiblio/bvirtual/021172/Protocoloparaelmonitore
oyseguimientodelagua.pdf
IDEAM. (2013). Zonificacion y Codificación de Cuencas Hidrográficas. Bogotá. Obtenido de
http://documentacion.ideam.gov.co/openbiblio/bvirtual/022655/MEMORIASMAPAZON
IFICACIONHIDROGRAFICA.pdf
 126

Jacobs, M. (1996). La economía verde. Medio ambiente, desarrollo sostenible y la política del
futuro. (Vol. XII). Barcelona: Icaria. Obtenido de
https://books.google.com.co/books?id=-
Ag1wbUrYj8C&pg=PA169&lpg=PA169&dq=La+capacidad+de+autoregeneraci%C3%
B3n+de+un+r%C3%ADo+depende&source=bl&ots=EvOF6oEfQY&sig=PQVSTCdbiuh
imKQ0bcZfxk5KVgo&hl=es&sa=X&ved=0ahUKEwil9b7B05TZAhWCxFkKHb8JBC0
Q6AEIRjAF#v=onepag
Jiménez Beltrán, A. A., & Cortes Naranjo, D. A. (2016). Modelación dinámica de la calidad del
agua del río Chicamocha teniendo en cuenta el vertimiento de agua residual de la ciudad
de Tunja. Tesis de grado, Universidad de Santo Tomás, División de Ingenierías, Bogotá.
Recuperado el 23 de Agosto de 2017, de
http://repository.usta.edu.co/bitstream/handle/11634/2843/Jimenezangie2016.pdf?sequen
ce=1&isAllowed=y
Jiménez, B. E. (2001). La contaminación ambiental en México.
Jørgensen, S. E., & Bendoricchio, G. (2001). Fundamentals of Ecological Modelling (Tercera
ed., Vol. XXI). Elsevier. Recuperado el 31 de Enero de 2018, de
https://books.google.es/books?hl=es&lr=&id=hd7aPFIdvTcC&oi=fnd&pg=PP2&dq=Fun
damentals+of++%22ecological+modelling%22+&ots=O5ehPde9eh&sig=DLIf8zYtvSHP
M6SlZD0SSS6xmVs#v=onepage&q&f=false
Ley N° 99. (22 de Diciembre de 1993). Diario Oficial 41146 de la República de Colombia.
Bogotá, Colombia.
Llavador Colomer, F. (2002). Modelos matemáticos de sistemas acuáticos dinámicos.
Publicaciones Universidad de Alicante.
López del Pino, S. J., & Martín Calderón, S. (2017). Depuración de aguas resiudales. Editorial
Elearning.
Lozano, G., Zapata, M., & Peña, L. (2003). Modelación de corrientes hídricas superficiales en el
departamento del Quindío. Seminario internacional la hidroinformática en la gestión
integrada de los recursos hídricos, 57-64. Recuperado el 14 de Octubre de 2017, de
http://www.bvsde.paho.org/bvsacd/agua2003/model.pdf
Lux Cardona, B. (2013). Conceptos básicos de Morfometría de Cuencas Hidrográficas. Maestría
en Energía y Ambiente. Obtenido de
http://www.repositorio.usac.edu.gt/4482/1/Conceptos%20b%C3%A1sicos%20de%20Mo
rfometr%C3%ADa%20de%20Cuencas%20Hidrogr%C3%A1ficas.pdf
Mancera Pineda, J. E., Peña Salamanca, E. J., Giraldo Henao, R., & Santos Martínez, A. (2003).
Introducción a la modelación ecológica. Principios y aplicaciones. Santafé de Bogotá.
Recuperado el 31 de Octubre de 2017, de https://books.google.com.co/books?id=-AfcTs-
AtvUC&printsec=frontcover&hl=es&source=gbs_ge_summary_r&cad=0#v=onepage&q
&f=false
 127

Menéndez, A. N. (2010). Transporte de contaminantes en el medio acuático. Maestría en

Ingeniería Ambiental, Universidad Tecnológica Nacional. Recuperado el 2018 de Enero
de 25, de
http://www.redmodelacion.com.ar/~modelacion/images/pdf/docencia/utn/Apunte_Conta
mAguas.pdf
Min Ambiente. (s.f.). Ministerio de Ambiente y Desarrollo Sostenible. Recuperado el 30 de
Octubre de 2017, de Administración del recurso hídrico:
http://www.minambiente.gov.co/index.php/gestion-integral-del-recurso-
hidrico/administracion-del-recurso-hidrico
Ministerio de Ambiente, Vivienda y Desarrollo Territorial. (2004). Plan Nacional de Manejo de
Aguas Residuales Municipales en Colombia. Bogotá. Obtenido de
http://portalterritorial.gov.co/apc-aa-
files/7515a587f637c2c66d45f01f9c4f315c/PLAN_NACIONAL_DE_MANEJO_DE_AG
UAS_RESIDUALES_MUNICIPALES_EN_COLOMBIA.pdf
Molero Fernandez, J., Saez Mercader, J., & Soler Andres, A. (1981). La autodepuracion en las
corrientes de agua. Recuperado el 16 de Julio de 2018, de Red de información educativa:
http://redined.mecd.gob.es/xmlui/bitstream/handle/11162/72890/00820073007882.pdf?se
quence=1
Montes, R. T., Navarro, I., Domínguez Mora, R., & Jiménez, B. E. (Noviembre-Diciembre de
2013). Changes in the self-purification capacity of the Magdalena River due to climate
change. Tecnologia y Ciencias del Agua, IV(5), 71-83. Recuperado el 3 de Marzo de
2018, de https://www.scopus.com/inward/record.uri?eid=2-s2.0-
84902188516&partnerID=40&md5=ab77fd08ec5b1da0b25b892cb23781e1
Navarro R., M. O. (2007). Demanda bioquímica de oxígeno 5 días, incubación y electrometría.
Obtenido de
http://www.ideam.gov.co/documents/14691/38155/Demanda+Bioqu%C3%ADmica+de+
Ox%C3%ADgeno..pdf/ca6e1594-4217-4aa3-9627-d60e5c077dfa
Organización de las naciones unidas para la agricultura y la alimentación (FAO). (2004).
Caracterización de peligros de patógenos en los alimentos y el agua. Roma. Recuperado
el 31 de Octubre de 2017, de http://www.fao.org/3/a-y4666s.pdf
PanaChlor. (2018). Interpretación de resultado de análisis de aguas. Obtenido de
http://panachlor.com/?p=593
Pelletier, G., & Chapra, S. (2008). A Modeling framework for simulating river and stream water
quality. Manual de Usuario (Versión 5..1).
Pérez Carreño, J. (2015). Modelación de la calidad del agua quebrada La Florida, municipio de
Armenia departamento del Quindío. CORPORACIÓN AUTÓNOMA REGIONAL DEL
QUINDÍO CRQ.
 128

Pérez, A., & Torres, P. (Septiembre de 2008). Evaluación del comportamiento hidrodinámico
como herramienta para optimización de reactores anaerobios de crecimiento en medio
fijo. Revista facultad de Ingeniería(45), 27-40. Recuperado el 1 de Septiembre de 2017,
de http://rccp.udea.edu.co/index.php/ingenieria/article/view/17963
POMCHRA. (2010). Plan de Ordenamiento y Manejo de la Cuenca Hidrográfica del Río
Algodonal. Ocaña. Obtenido de
http://corponor.gov.co/publica_recursos/pomca/algodonal/PLAN_DE_ORDENACION_
Y_MANEJO_CUENCA_RIO_ALGODONAL.pdf
Resolución N° 631. (18 de Abril de 2015). Diario Oficial 49486 de la República de Colombia.
Bogotá, Colombia.
Resolución N°330. (17 de Junio de 2017). Diario Oficial No. 50267 de la República de
Colombia. Bogotá D.C., Colombia.
Revista Ambientum. (Marzo de 2002). Revista Ambientum. Obtenido de Revista Ambientum:
https://www.ambientum.com/revista/2002_11/AUTDPRCNGS3.asp
Romo Saltos, L. A. (2003). AGUA: Patrimonio de la Humanidad. Quito: Editorial Universitaria.
Ronzano, E., & Dapena, J. L. (2002). Tratamiento Biológico de las aguas Residuales.
Recuperado el 2018 de Enero de 28, de
https://books.google.es/books?hl=es&lr=&id=M0ClCgAAQBAJ&oi=fnd&pg=PA5&dq=
Medidad+de+la+contaminaci%C3%B3n+Org%C3%A1nica+ronzano&ots=g2LpQdwew
1&sig=PewjslX_Cry3N2CBAy2MWC7g9ys#v=onepage&q&f=false
Rueda Valdivia, F. J. (2005). Autodepuración y vertidos en cursos fluviales. Universidad de
Granada. Recuperado el 15 de Octubre de 2017, de
http://www.ugr.es/~iagua/LICOM_archivos/Tema_AD2DOC.pdf
Santamaría Alzate, F. (2013). Análisis simplificado de oxígeno disuelto en el río Ubaté por el
modelo QUAL2K. Tesis de maestría, Escuela Colombiana de Ingeniería Julio Garavito,
Bogotá. Recuperado el 3 de Marzo de 2018, de
https://repositorio.escuelaing.edu.co/handle/001/120
Streeter, H., & Phelps, E. B. (1958). A Study of the Pollution and Natural Purification of the
Ohio River. (E. &. U.S. Department of Health, Ed.) Washington, D.C. Obtenido de
http://udspace.udel.edu/handle/19716/1590
Suárez, J. (Febrero de 2008). Calidad de aguas en ríos. Autodepuración. Recuperado el 15 de
Abril de 2018, de Ceres:
ftp://ceres.udc.es/master_en_ingenieria_del_agua/master%20antiguo_antes%20del%2020
12/Segundo_Curso/Modelos_de_Calidad_de_Aguas/material%202010-
2011/MCA_Tema_5__MODELOS_DE_CALIDAD_DE_AGUAS_EN_RIOS.pdf
 129

Terry Berro, C. (2007). Gestión de residuales líquidos desde la perspectiva del consumo
sustentable. Revista electrónica de la Agenda de Medio Ambiente(12), 1-11. Recuperado
el 10 de Febrero de 2018, de http://ama.redciencia.cu/articulos/12.02.pdf
Torrado Ascanio, D. (6 de Julio de 2018). Comunicación personal. (J. Pérez Rincón,
Entrevistador)
Torres Muñiz, A. E. (2013). Monitoreo de la comunidad de insectos para evaluar la calidad de
insectos del agua del río Pesquería, desde el municipio de García hasta Pesquería,
Nuevo León, México. Tesis de maestría, Universidad Autónoma de Nuevo León.
Obtenido de http://eprints.uanl.mx/7255/1/1080259443.pdf
Trapote Jaume, A. (2013). Infraestructuras Hidráulico-Sanitarias II. Saneamiento y drenaje
urbano.
Trento, A. (Enero de 2001). Reseña sobre modelos unidireccionales de calidad de aguas en rios.
Revista Brasilera de Recursos Hídricos, VI(4). doi:DOI: 10.21168/rbrh.v6n4.p175-183
UN Water. (24 de Noviembre de 2014). El Agua, Fuente de Vida. Recuperado el 07 de
Septiembre de 2017, de Departamento de Asuntos econòmicos y Sociales de Naciones
Unidas:
http://www.un.org/spanish/waterforlifedecade/water_and_sustainable_development.shtml
UNDP-SIWI Water Governance Facility, WIN, Cap-Net, & WaterNet. (Enero de 2009).
Training manual on water integrity. Obtenido de Water Governance Facility:
http://watergovernance.org/resources/training-manual-on-water-integrity/
United States Environmental Protection Agency (EPA). (2015). Watershed & Water Quality
Modeling Watershed & Water Quality Modeling. Obtenido de United States
Environmental Protection Agency: https://www.epa.gov/sites/production/files/2015-
10/documents/tech_center_fact_sheet.pdf
Universidad del Cauca. (s.f.). Verifacación y validación. Todo alumnos. Recuperado el 31 de
Octubre de 2017, de
univirtual.unicauca.edu.co/moodle/pluginfile.php/31641/mod.../0/.../tema6.pdf?...
Universidad Francisco de Paula Santander (Ocaña). (s.f.). Objetivos. Obtenido de
https://ufpso.edu.co/Objetivos
Universidade da Coruña. (2010-2011). Curso: Calidad de aguas en ríos. Autodepuración.
Obtenido de Ceres:
ftp://ceres.udc.es/Master_en_Ingenieria_del_Agua/master%20antiguo_antes%20del%202
012/Segundo_Curso/Modelos_de_Calidad_de_Aguas/
Vanessa Callejas, G., & Lopez, L. A. (2015). Determinación de alcalinidad del río Cali.
Academia.
Vilanova, R., Santìn, I., & Pedret, C. (2017). Control y operaciòn de estaciones depuradoras de
aguas residuales: modelo y simulaciòn. Revista Ibenoamericana de Automàtica e
 130

Informàtica Industrial, 14, 217-233. Recuperado el 3 de Diciembre de 2017, de

https://doi.org/10.1016/j.riai.2017.05.004
Walter J., W. J. (1979). Control de la calidad del agua. Procesos fisicoquímicos. Barcelona:
Editorial Reverté, S.A. Recuperado el 3 de Marzo de 2018, de
https://books.google.com.co/books?id=TLpzh5HQYvgC&printsec=frontcover&hl=es#v=
onepage&q&f=false
Weather Underground. (15 de Junio de 2018). Tiempo histórico en Ábrego. Obtenido de WU:
https://www.wunderground.com/weather/co/%C3%A1brego
World Economic Forum. (2014). Global Risks. Recuperado el 25 de Enero de 2018, de
http://www3.weforum.org/docs/WEF_GlobalRisks_Report_2014.pdf
World Health Organization. (s.f.). Organización Mundial de la Salud. Recuperado el 30 de
Enero de 2018, de WHO Web site: http://www.who.int/water_sanitation_health/dwq/es/
WWAP (Programa Mundial de Evaluación de los Recursos Hídricos de las Naciones Unidas).
(2017). Informe Mundial de las Naciones Unidas sobre el Desarrollo de los Recursos
Hídricos. París, UNESCO. Obtenido de
http://unesdoc.unesco.org/images/0024/002476/247647S.pdf
WWDR-UNESCO. (2015). Agua para un Mundo Sostenible. Datos y Cifras. Informe de las
Naciones Unidas sobre los recursos hídricos en el mundo. Recuperado el 7 de Marzo de
2018, de
http://www.unesco.org/new/fileadmin/MULTIMEDIA/HQ/SC/images/WWDR2015Facts
_Figures_SPA_web.pdf
Zamudio Rodriguez, C. (2012). Gobernabilidad sobre el recurso hídrico en Colombia: entre
avances y retos. Gestión y Ambiente, XV(3), 99-112. Recuperado el 2018 de Mayo de 17,
de https://revistas.unal.edu.co/index.php/gestion/article/view/36284/42930
 131

Capítulo 8 Apéndice

Apéndice A. Protocolo de monitoreo del río algodonal

Basado en el Protocolo para el monitoreo y seguimiento del agua (IDEAM, 2007).

REQUISITOS DEL MONITOREO.

El encargado debe tener conocimientos básicos de hidrometría, orientados a definir cómo

evolucionan las variables hidrológicas en función del tiempo y el espacio, teniendo en cuenta que

la información recolectada será la base técnica del estudio ambiental.

La información hidrológica se obtendrá de las estaciones de monitoreo ubicadas a lo

largo del cauce del río Algodonal.

Para la selección de las estaciones de monitoreo se tendrá en cuenta la geometría del

perfil de la sección de aforo. Buscando que la geometría sea constante, con el fin de que la

velocidad mantenga su comportamiento horizontal y en profundidad; si por el contrario la

geometría varía, la velocidad puede presentar alteraciones debido a cambios por sedimentación o

socavación del lecho y/o los taludes.

Instrumentos:

 Cartera de aforo.

 Planillero.

 Personal: dos técnicos (inspector y aforador).

 Para medir el nivel del agua se utilizará una estaca afora cada 10 cm o

jalones de colores vivos (rojo, blanco).

 Para la medición de las velocidades se empleará un correntómetro.

 Cuerda.
 132

 Cinta métrica.

 Cavas.

 Gel refrigerante.

 Marcador permanente.

 Etiquetas.

 Cinta adhesiva.

 Reactivos.

 Termómetro ambiental y para medir la temperatura del agua.

 6 Recipientes de plástico, tapa rosca de capacidad de 2 litros.

 6 Recipientes winkler.

 6 Recipientes de vidrio para análisis microbiológico.

 Papel Kraft o papel aluminio.

Método de aforo.

Aforo por vadeo: Se utiliza cuando la profundidad es menor de un metro (< 1 m) y la

velocidad de la corriente menor de un metro por segundo (< 1 m/s). Estas condiciones permiten

que los operarios y los equipos se metan al cauce con seguridad, garantizando de esta manera

que la medición se realice con comodidad y sin riesgo.

Medición del caudal: El caudal corresponde al volumen de agua que pasa en ese instante

por la sección aforada en una unidad de tiempo, y se expresa en litros por segundo (l/s) cuando

se manejan pequeñas magnitudes.

 133

Caudal es una corriente de agua en función del producto del área (A) de la sección de

aforo y la velocidad media del flujo (V).

𝑄 =𝑉∗𝐴

Para la medición del caudal en una corriente se han desarrollado diversos procedimientos

que se aplican según el tamaño del cauce, la magnitud del caudal, las características hidráulicas

del flujo, la necesidad de contar con datos inmediatos o a corto plazo, la viabilidad de visitar los

sitios de medición con mayor o menor frecuencia, dependiendo de las distancias y recursos

logísticos, entre otros.

Tiempo de viaje 𝒕̅: tiempo que tarda una sustancia disuelta en ser conducida por el

movimiento del medio y ser dispersada. Se puede estimar a partir de ensayos con trazadores o

aproximadamente como la razón entre la longitud del tramo 𝐿 y la velocidad media del agua 𝑣̅

𝐿
𝑡̅ =
𝑣̅

AFORO ÁREA- VELOCIDAD

Dado que el caudal es función del área de la sección y la velocidad media del flujo, este

procedimiento se basa en la determinación de estas variables. Este sistema de aforo es el de

mayor uso y requiere que el flujo tenga un comportamiento laminar y que las líneas de flujo sean

normales a la sección transversal de aforo.

La precisión de las mediciones del caudal depende en gran parte del número de verticales

que se tomen para la ejecución de las mediciones para el aforo, profundidad, velocidad, toma de

muestras de agua, etc.

 134

En las primeras mediciones, que se realizan en una nueva estación hidrométrica, se debe

aumentar el número de verticales, con el fi n de ir Instituto de Hidrología, Meteorología y

Estudios Ambientales 44 conociendo en detalle la distribución de profundidades y velocidades a

lo ancho de la sección de aforos, para posteriormente ajustar (reducir) el abscisado, si las

condiciones lo permiten.

El ancho de cada sección parcial se denomina ancho parcial y corresponde a la distancia

existente entre dos verticales de medición consecutivas. Normalmente son iguales a lo ancho del

cauce, sin embargo, si el fondo del cauce es irregular deberá reducirse en los sectores más

profundos para cumplir con la norma que establece que los caudales parciales deben ser

inferiores al 10% del caudal total.

Apéndice B. Etapas de la campaña de aforo y muestreo.

1. Reconocimiento de las estaciones de monitoreo

Verificar que la zona cumpla con los requisitos mencionados anteriormente.

2. Aforo
Se inicia con la ubicación de dos estacas a lado y lado del río Algodonal. Se ata la cuerda y esta se

afora dependiendo del ancho del río, se recomienda cada metro.

Se miden las alturas de cada sección del aforo, para crear el perfil transversal del río.

Luego se toman las mediciones de velocidad en medio de cada sección de aforo, y sumergido a la

mitad de su altura. En caso de tener más de un metro de profundidad se doman varias mediciones

en cada sección del aforo.

 135

Los cálculos de caudal y la ilustración del perfil transversal se realizaran en oficina a través de un

formato de Excel.

3. Estimación de los tiempos de viaje

El tiempo de viaje se debe calcular in situ, con la longitud del tramo 𝐿 (medida de la estación X1 a

la estación X2 previamente desde Google Earth) y la velocidad media del agua 𝑣̅ (determinada por

el promedio de las velocidades obtenidas).

𝐿 𝑚
𝑡̅ = → 𝑡̅ = →𝑠
𝑣̅ 𝑚⁄
̅̅̅̅̅
𝑠

4. Toma de muestras
Se toman la temperatura ambiental, la humedad y la temperatura del agua.

Se llena un recipiente de plástico de 2 litros en contra de la corriente.

Se llena un Winkler en contra de la corriente, se le agrega 30 gotas del reactivo #1 30 gotas del

reactivo #2, se mezcla por inversión. Luego se le agregan 30 gotas del reactivo #3, se mezcla por

inmersión y se envuelve en papel kraft.

Se llena un recipiente de vidrio en contra de la corriente para el análisis microbiológico.

 136

Se etiquetan y se llevan a refrigeración.

Apéndice C. Variables fisicoquímicas y microbiológicas y sus métodos de evaluación

Parámetro Método Referencia

Alcalinidad Titulación por técnica (Instituo Nacional de Salud INS,
volumétrica (Standard Methods) 2011)
2320 B
Conductividad Método electrométrico (Instituo Nacional de Salud INS,
(Standard Methods) 2510 B 2011)
DBO5 Método yodométrico (IDEAM, 2007)
DBO20 modificación de azida
DBO5 Filtrada
Detritos (POM) 𝑀𝑂𝑃 = 𝐷𝐵𝑂5 𝑇𝑜𝑡𝑎𝑙 − (Jiménez & Cortes, 2016, pág. 36)
𝐷𝐵𝑂𝑓𝑖𝑙𝑡𝑟𝑎𝑑𝑎
Fósforo inorgánico (SRP) Cloruro estannoso por (Instituo Nacional de Salud INS,
Fosfatos colorimetría (Standard Methods) 2011)
4500 D (RECINTO UNIVERSITARIO DE
MAYAGÜEZ, s.f.)
(EcuRed, 2018)
NO3-nitrógeno Espectrofotométrico (Standard (Instituo Nacional de Salud INS,
(Nitrógeno nitrato) Methods) 4500 NO3 B 2011)
NH4- nitrógeno Método fotométrico
Oxígeno disuelto Método Winkler (IDEAM, 2004)
Patógenos Coliformes totales Filtración por membrana
Coliformes fecales
pH Electrométrico (Standard (Instituo Nacional de Salud INS,
Methods) 4500 H+ B 2011)
Sólidos inorgánicos (SSI) Solidos suspendidos (Ramalho, 1990, pág. 256)
Inorgánicos en función de la
volatilidad
Temperatura del agua (T) Termómetro portátil con sonda Autores del proyecto.
flexible Hanna CheckTemp (HI
98509)

Apéndice D. Plan de muestreo para el análisis físico-químico y microbiológico.

Volumen mínimo Tipo de Tiempo maximo de

Parámetro Preservación
de la muestra recipiente almacenamiento
Alcalinidad Refrigerada 200 P, V 24 h/ 14 días
Conductividad Refrigerar 500 P, V 28 d / 28 días
DBO20 Refrigerada P, V
DBO5 Refrigerada 1000 P, V 6 h/ 48 h
DBO5 Final Refrigerada P, V
Analizar
Fosfatos 100 V(A) 48 h / N.S
inmediatamente.
 137

Analizar tan pronto

como
sea posible o
adicionar H2 48 h / 48 h (28 dias para
Nitratos S04 a pH< 2 ; 100 P, V muestras cloradas
Analizar tan pronto ninguna / 28 dias
como
sea posible o
refrigerar
Nitrógeno Amoniacal Refrigerada 500 Plástico 7 d / 28 dias
Oxígeno disuelto (In
Acidificación 300 Winkler 8h / 8 h
situ)
Frasco de
Patógenos Refrigerada 200 6h
vidrio estéril
Analizar
pH P, V 2 h / N.S
inmediatamente.
Sólidos Suspendidos
Refrigerar P, V 7 d/2 - 7 dias
inorgánicos
Analizar
Temperatura P, V stat / stat
inmediatamente

Apéndice E. Recursos humanos.

El ingeniero ambiental Wilson Angarita Castilla, especialista en práctica docente

universitaria, quien dirigió el desarrollo de la investigación como director.
La ingeniera ambiental Rocío Andrea Miranda Sanguino, especialista en práctica docente
universitaria, quién fue colaboradora en el proceso de investigación.
Los estudiantes de pregrado Heiner Crisanto Vergel Cañizares, Diego Andrés Díaz
Jiménez, Arley Arenas Sánchez, María Fernanda Mora Sanjuán, Camila Fernanda Sánchez
Hamilton, estudiantes de la Universidad Francisco de Paula Santander Ocaña, del plan de
estudios de ingeniería ambiental.
El grupo de investigación MINDALA, quienes apoyaron el monitoreo y toma de
muestras en la primera campaña.
Apéndice F. Recursos institucionales.

Los datos hidrometeorológicos fueron obtenidos del Instituto de Hidrología,

Meteorología y Estudios Ambientales (IDEAM).
Información de concesiones y permisos de vertimientos relacionados con el río algodonal
se obtuvieron de la Corporación Autónoma Regional de la frontera Nororiental (CORPONOR).
Datos demográficos fueron obtenidos del Departamento Administrativo Nacional de
estadísticas (DANE).
 138

Información acerca del plan básico de ordenamiento territorial (PBOT), planes de gestión
y planes de desarrollo se adquirieron a través de la Alcaldía Municipal de Ocaña y la Alcaldía
Municipal de Ábrego.
Apéndice G. Recursos financieros.

FUENTE
AUTORES DEL
RUBRO UFPSO EXTERNA TOTAL
PROYECTO
EFECTIVO ESPECIE EFECTIVO ESPECIE EFECTIVO ESPECIE
Implementos
200000 200000
de muestreo
Salidas de
150000 150000
campo
Análisis de
calidad del 1120000 1120000
agua
Fotografías
80000 80000
aéreas
TOTAL 1550000
 139

Apéndice H. Trabajo en campo.

 140

Apéndice I. Ficha técnica y etiqueta.

CAMPAÑA FECHA:
ESTACIÓN: CORRIENTE: río Algodonal COORDENADAS ALTURA

HORA INICIAL: HORA FINAL:

TEMPERATURA AMBIENTE: TEMPERATURA DEL AGUA: HUMEDAD: OBSERVACIONES:

Sección Velocidad media Profundidad (m) Ancho Ancho Área Caudal

del caudal (m/s) (m) acumulado (m) (m2) (m3/s)

0
1
2
 141

3
4
5
6
7
8
9
10
11
12
13
14
Velocidad
promedio
Tiempo de viaje de la masa de agua
En segundos En minutos En horas Hora
 142

Apéndice J. Laboratorio.
143
144