CN109164396B - 一种实时磷酸铁锂电池寿命损耗评估方法 - Google Patents

一种实时磷酸铁锂电池寿命损耗评估方法 Download PDF

Info

Publication number
CN109164396B
CN109164396B CN201811107902.XA CN201811107902A CN109164396B CN 109164396 B CN109164396 B CN 109164396B CN 201811107902 A CN201811107902 A CN 201811107902A CN 109164396 B CN109164396 B CN 109164396B
Authority
CN
China
Prior art keywords
iron phosphate
lithium iron
phosphate battery
soc
total
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
CN201811107902.XA
Other languages
English (en)
Other versions
CN109164396A (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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201811107902.XA priority Critical patent/CN109164396B/zh
Publication of CN109164396A publication Critical patent/CN109164396A/zh
Application granted granted Critical
Publication of CN109164396B publication Critical patent/CN109164396B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Secondary Cells (AREA)

Abstract

本发明公开了属于电力***领域的一种实时磷酸铁锂电池寿命损耗评估方法，该实时评估磷酸铁锂电池寿命损耗包括：进行磷酸铁锂电池循环寿命实验、确定磷酸铁锂电池的AH总量、确定某周期内磷酸铁锂电池的充电量、确定某周期内磷酸铁锂电池的充电量及确定磷酸铁锂电池的寿命损耗等步骤，本磷酸铁锂电池寿命损耗评估方法不需要一次性得到整个循环周期磷酸铁锂电池的充电电流I变化曲线和工作环境温度T变化曲线，在考虑充电速率和环境温度影响的同时，可以根据磷酸铁锂电池每一次充放电过程在线进行磷酸铁锂电池寿命损耗计算。该方法基本结构明确、计算速度快，满足了磷酸铁锂电池寿命损耗评估的实时性要求。

Description

一种实时磷酸铁锂电池寿命损耗评估方法
技术领域
本发明属于电力***领域，特别涉及一种实时磷酸铁锂电池寿命评估方法。
背景技术
由于自然资源的固有特性，以风电和光伏发电为代表的可再生资源发电普遍具有随机波动和不可预测的特性。这两个特性给可再生资源发电并网带来了极大的挑战。储能技术的应用可以有效地缓解可再生资源发电的随机波动和不可预测性，从而实现可再生资源发电的大规模并网，代替部分常规发电机组承担一定基本负荷。
储能技术包括抽水蓄能、压缩空气储能及各种电池储能等。其中，电池储能不受地形、气候等影响，也不需要专门设备，更适合在可再生资源发电中应用。各类电池中，磷酸铁锂电池具有高功率密度、长循环寿命和无记忆等优点，且具有良好的热稳定和化学稳定性，尤其适合于大容量高功率储能场合应用。
在可再生能源发电消纳、智能电网规划设计、优化协调控制中，定量的储能寿命损耗评估是考虑储能寿命成本和进行可靠性分析的关键技术。磷酸铁锂电池寿命的损耗与充放电速率、充放电次数、充放电深度、环境温度、电解液集中度、电解液浓度、高频扩散等因素相关，考虑所有相关因素的寿命预测模型的建立是非常困难的，而且计算速度无法满足电池寿命实时评估的要求。
目前，基于物理化学老化模型的电池寿命评估是最复杂和最完善的方法，但该模型需要大量的数据信息，包括详细的几何数据(电极的厚度和孔数)、电解液的体积和浓度等，这次参数在实验室可以很容易获得，但在电池的正常使用过程中是无法取得的。基于事件导向的老化模型存在事件定义和分类困难问题。本发明基于磷酸铁锂电池的AH总量不随充放电深度而变化的特征，提出了一种实时磷酸铁锂电池寿命损耗评估方法。
发明内容
本发明的目的是提出了一种实时磷酸铁锂电池寿命损耗评估方法，其特征在于，包括如下步骤：
1).进行磷酸铁锂电池循环寿命实验
随机抽取10个同型号的磷酸铁锂电池样本，放电深度DoD分别取10％、20％、30％、40％、50％、60％、70％、80％、90％和100％进行与10次DoD1-DoD10对应的CtF1-CtF10循环寿命10次实验；循环寿命CtF的单位为次，其中，CtF1-CtF10分别为磷酸铁锂电池样本相应放电深度DoD1-DoD10下的测试值，该测试值与磷酸铁锂电池规格和类型密切相关，从而得到磷酸铁锂电池循环寿命与放电深度的对应关系，
2).确定磷酸铁锂电池的AH总量AHT，
根据理论研究和实验表明，不同放电深度下，磷酸铁锂电池的AHT是保持不变的，因此，采用不同放电深度下的AHT平均值作为磷酸铁锂电池的AHT，单位为A·h；计算公式如下：
Figure BDA0001808314210000021
其中，Enom为额定电池容量，单位为A·h；i＝1,2,…,10；DoDi为10％*i的放电深度；CtFi为DoDi对应的循环寿命,单位为次；
磷酸铁锂电池的AH总量只需计算一次，如果规格型号、电池类型等不变，该值一旦确定，可以一直采用，直到磷酸铁锂电池报废。
3).确定某周期内磷酸铁锂电池的充放电量
从磷酸铁锂电池开始投入实际工程应用开始，以24h为循环周期，确定每一周期内磷酸铁锂电池的充放电量，具体内容如下：
(1)每隔1min记录磷酸铁锂电池荷电状态(State of Charge，SoC)，荷电状态SoC的单位为A·h，得到磷酸铁锂电池荷电状态SoC变化曲线，
(2)判断磷酸铁锂电池荷电状态SoC的变化趋势，当磷酸铁锂电池荷电状态SoC逐步增加时，此时磷酸铁锂电池处于充电状态；当磷酸铁锂电池荷电状态SoC逐步减小时，此时磷酸铁锂电池处于放电状态；当相邻时刻磷酸铁锂电池荷电状态SoC大小相等时，此时磷酸铁锂电池处于待机状态，既不充电也不放电；所有充电状态SoC增加量之和为本周期内磷酸铁锂电池的充电量，记为SoC+；所有放电状态SoC减小量之和为本周期内磷酸铁锂电池的放电量SoC-；
4).确定磷酸铁锂电池的总充放电量
从磷酸铁锂电池开始投入实际工程应用开始，将所有循环周期内磷酸铁锂电池的充电量相加，得到磷酸铁锂电池的总充电量，记为SoC+total，计算公式如下：
Figure BDA0001808314210000031
其中，k为循环周期数，k＝1,2,…,K；K为截至计算时的循环周期数；SoC+k为第k次循环周期的磷酸铁锂电池充电量。
从磷酸铁锂电池开始投入实际工程应用开始，将所有循环周期内磷酸铁锂电池的放电量相加，得到磷酸铁锂电池的总放电量，记为SoC-total，计算公式如下：
Figure BDA0001808314210000032
其中，k为循环周期数，k＝1,2,…,K；K为截至计算时的循环周期数；SoC-k为第k次循环周期的磷酸铁锂电池放电量；
5).确定磷酸铁锂电池的寿命损耗LoL，
当磷酸铁锂电池初始运行时的荷电状态SoC与磷酸铁锂电池运行结束时的荷电状态SoC大小相等时，磷酸铁锂电池的总放电量SoC-total等于总充电量SoC+total，磷酸铁锂电池的寿命损耗LoL为：
Figure BDA0001808314210000041
当磷酸铁锂电池初始运行时的荷电状态SoC与磷酸铁锂电池运行结束时的荷电状态SoC大小不相等时，磷酸铁锂电池的总放电量SoC-total与总充电量SoC+total不相等；考虑最恶劣运行状况，磷酸铁锂电池的寿命损耗LoL为：
Figure BDA0001808314210000042
当LoL＝1时，说明磷酸铁锂电池的寿命已经耗尽，此时磷酸铁锂电池需要进行更换。
所述磷酸铁锂电池的AHT只需计算一次，如果磷酸铁锂电池规格型号、电池类型都不变，该值一旦确定，能够一直采用，直到磷酸铁锂电池报废。
本发明的有益效果是提出的磷酸铁锂电池寿命损耗评估方法不需要一次性得到整个循环周期磷酸铁锂电池的荷电状态SoC变化曲线，也无需对磷酸铁锂电池荷电状态SoC时间序列进行预处理，可以根据磷酸铁锂电池荷电状态SoC每一次变化在线进行磷酸铁锂电池寿命损耗计算，满足了磷酸铁锂电池寿命损耗评估的实时性要求
附图说明
图1为磷酸铁锂电池荷电状态24h变化曲线。
图2确定某周期内磷酸铁锂电池充放电量的计算流程。
具体实施方式
本发明提出了一种实时磷酸铁锂电池寿命损耗评估方法，下面结合附图与实施例对本发明进一步说明。
本发明所述对磷酸铁锂电池寿命损耗的实时评估包括如下步骤：
1).进行磷酸铁锂电池循环寿命实验
随机抽取10个同型号的磷酸铁锂电池样本，放电深度DoD分别取10％、20％、30％、40％、50％、60％、70％、80％、90％和100％进行与10次DoD1-DoD10对应的CtF1-CtF10循环寿命10次实验；循环寿命CtF的单位为次，其中，CtF1-CtF10分别为磷酸铁锂电池样本相应放电深度DoD1-DoD10下的测试值，该测试值与磷酸铁锂电池规格和类型密切相关，从而得到磷酸铁锂电池循环寿命与放电深度的对应关系，如表1所示：
表1磷酸铁锂电池循环寿命与放电深度的关系
Figure BDA0001808314210000051
2).确定磷酸铁锂电池的AH总量(Amp Hours Throughput,AHT)
根据理论研究和实验表明，不同放电深度下，磷酸铁锂电池的AHT是保持不变的。因此，本发明采用不同放电深度下AHT平均值作为磷酸铁锂电池的AH T，AH T的单位为A·h，计算公式如下：
Figure BDA0001808314210000052
其中，Enom为额定电池容量，单位为A·h；i＝1,2,…,10；DoDi为10％*i的放电深度；CtFi为DoDi对应的循环寿命,单位为次。
磷酸铁锂电池的AH总量只需计算一次，如果规格型号、电池类型等不变，该值一旦确定，可以一直采用，直到磷酸铁锂电池报废。
3).确定某周期内磷酸铁锂电池的充放电量
从磷酸铁锂电池开始投入实际工程应用开始，以24h为循环周期，确定每一周期内磷酸铁锂电池的充放电量，具体内容如下：
(1)每隔1min记录磷酸铁锂电池荷电状态(State of Charge，SoC)，荷电状态SoC的单位为A·h，得到磷酸铁锂电池荷电状态SoC变化曲线如图1所示)。
(2)判断磷酸铁锂电池荷电状态SoC的变化趋势，当磷酸铁锂电池荷电状态SoC逐步增加时，此时磷酸铁锂电池处于充电状态；当磷酸铁锂电池荷电状态SoC逐步减小时，此时磷酸铁锂电池处于放电状态；当相邻时刻磷酸铁锂电池荷电状态SoC大小相等时，此时磷酸铁锂电池处于待机状态，既不充电也不放电；所有充电状态SoC增加量之和为本周期内磷酸铁锂电池的充电量，记为SoC+；所有放电状态SoC减小量之和为本周期内磷酸铁锂电池的放电量SoC：,其具体计算流程(如图2所示)如下：令SoCs＝SoC0，锂电池的充电量SoC+＝0，锂电池的放电量SoC-＝0；当SoCi>SoCi-1时，SoC+＝SoCi-SoCi-1+SoC+；当SoCi<SoCi-1时，SoC-＝SoCi-1-SoCi+SoC-；当SoCi＝SoCi-1时，不做处理。
4).确定磷酸铁锂电池的总充放电量
从磷酸铁锂电池开始投入实际工程应用开始，将所有循环周期内磷酸铁锂电池的充电量相加，得到磷酸铁锂电池的总充电量，记为SoC+total，计算公式如下：
Figure BDA0001808314210000061
其中，k为循环周期数，k＝1,2,…,K；K为截至计算时的循环周期数；SoC+k为第k次循环周期的磷酸铁锂电池充电量。
从磷酸铁锂电池开始投入实际工程应用开始，将所有循环周期内磷酸铁锂电池的放电量相加，得到磷酸铁锂电池的总放电量，记为SoC-total，计算公式如下：
Figure BDA0001808314210000071
其中，k为循环周期数，k＝1,2,…,K；K为截至计算时的循环周期数；SoC-k为第k次循环周期的磷酸铁锂电池放电量。
5).确定磷酸铁锂电池的寿命损耗(Loss of Lifetime，LoL)
当磷酸铁锂电池初始运行时的荷电状态SoC与磷酸铁锂电池运行结束时的荷电状态SoC大小相等时，磷酸铁锂电池的总放电量SoC-total等于总充电量SoC+total，磷酸铁锂电池的寿命损耗LoL为：
Figure BDA0001808314210000072
当磷酸铁锂电池初始运行时的荷电状态SoC与磷酸铁锂电池运行结束时的荷电状态SoC大小不相等时，磷酸铁锂电池的总放电量SoC-total与总充电量SoC+total不相等，考虑最恶劣运行状况，磷酸铁锂电池的寿命损耗LoL为：
Figure BDA0001808314210000073
当LoL＝1时，说明磷酸铁锂电池的寿命已经耗尽，此时磷酸铁锂电池需要进行更换。
本发明提出的磷酸铁锂电池寿命损耗评估方法不需要一次性得到整个循环周期磷酸铁锂电池的荷电状态SoC变化曲线，也无需对磷酸铁锂电池荷电状态SoC时间序列进行预处理，可以根据磷酸铁锂电池荷电状态SoC每一次变化在线进行磷酸铁锂电池寿命损耗计算，满足了磷酸铁锂电池寿命损耗评估的实时性要求。

