Universidad de Chile

Departamento de Ingeniería Civil

Otoño 20012

TEMA 2.10
MODELOS DE CALIDAD DEL AGUA

CI41B Ingeniería Ambiental

Profesores C. Espinoza, D. Rodríguez y R. Muñoz
Semestre Otoño 2005

1. INTRODUCCION

Un importante número de problemas relacionados con el control de la contaminación de agua

superficial involucra la descarga de residuos líquidos que para todos los efectos prácticos no
decaen en el tiempo o no desaparecen del curso de agua por sedimentación., adsorción u otro
proceso. Tales residuos se conocen como conservativos. Ejemplos de compuestos
conservativos son los sólidos disueltos y los cloruros. Algunos contaminantes como pesticidas y
herbicidas, los cuales decaen muy lentamente en el tiempo, también se consideran compuestos
conservativos.

Otros problemas de control de contaminación se refieren a aquellas sustancias que decaen con
el tiempo, entre lo que se puede mencionar a la demanda bioquímica de oxígeno y las bacterias
coliformes. Estos compuestos contaminantes se conocen como no conservativos.

Para formular una relación determinística entre los residuos descargados al ambiente y la
calidad del agua resultante, ambas en términos de un mismo compuesto, podemos considerar
un pequeño volumen de río donde la coordenada x se alinea con la longitud del río, la
coordenada y con su ancho, y la coordenada z con la profundidad del río, tal como se muestra
en la Figura 1.

Figura 1
Volumen Infinitesimal en un Río

En muchos ríos y esteros, es posible considerar al cuerpo de agua como homogéneo con
respecto a las variables de calidad del agua en las direcciones transversal y vertical. El

CI41B INGENIERIA AMBIENTAL 1

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

compuesto contaminante varía o "tiene un gradiente" sólo en la dirección del río. Este tipo de
sistema se conoce como unidimensional, en contraste con cuerpos de agua con variaciones
importantes en las tres direcciones o sistemas tridimensionales. La Figura 2 muestra un sistema
unidimensional en el cual el área de escurrimiento es A y el espesor del elemento de control es
∆x .
Figura 2
Elemento Infinitesimal en un Sistema Unidimensional

Para muchas situaciones relacionadas con la calidad del agua en ríos y esteros es posible
asumir que el medio es no dispersivo, es decir, cada elemento del material contaminante fluye
en la dirección de la corriente (agua abajo). No existe mezcla debido a difusión o dispersión.
Esto se conoce como flujo pistón y se ilustra en la Figura 3.

Figura 3
Ilustración de Flujo Pistón en un Río

CI41B INGENIERIA AMBIENTAL 2

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

2. MODELOS DE CALIDAD DE AGUAS - COMPUESTOS INDIVIDUALES

2.1 Balance de Masas

La masa total del compuesto C [M/L3] que entra a través de la sección de control en la Figura 2,
en un intervalo de tiempo ∆t es:

Q ⋅ C ⋅ ∆t (1)

donde Q es el caudal en el río. De manera similar, la masa de contaminante que sale del
volumen de control a través de una cara situada a una distancia ∆x de la zona de entrada
estaádada por:

(Q + ∆Q) ⋅  C + ∂C 
∆x  ⋅ ∆t (2)
 ∂x 

A partir de la definiciones anteriores, la variación de la masa del compuesto C dentro del

volumen de control debido a la entrada y salida de contaminante, así como también debido a
reacciones simples, sobre el intervalo de tiempo ∆t es:

 ∂C 
V ⋅ ∆C = Q ⋅ C ⋅ ∆t − (Q + ∆Q ) ⋅  C + ∆x  ⋅ ∆t ± K ⋅ V ⋅ C ⋅ ∆t (3)
 ∂x 

donde K [1/T] representa el coeficiente de una reacción de primer orden que describe la pérdida
o ganancia del compuesto C.

Al expandir la ecuación (3) y dividir por el intervalo de tiempo, ∆t , y el volumen V ( = A ⋅ ∆x ), se

tiene:

∆C Q ∂C C ∆ Q ∆ Q ∂C
=− ⋅ − ⋅ − ⋅ ± K ⋅C (4)
∆t A ∂x A ∆ x A ∂x

Si ahora tomamos el límite cuando el tamaño de la separación, ∆x , y el intervalo de tiempo,

∆t , se aproximan a cero podemos escribir:

∂C Q ∂C C ∂Q
=− ⋅ − ⋅ ± K ⋅C (5a)
∂t A ∂x A ∂x

o escrito en forma más reducida:

∂C 1 ∂(Q ⋅ C )
=− ⋅ ± K ⋅C (5b)
∂t A ∂x

Las ecuaciones (5a) y (5b) representan la ecuación diferencial básica para el estudio de un
sistema simple y no dispersivo. Como se indicó anteriormente, este modelo asume que no
existe mezcla (o dispersión) y se conoce comúnmente como modelo advectivo o modelo de
máximo gradiente.

CI41B INGENIERIA AMBIENTAL 3

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

Si la variable C representa un compuesto conservativo, K=0. El flujo, el área perpendicular, así

como el coeficiente de decaimiento pueden ser funciones de la distancia y tiempo.

Comúnmente se considera el origen del eje x en la ubicación de una descarga. De esta manera,
una condición de borde característica es:

C = C 0 (t ) en x = 0 (6)

donde C0 (t ) es la concentración del compuesto C en la localización de la descarga.

2.2 Régimen Permanente, Compuesto Conservativo

En esta situación, K es igual a cero, y los otros dos parámetros de la ecuación (5), así como la
concentración de la descarga se suponen constantes. Para efectos de simplificar aún más el
análisis supongamos que el resto de los parámetros no varían con el espacio.

De esta manera, el supuesto de régimen permanente permite escribir:

∂C
=0 y C = C0 en x=0
∂t

El supuesto de coeficientes constantes en el espacio permite escribir:

∂Q
=0 y A = A0
∂x

De esta manera, la ecuación (5) se puede escribir como:

Q dC dC
0=− ⋅ = −u ⋅ (7)
A0 dx dx

donde u es la velocidad de la corriente.

La solución de la ecuación (7) es obviamente una constante, el valor de la cual se puede

determinar a partir de la condición de borde en la descarga. La concentración C0 se puede
calcular con un balance de masas en la descarga tal que:

Q R ⋅ C R + QD ⋅ CD
C0 = (8)
QR + Q D

donde QR y C R son el caudal y la concentración en el río antes de la descarga, mientras que

QD y C D representan el caudal y la concentración de la descarga misma.

La Figura 4 ilustra en forma gráfica la solución anterior. De esta manera, en la Figura 4a se

tiene un río en el cual existen dos descargas y un tributario que se ubica aguas abajo de las
descargas. La Figura 4b muestra el aumento de flujo debido a las descargas y el tributario,
mientras que la Figura 4c muestra la disminución de la concentración en el río, en línea sólida,

CI41B INGENIERIA AMBIENTAL 4

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

por efecto de dilución debido a las descargas y el tributario. Las Figuras 4d y 4e muestran los
efectos individuales de las descargas 1 y 2, mientras que la Figura 4f muestra el efecto
combinado de las distintas descargas y del tributario sobre la calidad del agua en el río. En
línea punteada se observa la variación de la concentración en el caso de un compuesto no
conservativo.

2.3 Régimen Permanente, Compuesto No Conservativo

En este caso la ecuación (5) se puede escribir como:

Q dC C dQ
0=− ⋅ − ⋅ − K ⋅C (9)
A dx A dx

Si nuevamente se supone que el caudal y el área transversal son constantes, y sólo se

considera descargas puntuales, la ecuación (9) se puede escribir como:

Q dC dC
0=− ⋅ − K ⋅ C = −u ⋅ − K ⋅C (10)
A0 dx dx

con la condición de borde:

C = C0 en x = 0

La integración de la ecuación (10) permite escribir:

 K 
C (x ) = C 0 ⋅ exp − ⋅ x (11)
 u 

Al sustituir la expresión para la concentración C0 que se presenta en la ecuación (8) se puede

obtener:

 Q ⋅ C + QD ⋅ C D   K 
C (x ) =  R R  ⋅ exp  − ⋅ x  (12)
 QR + QD   u 

En contraste con la ecuación (8) se observa que para el caso de un compuesto no conservativo,
su concentración decae en forma exponencial con la distancia en el río. La ecuación (12) es
válida en todo un tramo de río hasta que una nueva descarga cambia las condiciones de borde.
Los esquemas de la Figura 4 muestra, con línea punteada, la variación de la concentración de
un compuesto no conservativo.

CI41B INGENIERIA AMBIENTAL 5

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

Figura 4
Variación de la Calidad del Agua en un Río

CI41B INGENIERIA AMBIENTAL 6

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

2.4 Régimen Impermanente o Transiente

Para analizar la situación de régimen transiente es conveniente volver a la ecuación general

obtenida a partir de un balance de masas:

∂C Q ∂C C ∂Q
=− ⋅ − ⋅ ± K ⋅C (13)
∂t A ∂x A ∂x

Para efectos de su solución analítica consideremos una situación más simple en la que los
parámetros Q, A y K son constantes en el tiempo, pero pueden variar en el espacio. La
descarga de residuos líquidos varía en el tiempo según una relación del tipo W (t ) . De esta
manera, la ecuación (13) se puede escribir como:

∂C Q ∂C
=− ⋅ − (K + ν ) ⋅ C (14)
∂t A ∂x

donde
1 ∂Q
ν (x) = ⋅
A ∂x

C = C 0 (t ) en x = 0

W (t ) W (t )
C0 (t ) = =
QR + QD QM

suponiendo que CR , la concentración en el río, es nula.

La solución de la ecuación (17) se puede obtener mediante el método de las características con
lo que se tiene:

( )
C (x , t ) = C0 t − t * ⋅ exp [− Φ 1 ( x )] (15a)

C (x , t ) =
(
W t − t* )
⋅ exp [− Φ 1 ( x )] (15b)
QM

donde

x
K +ν
Φ1 (x) = ∫ dξ (15c)
0
u

CI41B INGENIERIA AMBIENTAL 7

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

x
dx
t * = tiempo de viaje a x = ∫ (15d)
0
u

3. MODELOS DE CALIDAD DE AGUAS - SISTEMAS ACOPLADOS

3.1 Balance de Masas

Sistemas acoplados se refiere a aquellos fenómenos de calidad de aguas en los que un

parámetro específico está forzado por el valor de otro compuesto que reacciona
simultáneamente. Ejemplos de estos sistemas son la demanda bioquímica de oxígeno (DBO,
medida de la materia orgánica en un cuerpo de agua) y el oxígeno disponible para
estabilización de la materia orgánica. Asimismo se tiene a interacción entre polifosfatos y
ortofosfatos en un río.

El procedimiento seguido para construir un modelo matemático de un sistema acoplado es

idéntico a aquel correspondiente a un sistema individual. Se desarrolla un balance de masas
sobre un volumen elemental de río para cada uno de los compuestos considerados en este
análisis. Asimismo, no existe dilución y las reacciones que existen son de primer orden. Todo lo
anterior permite escribir las siguientes ecuaciones, donde C1 y C2 son las concentraciones de
los dos compuestos estudiados:

∂C1 1 ∂(Q ⋅ C1 )
=− ⋅ − K11 ⋅ C1 (16a)
∂t A ∂x

∂C2 1 ∂ (Q ⋅ C2 )
=− ⋅ − K 22 ⋅ C2 + K12 ⋅ C1 (16b)
∂t A ∂x

con:

C1 (0, t ) = C01 (t ) en x = 0

C2 (0, t ) = C 02 (t ) en x = 0

La solución de la ecuación (16a) aparece como una función forzante para definir el
comportamiento de la variable o parámetro C2.

3.2 Régimen Permanente

Bajo el supuesto de una situación en régimen permanente, así como parámetros constantes,
las ecuaciones (16a) y (16b) se pueden escribir como:

dC1
−u⋅ − K11 ⋅ C1 = 0 (17a)
dx

CI41B INGENIERIA AMBIENTAL 8

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

dC2
−u⋅ − K 22 ⋅ C 2 + K12 ⋅ C1 = 0 (17b)
dx

La solución de la ecuación (17a) con una condición de borde dada por:

C1 = C01 en x = 0

nos permite obtener:

 K 
C1 (x ) = C01 ⋅ exp  − 11 ⋅ x (18)
 u 

Si las concentraciones aguas arriba de la descarga son despreciables, entonces:

W1
C01 =
QM
y,

 K 
C1 (x ) =
W1
⋅ exp  − 11 ⋅ x  (19)
QM  u 

La ecuación (19) actúa como una función forzante para la ecuación (17b), lo que permite
escribir:

dC2 W  K 
u⋅ + K 22 ⋅ C 2 = K 12 ⋅ 1 ⋅ exp − 11 ⋅ x  (20)
dx QM  u 

sujeta a la condición de borde:

C2 = C02 en x = 0

Esta ecuación puede ser integrada para obtener:

  K   K  W  K 
C2 ( x ) =
K 12
⋅ exp  − 11 ⋅ x  − exp  − 22 ⋅ x  ⋅ 1 + C 02 ⋅ exp  − 22 ⋅ x  (21)
K 22 − K 11   u   u   QM  u 

Esta última ecuación, cuando se aplica al problema de consumo de oxígeno durante la

estabilización de materia orgánica en un río, se denomina Ecuación de Streeter y Phelps.

3.3 Régimen Impermanente o Transiente

Al igual que para el caso de un compuesto individual supongamos que Q, A, Kij son constantes
en el tiempo pero variables en el espacio. Considerando lo anterior podemos escribir:

CI41B INGENIERIA AMBIENTAL 9

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

∂C1 Q ∂C
= − ⋅ 1 − (K 11 + ν ) ⋅ C1 (22a)
∂t A ∂x

∂C 2 Q ∂C
= − ⋅ 2 − ( K 22 +ν ) ⋅ C 2 + K12 ⋅ C1 (22b)
∂t A ∂x
y

C1 = C 01 en x = 0 y t =0

C2 = C02 en x = 0 y t =0

La solución para C2 puede ser obtenida con la aplicación de Transformadas de Laplace y está
dada por:

 K12 
(
C2 ( x, t ) = C 01 t − t * ⋅  ) ( )
 ⋅ {exp [− Φ 1 ( x )] − exp [− Φ 2 ( x )]} + C01 t − t * ⋅ exp [− Φ 2 (x )] (24)
 K 22 − K11 

donde:

(
C01 t − t * = ) (
W1 t − t * ) (24)
QM

y t* es el tiempo de viaje.

CI41B INGENIERIA AMBIENTAL 10

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

EJEMPLO 1
Afloramiento Flujo Subterráneo

Supongamos un río en el cual su caudal se incrementa en forma exponencial debido a un flujo

subterráneo que aflora a lo largo del río, pero los otros parámetros no cambian
substancialmente.

Q( x ) = Q0 ⋅ exp (q ⋅ x )

Determine la variación de la concentración de un compuesto conservativo a lo largo del río.

EJEMPLO 2
Descarga hacia un Río

Una planta de tratamiento se encuentra descargando hacia un río un efluente con una elevada
concentración de bacterias coliformes. Aguas abajo del punto de descarga se encuentran
localizados una devolución de riego, la extracción para un canal de regadío y una zona de
recreación. En esta última es necesario cumplir con la normativa ambiental que indica que la
concentración de coliformes no debe exceder los 1000 coliformes por 100 ml.

Una vez en el río las bacterias coliformes son eliminadas a una tasa de primer orden (tasa de
consumo), mientras que su crecimiento ocurre a una tasa de orden cero:

dX
rC = = − k1 ⋅ X
dt CONSUMO

dX
rX = =b
dt CRECIMIENTO

donde k 1 y b son dos constantes.

Si cada tramo de río (ver Figura E2.1) puede ser modelado como un reactor de flujo pistón
determine la concentración de coliformes a la salida de un tramo de río en función de variables
como la concentración de entrada (XO), el volumen del tramo (V), el área de escurrimiento (A),
el caudal transportado (Q), entre otros.

Utilizando los datos entregados en la Tabla E2.1 determine la concentración de coliformes

esperada al inicio de la zona de recreación. ¿Es posible autorizar su funcionamiento? En caso
que su respuesta sea negativa, ¿Cuál sería el porcentaje de remoción de bacterias (mediante
desinfección) que debería utilizarse para lograr la autorización de la zona recreacional?

CI41B INGENIERIA AMBIENTAL 11

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005
Universidad de Chile
Departamento de Ingeniería Civil
Otoño 20012

Tabla E2.1
Resumen Información Básica
Parámetro Unidad Valor
QRIO m 3/s 10
XRIO Coliformes/100 ml 0
QPT m 3/día 2.000
XPT Coliformes/100 ml 10.000.000
QR m 3/día 500
XR Coliformes/100 ml 1.000
QCR m 3/día 2.500
k1 1/día 1.3
b Coliformes/100 ml/día 1.000

Tabla E2.2
Datos por Tramo de Río
Area Longitud Volumen
Tramo
(m2) (m) (m3)
1 3.0 700 2.100
2 3.5 1.000 3.500
3 3.2 500 1.600

Figura E2.1
Esquema de Sistema Analizado

CI41B INGENIERIA AMBIENTAL 12

TEMA 2.5 MODELOS DE CALIDAD DEL AGUA
SEMESTRE OTOÑO 2005