$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Esta seção descreve o pipeline completo de aprendizado profundo restrito pela física para avaliação de condições térmicas usando imagens termográficas infravermelhas. O framework integra (i) extração de campo de temperatura, (ii) regressão de CNN de linha de base, (iii) regularização baseada em física baseada em difusão, (iv) aprendizado adaptativo de parâmetros e (v) quantificação da incerteza. Todas as ferramentas usadas neste estudo são descritas na Tabela de Materiais.
Aquisição de imagem térmica e construção de campos de temperatura
Imagens infravermelhas (IR) de equipamentos de alta tensão foram adquiridas usando uma câmera termogràfica calibrada. Cada imagem era convertida em uma matriz de temperatura espacial T(x, y) ∈R H×W mapeando intensidades de pixels para valores físicos de temperatura (°C) usando parâmetros de calibração radiométrica da câmera. A Figura 1A–C mostra a visualização da temperatura. A imagem IR original é convertida em um mapa de temperatura fisicamente significativo. O histograma mostra uma distribuição de temperatura enviesada com concentração entre 50 e 90 °C. Zonas localizadas de alta temperatura (hotspots) são claramente visíveis. A Figura 2 mostra uma representação 3D térmica da elevação da superfície do campo de temperatura. O gráfico 3D de elevação da superfície revela inclinações acentuadas ao redor das regiões de pontos quentes. Os picos de temperatura se aproximam de ~120 °C. O campo térmico é espacialmente suave, exceto perto dos condutores energizados e das juntas de conexão. Essas características justificam o uso de um termo de regularização baseado em difusão para impor a suavidade física enquanto preserva fenômenos localizados de alto gradiente.

Figura 1: Visualização de Temperatura. (A) Imagem termográfica original infravermelha (IR) do equipamento de alta voltagem inspecionado mostrando o padrão térmico observado. (B) Mapa de temperatura extraído obtido da imagem termográfica, ilustrando a distribuição espacial das temperaturas superficiais e os gradientes térmicos associados. (C) Histograma da distribuição de temperatura mostrando a distribuição de frequência dos valores de temperatura dentro do campo de temperatura extraído. Esta figura demonstra a conversão de uma imagem termográfica em uma representação quantitativa da temperatura e resume a distribuição subjacente de temperatura usada para o desenvolvimento e análise subsequentes do modelo. Abreviações; IR = infravermelho. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 2: Representação 3D térmica da elevação da superfície do campo de temperatura. Visualização tridimensional da distribuição de temperatura, onde a elevação da superfície corresponde à magnitude da temperatura e destaca pontos quentes térmicos localizados. Por favor, clique aqui para ver uma versão ampliada desta figura.
Formulação do problema
Seja

denotam a imagem termográfica infravermelha RGB, onde H e W representam as dimensões espaciais e C = 3 corresponde aos canais de cor associados à paleta térmica arco-íris. Embora a imagem infravermelha represente visualmente informações térmicas, ela não é uma matriz numérica direta de temperatura porque as cores exibidas dependem do mapeamento da paleta, escala, interpolação e condições de imagem.
Seja

denotar o campo de temperatura correspondente de precisão no solo obtido por calibração radiométrica e extração de temperatura do sistema de imagem infravermelha.
O objetivo de aprendizagem é aproximar o mapeamento

onde θ denota os parâmetros treináveis da CNN e T̂ representa o campo de temperatura previsto. O objetivo desse mapeamento supervisionado não é simplesmente reproduzir a conversão radiométrica direta, mas aprender um modelo de regressão de temperatura espacialmente consistente a partir de imagens termográficas renderizadas. Ao contrário da extração direta de cor a temperatura pixel a pixel, a CNN pode explorar a estrutura térmica contextual em regiões vizinhas, reduzir a sensibilidade ao ruído, preservar características de hotspots e gerar campos de temperatura fisicamente plausíveis mais suaves sob as restrições de difusão impostas. A função objetivo de linha de base minimiza o erro de regressão pixel a pixel entre os campos de temperatura previsto e calibrado:

onde Tij e T̂ij denotam a verdade fundamental e as temperaturas previstas na localização dos pixels (i, j), respectivamente. Assim, o arcabouço deve ser interpretado como um modelo de regressão termográfica regularizado pela física, em vez de um substituto da calibração radiométrica direta.
Condições de aquisição de imagem infravermelha
Imagens termográficos infravermelhas foram obtidas usando uma câmera de imagem térmica FLIR sob condições de operação externas de subestações. O conjunto de dados consistia em imagens termográficas RGB com resolução espacial de 640 × 480 pixels, registradas usando a paleta de cores arco-íris FLIR. Cada imagem incluía uma barra de escala térmica embutida indicando a faixa de temperatura associada ao quadro de aquisição. A inspeção termográfica abrangeu múltiplas categorias de equipamentos de alta tensão, incluindo transformadores, buchas, isoladores, conectores e estruturas de terminais. Imagens foram capturadas durante condições normais de operação para incluir tanto o comportamento térmico padrão quanto padrões localizados de hotspots relevantes para aplicações de monitoramento de condições. A aquisição de imagem foi realizada a partir de pontos de vista práticos de inspeção, tipicamente usados na termografia de campo. A orientação da câmera e o ângulo de visão variavam dependendo do equipamento, acessibilidade e disposição da subestação. Imagens térmicas foram visualmente analisadas para garantir que o equipamento alvo permanecesse claramente observável e que as regiões de pontos quentes fossem distinguíveis do fundo. O conjunto de dados foi organizado em pastas rotuladas de acordo com a categoria do equipamento e condição térmica. Antes do desenvolvimento do modelo, as imagens eram segmentadas e redimensionadas em um formato espacial consistente para extração de campo de temperatura e treinamento CNN. As imagens termográficas contêm distribuições térmicas adequadas para visão computacional e estudos de aprendizado profundo envolvendo detecção de hotspots, identificação de falhas, localização de anomalias e previsão de campo de temperatura. No presente trabalho, os perfis térmicos extraídos das imagens infravermelhas foram usados para estimativa de temperatura limitada pela física e avaliação térmica consciente da incerteza.
Conjunto de dados termográficos e particionamento de dados
O conjunto de dados utilizado neste estudo consistiu em 174 imagens termográficas infravermelhas adquiridas de equipamentos de subestação de alta tensão usando um sistema de imagem térmicaFLIR 15. O conjunto de dados incluía vários tipos de componentes do sistema de energia, incluindo transformadores, buchas, isoladores, conectores e estruturas de terminais de alta tensão operando em condições práticas ao ar livre. Cada imagem termográfica era convertida em um campo de temperatura espacial usando as informações radiométricas fornecidas pela câmera infravermelha. O conjunto de dados contém variações térmicas que vão desde condições normais de operação até regiões localizadas de hotspots associadas a estresse térmico elevado. Os dados foram divididos em subconjuntos de treinamento, validação e teste com uma divisão de 70:15:15, resultando em 121, 26 e 27 imagens, respectivamente. A partição foi realizada no nível da imagem porque o conjunto de dados consistia em cenas termográficas capturadas de forma independente, coletadas sob diferentes condições operacionais e pontos de vista. Foi tomado cuidado para evitar quadros quase duplicados entre subconjuntos a fim de reduzir o risco de vazamento de dados e melhorar a justiça da avaliação. Todas as imagens foram redimensionadas para uma resolução espacial uniforme antes do treinamento. A normalização de temperatura foi aplicada para estabilizar a otimização, preservando os gradientes térmicos relativos e as características dos pontos quentes importantes para a avaliação da condição.
Arquitetura CNN básica
Uma arquitetura de regressão codificador-decodificador totalmente convolucional foi adotada para preservar a resolução espacial durante a previsão de campo de temperatura a partir de imagens termográficas infravermelhas. Todas as imagens infravermelhas de entrada foram redimensionadas para uma resolução uniforme de 256 × 256 × 3 antes do treinamento. A seção do codificador consistia em quatro blocos convolucionais com profundidade de características progressivamente crescente. Cada bloco continha uma camada de convolução 3 × 3, seguida de normalização em lote e ativação ReLU. A amostragem espacial foi realizada usando convoluções strided em vez de max pooling para melhor preservar estruturas térmicas localizadas e informações de hotspots. A progressão do canal do codificador foi:
32 → 64 → 128 → 256
A seção do decodificador reconstruiu o campo de temperatura espacial usando camadas de convolução transpostas para upsampling com a progressão reversa das características:
256 → 128 → 64 → 32
Cada bloco decodificador também utilizava convolução transposta, normalização em lote e ativação do ReLU. Uma camada final de convolução 1 × 1 gerou o mapa de temperatura previsto de canal único:
T̂ ∈ R256×256
A rede, portanto, gera um campo espacial de temperatura em resolução total correspondente à imagem de entrada no infravermelho. Para suportar a quantificação da incerteza, camadas dropout com probabilidade p = 0,2 foram inseridas próximas aos estágios de gargalo e decodificador. Durante a inferência, o dropout permaneceu ativo para permitir a amostragem de Monte Carlo para estimativa preditiva da incerteza. Todos os pesos convolucionais foram inicializados usando inicialização He-normal, enquanto termos de viés foram inicializados a zero. O modelo foi treinado usando o otimizador Adam com uma taxa inicial de aprendizado de 1 × 10−3, tamanho de lote de 8 e decaimento de peso de 1 × 10−5. O treinamento foi realizado por 100 épocas usando decaimento da taxa de aprendizado do cosseno e paradas precoces com base na validação RMSE.
A função objetivo de base é definida como:

onde dados L denotam o erro quadrático médio pixel a pixel entre os campos de temperatura previsto e o campo de verdade no solo.
O desempenho foi avaliado usando: R2, RMSE e MAE.
Embora a CNN de base tenha capturado com sucesso a estrutura térmica geral e a distribuição dos pontos quentes, ocasionalmente produziu artefatos não físicos em regiões de baixo gradiente e um leve suavizamento excessivo próximo a limites térmicos nítidos. Essas limitações motivaram a introdução da regularização da física baseada em difusão na formulação subsequente do modelo.
Aprendizado limitado por física via difusão
Equação governante do calor
Para impor consistência térmica no campo de temperatura previsto, um prior baseado em física baseada em difusão foi incorporado na equação do calor. Sob condições quase-estacionárias, a temperatura superficial satisfaz aproximadamente

onde T denota temperatura, α é a difusividade térmica e S representa a geração local de calor. Como as amostras termográficas disponíveis são instantâneos de inspeção únicos em vez de sequências resolvidas no tempo, a derivada temporal não pode ser avaliada diretamente. Assim, o modelo assume condução local quase-estacionária e impõe

longe de regiões de aquecimento ativo.
O Laplaciano 2D foi discretizado usando o estêncil padrão de diferenças finitas:

O resíduo físico era então calculado sobre o mapa de temperatura previsto como

Para evitar ambiguidade nas fronteiras da imagem, o resíduo Laplaciano foi avaliado apenas em pixels internos válidos, onde todos os valores vizinhos exigidos pelo estêncil estão disponíveis. Nenhuma máscara explícita foi aplicada às regiões de pontos quentes, já que o objetivo do termo físico era regularizar todo o campo espacial em vez de suprimir zonas termicamente ativas. A estrutura do hotspot foi, portanto, aprendida por meio do termo de dados, enquanto a penalidade laplaciana incentivou a suavidade consistente em difusão ao longo do mapa previsto.
A função objetivo final combinava a perda de dados e a perda física como

onde λphys controla a força da regularização física. No estudo de ablação, foi examinado λphys ∈ {0, 10−4, 10−3, 10−2, 10−1} para avaliar o efeito do aumento da força física de restrição. Valores moderados melhoraram a generalização, enquanto ponderação excessivamente forte tendia a gradientes de hotspot muito suaves e acentuados. Esse comportamento é consistente com as estruturas térmicas observadas nas imagens infravermelhas, onde regiões fisicamente lisas coexistem com concentrações localizadas de alta temperatura.
Deve-se enfatizar que a formulação de difusão em regime estacionário usada neste trabalho não pretende representar o mecanismo completo de transferência de calor que rege a formação de pontos quentes em equipamentos de alta tensão. Na prática, o superaquecimento localizado ocorre de processos internos de geração de calor, como aquecimento de Joule, resistência de contato e degradação do isolamento. Como o conjunto de dados disponível consiste em imagens termográficas instantâneas sem informações temporais ou medições de fonte, a restrição de difusão é empregada como um termo de regularização espacial inspirado na física. Seu papel é incentivar campos de temperatura localmente suaves e fisicamente plausíveis, permitindo que as características dos pontos quentes sejam aprendidas a partir dos dados. Portanto, a difusão anterior deve ser interpretada como uma restrição de suavidade espacial, e não como um modelo abrangente da física térmica subjacente.
Difusividade térmica aprendível e aumento de termo fonte
Para melhorar a interpretabilidade física, a difusividade térmica efetiva foi tratada como um parâmetro escalar treinável. Seja T̂(x, y) o campo de temperatura previsto e α > 0 denotar a difusividade térmica aprendível. Para impor positividade, α foi parametrizado por meio de uma variável não restrita β como

Nessa formulação, o prior físico foi escrito em termos de um campo de forçamento residual efetivo R(x, y), onde R(x, y) representa a desviação local da difusão pura induzida por efeitos de aquecimento não modelados. A relação adaptativo-difusividade é, portanto, expressa como

Aqui, ∇2T̂(x, y) é o Laplaciano do campo de temperatura previsto, e R(x, y) é o termo fonte residual efetiva inferido implicitamente a partir dos dados termográficos. Como as amostras disponíveis são imagens instantâneas em vez de sequências térmicas resolvidas no tempo, R(x, y) não é supervisionado separadamente e serve apenas como uma quantidade física latente durante a otimização. A perda física correspondente foi definida como

Para o experimento aumentado por termo fonte, uma cabeça auxiliar de CNN foi adicionada para estimar um mapa espacial de fonte Ŝ(x, y). Neste caso, Ŝ(x, y) representa o campo local de geração de calor previsto, e a relação governante torna-se

A cabeça fonte era acoplada ao gargalo do codificador e consistia em um pequeno ramo convolucional que mapeava características latentes compartilhadas para um mapa de saída de canal único Ŝ(x, y). Não havia rótulos de origem com base no terreno; portanto, Ŝ(x, y) era aprendido apenas através do objetivo conjunto

onde L é a perda por regressão de temperatura, λphys controla a restrição de difusão, e λs regulariza a previsão da fonte para desencorajar mapas triviais de alta magnitude. No contexto baseado em snapshots considerado aqui, a formulação do termo fonte permaneceu subdeterminada e não melhorou a precisão da previsão.
Abandono de Monte Carlo para quantificação da incerteza
Para estimar a incerteza epistêmica no campo de temperatura previsto, o dropout de Monte Carlo foi aplicado durante a inferência. Camadas dropout com taxa de 0,2 foram colocadas após o bloco de gargalo e dentro dos estágios decodificadores da CNN. Durante os testes, essas camadas de dropout foram mantidas ativas, e 30 passagens estocásticas para frente foram realizadas para cada imagem infravermelha. Para a passagem estocástica k-ésima, o campo de temperatura previsto é denotado por

onde X é a imagem infravermelha de entrada e θk representa os parâmetros da rede amostrada induzidos pelo dropout. A média preditiva foi calculada como

e a variância preditiva pixel a pixel foi obtida como

com K = 30. O mapa de variância resultante foi usado para identificar regiões espaciais com menor confiança preditiva. Como rótulos de incerteza de verdade no terreno não estão disponíveis para regressão termográfica, a incerteza foi avaliada usando mapas de variância pixel a pixel e estatísticas resumidas calculadas sobre regiões de hotspot e fundo, em vez de métricas baseadas em calibração. A análise quantitativa mostrou que o desvio padrão preditivo médio foi de 13,046 °C, enquanto a região do ponto quente apresentou um desvio padrão maior de 15,934 °C em comparação com 12,743 °C na região de fundo. Isso corresponde a um aumento de aproximadamente 25% na incerteza nas áreas de hotspot em relação ao fundo. Esses resultados indicam que o modelo expressa maior incerteza em regiões termicamente complexas, particularmente em torno dos limites de pontos quentes e gradientes acentuados de temperatura, o que é consistente com a estrutura observada nas visualizações da superfície térmica. O fluxo de trabalho é resumido nas Figuras 3 e 4 , que mostram o pipeline geral de processamento baseado em CNN.

Figura 3: Fluxo de trabalho. Visão geral da metodologia proposta, incluindo aquisição de imagens termográficas, pré-processamento, extração de campo de temperatura, treinamento CNN, aprendizado limitado por física e quantificação de incertezas. Abreviações; CNN = rede neural convolucional. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 4: Pipeline geral de processamento baseado em CNN. Arquitetura do framework proposto para CNN que mostra a transformação de imagens termográficas infravermelhas em previsões espaciais de campo de temperatura. Abreviações; CNN = rede neural convolucional. Por favor, clique aqui para ver uma versão ampliada desta figura.