CN106940429A - 基于发电机控制装置的低频振荡源定位及识别方法 - Google Patents

基于发电机控制装置的低频振荡源定位及识别方法 Download PDF

Info

Publication number
CN106940429A
CN106940429A CN201710103709.8A CN201710103709A CN106940429A CN 106940429 A CN106940429 A CN 106940429A CN 201710103709 A CN201710103709 A CN 201710103709A CN 106940429 A CN106940429 A CN 106940429A
Authority
CN
China
Prior art keywords
low
oscillation
generator
frequency oscillation
positioning
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
CN201710103709.8A
Other languages
English (en)
Other versions
CN106940429B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201710103709.8A priority Critical patent/CN106940429B/zh
Publication of CN106940429A publication Critical patent/CN106940429A/zh
Application granted granted Critical
Publication of CN106940429B publication Critical patent/CN106940429B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/34Testing dynamo-electric machines
    • G01R31/343Testing dynamo-electric machines in operation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/34Testing dynamo-electric machines
    • G01R31/346Testing of armature or field windings

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种基于发电机控制装置的低频振荡源定位及识别方法,包括以下步骤:S1:当电网发生低频振荡时,根据割集能量法锁定振荡区域内的机组群,并且获取发电机控制装置的量测电气数据;所述发电机控制装置包括励磁系统和调速系统;S2:采用经验模态分解算法提取已获取的量测电气数据在主导振荡模式下的电气变化量,并进行归一化处理;S3:分别计算励磁系统和调速系统在发生低频振荡时的输出能量变化率,进一步基于经验模态分解算法对输出能量变化率进行分解过滤并利用最小二乘法进行曲线拟合获得量化指标;S4:根据指标进行低频振荡源的定位以及识别,判断发电机的励磁系统及调速系统是否为振荡源并指出该振荡源的类型。

Description

基于发电机控制装置的低频振荡源定位及识别方法
技术领域
本发明涉及电力系统稳定分析与控制技术,特别是涉及基于发电机控制装置的低频振荡源定位及识别方法。
背景技术
随着现代电力系统互联规模的日益扩大,系统的动态过程越来越复杂,低频振荡已经成为影响电力系统安全稳定运行的重要问题。引起电力系统低频振荡的因素主要包括弱(负)阻尼理论以及共振机理强迫振荡理论这两类。其中,弱阻尼振荡与强迫功率振荡具有非常相似的表现形式,但其发生机理和抑制措施有本质的区别。如何能够准确的定位振荡源,并且识别振荡源类型对抑制振荡具有重要的意义。
对同步发电机而言,原动机及其调速系统是发电机的输入能量源,能够引发强迫功率振荡的可能性和危害性非常大;而励磁系统不能控制能量注入,不能单独作为强迫扰动源,其与调速系统特性的好坏能够影响系统的阻尼。
发明内容
发明目的:本发明的目的是提供一种能够解决现有技术中存在的缺陷的基于发电机控制装置的低频振荡源定位及识别方法。
技术方案:为达到此目的,本发明采用以下技术方案:
本发明所述的基于发电机控制装置的低频振荡源定位及识别方法,该方法包括以下步骤:
S1:当电网发生低频振荡时,根据割集能量法锁定振荡区域内的机组群,并且获取发电机控制装置的量测电气数据;所述发电机控制装置包括励磁系统和调速系统;
S2:采用经验模态分解算法提取已获取的量测电气数据在主导振荡模式下的电气变化量,并进行归一化处理;
S3:分别计算励磁系统和调速系统在发生低频振荡时的输出能量变化率,进一步基于经验模态分解算法对输出能量变化率进行分解过滤并利用最小二乘法进行曲线拟合获得量化指标;
S4:根据指标进行低频振荡源的定位以及识别,判断发电机的励磁系统及调速系统是否为振荡源并指出该振荡源的类型。
进一步,所述步骤S1中的量测电气数据包括励磁系统的励磁电动势Ef和励磁电流if、原动机输出机械功率Pm和发电机转速ω。
进一步,所述步骤S2包括以下步骤:
S2.1:采用经验模态分解算法过滤掉低频振荡时量测电气数据的干扰信号和衰减信号,提取出主导振荡模式下的电气变化量ΔEf、Δif、ΔPm和Δω;
S2.2:求取序列ΔEf、Δif、ΔPm和Δω的上下包络线,进一步求出其幅值,分别进行归一化处理,得到ΔE1 f,Δi1 f,ΔP1 m,Δω1
进一步,所述步骤S3包括以下步骤:
S3.1:根据式(1)计算励磁系统在发生低频振荡时的输出能量变化率根据式(2)计算调速系统在发生低频振荡时的输出能量变化率Em D
式(1)中,ΔE1 f为为归一化后的发电机励磁电动势的主导振荡模式下的变化量,Δi1 f为归一化后的发电机励磁电流的主导振荡模式下的变化量;
Em D=ΔP1 m*Δω1 (2)
式(2)中,ΔP1 m为归一化后的发电机机械功率的主导振荡模式下的变化量,Δω1为归一化后的发电机转子转速的主导振荡模式下的变化量;
S3.2:对和Em D进行EMD分解,减去经验模态分量,剩下恒定分量;对恒定分量进行最小二乘曲线拟合,获得振荡源定位指标及识别指标Em D1D
进一步,所述步骤S4中的判断方法包括以下步骤:
S4.1:判断是否成立:如果是,则进行步骤S4.2;否则,进行步骤S4.3;
S4.2:判断Em D1D>0是否成立:如果是,则进行步骤S4.4;否则,则进行步骤S4.5;
S4.3:判断Em D1D>0是否成立:如果是,则进行步骤S4.6;否则,则进行步骤S4.7;
S4.4:判断Em D1D≈1/2是否成立:如果是,则判定该机组励磁系统提供负阻尼且调速系统存在强迫扰动源;否则,则判定该机组调速系统和励磁系统同时提供负阻尼;
S4.5:判定该机组励磁系统提供负阻尼;
S4.6:判断Em D1D≈1/2是否成立:如果是,则判定该机组调速系统存在强迫扰动源;否则,则判定该机组调速系统提供负阻尼;
S4.7:判定该机组不是振荡源。
有益效果:本发明从发电机引发低频振荡的发生机理的角度出发,深入到发电机的励磁系统和调速系统进行分析,能够准确地判断出该发电机的控制装置是否为振荡源并且可以识别该振荡源的类型;同时本发明利用经验模态分解方法的滤波功能将振荡的主导分量提取出来,且进一步制定振荡源定位及识别的量化指标,既能提高准确性,又能实现自动功能。此外,本发明所需数据均来自发电机组,不需要电网中完善的PMU配置,可实现分散化监测。
附图说明
图1为本发明具体实施方式的方法流程图;
图2为本发明具体实施方式的简化原动机调节系统模型图;
图3为本发明具体实施方式的四机两区系统接线图。
具体实施方式
下面结合附图和具体实施方式对本发明的技术方案作进一步的介绍。
本具体实施方式以四机两区系统来详细说明。
如图2所示,发电机的励磁系统采用简化的汽轮机调节系统模型,其中,Ka为放大倍数(调差系数的倒数),Tg为调速器液压系统时间常数;TCH为汽轮机高压缸蒸汽容积时间常数。
四机两区系统如图3所示,根据小干扰稳定计算分析,该系统存在区域内弱阻尼振荡模式,振荡频率为1.145Hz,G1、G2相关因子分别为0.15和0.32;其中G1、G2调速系统参数设置一致,分界频率为0.411Hz。首先仿真系统弱阻尼振荡,在1s时对母线6施加一个瞬时三相短路故障。接着调整四台发电机的有功出力,使系统呈现正阻尼。根据小干扰稳定计算分析,该系统存在区域间振荡模式,振荡频率为0.708Hz。在机组G2的原动机侧添加频率为0.708Hz,大小为原机组输出有功功率10%的波动性功率,激发强迫功率振荡。其中,G1调速系统的分界频率为0.411Hz,G3、G4调速系统分界频率一致,为0.919Hz。
本具体实施方式公开了一种基于发电机控制装置的低频振荡源定位及识别方法,如图1所示,该方法包括以下步骤:
S1:当电网发生低频振荡时,根据割集能量法锁定振荡区域内的机组群G1-G4,并且获取发电机控制装置的量测电气数据;其中,发电机控制装置包括励磁系统和调速系统;量测电气数据包括励磁系统的励磁电动势Ef和励磁电流if、原动机输出机械功率Pm和发电机转速ω;
S2:采用经验模态分解算法提取已获取的量测电气数据在主导振荡模式下的电气变化量,并进行归一化处理;
S3:分别计算励磁系统和调速系统在发生低频振荡时的输出能量变化率,进一步基于经验模态分解算法对输出能量变化率进行分解过滤并利用最小二乘法进行曲线拟合获得量化指标;
S4:根据指标进行低频振荡源的定位以及识别,判断发电机的励磁系统及调速系统是否为振荡源并指出该振荡源的类型。
步骤S2包括以下步骤:
S2.1:采用经验模态分解算法过滤掉低频振荡时量测电气数据的干扰信号和衰减信号,提取出主导振荡模式下的电气变化量ΔEf、Δif、ΔPm和Δω;
S2.2:求取序列ΔEf、Δif、ΔPm和Δω的上下包络线,进一步求出其幅值,分别进行归一化处理,得到ΔE1 f,Δi1 f,ΔP1 m,Δω1
步骤S3包括以下步骤:
S3.1:根据式(1)计算励磁系统在发生低频振荡时的输出能量变化率根据式(2)计算调速系统在发生低频振荡时的输出能量变化率Em D
式(1)中,ΔE1 f为归一化后的发电机励磁电动势的主导振荡模式下的变化量,Δi1 f为归一化后的发电机励磁电流的主导振荡模式下的变化量;
Em D=ΔP1 m*Δω1 (2)
式(2)中,ΔP1 m为归一化后的发电机机械功率的主导振荡模式下的变化量,Δω1为归一化后的发电机转子转速的主导振荡模式下的变化量;
S3.2:对和Em D进行EMD分解,减去经验模态分量,剩下恒定分量;对恒定分量进行最小二乘曲线拟合,获得振荡源定位指标及识别指标Em D1D,分别如表1和表2所示。
表1振荡源识别分析结果(弱阻尼振荡)
表2振荡源识别分析结果(强迫功率振荡)
步骤S4中的判断方法包括以下步骤:
S4.1:判断是否成立:如果是,则进行步骤S4.2;否则,进行步骤S4.3;
S4.2:判断Em D1D>0是否成立:如果是,则进行步骤S4.4;否则,则进行步骤S4.5;
S4.3:判断Em D1D>0是否成立:如果是,则进行步骤S4.6;否则,则进行步骤S4.7;
S4.4:判断Em D1D≈1/2是否成立:如果是,则判定该机组励磁系统提供负阻尼且调速系统存在强迫扰动源;否则,则判定该机组调速系统和励磁系统同时提供负阻尼;
S4.5:判定该机组励磁系统提供负阻尼;
S4.6:判断Em D1D≈1/2是否成立:如果是,则判定该机组调速系统存在强迫扰动源;否则,则判定该机组调速系统提供负阻尼;
S4.7:判定该机组不是振荡源。
由表1可知,由于G2负荷较重,其励磁系统提供负阻尼;且G1、G2调速系统分界频率低于低频振荡频率,皆提供负阻尼,导致系统呈弱阻尼。G3、G4在该低频振荡模式下的参与度较小,提供的阻尼为正。振荡源识别分析结果与实际情况一致。
由表2可知,G2指标Em D1D为0.4842,接近于0.5,故强迫扰动源位于G2原动机侧;且G1调速系统分界频率低于低频振荡频率,提供负阻尼;而G3、G4调速系统分界频率高于低频振荡频率,提供阻尼为正。振荡源识别分析结果与实际情况一致。

Claims (5)

1.基于发电机控制装置的低频振荡源定位及识别方法,其特征在于:该方法包括以下步骤:
S1:当电网发生低频振荡时,根据割集能量法锁定振荡区域内的机组群,并且获取发电机控制装置的量测电气数据;所述发电机控制装置包括励磁系统和调速系统;
S2:采用经验模态分解算法提取已获取的量测电气数据在主导振荡模式下的电气变化量,并进行归一化处理;
S3:分别计算励磁系统和调速系统在发生低频振荡时的输出能量变化率,进一步基于经验模态分解算法对输出能量变化率进行分解过滤并利用最小二乘法进行曲线拟合获得量化指标;
S4:根据指标进行低频振荡源的定位以及识别,判断发电机的励磁系统及调速系统是否为振荡源并指出该振荡源的类型。
2.根据权利要求1所述的基于发电机控制装置的低频振荡源定位及识别方法,其特征在于:所述步骤S1中的量测电气数据包括励磁系统的励磁电动势Ef和励磁电流if、原动机输出机械功率Pm和发电机转速ω。
3.根据权利要求1所述的基于发电机控制装置的低频振荡源定位及识别方法,其特征在于:所述步骤S2包括以下步骤:
S2.1:采用经验模态分解算法过滤掉低频振荡时量测电气数据的干扰信号和衰减信号,提取出主导振荡模式下的电气变化量ΔEf、Δif、ΔPm和Δω;
S2.2:求取序列ΔEf、Δif、ΔPm和Δω的上下包络线,进一步求出其幅值,分别进行归一化处理,得到ΔE1 f,Δi1 f,ΔP1 m,Δω1
4.根据权利要求1所述的基于发电机控制装置的低频振荡源定位及识别方法,其特征在于:所述步骤S3包括以下步骤:
S3.1:根据式(1)计算励磁系统在发生低频振荡时的输出能量变化率根据式(2)计算调速系统在发生低频振荡时的输出能量变化率Em D
E f D = ΔE 1 f * Δi 1 f - - - ( 1 )
式(1)中,ΔE1 f为归一化后的发电机励磁电动势的主导振荡模式下的变化量,Δi1 f为归一化后的发电机励磁电流的主导振荡模式下的变化量;
Em D=ΔP1 m*Δω1 (2)
式(2)中,ΔP1 m为归一化后的发电机机械功率的主导振荡模式下的变化量,Δω1为归一化后的发电机转子转速的主导振荡模式下的变化量;
S3.2:对和Em D进行EMD分解,减去经验模态分量,剩下恒定分量;对恒定分量进行最小二乘曲线拟合,获得振荡源定位指标及识别指标Em D1D
5.根据权利要求4所述的基于发电机控制装置的低频振荡源定位及识别方法,其特征在于:所述步骤S4中的判断方法包括以下步骤:
S4.1:判断是否成立:如果是,则进行步骤S4.2;否则,进行步骤S4.3;
S4.2:判断Em D1D>0是否成立:如果是,则进行步骤S4.4;否则,则进行步骤S4.5;
S4.3:判断Em D1D>0是否成立:如果是,则进行步骤S4.6;否则,则进行步骤S4.7;
S4.4:判断Em D1D≈1/2是否成立:如果是,则判定该机组励磁系统提供负阻尼且调速系统存在强迫扰动源;否则,则判定该机组调速系统和励磁系统同时提供负阻尼;
S4.5:判定该机组励磁系统提供负阻尼;
S4.6:判断Em D1D≈1/2是否成立:如果是,则判定该机组调速系统存在强迫扰动源;否则,则判定该机组调速系统提供负阻尼;
S4.7:判定该机组不是振荡源。
CN201710103709.8A 2017-02-24 2017-02-24 基于发电机控制装置的低频振荡源定位及识别方法 Active CN106940429B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710103709.8A CN106940429B (zh) 2017-02-24 2017-02-24 基于发电机控制装置的低频振荡源定位及识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710103709.8A CN106940429B (zh) 2017-02-24 2017-02-24 基于发电机控制装置的低频振荡源定位及识别方法

Publications (2)

Publication Number Publication Date
CN106940429A true CN106940429A (zh) 2017-07-11
CN106940429B CN106940429B (zh) 2019-03-12

Family

ID=59469059

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710103709.8A Active CN106940429B (zh) 2017-02-24 2017-02-24 基于发电机控制装置的低频振荡源定位及识别方法

Country Status (1)

Country Link
CN (1) CN106940429B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109830972A (zh) * 2019-02-15 2019-05-31 湖南大学 一种新能源场站振荡源快速识别系统及方法
CN110323741A (zh) * 2019-07-10 2019-10-11 国网福建省电力有限公司 基于录波曲线的电力系统低频振荡分析与抑制方法
CN110879355A (zh) * 2019-12-03 2020-03-13 国网湖南省电力有限公司 一种基于fft的发电机调速系统引发低频振荡辨识方法
CN110879354A (zh) * 2019-12-02 2020-03-13 国网湖南省电力有限公司 基于fft算法的发电机组扰动源定位方法
CN110879324A (zh) * 2019-12-02 2020-03-13 国网湖南省电力有限公司 一种励磁系统在功率振荡中阻尼特性的辨别方法
CN113471960A (zh) * 2021-05-31 2021-10-01 国网江苏省电力有限公司 一种基于多源数据融合的低频振荡类型判断方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102983569A (zh) * 2012-12-12 2013-03-20 南京南瑞继保电气有限公司 一种电力系统低频振荡的振荡源定位方法
CN104007710A (zh) * 2013-02-25 2014-08-27 中国石油化工股份有限公司 化工振荡回路振荡现象的自动溯源方法
CN104330667A (zh) * 2014-10-29 2015-02-04 华中电网有限公司 一种用于判断调速系统对低频振荡阻尼影响的监测方法
CN104898649A (zh) * 2015-05-08 2015-09-09 东南大学 流程工业的厂级控制回路振荡源定位方法
CN104950843A (zh) * 2015-05-08 2015-09-30 东南大学 基于模型识别方法的流程工业厂级控制回路振荡源定位方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102983569A (zh) * 2012-12-12 2013-03-20 南京南瑞继保电气有限公司 一种电力系统低频振荡的振荡源定位方法
CN104007710A (zh) * 2013-02-25 2014-08-27 中国石油化工股份有限公司 化工振荡回路振荡现象的自动溯源方法
CN104330667A (zh) * 2014-10-29 2015-02-04 华中电网有限公司 一种用于判断调速系统对低频振荡阻尼影响的监测方法
CN104898649A (zh) * 2015-05-08 2015-09-09 东南大学 流程工业的厂级控制回路振荡源定位方法
CN104950843A (zh) * 2015-05-08 2015-09-30 东南大学 基于模型识别方法的流程工业厂级控制回路振荡源定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
蒋平 等: "基于功率及频率波动相位的强迫扰动源定位", 《电力自动化设备》 *
郑斌青 等: "基于WAMS的强迫振荡源在线定位方法", 《江苏电机工程》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109830972A (zh) * 2019-02-15 2019-05-31 湖南大学 一种新能源场站振荡源快速识别系统及方法
CN109830972B (zh) * 2019-02-15 2022-11-18 湖南大学 一种新能源场站振荡源快速识别系统及方法
CN110323741A (zh) * 2019-07-10 2019-10-11 国网福建省电力有限公司 基于录波曲线的电力系统低频振荡分析与抑制方法
CN110323741B (zh) * 2019-07-10 2022-05-13 国网福建省电力有限公司 基于录波曲线的电力系统低频振荡分析与抑制方法
CN110879354A (zh) * 2019-12-02 2020-03-13 国网湖南省电力有限公司 基于fft算法的发电机组扰动源定位方法
CN110879324A (zh) * 2019-12-02 2020-03-13 国网湖南省电力有限公司 一种励磁系统在功率振荡中阻尼特性的辨别方法
CN110879324B (zh) * 2019-12-02 2022-05-17 国网湖南省电力有限公司 一种励磁系统在功率振荡中阻尼特性的辨别方法
CN110879355A (zh) * 2019-12-03 2020-03-13 国网湖南省电力有限公司 一种基于fft的发电机调速系统引发低频振荡辨识方法
CN110879355B (zh) * 2019-12-03 2022-03-25 国网湖南省电力有限公司 一种基于fft的发电机调速系统引发低频振荡辨识方法
CN113471960A (zh) * 2021-05-31 2021-10-01 国网江苏省电力有限公司 一种基于多源数据融合的低频振荡类型判断方法和装置
CN113471960B (zh) * 2021-05-31 2024-02-23 国网江苏省电力有限公司 一种基于多源数据融合的低频振荡类型判断方法和装置

Also Published As

Publication number Publication date
CN106940429B (zh) 2019-03-12

Similar Documents

Publication Publication Date Title
CN106940429A (zh) 基于发电机控制装置的低频振荡源定位及识别方法
Delgado-Arredondo et al. Methodology for fault detection in induction motors via sound and vibration signals
Zaker et al. Improving synchronous generator parameters estimation using $ d-q $ axes tests and considering saturation effect
CN103036498A (zh) 一种基于pmu的同步发电机实用模型参数校核与辨识方法
CN109830972B (zh) 一种新能源场站振荡源快速识别系统及方法
CN109753689A (zh) 一种电力系统机电振荡模态特征参数在线辩识方法
CN105259414A (zh) 一种基于逆变器的电网阻抗在线检测方法
CN106526384A (zh) 一种用于大规模电力系统的振荡源定位方法
CN106058897B (zh) 一种基于相量的发电机强迫扰动源定位方法
US20150241869A1 (en) Numerical control device having function of calculating frequency characteristic of control loop
CN109605128A (zh) 一种基于功率谱熵差的铣削颤振在线检测方法
CN109283427A (zh) 含双馈风电机组的互联电网低频扰动源定位系统及方法
CN104392141A (zh) 电力系统低频振荡扰动源定位方法和装置
TW201545833A (zh) 切削加工參數估測裝置及其方法
Zhai et al. Influence of cutting parameters on force coefficients and stability in plunge milling
CN104101805B (zh) 一种基于振荡能量注入的励磁系统负阻尼检测方法
CN108535558B (zh) 一种电力稳定器的测试方法及系统
CN105808886B (zh) 一种主轴回转误差辨识方法
CN105158687B (zh) 一种发电机次同步振荡模态信号的测量设备和方法
CN105429157B (zh) 一种基于Prony分析的智能振荡模式识别方法
Thomsen et al. Online parameter identification methods for doubly fed induction generators
CN106066440B (zh) 一种pss负阻尼检测方法
CN110879355B (zh) 一种基于fft的发电机调速系统引发低频振荡辨识方法
CN107167733A (zh) 一种同步发电机励磁系统性能的基础评估数据的获取方法
CN105048920A (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