CONTENIDO

Pág.

1. PREDICCIÓN A TRAVÉS DE MODELOS DE SIMULACIÓN DE LOS IMPACTOS

QUE CAUSE EL VERTIMIENTO AL CUERPO DE AGUA, EN FUNCIÓN DE LA CAPACIDAD
DE ASIMILACIÓN Y DILUCIÓN DEL CUERPO DE AGUA RECEPTOR Y DE LOS USOS Y
CRITERIOS DE CALIDAD ESTABLECIDOS EN EL PLAN DE ORDENAMIENTO DEL
RECURSO HÍDRICO 8
1.1 GENERALIDADES DE LA MODELACIÓN 8
1.2 DESCRIPCIÓN DEL MODELO 8
1.3 ESQUEMATIZACIÓN DEL MODELO 10
1.4 BALANCE DE FLUJO 11
1.5 BALANCE TÉRMICO 12
1.6 BALANCE DE MASA 12
1.7 REAIREACIÓN 13
1.8 OTROS ASPECTOS DEL MODELO 15
1.9 ENTRADAS Y SALIDAS DEL MODELO 15
1.10 MONITOREO DE LA CALIDAD DEL AGUA15
1.11 FUENTES PUNTUALES Y DIFUSAS 17
1.12 CONSTANTES HIDRÁULICAS 17
1.13 CONSTANTES CINÉTICAS 18
1.14 INFORMACIÓN METEOROLÓGICA 19
1.15 DESCRIPCIÓN DE LA ZONA DE ESTUDIO 20
1.16 BALANCE DE FLUJO 22
1.16.1 Cuerpo receptor: Río Iza 22
1.16.2 Aguas arriba del primer vertimiento municipal 36
1.16.3 Vertimiento colegio 37
1.16.4 Vertimiento principal municipio de Iza 37
1.17 MODELACIÓN PARA EL SISTEMA ACTUAL CON CAUDAL MEDIDO 39
1.17.1 Zona de mezcla 50
1.18 ESCENARIOS DE SIMULACIÓN 54
1.18.1 ESCENARIO I: Caudal medido sin vertimientos: 54
1.18.2 ESCENARIO II: Caudal medido con sistema de tratamiento de aguas residuales
62
1.18.3 ESCENARIO III: Caudal mínimo con los tres vertimientos 71
1.18.4 ESCENARIO IV: Caudal mínimo sin vertimiento 78
1.18.5 ESCENARIO V: Modelación con caudal mínimo y sistema de tratamiento de aguas
residuales 85
1.19 ANÁLISIS FINAL DEL EJERCICIO DE MODELACIÓN EN EL SISTEMA DEFINIDO
PARA EL RIO IZA 93
1.20 EVALUACIÓN OBJETIVOS DE CALIDAD DE LA FUENTE DENOMINADA RIO
TOTA, RECEPTOR DE LOS VERTIIENTOS MUNICIPALES, UBICADOS EN LA
JURISDICCIÓN DEL MUNICIPIO DE IZA. 96
1.21 REFERENCIAS BIBLIOGRAFICAS 99
 TABLA DE FIGURAS
Pág.

Figura 1. Esquema general QUAL2Kw 11

Figura 2. El esquema del balance de flujo se puede ver de la siguiente manera: 11
Figura 3. Esquema de balance térmico 12
Figura 4. Esquema de balance térmico 13
Figura [Link]ón de la zona de estudio Rio Iza 20
Figura 6. Localización vertimientos al Rio Iza 21
Figura 7. Balance de flujo de la zona Rio Iza 22
Figura 8. Método de aforo al Rio Iza en los dos puntos de estudio 23
Figura [Link]ón horaria del caudal del vertimiento Colegio 29
Figura [Link]ón horaria del caudal del vertimiento Principal 31
Figura 11. Vertimientos principal municipio de Iza 31
Figura [Link]ón horaria del caudal del vertimiento Centro de Acopio de Leche 34
Figura 13. Vertimientos Centro de Acopio de Leche 34
Figura 14. Diagrama lineal zona de estudio Rio Iza 36
Figura 15. Zona de mezcla en cuerpos naturales de agua 50
Figura 16. Localización vertimientos al Rio Iza 62
Figura 17. Localización vertimientos al Rio Iza 85
 LISTA DE TABLAS
Pág.
.
Tabla 1. Ecuaciones matemáticas para la determinación de la Ka.............................................14
Tabla 2. Constantes cinéticas para calibración..............................................................................19
Tabla 3. Caudal Rio Iza aguas arriba del vertimiento Colegio.....................................................24
Tabla 4. Calidad del Rio Iza aguas arriba del vertimiento colegio...............................................25
Tabla 5. Caudal Rio Iza aguas abajo del vertimiento del centro de acopio de leche................25
Tabla 6. Calidad del Rio Iza aguas abajo del vertimiento centro de acopio...............................27
Tabla 7. Caudales horarios del vertimiento (L/s)...........................................................................28
Tabla 8. Caudales horarios del vertimiento (L/s)...........................................................................30
Tabla 9. Resultados calidad del vertimiento principal...................................................................32
Tabla 10. Caudales horarios del vertimiento (L/s).........................................................................33
Tabla 11. Calidad del vertimiento del Centro de Acopio de Leche..............................................35
Tabla 12. Localización puntos a evaluar.........................................................................................38
Tabla 13. Resultados caracterización de puntos de evaluación..................................................39
Tabla 14. Presentación Modelo QUAL2KW6 – Río Iza.................................................................40
Tabla 15. Reporte WQ OUTPUT para corrida preliminar con caudal medido...........................41
Tabla 16. Resultados de la DBO.............................................................................................................42
Tabla 17. Resultados para el pH, SST, DQO, Conductividad, NO3 y Coliformes Totales........................44
Tabla [Link] de la modelación para el OD.......................................................................49
Tabla 19. Reporte WQ OUTPUT para escenario I Rio Iza...........................................................55
Tabla 20. Resultados de la modelación para la DBO.............................................................................56
Tabla 21. Resultados para el pH, SST, DQO y Coliformes Totales..........................................................57
Tabla [Link] de la modelación para el OD.............................................................................61
Tabla 23. Reporte WQ OUTPUT para escenario I Rio Iza...........................................................64
Tabla 24. Resultados de la modelación para la DBO.............................................................................65
Tabla 25. Resultados para el pH, SST, DQO y Coliformes Totales..........................................................66
Tabla [Link] de la modelación para el OD.............................................................................70
Tabla 27. Reporte WQ OUTPUT para escenario I Rio Iza...........................................................71
Tabla 28. Resultados de la DBO.............................................................................................................72
Tabla 29. Resultados para el pH, SST, DQO y Coliformes Totales..........................................................73
Tabla [Link] de OD.......................................................................................................................77
Tabla 31. Reporte WQ OUTPUT para escenario IV- Rio Iza.......................................................78
Tabla 32. Resultados de la DBO.............................................................................................................80
Tabla 33. Resultados para el pH, SST, DQO y Coliformes Totales..........................................................81
Tabla [Link] de resultados para el OD...................................................................................84
Tabla 35. Reporte WQ OUTPUT para escenario V - Rio Iza.......................................................86
Tabla 36. Resultados de la DBO.............................................................................................................87
Tabla 37. Resultados para el pH, SST, DQO y Coliformes Totales..........................................................88
Tabla [Link] de OD.......................................................................................................................92
Tabla 39. Indicadores y parámetros de calidad para el sector del Rio Tota..............................96
Tabla 40. Criterios, normas o estándares de calidad del cuerpo receptor.................................97
 LISTA DE GRÁFICAS
Pág.

Gráfica 1. Flujo en Rio Iza – Área de Influencia............................................................................41

Gráfica 2. Resultados temperatura.......................................................................................................42
Gráfica 3. Resultados de la DBO...........................................................................................................43
Gráfica 4. Comportamiento de la DQO en la distancia.........................................................................44
Gráfica 5. Comportamiento del pH en la distancia...............................................................................45
Gráfica 6. Sólidos Suspendidos Inorgánicos..........................................................................................46
Gráfica 7. Conductividad en el sistema.................................................................................................47
Gráfica 8. Nitritos y Nitratos en el sistema.....................................................................................48
Gráfica 9. Coliformes Totales................................................................................................................48
Gráfica 10. Perfil de OD en la distancia.........................................................................................49
Gráfica 11. Flujo en Rio Iza – Área de Influencia.........................................................................55
Gráfica 12. Resultados temperatura.....................................................................................................56
Gráfica 13. Resultados de la DBO.........................................................................................................57
Gráfica 14. Resultados DQO.................................................................................................................58
Gráfica 15. Sólidos Suspendidos Totales...............................................................................................59
Gráfica 16. Conductividad en el sistema...............................................................................................59
Gráfica 17. Nitritos y Nitratos para el sistema evaludos........................................................................60
Gráfica 18. Coliformes Totales..............................................................................................................60
Gráfica 19. Perfil de OD en la distancia.........................................................................................61
Gráfica 20. Flujo en Rio Iza – Área de Influencia.........................................................................64
Gráfica 21. Resultados temperatura.....................................................................................................65
Gráfica 22. Resultados de la DBO.........................................................................................................66
Gráfica 23. Perfil de DBO en la distancia..............................................................................................67
Gráfica 24. Perfil de DQO en la distancia..............................................................................................67
Gráfica 25. Sólidos Suspendidos Totales...............................................................................................68
Gráfica 26. Conductividad en el sistema evaluado................................................................................68
Gráfica 27. Nitritos y Nitratos para el sistema evaluado.......................................................................69
Gráfica 28. Coliformes Totales..............................................................................................................69
Gráfica 29. Perfil de OD en la distancia.........................................................................................70
Gráfica 30. Flujo en Rio Iza – Área de Influencia.........................................................................71
Gráfica 31. Resultados temperatura.....................................................................................................72
Gráfica 32. Resultados de la DBO.........................................................................................................73
Gráfica 33. Perfil de DQO en la distancia..............................................................................................74
Gráfica 34. Comportamiento del pH en el sistema evaluado...............................................................74
Gráfica 35. Sólidos Suspendidos Inorgánicos........................................................................................75
Gráfica 36. Conductividad en el sistema...............................................................................................76
Gráfica 37. Nitritos y Nitratos para el sistema evaluado.......................................................................76
Gráfica 38. Coliformes Totales..............................................................................................................77
Gráfica 39. Perfil de OD en la distancia.........................................................................................78
Gráfica 40. Flujo en Rio Iza – Área de Influencia...........................................................................79
Gráfica 41. Resultados temperatura.....................................................................................................79
Gráfica 42. Resultados de la DBO.........................................................................................................80
Gráfica 43. Resultados DQO.................................................................................................................81
Gráfica 44. Comportamiento del pH en la distancia..............................................................................82
Gráfica 45. Sólidos Suspendidos Totales..............................................................................................82
Gráfica 46. Conductividad en el sistema...............................................................................................83
Gráfica 47. Nitritos y Nitratos en el sistema.........................................................................................83
Gráfica 48. Coliformes Totales..............................................................................................................84
Gráfica 49. Perfil de OD en la distancia..........................................................................................85
Gráfica 50. Flujo en Rio Iza – Área de Influencia.........................................................................86
Gráfica 51. Resultados temperatura.....................................................................................................87
Gráfica 52. Resultados de la DBO.........................................................................................................88
Gráfica 53. Perfil de DQO en la distancia..............................................................................................89
Gráfica 54. pH en la distancia...............................................................................................................89
Gráfica 55. Sólidos Suspendidos Inorgánicos........................................................................................90
Gráfica 56. Conductividad reportada en la distancia............................................................................90
Gráfica 57. Nitritos y Nitritos en la distancia........................................................................................91
Gráfica 58. Coliformes Totales..............................................................................................................92
Gráfica 59. Perfil de OD en la distancia.........................................................................................93
1. PREDICCIÓN A TRAVÉS DE MODELOS DE SIMULACIÓN DE LOS IMPACTOS
QUE CAUSE EL VERTIMIENTO AL CUERPO DE AGUA, EN FUNCIÓN DE LA
CAPACIDAD DE ASIMILACIÓN Y DILUCIÓN DEL CUERPO DE AGUA RECEPTOR
Y DE LOS USOS Y CRITERIOS DE CALIDAD ESTABLECIDOS EN EL PLAN DE
ORDENAMIENTO DEL RECURSO HÍDRICO

1.1 GENERALIDADES DE LA MODELACIÓN

Un modelo es una representación de un sistema, lo que lleva a obtener respuestas a

través de experimentos. No siempre se requiere el uso de programas de computador,
sin embargo, hoy en día, los modelos se han llevado a codificaciones en lenguajes
de programación. Modelar es de gran relevancia para el conocimiento de un sistema
y en este caso de los sistemas hídricos. Normalmente los fenómenos naturales se
representan por ecuaciones matemáticas, cuya mejor manera de resolverlas, es
elaborando un programa de computador.

Por otro lado, se cuenta con los procesos de simulación, los cuales son experimentos
que se hacen con un modelo. Las predicciones de calidad de agua son utilizadas
para conocer la variación de los contaminantes en el tiempo, es decir cuando se
necesita realizar simulaciones dinámicas.

Las normas han definido dos aspectos claves en los procesos de modelación, los
cuales están asociados con los niveles de asimilación de sustancias biodegradables
o acumulativas y la capacidad de dilución de sustancias no biodegradables en el
contexto de la calidad del agua.

Por lo anterior se define el uso del modelo de calidad del agua QUAL2Kw 1, software
de actualidad y uso en el contexto internacional para modelar la calidad del agua en
ríos y corrientes de agua superficiales.

1.2 DESCRIPCIÓN DEL MODELO

El modelo de calidad del agua QUAL2K fue desarrollado para la Agencia de

Protección Ambiental de Estados Unidos por Chapra y Pelletier en el año 2003,
desde ese momento ha venido mejorándose, recibiendo el nombre de QUAL2Kw en
el año 2008.

La nueva versión del modelo tiene la capacidad de simular una corriente principal y
tres corrientes secundarias, las cuales pueden ser manejadas de manera

1
Una de las razones por las cuales QUAL2Kw es uno de los programas de más uso en nuestro medio, es porque se encuentra
desarrollado en formato Excel® lo cual le da gran flexibilidad en el manejo de la información, amigable con los usuarios y
facilidad para entregar la información de resultados de forma gráfica y por medio de tablas.
independiente o integrarse a la corriente principal dependiendo de las necesidades
del usuario.

De manera general se pueden establecer las características más representativas del

modelo:

 El modelo define tramos y longitud de los mismos para modelar, por lo cual se
pueden incluir múltiples fuentes puntuales o dispersas en cualquier tramo.

 Modelo unidimensional en la dirección del flujo 2. El sistema hídrico, se define bien

mezclado tanto vertical como horizontalmente. Condiciones hidráulicas estables 3.

 El modelo permite realizar simulaciones en casos donde la simulación se extiende

no solo a la corriente principal, sino también a tributarios.

 El calor del volumen de agua y la temperatura son simulados como función de

aspectos meteorológicos en una escala de variación diurna.

 Todas las variables de la calidad del agua son simuladas en escala de tiempo
diurna, aunque el modelo funciona en flujo permanente, tiene la opción de trabajo en
la modalidad de un modelo cuasi-dinámico.

 Se simulan vertimientos y extracciones tanto puntuales como dispersas.

 El modelo QUAL2Kw usa dos formas de DBO carbonácea para representar el

carbono orgánico. Estas formas son: la DBO de rápida oxidación (DBO fast) y la
DBO de lenta oxidación (DBO slow)4.

 En la práctica se puede asumir lo siguiente: La DBO rápida corresponde a la

descomposición de la materia orgánica disuelta fácil de degradar. Por esta razón se
puede asumir que la DBOfast es igual a la DBO filtrada. Para el caso de la DBOslow
o lenta, corresponde a aquella materia orgánica que se demora mucho más tiempo
en degradarse

 Se modela un patógeno genérico (coliformes fecales o coliformes totales). La

remoción de patógenos está determinada como función de la temperatura, la luz y la
sedimentación.

 Los flujos en la interfase agua-sedimento de OD y nutrientes se simulan

internamente en el modelo. De esta forma el flujo de demanda béntica (SOD) y de
2
En la mayoría de los ríos las dimensiones de la longitud son mucho mayores a las profundidades y
anchos.
3
Flujo permanente, periodos de caudal y cargas constantes en el tiempo.
4
Adicionalmente, se simula la materia orgánica particulada (Detritus), este material de detritus está
compuesto por partículas de carbón, nitrógeno y fósforo en una estequiometria fija.
nutrientes son simulados como función de la sedimentación de materia orgánica
particulada, reacciones dentro de los sedimentos, y de concentraciones de las
formas solubles en la columna de agua.

 El modelo establece limitaciones relacionadas con la flexibilidad, por lo cual no

permite simular algunos parámetros, como los vertidos de derivados del petróleo. Por
otro lado, el modelo es unidimensional y por tanto asume que el río tiene una mezcla
perfecta lateral y vertical5.

1.3 ESQUEMATIZACIÓN DEL MODELO

La esquematización del modelo matemático de simulación QUAL2Kw corresponde a

la representación de las características del sistema que se requiere modelar. En el
caso de la modelación en corrientes superficiales dicha representación se refiere a
las características hidrogeométricas y la definición de las fronteras del modelo
(fronteras internas y externas6).

En el momento de seleccionar los tramos entre las estaciones se debe tener en

cuenta que las propiedades físicas e hidráulicas 7, químicas y biológicas permanecen
constantes a lo largo de cada tramo.

El tramo de la corriente que se estudia se divide en una sucesión de pequeños

elementos computacionales que operan como reactores completamente mezclados.
Para cada elemento computacional se realiza un balance hidrológico en términos del
caudal, un balance térmico en términos de la temperatura y un balance de masa en
términos de la concentración de cada compuesto. Se tienen en cuenta procesos de
transporte8, fuentes externas9 y fuentes internas10.

En la Figura siguiente se puede observar el esquema general del modelo QUAL2Kw

con sus fronteras externas e internas.

Figura 1. Esquema general QUAL2Kw

5
Esta es una asunción razonable para la mayoría de los ríos, excepto cuando los ríos son
particularmente anchos, profundos o tienen un movimiento lento. En estos casos debería considerarse
un modelo bidimensional.
6
Las fronteras externas corresponden a las estaciones de monitoreo ubicadas aguas arriba y aguas
abajo del tramo en estudio, mientras que las fronteras internas corresponden a los ríos tributarios, las
extracciones y los vertimientos para los cuales se dispone de información de caudales y calidad del
agua.
7
Sección transversal, pendiente de la corriente y rugosidad.
8
Advección y dispersión,
9
Descarga de agua residual o extracciones de agua.
10
Demanda béntica o transformaciones bioquímicas.
 Fuente: Chapra [Link]., 2008.

Las ecuaciones de balance resultante se resuelven para el caso de flujo permanente

con un método clásico de diferencias finitas implícito. Los resultados finales se
traducen en curvas que muestran la variación de los parámetros modelados a lo
largo de la corriente.

1.4 BALANCE DE FLUJO

El modelo trabaja bajo régimen estacionario, o sea que el flujo entrante y saliente en
cada uno de los elementos computacionales en análisis se mantiene constante a
través del tiempo.

Figura 2. El esquema del balance de flujo se puede ver de la siguiente manera:

Fuente: Chapra [Link]., 2008.

En donde:
Qi: Caudal de salida del elemento computacional
Qi-1: Caudal de entrada del elemento computacional
Qin,i: Caudal de entrada neto lateral por fuentes puntuales y difusas
Qout,i: Caudal de salida neto lateral por fuentes puntuales y difusas.

1.5 BALANCE TÉRMICO

El balance de temperatura en cada uno de los elementos que conforman el sistema

es función de la transferencia de calor de los elementos adyacentes, de la
transferencia de calor entre los sedimentos, el agua y la atmósfera y de las entradas
y salidas externas al sistema.

Figura 3. Esquema de balance térmico

Fuente: Chapra [Link]., 2008.

1.6 BALANCE DE MASA

El modelo considera por cada constituyente 11 un balance general de masa en cada

elemento computacional, el cual involucra los efectos de dispersión, dilución,
consumo o generación de constituyentes por reacciones químicas o bioquímicas y la
generación o pérdida de nutrientes por fuentes externas o internas 12
11
Parámetro fisicoquímico o biológico.
12
Descargas puntuales y captaciones.
Dónde:

Ci: Concentración
Qi: Flujo en el enésimo segmento
Qab: Salida en el elemento i debido a extracción puntual o difusa
Vi: Volumen incremental del enésimo segmento
t: Tiempo
Si: Fuentes y sumideros del constituyente debido a reacciones y mecanismos de
transferencia de masa
Ei: Coeficiente de dispersión longitudinal entre segmentos i y i+1
Si: Carga externa del constituyente en el elemento i.

Figura 4. Esquema de balance térmico

Fuente: Chapra [Link]., 2008

1.7 REAIREACIÓN

En este contexto se busca establecer la velocidad a la cual el oxígeno es transferido

a la masa de agua, para lo cual a través de la tasa de reaireación o también conocida
como constante de reaireación (Ka), se puede establecer esta situación.

QUAL2Kw, utiliza varias ecuaciones matemáticas para hacer la predicción de esta

constante, entre las cuales están: la ecuación de O´Connor y Dobbins, la ecuación
de Churchill, la ecuación de Owens, la ecuación de Tsivoglou, la ecuación de
Thanckston y Dawsan y las ecuaciones presentadas por USGS.
En la siguiente tabla se pueden observar las diferentes ecuaciones matemáticas que
permiten llegar a la determinación de la constante de reaireación Ka.

Tabla 1. Ecuaciones matemáticas para la determinación de la Ka

Fuente. Sierra Ramírez, Carlos. Calidad del Agua. Evaluación y Diagnóstico.

Universidad de Medellín.2011

1.8 OTROS ASPECTOS DEL MODELO

A continuación, se presentan otros aspectos relacionados con la implementación del
modelo de calidad del agua QUAL2Kw, entre los cuales están las entradas y salidas,
la información de calidad del agua, las fuentes definidas como descargas, los
aspectos relacionados a la hidráulica del sistema y finalmente lo relacionado a las
constantes cinéticas.

1.9 ENTRADAS Y SALIDAS DEL MODELO

 El libro de Excel contiene varias hojas de trabajo en las cuales se ingresan los
diferentes datos de entrada del modelo y desde donde éste es ejecutado.

 Cada hoja de trabajo tiene un nombre que identifica el tipo de información que
contiene y que puede resumirse de la siguiente forma: parámetros fisicoquímicos y
microbiológicos tomados en campo en las campañas de monitoreo, descripción de
los tramos del río, constantes hidráulicas, los valores de las constantes cinéticas de
calibración del modelo, condiciones meteorológicas (temperatura del aire,
temperatura del punto de rocío, velocidad del viento, nubosidad y sombra), y la
información fisicoquímica correspondiente a fuentes puntuales y difusas como
tributarios, captaciones y vertimientos

 Después de que el modelo QUAL2Kw ha efectuado los cálculos necesarios para

poder realizar las corridas de simulación, se generan diferentes hojas electrónicas y
gráficas con los resultados (datos de salida).

 En las gráficas resultantes se pueden comparar las curvas correspondientes al

perfil de concentraciones de los diferentes parámetros estimados por el modelo a lo
largo del río, con los datos medidos en campo en las estaciones de monitoreo

1.10 MONITOREO DE LA CALIDAD DEL AGUA

 La ejecución del modelo implica contar con datos históricos o recientes de

campañas de monitoreo de calidad del agua en estaciones ubicadas sobre la
corriente que se requiere modelar, de tal forma que sea posible realizar la calibración
de las constantes cinéticas que mejor representan las características de la corriente.

 Para el modelo, las fronteras externas corresponden a las estaciones de

monitoreo ubicadas aguas arriba y aguas abajo del tramo de estudio.

 El modelo permite ingresar los valores de calidad del agua medidos en la corriente
de forma puntual, el valor medio de los datos medidos en un tiempo determinado y
los valores mínimos y máximos observados para cada parámetro de calidad de agua
en las estaciones monitoreadas.
 El QUAL2Kw permite modelar los siguientes parámetros fisicoquímicos y
microbiológicos: Conductividad, sólidos suspendidos inorgánicos, oxígeno disuelto,
DBO rápida, DBO lenta, nitrógeno orgánico disuelto, nitrógeno amoniacal, nitratos,
fósforo orgánico disuelto, fósforo inorgánico, fitoplancton, detritus, patógenos,
alcalinidad, carbono orgánico total, algas de fondo, temperatura y caudal.

 La demanda bioquímica de oxígeno carbonácea (DBOC), estas dos formas 13 son

la DBO rápida (BOD fast) que es la materia orgánica fácilmente oxidada por los
microorganismos y la DBO lenta (BOD slow) que es la materia orgánica que es
difícilmente biodegradable por la población bacteriana.

 La DBO última o total es la medida hasta el agotamiento total de la materia

orgánica total, lo que usualmente ocurre en un periodo de 20 días 14.

 La cantidad de oxigeno remanente en un tiempo determinado es la Demanda

Bioquímica de Oxigeno a ese tiempo, esta medida de la DBO es proporcional a la
cantidad de material orgánico biodegradable o sustrato.

 La velocidad de consumo de sustrato sigue una cinética de primer orden, lo cual

matemáticamente se puede establecer de la siguiente manera 15:

Dónde:

yt: demanda bioquímica de oxígeno en un tiempo t

L0: DBO última
K16: constante cinética de reacción
t : tiempo.

1.11 FUENTES PUNTUALES Y DIFUSAS

13
Al modelo es necesario ingresar los valores de DBO rápida y DBO lenta en términos de DBO última,
esto supone una dificultad teniendo en cuenta que actualmente la materia orgánica es medida en el
laboratorio a 5 días (DBO5) por cuestiones de tiempo y presupuesto
14
esta oxidación se alcanza en dos etapas: en la primera se estabiliza la materia carbonácea y
corresponde a la primera fase de la desoxigenación; mientras que en la segunda fase se oxidan los
compuestos nitrogenados carbonáceos (Metcalf [Link]., 1998).
15
Pelletier [Link]., 2004):
16
Experimentalmente se ha encontrado que la constante cinética de reacción para la DBO en aguas
superficiales y residuales varía entre 0.05 y 0.3 d-1, estableciéndose que el valor más común es de
K≈0.23 d-1. Para aguas de otro tipo, como por ejemplo las aguas producidas por industrias K varía de
0.12 a 0.7 d-1 (Pelletier [Link]., 2004).
Las fuentes de contaminación del agua se dividen en fuentes puntuales 17 y fuentes
difusas18. En el QUAL2Kw las fronteras internas del modelo corresponden a las
fuentes puntuales (ríos, extracciones y vertimientos) y a las fuentes difusas
(vertimientos no puntuales) para los cuales se dispone de información de caudales y
calidad del agua ya que han sido previamente identificados en campo.

1.12 CONSTANTES HIDRÁULICAS

Este tipo de modelos implican un trabajo adecuado y pertinente en los relacionado a

la hidráulica, ya que la velocidad media de la corriente, la pendiente del tramo, el
ancho y la profundidad de la sección transversal, determinan la tasa de reaireación
de la corriente y por lo tanto su capacidad de autodepuración, así mismo la velocidad
del flujo determinan los tiempos de viaje o de residencia de los contaminantes en el
tramo. (Chapra [Link]., 2008).

El modelo define dos opciones para las constantes hidráulicas. La primera opción
está relacionada con Rating Curves y el segundo utiliza la fórmula de Manning. La
elección de cualquiera de estos métodos depende de la información con la que se
cuente y la confiabilidad de la misma.

Rating Curves

Donde:

Q: caudal
U: velocidad
H: profundidad
W: ancho
k, c, a: constantes
f, b, m: coeficientes

Ecuación de Manning:

17
la primera consiste en fuentes que entran al medio ambiente por un solo punto de entrada que se
identifica fácilmente Las fuentes puntuales tienden a ser aquellas directamente identificadas por sus
orígenes en actividades humanas y son descargadas por tuberías y alcantarillas a los cuerpos de
agua (fábricas, plantas de tratamiento de aguas residuales, minas, pozos petroleros, etc).
18
mientras que las fuentes difusas son aquellas que se reparten en áreas más extensas. Las fuentes
difusas por el contrario son relativamente más difíciles de identificar y controlar (Vertimiento de
sustancias químicas, tierras de cultivo, lotes para pastar ganado, construcciones, tanques sépticos)
(Manahan, 2006).
 Donde:

Q: caudal
So: pendiente del fondo
Ac: área transversal
P: perímetro mojado
n: coeficiente de Manning.

1.13 CONSTANTES CINÉTICAS

El QUAL2Kw permite definir una gran cantidad de variables cinéticas para describir
todas las reacciones que tienen lugar en la dispersión y la dilución de los
contaminantes en una corriente de agua, sin embargo, aunque todas las variables se
tienen en cuenta en la aplicación del modelo, es elección del usuario cuáles de estas
variables desea calibrar.

En la siguiente tabla se presentan las constantes cinéticas que permite calibrar el

QUAL2Kw, la descripción de cada una de ellas y el rango de valores de calibración
encontrados en la bibliografía a través de múltiples ensayos para cada constante 19.

Los valores iniciales de las constantes cinéticas a calibrar pueden ser definidos por el
usuario de acuerdo a experiencias en la región o lo puede asignar arbitrariamente
siempre y cuando estén entre los rangos establecidos.

Tabla 2. Constantes cinéticas para calibración

19
Chapra, 1997; Kannel [Link]., 2007; Pelletier [Link]., 2003.
 Fuente. Chapra, 1997; Kannel [Link]. 2007; Pelletier [Link]. 2003.

1.14 INFORMACIÓN METEOROLÓGICA

La información meteorológica necesaria para la implementación del modelo consiste

en datos de temperatura, temperatura de punto de rocío, velocidad del viento,
nubosidad, sombra en el río producida por la topografía o vegetación y radiación
solar.

Estos datos y la temperatura del agua del río permiten realizar el balance de energía
térmica para el cálculo de la temperatura del agua. Esta información puede ser
tomada en campo, se puede consultar la base de datos de la National Oceanic &
Atmospheric Administration (NOAA), o se pueden consultar estaciones
meteorológicas existentes en la zona para obtener los valores correspondientes al día
en el que se realizó el monitoreo de calidad del agua

1.15 DESCRIPCIÓN DE LA ZONA DE ESTUDIO

La red hidrográfica del municipio de Iza, forma parte de la cuenca alta del Río
Chicamocha, vertiente sur y suroriental. En dichas vertientes las corrientes que
fluyen al río, atraviesan y conforman variados paisajes con características
morfológicas bien definidas: La planicie aluvial y lacustre formada por los ríos Iza,
Pesca y Chiquito; las vertientes o montañas al oriente del municipio y las colinas al
noroccidente, con drenajes de orden inferior como las Quebradas Carichana,
Chiguatá, Toquilla, La Fragua y Agua Caliente, entre otras 20.

Figura [Link]ón de la zona de estudio Rio Iza

Fuente. INCO A&J S.A.S – Google Earth

Rio Tota: es la cuenca más importante del municipio; tiene su nacimiento en el
Páramo Las Alfombras a 3.600 m.s.n.m. En la parte alta, su curso recoge gran parte
del caudal producto de la abundante precipitación en la región; a medida que
desciende, el clima local cambia y las vertientes ahora secas, muestran marcados

20
EOT, 2000.
fenómenos de erosión. En la parte baja, conforma un pequeño valle, donde se
localiza gran parte del área urbana municipal.

La corriente del río Tota más conocido en el área puntual, como Río Iza, cruza el
municipio en dirección SE-NW, en una longitud de 5Km aproximadamente. El rio
Tota o Iza tiene una longitud total de 29 Km presenta como afluentes principales por
su margen izquierda las Quebradas Agua Caliente, Seca y La Chorrera, por su
margen derecho las Quebradas Toquilla, Chiguatá, Cruz Colorada y La Fragua,
además de otras corrientes menores y luego de pasar por el área urbana,
desemboca en el rio Chicamocha.

El objeto de la modelación de la calidad del agua para el caso, corresponde a un

sector del Rio Iza, que cumple la función de recibir el vertimiento principal del
municipio y los vertimientos secundarios del colegio y del centro de acopio de leche.

Figura 6. Localización vertimientos al Rio Iza

Fuente. INCO A&J S.A.S – Google Earth

La figura representa el sistema a modelar definido para sus límites con un punto
aguas arriba y otra agua abajo del sector. El inicio se define como aguas arriba del
primer vertimiento en Puente Colegio y su final en un punto aguas abajo del
vertimiento del centro de acopio de leche. También se pueden observar los
 vertimientos secundarios que corresponderían al Colegio y al Centro de acopio y
finalmente estaría el vertimiento principal del municipio de Iza.

1.16 BALANCE DE FLUJO

A continuación, se presenta el balance de flujo para el Rio Iza en el sector evaluado.

La figura, representa el diagrama conceptual de los flujos en la zona de estudio,
ubicando el puto aguas arriba a 284 metros del vertimiento del colegio y para el
punto aguas abajo se midieron 103 metros a partir del vertimiento del centro de
acopio de leche. Del vertimiento del colegio al vertimiento principal del municipio, se
tienen 1.378 m y de este al vertimiento de Centro de Acopio se tienen 0.725 Km.

Figura 7. Balance de flujo de la zona Rio Iza

Fuente: INCO A&J S.A.S.

1.16.1 Cuerpo receptor: Río Iza

Se llevó a cabo un aforo por el método área – velocidad, para cada uno de los puntos
de interés dentro del sistema evaluado. Bajo este método, el caudal corresponde al
volumen de agua que pasa instantáneamente por la sección de aforos en la unidad
de tiempo y se expresa en m3/s o en L/s; es un método utilizado a escalas menores.

El caudal en una corriente de agua es función del área de la sección de aforos y de

la velocidad media del flujo y se obtiene mediante el producto de estas dos variables.
 Figura 8. Método de aforo al Rio Iza en los dos puntos de estudio

Fuente: Hydrology and hydraulic systems (Gupta 1989).

En los puntos de interés para el estudio, también se hizo el muestreo para la

determinación de la calidad del recurso. Los parámetros analizados, son los
contemplados en el Artículo 8 de la Resolución 631 de 2015, en lo relativo a aguas
residuales domésticas de los prestadores del servicio público de alcantarillado a
cuerpos de agua superficial, con una carga menor a 625 kg/día de DBO 5.

La información del Rio Iza para los puntos aguas arriba y aguas abajo, se presenta
inicialmente para establecer los límites del sistema; posteriormente se presenta la
información de los vertimientos, sin embargo, se aclara que para el vertimiento del
colegio no se realizó la caracterización de la muestra, tomando como base las
condiciones del vertimiento principal del municipio, además que el aforo en este
punto estableció menos de un litro por segundo como caudal promedio del aporte.

A continuación, se presenta la información para cada uno de los puntos de interés en

lo concerniente a aforos y calidad.

Rio Iza aguas arriba del vertimiento Colegio

Este punto está definido a través de una longitud de 0.284 Km aguas arriba del
vertimiento colegio, sus coordenas geográficas corresponden a longitud 75°58´51.73”
W y latitud de 5°36´21.39” N. La altura sobre el nivel del mar es de 2529. En lo
relacionado al flujo, se pudo establecer un caudal de 416.21 L/s o 0.42 m 3/s.

Tabla 3. Caudal Rio Iza aguas arriba del vertimiento Colegio

Lugar de Muestreo Aguas arriba de vertimiento Colegio Fecha 25/09/2018
Este 1121486 Norte 1111774
Elevación 2529 Ancho del flujo (m) 5.4
Ancho parcial por sección 0.54 No. Secciones 10
Área trasversal Velocidad media en la
Sección Caudal (m3/s) Caudal (L/s)
(m2) sección (m/s)
1 0.057 0.000 0.0000 0.00
2 0.147 0.188 0.0276 27.61
3 0.183 0.343 0.0627 62.73
4 0.195 0.349 0.0679 67.88
5 0.170 0.373 0.0634 63.38
6 0.177 0.321 0.0567 56.71
7 0.174 0.284 0.0494 49.39
8 0.171 0.291 0.0499 49.87
9 0.185 0.202 0.0374 37.41
10 0.032 0.038 0.0012 1.22
Caudal total 0.42 416.21

Perfil hidraulico del río Iza aguas abajo de los vertimientos (m)

0
0 0.54 1.08 1.62 2.16 2.7 3.24 3.78 4.32 4.86 5.4
Profundidad (m)

-0.1

-0.2

-0.3

-0.4

Sección transversal aguas arriba Rio Iza

Fuente: INCO A&J S.A.S.

Tabla 4. Calidad del Rio Iza aguas arriba del vertimiento colegio
 Fuente: Analizar Laboratorio Ltda

Rio Iza aguas abajo del vertimiento del Centro de acopio de leche

Este punto está definido a través de una longitud de 0.103 Km aguas abajo del
vertimiento del cetro de acopio de leche, sus coordenas geográficas corresponden a
longitud 72°59´8.96” W y latitud de 5°37´28.05” N. La altura sobre el nivel del mar es
de 2530. El aforo estableció un caudal de 207.17 L/s o 0.21 m 3/s.

Tabla 5. Caudal Rio Iza aguas abajo del vertimiento del centro de acopio de
leche
Lugar de Muestreo Aguas abajo del vertimiento Fecha 25/09/2018
Este  1120935 Norte 1113915
Elevación 2515 Ancho del flujo (m) 4.79
Ancho parcial por sección 0.479 No. Secciones 10
Sección Área trasversal Velocidad media en Caudal (m3/s) Caudal (L/s)
 (m2) la sección (m/s)
1 0.023 0.015 0.0003 0.35
2 0.071 0.020 0.0014 1.42
3 0.108 0.405 0.0437 43.70
4 0.074 0.189 0.0140 14.02
5 0.081 0.233 0.0189 18.94
6 0.091 0.369 0.0336 33.60
7 0.108 0.369 0.0397 39.72
8 0.125 0.280 0.0349 34.88
9 0.105 0.183 0.0193 19.29
10 0.062 0.020 0.0013 1.25
Caudal total 0.21 207.17

Sección transversal aguas abajo Rio Iza

Fuente: INCO A&J S.A.S.

Tabla 6. Calidad del Rio Iza aguas abajo del vertimiento centro de acopio
 Fuente: Analizar Laboratorio Ltda.

Vertimiento Colegio

Este punto está definido a través de una longitud de 0.284Km del punto de inicio del
sistema, sus coordenas geográficas corresponden a longitud 75°58'51.50" W y latitud
de 5°36'27.12.51" N. La altura sobre el nivel del mar es de 2527. El aforo volumétrico
estableció un caudal de 0.966 L/s.

Coordenadas A.S.N.M
PUNTO Planas Coordenadas Geográficas .
X Y LONGITUD LATITUD
Vertimiento 75°58'51.50" 5°36'27.12.51" 2527
Colegio 1121493 1111950 E N
Tabla 7. Caudales horarios del vertimiento (L/s)
Caudal instantáneo
No. Toma Hora
(L/s)
1 9:00 0.847
2 10:00 2.978
3 11:00 2.444
4 12:00 1.153
5 13:00 2.556
6 14:00 1.889
7 15:00 0.931
8 16:00 0.500
9 17:00 0.506
10 18:00 2.056
11 19:00 1.153
12 20:00 0.431
13 21:00 0.417
14 22:00 0.700
15 23:00 0.000
16 0:00 0.000
17 1:00 0.000
18 2:00 0.000
19 3:00 0.000
20 4:00 0.000
21 5:00 0.000
22 6:00 2.667
23 7:00 2.083
24 8:00 0.600

Caudal Promedio (L/s) 0.966

Fuente: INCO A&J S.A.S.

 Figura [Link]ón horaria del caudal del vertimiento Colegio

Fuente: INCO A&J S.A.S.

El aforo realizado en un lapso de 24 horas, definió valores máximos de 2.978 L/s a

las 10:00am y caudales nulos en el lapso comprendidos entre las 11:00pm hasta las
5:00am. Este comportamiento es natural por las actividades propias que se
desarrollan en el colegio. El promedio de caudal se estableció en < 1.0 L/s.

Vertimiento principal del municipio de Iza

Este punto está definido a través de una longitud de 1.378 km desde el vertimiento
colegio. Sus coordenas geográficas corresponden a longitud 72°58'57.78" W y latitud
de 5°37'10.94" N. La altura sobre el nivel del mar es de 2531 y el aforo volumétrico
estableció un caudal de 2.216 L/s.

Coordenadas A.S.N.M
PUNTO Planas Coordenadas Geográficas .
X Y LONGITUD LATITUD
Vertimiento 72°58'57.78" 2531
Principal 1121297 1113296 W 5°37'10.94" N
Tabla 8. Caudales horarios del vertimiento (L/s)
RESUMEN DE RESULTADOS
Caudal instantáneo
No. Toma Hora
(L/s)
1 9:00 2.128
2 10:00 1.726
3 11:00 2.119
4 12:00 1.490
5 13:00 2.202
6 14:00 2.590
7 15:00 1.982
8 16:00 2.273
9 17:00 2.101
10 18:00 2.129
11 19:00 1.835
12 20:00 2.419
13 21:00 2.295
14 22:00 1.984
15 23:00 2.269
16 0:00 2.345
17 1:00 2.387
18 2:00 2.393
19 3:00 2.499
20 4:00 2.367
21 5:00 2.452
22 6:00 2.077
23 7:00 2.714
24 8:00 2.404

Caudal Promedio (L/s) 2.216

Fuente: INCO A&J S.A.S

 Figura [Link]ón horaria del caudal del vertimiento Principal

Fuente. INCO A&J S.A.S

El aforo realizado en un lapso de 24 horas, definió valores máximos de 2.714 L/s a

las 7:00am y mínimos a los 12:00 del mediodía, con un valor de 1.490 L/s. El
promedio de caudales en ese lapso de 24 horas fue de 2.216 L/s.

Figura 11. Vertimientos principal municipio de Iza

Localización punto de muestreo vertimiento principal

Fuente: INCO A&J S.A.S.

 Tabla 9. Resultados calidad del vertimiento principal

Fuente: Analizar Laboratorio Ltda.

Vertimiento del Centro de acopio de leche

Este punto está definido a través de una longitud de 0.103 km hasta el punto final del
sistema, sus coordenas geográficas corresponden a longitud 72°59'8.96" W y latitud
de 5°37'28.05" N. La altura sobre el nivel del mar es de 2530. El aforo volumétrico
estableció un caudal de 1.188 L/s.

Coordenadas A.S.N.M
PUNTO Planas Coordenadas Geográficas .
X Y LONGITUD LATITUD
Vertimiento 72°59'8.96" 2530
Principal 1120952 1113821 W 5°37'28.05" N
Tabla 10. Caudales horarios del vertimiento (L/s)
RESUMEN DE RESULTADOS
Caudal instantáneo
No. Toma Hora
(L/s)
1 10:00 1.140
2 11:00 0.949
3 12:00 1.108
4 13:00 1.212
5 14:00 1.206
6 15:00 1.040
7 16:00 1.124
8 17:00 1.121
9 18:00 1.133
10 19:00 1.057
11 20:00 1.022
12 21:00 1.022
13 22:00 0.924
14 23:00 0.916
15 0:00 0.918
16 1:00 1.161
17 2:00 1.014
18 3:00 1.409
19 4:00 1.437
20 5:00 1.380
21 6:00 1.238
22 7:00 1.584
23 8:00 1.763
24 9:00 1.633

Caudal promedio (Lps) 1.188

 Figura [Link]ón horaria del caudal del vertimiento Centro de Acopio de
Leche

Fuente: INCO A&J S.A.S.

El aforo realizado en un lapso de 24 horas, definió valores máximos de 1.763 L/s a

las 8:00am y mínimos a las 11:00pm, con un valor de 0.916 L/s. El promedio de
caudales en ese lapso de 24 horas fue de 1.188 L/s.

Figura 13. Vertimientos Centro de Acopio de Leche

Localización punto de muestreo vertimiento Lechería

Fuente: INCO A&J S.A.S.

Tabla 11. Calidad del vertimiento del Centro de Acopio de Leche

Fuente: Analizar Laboratorio Ltda.

 Figura 14. Diagrama lineal zona de estudio Rio Iza

Fuente: INCO A&J S.A.S.

1.16.2 Aguas arriba del primer vertimiento municipal

Para efectos del modelo la cabecera está ubicada en el kilómetro 2.49 del sistema
evaluado y a 284 metros del vertimiento del colegio. El caudal medido fue de 416.21
 L/s, dato reportado por medio del aforo por el método área – velocidad. En cuento a
la calidad de la fuente se reportan valores de pH de 7.37, DQO de 28 mg O 2 /L, para
la DBO5 total <2 mg O2 /L, para el OD de 5.40 mg O 2 /L y para coliformes totales de
680 NMP / 100ml.

1.16.3 Vertimiento colegio

En este punto de descarga y que corresponde al colegio, se realizó aforo volumétrico

durante 24 horas. Se pueden apreciar variaciones durante el transcurso del día, con
un máximo de 2.978 L/s a las 10:00am y un mínimo de 0.0 L/s en el periodo del día
comprendido entre las 11:00pm hasta las 5:00 am. Luego de hacer el promedio de
caudales durante el día, se obtuvo un valor de 0.996 L/s, el cual será utilizado para el
balance de flujo y la modelación.

Teniendo en cuenta que el promedio de caudales es muy bajo, aproximado a un litro

por segundo, no se realizó la caracterización del vertimiento, tomándose como
referencia la caracterización que se realizó para el vertimiento principal del municipio.

1.16.4 Vertimiento principal municipio de Iza

En el caso del vertimiento principal del municipio, se realizó aforo volumétrico

durante 24 horas, así como la respectiva caracterización, tomando alícuotas cada
hora. La variación de caudales no fue tan alta, reportando un valor promedio de
2.216 L/s, con máximos de 2.590 L/s a las 2:00pm y un mínimo de 1.490 L/s al medio
día.

En cuento a la calidad del vertimiento se reportan valores de pH de 7.14, DQO de

109 mg O2 /L, para la DBO5 total es de 55 mg O 2 /L, para el OD <0.08 mg O 2 /L y
para coliformes totales de 12 x 10 3 NMP / 100ml, entre otros parámetros reportados
por el Laboratorio Analizar Ltda.

1.16.5. Vertimiento secundario Centro de acopio lechería

Para el vertimiento secundario del centro de acopio de leche, se identificó un caudal

promedio de 1.188 L/s en una jornada de 24 horas, tomando alícuotas cada hora. El
flujo varía de 1.763 L/s a las 8:00pm hasta 0.916 L/s a las 11:00pm. Los mayores
flujos se reportan a tempranas horas de la madrugada hasta las 9:00am y
posteriormente en la tarde.

En cuento a la calidad del vertimiento se reportan valores de pH de 7.06, DQO de 77

mg O2 /L, para la DBO5 total es de 19 mg O2 /L, para el OD es de 4.44 mg O 2 /L y
para coliformes totales de 17 x 10 3 NMP / 100ml, entre otros parámetros reportados
por el Laboratorio Analizar Ltda.
1.16. 6. Aguas abajo de Vertimiento de centro de acopio de Leche

Las condiciones de flujo y de calidad definidas para el punto aguas abajo del sistema
o lo que se considera como el límite del sistema, identifica un caudal promedio de
207.17 L/s, el cual fue evaluado por el método de área velocidad, en una sección de
4.79 m de ancho, identificando 10 secciones.

En cuento a la calidad del vertimiento se reportan valores de pH de 7.27, DQO de 28

mg O2 /L, para la DBO5 total es de 2 mg O2 /L, para el OD es de 6.74 mg O 2 /L y para
coliformes totales de 22 x 10 2 NMP / 100ml, entre otros parámetros que fueron
evaluados por el laboratorio acreditado.

A continuación, se presenta la información relacionada con cada uno de los puntos

que definirán el límite del sistema trabajado, estableciendo la localización, las
coordenadas y la altura sobre el nivel del mar.

Tabla 12. Localización puntos a evaluar

DESCRIPCIÓN LOCALIZACIÓN COORDENADAS COORDENADAS ASNM
PLANAS GEOGRAFICAS
ESTE NORTE LATITUD LONGITUD
Aguas arriba
del primer
K0+00 1121486 1111774 5°36´21.39”N 75°58´51.73”W 2529
vertimiento
municipal
Vertimiento
K0+0.284 1121493 1111950 5°36´27.12”N 75°58´51.50”W 2527
Colegio
Vertimiento
principal K0+1.662 1121297 1113296 5°36´27.12”N 75°58´57.78”W 2531
municipio de Iza
Vertimiento
secundario
K0+2.387 1120952 1113821 5°37´28.05” N 72°58´57.78”W 2530
Centro de
Acopio Lechería
Aguas abajo del
Cetro de acopio K0+2.490 1120935 1113915 5°37´31.11” N 72°59´9.50”W 2503
Lechería
Fuente. INCO A&J S.A.S.

En la tabla se puede observar la información de calidad para el modelo, con soporte

en el trabajo de campo y laboratorio realizado por el laboratorio especializado.
 Tabla 13. Resultados caracterización de puntos de evaluación
Aguas Vertimiento Vertimiento Vertimiento Aguas
Arriba del Colegio principal secundario abajo del
primer municipio Centro de Cetro de
PARAMETRO UNIDADES Acopio acopio
Vertimiento de Iza
Lechería Lechería
municipal
TEMPERATURA (°C)
16 16 16 16 16
AMBIENTE **
TEMPERATURA (°C)
20.9 19.7 19.7 19.9 18.7
MUESTRA *
pH* Unidades
7.37 7.14 7.14 7.06 7.23
DQO* (mg O2 / L)
28 109 109 77 28
DBO5* (mg O2 / L)
<2 55 55 19 <2
DBO5 FAST** (mg O2 / L)
2 55 55 19 2
DBO5 SLOW ** (mg O2 / L)
26 54 54 58 26
OD* (mg O2 / L)
5.40 0.08 0.08 4.44 6.74
SST* (mg / L)
<11.51 14 14 19 <11.51
Conductividad * µSiemens/cm
150.2 429 429 211 187.4
Coliformes (NPM /100ml)
680 12 x 103 12 x 103 17 x 103 22 x 102
Totales*
Coliformes (NPM /100ml)
450 68 x 102 68 x 102 78 x 102 14 x 102
Fecales*
Fuente: Analizar Laboratorio Ltda.
*Analizar Laboratorio Ltda.
** El estudio.
1.17 MODELACIÓN PARA EL SISTEMA ACTUAL CON CAUDAL MEDIDO

Las aplicaciones preliminares son útiles para la identificación de datos importantes

que no existen y deficiencias teóricas no previstas. Estas corridas son utilizadas para
diseñar los trabajos y muestreos de campo y posteriormente se pueden calibrar para
lograr que los resultados y predicciones arrojados por el modelo estén lo bastante
cerca de los datos de calidad que se han obtenido previamente en registros de
campo.

La información de este ejercicio se presenta en la siguiente tabla y corresponde a los

datos WQ output o datos de calidad del agua de salida en el sistema. Para lograr la
estabilidad en el sistema modelado desde el contexto matemático, se definieron
cuatro tramos en una longitud total de 2.49 Km, definidos en los puntos donde se
recibe el vertimiento específico, es decir en el vertimiento del colegio, en el
vertimiento principal del municipio y en el vertimiento del centro de acopio de leche.
Los resultados del modelo permitieron observar un ajuste moderado para las
variables calculadas. Para efectos de la calibración, el modelo incluye un algoritmo
genético, por medio del cual se hacen los procesos de auto-calibrado, con el fin de
maximizar el ajuste en los datos registrados en las estaciones de monitoreo definidas
con respecto a los resultados del modelo.

El ajuste se realiza a partir del promedio de los coeficientes de regresión R 2 de las

diferencias entre las predicciones del modelo y los datos de calidad del agua
observados para los diferentes parámetros del agua.

Tabla 14. Presentación Modelo QUAL2KW6 – Río Iza

Fuente. INCO A&J S.A.S.

Para el caso del caudal se tomó en este primer ejercicio de modelación el caudal
medido correspondiente a 0.42 m 3/s, lo que permitió definir el comportamiento del
flujo a lo largo del sistema propuesto.

 Caudal

Los valores de caudal para el Rio Iza en el área de estudio fluctúan entre 520 L/s a
420 L/s, presentado un incremento debido a las descargas internas en el sistema
evaluado. En la tabla y gráfica, se puede apreciar las variaciones frente al
comportamiento del flujo, para la microcuenca analizada.
 Tabla 15. Reporte WQ OUTPUT para corrida preliminar con caudal medido
ELEMENTO Up Dist Down Dist Q T DBOu OD
X (Km) X(Km) ( m3/s) (°C) (mg/l) (mg/l)
Cabecera: Aguas Arriba del primer -
Vertimiento municipal 2.490 0.42 20.90 4.00 5.40
Vertimiento Colegio 2.348 1.517 0.43 20.72 3.87 5.31
Vertimiento principal municipio de 0.466
Iza 1.517 0.48 20.03 3.87 5.29
Vertimiento secundario Centro de 0.052
Acopio Lechería 0.466 0.51 19.74 4.21 5.30
Aguas abajo del Cetro de acopio 0.0
Lechería 0.052 0.52 19.70 4.27 5.29
Terminus: Fin del sistema - 0.0 19.70 4.27 5.29
Fuente. INCO A&J S.A.S.

Gráfica 1. Flujo en Rio Iza – Área de Influencia

Fuente. INCO A&J S.A.S.

 Temperatura

La temperatura del agua tiene una gran importancia en el desarrollo de los diversos
procesos que en ella se realizan, de forma que un aumento de la temperatura
modifica la solubilidad de las sustancias, aumentando la de los sólidos disueltos y
disminuyendo la de los gases21. Los resultados para la temperatura se pueden
apreciar en la siguiente gráfica. De acuerdo a los datos de temperatura, se puede

21
Jiménez, A. A. (2000). Determinación de los parámetros físico-químicos de calidad de las aguas.
Revista interdisciplinar de gestión ambiental, 2(23), 12-19.
 establecer una temperatura de 20.90°C en la cabecera hasta llegar a un valor de
19.70°C aguas abajo del vertimiento del centro de acopio de leche. Durante el
recorrido la disminución en la temperatura fue de 1.2°C para la longitud definida.

Gráfica 2. Resultados temperatura

Fuente. INCO A&J S.A.S.

 DBO

El modelo QUAL2Kw, requiere información para la DBO, diferenciándola en lo que

corresponde a la DBO rápida (fast) y a la DBO lenta (slow). Para el caso de la
descomposición lenta, la reacción se da por el incremento de la disolución de los
detritos, lo que genera una perdida a través del proceso de hidrolisis. En el caso
contrario, es decir para la DBO rápida, esta gana por vía del hidrolisis de la DBO lenta, la
pérdida se da por la oxidación y desnitrificación, finalmente.

Tabla 16. Resultados de la DBO

REACH X DBOslow DBOslow DBOslow DBOfast DBOfast DBOfast DBOu DBOu DBOu
(km) MIN (mg/l) MAX MIN (mg/l) MAX MIN (mg/l) MAX
(mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Headwate
r 2.490 2.00 2.00 2.00 2.00 2.00 2.00 4.00 4.00 4.00
1 2.348 1.97 1.97 1.97 1.89 1.89 1.89 3.87 3.87 3.87
2 1.517 1.97 1.97 1.97 1.90 1.90 1.90 3.87 3.87 3.87
3 0.5 2.14 2.14 2.14 2.07 2.07 2.07 4.21 4.21 4.21
4 0.05 2.17 2.17 2.17 2.10 2.10 2.10 4.27 4.27 4.27
Terminus 0.000 2.17 2.17 2.17 2.10 2.10 2.10 4.27 4.27 4.27
Fuente. INCO A&J S.A.S.
La demanda bioquímica de oxígeno se define de forma general como la cantidad de
oxígeno requerida por las bacterias en el proceso de estabilización de la materia
orgánica descomponible bajo condiciones aeróbicas.

Observando la información relacionada con la Demanda Bioquímica de Oxígeno

última o DBOu, se puede establecer una fluctuación de los valores entre 4.00 mg/L
en la cabecera, es decir en el kilómetro 2.49, para finalizar en el sistema con una
concentración de 4.27 mg/L. El perfil de DBO en la distancia establece un leve
aumento en 0.27 mg/L en este parámetro, situación generada por las condiciones de
reaireación en la sección del cauce.

Gráfica 3. Resultados de la DBO

Fuente. INCO A&J S.A.S.

La gráfica anterior, establece valores de DBO 22 para el Rio Iza por debajo de 5 mg/L,
lo que permite establecer una condición favorable para la corriente, debido a que se
considera que los valores de DBO mayores a 5 mg/L son para corrientes
contaminadas.

22
La DBO es el parámetro que más dificultades presenta para ser simulado, posiblemente por tratarse de un
bioensayo donde son muchos los factores que afectan la medición, como por ejemplo la relación de la materia
orgánica soluble a la suspendida, los sólidos sedimentables y los flotables, la presencia de hierro en su forma
oxidada o reducida, la presencia de compuestos azufrados entre otros (Corporación Autónoma Regional del
Cauca CAR, 2012).
  pH, SST, DQO, conductividad, Nitritos – Nitratos y Coliformes Totales

La siguiente tabla hace un reporte de los resultados del modelo en lo relacionado a

varios parámetros de interés para el análisis del sistema evaluado.

Tabla 17. Resultados para el pH, SST, DQO, Conductividad, NO3 y Coliformes
Totales
X (km) pH SST DQO Conductividad Nitrate + nitrite N Coliformes
(mg/l) (mg/L) (µS/cm 25C) Totales
(µgN/L) NMP/100ml

2.490 7.37 11.51 28.00 150.20 380.0 680

2.348 7.30 11.11 27.32 161.11 420.0 600
1.517 7.18 9.95 24.62 207.79 580.0 510
0.50 7.16 9.45 23.70 229.20 640.0 460
0.05 7.16 9.40 23.65 231.88 659.0 500
0.00 7.16 9.40 23.65 231.88 660.0 490
Fuente. INCO A&J S.A.S.

 DQO

La Demanda Química de Oxígeno o DQO presenta valores en un intervalo de 28

mg/l a 23.65 mg/l, con el máximo en la cabecera y el mínimo al final del sistema. A
través del trayecto se puede establecer que la DQO se reduce en 4.35 mg/l, situación
que favorece bastante el comportamiento del cauce y por su puesto identifica la
asimilación y dilución de los vertimientos que recibe.

Gráfica 4. Comportamiento de la DQO en la distancia

Fuente. INCO A&J S.A.S

  pH

Para el caso del pH, a través de la distancia se puede observar que los valores
fluctúan entre 7.37 a 7.16, mantiene la condición neutra a través de la distancia
evaluada.
Gráfica 5. Comportamiento del pH en la distancia

Fuente. INCO A&J S.A.S

 SST

La concentración de sólidos suspendidos totales para el monitoreo presenta valores

inferiores entre 11.51 mg/l hasta 9.40 mg/l. En la cabecera del sistema se puede
apreciar el mayor valor, disminuyendo gradualmente en 2.11 mg/L, hasta llegar a un
valor mínimo al final del sistema. La hidráulica de flujo en el cauce, así como el
efecto de dilución por la condición de caudal que se tiene en el inicio del sistema
evaluado, pueden ser las razones por las cuales se observa esta variación.
 Gráfica 6. Sólidos Suspendidos Inorgánicos

Fuente. INCO A&J S.A.S

 CONDUCTIVIDAD

La conductividad del agua es una expresión numérica de su habilidad para

transportar una corriente eléctrica. Depende de la concentración de sustancias
disueltas ionizables en el agua y de la temperatura a la cual se haga la
determinación. Por lo tanto, cualquier cambio en la cantidad de sustancias disuelta,
en la movilidad de los iones disueltos y en su valencia, implica un cambio en la
conductividad.

Los resultados obtenidos en el modelo establecen valores de conductividad que

oscilan entre 150.20 µS/cm 25ºC en la cabecera hasta 231.88 µS/cm 25ºC al final del
sistema. En el cauce se observa el incremento de este parámetro en 81.88 µS/cm
25ºC. Si se toma como referente el hecho de que la conductividad es una medida
indirecta de los sólidos disueltos, luego de multiplicarla por un factor que oscila entre
0.55 – 0.7, se puede observar el incremento de los sólidos disueltos en la longitud
evaluada.
 Gráfica 7. Conductividad en el sistema

Fuente. INCO A&J S.A.S

 NITRITOS Y NITRATOS

Los compuestos del nitrógeno se consideran de gran interés en el área de la

ingeniería de calidad de aguas. La química del nitrógeno es compleja, debido a los
varios estados de valencia efectuada que puede asumir este elemento y al hecho de
que los cambios en la valencia pueden ser efectuados por organismos vivos.
Además, se vuelven más interesantes si los cambios de valencia que realizan las
bacterias son positivos o negativos, según si las condiciones son anaerobias o
aerobias.

En aguas de polución reciente, contiene el nitrógeno en su mayoría como orgánico y

amoniacal, pero si el agua contiene en su mayoría nitratos, se considera que la
polución ocurrió en un amplio tiempo antes de que se realizara el análisis.

El reporte del modelo establece valores para la suma de nitritos y nitratos entre 380
µgN/L hasta 660 µgN/L, con los valores mayores al final del sistema, el
comportamiento es incremental y registra un valor de 280 µgN/L en los límites
definidos para el sistema.
 Gráfica 8. Nitritos y Nitratos en el sistema

Fuente. INCO A&J S.A.S

 Coliformes Totales

El modelo QUAL2Kw modela un patógeno genérico (coliformes fecales o coliformes

totales). La remoción de patógenos está determinada como función de la temperatura, la
luz y la sedimentación. Para el ejercicio se trabajó con coliformes totales, evidenciándose
una disminución en los valores. En la cabecera se reporta un valor de 680 NMP / 100ml,
disminuyendo gradualmente hasta el punto aguas abajo del sistema a 103 metros del
centro de acopio lechero, con un valor de 490 NMP / 100ml.

Gráfica 9. Coliformes Totales

Fuente. INCO A&J S.A.S

En la gráfica y la tabla siguiente, se presentan los resultados del oxígeno disuelto en el
Río Iza, en este caso, la línea azul representa la concentración de saturación, es decir el
máximo oxígeno disuelto que se puede alcanzar a lo largo de la corriente. La línea
continua representa los resultados de la simulación y los puntos los datos medidos en
campo.

Tabla [Link] de la modelación para el OD

X (km) T OD OD OD DO SAT DEFICIT
(°C) MIN (mg/l) MAX (mg/l) DE OD
(mg/l) (mg/l) (mg/l)
2.490 20.90 5.40 5.40 5.40 6.34 0.94
2.348 20.72 5.31 5.31 5.31 6.36 1.05
1.517 20.03 5.26 5.29 5.30 6.45 1.16
0.500 19.74 5.24 5.30 5.33 6.49 1.19
0.050 19.70 5.23 5.29 5.32 6.50 1.21
0.000 19.70 5.23 5.29 5.32 6.50 1.21
Fuente. INCO A&J S.A.S

En la modelación del oxígeno disuelto 23, se puede observar un valor de 5.40 mg/l en
la cabecera del sistema, avanzando a través de la distancia hacia un valor de 5.31
mg/l, para seguir avanzando hasta un valor de 5.29 mg/l al final del sistema.

En el reporte del modelo se puede establecer que la simulación plantea a través de

la longitud máximos de 5.40 mg/l y mínimos de 5.23 mg/l. El déficit de oxígeno
establecido a partir de la ecuación de Streeter and Phelps, en promedio sería de 1.13
mg/l, lo que establece una condición favorable, además es importante tener en
cuenta las condiciones morfométricas de las secciones del cauce, lo que beneficia la
hidráulica del cauce en la medida de los procesos de reaireación frente a la situación
de aportar o consumir oxigeno desde la base de la concentración de oxigeno de
saturación.
Gráfica 10. Perfil de OD en la distancia

23
La cinética del OD, establece su incrementa debido al proceso de fotosíntesis de las plantas, la biomasa de las algas de
fondo y el fitoplancton y por el contrario decrece debido al proceso de oxidación de la DBO rápida, por la respiración de las
plantas y por los procesos de nitrificación.
 Fuente. INCO A&J S.A.S

1.17.1 Zona de mezcla

La zona de mezcla es aquel volumen de agua en el cuerpo receptor en el que se

logra la dilución del vertimiento por procesos hidrodinámicos y de dispersión, sin
considerar otros factores además del decaimiento bacteriano, la sedimentación, la
asimilación en materia orgánica y la precipitación química. El propósito de la zona de
mezcla es asignar una región limitada para la mezcla completa del efluente con el
agua del cuerpo receptor, y así utilizar la capacidad de dilución.

Figura 15. Zona de mezcla en cuerpos naturales de agua

Fuente.[Link]
Para el cálculo de la extensión de la zona de mezcla aguas abajo del vertimiento se
puede determinar:
 - Asumiendo la mezcla instantánea
- Averiguar la longitud de mezcla

LONGITUD DE MEZCLA

Existen distintas ecuaciones propuestas para determinar esa longitud de mezcla, no

obstante, debido a la complejidad de las variaciones de profundidad en los ríos, su
distribución de velocidades y la tasa de mezcla, es difícil describir el fenómeno de
mezcla en fuentes hídricas naturales.

Para la determinación de la zona de mezcla aguas abajo del vertimiento se propone

el Método Simplificado desarrollado por la Agencia de Protección Ambiental de los
Estados Unidos de América (USEPA), que puede ser adoptado para descargas
superficiales en la orilla o en el centro del cuerpo de agua lotico.

La longitud necesaria para alcanzar la mezcla completa se determina con la siguiente

ecuación:

Donde:

= Ancho medio del cuerpo de agua, aguas abajo del vertimiento. (m)
U = Velocidad de flujo media del rio en la ubicación del vertimiento (m/s)
Dy = Coeficiente de dispersión lateral aguas abajo del vertimiento.

El coeficiente de dispersión se puede calcular con la siguiente ecuación:

Donde:

C = Coeficiente para cauces naturales con serpentear moderado (0.6)

d = Profundidad media del río aguas abajo del vertimiento (m)

=Velocidad de corte en metros por segundo (m/s).

Para determinar la velocidad se utiliza la siguiente ecuación:

 =

Dónde:

g = 9.81m/s2
d = Profundidad media del río aguas abajo del vertimiento (m)
S= Pendiente del cauce aguas abajo del vertimiento (m/m)

A continuación, se presenta la aplicación de esta metodología para cada uno de los

vertimientos evaluados en el sistema analizado para el Rio Iza.

 Vertimiento Colegio:

Para determinar la pendiente se tomó una distancia de 198 metros a partir del
vertimiento del colegio, para lo cual se obtiene una pendiente de 0.005 m/m. Los
datos tomados se registran a continuación:

COTA INICIAL COTA FINAL DISTANCIA PENDIENTE

2528m.s.n.m. 2527m.s.n.m. 198m 0.005

Con respecto a los insumos para la determinación de la zona de mezcla en relación

con la formulación planteada anteriormente, se toma como base una profundidad de
0.2827 metros, un ancho en la sección de 5.4 metros, velocidad de 0.1178 metros
por segundo, lo que reporta un coeficiente de dispersión de 0.01998 y una longitud
de zona de mezcla de 27.36 metro.

PROFUNDIDAD ANCHO VELOCIDAD COEFICIENTE DISPERSIÓN LONGITUD ZONA DE

d (m) W (m) Dy MEZCLA
(LzdM) (m)
(m/s)

0.2827 5.4 0.1178 0.01998 27.36

 Vertimiento principal del municipio de Iza

En este caso la longitud fue de 725 metros, por lo cual se obtiene una pendiente de
0.00138 m/m. Los datos tomados se registran a continuación:
COTA INICIAL COTA FINAL DISTANCIA PENDIENTE

2531m.s.n.m. 2530 m.s.n.m. 725m 0.00138

En este caso se tienen una profundidad promedio de 0.2827 metros, el ancho de 5.4
metros y velocidad de 0.062 metros por segundo, para un coeficiente de dispersión
de 0.01052 y la zona de mezcla de 27.35 metros.

PROFUNDIDAD ANCHO VELOCIDAD COEFICIENTE DISPERSIÓN LONGITUD ZONA DE

d (m) W (m) Dy MEZCLA
(LzdM) (m)
(m/s)

0.2827 5.4 0.062 0.01052 27.35

 Vertimiento centro de acopio de leche

En este caso la longitud fue de 103 metros, por lo cual se obtiene una pendiente de
0.1456 m/m. Los datos tomados se registran a continuación:

COTA INICIAL COTA FINAL DISTANCIA PENDIENTE

2530m.s.n.m. 2515 m.s.n.m. 103m 0.1456

La longitud de mezcla se define en 27.36 metros, con un coeficiente de dispersión de

0.1077 y velocidad de 0.625 metros por segundo.

PROFUNDIDAD ANCHO VELOCIDAD COEFICIENTE DISPERSIÓN LONGITUD ZONA DE

d (m) W (m) Dy MEZCLA
(LzdM) (m)
(m/s)

0.2827 5.4 0.625 0.1077 27.36

Para los tres vertimientos se puede establecer que la longitud de mezcla está
alrededor de 25 a 30 metros, manteniéndose un valor base, sin embargo, es claro
que se diferencian con respecto a cada una de las secciones evaluadas en lo
 relacionado a la pendiente, porque para el caso de la profundidad y ancho se
mantiene una sección promedio.

VERTIMIENTO DISTANCIA PENDIENTE VELOCIDAD COEFICIENTE LONGITUD

(m) (m/m) DISPERSIÓN ZONA MEZCLA
(m/s) Dy (LzdM) (m)

Vertimiento
Colegio 198m 0.005 0.1178 0.01998 27.36

Vertimiento
principal 725m 0.00138 0.062 0.01052 27.35
municipio
Vertimiento
centro de 103m 0.1456 0.625 0.1077 27.36
acopio leche

Para los tres casos se puede establecer que la zona de mezcla se encuentra dentro
de la longitud definida para cada uno de los tramos preestablecidos para el modelo.

1.18 ESCENARIOS DE SIMULACIÓN

Los escenarios que han sido planeados a nivel conceptual para lo relacionado al
saneamiento hídrico del Rio Iza y que han sido modelados con el objeto de evaluar
su impacto en la calidad de esta fuente hídrica y posteriormente apoyar la toma de
decisiones, corresponden a una base de caudales, que para el caso son el medido y
el mínimo, el primero preestablecido a través del trabajo de campo y el segundo a
partir de una disminución del medido en un 20%, como referencia frente al
comportamiento normal del cauce. Adicionalmente se plantea para cada uno de
estos flujos un escenario en el cual se tienen los vertimientos reales y normales de la
fuente, luego la fuente sin vertimientos y finalmente la fuente con la planta de
tratamiento.

ESCENARIO CAUDAL VERTIMIENTOS PTAR

I Medido Sin vertimientos -
II Medido Con vertimientos Con tratamiento
III Mínimo Con vertimientos -
IV Mínimo Sin vertimientos -
V Mínimo Con vertimientos Con tratamiento

1.18.1 ESCENARIO I: Caudal medido sin vertimientos:

Se plantea un primer escenario con el fin de poder tener una idea clara del
comportamiento del sistema del Rio Iza en la zona de evaluación sin el aporte de
cada uno de los vertimientos.
  Caudal

El registro de caudal medido fue de 0.42 m 3/s, en el inicio del sistema, fluctuando
entre valores de 0.43 m3/s, 0.48 m3/s y 0.51 m3/s. se observa un leve incremento a
través de la longitud, en 0.08 m 3/s, lo que puede ser definido por aportes naturales
y/o artificiales a través de la longitud del cauce.

En la tabla y gráfica siguiente se puede apreciar el comportamiento del flujo para el

sistema evaluado, en esta condición.

Tabla 19. Reporte WQ OUTPUT para escenario I Rio Iza

ELEMENTO Up Dist Down Dist Q T DBOu OD
X (Km) X(Km) ( m3/s) (°C) (mg/l) (mg/l)
Cabecera: 2.490 - 0.42 20.90 4.00 5.40
Tramo 1 2.206 0.828 0.43 20.54 3.87 5.31
Tramo 2 0.828 0.725 0.48 19.12 3.65 5.30
Tramo 3 0.725 0.103 0.51 18.52 3.55 5.34
Tramo 4 0.103 0.0 0.51 18.44 3.53 5.33
Terminus - 0.0 0.51 18.44 3.53 5.33
Fuente. INCO A&J S.A.S

Gráfica 11. Flujo en Rio Iza – Área de Influencia

Fuente. INCO A&J S.A.S.

 Temperatura

Los resultados para la temperatura se pueden apreciar en la siguiente gráfica. De

acuerdo a los datos de temperatura, se puede establecer una temperatura de
20.90°C en la cabecera hasta llegar a un valor de 18.44°C aguas abajo del
 vertimiento de centro de acopio de leche. Durante el recorrido el incremento 24 en la
temperatura es de 2.46°C para la longitud definida.

Gráfica 12. Resultados temperatura

Fuente. INCO A&J S.A.S.

 DBO

La demanda bioquímica de oxígeno se define de forma general como la cantidad de

oxígeno requerida por las bacterias en el proceso de estabilización de la materia
orgánica descomponible bajo condiciones aeróbicas.

Tabla 20. Resultados de la modelación para la DBO

REACH X DBOslow DBOslow DBOslow DBOfast DBOfast DBOfast DBOu DBOu DBOu
(km) MIN (mg/l) MAX MIN (mg/l) MAX MIN (mg/l) MAX
(mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Headwater 2.490 2.00 2.00 2.00 2.00 2.00 2.00 4.00 4.00 4.00
1 2.348 1.97 1.97 1.97 1.89 1.89 1.89 3.86 3.87 3.87
2 1.517 1.86 1.86 1.86 1.79 1.79 1.79 3.65 3.65 3.65
3 0.5 1.81 1.81 1.81 1.73 1.74 1.74 3.54 3.55 3.55
4 0.05 1.80 1.80 1.80 1.73 1.73 1.73 3.53 3.53 3.54
Terminus 0.000 1.80 1.80 1.80 1.73 1.73 1.73 3.53 3.53 3.54
Fuente. INCO A&J S.A.S.

Para el caso de la demanda bioquímica de oxígeno o DBOu, se puede establecer

valores máximos de 4.00 mg/l y mínimos de 3.53 mg/l. En el lapso de la distancia va
disminuyendo la concentración en 0.47 mg/l, si se toma como referencia la cabecera
y el final del sistema. La situación se puede generar por las condiciones de

24
La temperatura va aumentando progresivamente debido a la transferencia de calor de los elementos adyacentes, de la
transferencia de calor entre los sedimentos, el agua y la atmósfera, y de las entradas y salidas externas al sistema (Chapra
[Link]., 2008).
reaireación, favorecidas por las secciones del cauce y por otro lado se identifica el
proceso de dilución, debido a que en la cabecera se cuenta con una caudal de 0.42
m3/s, lo que favorece significativamente la condición de mezcla a través de la
longitud del sistema.

Gráfica 13. Resultados de la DBO

Fuente. INCO A&J S.A.S.

La gráfica anterior, establece valores de DBO para el Rio Iza, los cuales fluctúan
entre 4.0 mg/l, disminuyendo gradualmente a través de cada uno de los tramos
evaluados hasta llegar a 3.53 mg/l. Si se tiene como referencia que los valores por
encima de 5 mg/l corresponden a corrientes con algún grado de contaminación, se
puede decir que, para este escenario, la condición es favorable desde lo que
corresponde a la concentración de materia orgánica, en términos de DBOu.

 pH, SST, DQO y Coliformes Totales

La siguiente información presentada en el cuadro, corresponde a otros parámetros

modelados en este escenario, entre los cuales se tiene el pH, los sólidos
suspendidos totales, la DQO, la conductividad, los nitritos y nitratos y las coliformes
totales.

Tabla 21. Resultados para el pH, SST, DQO y Coliformes Totales

X (km) pH SST DQO Conductividad Nitrate + nitrite Coliformes
(mg/l) (mg/L) (µS/cm 25°C) N (µgN/L) Totales
NMP/100ml
2.490 7.37 11.51 28.00 150.20 0.38 47.96
 2.348 7.30 11.11 27.32 161.11 48.67 43.20
1.517 7.18 9.94 24.44 207.33 255.15 33.80
0.50 7.16 9.42 23.16 227.94 347.88 29.26
0.05 7.16 9.35 22.99 230.69 360.15 28.95
0.00 7.16 9.35 22.99 230.69 360.15 28.95
Fuente. INCO A&J S.A.S.

 pH

El reporte establece una variación de pH desde 7.37 en la cabecera variando en el

avance del cauce a 7.30, luego 7.18 y finalmente queda en 7.16. La disminución no
es muy significativa, alrededor de 0.21 unidades, sin embargo, la condición de pH en
todo e sistema es bastante buena.

 DQO:

En la tabla y la gráfica siguiente, se puede observar el comportamiento de la

Demanda Química de Oxigeno (DQO), la cual se asume en el modelo como un
constituyente genérico. Para el caso se pueden apreciar valores que fluctúan entre
28mg/L hasta 22.99mg/L, con una disminución gradual de 5.01 mg/L en toda la
extensión del sistema.

Gráfica 14. Resultados DQO

Fuente. INCO A&J S.A.S

 Sólidos Suspendidos Totales

La concentración de sólidos suspendidos totales para el monitoreo presentan valores
de 11.51 mg/l hasta 9.35 mg/l al final del sistema. La disminución se aprecia en el
recorrido por el cauce aguas abajo, con un valor de 2.16mg/l, desde la cabecera
hasta el punto aguas abajo del Rio Iza.

Gráfica 15. Sólidos Suspendidos Totales

Fuente. INCO A&J S.A.S

 Conductividad

Los resultados obtenidos en este escenario de simulación establecen valores de

conductividad que oscilan entre 150.20 µS/cm 25ºC en la cabecera hasta 230.69
µS/cm 25ºC al final del sistema. En el cauce se observa el incremento de este
parámetro en 81.49 µS/cm 25ºC.

Gráfica 16. Conductividad en el sistema

 Fuente. INCO A&J S.A.S

 NITRITOS Y NITRATOS

El reporte del modelo establece valores para la suma de nitritos y nitratos entre 0.38
µgN/L hasta 360.15 µgN/L, con los valores mayores al final del sistema, el
comportamiento es incremental y registra un valor de 359.77 µgN/L en los límites
definidos para el sistema.

Gráfica 17. Nitritos y Nitratos para el sistema evaludos

Fuente. INCO A&J S.A.S

 Coliformes totales
En la cabecera se reporta un valor de 47.96 NMP/100ml y al final de 28.95
NMP/100ml, lo que establece una disminución gradual para los 2.49 Km evaluados
con el modelo. La disminución se reporta en 19.01 NMP/ml.

Gráfica 18. Coliformes Totales

Fuente. INCO A&J S.A.S

 Oxígeno disuelto OD

En la gráfica y tabla siguiente, se presentan los resultados del oxígeno disuelto en el

Río Iza, en este caso, la línea azul representa el máximo oxígeno disuelto que se
puede alcanzar a lo largo de la corriente o lo que corresponde a la Concentración de
Saturación, con un valor inicial de 6.34 mg/L y luego a partir del kilómetro 2.348 de
6.36 mg/L, para finalizar con 6.50 mg/l. En todos los tramos evaluados para el
modelo, se puede establecer que la condición de oxígeno disuelto es óptima y con
muy buenos resultados.

Tabla [Link] de la modelación para el OD

X T OD OD OD DO SAT DEFICIT
(km) (°C) MIN (mg/l) MAX (mg/l) DE OD
(mg/l) (mg/l) (mg/l)
2.490 20.90 5.40 5.40 5.40 6.34 0.94
2.348 20.72 5.31 5.31 5.31 6.36 1.05
1.517 20.03 5.27 5.30 5.31 6.45 1.15
0.500 19.74 5.28 5.34 5.36 6.49 1.15
0.050 19.70 5.27 5.33 5.36 6.50 1.17
0.000 19.70 5.27 5.33 5.36 6.50 1.17
Fuente. INCO A&J S.A.S

El OD presenta valores de 5.40 mg/l en la cabecera del sistema hasta llegar a 5.36
mg/l al final del entorno de evaluación del escenario. Se observa una leve
 disminución en los 2.49 Km, lo que a su vez establece un déficit promedio de 1.105
mg/l.

Es importante observar las fronteras para este valor reportado del oxígeno,
estableciéndose un valor máximo y mínimo respectivamente en un rango entre 5.30
mg/L y 5.35 mg/L, condición bastante favorable para la respuesta del sistema hídrico.

Gráfica 19. Perfil de OD en la distancia

Fuente. INCO A&J S.A.S

1.18.2 ESCENARIO II: Caudal medido con sistema de tratamiento de aguas

residuales

Para el caso de este escenario se tiene en cuenta la ubicación de la planta de

tratamiento, por lo cual se amplía la longitud de la zona de estudio, atendiendo a que
la PTAR se localiza 180 metros más abajo del punto considerado como final del
sistema para las modelaciones anteriores. En ese sentido se hace una nueva
distribución de tramos, en una longitud de 2.87 Km desde la cabecera del sistema.

En la siguiente figura se puede establecer la ubicación de la planta de tratamiento de

aguas residuales para el municipio de Iza.

Figura 16. Localización vertimientos al Rio Iza

 Fuente. INCO A&J S.A.S – Google Earth

 Caudal

El registro de caudal medido fue de 0.42 m 3/s, en el inicio del sistema, fluctuando
entre valores de 0.52 m3/s y finalizando en 0.54 m3/s. Se puede apreciar un leve
incremento a través de la longitud de 0.12m 3/s, lo que puede ser definido por aportes
naturales y/o artificiales a través de la longitud del cauce.

En la tabla y gráfica siguiente se puede apreciar el comportamiento del flujo para el

sistema evaluado, en este escenario de simulación.
 Tabla 23. Reporte WQ OUTPUT para escenario I Rio Iza
ELEMENTO Up Dist Down Dist Q T DBOu OD
3
X (Km) X(Km) ( m /s) (°C) (mg/l) (mg/l)
Headwater 2.87 - 0.42 20.90 4.00 5.40
TRAMO 1 2.87 0.20 0.52 19.68 3.50 5.42
TRAMO 2 0.20 0.18 0.52 19.67 3.50 5.42
TRAMO 3 0.18 0.0 0.54 19.60 3.57 5.32
Terminus 0.0 0.54 19.60 3.57 5.32
Fuente. INCO A&J S.A.S

Gráfica 20. Flujo en Rio Iza – Área de Influencia

Fuente. INCO A&J S.A.S.

 Temperatura

Los resultados para la temperatura se pueden apreciar en la siguiente gráfica. De

acuerdo a los datos de temperatura, se puede establecer una temperatura de
20.90°C en la cabecera hasta llegar a un valor de 19.60°C hasta el punto final del
sistema. Durante el recorrido el incremento es de 1.3°C para la longitud definida.
 Gráfica 21. Resultados temperatura

Fuente. INCO A&J S.A.S.

 DBO

Para el caso de la demanda bioquímica de oxígeno o DBOu, se pueden establecer

valores máximos de 4.00 mg/l y mínimos de 3.57 mg/l. En el lapso de la distancia va
disminuyendo la concentración en 0.43 mg/l, si se toma como referencia la cabecera
y el final del sistema.

Tabla 24. Resultados de la modelación para la DBO

REACH X DBOslow DBOslow DBOslow DBOfast DBOfast DBOfast DBOu DBOu DBOu
(km) MIN (mg/l) MAX MIN (mg/l) MAX MIN (mg/l) MAX
(mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Headwate
r 2.490 2.00 2.00 2.00 2.00 2.00 2.00 4.00 4.00 4.00
1 2.348 1.78 1.79 1.79 1.69 1.70 1.70 3.48 3.50 3.51
2 1.517 1.79 1.79 1.79 1.69 1.70 1.70 3.48 3.50 3.51
3 0.5 1.84 1.84 1.84 1.71 1.71 1.71 3.55 3.57 3.58
4 0.05 1.84 1.84 1.84 1.71 1.71 1.71 3.55 3.57 3.58
Terminus 0.000 1.80 1.84 1.80 1.73 1.71 1.73 3.55 3.57 3.58
Fuente. INCO A&J S.A.S.

Los valores se limitan a máximos promedio de 3.63 mg/l y mínimos de 3.60 mg/l. Son
valores bajos que pueden ser asimilados a partir de las condiciones de calidad y en
especial de oxígeno disuelto de la fuente receptora.
 Gráfica 22. Resultados de la DBO

Fuente. INCO A&J S.A.S.

La gráfica permite observar el perfil de DBO, el cual presenta una disminución en la

distancia.

 pH, SST, DQO y Coliformes Totales

La siguiente información presentada en el cuadro, corresponde a otros parámetros

modelados en este escenario, entre los cuales se tiene el pH, los sólidos
suspendidos totales, la DQO, la conductividad, los nitritos y nitratos y las coliformes
totales.

Tabla 25. Resultados para el pH, SST, DQO y Coliformes Totales

X (km) pH SST DQO Conductividad Nitrate + nitrite Coliformes
(mg/l) (mg/L) (µS/cm 25°C) N (µgN/L) Totales
NMP/100ml
2.490 7.37 11.51 28.00 150.20 0.38 47.96
2.348 7.17 8.63 22.70 235.41 362.78 24.34
1.517 7.16 8.61 22.65 236.18 365.77 31.31
0.50 7.16 8.37 22.08 244.42 387.08 262.39
0.05 7.16 8.37 22.08 244.42 387.08 262.39
0.00 7.16 8.45 22.99 230.69 360.15 262.39
Fuente. INCO A&J S.A.S.

 pH
El reporte establece una variación de pH desde 7.37 en la cabecera variando en el
avance del cauce a 7.16. Se puede decir que para el sistema descrito se mantienen
las condiciones de neutralidad.

Gráfica 23. Perfil de DBO en la distancia

Fuente. INCO A&J S.A.S.

 DQO:

En la tabla y la gráfica siguiente, se puede observar el comportamiento de la

Demanda Química de Oxigeno (DQO), la cual se asume en el modelo como un
constituyente genérico. Para el caso se pueden apreciar valores que fluctúan entre
28mg/L hasta 22.99mg/L, con una disminución gradual de 5.01 mg/L, sin embargo,
se pueden observar algunos cambios leves al principio del sistema que implican
disminución en la concentración.

Gráfica 24. Perfil de DQO en la distancia

 Fuente. INCO A&J S.A.S

 Sólidos Suspendidos Totales

La concentración de sólidos suspendidos totales para el monitoreo presenta valores

de 11.51 mg/l al inicio, para ir disminuyendo hasta valores de 8.37 mg/l, sin embargo,
al final se incrementa a 8.45 mg/l, lo que significaría un aumento en 0.98mg/l.

Gráfica 25. Sólidos Suspendidos Totales

Fuente. INCO A&J S.A.S

 Conductividad

La conductividad oscila entre 150 y 230 µS/cm 25ºC, se puede decir que el promedio
en los 2.87 Km es de 223.55 µS/cm 25ºC, además se evidencia el incremento
gradual en el sentido del flujo.

Gráfica 26. Conductividad en el sistema evaluado

 Fuente. INCO A&J S.A.

 NITRITOS Y NITRATOS

El reporte del modelo establece valores para la suma de nitritos y nitratos en

promedio de 310.54 µgN/L, en todo caso en la distancia se observa el incremento en
este parámetro.

Gráfica 27. Nitritos y Nitratos para el sistema evaluado

Fuente. INCO A&J S.A.S

 Coliformes totales

Las coliformes totales presentan gran variación en a distancia, teniendo valores bajos
en la primera parte del sistema en promedio de 34.53 NMP/100ml y en la parte final
se evidencian valores de 262.39 NMP/100ml. Es claro que se beneficia en este
sector el crecimiento bacteriano, por lo cual el aumento es significativo.
 Gráfica 28. Coliformes Totales

Fuente. INCO A&J S.A.S

 Oxígeno disuelto OD

En la gráfica y tabla siguiente, se presentan los resultados del oxígeno disuelto en el

Río Iza. La concentración de saturación está alrededor de los 6.50 mg/l, por lo cual
los déficits calculados en promedio están alrededor de 1.12 mg/l. Se consideran
valores por encima de 5.0 mg/l bastante buenos y sobretodo frente a las condiciones
de asimilación de los vertimientos.

Se pueden analizar los valores de DBOu los cuales son bajos, garantizando las
condiciones de oxígenos disuelto en lo que se consideraría como una relación
inversa en su comportamiento, además se puede destacar el tema de las
condiciones hidráulicas que benefician en gran medida las condiciones de
reaireación de la fuente receptora.

Tabla [Link] de la modelación para el OD

X T OD OD OD DO SAT DEFICIT DE
(km) (°C) MIN (mg/l) MAX (mg/l) OD
(mg/l) (mg/l) (mg/l)
2.490 20.90 5.40 5.40 5.40 6.34 0.94
2.348 22.81 5.23 5.42 5.84 6.51 1.09
1.517 22.82 5.23 5.42 5.84 6.52 1.10
0.500 22.85 5.14 5.32 5.71 6.53 1.21
0.050 22.85 5.14 5.32 5.71 6.53 1.21
0.000 19.70 5.27 5.33 5.36 6.50 1.17
Fuente. INCO A&J S.A.S

Las condiciones de frontera definen un máximo de OD de 5.64 mg/l y mínimo de 5.24

mg/l. El máximo déficit de oxígeno calculado es de 1.21 mg/l, lo que favorece en gran
medida las condiciones de calidad de la fuente.

Gráfica 29. Perfil de OD en la distancia

 Fuente. INCO A&J S.A.S
1.18.3 ESCENARIO III: Caudal mínimo con los tres vertimientos

El escenario define la modelación con el caudal mínimo establecido en 0.336 m 3/s, el

cual disminuye en un 20% con respecto al caudal que fue medido en el trabajo de
campo. Esta situación se convertiría en una situación crítica del sistema para poder
evaluar el comportamiento al recibir las tres descargas actuales en la fuente
receptora.

 Caudal

En este escenario la modelación entre 0.34 m3/ s y 0.43 m3/s, lo que identifica el
incremento del flujo por los aportes de los vertimientos. Se puede establecer que en
el sistema el aumento es de 0.09 m 3/s. En la tabla y gráfica siguiente se puede
apreciar el comportamiento del flujo para el sistema evaluado, en esta condición.

Tabla 27. Reporte WQ OUTPUT para escenario I Rio Iza

ELEMENTO Up Dist Down Dist Q T DBOu OD
3
X (Km) X(Km) ( m /s) (°C) (mg/l) (mg/l)
Cabecera: 2.490 - 0.34 20.90 4.00 5.40
Tramo 1 2.206 0.828 0.35 20.68 3.85 5.30
Tramo 2 0.828 0.725 0.40 19.87 3.86 5.26
Tramo 3 0.725 0.103 0.43 19.53 4.27 5.27
Tramo 4 0.103 0.0 0.43 19.49 4.34 5.26
Terminus - 0.0 0.43 19.49 4.34 5.26
Fuente. INCO A&J S.A.S

Gráfica 30. Flujo en Rio Iza – Área de Influencia

 Fuente. INCO A&J S.A.S.

 Temperatura

Los resultados para la temperatura se pueden apreciar en la siguiente gráfica. De

acuerdo a los datos de temperatura, se puede establecer una temperatura de 20.90°C en
la cabecera hasta llegar a un valor de 19.49°C en el punto aguas abajo definido como
límite final del sistema. La disminución en la distancia es de 1.41 grados.

Gráfica 31. Resultados temperatura

Fuente. INCO A&J S.A.S.

 DBO
 La demanda bioquímica de oxígeno se define de forma general como la cantidad de
oxígeno requerida por las bacterias en el proceso de estabilización de la materia
orgánica descomponible bajo condiciones aeróbicas.

Tabla 28. Resultados de la DBO

REACH X DBOslow DBOslow DBOslow DBOfast DBOfast DBOfast DBOu DBOu DBOu
(km) MIN (mg/l) MAX MIN (mg/l) MAX MIN (mg/l) MAX
(mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Headwate
r 2.490 2.00 2.00 2.00 2.00 2.00 2.00 4.00 4.00 4.00
1 2.348 1.97 1.97 1.97 1.88 1.88 1.88 3.86 3.85 3.87
2 1.517 1.96 1.97 1.97 1.89 1.89 1.89 3.65 3.86 3.65
3 0.500 2.17 2.17 2.17 2.10 2.10 2.10 3.54 4.27 3.55
4 0.050 2.21 2.21 2.21 2.13 2.14 2.14 3.53 4.34 3.54
Terminus 0.000 2.21 2.21 2.21 2.13 2.14 2.14 3.53 4.34 3.54
Fuente. INCO A&J S.A.S.

Observando la información relacionada con la DBOu, se puede establecer una

fluctuación de los valores entre 4.0 mg/L en la cabecera, es decir en el kilómetro
2.49, para finalizar en el sistema con una concentración de 4.34 mg/L. El perfil de
DBO en la distancia establece un aumento de 0.34 mg/L, situación que se puede
presentar con la disminución en los flujos, disminuyendo la condición de dilución a
través del tiempo de viaje.

Gráfica 32. Resultados de la DBO

Fuente. INCO A&J S.A.S.

La gráfica anterior, establece valores de DBO para el Rio Iza, con un incremento a partr
de kilómetro 1.5 aumentando gradualmente hasta alcanzar un valor de 4.34 mg/l. Si se
tiene como referencia que los valores por encima de 5 mg/l corresponden a corrientes
con algún grado de contaminación, se puede decir que, para este escenario, la condición
no sería muy favorable desde lo que corresponde a la concentración de materia
orgánica, en términos de DBOu.

 pH, SST, DQO, Conductividad, Nitritos- Nitratos y Coliformes Totales

La tabla muestra el reporte de otros parámetros como resultado de la modelación

para este escenario.

Tabla 29. Resultados para el pH, SST, DQO y Coliformes Totales

X (km) pH SST DQO Conductividad Nitrate + nitrite N Coliformes
(mg/l) (mg/L) (µS/cm 25C) (µgN/L) Totales
NMP/100ml
2.490 7.37 11.51 28.00 150.20 380.0 680
2.348 7.28 11.02 27.16 163.76 420.0 570
1.517 7.16 9.62 23.91 219.94 620.0 490
0.500 7.15 9.04 22.85 244.75 720.0 460
0.050 7.14 8.99 22.80 247.77 730.0 590
0.000 7.14 8.99 22.80 247.77 730.0 590
Fuente. INCO A&J S.A.S.
 DQO

En la tabla y la gráfica siguiente, representa los valores relacionados con la

Demanda Química de Oxigeno (DQO), asumida en el modelo como un constituyente
genérico. Para el caso se pueden apreciar valores que fluctúan entre 28 mg/L hasta
22.80 mg/L. En la distancia se puede apreciar la disminución gradual de este
parámetro, la cual se puede generalizar a 5.2 mg/l en toda la distancia del sistema
evaluado.

Gráfica 33. Perfil de DQO en la distancia

Fuente. INCO A&J S.A.S

  pH

Para el caso del pH, a través de la distancia se puede observar que los valores
fluctúan entre 7.37 hasta 7.14, sin embargo, se mantiene la condición básica a través
de los 2.49 kilómetros. Es evidente la disminución en este parámetro, sin embargo,
se mantienen en los valores permisibles.
Gráfica 34. Comportamiento del pH en el sistema evaluado

Fuente. INCO A&J S.A.S

 Sólidos Suspendidos Totales

La concentración de sólidos suspendidos totales para el monitoreo presenta valores

inferiores a 15 mg/L. En la cabecera del sistema se puede apreciar un valor de 11.51
mg/L, disminuyendo gradualmente hasta 8.99 mg/L. Lo anterior define un diferencial
de 2.52 mg/l en la distancia. Es bastante interesante la disminución gradual en este
parámetro, sobre todo por las condiciones del cauce y en relación a que se modela
con el caudal mínimo.

Gráfica 35. Sólidos Suspendidos Inorgánicos

 Fuente. INCO A&J S.A.S

 Conductividad

La conductividad oscila entre 150 y 250 µS/cm 25ºC, se puede decir que el promedio
en los 2.49 Km es de 212.3 µS/cm 25ºC, además se evidencia el incremento gradual
en el sentido del flujo.

Gráfica 36. Conductividad en el sistema

Fuente. INCO A&J S.A.S

 NITRITOS Y NITRATOS
El reporte del modelo establece valores para la suma de nitritos y nitratos en
promedio de 600 µgN/L, con valores máximos entre 730 µgN/L al final del sistema y
mínimos de 380 µgN/L en la cabecera.

Gráfica 37. Nitritos y Nitratos para el sistema evaluado

Fuente. INCO A&J S.A.S

 Coliformes totales

En la cabecera se reporta un valor de 680 NMP / 100ml, disminuyendo gradualmente

hasta 460 NMP/100ml en el kilómetro 0.50 y en la parte final del sistema nuevamente se
incrementa a 590 NMP/100ml.

Gráfica 38. Coliformes Totales

Fuente. INCO A&J S.A.S

 Oxígeno disuelto OD
 En la gráfica y la tabla siguiente, se presentan los resultados del oxígeno disuelto en el
Río Iza. Los valores de oxígeno disuelto máximo representados por la concentración de
saturación se encuentran en valores que fluctúan entre 6.34 mg/l y 6.53 mg/l.

Tabla [Link] de OD
X T OD OD OD DO DEFICIT
(km) (°C) MIN (mg/l) MAX SAT DE OD
(mg/l) (mg/l) (mg/l) (mg/l)
2.490 20.90 5.40 5.40 5.40 6.34 0.94
2.348 20.68 5.30 5.30 5.30 6.37 1.07
1.517 19.87 5.22 5.26 5.28 6.47 1.21
0.500 19.53 5.20 5.27 5.31 6.52 1.25
0.050 19.49 5.19 5.26 5.30 6.53 1.27
0.000 19.49 5.19 5.26 5.30 6.53 1.27
Fuente. INCO A&J S.A.S

El resultado de la modelación permite establecer una condición bastante favorable para

el oxígeno a través de la distancia, con valores de 5.40 mg/l en la cabecera hasta
5.26mg/l en el límite final. Por lo anterior el déficit de oxígeno es bastante bajo del orden
de 1.2 mg/l.

Gráfica 39. Perfil de OD en la distancia

Fuente. INCO A&J S.A.S

1.18.4 ESCENARIO IV: Caudal mínimo sin vertimiento

Este escenario define un ejercicio de modelación e situación crítica en donde el

caudal es el mínimo para la fuente y además se excluyen los vertimientos principales
y secundarios del sistema.
  Caudal

Los valores de caudal para el Rio Iza en el área de estudio fluctúan entre 0.34 m 3/ s
hasta 0.43 m3/s presentando un incremento en la distancia de 0.09 m 3/s, el cual
puede ser generado por aportes naturales o artificiales en la corriente en su viaje
hasta el punto aguas abajo. En la tabla y gráfica siguiente se puede apreciar el
comportamiento del flujo para el sistema evaluado.

Tabla 31. Reporte WQ OUTPUT para escenario IV- Rio Iza

ELEMENTO Up Dist Down Dist Q T DBOu OD
X (Km) X(Km) ( m3/s) (°C) (mg/l) (mg/l)
Cabecera: 2.490 - 0.34 20.90 4.00 5.40
Tramo 1 2.206 0.828 0.35 20.68 4.02 5.70
Tramo 2 0.828 0.725 0.40 19.87 3.74 5.56
Tramo 3 0.725 0.103 0.42 19.53 3.61 5.53
Tramo 4 0.103 0.0 0.43 19.49 3.60 5.52
Terminus - 0.0 0.43 19.49 4.34 5.26
Fuente. INCO A&J S.A.S

Gráfica 40. Flujo en Rio Iza – Área de Influencia

Fuente. INCO A&J S.A.S

 Temperatura

Los resultados para la temperatura se pueden apreciar en la siguiente gráfica. De

acuerdo a los datos de temperatura, se puede establecer una temperatura de 20.90°C en
 la cabecera hasta llegar a un valor de 19.49°C aguas abajo. Se puede establecer una
disminución de 1.41 grados aproximadamente en todo el sistema.

Gráfica 41. Resultados temperatura

Fuente. INCO A&J S.A.S.

 DBO

La demanda bioquímica de oxígeno se define de forma general como la cantidad de

oxígeno requerida por las bacterias en el proceso de estabilización de la materia
orgánica descomponible bajo condiciones aeróbicas.

Tabla 32. Resultados de la DBO

REACH X DBOslow DBOslow DBOslow DBOfast DBOfast DBOfast DBOu DBOu DBOu
(km) MIN (mg/l) MAX MIN (mg/l) MAX MIN (mg/l) MAX
(mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Headwate
r 2.490 2.00 2.00 2.00 2.00 2.00 2.00 4.00 4.00 4.00
1 2.348 1.97 1.97 1.97 1.87 1.87 1.87 3.95 4.02 4.08
2 1.517 1.83 1.83 1.83 1.75 1.75 1.75 3.68 3.74 3.79
3 0.500 1.77 1.77 1.78 1.69 1.69 1.69 3.55 3.61 3.66
4 0.050 1.77 1.77 1.77 1.68 1.69 1.69 3.54 3.60 3.65
Terminus 0.000 1.77 1.77 1.77 1.68 1.69 1.69 3.54 3.60 3.65
Fuente. INCO A&J S.A.S.

Los valores de DBOu, para este escenario de simulación está entre 4.00 mg/l hasta
3.60 mg/l, con una disminución gradual a través de la distancia. La variación no es
tan significativa, por lo cual se puede definir como un valor constante a lo largo de la
microcuenca del Rio Iza evaluada.
 Gráfica 42. Resultados de la DBO

Fuente. INCO A&J S.A.S.

 pH, SST, DQO, Conductividad, Nitratos- Nitritos y Coliformes Totales

La siguiente tabla presenta los valores para cada uno de los parámetros definidos y
se complementan con las gráficas reportadas en el ejercicio de modelación.

Tabla 33. Resultados para el pH, SST, DQO y Coliformes Totales

X (km) pH SST DQO Conductividad Nitrate + nitrite N Coliformes
(mg/l) (mg/L) (µS/cm 25°C) (µgN/L) Totales
NMP/100ml
2.490 7.37 11.51 28.00 150.20 380.00 47.96
2.348 7.36 11.18 27.16 163.76 425.79 42.40
1.517 7.18 9.76 23.69 219.41 634.31 32.13
0.50 7.16 9.14 22.20 243.35 726.70 27.42
0.05 7.16 9.06 22.00 246.51 738.20 27.08
0.00 7.16 9.06 22.00 246.51 738.20 27.08
Fuente. INCO A&J S.A.S.

 DQO:

Los valores para la DQO están definidos entre 28 mg/l y 22 mg/l. El mayor valor se
encentra en la cabecera como parte de la valoración del sistema, sin embargo a
través de la distancia va disminuyendo hasta alcanzar el valor de 22 mg/l.

Gráfica 43. Resultados DQO

 Fuente. INCO A&J S.A.S

 pH

Para el caso del pH, a través de la distancia se puede observar que los valores
fluctúan entre 7.37 hasta 7.16, sin embargo, se mantiene la condición neutra en la
longitud a evaluar.

Gráfica 44. Comportamiento del pH en la distancia

Fuente. INCO A&J S.A.S

 Sólidos Suspendidos Totales

La concentración de sólidos suspendidos totales para la modelación de este escenario
reporta valores mínimos de 9.06 mg/l en el final del sistema y máximos de 11.51 mg/l en
la cabecera.

Gráfica 45. Sólidos Suspendidos Totales

Fuente. INCO A&J S.A.S

 CONDUCTIVIDAD

Los valores de conductividad están alrededor de 150 a 246 µS/cm 25°C, en

este proceso de modelación se evidencia el incremento en el sentido del flujo,
con incrementos graduales a través de la distancia.

Gráfica 46. Conductividad en el sistema

Fuente. INCO A&J S.A.S

  NITRITOS Y NITRATOS:

Para el caso de los nitritos y nitratos se pueden observar valores promedio de 607.2
µgN/L, con máximos de 738.20 µgN/L al final del sistema y mínimos de 380 µgN/L.

Gráfica 47. Nitritos y Nitratos en el sistema

Fuente. INCO A&J S.A.S

 Coliformes totales

En el ejercicio de modelación se presentan valores de 47.96 NMP/100ml en la

cabecera y a través de la distancia disminuye hasta un valor de 27.08 NMP/ml. Es
posible que no se incrementen las coliformes por que en la zona no se beneficia el
crecimiento de las mismas.
Gráfica 48. Coliformes Totales

Fuente. INCO A&J S.A.S

 Oxígeno disuelto OD
Para este escenario el OD fluctúa en valores mayores a 5.0 mg/l, con máximos de
5.70 mg/l y mínimos de 5.52 mg/l. El incremento se va dando en la distancia y se
definen unos límites para este parámetro entre 5.01mg/l como mínimo y un máximo
de 6.3 mg/l.

Teniendo como base la concentración de saturación promedio de 6.58 mg/l, los

déficits de oxígeno en toda la extensión del sistema evaluado oscilan alrededor de
1.07 mg/l. Condición muy favorable para la dinámica del sistema.

Tabla [Link] de resultados para el OD

X T OD OD OD DO SAT DEFICIT
(km) (°C) MIN (mg/l) MAX (mg/l) DE OD
(mg/l) (mg/l) (mg/l)
2.490 20.90 5.40 5.40 5.40 6.34 0.94
2.348 20.68 4.81 5.70 7.01 6.40 0.70
1.517 19.87 4.89 5.56 6.52 6.61 1.05
0.500 19.53 4.99 5.53 6.30 6.70 1.17
0.050 19.49 4.99 5.52 6.27 6.72 1.20
0.000 19.49 4.99 5.52 6.27 6.72 1.20
Fuente. INCO A&J S.A.S

La siguiente gráfica presenta el perfil de oxígeno disuelto en la distancia en una

longitud evaluada de 2.49 kilómetros, teniendo como referencia la concentración de
saturación, para poder establecer los faltantes de oxígenos en cada tramo.

Gráfica 49. Perfil de OD en la distancia

Fuente. INCO A&J S.A.S

1.18.5 ESCENARIO V: Modelación con caudal mínimo y sistema de tratamiento
de aguas residuales

Para el caso de este escenario se tiene en cuenta la ubicación de la planta de

tratamiento, por lo cual se amplía la longitud de la zona de estudio, atendiendo a que
la PTAR se localiza 180 metros más abajo del punto considerado como final del
sistema para las modelaciones anteriores. En ese sentido se hace una nueva
distribución de tramos, en una longitud de 2.87 Km desde la cabecera del sistema.

En la siguiente figura se puede establecer la ubicación de la planta de tratamiento de

aguas residuales para el municipio de Iza.

Figura 17. Localización vertimientos al Rio Iza

Fuente. INCO A&J S.A.S – Google Earth

 Caudal

El caudal presenta valores de 0.34 m 3/s y se incrementa a través de la longitud en

0.45 m3/s. Se estaría hablando de una diferencia de 0.11 m 3/ s, en toda la extensión
del sistema evaluado. Es importante aclarar que el escenario se trabaja con caudal
mínimo.

Tabla 35. Reporte WQ OUTPUT para escenario V - Rio Iza

ELEMENTO Up Dist Down Dist Q T DBOu OD
3
X (Km) X(Km) ( m /s) (°C) (mg/l) (mg/l)
Headwater 2.87 - 0.34 20.90 4.00 5.40
TRAMO 1 2.87 0.20 0.43 19.47 3.45 5.71
TRAMO 2 0.20 0.18 0.43 19.46 3.45 5.70
TRAMO 3 0.18 0.0 0.45 19.40 3.53 5.56
Terminus 0.0 0.45 19.40 3.53 5.56
Fuente. INCO A&J S.A.S

Gráfica 50. Flujo en Rio Iza – Área de Influencia

 Fuente. INCO A&J S.A.S.

 Temperatura

La temperatura promedio en los 2.87 kilómetros, está en 19.73°C, con valores máximos
de 20.90°C y mínimos de 19.40°C. En la distancia se aprecia una disminución de la
temperatura por las condiciones del mismo balance de calor que se realiza en el sistema
hídrico.

Gráfica 51. Resultados temperatura

Fuente. INCO A&J S.A.S.

  DBO

Los reportes de DBO en la modelación de este escenario definen un máximo de 4.00

mg/l e la cabecera y un mínimo de 3.45 mg/l, en el tramo 1 y 2, para incrementarse
nuevamente en 3.53 mg/l en la parte final del sistema.

Tabla 36. Resultados de la DBO

REACH X DBOslow DBOslow DBOslow DBOfast DBOfast DBOfast DBOu DBOu DBOu
(km) MIN (mg/l) MAX MIN (mg/l) MAX MIN (mg/l) MAX
(mg/l) (mg/l) (mg/l) (mg/l) (mg/l) (mg/l)
Headwater 2.87 2.00 2.00 2.00 2.00 2.00 2.00 3.00 4.00 4.00
1 1.54 1.75 1.75 1.75 1.65 1.66 1.66 2.40 3.45 3.49
2 0.19 1.75 1.75 1.75 1.65 1.66 1.66 3.40 3.45 3.49
3 0.09 1.81 1.81 1.82 1.67 1.68 1.68 3.49 3.53 3.58
Terminus 0.00 1.81 1.81 1.82 1.67 1.68 1.68 3.49 3.53 3.58
Fuente. INCO A&J S.A.S.

El sistema plantea límites de frontera promedio entre 3.55 mg/l a 3.63 mg/l, situación
que no demuestra mucha fluctuación en los valores, además que son bajos para las
condiciones analizadas. En la cabecera se tienen 4.0 mg/l y en el final del sistema se
definen valores entre 3.49 mg/l a 3.58 mg/l.

La gráfica permite observar el perfil de DBOu a través de la distancia definida para la

evaluación y que corresponde a 2.87 Km.

Gráfica 52. Resultados de la DBO

Fuente. INCO A&J S.A.S.

  pH, SST, DQO, Conductividad, Nitritos – Nitratos y Coliformes Totales

Para los parámetros enunciados, se presenta la siguiente tabla que permite

establecer el comportamiento a través d ela distancia de cada uno de ellos.

Tabla 37. Resultados para el pH, SST, DQO y Coliformes Totales

X (km) pH SST DQO Conductividad Nitrate + nitrite N Coliformes
(mg/l) (mg/L) (µS/cm 25°C) Totales
(µgN/L)
NMP/100ml

2.87 7.37 11.51 28.00 150.20 380.0 47.96

1.54 7.16 8.14 21.67 251.90 422.95 22.11
0.19 7.16 8.12 21.61 252.87 426.62 33.62
0.09 7.16 7.86 20.98 261.91 449.35 303.86
0.00 7.16 7.86 20.98 261.91 449.35 303.86
Fuente. INCO A&J S.A.S.

 DQO

Los valores de la Demanda Química de Oxigeno (DQO), para el caso se presentan

en promedio con un valor de 24.49mg/l, con un máximo de 28 mg/l en la cabecera y
un mínimo de 20.98 mg/l al final del sistema, es decir aguas abajo 0.18 Km de la
PTAR proyectada.

Gráfica 53. Perfil de DQO en la distancia

Fuente. INCO A&J S.A.S

  pH

Se mantiene en toda su extensión con un valor promedio de 7.27. Se observa una

leve disminución, aunque manteniendo las condiciones deseadas.

Gráfica 54. pH en la distancia

Fuente. INCO A&J S.A.S

 Sólidos Suspendidos Totales

La concentración de sólidos suspendidos totales para el monitoreo presenta valores

inferiores a 12 mg/L. En la cabecera del sistema se puede apreciar un valor de 11.51
mg/L, aumentando gradualmente en 3.65 mg/L, hasta llegar a un valor de 7.86 mg/L.

Gráfica 55. Sólidos Suspendidos Inorgánicos

 Fuente. INCO A&J S.A.S

 CONDUCTIVIDAD

Los valores de conductividad reportados para el modelo en este escenario presentan

máximos de 261.91 µS/cm 25°C en el final del sistema y mínimos de 150.20 µS/cm
25°C. Se establece un incremento a través de la distancia.

Gráfica 56. Conductividad reportada en la distancia

Fuente. INCO A&J S.A.S

 NITRITOS Y NITRATOS:
Los valores reportados en la modelación para el caso de los nitritos y nitratos
fluctúan entre 380 µgN/L y 449,35 µgN/L, con los valores máximos al final del
sistema evaluado. El incremento se acerca a 69.35 µgN/L, es decir unos 0.06935
mg/l.

Gráfica 57. Nitritos y Nitritos en la distancia

Fuente. INCO A&J S.A.S

 Coliformes totales

Las coliformes fecales en este escenario de simulación varían entre 47.96 NMP/100ml
en la cabecera y se incrementan a través de la distancia hasta alcanzar un valor máximo
de 303.86 NMP/100ml. El ejercicio plantea condiciones favorables para el crecimiento
bacteriano.
Gráfica 58. Coliformes Totales

Fuente. INCO A&J S.A.S

 Oxígeno disuelto OD
Las concentraciones de OD en la distancia definen una condición muy favorable para la
dinámica del cuerpo receptor. Se puede evidenciar la relación inversa de este parámetro
con respecto a la DBOu. Los valores presentan un promedio de 5.56 mg/l, con un
máximo de 5.71 mg/l y un mínimo de 5.40 mg/l en la cabecera. En la gráfica y la tabla
siguiente, se presentan los resultados del oxígeno disuelto en el Río Iza.

Tabla [Link] de OD
X (km) T OD OD OD DO SAT DEFICIT OD
(°C) MIN (mg/l) MAX (mg/l) (mg/l)
(mg/l) (mg/l)
2.87 20.90 5.40 5.40 5.40 6.34 0.94
1.54 19.47 5.09 5.71 6.77 6.54 0.83
0.19 19.46 5.09 5.70 6.76 6.55 0.85
0.09 19.40 4.99 5.56 6.54 6.56 1.00
0.00 19.40 4.99 5.56 6.54 6.56 1.00
Fuente. INCO A&J S.A.S

Las concentraciones de saturación presentan mínimos de 6.34 mg/l en el inicio y se

va incrementado gradualmente hasta llegar a un valor de 6.56 mg/l. Es bastante
favorable el comportamiento hidráulico para beneficiar la oxigenación del cauce
receptor. Por lo anterior, los déficits de oxígenos registrados a través de Streeter
and Phleps definen un mínimo de 0.83 mg/l y un máximo de 1.00mg/l.

Gráfica 59. Perfil de OD en la distancia

Fuente. INCO A&J S.A.S

1.19 ANÁLISIS FINAL DEL EJERCICIO DE MODELACIÓN EN EL SISTEMA
DEFINIDO PARA EL RIO IZA

El ejercicio de modelación se realizó con el modelo QUAL2Kw 6.0, realizando el

ingreso de datos a partir de los registros de campo y laboratorio realizados por el
laboratorio contratado, además se complementó alguna información con informes
oficiales de la cuenca del río y que son de consulta pública, lo que a su vez permitió
hacer los ajustes complementarios al proceso de calibración.

En un contexto general se plantea un sistema con 2.49 kilómetros en los cuales se

presenta la descarga del Colegio, el vertimiento principal del municipio de Iza y el
vertimiento del Centro de acopio de leche, los cuales se van mezclando en la
longitud hasta llegar a la parte final del sistema descrito, es decir 103 metros
después del vertimiento del centro de acopio de leche.

El ejercicio de modelación definió valores de DBOu y OD en la cabecera del sistema

del orden de 4.0 mg/l y 5.40 mg/l < de 2.0 mg/L respectivamente, sumado a un
caudal de 00.42 m3/s, lo que permite establecer una óptima condición de oxígeno,
además de la condición de reaireación favorecida por la morfometría del cauce. En el
caso del flujo, la cantidad establece la efectividad en los procesos de dilución frente a
las descargas de materia orgánica que se puedan aportar en la distancia analizada.

La evaluación en el punto aguas abajo, permite identificar un incremento en el

caudal, acercándose al orden de los 0.52 m 3/s, con una DBOu de 4.27 mg/l y un OD
de 5.29 mg/l. Para el caso la condición es bastante buena y la información define las
mejores capacidades del sistema hídrico para asimilar los aportes de los tributarios
evaluados.

Los resultados de laboratorio definieron una descarga de DBO 5 de menos de 2.0 mg/l
en la cabecera, 55 mg/l en el vertimiento principal del municipio y de 19mg/l en el
vertimiento del centro de acopio de leche. En lo relacionado al OD los valores aguas
arriba oscilan entre 5.40 mg/l y para el caso de los vertimientos en el principal del
municipio fue de 0.08 mg/l y en el vertimiento del centro de acopio fue de 4.44 mg/l.

Para efectos del modelo, se definieron 4 tramos para la modelación incial con los
reportes en campo y laboratorio, lo que definió una longitud de 2.49 kilómetros.
Posteriormente para los escenarios con el sistema de tratamiento, fue necesario
ampliar dicha longitud, debido a la ubicación de la PTAR, por lo cual la evaluación se
realizó en 2.87 kilómetros. Lo anterior permitió realizar los respectivos balances a
nivel de flujo, calidad y temperatura, sin embargo, en los reportes gráficos se pueden
observar los resultados en un contexto de mayor detalle.

Los resultados de la DBOu se pueden apreciar en el perfil, con un valor de 4.00 mg/l
en la cabecera y de 4.27 mg/l al final del sistema. Si se tiene en cuenta que este
parámetro evalúa la cantidad de oxígeno absorbido por un residuo en su proceso de
descomposición, se puede establecer que la corriente ejerce un proceso en las
mejores condiciones, lo que se evidencia a través de los resultados de OD
posteriores, los cuales están entre 5.40 mg/l hasta 5.29 mg/l. Para el caso del pH se
presentan valores alrededor de 7.37 a 7.16 y los sólidos suspendidos totales fluctúan
entre 11.51 mg/l en la cabecera y finaliza en una concentración de 9.40 mg/l.

El OD es la cantidad de oxígeno que está disuelta en el agua y que es esencial para

los sistemas hídricos. En ese sentido se convierte en un indicador de los sistemas
hídricos, en cuanto a los niveles de contaminación y puede dar soporte a las
condiciones bióticas de este tipo de ecosistemas. Generalmente, un nivel más alto de
oxígeno disuelto indica agua de mejor calidad., en caso contrario, si los niveles de
oxígeno disuelto son demasiado bajos, algunos peces y otros organismos no pueden
sobrevivir. Gran parte del oxígeno disuelto en el agua proviene del oxígeno en el aire
que se ha disuelto en el agua. Parte del oxígeno disuelto en el agua es el resultado
de la fotosíntesis de las plantas acuáticas. En el caso del OD, se define en 5.40 mg/l
y termina aguas abajo con un valor de 5.29 mg/l, con una concentración de
saturación alrededor de 6.50mg/l. La información de OD, permite establecer que el
déficit de oxígeno es bajo y obedece a un promedio de 1.075 mg/l, condición muy
favorable para el comportamiento del sistema, en lo relacionado a su capacidad de
autodepuración y dilución.

De manera general se puede establecer según criterio que niveles de oxígeno entre
0 y 4.00 mg/l, identifican una calidad del agua mala (algunas poblaciones de peces y
macro invertebrados empezarán a bajar); valores entre 4.1 a 7.9 mg/l, identifican
calidad aceptable; valores entre 8.00 y 12 mg/l la calidad del agua es buena.

En lo relacionado a las coliformes totales, se pueden observar valores en la cabecera

de 680 NMP/ 100ml, disminuyendo en el trayecto hasta un valor final de 490
NMP/100ml. Las coliformes son una familia de bacterias que se encuentran
comúnmente en las plantas, suelo y animales, incluyendo los humanos. Su presencia
en los cuerpos hídricos indica presencia de aguas residuales.

La asimilación es la capacidad del agua natural de recibir aguas residuales o

materiales tóxicos sin que tengan efectos negativos y sin daño para la vida acuática
o también puede entenderse como la capacidad del agua de purificarse de agentes
contaminantes. En otro contexto, pero de manera integrada, la dilución, relacionada
con disminución de la concentración de una disolución añadiendo disolvente, se
puede llegar a conocer un sistema frente a su comportamiento.

Una fuente hídrica se puede evaluar a través de su capacidad para asimilar y diluir
los contaminantes que recibe y por ello las principales transformaciones que se
presentan durante el proceso de autodepuración son de tipo bioquímico, además los
compuestos orgánicos más variados experimentan la acción de enzimas bacterianas
y de otros organismos que los transforman en células vivas o son oxidadas vía la
respiración.

La capacidad de asimilación y de dilución de un contaminante o también conocida

como la autopurificación natural del sistema hídrico, está gobernada por varios
factores, entre los cuales están: La cantidad del agua que lleva la corriente, la calidad
del agua de la corriente, la cantidad y calidad de las descargas de agua residuales,
otras fuentes contaminantes, la presencia de sustancias tóxicas a los
microorganismos y el régimen de flujo de la corriente. En ese conjunto, el modelo
integra estos aspectos y presenta una condición muy favorable para el Rio Iza en lo
que respecta a asimilación, por la condición de oxígeno disuelto y los resultados
encontrados, además las condiciones de reaireación definen una dinámica ajustada
por la misma morfometría y los aspectos hidráulicos asociados. La dilución se
acentúa en la medida en que la fuente en su cabecera cuenta con 0.42 m 3/s, que en
su viaje hacia el punto aguas abajo y luego de recibir los vertimientos permite varias
favorablemente las concentraciones de materia orgánica en términos de DBOu, sin
embargo, también se puede considerar que valores de DBO > 5 mg/l, identifican
condiciones de contaminación en la fuente.

Posteriormente se ha querido trabajar con algunos escenarios que permitan

establecer la dinámica del sistema, frente a las posibles modificaciones. En ese
orden de ideas se hace un modelado con el caudal medido sin los vertimientos, luego
con la PTAR en relación a las remociones presentadas para el diseño. En el análisis
realizado, se puede establecer que, aunque existe variación y por su puesto
modificación en los valores arrojados, se sigue manteniendo la condición favorable
de asimilación, en lo que respecta a la dinámica de reaireación, remoción y
oxigenación, aspectos que son gobernados a partir de las constantes cinéticas y por
otro lado la dilución es satisfactoria, con el caudal medido trabajado de 0.42 m 3/s.
Luego en el escenario con la PTAR en donde se definen las respectivas remociones
para el vertimiento, es evidente el mejoramiento de las condiciones en el rio, debido
a que el OD se mantiene en concentraciones bastante favorables para la fuente,
además que los déficits de oxígeno son bajos, lo que hace que las concentraciones
de saturación se mantengan y no varíe significativamente.

Los otros escenarios propuestos plantean la misma situación anterior, sistema sin
vertimientos y sistema con tratamiento, pero en una condición de flujo mínimo,
definido como el 20% del caudal medido, es decir con un flujo de 0.336 m 3/s. Esta
situación se puede tomar como condición extrema. La modelación permite establecer
la misma favorabilidad de la situación anterior, debido a las condiciones hidráulicas,
morfométricas y lo que implica los diferentes balances para garantizar las óptimas
condiciones de calidad necesarias, para que la fuente cumpla con sus funciones de
autopurificación, remoción y reaireación, para que el sistema hídrico no se altere y
por el contrario se mantenga en sus procesos de transformación.
 1.20 EVALUACIÓN OBJETIVOS DE CALIDAD DE LA FUENTE DENOMINADA
RIO TOTA, RECEPTOR DE LOS VERTIIENTOS MUNICIPALES, UBICADOS EN
LA JURISDICCIÓN DEL MUNICIPIO DE IZA.

Mediante resolución 0830 del 04 de septiembre de 2008, CORPOBOYACA

estableció los objetivos de calidad para la fuente denominada Rio Tota, receptor de
los vertimientos municipales, ubicados en la jurisdicción del municipio de Iza, por lo
cual se define el sector comprendido desde la acequia la Compañía con el Rio Tota
hasta la confluencia del Rio Pesca con el Rio Tota. En este caso los usos del recurso
están definidos para agrícola restringido, pecuario y consumo humano.

El siguiente cuadro presenta la información relacionada con los indicadores y

parámetros de calidad, en el contexto de los objetivos de calidad preestablecidos.

Tabla 39. Indicadores y parámetros de calidad para el sector del Rio Tota
INDICADOR Y PARAMETROS DE CALIDAD
(NMP/100ml)COLIFORMES FECALES

(mg/l N –NH3)AMONIACAL

(mg/l N-NO3)NITRATOS

(mg/l N-NO2)NITRITOS
USOS DEL RECURSO

COLIFORMES TOTALES

(mg/l)DBO

OD (mg/l)
SECTOR
TRAMO

(NMP/100ml)

pH
OBJETIVOS DE CALIDAD
1 Sector Rio Tota.(Desde DyA
la confluencia de la Agrícola
acequia la compañía restringido,
100-1000 - 1-3 4 0.5 10 0.1 6.5-8.5
con el rio tota hasta la pecuario,
confluencia del rio consumo
pesca con el rio tota) humano
Fuente. Resolución 0830 del 04 de septiembre de 2008. Objetivos de Calidad
CORPOBOYACA

En la tabla siguiente se presenta los criterios, normas o estándares de calidad del

cuerpo receptor

Tabla 40. Criterios, normas o estándares de calidad del cuerpo receptor

PARÁMETROS USOS DEL RECURSO: Resolución 0830 del RIO IZA CUMPLE RESULTADOS CUMPLE
 04 de septiembre de 2008

RESTRINGIDO
AGRÍCOLA

PECUARIO

CONSUMO
HUMANO
MODELACIÓN
CABECERA
SI NO CON PTAR SI NO

Coliformes
Totales 100 - 1000 100 - 1000 100 - 1000 680 X 269.4 X
(NMP/100ml)
Coliformes
Fecales - − - 450 X - X
(NMP/100ml)
DBO (mg/l) 1-3 1-3 1-3 <2 X 3.57 X
OD (mg/l) 4 4 4 5.40 X 5.32 X
Nitrógeno
Amoniacal 0.5 0.5 0.5 <0.42 X 0.33 X
(mg NH3 – N/ L)
Nitratos – N
10 10 10 <0.38 X 0.36 X
(mg NO3 –N/L)
Nitritos - N
(mg NO2 –N/L)
0.1 0.1 0.1 <0.018 X 0.36 X
pH 6.5 – 8.5 6.5 – 8.5 6.5 – 8.5 7.37 X 7.16 X
Fuente. INCO A&J S.A.S

Al analizar la información de la tabla anterior se puede establecer que la fuente Rio

Iza en su condición aguas arriba cumple con la mayoría de parámetros para los usos
definido, sin embargo se presenta un inconveniente con el tema de las coliformes
fecales. Se aclara que esta verificación se hace con los datos reportados para la
cabecera o lo que se estableció como aguas arriba del vertimiento colegio.

Para el caso de los resultados de la modelación con el caudal medido e

implementando el tratamiento, se establece que no cumple en el caso de los
coliformes fecales y la demanda bioquímica de oxígeno (DBO), sin embargo los
otros parámetros cumplen, en esta condición. Los datos obtenidos de la simulación
permiten determinar que la fuente receptora Río Iza tiene alta capacidad de dilución y
de autodepuración.
 1.21 REFERENCIAS BIBLIOGRÁFICAS

Chapra, S. C. (1997). Surface Water-Quality Modelling, the McGraw-Hill Companies. Inc.,

New York.

Chapra Steve and Pelletier Greg. (2008). QUAL2Kw theory and documentation (version
5.1). A modeling framework for simulating river and stream water quality.

Colombia, Ministerio de Ambiente, Vivienda y Desarrollo Territorial. (2010). Decreto 3930

de 2010

Colombia. Ministerio de Ambiente, Vivienda y Desarrollo Territorial. (2010). Política

Nacional para la Gestión Integral del Recurso Hídrico. Bogotá, D.C.: Colombia,
Garcia T., Pelletier G., Diaz J. (2009). Water quality simulation of the Chicamocha River,
Colombia. An application of the QUAL2Kw Model. Ecological Modeling.

Instituto de Hidrología, Meteorología y Estudios Ambientales – IDEAM. (2004). Guía Para

El Monitoreo y Seguimiento Del Agua. Bogotá.

Instituto de Hidrología, Meteorología y Estudios Ambientales – IDEAM. (2010). Estudio

Nacional del Agua 2010.

Jiménez, A. A. (2000). Determinación de los parámetros físico-químicos de calidad de las

aguas. Revista interdisciplinar de gestión ambiental, 2(23), 12-19.

Pelletier, G. J., Chapra, S. C. (2003). QUAL2K: A Modeling Framework for Simulating

River and Stream Water Quality. Documentation and User’s Manual. Civil and
Environmental Engineering Dept., Tufts University, Medford, MA. 2003.

Pelletier, G.j. Chapra, S. (2004). QUAL2Kw theory and documentation (version 5.1). A
modeling framework for emulating river and stream water quality. Civil and Environmental
Engineering Dept., Tufts University, Medford, MA. 104 pp. 2004.

Pelletier, G. J., Chapra, S., TAO, H. (2005). QUAL2Kw - A framework for modeling water
quality in streams and rivers using a genetic algorithm for calibration.

Vera Puerto, I. L. (2007). Aplicación de técnica de optimización mediante algoritmos

genéticos para calibración de modelo Qual2k como una aproximación a la modelación de
la calidad del agua de los principales ríos de la zona urbana de Bogotá DC.