Claims (2)

1.一种实时磷酸铁锂电池寿命损耗评估方法，其特征在于，包括如下步骤：
1).进行磷酸铁锂电池循环寿命实验
随机抽取10个同型号的磷酸铁锂电池样本，放电深度DoD分别取10％、20％、30％、40％、50％、60％、70％、80％、90％和100％进行与10次DoD1-DoD10对应的CtF1-CtF10循环寿命10次实验；循环寿命CtF的单位为次，其中，CtF1-CtF10分别为磷酸铁锂电池样本相应放电深度DoD1-DoD10下的测试值，该测试值与磷酸铁锂电池规格和类型密切相关，从而得到磷酸铁锂电池循环寿命与放电深度的对应关系；
2).确定磷酸铁锂电池的AH总量AHT
根据理论研究和实验表明，不同放电深度下，磷酸铁锂电池的AHT是保持不变的，因此，采用不同放电深度下的AHT平均值作为磷酸铁锂电池的AHT，单位为A·h；计算公式如下：
Figure FDA0002618731340000011
其中，Enom为额定电池容量，单位为A·h；i＝1,2,…,10；DoDi为10％*i的放电深度；CtFi为DoDi对应的循环寿命,单位为次；
磷酸铁锂电池的AH总量只需计算一次，如果规格型号、电池类型不变，该值一旦确定，可以一直采用，直到磷酸铁锂电池报废；
3).确定某周期内磷酸铁锂电池的充放电量
从磷酸铁锂电池开始投入实际工程应用开始，以24h为循环周期，确定每一周期内磷酸铁锂电池的充放电量，具体内容如下：
(1)每隔1min记录磷酸铁锂电池荷电状态(State of Charge，SoC)，荷电状态SoC的单位为A·h，得到磷酸铁锂电池荷电状态SoC变化曲线；
(2)判断磷酸铁锂电池荷电状态SoC的变化趋势，当磷酸铁锂电池荷电状态SoC逐步增加时，此时磷酸铁锂电池处于充电状态；当磷酸铁锂电池荷电状态SoC逐步减小时，此时磷酸铁锂电池处于放电状态；当相邻时刻磷酸铁锂电池荷电状态SoC大小相等时，此时磷酸铁锂电池处于待机状态，既不充电也不放电；所有充电状态SoC增加量之和为本周期内磷酸铁锂电池的充电量，记为SoC+；所有放电状态SoC减小量之和为本周期内磷酸铁锂电池的放电量SoC-；
4).确定磷酸铁锂电池的总充放电量
从磷酸铁锂电池开始投入实际工程应用开始，将所有循环周期内磷酸铁锂电池的充电量相加，得到磷酸铁锂电池的总充电量，记为SoC+total，计算公式如下：
Figure FDA0002618731340000021
其中，k为循环周期数，k＝1,2,…,K；K为截至计算时的循环周期数；SoC+k为第k次循环周期的磷酸铁锂电池充电量；
从磷酸铁锂电池开始投入实际工程应用开始，将所有循环周期内磷酸铁锂电池的放电量相加，得到磷酸铁锂电池的总放电量，记为SoC-total，计算公式如下：
Figure FDA0002618731340000022
其中，k为循环周期数，k＝1,2,…,K；K为截至计算时的循环周期数；SoC-k为第k次循环周期的磷酸铁锂电池放电量；
5).确定磷酸铁锂电池的寿命损耗LoL
当磷酸铁锂电池初始运行时的荷电状态SoC与磷酸铁锂电池运行结束时的荷电状态SoC大小相等时，磷酸铁锂电池的总放电量SoC-total等于总充电量SoC+total，磷酸铁锂电池的寿命损耗LoL为：
Figure FDA0002618731340000023
当磷酸铁锂电池初始运行时的荷电状态SoC与磷酸铁锂电池运行结束时的荷电状态SoC大小不相等时，磷酸铁锂电池的总放电量SoC-total与总充电量SoC+total不相等；考虑最恶劣运行状况，磷酸铁锂电池的寿命损耗LoL为：
Figure FDA0002618731340000031
当LoL＝1时，说明磷酸铁锂电池的寿命已经耗尽，此时磷酸铁锂电池需要进行更换。
2.根据权利要求1所述一种实时磷酸铁锂电池寿命损耗评估方法，其特征在于，所述磷酸铁锂电池的AHT只需计算一次，如果磷酸铁锂电池规格型号、电池类型都不变，该值一旦确定，能够一直采用，直到磷酸铁锂电池报废。
CN201811107902.XA 2018-09-21 2018-09-21 一种实时磷酸铁锂电池寿命损耗评估方法 Active CN109164396B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811107902.XA CN109164396B (zh) 2018-09-21 2018-09-21 一种实时磷酸铁锂电池寿命损耗评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811107902.XA CN109164396B (zh) 2018-09-21 2018-09-21 一种实时磷酸铁锂电池寿命损耗评估方法

Publications (2)

Publication Number Publication Date
CN109164396A CN109164396A (zh) 2019-01-08
CN109164396B true CN109164396B (zh) 2020-10-09

Family

ID=64879923

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811107902.XA Active CN109164396B (zh) 2018-09-21 2018-09-21 一种实时磷酸铁锂电池寿命损耗评估方法

Country Status (1)

Country Link
CN (1) CN109164396B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110673036B (zh) * 2019-07-24 2021-09-28 广州中国科学院计算机网络信息中心 基于obd计算技术的车辆故障预警方法及***
CN111707955B (zh) 2020-08-11 2021-01-12 江苏时代新能源科技有限公司 电池剩余寿命的估算方法、装置和介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08185890A (ja) * 1994-12-28 1996-07-16 Tokyo Gas Co Ltd 移動体の駆動源に用いられるバッテリの管理方法
CN1760691A (zh) * 2004-10-12 2006-04-19 三洋电机株式会社 电池的残余容量检测方法及电源装置
CN102520367A (zh) * 2011-12-26 2012-06-27 中国电子科技集团公司第十八研究所 一种空间用氢镍蓄电池寿命评估方法
CN106291366A (zh) * 2015-05-22 2017-01-04 中国电力科学研究院 一种锂离子电池等效循环寿命计算方法
CN107122603A (zh) * 2017-04-24 2017-09-01 国网山东省电力公司潍坊供电公司 一种微电网中蓄电池储能***寿命测算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08185890A (ja) * 1994-12-28 1996-07-16 Tokyo Gas Co Ltd 移動体の駆動源に用いられるバッテリの管理方法
CN1760691A (zh) * 2004-10-12 2006-04-19 三洋电机株式会社 电池的残余容量检测方法及电源装置
CN102520367A (zh) * 2011-12-26 2012-06-27 中国电子科技集团公司第十八研究所 一种空间用氢镍蓄电池寿命评估方法
CN106291366A (zh) * 2015-05-22 2017-01-04 中国电力科学研究院 一种锂离子电池等效循环寿命计算方法
CN107122603A (zh) * 2017-04-24 2017-09-01 国网山东省电力公司潍坊供电公司 一种微电网中蓄电池储能***寿命测算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Capacity loss mechanism of space lithium-ion cells and its life estimation method;Hiroaki YOSHIDA;《ELECTROCHEMISTRY》;20031231;第71卷(第12期);第1018-1024页 *
储能用磷酸铁锂电池循环寿命的能量分析;高飞;《中国电机工程学报》;20130215;第33卷(第5期);第41-45页 *

Also Published As

Publication number Publication date
CN109164396A (zh) 2019-01-08

Similar Documents

Publication Publication Date Title
CN107145628B (zh) 基于电化学-热耦合模型的预测锂电池循环寿命的方法
Stroe et al. Accelerated lifetime testing methodology for lifetime estimation of lithium-ion batteries used in augmented wind power plants
CN110764009B (zh) 一种电动汽车的电池组循环容量及健康状态的估算方法
CN111505503B (zh) 基于微观机理的锂离子电池老化试验方法和试验装置
Zhou et al. Battery behavior prediction and battery working states analysis of a hybrid solar–wind power generation system
CN111505502B (zh) 基于微观机理的时变循环工况下锂离子电池老化试验方法
Svoboda et al. Operating conditions of batteries in off-grid renewable energy systems
Pugach et al. Energy efficiency analysis for a kilo-watt class vanadium redox flow battery system
CN112464571B (zh) 基于多约束条件粒子群优化算法的锂电池组参数辨识方法
CN106597305A (zh) 一种锂离子电池的循环寿命预测方法
CN103412264B (zh) 蓄电池组内单体电池一致性的评价方法
CN109164397B (zh) 考虑充电速率和环境温度的锂电池寿命损耗评估方法
CN112363075A (zh) 一种锂离子电池老化的评估方法
CN108445414B (zh) 一种三元锂离子电池循环寿命的快速测试方法
CN117236264B (zh) 一种预测应力诱导锂离子电池容量衰减的方法
CN109164396B (zh) 一种实时磷酸铁锂电池寿命损耗评估方法
CN112415391A (zh) 基于改进粒子群优化的船用锂离子电池参数辨识方法
CN107066713B (zh) 一种预测锂离子电池材料电化学性能的仿真方法
CN109061478A (zh) 一种利用eis测试进行锂离子电池寿命定性预测的方法
You et al. Charging strategy optimization at low temperatures for li-ion batteries based on multi-factor coupling aging model
CN111999665B (zh) 一种基于微观机理汽车驾驶工况锂离子电池老化试验方法
CN115291131A (zh) 锂离子电池的循环寿命和使用温度预测方法及***
CN109412250B (zh) 一种离网发电***储能电池最佳充电速率的确定方法
CN113447817B (zh) 一种基于两点寿命特征的锂电池在线寿命预测方法
Karoui et al. Diagnosis and prognosis of complex energy storage systems: tools development and feedback on four installed systems

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