CN103869256B - 一种基于交流阻抗测试动力锂离子电池soh的估算方法 - Google Patents

一种基于交流阻抗测试动力锂离子电池soh的估算方法 Download PDF

Info

Publication number
CN103869256B
CN103869256B CN201410119748.3A CN201410119748A CN103869256B CN 103869256 B CN103869256 B CN 103869256B CN 201410119748 A CN201410119748 A CN 201410119748A CN 103869256 B CN103869256 B CN 103869256B
Authority
CN
China
Prior art keywords
battery
soh
curve
impedance
charge transfer
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
Application number
CN201410119748.3A
Other languages
English (en)
Other versions
CN103869256A (zh
Inventor
李伟善
张远明
胡佳娜
蔡燕凤
邱显焕
熊恋文
陈彦霖
王英博
郭均盛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China Normal University
Original Assignee
South China Normal University
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 South China Normal University filed Critical South China Normal University
Priority to CN201410119748.3A priority Critical patent/CN103869256B/zh
Publication of CN103869256A publication Critical patent/CN103869256A/zh
Application granted granted Critical
Publication of CN103869256B publication Critical patent/CN103869256B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Secondary Cells (AREA)

Abstract

本发明公开了一种基于交流阻抗测试动力锂离子电池SOH的估算方法,包含以下步骤:以1C到0.5C之间任一恒定电流对电池进行充电,再以0.1C到0.05C之间任一恒定电流对电池进行放电;在电池放电到0电荷状态时对电池进行交流阻抗测试,并拟合出的电池的电荷转移阻抗;在Origin9中拟合出容量衰减曲线以及电荷转移阻抗曲线;根据电荷转移阻抗递增曲线对出容量衰减曲线进行修正,建立出电池健康状况SOH随循环圈数变化的估算函数,计算出可供所述同类型电池健康状况的判断依据。本发明的方法,加入电荷转移阻抗对SOH进行估算,估算结果更加精确,且不需要拆解电池。

Description

一种基于交流阻抗测试动力锂离子电池SOH的估算方法
技术领域
本发明涉及锂离子电池健康状态估算领域,具体涉及一种基于交流阻抗测试动力锂离子电池SOH的估算方法。
背景技术
随着电动汽车(EDVs)的迅猛发展,包括纯电动汽车(BEV)、混合动力汽车(HEV)和插电式混合动力汽车(PHEV),电池技术得到世界范围的重视和关注。因此,基于现有的电池技术,一个重要的研究领域是提高电池在电动汽车中的使用效率和寿命。目前动力锂离子电池作为纯电动汽车最优的解决方案,其仍然有潜在过充电过放电引发安全事故的风险。准确的估算电池的健康状态SOH(State of Charge)是一个保证电动车安全运行的关键技术。伴随电池循环圈数增加,电池容量衰减,内阻增大,电池的健康状态SOH能代表电池的老化程度。
现有电池的健康状态SOH有多种定义方法,例如通过电池当前容量和初始容量的比值定义,即其中C为电池当前容量,C0为电池初始容量。目前也有根据电池内阻、功率密度和能量密度来定义SOH的算法。
电池容量衰减的主要原因在于电池正极活性材料的损失,电池负极材料的损失,电池电解液的损失和电池隔膜的老化所导致。不同类型电池的在不同的循环充放电制式的情况下电池的容量衰减的机制也不相同。目前阶段在对电池健康状态SOH估算的时候,仅仅估算电池容量衰减的幅度,并无法从电池容量衰减的机理来判断电池内部阻抗的变化。
现阶段对电池容量衰减机理的深入研究,通常将电池进行拆解,将正负极活性材料进行SEM(Scanning Electron Microscope,扫描电子显微镜),EDS(Energy DispersiveSpectroscopy)能量色散谱测试,分析电极活性材料的衰减情况,从而推断电池的衰减机理。这种方法在实际电动汽车所使用的动力锂离子电池并没有可行性,拆解电池的方法使得电池无法使用。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种基于交流阻抗测试动力锂离子电池SOH的估算方法。
本发明的目的通过以下的技术方案实现:
一种基于交流阻抗测试动力锂离子电池SOH的估算方法,包含以下顺序的步骤:
S1.以1C到0.5C之间任一恒定电流对电池进行充电,再以0.1C到0.05C之间任一恒定电流对电池进行放电,并记录每圈放电过程放出的电量;
S2.在电池放电到0电荷状态时对电池进行交流阻抗测试,取中频0.25Hz~1.2KHz交流阻抗曲线在NOVA1.1中拟合出的电池的电荷转移阻抗;
S3.根据每圈放电电量、电荷转移阻抗、圈数作为估算函数的输入参数,在0rigin9中拟合出容量衰减曲线以及电荷转移阻抗曲线;
S4.根据电荷转移阻抗递增曲线对出容量衰减曲线进行修正,建立出电池健康状况SOH随循环圈数变化的估算函数,计算出可供所述同类型电池健康状况的判断依据。
所述的估算函数,定义为SOH=1-C·Rct,其中C为容量衰减曲线,Rct为电池的电荷转移阻抗。
所述的基于交流阻抗测试动力锂离子电池SOH的估算方法,还包括步骤S5,验证估算方法的准确性:
所述电池的剩余容量和所述真实值之间的标准误差符合以下回归方程的显著性检验且误差小于7%:
其中,为所述电池的剩余容量的估算值,CAhi为所述电池的剩余容量的真实值,为所述电池的剩余容量的估算值,i为循环圈数,表示电池循环圈数从1到i的i个采集点;
误差技术满足方程:
本发明与现有技术相比,具有如下优点和有益效果:
本发明在传统电池健康状态SOH估算中加入电荷转移阻抗进行修正。电荷转移阻抗有别于目前的对于电池内阻的修正方法,将内阻中的随循环圈数变化不大的纯欧姆阻抗剔除,对电池健康状态SOH估算更加精确,同时不需要拆解电池,保持电池的完整性。
附图说明
图1为本发明所述的一种基于交流阻抗测试动力锂离子电池SOH的估算方法的流程图;
图2为图1所述方法中的电池容量变化曲线图;
图3为图1所述方法中电池电荷转移阻抗与电池循环次数关系曲线;
图4为图1所述方法中经过1700次充放电循环后点电池电荷转移阻抗的估计值和真实值的对比图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
如图1,一种基于交流阻抗测试动力锂离子电池SOH的估算方法,包含以下顺序的步骤:
S1.以1C到0.5C之间任一恒定电流对电池进行充电,再以0.1C到0.05C之间任一恒定电流对电池进行放电,并记录每圈放电过程放出的电量;
S2.在电池放电到0电荷状态时对电池进行交流阻抗测试,取中频0.25Hz~1.2KHz交流阻抗曲线在NOVA1.1中拟合出的电池的电荷转移阻抗;
S3.根据每圈放电电量、电荷转移阻抗、圈数作为估算函数的输入参数,在Origin9中拟合出容量衰减曲线以及电荷转移阻抗曲线;
S4.根据电荷转移阻抗递增曲线对出容量衰减曲线进行修正,建立出电池健康状况SOH随循环圈数变化的估算函数,计算出可供所述同类型电池健康状况的判断依据;所述的估算函数,定义为SOH=1-C·Rct,其中C为容量衰减曲线,Rct为电池的电荷转移阻抗。
S5.验证估算方法的准确性:
所述电池的剩余容量和所述真实值之间的标准误差符合以下回归方程的显著性检验且误差小于7%:
其中,为所述电池的剩余容量的估算值,CAhi为所述电池的剩余容量的真实值,为所述电池的剩余容量的估算值,i为循环圈数,表示电池循环圈数从1到i的i个采集点;
误差技术满足方程:
通常实施使用中的电动车动力型锂离子电池,其放电制式由实际电动车运行方式和路况决定,从而非常复杂没有规律。而放电到0电量时静置一段时间电池内部电化学系统区域稳定,因此,通过对电池交流阻抗测试可以非侵入式的检测电池的内部信息,进而分析电池容量的变化情况。
下面进一步具体实施例对技术方案进行说明,如图1、2、3、4,一种基于交流阻抗测试动力锂离子电池SOH的估算方法,包含以下顺序的步骤:
步骤S1:以恒定0.5C电流对电池进行充电,电池电压稳定在4.2V时,充电电流减小到0.01C视为电池充电充满;以0.05C对电池进行放电,电池电压降低到2.5V视为电池放电到0电量,同时记录放出的电量;
考虑到电池正极和负极材料的多样性,在本发明的实施例中以应用最广泛的钴酸锂为正极,石墨为负极的18650锂离子电池为例进行说明;需要说明的是,本发明实施例亦可以使用其他不同的正极和负极材料的电池;
图2为的电池容量变化曲线图;根据本发明一个实施例的钴酸锂为正极,石墨为负极的18650动力锂离子电池不同衰减情况下的容量变化曲线:其中一条曲线为电池经历1700次充放电循环过程记录的电池容量;另外一条曲线为电池容量变化曲线在Origin9中根据玻尔兹曼函数进行修正,得到电池容量衰减公式:
本示例公式为
实际不同型号、品牌电池得出电池容量衰减公式不同,本公式旨在对测试中同型号电池负责;其中,C为电池实际容量,x为循环圈数;
步骤S2:获取1,5,10,15,20,25……循环圈数中电池0电量时的中频0.25Hz~1.2KHz交流阻抗曲线,通过NOVA1.1电化学软件拟合出电荷转移阻抗,无损的非侵入性的得到电池内部的变化情况,同时在Origin9中拟合出电荷转移阻抗与循环圈数关系曲线;如图3,其中一条曲线为电荷转移阻抗与循环圈数关系曲线在Origin9中根据玻尔兹曼函数进行修正,得到电池电荷转移阻抗递增公式;即实线为电池经历1700次充放电循环过程记录的电荷转移阻抗,虚线为Origin9软件中拟合出电池容量衰减曲线,得出电池电荷转移阻抗递增公式
本示例公式为
实际不同型号、品牌电池得出电池容量衰减公式不同,本公式旨在对测试中同型号电池负责。其中,Rct为电池的电荷转移阻抗,x为循环圈数;
步骤S3:将电荷转移阻抗与循环圈数关系曲线与容量衰减曲线结合分析,在Origin9中定义电池健康状态SOH与循环圈数的曲线关系;所述曲线关系定义公式为
SOH=1-C·Rct (3)
可见图4,红色为电池SOH预测曲线,黑色为电池实际SOH值;图4中真实值为对应圈数的电池容量与电荷转移阻抗通过定义公式(3)计算值,虚线为拟合曲线经过定义公式(3)得到的预测值。预测值满足验证步骤S5。
步骤S4:进行验证,将同种型号18650动力锂离子电池在上述相同充放电制式条件下测试,由所述模型预测该测试电池的估算健康状态SOH值,所述电池的剩余容量和所述真实值之间的标准误差符合以下回归方程的显著性检验且误差小于7%:
其中,为所述电池的剩余容量的估算值,CAhi为所述电池的剩余容量的真实值,为所述电池的剩余容量的估算值,i为循环圈数,表示电池循环圈数从1到i的i个采集点;
误差技术满足方程:
从上述实施例中,可以发现对锂离子动力电池进行“充电-放电”的模式循环,随着循环次数的增加,SOH值线性递减,说明本发明方法估算的SOH能够得到相对准确的电池健康状态,相比现有技术的评估方法更准确,这是由于循环次数不同(充放电电流、电压参数设置一致),之后又以相同的倍率放电,再静置,中频0.25Hz~1.2KHz交流阻抗曲线通过NOVA1.1电化学软件拟合出电荷转移阻抗,并且电荷转移阻抗与循环次数呈现线性关系,即循环次数越多,电荷转移阻抗随循环圈数增加而线性递增,这与电池的内部元件的老化损失有关;锂离子动力电池的老化损失原因主要体现在电池内阻、充放电电流大小、极片制作工艺、涂料是否均匀、压实密度有关。在采用同一个电池的情况下,可以忽略充放电电流大小、极片制作工艺、涂布是否均匀、压实密度对电池极化的影响,所以此状态下电池的老化损失主要体现在电池内阻的增大。电池的内阻包括欧姆电阻和电荷转移阻抗,电池的欧姆电阻主要由电极材料、电解液、隔膜电阻及各部分零件的接触电阻组成,与电池的尺寸、结构、装配等有关,内阻小的电池的放电能力强,内阻大的电池放电能力弱,说明电池健康状态差,电池电荷转移阻抗大,说明电池内部产生较大阻抗。而Rb纯欧姆阻抗在电池循环寿命中变化不大,通过NOVA1.1电化学软件对中频0.25Hz~1.2KHz交流阻抗曲线中的Rb纯欧姆阻抗进行剔除,单独找出电池电荷转移阻抗,对电池的内部阻抗值作为加以修正。
Nova软件是瑞士万通公司Autolab电化学工作站配套软件,具有拟合交流阻抗曲线功能,OriginOrigin为OriginLab公司出品的较流行的专业函数绘图软件,是公认的简单易学、操作灵活、功能强大的软件,既可以满足一般用户的制图需要,也可以满足高级用户数据分析、函数拟合的需要。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (1)

1.一种基于交流阻抗测试动力锂离子电池SOH的估算方法,其特征在于,包含以下顺序的步骤:
S1.以1C到0.5C之间任一恒定电流对电池进行充电,再以0.1C到0.05C之间任一恒定电流对电池进行放电,并记录每圈放电过程放出的电量;
S2.在电池放电到0电荷状态时对电池进行交流阻抗测试,取中频0.25Hz~1.2KHz交流阻抗曲线在NOVA1.1中拟合出的电池的电荷转移阻抗;
S3.根据每圈放电电量、电荷转移阻抗、圈数作为估算函数的输入参数,在Origin9中拟合出容量衰减曲线以及电荷转移阻抗曲线;
S4.根据电荷转移阻抗递增曲线对容量衰减曲线进行修正,建立电池健康状况SOH随循环圈数变化的估算函数,计算出可供同类型电池健康状况的判断依据;
所述的估算函数,定义为SOH=1-C·Rct,其中C为容量衰减曲线,Rct为电池的电荷转移阻抗。
CN201410119748.3A 2014-03-27 2014-03-27 一种基于交流阻抗测试动力锂离子电池soh的估算方法 Active CN103869256B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410119748.3A CN103869256B (zh) 2014-03-27 2014-03-27 一种基于交流阻抗测试动力锂离子电池soh的估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410119748.3A CN103869256B (zh) 2014-03-27 2014-03-27 一种基于交流阻抗测试动力锂离子电池soh的估算方法

Publications (2)

Publication Number Publication Date
CN103869256A CN103869256A (zh) 2014-06-18
CN103869256B true CN103869256B (zh) 2017-02-15

Family

ID=50907997

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410119748.3A Active CN103869256B (zh) 2014-03-27 2014-03-27 一种基于交流阻抗测试动力锂离子电池soh的估算方法

Country Status (1)

Country Link
CN (1) CN103869256B (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502844A (zh) * 2014-11-18 2015-04-08 国家电网公司 一种基于交流阻抗的动力锂电池劣化程度诊断方法
CN104569841A (zh) * 2014-12-26 2015-04-29 国家电网公司 电池组的老化检测方法和装置
CN104678317A (zh) * 2015-03-05 2015-06-03 普天新能源车辆技术有限公司 动力锂电池容量衰减的检测方法和装置
CN107144790A (zh) * 2016-03-01 2017-09-08 上海卡耐新能源有限公司 一种预测锂离子电池循环寿命的方法
CN105652215A (zh) * 2016-04-07 2016-06-08 苏州协鑫集成科技工业应用研究院有限公司 蓄电池健康状态的检测方法及装置
CN105929336B (zh) * 2016-05-04 2019-03-26 合肥国轩高科动力能源有限公司 一种动力锂离子电池健康状态估算方法
CN106855610B (zh) * 2016-12-30 2021-09-10 中国电力科学研究院 钛酸锂电池健康状态估算方法
JP7020273B2 (ja) * 2018-04-26 2022-02-16 トヨタ自動車株式会社 電池情報処理システム、組電池、電池モジュールの容量算出方法、および、組電池の製造方法
CN109633452A (zh) * 2018-12-24 2019-04-16 广东省智能制造研究所 一种电池健康度检测方法及检测装置
JP7169519B2 (ja) * 2018-12-26 2022-11-11 トヨタ自動車株式会社 電池容量推定方法、および電池容量推定システム
CN110426639B (zh) * 2019-07-24 2022-09-23 中国电力科学研究院有限公司 一种基于动态阻抗谱的锂离子电池寿命预测方法及系统
CN110398693A (zh) * 2019-07-29 2019-11-01 国网河南省电力公司电力科学研究院 一种退役磷酸铁锂单体电池状态快速评价方法
CN110927585A (zh) * 2019-12-31 2020-03-27 四川长虹电器股份有限公司 一种基于自循环校正的锂电池soh估算系统及方法
CN113138340B (zh) * 2020-01-17 2022-11-11 华为技术有限公司 电池等效电路模型的建立方法、健康状态估算方法及装置
CN112363074A (zh) * 2020-11-10 2021-02-12 中国电力科学研究院有限公司 电池储能适用性的检测评估方法、系统、设备和存储介质
CN112611973A (zh) * 2020-11-27 2021-04-06 中国电力科学研究院有限公司 一种电池健康状态估算方法、系统、设备和存储介质
CN112698230A (zh) * 2020-12-02 2021-04-23 国网上海市电力公司 一种锂离子电池健康状态的动态阻抗快速测量方法
CN113687253B (zh) * 2021-08-23 2023-06-30 蜂巢能源科技有限公司 分析电芯内部组件阻抗的方法
CN116252682B (zh) * 2023-05-16 2023-07-21 车百中汽科技(北京)有限公司 一种新能源汽车动力电池健康度的确定方法及系统
CN118068204B (zh) * 2024-04-17 2024-07-09 贲安能源科技(上海)有限公司 一种利用阻抗分析钠离子电池衰减原因及寿命预测的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101208847A (zh) * 2005-06-30 2008-06-25 Lg化学株式会社 估算电池健康状态的方法和装置
CN102156265A (zh) * 2011-03-16 2011-08-17 深圳市派司德科技有限公司 一种电池健康状态测试装置及其方法
CN102445663A (zh) * 2011-09-28 2012-05-09 哈尔滨工业大学 一种电动汽车电池健康状态估算的方法
CN102854394A (zh) * 2012-08-31 2013-01-02 华南师范大学 锂离子电池健康状态估算系统及使用该系统对锂离子电池健康状态估算的方法
CN102866361A (zh) * 2012-08-31 2013-01-09 惠州市亿能电子有限公司 一种电池组soh在线估算方法
CN103163480A (zh) * 2013-03-29 2013-06-19 长城汽车股份有限公司 锂电池健康状态的评估方法
CN103344921A (zh) * 2013-07-08 2013-10-09 华南师范大学 锂离子动力电池健康状态评估系统及方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2943794B1 (fr) * 2009-03-24 2011-05-06 Saft Groupe Sa Procede de determination de l'etat de sante d'une batterie
US8738310B2 (en) * 2010-12-08 2014-05-27 Paul Swanton Automatic determination of baselines for battery testing
ES2629192T3 (es) * 2011-08-15 2017-08-07 Shawn P. Kelly Aparato y método para supervisión precisa de estado de salud (SoH) de dispositivo de energía

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101208847A (zh) * 2005-06-30 2008-06-25 Lg化学株式会社 估算电池健康状态的方法和装置
CN102156265A (zh) * 2011-03-16 2011-08-17 深圳市派司德科技有限公司 一种电池健康状态测试装置及其方法
CN102445663A (zh) * 2011-09-28 2012-05-09 哈尔滨工业大学 一种电动汽车电池健康状态估算的方法
CN102854394A (zh) * 2012-08-31 2013-01-02 华南师范大学 锂离子电池健康状态估算系统及使用该系统对锂离子电池健康状态估算的方法
CN102866361A (zh) * 2012-08-31 2013-01-09 惠州市亿能电子有限公司 一种电池组soh在线估算方法
CN103163480A (zh) * 2013-03-29 2013-06-19 长城汽车股份有限公司 锂电池健康状态的评估方法
CN103344921A (zh) * 2013-07-08 2013-10-09 华南师范大学 锂离子动力电池健康状态评估系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于开路电压回升速率和交流阻抗相结合的一种锂离子电池SOH算法;黄伟昭等;《新能源进展》;20140228;43-47 *
锂离子电池健康预测方法;高安同等;《重庆理工大学学报(自然科学)》;20140228;102-105 *

Also Published As

Publication number Publication date
CN103869256A (zh) 2014-06-18

Similar Documents

Publication Publication Date Title
CN103869256B (zh) 一种基于交流阻抗测试动力锂离子电池soh的估算方法
Fly et al. Rate dependency of incremental capacity analysis (dQ/dV) as a diagnostic tool for lithium-ion batteries
Stroe et al. Lithium-ion battery state-of-health estimation using the incremental capacity analysis technique
Abdel-Monem et al. Influence analysis of static and dynamic fast-charging current profiles on ageing performance of commercial lithium-ion batteries
Xiong et al. A systematic model-based degradation behavior recognition and health monitoring method for lithium-ion batteries
Tian et al. Fractional-order model-based incremental capacity analysis for degradation state recognition of lithium-ion batteries
Galeotti et al. Performance analysis and SOH (state of health) evaluation of lithium polymer batteries through electrochemical impedance spectroscopy
Wang et al. State of charge-dependent polynomial equivalent circuit modeling for electrochemical impedance spectroscopy of lithium-ion batteries
Barcellona et al. Aging effect on the variation of Li-ion battery resistance as function of temperature and state of charge
Monem et al. Lithium-ion batteries: Evaluation study of different charging methodologies based on aging process
CN107145628B (zh) 基于电化学-热耦合模型的预测锂电池循环寿命的方法
Eddahech et al. Behavior and state-of-health monitoring of Li-ion batteries using impedance spectroscopy and recurrent neural networks
CN109143108A (zh) 一种基于电化学阻抗谱的锂离子电池soh的估计方法
Leng et al. Hierarchical degradation processes in lithium-ion batteries during ageing
CN104076293B (zh) 基于观测器的锂电池soc估算误差的定量分析方法
CN109154636A (zh) 一种储能装置监测技术
JP6337233B2 (ja) 電池の評価方法及び電池特性評価装置
CN108614221B (zh) 一种锂离子电池化成工序的评价方法
CN114865117B (zh) 锂离子电池电极嵌锂量检测方法、装置及电池管理系统
Qu et al. A joint grey relational analysis based state of health estimation for lithium ion batteries considering temperature effects
CN114609523A (zh) 一种电池容量的在线检测方法、电子设备及存储介质
CN111722120A (zh) 一种锂离子电池可逆锂消耗量的评价方法及系统
CN115219918A (zh) 一种基于容量衰退组合模型的锂离子电池寿命预测方法
Cui et al. Prediction model and principle of end-of-life threshold for lithium ion batteries based on open circuit voltage drifts
Wang et al. Parameter sensitivity analysis and parameter identifiability analysis of electrochemical model under wide discharge rate

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant