María Fernanda QUAL2K

MODELACIÓN DE LA CALIDAD DEL AGUA DE LA
MICROCUENCA DEL RÍO BAHAMÓN CON EL PROGRAMA

DE LA AGENCIA DE PROTECCIÓN AMBIENTAL (EPA) QUAL2K
María Fernanda Gaitán Cárdenas
Informe de práctica empresarial para obtener el título de

Ingeniera Ambiental
Director
Ingeniero William Andrade
UNIVERSIDAD EL BOSQUE
FACULTAD DE INGENIERÍA AMBIENTAL
LÍNEA DE INVESTIGACIÓN: “MANEJO INTEGRADO DEL
RECURSO HÍDRICO”
BOGOTÁ D.C.
2006
1
“La Universidad El Bosque, no se hace responsable
de los conceptos emitidos por los investigadores en su
trabajo, solo velará por el rigor científico,
metodológico y ético del mismo en aras de la
búsqueda de la verdad y la justicia”
2
AGRADECIMIENTOS
Agradezco a Dios por darme lo que tengo y ser lo que soy, a la universidad El
Bosque, a la facultad de Ingeniería Ambiental, al profesor William Andrade por guiarme
y colaborarme con el proyecto, a mis compañeros durante esta etapa de mi vida y a mi
familia por brindarme no solo su apoyo incondicional sino las herramientas necesarias
para ser profesional.
A la familia Gaitán Medina por apoyarme y creer en mí siempre, a Patricia

Esteban por empujarme cuando lo necesite y a Jhon por estar cada día a mi
lado…GRACIAS POR TODO….
3
DEDICATORIA
Esto es para quienes me han dado todo…MIS PADRES…Carlos y Martha, para

mis hermanos quienes han estado siempre conmigo, para los amigos que tengo gracias a
estos años de carrera y para la persona que ha estado conmigo cada día durante estos
años…..JHON.
4
GUÍA DE CONTENIDO
Pág.
CAPÍTULO 1
INTRODUCCIÓN 13
1.1. Problemática 13
1.2. Justificación 14
1.3. Objetivos
1.3.1. General 15
1.3.2. Específicos 15
CAPÍTULO 2
DESCRIPCIÓN ZONA DE ESTUDIO
2.1. Localización 16
2.2. Topografía 19
2.3. Climatología 19
2.4. Hidrografía 20
CAPÍTULO 3
REVISIÓN BIBLIOGRÁFICA
3.1. Calidad de agua 23
3.2. Características del agua 23
3.2.1. Características físicas 23
3.2.2. Características químicas 24
3.2.3. Características biológicas 24
3.3. Modelación 25
3.4. Modelos de calidad de agua 26
3.4.1. Modelo QUAL2-E 26
3.4.2. Modelo QUAL-II 26
3.4.3. Modelo Quasar 27
3.4.4. Modelo ADZ 27
5
3.4.5. Modelo WASP 6 27
3.4.6. Modelo de dilución probalístico (PDM) 28
3.5. Modelo QUAL2K 28
3.6. Normatividad
3.6.1. Normatividad de agua en Colombia 34
3.6.2. Normatividad de agua según EPA 37
CAPÍTULO 4
METODOLOGÍA
4. Metodología
4.1. Revisión de literatura 39
4.2. Trabajo de campo 39
4.3. Calibración del modelo QUAL2K 45
4.4. Implementación del modelo QUAL2K 46
4.5. Análisis y conclusiones 50
4.6. Digitalización mapa final 51
4.7. Datos utilizados 51
CAPÍTULO 5
RESULTADOS

5.2. Calibración del modelo 67
5.3. Modelación final 70
CAPÍTULO 6
ANÁLISIS Y CONCLUSIONES 77
CAPÍTULO 7
RECOMENDACIONES 80
6
TERMINOLOGÍA BÁSICA 82
REFERENCIAS BIBLIOGRÁFICAS 85
ANEXOS 87
7
ANEXOS
Pág.
Anexo 1. Datos utilizados 87
Anexo 2. Modelo de Streeter – Phelps 91
Anexo 3. Resultados de Modelo de Streeter – Phelps 99
Anexo 4. Fotografías de la salida de campo 104
Anexo 5. Mapa de la microcuenca del río Bahamón con los puntos de muestreo 110
Anexo 6. Mapa de la contaminación de la microcuenca del río Bahamón 111
8
LISTADO DE TABLAS
Pág.
Tabla 1. Variables de estado “QUAL2K” 32
Tabla 2. Criterios calidad de agua Colombia 36
Tabla 3. Criterios calidad de agua EPA 38
Tabla 4. Georreferenciación puntos de muestreo 40
Tabla 5. Tipos de celdas “QUAL2K” 46
Tabla 6. Datos utilizados para la calibración 67
Tabla 7. Datos calibrados de temperatura 67
Tabla 8. Datos calibrados de oxígeno disuelto 68
Tabla 9. Datos calibrados de pH 69
Tabla 10. Datos calibrados de DBO 69
Tabla 11. Datos fisicoquímicos 88
Tabla 17. Datos de campo 90
Tabla 18. Balance de materia 94
Tabla 19. Variación de velocidades y profundidades usadas en las ecuaciones
12, 13 y 14 para flujos aireados a 20 °C 95
Tabla 20. Coeficiente de aireación de 1er orden para distintos ríos a 20 °C 96
Tabla 21. Coeficientes de descomposición a 1er orden para distintos ríos a 20 °C 96
Tabla 22. Factores Ө, para la obtención de las constantes cinéticas reales 97
Tabla 23. Tramos río Bahamón 100
Tabla 24. Resultados cálculos Streeter – Phelps 100
9
LISTADO DE FIGURAS
Pág.
Figura 1. Mapa Colombia/Cundinamarca 17
Figura 2. Mapa microcuenca del Río Bahamón 18
Figura 3. Diagrama de metodología 52
Figura 4. Hoja “QUAL2K” 53
Figura 5. Hoja “Headwater” 55
Figura 6. Hoja “Reach” 56
Figura 7. Hoja “Light” 57
Figura 8. Hoja “Meteorology” 57
Figura 9. Hoja “Point source” 58
Figura 10. Hoja “Rates” 61
Figura 11. Hoja “WQ Data” 61
Figura 12. Hoja “Temperature” 62
Figura 13. Gráfica QUAL2K “Temperatura” 63
Figura 14. Gráfica QUAL2K “Oxígeno disuelto” 64
Figura 15. Gráfica QUAL2K “pH” 65
Figura 16. Gráfica QUAL2K “DBO” 66
Figura 17. Gráfica calibración de temperatura 68
Figura 18. Gráfica calibración de OD 68
Figura 19. Gráfica calibración de pH 69
Figura 20. Gráfica calibración de DBO 70
Figura 21. Gráfica QUAL2K “pH” 71
Figura 22. Gráfica QUAL2K “OD” 72
Figura 23. Gráfica QUAL2K “Temperatura” 73
Figura 24. Gráfica QUAL2K “DBO” 74
Figura 25. Gráfica Streeter –Phelps “Oxígeno disuelto” 101
Figura 26. Gráfica Streeter – Phelps “DBO” 102
10
LISTADO DE FOTOGRAFÍAS
Pág.
Fotografía 1. Río Bahamón 41
Fotografía 2. Quebrada La Miquita 41
Fotografía 3. Quebrada Las Mercedes 42
Fotografía 5. Quebrada Ocotá 43
Fotografía 6. Quebrada Ocotá 43
Fotografía 7. Quebrada Tres Esquinas 44
Fotografía 10. Quebrada Tres Esquinas 105
Fotografía 11. Quebrada La Miquita 106
Fotografía 15. Zona aledaña a quebrada Ocotá 108
Fotografía16. Vivienda aledaña a quebrada Ocotá 108
Fotografía 17. Vivienda aledaña a quebrada Ocotá 109
Fotografía 18. Quebrada Las Mercedes 109
11
En este estudio se pretende aplicar un modelo matemático
para explorar la calidad del agua de un río de montaña. En
este se aplica un modelo matemático desarrollado por
Steven Chapra 2000 para la Agencia de Protección
Ambiental (EPA) y además se utiliza el modelo de las
ecuaciones de Streeter – Phelps para evaluar y comparar el
comportamiento del oxígeno disuelto y la demanda
biológica de oxígeno en un río Colombiano (río Bahamón) a
partir de datos obtenidos en campo.
Los resultados del análisis de la modelación muestran el

grado de contaminación en las diferentes partes del río,
exponiendo el punto de mayor contaminación del mismo, es
decir el punto con menor concentración de oxígeno disuelto
y de mayor concentración de DBO la cual se encuentra
ubicada en el nacimiento de la quebrada Ocotá.
12
CAPÍTULO 1
Introducción
La calidad del agua de las cuencas que recorren el país se ve cada vez más afectada por el
uso actual que se le da a este recurso y por el manejo de los residuos como aguas residuales
domésticas e industriales y residuos sólidos domésticos. Por esto es importante conocer la
calidad del agua de las cuencas y poder llegar a conclusiones que ayuden a dar soluciones
viables que disminuyan la posibilidad de contaminación del recurso.
La Universidad El Bosque creo la línea de investigación: “Manejo integrado del recurso

hídrico”, donde se esta estudiando en diferentes aspectos ambientales la microcuenca del
río Bahamón para con esto dar resultados a la población de los municipios aledaños al río,
con miras a que se preserve la cuenca y tomen medidas de planeación ambiental en torno a
esta.
Para poder estudiar la calidad del agua del río se utilizó el modelo de la agencia de
protección ambiental QUAL2K (2000) y el modelo de Streeter-Phelps, los cuales permiten
analizar las condiciones actuales de la microcuenca y conocer el estado de la calidad de
agua, para con esto poder llegar a planear estrategias y alternativas de saneamiento
ambiental que eviten riesgos a la salud pública de la zona.
1.1. Problemática
La microcuenca del río Bahamón esta ubicada en el suroeste del departamento de

Cundinamarca y en el noreste de la provincia del Tequendama limitando con los
municipios de Anolaima, Cachipay y las veredas La Florida y Peña Negra, sus afluentes
son las quebradas Ocotá, Las Mercedes, Tres Esquinas y La Miquita.
13
Según la ubicación de la microcuenca y de la actividad económica de la zona, la cual
corresponde al sector primario de la economía, es decir, al sector agrícola y al sector
pecuario, la problemática se presenta por el mal manejo que se da al agua.
A la microcuenca se vierten los desechos del municipio de Cachipay los cuales son
vertidos en la quebrada Ocotá, esta presenta la mayor contaminación de los sectores de
cuenca de la zona ya que soporta las aguas residuales domesticas del municipio, los
desechos del matadero y de la plaza de toros; en la quebrada La Miquita se presenta una
alta concentración de fosfatos ya que esta es la zona donde hay mayores cantidades de
cultivos para exportación, y además en la parte alta de la quebrada se encuentra una
bocatoma que abastece el municipio de Anolaima con el agua sin potabilizar; la quebrada
Tres Esquinas tiene altas concentraciones de coliformes fecales ya que recibe las aguas
residuales de la vereda Peña Negra.
Estas altas concentraciones de contaminantes en la microcuenca trae como consecuencia

la mala calidad de agua que afecta a la población de la zona, ya que esta agua contaminada
abastece el acueducto de Cayundá que provee a 350 usuarios con agua de La Miquita y el
acueducto de Peña Negra que toma el agua de la quebrada Ocotá que recibe contaminación
de La Miquita y Las Mercedes.
1.2. Justificación
El decreto 1594 de 1984, capítulo III, artículo 29, dice que los usos del agua son:
consumo humano y doméstico, preservación de flora y fauna, agrícola, pecuario, recreativo,
industrial y transporte, algunas de estas actividades son realizadas en los asentamientos
aledaños a la microcuenca del río Bahamón, pero los desechos que se producen de estas no
son manejados adecuadamente ya que son vertidos en la corriente sin ningún tratamiento.
En meses anteriores fue realizada una visita a la microcuenca donde se pudo observar
que la falta de un plan de manejo adecuado para la disposición final de los residuos sólidos
14
puede originar problemas de contaminación de las aguas ya que estos son arrojados a la
microcuenca sin trato alguno.
Por esta razón es importante conocer la calidad de agua realizando un modelamiento de

esta con el programa QUAL2K de la EPA, el cual permite pronosticar la calidad del agua
para planear estrategias y alternativas de saneamiento ambiental que eviten riesgos a la
salud pública de la zona.
1.3. Objetivos
1.3.1. General
Calibrar e implementar el modelo de calidad de agua de la EPA QUAL2K en la
microcuenca del río Bahamón.
1.3.2. Específicos
• Calibrar el modelo de calidad de agua QUAL2K.
• Realizar los estudios necesarios y requeridos por el modelo de calidad de
agua QUAL2K.
• Espacializar los datos reales obtenidos con el modelo en un mapa de
contaminación del cuerpo de agua de la microcuenca del río Bahamón.
15
CAPÍTULO 2
Descripción zona de estudio
2.1. Localización
La microcuenca del río Bahamón esta ubicada en el sureste del departamento de

Cundinamarca, limitando con los municipios de La mesa, Cachipay y Anolaima en la
provincia del Tequendama; para el estudio de calidad de agua las muestras que se
analizaron fueron realizadas en las veredas y corregimientos de La Florida, Peña Negra y
el municipio de Cachipay.
El municipio de Cachipay esta ubicado al occidente de Bogotá, limita al norte con el

