FR 3049653 A1 20171006 FR 1652924 A 20160404 1 La présente invention concerne le domaine de la modélisation des émissions polluantes d'un véhicule, en particulier d'un véhicule automobile. L'estimation des émissions polluantes des véhicules en usage réel est un enjeu majeur pour faire face aux challenges de santé publique et environnementaux rencontrés dans les grandes métropoles. Pour y parvenir, l'approche principalement utilisée actuellement consiste à instrumenter un véhicule à l'aide d'un Système de Mesure des Emissions Portable (PEMS, de l'anglais « Portable emissions measurement system »). Cependant, le coût de ces systèmes reste prohibitif et ne permet pas leur diffusion à grande échelle. En conséquence, ils ne permettent pas de mesurer les émissions polluantes en usage réel à grande échelle. Pour pallier ce problème lié à cette instrumentation, il existe différentes approches de modélisations. L'état de l'art des modèles de polluants est constitué de deux grandes familles, les modèles macroscopiques et les modèles microscopiques. Le domaine de l'invention se place ici à l'échelle du véhicule, c'est-à-dire à l'échelle macroscopique. Les modèles haute fréquence (échelle microscopique) modélisant les phénomènes physiques en jeu dans la combustion de manière extrêmement détaillée, comme par exemple décrit dans la demande de brevet FR2984557 (US 2013/0158967), ne sont pas adaptés pour estimer les émissions polluantes en usage réel et à grande échelle, à cause de leur temps de calcul prohibitif, ainsi que par le grand nombre de paramètres requis. Du côté des approches macroscopiques, le modèle le plus répandu utilise des facteurs d'émissions, c'est-à-dire des coefficients moyens d'émissions par kilomètre parcouru. Ces coefficients dépendent du polluant considéré, de la vitesse moyenne sur le trajet, du type de motorisation (essence/Diesel) et de la norme euro du véhicule. On peut citer notamment la demande de brevet CN102054222A qui couple des acquisitions GPS à des facteurs d'émissions. Cependant, ces approches ne considèrent qu'un comportement moyen du véhicule et du conducteur, et ne permettent donc pas de modéliser leur usage réel. En particulier, l'impact du style de conduite sur les émissions y est négligé. Ils permettent d'estimer les émissions moyennes totales sur un roulage complet d'une durée suffisamment longue. En outre, de nombreuses études se sont appuyées sur les facteurs d'émissions pour estimer les émissions polluantes. On peut citer notamment le document Liu, H., Chen, X., Wang, Y., & Han, S. (2013). Vehicle Emission and Near-Road Air Quality Modeling for Shanghai, China: Based on Global Positioning System Data from Taxis and Revised MOVES Emission Inventory.Transportation Research Record: Journal of the Transportation Research 3049653 2 Board, (2340), 38-48, dans lequel des données GPS connectées sur des taxis sont utilisées en entrée des facteurs d'émission, et le document Sentoff, K. M., Aultman-Hall, L., & Holmén, B. A. (2015). Implications of driving style and road grade for accurate vehicle activity data and emissions estimates. Transportation Research Part D: Transport and Environment, 5 35, 175-188, qui étudie la précision de l'estimation par facteur d'émission. La conclusion de cette dernière étude montre que la principale source d'erreur provient de la non prise en compte de l'impact du style de conduite et de la pente. Or, pour capter ces phénomènes, il est nécessaire d'utiliser un niveau de modèle plus fin, appelé modèle de polluant microscopique, dont l'entrée est un signal de vitesse véhicule de 1 Hertz. 10 De plus, il existe déjà de nombreux modèles microscopiques, prévus pour être utilisés dans le cadre d'étude d'ingénierie, inadaptés à une implémentation embarquée, par exemple sur téléphone intelligent (« smartphone ») ou sur un serveur. Les modèles microscopiques les plus connus sont le modèle CMEM (de l'anglais « comprehensive modal emissions » 15 pouvant être traduit par modèle modal approfondi des émissions) et le modèle PHEM (de l'anglais « Passenger Car and Heavy Duty Emission Model » pouvant être traduit par modèle d'émissions de véhicule de tourisme et de véhicule lourd) qui prennent comme données d'entrée des profils de vitesse instantanée. Le modèle CMEM, développé dans les années 1990 par University of California (Riverside), a représenté une contribution significative lors 20 de sa sortie, mais présente aujourd'hui deux inconvénients majeurs : des résultats imprécis dus à des simplifications excessives, et une mauvaise prise en compte des motorisations Diesel. Quant à lui, le modèle PHEM présente la limitation de ne proposer qu'un modèle type pour de nombreux véhicules relativement différents. Plusieurs documents ont déjà étudié le 25 couplage de ces modèles microscopiques avec des données GPS, fournies par des boitiers connectés installés dans le véhicule. On peut citer notamment le document Pluvinet, P., Gonzalez-Feliu, J., & Ambrosini, C. (2012). GPS data analysis for understanding urban goods movement. Procedia-Social and Behavioral Sciences, 39, 450-462 qui détaille une étude où des données GPS (pour « global positioning system » signifiant géo- 30 positionnement par satellite) réelles ont été utilisées pour alimenter un modèle d'émission microscopique. Cependant cette étude présente une limitation majeure puisque le type de véhicule et le type de moteur ne sont pas pris en compte. La prise en compte du type de véhicule et de motorisation est une problématique majeure pour la modélisation des émissions polluantes en usage réel à grande échelle. En effet, s'il est possible de récupérer 35 les mesures GPS, il est autrement plus compliqué de connaître les spécifications techniques de son moteur et de paramétrer le modèle en conséquence. Toutefois, ces spécifications sont indispensables pour paramétrer les modèles microscopiques précédemment cités. Il 3049653 3 s'agit d'une limitation majeure des modèles microscopiques existants qui ne peuvent pas être utilisés pour un déploiement à grande échelle puisqu'ils nécessitent un travail fastidieux de paramétrage manuel pour chaque véhicule considéré. 5 Pour pallier ces inconvénients, la présente invention concerne un procédé de détermination des émissions de polluants d'un véhicule, le procédé est basé sur l'utilisation de mesures de la position et/ou de l'altitude et/ou de la vitesse du véhicule, au moyen de modèles du véhicule, du moteur et du système de post-traitement construits avec des paramètres macroscopiques. Ainsi, il est possible de déterminer les émissions de polluants 10 du véhicule sans instrumentation spécifique, pour un usage réel (grâce aux mesures), avec un temps de calcul limité (grâce aux modèles), et de manière précise (grâce aux modèles construits à l'aide de paramètres macroscopiques). Le procédé selon l'invention 15 La présente invention concerne un procédé de détermination des émissions de polluants d'un véhicule, ledit véhicule comportant un moteur à combustion interne et un système de post-traitement des gaz d'échappement dudit moteur, dans lequel on acquiert au moins un paramètre macroscopique relatif à la conception dudit véhicule, et dans lequel on construit pour ledit véhicule : 20 i) un modèle dudit véhicule qui relie ladite position et/ou l'altitude et/ou la vitesse dudit véhicule au couple et au régime dudit moteur au moyen d'au moins un paramètre macroscopique ; ii) un modèle dudit moteur qui relie ledit couple et ledit régime dudit moteur aux émissions de polluants en sortie dudit moteur au moyen d'au moins un paramètre 25 macroscopique ; et iii) un modèle dudit système de post-traitement qui relie lesdites émissions de polluants en sortie dudit moteur au moyen aux émissions de polluants en sortie dudit système de post-traitement au moyen d'au moins un paramètre macroscopique. 30 Pour ce procédé, on réalise les étapes suivantes : a) on mesure la position, et/ou l'altitude et/ou la vitesse dudit véhicule ; b) on détermine ledit couple et ledit régime dudit moteur au moyen dudit modèle de véhicule et desdites mesures ; c) on détermine les émissions de polluants en sortie dudit moteur au moyen dudit 35 modèle du moteur et dudit couple et dudit régime dudit moteur ; et 3049653 4 d) on détermine les émissions de polluants du véhicule au moyen dudit modèle du système de post-traitement et desdites émissions de polluants en sortie dudit moteur. Selon un mode de réalisation, on réalise un prétraitement desdites mesures avant l'étape de détermination dudit couple et dudit régime du moteur. 5 Avantageusement, on réalise ledit prétraitement au moyen d'un sur-échantillonnage et d'un filtrage. De préférence, on acquiert au moins un paramètre macroscopique choisi parmi : le type de motorisation, la norme d'homologation, la cylindrée du moteur, le couple maximal et le régime du moteur associé, la puissance maximale et le régime du moteur associé, la 10 masse du véhicule, le type de transmission du véhicule, le type de système de post- traitement, le type du système d'injection, l'architecture de la boucle d'air. Conformément à une mise en oeuvre, on acquiert lesdits paramètres macroscopiques depuis une base de données et/ou au moyen d'une interface avec un utilisateur. Selon une caractéristique, on construit ledit modèle du véhicule au moyen d'un modèle 15 de dynamique du véhicule reliant ladite position et/ou l'altitude et/ou ladite vitesse du véhicule à la puissance estimée dudit moteur au moyen d'au moins un paramètre macroscopique, et au moyen d'un modèle de transmission reliant ladite puissance moteur au régime et au couple du moteur au moyen d'au moins un paramètre macroscopique. Selon une variante, on construit ledit modèle du moteur au moyen d'un modèle 20 énergétique reliant ledit régime et ledit couple du moteur à des températures et débits de fluides mis en oeuvre par la combustion au sein dudit moteur au moyen d'au moins un paramètre macroscopique, et au moyen d'un modèle de polluants en sortie dudit moteur reliant lesdites températures et lesdits débits de fluides aux émissions de polluants en sortie dudit moteur au moyen d'au moins un paramètre macroscopique. 25 Conformément à une option de réalisation, on construit ledit modèle de polluants en sortie dudit moteur, au moyen d'un modèle quasi-statique et d'un modèle transitoire, ces deux modèles reliant lesdites températures et lesdits débits de fluides auxdites émissions de polluants en sortie dudit moteur au moyen d'au moins un paramètre macroscopique. De manière avantageuse, ledit modèle transitoire correspond à un coefficient correctif 30 dudit modèle quasi-statique. Conformément à un mode de réalisation, on construit ledit modèle de post-traitement par discrétisation dudit système de post-traitement en plusieurs tranches et au moyen de l'efficacité de chaque tranche discrétisée. Selon une mise en oeuvre, lesdites émissions du véhicule déterminées sont stockées 35 dans une base de données. Avantageusement, lesdits polluants sont choisis parmi les oxydes d'azote, les particules, les monoxydes de carbone et/ou les hydrocarbures non brûlés. 3049653 5 Conformément à une caractéristique, on mesure la position et/ou l'altitude et/ou la vitesse du véhicule au moyen d'un système de géolocalisation ou d'un téléphone portable. De préférence, on affiche les émissions de polluants déterminées sur un écran d'un 5 système de géolocalisation, d'un téléphone portable, d'un tableau de bord du véhicule ou sur un site internet. De plus, l'invention concerne un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou 10 exécutable par un processeur ou un serveur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des caractéristiques précédentes, lorsque ledit programme est exécuté sur un ordinateur ou sur un téléphone portable. En outre, l'invention concerne l'utilisation du procédé selon l'une des caractéristiques 15 précédentes, pour estimer l'efficacité écologique d'une infrastructure routière ou une réglementation de la circulation routière. Présentation succincte des figures D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la 20 lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. La figure 1 illustre les étapes du procédé selon un mode de réalisation de l'invention. La figure 2 est un histogramme comparant les émissions de NOx pour plusieurs trajets pour une référence et le procédé selon l'invention. 25 Description détaillée de l'invention La présente invention concerne la détermination des émissions de polluants d'un véhicule. Le terme « polluants » désigne les oxydes d'azote (NOx), les particules, les monoxydes de carbone (CO), les hydrocarbures non brûlés (HC). Le procédé selon 30 l'invention permet de déterminer les émissions d'au moins un, avantageusement plusieurs, et de préférence tous, de ces polluants. Le véhicule comporte un moteur à combustion interne (ci-après désigné par le terme « moteur ») et un système de post-traitement des gaz d'échappement du moteur. Le moteur à combustion interne peut être un moteur essence ou un moteur Diesel. Le moteur peut 35 réaliser à lui seul le déplacement du véhicule, ou peut faire partie d'un système de 3049653 6 propulsion hybride. Le système de post-traitement permet de traiter les émissions de polluants en sortie du moteur, diminuant ainsi les émissions de polluants du véhicule. Le système de post-traitement peut comporter un catalyseur d'oxydation pour traiter les hydrocarbures imbrûlés et le monoxyde de carbone, et/ou des catalyseurs DeNOx pour 5 réduire les oxydes d'azote en présence d'oxygène, et/ou différents filtres pour éliminer les particules solides. Pour le procédé selon l'invention, on acquiert au moins un paramètre macroscopique relatif à la conception du véhicule. On appelle paramètre macroscopique, une caractéristique 10 générale relative au véhicule, son moteur ou son système de post-traitement. Il s'agit de paramètre constant pour un véhicule, correspondant à des données du constructeur du véhicule. Le paramètre est dit macroscopique car il est déterminé à l'échelle du véhicule, et il ne s'agit pas de paramètre microscopique qui peut être déterminé, comme par exemple dans la demande de brevet FR2984557 (US 2013/0158967), à l'échelle d'une maille représentant 15 une petite portion de la chambre de combustion. Les paramètres macroscopiques permettent la construction de modèles macroscopiques représentatifs du véhicule. Les paramètres macroscopiques peuvent être de trois types : - des paramètres liés à construction générale du véhicule (par exemple : masse du véhicule, transmission, ...) 20 - des paramètres liés à la motorisation (par exemple : le type d'injection, la cylindrée, le type de motorisation...), et - des paramètres liés au système de post-traitement (par exemple : type de post-traitement). Selon un mode de réalisation de l'invention, on peut acquérir au moins un paramètre 25 macroscopique choisi parmi : - le type de motorisation (essence, Diesel...), - le niveau de norme d'homologation (Euro 1, Euro 2, ...) , - la cylindrée du moteur, - le couple maximal et le régime du moteur associé, 30 - la puissance maximale et le régime du moteur associé, - la masse du véhicule, - le type de transmission du véhicule (type et constitution de boîte de vitesse,...), - le type de système de post-traitement, - le type du système d'injection, 35 - l'architecture de la boucle d'air (présence/absence d'une recirculation des gaz brûlés noté EGR, utilisation d'un turbocompresseur, d'une suralimentation...), - les dimensions des roues, etc. 3049653 7 Selon une variante de réalisation, les paramètres macroscopiques peuvent être obtenus depuis une base de données, qui répertorie les différents véhicules en circulation. Par exemple, les paramètres macroscopiques peuvent être obtenus en indiquant le numéro d'immatriculation du véhicule, la base de données associant le numéro de plaque à sa 5 conception (marque, modèle, motorisation...), et comprenant les paramètres macroscopiques du véhicule. Alternativement, les paramètres macroscopiques peuvent être des données constructeurs renseignés par l'utilisateur, en particulier au moyen d'une interface (par exemple un téléphone intelligent ou un système de géolocalisation). 10 Notations Dans la suite de la description, on utilise les notations suivantes : posGps Coordonnées mesurées par géolocalisation dans le repère de Lambert [m] a/ tGps Altitude mesurée par géolocalisation 15 v - GPS Vitesse du véhicule mesurée par géolocalisation [m/s] Vitesse du véhicule [m/s] Ne Régime du moteur estimé [tr/min] Cme Couple du moteur estimé [Nm] m Masse du véhicule [kg] 20 FT Effort de traction du véhicule à la roue [N] Fres Résultante des efforts de frottement subis par le véhicule [N] Fslope Effort normal subi par le véhicule (gravité) [N] Fbrk Force de freinage mécanique [N] a Angle d'inclinaison de la route [rad] 25 a, b, c Paramètres du véhicule [-] Pe Puissance moteur estimée [kW] Tltrans Rendement de la transmission [-] RMTH_V Rapport de réduction entre le régime de rotation du moteur et la vitesse du véhicule [rpm/km/h] 30 PSME Emissions polluantes en sortie du moteur [g/s] PSME, Emissions du polluant i en sortie du moteur [g/s] PSMEi_os Emissions du polluant i en sortie du moteur pour un régime quasi- statique [g/h] NOx0s Masse de Nox par unité de masse de carburant [g/kg de carburant] 35 COC Angle de vilebrequin au centre de combustion (50% de conversion d'énergie) [°CA] 3049653 8 rricyl m02 ao, al, a2, a3 Sootos 5 BGR BGRdyn AFratio AFratio-dyn CO ri-QS2TR 10 polluant i PSEE PSEE; Convi polluant i 15 Téch Qéch Masse d'air enfermée dans le cylindre par cycle [gti o3cm3] Masse d'oxygène enfermée dans le cylindre par cycle [gti o3cm3] Coefficients Emissions de particules en sortie du moteur en régime quasi-statique [ g/s] Fraction de gaz brûlés dans le cylindre Fraction dynamique de gaz brûlés dans le cylindre Richesse du mélange dans le cylindre Richesse dynamique du mélange dans le cylindre Coefficient de correction pour l'impact des phénomènes transitoires pour le [-] Emissions polluantes en sortie du système de post-traitement [g/s] Emissions du polluant i en sortie du système de post-traitement [g/s] Efficacité de conversion de la tranche j du système de post-traitement pour le [-] Température des gaz d'échappement [K] Débit des gaz d'échappement [g/s] Pour ces notations, la dérivée par rapport au temps est notée cl. Dans la présente demande, le terme f désigne de manière générale une fonction, qui 20 peut être de tout type. Selon l'invention, on construit trois modèles macroscopiques afin de déterminer les émissions de polluants. Ces modèles sont paramétrés avec les paramètres macroscopiques acquis, permettant ainsi de rendre le procédé selon l'invention représentatif du véhicule. Il 25 s'agit d'un modèle du véhicule, d'un modèle du moteur et d'un modèle de post-traitement (c'est-à-dire un modèle du système de post-traitement). Ces modèles ne nécessitent pas d'instrumentation ou de matériel spécifiques sur le véhicule. Modèle de véhicule 30 Le modèle du véhicule relie la position et/ou l'altitude et/ou la vitesse du véhicule au couple et au régime du moteur, au moyen d'au moins un paramètre macroscopique. Selon une mise en oeuvre de l'invention, pour construire le modèle du véhicule on peut utiliser au moins un des paramètres macroscopiques suivants : masse du véhicule, puissance maximale et le régime du moteur associé, vitesse maximale, type de transmission.... 35 Selon un mode de réalisation de l'invention, on peut construire le modèle du véhicule par association d'un modèle de dynamique du véhicule et d'un modèle de transmission du 3049653 9 véhicule. Le modèle de dynamique du véhicule relie la position et/ou la vitesse et/ou l'altitude du véhicule à la puissance estimée du véhicule au moyen d'au moins un paramètre macroscopique, par exemple la masse du véhicule, le type de transmission, les dimensions des roues. Le modèle de transmission du véhicule relie la puissance du véhicule au régime 5 et au couple du moteur, au moyen d'au moins un paramètre macroscopique, par exemple le type de transmission, la puissance maximale et le régime moteur associé. Le modèle de dynamique du véhicule prend en compte la dynamique du véhicule. Il peut être construit à partir de l'application du principe fondamental de la dynamique du 10 véhicule appliqué sur son axe longitudinal, et peut s'écrire sous la forme suivante : dv m dt = FT - F's - Fsiope - Fbrk avec : v est la vitesse du véhicule, une entrée du modèle de la dynamique du véhicule. Fres peut être exprimée en fonction de la vitesse sous la forme F's = a + bv + cv2 avec 15 a, b, c des paramètres du véhicule à identifier en fonction des caractéristiques générales du véhicule (paramètres macroscopiques du véhicule). F dope peut être exprimée en fonction de la masse du véhicule et de l'inclinaison a de la route : Fsiope = mg sin(a). L'angle d'inclinaison a est une donnée d'entrée du modèle de la dynamique du véhicule. En effet, l'inclinaison a peut être calculée à partir de l'altitude et de la 20 distance parcourue, elle dépend donc de l'altitude et de la position. Selon un mode de réalisation, l'angle d'inclinaison a peut être déterminé par une formule du type : (Aaltitude a = arctan Adistance) Ces équations permettent d'écrire une formule qui relie la puissance estimée du moteur à la vitesse du véhicule et d'autres paramètres macroscopiques connus ou déterminables. En effet, on peut écrire l'équation : Pe = FTin *v, . trans 25 Ainsi, en combinant les différentes équations, il est possible de déterminer une formule qui relie la puissance du moteur à la vitesse et l'altitude du véhicule, au moyen de paramètres macroscopiques connus et constants. Le modèle de transmission estime le rapport de réduction entre le régime de rotation 30 du moteur thermique et la vitesse du véhicule. Il peut être paramétré en fonction des caractéristiques générales (paramètres macroscopiques) du véhicule, notamment la masse du véhicule, la puissance maximale, le type de transmission, en particulier le nombre de 3049653 10 rapports. Ce modèle de transmission utilise seulement la vitesse du véhicule en tant que donnée d'entrée, pour estimer le rapport de réduction : RMTH-v = f (v) La fonction f peut être obtenue notamment à partir d'abaques données par le constructeur. 5 Ce rapport de réduction peut ensuite être utilisé pour déterminer le régime du moteur. En effet, on peut écrire les relations suivantes : Ne = RMTH-v * V Ensuite, le couple du moteur peut être déterminé en fonction de la puissance (estimée à l'aide du modèle de dynamique du véhicule) et du régime du moteur : Cme = f(Ne,Pe) La fonction f peut être obtenue par des cartographies données par le constructeur. 10 Modèle du moteur Le modèle du moteur relie le régime et le couple du moteur aux émissions de polluants en sortie du moteur (c'est-à-dire avant le système de post-traitement), au moyen d'au moins un paramètre macroscopique. Selon une mise en oeuvre de l'invention, pour construire le 15 modèle du moteur on peut utiliser au moins un des paramètres macroscopiques suivants : la cylindrée, le type de motorisation, le couple et la puissance , l'architecture de la boucle d'air, la norme d'homologation du véhicule, etc. Selon un mode de réalisation de l'invention, on peut construire le modèle du moteur par association d'un modèle énergétique et d'un modèle de polluants en sortie du moteur. Le 20 modèle énergétique relie le couple et le régime du moteur à des débits et des températures de fluides mis en oeuvre dans le moteur à combustion (carburants, gaz à l'admission, gaz à l'échappement, éventuellement recirculation des gaz brûlés) au moyen d'au moins un paramètre macroscopique, par exemple la cylindrée, le type de motorisation, le couple et la puissance maximaux, l'architecture de la boucle d'air. Le modèle de polluants en sortie 25 moteur relie des débits et des températures de fluides mis en oeuvre dans le moteur à combustion interne aux émissions de polluants en sortie du moteur, au moyen d'au moins un paramètre macroscopique, par exemple la norme d'homologation du véhicule, le type de motorisation, l'architecture de la boucle d'air. 30 Le modèle énergétique permet d'estimer les grandeurs physiques sur le point de fonctionnement courant (régime, couple). Il est paramétré en fonction de paramètres macroscopiques. Les grandeurs physiques estimées sont les débits et températures des 3049653 11 fluides mis en oeuvre dans le moteur à combustion (carburants, gaz à l'admission, gaz à l'échappement, éventuellement recirculation des gaz brûlés). Le modèle de polluants en sortie du moteur permet à partir de l'information du régime 5 et du couple du moteur, et des estimations issues du modèle énergétique, d'estimer les émissions polluantes en sortie du moteur. Il peut être paramétré en fonction des caractéristiques générales du véhicule et du moteur : la norme d'homologation du véhicule, le type de motorisation, l'architecture de la boucle d'air, etc. L'estimation des polluants en sortie moteur peut se faire en deux étapes : 10 estimation des émissions en quasi-statique au moyen d'un modèle quasi-statique, et estimation de l'impact des phénomènes transitoires au moyen d'un modèle transitoire. Alternativement, l'estimation des polluants en sortie moteur peut se faire uniquement 15 en une seule étape au moyen du modèle quasi-statique. Estimer les émissions en quasi-statique d'un moteur sur un point de fonctionnement à un instant donné revient à considérer que ce moteur est en fonctionnement stabilisé sur ce point de fonctionnement. L'estimation de l'impact des phénomènes transitoires (fonctionnement non stabilisé) 20 permet de prendre en compte les phénomènes transitoires, qui engendrent généralement un surplus d'émissions polluantes. Les modèles quasi-statiques de polluants peuvent être paramétrés au moyen de paramètres macroscopiques du véhicule et du moteur. Ils permettent à chaque instant d'estimer les émissions polluantes quasi-statiquesen sortie moteur, à partir des estimations 25 de régime et de couple du moteur thermique et des sorties du modèle énergétique. Les modèles quasi-statiques peuvent s'écrire sous la forme : PSMEi_Qs = f (Ne, Cme) La fonction f peut être de différent type, en fonction du type de polluant étudié. Par exemple, le modèle quasi-statique de NOx peut être issu des travaux de Gartner, (U. Gartner, G. Hohenberg, H. Daudel and H. Oelschlegel, Development and Application of a 30 Semi-Empirical NOx Model to Various HD Diesel Engines), et peut s'écrire sous la forme : log(NOxQs) = ao + a1 * COC + a2 * mcy1 + a3 * mo2 Les coefficients ao, al, a2, a3 sont obtenus à partir de données expérimentales. Un des avantages de ce modèle est que ces coefficients varient peu d'un moteur à l'autre. Ce point est démontré dans l'article de Gartner précité. 3049653 12 Les particules en sortie du moteur est la combinaison de deux phénomènes : la formation et la post-oxydation dans la chambre de combustion. Ces phénomènes sont au premier ordre influencés par la richesse, le régime, la quantité de carburant, et le taux de gaz brûlés. Ainsi, le modèle statique de particules en sortie du moteur peut s'écrire sous une 5 équation de la forme : SootQs = f (AFratio,N e, Fuel, BGR) La fonction f peut être déterminée par corrélation avec des données expérimentales. Des modèles similaires peuvent être construits pour les autres polluants. Pour le mode de réalisation, pour lequel on détermine l'impact des phénomènes transitoires, on peut mettre en oeuvre en outre les moyens décrits ci-dessous. Les 10 phénomènes de dynamique de boucle d'air génèrent un écart sur les taux de BGR (fraction de gaz brûlés, liée à la recirculation des gaz d'échappement) et la richesse par rapport au point de fonctionnement stabilisé, ce qui a un impact fort sur les polluants, en particuliers les hydrocarbures HC, le monoxyde de carbone CO et particules. Les modèles d'impact du transitoire sont paramétrés en fonction de paramètres macroscopiques du moteur, en 15 particulier des caractéristiques de la boucle d'air récupérées (atmosphérique/suralimentée, recirculation des gaz brûlés haute pression EGRHP / recirculation des gaz brûlés basse pression EGRBP). Ces modèles permettent d'estimer les fractions de gaz brulés et richesses dynamiques à partir des estimations quasi-statiques et de la variation du couple estimé : BGRciyi, = f (BGR, Cme, dCme / dt) AFratio-dyn = f (AFratio,Cme, dCme / dt) 20 Un coefficient de correction pour chaque polluant peut être calculé en fonction de ces grandeurs dynamiques : COri_QS2TR = f (B GRdyn, BGR, AF - - ) ratto-dyn, AF rano Ces coefficients de correction permettent d'estimer les émissions polluantes en sortie moteur en prenant compte les phénomènes transitoires. Pour cela, les émissions de polluants en sortie du moteur peuvent s'écrire par une formule du type : PSMEi = Cori_ 25 Modèle du post-traitement Le modèle du post-traitement relie les émissions de polluants en sortie du moteur (c'est-à-dire avant le système de post-traitement) aux émissions de polluants en sortie du système de post-traitement, au moyen d'au moins un paramètre macroscopique. Selon une 30 mise en oeuvre de l'invention, pour construire le modèle du moteur on peut utiliser au moins QS2TR * PSMEi_QS 3049653 13 un des paramètres macroscopiques suivants : la cylindrée, la norme d'homologation du véhicule, etc. Le modèle de post-traitement peut comprendre des sous-modèles pour chaque technologie de dépollution, qui sont associés en fonction de l'architecture du système de 5 dépollution du véhicule. Ces sous-modèles peuvent être paramétrés en fonction de paramètres macroscopiques du véhicule tels que la norme d'homologation, la cylindrée, etc. Par exemple, les différentes technologies de dépollution peuvent être : TWC, de l'anglais « Three-way catalytic converters » signifiant catalyseur trois voies, 10 GPF (pour moteur essence), de l'anglais « gasoline particle filter » signifiant filtre à particules pour essence, DOC (pour moteur diesel), de l'anglais « Diesel oxidation catalyst », signifiant catalyseur à oxydation pour Diesel, DPF (pour moteur diesel), de l'anglais « Diesel particle filter », signifiant filtre à 15 particules pour Diesel, LNT (pour moteur Diesel), de l'anglais « lean Nox trap », signifiant piège à Nox, SCR (pour moteur Diesel), de l'anglais « Selective catalytic reduction », signifiant réduction catalytique sélective. Le modèle de post-traitement permet d'estimer les émissions de polluants en sortie du 20 système de post-traitement à partir des estimations de température, débits, et émissions polluantes en sortie du moteur. Le modèle de post-traitement peut être construit en discrétisant le système de post-traitement en plusieurs tranches (ou couches), et par association de l'efficacité Convij de chaque tranche discrétisée. Selon un exemple, le modèle de post-traitement peut s'écrire : Nb pain PSEEi = Convii(Téch,Qécn) * PSMEi Î=1 25 L'efficacité des tranches du système de post-traitement peut être déterminée à partir de cartographies du constructeur. Etapes du procédé de détermination des émissions de polluant Une fois les modèles construits pour le véhicule considéré, le procédé selon l'invention 30 comporte les étapes suivantes : 1) Mesure de géolocalisation 2) (étape facultative) Prétraitement des mesures 3) Détermination du couple et du régime moteur 4) Détermination des émissions de polluants en sortie du moteur 3049653 14 5) Détermination des émissions de polluants du véhicule 6) (étape facultative) Stockage des données La figure 1 illustre, de manière schématique et non limitative, les étapes du procédé 5 selon un mode de réalisation de l'invention. Sur cette figure, les traits pointillés indiquent des éléments facultatifs du procédé. Préalablement aux étapes du procédé, les différents modèles (modèle du véhicule MOD VEH, modèle du moteur MOD MOT et modèle du post-traitement MOD POT) sont construits. Ces modèles sont construits à partir de paramètres macroscopiques PAR. De 10 manière facultative, les paramètres macroscopiques PAR peuvent être obtenus depuis une base de données BDD, qui répertorie les différents véhicules en circulation. Par exemple, les paramètres macroscopiques PAR peuvent être obtenus en indiquant le numéro d'immatriculation du véhicule, la base de données BDD associant le numéro de plaque à la conception du véhicule (marque, modèle, motorisation...), et comprenant les paramètres 15 macroscopiques du véhicule. Une première série de paramètres macroscopiques PARI est utilisée pour la construction du modèle du véhicule MOD VEH. Cette première série de paramètres macroscopiques PARI peut comprendre les paramètres suivants : la masse du véhicule, la puissance maximale et le régime du moteur associé, la vitesse maximale, le type de 20 transmission (liste non limitative). Une deuxième série de paramètres macroscopiques PAR2 est utilisée pour la construction du modèle du moteur MOD MOT. Cette deuxième série de paramètres macroscopiques PAR2 peut comprendre les paramètres suivants : la cylindrée, le type de motorisation, le couple et la puissance maximums, l'architecture de la boucle d'air, la norme 25 d'homologation du véhicule (liste non limitative). Une troisième série de paramètres macroscopiques PAR3 est utilisée pour la construction du modèle du post-traitement MOD POT. Cette troisième série de paramètres macroscopiques PAR3 peut comprendre les paramètres suivants : la cylindrée, la norme d'homologation du véhicule (liste non limitative). 30 Ces trois modèles peuvent être construits selon l'une des variantes de réalisation décrites ci-dessus. La première étape consiste en une étape de mesure MES de géolocalisation. Lors de cette étape, on peut mesurer la position posGps, et/ou l'altitude aItGPS et/ou la vitesse vGps du 35 véhicule. La prise en compte de l'altitude aItGPS permet notamment de prendre en compte la pente de la route. De préférence, on réalise les trois mesures, de manière à avoir l'information la plus précise possible concernant la géolocalisation du véhicule, car on peut 3049653 15 prendre alors en compte le style de conduite, et l'accélération du véhicule. Cette mesure peut être réalisée à l'aide d'un système de géolocalisation, par exemple de type GPS (pour « global positioning system » signifiant géo-positionnement par satellite) ou de type Galileo, ou au moyen d'un téléphone intelligent (« smartphone »), etc. Dans le cas d'un téléphone 5 intelligent, celui-ci peut être équipé d'un système de géolocalisation, alternativement les mesures peuvent être mises en oeuvre par d'autres moyens, notamment par triangulation. La deuxième étape, qui est une étape facultative, est une étape de prétraitement PRT des signaux de mesure. Cette étape permet d'améliorer la qualité des signaux mesurés 10 avant de les utiliser. Cette étape peut notamment être intéressante, si les mesures sont réalisées à partir d'un téléphone intelligent, car les mesures issues d'un tel dispositif peuvent être quelque peu imprécises. Ce prétraitement peut être variable, car il est dépendant de la qualité des données d'entrée. Selon un mode de réalisation de l'invention, le prétraitement PRT peut comprendre un sur-échantillonnage des signaux, puis un filtrage. En sortie de 15 cette étape, on retrouve donc des signaux relatifs à la position posGps, et/ou l'altitude altGps et/ou la vitesse vGps du véhicule, ces signaux ayant été prétraités. La troisième étape concerne la détermination du couple et du régime moteur. Cette étape est mise en oeuvre au moyen du modèle du véhicule MOD VEH, qui détermine le 20 couple Cme et le régime Ne du moteur, en fonction des données de géolocalisation : la position posGps, et/ou l'altitude aItGPS et/ou la vitesse vGps du véhicule. La quatrième étape concerne la détermination des émissions polluantes en sortie du moteur, cette étape est mise en oeuvre au moyen du modèle du moteur MOD MOT, qui 25 détermine les émissions de polluants en sortie du moteur PSME, en fonction du couple Cme et du régime Ne du moteur. La cinquième étape concerne la détermination des émissions de polluants du véhicule, c'est-à-dire en sortie du système de post-traitement. La détermination des émissions de 30 polluants peut être réalisée à chaque instant, par exemple à une fréquence de 1Hz. De plus, il est également possible de déterminer le cumul des émissions de polluants sur un trajet réalisé. Cette étape est mise en oeuvre au moyen du modèle de post-traitement MOD POT, qui détermine les émissions de polluants en sortie du système de post-traitement PSEE, en fonction des émissions de polluants en sortie du moteur PSME. 35 La sixième étape, facultative, concerne le stockage des données. Une fois que les émissions de polluants du véhicule PSEE sont déterminées, celles-ci peuvent être stockées 3049653 16 STO (enregistrées), en particulier dans une base de données (différente de la base de données qui comportent les paramètres macroscopiques). Ce stockage STO peut concerner uniquement les émissions de polluants du véhicule PSEE, mais peut aussi concerner les données déterminées après chaque étape du procédé : les mesures prétraitées et/ou le 5 couple Cme et le régime Ne du moteur et/ou les émissions de polluants en sortie du moteur PSME. Ces informations permettent un suivi des usages réels et des émissions associées avec une bonne résolution spatiale et temporelle. Ces informations peuvent par exemple permettre d'évaluer la pertinence environnementale des infrastructures routières à l'échelle d'une rue, identifier des pics d'émissions localisés, identifier l'impact du style de la conduite 10 sur les émissions, etc. Cette base de données peut associer les émissions instantanées aux données cartographiques pour constituer une carte des émissions polluantes réelles. Il est donc possible de tirer des conclusions à l'échelle d'une portion de route, d'un trajet complet ou même d'une zone géographique selon les besoins. Lors de cette étape, il est également 15 possible d'afficher les émissions de polluants en sortie du véhicule PSEE, par exemple sur un écran d'un système de géolocalisation (GPS, Galiléo), d'un téléphone intelligent, sur le tableau de bord du véhicule, sur un site internet, etc. Ainsi, il est possible d'informer l'utilisateur ou tout autre personne (par exemple un gestionnaire de flotte de véhicules, un responsable d'infrastructure routière...) de la pollution émise sur un trajet ou sur une route. 20 Le procédé selon l'invention peut être utilisé pour les véhicules automobiles. Toutefois, il peut être utilisé dans le domaine du transport routier, le domaine des deux-roues, le domaine ferroviaire, le domaine naval, le domaine aéronautique, le domaine des aéroglisseurs, et le domaine des véhicules amphibies... 25 L'invention concerne par ailleurs un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur ou un serveur. Ce programme comprend des instructions de code de programme pour la mise en oeuvre du procédé tel que décrit ci- 30 dessus, lorsque le programme est exécuté sur un ordinateur ou un téléphone portable. Le procédé selon l'invention présente donc les avantages suivants : Le procédé selon l'invention permet la prise en compte de la pente, de l'accélération et du style de conduite, pour déterminer de manière précise les 35 émissions polluantes, 3049653 17 La modélisation des polluants selon l'invention est réalisée par une approche physique nécessitant uniquement la connaissance de paramètres macroscopiques, disponibles facilement, par exemple à partir de la plaque d'immatriculation d'un véhicule, 5 Le procédé selon l'invention est adapté aux paramètres réels de chacun des véhicules, Le procédé selon l'invention ne nécessite pas de données mesurées au banc moteur, qui seraient spécifiques à un véhicule, Le procédé selon l'invention permet de modéliser une classe de véhicule, 10 Le procédé selon l'invention ne nécessite pas d'instrumentation ou de matériel spécifique, en effet, les émissions de polluants peuvent être déterminées au moyen d'un dispositif de géolocalisation ou d'un téléphone intelligent, Le procédé selon l'invention permet de placer les émissions polluantes sur les routes traversées, 15 L'invention permet d'estimer l'efficacité écologique d'une infrastructure routière existante, sans nécessiter l'installation coûteuse de capteurs de polluants spécifiques, Elle permet également d'observer l'impact d'une réglementation portant sur la circulation ou sur les véhicules sur les émissions polluantes en usage réel. 20 Exemple Afin de montrer l'intérêt du procédé selon l'invention, on réalise un exemple pour comparer les émissions de NOx mesurées et les émissions estimées au moyen des modèles mis en oeuvre dans le procédé selon l'invention. La figure 2 illustre une comparaison des 25 résultats obtenus, en termes d'émissions de NOx, pour 26 trajets différents, indiqués par leur numéro (N° trajet). Sur cet histogramme, la référerce REF est issue d'une campagne d'essais, et les valeurs estimées sont obtenues par le procédé selon l'invention INV (en utilisant les variantes de réalisation des modèles détaillés ci-dessus). On observe une bonne cohérence des résultats. 30 REVENDICATIONS1) Procédé de détermination des émissions de polluants d'un véhicule, ledit véhicule comportant un moteur à combustion interne et un système de post-traitement des gaz d'échappement dudit moteur, dans lequel on acquiert au moins un paramètre macroscopique (PAR) relatif à la conception dudit véhicule, et dans lequel on construit pour ledit véhicule : i) un modèle dudit véhicule (MOD VEH) qui relie ladite position et/ou l'altitude et/ou la vitesse dudit véhicule au couple et au régime dudit moteur au moyen d'au moins un paramètre macroscopique (PAR) ; ii) un modèle dudit moteur (MOD MOT) qui relie ledit couple et ledit régime dudit moteur aux émissions de polluants en sortie dudit moteur au moyen d'au moins un paramètre macroscopique (PAR) ; et iii) un modèle dudit système de post-traitement (MOD POT) qui relie lesdites émissions de polluants en sortie dudit moteur au moyen aux émissions de polluants en sortie dudit système de post-traitement au moyen d'au moins un paramètre macroscopique (PAR) ; caractérisé en ce qu'on réalise les étapes suivantes : a) on mesure la position (posGps), et/ou l'altitude (aItGPS) et/ou la vitesse (vGps) dudit véhicule ; b) on détermine ledit couple (Cme) et ledit régime (Ne) dudit moteur au moyen dudit modèle de véhicule (MOD VEH) et desdites mesures ; c) on détermine les émissions de polluants en sortie dudit moteur (PSME) au moyen dudit modèle du moteur (MOD MOT) et dudit couple (Cme) et dudit régime (Ne) dudit moteur ; et d) on détermine les émissions de polluants du véhicule (PSEE) au moyen dudit modèle du système de post-traitement (MOD POT) et desdites émissions de polluants en sortie dudit moteur (PSME). 2) Procédé selon la revendication 1, dans lequel on réalise un prétraitement (PRT) desdites mesures avant l'étape de détermination dudit couple (Cme) et dudit régime (Ne) du moteur. 3) Procédé selon la revendication 2, dans lequel on réalise ledit prétraitement (PRT) au moyen d'un sur-échantillonnage et d'un filtrage. 4) Procédé selon l'une des revendications précédentes, dans lequel on acquiert au moins un paramètre macroscopique (PAR) choisi parmi : le type de motorisation, la norme 3049653 19 d'homologation, la cylindrée du moteur, le couple maximal et le régime du moteur associé, la puissance maximale et le régime du moteur associé, la masse du véhicule, le type de transmission du véhicule, le type de système de post-traitement, le type du système d'injection, l'architecture de la boucle d'air. 5 5) Procédé selon l'une des revendications précédentes, dans lequel on acquiert lesdits paramètres macroscopiques (PAR) depuis une base de données et/ou au moyen d'une interface avec un utilisateur. 10 6) Procédé selon l'une des revendications précédentes, dans lequel on construit ledit modèle du véhicule (MOD VEH) au moyen d'un modèle de dynamique du véhicule reliant ladite position et/ou l'altitude et/ou ladite vitesse du véhicule à la puissance estimée dudit moteur au moyen d'au moins un paramètre macroscopique (PAR), et au moyen d'un modèle de transmission reliant ladite puissance moteur au régime et au couple du 15 moteur au moyen d'au moins un paramètre macroscopique (PAR). 7) Procédé selon l'une des revendications précédentes, dans lequel on construit ledit modèle du moteur (MOD MOT) au moyen d'un modèle énergétique reliant ledit régime et ledit couple du moteur à des températures et débits de fluides mis en oeuvre par la 20 combustion au sein dudit moteur au moyen d'au moins un paramètre macroscopique (PAR), et au moyen d'un modèle de polluants en sortie dudit moteur reliant lesdites températures et lesdits débits de fluides aux émissions de polluants en sortie dudit moteur au moyen d'au moins un paramètre macroscopique (PAR). 25 8) Procédé selon la revendication 7, dans lequel on construit ledit modèle de polluants en sortie dudit moteur, au moyen d'un modèle quasi-statique et d'un modèle transitoire, ces deux modèles reliant lesdites températures et lesdits débits de fluides auxdites émissions de polluants en sortie dudit moteur au moyen d'au moins un paramètre macroscopique (PAR). 30 9) Procédé selon la revendication 8, dans lequel ledit modèle transitoire correspond à un coefficient correctif dudit modèle quasi-statique. 10) Procédé selon l'une des revendications précédentes, dans lequel on construit ledit 35 modèle de post-traitement (MOD POT) par discrétisation dudit système de post- traitement en plusieurs tranches et au moyen de l'efficacité de chaque tranche discrétisée. 3049653 20 11) Procédé selon l'une des revendications précédentes, dans lequel lesdites émissions du véhicule déterminées sont stockées (STO) dans une base de données. 5 12) Procédé selon l'une des revendications précédentes, dans lequel lesdits polluants sont choisis parmi les oxydes d'azote, les particules, les monoxydes de carbone et/ou les hydrocarbures non brûlés. 13) Procédé selon l'une des revendications précédentes, dans lequel on mesure la position 10 (posGps) et/ou l'altitude (aItGPS) et/ou la vitesse (vGps) du véhicule au moyen d'un système de géolocalisation ou d'un téléphone portable. 14) Procédé selon l'une des revendications précédentes, dans lequel on affiche les émissions de polluants (PSEE) déterminées sur un écran d'un système de 15 géolocalisation, d'un téléphone portable, d'un tableau de bord du véhicule ou sur un site internet. 15) Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur 20 ou un serveur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur ou sur un téléphone portable. 16) Utilisation du procédé selon l'une des revendications 1 à 14 pour estimer l'efficacité 25 écologique d'une infrastructure routière ou une réglementation de la circulation routière.