Oxígeno disuelto

Profesora: Ing. Adriana Reyes Serrano

1
 3.1 Indicadores de la calidad del agua

3.2 Demanda bioquímica de oxígeno

3.3 Reaireación atmosférica

Contenido 3.4 Fotosíntesis y respiración

3.5 Demanda por sedimentos

3.6 Modelo para el oxígeno disuelto en ríos, estuarios y lagos.

3.7 Enfoque del volumen de control

3.8 Ingeniería de control

2
3.1 Indicadores de la calidad del agua

De acuerdo a la
SEMARNAT, el ICA (Índice
de calidad del agua) se
define como el grado de
contaminación existente en
el agua a la fecha de un
muestreo, expresado como
un porcentaje de agua pura.

3
• Dicho índice consideró 18 parámetros para su cálculo con distintos pesos relativos (Wi):

Tabla 1. Parámetros de calidad del agua

Parámetro Peso (Wi) Parámetro Peso (Wi)
Demanda bioquímica de 5.0 Nitrógeno en 2.0
oxígeno nitratos
Oxígeno disuelto 5.0 Alcalinidad 1.0
Coliformes fecales 4.0 Color 1.0
Coliformes totales 3.0 Dureza total 1.0
Sustancias activas al azul de 3.0 Potencial de 1.0
metileno (Detergentes) Hidrógeno (pH)
Conductividad eléctrica 2.0 Sólidos suspendidos 1.0
Fosfatos totales 2.0 Cloruros 0.5
Grasas y aceites 2.0 Sólidos disueltos 0.5
Nitrógeno amoniacal 2.0 Turbiedad 0.5

4
 Parámetros de calidad del agua
• Los parámetros de calidad del agua residual que tienen importancia en los vertidos de aguas
residuales son:

DBO5 DQO Sólidos Bacterias

suspendidos

Metales Pesticidas Nutrientes

5
Efectos
Oxígeno disuelto
• Es importante para la vida acuática.
• A niveles por debajo de 4-5 mg/L, se
producen efectos nocivos en
determinadas especies.
Sólidos suspendidos
• Afectan la turbiedad de la columna de
agua y acaban por sedimentar en el
fondo.
• Origina toxicidad y demanda de
oxígeno de los sedimentos.

6
Bacterias coliformes

• Indicador de la presencia de otros

organismos patógenos de origen fecal.
Se utiliza para determinar la seguridad
de la utilización de un agua para uso
recreativo.

Nutrientes

• Ocasionan eutrofización y descensos del

nivel de OD.

Compuestos químicos tóxicos

• En diferentes concentraciones provocan

efectos nocivos para la vida acuática y
para el hombre en caso de que se ingiera
el agua y/o la vida animal acuática.

7
 Tabla 2. Contaminantes de importancia en aguas residuales
Contaminante Parámetro típico de Impacto ambiental
medida
Materia orgánica DBO,DQO Desoxigenación del agua, generación de olores
degradable indeseables.
Materia suspendida SST, SSV Causa turbiedad en el agua, deposita lodos.
Patógenos CF Hace el agua insegura para consumo y
recreación.
Amoniaco N𝐻4− - N Desoxigena el agua, es tóxico para organismos
acuáticos y puede estimular el crecimiento de
algas.
Fósforo Ortofosfatos Puede estimular el crecimiento de algas.
Materiales tóxicos Como cada material Peligroso para la vida vegetal y animal.
tóxico específico
Sales inorgánicas SDT Limita los usos agrícolas e industriales del
agua.
Energía térmica Temperatura Reduce la concentración de saturación de
oxígeno en el agua, acelera el crecimiento de
organismos acuáticos.
Iones hidrógeno pH Riesgo potencial para organismos acuáticos 8
• Tabla 3. Elementos y compuestos tóxicos
Normas y criterios

Las normas de calidad

Son establecidas para
del agua son conjuntos
mantener o mejorar la
de limitaciones, tanto
calidad de los cuerpos
cuantitativos como
receptores.
cualitativos.

10
El sistema de control de vertidos de aguas
residuales más común

consiste en establecer limitaciones de

calidad de los efluentes y conceder
permisos de vertido.

En Estados Unidos, el control de los

vertidos de aguas residuales corre a cargo
de la (NPDES), administrada por la EPA.

11
Tabla 4. Ejemplo de clasificación de los cuerpos de agua

• Los cuerpos de agua receptores se

clasifican en diversas categorías en
función de su uso, actual o previsto,
estableciéndose, para cada uno de
ellos, un conjunto de limitaciones para
la protección de dichos usos.

12
En el caso de los compuestos tóxicos, el enfoque puede ser a nivel
global del efluente, o particularizado para cada compuesto
químico específico.

En este segundo caso, se utilizan limitaciones individuales

para cada uno de los compuestos tóxicos específicos presentes
en el agua residual.

Para proteger la vida acuática de los efectos tóxicos agudos o

crónicos, y para proteger al ser humano de posibles problemas
nocivos para la salud, entre ellos el cáncer.

13
 Los niveles de toxicidad aguda para la vida acuática no se
deben superar, ni siquiera durante cortos periodos de tiempo.

Niveles Los niveles de toxicidad crónica se pueden tolerar durante

periodos de tiempo más largos, que actualmente se
de establecen en torno a los cuatro días.
toxicidad

Las limitaciones que hacen referencia la protección de la

salud humana se suelen basar en periodos de exposición de
muy larga duración, que pueden llegar a ser de toda la vida,
de modo que los límites de concentración correspondientes se
pueden aplicar a las concentraciones medias.
14
Tarea sesión 11.
• Leer las normas oficiales mexicanas:
• NOM-001-SEMARNAT-1996
• NOM-002-SEMARNAT-1996
• NOM-003-SEMARNAT-1997
• Elaborar en un documento Word lo siguiente:
• Colocar el título de la norma, el objetivo y campo de aplicación, así como
indicar cuáles son los límites máximos permisibles para cada contaminante.
Lo anterior se debe realizar para cada NOM.
• Referencia: en formato APA sexta edición.
• La fecha límite de entrega es: viernes 22 de mayo a las 15:00 h,
después de este tiempo no se recibe.

15
3.2 Demanda bioquímica de oxígeno

Es la cantidad de oxígeno que requieren los microorganismos para oxidar

(estabilizar) la materia orgánica biodegradable en condiciones aerobias.

Cuando se refiere a la DBO necesaria para oxidar todo el material orgánico

carbonáceo biodegradable, se denomina demanda bioquímica última de
oxígeno carbonácea (DBOUC).

En condiciones normales de laboratorio, esta demanda se cuantifica a 20°C, el

ensayo estándar se realiza a cinco días de incubación y se conoce
convencionalmente como DBO, con valores numéricos expresados en mg/L –O2.

16
 DBO
Parámetro más usado para medir la calidad de
aguas residuales y superficiales.

Determinar la cantidad de oxígeno requerido

para estabilizar biológicamente la materia
orgánica del agua.

Diseñar unidades de tratamiento biológico.

Para evaluar la eficiencia de los procesos de

tratamiento.

Fijar las cargas orgánicas permisibles en fuentes

receptoras.
17
• La formulación matemática de la DBO carbonácea fue hecha por Steeter y Phelps con base en la
ley empírica de Theriault de que la tasa de oxidación bioquímica de la materia orgánica es
directamente proporcional a la cantidad de materia orgánica biodegradable presente; es decir,
obedece una ecuación de primer orden. Matemáticamente se tiene:

𝑑𝐿𝑡
• = −𝐾 𝐿𝑡 (1)
𝑑𝑡
• Donde: Lt = DBO remanente en el agua para el tiempo t, mg/L

• K = constante que expresa la tasa de oxidación, d-1

• t = tiempo de oxidación, d
𝑑 𝐿𝑡
• = tasa de oxidación de la materia orgánica carbonácea, mg/L. d
𝑑𝑡

• Integrando la ecuación anterior se obtiene:

𝐿𝑡
• = 𝑒 −𝐾𝑡 = 10−𝑘𝑡 (2)
𝐿
18
• Donde L = DBO remanente en el agua para el tiempo t= 0 o DBOUC

K= 2,303 k

• La DBO remanente para cualquier tiempo t es igual a:

𝐿𝑡 = 𝐿 𝑒 −𝐾𝑡 = 𝐿 𝑋 10−𝑘𝑡 (3)

• La DBO ejercida en cualquier tiempo t es:

𝑦 = 𝐿 − 𝐿𝑡 = 𝐿 − 𝐿 𝑒 −𝐾𝑡 = 𝐿 − 𝐿 𝑥 10−𝑘𝑡 (4)

𝑦 = 𝐿 1 − 𝑒 −𝐾𝑡 = 𝐿 1 − 10−𝑘𝑡 (5)

• La ecuación (5) es la expresión usual de la DBO, donde:

• y= DBO ejercida por el agua después de un tiempo t, mg/L

• L= DBOUC, mg/L, o DBO remanente en el agua para t=0

• K= constante de velocidad de reacción de la DBO, base natural, d -1

• k= constante de velocidad de reacción de la DBO, base decimal, d -1

• t= tiempo de reacción, d

19
• La DBO estándar, DBO5, o simplemente
DBO, es igual a:

𝑦 = 𝐿 1 − 𝑒 −5𝐾 = 𝐿 1 − 10−5𝐾 (6)

• Evidentemente, la DBOUC es igual a la

DBO satisfecha más la DBO remanente:

𝐿 = 𝑦 + 𝐿𝑡 7

La fracción remanente para cualquier

𝐿
tiempo t será: 1 y la fracción oxidada será:
𝐿
𝐿𝑡
1−
𝐿

Figura 1. Formulación de la curva de DBO carbonácea.

• Gráficamente, la ecuación (7) se ilustra
en la figura 1.

20
• La oxidación bioquímica es un proceso lento que requiere, matemáticamente, un tiempo infinito
para su culminación. A 20°C, valores típicos de K y k son respectivamente 0.23 y 0.10 d-1, para
aguas residuales domésticas. Sin embargo, los valores de K pueden variar entre 0.05 y 0.3 d-1. En
la tabla 4 se incluyen otros valores de K citados por la literatura.

• Tabla 4. Valores típicos de k, K y L

• Para determinar el valor de la constante de reacción K a una temperatura diferente a 20°C se

utiliza la ecuación deducida de la relación clásica de Van´t Hoff Arrhenius:

21
Para una misma DBOUC, el
consumo de oxígeno es función
del tiempo de reacción y del
valor de K, como puede verse en
la figura 2.

Figura 2. Efecto del valor de K sobre la DBO (L= constante)

22
 DBO nitrogenácea
La descomposición de la materia orgánica, especialmente la
hidrólisis de las proteínas, produce material no carbonáceo
como el amoníaco.

Este material, nitrógeno amoniacal, es oxidado por las

bacterias nitrificantes en nitrito y nitrato, causando una
Características demanda de oxígeno conocida como demanda bioquímica de
oxígeno nitrogenácea (DBON).

En un cultivo mixto, como el usado normalmente para

determinar DBO, se tiene una población de bacterias
nitrificantes baja, porque la tasa de reproducción de las
bacterias nitrificantes es lenta, normalmente de 6 a 10 días.

23
 • En general, se ha afirmado
que las bacterias comunes
del agua residual tienen un
tiempo de generación de 7
horas y que las bacterias
nitrificantes, por otra parte,
tienen un crecimiento
mucho más lento y su
tiempo de generación es
mayor de 8 días. El efecto
de la DBON se observa
gráficamente en la figura 3.

Figura 3. Efecto de la DBON.

24
• Los recursos de oxígeno requeridos para satisfacer la DBON pueden determinarse utilizando los
análisis de las diferentes formas del nitrógeno y usando las relaciones siguientes:

• Mediante las reacciones anteriores, las bacterias nitrificantes obtienen la energía para su
sostenimiento y crecimiento celular: Globalmente:

• De acuerdo con las relaciones estequiométricas anteriores se requieren 4.57 mg/L-O2 para oxidar un
mg/L de nitrógeno amoniacal expresado como N y se consumen 7.14 mg/L de alcalinidad como CaCO3.

25
 Ejemplos
• 1. La DBO estándar de un agua residual (AR) es de 200 mg/L y la constante de
reacción base 10, la k= 0.1 d-1. Determinar la DBO a un día y la DBOUC.

• Solución:

Datos Fórmula

k=0.1 d-1 y= L (1-10-kt)

y= 200 mg/L donde:
y=DBO ejercida por el agua después de un tiempo t, mg/L
L= DBOUC, mg/L, o DBO remanente en el agua para t=0
k= constante de velocidad de reacción de la DBO, base
decimal, d-1
t= tiempo de reacción, d

26
• Para la DBOUC:
𝑦 200 𝑚𝑔/𝐿
• 𝐿= = = 293 𝑚𝑔/𝐿
(1− 10−𝑘𝑡 ) (1−10−0.1∗5 )

• La DBO a un día será:

• 𝑦1 = 𝐿 1 − 10−𝑘𝑡 = 293 1 − 10−0.1∗1 = 60 𝑚𝑔/𝐿

27
2. La DBO estándar de un agua residual es de 200 mg/L y el contenido de NTK de 20
mg/L-N. La constante de reacción de la DBO, base 10, es igual a 0.1 d -1.

Calcular la DBOUC, la DBON y la DBO total del agua residual.

Solución:

1. La DBOUC del residuo será:

𝑦 200 𝑚𝑔/𝐿
𝐿= = = 293 𝑚𝑔/𝐿
(1 − 10−𝑘𝑡 ) (1 − 10−0.1∗5 )

2. La DBON del residuo será:

𝐷𝐵𝑂𝑁 = 20 𝑋 4.57 = 91 𝑚𝑔/𝐿

3. La DBO total del agua residual será:

𝐷𝐵𝑂𝑇 = 293 + 91 = 384 𝑚𝑔/𝐿

28
Existen otras fuentes y sumideros de oxígeno que es necesario
tener en cuenta. Entre ellos, se incluyen la reaireación superficial,
la demanda de oxígeno de los sedimentos, la fotosíntesis y la
respiración, todos ellos analizados en lo que sigue.

29
Oxígeno disuelto
Factor más importante que debe
ser medido en el agua.

Sólo tiene valor si se mide con

la temperatura.

Oxígeno

Figura 4. Variación del

Para establecer el porcentaje porcentaje de saturación de
de saturación. oxígeno al disminuir la
temperatura con la profundidad
en la presa del embalse La Fe,
Antioquía (Colombia).
30
 Sirve como variable indicadora de la
salud global del río y de su habilidad
para mantener y propagar un
sistema ecológicamente equilibrado.
Oxígeno
disuelto Todos los modelos de OD pretenden
cuantificar el OD resultante de la
interacción de las funciones que
afectan su concentración en el río.

31
 Reaireación
atmosférica

Fuentes
de oxígeno
disuelto
Inyección
artificial de en el río
OD por Fotosíntesis
equipos de
aireación

32
 Consumo de
oxígeno

Respiración de Demanda Bentónica de

plantas y animales bioquímica oxígeno

33
El principal factor de
consumo de oxígeno
libre es la oxidación
de materia orgánica
por respiración.
34
 3.3 Reaireación atmosférica

• Con el fin de comprender el

mecanismo básico de la transferencia
de O2 entre los cuerpos de agua y la
atmósfera, se consideró un sistema
de dos fases cuyo esquema es
mostrado en la figura 5.

Figura 5. Interfaz líquido-gas

35
 Ley de Henry
• El peso de cualquier gas que se disuelve en un volumen dado de líquido, a
una temperatura constante, es directamente proporcional a la presión que el
gas ejerce sobre el líquido.

• El OD en el agua se comporta de acuerdo a la ley de Henry, es decir:

𝐶𝑠 = 𝐻𝑒 𝑝
Donde:

Cs= concentración de saturación de OD en mg/L

𝑚𝑔
He= constante de Henry en ( ∗ 𝑚𝑚𝐻𝑔)
𝐿

p= presión parcial de O2 en mmHg

36
 Es el proceso de
absorción de O2
disponible en el
aire, por un cuerpo
de agua.

Reaireación atmosférica

37
Este es uno de los
principales mecanismos de
ingreso de O2 en los
sistemas acuáticos, se
caracteriza por su ka.

La absorción del O2 de la
atmósfera por el agua se
asume normalmente como
un proceso de primer
orden, es decir, la tasa de
absorción es directamente
proporcional al OD.

38
• Se sabe que:

𝑂𝐷 = 𝐶𝑆 − 𝐶0
• Donde:
• OD= El déficit de oxígeno disuelto, en (mg/L)
• 𝐶𝑆 = La concentración de saturación de OD a una temperatura T, en (mg/L)
• 𝐶0 = La concentración inicial de OD, en (mg/L)

• Utilizando un modelo simplificado, se tiene la siguiente relación:

𝐷𝑜𝑑 𝑑𝐶𝑠
= − 𝑘𝑎 𝑂𝐷 −
𝐷𝑇 𝑑𝑡
• Donde:
• ka= coeficiente de reaireación de superficie
t= tiempo
39
•
 𝐷𝐶𝑠
• Suponiendo que no hay variación de la temperatura en el tiempo = 0, por lo tanto:
𝑑𝑡

𝐷𝑑𝑜
= − 𝑘𝑎 𝑂𝐷
𝐷𝑇
Para OD = OD0 en un tiempo t= 0, la solución es

𝑂𝐷 = 𝑂𝐷0 𝑒 −𝑘𝑎 𝑡

También puede escribirse en términos de concentración, es decir:

𝐶 = 𝐶𝑆 − 𝐶𝑆 − 𝐶0 𝑒 −𝑘𝑎 𝑡

40
 • En la figura 6 se puede
observar dos casos extremos
de la recuperación del OD
debido a los efectos que
genera el ka después de un
disturbio con el cual se
presentó un OD de acuerdo a
las ecuaciones antes
presentadas. En el caso de un
bajo ka se presenta una
recuperación muy lenta en el
tiempo, incluso podría no
recuperar su concentración de
OD. Por otro lado, un Ka alto
devuelve la concentración de
OD en poco tiempo, lo que
permite al río autodepurarse.

Figura 6. Efecto del ka sobre la concentración de OD.

41
 Velocidad de la corriente: cuanto
mayor sea la velocidad del agua en el
sistema natural, mayor será su
capacidad de reaireación.
Factores
que Profundidad del lecho: a mayor
afectan la profundidad menor será la
reaireación concentración de O2 transferido a la
superficial
columna de agua.

Velocidad del viento: cuanto mayor

sea la velocidad del viento, mayor será
también la reaireación. Factor
importante en lagos.
42
 •Producen una
Existencia
pequeñas presas
de
y
pequeña reaireación
cascadas extra.

Surfactantes (p.e.,
sustancia que reduce •Impiden la
la tensión superficial
de un líquido, y que
transferencia de O2 a
sirve como agente través de la interfase
humectante
detergente).
o aire-agua.

43
• En los casos en los que la concentración de oxígeno disuelto en un cuerpo de agua con una
superficie libre expuesta a la atmósfera se halla por debajo de la concentración de saturación, se
produce un flujo neto de oxígeno entre la atmósfera y el agua.
• Este flujo (masa por unidad de tiempo y de superficie de contacto agua-atmósfera), es proporcional
a la diferencia entre la concentración de OD presente y la concentraci6n de saturación.
• Por lo tanto, para un volumen de control con una superficie libre de contacto, la velocidad del
aumento de la cantidad de OD debido ala reaireación superficial se puede escribir como:

𝐴 𝑘𝑅
𝑟𝑅 = 𝑘𝑅 𝑐𝑆 − 𝑐 = 𝑐 − 𝑐 = 𝐾2 𝑐𝑆 − 𝑐
𝑉 𝐻 𝑆
donde
rR = ganancia de oxígeno por unidad de tiempo y de volumen de agua debida a la reaireación
superficial, M/TL3
kR = velocidad del flujo de reaireación, L/T
A= superficie de contacto con la atmósfera del volumen de control, L 2
V = volumen del volumen de control, L3

44
Cs = concentración de saturación del oxígeno disuelto, M/L 3

C = concentración de oxígeno disuelto, M/L3

H = profundidad del volumen de control, L.

K2= velocidad de reaireación superficial, 1/T.

45
3.4 Demanda por sedimentos
• También conocida como: demanda bentónica de oxígeno

Es el consumo de oxígeno
necesario para estabilizar la
materia orgánica contenida en
los sedimentos e involucra
procesos aeróbicos y anaeróbicos.

46
 Parte de los sólidos descargados con el agua residual
son sólidos orgánicos. Una vez que sedimentan en el
fondo se descomponen, dependiendo de las
condiciones, mediante procesos aerobios y anaerobios.

El oxígeno consumido en el proceso de descomposición

aerobia representa otro sumidero de oxígeno disuelto
Características del cuerpo de agua.

En el caso de un volumen de control en contacto con

el fondo del cuerpo de agua, la tasa de reducción del
oxígeno disuelto debida a la demanda de oxígeno de
los sedimentos viene dada por:
47
 𝒌𝑺
𝒓𝑺 =
𝑯
Donde

rS = consumo de oxígeno debido a la OD por unidad de tiempo y de volumen de agua, M/TL 3

kS = tasa de consumo de oxígeno de los sedimentos, M/L 2T.

H = profundidad del volumen de control, L.

• Los factores que influyen, principalmente, en el valor de Ks:

❖ El contenido orgánico de los sedimentos

❖ La temperatura

❖ La concentración de oxígeno disuelto en la interfase sedimentos-agua

❖ El crecimiento de comunidades biológicas, y

❖ La velocidad de la corriente.
48
 Criterios para determinar ks
• La medición del valor de kS se puede llevar a cabo utilizando una cámara de
flujo para aislar los sedimentos del agua situada por encima de ellos.

• En la cámara de flujo, se registra la variación con el tiempo de la

concentración de oxígeno disuelto, lo cual permite determinar el valor de kS.
Es preferible realizar mediciones in situ, pero su fiabilidad suele quedar en
entredicho debido a la variabilidad espacial y a que las tensiones
tangenciales que el flujo provoca sobre el fondo del cuerpo de agua son
difíciles de reproducir en la cámara de flujo.

• En estudios previos, se puede emplear un valor de ks del siguiente orden de

magnitud: 2-10 g/m2*d en zonas próximas a los puntos de vertido de las
aguas residuales municipales; 1-2 g/m2*d en zonas aguas abajo de los puntos
de vertido de aguas residuales municipales y para el fango natural de
estuarios; 0.2-1 g/m2*d para fondos arenosos, y 0.05-0.1 g/m2*d para fondos
minerales.

49
La tasa de consumo ks también se puede estimar a partir de las velocidades de deposición de
sólidos, suponiendo que la velocidad de descomposición es igual a la velocidad de deposición:

𝑘𝑆 = 𝑟0 𝑎𝑅𝑑

Donde:

r0 = contenido orgánico oxidable de los sólidos descargados, normalmente entre 0,5 y 0,6 para
efluentes secundarios y 0,8 para efluentes primarios.

a = relación estequiométrica entre oxígeno y sedimentos = 1,07.

Rd = velocidad de deposición de sólidos, M/UT.

50
3.5 Fotosíntesis y respiración
La incidencia de la
luz solar en los
sistemas acuáticos o
terrestres se
captura por los
pigmentos de las
plantas,
principalmente la
clorofila, y es
energía luminosa se
convierte en energía
química mediante la
Fotosíntesis fotosíntesis.

51
Proceso por el cual la energía
química que se almacenó
mediante la fotosíntesis es
finalmente liberada en las
plantas y en otros
organismos (desde bacterias
hasta plantas y animales).

Respiración

52
Los niveles de oxígeno
disuelto en el entorno se
pueden ver afectados

por el crecimiento de
algas (fitoplancton), y
hierbas (macrofitos)
que se alimentan a base
de amoníaco y nitratos.

53
Efectos en el oxígeno disuelto

Fotosíntesis Respiración

Las algas y las

hierbas Sumidero de
constituyen una oxígeno.
fuente de oxígeno.

Sobresaturación
Niveles bajos de
de oxígeno
oxígeno disuelto.
disuelto.

54
Figura 7. Ciclo del nitrógeno en medio acuático
55
Ecuación de conservación de la masa
✓ Laincorporación de los procesos de transporte y transformación, así como la
existencia de fuentes externas, permite obtener la forma general de la
ecuación de conservación de la masa:

𝜕𝐶 𝜕𝐶 𝜕𝐶 𝜕𝐶 𝜕 𝜕𝐶 𝜕 𝜕𝐶 𝜕 𝜕𝐶
= −𝑈 −𝑉 −𝑊 + 𝐸 + 𝐸 + 𝐸 + ෍ 𝑟𝑖 + ෍ 𝐼𝑗 (1)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝑥 𝜕𝑥 𝜕𝑦 𝑦 𝜕𝑦 𝜕𝑧 𝑧 𝜕𝑧

Donde:

I = velocidad de aportación de las fuentes externas: masa introducida por

unidad de tiempo y de volumen de agua, M/L3T.

i = índice del proceso de transformación.

j = índice de identificación de la fuente externa.

56
Como ejemplo, en el caso en el que el constituyente considerado es el oxígeno
disuelto, el termino relativo a las transformaciones es:

෍ 𝑟𝑖 = −𝑟0 − 𝑟𝑅 − 𝑟𝑆 − 𝑟𝑃 − 𝑟𝑃𝑝 (2)

Donde:

rP = velocidad de producción de oxígeno por fotosíntesis por unidad de tiempo

y de volumen de agua, M/L3T.

RPp= velocidad de consumo de oxígeno por respiración por unidad de tiempo y

de volumen de agua, M/L3T.

57
 3.6 Modelo para el oxígeno
disuelto en ríos, estuarios y lagos

58
Vertido a lagos y embalses

En muchos lugares en
los que no existen cursos
de agua en las
proximidades, puede ser
necesario verter las
aguas residuales
tratadas en lagos o
embalses.

59
Análisis de mezcla completa

El agua de lagos y embalses pequeños tiende estar bien mezclada

como consecuencia de la turbulencia inducida por el viento.

Los lagos de mayor profundidad tienden a estratificarse durante

el verano; sin embargo, en muchos de ellos se producen dos
inversiones anuales de la estratificación.

En este tipo de análisis, se supone que las concentraciones de

constituyentes son aproximadamente homogéneas en todo el lago o
embalse.

60
 𝑑𝐶 1 𝑑𝑉
= −𝐾𝐶 + ⅀𝐼𝐽 + 𝐶
𝑑𝑡 𝑉 𝑑𝑡
El último término de esta ecuación tiene en cuenta las variaciones del volumen del
embalse producidas por la irregularidad de las aportaciones, los caudales aliviados, las
aportaciones pluviométricas y la evaporación.
Las aportaciones externas que hay que tener en cuenta son las precipitaciones, los cursos
de agua tributarios, la escorrentía superficial, la descarga de aguas subterráneas y los
caudales de aguas residuales (con caudales Qp, Qs, Qr, Qg, y Qw, y concentraciones Cp,
Cs, Cr, Cg, y Cw para cualquier constituyente de interés respectivamente.

Figura 8. Modelo de mezcla

completa para lagos y embalses. 61
• Los caudales extraídos del lago se expresan como una aportación con signo negativo –Q0 y una
concentración C igual a la concentración en el lago.

• La evaporación de la superficie del lago no suele representar una pérdida importante. Por lo
tanto, el término correspondiente a la aportación externa neta es:

𝑄𝑝 𝐶𝑝 + 𝑄𝑟 𝐶𝑟 + 𝑄𝑔 𝐶𝑔 + 𝑄𝑤 𝐶𝑤 − 𝑄0 𝐶 𝑀´ − 𝑄0 𝐶
⅀𝐼𝑗 = = (3)
𝑉 𝑉

Donde:
V= volumen del lago, L3.
M´= carga de constituyente al lago, M/T
= Qp Cp + Qs Cs + Qr Cr + Qg Cg + Qw Cw

62
Sustituyendo en la ecuación de balance de masa general (1) e
integrando, se obtiene:

𝑀′ 𝑡
− ‫׬‬0 𝛽𝑑𝜏
𝑡
− ‫׬‬0 𝛽𝑑𝜏 (4)
𝐶= 1−𝑒 + 𝐶0 𝑒
𝑉𝛽

Donde:
1 𝑑𝑉
𝛽=𝐾+ 𝑄0 +
𝑉 𝑑𝑡

C0 = concentración inicial (en el instante t = 0), M/L3

En los casos en los que el volumen del lago permanece constante,

dV/dt = 0, de forma que (J es constante y las integrales de los
𝑡
exponentes pasan a ser ‫׬‬0 𝛽𝑑𝜏 = 𝛽𝑡 en la figura 8, se representa la Figura 8. Modelo de
ecuación (4) para dos valores diferentes de la concentración inicial. mezcla completa para
La concentración converge asintóticamente al valor de la lagos y embalses.
concentración de equilibrio Ce = M'/βV, que se obtiene haciendo
tender t a infinito en la ecuación (4).

63
 Casi la totalidad de los lagos y embalses con
profundidades de 5 m o superiores presentan
estratificaciones durante una parte importante
del año.

Estratificación La estratificación se desarrolla durante la

de lagos
primavera, como consecuencia del aumento de
temperatura que se produce por efecto de la
radiación solar y atmosférica.

Debido a que la densidad del agua desciende

con la temperatura, se desarrolla una situación
hidrodinámicamente estable en la que el líquido
más ligero se sitúa por encima del más pesado.
64
• Por ejemplo, si se hace descender un
estrato de fluido hacia una zona
ocupada por un fluido más denso,
aparecerá una fuerza boyante
ascensional que se opondrá al
mezclado vertical restaurando el
equilibrio.
• Por lo tanto, se desarrolla una
estructura térmica vertical con un
estrato superior bien mezclado
(epilimnion), que se sitúa por
encima de una zona de rápido Figura 9. Ciclo de estratificación en un lago
descenso de la temperatura profundo.
(termoclina), situada a su vez sobre
un estrato de agua más fría y densa
(hipolimnion).
• Las corrientes de densidad y la
extracción selectiva son otros
fenómenos asociados a la
estratificación de los lagos.
65
 Ríos y estuarios

Los ríos y estuarios suelen tener una dimensión

preponderante, constituyendo cuerpos de agua mucho
más largos que anchos.

Modelo Las aportaciones de efluentes de plantas de tratamiento

unidimensional de aguas residuales u otras fuentes se mezclan
rápidamente en la sección transversal, lo cual justifica
la adopción de modelos unidimensionales.

Sólo se pretende calcular las variaciones longitudinales

de las concentraciones de los constituyentes, tomando
como punto de partida los valores medios en cada
sección.

66
 𝜕𝐶 𝜕𝐶 𝜕 𝜕𝑐
= −𝑈 + 𝐸𝑥 + 𝐸𝐿 − 𝐾𝐶 + 𝛴𝐼 (5)
𝜕𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑥

Donde:
x= distancia longitudinal a lo largo del río o estuario, L.
EL= coeficiente de dispersión longitudinal, L2/T.

En los cursos de agua naturales, la dispersión suele estar provocada por las variaciones laterales de la
velocidad. Para estimar el valor del coeficiente de dispersión se puede emplear la siguiente ecuación:
𝑈2 𝐵2
𝐸𝐿 = 0.011
𝐻𝑢 ∗
Donde:

EL= coeficiente de dispersión longitudinal, L2/T u* = velocidad tangencial, L/T= 𝑔 𝐻𝑠

U= velocidad media transversal, L/T g= aceleración de la gravedad, L/T 2

B= anchura del curso de agua, L. s= pendiente del curso de agua, L/L
H= profundidad del curso de agua, L
67
✓ En estuarios, las inversiones de flujo provocadas por la marea, así como las corrientes
secundarias provocadas por los gradientes de salinidad, tienden a aumentar la dispersión.
Normalmente, la dispersión es mucho mayor que la difusión turbulenta, de modo que se puede
despreciar el valor de Ex frente al de EL.

68
 -Desoxigenación carbonácea
Modelos -Reaireación atmosférica
más
usados -Fotosíntesis y respiración vegetal
para
simular
OD en
ríos -Demanda béntica de oxígeno
-Demanda de oxígeno por nitrificación

69
Aplicaciones de modelos de
simulación de OD
Evaluación de alternativas de control de contaminación para
visualizar los grados de tratamiento requerido.

Evaluación de alternativas de uso del agua para formular

normas de calidad del agua.

Cuantificación de caudales de dilución para control de

polución, y

Análisis de la localización óptima de plantas de tratamiento de

aguas residuales

70
Modelo de Streeter y Phelps

Publicado en 1944

Aplicable a ríos con condiciones aerobias

Flujo con superficie libre

Unidimensional

71
✓ Latasa de cambio del déficit de OD es igual a la diferencia entre la desoxigenación del agua, debida a
la remoción de DBO y materia orgánica carbonácea por oxidación biológica, y la reaireación o
suministro de OD proveniente de la atmósfera y causada por el déficit de OD y por la turbulencia.

𝑑𝐷
= 𝐾1 𝐿 − 𝐾2 𝐷
𝐷𝑇

𝐾1 𝐿0
𝐷= 𝑒 −𝐾1 𝑡 − 𝑒 −𝐾2𝑡 + 𝐷0 𝑒 −𝐾2𝑡
𝐾2 − 𝐾1

Donde: D: déficit de OD para el tiempo t, mg/L

L0= DBOUC inicial, t=0, mg/L
L= DBOUC remanente, t= t, mg/L
t= tiempo de flujo, d
K1= constante de desoxigenación, base e, d-1
K2= constante de reaireación, base e, d-1
D0= déficit inicial de OD, mg/L
72
• El punto crítico (punto de déficit máximo de OD, punto de concentración mínima de OD) ocurrirá
cuando el tiempo de flujo sea t= tc. Para este punto, las tasas de desoxigenación y reaireación son
iguales, el déficit de OD no cambia y:

𝑑𝐷
= 0 = 𝐾1 𝐿 − 𝐾2 𝐷𝑐
𝑑𝑡

𝐾1
𝐷𝑐 = 𝐿
𝐾2

𝐾1
𝐷𝑐 = 𝐿0 𝑒 −𝐾1 𝑡𝑐
𝐾2

Donde:

Dc= déficit crítico de OD, mg/L

tc= tiempo para el cual se alcanza el déficit crítico, d

73
• El valor de tc puede determinarse obteniendo dD/dt e igualándolo a cero:

1 𝐾2 𝐷0 𝐾2 − 𝐾1
𝑡𝑐 = 𝐼𝑛 [ 1− ]
𝐾2 − 𝐾1 𝐾1 𝐾1 𝐿0

La distancia al punto crítico (XC):

𝑋𝑐 = 𝑣𝑡𝑐

El valor de las constantes se puede obtener, para diferentes temperaturas por las ecuaciones:

74
Tabla 1. Valores de las constantes para el modelo de Streeter y Phelps

Tabla 2. Valores de K2 y f para el modelo de oxígeno disuelto

75
 Carga máxima permisible de DBO
Thomas, suponiendo que no existe deposición de lodos, dedujo una expresión para calcular
aproximadamente el valor de la carga máxima permisible de DBO, que puede introducirse en un
río, para mantener una concentración de oxígeno disuelto por encima de un estándar específico.

El cálculo por la fórmula de Thomas es aceptable en ríos con una sola carga puntual de aguas
residuales y características físicas muy constantes.
76
Ejemplo
• Un río con un caudal de 25 m3/s, a una temperatura de 14°C, tiene un OD de
7.0 mg/L y una DBOUC de 4.0 mg/L aguas arriba de la descarga de un agua
residual. El caudal de aguas residuales es de 5 m3/s, con DBO de 150 mg/L,
temperatura de 18°C y OD de 2.0 mg/L. La velocidad del río es de 0.5 m/s, la
profundidad media del flujo 1.0 m y la ecuación empírica para calcular K2 es:

𝑉
𝐾2 = 2.2
𝐻1.33

• Donde: K2= constante de reaireación, base neperiana, d-1

• V= velocidad de flujo, m/s

• H= profundidad del flujo, m

77
Los ensayos de DBO en el laboratorio indican un valor de la constante de desoxigenación, K1= 0.25
d-1 (base neperiana). El valor de la saturación de OD es de 7.5 mg/L

• Determinar:

• Déficit inicial de OD

• Tiempo de flujo al punto crítico

• Distancia al punto crítico

• Solución:

• Elaborar el diagrama de flujo Descarga aguas

residuales (2)

Río (1)

78
• Se calcula la temperatura de la mezcla:

𝑇 1 ∗ 𝑄1 + 𝑇2 ∗ 𝑄2 14 ∗ 25 + 18 ∗ 5
𝑇= = = 15°𝐶
𝑄1 + 𝑄2 25 + 5

• Se calculan las constantes:

• 𝐾1, 𝑇 = 𝐾1.20(1.047)𝑇−20 y 𝐾2 = 𝐾2.20(1.022)𝑇−20

• 𝐾1 = 0.25 (1.047)15−20 = 0.20 𝑑 −1

𝑉 2.2∗0.5
• 𝐾2 = 2.2 = = 1.1 𝑑 −1
𝐻 1.33 (1)1.33

• 𝐾2 = 1.10 (1.022)15−20 = 0.99 𝑑 −1

79
• Calcular el déficit inicial de OD

𝐷𝑜 = 𝐶𝑠 − 𝐶 = 7.5 − 𝐶

𝑂𝐷1 𝑄1 + 𝑂𝐷2 𝑄2
𝐶=
𝑄1 + 𝑄2

7 ∗ 25 + 2 ∗ 5
𝐶= = 6.17 𝑚𝑔/𝐿
25 + 5

𝐷𝑜 = 7.5 − 6.17 = 𝟏. 𝟑𝟑 𝒎𝒈/𝑳

▪ Calcular la DBOUC del agua residual

𝑌
𝐿𝐴𝑅 =
1 − 𝑒 −𝐾𝑡

80
 150
𝐿𝐴𝑅 = = 210.2 𝑚𝑔/𝐿
1 − 𝑒 −0.25∗5

▪ Calcular DBOUC de la mezcla

𝑌1 ∗ 𝑄1 + 𝑌1 𝑄1
𝐿0 =
𝑄1 + 𝑄2

4 ∗ 25 + 210.2 ∗ 5
𝐿0 = = 38.3 𝑚𝑔/𝐿
25 + 5

81
• Calcular el tiempo de flujo al punto crítico

1 𝐾2 𝐷0 𝐾2 − 𝐾1
𝑡𝑐 = 𝐼𝑛 [ 1− ]
𝐾2 − 𝐾1 𝐾1 𝐾1 𝐿0

1 0.99 1.33 0.99 − 0.20

𝑡𝑐 = 𝐼𝑛 1− = 𝟏. 𝟖𝟒 𝒅
0.99 − 0.20 0.20 0.20 ∗ 38.3

• Calcular la distancia al punto crítico

𝑋𝑐 = 𝑣𝑡𝑐

𝑚 𝑠
𝑋𝑐 = 0.5 1.84 𝑑 86400 = 79488𝑚 = 𝟕𝟗. 𝟓 𝒌𝒎
𝑠 𝑑

82
Tarea sesión 14
• Un río con un caudal de 30 m3/s, a una temperatura de 12°C, tiene un OD de
8.0 mg/L y una DBOUC de 6.0 mg/L aguas arriba de la descarga de un agua
residual. El caudal de aguas residuales es de 7 m3/s, con DBO de 180 mg/L,
temperatura de 20°C y OD de 3.0 mg/L. La velocidad del río es de 1 m/s, la
profundidad media del flujo 1.0 m y la ecuación empírica para calcular K2 es:

𝑉
𝐾2 = 2.2
𝐻1.33

• Donde: K2= constante de reaireación, base neperiana, d-1

• V= velocidad de flujo, m/s

• H= profundidad del flujo, m

83
Los ensayos de DBO en el laboratorio indican un valor de la constante de desoxigenación, K1= 0.25
d-1 (base neperiana). El valor de la saturación de OD es de 7.5 mg/L

• Determinar:

• Dibujar el diagrama de flujo

• Déficit inicial de OD

• Tiempo de flujo al punto crítico

• Distancia al punto crítico

• NOTA: La fecha límite de entrega es: viernes 29 de mayo a las 16:00 h, después
de este tiempo no se recibe.

84
 Agua residual Q= 7 m /s
3

T= 20 °C
OD= 3.0 mg/L
DBO= 180.0 mg/L

Río

Q= 30 m3/s V= 1 m/s
T= 12°C H= 1.0 m
OD= 8.0 mg/L ODS= 7.5 mg/L
DBOUC= 6.0 mg/L

Solución
Elaborar el diagrama
85
• Cálculo de la temperatura de la mezcla:

𝑇 1 ∗ 𝑄1 + 𝑇 2 ∗ 𝑄2 12 ∗ 30 + 20 ∗ 7
𝑇= = = 14°𝐶
𝑄1 + 𝑄 2 30 + 7

• Cálculo de las constantes:

• 𝐾1 = 0.25 (1.047)14−20 = 0.19 𝑑 −1

𝑉 2.2∗1
• 𝐾2 = 2.2 = = 2.2 𝑑 −1
𝐻1.33 (1)1.33

• 𝐾2 = 2.2 (1.022)14−20 = 1.93 𝑑 −1

86
• Cálculo del déficit inicial de OD

𝑂𝐷1 𝑄1 + 𝑂𝐷2 𝑄2
𝐶=
𝑄1 + 𝑄2

8 ∗ 30 + 3 ∗ 7
𝐶= = 7.05 𝑚𝑔/𝐿
30 + 7

𝐷𝑜 = 7.5 − 7.05 = 𝟎. 𝟒𝟓 𝒎𝒈/𝑳

▪ Cálculo del DBOUC del agua residual

180
𝐿𝐴𝑅 = = 252.3 𝑚𝑔/𝐿
1 − 𝑒 −0.25∗5

87
▪ Cálculo de DBOUC de la mezcla

6 ∗ 30 + 252.3 ∗ 7 𝑚𝑔
𝐿0 = = 53
30 + 7 𝐿

• Cálculo del tiempo de flujo al punto crítico

1 1.93 0.45 1.93 − 0.19

𝑡𝑐 = 𝐼𝑛 1− = 𝟏. 𝟐𝟖 𝒅
1.93 − 0.19 0.19 0.19 ∗ 53

• Cálculo de la distancia al punto crítico

𝑚 𝑠
𝑋𝑐 = 1.0 1.28 𝑑 86400 = 110,599 𝑚 = 𝟏𝟏𝟎. 𝟓 𝒌𝒎
𝑠 𝑑

88
 Referencias
• http://dgeiawf.semarnat.gob.mx:8080/ibi_apps/WFServlet?IBIF_ex=D3_R_AGUA0
5_01&IBIC_user=dgeia_mce&IBIC_pass=dgeia_mce

• Eddy, & Metcalf. (1995). Ingeniería de aguas residuales. España: McGRAW-HILL

• Romero-Rojas, J.A. (2000). Tratamiento de aguas residuales: Teoría y principios de

diseño. Colombia: Editorial Escuela Colombiana de Ingeniería.

• Anguiano Valencia, J. R. (2012). Determinación del coeficiente de reaireación en un

tramo corto de un río caudaloso con flujo turbulento, mediante el uso de los
trazadores hexafluoruro de azufre (SF6) y Rodamina WT. Morelia, Michoacán:
Universidad Michoacana de San Nicolás de Hidalgo.

• Roldán Pérez, G., & Ramírez Restrepo, J. J. (2008). Fundamentos de limnología

neotropical. Medellín, Colombia: Universidad de Antioquía .

89