UNIVERSIDAD DE

GUADALAJARA

DEPARTAMENTO DE FÍSICA

MODELACIÓN DE LA CONTAMINACIÓN
EN RÍOS

PROYECTO INTEGRADOR ETAPA III

Equipo: Happy Ganesh Chaturthi

José Alfredo Quiñonez Sandoval
Edilson Josafat Rodríguez Delgadillo
Daniel Zanella Alvarado
1. INTRODUCCIÓN
Desde los inicios de la humanidad hasta los días actuales, los humanos tenemos una
gran dependencia con el agua, dado que los sedentarios hacían sus campamentos cerca
de grandes reservas de agua, para poder realizar actividades como la agricultura y
el ganado, además es utilizada como fuente de energía, gracias a la gran fuerza y
potencia que esta posee, siendo de gran importancia para la existencia de la vida
misma.
Con la llegada de las fábricas, parte de los ríos y lagos se han visto afectados por
los desechos industriales que estos producen y desechan en dichos lugares, además
de los desechos que los humanos producimos, un ejemplo de esto serían los desechos
sanitarios, generando que estas aguas ya no sean seguras no solo para los humanos,
sino para los ecosistemas que se verán afectados por el paso de estas aguas, dado que
son su principal abastecimiento de agua.
Dado que el análisis del comportamiento de los contaminantes suele ser complicado
por los distintos factores que están implicados, podemos analizar este problema con
modelados matemáticos, pues puede ser de mayor utilidad al manipular las distintas
variables que están implicadas en el comportamiento, y de esta manera tener una
solución a dicho problema. Podemos hacer una interpretación del transporte de los
contaminantes con un modelo matemático, siendo expresado en EDP’s, podemos
generar un análisis más fácil proponiendo algunas condiciones iniciales y de frontera.

2. OBJETIVOS
En el siguiente documento se buscará abordar de una manera clara y concisa las pro-
blemáticas que se producen por la contaminación de las aguas de ríos y lagos, además
del contaminante a tratar y su fuente. Así mismo, exponer las consecuencias que este
tiene para el ser humano y el ecosistema, además de presentar algunos procesos para
mitigar su impacto. Se hará uso de un modelado matemático para el análisis de este
problema y una comparación con los casos reales, donde utilizaremos algunas ex-
presiones matemáticas que simulen los procesos de transporte de los contaminantes,
y de esta manera generar un problema matemático con condiciones iniciales y de
frontera. Se establecerá un área de estudio,es decir, una sección en particular del río,
así como también la localización de las distintas fuentes de polución y la velocidad
del río.

1
3. METODOLOGÍA
Durante el desarrollo del proyecto se recaudó información de diversas fuentes biblio-
gráficas y artículos acerca de las causas de la contaminación del agua más específica-
mente de ríos y lagos, así como su modelación matemática asociada y las condiciones
iniciales y de frontera.
Si tomamos un río como objeto de estudio, debemos establecer la velocidad del
mismo, tomando en cuenta la velocidad inicial, lo que genera diversas preguntas
como la dirección del río, la contaminación del mismo en ese punto, entre otras
cosas.
Si un material es desechado al agua, existen varios procesos que están inmiscuidos en
su traslado, así como su mezcla en la misma. Es decir, si el río tiene movimiento, esto
generará en el desecho un traslado a través de distintos puntos del río. En estos casos,
existen tres tipos de transporte: la advección, que es el trasporte que tiene el desecho
por la corriente del agua, la convección, el cual es el transporte que se genera por la
inestabilidad hidrostática, y por último, el corte, que es el transporte que depende de
variaciones espaciales de la velocidad del agua, la cual sufre de cambios dependiendo
de la profundidad.
Si hablamos únicamente del transporte por advección, podemos decir que es el más
importante por la naturaleza del ciclo del agua y de los procesos asociados que este
tiene, con la circulación entre los diferentes depósitos de la hidrosfera, el agua puede
ejercer este transporte por tres vías distintas: arrastre, suspensión y disolución. [1]
Teniendo esto, podemos proponer el modelo matemático con sus condiciones iniciales
o de frontera.
Por lo antes mencionado, la propuesta para el modelo puede ser de manera general,
el cual es el siguiente:
 
