CN105677952B - 一种考虑失稳位置不确定的输电线路故障概率评估方法 - Google Patents
一种考虑失稳位置不确定的输电线路故障概率评估方法 Download PDFInfo
- Publication number
- CN105677952B CN105677952B CN201511023780.2A CN201511023780A CN105677952B CN 105677952 B CN105677952 B CN 105677952B CN 201511023780 A CN201511023780 A CN 201511023780A CN 105677952 B CN105677952 B CN 105677952B
- Authority
- CN
- China
- Prior art keywords
- probability
- route
- failure
- unstability
- fault
- 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.)
- Expired - Fee Related
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 12
- 238000011156 evaluation Methods 0.000 title claims abstract description 9
- 230000007257 malfunction Effects 0.000 claims abstract description 9
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 230000001052 transient effect Effects 0.000 claims description 33
- 230000008859 change Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims 1
- 240000002853 Nelumbo nucifera Species 0.000 claims 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 25
- 238000011438 discrete method Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 2
- 238000012854 evaluation process Methods 0.000 abstract 1
- 238000007796 conventional method Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012502 risk assessment Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003121 nonmonotonic effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- APTZNLHMIGJTEW-UHFFFAOYSA-N pyraflufen-ethyl Chemical compound C1=C(Cl)C(OCC(=O)OCC)=CC(C=2C(=C(OC(F)F)N(C)N=2)Cl)=C1F APTZNLHMIGJTEW-UHFFFAOYSA-N 0.000 description 1
- 238000013468 resource allocation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种输电线路故障概率精确评估方法。该方法针对故障概率评估过程中失稳位置不确定情况，对故障失稳区间动态评估，得出不同场景下失稳区间的分布，从而使得故障概率的评估结果更为精确。与传统二分法、分段法等相比，该方法避开了用数值仿真得到摇摆曲线后，评估系统稳定程度或分析系统稳定机理的困难性；针对各种工况都能够准确的求出线路的失稳部分和失稳占比。所得结果用于系统失稳概率的求取中能够提高准确性，也降低特殊工况运行时计算结果与实际结果相差较大的可能性。
Description
技术领域
本发明涉及一种输电线路故障概率计算方法，特别是一种考虑失稳位置不确定的输电线路故障概率评估方法。
背景技术
近年来，随着电力工业的不断发展壮大，电力系统的规模不断增加，互联电网的规模也越来越大。一旦发生大停电事故将会在长时间、广范围内造成很严重的后果和损失。例如，09年的巴西大停电导致巴西国土面积一半的范围陷入停电状态，造成巨大损失和危害。这种情况下，系统可靠性风险评估对于保障系统安全运行具有重要意义。
由于实际运行中的电力系统深受所处环境影响，自然灾害也是导致电力系统发生故障甚至失稳崩溃的重要原因，通常对电力系统进行可靠性风险评估时不仅要考虑系统运行工况，还要结合外部环境因素。
现阶段大量使用的电力系统安全评估方法有传统的确定性方法和不确定性方法，确定性方法考虑在一定拓扑状态下会导致最严重系统故障的表现行为，并以此来确定系统的安全水平。虽然此方法能提供较高的可靠等级，但过于保守，使用过程势必造成设备的不充分利用或系统的不必要重复建设以及资源分配不合理等较大使用代价，同时不能全面考虑环境等诸多不确定因素。不确定性方法中利用风险的概念，从故障导致系统暂态失稳的概率和失稳造成的损失后果两个方面去分析去评估电力系统的安全性。对于故障导致系统暂态失稳概率部分，先基于全概率和条件概率理论建立暂态失稳模型，将系统发生暂态失稳的概率定义为求和系统各条线路上发生各类故障的概率和该故障发生导致系统失稳的概率之积。其中故障发生导致系统失稳的概率又由各类故障发生的概率和该类故障下线路的失稳区间在整条线路的占比之积求和得到。现有针对线路故障失稳占比的求取常用二分法和分段法等。分段法根据线路的长度将线路划分为多段，并以段内某点的稳定性反应该段的稳定性；二分法通过二分逼近原理确定线路失稳临界点来处理失稳区间问题。对于以上方法，当系统稳定性随故障位置变化出现不单调或异常时，或当故障发生导致线路上存在多段失稳域时，计算结果可能存在较大偏差；与此同时，当系统运行考虑外部灾害，例如雷击，其发生的不确定性导致的雷区分布不均可能会加重计算结果的偏差，导致最后评估结果不理想。并且，上述方法通过数值仿真得到摇摆曲线，评估系统稳定程度和分析系统稳定机理都具有一定的困难性。
由上可知，现有求取系统失稳概率方法在线路失稳占比和系统稳定裕度的确定及评价方面还存在一定缺陷。
发明内容
本发明的目的在于提供一种求取系统失稳概率中确定线路失稳占比方法。
实现本发明目的的技术解决方案为：一种考虑失稳位置不确定的输电线路故障概率评估方法，包括如下步骤：
步骤1、确定离散型随机变量和连续型随机变量的参数，所述离散变量的参数包括电网拓扑、线路参数、发电机参数，连续型变量的参数包括负荷水平、故障持续时间、故障位置；之后根据上述参数填写FASTEST潮流数据卡和稳定数据卡；
步骤2、确定故障类型和选择故障线路，生成FASTEST故障操作卡；所述故障类型包括：三相短路故障、双相接地故障、相间短路故障和单相接地故障；
步骤3、选择待计算线路和故障类型，根据步骤1和步骤2中的初始化数据和操作卡，利用FASTEST仿真计算得出系统稳定裕度随线路故障位置的变化曲线，确定该故障下线路失稳部分，得出该故障在线路上的失稳占比；
步骤4、重复步骤3，得出各线路在各种故障类型下的线路失稳占比；
步骤5、基于全概率和条件概率的理论建立暂态失稳概率模型，系统发生暂态失稳的概率为求和系统各条线路上发生故障的概率和该故障发生导致系统失稳的概率之积；
所述系统暂态失稳概率模型P(K)为：
其中：
K：系统发生暂态失稳；P(K)：系统发生暂态失稳的概率；Ki：系统由于线路i上的故障事故发生暂态失稳；P(Ki)：系统由于线路i上故障事故发生暂态失稳的概率；Kij：系统由于线路i上的j类故障事故失稳；Eij：在线路i上发生j类故障事故；P(Eij)为线路i上发生j类故障事故的概率，其公式为：
式中，k为tp时间段内元件的故障次数；λij为线路i上j类故障的故障率，系统状态和负荷水平在参数tp时间段内视为不变；
P(Kij|Eij)为故障导致系统暂态失稳的概率模型：
其中：lij为线路i上j类故障可以导致系统暂态失稳部分长度，Li为线路i的长度。
时间段tp的选取满足如下条件：
其中，0＜λijtp≤1，使得上式结果大于零，P(Eij)单调递增。
步骤6、整合上述计算数据，带入步骤5的模型求取输电线路故障概率。
本发明与现有方法相比，其显著优点为：本发明通过利用FASTEST，在确定系统参数的情况下，通过填卡仿真得出系统稳定裕度随线路上故障位置的变化曲线，统计得出线路失稳占比；与分段法和二分法相比，可以直观判断系统的稳定程度随故障位置的变化趋势，简化复杂的数值仿真计算，同时针对稳定裕度随故障位置非单调变化情况导致的线路多段失稳域现象，能够解决传统方法计算结果误差较大问题，得到准确结果；当结合外部环境因素考虑系统暂态稳定性时，对于环境因素的不确定性造成的线路分段故障率同时存在问题，也能很好的适用，给出准确的失稳故障概率结果。
附图说明
图1为本发明的考虑失稳位置不确定的输电线路故障概率评估流程图。
图2为计算6-11线路上发生故障系统暂态失稳概率线路状态图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
如图1所示，一种考虑失稳位置不确定的输电线路故障概率评估方法，包括以下步骤：
步骤1、确定离散型随机变量和连续型随机变量参数，包括负荷水平、电网拓扑、发电机参数、线路参数、故障持续时间和故障位置，根据参数初始化FASTEST潮流和稳定数据卡；
步骤2、选择待计算线路和故障类型，取线路1％占比为精度单位，从首端到末端分别放置100个故障位置点，逐一模拟系统暂态稳定性，生成FASTEST故障操作卡；
步骤3、根据步骤1和步骤2中的初始化数据和操作卡，利用FASTEST仿真计算得出系统稳定裕度随线路故障位置的变化曲线，确定该故障下线路失稳部分，得出该故障在线路上的失稳占比；
步骤4、重复步骤3，遍历每条线路在各种故障类型下的线路失稳占比；
步骤5、基于全概率和条件概率的理论构建暂态失稳概率模型，将系统发生暂态失稳的概率定义为求和系统各条线路上发生故障的概率和该故障发生导致系统失稳的概率之积；
所述系统各条线路上发生故障的概率由所研究时段内线路的故障率得到，采用泊松分布模型，表达如下式：
式中，k为tp时间段内元件的故障次数；λ为元件的故障率。
当线路受外部环境影响时，可以通过故障率的大小变化来直观的反应，即λ根据外部环境的同步变化。在不考虑可修复性情况下，一般我们认为设备元件发生一次故障的概率即为该设备元件的故障概率，则P(E)表达式如下：
所述故障发生导致系统失稳的概率由各类故障发生的概率和该类故障下线路的失稳区间在整条线路的占比之积求和得到，本文考虑的故障类型为三相短路故障、两相接地故障、相间短路故障以及单相接地故障，其他类型的故障由于发生的概率很小或者造成的影响很小或两者兼有而忽略不计。
根据贝叶斯——拉普拉斯准则，采用统一分布，认为线路上每一个位置发生故障的概率是公平的，则故障在线路上出现的密度分布函数为：
则线路i上的j类故障引发系统暂态失稳的概率为：
lij为线路i上j类故障可以导致系统暂态失稳部分长度，Li为线路i的长度。
各类型故障发生的频率为fj，P(j|Ei)为线路i上的故障为j类故障的概率，则有下式：
利用全概率原理，则系统暂态失稳概率为：
步骤6、整合仿真结果，带入步骤5中的暂态失稳概率模型中求得该时段下暂态失稳的概率。
为了进一步对本方法进行说明，以分析10机39节点系统故障失稳概率为例，步骤如下：
步骤一：确定负荷水平、电网拓扑、发电机参数以及线路参数如表一、二、三所示，选择考虑的故障类型包括三相短路故障、两相接地故障、相间短路故障以及单相接地故障，切除时间为0.14s，根据参数初始化潮流和稳定参数；
表1.IEEE39节点系统节点参数
表2.IEEE39节点系统支路数据
表3.IEEE39节点系统发电机参数
步骤二：选定故障发生线路及故障类型，取线路1％占比为精度单位，从首端到末端分别放置100个故障位置点，生产故障操作文件。
步骤三：根据步骤一和步骤二中的数据，利用FASTEST对该工况下系统潮流计算并仿真出系统稳定裕度随线路故障位置的变化曲线，确定该线路在选定故障类型下的失稳占比。
步骤四：重复步骤三，遍历每条线路在各种故障类型下的线路失稳占比。为了直观体现本文方法的优越性，在同样的系统参数、网络拓扑和负荷水平等条件下，运用二分法计算失稳占比，进行结果比较，取线路6-11为例，计算结果如表4所示：
表4失稳占比结果对比
| 线路 | FASTEST法失稳占比 | 二分法失稳占比 |
| 6-11 | 0.44 | 1 |
步骤五：构建暂态失稳概率模型；
步骤六：结合实际运行情况，考虑外部环境的影响，将线路失稳区间和雷区分布结合，假设此时线路6-11的BC段位于雷电中，雷电面积50km2，实际落雷数量为300，如图2所示，其中阴影部分为线路6-11故障失稳区间，则在此情况下，针对不同方法计算出该线路上发生故障时系统失稳概率对比如表5。
表5.失稳概率对比结果
| 线路 | FASTEST法 | 二分法 |
| 6-11 | 4.88×10<sup>-4</sup> | 4.64×10<sup>-2</sup> |
以上计算结果可以看出本方法的优越性，本发明的方法可以直观判断系统的稳定程度随故障位置的变化趋势，简化复杂的数值仿真计算，同时针对稳定裕度随故障位置非单调变化情况导致的线路多段失稳域现象，能够解决传统方法计算结果误差较大问题，得到准确结果。
Claims (2)
1.一种考虑失稳位置不确定的输电线路故障概率评估方法，其特征在于，包括如下步骤：
步骤1、确定离散型随机变量和连续型随机变量的参数，所述离散变量的参数包括电网拓扑、线路参数和发电机参数，连续型变量的参数包括负荷水平、故障持续时间和故障位置；之后根据上述参数填写FASTEST潮流数据卡和稳定数据卡；
步骤2、确定故障类型和选择故障线路，生成FASTEST故障操作卡；所述故障类型包括：三相短路故障、双相接地故障、相间短路故障和单相接地故障；
步骤3、选择待计算线路和故障类型，根据步骤1和步骤2中的初始化数据和操作卡，利用FASTEST仿真计算得出系统稳定裕度随线路故障位置的变化曲线，确定该故障下线路失稳部分，得出该故障在线路上的失稳占比；
步骤4、重复步骤3，得出各线路在各种故障类型下的线路失稳占比；
步骤5、基于全概率和条件概率的理论建立暂态失稳概率模型，系统发生暂态失稳的概率为求和系统各条线路上发生故障的概率和该故障发生导致系统失稳的概率之积；其中所述系统暂态失稳概率模型P(K)为：
其中：
K：系统发生暂态失稳；P(K)：系统发生暂态失稳的概率；Ki：系统由于线路i上的故障事故发生暂态失稳；P(Ki)：系统由于线路i上故障事故发生暂态失稳的概率；Kij：系统由于线路i上的j类故障事故失稳；Eij：在线路i上发生j类故障事故；P(Eij)为线路i上发生j类故障事故的概率，其公式为：
式中，k为tp时间段内元件的故障次数；λij为线路i上j类故障的故障率，系统状态和负荷水平在参数tp时间段内视为不变；
P(Kij|Eij)为故障导致系统暂态失稳的概率模型：
其中：lij为线路i上j类故障可以导致系统暂态失稳部分长度，Li为线路i的长度；
步骤6、整合上述计算数据，带入步骤5的模型求取输电线路故障概率。
2.根据权利要求1所述的考虑失稳位置不确定的输电线路故障概率评估方法，其特征在于，时间段tp的选取满足如下条件：
其中，0＜λijtp≤1，P(Eij)单调递增。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511023780.2A CN105677952B (zh) | 2015-12-30 | 2015-12-30 | 一种考虑失稳位置不确定的输电线路故障概率评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511023780.2A CN105677952B (zh) | 2015-12-30 | 2015-12-30 | 一种考虑失稳位置不确定的输电线路故障概率评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105677952A CN105677952A (zh) | 2016-06-15 |
| CN105677952B true CN105677952B (zh) | 2019-01-15 |
Family
ID=56189835
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511023780.2A Expired - Fee Related CN105677952B (zh) | 2015-12-30 | 2015-12-30 | 一种考虑失稳位置不确定的输电线路故障概率评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105677952B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107886220B (zh) * | 2017-10-23 | 2022-01-25 | 广西电网有限责任公司南宁供电局 | 一种基于历史因素分析的配电线路风险概率评估方法 |
| CN110309873B (zh) * | 2019-06-28 | 2021-05-14 | 清华大学 | 基于潮流向量贝叶斯分类的集成暂态稳定评估方法及系统 |
| CN111242459A (zh) * | 2020-01-07 | 2020-06-05 | 中国南方电网有限责任公司 | 全网设备参数异常值辨识方法及系统 |
| CN111950765B (zh) * | 2020-07-06 | 2024-04-19 | 四川大川云能科技有限公司 | 一种基于堆叠降噪自编码器的概率性暂态稳定预测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090319093A1 (en) * | 2008-03-31 | 2009-12-24 | The Royal Institution For Advancement Of Learning/ Mcgill University | Methods and processes relating to electricity power generation and distribution networks |
| CN102521667A (zh) * | 2011-12-26 | 2012-06-27 | 华北电力大学（保定） | 电力系统阶段式保护运行风险的概率评估方法 |
| CN104993464A (zh) * | 2015-07-15 | 2015-10-21 | 国网河南省电力公司漯河供电公司 | 一种电力系统继电保护故障分析方法 |
-
2015
- 2015-12-30 CN CN201511023780.2A patent/CN105677952B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090319093A1 (en) * | 2008-03-31 | 2009-12-24 | The Royal Institution For Advancement Of Learning/ Mcgill University | Methods and processes relating to electricity power generation and distribution networks |
| CN102521667A (zh) * | 2011-12-26 | 2012-06-27 | 华北电力大学（保定） | 电力系统阶段式保护运行风险的概率评估方法 |
| CN104993464A (zh) * | 2015-07-15 | 2015-10-21 | 国网河南省电力公司漯河供电公司 | 一种电力系统继电保护故障分析方法 |
Non-Patent Citations (4)
| Title |
|---|
| Probabilistic transient stability assessment for on-line applications;Mehdi Abapour等;《International Journal of Electrical Power & Energy Systems》;20120606;第42卷(第1期);第627-634页 |
| 冰灾对输电线故障率影响的时空评估;谢云云等;《电力系统自动化》;20130925;第37卷(第18期);第32-42页 |
| 基于稳定域及条件概率的暂态稳定不确定性分析;刘强等;《电力系统自动化》;20071010;第31卷(第19期);第1-6页 |
| 电网雷击故障概率的时空在线预警;谢云云等;《电力系统自动化》;20130910;第37卷(第7期);第44-51页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105677952A (zh) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105677952B (zh) | 一种考虑失稳位置不确定的输电线路故障概率评估方法 | |
| CN106709651B (zh) | 一种基于风险理论的电力系统安全性评估系统 | |
| CN105976257A (zh) | 基于隶属度函数的模糊综合评价法的电网脆弱性评估方法 | |
| CN105098763B (zh) | 一种基于wams和ems的暂态电压稳定在线评估方法 | |
| CN111327047B (zh) | 一种确定级联式电网多交流断面输电能力的方法及系统 | |
| CN109242239A (zh) | 一种交直流混合电网运行风险评估方法和装置 | |
| CN107316113A (zh) | 一种输电网规划方法及系统 | |
| CN103995982A (zh) | 一种考虑机组随机故障的概率潮流计算方法 | |
| CN109494719A (zh) | 一种中低压混合配电网层次化阻抗分析方法 | |
| CN109919398A (zh) | 基于图分割算法的含风电电力系统分区备用优化配置方法 | |
| CN103473602A (zh) | 一种电网理论线损数据预测系统及预测方法 | |
| CN115062534A (zh) | 天然气管道系统的供气可靠性计算方法及装置 | |
| CN104934982B (zh) | 一种多直流馈入系统的直流分区方法 | |
| CN107147108B (zh) | 一种配电系统安全域维度的计算方法 | |
| CN114996635A (zh) | 一种配电网参数确定方法、装置、设备及存储介质 | |
| CN108768748B (zh) | 一种用于电力通信服务的故障诊断方法、装置及存储介质 | |
| CN106410862B (zh) | 基于有功恢复斜率校正的风电场单机等值方法 | |
| CN110348540A (zh) | 基于聚类的电力系统暂态功角稳定故障筛选方法及装置 | |
| Zhou et al. | Robustness analysis of power system dynamic process and repair strategy | |
| Papadopoulos et al. | Probabilistic assessment of transient stability in reduced inertia systems | |
| CN107146166B (zh) | 一种配网规划辅助系统 | |
| CN108877969A (zh) | 核功率理论模型建立及验证方法、系统及终端设备 | |
| CN108183481B (zh) | 一种基于深度学习电网快速判稳方法和系统 | |
| Alzubaidi et al. | Probabilistic voltage stability assessment considering load and wind uncertainties | |
| Zhu et al. | Efficient identification of critical load model parameters affecting power system voltage stability |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190115 Termination date: 20191230 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |