IA-5037 INGENIERIA DEL AGUA

Clase 8: TRANSPORTE DE MASA DE

CONTAMINANTES ENLOS ACUIFEROS

Ing. Mg.Sc. Ricardo Apaclla Nalvarte

El transporte de masa de contaminantes en los
acuíferos es consecuencia de tres procesos:

◼ Difusión molecular.
◼ Dispersión mecánica.

◼ Advección.
 DIFUSION MOLECULAR
Es un proceso mediante el cual un sólido se mueve en
el agua de las zonas de mayor concentración a las
zonas de menor concentración. El movimiento de
solutos está condicionado por el gradiente de
concentraciones. Está regulado por las leyes de Fick
adaptada a un medio poroso:

Primera Ley de Fick para régimen

C
permanente de concentraciones: F = −D *

x

Segunda Ley de Fick para C *  C

2
=D
régimen transitorio de concentraciones: t x 2
 C C *  C
2
F = −D*
=D
x t x 2

F= Flujo de masa por unidad de área y unidad de

tiempo, ML-2T-1
C=concentración de contaminantes, ML-3
t = Tiempo, T
D*=D.ω, coeficiente de difusión para medio poroso
L2T-1
ω =Coeficiente dependiente del inverso del cuadrado
de la tortuosidad.
También
D*=D.ω, coeficiente de difusión para medio poroso
L2T-1
1
= f 2
 

D=Coeficiente de difusión en agua libre, L2T-1

ω =Coeficiente dependiente del inverso del cuadrado
de la tortuosidad.
τ = Tortuosidad
τ = Tortuosidad
 Lt 
 =  
 Lr 
Lt= Longitud de la trayectoria entre dos puntos
Lr=Distancia en línea recta entre esos dos puntos.
El coeficiente de difusión, D* se determina en
laboratorio para cada caso concreto. Es directamente
proporcional al coeficiente de difusión molecular en
agua libre e inversamente proporcional al cuadrado
de la tortuosidad del medio.

Su valor suele estar entre 0.01 y 0.5.

 DISPERSION MECANICA

Es un proceso producido porque todas las partículas

de agua subterránea no se mueven uniformemente en
el interior del medio poroso.

Como consecuencia, en el frente de avance, se

produce una dispersión que implica una dilución del
contaminante.
El fenómeno tiene lugar tanto en la dirección principal
del flujo (dispersión longitudinal) como transversal
(dispersión transversal). Esta última se produce tanto
en la horizontal como en la vertical.
El fenómeno de la dispersión posee una naturaleza
probabilística, la elección del fluido de ir hacia la
izquierda o derecha cuando se encuentra con una
partícula de suelo conlleva a interpretación
probabilística.

Las causas que producen la dispersión mecánica son:

◼ Fricción de las partículas de agua con las paredes de

los poros: viajan más de prisa las partículas
sometidas a menor fricción.
◼ Tortuosidad de las trayectorias que condicionan el
tiempo de tránsito.
◼ Tamaño de los poros que condiciona la velocidad del
flujo.
El coeficiente de dispersión mecánica se expresa:

◼ Coeficiente de dispersión mecánica longitudinal:

DML =  i  vi

◼ Coeficiente de dispersión mecánica transversal:

DMT =  j  vi
α= Es la dispersividad en la dirección que indica el
subíndice, [L]
v = Es la velocidad de flujo en la dirección que indica
el subíndice, [LT-1]
La Difusión molecular y la Dispersión mecánica se
consideran conjuntamente bajo la denominación de
Dispersión Hidrodinámica.
Los coeficientes Dispersión Hidrodinámica, son:

Coeficiente de dispersión DHL =  L  vi + D*

hidrodinámica en la dirección
del flujo [L2T-1]

Coeficiente de dispersión
hidrodinámica en la DHT =  T  vi + D*
dirección transversal al flujo
[L2T-1]

Estos coeficientes han de ser determinados para

cada caso en ensayos de laboratorio.
Para saber si en la dispersión hidrodinámica
predomina la difusión molecular o la dispersión
mecánica, se utiliza un parámetro denominado
número de Peclet:
vL
P=
DHL
Donde:
V= velocidad real media del flujo, [LT-1]
L= longitud de la trayectoria recorrida, [L]
DHL= coeficiente de dispersión hidrodinámica en la
dirección del flujo.
Si P>6 predomina la dispersión mecánica
Si P<0.4 predomina la dispersión molecular.
Para valores intermedios no está claro el predominio
de uno u otro fenómeno.
 TRANSPORTE ADVECTIVO
Se realiza según las líneas de corriente del flujo
subterráneo. Su causa es el gradiente hidráulico y
está gobernado por la Ley de Darcy. El flujo de masa
vendría dado por:
dM
= C  dQ = C  K  i  dA
dt
dM=Masa de contaminante que es transportada en un
intervalo de tiempo dt, [MT-1]
C= Concentración del contaminante en el agua
subterránea, [ML-3]
dQ= Volumen de agua subterránea que circula a
través de una área dA en un dt, [L3T-1]
i= gradiente hidráulico
K= Conductividad hidráulica del acuífero, [LT-1]
Si se realiza un balance en un volumen elemental de
acuífero para un intervalo de tiempo considerado, en
el que no hay fuentes ni sumideros de masa, se llega
a la expresión:

( ) ( ) ( )
 v*C x  v*C y  v*C z C
+ + =
x y z t

V*= es la velocidad real del flujo subterráneo, [LT-1]

x,y,z= son las tres coordenadas correspondientes a
un sistema cartesiano, [L]
C = es la concentración de soluto en las aguas
subterráneas, [ML3]
t = es el tiempo
Si se supone la velocidad constante y que el flujo es
unidireccional, quedaría:

C C
= −v x
t x

Esta es la expresión más simple del transporte de

masa por advección.
 ECUACION GENERAL DEL TRANSPORTE DE MASA

La ecuación general de
C  2C  2C C
transporte de masa puede = DHL 2 + DHT 2 − vx
t x y x
expresarse como:

El miembro izquierdo da la evolución de la

concentración a lo largo del tiempo. En el miembro
derecho están los términos correspondientes a la
dispersión hidrodinámica longitudinal, a la dispersión
hidrodinámica transversal (en la horizontal) y al
transporte advectivo.

La ecuación puede complicarse más si el

contaminante reacciona con el terreno o con otros
elementos en solución.
 REACCIONES

El contaminante como elemento disuelto en el que se

producen en el medio hidrogeológico, agua
subterránea puede participar en reacciones químicas
tales como las reacciones que ocurren en la fase
líquida.
 REACCIONES

En el caso de los contaminantes son particularmente

interesantes las reacciones de:
• Sorción.-entre la fase sólida-fase líquida.

• Desintegración Radiactiva.

• Reacciones de Biodegradación.

Todas ellas llevan consigo la eliminación, al menos en

parte, del contaminante presente en la fase líquida.
 REACCIONES DE SORCION

Son reacciones por las que las sustancias disueltas en

las aguas subterráneas pasan a ser retenidas por el
terreno, disminuyendo por lo tanto su concentración
en la fase líquida. Hay varios tipos de sorción (Wood,
Kramer y Hem, 1990):

Adsorción
Quimisorción
Absorción:
Intercambio iónico
Adsorción: el soluto se fija en la superficie de la fase
sólida por la acción de las fuerzas electrostáticas.

Quimisorción: el soluto pasa a la fase sólida como

consecuencia de una reacción química.

Absorción: el soluto se incorpora a la estructura

mineral de la fase sólida

Intercambio iónico: aniones y cationes de la fase

líquida son intercambiados con los de la fase sólida
tendiendo a alcanzar el estado de equilibrio iónico.
El efecto de la sorción implica un retraso en el
movimiento del contaminante, denominado efecto de
retardo. Considerando este tipo de reacciones, la
ecuación general del transporte puede escribirse:

C  2C  2C C  a C *
= DHL 2 + DHT 2 − vx −
t x y x  t

El último término del miembro derecho de la ecuación

corresponde a la reacción de sorción.

ρa=densidad aparente del terreno, [ML-3]

