ES2985237T3 - Procedimiento y dispositivo para la medición en tiempo real e in situ de los datos termodinámicos de una batería (entalpía y entropía) - Google Patents
Procedimiento y dispositivo para la medición en tiempo real e in situ de los datos termodinámicos de una batería (entalpía y entropía) Download PDFInfo
- Publication number
- ES2985237T3 ES2985237T3 ES19782947T ES19782947T ES2985237T3 ES 2985237 T3 ES2985237 T3 ES 2985237T3 ES 19782947 T ES19782947 T ES 19782947T ES 19782947 T ES19782947 T ES 19782947T ES 2985237 T3 ES2985237 T3 ES 2985237T3
- Authority
- ES
- Spain
- Prior art keywords
- battery
- electrical
- parameters
- phase
- estimate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/42—Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
- H01M10/48—Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/42—Methods or arrangements for servicing or maintenance of secondary cells or secondary half-cells
- H01M10/48—Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte
- H01M10/482—Accumulators combined with arrangements for measuring, testing or indicating the condition of cells, e.g. the level or density of the electrolyte for several batteries or cells simultaneously or sequentially
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
- H01M10/63—Control systems
- H01M10/633—Control systems characterised by algorithms, flow charts, software details or the like
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or discharging batteries or for supplying loads from batteries
- H02J7/80—Circuit arrangements for charging or discharging batteries or for supplying loads from batteries including monitoring or indicating arrangements
- H02J7/82—Control of state of charge [SOC]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or discharging batteries or for supplying loads from batteries
- H02J7/80—Circuit arrangements for charging or discharging batteries or for supplying loads from batteries including monitoring or indicating arrangements
- H02J7/84—Control of state of health [SOH]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J7/00—Circuit arrangements for charging or discharging batteries or for supplying loads from batteries
- H02J7/90—Regulation of charging or discharging current or voltage
- H02J7/971—Regulation of charging or discharging current or voltage the charge cycle being controlled or terminated in response to non-electric parameters
- H02J7/975—Regulation of charging or discharging current or voltage the charge cycle being controlled or terminated in response to non-electric parameters in response to temperature
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Engineering & Computer Science (AREA)
- General Chemical & Material Sciences (AREA)
- Manufacturing & Machinery (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Automation & Control Theory (AREA)
- Secondary Cells (AREA)
- Tests Of Electric Status Of Batteries (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
La invención se refiere a un procedimiento y a un dispositivo para la medida, en tiempo real e in situ, de datos termodinámicos de una batería (entalpía y entropía). El objeto de la invención es proporcionar un procedimiento fiable para la medida, in situ, online y en tiempo real, de la variación de entropía (ΔS) de una batería. Para ello, el procedimiento se caracteriza porque consiste principalmente en: (fase I) realizar un modelo previo de la batería: (a) cargar la batería; (b) y/o descargar la batería; (c) medir variables reales; (d) modelar el comportamiento eléctrico de la batería durante la carga (a) y/o descarga (b) para estimar los parámetros eléctricos de la batería; (e) estimar los parámetros eléctricos de la batería; (f) modelar el comportamiento térmico de la batería durante la carga (a) y/o descarga (b) para estimar ΔS in situ, online y en tiempo real; (g) estimar ΔS, mediante la utilización de al menos uno de los parámetros eléctricos estimados en la etapa (d); (fase II) medir ΔS de la batería durante su uso en cualquier aplicación y con cualquier estado de carga mediante la realización de la etapa (d) y la etapa (f) de la fase I, la etapa (e) y la etapa (g); (fase III) almacenar opcionalmente los datos medidos/calculados en la fase II y/o en la fase I. La invención se refiere también a un procedimiento para determinar el estado de carga y el estado de salud de una batería a partir de estos datos termodinámicos. Otro objeto de la invención es un dispositivo para la puesta en práctica de este procedimiento. (Traducción automática con Google Translate, sin valor legal)
Description
DESCRIPCIÓN
Procedimiento y dispositivo para la medición en tiempo real e in situ de los datos termodinámicos de una batería (entalpia y entropía)
Campo técnico
El campo de la invención es el de las baterías eléctricas.
Más concretamente, la invención se refiere a un procedimiento para medir en tiempo real e in situ los datos termodinámicos de una batería (entalpía y entropía).
La invención también se refiere a un dispositivo para implementar este procedimiento.
Estado de la técnica
El término “batería” designa un conjunto de acumuladores recargables conectados entre sí para obtener la capacidad y la tensión deseadas. Un acumulador se refiere a un dispositivo electroquímico unitario (celda) que comprende dos electrodos separados por un electrolito. En la presente discusión, el término “batería” se referirá tanto a un conjunto de celdas electroquímicas como a una celda unitaria.
Las baterías recargables, y las de ion Li en particular, ya se utilizan en una multitud de aplicaciones, como el almacenamiento de energía eólica y solar, los teléfonos móviles, los aviones y los coches. En este último sector, donde la cuota de mercado de los vehículos híbridos o eléctricos va a pasar del 1 % al 30 % en los próximos años, las baterías van a experimentar un crecimiento sin precedentes.
En todas estas aplicaciones, el diagnóstico de las baterías representa un reto importante en términos de rendimiento y seguridad. Las noticias están llenas de ejemplos de problemas de sobrecalentamiento, incendio o explosión de baterías.
Por lo tanto, es primordial desarrollar herramientas que permitan diagnosticar mejor el estado de las baterías y, en particular, el estado de carga y la salud de las baterías.
Con las perspectivas de desarrollo expuestas, optimizar el funcionamiento y la seguridad de las baterías, por poco que sea, tendría un impacto considerable. Esta optimización pasa por la mejora de los medios de diagnóstico del estado interno de las baterías.
Sabemos que existe un vínculo entre el estado interno de una batería y las magnitudes termodinámicas asociadas a él, en particular la variación de entalpía y la variación de entropía. Se ha observado que los perfiles termodinámicos dependen de determinadas características de la batería. Estos perfiles dependen del envejecimiento térmico, del envejecimiento por carga-descarga, del envejecimiento por sobrecarga y también del estado de carga de la batería. También se ha observado una modificación de los perfiles cuando la batería se sobrecarga. El conocimiento de los datos termodinámicos es, por lo tanto, una rica fuente de información sobre el estado y la historia de una batería.
Además, el conocimiento de la termodinámica de las reacciones en los electrodos y de las transformaciones físicas que tienen lugar en las baterías es esencial para predecir su rendimiento y estabilidad.
Así, la densidad de energía refleja las cantidades totales de cargas intercambiadas reversiblemente y el potencial al que se producen estos intercambios. Además, el ciclo de vida se refiere a la estabilidad de los estados o fases que resultan de las transformaciones en los electrodos en los procesos de carga y descarga. Estos procesos están controlados por la termodinámica.
Por lo tanto, existe un vínculo entre el estado interno de la batería y las magnitudes termodinámicas asociadas a él (variación de entalpía y variación de entropía).
El estado interno de una batería está dado en particular por las siguientes características: envejecimiento térmico, envejecimiento por carga-descarga, envejecimiento por sobrecarga, estado de carga de la batería, sobrecarga. El conocimiento de los datos termodinámicos es, por lo tanto, una rica fuente de información sobre el estado y la historia de una batería.
Se han desarrollado y aplicado varias técnicas para evaluar la cinética termoquímica de las reacciones en los electrodos, en particular procedimientos electroanalíticos (voltamperometría cíclica, potenciometría, etc.) y técnicas espectroscópicas (difracción de rayos X, RMN, LEED, etc.).
Los documentos de patente EP1924849, US8446127, US9599584, US20160146895, WO2017204750, entre otros, describen procedimientos para evaluar el estado de carga y el estado de salud de las baterías mediante perfiles termodinámicos medidos.
En este estado de la técnica, el procedimiento para medir los perfiles termodinámicos es un procedimiento estándar que requiere condiciones de laboratorio y una cantidad significativa de tiempo del orden de varios días a varias semanas. Este procedimiento estándar, también conocido como procedimiento ETM “Electro-chemical Thermodynamic Measurement”, consiste esencialmente en utilizar un dispositivo que comprende:
1 Un dispositivo para controlar la temperatura de las baterías
2 Un potenciómetro-galvanómetro para medir la tensión y la corriente con gran precisión
3 Un ordenador para controlar el proceso de medición termodinámica y recopilar los datos
Para medir AS (la variación de entropía) con este dispositivo, el proceso ETM se basa en la siguiente fórmula:
en donde
• AS es la variación de entropía
• F es la constante de Faraday
• OCV es la tensión en circuito abierto
• Tbates la temperatura de la batería
• x es el estado de carga
En la práctica, para medir AS, se coloca la batería en un estado de carga determinado, se desconecta la corriente y se deja que la tensión en los terminales de la batería se relaje hacia OCV durante varias horas. A continuación, se varía la temperatura de la batería, lo que tiene el efecto de variar la tensión. Existe una relación lineal entre Tbaty OCV. Solo queda determinar el coeficiente que une a ambos y multiplicarlo por la constante de Faraday para obtener AS. Para obtener un perfil completo, la medición debe repetirse en diferentes estados de carga.
OCV se obtiene directamente midiendo la tensión después de la relajación y AH (la variación de entalpía) se obtiene mediante la ecuación (4)
Por lo tanto, parece que esta medición ETM estándar de los perfiles termodinámicos es una técnica de laboratorio fastidiosa y lenta, que requiere cierto dominio de las condiciones de funcionamiento de la batería. Este procedimiento ETM estándar no puede utilizarse in situ, en línea y en tiempo real.
En estas circunstancias, la presente invención pretende satisfacer al menos uno de los objetivos que se exponen a continuación.
CjsrUno de los objetivos esenciales de la presente invención es proporcionar un procedimiento para medir, in situ, en línea y en tiempo real, la variación de entropía de una batería.
^"Uno de los objetivos esenciales de la presente invención es proporcionar un procedimiento sencillo y fiable para medir, in situ, en línea y en tiempo real, la variación de entropía de una batería.
^"Uno de los objetivos esenciales de la presente invención es proporcionar un procedimiento eficaz para medir, in situ, en línea y en tiempo real, la variación de entropía de una batería.
^"Uno de los objetivos esenciales de la presente invención es proporcionar un procedimiento rápido para medir, in situ, en línea y en tiempo real, la variación de entropía de una batería.
^"Uno de los objetivos esenciales de la presente invención es proporcionar un procedimiento rentable para medir, in situ, en línea y en tiempo real, la variación de entropía de una batería.
^"Uno de los objetivos esenciales de la presente invención es proporcionar un procedimiento económico de medir, in situ, en línea y en tiempo real, la variación de entropía de una batería, es decir, en condiciones más próximas al uso normal de la batería.
^"Uno de los objetivos esenciales de la presente invención es proporcionar un procedimiento simple y eficaz de diagnóstico del estado de carga y/o del estado de salud de las baterías, sobre la base de la variación de entropía medida por el procedimiento mencionado en los objetivos anteriores.
^ U n o de los objetivos esenciales de la presente invención es proporcionar un dispositivo sencillo y de alto rendimiento para llevar a cabo el procedimiento mencionado en los objetivos anteriores.
Breve descripción
Estos objetivos, entre otros, se consiguen mediante la presente invención, que se refiere en primer lugar a un procedimiento de medición en línea e in situ de la variación de entropía AS de al menos una batería,caracterizado porqueconsiste esencialmente en:
(Fase I) Realizar una modelización previa de la batería cuyo estado de carga (SOC: “State Of Charge”) está comprendido entre el 0 y el 100 % mediante la ejecución de las siguientes etapas:
(a) cargar al menos parcialmente la batería con una señal Sc de corriente de carga;
(b) y/o descargar al menos parcialmente la batería con una señal Sd de corriente de descarga;
(c) medir variables reales útiles en las etapas siguientes;
(d) modelizar el comportamiento eléctrico de la batería durante la carga (a) con una señal Sc de corriente de carga y/o la descarga (b) con una señal Sd de corriente de descarga, a fin de estimar los parámetros eléctricos de la batería; (e) estimar periódicamente los parámetros eléctricos de la batería con una frecuencia Fe;
(f) modelizar el comportamiento térmico de la batería durante la carga (a) con una señal Sc de corriente de carga Sc y/o durante la descarga (b) con una señal Sd de corriente de descarga, a fin de estimar in situ, en línea y en tiempo real, al menos uno de los parámetros del modelo térmico, a saber AS;
(g) estimar periódicamente, a una frecuencia Fg, al menos uno de los parámetros del modelo térmico, en particular AS, utilizando al menos uno de los parámetros eléctricos estimados en la etapa (e);
(Fase II) Medir los datos termodinámicos, en particular AS, de la batería en uso en una aplicación y en cualquier estado de carga, implementando los modelos eléctricos [etapa (d)] y térmico [etapa (f)] de la fase I, la estimación de los parámetros eléctricos [etapa (e)] y la estimación de al menos uno de los parámetros del modelo térmico, en particular AS [etapa (g)];
(Fase III) Posiblemente almacenar los datos medidos/calculados durante la fase II y/o la fase I.
El procedimiento es esencialmente tal que:
^ La modelización de la etapa d) consiste en considerar que la batería es un circuito eléctrico o modelo eléctrico que comprende una resistencia R<0>, una tensión de circuito abierto OCV, y un circuito R<1>C<1>en serie;
el comportamiento eléctrico de la batería se describe, en este modelo, mediante las siguientes ecuaciones;
en donde U<1>es la tensión a través del circuito RC, / es la corriente a través de la batería y Vbat es la tensión a través de los bornes de la batería.
La ecuación (2) se discretiza como sigue:
Vbat.k ~ ^k^O.k3”^k-1 1^,k 3”al , k ^ ^ ^ k - l ^bat,/c-l) 3"® ^ k(2') y se reescribe como sigue:
con:
en donde
Ts es el período de muestreo de la señal Se;
u kes el vector de parámetros;
^ La modelización de la etapa (f) consiste en considerar la batería como un modelo térmico en el que la batería, por un lado, está sometida a una corriente Sc de carga susceptible de ser sometida a muestreo o a una corriente Sd de descarga susceptible de ser sometida a muestreo, y, por otro lado, es un intercambiador de calor con su entorno y en el que el comportamiento térmico de la batería se describe mediante la siguiente ecuación:
en donde:
• m es la masa de la batería
• Cp es la capacidad calorífica de la batería
• Tbates la temperatura de la batería
• t es la variable temporal
• / es la corriente que circula por la batería
• Vbat es la tensión en los bornes de la batería
• OCV es la tensión en circuito abierto de la batería
• AS es la variación de entropía de la batería
• F es la constante de Faraday
• h es el coeficiente de intercambio de calor con el exterior
• A es el área de la batería en contacto con el exterior
• Tamb es la temperatura del ambiente exterior
en donde la ecuación (3) se discretiza como sigue:
y se reescribe como sigue:
con:
en donde:
Ts es el período de muestreo de la medición.
Por lo tanto, este procedimiento puede utilizarse para estimar un gran número de parámetros de la batería en cuestión, en particular la fuerza electromotriz (“Open Circuit Voltage-OCV”), la variación de entropía AS y, por lo tanto, la variación de entalpia AH. Además de estos parámetros, el procedimiento también puede utilizarse para estimar la resistencia interna y la capacidad calorífica de la batería.
Este procedimiento es especialmente ventajoso porque permite medir en línea, en tiempo real e in situ los parámetros de una batería en uso, por ejemplo, en un vehículo eléctrico.
Una vez estimada, por ejemplo, la variación de entropía, es posible estimar el estado de salud de la batería, conocimiento fundamental para una buena gestión de la misma.
Así, en otro de sus aspectos, la presente invención se refiere a un procedimiento para determinar el estado de carga (SOC) y/o el estado de salud (SOH: “State Of Health”) de una batería a partir de AS y/o AH medidos por el procedimiento según la invención.
La presente invención también se refiere a un dispositivo para implementar el procedimiento de medición de variación de entropía según la invención.
Definiciones
A lo largo del presente documento, el singular se refiere tanto al singular como al plural. Las definiciones que se dan a continuación, a modo de ejemplo, pueden utilizarse en la interpretación de este documento:
• “medición de datos termodinámicos in situ”: La medición in situ se refiere a una medición durante el uso de la batería en una aplicación, por ejemplo: batería utilizada en un vehículo eléctrico. Esto contrasta con una medición que tendría lugar en un laboratorio.
• “medición de datos termodinámicos en línea”: La medición en línea es una medición que estima los parámetros o estados de un modelo (etapa (d) modelo eléctrico o etapa (f) modelo térmico), cuando se dispone de nuevos datos durante el uso de la batería en una aplicación.
• “batería”: un dispositivo electroquímico unitario (celda) que comprende dos electrodos separados por un electrolito, o un conjunto de acumuladores conectados entre sí para obtener la capacidad y la tensión deseadas. • aproximadamente” o “sustancialmente” significa dentro de más o menos el 10 %, o incluso dentro de más o menos el 5 %, de la unidad de medición utilizada.
• “ incluido entre Z1 y ZZ” significa que uno y/u otro de los terminales Z1, Z2 está incluido o no en el intervalo [Z1, Z2].
Descripción detallada
Esta descripción se hace con referencia a las figuras adjuntas en las que:
- La figura 1 es un diagrama esquemático del modelo eléctrico implementado en la etapa (d) del procedimiento según la invención.
- La figura 2 es un diagrama esquemático del uso de 2 algoritmos de mínimos cuadrados recursivos RLS en las etapas (e) y (g) de la estimación periódica, respectivamente.
- La figura 3 es un diagrama del banco de pruebas utilizado en el ejemplo de implementación del procedimiento según la invención.
- La figura 4 en la parte superior es una curva que muestra la respuesta de la temperatura de una batería 18650 a la corriente, en el ejemplo según la invención, con vistas a estimar mCp y hA.
- La figura 4 en la parte inferior es una curva que muestra la respuesta de la tensión de una batería 18650 a la corriente, en el ejemplo según la invención, con vistas a estimar mCp y hA.
- La figura 5 muestra una curva que estima el parámetro mCp en función del tiempo.
- La figura 6 es una curva de estimación del parámetro hA en función del tiempo.
- La figura 7 muestra la corriente de entrada Se o corriente de carga Sc constituida por una señal SBPA en A en función del tiempo (s).
- La figura 8 muestra la respuesta de la tensión (V) de una batería 18650 a la corriente, en el ejemplo, con vistas a estimar R<0>y OCV.
- La figura 9 muestra una estimación de R<0>(ohmios) en función del SOC (%) durante la descarga.
- La figura 10 muestra una estimación de OCV (V) en función del estado de carga SOC (%) durante la descarga. - La figura 11 es una curva de AS (J/K/mol) en función del estado de carga SOC (%) durante la descarga.
- La figura 12 es una curva de AH (kJ/mol) en función del estado de carga SOC (%) durante la descarga.
- La figura 13 muestra el banco de pruebas para una medición estándar de los datos termodinámicos del ejemplo comparativo.
- La figura 14 muestra curvas comparativas de los perfiles OCV(V) en función del estado de carga SOC (%) obtenidos por el procedimiento según invención (“procedimiento online”) y por ETM en el ejemplo comparativo (procedimiento estándar) en descarga.
- La figura 15 muestra curvas comparativas de los perfiles de AS (J/K/mol) en función del estado de carga SOC (%) obtenidos por el proceso según invención (procedimiento online) y por ETM en el ejemplo comparativo (procedimiento estándar) en descarga.
- La figura 16 muestra curvas comparativas de los perfiles de AH (kJ/mol) en función del estado de carga de SOC (%) obtenidos por el procedimiento según la invención (procedimiento en línea) y por ETM en el ejemplo comparativo (procedimiento estándar) en descarga.
- La figura 17 muestra las curvas OCV(V) en función del estado de carga SOC (%) obtenido por el procedimiento según la invención, medidas varias veces en las mismas condiciones.
- La figura 18 muestra curvas de AS (J/K/mol) en función del OCV(V) obtenido por el procedimiento según la invención, medido varias veces en las mismas condiciones.
- La figura 19 muestra curvas de AH (kJ/mol) en función del OCV (V) obtenido por el procedimiento según la invención, medido varias veces en las mismas condiciones.
Procedimiento para medir el cambio de entropía AS
Fase I de modelización
La fase I de modelización del procedimiento según la invención comprende ventajosamente las siguientes etapas: Etapa (ao) opcional pero no obstante preferida: estimación previa de los parámetros mCp y hA
Etapa (a) carga de la batería con una señal Sc de corriente de carga
Etapa b) y/o descarga de la batería con una señal Sd de corriente de carga
Etapa c) medición de las variables reales
Etapa d) modelización eléctrica
Etapa e) estimación de los parámetros eléctricos
Etapa f) modelización térmica
Etapa g) estimación de los parámetros termodinámicos
Para las modelizaciones (d) y (f), la etapa de carga (a) y/o la etapa de descarga (b) se realizan en paralelo, es decir, simultáneamente, con las etapas (c) a (g).
En un modo de realización que comprende una etapa (b) de descarga al menos parcial durante la cual se llevan a cabo las etapas (c) a (g), se prevé previamente una etapa (a) para cargar completamente la batería, de acuerdo con la invención.
De acuerdo con una variante de este modo de realización, la etapa (a) previa es una carga parcial de la batería. En otro modo de realización que comprende una etapa (a) de carga, durante la cual se ejecutan las etapas (c) a (g), se prevé previamente una etapa (b) de descarga total de la batería, de acuerdo con la invención.
De acuerdo con una variante de este modo de realización, la etapa (b) previa es una descarga parcial de la batería. Independientemente de su creciente referenciación alfabética, las etapas (c) a (g) no son necesariamente sucesivas. Modelo termoeléctrico [Etapas (d) y (f)]
El enfoque adoptado por la invención para estimar los datos termodinámicos de una batería consiste en elegir una representación termoeléctrica simple del sistema que comprende la batería.
Esta representación termoeléctrica que tiene, además, la ventaja de poder aplicarse a una gran variedad de baterías, consta de un modelo eléctrico [Etapa (d)] y de un modelo térmico [Etapa (f)].
• El modelo eléctrico se utiliza para estimar los parámetros eléctricos de la batería, como la resistencia interna y la fuerza electromotriz. Se basa ventajosamente en el llamado “modelo de Thévenin”, que considera la batería como una resistencia en serie con un circuito RC. Como se muestra en la figura 1, este modelo consiste en una resistencia Ro, una tensión de circuito abierto OCV y un circuito R<1>C<1>en serie.
En este modelo, el comportamiento eléctrico de la batería también se describe mediante las siguientes ecuaciones:
en donde Ui es la tensión terminal del circuito RC, / es la corriente a través de la batería y Vbat es la tensión terminal de la batería.
A continuación, las ecuaciones se discretizan del mismo modo que en “V.-H. Duong, H. A. Bastawrous, K. Lim, K. W. See, P. Zhang, and S. X. Dou, “Online State of charge and model parameters estimation of the LiFePO<4>battery in electric vehicles using multiple adaptive forgetting factors recursive least-squares”, JOURNAL Of POWER SOURCES, vol. 296, pp. 215-224, Nov 20, 2015”. Se obtiene:
que se puede reescribir:
con:
en donde:
Ts es el período de muestreo.f í7kes el vector de parámetros, identificando este vector es posible estimar un valorpara Rq y OCV.
De acuerdo con una disposición preferida de realización de la invención, OCV se utiliza en el modelo térmico.
• El modelo térmico permite estimar los parámetros termodinámicos de la batería, como la variación de entropía, de forma precisa, continua y sin control específico de las condiciones de funcionamiento.
Este modelo térmico se basa en particular en el hecho de que el calor liberado o absorbido por la batería procede principalmente de 2 fenómenos:
- El efecto Joule y la variación de entropía. De hecho, como la batería tiene una resistencia interna, produce calor cuando fluye la corriente, durante la carga o la descarga.
- Además, el calor puede ser absorbido o emitido en función del signo de la variación de entropía y de la corriente que circula por la batería. También se intercambia calor entre la batería y su entorno.
La modelización de la etapa (f) consiste en considerar la batería como un modelo térmico en el que la batería, por un lado, está sometida a una corriente Sc de carga susceptible de ser sometida a muestreo o a una corriente Sd de descarga susceptible de ser sometida a muestreo, y, por otro lado, es un intercambiador de calor con su entorno y en el que el comportamiento térmico de la batería se describe mediante la siguiente ecuación:
en donde:
• m es la masa de la batería
• Cp es la capacidad calorífica de la batería
• Tbates la temperatura de la batería
• t es la variable temporal
• / es la corriente que circula por la batería
• Vbat es la tensión en los bornes de la batería
• OCV es la tensión en circuito abierto de la batería
• AS es la variación de entropía de la batería
• F es la constante de Faraday
• h es el coeficiente de intercambio de calor con el exterior
• A es el área de la batería en contacto con el exterior
• Tamb es la temperatura del ambiente exterior
La ecuación (3) se discretiza entonces:
que puede reescribirse:
con:
en donde:
Ts es el período de muestreo de la medición.
En una variante, R puede utilizarse en lugar o además de OCV en el modelo térmico.
Los dos modelos propuestos pueden escribirse como una combinación lineal de magnitudes medidas. Por lo tanto, pueden aplicarse algoritmos como el de mínimos cuadrados para estimar los coeficientes de esta combinación lineal y, por lo tanto, los parámetros.
Este procedimiento permite, por lo tanto, estimar un gran número de parámetros para la batería considerada, y puede utilizarse para estimar la fuerza electromotriz, la variación de entropía y, por lo tanto, la variación de entalpía. Además de estos parámetros, esta metodología también puede utilizarse para estimar, entre otras cosas, la resistencia interna, la capacidad calorífica de la batería, el estado de carga y el estado de salud de la batería.
Etapa (ao) opcional pero preferida: estimación previa de los parámetros mCp y hA
En un modo de realización preferido de la invención, los productos mCp y hA se consideran constantes con respecto al estado de carga y al estado de salud de la batería.
Por lo tanto, es ventajoso estimar estos dos parámetros de una vez por todas y utilizar estos valores en el modelo eléctrico según la invención, para una batería dada.
Esta estimación previa consiste esencialmente en:
- relajar la batería durante un período de tiempo (que varía según las baterías) comprendido entre 10 min y 60 min para que el parámetro OCV de la ecuación (0) del modelo eléctrico tenga un valor preciso dado;
- aplicar una señal eléctrica de entrada Se periódica cuyo período se elija de manera que la media del calor generado por AS durante un período sea del orden de 0, siendo Se preferiblemente una señal de onda cuadrada con un período comprendido entre 10 y 30 segundos, idealmente del orden de 20 segundos, teniendo esta señal de onda cuadrada más preferiblemente aún una media de cero;
- estimar los productos mCpy hA, preferiblemente mediante un algoritmo de mínimos cuadrados recursivos, utilizando la ecuación (3) que se convierte en la ecuación (3') siguiente:
y midiendo las variables reales correspondientes a los parámetros Vbat, I, Tbat, Tamb,
- integrar esta estimación de los productos mCp y hA en el modelo térmico para las etapas (f) y (g).
El período de relajación de la tensión de la batería puede ser de unos diez minutos, con lo que se obtiene un valor de OCV preciso.
A continuación, se utiliza preferiblemente un algoritmo de mínimos cuadrados recursivos (algoritmo RLS) para estimar estos dos parámetros.
Etapa (a) carga de la batería con una señal Sc de corriente de carga
Para llevar a cabo la etapa (a) de carga, se aplica repetidamente una corriente de carga Sc, cuyo intervalo de frecuencia está comprendido entre 0 y 1 Hz, preferiblemente una señal correspondiente a una Secuencia Binaria PseudoAleatoria -SBPA-, elegida de forma que la “tasa C” (C: capacidad de la batería) de la batería esté comprendida entre 0,01. C y 3.C, preferiblemente entre 0.1.C y 2.5.C, y aún más preferiblemente entre 0.2.C y 2.C.
Esto puede corresponder, por ejemplo, para una batería de litio dada, a una corriente de carga Sc cuya intensidad está comprendida, por ejemplo, entre 0 y 1A.
De acuerdo con un modo ventajoso de la invención, la corriente de carga Sc se muestrea cada segundo.
Etapa (b) descarga de la batería con una señal Sd de corriente de descarga
Para llevar a cabo la etapa (b) de descarga, se aplica repetidamente una corriente de descarga Sd, cuyo intervalo de frecuencia está comprendido entre 0 y 1 Hz, preferiblemente una señal correspondiente a una Secuencia Binaria Pseudoaleatoria -SBPA, elegida de modo que la “tasa C” (C: capacidad de la batería) de la batería esté comprendida entre 0.01.C y 3.C, preferiblemente entre 0.1.C y 2.5.C, y aún más preferiblemente entre 0.2.C y 2.C.
Esto puede corresponder, por ejemplo, para una batería de litio dada, a una corriente de descarga Sd cuya intensidad está comprendida, por ejemplo, entre 0 y -1 A.
De acuerdo con un modo ventajoso de la invención, la corriente de descarga Sd se muestrea cada segundo.
Etapa c) medición de las variables reales
Las variables reales medidas en la etapa (c) corresponden ventajosamente a los parámetros Vsat, I, Tbat, Tamb de los modelos eléctrico y térmico.
Etapa (d) modelización eléctrica / Etapa (f) modelización térmica
Véase más arriba.
Etapa (e) estimación de los parámetros eléctricos / Etapa (g) estimación de los parámetros termodinámicos
De acuerdo con una característica destacable de la invención, la etapa (e) se lleva a cabo mediante un algoritmo de mínimos cuadrados recursivos y la estimación según la etapa (g) se lleva a cabo mediante un algoritmo de mínimos cuadrados recursivos.
Para la estimación de los parámetros eléctricos en la etapa (e), la batería se carga y descarga de acuerdo con las etapas (a) y (b).
Los parámetros del modelo eléctrico que se tienen en cuenta preferiblemente son R<0>y OCV. R<0>es un parámetro utilizado para la estimación del envejecimiento de la batería. OCV es útil para el modelo térmico, para estimar AS.
El algoritmo RLS se aplica a los datos para estimar los parámetros del modelo. Los parámetros se actualizan cada 0,2 segundos.
Para estimar AS, el algoritmo RLS se aplica preferiblemente a los datos vinculados al modelo térmico.
Como se muestra en la Figura 2, dos algoritmos RLS trabajan simultáneamente. El primero, como se ha mencionado con anterioridad, estimará OCV a partir de la corriente y la tensión. Esta estimación será utilizada por un segundo algoritmo RLS basado en el modelo térmico y utilizando la temperatura para estimar AS. Hay que tener en cuenta que, como la dinámica eléctrica es mucho más rápida que la térmica, la frecuencia utilizada para actualizar los parámetros eléctricos es mayor que para los parámetros térmicos.
También se puede obtener una estimación de la entalpía AH a partir de OCV y AS de la siguiente manera:
AH = -F .O C V -T batAS(6)
El parámetro eléctrico OCV y los parámetros termodinámicos (AS y AH) se estiman así en tiempo real mediante un algoritmo recursivo.
En otro de sus aspectos, la invención se refiere a un procedimiento para determinar el estado de carga y/o el estado de salud de una batería a partir de AS y/o AH.
• Para estimar el SOC: Se ha demostrado que, para un valor dado de SOC, corresponde un único par de valores (AS y AH). Si estos se estiman, entonces el SOC también puede estimarse, de ahí la importancia de la estimación en tiempo real. Una implementación sencilla consiste en estimar el SOC considerándolo como una combinación lineal de AS y AH. El SOC puede escribirse SOC(%)= a AS p AH<y>con a, p y<y>parámetros por determinar para una batería dada. Para determinarlos, se miden las magnitudes termodinámicas AS y AH para una batería dada del 0 al 100 %. A continuación, se determina la relación entre AS, AH y SOC de forma que SOC(%)= a AS p AH<y>, a, p y<y>se determinen preferiblemente por el procedimiento de los mínimos cuadrados. Para un ejemplo de dicho procedimiento, puede hacerse referencia a la solicitud de patente US2016146895A1.
• Para la estimación del SOH, es posible estimar un indicador de SOH basado en la energía de una batería a partir de sus perfiles termodinámicos utilizando herramientas de “machine learning”.
Por ejemplo, las celdas electroquímicas de iones litio de tipo 18650 (65 mm de altura y 18 mm de diámetro) se someten a envejecimiento por ciclos (cargas y descargas sucesivas). Las células se envejecieron a una velocidad de 1,5 C y a 55 °C. En diferentes etapas del envejecimiento, se midieron un indicador de SOH basado en la energía de la batería y varios parámetros termodinámicos. Estos últimos se midieron mediante el procedimiento estándar ETM. A continuación, fue posible establecer un vínculo entre una de las variables termodinámicas, AS en este caso, y el indicador SOH basado en la energía utilizando herramientas de “machine learning”. Un ejemplo de algoritmo de “machine learning” utilizado para vincularlas es la regresión lineal múltiple. En este caso, el indicador SOH se estima como una combinación lineal de AS en determinados valores de SOC u OCV.
De acuerdo con un modo de realización, la ejecución de la etapa (e) permite estimar, a una frecuencia Fe comprendida entre 0 y 10 Hz, preferiblemente del orden de 1 Hz, R<0>y OCV en función del estado de carga de la batería.
De acuerdo con un modo de realización y en el que la realización de la etapa (g) permite estimar, a una frecuencia Fg comprendida entre 0 y 1 Hz, preferiblemente del orden de 0,2 Hz, AS en función del estado de carga de la batería.
Fase II de medición
Esta fase II de medición in situ, durante la utilización de la batería en una aplicación determinada, por ejemplo, la alimentación de un vehículo eléctrico, tiene lugar después de la fase I de modelización teórica.
Por lo general, el procedimiento se interrumpe durante cierto tiempo entre la fase I y la fase II.
Fase III de almacenamiento
Esta fase III opcional de almacenamiento de los datos medidos/calculados durante la fase II se realiza preferiblemente en forma simultánea a esta fase II.
De acuerdo con una variante, estos datos pueden (tele)transmitirse por cualquier medio adecuado a un centro de análisis y tratamiento.
Los datos recogidos in situ pueden utilizarse para mejorar la modelización.
Dispositivo de aplicación
i. Cargador/descargador: se trata de un elemento capaz de cargar y descargar la batería de forma programable. En un ejemplo, este cargador/descargador es un módulo de evaluación FRDM-BC3770-EVM de la empresa NXP. La corriente máxima de carga es de 2 A y la corriente máxima de descarga es de 1 A. El cargador está equipado con un microcontrolador que controla el perfil de corriente. Además, el microcontrolador dispone de convertidores analógicodigitales de 16 bits, que permiten realizar mediciones de alta resolución.
ii. Sensores de variables reales.
Tbat: el sensor utilizado para medir la temperatura de la batería es ventajosamente un termistor que permite medir la temperatura con alta resolución. Por ejemplo, podría ser un termistor de 10 kü conectado al módulo de evaluación BC3770. Este termistor tiene un convertidor analógico-digital de 16 bits y ofrece una resolución de temperatura de 0,01 °C con un ruido aceptable.
El módulo de evaluación también mide Vbat, así como la corriente de las baterías. Se trata de funciones nativas del módulo. Tamb se mide utilizando un termistor conectado a los convertidores analógico-digitales del microcontrolador. iii. Registrador de datos
iv. Generador de señales de corriente de carga
v. Unidad central de control y cálculo
Estos 3 elementos iii. iv. y v. pueden reunirse en un nano-ordenador monoplaca con un procesador ARM, como la raspberry pi 3.
Ejemplo
El siguiente ejemplo ilustra una forma de realización preferida del procedimiento según la invención en una batería de iones litio.
Batería: batería cilíndrica de iones litio 18650 (diámetro: 18 mm, altura: 65 mm), de capacidad nominal de 3070 mAh. El dispositivo según la invención utilizado en este ejemplo es el descrito con anterioridad.
La figura 3 muestra el banco de pruebas que comprende este dispositivo y la batería.
Este banco de pruebas está colocado en una cámara ambiental para poder controlar la temperatura del aire alrededor de la batería y evitar perturbaciones innecesarias.
El termistor se pone en contacto con la batería utilizando un fluido que favorece el contacto térmico, en este ejemplo, grasa de silicona comercializada por la sociedad VELLEMAN. También se aísla térmicamente del entorno exterior mediante caucho sintético.
Por lo tanto, se considera que la temperatura medida es la temperatura superficial de la batería.
Etapa (aü) estimación previa de los parámetros mCp y hA
En nuestro modelo termoeléctrico, los dos parámetros mCp y hA se consideran constantes con respecto al estado de carga y al estado de salud de la batería.
Estos dos parámetros se determinan de una vez por todas y se utilizan en los modelos eléctrico y térmico de una batería determinada.
Para estimar estos parámetros, se aplica una señal de onda cuadrada y media cero de ± 1 A con un período de 20 segundos. Esto elimina la influencia del calor generado por la variación de entropía, ya que a lo largo de un período el calor medio generado por AS es igual a cero.
La señal de onda cuadrada se aplica tras un período de relajación de la tensión de la batería de unos diez minutos. De este modo, se obtiene un valor OCV preciso. Esto deja solo dos parámetros desconocidos en la ecuación (3'): mCp y hA.
A continuación, se utiliza el banco de pruebas descrito con anterioridad y el algoritmo RLS para estimar estos dos parámetros.
El gráfico superior de la figura 4 muestra la respuesta de la temperatura de la batería. El gráfico inferior muestra la respuesta de la tensión. Estas señales y el algoritmo RLS se pueden utilizar para estimar los parámetros deseados.
Las figuras 5 y 6 muestran la convergencia de estos dos parámetros. Se obtienen los siguientes valores: ■ mCp = 110 JK-1 ■ hA = 0,09 WK-1
Esta estimación se incorpora al modelo térmico. A continuación, se pueden estimar los demás parámetros.
Etapas (a) de carga y (b) de descarga
Se aplica una secuencia binaria pseudoaleatoria (SBPA) a la corriente de entrada de la batería para cargarla. Por ejemplo, se genera una señal SBpA entre 0A y 1A, como se muestra en la figura 7, y se repite. La batería alcanza una tasa C de 0,3.C. Se utiliza el mismo procedimiento para la descarga, y la señal Sd se sitúa entre 0A y -1A.
Etapa (e) de estimación periódica del modelo eléctrico
Se aplica una corriente como la descrita en el apartado anterior a la entrada de nuestro sistema de baterías. Se comienza descargando la batería, que inicialmente está cargada por completo. La respuesta de la tensión puede verse en la figura 8.
Los parámetros de interés en el modelo eléctrico gobernado por la ecuación 1 son Ro y OCV. Ro es un parámetro útil y se utiliza a menudo para estimar el envejecimiento de la batería. OCV es útil para el modelo térmico, ya que nos permitirá estimar mejor AS.
A continuación, se aplica el algoritmo RLS a los datos para estimar los parámetros del modelo. Los parámetros se actualizan cada 0,2 segundos. La figura 9 muestra la estimación de Ro en función del estado de carga. El estado de carga se define mediante el procedimiento de recuento Coulomb Se trata de seguir la corriente integrándola durante el uso de las baterías. Esta integración proporciona una indicación directa de la cantidad de carga eléctrica inyectada o extraída de la batería, lo que permite cuantificar con precisión el estado de carga de la batería. La figura 10 muestra también la estimación del OCV.
Etapa (g) de estimación periódica del modelo térmico
Como se muestra en la figura 2, el algoritmo RLS para la estimación periódica del modelo eléctrico funciona al mismo tiempo que el algoritmo RLS para la estimación periódica del modelo térmico. El primero, como ya se ha mencionado, estimará OCV a partir de la corriente y la tensión. Esta estimación es utilizada por un segundo algoritmo RLS basado en el modelo térmico y que utiliza la temperatura para estimar AS. Cabe señalar que, dado que la dinámica eléctrica es mucho más rápida que la dinámica térmica, la frecuencia utilizada para actualizar los parámetros eléctricos es mayor que para los parámetros térmicos. Así se obtiene la estimación de AS en la figura 11 en función del estado de carga de la batería. Un dato importante que hay que añadir es el tiempo que se tarda en obtener el perfil completo. Dado que la corriente media es de 0,5 A, se necesitaron 6 horas.
También puede obtenerse una estimación de la entalpía AH (véase la figura 12), ya que puede calcularse a partir de OCV y<a>S de la siguiente manera: AH = -F.OCV- TbatAS
Así pues, los datos termodinámicos (OCV, AS y AH) acaban de estimarse en tiempo real mediante un algoritmo recursivo. Ahora es necesario comparar estos perfiles estimados con los perfiles medidos por el procedimiento estándar, ETM (“electo-chemical thermodynamic measurement”).
Ejemplo de comparación con el procedimiento ETM
Banco de pruebas
Para realizar la comparación con el procedimiento estándar, se construye un banco de pruebas para llevar a cabo la medición. En la figura 13, se muestra un diagrama de este banco. Se compone principalmente de 9 elementos:
1. Una batería idéntica a la utilizada en el ejemplo según la invención.
2. Un cargador de batería idéntico al utilizado en el ejemplo según la invención.
3. Un termistor de 10 kü para medir la temperatura de la batería.
4. Una caja metálica que contiene la batería.
5. Un elemento Peltier en contacto con la carcasa para modificar la temperatura de la misma.
6. Un radiador para disipar el calor.
7. Un ventilador para enfriar el radiador.
8. Un controlador Peltier para controlar la temperatura de la batería.
9. Un Raspberry Pi 3 para controlar el cargador y el controlador Peltier (es decir, para controlar la temperatura de la batería). También se utiliza para registrar los datos de corriente, temperatura y tensión de la batería.
Medición de datos termodinámicos
Para medir AS, se utiliza la siguiente fórmula:
• AS es la variación de entropía
• F es la constante de Faraday
• OCV es la tensión en circuito abierto
• Tbates la temperatura de la batería
• x es el estado de carga
En la práctica, para medir AS, se coloca la batería en un estado de carga determinado, se desconecta la corriente y se esperan varias horas para que la tensión se relaje y tienda hacia OCV. A continuación, variamos la temperatura de la batería, lo que tendrá como efecto la variación de la tensión (esta variación es muy pequeña, por lo que la medida debe realizarse a alta resolución). Si no hay cambio de fase (cambio de la estructura cristalina en el interior del ánodo o del cátodo), existe una relación lineal entre Tbat y OCV. Solo queda determinar el coeficiente que une a ambos y multiplicarlo por la constante de Faraday para obtener AS. Para obtener un perfil completo, la medición debe repetirse en diferentes estados de carga.
OCV se obtiene directamente midiendo la tensión después de la relajación y AH se obtiene utilizando la ecuación: AH = -F.OCV -TbatAS.
Comparación
Ahora se comparan los perfiles obtenidos con el procedimiento estándar y el procedimiento en línea. En las figuras 14, 15 y 16, se comparan estos dos procedimientos. Se puede observar que existe un sesgo entre la OCV medida en línea según la invención y la medida después de la relajación según el ejemplo comparativo. Este sesgo se refleja entonces en la estimación de AH y AS.
Otro aspecto importante de la comparación es el tiempo de medición para obtener un perfil completo. Con el procedimiento ETM en el ejemplo comparativo, se requieren 5 días. Con el procedimiento según la invención: 6 horas. Esto supone un ahorro de más de un orden de magnitud.
Repetibilidad
Una característica importante a determinar es la repetibilidad del procedimiento según la invención. Para la misma batería, los perfiles termodinámicos de descarga se estimaron varias veces (7 veces). A continuación, se trazaron los perfiles medios y sus intervalos de confianza del 95 %. Esto puede verse en las figuras 17, 18 y 19.
Más que una compatibilidad perfecta con el perfil medido con el procedimiento estándar, lo importante es garantizar que la medición del perfil en línea según la invención sea repetible. En efecto, este perfil puede entonces utilizarse para estimar el estado de salud de la batería, y es su variación en el tiempo lo que es importante en este sentido. Del mismo modo, para estimar el estado de carga, lo importante es su repetibilidad.
Claims (9)
- REIVINDICACIONES 1. Procedimiento para medir, in situ, en línea y en tiempo real, los datos termodinámicos, en particular la variación de entropía AS, de al menos una batería, que consiste en: (Fase I) Realización de una modelización preliminar de la batería cuyo estado de carga (SOC: “State Of Charge”) está comprendido entre el 0 y el 100 % mediante la ejecución de las siguientes etapas: (a) cargar al menos parcialmente la batería con una señal Sc de corriente de carga; (b) y/o descargar al menos parcialmente la batería con una señal Sd de corriente de descarga; (c) medir variables reales útiles en las etapas siguientes; (d) modelizar el comportamiento eléctrico de la batería durante la carga (a) con una señal Sc de corriente de carga y/o la descarga (b) con una señal Sd de corriente de descarga, a fin de estimar los parámetros eléctricos de la batería; (e) estimar periódicamente los parámetros eléctricos de la batería con una frecuencia Fe; (f) modelizar el comportamiento térmico de la batería durante la carga (a) con una señal Sc de corriente de carga y/o la descarga (b) con una señal Sd de corriente de descarga, con el fin de estimar in situ, en línea y en tiempo real, al menos uno de los parámetros del modelo térmico, a saber, AS; (g) estimar periódicamente, a una frecuencia Fg, al menos uno de los parámetros del modelo térmico, en particular AS, utilizando al menos uno de los parámetros eléctricos estimados en la etapa (e); (Fase II) Medir los datos termodinámicos, en particular AS, de la batería en uso en una aplicación y en cualquier estado de carga, implementando los modelos eléctrico [etapa (d)] y térmico [etapa (f)] de la fase I, la estimación de los parámetros eléctricos [etapa (e)] y la estimación de al menos uno de los parámetros del modelo térmico, en particular AS [etapa (g)]; (Fase III) Almacenar opcionalmente los datos medidos/calculados durante la fase II y/o la fase I, en donde el procedimiento se caracteriza porque ^ la modelización de la etapa (d) consiste en considerar que la batería es un circuito eléctrico o modelo eléctrico que comprende una resistencia Ro, una tensión de circuito abierto OCV, y un circuito R<1>C<1>en serie; el comportamiento eléctrico de la batería se describe, en este modelo, mediante las siguientes ecuaciones:en donde Ui es la tensión a través del circuito RC, / es la corriente a través de la batería y Vbat es la tensión en los bornes de la batería. La ecuación (2) se discretiza como sigue: Vbat,k ~ Ikb0,k E¡k -l^ l.kE &i,k(^£Vfc-l —Vbat.k-1) EOCV^(2 )y se reescribe como: con:en donde:Ts es el período de muestreo de la señal Se; es el vector de parámetros; ^ La modelización de la etapa f) consiste en considerar la batería como un modelo térmico en el que la batería, por una parte, está sometida a una corriente de carga Sc susceptible de ser sometida a muestreo o a una corriente de descarga Sd susceptible de ser sometida a muestreo y, por otra parte, es un intercambiador de calor con su entorno y en el que el comportamiento térmico de la batería se describe mediante la siguiente ecuación:en donde: • m es la masa de la batería • Cp es la capacidad calorífica de la batería • Tbates la temperatura de la batería • t es la variable temporal • I es la corriente que circula por la batería • Vbat es la tensión en bornes de la batería • OCV es la tensión en circuito abierto de la batería • AS es la variación de entropía de la batería • F es la constante de Faraday • h es el coeficiente de intercambio de calor con el exterior • A es el área de la batería en contacto con el exterior • Tamb es la temperatura del ambiente exterior en donde la ecuación (3) se discretiza como sigue:y se reescribe como: con:en donde:Ts es el período de muestreo de la medición.
- 2. Procedimiento de acuerdo con la reivindicación 1,caracterizado porquelos productos mCp y hA son constantes y porque se estiman en una etapa (aü) previa a la etapa (a) que consiste esencialmente en: - aplicar una relajación de la batería para que el parámetro OCV de ecuación (2') del modelo eléctrico tenga un valor dado preciso; - aplicar una señal eléctrica de entrada Se periódica, cuyo período se elige de manera que la media del calor generado por AS durante un período sea del orden de ü, siendo Se preferiblemente una señal de onda cuadrada con un período de entre 10 y 30 segundos, idealmente del orden de 20 segundos, siendo esta señal de onda cuadrada más preferiblemente aún de media cero; - estimar los productos mCp y hA, preferiblemente mediante un algoritmo de mínimos cuadrados recursivos, utilizando la ecuación (3) que se convierte en la ecuación (3') siguiente: ( y midiendo las variables reales correspondientes a los parámetros VBat, I, Tbat, Tamb, - integrar esta estimación de los productos mCp y hA en el modelo térmico para las etapas (f) y (g).
- 3. Procedimiento de acuerdo con al menos una de las reivindicaciones anteriores,caracterizado porquelas variables reales medidas en la etapa (c) corresponden a los parámetros VBat, I, Tbat, Tamb de los modelos eléctrico y térmico.
- 4. Procedimiento de acuerdo con al menos una de las reivindicaciones anteriores,caracterizado porquela estimación según la etapa (e) se realiza mediante un algoritmo de mínimos cuadrados recursivos y porque la estimación según la etapa (g) se realiza mediante un algoritmo de mínimos cuadrados recursivos.
- 5. Procedimiento de acuerdo con la reivindicación 1, opcionalmente de acuerdo con al menos una de las reivindicaciones 2 a 4,caracterizado porque, para implementar la etapa de carga (a), se aplica repetidamente una corriente de carga Sc, cuyo intervalo de frecuencia está comprendido entre 0 y 1 Hz, preferiblemente una señal correspondiente a una Secuencia Binaria Pseudoaleatoria -SBPA-, eligiéndose dicha señal de forma que la “tasa C” de la batería esté comprendida entre 0,01.C y 3.C, preferiblemente entre 0,1.C y 2,5.C, y aún más preferiblemente entre 0,2.C y 2.C.
- 6. Procedimiento de acuerdo con la reivindicación 1, opcionalmente de acuerdo con al menos una de las reivindicaciones 2 a 5,caracterizado porque, para llevar a cabo la etapa de descarga b), se aplica repetidamente una corriente de descarga Sd, que es la corriente de entrada Se, cuyo intervalo de frecuencia está comprendido entre 0 y 1 Hz, preferiblemente una señal correspondiente a una Secuencia Binaria Pseudoaleatoria -SBPA-, eligiéndose dicha señal de forma que la “tasa C” de la batería esté comprendida entre 0,01.C y 3.C, preferiblemente entre 0,1.C y 2,5.C, y aún más preferiblemente entre 0,2.C y 2.C.
- 7. Procedimiento de acuerdo con la reivindicación 1, opcionalmente de acuerdo con al menos una de las reivindicaciones 2 a 6,caracterizado porquela variación de entalpia AH se estima a partir de OCV y AS mediante la siguiente ecuación: AH = —F.OCV — TbatAS(6 )
- 8. Procedimiento para determinar el estado de carga y/o el estado de salud de una batería a partir de AS y/o AH medido por el procedimiento de acuerdo con al menos una de las reivindicaciones anteriores.
- 9. Dispositivo configurado para implementar el procedimiento de acuerdo con al menos una de las reivindicaciones 2 a 8,caracterizado porquecomprende: i. al menos un cargador/descargador programable; ii. al menos un sensor de al menos una de las variables reales que corresponden a los siguientes parámetros: Veat, L Tbat Tamb; iii. al menos un registrador de datos; iv. al menos un generador de señales de corriente de carga; v. al menos una unidad central de control e informática capaz de controlar el o los elementos (i), (ii), (iii), (iv), para recoger y procesar los datos estimados en las etapas (ao), (e) y (g), preferiblemente mediante la aplicación de algoritmos de mínimos cuadrados recursivos.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1858919A FR3086802B1 (fr) | 2018-09-27 | 2018-09-27 | Procede et dispositif de mesure en temps reel et in situ des donnees thermodynamiques d'une batterie (enthalpie et entropie) |
| PCT/EP2019/076102 WO2020064959A1 (fr) | 2018-09-27 | 2019-09-26 | Procede et dispositif de mesure en temps reel et in situ des donnees thermodynamiques d'une batterie (enthalpie et entropie) |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2985237T3 true ES2985237T3 (es) | 2024-11-04 |
Family
ID=65243956
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES19782947T Active ES2985237T3 (es) | 2018-09-27 | 2019-09-26 | Procedimiento y dispositivo para la medición en tiempo real e in situ de los datos termodinámicos de una batería (entalpía y entropía) |
Country Status (11)
| Country | Link |
|---|---|
| US (1) | US12123915B2 (es) |
| EP (1) | EP3857639B1 (es) |
| JP (1) | JP7497345B2 (es) |
| KR (1) | KR102871414B1 (es) |
| CN (1) | CN112840496B (es) |
| CA (1) | CA3114194A1 (es) |
| ES (1) | ES2985237T3 (es) |
| FR (1) | FR3086802B1 (es) |
| HU (1) | HUE067374T2 (es) |
| PL (1) | PL3857639T3 (es) |
| WO (1) | WO2020064959A1 (es) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220344734A1 (en) * | 2021-04-14 | 2022-10-27 | Analog Devices, Inc. | Technique for estimation of internal battery temperature |
| CN113420471B (zh) * | 2021-06-01 | 2022-02-18 | 上海交通大学 | 基于电化学机理的动力锂电池热模型构建建立方法及系统 |
| US20250198958A1 (en) * | 2022-03-29 | 2025-06-19 | Temple University-Of The Commonwealth System Of Higher Education | System and method to detect mechanically damaged energy storage cells using electrical signals |
| CN114583301B (zh) * | 2022-04-29 | 2022-08-23 | 国网浙江省电力有限公司电力科学研究院 | 基于安全特征参量表征体系的电站热失控预警方法及系统 |
| CN115219939B (zh) * | 2022-09-15 | 2023-01-03 | 小米汽车科技有限公司 | 电池温度预测方法、装置、车辆和存储介质 |
| CN115774204B (zh) * | 2022-12-01 | 2026-04-07 | 章鱼博士智能技术(上海)有限公司 | 电池热失控异常的检测方法及电池管理单元 |
| FR3144868B1 (fr) | 2023-01-05 | 2024-11-29 | Entroview | Procédé automatique d'estimation de l'état de charge d'une cellule d'une batterie |
| CN115952624B (zh) * | 2023-03-10 | 2023-05-23 | 陕西空天信息技术有限公司 | 一种叶轮机械真实流场的cfd加速分析方法及系统 |
| KR102888925B1 (ko) * | 2024-01-17 | 2025-11-19 | 주식회사 엘지에너지솔루션 | 배터리 진단 장치 및 방법 |
| FR3166213A1 (fr) | 2024-09-10 | 2026-03-13 | Entroview | Procédé automatique d’estimation de la variation d’entropie d’une cellule d’une batterie |
| CN119936714B (zh) * | 2025-01-07 | 2025-11-04 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种用于锂离子电池不同快充工况下的健康状态估计方法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8446127B2 (en) | 2005-08-03 | 2013-05-21 | California Institute Of Technology | Methods for thermodynamic evaluation of battery state of health |
| US7595611B2 (en) | 2005-08-03 | 2009-09-29 | California Institute Of Technology | Electrochemical thermodynamic measurement system |
| US10556510B2 (en) | 2012-04-27 | 2020-02-11 | California Institute Of Technology | Accurate assessment of the state of charge of electrochemical cells |
| CN104471415B (zh) | 2012-04-27 | 2017-09-01 | 加州理工学院 | 用于电池应用的嵌入式芯片 |
| CN105301509B (zh) * | 2015-11-12 | 2019-03-29 | 清华大学 | 锂离子电池荷电状态、健康状态与功率状态的联合估计方法 |
| KR101805514B1 (ko) * | 2015-11-20 | 2017-12-07 | 한국과학기술원 | 배터리의 동적 엔트로피 추정 방법 |
| US10243385B2 (en) * | 2016-01-29 | 2019-03-26 | Robert Bosch Gmbh | Secondary battery management system |
| WO2017204750A1 (en) | 2016-05-27 | 2017-11-30 | Nanyang Technological University | Method of assessing a performance of an electrochemical cell, and apparatus thereof |
| US20180292465A1 (en) * | 2017-04-07 | 2018-10-11 | Board Of Regents, The University Of Texas System | Systems and methods for degradation analysis |
-
2018
- 2018-09-27 FR FR1858919A patent/FR3086802B1/fr active Active
-
2019
- 2019-09-26 KR KR1020217012446A patent/KR102871414B1/ko active Active
- 2019-09-26 CA CA3114194A patent/CA3114194A1/fr active Pending
- 2019-09-26 EP EP19782947.6A patent/EP3857639B1/fr active Active
- 2019-09-26 CN CN201980063847.5A patent/CN112840496B/zh active Active
- 2019-09-26 PL PL19782947.6T patent/PL3857639T3/pl unknown
- 2019-09-26 JP JP2021517449A patent/JP7497345B2/ja active Active
- 2019-09-26 WO PCT/EP2019/076102 patent/WO2020064959A1/fr not_active Ceased
- 2019-09-26 US US17/280,011 patent/US12123915B2/en active Active
- 2019-09-26 HU HUE19782947A patent/HUE067374T2/hu unknown
- 2019-09-26 ES ES19782947T patent/ES2985237T3/es active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112840496A (zh) | 2021-05-25 |
| PL3857639T3 (pl) | 2024-07-08 |
| JP7497345B2 (ja) | 2024-06-10 |
| KR102871414B1 (ko) | 2025-10-16 |
| KR20210096076A (ko) | 2021-08-04 |
| EP3857639A1 (fr) | 2021-08-04 |
| US20220034973A1 (en) | 2022-02-03 |
| JP2022503861A (ja) | 2022-01-12 |
| EP3857639B1 (fr) | 2024-02-21 |
| HUE067374T2 (hu) | 2024-10-28 |
| FR3086802A1 (fr) | 2020-04-03 |
| US12123915B2 (en) | 2024-10-22 |
| WO2020064959A1 (fr) | 2020-04-02 |
| FR3086802B1 (fr) | 2021-12-24 |
| CN112840496B (zh) | 2024-05-31 |
| CA3114194A1 (fr) | 2020-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| ES2985237T3 (es) | Procedimiento y dispositivo para la medición en tiempo real e in situ de los datos termodinámicos de una batería (entalpía y entropía) | |
| CN103502829B (zh) | 用于电化学存储系统热管理的优化方法 | |
| US9660299B2 (en) | Strain measurement based battery testing | |
| Troxler et al. | The effect of thermal gradients on the performance of lithium-ion batteries | |
| KR102501672B1 (ko) | 배터리 팩에서의 벌크 힘 및 충전 상태 추정에 대한 벌크 힘의 적용 | |
| US10556510B2 (en) | Accurate assessment of the state of charge of electrochemical cells | |
| Williard et al. | Comparative analysis of features for determining state of health in lithium-ion batteries | |
| Haifeng et al. | A new SOH prediction concept for the power lithium-ion battery used on HEVs | |
| KR101805514B1 (ko) | 배터리의 동적 엔트로피 추정 방법 | |
| KR102527327B1 (ko) | 배터리의 온도 추정 장치 및 방법, 배터리 관리 장치 및 방법 | |
| US9599584B2 (en) | Imbedded chip for battery applications | |
| Mingant et al. | EIS measurements for determining the SOC and SOH of Li-ion batteries | |
| KR102654635B1 (ko) | 배터리의 용량상태 결정방법 및 그를 이용한 용량감소 측정방법 | |
| KR101329915B1 (ko) | 하이브리드 전기 자동차용 배터리의 주행 중 최대 출력파워 예측 방법 및 장치 | |
| Zhang et al. | Investigation of current sharing and heat dissipation in parallel-connected lithium-ion battery packs | |
| Millet et al. | Isothermal calorimeter heat measurements of a 20Ah lithium iron phosphate battery cell | |
| Eddahech et al. | Lithium-ion battery heat generation investigation based on calorimetric entropy measurements | |
| Tripathy et al. | Internal temperature prediction of Lithium-ion cell using differential voltage technique | |
| Gorny et al. | Monitoring and Health prognosis of Lithium-Ion battery system | |
| Stanciu et al. | Performance degradation of thermal parameters during cycle ageing of high energy density ni-mn-co based lithium-ion battery cells | |
| Stroe et al. | Electric circuit modeling of lithium-sulfur batteries during discharging state | |
| Kellner et al. | Electrical and thermal behavior of pouch-format lithium ion battery cells under high-performance and standard automotive duty-cycles | |
| Ovejas et al. | SoH evaluation of LiFePO 4 cells using impedance and thermal measurements | |
| Paramadayalan et al. | Experimental and Numerical Studies on Battery Temperature Distribution | |
| Novak et al. | State of Health Estimation of Lithium-Ion Batteries Based on Incremental Capacity and Pulse Analysis |