municipio de Anolaima, al sur con La Mesa, al oriente con Zipacón y al occidente con
Quipile; se encuentra en la vertiente de la cordillera oriental, la cabecera municipal se
localiza a 04° 05’ 05’’de latitud norte y a los 74° 31’ 05’’ de longitud oeste de Greenwich,
su extensión territorial es de 8.417.781 m2. La vereda de Peña Negra esta ubicada al
suroeste del municipio de Cachipay, su extensión territorial es de 10.455.187 m2 y la vereda
de La Florida esta ubicada la norte del municipio, su extensión es de 4.517.536 m2. (EOT
Cachipay, 1999).
16
Figura 1. Mapa Colombia / Cundinamarca
17
Figura 2. Mapa Microcuenca Río Bahamón
18
2.2. Topografía
El municipio posee una topografía quebrada, entre los accidentes geográficos más
importantes de la zona de estudio, están los cerros de Tablanca en donde se encuentran las
principales formaciones hídricas como la microcuenca del río Bahamón. (EOT Cachipay,
1999).
2.3. Climatología
2.3.1. Precipitación
La precipitación promedia anual en la zona es de 1472.1 mm. con dos periodos secos en
los meses de enero y julio y dos periodos de alta precipitación en abril y octubre, para el
municipio de Cachipay la precipitación media anual es de 1492 m.m. (EOT Cachipay,
1999).
2.3.2. Temperatura y altura m.s.n.m

El municipio se encuentra en los pisos térmicos templado y subpáramo cuya temperatura
oscila entre 8° y 24°C, la media se conserva en 16.8 °C. Su altura sobre el mar oscila entre
650 y 2400 m.s.n.m, (EOT Cachipay, 1999).
2.3.3. Humedad relativa

Esta supera el 80%, los valores más bajos se presentan en los meses de julio y agosto.
(EOT Cachipay, 1999).
2.3.4. Evaporación
Su punto mínimo esta registrado con 49.3 m.m. en febrero y su punto máximo es 100.8
m.m. en el mes de mayo. (EOT Cachipay, 1999).
19
2.3.5. Velocidad del viento
La velocidad máxima alcanzada en el año es de 2.00m/seg., los vientos que predominan
en la zona son los de ladera, además las barreras orográficas que esta presenta no permiten
velocidades catastróficas y no admiten vientos foráneos de alto impacto. (EOT Cachipay,
1999).
2.3.6 Brillo solar

La zona presenta 64.9 horas luz/mes mínimo en el mes de noviembre y 109 horas
luz/mes como máximo en el mes de agosto. (EOT Cachipay, 1999).
2.4. Hidrografía
La ubicación del municipio sobre la cordillera oriental determina la hidrografía del

municipio, ya que esta es caracterizada por cadenas montañosas separadas por profundos
cañones de ríos y quebradas que forman pequeños valles. (EOT Cachipay, 1999).
El eje hidrográfico es conformado por el río Bahamón y la Quebrada Doña Juana, las
cuales forman pequeñas hoyas hidrográficas que resultan ser propicios para la agricultura y
la ganadería. Los ríos y quebradas que se hallan son de lecho profundo y en forma de V.
2.4.1. Microcuenca Río Bahamón

Esta microcuenca presenta una red dendrítica del orden 5. Se origina del río Cachipay,
quebrada Salitre, quebrada Malabrigo, quebrada Chay y quebrada La miquita. Las aguas
recorren el municipio de Cachipay de noreste a suroeste. Esta desemboca en el río Curí. Es
considerada una arteria fluvial en la zona. Nace en el Alto del Roble y tiene una longitud de
15.65 Km. (EOT Cachipay, 1999).
Beneficia las veredas de Tocarema, Petaluma, Cayundá, Mesitas de Santa Inés, Peña
Negra, Uchuta, entre otras; es utilizada en el acueducto de Cachipay, Anolaima y la Mesa,
además para consumo humano, animal y agrícola. (EOT Cachipay, 1999).
20
2.4.2. Río Cachipay
Este nace en el norte del municipio en la zona montañosa de Tablanca. Tiene una
extensión de 9.5 Km. Es utilizado para acueductos municipales, consumo humano, animal,
agrícola y regadíos de hortalizas y flores. (EOT Cachipay, 1999).
2.4.3. Quebrada el Salitre

Nace en el alto de la Pava en limites con el municipio de Zipacón, tiene una extensión de
3.5 Km. Vierte sus aguas al río Cachipay. (EOT Cachipay, 1999).
2.4.4. Quebrada Agua Regada

Nace en la vereda Agua Regada en el municipio de Zipacón, tiene una extensión de 2.5
Km. Es el límite entre Cachipay y Zipacón. (EOT Cachipay, 1999). (EOT Cachipay, 1999).
2.4.5. Quebrada Petaluma

Nace en la vereda Petaluma, tiene una longitud de 2.1 Km. y vierte sus aguas al río
Cachipay. (EOT Cachipay, 1999).
2.4.6. Quebrada Canta Gallo

Nace en Tres Esquinas en límites entre Anolaima y Cachipay. Tiene una longitud de 2
Km. Recoge aguas de las quebradas Bisalope, San Juanito y Los monos (Anolaima). (EOT
Cachipay, 1999).
2.4.7. Quebrada Chay

Nace entre la parte baja del municipio de Anolaima y el norte del municipio de Cachipay,
tiene una longitud de 3.2. Km., vierte sus aguas al río Cachipay. (EOT Cachipay, 1999).
2.4.8. Quebrada La Pedregosa

Nace en la vereda Naranjal, tiene una longitud de 1.9 Km., vierte sus aguas en el río
Bahamón. (EOT Cachipay, 1999).
21
2.4.9. Quebrada Ocotá
Nace cerca del casco urbano de Cachipay, su longitud es de 4.1 Km. y vierte sus aguas al
río Bahamón. (EOT Cachipay, 1999).
2.4.10. Quebrada Las Mercedes

Nace en el sector suburbano del municipio de Cachipay, tiene una longitud de 1.5 Km.,
vierte sus aguas a la quebrada los Micos. (EOT Cachipay, 1999).
2.4.11. Quebrada el Progreso

Nace en la vereda Petaluma, tiene una extensión de 1.7 Km., vierte sus aguas a la
quebrada El Hueso. (EOT Cachipay, 1999).
2.4.12. Quebrada los Micos

Nace de las quebradas el Hueso y el Progreso, tiene una longitud de 2.00 Km. y vierte
sus aguas al río Bahamón. (EOT Cachipay, 1999).
2.4.13. Quebrada Tres Esquinas

Nace en el cerro de Pan de Azúcar, tiene una longitud de 3.5 Km., vierte sus aguas al río
Bahamón. (EOT Cachipay, 1999).
2.4.14. Quebrada La Palestina

Nace en Peña Negra, en una longitud de 1.1 Km. y vierte sus aguas al río Bahamón.
22
CAPITULO 3
Revisión bibliográfica
3.1. Calidad de agua
La calidad del agua esta relacionada directamente con el manejo que se le da a esta, lo
que indica que es importante conocer las actividades que se le dan al agua para poder
establecer parámetros que se puedan cumplir con respecto a este tema. Si en dado caso el
agua no cumple con los parámetros establecidos se afirma que el agua esta contaminada y
que puede traer problemas si es usada en algunas actividades.
3.2. Características del agua
El agua tiene características físicas, químicas y biológicas que están ligadas con la
contaminación de esta ya que son las que muestran las variaciones causadas a las corrientes
de agua debido a vertimientos de residuos sólidos o líquidos sin ninguna clase de
tratamiento.
Las propiedades que están más relacionadas con las condiciones climáticas,
geomorfológicas y geoquímicas son las físicas y las químicas. Las propiedades químicas
varían con el clima, los usos de la tierra, entre otros, al alterarse estas el agua ya no es apta
para consumo humano. (Chapman, 1996).
3.2.1. Características físicas

• Temperatura
La temperatura esta relacionada con la actividad biológica, la adsorción de oxígeno y
demás procesos en un ecosistema acuático. Afecta todos los procesos que allí se dan, es
decir, la temperatura varía las concentraciones dependiendo de su aumento o
23
disminución, se influencia por la latitud y la longitud, época del año, hora del día,
nubosidad, caudal y profundidad del río.
3.2.2. Características químicas

• pH
Es el coeficiente que caracteriza la acidez o la basicidad de una sustancia; es el
cologaritmo decimal de la concentración en iones H+ presentes en el agua. En aguas
ácidas la concentración de iones hidronios aumenta, por lo tanto, el pH disminuye. El
pH es de gran importancia en la calidad de agua ya que muchos procesos no se pueden
llevar a cabo si el pH esta muy alto o muy bajo. (Sawyer, McCarty y Perkin, 2001).
3.2.3. Características biológicas

• Oxigeno disuelto
Es uno de los parámetros más relevantes en la calidad de agua ya que tiene relación
con la cantidad de materia orgánica presente en el agua. Mantiene los procesos que se
dan en los cuerpos de agua y además es un indicador de la destrucción de sustancias
orgánicas y del nivel de auto purificación del agua. Los niveles de oxigeno disuelto
dependen de la temperatura, salinidad, presión atmosférica y actividad fotosintética de
las plantas. La solubilidad del oxigeno disminuye cuando la temperatura y la salinidad
aumentan. Las concentraciones típicas al nivel del mar en el agua dulce son: 15 mg/l a
0°C – 8 mg/l a 25°C. El oxigeno disuelto en el agua se reduce con la respiración
biológica, vertimientos industriales y municipales con altos contenidos de materia
orgánica y nutrientes, aumento de la actividad microbiana y producción de biomasa por
medio de las plantas. Cuando las concentraciones son menores de 2 mg/l los peces
mueren. (Chapman, 1996).
• DBO
El DBO es importante para conocer el ingreso y uso del oxigeno en un cuerpo de
agua. Esta relacionado con la oxidación de la materia orgánica, es decir, la cantidad de
oxigeno que necesitan los microorganismos para oxidar la materia orgánica
convirtiéndola en materia estable.
24
Para la calidad de agua el DBO es esencial ya que ayuda a medir el grado de
contaminación, debido a que de forma indirecta mide la cantidad de materia orgánica
que se puede degradar biológicamente.
3.3. Modelación
Es una herramienta matemática utilizada para relacionar variables de calidad de agua y

puede ser usada para la predicción y planificación del saneamiento básico ambiental
(Chapra 1997). Es además una herramienta útil para la investigación de las corrientes de
agua.
A finales del siglo XIX, los modelos matemáticos fueron usados por los ingenieros
civiles para planear la creación de acueductos para agua de consumo humano;
posteriormente comenzaron a involucrarse en el desarrollo de plantas de tratamiento de
agua, redes de distribución y sistemas de recolección. El planteamiento de estos proyectos
fue directo por que sus metas estuvieron definidas, las cuales eran distribuir una cantidad
adecuada de agua potable para la ciudad y llevar de manera seguro a su consumo. En un
principio, descargaban su drenaje directamente a las aguas, lo cual fue observado como una
acción que iría transformando los ríos, lagos y estuarios, en grandes cloacas, en
consecuencia se empezaron a construir las plantas de tratamiento de agua, que manejaban
tasas de sedimentación simple y realizaban un tratamiento físico-químico, lo que podría
producir un efluente que fuera más original y que pudiera recibir agua, el problema era que
los costos eran muy elevados. Se observó entonces que ninguna de las opciones era viable
lo que llevo a replantear el proyecto y decidir que se debía proteger el ambiente de forma
adecuada pero económica. Entonces se decidió que el tratamiento de agua para consumo
humano se debería basar en la producción de un efluente que indujera a un nivel de calidad
aceptable en las reservas de agua. (Chapra, 1997).
Los modelos de calidad del agua ayudan a analizar los procesos físicos, químicos y
biológicos del agua de una corriente. Estos sirven como herramienta para determinar el
25
nivel de contaminación que presenta el agua de un río, el cual puede estar expuesto a
vertimientos industriales constantes.
La modelación busca simular los procesos físicos, químicos y biológicos referentes a la

calidad del agua. Además es un instrumento muy útil para determinar alternativas de
saneamiento y así poder implementar mejores prácticas en las industrias y sus vertimientos.
De esta forma se controlan los niveles de calidad físicos, químicos y biológicos del agua.
3.4. Modelos de calidad de agua
Existen diferentes modelos comerciales para la calidad del agua tales como QUAL2-E
(EPA), Modelo QUAL-II, QUASAR (Whitehead, et al., 1997), ADZ (Beer y Young,
1983), WASP 6 (Blanco, 2004), Modelo de Dilución Probabilística (PDM) y QUAL2-K
(Chapra, 2001).
3.4.1. Modelo QUAL2-E

El QUAL2E es un programa comprensivo y versátil que contiene un modelo de calidad
de agua fluvial, es una versión actualizada del modelo QUAL-II. Este modelo es
unidimensional, que puede usarse en condiciones de régimen permanente, como
impermanente. Permite, de acuerdo a su estructura, múltiples descargas, extracciones y
afluentes, considerados puntuales, como también caudales distribuidos. Puede simular,
numerosos parámetros de calidad de aguas, como Oxígeno Disuelto, Demanda Bioquímica
de Oxígeno, temperatura, Nitrógeno orgánico, amonio, nitrito, nitrato, Fósforo orgánico,
fósforo disuelto, clorofila a, coliformes etc. (Brown, 1987).
3.4.2. Modelo QUAL-II

Este modelo fue a su vez el modelo de Streeter-Phelps, es uno de los primeros modelos
para los cuerpos fluviales presentado por Streeter y Phelps en 1925. Permite predecir en el
caso de las descargas de material orgánico biodegradable su efecto sobre la concentración
de oxígeno disuelto. Dicho modelo considera el río como un reactor de flujo pistón, y
establece un balance para el oxígeno disuelto y el material biodegradable, bajo condiciones
26
de estado estacionario. En el caso del balance de oxígeno, se incluye la entrada de oxígeno
por transferencia desde el aire, también expresada por un modelo cinético de primer orden.
En versiones posteriores, se toma en cuenta, además, el aporte de oxígeno por vía
fotosintética. Por lo tanto la solución analítica del modelo de Streeter y Phelps para una
descarga puntual y continua, permite estimar la concentración de oxígeno disuelto a lo
largo de un río (Streeter & Phelps, 1925). Siendo este modelo el mejor adaptado para la
modelación de Oxígeno Disuelto y Demanda Bioquímica de Oxígeno, ya que no solo se
basa sobre la premisa estadística.
3.4.3. Modelo Quasar

Quality simulation along river systems (Blanco, 2004) es un modelo dinámico de la
calidad del agua de un río. Con este se modela nitrógeno, OD, temperatura, amoniaco, pH,
DBO y otros parámetros de calidad. El modelo involucra los efectos de la advección y
dispersión de la corriente en el transporte de la corriente, además del efecto de descargas
puntuales y difusas en el tiempo. El usuario tiene la posibilidad de seleccionar las
subrutinas cinéticas más representativas del proceso de contaminación de mayor interés
para su estudio
3.4.4. Modelo ADZ

El modelo ADZ es un modelo de transporte de solutos, que considera zonas muertas y
almacenamiento, explicando los procesos de advección y dispersión. El modelo toma el
tramo del río como un sistema imperfectamente mezclado donde el contaminante, este se ve
afectado por procesos de advección y dispersión.
3.4.5. Modelo WASP 6

El WASP (Water Quality Simulation Program, Di toro et al 1983, Ambrose R. B. et al,
1988), es un modelo que permite al usuario interpretar y predecir la respuesta de la calidad
del agua a un fenómeno natural o un evento de contaminación de origen antropogénico,
permitiendo analizar las efecto de las medidas y gestiones involucradas para controlar y
mitigar el impacto. El WASP es un compendio de modelación dinámico de sistemas
acuáticos, que involucra los efectos tanto sobre la columna de agua como sobre la zona
béntica.
27
Una importante ventaja del modelo es que el análisis dinámico lo convierte en una
importante herramienta para el seguimiento del ecosistema en diferentes épocas del año,
involucrando la variación de las condiciones climáticas de que caracterizan una
determinada temporada del año, (lluvia, verano, crecientes), además de permitir evaluar en
el tiempo la eficiencia de operación de cualquier medida de mitigación de la contaminación
en la corriente.
3.4.6. Modelo de Dilución Probabilística (PDM)
El modelo de Dilución Probabilística (PDM), originalmente desarrollado por Office of

Water, estima cuan frecuentemente se excede una concentración dada en un río o arroyo
que recibe contaminantes. Las estimaciones se basan en la variabilidad natural de los flujos
de los cursos de agua y en los flujos de efluentes, y en una ecuación de balance de masa y
dilución. Este modelo es una herramienta de evaluación de riesgo ecológico acuático usada
para predecir el número de días por año que una concentración química de interés será
excedida después de la descarga desde una instalación industrial (Río Damas, 1997).
El modelo QUAL2K fue desarrollado con el fin de mejorar las características de su

antecesor QUAL2-E. Este no tiene los limitantes del modelo Q2E como lo son un número
limitado de tramos y su interfase poco amigable. Este modelo presenta características útiles
para los ríos de montaña como son permitir la prescripción de las constantes de reaireción
por el usuario o calculo interno bajo las condiciones hidráulicas, modelación de
fitoplancton, etc. Este modelo fue desarrollado por el Profesor Steven Chapra, calcula
endógenamente la degradación de la materia orgánica en la metanogénesis bajo condiciones
anaeróbicas y representa adecuadamente interacciones agua sedimentos.
1
CHAPRA, PELLETIER G.J. 2003. QUAL2K: A modeling framework for simulating river and stream water
quality: Documentation and user manual.
28
El modelo de la EPA QUAL2K (o Q2K) es un modelo de la calidad del agua del río y de
la corriente que se piensa para representar una versión modernizada del modelo de QUALÈ
(o QÈ) (Marrón y Barnwell 1987). Q2K es similar a QÈ en los respectos siguientes:
• Una dimensión. El canal es bien mezclado vertical y lateralmente. Los ríos se

identifican por tener una corriente unidireccional con una altura relativa, el flujo es
muy variable en el tiempo, dependiendo de los cambios climáticos y patrones de
drenaje. En general los ríos alcanzan una mezcla vertical continua debido a la
corriente y turbulencia presentes en estos.
• Hidráulica. Simula flujo no uniforme, constante.
• La temperatura se simula en función de meteorología en escala de tiempo diurno.
• Todas las variables de calidad de agua se simulan en escala de tiempo diurno.
• Para las entradas de calor y masa se simulan las cargas en el punto de la descarga y
en el punto del río.
Este modelo incluye:
• Es funcional en hoja de cálculo de excel.

• El modelo QUAL2K tiene dos versiones una esta en estado estable y la siguiente
versión es en estado dinámico. Para nuestro estudio se utilizará el modelo en estado
estable ya que los ríos colombianos por ser ríos de montaña que están ubicados en
el trópico que no requieren de un estudio en las diferentes épocas para conocerse el
comportamiento de los ecosistemas se consideran flujos estables; además de
permitir evaluar en el tiempo la eficiencia de operación de cualquier medida de
mitigación de la contaminación en la corriente.
• El QUAL2K en el sistema de segmentos en los tramos del río abarca igualmente
los elementos espaciados. En contraste Q2K usa alcances espaciados desiguales.
Además, los cargamentos múltiples y las abstracciones se pueden entrar a cualquier
tramo.
• Q2K usa 2 formas de DBO que representan el carbono orgánico, estas formas son
una oxidación lenta y rápida; además se simula la materia orgánica de partículas no
29
vivas (Detritos). Este material detrítico es compuesto por partículas de carbono,
nitrógeno y fósforo en una estequiometría fija.
• Anoxia. Q2K acomoda la anoxia reduciendo las reacciones de oxidación a cero en
los niveles bajos de oxigeno, además la desnitrificación es modelada como una
reacción de primer orden que llega a ser pronunciado en las concentraciones de
oxigeno bajas.
• Interacciones sedimentos-agua. Los flujos del sedimento-agua y de los nutrientes
disueltos se simulan internamente en vez de prescribirse, es decir, los flujos de
oxigeno y nutrientes se simulan en función de fijar materia orgánica de reacciones
dentro de los sedimentos y las concentraciones de formas solubles en las aguas
superficiales.
• Algas inferiores. El modelo simula explícitamente algas inferiores unidas.
• Extensión ligera. Se calcula en función de algas, de detritos y de sólidos orgánicos.
• pH. La alcalinidad y el carbono inorgánico total se simularon. El pH del río se
simula basado en las 2 cantidades.
• Patógenos. Se simula un patógeno genérico. El retiro del patógeno se determina en
función de la temperatura, la luz y la posición.
El modelo de calidad de agua Q2K maneja matemática, las ecuaciones utilizadas

describen las cualidades del mismo y además permite ver como el modelo realiza las
simulaciones del cuerpo del agua, a continuación se muestran algunas de las ecuaciones:
Segmentación e Hidráulica:
• Balance de flujo:
Qi = Qi − 1 + Qin, i − Qab, I (1)
Donde: Qi: salida del tramo i en el tramo i+1; Qi-1: entrada del tramo por aguas arriba;
Qin,i: entrada total del tramo del punto de nacimiento y no punto de nacimiento; Qab,i:
salida total del tramo debido a l punto y no punto de abstracción.
30
• Ecuación de Manning:
Q=
( So1 / 2 × Ac 5 / 3 )
(2)
n × P2/3
Donde: Q: caudal; So: pendiente longitudinal del fondo; Ac: área mojada; n: número de
Manning; P: perímetro mojado.
• Tiempo medio de viaje:
Tk = Vk ÷ Ok (3)
Donde: Tk: tiempo medio de viaje; Vk: velocidad media; Qk: caudal.
Modelo temperatura
• Flujo superficial de calor:
Jh = I (0 ) + Jan − Jbr − Jc − Je (4)
Donde: J: flux; I(0): radiación solar neta de onda corta en la superficie del agua; Jan:
radiación solar neta de onda larga, Jbr: onda larga reflejada en el agua; Jc: conducción;
Je: evaporación.
• Transferencia de calor sedimento – agua:
∞s
Js, i = ρs × Cps × (Tsi − Ti ) (5)
Hsed / 2
Donde: ∞s: difusión térmica del sedimento; Tsi: temperatura de sedimentos del fondo
debajo del tramo; Cps: calor especifico; Js,i: flujo de calor sedimento – agua; ρs:
densidad de sedimentos; Hsed: volumen de la capa del sedimento.
31
El modelo utiliza los siguientes componentes o variables de estado, estas variables son los
parámetros que el modelo QUAL2K utiliza y con los que funciona para arrojar resultados
acerca de la calidad del agua:
VARIABLE SIMBOLO UNIDADES

Conductividad S µmhos
Sólidos suspendidos Mi mgD/l
inorgánicos
Oxigeno disuelto O mgO2/l
Reacción lenta CDBO Cs mgO2/l
Reacción rápida CDBO Cf mgO2/l
Nitrógeno orgánico No ugN/l
disuelto
Nitrógeno del amoniaco Na ugN/l
Nitrógeno del nitrato Nn ugN/l
Fósforo orgánico disuelto Po ugP/l
Fósforo inorgánico Pi ugP/l
Fitoplancton Ap ugA/l
Detritus Mo ugD/l
Patógenos X Cfu/100ml
Alcalinidad Alk mg CaCO3/l
Carbono inorgánico total Ct Mole/l
Algas inferiores ab gD/m2
Tomado de: Chapra, 2003
Tabla 1. Variables de estado “QUAL2K”
32
Además de las ecuaciones mostradas anteriormente maneja también las siguientes:
• Balance de masa:
dCi Qi −1 Qi Qab ,i E' E' Wi

= Ci −1 − Ci − Ci + i −1 (C i −1 − Ci ) + i (C i +1 − Ci ) + + Si (6)
dt Vi Vi Vi Vi Vi Vi
Donde: Wi: carga externa del componente para tramo i; Si: fuentes y sumideros del
componente debido a las reacciones y a los mecanismos de transferencia.
• Degradación lenta de carbono:
Scs = rro DetrDiss − SlowCHydr (7)
Donde: Scs: degradación lenta de carbono; rroDetrDiss: cociente de generación y

consumo de oxígeno en la disolución de detritos; SlowCHydr: reacción lenta de
hidrólisis;
• Degradación rápida de carbono:
Scf = SlowCHydr − FastCOxid − rondn Denitr (8)
Donde: Scf: reacción rápida de CBOD; SlowCHydr: hidrólisis lenta de carbono;

FastCOxid: oxidación rápida de carbono; rondnDenitr: cociente de oxígeno equivalente
perdido por el nitrato que es desnitrificado.
FastCOxid = Foxcf k dc (T )Cf (9)
Donde: Foxdf: atenuación debido a baja de oxígeno; Kdc(T)Cf: tarifa de temperatura

dependiente de oxidación rápida de DBO.
33
• Oxígeno disuelto:
So = roa PhytoGrowth + rod BotAgGrowth + roc FastCOxid − ron NH 4 Nitr

(10)
− ros PhytoRsp − rod BotAgRsp + OxRaer
Donde: So: oxígeno disuelto; roaPhytoGrowth: cociente de oxígeno debido al

crecimiento de el fitoplancton; rodBotAgGrowth: cociente de oxígeno debido a
crecimiento de algas inferiores; rocFastCOxid: cociente de oxidación rápida de oxígeno;
ronNH4Nitr: cociente de nitrificación NH4; rosPhytoRsp: cociente de oxígeno debido a la
respiración de fitoplancton; rodBotAgRsp: cociente de oxígeno debido a la respiración
de la algas inferiores.
OxRaer = Ka (T )(Os (T , elev ) − 0 ) (11)
Donde: Ka(T): coeficiente de oxígeno reaireación dependiente de la temperatura;

Os[T,elev]: concentración de saturación de oxígeno a temperatura T y elevación sobre
el nivel del mar.
3.6. Normatividad
3.6.1. Normatividad de agua en Colombia 2

En Colombia existen diversas leyes, decretos y resoluciones que indican la
normatividad adecuada para la calidad de agua de uso doméstico, humano, agrícola,
pecuario, recreativo, industrial y transporte. Es importante nombrar los artículos de la
legislación ambiental más relevantes que reglamentan la calidad del agua en nuestro
país, estos son:
2
Tomado de: http://www.ideam.gov.co:8080/legal/decretos.shtml. y
http://www.ideam.gov.co:8080/legal/resoluciones.shtml, el 17 de febrero de 2006.
34
• Decreto 1594 de 1984:
Articulo 24. Para el establecimiento de los modelos de simulación de calidad se
deberá realizar periódicamente a partir de la vigencia del presente decreto los análisis
pertinentes para obtener, por lo menos la siguiente información:
a. DBO
b. DQO
c. SS: sólidos suspendidos
d. pH
e. Temperatura
f. OD: oxigeno disuelto
g. Caudal
h. Datos hiodrobiológicos
i. Coliformes (NPM)
Articulo 29. Para los efectos del presente decreto se tendrán en cuenta los siguientes
usos del agua, sin que su enunciado indique orden de prioridad:
a. consumo humano y domestico

b. preservación de flora u fauna
c. agrícola
d. pecuario
e. recreativo
f. industrial
g. transporte
Articulo 39. Los criterios de calidad admisibles para la destinación del recurso para
consumo humano y doméstico son los que relacionan a continuación, e indican que para
su potabilización se requiere solo desinfección:
35
Nitritos 1.0
Nitratos 10.0
pH 6.5-8.5
Turbiedad UJT 10 unidades
Coliformes totales NPM 1000 microorganismos/100 ml
Artículo 38. Los criterios de calidad admisibles para la destinación del recurso para
consumo humano y doméstico son los que se relacionan a continuación, e indican que
para su potabilización se requiere solamente tratamiento convencional:
Referencia Expresado como Valor
Amoníaco N 1.0
Arsénico As 0.05
Bario Ba 1.0
Cadmio Cd 0.01
Cianuro Cn- 0.2
Cinc Zn 15.0
Cloruros Cl- 250.0
Cobre Cu 1.0
Compuestos
Fenólicos Fenol 0.002
Cromo Cr+6 0.05
Mercurio Hg 0.002
Nitratos N 10.0
Nitritos N 1.0
pH unidades 5.0-9.0
Plata Ag 0.05
Plomo Pb 0.05
=4
Sulfatos SO 400.0
Coliformes totales NPM 20.000 micro/100 ml
Coliformes fecales NPM 2.000 micro/ 100 ml
Tabla 2. Criterios calidad de agua Colombia
36
3.6.2. Normatividad de agua según EPA 3
Los criterios de calidad de agua, se han convertido en ítems importantes y algunas veces
controversiales en el campo del manejo adecuado del agua. Muchos de estos énfasis
contienen los descubrimientos de las últimas décadas vinculando los niveles de
contaminación en la alta incidencia de enfermedades como el cáncer.
Una breve revisión del proceso regulador destaca los mecanismos detrás de ella. Se
discuten los criterios, los estándares, y las metas de la calidad del agua aunque están
utilizados a menudo alternativamente en discurso diario. Todos cupieron bajo categoría
general de la calidad del agua "regulación". La regulación de la calidad del agua procede
típicamente en una manera de pasos aceptados lógica:
1. Se señalan las aplicaciones beneficiosas

2. Se desarrollan los criterios
3. Se promulgan los estándares
4. Se fijan las metas
Un sistema seleccionado de regulaciones de la EPA para el agua potable, representando

una gama de componentes con salud o impactos estéticos, se demuestra en la tabla
siguiente:
3
Montgomery, James m. & consulting engineers, Inc. 1985. Water treatment, principles and design, John
Wiley & sons, Inc. 67-77 pág.
37
Regulaciones internacionales seleccionadas del agua potable.
Parámetro Unidades U.S. EPA
Recomendado Nivel máximo
contaminación
Arsénico mg/l - 0.05
Cloruro mg/l 250 -
Cobre mg/l 1 -
Fluoruro mg/l - 1.4-2.4
Hierro mg/l 0.3 -
Plomo mg/l - 0.05
Manganeso mg/l 0.05 -
Mercurio mg/l - 0.002
Nitrato mg/l-N - 10
Sodio mg/l 20 -
Sulfatos mg/l 250 -
Zinc mg/l 5 -
TDS mg/l 500 -
Color mg pt/l 15 -
Turbides TU - 1
pH 6.5-8.5 -
Coliformes org/100ml - 1
Elementos espumo. mg/l 0.5 -
Tabla 3. Criterios calidad de agua EPA
38
CAPÍTULO 4
Metodología
4.1. Revisión de literatura
La revisión de literatura es el primer y más importante paso a seguir cuando se desarrolla

un proyecto ya que este brinda la información necesaria para enfocar y guiar a la
realización óptima de este, lo que traerá como consecuencia resultados viables a la solución
del problema.
Se realizó una búsqueda de proyectos similares realizados en Colombia para su lectura y

revisión, luego de tener un conocimiento acerca del tema se comenzó la lectura del manual
del programa de modelación de calidad de agua QUAL2K, para luego revisar la
normatividad que este utiliza como base (normatividad agua potable EPA), y la
normatividad Colombiana que es la que rige a la microcuenca del río.
De esta revisión surgió la redacción del marco conceptual, la estandarización e

implementación del modelo Q2K.
4.2. Trabajo de campo
La visita a la zona de estudio tiene como primer objetivo hacer un reconocimiento de la

microcuenca para luego corregir la georreferenciación de los puntos de muestreo y además
para tomar mediciones de pH, temperatura y oxígeno disuelto con ayuda del
multiparámetro, el oxímetro y el GPS.
39
Esto con el fin de obtener información acerca de los parámetros y poder ingresar los
datos correctos a la estandarización e implementación del modelo, para que este arroje
resultados confiables acerca de los cuatro parámetros a analizar.
En la zona de estudio se tomaron siete puntos de muestreo, la georreferenciación de estos

se presenta a continuación:
Punto de toma Coordenadas Altura

Río Bahamón después de Quebrada X 954.300.27
Tres Esquinas Y 1.010.822.23 1265msnm
X 955.244.24
Quebrada tres Esquinas Y 1.011.442.55 1775 msnm
X 956.646.72
Quebrada Bahamón antes de Ocotá Y 1.012.386.52 1600 msnm
X 959.748.34
Quebrada Ocotá Y 1.014.409.32 1265 msnm
X 960.125.93
Quebrada Las Mercedes Y 1.015.083.59 1755 msnm
Río Bahamón antes de Quebrada La X 958.750.43
Miquita Y 1.015.056.62 1265 msnm
X 958.966.19
Quebrada La Miquita Y 1.016.054.53 775 msnm
Tabla 4. Georreferenciación puntos de muestreo
40
1. Río Bahamón antes de la Quebrada La Miquita
Fotografía 1. Río Bahamón
2. Quebrada La Miquita
Fotografía 2. Q. La Miquita
41
3. Quebrada Las Mercedes
Fotografía 3. Q. Las Mercedes
4. Río Bahamón después de Quebrada Ocotá
42
5. Quebrada Ocotá
Fotografía 5. Q. Ocotá.
Fotografía 6. Q. Ocotá
43
6. Quebrada Tres Esquinas
Fotografía 7. Q. Tres Esquinas
7. Río Bahamón después de Tres Esquinas
44
4.3. Calibración del modelo QUAL2K
Hay que aclarar que se debe realizar una modelación inicial antes de la calibración del
modelo, la cual puede tener datos no confiables a causa de la fluctuación que estos pueden
tener por lo que se considera prudente calibrar los resultados iniciales en la simulación en
busca de disminuir la incertidumbre de los datos para luego implementar la modelación
final.
Para la calibración del modelo Q2K se puede utilizar una metodología de análisis de
sensibilidad y cuantificación de la incertidumbre con ayuda del software MATLAB, el cual
ejecuta herramientas de análisis para generar constantes que permitan la implementación
del modelo definitivo que es en si el resultado de este proyecto.
Esta herramienta de estandarización genera gráficas que muestran los percentiles

escogidos como por ejemplo 5 y 95%, a partir de una matriz; además puede generar
diagramas de dispersión que muestran la condición de los datos, esto con el fin de observar
el comportamiento de cada parámetro y de estimar los coeficientes de correlación para
obtener el dato óptimo para cada parámetro.
Cabe anotar que el método sugerido para la estandarización es el análisis de simulaciones

de Montecarlo o rutina de Montecarlo, el cual no es pertinente usar en este caso ya que los
datos del análisis fisicoquímico que se van a utilizar son muy pocos (entre 2 y 3 por
parámetro), y para implementar este método se necesita una cantidad mayor de datos por
parámetro y mayor cantidad de campañas.
Por esta razón se utilizó el método de ajuste de las gráficas por medio del R2 en excel
para linealizar las ecuaciones y poder llevar a cabo los pasos siguientes de la metodología
del proyecto.
45
4.4. Implementación del modelo Q2K
Para la correcta implementación del modelo Q2K hay que tener en cuenta los colores que
se utilizan para identificar la información de entrada (usuario), o de salida (programa), las
cuales se muestran en la tabla siguiente:
Datos usuario (campo)

Datos usuario
Datos calculados por el modelo
Parámetros a calibrar
Tabla 5. Tipos de celdas QUAL2K.
La primera hoja del modelo se utiliza para ingresar la información general:
• Nombre del río

• Fecha de simulación
• Tiempo luz en el día
• Tiempo de cálculo (horas)
• Tiempo final: duración total del cálculo en días.
El manual tiene diferentes hojas de trabajo, las cuales necesitan información; a

continuación se describen brevemente:
• Hoja de trabajo cabecera:

Es usada para incorporar el flujo y la concentración para los límites del sistema.
Requiere de los siguientes datos:
o Flujo (m3/s)
o Altitud y latitud de la cabecera
o Elevación
46
o Temperatura
o Oxigeno
o pH
o alcalinidad
• Hoja de trabajo del tramo:

Se utiliza para incorporar la información de la cabecera del río y del tramo
medido:
o Nombre del segmento

o Número de tramos: el modelo numera automáticamente en forma
ascendente.
o Longitud del tramo: el modelo computa y exhibe la longitud de cada tramo.
o Latitud y longitud de cada tramo: el modelo computa y exhibe
automáticamente la latitud y la longitud de cada tramo.
o Localización: el usuario incorpora el km del río para cada tramo final.
o Elevación: el usuario incorpora la elevación en m.s.n.m.
o Latitud y longitud aguas arriba: el usuario debe ingresar en horas, minutos y
segundos o en grados decimales este dato.
o Modelo hidráulico: hay dos opciones para computar la velocidad y
profundidad del flujo; curvas de rango y fórmula de Manning, el usuario
decide cual de las dos incorporando los datos requeridos.
o Dispersión prescrita: el modelo automáticamente lo calcula.
o Altura del vertedero
o Reaireación prescrita
o Cobertura de algas inferiores: usuario incorpora el porcentaje que cree que
tiene el río.
o SOD prescrito: Q2K simula este dato a partir de la cantidad de biomasa de
detritos y fitoplancton.
o CH4 prescrito: Q2K simula este dato.
o NH4 prescrito: Q2K simula este dato.
47
o Fósforo inorgánico prescrito: el usuario puede o no ingresar este dato.
• Hoja de trabajo de luz:

El usuario puede ingresar datos de luz ya que en el modelo las reacciones se realizan
en el foto periodo de la zona lo que afecta las reacciones causadas por la luz, estos
valores están por defecto en la hoja de trabajo y se recomienda dejarlos así.
• Hoja de trabajo de meteorología:
o Temperatura del aire: etiquetas y distancias y temperatura del aire.

o Temperatura punto de condensación: identificadores del tramo y temperatura
punto de condensación.
o Velocidad del viento: identificadores del tramo y velocidad del viento.
o Cubierta de nubes: identificadores del tramo y cubierta de nube.
o Sombra.
• Hoja de trabajo punto de fuente

Esta hoja de trabajo se utiliza para incorporar la información relacionando los
puntos de fuentes del sistema.
o Nombre: el usuario especifico en la etiqueta para identificar la abstracción o

afluencia particular.
o Localización: el kilómetro donde el punto de la fuente que entra o sale del
río.
o Fuentes de flujo de entrada o de salida: una fuente puede ser un afluente
(cargamento) o efluente (abstracción), pero no ambos.
o Componentes: la temperatura, calidad del agua, concentraciones de
afluentes.
48
• Hoja de trabajo de fuentes difusas:
Esta hoja de trabajo se utiliza para incorporar información que relaciona las
fuentes difusas del sistema.
o Nombre: el usuario especifico en la etiqueta para identificar la abstracción o

afluencia particular.
o Localización: kilómetros aguas arriba y aguas abajo que exceden la fuente o
la abstracción difusa que entra o sale del río.
o Afluencias y salidas de la fuente.
o Componentes: la temperatura, calidad del agua, concentraciones de
afluentes.
• Hoja de trabajo de constantes:

Se utiliza para incorporar las constantes de reacción para todos los procesos que
se dan en el modelo.
• Hoja de trabajo datos hidráulicos

Se utiliza para incorporar los datos hidráulicos requeridos:
o Distancia: distancia en la cual se trazan los datos de hidráulica (km),

o Datos – Q: datos de flujo en m3/s.
o Datos – H: datos profundidad en m.
o Datos – V: datos velocidad en m/s.
o Datos tiempo – viaje: tiempo viaje en días.
• Hoja de trabajo de datos de temperatura

Se usa para incorporar datos de temperatura:
o Distancia: distancia en la cual se trazan los datos de temperatura (km).

o Dato de temperatura en °C
49
• Hoja de trabajo de datos de calidad de agua (WQ)
o Distancia: distancia en la cual se trazan los datos de calidad de agua (km).

o Componentes: concentraciones medias de calidad de agua
o Otras concentraciones y fluyes: algas inferiores, nitrógeno total, fósforo
total, sólidos suspendidos totales, NH3, porcentaje de saturación, SOD,
sedimentos amonio, sedimentos metano, sedimentos fósforo inorgánico,
BOD, carbono orgánico total.
Luego de haber incorporado los datos requeridos por el modelo y de haber completado en
su totalidad las hojas de trabajo del programa, las últimas hojas de trabajo son los datos
arrojados por este y los diagramas resultantes del mismo con los datos que se deben
analizar.
• Hoja de trabajo de salida

Estas son series de hojas de trabajo que las tablas actuales de salida generadas por
el modelo.
o Resumen de fuentes: esta hoja de trabajo resume la carga total para cada
tramo del modelo.
o Temperatura de salida: resume la temperatura de salida para cada tramo del
modelo.
o WQ de salida: resume las concentraciones medias de salida para cada tramo
del modelo.
o Fluxes de sedimentos: resume los fluxes de oxigeno y nutrientes entre el
agua y el compartimiento subyacente del sedimento para cada tramo del
modelo.
4.5. Análisis y conclusiones
Después de haber ingresado los datos y de haber culminado la modelación con el

