MODELIZACIN
AMBIENTAL
Modelado de calidad de agua: Modelo de corrientes
Introduccin
Modelo de calidad de agua de ros y
corrientes
Unidimensional
Mezclado perfecto vertical y lateral
Corriente principal y tributarios
Hidrulica en estado estacionario
Balance de calor y de cintica del agua
diarios
Entradas de masa y calor
Especial
Excel como interface
Especiacin de DBO orgnico
Interacciones agua-sedimento
Extincin de la luz, pH, patgenos
Presas, cascadas
Obtencin
[Link]
Q2KDocv2_07.pdf: Manual de usuario en versin pdf
Q2KFortran2_07.exe: ejecutable Fortran
Q2KMasterv2_07.xls: planilla Excel del modelo
Segmentacin
Ro: suma de tramos
Caractersticas
hidrulicas constantes
Numeracin desde
naciente a
desembocadura
Numeracin de afluentes
Tramo: dividido en
elementos
Resumen
Organizacin de la topologa del ro
Tramo (reach): longitud de ro con caractersticas
hidrulicas constantes
Elemento (element): unidad computacional
fundamental
Segmento (segment): conjunto de tramos que
constituye una rama del sistema
Cabecera (headwater): lmite superior de un
segmento
Planilla
Hojas azules: valores de variables y parmetros dados
por el usuario (datos de cabecera del ro, parmetros
hidrulicos y cinticos, parmetros de cargas puntuales
y difusas, etc.)
Hojas amarillas: entrada de datos de campo a lo largo
del ro. Se utilizan para comparacin con los valores
simulados, auxiliando la calibracin. Son graficados
junto a los simulados. No interfieren en los resultados
de la simulacin
Planilla
Hojas verdes: los valores simulados a lo largo del ro,
generados cuando se ejecuta el modelo (no deben ser
modificadas por el usuario);
Hojas rosa y hojas violeta presentan los resultados en
forma de grficos, que pueden ser formateados por el
usuario. Las curvas ilustradas son referenciadas de las
planillas amarillas (curva de datos medidos) y verdes
(curva de datos simulados).
Hidrulica
Balance de flujo en cada elemento
Velocidad y profundidad: vertederos,
Manning, cascadas.
Modelo de temperatura: atmsfera,
sedimentos, carga, descarga.
Balance de constituyentes: fuentes,
sumideros, sedimentos, algas, transferencia
atmosfrica.
Variables
Variable
Conductivity
Inorganic suspended
solids
Dissolved oxygen
Slowly reacting CBOD
Fast reacting CBOD
Organic nitrogen
Ammonia nitrogen
Nitrate nitrogen
Organic phosphorus
Inorganic phosphorus
Phytoplankton
Detritus
Pathogen
Alkalinity
Total inorganic carbon
Bottom algae biomass
Bottom algae nitrogen
Bottom algae
phosphorus
Smbolo
s
mi
Unidades
mhos
mgD/L
o
cs
cf
no
na
nn
po
pi
ap
mo
X
Alk
cT
ab
INb
IPb
mgO2/L
mgO2/L
mgO2/L
gN/L
gN/L
gN/L
gP/L
gP/L
gA/L
mgD/L
cfu/100 mL
mgCaCO3/L
mole/L
mgA/m2
mgN/m2
mgP/m2
cfu= unidades formadoras de colonias
mg/L g/m3
Reacciones
Bioqumicas: fotosntesis/respiracin,
nitrificacin/denitrificacin
Estequiometra de materia orgnica: generacin y
consumo de oxgeno
Efecto de temperatura en reacciones
Reacciones constitutivas: fitoplacton, algas de
fondo
Balance
Implementado para cada elemento del modelo:
Qi Qi 1 Qin ,i Qout ,i
Qin,i
i1
Qi1
Qout,i
Qi
i+1
Fuentes puntuales
Flujo total entrante a partir de fuentes
Qin,i
psi
npsi
j 1
j 1
Q ps,i, j Qnps,i, j
Flujo total saliente a partir de derivaciones
Qout,i
pai
pa ,i , j
j 1
npai
Q
j 1
13
npa ,i , j
Fuentes difusas
Las fuentes y derivaciones son modeladas como lneas (Flujo distribuido en
cada elemento ponderado por la longitud)
Qnpt
25%
25%
50%
start
end
Flujo
Calculado el flujo saliente de cada elemento, la
profundidad y velocidad se calcula por: vertedero,
curva h-q o Ecuacin de Manning.
Decisin:
Se implementa la opcin vertedero, si se ingresa altura y
ancho del vertedero
Si estas variables son nulas y se ingresan los coeficientes
(a y b) de una curva, se implementa la opcin curva h-q
Si no se satisfacen las condiciones anteriores, Q2K usa la
ecuacin de Manning
Manning
S 01 / 2 Ac5 / 3
Q
n P2/3
S0
B1
1
H
ss1
B0
ss2
1
Q, U
Cada elemento puede idealizarse como un canal trapezoidal. Bajo
condiciones de flujo estacionario, la ecuacin de Manning expresa la
relacin entre caudal y profundidad
MATERIAL
Man-made channels
Concrete
0.012
Gravel bottom with sides:
Concrete
0.020
mortared stone
0.023
Riprap
0.033
Natural stream channels
Clean, straight
0.025-0.04
Clean, winding and some weeds
0.03-0.05
Weeds and pools, winding
0.05
Mountain streams with boulders
0.04-0.10
Heavy brush, timber
0.05-0.20
Coeficiente de rugosidad de Manning para diferentes superficies (Chow et al. 1988)
Calor
Transferencia atmosfrica
Carga de calor
Flujo
entrante
dispersin
Salida de calor
Flujo
saliente
dispersin
Transferencia
sedimento-agua
Sedimento
Radiacin
El modelo calcula la cantidad de radiacin solar que entra al agua a una
latitud y longitud particular en la superficie terrestre (Lat y Llm), funcin de la
radiacin en el tope de la atmsfera atenuada por la transmisin atmosfrica, la
cobertura de nubes, reflexin y sombras
I (0)
I0
at
radiacin
atenuacin
extraterrestre
atmosfric a
ac
nubes
(1 Rs )
reflexin
(1 S f )
sombras
donde I(0) = radiacin solar en la superficie de agua [cal/cm2/d], I0 =
radiacin extraterrestre ( en el tope de la atmsfera) [cal/cm2/d], at = atenuacin
atmosfrica, ac = atenuacin por nubes, Rs = albedo (fraccin reflejada), y Sf =
sombra efectiva (fraccon bloqueada por la vegetacin y topografa).
Calor
Trminos
de no-radiacin
Trminos de radiacin
Radiacin
solar de
onda corta
Radiacin
atmosfrica
onda larga
Radiacin
onda larga
del agua
Radiacin neta absorbida
conduccin
y
conveccin
Interface
Aireagua
evaporacin
y
condensacin
trminos dependientes del agua
I(0) = radiacin solar neta de onda corta en la superficie del agua,
Jan =radiacin atmosfrica neta de onda larga,
Jbr = radiacin de onda larga desde la superficie del agua,
Jc = conduccin, and
Je = evaporacin.
Todos los flujos expresados en cal/cm2/d.
J I ( 0) J J J J e
h
an
br
c
Balance
Transferencia atmosfrica
Carga de masa
Sumidero de masa
Flujo entrante
Flujo
saliente
i
dispersin
dispersin
Algas de
fondo
sedimentos
Qout ,i
dci Qi 1
Qi
E i'1
E i'
W
ci 1 ci ci 1 ci i S i
c i 1
ci
ci
dt
Vi
Vi
Vi
Vi
Vi
Vi
QUALK2
Se elige el mtodo de solucin
general para las variables
(Solution method integration) y
el mtodo de solucin para el
pH (Solution method pH)
La opcin Final time define a
duracin de la simulacin, que
debe ser mayor o igual a dos
das, (tiempo mnimo necesario
para que el modelo alcance el
estado estacionario).
Afluentes
Headwater
Parmetros fsicoqumicos e hidrulicos
de cabecera del cuerpo
de agua superficial
principal, y de los
afluentes a ste.
Reach
Parmetros hidrulicos de los cursos. Se
representa un ro como una serie de
tramos de caractersticas hidrulicas
(So, B, n, z, etc.) y cinticas (kr, ks, ka,
etc.) constantes.
Los tramos se numeran en forma
ascendente comenzando desde la
cabecera
Las fuentes puntuales y no puntuales
pueden ubicarse en cualquier lugar a lo
largo del curso
Cuando se llega a una confluencia, la
numeracin contina en la cabecera del
tributario
El elemento computacional es la unidad
computacional fundamental del modelo,
que consiste en una subdivisin de
igual longitud en cada tramo
Reach rates
Tasas de reacciones
fsico-qumicas y
biolgicas simuladas
discriminadas por
tramo.
Rates
Valores globales de las tasas de
reacciones del cuerpo simulado (slo
aquellas que no varan a lo largo del
curso)
Modelo de reaereacin y coeficientes
de correccin de temperatura. Se
aconseja elegir el modelo Internal,
(seleccin automtica entre OConnorDobbins, Owens-Gibbs e Churchill de
acuerdo
a
las
caractersticas
hidrulicas de cada tramo.
Si se tienen valores calibrados de
reaireacin
los
valores
deben
insertarse en Reach Rates Prescribed
reareation , e indicar User specified
en la planilla azul Rates Reaeration
model
Sources
Se indica ubicacin de
la fuente y valores de
parmetros fsico
qumicos
WQ Data
Datos medidos a lo
largo del ro, para la
comparacin con los
simulados para la
calibracin del modelo.
Salidas
Oxgeno Disuelto
CBOD fast
