Etude d’impact du Centre de Production Thermique de Porcheville (78)
description
Transcript of Etude d’impact du Centre de Production Thermique de Porcheville (78)
Etude d’impact du Centre deProduction Thermique
de Porcheville (78)
CEREA
Centre d’Enseignement et de Recherche Environnement Atmosphérique
Laboratoire Commun ENPC-EDF R&D
Objectifs
• évaluer l’incidence de divers scénarios de fonctionnement du CPT de Porcheville sur la qualité de l’air en région parisienne dans un rayon d’une centaine de kilomètres.
• Polluants : NOx, SO2, PM10 (PM2.5)
• Deux axes d’étude:– Vision générale : cartes de moyennes, …
– Vision réglementaire : autour des normes réglementaires (percentiles, seuils d’information, seuils d’alerte)
Présentation de Polyphemus
• Polyphemus est une plateforme de simulation de la qualité de l’air :
– Plusieurs modélisations : eulérienne, gaussienne
– Plusieurs paramétrisations pour les phénomènes physiques (dépôts, lessivages, diffusion turbulente, schémas chimiques, émissions, …)
– Mallet et al., Technical Note: The air quality modeling system Polyphemus, Atmos. Chem. Phys., 7, 5479-5487, 2007.
• Modélisation eulérienne offline :– Preprocessing : librairie C++ AtmoData pour les paramétrisations
des phénomènes physiques
– Simulation chimie-transport : Polair3D
– Postprocessing : librairie python AtmoPy pour la manipulation des résultats (visualisation, validation, …)
Résultats proposés
• Vision d’ensemble :– Cartes de moyennes de concentrations, de moyennes de
différences.
– Cartes de maximums de concentrations, de maximums de différences.
• Vision réglementaire :– En des points précis (stations d’observations et de
surveillance), impact des scénarios de fonctionnement de la centrale sur les valeurs suivies réglementairement pour chaque polluant (percentiles, nombre d’heures de dépassement de seuils).
• Éventuellement combiner les deux visions : cartes de percentiles ou cartes d’heures de dépassement.
Configurations de l’étude
• Les domaines d’étude :– « régional large » :
rectangle de 80 km x 90 km (mailles : 5-6 km)
– « régional fin » : rectangle de 16 km x 11 km (mailles : 1 km)
• Périodes de simulation : une année, à terme
cinq.
Champs météorologiquesConditions aux limites
• Modèle méso-échelle MM5 :– Maîtrise des résolutions spatiales (suffisamment fines : 1-2
km) et temporelles (1h).
– Champs consistants : les données d’entrée sont des analyses NCAR-NCEP + nudging à partir d’observations NCAR.
– Liberté sur les droits d’utilisation.
• Les concentrations aux limites des domaines sont issues de simulations à plus grande échelle (niveau européen) par une méthode de nesting one-way (interpolations aux bords).
Cadastres d’émissions de polluants
• Contrairement aux outils habituellement utilisés pour des études d’impact (gaussiens), Polair3D peut prendre en compte les sources d’émissions autres que le CPT et les intégrer au calcul de chimie-transport.
• Nécessite des inventaires d’émissions :– EMEP (résolution 50 km x 50 km).
– Airparif (1 km x 1 km), uniquement sur l’Ile-de-France.
– Dans un futur proche, 1 km x 1 km sur toute la France.
Scénarios d’émissions du CPT
• Deux stratégies possibles :– Fonctionnement continu. Avantage : on balaie toutes les
conditions météorologiques.– Fonctionnement épisodique. Avantage : on simule un
comportement plus proche de la réalité.
• Pour cette étude, 4 scénarios :– 450 h PCN, fioul TTBTS– 450 h PCN, fioul TTTBTS– 450 h PCN, fioul TTTBTS + SCR– 600 h PCN, fioul TTTBTS + 4 x SCR
Le modèle Polair3D
• Advection (transport par le vent)
• Diffusion (transport par la turbulence)
• Chimie (RACM : 72 espèces gazeuses
+ SIREAM : 17 espèces d’aérosols)
• Temps de calcul : – Environ 24 h de calcul pour le preprocessing
– Environ six jours de calcul par scénario pour un an simulé avec polair3d-siream
– Depuis peu, polair3d-siream est parallélisé, on peut donc envisager dans un futur proche des temps de calculs inférieurs à condition d’avoir les capacités de calculs adéquates.
Résultats (1)vision globale : NO2 en moyenne
Concentration moyenne
Scénario « sans CPT »
Différence de concentrations moyennes entre le scénario « 450h + TTBTS » et le scénario « sans CPT »
Différence de concentrations moyennes entre le scénario « 600h + fioul TTTBTS + 4 x SCR » et le scénario « sans CPT »
Résultats (2):vision globale : NO2 en max
Concentration moyenne
Scénario « sans CPT »
Différences maximales des concentrations horaires entre le scénario « 450h + TTBTS » et le scénario « sans CPT »
Différences maximales des concentrations horaires entre le scénario « 600h + fioul TTTBTS + 4 x SCR » et le scénario « sans CPT »
Résultats (3) :vision globale : SO2 en moyenne
Concentration moyenne
Scénario « sans CPT »
Différence de concentrations moyennes entre le scénario « 450h + TTBTS » et le scénario « sans CPT »
Différence de concentrations moyennes entre le scénario « 450h + fioul TTTBTS » et le scénario « sans CPT »
Résultats (4):vision globale : SO2 en max
Concentration moyenne
Scénario « sans CPT »
Différences maximales des concentrations horaires entre le scénario « 450h + TTBTS » et le scénario « sans CPT »
Différences maximales des concentrations horaires entre le scénario « 450h + fioul TTTBTS » et le scénario « sans CPT »
Résultats (5) :vision réglementaire : NO2
Station :
Porcheville
Scénario de
référence
« sans CPT »
450 h
fioul TTBTS
450 h
fioul TTTBTS
450 h
fioul TTTBTS
SCR
600 h
fioul TTTBTS
4 x SCR
Moyenne 32.09 32.12 32.12 32.12 32.11
P100 horaires 127.88 127.73 127.73 127.73 127.73
P99.8 horaires 100.58 100.58 100.58 100.58 100.58
Seuil d’information (h)
0 0 0 0 0
Seuil d’alerte (h) 0 0 0 0 0
Maximum des différences horaires
- 51.31 51.3 40.91 34.91
Maximum des différences journalières
- 5.45 5.45 4.21 3.56
Résultats (6) :vision réglementaire : NO2
Station :
Boissy-Mauvoisin
Scénario de référence
« sans CPT »
450 h
fioul TTBTS
450 h
fioul TTTBTS
450 h
fioul TTTBTS
SCR
600 h
fioul TTTBTS
4 x SCR
Moyenne 24.16 24.21 24.21 24.19 24.2
P100 horaires 166.57 169.52 169.36 168.12 167.9
P99.8 horaires 98.4 103.5 103.5 100.29 99.83
Seuil d’information (h)
0 0 0 0 0
Seuil d’alerte (h) 0 0 0 0 0
Maximum des différences horaires
- 60.8 60.75 47.92 39.46
Maximum des différences journalières
- 8.41 8.4 6.51 5.32
Conclusions
• En moyenne, l’impact se révèle faible quel que soit le polluant considéré, en accord avec les niveaux d’émissions.
• L’impact est également faible sur les valeurs réglementaires (percentiles, …) quel que soit le polluant et quelle que soit la station d’observation considérée.
• Épisodiquement et localement, l’impact peut éventuellement être important mais sans dépasser les seuils d’alerte en général.
Perspectives
• Sur-hauteurs : prendre en compte les vitesses et les températures d’éjection en sortie de cheminées pour évaluer une hauteur d’émissions (fin 2008)
• Interface graphique (fin 2008)
• Métaux lourds (fin 2008 – début 2009)
• Plume-in-grid : mettre en place et valider un modèle sous-maille, capable de modéliser les centrales thermiques comme sources ponctuelles avec un panache gaussien (courant 2009)
• Nesting : automatiser les liens entre les simulations sur des domaines imbriqués (courant 2009)