Seminario: AQUASIM

Aplicación del Modelo de Fangos

Activos Nº 1 (ASM1)

1.- MODELO TRADICIONAL: ESTADO ESTACIONARIO

ACUMULACIÓN = ENTRA – SALE  REACCIONA

En estado estacionario: ACUMULACIÓN = 0

0  Q S 0  Qr S f   Q  Qr  S f  V
dS f
dt

S0 S0
Sf  
K1 V X K
1 1 1
Q S0 CM

1
  dX   dX 
0  Q X 0  Q X e  Qw X r  V   V ' 
 dt  REACTOR  dt  DEC RECIRC

Pf  Qw X r  Y Q S 0  S f   k d X V

NcOx

0  Q OD0  Qr OD f   Q  Qr  OD f  NcOx  V
dOD
dt

NcOx  a Q S0  S f   bV X

2.- Ventajas y limitaciones de la modelación dinámica

Predecir la calidad del efluente, la demanda de oxígeno y la

producción de lodos como respuesta a las fluctuaciones en tiempo
real de la carga y del caudal afluente.

El modelo calibrado de un proceso se puede utilizar para:

diagnóstico y evaluación de las prácticas actuales, proyecto,
comparación de variantes, pruebas de cambios operacionales y
optimización de la explotación, así como para capacitación personal.

En la fase de proyecto, se puede probar configuraciones

innovadoras y comparar variantes.

Se puede optimizar los ensayos a escala piloto (tiempo e

inversión requerida), realizando primero la aplicación de un
modelo pre-calibrado y sólo probar a escala piloto algunas de las
soluciones propuestas.

2
Ventajas y limitaciones de la modelación dinámica

Se puede optimizar la configuración de una EDAR, dando respuesta a

cuestiones sobre su capacidad y los impactos de modificaciones
proyectadas.

Herramienta de ayuda a la toma de decisión. Medio gráfico visual.

Se pueden realizar estudios sin perturbar el funcionamiento de la EDAR, ni

poner en riesgo sus equipos.

Herramienta de formación y de toma de conciencia de los operadores

sobre aspectos críticos (recirculación, TRC, DO, Pf, …).

Como principal limitación: necesidad de realizar estudios previos para pre-

calibrar los modelos. Los datos que se requieren, en general, difieren de
los datos rutinarios de caracterización que se tiene en las EDAR.

En el caso de aguas residuales industriales falta más investigación para

evaluar la aplicabilidad del modelo ASM1.

3.- Notación general del modelo ASM1

Matriz de un modelo que representa un crecimiento aerobio heterótrofo

Componente i  1 2 3 4 Cinética
Proceso j  XH SDQO SO2 SNH4
1 Crecimiento +1
1 1  YH - iXB
S DQO
heterotrofo   H XH
aerobio YH YH K DQO  S DQO
2 R. endógena -1 -1
SO 2
heterotrofa bH' XH
K O 2, H  S O 2
Parámetros Velocidad de reacción observada (M iL-3T-1) Parámetros cinéticos:
estequiométricos: ri  v ji Pj H, KDQO, b’H, KO2,H
Y j
Biomasa Sustrato Oxígeno Amonio
(g DQO/m3) (g DQO/m3) (g O2/m3) (g N/m3)

Se puede calcular la velocidad neta de reacción de cualquier componente i como

efecto de los diferentes procesos j: ri neta = vji pj

1  YH  H S DQO SO 2
rO 2   X H  bH' XH
YH K DQO  S DQO KO 2, H  SO 2

3
Funciones de activación y desactivación (ON/OFF)

Se toma como ejemplo la expresión cinética utilizada en el ASM1 para

el proceso de crecimiento heterotrofo en medio anóxico.

dX H SS K O 2, H S NO
  H , max g X H
dt K S  S S KO 2, H  SO 2 K NO  S NO

Donde:
XH: biomasa heterótrofa activa
SS: sustrato orgánico rápidamente degradable
SO2: concentración de oxígeno disuelto
SNO: concentración de nitrato
g: corrector de la velocidad de crecimiento en medio anóxico.
t : tiempo
KS, KO2H, KNO: constantes de semi-saturación o de activación /
desactivación

Caracterización del agua residual

• Biodegradable:
… rápida
… lenta (hidrólisis)

• Inerte

4
Compuestos (reactivos y productos) del ASM 1

DQO total

DQO DQO DQO

biodegradable inerte biomasa activa

SS XS XBH XBA
soluble en suspens. heterotrofa nitrificantes

SI XI XP
soluble inerte inerte producto

NTK SNO
nitrato + nitrito

