专利摘要:
Méthode pour évaluer un niveau de bouchage d'un filtre à air (20) d'une unité HVAC (5), comprenant les étapes suivantes :- acquisition d'une première pluralité de paramètres (ΔP) représentatifs d'une perte de charge entre deux côtés (28A, 28B) du filtre à air à une pluralité d'instants,- calcul d'une seconde pluralité de paramètres, chacun étant représentatif d'une moyenne d'au moins certains des paramètres sur un segment de temps passé,- obtention d'une variété de scénarios représentatifs de conditions de fonctionnement futures de l'unité HVAC et ayant respectivement une probabilité p<sub>i</sub>,- obtention d'un paramètre représentatif d'une perte de charge maximale autorisée entre les deux côtés du filtre à air,- pour chacun des scénarios, calcul d'une troisième pluralité de paramètres représentatifs d'une évolution future de la perte de charge sur des segments de temps futurs, et- à l'aide de la troisième pluralité de paramètres, calcul du niveau de bouchage du filtre à air sous forme, soit d'une durée de vie utile résiduelle probable du filtre à air, soit d'une probabilité pour que, au bout d'un laps de temps donné, la perte de charge ait atteint la perte de charge maximale autorisée.Installation correspondante.
公开号:EP3680734A1
申请号:EP20151004.7
申请日:2020-01-09
公开日:2020-07-15
发明作者:Andrea Staino;Rami ABOU-EID
申请人:Alstom Transport Technologies SAS;
IPC主号:F24F11-00
专利说明:
[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] Le 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 indications 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 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.
[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 : acquisition d'une première pluralité de paramètres ΔPraw représentatifs d'une perte de charge entre deux côtés du filtre à air à une pluralité d'instants ; calcul d'une seconde pluralité de paramètres ΔPseg à partir de la première pluralité de paramètres ΔPraw, chacun des paramètres ΔPseg étant respectivement représentatif d'une moyenne d'au moins certains des paramètres ΔPraw sur un segment de temps passé ; obtention d'une variété de scénarios Si représentatifs de conditions de fonctionnement futures de l'unité HVAC sur une pluralité de segments de temps futurs, chacun des scénarios ayant respectivement une probabilité pi ; obtention d'un paramètre ΔPmax représentatif d'une perte de charge maximale autorisée entre les deux côtés du filtre à air ; à l'aide d'un modèle, pour chacun des scénarios Si, calcul d'une troisième pluralité de paramètres (ΔP'seg)i 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 à l'aide de la troisième pluralité de paramètres (ΔP'seg)i, calcul du niveau de bouchage du filtre à air sous forme : 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 (ΔP'seg)i auront atteint le paramètre ΔPmax, soit d'une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax.
[0007] 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 : la méthode comporte en outre les étapes suivantes : 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 ΔPseg 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 ΔPseg divisé par le carré d'un débit d'air Qair dans le filtre à air, le débit d'air Qair étant un débit d'air capable de créer dans le filtre à air une perte de charge égale audit paramètre ΔPseg ; et à 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 ; 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 (ΔP'seg)i aient atteint le paramètre ΔPmax 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 D0 ; 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 ΔPseg, certains des paramètres ΔPraw sont sélectionnés respectivement pour obtenir chacun des paramètres ΔPseg, les paramètres ΔPraw sélectionnés étant tels que la tension électrique U mesurée au même instant était supérieure à un seuil de tension Umax ; la méthode comporte en outre une acquisition d'un angle d'ouverture θ à ladite pluralité d'instants, l'unité HVAC comprenant une enceinte, et une prise d'air ayant un clapet définissant l'angle d'ouverture θ par rapport à l'enceinte, et, à l'étape de calcul de la seconde pluralité de paramètres ΔPseg, certains des paramètres ΔPraw sont sélectionnés respectivement pour obtenir chacun des paramètres ΔPseg, les paramètres ΔPraw sélectionnés étant tels que l'angle d'ouverture θ est sensiblement égal à une même valeur prédéterminée, de préférence nulle ; chacun des segments de temps passés et chacun des segments de temps futurs possède une durée égale à 1 jour ; l'étape de calcul d'une troisième pluralité de paramètres (ΔP'seg)i comporte un calage du modèle sur les segments de temps passé, le calage utilisant les paramètres ΔPseg ; 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 (ΔP'seg)i auront atteint le paramètre ΔPmax dans l'un des scénarios Si, ladite moyenne étant pondérée par les probabilités pi ; et l'obtention du paramètre ΔPmax comprend les sous-étapes suivantes : obtention d'un paramètre Qair,min à partir d'un débit nominal d'air de l'unité HVAC, le paramètre Qair,min étant représentatif d'un débit minimal d'air conditionné souhaité pour l'unité HVAC, et calcul du paramètre ΔPmax en utilisant le paramètre Qair,min ; mesure de la première pluralité de paramètres ΔPraw ; et mise à disposition du niveau de bouchage calculé, par affichage ou transmission du niveau de bouchage.
[0008] L'invention a également pour objet un procédé de maintenance de l'unité HVAC comprenant les étapes suivantes : mise en œuvre d'une méthode telle que décrite ci-dessus, en fonction du niveau de bouchage, réalisation d'un acte de maintenance tel qu'un remplacement du filtre à air.
[0009] 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 : un module d'acquisition adapté pour acquérir une première pluralité de paramètres ΔPraw représentatifs d'une perte de charge entre deux côtés du filtre à air à une pluralité d'instants, une mémoire dans laquelle sont enregistrés : une variété de scénarios Si représentatifs de conditions de fonctionnement futures de l'unité HVAC sur une pluralité de segments de temps futurs, et un paramètre ΔPmax représentatif d'une perte de charge maximale autorisée entre les deux côtés du filtre à air, un module de calcul adapté pour calculer : une seconde pluralité de paramètres ΔPseg à partir de la première pluralité de paramètres ΔPraw, chacun des paramètres ΔPseg étant respectivement représentatif d'une moyenne d'au moins certains des paramètres ΔPraw sur un segment de temps passé, pour chacun des scénarios Si, une troisième pluralité de paramètres (ΔP'seg)i 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, à l'aide de la troisième pluralité de paramètres (ΔP'seg)i, le niveau de bouchage du filtre à air sous forme : 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 (ΔP'seg)i auront atteint le paramètre ΔPmax, soit d'une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax.
[0010] 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é.
[0011] L'invention a également pour objet une installation comportant : une unité HVAC, et un calculateur tel que décrit ci-dessus.
[0012] 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 : [Fig 1] la figure 1 est une vue schématique d'une installation selon l'invention, [Fig 2] la figure 2 est un graphique représentant des pertes 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, [Fig 3] la figure 3 est un graphique illustrant une phase de calage d'une méthode selon l'invention, [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 à air, et [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 à air ait atteint une valeur prédéterminée.
[0013] En référence à la figure 1, on décrit une installation 1 selon l'invention.
[0014] L'installation 1 comprend une unité HVAC 5, et un calculateur 10.
[0015] 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é).
[0016] 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.
[0017] L'unité HVAC 5 est conçue pour fournir un débit d'air conditionné Qair é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é Qair que l'unité HVAC 5 est capable de fournir diminue.
[0018] 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.
[0019] 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 ΔP, c'est-à-dire une différence de pression entre le côté amont et le côté aval. La perte de charge ΔP 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.
[0020] Le filtre à air 20 est équipé d'un capteur de pression 30 adapté pour mesurer la perte de charge ΔP à une pluralité d'instants, par exemple tous les 1 seconde lorsque l'unité HVAC 5 est en service.
[0021] La prise d'air 22 comporte un clapet 32 définissant un angle d'ouverture θ par rapport à l'enceinte 24 de l'unité HVAC 5, et avantageusement un capteur d'ouverture 33 adapté pour mesurer l'angle d'ouverture.
[0022] L'angle d'ouverture θ conditionne la quantité d'air extérieur admise dans l'unité HVAC 5.
[0023] 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é).
[0024] En variante (non représentée), le calculateur 10 est embarqué à bord du véhicule ferroviaire 15.
[0025] 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.
[0026] Le module de mise à disposition 40 est par exemple un écran.
[0027] Le module d'acquisition 34 est adapté pour acquérir une première pluralité de paramètres ΔPraw 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 θ à la même pluralité d'instants.
[0028] Par exemple, les paramètres ΔPraw sont des pertes de charges ΔP 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.
[0029] 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 ΔPmax représentatif d'une perte de charge maximale autorisée entre les deux côtés 28A, 28B du filtre à air 20.
[0030] 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.
[0031] Par définition, la somme des 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).
[0032] Le paramètre ΔPmax est par exemple une perte de charge maximale, avantageusement calculée par le calculateur 38 à partir d'un paramètre Qair,min représentatif d'un débit minimal d'air conditionné souhaité pour l'unité HVAC 5.
[0033] En variante, le paramètre ΔPmax est directement fourni et inscrit dans la mémoire 36.
[0034] Le paramètre Qair,min 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.
[0035] En variante, le paramètre Qair,min est directement fourni et inscrit dans la mémoire 36.
[0036] Le module de calcul 38 est adapté pour calculer une seconde pluralité de paramètres ΔPseg à partir de la première pluralité de paramètres ΔPraw, et, pour chacun des scénarios Si, une troisième pluralité de paramètres (ΔP'seg)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. Le module de calcul 38 est en outre adapté pour calculer, à l'aide de la troisième pluralité de paramètres (ΔP'seg)i et des probabilités pi, le niveau de bouchage du filtre à air 20.
[0037] 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.
[0038] Chacun des paramètres ΔPseg est respectivement représentatif d'une moyenne d'au moins certains des paramètres ΔPraw sur un segment de temps passé.
[0039] Par exemple, les segments de temps passés et les segments de temps futurs ont des durées égales à 1 jour.
[0040] Le niveau de bouchage du filtre à air 20 est par exemple : 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 (ΔP'seg)i auront atteint le paramètre ΔPmax, soit une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax.
[0041] 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.
[0042] 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 ΔP dans le filtre à air 20 en pascals (Pa).
[0043] Le module d'acquisition 34 réalise tout d'abord l'acquisition de la première pluralité de paramètres ΔPraw, 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 θ 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.
[0044] Sur la figure 2, une barre verticale T1 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.
[0045] Puis, le module de calcul 38 calcule la seconde pluralité de paramètres ΔPseg à partir de la première pluralité de paramètres ΔPraw. Les paramètres ΔPseg forment une courbe C1 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.
[0046] Chaque paramètre ΔPseg est une moyenne de certains des paramètres ΔPraw mesurés pendant l'un des segments de temps passés, c'est-à-dire une moyenne journalière éventuellement conditionnelle.
[0047] Avantageusement, les paramètres ΔPraw sélectionnés par le module de calcul 38 pour obtenir chacun des paramètres ΔPseg sont tels que la tension électrique U mesurée au même instant était supérieure à un seuil de tension Umax, par exemple 470 V. Dit autrement, on ne tient pas compte des paramètres ΔPraw obtenus lorsque la tension électrique U est trop basse.
[0048] 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 tension électrique U. Il a été observé que le fait de ne pas prendre en compte dans les moyennes ΔPseg les paramètres ΔPraw 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 ΔPseg.
[0049] Avantageusement, les paramètres ΔPraw sélectionnés pour obtenir chacun des paramètres ΔPseg sont tels que l'angle d'ouverture θ 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 ΔPraw obtenus lorsque l'angle d'ouverture θ n'est pas sensiblement égal à 0°, c'est-à-dire que le clapet n'est pas fermé.
[0050] Par « sensiblement égal », on entend en pratique, par exemple, que la différence est inférieure à 5°, par exemple inférieure à 1°.
[0051] 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 θ est alors non nul. Ensuite, le clapet 32 est refermé et l'angle d'ouverture θ redevient sensiblement nul. Or, la valeur de l'angle d'ouverture θ 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 ΔPraw mesurés lorsque l'angle d'ouverture θ 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 ΔPseg.
[0052] Ainsi, ΔPseg = f1(ΔPraw,U,θ) oùF1 est une moyenne conditionnelle.
[0053] Selon une variante, tous les paramètres ΔPraw sont utilisés pour calculer les moyennes ΔPseg. Dans ce cas, Pseg = f1(ΔPraw) et la moyenne f1 n'est pas conditionnelle.
[0054] 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 ΔPseg correspondant audit segment de temps passé.
[0055] La résistance hydraulique est définie comme le ratio dudit paramètre ΔPseg divisé par le carré d'un débit d'air Qair capable de créer dans le filtre à air 20 une perte de charge égale audit paramètre ΔPseg :K = ΔP seg /Q air 2
[0056] 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 à air 20 et la résistance hydraulique K0 d'un filtre à air neuf (non représenté). Le paramètre D est un indicateur de santé du filtre à air 20.Ainsi, D = f3(K, K0)Si K0 est un nombre, alors, par exemple, D = K - K0Si K0 est exprimée sous la forme d'une distribution de moyenne µ et d'écart type σ, la fonction f3 est par exemple la distance de Mahalanobis :D = K − σ / µ
[0057] La perte de charge maximale ΔPmax est par exemple donnée.
[0058] En variante, elle est calculée à partir du paramètre Qair,min représentatif d'un débit minimal d'air conditionné (ou de froid) souhaité pour l'unité HVAC 5 : ΔP max = f 4 Q air , min= f 2 − 1 Q air , min
[0059] Dans l'exemple représenté sur la figure 2, ΔPmax vaut environ 130 Pa.
[0060] Chacun des scénarios Si 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).
[0061] 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 (ΔP'seg)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.
[0062] Sur la figure 2, on a symbolisé les paramètres (ΔP'seg)i par des courbes C2, C3, C4.
[0063] 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.
[0064] Les courbes C3 et C4 encadrent la courbe C2 et représentent respectivement et symboliquement 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.
[0065] Avantageusement, le module de calcul 38 réalise un calage de certains de ses paramè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 ΔPseg réellement obtenus à partir des mesures sur cette période.
[0066] Avantageusement, le module de calcul 38 ne calcule les paramètres (ΔP'seg)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é.
[0067] Les scénarios Si fournissent par principe tous les paramètres dont le modèle utilisé par le module de calcul 38 a besoin.
[0068] Pour chacune des pluralités de paramètres (ΔP'seg)i, 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 (ΔP'seg)i atteignent la perte de charge maximale ΔPmax.
[0069] Sur la figure 4 est représenté un histogramme de fréquence des durées de vie utiles résiduelles RUL 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.
[0070] 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 (ΔP'seg)i auront atteint le paramètre ΔPmax. 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 ΔPseg dans le filtre à air 20 aura atteint la perte de charge maximale ΔPmax souhaitée.
[0071] 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.
[0072] En variante, à partir des paramètres (ΔP'seg)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.
[0073] La figure 5 est un histogramme de fréquence de ces pertes de charge. La barre verticale C8 symbolise le paramètre ΔPmax. Le module de calcul 38 estime la probabilité pour que, au bout du laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax. Dans l'exemple, cette probabilité est d'environ 45%.
[0074] 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.
[0075] Optionnellement, en fonction du niveau de bouchage, un remplacement du filtre à air 20 est effectué.
[0076] 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, 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.
[0077] 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).
权利要求:
Claims (10)
[0001] Méthode pour évaluer un niveau de bouchage d'un filtre à air (20) d'une unité HVAC (5), la méthode comprenant les étapes suivantes :
- acquisition d'une première pluralité de paramètres ΔPraw représentatifs d'une perte de charge entre deux côtés (28A, 28B) du filtre à air (20) à une pluralité d'instants,
- calcul d'une seconde pluralité de paramètres ΔPseg à partir de la première pluralité de paramètres ΔPraw, chacun des paramètres ΔPseg étant respectivement représentatif d'une moyenne d'au moins certains des paramètres ΔPraw sur un segment de temps passé,
- obtention d'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, chacun des scénarios ayant respectivement une probabilité pi,
- obtention d'un paramètre ΔPmax représentatif d'une perte de charge maximale autorisée entre les deux côtés (28A, 28B) du filtre à air (20),
- à l'aide d'un modèle, pour chacun des scénarios Si, calcul d'une troisième pluralité de paramètres (ΔP'seg)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, et
- à l'aide de la troisième pluralité de paramètres (ΔP'seg)i, calcul du niveau de bouchage du filtre à air (20) sous forme :
- soit d'une durée de vie utile résiduelle probable du filtre à air (20), représentative d'un laps de temps probable au bout duquel les paramètres (ΔP'seg)i auront atteint le paramètre ΔPmax,
- soit d'une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax.
[0002] Méthode selon la revendication 1, comportant en outre les étapes suivantes :
- 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 ΔPseg correspondant audit segment de temps passé, le paramètre K étant représentatif d'une résistance hydraulique du filtre à air (20), la résistance hydraulique étant définie comme le ratio dudit paramètre ΔPseg divisé par le carré d'un débit d'air Qair dans le filtre à air (20), le débit d'air Qair étant un débit d'air capable de créer dans le filtre à air (20) une perte de charge égale audit paramètre ΔPseg, et
- à 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 (20) et un filtre à air (20) neuf.
[0003] Méthode selon la revendication 2, dans laquelle l'étape de calcul de la durée de vie utile résiduelle probable du filtre à air (20) ou de la probabilité pour que, au bout d'un laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax est réalisée uniquement si le paramètre D représentatif d'une différence de résistance hydraulique entre le filtre à air (20) et un filtre à air (20) neuf a dépassé un seuil prédéterminée D0.
[0004] Méthode selon l'une quelconque des revendications 1 à 4, comportant en outre une acquisition d'une tension électrique U d'alimentation de l'unité HVAC (5) à ladite pluralité d'instants, etdans laquelle, à l'étape de calcul de la seconde pluralité de paramètres ΔPseg, certains des paramètres ΔPraw sont sélectionnés respectivement pour obtenir chacun des paramètres ΔPseg, les paramètres ΔPraw sélectionnés étant tels que la tension électrique U mesurée au même instant était supérieure à un seuil de tension Umax.
[0005] Méthode selon l'une quelconque des revendications 1 à 4, comportant en outre une acquisition d'un angle d'ouverture θ à ladite pluralité d'instants, l'unité HVAC (5) comprenant une enceinte (24), et une prise d'air (22) ayant un clapet (32) définissant l'angle d'ouverture θ par rapport à l'enceinte (24), etdans laquelle, à l'étape de calcul de la seconde pluralité de paramètres ΔPseg, certains des paramètres ΔPraw sont sélectionnés respectivement pour obtenir chacun des paramètres ΔPseg, les paramètres ΔPraw sélectionnés étant tels que l'angle d'ouverture θ est sensiblement égal à une même valeur prédéterminée, de préférence nulle.
[0006] Méthode selon l'une quelconque des revendications 1 à 5, dans laquelle chacun des segments de temps passés et chacun des segments de temps futurs possède une durée égale à 1 jour.
[0007] Méthode selon l'une quelconque des revendications 1 à 6, dans laquelle l'étape de calcul d'une troisième pluralité de paramètres (ΔP'seg)i comporte un calage du modèle sur les segments de temps passé, le calage utilisant les paramètres ΔPseg.
[0008] Méthode selon l'une quelconque des revendications 1 à 7, dans laquelle 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 (ΔP'seg)i auront atteint le paramètre ΔPmax dans l'un des scénarios Si, ladite moyenne étant pondérée par les probabilités pi.
[0009] Méthode selon l'une quelconque des revendications 1 à 8, dans laquelle l'obtention du paramètre ΔPmax comprend les sous-étapes suivantes :
- obtention d'un paramètre Qair,min à partir d'un débit nominal d'air de l'unité HVAC (5), le paramètre Qair,min étant représentatif d'un débit minimal d'air conditionné souhaité pour l'unité HVAC (5), et
- calcul du paramètre ΔPmax en utilisant le paramètre Qair,min.
[0010] Calculateur (10) configuré pour évaluer un niveau de bouchage d'un filtre à air (20) d'une unité HVAC (5), le calculateur (10) comportant :
- un module d'acquisition (34) adapté pour acquérir une première pluralité de paramètres ΔPraw représentatifs d'une perte de charge entre deux côtés (28A, 28B) du filtre à air (20) à une pluralité d'instants,
- une mémoire (36) dans laquelle 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, et
- un paramètre ΔPmax représentatif d'une perte de charge maximale autorisée entre les deux côtés (28A, 28B) du filtre à air (20),
- un module de calcul adapté pour calculer :
- une seconde pluralité de paramètres ΔPseg à partir de la première pluralité de paramètres ΔPraw, chacun des paramètres ΔPseg étant respectivement représentatif d'une moyenne d'au moins certains des paramètres ΔPraw sur un segment de temps passé,
- pour chacun des scénarios Si, une troisième pluralité de paramètres (ΔP'seg)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,
- à l'aide de la troisième pluralité de paramètres (ΔP'seg)i, le niveau de bouchage du filtre à air (20) sous forme :
- soit d'une durée de vie utile résiduelle probable du filtre à air (20), représentative d'un laps de temps probable au bout duquel les paramètres (ΔP'seg)i auront atteint le paramètre ΔPmax,
- soit d'une probabilité pour que, au bout d'un laps de temps donné H, les paramètres (ΔP'seg)i aient atteint le paramètre ΔPmax.
类似技术:
公开号 | 公开日 | 专利标题
US9911147B2|2018-03-06|Systems and methods for implementing automated intelligence-based bidding for repair services for environmental control systems in monitored buildings
US10202980B2|2019-02-12|Centrifugal compressor apparatus and method for preventing surge therein
EP2209976B1|2017-10-11|Dispositif de surveillance d&#39;une pompe haute-pression dans un circuit d&#39;alimentation en carburant d&#39;une turbomachine
US8150720B2|2012-04-03|Dispatch management model
US10675579B2|2020-06-09|Method and system for predicting residual useful life of an air filter
EP2568295B1|2014-01-15|Procédé et dispositif d&#39;estimation automatique d&#39;une vitesse air d&#39;un aéronef
CN102862685B|2015-06-03|检测飞机部件的性能进入衰退期的方法及飞机的维修方法
US6644095B2|2003-11-11|Continuous on-board diagnostic lubricant monitoring system and method
EP1604170B1|2015-05-06|Procede et systeme pour surveiller le comportement d&#39;une tuyauterie contenant un fluide sous pression
US20050075769A1|2005-04-07|Aircraft accessory monitor
EP1456045B1|2007-05-09|Installation de climatisation de vehicule munie d&#39;un dispositif electronique de controle
US4852693A|1989-08-01|Lubrication monitoring apparatus for machine
EP2078995B1|2016-08-17|Procédé de surveillance d&#39;un moteur à turbine à gaz
EP1689986B1|2007-03-21|Systeme d&#39;aide a la maintenance d&#39;un filtre a particules integre dans une ligne d&#39;echappement d&#39;un moteur de vehicule automobile
US10371406B2|2019-08-06|Maintenance and diagnostics for refrigeration systems
CN107436014A|2017-12-05|用于hvac安装和诊断的移动应用的系统及方法
KR20000011630A|2000-02-25|스팀트랩을검사하는방법및장치,및스팀트랩을관리하기위한관리시스템
US8131400B2|2012-03-06|Adaptive on-tool mass flow controller tuning
EP2488833B1|2020-08-05|Dispositif de détection du blocage d&#39;un compteur mécanique de fluide, et compteur mécanique de fluide avec détection de blocage
EP1831663B1|2008-09-10|Dispositif de mesure a resonateur et procede mettant en oeuvre le dispositif
US10775085B2|2020-09-15|Energy management for refrigeration systems
EP1784561A1|2007-05-16|Système d&#39;evaluation de l&#39;état de charge en suies de moyens de dépollution
US20040165185A1|2004-08-26|Fluid particle monitor and methods related thereto
AU2014206171B2|2017-11-02|System and method for monitoring lubricant of an engine
JP6244276B2|2017-12-06|航空機翼防氷バルブの監視システム及び方法
同族专利:
公开号 | 公开日
ZA202000136B|2021-01-27|
FR3091576B1|2021-03-12|
FR3091576A1|2020-07-10|
引用文献:
公开号 | 申请日 | 公开日 | 申请人 | 专利标题
法律状态:
2020-06-12| PUAI| Public reference made under article 153(3) epc to a published international application that has entered the european phase|Free format text: ORIGINAL CODE: 0009012 |
2020-06-12| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED |
2020-06-26| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
2020-07-15| AX| Request for extension of the european patent|Extension state: BA ME |
2020-07-15| AK| Designated contracting states|Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
2020-07-29| 17P| Request for examination filed|Effective date: 20200624 |
2020-07-29| RIN1| Information on inventor provided before grant (corrected)|Inventor name: ABOU-EID, RAMI Inventor name: STAINO, ANDREA |
2020-07-29| RBV| Designated contracting states (corrected)|Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
2021-10-25| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: GRANT OF PATENT IS INTENDED |
2021-10-25| GRAP| Despatch of communication of intention to grant a patent|Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
2021-10-28| GRAJ| Information related to disapproval of communication of intention to grant by the applicant or resumption of examination proceedings by the epo deleted|Free format text: ORIGINAL CODE: EPIDOSDIGR1 |
2021-10-28| GRAP| Despatch of communication of intention to grant a patent|Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
2021-10-28| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
2021-11-03| RIC1| Information provided on ipc code assigned before grant|Ipc: G05B 23/02 20060101ALI20210930BHEP Ipc: G05B 17/02 20060101AFI20210930BHEP |
2021-11-08| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: GRANT OF PATENT IS INTENDED |
2021-11-08| GRAP| Despatch of communication of intention to grant a patent|Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
2021-11-24| INTG| Intention to grant announced|Effective date: 20211026 |
2021-12-01| INTC| Intention to grant announced (deleted)|
2021-12-08| INTG| Intention to grant announced|Effective date: 20211109 |
优先权:
申请号 | 申请日 | 专利标题
[返回顶部]