programa de calidad de agua Q2K, se analizan los datos y diagramas que este arrojo
50
después de compararlos con los parámetros de calidad de agua de la EPA, y de esto surgen
los análisis y las conclusiones finales del proyecto.
4.6. Digitalización de mapa final
Luego de tener los resultados finales es importante espacializar el proyecto en un mapa

que presente las zonas de la microcuenca con los diferentes puntos de muestreo y las zonas
mas contaminadas. Esta actividad se realizará en las instalaciones de la universidad
haciendo uso de la mesa digitalizadora y del software ArcView.
4.7. Datos utilizados
Los datos utilizados para este proyecto fueron los análisis fisicoquímicos tomados en
meses anteriores en siete puntos de muestreo en seis fechas distintas de DBO, oxigeno
disuelto, y además de esto los datos recogidos en la salida de campo como pH y
temperatura del agua en los mismos puntos de muestreo. (Anexo 1)
51
Revisión bibliográfica
Proyectos similares Manual Normatividad

programa Q2K
EPA Colombia
Salida de campo
Georreferenciación de puntos Toma de datos fisicoquímicos
Modelación inicial
Calibración del modelo
Ajuste estadístico
Implementación del modelo
Incorporar datos generales Incorporar datos necesarios
Se obtiene
Datos finales Gráficas
Análisis y conclusiones Digitalización mapa final
Figura 3. Diagrama de metodología
52
CAPÍTULO 5
Resultados
5.1. Modelo QUAL2K

Se realizaron los pasos descritos anteriormente en la metodología donde se realizo la
modelación de los parámetros de pH, temperatura y oxígeno disuelto con el modelo. Se
ingresaron los datos a las hojas de cálculo donde van los datos requeridos por el modelo y
los que este arroja, así:
• Hoja QUAL2K: En esta se ingresan los datos básicos y generales de la cuenca como
su nombre y la fecha de la campaña a modelar.
QUAL2K
Stream Water Quality Model
Steve Chapra
System ID:
River Name Río Bahamón
Month 3
Day 18
Year 2006
Sunrise 6.06
pCO2 0.000355 atm
Julian day 77
Calculation:
User calc step 0.009 day
Final time 3 day
Program determined calc step 0.007813 day
Time of last calculation 0.04 seconds
Figura 4. “Hoja QUAL2K”
53
• Hoja de cálculo “Headwater”, en ella se introducen datos importantes de la cuenca
como la longitud y latitud, datos hidráulicos y algunos datos fisicoquímicos como
pH y temperatura de la cabecera del río:
QUAL2K
Río Bahamón (3/18/2006)
Headwater Data:
Constant
characteristics:
Location 10.85 km
Headwater Latitude 4 Degrees
44 Minutes
7 Seconds
Headwater Longitude 74 Degrees
27 Minutes
36 Seconds
Elevation 1265 m
Flow cms
Bottom Width m
Side slope m/m
Slope m/m
Manning n 0.045
Alpha-Velocity
Beta-Velocity
Alpha-Depth
Beta-Depth
Alpha-Area
54
Beta-Area
Prescribed Dispersion m2/d
Upstream Conditions:
Temp (C) 21.90 C
Cons 0.00 units
Inorganic Solids 0.00 mgD/L
Oxygen 7.48 mg/L
Slow C 0.00 mgC/L
Fast C 0.00 mgC/L
Norg 0.00 ugN/L
NH4-N 0.00 ugN/L
NO3-N 0.00 ugN/L
Porg 0.00 ugP/L
SRP 0.00 ugP/L
Phyto 0.00 ugA/L
POM 0.00 mgD/L
Alkalinity 0.00 mgCaCO3/L
pH 5.89
Figura 5. “Hoja Headwater”
• Hoja de calculo “Reach”: es una de las mas importantes del modelo ya que en esta
se ingresan las características de cada tramo de estudio. Se deben ingresar datos
reales de la zona de estudio sino se pueden dejar en blanco.
55
QUAL2K
Reach Data:
Reach
length Downstream
Downstream Number (km) Latitude Longitude
R. Bahamón antes de Q. La
Miquita 1.00 0.00 4.74 74.46
Q. La Miquita 2.00 1.38 4.75 74.45
Q. Las Mercedes 3.00 4.50 4.73 74.44
R. Bahamón antes de Q. Ocotá 4.00 5.30 4.74 74.45
Q. Ocotá 5.00 4.95 4.73 74.44
Q. 3 esquinas 6.00 8.18 4.70 74.49
R. Bahamón después de 3
Esquinas 7.00 1.42 4.70 74.49
Figura 6. “Hoja Reach”
• Hoja de cálculo “Light”: el modelo realiza las reacciones durante el foto periodo de
la zona, esta no es conveniente modificar ya que son valores que se tienen por
defecto.
56
QUAL2K
Light Data:
Parameter Value Unit

Photosynthetically Avail Rad 0.45
Background light extinction 0.5 /m
Linear Chl light extinction 0.0088 1/m-(ugA/L)
Nonlinear Chl light extinction 0.054 1/m-(ugA/L)^2/3
SS light extinction 0.052 1/m-(mgD/L)
POM light extinction 0.174 1/m-(mgD/L)
Figura 7. “Hoja Light”
• Hoja de cálculo “Meteorology”: en esta se agregan los datos del ambiente que rodea
la cuenca ya que estos son de importancia por el intercambio térmico entre agua y
aire que se presenta en este.
QUAL2K

Met Data:
Air Dew-Point
Reach Reach Kilometer Temp Temp
Downstream Number Start End Ta(C) Td(C)
Miquita 1.00 10.85 10.85 16.000 0.000
Q. La Miquita 2.00 10.85 12.23 14.800 0.000
Q. Las Mercedes 3.00 12.23 16.73 15.300 0.000
R. Bahamón antes de Q. Ocotá 4.00 16.73 22.02 14.400 0.000
Q. Ocotá 5.00 22.02 26.98 13.800 0.000
Q. 3 esquinas 6.00 26.98 35.15 15.700 0.000
Esquinas 7.00 35.15 36.58 16.200 0.000
Figura 8. “Hoja Meteorology”
57
• Hoja de cálculo “Point Source”: en esta se ingresan los datos de cada uno de los
puntos, se especifican parámetros de calidad de agua como DBO, DQO, OD,
temperatura y pH y además el caudal en cada uno de los puntos.
QUAL2K
Point Source Data:
Point Point
Outflow Inflow Temp
Location
Name (km) (m3/s) (m3/s) (C)
Miquita 10.850 0.65 0 20.3
Q. La Miquita 12.225 0 0.036 19.8
Q. Las Mercedes 16.725 0 0.014 21.9
R. Bahamón antes de Q. Ocotá 22.025 0.7 0 19
Q. Ocotá 26.975 0.04 19.9
Q. 3 esquinas 35.150 0 0.008 22
Esquinas 36.575 0.748 0 19.9
Figura 9. “Hoja Point Source”
• Hoja de cálculo “Rates”: estos datos están por defecto y si el usuario no los obtiene
exactos puede dejar estos pero no se pueden modificar erróneamente ya que esto
puede agrandar las posibilidades de errores en los resultados.
58
QUAL2K
Water Column Rates
Parameter Value Units

Inorganic suspended solids
Settling velocity 1 m/d
Stoichiometry:
Carbon 40 mgC
Nitrogen 7.2 mgN
Phosphorus 1 mgP
Dry weight 100 mgd
Chlorophyll 1 mgA
Oxygen
Temp Corr (Reaeration) 1.024
O2 for Carbon oxidation 2.67 gO/gC
O2 for Nitrification 4.57 gO/gN
Oxygen inhib C Decomp 0.5 mgO/L
Oxygen inhib Nitrif 0.5 mgO/L
Slow C:
Hydrolysis 0.3 /d
Temp Corr 1.07
Fast C:
Decomposition 0.1 /d
Temp Corr 1.047
Organic N:
Hydrolysis 0.3 /d
Temp Corr 1.07
59
Ammonium:
Nitrification 0.5 /d
Temp Corr 1.07
Nitrate
Denitrification 0 /d
Temp Corr 1.07
Oxygen level denitr 0.5 mgO/L
Organic P:
Hydrolysis 0.1 /d
Temp Corr 1.07
Floating Plants (Phytoplankton)
Max Growth 2 /d
Temp Corr 1.07
Respiration 0.1 /d
Temp Corr 1.07
Death 0.1 /d
Temp Corr 1.07
Nitr Half Sat Constant 15 ugN/L
Phos Half Sat Constant 3 ugP/L
Light Half Sat Constant 35 cal/cm2/d
Light Model
1(HalfSat);2(Smith);3(Steele) 1
Ammonia preference 15 ugN/L
Settling velocity 0.1 m/d
Stationary Plants
(Periphyton/Macrophytes)
Max Growth 15 gD/m2/d
Temp Corr 1.07
Respiration 0.1 /d
Temp Corr 1.07
60
Death 0.05 /d
Temp Corr 1.07
Nitr Half Sat Constant 15 ugN/L
Phos Half Sat Constant 2 ugP/L
Light Half Sat Constant 50 cal/cm2/d
Ammonia preference 15 ugN/L
POM:
Hydrolysis 0.1 /d
Temp Corr 1.07
Settling Velocity 2 m/d
Figura 10. “Hoja Rates”
• Hoja de cálculo “temperatura Data”: esta hoja es para ingresar los datos de cada
punto que se tomaron en el agua. Estos datos no son obligatorios pero es
recomendable agregarlos para que el modelo se ajuste mejor.
Distance Data
x(km) Temp-data
10.850 20.3
12.225 19.8
16.725 21.9
22.025 19
26.975 19.9
35.150 22
36.575 19.9
Figura 11. “Hoja Temperature”
• Hoja de cálculo “WQ Data”: en esta hoja se agregan los datos de calidad de agua
medidos en campo.
61
QUAL2K
Boulder Creek, CO (8/21/1987)
Constituent Summary
Data:
Distance (km) cond-data SS-data DO-data
10.85 5.5
12.23 5.55
16.73 3.73
22.02 6.97
26.98 0.67
35.15 4.09
36.58 5.68
Figura 12. “Hoja WQ data”
Después del ingreso de los datos en las hojas de cálculo presentadas anteriormente, se
obtienen hojas con los resultados y gráficas de cada uno de los parámetros que se
corrieron en la primera modelación:
62
25
20
Temperatura (°C)
15
10
0
0 5 10 15 20 25 30 35 40
Distancia (km)
Temp(C) Temp-data
Figura 13. Gráfica QUAL2K “Temperatura”
63
5
OD (mg/L)
0
0 5 10 15 20 25 30 35 40
Distancia (km)
DO(mg/L) DO-data Dos(mg/L)
Figura 14. Gráfica QUAL2K “Oxígeno Disuelto”
64
5
pH
0
0 5 10 15 20 25 30 35 40
Distancia (km)
pH-data pH
Figura 15. Gráfica QUAL2K “pH”
65
300
250
200
DBO (mg/l)
150
100
50
0
0 5 10 15 20 25 30 35 40
-50
DISTANCIA (km)
DBO-data DBO
Figura 16. Gráfica QUAL2K “DBO”
66
5.2. Calibración del modelo
Luego de tener la primera modelación es necesario realizar una linealización de

las ecuaciones con la hoja de cálculo excel y las ecuaciones de tendencia con sus
ecuaciones y R2 para luego poder realizar la modelación final como ya se explico en
la metodología en el capítulo anterior. Para esto se utilizaron los datos que se
presentan a continuación:
Temperatura OD DBO pH
20.3 5.99 1.00 1
19.8 6.21 1.00 1
21.9 4.77 38.00 13.8
19 6.59 64.00 1.33
19.9 1.24 98.00 137.5
22 5.2 10.00 1
19.9 6.9 1.00 1
Tabla 6. Datos utilizados para la calibración
A continuación se presentan las ecuaciones linealizadas en excel, las cuales

arrojan los datos a utilizar en la modelación final.
• Temperatura: para este parámetro fue necesario realizar la calibración con la

línea de tendencia lineal ya que los datos son cercanos y la fluctuación no es
considerable.
Distancia Temperatura
10.850 20.73
12.225 20.7
16.725 20.8
22.025 20.55
26.975 20.48
35.150 20.32
36.575 20.31
Tabla 7. Datos calibrados de temperatura
67
y = -0.018x + 20.967
Río Bahamón R2 = 0.9031
20.9
20.8
Temperatura (°C)
20.7
20.6
20.5
20.4
20.3
20.2
0.000 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
Distancia (m)
Figura 17. Gráfica calibración de temperatura
• Oxígeno disuelto: en este caso se utilizó la línea de tendencia polinomial ya que

se deben considerar los picos altos y bajos que este parámetro presento en sus
mediciones.
Distancia OD
10.850 6.7
12.225 6.54
16.725 6.3
22.025 6.2
26.975 5.75
35.150 6.53
36.575 6.9
Tabla 8. Datos calibrados de OD
y = 0.0047x 2 - 0.2242x + 8.6193

Río Bahamón
R2 = 0.8105
8.00
7.00
6.00
OD (mg/L)
5.00
4.00
3.00
2.00
1.00
0.00
0.000 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
Distancia (m)
Figura 18. Gráfica calibración de OD
68
• pH: para este parámetro se utilizó la línea de tendencia lineal ya que los datos no
presentan una fluctuación con picos altos o bajos que se deban considerar.
Distancia pH
10.850 6.9
12.225 6.5
16.725 6.2
22.025 5.9
26.975 5.74
35.150 5.59
36.575 5.55
Tabla 9. Datos calibrados de pH
y = -0.0454x + 7.095
Río Bahamón R2 = 0.8785
8.00
7.00
6.00
5.00
pH
4.00
3.00
2.00
1.00
0.00
0.000 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
Distancia (m)
Figura 19. Gráfica calibración de pH
• DBO: para este parámetro se utilizó la línea de tendencia polinomial ya que su

