In this work we are interested in the numerical solution of the steady incompressible Navier-Stokes equations for fluid flow in pipes with varying curvatures and cross-sections. We intend to compute a reduced basis approximation of the solution, employing the geometry as a parameter in the reduced basis method. This has previously been done in a spectral element PN {-} PN {-}2 setting in two dimensions for the steady Stokes equations. To compute the necessary basis-functions in the reduced basis method, we propose to use a stabilized P1 {-} P1 finite element method for solving the Navier-Stokes equations on different geometries. By employing the same stabilization in the reduced basis approximation, we avoid having to enrich the velocity basis in order to satisfy the inf-sup condition. This reduces the complexity of the reduced basis method for the Navier-Stokes problem, while keeping its good approximation properties. We prove the well posedness of the reduced problem and present numerical results for selected parameter dependent three dimensional pipes.