Modelación OD y DBO

Curso – Taller
Modelación de la calidad del agua en corrientes superficiales

utilizando el modelo QUAL2Kw
Tema 4.
Ing. Javier E. Holguín G. MSc. Dr.Sc.
Curso Mod. Calid. Agua con el QUAL2Kw Prof. Javier E. Holguin Tema 4. Mod. OD y DBO
Curso – Taller
Modelación de la calidad del agua en corrientes
superficiales utilizando el modelo QUAL2Kw
Tema 4. Modelación simplificada del OD y

la DBO en corrientes superficiales
Ing. Javier E. Holguín G. MSc. Dr.Sc.

e-mail: jaholgu@hotmail.com
1
Tema 4. Modelación simplificada del OD y la DBO en
corrientes superficiales
Contenido:
I. Modelación del OD y la DBO

1. Modelación del OD
2. Modelación de la DBO
3. Modelación conjunta del OD y la DBO
II. Streeter-Phelps multiples fuentes

puntuales
III. Modelo completo de Streeter-Phelps
Primera Parte
1. Modelación del
Oxígeno Disuelto
2
Máxima carga de DBO5 y
Oxígeno disuelto mínimo
Modelo de Calidad del Agua
C = f (W, física, química y biología)
Ciudad
Planta Tratamiento Planta Tratamiento

Carga removida
Agua Potable Aguas Residuales
Carga contaminante, W [M/T]
Río
Concentración crítica, C [M/L3]
Fuente: Camacho (2002).
Balance del
Reaireación
atmosférica Oxígeno
DBO
Carbonácea
N. Orgánico
O Resuspención
X Sedimentación
Amonio Í
G Demanda
E Béntica
Nitritos N
Fósforo
O Organico
Nitratos Respiración Fotosíntesis

Fósforo
Disuelto
Algas
Clorofila a
3
Saturación del Oxígeno en agua
Efecto de la Temperatura y Salinidad
10
Saturación del Oxígeno en agua

Efecto de la Presión (Elevación)
4
Efectos de Temperatura, Salinidad y Elevación
sobre la Saturación del Oxígeno
- Efecto de la Temperatura: Ta (K )  T ( o C )  273.15
1.575701* 10 5 6.642308* 10 7 1.2438* 1010 8.621949* 1011

ln Os,T  -139.34411    
(Ta ) (Ta )2 (Ta )3 (Ta )4
- Corrección por Salinidad (a mayor salinidad menor OD saturación)

 1.0754 * 101 2.1407 * 10 3 
ln Oss  ln OS,T  S1.7674 * 10 2   2


 Ta Ta 
- Corrección por altura: (a mayor altitud menor OD saturación)
O sf  O s 1  0.1148 * Alt (km)
FENÓMENO DE REAIREACIÓN EN
CORRIENTES SUPERFICIALES
Antecedentes
I. Streeter and Phelps en 1925
II. Método semi-empírico:
 O’Connor and Dobbins en 1956
Ux
Ka  c y
 Krenkel y Orlob en 1962 H
 Tsivoglou and Wallace en 1972
 Fok and Stall en 1972.
III. Métodos empíricos: H

 Churchill et al. en 1962 Ka  c
 Langbein and Durum en 1967 t
 Owens and Gibbs et al. en 1964
5
Antecedentes
IV. Reaireación medida directamente en campo:
1. Balance de oxígeno en estado estable

2. Desoxigenación con Sulfito de Sodio
3. Oscilación diurna del oxígeno
4. Inyección de sustancias volátiles (trazadores)
IV. Reaireación medida directamente en campo
4. Inyección de sustancias volátiles (trazadores):
Trazador R Referencia
Krypton 1.2 Tsivogou y Wallace 0.5

(1972)  DO 
Ka   2  Kg
Etileno 1.15 Rathbun et al. (1978) D 
 g
Propano 1.39 Rathbun et al. (1978)
Metil 1.4 Wilcox (1984) K a  RK g

cloruro
SF6 1.38 Canale et al. (1995)
6
Estimación de las Tasas de Reaireación
Opción 1. Método del área bajo la curva: uso de trazadores
1 M
kg  ln 1 K a  RK g
t 2  t1 M 2
Concentración
Distribución del Distribución del

trazador aguas arriba trazador aguas abajo
t1
t
t2
c0 Tiempo
Estimación de las Tasas de Reaireación

