CI4102 - Ingenierı́a Ambiental Prof.

Ana Lucı́a Prieto

Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Pauta Taller 2
P1 Reactores

En el norte de Chile el mercurio es un compuesto comúnmente usado en minerı́a para solubilizar o

separar ciertos metales. Los efluentes producidos en el proceso son altamente corrosivos y difı́ciles
de tratar. La disposición de estos efluentes en pozos de retención o almacenamiento es una práctica
común para el manejo de ellos. En el fondo de uno de estos pozos se acumulan grandes cantidades
de sulfuro de mercurio (HgS), que se mantiene en equilibrio con el agua en contacto con el fondo
del pozo.

HgS ↔ Hg 2+ + S 2− Ksp = 4 · 10−23 M 2

Algunos microorganismos que proliferan en estos pozos consumen los iones de mercurio que cons-
tantemente se disuelven en el agua, manteniendo una concentración estable de 10−5 M .

S 2− + 2H + ↔ H2 S K = 1021,1 M −2
A su vez, el S 2− reacciona con los protones en el agua para formar ácido sulfúrico gaseoso (H2 S).
El área superficial del pozo es de 10 m2 y la profundidad es de 2 m.

a) Elabore un balance de masa para el H2 S en el pozo y derive una ecuación que describa el cambio
de la concentración de H2 S con la profundidad (Representa la concentración de H2 S con C y
la profundidad con z).
Hint: Considere una contribución de azufre desde el fondo que se distribuye en la columna de
agua por flujo difusivo, el cual se representa como Fz = −D · dCdz , donde D es la constante de
2
difusión (m /s).

Solución:

Balance de masa para H2 S asumiendo estado permanente:

dM
= In − Out =⇒ 0 = Fz · A − Fz+∆z · A (1)
dt
Dividiendo por V = ∆Z · A y tomando el lı́mite cuando ∆z ⇒ 0:

dFz
− =0 (2)
dz
Utilizando el hint:

d2 C dC
=⇒ D 2
= 0 =⇒ =0 (3)
dz dz
E integrando:

Departamento de Ingenierı́a Civil 1 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

∴ C(z) = A · z + B (4)

b) Derive una expresión que describa la concentración de H2 S como función de la profundidad

(C(z)) en términos de [H + ].
Hint: tenga en cuenta las condiciones de borde para su cálculo.

Solución:

Evaluando en las condiciones de borde:

Para z = 2 m se tiene que [H2 S]=0

Para z = 0 se tiene que [H2 S] depende del equilibrio quı́mico del sistema.

El sedimento de HgS contribuye a S 2− al agua que luego forma H2 S.

Por equilibrio quı́mico:

[S 2− ] = Ksp [Hg 2+ ] (5)

[H2 S] = K[S 2− ][H + ]2 (6)

Reemplazando (5) en (6) y usando los valores de Ksp y K, se tiene que:

C(0) = 5, 04 · 103 M −1 [H + ]2 (7)

Usando los valores de C(0) y C(2) en (4):

z
C(z) = 5, 04 · 103 M −1 [H + ]2 (1 − ) (8)
2
c) Si se sabe que el pH del agua es cercano a 2, ¿Cuál es la concentración de H2 S en la mitad del
pozo (z=1 m)?

Solución:

Para z=1 y [H + ]=10 −2 M, se tiene que C(z=1[m]) = 0,25 M

d) Si asumiera por un momento condiciones de mezcla completa en el pozo, ¿Cuál es la concentra-

ción de H2 S en el aire justo sobre la superficie? Exprese la concentración en ppm (KH =0,115
M/atm)

Departamento de Ingenierı́a Civil 2 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Solución:
Mezcla completa permite asumir condiciones de concentración uniformes en el tanque.

0, 25M
CH2 S,aq = KH · PH2 S =⇒ PH2 S = = 2, 19[atm] ≈ 2, 2[atm] (9)
0, 115M/atm
Expresando en ppm y asumiendo la presión atmosférica igual a 1 [atm]:

Cppm
2, 2 = =⇒ Cppm = 2, 2 · 106 [ppm] (10)
106

Departamento de Ingenierı́a Civil 3 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

P2 Limnologı́a
Usted se encuentra estudiando un lago de la zona sur de Chile. Un esquema de éste se presenta en
la figura 1.

Figura 1: Esquema lago con flujo subterráneo y descarga de contaminación.

El lago se encuentra estratificado por temperatura, separando epilimnio e hipolimnio cada uno
con volúmenes Vepi y Vhipo . Del hipolimnio escurre agua con Qesc hacia un acuı́fero subterráneo de
volumen desconocido. De este acuı́fero, un 90 % del flujo que ingresa sale nuevamente a la superficie,
para juntarse con un rı́o de caudal Qnat e ingresar al lago hacia el epilimnio. El rı́o trae una carga
de materia orgánica Cnat .
La materia orgánica, se degrada SOLO en el lago a una tasa kepi en el epilimnio y khipo en el
hipolimnio. Además, en el epilimnio y en el hipolimnio, la materia orgánica presenta una tasa de
sedimentación ksed . Por otro lado, la materia orgánica en el hipolimnio es capaz de resuspenderse
al epilimnio debido a turbulencia a una tasa kres .
Del lago sale un caudal Qef . Una empresa aguas abajo, extrae un 0,5 % del caudal del rı́o, y lo
devuelve al lago con una carga de Ci .

a) Si la materia orgánica NO presenta degradación en las etapas subterráneas (escurrimiento,

acuı́fero, ascensión). Encuentre las concentraciones de la materia orgánica en el epilimnio e
hipolimnio. Hint: Plantee los balances de masa en el punto de mezcla, en el epilimnio y en el
hipolimnio. Asuma estado permanente.

b) Compare su resultado con los valores estipulados en el decreto 1, adjunto con el enunciado.
¿Cómo los supuestos tomados anteriormente podrı́an afectar la diferencia entre el valor calculado
y el valor del decreto?

c) Si la materia orgánica SÍ presenta degradación en las etapas subterráneas (escurrimiento, acuı́fe-
ro, ascensión) y desde el punto de mezcla a la entrada del lago. Plantee un modelo en estado
permanente que permita encontrar la concentración de materia orgánica en el acuı́fero. Para

Departamento de Ingenierı́a Civil 4 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

este modelo puede incluir los parámetros: Vhipo−acuı́f ero : volumen desde el hipolimnio al acuı́fe-
ro; Vacuı́f ero : volumen del acuı́fero; Vacuı́f ero−superf icie : volumen desde el acuı́fero a la superficie;
ksub : constante de decaimiento de la materia orgánica de primer orden; Vmezcla−lago : volumen
desde el punto de mezcla al lago.

Recuerde mencionar TODOS LOS SUPUESTOS que estime convenientes para su modelo.

Departamento de Ingenierı́a Civil 5 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Datos:

Parámetro Valor Unidad

Qnat 8 m3 /s
Qesc 0,5 m3 /s
Cnat 20 mg/L
Ci 80 mg/L
Vepi 109 m3
Vhipo 107 m3
kepi 0,06 1/dı́a
khipo 0,03 1/dı́a
kres 0,0002 1/dı́a
ksed 0,0015 1/dı́a

Departamento de Ingenierı́a Civil 6 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Solución:

a) En primer lugar se plantea el balance de volumen en el lago con el fin de obtener el valor del
caudal de salida Qef :

dVlago
= Qnat + Qi + 0, 9Qesc − Qef − Qesc
dt
Bajo el supuesto de estado permanente, se tiene que el lago conserva su volumen, lo que en
ecuaciones se traduce a lo siguiente:

Qnat + Qi + 0, 9Qesc = Qef + Qesc

Qi = 0, 005Qef
Qef = 690332[m3 /dı́a]
Qi = 3451[m3 /dı́a]
Siguiendo las indicaciones, se realizan los balances de masa en los distintos sectores:
Balance punto de mezcla Se comienza con el balance global:
dM
= Qnat Cnat + 0, 9Qesc Chipo − QM CM
dt
Asumiendo estado permanente (puntos en donde se juntan caudales como RMC):

Qnat Cnat + 0, 9Qesc Chipo = QM CM

Balance en epilimnio Se comienza con el balance global:

