CN104008279A - 一种求解电力网络静态安全域的方法 - Google Patents
一种求解电力网络静态安全域的方法 Download PDFInfo
- Publication number
- CN104008279A CN104008279A CN201410201937.5A CN201410201937A CN104008279A CN 104008279 A CN104008279 A CN 104008279A CN 201410201937 A CN201410201937 A CN 201410201937A CN 104008279 A CN104008279 A CN 104008279A
- Authority
- CN
- China
- Prior art keywords
- node
- electric power
- security
- power networks
- security domain
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000003068 static effect Effects 0.000 title claims abstract description 16
- 238000002347 injection Methods 0.000 claims abstract description 22
- 239000007924 injection Substances 0.000 claims abstract description 22
- 239000000243 solution Substances 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 9
- 235000013350 formula milk Nutrition 0.000 description 13
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种求解电力网络静态安全域的方法。该方法线性化电力网络节点功率注入方程,建立节点电压幅值与相角和节点功率注入之间的线性关系,进而得到代表所有对应于静态安全约束的子边界的超平面组。通过***电力网络,把断面潮流安全域等价为***后两个子电力网络***节点安全域的交集,间接地获得了断面潮流安全域边界超平面组。所提方法采用了线性化方法,解决了电力网络节点功率注入方程超高维数和强非线性问题,使得求解安全域边界的计算量小,所需计算时间少,适合于实际在线应用。本发明可广泛应用于电力网络的静态安全性分析。
Description
技术领域
本发明涉及一种求解电力网络静态安全域的方法,属于电力***静态安全分析领域。
背景技术
电力***安全直接影响到人们的日常生活和社会稳定。近年来世界范围内大停电事故频发,导致了巨大的经济损失和社会负面影响,因此研究能够实时在线应用的电力***静态安全分析方法十分必要。
电力网络静态安全分析是保障电力***安全十分有效的一种工具。传统静态安全分析方法“逐点法”针对给定功率注入和故障条件,采用数值仿真、特征值计算或能量函数法等方法来分析电力***的安全性。这类方法是当前工业界进行静态安全分析最主要和最实用的方法,但其缺点亦很多,比如功率注入给定的人为主观性、选择故障集的历史经验性和计算过程的耗时低效。另外这类方法只针对一个或多个事件,无法对***的安全性进行全局性分析。
电力网络静态安全域能提供丰富的安全信息。目前拟合法是常有的计算静态安全域的方法。该方法首先采用仿真方法获得大量安全域边界上的运行点,基于获得的临界运行点采用最小二乘法拟合来安全域边界。然而对于大规模的电力***,这种方法由于需要大量数值仿真造成计算负担繁重而难以实际在线应用。而且采用仿真方法来获得安全域边界上的运行点,在操作上依然带有人为主观性的缺点。
另外,当前电力网络静态安全域方法都是针对电力网络节点功率注入的,不获得电力网络中由几条线路组成的断面的安全域。而电力***实际运行中,由几条线路构成的某些关键断面的安全性可能是调度员非常关心的。因此,研究断面线路潮流的安全域具有十分重要的实用价值。
发明内容
本发明所要解决的技术问题是提供一种求解电力网络静态安全域的方法,解决现有电力网络静态安全域分析方法的不足和不能求取电力网络断面线路潮流安全域的缺点。该方法解决了电力网络节点功率注入方程超高维数和强非线性问题,使得求解安全域边界的计算量小,所需计算时间少,适合于实际在线应用。
本发明为解决上述技术问题采用以下技术方案:
本发明提供一种求解电力网络静态安全域的方法,首先通过线性化电力网络节点功率注入方程,建立节点电压幅值与相角和节点功率注入之间的线性关系,得到代表静态安全域边界的超平面组;然后通过***电力网络,将求解断面潮流安全域转化为求解***节点安全域,间接地获得了断面潮流安全域边界超平面组。
本发明的具体实施步骤如下:
步骤1,输入电力网络数据;所述电力网络数据包括电力网络的结构数据、电气数据、静态安全约束数据;所述结构数据包括电力网络中线路的连接关系、线路与变压器的连接关系、发电机和负荷的位置;所述电气数据包括线路及变压器的电阻和电抗、发电机节点和负荷节点的功率注入;所述静态安全约束数据包括节点电压幅值安全约束、线路潮流安全约束;
步骤2,根据电力网络的拓扑结构,建立节点功率注入方程,即节点有功、节点无功注入方程;
步骤3,在理想运行点(V=1,θ=0)的小邻域内,线性化步骤2中建立的节点功率注入方程,得到节点电压幅值与相角和节点功率注入之间的线性关系;其中,V是节点电压幅值向量;θ是节点电压相角向量;
步骤4,将静态安全约束数据代入步骤3中的节点电压幅值与相角和节点功率注入之间的线性关系方程,得到由线性方程表达的电力网络节点功率注入静态安全域边界;
步骤5,将电力网络从断面***开,***后所得的两个子电力网络都保留联络线和***节点;
步骤6,求解步骤5中每个子电力网络的***节点的功率注入安全域,则原电力网络断面潮流安全域等价于两个子网络的***节点的功率注入安全域的交集。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)本发明通过线性化电力网络节点功率注入方程,得到节点电压幅值与相角和节点功率注入之间的线性关系,将电力网络节点功率注入静态安全域边界分解为节点电压幅值约束子边界和线路潮流约束子边界,解决了求解电力网络节点功率注入静态安全域边界问题的超高维数和强非线性问题;
(2)本发明把电力网络从断面裂开,把求解电力网络断面潮流静态安全域边界问题转化为求解电力网络节点功率注入静态安全域边界问题,解决了求解电力网络断面潮流静态安全域边界问题的超高维数和强非线性问题;
(3)本发明通过线性化节点功率注入方程解析地求取了安全域,方法通用性强,实现方法简单,便于推广应用;
(4)本发明中安全域边界计算量小速度快,适合大电力***的实时在线静态安全分析。
附图说明
图1是本发明的方法流程图。
图2是IEEE14节点标准电力网络。
图3是IEEE14节点标准电力网络中节点3和节点8的电压幅值约束子安全域和线路潮流约束子安全域,其中,(a)为考虑节点3和节点8的电压幅值约束子安全域超平面组;(b)为考虑节点3和节点8的线路潮流约束子安全域超平面组。
图4是IEEE14节点标准电力网络中节点3和节点8的有功功率注入空间安全域。
图5是电力网络***方法示意图,其中,(a)为两区域互联电力网络;(b)为***后的子网络示意图。
图6是IEEE18节点标准电力网络。
图7是IEEE18节点标准电力网络断面中线路23至24和线路30至38潮流安全域。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解的是,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明设计一种求解电力网络静态安全域的方法,线性化电力网络节点功率注入方程,建立节点电压幅值与相角和节点功率注入之间的线性关系,得到代表静态安全域边界的超平面组;通过***电力网络,将求解断面潮流安全域转化为求解***节点安全域,间接地获得了断面潮流安全域边界超平面组。
如图1所示,本发明一种求解电力网络静态安全域的方法的具体实施步骤如下:
步骤1,输入电力网络数据;所述电力网络数据包括电力网络的结构数据、电气数据、静态安全约束数据;所述结构数据包括电力网络中线路的连接关系、线路与变压器的连接关系、发电机和负荷的位置;所述电气数据包括线路及变压器的电阻和电抗、发电机节点和负荷节点的功率注入;所述静态安全约束数据包括节点电压幅值安全约束、线路潮流安全约束。
其中,电压幅值安全约束和线路潮流安全约束为:
Vmin≤V≤Vmax (1)
-δ≤ATθ≤δ (2)
式中,V是节点电压幅值向量;θ是节点电压相角向量;Vmin为节点电压幅值下限约束向量;Vmax为节点电压幅值上限约束向量;δ为线路最大允许相角差;A为电力网络节点支路关联矩阵。
步骤2,根据电力网络的拓扑结构,建立节点功率注入方程,即节点有功、节点无功注入方程,具体如下:
式中,Pi为节点i的有功注入功率;Qi为节点i的无功注入功率;Vi为节点i的电压幅值;Vj为节点j的电压幅值;θij为节点i和节点j之间的电压相角差;G和B为n×n阶矩阵,其中G和B的对角元素和非对角元素分别为 和 rij为节点i和节点j之间支路ij的电阻;xij为节点i和节点j之间支路ij的电抗;i=1,2...,n;j=1,2...,n;n为电力网络节点数。
步骤3,线性化节点功率注入方程,具体如下:
将公式(3)改写成如下的紧凑形式:
线路潮流的计算式为
式中,Pij和Qij分别为支路ij的有功潮流和无功潮流;
重写公式(5)为
其中,Γ1、Γ2为中间变量,Γ1=Vjcosθij-Vi,Γ2=Vjsinθij。
在理想运行点(V=1,θ=0)的小邻域内,Vi=1,则得到
进一步,利用cosθij=1,得到Γ1=Vj-Vi;利用Vj=1和sinθij=θi-θj,得到Γ2=Vjsinθij;
最终可得
对节点i应用基尔霍夫定律可以得到
最终的得到
公式(10)写成矩阵形式为
公式(11)是一组线性方程,则容易求得以下表达式
公式(12)给出了以节点功率注入表达的电压幅值和相角方程,是经线性化后得到的状态量空间到节点功率注入空间的线性映射。
步骤4,求解电力网络节点功率注入静态安全域边界,具体如下:
重写公式(12)为如下形式
式中,αij和βij是安全域平面方程系数,求解方法为
线路潮流约束是线路相角差函数,变换公式(14)得到
将安全约束限值Vmax、Vmin、-δ、δ代入公式(14),便可得到近似真实安全域的子边界的超平面组。
下面以图2、图3和图4所示为例介绍近似真实安全域的超平面组,具体如下:
如图2所示的IEEE14节点标准电力网络,该IEEE14节点标准电力网络共有5台发电机和11个负荷,节点1为平衡节点,安全约束为Vi∈[0.9,1.1],θij∈[-10°,10°],其中i,j∈(2,3,…,14)。图3中(a)给出了考虑节点3和节点8的电压幅值约束子安全域超平面组,图3中(b)给出了考虑节点3和节点8的线路潮流约束子安全域超平面组,图4给出了考虑节点3和节点8的有功功率注入空间安全域超平面组。
步骤5,电力网络解裂。
下面以图5为例介绍一种是用于求取线路潮流安全域的电力网络解裂方法。对后续步骤进行具体说明:
***如图5中(a)所示的两区域互联电力网络后,可得到如图5中(b)所示的子网络-A和子网络-B。其中,子网络-A中***节点为节点1B’,2B’,和3B’,子网络-B中***节点为节点1A’,2A’和3A’。
步骤6,求取电力网络断面潮流静态安全域边界,具体如下:
对步骤5中的每个子电力网络而言,其***节点的功率注入安全域是只考虑该子电力网络安全约束时的安全域,因此,原电力网络断面潮流安全域等价于***后两个子电力网络的***节点的功率注入安全域的交集。
下面图6和图7为例介绍线路潮流安全域的求取方法及其超平面组,对后续步骤进行具体说明。
如图6所示的IEEE118节点标准电力网络中包含由线路15-33、19-34、30-38和23-24组成的断面,将该******成由线路15-33’、19-34’、30-38’、23-24’和区域-A构成的子电力网络-A、由线路15’-33、19’-34、30’-38、23’-24和区域-B构成的子电力网络-B,求取每个子电力网络的节点功率注入安全域。节点33’,34’,38’和节点24’的功率注入构成的子安全域即为原电力网络断面潮流只考虑区域-A和断面的安全约束下的安全域,节点15’,19’,30’和23’的功率注入构成的子安全域即为原电力网络断面潮流只考虑区域-B和断面的安全约束下的安全域。线路15-33,19-34,30-38和23-24组成的断面潮流安全域为由节点33’,34’,38’和节点24’的功率注入构成的子安全域和由节点15’,19’,30’和23’的功率注入构成的子安全域的交集。如图7所示,给出了线路23-24和线路30-38潮流的子安全域的可视化图形。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (4)
1.一种求解电力网络静态安全域的方法,其特征在于,首先通过线性化电力网络节点功率注入方程,建立节点电压幅值与相角和节点功率注入之间的线性关系,得到代表静态安全域边界的超平面组;然后通过***电力网络,将求解断面潮流安全域转化为求解***节点安全域,间接地获得了断面潮流安全域边界超平面组;具体实施步骤如下:
步骤1,输入电力网络数据;
步骤2,根据电力网络的拓扑结构,建立节点功率注入方程,即节点有功、节点无功注入方程;
步骤3,在理想运行点的小邻域内,线性化步骤2中建立的节点功率注入方程,得到节点电压幅值与相角和节点功率注入之间的线性关系;
步骤4,将静态安全约束数据代入步骤3中的节点电压幅值与相角和节点功率注入之间的线性关系方程,得到由线性方程表达的电力网络节点功率注入静态安全域边界;
步骤5,将电力网络从断面***开,***后所得的两个子电力网络都保留联络线和***节点;
步骤6,求解步骤5中每个子电力网络的***节点的功率注入安全域,则原电力网络断面潮流安全域等价于两个子网络的***节点的功率注入安全域的交集。
2.根据权利要求1所述的一种求解电力网络静态安全域的方法,其特征在于,步骤1中所述电力网络数据包括电力网络的结构数据、电气数据、静态安全约束数据;所述结构数据包括电力网络中线路的连接关系、线路与变压器的连接关系、发电机和负荷的位置;所述电气数据包括线路及变压器的电阻和电抗、发电机节点和负荷节点的功率注入;所述静态安全约束数据包括节点电压幅值安全约束、线路潮流安全约束。
3.根据权利要求1所述的一种求解电力网络静态安全域的方法,其特征在于,步骤3中在理想运行点的小邻域内,对电力网络节点功率注入方程在进行线性化,进而可采用超平面组来近似节点功率注入静态安全域边界。
4.根据权利要求1所述的一种求解电力网络静态安全域的方法,其特征在于,步骤3中所述理想运行点是(V=1,θ=0),其中,V是节点电压幅值向量;θ是节点电压相角向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410201937.5A CN104008279B (zh) | 2014-05-13 | 2014-05-13 | 一种求解电力网络静态安全域的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410201937.5A CN104008279B (zh) | 2014-05-13 | 2014-05-13 | 一种求解电力网络静态安全域的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104008279A true CN104008279A (zh) | 2014-08-27 |
CN104008279B CN104008279B (zh) | 2017-03-22 |
Family
ID=51368931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410201937.5A Expired - Fee Related CN104008279B (zh) | 2014-05-13 | 2014-05-13 | 一种求解电力网络静态安全域的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104008279B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105514999A (zh) * | 2016-01-29 | 2016-04-20 | 江苏省电力公司电力经济技术研究院 | 一种电网静态电压稳定预防控制模型及其算法 |
CN107039981A (zh) * | 2017-04-20 | 2017-08-11 | 河海大学 | 一种拟直流线性化概率最优潮流计算方法 |
CN107947199A (zh) * | 2017-12-01 | 2018-04-20 | 东北电力大学 | 基于优化模型的电力***热稳定安全域边界的搜索方法 |
CN108493921A (zh) * | 2018-01-19 | 2018-09-04 | 中国南方电网有限责任公司 | 基于关键节点注入功率空间的热稳定安全域快速构建方法 |
CN108521128A (zh) * | 2018-04-25 | 2018-09-11 | 东北电力大学 | 电力***静态电压安全域边界的快速搜索方法 |
CN111416348A (zh) * | 2020-04-21 | 2020-07-14 | 竺炜 | 全网静态电压安全及功角稳定的线路负荷安全域构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102983573B (zh) * | 2012-11-09 | 2014-10-15 | 天津大学 | 基于安全域的安全约束经济调度方法 |
-
2014
- 2014-05-13 CN CN201410201937.5A patent/CN104008279B/zh not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
哈比比等: "基于安全域的电力***有功及无功优化", 《中国电机工程学报》 * |
李晓珺: "基于静态安全域的大型互联电网在线安全监视方法研究", 《中国优秀硕士论文全文数据库》 * |
罗春雷等: "人工神经网络在电力***静态安全域分析中的应用", 《人工神经网络在电力***静态安全域分析中的应用》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105514999A (zh) * | 2016-01-29 | 2016-04-20 | 江苏省电力公司电力经济技术研究院 | 一种电网静态电压稳定预防控制模型及其算法 |
CN107039981A (zh) * | 2017-04-20 | 2017-08-11 | 河海大学 | 一种拟直流线性化概率最优潮流计算方法 |
CN107947199A (zh) * | 2017-12-01 | 2018-04-20 | 东北电力大学 | 基于优化模型的电力***热稳定安全域边界的搜索方法 |
CN108493921A (zh) * | 2018-01-19 | 2018-09-04 | 中国南方电网有限责任公司 | 基于关键节点注入功率空间的热稳定安全域快速构建方法 |
CN108493921B (zh) * | 2018-01-19 | 2021-04-23 | 中国南方电网有限责任公司 | 基于关键节点注入功率空间的热稳定安全域快速构建方法 |
CN108521128A (zh) * | 2018-04-25 | 2018-09-11 | 东北电力大学 | 电力***静态电压安全域边界的快速搜索方法 |
WO2019205222A1 (zh) * | 2018-04-25 | 2019-10-31 | 东北电力大学 | 电力***静态电压安全域边界的快速搜索方法 |
CN108521128B (zh) * | 2018-04-25 | 2020-04-28 | 东北电力大学 | 电力***静态电压安全域边界的快速搜索方法 |
CN111416348A (zh) * | 2020-04-21 | 2020-07-14 | 竺炜 | 全网静态电压安全及功角稳定的线路负荷安全域构建方法 |
CN111416348B (zh) * | 2020-04-21 | 2024-01-16 | 竺炜 | 全网静态电压安全及功角稳定的线路负荷安全域构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104008279B (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104008279A (zh) | 一种求解电力网络静态安全域的方法 | |
Medicherla et al. | Generation rescheduling and load shedding to alleviate line overloads-analysis | |
CN101860037B (zh) | 一种高压直流输电***网侧谐波电流的确定方法 | |
CN112531694A (zh) | 基于数字孪生技术的交直流混合电网全域实时模拟方法 | |
CN103973203B (zh) | 一种适用于安全稳定分析的大型光伏电站在线等值建模方法 | |
CN102831315B (zh) | 一种电力***状态估计量测方程的精确线性化方法 | |
CN101741092A (zh) | 电力***电压无功分区控制方法 | |
CN103326392B (zh) | 一种特高压直流输电换流阀组***可靠性计算方法 | |
CN105006839B (zh) | 并网发电***的弱电网网源荷阻抗模型分析方法 | |
WO2016061745A1 (zh) | 交流电力网的最小相位线性有功潮流的获取方法 | |
CN103714387A (zh) | 一种精细化电网风险评估方法 | |
CN101242101B (zh) | 电力***负荷谐波稳态模型的建立方法及仿真方法 | |
CN104104081A (zh) | 一种基于优化方法的非迭代不确定潮流分析方法 | |
Wang et al. | Resilience enhancement strategy using microgrids in distribution network | |
CN104993491A (zh) | 一种计及电压和无功的线性化潮流计算方法 | |
Liu et al. | Fast power system dynamic simulation using continued fractions | |
CN103618325B (zh) | 计及电力***运行方式变化的广域阻尼控制器设计方法 | |
CN106877309A (zh) | 主从控制模式下直流微电网的稳定性判定方法 | |
CN105790308A (zh) | 一种太阳能光伏发电***并网运行的控制方案 | |
Li et al. | Power system reliability analysis system based on PSASP and fault enumeration method and applications | |
CN105720579A (zh) | 一种基于lmi时滞电力***动态输出反馈控制器 | |
CN104600743A (zh) | 考虑风电集群功率波动的***关键变量提取方法 | |
CN110048428A (zh) | 基于概率守恒原理的电力***概率潮流计算方法 | |
Yesil et al. | Preliminary Studies on Dynamic Reduction of the Turkish Transmission Network | |
Wang et al. | Simulation case studies on period optimal power flow |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170322 |