Brevet FR3091576A1 — publié le 20200710

Description Titre de l'invention : Méthode pour évaluer un niveau de bouchage d'un filtre à air dans une unité HVAC [0001] La présente invention concerne une méthode pour évaluer un niveau de bouchage d'un filtre à air dans une unité HVAC (abréviation de l'anglais Heating, Ventilation and Air-Conditioning, c'est-à-dire chauffage, ventilation et climatisation), en particulier une unité embarquée dans un véhicule ferroviaire. [0002] L'invention concerne également un calculateur correspondant, et une installation comportant l'unité HVAC et un tel calculateur. [0003] Lc filtre à air d'une unité HVAC est obstrué progressivement par l'accumulation de particules telles que la poussière, la saleté, l'herbe, le pollen, etc. Un filtre à air bouché a tendance à prolonger les plages de fonctionnement du compresseur, réduire le rendement énergétique de l'unité, et réduire la capacité de refroidissement et de chauffage. Ceci peut conduire à l'apparition de givre sur les serpentins d'évaporateur et éventuellement causer une défaillance de l'unité. Par conséquent, le filtre doit être changé périodiquement afin d'éviter ces problèmes. [0004] Généralement, les filtres à air sont remplacés périodiquement en fonction des in- dications du fabricant. Le remplacement est souvent effectué sans évaluer la capacité de rétention résiduelle réelle du filtre. De plus, l'inspection des filtres à air pour un véhicule ferroviaire nécessite habituellement l'arrêt du véhicule ferroviaire dans un dépôt et l'accès au toit. Cela a une incidence sur la disponibilité des véhicules ferroviaires. Un remplacement inutile du filtre est coûteux et prend du temps. Il n'est donc pas souhaitable. [0005] Un but de l'invention est de permettre unc estimation précise du niveau de bouchage du filtre à air, par exemple dans le but d'éviter un remplacement trop précoce du filtre à air. [0006] A cet effet, l'invention concerne une méthode pour évaluer un niveau de bouchage d'un filtre à air d'une unité HVAC, la méthode comprenant les étapes suivantes : [0007] - acquisition d'une première pluralité de paramètres AP,,,' représentatifs d'une perte de charge entre deux côtés du filtre à air à une pluralité d'instants ; [0008] - calcul d'une seconde pluralité de paramètres AP,, à partir de la première pluralité de paramètres chacun des paramètres AP,, étant respectivement représentatif d'une moyenne d'au moins certains des paramètres AP,' sur un segment de temps passé ; [0009] - obtention d'une variété de scénarios S; représentatifs de conditions de fonc- tionnement futures de l'unité HVAC sur une pluralité de segments de temps futurs, 2 chacun des scénarios ayant respectivement une probabilité pi ; [0010] - obtention d'un paramètre AP,'', représentatif d'une perte de charge maximale autorisée entre les deux côtés du filtre à air ; [0011] - à l'aide d'un modèle, pour chacun des scénarios Si, calcul d'une troisième pluralité de paramètres (AP''di représentatifs d'une évolution future de la perte de charge entre les deux côtés du filtre à air sur au moins certains des segments de temps futurs ; et [0012] - à l'aide de la troisième pluralité de paramètres (AP's',)i, calcul du niveau de bouchage du filtre à air sous forme : [0013] - soit d'une durée de vie utile résiduelle probable du filtre à air, représentative d'un laps de temps probable au bout duquel les paramètres (AP's')i auront atteint le paramètre AP'',, [0014] - soit d'une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (AP''di aient atteint le paramètre AP'm'.. [0015] Selon des modes particuliers de réalisation, la méthode comprend l'une ou plusieurs des caractéristiques suivantes, prise(s) seule ou selon toutes les combinaisons techniquement possibles : [0016] - la méthode comporte en outre les étapes suivantes : [0017] - pour l'un des segments de temps passés, de préférence le dernier temporellement, calcul d'un paramètre K en utilisant le paramètre AP,' correspondant audit segment de temps passé, le paramètre K étant représentatif d'une résistance hydraulique du filtre à air, la résistance hydraulique étant définie comme le ratio dudit paramètre AP,' divisé par le carré d'un débit d'air (lir dans le filtre à air, le débit d'air Q';, étant un débit d'air capable de créer dans le filtre à air une perte de charge égale audit paramètre AP,'; et [0018] - à partir du paramètre K, calcul d'un paramètre D représentatif d'une différence de résistance hydraulique entre le filtre à air et un filtre à air neuf ; [0019] - l'étape de calcul de la durée de vie utile résiduelle probable du filtre à air ou de la probabilité pour que, au bout d'un laps de temps donné H, les paramètres (AP's'); aient atteint le paramètre AP,'', est réalisée uniquement si le paramètre D représentatif d'une différence de résistance hydraulique entre le filtre à air et un filtre à air neuf a dépassé un seuil prédéterminée Do; [0020] - la méthode comporte en outre une acquisition d'une tension électrique U d'alimentation de l'unité HVAC à ladite pluralité d'instants, et, à l'étape de calcul de la seconde pluralité de paramètres AP,', certains des paramètres APraw sont sélectionnés respectivement pour obtenir chacun des paramètres AP,,, les paramètres APraw sélectionnés étant tels que la tension électrique U mesurée au même instant était supérieure à un seuil de tension Umax ; [0021] - la méthode comporte en outre une acquisition d'un angle d'ouverture 0 à ladite pluralité d'instants, l'unité HVAC comprenant une enceinte, et une prise d'air ayant un 3 clapet définissant l'angle d'ouverture 0 par rapport à l'enceinte, et, à l'étape de calcul de la seconde pluralité de paramètres AP',', certains des paramètres AP,,' sont sélectionnés respectivement pour obtenir chacun des paramètres AP',,, les paramètres AP 'an, sélectionnés étant tels que l'angle d'ouverture 0 est sensiblement égal à une même valeur prédéterminée, de préférence nulle ; [0022] - chacun des segments de temps passés et chacun des segments de temps futurs possède une durée égale à 1 jour ; [0023] - l'étape de calcul d'une troisième pluralité de paramètres comporte un calage du modèle sur les segments de temps passé, le calage utilisant les paramètres AP 'g; [0024] - la durée de vie résiduelle probable est calculée comme une moyenne de durées de vie résiduelles, chaque durée de vie résiduelle étant respectivement représentative d'un laps de temps au bout duquel les paramètres (AP'''), auront atteint le paramètre AP''' dans l'un des scénarios Si, ladite moyenne étant pondérée par les probabilités p, ; et [0025] - l'obtention du paramètre AP,''' comprend les sous-étapes suivantes : [0026] - obtention d'un paramètre Q'i,'',,, à partir d'un débit nominal d'air de l'unité HVAC, le paramètre Qd,-.,'in étant représentatif d'un débit minimal d'air conditionné souhaité pour l'unité HVAC, et [0027] - calcul du paramètre AP',,,' en utilisant le paramètre ; [0028] - mesure de la première pluralité de paramètres AP,',''; et [0029] - mise à disposition du niveau de bouchage calculé, par affichage ou transmission du niveau de bouchage. [0030] L'invention a également pour objet un procédé de maintenance de l'unité HVAC comprenant les étapes suivantes : [0031] - mise en oeuvre d'une méthode telle que décrite ci-dessus, [0032] - en fonction du niveau de bouchage, réalisation d'un acte de maintenance tel qu'un remplacement du filtre à air. [0033] L'invention a également pour objet un calculateur configuré pour évaluer un niveau de bouchage d'un filtre à air dans une unité HVAC, le calculateur comportant : [0034] - un module d'acquisition adapté pour acquérir une première pluralité de paramètres AP,'', représentatifs d'une perte de charge entre deux côtés du filtre à air à une pluralité d'instants, [0035] - une mémoire dans laquelle sont enregistrés [0036] - une variété de scénarios S, représentatifs de conditions de fonctionnement futures de l'unité HVAC sur une pluralité de segments de temps futurs, et [0037] - un paramètre AP',' représentatif d'une perte de charge maximale autorisée entre les deux côtés du filtre à air, [0038] - un module de calcul adapté pour calculer : [0039] - une seconde pluralité de paramètres AP,' à partir de la première pluralité de pa- ramètres AP',', chacun des paramètres AP', étant respectivement représentatif d'une moyenne d'au moins certains des paramètres AP,'', sur un segment de temps passé, [0040] - pour chacun des scénarios Si, une troisième pluralité de paramètres (AP'',), repré- sentatifs d'une évolution future de la perte de charge entre les deux côtés du filtre à air sur au moins certains des segments de temps futurs, [0041] - à l'aide de la troisième pluralité de paramètres (AP''), le niveau de bouchage du filtre à air sous forme : [0042] - soit d'une durée de vie utile résiduelle probable du filtre à air, représentative d'un laps de temps probable au bout duquel les paramètres (AP''), auront atteint le paramètre AP'', [0043] - soit d'une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (AP'seg)i aient atteint le paramètre AP',,,. [0044] Selon un mode de réalisation particulier, le calculateur comprend un module de mise à disposition du niveau de bouchage adapté pour afficher ou transmettre le niveau de bouchage calculé. [0045] L'invention a également pour objet une installation comportant : [0046] - une unité HVAC, et [0047] - un calculateur tel que décrit ci-dessus. [0048] L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés, sur lesquels : [0049] [fig.1] la figure 1 est une vue schématique d'une installation selon l'invention, [0050] [fig.2] la figure 2 est un graphique représentant des pelles de charge dans le filtre à air de l'installation représentée sur la figure 1, à différentes phases d'une méthode selon l'invention, [0051] [fig.3] la figure 3 est un graphique illustrant une phase de calage d'une méthode selon l'invention, [0052] [fig.4] la figure 4 est un histogramme illustrant une phase de calcul d'une méthode selon l'invention pour obtenir la durée de vie utile résiduelle probable du filtre à ah-, et [0053] [fig.5] la figure 5 est un histogramme illustrant une phase de calcul d'une méthode selon l'invention pour obtenir une probabilité pour que, au bout d'un laps de temps donné H de deux semaines, la perte de charge dans le filtre à ah- ait atteint une valeur prédéterminée. [0054] En référence à la figure 1, on décrit une installation 1 selon l'invention. [0055] L'installation 1 comprend une unité HVAC 5, et un calculateur 10. [0056] L'unité HVAC 5 est par exemple embarquée à bord d'un véhicule ferroviaire 15, tandis que le calculateur 10 est distant, par exemple situé dans un centre de maintenance (non représenté). [0057] L'unité HVAC 5 comprend un filtre à air 20, une prise d'air 22, et une enceinte 24 adaptée pour abriter les autres éléments classiques d'une unité HVAC qui ne seront pas décrits ici. [0058] L'unité HVAC 5 est conçue pour fournir un débit d'air conditionné Qui, égal à une valeur nominale lorsque l'unité HVAC 5 fonctionne à 100% de sa capacité nominale. Au fur et à mesure que le filtre à air 20 se bouche, le débit d'air conditionné Qu que l'unité HVAC 5 est capable de fournir diminue. [0059] L'unité HVAC 5 est par exemple connectée électriquement à un convertisseur statique 26 situé dans le véhicule ferroviaire 15 et adapté pour délivrer à l'unité HVAC une tension électrique U ajustable. [0060] Le filtre à air 20 possède deux côtés 28A, 28B entre lesquels le passage d'un flux d'air (non représenté) crée une perte de charge AP, c'est-à-dire une différence de pression entre le côté amont et le côté aval. La perte de charge AP ainsi définie est le plus souvent positive, mais peut parfois s'inverser et devenir négative selon les conditions de fonctionnement de l'unité HVAC 5. [0061] Le filtre à ah- 20 est équipé d'un capteur de pression 30 adapté pour mesurer la perte de charge AP à une pluralité d'instants, par exemple tous les 1 seconde lorsque l'unité HVAC 5 est en service. [0062] La prise d'air 22 comporte un clapet 32 définissant un angle d'ouverture 0 par rapport à l'enceinte 24 de l'unité HVAC 5, et avantageusement un capteur d'ouverture 33 adapté pour mesurer l'angle d'ouverture. [0063] L'angle d'ouverture 0 conditionne la quantité d'air extérieur admise dans l'unité HVAC 5. [0064] Le calculateur 10 est configuré pour évaluer un niveau de bouchage du filtre à air 20. Le calculateur 10 est par exemple distant, c'est-à-dire non embarqué dans le véhicule ferroviaire 15. Le calculateur 10 est avantageusement hébergé dans un centre de maintenance (non représenté). [0065] En variante (non représentée), le calculateur 10 est embarqué à bord du véhicule fer- roviaire 15. [0066] Le calculateur 10 comprend un module d'acquisition 34, une mémoire 36, et un module de calcul 38. Avantageusement, le calculateur comprend aussi un module de mise à disposition 40 pour afficher ou transmettre le niveau de bouchage. [0067] Le module de mise à disposition 40 est par exemple un écran. 6 [0068] Le module d'acquisition 34 est adapté pour acquérir une première pluralité de pa- ramètres AP,,' représentatifs de la perte de charge entre les deux côtés 28A, 28B du filtre à air 20 à une pluralité d'instants, par exemple toutes les secondes lorsque l'unité HVAC 5 est en service. Avantageusement, le module d'acquisition 34 est aussi adapté pour acquérir la tension électrique U et/ou l'angle d'ouverture 0 à la même pluralité d'instants. [0069] Par exemple, les paramètres AP'' sont des pertes de charges AP entre les deux côtés 28A, 28B, avantageusement mesurées par le capteur de pression 30 ou obtenues à partir de mesures faites par le capteur de pression. [0070] Dans la mémoire 36 sont enregistrés une variété de scénarios Si représentatifs de conditions de fonctionnement futures de l'unité HVAC 5 sur une pluralité de segments de temps futurs, des probabilités pi que les scénarios Si se produisent effectivement, et un paramètre AP''', représentatif d'une perte de charge maximale autorisée entre les deux côtés 28A, 28B du filtre à air 20. [0071] Par exemple, chacun des scénarios Si définit la distribution de diamètre de particule, la concentration de particule et temps de fonctionnement du système. [0072] Par définition, la somme dcs probabilités pi est égale à 1. Selon un cas particulier, les scénarios sont générés de manière à être équiprobables (méthode de Monte Carlo). [0073] Le paramètre AP,''' est par exemple une perte de charge maximale, avantageusement calculée par le calculateur 38 à partir d'un paramètre Qi',' représentatif d'un débit minimal d'air conditionné souhaité pour l'unité HVAC 5. [0074] En variante, le paramètre AP',' est directement fourni et inscrit dans la mémoire 36. [0075] Le paramètre est avantageusement calculé par le calculateur 38 comme étant une fraction, par exemple 90%, du débit nominal d'air conditionné de l'unité HVAC 5. [0076] En variante, le paramètre Q',', est directement fourni et inscrit dans la mémoire 36. [0077] Le module de calcul 38 est adapté pour calculer une seconde pluralité de paramètres AP', à partir de la première pluralité de paramètres AP',, et, pour chacun des scénarios Si, une troisième pluralité de paramètres (AP''di représentatifs d'une évolution future de la perte de charge entre les deux côtés 28A, 28B du filtre à air 20 sur au moins certains des segments de temps futurs. Le module de calcul 38 est en outre adapté pour calculer, à l'aide de la troisième pluralité de paramètres (AP''g)i et des probabilités pi, le niveau de bouchage du filtre à air 20. [0078] Avantageusement, le module de calcul 38 est également adapté pour calculer un paramètre K représentatif d'une résistance hydraulique du filtre à air 20. [0079] Chacun des paramètres APs, est respectivement représentatif d'une moyenne d'au moins certains des paramètres AP', sur un segment de temps passé. [0080] Par exemple, les segments de temps passés et les segments de temps futurs ont des durées égales à 1 jour. 7 [0081] Le niveau de bouchage du filtre à air 20 est par exemple : [0082] - soit une durée de vie utile résiduelle probable du filtre à air, représentative d'un laps de temps probable au bout duquel les paramètres (AP'',), auront atteint le paramètre AP,',,,, [0083] - soit une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (AP''), aient atteint le paramètre AP',,,. [0084] Le fonctionnement de l'installation 1 et en particulier celui du calculateur 10 vont maintenant être expliqués pour illustrer une méthode selon l'invention. [0085] Sur la figure 2 est représenté un diagramme avec, en abscisse, le temps sous forme de date (du 6 mars au 22 septembre dans l'exemple), et, en ordonnée, la perte de charge AP dans le filtre à air 20 en pascals (Pa). [0086] Le module d'acquisition 34 réalise tout d'abord l'acquisition de la première pluralité de paramètres et avantageusement, d'une part, de la tension électrique U fournie par le convertisseur statique 26 et, d'autre part, de l'angle d'ouverture 0 du clapet 32. Ces paramètres ne sont pas représentés sur la figure 2. L'acquisition est réalisée dans l'exemple à une pluralité d'instants s'étalant sur une période allant du 6 mars au 25 juillet. [0087] Sur la figure 2, une barre verticale Tl représente le moment où une opération de maintenance du filtre à air 20, par exemple son remplacement, était initialement prévue. Une barre verticale T2 représente le temps présent, ici le 26 juillet. Grâce à la méthode selon l'invention, le remplacement prévu a pu être évité, et n'a toujours pas été jugé nécessaire au temps présent correspondant à la barre verticale T2. [0088] Puis, le module de calcul 38 calcule la seconde pluralité de paramètres AP,, à partir de la première pluralité de paramètres AP','. Les paramètres AP,, forment une courbe Cl sur la figure 2, montrant une augmentation à peu près régulière de la perte de charge dans le filtre à air 20 jusqu'au temps présent. [0089] Chaque paramètre AP', est une moyenne de certains des paramètres AP,,,' mesurés pendant l'un des segments de temps passés, c'est-à-dire une moyenne journalière éventuellement conditionnelle. [0090] Avantageusement, les paramètres AP,' sélectionnés par le module de calcul 38 pour obtenir chacun des paramètres AP,, sont tels que la tension électrique U mesurée au même instant était supérieure à un seuil de tension Unt,' par exemple 470 V. Dit autrement, on ne tient pas compte des paramètres AP,'' obtenus lorsque la tension électrique U est trop basse. [0091] En effet, la valeur de la perte de charge à travers le filtre à air 20 change lorsque la fréquence des ventilateurs (non représentés) de l'unité HVAC 5 change. Si l'unité HVAC 5 est équipée d'un variateur de fréquence (non représenté), les ventilateurs sont contrôlés en ajustant leur fréquence de fonctionnement, qui est proportionnelle à la 8 tension électrique U. Il a été observé que le fait de ne pas prendre en compte dans les moyennes AP'.' les paramètres AP'' mesurés lorsque la tension électrique U est inférieure à un certain seuil augmente la précision de la méthode en améliorant la cohérence interne des paramètres AP,''. [0092] Avantageusement, les paramètres AP,,,' sélectionnés pour obtenir chacun des pa- ramètres APseg sont tels que l'angle d'ouverture 0 est sensiblement égal à une même valeur prédéterminée, de préférence la plus petite possible, de préférence nulle. Dit autrement, on ne tient pas compte des paramètres AP,'' obtenus lorsque l'angle d'ouverture 0 n'est pas sensiblement égal à 0e, c'est-à-dire que le clapet n'est pas fermé. [0093] Par « sensiblement égal », on entend en pratique, par exemple, que la différence est inférieure à 5', par exemple inférieure à 1'. [0094] Lorsque l'unité HVAC 5 a besoin d'air frais, le clapet 32 est ouvert pour permettre à de l'air extérieur à l'enceinte 15 de pénétrer dans l'unité HVAC. L'angle d'ouverture 0 est alors non nul. Ensuite, le clapet 32 est refermé et l'angle d'ouverture 0 redevient sensiblement nul. Or, la valeur de l'angle d'ouverture 0 est susceptible d'affecter la perte de charge mesurée dans le filtre à air 20. Il a été observé que ne sélectionner que les paramètres AP,', mesurés lorsque l'angle d'ouverture 0 est sensiblement égal à une même valeur prédéterminée augmente la précision de la méthode, en améliorant la cohérence interne des paramètres 4P,. [0095] Ainsi, AP', = fl(AP,',U,O) où [0096] F1 est une moyenne conditionnelle. [0097] Selon une variante, tous les paramètres APraw sont utilisés pour calculer les moyennes AP',. Dans ce cas, P,ae = fl(AP,'',) et la moyenne fl n'est pas conditionnelle. [0098] Avantageusement, pour l'un des segments de temps passés, de préférence le dernier temporellement, le module de calcul 38 calcule le paramètre K, dans l'exemple une résistance hydraulique du filtre à air 20, en utilisant le paramètre AP,, correspondant audit segment de temps passé. [0099] La résistance hydraulique est définie comme le ratio dudit paramètre AP,, divisé par le carré d'un débit d'air (2, capable de créer dans le filtre à air 20 une perte de charge égale audit paramètre AP', : [0100] K = AP', / Q,2 où [0101] Qin = f2 (AP,,,), [0102] f2 étant une fonction connue de l'homme du métier et issue des principes de fonc- tionnement des turbomachines et de la mécanique des fluides. La fonction f2 est par exemple polynômiale et est avantageusement fournie par le fabricant du filtre. [0103] A partir du paramètre K, le module calcule avantageusement d'un paramètre D repré- sentatif d'une différence entre la résistance hydraulique du filtre à ah- 20 et la ré- 9 sistance hydraulique KO d'un filtre à air neuf (non représenté). Le paramètre D est un indicateur de santé du filtre à air 20. [0104] Ainsi, D = f3(K, KO) [0105] Si KO est un nombre, alors, par exemple, D = K - KO [0106] Si KO est exprimée sous la forme d'une distribution de moyenne Ft et d'écart type o, la fonction f3 est par exemple la distance de Mahalanohis : [0107] D= 1K-G /Ft [0108] La perte de charge maximale AP',', est par exemple donnée. [0109] En variante, elle est calculée à partir du paramètre représentatif d'un débit minimal d'air conditionné (ou de froid) souhaité pour l'unité HVAC 5 : [0110] AP'',, = f4(Q,', ''') = f21 (Q','''') [0111] Dans l'exemple représenté sur la figure 2, AP',' vaut environ 130 Pa. [0112] Chacun des scénarios S, représente par exemple les conditions de fonctionnement de l'unité HVAC 5 sur la période du 26 juillet au 22 septembre (c'est-à-dire le futur, pour un temps présent arbitrairement fixé au 26 juillet dans l'exemple). [0113] A l'aide d'un modèle connu de l'homme du métier, tel que les modèles de colmatage de milieu poreux basés sur la loi de Darcy, le module de calcul 38 calcule, pour chacun des scénarios Si, une troisième pluralité de paramètres (AP',,g)i représentatifs d'une évolution future de la perte de charge entre les deux côtés 28A, 28B du filtre à air 20 sur au moins certains des segments de temps futurs, couvrant dans l'exemple la période du 26 juillet au 22 septembre. [0114] Sur la figure 2, on a symbolisé les paramètres (AP's')i par des courbes C2, C3, C4. [0115] La courbe C2 représente symboliquement l'évolution la plus probable de la perte de charge dans le filtre à air 20. La courbe C2 est obtenue en calculant, pour chacun des segments de temps futur (ici chaque jour), la perte de charge la plus probable. [0116] Les courbes C3 et C4 encadrent la courbe C2 et représentent respectivement et sym- boliquement les valeurs entre lesquelles la perte de charge a 90% de chances de se situer. Dit autrement, pour chaque segment de temps futur, la perte de charge a 95% de chances de se situer au-dessus de la courbe C3 et 95% de chances de se situer au-dessous de la courbe C4. [0117] Avantageusement, le module de calcul 38 réalise un calage de certains de ses pa- ramètres internes non fournis par les scénarios Si. Par exemple, comme représenté sur la figure 3, sur la période passée allant du 6 mars au 26 juillet, le module de calcul 38 utilise l'état du filtre à 20 au 6 mars et le scénario qui s'est réellement produit entre le 6 mars et le 26 juillet pour calculer, à l'aide du modèle, l'évolution de perte de charge journalière sur cette période représentée par une courbe C7. Puis, le module de calcul 38 ajuste lesdits paramètres à caler pour que la courbe C7 soit la plus proche possible d'une courbe C6 représentant les valeurs de AP,' réellement obtenus à partir des mesures sur cette période. [0118] Avantageusement, le module de calcul 38 ne calcule les paramètres (AP',')i que si la résistance hydraulique K a dépassé un seuil haut prédéterminé, ou si le paramètre D (indicateur de santé du filtre à air 20) a franchi un seuil bas prédéterminé. [0119] Les scénarios S, fournissent par principe tous les paramètres dont le modèle utilisé par le module de calcul 38 a besoin. [0120] Pour chacune des pluralités de paramètres (AP',')' c'est-à-dire chacune des évolutions futures de la perte de charge associée à un des scénarios Si, le module de calcul 38 calcule une durée de vie utile résiduelle RUL,, en jour, représentative du laps de temps au bout duquel les paramètres (AP'''), atteignent la perte de charge maximale [0121] Sur la figure 4 est représenté un histogramme de fréquence des durées de vie utiles résiduelles RULi pour tous les scénarios Si. Dans le cas particulier où les scénarios Si sont équiprobables, ces fréquences représentent directement des probabilités. [0122] Le module de calcul 38 calcule la durée de vie utile résiduelle la plus probable du filtre à air 20, représentative d'un laps de temps probable au bout duquel les paramètres (AP'',), auront atteint le paramètre AP,'''. Dans l'exemple, la durée de vie utile résiduelle la plus probable est d'environ 14 ou 15 jours. Dit autrement, il est probable que, 14 jours après le 26 juillet, la perte de charge AP,,,g dans le filtre à air 20 aura atteint la perte de charge maximale AP,''' souhaitée. [0123] La durée de vie utile résiduelle la plus probable est représentative du niveau de bouchage du filtre à air 20 en ce sens que, toutes choses égales par ailleurs, elle est d'autant plus courte que le filtre à air est déjà bouché au 26 juillet. [0124] En variante, à partir des paramètres (AP'',)i, le module de calcul 38 évalue, pour chacun des scénarios Si, la perte de charge au bout d'un laps de temps donné H, par exemple deux semaines. [0125] La figure 5 est un histogramme de fréquence de ces pertes de charge. La barre verticale C8 symbolise le paramètre AP,'''. Le module de calcul 38 estime la probabilité pour que, au bout du laps de temps donné H, les paramètres (AP''), aient atteint le paramètre AP,'',. Dans l'exemple, cette probabilité est d'environ 45%. [0126] Cette probabilité est représentative du niveau de bouchage du filtre à air 20 en ce sens que, toutes choses égales par ailleurs, elle est d'autant plus élevée que le filtre à air est déjà bouché au 26 juillet. [0127] Optionnellement, en fonction du niveau de bouchage, un remplacement du filtre à air 20 est effectué. [0128] Ainsi, grâce aux caractéristiques décrites ci-dessus, le calculateur et la méthode ci- dessus permettent une estimation précise du niveau de bouchage du filtre à air, par exemple dans le but d'éviter un remplacement trop précoce du filtre à air 20. En effet, 11 alors que la date normale programmée pour ce remplacement était vers le 25 mai, la méthode permet d'envisager un remplacement peu avant le 25 août. [0129] La méthode fournit une prédiction probabiliste de la durée de vie utile restante du filtre. L'état bouché du filtre est anticipé et la durée de vie résiduelle probable est avantageusement communiquée au mainteneur. Les inspections périodiques deviennent moins utiles voire inutiles. Les pièces de rechange requises peuvent être commandées à l'avance et la disponibilité des ressources (personnel, installations, etc.) peut être planifiée à l'avance. La détection des filtres et le système d'estimation de la durée de vie restante peuvent donc améliorer la disponibilité de la flotte et réduire les coûts de maintenance (le nombre de pièces de rechange utilisées est considérablement réduit). 12 [Revendication I] [Revendication 2]