CN112098846B - 一种退役动力锂电池的混合建模方法 - Google Patents
一种退役动力锂电池的混合建模方法 Download PDFInfo
- Publication number
- CN112098846B CN112098846B CN202010827778.5A CN202010827778A CN112098846B CN 112098846 B CN112098846 B CN 112098846B CN 202010827778 A CN202010827778 A CN 202010827778A CN 112098846 B CN112098846 B CN 112098846B
- Authority
- CN
- China
- Prior art keywords
- battery
- model
- equivalent circuit
- voltage
- retired power
- 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
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/367—Software therefor, e.g. for battery testing using modelling or look-up tables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/385—Arrangements for measuring battery or accumulator variables
- G01R31/387—Determining ampere-hour charge capacity or 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/389—Measuring internal impedance, internal conductance or related variables
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种退役动力锂电池的混合建模方法,属于动力锂电池应用领域。该方法包括建立退役动力锂电池的一阶等效电路模型,并获取模型参数;采用偏最小二乘法对一阶等效电路模型形成的误差进行二次拟合,建立退役动力锂电池的偏差补偿模型;联合一阶等效电路模型和偏差补偿模型,建立退役动力锂电池的混合模型。本发明对外特性已经发生一定程度变化的退役电池,通过一阶等效电路模型描述电池的基本外特性,建立基于偏最小二乘的偏差补偿模型,能够有效提高退役动力电池建模的精度,便于后续电池管理系统实现电池状态参数的准确估计及电池组的高效能量管理,有助于退役动力电池的梯次利用。
Description
技术领域
本发明涉及动力锂电池应用领域,具体涉及一种退役动力锂电池的混合建模方法。
背景技术
锂电池已经成为了当前新能源汽车最主要的储能元件。然而,伴随着新能源汽车的日常使用,动力锂电池会不可避免地逐渐老化。从车辆上退役下来的锂电池,仍然具备在性能要求较低的场合继续使用的潜在价值,即退役电池仍然具备梯次利用价值。电池模型是梯次利用环节中进行系统能量管理的重要基础,电池核心状态参数(如荷电状态(Stateof Charge,SOC)、健康状态(State of Health,SOH)等)的获取也需要精确的电池模型。
相较于全新电池,退役锂电池经过了新能源汽车的一次使用,其内部已经发生一定程度老化。但是,锂电池初次使用过程中历史数据是缺失的,这就导致无法直接使用大量数据建立基于数据驱动的电池模型。传统的经验模型仅通过数学函数来描述电池的外特性,显然无法满足退役电池建模的复杂需求。电化学模型的计算量较大,参数多且难以获取。常用的等效电路模型为提高建模精度,一般需要增加并联的电阻电容网络数目,应用于退役动力锂电池时,很可能会造成模型结构复杂、参数在线辨识困难等问题。
发明内容
针对现有技术中的上述不足,本发明提供了一种退役动力锂电池的混合建模方法。
为了达到上述发明目的,本发明采用的技术方案为:
一种退役动力锂电池的混合建模方法,包括以下步骤:
S1、建立退役动力锂电池的一阶等效电路模型,并获取模型参数;
S2、采用偏最小二乘法对一阶等效电路模型形成的误差进行二次拟合,建立退役动力锂电池的偏差补偿模型;
S3、联合一阶等效电路模型和偏差补偿模型,建立退役动力锂电池的混合模型。
优选地,所述退役动力锂电池的一阶等效电路模型表示为:
Ut=Uocv-U1-It·R0
其中,Ut为电池的端电压,It为电池的充放电电流,Uocv为电池开路电压,U1为电容电阻网络的端电压,R0为电池内阻,R1、C1为电容电阻网络对应的电阻、电容值。
优选地,所述步骤S1中获取一阶等效电路模型参数的方法为:
S11、对一阶等效电路模型进行离线电流脉冲测试,得到电压响应曲线;
S12、计算一阶等效电路模型在测试条件下的充放电内阻;
S13、分别提取电流脉冲充放电后电池静置时间段的电压测量值,并拟合电压曲线;
S14、采用非线性最小二乘法计算电压曲线中的拟合参数;
S15、根据步骤S14得到的拟合参数计算一阶等效电路模型参数。
优选地,所述一阶等效电路模型在测试条件下的充放电内阻的计算公式为:
其中,R0_dis为放电内阻,R0_ch为充电内阻,分别为电压响应曲线对应电流脉冲放电时刻t1、t0的电压测量值,分别为电压响应曲线对应电流脉冲充电时刻t5、t4的电压测量值,Ich、Idis分别为测试过程的充放电电流倍率。
优选地,所述一阶等效电路模型参数的计算公式为:
其中,τ为电容电阻网络的时间常数,R1_dis为电容电阻网络对应的放电电阻,R1_ch为电容电阻网络对应的充电电阻,α、B为电压曲线的拟合参数,tdis为电压响应曲线对应电流脉冲放电时刻t2、t1之间的时间段,tch为电压响应曲线对应电流脉冲放电时刻t6、t5之间的时间段。
优选地,所述步骤S2具体为:
根据一阶等效电路模型的输出与电池端电压测量值之间的偏差,采用偏最小二乘法进行二次拟合,计算得到偏最小二乘偏差补偿模型,建立基于偏最小二乘的退役动力锂电池的偏差补偿模型。
优选地,所述退役动力锂电池的偏差补偿模型表示为:
ΔU=b1+b2·SOC+b3·It
其中,SOC为偏最小二乘法训练窗口内电池荷电状态的估计值。
优选地,所述退役动力锂电池的混合模型表示为:
Uhybrid=Ut+ΔU=Uocv-U1-It·R0+b1+b2·SOC+b3·It
其中,Uhybrid表示混合模型的输出电压值。
本发明具有以下有益效果:
本发明对外特性已经发生一定程度变化的退役电池,通过一阶等效电路模型描述电池的基本外特性,建立基于偏最小二乘的偏差补偿模型,能够有效提高退役动力电池建模的精度,便于后续电池管理系统实现电池状态参数的准确估计及电池组的高效能量管理,有助于退役动力电池的梯次利用。
附图说明
图1为本发明的退役动力锂电池的混合建模方法流程示意图;
图2为本发明实施例中一阶等效电路模型结构示意图;
图3为本发明实施例中离线电池参数辨识工况示意图;其中图(a)为电流脉冲曲线,图(b)为电压响应曲线;
图4为本发明实施例中偏最小二乘偏差补偿模型的建立示意图;
图5为本发明实施例中退役动力锂电池混合模型示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
由于退役电池内部已经产生了不同程度的老化,外特性与全新电池相比不尽相同,为此,传统锂电池建模方法应对退役电池时有一定难度。电池模型的参数在线辨识方法对模型精度的提升非常有限,且易受测量噪声的干扰。而精确的电池模型一直是进行电池组能量管理所需的关键技术。
为此,本发明为提高梯次利用中退役动力锂电池的建模精度,方便退役动力锂电池的梯次利用,针对电池状态估计、电池组能量管理所需的精确电池模型展开分析,提出了一种联合不同电池模型的混合建模方法。该方法综合利用了锂电池等效电路模型与数据驱动建模方法的优势,选取一阶等效电路模型作为混合模型的基本框架,而等效电路模型的参数可由不同条件下的电流脉冲离线测量获得;在此基础上,计算特定窗口内一阶等效电路模型的误差,将该误差、电流、SOC估计值作为训练样本,建立退役电池的偏差补偿模型;最后,通过联合一阶等效电路模型与基于最小二乘法的偏差补偿模型,得到了退役动力锂电池混合模型。本发明提供的混合建模方法,对于促进梯次利用中退役动力锂电池的状态获取及电池组能量管理均具有重要意义。
如图1所示,本发明实施例提供的一种退役动力锂电池的混合建模方法,包括以下步骤S1至S3:
S1、建立退役动力锂电池的一阶等效电路模型,并获取模型参数;
在本实施例中,本发明综合考虑电池模型的复杂度与精度,选取一阶等效电路模型作为电池建模的主要框架,利用离线电流脉冲测试,建立初步的退役电池等效电路模型。
本发明建立的一阶等效电路模型主要由电压源、电阻及电阻-电容网络组成,如图2所示。根据一阶等效电路模型的结构,得到退役动力锂电池的一阶等效电路模型表示为:
Ut=Uocv-U1-It·R0
其中,Ut为电池的端电压,It为电池的充放电电流,Uocv为电池开路电压,U1为电容电阻网络的端电压,R0为电池内阻,R1、C1为电容电阻网络对应的电阻、电容值。式中参数Uocv、R0、R1、C1均可以离线测量计算得到,离线电池参数测试工况如图3所示。
本发明获取一阶等效电路模型参数的方法为:
S11、对一阶等效电路模型进行离线电流脉冲测试,得到电压响应曲线;
如图3所示工况为电流脉冲测试及相应的电压响应曲线,t0为电流脉冲跳变起始时刻,t1为电流脉冲跳变至预定放电电流幅值的时刻,t2为电流脉冲结束时刻,t3为电流脉冲跳变为0的时刻,t4为电池静置1小时后的结束时刻,t5为电流跳变至预定充电电流幅值的时刻,t6为电流脉冲结束时刻,t7为电流脉冲跳变为0的时刻,t8为电池静置1小时后的结束时刻,其中:t0~t3为持续时间20秒的电流脉冲放电,t3~t4为持续时间1小时的电池静置,t4~t7为持续时间20秒的电流脉冲充电,t7~t8为持续时间1小时的电池静置。为了获得不同SOC、温度、电流倍率下的电池模型参数,可以在多个不同条件下施加图3所示电流脉冲测试工况,收集电压的响应数据。
S12、计算一阶等效电路模型在测试条件下的充放电内阻;
充放电内阻的计算公式为:
其中,R0_dis为放电内阻,R0_ch为充电内阻,分别为电压响应曲线对应电流脉冲放电时刻t1、t0的电压测量值,分别为电压响应曲线对应电流脉冲充电时刻t5、t4的电压测量值,Ich、Idis分别为测试过程的充放电电流倍率。
S13、分别提取电流脉冲充放电后电池静置时间段的电压测量值,并拟合电压曲线;
本发明分别提取电流脉冲充放电后电池静置时间段t3~t4及t7~t8内的电压测量值,并采用如下拟合方程拟合电压曲线。
y(t)=A+B·e-αt
其中,y(t)为曲线拟合得到的电压值,A、B及α为待拟合的参数。
S14、采用非线性最小二乘法计算电压曲线中的拟合参数;
S15、根据步骤S14得到的拟合参数计算一阶等效电路模型参数。
一阶等效电路模型参数的计算公式为:
其中,τ为电容电阻网络的时间常数,R1_dis为电容电阻网络对应的放电电阻,R1_ch为电容电阻网络对应的充电电阻,α、B为电压曲线的拟合参数,tdis为电压响应曲线对应电流脉冲放电时刻t2、t1之间的时间段,即tdis=t2-t1,tch为电压响应曲线对应电流脉冲放电时刻t6、t5之间的时间段,即tch=t6-t5。
在通过以上过程计算出退役动力锂电池的一阶等效电路模型的参数后,即可通过查表运算,根据电池当前的充放电倍率、温度及SOC估计值在线更新模型的参数。
S2、采用偏最小二乘法对一阶等效电路模型形成的误差进行二次拟合,建立退役动力锂电池的偏差补偿模型;
在本实施例中,根据一阶等效电路模型的输出与电池端电压测量值之间的偏差,采用偏最小二乘法进行二次拟合,计算得到偏最小二乘偏差补偿模型,建立基于偏最小二乘的退役动力锂电池的偏差补偿模型。
如图4所示,建立基于偏最小二乘的退役动力锂电池的偏差补偿模型的方法为:
S21、设定训练窗口宽度为L,在窗口范围内计算一阶等效电路模型的输出Ut与电池端电压测量值U之间的偏差,获得ΔU=[ΔU1,ΔU2,...,ΔUL]。同时,根据窗口内电池SOC的估计值SOC=[SOC1,SOC2,...,SOCL]及电流测量值It=[It1,It2,...,ItL],建立偏最小二乘的训练样本库,并定义XP=[It,SOC],YP=ΔU;
经过以上步骤,计算得到的偏最小二乘偏差补偿模型为YP=XP·BP,其中,BP=WP·QP,WP=[w1,w2,...,wn],QP=[q1,q2,...,qn]。令由偏最小二乘计算得到的参数向量为BP=[b1,b2,b3],b1、b2、b3为参数向量中包含的元素,则可将基于偏最小二乘的退役动力锂电池偏差补偿模型表示为:
ΔU=b1+b2·SOC+b3·It
其中,SOC为偏最小二乘法训练窗口内电池荷电状态的估计值。
S3、联合一阶等效电路模型和偏差补偿模型,建立退役动力锂电池的混合模型。
在本实施例中,退役动力锂电池的混合模型表示为:
Uhybrid=Ut+ΔU=Uocv-U1-It·R0+b1+b2·SOC+b3·It
其中,Uhybrid表示混合模型的输出电压值。
如图5所示,混合模型可以采用移动加窗的方式更新,定义训练窗口宽度为L,每间隔L个采样周期,根据模型偏差的计算结果,重新训练偏最小二乘偏差补偿模型,提高退役动力电池混合模型的精度。
经过以上所有步骤,即可建立退役动力锂电池的混合模型,并完成模型的在线动态更新。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (5)
1.一种退役动力锂电池的混合建模方法,其特征在于,包括以下步骤:
S1、建立退役动力锂电池的一阶等效电路模型,并获取模型参数;
S2、采用偏最小二乘法对一阶等效电路模型形成的误差进行二次拟合,建立退役动力锂电池的偏差补偿模型,具体为:
根据一阶等效电路模型的输出与电池端电压测量值之间的偏差,采用偏最小二乘法进行二次拟合,计算得到偏最小二乘偏差补偿模型,建立基于偏最小二乘的退役动力锂电池的偏差补偿模型,表示为:
ΔU=b1+b2·SOC+b3·It
其中,ΔU为一阶等效电路模型的输出与电池端电压测量值之间的偏差,SOC为偏最小二乘法训练窗口内电池荷电状态的估计值,It为电池的充放电电流,b1、b2、b3为由偏最小二乘计算的参数;
S3、联合一阶等效电路模型和偏差补偿模型,建立退役动力锂电池的混合模型,表示为:
Uhybrid=Ut+ΔU=Uocv-U1-It·R0+b1+b2·SOC+b3·It
其中,Uhybrid表示混合模型的输出电压值,Ut为电池的端电压,Uocv为电池开路电压,U1为电容电阻网络的端电压,R0为电池内阻。
3.根据权利要求2所述的退役动力锂电池的混合建模方法,其特征在于,所述步骤S1中获取一阶等效电路模型参数的方法为:
S11、对一阶等效电路模型进行离线电流脉冲测试,得到电压响应曲线;
S12、计算一阶等效电路模型在测试条件下的充放电内阻;
S13、分别提取电流脉冲充放电后电池静置时间段的电压测量值,并拟合电压曲线,电压曲线的拟合方程表示为
y(t)=A+B·e-αt
其中,y(t)为电压曲线拟合得到的电压值,A、B及α为电压曲线的拟合参数;
S14、采用非线性最小二乘法计算电压曲线中的拟合参数;
S15、根据步骤S14得到的拟合参数计算一阶等效电路模型参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010827778.5A CN112098846B (zh) | 2020-08-17 | 2020-08-17 | 一种退役动力锂电池的混合建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010827778.5A CN112098846B (zh) | 2020-08-17 | 2020-08-17 | 一种退役动力锂电池的混合建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112098846A CN112098846A (zh) | 2020-12-18 |
CN112098846B true CN112098846B (zh) | 2021-10-22 |
Family
ID=73753844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010827778.5A Active CN112098846B (zh) | 2020-08-17 | 2020-08-17 | 一种退役动力锂电池的混合建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112098846B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104849672A (zh) * | 2015-05-27 | 2015-08-19 | 中国人民解放军国防科学技术大学 | 基于等效电路模型的锂电池动态阻抗参数识别方法 |
EP3224637B1 (fr) * | 2014-11-28 | 2018-10-31 | Renault S.A.S. | Procede automatique d'estimation de la capacite d'une cellule d'une batterie |
CN109633472A (zh) * | 2019-01-23 | 2019-04-16 | 刘平 | 一种锂电池单体的荷电状态估计算法 |
CN110095995A (zh) * | 2019-04-16 | 2019-08-06 | 杭州坤天自动化系统有限公司 | 基于工业大数据和机理混合建模的仿真方法及装置 |
CN110320472A (zh) * | 2019-05-17 | 2019-10-11 | 枣庄学院 | 一种用于矿用锂电池的自修正soc估计方法 |
CN111323705A (zh) * | 2020-03-19 | 2020-06-23 | 山东大学 | 基于鲁棒递归最小二乘的电池参数辨识方法及系统 |
CN111366855A (zh) * | 2020-03-19 | 2020-07-03 | 北京理工大学 | 一种电池等效电路模型抗扰动参数化方法 |
CN111487535A (zh) * | 2020-04-29 | 2020-08-04 | 华中科技大学 | 一种液态金属电池双等效电路模型的参数获取及切换方法 |
-
2020
- 2020-08-17 CN CN202010827778.5A patent/CN112098846B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3224637B1 (fr) * | 2014-11-28 | 2018-10-31 | Renault S.A.S. | Procede automatique d'estimation de la capacite d'une cellule d'une batterie |
CN104849672A (zh) * | 2015-05-27 | 2015-08-19 | 中国人民解放军国防科学技术大学 | 基于等效电路模型的锂电池动态阻抗参数识别方法 |
CN109633472A (zh) * | 2019-01-23 | 2019-04-16 | 刘平 | 一种锂电池单体的荷电状态估计算法 |
CN110095995A (zh) * | 2019-04-16 | 2019-08-06 | 杭州坤天自动化系统有限公司 | 基于工业大数据和机理混合建模的仿真方法及装置 |
CN110320472A (zh) * | 2019-05-17 | 2019-10-11 | 枣庄学院 | 一种用于矿用锂电池的自修正soc估计方法 |
CN111323705A (zh) * | 2020-03-19 | 2020-06-23 | 山东大学 | 基于鲁棒递归最小二乘的电池参数辨识方法及系统 |
CN111366855A (zh) * | 2020-03-19 | 2020-07-03 | 北京理工大学 | 一种电池等效电路模型抗扰动参数化方法 |
CN111487535A (zh) * | 2020-04-29 | 2020-08-04 | 华中科技大学 | 一种液态金属电池双等效电路模型的参数获取及切换方法 |
Non-Patent Citations (1)
Title |
---|
观测模型误差不确定的锂电池SOC估计方法;谈发明 等;《电测与仪表》;20200210;第57卷(第3期);第32-38页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112098846A (zh) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112946499B (zh) | 基于机器学习的锂电池健康状态及荷电状态联合估算方法 | |
CN111581904B (zh) | 考虑循环次数影响的锂电池soc及soh协同估算方法 | |
CN106443473B (zh) | 一种动力锂离子电池组soc估算方法 | |
CN111007400A (zh) | 基于自适应双扩展卡尔曼滤波法的锂电池soc估算方法 | |
CN110261778B (zh) | 一种锂离子电池soc估测算法 | |
CN111537894A (zh) | 一种用于估计锂电池soc和sop的方法 | |
CN110133525B (zh) | 一种应用于电池管理系统的锂离子电池健康状态估计方法 | |
CN110703114B (zh) | 一种基于电-热-神经网络耦合模型的动力电池soc和sot联合状态估计方法 | |
WO2023274194A1 (zh) | 一种适用于富锂锰基电池的高阶模型参数辨识方法和系统 | |
CN109239602B (zh) | 一种动力电池的欧姆内阻的估算方法 | |
CN112557907A (zh) | 一种基于gru-rnn的电动汽车锂离子电池的soc估计方法 | |
CN111142025A (zh) | 一种电池soc估算方法、装置、存储介质及电动汽车 | |
CN110554321B (zh) | 一种实时检测退役动力电池soc的方法 | |
CN110673037B (zh) | 基于改进模拟退火算法的电池soc估算方法及系统 | |
CN113156321A (zh) | 一种锂离子电池荷电状态soc的估算方法 | |
CN112433154A (zh) | 基于ffrls和ekf的锂离子电池soc估测算法 | |
CN112989690A (zh) | 一种混合动力汽车锂电池多时间尺度荷电状态估计方法 | |
CN115389936A (zh) | 一种数模混合驱动的锂电池持续峰值功率能力的在线预测方法 | |
CN109901072B (zh) | 基于历史数据与实验室测试数据的退役电池参数检测方法 | |
CN115327415A (zh) | 基于限定记忆递推最小二乘算法的锂电池soc估算方法 | |
CN110716141A (zh) | 一种基于扩展卡尔曼滤波的电池荷电状态估计方法 | |
CN112098846B (zh) | 一种退役动力锂电池的混合建模方法 | |
WO2021035500A1 (zh) | 用于48v轻混汽车锂离子电池的荷电状态在线估算系统 | |
CN112462281A (zh) | 一种基于气液动力学模型带参数修正的soc估计方法及系统 | |
CN113255215B (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 |