Abstract: A boundary element method(BEM) for calculating groundwater flow to a multiple well system in a leaky aquifer is presented.This method treats the well walls as internal boundaries and avoids discretizing the boundaries of the well walls.At first, the boundary integral equations of the water head H and its normal derivative ∂/∂n are derived.After discretization, a system of linear algebraic.equations forsolving the integral equations is. obtained To illustrate application of the method,some type well problems are solved and the results show that the method is accurate.