CN113406522B - 一种电动车锂电池系统寿命预测评估方法及系统 - Google Patents
一种电动车锂电池系统寿命预测评估方法及系统 Download PDFInfo
- Publication number
- CN113406522B CN113406522B CN202110680465.6A CN202110680465A CN113406522B CN 113406522 B CN113406522 B CN 113406522B CN 202110680465 A CN202110680465 A CN 202110680465A CN 113406522 B CN113406522 B CN 113406522B
- Authority
- CN
- China
- Prior art keywords
- working condition
- sys
- life
- distribution
- vehicle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 title claims abstract description 30
- 229910052744 lithium Inorganic materials 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 28
- 238000005265 energy consumption Methods 0.000 claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims description 23
- 238000007599 discharging Methods 0.000 claims description 19
- 230000000284 resting effect Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004378 air conditioning Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 7
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000012804 iterative process Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 229910001416 lithium ion Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/392—Determining battery ageing or deterioration, e.g. state of health
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/382—Arrangements for monitoring battery or accumulator variables, e.g. SoC
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/396—Acquisition or processing of data for testing or for monitoring individual cells or groups of cells within a battery
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Secondary Cells (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Tests Of Electric Status Of Batteries (AREA)
Abstract
本发明公开了一种电动车锂电池系统寿命预测评估方法,根据系统参数和工况参数,结合设定的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布。电芯寿命参数和统计的工况分布输入至电芯寿命模型内,通过加权算法可得系统在此工况下的日历寿命曲线和循环寿命曲线。整车参数和统计的工况分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正。最终,系统寿命迭代模型整合所有参数和数据,对系统的寿命进行预测并得到置信区间本发明在整车寿命预测中结合单体电芯寿命参数,考虑整车运行工况对系统寿命的影响,对续航里程进行修正,实现了锂电池系统使用寿命的预测评估。
Description
技术领域
本发明涉及新能源汽车技术领域,具体来说是一种电动车锂电池系统寿命预测评估方法及系统。
背景技术
锂离子电池作为新一代的绿色能源,因其具有能量密度大、工作电压高、循环寿命长等优势,被广泛应用在电动汽车,储能系统领域。作为电动汽车的核心组成部分,锂电池系统的性能对整车使用有着至关重要的作用,需要对其有充分的研究来确保电池包的经济寿命和寿命期内的可靠性。电池系统最为大家关注的就是电池容量的退化。
其中,锂电池的衰减可分为日历衰减和循环衰减两个部分。日历衰减受环境温度和搁置SOC状态影响,循环衰减受环境温度、充放电倍率、放电深度等因素影响。目前,对锂电池系统进行寿命评估只是简单的将单体电芯的日历衰减和循环衰减相加得到系统的衰减情况。不考虑整车工况、日历衰减和循环衰减耦合作用的影响。因此,无法较为精准地评估锂电池系统的寿命趋势。本发明在整车寿命预测中结合单体电芯寿命数据,考虑整车运行工况对系统寿命的影响,对续航里程进行修正,采用迭代算法对日历衰减和循环衰减进行耦合计算,实现了锂电池系统使用寿命的预测评估。
如申请号为CN201910023616.3公开的一种锂离子电池的寿命预测方法,该方法,首先建立表征电池存储寿命的第一容量衰减模型,根据存储实验的存储实验数据求解第一容量衰减模型的未知参数,然后,根据确定的第一容量衰减模型,结合电池的在设定时间内的充放电循环次数,建立第二容量衰减模型,根据充放电循环实验的循环实验数据,求解第二容量衰减模型的未知参数,从而最终确定表征电池存储寿命和循环寿命的第二容量衰减模型。最后,将电池实际使用的工况信息代入第二容量衰减模型,从而求得电池的预测寿命。相对于现有技术中单纯预测存储寿命或循环寿命的方法更符合电池的实际使用情况,能够快速预测电池的实际使用寿命。该方法针对单体电芯进行寿命预测,没有考虑整车公开,无法直接得到电池系统的寿命预测结果。
发明内容
本发明所要解决的技术问题在于提供一种考虑整车运行工况对电池系统寿命影响的锂电池系统寿命预测评估方法。
本发明通过以下技术手段实现解决上述技术问题的:
一种电动车锂电池系统寿命预测评估方法,包括以下步骤:
S1、统计具体车辆的工况数据,至少包括每日里程、平均速度和运营区域作为寿命预测的工况参数,将总电量和续航里程作为寿命预测的整车参数;
S2、对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;
S3、根据整车参数和工况参数,结合对应的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布;其中,工况分布概率模型包括私家车工况模型、运营车工况模型;工况分布概率模型为:
将车辆起始出行时刻的分布主要划分为n个时间区间,μn表示第n个时间区间,f(x)表示当前区间某一时刻的概率,wn为当前时间区间的权重,σn为当前时间区间的概率参数。
S4、电芯寿命数据和各工况的特征分布输入至电芯寿命模型内,通过加权算法可得电池系统在此工况下的日历寿命曲线和循环寿命曲线;
S5、整车参数和统计的工况特征分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正;
S6、电池系统寿命迭代整合整车参数和工况参数,对电池系统的寿命进行预测。
本发明在整车寿命预测中结合单体电芯寿命数据,考虑整车运行工况对系统寿命的影响,对续航里程进行修正,采用迭代算法对日历衰减和循环衰减进行耦合计算,实现了锂电池系统使用寿命的预测评估。
进一步的,所述步骤S3中将车辆起始出行时刻的分布主要划分为四个时间区间,分别为[6,10)、[10,14)、[14,20)、[20,24]。
进一步的,所述私家车工况模型和运营车工况模型中的权重和概率参数分别设置为:
1)私家车工况分布概率模型对应于四个时间区间的概率参数σn分别为:1、1.5、1.5、0.5;四个时间区间对应的权重wn分别为:0.35、0.15、0.45、0.1。
2)运营工况分布概率模型对应于四个时间区间的概率参数σn分别为:2,2,2,2;四个时间区间对应的权重wn分别为:0.3、0.2、0.3、0.2。
进一步的,所述步骤S3中工况包括充电、放电、搁置;工况特征分布包括温度分布、SOC分布、电流倍率分布;所述工况特征分布具体为:
其中,cha、dch、sta分别表示充电状态、放电状态、搁置状态。SOC每隔20%划分,电流倍率每隔0.2C划分,温度每隔10℃划分;因此,wsta/sic_50表示搁置状态下,SOC位于[40%,60%]的时间占比,wdch/temp_25表示放电状态下,温度位于[20℃,30℃]的时间占比,wcha/rate_0.5表示放充电状态下,充电倍率位于[0.4C,0.6C]的时间占比。
进一步的,所述步骤S4中加权算法如下:
Calsys(τ)=∑,Calcell(T,s,τ)*wsta/soc_s*wsta/temp_t
Cycsys(n)=∑Cyccell(T,s,n)*wcha/rate_r*wcha/temp_t
其中,Calsys(τ)、Cycsys(τ)、Calcell(T,s,τ)、Cyccell(T,s,τ)分别表示电池系统的日历寿命曲线和循环寿命曲线、单体电芯的日历寿命曲线和循环寿命曲线;单体电芯的日历寿命曲线是温度(T)、SOC(s)和搁置时间(τ)的衰减曲线;单体电芯的循环寿命曲线是温度(T)、充电倍率(r)和循环次数(n)的衰减曲线。
进一步的,所述步骤S5中修正因子如下:
f=fe/(fa*fm)
fe=∑Et*wdch/temp_t
其中,f、fe、fa、fm分别表示总修正因子、能量效率因子、配件能耗因子和载重影响因子;能量效率因子fe与单体电芯在不同温度下的放电能量效率相关,Et表示特定温度下单体电芯的放电能量效率;wt表示空调开启下温度区间的占比,et表示该温度条件的空调小时能耗;u和ec分别表示平均车速和整车的百公里能耗;load表示载重与整车重量的比值。
进一步的,所述步骤S6中迭代算法为:
Qcal,sys=Calsys[>=Qloss,sys](Δτ)
Qcyc,sys=Cycsys[>=Qloss,sys](Δn)
Qloss,sys+=Qcal,sys+Qcyc,sys
m+=Δm
其中,m、Qloss,sys分别表示迭代过程中累积行驶里程、累积容量衰减率,Δm、Δτ、Δn分别表示一次迭代周期内里程增长量、搁置时间、循环次数;Qcal,sys、Qcyc,sys分别表示一次迭代周期内日历衰减和循环衰减。
与上述方法对应的,本发明还提供一种电动车锂电池系统寿命预测评估系统,包括:
参数设置模块,统计具体车辆的工况数据,至少包括每日里程、平均速度和运营区域作为寿命预测的工况参数,将总电量和续航里程作为寿命预测的整车参数;
电芯寿命参数预测模块,对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;
工况特征分布统计模块,根据整车参数和工况参数,结合对应的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布;其中,工况分布概率模型包括私家车工况模型、运营车工况模型;
电池系统寿命计算模块,电芯寿命数据和各工况的特征分布输入至电芯寿命模型内,通过加权算法可得电池系统在此工况下的日历寿命曲线和循环寿命曲线;
修正模块,整车参数和统计的工况特征分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正;
电池系统寿命预测模块,电池系统寿命迭代整合所有参数和数据,对电池系统的寿命进行预测。
进一步的,所述工况特征分布统计模块中工况分布概率模型为:
将车辆起始出行时刻的分布主要划分为四个时间区间,分别为[6,10)、[10,14)、[14,20)、[20,24],μn表示第n个时间区间,f(x)表示当前区间某一时刻的概率,wn为当前时间区间的权重,σn为当前时间区间的概率参数。
进一步的,所述私家车工况模型和运营车工况模型中的权重和概率参数分别设置为:
1)私家车工况分布概率模型对应于四个时间区间的概率参数σn分别为:1、1.5、1.5、0.5;四个时间区间对应的权重wn分别为:0.35、0.15、0.45、0.1。
2)运营工况分布概率模型对应于四个时间区间的概率参数σn分别为:2,2,2,2;四个时间区间对应的权重wn分别为:0.3、0.2、0.3、0.2。
进一步的,所述工况特征分布统计模块中工况包括充电、放电、搁置;工况特征分布包括温度分布、SOC分布、电流倍率分布;所述工况特征分布具体为:
其中,cha、dch、sta分别表示充电状态、放电状态、搁置状态。SOC每隔20%划分,电流倍率每隔0.2C划分,温度每隔10℃划分;因此,wsta/soc_50表示搁置状态下,SOC位于[40%,60%]的时间占比,wdch/temp_25表示放电状态下,温度位于[20℃,30℃]的时间占比,wcha/rate_0.5表示放充电状态下,充电倍率位于[0.4C,0.6C]的时间占比。
进一步的,所述步骤S4中加权算法如下:
Calsys(τ)=∑Calcell(T,s,τ)*wsta/soc_s*wsta/temp_t
Cycsys(n)=∑Cyccell(T,s,n)*wcha/rate_r*wcha/temp_t
其中,Calsys(τ)、Cycsys(τ)、Calcell(T,s,τ)、Cyccell(T,s,τ)分别表示电池系统的日历寿命曲线和循环寿命曲线、单体电芯的日历寿命曲线和循环寿命曲线;单体电芯的日历寿命曲线是温度(T)、SOC(s)和搁置时间(τ)的衰减曲线;单体电芯的循环寿命曲线是温度(T)、充电倍率(r)和循环次数(n)的衰减曲线。
进一步的,所述修正模块中修正因子如下:
f=fe/(fa*fm)
fe=ΣEt*wdch/temp_t
其中,f、fe、fa、fm分别表示总修正因子、能量效率因子、配件能耗因子和载重影响因子;能量效率因子fe与单体电芯在不同温度下的放电能量效率相关,Et表示特定温度下单体电芯的放电能量效率;wt表示空调开启下温度区间的占比,et表示该温度条件的空调小时能耗;u和ec分别表示平均车速和整车的百公里能耗;load表示载重与整车重量的比值。
进一步的,所述电池系统寿命预测模块中迭代算法为:
Qcal,sys=Calsys[>=Qloss,sys](Δτ)
Qcyc,sys=Cycsys[>=Qloss,sys](Δn)
Qloss,sys+=Qcal,sys+Qcyc,sys
m+=Δm
其中,m、Qloss,sys分别表示迭代过程中累积行驶里程、累积容量衰减率,Δm、Δτ、Δn分别表示一次迭代周期内里程增长量、搁置时间、循环次数;Qcal,sys、Qcyc,sys分别表示一次迭代周期内日历衰减和循环衰减。
本发明的优点在于:
本发明根据整车参数和工况参数,结合设定的工况分布概率模型,随机模拟整车的具体工况,并统计各工况(放电工况、充电工况和搁置工况)的特征(SOC、温度、电流倍率)分布。电芯寿命数据和统计的工况分布输入至电芯寿命模型内,通过加权算法可得系统在此工况下的日历寿命曲线和循环寿命曲线。整车参数和统计的工况分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正。最终,系统寿命迭代模型整合所有参数和数据,对系统的寿命进行预测。本发明在整车寿命预测中结合单体电芯寿命数据,考虑整车运行工况对系统寿命的影响,对续航里程进行修正,采用迭代算法对日历衰减和循环衰减进行耦合计算,实现了锂电池系统使用寿命的预测评估。
附图说明
图1为本发明实施例电动车锂电池系统寿命预测评估方法逻辑结构框图;
图2为本发明实施例电动车锂电池系统寿命预测评估方法中迭代流程图;
图3为本发明实施例中路试车拟合得到的日历寿命曲线和循环寿命曲线。
图4为本发明实施例中路试车的寿命预测结果。
图5为本发明实施例中路试车的寿命预测和实测结果对比。
图6为本发明实施例中市场运营车辆容量衰减和模型预测的对比。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例公开一种电动车锂电池系统寿命预测评估方法,如图1所示,该方法包括三个维度参数输入、四个层次评估模型和1个指标结果输出。三个维度参数输入包括:系统参数(总电量和续航里程)、工况参数(每日里程、平均速度和运营区域)、电芯寿命数据(日历寿命和循环寿命数据)。四个层次评估模型包括:工况模拟模型、电芯寿命模型、续航里程修正模型和系统寿命迭代模型。1个指标结果输出包括:寿命预测。具体方法包括以下步骤:
步骤1、统计具体车辆的工况数据,至少包括每日里程、平均速度和运营区域作为寿命预测的工况参数,将总电量和续航里程作为寿命预测的整车参数;
步骤2、对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;
步骤3、根据整车参数和工况参数,结合对应的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布;其中,工况分布概率模型包括私家车工况模型、运营车工况模型;
工况分布概率模型为:
将车辆起始出行时刻的分布主要划分为四个时间区间,分别为[6,10)、[10,14)、[14,20)、[20,24],μn表示第n个时间区间,f(x)表示当前区间某一时刻的概率,wn为当前时间区间的权重,σn为当前时间区间的概率参数。
所述私家车工况模型和运营车工况模型中的权重和概率参数分别设置为:
1)私家车工况分布概率模型对应于四个时间区间的概率参数σn分别为:1、1.5、1.5、0.5;四个时间区间对应的权重wn分别为:0.35、0.15、0.45、0.1。
2)运营工况分布概率模型对应于四个时间区间的概率参数σn分别为:2,2,2,2;四个时间区间对应的权重wn分别为:0.3、0.2、0.3、0.2。
工况包括充电、放电、搁置;工况特征分布包括温度分布、SOC分布、电流倍率分布;所述工况特征分布具体为:
其中,cha、dch、sta分别表示充电状态、放电状态、搁置状态。SOC每隔20%划分,电流倍率每隔0.2C划分,温度每隔10℃划分;因此,wsta/soc_50表示搁置状态下,SOC位于[40%,60%]的时间占比,wdch/temp_25表示放电状态下,温度位于[20℃,30℃]的时间占比,wcha/rate_0.5表示放充电状态下,充电倍率位于[0.4C,0.6C]的时间占比。
步骤4、电芯寿命数据和各工况的特征分布输入至电芯寿命模型内,通过加权算法可得电池系统在此工况下的日历寿命曲线和循环寿命曲线;
权算法如下:
Calsys(τ)=∑Calcell(T,s,τ)*wsta/soc_s*wsta/temp_t
Cycsys(n)=∑Cyccell(T,s,n)*wcha/rate_r*wcha/temp_t
其中,Calsys(τ)、Cycsys(τ)、Calcell(T,s,τ)、Cyccell(T,s,τ)分别表示电池系统的日历寿命曲线和循环寿命曲线、单体电芯的日历寿命曲线和循环寿命曲线;单体电芯的日历寿命曲线是温度(T)、SOC(s)和搁置时间(τ)的衰减曲线;单体电芯的循环寿命曲线是温度(T)、充电倍率(r)和循环次数(n)的衰减曲线。
步骤5、整车参数和统计的工况特征分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正;
修正因子如下:
f=fe/(fa*fm)
fe=∑Et*wdch/temp_t
其中,f、fe、fa、fm分别表示总修正因子、能量效率因子、配件能耗因子和载重影响因子;能量效率因子fe与单体电芯在不同温度下的放电能量效率相关,Et表示特定温度下单体电芯的放电能量效率;wt表示空调开启下温度区间的占比,et表示该温度条件的空调小时能耗;u和ec分别表示平均车速和整车的百公里能耗;load表示载重与整车重量的比值。
步骤6、电池系统寿命迭代整合整车参数和工况参数,对电池系统的寿命进行预测。
代算法为:
Qcal,sys=Calsys[>=Qloss,sys](Δτ)
Qcyc,sys=Cycsys[>=Qloss,sys](Δn)
Qloss,sys+=Qcal,sys+Qcyc,sys
m+=Δm
迭代过程如图2所示,其中,m、day、Qloss,sys分别表示迭代过程中累积行驶里程、累积行驶天数和累积容量衰减率。d表示迭代周期,Δm、Δτ、Δn分别表示一次迭代周期内里程增长量、搁置时间、循环次数。Qcal,sys、Qcyc,sys分别表示一次迭代周期内日历衰减和循环衰减。M、Day、Qloss分别表示设定的迭代截止条件:总里程、总时间、总衰减量。
本实施例根据整车参数和工况参数,结合设定的工况分布概率模型,随机模拟整车的具体工况,并统计各工况(放电工况、充电工况和搁置工况)的特征(SOC、温度、电流倍率)分布。电芯寿命数据和统计的工况分布输入至电芯寿命模型内,通过加权算法可得系统在此工况下的日历寿命曲线和循环寿命曲线。整车参数和统计的工况分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正。最终,系统寿命迭代模型整合所有参数和数据,对系统的寿命进行预测。本实施例在整车寿命预测中结合单体电芯寿命数据,考虑整车运行工况对系统寿命的影响,对续航里程进行修正,采用迭代算法对日历衰减和循环衰减进行耦合计算,实现了锂电池系统使用寿命的预测评估。
下面对某款搭载锂电池的电动车进行长里程路试,并在市场运营的车辆中随机抽取100辆进行分析。数据分析处理步骤如下:
(1)统计路试车辆或100辆市场运营车的工况数据,统计数据中每日里程、平均速度和运营区域作为寿命预测模型的工况参数。另外,总电量和续航里程作为寿命预测模型的系统参数;
(2)对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;
(3)根据整车参数和工况参数,结合特定的工况分布概率模型,随机模拟整车的具体工况,并统计各工况(放电工况、充电工况和搁置工况)的特征(SOC、温度、电流倍率)分布。其中,特定的工况分布概率模型包括私家车工况模型、运营车工况模型;
(4)电芯寿命数据和上述统计的工况分布输入至电芯寿命模型内,通过加权算法可得系统在此工况下的日历寿命曲线和循环寿命曲线。
如图2所示,为路试车拟合得到的日历寿命曲线和循环寿命曲线。图3中日历寿命曲线和循环寿命曲线分别表示日历寿命随搁置时间的衰减量和循环寿命随循环次数的衰减量。
(5)整车参数和统计的工况分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正。
(6)系统寿命迭代模型整合所有参数和数据,对系统的寿命进行预测并得到置信区间。如图4所示,为路试车的寿命预测结果。图4中的日历衰减曲线和循环衰减曲线分别表示日历寿命和循环寿命在整车工况下随时间的衰减量。
将路试车的寿命预测结果和实测数据进行对比,如图5所示表示路试车预测结果与实测结果基本相吻合。
(7)在市场运营的车辆中随机抽取100辆进行分析,将市场运营车辆容量衰减和模型预测的对比,如图6所示表示市场运行车辆实测容量衰减在预测曲线周围波动,误差由工况差异导致。
本实施例还提供一种电动车锂电池系统寿命预测评估系统,包括:
参数设置模块,统计具体车辆的工况数据,至少包括每日里程、平均速度和运营区域作为寿命预测的工况参数,将总电量和续航里程作为寿命预测的整车参数;具体执行上述步骤1。
电芯寿命参数预测模块,对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;具体执行上述步骤2。
工况特征分布统计模块,根据整车参数和工况参数,结合对应的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布;其中,工况分布概率模型包括私家车工况模型、运营车工况模型;工况分布概率模型为:
将车辆起始出行时刻的分布主要划分为四个时间区间,分别为[6,10)、[10,14)、[14,20)、[20,24],μn表示第n个时间区间,f(x)表示当前区间某一时刻的概率,wn为当前时间区间的权重,σn为当前时间区间的概率参数。
所述私家车工况模型和运营车工况模型中的权重和概率参数分别设置为:
1)私家车工况分布概率模型对应于四个时间区间的概率参数σn分别为:1、1.5、1.5、0.5;四个时间区间对应的权重wn分别为:0.35、0.15、0.45、0.1。
2)运营工况分布概率模型对应于四个时间区间的概率参数σn分别为:2,2,2,2;四个时间区间对应的权重wn分别为:0.3、0.2、0.3、0.2。
工况包括充电、放电、搁置;工况特征分布包括温度分布、SOC分布、电流倍率分布;所述工况特征分布具体为:
其中,cha、dch、sta分别表示充电状态、放电状态、搁置状态。SOC每隔20%划分,电流倍率每隔0.2C划分,温度每隔10℃划分;因此,wsta/soc_50表示搁置状态下,SOC位于[40%,60%]的时间占比,wdch/temp_25表示放电状态下,温度位于[20℃,30℃]的时间占比,wcha/rate_0.5表示放充电状态下,充电倍率位于[0.4C,0.6C]的时间占比。
电池系统寿命计算模块,电芯寿命数据和各工况的特征分布输入至电芯寿命模型内,通过加权算法可得电池系统在此工况下的日历寿命曲线和循环寿命曲线;权算法如下:
Calsys(τ)=∑Calcell(T,s,τ)*wsta/soc_s*wsta/temp_t
Cycsys(n)=∑Cyccell(T,s,n)*wcha/rate_r*wcha/temp_t
其中,Calsys(τ)、Cycsys(τ)、Calcell(T,s,τ)、Cyccell(T,s,τ)分别表示电池系统的日历寿命曲线和循环寿命曲线、单体电芯的日历寿命曲线和循环寿命曲线;单体电芯的日历寿命曲线是温度(T)、SOC(s)和搁置时间(τ)的衰减曲线;单体电芯的循环寿命曲线是温度(T)、充电倍率(r)和循环次数(n)的衰减曲线。
修正模块,整车参数和统计的工况特征分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正;修正因子如下:
f=fe/(fa*fm)
fe=∑Et*wdch/temp_t
其中,f、fe、fa、fm分别表示总修正因子、能量效率因子、配件能耗因子和载重影响因子;能量效率因子fe与单体电芯在不同温度下的放电能量效率相关,Et表示特定温度下单体电芯的放电能量效率;wt表示空调开启下温度区间的占比,et表示该温度条件的空调小时能耗;u和ec分别表示平均车速和整车的百公里能耗;load表示载重与整车重量的比值。
电池系统寿命预测模块,电池系统寿命迭代整合所有参数和数据,对电池系统的寿命进行预测。代算法为:
Qcal,sys=Calsys[>=Qloss,sys](Δτ)
Qcyc,sys=Cycsys[>=Qloss,sys](Δn)
Qloss,sys+=Qcal,sys+Qcyc,sys
m+=Δm
迭代过程如图2所示,其中,m、day、Qloss,sys分别表示迭代过程中累积行驶里程、累积行驶天数和累积容量衰减率。d表示迭代周期,本实施例中迭代周期为1~30天,Δm、Δτ、Δn分别表示一次迭代周期内里程增长量、搁置时间、循环次数。Qcal,sys、Qcyc,sys分别表示一次迭代周期内日历衰减和循环衰减。M、Day、Qloss分别表示设定的迭代截止条件:总里程、总时间、总衰减量。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种电动车锂电池系统寿命预测评估方法,其特征在于,包括以下步骤:
S1、统计具体车辆的工况数据,至少包括每日里程、平均速度和运营区域作为寿命预测的工况参数,将总电量和续航里程作为寿命预测的整车参数;
S2、对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;
S3、根据整车参数和工况参数,结合对应的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布;其中,工况分布概率模型包括私家车工况模型、运营车工况模型;当前某一时刻工况分布概率为:
将车辆起始出行时刻的分布主要划分为n个时间区间,μn表示第n个时间区间,f(x)表示当前区间某一时刻的概率,wn为当前时间区间的权重,σn为当前时间区间的概率参数;
S4、电芯寿命参数和各工况的特征分布输入至电芯寿命模型内,通过加权算法可得电池系统在此工况下的日历寿命曲线和循环寿命曲线;
S5、整车参数和统计的工况特征分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正;
所述加权算法如下:
Calsys(τ)=∑Calcell(T,s,τ)*wsta/soc_s*wsta/temp_t
Cycsys(n)=∑Cyccell(T,s,n)*wcha/rate_r*wcha/temp_t
其中,Calsys(τ)、Cycsys(τ)、Calcell(T,s,τ)、Cyccell(T,s,τ)分别表示电池系统的日历寿命曲线和循环寿命曲线、单体电芯的日历寿命曲线和循环寿命曲线;单体电芯的日历寿命曲线是温度(T)、SOC(s)和搁置时间(τ)的衰减曲线;单体电芯的循环寿命曲线是温度(T)、充电倍率(r)和循环次数(n)的衰减曲线;cha、dch、sta分别表示充电状态、放电状态、搁置状态;wsta/soc_s表示搁置状态下,SOC位于s的时间占比,wdch/temp_t表示放电状态下,温度位于t的时间占比,wcha/rate_r表示放充电状态下,充电倍率位于r的时间占比;
S6、对电池系统的整车参数和工况参数进行迭代运算,对电池系统的寿命进行预测;
迭代算法为:
Qcal,sys=Calsys[>=Qloss,sys](Δτ)
Qcyc,sys=Cycsys[>=Qloss,sys](Δn)
Qloss,sys+=Qcal,sys+Qcyc,sys
m+=Δm
其中,m、Qloss,sys分别表示迭代过程中累积行驶里程、累积容量衰减率,Δm、Δτ、Δn分别表示一次迭代周期内里程增长量、搁置时间、循环次数;Qcal,sys、Qcyc,sys分别表示一次迭代周期内日历衰减和循环衰减。
2.根据权利要求1所述的一种电动车锂电池系统寿命预测评估方法,其特征在于,所述步骤S3将车辆起始出行时刻的分布主要划分为四个时间区间,分别为[6,10)、[10,14)、[14,20)、[20,24]。
3.根据权利要求2所述的一种电动车锂电池系统寿命预测评估方法,其特征在于,所述私家车工况模型和运营车工况模型中的权重和概率参数分别设置为:
1)私家车工况分布概率模型对应于四个时间区间的概率参数σn分别为:1、1.5、1.5、0.5;四个时间区间对应的权重wn分别为:0.35、0.15、0.45、0.1;
2)运营工况分布概率模型对应于四个时间区间的概率参数σn分别为:2,2,2,2;四个时间区间对应的权重wn分别为:0.3、0.2、0.3、0.2。
4.根据权利要求3所述的一种电动车锂电池系统寿命预测评估方法,其特征在于,所述步骤S3中工况包括充电、放电、搁置;工况特征分布包括温度分布、SOC分布、电流倍率分布;所述工况特征分布具体为:
SOC每隔20%划分,电流倍率每隔0.2C划分,温度每隔10℃划分;因此,wsta/soc_s表示搁置状态下,SOC位于s±20%的时间占比,wdch/temp_t表示放电状态下,温度位于t±10℃的时间占比,wcha/rate_r表示放充电状态下,充电倍率位于r±0.2C的时间占比。
6.一种电动车锂电池系统寿命预测评估系统,其特征在于,包括:
参数设置模块,统计具体车辆的工况数据,至少包括每日里程、平均速度和运营区域作为寿命预测的工况参数,将总电量和续航里程作为寿命预测的整车参数;
电芯寿命参数预测模块,对单体电芯进行日历寿命和循环寿命测试,测试结果作为寿命预测模型的电芯寿命参数;
工况特征分布统计模块,根据整车参数和工况参数,结合对应的工况分布概率模型,随机模拟整车的具体工况,并统计各工况的特征分布;其中,工况分布概率模型包括私家车工况模型、运营车工况模型;当前某一时刻工况分布概率为:
将车辆起始出行时刻的分布主要划分为n个时间区间,μn表示第n个时间区间,f(x)表示当前区间某一时刻的概率,wn为当前时间区间的权重,σn为当前时间区间的概率参数;
电池系统寿命计算模块,电芯寿命参数和各工况的特征分布输入至电芯寿命模型内,通过加权算法可得电池系统在此工况下的日历寿命曲线和循环寿命曲线;
修正模块,整车参数和统计的工况特征分布输入至续航里程修正模型,在能量效率、配件能耗和载重影响三个方面对续航里程进行修正;
所述加权算法如下:
Calsys(τ)=∑Calcell(T,s,τ)*wsta/soc_s*wsta/temp_t
Cycsys(n)=∑Cyccell(T,s,n)*wcha/rate_r*wcha/temp_t
其中,Calsys(τ)、Cycsys(τ)、Calcell(T,s,τ)、Cyccell(T,s,τ)分别表示电池系统的日历寿命曲线和循环寿命曲线、单体电芯的日历寿命曲线和循环寿命曲线;单体电芯的日历寿命曲线是温度(T)、SOC(s)和搁置时间(τ)的衰减曲线;单体电芯的循环寿命曲线是温度(T)、充电倍率(r)和循环次数(n)的衰减曲线;cha、dch、sta分别表示充电状态、放电状态、搁置状态;wsta/soc_s表示搁置状态下,SOC位于s的时间占比,wdch/temp_t表示放电状态下,温度位于t的时间占比,wcha/rate_r表示放充电状态下,充电倍率位于r的时间占比;
电池系统寿命预测模块,对电池系统的整车参数和工况参数进行迭代运算,对电池系统的寿命进行预测;迭代算法为:
Qcal,sys=Calsys[>=Qloss,sys](Δτ)
Qcyc,sys=Cycsys[>=Qloss,sys](Δn)
Qloss,sys+=Qcal,sys+Qcyc,sys
m+=Δm
其中,m、Qloss,sys分别表示迭代过程中累积行驶里程、累积容量衰减率,Δm、Δτ、Δn分别表示一次迭代周期内里程增长量、搁置时间、循环次数;Qcal,sys、Qcyc,sys分别表示一次迭代周期内日历衰减和循环衰减。
7.根据权利要求6所述的一种电动车锂电池系统寿命预测评估系统,其特征在于,所述工况特征分布统计模块中将车辆起始出行时刻的分布主要划分为四个时间区间,分别为[6,10)、[10,14)、[14,20)、[20,24]。
8.根据权利要求7所述的一种电动车锂电池系统寿命预测评估系统,其特征在于,所述私家车工况模型和运营车工况模型中的权重和概率参数分别设置为:
1)私家车工况分布概率模型对应于四个时间区间的概率参数σn分别为:1、1.5、1.5、0.5;四个时间区间对应的权重wn分别为:0.35、0.15、0.45、0.1;
2)运营工况分布概率模型对应于四个时间区间的概率参数σn分别为:2,2,2,2;四个时间区间对应的权重wn分别为:0.3、0.2、0.3、0.2。
9.根据权利要求8所述的一种电动车锂电池系统寿命预测评估系统,其特征在于,所述工况特征分布统计模块中工况包括充电、放电、搁置;工况特征分布包括温度分布、SOC分布、电流倍率分布;所述工况特征分布具体为:
SOC每隔20%划分,电流倍率每隔0.2C划分,温度每隔10℃划分;因此,wsta/soc_s表示搁置状态下,SOC位于s±20%的时间占比,wdch/temp_t表示放电状态下,温度位于t±10℃的时间占比,wchs/rate_r表示放充电状态下,充电倍率位于r±0.2C的时间占比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110680465.6A CN113406522B (zh) | 2021-06-18 | 2021-06-18 | 一种电动车锂电池系统寿命预测评估方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110680465.6A CN113406522B (zh) | 2021-06-18 | 2021-06-18 | 一种电动车锂电池系统寿命预测评估方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113406522A CN113406522A (zh) | 2021-09-17 |
CN113406522B true CN113406522B (zh) | 2022-05-03 |
Family
ID=77681710
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110680465.6A Active CN113406522B (zh) | 2021-06-18 | 2021-06-18 | 一种电动车锂电池系统寿命预测评估方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113406522B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114047452B (zh) * | 2022-01-13 | 2022-05-13 | 浙江玥视科技有限公司 | 一种确定电池循环寿命的方法及装置 |
CN116050662A (zh) * | 2023-03-07 | 2023-05-02 | 中环洁集团股份有限公司 | 一种环卫设备报废预测方法、系统及电子设备 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012127938A (ja) * | 2010-11-24 | 2012-07-05 | Mitsubishi Electric Corp | 蓄電装置の劣化監視方法、及びその劣化監視装置 |
CN105904981A (zh) * | 2016-04-07 | 2016-08-31 | 北京现代汽车有限公司 | 一种电动汽车续航里程估计控制方法、装置及整车控制器 |
JP2017502820A (ja) * | 2014-01-16 | 2017-01-26 | ボストン サイエンティフィック ニューロモデュレイション コーポレイション | 充電式バッテリを有する埋め込み可能医療デバイスに対する寿命末期の決定及び予想 |
CN107064800A (zh) * | 2016-11-29 | 2017-08-18 | 北京交通大学 | 锂离子电池剩余使用寿命的实时预测方法 |
CN107317057A (zh) * | 2016-04-27 | 2017-11-03 | 宝沃汽车(中国)有限公司 | 一种动力电池运行寿命预测和预测装置 |
CN108710084A (zh) * | 2018-06-01 | 2018-10-26 | 王泽霖 | 一种电动汽车电源监控及能源管理系统 |
CN109747428A (zh) * | 2019-02-27 | 2019-05-14 | 合肥国轩高科动力能源有限公司 | 一种电动汽车电池包剩余容量估算方法 |
CN110470991A (zh) * | 2019-08-23 | 2019-11-19 | 江西优特汽车技术有限公司 | 一种动力电池循环寿命评测方法 |
CN110901470A (zh) * | 2019-11-29 | 2020-03-24 | 安徽江淮汽车集团股份有限公司 | 电动汽车电池寿命预测方法、装置、设备及存储介质 |
CN111426952A (zh) * | 2019-01-10 | 2020-07-17 | 郑州宇通客车股份有限公司 | 一种锂离子电池的寿命预测方法 |
CN111999665A (zh) * | 2020-10-29 | 2020-11-27 | 北京航空航天大学 | 一种基于微观机理汽车驾驶工况锂离子电池老化试验方法 |
CN112595980A (zh) * | 2020-12-17 | 2021-04-02 | 北京海博思创科技股份有限公司 | 电池储能系统寿命的预测方法、装置及设备 |
CN112731164A (zh) * | 2020-12-21 | 2021-04-30 | 惠州亿纬锂能股份有限公司 | 一种电池寿命评估方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9707395B2 (en) * | 2014-01-16 | 2017-07-18 | Boston Scientific Neuromodulation Corporation | Circuitry for charging a battery in an implantable medical device in accordance with historical parameters impacting battery capacity |
US9864016B2 (en) * | 2014-10-31 | 2018-01-09 | GM Global Technology Operations LLC | Battery system pack life estimation systems and methods |
US10422835B2 (en) * | 2015-10-27 | 2019-09-24 | Nec Corporation | Innovative framework combining cycling and calendar aging models |
CN106515478B (zh) * | 2016-10-27 | 2017-10-20 | 合肥工业大学 | 电动汽车剩余行驶里程在线预测方法及装置 |
-
2021
- 2021-06-18 CN CN202110680465.6A patent/CN113406522B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012127938A (ja) * | 2010-11-24 | 2012-07-05 | Mitsubishi Electric Corp | 蓄電装置の劣化監視方法、及びその劣化監視装置 |
JP2017502820A (ja) * | 2014-01-16 | 2017-01-26 | ボストン サイエンティフィック ニューロモデュレイション コーポレイション | 充電式バッテリを有する埋め込み可能医療デバイスに対する寿命末期の決定及び予想 |
CN105904981A (zh) * | 2016-04-07 | 2016-08-31 | 北京现代汽车有限公司 | 一种电动汽车续航里程估计控制方法、装置及整车控制器 |
CN107317057A (zh) * | 2016-04-27 | 2017-11-03 | 宝沃汽车(中国)有限公司 | 一种动力电池运行寿命预测和预测装置 |
CN107064800A (zh) * | 2016-11-29 | 2017-08-18 | 北京交通大学 | 锂离子电池剩余使用寿命的实时预测方法 |
CN108710084A (zh) * | 2018-06-01 | 2018-10-26 | 王泽霖 | 一种电动汽车电源监控及能源管理系统 |
CN111426952A (zh) * | 2019-01-10 | 2020-07-17 | 郑州宇通客车股份有限公司 | 一种锂离子电池的寿命预测方法 |
CN109747428A (zh) * | 2019-02-27 | 2019-05-14 | 合肥国轩高科动力能源有限公司 | 一种电动汽车电池包剩余容量估算方法 |
CN110470991A (zh) * | 2019-08-23 | 2019-11-19 | 江西优特汽车技术有限公司 | 一种动力电池循环寿命评测方法 |
CN110901470A (zh) * | 2019-11-29 | 2020-03-24 | 安徽江淮汽车集团股份有限公司 | 电动汽车电池寿命预测方法、装置、设备及存储介质 |
CN111999665A (zh) * | 2020-10-29 | 2020-11-27 | 北京航空航天大学 | 一种基于微观机理汽车驾驶工况锂离子电池老化试验方法 |
CN112595980A (zh) * | 2020-12-17 | 2021-04-02 | 北京海博思创科技股份有限公司 | 电池储能系统寿命的预测方法、装置及设备 |
CN112731164A (zh) * | 2020-12-21 | 2021-04-30 | 惠州亿纬锂能股份有限公司 | 一种电池寿命评估方法 |
Non-Patent Citations (4)
Title |
---|
Lithium Iron Phosphate - Assessment of Calendar Life and Change of Battery Parameters;N. Omar等;《2014 IEEE Vehicle Power and Propulsion Conference (VPPC)》;20150209;全文 * |
Remaining useful life estimation of Lithium-ion battery based on interacting multiple model particle filter and support vector regression;Sai Li等;《Reliability Engineering and System Safety》;20210225;全文 * |
磷酸铁锂动力电池日历寿命加速测试与拟合;王芳等;《电源技术》;20150120(第01期);全文 * |
锂离子电池剩余寿命预测方法研究综述;郑文芳等;《计算机测量与控制》;20201231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113406522A (zh) | 2021-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schaltz et al. | Incremental capacity analysis applied on electric vehicles for battery state-of-health estimation | |
CN108919129B (zh) | 一种时变工况下动力电池的寿命预测方法 | |
CN113406522B (zh) | 一种电动车锂电池系统寿命预测评估方法及系统 | |
CN111239630A (zh) | 一种储能电池寿命预测方法及管理系统 | |
CN110795851B (zh) | 一种考虑环境温度影响的锂离子电池建模方法 | |
CN112014735B (zh) | 一种基于全寿命周期的电芯老化寿命预测方法及装置 | |
CN109358293B (zh) | 基于ipf的锂离子电池soc估计方法 | |
CN112198434B (zh) | 电池容量衰减模型参数的辨识方法、系统、设备、装置及介质 | |
Soto et al. | Analysis of the main battery characterization techniques and experimental comparison of commercial 18650 Li-ion cells | |
CN111562509B (zh) | 一种用于确定退役动力电池剩余寿命的方法及系统 | |
CN113109729B (zh) | 基于加速老化试验与实车工况的车用动力电池soh评估方法 | |
CN110045291A (zh) | 一种锂电池容量估计方法 | |
Kataoka et al. | Battery state estimation system for automobiles | |
CN111257770B (zh) | 一种电池包功率估算方法 | |
CN113406520A (zh) | 针对现实新能源汽车的电池健康状态估计方法 | |
CN115407206A (zh) | 一种基于容量累积的soh自适应估算方法 | |
CN115754372A (zh) | 一种变温度循环的锂离子电池老化测试方法 | |
CN115219918A (zh) | 一种基于容量衰退组合模型的锂离子电池寿命预测方法 | |
CN117388737A (zh) | 电池健康状态的评估方法、装置、设备及存储介质 | |
US20230314528A1 (en) | Method of Estimation of Battery Degradation | |
Chen et al. | Online estimation of state of power for lithium-ion battery considering the battery aging | |
CN116679213A (zh) | 基于集成深度学习的电动汽车动力电池soh估算方法 | |
CN113687256B (zh) | 一种评估单体自放电率对电池系统一致性影响的方法 | |
Liu et al. | Economic assessment of Lithium-Ion batteries in terms of V2G utilisation | |
CN114994559A (zh) | 一种针对运营车辆的动力电池循环寿命的试验方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |