Le système de collaboration cyber-physique intégrant une IA explicable a été testé expérimentalement sur environ 12 000 enregistrements audio provenant du jeu de données MIMII, qui contient des conditions de fonctionnement normales et anormales de machines. Le jeu de données d'apprentissage comprenait environ 9 600 enregistrements, tandis que le jeu de données de test comprenait environ 2 400 échantillons audio. Selon les résultats, le modèle CNN-transformeur proposé a atteint une précision de classification de 98,5 % sur le jeu de données de test, entre 0,97 et 0,98. Ces résultats correspondent à ceux du Tableau 3, de la matrice de confusion et de l'analyse de stabilité statistique. L'intégration de l'encodeur transformeur améliore les capacités de modélisation temporelle de l'architecture, tandis que le CNN extrait des caractéristiques spatiales discriminantes à partir du spectrogramme d'entrée.
Outre les performances de classification, d'autres caractéristiques liées au système ont été évaluées. Premièrement, la latence d'inférence pour une implémentation sur appareil embarqué n'a pas dépassé 20 ms, permettant une détection en temps réel des anomalies. Deuxièmement, le système présenté contribue à améliorer l'efficacité énergétique en réduisant les besoins en communication. Enfin, le cadre d'intelligence artificielle explicative proposé fournit une interprétabilité SHAP, facilitant l'analyse du modèle du point de vue d'un expert du domaine. Les données MIMII ont été converties en spectrogrammes à l'aide de la transformation de Fourier à court terme (STFT), en utilisant une taille de fenêtre de 1024, une taille de saut de 512 et une taille de transformation de Fourier rapide (FFT) de 1024. Les log-mel-spectrogrammes obtenus ont été redimensionnés à 224 × 224 pixels, puis normalisés selon la méthode z-score. Des méthodes d'augmentation des données, incluant le masquage temporel, le masquage fréquentiel et l'ajout de bruit gaussien, ont été utilisées afin d'accroître la robustesse du modèle entraîné. L'architecture du CNN comprenait quatre couches de convolution avec respectivement 32, 64, 128 et 256 filtres, ainsi que des couches de normalisation par lot, de ReLU et de max-pooling. Les cartes de caractéristiques extraites ont été transformées en séquences et utilisées comme entrées pour l'encodeur de transformeur, composé de quatre couches d'encodeur, huit têtes d'attention, un vecteur d'incorporation de dimension 256 et une dimension de couche entièrement connectée de 1024. Une probabilité de dropout de 0,3 a été appliquée afin de minimiser le surapprentissage. L'optimiseur Adam, avec un taux d'apprentissage de 0,0001 et une décroissance du poids de 1 × 10⁻5, a été utilisé pour l'entraînement. Le nombre d'époques et la taille du lot lors de l'entraînement étaient respectivement de 100 et 32. La graine aléatoire a été fixée à 42 afin de garantir la reproductibilité de l'expérience. Les méthodes suivantes ont été utilisées pour vérifier la validité statistique des résultats : toutes les expériences ont été exécutées indépendamment 10 fois avec différentes graines aléatoires et différents réarrangements des données. Les valeurs moyennes et les écarts-types de l'exactitude (Accuracy), de la précision (Precision), du rappel (Recall), du score F1, de l'aire sous la courbe caractéristique de fonctionnement du récepteur (ROC-AUC) et de l'aire sous la courbe précision-rappel (PR-AUC) ont été calculés. De plus, des intervalles de confiance à 95 % ont été calculés afin d'estimer l'incertitude des performances. Enfin, un test t apparié a été effectué pour comparer les résultats de notre approche CNN-transformer avec ceux du modèle de référence offrant les meilleures performances. La méthode d'arrêt précoce a été appliquée en utilisant 10 époques comme paramètre de patience. Tableau 4 présente la configuration environnementale et des hyperparamètres, et Tableau 5 fournit les détails relatifs à l'environnement de simulation.
L'expérience d'évaluation menée dans cette étude a été réalisée à l'aide du jeu de données MIMII, qui inclut diverses conditions de fonctionnement des machines et des scénarios acoustiques de défauts. Bien que les résultats du cadre proposé valident le modèle, une validation supplémentaire sur divers jeux de données industriels et dans différents environnements de fabrication est nécessaire.
Le cadre proposé a été mis en œuvre sur une plateforme de calcul en périphérie équipée d'un processeur multicœur, d'une mémoire système de 16 Go et d'un accélérateur périphérique doté d'une unité de traitement graphique (GPU) afin de prendre en charge l'inférence en temps réel et l'explicabilité. Des capteurs acoustiques industriels transmettaient des flux audio à la passerelle périphérique via la couche de communication IIoT. Les nœuds périphériques étaient utilisés pour effectuer des tâches telles que la création de spectrogrammes, l'inférence d'anomalies et la génération d'explications SHAP, tandis que le serveur cloud servait à stocker les données à long terme et à mettre à jour le modèle. Une telle configuration réduit non seulement la quantité de communication requise, mais permet également une interaction en temps réel entre les dispositifs de détection, les modules d'intelligence artificielle et les opérateurs humains dans les systèmes cyber-physiques.
L'approche CNN-transformer présentée ci-dessus est comparée à plusieurs méthodes existantes, notamment des modèles classiques d'apprentissage automatique tels que la machine à vecteurs de support (SVM) et la forêt aléatoire, des modèles CNN sans transformeur, ainsi que des méthodes basées sur la mémoire à long court terme (LSTM). Alors que les modèles traditionnels d'apprentissage automatique n'obtiennent pas de résultats satisfaisants en raison de leur incapacité à traiter les corrélations spatiales et temporelles dans les données, les modèles CNN peuvent gérer efficacement les caractéristiques spatiales, tandis que les modèles LSTM sont capables de traiter les aspects temporels des signaux acoustiques. Néanmoins, les premiers sont inférieurs aux seconds en termes de coût de calcul et de capacité de traitement parallèle.
Afin d'assurer une comparaison équitable entre le modèle CNN-transformeur proposé et les modèles de référence, ce dernier a été entraîné et évalué à l'aide de la même répartition du jeu de données MIMII (80 % d'entraînement et 20 % de test), des mêmes étapes de prétraitement, du même processus de création des spectrogrammes et des mêmes mesures d'évaluation. Dans la méthode SVM, un noyau à fonction de base radiale (RBF) a été utilisé, avec C = 10 et γ = 0,01. La méthode forêt aléatoire comprenait 200 arbres, avec une profondeur maximale de 20. En ce qui concerne le CNN, il comportait quatre couches de convolution (32, 64, 128 et 256 filtres), suivies de couches entièrement connectées sans module Transformeur. Enfin, pour la méthode LSTM, deux couches LSTM empilées ont été utilisées avec 128 unités cachées et un taux de dropout de 0,3.
Au contraire, l'utilisation du cadre CNN-transformer permet d'obtenir les meilleurs résultats en exploitant efficacement les caractéristiques spatiales et temporelles des signaux acoustiques. Par ailleurs, le recours à l'informatique en périphérie permet de réduire les coûts de calcul, d'améliorer la latence et de réaliser des économies d'énergie. En outre, l'explicabilité des prédictions mise en œuvre à l'aide de SHAP peut être considérée comme un autre avantage qui distingue la méthode proposée des techniques existantes.
Malgré les récents progrès réalisés dans le développement de nouveaux modèles pour la détection d'anomalies acoustiques, tels que le Vision Transformer (ViT), l'Audio Spectrogram Transformer (AST), le Conformer et les modèles basés sur l'apprentissage auto-supervisé, une comparaison expérimentale de ces architectures n'a pas été menée dans la présente étude. À l'avenir, ces modèles devraient être utilisés comme références pour évaluer le cadre proposé.
Les métriques, notamment la précision, l'exactitude, le rappel et le score F1, sont utilisées pour évaluer les performances du modèle proposé. Alors que la précision mesure le rapport entre les vrais positifs et le nombre total de vrais positifs et de faux positifs, l'exactitude évalue la justesse globale des prédictions. Le rappel fournit une estimation de la capacité du modèle à prévoir les anomalies dans un système donné. Ceci est crucial, car l'incapacité à prévoir un défaut pourrait entraîner des pannes du système. Voici comment ces métriques sont calculées mathématiquement :
(25)
(26)
(27)
(28)
Outre ceux-ci, des métriques supplémentaires liées aux performances du système, telles que la latence, l'utilisation de l'énergie et la scalabilité, sont prises en compte lors de l'évaluation des performances dans des conditions réelles du modèle CPS proposé fondé sur l'informatique en périphérie.
Outre les techniques classiques de mesure de performance, il existe de nombreuses autres mesures pouvant aider à évaluer la solution proposée selon une perspective plus globale, particulièrement en raison de son applicabilité à des études de cas réelles portant sur des données déséquilibrées et en flux continu. Par exemple, la capacité de l'algorithme à distinguer entre les deux groupes selon différents seuils est mesurée par l'aire sous la courbe caractéristique de fonctionnement du récepteur (ROC-AUC). Étant donné sa robustesse et sa séparabilité, la valeur de ROC-AUC devrait être élevée. De plus, la métrique aire sous la courbe précision-rappel (PR-AUC) est importante pour l'étude de cas de détection d'anomalies, car les classes anormales sont rares.
Une mesure importante à prendre en compte est le coefficient de corrélation de Matthews (MCC). Il tient compte des quatre éléments de la matrice de confusion et fournit une évaluation précise, même en présence d'un déséquilibre entre les classes. Cette mesure est calculée comme suit :
(29)
Les valeurs du MCC varient de -1 à +1, une valeur de +1 indiquant une prédiction parfaite. L'autre métrique qui pourrait être appliquée est la Spécificité (taux de vrais négatifs), qui serait utile dans une application industrielle pour évaluer s'il existe des prédictions faussement positives :
(30)
Afin de s'assurer qu'il n'existait aucune incohérence entre la matrice de confusion et les métriques de performance, tous les critères d'évaluation ont été calculés à partir du jeu de données de test final. Les notations vrais positifs (VP), vrais négatifs (VN), faux positifs (FP) et faux négatifs (FN) sont utilisées dans les équations suivantes. L'exactitude, la précision, le rappel, la spécificité, le score F1 et le MCC ont été calculés selon les équations (25) à (30). De plus, la ROC-AUC et la PR-AUC ont été déterminées à partir des sorties de probabilité du modèle CNN-transformer par évaluation indépendante des seuils.
Les métriques clés à prendre en compte au niveau du système incluent la latence, le débit et le coût de calcul. La latence correspond au temps nécessaire pour effectuer une prédiction. Le débit représente la quantité d'échantillons de données traités par seconde. Le coût de calcul, qui reflète l'efficacité énergétique, est important pour garantir une efficacité informatique. En outre, des métriques d'explicabilité, telles que la cohérence SHAP et la stabilité de l'importance des caractéristiques, peuvent également être prises en compte lors de l'évaluation du modèle.
Il ressort clairement des résultats du Tableau 3 pour le jeu de données MIMII que le modèle proposé combinant CNN et transformeur surpasse toutes les autres méthodes d'apprentissage automatique et d'apprentissage profond. Les performances des modèles SVM et forêts aléatoires sont acceptables, bien que limitées, avec des précisions respectives de 88,6 % et 91,2 %. Ceci pourrait s'expliquer par l'incapacité de ces modèles à extraire des caractéristiques temporelles et spatiales complexes à partir des signaux acoustiques. Toutefois, le LSTM surpasse les modèles SVM et forêts aléatoires, atteignant une précision de 94,5 % en modélisant les données de séries temporelles. Néanmoins, l'approche proposée combinant CNN et transformeur présente les meilleures performances, avec une précision améliorée (98,5 %), une justesse (98,06 %), un rappel (99,02 %) et un score-F1 (98,54 %) élevés. Ainsi, l'algorithme proposé permet de classer avec précision les états normaux et anormaux des machines. La spécificité améliorée (97,96 %) démontre l'efficacité de cette approche pour réduire les faux positifs, ce qui est crucial pour l'industrie. Il convient de noter que la performance du CNN-transformeur s'explique par la combinaison du CNN et du transformeur, le premier extrayant des caractéristiques discriminantes et le second capturant les relations temporelles dans de longues séquences.
La courbe ROC de la Figure 6 montre la performance comparative des approches SVM, forêt aléatoire, LSTM et CNN-transformeur proposée en ce qui concerne leurs capacités de classification sur le jeu de données MIMII. Comme indiqué dans la Figure 6, la plus grande aire sous la courbe (AUC) a été observée pour le CNN-transformeur proposé, qui atteint 0,994. Il est clair que la méthode proposée offre une meilleure performance discriminante pour distinguer les catégories normales des catégories anormales. En revanche, l’LSTM donne des résultats légèrement inférieurs à ceux de l’approche CNN-transformeur proposée, son AUC étant d’environ 0,97. La forêt aléatoire et le SVM ont tous deux des AUC relativement faibles, d’environ 0,958 et 0,913 respectivement. Ceci s’attribue à l’incapacité des deux approches à apprendre des motifs acoustiques complexes à partir des données fournies.
Dans la détection d'anomalies, où les échantillons anormaux sont peu nombreux, la courbe précision-rappel (PR) permet une analyse plus pertinente. Sur la courbe PR, le modèle CNN-transformeur atteint la plus haute précision moyenne (AP), d’environ 0,97 à 0,98, ce qui indique une capacité supérieure à détecter les instances anormales avec moins de fausses alarmes, comme illustré dans la Figure 7. En revanche, le classifieur LSTM arrive en deuxième position, avec une AP comprise entre 0,92 et 0,94. D'autre part, les classifieurs forêt aléatoire et SVM se classent respectivement troisième et quatrième, avec des AP de 0,88 et 0,84. Le classifieur hybride proposé atteint une précision plus élevée à tous les niveaux de rappel, ce qui est important pour les applications industrielles pratiques, car celles-ci exigent à la fois une grande précision et un faible taux de fausses alarmes.
Le modèle CNN-transformer a été comparé aux modèles existants sur le jeu de données MIMII à l'aide des métriques de performance avancées présentées dans le Tableau 6. La valeur ROC-AUC indique la capacité discriminante du modèle entre les deux classes à différents seuils. Les classificateurs SVM et forêt aléatoire atteignent des valeurs de ROC-AUC de 0,913 et 0,958 respectivement, tandis que le modèle LSTM obtient un ROC-AUC de 0,970, et le CNN obtient 0,98. Le cadre CNN-transformer proposé atteint le ROC-AUC le plus élevé, de 0,994, démontrant une discrimination supérieure entre les conditions normales et anormales des machines. De même, les valeurs PR-AUC indiquent que le modèle proposé peut gérer efficacement le déséquilibre des classes, un facteur essentiel dans la détection d'anomalies. Le CNN-transformer proposé atteint le PR-AUC (précision moyenne) le plus élevé, de 0,982, suivi par le CNN (0,950), LSTM (0,912), forêt aléatoire (0,891) et SVM (0,765), ce qui témoigne d'une performance supérieure en termes de précision et de rappel pour la détection d'anomalies acoustiques industrielles. En outre, le coefficient de corrélation de Matthews (MCC), qui mesure l'équilibre de la performance du modèle, présente une valeur relativement élevée de 0,97 pour le modèle proposé. En revanche, les méthodes antérieures telles que le SVM (0,73) et la forêt aléatoire (0,78) présentaient une précision prédictive inférieure.
Afin de tester la stabilité de l'approche présentée, la même expérience a été réalisée 10 fois avec des graines d'initialisation et des schémas de mélange des données variables. Les résultats du tableau 7 pour l'architecture CNN-transformer étaient les suivants : exactitude = 98,50 % ± 0,19 %, précision = 98,06 % ± 0,23 %, rappel = 99,02 % ± 0,22 % et score F1 = 98,54 % ± 0,24 %. L'écart-type relativement faible indique la stabilité du modèle. Des intervalles de confiance à 95 % ont été calculés pour toutes les mesures d'évaluation. Il a été déterminé que l'intervalle de confiance pour la métrique d'exactitude était [98,1 % ; 98,9 %], ce qui indique que le modèle présente des performances constantes. Des intervalles de confiance étroits ont également été obtenus pour la ROC-AUC, la PR-AUC et le MCC. Le test t apparié comparant le modèle CNN-transformer proposé et le modèle de référence ayant les meilleures performances (LSTM) a donné une valeur p < 0,05, indiquant que le gain de performance observé est significatif et non dû au hasard.
Afin de tester la robustesse de l'architecture CNN-transformeur proposée, les expériences ont été répétées 10 fois avec différentes graines d'initialisation aléatoire et configurations de mélange, comme illustré dans la Figure 8. La précision moyenne, l'exactitude, le rappel et le score F1 obtenus par le modèle étaient respectivement de 98,50 % ± 0,19 %, 98,06 % ± 0,23 %, 99,02 % ± 0,22 % et 98,54 % ± 0,24 %. D'après l'analyse des boîtes à moustaches, les variations des métriques de performance du modèle sont minimes, ce qui indique sa stabilité et sa robustesse. Ainsi, l'approche proposée présente une forte stabilité au cours de plusieurs expériences, sans variation significative, ce qui témoigne des bonnes capacités de généralisation du modèle.
De plus, un test statistique a été effectué en se basant sur les résultats de plusieurs exécutions de l'expérience. Les faibles écarts types observés pour toutes les mesures d'évaluation suggèrent une faible sensibilité aux variations aléatoires ou à l'apprentissage. Cela confirme que le modèle CNN-transformer proposé peut offrir des performances stables et reproductibles dans des systèmes cyber-physiques de fabrication en conditions réelles.
La carte thermique SHAP donne une indication sur l'importance des régions temps-fréquence dans le spectrogramme pour prédire si l'état de la machine est normal ou anormal. Les valeurs SHAP indiquent la contribution de chaque prédicteur à la prédiction finale.
Dans la carte thermique des valeurs SHAP présentée dans la Figure 9A, les valeurs SHAP positives (par exemple, comprises entre +0,6 et +1,2) représentent une contribution en faveur de l’anormalité, tandis que les valeurs SHAP négatives (c’est-à-dire comprises entre −0,3 et −0,8) indiquent une contribution en faveur de la catégorie normale. Le score SHAP le plus élevé a été observé dans la gamme de fréquences moyennes (40 à 80 intervalles de fréquence) aux instants 50 à 100. Ces régions constituent les zones clés où les signaux de défaut critiques sont détectés dans les mesures acoustiques. De plus, la bande de basses fréquences (moins de 20 intervalles) présente des scores SHAP presque nuls (−0,1 à +0,1), indiquant un impact négligeable sur la prise de décision du modèle. Cela montre que le modèle ignore le bruit de fond non pertinent pour la prédiction des défauts. La cohérence temporelle des scores SHAP élevés au sein de bandes temporelles spécifiques met également en évidence la capacité du module Transformeur à capturer les dépendances à longue portée dans les signaux acoustiques. En résumé, la carte thermique SHAP montre clairement que le modèle CNN-Transformeur se concentre sur des bandes temps-fréquence spécifiques possédant un fort pouvoir discriminant. Cela améliore l’interprétabilité du modèle et aide l’opérateur humain à détecter visuellement les défauts.
On montre que la méthode SHAP est utilisée uniquement comme une technique d'interprétabilité et ne fait pas partie du processus d'entraînement du modèle CNN-transformeur. La précision de classification est déterminée uniquement par les paramètres appris par les composants CNN et Transformeur. SHAP est utilisée après l'étape d'inférence du modèle pour interpréter les prédictions en identifiant les fenêtres temps-fréquence les plus influentes qui orientent les prédictions vers les catégories normales ou anormales. Ainsi, SHAP améliore la transparence et la compréhension humaines sans affecter la précision de classification. Les valeurs SHAP fournies par le cadre proposé permettent une visualisation claire qui aidera les opérateurs de maintenance à vérifier les prédictions du modèle en comparant les régions détectées avec le comportement réel de la machine et les historiques de maintenance. Par conséquent, il sera plus facile pour les humains de collaborer avec les machines lors d'activités de maintenance prédictive.
Afin d'évaluer le niveau d'explicabilité, une analyse globale SHAP sur l'importance des caractéristiques a été réalisée, comme illustré dans la Figure 9B. Les résultats montrent que certaines composantes de fréquence moyenne jouent un rôle important dans la détection d'anomalies. De plus, une étude de cas illustrée dans la Figure 9C montre que SHAP identifie correctement les régions pertinentes aux défauts dans les spectrogrammes anormaux des machines. Cette étude a impliqué une étude de cas qualitative utilisant des échantillons de données anormales provenant de pompes et de vannes du jeu de données MIMII. En ce qui concerne les pompes défectueuses, SHAP a identifié des plages de fréquence de 2–4 kHz associées à la cavitation et à l'usure. Dans le cas des vannes défectueuses, l'activation de SHAP était dominante dans les zones présentant des signaux harmoniques répétés, indicatifs de fuites. Sur des échantillons tests représentatifs, les visualisations SHAP ont identifié de manière fiable des zones temporelles et fréquentielles discriminantes associées à des états anormaux des machines. En mettant en évidence les régions spectrales qui contribuent le plus significativement à chaque prédiction, ces explications permettent de mieux comprendre comment le modèle prend ses décisions. L'analyse SHAP est proposée comme un outil d'interprétabilité permettant de comprendre le comportement du modèle, et non comme une validation experte des défauts, car une validation formelle par des professionnels de la maintenance sortait du cadre de cette étude.
La matrice de confusion met en évidence la performance du modèle CNN-transformer proposé pour la classification sur le jeu de données MIMII illustré dans la Figure 10. Au total, 960 échantillons normaux ont été correctement identifiés comme normaux (vrais négatifs), tandis que 20 échantillons normaux ont été incorrectement identifiés comme anormaux (faux positifs). En outre, parmi tous les échantillons, 1010 échantillons anormaux ont été détectés comme tels (vrais positifs), alors que 10 échantillons anormaux ont été identifiés comme normaux (faux négatifs).
Ainsi, une performance élevée de détection a été obtenue, en particulier pour les échantillons anormaux, qui sont cruciaux dans la maintenance prédictive. Comme indiqué dans la Figure 10, une erreur de faux négatif est observée très rarement (10 échantillons) – ce qui est essentiel pour garantir la sécurité et la fiabilité dans tout secteur industriel. En revanche, l'erreur de faux positif (20 échantillons) se produit légèrement plus fréquemment, mais reste à un niveau relativement faible. La matrice de confusion présentée dans la Figure 10 a été générée à l'aide d'un sous-ensemble équilibré de 2 000 enregistrements de test (nombres égaux d'échantillons normaux et anormaux) afin de fournir une visualisation claire des performances du classificateur. Toutes les métriques d'évaluation rapportées, y compris l'exactitude, la précision, le rappel, la spécificité et le score F1, ont été calculées à l'aide de l'ensemble complet de test (environ 2 400 enregistrements).
En tenant compte de la matrice de confusion illustrée dans la Figure 10 et dans le Tableau 8 (VN = 960, FP = 20, VP = 1010, FN = 10), les calculs suivants peuvent être effectués : Précision globale = (VP + VN) / (VP + VN + FP + FN) = 98,50 %, Précision = VP / (VP + FP) = 98,06 %, Rappel = VP / (VP + FN) = 99,02 %, Spécificité = VN / (VN + FP) = 97,96 %, Score F1 = 98,54 %, et MCC = 0,97. De plus, les valeurs de ROC-AUC et de PR-AUC étaient respectivement de 0,994 et de 0,982. La matrice de confusion présentée dans la Figure 10 a été générée à l’aide d’un sous-ensemble équilibré et représentatif de 2 000 enregistrements provenant des données de test, uniquement à des fins de visualisation. Toutes les métriques quantitatives de performance rapportées dans cette étude ont été calculées en utilisant l’ensemble complet des données de test (environ 2 400 enregistrements). Dans l’ensemble, les résultats démontrent que l’architecture CNN-transformeur proposée atteint une précision de classification de 98,5 %, en accord avec les résultats indiqués dans le Tableau 8 et l’analyse d’ablation.
Afin de vérifier l'efficacité de l'implémentation basée sur les nœuds périphériques, des mesures supplémentaires ont été effectuées. La latence moyenne d'inférence du modèle proposé est de 18,7 ms/échantillon, et il peut traiter environ 53 échantillons par seconde, comme indiqué dans Tableau 9. La consommation d'énergie pendant l'inférence était de 8,9 W, ce qui correspond à une consommation énergétique estimée de 0,167 J/prédiction. Pour évaluer l'impact environnemental de l'approche proposée dans Tableau 10, la surcharge de communication, la consommation d'énergie et l'utilisation des ressources ont été prises en compte. Étant donné que le processus d'inférence s'effectue localement sur les nœuds périphériques, seules les données diagnostiques résumées sont transmises vers l'environnement cloud. Les résultats expérimentaux ont montré une réduction de 73 % de la communication par rapport aux systèmes CPS traditionnels basés sur le cloud. De plus, l'inférence locale a permis de réduire la consommation d'énergie de 32 %, tandis que le temps de communication est passé de 75 ms à 20 ms.
Étude d'ablation
Le jeu de données MIMII est soumis à une étude d'ablation afin d'évaluer l'impact de chaque composant de l'architecture CPHS basée sur l'intelligence artificielle explicable. Cette étude d'ablation examine les effets de l'utilisation d'un CNN pour apprendre les caractéristiques spatiales, d'un encodeur de type transformeur pour apprendre les caractéristiques temporelles, ainsi que de l'algorithme d'explicabilité SHAP lui-même. Plusieurs variantes du modèle proposé sont envisagées, à savoir : (i) une version CNN uniquement, (ii) une version transformeur uniquement, (iii) le modèle combinant CNN et transformeur mais sans l'algorithme SHAP, et (iv) le modèle complet proposé combinant CNN, transformeur et SHAP. Pour l'évaluation, des métriques classiques telles que la précision, le rappel, la justesse et le score F1 sont utilisées.
Ces résultats d'étude d'ablation démontrent l'importance de chaque élément du cadre décrit dans le Tableau 11. Le modèle reposant uniquement sur un CNN atteint une précision de 93,8 %, ce qui indique que l'extraction de caractéristiques spatiales à partir des spectrogrammes est efficace, bien qu'il ne puisse pas capturer les dépendances temporelles présentes dans les signaux acoustiques. Il en va de même pour le modèle utilisant uniquement un transformeur, qui atteint une précision légèrement inférieure (92,6 %), car il extrait des caractéristiques temporelles tout en négligeant l'information spatiale. Enfin, le modèle CNN-transformeur proposé atteint une précision de 98,5 %, une exactitude de 98,06 %, un rappel de 99,02 % et un score F1 de 98,54 %, en accord avec les résultats de l'expérience principale. Étant donné que SHAP est utilisé uniquement pour l'interprétation a posteriori et non pour l'apprentissage du modèle, la précision de classification est identique à celle du modèle CNN-transformeur entraîné. SHAP est ensuite utilisé pour produire des explications fondées sur l'attribution des caractéristiques. Le rappel amélioré (99,02 %) indique que le modèle est efficace pour détecter les conditions anormales des machines, garantissant ainsi qu'aucune pièce défectueuse n'est omise, ce qui est essentiel pour la mise en œuvre d'un système de maintenance prédictive.
DISPONIBILITÉ DES DONNÉES :
Le jeu de données MIMII (Malfunctioning Industrial Machine Investigation and Inspection) utilisé dans cette étude est publiquement disponible sur le dépôt officiel Zenodo. Les expériences ont été menées à l'aide des enregistrements audio et des annotations publiquement accessibles fournis par les auteurs du jeu de données. Le jeu de données peut être consulté à l'adresse suivante : https://zenodo.org/records/3384388. Le code de mise en œuvre accompagnant cette étude est publiquement disponible via le dépôt Zenodo à l'adresse https://zenodo.org/records/21405292. Le dépôt inclut le code de mise en œuvre du cadre proposé, comprenant le pipeline de prétraitement des données, l'architecture du modèle, le flux d'entraînement, les scripts d'évaluation, les fichiers de configuration et les utilitaires complémentaires afin de faciliter la compréhension et la reproduction indépendante de la méthodologie proposée. Le jeu de données MIMII n'est pas redistribué avec le code de mise en œuvre et doit être obtenu séparément depuis son dépôt officiel.

Figure 1 : Architecture du travail proposé. Conception architecturale d'un système cyber-physique humain piloté par une IA explicative combinant un prétraitement, un apprentissage spatial basé sur les CNN, une modélisation temporelle basée sur les transformeurs, une explicabilité SHAP et un calcul en périphérie pour une fabrication durable. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 2 : Pipeline de prétraitement. Pipeline de prétraitement du signal acoustique comprenant la suppression du bruit, la normalisation, la segmentation, la transformation STFT et l'extraction du spectrogramme. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 3 : Extraction de caractéristiques spatiales à l'aide de CNN. Architecture permettant l'extraction de caractéristiques spatiales par des réseaux de neurones convolutifs (CNN) utilisant des spectrogrammes, au moyen des étapes de convolution, d'activation et de sous-échantillonnage. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 4 : Modélisation temporelle à l'aide d'un encodeur de type transformeur. Architecture d'encodeur de type transformeur permettant de capturer les dépendances temporelles et les dépendances à longue portée au moyen d'une représentation sous forme de séquence de spectrogramme. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 5 : Couche de classification. Une couche de classification qui associe les caractéristiques temporelles encodées à des valeurs de probabilité pour prédire l'état de la machine. Veuillez cliquer ici pour afficher une version agrandie de cette figure.

Figure 6 : Résultat de la courbe ROC. Comparaison des courbes ROC pour les algorithmes SVM, forêt aléatoire, LSTM, CNN et CNN-transformeur proposé sur le jeu de données MIMII. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 7 : Résultat de la courbe précision-rappel. Comparaison des courbes précision-rappel illustrant l'efficacité de détection des anomalies de diverses méthodes d'apprentissage automatique et d'apprentissage profond. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 8 : Analyse de la stabilité des performances sur plusieurs exécutions. Diagramme en boîte montrant la stabilité des performances en termes de précision, exactitude, rappel et score F1 obtenues lors de dix expériences indépendantes menées avec l'architecture CNN-transformeur proposée. Les faibles variations indiquent la stabilité et la robustesse du modèle. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 9 : Analyse explicative basée sur SHAP du cadre proposé pour la détection d'anomalies. (A) Analyse explicative basée sur SHAP. Carte thermique utilisant la méthode SHAP pour montrer les zones temporelles et fréquentielles les plus pertinentes menant à l'identification d'un comportement anormal. (B) Analyse de l'importance des caractéristiques à l'aide de l'approche globale SHAP pour identifier les caractéristiques les plus influentes. (C) Visualisation de l'explicabilité SHAP pour une machine anormale. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 10 : Matrice de confusion pour le jeu de données MIMII. La matrice de confusion illustre les performances de classification de l'algorithme CNN-transformer proposé. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
| Paramètre | Description |
| Nom du jeu de données | Jeu de données MIMII (Malfunctioning Industrial Machine Investigation and Inspection) |
| Types de machines | Vannes, pompes, ventilateurs, rails coulissants |
| Nombre total d'échantillons utilisés | Environ 12 000 enregistrements audio |
| Classes | Normale, Anormale |
| Fréquence d'échantillonnage | 16 kHz |
| Durée par échantillon | 10 secondes |
| Format audio | WAV |
| Représentation des caractéristiques | Images de spectrogrammes basées sur la TFCT |
| Étapes de prétraitement | Filtrage du bruit, normalisation, segmentation, transformation par TFCT |
| Attribution des étiquettes | Les enregistrements normaux sont étiquetés comme la classe 0 ; les enregistrements anormaux sont étiquetés comme la classe 1 selon les annotations MIMII |
| Répartition entraînement-test | 80 % pour l'entraînement (≈9 600 échantillons), 20 % pour le test (≈2 400 échantillons) |
| Stratégie de partitionnement des données | Classement au niveau des fichiers audio effectué avant la segmentation afin d'éviter les fuites de données |
| Conditions de bruit industriel | Environnements acoustiques industriels réalistes inclus dans les enregistrements MIMII |
| Domaine d'application | Maintenance prédictive et détection d'anomalies des machines |
| Objectif ciblé de l'industrie 5.0 | Détection de défauts explicative, centrée sur l'humain et durable |
Tableau 1 : Description du jeu de données. Description du jeu de données MIMII en fonction du type de machine, de la distribution des échantillons, des caractéristiques du signal, de la représentation des caractéristiques et des aspects applicatifs.
| Type de machine | Identifiants de machine | Enregistrements normaux | Enregistrements anormaux | Conditions de rapport signal/bruit (dB) | Répartition expérimentale |
| Ventilateur | ID00–ID06 | 8 400 | 1 600 | −12 à −9 | 80 % apprentissage / 20 % test |
| Boîte de vitesses | ID00–ID06 | 7 124 | 1 147 | −17 à −15 | 80 % apprentissage / 20 % test |
| Pompe | ID00–ID06 | 7 200 | 1 800 | −18 à −9 | 80 % apprentissage / 20 % test |
| Glissière | ID00–ID05 | 7 000 | 1 800 | −14 à −12 | 80 % apprentissage / 20 % test |
| Vanne | ID00–ID05 | 7 000 | 1 800 | −12 à −9 | 80 % apprentissage / 20 % test |
Tableau 2 : Description du jeu de données MIMII utilisé dans cette étude. Types de machines industrielles, identifiants des machines, nombres d'enregistrements normaux et anormaux, conditions de rapport signal sur bruit (SNR), et la répartition entraînement/test utilisée pour la classification binaire supervisée.
| Modèle | Précision (%) | Exactitude (%) | Rappel (%) | Score F1 (%) | Spécificité (%) |
| SVM | 88.6 | 86.9 | 85.4 | 86.1 | 90.2 |
| Forêt aléatoire (RF) | 91.2 | 89.8 | 88.5 | 89.1 | 92.6 |
| LSTM | 94.5 | 93.2 | 92.8 | 93 | 95.1 |
| CNN | 96.3 | 95.4 | 94.9 | 95.1 | 96.8 |
| CNN–Transformeur (proposé) | 98.5 | 98,06 | 99,02 | 98,54 | 97,96 |
Tableau 3 : Résultat de classification pour le jeu de données MIMII. Comparaison des performances des classificateurs d'apprentissage automatique et d'apprentissage profond à l'aide de l'exactitude, de la précision, du rappel, du score F1 et de la spécificité.
| Paramètre | Valeur |
| Jeu de données | MIMII Dataset |
| Nombre total d'échantillons | ~12 000 |
| Répartition entraînement-test | 80 % / 20 % |
| Type d'entrée | Images de spectrogrammes |
| Couches CNN | 3 à 5 couches de convolution |
| Couches de transformeur | 2 à 4 couches d'encodeur |
| Taille du lot | 32 |
| Taux d'apprentissage | 0,001 |
| Optimiseur | Adam |
| Époques | 100 |
| Matériel | Dispositif embarqué + GPU (pour l'entraînement) |
| Longueur de fenêtre STFT | 1024 |
| Taille du pas | 512 |
| Taille de la FFT | 1024 |
| Résolution du spectrogramme | 224 × 224 |
| Normalisation | Z-score |
| Augmentation des données | Masquage temporel, masquage fréquentiel, bruit gaussien |
| Couches CNN | 4 |
| Filtres CNN | 32, 64, 128, 256 |
| Couches de transformeur | 4 |
| Têtes d'attention | 8 |
| Dimension de l'incorporation | 256 |
| Dimension du réseau feed-forward | 1024 |
| Dropout | 0,3 |
| Graine aléatoire | 42 |
| Patience pour l'arrêt précoce | 10 |
Tableau 4 : Configuration de l'environnement et des hyperparamètres. Paramètres d'apprentissage et d'expériences utilisés pour mettre en œuvre l'architecture CNN-transformeur proposée.
| Composant | Spécification |
| Langage de programmation | Python |
| Framework | TensorFlow / PyTorch |
| Processeur | Intel i7 / Ryzen 7 |
| GPU | NVIDIA GTX 1660 / RTX Series |
| Dispositif embarqué | Raspberry Pi / NVIDIA Jetson |
| Visualisation | Matplotlib, SHAP |
| Dispositif embarqué | NVIDIA Jetson Xavier NX |
| UC | 6 cœurs ARM v8.2 |
| RAM | 16 Go |
| Système d'exploitation | Ubuntu 20.04 |
| Framework d'apprentissage profond | TensorFlow/PyTorch |
| Taille du jeu de données | 12 000 échantillons MIMII |
| Connectivité | Ethernet/Wi-Fi |
| Paramètres du modèle | 8,4 millions |
| Taille du modèle | 32,7 Mo |
| Framework | PyTorch + TensorRT |
| Taille du lot | 1 |
| Quantification | FP16 |
| Élagage | Non |
| Latence moyenne | 17,8 ms |
| Latence médiane | 17,2 ms |
| Latence au 95e centile | 19,6 ms |
| Débit | 56 échantillons/s |
| Utilisation de la mémoire | 1,4 Go |
| Consommation électrique | 11,3 W |
Tableau 5 : Environnement de simulation. Du matériel et des logiciels ont été utilisés pour entraîner, déployer et évaluer le modèle proposé.
| Modèle | ROC-AUC | PR-AUC | MCC |
| SVM | 0.913 | 0.765 | 0.73 |
| Forêt aléatoire (RF) | 0.958 | 0.891 | 0.78 |
| LSTM | 0.97 | 0.912 | 0.86 |
| CNN | 0.98 | 0.95 | 0.9 |
| CNN–Transformeur (proposé) | 0.994 | 0.982 | 0.97 |
Tableau 6 : Comparaison des résultats expérimentaux pour ROC-AUC, PR-AUC et MCC. Comparaison entre différents modèles selon les mesures ROC-AUC, PR-AUC et coefficient de corrélation de Matthews.
| Métrique | Moyenne (%) + Écart type (%) | Intervalle de confiance à 95 % |
| Exactitude | 98.50 ± 0.19 | [98.1, 98.9] |
| Précision | 98.06 ± 0.23 | [98.0, 98.2] |
| Rappel | 99.22 ± 0.22 | [98.8, 99.4] |
| Score-F1 | 98.54 ± 0.24 | [98.1, 98.9] |
| ROC-AUC | 0.984 | 0.004 |
| PR-AUC | 0.982 | 0.005 |
| MCC | 0.97 | 0.007 |
Tableau 7 : Résultat de l'analyse statistique sur plusieurs exécutions. Statistiques de robustesse de l'architecture CNN-transformer proposée sur 10 expériences, incluant la moyenne, l'écart type, l'intervalle de confiance à 95 % et la significativité pour diverses mesures de performance.
| Métrique | Valeur |
| Exactitude | 98,50 % |
| Précision | 98,06 % |
| Rappel | 99,02 % |
| Spécificité | 97,96 % |
| Score-F1 | 98,54 % |
| MCC | 0,97 |
| ROC-AUC | 0,994 |
| PR-AUC | 0,982 |
Tableau 8 : Ensemble final de test de la matrice de confusion. Matrice de confusion sur l'ensemble de test pour l'architecture proposée, montrant la classification des états normaux et anormaux des machines à partir desquels les métriques d'évaluation ont été déduites.
| Métrique | Déploiement dans le cloud | Déploiement en périphérie (edge) |
| Latence d'inférence (ms) | 85,6 | 18,7 |
| Débit (échantillons/s) | 22 | 53 |
| Utilisation du réseau (Mo/min) | 120 | 18 |
| Énergie par prédiction (J) | 0,52 | 0,17 |
| Capacité en temps réel | Modérée | Élevée |
Tableau 9 : Résultat de la latence d'inférence. Les performances en termes de latence du déploiement en périphérie du cadre CPHS doté d'une IA explicable, incluant le temps de traitement, le débit, la consommation mémoire et d'autres caractéristiques en temps réel.
| Métrique | CPS basé sur le cloud | CPHS proposé avec fonctionnalités Edge | Amélioration |
| Transmission de données par heure | 100% | 35% | Réduction de 65 % |
| Latence de communication | 75 ms | 20 ms | Réduction de 73,3 % |
| Consommation d'énergie | 100% | 68% | Réduction de 32 % |
| Émission de carbone estimée | 100% | 70% | Réduction de 30 % |
Tableau 10 : Évaluation de la durabilité du CPHS doté de capacités de périphérie (edge). Des mesures de la surcharge de communication, de la consommation d'énergie, de l'utilisation des ressources et de la latence sont utilisées pour évaluer la durabilité du CPHS doté de capacités de périphérie.
| Variante du modèle | Précision (%) | Exactitude (%) | Rappel (%) | Score F1 (%) |
| CNN uniquement | 93.8 | 92.5 | 94.2 | 93.3 |
| Transformeur uniquement | 92.6 | 91.2 | 93.5 | 92.3 |
| CNN + Transformeur | 98.5 | 98.06 | 99.02 | 98.54 |
Tableau 11 : Résultat de l'étude d'ablation. Résultat de l'étude d'ablation mettant en évidence l'effet des modules CNN, Transformeur et SHAP sur les performances du modèle proposé.