Simulación del Flujo Subterráneo Regional del Sistema Acuífero Guaraní en Régimen Estacionario

Abstract

En este trabajo se presentan los resultados del modelo numérico del flujo subterráneo regional del Sistema Acuífero Guaraní (SAG), desarrollado en el marco del Proyecto para la Protección Ambiental y el Desarrollo Sostenible del Sistema Acuífero Guaraní financiado por el Banco Mundial. Se simuló el funcionamiento del denominado “SAG operativo”, acuífero que no incluye las formaciones geológicas que lo confinan en su parte superior e inferior, de menores permeabilidades. Además se lo consideró continuo en toda su extensión. El acuífero posee un funcionamiento hidráulico complejo, se comporta como medio transmisor de la recarga proveniente de las precipitaciones en regiones donde aflora, transfiriéndola a otras zonas de afloramiento donde descarga, sea a través de sus bordes, sea a través del contacto con escurrimientos superficiales. Otra fuente de salida la constituyen los pozos de bombeo que extraen en el orden de 1040 hm3 anuales. El flujo subterráneo del SAG fue simulado mediante el código numérico TRANSIN, un modelo que resuelve las ecuaciones del flujo de agua subterránea y el transporte de solutos el método de elementos finitos (EF), y que cuenta con un potente algoritmo de estimación automática de parámetros para la resolución del Problema Inverso. El área modelada cubrió 1.200.098 km2 abarcando toda la extensión del SAG operativo. El dominio computacional bidimensional se discretizó mediante 46862 elementos triangulares y 23890 nodos, en contraste con el modelo anterior del SAG que cubría sólo un 75 % del área, discretizada mediante 5423 elementos y 2846 nodos, es decir se realizó un refinamiento notable de la malla. El régimen temporal de las simulaciones fue estacionario. Se utilizó una combinación de condiciones de borde externas e internas de primer y segundo tipo y mixtas, además de representar la explotación subterránea y zonas de recarga/descarga. Durante el proceso de la zonificación/calibración de la conductividad hidráulica K, uno de los principales parámetros, se emplearon varias distribuciones areales con el fin de conseguir un ajuste aceptable con las medidas disponibles de carga hidráulica en puntos de observación. Los resultados de la modelación fueron razonables y coherentes con la información disponible, aunque evidenciaron una importante incertidumbre respecto de algunos procesos, por ejemplo la localización de las zonas de descarga del sistema regional a lo largo de tramos del río Uruguay y Paraná y, en general, la magnitud de los flujos de agua en el sistema. Las hipótesis de las posibles alternativas de descarga son variadas, aunque deben analizarse en mayor profundidad con metodologías alternativas para reducir las incertidumbres.