In the framework of the theory of approximate transformation groups proposed by Baikov, Gaziziv and Ibragimov [1], the first-order approximate symmetry operator is calculated for the Navier-Stokes equations. The symmetries of the coupled system obtained by expanding the dependent variables of the Navier-Stokes equations in the perturbation series with respect to a small parameter (viscosity) are used to derive approximate symmetries in the sense by Baikov et al.