La formation et l'évolution du disque épais de la Voie Lactée restent controversées. Nous avons utilisé un modèle de synthèse de la population de la Galaxie, le Modèle de la Galaxie de Besançon (Robin et al., 2003), qui peut être utilisé pour l'interprétation des données, étudier la structure galactique et tester différents scénarios de formation et évolution Galactique. Nous avons examiné ces questions en étudiant la forme et la distribution de métallicité du disque mince et du disque épais en utilisant l'approche de synthèse de la population. Nous avons imposé sur des simulations les erreurs d'observation et les biais afin de les rendre directement comparables aux observations. Nous avons corrigé les magnitudes et les couleurs des étoiles de la simulation, en utilisant un modèle d'extinction. Les modèles d'extinction disponibles ne reproduisent pas toujours la quantité exacte d'extinction le long de la ligne de visée. Un programme a été développé pour corriger la distribution de l'extinction en fonction de la distance le long de ces lignes. Les extinctions correctes ont ensuite été appliquées sur les simulations du modèle. Nous avons étudié la forme du disque mince en utilisant des données photométriques aux basses latitudes du sondage SDSS-SEGUE. Nous avons comparé qualitativement et quantitativement les observations et les simulations et nous avons essayé de contraindre la fonction de masse initiale. En utilisant la spectroscopie du relevé SEGUE, nous avons sélectionné les étoiles du turn-off de la séquence principale (MSTO) (Cheng et al 2012) et des géantes K pour étudier la distribution de métallicité du disque mince et du disque épais. Nous avons calculé une estimation de distance pour chaque étoile à partir de la relation entre les températures effectives et magnitudes absolues pour les catalogues observés et simulés. Ces deux catalogues ont les mêmes biais sur les distances, elles sont donc comparables. Nous avons développé un outil basé sur une méthode MCMC-ABC pour déterminer la distribution de la métallicité et étudier les corrélations entre les paramètres ajustés. Nous avons confirmé la présence d'un gradient de métallicité radiale de -0.079 ± 0.015 dex kpc−1 pour le disque mince. Nous avons obtenu une métallicité du disque épais au voisinage solaire de -0.47 ± 0.03 dex, compatible avec les résultats obtenus par les études précédentes. De plus, le disque épais ne montre pas de gradient, mais les données sont compatibles avec un gradient positif intérieur suivi d'un négatif extérieur. Nous avons ensuite appliqué les outils développés au relevé spectroscopique Gaia-ESO et calculé la distribution de métallicité des étoiles F/G/K dans le disque mince et épais en supposant une formation en deux époques du disque épais de la Voie Lactée. Nous avons obtenu une métallicité locale dans le disque épais de -0.23 ± 0.04 dex légèrement plus élevée que celle obtenue avec SEGUE mais en accord avec Adibekyan et al. (2013) et un gradient de métallicité radiale du disque épais en accord avec notre analyse précédente des données de SEGUE et la littérature. La métallicité locale est en accord avec la littérature au niveau de 3σ mais parce que les données GES sont préliminaires, une analyse plus approfondie avec plus de données et de meilleurs calibrations doit être faite. L'existence d'un gradient plat dans le disque épais peut être une conséquence d'une formation à partir d’un gaz turbulent et bien homogène, ou bien un fort mélange radial a brassé après coup les étoiles. / The formation and evolution of the thick disc of the Milky Way remain controversial. We made use of a population synthesis model of the Galaxy, the Besançon Galaxy Model (Robin et al. 2003), which can be used for data interpretation, study the Galactic structure and test different scenarios of Galaxy formation and evolution. We examined these questions by studying the shape and the metallicity distribution of the thin and thick disc using the population synthesis approach. We imposed on simulations observational errors and biases to make them directly comparable to observations. We corrected magnitudes and colors of stars, from the simulation, using an extinction model. The available extinction models do not always reproduce the exact quantity of extinction along the line of sight. A code to correct the distribution of extinction in distance along these lines have been developed and the corrected extinctions have been applied on model simulations. We studied the shape of the thin disc using photometric data at low latitudes from the SDSS-SEGUE survey. We compared qualitatively and quantitatively observations and simulations and try to constrain the Initial Mass Function. Using the spectroscopic survey SEGUE we selected Main Sequence Turnoff (MSTO) stars (Cheng et al 2012) and K giants to study the metallicity distribution of the thin and thick discs. We computed a distance for each star from the relation between effective temperatures and absolute magnitudes for the observed and simulated catalogs. These two catalogues have the same biases in distances, therefore are comparable. We developed a tool based on a MCMC-ABC method to determine the metallicity distribution and study the correlations between the fitted parameters. We confirmed a radial metallicity gradient of -0.079 ± 0.015 dex kpc−1 for the thin disc. We obtained a solar neighborhood metallicity of the thick disc of -0.47 ± 0.03 dex similar to previous studies and the thick disc shows no gradient but the data are compatible with an inner positive gradient followed by a outer negative one. Furthermore, we have applied the developed tools to the Gaia-ESO spectroscopic survey and computed the metallicity distribution of F/G/K stars in the thin and thick disc assuming a two epoch formation for the thick disc of the Milky Way. We obtained a local metallicity in the thick disc of -0.23 ± 0.04 dex slightly higher than the one obtained with SEGUE but in agreement with Adibekyan et al. (2013) and a radial metallicity gradient for the thick disc in agreement with our previous analysis of SEGUE data and the literature. The local metallicity is in fair agreement with literature at the 3σ level but because the GES data is an internal release under testing further analysis with more data and better calibrations have to be done. The existence of a flat gradient in the thick disc can be a consequence of an early formation from a highly turbulent homogeneous well mixed gas, unless it has suffered heavy radial mixing later on.