CN116184246A - SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment - Google Patents

SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment Download PDF

Info

Publication number
CN116184246A
CN116184246A CN202310162497.6A CN202310162497A CN116184246A CN 116184246 A CN116184246 A CN 116184246A CN 202310162497 A CN202310162497 A CN 202310162497A CN 116184246 A CN116184246 A CN 116184246A
Authority
CN
China
Prior art keywords
redox flow
vanadium redox
flow battery
soc
soh
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.)
Pending
Application number
CN202310162497.6A
Other languages
Chinese (zh)
Inventor
高翔
宋浩
郑成航
吴卫红
韩爱
余鸿敏
王丁振
刘少俊
张霄
徐甸
宋思慧
张悠
张成健
岑可法
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202310162497.6A priority Critical patent/CN116184246A/en
Publication of CN116184246A publication Critical patent/CN116184246A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/39—Circuit design at the physical level
    • G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
    • 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/378—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC] specially adapted for the type of battery or accumulator
    • 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
    • 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/392—Determining battery ageing or deterioration, e.g. state of health
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2115/00—Details relating to the type of the circuit
    • G06F2115/02—System on chip [SoC] design
    • 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/30—Hydrogen technology
    • Y02E60/50—Fuel cells

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Fuel Cell (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)

Abstract

The invention provides a method for jointly estimating SOC and SOH of an all-vanadium redox flow battery, which comprises the following steps: s1, constructing an equivalent circuit model, and determining parameters of the equivalent circuit model as activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltage V s The method comprises the steps of carrying out a first treatment on the surface of the S2, constructing parameters based on multiple physical fields for calculating parameters in equivalent circuit modelCalculating a model; s3, substituting the material performance parameter values of the all-vanadium redox flow battery into a parameter calculation model to obtain an all-vanadium redox flow battery simulation model with the same or similar material performance; s4, inputting specific size parameter values of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery. The invention improves the universality of the equivalent circuit model on all-vanadium redox flow batteries of various types, can estimate the SOC and SOH of the positive half battery and the negative half battery respectively, and can monitor the actual running condition of the all-vanadium redox flow battery more comprehensively.

Description

SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment
Technical Field
The invention relates to a state estimation method of an all-vanadium redox flow battery, in particular to a method and a device for jointly estimating SOC and SOH of the all-vanadium redox flow battery and terminal equipment, and belongs to the technical field of all-vanadium redox flow batteries.
Background
The all-vanadium redox flow battery (VRB) is an efficient electrochemical energy storage device, has the advantages of large energy storage scale, long service life and the like, and can effectively inhibit the influence of new energy power generation fluctuation on a power grid.
The electrical energy storage medium of the VRB is an electrolyte containing vanadium ions, which is stored in a reservoir external to the battery. Parameters such as state of charge (SOC) and state of health (SOH) of the VRB are important indicators during battery operation. The traditional SOC estimation method is mainly realized by constructing an equivalent circuit model and adopting an ampere-hour integration method. The method has two defects, namely, element parameters of an equivalent circuit model are required to be obtained through a large number of VRB experimental tests, new experimental test results are required to be supplemented for updating the element parameters when the VRB model is replaced, and therefore the method is poor in universality; secondly, the equivalent circuit model can only estimate the state parameter of the SOC and can not estimate the SOH of the battery, and the SOH mainly provides reference for operations such as mixing of positive and negative poles and recovery of battery capacity, and the like, and is also one of important indexes in the battery operation process, and the SOH can not be estimated only but also can not comprehensively reflect real-time monitoring of the VRB operation state. Therefore, aiming at all-vanadium redox flow batteries in different scenes, a general method is developed for monitoring the running states of VRB of various models in real time, and the method has very important significance for prolonging the service life of the VRB and promoting the large-scale commercial use of the VRB.
Disclosure of Invention
Based on the background, the invention aims to provide the SOC and SOH combined estimation method of the all-vanadium redox flow battery, which can estimate and calculate the SOC and SOH of the positive half battery and the negative half battery of all-vanadium redox flow batteries of various types respectively without fitting parameters of an equivalent circuit model through experimental data, and can monitor the running condition of the all-vanadium redox flow battery more comprehensively.
In order to achieve the above object, the present invention provides the following technical solutions:
the combined estimation method of the SOC and the SOH of the all-vanadium redox flow battery comprises the following steps:
s1, constructing an equivalent circuit model, and determining parameters of the equivalent circuit model as activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltageV s ；
S2, constructing a parameter calculation model based on multiple physical fields for calculating parameters in an equivalent circuit model, and correspondingly calculating continuous equation expression of the change of the concentration of each vanadium ion in a galvanic pile and a storage tank, and calculating the activation polarization overpotential V act And concentration polarization overpotential V con Equation expression of the reaction kinetics of (2) for calculating ohmic loss R ohm Equation expression of (2), calculate open circuit voltage V s Equation expression of (2) and calculating pump loss current I pump Momentum equation expression of (a);
s3, substituting the material performance parameter value of the all-vanadium redox flow battery into a parameter calculation model;
s4, inputting specific size parameter values of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery.
Preferably, in step S1, the equivalent circuit model includes a voltage source, a variable resistor, a controlled current source and two controlled voltage sources, wherein the voltage source, the variable resistor and the two controlled voltage sources are connected in series and then connected in parallel with the controlled current source, and the controlled voltage source is used for representing an activation polarization overpotential V caused by electrochemical reaction during operation of the vanadium redox flow battery act Another controlled voltage source is used for representing concentration polarization overpotential V caused by electrochemical reaction in the operation process of the all-vanadium redox flow battery con The controlled current source is used for representing pump loss current I pump The variable resistor is used for representing ohmic loss R caused by a pile and electrolyte ohm The voltage source is used for representing the open circuit voltage V of the all-vanadium redox flow battery s 。
Preferably, in step S2, the continuous equation expression for calculating the change in the concentration of each valence vanadium ion in the galvanic pile and the storage tank is,
Figure SMS_1
Figure SMS_2
in the formula (1) and the formula (2),
Figure SMS_3
respectively representing the concentration of the i-valent vanadium ions in the galvanic pile and the storage tank; m is the serial number of single cells; v (V) cell ,V tank Respectively representing the electrolyte volumes in the single cell and the single storage tank; d (D) i Represents the transmembrane diffusion coefficient of the i-valent ion; thick (thick) ed ,thick mem The thickness of the electrode and the film, respectively; z is the electron transfer number; f is Faraday constant; i is the charge and discharge current of the battery, and Q is the electrolyte flow of the battery;
calculating the activation polarization overpotential V act And concentration polarization overpotential V con The expression of the reaction kinetic equation of (c) is as follows,
Figure SMS_4
Figure SMS_5
in the formula (3) and the formula (4), gamma 0p ,γ 0n Respectively positive and negative standard reaction rate constants; e (E) 0p ,E 0n Respectively positive and negative half-cell reaction standard potentials; a is that ed The specific surface area of the electrode;
Figure SMS_6
represents the electrode surface diffusion coefficient of I or j vanadium ions, and when charged (I>0) I4,j takes 3 and when discharging (I<0) i, 5,j, 2; ρ, μ are electrolyte density and viscosity, respectively; d, d ed Is the average fiber diameter of the electrode; a is that cell Is the sectional area of the pile;
calculation of ohmic loss R ohm The equation expression of (2) is that,
Figure SMS_7
in formula (5), the thick et Representing half cell thickness; sigma (sigma) ed ,σ mem ,σ et Respectively representing the conductivities of the electrode, the membrane and the electrolyte;
calculating open circuit voltage V s The equation expression of (2) is that,
Figure SMS_8
calculating pump loss current I pump The expression of the momentum equation of (c) is that,
Figure SMS_9
in the formula (7), epsilon is the porosity of the electrode; lambda (lambda) CK A Carman-Kozeny constant that is a fibrous material; l (L) pipe ,d pipe ,h pipe The length, the diameter and the height of the half-cell pipeline are respectively; f is the total secondary loss coefficient; v pipe Is the flow rate of the electrolyte in the pipeline; η (eta) pump Is the efficiency of the pump; u (U) d Is the VRB system terminal voltage.
When the traditional all-vanadium redox flow battery equivalent circuit model is used for simulation, the circuit element parameters of the traditional all-vanadium redox flow battery equivalent circuit model need to be identified on the basis of a large number of experimental tests, and the element parameters of the equivalent circuit model need to be identified again in the face of all-vanadium redox flow batteries of different models, and the element parameters in the equivalent circuit model are calculated by utilizing the parameter calculation model constructed based on the expression, so that the identification of the circuit element parameters based on the experimental tests can be avoided, and the concentration condition of vanadium ions of various valences can be calculated in real time.
Preferably, in step S3, after substituting the material performance parameter value of the all-vanadium redox flow battery into the parameter calculation model, an equation expression with a general fixed parameter value with the same or similar material performance is obtained.
Preferably, in step S4, estimating the SOC and SOH of the monitored vanadium redox flow battery includes the steps of:
equation expressions for calculating SOC and SOH are constructed,
Figure SMS_10
Figure SMS_11
in the formula (14) and the formula (15), Q practical And Q theoretical Respectively representing the actual available capacity and the theoretical available capacity; c v,p And c v,n Respectively represents the initial vanadium ion concentration of the positive and negative electrolyte, U d The end voltage of the vanadium redox flow battery is obtained.
An SOC and SOH joint estimation device of an all-vanadium redox flow battery, comprising:
a first construction module for constructing an equivalent circuit model, determining the parameters of the equivalent circuit model as the activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltage V s ；
A second construction module for constructing a parameter calculation model based on multiple physical fields for calculating parameters in the equivalent circuit model, and a continuous equation expression for correspondingly calculating the change of the concentration of each valence vanadium ion in the galvanic pile and the storage tank, and calculating the activation polarization overpotential V act And concentration polarization overpotential V con Equation expression of the reaction kinetics of (2) for calculating ohmic loss R ohm Equation expression of (2), calculate open circuit voltage V s Equation expression of (2) and calculating pump loss current I pump Momentum equation expression of (a);
the first calculation module is used for substituting the general fixed parameter value of the all-vanadium redox flow battery into the parameter calculation model;
and the second calculation module is used for inputting the specific size parameter value of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery.
A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method as claimed in any one of the preceding claims when the computer program is executed.
A computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method as claimed in any one of the preceding claims.
Compared with the prior art, the invention has the following advantages:
according to the SOC and SOH joint estimation method of the all-vanadium redox flow battery, disclosed by the invention, parameters in an equivalent circuit model are calculated through a parameter calculation model based on multiple physical fields and corresponding equation expressions such as a reaction dynamics equation, a continuous equation and a momentum equation, and the like, so that the parameters are not required to be simulated through a large amount of experimental data, and the universality of the equivalent circuit model on all-vanadium redox flow batteries of various types is improved; the method can estimate the SOC and SOH of the positive half-cell and the negative half-cell respectively, and can monitor the actual running condition of the all-vanadium redox flow battery more comprehensively.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present invention, and that other drawings can be obtained according to the provided drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of a method for joint estimation of SOC and SOH of an all-vanadium redox flow battery of the present invention;
FIG. 2 is a schematic diagram of an equivalent circuit model and a parameter calculation model in accordance with the present invention;
FIG. 3 is a charge-discharge curve obtained by a simulation test and an experimental test of the combined estimation method of the SOC and the SOH of the all-vanadium redox flow battery;
fig. 4 is real-time estimation data of the SOC of the all-vanadium redox flow battery obtained by the simulation test of the SOC and SOH joint estimation method of the all-vanadium redox flow battery.
Detailed Description
The technical scheme of the invention is further specifically described below through specific embodiments and with reference to the accompanying drawings. It should be understood that the practice of the invention is not limited to the following examples, but is intended to be within the scope of the invention in any form and/or modification thereof.
In the present invention, unless otherwise specified, all parts and percentages are by weight, and the equipment, materials, etc. used are commercially available or are conventional in the art. The methods in the following examples are conventional in the art unless otherwise specified. The components and devices in the following examples are, unless otherwise indicated, all those components and devices known to those skilled in the art, and their structures and principles are known to those skilled in the art from technical manuals or by routine experimentation.
The embodiment of the invention discloses a combined estimation method of SOC and SOH of an all-vanadium redox flow battery, which utilizes a parameter calculation model based on multiple physical fields to calculate parameters of an equivalent circuit model, and forms a simulation model aiming at the monitored all-vanadium redox flow battery with a specific model by inputting specific dimension parameter values of the monitored all-vanadium redox flow battery, so as to simulate a charging and discharging curve of the actual monitored all-vanadium redox flow battery and estimate real-time SOC and SOH, thereby realizing state monitoring of the monitored all-vanadium redox flow battery. According to the method, parameters of an equivalent circuit model are prevented from being obtained through a large number of experiments, and the SOC and SOH of the positive half battery and the negative half battery of the all-vanadium redox flow battery are estimated respectively.
In the following detailed description, which is given by way of example of an all-vanadium flow battery of 5kW/30kWh, and with reference to the accompanying drawings, numerous specific details are set forth in order to provide a thorough understanding of embodiments of the invention. However, one or more embodiments may be practiced by one of ordinary skill in the art without these specific details.
The method for jointly estimating the SOC and the SOH of the all-vanadium redox flow battery shown in FIG. 1 comprises the following steps:
s1, constructing an equivalent circuit model, and determining the equivalent circuit modelThe parameter being the activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltage V s ；
S2, constructing a parameter calculation model based on multiple physical fields for calculating parameters in an equivalent circuit model, and correspondingly calculating continuous equation expression of the change of the concentration of each vanadium ion in a galvanic pile and a storage tank, and calculating the activation polarization overpotential V act And concentration polarization overpotential V con Equation expression of the reaction kinetics of (2) for calculating ohmic loss R ohm Equation expression of (2), calculate open circuit voltage V s Equation expression of (2) and calculating pump loss current I pump Momentum equation expression of (a);
s3, substituting the general fixed parameter value of the all-vanadium redox flow battery into a parameter calculation model;
s4, inputting specific size parameter values of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery.
The equivalent circuit model includes a voltage source, a variable resistor, a controlled current source and two controlled voltage sources, wherein the voltage source, the variable resistor and the two controlled voltage sources are connected in series and then connected in parallel with the controlled current source, as shown in fig. 2. A controlled voltage source is used for representing activation polarization overpotential V caused by electrochemical reaction in the operation process of the all-vanadium redox flow battery act Another controlled voltage source is used for representing concentration polarization overpotential V caused by electrochemical reaction in the operation process of the all-vanadium redox flow battery con The controlled current source is used for representing the pump loss current I pump Variable resistor for characterizing ohmic losses R caused by galvanic pile and electrolyte ohm The voltage source is used for representing open circuit voltage V of all-vanadium redox flow battery s 。
Each parameter in the equivalent circuit model constructed in step S1 is obtained by calculation of a parameter calculation model based on multiple physical fields, and the expression of the correlation equation used for the calculation is specifically as follows.
The continuous equation expression for calculating the change of the concentration of each valence vanadium ion in the galvanic pile and the storage tank is as follows,
Figure SMS_12
Figure SMS_13
in the formula (1) and the formula (2),
Figure SMS_14
respectively representing the concentration of the i-valent vanadium ions in the galvanic pile and the storage tank; m is the serial number of single cells; v (V) cell ,V tank Respectively representing the electrolyte volumes in the single cell and the single storage tank; d (D) i Represents the transmembrane diffusion coefficient of the i-valent ion; thick (thick) ed ,thick mem The thickness of the electrode and the film, respectively; z is the electron transfer number; f is Faraday constant; i is the charge and discharge current of the battery, and Q is the electrolyte flow of the battery;
calculating the activation polarization overpotential V act And concentration polarization overpotential V con The expression of the reaction kinetic equation of (c) is as follows,
Figure SMS_15
Figure SMS_16
in the formula (3) and the formula (4), gamma 0p ,γ 0n Respectively positive and negative standard reaction rate constants; e (E) 0p ,E 0n Respectively positive and negative half-cell reaction standard potentials; a is that ed The specific surface area of the electrode;
Figure SMS_17
represents the electrode surface diffusion coefficient of I or j vanadium ions, and when charged (I>0) I4,j takes 3 and when discharging (I<0) i, 5,j, 2; ρ, μ are electrolyte density and viscosity, respectively; d, d ed Is the average fiber diameter of the electrode; a is that cell Is cut off for galvanic pileAn area;
calculation of ohmic loss R ohm The equation expression of (2) is that,
Figure SMS_18
in formula (5), the thick et Representing half cell thickness; sigma (sigma) ed ,σ mem ,σ et Respectively representing the conductivities of the electrode, the membrane and the electrolyte;
calculating open circuit voltage V s The equation expression of (2) is that,
Figure SMS_19
calculating pump loss current I pump The expression of the momentum equation of (c) is that,
Figure SMS_20
in the formula (7), epsilon is the porosity of the electrode; lambda (lambda) CK A Carman-Kozeny constant that is a fibrous material; l (L) pipe ,d pipe ,h pipe The length, the diameter and the height of the half-cell pipeline are respectively; f is the total secondary loss coefficient; v pipe Is the flow rate of the electrolyte in the pipeline; η (eta) pump Is the efficiency of the pump; u (U) d Is the VRB system terminal voltage.
The material performance parameter values of the all-vanadium redox flow battery are shown in table 1.
TABLE 1 Material Performance parameter values for vanadium redox flow batteries
Figure SMS_21
Figure SMS_22
After substituting the material performance parameter values of the all-vanadium redox flow battery in table 1 into the parameter calculation model, the following equation expression of the general fixed parameter values with the same or similar material performance is obtained:
the continuous equation expression for calculating the change of the concentration of each valence vanadium ion in the galvanic pile and the storage tank is as follows,
Figure SMS_23
calculating the activation polarization overpotential V act And concentration polarization overpotential V con The expression of the reaction kinetic equation of (c) is as follows,
Figure SMS_24
Figure SMS_25
calculation of ohmic loss R ohm The equation expression of (2) is that,
R ohm ＝4.17×10 -4 M(29)
calculating open circuit voltage V s The equation expression of (2) is that,
Figure SMS_26
calculating pump loss current I pump The expression of the momentum equation of (c) is that,
Figure SMS_27
the specific size parameter values of the monitored all-vanadium redox flow battery are input to a parameter calculation model with general fixed parameter values of the same or similar material properties, and in this embodiment, the specific size parameter values of the monitored all-vanadium redox flow battery are shown in table 2.
TABLE 2 specific dimensional parameter values for all-vanadium redox flow batteries
Figure SMS_28
Figure SMS_29
Inputting the parameter values in the table 2 into a parameter calculation model with general fixed parameter values, and obtaining the model of the all-vanadium redox flow battery suitable for the monitored model.
Estimating the SOC and SOH of the monitored all-vanadium redox flow battery comprises the steps of constructing an equation expression for calculating the SOC and SOH,
Figure SMS_30
Figure SMS_31
in the formula (14) and the formula (15), Q practical And Q theoretical Respectively representing the actual available capacity and the theoretical available capacity; c v, p and c v, n represents the initial vanadium ion concentration of the positive and negative electrode electrolyte, U d The end voltage of the vanadium redox flow battery is obtained.
Simulation and test are carried out under the conditions that the temperature T is 25 ℃, the charge-discharge current I is 80A, the electrolyte flow Q is 32L/min, and the effectiveness of the method is verified by comparing simulation data with experimental data. The charge-discharge curve of the monitored vanadium redox flow battery drawn according to the simulation data and the experimental data is shown in fig. 3, the SOC of the monitored vanadium redox flow battery obtained through simulation changes with time as shown in fig. 4, and the SOH of the monitored vanadium redox flow battery obtained through simulation is 92.399%.
According to the test, the simulation estimation method of the SOC and SOH of the all-vanadium redox flow battery is used for carrying out simulation estimation on the monitored all-vanadium redox flow battery, and compared with an actual experiment, the method is basically consistent with the charge-discharge curve of the monitored all-vanadium redox flow battery, and the accuracy of the method is verified. The method can estimate the SOC and SOH of the positive half battery and the negative half battery respectively, so that the condition of the vanadium redox flow battery in the actual running process can be monitored more comprehensively.
The embodiment of the invention also discloses a device for jointly estimating the SOC and the SOH of the all-vanadium redox flow battery, which comprises a first construction module, a second construction module, a first calculation module and a second calculation module. The first construction module is used for constructing an equivalent circuit model, and determining the parameters of the equivalent circuit model as the activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltage V s . The second construction module is used for constructing a parameter calculation model based on multiple physical fields for calculating parameters in the equivalent circuit model, and a continuous equation expression for correspondingly calculating the change of the concentration of each valence vanadium ion in the galvanic pile and the storage tank, and calculating the activation polarization overpotential V act And concentration polarization overpotential V con Equation expression of the reaction kinetics of (2) for calculating ohmic loss R ohm Equation expression of (2), calculate open circuit voltage V s Equation expression of (2) and calculating pump loss current I pump Momentum equation expression of (c). The first calculation module is used for substituting the general fixed parameter value of the all-vanadium redox flow battery into the parameter calculation model. The second calculation module is used for inputting specific size parameter values of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery.
The embodiment of the invention also discloses a terminal device which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor executes the computer program to realize the steps of the method.
Embodiments of the present invention also disclose a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method as described in any of the above.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to facilitate an understanding of the method of the present invention and its core ideas. It should be noted that it will be apparent to those skilled in the art that various modifications and adaptations of the invention can be made without departing from the principles of the invention and these modifications and adaptations are intended to be within the scope of the invention as defined in the following claims.

Claims (8)

1. A combined estimation method of SOC and SOH of an all-vanadium redox flow battery is characterized by comprising the following steps of: the method comprises the following steps:
s1, constructing an equivalent circuit model, and determining parameters of the equivalent circuit model as activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltage V s ；
S2, constructing a parameter calculation model based on multiple physical fields for calculating parameters in an equivalent circuit model, and correspondingly calculating continuous equation expression of the change of the concentration of each vanadium ion in a galvanic pile and a storage tank, and calculating the activation polarization overpotential V act And concentration polarization overpotential V con Equation expression of the reaction kinetics of (2) for calculating ohmic loss R ohm Equation expression of (2), calculate open circuit voltage V s Equation expression of (2) and calculating pump loss current I pump Momentum equation expression of (a);
s3, substituting the material performance parameter value of the all-vanadium redox flow battery into a parameter calculation model;
s4, inputting specific size parameter values of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery.
2. The method for jointly estimating the SOC and the SOH of the all-vanadium redox flow battery according to claim 1, which is characterized by comprising the following steps of: in step S1, the equivalent circuit model includes a voltage source, a variable resistor, a controlled current source, and two controlled voltage sources, where the voltage source, the variable resistor, and the two controlled voltage sources are connected in series and then connected in parallel to the controlled current source, and the controlled voltage source is used to characterize an activated polarization overpotential V caused by an electrochemical reaction during operation of the all-vanadium redox flow battery act Another one of the controlled voltage sources is used for representing the all-vanadium redox flow batteryConcentration polarization overpotential V caused by electrochemical reaction during operation con The controlled current source is used for representing pump loss current I pump The variable resistor is used for representing ohmic loss R caused by a pile and electrolyte ohm The voltage source is used for representing the open circuit voltage V of the all-vanadium redox flow battery s 。
3. The method for jointly estimating the SOC and the SOH of the all-vanadium redox flow battery according to claim 1, which is characterized by comprising the following steps of: in step S2, the continuous equation expression for calculating the change in the concentration of each valence vanadium ion in the electric pile and the storage tank is,
Figure FDA0004094660550000011
Figure FDA0004094660550000012
in the formula (1) and the formula (2),
Figure FDA0004094660550000013
respectively representing the concentration of the i-valent vanadium ions in the galvanic pile and the storage tank; m is the serial number of single cells; v (V) cell ,V tank Respectively representing the electrolyte volumes in the single cell and the single storage tank; d (D) i Represents the transmembrane diffusion coefficient of the i-valent ion; thick (thick) ed ,thick mem The thickness of the electrode and the film, respectively; z is the electron transfer number; f is Faraday constant; i is the charge and discharge current of the battery, and Q is the electrolyte flow of the battery;
calculating the activation polarization overpotential V act And concentration polarization overpotential V con The expression of the reaction kinetic equation of (c) is as follows,
Figure FDA0004094660550000021
Figure FDA0004094660550000022
in the formula (3) and the formula (4), gamma 0p ,γ 0n Respectively positive and negative standard reaction rate constants; e (E) 0p ,E 0n Respectively positive and negative half-cell reaction standard potentials; a is that ed The specific surface area of the electrode;
Figure FDA0004094660550000023
represents the electrode surface diffusion coefficient of I or j vanadium ions, and when charged (I>0) I4,j takes 3 and when discharging (I<0) i, 5,j, 2; ρ, μ are electrolyte density and viscosity, respectively; d, d ed Is the average fiber diameter of the electrode; a is that cell Is the sectional area of the pile;
calculation of ohmic loss R ohm The equation expression of (2) is that,
Figure FDA0004094660550000024
in formula (5), the thick et Representing half cell thickness; sigma (sigma) ed ,σ mem ,σ et Respectively representing the conductivities of the electrode, the membrane and the electrolyte;
calculating open circuit voltage V s The equation expression of (2) is that,
Figure FDA0004094660550000025
calculating pump loss current I pump The expression of the momentum equation of (c) is that,
Figure FDA0004094660550000026
in the formula (7), ζ is an electrode porosity; lambda (lambda) CK A Carman-Kozeny constant that is a fibrous material; l (L) pipe ,d pipe ,h pipe Respectively half-cellsLength, diameter and height of the pipe; f is the total secondary loss coefficient; v pipe Is the flow rate of the electrolyte in the pipeline; η (eta) pump Is the efficiency of the pump; u (U) d Is the VRB system terminal voltage.
4. The method for jointly estimating the SOC and the SOH of the all-vanadium redox flow battery according to claim 1, which is characterized by comprising the following steps of: in step S3, after substituting the material performance parameter value of the all-vanadium redox flow battery into the parameter calculation model, an equation expression with the same or similar material performance and general fixed parameter value is obtained.
5. The method for jointly estimating the SOC and the SOH of the all-vanadium redox flow battery according to claim 1, which is characterized by comprising the following steps of: in step S4, estimating the SOC and SOH of the monitored all-vanadium redox flow battery includes the steps of:
equation expressions for calculating SOC and SOH are constructed,
Figure FDA0004094660550000031
Figure FDA0004094660550000032
in the formula (14) and the formula (15), Q practical And Q theoretical Respectively representing the actual available capacity and the theoretical available capacity; c v,p And c v,n Respectively represents the initial vanadium ion concentration of the positive and negative electrolyte, U d The end voltage of the vanadium redox flow battery is obtained.
6. The utility model provides an SOC and SOH joint estimation device of all vanadium redox flow battery which characterized in that: the SOC and SOH joint estimation device of the all-vanadium redox flow battery comprises:
a first construction module for constructing an equivalent circuit model, determining the parameters of the equivalent circuit model as the activation polarization overpotential V act Concentration polarization overpotential V con Pump loss current I pump Ohmic loss R ohm And open circuit voltage V s ；
A second construction module for constructing a parameter calculation model based on multiple physical fields for calculating parameters in the equivalent circuit model, and a continuous equation expression for correspondingly calculating the change of the concentration of each valence vanadium ion in the galvanic pile and the storage tank, and calculating the activation polarization overpotential V act And concentration polarization overpotential V con Equation expression of the reaction kinetics of (2) for calculating ohmic loss R ohm Equation expression of (2), calculate open circuit voltage V s Equation expression of (2) and calculating pump loss current I pump Momentum equation expression of (a);
the first calculation module is used for substituting the general fixed parameter value of the all-vanadium redox flow battery into the parameter calculation model;
and the second calculation module is used for inputting the specific size parameter value of the monitored all-vanadium redox flow battery into the parameter calculation model in the step S3, and estimating the SOC and SOH of the monitored all-vanadium redox flow battery.
7. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: the processor, when executing the computer program, implements the steps of the method according to any one of claims 1-5.
8. A computer-readable storage medium storing a computer program, characterized in that: the computer program implementing the steps of the method according to any of claims 1-5 when executed by a processor.
CN202310162497.6A 2023-02-24 2023-02-24 SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment Pending CN116184246A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310162497.6A CN116184246A (en) 2023-02-24 2023-02-24 SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310162497.6A CN116184246A (en) 2023-02-24 2023-02-24 SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment

Publications (1)

Publication Number Publication Date
CN116184246A true CN116184246A (en) 2023-05-30

Family

ID=86444121

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310162497.6A Pending CN116184246A (en) 2023-02-24 2023-02-24 SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment

Country Status (1)

Country Link
CN (1) CN116184246A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117239186A (en) * 2023-11-10 2023-12-15 杭州德海艾科能源科技有限公司 SOH (solid-state oxygen) online evaluation device and method for all-vanadium redox flow battery system
CN118070714A (en) * 2024-04-19 2024-05-24 艾博特瑞能源科技(苏州)有限公司 Simulation method of all-vanadium redox flow battery based on different charge-discharge strategies

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117239186A (en) * 2023-11-10 2023-12-15 杭州德海艾科能源科技有限公司 SOH (solid-state oxygen) online evaluation device and method for all-vanadium redox flow battery system
CN117239186B (en) * 2023-11-10 2024-02-06 杭州德海艾科能源科技有限公司 SOH (solid-state oxygen) online evaluation device and method for all-vanadium redox flow battery system
CN118070714A (en) * 2024-04-19 2024-05-24 艾博特瑞能源科技(苏州)有限公司 Simulation method of all-vanadium redox flow battery based on different charge-discharge strategies
CN118070714B (en) * 2024-04-19 2024-06-28 艾博特瑞能源科技(苏州)有限公司 Simulation method of all-vanadium redox flow battery based on different charge-discharge strategies

Similar Documents

Publication Publication Date Title
CN107145628B (en) Method for predicting cycle life of lithium battery based on electrochemical-thermal coupling model
Lyu et al. A lead-acid battery's remaining useful life prediction by using electrochemical model in the Particle Filtering framework
CN116184246A (en) SOC and SOH joint estimation method and device of all-vanadium redox flow battery and terminal equipment
CN107066722B (en) Electrochemical model-based combined estimation method for state of charge and state of health of power battery system
Lamorgese et al. Electrochemical-thermal P2D aging model of a LiCoO2/graphite cell: Capacity fade simulations
CN105932349B (en) A kind of lithium ion battery long-life fast charge method
CN108761341A (en) A kind of lithium ion battery battery chemical modeling parameter acquisition methods
Leng et al. A practical framework of electrical based online state-of-charge estimation of lithium ion batteries
CN115632179B (en) Intelligent quick charging method and system for lithium ion battery
KR101282687B1 (en) Method for estimating battery SOC of vehicle
CN106198293B (en) Method for testing moisture content in power lithium battery
Merei et al. Multi-physics model for a vanadium redox flow battery
Khaki et al. Voltage loss and capacity fade reduction in vanadium redox battery by electrolyte flow control
CN105891715A (en) Lithium ion battery health state estimation method
Khaki et al. Multi-objective optimal charging current and flow management of Vanadium Redox Flow Batteries for fast charging and energy-efficient operation
CN114609523A (en) Online battery capacity detection method, electronic equipment and storage medium
CN109164397B (en) Lithium battery life loss assessment method considering charging rate and ambient temperature
Zhang et al. The state-of-charge estimation of supercapacitor with kalman filtering algorithm
Zhou et al. Tab design based on the internal distributed properties in a Zinc–Nickel single-flow battery
US20170207475A1 (en) Electrode for redox flow battery, redox flow battery, and electrode characteristics evaluation method
CN116721706A (en) Battery simulation method, storage medium and electronic equipment
CN112447266B (en) All-vanadium redox flow battery capacity attenuation modeling method considering water molecule migration
Guo The application of genetic algorithms to parameter estimation in lead-acid battery equivalent circuit models
Pranata et al. Static and dynamic characteristic lead acid flow battery
Badrinarayanan et al. Temperature dependence of capacity decay due to ion diffusion in vanadium redox flow battery

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination