CN110007198B - 一种单相接地故障启动方法 - Google Patents
一种单相接地故障启动方法 Download PDFInfo
- Publication number
- CN110007198B CN110007198B CN201910404624.2A CN201910404624A CN110007198B CN 110007198 B CN110007198 B CN 110007198B CN 201910404624 A CN201910404624 A CN 201910404624A CN 110007198 B CN110007198 B CN 110007198B
- Authority
- CN
- China
- Prior art keywords
- adaptive threshold
- frequency
- component
- low
- frequency component
- 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/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- 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/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Emergency Protection Circuit Devices (AREA)
Abstract
本发明涉及一种新型的单相接地故障启动方法。包括对母线零序电压信号进行小波包变换获得特定频带分量,以及通过对高、低频分量的历史数据进行进一步处理,得到用于诊断系统是否发生扰动的自适应阈值,根据系统是否发生扰动以及扰动持续的时间,实现对单相接地故障的快速检测。本发明能够快速、准确地识别单相接地故障,此外,本发明可以有效区分单相接地故障和投切带载线路、投切电容器组等非故障行为;本发明具有较强的适应能力,在噪声干扰的工况下仍具有较高的识别正确率。
Description
技术领域
本发明涉及一种配电网单相接地故障启动方法。
背景技术
中低压配电网通常指35kV及以下电压等级的电力网络,其分布范围广,结构复杂,与用户联系紧密。我国6~35kV配电网普遍采用小电流接地方式,在小电流接地方式下单相接地故障发生率最高,占配电网故障的80%以上。随着电网规模的扩大,特别是电缆的使用量增多,线路距离增长,使得电网对地电容增大,一旦发生接地故障,由此产生的电弧很难熄灭,将对电力系统的正常运行造成严重影响。因此,研究单相接地故障的快速识别方法,有利于快速定位和排除故障,对保证电网的安全运行具有实际应用价值。
本申请提出一种基于离散小波变换和自适应阈值的单相接地故障启动算法,先对母线零序电压进行小波包变换得到特定频带分量,再利用所得到的高、低频分量构造自适应阈值,将当前时刻的高、低频分量变化量与由历史数据计算得到的自适应阈值进行比较,若变化量高于自适应阈值,则认为系统发生了扰动。若扰动持续时间超过预设阈值,则认为配电网中发生了单相接地故障。
发明内容
本发明的目的在于提供一种配电网单相接地故障启动方法,该方法能够快速、准确地识别单相接地故障,此外,可以有效区分单相接地故障和投切带载线路、投切电容器组等非故障行为;本发明方法具有较强的适应能力,在噪声干扰的工况下仍具有较高的识别正确率。
为实现上述目的,本发明的技术方案是:一种配电网单相接地故障启动方法,包括如下步骤:
步骤S1、对实时采集的零序电压信号进行Mallat小波包分解,获取零序电压信号的频带分量,在此基础上,计算选定的高、低频分量的周期变化量;
步骤S2、根据选定的高、低频分量及其周期变化量的历史数据,分别构造高频自适应阈值以及低频自适应阈值;
步骤S3、比较高、低频分量周期变化量与其自适应阈值之间的大小关系,若高、低频分量周期变化量中任意一者超过其自适应阈值,则认为系统发生了扰动;
步骤S4、监视扰动持续时间,根据持续时间判断是否发生了单相接地故障。
在本发明一实施例中,所述步骤S1的具体过程如下:
根据实时采集的零序电压信号,对其进行Mallat小波包分解:
(p为偶数,j>0)
(p为奇数,j>0)
其中,Sj和Wj为第j尺度的近似分量和细节分量,hk和gk分别为低通滤波器系数和高通滤波器系数;
从分解得到的各分量中选择由第三尺度细节分量进一步分解得到的细节分量(dd4)和由第四尺度近似分量(cc4)作为特征信号,并进一步计算它们的周期变化量,计算式如式(3)所示:
其中,n代表当前点与参考点之间相隔的周期数,L代表每个工频周期的采样点数。
在本发明一实施例中,所述步骤S2的具体过程如下:
自适应阈值由两种分量组成,其中一个根据高、低频分量周期变化量计算得到,计算式如式(4)所示:
其中,Nl为每一片历史数据片段的长度,Np为分隔区的长度;
另一个分量根据当前周期之前数个周期的小波系数的最大值计算得到,其计算式如式(5) 所示:
其中,M为计算一个阈值所需的历史片段的个数;
将上述得到的数值乘以相应的系数后再叠加,即可得到算法所需的自适应阈值,其计算式如式(6)所示:
其中,Kε_d和Kε_c为高低频分量的极大值的利用系数,Kε_c为可靠度;
计算得到的λs_d和λs_c两种自适应阈值分别用于高频变化量Δd和低频变化量Δc的监测。
在本发明一实施例中,所述步骤S3的具体实现方式为:计算得到的λs_d和λs_c两种自适应阈值分别用于高频变化量Δd和低频变化量Δc的监测,当高频变化量Δd和低频变化量Δc中任意一种变化量超过了所对应的自适应阈值,则认为出现了扰动,所对应的时刻被记录下来。
在本发明一实施例中,所述步骤S4的具体实现方式为:当算法认为出现扰动时,自适应阈值的计算工作将被暂停,直到判断流程结束;在等待t1秒后,算法将再次判断变化量是否还高于自适应阈值;若在后续的t2秒内,变化量超过自适应阈值的次数超过δ,算法将输出启动信号以及故障时刻;反之,所记录的故障时刻将被清除,算法回到正常监测状态。
相较于现有技术,本发明具有以下有益效果:
1、本发明利用Mallat小波包变换算法对采集到的零序电压信号进行实时分解计算,可在仅仅已知低通和高通滤波器系数的情况下实现对离散信号的小波分解。该算法既保留了小波变换强大的时频分解能力,也简化了小波变化的复杂的计算函数,可在极大程度上降低启动算法的计算负担和应用难度;
2、本发明利用由历史数据构造的自适应阈值来检测单相接地故障,可提高算法的灵敏性,并可降低噪声、采样率、电力网络结构变化等因素对算法正确率的影响;
3、本发明将扰动持续时间作为启动算法的判据之一,可以极大提高算法的抗干扰能力。
附图说明
图1为本发明的流程图。
图2为本发明实施例中所应用的10kV配电网模型。
图3为本发明的算法效果示意图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本实施例提供一种配电网单相接地故障启动方法,如图1所示,包括如下步骤:
步骤S1:对实时采集的零序电压信号进行Mallat小波包分解,获取零序电压信号的特定频带分量,在此基础上,计算选定的高、低频分量的周期变化量;本步骤具体包括如下步骤:
根据实时采集的零序电压数据,对其进行Mallat小波包分解,Mallat算法的递推公式如式 (1)和式(2)所示:
(p为偶数,j>0)
(p为奇数,j>0)
其中,Sj和Wj为第j尺度的近似分量和细节分量,hk和gk分别为低通滤波器系数和高通滤波器系数;
从分解得到的各分量中选择由第三尺度细节分量进一步分解得到的细节分量(dd4)和由第四尺度近似分量(cc4)作为特征信号,并进一步计算它们的周期变化量,计算式如式(3)所示:
其中,n代表当前点与参考点之间相隔的周期数,L代表每个工频周期的采样点数。
步骤S2:根据选定的高、低频分量及其周期变化量的历史数据,分别构造高频自适应阈值以及低频自适应阈值;本步骤具体包括如下步骤:
自适应阈值由两种分量组成,其中一个根据高、低频分量周期变化量计算得到,计算式如式(4)所示:
其中,Nl为每一片历史数据片段的长度,Np为分隔区的长度;
另一个分量根据当前周期之前数个周期的小波系数的最大值计算得到,其计算式如式(5) 所示:
其中,M为计算一个阈值所需的历史片段的个数;
将上述得到的数值乘以相应的系数后再叠加,即可得到算法所需的自适应阈值,其计算式如式(6)所示:
其中,Kε_d和Kε_c为高低频分量的极大值的利用系数,Kε_c为可靠度;
步骤S3:比较高、低频分量变化量与其自适应阈值之间的大小关系,若其中任意一者超过其自适应阈值,则认为系统发生了扰动;本步骤具体包括如下步骤:
计算得到的λs_d和λs_c两种自适应阈值分别用于高频变化量Δd和低频变化量Δc的监测,当其中任意一种变化量超过了所对应的自适应阈值,则认为出现了扰动,所对应的时刻被记录下来。
步骤S4:监视扰动持续时间,根据持续时间判断是否发生了单相接地故障。本步骤具体包括如下步骤:
当算法认为出现了扰动时,自适应阈值的计算工作将被暂停,直到判断流程结束。在等待 t1秒后,算法将再次判断变化量是否还高于自适应阈值。若在后续的t2秒内,变化量超过自适应阈值的次数超过δ,算法将输出启动信号以及故障时刻;反之,所记录的故障时刻将被清除,算法回到正常监测状态。
在本实施例中,如图2所示,利用PSCAD/EMTDC仿真软件搭建10kV配电网模型用于获取扰动数据,测试结果表明,该方法对不同的故障点、初相角、接地电阻所发生的配电网单相接地故障均能快速准确识别,且在噪声干扰下有较好的适应性,在此基础上进行多种类型的单相接地故障、投切带载线路行为和投切电容器组行为的模拟实验,并采集零序电压波形。仿真模型中,110kV高压线路用三相电源代替,仿真元件主要包括:电压互感器、系统电源、变压器、线路等。110kV/10kV变压器连接组别为YNd11,一次侧与二次侧的电阻标幺值为0.0019,电感标幺值为0.75,励磁铁芯电阻标幺值为1615.12,励磁铁芯电感标幺值为833.23;10kV/0.4kV 变压器连接组别为Dyn11,一次侧与二次侧的电阻标幺值为0.00501,电感标幺值为0.0223,励磁铁芯电阻标幺值为869.27,励磁铁芯电感标幺值为142.35;电磁式PT的励磁参数为:电压标幺值为1、1.328、1.501、1.79、1.963,对应的电流标幺值为1、1.733、3.067、7.33、11.93;10kV 线路模块选用三相π型等值电路模块,架空线路的正序电阻0.17Ω/km,正序电容0.0097μF/km,正序电感1.21mH/km,零序电阻0.23Ω/km,零序电容0.008μF/km,零序电感5.478mH/km;电缆线路的正序电阻0.27Ω/km,正序电容0.339μF/km,正序电感0.255mH/km,零序电阻2.7Ω/km,零序电容0.28μF/km,零序电感1.019mH/km。
在仿真样本数据中加入30dB的高斯白噪声,再利用所提出的启动算法对仿真数据进行检测,验证本启动方法的性能。图3展示了所提出的启动方法针对单相接地故障、投切带载线路行为和投切电容器组行为的算法效果示意图,其中图(a)从左至右四张小图分别为单相接地故障时的零序电压波形、高频变化量Δd及其自适应阈值λs_d波形、低频变化量Δc及其自适应阈值λs_c波形、启动信号波形,从图中可以看出,当故障发生后,高低频变化量可迅速越过自适应阈值且衰减速度较慢,算法可以准确识别故障,并输出启动信号。从图3的(b)和(c)两张子图可以看出,高低频变化量基本不会超过自适应阈值,算法不会对此类扰动做出反应。针对所有仿真数据的识别结果如表1所示,识别准确率高达100%,说明了所提方法具有良好的可靠性和适应性。
表1启动算法识别效果
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (3)
1.一种单相接地故障启动方法,其特征在于,包括如下步骤:
步骤S1、对实时采集的零序电压信号进行Mallat小波包分解,获取零序电压信号的频带分量,在此基础上,计算选定的高、低频分量的周期变化量;具体实现如下:
根据实时采集的零序电压信号,对其进行Mallat小波包分解:
p为偶数,j>0;
p为奇数,j>0 ;
其中,Sj和Wj为第j尺度的近似分量和细节分量,hk和gk分别为低通滤波器系数和高通滤波器系数;
从分解得到的各分量中选择由第三尺度细节分量进一步分解得到的细节分量(dd4)和由第四尺度近似分量(cc4)作为特征信号,并进一步计算它们的周期变化量,计算式如式(3)所示:
其中,n代表当前点与参考点之间相隔的周期数,L代表每个工频周期的采样点数;
步骤S2、根据选定的高、低频分量及其周期变化量的历史数据,分别构造高频自适应阈值以及低频自适应阈值;具体实现如下:
自适应阈值由两种分量组成,其中一个根据高、低频分量周期变化量计算得到,计算式如式(4)所示:
其中,Nl为每一片历史数据片段的长度,Np为分隔区的长度,Δd、Δc分别为高、低频分量周期变化量;
另一个分量根据当前周期之前数个周期的小波系数的最大值计算得到,其计算式如式(5)所示:
其中,M为计算一个阈值所需的历史片段的个数;
将上述得到的数值乘以相应的系数后再叠加,即可得到算法所需的自适应阈值,其计算式如式(6)所示:
其中,Kε_d和Kε_c为高低频分量的极大值的利用系数;
计算得到的λs_d和λs_c两种自适应阈值分别用于高频变化量Δd和低频变化量Δc的监测;
步骤S3、比较高、低频分量周期变化量与其自适应阈值之间的大小关系,若高、低频分量周期变化量中任意一者超过其自适应阈值,则认为系统发生了扰动;
步骤S4、监视扰动持续时间,根据持续时间判断是否发生了单相接地故障。
2.根据权利要求1所述的一种单相接地故障启动方法,其特征在于,所述步骤S3的具体实现方式为:计算得到的λs_d和λs_c两种自适应阈值分别用于高频变化量Δd和低频变化量Δc的监测,当高频变化量Δd和低频变化量Δc中任意一种变化量超过了所对应的自适应阈值,则认为出现了扰动,所对应的时刻被记录下来。
3.根据权利要求1或2所述的一种单相接地故障启动方法,其特征在于,所述步骤S4的具体实现方式为:当算法认为出现扰动时,自适应阈值的计算工作将被暂停,直到判断流程结束;在等待t1秒后,算法将再次判断变化量是否还高于自适应阈值;若在后续的t2秒内,变化量超过自适应阈值的次数超过δ,算法将输出启动信号以及故障时刻;反之,所记录的故障时刻将被清除,算法回到正常监测状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910404624.2A CN110007198B (zh) | 2019-05-16 | 2019-05-16 | 一种单相接地故障启动方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910404624.2A CN110007198B (zh) | 2019-05-16 | 2019-05-16 | 一种单相接地故障启动方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110007198A CN110007198A (zh) | 2019-07-12 |
CN110007198B true CN110007198B (zh) | 2021-03-02 |
Family
ID=67177054
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910404624.2A Active CN110007198B (zh) | 2019-05-16 | 2019-05-16 | 一种单相接地故障启动方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110007198B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111487507A (zh) * | 2020-06-01 | 2020-08-04 | 南方电网调峰调频发电有限公司西部检修试验分公司 | 基于小波包能量比值的配电网高阻故障辨识方法 |
CN112067961B (zh) * | 2020-10-13 | 2023-08-15 | 哈尔滨工业大学(深圳) | 电弧故障检测方法、系统及存储介质 |
CN112684282B (zh) * | 2020-11-12 | 2022-07-19 | 国网河北省电力有限公司电力科学研究院 | 配电网单相接地故障识别方法、装置及终端设备 |
CN113452002B (zh) * | 2021-08-30 | 2021-12-31 | 国网江西省电力有限公司电力科学研究院 | 基于垂段最小距离拟合的高压送出线距离保护方法及装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858948A (zh) * | 2009-04-10 | 2010-10-13 | 阿海珐输配电英国有限公司 | 用于在三相中压配电系统中进行暂态和间歇性接地故障检测和方向确定的方法和系统 |
CN101865967A (zh) * | 2010-06-18 | 2010-10-20 | 广州智光电气股份有限公司 | 单相接地故障选线的装置及其方法 |
CN101943737A (zh) * | 2010-08-04 | 2011-01-12 | 清华大学 | 单相接地故障诊断方法和装置 |
CN103323728A (zh) * | 2013-06-05 | 2013-09-25 | 昆明理工大学 | 基于全周波能量比的单相接地故障与虚幻接地识别方法 |
RU2013142540A (ru) * | 2013-09-17 | 2015-03-27 | Общество с ограниченной ответственностью "Исследовательский центр "Бреслер" | Способ определения места однофазного замыкания фидера на землю |
CN104991162A (zh) * | 2015-06-21 | 2015-10-21 | 云南电力试验研究院(集团)有限公司 | 一种小电流接地系统单相接地故障的定位装置及使用方法 |
CN106771870A (zh) * | 2016-12-26 | 2017-05-31 | 北京国电通网络技术有限公司 | 一种配电网接地故障定位方法及系统 |
CN108279364A (zh) * | 2018-01-30 | 2018-07-13 | 福州大学 | 基于卷积神经网络的配电网单相接地故障选线方法 |
CN108802557A (zh) * | 2018-06-28 | 2018-11-13 | 安徽康能电气有限公司 | 一种基于扰动信号的故障定位装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508076B (zh) * | 2011-11-09 | 2013-09-25 | 东北大学 | 一种基于多智能体系统和小波分析的故障诊断装置及方法 |
-
2019
- 2019-05-16 CN CN201910404624.2A patent/CN110007198B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858948A (zh) * | 2009-04-10 | 2010-10-13 | 阿海珐输配电英国有限公司 | 用于在三相中压配电系统中进行暂态和间歇性接地故障检测和方向确定的方法和系统 |
CN101865967A (zh) * | 2010-06-18 | 2010-10-20 | 广州智光电气股份有限公司 | 单相接地故障选线的装置及其方法 |
CN101943737A (zh) * | 2010-08-04 | 2011-01-12 | 清华大学 | 单相接地故障诊断方法和装置 |
CN103323728A (zh) * | 2013-06-05 | 2013-09-25 | 昆明理工大学 | 基于全周波能量比的单相接地故障与虚幻接地识别方法 |
RU2013142540A (ru) * | 2013-09-17 | 2015-03-27 | Общество с ограниченной ответственностью "Исследовательский центр "Бреслер" | Способ определения места однофазного замыкания фидера на землю |
CN104991162A (zh) * | 2015-06-21 | 2015-10-21 | 云南电力试验研究院(集团)有限公司 | 一种小电流接地系统单相接地故障的定位装置及使用方法 |
CN106771870A (zh) * | 2016-12-26 | 2017-05-31 | 北京国电通网络技术有限公司 | 一种配电网接地故障定位方法及系统 |
CN108279364A (zh) * | 2018-01-30 | 2018-07-13 | 福州大学 | 基于卷积神经网络的配电网单相接地故障选线方法 |
CN108802557A (zh) * | 2018-06-28 | 2018-11-13 | 安徽康能电气有限公司 | 一种基于扰动信号的故障定位装置 |
Non-Patent Citations (2)
Title |
---|
Features-clustering-based earth fault detection using singular-value decomposition and fuzzy c-means in resonant grounding distribution systems;Mou-Fa Guo, Nien-Che Yang;《Electrical Power and Energy Systems》;20170529;第97-108页 * |
基于HHT及信号注入的配电网谐振与单相接地故障识别;蔡晋等;《电气技术》;20151215;第31-35、51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110007198A (zh) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110007198B (zh) | 一种单相接地故障启动方法 | |
CN103018627B (zh) | 非有效接地系统故障类型自适应接地选线方法 | |
CN103675605B (zh) | 一种基于故障信号暂态相关分析的配电网小电流接地故障选线方法 | |
CN103344875B (zh) | 谐振接地系统单相接地故障分类选线方法 | |
CN102510044A (zh) | 基于小波变换和概率神经网络的励磁涌流鉴别方法 | |
CN111308272B (zh) | 一种小电流接地故障区段定位方法 | |
CN111308263B (zh) | 一种配电网高阻接地故障检测方法 | |
CN105738764A (zh) | 基于暂态信息全频带的配电网故障区段定位方法 | |
CN111007364B (zh) | 一种电缆早期自恢复故障的识别方法 | |
CN107329044A (zh) | 一种基于电弧暂态分量的配电网单相接地故障选线方法 | |
CN111381129B (zh) | 基于超低频信号的接地故障线路及类型识别方法和装置 | |
CN112014773B (zh) | 一种用于检测小电流接地系统电缆早期故障的方法 | |
CN102435896B (zh) | 船舶中压电力系统中间歇性接地故障的快速识别方法 | |
CN115128400A (zh) | 一种配网故障类型识别与故障选线综合研判方法及系统 | |
CN116298665A (zh) | 一种配电电缆弧光接地故障判定方法及系统 | |
CN117741505A (zh) | 基于峰度和Teager能量算子的配电网相继故障检测方法及系统 | |
CN113376477B (zh) | 基于行波能谱矩阵相似度的柔性直流电网单端保护方法 | |
CN112782528B (zh) | 一种利用pmu的配电网故障区段定位方法 | |
CN112748362B (zh) | 基于vmd和灰色关联度相结合的小电流接地故障检测方法 | |
CN109324268B (zh) | 基于贝叶斯推理的配电网早期故障检测方法及装置 | |
CN114910744B (zh) | 基于s变换和自适应平均奇异熵的高阻接地故障检测方法 | |
CN116660677A (zh) | 一种基于vmd分解的配电网单相接地故障的多判据选线方法 | |
CN113552441A (zh) | 一种单相接地故障检测方法及装置 | |
CN116609610A (zh) | 一种基于自适应小波包变换的单相接地故障选相方法 | |
CN112731063B (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 |