Brevet FR3047272A1 — publié le 20170804

Domaine de l'invention La présente invention se rapporte à un procédé de déter- mination d'un débit massique d'ammoniac entre deux catalyseurs SCR installés l'un derrière l'autre dans une conduite de gaz d'échappement faisant partie d'un système de catalyseurs SCR qui ne comporte qu'une unité de dosage d'agent réducteur en amont du premier catalyseur SCR. L'invention se rapporte également à un programme d'ordinateur exécutant les étapes du procédé ainsi qu'à un support de mémoire lisible par une machine et contenant le programme de l'ordinateur. L'invention a également pour objet un appareil de commande électronique pour exécuter le procédé. Etat de la technique La réduction catalytique sélective SCR constitue un pro- cédé promoteur pour réduire la teneur en oxydes d'azote des gaz d'échappement riches en oxygène avec de l'ammoniac ou des réactifs décomposant l'ammoniac. Le rendement d'un catalyseur SCR dépend de sa température, de la vitesse des gaz d'échappement et d'une manière très importante, du niveau de remplissage de l'ammoniac adsorbé à sa surface. Comme pour réduire les oxydes d'azote, on dispose à côté de l'ammoniac dosé directement également de l'ammoniac adsorbé, le rendement d'un catalyseur SCR chargé est supérieur à celui d'un catalyseur vide. Le comportement d'accumulation dépend de la température de fonctionnement du catalyseur. Plus cette température est faible, et plus grande est sa capacité d'accumulation. Lorsque l'accumulateur d'un catalyseur SCR est complètement rempli, on pourra le réduire par des variations brusques de la charge du moteur à combustion interne dont les gaz d'échappement seront réduits par le catalyseur SCR, même avec un glissement d'ammoniac, si l'on ne dose plus d'ammoniac ou pas de réactif donnant de l'ammoniac dans la conduite des gaz d'échappement. Pour arriver à des conversions d'oxydes d'azote aussi élevées que possible, il est toutefois indispensable de faire fonctionner le système SCR avec un niveau de remplissage élevé en ammoniac. Mais, si alors, du fait d'une varia- tion brusque de charge du moteur à combustion interne, la température 3047272 2 du catalyseur SCR complètement chargée augmente, sa capacité d'accumulation de l'ammoniac diminue ce qui se traduit par un glissement d'ammoniac. Cet effet est particulièrement accentué du fait que les ca- 5 talyseurs SCR sont installés à proximité du moteur à combustion in- terne et qu'ainsi le catalyseur SCR atteint rapidement sa température de fonctionnement après un démarrage à froid du moteur à combustion interne. Un second catalyseur SCR installé en aval du premier catalyseur SCR peut ainsi être prévu dans la conduite des gaz d'échappement 10 pour absorber le glissement d'ammoniac à travers le premier catalyseur et finalement convertir cet ammoniac. Les directives relatives au diagnostique embarqué (dia- gnostique OBD) exigent la surveillance des deux catalyseurs SCR. Pour cela, en général, à l'aval des deux catalyseurs SCR, on a un détecteur 15 d'oxydes d'azote. Pour des raisons de coûts, on n'installe habituellement qu'une soupape de dosage en amont du premier catalyseur SCR, pour doser la solution d'agent réactif décomposant l'ammoniac dans la conduite des gaz d'échappement. Le remplissage en ammoniac du second catalyseur SCR ne se fait, dans ces conditions, que par le glissement 20 d'ammoniac à travers le premier catalyseur SCR. On peut utiliser les données de tels catalyseurs pour modéliser le niveau de remplissage des deux catalyseurs SCR. Toutefois, en cas d'écart pour modéliser le vieillissement des catalyseurs SCR, les niveaux de remplissage réels peuvent différer de façon significative des niveaux de remplissage 25 modélisés. Cela se traduit par des modifications du rendement de la ré- duction des oxydes d'azote et ainsi le cas échéant au dépassement des limites d'émission. Exposé et avantages de l'invention Le procédé selon l'invention sert à déterminer un débit 30 massique d'ammoniac entre deux catalyseurs SCR installés l'un à la suite de l'autre dans la conduite des gaz d'échappement et faisant partie d'un système de catalyseurs SCR ne comportant qu'une unité de dosage d'agent réducteur en amont du premier catalyseur SCR. On détermine à partir d'un signal fourni par un détecteur d'oxydes d'azote NOx installé 35 entre les deux catalyseurs SCR et d'un signal d'un détecteur d'oxydes 3047272 3 d'azote NOx installé en amont du second catalyseur SCR. On utilise pour cela le fait que le détecteur d'oxydes d'azote installé dans un tel système de catalyseur SCR ainsi qu'entre les deux catalyseurs SCR réagit en sensibilité transversale à l'ammoniac et c'est pourquoi il mesure 5 non seulement les oxydes d'azote non réduits dans le premier cataly- seur SCR, mais également l'ammoniac provenant du glissement d'ammoniac sur le premier catalyseur SCR. Un procédé de calcul approprié répartit le signal fourni par le capteur d'oxydes d'azote NOx installé entre les deux catalyseurs SCR, en ammoniac et en oxydes d'azote 10 pour déterminer ainsi le débit massique d'oxydes d'azote. Selon un développement simple du procédé, on utilise la différence entre les signaux de concentration fournis par les deux capteurs comme base pour déterminer le débit massique d'ammoniac. Ensuite, il suffit de convertir les concentrations massiques. 15 Selon un développement plus complexe du procédé, à partir du signal du détecteur d'oxydes d'azote NOx installé entre les deux catalyseurs SCR et le signal du second détecteur d'oxydes d'azote NOx installé en aval du second catalyseur SCR, on obtient une demande d'agent de dosage spécifique du système de catalyseur SCR et 20 une demande d'agent de dosage spécifique du premier catalyseur SCR. Le débit massique d'ammoniac peut s'obtenir par la comparaison de l'intégrale des deux demandes d'agent de dosage en fonction du temps. La surface comprise entre les deux intégrales est ainsi proportionnelle au glissement d'ammoniac en aval du premier catalyseur SCR. Le 25 moyen de dosage spécifique est le quotient de la masse d'ammoniac do- sée dans le premier catalyseur SCR et de la masse d'oxydes d'azote convertie dans le premier catalyseur SCR ou dans l'ensemble du système de catalyseur SCR. Pour le calcul de la masse d'oxydes d'azote convertie, on calcule toujours avec la masse molaire du dioxyde d'azote, c'est- 30 à-dire 46g/mol. Si la masse d'ammoniac effectivement dosée diffère de cette valeur à cause des actions de régulation de la masse d'ammoniac à doser selon le modèle, cela suppose que la quantité d'oxydes d'azote convertie était élevée. Pour le calcul des deux besoins en agent de dosage on 35 tient compte de préférence de la masse d'ammoniac dégagée de par 3047272 4 l'agent réducteur et qui est dosée par l'unité de dosage d'agent réducteur dans la conduite des gaz d'échappement et de la masse des oxydes d'azote en amont du premier catalyseur SCR. La masse d'ammoniac dosée se détermine par la commande de l'unité de dosage d'agent réduc- 5 teur. La masse des oxydes d'azote en amont du premier catalyseur SCR se détermine avec le capteur d'oxydes d'azote NOx installé à cet endroit ou avec un modèle. Dans les formes de réalisation du procédé, décrites ci- dessus, la détermination se fait de préférence à un point de fonction- 10 nement du système de catalyseur SCR, point où une variation des ni- veaux de remplissage en ammoniac des deux catalyseurs SCR ne dépasse pas un seuil de variation. Cela augmente la fiabilité du débit massique d'ammoniac obtenu si approximativement on peut supposer que dans les deux catalyseurs SCR il n'y a pas eu absorption ni désorp- 15 tion d'ammoniac et que cela ne participe pas au débit massique d'ammoniac dans le système de catalyseur SCR. Suivant une autre caractéristique, dans ces formes de réalisation du procédé, de manière préférentielle, on détermine un point de fonctionnement du système de catalyseur SCR là où le glissement 20 d'ammoniac par le second catalyseur SCR ne dépasse pas un seuil de glissement. D'une manière particulièrement préférentielle, ce seuil de glissement est égal à zéro. En principe, on utilise des systèmes de catalyseur SCR à deux catalyseurs SCR de façon à ce que le glissement d'ammoniac en aval du second catalyseur SCR n'existe pas car cet am- 25 moniac s'échapperait à l'air ambiant aux points de fonctionnement pour lesquels on peut estimer que l'objectif est atteint ; pour déterminer le débit massique d'ammoniac on peut estimer qu'entre les deux catalyseurs SCR, le débit massique d'ammoniac en aval du second catalyseur SCR est pratiquement égal à zéro. Si cette condition n'est pas remplie, 30 cela détériore la fiabilité du débit massique d'ammoniac obtenu. Selon un autre exemple de réalisation du procédé, on détermine le débit massique d'ammoniac d'une manière particulièrement précise pour déterminer le débit massique d'ammoniac du rendement du second catalyseur SCR. Ce rendement peut ainsi être pris comme 35 modèle. 3047272 5 Les niveaux de remplissage en ammoniac des deux catalyseurs SCR peuvent être corrigés par le débit massique d'ammoniac que l'on détermine selon toutes les formes de réalisation du procédé, pour minimiser les écarts entre les niveaux de remplissage physiques et 5 les niveaux de remplissage modélisés d'ammoniac. Le programme d'ordinateur selon l'invention est conçu pour exécuter chaque étape du procédé, notamment lorsque le programme est exécuté par un calculateur ou un appareil de commande électronique. Il permet d'implémenter le procédé avec un appareil de 10 commande électronique usuel sans avoir à faire pour cela des modifica- tions constructives. Le programme est enregistré sur une mémoire lisible par une machine. L'exécution du programme par un appareil de commande électronique usuel consiste à déterminer par le procédé, un débit massique d'ammoniac entre deux catalyseurs SCR installés l'un 15 après l'autre dans la conduite des gaz d'échappement et faisant partie d'un système de catalyseurs SCR, et qui ne comportent qu'une unité de dosage d'agent réducteur en amont du premier catalyseur SCR. Dessins La présente invention sera décrite ci-après, de manière 20 plus détaillée à l'aide d'exemples de procédés de détermination d'un dé- bit massique d'ammoniac dans une conduite de gaz d'échappement d'un moteur, représentés dans les dessins annexés dans lesquels : la figure 1 est un schéma d'un système de catalyseurs SCR comportant deux catalyseurs SCR dont on détermine le débit massique 25 d'ammoniac entre les deux catalyseurs SCR à l'aide d'un procédé selon l'un des modes de réalisation de l'invention, la figure 2 montre le chronogramme de la demande spécifique d'agent de dosage d'un système de catalyseurs SCR comportant deux catalyseurs SCR et du moyen de dosage spécifique à une ap- 30 plication de l'un des deux catalyseurs SCR, selon un procédé cor- respondant à l'exemple de réalisation de l'invention, la figure 3 montre le chronogramme de la concentration en ammoniac entre deux catalyseurs SCR selon un exemple de réalisation du procédé de l'invention. 35 3047272 6 Description de modes de réalisation La conduite des gaz d'échappement 11 d'un moteur à combustion interne 10 comporte un système de catalyseurs SCR 20 comme le montre la figure 1. Il dispose d'une unité de dosage d'agent 5 réducteur 21 qui injecte une solution aqueuse d'urée dans la conduite des gaz d'échappement 11. Aux températures élevées des gaz d'échappement, cette solution d'urée libère de l'ammoniac. En aval de l'unité de dosage d'agent réducteur 21, il y a un premier catalyseur SCR 22 et un second catalyseur SCR 23. Un premier détecteur d'oxydes 10 d'azote NOx 31 est installé en amont de l'unité de dosage d'agent réduc- teur 21 dans la conduite des gaz d'échappement 11. Un second détecteur d'oxydes d'azote NOx 32 est installé entre les deux catalyseurs SCR 22, 23. Un troisième détecteur d'oxydes d'azote NOx est installé en aval du second catalyseur SCR 23. Tous les détecteurs d'oxydes d'azote NOx 15 31, 32, 33 alimentent en signaux un appareil de commande électro- nique 40. Comme les détecteurs d'oxydes d'azote NOx 31, 32, 33 ont une sensibilité transversale vis-à-vis de l'ammoniac, les signaux qu'ils fournissent sont des signaux somme correspondant à la somme des oxydes d'azote et de l'ammoniac. Le premier détecteur d'oxydes d'azote 20 NOx est toutefois installé en amont de l'unité de dosage d'agent réduc- teur 21 de sorte qu'il mesure de façon fiable la quantité d'oxydes d'azote contenue dans les gaz d'échappement. Si l'on fait fonctionner le système de catalyseur SCR 20 pour ne pas avoir de glissement d'ammoniac sur le second catalyseur SCR 23, on peut supposer que le signal du troi- 25 sième détecteur d'oxydes NOx mesure exclusivement les oxydes d'azote. Comme il est prévu un glissement d'ammoniac sur le premier catalyseur SCR 22, pour alimenter le second catalyseur SCR 23 en ammoniac, le second détecteur d'oxydes d'azote NOx fournit toujours un signal somme pour l'ammoniac et les oxydes d'azote. L'unité de dosage 30 d'agent réducteur 21 signale la quantité d'ammoniac dosée dans la con- duite des gaz d'échappement 11 également à l'appareil de commande 40. Un premier exemple du procédé selon l'invention repose sur la connaissance de la différence entre la concentration r(NH3) 21 de 35 l'ammoniac qui a été libéré à partir de la solution d'agent réducteur do- 3047272 7 sée par l'unité de dosage d'agent réducteur 21 et de la concentration en ammoniac r(NH3) 32 sur le second détecteur d'oxydes d'azote NOx 32 selon la formule 1 de la différence entre la concentration en oxydes d'azote r(NOx)_31 en amont du premier catalyseur SCR 22 et de la con- 5 centration en oxydes d'azote r(NOx)_32 entre les deux catalyseurs SCR 22 et 23 : r(NH3) 21 - r(NH3) 32 = r(NOx)_31 - r(NOx)_32 (Formule 1) 10 La valeur de la concentration r(sens) 32 transmise par le second détecteur d'oxydes d'azote NOx 32 à l'appareil de commande 40 correspond ainsi, selon la formule 2 à la somme des concentrations en oxydes d'azote et en ammoniac entre les deux catalyseurs SCR 22, 23 : 15 r(sens) 32 = r(NOx)_32 + r(NH3) 32 (Formule 2) Selon la formule 3, la concentration d'oxydes d'azote r(NH3) 21 dosée correspond à la différence entre la concentration en oxydes d'azote en amont du premier catalyseur 22 et la concentration 20 en oxydes d'azote (NOx) _33 en aval du second catalyseur SCR 23 : r(NH3) 21 = r(NOx)_31 - r(NOx)_33 (Formule 3) A partir des formules 1 à 3 on obtient la formule 4 qui 25 permet de calculer la concentration en ammoniac r(NH3) 32 entre les deux catalyseur SCR 22, 23 à partir du signal de détecteur fourni par le second détecteur NOx 32 et du troisième détecteur NOx 33 : 30 r(NH 3) _32 - r(sens) _ 32 - r(NOx) _33 (Formule 4) 2 S'il n'y a pas de glissement d'ammoniac sur le second catalyseur SCR 23, on pourra prendre la concentration en oxydes d'azote r(NOx)_33 en aval du second catalyseur SCR 23 à partir du signal du troisième détecteur NOx 33. En déterminant les deux valeurs de con- 3047272 8 centration avec le temps, on peut calculer le débit massique d'ammoniac. Selon un second exemple de réalisation du procédé de l'invention, en appliquant la formule 5, on calcule la demande spéci- 5 figue d'agent de dosage facDos_22 du premier catalyseur SCR 22. Dans cette formule m(NH3)_21 désigne la masse d'ammoniac introduite par dosage et la différence entre la masse d'oxydes d'azote m(NOx)_31 en amont du premier catalyseur SCR 22 et la masse d'oxydes d'azote m(NOx)_32 en aval du premier catalyseur SCR correspond à la masse 10 d'oxydes d'azote convertie dans le premier catalyseur SCR : Dos 22 = m(NH3)_21 (Formule 5) m(NOx)_31 -(m(N0x)_32 + m(N 1-13 )_ 32) La masse d'oxydes d'azote m(NOx)_32 et la masse 15 d'ammoniac m(NH3)_32 ne peuvent être mesurées directement. Comme le signal de masse m(sens)_32 fourni par le second détecteur d'oxydes d'azote 32 correspond à la somme de la masse d'oxydes d'azote m(NOx)_32 et la masse d'ammoniac m(NH3)_32 entre les deux catalyseurs SCR 22, 23, on peut remplacer la formule 5 par la formule 6 sui- 20 vante : Dos 22 = m(NH3)_21 (Formule 6) m(NOx)_31 - m(sens)_32 Pour calculer la demande spécifique de l'agent de dosage 25 facDos_20 du système de catalyseur SCR 20 il faut, selon la formule 7, la masse d'oxydes d'azote mesurée par le troisième détecteur d'oxydes d'azote 33 en aval du second catalyseur SCR 22 : 30 m(NH3)_21 (Formule 7) Dos 20- m(NOx)_31 - m(NOx)_33 On suppose qu'en aval du second catalyseur SCR 23 il n'y a pas de glissement d'ammoniac de sorte que le signal fourni par le troisième détecteur d'oxydes d'azote NOx 33 n'est pas faussé par la présence d'ammoniac. De plus, partant le la réduction des oxydes d'azote 3047272 9 dans le premier catalyseur SCR 22 on peut établir la formule 8 suivante pour l'ensemble du système de catalyseur SCR 20 : 5 Dos-20 = (Formule 8) m(NH3)_21 - m(NH3)_32 m(NOx)_31 - m(NOx)_32 Avec les formules 6 à 8 et avec les données des capteurs on obtient les demandes d'agent de dosage facDos 20, facDos 22. Pour cela, il faut tenir compte de ce que les masses molaires d'ammoniac et des oxydes d'azote sont différentes. Pour la masse d'oxydes d'azote on 10 utilise toujours pour le calcul, la masse molaire du dioxyde d'azote. La figure 2 montre son évolution dans un exemple d'état de fonctionnement du moteur à combustion interne 10 en fonction du temps t. La différence entre l'intégrale en fonction du temps des deux demandes d'agent de dosage facDos 20, facDos 22 est proportionnelle au débit 15 massique d'ammoniac entre les deux catalyseurs SCR 22, 23. On a un débit massique d'ammoniac élevé et ainsi un glissement d'ammoniac élevé sur le premier catalyseur SCR 22 là où les deux chronogrammes des deux demandes spécifiques d'agent de dosage facDos 20, facDos 22 diffèrent de manière significative. Là où ces courbes se rapprochent, le 20 débit massique d'ammoniac est faible. Cela a été vérifié en ce que l'on a installé un détecteur d'ammoniac non représenté, entre les deux catalyseurs SCR 22, 23 dans la conduite des gaz d'échappement 11. L'évolution chronologique de la concentration en ammoniac r(NH3) 32 mesurée est représentée à la figure 3. Le calcul de la masse d'ammoniac 25 m(NH3) 32 s'obtient à partir des valeurs mesurées selon la formule 9 suivante : 30 m(NH3)_21 m(NH3)_21 (Formule 9) facDos_20 facDos_22 m(NH3)_32 - 1 +1 facDos _20 Dans un troisième mode de réalisation du procédé de l'invention on utilise le fait que, selon la formule 10, le débit massique 3047272 10 d'oxydes d'azote dm(NOx)_33 en aval du second catalyseur SCR 23 correspond à la différence entre le débit massique d'oxydes d'azote dm(NOx)_32 entre les deux catalyseurs SCR 22, 23 et le débit massique du glissement d'ammoniac dm(NH3) 32 du premier catalyseur SCR 22, 5 différence à laquelle on ajoute la différence Am(NH3) entre la masse d'ammoniac fournie au second catalyseur SCR 23 et celle qu'il consomme : dm(NOx)_33 = dm(NOx)_32 - dm(NH3) 32 + 4m(NH3) (Formule 10) 10 Dans ce cas il faut également tenir compte de la différence des masses molaires de l'ammoniac et des oxydes d'azote. La différence Am(NH3) correspond, selon la formule 11, par ailleurs au débit massique d'oxydes d'azote dm(NOx)_32 entre les deux catalyseurs SCR 15 22, 23 et du débit massique d'oxydes d'azote dm(NOx)_32 multiplié par le rendement ETA(23) du second catalyseur SCR, entre les deux catalyseurs SCR 22, 23 : 4m(NH3) = dm(NOx)_32 - ETA(23) x dm(NOx)_32 (Formule 11) 20 A partir des formules 10 et 11 on obtient la formule sui- vante : dm(NOx)_33 = (1 - ETA(23)) x dm(NOx)_32 (Formule12) 25 De façon analogue à la formule 2, la formule 13 correspond au débit massique dm(sens) 32 mesuré à partir du second détecteur d'oxydes d'azote NOx à la somme du débit massique d'oxydes d'azote dm(NOx)_32 et du débit massique d'ammoniac dm(NH3) 32 30 entre les deux catalyseurs SCR : dm(sens) 32 = dm(NOx)_32 + dm(NH3) 32 (Formule 13) A partir des formules 11 et 12, et avec la formule 14 on 35 peut calculer le débit massique d'ammoniac dm(NH3) 32 : 3047272 11 (Formule14) dm(NH3)_32 = dm(sens) 32 x dm(NOx)_33 1 - ETA(23) On utilise pour cela le rendement ETA(23) d'un modèle du second catalyseur SCR 23. Le débit massique d'oxydes d'azote 5 dm(N0x) 33 en aval du second catalyseur SCR 23 peut être obtenu à partir du signal du troisième détecteur d'oxydes d'azote NOx 33 s'il n'y a pas de glissement d'ammoniac au niveau du second catalyseur SCR 23. 10 3047272 12 10 NOMENCLATURE DES ELEMENTS PRINCIPAUX 11 Moteur à combustion interne Conduite des gaz d'échappement 5 20 Système de catalyseurs 21 Unité de dosage d'agent réducteur 22, 23, 24 Catalyseurs SCR 31, 32, 33 Détecteurs d'oxydes d'azote NOx Appareil de commande électronique 10 REVENDICATIONS1°) Procédé de détermination du débit massique d'ammoniac (dm(NH2)_32) entre deux catalyseurs SCR (22, 23) installés l'un derrière l'autre dans une conduite de gaz d'échappement (11) faisant partie d'un système de catalyseurs SCR (20) qui ne comporte qu'une unité de do- sage d'agent réducteur (21) en amont du premier catalyseur SCR (22), procédé caractérisé en ce qu' on détermine le débit massique d'ammoniac à partir du signal d'un dé- tecteur d'oxydes d'azote (22) installé entre les deux catalyseurs SCR (22, 23) et le signal d'un détecteur d'oxydes d'azote (33) installé en aval du second catalyseur SCR (23). 2°) Procédé selon la revendication 1, caractérisé en ce que pour la détermination, on effectue la différence du signal du détecteur d'oxydes d'azote NOx (32) installé entre deux catalyseurs SCR et du signal du détecteur d'oxydes d'azote NOx (33) installé en aval du second catalyseur SCR (23). 3°) Procédé selon la revendication 1, caractérisé en ce qu' à partir du signal du détecteur d'oxydes d'azote NOx (32) installé entre les deux catalyseurs SCR (22, 23) et du signal du détecteur d'oxydes d'azote NOx (33) installé en aval du second catalyseur SCR, on calcule le besoin spécifique d'agent de dosage (facDos_20) du système de cata- lyseur SCR(20) et le besoin spécifique d'agent de dosage (facDos_22) du premier catalyseur SCR et on détermine le débit massique d'ammoniac (dm(NH3)_32) en comparant l'intégrale dans le temps des deux besoins d'agent de dosage (facDos_20, facDos_22). 4°) Procédé selon la revendication 3, caractérisé en ce que pour le calcul des deux besoins d'agent de dosage (facDos_20, facDos_22) on tient compte de la masse d'ammoniac (m(NH3)_21) qui a été dégagée de l'agent réducteur et qui a été dosée par l'unité de dosage 3047272 14 d'agent réducteur (21) dans la conduite des gaz d'échappement et de la masse (m(N0x)_31) des oxydes d'azote en amont du premier catalyseur SCR (22). 5 5°) Procédé selon l'une des revendications 1 à 4, caractérisé en ce qu' on fait la détermination en un point de fonctionnement du système de catalyseur SCR (20) où la variation des niveaux de remplissage d'ammoniac des deux catalyseurs SCR (22, 23) ne dépasse pas un seuil 10 de variation. 6°) Procédé selon l'une des revendications 1 à 5, caractérisé en ce qu' on détermine en un point de fonctionnement du système de catalyseur 15 SCR (20) où le glissement d'ammoniac sur le second catalyseur SCR (23) ne dépasse pas un seuil de glissement. 7°) Procédé selon la revendication 1, caractérisé en ce que 20 pour déterminer le débit massique d'ammoniac (dm(NH2)_32) on tient compte du rendement (ETA(23)) du second catalyseur SCR (23). 8°) Procédé selon l'une des revendications 1 à 7, caractérisé en ce que 25 l'on corrige les niveaux de remplissage d'ammoniac, modélisés des deux catalyseurs SCR (22, 23) à l'aide du débit massique d'ammoniac (dm(NH2)_32) obtenu. 9°) Programme d'ordinateur conçu pour exécuter les étapes du procédé 30 selon l'une des revendications 1 à 8 et support de mémoire lisible par une machine et contenant le programme d'ordinateur. 10°) Appareil de commande électronique (40) conçu pour déterminer à l'aide du procédé selon l'une des revendications 1 à 8, un débit mas- 35 sique d'ammoniac (dm(NH2)_32) entre deux catalyseurs SCR (22, 23) 3047272 15 installés l'un derrière l'autre dans la conduite des gaz d'échappement (11) dans un système de catalyseurs SCR (20) et qui ne comporte qu'une unité de dosage d'agent réducteur (21) en amont du premier catalyseur SCR (22). 5