Evaluación de La Calidad de Las Aguas Del Río Katari, La Paz, Bolivia, Mediante Un Modelo Matemático
Evaluación de La Calidad de Las Aguas Del Río Katari, La Paz, Bolivia, Mediante Un Modelo Matemático
Evaluación de La Calidad de Las Aguas Del Río Katari, La Paz, Bolivia, Mediante Un Modelo Matemático
FACULTAD DE INGENIERIA
CARRERA DE INGENIERIA AMBIENTAL
IIDEPROQ
INSTITUTO DE INVESTIGACIN Y DESARROLLO DE
PROCESOS QUIMICOS
PROYECTO DE GRADO
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ,
BOLIVIA, MEDIANTE UN MODELO MATEMTICO
POSTULANTE:
Michael Russel Osina Torrez
ASESOR:
Ing. Waldo P. Vargas Ballester
La Paz Bolivia
2011
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
DEDICATORIA
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
AGRADECIMIENTOS
El presente proyecto fue realizado con mucho esfuerzo, dedicacin y amor gracias a todas las
personas e instituciones que directa e indirectamente hicieron posible su culminacin e
indudablemente gracias a Dios.
Quisiera expresar mis ms sinceros y profundos agradecimientos a mi asesor M.sc Waldo P.
Vargas Ballester por el asesoramiento del presente proyecto de grado, sus consejos tcnicos,
constante apoyo, enseanzas y su amistad.
A toda mi familia, que confiaron en m en todo momento, especialmente a mi madre y
abuelita por su paciencia, comprensin y apoyo en todo momento a lo largo de mi vida.
A Paola por estar siempre a mi lado en los buenos y malos momentos, por su apoyo y
comprensin.
A mi primo Emir, por proveer tan generosamente de su tiempo, para las actividades de
trabajo de campo.
Al Instituto de Boliviano de Ciencia y Tecnologa Nuclear (IBTEN) por su apoyo en los
anlisis e informacin proporcionada.
Al Ing. M. Sc. Luis Romero Bolaos, por su apoyo desinteresado para el trabajo de campo y
por su amistad, palabras de nimo y sabios consejos.
Al Ing. Jorge Chungara del Instituto Boliviano de Ciencia y Tecnologa Nuclear (IBTEN),
por el apoyo para el trabajo de campo, informacin y por su amistad.
Al Tc. Sup. Wilfredo Tarquino, Tc. Rocio Choque, Lic. Samuel Fernndez, Tc. Alex
Chiara y Sr. Eliseo Rojas P. (IBTEN), por el apoyo para el trabajo de campo y por su
amistad.
Al Ing. Luis Noriega y Tc. Franklin Mamani del Servicio Nacional de Meteorologa e
Hidrologa (SENAMHI), por el apoyo decidido de su institucin para el trabajo de campo y
por la informacin proporcionada.
Al Instituto de Investigacin y Desarrollo de Procesos Qumicos (IIDEPROQ), por el apoyo
en el anlisis de muestras.
A todos los miembros investigadores del IIDEPROQ, por su amistad y su cario
desinteresado.
Al Ing. Edwin Astorga del Instituto de Ingeniera Sanitaria (IIS), por el apoyo en el anlisis
de muestras.
A Lic. Joaqun Gutirrez Daz del Centro de Informacin, Gestin y Educacin Ambiental
(CIGEA) de Cuba, por sus consejos tcnicos, enseanzas y amistad.
Michael Russel Osina Torrez
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
RESUMEN
Actualmente la contaminacin hdrica es uno de los problemas ambientales ms grandes a
nivel mundial, ya que la escasez del agua dulce y la creciente contaminacin de sta, estn
haciendo que su uso sea cada vez ms dificultoso. Bolivia no es ajena a esta problemtica y
los casos de contaminacin en los cuerpos naturales de agua en su territorio, son cada vez
ms significativos y frecuentes.
El establecimiento de objetivos de calidad ha de realizarse teniendo en cuenta tanto la
cantidad de contaminante vertido como su comportamiento posterior en el medio en funcin
del tiempo y cmo puede interferir con el resto de sustancias presentes en el medio. En este
sentido, se hace necesario utilizar modelos matemticos de simulacin conocidos como
"modelos de calidad" o "modelos de transporte de contaminantes que permitan predecir la
evolucin a largo plazo del sistema en su conjunto ante una accin determinada. La validez
de la prediccin vendr determinada por la capacidad del modelo propuesto de reflejar el
comportamiento del sistema. Por ello, es necesario, por una parte, conocer todos los
fenmenos que afectan al comportamiento de las distintas sustancias presentes en el sistema y
por otra conocer muy bien las caractersticas del sistema, de tal forma desde el punto de vista
del enfoque sistmico realizar un seccionamiento ambiental representativo.
En este trabajo se presenta el estudio de la calidad del agua del Ro Katari el cual comprende
la confluencia del Ro Pallina y Ro Colorado hasta la descarga a la Baha de Cohana
utilizando modelos matemticos. En la segmentacin ambiental se considero cuatro tramos
divididos por cinco puntos de muestreo a lo largo de la corriente. Se utilizaron como datos de
entrada a los modelos: demanda biolgica de oxigeno (DBO 5 ), oxigeno disuelto (OD),
temperatura (T), altitud sobre nivel del mar, velocidad del ro (V), caudal (Q), profundidad
(H), ancho del ro (W), rea transversal (A), constante de desoxigenacin (K1 ), constante de
reaireacin (K 2 ), longitud de cada tramo (L). Tambin se considero en el anlisis de calidad
de agua parmetros como el pH, conductividad elctrica (CE), Concentracin de saturacin
de oxigeno disuelto (Cs). Previamente se realizo varias precampaas para conocer las
caractersticas del lugar de estudio, para posteriormente realizar la segmentacin ambiental
requerida por los modelos y georeferenciar los puntos de muestreo para luego poder
utilizarlos como referencia en futuras campaas de muestreo. Se realizaron tres campaas de
muestreo en poca seca el 2010 y una en poca de lluvia el 2011, tambin se utilizo datos de
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
ABSTRACT
Nowadays, water pollution is one of the largest environmental problems worldwide, the
scarcity of fresh water and increased pollution of it are making its use more difficult than
used to be. Bolivia is not apart to this problem and, cases of contamination of natural water
bodies in its territory are becoming more significant and frequent.
The establishment of quality objectives must be carried out taking into account both the
amount of pollutants and their subsequent behavior in the environment in terms of time, and
how it can interfere with other substances in the environment. In this sense, it is necessary to
use mathematical simulation models known as "quality models" or "contaminant transport
models" that predict the long-term system as a whole to a particular action. The validity of
the prediction will be determined by the ability of the proposed model to reflect the system
behavior. It is necessary; know all the phenomena that affect the behavior of the substances
presented in the system and other well known system features, so from the systemic approach
point of view we carry out a representative environmental segmentation.
This thesis presents a study of water quality in the Katari River which includes the
confluence of the Colorado River and Pallina River until discharge to the Cohana Bay using
mathematical models. For the environmental segmentation was considered four sections
divided by five sampling points along the stream. Were used as input data to models:
biological oxygen demand (BOD 5 ), dissolved oxygen (DO), temperature (T), altitude above
sea level, river velocity (V), flow (Q), depth (H ), river width (W), cross-sectional area (A),
deoxygenation constant (K1 ), reaeration constant (K2 ), length of each section (L). Also was
considered in the water quality analysis parameters such as pH, electrical conductivity (EC),
dissolved oxygen concentration saturation (Cs). Previously we visited the place to determine
the characteristics of the study area to later perform the required environmental segmentation
required by the models and georeferenced sampling points to be can use them as reference for
future sampling collections. We carried out three sample collections in dry season in the 2010
and one in rainy season in 2011, also; we use data from previous sample collections carried
out by the National Institute of Nuclear Science and Technology (IBTEN) and National
Service of Meteorology and Hydrology (SENAMHI).
We used three computational models which are RIOSep, SIMOD and WASP 7.4. The
information obtained was prepared to introduce to the models and then proceed to the
Michael Russel Osina Torrez
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
calibration. The models were calibrated with 5 sample collections and validated with 4
sample collections. RIOSep and WASP 7.4 model does not have the option of calculating the
constant of deoxygenation, which was determined from the equation of Chapra (1997) and
we get a correlation from the data obtained by the model SIMOD between the constant of
reaeration (K2 ) and deoxygenation (K1 ). Once the models were calibrated and validated, the
evaluation of the models was done and we verifies that the WASP 7.4 does not calibrated
because it requires a lot of input information which is not available. However, the model
SIMOD and RIOSep were well calibrated. It was concluded that the model SIMOD is the
most optimal for the study of water quality because it had the lowest percentage of error,
correlation coefficient close to 1 and a mean difference close to zero showing better
correlation of the analyzed values than RIOSep model, plus its easy handling, low input data
and calculation of the constant deoxygenation (K 1 ), which saves time and money. Finally,
scenario analysis performed with the model that presented better results in terms of
calibration, validation and sensitivity analysis being the model SIMOD. As a result of
prediction scenarios we have:
-
All additions and improvements to be made in the WWTP Puchukollo is not enough
to solve the pollution in Cohana Bay; unless all the wastewater generated in El Alto,
Viacha and Laja is collected and treated.
If there is un increased discharge of pollutants from the city of El Alto, Laja, Viacha,
this would cause more pollution to the Katari River because the river could not debug
it naturally.
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
TABLA DE CONTENIDO
TABLAS ........................................................................................................................ VII
GRFICO S.................................................................................................................... IX
CAPTULO 1 INTRODUCCIN .................................................................................. 1
1.1. JUSTIFICACIN ........................................................................................................ 2
1.2. DEFINICIN DE OBJETIVOS .................................................................................. 5
1.2.1. OBJETIVO GENERAL ............................................................................................................ 5
1.2.2.
Pgina | I
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
3.4.3.2.2.
Pgina | II
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
4.2.
4.3.
LIMITACIONES.................................................................................................. 74
4.3.1. SIMOD 1.01 ................................................................................................... 74
4.3.2. RIOSep ........................................................................................................... 74
4.3.3. WASP 7.4 ....................................................................................................... 74
4.4.
Pgina | III
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
4.5.
4.6.
4.7.
4.8.
4.9.
Pgina | IV
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
7.1.1.1.1.
7.1.1.2.
7.1.1.2.1.
7.1.1.2.2.
7.1.1.2.3.
7.1.1.3.
7.1.1.3.1.
7.1.1.3.2.
7.1.1.3.3.
7.2.
7.1.2.1.
7.1.2.2.
7.3.
Pgina | V
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
9.2.
9.3.
Pgina | VI
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
TABLAS
Tabla 1.1. Valores mximos admisibles de parmetros en cuerpos receptores....................... 4
Tabla 2.1. Calendario Agrcola, Cuenca del Rio Katari, La Paz, Bolivia ............................. 14
Tabla 3.1. Clasificacin de los cuerpos de agua segn su aptitud de uso ............................. 17
Tabla 3.2. Formulas para el clculo de la constante de reaireacin ...................................... 67
Tabla 4.1. Coeficientes de reacciones corregidas por la temperatura ................................... 87
Tabla 5.1. Estaciones de monitoreo: ubicacin y tipo de muestra tomada (W = Agua, IS = in
situ = Laboratorio, S = Simple) .......................................................................................... 99
Tabla 5.2. Parmetros fsicos y qumicos medidos en laboratorio, equipos utilizados,
periodicidad de la medicin y protocolo de medida seguido ................................................ 99
Tabla 5.3. Detalle de los parmetros medidos in situ y en laboratorio en las campaas
realizadas .......................................................................................................................... 101
Tabla 6.1. Variacin de la concentracin de saturacin en funcin de la CE ..................... 106
Tabla 7.1. Combinacin de corridas con el modelo SIMOD v1.01.................................... 117
Tabla 7.2. Rango de valores para las constantes de reaireacin y desoxigenacin ............. 117
Tabla 7.3. Rango de valores para las constantes de reaireacin y desoxigenacin ............. 118
Tabla 7.4. Resultados Calibrado con el modelo SIMOD v1.01 (3ra campaa) .................. 119
Tabla 7.5. Resultados Calibrado con el modelo SIMOD v1.01 (5ta campaa) .................. 120
Tabla 7.6. Combinacin de corridas con el modelo RIOSep ............................................. 122
Tabla 7.7. Constantes de reaireacin para cada tramo de cada campaa por cinco diferentes
mtodos ............................................................................................................................ 123
Tabla 7.8. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la
ecuacin ajustada .............................................................................................................. 124
Tabla 7.9. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la
ecuacin Chapra................................................................................................................ 124
Tabla 7.10. Resultados Calibrado con el modelo RIOSep (3ra campaa) .......................... 125
Tabla 7.11. Resultados Calibrado con el modelo RIOSep (5ta campaa) .......................... 127
Tabla 7.12. Combinacin de corridas con el modelo WASP 7.4 ....................................... 130
Tabla 7.13. Constantes de reaireacin para cada tramo de cada campaa por cinco diferentes
mtodos ............................................................................................................................ 131
Tabla 7.14. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la
ecuacin ajustada ............................................................................................................. 132
Tabla 7.15. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la
ecuacin Chapra ............................................................................................................... 132
Tabla 7.16. Resultados Calibrado con el modelo WASP 7.4 (1ra campaa) ...................... 133
Tabla 7.17. Resultados Validacin del modelo SIMOD (1ra campaa) ............................ 135
Tabla 7.18. Resultados Validacin del modelo SIMOD (8va campaa) ........................... 136
Tabla 7.19. Resultados Validacin del modelo RIOSep (1ra campaa) ............................ 137
Tabla 7.20. Resultados Validacin del modelo RIOSep (8va campaa) ........................... 139
Tabla 7.21. Coeficientes de correlacin para las corridas mejor calibradas de DBO 5 ........ 141
Tabla 7.22. Coeficientes de correlacin para las corridas mejor calibradas de OD ............ 142
Tabla 7.23. Errores e incertidumbre para las corridas mejor calibradas de DBO 5 , Modelo
SIMOD ............................................................................................................................. 143
Tabla 7.24. Errores e incertidumbre para las corridas mejor calibradas de DBO 5 , Modelo
RIOSep ............................................................................................................................. 144
Tabla 7.25. Errores e incertidumbre para las corridas mejor calibradas de OD, Modelo
SIMOD ............................................................................................................................. 144
Pgina | VII
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tabla 7.26. Errores e incertidumbre para las corridas mejor calibradas de OD, Modelo
RIOSep ............................................................................................................................ 145
Tabla 7.27. Diferencia promedio para las corridas mejor calibradas de DBO 5 .................. 146
Tabla 7.28. Diferencia promedio para las corridas mejor calibradas de OD ..................... 147
Tabla 7.29. Caudales mximos estimados de las descargas de la ciudad de El Alto, Viacha y
Laja a partir de datos tomados por el IBTEN ..................................................................... 149
Tabla 7.30. Datos de proyeccin de la carga que llegara a la PTAR Puchukollo hasta el 2035
......................................................................................................................................... 153
Tabla 7.31. Datos Promedios de las cargas que llegaran a la PTAR Puchukollo................ 153
Tabla 7.32. Ampliacin de la PTAR Puchukollo .............................................................. 155
Tabla 7.33. Potencial de la eliminacin de la DBO ........................................................... 156
Tabla 7.34. Eficiencia de tratamientos .............................................................................. 156
Tabla 7.35. Caudales estimados, concentraciones de OD y DBO5 de descarga de las
principales fuentes de contaminacin de la cuenca del Ro Katari ..................................... 157
Tabla 7.36. Datos de ingreso al modelo para prediccin del primer escenario (descarga del
Ro Pallina) ....................................................................................................................... 158
Tabla 7.37. Resultados Escenario de Prediccin 1 con el modelo SIMOD v1.01 .............. 158
Tabla 7.38. Datos de entrada al modelo segn la aplicacin de la eficiencia de los distintos
tratamientos ..................................................................................................................... 159
Tabla 7.39. Resultados Escenario de Prediccin 2 (Aplicando Tratamiento Primario) con el
modelo SIMOD v1.01 ....................................................................................................... 159
Tabla 7.40. Resultados Escenario de Prediccin 2 (Aplicando Tratamiento Secundario) con
el modelo SIMOD v1.01 ................................................................................................... 160
Tabla 7.41. Resultados Escenario de Prediccin 2 (Aplicando Tratamiento Terciario) con el
modelo SIMOD v1.01 ....................................................................................................... 161
Tabla 7.42. Datos de entrada al modelo sin la aplicacin de tratamientos a los principales
efluentes de la cuenca del Ro Katari................................................................................. 162
Tabla 7.43. Resultados Escenario de Prediccin 3 con el modelo SIMOD v1.01 .............. 162
Tabla 9.1. Rangos de las constantes de reaireacin y desoxigenacin ............................... 168
Pgina | VIII
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
GRFICOS
Grfico 1.1. Resultados de ndices de calidad de agua (ICA) de la cuenca del Ro Katari ... 3
Grfico 2.1. Precipitaciones en la Estacin Tambillo, Rio Katari ........................................ 9
Grfico 2.2. Puntos de muestreo en el Rio Katari.............................................................. 15
Grfico 3.1. Variacin en el tiempo de la DBO: a) Remanente, b) Satisfecha ................... 21
Grfico 3.2. Relacin de la DBO Carboncea (DBOC) y la DBO Nitrogencea (DBON) . 22
Grfico 3.3. Concepto de los modelos de calidad de agua ................................................. 29
Grfico 3.4. Concepto de transporte advectivo.................................................................. 32
Grfico 3.5. Masa de un constituyente que entra y sale de un volumen de control debido a
la adveccin en direccin x ................................................................................................. 33
Grfico 3.6. Representacin esquemtica de la transferencia de O2 en cuatro pasos desde la
fase vapor a la fase solucin ................................................................................................ 35
Grfico 3.7. Relaciones entre la presin parcial del O2 en la fase gaseosa y la concentracin
en la fase disolucin ............................................................................................................ 37
Grfico 3.8. Esquema del modelo de la pelcula ............................................................... 40
Grfico 3.9. Esquema del modelo de la penetracin.......................................................... 40
Grfico 3.10. Fases de Crecimiento de una Poblacin Bacteriana .................................... 43
Grfico 3.11. Tasa de Crecimiento Especfico versus Concentracin de Sustrato .............. 45
Grfico 3.12. Capacidad de Crecimiento para una serie de Cuatro Experiencias ............... 46
Grfico 3.13. Transmisin longitudinal de contaminacin ................................................ 51
Grfico 3.14. Elemento del cuerpo de agua....................................................................... 52
Grfico 3.15. Esquema de un sistema ............................................................................... 55
Grfico 3.16. Esquema de un ro con compuestos conservativos ....................................... 56
Grfico 3.17. Volumen de control de un ro ...................................................................... 57
Grfico 3.18. Curva caracterstica del agotamiento de OD obtenida mediante la aplicacin
de Streeter and Phelps ......................................................................................................... 59
Grfico 3.19. Comparacin del comportamiento de la DBO en el rio y laboratorio a)
decaimiento tpico de DBO en el rio y b) determinacin de la DBO en el laboratorio ......... 60
Grfico 3.20. Materia orgnica vs. Tiempo ....................................................................... 61
Grfico 3.21. Relacin entre la DBO y la materia orgnica ............................................... 62
Grfico 3.22. Esquema de la configuracin de un modelo, como sistema simplificado ..... 69
Grfico 4.1. Esquema del Seccionamiento Ambiental para el modelo SIMOD ................. 71
Grfico 4.2. Esquema del Seccionamiento Ambiental para el modelo RIOSep ................. 72
Grfico 4.3. Esquema del Seccionamiento Ambiental para el modelo WASP 7.4 ............. 72
Grfico 4.4. Sistema de coordenadas para la ecuacin de balance de masa ....................... 79
Grfico 5.1. Maniobra tpica de toma de muestra en dos estaciones (PC1 y PC5) situada en
el medio del ro y sobre un puente del Ro Colorado y Ro Katari .................................... 102
Grfico 6.1. Grafico del Oxigeno Disuelto (OD) de los diferentes puntos en las distintas
campaas realizadas .......................................................................................................... 103
Grfico 6.2. Grafico de la Salinidad (S) de los diferentes puntos en las distintas campaas
realizadas .......................................................................................................................... 104
Grfico 6.3. Grafico de la Concentracin de Saturacin del OD (Cs) de los diferentes
puntos en las distintas campaas realizadas (Ec. 3-12 y 3-14)............................................ 105
Grfico 6.4. Grafico de la Concentracin de Saturacin del OD (Cs) de los diferentes
puntos en las distintas campaas realizadas (Ec. 4-56 y 4-57)............................................ 105
Grfico 6.5. Grafico de la Demanda Bioqumica de Oxigeno (DBO) de los diferentes
puntos en las distintas campaas realizadas ....................................................................... 107
Pgina | IX
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | X
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | XI
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
1. INTRODUCCIN
Actualmente la contaminacin hdrica es uno de los problemas ambientales ms grandes a
nivel mundial, ya que la escasez del agua dulce y la creciente contaminacin de sta, estn
haciendo que su uso sea cada vez ms dificultoso. Bolivia no es ajena a esta problemtica y
los casos de contaminacin en los cuerpos naturales de agua en su territorio, son cada vez
ms significativos y frecuentes.
En la actualidad est rebrotando en el pas, una importante corriente de conservacin
ambiental que promueve en el caso especfico del Departamento de La Paz, la necesidad de
implantar medidas de mitigacin que puedan frenar y hasta revertir los niveles de
contaminacin que afectan las aguas de distintos cuerpos naturales, como por ejemplo el Lago
Titicaca, entre las cuales se destaca la necesidad perentoria de tratar las aguas residuales que
se descargan a dicho cuerpo de agua, ya que el Lago Titicaca es uno de los lagos navegables
de altura con mayor extensin en el mundo y tambin es un importante destino turstico, que
en la ltima dcada se ha visto significativamente afectado por la creciente contaminacin
hdrica, la degradacin del paisaje, prdida del hbitat, y la aparicin de procesos eutrficos,
localizados en numerosos sitios.
La Baha de Cohana es considerada como el sitio ms contaminado de la parte boliviana del
lago Titicaca, y el segundo ms contaminado de todo el lago despus de la Baha de Puno en
Per. Las problemticas ambientales, y en especial las que se refieren a la contaminacin
hdrica, estn relacionadas con variables causales internas y externas. La carencia de una
gestin integral de residuos slidos y la prctica de la ganadera intensiva en las poblaciones
cercanas a los ros que desembocan en el Lago Titicaca, crticamente sustentada por las
macrfitas acuticas (especialmente Myriophyllum qutense conocida comnmente como
Llacho, Chchinqui y Schoenoplectus californicus ssp. conocida comnmente como totora),
representan las variables internas, mientras que la variedad de contaminantes y materiales
residuales arrastrados por los ros que desembocan en la Baha, desde la ciudad de El Alto y
diversas poblaciones, constituyen las variables externas.
Esta situacin afecta significativamente el nivel de la calidad ambiental de los cuerpos de
agua de todo el sistema TiticacaDesaguaderoSalar de PoopSalar de Coipasa (TDPS), el
uso y consumo de los mismos, y tambin en todo lo concerniente a la salud pblica y la
epidemiologa. Las actividades productivas del lugar y, principalmente, el cmo se llevan a
Michael Russel Osina Torrez
Pgina | 1
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 2
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
por lo que los nutrientes no se mantienen en el fondo, sino que estn disponibles para el
aprovechamiento de los productores primarios (algas y macrfitas). Su presencia causa un
desarrollo masivo de stas, especialmente la Lemna (o lenteja de agua) y otras plantas no
enraizadas que viven en la superficie del agua.
Debido al exceso de nutrientes, las plantas crecen en abundancia y colonizados por otros
organismos, propios de aguas contaminadas. Ms tarde, cuando mueren, se pudren y llenan el
agua de malos olores y le dan un aspecto nauseabundo, disminuyendo drsticamente su
calidad. El proceso de degradacin de la materia orgnica consume una gran cantidad del
oxgeno disuelto y las aguas dejan de ser aptas para la mayor parte de los seres vivos.
Pgina | 3
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
PARAMETRO
UNIDAD
pH
Temperatura
Slidos disueltos
CLASE A
6.0 a 8.5
3 C de c.
receptor
CLASE B
CLASE C
CLASE D
6.0 a 9.0
6.0 a 9.0
6.0 a 9.0
3 C de c. receptor
3 C de c. receptor
3 C de c. receptor
mg/l
1000
1000
1500
1500
Aceites y Grasas
mg/l
Ausentes
Ausentes
0.3
DBO 5
mg/l
<2
<5
<20
<30
DQO
mg/l
<5
<10
<40
<60
NMP Colifecales
N/100ml
<1000 y <200 en
<5000 y <1000 en
<50000 y <5000 en
de muestras
80% de muestras
80% de muestras
80% de muestras
Parsitos
N/l
<1
<1
<1
<1
Color mg Pt/l
mg/l
<10
<50
<100
<200
Oxigeno disuelto
mg/l
>70% sat.
>60% sat.
50% sat.
Turbidez
UNT
<10
<50
<100 - <2000***
<200 10000***
Slidos
mg/l -
Sedimentables
ml/l
<10 mg/l
Aluminio
mg/l
0.2 c. Al
0.5 c. Al
1.0 c. Al
1.0 c. Al
Amoniaco
mg/l
0.05 c. NH
1.0 c. NH
2 c. NH
4 c. NH
Antimonio
mg/l
0.01 c. Sb
0.01 c. Sb
0.01 c. Sb
0.01 c. Sb
Arsnico total
mg/l
0.05 As
0.05 c. As
0.05 c. As
0.1 c. As
Benceno
g/l
2.0 c. Be
6.0 c: Be
10.0 c. Be
10.0
Bario
mg/l
1.0 0.05 c. Ba
1.0 c. Ba
2.0 c. Ba
5.0 c. Ba
Berilio
mg/l
0.001 c. Be
0.001 c. Be
0.001 c. Be
0.001 c. Be
Boro
mg/l
1.0 c. B
1.0 c. B
1.0 c. B
1.0 c. B
Calcio
mg/l
200
300
300
400
Cadmio
mg/l
0.005
0.005
0.005
0.005
Cianuros
mg/l
0.002
0.1
0.2
0.2
Cloruros
mg/l
250 c. Cl
300 c. Cl
400 c. Cl
500 c. Cl
Cobre
mg/l
0.05 c. Cu
1.0 c. Cu
1.0 c. Cu
1.0 c. Cu
Cobalto
mg/l
0.1 c. Co
0.2 c. Co
totales
0.05 c. Cr total
0.2 c. Co
+6
0.2 c. Co
+6
0.05 c. Cr
0.05 c. Cr
0.05 c. Cr+6
0.6c. Cr+3
0.6c. Cr+3
1.1 c. Cr+3
Cromo Trivalente
mg/l
Hierro soluble
mg/l
0.3 c. Fe
0.3 c. Fe
0.1 c. Fe
1.0 c. Fe
Manganeso
mg/l
0.5 c. Mn
1.0 c. Mn
1.0 c. Mn
1.0 c. Mn
Mercurio
mg/l
0.001 Hg
0.001 Hg
0.001 Hg
0.001 Hg
Nquel
mg/l
0.05 c. Ni
0.05 c. Ni
0.5 c. Ni
0.5 c. Ni
Nitrato
mg/l
20.0 c. NO 3
50.0 c. NO 3
50.0 c. NO 3
50.0 c. NO 3
Nitrito
mg/l
<1.0 c. N
<1.0 c. N
<1.0 c. N
<1.0 c. N
Nitrgeno Total
mg/l
5 c. N
12 c. N
12 c. N
12 c. N
Plomo
mg/l
0.05 c. Pb
0.05 c. Pb
0.05 c. Pb
0.1 c. Pb
Pgina | 4
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Desde hace algn tiempo, las comunidades pesqueras de esta regin han venido denunciando
la muerte masiva de los peces del sector de la Baha de Cohana. Por otra parte, la prdida de
la calidad del pescado ha generado un rechazo del consumidor a este producto, especialmente
a las especies pisccolas nativas. Estas dos situaciones han creado un problema social entre las
comunidades pesqueras, que ya no pueden contar con el ingreso de la pesca.
En este sentido, es importante conocer y evaluar el estado del recurso hdrico del Ro Katari
ya que este ro desfoga sus aguas al lago Titicaca en la Baha de Cohana, para poder proponer
medidas de mitigacin y control de la contaminacin. El uso de un modelo matemtico para
simular el Oxigeno Disuelto y la Demanda Bioqumica de Oxigeno nos permitir alcanzar lo
anteriormente mencionado de una forma econmica. (Ver Anexo E)
1.2. DEFINICIN DE OBJETIVOS
Los objetivos para el presente trabajo de investigacin son:
1.2.1. OBJETIVO GENERAL
o Aplicar el modelo SIMOD para estudiar la calidad del agua del Ro Katari en el tramo
conformado por la confluencia del Ro Pallina y Colorado, hasta la descarga a la
Baha de Cohana, Lago Titicaca.
1.2.2. OBJETIVOS ESPECFICOS
o Definir la segmentacin ambiental de la cuenca del Ro Katari.
o Identificar el rea de Influencia vinculado con la Calidad del Agua del Ro Katari.
o Determinar la calidad de agua del Ro Katari en el tramo final.
o Estudiar la capacidad de Autopurificacin del Ro Katari en el tramo final.
o Desarrollar las lneas maestras de un Plan de Accin Ambiental de control de la
contaminacin de las aguas del Ro Katari en trminos del tratamiento de las
descargas.
o Uso de los Modelos WASP 7.4 y RIOSep como comparacin con el Modelo SIMOD.
Pgina | 5
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 6
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Cantn Cascachi.
Cantn Augachi.
Municipio Pucarani:
Cantn Pucarani.
Cantn Lacaya.
Cantn Catavi.
Pgina | 7
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Municipio Laja:
Cantn Laja.
Cantn Tambillo.
o Provincia Ingavi
Municipio Viacha
o Provincia Murillo
Clima Fro, Subhmedo, el cual se presenta en las partes este y sur del rea. Su
particularidad principal es la fluctuacin de la precipitacin pluvial (entre 300 a
600 mm al ao) acompaada de bajas temperaturas.
2.4. TEMPERATURA
Las temperaturas medias anuales determinadas en las estaciones cercanas al rea de estudio
son de 11C en la estacin de El Alto y 7.7C en la estacin Tiwanaku (ubicada al Sur), se
debe al efecto trmico que ocasiona la presencia del Lago Titicaca.
Pgina | 8
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
2.7. SUELOS
Las caractersticas de los suelos son muy variables, dependiendo de la zona fisiogrfica en
que se encuentran, aunque, en trminos generales predominan suelos de transicin no bien
desarrollados.
En las llanuras fluvio-lacustres los suelos son de origen volcnico de textura
predominantemente franco-arenosa, con abundantes sedimentos, razn por la que la zona es
Michael Russel Osina Torrez
Pgina | 9
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
prcticamente de baja fertilidad, por lo que no se encuentra mucha diversidad tanto en flora
como en fauna y mucho menos en especies forestales. La erosin en este tipo de suelos es de
tipo laminar y progresiva.
En las cimas de las colinas y serranas, los suelos presentan abundante pedregosidad y soltura
superficial, en cambio en las montaas, los suelos son superficiales de textura franco
arenosa, estos mismos evitan el proceso erosivo acelerado en las pendientes.
2.8. HIDROGRAFA E HIDROGEOLOGA
2.8.1.
HIDROGRAFA
El Ro Katari es el principal cuerpo de agua de la cuenca, y discurre en direccin SuresteNoroeste; el Ro Seco confluye en el Ro Pallina, hasta la confluencia con el Ro Colorado,
donde cambia de nombre a Ro Katari. Tiene una longitud de 86 km, desde su naciente hasta
su desembocadura en el Lago Menor. Sus principales afluentes son: Ro Seco (que atraviesa la
ciudad de El Alto) y el Ro Pallina (que bordea la ciudad de Viacha y la poblacin de Laja).
2.8.2. RO SECO RO PALLINA
El Alto pertenece al rea de influencia del sistema hidrogrfico del Altiplano Norte;
conformado por los siguientes ros: Seco, Seque, Kantutani, San Roque y Hernani, que nacen
en la Cordillera de La Paz y/o en sus estribaciones.
2.9. MEDIO BITICO
La zona de estudio, pertenece a la ecoregin del Altiplano denominada puna semihmeda,
con una vegetacin caracterizada por una cubierta de gramneas duras dispuestas en macollos
y arbustos resinosos bajos.
El Altiplano norte puede dividirse en cuatro unidades fitofisiogrficas:
o Zona montaosa
o Zona Colinosa
o Zona de pie de monte
o Zona de llanura aluvial
Pgina | 10
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
FAUNA ACUTICA
Entre la fauna acutica presente en la zona, se pueden diferenciar diferentes especies de aves,
siendo la principal la choka que se caracteriza por presentar importantes poblaciones en
cercana a los totorales.
La pesca es una actividad practicada; siendo la fauna ictica ms comn, las diferentes
especies del gnero Orestia, especies tpicas del lago que conviven junto a otras dos especies
introducidas: la trucha y el pejerrey.
Pgina | 11
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 12
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 13
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
bajo riego, son la haba, arveja y otras semillas de hortalizas, de diversas variedades adquiridas
en el mercado local; y en reas muy pequeas, se est introduciendo alfalfa (variedades:
Ranger, Pampeana y Bolivia 2000). Entre otros cultivos podemos mencionar al tarwi,
papaliza, isao, y otros.
o CALENDARIO AGRCOLA
Las principales actividades dentro el calendario agrcola por tipo de cultivo se detalla en el
siguiente cuadro: (35)
Cultivo
Actividades
Meses
J
J
Preparacin del
terreno
Abonado
Siembra
Papa, cebolla y
Deshierbe
oca
Aporque
Tratamiento
sanitario
Cosecha
Siembra
Deshierbe
Haba y arveja
Riego
Cosecha
Siembra
Cebada en grano
Cosecha
Cebada y avena Siembra
para forraje
Cosecha
Tabla2-1. Calendario Agrcola, Cuenca del Ro Katari, La Paz, Bolivia.
Fuente: Prefectura del Departamento de La Paz Direccin de Recursos Naturales y Medio Ambiente. Primer Informe del Estudio Ambiental de la
Contaminacin Baha Cohana Lago Titicaca, PSID S.R.L. Noviembre de 2004.
Pgina | 14
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 15
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
3. FUNDAMENTO TERICO
3.1. CLASIFICACIN DE LOS CUERPOS DE AGUA
ORDEN
USOS
CLASE
A
CLASE
B
CLASE
C
CLASE
D
SI
NO
NO
NO
No
necesario
SI
NO
NO
No
necesario
No
necesario
SI
NO
No
necesario
No
necesario
No
necesario
SI
Pgina | 16
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
2
3
SI
SI
SI
NO
SI
SI
SI
NO
SI
SI
NO
NO
SI
SI
SI
SI
SI
SI
SI
NO
NO (*)
SI
SI
NO
NO (**)
SI
SI
SI
(SI) Es aplicable, puede tener todos los usos indicados en las clases correspondientes
(*) No en represas usadas para abastecimiento de agua potable
(**) No a navegacin a motor
(***) No aplicable a acuferos
Tabla3-1. Clasificacin de los cuerpos de agua segn su aptitud de uso.
Fuente: Cuadro N 1 - Anexo A del Reglamento en Materia de Contaminacin Hdrica.
Pgina | 17
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
El origen de la contaminacin de las aguas hay que buscarlo en procesos naturales pero, sobre
todo, en la actividad humana. La contaminacin inducida por el hombre es un fenmeno
antiguo si bien, la industrializacin y el crecimiento demogrfico y urbanstico han agravado
considerablemente el problema en muchas regiones. A pesar de su origen comn, se puede
establecer una distincin entre los principales problemas de contaminacin en aguas
superficiales y en aguas subterrneas. (27)
3.2.1. PRINCIPALES PROBLEMAS EN AGUAS SUPERFICIALES
o
Salinizacin. Se produce cuando la tasa de deposicin de las sales que transportan las
aguas de riego supera la velocidad de arrastre de la sal por el curso del agua. Es un problema
que se puede plantear debido a la reutilizacin de aguas residuales parcialmente tratadas,
puesto que estas aguas han sufrido un pequeo incremento de sales disueltas durante su uso
domstico o industrial. Es ms habitual que se produzca debido al riego con aguas de
acuferos sobreexplotados, con altas concentraciones de sales disueltas.
o
reaccin del vapor de agua con xidos de nitrgeno y azufre provoca la acidificacin de los
lagos. Tambin puede tener lugar por deposicin seca en las proximidades de los focos de
emisin de los contaminantes atmosfricos citados. El fenmeno se da en pases
industrializados y su correccin requiere el control y supresin de las emisiones a la atmsfera
de xidos de nitrgeno y azufre.
o
Pgina | 18
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
tambin provocar efectos altamente indeseables. Por ejemplo, la prdida de transparencia del
agua del mar en zonas costeras hace desplazarse a cierto tipo de algas desde su profundidad
normal a otras inferiores. De este modo quedan ms expuestas a la accin del oleaje y pueden
ser arrancadas de forma masiva, acumulndose en las playas. Su descomposicin puede
generar problemas de anoxia en las aguas costeras.
o
metales pesados son ingeridos por seres vivos a travs de la cadena alimenticia, acumulndose
en los tejidos grasos y musculares. Estos compuestos debilitan las defensas de los organismos
y provocan numerosas enfermedades, como ltimo eslabn, pueden alcanzar al ser humano
produciendo lesiones y hasta la muerte.
o
son vehculo de agentes infecciosos (protozoos, bacterias y virus) que pueden provocar un
gran nmero de enfermedades en los seres vivos (conjuntivitis, poliomielitis, hepatitis, tifus,
salmonelosis, clera, etc.). Este tipo de enfermedades de transmisin hdrica son todava muy
frecuentes en los pases del Tercer Mundo debido a la falta de condiciones de salubridad del
agua. La reutilizacin de aguas residuales tratadas para riego debe tener en cuenta la posible
presencia de estos agentes infecciosos. El vertido de aguas residuales a travs de emisarios
submarinos tambin puede ser otro foco potencial de enfermedades.
o
Pgina | 19
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
3
3
2 2 + 2 + 3
4 2 4
2 2
. (3 1)
La diferencia en el contenido de oxgeno al inicio y al final del anlisis, por litro de muestra
utilizada constituye la DBO de la muestra expresada en mg de O 2 L-1.
Lgicamente la velocidad y extensin a la cual ocurre esta degradacin depende de la
temperatura y del tiempo, de ah que se ha aceptado para estos dos parmetros 20 C y cinco
das respectivamente, tal como ya fue expresado.
Pgina | 20
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
En primer lugar se oxida la DBO carboncea (DBO C ), proceso al que sigue entre 8 y 12 das
despus, la oxidacin de la DBO nitrogenada (DBO N ).
Tanto la oxidacin de la DBO carboncea, como de la nitrogencea, son procesos de primer
orden, siendo la velocidad de oxidacin (equivalente a la velocidad de remocin de la DBO),
proporcional a la DBO presente: (23)
Donde:
Pgina | 21
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
como las velocidades, deben ser homogneas, lo cual puede hacer necesaria la conversin de
unidades de algunos parmetros.
En los anlisis de OD, la oxidacin de la DBO representa una prdida del OD que se produce
a la misma velocidad que la remocin de la DBO. La suma de las componentes asociadas a la
DBO carboncea y nitrogencea, da lugar a:
= +
Donde;
. (3 2)
= 1
. (3 3)
Pgina | 22
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 4)
. (3 5)
. (3 6)
. (3 7)
. (3 8)
. (3 9)
4.57( + 1 ) + 2
. (3 10)
Donde:
Pgina | 23
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
(3 11)
2
3
4
10.754 2.1407 103
1.7674 102
+
. (3 12)
[2 ] = 139.34411 +
= 1.80655 [ ]
Donde:
. (3 13)
. (3 14)
Pgina | 24
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Cl = Cloruros (mg/m3)
SC = Conductividad elctrica (S/cm) = CE
La disminucin de la concentracin de saturacin con la altitud es debida a la disminucin de
la presin atmosfrica y, por lo tanto, de la presin parcial de oxgeno en el aire. El clculo de
la concentracin se puede efectuar considerando que la presin atmosfrica desciende un 1%
por cada 82 metros que se asciende. Esta aproximacin da lugar a la siguiente expresin:
[2 ], = [2 ],0 8200
. (3 15)
. (3 16)
Por ltimo, la humedad relativa tambin, aunque en pequea medida, puede modificar la
concentracin de saturacin, dado que modifica las presiones parciales de los gases en la
atmsfera. (5)
Las principales fuentes de oxgeno en un sistema natural son las siguientes:
o Aportes desde afluentes. Se trata del oxgeno disuelto que entra en el sistema a travs
de diversos aportes. Se obtiene a partir de mediciones in situ de oxgeno disuelto y del
caudal.
o Reaireacin superficial. Cuando la concentracin de oxgeno disuelto en el agua es
inferior a la de saturacin, se produce una transferencia de oxgeno desde la atmsfera
hacia el seno del agua. En el caso, poco habitual pero posible, de que la concentracin
de oxgeno disuelto en el agua sea superior a la de saturacin, el flujo de oxgeno es
desde el agua hacia la atmsfera.
o Fotosntesis. La sntesis de materia orgnica a partir de sales minerales y absorcin de
energa por parte de algas y plantas verdes genera oxgeno en el sistema. (27)
Como principales consumos de oxgeno en un sistema natural se pueden considerar los
siguientes:
Michael Russel Osina Torrez
Pgina | 25
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 26
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 17)
Donde: (T) es la temperatura, (R) la constante de los gases ideales (3,314 J/mol/K) y (E) la
llamada energa de activacin.
Al integrar la ecuacin (3-17) entre los lmites de temperatura (T 1 ) y (T 2 ) se obtiene:
2 (2 1 )
=
1
(1 2)
. (3 18)
donde es posible obtener el valor de (K 2 ) deseado a partir del conocimiento de (K1 ) para una
temperatura dada y (E). Asumiendo que el factor (E/R T 1 T 2 ) es relativamente constante la
ecuacin (3-18) tomara la forma de:
2
= (2 1 )
1
2
= (21)
1
. (3 19)
. (3 20)
Pgina | 27
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 21)
Debe de considerarse que en algunos casos () vara con la temperatura lo que hace necesario
tomar valores de () vlidos para el rango de temperatura a utilizar en la ecuacin (3-21). (39)
3.4. MODELACIN MATEMTICA DE LA CALIDAD DE LAS AGUAS
3.4.1. DEFINICIN
Un modelo es una representacin simplificada o abstraccin de la realidad, cuyo objetivo es el
pronosticar alguna faceta del comportamiento del sistema sometido a estudio. En nuestro caso
se trata de la representacin de hechos en la naturaleza mediante el desarrollo de tcnicas y
herramientas para expresar en trminos matemticos los diversos aspectos fsicos, qumicos y
biolgicos de los cambios experimentados por sustancias agregadas al recurso hdrico.
Los modelos matemticos pueden representar muchos problemas complejos de modo
eficiente y conciso, y en muchos casos, pueden ser la manera ms barata de analizar los
problemas. Facilitan la experimentacin por que todas las variables dependientes, las
variables independientes, las constantes y los parmetros estn explcitamente relacionados
por el lenguaje de la matemtica. (7)
3.4.2. FINALIDAD DE LOS MODELOS DE CALIDAD DE AGUAS
Los modelos de calidad de agua tienen por finalidad:
o Determinar las nuevas concentraciones del cuerpo de agua en cada punto y a lo largo
del lapso de inters cuando las condiciones de modificacin y el estado primitivo son
conocidos.
o Determinar el nivel de impacto que ejerce sobre un cuerpo de agua por descargas.
o Simular situaciones extremas acerca de los cuales se tiene poca informacin con el
principal objetivo de estar preparados para la eventualidad.
Pgina | 28
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
o Fijar los puntos de control y niveles de calidad a fin de evitar que se produzcan
condiciones anaerobias en el futuro y en su caso determinar las acciones a tomarse.
La principal ventaja en el empleo de modelos es que pueden ser aplicados repetidas veces con
diferentes conjuntos de datos de entrada, para simular diferentes condiciones, en el sistema
estudiado, lo cual permite predecir los efectos de varios planes y polticas en el sentido
espacial y temporal.
Los modelos de calidad de agua determinan las nuevas concentraciones del cuerpo de agua en
cada punto y a lo largo del lapso de inters cuando las condiciones de modificacin y el
estado primitivo son conocidos. El Grafico 3-3 precisa este concepto, suponiendo un punto
cualquiera. Se busca el valor de C para un incremento t siendo t el tiempo y C la
concentracin. Como se ver t puede ser cualquiera, desde un segundo a un ao. (6)
PROCESOS ACTUANTES
Pgina | 29
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
generacion de masa
perdida de masa
aumento de masa
entrada de masa
salida de masa
Ec. (3 22)
Todos los trminos de la ecuacin (3-22), se expresan en dimensiones de, masa por unidad de
tiempo [MT-1]. La ecuacin de conservacin de la masa puede aplicar, tanto si el vertido es a
un lago, a una corriente de agua, o a una zona costera. Sin embrago, las diferentes
caractersticas fsicas de estos cuerpos de aguas receptores obligan a adoptar enfoques y
aproximaciones diferentes en cada caso para obtener las concentraciones de constituyentes. (31)
3.4.3.1.1.
PROCESOS DE TRANSPORTE
DISPERSIN
El trmino dispersin hace referencia a los procesos de transporte que sufren las sustancias
disueltas y/o particuladas y que no estn directamente relacionados con el desplazamiento de
las masas de agua. Hay tres procesos que contribuyen a dicho transporte:
o Difusin molecular. La difusin molecular es el desplazamiento aleatorio de las molculas
de solutos entre las molculas de agua debido al movimiento browniano de las partculas.
En la prctica se observa una homogenizacin de las partculas en la solucin de manera
que el efecto es como si las molculas se desplazaran desde las zonas de mayor
concentracin a las de menor, por lo que siempre que exista un gradiente de concentracin,
se dar este fenmeno. Es un proceso extraordinariamente lento que no suele ser
importante a efectos prcticos en el estudio de sistemas naturales. (30)
o Difusin turbulenta. La difusin turbulenta hace referencia al transporte de sustancias
disueltas y pequeas partculas causada por turbulencias a pequea (micro) escala. Se trata
Michael Russel Osina Torrez
Pgina | 30
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 23)
. (3 24)
Donde:
. (3 25)
(3 26)
Pgina | 31
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Haciendo coincidir los ejes coordenadas con los ejes principales del tensor de dispersin, solo
quedan tres componentes no nulas de este, las correspondientes a E xx , E yy y E zz o bien E x, E y
y E z, que son los trminos de la dispersin que aparecen en la ecuacin general del balance de
materia.
=
+
+
(3 27)
ADVECCIN
Utilizando el volumen de control cubico infinitesimal del Grafico 3-5, los trminos de la
ecuacin de conservacin de la masa relativos a la adveccin en la direccin x, se pueden
escribir como: (6)
Velocidad de aumento de masa en el volumen de control =
. (3 28)
. (3 29)
Pgina | 32
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
. (3 30)
Grafico3-5. Masa de un constituyente que entra y sale de un volumen de control debido a la adveccin en direccin x.
Fuente: Metcalf&Eddy. 1995. Ingeniera de aguas residuales: Tratamiento vertido y reutilizacin, Tercera Edicin, Ed. McGraw-Hill, Inc., Espaa.
. (3 31)
. (3 32)
Pgina | 33
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
=0
en la ecuacin (3-32), permite llevar a cabo una simplificacin final. En la ecuacin que
resulta, se define el defecto de la adveccin en los cambios de las concentraciones con el paso
del tiempo:
3.4.3.1.2.
PROCESOS DE TRANSFORMACIN
. (3 33)
PASOS DE TRANSPORTE
Pgina | 34
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Cada uno de estos pasos o de transferencia implica un incremento finito de tiempo. Sin
embargo, las condiciones sern tales que un paso particular ser tan largo con relacin a los
otros que las velocidades de los pasos restantes sern insignificantes en el proceso global.
Para unas condiciones dadas, el paso ms largo se llama el paso que determina la velocidad
y los clculos de velocidad de transferencia global pueden basarse en este paso solamente.
Grafico3-6. Representacin esquemtica de la transferencia de O2 en cuatro pasos desde la fase vapor a la fase solucin.
Fuente: Weber, W. J. 1979. Control de calidad del agua: Procesos Fisicoqumicos. Universidad de Michigan. Editorial Reverte, S. A. 1979. Edicin en espaol.
Capitulo 11. Espaa
El aire es una mezcla de gases, compuesto principalmente por nitrgeno y oxigeno, con
menores porcentajes de anhdrido carbnico y gases nobles. Ninguno de los gases, con
excepcin del anhdrido carbnico, reacciona qumicamente en el agua; as sus solubilidades
en el agua se definen a partir de la ley de Henry, la que establece para soluciones diluidas en
Michael Russel Osina Torrez
Pgina | 35
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 34)
Donde:
= 298.15
1
1
298.15
. (3 35)
. (3 36)
Pgina | 36
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Grafico3-7. Relaciones entre la presin parcial del O 2 en la fase gaseosa y la concentracin en la fase disolucin.
Fuente: Weber, W. J. 1979. Control de calidad del agua: Procesos Fisicoqumicos. Universidad de Michigan. Editorial Reverte, S. A. 1979. Edicin en espaol.
Capitulo 11. Espaa
El proceso de transferencia de O 2 , desde la fase gaseosa a la fase disolucin esta descrito por
la simple reaccin: (48)
2() 12 2()
. (3 37)
. (3 38)
Pgina | 37
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 39)
= 2,
. (3 40)
. (3 41)
. (3 42)
+
=
+ 1
Entonces:
1
=
+
1
1
1
=
+
. (3 43)
. (3 44)
. (3 45)
Las ecuaciones (3-44) y (3-45) expresan los coeficientes de transferencia globales en funcin
de las resistencias de fase individuales y la constante de la ley de Henry.
Pgina | 38
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Las velocidades de disolucin para gases altamente solubles en agua, como el amoniaco, son
controladas por la velocidad de transferencia a travs de la pelcula de gas porque, para estos
sistemas KH es grande. La ecuacin (3-45) indica que el coeficiente de transferencia global es
entonces prcticamente igual al coeficiente de transferencia de la fase gaseosa. En este caso,
puede aumentarse la transferencia por agitacin de la fase gas con el fin de reducir el espesor
de la pelcula de gas. Por el contrario, la transferencia de gases poco solubles en el agua,
tales como el oxigeno, nitrgeno y dixido de carbono, se controlan generalmente por la
velocidad de transferencia a travs del lado liquido de la interface, ya que para estos sistemas
K H es pequea y como indica la ecuacin (3-44), la mayor resistencia est en la fase liquida.
Si el espesor de la pelcula liquida se reduce, en este caso por agitacin de la fase liquida, la
velocidad de transferencia global del gas es mayor. La transferencia de gases de solubilidad
intermedia en el agua, tal como el sulfuro de hidrogeno, viene controlada por una
combinacin de las resistencias intefaciales del gas y liquido y es preciso una agitacin de
ambas fases con el fin de aumentar la transferencia gaseosa.
Realmente la velocidad global de transferencia de un gas a una disolucin viene controlada
por la velocidad de desplazamiento del gas a travs de uno de los diversos planos de
transferencia hipotticos. Independientemente de que plano ofrece la mayor resistencia al
desplazamiento del gas, y por tanto limita el proceso de transferencia global, la velocidad ser
siempre una funcin de la diferencia de contraccin gaseosa a travs de este plano, es decir el
gradiente de concentracin. Adems, el valor limitante del gradiente de concentracin es una
funcin de la solubilidad o limite de saturacin para un gas determinado.
3.4.3.1.2.3.
En forma ms detallada los conceptos relacionados con los temas indicados anteriormente se
presentan a continuacin.
Modelo de la pelcula: fue propuesto por Nernst en 1904 y desarrollado posteriormente por
Withman en 1924, y supone una relacin entre el coeficiente de masa y la difusividad
molecular de la especie que se trasfiere (D m l), como la que sigue:
M =
. (3 46)
Donde es el espesor efectivo de la capa estancada en la fase liquida. Este espesor constituye
el nico parmetro que debe estimarse en base a datos experimentales para calcular el
Michael Russel Osina Torrez
Pgina | 39
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Si consideramos un paquete de fluido que llega a la interface desde el seno del lquido,
como se muestra esquemticamente en la figura, el elemento que originalmente tena una
Michael Russel Osina Torrez
Pgina | 40
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
= 2
0.5
. (3 47)
. (3 48)
0.5
. (3 49)
0.5
= 0.34
. (3 50)
Donde f es la frecuencia de los remolinos en la regin cercana a la interface gas lquido. (45)
3.4.3.1.2.4.
REAIREACIN NATURAL
En los casos en los que la concentracin de Oxigeno Disuelto (OD) en un cuerpo de agua con
una superficie libre expuesta a la atmosfera, se encuentra por debajo de la concentracin de
saturacin, se produce un flujo neto de oxigeno molecular entre la atmosfera y el agua. Este
Michael Russel Osina Torrez
Pgina | 41
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
( ) = ( ) = 2 ( )
. (3 51)
Pgina | 42
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
2 =
Donde:
( )0.5
1.5
. (3 52)
Donde:
. (3 53)
(T-
20C)
subtitulo 3.4.8.2.
3.4.3.1.2.5.
Pgina | 43
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Despus de un corto perodo de adaptacin al nuevo ambiente (fase lag) las bacterias se
comienzan a reproducir por fisin binaria, incrementando su nmero y su biomasa en forma
exponencial. La existencia de alimento (sustrato) en cantidad suficiente proporciona las
condiciones adecuadas para una tasa de crecimiento mxima. La tasa de metabolismo en
la fase de crecimiento exponencial (exponential growth phase) est limitada nicamente por
la capacidad de las bacterias de procesar el alimento. Si X representa la concentracin de
biomasa y m es una constante de proporcionalidad, la tasa de crecimiento de la biomasa se
puede escribir como:
Donde
. (3 54)
La fase de crecimiento en declinacin (declining growth phase) es causada por una reduccin
sostenida del alimento o sustrato. La tasa de reproduccin decrece hasta que el nmero de
bacterias existentes es estacionario, lo que ocurre cuando la tasa de reproduccin iguala a
la tasa de mortalidad. En experimentos de laboratorio se estudi el crecimiento de
bacterias en cultivos artificiales. Como resultado de estas
el
experiencias se
encontr que
Donde:
. (3 55)
Pgina | 44
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
==
+
. (3 56)
. (3 57)
Pgina | 45
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Dado que el crecimiento est limitado por el consumo de todo el sustrato, se puede asumir que
Sm=0, con lo que se tiene:
=
. (3 58)
. (3 59)
Donde:
=
( + )
. (3 60)
Pgina | 46
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde
. (3 61)
=
=
+
+
. (3 62)
> (3 63)
=
= = (3 64)
Pgina | 47
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
L= Concentracin del sustrato o materia orgnica = S (mg/l)
K = tasa de crecimiento mximo = m
K ox = Constante de oxidacin o desoxigenacin = K d
3.4.3.2.
+ +
(3 65)
Donde:
I = velocidad de exportacin de las fuentes externas; masa introducida por unidad de tiempo y
de volumen de agua MT-1L-3
i = ndice del proceso de transformacin.
J = ndice de identificacin de la fuente externa
Como ejemplo, en el caso en que el constituyente considerado es el OD, el trmino relativo a
las transformaciones es:
= + + (3 66)
Michael Russel Osina Torrez
Pgina | 48
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
r P = velocidad de produccin de oxigeno por fotosntesis por unidad de tiempo y de volumen
de agua, MT-1L-3
r RP = velocidad de consumo de oxigeno en la respiracin por unidad de tiempo y de volumen
de agua, MT-1L-3
Muchos de los constituyentes presentes en el agua residual, como la DBO, los slidos
suspendidos, o las bacterias, solo estn sometidos a un proceso de transformacin, que sigue
principios
cinticos
de
primer
orden.
En
este
caso,
Por lo tanto el anlisis de varios de los constituyentes se puede realizar con nica forma de la
ecuacin de conservacin de la masa, para lo cual es necesario ajustar adecuadamente el valor
de K.
La ecuacin de conservacin de la masa es la base de prcticamente, la totalidad de los
modelos de calidad de agua. Es una ecuacin diferencial en diversas parciales de segundo
orden, cuya resolucin en el caso general, resulta complicada. Sin embargo se pueden realizar
numerosas simplificaciones para las cuales si existe una solucin analtica. En primer lugar, se
trata de una ecuacin lineal, de modo que la suma de soluciones (Como funcin de x, y, z y t)
tambin es solucin. Por ejemplo, las concentraciones de constituyentes correspondientes a
dos descargas independientes que tienen lugar simultneamente, son las suma de las
concentraciones correspondientes a cada una de las descargas suponiendo que no se producen
de forma simultnea. Adems, si la concentracin de un constituyente en el efluente se
aumenta o reduce en un determinado factor, las concentraciones en cualquier punto del cuerpo
de aguas receptor aumentan o disminuyen en la misma relacin. Ntese que, si existe una
concentracin por encima de la concentracin de base. Segundo, el campo de velocidades en
toda la zona afectada en funcin del tiempo, la determinacin del campo de velocidades se
puede llevar a cabo mediante mediciones o modelos.
3.4.3.2.1.
Pgina | 49
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Ec.(3-67)
Pgina | 50
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Una vez que se formula el modelo matemtico que es una ecuacin diferencial o un
sistema de ecuaciones diferenciales. Por supuesto al incrementar la resolucin, se
agrega complejidad al modelo matemtico y se incrementa la probabilidad de que no
se obtenga complejidad al modelo matemtico y se incrementa la probabilidad de que
no se obtenga una solucin explicita.
En un modelo matemtico de un sistema fsico, suele intervenir la variable tiempo t.
entonces una solucin del modelo da el estado de sistema; en otras palabras los valores
de la variable dependiente (o variables) para valores apropiados de t describen el
sistema en el pasado, presente y futuro. (7)
3.4.3.2.2.
EXACTITUD Y PRECISIN
Los errores asociados con los clculos y medidas se caracterizan observando su exactitud y
precisin. La exactitud se refiere a que tan cercano esta el valor calculado o medido con el
valor verdadero. La precisin de refiere que tan cercano esta un valor individual medido o
calculado con respecto a los otros
La inexactitud (conocida tambin como sesgo) se define como un alejamiento sistemtico de
la verdad. La imprecisin (tambin llamada incertidumbre) se refiere a la magnitud del
esparcimiento de los datos. (22)
3.4.4. DEFINICIN DE TRANSMISIN LONGITUDINAL DE CONTAMINACIN
Pgina | 51
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Efecto advectivo:
Entrada:
. (3 68)
Salida: ( + )( + )
( + ) +
Donde:
. (3 69)
. (3 70)
Q = Caudal (m3/s)
C = Concentracin del contaminante (g/m3)
t = Diferencial de tiempo (s)
o
negativo.
Entrada:
. (3 71)
Pgina | 52
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Salida: ( + )( + )
. (3 72)
Donde:
E = Efecto dispersivo en (L2/T)
A = Superficie (m2)
o
+( + )( + )
Donde:
+
+ () +
. (3 73)
2
+
2
. (3 74)
+ 2 +
. (3 75)
Pgina | 53
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
1
2 1
1
=
+ 2
+
+
. (3 76)
Simplificando:
1
2
1 1
=
+ 2 +
+ 1
. (3 77)
2
1
=
+ 2 +
+ 1
. (3 78)
Si se asume que:
Donde:
2
=
+ 2 +
Ordenando:
. (3 79)
. (3 80)
. (3 81)
Influencia
de laadvectivo
reaccin
Influencia del
transporte
Influencia de la dispersin
3.4.5. SISTEMAS ACOPLADOS
Pgina | 54
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
sistema sus elementos interactan, ms all de la frontera existen factores exgenos que
incluyen en los elementos del sistema pudiendo llegar a producir perturbaciones. (42)
Factores Exgenos
Frontera
Entrada
Salida
Ncleo
Pgina | 55
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
plomo, mercurio, etc. Estos componentes (dando por sentado que no existen fuentes
antropognicas) se introducen en el ro a travs de las aguas subterrneas. Se da por supuesto
que la contribucin de las aguas subterrneas al flujo del ro no cambia de forma considerable
(incluso en pocas de riada) estas sales disueltas se mantendrn, por lo tanto, constantes a lo
largo del ao.
Como los parmetros analizados es en te trabajo son no conservativos, no se proceder a
detallar la modelacin de los componentes conservativos, sin embargo la explicacin de la
misma se puede encontrar en la bibliografa presentada.
Los contaminantes conservativos se calculan aguas abajo en condiciones de mezcla completa
de acuerdo a la ecuacin:
(1 1 ) + (2 2 )
(1 + 2 )
. (3 82)
La ecuacin de mezcla se utiliza siempre como la primera opcin del balance de masa, su
estructura permite por un simple despeje estimar un caudal si conocemos los gastos, y una
concentracin. (39)
3.4.7. MODELACIN DE LOS COMPONENTES NO CONSERVATIVOS
Los componentes no conservativos son aquellos qumicamente o bioqumicamente reactivos.
La concentracin vara durante el transporte. Usualmente se asume para el cambio una
cintica de primer orden. (39)
3.4.7.1.
MODELACIN DEL OD
Consideramos un volumen de control del ro con entradas y salidas como las mostradas en la
Grafica 3-17.
Michael Russel Osina Torrez
Pgina | 56
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
+ + 01 = 02 =
. (3 83)
Pgina | 57
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
R 02 = constante de reaireacin
Tomando el limite como x0 y para el caso en rgimen estacionario
1 + 2 = 0
= 0
. (3 84)
Sin embargo, t= V/Q por lo que la anterior ecuacin se convierte en la ecuacin diferencial
para el modelo simple de la curva de dficit de oxigeno:
= 1 + 2
. (3 85)
Donde:
C = Concentracin de Oxigeno Disuelto
Cs = Concentracin de Oxigeno Disuelto de saturacin
Por tanto:
= 1 + 2
Po tanto;
. (3 86)
. (3 87)
1 = 1 1
. (3 88)
+ 2 = 1 1
. (3 89)
1
( 1 2 ) + 2
2 1
. (3 90)
Pgina | 58
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
D(t) = Dficit de oxigeno en cualquier momento D = Cs-C(t).
Lo = DBO ultimo en el lugar de descarga.
Do = Dficit inicial de oxigeno DBO en el lugar de descarga.
La anterior ecuacin es la ecuacin del dficit de oxigeno de Streeter and Phelps, y se
representa en el Grafico 3-18.
Grafica3-18: Curva caracterstica del agotamiento de OD obtenida mediante la aplicacin de Streeter and Phelps.
Fuente: River Pollution and the Dissolved Oxygen Sag Curve, Lecture 9, Spring 2000.
El dficit mximo de oxigeno Dc, que tiene lugar a una distancia Xc aguas del punto de
descarga se puede calcular estableciendo que dD/dt = 0, es decir;
= 1 1 2 = 0
1
1
2
. (3 91)
. (3 92)
Para determinar el tiempo t, en el que esto se produce la ecuacin principal se despeja con
respecto al tiempo y se resuelve para t:
()
1
( 1 1 2 ) 2 2 = 0
=
2 1 2
=
1
2
(2 1 )
1
2 1
1
1
. (3 93)
. (3 94)
Pgina | 59
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
3.4.7.2.
MODELACIN DE LA DBO
Grafica3-19: Comparacin del comportamiento de la DBO en el ro y laboratorio a) decaimiento tpico de DBO en el ro y b) determinacin de la DBO en el
laboratorio
Fuente: U. S. Environmental Protection Agency (EPA): Streams and Rivers, Technical Guidance Manual for Developing Total Maximum Daily Loads, 1995,
20-45.
La DBO, ser removida como resultado de la biodegradacin aerbica y anaerbica, segn las
caractersticas hidrodinmicas del flujo.
Como el oxigeno es equivalente a la materia orgnica se tiene:
Donde:
. (3 95)
L = Materia orgnica
T = Tiempo de incubacin
Pgina | 60
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 96)
=0
Donde:
. (3 97)
Durante la prueba en el laboratorio se mide cuanto oxigeno remanente queda y por diferencia
se tiene el oxigeno que se consumi. Es decir no se realiza medicin de materia orgnica en
laboratorio sino de oxigeno consumido durante el proceso.
Pgina | 61
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (3 98)
= (1 )
Donde:
. (3 99)
6
5
Conc. (mg/L)
4
3
2
L(t) remanente
BOD(t) ejercida
0
0
10
15
20
25
30
Tiempo (das)
En trminos de DBO:
= +
= +
Michael Russel Osina Torrez
Pgina | 62
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
DBO ultima : es la cantidad de oxigeno necesaria para oxidar toda la materia orgnica, presente
en un agua, por medio de procesos bioqumicos; por lo tanto es la cantidad de oxigeno
equivalente a toda la materia orgnica.
Es conveniente expresar los resultados de la DBO 5 a la DBO total en cualquier tiempo dado
por extrapolacin, de acuerdo a la siguiente relacin: (6)
5 = 1 5
. (3 100)
TASA DE REAIREACIN (K 2 )
Pgina | 63
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
= 3
. (3 101)
Pgina | 64
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 65
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 66
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
3.4.9. CALIBRACIN
Para que un modelo pueda ser aplicado en predecir condiciones o escenarios futuros de
planificacin, debe ser capaz de reproducir aceptablemente las condiciones actuales, siendo
este entonces el objetivo de la calibracin
La calibracin del modelo se efecta usando uno o ms conjuntos de datos tanto de entrada
como de salida. Los parmetros del modelo y de hecho el modelo en si se ajustan y modifican
de forma que se produzca un resultado que sea tan cercano posible a la calidad de agua
observada en el campo. La evaluacin de los parmetros del modelo, que suele realizarse por
un procedimiento subjetivo de prueba y error, pueden ser ayudada por medio de mtodos de
optimizacin, y el proceso es llamado identificacin de los parmetros. En general se usan
como puntos de partida para el ajuste de los parmetros los valores de las constantes de
Pgina | 67
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 68
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 69
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
= ( )
( )
( )
Modelos estocsticos: Los modelos estocsticos se caracterizan por tener alguna componente
desconocida, de la cual slo se conoce su probabilidad de ocurrencia. Una forma clsica de
representar estas variables es la siguiente:
= +
Pgina | 70
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
de la Facultad
de Ingeniera
(DEPFI
herramienta tcnica, a la Primera Fase del Estudio de Calidad del Agua en el Ro Bravo,
patrocinado por el Instituto Mexicano de la Tecnologa del Agua (IMTA) de la Secretara
de Agricultura y Recursos Hidrulicos
4.1.2. RIOSep
El programa RIOSep es un modelo matemtico de calidad de agua desarrollado por el Centro
Panamericano de Ingeniera Sanitaria y Ciencias del Ambiente (CEPIS) reemplazando a
RIOSII, RIOSIII, RIOSIV, SIMOX.
Realiza el anlisis en estado permanente de oxigeno disuelto en ros, considerando que la
dispersin no es significativa en comparacin con el transporte advectivo. Evala fuentes
puntuales o dispersas de contaminantes. (8)
Pgina | 71
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 72
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Dficit de oxigeno.
o Procesos de separacin
Cintica
(precipitacin-disolucin,
adsorcin-desorcin,
sedimentacin-
resuspensin)
o Mdulos
M
Michael Russel Osina Torrez
Pgina | 73
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Euler
COSMIC
4.3.LIMITACIONES
4.3.1. SIMOD 1.01
o El modelo SIMOD ha sido desarrollado para modelar condiciones aerbicas y no as
condiciones anaerbicas. (43)
o El empleo del modelo SIMOD no tiene limitaciones en rango de valores permitidos.
4.3.2. RIOSep
o El nmero de tramos limite es de 100 segmentos
o En cada tramo analizado solo se puede considerar una descarga, a mayor numero de
descargas en el ro (domesticas, industriales y/o tributarios) se tiene mayor nmero de
tramos.
o No se puede considerar caudales negativos o extracciones
o Se debe calcular manualmente los valores de K 1 (8)
Pgina | 74
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 75
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
( )
=
+ ()
Donde:
. (4 1)
M = Masa (M)
X = Distancia (L)
t = Tiempo (T)
C = Concentracin (M/L3)
A = rea de la seccin transversal (L2)
DL = Coeficiente de dispersin (L2/T)
U = velocidad media (L/T)
Michael Russel Osina Torrez
Pgina | 76
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
=
=
+
. (4 2)
Donde: =
= 0, por tanto:
= 0 y la ec. (4 2) queda:
(4 3)
( )
=
(4 4)
4.4.2.2.RIOSep
(4 5)
(4 6)
= 0 =
+
. (4 7)
Pgina | 77
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
= 0 =
+
. (4 8)
= 0 =
+ + + ( )
Donde:
D=
L=
N=
Ka =
Kr =
Kd =
Krn =
Kn =
WL =
. (4 9)
Pgina | 78
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
4.4.2.3.WASP 7.4
Una ecuacin de balance de masa para constituyentes disueltos en un cuerpo de agua, debe
tomar en cuenta toda la materia que entra y sale a travs de la carga directa y difusa;
transporte advectivo y difusivo; transformacin fsico, qumico, y biolgico. Considera el
sistema de coordenadas mostrada en la ecuacin 5-1, donde x y y son coordenadas en el
plano horizontal y z es coordenada en el plano vertical.
( ) +
= ( )
+ + + + (5 1)
. (4 10)
Donde:
dia-1)
Pgina | 79
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
( ) =
+ + ( + ) +
. (4 11)
Profundidad media: = 2 2
. (4 12)
. (4 13)
Donde A1 , B 1 , A2 , B 2 , son constantes que son calculadas a partir de datos de escala, gasto y
. (4 14)
Pgina | 80
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
Q = caudal
V = Velocidad
A = rea de la seccin transversal
4.5.3. WASP 7.4
La ecuacin de continuidad para una dimensin es:
Donde:
Q = es el caudal (m3/s)
+
=0
. (4 15)
Donde:
+
=0
. (4 16)
Donde:
= 0
. (4 17)
Pgina | 81
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
S f = Friccin de la pendiente
La ecuacin de Manning expresa la fuerza de friccin como funcin de la velocidad del agua y su
radio hidrulico:
Donde:
2 2
. (4 18)
n = Factor de friccin
v = Velocidad del agua (m/s)
R = Radio hidrulico (m) equivalente a la profundidad media transversal, D
De la ecuacin de momento simplificada, S O puede igualarse con S f. El radio hidrulico
puede
Puede ser expresado como el rea transversal dividido por el ancho, B (m). Sustituyendo esto
en la ecuacin de Manning:
5
1 3 12
=
23
. (4 19)
+ 1
=0
. (4 20)
Las ecuaciones que relacionan ancho, profundidad y velocidad en el flujo de agua son:
=
Donde:
. (4 21)
. (4 22)
. (4 23)
Pgina | 82
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (4 24)
= = = ( ) ( ) ( ) = ( ) ++
=1
++ =1
. (4 26)
Ec. (4 25)
. (4 27)
El modelo WASP 7.4 solo requiere especificaciones para las relaciones de velocidad en la
ecuacin (4-21), y profundidad (4-22) los coeficientes para la ecuacin (4-23) esta
implcitamente especificado por las ecuaciones (4-26) y (4-27). (38)
4.6.CONSTANTES O TASAS
4.6.1. TASA DE REAIREACIN
4.6.1.1.SIMOD 1.01
La tasa de reaireacin se denomina en este modelo Ka presentando las siguientes opciones:
o Churchill et al (1962)
o O`Connor et al (1958)
o Bennett & Rathbun (1972)
3
Pgina | 83
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
. (4 28)
. (4 29)
Donde:
. (4 30)
Donde:
86400 3 2 3
=
(100 )
100
86400
1 1
(1 100 )1 + 2100
100
2
3 2 3
1 =
u
1
2
2 =
0
2
86400
=
100
100
. (4 31)
. (4 32)
. (4 33)
. (4 34)
. (4 35)
Pgina | 84
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
K wj = Constante de reaireacin inducida por el viento (dia-1)
W = tiempo de variacin de la velocidad del viento a 10 cm sobre la superficie (ms-1)
Ta = Temperatura del aire (C)
a = Densidad del aire en funcin de Ta (g*cm-3)
w = Densidad del agua (g*cm-3)
va = Viscosidad del aire en funcin de Ta (cm2s-1)
v W = Viscosidad del aire en funcin de Ta (cm2s-1)
DOW = Difusividad del oxigeno en el agua en funcin de T (cm2s-1)
D j = Coeficiente de Von Karman 0.4
vT = Transicin del esfuerzo cortante, establecido 9, 10, y 10 para pequeas, medianas y
grandes escalas (cm*s-1)
v c = Transicin del esfuerzo cortante, establecido 22, 11, y 11 para pequeas, medianas y
grandes escalas (cm*s-1)
z e = rugosidad equivalente, establecido 0.25, 0.35, y 0.35 para pequeas, medianas y grandes
escalas (cm)
z 0 = rugosidad efectiva, en funcin de z e , , C d , V T , va , W (cm)
= inverso del numero de Reynold, establecido 10, 3 y 3 5 para pequeas, medianas y
grandes escalas
= Coeficiente adimensional, establecido 10, 6.5, y 5 para pequeas, medianas y grandes
escalas
u = Coeficiente adimensional, en funcin de , vc . C d
C d = Coeficiente de arrastre en funcin de z e , , va , D j , vT , W
Michael Russel Osina Torrez
Pgina | 85
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
(4 36)
Pgina | 86
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
o Constante de reaireacin
= 20 1.02420
(4 37)
Los valores obtenidos con las ecuaciones 4-36 y 4-37 estn dados en (1/da) con base e.
4.7.2. RIOSep
Los coeficientes de reacciones se corrigen para temperatura de la siguiente manera: (8)
Donde:
= 20 20
. (4 38)
Ka
1.024
Kd
1.047
Kn
1.047
SB
1.065
Donde:
() = (20) 20
. (4 39)
Pgina | 87
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
4.8.REACCIONES O INTERRELACIONES
4.8.1. SIMOD 1.01
La variacin de los constituyentes
=
=
(4 40)
= ( )
(4 41)
Para resolver las ecuaciones (4-40) y (4-41) se tomo las siguientes consideraciones de
frontera:
X=0
L = Lo
C = Co
Pgina | 88
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
=
0
= +
. (4 42)
Donde:
D = C S -C
(Dficit de OD)
D O = C S -C O
. (4 43)
. (4 44)
Pgina | 89
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
4.8.2. RIOSep
Las soluciones generales son: (8)
=
=
. (4 45)
. (4 46)
. (4 47)
1
+
. (4 47)
Donde:
. (4 47)
(4 47 )
. (4 47)
. (4 47)
. (4 47)
. (4 48)
5 = 20 5
1 = 20 1
Michael Russel Osina Torrez
3
(1 5 )5
3
(1 1 )1
. (4 49)
. (4 50)
Pgina | 90
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
6 = 2 220 ( 6 ) 20 5
Donde:
64
20
20 1
14
. (4 51)
. (4 52)
Pgina | 91
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
DBO: Lm =
O.D.: Cm =
(LrQr + LdQd)
Qm
(CrQr + CdQd)
Temperatura: Tm =
Qm
(TrQr + TdQd)
Qm
. (4 53)
. (4 54)
. (4 55)
Donde:
Qm = Gasto despus de la mezcla (L3 /T)
Qr = Gasto del ro (L3 /T)
Qd = Gasto de la descarga (L3 /T)
Lm = DBO despus de la mezcla (M / L3)
Lr = DBO del ro (M / L3)
Ld = DBO de la descarga (M / L3)
Cm = O.D. despus de la mezcla (M / L3)
Cr = O.D. del ro (M / L3)
Cd = O.D. de la descarga (M / L3)
Tm = Temperatura despus de la mezcla (C)
Tr = Temperatura del ro (C)
Td = Temperatura de la descarga (C)
o Concentracin de saturacin del OD: varias en funcin de la temperatura, presin
atmosfrica, concentracin de cloruros e impurezas.
El modelo emplea la expresin desarrollada para el clculo de la concentracin del OD
a presin por los investigadores Elmore y Hayes (1960):
, = 14.62 0.3898 + 0.0069692 0.00005893
. (4 56)
Pgina | 92
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
= (1 0.00002287675)5.167 ,
. (4 57)
En casos en que debido a la descarga el Oxigeno Disuelto sea igual a cero, automticamente
RIOSep estima los niveles de DBO C y DBO N en condiciones anaerbicas y determina la
probable ubicacin de la recuperacin del Oxigeno Disuelto.
Las ecuaciones bsicas son las siguientes:
0 =
0 =
0 =
Donde:
. (4 58)
. (4 59)
. (4 60)
Ks = Kr - Kd + Kan
Kan = Descomposicin anaerbica de DBO C
Ksn = Krn - Kn
Michael Russel Osina Torrez
Pgina | 93
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
= +
. (4 61)
. (4 62)
. (4 63)
= +
. (4 64)
. (4 65)
. (4 66)
Se tiene que notar que dentro de tramos anaerbicos se aplican las siguientes condiciones y/o
limitaciones:
1. Ninguna fotosntesis/respiracin de algas
2. Ninguna demanda bntica de oxigeno disuelto
Hay varios fenmenos que pueden modificar el fenmeno de nitrificacin tal como se
presenta a continuacin:
o Los niveles de oxigeno disuelto bajos (<1.5 mg/l) suprimen la nitrificacin
o Las temperaturas bajas de agua (<10C) suprimen la nitrificacin
o El uso de amonio como fuente de nutrientes por algas
o El pH optimo para nitrificacin es de 7.0 a 8.5
Pgina | 94
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
= ( ) > 0
. (4 67)
Reordenando:
= ( ) < 0
. (4 68)
. (4 69)
Para la solucin de este caso se mantiene constante el nivel; de oxigeno disuelto permitiendo
la oxidacin completa de la DBO C , y, una vez suprimidos los trminos de demanda bntica,
fotosntesis y respiracin, se permite la oxidacin parcial de la DBO N , hasta que la
disminucin o aumento de esta conduzca a un punto en que dos tasas de variacin de oxigeno
disuelto tengan el mismo signo. (8)
= 0 = ( ) +
+
Entonces:
= ( ) +
. (4 70)
. (4 71)
Pgina | 95
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Para Ksn = 0
= +
1 +
1
(1 )
( ) ( )
+
+
. (4 72)
. (4 73)
. (4 74)
Oxgeno Disuelto
DBO
Fitoplancton
Periphyton
Detritus (C, N, P)
Amonio
Nitrato
Ortofosfato
Salinidad
Slidos
Pgina | 96
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Sedimentos
o CALOR
Temperatura
Salinidad
Coliformes
Conservativos 1 y 2
1 Txico
Limo y Arcilla
Arenas
Slidos orgnicos
o TXICOS ORGNICOS
Qumico 1
Qumico 2
Qumico 3
Limo y Arcilla
Arenas
Slidos orgnicos
o MERCURIO
Elemental, Hg0
Divalente, HgII
Metil, MeHg
Limo y Arcilla
Arena
Slidos orgnicos
Pgina | 97
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 98
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Posicin Geogrfica
Estacin
Longitud
Latitud
Tipo de muestra
PC1
68 27.791' W
16 33.638'
W, IS, L, S
PI2
68 28.138' W
16 32.304'
W, IS, L, S
PC3
68 29.133' W
16 32.159'
W, IS, L, S
PC4
68 35.799' W
16 25.373'
W, IS, L, S
PC5
68 37.801' W
16 23.579'
W, IS, L, S
Tabla5-1. Estaciones de monitoreo: ubicacin y tipo de muestra tomada (W = Agua, IS = in situ = Laboratorio, S = Simple).
Fuente: Apuntes de clases de Monitoreo Ambiental (IMA-901).
Parmetro
Equipo
Longitud
GPS
Latitud
Ancho del ro
Profundidad
media del ro
Frecuencia
Tcnica
Medicin
in situ
GPS
in situ
Cinta mtrica
in situ
Cinta mtrica
in situ
Anual
Aforo con
Caudal
Molinete/Sal y Conductimetro
molinete/aforo
in situ
qumico
Aforo con
Velocidad
molinete/aforo
in situ
qumico
Temperatura
Probe
in situ
Pgina | 99
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Conductividad
Electrodo
in situ
pH
Electrodo
in situ
in situ/Mtodo
in situ /
Wrinkler
Laboratorio
Oxigeno
Electrodo de O 2 , Hach
Dissolved Oxygen
Disuelto
Probe/Bureta, matraz
Turbidez
DBO 5
Incubadora
DQO
Nitrgeno-NH 4
Nitrgeno-NO 3
Nitrgeno-NO 2
Recirculacin/Espectrofotmetro
portable DR 2800
Espectrofotmetro portable DR
2800
Espectrofotmetro portable DR
2800
Espectrofotmetro portable DR
2800
Nitrgeno
Digestor de Semi-micro
Kjeldahl
Kjeldahl
in situ /
laboratorio
Incubacin 5
das 20C
Laboratorio
Reflujo
2
abierto/ Hach
Laboratorio
8000
2
Laboratorio
Laboratorio
Laboratorio
Semi-micro
Kjeldahl
Laboratorio
Tabla5-2. Parmetros fsicos y qumicos medidos en laboratorio, equipos utilizados, periodicidad de la medicin y protocolo de medida seguido.
Fuente: Apuntes de clases de Monitoreo Ambiental (IMA-901).
En la Tabla 5-3 se presentan los parmetros (in situ y de laboratorio) determinados en cada
campaa de muestreo.
Pgina | 100
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
1 Campaa
2 Campaa
3 Campaa
4 Campaa
5 Campaa
6 Campaa
7 Campaa
8 Campaa
9 Campaa
Fecha de realizacin
Abril/2007
Noviembre/2007
Abril/2008
Abril/2009
Octubre/2009
Junio/2010
Agosto/2010
Noviembre/2010
Marzo/2011
poca del ao
Lluviosa
Seca
Lluviosa
Lluviosa
Seca
Seca
Seca
Seca
Lluviosa
N Puntos muestreo
Longitud
Si
No
No
No
No
Si
Si
Si
Si
Latitud
Si
No
No
No
No
Si
Si
Si
Si
Ancho del ro
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Caudal
Si
Si
Si
Si
Si
Si
Si
Si
Si
Velocidad
No
No
No
No
No
Si
Si
Si
Si
Temperatura
Si
Si
Si
Si
Si
Si
Si
Si
Si
Conductividad
Si
Si
Si
Si
Si
Si
Si
Si
Si
pH
Si
Si
Si
Si
Si
Si
Si
Si
Si
Oxigeno Disuelto
Si
Si
Si
Si
Si
Si
Si
Si
Si
Turbidez
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Si
Parmetros in situ
ro
Parmetros en
laboratorio
DBO 5
Si
Si
Si
Si
Si
Si
Si
Si
Si
Nitrgeno-NH
Si
No
Si
No
No
Si
No
Si
No
Nitrgeno-NO 3
No
No
Si
Si
No
Si
Si
Si
No
Nitrgeno-NO 2
No
No
No
No
No
No
No
Si
No
Nitrgeno Kjeldal
No
Si
Si
Si
Si
Si
Si
Si
No
DQO
Tabla5-3. Detalle de los parmetros medidos in situ y en laboratorio en las campaas realizadas.
Fuente: Elaboracin Propia.
Pgina | 101
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Idealmente, todas las estaciones de monitoreo, deberan poseer caractersticas similares con
respecto al punto de recoleccin de muestra, ya sea que el mismo se encuentre en la orilla o en
la mitad del ro. Durante el trabajo de seleccin de las mismas se prioriz que se encuentren a
la mitad del ro y a una profundidad adecuada, siendo stos lugares ptimos para el muestreo
de todos los parmetros y en donde las muestras son ms representativas del mezclado del ro.
Tambin se tomo en cuenta el tiempo de paso en las ultimas 2 campaas que sirvieron para
validar el modelo SIMOD, ya que la principal caracterstica de este modelo es aplicar el
muestreo intensivo para poder obtener un modelado preciso.
El Grfico 19 muestra las diferentes maniobras de muestreo desarrolladas, en una estacin
situada en el medio del ro y en otra situada en un puente respectivamente (debido al riesgo de
entrar al ro por la fangosidad). Ver Anexos A y C
Grafico5-1. Maniobra tpica de toma de muestra en dos estaciones (PC1 y PC5) situada en el medio del ro y sobre un puente del Ro Colorado y Ro Katari.
Fuente: Elaboracin Propia.
Pgina | 102
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
OD (mg/l)
8
7
OD (mg/l)
6
5
4
3
2
1
0
ABRIL 2007
NOVIEMBRE 2007
OD - PC1
ABRIL 2008
ABRIL 2009
OD - PI2
OCTUBRE 2009
OD - PC3
JUNIO 2010
AGOSTO 2010
OD - PC4
NOVIEMBRE 2010
MARZO 2011
OD - PC5
Grafico6-1. Grafico del Oxigeno Disuelto (OD) de los diferentes puntos en las distintas campaas realizadas.
Fuente: Elaboracin Propia.
Los perfiles longitudinales de la concentracin de OD muestran que los valores ms altos son
obtenidos del Ro Pallina siento este la principal fuente de contaminacin del Ro Katari.
La concentracin de OD en el agua del Ro Katari, est vinculada al proceso de reaireacin
natural (fuente) y oxidacin de la materia orgnica biodegradable (sumidero). En el ro se
tiene apropiadas condiciones fsicas para el proceso de reaireacin natural por el contacto con
el aire atmosfrico, de esta manera se favorece la actividad microbiolgica aerbica.
En el reglamento en Materia de Contaminacin Hdrica (RMCH) que pertenece a la Ley del
Medio Ambiente N1333, el limite permisible para la concentracin de OD es 4[mg/l], en los
datos obtenidos podemos observar que los valores que no cumplen son lo de la descarga del
Ro Pallina y los valores despus del punto de confluencia.
Michael Russel Osina Torrez
Pgina | 103
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Salinidad
1.4
1.2
S (mg/l)
1.0
0.8
0.6
0.4
0.2
0.0
ABRIL 2007
NOVIEMBRE 2007
ABRIL 2008
PC1
ABRIL 2009
PI2
OCTUBRE 2009
PC3
JUNIO 2010
PC4
AGOSTO 2010
NOVIEMBRE 2010
MARZO 2011
PC5
Grafico6-2. Grafico de la Salinidad (S) de los diferentes puntos en las distintas campaas realizadas.
Fuente: Elaboracin Propia.
Las sales disueltas en agua reducen los espacios intermoleculares disponibles para la
disolucin de O 2 . La salinidad se calcula a partir de la ecuacin de Weiss (1970), utilizando el
valor de CE. La salinidad sirve tambin para corregir la Concentracin de Saturacin del
Oxigeno Disuelto (Cs), a altos valores de S menor es el valor de Cs.
= 5.572 104 [ ] + 2.02 109 [ ]2
Pgina | 104
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Cs (mg/l)
7.0
6.5
6.0
5.5
5.0
0
5.33
5.62
28.67
34.43
Distancia (Km)
1C
2C
3C
4C
5C
6C
7C
8C
9C
Grafico6-3. Grafico de la Concentracin de Saturacin del OD (Cs) de los diferentes puntos en las distintas campaas realizadas (Ec. 3-12 y 3-14).
Fuente: Elaboracin Propia.
Cs (mg/l)
7.0
6.5
6.0
5.5
5.0
0
5.33
5.62
28.67
34.43
Distancia (Km)
1C
2C
3C
4C
5C
6C
7C
8C
9C
Grafico6-4. Grafico de la Concentracin de Saturacin del OD (Cs) de los diferentes puntos en las distintas campaas realizadas (Ec. 4-56 y 4-57).
Fuente: Elaboracin Propia.
Pgina | 105
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
La Concentracin de saturacin (Cs) de OD, se calcula utilizando las ecuaciones 3-12, 3-14
donde toma en cuenta la correccin de la salinidad, temperatura y presin, y las ecuaciones 456, 4-57 que son las ecuaciones con las cuales calcula el modelo SIMOD donde toma en
cuenta la temperatura y la presin.
Los valores calculados muestran valores mximos de 7.063 y 7.049 (mg/l) y mnimos de
5.198 y 5.185 (mg/l).
Como se puede observar en las graficas la importancia de la salinidad es poco significativa ya
que no existe una gran variacin entre los valores obtenidos. Prcticamente son similares.
Las descargas de aguas residuales domesticas, comerciales e industriales, influyen en los
niveles de temperatura del agua en los ros, porque despus de la descarga y punto de mezcla,
provocan su incremento disminuyendo el valor de la Cs.
El aumentar el valor de S, se puede observar que su efecto es la disminucin de los valores de
Cs de OD. Es importante considerar esta correccin cuando se tengan valores muy elevados
de CE.
A continuacin se presenta algunos ejemplos de la variacin de Cs a valores altos de
Conductividad especfica.
Ejemplo
Conductividad
Especfica [S/cm]
Salinidad [g/Kg]
Descenso en la
Cs de OD [mg/l]
1300
0.73
0.05
6000
3.42
0.23
20000
11.95
0.77
30000
18.53
1.18
40000
25.52
1.58
55000
36.76
2.21
1
2
3
4
5
6
Pgina | 106
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
DBO5 (mg/l)
100
90
80
DBO (mg/l)
70
60
50
40
30
20
10
0
ABRIL 2007
NOVIEMBRE 2007
ABRIL 2008
PC1
ABRIL 2009
OCTUBRE 2009
PI2
PC3
JUNIO 2010
PC4
AGOSTO 2010
NOVIEMBRE 2010
MARZO 2011
PC5
Grafico6-5. Grafico de la Demanda Bioqumica de Oxigeno (DBO) de los diferentes puntos en las distintas campaas realizadas.
Fuente: Elaboracin Propia.
Los perfiles de la Demanda Bioqumica de oxigeno a los 5 das (DBO 5 ), muestran que los
valores ms altos son de las aguas del Ro Pallina debido a las descargas que son vertidas a
este ro en su trayecto y las ms bajas son del Ro Colorado.
La concentracin de la DBO 5 se relaciona con la materia orgnica biodegradable soluble y
particulada, que se encuentra en las descargas de aguas residuales domesticas, comerciales e
industriales y su valor puede incrementar por la influencias de los residuos slidos que son
arrojados indiscriminadamente el los cuerpos receptores, que tienen su origen en la
actividades domsticas, comerciales e industriales.
Como se puede observar en la grafica los valores que se encuentran en las campaas
realizadas en poca de lluvia (abril y marzo) son relativamente ms bajos comparados a los
dems valores que se encuentran en poca seca, esto es debido a la dilucin que existe con el
agua de lluvia. Por este motivo es recomendable y ms preciso modelar con datos de poca
Pgina | 107
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
seca. En poca seca se presenta el escenario, ms desfavorable debido a que no existe dilucin
de los contaminantes.
Tro
25
20
T(C)
15
10
ABRIL 2007
NOVIEMBRE 2007
ABRIL 2008
PC1
ABRIL 2009
PI2
OCTUBRE 2009
PC3
JUNIO 2010
PC4
AGOSTO 2010
NOVIEMBRE 2010
MARZO 2011
PC5
Grafico6-6. Grafico de la temperatura del ro en grados centgrados de los diferentes puntos en las distintas campaas realizadas.
Fuente: Elaboracin Propia.
Pgina | 108
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
CE
2.5
CE (mS/cm)
1.5
0.5
ABRIL 2007
NOVIEMBRE 2007
ABRIL 2008
PC1
ABRIL 2009
PI2
OCTUBRE 2009
PC3
JUNIO 2010
PC4
AGOSTO 2010
NOVIEMBRE 2010
MARZO 2011
PC5
Grafico6-7. Grafico de la Conductividad elctrica de los diferentes puntos en las distintas campaas realizadas.
Fuente: Elaboracin Propia.
Pgina | 109
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
pH
11
10
pH
ABRIL 2007
NOVIEMBRE 2007
ABRIL 2008
PC1
ABRIL 2009
PI2
OCTUBRE 2009
PC3
JUNIO 2010
PC4
AGOSTO 2010
NOVIEMBRE 2010
MARZO 2011
PC5
Grafico6-8. Grafico del Potencial de hidrogeno (pH) de los diferentes puntos en las distintas campaas realizadas.
Fuente: Elaboracin Propia.
Pgina | 110
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 111
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tramo N2
Tramo N3
Pgina | 112
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tramo N4
Perfil Longitudinal
3840
Altitud (m.s.n.m.)
3838
3836
3834
3832
3830
3828
3826
3824
0
10
15
20
25
30
35
40
Distancia (Km)
Pgina | 113
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Cada tramo presenta condiciones ambientales propias. Toda la trayectoria del ro tiene una
longitud de 34430 [m] y diferencias de nivel de 13 [m] respectivamente. En el grafico 7-5 se
aprecia la diferencia de nivel para cada tramo.
7.1.1. CALIBRACIN
Para la Calibracin de los modelos se utilizo el siguiente algoritmo de decisin:
K1
K1
Resultados
Resultados
Rio
Katari
Rio Katari
Mtodos
Mtodos
K2
K2
No
Calibrado?
Si
Metodoi,j
Metodoi,j
monitoreo de la calidad del agua en ros. Para su aplicacin, utilizando el concepto de enfoque
sistmico, previamente es preciso definir la configuracin del sistema ambiental (aspectos
hdricos y de borde, el ncleo, establecer el seccionamiento ambiental del ro y la interseccin
de todos los componentes del indicado sistema.
Pgina | 114
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 115
5.40
24.00
19.59
130.00
11.1
14.9
14.5
12.7
3826
130.00
5.40
Temp. (C)
13.0
DBO 5 (mg/l)
1.15 3.750
1.380
1.15 3.745
1.380
0.55
16 23.574'S,
68 37.796'O
3828
0.45
16 25.610'S,
68 36.445'O
3835
16 32.159'S,
68 29.133'O
0.35
3836
16 32.311'S,
68 29.161'O
1.119
34+43 PC5
28+67 PC4
PC3
P2
5+33
5+62
PC1
5+33
3836
0.319
5.76
11.1
OD (mg/l)
7.30 0.574
Q (m^3/s)
0.319
23.05
3.15 3.745
Vm (m/s)
0.29
0.35
Prof. (m)
0.001
2.40 3.171
16 33.638'S,
68 27.791'O
3839
Coordenadas
(, min dec.)
0.30
Elev. (m.s.n.m.)
5.33
16 32.311'S,
68 29.161'O
Distancia (Km)
0.30
P0
Progresiva
7.30 0.574
Punto
0+00
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Grafico7-7. Seccionamiento Ambiental del Ro Katari (Datos 1 campaa) para el modelo SIMOD.
Fuente: Elaboracin Propia.
Pgina | 116
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Para calibrara el modelo SIMOD se aplica el algoritmo de decisin (grafico 7-6), tomando en
cuenta las caractersticas mencionadas en el Captulo 4.
Introduciendo la informacin del Anexo B se realizaron 9 combinaciones para las campaas
Empleando el modelo SIMOD, en cada corrida se modificaron las constantes de
desoxigenacin y reaireacin mediante las opciones que ofrece el modelo, como se muestra
en la Tabla 7-1.
N de Corrida
Combinacin
Mtodos K 1 (1/da)
Mtodos K 2 (1/da)
11
12
13
21
Bansal (2)
22
Bansal (2)
23
Bansal (2)
31
Bosko (3)
32
Bosko (3)
33
Bosko (3)
Unidad
Rango de valores
Reaireacin
dia-1
0.00 100.00
Desoxigenacin
dia-1
0.00 2.00
Pgina | 117
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
El rango de valores obtenido por el modelo SIMOD v1.01 en el Ro Katari para las diferentes
campanas se muestran a continuacin:
Constante
Unidad
Rango de valores
Reaireacin
dia-1
0.65 30.00
Desoxigenacin
dia-1
0.12 1.95
A continuacin, a manera de ejemplo se muestran dos de los nueve resultados obtenidos por
la modelacin con SIMOD que corresponden a la calibracin de dicho modelo. Los resultados
mencionados estn expresados en dos graficas de OD y DBO y una tabla. Estas corridas se
calibraron con el mtodo de Bansal para la determinacin de K1 y OConnor et al para la
determinacin de K 2 . Las restantes corridas se encuentran en el Anexo D.
Pgina | 118
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
POCA DE LLUVIA
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O'Connor et al)
8
7
5
DBO (mg/L)
OD (mg/L)
6
4
3
2
1
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
O.D. Observado (mg/L)
Distancia Km
SIMULACION DBO5
(Kd: M. Bansal y K2: O'Connor et al)
40
36
32
28
24
20
16
12
8
4
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
Abril
Distancia
Segmento
Acumulada(Km)
2008
Temperatura
O.D.
DBO
O.D.
del Agua
SIMOD
SIMOD
Saturacin
(C)
(mg/l)
(mg/L)
(mg/l)
O.D.
K1
K2
Mnimo
(1/da)
(1/da)
Permitido
(mg/L)
O.D.
DBO 5
Observado
Observada
(mg/L)
(mg/L)
5.45
16.00
0.00
22.20
5.45
16.00
5.45
0.28
26.16
4.00
1.57
22.20
5.24
15.49
5.45
0.28
26.16
4.00
3.13
22.20
5.23
14.99
5.45
0.28
26.16
4.00
4.70
22.20
5.24
14.51
5.45
0.28
26.16
4.00
5.33
22.20
5.24
14.32
5.45
0.28
26.16
4.00
5.33
20.87
5.59
34.07
5.59
0.24
9.15
4.00
5.33
20.87
5.59
34.07
5.59
0.24
9.15
4.00
5.33
20.87
5.59
34.07
5.59
0.24
13.06
4.00
5.62
20.87
5.47
33.98
5.59
0.24
13.06
4.00
5.62
20.87
5.47
34.52
5.60
0.25
8.04
4.00
15.62
20.87
4.30
32.20
5.60
0.25
8.04
4.00
25.62
20.87
4.26
30.03
5.60
0.25
8.04
4.00
4.5
25.00
34.43
20.87
4.33
28.25
5.60
0.25
8.04
4.00
4.5
24.00
Pgina | 119
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
POCA SECA
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O'Connor et al)
8
7
6
OD (mg/L)
DBO (mg/L)
4
3
2
1
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia Km
SIMULACION DBO5
(Kd: M. Bansal y K2: O'Connor et al)
60
55
50
45
40
35
30
25
20
15
10
5
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia enDBO
Km
SIMOD (mg/L)
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
N
Segmento
Distancia
Acumulada(Km)
Temperatura
del Agua
(C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D.
Mnimo
Permitido
(mg/L)
1
1
1
1
1
1
1
2
2
3
3
4
4
4
4
0.00
1.00
2.00
3.00
4.00
5.00
5.33
5.33
5.33
5.33
5.62
5.62
15.62
25.62
34.43
17.00
17.00
17.00
17.00
17.00
17.00
17.00
20.08
20.08
20.08
20.08
20.08
20.08
20.08
20.08
5.00
5.50
5.76
5.89
5.97
6.00
6.01
1.34
1.34
1.34
1.42
1.42
2.31
3.02
3.53
12.00
6.29
3.30
1.73
0.91
0.48
0.39
37.71
37.69
53.99
53.00
51.96
40.99
32.33
26.24
6.05
6.05
6.05
6.05
6.05
6.05
6.05
5.68
5.68
5.68
5.68
5.68
5.68
5.68
5.68
0.28
0.28
0.28
0.28
0.28
0.28
0.28
0.23
0.23
0.80
0.80
0.66
0.66
0.66
0.66
4.52
4.52
4.52
4.52
4.52
4.52
4.52
3.16
3.16
10.92
10.92
8.94
8.94
8.94
8.94
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
O.D.
Observado
(mg/L)
DBO 5
Observada
(mg/L)
12.00
3.98
4.15
39.00
38.00
Pgina | 120
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
7.1.1.2.MODELO RIOSep
La introduccin de datos al modelo RIOSep se realiza por tramos, y de acuerdo a los
requerimientos del modelo solo es posible considerar una descarga por tramo.
7.1.1.2.1. SEGMENTACIN AMBIENTAL DEL RO
Para la recopilacin de los datos necesarios para la calibracin del modelo RIOSep se realiza
el siguiente esquema (grafico 7-12) en el que se presenta el seccionamiento del Ro Katari.
ODobs = v
DBOobs = v
Rio Katari
Rio Katari
ODobs = v
DBOobs = v
Rio Katari
Descarga
Rio Pallina
Condiciones iniciales
Rio Colorado
ODobs = v
DBOobs = v
Qo = v
Vo = v
Ho = v
ODo = v
DBOo = v
To = v
Eo = v
QD = v
VD = v
HD = v
ODD = v
DBOD = v
TD = v
ED = v
Los datos que ingresan al modelo son los datos iniciales y los del afluente (Ro Pallina)
considerando que los caudales son constantes a lo largo de la trayectoria, ya que no existen ni
extracciones ni otras descargas hasta llegar a la Baha de Cohana. Para los datos observados
se considera los valores obtenidos de DBO 5 y OD.
Para emplear el modelo RIOSep se utilizo informacin que se encuentra en el Anexo B.
Se realizaron diferentes corridas con cada campaa utilizada para calibrar el modelo. Se
emplearon las diferentes opciones que el modelo ofrece para el valor de la constante de
reaireacin como se muestra en la Tabla 7-6. Sin embargo el modelo RIOSep no ofrece
mtodos de clculo de la constante de desoxigenacin, solamente requiere la introduccin del
valor.
Michael Russel Osina Torrez
Pgina | 121
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
(7 1)
Donde:
Combinacin
Mtodos K 1 (1/da)
Mtodos K 2 (1/da)
11
O'Connor y Dobbins
12
13
14
Bennett y Rathbun
15
Tsivoglou (5)
21
O'Connor y Dobbins
22
23
24
Bennett y Rathbun
10
25
Tsivoglou (5)
Pgina | 122
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tra
mo
1
Campa
a
2
Campa
a
3
Campa
a
4
Campa
a
5
Campa
a
6
Campa
a
No
cumple
7
Campa
a
8
Campa
a
9
Campa
a
13.51
5.93
22.08
15.41
5.93
0.71
5.93
78.60
14.56
22.15
8.15
13.36
12.78
5.52
15.40
7.41
12.53
12.45
0.90
16.00
12.82
0.90
67.67
0.07
0.90
No
cumple
23.71
27.19
7.49
15.78
10.26
4.12
18.03
3.42
22.72
0.48
6.69
No
cumple
Churchill et al.
1
23.03
6.69
40.62
26.68
6.69
No
cumple
27.16
43.67
12.44
23.46
21.15
7.52
27.86
10.13
22.94
0.94
6.68
61.59
13.72
7.76
10.62
Owens et al.
Bennett y
Rathbun
11.45
6.68
18.91
13.16
6.68
No
cumple
12.43
19.29
7.73
11.97
12.10
10.28
2.75
0.05
1.34
2.07
0.05
0.72
0.03
0.05
83.98
6.56
2.94
1.93
3.08
1.22
12.65
2.90
0.37
2.64
Tsivoglou
Tabla7-7. Constantes de reaireacin para cada tramo de cada campaa por cinco diferentes mtodos.
Fuente: Elaboracin propia.
2 = 0.9862
(7 2)
K1 vs. K2
1.0
0.8
0.6
0.4
0.2
0.0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
K1 vs. K2
Pgina | 123
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Ajuste
de K 1
y K2 SIMO
D
Tram
o
1
Campaa
2
Campaa
3
Campaa
4
Campaa
5
Campaa
6
Campaa
No
cumple
7
Campaa
8
Campaa
0.65
0.32
1.00
0.74
0.32
0.70
1.03
0.42
0.65
9
Campaa
No
cumple
0.09
0.32
0.62
0.30
0.74
0.38
0.61
0.06
0.10
No
cumple
0.61
0.10
0.76
0.62
0.10
No
cumple
1.10
1.25
0.39
0.75
0.51
0.24
0.85
0.21
1.00
0.08
0.35
No
cumple
1.07
0.35
1.85
1.23
0.35
No
cumple
1.25
1.98
0.61
1.09
0.99
0.39
1.28
0.50
1.00
0.10
0.35
2.77
0.56
0.35
0.89
0.64
0.35
No
cumple
0.60
0.91
0.40
0.58
0.59
0.51
0.66
0.40
0.53
0.18
0.06
0.12
0.15
0.06
0.09
0.06
0.06
No
cumple
0.35
0.19
0.14
0.19
0.11
0.61
0.19
0.07
0.17
Tabla7-8. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la ecuacin ajustada.
Fuente: Elaboracin propia.
Calculando los valores de K1 con la ecuacin propuesta por Chapra (1997). (10)
0.434
1 = 0.3
8
Constante de desoxigenacin (K 1 )
Mtodo
Tramo
1 Campaa
2 Campaa
3 Campaa
4 Campaa
5 Campaa
6 Campaa
7 Campaa
8 Campaa
9 Campaa
Chapra
(1997)
1.25
1.79
1.60
1.35
1.79
2.71
1.05
1.79
1.27
1.05
1.33
1.06
1.14
1.29
1.00
1.20
1.31
0.96
Tabla7-9. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la ecuacin Chapra.
Fuente: Elaboracin propia.
A continuacin se muestra dos de las mejores corridas calibradas por el modelo utilizando la
combinacin 11 de la tabla 7-6. (Ver Anexo F para el clculo de la K 1 en laboratorio con
datos obtenidos en la octava campaa).
Pgina | 124
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
POCA DE LLUVIA
SIMULACION DE O.D.
(Kd: Ec. ajustada y K2: O'Connor et al)
SIMULACION DBO5
(Kd: Ec. ajustada y K2: O'Connor et al)
40
35
30
25
DBO (mg/L)
OD (mg/L)
6
5
20
15
10
2
0
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
OD RIOSep (mg/l)
t (h)
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
X (Km)
DBO 5 (mg/l)
DOD (mg/l)
OD (mg/l)
OD saturacin
16
-0.5786
5.45
4.8714
0.5
0.8956
15.3339
0.2014
4.6701
4.8715
1.7912
14.6955
0.5082
4.3633
4.8715
1.5
2.6868
14.0837
0.6184
4.253
4.8714
3.5824
13.4974
0.6475
4.224
4.8715
2.5
4.478
12.9355
0.6434
4.2281
4.8715
5.3736
12.397
0.6261
4.2453
4.8714
3.5
6.2692
11.8809
0.604
4.2674
4.8714
7.1648
11.3863
0.5805
4.2909
4.8714
4.5
8.0603
10.9122
0.5571
4.3144
4.8715
8.9559
10.458
0.5342
4.3373
4.8715
5.33
9.547
10.1685
0.5195
4.352
4.8715
5.33
33.3827
5.0387
5.0387
5.83
0.4044
33.1407
0.2262
4.8124
5.0386
6.33
0.8089
32.9004
0.4221
4.6166
5.0387
6.83
1.2133
32.6619
0.5914
4.4473
5.0387
7.33
1.6177
32.4252
0.7376
4.301
5.0386
7.83
2.0222
32.1901
0.8637
4.175
5.0387
8.33
2.4266
31.9568
0.9721
4.0666
5.0387
8.83
2.831
31.7251
1.0651
3.9736
5.0387
9.33
3.2354
31.4952
1.1448
3.8939
5.0387
9.83
3.6399
31.2669
1.2128
3.8259
5.0387
10.33
4.0443
31.0402
1.2706
3.7681
5.0387
10.83
4.4487
30.8152
1.3195
3.7192
5.0387
11.33
4.8532
30.5918
1.3606
3.678
5.0386
11.83
5.2576
30.3701
1.3951
3.6436
5.0387
12.33
5.662
30.1499
1.4236
3.6151
5.0387
12.83
6.0665
29.9314
1.447
3.5916
5.0386
13.33
6.4709
29.7144
1.466
3.5727
5.0387
13.83
6.8753
29.499
1.4811
3.5576
5.0387
Pgina | 125
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
14.33
7.2797
29.2852
1.4928
3.5459
5.0387
14.83
7.6842
29.0729
1.5016
3.5371
5.0387
15.33
8.0886
28.8621
1.5078
3.5309
5.0387
15.83
8.493
28.6529
1.5118
3.5269
5.0387
16.33
8.8975
28.4452
1.5139
3.5248
5.0387
16.83
9.3019
28.239
1.5143
3.5244
5.0387
17.33
9.7063
28.0343
1.5132
3.5254
5.0386
17.83
10.1108
27.8311
1.5109
3.5277
5.0386
18.33
10.5152
27.6294
1.5076
3.5311
5.0387
18.83
10.9196
27.4291
1.5032
3.5355
5.0387
19.33
11.3241
27.2302
1.4981
3.5406
5.0387
19.83
11.7285
27.0329
1.4923
3.5464
5.0387
20.33
12.1329
26.8369
1.4859
3.5528
5.0387
20.83
12.5373
26.6424
1.4789
3.5598
5.0387
21.33
12.9418
26.4492
1.4716
3.5671
5.0387
21.83
13.3462
26.2575
1.4638
3.5749
5.0387
22.33
13.7506
26.0672
1.4558
3.5829
5.0387
22.83
14.1551
25.8782
1.4474
3.5913
5.0387
23.33
14.5595
25.6906
1.4389
3.5998
5.0387
23.83
14.9639
25.5044
1.4301
3.6085
5.0386
24.33
15.3684
25.3195
1.4213
3.6174
5.0387
24.83
15.7728
25.136
1.4122
3.6264
5.0386
25.33
16.1772
24.9538
1.4031
3.6355
5.0386
25.83
16.5816
24.7729
1.394
3.6447
5.0387
26.33
16.9861
24.5933
1.3847
3.654
5.0387
26.83
17.3905
24.415
1.3754
3.6633
5.0387
27.33
17.7949
24.2381
1.3661
3.6726
5.0387
27.83
18.1994
24.0624
1.3568
3.6819
5.0387
28.33
18.6038
23.8879
1.3474
3.6912
5.0386
28.83
19.0082
23.7148
1.3381
3.7006
5.0387
29.33
19.4127
23.5429
1.3288
3.7099
5.0387
29.83
19.8171
23.3722
1.3195
3.7192
5.0387
30.33
20.2215
23.2028
1.3102
3.7285
5.0387
30.83
20.626
23.0346
1.301
3.7377
5.0387
31.33
21.0304
22.8676
1.2918
3.7469
5.0387
31.83
21.4348
22.7018
1.2826
3.7561
5.0387
32.33
21.8392
22.5373
1.2735
3.7652
5.0387
32.83
22.2437
22.3739
1.2644
3.7743
5.0387
33.33
22.6481
22.2117
1.2553
3.7833
5.0386
33.83
23.0525
22.0507
1.2464
3.7923
5.0387
34.33
23.457
21.8909
1.2374
3.8013
5.0387
34.43
23.5379
21.859
1.2356
3.8031
5.0387
Pgina | 126
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
POCA SECA
SIMULACION DE O.D.
(Kd: Ec. ajustada y K2: O'Connor et al)
SIMULACION DBO5
(Kd: Ec. ajustada y K2: O'Connor et al)
55
50
45
40
30
35
DBO (mg/L)
OD (mg/L)
2
1
25
20
15
10
5
0
0
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
OD RIOSep (mg/l)
X (Km)
t (h)
DBO 5 (mg/l)
DOD (mg/l)
OD (mg/l)
OD saturacin
12
0.4054
5.4054
0.5
25.2778
8.8647
0.4592
4.9462
5.4054
50.5556
6.5485
0.3396
5.0658
5.4054
1.5
75.8333
4.8376
0.2508
5.1546
5.4054
101.1111
3.5736
0.1853
5.2201
5.4054
2.5
126.3889
2.6399
0.1369
5.2685
5.4054
151.6667
1.9502
0.1011
5.3043
5.4054
3.5
176.9444
1.4406
0.0747
5.3307
5.4054
202.2222
1.0642
0.0552
5.3502
5.4054
4.5
227.5
0.7862
0.0408
5.3646
5.4054
252.7778
0.5808
0.0301
5.3753
5.4054
5.33
269.4611
0.4756
0.0247
5.3807
5.4054
5.33
37.7114
3.7573
1.3314
5.0887
5.83
1.7083
36.6742
2.9804
2.1083
5.0887
6.33
3.4167
35.6655
2.5074
2.5813
5.0887
6.83
5.125
34.6846
2.2114
2.8773
5.0887
7.33
6.8333
33.7307
2.0188
3.0699
5.0887
7.83
8.5417
32.803
1.8867
3.202
5.0887
8.33
10.25
31.9008
1.7904
3.2983
5.0887
8.83
11.9583
31.0235
1.7154
3.3733
5.0887
9.33
13.6667
30.1702
1.6532
3.4355
5.0887
9.83
15.375
29.3405
1.5991
3.4896
5.0887
10.33
17.0833
28.5335
1.55
3.5387
5.0887
10.83
18.7917
27.7488
1.5045
3.5842
5.0887
11.33
20.5
26.9856
1.4614
3.6273
5.0887
11.83
22.2083
26.2434
1.4202
3.6685
5.0887
12.33
23.9167
25.5216
1.3806
3.7081
5.0887
12.83
25.625
24.8197
1.3423
3.7464
5.0887
Pgina | 127
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
13.33
27.3333
24.1371
1.3052
3.7835
5.0887
13.83
29.0417
23.4733
1.2692
3.8195
5.0887
14.33
30.75
22.8277
1.2342
3.8545
5.0887
14.83
32.4583
22.1998
1.2002
3.8885
5.0887
15.33
34.1667
21.5893
1.1672
3.9215
5.0887
15.83
35.875
20.9955
1.1351
3.9536
5.0887
16.33
37.5833
20.4181
1.1038
3.9849
5.0887
16.83
39.2917
19.8565
1.0735
4.0152
5.0887
17.33
41
19.3104
1.044
4.0447
5.0887
17.83
42.7083
18.7793
1.0152
4.0735
5.0887
18.33
44.4167
18.2628
0.9873
4.1014
5.0887
18.83
46.125
17.7605
0.9602
4.1285
5.0887
19.33
47.8333
17.2721
0.9338
4.1549
5.0887
19.83
49.5417
16.797
0.9081
4.1806
5.0887
20.33
51.25
16.3351
0.8831
4.2056
5.0887
20.83
52.9583
15.8858
0.8588
4.2299
5.0887
21.33
54.6667
15.4489
0.8352
4.2535
5.0887
21.83
56.375
15.024
0.8122
4.2765
5.0887
22.33
58.0833
14.6108
0.7899
4.2988
5.0887
22.83
59.7917
14.209
0.7682
4.3205
5.0887
23.33
61.5
13.8182
0.747
4.3417
5.0887
23.83
63.2083
13.4381
0.7265
4.3622
5.0887
24.33
64.9167
13.0686
0.7065
4.3822
5.0887
24.83
66.625
12.7091
0.6871
4.4016
5.0887
25.33
68.3333
12.3596
0.6682
4.4205
5.0887
25.83
70.0417
12.0197
0.6498
4.4389
5.0887
26.33
71.75
11.6891
0.6319
4.4568
5.0887
26.83
73.4583
11.3676
0.6146
4.4741
5.0887
27.33
75.1667
11.055
0.5977
4.491
5.0887
27.83
76.875
10.7509
0.5812
4.5075
5.0887
28.33
78.5833
10.4552
0.5652
4.5235
5.0887
28.83
80.2917
10.1677
0.5497
4.539
5.0887
29.33
82
9.8881
0.5346
4.5541
5.0887
29.83
83.7083
9.6161
0.5199
4.5688
5.0887
30.33
85.4167
9.3516
0.5056
4.5831
5.0887
30.83
87.125
9.0944
0.4917
4.597
5.0887
31.33
88.8333
8.8443
0.4781
4.6106
5.0887
31.83
90.5417
8.6011
0.465
4.6237
5.0887
32.33
92.25
8.3645
0.4522
4.6365
5.0887
32.83
93.9583
8.1345
0.4398
4.6489
5.0887
33.33
95.6667
7.9107
0.4277
4.661
5.0887
33.83
97.375
7.6932
0.4159
4.6728
5.0887
34.33
99.0833
7.4816
0.4045
4.6842
5.0887
34.43
99.425
7.44
0.4022
4.6865
5.0887
Pgina | 128
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 129
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
El modelo requiere que se introduzca datos de frontera los cuales se muestran en la grafica 718. En los cuales se introduce las concentraciones de contaminante que ingresan y los datos
hidrodinmicos. Tambin se necesita introducir los datos hidrodinmicos de cada segmento,
en los cuales se asume que tienen las mismas caractersticas dentro de cada segmento.
Para emplear el modelo WASP 7.4 se utilizo informacin que se encuentra en el Anexo B.
Se realizaron diferentes corridas con cada campaa utilizada para calibrar el modelo. Se
emplearon las diferentes opciones que el modelo ofrece para el valor de la constante de
reaireacin como se muestra en la Tabla 7-12. Sin embargo el modelo WASP 7.4 no ofrece
mtodos de clculo de la constante de desoxigenacin, solamente requiere la introduccin del
valor. Tambin tiene la facilidad de introducir el valor de para la correccin por
temperatura de las constantes de desoxigenacin y reaireacin.
Para determinar el valor de la constante de desoxigenacin se utilizo la correlacin obtenida
de los valores promedio de K1 y K 2 por el modelo SIMOD (ver subtitulo 7.1.1.2.3.) y tambin
se calculo con la formula de Chapra [1997] (ver subtitulo 7.1.1.2.1.).
N de Corrida
Combinacin
Mtodos K 1 (1/da)
Mtodos K 2 (1/da)
11
O'Connor y Dobbins
12
13
14
Bennett y Rathbun
15
Tsivoglou (5)
21
O'Connor y Dobbins
22
23
24
Bennett y Rathbun
10
25
Tsivoglou (5)
Pgina | 130
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
O'Connor y
Dobbins
Churchill et
al.
Tramo
1
Campaa
2
Campaa
3
Campaa
4
Campaa
5
Campaa
6
Campaa
No
cumple
7
Campaa
8
Campaa
9
Campaa
13.51
6.22
22.08
15.41
6.22
0.71
5.93
78.60
20.09
19.97
11.66
17.78
16.68
5.52
12.22
4.79
12.53
15.30
26.58
13.39
14.93
23.11
5.95
17.41
9.95
16.20
11.32
25.06
7.18
12.22
12.29
5.95
17.41
7.57
16.20
12.45
0.99
16.00
12.82
0.99
67.67
0.07
0.90
No
cumple
32.46
21.18
10.52
21.18
11.35
4.12
10.25
1.66
22.72
26.12
36.87
13.76
18.26
22.68
4.76
24.61
5.49
47.69
18.67
34.53
6.87
14.60
11.21
4.76
24.61
4.05
47.69
0.48
6.69
No
cumple
23.03
7.13
40.62
26.68
7.13
No
cumple
40.16
37.75
19.19
33.27
28.63
7.52
20.14
5.73
22.94
29.03
55.38
23.11
26.98
44.66
8.32
33.17
14.84
33.41
20.03
51.51
10.72
21.07
20.49
8.32
33.17
10.59
33.41
1.18
8.34
61.59
Owens et al.
Bennett y
Rathbun
14.30
8.68
23.62
16.44
8.68
No
cumple
31.76
33.45
20.47
29.47
29.94
10.28
21.63
9.97
19.93
12.52
21.83
11.82
12.80
20.03
5.64
14.54
9.69
12.28
9.81
21.46
6.85
11.03
11.50
5.86
15.11
7.79
12.76
2.75
0.05
1.34
2.07
0.05
0.72
0.03
0.04
83.98
59.06
22.15
16.11
29.53
7.70
12.65
13.63
2.11
89.20
6.43
3.32
1.94
2.87
1.49
1.33
3.91
0.46
23.85
7.38
3.81
2.21
3.28
1.70
1.52
4.49
0.53
27.26
Tsivoglou
Tabla7-13. Constantes de reaireacin para cada tramo de cada campaa por cinco diferentes mtodos.
Fuente: Elaboracin propia.
Pgina | 131
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Constante de desoxigenacin (K 1 )
Mtodo
Tramo
1
Campaa
2
Campaa
3
Campaa
4
Campaa
5
Campaa
6
Campaa
7
Campaa
8
Campaa
9
Campaa
0.65
0.33
1.03
0.74
0.33
No cumple
0.09
0.32
No cumple
0.94
0.94
0.57
0.84
0.79
0.30
0.60
0.27
0.61
0.73
1.23
0.65
0.71
1.08
0.32
0.82
0.50
0.77
0.56
1.16
0.37
0.60
0.60
0.32
0.82
0.39
0.77
0.61
0.10
0.76
0.62
0.10
3.04
0.06
0.10
No cumple
1.49
0.99
0.52
0.99
0.56
0.24
0.51
0.13
1.06
1.21
1.68
0.66
0.86
1.06
0.27
1.14
0.30
2.16
0.88
1.58
0.36
0.70
0.55
0.27
1.14
0.24
2.16
1.07
0.37
1.85
1.23
0.37
No cumple
0.08
0.35
No cumple
1.83
1.72
0.90
1.52
1.32
0.39
0.94
0.31
1.07
1.34
2.50
1.08
1.25
2.02
0.42
1.52
0.71
1.53
0.94
2.32
0.53
0.99
0.96
0.42
1.52
0.52
1.53
0.69
0.44
1.10
0.78
0.44
No cumple
0.11
0.43
No cumple
1.46
1.53
0.96
1.35
1.38
0.51
1.01
0.50
0.93
0.61
1.02
0.58
0.62
0.94
0.31
0.70
0.48
0.60
0.49
1.00
0.36
0.54
0.56
0.32
0.72
0.40
0.62
0.18
0.06
0.12
0.15
0.06
0.09
0.06
0.06
No cumple
2.66
1.03
0.77
1.36
0.40
0.61
0.66
0.15
No cumple
0.34
0.20
0.14
0.18
0.12
0.12
0.23
0.08
1.11
0.38
0.23
0.16
0.20
0.13
0.13
0.26
0.08
1.26
Ec.
ajustada
Tabla7-14. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la ecuacin ajustada.
Fuente: Elaboracin propia.
Calculando los valores de K1 con la ecuacin propuesta por Chapra (1997). (10)
0.434
1 = 0.3
8
Constante de desoxigenacin (K 1 )
Mtodo
Tramo
1 Campaa
2 Campaa
3 Campaa
4 Campaa
5 Campaa
6 Campaa
7 Campaa
8 Campaa
9 Campaa
1.25
1.79
1.60
1.35
1.79
No cumple
1.05
1.79
1.27
Chapra
1.17
1.35
1.20
1.25
1.49
1.00
1.25
1.25
0.96
(1997)
1.05
1.35
1.20
1.17
1.46
1.00
1.17
1.49
0.88
0.96
1.33
1.00
1.10
1.21
1.00
1.17
1.25
0.88
Tabla7-15. Constantes de desoxigenacin para cada tramo de cada campaa, a partir de la ecuacin Chapra.
Fuente: Elaboracin propia.
A continuacin se presenta la grafica de la campaa que tuvo una mejor corrida, sin embargo
ninguna corrida calibro correctamente debido al gran demanda de datos necesarios para
utilizar el modelo.
Pgina | 132
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
SIMULACION DE O.D.
(Kd: Ec. ajustada y K2: O`Connor et al)
12
11
10
9
8
7
6
5
4
3
2
1
0
SIMULACION DBO5
(Kd: Ec. ajustada y K2: O`Connor et al)
28
25
23
20
OD (mg/L)
18
DBO (mg/L)
15
13
10
8
5
3
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
O.D. Observado (mg/L)
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
X (Km)
t (d)
DBO (mg/l)
OD (mg/l)
Sat OD (mg/l)
0.00
5.4
7.3
11.0395
0.5
0.02
5.3515
7.6779
11.0395
0.04
5.3035
7.6928
11.0395
1.5
0.06
5.2559
7.7596
11.0395
0.08
5.2087
7.9542
11.0395
2.5
0.10
5.1853
8.5321
11.0395
0.12
5.1387
8.6734
11.0395
3.5
0.14
5.1156
8.7453
11.0395
0.16
5.0926
8.7666
11.0395
4.5
0.18
5.0697
8.9061
11.0395
0.20
5.0469
9.3004
11.0395
5.33
0.22
5.0279
9.4972
11.0395
5.33
0.24
24
9.7003
11.0395
5.83
0.26
23.3438
9.9545
11.0395
6.33
0.30
22.7053
9.9869
11.0395
6.83
0.32
22.0844
10.1166
11.0395
7.33
0.34
21.4811
10.2192
11.0395
7.83
0.38
20.8935
10.2844
10.1381
8.33
0.40
20.3217
10.3145
10.1381
8.83
0.42
19.7666
10.3258
10.1381
9.33
0.46
19.2253
10.3522
10.1381
9.83
0.48
18.6999
10.5346
10.1381
10.33
0.50
18.1884
10.6819
10.1381
10.83
0.54
17.6903
10.7807
10.1381
11.33
0.56
17.2072
10.8471
10.1381
11.83
0.58
16.7365
10.8915
10.1381
12.33
0.60
16.2784
10.9214
10.1381
12.83
0.64
15.8342
10.9415
10.1381
13.33
0.66
15.4712
10.955
10.1381
13.83
0.68
16.6669
10.9629
10.1381
14.33
0.70
16.4578
3.4479
10.1381
Pgina | 133
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
14.83
0.72
16.1008
5.5343
10.1381
15.33
0.76
15.7481
6.8351
10.1381
15.83
0.78
15.4003
7.6428
10.2267
16.33
0.80
15.0578
8.15
10.2267
16.83
0.82
14.7218
8.4727
10.2267
17.33
0.84
14.3907
8.5321
10.2267
17.83
0.88
14.0664
8.6819
10.2267
18.33
0.90
13.7478
8.822
10.2267
18.83
0.92
13.4349
8.919
10.2267
19.33
0.94
13.1291
8.9061
10.2267
19.83
0.96
12.8289
8.9898
10.2267
20.33
1.00
12.5346
9.044
10.2267
20.83
1.02
12.2474
9.0878
10.2267
21.33
1.04
12.0114
9.1251
10.2267
21.83
1.06
0.6597
9.1501
10.2267
22.33
1.08
2.5991
9.158
10.2267
22.83
1.10
2.8673
9.1879
10.2267
23.33
1.12
3.0807
9.2114
10.6431
23.83
1.14
3.2502
9.3123
10.6431
24.33
1.16
3.3827
9.3123
10.6431
24.83
1.18
3.485
9.4216
10.6431
25.33
1.20
3.562
9.4974
10.6431
25.83
1.22
3.6175
9.5516
10.6431
26.33
1.24
3.6552
9.5919
10.6431
26.83
1.26
3.6777
9.6232
10.6431
27.33
1.28
3.6874
9.6487
10.6431
27.83
1.30
3.687
9.7003
10.6431
28.33
1.32
0.0507
9.9545
10.6431
28.83
1.34
0.2078
10.1166
10.6431
29.33
1.36
0.4272
10.1166
10.5716
29.83
1.38
1.1834
10.2192
10.5716
30.33
1.40
1.4207
10.2844
10.5716
30.83
1.42
1.6407
10.3258
10.5716
31.33
1.44
1.8403
10.3522
10.5716
31.83
1.46
2.0192
10.3675
10.5716
32.33
1.48
2.1778
10.3876
10.5716
32.83
1.50
2.3162
10.3876
10.5716
33.33
1.52
2.4362
10.3985
10.5716
33.83
1.54
2.5389
10.4039
10.5716
34.33
1.56
2.6252
10.406
10.5716
34.43
1.58
2.6863
10.4063
10.5716
Pgina | 134
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
7.1.2. VALIDACIN
Para validar los modelos, se realizan corridas con los modelos ya calibrados con distintas
campaas utilizadas en la calibracin. Los resultados son presentados en los grficos y tablas
presentados en los subttulos a continuacin.
El modelo que represente mejores resultados se utilizara para realizar los escenarios de
prediccin.
7.1.2.1.MODELO SIMOD
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O`Connor et al)
25
23
20
18
15
DBO (mg/L)
OD (mg/L)
SIMULACION DBO5
(Kd: M. Bansal y K2: O`Connor et al)
13
10
2
1
8
5
3
0
0
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia Km
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
N
Segmento
Distancia
Acumulada(Km)
Temperatura
del Agua
(C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D.
Mnimo
Permitido
(mg/L)
1
1
1
1
1
1
1
2
2
3
3
4
4
4
4
0
1.00
2.00
3.00
4.00
5.00
5.33
5.33
5.33
5.33
5.62
5.62
15.62
25.62
34.43
11.10
11.10
11.10
11.10
11.10
11.10
11.10
14.32
14.32
14.32
14.32
14.32
14.32
14.32
14.32
6.89
6.84
6.81
6.79
6.78
6.77
6.77
3.07
3.07
3.07
3.08
3.08
3.25
3.41
3.55
5.40
5.38
5.35
5.33
5.31
5.28
5.28
21.13
21.13
22.11
22.04
21.23
20.13
19.08
18.21
6.89
6.89
6.89
6.89
6.89
6.89
6.89
6.40
6.40
6.41
6.41
6.41
6.41
6.41
6.41
0.12
0.12
0.12
0.12
0.12
0.12
0.12
0.60
0.60
1.05
1.05
0.64
0.64
0.64
0.64
10.06
10.06
10.06
10.06
10.06
10.06
10.06
4.61
4.61
8.06
8.06
4.84
4.84
4.84
4.84
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
O.D.
Observado
(mg/L)
DBO 5
Observada
(mg/L)
7.3
5.40
3.15
3.5
18.00
15.54
Pgina | 135
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O'Connor et al)
SIMULACION DBO5
(Kd: M. Bansal y K2: O'Connor et al)
7
6
DBO (mg/L)
OD (mg/L)
5
4
3
2
1
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
O.D. Observado (mg/L)
Distancia Km
55
50
45
40
35
30
25
20
15
10
5
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
2010
N
Segmento
Distancia
Acumulada(Km)
Temperatura
del Agua
(C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D.
Mnimo
Permitido
(mg/L)
1
1
1
1
1
1
1
2
2
3
3
4
4
4
4
0.00
1.00
2.00
3.00
4.00
5.00
5.33
5.33
5.33
5.33
5.62
5.62
15.62
25.62
34.43
16.00
16.00
16.00
16.00
16.00
16.00
16.00
17.08
17.08
17.08
17.08
17.08
17.08
17.08
17.08
5.10
5.75
6.01
6.11
6.15
6.16
6.17
2.96
2.95
2.95
3.04
3.04
4.01
4.67
5.07
6.00
2.37
0.94
0.37
0.15
0.06
0.04
28.48
28.47
41.84
40.61
38.34
25.92
17.52
12.41
6.18
6.18
6.18
6.18
6.18
6.18
6.18
6.04
6.04
6.04
6.04
6.04
6.04
6.04
6.04
0.40
0.40
0.40
0.40
0.40
0.40
0.40
0.18
0.18
0.40
0.40
0.34
0.34
0.34
0.34
2.97
2.97
2.97
2.97
2.97
2.97
2.97
1.62
1.62
6.64
6.64
5.62
5.62
5.62
5.62
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
4.00
O.D.
Observado
(mg/L)
DBO 5
Observada
(mg/L)
5.1
6.00
2.9
43.00
4
4.3
17.00
15.00
Pgina | 136
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
7.1.2.2.MODELO RIOSep
SIMULACION DE O.D.
(Kd: Ec. ajustada y K2: O`Connor et al)
25
23
20
18
15
DBO (mg/L)
OD (mg/L)
SIMULACION DBO5
(Kd: Ec. ajustada y K2: O`Connor et al)
4
3
13
10
8
5
3
2
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
O.D. Observado (mg/L)
Distancia Km
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en DBO
KmRIOSep (mg/l)
X (Km)
t (h)
DBO 5 (mg/l)
DOD (mg/l)
OD (mg/l)
OD saturacin
5.4
0.0031
6.15
6.1531
0.5
0.4355
5.3578
0.0407
6.1124
6.1531
0.8711
5.316
0.0712
6.0818
6.153
1.5
1.3066
5.2745
0.0959
6.0571
6.153
1.7422
5.2333
0.1159
6.0372
6.1531
2.5
2.1777
5.1925
0.1319
6.0212
6.1531
2.6132
5.1519
0.1447
6.0083
6.153
3.5
3.0488
5.1117
0.155
5.9981
6.1531
3.4843
5.0718
0.163
5.99
6.153
4.5
3.9199
5.0322
0.1694
5.9837
6.1531
4.3554
4.9929
0.1743
5.9788
6.1531
5.33
4.6429
4.9671
0.1769
5.9762
6.1531
5.33
21.0828
2.7189
2.9481
5.667
5.83
0.2016
20.9849
2.6012
3.0658
5.667
6.33
0.4033
20.8875
2.4923
3.1747
5.667
6.83
0.6049
20.7906
2.3914
3.2756
5.667
7.33
0.8066
20.6941
2.2979
3.3691
5.667
7.83
1.0082
20.598
2.2113
3.4557
5.667
8.33
1.2098
20.5024
2.131
3.536
5.667
8.83
1.4115
20.4072
2.0565
3.6105
5.667
9.33
1.6131
20.3125
1.9874
3.6796
5.667
9.83
1.8148
20.2182
1.9232
3.7438
5.667
10.33
2.0164
20.1243
1.8636
3.8034
5.667
10.83
2.218
20.0309
1.8083
3.8587
5.667
11.33
2.4197
19.9379
1.7568
3.9102
5.667
11.83
2.6213
19.8453
1.7089
3.9581
5.667
12.33
2.8229
19.7532
1.6644
4.0026
5.667
12.83
3.0246
19.6615
1.6229
4.0441
5.667
13.33
3.2262
19.5702
1.5842
4.0828
5.667
13.83
3.4279
19.4794
1.5481
4.1189
5.667
14.33
3.6295
19.389
1.5145
4.1525
5.667
Pgina | 137
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
14.83
3.8311
19.299
1.483
4.184
5.667
15.33
4.0328
19.2094
1.4536
4.2134
5.667
15.83
4.2344
19.1202
1.4261
4.2409
5.667
16.33
4.4361
19.0314
1.4003
4.2667
5.667
16.83
4.6377
18.9431
1.3762
4.2908
5.667
17.33
4.8393
18.8552
1.3535
4.3135
5.667
17.83
5.041
18.7676
1.3322
4.3348
5.667
18.33
5.2426
18.6805
1.3122
4.3548
5.667
18.83
5.4443
18.5938
1.2933
4.3737
5.667
19.33
5.6459
18.5075
1.2755
4.3915
5.667
19.83
5.8475
18.4215
1.2587
4.4083
5.667
20.33
6.0492
18.336
1.2429
4.4241
5.667
20.83
6.2508
18.2509
1.2279
4.4391
5.667
21.33
6.4525
18.1662
1.2137
4.4533
5.667
21.83
6.6541
18.0819
1.2002
4.4668
5.667
22.33
6.8557
17.9979
1.1874
4.4796
5.667
22.83
7.0574
17.9144
1.1752
4.4918
5.667
23.33
7.259
17.8312
1.1636
4.5034
5.667
23.83
7.4607
17.7484
1.1525
4.5145
5.667
24.33
7.6623
17.666
1.1419
4.5251
5.667
24.83
7.8639
17.584
1.1318
4.5352
5.667
25.33
8.0656
17.5024
1.1221
4.5449
5.667
25.83
8.2672
17.4211
1.1128
4.5543
5.6671
26.33
8.4688
17.3403
1.1038
4.5632
5.667
26.83
8.6705
17.2598
1.0952
4.5718
5.667
27.33
8.8721
17.1797
1.0869
4.5801
5.667
27.83
9.0738
17.0999
1.0789
4.5881
5.667
28.33
9.2754
17.0205
1.0711
4.5959
5.667
28.83
9.477
16.9415
1.0636
4.6034
5.667
29.33
9.6787
16.8629
1.0564
4.6106
5.667
29.83
9.8803
16.7846
1.0493
4.6177
5.667
30.33
10.082
16.7067
1.0425
4.6246
5.6671
30.83
10.2836
16.6291
1.0358
4.6312
5.667
31.33
10.4852
16.5519
1.0293
4.6377
5.667
31.83
10.6869
16.4751
1.023
4.644
5.667
32.33
10.8885
16.3986
1.0168
4.6502
5.667
32.83
11.0902
16.3225
1.0107
4.6563
5.667
33.33
11.2918
16.2467
1.0048
4.6622
5.667
33.83
11.4934
16.1713
0.9991
4.668
5.6671
34.33
11.6951
16.0962
0.9934
4.6736
5.667
34.43
11.7354
16.0812
0.9923
4.6747
5.667
Pgina | 138
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
SIMULACION DE O.D.
(Kd: Ec. ajustada y K2: O'Connor et al)
SIMULACION DBO5
(Kd: Ec. ajustada y K2: O'Connor et al)
55
50
45
40
35
30
DBO (mg/L)
OD (mg/L)
3
2
1
25
20
15
10
5
0
0
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia Km
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
OD RIOSep (mg/l)
X (Km)
t (h)
DBO 5 (mg/l)
DOD (mg/l)
OD (mg/l)
OD saturacin
0.4203
5.1
5.5203
0.5
25.2778
5.6926
0.6767
4.8436
5.5203
50.5556
5.401
0.8955
4.6249
5.5204
1.5
75.8333
5.1243
1.0807
4.4397
5.5204
101.1111
4.8618
1.236
4.2844
5.5204
2.5
126.3889
4.6127
1.3647
4.1556
5.5203
151.6667
4.3764
1.4699
4.0504
5.5203
3.5
176.9444
4.1522
1.5543
3.966
5.5203
202.2222
3.9395
1.6203
3.9001
5.5204
4.5
227.5
3.7377
1.67
3.8503
5.5203
252.7778
3.5462
1.7055
3.8148
5.5203
5.33
269.4611
3.4252
1.7219
3.7984
5.5203
5.33
28.5408
2.4924
2.9161
5.4085
5.83
6.8572
28.0456
2.4893
2.9193
5.4086
6.33
13.7143
27.5589
2.4789
2.9297
5.4086
6.83
20.5715
27.0806
2.4627
2.9458
5.4085
7.33
27.4287
26.6107
2.442
2.9666
5.4086
7.83
34.2859
26.1489
2.4176
2.9909
5.4085
8.33
41.143
25.6951
2.3904
3.0181
5.4085
8.83
48.0002
25.2492
2.361
3.0476
5.4086
9.33
54.8574
24.811
2.3299
3.0786
5.4085
9.83
61.7145
24.3805
2.2976
3.111
5.4086
10.33
68.5717
23.9574
2.2643
3.1442
5.4085
10.83
75.4289
23.5416
2.2304
3.1781
5.4085
11.33
82.2861
23.1331
2.1962
3.2124
5.4086
11.83
89.1432
22.7317
2.1617
3.2468
5.4085
12.33
96.0004
22.3372
2.1272
3.2814
5.4086
12.83
102.8576
21.9495
2.0927
3.3159
5.4086
13.33
109.7147
21.5686
2.0584
3.3502
5.4086
Pgina | 139
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
13.83
116.5719
21.1944
2.0243
3.3843
5.4086
14.33
123.4291
20.8266
1.9905
3.4181
5.4086
14.83
130.2863
20.4651
1.957
3.4515
5.4085
15.33
137.1434
20.11
1.924
3.4846
5.4086
15.83
144.0006
19.761
1.8913
3.5172
5.4085
16.33
150.8578
19.4181
1.8591
3.5495
5.4086
16.83
157.7149
19.0811
1.8273
3.5812
5.4085
17.33
164.5721
18.75
1.796
3.6125
5.4085
17.83
171.4293
18.4246
1.7652
3.6434
5.4086
18.33
178.2865
18.1049
1.7348
3.6737
5.4085
18.83
185.1436
17.7907
1.7049
3.7036
5.4085
19.33
192.0008
17.4819
1.6755
3.733
5.4085
19.83
198.858
17.1786
1.6466
3.762
5.4086
20.33
205.7151
16.8805
1.6181
3.7904
5.4085
20.83
212.5723
16.5875
1.5902
3.8184
5.4086
21.33
219.4295
16.2997
1.5626
3.8459
5.4085
21.83
226.2867
16.0168
1.5356
3.873
5.4086
22.33
233.1438
15.7389
1.509
3.8995
5.4085
22.83
240.001
15.4657
1.4829
3.9257
5.4086
23.33
246.8582
15.1973
1.4572
3.9514
5.4086
23.83
253.7153
14.9336
1.4319
3.9766
5.4085
24.33
260.5725
14.6745
1.4071
4.0015
5.4086
24.83
267.4297
14.4198
1.3827
4.0259
5.4086
25.33
274.2869
14.1696
1.3587
4.0498
5.4085
25.83
281.144
13.9237
1.3351
4.0734
5.4085
26.33
288.0012
13.6821
1.312
4.0966
5.4086
26.83
294.8584
13.4446
1.2892
4.1193
5.4085
27.33
301.7155
13.2113
1.2669
4.1417
5.4086
27.83
308.5727
12.982
1.2449
4.1637
5.4086
28.33
315.4299
12.7568
1.2233
4.1853
5.4086
28.83
322.2871
12.5354
1.2021
4.2065
5.4086
29.33
329.1442
12.3178
1.1812
4.2273
5.4085
29.83
336.0014
12.1041
1.1607
4.2478
5.4085
30.33
342.8586
11.894
1.1406
4.268
5.4086
30.83
349.7157
11.6876
1.1208
4.2878
5.4086
31.33
356.5729
11.4848
1.1013
4.3072
5.4085
31.83
363.4301
11.2855
1.0822
4.3263
5.4085
32.33
370.2873
11.0897
1.0634
4.3451
5.4085
32.83
377.1444
10.8972
1.045
4.3636
5.4086
33.33
384.0016
10.7081
1.0268
4.3817
5.4085
33.83
390.8588
10.5223
1.009
4.3995
5.4085
34.33
397.7159
10.3397
0.9915
4.417
5.4085
34.43
399.0874
10.3035
0.9881
4.4205
5.4086
Pgina | 140
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Distancia
acumulada (Km)
DBO 5 (mg/l)
SIMOD K 1 : M.
Bansal y K 2 :
O'Connor et al
DBO 5 (mg/l)
RIOSep K 1 : Ec.
Ajustada y K 2 :
O'Connor y
Dobbins
DBO 5 Observado
(mg/l)
0.00
16.00
16.00
16.00
28.67
30.03
23.71
25.00
34.43
28.25
21.86
24.00
0.00
22.00
22.00
22.00
28.67
41.80
33.95
34.00
34.43
38.79
31.18
30.00
0.00
16.00
16.00
16.00
28.67
48.88
23.64
39.00
34.43
43.39
19.98
33.00
0.00
12.00
12.00
12.00
28.67
32.33
10.17
39.00
34.43
26.24
7.44
38.00
0.00
15.90
15.90
15.92
5.33
65.15
53.61
70.00
28.67
60.00
40.16
46.84
34.43
57.86
37.48
44.13
PROMEDIO
Coeficiente de
correlacin entre
los valores del
modelo SIMOD y
datos Observados
Coeficiente de
correlacin entre
los valores del
modelo RIOSep y
datos Observados
0.9998
0.9829
0.9636
0.9878
0.9905
0.9410
0.9320
0.6152
0.8337
0.9924
0.9439
0.9039
Pgina | 141
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Distancia
acumulada (Km)
OD (mg/l) SIMOD
K 1 : M. Bansal y
K 2 : O'Connor et al
OD (mg/l) RIOSep
K 1 : Ec. Ajustada y
K 2 : O'Connor y
Dobbins
OD Observado
(mg/l)
0.00
5.45
5.45
5.45
28.67
4.26
3.70
4.50
34.43
4.33
3.80
4.50
0.00
6.02
6.20
6.20
28.67
3.45
3.82
4.30
34.43
3.66
3.97
5.30
0.00
6.00
6.00
5.97
28.67
3.67
3.74
3.50
34.43
3.88
3.92
4.00
0.00
5.00
5.00
5.00
28.67
3.02
4.54
3.98
34.43
3.53
4.69
4.15
0.00
5.40
5.40
5.40
5.33
4.36
4.34
4.50
28.67
4.43
4.22
4.12
34.43
4.48
4.28
4.50
PROMEDIO
Coeficiente de
correlacin entre
los valores del
modelo SIMOD y
datos Observados
Coeficiente de
correlacin entre
los valores del
modelo RIOSep y
datos Observados
0.9975
0.9973
0.7844
0.7734
0.9875
0.9861
0.9918
0.9738
0.8794
0.9331
0.9281
0.9327
Comparando los valores de los coeficientes de correlacin en las corridas realizadas para el
calibrado podemos observar que los valores ms cercanos a 1 son los del modelo SIMOD en
los que respecta al modelado de la DBO 5 (tabla 7-27), el que tiene mejor correlacin del
modelado de OD es el modelo RIOSep, ya que sus valores son ms cercanos a 1 (tabla 7-28).
Los valores del modelo WASP 7.4 no se tomaron en cuenta debido a la complejidad del
modelo no se puede calibrar con los datos que se posee, por tal motivo no se considera
necesario analizar el coeficiente de correlacin de este modelo.
7.2.2. ERRORES E INCERTIDUMBRE
Se conoce como valor verdadero () de una magnitud, a la cantidad que expresa el valor
absoluto de la misma. Su determinacin por medicin en la mayora de la generalidad de los
casos es muy difcil y solo se conocen aproximaciones.
Si se compara los valores obtenidos con los dos modelos con las diferentes campaas con los
valores observados, considerando a estos ltimos como valores verdaderos () se puede
Pgina | 142
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Donde:
. (7 4)
. (7 5)
. (7 6)
% = 100
. (7 7)
E = Error absoluto
xi = Valor generado por el modelo
= Valor observado
= Error relativo
% = Error relativo porcentual
Distancia
acumulada
(Km)
DBO 5 (mg/l)
SIMOD K 1 : M.
Bansal y K 2 :
O'Connor et al
DBO 5
Observado
(mg/l)
Error
Absoluto
Error
Relativo
0.00
16.00
16.00
0.00
0.00
16.00
+/-
0.00
0.00
28.67
30.03
25.00
5.03
0.17
30.03
.+/-
0.17
16.76
34.43
28.25
24.00
4.25
0.15
28.25
.+/-
0.15
15.04
0.00
22.00
22.00
0.00
0.00
22.00
.+/-
0.00
0.00
28.67
41.80
34.00
7.80
0.19
41.80
.+/-
0.19
18.67
34.43
38.79
30.00
8.79
0.23
38.79
.+/-
0.23
22.67
0.00
16.00
16.00
0.00
0.00
16.00
.+/-
0.00
0.00
28.67
48.88
39.00
9.88
0.20
48.88
.+/-
0.20
20.21
34.43
43.39
33.00
10.39
0.24
43.39
.+/-
0.24
23.94
0.00
12.00
12.00
0.00
0.00
12.00
.+/-
0.00
0.00
28.67
32.33
39.00
6.67
0.21
32.33
.+/-
0.21
20.61
34.43
26.24
38.00
11.76
0.45
26.24
.+/-
0.45
44.83
0.00
15.90
15.92
0.02
0.00
15.90
.+/-
0.00
0.13
5.33
53.54
70.00
16.46
0.31
53.54
.+/-
0.31
30.74
28.67
60.00
46.84
13.16
0.22
60.00
.+/-
0.22
21.93
34.43
57.86
44.13
13.73
0.24
57.86
.+/-
0.24
23.73
Grado de
Incertidumbre
Error
porcentual
%
Promedio de
Errores
porcentuales
%
10.60
13.78
14.72
21.81
19.13
PROMEDIO
16.01
Tabla7-23. Errores e incertidumbre para las corridas mejor calibradas de DBO 5 , Modelo SIMOD.
Fuente: Elaboracin Propia.
Pgina | 143
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Distancia
acumulada
(Km)
DBO 5 (mg/l)
RIOSep K 1 : Ec.
Ajustada y K 2 :
O'Connor y
Dobbins
DBO 5
Observado
(mg/l)
Error
Absoluto
Error
Relativo
0.00
16.00
16.00
0.00
0.00
28.67
34.43
0.00
28.67
23.71
21.86
22.00
33.95
25.00
24.00
22.00
34.00
1.29
2.14
0.00
0.05
0.05
0.10
0.00
0.00
Grado de
Incertidumbre
Error
porcentual
%
16.00
+/-
0.00
0.00
23.71
+/-
0.05
5.42
21.86
+/-
0.10
9.79
22.00
+/-
0.00
0.00
33.95
+/-
0.00
0.15
5.07
1.31
34.43
31.18
30.00
1.18
0.04
31.18
+/-
0.04
3.79
0.00
16.00
16.00
0.00
0.00
16.00
+/-
0.00
0.00
28.67
23.64
39.00
15.36
0.65
23.64
+/-
0.65
64.97
19.98
+/-
0.65
65.17
12.00
+/-
0.00
0.00
10.17
+/-
2.84
283.57
7.44
+/-
4.11
410.75
15.90
+/-
0.00
0.13
34.43
0.00
28.67
34.43
0.00
19.98
12.00
10.17
7.44
15.90
33.00
12.00
39.00
38.00
15.92
13.02
0.00
28.83
30.56
0.02
0.65
0.00
2.84
4.11
0.00
5.33
53.61
70.00
16.39
0.31
53.61
+/-
0.31
30.58
28.67
40.16
46.84
6.68
0.17
40.16
+/-
0.17
16.65
34.43
37.48
44.13
6.65
0.18
37.48
+/-
0.18
17.73
Promedio de
Errores
porcentuales
%
43.38
231.44
16.27
PROMEDIO
59.49
Tabla7-24. Errores e incertidumbre para las corridas mejor calibradas de DBO 5 , Modelo RIOSep.
Fuente: Elaboracin Propia.
Distancia
acumulada
(Km)
OD (mg/l)
SIMOD K 1 : M.
Bansal y K 2 :
O'Connor et al
OD
Observado
(mg/l)
Error
Absoluto
Error
Relativo
0.00
5.45
5.45
0.00
0.00
5.45
.+/-
0.00
0.05
28.67
4.26
4.50
0.24
0.06
4.26
.+/-
0.06
5.60
34.43
4.33
4.50
0.17
0.04
4.33
.+/-
0.04
3.99
Grado de
Incertidumbre
Error
porcentual
%
0.00
6.02
6.20
0.18
0.03
6.02
.+/-
0.03
2.98
28.67
3.45
4.30
0.85
0.25
3.45
.+/-
0.25
24.53
34.43
3.66
5.30
1.64
0.45
3.66
.+/-
0.45
44.97
0.00
6.00
5.97
0.03
0.01
6.00
.+/-
0.01
0.50
28.67
3.67
3.50
0.17
0.05
3.67
.+/-
0.05
4.59
34.43
3.88
4.00
0.12
0.03
3.88
.+/-
0.03
3.19
0.00
5.00
5.00
0.00
0.00
5.00
.+/-
0.00
0.00
28.67
3.02
3.98
0.96
0.32
3.02
.+/-
0.32
31.60
34.43
3.53
4.15
0.62
0.18
3.53
.+/-
0.18
17.70
0.00
5.40
5.40
0.00
0.00
5.40
.+/-
0.00
0.00
5.33
4.36
4.50
0.14
0.03
4.36
.+/-
0.03
3.11
28.67
4.43
4.12
0.31
0.07
4.43
.+/-
0.07
6.95
34.43
4.48
4.50
0.02
0.01
4.48
.+/-
0.01
0.53
Promedio de
Errores
porcentuales
%
3.21
24.16
2.76
16.43
2.65
PROMEDIO
9.84
Tabla7-25. Errores e incertidumbre para las corridas mejor calibradas de OD, Modelo SIMOD.
Fuente: Elaboracin Propia.
Pgina | 144
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Distancia
acumulada
(Km)
OD (mg/l)
RIOSep K 1 : Ec.
Ajustada y K 2 :
O'Connor y
Dobbins
OD
Observado
(mg/l)
Error
Absoluto
Error
Relativo
0.00
5.45
5.45
0.00
0.00
5.45
+/-
0.00
0.00
28.67
3.70
4.50
0.80
0.22
3.70
+/-
0.22
21.60
34.43
3.80
4.50
0.70
0.18
3.80
+/-
0.18
18.32
0.00
6.20
6.20
0.00
0.00
6.20
+/-
0.00
0.00
28.67
3.82
4.30
0.48
0.13
3.82
+/-
0.13
12.70
34.43
3.97
5.30
1.33
0.34
3.97
+/-
0.34
33.51
0.00
6.00
5.97
0.03
0.01
6.00
+/-
0.01
0.50
28.67
3.74
3.50
0.24
0.06
3.74
+/-
0.06
6.33
34.43
3.92
4.00
0.08
0.02
3.92
+/-
0.02
1.96
0.00
5.00
5.00
0.00
0.00
5.00
+/-
0.00
0.00
28.67
4.54
3.98
0.56
0.12
4.54
+/-
0.12
12.32
34.43
4.69
4.15
0.54
0.11
4.69
+/-
0.11
11.45
0.00
5.40
5.40
0.00
0.00
5.40
+/-
0.00
0.00
Grado de
Incertidumbre
Error
porcentual
%
5.33
4.34
4.50
0.16
0.04
4.34
+/-
0.04
3.75
28.67
4.22
4.12
0.10
0.02
4.22
+/-
0.02
2.30
34.43
4.28
4.50
0.22
0.05
4.28
+/-
0.05
5.11
Promedio de
Errores
porcentuales
%
13.31
15.40
2.93
7.92
2.79
PROMEDIO
8.47
Tabla7-26. Errores e incertidumbre para las corridas mejor calibradas de OD, Modelo RIOSep.
Fuente: Elaboracin Propia.
El modelo SIMOD es el que presenta un porcentaje de error mucho menor que el modelo
RIOSep en lo que refiere al modelado de al DBO 5 , pero el modelo RIOSep muestra un menor
porcentaje de error en la modelacin de OD no siendo mucha la diferencia con el modelo
SIMOD.
Los valores del modelo WASP 7.4 no se tomaron en cuenta debido a la complejidad del
modelo no se puede calibrar con los datos que se posee, por tal motivo no se considera
necesario analizar el coeficiente de correlacin de este modelo.
Pgina | 145
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
=
Donde:
. (7 8)
x o = Valor observado
x g = Valor generado por el modelo
n = Numero de puntos observados
Los clculos de la diferencia promedio se presentan en las siguientes tablas.
Distancia
acumulada
(Km)
DBO 5
Observado
(mg/l)
0.00
16.00
16.00
16.00
28.67
30.03
23.71
25.00
34.43
28.25
21.86
24.00
0.00
22.00
22.00
22.00
28.67
41.80
33.95
34.00
34.43
38.79
31.18
30.00
0.00
16.00
16.00
16.00
28.67
48.88
23.64
39.00
34.43
43.39
19.98
33.00
0.00
12.00
12.00
12.00
28.67
32.33
10.17
39.00
34.43
26.24
7.44
38.00
0.00
15.90
15.90
15.92
5.33
53.54
53.61
70.00
28.67
60.00
40.16
46.84
34.43
57.86
37.48
44.13
PROMEDIO
3.094
1.142
5.533
0.411
6.755
9.460
6.143
19.797
10.843
7.436
6.473
7.649
Pgina | 146
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Distancia
acumulada
(Km)
OD (mg/l) SIMOD
K 1 : M. Bansal y K 2 :
O'Connor et al
OD (mg/l) RIOSep K 1 :
Ec. Ajustada y K 2 :
O'Connor y Dobbins
OD
Observado
(mg/l)
0.00
5.45
5.45
5.45
28.67
4.26
3.70
4.50
34.43
4.33
3.80
4.50
0.00
6.02
6.20
6.20
28.67
3.45
3.82
4.30
34.43
3.66
3.97
5.30
0.00
6.00
6.00
5.97
28.67
3.67
3.74
3.50
34.43
3.88
3.92
4.00
0.00
5.00
5.00
5.00
28.67
3.02
4.54
3.98
34.43
3.53
4.69
4.15
0.00
5.40
5.40
5.40
5.33
4.36
4.34
4.50
28.67
4.43
4.22
4.12
34.43
4.48
4.28
4.50
PROMEDIO
0.138
0.499
0.890
0.605
0.107
0.114
0.527
0.365
0.117
0.120
0.356
0.341
El criterio de la Diferencia Promedio indica que el valor ms prximo a cero tiene una mejor
correlacin de los valores analizados, de acuerdo a esto los valores generados por el modelo
SIMOD presenta una mejor correlacin en el modelado de la DBO 5 comparado con el modelo
RIOSep, sin embargo el modelo RIOSep tiene una mejor correlacin en el modelado del OD
respecto al modelo SIMOD, sin embargo las diferencias no son muy grandes.
Los resultados estadsticamente obtenidos podran tender a ubicar los valores generados del
modelo SIMOD en primer lugar en comparacin con el modelo RIOSep, sin embargo, hay
que tener en cuenta que se asumi como valores verdaderos aquellos datos obtenidos en los
muestreos realizados sin considerar el margen de error que estos datos podran tener ya que
solo representan una medida del Ro Katari.
Pgina | 147
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Q M = Q1 + Q 2
Donde:
Q 2 = Q M Q1
Ec. (7 9)
Pgina | 148
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Q 2 = Caudal de descarga
Q M = Caudal aguas abajo
Qmax (m^3/s)
El alto
Viacha
Laja
Ro Seke
Ro Seco
Ro Pallina
Ro Pallina
0.417
1.103
1.002
1.038
Tabla7-29. Caudales mximos estimados de las descargas de la ciudad de El Alto, Viacha y Laja.
Fuente: Elaboracin Propia.
Para ver la relacin entre caudal, DBO 5 y OD en los Ro Seco y Pallina se tiene las siguientes
graficas:
DBO5 vs. Q
250
DBO5 (mg/l)
200
150
100
50
0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
Q (m^3/s)
Grafico7-30. Relacin de la DBO 5 vs. Q del Ro Seco con datos tomados por el IBTEN.
Fuente: Elaboracin Propia.
Obteniendo la siguiente ecuacin del ajuste, de datos obtenidos en el Ro Seco por el IBTEN:
DBO5 = 459.21 Q + 10.169
Michael Russel Osina Torrez
2 = 0.8360
(7 10)
Pgina | 149
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Como podemos observar en el grafico a medida que aumenta el caudal, aumenta la DBO 5
esto se debe a que el Ro Seco empieza en una vertiente con poco caudal y a medida que va
avanzando va recibiendo descargas contaminadas, lo cual aumenta el caudal y la DBO 5 .
DBO5 vs. OD
350
DBO5 (mg/l)
300
250
200
150
100
50
0
0
10
12
OD (mg/l)
Q<0.1
Q>0.1
Grafico7-31. Relacin de la DBO 5 vs. OD del Ro Seco con datos tomados por el IBTEN.
Fuente: Elaboracin Propia.
. (7 11)
. (7 12)
r 2 = 0.7112
r 2 = 0.6794
En la siguiente grafica es tiene la relacin de DBO 5 vs. OD la cual nos muestra que a medida
que aumenta la DBO 5 el OD disminuye. Esto se debe a que los organismos consumen OD
para poder degradar la materia orgnica y as dan valores elevados de DBO 5 .
Lo mismo se ve en los siguientes grficos del Ro Pallina el cual recoge las aguas del Ro
Seco, las descargas de Viacha y Laja que desemboca en el Ro Katari.
Pgina | 150
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
DBO5 vs. Q
100
DBO5 (mg/l)
80
60
40
20
0
0.00
0.02
0.04
0.06
0.08
0.10
Q (m^3/s)
Q<0.1
Lineal (Q<0.1)
Grafico7-32. Relacin de la DBO 5 vs. Q del Ro Pallina a partir de datos tomados por el IBTEN, para caudales menor a 0.1 (m^3/s).
Fuente: Elaboracin Propia.
Q < 0.1
r 2 = 0.9142
. (7 13)
DBO5 vs. Q
DBO5 (mg/l)
100
80
60
40
20
0
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0.70
0.80
Q (m^3/s)
0.1<Q>1
Lineal (0.1<Q>1)
Grafico7-33. Relacin de la DBO 5 vs. Q del Ro Pallina a partir de datos tomados por el IBTEN, para caudales mayor a 0.1 y menor a 1(m^3/s).
Fuente: Elaboracin Propia.
r 2 = 0.7835
. (7 14)
A medida que aumenta el caudal la carga de la DBO 5 va aumentando. Debido a que a medida
que se descarga las aguas residuales aumenta la DBO 5 y el caudal.
Michael Russel Osina Torrez
Pgina | 151
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Al ver que no exista una tendencia en los datos graficados de DBO 5 vs. Q del Ro Pallina,
resulto conveniente adoptar el criterio de relacionar DBO 5 y Q, graficando por rangos de
caudales.
DBO5 vs. OD
DBO5 (mg/l)
200
150
100
50
0
0.00
1.00
2.00
3.00
4.00
5.00
6.00
7.00
8.00
9.00
10.00
OD (mg/l)
Q<0.1
0.1<Q>1
Q>1
Polinmica (Q<0.1)
Polinmica (0.1<Q>1)
Polinmica (Q>1)
Grafico7-34. Relacin de la DBO 5 vs. OD del Ro Pallina a partir de datos tomados por el IBTEN, para caudales menor a 0.1, entre 0.11 y mayor a 1 (m^3/s).
Fuente: Elaboracin Propia.
Q < 0.1
Q>1
. (7 15)
. (7 16)
. (7 17)
Pgina | 152
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
AO
HABITANTES
CONECTADOS
CAUDALES
(m3/d)
DBO (mg/l)
OD (mg/l)
2009
2010
2011
2012
2013
2014
2015
2016
2017
2018
2019
2020
2021
2022
2023
2024
2025
2026
2027
2028
2029
2030
2031
2032
2033
2034
2035
571.868
617.818
665.579
714.904
765.497
817.027
869.125
921.395
973.423
1.024.790
1.075.078
1.123.888
1.170.845
1.215.616
1.257.908
1.297.489
1.334.182
1.367.874
1.419.430
1.472.963
1.528.551
1.584.572
1.642.643
1.702.840
1.765.243
1.829.932
1.896.996
39.476
43.075
46.869
50.845
54.988
59.276
63.687
68.192
72.763
77.369
81.977
86.556
91.074
95.502
99.813
103.983
107.993
111.827
117.203
122.839
128.75
134.803
141.141
147.776
154.724
161.997
169.614
478.063
473.314
468.625
463.998
459.391
454.855
450.343
445.888
441.474
437.100
432.780
428.486
424.248
420.044
415.888
411.769
407.693
403.659
399.657
395.705
391.783
387.907
384.063
380.265
376.496
372.772
369.079
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Tabla7-30. Datos de proyeccin de la carga que llegara a la PTAR Puchukollo hasta el 2035.
Fuente: EPSAS.
Pgina | 153
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
H vs. Q
0.6
0.5
H (m)
0.4
0.3
0.2
0.1
0.0
0.0
0.4
0.8
1.2
1.6
2.0
2.4
2.8
3.2
Q (m^3/s)
Grafico7-35. Relacin de la Profundidad vs. Caudal del Ro Pallina a partir de datos tomados por el IBTEN.
Fuente: Elaboracin Propia.
r 2 = 0.7279
. (7 18)
V vs. Q
1.2
V (m/s)
1.0
0.8
0.6
0.4
0.2
0.0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Q (m^3/s)
Grafico7-36. Relacin de la Velocidad vs. Caudal del Ro Pallina a partir de datos tomados por el IBTEN.
Fuente: Elaboracin Propia.
Pgina | 154
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
V = 0.5184 Q0.4007
r 2 = 0.8717
. (7 19)
Puesta en marcha
Hasta
2009
Etapa 1
2013
Etapa 2
2017
Etapa 3
2026
Etapa 4
2035
Inversin
3 Filtros percoladores.
2 Filtros percoladores,
desinfeccin.
2 Filtros percoladores, 1 Laguna
anaerbica, desinfeccin.
2 Filtros percoladores, 1 Laguna
anaerbica.
4 Filtros percoladores, 1 Laguna
anaerbica.
Colector
II=6
III=4
II=5
III=3
III=5
III=6
Filtros
percoladores
Serie III
III=1
II=4
Desinfeccion
Bombeo
III=2
0,25 Q
Recirculacion
Serie IIIII=3
Tanque de
ecualizacion
bombeo
0,25 Q
II=2
II=1
Lagunas
anaerobicas (nuevas)
Planta existente
Parschall 1
Canal de descarga
Parschall 2
Distribudor
Tamices
Afluente
Pgina | 155
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
El tercer escenario, ningn tipo de tratamiento. En este escenario se ver como incrementa la
contaminacin con el lapso de tiempo si es que no se realiza un plan de accin ambiental, que
pueda frenar esta constante contaminacin de los ros.
Michael Russel Osina Torrez
Pgina | 156
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Q (m^3/s)
DBO 5 (mg/l)
OD (mg/l)
Descripcin
Prom
Max
Min
Prom
Max
Min
Prom
Max
Min
0.128
0.474
0.004
3.45
5.00
0.00
56.12
176.00
5.00
0.227
1.151
0.021
4.00
5.75
0.10
43.55
185.00
1.00
0.322
0.435
0.245
3.20
4.70
0.30
140.00
235.00
71.00
0.457
1.117
0.162
1.45
2.85
0.00
113.31
198.00
6.00
0.083
0.174
0.028
1.49
6.60
0.00
68.16
107.13
30.00
1.207
3.237
0.129
2.61
6.25
1.00
44.98
89.00
7.00
Tabla7-35. Caudales estimados, concentraciones de OD y DBO 5 de descarga de las principales fuentes de contaminacin de la cuenca del Ro Katari.
Fuente: Elaboracin propia.
Como se puede observar, no se tiene datos de las descargas de la ciudad de Laja, esto se debe
a que no tiene una sola descarga si no varias pequeas descargas puntuales, las cuales no
estn bien identificadas, sin embargo en el muestreo realizado a fines de marzo del 2011 se
muestreo una de las descargas, y presentaba las siguientes caractersticas: 0 (mg/l) de OD,
0.01 (m^3/s) de Caudal y 45 (mg/l) de DBO 5 .
Por balance de masa podemos inferir que el caudal promedio que aporta el la ciudad de Laja
es de 0.249 (m^3/s), y que el valor de la DBO 5 es aproximadamente de 45 (mg/l) con un OD
de 0 (mg/l). El Ro Pallina llega a la confluencia con el Ro colorado con un valor promedio
de DBO 5 de 54.81 (mg/l) y 2.16 (mg/l) de OD. Con lo cual podemos estimar una remocin
natural de 32.35 %.
El Ro Seco descarga un caudal bastante representativo Clase D. Sabiendo que los valores de
DBO son altos antes de la descarga de la PTAR Puchukollo, entonces si existe una ampliacin
en la Planta no afectara de gran manera a mitigar la contaminacin de los Ros dentro de la
cuenca. Lo cual si los filtros percoladores funcionaran correctamente conjuntamente con la
construccin de las lagunas anaerobias la concentracin de la descarga tendra que disminuir
a 80 (mg/l) aproximadamente.
El Ro seco descargara sus aguas en el peor de los casos con una concentracin de DBO
aproximada de 132.5 (mg/l), ms las descargas del Ro Seke que si bien estn contaminadas
no tienen mucho caudal lo cual podra ascender a un mximo de 145.485 (mg/l). Aplicando el
Michael Russel Osina Torrez
Pgina | 157
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
OD (mg/l)
Q (m^3/s)
H (m)
V (m/s)
T (C)
98.42
0.57
1.248
0.35
0.57
14
Tabla7-36. Datos de ingreso al modelo para prediccin del primer escenario (descarga del Ro Pallina).
Fuente: Elaboracin propia.
90
80
70
60
50
DBO (mg/L)
OD (mg/L)
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O`Connor et al)
3
2
1
0
40
30
20
10
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
OD SIMOD (mg/L)
Distancia Km
SIMULACION DBO5
(Kd: M. Bansal y K2: O`Connor et al)
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
Escenario
Prediccin 1
Distancia
Acumulada(Km)
Temperatura
del Agua (C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D.
Mnimo
Permitido
(mg/L)
15.73
5.30
12.00
6.21
1.72
24.33
4.00
1.00
15.73
5.32
11.80
6.21
1.72
24.33
4.00
2.00
15.73
5.33
11.60
6.21
1.72
24.33
4.00
3.00
15.73
5.35
11.40
6.21
1.72
24.33
4.00
4.00
15.73
5.36
11.21
6.21
1.72
24.33
4.00
5.00
15.73
5.37
11.02
6.21
1.72
24.33
4.00
5.33
15.73
5.38
10.96
6.21
1.72
24.33
4.00
5.33
14.36
1.61
79.98
6.40
0.67
12.30
4.00
5.33
14.36
1.61
79.98
6.40
0.67
12.30
4.00
5.33
14.36
1.61
78.86
6.41
0.61
11.08
4.00
15.33
14.36
1.70
77.25
6.41
0.61
11.08
4.00
25.33
14.36
1.80
75.67
6.41
0.61
11.08
4.00
34.43
14.36
1.89
74.26
6.41
0.61
11.08
4.00
N
Segmento
Pgina | 158
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Segundo escenario: proponiendo tres tipos de tratamiento tendremos los siguientes datos del
Ro Pallina:
Tratamiento
% Remocin
DBO 5 (mg/l)
OD (mg/l)
Q (m^3/s)
H (m)
V (m/s)
T (C)
Primario
50
44.5
1.41
3.237
0.44
0.82
14
Secundario
80
17.8
3.25
3.237
0.44
0.82
14
Terciario
95
4.45
5.40
3.237
0.44
0.82
14
Tabla7-38. Datos de entrada al modelo segn la aplicacin de la eficiencia de los distintos tratamientos.
Fuente: Elaboracin Propia.
Tratamiento Primario
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O`Connor et al)
8
7
6
4
DBO (mg/L)
OD (mg/L)
5
3
2
1
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia Km
SIMULACION DBO5
(Kd: M. Bansal y K2: O`Connor et al)
48
44
40
36
32
28
24
20
16
12
8
4
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
Escenario
Distancia
Acumulada(Km)
Temperatura
del Agua (C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D. Mnimo
Permitido
(mg/L)
15.73
5.30
12.00
6.21
1.72
24.33
4.00
1.00
15.73
5.32
11.80
6.21
1.72
24.33
4.00
2.00
15.73
5.33
11.60
6.21
1.72
24.33
4.00
3.00
15.73
5.35
11.40
6.21
1.72
24.33
4.00
4.00
15.73
5.36
11.21
6.21
1.72
24.33
4.00
5.00
15.73
5.37
11.02
6.21
1.72
24.33
4.00
5.33
15.73
5.38
10.96
6.21
1.72
24.33
4.00
5.33
14.16
1.77
41.37
6.43
0.69
6.98
4.00
5.33
14.16
1.77
41.37
6.43
0.69
6.98
4.00
5.33
14.16
1.77
42.12
6.44
0.84
8.56
4.00
15.33
14.16
1.90
40.93
6.44
0.84
8.56
4.00
25.33
14.16
2.03
39.77
6.44
0.84
8.56
4.00
34.43
14.16
2.14
38.75
6.44
0.84
8.56
4.00
N
Segmento
Tabla7-39. Resultados Escenario de Prediccin 2 (Aplicando Tratamiento Primario) con el modelo SIMOD v1.01.
Pgina | 159
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tratamiento Secundario
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O`Connor et al)
20
18
16
14
12
DBO (mg/L)
OD (mg/L)
SIMULACION DBO5
(Kd: M. Bansal y K2: O`Connor et al)
3
2
1
10
8
6
4
2
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
OD SIMOD (mg/L)
Distancia Km
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
DBO SIMOD (mg/L)
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
Escenario
N
Segmento
Distancia
Acumulada(Km)
Temperatura
del Agua (C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D. Mnimo
Permitido
(mg/L)
15.73
5.30
12.00
6.21
1.72
24.33
4.00
1.00
15.73
5.32
11.80
6.21
1.72
24.33
4.00
2.00
15.73
5.33
11.60
6.21
1.72
24.33
4.00
3.00
15.73
5.35
11.40
6.21
1.72
24.33
4.00
4.00
15.73
5.36
11.21
6.21
1.72
24.33
4.00
5.00
15.73
5.37
11.02
6.21
1.72
24.33
4.00
5.33
15.73
5.38
10.96
6.21
1.72
24.33
4.00
5.33
14.16
3.49
17.16
6.43
0.82
5.67
4.00
5.33
14.16
3.49
17.16
6.43
0.82
5.67
4.00
5.33
14.16
3.49
17.34
6.44
1.00
6.95
4.00
15.33
14.16
3.59
16.76
6.44
1.00
6.95
4.00
25.33
14.16
3.69
16.19
6.44
1.00
6.95
4.00
34.43
14.16
3.77
15.70
6.44
1.00
6.95
4.00
Tabla7-40. Resultados Escenario de Prediccin 2 (Aplicando Tratamiento Secundario) con el modelo SIMOD v1.01.
Fuente: Modelo SIMOD v1.01.
Pgina | 160
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tratamiento Terciario
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O`Connor et al)
14
12
10
DBO (mg/L)
OD (mg/L)
SIMULACION DBO5
(Kd: M. Bansal y K2: O`Connor et al)
4
3
2
1
8
6
4
2
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
OD SIMOD (mg/L)
Distancia Km
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
N
Segmento
Temperatura
del Agua (C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D.
Mnimo
Permitido
(mg/L)
15.73
5.30
12.00
6.21
1.72
24.33
4.00
1.00
15.73
5.32
11.80
6.21
1.72
24.33
4.00
2.00
15.73
5.33
11.60
6.21
1.72
24.33
4.00
3.00
15.73
5.35
11.40
6.21
1.72
24.33
4.00
4.00
15.73
5.36
11.21
6.21
1.72
24.33
4.00
5.00
15.73
5.37
11.02
6.21
1.72
24.33
4.00
5.33
15.73
5.38
10.96
6.21
1.72
24.33
4.00
5.33
14.16
5.40
5.10
6.43
0.86
5.20
4.00
5.33
14.16
5.40
5.10
6.43
0.86
5.20
4.00
5.33
14.16
5.40
5.15
6.44
1.06
6.35
4.00
15.33
14.16
5.43
4.96
6.44
1.06
6.35
4.00
25.33
14.16
5.47
4.79
6.44
1.06
6.35
4.00
34.43
14.16
5.50
4.63
6.44
1.06
6.35
4.00
Tabla7-41. Resultados Escenario de Prediccin 2 (Aplicando Tratamiento Terciario) con el modelo SIMOD v1.01.
Fuente: Modelo SIMOD v1.01.
Pgina | 161
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Tercer escenario: sin ningn tratamiento se llegaran a valores mximos de DBO 5 y valores
bajos de OD. Por lo cual se tiene los siguientes datos para el Ro Pallina:
DBO 5 (mg/l)
OD (mg/l)
Q (m^3/s)
H (m)
V (m/s)
T (C)
125
3.237
0.44
0.82
14
Tabla7-42. Datos de entrada al modelo sin la aplicacin de tratamientos a los principales efluentes de la cuenca del Ro Katari.
Fuente: Elaboracin Propia.
SIMULACION DE O.D.
(Kd: M. Bansal y K2: O`Connor et al)
8
7
6
4
DBO (mg/L)
OD (mg/L)
5
3
2
1
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
O.D. Saturacin (mg/l)
OD SIMOD (mg/L)
Distancia Km
SIMULACION DBO5
(Kd: M. Bansal y K2: O`Connor et al)
132
120
108
96
84
72
60
48
36
24
12
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
Distancia en Km
TABLA DE RESULTADOS
PROYECTO:
PERIODO:
Escenario
Prediccin 3
Distancia
Acumulada(Km)
Temperatura
del Agua (C)
O.D.
SIMOD
(mg/l)
DBO
SIMOD
(mg/L)
O.D.
Saturacin
(mg/l)
K1
(1/da)
K2
(1/da)
O.D.
Mnimo
Permitido
(mg/L)
15.73
5.30
12.00
6.21
1.72
24.33
4.00
1.00
15.73
5.32
11.80
6.21
1.72
24.33
4.00
2.00
15.73
5.33
11.60
6.21
1.72
24.33
4.00
3.00
15.73
5.35
11.40
6.21
1.72
24.33
4.00
4.00
15.73
5.36
11.21
6.21
1.72
24.33
4.00
5.00
15.73
5.37
11.02
6.21
1.72
24.33
4.00
5.33
15.73
5.38
10.96
6.21
1.72
24.33
4.00
5.33
14.16
0.59
114.36
6.43
0.41
9.70
4.00
5.33
14.16
0.59
114.36
6.43
0.41
9.70
4.00
5.33
14.16
0.59
120.61
6.44
0.51
11.89
4.00
15.33
14.16
0.69
118.54
6.44
0.51
11.89
4.00
25.33
14.16
0.78
116.50
6.44
0.51
11.89
4.00
34.43
14.16
0.87
114.67
6.44
0.51
11.89
4.00
N
Segmento
Pgina | 162
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 163
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
La combinacin 11 en el modelo RIOSep (ver Tabla 7-6) fue la que dio mejores resultados.
La validacin del modelo mostro una tendencia similar a la de los datos observados,
confirmando as la combinacin de mtodos de la constante de reaireacin y desoxigenacin.
El modelo WASP 7.4 no calibro con los datos obtenidos en campo. El modelo necesita mayor
cantidad y calidad de informacin de entrada.
Los escenarios de prediccin se realizaron con el modelo SIMOD por dar mejores resultados.
Valor de correlacin cercano a 1 y menor porcentaje de error comparado a los otros dos
modelos.
Se plantearon tres tipos de escenarios, los cuales se obtuvieron de informacin recopilada y
experiencia por las campaas realizadas. En los escenarios de prediccin se tiene los
siguientes resultados:
- Escenario uno, el Ro Seco tiene valores de DBO 5 muy altos y OD muy bajos
clasificndolo en aguas de clase D segn el RMCH. Este ro recoge aguas pluviales,
residuales domsticas, industriales, lixiviados de basura, residuos slidos y hasta
hospitalarios. Pasa cerca de la planta de Puchukollo, absolutamente contaminado, donde se
descargan los efluentes de la planta depuradora. Para la realidad del Ro Katari y el Lago
menor, mientras no se traten todas las aguas residuales de El Alto la PTAR Puchukollo es un
esfuerzo que a pesar de su valor, no tiene significancia real en la problemtica de la Baha de
Cohana.
La PTAR Puchukollo no ha sido prevista ni tendra capacidad para interceptar y tratar las
aguas del Ro Seco, uno de los principales emisores de aguas contaminadas al Ro Pallina. Por
tal motivo se debe proponer una ampliacin para captar todas las aguas de la ciudad de El
Alto o construir otra planta de tratamiento.
- Escenario dos, en este escenario se supone tratamiento en las diferentes fuentes de descarga,
se proponen tratamientos primarios, secundarios y terciarios.
Aplicando un Tratamiento Primario disminuye la concentracin de la DBO 5 y mejora la
calidad del agua, llega a concentraciones entre 30 y 40 (mg/l) Clase D, mucho mejor
comparada a concentraciones de 80 (mg/l), no obstante se debe proponer otro tipos de
tratamiento para poder mejorar la calidad al menos un agua de Clase B. A pesar que el
Michael Russel Osina Torrez
Pgina | 164
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
oxigeno disuelto es afectado por las descargas, aumenta su concentracin en funcin del
tiempo y distancia recorrida, aun as las concentraciones son muy bajas, menores a 3 (mg/l).
(Ver grficos 7-40, 7-41)
Aplicando un Tratamiento Secundario la DBO 5 disminuye su concentracin clasificando la
calidad de agua Clase C y el oxigeno se encuentra muy cercano a 4 (mg/l) segn exige la
norma RMCH. Se puede considerar una concentracin aceptable (ver Grficos 7-42, 7-43).
Aplicando un Tratamiento Terciario la calidad del Ro Katari se encuentra dentro de la Clase
A, con concentraciones de DBO 5 menores a 5 (mg/l) y concentracin del OD de 5.50 (mg/l)
cumpliendo con la norma. Este tratamiento es el mejor de los tres propuestos, pero tiene la
desventaja de tener un gran costo, por lo que se ve como mejor solucin proponer un
tratamiento secundario a las descargas principales que aportan sus aguas a la cuenca del Ro
Katari, para mejorar la calidad de agua que se tiene a un agua de Clase B y no tener un costo
muy exagerado de inversin, el cual puede ser posiblemente financiado (ver Grficos 7-44, 745).
- Escenario tres, Se analiza la situacin sin la existencia de algn tratamiento en las descargas
a los ros principales de la cuenca del Ro Katari. La contaminacin en la Baha de Cohana
puede empeorar llegando a concentraciones muy elevadas de DBO 5 mayores a 100 (mg/l) y
oxigeno disuelto con concentraciones menores a 1 (mg/l) lo cual perjudicara al proceso de
degradacin de la materia contaminante (ver Grficos 7-46, 7-47) e incrementara la
eutrofizacin en el lugar de descarga.
Pgina | 165
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
9. CONCLUSIONES Y RECOMENDACIONES
9.1. CONCLUSIONES GENERALES
Los resultados del trabajo de campo y del modelado, fueron comparados con los
lmites permisibles del Reglamento en Materia de Contaminacin Hdrica (RMCH) de
la ley de Medio Ambiente N 1333, evidencindose que el Ro Pallina (afluente
principal del tramo de estudio), presenta en la mayora de las campaas,
concentraciones de oxigeno disuelto (OD) por debajo del valor mnimo establecido y a
su vez, elevados niveles de contaminacin por materia orgnica biodegradable
(DBO 5 ), con valores superiores al Limite Permisible correspondiente a la Clase D que
implica que la calidad del agua no es apta para consumo humano. La descarga del Ro
Katari a la Baha de Cohana vara de Clase C a Clase D, segn la poca del ao,
ocasionando indistintamente su eutrofizacin.
Ninguno de los tres municipios (El Alto, Viacha y Laja) que descargan sus aguas al
Ro Pallina, principal afluente del Ro Katari, no efectan el tratamiento de las aguas
residuales generadas en su mbito de influencia; tampoco, impiden el lanzamiento
indiscriminado de residuos slidos en su cauce o el de sus afluentes.
Los tres modelos SIMOD, RIOSep y WASP 7.4 que fueron utilizados, constituyen
herramientas tcnicas de gran utilidad para la simulacin del comportamiento de la
calidad de agua en trminos de DBO 5 y OD. Desde este punto de vista, son
consideradas herramientas apropiadas de prediccin, planificacin, anlisis de riesgos,
validacin y procesamiento de informacin de calidad del agua. Su uso eficiente
depender del escenario y las condiciones del medio en las que se apliquen.
Pgina | 166
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Las relaciones matemticas que permiten calcular con mayor precisin los valores
ejercidos en el Ro Katari, son el mtodo de OConnor et al para la constante de
reaireacin (K 2 ) y el mtodo de M. Bansal para la constante de desoxigenacin (K 1 ).
Pgina | 167
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Constante de reaireacin
K 1 (1/da)
K 2 (1/da)
SIMOD
0.12 1.95
0.65 30.00
0.09 1.03
0.71 22.15
EPA
0.00 2.00
0.00 100.00
Para el caso presente no result til el empleo del modelo WASP 7.4 para simular la
calidad del agua en trminos de DBO 5 y OD. Si bien el modelo ofrece mayor cantidad
de parmetros a simular, necesita excesiva cantidad y calidad de informacin de
entrada, eso lo convierte en un modelo poco aplicable en nuestro medio, debido a la
carencia de informacin existente para nuestros ros.
Pgina | 168
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
9.3. RECOMENDACIONES
Inventariar las industrias para determinar con exactitud el efecto o grado contaminante
del sector industrial.
El servicio de Aseo Urbano debe mejorar la limpieza sobre los sumideros y bocas de
tormentas, que se ubican en diferentes puntos del rea urbana de El Alto, Viacha y
Laja.
Pgina | 169
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 170
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
10.
BIBLIOGRAFA Y REFERENCIAS
Pgina | 171
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
15. Fontrbel Rada, Francisco. (n.d.): Evaluacin preliminar de la calidad hdrica, mediante
indicadores fisicoqumicos y biolgicos, en la Baha de Cohana, Lago Titicaca del
departamento de La Paz, Bolivia. La Paz, Bolivia, Centro de Estudios de Postgrado e
Investigacin (CEPI) de la Universidad Mayor Real y Pontificia de San Francisco
Xavier de Chuquisaca, Sucre (Bolivia).
16. Humerez, Susana E.: Evaluacin de la ecotoxicidad de las aguas del Ro Pallina, Tesis
de Ingeniera Ambiental, Universidad Mayor de San Andres, La Paz, Bolivia. 2011.
17. Husnain Haider and Waris Ali: Effect of Wastewater Treatment on Bio-kinetics of
Dissolved Oxygen in River Ravi, (6), 2010.
18. IBTEN: Datos de campaas, 2007-2011
19. Lam, E.; Demergasso, C.; Galleguillos, P.; Escudero, L.: Diseo de un Biorreactor
Anaerbico para el tratamiento de aguas residuales en la produccin de agua potable,
XXVII Congresso Interamericano de Engenharia Sanitria e Ambiental, Porto Alegre,
2000.
20. Llavador Colomer Fernando: Modelos Matemticos De Sistemas Variables En El
Tiempo, Entidad Pblica de Saneamiento de Aguas Residuales de la Comunidad
Valenciana, 2008.
21. Manhattan College: Quality models of natural water systems, Course notes for the
Summer Institute in Water Pollution Control, Manhattan College, Bronx, NY, 1983.
22. Mattos, S: Evaluacin del Comportamiento de modelos matemticos de calidad de agua,
Tesis de Ingeniera Civil, Universidad Mayor de San Andres, La Paz, Bolivia. 2003.
23. Menendez Gutierrez Carlos y. Prez Olmo Jess M: Procesos para el tratamiento de
aguas residuales, 2007, 3-13.
24. Miller J. C.: Estadstica para qumica analtica, (2), 1993, 87-114
25. Ministerio del Medio Ambiente. 1997. Reglamento en Materia de Contaminacin
Hdrica. Articulo 4.
26. Metcalf&Eddy. Ingeniera de aguas residuales; tratamiento, vertido y reutilizacin,
Tercera Edicin, Ed. McGraw-Hill, Inc., Espaa, 1995.
27. Monerris Miguel Martn y Marzal Domnech Paula: Modelacin de la Calidad del Agua.
Universidad Politcnica de Valencia, Departamento de Ingeniera Hidrulica y Medio
Ambiente, Escuela Tcnica Superior de Ingenieros de Caminos, Canales y Puertos.
28. M.Sc. Chuco,Tolessa Deksissa: Dynamic integrated modeling of basic water quality and
fate and effect of organic contaminants in rivers, 2004.
Pgina | 172
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 173
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
Pgina | 174
EVALUACIN DE LA CALIDAD DE LAS AGUAS DEL RO KATARI, LA PAZ, BOLIVIA, MEDIANTE UN MODELO MATEMTICO
11.
ANEXOS
ANEXO A
FOTOGRAFAS DE LOS PUNTOS DE MUESTREO, TOMA DE MUESTRA,
MEDICIONES IN SITU, EQUIPOS DE CAMPO Y MEDICIONES EN LABORATORIO
ANEXO B
DATOS OBTENIDOS EN EL TRABAJO DE CAMPO
ANEXO C
PLAN DE MONITOREO
ANEXO D
CORRIDAS DE LOS MODELOS
ANEXO E
MATRIZ DE IDENTIFICACIN DE IMPACTOS AMBIENTALES (M1)
ANEXO F
CALCULO DE LA K1 CON DATOS DE LABORATORIO PARA LA OCTAVA
CAMPAA
ANEXO G
RESULTADOS DE ANALISIS QUIMICO DE LA OCTAVA CAMPAA
Pgina | 175