SNH Nitrógeno XNB = iXB XH

amonio orgánico nitrógeno en la biomasa

N-Orgánico N-Orgánico
soluble en suspensión

No biodegradable No biodegradable
SNI Biodegradable XNI & XNP
Soluble en susp.
SND XND

Características de aguas residuales domésticas decantadas

Símbolo Unidades Dinamarca Suiza Hungría

3
SS g DQO/m 125 70 100
3
SI g DQO/m 40 25 30
XS g DQO/m3 250 100 150
3
XI g DQO/m 100 25 70
3
SND g N/m 8 5 10
3
XND g N/m 10 10 15
+ 3
SNH g N-NH4 /m 30 10 30
SNI g N/m3 2 2 3
- 3
SNO g N-NO3 /m 0.5 1 1

5
 Tabla.- Fracciones de la DQO en aguas residuales municipales

Fracciones (%)
fSI fXI fSS fXS fXH
Aguas Residuales Brutas
Sudáfrica 5 13 20 62 -
Suiza 11 11 32 45 -
Hungría 9 20 29 43 -
Dinamarca 8 19 24 49 -
Suiza
Flawill (22 ºC) 20 9 11 53 7
Tuffenwies (13 ºC) 10 8 7 60 15
Dietikon (15 ºC) 12 10 8 55 15
Dinamarca-Lundtofte 2 18 20 40 20
Aguas Residuales Primarias
España (*) 5 - 10 10 - 15 10 - 25 30 - 60 5 - 15
Dinamarca-Lundtofte 3 11 29 43 14
Suiza-Zurich 10 9 16 40 25
Francia-Pilot 10 13 33 44 -
Francia-Valenton 6 8 25 41 -
Sudáfrica 8 4 28 60 -
Aguas Residuales Pluviales
Dinamarca-Lundtofte 5 5 56 26 8

Composición media del agua residual decantada del área metropolitana de Coruña.

DQOT DQOS NTK NH4 NO3 SST Alcalinidad pH CE SSV P total DBOC
(mg/L) (mg/L) (mg N/L) (mg N/L) (mg N/L) (ppm) (mg/L CaCO3) (S/cm) (ppm) (mg P/L) (mg/L)
298 168 43 33 1,5 103 211 7,1 2931 75 9,8 250

Fracciones orgánicas del agua residual decantada Coruña

PARÁMETRO DESCRIPCIÓN VALOR MEDIO
fXI Fracción inerte en suspensión 0.03
fSI Fracción inerte soluble 0.13
fSS Fracción rápida degradación 0.47
fXS Fracción lenta degradación 0.37

Coeficientes cinéticos y estequiométrico a 20 ºC del agua residual decantada Coruña

PARÁMETRO DESCRIPCIÓN PROMEDIO RANGO
YH (g DQO/g DQO) Producción heterotrofa 0.73 0.5 – 0.84
H, max (d-1) Crecimiento específico máximo 5.9 4.5 – 7.1
de heterotrofas
bH (d-1) (*) Respiración endógena 0.92 0.67 – 1.1
de heterotrofas
A, max (d-1) Crecimiento específico máximo 0.42 0.31 – 0.53
de nitrificantes
(*) Modelo muerte - regeneración

6
 SO

SS XH XP
Crecimiento
heterótrofas
Hidrólisis Decaimiento
orgánica bacteriano
XS

Crecimiento
nitrificantes
SNH XA
Amonificación

SND SO SALC SNO

Hidrólisis
N-Org
XND Procesos aerobios
ASM1

SNO SALC

SS XH XP
Crecimiento
heterótrofas
Hidrólisis Decaimiento
orgánica bacteriano
XS

SNH XA
Amonificación

SND
Hidrólisis
N-Org
XND Procesos anóxicos
ASM1

7
 Matriz para oxidación carbonosa, nitrificación y desnitrificación

Componente i 1 2 3 4 5 6 7 8 9 10 11 12 13
SI SS XI XS XH XA XP SO2 SNO SNH4 SND XND SALK
Cinética, pj (ML-3T-1)
j Proceso
1 Crecimiento 1 +1 1 YH -iXB i SS SO
heterotrofo    XB  H ,max XH
YH
aerobio YH 14 K S  S S K O, H  SO

2 Crecimiento 1 +1

1 YH - iXB 1  YH i SS S NO K O, H
heterotrofo  2.86 YH  XB   XH
anóxico YH 14·2,86 YH 14 H , max K S  S S K NO  S NO K O, H  SO g

3 Crecimiento +1 4.57  YA 1 1 i XB 1
  i XB    S NH SO
