FR 3067124 A1 20181207 FR 1754949 A 20170602 - 1 - « Procédé et système pour diagnostiquer en temps réel l'état de fonctionnement d'un système électrochimique, et système électrochimique intégrant ce système de diagnostic » Domaine technique La présente invention concerne un procédé pour diagnostiquer en temps réel l'état de fonctionnement d'un système électrochimique. Elle vise également un système de diagnostic mettant en oeuvre ce procédé, ainsi qu'un système électrochimique, notamment une pile à combustible, intégrant ce système de diagnostic. Etat de la technique antérieure Les piles à combustible actuelles, qu'elles soient embarquées ou stationnaires, doivent faire l'objet d'une surveillance de leur état, afin 15 notamment de prévenir une défaillance ou une dégradation de leur fonctionnement. Il existe déjà des procédés de détection de défauts dans une pile à combustible, qui ne peuvent être mis en oeuvre que hors fonctionnement de cette pile et dans un contexte de laboratoire. Ces méthodes mettent parfois en 20 oeuvre des capteurs intrusifs. Par ailleurs, pour des piles à combustible à forte valeur ajoutée, utilisées par exemple en milieu nécessitant un fonctionnement certain à tout moment, tel qu'en milieu hospitalier par exemple, des méthodes de détection impliquant une mise hors fonctionnement temporaire ne sont pas 25 envisageables. Le but de la présente invention est de proposer un nouveau procédé de détection de défauts dans un système électrochimique, notamment une pile à combustible, qui permette de connaître l'état de ce système en fonctionnement et sans capteur intrusif et/ou complexe ou coûteux. 30 Exposé de l'invention Cet objectif est atteint avec un procédé pour diagnostiquer en temps réel l'état de fonctionnement d'un système électrochimique comprenant un empilement de cellules, ce procédé comprenant des étapes pour réaliser des 3067124 - 2 - mesures de tension desdites cellules, caractérisé en ce qu'il comprend en outre : - un traitement des mesures de tension ainsi réalisées en temps réel pour en extraire des formes d'onde spécifiques, 5 - une transformation desdites formes d'onde spécifiques pour en générer des points spécifiques du fonctionnement temps-réel dudit système électrochimique, et - une comparaison desdits points spécifiques du fonctionnement temps-réel et de points spécifiques d'un fonctionnement hors-ligne dudit 10 système électrochimique qui sont issus d'une transformation de formes d'onde spécifiques extraites de mesures de tension effectuées hors ligne alors que ledit système électrochimique est placé dans des états de fonctionnement connus incluant des états de défaut, de façon à produire des informations sur l'état de fonctionnement temps-réel dudit système électrochimique. Les étapes de mesure de tension de cellule peuvent être effectuées périodiquement. Dans un mode de réalisation préféré de l'invention, les formes d'onde spécifiques extraites des mesures de tension effectuées hors-ligne et/ou en temps réel sont soumises à une transformation mathématique - appelée « Shapelet Transform » - pour créer des points spécifiques à partir de formes d'onde spécifiques. Ces points spécifiques peuvent être identifiés dans un espace bidimensionnel (2D) puis classés dans un espace tridimensionnel (3D) Dans un premier mode particulier de réalisation du procédé selon l'invention, des points spécifiques sont classés au moyen d'une fonction de type SSM-SVM (pour « Sphere Shaped Multi-Class - Support Vector Machine »). Dans un second mode particulier de réalisation du procédé selon l'invention, des points spécifiques sont classés au moyen d'une fonction de type kNN (pour « K Nearest Neighbor »). Dans un troisième mode particulier de réalisation du procédé selon l'invention, des points spécifiques sont classés au moyen d'une fonction de type GMM (pour « Gaussian Mixture Model »). 3067124 - 3 - Le procédé de diagnostic selon l'invention peut en outre avantageusement comprendre une constitution initiale d'une base de données par apprentissage sur des états de fonctionnement connus du système électrochimique, qui sont testés hors ligne. 5 Lorsque le procédé selon l'invention est mis en oeuvre pour diagnostiquer l'état de fonctionnement d'une pile à combustible, les états de fonctionnement connus qui sont testés hors ligne comprennent tout ou partie des états suivants : Fonctionnement normal, Pression basse hydrogène, Stoechiométrie basse du côté air, Asséchement, Noyage. 10 Suivant un autre aspect de l'invention, il est proposé un système pour diagnostiquer l'état de fonctionnement en temps réel d'un système électrochimique comprenant un empilement de cellules, mettant en oeuvre le procédé de diagnostic selon l'invention, ce système de diagnostic comprenant des moyens pour mesurer les tensions aux bornes desdites cellules, 15 caractérisé en ce qu'il comprend en outre : - des moyens pour traiter les tensions ainsi mesurées pour en extraire des formes d'onde spécifiques, des moyens pour transformer lesdites fonctions d'onde spécifiques pour en générer des points spécifiques du fonctionnement temps-réel dudit système électrochimique, et 20 - des moyens pour comparer lesdits points spécifiques du fonctionnement temps-réel et des points spécifiques d'un fonctionnement hors-ligne dudit système électrochimique qui sont issus d'une transformation de formes d'onde spécifiques extraites de mesures de tension effectuées hors ligne - alors que ledit système électrochimique est placé dans des 25 états de fonctionnement connus incluant des états de défaut, de façon à produire des informations de diagnostic de l'état de fonctionnement temps-réel dudit système électrochimique. Le système de diagnostic selon l'invention peut en outre avantageusement comprendre des moyens pour transformer les formes 30 d'onde spécifiques extraites des mesures de tension effectuées hors-ligne et/ou en temps réel, de façon à créer des points spécifiques. Ce système de diagnostic peut en outre comprendre des moyens pour identifier les points spécifiques dans un espace bidimensionnel (2D) et des moyens pour classer lesdits points spécifiques ainsi identifiés dans un espace 3067124 - 4 - tridimensionnel (3D), ainsi qu'une base de données par apprentissage sur des états de fonctionnement connus du système électrochimique. Suivant encore un autre aspect de l'invention, il est proposé un système électrochimique intégrant un système de diagnostic selon l'invention. Ce 5 système électrochimique peut notamment comporter une pile à combustible. Description des figures D'autres avantages et particularités de l'invention apparaîtront à la lecture de la description détaillée de mises en oeuvre et de modes de 10 réalisation nullement limitatifs, et des dessins annexés suivants : - la figure 1 illustre schématiquement le principe de mesure de tension mis en oeuvre dans le procédé de détection de défauts selon l'invention ; - la figure 2 illustre schématiquement une mise en oeuvre d'un microcontrôleur pour des mesures analogiques ; 15 - la figure 3 représente un dispositif de mesure de tension mis en oeuvre dans un exemple de réalisation d'un système de détection de défauts selon l'invention ; - la figure 4 représente des formes d'onde spécifiques (ou « shapelets ») produits dans le procédé de détection de diagnostic selon l'invention ; 20 - la figure 5 illustre une classification des formes d'onde spécifiques avec la méthode SSM-SVM ; et - la figure 6 illustre schématiquement un exemple de réalisation du procédé de diagnostic selon l'invention. 25 Modes de réalisation détaillés Ces modes de réalisation étant nullement limitatifs, on pourra notamment considérer des variantes de l'invention ne comprenant qu'une sélection de caractéristiques décrites ou illustrées par la suite isolées des autres caractéristiques décrites ou illustrées (même si cette sélection est 30 isolée au sein d'une phase comprenant ces autres caractéristiques), si cette sélection de caractéristiques est suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique antérieure. Cette sélection comprend au moins une caractéristique de préférence fonctionnelle sans détails structurels, et/ou avec seulement une partie des détails structurels si cette partie uniquement est suffisante pour 3067124 - 5 - conférer un avantage technique ou à différencier l'invention par rapport à l'état de la technique antérieure. On va tout d'abord décrire, en référence à la figure 1, un dispositif 10 de mesure des tensions aux bornes des différentes cellules d'un système 5 électrochimique 1 tel qu'une batterie, un supercondensateur ou une pile à combustible. Les différentes tensions Vcell 1,.., k+1 sont mesurées puis traitées par des composants programmables de type DSP (pour « Digital Signal Processor ») ou microcontrôleur 20, 70 et les données sont transmises par un 10 bus de communication CAN 6. Le dispositif de mesure 10 se présente sous forme modulaire. Chaque module 2, 7 isolé électriquement est chargé de mesurer les tensions Vcell_1,Vcell_2,Vcell_3 ; Vcell_k,Vcell_k+1 d'un nombre de cellules en respectant l'ordre défini par la connexion en série de celles-ci. Il est donc 15 possible de connecter autant de modules 2, 7 que nécessaire pour répondre aux besoins de l'utilisateur, la limite étant fixée par le niveau d'isolation en tension des composants électroniques (environ 1000V dans le cas le plus critique). Chaque module de mesure 2, 7 comporte un circuit électronique basé 20 autour d'un microcontrôleur 20, 70, son potentiel de référence (masse) correspondant au potentiel le plus faible du paquet d'éléments (Vréf_1, Vréf_2). Ces modules 2, 7 sont alimentés par l'intermédiaire d'un convertisseur DC/DC isolé 3, 8 connecté sur l'alimentation générale du dispositif 10. Les différentes données issues des mesures sont transmises par 25 l'intermédiaire d'un bus de communication. Pour cela, des drivers CAN isolés 5, 9 permettent des liaisons entre les microcontrôleurs 20, 70 et le bus CAN 6 du dispositif de mesure 10. Les module de mesure 2, 7 sont dupliqués autant de fois que nécessaire pour mesurer les différentes tensions aux bornes de tous les éléments du 30 système électrochimique 1 étudié. En référence à la figure 2, le potentiel de référence des mesures analogiques du microcontrôleur 30 est connecté à la tension la plus faible de l'ensemble des éléments concernés n, n+1, n+2. Les différentes mesures de tension Vn, Vn+1, Vn+2 sont connectées aux entrées analogiques AN_1, 35 AN_2, AN_3 du microcontrôleur 30 en respectant l'ordre établi par les 3067124 - 6 - potentiels : du plus faible au plus fort. La tension aux bornes d'un élément est alors déterminée par la différence de deux tensions consécutives mesurées. En fonction des niveaux de tension mesurés, les entrées analogiques du microcontrôleur 30 peuvent être connectées directement sur les éléments à 5 mesurer. Pour des tensions plus élevées, un circuit d'adaptation permet d'ajuster la tension dans une plage compatible avec celle du microcontrôleur. En fonction du microcontrôleur ou du circuit DSP utilisé, le nombre de voies analogiques est variable et limite ainsi le nombre de tensions mesurées sur un sous-ensemble de cellules. Les mesures analogiques peuvent être 10 effectuées de manière simultanée (mesures synchrones) ou séquentielle (mesures asynchrones). Dans l'exemple de réalisation illustré par la figure 3, le dispositif de mesure 100 est utilisé pour la mesure des tensions des cellules d'une pile à combustible 30. En pratique, la tension des cellules peut atteindre environ 15 1.43V lors de leur activation, mais la plage de mesure utile est comprise entre 0 et 1V. Le dispositif de mesure 100 comprend un microcontrôleur 50 comportant 11 entrées analogiques et un seul convertisseur analogique/numérique intégré. On utilise alors le principe de mesures 20 séquentielles pour les différentes voies analogiques. La durée globale de conversion pour les 11 voies reste inférieure à 2ms (150ps par voie, soit 1,6 ms). Les différentes tensions des cellules de la pile à combustible 40 sont connectées sur les entrées analogiques, de manière directe pour les quatre 25 premières voies (tensions inférieures à 4V) et par l'intermédiaire d'un circuit d'adaptation pour les voies AN_4 à AN_10. L'alimentation du dispositif de mesure 100 est réalisée à partir d'une source de tension de 5V continue (+5v_1). Un convertisseur DC/DC 51 permet d'obtenir une tension de 5V (+5V_2) isolée pour l'alimentation de l'ensemble de 30 l'électronique. Le circuit utilisé assure une tension d'isolement de 1000V. L'ensemble du dispositif de mesure 100 se trouve ainsi isolé avec un potentiel de référence Vref_1. Pour assurer cet isolement, les communications entre le microcontrôleur 50 et le reste du système sont effectuées par un bus CAN 52 également isolé. 35 Un circuit driver CAN 53 isole le microcontrôleur 50 en assurant une tension 3067124 - 7 - d'isolement de 2500V. On va maintenant décrire, en référence aux figures 4 à 6, des étapes de traitement qui sont implémentées dans le procédé de diagnostic selon l'invention. 5 Le procédé de diagnostic selon l'invention est construit autour d'un apprentissage basé sur des données mesurées et connues. Il est ici mis en oeuvre dans le but de détecter des défauts d'un système de pile à combustible (PAC) de type PEMFC (pour « Polymer Electrolyte Membrane Fuel Cell »). Il est entendu que le procédé de diagnostic selon l'invention peut être mis en oeuvre 10 sans limitation dans d'autres types de piles à combustible. Ce procédé de diagnostic met en oeuvre un dispositif de mesure de tension individuelle des cellules unitaires de la PAC tel que décrit précédemment. Afin de tenir compte de la dynamique du système et de l'hétérogénéité spatiale des différentes cellules qui le composent, le procédé de 15 diagnostic selon l'invention fait appel à un outil d'analyse de série temporelle, appelé dans la suite « Shapelet transform », qui réalise une caractérisation de formes d'onde spécifiques extraites des mesures de tension aux bornes des cellules du système électrochimique. En référence à la figure 6, dans une forme particulière de réalisation du 20 procédé de diagnostic selon l'invention, celui-ci comprend un processus 6A de traitement en temps réel de mesures de tension effectuées sur les cellules de la pile, et un processus 6B de traitement hors ligne de mesures de tension effectuées sur les mêmes cellules mais dans des états de fonctionnement connus incluant des états de dysfonctionnement. 25 L'outil d'analyse temporelle est utilisé pour extraire des formes d'onde spécifiques ou « shapelets » 40 (Figure 4) issues des données d'une base de données 61 (Figure 6) servant au diagnostic en temps-réel. L'outil d'analyse temporelle est utilisé dans une fenêtre 41 dont la durée est connue et définie préalablement en fonction des défauts considérés. Dans un second temps, un 30 outil de classification SSM-SVM est appliqué dans l'espace dimensionnel considéré (fonction du nombre de défauts). Les tensions des cellules du système de pile à combustible sont par nature non homogènes de par leurs états physiques au sein d'une pile (ou « stack »), tels que la température, la distribution des gaz, le taux d'humidité 35 relative, la distribution du catalyseur sur la couche d'activation, et on doit 3067124 - 8 - considérer la corrélation temporelle entre les échantillons mesurés. Les tensions de cellules échantillonnées dans un intervalle de temps (fenêtre 41) sont considérées comme des variables pour le procédé de diagnostic selon l'invention. 5 Cette fenêtre d'acquisition est définie mathématiquement comme suit : T = [T(1), ... , T(Iw)] = [v1, ... , vlw] (1) OÙ y; est le vecteur échantillonné à l'index i et composé de l'ensemble des tensions d'un nombre ncell decellules. 10 Ce vecteur vi = [vi (1), vi (2), ... , vi (nceii)]. Lest la largeurde la fenêtre d'acquisition. Chaque observation T est une série chronologique dynamique qui peut être représentée mathématiquement sous la forme d'une matrice, dont la dimension esr définie d'une part par le nombre de cellules et d'autre part par la 15 durée de l'observation lw et dont le nombre de points est défini par la période d'échantillonnage (ou fréquence d'acquisition) du dispositif de mesure. Dans le processus hors-ligne 6B, une base de données d'apprentissage 71 est constituée à partir de N mesures hors-ligne obtenues au moyen du dispositif de mesure décrit précédemment. Elle est répertoriée par classes 20 notées S20, S21, S22,..., QC qui correspondent aux différents défauts. La classe d'indice 0 est définie comme l'état de fonctionnement normal (sans défauts). L'espace des classes sur l'ensemble des observations Ti est connu au préalable, soit g; e {0,1, ... , C}. Les nombres d'observations 78 qui 25 constituent la base de donnée 71 dans l'ensemble de ces classes sont respectivement NO, Ni, ...,NC et satisfont l'équation NO + Ni + ... + NC = N. On va maintenant plus particulièrement décrire la génération de formes d'onde spécifiques (« Shapelets ») dans les deux processus temps-réel 6A et hors-ligne 6B. 30 Pour identifier une forme d'onde spécifique, il est nécessaire de produire une quantité de signaux éligibles que l'on appelle candidats et qui sont généré dans le vecteur SCji. En faisant l'hypothèse que la forme d'onde spécifique est d'une durée ls, les candidats sont sélectionnés en tenant compte de la durée 3067124 - 9 - d'observation lw. Par exemple, on considère des candidats jusqu'à lw - ls + 1 pour une observation Ti suivant la formule ci-après : SCi = {Ti (.1), Ti (.1 + 1), , Ti (j + Is - 1)}, j = 1, , 1w-1s + 1 (2) 5 OÙ SCJi représente le jiéme candidat de la matrice Ti. Comme SCji est issu de Ti, l'espace dans lequel il se situe est également gi. Le nombre de candidats générés à partir de la base d'apprentissage 71 est donc (lw - Is + 1)N. Dans ce procédé, les durées Is et lw sont définies préalablement pour 10 chaque classe. Cette définition est faite de manière empirique lors de la constitution de la base de données d'apprentissage. Pour pouvoir classifier et caractériser les candidats, trois opérations sont nécessaires. La première consiste à calculer la distance entre un candidat SCij et l'observation Tk, qui est égale à la valeur minimum entre le candidat SCij et 15 l'ensemble de candidats générés par cette observation Tk comme suit : ) = min ILc1 - SCk,i (3) Pour chaque candidat cette distance est calculée, puis la distance moyenne de cette même classe est également définie. pk=gi CIIS(SCI,d) = 20 Grâce aux deux valeurs précédentes, on calcule le facteur de qualité de chaque candidat d'une même classe suivante la formule de discrimination quadratique suivante : Qu SC) (5) k =y 25 Lorsque tous les candidats de chaque classe ont été évalués, ils sont classés par ordre décroissant suivant leurs facteurs de qualité. Dans le cas où deux formes d'onde spécifiques candidats de même la base d'observation obtiennent des indices redondants, c'est alors celui qui présente le facteur de qualité le plus élevé qui est conservé. Dans les candidats restants, ce sont 30 ceux qui ont le meilleur facteur de qualité qui sont sélectionnés. La procédure d'extraction des formes d'onde spécifiques qui vient d'être décrite peut aisément être mise en oeuvre sous la forme d'un algorithme avec Nt (4) 3067124 - 10 - des techniques logicielles connues. On peut ainsi générer des formes d'onde spécifiques 74 dans le processus hors-ligne 6B et des formes d'onde spécifiques 40 dans le processus temps-réel 6A. On va maintenant plus particulièrement décrire un exemple d'extraction 5 des formes d'onde spécifiques lors d'un diagnostic temps réel. Pour une fenêtre de diagnostic T (Figure 4) prise arbitrairement, les caractéristiques des formes d'onde spécifiques de cette base arbitraire issue de mesures en temps réel 64 (voir figure 6) sont déterminées en exécutant le calcul de distance avec ceux issus de la base d'apprentissage 71 et sélectionnés comme 10 pertinents. Ce calcul permet également de connaître la dimension du domaine de caractérisation qui servira à l'initialisation de l'outil 63 de classification SSMSVM en référence à la figure 6. On va maintenant décrire un exemple de mise en oeuvre de cet outil de 15 classification SSM-SVM dans le procédé de diagnostic selon l'invention. Pour solutionner le problème de classification non linéaire issue des données d'apprentissage (notamment dû au fait que l'on considère plus de deux classes), la méthode utilisée ci-après est connue sous le nom « kernel trick ». 20 Dans un premier temps, on réalise une projection des vecteurs dans un espace multidimensionnel via la fonction 4, . Les vecteurs sont définis comme suit : e n Par exemple, pour une classe i, la sphère est calculée comme suit : (5) 25 Où Rî et ai sont respectivement le rayon et le centre du cercle de la sphère i. La fonction ai = (11) r - D et D (12) Le rayon de la sphère i, R' peut être calculé d'après les conditions de 25 Karush-Kuhn-Tucker ( « KKT conditions ») décrites ci-après. 1 (1195 (4112 R2) e ) 0 - a Donc 110 pour tout les zn avec 30 A partir de ces calculs analytiques, l'algorithme SSM-SVM peut être codé. Dans l'exemple de mise en oeuvre décrit ici, on opère sur un nombre de classe supérieur à 2 et dans un espace de dimensions supérieures à 2. 3067124 - 12 - Cet algorithme permet de décomposer le problème quadratique en une somme de problèmes linéaires qui sont traités un à un. Chaque multiplicateur de Lagrange est calculé individuellement et les vecteurs supports décrivant les sphères englobant les classes sont mis à jours à chaque itération. 5 La méthode SSM-SVM qui vient d'être décrite permet ainsi de classifier des formes d'onde spécifiques, en référence à la figure 5. Par ailleurs, l'association des formes d'onde spécifiques (ou « shapelets ») et de la classification SSMSVM permet d'avoir des performances en terme de taux d'erreur supérieures à celles habituellement observées avec les procédés actuels de diagnostic. 10 L'utilisation de la méthode SSM-SVM permet d'obtenir à la fois des résultats de classification 77 des formes d'onde spécifiques 75 issues d'observations hors-ligne dans un mode d'apprentissage, et des résultats de classification 66 des formes d'onde spécifiques 65 issues de mesures de tension effectuées en temps réel. 15 Les résultats de classification issus du mode d'apprentissage SSM-SVM 76 sont exploités dans le traitement SSM-SVM 63 des points spécifiques issues de la transformation 62 des formes d'onde spécifiques issues des mesures temps-réel. On va maintenant décrire, en référence à la figure 6, une stratégie de 20 diagnostic en temps réel mise en oeuvre dans le procédé de diagnostic selon l'invention. Les tensions aux bornes des cellules du système électrochimique sont mesurées au moyen du dispositif de mesure décrit précédemment en référence aux figures 1 à 3. Ces mesures sont numérisées à une fréquence d'échantillonnage de 1hz pendant un temps Lw puis stockées sous forme d'un 25 tableau de données temps-réel 61 qui sera ensuite utilisé pour les calculs. Les données en temps réel sont acquises selon une fenêtre glissante d'acquisition. Une fois le tableau stocké, les distances minimum sont alors calculées via la transformation 62 en forme d'onde spécifique. Lorsque la forme d'onde spécifique est caractérisée et représentée 30 graphiquement par un point spécifique 65, il est ensuite analysé par l'algorithme SSM-SVM 63 qui va lui affecter une classe en fonction de sa position dans les volumes définis par les sphères, comme l'illustre la figure 5. Chaque classe est associée à un état de fonctionnement (Normal ou décrivant un défaut connu). 3067124 - 13 - Des règles de diagnostic 67 sont ensuite appliquées aux résultats de classification 66 de façon à générer des résultats de diagnostic 68 qui vont permettre de poser un diagnostic de fonctionnement 69 (normal ou défaut, type de défaut,..) qui sera communiqué à la gestion du système 5 électrochimique et possiblement à l'utilisateur ou au personnel chargé de la maintenance de ce système. Ces opérations sont répétées en boucle à chaque nouvelle acquisition (fenêtre) et ce à chaque mise en route du système. Il est à noter que d'autres méthodes mathématiques, faisant notamment appel aux réseaux de neurones ou encore de type kNN (K Nearest 10 Neighbor) ou GMM (Gaussian Mixture Model), pourraient être utilisées pour réaliser des classifications de formes d'onde spécifiques. Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. Bien entendu, les différentes 15 caractéristiques, formes, variantes et modes de réalisation de l'invention peuvent être associées les unes avec les autres selon diverses combinaisons dans la mesure où elles ne sont pas incompatibles ou exclusives les unes des autres. En particulier toutes les variantes et modes de réalisation décrits précédemment sont combinables entre eux. 20 REVENDICATIONS1. Procédé pour diagnostiquer en temps réel l'état de fonctionnement d'un système électrochimique (1, 30) comprenant un empilement de cellules, ce procédé comprenant des étapes pour réaliser des mesures de tension aux bornes desdites cellules, caractérisé en ce qu'il comprend en outre : - un traitement des mesures de tension ainsi réalisées en temps réel pour en extraire des formes d'onde spécifiques (40), - une transformation (62) desdites formes d'onde spécifiques (40) pour en générer des points (65) spécifiques du fonctionnement temps-réel dudit système électrochimique (1, 30), et - une comparaison desdits points (65) spécifiques du fonctionnement temps-réel et de points (75) spécifiques d'un fonctionnement hors-ligne dudit système électrochimique (1, 30) qui sont issus d'une transformation (73) de formes d'onde spécifiques (74) extraites de mesures de tension (78) effectuées hors ligne alors que ledit système électrochimique (1, 30) est placé dans des états de fonctionnement connus incluant des états de défaut, de façon à produire des informations (69) sur l'état de fonctionnement temps-réel dudit système électrochimique (1, 30). 2. Procédé selon la revendication 1, caractérisé en ce que les étapes de mesure de tension de cellule sont effectuées périodiquement. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que les formes d'onde spécifiques (40, 74) extraites des mesures de tension effectuées hors-ligne et/ou en temps réel sont soumises à une transformation mathématique (62, 73) - appelée « Shapelet Transform » - pour créer des points spécifiques (65, 75). 4. Procédé selon la revendication 3, caractérisé en ce que les points spécifiques sont identifiés dans un espace bidimensionnel (2D) puis classés dans un espace tridimensionnel (3D). 3067124 - 15 - 5. Procédé selon la revendication 4, caractérisé en ce que des points spécifiques sont classés au moyen d'une fonction (63, 76) de type SSM- SVM (pour « Sphere Shaped Multi-Class - Support Vector Machine »). 5 6. Procédé selon l'une des revendications 4 ou 5, caractérisé en ce que des points spécifiques sont classés au moyen d'une fonction de type kNN (pour « K Nearest Neighbor »). 7. Procédé selon l'une des revendications 4 à 6, caractérisé en que des points 10 spécifiques sont classés au moyen d'une fonction de type GMM (pour « Gaussian Mixture Model »). 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre une constitution initiale d'une 15 base de données par apprentissage (71) sur des états de fonctionnement connus du système électrochimique (1, 30), qui sont testés hors ligne. 9. Procédé selon l'une quelconque des revendications précédentes, mis en oeuvre pour diagnostiquer l'état de fonctionnement d'une pile à 20 combustible (30), caractérisé en ce que les états de fonctionnement connus testés hors ligne comprennent tout ou partie des états suivants : Fonctionnement Normal, Pression basse hydrogène, Stoechiométrie basse du côté Air, Asséchement, Noyage, Pression élevée. 25 10.Système pour diagnostiquer l'état de fonctionnement en temps réel d'un système électrochimique (1, 30) comprenant un empilement de cellules, mettant en oeuvre le procédé de diagnostic selon l'une quelconque des revendications précédentes, comprenant des moyens pour mesurer les tensions aux bornes desdites cellules, caractérisé en ce qu'il comprend en 30 outre : - des moyens pour traiter les tensions ainsi mesurées pour en extraire des formes d'onde spécifiques, - des moyens pour transformer lesdites formes d'onde spécifiques pour en générer des points spécifiques du fonctionnement temps-réel dudit 35 système électrochimique, et 3067124 - 16 - - des moyens pour comparer lesdits points spécifiques du fonctionnement temps-réel et des points spécifiques d'un fonctionnement hors-ligne dudit système électrochimique qui sont issus d'une transformation de formes d'onde spécifiques extraites de mesures de tension effectuées 5 hors ligne alors que ledit système électrochimique est placé dans des états de fonctionnement connus incluant des états de défaut, de façon à extraire des informations sur l'état de fonctionnement temps-réel dudit système électrochimique. 10 11. Système de diagnostic selon la revendication 10, caractérisé en ce qu'il comprend en outre des moyens pour transformer les formes d'onde spécifiques extraites des mesures de tension effectuées hors-ligne et/ou en temps réel, de façon à créer des points spécifiques. 15 12.Système de diagnostic selon la revendication 11, caractérisé en ce qu'il comprend en outre des moyens pour identifier les points spécifiques dans un espace bidimensionnel (2D) et des moyens pour classer lesdits points spécifiques ainsi identifiés dans un espace tridimensionnel (3D). 20 13.Système de diagnostic selon l'une quelconque des revendications 10 à 12, caractérisé en ce qu'il comprend en outre une base de données par apprentissage sur des états de fonctionnement connus dudit système électrochimique. 25 14.Système électrochimique comportant un empilement de cellules, intégrant un système de diagnostic selon l'une quelconque des revendications 10 à 13. 15.Système électrochimique selon la revendication 14, caractérisé en ce qu'il 30 comporte une pile à combustible.