forma presenta curvas pronunciadas y un pico alto.
Distancia DBO
10.850 1.00
12.225 1.00
16.725 38.00
22.025 64.00
26.975 98.00
35.150 10.00
36.575 1.00
Tabla 10. Datos calibrados de DBO
69
y = -0.5014x 2 + 24.249x - 215.29
Río Bahamón
R2 = 0.8758
120.00
100.00
80.00
DBO (mg/L)
60.00
40.00
20.00
0.00
-20.000.000 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
Distancia (m)
Figura 20. Gráfica calibración de DBO
5.3.Modelación final
Luego de realizar la calibración del modelo con el método de ajuste de gráficas

por medio del R2 en una hoja de cálculo, los datos que obtenidos se ingresan al
programa para la modelación final, lo que nos arroja como resultado las siguientes
gráficas:
70
5
pH
0
0 5 10 15 20 25 30 35 40
Distancia (km)
pH-data pH
Figura 21. Gráfica QUAL2K “pH”
71
5
OD (mg/L)
0
0 5 10 15 20 25 30 35 40
Distancia (km)
DO(mg/L) DO-data Dos(mg/L)
Figura 22. Gráfica QUAL2K “OD”
72
25
20
Temperatura (°C)
15
10
0
0 5 10 15 20 25 30 35 40
Distancia (km)
Temp(C) Temp-data
Figura 23. Gráfica QUAL2K “Temperatura”
73
120.00
100.00
DBO (mg/L)
80.00
60.00
40.00
20.00
0.00
0.000 5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000
Distancia (km)
DBO-data DBO
Figura 24. Gráfica QUAL2K “DBO”
74
Con estas gráficas arrojadas por el modelo se puede conocer el comportamiento de los
parámetros fisicoquímicos que se modelaron como lo son pH, OD y temperatura en el río
Bahamón en función de la distancia.
Además se conoce la capacidad de asimilación del río ocasionado por la carga producida
por las aguas residuales y las condiciones de calidad del agua del río constituyendo el
impacto negativo que tienen los vertimientos.
La gráfica de temperatura muestra que esta es normal y esta en los parámetros

establecidos para su calidad entre 15 – 20 °C, por lo que no es un factor que favorezca al
aumento de la contaminación en la microcuenca, es importante que este parámetro sea apto
ya que este tiene gran importancia en la supervivencia de organismos que requieren
condiciones específicas para vivir, además la temperatura influye directamente en el
comportamiento de parámetros como el pH, el déficit de oxígeno, la conductividad
eléctrica, entre otros.
La grafica de oxígeno disuelto muestra que no hay una homogeneidad en este dato en los
diferentes puntos de muestreo por lo que se puede asegurar que hay problemas de
autopurificación del agua en algunos puntos, pero esto no es de mayor importancia aguas
abajo ya que se presenta en los puntos de las quebradas y la dilución disminuye
considerablemente este impacto. Hay que mencionar que en el punto cinco, es decir, el
punto de la Quebrada Ocotá es donde hay menor presencia de oxígeno disuelto esto es
debido al alto grado de contaminación con residuos sólidos en el nacimiento de este cuerpo
de agua, esto es causado por el mal manejo de los residuos y por el vertimiento de aguas
residuales sin tratamiento, lo que puede traer como consecuencia problemas en la salud de
los pobladores de la zona, específicamente los habitantes aledaños a la quebrada en el
municipio de Cachipay.
La cantidad de oxígeno disuelto en el agua es de gran importancia ya que es esencial para

que los cuerpos de agua sean saludables y que pueda dar soporte a la vida animal y vegetal
75
ya que si los niveles de oxigeno disuelto en el agua son bajos algunas especies de peces y
otros organismos no pueden sobrevivir.
El pH de la microcuenca esta dentro de los parámetros establecidos por la ley

colombiana (Ley 1594 de 1984) que dice que debe estar entre 6.5-8.5 unidades, su valor
define la capacidad de autodepuración del río, además el contenido de materia orgánica y la
presencia de metales pesados que este pueda tener indicando contaminación.
En la microcuenca solo los puntos de la Quebrada Las Mercedes y del río Bahamón
después de la Quebrada Tres Esquinas están un poco por debajo de este parámetro sin ser
aguas lo suficientemente ácidas para traer problemas de consideración al cuerpo de agua en
ese punto, como la disminución de la vida acuática a causa de bajas en la tasa de
reproducción de los animales por la calidad del agua.
Para el DBO se observa claramente que el punto de mayor DBO es el punto de la

quebrada Ocotá, esto es complementario a la relación inversamente proporcional que existe
entre el oxígeno disuelto y el DBO, es decir, entre mayor sea el DBO menor será la
concentración de oxígeno disuelto ya que las bacterias que descomponen la materia
orgánica que se encuentra en el agua consumen en grandes cantidades y a gran velocidad el
oxígeno que esta presente en el agua disminuyendo los niveles de este hasta limitar la vida
acuática.
Se llegó a la conclusión de modelar el oxígeno disuelto y la demanda biológica de

oxígeno con el modelo de Streeter – Phelps, el cual brinda una opción más simple que se
podía llevar a cabo con la información obtenida, para luego con los resultados de este poder
realizar una comparación de los datos arrojados por los dos modelos en los parámetros
anteriormente mencionados. (Anexos 2 y 3).
76
CAPÍTULO 6
Análisis y conclusiones
• En este presente informe se presenta la calibración del modelo QUAL2K de

calidad de agua utilizando mediciones de campo. El procedimiento de
calibración utiliza la metodología de ajuste de ecuaciones y gráficas con el R2,
con el coeficiente de correlación. Los valores “óptimos” de los parámetros
calibrados se utilizaron para la modelación final. Hay que anotar que hacen falta
condiciones hidrológicas para que el modelo de calidad del agua sea validado y
las alternativas de saneamiento propuestas sean firmes.
• El oxígeno disuelto presenta un impacto significativo en el punto de la Quebrada

Ocotá ya que la cantidad de oxígeno es baja, esto puede ser causado por los
vertimientos de aguas residuales y residuos sólidos a esta sin ningún control ni
tratamiento. Entonces es necesario considerar alternativas novedosas de
saneamiento en la quebrada, por lo que se propone extender realizar de manera
más exhaustiva la modelación de la calidad del agua en cada quebrada del río.
• Según la gráfica no hay una homogeneidad en las cantidades de oxígeno disuelto

en los demás tramos del río lo que indica que no se esta presentando una
eficiente reaireación lo que disminuye la autopurificación del agua.
• Es de gran importancia tener condiciones aptas en la cabecera del río ya que

entre mejor sea el escenario inicial los vertimientos no afectaran de manera
considerable la dinámica de la calidad del agua, es decir, que si en las
condiciones iniciales del río ya hay problemas de contaminación se afectara el
ambiente del mismo aguas abajo.
• Si se tienen problemas de contaminación a lo largo de la longitud del río este

puede traer consecuencias en la salud de las poblaciones aledañas a la cuenca y
puede conllevar a problemas en las actividades económicas de la zona como la
77
agricultura y la ganadería, ya que las actividades económicas de la zona se
abastecen del agua del río Bahamón.
• Es necesario llevar a cabo un tratamiento de las aguas residuales domesticas y

del matadero que son vertidas al río sin ningún tratamiento, para con esto buscar
minimizar los impactos causados debido a la contaminación que se puede
presentar en la microcuenca.
• El pH del agua esta por dentro de los parámetros establecidos por la Agencia de
Protección Ambiental (EPA) y la normatividad Colombiana, esto indica que
teniendo en cuenta este parámetro es un agua apta para consumo humano,
siempre y cuando se realicen estudios microbiológicos.
• Es importante comparar las gráficas de oxígeno disuelto y el DBO para así

constatar que los datos arrojados por el modelo no son erróneos y muestran que
estos parámetros son inversamente proporcionales uno del otro, aunque esto sea
algo consecuente pero no suficiente para mostrar la relación existente entre estos
dos parámetros.
• Con el DBO se confirma que el punto más contaminado de la microcuenca del

río Bahamón, desde el punto de vista de la materia orgánica soluble y
biodegradable, es la quebrada Ocotá en su nacimiento ya que es la que mayor
concentración de este parámetro posee a la largo del cuerpo de agua.
• En el resto del río se presentan concentraciones bajas de DBO lo que indica poca
presencia de bacterias anaerobias descomponedoras de materia orgánica, lo que
a su vez muestra poca contaminación por vertimientos en estas zonas.
• Es importante concientizar a la población del municipio de Cachipay de la

importancia que tiene el preservar la microcuenca del río Bahamón, ya que esta
es la encargada de abastecerlos de agua en la zona tanto para su uso doméstico
como para sus actividades económicas.
78
• Se debe educar a la población en el tema de manejo adecuado de los residuos
sólidos para que se disminuya la cantidad de desechos arrojados al río
especialmente al nacimiento de la quebrada Ocotá, con talleres de educación
ambiental que les indiquen los lugares de disposición de los residuos ya
establecidos por el gobierno municipal de Cachipay.
79
CAPÍTULO 7
Recomendaciones
• Para realizar correctamente la modelación del río con el programa de la EPA

QUAL2K es importante tener en cuenta que se necesitan al menos dos campañas
de muestreo en cada época del año (seca y húmeda), además una campaña extra
para la calibración y el cálculo real de la reaireción.
• Es importante realizar más campañas de muestreo con el equipo humano
multidisciplinario necesario para que los datos obtenidos de la modelación sean
más confiables.
• Se debe tomar más datos para poder realizar una calibración estadística del
modelo con software como por ejemplo, Montecarlo en MATLAB o en Excel o
en C++, esto depende del investigador, para que los resultados sean más
confiables.
• Es necesario tener en cuenta el tiempo exacto de toma de cada muestra para
poder ingresar estos datos en el modelo y obtener las graficas que muestran el
tiempo de viaje, tiempo de muestreo y la exactitud de la modelación.
• Se deben hacer más estudios fisicoquímicos como nitritos, nitratos y fosfatos
para conocer en forma detallada la calidad del agua.
• Es necesario incluir dentro del presupuesto un aforo que permita conocer los
datos de caudal, profundidad, ancho del río y velocidad de este en los diferentes
puntos de muestreo.
• Es necesario para dar recomendaciones de saneamiento ambiental que se realice
adecuadamente la modelación de la microcuenca para que se pueda plantear la
posibilidad de una planta de tratamiento de agua sobre todo en la parte mas
contaminada del río que es la quebrada Ocotá.
• Se deben realizar campañas de educación ambiental que concienticen a la
población acerca de la importancia de la conservación de la microcuenca y
además para que conozcan el método del reciclaje para colaborar con el manejo
de los residuos sólidos del municipio.
80
• Realizar el Plan de Gestión Integral de Residuos Sólidos (PGIRS), para los
municipios aledaños a la microcuenca, para así mitigar el impacto causado en
esta por el vertimiento de los residuos sólidos domiciliarios.
• Se recomienda realizar estudios sectorizados de cada una de la quebradas de la

microcuenca con este modelo de calidad de agua QUAL2K.
81
TERMINOLOGÍA BÁSICA
Análisis de sensibilidad: este concierne el estudio de posibles cambios en la solución

óptima disponible como resultado de hacer cambios en el modelo original.
Calibración: es el procedimiento de comparación entre lo que indica un instrumento y

lo que "debiera indicar" de acuerdo a un patrón de referencia con valor conocido.
Cuenca: se denomina cuenca a un área física debidamente delimitada, en donde las

aguas superficiales y subterráneas vierten a una red natural mediante uno o varios
cauces de caudal continuo o intermitente, que confluyen a un cauce mayor que
desemboca o puede desembocar en un río principal, un depósito natural de aguas, un
pantano o directamente en el mar. (Código Nacional de Recursos Naturales. 1995)
Demanda Biológica de Oxígeno (DBO): es una medida de oxígenos que usan los
microorganismos para decomponer esta agua. Si hay una gran cantidad de desechos
orgánicos en el suministro de agua, también habrá muchas bacterias presentes
trabajando para descomponer este desecho. (Tomado de:
http://www.k12science.org/curriculum/depproj2/es/fieldbook/oxigeno.shtml. abril 8 de
2006).
Mesocuenca: es una subdivisión de la subcuenca, la cual se caracteriza por su eje axial

o cauce central que cae directamente al eje central de la subcuenca y su parte alta se
encuentra definida por su perímetro, filos o líneas divisorias de agua. (IGAC, 1999).
Microcuenca: es la unidad hidrográfica más pequeña correspondiente a los sectores

altos de las cuencas, está compuesta esencialmente por quebradas, riachuelos y
nacimientos, entre otros. (IGAC).
Modelación: Es una herramienta matemática utilizada para relacionar variables de

calidad de agua y puede ser usada para la predicción y planificación del saneamiento
básico ambiental (Chapra 1997).
82
Oxígeno Disuelto (OD): es la cantidad de oxígeno que está disuelta en el agua y que es
esencial para los riachuelos y lagos saludables. El nivel de oxígeno disuelto puede ser
un indicador de cuán contaminada está el agua y cuán bien puede dar soporte esta agua
a la vida vegetal y animal. Generalmente, un nivel más alto de oxígeno disuelto indica
agua de mejor calidad. Si los niveles de oxígeno disuelto son demasiado bajos, algunos
peces y otros organismos no pueden sobrevivir. (Tomado de:
http://www.k12science.org/curriculum/depproj2/es/fieldbook/oxigeno.shtml. Abril 8 de
2006).
Parámetros: Denominación de una cantidad o elementos de información que se utiliza

