CN115308638A - 一种适用于多端直流输电系统的故障线路识别方法 - Google Patents
一种适用于多端直流输电系统的故障线路识别方法 Download PDFInfo
- Publication number
- CN115308638A CN115308638A CN202211050052.0A CN202211050052A CN115308638A CN 115308638 A CN115308638 A CN 115308638A CN 202211050052 A CN202211050052 A CN 202211050052A CN 115308638 A CN115308638 A CN 115308638A
- Authority
- CN
- China
- Prior art keywords
- fault
- line
- traveling wave
- time
- reverse
- 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
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000001174 ascending effect Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000000630 rising effect Effects 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 230000007704 transition Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 6
- 230000004913 activation Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/52—Testing for short-circuits, leakage current or ground faults
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Locating Faults (AREA)
Abstract
本发明涉及一种适用于多端直流输电系统的故障线路识别方法,所述多端直流输电系统的线路L2发生故障后LCC站首端m处反向线模行波呈先下降后上升的趋势,而多端直流输电系统的线路L1故障后反向线模行波呈先快速下降然后趋于恒定的趋势,因此通过检测行波保护启动后线模反向行波的上升速率来进行故障线路的识别。本发明仅利用单端数据信息即可完成故障线路识别,识别速度快,在阈值整定方面不依赖于仿真,可独立完成;且本发明并不依赖于线路边界元件,具有广泛的适用性。
Description
技术领域
本发明属于电力系统及自动化技术领域,具体涉及一种适用于多端直流输电系统的故障线路识别方法。
背景技术
混合多端直流输电系统送端采用电网式换相换流器(LCC),受端采用模块化多电平换流器(MMC),可有效避免逆变侧的换相失败问题。相比其他直流输电方式,混合多端直流输电更是兼具经济性和灵活性的优势,在分布式新能源接入,多落点送电等方面具有广阔的应用前景。
目前,我国已建成昆柳龙混合三端直流输电工程,在支撑国家西电东送战略起到了重要作用。较长的输电距离和复杂的工作环境使得输电线路成为高压直流输电系统中故障率最高的元件。在混合多端直流输电系统中,输电线路发生故障后,LCC换流站强制移相,MMC换流站闭锁或投入负压,以实现故障电流的清除。由于混合多端直流输电系统不止一条直流线路,因此,准确识别并切除故障线路,最小范围地缩小故障隔离区域,是实现非故障区域快速恢复供电的基本前提。
目前已有基于故障后汇流母线两侧电气量特征构造故障线路识别方案的研究成果,但上述方案在识别出故障线路后,还需依靠通信系统将故障线路信息发送到其他换流站,并且在相关阈值整定方面依赖于仿真。
发明内容
本发明的目的在于克服现有技术的不足,提供一种适用于多端直流输电系统的故障线路识别方法,该方法能够实现混合多端直流输电系统中故障线路的快速识别。
本发明解决其技术问题是通过以下技术方案实现的:
一种适用于多端直流输电系统的故障线路识别方法,其特征在于:所述方法的步骤为:
S1、线路L2故障后反向行波由下降转变为上升时刻t0的确定
所述多端直流输电系统的线路L2发生故障后LCC站首端m处反向线模行波呈先下降后上升的趋势,而多端直流输电系统的线路L1故障后反向线模行波呈先快速下降然后趋于恒定的趋势,因此通过检测行波保护启动后线模反向行波的上升速率来进行故障线路的识别;
令线路L2故障后首端m处反向行波导数表达式(1)为0,可以求解得到线路L2故障后理论上反向行波由下降开始转变为上升的时刻t0,t0表达式如式(2)所示:
根据式(1)同时结合实际工程参数特点,可得:
根据式(2)可以将t0表达式简化,结果如式(3)所示。
由式(3)可得t0与故障类型、过渡电阻和故障位置等因素基本无关,根据线路L1参数以及MMC1换流站参数,计算得到线路L2发生故障后,t0约为75μs,因此通过检测保护启动元件启动t0时间后反向行波的上升速率来实现故障线路识别;
S2、故障线路识别判据确定
为去除高频噪声,首先采用巴特沃斯低通滤波器对信号进行处理,在保护元件启动后,开始计算反向行波的差分值,如果当前时刻反向行波差分值大于设定阈值,并且该时刻满足时间范围要求,则判断为线路L2故障,否则判断为线路L1故障;
线路L2故障的识别判据为:
其中:Bm为反向行波值;
t为当前采样时刻;
tst为判断开始时刻,可以设置为保护启动时刻;
tend为判断结束时刻,为在设定时间范围内可以检测到线路L2故障后反向行波上升,要求tend>t0;
Kset为反向行波变化率阈值,可以根据式(4)来进行整定;
t0为理论计算中得到的反向行波开始上升时刻,可由式(3)计算得到;
Krel为可靠系数,取1.2,结合线路L2故障后反向线模行波表达式,可以通过解析计算的方式带入实际运行数据对阈值Kset进行整定;
从反向行波表达式和判据阈值计算式(4)可知,Kset与故障类型和过渡电阻有关,按照单极接地故障且过渡电阻最大的情况进行整定计算;
为说明阈值Kset计算方式的有效性,对线路L2故障后反向行波表达式求二阶导数,可以得到二阶导数为0的时刻t1表达式如式(5)所示:
根据反向行波的二阶导数表达式,可以得到在t0~t1时间段上二阶导数恒大于0,这说明反向行波在t0~t1时间段是凸函数;在t0~t1时间段内对任意两个反向行波采样点做差分运算,其差分值一定大于反向行波极小值点与之后第一个采样点间的差分值,即Kset;据此可以得到判据结束时刻tend应满足tend≤2t0;
在故障点距换流站较近的情况下,故障后第二次折反射行波在判据识别时间范围内可能再次到达首端m处,对故障线路识别造成干扰,因此tend应尽可能接近t0,文中tend取1.2t0;
行波从故障点传输到线路首端需要一定时间,故障点距首端越近,传输时间也越短,通过计算故障行波到达线路首端位置tend时间内反向行波的变化率来实现故障识别,如果在首行波到达线路首端tend时间内,故障点处第二次折反射行波到达,则会对故障识别产生影响,从而形成识别死区;
只有当线路上故障点距换流站的距离lf满足式(6)所示的条件后,该方法才可以区分故障线路,式中v为行波传输速度。
根据式(6)可以计算得到故障距离lf应大于14km,因此该故障线路识别判据在线路L1和L2首端处存在14km的死区。
本发明的优点和有益效果为:
本发明适用于多端直流输电系统的故障线路识别方法,与现有技术相比,本发明仅利用单端数据信息即可完成故障线路识别,识别速度快,在阈值整定方面不依赖于仿真,可独立完成;且本发明并不依赖于线路边界元件,具有广泛的适用性。
附图说明
图1为本发明多端直流输电系统的示意图;
图2为本发明行波传输示意图。
具体实施方式
下面通过具体实施例对本发明作进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
一种适用于多端直流输电系统的故障线路识别方法,其创新之处在于:所述方法的步骤为:
图1中线路L2发生故障后LCC站首端m处反向线模行波呈先下降后上升的趋势,而线路L1故障后反向线模行波呈先快速下降然后趋于恒定的趋势,因此可以通过检测行波保护启动后线模反向行波的上升速率来实现故障线路的识别。
令线路L2故障后首端m处反向行波导数表达式(1)为0,可以求解得到线路L2故障后理论上反向行波由下降开始转变为上升的时刻t0,t0表达式如式(2)所示。
根据式(1)同时结合实际工程参数特点,可得:
根据式(2)可以将t0表达式简化,结果如式(3)所示。
由式(3)可得t0与故障类型、过渡电阻和故障位置等因素基本无关。根据传输线L1参数以及MMC1换流站参数,计算得到线路L2发生故障后,t0约为75μs。因此可以通过检测保护启动元件启动t0时间后反向行波的上升速率来实现故障线路识别。
为去除高频噪声,首先采用巴特沃斯低通滤波器对信号进行处理。在保护元件启动后,开始计算反向行波的差分值。如果当前时刻反向行波差分值大于设定阈值,并且该时刻满足时间范围要求,则判断为线路L2故障,否则判断为线路L1故障。
线路L2故障的识别判据为:
Bm为反向行波值,t为当前采样时刻,tst为判断开始时刻,可以设置为保护启动时刻,tend为判断结束时刻。为在设定时间范围内可以检测到线路L2故障后反向行波上升,要求tend>t0。
Kset为反向行波变化率阈值,可以根据式(4)来进行整定。式中t0为理论计算中得到的反向行波开始上升时刻,可由式(3)计算得到。Bm代表反向行波,Krel为可靠系数,可取1.2。结合线路L2故障后反向线模行波表达式,可以通过解析计算的方式带入实际运行数据对阈值Kset进行整定。
式中t0为理论计算中得到的反向行波开始上升时刻,可由式(3)计算得到。从反向行波表达式和判据阈值计算式(4)可知,Kset与故障类型和过渡电阻有关,因此应按照单极接地故障且过渡电阻最大的情况进行整定计算。
为说明阈值Kset计算方式的有效性,对线路L2故障后反向行波表达式求二阶导数。可以得到二阶导数为0的时刻t1表达式如式(5)所示。
根据反向行波的二阶导数表达式,可以得到在t0~t1时间段上二阶导数恒大于0,这说明反向行波在t0~t1时间段是凸函数。因此理论上在t0~t1时间段内对任意两个反向行波采样点做差分运算,其差分值一定大于反向行波极小值点与之后第一个采样点间的差分值,即Kset。据此可以得到判据结束时刻tend应满足tend≤2t0。在故障点距换流站较近的情况下,故障后第二次折反射行波在判据识别时间范围内可能再次到达首端m处,对故障线路识别造成干扰。因此tend应尽可能接近t0,文中tend取1.2t0。
行波从故障点传输到线路首端需要一定时间,故障点距首端越近,传输时间也越短。本文提出的故障线路识别方法是通过计算故障行波到达线路首端位置tend时间内反向行波的变化率来实现的。如果在首行波到达线路首端tend时间内,故障点处第二次折反射行波到达,则会对故障识别产生影响,从而形成识别死区,如图2所示。
因此只有当线路上故障点距换流站的距离lf满足式(6)所示的条件后,该方法才可以区分故障线路。式中v为行波传输速度。
根据式(6)可以计算得到故障距离lf应大于14km,因此本文提出的故障线路识别判据在线路L1和L2首端处存在14km的死区。
尽管为说明目的公开了本发明的实施例和附图,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内,各种替换、变化和修改都是可能的,因此,本发明的范围不局限于实施例和附图所公开的内容。
Claims (1)
1.一种适用于多端直流输电系统的故障线路识别方法,其特征在于:所述方法的步骤为:
S1、线路L2故障后反向行波由下降转变为上升时刻t0的确定
所述多端直流输电系统的线路L2发生故障后LCC站首端m处反向线模行波呈先下降后上升的趋势,而多端直流输电系统的线路L1故障后反向线模行波呈先快速下降然后趋于恒定的趋势,因此通过检测行波保护启动后线模反向行波的上升速率来进行故障线路的识别;
令线路L2故障后首端m处反向行波导数表达式(1)为0,可以求解得到线路L2故障后理论上反向行波由下降开始转变为上升的时刻t0,t0表达式如式(2)所示:
根据式(1)同时结合实际工程参数特点,可得:
根据式(2)可以将t0表达式简化,结果如式(3)所示。
由式(3)可得t0与故障类型、过渡电阻和故障位置等因素基本无关,根据线路L1参数以及MMC1换流站参数,计算得到线路L2发生故障后,t0约为75μs,因此通过检测保护启动元件启动t0时间后反向行波的上升速率来实现故障线路识别;
S2、故障线路识别判据确定
为去除高频噪声,首先采用巴特沃斯低通滤波器对信号进行处理,在保护元件启动后,开始计算反向行波的差分值,如果当前时刻反向行波差分值大于设定阈值,并且该时刻满足时间范围要求,则判断为线路L2故障,否则判断为线路L1故障;
线路L2故障的识别判据为:
其中:Bm为反向行波值;
t为当前采样时刻;
tst为判断开始时刻,可以设置为保护启动时刻;
tend为判断结束时刻,为在设定时间范围内可以检测到线路L2故障后反向行波上升,要求tend>t0;
Kset为反向行波变化率阈值,可以根据式(4)来进行整定;
t0为理论计算中得到的反向行波开始上升时刻,可由式(3)计算得到;
Krel为可靠系数,取1.2,结合线路L2故障后反向线模行波表达式,可以通过解析计算的方式带入实际运行数据对阈值Kset进行整定;
从反向行波表达式和判据阈值计算式(4)可知,Kset与故障类型和过渡电阻有关,按照单极接地故障且过渡电阻最大的情况进行整定计算;
为说明阈值Kset计算方式的有效性,对线路L2故障后反向行波表达式求二阶导数,可以得到二阶导数为0的时刻t1表达式如式(5)所示:
根据反向行波的二阶导数表达式,可以得到在t0~t1时间段上二阶导数恒大于0,这说明反向行波在t0~t1时间段是凸函数;在t0~t1时间段内对任意两个反向行波采样点做差分运算,其差分值一定大于反向行波极小值点与之后第一个采样点间的差分值,即Kset;据此可以得到判据结束时刻tend应满足tend≤2t0;
在故障点距换流站较近的情况下,故障后第二次折反射行波在判据识别时间范围内可能再次到达首端m处,对故障线路识别造成干扰,因此tend应尽可能接近t0,文中tend取1.2t0;
行波从故障点传输到线路首端需要一定时间,故障点距首端越近,传输时间也越短,通过计算故障行波到达线路首端位置tend时间内反向行波的变化率来实现故障识别,如果在首行波到达线路首端tend时间内,故障点处第二次折反射行波到达,则会对故障识别产生影响,从而形成识别死区;
只有当线路上故障点距换流站的距离lf满足式(6)所示的条件后,该方法才可以区分故障线路,式中v为行波传输速度。
根据式(6)可以计算得到故障距离lf应大于14km,因此该故障线路识别判据在线路L1和L2首端处存在14km的死区。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211050052.0A CN115308638B (zh) | 2022-08-30 | 2022-08-30 | 一种适用于多端直流输电系统的故障线路识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211050052.0A CN115308638B (zh) | 2022-08-30 | 2022-08-30 | 一种适用于多端直流输电系统的故障线路识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115308638A true CN115308638A (zh) | 2022-11-08 |
CN115308638B CN115308638B (zh) | 2024-08-20 |
Family
ID=83865097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211050052.0A Active CN115308638B (zh) | 2022-08-30 | 2022-08-30 | 一种适用于多端直流输电系统的故障线路识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115308638B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117277234A (zh) * | 2023-08-29 | 2023-12-22 | 天津大学 | 适用于并联型多端高压直流输电系统的单端量保护方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170082675A1 (en) * | 2015-09-18 | 2017-03-23 | Schweitzer Engineering Laboratories, Inc. | Time-domain directional line protection of electric power delivery systems |
CN106771700A (zh) * | 2016-11-11 | 2017-05-31 | 清华大学 | 柔性直流输电线路雷击干扰的快速识别方法及装置 |
CN108767820A (zh) * | 2018-06-09 | 2018-11-06 | 天津大学 | 一种适用于多端直流输电系统的新型方向保护方法 |
CN110391645A (zh) * | 2019-06-25 | 2019-10-29 | 北京四方继保自动化股份有限公司 | 基于故障行波突变量的直流配电网故障选线方法及系统 |
CN111257693A (zh) * | 2020-02-25 | 2020-06-09 | 杭州柯林电气股份有限公司 | 小电流接地系统单相接地故障定位方法、装置 |
CN111342435A (zh) * | 2020-03-26 | 2020-06-26 | 四川大学 | 适用于柔性直流输电系统输电线路的单端量行波保护方法 |
CN111812453A (zh) * | 2020-06-23 | 2020-10-23 | 华南理工大学 | 高压直流输电线路单端量保护方法、系统及存储介质 |
CN112363013A (zh) * | 2020-10-30 | 2021-02-12 | 西安理工大学 | 一种混合直流输电系统直流线路区内外故障判别方法 |
CN112886547A (zh) * | 2020-12-31 | 2021-06-01 | 山东大学 | 基于行波波形距离的高压输电线路行波保护方法及系统 |
CN113702757A (zh) * | 2021-08-13 | 2021-11-26 | 天津大学 | 基于电流/电压高频行波比的直流配网快速方向保护方法 |
CN114089104A (zh) * | 2021-11-12 | 2022-02-25 | 西安交通大学 | 一种基于线路边界线路侧电压的故障方向判别方法及系统 |
CN114629089A (zh) * | 2022-03-14 | 2022-06-14 | 四川大学 | 适用于柔性直流输电线路的单端量波形相似度保护方法 |
-
2022
- 2022-08-30 CN CN202211050052.0A patent/CN115308638B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170082675A1 (en) * | 2015-09-18 | 2017-03-23 | Schweitzer Engineering Laboratories, Inc. | Time-domain directional line protection of electric power delivery systems |
CN106771700A (zh) * | 2016-11-11 | 2017-05-31 | 清华大学 | 柔性直流输电线路雷击干扰的快速识别方法及装置 |
CN108767820A (zh) * | 2018-06-09 | 2018-11-06 | 天津大学 | 一种适用于多端直流输电系统的新型方向保护方法 |
CN110391645A (zh) * | 2019-06-25 | 2019-10-29 | 北京四方继保自动化股份有限公司 | 基于故障行波突变量的直流配电网故障选线方法及系统 |
CN111257693A (zh) * | 2020-02-25 | 2020-06-09 | 杭州柯林电气股份有限公司 | 小电流接地系统单相接地故障定位方法、装置 |
CN111342435A (zh) * | 2020-03-26 | 2020-06-26 | 四川大学 | 适用于柔性直流输电系统输电线路的单端量行波保护方法 |
CN111812453A (zh) * | 2020-06-23 | 2020-10-23 | 华南理工大学 | 高压直流输电线路单端量保护方法、系统及存储介质 |
CN112363013A (zh) * | 2020-10-30 | 2021-02-12 | 西安理工大学 | 一种混合直流输电系统直流线路区内外故障判别方法 |
CN112886547A (zh) * | 2020-12-31 | 2021-06-01 | 山东大学 | 基于行波波形距离的高压输电线路行波保护方法及系统 |
CN113702757A (zh) * | 2021-08-13 | 2021-11-26 | 天津大学 | 基于电流/电压高频行波比的直流配网快速方向保护方法 |
CN114089104A (zh) * | 2021-11-12 | 2022-02-25 | 西安交通大学 | 一种基于线路边界线路侧电压的故障方向判别方法及系统 |
CN114629089A (zh) * | 2022-03-14 | 2022-06-14 | 四川大学 | 适用于柔性直流输电线路的单端量波形相似度保护方法 |
Non-Patent Citations (1)
Title |
---|
白嘉;徐玉琴;王增平;: "基于形态学小波综合算法的超高压输电线路单端行波保护新原理", 电网技术, no. 02, 30 January 2006 (2006-01-30) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117277234A (zh) * | 2023-08-29 | 2023-12-22 | 天津大学 | 适用于并联型多端高压直流输电系统的单端量保护方法 |
CN117277234B (zh) * | 2023-08-29 | 2024-06-14 | 天津大学 | 适用于并联型多端高压直流输电系统的单端量保护方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115308638B (zh) | 2024-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107359602B (zh) | 一种适用于输电线路的自同步纵联电流差动保护方法 | |
CN103176107B (zh) | 高压直流输电线路混合故障测距方法 | |
CN107422196B (zh) | 一种基于谐波畸变率和频率的孤岛检测方法 | |
CN115308638A (zh) | 一种适用于多端直流输电系统的故障线路识别方法 | |
CN103412240B (zh) | 一种不依赖波头辨识的同塔双回输电线路单端行波故障测距方法 | |
CN110932264B (zh) | 基于台区末端智能感知的拓扑识别方法 | |
CN108599114A (zh) | 一种高压交直流混联系统交流线路暂态方向保护方法 | |
CN107390084B (zh) | 故障方向检测方法、装置、继电器和计算机可读存储介质 | |
CN112886547B (zh) | 基于行波波形距离的高压输电线路行波保护方法及系统 | |
CN104820168A (zh) | 一种基于波形差异度和雷击故障样本库的雷击故障判别方法 | |
CN108199356B (zh) | 基于波前信息的直流输电线路超高速保护方法 | |
CN113655339B (zh) | 一种直流输电线路保护系统故障定位方法及装置 | |
CN112083280B (zh) | 一种识别混合多端直流输电系统故障区间的方法 | |
CN102313858B (zh) | 一种初始反极性方向行波的辨识方法 | |
CN112098779A (zh) | 一种架空线-电缆混合线路故障点定位方法和系统 | |
CN116345418A (zh) | 一种基于标准反行波面积的柔性直流输电线路距离保护方法 | |
CN109613398B (zh) | 一种用于高压交直流混联电网的故障方向判别方法、判别元件 | |
CN114142442B (zh) | 一种直流输电线路保护方法及系统 | |
CN113702757B (zh) | 基于电流/电压高频行波比的直流配网快速方向保护方法 | |
CN114465215A (zh) | 一种海上风电直流送出线路纵联保护方法及系统 | |
CN113644635B (zh) | 一种柔性直流输电系统直流线路区内外故障判别方法 | |
CN110943447B (zh) | 一种电力系统失步振荡中心识别方法及装置 | |
CN111766477B (zh) | 直流无断路器配电网线路故障快速检测与识别方法及装置 | |
CN114629089A (zh) | 适用于柔性直流输电线路的单端量波形相似度保护方法 | |
CN113820563A (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 |