CN101793918A - 一种电压暂降检测方法 - Google Patents

一种电压暂降检测方法 Download PDF

Info

Publication number
CN101793918A
CN101793918A CN 200910237554 CN200910237554A CN101793918A CN 101793918 A CN101793918 A CN 101793918A CN 200910237554 CN200910237554 CN 200910237554 CN 200910237554 A CN200910237554 A CN 200910237554A CN 101793918 A CN101793918 A CN 101793918A
Authority
CN
China
Prior art keywords
component
voltage
sequence
amplitude
positive
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
CN 200910237554
Other languages
English (en)
Other versions
CN101793918B (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 Electric Power Research Institute Co Ltd CEPRI
China EPRI Science and Technology Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
China EPRI Science and Technology 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 China Electric Power Research Institute Co Ltd CEPRI, China EPRI Science and Technology Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN 200910237554 priority Critical patent/CN101793918B/zh
Publication of CN101793918A publication Critical patent/CN101793918A/zh
Application granted granted Critical
Publication of CN101793918B publication Critical patent/CN101793918B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Current Or Voltage (AREA)

Abstract

本发明涉及一种电压暂降检测方法,(1)若发生电压跌落故障时,将包含正序、负序和五次谐波分量的系统三相电压送入dq坐标变换模块;(2)经坐标变换后的正序分量为直流分量、负序分量为100Hz分量、五次谐波分量为300Hz分量,通过等值滤波网络将100Hz和300Hz分量滤除,得到d轴和q轴的直流分量,对两直流分量进行平方和相加、再取均方根运算,得到系统正序电压幅值;(3)选取额定电压幅值的90%为电压跌落的阈值,将系统正序电压幅值与该电压跌落的阈值比较,当正序电压幅值小于电压跌落的阈值时,电压发生跌落;当正序电压幅值大于等于电压跌落的阈值时,则没有发生电压跌落。本发明由于采用了等值滤波网络,因此解决了谐波放大问题。本发明可以广泛应用于各种电力系统的电压故障检测中。

Description

一种电压暂降检测方法
技术领域
本发明涉及一种电力系统中的检测方法,特别是关于一种电压暂降检测方法。
背景技术
随着现代科学技术的发展,越来越多的工业生产过程和流水线(如半导体制造、计算机集成制造、造纸业等)依赖于对电能质量十分敏感的以微处理器芯片为核心的设备,即使几个周期的供电中断或电压骤降都将影响这些设备正常工作。因此高质量电力供应已成为现代社会生产、生活正常进行的基本条件。目前,在对电能质量问题的研究中,电压骤降已被认为是影响许多用电设备正常、安全运行的最严重的动态电能质量问题。统计表明,大型电力用户,幅度超过20%的骤降年发生率在10到20次左右,许多高度自动化连续生产过程,每次电压骤降造成的经济损失达数十万至数百万美元之多。考虑到电压跌落发生的随机性和快速性,要使动态电能质量调节装置具有良好的实时控制效果,首先要解决的是在保证能对装置的控制信号(通常为电压、电流)在一定检测准确度的前提下实现快速跟踪检测问题。
目前研究较多的方法,主要有有效值法、基于瞬时无功功率理论的dq0变换方法和小波分析法等。但这些方法要么实时性差,要么运算复杂,因此都仅限用于软件仿真和试验样机模拟。其中,基于“abc-dq”变换的检测算法是目前被研究最多的算法,该方法的常规算法需要通过低通滤波器或滑动窗来分离直流分量,存在延时,难以保证实时性。随后,有人对“abc-dq”变换的检测算法进行了改进,通过微分运算代替滤波器来分离直流分量,提高了实时性,但是该改进算法存在谐波放大的缺点。
发明内容
针对上述问题,本发明的目的是提供一种实时性较高、并能避免谐波放大的电压暂降检测方法。
为实现上述目的,本发明采取以下技术方案:一种电压暂降检测方法,其步骤如下:(1)假定发生电压跌落故障时,将包含正序分量、负序分量和五次谐波分量的系统三相电压送入dq坐标变换模块内;(2)经坐标变换后得到的正序分量为直流分量、负序分量为100Hz分量、五次谐波分量为300Hz分量,通过等值滤波网络将100Hz分量和300Hz分量滤除,得到d轴和q轴的直流分量,并对两直流分量进行平方和相加、再取均方根运算,得到系统正序电压的幅值;(3)选取额定电压幅值的90%为电压跌落的阈值,将系统正序电压幅值与该电压跌落的阈值进行比较,当正序电压幅值小于电压跌落的阈值时,则电压发生跌落;当正序电压幅值大于等于电压跌落的阈值时,则没有发生电压跌落。
所述步骤(2)中,所述d轴直流分量和所述q轴直流分量分别为:Vd=Vd+Vq_div=V1cos(θ1)+V′h5cos(6ωt-θ′h5),Vq=Vq-Vd_div=V1sin(θ1)-V′h5sin(6ωt-θ′h5),式中,Vd为所述d轴直流电压分量;Vq为所述q轴直流电压分量;Vd_div和Vq_div分别为d轴电压幅值Vd、q轴电压幅值Vq经所述等值滤波网络后的输出电压;V′h5为运算后剩余五次谐波幅值;θ′h5为运算后剩余五次谐波分量相角。
所述等值滤波网络在中频段时,其增益为1,相移为90度;在大于100Hz的高频段时,谐波分量稳定或减小。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于采用dq变换模块,将系统发生电压跌落故障时的三相电压经dq变换后,得到的直流分量、100Hz分量和300Hz分量经过等值滤波网络后,输出的直流分量与电压跌落的阈值进行比较,进而判断是否发生电压跌落,因此保证了检测算法的快速性,避免了谐波放大的缺点,并且实用性较高。2、本发明由于采用了等值滤波网络,该等值滤波网络在小于100Hz的低频段,其衰减特性与现有技术采用微分环节在低频段的衰减特性相同;在中频段,其增益为1,相移为90度;在大于100Hz的高频段,谐波分量不被放大。因此保证了高频段不再被放大,也就是解决了谐波放大问题。本发明可以广泛应用于各种电力系统的电压故障检测中。
附图说明
图1是本发明的电压跌落检测原理框图
图2是本发明的等值滤波网络的频率特性伯德图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,在系统发生电压跌落故障时,本发明的电压暂降检测方法,其步骤如下:
1)假定发生电压跌落故障时,系统三相电压包含正序分量、负序分量和五次谐波分量,并将该系统三相电压送入dq变换模块内,进行坐标变换;
2)经坐标变换后得到的正序分量变为直流分量、负序分量为100Hz分量、五次谐波分量为300Hz分量,通过等值滤波网络将100Hz分量和300Hz分量滤除,得到d轴和q轴的直流分量,并对d轴和q轴直流分量进行平方和相加、再取均方根运算,得到系统正序电压的幅值Vrms
3)选取额定电压幅值的90%为电压跌落的阈值,将系统正序电压幅值Vrms与该电压跌落的阈值进行比较,当正序电压幅值Vrms小于电压跌落的阈值时,则电压发生跌落;当正序电压幅值Vrms大于等于电压跌落的阈值时,则没有发生电压跌落。
上述步骤1)中,三相电压为:
Va=V1cos(ωt+θ1)+V2cos(ωt+θ2)+Vh5(5ωt+θh5)
Vb=V1cos(ωt+θ1-120°)+V2cos(ωt+θ2+120°)+Vh5(5ωt+θh5+120°)    (1)
Vc=V1cos(ωt+θ1+120°)+V2cos(ωt+θ2-120°)+Vh5(5ωt+θh5-120°)
公式(1)中的三相电压经32变换后为:
V α V β = 2 / 3 1 - 1 / 2 - 1 / 2 0 3 / 2 - 3 / 2 V a V b V c - - - ( 2 )
公式(2)经dq变换后为:
V d V q = cos ωt sin ωt - sin ωt cos ωt V α V β = V 1 cos ( θ 1 ) + V 2 cos ( 2 ωt - θ 2 ) + V h 5 cos ( 6 ωt - θ h 5 ) V 1 sin ( θ 1 ) - V 2 sin ( 2 ωt - θ 2 ) - V h 5 sin ( 6 ωt - θ h 5 ) - - - ( 3 )
式中,V1为系统电压正序分量幅值;θ1为系统电压正序分量初始相角;V2为系统电压负序分量幅值;θ2为系统电压负序分量初始相角;Vh5为系统电压五次谐波幅值;θh5为系统电压五次谐波分量初始相角;ω为角频率;
上述步骤2)中,经等值滤波网络输出的电压与经公式(3)运算后得到的的d轴和q轴电压相加后,得到的d轴直流分量Vd和q轴直流分量Vq分别为:
Vd=Vd+Vq_div=V1cos(θ1)+V′h5cos(6ωt-θ′h5)    (4)
Vq=Vq-Vd_div=V1sin(θ1)-V′h5sin(6ωt-θ′h5)    (5)
根据公式(4)和公式(5)得到系统正序电压的幅值Vrms为:
V rms = V ‾ d 2 + V ‾ q 2
= [ V 1 cos ( θ 1 ) + V h 5 ′ cos ( 6 ωt - θ h 5 ′ ) ] 2 + [ V 1 sin ( θ 1 ) - V h 5 ′ sin ( 6 ωt - θ h 5 ′ ) ] 2 - - - ( 6 )
= V 1 2 + V h 5 ′ 2 + cos ( 6 ωt + θ 1 - θ h 5 ′ )
式中,Vd_div和Vq_div分别为d轴电压幅值Vd、q轴电压幅值Vq经等值滤波网络后的输出电压;V′h5为运算后剩余五次谐波幅值;θ′h5为运算后剩余五次谐波分量相角。
由公式(4)和公式(5)可知V′h5和θ′h5与等值滤波网络的高频特性相关。本发明采用的等值滤波网络的特性为(如图2所示):其中A点处频率为101Hz,相移为89.9度;B点处频率为101Hz,幅值衰减为0.0424dB;C点处频率为301Hz,幅值衰减为7.94dB。因此可以看出本发明的等值滤波网络在小于100Hz的低频段时,其衰减特性与现有技术采用微分环节在低频段的衰减特性相同;在中频段,其增益为1,相移为90度;在大于100Hz的高频段,其衰减特性可以通过改变等值滤波网络相应系数进行控制,以使得运算后剩余五次谐波幅值V′h5值很小,保证了谐波分量不被放大。
上述各实施例仅为本发明的应用,并非用于限定本发明的实施范围。凡基于本发明技术方案上的变化和改进,不应排除在本发明的保护范围之外。

Claims (3)

1.一种电压暂降检测方法,其步骤如下:
(1)假定发生电压跌落故障时,将包含正序分量、负序分量和五次谐波分量的系统三相电压送入dq坐标变换模块内;
(2)经坐标变换后得到的正序分量为直流分量、负序分量为100Hz分量、五次谐波分量为300Hz分量,通过等值滤波网络将100Hz分量和300Hz分量滤除,得到d轴和q轴的直流分量,并对两直流分量进行平方和相加、再取均方根运算,得到系统正序电压的幅值;
(3)选取额定电压幅值的90%为电压跌落的阈值,将系统正序电压幅值与该电压跌落的阈值进行比较,当正序电压幅值小于电压跌落的阈值时,则电压发生跌落;当正序电压幅值大于等于电压跌落的阈值时,则没有发生电压跌落。
2.如权利要求1所述的一种电压暂降检测方法,其特征在于:所述步骤(2)中,所述d轴直流分量和所述q轴直流分量分别为:
Vd=Vd+Vq_div=V1cos(θ1)+V′h5cos(6ωt-θ′h5),
Vq=Vq-Vd_div=V1sin(θ1)-V′h5sin(6ωt-θ′h5),
式中,Vd为所述d轴直流电压分量;Vq为所述q轴直流电压分量;Vd_div和Vq_div分别为d轴电压幅值Vd、q轴电压幅值Vq经所述等值滤波网络后的输出电压;V′h5为运算后剩余五次谐波幅值;θ′h5为运算后剩余五次谐波分量相角。
3.如权利要求1或2所述的一种电压暂降检测方法,其特征在于:所述等值滤波网络在中频段时,其增益为1,相移为90度;在大于100Hz的高频段时,谐波分量稳定或减小。
CN 200910237554 2009-11-18 2009-11-18 一种电压暂降检测方法 Active CN101793918B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200910237554 CN101793918B (zh) 2009-11-18 2009-11-18 一种电压暂降检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200910237554 CN101793918B (zh) 2009-11-18 2009-11-18 一种电压暂降检测方法

Publications (2)

Publication Number Publication Date
CN101793918A true CN101793918A (zh) 2010-08-04
CN101793918B CN101793918B (zh) 2012-10-03

Family

ID=42586724

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200910237554 Active CN101793918B (zh) 2009-11-18 2009-11-18 一种电压暂降检测方法

Country Status (1)

Country Link
CN (1) CN101793918B (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072983A (zh) * 2010-11-22 2011-05-25 华北电力大学(保定) 一种电压暂降起因监判方法
CN102095915A (zh) * 2010-11-24 2011-06-15 重庆大学 一种采用多同步参考坐标系变换的电压信号检测装置
CN102156217A (zh) * 2011-03-31 2011-08-17 山亿新能源股份有限公司 一种基于平均值滤波的三相电网电压跌落检测方法
CN102854421A (zh) * 2012-09-11 2013-01-02 江苏旭源科技有限公司 一种光伏逆变器低电压穿越的快速判定方法
CN103018532A (zh) * 2012-12-13 2013-04-03 南京安炤电力电子有限公司 高速电压跌落检测方法及系统
CN103163363A (zh) * 2013-03-11 2013-06-19 西安交通大学 用于动态电压恢复器的电网电压跌落检测算法
CN104360137A (zh) * 2014-11-17 2015-02-18 华北电力大学 一种适用于动态电压恢复器的电压暂降检测方法
CN104749453A (zh) * 2013-12-30 2015-07-01 上海宝钢工业技术服务有限公司 降低外网单相接地故障对用户电压暂降影响的方法
CN104993711A (zh) * 2015-05-22 2015-10-21 国网河南省电力公司电力科学研究院 一种电压暂降过渡过程模拟装置及方法
CN105116195A (zh) * 2015-07-13 2015-12-02 中国人民解放军海军工程大学 适用于电网电压谐波含量大的电压跌落检测方法
CN105137164A (zh) * 2015-08-06 2015-12-09 江苏省电力公司苏州供电公司 应用于电力系统中的电压暂降在线监测装置
CN105388396A (zh) * 2015-11-04 2016-03-09 中国矿业大学 一种用序有功增量电流方向追溯电压暂降源的方法
CN107085136A (zh) * 2017-03-15 2017-08-22 中国电力科学研究院 一种三相交流电压跌落状态的判断方法及系统
CN108169544A (zh) * 2017-12-25 2018-06-15 国网冀北电力有限公司秦皇岛供电公司 电压暂降展示方法及装置
CN110045175A (zh) * 2019-01-29 2019-07-23 广东电网有限责任公司 一种单相配电系统电压跌落检测方法
CN110749769A (zh) * 2019-11-20 2020-02-04 山东泰开自动化有限公司 一种三相电压跌落的快速检测方法
CN112510717A (zh) * 2020-11-04 2021-03-16 天合光能股份有限公司 大功率储能双向变流器零电压穿越控制方法
CN112986680A (zh) * 2021-02-10 2021-06-18 广东电网有限责任公司 一种基于坐标变换的电压暂降快速检测方法和系统
CN113376424A (zh) * 2021-05-20 2021-09-10 国网湖北省电力有限公司电力科学研究院 用于级联式电能质量综合治理装置的电压暂降检测方法
CN113991835A (zh) * 2021-11-29 2022-01-28 国网江苏省电力有限公司扬州市江都区供电分公司 一种基于备自投装置的双电源无缝切换的控制方法
CN114325070A (zh) * 2021-12-13 2022-04-12 广西电网有限责任公司电力科学研究院 一种极端电网工况下高鲁棒电压暂降检测方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101452016A (zh) * 2007-11-30 2009-06-10 北京市电力公司 电压暂降检测方法和装置
CN101339208B (zh) * 2008-08-12 2010-07-28 中国矿业大学 一种基于时域分析的电压质量监测与扰动自动分类方法
CN101487861B (zh) * 2009-02-27 2011-08-31 国电南瑞科技股份有限公司 电网电压跌落时刻电压相角跳变检测方法

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072983A (zh) * 2010-11-22 2011-05-25 华北电力大学(保定) 一种电压暂降起因监判方法
CN102095915A (zh) * 2010-11-24 2011-06-15 重庆大学 一种采用多同步参考坐标系变换的电压信号检测装置
CN102095915B (zh) * 2010-11-24 2013-04-17 重庆大学 一种采用多同步参考坐标系变换的电压信号检测装置
CN102156217A (zh) * 2011-03-31 2011-08-17 山亿新能源股份有限公司 一种基于平均值滤波的三相电网电压跌落检测方法
CN102854421A (zh) * 2012-09-11 2013-01-02 江苏旭源科技有限公司 一种光伏逆变器低电压穿越的快速判定方法
CN103018532A (zh) * 2012-12-13 2013-04-03 南京安炤电力电子有限公司 高速电压跌落检测方法及系统
CN103018532B (zh) * 2012-12-13 2014-12-31 南京安炤电力电子有限公司 高速电压跌落检测方法及系统
CN103163363B (zh) * 2013-03-11 2015-05-27 西安交通大学 用于动态电压恢复器的电网电压跌落检测算法
CN103163363A (zh) * 2013-03-11 2013-06-19 西安交通大学 用于动态电压恢复器的电网电压跌落检测算法
CN104749453A (zh) * 2013-12-30 2015-07-01 上海宝钢工业技术服务有限公司 降低外网单相接地故障对用户电压暂降影响的方法
CN104749453B (zh) * 2013-12-30 2019-02-15 上海宝钢工业技术服务有限公司 降低外网单相接地故障对用户电压暂降影响的方法
CN104360137B (zh) * 2014-11-17 2017-06-09 华北电力大学 一种适用于动态电压恢复器的电压暂降检测方法
CN104360137A (zh) * 2014-11-17 2015-02-18 华北电力大学 一种适用于动态电压恢复器的电压暂降检测方法
CN104993711A (zh) * 2015-05-22 2015-10-21 国网河南省电力公司电力科学研究院 一种电压暂降过渡过程模拟装置及方法
CN104993711B (zh) * 2015-05-22 2018-01-30 国网河南省电力公司电力科学研究院 一种电压暂降过渡过程模拟装置及方法
CN105116195A (zh) * 2015-07-13 2015-12-02 中国人民解放军海军工程大学 适用于电网电压谐波含量大的电压跌落检测方法
CN105116195B (zh) * 2015-07-13 2018-07-06 中国人民解放军海军工程大学 适用于电网电压谐波含量大的电压跌落检测方法
CN105137164A (zh) * 2015-08-06 2015-12-09 江苏省电力公司苏州供电公司 应用于电力系统中的电压暂降在线监测装置
CN105388396B (zh) * 2015-11-04 2018-10-26 中国矿业大学 一种用序有功增量电流方向追溯电压暂降源的方法
CN105388396A (zh) * 2015-11-04 2016-03-09 中国矿业大学 一种用序有功增量电流方向追溯电压暂降源的方法
CN107085136A (zh) * 2017-03-15 2017-08-22 中国电力科学研究院 一种三相交流电压跌落状态的判断方法及系统
CN108169544A (zh) * 2017-12-25 2018-06-15 国网冀北电力有限公司秦皇岛供电公司 电压暂降展示方法及装置
CN110045175A (zh) * 2019-01-29 2019-07-23 广东电网有限责任公司 一种单相配电系统电压跌落检测方法
CN110749769A (zh) * 2019-11-20 2020-02-04 山东泰开自动化有限公司 一种三相电压跌落的快速检测方法
CN112510717A (zh) * 2020-11-04 2021-03-16 天合光能股份有限公司 大功率储能双向变流器零电压穿越控制方法
CN112510717B (zh) * 2020-11-04 2023-07-18 江苏天合储能有限公司 大功率储能双向变流器零电压穿越控制方法
CN112986680A (zh) * 2021-02-10 2021-06-18 广东电网有限责任公司 一种基于坐标变换的电压暂降快速检测方法和系统
CN113376424A (zh) * 2021-05-20 2021-09-10 国网湖北省电力有限公司电力科学研究院 用于级联式电能质量综合治理装置的电压暂降检测方法
CN113376424B (zh) * 2021-05-20 2022-04-22 国网湖北省电力有限公司电力科学研究院 用于级联式电能质量综合治理装置的电压暂降检测方法
CN113991835A (zh) * 2021-11-29 2022-01-28 国网江苏省电力有限公司扬州市江都区供电分公司 一种基于备自投装置的双电源无缝切换的控制方法
CN114325070A (zh) * 2021-12-13 2022-04-12 广西电网有限责任公司电力科学研究院 一种极端电网工况下高鲁棒电压暂降检测方法及系统
CN114325070B (zh) * 2021-12-13 2023-09-05 广西电网有限责任公司电力科学研究院 一种极端电网工况下高鲁棒电压暂降检测方法及系统

Also Published As

Publication number Publication date
CN101793918B (zh) 2012-10-03

Similar Documents

Publication Publication Date Title
CN101793918B (zh) 一种电压暂降检测方法
CN102305886B (zh) 电网电压谐波畸变及不平衡时基波电压同步信号检测方法
CN102385003B (zh) 检测电力系统中电压信号或电流信号频率和相位的方法
CN103163363B (zh) 用于动态电压恢复器的电网电压跌落检测算法
CN103269086B (zh) 光伏并网逆变器低电压穿越控制的正负序分量分离方法
CN103149494B (zh) 三相电网电压短路故障诊断方法
CN103399202B (zh) 可变速抽水蓄能机组控制系统的锁相方法
CN104811188A (zh) 基于滑动滤波器的锁相环动态性能改进方法
CN107632200B (zh) 基于离散傅里叶变换的谐波检测方法
CN108627731A (zh) 一种单相断电的快速检测方法
CN102868183A (zh) 单相并网逆变器的基于多谐振滑模面的滑模变结构控制方法
CN105425171A (zh) 一种变频器低电压穿越电源电压跌落检测方法
CN102832643A (zh) 一种基于逆系统的三相光伏并网逆变器的控制方法
CN103872690A (zh) 一种基于hht 检测法和pfc 的动态电压恢复器控制方法
CN103117562A (zh) 一种高压级联能量回馈变频器功率模块的控制方法
CN103647550A (zh) 一种用于动态电压无功补偿的锁相环方法
CN110361629A (zh) 一种基于动态电压恢复器的电压暂降检测系统
CN105024393A (zh) 一种直流输电系统换相面积计算方法及换相失败检测方法
CN102868163B (zh) 一种自适应并网变流器三相软锁相环
CN106374917A (zh) 一种适用于电压暂降情况下的锁相环实现方法
CN104459389A (zh) 一种输入掉电晃电的判断方法与系统
CN103995180A (zh) 一种考虑不等式约束的电力系统频率估计方法
CN205389108U (zh) 一种微电网孤岛运行检测装置
CN107703358A (zh) 一种基于改进二阶广义积分器的锁相算法
CN103560509B (zh) 一种基于小波分析的电压下陷检测装置及该装置的控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant