$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Segnali di ciclatura multicanale e distribuzione dei dati
I segnali ciclici multicanale rappresentativi e le distribuzioni di prestazione sono stati analizzati per caratterizzare il comportamento di risposta ciclica delle batterie agli ioni di litio. Lo studio ha utilizzato cinque dataset di batterie agli ioni di litio rappresentanti diversi sistemi chimici, inclusi l'ossido di litio e cobalto (LCO), l'ossido di nichel, cobalto e alluminio (NCA), l'ossido di nichel, manganese e cobalto (NMC) e il LFP (Tabella 1). Le variabili fondamentali del ciclo, tra cui tensione, corrente, capacità, temperatura, resistenza interna ed efficienza coulombica o energetica, sono state integrate con caratteristiche derivate, come le caratteristiche di capacità incrementale e di tensione differenziale, al fine di stabilire un quadro di caratterizzazione multicanale. I profili capacità-tensione e le distribuzioni di SOH per il dataset CALCE (Figura 3A), i profili medi di tensione e le distribuzioni di SOH per il dataset SANYO (Figura 3B), le traiettorie di invecchiamento per il dataset PANASONIC (Figura 3C), le mappe termiche di evoluzione della tensione per il dataset KOKAM (Figura 3D) e i profili di distribuzione della tensione per il dataset GOTION (Figura 3E) hanno evidenziato comportamenti di degrado distinti tra le diverse chimiche di batteria valutate. Ad esempio, il dataset GOTION IFP20100140A ha mostrato caratteristiche di tensione a gradino tipicamente associate ai sistemi LFP, mentre il dataset PANASONIC NCR18650BD ha mostrato un degrado progressivo della capacità accompagnato da spostamenti del profilo di tensione durante il ciclo a lungo termine.
| ID Dataset | Tipo/Modello Batteria | Chimica | Formato | Numero di Celle | Cicli Totali Analizzati | Finestre dei Primi Cicli Utilizzate | Punti di Campionamento per Ciclo | Canali Principali del Ciclo | Canali Aggiuntivi di Stato | Caratteristiche Derivate | Obiettivi di Predizione |
| Dataset #1 | CALCE LCO battery | LCO | Pouch | 38 | 4.216 | Primi 50, primi 100, primi il 20% della vita | 128 | Voltage, current, time, capacity | Temperature, internal resistance | dQ/dV, dV/dQ, charge duration, discharge energy | SOH, RUL, dischargeable energy |
| Dataset #2 | SANYO UR18650E | NCA | Cilindrico 18650 | 24 | 2.487 | Primi 50, primi 100, primi il 20% della vita | 128 | Voltage, current, time, capacity | Temperature, coulombic efficiency | dQ/dV, voltage plateau slope, energy efficiency | SOH, RUL, dischargeable energy |
| Dataset #3 | PANASONIC NCR18650BD | NCA | Cilindrico 18650 | 42 | 5.134 | Primi 50, primi 100, primi il 20% della vita | 128 | Voltage, current, time, capacity | Temperature, internal resistance, coulombic efficiency | dQ/dV, dV/dQ, charge duration, energy throughput | SOH, RUL, dischargeable energy |
| Dataset #4 | KOKAM SLPB533459H4 | NMC | Pouch | 19 | 1.968 | Primi 50, primi 100, primi il 20% della vita | 128 | Voltage, current, time, capacity | Temperature, internal resistance | dQ/dV, discharge duration, voltage slope | SOH, RUL, dischargeable energy |
| Dataset #5 | GOTION IFP20100140A | LFP | Prismatic | 57 | 6.482 | Primi 50, primi 100, primi il 20% della vita | 128 | Voltage, current, time, capacity | Temperature, internal resistance, energy efficiency | dQ/dV, plateau duration, charge/discharge energy | SOH, RUL, dischargeable energy |
Tabella 1: Caratteristiche dei dataset delle batterie agli ioni di litio, delle variabili di ciclaggio multicanale, delle caratteristiche ingegnerizzate e degli obiettivi di previsione utilizzati nel flusso di lavoro DL riproducibile. La tabella riassume i dataset delle batterie agli ioni di litio inclusi in questo studio, comprendenti la chimica della batteria, il formato della cella, il numero di celle, i cicli totali analizzati, le finestre di osservazione, la densità di campionamento, le variabili di ciclaggio multicanale, le caratteristiche elettrochimiche ingegnerizzate e gli obiettivi di previsione utilizzati per lo sviluppo e la valutazione del modello. Abbreviazioni: SOH = stato di salute; RUL = vita utile residua; LCO = ossido di cobalto e litio; NCA = ossido di nichel, cobalto e alluminio; NMC = ossido di nichel, manganese e cobalto; LFP = fosfato di ferro e litio.

Figura 3. Profili rappresentativi di ciclatura multicanale e distribuzioni dello stato di salute delle batterie in diversi dataset di batterie agli ioni di litio. (A) Profili capacità-tensione e distribuzioni dello stato di salute (SOH) per il dataset CALCE relativo all'ossido di cobalto al litio. (B) Profili medi di tensione e distribuzioni dello stato di salute (SOH) per il dataset SANYO relativo all'ossido di nichel, cobalto e alluminio. (C) Traiettorie di invecchiamento e distribuzioni cumulative dello stato di salute (SOH) per il dataset PANASONIC relativo all'ossido di nichel, cobalto e alluminio. (D) Mappe termiche di evoluzione della tensione e distribuzioni campione dello stato di salute (SOH) per il dataset KOKAM relativo all'ossido di nichel, manganese e cobalto. (E) Profili di distribuzione della tensione e relazioni tra stato di salute (SOH) e resistenza interna per il dataset GOTION relativo al fosfato di ferro e litio. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Prestazioni Predittive Complessive per l'Immaganazzinamento dell'Energia
Il flusso di lavoro proposto ha dimostrato errori di previsione ridotti e prestazioni predittive costanti in tutti i compiti di valutazione (Tabella 2). Per la stima dello stato di salute (SOH), il modello ha raggiunto un errore assoluto medio (MAE) di 0,021 nel dataset di validazione e di 0,024 nel dataset di test, con valori corrispondenti di R2 pari a 0,962 e 0,955. Questi risultati indicano prestazioni predittive stabili per la stima continua dello stato di salute della batteria nelle condizioni valutate. Il flusso di lavoro ha inoltre dimostrato una capacità predittiva affidabile per la stima della vita utile residua (RUL) e dell'energia scaricabile, ottenendo valori di R2 nel dataset di test rispettivamente pari a 0,927 e 0,941. Per valutare ulteriormente la stabilità predittiva, tutte le metriche di prestazione sono state calcolate su cinque inizializzazioni indipendenti del modello utilizzando partizioni identiche dei dataset. Gli scarti tipo risultanti e i coefficienti di variazione sono rimasti bassi nelle valutazioni ripetute, indicando una limitata sensibilità all'inizializzazione stocastica e un comportamento predittivo stabile in condizioni ripetute di addestramento. Inoltre, gli intervalli di confidenza al 95% hanno mostrato una variabilità ridotta tra le esecuzioni ripetute, confermando la robustezza e la riproducibilità del flusso di lavoro multicanale proposto.
| Attività di previsione | Divisione del dataset | MAE | RMSE | MAPE (%) | R² |
| Stima dell'SOH | Validazione | 0.021 | 0.029 | 2.47 | 0.962 |
| Stima dell'SOH | Test | 0.024 | 0.033 | 2.83 | 0.955 |
| Previsione dell'RUL | Validazione | 18.6 | 25.9 | 8.74 | 0.938 |
| Previsione dell'RUL | Test | 21.4 | 29.7 | 9.63 | 0.927 |
| Previsione dell'energia scaricabile | Validazione | 0.087 | 0.121 | 3.18 | 0.948 |
| Previsione dell'energia scaricabile | Test | 0.094 | 0.129 | 3.56 | 0.941 |
Tabella 2: Prestazioni predittive del flusso di lavoro DL riproducibile per la stima delle prestazioni delle batterie agli ioni di litio su dataset di validazione e di test. La tabella riassume le prestazioni predittive del flusso di lavoro proposto per la stima dello stato di salute (SOH), la previsione della vita residua utile (RUL) e la previsione dell'energia scaricabile, utilizzando dataset di validazione e di test. Le metriche di prestazione includono l'errore assoluto medio (MAE), l'errore quadratico medio (RMSE), l'errore percentuale assoluto medio (MAPE) e il coefficiente di determinazione (R2). Abbreviazioni: SOH = stato di salute; RUL = vita utile residua; MAE = errore assoluto medio; RMSE = errore quadratico medio; MAPE = errore percentuale assoluto medio.
Confronto delle prestazioni tra diverse architetture di apprendimento profondo
La valutazione su diverse architetture di rete e configurazioni di input ha dimostrato una notevole variabilità nelle prestazioni predittive (Tabella 3). Le distribuzioni dell'RMSE prima e dopo la selezione dell'architettura basata sulla riproducibilità sono mostrate in Figura 4A. L'architettura della rete neurale profonda ha prodotto valori di MAE e RMSE relativamente più elevati rispetto alle architetture basate su sequenze e a quelle ibride. I confronti tra diverse configurazioni di input multicanale sono riassunti in Figura 4B. I modelli che integravano l'estrazione di caratteristiche locali e la modellazione temporale incrociata tra cicli hanno mostrato un'accuratezza predittiva migliore nelle condizioni valutate. I modelli con minore variabilità predittiva si sono generalmente raggruppati in corrispondenza di punteggi medi di prestazione più elevati e con deviazioni standard ridotte (Figura 4C). Le relazioni tra prestazione media e deviazione standard nelle diverse configurazioni di input valutate sono illustrate in Figura 4D. Tra le configurazioni testate, l'architettura CNN–LSTM ha ottenuto il valore di MAE più basso (0.024) e il valore di R2 più alto (0.955) all'interno dei dataset e delle impostazioni sperimentali valutati. Inoltre, l'integrazione di caratteristiche progettate con input multicanale completi ha ridotto il MAE a 0.022 con un valore di R2 pari a 0.961, mentre gli input basati esclusivamente sulla capacità hanno mostrato prestazioni predittive ridotte.
| Categoria | Modello / Configurazione dell'input | MAE | RMSE | MAPE (%) | R² |
| Architettura del deep learning | DNN | 0.031 | 0.042 | 3.58 | 0.928 |
| CNN | 0.027 | 0.037 | 3.14 | 0.941 |
| LSTM | 0.026 | 0.036 | 3.05 | 0.944 |
| CNN–LSTM | 0.024 | 0.033 | 2.81 | 0.955 |
| Transformer | 0.025 | 0.034 | 2.89 | 0.951 |
| Configurazione dell'input | Solo capacità | 0.034 | 0.046 | 3.92 | 0.916 |
| Tensione + corrente | 0.029 | 0.039 | 3.29 | 0.936 |
| Input multicanale completo | 0.024 | 0.033 | 2.83 | 0.955 |
| Multicanale + caratteristiche progettate | 0.022 | 0.031 | 2.61 | 0.961 |
Tabella 3: Confronto delle prestazioni predittive tra le diverse architetture di DL e configurazioni di input. La tabella confronta le prestazioni predittive di diverse architetture di apprendimento profondo (DL) e configurazioni di input utilizzate nel flusso di lavoro multicanale riproducibile. Le metriche di prestazione includono l'errore assoluto medio (MAE), l'errore quadratico medio (RMSE), l'errore percentuale assoluto medio (MAPE) e il coefficiente di determinazione (R2). Abbreviazioni: DNN = deep neural network; CNN = convolutional neural network; LSTM = long short-term memory; MAE = mean absolute error; RMSE = root mean square error; MAPE = mean absolute percentage error.

Figura 4. Confronto tra architetture di modelli e configurazioni di input guidato dalla riproducibilità. (A) Distribuzioni dell'errore quadratico medio prima e dopo la selezione dell'architettura basata sulla riproducibilità. (B) Distribuzioni dell'errore quadratico medio tra diverse configurazioni di input prima e dopo la selezione del modello stabile. (C) Confronto tra prestazione media e deviazione standard delle architetture di modelli candidate, con colori che indicano l'errore quadratico medio e cerchi vuoti che indicano i modelli stabili selezionati. (D) Confronto tra prestazione media e deviazione standard delle configurazioni di input, con colori che indicano l'errore quadratico medio e cerchi vuoti che indicano le configurazioni ottimali selezionate. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Valutazione della Distribuzione degli Errori e della Coerenza
Le distribuzioni degli errori assoluti sono state analizzate per valutare la stabilità delle previsioni nei confronti di basi di riferimento tradizionali di ML e architetture di DL. Le basi di riferimento tradizionali di ML hanno mostrato distribuzioni di errore più ampie e code di errore elevato prolungate, indicando una stabilità ridotta per compiti di previsione di modelli di degrado complessi. Tra i modelli di DL valutati, l'architettura CNN–LSTM ha dimostrato distribuzioni di errore relativamente più ristrette e code di errore elevato più corte, suggerendo un'integrazione migliore dell'estrazione di caratteristiche locali e della modellazione di sequenze temporali (Figura 5). L'inclusione di variabili relative a temperatura, resistenza interna ed efficienza è stata associata a riduzioni progressive della dispersione dell'errore di previsione nei diversi dataset valutati.

Figura 5. Distribuzioni di errore assoluto per basi di confronto convenzionali di apprendimento automatico (ML) e modelli DL multicanale. Diagrammi a violino che confrontano la distribuzione degli errori assoluti di previsione tra basi di confronto convenzionali di ML, architetture DL e configurazioni di ingresso multicanale. Le linee nere indicano l'errore medio assoluto, mentre le linee colorate indicano la media ± deviazione standard. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Riproducibilità e robustezza in diversi scenari di validazione
Il flusso di lavoro ha mantenuto prestazioni predittive costanti attraverso ripetuti cicli di addestramento e diverse strategie di validazione (Tabella 4). Le ripetizioni effettuate utilizzando partizioni fisse dei dati e inizializzazione fissa del seme casuale (random-seed) hanno prodotto metriche di prestazione stabili, con valori di R2 costantemente superiori a 0,95 e basse deviazioni standard nell'errore assoluto medio (MAE). Un’ulteriore valutazione mediante validazione incrociata a cinque ripiegamenti (five-fold cross-validation) e dataset di validazione esterna ha mostrato una limitata variazione delle prestazioni in condizioni sperimentali ripetute. La fase di validazione esterna ha utilizzato esclusivamente il dataset GOTION, che è stato completamente escluso dalle procedure di addestramento del modello e di ottimizzazione degli iperparametri. Poiché questo dataset rappresenta una chimica LFP distinta e una configurazione a celle prismatiche, la valutazione esterna ha fornito una rigorosa valutazione della capacità di generalizzazione tra dataset in condizioni elettrochimiche e strutturali eterogenee. Inoltre, le analisi del coefficiente di variazione e le valutazioni su esecuzioni ripetute hanno dimostrato una bassa dispersione predittiva tra esperimenti indipendenti, confermando la robustezza e la riproducibilità del flusso di lavoro multicanale proposto.
| Impostazione della validazione | Metrica di prestazione | Media | DS | Minimo | Massimo | Coefficiente di variazione (%) |
| Esecuzioni ripetute con suddivisione fissa dei dati (n = 5) | MAE | 0.023 | 0.0021 | 0.021 | 0.026 | 9.13 |
| RMSE | 0.032 | 0.0028 | 0.029 | 0.036 | 8.75 |
| MAPE (%) | 2.76 | 0.24 | 2.48 | 3.11 | 8.7 |
| R² | 0.957 | 0.006 | 0.949 | 0.964 | 0.63 |
| Esecuzioni ripetute con diversi semi casuali (n = 5) | MAE | 0.024 | 0.0024 | 0.021 | 0.027 | 10 |
| RMSE | 0.033 | 0.0031 | 0.029 | 0.037 | 9.39 |
| MAPE (%) | 2.84 | 0.27 | 2.52 | 3.19 | 9.51 |
| R² | 0.954 | 0.007 | 0.946 | 0.962 | 0.73 |
| Validazione incrociata a cinque ripiegamenti | MAE | 0.025 | 0.003 | 0.022 | 0.029 | 12 |
| RMSE | 0.034 | 0.0034 | 0.03 | 0.039 | 10 |
| MAPE (%) | 2.95 | 0.31 | 2.57 | 3.41 | 10.51 |
| R² | 0.951 | 0.009 | 0.94 | 0.961 | 0.95 |
| Validazione esterna | MAE | 0.027 | 0.0032 | 0.023 | 0.031 | 11.85 |
| RMSE | 0.037 | 0.0038 | 0.032 | 0.042 | 10.27 |
| MAPE (%) | 3.18 | 0.35 | 2.76 | 3.67 | 11.01 |
| R² | 0.944 | 0.011 | 0.931 | 0.956 | 1.17 |
Tabella 4: Analisi della riproducibilità e robustezza del flusso di lavoro multicanale basato su DL proposto attraverso ripetuti cicli di addestramento e strategie di validazione. La tabella riassume le prestazioni in termini di riproducibilità e robustezza del flusso di lavoro proposto in ripetuti cicli con partizioni fisse dei dati, ripetuti cicli con diversi semi casuali, validazione incrociata a cinque ripiegamenti e validazione esterna. Le metriche di prestazione includono l'errore assoluto medio (MAE), l'errore quadratico medio (RMSE), l'errore percentuale assoluto medio (MAPE), il coefficiente di determinazione (R2), la deviazione standard (SD) e il coefficiente di variazione. Abbreviazioni: MAE = errore assoluto medio; RMSE = errore quadratico medio; MAPE = errore percentuale assoluto medio; SD = deviazione standard.
Sensibilità agli Iperparametri e Ablazione delle Caratteristiche
L'analisi della sensibilità agli iperparametri ha dimostrato che le prestazioni di previsione variavano in base alle dimensioni della finestra di osservazione, alle funzioni di attivazione, alle dimensioni della dimensione nascosta e alla profondità del modello sequenziale. L'effetto delle dimensioni della finestra di osservazione sulle distribuzioni dell'errore di previsione è illustrato nella Figura 6A. Nelle condizioni valutate, finestre di osservazione comprese tra 100 e 200 cicli erano associate a errori di previsione inferiori, mentre finestre di osservazione più lunghe mostravano miglioramenti decrescenti nell'accuratezza predittiva. La variabilità della previsione dipendente dalla funzione di attivazione è riassunta nella Figura 6B. Le funzioni di attivazione rectified linear unit (ReLU) e Gaussian error linear unit (GELU) hanno mostrato errori di previsione inferiori rispetto ad altre funzioni di attivazione nelle condizioni valutate. La sensibilità delle prestazioni di previsione alla dimensione della dimensione nascosta e alla profondità del modello sequenziale è illustrata nella Figura 6C. Dimensioni nascoste pari a 128 e da due a tre strati di modellazione sequenziale erano associate a prestazioni predittive migliorate, mentre architetture poco profonde e dimensioni nascoste ridotte erano associate a un aumento dell'errore di previsione. L'analisi di ablazione delle caratteristiche ha inoltre dimostrato che la rimozione di caratteristiche ingegnerizzate o di variabili relative all'efficienza aumentava l'errore assoluto medio (MAE) rispetto alla configurazione multicanale completa (Tabella 5), suggerendo che l'uso combinato di ingressi multicanale e di caratteristiche elettrochimiche ingegnerizzate contribuiva a migliorare le prestazioni predittive.

Figura 6. Sensibilità delle prestazioni predittive alla variazione degli iperparametri. (A) Effetto della dimensione della finestra di input sulla distribuzione dell'errore di previsione. (B) Effetto della scelta della funzione di attivazione sulla distribuzione dell'errore di previsione. (C) Mappa termica che mostra la sensibilità della previsione alla dimensione dello strato nascosto e al numero di strati del modello sequenziale. Le linee nere indicano l'errore assoluto medio, mentre le linee colorate indicano la media ± deviazione standard, ove applicabile. Cliccare qui per visualizzare una versione ingrandita di questa figura.
| Categoria | Impostazione | MAE | RMSE | MAPE (%) | R² |
| Dimensione della finestra di input | 20 cicli | 0.031 | 0.043 | 3,74 | 0,927 |
| 50 cicli | 0,026 | 0,036 | 3,08 | 0,944 |
| 100 cicli | 0,023 | 0,032 | 2,71 | 0,957 |
| 150 cicli | 0,022 | 0,031 | 2,63 | 0,96 |
| 200 cicli | 0,022 | 0,031 | 2,59 | 0,961 |
| 250 cicli | 0,023 | 0,032 | 2,67 | 0,958 |
| 300 cicli | 0,024 | 0,033 | 2,79 | 0,955 |
| Funzione di attivazione | ReLU | 0,023 | 0,032 | 2,74 | 0,957 |
| GELU | 0,022 | 0,031 | 2,61 | 0,961 |
| Tanh | 0,026 | 0,036 | 3,05 | 0,946 |
| Sigmoid | 0,034 | 0,047 | 3,96 | 0,918 |
| Dimensione nascosta / Canali per livello | 32 | 0,028 | 0,039 | 3,22 | 0,939 |
| 64 | 0,024 | 0,033 | 2,81 | 0,954 |
| 128 | 0,022 | 0,031 | 2,58 | 0,961 |
| 256 | 0,023 | 0,032 | 2,69 | 0,958 |
| Numero di livelli di modellazione della sequenza | 1 | 0,027 | 0,038 | 3,14 | 0,942 |
| 2 | 0,023 | 0,032 | 2,73 | 0,956 |
| 3 | 0,022 | 0,031 | 2,6 | 0,961 |
| 4 | 0,024 | 0,034 | 2,86 | 0,952 |
| Ablazione delle caratteristiche | Caratteristiche multicanale complete + caratteristiche ingegnerizzate | 0,022 | 0,031 | 2,61 | 0,961 |
| Senza caratteristiche ingegnerizzate | 0,024 | 0,033 | 2,84 | 0,955 |
| Senza temperatura | 0,025 | 0,034 | 2,93 | 0,951 |
| Senza resistenza interna | 0,026 | 0,035 | 3,01 | 0,949 |
| Senza variabili relative all'efficienza | 0,025 | 0,034 | 2,95 | 0,95 |
| Solo tensione e corrente | 0,029 | 0,039 | 3,29 | 0,936 |
| Solo capacità | 0,034 | 0,046 | 3,92 | 0,916 |
Tabella 5: Analisi di sensibilità delle prestazioni del modello in diverse configurazioni di iperparametri e di ablazione delle caratteristiche. La tabella riassume gli effetti della dimensione della finestra di input, della funzione di attivazione, della dimensione della dimensione nascosta, della profondità del modello sequenziale e delle strategie di ablazione delle caratteristiche sulle prestazioni predittive all'interno del flusso di lavoro DL multicanale proposto. Le metriche di prestazione includono l'errore assoluto medio (MAE), l'errore quadratico medio (RMSE), l'errore percentuale assoluto medio (MAPE) e il coefficiente di determinazione (R2). Abbreviazioni: MAE = errore assoluto medio; RMSE = errore quadratico medio; MAPE = errore percentuale assoluto medio; ReLU = unità lineare rettificata; GELU = unità lineare di errore gaussiano.