como valor de comparación en una rutina, subrutina, programa o cálculo matemático; se
le puede asignar un valor diferente cada vez que el proceso se repite. (Tomado de:
http://www.google.com.co/url?sa=X&start=4&oi=define&q=http://www.lorenzoservido
r.com.ar/info01/diccio-p-r.htm. Mayo 13 de 2006)
pH: es una expresión de la intensidad de las condiciones ácidas o básicas de un líquido

(o una suspensión, en el caso de los suelos); matemáticamente es el logaritmo en base
10 del recíproco de la concentración iónica de hidrógeno en moles por litro de solución,
y puede variar entre 0 y 14, donde 0 es el más ácido y 7 es neutro. Las aguas naturales
usualmente tienen un pH entre 6,5 y 8,5. (Tomado de:
http://www.ideam.gov.co/indicadores/calidad2.htm. Abril 5 de 2006)
Río de montaña: Un río de montaña o torrente, es un curso de aguas superficiales con

pendientes pronunciadas que presenta gran diferencia entre los caudales de crecida y de
estiaje. Durante la época húmeda conducen caudales grandes y en época seca los
caudales de escurrimiento son muy pequeños, llegando en muchos casos a no tener
escurrimiento superficial.
Tienen en común su gran pendiente y los importantes volúmenes de suelo que mueven,
principalmente en época húmeda. Durante los meses de estiaje el proceso de transporte
de sedimentos se detiene. (Tomado de:
http://www.eia.edu.co/sitios/webalumnos/laderas%20andinas/index.htm. Marzo 1 de
2006)
83
Simulación: es la experimentación con un modelo de una hipótesis de trabajo. La
experimentación puede ser un trabajo de campo o de laboratorio. El modelo de método
usado para la simulación seria teórico, conceptual o sistémico. (Tomado de:
http://es.wikipedia.org/wiki/Simulaci%C3%B3n. Mayo 13 de 2006).
Subcuenca: la primera división de las cuencas nacionales hidrográficas mayores se

denominan subcuenca. (IGAC, 1999).
Variables de estado: Entre las formas de modelar un sistema matemáticamente se
encuentra la de describir al sistema mediante la representación de variables de estado en
el modelo QUAL2K estas variables son los datos con los cuales el funciona. (Tomado
de: http://ib.cnea.gov.ar/~dsc/capitulo2/Capitulo2.htm. Mayo 13 de 2006).
Verificación: La constatación ocular o comprobación mediante muestreo, medición,

pruebas de laboratorio, o examen de documentos que se realiza para evaluar la
conformidad en un momento determinado. (Tomado
de:http://www.google.com.co/url?sa=X&start=2&oi=define&q=http://www.mejorambi
ente.com.mx/informacion/glosarios/glosario_turistico.htm. Mayo 13 de 2006).
84
REFERENCIAS BIBLIOGRAFICAS
• Alcaldía Municipal de Cachipay. Esquema de Ordenamiento Territorial, EOT

1999.
• Arenas F. Gustavo A, 2004. Modelación de la calidad del agua en un río de
montaña colombiano (Quebrada La Lejía). Universidad de los Andes, Facultad
de Ingeniería Ambiental. Departamento de Ingeniería Civil y Ambiental. Bogotá
Colombia.
• Borrero Mutis Santiago, et al. Estudio general de suelos y zonificación de tierras
del departamento de Cundinamarca. Instituto geográfico Agustín Codazzi.
Subdirección de agrología. Tomo 3. Bogotá 2000.
• Chapman D, 1996. “Water Quality Assessments, London & New York”.
• Chapra, S. C, 1997. Surface Water-Quality modelling, McGraw-Hill.
• Esteban mancera Patricia, 2005. Determinación experimental de parámetros de
transporte de solutos en la subcuenca del río Bahamón. Universidad El Bosque.
Facultad de Ingeniería Ambiental. Bogotá Colombia.
• http://www.car.gov.co/documentos/10_6_2004_11_46_05_AM_MapaRegionale
sCar.jpg. Tomado el 25 febrero de 2006.
• http://www.cra.gov.co/gc/www/section-10.jsp. Tomado el 20 febrero de 2006.
• http://www.ideam.gov.co:8080/legal/resoluciones.shtml. Tomado el 17 febrero
de 2006.
• http://www.ideam.gov.co:8080/legal/decretos.shtml. Tomado el 17 febrero de
2006.
• http://www.car.gov.co/documentos/10_6_2004_11_47_26_AM_Mapa_Politico_
Colombia.jpg. Tomado el 25 febrero de 2006.
• Jorge arboleda Valencia. Teoría practica de la purificación del agua. Tomo 1.
McGrawHill.
• L.A. Camacho B., J.A. Díaz C. & F.E. Díaz l. Modelación del transporte y de la
calidad del agua en la cuenca baja del río Bogotá. Universidad de los Andes.
Bogotá Colombia.
• L.A. Camacho B, Mario A. Díaz G. & Eugenio Giraldo. Contribución al
desarrollo de un modelo de calidad del agua apropiado para evaluar
85
alternativas de saneamiento ambiental en el río Bogotá. Departamento de
ingeniería civil y ambiental de la Universidad de los Andes, Bogotá Colombia.
• Montgomery, James m. & consulting engineers, Inc. 1985. Water treatment,
principles and design, John Wiley & sons, Inc. 67-77 pág.
• Oyarzún G. Jessica P, 2004. Modelación y simulación de oxígeno disuelto,

materia orgánica y relación distribucional de macroinvertebrados en la
subcuenca del río Traiguén. Universidad Católica de Temuco, Facultad de
ciencias. Temuco.
• Sabrina, Porras Lina, 2004. Modelación dinámica de la calidad del agua del río
Bogotá Villapinzón. Universidad de los Andes, Facultad de Ingeniería
Ambiental. Departamento de Ingeniería Civil y Ambiental. Bogotá Colombia.
86
ANEXO 1
Datos utilizados
87
Informe de laboratorio Daphnia Ltda.
Fecha: 20/08/05
PARAMETRO LUGAR
Bahamón antes
de Ocotá y Quebrada La Quebrada las
Quebrada Ocotá Miquita Miquita Mercedes
DBO5 120 mg O2/L < 1 mg O2/L <1 mg O2/L 12 mg O2/L
DQO 274 mg O2/L 5 mg O2/L 6 mg O2/L 51 mg O2/L
Oxigeno
disuelto 0,67 mg O2/L 5,50 mg O2/L 5,55 mg O2/L 3,73 mg O2/L
Tabla 11. Datos fisicoquímicos
Fecha: 17/19/05
PARAMETRO LUGAR
Bahamón Quebrada
antes de Tres Tres Quebrada La Quebrada las
esquinas Esquinas Miquita Mercedes Río Bahamón
DBO5 <1 mg O2/L <1 mg O2/L <1 mg O2/L 28 mg O2/L <1 mg O2/L
DQO 5 mg O2/L 19 mg O2/L 5 mg O2/L 56 mg O2/L 6 mg O2/L
Oxigeno
disuelto 5,68 mg O2/L 4,09 mg O2/L 5,04 mg O2/L 3,18 mg O2/L 5,57 mg O2/L
Fecha: 24/09/05
PARAMETRO LUGAR
Bahamón Quebrada
antes de Tres Tres Quebrada La Quebrada las
esquinas Esquinas Miquita Mercedes Río Bahamón
DBO5 1 mg O2/L 1 mg O2/L 1 mg O2/L 5 mg O2/L 1 mg O2/L
DQO 3 mg O2/L 25 mg O2/L 16 mg O2/L 44 mg O2/L 5 mg O2/L
Oxigeno
disuelto 8,29 mg O2/L 3,76 mg O2/L 5,66 mg O2/L 4,20 mg O2/L 8,67 mg O2/L
88
Fecha: 30/09/05
PARAMETRO LUGAR
Bahamón
Quebrada Quebrada antes de
Quebrada Tres Quebrada las Río Ocotá y
Ocotá Esquinas La Miquita Mercedes Bahamón Miquita
DBO5 166 mg O2/L 1 mg O2/L 1 mg O2/L 12 mg O2/L 1 mg O2/L 1 mg O2/L
13 mg
DQO 211 mg O2/L 5 mg O2/L 5 mg O2/L 23 mg O2/L O2/L 11 mg O2/L
Oxigeno 6,78 mg 7,35 mg 4,78 mg 6,43 mg 6,53 mg
disuelto 1,03 mg O2/L O2/L O2/L O2/L O2/L O2/L
Fecha: 9/10/05
PARAMETRO LUGAR
Bahamón Bahamón
antes de Quebrada Quebrada antes de Quebrada
Quebrada Ocotá y La las Tres Río Tres
Ocotá Miquita Miquita Mercedes Esquinas Bahamón Esquinas
202 mg < 1 mg <1 mg 12 mg < 1mg 2 mg <1 mg
DBO5 O2/L O2/L O2/L O2/L O2/L O2/L O2/L
364 mg 13 mg 17 mg 38 mg 6 mg 12 mg 13 mg
DQO O2/L O2/L O2/L O2/L O2/L O2/L O2/L
Oxigeno 0,68 mg 5,95 mg 7,48 mg 7,99 mg 8,33 mg 6,97 mg 6,12 mg
disuelto O2/L O2/L O2/L O2/L O2/L O2/L O2/L
89
Fecha: 17/10/05
PARAMETRO LUGAR
Bahamón antes
de Tres
Quebrada Ocotá Esquinas Tres Esquinas Río Bahamón
DBO5 62 mg O2/L < 1mg O2/L 1 mg O2/L <1 mg O2/L
DQO 174 mg O2/L 15 mg O2/L 18 mg O2/L 15 mg O2/L
Oxigeno disuelto 2,60 mg O2/L 5,31mg O2/L 5,26 mg O2/L 5,31 mg O2/L
Datos tomados en campo
Fecha: 18/03/06
Temperatura
NUMERO NOMBRE PUNTO pH (°C)
RIO BAHAMON EN PEÑA NEGRA DESPUES DE Q.

1 TRES ESQUINAS 7.85 19.9
2 QUEBRADA TRES ESQUINAS 5.9 22
3 RIO BAHAMON ANTES DE Q. OCOTA 6.1 20.3
4 QUEBRADA LA MIQUITA 6.7 19.8
5 RIO BAHAMON ANTES DE Q. LA MIQUITA 7.14 19
6 QUEBRADA OCOTA 6.8 19.9
7 QUEBRADA LAS MERCEDES 5.89 21.9
Tabla 17. Datos campo
90
ANEXO 2
Modelo de Streeter – Phelps
91
MODELO DE STREETER – PHELPS
Modelación de oxígeno disuelto.
En general, el balance de masa para una especie i considera los siguientes flujos:
(Entrada)i – (Salida)i – (Consumo)i + (Generación)i = (Acumulación)i
En el caso de un recipiente abierto a la atmósfera el balance de masa para la DBO y

Oxígeno disuelto se puede escribir como: (Chapra, 1997)
dDBO
V = kdVDBO (1)
dt
y
VdO = − kdVDBO + kaV (Os − O ) (2)
Ahora antes de proceder con la solución, se realiza una transformación que simplifica
el equilibrio del oxígeno. Para hacer esto se introduce una nueva variable,
D = Os − O (3)
Donde D es llamado “Déficit de Oxígeno Disuelto”, y la Ecuación 3 puede ser

diferenciada para dar,
dD dO
= (4)
dt dt
La Ecuación 3 y la 4 puede ser sustituida en la Ecuación 2 para dar,
dD
V = kdVDBO − kaVd (5)
dt
92
Así el uso del déficit simplifica la ecuación diferencial (Chapra, 1997). Supuestos,
según Zaror 1998:
1.- Contaminante sin difusión D = 0.
2.- Estado Estacionario dC/dt = 0.
3.- Decaimiento de 1er orden (Ecuación 1).
Si DBO = DBOo y D = 0 a t = 0, las ecuaciones 1 y 5 pueden ser resueltas por,
DBO = DBOoe− kdt (6)
D=
ka − kd
(
kdDBOo − kdt
e − e kat ) (7)
Ahora se puede pasar la teoría anterior al agua natural. Específicamente a la

modelación de un río con una sola fuente puntual y continua de DBO. Se caracteriza a
través de un reactor flujo pistón con hidrología y geometría constante. Ésta es la
manifestación más simple del modelo clásico de Streeter-Phelps 1925. Donde el balance
de masa se puede escribir como:
dDBO
−U − kdDBO = 0 (8)
dx
y
dU
−U + kdDBO − kaD = 0 (9)
dx
Si DBO = DBOo y D = Do a t = 0, estas ecuaciones pueden ser resueltas por (Chapra,

1997):
DBO = DBOoe− kdx / U (10)
93
(
D = Doe− kax / U + kdDBOo e− kdx / U − e− kax / U ) (11)
Ahora, se puede aplicar secuencialmente a una serie de tramos comprendidos entre

fuentes o cursos de agua tributarios, adoptando los valores de DBOo y Do adecuado
para cada tramo. En la siguiente tabla, se indican las condiciones aguas debajo de un
punto de vertido, en función de las condiciones de aguas arriba y de las características
del vertido.
AGUAS VERTIDO AGUAS ABAJO DEL

ARRIBA DEL VERTIDO
VERTIDO
CAUDAL Qx QD Q = Qx + QD
DBO DBOx DBOD DBOO = (QxDBOx +
QDDBOD)/Q
DEFICIT DE Dx DD DO = (QxDx + QDDD)/Q
OD
Fuente: Metcalf, 1995
Tabla 18. Balances de Materia
Constantes Cinéticas.
La estimación de la constante cinética de primer orden para la transferencia de

oxígeno desde el aire ka, se realiza a partir de varios modelos a 20°C:
O’Connor-Dobbins. (Schnoor et al., 1996); usada preferentemente para velocidades

moderadas y elevadas profundidades.
U 1.5
H = 3.93 (12)
H 1.5
Churchill. (Churchill et al., 1962); usada preferentemente para flujos muy rápidos.
U
Ka = 5.026 (13)
H 1.67
94
Owens and Gibbs. (Owens et al., 1964); usada preferentemente para velocidades
moderadas y sistemas más superficiales.
U 0.67
Ka = 5.32 (14)
H 1.85
O’CONNOR- CHURCHILL OWEN-GIBBS

PARAMETROS DOBBINS
PROFUNDIDAD 0.30-9.9 0.61-3.35 0.12-0.73
(m)
VELOCIDAD (m/s) 0.15-0.49 0.55-1.52 0.03-0.55
Fuente: Chapra, 1997.
Tabla 19. Variación de velocidades y profundidades usadas en las

ecuaciones 12, 13 y 14, para flujos aireados a 20°C.
La cuantificación de los valores de los parámetros del modelo se realiza durante una
calibración a priori, usando los valores de literatura (teóricos y experimentales), que
corresponden con la situación específica del sistema en estudio.
COEFICIENTES ka VALOR UNIDADES

Ríos lentos, profundos 0.1-0.4 1/d
Ríos condiciones típicas 0.4-1.5 1/d
Ríos veloces, profundos 1.5-4.0 1/d
Ríos veloces, poco 4.0-10.0 1/d
profundos
Fuente: Schnoor, 1996.
Tabla 20. Coeficientes de aireación de 1er orden, para distintos ríos a 20°C.
La estimación de la constante cinética de primer orden para el consumo de oxígeno por

acción metabólica kb, se realiza a partir de la tendencia expresada por Hydroscience,
1971 a 20°C (Chapra, 1997):
−0.434
⎛H⎞
Kd = 0.3⎜ ⎟ (15)
⎝8⎠
95
0 < H < 2.4 m
kd = 0.3 (16)
H > 2.4 m
COEFICIENTE kd VALOR UNIDADES

Sin previo tratamiento 0.35 (0.20-0.50) 1/d
Tratamiento primario 0.20 (0.10-0.30) 1/d
Lodo activado 0.075 (0.05-0.10) 1/d
Fuente: Chapra, 1997.
Tabla 21. Coeficientes de descomposición de 1er orden, para
distintos ríos a 20°C.
Efecto de la temperatura.
Mediante un balance de calor para la temperatura de mezcla, se obtiene la temperatura
mezcla del río.
T =
[(Qrío × T °río ) + (Q fp × T ° fp )] (17)
[Q río + Q fp ]
Tanto la transferencia de oxígeno como la tasa de eliminación de la DBO, son

influenciadas por la temperatura del medio. Esto se estudia con la relación de Van’t
Hoff-Arrhenius (Metcalf, 1995).
Ka = ka, 20°Cθ T − 20 (18)
y
Kd = kd , 20°C θ T − 20 (19)
96
VALOR REFERENCIA
Өa 1.005-1.030 Thomann and Mueller, 1987
1.024 Chapra, 1997
1.024 Schnoor, 1996
1.024-1.037 Schnoor, 1996
Өd 1.048 Schnoor, 1996
1.047 Schnoor, 1996
1.04 Thomann and Mueller, 1987
Tabla 22: Factores Ө, para la obtención de las constantes cinéticas reales.
Concentración de Saturación de Oxígeno Disuelto [O2]s en el río.
La [O2]s, se determina mediante un modelo analítico, en función de la temperatura

media del agua, obtenido por un polinomio de segundo orden (Chapra, 1996). Sin
embargo para comprobar se usan otras metodologías como la Ley de Henry o el uso de
una tabla que presente la concentración de OD en función de la temperatura, la salinidad
y de la presión barométrica (Metcalf, 1995).
Modelo Analítico (Chapra, 1996):
[O2 ]s = 0.0035T 2 − 0.3369T + 14.407 (20)
Concentración de Oxígeno Disuelto (OD) en el río a x (m) aguas abajo de la

descarga.
Se estima a partir de la diferencia entre la Concentración de Saturación del OD y el

Déficit de OD.
[O2 ] = [O2 ]s − [D ] (21)
97
ANEXO 3
Resultados Modelo Streeter – Phelps
98
RESULTADOS MODELO STREETER – PHELPS
Para el cálculo del modelo de Streeter – Phelps se requería de tramos que cumplieran
el criterio de aguas arriba o abajo de una quebrada que tenga vertimientos
contaminantes por lo que los siete puntos de muestreo que se tienen en el río se
agruparon como se muestra a continuación:
Tramo Río Quebrada
Río Bahamón antes de Q. La
1 Quebrada La Miquita
Miquita
Río Bahamón antes de Q. La
2 Quebrada Las Mercedes
Miquita
3 Río Bahamón antes de Q. Ocotá Quebrada Ocotá
Río Bahamón después de Q. Tres
4 Quebrada Tres Esquinas
Esquinas
Tabla 23. Tramos río Bahamón
Además el modelo requiere de datos de la zona que fueron aproximados a los reales y
de los datos utilizados para la modelación con el programa QUAL2K. Es decir que el
modelo se realizó con datos reales y teóricos.
A continuación se presenta una tabla de los resultados obtenidos a través de los
cálculos realizados con las ecuaciones presentadas anteriormente:
ka kd To Do DBOo DBO D [O2] [O2]s

(1/día) (1/día) (°C) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Tramo
1 8.06 0.34 19.92 -5.365 9.667 3.94 0.17 22.217 22.497
Tramo 1.1099 2.47x10-
2 1.81 0.33 19.06 -5.37 10.35 x 10-15 16 9.25 9.25
Tramo
3 0.55 0.35 20.27 -5.89 23.47 2.1 2.6265 6.38 9.01
Tramo -
4 2.96 0.35 19.9 6.1525 7.34 0.007 0.0009 9.08 9.08
Tabla 24. Resultados cálculos Streeter - Phelps
99
A partir de estos datos se realizaron gráficas de déficit de oxígeno disuelto,
concentración de oxígeno disuelto y decaimiento de DBO:
Oxígeno disuelto vs. distancia:
Oxígeno Disuelto
25
20
15
[O2]
10
0
0 2 4 6 8 10
Distancia
Figura 25. Gráfica Streeter – Phelps “Oxígeno disuelto”
Esta gráfica de oxígeno disuelto muestra que el tramo 1 es el que mayor cantidad de
oxígeno disuelto posee lo que indica buena calidad de agua, ya que hay una buena
autopurificación debido a la reaireación y que en esta zona el río posee vida animal y
vegetal. Además se puede decir que el agua de esta zona es óptima para consumo
humano y para su uso en las diferentes actividades económicas de la zona.
Al comparar esta gráfica con la obtenida con el modelo QUAL2K podemos ver que hay
una coincidencia en que el punto de menor cantidad de oxígeno disuelto es el de la
Quebrada Ocotá (6.38 mg/l) que es un lugar donde se vierten los desechos líquidos y
sólidos del municipio de Cachipay sin tratamiento alguno.
100
El agua de este punto posiblemente no sea apta para consumo humano especialmente en
el nacimiento de la quebrada Ocotá que es la zona más cercana al municipio y donde en
sus laderas hay habitantes, lo que ocasiona mayor problemas de contaminación al
recurso ya que a simple vista se observe el deterioro en esta zona, con los malos olores y
la cantidad de residuos sólidos que corren en el río.
En los tramos 2 y 4 la cantidad de oxígeno disuelto es mayor que en este punto con
valores de 9.25 y 9.08 mg/l respectivamente, indicando que hay problemas de
contaminación a lo largo del río debido a la mala reaireación que este tiene, pero que
esta por encima de lo mínimo necesario para que halla vida animal y vegetal en un río
que es de 5 mg/l.
Se observa además que el oxígeno disuelto disminuyo más del 50% de su saturación
entre los tramos 1 (22.497 mg/l) y 2 (9.25 mg/l), lo que confirma que en el tramo 2 ya
se presentan problemas de contaminación del agua.
DBO vs. distancia
DBO
25
20
DBO (mg/l)
15
10
0
0 2 4 6 8 10
Distancia
Figura 26. Gráfica Streeter – Phelps “DBO”
101
Dado que la demanda bioquímica de oxígeno es la cantidad de oxígeno disuelto
requerido por los microorganismos para la oxidación aerobia de la materia orgánica
biodegradable presente en el agua (Metcalf, 1995), la relación que existe entre estos dos
parámetros de la calidad del agua es inversa. Para este caso esta analogía se cumple
completamente ya que como se puede observar en la gráfica el punto de mayor DBO es
el de menor OD.
Esto es causado por la cantidad de microorganismos como las bacterias aeróbicas que
están presentes en el agua a causa de la contaminación, estas descomponen la materia
orgánica causando disminución del oxígeno disuelto en el agua, robándoles la
posibilidad a otros organismos acuáticos de vivir.
En el punto de la quebrada Ocotá hay una gran cantidad de desechos orgánicos lo que
trae como consecuencia también una gran cantidad de bacterias descomponiendo estos
desechos. Lo que indica una alta concentración de DBO que disminuirá con el consumo
o dispersión de los desechos.
El DBO y la temperatura están relacionadas ya que si uno de estos parámetros aumenta

el otro también, es decir, si el agua es tibia tendrá un nivel mas alto de DBO que el agua
fría, esto es debido a que cuando la temperatura del agua aumenta la velocidad de la
fotosíntesis también, lo que ocasiona el desarrollo más rápido y por consiguiente la
muerte más rápida de las plantas. Cuando esto sucede están son descompuestas en el
fondo por las bacterias que requieren de oxígeno para realizar esta labor, aumentando
las concentraciones de DBO en el cuerpo de agua.
102
ANEXO 4
Fotografías de la microcuenca del río Bahamón
103
Fotografía 10. Quebrada Tres Esquinas
104
Fotografía 11. Quebrada La Miquita
105
106
Fotografía 15. Zona aledaña a Quebrada Ocotá
Fotografía 16. Vivienda aledaña a Quebrada Ocotá
107
Fotografía 17. Vivienda aledaña a Quebrada Ocotá
Fotografía 18. Quebrada Las Mercedes
108
ANEXO 5
Mapa de microcuenca del río Bahamón con puntos de muestreo
109
110
ANEXO 6
Mapa de contaminación de la microcuenca del río Bahamón
111
112

María Fernanda QUAL2K

Cargado por

Copyright:

Formatos disponibles

María Fernanda QUAL2K

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

María Fernanda QUAL2K

Cargado por

Copyright:

Formatos disponibles

MODELACIÓN DE LA CALIDAD DEL AGUA DE LA

MICROCUENCA DEL RÍO BAHAMÓN CON EL PROGRAMA

María Fernanda Gaitán Cárdenas

Informe de práctica empresarial para obtener el título de

A la familia Gaitán Medina por apoyarme y creer en mí siempre, a Patricia

Esto es para quienes me han dado todo…MIS PADRES…Carlos y Martha, para

5.1. Modelo QUAL2K 53

Los resultados del análisis de la modelación muestran el

La Universidad El Bosque creo la línea de investigación: “Manejo integrado del recurso

La microcuenca del río Bahamón esta ubicada en el suroeste del departamento de

Estas altas concentraciones de contaminantes en la microcuenca trae como consecuencia

Por esta razón es importante conocer la calidad de agua realizando un modelamiento de

Descripción zona de estudio

La microcuenca del río Bahamón esta ubicada en el sureste del departamento de

El municipio de Cachipay esta ubicado al occidente de Bogotá, limita al norte con el

2.3.2. Temperatura y altura m.s.n.m

2.3.3. Humedad relativa

2.3.6 Brillo solar

La ubicación del municipio sobre la cordillera oriental determina la hidrografía del

2.4.1. Microcuenca Río Bahamón

2.4.3. Quebrada el Salitre

2.4.4. Quebrada Agua Regada

2.4.5. Quebrada Petaluma

2.4.6. Quebrada Canta Gallo

2.4.7. Quebrada Chay

2.4.8. Quebrada La Pedregosa

2.4.10. Quebrada Las Mercedes

2.4.11. Quebrada el Progreso

2.4.12. Quebrada los Micos

2.4.13. Quebrada Tres Esquinas

2.4.14. Quebrada La Palestina

3.1. Calidad de agua

3.2. Características del agua

3.2.1. Características físicas

3.2.2. Características químicas

3.2.3. Características biológicas

Es una herramienta matemática utilizada para relacionar variables de calidad de agua y

La modelación busca simular los procesos físicos, químicos y biológicos referentes a la

3.4. Modelos de calidad de agua

3.4.1. Modelo QUAL2-E

3.4.2. Modelo QUAL-II

3.4.3. Modelo Quasar

3.4.4. Modelo ADZ

3.4.5. Modelo WASP 6

3.4.6. Modelo de Dilución Probabilística (PDM)

El modelo de Dilución Probabilística (PDM), originalmente desarrollado por Office of

3.5. Modelo QUAL2K 1

El modelo QUAL2K fue desarrollado con el fin de mejorar las características de su

• Una dimensión. El canal es bien mezclado vertical y lateralmente. Los ríos se

Este modelo incluye:

• Es funcional en hoja de cálculo de excel.

El modelo de calidad de agua Q2K maneja matemática, las ecuaciones utilizadas

Qi = Qi − 1 + Qin, i − Qab, I (1)

• Tiempo medio de viaje:

• Flujo superficial de calor:

Jh = I (0 ) + Jan − Jbr − Jc − Je (4)

• Transferencia de calor sedimento – agua:

VARIABLE SIMBOLO UNIDADES

Tomado de: Chapra, 2003

Tabla 1. Variables de estado “QUAL2K”