dMepi
= QM CM + Qi Ci + kres Chipo Vhipo − ksed Cepi Vepi − kepi Cepi Vepi − Qef Cepi
dt
Analizando en estado permanente (sin acumulación), se llega a lo siguiente:

0 = QM CM + Qi Ci + kres Chipo Vhipo − ksed Cepi Vepi − kepi Cepi Vepi − Qef Cepi

Balance hipolimnio Se comienza con el balance global:

dMhipo
= ksed Cepi Vepi − kres Chipo Vhipo − ksed Chipo Vhipo − khipo Chipo Vhipo − Qesc Chipo
dt
Analizando en estado permanente (sin acumulación), se llega a lo siguiente:

0 = ksed Cepi Vepi − kres Chipo Vhipo − ksed Chipo Vhipo − khipo Chipo Vhipo − Qesc Chipo

Del balance en el hipolimnio, se encuentra la concentración en éste en función de la concentración

Departamento de Ingenierı́a Civil 7 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

en el epilimnio y constantes conocidas:

ksed Cepi Vepi
Chipo =
kres Vhipo + khipo Vhipo + Qes + ksed Vhipo

Uniendo las ecuaciones, se obtiene:

ksed Cepi Vepi
0 = Qnat Cnat + 0, 9Qesc + Qi Ci − kepi Cepi Vepi − Qef Cepi
kres Vhipo + khipo Vhipo + Qes + ksed Vhipo
ksed Cepi Vepi
+kres Vhipo − ksed Cepi Vepi
kres Vhipo + khipo Vhipo + Qes + ksed Vhipo

Despejando Cepi y reemplazando, se obtiene:

Cepi = 0, 227[mg/L]

b) Sı́ se cumple la norma.

¿Cómo influyen los supuestos?

No se está considerando degradación en las etapas subterráneas, lo cual, si bien podrı́a

significar la disminución de materia orgánica, puede dar paso a la generación de otros
contaminantes.
El supuesto de mezcla perfecta podrı́a no modelar de manera adecuada el sector en donde
está la descarga industrial, lo cual podrı́a significar que hayan zonas muy contaminadas, lo
cual podrı́a matar organismos y contaminar en cadena el resto del lago.
No se considera la resuspensión de materia orgánica desde el fondo del lago, lo cual podrı́a
ser una fuente importante de ésta.
Al hacer balance de volumen en el epilimnio, se podrı́a observar que éste va aumentando
su tamaño con el tiempo, esto podrı́a alterar los resultados.

Se pueden nombrar muchas implicancias. Se requieren al menos 2, debido a lo conversado en las

clases auxiliares y los balances de masa.

c) El balance en el epilimnio y en el hipolimnio se mantienen. La concentración a la entrada del

acuı́fero se puede modelar como un PFR:
Z Cin,acuif ero
Vhipo−acuif ero dC
=
Qesc Chipo −ksub C

Luego:
Cin,acuı́f ero=Chipo ·exp(−ksub ·θH1 )
El acuı́fero se puede modelar como un reactor de mezcla perfecta.
Qesc Cin,acuif ero
Cout,acuif ero =
Qout + ksub Vacuif ero

Departamento de Ingenierı́a Civil 8 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

En el punto de mezcla, el agua subterránea vendrá con una concentración que se rige por:

Csup = Cout,acuif ero · exp(−ksub θH2 )

De esta manera, en el punto de mezcla:

Qnat Cnat + 0, 9Qesc Csup = QM CM

Hay que notar que el flujo de entrada y salida del acuı́fero no son iguales, por lo que el volumen
del acuı́fero y la concentración de materia orgánica en el mismo serán dependientes del tiempo.

Departamento de Ingenierı́a Civil 9 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

P3 Ecologı́a

Cuando se tienen una serie de datos de DBO y se quiere calcular el valor L0 , se puede seguir el
siguiente método:
Realizar un gráfico con DBOt en el eje x, y DBOt+1 en el eje y.

Ajustar los datos con una regresión lineal.

Intersectar la curva ajustada con el curva DBOt+1 =DBOt (y=x).