∂C ∂C ∂ ∂C
+u − µ + κC = S(x, t) (1)
∂t ∂x ∂x ∂x

Donde:
∂C
∂t
Es el cambio en la concentración de contaminantes con el tiempo.

u ∂C
∂x
Es el termino que gobierna el proceso de advección en el rió.
∂
µ ∂C


− ∂x ∂x
Es el termino que gobierna la difusión.

κC Es el termino que gobierna la reacción de los contaminantes en el río.

2
 S(x, t) Son las fuentes externas de contaminante, las cuales pueden ser: constan-
tes, depender de posición y tiempo, o depender de la presencia de contaminantes
o no estar presentes directamente, si no en las condiciones de frontera.

Las variables, u, µ y κ pueden ser valores constantes o funciones en el caso de

µ depende del tipo de contaminante a analizar.

4. PLANTEAMIENTO DEL PROBLEMA

Nuestro objeto de estudio es el Río Santiago, recorre parte de seis estados de la
República Mexicana: Aguascalientes, Durango, Guanajuato, Jalisco, Nayarit y Zaca-
tecas, dicho río es uno de los más contaminados del país, dado que ahí terminan los
desechos tóxicos de diversas fábricas que se encuentran en nuestro estado, además de
las aguas negras de la ZMG y las diversas poblaciones que se encuentran a pie del
río.
El grave deterioro ambiental de la cuenca del río Santiago -que para algunos se inicia
en los años setenta – sufre el consecuente pasmo de las autoridades públicas esta-
tales, la indiferencia de las federales electas en el año 2006 y la impotencia de las
autoridades municipales. El panorama resulta no sólo desolador en cuanto al cuida-
do y preservación de la riqueza natural, de la vida vegetal y animal que anima este
ecosistema, sino ominoso en cuanto se le relaciona con la población que habita en
sus márgenes o inmediaciones.
Habiendo hecho una breve descripción de este río, comenzaremos con el análisis de
concentración de contaminantes en el río Santiago. Tomamos una parte del río para
facilitar el estudio, Nuestro dominio comenzará en el inicio del río, justo en el lago
de Chapala y culminará unos cuantos metros después del municipio de El Salto, pa-
sando por los municipios de Ocotlán, Poncitlán, Juanacatlán y El Salto, hay algunos
otros municipios, los cuales no tomaremos en cuenta dado que son comunidades muy
pequeñas. Estos municipios serán nuestras fuentes de polución, dado que su pobla-
ción es considerable, significa que existe una alta cantidad de desechos vertidos al río
en esos puntos.
Para determinar el dominio de nuestra área de estudio se hace uso de Google Earth
y junto a su herramienta de medición entre dos puntos se fue trazando la trayectoria
que sigue el río para poder determinar a que distancia se encontraba el inicio (Lago
de Chapala) de los distintos municipios, hasta llegar al final del dominio.

3
Figura 1: Vista de satélite de los lugares de descarga haciendo uso de Google Earth.

Una vez seleccionandos los lugares de descarga, calculamos la distancia a la cual se

encuentran dichos puntos, teniendo como resultado la siguiente tabla:

Localizaciones Dominio
Chapala 0 km
Ocotlán 3.5 km
Poncitlán 22.63 km
El Salto 66.46 km
Juanacatlán 66.46 km
Fin del dominio 69 km

Cuadro 1: Dominio y localización de las fuentes de polución.

Lo siguiente es establecer la velocidad del Río, para esto haremos uso de la fórmula
de Manning:
1 2 1
U= R3 S 2 (2)
n

4
Donde U es la velocidad del río, n el coeficiente de rugosidad, R la anchura del río y
S la pendiente del mismo.
Para eso, nuevamente haremos uso de la herramienta Google Earths. Se mide la
altitud del Río aproximadamente cada kilometro y lo mismo con la anchura, con el
fin de promediar estas cantidades y dividir el río, es decir, tendremos cuatro secciones.

Figura 2: Altitud.

Figura 3: Anchura.

5
 Región [km] Pendiente Ancho [m]
0 a 3.5 0.00102 30
3.5 a 22.63 0.000105 51.05
22.63 a 66.46 0.000326 62.04
66.46 a 69 0.0165 30

