CN109444657A - 一种配电网高阻接地故障区段定位方法 - Google Patents

一种配电网高阻接地故障区段定位方法 Download PDF

Info

Publication number
CN109444657A
CN109444657A CN201811234009.3A CN201811234009A CN109444657A CN 109444657 A CN109444657 A CN 109444657A CN 201811234009 A CN201811234009 A CN 201811234009A CN 109444657 A CN109444657 A CN 109444657A
Authority
CN
China
Prior art keywords
periodic component
zero
residual voltage
power distribution
sequence current
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.)
Granted
Application number
CN201811234009.3A
Other languages
English (en)
Other versions
CN109444657B (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.)
China Southern Power Grid Co Ltd
Original Assignee
Research Institute of Southern Power Grid Co Ltd
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 Research Institute of Southern Power Grid Co Ltd filed Critical Research Institute of Southern Power Grid Co Ltd
Priority to CN201811234009.3A priority Critical patent/CN109444657B/zh
Publication of CN109444657A publication Critical patent/CN109444657A/zh
Application granted granted Critical
Publication of CN109444657B publication Critical patent/CN109444657B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Abstract

本发明公开了一种配电网高阻接地故障区段定位方法,包含以下步骤:提取零序电流和零序电压及其衰减周期分量;计算各相邻终端乘积为正的个数差;估算故障点过渡电阻;计算故障定段保护定值;判断是否满足保护判据;判断相位差的最大值。本发明所提故障区段定位方法适应性强,对于不同故障位置、不同过渡电阻,本方法均能准确地找出故障区段,并具有投入成本低,实用性强的优点。

Description

一种配电网高阻接地故障区段定位方法
技术领域
本发明涉及电力系统故障检测领域,特别涉及一种基于零序电流衰减周期分量的配电网高阻接地故障区段定位方法。
背景技术
由于电力需求的日益增长以及土地资源的日渐稀缺,电缆线路在我国大中型城市得到了广泛使用。电缆线路本身存在着较大的对地电容,因此一般采用消弧线圈接地方式对电容电流进行补偿,从而减小单相接地故障发生时的故障电流。
消弧线圈接地(谐振接地)配电网发生高阻接地故障的概率较大,现有的故障检测和定段方法存在准确度低,实现难度大等问题。当其发生高阻接地故障时,如树障或导线坠落在草地,马路,水塘等;由于故障电气量特征不明显,极大增加了故障定位和排除的困难,若系统长时间带故障运行,则可能发展为相间故障,扩大故障范围。
为了减小过渡电阻的影响,提高定位精度的影响,高压电网通常采用双端行波法进行故障定位。然而双端行波故障定位方法对线路两端的同步采样和采样频率要求较高,设备投入较大,且配电网中线路分支众多,因此在配电网中实现双端行波故障定位成本较高。
随着无线通信技术的不断发展和配电网自动化程度的不断提高,配电自动化终端被越来越广泛的应用,配电网也开始利用配电自动化终端进行故障定段。但并未考虑高阻接地的情况,因此,有必要进一步研究针对消弧线圈接地配电网的高阻接地故障检测和定位技术。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种配电网高阻接地故障区段定位方法,此方法可准确地判断故障区段,适应不同故障位置和过渡电阻,具有准确性高、适应性强等优点。
本发明的目的通过以下的技术方案实现:一种配电网高阻接地故障区段定位方法,包含以下步骤:
S1、故障定段终端分别提取零序电流和零序电压及其衰减周期分量;
S2、将各相邻终端零序电流和零序电压的衰减周期分量相乘,并求得乘积为正的个数差;
S3、配电自动化主站利用母线零序电压估算故障过渡电阻,确保所提故障定段方法的灵敏度,得到更精确的定段结果;
S4、计算故障定段保护定值,配电自动化主站利用母线零序电压估算的故障过渡电阻计算故障定段保护定值;
S5、判断是否满足保护判据,配电自动化主站利用各终端零序电流衰减周期分量的相位差差值判断故障区段,只要满足保护判据,则可以判断故障在线路末端,若不满足则可以判断故障在线路某一区段,进行下一步;
S6、判断故障定段终端相位差的最大值的区段,则故障点在该相位差最大的两测量区段之间。
优选的,所述步骤S1中,t0时刻零序电流和零序电压的衰减周期分量,均采用下式计算:
其中,表示t0时刻零序电流或零序电压的衰减周期分量,AQk(t0)表示线路在t0时刻的零序电流值或零序电压值,T表示一个周波时间,单位为ms。
优选的,所述步骤S2中,取两个周波时间的数据窗,计算各终端零序电流衰减周期分量相对于零序电压的相位,将所采集的零序电压和零序电流衰减周期分量的瞬时值相乘,从乘积的正负情况可得出零序电流衰减周期分量和零序电压衰减周期分量的相对相位关系,并将计算的相位上传至配电自动化主站。
更进一步的,令所述各相邻终端乘积为正的个数差为Dpk,采用下式计算:
DPk=|SPk-SP(k+1)|
其中,在一个周波时间T的数据窗内,每个终端有T·n个乘积结果,n为每个终端的采样频率,单位为kHz;设一个周波时间内终端k中乘积为正的个数为SPk(k=1,2,…m),则该终端的零序电流衰减周期分量相对于零序电压的衰减周期分量的相位为360为角度大小。
优选的,所述的故障过渡电阻Rf,采用下式计算:
其中,ω0为角频率,LC为消弧线圈等值电抗,C为系统对地分布电容之和,Uf为故障点虚拟电源的幅值,U0为母线零序电压的幅值。
优选的,步骤S4中,所述的故障保护定值采用下式计算:
DSETO=Krel0·|ks·θmax|
其中,可靠系数Krel0=1.2;采样转换系数ks表示零序电流衰减周期分量相对于零序电压的衰减周期分量的相位与终端采样点数之间的转换关系,即θmax为线路末端发生故障时,不同过渡电阻下各相邻终端相位差之间存在的最大差值。
优选的,步骤S5中,保护的判据为:
max|DPk-DPk+1|≤DSET0
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明随着过渡电阻的增大,线路参数对母线零序电压幅值的影响比重减小,估算误差逐渐减小。定值计算时,过渡电阻误差以每100欧姆为一档来确定,由于计算误差均在可接受的范围内,故其不会对所提定段方法的准确性产生影响。
2、本发明适应性强,对于不同故障位置、不同过渡电阻,所提故障定段方法均能准确地找出故障区段。
3、本发明利用自动化终端获取较高的采样频率,结合零序电流衰减周期分量,降低设备投入,成本较低,实用性较好。
附图说明
图1是本发明一种配电网高阻接地故障区段定位方法流程图。
图2是本发明实施例中消弧线圈接地系统模型图。
具体实施方式
为了更好的理解本发明的技术方案,下面结合附图详细描述本发明提供的实施例,但本发明的实施方式不限于此。
实施例
在RTDS(实时数字仿真仪)中搭建如图2所示的消弧线圈接地系统模型。图2为10.5kV消弧线圈接地配电系统。系统共有三回出线,L1、L2和L3均为电缆线路,长度均为10km,单位长度线路电容为CL=0.37μF,每条线路的负荷均为6MW,功率因数0.95;系统采用过补偿方式,脱谐度为-6%,消弧线圈等值电抗为LC=0.3075H;L1为故障发生线路,n1—n4为不同的故障点(每个故障点均在该区段的中点);L1上等距安装有4个故障定段终端,用以测量零序电流和零序电压。
如图1,测量故障区段的定位方法包含以下顺序的步骤:
S1、提取零序电流和零序电压及其衰减周期分量。
故障定段终端分别提取零序电流和零序电压及其衰减周期分量,所述的t0时刻零序电流和零序电压的衰减周期分量均采用下式计算:
其中,表示t0时刻零序电流或零序电压的衰减周期分量,表示线路在t0时刻的零序电流值或零序电压值,T表示一个周波时间,单位为ms。
S2、将各相邻终端零序电流和零序电压的衰减周期分量相乘,并求得乘积为正的个数差;
取两个周波时间的数据窗,计算各终端零序电流衰减周期分量相对于零序电压的相位,将所采集的零序电压和零序电流衰减周期分量的瞬时值相乘,从乘积的正负情况可得出零序电流衰减周期分量和零序电压衰减周期分量的相对相位关系,并将计算的相位上传至配电自动化主站。令各相邻终端乘积为正的个数差Dpk,采用下式计算:
DPk=|SPk-SP(k+1)|
其中,在一个周波时间即20ms的数据窗内,每个终端有20n个乘积结果,n为每个终端的采样频率,单位为kHz。设一个周波时间内终端k中乘积为正的个数为SPk(k=1,2,…m),则该终端的零序电流衰减周期分量相对于零序电压的衰减周期分量的相位为360为角度大小。
S3、估算故障点过渡电阻。
当发生高阻接地时,各终端测量的与过渡电阻相关,且非故障区段的较低电阻和故障区段的较高电阻可能存在重叠部分。为了确保所提故障定段方法的灵敏度,得到更精确的定段结果,配电自动化主站利用母线零序电压估算故障过渡电阻。所述的故障过渡电阻Rf,采用下式计算:
其中,ω0为角频率,LC为消弧线圈等值电抗,C为系统对地分布电容之和,Uf为故障点虚拟电源的幅值,U0为母线零序电压的幅值。
设线路f2处发生单相接地故障,通过MATLAB计算出故障点的过渡电阻,并与RTDS中实际过渡电阻相比较得到其估算误差,如表1所示。
表1
S4、计算故障定段保护定值。
配电自动化主站利用母线零序电压估算的故障过渡电阻计算故障定段保护定值,具体步骤如下:
DSETO=Krel0·|ks·θmax|
其中,可靠系数Krel0=1.2;采样转换系数ks表示零序电流衰减周期分量相对于零序电压的衰减周期分量的相位与终端采样点数之间的转换关系,即θmax为线路末端发生故障时,不同过渡电阻下各相邻终端相位差之间存在的最大差值。
S5、判断是否满足保护判据。
配电自动化主站利用各终端零序电流衰减周期分量的相位差差值判断故障区段。只要满足保护判据,则可以判断故障在线路末端,若不满足则可以判断故障在线路某一区段,进行下一步。保护的判据为:
max|DPk-DPk+1|≤DSET0
S6、判断故障定段终端相位差的最大值的区段。
对于高阻接地故障发生在线路中某一区段时,基于故障点相邻的两个测量终端测量的相差较大,可依据各个DPk中的最大值max(DPk)来确定故障区段,即故障点在k段与k+1段之间。
基于RTDS对所提的故障定段方法进行验证。假设每个终端的采样频率均为10kHz,取两个周波时间的数据窗进行故障区段判别。表2为不同过渡电阻和不同故障位置下各检测点的SPk、故障定段定值以及故障定段结果。
表2
由表中的计算结果可知,本发明一种配电网高阻接地故障区段定位方法可针对不同故障位置、不同过渡电阻准确地找出故障区段,适应性强、准确率高。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (7)

1.一种配电网高阻接地故障区段定位方法,其特征在于,包含以下步骤:
S1、故障定段终端分别提取零序电流和零序电压及其衰减周期分量;
S2、将各相邻终端零序电流和零序电压的衰减周期分量相乘,并求得乘积为正的个数差;
S3、配电自动化主站利用母线零序电压估算故障过渡电阻;
S4、计算故障定段保护定值,配电自动化主站利用母线零序电压估算的故障过渡电阻计算故障定段保护定值;
S5、判断是否满足保护判据,配电自动化主站利用各终端零序电流衰减周期分量的相位差差值判断故障区段,只要满足保护判据,则可以判断故障在线路末端,若不满足则可以判断故障在线路某一区段,进行下一步;
S6、判断故障定段终端相位差的最大值的区段,则故障点在该相位差最大的两测量区段之间。
2.根据权利要求1所述的配电网高阻接地故障区段定位方法,其特征在于,所述步骤S1中,t0时刻零序电流和零序电压的衰减周期分量均采用下式计算:
其中,表示t0时刻零序电流或零序电压的衰减周期分量,表示线路在t0时刻的零序电流值或零序电压值,T表示一个周波时间,单位为ms。
3.根据权利要求1所述的配电网高阻接地故障区段定位方法,其特征在于,所述步骤S2中,取两个周波时间的数据窗,计算各终端零序电流衰减周期分量相对于零序电压的相位,将所采集的零序电压和零序电流衰减周期分量的瞬时值相乘,从乘积的正负情况可得出零序电流衰减周期分量和零序电压衰减周期分量的相对相位关系,并将计算的相位上传至配电自动化主站。
4.根据权利要求3所述的配电网高阻接地故障区段定位方法,其特征在于,令所述各相邻终端乘积为正的个数差为Dpk,采用下式计算:
DPk=|SPk-SP(k+1)|
其中,在一个周波时间T的数据窗内,每个终端有T·n个乘积结果,n为每个终端的采样频率,单位为kHz;设一个周波时间内终端k中乘积为正的个数为SPk(k=1,2,…m),则该终端的零序电流衰减周期分量相对于零序电压的衰减周期分量的相位为360为角度大小。
5.根据权利要求1所述的配电网高阻接地故障区段定位方法,其特征在于,所述的故障过渡电阻Rf,采用下式计算:
其中,ω0为角频率,LC为消弧线圈等值电抗,C为系统对地分布电容之和,Uf为故障点虚拟电源的幅值,U0为母线零序电压的幅值。
6.根据权利要求1所述的配电网高阻接地故障区段定位方法,其特征在于,步骤S4中,所述的故障保护定值采用下式计算:
DSETO=Krel0·|ks·max|
其中,可靠系数Krel0=1.2;采样转换系数ks表示零序电流衰减周期分量相对于零序电压的衰减周期分量的相位与终端采样点数之间的转换关系,即θmax为线路末端发生故障时,不同过渡电阻下各相邻终端相位差之间存在的最大差值。
7.根据权利要求1所述的配电网高阻接地故障区段定位方法,其特征在于,步骤S5中,保护的判据为:
max|DPk-DPk+1|≤DSET0
CN201811234009.3A 2018-10-23 2018-10-23 一种配电网高阻接地故障区段定位方法 Active CN109444657B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811234009.3A CN109444657B (zh) 2018-10-23 2018-10-23 一种配电网高阻接地故障区段定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811234009.3A CN109444657B (zh) 2018-10-23 2018-10-23 一种配电网高阻接地故障区段定位方法

Publications (2)

Publication Number Publication Date
CN109444657A true CN109444657A (zh) 2019-03-08
CN109444657B CN109444657B (zh) 2021-06-15

Family

ID=65547972

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811234009.3A Active CN109444657B (zh) 2018-10-23 2018-10-23 一种配电网高阻接地故障区段定位方法

Country Status (1)

Country Link
CN (1) CN109444657B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110311380A (zh) * 2019-07-19 2019-10-08 云南电网有限责任公司电力科学研究院 一种有源接地补偿器控制方法
CN110542832A (zh) * 2019-09-18 2019-12-06 南方电网科学研究院有限责任公司 一种配电网高阻接地故障区段定位方法、装置及存储介质
CN111725794A (zh) * 2020-06-03 2020-09-29 南方电网科学研究院有限责任公司 一种消弧线圈并联低电阻接地装置的控制方法及控制装置
CN111896842A (zh) * 2020-07-27 2020-11-06 国网上海市电力公司 基于区间斜率的配电网弧光高阻故障区段定位方法
CN112083283A (zh) * 2020-09-07 2020-12-15 昆明理工大学 一种基于故障分量频带分布差异的配电网故障区段定位方法
CN112782528A (zh) * 2020-12-31 2021-05-11 西安理工大学 一种利用pmu的配电网故障区段定位方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101242097A (zh) * 2008-03-11 2008-08-13 昆明理工大学 利用零模电流测后模拟的小电流接地系统故障选线方法
CN102928731A (zh) * 2012-11-06 2013-02-13 昆明理工大学 一种利用零序电流全量Hough变换的配电网故障选线方法
CN104764969A (zh) * 2015-03-04 2015-07-08 国家电网公司 基于接地电阻实测的双回线路非同名相跨线高阻接地故障定位方法
CN106443349A (zh) * 2016-10-17 2017-02-22 国网福建省电力有限公司 一种高阻接地故障定位方法及系统
CN106980069A (zh) * 2017-05-05 2017-07-25 国网山东省电力公司电力科学研究院 基于暂态电流投影系数差值比较的高阻接地故障定位方法
CN107508265A (zh) * 2017-08-23 2017-12-22 广州供电局有限公司 小电阻接地系统高阻接地保护方法和系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101242097A (zh) * 2008-03-11 2008-08-13 昆明理工大学 利用零模电流测后模拟的小电流接地系统故障选线方法
CN102928731A (zh) * 2012-11-06 2013-02-13 昆明理工大学 一种利用零序电流全量Hough变换的配电网故障选线方法
CN104764969A (zh) * 2015-03-04 2015-07-08 国家电网公司 基于接地电阻实测的双回线路非同名相跨线高阻接地故障定位方法
CN106443349A (zh) * 2016-10-17 2017-02-22 国网福建省电力有限公司 一种高阻接地故障定位方法及系统
CN106980069A (zh) * 2017-05-05 2017-07-25 国网山东省电力公司电力科学研究院 基于暂态电流投影系数差值比较的高阻接地故障定位方法
CN107508265A (zh) * 2017-08-23 2017-12-22 广州供电局有限公司 小电阻接地系统高阻接地保护方法和系统

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110311380A (zh) * 2019-07-19 2019-10-08 云南电网有限责任公司电力科学研究院 一种有源接地补偿器控制方法
CN110311380B (zh) * 2019-07-19 2023-05-23 云南电网有限责任公司电力科学研究院 一种有源接地补偿器控制方法
CN110542832A (zh) * 2019-09-18 2019-12-06 南方电网科学研究院有限责任公司 一种配电网高阻接地故障区段定位方法、装置及存储介质
CN111725794A (zh) * 2020-06-03 2020-09-29 南方电网科学研究院有限责任公司 一种消弧线圈并联低电阻接地装置的控制方法及控制装置
CN111896842A (zh) * 2020-07-27 2020-11-06 国网上海市电力公司 基于区间斜率的配电网弧光高阻故障区段定位方法
CN112083283A (zh) * 2020-09-07 2020-12-15 昆明理工大学 一种基于故障分量频带分布差异的配电网故障区段定位方法
CN112782528A (zh) * 2020-12-31 2021-05-11 西安理工大学 一种利用pmu的配电网故障区段定位方法
CN112782528B (zh) * 2020-12-31 2023-07-18 西安理工大学 一种利用pmu的配电网故障区段定位方法

Also Published As

Publication number Publication date
CN109444657B (zh) 2021-06-15

Similar Documents

Publication Publication Date Title
CN109444657A (zh) 一种配电网高阻接地故障区段定位方法
CN102129009B (zh) 基于双端测量信息的特高压输电线路正序参数测量方法
CN102645613B (zh) 一种基于非接触式磁场测量的输电线路故障定位方法
CN101232176B (zh) 基于暂态行波的中性点非有效接地配电系统故障定位方法
RU2571629C1 (ru) Способ и устройство для оценки угла напряжения нулевой последовательности при однофазном замыкании на землю
CN102135571B (zh) 超高压/特高压多回路输电线路零序阻抗抗干扰测量方法
CN110082636A (zh) 一种电力电缆故障定位方法及系统
CN106771869A (zh) 基于阻性电流分离法的长距离电力电缆绝缘在线监测方法
CN101943737A (zh) 单相接地故障诊断方法和装置
CN105182184B (zh) 一种基于沿线方向行波分解的配网线缆混合直配线路故障测距方法
CN105738759B (zh) 一种基于暂态录波数据的直流输电线路故障定位方法
CN106019080B (zh) 一种基于沿线能量突变的同塔双回直流线路单端行波测距方法
CN103399209A (zh) 一种特高压双极直流输电线路工频参数测量方法
CN103954885A (zh) 基于分布参数的双回线故障单端定位系统及定位方法
CN112540259A (zh) 一种适用于智能配电终端的配网断线故障识别方法和系统
CN108845233A (zh) 配电网架空线路单相接地双端检测定位方法
CN107508265B (zh) 小电阻接地系统高阻接地保护方法和系统
WO2019166903A1 (en) Method and device for fault location in a two-terminal transmission system
CN105842582B (zh) 基于emtr的柔性直流线路故障测距方法
CN103424627B (zh) 双端测量平行电网线路零序阻抗的方法
CN103487725B (zh) 一种基于零序分量法的架空配电线路接地故障指示装置
CN105388355A (zh) 一种运用gps同步的地网分流矢量测试系统及测试方法
CN103884959A (zh) 供电线路故障定位仪及故障定位方法
CN101943738A (zh) 单相接地测距方法和测距装置
CN203825142U (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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20191213

Address after: 510530 Guangdong city of Guangzhou province Luogang District Science City Kexiang Road No. 11

Applicant after: China Southern Power Grid Co., Ltd.

Applicant after: Research Institute of Southern Power Grid Co., Ltd.

Address before: 510623 Guangdong city of Guangzhou province Tianhe District Pearl River Metro Chinese Sui Road No. 6

Applicant before: Research Institute of Southern Power Grid Co., Ltd.

GR01 Patent grant
GR01 Patent grant