FR 3011946 A1 20150417 FR 1359903 A 20131011 SLEWEILLANCE N'OTT:Un D ,f.:AONEF POUR ANTICIKR LES OPÉRATIONS DL MANTENANCE 1171.,.F VON La présente invention se rapporte au domaine de surveillance d'un moteur d'aéronef. En particulier, l'invention concerne un procédé et un système de surveillance d'un moteur d'aéronef pour pronostiquer les opérations de maintenance en estimant une marge résiduelle de température de sortie des gaz d'échappement du moteur. Durant chaque vol, un aéronef procède à l'enregistrement et l'envoi d'informations sur son fonctionnement selon un système de message nommé ACARS (Aircraft Communication Addressing and Reporting System). Ces données sont récupérées par les stations au sol en temps réel pour être traitées imméd'i,tement dans le cas d'anomalies évidentes et sinon pour être archivées avec toutes les données rle la flotte. Les données relatives au moteur sont ensuite consultées par des experts pour surveiller le bon fonctionnement du moteur. Entre autres, les experts analysent les données concernant la température des gaz d'échappement EGT (Exhaust Gas Temperature) issus des vols successifs pour c:iagno:»:quer l'état d'un moteur et 20 pronostiquer les ,2éiaillances. En particulier, les eyperts estiment la marge résidulle de température de sortie EGT qui représente la différence entre un seuil de température limite et la température EGT mesurée. Ainsi, il est important de sun'eil'i:,r l'évolution de cette marge de manière à pouvoir anticiper le franchissement d'un seuil au-delà duquel le moteur doit être c:-'posé pour restauration de ses performances. 25 On notera qu'en début de vie, la marge est élevée et diminue progressivement avec le nombre de cycles. Ainsi, le signal temporel de la marge résiduelle de température évolue normalement selon une courbe décroissante représentative c'R l'usure naturelle du moteur. Toutefois, au lieu de décroitre continuellement, le signal temporel de la 30 marge qu'on observe habituellement présente ies oscillations aléatoires. Ces oscillations sont principalement dues aux interventions effectuées sur le moteur et à moindre mesure aux conditions j'acuisition des mesures ck, t,;:mpérature. Des augmentations artificielles de la marge résiduelle de température peuvent être provoquées par des nettoyages plus ou moins fréquents des moteurs réalisés pour éliminer les saletés accumulées de vol en vol. Ces oscillations du signal temporel de la marge résiduelle de température compliquent l'analyse et ne permettent pas de prévoir convenablement à l'avance et avec précision la date de défaillance du moteur. L'objet de la présente invention est par conséquent, de remédier aux inconvénients précités en proposant un procédé et un système de surveillance de la marge résiduelle de température de sortie des gaz d'échappement d'un moteur d'aéronef permettant de pronostiquer de manière simple, précise et suf:hamment à l'avance la défaillance rr,,.teur. OFPFT ET RÉS N La préserV.~ invention concerne un procédé de surveillance d'un moteur d'aéronef, comportant les étapes suivantes : - acquisition d'un s nal temporel de marge résiduelle de température de sortie des gaz d'échappement dudit moteur d'aéronef, - lissage dudit signal temporel pour former une première courbe représentative de ladite marI,,c résiduelle de température, - identification des morceaux descendants dans ladite première courbe, - construction d'une deuxième courbe par con,:aténation desdits morceaux descendants, ladite deuxième courbe étant continue tout en étant restreinte auxdits morceaux descendants de ladite première courbe, et - construction d'un modèle de prédiction à partir de ladite deuxième courbe pour déterminer au moins un indicateur de pronostic de déf.:311.,nce. Ce procédé modélise de manière simple et rapide l'usure réelle du moteur en éliminant les augmentations artificielles de la marge de température et permettant ainsi de pronostiquer avec précision la défaillance du moteur. Avantageusement, l'identification des morceaux descendants dans ladite première courbe comporte les étapes suivantes : - application sur la première courbe d'un modèle statistis:vJe d'accroissement décomposable en deux parties imMpendantes composées d'une première fonction décroissante représentative d'une usure habituelle du moteur d'aéronef et d'une deuxième fonction échelon composée de sauts déclenchés de manière aléatoire et représentatifs d'interventions ponctuelles sur le moteur d'aéronef, - recherche de remontées correspondantes auxdits sauts, et - identification des morceaux descendants en supprimant les points desdites remontées sur la première courbe. Ceci permet d'éliminer les augmentations artificielles de la marge de température pour ne conserver que les intervalles correspondants aux morceaux décroissants représentatifs de l'usure du moteur. Avantageusement, la construction par concaténation de ladite deuxième courbe comporte un collage ees..iits morceaxe7,,endants en déplaçant chaque morceau précédent pour la faire joindre au morceau suivant en commençant par le dernier morceau et en remontant dans le temps de proche en proche. Ceci permet de modéliser la tendance décroissan:e L-1..jiquant l'usure réelle du moteur. Selon un premier mode de réalisation de l'invention, la construction dudit modèle de prédiction comporte les étapes suivantes : - construction d'un modèle autorégressif modélisant l'évolution de la marge résiduelle de température en utilisant l'historique ladite deuxième courbe, et - application d'un filtre dynamique sur ledit modèle autorégressif pour déterminer leiit au moins un indicateur de pronostic de défaillance. Cette construction d'un modèle ,-yrIér ligue permet de faire des pronostics précis sur tout type de moteur et à n'importe quelle période de la vie du moteur. Avantageusement, le filtre dynamique est sélectionné parmi l'ensemble des filtres particulaires suivants : filtre bayesien, filtre de Kalman, filtres de Kalman étendus. Selon un deuxième mode de réalisation de l'invention, la construction dudit modèle de prédiction comporte les étapes suivantes : - construction d'un modèle linéaire de l'évolution de la marge de température en utilisant l'historique de ladite deuxième courbe, et - application d'une technique de régression sur ledit modèle linéaire pour déterminer ledit au moins un indicateur de pronostic de défaillance. La construction d'un modèle linéaire est très simple et permet de faire des pronostics précis avec très peu de calcul. Avantageusement, ledit au moins un indicatrur pronostic de défaillance est sélectionné parmi un ensemble d'indicateurs comprenant : -un premier indicateur d'estimation de la probabilité de franchissement d'un seuil de défaillance avant un horizon temporel prédéterminé, et -un deuxième indicateur d'estimation de la date de franchissement d'un seuil de défaillance. Avantageusement, l'acquisition dudit signal temporel de marge résiduelle de température comporte les étapes suivantes : - acquisition au cours du temps des mesures de la température de sortie des gaz d'échappement du moteur d'aéronef, - normalisation desdites mesures de tempratf,re par rapport à une température de référence iso standard formant ainsi des mesures de température normalisées, - standardisation desdites nsures de température normalisées en prenant en compte des données de contexte formant ainsi des mesures de température standarc;sées, et -calcul des marges entre lesdites mesures de température standardisées et une valeur de température maximale pré.-14,-rminée (en fonction du moteur) pour former ledit signal temporel de marge. Ainsi, le signal temporel de marge est recueilli selon des conditions physiques standards et indépendamment du contexte. L'invention vise également un système de surveillance d'un moteur d'aéronef, comportant: - des moyens d'acquisition et de traitement configurés pour recueillir un signal temporel de marge résiduelle de température de sortie des gaz d'échappement dudit moteur d'aéronef, - des moyens d'acquisition et de traitement configurés pour lisser ledit signal 1,,mporel formant ainsi une première courbe représentative de ladite marge résiduelle Ce température, - des moyens d'acquisition et de traitement configurés pour identifier des morceaux descendants dans ladite première courbe, - des moyens d'acquisition et de traitement configurés pour construire une deuxième courbe par concaténation desdits liorceaux descendants, ladite deuxième courbe étant continue tout en étant restreinte ::,ixdits morceaux descendants de ladite première courbe, - des moyens d'acquisition et de traitement configurés pour construire un modèle de prédiction à partir 1-.1e I ditr deuxième ccrwb(?. pour déterminer au moins un indic de prcnostic de défaillance. L'invenion vise aussi un programme d'ordinateur susceptible d'être mis en oeuvre par des moyens cle traitement et comportant des instructions de code a- aptées à la mise en 7:uvre dl procédé de surveillance selon l'une quelconc,ue des caractéristiques ci-dessus. ,euge niFez S!NS D'autr,s particularités et avantages du système et du procédé selon l'invention ressortiront mieux à la lecture de L, description faite ci-après, à titre indicatif mais non limitatif, en référence aux dessins annexés sur lesqu ? - la Fig. 1 illustre de manière schématique des moyens matériels mis en oeuvre dans le système ou procédé qui Î-:euvent être employés pour la surveillance d'un moteur d'aéronef selon l'invention ; - la Fig. 2A illustre les mesures EGT d'un moteur d'aéronef pendant un ensemble de vols successifs ; - la F. 2B illustre les mesures EGT standardisées du moteur d'aéronef pendant un ensemble de vols successifs ; - la Fig. 3 illustre le signal temporel de marge 'ésiduelle de temnérature de sortie des gaz d'échappement ; - la Fig. 4 illustre une première courbe correspondant au lissage du signal temporel de la Fig. 3 sur une période prédéterminée, selon l'invention ; - la Fig. 5 illustre une représentation graphique des remontées correspondant aux sauts de la première courbe, selon l'invention ; - la Fig. 6 illustre une représentation graphique d'une deuxième courbe formée par la concaténation des morceaux descendants, selon l'invention ; - la Fig. 7 illustre l'application d'un filtre particuiJk sur le système dynamique de la deuxième courbe, selon l'invention ; - la Fig. 8 illustre l'application d'un filtre particulaire sur le système dynamique de la deuxième courbe relative à un moteur jeune, selon l'invention ; et - la Fig. 9 illustre de manière schématique deux indicateurs de pronostic de défaillance, selon l'invention. 15 DESCRIPTION i TAPI -F. F Le principe d:? c,)nsiste à supprimer des points correspondant à une augmentation artificielle c,c. la marge résiduelle de température pour ne conserver que les parties décroissantes représentatives de l'usure réelle du mote . f,;nsi, en 20 observant l'évolution de cette usure, il est possible de pronostiquer avec une grande précision la future défaillance du moteur et les opérations de maintenance à mener. La Fig. 1 illustre un exemple de moyen.; nq-1, riels mis en oeuvre dans le système ou procédé de surveillance d'un moteur d'aéronef selon l'invention. La présente invention peut également s'appliquer à un aéronef cornpor1ant plusieurs 25 moteurs, tel que représenté sur la Fig. 1. Au cours d'un vol, un aéron,4 1 procède à l'enregistrement d'informations sur son fonctionnement ainsi que sur di-iférents paramètres environnementaux. Ces données enregistrées par des calculateurs embarqués par l'aéronef (par exemple, FADEC, ACMS, etc.) sont issues des mesures fournies par des moyens de mesure ou 30 capteurs intégrés dans l'aéronef 1. Par exemple, le FADEC (qui contrôle le moteur 2 de l'aéronef 1) enregistre un certain nombre de données mesurées par des capteurs intégrés au moteur 2 permettant à la fois de contrôler le moteur 2 et de servir de base à une procédure prédictive de maintenance. Les calculateurs de l'aéronef 1 recueillent ainsi au cours du temps des données relatives au moteur 2 de l'aéronef et son environnement. A chaque acquisition, ces données comprennent des informations relatives à des paramètres endogènes décrivant le comportement du moteur 2 ainsi qu'à des paramètres exogènes décrivant le contexte d'acquisition. A titre d'exemple, les paramètres endogènes comprennent la température de sortie des gaz d'échappement EGT (Exhaust Gas Temperature), les vitesses de rotation des arbres, le dé_it 'lu carburant, les températures et pressions de fluides en différentes localisations Co moteur (par exemple, avant et/ou après compression), etc. Les paramètres exogènes peuvent comprendre la température extérieure, l'altitude, le poids de l'avion, la géométrie variable de vanne de décharge, les consignes des turbines haute pression et basse pression, vitesse de l'avion, etc. Par ailleurs, un aéronef 1 envoie régulièrement au sol de courts messa3es instantanés concernant les paramètres endogènes et ex: 'f. 2?.nes. Durant chaque v0., l'aéronef 1 envoie en général au moins deux messages au sol, l'un lors du décollage et l'autre pendant la phase de croisière. Ces messages sont notamment envoyés par satellite (protocole ACARS) grâce à un système de transmission do.? données numériques entre l'aéronef en vol et le sol (d'autres protocoles de communication sont possibles : PCMCIA, 3G, etc.). Les stations au sol 3 récupèrent les différents messages émis à différentes dates pour différents aéronefs 1 et pour différents moteurs 2 et les envoient ensuite par l'intermédiaire d'un moyen ue communication à un centre de gestion 5. Ce dernier comporte un système informatiue 7 comprenant de manière habituelle des moyens d'entrée 9, des moyens d'acquisition et de traitement 11, des moyens de stockage 13, et des moyens de sortie 15. On notera que d'autres données enregistrées durant le vol sur les calculateurs embarqués peuvent aussi être déchargées régulièrement pour enrichir le recueil des infbrmiJons relatives aux moteurs 2. Les différentus données issues des messages reçus directement des différents aéronefs 1 ou celles récupérées au sol depuis les mémoires internes des calculateurs embarqués, sont stockées dans les moyens de stockage 13 pour constituer une base de données 14 sur toute une flotte de moteurs 2. La présente invention s'intéresse surtout aux données relatives à la température de sortie des gaz d'échappement EGT des moteurs. La Fig. 2A illustre les mesures EGT d'un moteur d'aéronef pendant un ensemble de vols successifs. L'axe des ordonnées représente les valeurs EGT et l'axe des abscisses le nombre de cycles ou vols. Certains capteurs intégrés dans le moteur d'aéronef sont configurés pour acquérir au cours du temps des mesures de la tempà%.iture de sortie des gaz 10 d'échappement du moteur. Le contexte d'acquisition de ces données peut être très variable. Par exemple, les mesures concernant le décollage acquises lors d'un premier vol j.C. |a journée où le moteur est démarré à froid peuvent être différentes de celles acquises lors des autres vols de la journée. D'autres exemples concernent la variation des conditions météo (pluie, neige, gel, etc.), le changement de pilote, le lieu de survol 15 (au-dessus de la mer, du désert, de la terre, etc.). Ainsi, les mesures EGT sont très dépendantes des conditions extérieures. Avantageusement, les moyens de traitement 11 sont configurés pour réaliser une double normalisation sur ces mesures EGT par rapport à un référentiel standard ainsi que par rapport au contexte afin de supprimer l'influence des conditions 20 extérieures. Plus particulièrement, les mesures de température (EGT) sont normalisées par rapport à une température de référence iso standard formant ainsi des mesures de température normalisées. La température de référence est définie par rapport à une température de décollage mesurée au niveau de la mer. Ces mesures de tern' 25 normalisées sont en outre standardisées par rapport aux données de conte,J:e pour former des mesures EGT standardisées. La technique de standardisation est par exemple décrite dans le brevet la demanderesse EP2376988 et repose en particulier sur un modèle de régression en prenant éventuellement en considération des paramètres supplémentaires construits à partir de calculs utilisant des paramètres 30 exogènes initiaux. La Fig. 2B illustre les mesures EGT standardisées du moteur d'aéronef pendant un ensemble de vols successifs. La représentation graphique de la Fig. 2B correspond à la standardisation de celle de la Fig. 2A. Les deux représentations ont la même échelle et sont centrées sur la valeur moyenne des mesures EGT. Les traits discontinus horizontaux représentent les bornes +3o- et +6o- autour de la valeur moyenne. Enfin, les moyens de traitement 11 sont configurés pour calculer des marges entre les mesures EGT standardisées et une valeur de température maximale prédéterminée pour construire le signal temporel de marge résiduelle de température de sortie des gaz d'échappement comme illustré sur la Fig. 3. La ',empérure maximale correspond en général à une température limite requérant la dépose du moteur. en atelier. En particulier, la Fig. 3 montre que le signal temporel de marge résiduelle de température présente des oscillations aléatoires malgré le fait qu'il est basé sur des rrsur::'s EGT standardisées. Ces oscillations sont essentiellernent dues aux dii;',rentes interventions (par exemple, nettoyages à l'eau) effectuées sur le moteur. Les moyens de traitement 11 sont en outre configurés pour lisser le signal temporel de la Fig. 3 afin de former une première courbe Cl représentative de la marge résiduelle de température comme illustré sur la Fig. 4. Le lissage du signal temporel peut être réalisé selon une technique connue de filtrage (moyenne mobile, filtre gaussien, etc.). La première courbe Cl montre clairement des sauts montants correspondant à des augmentations artificielles de la marge résiduelle de température. Ces sauts montants semblent se présenter de manière plus ou moins régulière mais avec des effets différents. Conformément à l'invention, il est proposé de supprimer les sauts montants pour ne conserver que les parties descendantes représentatives de l'usure réelle du moteur. Ainsi, les moyens de traitement 11 sont configurés pour identifier automatiquement les morceaux descendants dans la première courbe. Plus particulièrement, l'identification des morceaux descendants peut par exemple être réalisée en recherchant les remontées correspondant aux sauts. Selon un mode de réalisation préféré de l'invention, les moyens de traitement 11 appliquent un modèle statistique correspondant à un processus statistique d'accroissement sur la première courbe Cl pour extraire les zones dépourvues de sauts montants. Ce processus statistique est décomposable en deux parties indépendantes de sorte qu'un accroissement dXt = Xt+1.- Xt de la marge résiduelle de température relatif à la première courbe Cl, est défini par la somme d'une première fonction décroissante Ut et d'une ,uxième fonction échelon Ht selon l'équation suivante: dXt = Ut + Ht La fonction décroissante Ut est représentative d'une usure habituelle du moteur d'aéronef et peut être considérée comme une variable aléatoire Ut qui suit une distribution gaussienne Ut-N(-u, au) paramétrée par une espérance négative -u représentant la décroissance naturelle de la marge résiduelle de température et par un écart-type au. La fonction échelon Ht est représentative des interventions sur le moteur d'aéronef et peut être définie par un produit Ht = ZtGt d'une loi boolienne aléatoire Zt et d'une variable aléatoire de saut positif Gt. La variable aléatoire de saut Gt peut aussi être décrite suivant une distribution gaussienne Gt N +g, crg) paramétrée par un saut moyen g positif et par un écart-type ag, le saut moyen g étant supérieur à la valeur absolue u de la décroissance naturelle -u de la marge résiduelle de température. La loi boolienne Zt déclenche le saut Gt ,!.e manière aléatoire avec une petite probabilité p et peut par exemple être définie par une loi binomiale Zt-B(p) de paramètre p prédéterminé. Ainsi, la probabilité d'une intervention sur le moteur (i.e., une remontée) est donnée par p = P(Zt = 1). Alors, pour identifier les morceaux descendants modélisés par la variable aléatoire Ut, il suffit de prendre l'ensemble des instants t où Zt = O en supprimant les points correspondants aux remontées, La Fig. 5 illustre une représentation graphique des remontées correspondant aux sauts. Cette représentation montre des points localisés autour des instants isolés correspondant à actions ponctuelles, séparés par des longs intervalles. Les instants isolés présentent des épaisseurs artificielles qui sont induites par les opérations de normalisation et de lissage. On notera que ces instants isolés sont beaucoup moins nombreux que les instants de décroissance et peuvent donc être ignorés sans incidence sur la valeur des paramètres -u et au la distribution gaussienne Ut-N(-u,o-u). Une fois que les instants isolés de remontées sont supprimés, les moyens de traitement 11 sont configurés pour construire une deuxième courbe C2 par concaténation des morceaux descendants, comme illustré sur la Fig. 6. La deuxième courbe est alors continue et décroissante tout en étant restreinte aux seuls morceaux descendants de première courbe Cl. La construction par concaténation de la deuxième courbe C2 consiste à déplacer verticalement et horizontalement les morceaux descendants pour coller les bouts entre eux manière à assurer une continuité entre les différents morceaux. En particulier, les moyens de traitements 11 sont configurés pour commencer par le dernier morceau (i.e., le morceau le plus récent) pour que les dernières valeurs aient un sens physique. Ensuite, les autres morceaux sont concaténés en remontant dans le temps de proche en proche. Ainsi, chaque morceau précédent (i.e., antérieur) est c '_Vacé pour le faire joindre au morceau suivant de manière successive dans le sens inverse du temps. Autrement dit, on ajoute progressivement les valeurs des variations dXF = Xt_i - Xt de droite à gauche pour former la deuxième courbe C2. En outre, les moyens de traitement 11 sont configurés pour construire à partir de la deuxième courbe C2 un modèle de prédiction de défaillance permettant de déterminer au moins un indicateur de pronostic de ,.::,faillance du moteur. On notera que la courbe illustrée dans l'exemple de la Fig. 6 est pratiquement linéaire. Ainsi, selon un premier mode de réalistiol, il est avantageux d'utiliser l'historique de la deuxième courbe C2 pour construire un modèle linéaire de l'évolution de la marge résiduelle de température et d')ppliquer ensuite sur ce modèle linéaire une technique de régression assez simple pour déterminer le ou les indicateur(s) de pronostic de défaillance. Selon ce mode de réalisation, il n'est pas nécessaire d'analyser la dynamique du signal et il suffit d'estimer la tendance générale de la décroissance qui suit une loi de type Student. Toutefois, pour des moteurs assez jeunes, la décroissance de la marge est lente et n'est pas strictement linéaire. Alors, selon un deuxième mode de réalisation, on utilise une méthode stochastique basée sur des filtres dynamiques ou particulaires. De manière générale, un bon moyen pour anticiper un processus continu consiste à modéliser son comportement avec un modèle autonAgrrssif. Ainsi, on construit un modèle autorégressif modélisant l'évolution de la marge résiduelle de température en utilisant l'historique de la deuxième courbe. Alors, il est possib!e d'extraire un processus d'états (Xt)t>0 dits cachés à partir des observations (Y)t En particulier, on suppose que le processus d'états (Xt)t>0 est une chaîne de Markov d'ort!re 1 et que le lien entre (Xt)t>0 et (Y3t>1 est régi par une hypothèse de canal sans ME moire. Alors, l'espace d'états du modèle autorégressif peut être défini de la manière suivante : Xt = FX_i+ et Yt - HXt_i et est un bruit blanc, l'opérateur F est la matrice de transition d'états et l'opérateur H est la matrice d'observation définissant le système dynamique. On notera que grâce à la suvression des sauts montants, le système dynamique qui capte la décroissance du signal ce la marge résiduelle de température est très simple à analyser et à mettre en oeuvre et consomme très peu d'étapes de calcul tot.:': on fournissant des résultats très précis. Ensuite, on applique un filtre dynamique sur le modèle autorégressif pour estimer de façon récursive le processus d'états cachés (Xt)t>0 à partir des observations (Y3t>1 selon une technique bayésienne. Le filtrage dynamique permet ainsi de déterminer pour tout instant k l'état caché Xk à partir des observations disponibles Y1, , Yk (i.e., jusqu'à l'instant k). On peut ainsi déterminer ensuite le ou les indicateur(s) de pronostic de .!éfaillance à long terme et avec précision. En effet, la Fig. 7 illustre l'application d'un filtre particulaire sur le système dynamique de la deuxième courbe. La courbe C2 jusqu'à la ligne verticale en traits discontinus dl correspond au processus d'observations (Y) représentant l'évolution observée de la marge. La ligne verticale dl correspond ainsi à l'instant où l'on fait une prédiction en fonction du processus d'observations (Y). A partir de cet instant, l'application du filtre dynamique simule une pluralité de trajectoires particulaires tl, la courbe en gras C3 à l'intérieur , el'i',,:rentes trajectoires tl représentant la trajectoire moyenne. La ligne horizontale d2 représente le seuil de défaillance. Il est à noter que la qualité de la prédiction est très bonne grâce à la connaissance précise du système dynamique et notamment grâce à la suppression des remontées correspondant aux interventions ponctueles s:ir le moteur. La Fig. 8 illustre d'un filtre particulaire sur le système dynamique relatif à un moteur jeune. Au départ, la diminution de la marge résiduelle de température est assez forte à cause d'un phénomène de rodage du jeune moteur. Ensuite, l'évolution s'adoucit et la décroissance devient assez lente. De la même manière que dans la Fig. 7, l'application du filtre dynamique simule une pluralité de .ectoires particulaires avec une grande qualité de prédiction. Ainsi, l'utilisation d'un modèle,.,,n;-,mique à filtre particulaire permet de faire des pronostics précis à long terme pour tous types de moteurs et à tout âge. Il est à noter que le filtre dynamique peut être un filtre bayesien linéaire ou non-linéaire. En variante, on peut utiliser un filtre de Kalman, ou une extension d'un filtre de Kalman (par exemple, un filtre de Kalman étendu). L'application quelconque de l'un de ces filtres dynamiques sur le modèle de prédiction permet d'estimer l'évolution de la marge résiduelle de température et par conséquent de déterminer des indicateurs de pronostic de défaillance peuvent être utilisés comme des indicateurs d'alerte. La Fig. 9 illustre de manière schématique deux indicateu ;1-onostic de défaillance, selon l'invention. Le premier indicateur 11 estime la probabilité de franchissement du seuil de défaillance d2 avant un horizon temporel hl prédéterminé. Ceci est schématisé par le calcul d'une probabilité de détection POD (Probability Of Detection) à l'instant t+h. A titre d'exemple, l'horizon temporel hl peu'. correspondre à une date programmée pour. l'inspection du moteur. En variante, on peut raisonner en termes d'estimation de la durée de vie restante RUL (Remaining Useful Life) du moteur. Dans ce cas, le deuxième indicateur estime la date de franchissement du seuil de c,:âfaillance d2. Ainsi, ces indicateurs de pronostic c,-.? ,..r.i11,-nce 11, 12 permettent de prévoir la probabilité de défaillance une fois arrivé à un horizon futur hl, ou de prédire la date de défaillance. L'invention vise aussi un programme d'ordinateur, susceptible d'être mis en oeuvre dans les moyens de traitement et comportant des instructions de code adaptées à la mise en oeuvre d'un procédé selon l'invention telle que décrite ci-dessus. REVENDICATIONS1. Procédé de surveillance d'un moteur d'aéronef, caractérisé en ce qu'il comporte les étapes suivantes : - acquisition d'un signal temporel de marge résiduelle de température de sortie des gaz d'échappement dudit moteur (2) d'aéronef, - lissage dudit signal temporel pour former une première courbe (Cl) représentative de ladite marge résiduelle de température, - identification des morceaux descendants dans ladite première courbe, - construction d'une deuxième courbe (C2) par concaténation desdits morceaux descendants, ladite deuxième courbe étant continue tout en étant restreinte auxdits morceaux descendants de ladite première courbe, et -construction d'un modèle de prédiction à partir de ladite deuxième courbe pour déterminer au moins un indicateur d'? pronostic d^ &eaillance (11,12). 15 2. Procédé selon la revendication 1, caractérisé en ce que l'identification des morceaux descendants dans ladite première courbe comporte les étapes suivantes : - application sur la première courbe (Cl) d'un modèle statistique d'accroissement décomposable en deux p::r.ties in;!é;:enciintes composées d'une 20 première fonction décroissante représentative d'une usure habituelle du moteur d'aéronef et d'une deuxième fonction échelon composée de sat.:i:s cklenchés de manière aléatoire et représenteq-`:i; d'interventions ponctuelles sur le moteur d'aéronef, - recherche (:e remontées correspondantes auxdts sauts, et 25 - identification des morceaux descendants en supprimant les points desdites remontées sur la première courbe. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que la construction par concaténation de ladite deuxième courbe (C2) comporte un collage desdits 30 morceaux descendants en déplaçant chaque morceau précédent, pour la faire joindreau morceau suivant en commençant par le dernier morceau et en remontant dans le temps de proche en proche. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la construction dudit modèle de prédiction comporte les étapes suivantes: - construction d'un modèle au.orégressif me,..1élisant l'évolution de la marge résiduelle de température en utilisant nistorique de ladite deuxième courbe, et - application d'un filtre dynamique sur ledit modèle autorégressif pour déterminer ledit au moins un indicateur de pronostic de défaillance. 5. Procédé selon la revendication 4, caractérisé en ce que le filtre dynamique est sélectionné parmi l'ensemble des filtres particulaires suivants : filtre bayesien, filtre de Kalman, filtres de Kalman étendus. 6. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que la construction dudit modèle de prédiction comporte les étapes suivantes : - construction d'un modèle linéaire de l'évolution de la marge résiduelle de température en utilisant l'historique de laJite deuxième courbe, et - application d'une technique de régression sur ledit modèle linéaire pour déterminer ledit au moins un indicateur de pronostic de défaillance. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que ledit au moins un indicateur de pronostic de défaillance est sélectionné parmi un ensemble d'indicateurs comprenant : -un premier indicateur (11) d'estimation de la probabilité de franchissement d'un seuil de défaillance avant un horizon temporel prédéterminé, et -un deuxième indicateur (12) d'estimation de la date de franchissement d'un seuil de défaillance.30 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce cit. e l'acquisition dudit signal temporel de marge re,siduelle de température comporte les étapes suivantes : - acquisition au cours du temps des mesures de la température de sortie des gaz d'échappement du moteur d'aéronef, - normalisation desdites mesures de température par rapport à une température de référence iso standard formant ainsi mesures de température normalisées, - standardisation desdites mesures de température normalisées en prenant en compte des c'»rl ées de contexte formant ainsi des mesures de température standardisées, et - calcul des marges entre lesdites mesures de terrp,rature standardisées et une valeur de température maximale prédéterminée (en fonction du moteur) pour former ledit signal temporel de marge. 9. Système de surveillance d'un moteur d'aéronef, caractérisé en ce qu'il comporte : - des moyens d'acquisition et de traitement (11) configurés pour recueillir un signal temporel de marge résiduelle de température de sortie des gaz d'échappement dudit moteur (2) d'aéronef, - des moyens d'acquisition et de traitement (11) configurés pour lisser ledit signal tem ,orel formant ainsi une première courbe (Cl) représentative de ladite marge résiduelle de température, - des moyens d'acquisition et de traitement (11) configurés pour identifier des 2 morceaux descendants dans ladite première courbe, - des moyens d'acquisition et traie.?.ment (11) configurés pour construire une deuxième courbe (C2) par concaténation desdits morceaux descendants, ladite deuxième courbe étant continue tout en étant restreinte auxdits morceaux descendants de ladite première courbe,- des moyens d'é uisition et de traitement (11) configurés pour construire un modèle de prédiction à partir ladite deuxième courbe pour déterminer au moins un indicateur de pronostic de défaillance (11,12). 10. Programme d'ordinateur susceptible d'être mis en oeuvre par les moyens de traitement et comportant des instructions de code adaptées à la mise en oeuvre du procédé de surveillance selon l'une quelconque des revendications 1 à 8.