La zona de mezcla de un vertimiento es el tramo del río, aguas abajo de la descarga, donde el efluente todavía no se ha mezclado con el agua del cuerpo receptor y donde la norma tolera que se superen los criterios de calidad. En Perú, el artículo 7 del D.S. 004-2017-MINAM dice que la Autoridad Nacional del Agua (ANA) verifica los Estándares de Calidad Ambiental para Agua fuera de esa zona. En Colombia, el Decreto 1076 de 2015 la define como el área "indispensable para que se produzca mezcla homogénea" del vertimiento con el cuerpo receptor, donde se permite sobrepasar los criterios de calidad siempre que se cumplan las normas de vertimiento. La diferencia práctica está en el cómo: Perú tiene desde 2017 una guía de la ANA con ecuaciones y restricciones explícitas; Colombia remite a la Guía Nacional de Modelación, que da el marco pero deja la elección del método al modelador.
Importa tu libro QUAL2Kw o abre el proyecto de ejemplo. Los resultados reproducen 15 modelos QUAL2Kw calibrados dentro del 5% o de una pequeña banda absoluta.
Ver el motor de predicciónLa Guía para la determinación de la zona de mezcla y la evaluación del impacto de un vertimiento de aguas residuales tratadas a un cuerpo natural de agua fue aprobada por la Resolución Jefatural 108-2017-ANA como herramienta "referencial", en tanto se apruebe una metodología vinculante. La propia guía aclara que no excluye otras metodologías debidamente sustentadas, pero recomienda aplicarla porque agiliza la evaluación de la ANA. Cubre ríos y quebradas, cuerpos lénticos y el medio marino-costero; lo que sigue se refiere a ríos.
La guía propone el método simplificado de la EPA de EE. UU. para descargas superficiales. Para un vertimiento en la orilla:
L_ZdM = W² · u / (2 · π · Dy)
y para un vertimiento en el centro del cauce el denominador es 8 · π · Dy, es decir, una longitud cuatro veces menor. W es el ancho medio del río en los 500 m aguas abajo del vertimiento, u la velocidad media en el punto de descarga y Dy el coeficiente de dispersión lateral:
Dy = c · d · u*, con u* = √(g · d · s)
donde d es la profundidad media aguas abajo, s la pendiente del cauce medida entre la descarga y un punto unos 500 m aguas abajo, y c un factor de irregularidad del cauce: 0,1 para ríos rectos de sección rectangular, 0,3 para ríos canalizados, 0,6 para cauces naturales con meandros moderados, 1,0 para meandros marcados y más de 1,0 con cambios bruscos de dirección de 90° o más. Las condiciones hidráulicas deben ser las del período crítico (estiaje); si el caudal del vertimiento cambia mucho durante el año, se calcula para estiaje y avenida y se usa el mayor valor. Como alternativa, la guía acepta modelos numéricos como EFDC, CORMIX o Visual Plumes, simulando las condiciones más desfavorables.
La longitud calculada no se acepta sin más. La guía la limita en cuatro casos:
La guía menciona el criterio usado en muchos estados de EE. UU. (el caudal mínimo medio de 30 días con período de retorno de 5 años, sobre series de diez años o más) y, como pocas veces hay series tan largas, propone dos métodos simplificados. Con una serie de al menos 5 años, el caudal crítico es el menor de los caudales medios mensuales de los últimos 5 años, sin considerar extremos como el fenómeno El Niño. Sin serie, se mide el caudal al final del mes de menor precipitación media multianual y se divide entre 2, un factor de seguridad por la variabilidad entre años y por depender de una sola medición. Si la zona de mezcla quedó restringida, el caudal disponible también se reduce, porque no toda la sección alcanza a diluir:
Q_crít,red = Q_crít · √(2 · π · Dy · L_ZdM / u) / W (descarga en la orilla; 8 · π para el centro)
Para la mayoría de parámetros la concentración más alta aparece justo al final de la zona de mezcla, y la guía evalúa el cumplimiento con un balance de masas: C_a = (C_RH · Q_crít + C_vert · Q_vert) / (Q_crít + Q_vert), que debe ser menor o igual al ECA-Agua. De ahí despeja la carga, el caudal o la concentración máxima admisibles del efluente. Dos parámetros son la excepción porque su peor valor aparece lejos de la zona de mezcla: los nitratos, que crecen con la nitrificación del amonio, y el oxígeno disuelto, cuyo mínimo suele estar varios kilómetros aguas abajo. Para los nitratos (en cuerpos de categoría 1, o de categoría 3 que tributan a un cuerpo léntico) la guía usa un balance de nitrógeno que supone nitrificación completa y aplica el caudal crítico total, no el reducido. Para el oxígeno propone un método simplificado del déficit máximo y, como alternativa, un modelo de autodepuración; en una nota al pie menciona el modelo de Streeter-Phelps aplicado en programas como RIOS-4, RIOS EP y QUAL2K.
Verifica el marco vigente: la guía se apoya en el artículo 5 del D.S. 023-2009-MINAM, que fue derogado menos de un mes después por el D.S. 004-2017-MINAM (publicado el 7 de junio de 2017). El nuevo decreto mantiene el encargo a la ANA de fijar la metodología de la zona de mezcla (artículo 7.3), y la guía sigue publicada como referencia. Antes de radicar, confirma con la ANA si ya existe una metodología vinculante que la reemplace.
En Colombia no hay una ecuación oficial. El artículo 2.2.3.3.3.5 del Decreto 1076 dice que la autoridad ambiental controla los criterios de calidad por fuera de la zona de mezcla, "la cual será determinada para cada situación específica por dicha autoridad", teniendo en cuenta la Guía Nacional de Modelación. El numeral 9 del artículo 2.2.3.3.5.3 exige en la Evaluación Ambiental del Vertimiento los estudios y diseños de la estructura de descarga para minimizar la extensión de la zona de mezcla.
La Guía Nacional de Modelación (sección 5.3) fija cuatro reglas:
La longitud de la zona de mezcla tiene además una consecuencia directa en el monitoreo de la EAV: la Guía pide un punto aguas abajo "a una distancia mayor a la longitud estimada para la zona de mezcla". Lo explicamos en la guía de modelación para la Evaluación Ambiental del Vertimiento.
Un río ilustrativo en estiaje: ancho medio 20 m, velocidad 0,4 m/s, profundidad 0,8 m y pendiente 0,002. La velocidad de corte es u* = √(9,81 · 0,8 · 0,002) ≈ 0,125 m/s.
| Método | Coeficiente lateral | Descarga en la orilla | Descarga en el centro |
|---|---|---|---|
| Guía ANA, meandros moderados (c = 0,6) | Dy ≈ 0,060 m²/s | ≈ 424 m | ≈ 106 m |
| Guía ANA, meandros marcados (c = 1,0) | Dy ≈ 0,100 m²/s | ≈ 254 m | ≈ 64 m |
| USGS (Kilpatrick y Wilson, 1989), 95% de mezcla | Ez ≈ 0,020 m²/s | ≈ 3.200 m (K = 0,4) | ≈ 800 m (K = 0,1) |
La tercera fila usa la ecuación del manual del USGS que cita la Guía colombiana: L = K · u · W² / Ez, con K = 0,1 para una inyección central y 0,4 para una lateral, que corresponde al 95% de mezcla, y el coeficiente lateral de la nota de su tabla 1 (Ez = 1,13 d^1,5 s^0,5, con d en pies y Ez en pies²/s; aquí convertido a metros). El mismo río da entre 64 m y más de 3 km según la fórmula y el coeficiente. Ninguna está "mal": cambian la definición de mezcla completa y el coeficiente de dispersión lateral. Por eso el informe debe decir qué ecuación usaste, con qué coeficiente y por qué, y, si el resultado es decisivo, contrastarlo con un ensayo con trazador o un modelo de campo cercano. En Perú, además, cualquier resultado por encima de 500 m choca con el tope de la guía, y la salida sería rediseñar la descarga, no discutir la fórmula.
Sigamos con el caso peruano. Supón una captación 200 m aguas abajo de una descarga en la orilla: la zona de mezcla queda limitada a 100 m y, con c = 0,6, el caudal crítico disponible baja a cerca de la mitad (factor √(2π · 0,060 · 100 / 0,4) / 20 ≈ 0,49). Con un caudal crítico de 1,5 m³/s, un vertimiento de 0,05 m³/s con 60 mg/L de un parámetro y 2 mg/L en el río, el balance de masas da unos 3,9 mg/L con el caudal completo y unos 5,7 mg/L con el caudal reducido. Esa diferencia puede decidir si el vertimiento cumple o no.
La zona de mezcla es un fenómeno de campo cercano: dilución física de una sustancia que se trata como conservativa. La longitud de influencia es otra pregunta: hasta dónde, aguas abajo, el vertimiento sigue cambiando la calidad del río de forma apreciable. Para la materia orgánica y el nitrógeno la respuesta está en el campo lejano, porque el mínimo de oxígeno disuelto y el máximo de nitratos ocurren lejos de la descarga, como reconoce la propia guía de la ANA.
El término "longitud de influencia" no está definido en el Decreto 1076 ni en la Guía Nacional de Modelación. La ANLA, por ejemplo, ha publicado un documento metodológico sobre la longitud de influencia de vertimientos, pero la norma nacional no fija qué variable mirar ni en qué punto se considera que el efecto terminó. Si tu evaluación lo pide, pregunta a la autoridad qué variable y qué criterio de corte usar, y escríbelo en el informe.
En la práctica ambas cosas se calculan con herramientas distintas y se complementan. Los modelos unidimensionales como QUAL2K y QUAL2Kw suponen el cauce mezclado en la vertical y en la transversal, así que describen bien el río aguas abajo de la zona de mezcla, pero no dentro de ella. La secuencia razonable es estimar la zona de mezcla con una ecuación o un modelo de campo cercano, ubicar ahí el punto de control y usar el modelo 1D calibrado para el perfil de oxígeno, amonio y nitratos, del que sale la longitud de influencia.
Para el balance de masas al final de la zona de mezcla sirve la calculadora de balance de masas, y para una primera curva de oxígeno, la calculadora Streeter-Phelps. Para el perfil completo aguas abajo, con varios vertimientos, nitrificación, algas y temperatura, el motor de Hydrolitica corre desde el navegador un port del solver de QUAL2Kw 5.1 e importa libros QUAL2Kw existentes. Ninguna de estas herramientas calcula por sí sola la zona de mezcla: esa estimación sigue siendo una decisión del modelador, con la ecuación o el modelo de campo cercano que la autoridad acepte.
Pon en práctica lo que acabas de leer con el motor de predicción QUAL2K completo
Prueba gratis el motor de predicción