El punto de la intersección de estas curvas se interpreta como el punto en donde la DBOt+1 es igual
a DBOt , es decir, que ya no hay variación en el valor del oxı́geno consumido, en otras palabras L0 .
Esto se ilustra en la figura 2.
Ejemplo de obtención de L0 = 192 [mg/L]
300

250

200
DBOt+1 [mg/L]

150

100

50

0
0 50 100 150 200 250 300
DBOt [mg/L]

Figura 2: Gráfico de metodologı́a descrita.

Usted ha encontrado los siguientes valores medidos de DBO para la muestra de una descarga
doméstica:

t [dı́a] DBO5 [mg/L]

0,5 28
1 53
1,5 74
2 92
2,5 109
3 125
3,5 137
4 149
4,5 157
5 167

a) Determine la concentración inicial de DBO (L0 ). ¿Esta es carbonatada o nitrogenada?. Argu-

mente.

Departamento de Ingenierı́a Civil 10 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

b) Como usted sabe, la demanda biológica de oxı́geno nitrogenada NDBO es posible observarla
aproximadamente después del dı́a 5 en una muestra. Una ecuación que puede modelar esta
demanda está dada por:

N DBO(t) = N DBOtotal · (1 − exp(−0, 22(t)))

Sabiendo que el contenido en nitrógeno de Kjeldahl es de 27 [mg/L]. Grafique la curvas CDBO,

NDBO y DBO para un horizonte de 30 dı́as.

c) Cree una función que represente la fracción CDBO/DBOtotal y grafı́quela en el mismo horizonte
anterior.

Departamento de Ingenierı́a Civil 11 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Solución

a) Gráficamente:

Ejemplo de obtención de L0 = 231 [mg/L]

300

250

200
DBOt+1 [mg/L]

150

100

50

0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300
DBOt [mg/L]

Figura 3: Gráfico de metodologı́a descrita.

El valor obtenido 231 [mg/L] permite calcular la constante de reacción según:

k = ln(1 − DBO5 /L0 )/ − 5

Esto permite plantear l modelo completo.

NB: Para utilizar la ecuación anterior, se está suponiendo que la curva de la DBO pasa exacta-
mente por el valor de la DBO5 en el dı́a 5. Para ser más exhaustivos se debe hacer un ajuste de
parámetros. De todas maneras, incluso haciendo un ajuste de parámetros habrá un sesgo debido
a la toma experimental de los datos.

Departamento de Ingenierı́a Civil 12 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

b) Con el modelo anterior, y teniendo en cuenta que: DBOT = NDBO + CDBO, se tiene la
siguiente salida gráfica.

Demanda de oxígeno en agua residual

400
CDBO
350 NDBO
DBOt

300

250
DBO [mg/L]

200

150

100

50

0
0 5 10 15 20 25 30
Tiempo [días]

Figura 4: Sistema DBO

c) La fracción se construye según:

CDBO
f=
DBOT

Departamento de Ingenierı́a Civil 13 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Cuya solución gráfica está dada por:

Fracción de DBO
1

0.95

0.9

0.85
Fracción

0.8

0.75

0.7

0.65
0 5 10 15 20 25 30
Tiempo [días]

Figura 5: Fracción

Departamento de Ingenierı́a Civil 14 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

P4 Modelo Streeter-Phelps

Una industria desea descargar sus residuos a un rı́o cercano. Como sus residuos poseen una DBO5
más alta que la permitida (máx 35 mg/l), ésta debe tratar sus efluentes. El tratamiento consiste en
un reactor de mezcla completa abierto a la atmósfera (laguna de estabilización) donde se consume
materia orgánica y además se produce una reaireación de la descarga. Antes de alcanzar la laguna,
el agua es conducida por un ducto CERRADO a la atmósfera. La laguna desagua a un canal (por
definición, abierto a la atmósfera), que finalmente descarga en un rı́o. En la Figura 6 se muestra
un esquema de la situaciónn y en la Tabla 1 los datos necesarios para resolver el problema.

Figura 6: Esquema problema 4.

Parámetro Simbolo Magnitud Unidad

Caudal Industria Q[ ind] 70 [m3 /s]
DBO Industria L[ ind] 50 [mg/L]
OD Industria OD[ ind] 3 [mg/L]
Área transversal del ducto AD 1 [m2 ]
Largo del ducto DD 500 [m]
Volumen de la Laguna V 5000 [m3 ]
Área transversal del canal AC 2 [m2 ]
Largo del canal LC 600 [m]
Tasa de decaimiento kd 0.3 [1/dı́a]
Tasa de reoxigenación kr 1.5 [1/dı́a]
Oxı́geno disuelto de saturación CS 9 [mg/L]

Tabla 1: Datos problema 4.

a) Determine una expresión para la materia orgánica L y oxı́geno disuelto OD en función de la

distancia en el ducto cerrado. Usando los datos de la Tabla 1, calcule la concentración de L y
OD al final del tramo. Hint recuerde que en el ducto no existe reaireacion).

Departamento de Ingenierı́a Civil 15 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Solución:

Para la materia orgánica -DBO- se obtiene una expresión del balance de masa que además
coincide con la expresión del modelo S-P:

L(t) = L0 · e−kd ·t

Para el caso del oxı́geno, como nos encontramos en un ducto cerrado, se tiene que kr = 0.
Reemplazando esto en la ecuación del modelo se obtiene la expresión:

D(t) = D0 − L0 · (e−kr ·t − 1)

Notando que podemos relacionar el tiempo con la distancia por medio de la velocidad, hacemos
el reemplazo de variables en la ecuación. Al reemplazar se obtiene:
z
L(t) = L0 · e−kd · v
z
D(t) = D0 − L0 · (e−kr · v − 1)
Al evaluar esto, se obtienen los siguientes valores:

D(x = Lducto ) = 6,00[mg/L]

L(x = Lducto ) = 49,997[mg/L]

b) Determine una expresión para la materia orgánica L y oxıgeno disuelto OD a la salida de la

laguna de estabilización. Considere reaireación. Evalúe para los datos de la Tabla 1.

Solución:

Se plantean balances de masa para el oxı́geno y la DBO en a laguna de la forma:

Caudales:
Qin = Qout

DBO:
Lin · Qin = Lout · Qout + V · kd · Lout

OD:
ODin · Qin + V · kr · (Cs − ODout ) = ODout · Qout + V · kd · Lout

Evaluando se obtiene los siguientes valores:

Llaguna = 49,987[mg/L]

ODlaguna = 5,9991[mg/L]

c) ¿Que expresión es valida para el canal? Evalúe la calidad del agua en la descarga al rı́o.

Departamento de Ingenierı́a Civil 16 Universidad de Chile

CI4102 - Ingenierı́a Ambiental Prof. Ana Lucı́a Prieto
Otoño 2020 Auxs. Beatriz Espinoza
Taller#2 Gino Vera

Solución:

Para este tramo, se utilizan las ecuaciones S-P, utilizando como condiciones iniciales las condi-
ciones de equilibrio de la laguna de estabilización.

Llaguna = L0

ODlaguna = OD0

d) ¿Se cumple el máximo permitido para la DBO?. Suponga también que se ha determinado que
en la descarga, los valores de oxı́geno disuelto deben ser superiores a los 3 mg/L. ¿Se cumple
esta condición?

Solución:

Con lo anterior, se obtienen los siguientes valores:

Lf inalcanal = 49,987[mg/L]

Df inalcanal = 5,991[mg/L]
ODf inalcanal = Cs − D = 9 − 5,991 = 3,009[mg/L]
Por lo tanto, no se cumple la condición impuesta para el OD, sin embargo, si cumple la condición
para la DBO.

e) Determine el volumen mı́nimo necesario del reactor para garantizar que la descarga al rı́o cumpla
con los niveles exigidos para la DBO.

Solución:

Para esto, se impone que el déficit y la DBO tengan las condiciones que se piden. Con esto, se
obtiene que los valores de DBO y de déficit al principio del canal sean:

L0 = 35,005[mg/L]

Para obtener este valor, se resuelve la ecuación de balance de masa de DBO anteriormente
utilizada, ahora sabiendo el L y teniendo como incognita el volumen. Al despejar se obtiene un
valor:
V olumen = 8, 634, 160[m3 ]

Departamento de Ingenierı́a Civil 17 Universidad de Chile