WO2023130776A1 - Method and system for predicting working condition health status of battery in energy storage power station - Google Patents

Method and system for predicting working condition health status of battery in energy storage power station Download PDF

Info

Publication number
WO2023130776A1
WO2023130776A1 PCT/CN2022/121940 CN2022121940W WO2023130776A1 WO 2023130776 A1 WO2023130776 A1 WO 2023130776A1 CN 2022121940 W CN2022121940 W CN 2022121940W WO 2023130776 A1 WO2023130776 A1 WO 2023130776A1
Authority
WO
WIPO (PCT)
Prior art keywords
battery
data
original
signal
energy storage
Prior art date
Application number
PCT/CN2022/121940
Other languages
French (fr)
Chinese (zh)
Inventor
张雪松
林达
钱平
赵波
杨帆
章雷其
马瑜涵
葛晓慧
刘敏
陈凌宇
陈哲
Original Assignee
国网浙江省电力有限公司电力科学研究院
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 国网浙江省电力有限公司电力科学研究院 filed Critical 国网浙江省电力有限公司电力科学研究院
Publication of WO2023130776A1 publication Critical patent/WO2023130776A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/392Determining battery ageing or deterioration, e.g. state of health
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/396Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02WCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
    • Y02W30/00Technologies for solid waste management
    • Y02W30/50Reuse, recycling or recovery technologies
    • Y02W30/84Recycling of batteries or fuel cells

Definitions

  • the invention relates to the technical field of energy storage batteries, in particular to a method and system for predicting the health state of batteries in energy storage power stations.
  • the current SOH prediction of batteries generally uses the historical data of battery decline to fit the empirical decline model or some commonly used data-driven models. This method can easily capture the law of battery decline, and some better results in the scene. This method requires a large amount of relatively complete battery decline historical data to fit the model, which is generally a linear or exponential decline process, and is not suitable for the superimposition of periodic factors (such as temperature, etc.). Higher, generally it is relatively complete charge and discharge data, otherwise it is necessary to estimate the SOH of the battery first and then use the estimated value to realize the SOH prediction, resulting in error superposition.
  • the present invention provides a method and system for predicting the health status of batteries in energy storage power stations, which approximates by intercepting the capacity (energy) of the platform area The full interval capacity (energy) of the battery, while ensuring a certain accuracy, reduces the demand for battery data for battery health status evaluation. It is suitable for SOH evaluation of energy storage systems working under peak load conditions, and solves the health status of energy storage batteries in complex environments.
  • the corresponding signal decomposition algorithm is used to realize the decoupling of the complex factors of the battery SOH sequence, so as to pay attention to the essential law of battery decline and improve the accuracy of battery health state prediction and evaluation.
  • a technical solution adopted by the present invention is as follows: a method for predicting the health status of batteries in energy storage power stations, which includes:
  • Step 1 Obtain the historical operation data of the battery of the energy storage power station, and preprocess the original data
  • Step 2 According to the characteristics of the battery, intercept the platform voltage range of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform range for a period of time from the charging period of historical operating data, and obtain the charging capacity value each time through the integration of ampere hours , to obtain the original battery health state change sequence;
  • Step 3 Use the empirical mode decomposition algorithm to separate the battery cycle change segment and the linear decay segment from the original battery health state change sequence
  • Step 4 Establish the data-driven models for the periodical change segment and the linear decline segment respectively, and realize the forecast of the change trend of the battery linear decline segment and the seasonal cycle segment.
  • the battery platform area contains about 80% of the capacity (energy) of the battery.
  • the platform area and the two fast polarization areas of the battery can be easily distinguished.
  • the annual power station operation data is processed to obtain the daily charging amount of the battery cell in the platform area, and then the signal decomposition method is used to obtain the linear decline period and season of the battery In order to realize the decoupling of battery health status and temperature influence, obtain the approximate linear decline path of the battery, so as to realize the prediction of battery SOH based on the field data of the power station.
  • the historical battery operation data of the energy storage power station includes: battery pack voltage data, battery pack current data, battery cell voltage data, battery cell current data, battery cell temperature data, battery module temperature data and the accumulated running time of the energy storage power station; among them, battery cell voltage data, battery cell current data and battery cell temperature data are necessary data;
  • the data preprocessing using data interpolation to unify the sampling interval of the data; removing outliers, interpolating and smoothing the data whose sampling interval has been unified.
  • step 2 the plateau voltage range of the battery is defined by the battery capacity incremental curve, and the voltage platform is required to include all peaks of the battery capacity incremental curve.
  • step 2 the battery capacity increment curve is obtained by the following method:
  • dV is the voltage interval
  • dQ is the capacity corresponding to the equal voltage interval of the battery
  • IC is the battery capacity increment.
  • step 2 based on the battery IC curve, the voltage interval of the battery platform area is quantitatively fixed to obtain the reference value of the platform area capacity;
  • the capacity sequence of the platform area of the charging section of the battery is obtained by using the method, which is the original battery health state change sequence.
  • step 3 is as follows:
  • the upper envelope and the lower envelope of the original signal are constructed by means of the curve fitting method; the original signal is taken from the original The battery health state change sequence;
  • the original battery health state change sequence is decomposed into a number of empirical modal components and residual sequences, in which the residual sequence is mostly the noise of the signal, which is ignored; some empirical modal component sequences include some periodic signal sequences and Approximate linearly decaying signal sequence.
  • step 4 Furthermore, the specific content of step 4 is as follows:
  • the modeling methods include linear regression, support vector machine and Gaussian process regression, and the square sine kernel function is used for the quasi-periodic signal sequence
  • the Gaussian process regression characterizes the periodic variation characteristics of periodic signals
  • the data of the linear decay process and the periodic process are respectively selected as the training set and cross-validation is turned on;
  • Another technical solution adopted by the present invention is: a system for predicting the health status of batteries in energy storage power stations, which includes:
  • Data acquisition and preprocessing unit acquire the historical operation data of the battery of the energy storage power station, and preprocess the original data
  • Original battery health state change sequence acquisition unit According to the characteristics of the battery, intercept the battery platform voltage range, use this as a reference to obtain a period of time from the historical operating data charging section voltage and current data of the battery charging platform interval, and obtain it through ampere-hour integration Get the original battery health state change sequence for each charging capacity value;
  • Separation unit Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state change sequence;
  • Battery change trend prediction unit establish data-driven models for the cycle change segment and the linear decline segment respectively, and realize the change trend prediction of the battery linear decline segment and the seasonal cycle segment.
  • the present invention approximates the capacity (energy) of the entire battery interval by intercepting the capacity (energy) of the platform area, which greatly reduces the demand for battery data for battery health status evaluation while ensuring a certain accuracy, and is suitable for storage batteries working under peak-shaving conditions.
  • Energy system SOH evaluation solve the problem of multi-modal superposition of energy storage battery health status in a complex environment, and use the corresponding signal decomposition algorithm to realize the decoupling of complex factors in the battery SOH sequence, so as to pay attention to the essential law of battery decline and improve battery health Accuracy of state prediction and assessment.
  • Fig. 1 is a schematic flow chart of the method for predicting the health state of the battery working condition of the energy storage power station according to the present invention
  • Fig. 2 is the battery IC graph in the application example of the present invention.
  • Fig. 3-4 is the comparison chart before and after working condition data processing in the application example of the present invention (Fig. 3 is the schematic diagram before working condition data processing, and Fig. 4 is the schematic diagram after working condition data processing);
  • Fig. 5 is the original sequence diagram of the platform capacity of all monomers in the battery cluster in the application example of the present invention.
  • Fig. 6 is a comparison chart between the predicted value and the actual value of the battery cell capacity in the application example of the present invention.
  • Fig. 7 is a structural diagram of the system for predicting the health state of the battery in the energy storage power station according to the present invention.
  • this embodiment provides a method for predicting the health state of a battery in an energy storage power station, the steps of which are as follows:
  • Step 1 Obtain the historical operation data of the battery of the energy storage power station, specific to the battery cell, with a sampling accuracy of 1min to 5min, and preprocess the original data, including data interpolation and data cleaning.
  • Step 2 According to the characteristics of the battery, intercept the platform voltage range of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform range for a period of time from the charging period of historical operating data, and obtain the charging capacity value each time through the integration of ampere hours , to obtain the original battery health state change sequence;
  • Step 3 Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state sequence
  • Step 4 Establish the data-driven models for the periodical change segment and the linear decline segment respectively, and realize the forecast of the change trend of the battery linear decline segment and the seasonal cycle segment.
  • Battery pack voltage data, battery pack current data, battery cell voltage data, battery cell current data, battery cell temperature data, battery module temperature data and the accumulated time of system operation are necessary data.
  • Data preprocessing due to poor sampling accuracy or short-term failure of the monitoring system, some sampled data may appear outliers or missing data to varying degrees, resulting in data distortion, thereby affecting subsequent data extraction. Therefore, it is necessary to delete outlier data and interpolate missing data in this step to obtain usable data that can be further processed.
  • the power station does not participate in peak shaving every day, so the battery of the power station does not work or only works in a small DOD interval at certain times. This part of the data is not representative and it is difficult to extract new information from it, so this part of the data will be Clean up and make empty.
  • the platform voltage range of the battery described in step 2 is obtained through the battery capacity increment curve, and the voltage platform is required to include all peaks of the capacity increment curve.
  • the capacity increment curve of the battery is obtained by the following method:
  • the accumulative charging ampere-hours of each sampling point in the charging process of the battery is obtained by the ampere-hour integration method, and the battery voltage-capacity curve is obtained.
  • dV is the voltage interval
  • dQ is the capacity corresponding to the equal voltage interval of the battery
  • IC is the battery capacity increment sequence.
  • the voltage interval of the plateau area of the battery can be fixed quantitatively, and the reference value for obtaining the capacity of the plateau area can be obtained.
  • the data of each battery charging section is intercepted, and the capacity sequence of the platform area of the charging section of the battery is obtained by the ampere-hour integration method, which is the reference sequence of the battery SOH.
  • the empirical mode decomposition (EMD) algorithm is used to decompose the original sequence, and the steps of the EMD algorithm are as follows:
  • the original sequence of battery decline is decomposed into several empirical mode components and residuals, among which the residual sequence is mostly the noise of the signal and can be ignored.
  • Several empirical mode components include local characteristic signals of different time scales, mainly including some quasi-periodic signals and approximate linear decay signals.
  • Gaussian process regression can well describe the periodic variation characteristics of periodic signals.
  • Gaussian process regression is used to model the two sequences:
  • this embodiment provides a system for predicting the state of health of a battery in an energy storage power station, which includes:
  • Data acquisition and preprocessing unit acquire the historical operation data of the battery of the energy storage power station, and preprocess the original data
  • Original battery health state change sequence acquisition unit According to the characteristics of the battery, intercept the battery platform voltage range, use this as a reference to obtain a period of time from the historical operating data charging section voltage and current data of the battery charging platform interval, and obtain it through ampere-hour integration Get the original battery health state change sequence for each charging capacity value;
  • Separation unit Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state change sequence;
  • Battery change trend prediction unit establish data-driven models for the cycle change segment and the linear decline segment respectively, and realize the change trend prediction of the battery linear decline segment and the seasonal cycle segment.
  • the historical battery operation data of the energy storage power station includes: battery pack voltage data, battery pack current data, battery cell voltage data, battery cell current data, battery cell temperature data, Battery module temperature data and the accumulated time of energy storage power station operation;
  • the data preprocessing using data interpolation to unify the sampling interval of the data; removing outliers, interpolating and smoothing the data whose sampling interval has been unified.
  • the plateau voltage range of the battery is defined by the battery capacity increment curve, and the voltage plateau is required to include all peaks of the battery capacity increment curve.
  • the battery capacity increment curve is obtained by the following method:
  • dV is the voltage interval
  • dQ is the capacity corresponding to the equal voltage interval of the battery
  • IC is the battery capacity increment.
  • the capacity sequence is the original battery health state change sequence.
  • the upper envelope and the lower envelope of the original signal are constructed by means of the curve fitting method; the original signal is taken from the original The battery health state change sequence;
  • the original battery health state change sequence is decomposed into a number of empirical modal components and residual sequences, in which the residual sequence is mostly the noise of the signal, which is ignored; some empirical modal component sequences include some periodic signal sequences and Approximate linearly decaying signal sequence.
  • the specific content of the battery change trend prediction unit is as follows:
  • the modeling methods include linear regression, support vector machine and Gaussian process regression, and the square sine kernel function is used for the quasi-periodic signal sequence
  • the Gaussian process regression characterizes the periodic variation characteristics of periodic signals
  • the data of the linear decay process and the periodic process are respectively selected as the training set and cross-validation is turned on;
  • Example 1 Refer to Example 1 for the parts not described in detail in Example 2.
  • the data of the application example comes from the operation data of a demonstration project using battery energy storage in an echelon.
  • the data recording period is 1 year.
  • the specific real-time steps are as follows:
  • Step 1 Obtain the historical operation data of the energy storage power station, specific to the battery cell, with a sampling accuracy of 1min to 5min, and perform preprocessing on the original data, including data interpolation and data cleaning.
  • the data sampling points are first unified, and the measure taken is to interpolate the sparsely sampled data and unify the data sampling interval to 1min. Due to the poor quality of the original voltage data due to insufficient sampling accuracy, outlier points, interpolation, and smoothing are performed on the data that has been uniformly sampled to obtain a smoother voltage curve of the battery.
  • the battery voltage curve before and after processing is shown in Figure 2 .
  • Step 2 according to the characteristics of the battery, intercept the platform voltage range of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform range for a period of time from the charging period of historical operating data, and obtain the charging capacity value each time through the ampere-hour integration , to obtain the battery health state change sequence.
  • the voltage interval of the battery platform area is defined by the capacity increment curve of the battery.
  • a lithium iron phosphate battery of the same batch as the power station battery is taken, and the constant current charging data of the battery at a certain rate is obtained to draw a battery capacity increment curve.
  • dV takes a fixed value of 0.002V
  • dQ corresponding to the voltage interval can be obtained by using the difference method
  • the capacity increment curve of the battery is obtained according to the definition of the capacity increment curve, as shown in Figure 3-4.
  • the range of the plateau area can be easily determined through the capacity increment curve, and the voltage range of the plateau area is limited to 3.25-3.4V in this embodiment.
  • Based on the obtained voltage interval intercept the daily charging section data and integrate the current to obtain the battery platform area capacity sequence.
  • the normalized battery capacity sequence operation is not performed to make the result more intuitive, and a certain battery cluster is obtained.
  • the original sequences of all the monomers in are shown in Figure 5.
  • Step 3 the empirical mode decomposition method is used to decompose the historical data of the battery, and the mode decomposition can be realized by using the EMD algorithm toolkit in the signal processing toolbox of Matlab software.
  • the maximum IMF The number is set to 3. At this point, the post-signal decomposition and decoupling work is completed, and the battery linear segment and cycle segment sequence are obtained.
  • Step 4 Establish data-driven models for the periodical variation segment and the linear recession segment respectively, so as to realize the prediction of the change trend of the battery linear recession segment and the seasonal cycle segment.
  • the selected kernel function is a square exponential kernel function that can characterize the battery's long-term decline trend.
  • the expression of the kernel function is:
  • ⁇ , l and ⁇ are hyperparameters. Since the analytical expression of the kernel function contains several hyperparameters, selecting appropriate parameter values is crucial for establishing an accurate Gaussian process regression model.
  • x represents the input of training data
  • y represents the output of training data
  • the posterior estimated value of the hyperparameters of the kernel function can be obtained according to the training data.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Secondary Cells (AREA)
  • Tests Of Electric Status Of Batteries (AREA)

Abstract

Disclosed are a method and system for predicting the working condition health status of a battery in an energy storage power station. The prediction method of the present invention comprises: acquiring historical operation data of a battery in an energy storage power station, and pre-processing the original data; intercepting a platform voltage interval of the battery according to the characteristics of the battery, using the platform voltage interval as a reference to acquire, from a historical operation data charging section, charging platform interval voltage and current data of the battery in a period of time, acquiring a charging capacity value for each instance by means of an ampere-hour integral, and acquiring an original battery health status change sequence; separating a periodic change section and a linear degradation section of the battery from the original battery health status change sequence by using an empirical mode decomposition algorithm; establishing data-driven models for the periodic change section and the linear degradation section, respectively, and predicting the change trend of the linear degradation section and a seasonal periodic section of the battery. In the present invention, the demand during battery health status evaluation on battery data is reduced while precision is ensured to a certain extent, and the precision of battery health status prediction and evaluation is increased.

Description

储能电站电池工况健康状态预测方法及系统Method and system for predicting the health status of batteries in energy storage power stations 技术领域technical field
本发明涉及储能电池技术领域,特别是一种储能电站电池工况健康状态预测方法及预测系统。The invention relates to the technical field of energy storage batteries, in particular to a method and system for predicting the health state of batteries in energy storage power stations.
背景技术Background technique
随着电化学储能技术的发展,锂离子电池的成本逐年下降,除动力电池外,以磷酸铁锂电池为代表的大量储能电池被投放到电化学储能领域。同时电动汽车产业发展也带来了动力电池到了退役潮,大量的退役动力电池一方面被直接拆解回收,一些性能较好的退役电池被分选后进入梯次利用阶段,这些梯次电池也往往用于用户侧储能特别是调峰场景。储能电池投运装机量一般较大,对电池健康状态进行评估和预测对系统的健康稳定运行至关重要,准确的状态评估和趋势预测有助于发现问题电池方便系统维护,保证了储能系统的健康稳定。With the development of electrochemical energy storage technology, the cost of lithium-ion batteries has been decreasing year by year. In addition to power batteries, a large number of energy storage batteries represented by lithium iron phosphate batteries have been put into the field of electrochemical energy storage. At the same time, the development of the electric vehicle industry has also brought about a wave of decommissioning of power batteries. On the one hand, a large number of decommissioned power batteries are directly dismantled and recycled, and some decommissioned batteries with better performance are sorted and enter the cascade utilization stage. These cascade batteries are often used It is used for energy storage on the user side, especially in peak shaving scenarios. The installed capacity of energy storage batteries is generally large, and evaluating and predicting the health status of batteries is crucial to the healthy and stable operation of the system. Accurate status assessment and trend prediction help to find problems with batteries, facilitate system maintenance, and ensure energy storage The health and stability of the system.
常规的数据驱动SOH预测方法:当前对电池的SOH预测一般采用直接利用电池衰退历史数据来拟合经验衰退模型或一些常用的数据驱动模型,该方法可以很容易捕捉到电池衰退的规律,某些场景下取得较好的效果。该方法需要大量的较为完整的电池衰退历史数据进行模型的拟合,一般为线性衰退或指数衰退过程,并不适用于周期影响因素(如温度等)叠加的情况,同时对电池衰退历史数据要求较高,一般是较为完整的充放电数据,否则需要先对电池进行SOH估计再利用估计值实现SOH预测,造成误差叠加。Conventional data-driven SOH prediction method: The current SOH prediction of batteries generally uses the historical data of battery decline to fit the empirical decline model or some commonly used data-driven models. This method can easily capture the law of battery decline, and some better results in the scene. This method requires a large amount of relatively complete battery decline historical data to fit the model, which is generally a linear or exponential decline process, and is not suitable for the superimposition of periodic factors (such as temperature, etc.). Higher, generally it is relatively complete charge and discharge data, otherwise it is necessary to estimate the SOH of the battery first and then use the estimated value to realize the SOH prediction, resulting in error superposition.
发明内容Contents of the invention
针对电池充放电数据不完整电池容量(能量)真值难以直接获取的问题,本发明提供一种储能电站电池工况健康状态预测方法及预测系统,其通过截取平台区容量(能量)来近似电池全区间容量(能量),在保证一定精度的同时降低电池健康状态评估对电池数据的需求,适用于工作在调峰工况下的储能系统SOH评估,解决复杂环境下储能电池健康状态多模态叠加的问题,采 用相应的信号分解算法实现电池SOH序列复杂因素的解耦合,从而关注电池衰退的本质规律,提高电池健康状态预测和评估的精度。Aiming at the problem that it is difficult to directly obtain the true value of battery capacity (energy) due to incomplete battery charging and discharging data, the present invention provides a method and system for predicting the health status of batteries in energy storage power stations, which approximates by intercepting the capacity (energy) of the platform area The full interval capacity (energy) of the battery, while ensuring a certain accuracy, reduces the demand for battery data for battery health status evaluation. It is suitable for SOH evaluation of energy storage systems working under peak load conditions, and solves the health status of energy storage batteries in complex environments. For the problem of multi-modal superposition, the corresponding signal decomposition algorithm is used to realize the decoupling of the complex factors of the battery SOH sequence, so as to pay attention to the essential law of battery decline and improve the accuracy of battery health state prediction and evaluation.
为此,本发明采用的一种技术方案如下:储能电站电池工况健康状态预测方法,其包括:For this reason, a technical solution adopted by the present invention is as follows: a method for predicting the health status of batteries in energy storage power stations, which includes:
步骤1:获取储能电站电池历史运行数据,对原始数据进行预处理;Step 1: Obtain the historical operation data of the battery of the energy storage power station, and preprocess the original data;
步骤2:根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取原始的电池健康状态变化序列;Step 2: According to the characteristics of the battery, intercept the platform voltage range of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform range for a period of time from the charging period of historical operating data, and obtain the charging capacity value each time through the integration of ampere hours , to obtain the original battery health state change sequence;
步骤3:使用经验模态分解算法从原始的电池健康状态变化序列中分离出电池周期变化段和线性衰退段;Step 3: Use the empirical mode decomposition algorithm to separate the battery cycle change segment and the linear decay segment from the original battery health state change sequence;
步骤4:分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Step 4: Establish the data-driven models for the periodical change segment and the linear decline segment respectively, and realize the forecast of the change trend of the battery linear decline segment and the seasonal cycle segment.
针对电池充放电数据不完整电池容量(能量)真值难以直接获取的问题,提出采用电池平台区的容量(能量)作为电池的健康状态参考值。其中电池平台区包含了电池的80%左右的容量(能量),对于储能磷酸铁锂电池曲线,平台区和电池的两段的快速极化区可以很容易进行区分。Aiming at the problem that it is difficult to directly obtain the true value of battery capacity (energy) due to incomplete battery charge and discharge data, it is proposed to use the capacity (energy) of the battery platform area as the reference value of the battery's health status. The battery platform area contains about 80% of the capacity (energy) of the battery. For the energy storage lithium iron phosphate battery curve, the platform area and the two fast polarization areas of the battery can be easily distinguished.
由于现场电池电压、电流采样频率较低的问题,常用的基于容量增量曲线特征、差分电压特征的方法由于精度问题基本不可用,为了尽量减小采样精度对特征计量值的影响,采用电流积分方法获取电池平台区累计安时数,以此作为电池的容量健康状态的参考依据。Due to the low frequency of on-site battery voltage and current sampling, the commonly used methods based on capacity increment curve characteristics and differential voltage characteristics are basically unavailable due to accuracy problems. In order to minimize the impact of sampling accuracy on characteristic measurement values, current integration is used. The method obtains the accumulative ampere-hours in the battery platform area, which is used as a reference for the capacity and health status of the battery.
针对在无温控或弱温控条件下的电池健康状态受温度影响较大,对年度的电站运行数据处理得到电芯每天的平台区充电量后再利用信号分解方法获取电池线性衰退段和季节性周期段,从而实现电池健康状态和温度影响的解耦,获取电池近似线性衰退路径,从而实现基于电站现场数据的电池SOH的预测。In view of the fact that the health status of the battery under the condition of no temperature control or weak temperature control is greatly affected by the temperature, the annual power station operation data is processed to obtain the daily charging amount of the battery cell in the platform area, and then the signal decomposition method is used to obtain the linear decline period and season of the battery In order to realize the decoupling of battery health status and temperature influence, obtain the approximate linear decline path of the battery, so as to realize the prediction of battery SOH based on the field data of the power station.
进一步地,步骤1中,所述储能电站电池历史运行数据包括:电池组电压数据、电池组电流数据、电池单体电压数据、电池单体电流数据、电池单体温度数据、电池模块温度数据及储能电站运行的累计时间;其中,电池单体电压数据、电池单体电流数据和电池单体温度数据是必要数据;Further, in step 1, the historical battery operation data of the energy storage power station includes: battery pack voltage data, battery pack current data, battery cell voltage data, battery cell current data, battery cell temperature data, battery module temperature data and the accumulated running time of the energy storage power station; among them, battery cell voltage data, battery cell current data and battery cell temperature data are necessary data;
所述的数据预处理:采用数据插值,将数据的采样间隔统一;对已经统一采样间隔的数据进行去除离群点、插补和平滑处理。The data preprocessing: using data interpolation to unify the sampling interval of the data; removing outliers, interpolating and smoothing the data whose sampling interval has been unified.
进一步地,步骤2中,电池的平台电压区间通过电池容量增量曲线进行界定,要求电压平台包括电池容量增量曲线的所有峰。Further, in step 2, the plateau voltage range of the battery is defined by the battery capacity incremental curve, and the voltage platform is required to include all peaks of the battery capacity incremental curve.
更进一步地,步骤2中,所述电池容量增量曲线由以下方法获取:Furthermore, in step 2, the battery capacity increment curve is obtained by the following method:
获取电池在某充电倍率下的完整充电段电压、电流数据;Obtain the voltage and current data of the complete charging section of the battery at a certain charging rate;
通过安时积分法获取电池在充电过程中各个采样点的累计充电安时数得到电池电压-容量曲线;Obtain the cumulative charging ampere-hours of each sampling point during the charging process by the ampere-hour integration method to obtain the battery voltage-capacity curve;
设置合理的电压间隔dV,通过线性插值方法获取按固定电压间隔增加的电压-容量曲线,随后通过差分方法获取电池IC曲线数据,如下式:Set a reasonable voltage interval dV, obtain the voltage-capacity curve increasing at a fixed voltage interval through the linear interpolation method, and then obtain the battery IC curve data through the differential method, as follows:
Figure PCTCN2022121940-appb-000001
Figure PCTCN2022121940-appb-000001
上式中,dV为电压间隔,dQ为电池等电压间隔对应的容量,IC为电池容量增量。In the above formula, dV is the voltage interval, dQ is the capacity corresponding to the equal voltage interval of the battery, and IC is the battery capacity increment.
再进一步地,步骤2中,基于电池IC曲线定量固定电池平台区的电压区间,得到获取平台区容量的基准值;基于平台区容量的基准值,截取每一次电池充电段数据,通过安时积分法获取电池的充电段平台区容量序列,即为原始的电池健康状态变化序列。Furthermore, in step 2, based on the battery IC curve, the voltage interval of the battery platform area is quantitatively fixed to obtain the reference value of the platform area capacity; The capacity sequence of the platform area of the charging section of the battery is obtained by using the method, which is the original battery health state change sequence.
进一步地,步骤3的具体步骤如下:Further, the specific steps of step 3 are as follows:
a.根据从原始信号中找出的全部局部极大值和极小值点,借助曲线拟合方法,构造出原始信号的上包络线和下包络线;所述的原始信号取自原始的电池健康状态变化序列;a. According to all the local maximum and minimum points found from the original signal, the upper envelope and the lower envelope of the original signal are constructed by means of the curve fitting method; the original signal is taken from the original The battery health state change sequence;
b.求上、下包络线的均值并构造均值曲线,然后用原始信号减去该曲线,得到信号的中间分量;b. Find the mean value of the upper and lower envelopes and construct the mean value curve, and then subtract the curve from the original signal to obtain the intermediate component of the signal;
c.判断所得的信号中间分量是否满足本征模函数IMF的约束条件:如果满足,则该分量就是一个IMF;如果不满足,则以该分量为待分解信号重复步骤a至步骤b,直到满足IMF的约束条件,或筛分门限值小于预设门限值;c. Determine whether the intermediate component of the obtained signal satisfies the constraints of the intrinsic mode function IMF: if it is satisfied, the component is an IMF; if it is not satisfied, repeat steps a to b with the component as the signal to be decomposed until it is satisfied The constraint condition of IMF, or the screening threshold value is less than the preset threshold value;
d.通过上述步骤得到第一个IMF后,用原始信号与该第一个IMF相减,作为新的原始信号重复步骤a至步骤c,迭代直到残差分量满足预设条件,剩余信号分量即为残差分量,迭代结束;d. After the first IMF is obtained through the above steps, subtract the original signal from the first IMF, repeat steps a to c as a new original signal, iterate until the residual component meets the preset condition, and the remaining signal component is is the residual component, and the iteration ends;
经过以上分解过程,原始的电池健康状态变化序列被分解为若干经验模态分量及残差序列,其中残差序列多是信号的噪声,忽略;若干经验模态分量序列包括一些类周期信号序列及近似线性衰退信号序列。After the above decomposition process, the original battery health state change sequence is decomposed into a number of empirical modal components and residual sequences, in which the residual sequence is mostly the noise of the signal, which is ignored; some empirical modal component sequences include some periodic signal sequences and Approximate linearly decaying signal sequence.
更进一步地,步骤4的具体内容如下:Furthermore, the specific content of step 4 is as follows:
基于获取的类周期信号序列和近似线性衰退信号序列分别建模,对近似线性衰退信号序列建模的方式包括线性回归、支持向量机和高斯过程回归,对类周期信号序列,采用平方正弦核函数的高斯过程回归刻画类周期信号的周期变化特性;Based on the obtained quasi-periodic signal sequence and approximate linear decay signal sequence, the modeling methods include linear regression, support vector machine and Gaussian process regression, and the square sine kernel function is used for the quasi-periodic signal sequence The Gaussian process regression characterizes the periodic variation characteristics of periodic signals;
采用高斯过程回归建模的步骤如下:The steps for modeling with Gaussian process regression are as follows:
a.根据经验模态分解算法分解结果,分别选取线性衰退过程和周期过程数据作为训练集并开启交叉验证;a. According to the decomposition results of the empirical mode decomposition algorithm, the data of the linear decay process and the periodic process are respectively selected as the training set and cross-validation is turned on;
b.根据不同的训练集数据特征,选取不同的核函数,并初始化超参数;b. According to different training set data characteristics, select different kernel functions and initialize hyperparameters;
c.使用训练集数据训练回归模型,采用极大似然估计获取超参数最优估计得到后验模型;c. Use the training set data to train the regression model, and use the maximum likelihood estimation to obtain the optimal estimate of the hyperparameters to obtain the posterior model;
d.利用训练后的模型对输入的测试集数据进行预测,并将预测值和验证集数据进行对比分析模型预测能力。d. Use the trained model to predict the input test set data, and compare the predicted value with the verification set data to analyze the predictive ability of the model.
本发明采用的另一种技术方案为:储能电站电池工况健康状态预测系统,其包括:Another technical solution adopted by the present invention is: a system for predicting the health status of batteries in energy storage power stations, which includes:
数据获取及预处理单元:获取储能电站电池历史运行数据,对原始数据进行预处理;Data acquisition and preprocessing unit: acquire the historical operation data of the battery of the energy storage power station, and preprocess the original data;
原始电池健康状态变化序列获取单元:根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取原始的电池健康状态变化序列;Original battery health state change sequence acquisition unit: According to the characteristics of the battery, intercept the battery platform voltage range, use this as a reference to obtain a period of time from the historical operating data charging section voltage and current data of the battery charging platform interval, and obtain it through ampere-hour integration Get the original battery health state change sequence for each charging capacity value;
分离单元:使用经验模态分解算法从原始的电池健康状态变化序列中分离出电池周期变化段和线性衰退段;Separation unit: Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state change sequence;
电池变化趋势预测单元:分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Battery change trend prediction unit: establish data-driven models for the cycle change segment and the linear decline segment respectively, and realize the change trend prediction of the battery linear decline segment and the seasonal cycle segment.
本发明具有的有益效果在于:The beneficial effects that the present invention has are:
本发明通过截取平台区容量(能量)来近似电池全区间容量(能量),在保证一定精度的同时大大降低了电池健康状态评估对电池数据的需求,适用于工作在调峰工况下的储能系统SOH评估;解决了复杂环境下储能电池健康状态多模态叠加的问题,采用相应的信号分解算法实现电池SOH序列复杂因素的解耦合,从而关注电池衰退的本质规律,提高了电池健康状态预测和评估的精度。The present invention approximates the capacity (energy) of the entire battery interval by intercepting the capacity (energy) of the platform area, which greatly reduces the demand for battery data for battery health status evaluation while ensuring a certain accuracy, and is suitable for storage batteries working under peak-shaving conditions. Energy system SOH evaluation; solve the problem of multi-modal superposition of energy storage battery health status in a complex environment, and use the corresponding signal decomposition algorithm to realize the decoupling of complex factors in the battery SOH sequence, so as to pay attention to the essential law of battery decline and improve battery health Accuracy of state prediction and assessment.
附图说明Description of drawings
图1为本发明储能电站电池工况健康状态预测方法的流程示意图;Fig. 1 is a schematic flow chart of the method for predicting the health state of the battery working condition of the energy storage power station according to the present invention;
图2为本发明应用例中电池IC曲线图;Fig. 2 is the battery IC graph in the application example of the present invention;
图3-4为本发明应用例中工况数据处理前后对比图(图3为工况数据处理前的示意图,图4为工况数据处理后的示意图);Fig. 3-4 is the comparison chart before and after working condition data processing in the application example of the present invention (Fig. 3 is the schematic diagram before working condition data processing, and Fig. 4 is the schematic diagram after working condition data processing);
图5为本发明应用例中电池簇内所有单体的平台容量原始序列图;Fig. 5 is the original sequence diagram of the platform capacity of all monomers in the battery cluster in the application example of the present invention;
图6为本发明应用例中电池单体容量的预测值和实际值对比图;Fig. 6 is a comparison chart between the predicted value and the actual value of the battery cell capacity in the application example of the present invention;
图7为本发明储能电站电池工况健康状态预测系统的结构图。Fig. 7 is a structural diagram of the system for predicting the health state of the battery in the energy storage power station according to the present invention.
具体实施方式Detailed ways
下面结合说明书附图与具体实施方式,对本发明做进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings and specific implementation methods.
实施例1Example 1
参考图1,本实施例提供一种储能电站电池工况健康状态预测方法,其步骤如下:Referring to FIG. 1 , this embodiment provides a method for predicting the health state of a battery in an energy storage power station, the steps of which are as follows:
步骤1:获取储能电站电池历史运行数据,具体到电芯,采样精度1min~5min,对原始数据进行预处理,包括数据插补、数据清洗。Step 1: Obtain the historical operation data of the battery of the energy storage power station, specific to the battery cell, with a sampling accuracy of 1min to 5min, and preprocess the original data, including data interpolation and data cleaning.
步骤2:根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取原始的电池健康状态变化序列;Step 2: According to the characteristics of the battery, intercept the platform voltage range of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform range for a period of time from the charging period of historical operating data, and obtain the charging capacity value each time through the integration of ampere hours , to obtain the original battery health state change sequence;
步骤3:使用经验模态分解算法从原始的电池健康状态序列中分离出电池周期变化段和线性衰退段;Step 3: Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state sequence;
步骤4:分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Step 4: Establish the data-driven models for the periodical change segment and the linear decline segment respectively, and realize the forecast of the change trend of the battery linear decline segment and the seasonal cycle segment.
在上述技术方案的基础上,所述储能电站历史运行数据,及数据的预处 理方法:On the basis of the above technical solution, the historical operation data of the energy storage power station, and the data preprocessing method:
电池组电压数据、电池组电流数据、电池单体电压数据、电池单体电流数据、电池单体温度数据、电池模块温度数据及系统运行的累计时间。其中电池单体的电压数据、单体电流数据及温度数据是必要数据。Battery pack voltage data, battery pack current data, battery cell voltage data, battery cell current data, battery cell temperature data, battery module temperature data and the accumulated time of system operation. Among them, the voltage data, current data and temperature data of the battery cells are necessary data.
数据预处理:由于采样精度较差或监测系统短时故障有可能会造成某些采样数据出现不同程度的离群或者缺失会造成数据失真,从而影响后续的数据提取。因此,对在该步骤需要删除离群数据并对缺失数据进行插补得到可以进一步处理的可用数据。同时,电站并非每一天都参与调峰,因此某些时间电站电池不工作或仅工作在很小的DOD区间,该部分数据不具有代表性也很难从中提取新信息,因此该部分数据会被清理并置空。Data preprocessing: due to poor sampling accuracy or short-term failure of the monitoring system, some sampled data may appear outliers or missing data to varying degrees, resulting in data distortion, thereby affecting subsequent data extraction. Therefore, it is necessary to delete outlier data and interpolate missing data in this step to obtain usable data that can be further processed. At the same time, the power station does not participate in peak shaving every day, so the battery of the power station does not work or only works in a small DOD interval at certain times. This part of the data is not representative and it is difficult to extract new information from it, so this part of the data will be Clean up and make empty.
具体地,步骤2中所述的电池的平台电压区间通过电池容量增量曲线获取,要求电压平台要包括容量增量曲线的所有峰。Specifically, the platform voltage range of the battery described in step 2 is obtained through the battery capacity increment curve, and the voltage platform is required to include all peaks of the capacity increment curve.
具体地,电池的容量增量曲线由以下方法获取:Specifically, the capacity increment curve of the battery is obtained by the following method:
获取电池在某充电倍率下的完整充电段电压、电流数据。Obtain the voltage and current data of the complete charging section of the battery at a certain charging rate.
通过安时积分法获取电池在充电过程中各个采样点的累计充电安时数,得到电池电压-容量曲线。The accumulative charging ampere-hours of each sampling point in the charging process of the battery is obtained by the ampere-hour integration method, and the battery voltage-capacity curve is obtained.
设置合理的电压间隔(dV),通过线性插值方法获取按固定电压间隔增加的电压-容量曲线,随后通过差分方法获取电池IC曲线数据(IC(V)),如下式:Set a reasonable voltage interval (dV), obtain the voltage-capacity curve increasing at a fixed voltage interval by the linear interpolation method, and then obtain the battery IC curve data (IC(V)) by the differential method, as follows:
Figure PCTCN2022121940-appb-000002
Figure PCTCN2022121940-appb-000002
上式中,dV为电压间隔,dQ为电池等电压间隔对应的容量,IC为电池容量增量序列。In the above formula, dV is the voltage interval, dQ is the capacity corresponding to the equal voltage interval of the battery, and IC is the battery capacity increment sequence.
基于得到的具体IC曲线可以定量固定电池平台区的电压区间,得到获取平台区容量的基准值。基于得到的平台区电压基准,截取每一次电池充电段数据,通过安时积分法获取电池的充电段平台区容量序列即为电池SOH的参考序列。Based on the obtained specific IC curve, the voltage interval of the plateau area of the battery can be fixed quantitatively, and the reference value for obtaining the capacity of the plateau area can be obtained. Based on the obtained voltage reference of the platform area, the data of each battery charging section is intercepted, and the capacity sequence of the platform area of the charging section of the battery is obtained by the ampere-hour integration method, which is the reference sequence of the battery SOH.
具体地,为了分离获取序列中的季节周期性信号和电池的近似线性衰退信号,采用经验模态分解(EMD)算法对原始序列进行分解,其中EMD算法的步骤如下:Specifically, in order to separate the seasonal periodic signal and the approximate linear decay signal of the battery in the acquisition sequence, the empirical mode decomposition (EMD) algorithm is used to decompose the original sequence, and the steps of the EMD algorithm are as follows:
1.根据从原始信号中找出的全部局部极大值和极小值点,借助曲线拟合方 法,构造出原始信号的上包络线和下包络线;1. Construct the upper envelope and lower envelope of the original signal according to all the local maximum and minimum points found from the original signal, with the help of curve fitting method;
2.求上、下包络线的均值并构造均值曲线,然后用原始信号减去该曲线,得到信号的中间分量;2. Find the mean value of the upper and lower envelopes and construct the mean value curve, and then subtract the curve from the original signal to obtain the intermediate component of the signal;
3.判断所得的信号中间分量是否满足本征模函数(IMF)的约束条件:如果满足,则该分量就是一个IMF;如果不满足,则以该分量为待分解信号重复步骤一至步骤三,直到满足IMF的约束条件,或筛分门限值小于预设门限值(通常取0.2~0.3);3. Judge whether the intermediate component of the obtained signal satisfies the constraints of the intrinsic mode function (IMF): if it is satisfied, then the component is an IMF; if not, repeat steps 1 to 3 with the component as the signal to be decomposed until Meet the constraints of the IMF, or the screening threshold is less than the preset threshold (usually 0.2 to 0.3);
4.通过上述步骤得到第一个IMF后,用原始信号与该IMF相减,作为新的原始信号重复步骤1至步骤4,迭代直到残差分量满足预设条件,剩余信号分量即为残差分量,迭代结束。4. After the first IMF is obtained through the above steps, subtract the original signal from the IMF, repeat steps 1 to 4 as a new original signal, iterate until the residual component meets the preset conditions, and the remaining signal component is the residual component, the iteration ends.
经过以上分解过程,电池衰退原始序列被分解为若干经验模态分量及残差,其中残差序列多是信号的噪声,可以忽略。若干经验模态分量中包括不同时间尺度的局部特征信号,主要包括一些类周期信号及近似线性衰退信号。After the above decomposition process, the original sequence of battery decline is decomposed into several empirical mode components and residuals, among which the residual sequence is mostly the noise of the signal and can be ignored. Several empirical mode components include local characteristic signals of different time scales, mainly including some quasi-periodic signals and approximate linear decay signals.
具体地,基于获取的周期信号和线性信号分别建模,对近似线性序列建模的方式有很多如线性回归、支持向量机、高斯过程回归等回归算法,对周期序列,采用平方正弦核函数的高斯过程回归可以很好地刻画周期信号的周期变化特性。在本实施例中,为了保持统一,均采用高斯过程回归对两种序列进行建模:Specifically, based on the obtained periodic signal and linear signal, there are many ways to model the approximate linear sequence, such as linear regression, support vector machine, Gaussian process regression and other regression algorithms. For the periodic sequence, the square sine kernel function is used Gaussian process regression can well describe the periodic variation characteristics of periodic signals. In this example, in order to maintain unity, Gaussian process regression is used to model the two sequences:
采用高斯过程回归训练模型的步骤如下:The steps to train a model using Gaussian process regression are as follows:
1.根据EMD算法分解结果,分别选取线性衰退过程和周期过程数据作为训练集并开启交叉验证;1. According to the decomposition results of the EMD algorithm, select the linear decay process and periodic process data as the training set and enable cross-validation;
2.根据不同的训练集数据特征,选取不同的核函数,并初始化超参数;2. According to the data characteristics of different training sets, select different kernel functions and initialize hyperparameters;
3.使用训练集数据训练回归模型,采用极大似然估计获取超参数最优估计得到后验模型;3. Use the training set data to train the regression model, and use the maximum likelihood estimation to obtain the optimal estimate of the hyperparameters to obtain the posterior model;
4.利用训练后的模型对输入的测试集数据进行预测,并将预测值和验证集数据进行对比分析模型预测能力。4. Use the trained model to predict the input test set data, and compare the predicted value with the verification set data to analyze the predictive ability of the model.
以上,完成夹杂季节周期信号的电池健康状态趋势预测。Above, the battery health status trend prediction with seasonal cycle signals is completed.
实施例2Example 2
参考图7,本实施例提供一种储能电站电池工况健康状态预测系统,其包括:Referring to Fig. 7, this embodiment provides a system for predicting the state of health of a battery in an energy storage power station, which includes:
数据获取及预处理单元:获取储能电站电池历史运行数据,对原始数据进行预处理;Data acquisition and preprocessing unit: acquire the historical operation data of the battery of the energy storage power station, and preprocess the original data;
原始电池健康状态变化序列获取单元:根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取原始的电池健康状态变化序列;Original battery health state change sequence acquisition unit: According to the characteristics of the battery, intercept the battery platform voltage range, use this as a reference to obtain a period of time from the historical operating data charging section voltage and current data of the battery charging platform interval, and obtain it through ampere-hour integration Get the original battery health state change sequence for each charging capacity value;
分离单元:使用经验模态分解算法从原始的电池健康状态变化序列中分离出电池周期变化段和线性衰退段;Separation unit: Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state change sequence;
电池变化趋势预测单元:分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Battery change trend prediction unit: establish data-driven models for the cycle change segment and the linear decline segment respectively, and realize the change trend prediction of the battery linear decline segment and the seasonal cycle segment.
具体地,数据获取及预处理单元中,所述储能电站电池历史运行数据包括:电池组电压数据、电池组电流数据、电池单体电压数据、电池单体电流数据、电池单体温度数据、电池模块温度数据及储能电站运行的累计时间;Specifically, in the data acquisition and preprocessing unit, the historical battery operation data of the energy storage power station includes: battery pack voltage data, battery pack current data, battery cell voltage data, battery cell current data, battery cell temperature data, Battery module temperature data and the accumulated time of energy storage power station operation;
所述的数据预处理:采用数据插值,将数据的采样间隔统一;对已经统一采样间隔的数据进行去除离群点、插补和平滑处理。The data preprocessing: using data interpolation to unify the sampling interval of the data; removing outliers, interpolating and smoothing the data whose sampling interval has been unified.
具体地,原始电池健康状态变化序列获取单元中,电池的平台电压区间通过电池容量增量曲线进行界定,要求电压平台包括电池容量增量曲线的所有峰。Specifically, in the acquisition unit of the original battery health state change sequence, the plateau voltage range of the battery is defined by the battery capacity increment curve, and the voltage plateau is required to include all peaks of the battery capacity increment curve.
所述电池容量增量曲线由以下方法获取:The battery capacity increment curve is obtained by the following method:
获取电池在某充电倍率下的完整充电段电压、电流数据;Obtain the voltage and current data of the complete charging section of the battery at a certain charging rate;
通过安时积分法获取电池在充电过程中各个采样点的累计充电安时数得到电池电压-容量曲线;Obtain the cumulative charging ampere-hours of each sampling point during the charging process by the ampere-hour integration method to obtain the battery voltage-capacity curve;
设置合理的电压间隔dV,通过线性插值方法获取按固定电压间隔增加的电压-容量曲线,随后通过差分方法获取电池IC曲线数据,如下式:Set a reasonable voltage interval dV, obtain the voltage-capacity curve increasing at a fixed voltage interval through the linear interpolation method, and then obtain the battery IC curve data through the differential method, as follows:
Figure PCTCN2022121940-appb-000003
Figure PCTCN2022121940-appb-000003
上式中,dV为电压间隔,dQ为电池等电压间隔对应的容量,IC为电池容量增量。In the above formula, dV is the voltage interval, dQ is the capacity corresponding to the equal voltage interval of the battery, and IC is the battery capacity increment.
基于电池IC曲线定量固定电池平台区的电压区间,得到获取平台区容量的基准值;基于平台区容量的基准值,截取每一次电池充电段数据,通过安时积分法获取电池的充电段平台区容量序列,即为原始的电池健康状态变化序列。Quantitatively fix the voltage interval of the battery platform area based on the battery IC curve, and obtain the benchmark value of the platform area capacity; based on the benchmark value of the platform area capacity, intercept the data of each battery charging section, and obtain the platform area of the battery charging section through the ampere-hour integral method The capacity sequence is the original battery health state change sequence.
具体地,所述分离单元的具体步骤如下:Specifically, the specific steps of the separation unit are as follows:
a.根据从原始信号中找出的全部局部极大值和极小值点,借助曲线拟合方法,构造出原始信号的上包络线和下包络线;所述的原始信号取自原始的电池健康状态变化序列;a. According to all the local maximum and minimum points found from the original signal, the upper envelope and the lower envelope of the original signal are constructed by means of the curve fitting method; the original signal is taken from the original The battery health state change sequence;
b.求上、下包络线的均值并构造均值曲线,然后用原始信号减去该曲线,得到信号的中间分量;b. Find the mean value of the upper and lower envelopes and construct the mean value curve, and then subtract the curve from the original signal to obtain the intermediate component of the signal;
c.判断所得的信号中间分量是否满足本征模函数IMF的约束条件:如果满足,则该分量就是一个IMF;如果不满足,则以该分量为待分解信号重复步骤a至步骤b,直到满足IMF的约束条件,或筛分门限值小于预设门限值;c. Determine whether the intermediate component of the obtained signal satisfies the constraints of the intrinsic mode function IMF: if it is satisfied, the component is an IMF; if it is not satisfied, repeat steps a to b with the component as the signal to be decomposed until it is satisfied The constraint condition of IMF, or the screening threshold value is less than the preset threshold value;
d.通过上述步骤得到第一个IMF后,用原始信号与该第一个IMF相减,作为新的原始信号重复步骤a至步骤c,迭代直到残差分量满足预设条件,剩余信号分量即为残差分量,迭代结束;d. After the first IMF is obtained through the above steps, subtract the original signal from the first IMF, repeat steps a to c as a new original signal, iterate until the residual component meets the preset condition, and the remaining signal component is is the residual component, and the iteration ends;
经过以上分解过程,原始的电池健康状态变化序列被分解为若干经验模态分量及残差序列,其中残差序列多是信号的噪声,忽略;若干经验模态分量序列包括一些类周期信号序列及近似线性衰退信号序列。After the above decomposition process, the original battery health state change sequence is decomposed into a number of empirical modal components and residual sequences, in which the residual sequence is mostly the noise of the signal, which is ignored; some empirical modal component sequences include some periodic signal sequences and Approximate linearly decaying signal sequence.
具体地,电池变化趋势预测单元的具体内容如下:Specifically, the specific content of the battery change trend prediction unit is as follows:
基于获取的类周期信号序列和近似线性衰退信号序列分别建模,对近似线性衰退信号序列建模的方式包括线性回归、支持向量机和高斯过程回归,对类周期信号序列,采用平方正弦核函数的高斯过程回归刻画类周期信号的周期变化特性;Based on the obtained quasi-periodic signal sequence and approximate linear decay signal sequence, the modeling methods include linear regression, support vector machine and Gaussian process regression, and the square sine kernel function is used for the quasi-periodic signal sequence The Gaussian process regression characterizes the periodic variation characteristics of periodic signals;
采用高斯过程回归建模的步骤如下:The steps for modeling with Gaussian process regression are as follows:
a.根据经验模态分解算法分解结果,分别选取线性衰退过程和周期过程数据作为训练集并开启交叉验证;a. According to the decomposition results of the empirical mode decomposition algorithm, the data of the linear decay process and the periodic process are respectively selected as the training set and cross-validation is turned on;
b.根据不同的训练集数据特征,选取不同的核函数,并初始化超参数;b. According to different training set data characteristics, select different kernel functions and initialize hyperparameters;
c.使用训练集数据训练回归模型,采用极大似然估计获取超参数最优估计 得到后验模型;c. Use the training set data to train the regression model, and use the maximum likelihood estimation to obtain the optimal estimate of the hyperparameters to obtain the posterior model;
d.利用训练后的模型对输入的测试集数据进行预测,并将预测值和验证集数据进行对比分析模型预测能力。d. Use the trained model to predict the input test set data, and compare the predicted value with the verification set data to analyze the predictive ability of the model.
实施例2中未详细说明的部分参见实施例1。Refer to Example 1 for the parts not described in detail in Example 2.
应用例Application example
应用例数据来自于某梯次利用电池储能示范工程工况运行数据,数据记录周期为1年,具体的实时步骤如下:The data of the application example comes from the operation data of a demonstration project using battery energy storage in an echelon. The data recording period is 1 year. The specific real-time steps are as follows:
步骤1,获取储能电站历史运行数据,具体到电芯,采样精度1min~5min,对原始数据进行预处理,包括数据插补、数据清洗。 Step 1. Obtain the historical operation data of the energy storage power station, specific to the battery cell, with a sampling accuracy of 1min to 5min, and perform preprocessing on the original data, including data interpolation and data cleaning.
由于不同时间的数据采样间隔不统一,为了便于分析计算首先对数据采样点统一,采取的措施是对采样较疏的数据进行插值,将数据的采样间隔统一为1min。由于原始电压数据因采样精度不足数据质量较差,对已经统一采样间隔的数据进行去除离群点、插补、平滑处理得到电池较平滑的电压曲线,处理前后的电池电压曲线如图2所示。Since the data sampling intervals at different times are not uniform, in order to facilitate analysis and calculation, the data sampling points are first unified, and the measure taken is to interpolate the sparsely sampled data and unify the data sampling interval to 1min. Due to the poor quality of the original voltage data due to insufficient sampling accuracy, outlier points, interpolation, and smoothing are performed on the data that has been uniformly sampled to obtain a smoother voltage curve of the battery. The battery voltage curve before and after processing is shown in Figure 2 .
步骤2,根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取电池健康状态变化序列。 Step 2, according to the characteristics of the battery, intercept the platform voltage range of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform range for a period of time from the charging period of historical operating data, and obtain the charging capacity value each time through the ampere-hour integration , to obtain the battery health state change sequence.
其中,电池平台区电压区间通过电池的容量增量曲线进行界定。在本实施例中,取和电站电池同一批次的某款磷酸铁锂电池,获取电池在一定倍率下恒流充电数据绘制电池容量增量曲线。在本实施例中,dV取固定值0.002V,采用差值方法可以获取对应电压间隔的dQ,根据容量增量曲线定义得到该电池的容量增量曲线如图3-4所示。Among them, the voltage interval of the battery platform area is defined by the capacity increment curve of the battery. In this embodiment, a lithium iron phosphate battery of the same batch as the power station battery is taken, and the constant current charging data of the battery at a certain rate is obtained to draw a battery capacity increment curve. In this embodiment, dV takes a fixed value of 0.002V, dQ corresponding to the voltage interval can be obtained by using the difference method, and the capacity increment curve of the battery is obtained according to the definition of the capacity increment curve, as shown in Figure 3-4.
Figure PCTCN2022121940-appb-000004
Figure PCTCN2022121940-appb-000004
通过容量增量曲线很容易确定平台区的范围,在本实施例中平台区电压范围限定在3.25~3.4V。基于得到的电压区间,截取每天的充电段数据并对电流作积分,获取电池平台区容量序列,在本实施例中为了使结果更直观未进行归一化的电池容量序列操作,得到某电池簇内所有单体的原始序列如图5 所示。The range of the plateau area can be easily determined through the capacity increment curve, and the voltage range of the plateau area is limited to 3.25-3.4V in this embodiment. Based on the obtained voltage interval, intercept the daily charging section data and integrate the current to obtain the battery platform area capacity sequence. In this embodiment, the normalized battery capacity sequence operation is not performed to make the result more intuitive, and a certain battery cluster is obtained. The original sequences of all the monomers in are shown in Figure 5.
通过原始序列可以发现,电池平台区容量变化趋势为多模态混叠,部分电池平台区容量受温度影响较明显,因此必须对原始数据进行信号分解。Through the original sequence, it can be found that the capacity change trend of the battery platform area is multi-modal aliasing, and the capacity of some battery platform areas is significantly affected by temperature, so the signal decomposition of the original data must be performed.
步骤3,在本应用例中,对电池历史数据采用经验模态分解方法进行分解,利用Matlab软件信号处理工具箱中的EMD算法工具包即可实现模态分解,在本实施例中,最大IMF个数设置为3,至此,完后信号分解解耦工作,获取电池线性段和周期段序列。 Step 3. In this application example, the empirical mode decomposition method is used to decompose the historical data of the battery, and the mode decomposition can be realized by using the EMD algorithm toolkit in the signal processing toolbox of Matlab software. In this embodiment, the maximum IMF The number is set to 3. At this point, the post-signal decomposition and decoupling work is completed, and the battery linear segment and cycle segment sequence are obtained.
步骤4,分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Step 4. Establish data-driven models for the periodical variation segment and the linear recession segment respectively, so as to realize the prediction of the change trend of the battery linear recession segment and the seasonal cycle segment.
对分解出的线性衰退过程,采用高斯过程回归对该过程建模时,选取的核函数为可以表征电池长时间衰退趋势的平方指数核函数,该核函数的表达式为:For the decomposed linear decay process, when using Gaussian process regression to model the process, the selected kernel function is a square exponential kernel function that can characterize the battery's long-term decline trend. The expression of the kernel function is:
Figure PCTCN2022121940-appb-000005
Figure PCTCN2022121940-appb-000005
其中,σ、l和α为超参数。由于核函数解析表达式中均包含数个超参数,选取适当的参数值对于建立准确的高斯过程回归模型至关重要。超参数估计方法一般采用极大似然估计法,极大似然估计是一种常用的模型参数估计方法:设超参数集θ={θ 1,...,θ n},通过贝叶斯定理可知超参数的后验分布为: Among them, σ, l and α are hyperparameters. Since the analytical expression of the kernel function contains several hyperparameters, selecting appropriate parameter values is crucial for establishing an accurate Gaussian process regression model. The hyperparameter estimation method generally adopts the maximum likelihood estimation method, which is a commonly used model parameter estimation method: set the hyperparameter set θ={θ 1 ,...,θ n }, through Bayesian Theorem shows that the posterior distribution of hyperparameters is:
Figure PCTCN2022121940-appb-000006
Figure PCTCN2022121940-appb-000006
其中,x表示训练数据的输入,y表示训练数据的输出。若并无关于超参数集θ={θ 1,...,θ n}的先验知识,那么θ的最大后验估计为p(y|x,θ)的极大似然估计,即: Among them, x represents the input of training data, and y represents the output of training data. If there is no prior knowledge about the hyperparameter set θ={θ 1 ,...,θ n }, then the maximum a posteriori estimate of θ is the maximum likelihood estimate of p(y|x,θ), namely:
Figure PCTCN2022121940-appb-000007
Figure PCTCN2022121940-appb-000007
由此,即可根据训练数据得到核函数超参数的后验估计值。In this way, the posterior estimated value of the hyperparameters of the kernel function can be obtained according to the training data.
对于季节周期信号,本应用例中采用平方正弦核函数来描述样本变量之间的关系,其形式如下式所示:For seasonal periodic signals, in this application example, the square sine kernel function is used to describe the relationship between sample variables, and its form is shown in the following formula:
Figure PCTCN2022121940-appb-000008
Figure PCTCN2022121940-appb-000008
在电池衰退模型构建完成后,为验证模型精度,选取电池簇中两只电池单体,以2020年6月-2021年1月的现场数据为训练数据,预测2021年2月到6月的电池衰退情况,得到的电池单体容量的预测值和实际值如图6所示,电池的线性衰退趋势被很好的拟合,预测最大误差小于5%。After the battery decay model is built, in order to verify the accuracy of the model, two battery cells in the battery cluster are selected, and the field data from June 2020 to January 2021 are used as training data to predict the battery life from February to June 2021. For the decline situation, the predicted value and actual value of the obtained battery cell capacity are shown in Figure 6. The linear decline trend of the battery is well fitted, and the maximum prediction error is less than 5%.
本发明中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。In the present invention, specific examples are used to illustrate the principle and implementation of the present invention, and the descriptions of the above embodiments are only used to help understand the method and core idea of the present invention. It should be pointed out that for those skilled in the art, without departing from the principle of the present invention, some improvements and modifications can be made to the present invention, and these improvements and modifications also fall within the protection scope of the claims of the present invention.

Claims (10)

  1. 储能电站电池工况健康状态预测方法,其特征在于,包括:The method for predicting the state of health of a battery in an energy storage power station is characterized in that it includes:
    步骤1:获取储能电站电池历史运行数据,对原始数据进行预处理;Step 1: Obtain the historical operation data of the battery of the energy storage power station, and preprocess the original data;
    步骤2:根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取原始的电池健康状态变化序列;Step 2: According to the characteristics of the battery, intercept the platform voltage interval of the battery, and use this as a reference to obtain the voltage and current data of the battery charging platform interval for a period of time from the historical operation data charging section, and obtain the charging capacity value each time through the integration of ampere hours , to obtain the original battery health state change sequence;
    步骤3:使用经验模态分解算法从原始的电池健康状态变化序列中分离出电池周期变化段和线性衰退段;Step 3: Use the empirical mode decomposition algorithm to separate the battery cycle change segment and the linear decay segment from the original battery health state change sequence;
    步骤4:分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Step 4: Establish the data-driven models for the periodical change segment and the linear decline segment respectively, and realize the forecast of the change trend of the battery linear decline segment and the seasonal cycle segment.
  2. 如权利要求1所述的储能电站电池工况健康状态预测方法,其特征在于,步骤1中,所述储能电站电池历史运行数据包括:电池组电压数据、电池组电流数据、电池单体电压数据、电池单体电流数据、电池单体温度数据、电池模块温度数据及储能电站运行的累计时间;The method for predicting the state of health of a battery in an energy storage power station according to claim 1, wherein in step 1, the historical battery operation data of the energy storage power station includes: battery pack voltage data, battery pack current data, battery cell Voltage data, battery cell current data, battery cell temperature data, battery module temperature data and the accumulated time of energy storage power station operation;
    所述的数据预处理:采用数据插值,将数据的采样间隔统一;对已经统一采样间隔的数据进行去除离群点、插补和平滑处理。The data preprocessing: using data interpolation to unify the sampling interval of the data; removing outliers, interpolating and smoothing the data whose sampling interval has been unified.
  3. 如权利要求1所述的储能电站电池工况健康状态预测方法,其特征在于,步骤2中,电池的平台电压区间通过电池容量增量曲线进行界定,要求电压平台包括电池容量增量曲线的所有峰。The method for predicting the state of health of a battery in an energy storage power station according to claim 1, wherein in step 2, the platform voltage range of the battery is defined by the battery capacity incremental curve, and the voltage platform is required to include the battery capacity incremental curve all peaks.
  4. 如权利要求3所述的储能电站电池工况健康状态预测方法,其特征在于,步骤2中,所述电池容量增量曲线由以下方法获取:The method for predicting the state of health of a battery in an energy storage power station according to claim 3, wherein in step 2, the battery capacity increment curve is obtained by the following method:
    获取电池在某充电倍率下的完整充电段电压、电流数据;Obtain the voltage and current data of the complete charging section of the battery at a certain charging rate;
    通过安时积分法获取电池在充电过程中各个采样点的累计充电安时数得到电池电压-容量曲线;Obtain the cumulative charging ampere-hours of each sampling point during the charging process by the ampere-hour integration method to obtain the battery voltage-capacity curve;
    设置合理的电压间隔dV,通过线性插值方法获取按固定电压间隔增加的电压-容量曲线,随后通过差分方法获取电池IC曲线数据,如下式:Set a reasonable voltage interval dV, obtain the voltage-capacity curve increasing at a fixed voltage interval through the linear interpolation method, and then obtain the battery IC curve data through the differential method, as follows:
    Figure PCTCN2022121940-appb-100001
    Figure PCTCN2022121940-appb-100001
    上式中,dV为电压间隔,dQ为电池等电压间隔对应的容量,IC为电池容量增量。In the above formula, dV is the voltage interval, dQ is the capacity corresponding to the equal voltage interval of the battery, and IC is the battery capacity increment.
  5. 如权利要求4所述的储能电站电池工况健康状态预测方法,其特征在于,步骤2中,基于电池IC曲线定量固定电池平台区的电压区间,得到获取平台区容量的基准值;基于平台区容量的基准值,截取每一次电池充电段数据,通过安时积分法获取电池的充电段平台区容量序列,即为原始的电池健康状态变化序列。The method for predicting the state of health of a battery in an energy storage power station according to claim 4, wherein in step 2, the voltage interval of the battery platform area is quantitatively fixed based on the battery IC curve, and the reference value for obtaining the capacity of the platform area is obtained; based on the platform Intercept the data of each charging section of the battery, and obtain the platform area capacity sequence of the charging section of the battery through the ampere-hour integration method, which is the original battery health state change sequence.
  6. 如权利要求1所述的储能电站电池工况健康状态预测方法,其特征在于,步骤3的具体步骤如下:The method for predicting the state of health of a battery in an energy storage power station according to claim 1, wherein the specific steps of step 3 are as follows:
    a.根据从原始信号中找出的全部局部极大值和极小值点,借助曲线拟合方法,构造出原始信号的上包络线和下包络线;所述的原始信号取自原始的电池健康状态变化序列;a. According to all the local maximum and minimum points found from the original signal, the upper envelope and the lower envelope of the original signal are constructed by means of the curve fitting method; the original signal is taken from the original The battery health state change sequence;
    b.求上、下包络线的均值并构造均值曲线,然后用原始信号减去该曲线,得到信号的中间分量;b. Find the mean value of the upper and lower envelopes and construct the mean value curve, and then subtract the curve from the original signal to obtain the intermediate component of the signal;
    c.判断所得的信号中间分量是否满足本征模函数IMF的约束条件:如果满足,则该分量就是一个IMF;如果不满足,则以该分量为待分解信号重复步骤a至步骤b,直到满足IMF的约束条件,或筛分门限值小于预设门限值;c. Judging whether the intermediate component of the obtained signal satisfies the constraints of the intrinsic mode function IMF: if it is satisfied, the component is an IMF; if it is not satisfied, then repeat steps a to b with this component as the signal to be decomposed until it is satisfied The constraint condition of IMF, or the screening threshold value is less than the preset threshold value;
    d.通过上述步骤得到第一个IMF后,用原始信号与该第一个IMF相减,作为新的原始信号重复步骤a至步骤c,迭代直到残差分量满足预设条件,剩余信号分量即为残差分量,迭代结束;d. After the first IMF is obtained through the above steps, subtract the original signal from the first IMF, repeat steps a to c as a new original signal, iterate until the residual component meets the preset condition, and the remaining signal component is is the residual component, and the iteration ends;
    经过以上分解过程,原始的电池健康状态变化序列被分解为若干经验模态分量及残差序列,其中残差序列多是信号的噪声,忽略;若干经验模态分量序列包括一些类周期信号序列及近似线性衰退信号序列。After the above decomposition process, the original battery health state change sequence is decomposed into a number of empirical modal components and residual sequences, in which the residual sequence is mostly the noise of the signal, which is ignored; some empirical modal component sequences include some periodic signal sequences and Approximate linearly decaying signal sequence.
  7. 如权利要求6所述的储能电站电池工况健康状态预测方法,其特征在于,步骤4的具体内容如下:The method for predicting the state of health of a battery in an energy storage power station according to claim 6, wherein the specific content of step 4 is as follows:
    基于获取的类周期信号序列和近似线性衰退信号序列分别建模,对近似线性衰退信号序列建模的方式包括线性回归、支持向量机和高斯过程回归,对类周期信号序列,采用平方正弦核函数的高斯过程回归刻画类周期信号的周期变化特性;Based on the obtained quasi-periodic signal sequence and approximate linear decay signal sequence, the modeling methods include linear regression, support vector machine and Gaussian process regression, and the square sine kernel function is used for the quasi-periodic signal sequence The Gaussian process regression characterizes the periodic variation characteristics of periodic signals;
    采用高斯过程回归建模的步骤如下:The steps for modeling with Gaussian process regression are as follows:
    a.根据经验模态分解算法分解结果,分别选取线性衰退过程和周期过程数 据作为训练集并开启交叉验证;a. According to the decomposition results of the empirical mode decomposition algorithm, the data of the linear decay process and the periodic process are selected as the training set and cross-validation is turned on;
    b.根据不同的训练集数据特征,选取不同的核函数,并初始化超参数;b. According to different training set data characteristics, select different kernel functions and initialize hyperparameters;
    c.使用训练集数据训练回归模型,采用极大似然估计获取超参数最优估计得到后验模型;c. Use the training set data to train the regression model, and use the maximum likelihood estimation to obtain the optimal estimate of the hyperparameters to obtain the posterior model;
    d.利用训练后的模型对输入的测试集数据进行预测,并将预测值和验证集数据进行对比分析模型预测能力。d. Use the trained model to predict the input test set data, and compare the predicted value with the verification set data to analyze the predictive ability of the model.
  8. 储能电站电池工况健康状态预测系统,其特征在于,包括:The system for predicting the state of health of a battery in an energy storage power station is characterized in that it includes:
    数据获取及预处理单元:获取储能电站电池历史运行数据,对原始数据进行预处理;Data acquisition and preprocessing unit: acquire the historical operation data of the battery of the energy storage power station, and preprocess the original data;
    原始电池健康状态变化序列获取单元:根据电池的特性,截取电池的平台电压区间,以此为参考从历史运行数据充电段中获取一段时间的电池充电平台区间电压电流数据,并通过安时积分获取每次的充电容量值,获取原始的电池健康状态变化序列;Original battery health state change sequence acquisition unit: According to the characteristics of the battery, intercept the battery platform voltage range, use this as a reference to obtain a period of time from the historical operating data charging section voltage and current data of the battery charging platform interval, and obtain it through ampere-hour integration Get the original battery health state change sequence for each charging capacity value;
    分离单元:使用经验模态分解算法从原始的电池健康状态变化序列中分离出电池周期变化段和线性衰退段;Separation unit: Use the empirical mode decomposition algorithm to separate the battery cycle change segment and linear decay segment from the original battery health state change sequence;
    电池变化趋势预测单元:分别建立周期变化段和线性衰退段的数据驱动模型,实现电池线性衰退段和季节性周期段的变化趋势预测。Battery change trend prediction unit: establish the data-driven models of the cycle change segment and the linear decline segment respectively, and realize the change trend prediction of the battery linear decline segment and the seasonal cycle segment.
  9. 根据权利要求8所述的储能电站电池工况健康状态预测系统,其特征在于,所述分离单元的具体步骤如下:The system for predicting the state of health of the battery working condition of the energy storage power station according to claim 8, wherein the specific steps of the separation unit are as follows:
    a.根据从原始信号中找出的全部局部极大值和极小值点,借助曲线拟合方法,构造出原始信号的上包络线和下包络线;所述的原始信号取自原始的电池健康状态变化序列;a. According to all the local maximum and minimum points found from the original signal, the upper envelope and the lower envelope of the original signal are constructed by means of the curve fitting method; the original signal is taken from the original The battery health state change sequence;
    b.求上、下包络线的均值并构造均值曲线,然后用原始信号减去该曲线,得到信号的中间分量;b. Find the mean value of the upper and lower envelopes and construct the mean value curve, and then subtract the curve from the original signal to obtain the intermediate component of the signal;
    c.判断所得的信号中间分量是否满足本征模函数IMF的约束条件:如果满足,则该分量就是一个IMF;如果不满足,则以该分量为待分解信号重复步骤a至步骤b,直到满足IMF的约束条件,或筛分门限值小于预设门限值;c. Determine whether the intermediate component of the obtained signal satisfies the constraints of the intrinsic mode function IMF: if it is satisfied, the component is an IMF; if it is not satisfied, repeat steps a to b with the component as the signal to be decomposed until it is satisfied The constraint condition of IMF, or the screening threshold value is less than the preset threshold value;
    d.通过上述步骤得到第一个IMF后,用原始信号与该第一个IMF相减,作为新的原始信号重复步骤a至步骤c,迭代直到残差分量满足预设条件,剩 余信号分量即为残差分量,迭代结束;d. After the first IMF is obtained through the above steps, subtract the original signal from the first IMF, repeat steps a to c as a new original signal, iterate until the residual component meets the preset condition, and the remaining signal component is is the residual component, and the iteration ends;
    经过以上分解过程,原始的电池健康状态变化序列被分解为若干经验模态分量及残差序列,其中残差序列多是信号的噪声,忽略;若干经验模态分量序列包括一些类周期信号序列及近似线性衰退信号序列。After the above decomposition process, the original battery health state change sequence is decomposed into a number of empirical modal components and residual sequences, in which the residual sequence is mostly the noise of the signal, which is ignored; some empirical modal component sequences include some periodic signal sequences and Approximate linearly decaying signal sequence.
  10. 如权利要求9所述的储能电站电池工况健康状态预测系统,其特征在于,电池变化趋势预测单元的具体内容如下:The system for predicting the state of health of a battery in an energy storage power station according to claim 9, wherein the specific content of the battery change trend prediction unit is as follows:
    基于获取的类周期信号序列和近似线性衰退信号序列分别建模,对近似线性衰退信号序列建模的方式包括线性回归、支持向量机和高斯过程回归,对类周期信号序列,采用平方正弦核函数的高斯过程回归刻画类周期信号的周期变化特性;Based on the obtained quasi-periodic signal sequence and approximate linear decay signal sequence, the modeling methods include linear regression, support vector machine and Gaussian process regression, and the square sine kernel function is used for the quasi-periodic signal sequence The Gaussian process regression characterizes the periodic variation characteristics of periodic signals;
    采用高斯过程回归建模的步骤如下:The steps for modeling with Gaussian process regression are as follows:
    a.根据经验模态分解算法分解结果,分别选取线性衰退过程和周期过程数据作为训练集并开启交叉验证;a. According to the decomposition results of the empirical mode decomposition algorithm, the data of the linear decay process and the periodic process are respectively selected as the training set and cross-validation is turned on;
    b.根据不同的训练集数据特征,选取不同的核函数,并初始化超参数;b. According to different training set data characteristics, select different kernel functions and initialize hyperparameters;
    c.使用训练集数据训练回归模型,采用极大似然估计获取超参数最优估计得到后验模型;c. Use the training set data to train the regression model, and use the maximum likelihood estimation to obtain the optimal estimate of the hyperparameters to obtain the posterior model;
    d.利用训练后的模型对输入的测试集数据进行预测,并将预测值和验证集数据进行对比分析模型预测能力。d. Use the trained model to predict the input test set data, and compare the predicted value with the verification set data to analyze the predictive ability of the model.
PCT/CN2022/121940 2022-01-07 2022-09-28 Method and system for predicting working condition health status of battery in energy storage power station WO2023130776A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210016460.8 2022-01-07
CN202210016460.8A CN114547849A (en) 2022-01-07 2022-01-07 Method and system for predicting working condition and health state of battery of energy storage power station

Publications (1)

Publication Number Publication Date
WO2023130776A1 true WO2023130776A1 (en) 2023-07-13

Family

ID=81670376

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/121940 WO2023130776A1 (en) 2022-01-07 2022-09-28 Method and system for predicting working condition health status of battery in energy storage power station

Country Status (2)

Country Link
CN (1) CN114547849A (en)
WO (1) WO2023130776A1 (en)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116699445A (en) * 2023-08-07 2023-09-05 江苏天合储能有限公司 Capacity prediction method and system for battery energy storage system
CN116736142A (en) * 2023-08-14 2023-09-12 新誉集团有限公司 Method, system and device for early warning of health condition of battery pack
CN116743180A (en) * 2023-08-08 2023-09-12 东莞市欧派奇电子科技有限公司 Intelligent storage method for energy storage power supply data
CN116736150A (en) * 2023-08-16 2023-09-12 杭州高特电子设备股份有限公司 Battery abnormality detection method, battery system and computer program
CN116774075A (en) * 2023-08-28 2023-09-19 清华四川能源互联网研究院 Lithium ion battery health state evaluation method and system
CN117007978A (en) * 2023-10-07 2023-11-07 中国华能集团清洁能源技术研究院有限公司 Battery voltage prediction method and device, electronic equipment and storage medium
CN117054892A (en) * 2023-10-11 2023-11-14 特变电工西安电气科技有限公司 Evaluation method, device and management method for battery state of energy storage power station
CN117151255A (en) * 2023-10-30 2023-12-01 速源芯(东莞)能源科技有限公司 Modularized energy storage battery charging and discharging scheduling management system
CN117235446A (en) * 2023-11-09 2023-12-15 北京豪能汇新能源科技有限公司 Automatic change power station data intelligent management system
CN117269766A (en) * 2023-08-28 2023-12-22 广东工业大学 Battery SOH prediction method for unbalanced use scene
CN117318255A (en) * 2023-11-30 2023-12-29 北京中铁建电气化设计研究院有限公司 Battery state analysis system and method based on big data visualization
CN117310506A (en) * 2023-11-23 2023-12-29 宁德时代新能源科技股份有限公司 Method, apparatus and computer readable storage medium for monitoring energy storage power station
CN117371726A (en) * 2023-10-20 2024-01-09 国网新疆电力有限公司 Energy storage information substation state sensing and performance evaluation system and method
CN117451113A (en) * 2023-12-22 2024-01-26 中国电建集团华东勘测设计研究院有限公司 Self-elevating platform spud leg structure health monitoring system based on optical fiber sensing
CN117494028A (en) * 2023-12-29 2024-02-02 山东国华科技发展有限公司 Energy storage data anomaly detection method and system for energy storage power station
CN117578664A (en) * 2023-11-28 2024-02-20 浙江智格科技有限公司 Intelligent Internet of things-based battery-changing cabinet safety control method and system
CN117686937A (en) * 2024-02-02 2024-03-12 河南科技学院 Method for estimating health state of single battery in battery system
CN117741450A (en) * 2024-02-21 2024-03-22 新风光电子科技股份有限公司 Energy storage battery detection method for electrical parameter analysis
CN117851179A (en) * 2024-01-15 2024-04-09 上海异工同智信息科技有限公司 Portable equipment system health state detection and evaluation method and system
CN117849644A (en) * 2024-01-10 2024-04-09 中国电力科学研究院有限公司 On-line charge health state assessment method for electric automobile power battery
CN117949831A (en) * 2024-03-27 2024-04-30 牡丹江师范学院 Adjustable physical similarity simulation experiment platform
CN118050644A (en) * 2024-02-04 2024-05-17 淮阴工学院 Hydrogen fuel cell performance degradation prediction method, device and medium
CN118068253A (en) * 2024-04-18 2024-05-24 山东德源电力科技股份有限公司 High-precision measurement transformer error measurement method
CN118155047A (en) * 2024-05-10 2024-06-07 中央民族大学 Small battery management method suitable for alpine region

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114547849A (en) * 2022-01-07 2022-05-27 国网浙江省电力有限公司电力科学研究院 Method and system for predicting working condition and health state of battery of energy storage power station
CN116819209B (en) * 2023-06-29 2023-12-26 广州市兆能有限公司 Method and system for testing network-related performance of energy storage power station
CN117289168B (en) * 2023-09-19 2024-05-28 北京翼新数智科技有限公司 Product quality monitoring system based on sensor data and intelligent algorithm
CN117148170B (en) * 2023-10-30 2024-01-09 深圳市普裕时代新能源科技有限公司 Battery energy storage system and energy storage test method thereof

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031153A (en) * 2018-10-16 2018-12-18 北京交通大学 A kind of health status On-line Estimation method of lithium ion battery
CN110221225A (en) * 2019-07-08 2019-09-10 中国人民解放军国防科技大学 Spacecraft lithium ion battery cycle life prediction method
CN110941929A (en) * 2019-12-06 2020-03-31 长沙理工大学 Battery health state assessment method based on ARMA and Elman neural network combined modeling
CN113075574A (en) * 2021-03-30 2021-07-06 上海交通大学 Battery health state prediction method and equipment based on self-adaptive information fusion
CN113627671A (en) * 2021-08-11 2021-11-09 万克能源科技有限公司 SOH prediction calculation method for single battery in energy storage scene
CN114547849A (en) * 2022-01-07 2022-05-27 国网浙江省电力有限公司电力科学研究院 Method and system for predicting working condition and health state of battery of energy storage power station

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031153A (en) * 2018-10-16 2018-12-18 北京交通大学 A kind of health status On-line Estimation method of lithium ion battery
CN110221225A (en) * 2019-07-08 2019-09-10 中国人民解放军国防科技大学 Spacecraft lithium ion battery cycle life prediction method
CN110941929A (en) * 2019-12-06 2020-03-31 长沙理工大学 Battery health state assessment method based on ARMA and Elman neural network combined modeling
CN113075574A (en) * 2021-03-30 2021-07-06 上海交通大学 Battery health state prediction method and equipment based on self-adaptive information fusion
CN113627671A (en) * 2021-08-11 2021-11-09 万克能源科技有限公司 SOH prediction calculation method for single battery in energy storage scene
CN114547849A (en) * 2022-01-07 2022-05-27 国网浙江省电力有限公司电力科学研究院 Method and system for predicting working condition and health state of battery of energy storage power station

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116699445A (en) * 2023-08-07 2023-09-05 江苏天合储能有限公司 Capacity prediction method and system for battery energy storage system
CN116699445B (en) * 2023-08-07 2023-10-20 江苏天合储能有限公司 Capacity prediction method and system for battery energy storage system
CN116743180A (en) * 2023-08-08 2023-09-12 东莞市欧派奇电子科技有限公司 Intelligent storage method for energy storage power supply data
CN116743180B (en) * 2023-08-08 2024-02-13 东莞市欧派奇电子科技有限公司 Intelligent storage method for energy storage power supply data
CN116736142B (en) * 2023-08-14 2023-10-24 新誉集团有限公司 Method, system and device for early warning of health condition of battery pack
CN116736142A (en) * 2023-08-14 2023-09-12 新誉集团有限公司 Method, system and device for early warning of health condition of battery pack
CN116736150B (en) * 2023-08-16 2023-11-03 杭州高特电子设备股份有限公司 Battery abnormality detection method, battery system and computer program
CN116736150A (en) * 2023-08-16 2023-09-12 杭州高特电子设备股份有限公司 Battery abnormality detection method, battery system and computer program
CN116774075A (en) * 2023-08-28 2023-09-19 清华四川能源互联网研究院 Lithium ion battery health state evaluation method and system
CN117269766A (en) * 2023-08-28 2023-12-22 广东工业大学 Battery SOH prediction method for unbalanced use scene
CN117269766B (en) * 2023-08-28 2024-05-14 广东工业大学 Battery SOH prediction method for unbalanced use scene
CN117007978A (en) * 2023-10-07 2023-11-07 中国华能集团清洁能源技术研究院有限公司 Battery voltage prediction method and device, electronic equipment and storage medium
CN117007978B (en) * 2023-10-07 2024-01-30 中国华能集团清洁能源技术研究院有限公司 Battery voltage prediction method and device, electronic equipment and storage medium
CN117054892B (en) * 2023-10-11 2024-02-27 特变电工西安电气科技有限公司 Evaluation method, device and management method for battery state of energy storage power station
CN117054892A (en) * 2023-10-11 2023-11-14 特变电工西安电气科技有限公司 Evaluation method, device and management method for battery state of energy storage power station
CN117371726A (en) * 2023-10-20 2024-01-09 国网新疆电力有限公司 Energy storage information substation state sensing and performance evaluation system and method
CN117371726B (en) * 2023-10-20 2024-05-24 国网新疆电力有限公司 Energy storage information substation state sensing and performance evaluation system and method
CN117151255B (en) * 2023-10-30 2024-01-26 速源芯(东莞)能源科技有限公司 Modularized energy storage battery charging and discharging scheduling management system
CN117151255A (en) * 2023-10-30 2023-12-01 速源芯(东莞)能源科技有限公司 Modularized energy storage battery charging and discharging scheduling management system
CN117235446B (en) * 2023-11-09 2024-01-26 北京豪能汇新能源科技有限公司 Automatic change power station data intelligent management system
CN117235446A (en) * 2023-11-09 2023-12-15 北京豪能汇新能源科技有限公司 Automatic change power station data intelligent management system
CN117310506A (en) * 2023-11-23 2023-12-29 宁德时代新能源科技股份有限公司 Method, apparatus and computer readable storage medium for monitoring energy storage power station
CN117310506B (en) * 2023-11-23 2024-04-12 宁德时代新能源科技股份有限公司 Method, apparatus and computer readable storage medium for monitoring energy storage power station
CN117578664A (en) * 2023-11-28 2024-02-20 浙江智格科技有限公司 Intelligent Internet of things-based battery-changing cabinet safety control method and system
CN117578664B (en) * 2023-11-28 2024-06-11 浙江智格科技有限公司 Intelligent Internet of things-based battery-changing cabinet safety control method and system
CN117318255B (en) * 2023-11-30 2024-02-20 北京中铁建电气化设计研究院有限公司 Battery state analysis system and method based on big data visualization
CN117318255A (en) * 2023-11-30 2023-12-29 北京中铁建电气化设计研究院有限公司 Battery state analysis system and method based on big data visualization
CN117451113A (en) * 2023-12-22 2024-01-26 中国电建集团华东勘测设计研究院有限公司 Self-elevating platform spud leg structure health monitoring system based on optical fiber sensing
CN117451113B (en) * 2023-12-22 2024-03-26 中国电建集团华东勘测设计研究院有限公司 Self-elevating platform spud leg structure health monitoring system based on optical fiber sensing
CN117494028A (en) * 2023-12-29 2024-02-02 山东国华科技发展有限公司 Energy storage data anomaly detection method and system for energy storage power station
CN117494028B (en) * 2023-12-29 2024-04-19 山东国华科技发展有限公司 Energy storage data anomaly detection method and system for energy storage power station
CN117849644A (en) * 2024-01-10 2024-04-09 中国电力科学研究院有限公司 On-line charge health state assessment method for electric automobile power battery
CN117851179B (en) * 2024-01-15 2024-06-04 上海异工同智信息科技有限公司 Portable equipment system health state detection and evaluation method and system
CN117851179A (en) * 2024-01-15 2024-04-09 上海异工同智信息科技有限公司 Portable equipment system health state detection and evaluation method and system
CN117686937B (en) * 2024-02-02 2024-04-12 河南科技学院 Method for estimating health state of single battery in battery system
CN117686937A (en) * 2024-02-02 2024-03-12 河南科技学院 Method for estimating health state of single battery in battery system
CN118050644A (en) * 2024-02-04 2024-05-17 淮阴工学院 Hydrogen fuel cell performance degradation prediction method, device and medium
CN117741450B (en) * 2024-02-21 2024-05-14 新风光电子科技股份有限公司 Energy storage battery detection method for electrical parameter analysis
CN117741450A (en) * 2024-02-21 2024-03-22 新风光电子科技股份有限公司 Energy storage battery detection method for electrical parameter analysis
CN117949831A (en) * 2024-03-27 2024-04-30 牡丹江师范学院 Adjustable physical similarity simulation experiment platform
CN117949831B (en) * 2024-03-27 2024-05-31 牡丹江师范学院 Adjustable physical similarity simulation experiment platform
CN118068253A (en) * 2024-04-18 2024-05-24 山东德源电力科技股份有限公司 High-precision measurement transformer error measurement method
CN118155047A (en) * 2024-05-10 2024-06-07 中央民族大学 Small battery management method suitable for alpine region

Also Published As

Publication number Publication date
CN114547849A (en) 2022-05-27

Similar Documents

Publication Publication Date Title
WO2023130776A1 (en) Method and system for predicting working condition health status of battery in energy storage power station
Deng et al. Data-driven state of charge estimation for lithium-ion battery packs based on Gaussian process regression
Liu et al. State-of-health estimation of lithium-ion batteries based on electrochemical impedance spectroscopy: a review
Wang et al. State of health estimation based on modified Gaussian process regression for lithium-ion batteries
CN110187290B (en) Lithium ion battery residual life prediction method based on fusion algorithm
Song et al. Remaining life prediction of lithium-ion batteries based on health management: A review
CN111398833A (en) Battery health state evaluation method and evaluation system
CN112415414A (en) Method for predicting remaining service life of lithium ion battery
CN113866655B (en) Power battery capacity assessment method integrating vehicle networking operation data and test data
Wang et al. Multiple indicators-based health diagnostics and prognostics for energy storage technologies using fuzzy comprehensive evaluation and improved multivariate grey model
CN115902647B (en) Intelligent battery state monitoring method
US20220170995A1 (en) Method and Apparatus for Predicting a State of Health of a Device Battery in a Battery-operated Device
CN111983474A (en) Lithium ion battery life prediction method and system based on capacity decline model
CN113406558A (en) Electric meter misalignment detection method and device based on linear regression and electronic equipment
Kharlamova et al. A digital twin of battery energy storage systems providing frequency regulation
Liu et al. Review on degradation mechanism and health state estimation methods of lithium-ion batteries
Lin et al. State of health estimation of lithium-ion batteries based on a novel indirect health indicator
CN115293257A (en) Detection method and system for abnormal electricity utilization user
CN117783745B (en) Data online monitoring method and system for battery replacement cabinet
Xiao et al. State of health estimation framework of li-on battery based on improved Gaussian process regression for real car data
Dong et al. State of health estimation and remaining useful life estimation for Li-ion batteries based on a hybrid kernel function relevance vector machine
Yang et al. Probabilistic revenue analysis of microgrid considering source-load and forecast uncertainties
Jiang et al. A battery capacity estimation method based on the equivalent circuit model and quantile regression using vehicle real-world operation data
CN116249186A (en) Data processing method and device of wireless network equipment, storage medium and electronic equipment
CN115640542A (en) Intelligent electric energy meter state evaluation method and evaluation device based on Bayesian theory

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22918230

Country of ref document: EP

Kind code of ref document: A1