CN106771594B - 一种电力系统的次/超同步谐波检测方法 - Google Patents

一种电力系统的次/超同步谐波检测方法 Download PDF

Info

Publication number
CN106771594B
CN106771594B CN201611122173.6A CN201611122173A CN106771594B CN 106771594 B CN106771594 B CN 106771594B CN 201611122173 A CN201611122173 A CN 201611122173A CN 106771594 B CN106771594 B CN 106771594B
Authority
CN
China
Prior art keywords
frequency
supersynchronous
group
harmonic
signal
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
CN201611122173.6A
Other languages
English (en)
Other versions
CN106771594A (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.)
Guo Wang Xinjiang Power Co
Tsinghua University
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
Guo Wang Xinjiang Power Co
Tsinghua University
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
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 Guo Wang Xinjiang Power Co, Tsinghua University, State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI filed Critical Guo Wang Xinjiang Power Co
Priority to CN201611122173.6A priority Critical patent/CN106771594B/zh
Publication of CN106771594A publication Critical patent/CN106771594A/zh
Application granted granted Critical
Publication of CN106771594B publication Critical patent/CN106771594B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16—Spectrum analysis; Fourier analysis
    • G01R23/165—Spectrum analysis; Fourier analysis using filters

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种电力系统的次/超同步谐波检测方法，包括：提取并分离工频分量；根据模态混叠条件进行分组；提取含有最高频率谐波分量的分组信号并判断是否发生模态混叠；发生模态混叠时进行频移调制，经验模态分解得到该组所含的各次/超同步谐波对应的固有模态函数；进行幅值相位和频率补偿；希尔伯特变换得到该组次/超同步谐波的瞬时频率和幅值；分离该组次/超同步谐波；重复上述步骤直至分离所有分组。本发明具有如下优点：基于信号调制的频移方法，通过增大两个次/超同步谐波分量的频率比，扩展经验模态分解的信号分解范围，适用于电力系统次/超同步谐波的实时精确检测。

Description

一种电力系统的次/超同步谐波检测方法
技术领域
本发明涉及电力系统动态监测与电能质量技术领域，具体涉及一种电力系统的次/超同步谐波检测方法。
背景技术
近年来，风电、光伏等新型可再生能源发电迅速扩展，其渗透率不断提高，由风电机组控制器与交流电网之间的相互作用引起的次同步振荡问题日益凸显，影响电网和设备的安全运行，甚至造成严重的稳定性事故或设备损伤，危害极大。
次同步振荡在电网中主要表现为幅值较大的持续增长(发散)或恒幅的次同步和/或超同步电流、电压和功率谐波，准确检测次/超同步谐波的幅值、频率等特征参数，有助于深入分析电力系统次同步振荡特性，从而有效抑制振荡，保障机组及电力系统稳定运行。
电力系统次/超同步谐波具有非线性、非平稳性特性，而传统的傅里叶方法不能准确得到振荡模式的时变特性，因此难以有效检测次/超同步谐波。1998年，美国宇航局Norden E.Huang等人提出了针对非线性、非平稳信号的自适应分析方法——希尔伯特·黄变换(Hilbert-Huang Transform，HHT)，HHT分为两步，第一步为经验模态分解(EmpiricalMode Decomposition,EMD)，第二步为希尔伯特变换。其中EMD是HHT的核心部分，将信号中多个组成分量分解成一系列固有模态函数(Intrinsic Mode Function,IMF)，使得每一瞬间每一个IMF只有一个频率，然后经过希尔伯特变换得到有实际物理意义的每一个IMF的瞬时幅值和瞬时频率。但EMD方法应用中存在端点效应、模态混叠等问题，尤其是当信号中两个组成分量的频率在2倍频内时，EMD无法将两者分解开，产生模态混叠现象，降低了频率分辨率。由于次/超同步谐波与基波或者次/超同步谐波之间的频率接近2倍频或在2倍频内，因此这个缺点制约了HHT在次/超同步谐波检测中的应用。
发明内容
本发明旨在至少解决上述技术问题之一。
为此，本发明的目的在于提出一种电力系统的次/超同步谐波检测方法，解决模态混叠问题。
为了实现上述目的，本发明的实施例公开了一种电力系统的次/超同步谐波检测方法，包括以下步骤：S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离；S2：通过对所述实测信号进行快速傅里叶变换(Fast Fourier Transform,FFT)，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，直接进行EMD会发生混叠，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组；S3：提取含有所述实测信号未分离部分最高频率谐波分量的分组信号，根据预设混叠设定条件判断所述分组信号是否使用信号频移方法进行调制；S4：如果所述分组信号不需要调制，EMD得到所述分组所含的各次/超同步谐波对应的IMF，如果所述分组信号需要调制，则对经过调制的信号进行EMD，得到所述分组所含的各次/超同步谐波经过调制后对应的IMF；S5：对所有IMF的幅值、相位和频率根据预设补偿条件进行补偿；S6：对补偿后的函数进行希尔伯特变换得到所述分组所含的各次/超同步谐波的瞬时频率和瞬时幅值；S7：将所述分组信号从所述实测信号中分离；S8：重复S3-S7，直至将所述实测信号中全部分组进行分离。
根据本发明实施例的电力系统的次/超同步谐波检测方法，通过增大两个次/超同步谐波分量的频率比，扩展EMD的信号分解范围，适用于电力系统次/超同步谐波的实时精确检测
另外，根据本发明上述实施例的电力系统的次/超同步谐波检测方法，还可以具有如下附加的技术特征：
进一步地，在步骤S2中，通过以下步骤对所述实测信号所含的次/超同步谐波分量进行分组：将所述实测信号所含的次/超同步谐波分量按照频率的高低进行降序排列，计算两个相邻谐波分量的频率比；将频率比小于预设频率比阈值的两个谐波分量分为一组，对得到的多个分组按照组成分量频率的高低进行降序排列；将频率位于两个相邻分组所含谐波频率之间的剩余谐波分量分为一组。
进一步地，在步骤S3中，使用高通滤波器得到含有最高频率谐波分量的分组信号。
进一步地，在步骤S4中，在得到所述多个IMF之后还包括：筛除包含混叠频率的IMF。
进一步地，所述预设倍数为5倍。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是本发明一个实施例的电力系统的次/超同步谐波检测方法的流程图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
参照下面的描述和附图，将清楚本发明的实施例的这些和其他方面。在这些描述和附图中，具体公开了本发明的实施例中的一些特定实施方式，来表示实施本发明的实施例的原理的一些方式，但是应当理解，本发明的实施例的范围不受此限制。相反，本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
以下结合附图描述本发明。
图1是本发明一个实施例的电力系统的次/超同步谐波检测方法的流程图。如图1所示，一种电力系统的次/超同步谐波检测方法，包括以下步骤：
S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离。
具体地，实测信号s(t)中，如果工频分量幅值是频率与其相近的次/超同步谐波幅值的5倍以上，为了避免工频分量对该次/超同步谐波的特征参数提取产生干扰，将工频分量从原信号中提取并分离，具体过程如下：
首先使用带通滤波器获得所述工频分量，经过EMD得到相应的IMF，记为g(t)，进行幅值相位补偿后，通过希尔伯特变换得到所述工频分量的瞬时幅值和瞬时频率，再将所述工频分量从实测信号中分离。
S2：通过对所述分离出工频分量的实测信号进行FFT，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，直接进行EMD会发生混叠，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组。
具体地，假设s(t)中含有N个次/超同步谐波分量，根据预设的频率分辨率，选择相应的数据窗长对s(t)进行FFT计算，得到N个次/超同步谐波分量的大致频率值fn(n＝1,2,…,N),其中f1>f2>f3>…>fN，将N个次/超同步谐波分量按照频率的高低进行降序排列，计算两个相邻谐波分量的频率比mn,mn＝fn/fn+1(n＝1,2,…,N-1)；根据混叠现象的产生条件，将所述实测信号所含的次/超同步谐波分量进行分组，具体过程如下：将mn(1≤n≤N-1)依次与2进行比较，若mn≤2，则将频率为fn和fn+1的组成分量分为一组，若此时mn+1≤2同时成立，由于每个次/超同步谐波分量都唯一存在于某一分组中，不再将频率为fn+1和fn+2的组成分量分为一组，而继续根据mn+2的大小考虑是否将频率为fn+2和fn+3的组成分量分为一组；假设共得到K个组，按照每组所含组成分量的频率高低，依次记为Hi(i＝1,2,…,K)；
将剩余谐波分量根据各自频率所在区间进行分组，将频率高于H1中频率最高谐波分量的剩余谐波分量分为一组；将频率低于Hi(1≤i≤K-1)中频率最低谐波分量且高于Hi+1(1≤i≤K-1)中频率最高谐波分量的剩余谐波分量分为一组；将频率低于HK中频率最低谐波分量的剩余谐波分量分为一组；假设共得到J个组，按照每组所含组成分量的频率高低，依次记为Wi(i＝1,2,…,J)。
至此，将实测信号所含的次/超同步谐波分量全部进行分组，其中，Hi(i＝1,2,…,K)分别由直接进行EMD会发生混叠的两个谐波分量组成，Wi(i＝1,2,…,J)分别由直接进行EMD互不发生混叠的若干个谐波分量组成。
S3：分离出含有实测信号未分离部分最高频率谐波分量的一组信号，根据预设混叠设定条件判断所述分组信号是否使用信号频移方法进行调制。
具体地，使用高通滤波器，得到含有最高频率谐波分量的分组信号，若该组为发生混叠的谐波分量组成的分组Hi，首先使用信号频移方法，将该组信号进行调制，具体过程如下：
将Hi中包含的两个次/超同步谐波频率分别记为fa和fb，其中fa>fb，设置调制信号为mi(t)＝Asin(2πfmit)，其中频率fmi需要满足2fb-fa<fmi<fb，幅值A为任意大于零的实数；
设调制后信号为li(t),li(t)＝Hi(t)mi(t),li(t)包含四个频率fa-fmi,fb-fmi,fa+fmi,fb+fmi，其中fa-fmi,fb-fmi为Hi的两个组成分量经过频移过程对应的新的频率，满足两个组成分量直接进行EMD不发生混叠的条件。
S4：如果所述分组信号不需要调制，EMD得到所述分组所含的各次/超同步谐波的IMF，如果所述分组需要调制，则对经过调制的信号进行EMD，得到所述分组所含的各次/超同步谐波经过调制后对应的IMF。
具体地，若该组次/超同步谐波不需要调制，直接进行EMD，若该组次/超同步谐波需要调制，则将调制后的信号进行EMD。在得到的若干个IMF中，经过筛选去除无效的包含混叠频率的IMF，得到该组各次/超同步谐波分别对应的仅含单一频率的IMF，记为fm(t)(m＝1,2,…)。
S5：对所有IMF的幅值、相位和频率根据预设补偿条件进行补偿。
具体地，将fm(t)(m＝1,2,…)进行幅值、相位和频率补偿，具体步骤如下：
若该组次/超同步谐波经过调制过程，则得到的fm(t)(m＝1,2)的幅值为对应的次/超同步谐波幅值与调制信号幅值A的乘积的一半，fm(t)的频率为对应的次/超同步谐波频率fm与调制信号频率fmi之差，因此根据设置的调制信号依次进行幅值补偿和频率补偿；
计算该组次/超同步谐波在步骤(3)的滤波环节中在各自频率fm(m＝1,2,…)上的幅频响应αm和相频响应将fm(t)的幅值除以αm，相位减去将补偿后的函数记为hm(t)(m＝1,2,…)。
S6：对补偿后的函数进行希尔伯特变换得到所述分组所含的各次/超同步谐波的瞬时频率和瞬时幅值。
具体地，将hm(t)(m＝1,2,…)进行希尔伯特变换，得到该组所含的次/超同步谐波的瞬时频率和瞬时幅值
S7：将所述分组信号从所述实测信号中分离。
具体地，从原信号中分离出该组次/超同步谐波，s(t)＝s(t)-∑hm(t)。
S8：重复S3-S7，直至将所述实测信号中全部分组进行分离。
另外，本发明实施例的电力系统的次/超同步谐波检测方法的其它构成以及作用对于本领域的技术人员而言都是已知的，为了减少冗余，不做赘述。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同限定。

Claims (4)

1.一种电力系统的次/超同步谐波检测方法，其特征在于，包括以下步骤：
S1：假设实测信号包含工频分量和多个次/超同步谐波分量，如果工频分量幅值超过频率与其相近的次/超同步谐波幅值预设倍数时，将所述工频分量从原信号中提取并分离；
S2：通过对所述实测信号进行快速傅里叶变换，得到所述实测信号所含的各次/超同步谐波的大致频率，将频率比小于频率比阈值的两个谐波分量分为一组，将剩余谐波分量根据各自频率所在区间进行分组；从而将所述实测信号所含的次/超同步谐波分量全部进行分组；其中，通过以下步骤对所述实测信号所含的次/超同步谐波分量进行分组：将所述实测信号所含的次/超同步谐波分量按照频率的高低进行降序排列，计算两个相邻谐波分量的频率比；将频率比小于等于预设频率比阈值的两个谐波分量分为一组，从而得到按照每组所含组成分量的频率自高而低降序排列的K个组，依次记为Hi，i＝1,2,…,K；将频率高于所述K个组的第一个组H1中频率最高谐波分量的剩余谐波分量分为一组，将频率位于所述K个组的两个相邻分组所含谐波频率之间的剩余谐波分量分为一组，将频率低于所述K个组的最后一个组HK中频率最低谐波分量的剩余谐波分量分为一组；
S3：提取含有所述实测信号未分离部分最高频率谐波分量的分组信号，根据预设混叠设定条件判断所述分组信号是否使用信号频移方法进行调制；
S4：如果所述分组信号不需要调制，经验模态分解得到所述分组所含的各次/超同步谐波对应的IMF，如果所述分组信号需要调制，则对经过调制的信号进行经验模态分解，得到所述分组所含的各次/超同步谐波经过调制后对应的固有模态函数；
S5：对所有固有模态函数的幅值、相位和频率根据预设补偿条件进行补偿；
S6：对补偿后的函数进行希尔伯特变换得到所述分组所含的各次/超同步谐波的瞬时频率和瞬时幅值；
S7：将所述分组信号从所述实测信号中分离；
S8：重复S3-S7，直至将所述实测信号中全部分组进行分离。
2.根据权利要求1所述的电力系统的次/超同步谐波检测方法，其特征在于，在步骤S3中，使用高通滤波器得到含有最高频率谐波分量的分组信号。
3.根据权利要求1所述的电力系统的次/超同步谐波检测方法，其特征在于，在步骤S4中，在得到所述多个固有模态函数之后还包括：
筛除包含混叠频率的固有模态函数。
4.根据权利要求1所述的电力系统的次/超同步谐波检测方法，其特征在于，所述预设倍数为5倍。
CN201611122173.6A 2016-12-08 2016-12-08 一种电力系统的次/超同步谐波检测方法 Active CN106771594B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611122173.6A CN106771594B (zh) 2016-12-08 2016-12-08 一种电力系统的次/超同步谐波检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611122173.6A CN106771594B (zh) 2016-12-08 2016-12-08 一种电力系统的次/超同步谐波检测方法

Publications (2)

Publication Number Publication Date
CN106771594A CN106771594A (zh) 2017-05-31
CN106771594B true CN106771594B (zh) 2019-08-09

Family

ID=58881488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611122173.6A Active CN106771594B (zh) 2016-12-08 2016-12-08 一种电力系统的次/超同步谐波检测方法

Country Status (1)

Country Link
CN (1) CN106771594B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966611B (zh) * 2017-11-24 2019-10-22 广东电网有限责任公司电力调度控制中心 一种基于矢量匹配的电力系统次超同步谐波检测方法
CN108957129A (zh) * 2018-07-27 2018-12-07 中国南方电网有限责任公司 一种新能源并网系统次/超同步振荡自适应监测方法
CN111413578B (zh) * 2019-05-29 2022-07-05 中国电力工程顾问集团华北电力设计院有限公司 一种次同步振荡的实时监测预警方法
CN110837001B (zh) * 2019-11-13 2021-10-01 河南工程学院 一种电力系统中谐波和间谐波的分析方法与装置
CN111257636B (zh) * 2020-02-25 2021-05-11 北京交通大学 一种电力系统次同步和超同步振荡频率辨识方法
CN111984920B (zh) * 2020-08-31 2022-03-18 广东电网有限责任公司广州供电局 次/超同步谐波参数识别方法、装置、设备和介质
CN113608018A (zh) * 2021-06-30 2021-11-05 中冶南方都市环保工程技术股份有限公司 提高谐波检测精度的自适应vmd检测方法、装置及存储介质
CN117076933B (zh) * 2023-10-16 2024-02-06 锦浪科技股份有限公司 电弧判断模型的训练、光伏直流电弧检测方法及计算设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323480A (zh) * 2011-05-19 2012-01-18 西南交通大学 一种基于Hilbert-Huang变换的电能质量分析方法
CN103163372A (zh) * 2013-03-26 2013-06-19 山西省电力公司长治供电分公司 采用hht的电力系统谐波分析方法
CN105510711A (zh) * 2015-12-24 2016-04-20 合肥工业大学 一种改进的经验模态分解的谐波分析法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323480A (zh) * 2011-05-19 2012-01-18 西南交通大学 一种基于Hilbert-Huang变换的电能质量分析方法
CN103163372A (zh) * 2013-03-26 2013-06-19 山西省电力公司长治供电分公司 采用hht的电力系统谐波分析方法
CN105510711A (zh) * 2015-12-24 2016-04-20 合肥工业大学 一种改进的经验模态分解的谐波分析法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
希尔伯特—黄变换（HHT）及其在结构动态性能分析中的应用研究;汪家慰;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20110415(第04期);正文第5页第2-3段、第18页第1段-第24页最后1段 *
改进的 Hilbert- Huang 变换在船舶电网谐波检测中的应用;苏玉香等;《电工电能新技术》;20110731;第30卷(第3期);摘要、第10页第2栏最后1段-第12页第2栏第2段 *
电力系统次同步和超同步谐波相量的检测方法;谢小荣等;《电力系统自动化》;20161110;第40卷(第21期);摘要、第189页第2栏最后1段-191页第2栏第4段 *

Also Published As

Publication number Publication date
CN106771594A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106771594B (zh) 一种电力系统的次/超同步谐波检测方法
Achlerkar et al. Variational mode decomposition and decision tree based detection and classification of power quality disturbances in grid-connected distributed generation system
Gao et al. Waveform-based method for fast and accurate identification of subsynchronous resonance events
CN104638664B (zh) 风电厂次同步振荡保护方法及保护装置
US20120303306A1 (en) Method of subsynchronous resonance detection
CN108957129A (zh) 一种新能源并网系统次/超同步振荡自适应监测方法
Langella et al. On the interharmonic emission of PV inverters under different operating conditions
CN107345984A (zh) 一种基于信号识别的自适应同步相量测量方法
CN108535613B (zh) 一种基于组合窗函数的电压闪变参数检测方法
WO2016065959A1 (zh) 中性点不接地的10kV系统中铁磁谐振的诊断方法
CN108169540A (zh) 一种风力发电机组电压闪变的测量方法
Prince et al. Targeted approach to apply masking signal-based empirical mode decomposition for mode identification from dynamic power system wide area measurement signal data
CN105891625B (zh) 基于能量流的电力系统次同步振荡扰动源辨别方法
CN107179440A (zh) 一种在线自适应频率变化的次同步振荡间谐波提取方法
Abdeen et al. Improved subsynchronous oscillation detection method in a DFIG-based wind farm interfaced with a series-compensated network
CN105223483A (zh) 一种直流输电工程绝缘子串缺陷在线检测方法及装置
Yang et al. Research on a power quality monitoring technique for individual wind turbines
CN104466958B (zh) 电网电压故障下自适应分数阶锁相环的获取方法
Vanfretti et al. A PMU-based fast real-time sub-synchronous oscillation detection application
CN105301391B (zh) 基于有源滤波及谐波补偿的孤岛检测装置和方法
Shuai et al. Application of type 4 wind turbine harmonic model for wind power plant harmonic study
AlAli et al. Impact of increasing wind power generation on the north-south inter-area oscillation mode in the European ENTSO-E system
CN108490254A (zh) 一种基于多算法融合的宽带多频信号测量方法
Amirat et al. Phasor estimation for power quality monitoring: Least square versus Kalman filter
CN106610446A (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