Cuadro 2: Pendientes y anchuras.

Utilizando como referencia el libro "Water Quality Modelling for Rivers and Streams"de
Benedini y Tsakiris, tomando los datos de la tabla de la página 43, extraemos el co-
eficiente de rugosidad.[4]

Figura 4: Coeficientes de rugosidad.

Suponemos que la superficie es rocosa e irregular en todo el río, siendo una superficie
en buen estado, por lo que ya tenemos nuestro coeficiente de rugosidad.
Tenemos todos los datos necesarios para poder obtener la velocidad con la ecuación
2, de este modo:

Región [km] Velocidad [km/h]

0 a 3.5 27.7516
3.5 a 22.63 12.6910
22.63 a 66.46 25.4628
66.46 a 69 111.6172

Cuadro 3: Pendientes y anchuras.

6
Nuestro coeficiente κ, según la referencia [4], está en el intervalo de 0.2 a 0.5 por día,
por lo que tomaremos un valor constante de 0.2, esto para todo el dominio.

4.1. Condiciones iniciales y de Frontera

Nuestra condición inicial es que el río comienza limpio por tanto:

C(x, t) = C(x, 0) = 0 (3)

Para las condiciones de frontera, consideramos que nuestro río comienza limpio, por
lo que el contaminante vendrá por las descargas que sucederán a medida que el río
fluya por diversas regiones y para el final del área de estudio vamos a considerar
que la contaminación seguirá fluyendo, lo que quiere decir los contaminantes no se
quedarán estancados en esa zona, dicho esto, vamos a aplicar las condiciones tipo
Neumann:
∂
C(0, t) = 0 (4)
∂x
∂
C(69km, t) = 0 (5)
∂x
Sea ec.4 la condición inicial y ec.5 la condición final.

Para calcular la magnitud de los desechos tomaremos únicamente los residuos fecales
de las regiones que colindan con el río Santiago, que de acuerdo a [4], en promedio,
cada persona produce 70 g/día de desecho, por lo que debemos calcular los desechos
que generan las poblaciones en total, para esto utilizaremos datos del Instituto de
Información Estadística y Geográfica de Jalisco (IIEG), los datos que nos propor-
cionan son del año 2020, por lo que quizá estén un poco alejados de la realidad,
actualmente.

Poblaciones de estudio en el año 2020

Municipio Habitantes
Ocotlán 106,050
Poncitlán 53,659
El Salto 232,852
Juanacatlán 30,855

Cuadro 4: Poblaciones en 2020

7
Solamente nos centraremos en las localidades que se encuentra pegadas al río, dado
que es muy probable que sus habitantes terminen desechando sus residuos en el río.
Por lo tanto, nos centraremos en las localidades de Ocotlán, Poncitlán, El Salto y
Juanacatlán, cuyas poblaciones son las más grandes y se encuentran al pie del río.
La distribución de la población de los municipios según el Instituto de Información
Estadística y Geográfica de Jalisco (IIEG) son las siguientes:

Figura 5: Población de Ocotlán 2020.

Figura 6: Población de Poncitlán 2020.

8
 Figura 7: Población de El Salto 2020.

Figura 8: Población de Juanacatlán 2020.

9
Por lo antes mencionado, tomamos solo los datos que necesitamos, teniendo como
resultado:

Poblaciones de interés en el año 2020

Localidad Habitantes
Ocotlán 94,978
Poncitlán 14,636
El Salto 42,576
Juanacatlán 9,626

Cuadro 5: Poblaciones en las localidades en 2020

De esta forma, el total de habitantes que generan desechos que terminan en el río es
de 161,816 habitantes, lo cual multiplicamos por el desecho en promedio que genera
una persona por día, de este modo:
 
kg
ST OT AL = 161, 816 × 70 g/día = 11, 327.12 kg/día (6)
1000g

La función que determina las descargas S(x,t) está propuesta para que exista una
descarga en cada punto de polución, con una forma de distribución normal, siendo que
mientras más cerca se encuentre del centro de cada localidad habrá mayor descarga
y mientras más alejado de ese sitio habrá menor descarga de desecho, dicha función
es la siguiente:
 
