Modeling an aquifer: Numerical solution to the groundwater flow equation

Abstract

We present a model of groundwater dynamics under stationary flow and, governed by Darcy's law of water motion through porous media, we apply it to study a 2D aquifer with water table of constant slope comprised of a homogeneous and isotropic media; the more realistic case of an homogeneous anisotropic soil is also considered. Taking into account some geophysical parameters we develop a computational routine, in the Finite Difference Method, which solves the resulting elliptic partial equation, both in a homogeneous isotropic and in a homogeneous anisotropic media. After calibration of the numerical model, this routine is used to begin a study of the Ayamonte-Huelva aquifer in Spain, a modest analysis of the system is given, and we compute the average discharge vector as well as its root mean square as a first predictive approximation of the flux in this system, providing us a signal of the location of best exploitation; long term goal is to develop a complete computational tool for the analysis of groundwater dynamics.