Revista Internacional de Contaminación Ambiental

Universidad Nacional Autónoma de México

rvp@atmosfera.unam.mx
ISSN (Versión impresa): 0188-4999
MÉXICO

2007
Franklin Torres Bejarano / Hermilo Ramírez León
EL PROGRAMA ANAITE PARA EL ESTUDIO DE LA DINÁMICA Y LA CALIDAD DEL
AGUA EN CAUCES NATURALES CON PENDIENTE SUAVE
Revista Internacional de Contaminación Ambiental, año/vol. 23, número 003
Universidad Nacional Autónoma de México
Distrito Federal, México
pp. 115-127

Red de Revistas Científicas de América Latina y el Caribe, España y Portugal

Universidad Autónoma del Estado de México

http://redalyc.uaemex.mx
 Rev. Int. Contam. Ambient. 23 (3) 115-127, 2007

El PROGRAMA ANAITE PARA EL ESTUDIO DE LA DINÁMICA Y LA CALIDAD DEL AGUA

EN CAUCES NATURALES CON PENDIENTE SUAVE

Franklin TORRES-BEJARANO1 y Hermilo RAMÍREZ-LEÓN2

1
Programa de Estudios de Posgrado. Instituto Mexicano del Petróleo. Eje Central Lázaro Cárdenas No. 152
México D.F. fmtorres@imp.mx
2
Instituto Mexicano del Petróleo. Programa de Investigación en Matemáticas Aplicadas y Computación. Eje
Central Lázaro Cárdenas No. 152 México D. F. hrleon@imp.mx

(Recibido enero 2007, aceptado junio 2007)

Palabras clave: hidrodinámica, transporte de contaminantes, modelación de la calidad del agua

RESUMEN

Se presenta la descripción, validación y aplicación de un modelo para el cálculo de

la calidad del agua en cauces naturales. La componente hidrodinámica se resuelve
en dos tiempos, primero se calcula el perfil hidráulico satisfaciendo la ecuación de la
energía a partir de una sección de control en donde se conoce el nivel y el caudal; en
segundo término, se calcula la hidrodinámica para el plano superficial, utilizando una
adecuación a la ley universal de velocidades para secciones irregulares. Esta compo-
nente hidrodinámica sirve de base para la solución del módulo de calidad del agua, el
cual está constituido por tres submódulos en las que se resuelven las variables físicas,
variables químicas y biológicas, y sustancias tóxicas, respectivamente. El transporte
de los parámetros modelados se realiza utilizando una solución numérica en diferen-
cias finitas de la ecuación de advección-difusión-reacción. Los diferentes módulos de
este programa se validan utilizando diferentes modelos existentes en la literatura. Los
resultados obtenidos con el modelo, llamado ANAITE, son satisfactorios.

Key words: hydrodynamic, pollutants transport, water quality modeling

ABSTRACT

In this paper, description, validation and application of a 2D water quality model for
natural rivers is presented. The hydrodynamic component is solved in two times, first,
the hydraulic profile is calculated satisfying the energy equation starting from a control
section where the level and the flow are know; then, the hydrodynamic is calculated
for the superficial plane, using an adaptation to the universal law of velocities for ir-
regular sections. The hydrodynamic component generated is the base for the solution
of the water quality module, constituted respectively by three sub-modules where
physical, chemical and biological, and toxicological substances variables are solved.
The transport of the modeling parameters is carried out using a numerical solution in
finite differences of advection-diffusion-reaction equation. The different modules of
this program are validated using different existent models in the literature. The results
obtained with the model, called ANAITE, are satisfactory.
116 F. Torres-Bejarano y H. Ramírez-León

INTRODUCCIÓN transporte de las sustancias disueltas y decaimiento

y transformación de las mismas a causa de los pro-
Un modelo de calidad del agua adecuado requie- cesos físicos, químicos, biológicos y toxicológicos.
re la especificación de una formulación apropiada de La representación funcional vincula la formulación
los procesos para tomar en cuenta aspectos del trans- de las características físicas, los procesos y condi-
porte longitudinal, lateral y vertical. La predicción ciones de frontera dentro del conjunto de ecuaciones
de la calidad del agua depende del procedimiento algebraicas y envuelve la definición precisa de cada
en el cual los procesos físico-químicos e hidrodiná- variable y su relación con los otros parámetros que
micos sean simulados (Maskell 1991, Calow 1994). caracterizan el modelo o las relaciones de entrada
Es importante que los métodos utilizados para re- y salida. En la figura 2 se presenta el diagrama de
presentar los diversos procesos, sean apropiados a bloques del modelo.
la aplicación del modelo. La finalidad de desarrollar
un modelo de calidad del agua es disponer de una
herramienta capaz de simular el comportamiento Lectura general
de datos
de los componentes hidrológicos y de calidad del
agua de un sistema de corrientes y realizar con ello Cálculo de los perfiles hidrodinámicos
estudios de diagnóstico y pronóstico del estado del Longitudinal y transversal

sistema en condiciones normales y por la presencia

de eventos extremos. Solución de los módulos de
calidad del agua
Shanahan et al. (1998) presentan un esquema de
t = t + ∆t
un modelo conceptual de ecosistemas de corrientes de n=n+1
T = Tmáx
ó Si
agua que consiste de elementos abióticos y bióticos Convergencia
No
enlazados dentro del ciclo hidrológico (Fig. 1). Los
procesos dentro del sistema y entre los elementos Si
Régimen hidrodinámico
permanente Resultados
son complejos y pueden ser descritos por una serie Salidas gráficas
de parámetros físico-químicos, hidro-morfológicos y
No
biológicos. Las estructuras abióticas y bióticas de la
corriente de agua son caracterizadas por gradientes
longitudinales, laterales, verticales y temporales. Fin

