CN106772042A - 一种发电机原动系统阻尼在线检测方法 - Google Patents

一种发电机原动系统阻尼在线检测方法 Download PDF

Info

Publication number
CN106772042A
CN106772042A CN201611209173.XA CN201611209173A CN106772042A CN 106772042 A CN106772042 A CN 106772042A CN 201611209173 A CN201611209173 A CN 201611209173A CN 106772042 A CN106772042 A CN 106772042A
Authority
CN
China
Prior art keywords
generator
energy stream
delta
dissipation energy
excitation
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
CN201611209173.XA
Other languages
English (en)
Other versions
CN106772042B (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.)
Tsinghua University
China Southern Power Grid Co Ltd
Original Assignee
Tsinghua University
China Southern Power Grid 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 Tsinghua University, China Southern Power Grid Co Ltd filed Critical Tsinghua University
Priority to CN201611209173.XA priority Critical patent/CN106772042B/zh
Publication of CN106772042A publication Critical patent/CN106772042A/zh
Application granted granted Critical
Publication of CN106772042B publication Critical patent/CN106772042B/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

Abstract

本发明提供一种发电机原动系统阻尼在线测量方法,包括以下步骤:S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流;S2、计算发电机流入励磁系统的励磁耗散能量流;S3、根据电网流入发电机的耗散能量流和励磁耗散能量流之差计算出发电机流入原动系统的原动耗散能量流;S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合,并计算两者斜率比,得到原动系统的阻尼转矩系数。通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流,并在线计算发电机原动系统的阻尼转矩系数。

Description

一种发电机原动系统阻尼在线检测方法
技术领域
本发明涉及电力系统分析技术领域,更具体地,涉及一种发电机原动系统阻尼在线检测方法。
背景技术
超低频振荡是威胁电力系统安全运行的一个突出问题,当发生超低频振荡时,系统内所有机组共同振荡,频率和功率也发生持续振荡,严重威胁电网安全。超低频振荡是一次调频过程中出现的稳定问题,与发电机的原动机和调速系统密切相关,原动系统产生负阻尼是导致超低频振荡的重要原因之一。目前没有对原动系统阻尼在线评估方法,使得超低频振荡的快速有效控制成为难题。本发明通过在线评估计算发电机原动系统的阻尼转矩系数,能够在线得到各发电机原动系统的阻尼情况,当系统发生超低频振荡时可以快速确定提供负阻尼的原动系统并将其调速器退出,使振荡快速平息,以满足电网安全需求。
发明内容
本发明提供一种克服上述问题或者至少部分地解决上述问题的发电机原动系统阻尼在线测量方法,通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流,并在线计算发电机原动系统的阻尼转矩系数。
根据本发明的一个方面,提供一种发电机原动系统阻尼在线测量方法,包括以下步骤:
S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流;
S2、计算发电机流入励磁系统的励磁耗散能量流;
S3、根据电网流入发电机的耗散能量流和励磁耗散能量流之差计算发电机流入原动系统的原动耗散能量流;
S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合,并计算两者斜率比,得到原动系统的阻尼转矩系数。
作为优选的,所述步骤S1具体包括:
S11、监测电力系统振荡后发电机的有功功率振幅,确定有功功率振幅稳定的一段时间窗口;
S12、获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f,并进行滤波预处理;
S13、计算该时间窗口内每个采样点电网流入发电机的耗散能量流和每个采样点单位阻尼耗散能量流。
作为优选的,所述步骤S12具体包括:获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f;计算电压U的自然对数lnU,对P、Q、lnU、f进行高通滤波,去除直流分量和低频分量,得到相应偏差△P、△Q、△lnU、△f。
作为优选的,所述步骤S13具体包括:
对时间窗口内的每个采集样点t,计算电网流入发电机的发电机耗散能量流:
式中,为时刻t电网流入发电机的的耗散能量流,ts为计算的起始时刻,初始值
对时间窗口内每个采样点t,计算单位阻尼耗散能量流:
式中,f0为额定频率,为时刻t单位阻尼消耗能量流,ts为计算的起始时刻,初始值
作为优选的,所述步骤S13进一步包括,当采样数据为离散数据时,采样周期为h,则:
电网流入发电机的发电机耗散能量流为:
式中,下标i表示对应变量在t=ts+ih处的值,初始值
单位阻尼消耗能量流:
式中,下标i表示对应变量在t=ts+ih处的值,初始值
作为优选的,所述步骤S2包括:
S21、采集该段时间窗口内发电机的励磁电压Ufd、励磁电流Ifd并进行滤波预处理,得到励磁电压、励磁电流的偏差量△Ufd、△Ifd
S22、并计算计算机流入励磁系统的励磁耗散能量流。
作为优选的,所述步骤S22具体包括:
对于每个时间窗口内的每个采样点t,计算发电机流入励磁系统的励磁耗散能量流:
式中,为时刻t发电机流入励磁系统的励磁耗散能量流,ts为计算的起始时刻,初始值
当采集数据为离散数据时,采集周期为h,发电机流入励磁系统的励磁耗散能量流为:
式中,下标i表示对应变量在t=ts+ih处的值,初始值
作为优选的,所述步骤S3具体包括:
根据电网流入发电机的耗散能量流和励磁耗散能量流得到原动耗散能量流:
当采集数据为离散数据时,原动耗散能量流为:
本申请提出一种发电机原动系统阻尼在线检测方法,通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流,并在线计算发电机原动系统的阻尼转矩系数,获取各个发电机原动系统的阻尼情况,当系统发生超低频振荡时可以快速确定提供负阻尼的原动系统,并将其一次调频退出,使振荡可以快速平息,保证电网的安全性。
附图说明
图1为本发明实施例中的原动系统阻尼在线测量方法流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1示出了一种发电机原动系统阻尼在线检测方法,利用能量流方法计算发电机流入原动系统的耗散能量流,进而通过耗散能量流分析原动系统的阻尼作用。由于计算发电机流入原动系统的耗散能量流公式中很多变量在实际电网中无法测量或测量精度较低,因此本发明采取了间接计算的方法,通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的耗散能量流。将得到的能量流进行线性拟合得到斜率即为流入原动系统的耗散功率,该功率与单位阻尼耗散功率的比值即为原动系统阻尼转矩系数,具体包括以下步骤:
S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流;
S2、计算发电机流入励磁系统的励磁耗散能量流;
S3、根据发电机耗散能量流和励磁耗散能量流之差计算发电机流入原动系统的原动耗散能量流;
S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合,并计算两者斜率比,得到原动系统的阻尼转矩系数,该值为正代表原动系统提供正阻尼,否则为负阻尼。
作为优选的,所述步骤S1具体包括:
S11、监测电力系统振荡后发电机的有功功率振幅,确定有功功率振幅稳定的一段时间窗口;
S12、获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f,并进行滤波预处理;其中P、Q分别为发电机发出的有功功率和无功功率,f为发电机所连母线处的频率,U为发电机所连母线处电压。
S13、计算该时间窗口内每个采样点电网流入发电机的发电机耗散能量流和每个采样点单位阻尼耗散能量流。
作为优选的,所述步骤S12具体包括:获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f;计算电压U的自然对数lnU,对P、Q、lnU、f进行高通滤波,去除直流分量和低频分量,得到相应偏差△P、△Q、△lnU、△f。
作为优选的,所述步骤S13具体包括:
对时间窗口内的每个采集样点t,计算电网流入发电机的发电机耗散能量流:
式中,为时刻t电网流入发电机的的耗散能量流,ts为计算的起始时刻,初始值
当采样数据为离散数据时,采样周期为h,则:
电网流入发电机的发电机耗散能量流为:
式中,下标i表示对应变量在t=ts+ih处的值,初始值
对时间窗口内每个采样点t,计算单位阻尼耗散能量流:
式中,f0为额定频率,为时刻t单位阻尼消耗能量流,ts为计算的起始时刻,初始值
当数据为离散数据时,采样周期为h时,单位阻尼消耗能量流:
式中,下标i表示对应变量在t=ts+ih处的值,初始值
作为优选的,所述步骤S2包括:
S21、采集该段时间窗口内发电机的励磁电压Ufd、励磁电流Ifd并进行滤波预处理,得到励磁电压、励磁电流的偏差量△Ufd、△Ifd;Ufd、Ifd分别为励磁电压、流入发电机励磁绕组的励磁电流。
S22、并计算计算机流入励磁系统的励磁耗散能量流。
作为优选的,所述步骤S22具体包括:
对于每个时间窗口内的每个采样点t,计算发电机流入励磁系统的励磁耗散能量流:
式中,为时刻t发电机流入励磁系统的励磁耗散能量流,ts为计算的起始时刻,初始值
当采集数据为离散数据时,采集周期为h,发电机流入励磁系统的励磁耗散能量流为:
式中,下标i表示对应变量在t=ts+ih处的值,初始值
作为优选的,所述步骤S3具体包括:
根据电网流入发电机的耗散能量流和励磁耗散能量流的得到原动耗散能量流:
当采集数据为离散数据时,原动耗散能量流为:
最后,对发电机流入原动系统的耗散能量流和单位阻尼耗散能量流分别进行线性拟合,得到拟合结果Wpm D≈kpmt+cpm使Wpm D=kpmt+cpm计算两者斜率比值即为发电机原动系统的阻尼转矩系数。
本申请提出一种发电机原动系统阻尼在线检测方法,通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流,并在线计算发电机原动系统的阻尼转矩系数,获取各个发电机原动系统的阻尼情况,当系统发生超低频振荡时可以快速确定提供负阻尼的原动系统,并将其一次调频退出,使振荡可以快速平息,保证电网的安全性。
最后,本申请的方法仅为较佳的实施方案,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种发电机原动系统阻尼在线测量方法,其特征在于,包括以下步骤:
S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流;
S2、计算发电机流入励磁系统的励磁耗散能量流;
S3、根据电网流入发电机的耗散能量流和励磁耗散能量流之差计算出发电机流入原动系统的原动耗散能量流;
S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合,并计算两者斜率比,得到原动系统的阻尼转矩系数。
2.根据权利要求1所述的发电机原动系统阻尼在线测量方法,其特征在于,所述步骤S1具体包括:
S11、监测电力系统振荡后发电机的有功功率振幅,确定有功功率振幅稳定的一段时间窗口;
S12、获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f,并进行滤波预处理;
S13、计算该时间窗口内每个采样点电网流入发电机的发电机耗散能量流和每个采样点单位阻尼耗散能量流。
3.根据权利要求2所述的发电机原动系统阻尼在线测量方法,其特征在于,所述步骤S12具体包括:
获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f;
计算电压U的自然对数lnU,对P、Q、lnU、f进行高通滤波,去除直流分量和低频分量,得到相应偏差△P、△Q、△lnU、△f。
4.根据权利要求3所述的发电机原动系统阻尼在线测量方法,其特征在于,所述步骤S13具体包括:
对时间窗口内的每个采集样点t,计算电网流入发电机的发电机耗散能量流:
W I N D ( t ) = - ∫ t s t ( 2 π Δ P Δ f d t + Δ Q d ( Δ ln U ) )
式中,为时刻t电网流入发电机的耗散能量流,ts为计算的起始时刻,初始值
对时间窗口内每个采样点t,计算单位阻尼耗散能量流:
W u n i t D ( t ) = 2 π f 0 ∫ t s t Δf 2 d t
式中,f0为额定频率,为时刻t单位阻尼消耗能量流,ts为计算的起始时刻,初始值
5.根据权利要求4所述的发电机原动系统阻尼在线测量方法,其特征在于,所述步骤S13进一步包括,当采样数据为离散数据时,采样周期为h,则:
电网流入发电机的发电机耗散能量流为:
W I N k D = - Σ i = 1 k π ( ΔP i - 1 Δf i - 1 + ΔP i Δf i ) h + 1 2 ( ΔQ i - 1 + ΔQ i ) ( Δ ln U i - Δ ln U i - 1 )
式中,下标i表示对应变量在t=ts+ih处的值,初始值
单位阻尼消耗能量流:
W u n i t k D = π f 0 Σ i = 1 k ( Δf i - 1 2 + Δf i 2 ) h
式中,下标i表示对应变量在t=ts+ih处的值,初始值
6.根据权利要求5所述的发电机原动系统阻尼在线检测方法,其特征在于,所述步骤S2包括:
S21、采集该段时间窗口内发电机的励磁电压Ufd、励磁电流Ifd并进行滤波预处理,得到励磁电压、励磁电流的偏差量△Ufd、△Ifd
S22、并计算计算机流入励磁系统的励磁耗散能量流。
7.根据权利要求6所述的发电机原动系统阻尼在线检测方法,其特征在于,所述步骤S22具体包括:
对于每个时间窗口内的每个采样点t,计算发电机流入励磁系统的励磁耗散能量流:
W f d D ( t ) = - ∫ t s t ( ΔU f d ΔI f d - ΔI f d 2 R f d ) d t
式中,为时刻t发电机流入励磁系统的励磁耗散能量流,ts为计算的起始时刻,初始值
当采集数据为离散数据时,采集周期为h,发电机流入励磁系统的励磁耗散能量流为:
W f d k D ( t ) = - 1 2 Σ i = 1 k ( ΔU f d i - 1 ΔI f d i - 1 + ΔU f d i ΔI f d i - R f d ( ΔI f d i - 1 2 + ΔI f d i 2 ) ) h
式中,下标i表示对应变量在t=ts+ih处的值,初始值
8.根据权利要求7所述的发电机原动系统阻尼在线检测方法,其特征在于,所述步骤S3具体包括:
根据电网流入发电机的耗散能量流和励磁耗散能量流的得到原动耗散能量流:
W p m D ( t ) = W I N D ( t ) - W f d D ( t )
当采集数据为离散数据时,原动耗散能量流为:
W p m k D = W I N k D - W f d k D .
CN201611209173.XA 2016-12-23 2016-12-23 一种发电机原动系统阻尼在线检测方法 Active CN106772042B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611209173.XA CN106772042B (zh) 2016-12-23 2016-12-23 一种发电机原动系统阻尼在线检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611209173.XA CN106772042B (zh) 2016-12-23 2016-12-23 一种发电机原动系统阻尼在线检测方法

Publications (2)

Publication Number Publication Date
CN106772042A true CN106772042A (zh) 2017-05-31
CN106772042B CN106772042B (zh) 2019-06-18

Family

ID=58920232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611209173.XA Active CN106772042B (zh) 2016-12-23 2016-12-23 一种发电机原动系统阻尼在线检测方法

Country Status (1)

Country Link
CN (1) CN106772042B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108896908A (zh) * 2018-04-03 2018-11-27 上海电器科学研究所(集团)有限公司 电动机有限元振动计算模型中阻尼参数的精确识别方法
CN110417054A (zh) * 2019-07-30 2019-11-05 华北电力大学 一种双馈风电并网系统稳定性调整系统及方法
CN110875599A (zh) * 2019-11-04 2020-03-10 云南电网有限责任公司 一种电网频率振荡的控制方法及系统
CN110879324A (zh) * 2019-12-02 2020-03-13 国网湖南省电力有限公司 一种励磁系统在功率振荡中阻尼特性的辨别方法
CN111555312A (zh) * 2020-05-27 2020-08-18 四川大学 一种适用于电力系统超低频振荡稳定性评估的方法
CN112861328A (zh) * 2021-01-22 2021-05-28 东北电力大学 一种基于随机响应信号的发电机阻尼评估装置及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201303308Y (zh) * 2008-11-19 2009-09-02 湖南大学 一种发电机原动系统动态仿真的装置
CN102944798A (zh) * 2012-11-29 2013-02-27 武汉华中电力电网技术有限公司 一种负阻尼低频振荡与强迫功率振荡判别方法
CN103117547A (zh) * 2013-03-04 2013-05-22 中国南方电网有限责任公司 一种电力系统低频振荡的振荡源定位方法
CN104111405A (zh) * 2014-06-23 2014-10-22 华中科技大学 一种基于阻尼转矩分析法的电力系统低频振荡源定位方法
CN106058897A (zh) * 2016-07-28 2016-10-26 东南大学 一种基于相量的发电机强迫扰动源定位方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201303308Y (zh) * 2008-11-19 2009-09-02 湖南大学 一种发电机原动系统动态仿真的装置
CN102944798A (zh) * 2012-11-29 2013-02-27 武汉华中电力电网技术有限公司 一种负阻尼低频振荡与强迫功率振荡判别方法
CN103117547A (zh) * 2013-03-04 2013-05-22 中国南方电网有限责任公司 一种电力系统低频振荡的振荡源定位方法
CN104111405A (zh) * 2014-06-23 2014-10-22 华中科技大学 一种基于阻尼转矩分析法的电力系统低频振荡源定位方法
CN106058897A (zh) * 2016-07-28 2016-10-26 东南大学 一种基于相量的发电机强迫扰动源定位方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LEI CHEN ET AL.: "An Energy-Based Method for Location of Power System Oscillation Source", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
陈磊 等: "基于振荡能量的低频振荡分析与振荡源定位(一)理论基础与能量流计算", 《电力系统自动化》 *
陈磊 等: "基于振荡能量的低频振荡分析与振荡源定位(二)振荡源定位方法与算例", 《电力系统自动化》 *
陈磊 等: "振荡能量分析和特征值分析的联系及发电机阻尼评估", 《电力系统自动化》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108896908A (zh) * 2018-04-03 2018-11-27 上海电器科学研究所(集团)有限公司 电动机有限元振动计算模型中阻尼参数的精确识别方法
CN108896908B (zh) * 2018-04-03 2021-07-09 上海电器科学研究所(集团)有限公司 电动机有限元振动计算模型中阻尼参数的精确识别方法
CN110417054A (zh) * 2019-07-30 2019-11-05 华北电力大学 一种双馈风电并网系统稳定性调整系统及方法
CN110417054B (zh) * 2019-07-30 2020-10-16 华北电力大学 一种双馈风电并网系统稳定性调整系统及方法
CN110875599A (zh) * 2019-11-04 2020-03-10 云南电网有限责任公司 一种电网频率振荡的控制方法及系统
CN110879324A (zh) * 2019-12-02 2020-03-13 国网湖南省电力有限公司 一种励磁系统在功率振荡中阻尼特性的辨别方法
CN110879324B (zh) * 2019-12-02 2022-05-17 国网湖南省电力有限公司 一种励磁系统在功率振荡中阻尼特性的辨别方法
CN111555312A (zh) * 2020-05-27 2020-08-18 四川大学 一种适用于电力系统超低频振荡稳定性评估的方法
CN112861328A (zh) * 2021-01-22 2021-05-28 东北电力大学 一种基于随机响应信号的发电机阻尼评估装置及方法
CN112861328B (zh) * 2021-01-22 2022-08-30 东北电力大学 一种基于随机响应信号的发电机阻尼评估装置及方法

Also Published As

Publication number Publication date
CN106772042B (zh) 2019-06-18

Similar Documents

Publication Publication Date Title
CN106772042A (zh) 一种发电机原动系统阻尼在线检测方法
DE112012007019B4 (de) Motorsteuergerät für ein elektrisches Fahrzeug
CN102288283B (zh) 汽轮发电机组振动故障的监测方法及装置
CN102545245B (zh) 基于端口供给能量的电力系统振荡源定位方法
CN109312716A (zh) 用于控制风力涡轮机以管理缘向叶片振动的方法和系统
CN103344914B (zh) 基于归一化的风电机组故障预警方法
CN104865474B (zh) 一种基于pmu数据实时监测低频振荡源的方法
CN102928697B (zh) 励磁调节器pss2a模型低频段阻尼检测方法及系统
CN106948796A (zh) 一种抽油机不停机间抽控制装置及方法
CN103117547A (zh) 一种电力系统低频振荡的振荡源定位方法
CN103618492A (zh) 一种基于时频变换的同步发电机参数辨识方法
CN104821579A (zh) 一种基于换流站电气信号的次同步振荡监测分析方法
DE102012201200A1 (de) Fahrzeug-Motorsteuervorrichtung
CN109283427B (zh) 含双馈风电机组的互联电网低频扰动源定位系统及方法
CN104280682B (zh) 一种基于磁场定向控制的电机转子故障诊断方法
CN107104448A (zh) 电力系统稳定器及方法
CN104777396A (zh) 一种综合考虑电机电流和电压的转子故障诊断方法
CN105891625A (zh) 基于能量流的电力系统次同步振荡扰动源辨别方法
CN106058897B (zh) 一种基于相量的发电机强迫扰动源定位方法
CN104101805B (zh) 一种基于振荡能量注入的励磁系统负阻尼检测方法
CN104638638A (zh) 一种用于大电网的在线安全稳定趋势分析方法
CN108988387A (zh) 基于转子电流移相平均的双馈风电机组次同步谐振抑制方法
CN103511553A (zh) 一种拍振噪声的控制设备、方法、系统及工程机械
CN101026302A (zh) 混合型有源滤波器的分频控制方法
CN109921439B (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