CN105067966B - 基于特征模态分量能量分析的低压交流故障电弧检测方法 - Google Patents

基于特征模态分量能量分析的低压交流故障电弧检测方法 Download PDF

Info

Publication number
CN105067966B
CN105067966B CN201510398287.2A CN201510398287A CN105067966B CN 105067966 B CN105067966 B CN 105067966B CN 201510398287 A CN201510398287 A CN 201510398287A CN 105067966 B CN105067966 B CN 105067966B
Authority
CN
China
Prior art keywords
current
fault arc
half period
detection method
fault
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
CN201510398287.2A
Other languages
English (en)
Other versions
CN105067966A (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201510398287.2A priority Critical patent/CN105067966B/zh
Publication of CN105067966A publication Critical patent/CN105067966A/zh
Application granted granted Critical
Publication of CN105067966B publication Critical patent/CN105067966B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Relating To Insulation (AREA)

Abstract

本发明涉及一种基于特征模态分量能量分析的低压交流故障电弧检测方法,该方法采用经验模态分解获得采样电流的本征模态函数,再对所述本征模态函数进行后分解得到原始序列的特征模态分量I1,经Hilbert变换提取所述特征模态分量I1的瞬时幅值分布,并计算特征模态分量I1的各个半周期能量,通过所述半周期能量与设定参考阈值的比较,实现对串联故障电弧半周波的识别。与现有技术相比,本发明具有识别度高、多种类型负载适用等优点。

Description

基于特征模态分量能量分析的低压交流故障电弧检测方法
技术领域
本发明涉及一种低压交流串联故障电弧的检测方法,尤其是涉及一种基于特征模态分量能量分析的低压交流故障电弧检测方法。
背景技术
建筑电气电路长期带载或过载运行,电气接触不良等会造成线路绝缘破损、老化、碳化,进而造成电弧隐患,电弧局部电阻增大,产生热量加剧线路恶化,甚至引燃周围物体。
国内外现有的低压故障电弧断路器(AFCI)的内置算法多是基于对电路中故障电弧电流、电压信号的时域、频域特征分析。
考察电气量时域特征时,有学者利用卡尔曼滤波器的学习能力,通过比较习得信号与时域实际信号的差值大小判断电弧是否发生,判断结果准确,但其学习过程耗时,不满足实时检测要求。另有学者把电流的时域信号特征值作为BP神经网络的输入之一,把时域法与其他智能算法相结合。还有学者通过比较正常时和故障时每半周期电流峰值或平均值变化量来做判断,但局限于某些具体负载,缺乏通用性。
研究电流信号的频域特征时,有学者提出用快速傅里叶分析法,取电流五次谐波的幅值大小变化作为电弧发生的判据,但通过实验测量和数据分析,对于电吹风高档、电容型启动风扇负载该判断方法适用,对于电水壶等阻性负载,五次谐波的幅值大小在正常运行和故障时区别不大,对于电脑等开关电源型负载,甚至正常时的五次谐波的幅值高于故障时的,可见该判断方法对不同负载的适用有局限性。有研究人员用离散小波分解的方法,提取经若干层分解后的电流信号细节,结合神经网络学习方法,用其幅值变化或与横轴围成的面积大小判断电弧发生,具有一定的负载通用性,但自学习和样本训练过程消耗线下时间,且不同待测对象需要分别训练学习。一些满足计算机仿真的算法有一定的计算复杂度而在微处理器中移植受限。有学者把电弧视为一种局部放电现象,用高频电流互感器、射频线圈提取电弧发生时的高频信号,并分析其所在频段,制作专用频带的天线来检测电弧,但该方法易受外界空间信号干扰,有一定的误判率。
现有的故障电弧判别算法受负载特性影响较大,对线性负载时的故障电弧识别准确,而对现代电器中应用较多的感性负载、开关电源负载发生电弧故障时识别有困难,缺乏通用性。
希尔伯特黄变换(Hilbert-Huang transform,HHT)包含经验模态分解(EmpiricalMode Decomposition,EMD)技术和Hilbert变换技术,处理数据时,先通过EMD将信号分解为若干个本征模态函数(Intrinsic Mode Functions,IMF),再经Hilbert变换得到时间序列的各个IMF的瞬时幅值、瞬时频率和Hilbert谱。通过研究HHT所得的特征与分布来辨识输入时间序列的参数和特性。HHT广泛应用于工程体系结构参数识别、重大工程(如:高层建筑、桥梁、道路等)变形监测、振动分析、机械故障诊断、信号谐波分析等诸多领域,能够准确捕捉线性平稳信号的特征,也适用于分析非线性、非稳态数据,具有识别信号的不同频率组分功能的同时,所得分量具有较明显的实际物理意义。
发明内容
本发明的目的就是为了克服上述现有技术存在的某些缺点而提供一种识别度高、多种类型负载适用的基于特征模态分量能量分析的低压交流故障电弧检测方法。
本发明的目的可以通过以下技术方案来实现:
一种基于特征模态分量能量分析的低压交流故障电弧检测方法,该方法采用经验模态分解获得采样电流的本征模态函数,再对所述本征模态函数进行后处理得到原始序列的特征模态分量I1,经Hilbert变换提取所述特征模态分量I1的瞬时幅值分布,并计算特征模态分量I1的半周期能量,通过所述半周期能量与设定参考阈值的比较,实现对串联故障电弧半周波的识别。
对采样电流进行经验模态分解前包括步骤:对所述采样电流进行归一化和软件降噪。
所述后处理具体为:
从本征模态函数中选取物理意义描述故障特征的多个固有模态分量,将其中波形趋势相似度最高的1~3个相邻分量进行组合,得到特征模态分量I1
通过所述半周期能量与设定参考阈值的比较,实现对串联故障电弧半周波的识别具体为:
1)初始化电弧故障发生标志位flag=0;
2)判断第j个半周期能量ej与设定参考阈值e0的大小关系,若ej>e0,则执行步骤3),若ej≤e0,则令j=j+1,返回步骤2);
3)电弧故障发生标志位flag=1,记为故障半周波。
所述设定参考阈值的计算方法是:
提取当前负载下正常运行时电流的特征模态分量,取该特征模态分量的p个半周期能量的平均值作为半周期能量的设定参考阈值,其中,erp为第p个半周期能量,p为选取的半周期的个数。
所述选取的半周期的个数p取为5。
所述半周期能量ej的计算公式为:
其中,N=f/100是半周期信号的采样点数,f为采样电流的采样频率,ai为特征模态分量I1第i个采样点处的瞬时幅值。
当第1个故障半周波与第8个故障半周波之间的时间差Δt满足Δt≤0.5s时,开断负载所在支路。
根据上述技术方案,基于特征模态分量能量分析的低压交流故障电弧检测方法的具体步骤可描述为:
a)按采样率采集负载所在电路的串联电流;
b)将采样电流x(t)进行归一化和软件降噪,得到y(t);
c)对y(t)进行经验模态分解得到各个本征模态分量IMF,对IMF进行后处理得到特征模态分量I1,作为故障特征分析对象;
d)对特征模态分量I1进行Hilbert变换,得到I1的瞬时幅值分布;
e)计算特征模态分量I1的半周期能量ej和半周期能量参考阈值e0
f)比较ej与参考阈值e0的大小关系,若ej≤e0,则令j=j+1,返回步骤e),若ej>e0,则故障电弧发生标志位flag=1,执行步骤g);
g)flag每置1一次,故障电弧个数count加1,记录该序列的序列号index,并把flag置0;
h)判断第1个故障半周波与第8个故障半周波之间的时间差Δt是否不大于0.5s,若是,转i),若否,转a);
i)开断负载所在支路。
与现有技术相比,本发明具有以下优点:
(1)本发明在进行EMD前对采样电流进行归一化、降噪等处理,减小电弧以外的因素对电流波形的影响,方便不同负载时的电流数据特征比较,提高对不同负载电流细节特点的有效识别度,有利于特征模态分量能量分析算法对串联故障电弧的准确识别。
(2)本发明方法适用实验条件下各类负载的串联故障电弧判断,设置适当的半周期能量参考阈值,可以实现多种类型负载下的串联故障电弧判断,具有一定的通用性。
(3)本发明方法对工频电流的半周期进行考察,满足UL1699标准对串联故障电弧半周波个数的计数要求,可以识别电弧发生初期,非连续状态的故障电弧,算法的逻辑复杂度小,不必区分前半周波和后半周波。
(4)本发明方法便于嵌入式系统的数据计算和高效传输,不需离线自学习或样本训练过程,缩短故障判断过程时间。
附图说明
图1是本发明检测方法的流程示意图;
图2是本发明检测方法在阻性负载1000W电水壶时串联故障电弧的判断过程示意图;
图3是本发明检测方法在感性负载500W电钻时串联故障电弧的判断过程示意图;
图4是本发明检测方法在容性负载50W电容启动型风扇时串联故障电弧的判断过程示意图;
图5是本发明检测方法在开关电源负载300W电脑时串联故障电弧的判断过程示意图;
图6是本发明检测方法在开关电源负载1200W复印机时串联故障电弧的判断过程示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
图1是本发明提出的基于特征模态分量能量分析的低压交流故障电弧检测方法的逻辑流程图,该方法的工作流程是:
开始→参数初始化,包括:电弧故障发生标志位flag,电弧故障个数count,跳闸标志位Trip均置0,设置电流采样率f→按既定采样率采样得到电流原始信号x(t)→采样电流x(t)的归一化和软件降噪,得到y(t)→对y(t)进行经验模态分解和后处理得到特征模态分量I1→通过Hilbert变换得到I1的瞬时幅值分布→计算特征模态分量I1的各个半周期能量→设置当前负载下半周期能量的参考阈值e0,即正常运行时电流的特征模态分量的5个半周期的能量的平均值→比较ej与参考阈值e0的大小关系,若ej≤e0,待测序列号j加1,转下一步骤,若ej>e0,则故障电弧发生标志位flag=1,转下一步骤→记录电弧故障个数,若flag值为1,故障电弧个数count加1,并记录该序列的序列号index,再把flag置0→判断第1个故障半周波与第8个故障半周波之间的时间差Δt是否不大于0.5s。若是,转下一步骤,若否,继续下一个上述采样计算过程→跳闸信号Trip=1,开断负载所在支路。
上述过程中,后处理具体为:从本征模态函数中选取物理意义描述故障特征的多个固有模态分量,将其中波形趋势相似度最高的1~3个相邻分量进行组合,得到特征模态分量I1
实施例1
本实施例将上述方法应用于负载1000W电水壶阻性的串联故障电弧识别中,具体步骤如下:
1)按采样率f=104Hz采集负载所在电路的串联电流。
2)将采样电流x(t)进行归一化和软件降噪,得到y(t)。
3)对y(t)进行经验模态分解得到各个本征模态分量IMF,对IMF进行后处理得到特征模态分量I1,作为故障特征分析对象。
4)对特征模态分量I1进行Hilbert变换,得到I1的瞬时幅值分布。
5)计算特征模态分量I1的半周期能量ej,其计算公式如下:
其中,N=f/100是半周期信号的采样点数,ai为特征模态分量I1第i个采样点处的瞬时幅值。
6)提取该负载下正常运行时电流的特征模态分量,取其5个半周期的能量的平均值作为半周期能量的参考阈值。
7)比较ej与参考阈值e0的大小关系,若ej≤e0,待测序列号j加1,若ej>e0故障电弧发生标志位flag=1。
8)flag每置1一次,故障电弧个数count加1,记录该序列的序列号index,并把flag置0。
9)判断第1个故障半周波与第8个故障半周波之间的时间差Δt是否不大于0.5s。若是,转10),若否,转1)。
10)跳闸信号标志位Trip=1,开断负载所在支路。
图2是本发明检测方法在阻性负载1000W电水壶的串联故障电弧的判断过程。图2中从上到下依次是电水壶负载归一化降噪后的电流y(t)、y(t)的特征模态分量I1、I1的归一化半周期能量Power、负载所在支路的跳闸信号标志位Trip。由图2可以看到,本发明的算法可以实现阻性负载电水壶时串联故障电弧(arc)的准确识别判断。
实施例2
本实施例将本发明检测方法应用于感性负载500W电钻的串联故障电弧识别中,识别的具体步骤同实施例1。图3显示了故障电弧的判断过程,从上到下依次是电钻负载归一化降噪后的电流y(t)、y(t)的特征模态分量I1、I1的归一化半周期能量Power、负载所在支路的跳闸信号标志位Trip。由图3可以看到,本发明的算法可以实现感性负载电钻时串联故障电弧的准确识别判断。
实施例3
本实施例将本发明检测方法应用于容性负载50W电容启动型风扇的串联故障电弧识别中,识别的具体步骤同实施例1。图4显示了故障电弧的判断过程,从上到下依次是电容启动型风扇负载归一化降噪后的电流y(t)、y(t)的特征模态分量I1、I1的归一化半周期能量Power、负载所在支路的跳闸信号标志位Trip。由图4可以看到,本发明的算法可以实现容性负载电容启动型风扇时串联故障电弧的准确识别判断。
实施例4
本实施例将本发明检测方法应用于开关电源负载300W电脑的串联故障电弧识别中,识别的具体步骤同实施例1。图5显示了故障电弧的判断过程,从上到下依次是电脑负载归一化降噪后的电流y(t)、y(t)的特征模态分量I1、I1的归一化半周期能量Power、负载所在支路的跳闸信号标志位Trip。由图5可以看到,本发明的算法可以实现开关电源负载电脑时串联故障电弧的准确识别判断。
实施例5
本实施例将本发明检测方法应用于开关电源负载1200W复印机的串联故障电弧识别中,识别的具体步骤同实施例1。图6显示了故障电弧的判断过程,从上到下依次是复印机负载归一化降噪后的电流y(t)、y(t)的特征模态分量I1、I1的归一化半周期能量Power、负载所在支路的跳闸信号标志位Trip。由图6可以看到,本发明的算法可以实现开关电源负载复印机时串联故障电弧的准确识别判断。
本领域技术人员应理解,本发明的保护范围不限于上述的实施例,任何通过考察采样电流与正常运行时电流的特征模态分量的半周期能量相对大小关系而做出故障电弧发生与否判断的做法均属于本发明检测方法之内。

Claims (7)

1.一种基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,该方法采用经验模态分解获得采样电流的本征模态函数,再对所述本征模态函数进行后处理得到原始序列的特征模态分量I1,经Hilbert变换提取所述特征模态分量I1的瞬时幅值分布,并计算特征模态分量I1的各个半周期能量,通过所述半周期能量与设定参考阈值的比较,实现对串联故障电弧半周波的识别。
2.根据权利要求1所述的基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,对采样电流进行经验模态分解前包括步骤:
对所述采样电流进行归一化和软件降噪。
3.根据权利要求1所述的基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,所述后处理具体为:
从本征模态函数中选取物理意义描述故障特征的多个固有模态分量,将其中波形趋势相似度最高的1~3个相邻分量进行组合,得到特征模态分量I1
4.根据权利要求1所述的基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,通过所述半周期能量与设定参考阈值的比较,实现对串联故障电弧半周波的识别具体为:
1)初始化电弧故障发生标志位flag=0;
2)判断第j个半周期能量ej与设定参考阈值e0的大小关系,若ej>e0,则执行步骤3),若ej≤e0,则令j=j+1,返回步骤2);
3)电弧故障发生标志位flag=1,记为串联故障电弧半周波。
5.根据权利要求1或4所述的基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,所述设定参考阈值的计算方法是:
提取当前负载下正常运行时电流的特征模态分量,取该特征模态分量的p个半周期能量的平均值作为半周期能量的设定参考阈值,其中,erp为第p个半周期能量,p为选取的半周期的个数。
6.根据权利要求5所述的基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,所述选取的半周期的个数p取为5。
7.根据权利要求4所述的基于特征模态分量能量分析的低压交流故障电弧检测方法,其特征在于,当第1个串联故障电弧半周波与第8个串联故障电弧半周波之间的时间差Δt满足Δt≤0.5s时,断开负载所在支路。
CN201510398287.2A 2015-07-08 2015-07-08 基于特征模态分量能量分析的低压交流故障电弧检测方法 Active CN105067966B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510398287.2A CN105067966B (zh) 2015-07-08 2015-07-08 基于特征模态分量能量分析的低压交流故障电弧检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510398287.2A CN105067966B (zh) 2015-07-08 2015-07-08 基于特征模态分量能量分析的低压交流故障电弧检测方法

Publications (2)

Publication Number Publication Date
CN105067966A CN105067966A (zh) 2015-11-18
CN105067966B true CN105067966B (zh) 2017-10-10

Family

ID=54497382

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510398287.2A Active CN105067966B (zh) 2015-07-08 2015-07-08 基于特征模态分量能量分析的低压交流故障电弧检测方法

Country Status (1)

Country Link
CN (1) CN105067966B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107783013A (zh) * 2016-08-30 2018-03-09 上海复旦微电子集团股份有限公司 一种电弧故障的检测方法及装置
CN107797031B (zh) * 2016-08-30 2021-01-19 上海复旦微电子集团股份有限公司 电弧故障的检测方法及装置
GB201617809D0 (en) * 2016-10-21 2016-12-07 Analog Devices Global A method of detecting arc events in a power system, and a power system including an arc detector
US10509067B2 (en) * 2017-07-06 2019-12-17 Mersen Usa Newburyport-Ma, Llc Method for AC arc fault detection using multidimensional energy points
CN109239560B (zh) * 2018-11-29 2021-04-13 浙江习羽智能科技有限公司 故障电弧检测方法、装置及存储介质
CN109975673B (zh) * 2019-04-23 2021-03-16 辽宁工程技术大学 一种光伏微电网直流侧故障电弧识别方法
CN110531142A (zh) * 2019-08-13 2019-12-03 珠海格力电器股份有限公司 一种空调电弧故障的检测方法
CN110646714A (zh) * 2019-09-26 2020-01-03 杭州电力设备制造有限公司 一种超声波局放信号的特征提取方法及相关设备
CN110740063B (zh) * 2019-10-25 2021-07-06 电子科技大学 基于信号分解和周期特性的网络流量特征指标预测方法
CN110854821B (zh) * 2019-11-15 2021-09-24 上海红檀智能科技有限公司 一种基于电流及电弧光变化趋势的自适应电弧光保护方法
CN113484658B (zh) * 2021-05-26 2022-06-28 上海红檀智能科技有限公司 电弧故障的诊断方法、系统、介质及电子设备
CN113589105B (zh) * 2021-08-09 2022-07-26 山东建筑大学 基于物联网的建筑电气火灾串联故障电弧识别方法及系统
CN113608084B (zh) * 2021-08-12 2024-08-16 山东理工大学 一种串联故障电弧检测方法
TWI837770B (zh) * 2022-08-17 2024-04-01 國立臺灣科技大學 電弧故障檢測模型的訓練及線上更新方法
CN116298735A (zh) * 2023-05-10 2023-06-23 广东电网有限责任公司湛江供电局 一种低压配电网交流电弧故障检测方法及相关装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5132629A (en) * 1991-02-25 1992-07-21 Clinton Instrument Company Apparatus for testing the insulation of an electrical conductor passing through an electrode
CN101551433A (zh) * 2009-05-05 2009-10-07 昆明理工大学 一种利用hht检测的配电网馈出线路故障选线方法
CN101814731A (zh) * 2010-04-15 2010-08-25 中国石油大学(华东) 一种应用非工频暂态分量的配电网故障选线方法
CN104635133A (zh) * 2015-03-14 2015-05-20 北京芯同汇科技有限公司 一种新型故障电弧检测方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI453432B (zh) * 2012-11-28 2014-09-21 Simplo Technology Co Ltd 絕緣偵測電路及其方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5132629A (en) * 1991-02-25 1992-07-21 Clinton Instrument Company Apparatus for testing the insulation of an electrical conductor passing through an electrode
CN101551433A (zh) * 2009-05-05 2009-10-07 昆明理工大学 一种利用hht检测的配电网馈出线路故障选线方法
CN101814731A (zh) * 2010-04-15 2010-08-25 中国石油大学(华东) 一种应用非工频暂态分量的配电网故障选线方法
CN104635133A (zh) * 2015-03-14 2015-05-20 北京芯同汇科技有限公司 一种新型故障电弧检测方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
关于检测低压故障电弧的初步研究;印永嘉 等;《实验室研究与探索》;20070815;第26卷(第8期);第20-23页 *