θ = es la relación entre el volumen de agua contenido en el
terreno y el volumen del terreno (relación volumétrica). Si el
medio está saturado coincide con la porosidad efectiva.
C*=masa de contaminante sorbida por unidad de peso de la fase
sólida. Generalmente se mide en mg/kg.
Reacciones de sorción en estado de equilibrio.
Cuando las reacciones de sorción son rápidas con relación al
flujo subterráneo llega un momento en que cesa la
transferencia de soluto y puede asumirse un estado de
equilibrio entre la fase líquida y la fase sólida de su entorno.

Reacciones cinéticas
A veces las reacciones son más lentas y este equilibrio no
se alcanza. Son las reacciones en las que los procesos de
sorción (paso de soluto a la fase sólida) y de desorción
(paso de soluto de la fase sólida a la fase líquida) son
continuos y simultáneos.
A. Reacciones en estado de Equilibrio.

Existen dos modelos para definir la incorporación del

soluto a la fase sólida en las reacciones en estado de
equilibrio:

◼ Modelo de Freundlich

◼ Modelo de Legmuir
 Modelo de Freundlich

El Modelo de Freundlich viene dado por: C = KdC

* N

La ecuación general del transporte quedaría:

C  2C  2C
= DHL 2 + DHT 2 − vx
C  a  K d C N
−
( )
t x y x  t
 a K d N C N −1
El factor de retardo es: rf = 1 +


Kd y N son constantes que dependen de las características

propias del soluto y de la fase sólida y que es necesario
determinar en laboratorio.
Si se toman logaritmos a la expresión de Freundlich:

Log C * = log K d + N log C

Esta expresión representada en papel doble logarítmico C*

en ordenadas y C en abscisas, es una recta de pendiente N
y de ordenada en el origen log Kd

Un caso particular de este modelo se da cuando N=1,

entonces:
C* = Kd C

Representada en papel aritmético, C* en ordenadas y C en

abscisas, se obtiene una recta que pasa por el origen con
pendiente Kd.
 Este modelo se conoce como modelo lineal. Tiene la
ventaja de ser sencillo y el inconveniente de que supone
ningún límite a la capacidad de sorción de la fase sólida.

A Kd se le llama coeficiente
de reparto, indica la
relación entre la masa de
soluto que se incorpora a la
fase sólida y la que queda
en solución.

El factor de retardo en este caso, es: a Kd

rf = 1 +

 El modelo de Lagmuir

Considera que la fase sólida tiene una capacidad

limitada de sorción y que en ella pueden existir varias
zonas de sorción con propiedades diferentes.

Para una sola zona de sorción seria:

C 1 C C
= + o C =
*

C *
  1+  C

α =constante de sorción [L3M-1], generalmente en L/mg

β = cantidad máxima de soluto que puede incorporarse a
la unidad de peso de fase sólida, generalmente en
mg/kg.
La ecuación general del transporte quedaría en este caso:

 C 
  
C  2C C
2
C  a  1 +  C 
= DHL 2 + DHT − v −
t x y 2 x  t
x

El factor de retardo es: a 

rf = 1 +
 (1 +  C )2

Representando en papel con escala aritmética C/C* en

función de C se obtiene una recta de pendiente 1/β y
de ordenada en el origen 1/αβ

Si existen varias zonas con distintas propiedades de

sorción:
1 1 C  2  2 C
C =
*
+ +
1 + 1 C 1 +  2 C
 B. Reacciones en estado de no Equilibrio o Cinéticas

Las reacciones cinéticas describen la continua transferencia

de soluto desde la fase líquida a la fase sólida y viceversa, al
no llegar a alcanzarse un estado de equilibrio.

Este tipo de reacciones se ajustan a un determinado tipo

de modelos conceptuales, que, al no alcanzarse el estado
de equilibrio han de tener en cuenta el factor tiempo

• Modelo irreversible
• Modelo reversible
Modelo irreversible: asume que la reacción se verifica
exclusivamente incorporándose el soluto a la fase sólida. En
el límite, la reacción acabaría al desaparecer el soluto de la
fase líquida.

C *
= K1 C
t

K1 es la constante de degradación. La ecuación general del

transporte quedaría:

C  2C  2C C  a
= DHL 2 + DHT 2 − vx − K1C
t x y x 
Modelo reversible: supone transferencia de soluto entre la
fase sólida y la líquida en los dos sentidos y sin llegar a
alcanzar el equilibrio. Los más usuales son el lineal, el no
lineal y el bilineal.

Modelo Cinético Reversible Lineal

C *
= K 2 C + K 3C *
t

K2 es la constante que regula la transferencia de soluto

fase líquida a fase sólida.
K3 es la constante que regula la transferencia de soluto
fase sólida a fase líquida.
C *
Si se alcanza el equilibrio =0
t

se llega al modelo de Freundlich, para N=1

 Modelo Cinético Reversible no Lineal
C *
= K 5 C N + K 6C *
t

K5 es la constante de sorción y K6 la de desorción

Modelo Cinético Reversible Bilineal

C *
t
( )
= K 7 C  − C * + K 8C *

K7 es la constante de sorción y K8 la de desorción

β = máxima cantidad de soluto que puede ser sorbida por
la fase sólida

En los casos de modelos cinéticos, C *

reversibles e irreversibles, el término: t

se debe introducir en la ecuación general del transporte

para su resolución en cada caso.
 REACCIONES DE DESINTEGRACION RADIACTIVA
Los contaminantes radiactivos, además de participar
en las reacciones de sorción, desaparecen de la fase
sólida y de la fase líquida siguiendo la ley de
desintegración radiactiva.
A = A0 e − kt

A0= actividad inicial del radioisótopo

A= es la actividad después de un tiempo t
λ = vida media del elemento
ln 2
k= , cons tan te de de sin tegración


La desintegración radiactiva ha de añadirse como un

término más a la ecuación general del transporte de
masa.
 REACCIONES DE BIODEGRADACION

Afectan a las moléculas orgánicas disueltas en aguas

subterráneas. Estas moléculas quedan adheridas al
terreno formando una película que sirve de base para
el desarrollo de la actividad biológica.

Los microbios y bacterias se alimentan de las

moléculas orgánicas haciéndolas desaparecer del
sistema hidrogeológico. Si la actividad biológica
necesita del oxígeno para su desarrollo, al proceso se
le llama biodegradación aerobia. Si se desarrolla en
ausencia de oxígeno, biodegradación anaerobia.
Los procesos de biodegradación son complejos y
dependen de parámetros difíciles de evaluar:
concentración de microbios,
◼ tasa máxima de consumo de materia orgánica por
parte de los microbios,
◼ tasa de degradación microbiana, etc.

Debido a esto su formulación matemática presenta

dificultades, aunque se han propuesto ecuaciones al
respecto (Fetter, 1999).
 SOLUCION DE LA ECUACION GENERAL DEL
TRANSPORTE

Requiere la solución de cada uno de los términos que

la componen:
◼ advección,

◼ dispersión hidrodinámica y

◼ reacciones del contaminante

Teniendo en cuenta, además, las posibles fuentes y

sumideros de masa.

Como cualquier ecuación diferencial la obtención de

una solución particular exige imponer las condiciones
de contorno propias del problema que se aborda.
La solución de la ecuación del transporte de masa por
métodos analíticos se limita frecuentemente a casos
sencillos que suelen considerar medio homogéneo y
transporte unidimensional (Domenico y Schwartz,
1990; Zheng y Gordon 1995; Fetter, 1999).

Actualmente los problemas de transporte de masa

suelen abordarse por métodos numéricos que
permiten resolver problemas complejos para sistemas
uni, bi o tridimensionales. Se trata de los llamados
modelos digitales de transporte de masa.
El objetivo de estos modelos es el cálculo de la
distribución espacial y temporal de la concentración
de un contaminante introducido en un acuífero a
partir de un foco de inyección que puede ser puntual,
lineal o areal.

En esencia, los modelos digitales discretizan el medio

físico en una serie de elementos y el tiempo en una
serie de pasos.

Transforman la ecuación diferencial en un sistema

de ecuaciones cuyas incógnitas son las
concentraciones de contaminante en
determinados puntos para cada paso de tiempo.
El sistema de ecuaciones se resuelve a partir de un
algoritmo matemático o modelo numérico implícito o
explícito. El esquema de resolución adoptado puede
ser el de las diferencias finitas o el de los elementos
finitos.

Existen modelos digitales de transporte que se limitan

a resolver el término advectivo de la ecuación general
del transporte. Otros modelos resuelven la ecuación
de transporte completa (advección, dispersión
hidrodinámica, fuentes y/o sumideros y reacciones
químicas).
Las reacciones que actualmente pueden simularse
son las de sorción en equilibrio, las de sorción
cinéticas, las de desintegración radiactiva,
degradación biológica y algunas relacionadas con el
comportamiento de los hidrocarburos solubles que se
incorporan al flujo subterráneo.
 TRANSPORTE DE MASA EN LA ZONA NO SATURADA

En el transporte de masa en la zona no saturada

intervienen los mismos procesos que en el transporte
de masa en la zona saturada:
◼difusión molecular

◼dispersión mecánica y

◼advección.

La diferencia fundamental estriba en que en la zona

no saturada todos estos procesos están
condicionados por el contenido de agua en los poros.
En el caso de la difusión molecular, la
primera ley de Fick, que expresa el flujo de
masa a través de una unidad de superficie C
Fdif = D* ( )
en régimen permanente, se expresa: z

La segunda ley de Fick, para régimen transitorio de

concentración:
C   * C   2C C D* ( )
= D ( )  = D ( ) 2 +
*

t z  z  z z z

Fdif= Flujo de masa por unidad de área y de tiempo

D*(θ)=coeficiente de difusión molecular en la zona no
saturada. Depende de la tortuosidad, de las
fuerzas electrostáticas y de la humedad
volumétrica θ.
C =concentración de contaminante.
t = Tiempo
z = Expresa la dirección del movimiento del soluto.
La dispersión mecánica en la zona no saturada, además
de depender de la textura y estructura del terreno,
depende de θ que condiciona el espacio por el que se
mueve el agua y por tanto la velocidad de movimiento.

De esta manera se define un coeficiente de dispersión

mecánica Dm dependiente de los mismos factores que en el
caso de la zona saturada y de θ

El coeficiente de difusión molecular y el de dispersión

mecánica se integran en uno solo denominado coeficiente
de dispersión hidrodinámica, Dh. Considerando este
coeficiente, las leyes de Fick, para la zona no saturada,
expresan el transporte de masa como consecuencia de la
dispersión hidrodinámica.
 C
En régimen permanente: Fdh = Dh ( )
z

En régimen transitorio:

C   C   2C C Dh ( )
=  Dh ( )  = Dh ( ) 2 +
t z  z  z z z

El flujo de masa por unidad de superficie consecuencia del

transporte advectivo, será el producto del caudal
circulante por unidad de superficie (velocidad de flujo en
la zona no saturada, dependiente de θ) y la concentración
de masa existente en el agua, C:

Fadv = v( )C
El flujo total de masa por unidad de superficie será la
suma del debido al transporte por dispersión
hidrodinámicaa más el debido al transporte advectivo.

C
FT = Dh ( ) + v( )C
z

Si a partir de esta expresión, se realiza un balance de masa

en un elemento diferencial de acuífero de dimensiones dx,
dy y dz, para un tiempo dt, considerando que el transporte
de masa se realiza principalmente en la dirección vertical,
se llega a:

(C )   C  v( )C 
=  Dh ( )  +
t z  z  z

Esta ecuación expresa el principio de conservación de la

masa en la fase líquida del elemento considerado.
 (C )   C  v( )C 
=  Dh ( )  +
t z  z  z

El lado izquierdo es la variación de la masa que se

produce en la fase líquida en el intervalo de tiempo
considerado.

El lado derecho representa la diferencia de entrada y

salidas de masa en el elemento para dicho intervalo de
tiempo.
Si se considera que parte de la masa puede ser retenida
(o aportada) por la fase sólida del elemento y que
pueden existir fuentes o sumideros de masa, la ecuación
se puede completar como:

+
(
(C )   aC *)   C  v( )C 
=  Dh ( )  +  FS
t t z  z  z

ρa: densidad aparente de la fase sólida

C*: concentración de masa en dicha fase.
FS: entradas o salidas de masa por fuentes o sumideros.