Opción 2. Método de la concentración máxima: uso de trazadores
Inyección puntual
____ Trazador Conservativo
de ambos
---- Trazador Gaseoso
trazadores
+ Centroide de la curva del
trazador conservativo  Cg 
 
 
tp2 1  Cc 1
tp1
K g  ln J
tp t p  Cg 
 
Concentración
Cp C 
 c 2
Cg
Cp
Ac
Ac
Ag Ag
Cg
K a  RK g
tiempo
tc1 tc
tc2
Fuente : KILPATRICK et al, 1989
7
Por que necesitamos usar las
Tasas de reaireación ?
dO2
 K a * (Cs  C )
dt
dO2
 tasa de intercambio (mg )
dt l *s
Ka  Coeficiente de reaireación(s-1 )
Cs  Concentración de saturación del oxígeno(mg/l)
C  Concentración de oxígeno disuelto (mg/l)
Como modelamos las Tasas de

Reaireación? Tres Ecuaciones Predictivas generales
Tasas de reaireación a una temperatura de 20 °C

Con las siguientes unidades U = [m/s] y H= [m]
U 0.67
Owens and Gibbs Ka  5.32
H 1.85
U
Churchill Ka  5.03 1.67
H
Oconor Dobbins U 0.5

Ka  3.93
H 1.5
8
TABLE 3-6. REAERATION COEFFICIENTS FOR RIVERS ANO STREAMS
Author(s) t 2, base e(1/d~ at 20°C) Un1ts Appltcabt 1tty
o•connor and Dobbtns (1958) U-fps Moderately deep to deep ch~nnels; 1~ttS30', 0.5SUS..6 fps;
H-feet 0.05Sk.st2.2/dl1. o•connor and Oobbtns also developed 1
secondffo~la for shallow stre..s but o•connor (1958)
showed the dtfferences bet~~~~en the two fo,.,las was
tnstgntf1cant, and recoaended that the ftrst fo~la be
used.
Churchtll .!! .!!· (1962) U-fps Based oa observed reaeratton rates below d.-s tro. Whtch
H-feet oxygen deftetent water •s released. 2'~11$11'¡ 1.8:91::!35 fps.
Church111 J!J!. also developed other fo,..las, but
reeo.endeTtlils foiWUla.
Owens .!! .!.!.· (1964) U-fps Oxygen reeovery •n1tored for six stre•s 1n England
H-feet follow1ng deoxygenat1on w1th sodtu. sulf1te.
o.~ fps: o.c•sHSJ.t•.
Owens .!! al. ( 1964) U-fps Thts 1s a second fo...,la developed by Owens et al., and
H-feet app11es for 0.1SU~.8 fps¡ 0.4'SH 11'.
__, Langbetn and o.m. (1967) 7.6U U-fps Basad on synthests of data fra. o•eonnor-Oobbtns (1958),
o
•(,.,..l ttr.D H-feet Church111 et al. (1962), Krenkel and Orlob (1963).,and
Streeter erar. (1936).
lsaacs and llud.Y (1968) 8.62U U-fps Developed using regression analyses fraa data collected
HD H-feet. us1ngA rectrculattng eyltndrteal tank.
o.&:s:USl.& fps; o.s·sHSt.5•.
Parkhurst II'ICI P~ro~ (1972) 48.4(1+0.17F2}(SU)3/8 U_.,s Oeveloped fra. data collected in 12 sewers lftd tn 111tural
A S-11/• stre•s.
H-meters
Negulescu lftd RojanSkt (1969) 10.9~t·85 U-fps Developed fra1 1 rectrculattng.flUie .tth depths le$s
H-feet than 0.5 feet.
Thackston and Krenkel (1969) 24.9~+~· 5)u. u.-fps

H-feet
Blsed on •uurelllftts ude tn a 2' wtde
deoxylenated waters • O.os•s H~0.23'.
n..- wtth
R
Lau (1972b) 2515/_uiu Based on reanalys1s of the data of Thackston and Krenkel
\UfA u.-fps
U -fps (1969}, Krentel (1960), a"d Church111 et !1· (1~2).
H-feet
(cont1nued)
''11
Ecuaciones Predictivas de las
Tasas de Reaireación
DETERMINACION EXPERIMENTAL
DE LA TASA DE REAIREACION CON TRAZADORES
CASO 1. UN RÍO DE MONTAÑA

CASO 2. RÍO CAUCA
25 6
Concentración Rodamina (ppb)
Sal
5
Concentración Sal (mg/l)
20 Rodamina
15
3
10
2
Fuente: Archivo personal

5
1
0 0
300
600
900
1200
1500
1800
2100
2400
2700
3000
3300
3600
3900
4200
4500
4800
5100
5400
0
Tiempo (S)
9
INVESTIGACIÓN TASA DE REAIREACION
OBJETIVO GENERAL
Determinar la tasa de reaireación

(Ka) en un río de montaña (Guavio) y
en un río de Planicie (Cauca) mediante
un estudio con trazadores.
OBJETIVOS ESPECIFICOS
Determinar cual es la ecuación predictiva
de las tasas de reaireación que más se
aproxima a los valores observados,
mediante comparación de los resultados
experimentales, con los calculados con las
ecuaciones predictivas
10
OBJETIVOS ESPECIFICOS
Presentar técnicas adecuadas de
utilización de trazadores e identificar las
potenciales fuentes de error, que
afectan la determinación de la tasa de
reaireación, mediante estudios con
trazadores.
METODOLOGIA
Selección del Río y tramo de estudio (río Montaña)
Río Guavio

11
METODOLOGIA
Selección del Río y tramo de estudio (río Planicie)
Río Cauca
Sector La Balsa
Fuente: Archivo Proyecto PMC
METODOLOGIA
Selección del Río y tramo de estudio (río Planicie)
Río Cauca
Sector Río Cali
Río Cali
Río Cauca
Fuente: http://www.publimetro.co/cali
12
METODOLOGIA
Diseño de experimentos
Manómetro Ent. agua
Ent. Sf6
Sf6(g)
Sf6(g)
Sf6(l) Rodamina
Sal. mezcla
Cilindro Sf6 Reactor Saturador
Fuente: Archivo Personal
METODOLOGIA
Diseño de experimentos (río Montaña)
13
METODOLOGIA
Diseño de experimentos (río Planicie)
Fuente: Archivo Proyecto PMC
METODOLOGIA
Recolección de muestras
14
METODOLOGIA
Recolección de muestras
METODOLOGIA
Análisis de muestras
15
METODOLOGIA – Río de Montaña

Modelación de los fenómenos de transporte de
solutos, Calibración y análisis de sensibilidad
a) Modelo Advección Dispersión (ADE)

b) Almacenamiento Temporal (TS)
Calibración mediante simulaciones de Montecarlo

Modelación de los fenómenos de reaireación
Método de la concentración máxima
Inyección puntual
____ Trazador Conservativo
de ambos
---- Trazador Gaseoso
trazadores
+ Centroide de la curva del
trazador conservativo  Cg 
 
 
tp2 1  Cc 1
tp1
K g  ln J
tp t p  Cg 
 
Concentración
Cp C 
 c 2
Cg
Cp
Ac
Ac 0.5
Ag
Cg Ag  DO 
Ka   2  Kg
D 
tc1 tc
tiempo  g
tc2
Fuente : KILPATRICK et al, 1989
16
RESULTADOS – Río de Montaña
Curvas de distribución de las concentraciones de ambos trazadores.
Sitio 2.
25 6
Concentración Rodamina (ppb)

Sal
5
Concentración Sal (mg/l)
20 Rodamina
4
15
3
10
2
5 1
0 0
0
300
600
900
1200
1500
1800
2100
2400
2700
3000
3300
3600
3900
4200
4500
4800
5100
5400
Tiempo (S)

Calibración y análisis de sensibilidad (Montecarlo)
Modelo TS Trazador Rodamina
2
17
Resultados Modelo TS - Trazador Rodamina
2

Selección de los Parámetros de Transporte
Exp A D As  AT U As/A
___ ___ _____ ______ __ ___ _____
m2 m2/s m2 d-1 m2 m/s
____________________________________________________________
No. 1 1.6 1.09 0.133 0.019 1.7 0.156 0.08
No. 2 1.5 1.04 0.190 0.011 1.7 0.146 0.13
____________________________________________________________
Parámetros hidráulicos utilizados para determinar Ka

Exp Q H W S H tc X
___ ___ ___ _____ ____ ____ _____
m3/s m m m/m m min m
___________________________________________________________
No. 1 0.23 0.19 8.5 0.04162 10 22.3 186.23
No. 2 0.22 0.17 8.5 0.04162 10 24.7 186.23
___________________________________________________________
18
Determinación de la Tasa Reaireación
0.5
 DO 
Ka   2  Kg Ka  1.37 Kg
D 
 g
Experimento T Kg Ka Ka 20°
____ ____ ___ _______
°C d-1 d-1 d-1
______________________________________________
No. 1 17.2 77.97 106.66 114
No. 2 16.3 72.73 99.49 108.65
_____________________________________________

Comparación de los resultados con la ecuaciones predictivas
Experimento No. 2 Experimento No. 3

03/01/03 04/01/03
Ecuación MME
Ka,obs =114 d-1 Ka,obs = 108.65 d-1
Tsivoglou y  H 
K a  0,1772   114.64 103.51 1.03
Neal (1976)  tc 
 H 
Grant (1978) K a  0 ,1748   113.09 102.11 1.04
 tc 
Shindala y  H 
K a  0 ,1748   113.09 102.11 1.04
Truax (1980)  tc 
Tsivoglou y  H 
K a  0,1573  101.77 91.89 1.15
Wallace (1972)  tc 
19
CONCLUSIONES
 Las Ka en el río de montaña, son mejor descritas por los modelos de
Disipación de Energía, ya que este tipo de modelos tienen mayor
habilidad de describir la reaireación en aguas turbulentas en
comparación con los modelos de Renovación Superficial (aguas
calmadas y de superficie libre).
 La ecuación predictiva de Ka que mejor se ajusta a los datos

observados en el río de montaña, fue la de Tsivoglou y Neal (1976.
 La ecuación predictiva de Ka que mejor se ajusta a los datos

observados en el río de planicie, fue la de O’Connor and Dobbins en 1956

Contenido:


puntuales
20
Primera Parte
Modelo de DBOc en un
reactor cerrado
L
dL
V  k dVL
dt
en t  0, L  0
L
L  L0e  k d t
L0
Disminución de k d
t
0
21
Descarga puntual a un río
en estado estable
Qwcw
Qr cr Qc0
x
0
Qw cw  Qr cr
c0 
Qw  Qr
Flujo a pistón en estado estable

0
c c  2c
 U  D 2  kc
t x x
x
0 L
22
Flujo a pistón en estado estable
Solución aguas abajo
c
0  U  kc
x
donde c = c0 en x = 0
c k
 c
x U
k 
  x
U 
c  c0 e
Carga de DBOc en un río
c Qw Lw  Qr Lr
c0 L0 
Qw  Qr
x
tx U
23
Carga de DBOc en un río y
tiempo de viaje
Debido a que la velocidad es constante, x = Ut,
donde, t = es el tiempo de viaje. Entonces,
k 
 x
L  L0 e U 
 L0 e  kt
C
C0 Decrecimiento k
Incremento U
x
txU

L L
 U  kr L
t x
kr  kd  ks
En estado estable: vs
ks 
L H
0  U  kr L
x
Condición de frontera:
Qw Lw  Qr Lr
L0 
Qw  Qr
Solución:
kr
 x
U
L  L0 e
24
kr
 x
L  L0 e U L  L0 e  k r t
c Qw Lw  Qr Lr
c0 L0 
Qw  Qr
x
t  x U
Contenido:
Modelación del oxígeno disuelto y la DBO
en corrientes superficiales

puntuales
25
Primera Parte
3. Modelación conjunta del

Oxígeno Disuelto y la DBO
Modelo de DBOc en un
reactor cerrado
L
dL
V  k dVL
dt
en t  0, L  0
L
L  L0e  k d t
L0
Disminución de k d
t
0
26
Modelo de Oxígeno en un
reactor cerrado
o dL
 k d L
L dt
do
 k d L
dt
Si L  L0 y o  o0 en t  0
L  L0 e  kd t

o  o0  L0 1  e k d t 
Modelo de Oxígeno / DBO

de un reactor abierto
Os
A que
velocidad
O ocurre el
proceso?
L
27
Modelo de transferencia de gas
de un reactor
d Os
A
D H

de un reactor
Balance de masas y Ley de Fick dc
V   JAc
dt
J  D
c dc  dc 
siendo
x
entonces V    D  Ac
dt  dx 
do DA
siendo
dc
dx
c  c1
 2
L
entonces V  os  o 
dt H
dc  do
dx  H
28
de un reactor
do DA
V  os  o 
dt H
D Coeficiente transferencia  L 
KL  
H Masa de oxígeno T 
 
do
V  K L Aos  o 
dt
do
 k a os  o 
dt
KL A Coeficiente 1
ka  de reaireación O2 T 
V  
Modelo de DBOc / OD de un
reactor abierto
L
dL
V   k dVL
dt
en t  0, L  0
L L  L0 e  k d t
L0
Disminución de k d
t
0
29
reactor abierto
dL
 k d L
dt
do
 k d L  k a os  o 
dt
Si L  L0 en t  0
L  L0 e k d t
Se define como déficit (D)
D  os  o dD do
 Integrando
Sustituyendo dt dt
dD k d L0
 k d L0 e  kd t  k a D D( t ) 
k a  kd

e kd t  e  kat 
dt
dD
 k a D  k d L0 e k d t
dt
reactor abierto
Streeter-Phelps Bolsa de Oxígeno
30
Contenido:


puntuales
II. Streeter-Phelps multiples

fuentes puntuales
31
Balance de Oxígeno Disuelto
Modelo Streeter-Phelps (Bolsa de Oxígeno)
Streeter-Phelps fuente puntual

D D dD do
 U  kd L  ka D recordar que 
t x dt dt
En estado estable sustituir la solución de DBO para obtener

k
D  rx
0  U  k d L0e U  k a D
x
La condición de frontera es:
Q o  Qr or D0  os0  o0
o0  w w ;
Qw  Qr
Y la solución es:
ka
 x k d L0   kUa x k
 r x
D( x )  D0e U
 e  e U 
with D( x )  os  o( x )
kr  ka  
 
k d L0
D( t )  D0e  k at 
kr  ka

e kat  e k r t 
32
Taller 2. Streeter-Phelps fuentes puntuales
Calcule el nivel del oxígeno disuelto en el río Negro en el tramo comprendido
entre la estación Puente Miravalle (K 0+0) hasta el Puente El Salto (K
51+875). La desembocadura del río Molinos se encuentra ubicada en el
kilómetro (K 5+625) y la del río San Pedro en el kilómetro (K 16+875). La
caracterización fisicoquímica de las descargas y las condiciones iniciales se
muestran en la Figura y las características hidráulicas de cada tramo se
presentan en el Cuadro.
Asuma ks= 0.1 d-1, kd = 0.25 d-1 a 20 ºC y Altitud = 2562 m.s.n.m

Río Molinos
Q=0.2 m3 /s
T=16.8 ºC
DBO5 = 83 mgl
Sitio Velocidad Profundidad OD= 0 mg/l
(ms) (m)
Puente Miravalle - Río Molinos 0.3 2.76
Río Molinos- Río San Pedro 0.3 2.78
Río San Pedro -Puente El Salto 0.34 3.55
Río Negro Pte. El Salto
Pte. Miravalle
Q=10 m3 /s Río San Pedro
T=15 ºC Q= 8 m3 /s
DBO5= 2 mgl T= 13.4 ºC
OD= 2.4 mg/l DBO5=36 mgl
OD= 6 mg/l
Taller 2 : Streeter-Phelps fuentes puntuales

I. Efectos de Temperatura y Elevación sobre el OS
1.575701* 10 5 6.642308* 10 7 1.2438* 1010 8.621949* 1011
ln Os,T  -139.34411    
(Ta ) (Ta )2 (Ta )3 (Ta )4
O sf  O s 1  0.1148 * Alt ( km) Ta (K )  T ( o C )  273.15
II. Efectos de la Temperatura en las tasas de Rx: k T   k 20  T  20

III. Modelo Streeter-Phelps:
kr
 x
U
L  L0 e
ka with D( x )  os  o( x )
 x k L   kUa x k
 r x
D( x )  D0 e U
 d 0 e  e U 
kr  ka   D0  os0  o0
 
k d L0
D( t )  D0 e  k at 
kr  ka

e kat  e k r t 
33
Taller 2 : Streeter-Phelps fuentes puntuales
Efectos de la Temperatura
Las tasas de reacción en los cuerpos de agua naturales se
incrementan con la temperatura (Ecuación de Arrhenius).
k T   k 20  T  20
Reacción 
Reaireación de oxígeno 1.024
Descomposición de DBO 1.047
Crecimiento de Phytoplankton 1.066
Biología “Regla de oro” 1.07
Demanda béntica de O2 1.08
Taller 2 : Streeter-Phelps fuente puntual

Ecuaciones Predictivas de las Tasas de Reaireación
34
Contenido:


puntuales
Streeter-Phelps – Carga distribuida
35
El “eslabón perdido” en el modelo
clásico de Streeter-Phelps
aire
ka
Materia org.
déficit
Amonio kd , kn
O2
(DBOC/DBON) agua
ks DOS
Materia ? déficit
O2 sedimentos
Org.
Modelo Mecanistico de
Demanda Bentica
COLUMNA AGUA
MOP FLUJO DOS, NH 4+, Y CH4

FLUJOS
ZONA AEROBICA
CAPA ACTIVA SEDIMENTOS
ZONA ANAEROBICA
PRODUCCION de
NH 4+, y CH 4
Chapra(1987)
36
Demanda de oxígeno de los
sedimentos (SOD)
dD S
0  U  ka D  B
dx H
Si D  0 en x  0,
k
S H   ax
D B 1  e U 
ka  
 
Demanda de oxígeno de los

sedimentos (SOD)
dD S
0  U  ka D  B
dx H
Si D  0 en x  0,
k
SB H   ax
D 1  e U 
ka  
 
D
x, t
37
Fotosíntesis (P) y respiración (R)
Variación de la fotosíntesis en el ciclo diurno debido
a la variación en la intensidad horaria de la luz solar
Fotosíntesis (P) y respiración (R)

Método de las botellas claras y oscuras
38
Streeter-Phelps – Modelo semi-completo
SL
L  L0 e  kr t 
kr

1  e kr t 
(puntual distribuido)
SB
PR
k d L0  k r t
D  D0 e  k a t 
ka  kr
e   e kat   ka

H 1  e kat 
déficit DBO déficit
puntual Punt. distribuido
kd S L kd S L

kr ka

1  e kat  
k r k a  k r 

e  k r t  e  k at 
(DBO distribuido)
CICLO DEL NITRÓGENO
Película-gas
controla Película-liq.
volatilización O2,s N2,s
controla
aire volatilización
agua
NH3 O2
alto
pH denitrificación
Bajo
amonificación pH nitrificación nitrificación
Norganico NH4+ NO2- NO3- N2
denitrificación
Producción neta
Fijación de nitrógeno
Fuente: adaptado de Chapra (1997)
39
PROCESO DE NITRIFICACIÓN
Agua residual Agua residual
nitrosomonas nitrobacterias
nitrógeno
amonio nitrito nitrato
orgánico
ammonificación oxígeno oxígeno
(hidrólisis)
NITRIFICACION
I. Modelo de nitrógeno total

DBOn
DBO DBOn = 4.57 NTK (Nitrógeno total Kjeldahl)
mg/L DBOn = 4.57 (Norg + Namoniacal)
DBOn
DBOc
tiempo
40
CARGA DE DBOn A UN RÍO
LN  4.57 * NTK (Nitrógeno total Kjeldahl)
kn x

U
LN  LN0 e
ka k
 x k n Lno   kUn x  ax
D  D0e U
 e  e U 
ka  kn 

k n : tasa de oxidación o descomposición de DBOn
0.1 – 0.5 d-1 para ríos profundos (e.g. río Cauca)
> 1 d-1 corrientes poco profundas (e.g. ríos de montaña)
I. Modelo de nitrógeno total DBOn

Deficiencias:
• Se agrega Nitrógeno orgánico y amonio en una sola entidad

(DBOn) y por lo tanto se ignora la tasa de conversión de uno a
otro (hidrólisis) y se omite el tiempo retraso.
• Si se utiliza una sola tasa de descomposición se omite el

tiempo de conversión nitrito – nitrato
• La perspectiva de DBOn no tiene en cuenta en forma directa

los cofactores requeridos para que ocurra la nitrificación
(presencia suficiente de bacterias (biopelicula en rocas),
valores de pH (optimo 8.5 unds) y suficiente oxígeno disuelto
(mínimo OD entre 1 y 2 mg/l)).
41
PROCESO DE NITRIFICACIÓN
Agua residual Agua residual
nitrosomonas nitrobacterias
nitrógeno
amonio nitrito nitrato
orgánico
ammonificación oxígeno oxígeno
(hidrólisis)
NITRIFICACION
II. MODELO DE ESPECIES DE

NITROGENO
N org NH 4 NO2 NO3

N0 N1 N2 N3
dN 0
  k01 N 0
dt
dN1
 k 01 N 0  k12 N1
dt
dN 2
 k12 N1  k 23 N 2
dt
dN 3
 k 23 N 2
dt
42
NITROGENO
N org NH 4 NO2 NO3
N0 N1 N2 N3
N o  N o 0 e  koat
k ai koa N o 0  koat
N a  N a 0 e  k ai t 
k ai  koa
e  e  kai t 
k ai N a 0  k ai t k k N  e  koat  e  kint e  k ait  e  kint 
Ni 
kin  k ai

e 
 e  kint  ai oa o 0
k ai  koa
  
 kin  koa kin  k ai 
Nn 

NITROGENO
43
INHIBICION DE LA NITRIFICACIÓN
f nitr  1  e  knitr o
f nitr 1,2
0,8
0,4
0
0 2 4 6 8 10
omg / L 
INHIBICION DE LA NITRIFICACIÓN
44
Modelo completo
de Streeter-Phelps
Modelo completo Streeter-Phelps
45
GRACIAS POR SU ATENCION
46

Modelación OD y DBO

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Modelación OD y DBO

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Modelación OD y DBO

Cargado por

Copyright:

Formatos disponibles

Curso – Taller

Modelación de la calidad del agua en corrientes superficiales

Ing. Javier E. Holguín G. MSc. Dr.Sc.

Tema 4. Modelación simplificada del OD y

Ing. Javier E. Holguín G. MSc. Dr.Sc.

I. Modelación del OD y la DBO

II. Streeter-Phelps multiples fuentes

III. Modelo completo de Streeter-Phelps

Planta Tratamiento Planta Tratamiento

Carga contaminante, W [M/T]

Nitratos Respiración Fotosíntesis

Saturación del Oxígeno en agua

- Efecto de la Temperatura: Ta (K )  T ( o C )  273.15

1.575701* 10 5 6.642308* 10 7 1.2438* 1010 8.621949* 1011

- Corrección por Salinidad (a mayor salinidad menor OD saturación)

- Corrección por altura: (a mayor altitud menor OD saturación)

O sf  O s 1  0.1148 * Alt (km)

III. Métodos empíricos: H

1. Balance de oxígeno en estado estable

Krypton 1.2 Tsivogou y Wallace 0.5

Metil 1.4 Wilcox (1984) K a  RK g

Distribución del Distribución del

Estimación de las Tasas de Reaireación

Como modelamos las Tasas de

Tasas de reaireación a una temperatura de 20 °C

Oconor Dobbins U 0.5

Thackston and Krenkel (1969) 24.9~+~· 5)u. u.-fps

CASO 1. UN RÍO DE MONTAÑA

Fuente: Archivo personal

Determinar la tasa de reaireación

INVESTIGACIÓN TASA DE REAIREACION

INVESTIGACIÓN TASA DE REAIREACION

Fuente: Archivo personal

Fuente: Archivo personal

Fuente: Archivo Proyecto PMC

INVESTIGACIÓN TASA DE REAIREACION

Manómetro Ent. agua

Cilindro Sf6 Reactor Saturador

Fuente: Archivo Personal

INVESTIGACIÓN TASA DE REAIREACION

Fuente: Archivo Personal

Fuente: Archivo Proyecto PMC

INVESTIGACIÓN TASA DE REAIREACION

Fuente: Archivo Personal

INVESTIGACIÓN TASA DE REAIREACION

METODOLOGIA – Río de Montaña

a) Modelo Advección Dispersión (ADE)

Calibración mediante simulaciones de Montecarlo

INVESTIGACIÓN TASA DE REAIREACION

Concentración Rodamina (ppb)

INVESTIGACIÓN TASA DE REAIREACION

INVESTIGACIÓN TASA DE REAIREACION

Parámetros hidráulicos utilizados para determinar Ka

INVESTIGACIÓN TASA DE REAIREACION

Experimento No. 2 Experimento No. 3

 La ecuación predictiva de Ka que mejor se ajusta a los datos

 La ecuación predictiva de Ka que mejor se ajusta a los datos

Tema 4. Modelación simplificada del OD y la DBO en