Also Published As

Publication number Publication date
CN105067966A (zh) 2015-11-18

Similar Documents

Publication Publication Date Title
CN105067966B (zh) 基于特征模态分量能量分析的低压交流故障电弧检测方法
CN107025365B (zh) 一种用于用户侧的非侵入负荷辨识方法
CN106786534B (zh) 一种非侵入式电力负荷暂态过程辨识方法及系统
CN104597344A (zh) 基于小波一层高频分量相关性的故障电弧在线检测方法
Zhang et al. A real-time classification method of power quality disturbances
Liu et al. Application of the variational mode decomposition-based time and time–frequency domain analysis on series DC arc fault detection of photovoltaic arrays
CN102331543A (zh) 基于支持向量机的故障电弧检测方法
CN106908671A (zh) 一种非侵入式家用负荷智能检测方法及系统
CN106250904B (zh) 基于改进s变换的电能扰动分析仪及分类方法
CN109447473B (zh) 一种电力负荷监测方法、装置、设备及可读存储介质
CN109470957B (zh) 一种基于弗雷歇距离的非侵入式电器识别方法
CN112505511A (zh) 一种非侵入式低压故障电弧检测与定位方法及系统
CN114252682B (zh) 基于非侵入式负荷监测的大型楼宇异常负荷控制方法
CN112183628B (zh) 基于多种线性时频变换的交流电弧故障检测方法及系统
CN113567794A (zh) 基于动态时间规整的电动自行车户内充电辨识方法及系统
CN112949563B (zh) 基于变点检测和改进knn算法的非侵入式负荷识别方法
Gong et al. Series arc fault identification method based on wavelet transform and feature values decomposition fusion DNN
CN117370777A (zh) 一种电动自行车充电负荷辨识特征的选择方法及系统
Lu et al. A new series arc fault identification method based on wavelet transform
Alaca et al. CNN-Based Phase Fault Classification in Real and Simulated Power Systems Data
Zeng et al. High impedance fault detection in distribution network based on S-transform and average singular entropy
CN110542855B (zh) 基于离散余弦变换的负荷开关事件检测方法和系统
Li et al. The investigation of residential load identification technologies
CN110749841A (zh) 一种利用转换空间因子的负荷开关事件检测方法和系统
Liu et al. Dc series arc fault detection based on random forest combined with entropy weight method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant