Numerical Solutions of Isothermal Pumping and Re-Injection in a Confined Aquifer

Abstract

Two numerical schemes are used for the computation of steady and transient pressure distribution in a porous medium. These two schemes are applied to two different types of boundary conditions. Numerical results for pressure distribution in a rectangular reservoir are compared with analytical solutions. It is found that the numerical results based on two different numerical methods are in close agreement with the analytical solutions for the region away from the sites of pumping and re-injection. Both the analytical solutions and the singular finite difference solutions are singular at the points of pumping and re-injection. The regular finite difference method, however, gives finite values of pressure at these points.