A Convolution Boundary Element Method for Unsteady State Groundwater Flow in Homogeneous Aquifers


J. Kikani, R.N. Horne


Stanford Geothermal Workshop




Reservoir Modelling



File Size:


View File:


In this paper, Boundary Element (BEM) solutions were obtained for the transient flow of fluids through homogeneous, anisotropic porous media. The Greenís function method with Euler method of forward time differencing and Laplace transform method have been used by previous authors. Unlike these methods, this paper uses the fundamental solution to the differential equation and the convolution behavior of the resulting integrals to obtain an implicit and stable solution. This allows large time steps to be taken without significant loss in accuracy. Comparison with the Laplace transform method and Greenís function method with discrete time stepping, for two test cases, show that the method is very accurate. The computations however, become quite storage intensive owing to the dynamic increase in the number of stored matrices. It has been shown elsewhere that for certain problems with both Dirichlet and Neumann boundary conditions, asymptotic expression generated from exact solution is needed for starting the computational procedure. The present formulation alleviates this requirement. These solutions are developed for use in the analysis of pressure transients in complex reservoir problems.

ec2-18-204-56-97.compute-1.amazonaws.com, you have accessed 0 records today.

Press the Back button in your browser, or search again.

Copyright 1988, Stanford Geothermal Program: Readers who download papers from this site should honor the copyright of the original authors and may not copy or distribute the work further without the permission of the original publisher.

Attend the nwxt Stanford Geothermal Workshop, click here for details.

Accessed by: ec2-18-204-56-97.compute-1.amazonaws.com (
Accessed: Thursday 29th of September 2022 03:18:22 AM