Modelización del flujo en acuífero confinado sin discretización temporal

Resumen

Se presenta un programa interactivo para el modelado del flujo en acuífero confinado, homogéneo o heterogéneo (MGAC). La ecuación a partir de la cual se elabora el programa, es la obtenida al aplicar la transformada de Laplace a la ecuación diferencial de flujo en las dos direcciones del plano. El acuífero se discretiza espacialmente en celdas cuadradas y se plantea el sistema de ecuaciones transformadas en forma de diferencias finitas, no estando presente la variable tiempo. Para la solución del sistema se utiliza el método de SOR. La solución que se obtiene es invertida numéricamente para obtener los valores de los descensos, s_i, en vez de sus transformadas, s̄_i, y para recuperar la variable tiempo en vez del parámetro de inversión p. El método de inversión utilizado es el de Stehfest (1970). La aplicación del tiempo es de una vez y en su totalidad, por lo que los modelos no requieren discretización temporal. El uso del método de Stehfest conlleva un error para valores de Tt/r²S pequeños. La validez de este método se aproxima para acuífero confinado homogéneo mediante la comparación de los valores a los que se llega con la inversión de la expresión analítica transformada de la ecuación de Theis y los de la propia expresión de Theis. También se comparan los resultados obtenidos analíticamente con los de un modelo de acuífero genérico. Se reseñan los errores por truncación y las limitaciones originadas al no poder considerar con precisión los potenciales iniciales o fijos reales en régimen transitorio. Se muestran los tiempos de computación y la memoria requerida para distintos ejemplos.