1 2 πt  − 2k1 (x−a)2 − 1 (x−b)2 − 1 (x−c)2

S(x, t) = sin s1 e 1 + s2 e 2k2 + s3 e 2k3 (7)
2 24

Donde a, b y c son las posiciones de las poblaciones, los valores ki son la mitad de la
distancia que el río pasa por la población y si son las descargas en cada población.

4.2. Esquema numérico

Tomando en cuenta las condiciones tanto iniciales como de frontera, se ha llegado a
un esquema numérico explicito de la ecuación general de transporte de contaminantes
en un río(vease ec. 1). Ahora bien, simplificando la expresión nos queda que:

Cin+1 = αCin + βCi−1

n n
+ γCi+1 + Sin ∆t (8)

10
Donde:
α = 1 − κ∆t − 2β + λ

∆t ∆t
β= , λ = u , γ = β − λ.
∆x2 ∆x

5. SOLUCIÓN DEL PROBLEMA

Figura 9: Resultado de la simulación

11
 Figura 10: Resultado de la simulación

Para obtener el modelo matemático se utilizó la ecuación 7 de donde se obtuvieron

los contaminantes por hora de cada municipio, de esta manera, introdujimos estos
datos en una matriz S, así mismo utilizamos una velocidad promedio como se men-
cionó con anterioridad para facilitar los calculos númericos.
De nuestra simulación se puede observar que la incersión de contaminantes es cons-
tante la mayoría del tiempo, no obstante en la región de segmentos alrededor de los
50 km, x=22, en la madrugada hay una gran cantidad de ingresos en contaminantes.
Se observa que esta es una simulación de los ingresos de los contaminantes respecto
al tiempo río abajo, por último también se observa que el nivel de deshechos es ma-
yor en el extremo del río, esto coincide con los poblados mas grandes y con mayor
cercanía al río.

12
6. CONCLUSIONES
En esta etapa, seleccionamos nuestra área de estudio (dominio), el cual comienza en
el lago de Chapala y culmina unos cuantos metros después del municipio de El Sal-
to, además calculamos los desechos que generan las comunidades que se encuentran
cerca del río, aunque esto es solamente un aproximado.
Teniendo el modelo matemático, pudimos calcular la contaminación del río Santiago,
aunque es cierto que utilizamos algunos datos que nos favorecieron a simplicar los
cálculos, así como también los datos de las poblaciones, pues los datos que encon-
tramos no estaban actualizados al año en curso, pero al ser un proyecto estudiantil
no tiene tanta importancia, cosa que estaría mal si se tratara de algo más profesio-
nal, pues en esos casos si son necesario datos más precisos para poder obtener una
aproximación más precisa y de acorde a lo que pasa en la realidad, para obtener
una mejor simulación habría que ir a tomar datos al área de estudio ya que al usar
una herramienta virtual pudo haber mucho errores de medición aunque esto llevaría
muchísimo tiempo por la longitud del Río, sin embargo nos acercaría mucho a un
modelo más realista del mismo.
Los problemas de contaminación de las reservas de aguas, como en este caso el río
Santiago son un tema de gran importancia, por lo que este tipo de proyectos pueden
ser de gran ayudar para encontrar un solución a esto, ya que podemos identificar los
lugares donde se generan la mayor cantidad de desechos y con esto encontrar una
solución que sea eficiente para poder limpiar un poco el río.

7. BIBLIOGRAFÍA Y REFERENCIAS
[1] Universidad De las Palmas De Gran Canaria. Transporte y difusión de contami-
nantes en las aguas naturales. url:https://www2.ulpgc.es/hege/almacen/ downloa-
d/42/42250/tema,v1.pdf.

[2] Jairo Escobar. La contaminación de los ríos y sus efectos en las áreas costeras y
el mar. CEPAL, 2002.

[3] MA Jurado and ID Mercado. Revisión sistemática de técnicas no convencionales

para la evaluación de la calidad del agua de ríos contaminados con plaguicidas. Entre
ciencia e ingeniería, 11(21):56–65, 2017.

[4] Marcello Benedini and George Tsakiris. Water quality modelling for rivers and
streams. Springer Science and Business Media, 2013.

13
14