The vertical averages of the incompressible Navier-Stokes equations are studied from the point of view of numerical analysis: existence of solution and convergence of algorithms. Three formulations are analysed; existence theorems are obtained when the Reynolds number is small. Convergence of a time implicit algorithm is shown, while the discretization in space is achieved with the finite element method.