Fig. 2. Diagrama de bloques del programa ANAITE

Hidrología
gua Eco
de a siste
iente ma
Corr En los resultados se presenta la validación del
modelo con ejemplos de aplicación del programa
Morfología Comunidades Interfase ANAITE comparados con soluciones analíticas y
biológicas tierra-agua
con otros modelos de uso común por la comunidad
tecnológica. Se presenta la verificación del módulo
hidrodinámico longitudinal en el plano (x, z), la va-
Físico-químicas
lidación del módulo hidrodinámico bidimensional
en el plano (x, y) y la verificación del módulo de
Fig. 1. Estructura de un ecosistema de corriente de agua calidad del agua.
Entre las aplicaciones prácticas que se están
llevando a cabo con el programa ANAITE están el
El modelo de calidad del agua desarrollado está tramo final del río Magdalena, Colombia, (Torres-
dividido en tres partes conceptuales: un módulo hi- Bejarano 2003) y la desembocadura del río Pánuco,
drodinámico longitudinal (x, z) en el que se resuelven México (Islao et al. 2003).
las características hidráulicas básicas del cauce natu-
ral en estudio, proporcionando velocidades y niveles
del agua; un módulo hidrodinámico bidimensional EL MÓDULO HIDRODINÁMICO
que genera en el plano horizontal (x, y) los campos de
velocidades del río para secciones preestablecidas; y El remanso hidráulico, perfiles longitudinales x-z
un módulo de calidad del agua para la resolución del Básicamente, el método utilizado consiste en
 DINÁMICA Y CALIDAD DEL AGUA EN CAUCES NATURALES 117

1 2 vaciones-radio hidráulico y elevaciones perímetro

Línea de energía total mojado (Fig. 5). El proceso de solución fue tomado
2 de Berezowsky y Ramírez (1985).
V1 Sf
hf
2g La pérdida de energía entre dos secciones es
V1
∆ evaluada como:
2
V2 2
y1
2g  nV 
Q h f = ∆ x⋅ S f =  2 / 3  (2)
R 
V2
y2
donde:
Z1

∆x: ���������������������������������������
es la distancia entre dos secciones [L]
Z2
∆X n: es el coeficiente de Manning [L1/6]
Plano horizontal de comparación R: es el radio hidráulico [L]
Fig. 3. Definición de variables en un tramo de río Sf: pendiente de fricción []
El coeficiente de rugosidad de Manningn, requiere
siempre de un proceso de calibración específica y es
la aplicación de la ecuación de energía entre dos particular de cada cauce natural.
secciones, tal y como se ilustra en la figura 3 y se
demuestra a partir de la ecuación (1):
y
V2 V2 A
z1 + y1 + 1 = z 2 + y2 + 2 + hf (1) RH
2g 2g
P

y
donde:
A, RH, P

z: elevación sobre una referencia del Fig. 5. Curvas de elevaciones-área hidráulica, elevaciones-radio
punto más bajo de la sección [L] hidráulico y elevaciones-perímetro mojado, para una
y: tirante [L] sección cualquiera
hf: perdida de energía entre la sección 1 y 2 [L]
V: velocidad media en cada sección [LT -1] El programa es aplicado en régimen subcrítico
g: aceleración debida a la gravedad [LT -2] debido a que en general en el escurrimiento en cauces
con pendientes suaves se presenta esta condición y
Las secciones marcadas con 1 y 2 en la figura 3, para perfiles hidráulicos definidos como de tipo M
representan dos secciones cualesquiera de un tramo (Fig. 6). Un perfil M2, por ejemplo, ocurre cuando
en estudio. En realidad el sistema en estudio puede el fondo del cauce en el extremo de aguas abajo se
dividirse en varios tramos como se indica en la figura sumerge en un embalse con una profundidad menor
4. El proceso de solución de la ecuación (1) se lleva que la normal, o también en la desembocadura de un
a cabo satisfaciendo ambas partes de la ecuación, río cuando el nivel medio del mar es superior al de
para lo cual se requiere construir, para cada una de la superficie libre del agua del río. Los parámetros
las secciones, las curvas de elevaciones-gastos, ele- que arroja el modelo hidrodinámico longitudinal

Puntos marcados en
el campo

Secciones
Eje del cauce
del río
levantadas
en campo

a) Planta del río b) Sección transversal del río

Fig. 4. Esquema de las secciones transversales en un cauce natural
118 F. Torres-Bejarano y H. Ramírez-León

Nivel de marea alta

M1
Línea d
e tirante
s norma
les
Línea
de tira
ntes c M2
ríticos

Nivel de marea baja

Pendie
ntes s
ubcríti
cas

Fig. 6. Los perfiles típicos M1 y M2 de acuerdo a la condición de la marea en la

sección de control

son el ancho de la superficie libre, velocidad media, El campo de velocidades superficial se obtiene a
tirante, y radio hidráulico. En el cuadro I se presenta partir de la velocidad media obtenida con el modelo
un extracto de los resultados arrojados por el modelo hidrodinámico unidimensional (Fig. 7), a la que se
para la última sección de cálculo del río Magdalena; le aplica la ley de distribución de velocidades para
en forma similar se obtienen estos parámetros para canales naturales rugosos (Chow 1994). La fórmula
cada sección que se considere. utilizada para este cálculo es la siguiente:

Los perfiles horizontales, perfiles x-y  R 

La segunda parte del módulo hidrodinámico ge- V = U * ⋅  A 0 + 5.75log e  (3)
 k 
nera un campo hidrodinámico en el plano superficial
e interpola linealmente a un ∆Xij determinado, más gU
refinado e independiente de los ∆X considerados U* = (4)
en el cálculo del perfil hidráulico longitudinal. De C
igual forma calcula un ∆Yij con base en la variación
del ancho del canal a lo largo de la distancia. Este C = 28.6 Re (1 / 8 ) (5)
módulo tiene como función primordial determinar
la hidrodinámica bidimensional del río, de la que se donde:
adquiere el campo de velocidades superficial con el
que se realizará el transporte de los contaminantes. velocidad [LT -1]
V: �������������

CUADRO I. PARÁMETROS GEOMÉTRICOS E HIDRÁULICOS CAL-

CULADOS POR EL MÓDULO HIDRODINÁMICO UNIDI-
MENSIONAL

Sección No. 29
Elev. Área Perímetro Radio Hid. Ancho Caudal, Qc Factor de
(m) (m2) (m) (m) (m) (m3/s) pendiente
S

-1.00 13.85 81.09 0.17 81.08 17.93 0.372E-04

-0.50 92.79 261.12 0.36 261.10 173.25 0.312E-06
0.00 266.48 362.46 0.74 362.38 715.72 0.143E-07
0.50 451.80 379.07 1.19 378.92 1545.18 0.262E-08
1.00 645.40 395.68 1.63 395.47 2582.38 0.845E-09
1.50 847.27 412.29 2.06 412.02 3805.49 0.360E-09
. . . . . . .
. . . . . . .
. . . . . . .
7.00 3613.99 595.03 6.07 594.06 27919.06 0.467E-11
7.50 3915.15 611.64 6.40 610.61 31051.14 0.371E-11
8.00 4224.59 628.25 6.72 627.15 34341.95 0.298E-11
8.50 4542.31 644.86 7.04 643.70 37792.61 0.243E-11
9.00 4868.30 661.47 7.36 660.25 41404.30 0.199E-11
 DINÁMICA Y CALIDAD DEL AGUA EN CAUCES NATURALES 119

Eje transversal, en m
1400
1200
1000
800
600 800
400 700
200 Velocidad Media
600
20 25 30
Eje longitudinal en, m 500

400
0 10 20 30

Fig. 7. Generación del perfil hidrodinámico en los planos superficiales

k: altura de rugosidad de Nikuradse, U: velocidad media de la corriente [LT -1]

1.55 es el valor típico para ríos naturales [] B: ancho medio [L]
A0: factor de forma, para este caso se utilizó H: profundidad media [L]
un valor de 6.25 (Keulegan, 1938) [] U*: velocidad al cortante [LT -1]
U*: velocidad de fricción en las paredes [LT-1]
U: velocidad media del flujo [LT -1] La velocidad al cortante es evaluada con:
C: coeficiente de Chezy [L1/2T -1/4]
Re: número de Reynolds [] U * = gHS 0 (7)

Los coeficientes de dispersión longitudinal y en la que

transversal
Coeficiente de dispersión longitudinal g: aceleración debida a la gravedad [LT -2]
Para la evaluación de este coeficiente se utiliza el S0: pendiente del fondo del canal []
siguiente criterio (Fischer 1968):
U 2B2 Coeficiente de dispersión transversal
E x = 0.011 (6) Para este coeficiente, se utilizó la siguiente for-
HU *
mulación (Elder 1959):

donde: E y = φHU * (8)

Ex: coeficiexnte de dispersión longitudinal donde:

dirección x [L2T -1]

CUADRO II. EXTRACTO DE LOS DATOS DE ENTRADA PARA EL MODULO HIDRODINÁ-

MICO BIDIMENSIONAL

Secc Ancho Longitud Tirante Radio hidráulico Vel. Ex Ey

(m) (m) (m2) (m) (m/s) (m2/s) (m2/s)
29 660.25 546.00 9.75 6.886 0.68 218.0634 0.1875259
28 1400.10 695.00 9.78 8.150 0.26 215.2130 0.1875074
27 1673.03 558.00 9.78 8.287 0.21 212.4338 0.1874889
26 812.04 1648.00 9.91 4.921 1.33 209.7230 0.1874703
25 620.31 853.00 10.02 5.256 1.56 142.0012 0.1888295
24 580.20 543.00 9.14 4.982 1.05 140.2038 0.1887898
. . . . . . . .
. . . . . . . .
. . . . . . . .
6 544.49 617.00 11.95 6.964 0.88 142.0012 0.1886705
5 496.23 1061.00 11.99 9.333 0.66 140.2038 0.1826394
4 504.90 460.00 12.00 10.124 0.58 138.4559 0.1822683
3 561.73 458.00 12.01 9.555 0.54 136.7555 0.1818973
2 599.16 626.00 11.91 8.102 0.61 135.1008 0.1815267
1 736.05 661.00 11.92 7.519 0.52 107.8016 0.1811562
120 F. Torres-Bejarano y H. Ramírez-León

Ey: coeficiente de mezclado lateral [L2T] este término generalmente se anula en los ríos.
φ: coeficiente particular que depende de Las variables C, U, V, Ex, Ey y �
ΓC se pueden mo-
cada río, se utiliza un valor de 0.23 dificar con el tiempo y la ubicación, en el plano (x,
recomendado por Elder (1959) [] y). Esta ecuación está sujeta a varias condiciones de
frontera e iniciales cuya formulación varía depen-
En el cuadro II se muestra otro grupo de datos diendo de la aplicación. Las variables U, V, Ex y Ey,
que son generados por el modelo para cada una de son independientes de la concentración C y son deter-
las secciones consideradas; estos parámetros son minadas externamente por el módulo hidrodinámico
utilizados por los módulos de calidad del agua. para el modelo de calidad del agua.

LOS MÓDULOS DE CALIDAD DEL AGUA Solución numérica de la ecuación de advección-

difusión
Las partes constitutivas del modelo de calidad del La ecuación de advección-difusión que se resuel-
agua están divididas en tres módulos identificados y ve es la (9); esta ecuación es resuelta utilizando dife-
compuestos de la siguiente manera: rencias finitas centradas, lo cual facilita el tratamiento
numérico que se le aplica. En particular, el término
MODCA 1. En el que se resuelven las variables físi- netamente advectivo es resuelto en el programa por
cas, como la temperatura y la salinidad. dos métodos, el de diferencias centradas o bien por
MODCA 2. Resolución de las variables químicas y un esquema completamente euleriano. Ambos pro-
biológicas; OD, DBO, coliformes fecales, nitrógeno cedimientos se describen a continuación.
orgánico, fósforo orgánico, etc.
MODCA 3. Correspondiente a la resolución de Método de diferencias finitas con esquemas cen-
sustancias tóxicas, propiamente metales pesados; trados
zinc, hierro y manganeso. Aunque también se tiene Si se considera que en cada incremento de tiempo
la capacidad de modelar hidrocarburos. (∆t), las velocidades U y V y los coeficientes de disper-
sión Ex y Ey son estacionarios, se puede plantear una
Las ecuaciones diferenciales de cada uno de los solución general como se indica en la ecuación (10),
componentes modelados están escritas con base en sin considerar por ahora el término de decaimiento:
la ecuación de advección-difusión-reacción:
Cin, +j 1 = A1Cin, j + A2Cin+1, j + A3Cin−1, j +
∂C ∂C ∂C ∂  ∂ C  2
(10)
+U +V =  Ex + A4Cin, j +1 + A5Cin, j −1 + f ( ET )
∂t ∂x ∂y ∂x  ∂x 2 
(9)
∂  ∂ 2C  El término f(ET) representa el error de trunca-
 E y 2  ± ΓC
∂y  ∂y  miento que contiene esencialmente errores de aproxi-
mación de alto orden; se ha comprobado el aporte de
en donde: este término a través de diferentes trabajos (Ekebjaerg
y Justesen 1991, Laurel 2004), en donde se verifica el
C: representa alguna variable de calidad nivel de estabilidad de estos términos. Particularmen-
del agua: e.g. OD, DBO, nitrógeno te, en Laurel (2004) se hace una revisión utilizando
total, etc. [ML-3] diferentes esquemas de solución fundamentalmente
U: velocidad media de la corriente para los términos advectivos.
longitudinal [LT -1] Si se consideran diferencias centradas para aproxi-
V: velocidad media de la corriente mar tanto los términos advectivos como los términos
transversal [LT -1] difusivos, se pueden encontrar los siguientes valores
Ex: coeficiente de dispersión longitudinal [L2T -1] de los coeficientes Ai:
Ey: coeficientes de dispersión transversal [L2T -1]
ΓC: los mecanismos de reacción de las Cr − x
variables de calidad del agua [ML-3T -1] A1 = 1 − 2λx − 2λ y ; A2 = λx + ;
2
Cr − x C
El término �(V ∂ C/∂ y)���������������������������
se incluye cuando las des- A3 = λx − ; A4 = λ y + r − y ; (11)
2 2
cargas laterales lleven una velocidad importante con Cr − y
respecto a la de la corriente principal; de otra forma A5 = λ y −
2
 DINÁMICA Y CALIDAD DEL AGUA EN CAUCES NATURALES 121

Los coeficientes Ai son funciones de números j+1

característicos de estabilidad numérica, que para la
dirección en x, por ejemplo, se escriben como:
Fy+(i,j)
U∆t
Cr − x = , número de Courant ;
∆x
∆tE x Cr − x Fx+(i,j)
λx = = , número difusivo ;
∆ x2 Pe− x (12) j
U∆t
Pe− x = , número de Peclet
Ex

El término que incluye los errores de truncamien-

to, resulta:

 ∆t  ∂ C  ∆t  ∂ C
2 2 3 j+1
ET =   2 +   3 + i-1 i i+1
 2  ∂t  6  ∂t
(13) Fig. 8. Definición del esquema numérico
 ∆ x 2  ∂ 3C  ∆y 2  ∂ 3C
U  3 +V   3 + TAO
 6  ∂x  6  ∂y
C p = Ci−a = Ci−s − p(Ci−s − Ci−s−1 ) (15)
En donde TAO se refiere a términos de trunca-
miento de alto orden, que no son considerados en donde a es el número de Courant en la dirección x, s
estos desarrollos. es la parte entera y p la fraccionaria del mismo, por
Sustituyendo las ecuaciones (11), (12) y (13) en la lo tanto a = s + p donde s es el número de nodos que
(10), puede obtenerse la formulación siguiente: recorre la curva característica y p es la fracción de
∆x entre los nodos (i-s) e (i-s-1).
Cin, +j 1 = Cin, j + {Fxi−, j + Fxi+, j + Fyi−, j + Fyi−, j } (14) Generalizando para el caso bidimensional, pero
ahora considerando la ecuación de advección pura,
En donde Fxi,j = [Cr-x,y, λx,y, y las concentraciones C podemos escribir la siguiente formulación:
alrededor del nodo (i, j)]. Los signos (-) se refieren a
la evaluación de la función antes del nodo (i, j) y el (+) [
Fi−a , j −b = (1 − r ) (1 − p)(Fi−l , j −m + Fi−1, j −m−1 ) +
(16)
] [ ]
después del nodo (Fig. 8). Las funciones F incluyen la
optimización de los términos en la ecuación (13). pFi−l−1, j −m + r (1 − p) Fi−1, j −m + pFi−l−1,i−m

Solución del término advectivo por el método donde (l, m) son las partes enteras de (a,b), res-
euleriano-lagrangeano pectivamente y (p,q) sus correspondientes partes
Otra implementación que tiene el modelo ANAITE decimales, de tal manera que a = l + p; b = m + q.
es la solución de los términos advectivos a través de La formulación F es válida para la parte advectiva
una formulación euleriano-lagrangeana. Este método de la ecuación (9), de tal forma que:
consiste en desplazar una partícula desde un punto P
viajando a una velocidad constante U (tomando el caso ∂C ∂C ∂C
unidimensional como ejemplo) y llegando a una nueva F= = −U −V (17)
∂t ∂x ∂y
posición (x + U ∆t) en un tiempo (t + ∆t), generalizan-
do este concepto para cualquier número de Courant.
Siendo p un factor de peso entre las concentraciones Con esto, la ecuación (9) es resuelta en dos
más próximas al punto P, es necesario adicionarle otra tiempos, primero para la parte advectiva (16) y en
variable que represente el número de nodos hacia atrás segunda instancia la parte difusiva-reactiva, para cada
en el que se encuentra el segmento de interpolación, ∆t (Laurel 2004).
designado en este caso como s.
Así, una concentración que incluye un esquema Requerimientos de estabilidad
flexible en cuanto a la restricción de Courant-Frede- Para garantizar la estabilidad de la solución numé-
rich-Lewy (Rodríguez 2003) para la interpolación rica es necesario cumplir con los siguientes criterios
lineal se escribe como: Abbot (1979):
122 F. Torres-Bejarano y H. Ramírez-León

6
1 1
Cγ < Pe ; Pe < 2 ; λ <
2 2 5 Anaite
HEC-RAS
4

Vel, (m/s)
En la práctica suelen encontrarse recomendacio-
3
nes usuales que establecen que:
2
0.5 ≤ Cr < 1 ; Pe ≈ 1.5 ; 1 / 6 < λ < 1 / 2
1

0
Con lo cual deben de definirse los tamaños de ∆t, 0.000 5.000 10.000 15.000 20.000 25.000
∆x, tomando además en cuenta el valor del coeficiente Distancia, (m)

de dispersión, el cual a su vez, es definido por las Fig. 9. Comparación de velocidades obtenidas con el modelo
características del flujo, tal y como se especificó en ANAITE y el modelo HEC-RAS
la sección 2.3.
del modelo nos garantiza la versatilidad del mismo.
El mecanismo de reacción En particular se hizo énfasis en las velocidades, ya
A la resolución del transporte de sustancias me- que éste parámetro es el que define el transporte
diante el método previamente descrito, se le sumó de los contaminantes a lo largo del río y es la base
la capacidad de analizar las interacciones del medio para el desarrollo y funcionamiento apropiado del
con los contaminantes y de diversos contaminantes módulo hidrodinámico bidimensional, así como de
entre sí. Esta capacidad está representada dentro de los módulos de calidad del agua.
la ecuación de transporte presentada en la sección an-
terior (ec. 9), por el término ΓC. Cada mecanismo de El modulo para el transporte de escalares
reacción es específico para cada una de las variables Se resuelve la ecuación de transporte de escalares
que cambian debido a procesos de conversión físicos, para lo cual se seleccionaron los siguientes ejemplos
químicos y biológicos en función de la concentración que se mencionan a continuación.
de cada variable (Torres-Bejarano 2003). Una primera verificación de la solución del mo-
delo es hecha considerando una aproximación con
Pruebas de calibración y ajuste la solución analítica de la ecuación (9) en su versión
El cálculo del perfil hidrodinámico longitudinal unidimensional, como sigue:
Para la validación del módulo hidrodinámico
unidimensional se llevó a cabo una comparación de ∂C ∂C ∂ 2C
los resultados obtenidos con el ANAITE y con los +U = Ex 2 (18)
∂t ∂x ∂x
obtenidos de una simulación realizada con el soft-
ware HEC-RAS (River Analysis system) de la U.S.
Army Corps of Engineers. Las simulaciones fueron Cuya solución analítica puede escribirse como
aplicadas a los últimos 21 Km del río Magdalena, (Fetter 1992):
Colombia (Torres-Bejarano 2003). En este tramo
C0   L −U ⋅t 
 + exp U ⋅ L  x
el río presenta una hidrodinámica compleja, ya que C= erfc   E 
está influenciado por las descargas de los arroyos 2   2 E ⋅t   x 
  x 
de la ciudad de Barranquilla y por el gran numero (19)
de industrias ubicadas en su ribera, además de la  L + U ⋅ t 
erfc   
influencia ocasionada por la marea del mar Caribe.  2 E ⋅ t 
 x 
Para este tramo se disponían de los levantamientos de
29 secciones transversales, se calibró el coeficiente de Los datos considerados son: velocidad del flujo, U
rugosidad de Manning para cada tramo en estudio, la = 0.5 m/s, coeficiente de dispersión, Ex = 50 m2/s, y
pendiente del fondo fue considerada como S0 = 4x10-5 longitud del cauce de L = 2,000 m, con una concentra-
constante en todo el tramo de estudio. Se eligió un ción inicial, C0 = 100 mg/L en X0 = 0 y para t ≥ t0. Se
gasto de 3,000 m3/s, y una elevación conocida en la trata de calcular la dispersión del contaminante en el
última sección aguas abajo de 8.0 m. tiempo. Los resultados de concentración obtenidos con
En la figura 9 se presentan los resultados obte- el ANAITE y la solución analítica, en toda la longitud
nidos para las velocidades promedio en cada una de del cauce se muestran en la figura 10. Asimismo, se
las secciones. Como puede verse, los resultados son compara la evolución de la concentración en el tiempo
bastante satisfactorios por lo que esta primera etapa en un punto fijo del tramo estudiado (Fig. 11).
 DINÁMICA Y CALIDAD DEL AGUA EN CAUCES NATURALES 123

Concentraciones (mg/L) Concentraciones en (µg/L)

120 120
110 110
Ec_Fetter (1992)
100 Numérica 100 C inicial
C 4100-Medida
90 90 C 4000-Calculada
80 C 4400-Calculada
80
70
70
60
50 60
40 50
30 40
20 30
0 500 1000 1500 2000
Distancia (m) 20
10
Fig. 10. Comparación de las concentraciones obtenidas con el 0
0 10000 20000 30000
ANAITE y la solución analítica a lo largo del cauce Tiempo, en (s)

Fig. 12. Comparación de concentraciones medidas con las cal-

Como se aprecia, los resultados muestran un com- culadas por ANAITE en su versión unidimensional
portamiento y un decaimiento de la concentración
semejante para ambas soluciones tanto en la distancia resultados obtenidos muestran buena concordancia
como en el tiempo. con los resultados medidos en la sección 4100.
El módulo MODCA-2, para el estudio de la calidad
Concentraciones (mg/L) del agua
120 Se presenta un ejemplo de aplicación y calibra-
100 ción del modelo de calidad de agua (MODCA-2)
80
en su versión unidimensional, comparado con una
60
Ec_Fetter (1992) simulación realizada con el modelo QUAL2K para
Numérica

40
las mismas condiciones.
QUAL2K v 2.04 (o Q2K) es un modelo de calidad
20
del agua para ríos y corrientes desarrollado por la En-
0
0 500 1000 1500 2000 2500 3000 3500 vironmental Protection Agency (Chapra et al. 2006)
Tiempo (s)
como una versión modernizada del bien conocido
Fig. 11. Comparación de las concentraciones obtenidas con el modelo QUAL2E (Brown and Barnwell 1987).
ANAITE y la solución analítica en el tiempo Así como QUAL2E, QUAL2K simula el trans-
porte y destino de contaminantes convencionales. El
modelo representa el río como un canal unidimen-
Se reproduce ahora un problema planteado en Mc- sional con flujo no uniforme en estado fijo y simula
Cutcheon (1989) –ejemplo 3.2– en donde se tienen el impacto de descargas puntuales y difusas. Incluye
mediciones en un río localizadas a 2,400 y 4,100 m de el transporte y destino de constituyentes tales como,
un cierto punto de emisión las mediciones se mues- temperatura, demanda bioquímica de oxígeno, oxí-
tran con línea punteada en la figura 12. La velocidad geno disuelto, nutrientes, etc.
del río es constante e igual a 0.586 m/s. Los parámetros que se modelaron son la demanda
En este caso se procedió a resolver la ecuación bioquímica de oxígeno (DBO) y nitrógeno orgánico
genérica (10), para lo cual se eligieron los siguientes (NO). El ejercicio considera un tramo de río con
parámetros, dadas las características del problema: longitud total de 4 km, que sufre una descarga con-
∆x = 400 m y ∆t = 600 s y Ex = 130 m2/s. Con estos tinua de contaminante en el principio del tramo, con
datos se procedió a satisfacer los requerimientos de las características que se indican en el cuadro III.
estabilidad:
Cr = 0.88 ; Pe = 0.525 ; λ = 0.485 CUADRO III. CONDICIONES DE DESCARGA

Los resultados obtenidos para este problema Parámetro Unidad Valor

también se presentan en la figura 12, en donde al U m/s 2.0
utilizar un incremento en x de 400 m, no se ajusta Ex m2/s 150
totalmente a la curva medida en la sección 4100, por DBO mg/L 13
ello se muestran dos curvas obtenidas por el programa N org mg/L 280
P org mg/L 16
ANAITE, una a los 4000 m y otra a los 4400 m. Los
124 F. Torres-Bejarano y H. Ramírez-León

Se trata de predecir con ambos modelos el efecto de y la salinidad. Se eligió un dominio de estudio de
la descarga del efluente a lo largo del cauce para un 1,600 m de ancho constante y de 20 km de longitud
tiempo de simulación de 2 días: (Fig. 15).
Cr = 0.60 ; λ = 0.45 ; Pe = 1.33
Dirección transversal, en (m)
1600
La revisión de los criterios de estabilidad dados Margen izquierda
1400
anteriormente se establecen como óptimos a un ∆t = 1200
30 s y un ∆x = 100 m, con lo cual se tiene que satis- 1000
Eje central
facen los requerimientos de estabilidad. 800
Los resultados obtenidos para este problema se 600
presentan en las figuras 13 y 14, en donde aparece 400

también la solución obtenida con el QUAL2K, para 200 Margen derecha

0
las mismas condiciones de simulación. Como se 0 5000 10000 15000 20000
Dirección longitudinal, en (m)
aprecia, estos resultados muestran una concordan-
cia aceptable entre ambos modelos. Los puntos de Fig. 15. Dominio de estudio y vectores de velocidad conside-
máxima concentración coinciden aproximadamente rados
al 99 %, en tiempo e intensidad.

En la primera simulación se presenta una descarga

Concentraciones (mg/L) puntual y al centro del dominio con las siguientes
14
características: temperatura total de descarga igual
12 Qual2k
Anaite
a 30 °C y en el centro de la frontera inicial del do-
10
minio, la velocidad media del flujo, U = 5.0 m/s, los
8
coeficientes de dispersión longitudinal y lateral, Ex =
6
150 m2/s y Ey = 10 m2/s, respectivamente.
4

2
La temperatura natural del agua establecida para
0
este ejemplo es de 20 °C. La temperatura máxima se
0 1000 2000 3000 4000
Distancia (m) encuentra elementalmente en el punto de la descar-
ga, con una rápida disminución en la dirección de
Fig. 13. Resultados de DBO obtenidos con ANAITE y
QUAL2K la corriente, predominando en la mayor parte del
dominio una temperatura que oscila entre 20.5 a
21.4 °C para la figura 16a.
Concentraciones (mg/L)
300
En la figura 16b se puede confirmar con claridad
lo ilustrado en la figura anterior, esto es, el decai-
250 Qual2k
Anaite miento normal de la temperatura con la distancia
200
tendiendo a la temperatura natural del cuerpo de
150
agua. En cuanto a la temperatura en las paredes
100 izquierda y derecha, se observa que estas comien-
50 zan a incrementarse a una distancia aproximada
0 de 5,000 m del punto de descarga; se nota también
0 1000 2000 3000 4000
Distancia (m) que ambas temperaturas a la pared, tienen idéntico
Fig. 14. Resultados para el nitrógeno orgánico obtenidos con comportamiento y valor, lo cual se debe a la forma
ANAITE y QUAL2K regular del dominio.
Se resuelve ahora el mismo problema, conside-
rando además una descarga lateral ubicada a 5 km
El modelo bidimensional para el transporte de sobre la margen derecha y una descarga difusa de
escalares 3 km de longitud a partir del kilómetro 12 sobre la
Como se discutió anteriormente, la ecuación de margen izquierda. Para este ejemplo, la descarga
advección-difusión-reacción para el transporte de al centro del dominio tiene un valor de 35 °C, la
escalares que se incluye en el modelo es la ecuación descarga lateral es de 30 °C y la temperatura para
(9). La cual fue probada para el módulo de calidad la descarga difusa es de 25 °C constante; la tempe-
del agua MODCA1, que resuelve la temperatura ratura natural del agua es de 20 °C. La velocidad
 DINÁMICA Y CALIDAD DEL AGUA EN CAUCES NATURALES 125

T. en (C)
1500 27 2234
24 1251
22 2517
1200 21 6543
21 3114
1000 20 9993
20 6753
20 5716
750 20 5643
20 4941
20 4756
500 20 4502
20 4155
250 20 3735
20 3521

0
5000 10000 15000 20000
(a)
0

Temperatura máxima, en (C) Temperatura en las márgenes, en (C)

30 21.0
20.9
28 T MAX 20.8
T P. DER.
T P. IZQ. 20.7
26 20.6
20.5
24 20.4
20.3
22 20.2
20.1
20 20.0 (b)
0 5000 10000 15000 20000

Fig. 16. Plano de concentración de la temperatura. Descarga puntual y al centro del dominio.
(a) Temperatura total (b) Evolución al centro del dominio y en las márgenes izquierda
y derecha

media del flujo, U = 5.0 m/s, los coeficientes de que la temperatura se distribuye en casi todo lo largo y
dispersión longitudinal y lateral, Ex = 150 m2/s y ancho del dominio de estudio manteniendo valores ma-
Ey = 10 m2/s, respectivamente. yores a los de la temperatura natural. Para la descarga
Una vez resuelto el problema se concluye que el difusa, se puede notar en este ejemplo, cómo una des-
comportamiento presentado por este escenario de carga aunque relativamente pequeña en comparación
descarga es muy singular, resaltándose principalmente con la descarga puntual al centro, mantiene valores

T. en (C)
1500 34.1280
30.3679
28.1134
1200 24.6071
24.0107
1000 23.0890
21.9723
21.8426
750 21.6210
21.4983
21.3789
500 21.1109
20.9158
250 20.8190
20.4475

0 (a)
0 5000 10000 15000 20000

Temperatura máxima, en (C) Temperatura en las márgenes, en (C)

40 40
38 38
36 T MAX 36
T P. DER.
34 T P. IZQ. 34
32 32
30 30
28 28
26 26
24 24
22 22
20 20
(b)
0 5000 10000 15000 20000

Fig. 17. Plano de concentración de la temperatura. Descargas puntual, al centro del dominio,
lateral, y descarga difusa. (a) Temperatura total (b) Perfiles al centro y en ambas
márgenes
126 F. Torres-Bejarano y H. Ramírez-León

importantes y aún mayores que si se considerara como químicos, biológicos, tóxicos, etc. Actualmente se ha
una descarga puntual, siendo ésta el tipo de descarga ganado conocimiento experimental al respecto que
más influyente en el aumento de temperatura en gran aunado al desarrollo de equipo informático, permiten
parte del dominio de estudio (Fig. 17a). incluir aspectos como el mezclado lateral que ha sido
Los perfiles de temperatura para este escenario de un tema de auge en los últimos años y son particular-
descarga presentan un comportamiento característico, mente útiles en ríos anchos en donde la aproximación
principalmente para los perfiles a las márgenes. La unidimensional puede resultar cuestionable. La ten-
temperatura máxima tiende a su disminución gradual dencia mundial de los modelos hidrodinámicos y de
encontrándose con incrementos en su valor en los calidad del agua de los ríos es incluirlos en sistemas
puntos de la descarga puntual y difusa, para luego de información geográfica, en modelos de tiempo real
tender a su decaimiento normal tratando de llegar a y para evaluar indicadores ambientales. El modelo
la temperatura normal. La temperatura para la pared ANAITE continua desarrollándose en este sentido
derecha inicia con un incremento súbito en el punto (Islao et al. 2003, Laurel 2004).
de la descarga lateral; la temperatura a la pared iz-
quierda se mantiene constante a 20 °C hasta llegar al
punto de inicio de la descarga difusa (km 12), la cual REFERENCIAS
hace que se incremente el valor de la temperatura a
aproximadamente 25 °C conservándose constante Abbott M.B. (1979). Computational hydraulics: Elements
en todo el tramo de la descarga para luego decaer of the theory of free surface flows. Pitman, Londres.
gradualmente (Fig. 17b). 326 pp.
Berezowsky M. y Ramírez H. (1985). Calibración de la
rugosidad a lo largo del río Colorado, Instituto de In-
CONCLUSIONES geniería, UNAM, México. Informe técnico , Proyecto
No 4309, pp 150.
El modelo ANAITE demuestra una gran ver- Brown L. y Barnwell T. (1987). The enhanced stream wa-
satilidad para reproducir la hidrodinámica de una ter quality models QUAL2E and QUAL2E-UNCAS:
corriente unidimensional y bidimensional y puede ser Documentation and User Manual, Report EPA/600/3-
aplicado con gran facilidad a cualquier sistema con 87/007, U.S. EPA, Athens, GA, USA.
variación irregular o de sección constante. Asimismo, Calow P. and G. Petts (1994). The Rivers Handbook:
en su módulo de calidad del agua puede manejar Hydrological and Ecological Principles. Volume 2.
diferentes tipos de descargas, ya sean éstas difusas, Boston, MA. Blackwell Scientific Publications.
puntuales o transversales en todo el dominio como Chapra S.C., Pelletier G.J. y Tao H. (2006). QUAL2K: A
se vio anteriormente. Modeling Framework for Simulating River and Stream
Requiere relativamente de poca información Water Quality, Version 2.04: Documentation and Users
para su ejecución: algunas secciones transversales Manual. Civil and Environmental Engineering Dept.,
del tramo o longitud en estudio; para cada una de Tufts University, Medford, MA.
las secciones transversales necesita un número de Chow V.T. (1994). Hidráulica de canales abiertos. Mc-
puntos dispuestos en coordenada (x, y), requiere de Graw-Hill, Trad. Saldarriaga, J., Bogotá, 667 p.
una sección de control en donde se identifique el Elder J. (1959). The dispersion of marked fluid in turbulent
gasto que pase en un instante así como la elevación shear flow. J. Fluid Mech. 5, 544–560.
de la superficie libre. El programa estima coeficien- Ekebjaerg L. y Justesen P. (1991). An explicit scheme for
tes de rugosidad de Manning, velocidades medias y advection–diffusion modelling in two dimensions.
elevaciones de la superficie libre. Una vez generada Comp. Method. Appl. M. 88, 287-297.
la hidrodinámica se pueden aplicar los módulos de Fetter C.W. (1992). Contaminant Hydrogeology. Macmil-
calidad del agua, para lo que se requiere conocer al- lan, Nueva York, 458 p.
gunas características del sitio y el nivel de descarga. Fischer H.B. (1968). Methods for predicting dispersion
Tiene además la capacidad de incluir los términos de coefficients in natural streams, with applications to
creación/destrucción de los compuestos selecciona- lower reaches of the Green and Duwamish rivers.
dos (Torres-Bejarano 2003). Washington, U.S. Geological Survey Professional
Aunque los modelos de calidad del agua bidi- Paper, 582-A.
mensionales aún hoy en día son menos usados en Islao J.P., Ramírez León H., Navarro J.M. y Torres F.M.
ríos que los modelos unidimensionales, es más que (2003). Modelación numérica e indicadores de sen-
nada por el desconocimiento de los procesos físicos, sibilidad ambiental. Aplicación a la desembocadura
 DINÁMICA Y CALIDAD DEL AGUA EN CAUCES NATURALES 127

del río panuco. IX Congreso de Dinámica de Fluidos. Maskell J. (1991). Hydrodynamic and physical consider-
Cd. México. ations for water quality modeling. IWEM Water Qual-
Keulegan G.H. (1938). Laws of turbulent flow in open ity Modeling Symposium, Harrogate. The Institution of
channels. Journal of Research, Vol. 21, U.S. National Water and Environmental Management, Londres.
Bureau of Standard, artículo de investigación RP 1151, Rodríguez C. (2003). Modelación hidrodinámica en 2D
707-741. y 3D para lagos. Aplicación a la laguna de Coyuca y
Laurel A. (2004). El sistema ANAITE para el estudio de al embalse de Valle de Bravo. Tesis de Maestría en
la dinámica, contaminación y el diagnostico ambiental Ingeniería Hidráulica, DEPFI-UNAM.
en ríos. Tesis de Licenciatura, Escuela Superior de Shanahan P., Borchardt D., Henze M., Rauch W., Reichert
Ingeniería y Arquitectura del IPN, México. P., Somlyódy L. y Vanrolleghem P. (1998). River water
Liu H. et al. (FALTAN) (1977). Predicting dispersion quality No 1: I. Modeling Approach. Wat. Sci. Tech,
coefficient of streams. J. Environ. Eng. Div-ASCE, 38, 245-252.
103, EE1, pp. 59-69. Torres-Bejarano F. (2003). El programa ANAITE para el
McCutcheon S. (1989). Water quality modeling. Vol. I, estudio de la calidad del agua en sistemas de cauces
CRC Press, Boca Raton, FL. EUA, 334 pp. naturales. Aplicación al río Magdalena, Colombia.
Tesis de Maestría en Ingeniería, DEPFI-UNAM.