Turbulent premixed methane-air combustion in an axisymmetric cylindrical chamber is numerically simulated by a second-order finite volume method. Combustion chemistry is modeled by a reduced mechanism consists of 16 species and 31 elementary reactions, which includes some C1 and C2-chain species. Reaction rates are considered as minimum rates between Arrhenius rates and eddy break-up rates. Turbulent modeling is done via standard k-epsilon model with modified coefficients in flame cases and empirical inlet boundary conditions. All the numerical results were compared with experimental data for 3 turbulent premixed methane- air flames. Since combustion is occurred in distributed reaction zone, the turbulent intensities play an important role in prediction of the flame characteristics such as temperature and species mass fractions, specially when the flame occurs near the broken reaction zone. The results show the important role of C2-chain species in prediction of turbulent premixed flames behavior, especially in CO and O-2 mass fractions and the temperature field.