Dans un contexte environnemental fort, la législation Européenne actuelle impose aux centrales thermiques à gaz des limites d'émission à 25 ppmvd pour le monoxyde de carbone (CO) et les oxydes d'azote (NOx). L'exploitant est alors contraint de régler finement ses machines pour éviter l'imposition de pénalités financières. Afin d'orienter les opérateurs vers les réglages optimaux, un outil numérique est développé. Il a pour but de prédire les émissions polluantes avec précision (quelques ppmvd), rapidité (utilisation sur site) et pour toutes sortes de configurations (géométrie, conditions de fonctionnement). Toutefois, modéliser et simuler la formation d'espèces minoritaires dans un écoulement turbulent réactif reste, aujourd'hui encore, un challenge scientifique et numérique en raison des larges spectres spatiaux et temporels mis en jeu. La solution proposée dans cette thèse consiste à réaliser un calcul CFD (Computational Fluid Dynamics) basé sur un mécanisme réactionnel simple en vue de construire un réseau de réacteurs 0D incluant une cinétique chimique détaillée. Une méthodologie générale est présentée pour le découpage de l'écoulement en zones statistiquement homogènes qui, une fois connectées, forment un réseau de réacteurs. Deux modèles de réseaux de réacteurs sont définis et prennent ou non en compte les fluctuations turbulentes (densité de probabilité) et une distribution des temps de séjour dans chaque réacteur. Des études de sensibilité des émissions de NOx à l'hygrométrie, à la variabilité combustible et aux changements de charge sont menées. La comparaison aux mesures sur site souligne la représentativité des modèles proposés. Enfin, une nouvelle modélisation du terme source chimique moyen est introduite (méthode de tabulation du terme source chimique et nouvelle densité de probabilité) pour permettre une prédiction affinée des émissions de NOx dès le calcul CFD.