A Three Dimensional Finite Element Method And Multigrid Solver For A Darcy-Stokes System And Applications to Vuggy Porous Media