Boundary integral equations for the computational modelling of three-dimensional steady grounwater flow problems