FR 3085423 A1 20200306 FR 1857768 A 20180829 PROCÉDÉ D'ESTIMATION DE CHARGE D'UN FILTRE À PARTICULES L'invention a trait au domaine de la dépollution des moteurs à combustion, plus particulièrement à un procédé d'estimation de la charge en suies d'un filtre à particules ainsi qu'un véhicule équipé d'une ligne d'échappement comprenant un filtre à particules et d'un calculateur apte à mettre en oeuvre le procédé. Le document de brevet publié DE 10 2017 100 500 Al divulgue un procédé de contrôle de la régénération d'un filtre à particules comprenant plusieurs zones de forme annulaire dans lesquelles la charge en suies et la température sont déterminées séparément. Le déclenchement de la régénération intervient lorsque l'anneau en périphérie atteint une masse critique. Le document a pour objet d'éviter une combustion incomplète des suies dans la périphérie du filtre. Le système connu a le défaut de ne pas prendre en compte l'ensemble du filtre à particules. Le document de brevet publié FR 2989422 divulgue un procédé d'estimation de la charge en suies d'un filtre à particules basé sur des mesures de pression de gaz d'échappement en entrée et sortie d'un filtre à particules ainsi que sur une mesure d'un débit des gaz d'échappement. Ce document ne prend pas en compte le problème d'inhomogénéité de la charge en suies, ce qui réduit la plausibilité de l'estimation de la charge en suies. L'invention a pour objectif de pallier au moins un des inconvénients de l'état de la technique susmentionné. Plus particulièrement, l'invention a pour objectif d'améliorer l'estimation de la charge en suies d'un filtre à particules. L'invention a pour objet un procédé pour estimer la charge en suies d'un filtre à particules dans un conduit de gaz d'échappement d'un moteur à combustion interne, comprenant les étapes suivantes : modéliser le filtre à particules en au moins deux zones, remarquable par déterminer pour chaque zone du filtre à particules, un paramètre représentatif d'une densité d'une charge en suies dans la zone en question et sélectionner un mode de détermination d'une masse de suies 2 accumulées dans le filtre à particules en dépendance d'au moins deux paramètres représentatifs de la densité de la charge en suies déterminés à l'étape précédente. Avantageusement, le mode de détermination de la masse de suies accumulées dans le filtre à particules est au moins une détermination de la masse de suies dans 5 le filtre à particules basée sur au moins une différence de pression entre l'entrée et la sortie du filtre à particules ou une détermination de la masse de suies dans le filtre à particules basée sur la somme des masses déterminées pour chaque zone. Avantageusement, le paramètre représentatif de la densité de la charge en suies pour une zone donnée est basé sur une masse de suies accumulées déterminée 10 pour la zone en question. Avantageusement, la sélection du mode de détermination de la masse de suies accumulées dans le filtre à particules est fonction d'une comparaison entre les au moins deux paramètres représentatifs de la densité de la charge en suies. Avantageusement, la détermination de la masse de suies dans le filtre à particules 15 est basée sur au moins une différence de pression entre l'entrée et la sortie du filtre à particules lorsque au moins deux valeurs correspondant aux au moins deux paramètres représentatifs de la densité de la charge en suies sont égales ou la valeur absolue d'un écart entre deux des au moins deux valeurs est inférieure ou égale à un premier seuil. 20 Avantageusement, la détermination de la masse de suies dans le filtre à particules est basée sur la somme des masses déterminées pour chaque zone du filtre à particules, lorsque les au moins deux valeurs correspondant aux au moins deux paramètres représentatifs de la densité de la charge en suies sont différentes ou la valeur absolue de l'écart entre deux des au moins deux valeurs est supérieure au 25 premier seuil. Avantageusement, chaque zone du filtre à particules est modélisée avec un modèle OD et au moins une des variables modélisées pour une des au moins deux zones est(sont) au moins une température dans la zone en question, une exotherme de 3 combustion dans la zone en question, un débit de suies brûlées dans la zone en question et/ou une masse de suies accumulées dans la zone en question. Avantageusement, au moins une des entrées du modèle OD pour la zone en question est(sont) au moins une pression partielle d'Oxygène en amont du filtre à 5 particules, un débit de gaz d'échappement , un débit de suies, et/ou une température des gaz d'échappement entrant dans la zone en question. Avantageusement, une pluralité de cylindres de délimitation des au moins deux zones est prévue, chacune des au moins deux zones étant délimitée en partie par au moins un cylindre parmi les cylindres de délimitation, les cylindres de délimitation 10 étant coaxiaux ou parallèles entre eux. L'invention concerne aussi un véhicule comportant un moteur à combustion interne à injection directe et/ou allumage commandé, une ligne d'échappement comprenant un filtre à particules et un calculateur adapté à mettre en oeuvre le procédé. Préférentiellement, le paramètre représentatif de la densité de la charge en suies 15 pour une zone donnée est le rapport entre la masse de suies accumulées dans la zone rapportée au volume effectif de la zone. L'invention peut également avoir pour objet un procédé pour estimer la charge en suies d'un filtre à particules dans un conduit de gaz d'échappement d'un moteur à combustion interne, comprenant les étapes suivantes : 20 - modéliser le filtre à particules en au moins deux zones ; - déterminer pour au moins une de ces zones du filtre à particules une masse de suies ; remarquable en ce que la au moins une de ces zones du filtre à particules est modélisée avec un modèle OD 25 et au moins une des variables modélisées pour la zone en question est(sont) au moins une température dans la zone en question, une exotherme de combustion dans la zone en question, un débit de suies brûlées dans la zone en question et/ou une masse de suies accumulées dans la zone en question. Les mesures de l'invention sont intéressantes en ce qu'elles permettent d'éviter de sous-estimer la quantité de suies dans le filtre à particules et éviter une libération 4 excessive d'énergie en phase de régénération, évitant une usure thermique prématurée du filtre à particules. Elles permettent aussi d'établir des données susceptibles de recalibrer le modèle d'estimation de charge en suies. D'autres caractéristiques et avantages de la présente invention seront mieux 5 compris à l'aide de la description et des dessins parmi lesquels : - La figure 1 illustre un logigramme du procédé selon l'invention ; - La figure 2 représente un filtre à particules avec deux zones modélisées ; - La figure 3 illustre deux modèles OD correspondant aux deux zones avec les grandeurs d'entrée et les variables modélisées associées ; 10 - La figure 4 représente deux configurations de filtres à particules comprenant deux zones. La figure 1 décrit les étapes qui permettent d'estimer la charge en suies d'un filtre à particules 10 (acronyme FAP) dans un conduit de gaz d'échappement d'un moteur à combustion interne selon l'invention. Des masses de suies accumulées msuiesl, 15 msuies2 sont déterminées dans deux zones 11, 12 du filtre à particules à d'une l'aide - d modélisation numérique, par exemple. La figure 2 représente les deux zones 11, 12 dont l'une occupe une position centrale tandis que l'autre est en périphérie du filtre à particules 10. Un paramètre dl , d2 représentatif d'une densité d'une charge en suies dans une zone donnée 11, 12 tel qu'une densité dl, d2 est ensuite évalué pour 20 chaque zone 11, 12. La densité de charge en suies pour une zone donnée 11, 12 peut être définie comme étant le rapport entre la masse de suies accumulée m -suies 1, msuies 2 dans la zone 11, 12 rapportée au volume effectif de la zone V1, V2. Le volume effectif V1, V2 d'une zone donnée 11, 12 peut être défini par le volume de la zone 11, 12 du filtre à particules 10 ou par le volume des canaux de la zone 11, 12 25 en question. Egalement, si les zones 11, 12 ont le même volume effectif V1,V2, le paramètre dl, d2 représentatif de la densité de la charge en suies peut être relié directement à la masse msuies 1, msuies 2 de suies accumulées dans la zone 11, 12 en question. 5 La connaissance des valeurs des densités dl , d2 permet d'évaluer l'homogénéité de la charge en suies dans le filtre à particules 10. Dans certaines conditions, il peut arriver que la périphérie du filtre à particules soit plus chargée en suies que la partie centrale, car la périphérie est relativement plus froide, subissant plus de pertes 5 thermiques que la partie centrale. En conséquence, la partie périphérique du filtre à particules 10 profite moins du phénomène d'autocombustion des suies, s'activant pour une température seuil. Un manque d'homogénéité de la charge en suies fausse l'estimation de la charge en suies. La présente invention propose d'améliorer la détermination de la charge en suies même lorsque le filtre à particules 10 n'est pas 10 chargé en suies de manière homogène. A cette fin, l'invention sélectionne une méthode d'estimation de charge en suies adaptée à l'état d'homogénéité du filtre à particules 10. Le mode de détermination de la masse m de suies accumulées dans le filtre à particules est établi en dépendance d'une comparaison entre les deux densités dl, d2 estimées à l'étape précédente. Si 15 les valeurs correspondant aux deux densités dl d2 sont égales ou la valeur absolue de leur écart est inférieure ou égale à un premier seul, il est établi que le chargement en suies du filtre est homogène, et une masse m de suies accumulée dans le filtre à particules peut être évaluée en fonction d'une différence de pression entre l'entrée et la sortie du filtre à particules 10 et optionnellement du débit volumique traversant le 20 filtre à particules 10. Ce mode de détermination de la charge en suies est avantageux, car il découle d'une mesure directe. En effet, la masse m de suies accumulée dans le filtre à particules basée sur des mesures peut servir dans la prise de décision, par exemple pour le déclenchement d'une régénération du filtre à particules 10. Cette masse m peut aussi servir à recalibrer les modèles utilisés pour 25 déterminer les masses msuiesl, msuiesz. Il se peut qu'un modèle OD du filtre à particules (complet) 10 soit aussi recalibré. Si les valeurs correspondant aux deux densités dl, d2 sont différentes ou la valeur absolue de leur écart est supérieur au premier seul, il est établi que le chargement en suies dans le filtre à particules est inhomogène. Le masse m de suies 30 accumulées dans le filtre à particules basée sur une différence de pression entre l'entrée et la sortie du filtre à particules 10 et optionnellement sur le débit volumique traversant le filtre à particules n'est pas plausible, car les gaz d'échappement 6 passent préférentiellement dans la partie centrale du filtre à particule 10 moins chargée en suies que la périphérie. Étant donné que la correspondance entre la différence de pression et la masse de suies est calibrée pour une masse homogène, le mode de détermination de la masse m de suies basée sur une différence de 5 pression sous-estime la masse réelle de suies accumulées dans le filtre à particules 10. La masse m de suies accumulées déterminée par une différence de pression ne peut donc pas servir dans la prise de décision, par exemple pour le déclenchement d'une régénération du filtre à particules 10. Dans ce cas, la masse m de suies accumulées basée sur la somme des masses msuiesl, msuies2 est préférée à une 10 valeur déduite d'une mesure (par exemple différence de pression), et ce jusqu'à ce qu'il soit estimé que le chargement en suies est redevenu homogène. En plus des deux modes décrits précédemment, il existe d'autres modes de détermination d'une masse m de suies accumulées dans un filtre à particules tels qu'une méthode de radiofréquence, un modèle OD du filtre à particules complet, ou 15 un modèle du filtre découpé en plus de deux zones. L'invention ne doit pas être limitée aux deux modes présentés en figure 1, c.-à-d. la somme des masses ml, m2 ou une masse de suies accumulée dans le filtre à particules basé sur au moins une différence de pression entre l'entrée et la sortie du filtre à particules. Le filtre à particules 10 pour les moteurs à allumage commandé, notamment le 20 moteur essence, est dénommé GPF (Acronyme anglo-saxon pour désigner Gasoline Particulate Filter). Dans la suite, afin de faciliter la lecture, il sera fait référence au filtre à particules ou son acronyme FAP, tout en considérant qu'il englobe aussi le concept GPF. Comme illustré en figure 3, le filtre à particules 10 peut être décomposé en deux 25 modèles 21, 22 (sous-modèles). Chaque modèle 21, 22 correspond à une zone 11, 12 du filtre à particules et peut être décrit par modèle du type OD. Un modèle OD est basé sur une ou des équation(s) différentielle(s) dépendant du temps uniquement, contrairement aux modèles 1 D, 2D et 3D qui intègrent une ou des composante(s) spatiales. Le modèle OD selon un mode de réalisation de l'invention est basé sur un 30 bilan de la quantité de suies (conservation de la masse) et un bilan de puissance thermique (conservation d'énergie) dans chaque zone 11, 12 comme illustrée ci- 7 dessous. Le modèle OD pour chaque zone du filtre à particules 10 peut suivre l'évolution de 4 variables : une température dans la zone en question TFAP1, TFAP2, une exotherme de combustion dans la zone en question ExoTi, ExoT2, un débit de suies brûlées dans ladite zone en question Qmsuiescombi, Qmsuiescomb2 et une 5 masse de suies dans la zone en question msuiesl, msuies2. Les grandeurs d'entrée du modèle OD pour chaque zone du filtre à particules sont une pression partielle d'Oxygène 02 en amont du filtre à particules, un débit de gaz d'échappement Qmgart Qmgaz2, un débit des suies Qmsuiest Qmsuies2 et une température des gaz d'échappement Tgazi, Tgaz2 entrant dans la zone en question. 10 Pour chaque zone 11, 12, les variables et les grandeurs d'entrée sont indicées i: 1, 2. La première équation du modèle se rapporte à la conservation de la masse. mairies (g) = Q,n GeCaornüC g s)dt avec o le débit de suies entrant dans la zone i = 1, 2 ; 15 - ninâe'sQrnb le débit de suies brûlées dans la zone i, défini par une loi de type Arrhenius : 13E Qnsuis = KG KITF AF tri Éss É avec la constante cinétique de la réaction suies + 02 4 CO2 ; 20 l'énergie d'activation de cette réaction ; O TFAP la température dans la zone i la masse de suies dans la zone i ; o Pa la pression partielle en 02. 25 Un bilan de puissance thermique pour déterminer une température de la zone i moyenne est donné par l'équation suivante : 8 .0 Ti .0 PPM yt-1 Qmgaz L Pgar At FAP Qmgetz - FAP ifitFAp At avec o Q=7= gai ale débit d'échappement entrant dans la zone i ; 5 O gré la température du débit d'échappement entrant dans la zone i ; o 611-72 et CI PFAF les capacités thermiques de la zone i ; la masse de la zone i ; if la contribution des échanges thermiques : 10 o hi] le coefficient d'échange conducto-convectif de la surface j de la zone i ; osi./ la surface d'échange j de la zone i ; o rij la température externe au niveau de surface d'échange j de la zone i ; 15 o ExoT l'exotherme de combustion des suies déduit de la vitesse de combustion définie précédemment et de l'enthalpie de la réaction de la zone i : ExoT I = 0 Thpl i la température de la zone i calculée au pas de temps 20 précédent ; at le pas de temps du calcul. La figure 3 représente aussi un module 20 d'estimation permettant de déterminer les grandeurs d'entrée, comme le débit des gaz d'échappement Qmgaz-i, Qmgaz2 et de 25 suies Qmsffiesi, Qmsffies2, la température des gaz d'échappement entrant Tgazi, Tgaz2 pour chaque zone 11,12 du modèle 21, 22 à partir de valeurs mesurées et/ou modélisées, par exemple par cartographie. Le débit des gaz d'échappement Qmgari, 9 Qmgaz2, le débit de suies Qmsffiesi , Qmsffies2 et la température des gaz d'échappement Tgazi, Tgaz2 entrant dans chaque zone 11, 12 peuvent être modélisés par des cartographies respectives qui sont fonction d'au moins un des paramètres suivants : régime moteur, couple, un débit des gaz d'échappement en 5 amont du FAP Qmgaz et/ou une température des gaz d'échappement en amont du FAP Tgaz. Alternativement, la température des gaz d'échappement Tgarl entrant dans la zone centrale 11 peut être estimée à partir d'une cartographie d'offset (terme anglo-saxon signifiant décalage) de température qui peut dépendre du débit des gaz d'échappement amont du FAP Qmgaz et/ou température des gaz d'échappement en 10 amont du FAP Tgaz, et/ou de la vitesse du véhicule. Étant donné que la zone centrale 11 du filtre à particules 10 n'est pas en contact direct avec l'air ambiant, la conduction thermique entre la zone centrale 11 et la zone périphérique 12 peut être considérée comme négligeable. La contribution aux échanges thermiques des autres surfaces telles que les surfaces frontales de la 15 zone centrale 11, est aussi négligeable. En conséquence, la zone centrale 11 peut être modélisée comme étant adiabatique. Ce qui se traduit par une surface nulle d'échange stj = o pour chaque surface j de la zone centrale 11. Les paramètres liés aux échanges thermiques de la zone centrale sont repris dans les équations suivantes : h. n, = 0 h =0 20 Pour la zone périphérique 12, seuls les échanges de la paroi externe sont pris en considération dans les équations suivantes : = h 5 T =hSTanih o h le coefficient d'échange conducto-convectif entre le FAP (c.-à-d. la 25 zone périphérique 12 et l'air ambiant) ; o 5 la surface d'échange entre le FAP et l'air ambiant , o T',', la température ambiante. 10 La figure 4 représente deux illustrations schématiques du filtre à particules 10 modélisé avec les deux zones 11, 12. Comme représentées en figure 4, les deux zones 11, 12 sont délimitées par une surface cylindrique commune. Les deux zones 11, 12 peuvent être coaxiales ou parallèles. Le modèle ne doit pas être limité à deux 5 zones 11, 12, il est possible d'en REVENDICATIONS1. Procédé pour estimer la charge en suies d'un filtre à particules (10) dans un conduit de gaz d'échappement d'un moteur à combustion interne, comprenant les étapes suivantes : - modéliser le filtre à particules en au moins deux zones (11, 12) ; caractérisé par - déterminer pour chaque zone (11, 12) dudit filtre à particules (10), un paramètre (dl, d2) représentatif d'une densité d'une charge en suies dans ladite zone (11, 12) ; - sélectionner un mode de détermination d'une masse (m) de suies accumulées dans le filtre à particules (10) en dépendance d'au moins deux desdits paramètres (dl, d2) déterminés à l'étape précédente. 2. Procédé selon la revendication 1, caractérisé en ce que le mode de détermination de la masse (m) de suies accumulées dans le filtre à particules (10) est au moins une détermination de la masse (m) de suies dans le filtre à particules basée sur au moins une différence de pression entre l'entrée et la sortie du filtre à particules ou une détermination de la masse (m) de suies dans le filtre à particules basée sur la somme des masses (meuieel, nuisez) déterminées pour chaque zone (11, 12). 3. Procédé selon l'une des revendications précédentes, caractérisé en ce que le paramètre (dl, d2) représentatif de la densité de la charge en suies pour une zone donnée du filtre à particules (10) est basé sur une masse de suies accumulées (msuiesi, msuies2) déterminée pour ladite zone (11 ,12). 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que la sélection du mode de détermination de la masse (m) de suies accumulées dans le filtre à particules (10) est fonction d'une comparaison entre lesdits au moins deux paramètres (d1, d2). 5. Procédé selon la revendication 4, caractérisé en ce que la détermination de la masse (m) de suies dans le filtre à particules (10) est basée sur une 12 différence de pression entre l'entrée et la sortie du filtre à particules lorsque au moins deux valeurs (dl, d2) correspondant auxdits au moins deux paramètres (dl, d2) sont égales ou la valeur absolue d'un écart entre deux desdites au moins deux valeurs est inférieure ou égale à un premier seuil. 5 6. Procédé selon l'une des revendications 4 et 5, caractérisé en ce que la détermination de la masse (m) de suies dans le filtre à particules (10) est basée sur la somme des masses (mstuesl, msuies2) déterminées pour chaque zone du filtre à particules (10), lorsque lesdites au moins deux valeurs (dl, d2) correspondant auxdits au moins deux paramètres (dl, d2) sont différentes 10 ou la valeur absolue de l'écart entre deux desdites au moins deux valeurs est supérieure au premier seuil. 7. Procédé selon l'une des revendications précédentes, caractérisé en ce que chaque zone (11, 12) du filtre à particules (10) est modélisée avec un modèle OD (21, 22) et au moins une des variables modélisées pour une desdites au 15 moins deux zones (11, 12) est(sont) au moins une température (TFAP1, TFAP2) dans ladite zone (11, 12), une exotherme de combustion (ExoTt, ExoT2) dans ladite zone (10, 12), un débit de suies brûlées (Qmsuiescombi, Qmsuiescomb2) dans ladite zone (10, 12) et/ou une masse de suies (msuiest, msuies2) accumulées dans ladite zone (10, 12). 20 8. Procédé selon la revendication 7, caractérisé en ce qu' au moins une des entrées dudit modèle OD (21, 22) pour ladite zone est(sont) au moins une pression partielle d'Oxygène (02) en amont du filtre à particules (10), un débit de gaz d'échappement (Qmgazt , Qmgazz), un débit de suies (Qmsuiest, Qmsuies2), et/ou une température des gaz d'échappement (Tgazi, Tgaz2) entrant 25 dans ladite zone (10, 12). 9. Procédé selon l'une des revendications précédentes, caractérisé par une pluralité de cylindres de délimitation desdites au moins deux zones (11, 12), chacune desdites au moins deux zones étant délimitée en partie par au moins un cylindre parmi lesdits cylindres de délimitation, lesdits cylindres de délimitation étant coaxiaux ou parallèles entre eux. 13 10.Un véhicule comportant un moteur à combustion interne à injection directe et/ou allumage commandé, une ligne d'échappement comprenant un filtre à particules (10) et un calculateur adapté à mettre en oeuvre le procédé selon l'une des revendications 1 à 9. 5