nitrificantes YA YA YA  A, max XA
14 7 YA K NH  S NH KO, A  SO

4 Decaimiento 1 - fp -1 fp iXB – fp iXP

heterotrofo bH X H

5 Decaimiento 1 - fp -1 fp iXB – fp iXP

nitrificantes bA X A
6 Amonificación 1 -1 1
N-Org. soluble
ka SND XH
14
7 Hidrólisis de 1 -1 
XS XH SO  K O, H   S NO  
mat. orgánica kh   h    XH
K X   X S X H   K O, H  S O  
 K O, H  S O   K NO  S NO

 


8 Hidrólisis de 1 -1
p7 (XND/XS)
N-Org. susp.

DQO DQO DQO DQO DQO DQO DQO DQO N N N N mol

Expresiones cinéticas del ASM1

j Proceso Cinética, pj (ML-3T-1)

1 Crecimiento SS SO
heterótrofo  H ,max XH
K S  S S K O, H  SO
aerobio
2 Crecimiento SS S NO K O, H
heterótrofo   XH
H , max K S  S S K NO  S NO K O, H  SO g
anóxico
3 Crecimiento S NH SO
nitrificantes  A, max XA
K NH  S NH KO, A  SO
4 Decaimiento
H H b X
heterótrofo
5 Decaimiento
bA X A
nitrificantes
6 Amonificación
k S X
N-Org. soluble a ND H
7 Hidrólisis de 
XS XH SO  K O, H  S NO  
mat. orgánica k h   h    XH
K X  X S X H   K O, H  S O  K O,H  S O

  K NO  S NO

 
 

8 Hidrólisis de
p7 (XND/XS)
N-Org. susp.

8
 Origen y cambios que sufren los componentes

Componentes orgánicos del agua

Símbolo Descripción Columna

SI DQO inertes soluble 1

SS DQO rápidamente biodegradable 2

XI DQO inerte en suspensión 3

XS DQO lentamente biodegradable en suspensión 4

Columnas 1 y 3: SI y XI no participan en ninguna conversión, pero SI

contribuye a la concentración del efluente, mientras que XI es parte de
los sólidos volátiles del licor mezcla y de los lodos purgados.

Columna 2: SS se utiliza en el crecimiento de la biomasa heterótrofa

tanto en condiciones aerobias como anóxicas. Además de la cantidad
inicial presente en el afluente, SS se genera continuamente como
producto de la hidrólisis de la materia orgánica particulada XS.

Columna 4: XS se consume (transforma) por hidrólisis pero se produce

en el decaimiento de la biomasa heterótrofa y nitrificante.

Componentes orgánicos de la biomasa

Símbolo Descripción Columna

XH Biomasa heterótrofa 5

XA Biomasa nitrificante 6

XP Residuo inerte de la biomasa muerta 7

Se suele despreciar el aporte de biomasa activa por el afluente. Por otro

lado, la DQO total y los sólidos en suspensión volátiles (SSV) del lodo
purgado y licor mezcla están formados por biomasa activa y biomasa
inactiva: XH + XA + XP + XI + XS.

Columnas 5 y 6. XH y XA, se forman en el crecimiento aerobio, por un

lado, y anóxico (sólo XH) por otro, mientras que se consumen en el
proceso de decaimiento.

Columna 7. Al decaer, XH y XA se forma, entre otros productos, un

residuo particulado inerte XP que no participa en ninguna otra conversión
(contribuye a la biomasa inactiva).

9
Oxígeno disuelto

Columna 8. El oxígeno se consume en el crecimiento de heterótrofas y

nitrificantes.

El modelo considera que no se consume oxígeno directamente en el

decaimiento (pero si indirectamente a través de SS que procede de la lisis).

El modelo permite determinar el consumo de oxígeno (r O2), pero no las

variaciones en la concentración de oxígeno, SO. Para poder predecir SO, se
necesita incluir el proceso de aireación en el modelo.

Compuestos nitrogenados

Símbolo Descripción Columna

XND N-Orgánico en suspensión biodegradable 12
SND N-Orgánico soluble biodegradable 11
SNH N-amoniacal 10
SNO N-nitrato + nitrito 9

Además de las formas contempladas en la matriz, existen otras tres formas

de nitrógeno orgánico incorporadas en:

• la biomasa (XNB = iXB · XB),

• los productos inertes de lisis (XNP = iXP · XP) y
• la materia orgánica inerte (XNI = iXI · XI).

Estas formas se pueden calcular pero no son necesarias en el modelo

puesto que no se puede hacer un control del balance de nitrógeno, a
menos que se capture y cuantifique el nitrógeno gaseoso del proceso (N2).

10
Columna 9 (SNO). El segundo aceptor de electrones que incluye el modelo es
el nitrato (y nitrito). Se produce en el crecimiento aerobio de nitrificantes y se
consume (elimina) en el crecimiento anóxico de las heterótrofas.

El factor 2.86 que aparece en uno de los coeficientes estequiométricos

corresponde al equivalente en O2 de la conversión de NO3 en N2.
1 1
mol NO3  mol O2  1eq  e 
5 4

Columna 10. El nitrógeno amoniacal (SNH) utilizado procede, además de lo

que traiga el agua residual, de la reacción de amonificación del nitrógeno
orgánico soluble y se consume principalmente cuando sirve de fuente de
energía en el crecimiento de la biomasa nitrificante.

Hay que contar también una fracción (iXB) de N-NH4+ que se incorpora a las
células durante la síntesis. Se considera que a pH  7.0, SNH está formado
principalmente del ion amonio.

Como simplificación, el ASM 1 considera que la formación de NO3- se hace

directamente sin otra etapa que pase por NO2-.

Columna 11. El nitrógeno orgánico soluble (SND) es producto de la hidrólisis

de su forma particulada (XND), mientras que es transformado en amonio por
amonificación.

Columna 12. XND se pierde por hidrólisis, pero se reacumula con la muerte de
los microorganismos, tanto autótrofos como heterótrofos, donde se recupera
el equivalente de iXB – (fP · iXP).

El producto (fP · iXP) es la fracción que no es recuperable, porque está asociada

a la fracción inerte de los productos de lisis.

Alcalinidad
Columna 13. Representa el balance de alcalinidad (Salk). Se optó por hacer el
balance de alcalinidad, pudiendo obtener de éste, señales sobre el pH.

Para fines prácticos, se considera que una caída de 50 mg/L de CaCO3 en la

alcalinidad puede causar inestabilidad y caída del pH a valores menores que 6,
lo que podría afectar de forma significativa el proceso de nitrificación.

11
Utilización del modelo ASM 1
Para simular una EDAR en particular, se debe realizar una pre-calibración del
modelo, actividad que implica varias etapas.

Caracterización de las aguas y determinación de los parámetros del modelo

Los parámetros que se pueden tomar como supuestamente fijos son: YA, fP, iXB,
iXP, KO,H, KNO, KO,A, bA.

Símbolo Unidades Valor a 10 º C Valor a 20 ºC Rango a 20 ºC

Coeficientes estequiométricos

YA g DQOC/g Nox 0.24 0.24 0.07 – 0.28

YH g DQOC/g DQOS 0.67 0.67 0.46 – 0.69
fp adimensional 0.08 0.08 ---
iXB g NC/g DQOC 0.086 0.086 ---
iXP g Np/g DQOP 0.06 0.06 0.02 – 0.1

Parámetros cinéticos Valor 10 º C Valor 20 ºC Rango 20 ºC

H,max d-1 3.0 6.0 3.0 – 13.2
KS g DQO/m3 20.0 20.0 10 – 180
KO,H g O2 /m3 0.20 0.20 0.01 – 0.2
KNO g N-NO3 -/m3 0.50 0.50 ---
bH d-1 0.20 0.62 0.05 – 1.6
g adimensional 0.8 0.8 --
h adimensionmal 0.4 0.4 0.4 – 1
kh g DQO [Link]./g DQOC/d 1.0 3.0
KX g DQO [Link]./g DQOC 0.01 0.03
A,max d-1 0.3 0.8 0.34 – 0.8
KNH g N-NH4+/m3 1.0 1.0 --
KO,A g O2 /m3 0.4 0.4 --
ka m3/g DQO/d 0.04 0.08 --
bA d-1 0.15 --

12
Protocolo de simulación de plantas de tratamiento
Un ejemplo de protocolo que se podría seguir es el de STOWA (Fundación
holandesa de investigación aplicada al agua) (Hulsbeek et al., 2002).

Establecimiento de objetivos Sugieren calibrar el modelo

directamente en la planta,
Descripción de procesos
empezando con los valores por
defecto y ajustando algunos
parámetros específicos:
Recopilación y verificación de datos

a) Medir y calibrar la
Estructura del modelo composición y producción de
lodo, ajustando XS y XI del
Caracterización de flujos afluente, así como iXB y iXI.

Calibración b) Medir y calibrar el amonio

efluente, ajustando KO,H, KNH
y bA.
Caracterización detallada

Validación
c) Medir y calibrar el nitrato
Estudios de optimización efluente, ajustando g, bH,
KO,A, KO,H.

13