$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Cette section décrit l’ensemble du pipeline d’apprentissage profond contraint en physique pour l’évaluation des conditions thermiques à l’aide d’images thermographiques infrarouges. Le cadre intègre (i) l’extraction du champ de température, (ii) la régression CNN de base, (iii) la régularisation de la physique basée sur la diffusion, (iv) l’apprentissage adaptatif des paramètres, et (v) la quantification de l’incertitude. Tous les outils utilisés dans cette étude sont décrits dans le Tableau des Matériaux.
Acquisition d’image thermique et construction du champ de température
Des images infrarouges (IR) des équipements de puissance haute tension étaient obtenues à l’aide d’une caméra thermographique calibrée. Chaque image était convertie en une matrice de température spatiale T(x, y) ∈R H×W en cartographiant les intensités des pixels en valeurs physiques de température (°C) à l’aide des paramètres d’étalonnage radiométrique de la caméra. La figure 1A–C montre la visualisation de la température. L’image IR originale est convertie en une carte de température physiquement significative. L’histogramme montre une distribution de température déséquilibrée avec une concentration comprise entre 50 et 90 °C. Les zones de haute température localisées (points chauds) sont clairement visibles. La figure 2 montre une représentation 3D de l’élévation thermique du champ de température. Le graphique d’élévation de surface en 3D révèle des pentes abruptes autour des zones chaudes. Les pics de température approchent les ~120 °C. Le champ thermique est spatialement lisse, sauf près des conducteurs sous tension et des joints de connexion. Ces caractéristiques justifient l’utilisation d’un terme de régularisation basé sur la diffusion pour imposer une lissabilité physique tout en préservant des phénomènes localisés à gradient élevé.

Figure 1 : Visualisation de la température. (A) Image thermographique infrarouge (IR) originale de l’équipement haute tension inspecté montrant le motif thermique observé. (B) Carte de température extraite à partir de l’image thermographique, illustrant la distribution spatiale des températures de surface et des gradients thermiques associés. (C) Histogramme de la distribution de la température montrant la distribution fréquentielle des valeurs de température dans le champ de température extrait. Cette figure montre la conversion d’une image thermographique en une représentation quantitative de la température et résume la distribution de température sous-jacente utilisée pour le développement et l’analyse ultérieurs du modèle. Abréviations ; IR = infrarouge. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 2 : Représentation 3D de l’élévation thermique de la surface du champ de température. Visualisation tridimensionnelle de la distribution de température, où l’élévation de la surface correspond à la magnitude de la température et met en évidence les points chauds thermiques localisés. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Formulation du problème
Soit

désignons l’image thermographique infrarouge RGB, où H et W représentent les dimensions spatiales et C = 3 correspond aux canaux de couleur associés à la palette thermique arc-en-ciel. Bien que l’image infrarouge représente visuellement l’information thermique, ce n’est pas une matrice numérique directe de température car les couleurs affichées dépendent de la cartographie de la palette, de l’échelle, de l’interpolation et des conditions d’imagerie.
Soit

désigner le champ de température de référence au sol correspondant obtenu par étalonnage radiométrique et extraction de température à partir du système d’imagerie infrarouge.
L’objectif d’apprentissage est d’approximer la correspondance

où θ désigne les paramètres entraînables du CNN et T̂ représente le champ de température prédit. Le but de cette cartographie supervisée n’est pas simplement de reproduire la conversion radiométrique directe, mais d’apprendre un modèle de régression de température spatialement cohérent à partir d’images thermographiques rendues. Contrairement à l’extraction directe couleur-température pixel par pixel, le CNN peut exploiter la structure thermique contextuelle à travers les régions voisines, réduire la sensibilité au bruit, préserver les caractéristiques des points chauds et générer des champs de température physiquement plausibles plus lisses sous les contraintes de diffusion imposées. La fonction objectif de base minimise l’erreur de régression pixel par pixel entre les champs de température prédit et calibré :

où Tij et T̂ij désignent respectivement la vérité fondamentale et les températures prédites à la position des pixels (i, j). En conséquence, le cadre doit être interprété comme un modèle de régression thermographique régularisé en physique plutôt que comme un remplacement de l’étalonnage radiométrique direct.
Conditions d’acquisition d’image infrarouge
Des images thermographiques infrarouges ont été prises à l’aide d’une caméra thermique FLIR dans des conditions de fonctionnement extérieures de la substation. Le jeu de données comprenait des images thermographiques RVB avec une résolution spatiale de 640 × 480 pixels, enregistrées à l’aide de la palette de couleurs arc-en-ciel FLIR. Chaque image comportait une barre d’échelle thermique intégrée indiquant la plage de température associée à l’image d’acquisition. L’inspection thermographique couvrait plusieurs catégories d’équipements haute tension, notamment les transformateurs, les bagues, les isolants, les connecteurs et les structures de bornes. Des images ont été capturées dans des conditions normales de fonctionnement pour inclure à la fois le comportement thermique standard et les motifs localisés de points chauds pertinents pour les applications de surveillance de l’état. L’acquisition d’images a été réalisée à partir de points de vue d’inspection pratiques généralement utilisés en thermographie de terrain. L’orientation de la caméra et l’angle de vue variaient selon l’accessibilité des équipements et la disposition des sous-stations. Les images thermiques étaient visualement contrôlées pour garantir que l’équipement cible restait clairement observable et que les zones chaudes étaient distinguables du fond. Le jeu de données était organisé en dossiers étiquetés selon la catégorie d’équipement et l’état thermique. Avant le développement du modèle, les images étaient segmentées et redimensionnées dans un format spatial cohérent pour l’extraction des champs de température et la formation CNN. Les images thermographiques contiennent des distributions thermiques adaptées à la vision par ordinateur et aux études d’apprentissage profond impliquant la détection de points chauds, l’identification de défauts, la localisation d’anomalies et la prévision des champs de température. Dans les présents travaux, les profils thermiques extraits des images infrarouges ont été utilisés pour l’estimation de la température limitée par la physique et l’évaluation thermique consciente de l’incertitude.
Jeu de données thermographiques et partitionnement des données
Le jeu de données utilisé dans cette étude comprenait 174 images thermographiques infrarouges acquises à partir d’équipements de sous-station haute tension utilisant un système d’imagerie thermiqueFLIR 15. Le jeu de données comprenait plusieurs types de composants du système électrique, notamment des transformateurs, des bagues, des isolants, des connecteurs et des structures de bornes haute tension fonctionnant dans des conditions extérieures pratiques. Chaque image thermographique était convertie en champ de température spatiale à l’aide des informations radiométriques fournies par la caméra infrarouge. Le jeu de données contient des variations thermiques allant de conditions normales de fonctionnement à des zones chaudes localisées associées à une contrainte thermique élevée. Les données ont été divisées en sous-ensembles d’entraînement, de validation et de test avec une répartition 70:15:15, donnant respectivement 121, 26 et 27 images. Le partitionnement a été effectué au niveau de l’image car le jeu de données consistait en des scènes thermographiques capturées indépendamment, collectées sous des conditions de fonctionnement et des points de vue variés. Des soins ont été pris pour éviter les quasi-doubles de trames entre sous-ensembles afin de réduire le risque de fuite de données et d’améliorer l’équité de l’évaluation. Toutes les images ont été redimensionnées à une résolution spatiale uniforme avant la formation. La normalisation de la température a été appliquée pour stabiliser l’optimisation tout en préservant les gradients thermiques relatifs et les caractéristiques des points chauds importantes pour l’évaluation de l’état.
Architecture CNN de base
Une architecture de régression encodeur-décodeur entièrement convolutionnelle a été adoptée pour préserver la résolution spatiale lors de la prédiction température-champ à partir d’images thermographiques infrarouges. Toutes les images infrarouges d’entrée étaient redimensionnées à une résolution uniforme de 256 × 256 × 3 avant l’entraînement. La section encodeur se composait de quatre blocs convolutionnels avec une profondeur de caractéristique croissante progressivement. Chaque bloc contenait une couche de convolution 3 × 3 suivie d’une normalisation par lots et d’une activation ReLU. Le sous-échantillonnage spatial a été réalisé à l’aide de convolutions stridées plutôt que de pooling maximal afin de mieux préserver les structures thermiques localisées et les informations sur les points chauds. La progression des canaux encodeurs était la suivante :
32 → 64 → 128 → 256
La section décodeur a reconstruit le champ de température spatiale en utilisant des couches de convolution transposées pour l’upsampling avec la progression inversée des caractéristiques :
256 → 128 → 64 → 32
Chaque bloc décodeur utilisait de manière similaire la convolution transposée, la normalisation par lots et l’activation de ReLU. Une dernière couche de convolution 1 × 1 a généré la carte de température à canal unique prédite :
T̂ ∈ R256×256
Le réseau produit donc un champ de température spatiale en pleine résolution correspondant à l’image d’entrée infrarouge. Pour soutenir la quantification de l’incertitude, des couches de dropout avec une probabilité p = 0,2 ont été insérées près des étapes goulot d’étranglement et décodeur. Pendant l’inférence, le dropout est resté actif pour permettre l’échantillonnage Monte Carlo pour l’estimation prédictive de l’incertitude. Tous les poids convolutionnels étaient initialisés par initialisation He-normale, tandis que les termes de polarisation étaient initialisés à zéro. Le modèle a été entraîné à l’aide de l’optimiseur Adam avec un taux d’apprentissage initial de 1 × 10−3, une taille de lot de 8, et une désintégration en poids de 1 × 10−5. La formation a été réalisée pendant 100 époques en utilisant la désintégration du taux d’apprentissage du cosinus et l’arrêt précoce basés sur la validation RMSE.
La fonction objectif de base est définie comme suit :

où Ldonnées désignent l’erreur quadratique moyenne pixel par pixel entre les champs de température prédit et de vérité du terrain.
La performance a été évaluée en utilisant : R2, RMSE et MAE.
Bien que le CNN de base ait réussi à capturer la structure thermique globale et la distribution des points chauds, il a parfois produit des artefacts non physiques dans des régions à faible gradient et un léger surlissage près de limites thermiques nettes. Ces limitations ont motivé l’introduction de la régularisation de la physique basée sur la diffusion dans la formulation ultérieure du modèle.
Apprentissage limité par la physique via diffusion
Équation de la chaleur régissant
Pour assurer la cohérence thermique dans le champ de température prévu, une hypophysique basée sur la diffusion a été intégrée dans l’équation de la chaleur. En conditions quasi-stables, la température de surface satisfait approximativement

où T désigne la température, α est la diffusivité thermique, et S représente la génération locale de chaleur. Puisque les échantillons thermographiques disponibles sont des instantanés d’inspection uniques plutôt que des séquences résolues dans le temps, la dérivée temporelle ne peut pas être évaluée directement. En conséquence, le modèle suppose une conduction quasi-stationnaire locale et impose

loin des zones de chauffage active.
Le laplacien 2D a été discrétisé à l’aide du pochoir standard aux différences finies :

Le résidu physique était alors calculé sur la carte de température prédite comme

Pour éviter toute ambiguïté aux limites de l’image, le résidu laplacien n’a été évalué que sur des pixels intérieurs valides, où toutes les valeurs voisines requises par le pochoir sont disponibles. Aucun masquage explicite n’a été appliqué aux régions chaudes, car l’objectif du terme physique était de régulariser l’ensemble du champ spatial plutôt que de supprimer les zones thermiquement actives. La structure du point chaud a donc été apprise via le terme des données, tandis que la pénalité laplacienne favorisait une fluidité cohérente en diffusion sur la carte prédite.
La fonction objectif finale combinait la perte de données et la perte physique comme

où λphys contrôle la force de la régularisation physique. Dans l’étude d’ablation, λphys ∈ {0, 10−4, 10−3, 10−2, 10−1} a été examinée afin d’évaluer l’effet de l’augmentation de la résistance physique à la contrainte. Des valeurs modérées ont amélioré la généralisation, tandis qu’une pondération trop forte tendait à des gradients de points chauds trop lisses. Ce comportement est cohérent avec les structures thermiques observées dans les images infrarouges, où des régions physiquement lisses coexistent avec des concentrations localisées à haute température.
Il convient de souligner que la formulation de diffusion en régime permanent utilisée dans ce travail n’a pas pour but de représenter le mécanisme complet de transfert de chaleur régissant la formation des points chauds dans les équipements haute tension. En pratique, la surchauffe localisée résulte de processus internes de production de chaleur tels que le chauffage en joule, la résistance de contact et la dégradation de l’isolation. Parce que le jeu de données disponible consiste en images thermographiques instantanées sans informations temporelles ni mesures de source, la contrainte de diffusion est utilisée comme terme de régularisation spatiale inspiré par la physique. Son rôle est de favoriser des champs de température localement lisses et physiquement plausibles, tout en permettant d’apprendre les caractéristiques du point chaud à partir des données. Par conséquent, la diffusion antérieure doit être interprétée comme une contrainte de fluidité spatiale plutôt que comme un modèle global de la physique thermique sous-jacente.
Diffusivité thermique apprenable et augmentation en termes sources
Pour améliorer l’interprétabilité physique, la diffusivité thermique effective a été traitée comme un paramètre scalaire entraînable. Soit T̂(x, y) le champ de température prévu et α > 0 la diffusivité thermique apprenable. Pour faire respecter la positivité, α était paramétrée par une variable non contrainte β comme

Dans cette formulation, le prior physique a été écrit en termes d’un champ de forçage résiduel effectif R(x, y), où R(x, y) représente l’écart local de la diffusion pure induit par des effets de chauffage non modélisés. La relation adaptative-diffusivité s’exprime donc comme suit

Ici, ∇2T̂(x, y) est le laplacien du champ de température prédit, et R(x, y) est le terme source résiduelle effective déduit implicitement des données thermographiques. Parce que les échantillons disponibles sont des images instantanées plutôt que des séquences thermiques résolues dans le temps, R(x, y) n’est pas supervisé séparément et ne sert que comme une grandeur physique latente lors de l’optimisation. La perte physique correspondante était définie comme suit

Pour l’expérience augmentée par terme source, une tête CNN auxiliaire a été ajoutée pour estimer une carte source spatiale Ŝ(x, y). Dans ce cas, Ŝ(x, y) représente le champ local de production de chaleur prédit, et la relation gouvernante devient

La tête source était attachée au goulot d’étranglement de l’encodeur et consistait en une petite branche convolutionnelle qui mappait des caractéristiques latentes partagées sur une carte de sortie à canal unique Ŝ(x, y). Aucune étiquette de source fondée sur le terrain n’était disponible ; par conséquent, Ŝ(x, y) n’était appris que par l’objectif commun

où Lest la perte de régression de température, λphys contrôle la contrainte de diffusion, et λs régularise la prédiction de la source pour décourager des cartes de sources triviales et de grande intensité. Dans le contexte basé sur les instantanés considéré ici, la formulation du terme source est restée sous-déterminée et n’a pas amélioré la précision de la prédiction.
Abandon de Monte Carlo pour la quantification de l’incertitude
Pour estimer l’incertitude épistémique dans le champ de température prévu, le dropout de Monte Carlo a été appliqué lors de l’inférence. Des couches de dropout à un taux de 0,2 étaient placées après le bloc de goulot d’étranglement et à l’intérieur des étages décodeurs du CNN. Lors des essais, ces couches de dropout ont été maintenues actives, et 30 passes stochastiques vers l’avant ont été effectuées pour chaque image infrarouge. Pour le k-ième passage stochastique, le champ de température prédit est noté par

où X est l’image infrarouge d’entrée et θk représente les paramètres du réseau échantillonnés induits par le dropout. La moyenne prédictive était calculée comme suit

et la variance prédictive pixel par pixel était obtenue comme

avec K = 30. La carte de variance résultante a été utilisée pour identifier les régions spatiales ayant une confiance prédictive plus faible. Comme les étiquettes d’incertitude basées sur le terrain ne sont pas disponibles pour la régression thermographique, l’incertitude a été évaluée à l’aide de cartes de variance pixel par pixel et de statistiques résumées calculées sur les régions chaudes et de fond, plutôt que des métriques basées sur l’étalonnage. L’analyse quantitative a montré que l’écart-type prédictif moyen était de 13,046 °C, tandis que la région du point chaud présentait un écart-type plus élevé de 15,934 °C contre 12,743 °C dans la région de fond. Cela correspond à une augmentation d’environ 25 % de l’incertitude dans les zones chaudes par rapport au fond. Ces résultats indiquent que le modèle exprime une incertitude plus grande près des régions thermiquement complexes, en particulier autour des frontières de points chauds et des gradients de température abrupts, ce qui est cohérent avec la structure observée dans les visualisations thermiques de surface. Le flux de travail est résumé dans les Figures 3 et 4 qui montrent le pipeline global de traitement basé sur CNN.

Figure 3 : Flux de travail. Aperçu de la méthodologie proposée, incluant l’acquisition d’images thermographiques, le prétraitement, l’extraction des champs de température, la formation CNN, l’apprentissage contraint par la physique et la quantification de l’incertitude. Abréviations ; CNN = réseau de neurones convolutionnel. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 4 : Pipeline global de traitement basé sur CNN. Architecture du cadre CNN proposé montrant la transformation des images thermographiques infrarouges en prédictions spatiales de champ de température. Abréviations ; CNN = réseau de neurones convolutionnel. Veuillez cliquer ici pour voir une version agrandie de cette figurine.