CN104333283A - 基于回路成形的双馈电机定子电流鲁棒控制方法 - Google Patents

基于回路成形的双馈电机定子电流鲁棒控制方法 Download PDF

Info

Publication number
CN104333283A
CN104333283A CN201410567661.2A CN201410567661A CN104333283A CN 104333283 A CN104333283 A CN 104333283A CN 201410567661 A CN201410567661 A CN 201410567661A CN 104333283 A CN104333283 A CN 104333283A
Authority
CN
China
Prior art keywords
stator
stator current
current
rotor
electric engine
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
CN201410567661.2A
Other languages
English (en)
Other versions
CN104333283B (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 South Power Grid International Co ltd
Original Assignee
China South Power Grid International 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 South Power Grid International Co ltd filed Critical China South Power Grid International Co ltd
Priority to CN201410567661.2A priority Critical patent/CN104333283B/zh
Publication of CN104333283A publication Critical patent/CN104333283A/zh
Application granted granted Critical
Publication of CN104333283B publication Critical patent/CN104333283B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Ac Motors In General (AREA)

Abstract

本发明公开了一种基于回路成形的双馈电机定子电流鲁棒控制方法。该方法主要包含H∞电流调节器设计与双馈电机定子电流控制。本发明通过建立双馈电机的数学模型,将双馈电机定子电流调节器的设计问题转变为求解H∞标准问题;考虑电机转速变化、参数不确定、定转子耦合等因素的影响,构造合理的权函数,并基于回路成形法对H∞电流调节器进行优化设计。本发明在定子静止坐标系下,通过采集双馈电机定子电压、电流和转速计算得到定子电流给定值,并与实际值进行比较后送入H∞电流调节器;H∞电流调节器输出信号为转子电压控制信号,由转子侧PWM逆变器实现并施加到转子绕组上,从而起到控制双馈电机定子电流以及输出功率的目的。

Description

基于回路成形的双馈电机定子电流鲁棒控制方法
技术领域
本发明涉及一种风力发电技术领域的控制方法,具体是一种基于回路成形的双馈电机定子电流鲁棒控制方法。
技术背景
随着全球能源危机以及气候环境问题的不断发展,可再生能源发电技术日益受到世界各国的重视。目前,风力发电技术较为成熟,并网风电装机容量上升较快,风电机组也正向着大型化、海上化的方向发展,这些因素在一定程度上要求风力发电机具有更好的控制性能。在各种风力发电技术中,双馈型感应发电机以其所需变频器容量小、系统效率高以及功率因数可调等优点,成为了当前变速恒频风力发电机的主流机型。双馈电机现有控制方法大致可以分为基于磁场/电压定向的矢量控制,以及直接功率/转矩控制;一般都是基于旋转坐标变换与PI较为(比例积分)调节器,通过控制双馈电机交流侧电流间接控制定子电流以及输出功率。
现有控制方法在一定程度上解决了双馈电机的有关控制问题。但PI调节器的参数设计往往依赖于工程调试经验,实际中由于电机参数与外界扰动的并不确定性,造成现有控制器参数设计存在一定盲目性。为了解决现有控制器鲁棒性差、参数整定繁琐等不足,专利CN101388637B提出一种带前馈补偿的双馈风力发电机鲁棒控制方法,利用加权函数与回路成形方法设计具有良好鲁棒性的双馈电机转子电流控制器。该方法的不足在于:1、在同步坐标系下进行控制,须附加前馈补偿,降低了鲁棒控制器的改善效果;2、通过转子电流间接控制定子电流,增加了参数不确定性的干扰。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提出一种基于回路成形的双馈电机定子电流鲁棒控制方法。本发明可以设计满足控制要求的H∞电流调节器,并可以直接对定子电流进行控制,提高双馈电机控制系统的鲁棒性能,且便于实现。
本发明的目的是通过以下技术方案来实现的。
本发明的一种基于回路成形的双馈电机定子电流鲁棒控制方法,包括以下步骤:
(1)在定子两相静止αβ坐标系下,建立双馈电机的数学模型,该模型的输入为定子电压us_αβ、转子电压ur_αβ、转子转速ωr,输出为定子电流is_αβ
(2)根据控制目标、参数摄动以及外界干扰构造一组权函数,包括性能权函数Wp、输出权函数Wy和控制权函数Wu,其中性能权函数用于满足定子电流无差跟踪调节的要求,输出权函数用于增强定子电流控制的鲁棒性,控制权函数用于限制电流调节器的输出幅值,权函数与双馈电机的数学模型一起构成广义控制对象模型P,广义控制对象模型的输入为定子电流参考信号电网电压us_αβ、转子转速ωr与转子电压ur_αβ,输出为性能加权zp、输出加权zy、控制加权zu与定子电流误差△is_αβ
(3)利用回路成形法,迭代地求得一个2输入2输出的参数化H∞电流调节器,使得广义控制对象模型P与H∞电流调节器构成的闭环回路具有最小的H∞范数;
(4)测量得到双馈电机的三相定子电压us_abc、三相定子电流is_abc,并转换到定子两相静止αβ坐标系下得到us_αβ和is_αβ,计算得到定子有功功率Ps和无功功率Qs,测量得到转子角速度ωr与位置信号θr
(5)将有功功率Ps与其参考信号比较后得到有功功率误差信号△Ps,将无功功率Qs与其参考小信号比较后得到无功功率误差信号△Qs
(6)将功率误差信号△Ps和△Qs通过比例积分调节,得到dq轴同步坐标系下的定子电流参考信号
(7)利用定子电压us_αβ将dq轴同步坐标系下的定子电流参考信号转换到静止坐标系下,得到静止坐标系下的定子电流参考信号
(8)将静止坐标系下的定子电流is_αβ与参考信号进行比较,得到电流误差信号△is_αβ
(9)将电流误差信号△is_αβ送入H∞电流调节器,得到定子两相静止坐标系下的转子电压参考信号
(10)利用转子位置信号θr,将转子电压参考信号转换到转子三相静止坐标系下,得到并由转子侧PWM逆变器施加到转子绕组,从而实现对双馈电机定子电流及输出功率的控制。
本发明的有益效果是:本发明基于回路成形的双馈电机定子电流鲁棒控制方法,在定子静止坐标系下,建立双馈电机数学模型与权函数,利用回路成形法设计满足控制要求的H∞电流调节器,能够在定子静止坐标系下对双馈电机的定子电流进行实时控制,并实现定子输出功率因数的连续可调;通过构造合理的性能权函数,实现了定子电流的无差跟踪控制;通过构造合理的输出权函数,能够针对电机参数和外部扰动的不确定性影响,提高定子电流控制的鲁棒稳定性。本发明提高了双馈电机控制系统的鲁棒性能,且便于实现。
附图说明
图1是本发明的H∞电流调节器设计原理图;
图2是本发明的双馈电机控制原理图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明。
本发明一种基于回路成形的双馈电机定子电流鲁棒控制方法,主要包含了两个部分。首先在在定子静止坐标系下建立双馈电机数学模型及一组权函数,并利用回路整形方法得到H∞电流调节器。其次,采用双闭环的矢量控制策略,通过H∞电流调节器直接对电机定子电流进行无差跟踪调节,从而控制电机输出功率。
图1是本发明的H∞电流调节器设计原理图,图中,方框1表示控制对象数学模型,方框2表示控制权函数Wu,方框3表示输出权函数Wy,方框4表示性能权函数Wp,方框5表示H∞电流调节器。
如图1所示,本发明的H∞电流调节器设计方法是,通过迭代寻找得到一个2输入2输出的参数化的H∞电流调节器,使得从输入向量w到输出向量z的闭环系统的H∞范数最小。本发明所述的H∞电流调节器设计过程如下:
1、在定子静止坐标系下,建立双馈电机的数学模型,数学模型的输入为定子电压us_αβ,转子电压ur_αβ,转速ωr;输出为定子电流is_αβ;中间变量有定子磁链ψs_αβ,转子磁链ψr_αβ;模型参数有定子绕组自感Ls,转子绕组自感Lr,定转子互感Lm,定子绕组电阻Rs,转子绕组电阻Rr。采用国际单位制,电机模型的数学表达式为:
dψ s _ αβ dt = - R s i s _ αβ + u s _ αβ dψ r _ αβ dt = - R r i r _ αβ + u r _ αβ + ω r ψ r _ αβ - - - ( 1 )
在公式(1)中,
ψ s _ αβ = L s i s _ αβ + L m i r _ αβ ψ r _ αβ = L m i s _ αβ + L r i r _ αβ - - - ( 2 )
公式(1)(2)可以简化地表示为:
is_αβ=f(us_αβ,ur_αβr)     (3)
2、将定子电流is_αβ与参考信号比较后得到误差信号△is_αβ,并与性能权函数Wp相乘得到广义输出zp;将定子电流is_αβ与输出权函数Wy相乘得到广义输出zp;将控制信号ur_αβ与控制权函数Wu相乘得到广义输出zu。图1所示的虚线框中从w、uc_αβ到z、y的部分即为广义控制对象P。其中性能权函数用于满足定子电流无差跟踪调节的要求,输出权函数用于增强定子电流控制的鲁棒性,控制权函数用于限制电流调节器的输出幅值。
性能权函数的数学表达式为:
W p = k ω s 2 s 2 + 2 ζ ω s s + ω s 2 - - - ( 4 )
其中,ωs为电网电压角频率;k为增益系数,用于控制跟踪精度;ζ为防止无穷大增益的阻尼系数,一般设为一较小常数,如0.001。
输出权函数的表达式为:
W y = ( s + ω b ) n ( s M 1 / n + ω b A 1 / n ) n - - - ( 5 )
其中,n为权函数阶数,一般选择1~2阶即可;A为低频段增益上界,典型值为2;M为高频段增益上界,典型值为0.001;ωb为穿越频率,可设置在100~500Hz。
控制权函数的表达式为:
Wu=ku     (6)
其中,ku为增益系数,当使用标幺参数时,ku可设为1。
H∞电流调节器的设计问题可以表述为寻找一个合适的电流调节器,使得广义控制对象P的H∞范数最小化,即寻找一个足够小的正数γmin,使得:
W p S W u KS W y T &infin; &le; &gamma; min < 1 - - - ( 7 )
其中S称为灵敏度函数,是从w的y的闭环传递函数。T称为互补灵敏度函数,是从w到z的闭环传递函数。γmin为H∞范数的峰值。
3、基于回路成形,迭代地求解公式(7)所表示的H∞问题,得到所需的2输入2输出的H∞电流调节器。
图2是本发明的定子电流矢量控制原理图,其中包括H∞电流调节器5、两相/三相静止坐标系变换模块6、转子侧PWM逆变器7、直流母线电容8、直流电压传感器9、双馈电机机械转轴10、光电码盘11、转子绕组12、定子绕组13、三相交流电流传感器14、三相交流电压传感器15、三相/两相静止坐标变换模块16、瞬时功率计算模块17、有功功率调节器18、无功功率调节器19、两相旋转/两相静止坐标变换模块20。
如图2所示,本发明提出的定子电流矢量控制策略的实施过程如下:
1、利用光电码盘11测算得到双馈电机转子角速度ωr与转子位置角θr。利用交流电流传感器14、交流电压传感器15测量得到双馈电机定子三相电压us_abc和三相电流is_abc,并通过模块16在静止坐标系下进行三相/两相变换,得到定子电流矢量is_αβ和电压矢量us_αβ,变换公式为:
i s _ &alpha; i s _ &beta; = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i s _ a i s _ b i s _ c - - - ( 8 )
u s _ &alpha; u s _ &beta; = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u s _ a u s _ b u s _ c - - - ( 9 )
将电压矢量us_αβ和电流矢量is_αβ送入瞬时功率计算模块17计算电机定子输出有功功率Ps和无功功率Qs,计算公式为:
P s Q s = u s _ &alpha; u s _ &beta; u s _ &beta; - u s _ &alpha; i s _ &alpha; i s _ &beta; - - - ( 10 )
2、将有功功率Ps与其参考信号比较后得到有功功率误差信号△Ps
&Delta;P s = P s * - P s - - - ( 11 )
将无功功率Qs与其参考小信号比较后得到无功功率误差信号△Qs
&Delta;Q s = Q s * - Q s - - - ( 12 )
3、将功率误差信号△Ps和△Qs通过PI(比例积分)调节器,得到dq轴同步坐标系下的定子电流参考信号PI调节器的数学表达式为:
i s _ d * = ( k p _ d + k i _ d s ) &Delta;P s - - - ( 13 )
i s _ q * = ( k p _ q + k i _ q s ) &Delta;Q s - - - ( 14 )
其中,kp_d、kp_q分别为dq轴比例系数,ki_d、ki_q分别为dq轴积分系数;
4、利用定子电压us_αβ将dq轴同步坐标系下的定子电流参考信号转换到静止坐标系下,得到静止坐标系下的定子电流参考信号计算公式为:
i s _ &alpha; * i s _ &beta; * = 1 u s _ &alpha; 2 + u s _ &beta; 2 u s _ &alpha; - u s _ &beta; u s _ &beta; u s _ &alpha; i s _ d * i s _ q * - - - ( 15 )
5、将静止坐标系下的定子电流is_αβ与参考信号进行比较,得到电流误差信号△is_αβ
&Delta;i s _ &alpha;&beta; = i s _ &alpha;&beta; * - i s _ &alpha;&beta; - - - ( 16 )
6、将电流误差信号△is_αβ送入H∞电流调节器5,得到定子两相静止坐标系下的转子电压参考信号
u r _ &alpha;&beta; * = K&Delta; i s _ &alpha;&beta; - - - ( 17 )
7、利用转子位置信号θr,将转子电压参考信号转换到转子三相静止坐标系下,得到计算公式为:
u r _ a * u r _ b * u r _ c * = 2 3 cos &theta; r - sin &theta; r cos ( &theta; r - 2 &pi; / 3 ) - sin ( &theta; r - 2 &pi; / 3 ) cos ( &theta; r - 4 &pi; / 3 ) - sin ( &theta; r - 4 &pi; / 3 ) u r _ &alpha; * u r _ &beta; * - - - ( 18 )
最后由转子侧PWM逆变器实现参考电压,并施加到转子绕组上,从而实现对双馈电机定子电流及输出功率的控制。

Claims (8)

1.一种基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于包括以下步骤:
(1)在定子两相静止αβ坐标系下,建立双馈电机的数学模型,该模型的输入为定子电压us_αβ、转子电压ur_αβ、转子转速ωr,输出为定子电流is_αβ
(2)根据控制目标、参数摄动以及外界干扰构造一组权函数,包括性能权函数Wp、输出权函数Wy和控制权函数Wu,其中性能权函数用于满足定子电流无差跟踪调节的要求,输出权函数用于增强定子电流控制的鲁棒性,控制权函数用于限制电流调节器的输出幅值,权函数与双馈电机的数学模型一起构成广义控制对象模型P,广义控制对象模型的输入为定子电流参考信号电网电压us_αβ、转子转速ωr与转子电压ur_αβ,输出为性能加权zp、输出加权zy、控制加权zu与定子电流误差△is_αβ
(3)利用回路成形法,迭代地求得一个2输入2输出的参数化H∞电流调节器,使得广义控制对象模型P与H∞电流调节器构成的闭环回路具有最小的H∞范数;
(4)测量得到双馈电机的三相定子电压us_abc、三相定子电流is_abc,并转换到定子两相静止αβ坐标系下得到us_αβ和is_αβ,计算得到定子有功功率Ps和无功功率Qs,测量得到转子角速度ωr与位置信号θr
(5)将有功功率Ps与其参考信号比较后得到有功功率误差信号△Ps,将无功功率Qs与其参考小信号比较后得到无功功率误差信号△Qs
(6)将功率误差信号△Ps和△Qs通过比例积分调节,得到dq轴同步坐标系下的定子电流参考信号
(7)利用定子电压us_αβ将dq轴同步坐标系下的定子电流参考信号转换到静止坐标系下,得到静止坐标系下的定子电流参考信号
(8)将静止坐标系下的定子电流is_αβ与参考信号进行比较,得到电流误差信号△is_αβ
(9)将电流误差信号△is_αβ送入H∞电流调节器,得到定子两相静止坐标系下的转子电压参考信号
(10)利用转子位置信号θr,将转子电压参考信号转换到转子三相静止坐标系下,得到并由转子侧PWM逆变器施加到转子绕组,从而实现对双馈电机定子电流及输出功率的控制。
2.根据权利要求1所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述步骤(2)所构造的性能权函数的数学表达式为:
W p = k&omega; s 2 s 2 + 2 &zeta; &omega; s s + &omega; s 2
其中,ωs为电网电压角频率;k为增益系数,用于控制跟踪精度;ζ为防止无穷大增益的阻尼系数。
3.根据权利要求1所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述输出权函数的表达式为:
W y = ( s + &omega; b ) n ( sM 1 / n + &omega; b A 1 / n ) n - - - ( 5 )
其中,n为权函数阶数,A为低频段增益上界,M为高频段增益上界,ωb为穿越频率。
4.根据权利要求1所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述控制权函数的表达式为:
Wu=ku                      (6)
其中,ku为增益系数。
5.根据权利要求1所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述H∞电流调节器在两相静止坐标系下对定子电流参考信号进行无差跟踪调节,从而实现对双馈电机输出功率的控制,且对电机参数变化以及外部扰动不敏感。
6.根据权利要求2所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述防止无穷大增益的阻尼系数ζ设为0.001。
7.根据权利要求3所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述权函数阶数n为1~2阶;低频段增益上界A典型值为2;高频段增益上界M典型值为0.001;穿越频率ωb设置在100~500Hz。
8.根据权利要求4所述的基于回路成形的双馈电机定子电流鲁棒控制方法,其特征在于上述增益系数ku当使用标幺参数时,ku设为1。
CN201410567661.2A 2014-10-22 2014-10-22 基于回路成形的双馈电机定子电流鲁棒控制方法 Active CN104333283B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410567661.2A CN104333283B (zh) 2014-10-22 2014-10-22 基于回路成形的双馈电机定子电流鲁棒控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410567661.2A CN104333283B (zh) 2014-10-22 2014-10-22 基于回路成形的双馈电机定子电流鲁棒控制方法

Publications (2)

Publication Number Publication Date
CN104333283A true CN104333283A (zh) 2015-02-04
CN104333283B CN104333283B (zh) 2017-02-08

Family

ID=52407954

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410567661.2A Active CN104333283B (zh) 2014-10-22 2014-10-22 基于回路成形的双馈电机定子电流鲁棒控制方法

Country Status (1)

Country Link
CN (1) CN104333283B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106982021A (zh) * 2017-03-31 2017-07-25 中国电力科学研究院 一种并网型双馈感应发电机定子电流控制方法及装置
CN107666220A (zh) * 2017-11-16 2018-02-06 江苏大中电机股份有限公司 一种超高效交流永磁同步电机的鲁棒性优化方法
CN112564565A (zh) * 2020-12-07 2021-03-26 北京航空航天大学 一种永磁同步电机伺服系统的全回路抗干扰鲁棒控制方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19919752C5 (de) * 1999-04-29 2010-12-16 Sew-Eurodrive Gmbh & Co. Kg Verfahren zum Zuschalten eines Umrichters an einen Asynchronmotor
CN102055205A (zh) * 2010-10-26 2011-05-11 天津理工大学 一种基于电压不平衡数学模型的无刷双馈电机网侧变换器控制策略
CN103414212B (zh) * 2013-08-09 2015-10-21 湖南大学 考虑双馈式风电机接入的配网系统的模拟方法
CN104065100B (zh) * 2014-06-13 2017-05-17 云南电力试验研究院(集团)有限公司电力研究院 一种提高不对称故障穿越能力的双馈风力机变流器控制装置
CN104065320B (zh) * 2014-06-17 2017-08-01 合肥学院 双馈风力发电机组无速度传感器控制方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106982021A (zh) * 2017-03-31 2017-07-25 中国电力科学研究院 一种并网型双馈感应发电机定子电流控制方法及装置
CN106982021B (zh) * 2017-03-31 2021-11-19 中国电力科学研究院 一种并网型双馈感应发电机定子电流控制方法及装置
CN107666220A (zh) * 2017-11-16 2018-02-06 江苏大中电机股份有限公司 一种超高效交流永磁同步电机的鲁棒性优化方法
CN112564565A (zh) * 2020-12-07 2021-03-26 北京航空航天大学 一种永磁同步电机伺服系统的全回路抗干扰鲁棒控制方法
CN112564565B (zh) * 2020-12-07 2022-07-05 北京航空航天大学 一种永磁同步电机伺服系统的全回路抗干扰鲁棒控制方法

Also Published As

Publication number Publication date
CN104333283B (zh) 2017-02-08

Similar Documents

Publication Publication Date Title
CN101478283B (zh) 不平衡电网电压下双馈异步风力发电机转子侧逆变器控制方法
CN108429284B (zh) 基于电压源输出双馈风电机组的谐波电压比例前馈补偿方法
CN101521481B (zh) 一种双馈异步风力发电系统不对称协调直接功率控制方法
CN101272119B (zh) 抑制双馈式风力发电机组定子电流不平衡及畸变的方法
CN105896600B (zh) 一种永磁同步直驱风力发电机并网系统控制方法
CN101557190A (zh) 不平衡电网电压下双馈异步风力发电机网侧整流器控制方法
CN107565599B (zh) 一种基于vsg的风电并网半实物仿真系统
CN104579060B (zh) 笼型转子无刷双馈风力发电机的间接功率控制方法
CN103762921B (zh) 一种基于粒子群算法的不平衡电网下dfig的多目标控制方法
CN104079226A (zh) 一种同步坐标系下无锁相环的dfig控制方法
CN102005782A (zh) 双馈风力发电机的并网控制方法
CN103117700A (zh) 一种不平衡电网下基于谐振反馈的dfig控制方法
CN107124126A (zh) 一种双馈感应发电机的无锁相环电流控制方法及装置
CN112542855A (zh) 一种双馈风力发电系统相量模型建模及仿真方法
CN104333283A (zh) 基于回路成形的双馈电机定子电流鲁棒控制方法
CN102790515B (zh) 基于可变开关频率的风力发电双馈变流器svpwm的控制方法
CN102355192B (zh) 双馈风力发电机无功功率的控制方法
CN102522942B (zh) 双馈风力发电机励磁控制方法
Abdelbaset et al. Wind driven doubly fed induction generator
CN101272118A (zh) 并网型双馈式风力发电机转子位置检测方法
Zhen et al. Study on control strategy of maximum power capture for DFIG in wind turbine system
He et al. Grey prediction pi control of direct drive permanent magnet synchronous wind turbine
CN105552951A (zh) 一种基于重复滑模的dfig系统控制方法
CN102332861B (zh) 双馈风力发电机有功功率的控制方法
CN105515040A (zh) 一种基于滑模+重复的dfig控制方法

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