Coeficiente de reaireacion en rios: como elegir entre O'Connor-Dobbins, Churchill, Owens-Gibbs y Tsivoglou
El coeficiente de reaireacion en rios (Ka, en unidades de 1/dia) es el parametro mas sensible en cualquier modelo de oxigeno disuelto, segun los analisis de sensibilidad globales tipo Sobol y Morris reportados por la EPA y por el equipo que mantiene QUAL2K. Si Ka cambia un 30%, el OD modelado en el punto critico aguas abajo de un vertimiento puede moverse mas de 2 mg/L. Existen al menos nueve formulas empiricas publicadas y ninguna funciona en todo el rango hidraulico. Elegir mal es la causa numero uno de modelos de calidad de agua que pasan la calibracion en escritorio pero fallan al ser auditados.
¿Listo para modelar tu río?
Prueba el motor de predicción QUAL2K con tus propios datos
Iniciar Motor de PredicciónQue es el coeficiente de reaireacion y por que importa tanto
La reaireacion atmosferica es el flujo neto de oxigeno desde la atmosfera hacia el agua. En un rio, ese flujo es proporcional al deficit de oxigeno disuelto y al coeficiente Ka, expresado como tasa de primer orden. Cuanto mas turbulento y delgado sea el flujo, mayor es la transferencia de oxigeno por unidad de volumen, y mayor sera Ka. En rios profundos y lentos, Ka cae a valores cercanos a 1 1/d; en rios someros y rapidos puede superar los 50 1/d.
La importancia practica es directa: en un balance de masa estacionario tipo Streeter-Phelps moderno, el deficit critico de OD aguas abajo de una descarga depende del cociente entre la tasa de deoxigenacion (Kd) y Ka. Si subestimas Ka, predices un deficit mayor del real y dimensionas en exceso un tratamiento; si lo sobrestimas, declaras cumplimiento donde el rio sigue estresado. Para un ingeniero ambiental que defiende un disenio de PTAR ante una autoridad ambiental con datos limitados, la eleccion de la formula de Ka es muchas veces la diferencia entre cumplir o exceder el OD minimo regulatorio.
La formula clasica de O'Connor-Dobbins (1958): cuando aplicarla
La formula de O'Connor y Dobbins (1958, ASCE) es la mas conocida y la que aparece por defecto en muchos textos. Su forma es:
Ka (1/d) = 3.93 · V^0.5 / H^1.5
donde V es la velocidad media (m/s) y H la profundidad media (m). El rango de validez derivado de los datos originales es H entre 0.30 y 9.10 m, con velocidades entre 0.15 y 0.49 m/s. Funciona bien para rios de llanura, profundos, con velocidad moderada. Su debilidad es sistematica: subestima Ka en rios rapidos, especialmente cuando el numero de Froude se acerca a 1, porque la teoria de penetracion no captura la turbulencia adicional inducida por la rugosidad del lecho.
Churchill, Owens-Gibbs y Tsivoglou-Neal: el resto del menu
La formula de Churchill (1962) fue derivada en rios del Tennessee Valley Authority aguas abajo de presas, donde el agua sale fria, sin carga de DBO y con velocidades altas. Su expresion:
Ka (1/d) = 5.026 · V^0.969 / H^1.673
Es la opcion preferida cuando la velocidad media supera los 0.55 m/s y la profundidad esta entre 0.6 y 3.4 m. Para rios con esas caracteristicas, Churchill suele dar valores de Ka entre dos y cuatro veces mas altos que O'Connor-Dobbins.
Owens-Gibbs (1964) fue desarrollada para arroyos pequenios y poco profundos del Reino Unido. Es la unica formula tradicional con suficiente respaldo empirico cuando H baja de 0.6 m, terreno donde O'Connor-Dobbins se vuelve numericamente inestable porque H^1.5 tiende a cero. Tsivoglou-Neal (1976) rompe el esquema: en lugar de velocidad y profundidad, usa la pendiente del lecho (S) y la velocidad. Su forma simplificada es Ka aproximadamente igual a c · S · V, donde c depende del rango de caudal. Es la opcion natural para rios de montania donde la pendiente domina la disipacion de energia y la entrada de aire por estructuras y rapidos.
Como elegir la formula correcta segun profundidad y velocidad
La USGS Professional Paper 737 y revisiones posteriores como Wilson y Macleod (1974) o Pelletier y Chapra (2008) coinciden en algo importante: ninguna formula domina estadisticamente en todo el rango. Suele existir un grupo de tres o cuatro formulas con desempenio equivalente para una hidraulica dada. La estrategia razonable es elegir el grupo correcto y luego promediar o usar el valor central.
| Hidraulica del tramo | Formula recomendada | Rango tipico de Ka (1/d a 20 grados C) |
|---|---|---|
| Rio profundo (H mayor a 3 m), V baja (menor a 0.5 m/s) | O'Connor-Dobbins | 0.5 a 4 |
| Rio intermedio (H 0.6 a 3 m), V alta (mayor a 0.55 m/s) | Churchill | 3 a 15 |
| Arroyo somero (H menor a 0.6 m) | Owens-Gibbs | 10 a 50 |
| Rio de montania, alta pendiente (S mayor a 0.005) | Tsivoglou-Neal | 10 a 60 |
Lo que casi nadie te cuenta: por que tu Ka calculada se desvia del observado
Aca aparece la limitacion practica que las formulas tradicionales no contemplan: las estructuras hidraulicas. Vertederos, saltos, presas pequenias y caidas naturales producen una reaireacion adicional que ninguna formula derivada de canales prismaticos captura. Si tu tramo de estudio cruza un vertedero o una pequenia presa, el Ka medio del tramo puede ser dos o tres veces mayor que el predicho por O'Connor-Dobbins solo por efecto local de esa estructura.
Pro tip de campo: usa la correccion de Butts y Evans (1983) para saltos. Para una caida de altura h (metros) a temperatura T (grados Celsius), el deficit residual de OD aguas abajo se relaciona con el de aguas arriba como Da/Db igual a 1 dividido entre (1 mas 0.38 · a · b · h · (1 menos 0.11 · h) · (1 mas 0.046 · T)), donde a es el factor de calidad del agua (1.0 agua limpia, 0.65 efluente de PTAR) y b es el factor del tipo de estructura (1.0 vertedero plano, 1.3 escalonado). Una caida de 1 m en agua a 20 grados C en un rio con efluente tipico (a igual a 0.8, b igual a 1.0) elimina aproximadamente el 23% del deficit de OD existente. La fuente original es el USGS Open-File Report 1981/0648. En ningun modelo numerico aparece esta correccion automaticamente: hay que aplicarla como condicion de contorno entre tramos.
El otro factor que desvia los valores predichos es la temperatura. Ka tiene su propia correccion theta, tipicamente theta igual a 1.024, lo que implica que un cambio de 10 grados C en el rio modifica Ka un 27%. En modelos donde la temperatura varia mucho a lo largo del dia, asumir un Ka constante puede generar errores de calibracion del orden del 0.5 a 1 mg/L en el ciclo diel.
Implementacion en QUAL2K: la opcion Internal y como validarla
QUAL2K, segun el manual oficial de Chapra, Pelletier y Tao, ofrece una opcion automatica llamada Internal que selecciona dinamicamente entre O'Connor-Dobbins, Churchill y Owens-Gibbs en funcion de la velocidad y la profundidad del tramo. Es la recomendacion para usuarios que no quieren forzar manualmente la formula tramo por tramo. La logica usa el diagrama de Covar (1976), que divide el plano V-H en zonas y asigna la formula del trio clasico que mejor se ajusta a esa region.
Para validar lo que la opcion Internal te dio, hay tres caminos. Primero, contrastar con un calculo manual usando la formula que aplica a tu hidraulica dominante. Segundo, comparar con el Ka derivado de un balance de masa cerrado entre dos puntos de monitoreo donde se midio OD y los flujos atmosfericos no son la unica explicacion del cambio observado. Tercero, ejecutar el simulador de Hydrolitica con dos opciones de Ka (Internal y formula forzada) y observar la sensibilidad del OD modelado en el punto critico: si la diferencia es menor a 0.3 mg/L, el modelo no es sensible a la eleccion en ese tramo y puedes seguir adelante. Si supera 1 mg/L, necesitas datos de campo para decidir.
Una alternativa pragmatica: para sitios sin datos de monitoreo, la calculadora online de Ponce (San Diego State University) permite estimar Ka con varias formulas en simultaneo y observar el rango. Si las formulas aplicables a tu hidraulica caen en un rango estrecho, ese es tu intervalo de confianza implicito; si el rango supera un orden de magnitud, hay que revisar los datos de entrada.
Como propagar incertidumbre cuando no estas seguro de la formula
Cuando el rango de Ka entre las formulas aplicables supera un factor de dos, la mejor practica no es elegir un valor unico sino propagar la incertidumbre. Una metodologia simple consiste en correr el modelo tres veces: con el Ka minimo del grupo de formulas aplicables, con el Ka medio y con el Ka maximo. La banda de OD resultante en el punto critico aguas abajo es la incertidumbre que debes reportar a la autoridad ambiental.
En la practica, esto produce graficos donde el OD minimo modelado se presenta como una linea central con una sombra de variabilidad. Si la linea central cumple el limite regulatorio pero el limite inferior de la banda no, la conclusion tecnica es que el sistema esta en el limite de cumplimiento y se recomienda monitoreo continuo o medidas de respaldo. Si toda la banda cumple, hay margen suficiente para validar el disenio. Si toda la banda incumple, hay un problema estructural en el sistema que ningun ajuste de Ka va a resolver.
Tres advertencias adicionales sobre Ka que rara vez se documentan en los manuales basicos. Primera, en rios con vegetacion sumergida densa (macrofitas), la velocidad efectiva en la columna de agua puede ser 30 a 60% menor que la velocidad media calculada a partir del caudal y el area mojada, y eso reduce Ka mas alla de lo que cualquier formula clasica anticipa. Segunda, en tramos con presencia de hielo o de espuma persistente (descargas industriales), la transferencia atmosferica se inhibe casi por completo y Ka tiende a cero localmente. Tercera, los rios de planicie con elevada turbidez no necesariamente tienen menor Ka, pero si pueden tener saturacion de OD reducida por la respiracion bacteriana en suspension, lo que hay que separar conceptualmente del termino de reaireacion.
El cierre numerico: cuanto mueve realmente Ka tu deficit critico
Volvamos al numero del lead: una variacion del 30% en Ka mueve el deficit critico de OD mas de 2 mg/L en escenarios tipicos. Eso significa que para un rio con OD de saturacion de 8.5 mg/L y un OD minimo regulatorio comun de 4 o 5 mg/L, la diferencia entre cumplimiento y exceso depende muchas veces de un parametro que se calcula con una formula de seis decimas de exactitud. La conclusion practica es doble. Primero, no aceptes el Ka por defecto sin verificar el rango hidraulico de tu rio. Segundo, declara siempre en tu informe tecnico que formula usaste, por que la elegiste y cual es el rango plausible de Ka segun el grupo de formulas equivalentes. Esa transparencia es lo que separa un modelo defendible de uno que sera devuelto en revision. Y si tienes el luxo de medir, medir es siempre mejor que estimar: una sola medicion in situ con tecnicas de gas trazador (radon, propano) puede ahorrarte semanas de discusion sobre cual formula aplica.