CN102522742A - 基于单点量测信息的外网戴维南等值参数的估计方法 - Google Patents

基于单点量测信息的外网戴维南等值参数的估计方法 Download PDF

Info

Publication number
CN102522742A
CN102522742A CN2011103278890A CN201110327889A CN102522742A CN 102522742 A CN102522742 A CN 102522742A CN 2011103278890 A CN2011103278890 A CN 2011103278890A CN 201110327889 A CN201110327889 A CN 201110327889A CN 102522742 A CN102522742 A CN 102522742A
Authority
CN
China
Prior art keywords
estimated value
equivalent
line end
dai weinan
built
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
CN2011103278890A
Other languages
English (en)
Other versions
CN102522742B (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.)
Chongqing University
Original Assignee
Chongqing 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 Chongqing University filed Critical Chongqing University
Priority to CN201110327889.0A priority Critical patent/CN102522742B/zh
Publication of CN102522742A publication Critical patent/CN102522742A/zh
Application granted granted Critical
Publication of CN102522742B publication Critical patent/CN102522742B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

一种基于单点量测信息的外网戴维南等值参数的估计方法,涉及电力系统外网戴维南等值参数估计方法技术领域。本发明利用计算机和线路末端的数据采集装置,通过程序,首先输入被监测线路的参数,然后根据被监测线路末端多个时段的SCADA量测数据,建立以残差平方和最小为目标的带约束的最小二乘等值模型,求解该模型,得出被监测线路以外网络的等值内电势和等值阻抗参数。本发明充分考虑了外网戴维南等值参数的物理意义,能够有效保证外网等值参数的合理性,能够有效避免量测误差对等值参数计算结果的影响,准确性高,方法简单,估计精度高,计算速度快,可广泛应用于任何末端安装有数据采集装置的电力线路外网戴维南等值参数的估计。

Description

基于单点量测信息的外网戴维南等值参数的估计方法
技术领域
本发明属于电力系统外网等值参数估计方法技术领域,具体涉及电力系统外网戴维南等值参数的估计方法。
背景技术
现代电力系统已经发展成为大规模、多区域的互联电网。由于对整个电网进行分析的计算量大,计算时间长,而且同步获取全网所有实时测量信息也有困难,电力系统在线计算和分析中通常需要将外网进行简化等值,其中,戴维南等值是针对单端外网的一种最常用的等值方法。外网等值参数的准确性直接关系到被研究局部网络各种计算和分析的准确性,特别是会对电压稳定等在线控制与决策的有效性产生重要影响,对于保证电力系统的安全稳定运行具有重大意义。
现有的电力系统外网戴维南等值参数的估计方法,如2009年第33卷第10期《电力系统自动化》中“基于PMU和改进戴维南等值模型的电压稳定在线监视”一文,公开的方法是假设戴维南等值电路中内电势幅值和等值阻抗恒定,考虑等值内电势相角的变化,根据被监测节点三个时间断面的量测数据,建立等值电源点与被监测节点间的电压方程,进而解得外网的戴维南等值参数。该方法的主要缺点是:不能避免量测误差对等值参数计算结果的影响,等值精度低。又如2011年第44卷第4期《中国电力》中“基于最小二乘外网等值的线路电压稳定性指标”一文,公开的方法是根据被监测线路末端节点三个或三个以上时段的SCADA(Supervisory Control and Data Acquisition,监控与数据采集装置)量测数据,建立以被监测线路末端节点电压幅值、线路有功功率和无功功率的残差平方和最小为目标的最小二乘模型,并用牛顿法求解得到戴维南等值电路中内电势幅值、等值电阻、等值电抗和每个时段的等值内电势相角。该方法的主要缺点是:没有考虑等值参数的物理意义,不能避免等值参数中出现电阻为负值等不合理现象,其等值结果不能直接应用于电力系统电压稳定的在线分析与控制决策。
发明内容
本发明的目的是针对现有电力系统外网戴维南等值参数估计方法的不足,提供一种基于单点量测信息的外网戴维南等值参数估计方法。本发明方法充分考虑了外网戴维南等值参数的物理意义,能够有效保证外网等值参数的合理性,能够有效避免量测误差对等值参数计算结果的影响,等值效果好,其外网等值结果能够直接应用于电力系统电压稳定性的在线分析与控制决策。
实现本发明目的之技术方案是:一种基于单点量测信息的外网戴维南等值参数的估计方法,利用计算机和被监测线路末端的数据采集装置,通过程序,首先输入被监测线路的电阻、电抗和对地电纳参数,然后根据被监测线路末端多个时段的SCADA量测数据,建立以残差平方和最小为目标的带不等式约束的最小二乘等值模型,求解该模型,得出被监测线路以外网络的等值内电势和等值阻抗参数。其具体方法步骤如下:
(1)   输入基础数据
首先输入被监测线路的基础数据和收敛精度ε及时段数n,要求n≥3。所述的线路的基础数据为线路电阻、电抗和电纳参数,线路的额定电压、功率基准以及线路末端n个时段的SCADA量测数据,即线路末端电压的幅值、末端电流的幅值以及线路末端的有功功率和无功功率;
(2) 计算各时段电压和电流量测量的实部和虚部
第(1)步完成后,根据被监测线路末端各时段的SCADA量测数据,以末端节点电压为参考相量,计算各时段被监测线路末端电压和电流量测量的实部和虚部,计算公式为:
Figure 137478DEST_PATH_IMAGE001
式中,
Figure 992302DEST_PATH_IMAGE002
Figure 772039DEST_PATH_IMAGE003
分别为第m个时段线路末端节点电压幅值、线路末端电流幅值、线路末端有功功率和无功功率的量测量,m=1,2,…,n;取末端节点电压相量的相角为0°,
Figure 167248DEST_PATH_IMAGE004
Figure 919303DEST_PATH_IMAGE005
分别为第m个时段线路末端节点电压量测量的实部和虚部,
Figure 741766DEST_PATH_IMAGE006
Figure 8799DEST_PATH_IMAGE007
分别为第m个时段线路末端电流量测量的实部和虚部,tg-1表示反正切函数;
  (3)  计算外网戴维南等值参数的估计值   
第(2)步完成后,按下式建立外网戴维南等值参数的最小二乘估计模型:
Figure 145382DEST_PATH_IMAGE008
式中:
Figure 814261DEST_PATH_IMAGE009
为估计模型的变量,其中,E eq R eq X eq 分别为戴维南等值电路中内电势幅值、等值电阻和等值电抗的估计值,
Figure 807625DEST_PATH_IMAGE010
分别为第1,2,…,n个时段戴维南等值电路中内电势相角的估计值,
Figure 496708DEST_PATH_IMAGE011
分别为第1,2,…,n个时段线路末端电压实部的估计值,
Figure 233720DEST_PATH_IMAGE012
分别为第1,2,…,n个时段线路末端电流实部的估计值,
Figure 757105DEST_PATH_IMAGE013
别为第1,2,…,n个时段线路末端电流虚部的估计值;
Figure 124632DEST_PATH_IMAGE014
Figure 366258DEST_PATH_IMAGE015
分别为第m个时段线路末端电压量测量的实部、电流量测量的实部和电流量测量的虚部;
Figure 222535DEST_PATH_IMAGE017
分别为第m个时段线路末端电压实部的估计值、电流实部的估计值和电流虚部的估计值;λ (m)ρ (m)分别为第m个时段虚拟量测量的估计值,其表达式为:
式中:RXb分别为被监测线路的电阻、电抗和对地电纳的参考值;E eq R eq X eq 分别为戴维南等值电路中内电势幅值、等值电阻和等值电抗的估计值;
Figure 224306DEST_PATH_IMAGE019
为第m个时段戴维南等值电路中内电势相角的估计值,
Figure 801098DEST_PATH_IMAGE017
分别为第m个时段线路末端电压实部的估计值、电流实部的估计值和电流虚部的估计值;
求解式(5)-(8)所示的最小二乘估计模型,得到外网戴维南等值参数的估计值,即得到戴维南等值内电势幅值的估计值E eq 、等值电阻的估计值R eq 、等值电抗的估计值X eq 和各时段等值内电势相角的估计值
Figure 510428DEST_PATH_IMAGE010
(4) 输出外网戴维南等值参数的估计值
 第(3)步完成后,输出外网戴维南等值参数的估计值,即输出戴维南等值内电势幅值的估计值E eq 、等值电阻的估计值R eq 、等值电抗的估计值X eq 和各时段等值内电势相角的估计值
Figure 461067DEST_PATH_IMAGE010
本发明采用上述技术方案后,主要有以下效果:
1. 本发明方法考虑了戴维南等值参数的物理意义和应该满足的约束条件,能够有效保证外网等值参数估计结果的合理性,等值效果好,等值结果能够直接应用于电力系统电压稳定性的在线分析与控制决策;
2. 本发明方法充分利用了最小二乘法的滤波作用,能够有效降低量测误差对等值参数计算结果的影响,等值精度高;
3. 本发明方法利用线路末端的SCADA量测数据进行外网等值参数的估计,对基础数据的要求不高,方法简单,计算速度快,实用性强。
本发明可广泛应用于任何末端安装有数据采集装置的电力线路的外网等值参数的估计。
附图说明
图1为本发明方法的程序流程框图;
图2为应用本发明方法进行外网戴维南等值参数估计的系统接线图。
图中:ij分别为被监测线路首端节点和末端节点的节点号;R、Xb分别为被监测线路的电阻、电抗和对地电纳;UI分别线路末端节点的电压幅值和电流幅值,PQ分别为线路末端的有功功率和无功功率。E eq R eq X eq 分别为被监测线路外网戴维南等值电路中内电势的幅值、等值电阻和等值电抗,
Figure 343572DEST_PATH_IMAGE019
为第m个时段等值内电势的相角。
具体实施方式
下面结合具体实施方式,进一步说明本发明。
实施例
如图1、2所示,一种基于单点量测信息的外网戴维南等值参数的估计方法的具体步骤如下:
(1) 输入基础数据
首先输入被监测线路的基础数据和收敛精度ε及时段数。线路的基础数据为线路电阻、电抗和电纳参数,线路的额定电压、功率基准以及线路末端n个时段的SCADA量测数据,即线路末端电压的幅值、末端电流的幅值以及末端的有功功率和无功功率;
如图2所示的被监测电力线路的参数为:电阻R= 1.1297Ω,电抗X= 3.3471Ω,对地电纳b = 2.45306×10-5 S;线路的额定电压为U B  = 110 kV,功率基准S = 100 MVA。时段数n=10;收敛精度ε= 0.00001;
以第一个时段为例,线路末端的SCADA量测数据如下表所示:
(2) 计算各时段电压和电流量测量的实部和虚部
第(1)步完成后,根据被监测线路末端各时段的SCADA量测数据,以末端节点电压为参考相量,计算各时段被监测线路末端电压和电流量测量的实部和虚部,计算公式为技术方案中的公式(1)-(4)。
对第一个时段举例,按技术方案的公式(1)-(4),可得:
Figure 576287DEST_PATH_IMAGE021
即计算得到的线路末端节点电压量测量的实部为
Figure 14222DEST_PATH_IMAGE022
(p.u.表示标么值);末端节点电压量测量的虚部为
Figure 372522DEST_PATH_IMAGE023
,末端电流量测量的实部为
Figure 845092DEST_PATH_IMAGE024
,末端电流量测量的虚部为
Figure 958541DEST_PATH_IMAGE025
(3) 计算外网戴维南等值参数的估计值
第(2)步完成后,按技术方案中的公式(5)-(8),建立外网戴维南等值参数的最小二乘估计模型,即:
Figure 821455DEST_PATH_IMAGE026
式中:
Figure 45763DEST_PATH_IMAGE027
求解该最小二乘估计模型,得到外网戴维南等值参数的估计值;
(4) 输出外网戴维南等值参数的估计值
第(3)步完成后,输出外网戴维南等值参数的估计值,即输出戴维南等值内电势幅值的估计值E eq 、等值电阻的估计值R eq 、等值电抗的估计值X eq 和各时段等值内电势相角的估计值
Figure 107260DEST_PATH_IMAGE010
附图2中的线路,外网戴维南等值参数的计算结果为:E eq =1.062p.u.,R eq =0.795×10-4p.u.,X eq =2.817×10-4p.u.,每个时段的等值内电势相角如下表所示:
Figure 349802DEST_PATH_IMAGE028
实验效果
以本实施例的被监测线路(如图2所示)为对象,通过仿真算例验证本发明方法的有效性。
当被监测电力线路的参数为:电阻R= 1.1297 Ω,电抗X= 3.3471 Ω,对地电纳b = 2.45306×10-5 S,分别应用本发明方法和现有技术中的第三种方法,估计被监测线路的外网戴维南等值参数,其中5个时段的估计结果如下表所示:
Figure 762328DEST_PATH_IMAGE030
从上述实验结果可知:应用现有技术估计外网的戴维南等值参数时,会出现电阻为负、电抗为负或输电网电阻远大于电抗三种不合理结果,而应用本发明方法,则能够有效保证等值参数估计值的合理性。

Claims (1)

1.一种基于单点量测信息的外网戴维南等值参数的估计方法,利用计算机和被监测线路末端的数据采集装置,通过程序,计算被监测线路以外网络的等值内电势和等值阻抗参数,其特征在于其具体的方法步骤如下:
(1)   输入基础数据
首先输入被监测线路的基础数据和收敛精度ε及时段数n,要求n≥3;所述的线路的基础数据为线路电阻、电抗和电纳参数,线路的额定电压、功率基准以及线路末端n个时段的SCADA量测数据,即线路末端电压的幅值、末端电流的幅值以及线路末端的有功功率和无功功率;
(2) 计算各时段电压和电流量测量的实部和虚部
第(1)步完成后,根据被监测线路末端各时段的SCADA量测数据,以末端节点电压为参考相量,计算各时段被监测线路末端电压和电流量测量的实部和虚部,计算公式为:
Figure 774770DEST_PATH_IMAGE001
式中,
Figure 127254DEST_PATH_IMAGE002
Figure 69802DEST_PATH_IMAGE003
分别为第m个时段线路末端节点电压幅值、线路末端电流幅值、线路末端有功功率和无功功率的量测量,m=1,2,…,n;取末端节点电压相量的相角为0°,
Figure 445420DEST_PATH_IMAGE004
分别为第m个时段线路末端节点电压量测量的实部和虚部,
Figure 327422DEST_PATH_IMAGE007
分别为第m个时段线路末端电流量测量的实部和虚部,tg-1表示反正切函数;
  (3)  计算外网戴维南等值参数的估计值   
第(2)步完成后,按下式建立外网戴维南等值参数的最小二乘估计模型:
Figure 252653DEST_PATH_IMAGE008
式中:
Figure 411715DEST_PATH_IMAGE009
为估计模型的变量,其中,E eq R eq X eq 分别为戴维南等值电路中内电势幅值、等值电阻和等值电抗的估计值,
Figure 473211DEST_PATH_IMAGE010
分别为第1,2,…,n个时段戴维南等值电路中内电势相角的估计值,
Figure 757562DEST_PATH_IMAGE011
分别为第1,2,…,n个时段线路末端电压实部的估计值,
Figure 107772DEST_PATH_IMAGE012
分别为第1,2,…,n个时段线路末端电流实部的估计值,
Figure 135771DEST_PATH_IMAGE013
别为第1,2,…,n个时段线路末端电流虚部的估计值;
Figure 317354DEST_PATH_IMAGE014
分别为第m个时段线路末端电压量测量的实部、电流量测量的实部和电流量测量的虚部;
Figure 344533DEST_PATH_IMAGE016
Figure 910643DEST_PATH_IMAGE017
分别为第m个时段线路末端电压实部的估计值、电流实部的估计值和电流虚部的估计值;λ (m)ρ (m)分别为第m个时段虚拟量测量的估计值,其表达式为:
式中:RXb分别为被监测线路的电阻、电抗和对地电纳的参考值;E eq R eq X eq 分别为戴维南等值电路中内电势幅值、等值电阻和等值电抗的估计值;
Figure 510569DEST_PATH_IMAGE020
为第m个时段戴维南等值电路中内电势相角的估计值,
Figure 897688DEST_PATH_IMAGE016
Figure 205172DEST_PATH_IMAGE017
分别为第m个时段线路末端电压实部的估计值、电流实部的估计值和电流虚部的估计值;
求解式(5)-(8)所示的最小二乘估计模型,得到外网戴维南等值参数的估计值,即得到戴维南等值内电势幅值的估计值E eq 、等值电阻的估计值R eq 、等值电抗的估计值X eq 和各时段等值内电势相角的估计值
Figure 95768DEST_PATH_IMAGE010
(4) 输出外网戴维南等值参数的估计值
 第(3)步完成后,输出外网戴维南等值参数的估计值,即输出戴维南等值内电势幅值的估计值E eq 、等值电阻的估计值R eq 、等值电抗的估计值X eq 和各时段等值内电势相角的估计值
Figure 892823DEST_PATH_IMAGE010
CN201110327889.0A 2011-10-26 2011-10-26 基于单点量测信息的外网戴维南等值参数的估计方法 Expired - Fee Related CN102522742B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110327889.0A CN102522742B (zh) 2011-10-26 2011-10-26 基于单点量测信息的外网戴维南等值参数的估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110327889.0A CN102522742B (zh) 2011-10-26 2011-10-26 基于单点量测信息的外网戴维南等值参数的估计方法

Publications (2)

Publication Number Publication Date
CN102522742A true CN102522742A (zh) 2012-06-27
CN102522742B CN102522742B (zh) 2014-02-26

Family

ID=46293566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110327889.0A Expired - Fee Related CN102522742B (zh) 2011-10-26 2011-10-26 基于单点量测信息的外网戴维南等值参数的估计方法

Country Status (1)

Country Link
CN (1) CN102522742B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103023015A (zh) * 2012-11-07 2013-04-03 中国电力科学研究院 一种配电线路阻抗在线虚拟量测方法
CN103116097A (zh) * 2013-01-25 2013-05-22 中国电力科学研究院 基于多断面混合量测信息的设备参数在线辨识方法
CN103606921A (zh) * 2013-12-10 2014-02-26 上海交通大学 基于实时广义戴维南等值的暂态电压稳定判断方法
CN109188045A (zh) * 2018-06-08 2019-01-11 中国电力科学研究院有限公司 一种配电网线路电压的平差方法和装置
CN111181153A (zh) * 2019-11-20 2020-05-19 国网湖南省电力有限公司经济技术研究院 一种基于支路响应信息确定戴维南等值参数的方法及系统

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104917177A (zh) * 2015-06-26 2015-09-16 四川大学 一种电力系统戴维南等值参数在线识别方法及其装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1912304A2 (en) * 2006-10-09 2008-04-16 Electric Power Research Institute, Inc. Method for voltage instability load shedding using local measurements
CN101350524A (zh) * 2008-09-11 2009-01-21 中国电力科学研究院 一种基于全微分方程的戴维南等值参数跟踪的计算方法
CN102055191A (zh) * 2010-12-31 2011-05-11 重庆大学 一种考虑外网等值的静态电压稳定广域切负荷控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1912304A2 (en) * 2006-10-09 2008-04-16 Electric Power Research Institute, Inc. Method for voltage instability load shedding using local measurements
CN101350524A (zh) * 2008-09-11 2009-01-21 中国电力科学研究院 一种基于全微分方程的戴维南等值参数跟踪的计算方法
CN102055191A (zh) * 2010-12-31 2011-05-11 重庆大学 一种考虑外网等值的静态电压稳定广域切负荷控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
汤涌等: "基于时域仿真的戴维南等值参数跟踪计算方法", 《中国电机工程学报》 *
王漪等: "基于戴维南等值的系统参数跟踪估计", 《电网技术》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103023015A (zh) * 2012-11-07 2013-04-03 中国电力科学研究院 一种配电线路阻抗在线虚拟量测方法
CN103023015B (zh) * 2012-11-07 2014-12-10 中国电力科学研究院 一种配电线路阻抗在线虚拟量测方法
CN103116097A (zh) * 2013-01-25 2013-05-22 中国电力科学研究院 基于多断面混合量测信息的设备参数在线辨识方法
CN103116097B (zh) * 2013-01-25 2015-06-24 中国电力科学研究院 基于多断面混合量测信息的设备参数在线辨识方法
CN103606921A (zh) * 2013-12-10 2014-02-26 上海交通大学 基于实时广义戴维南等值的暂态电压稳定判断方法
CN103606921B (zh) * 2013-12-10 2015-08-19 上海交通大学 基于实时广义戴维南等值的暂态电压稳定判断方法
CN109188045A (zh) * 2018-06-08 2019-01-11 中国电力科学研究院有限公司 一种配电网线路电压的平差方法和装置
CN109188045B (zh) * 2018-06-08 2023-07-18 中国电力科学研究院有限公司 一种配电网线路电压的平差方法和装置
CN111181153A (zh) * 2019-11-20 2020-05-19 国网湖南省电力有限公司经济技术研究院 一种基于支路响应信息确定戴维南等值参数的方法及系统

Also Published As

Publication number Publication date
CN102522742B (zh) 2014-02-26

Similar Documents

Publication Publication Date Title
CN102495288B (zh) 基于scada量测数据的输电线路参数的辨识与估计方法
CN102522742B (zh) 基于单点量测信息的外网戴维南等值参数的估计方法
CN102175922B (zh) 基于pmu量测数据的电力线路参数的辨识与估计方法
CN106841778B (zh) 基于pmu实现的次同步和超同步谐波参数的处理方法
TWI479159B (zh) 電壓穩定度即時估測方法
CN101441244B (zh) 一种基于就地电气量测量失步振荡中心位置的判据
CN101750562A (zh) 基于潮流方程灵敏度分析的非pmu测点动态过程估计方法
CN102420427A (zh) 一种考虑外网等值的区域电网电压稳定裕度计算方法
CN107064630B (zh) 一种电力系统频率测量方法及装置
CN102611102A (zh) 一种基于pmu实测的电力系统分区建模方法
CN103020726A (zh) 面向全pmu量测的抗差状态估计方法
CN104242267A (zh) 一种风力发电送出输电线路距离保护方法
CN103872706A (zh) 一种基于同步相量测量技术的配网合环方法
CN103473602A (zh) 一种电网理论线损数据预测系统及预测方法
CN104052053A (zh) 基于一元线性回归模型的配电网双端线路段阻抗分析方法
CN108896852B (zh) 一种公共接入点短路容量在线测量方法及系统
Goh et al. Dynamic estimation of power system stability in different Kalman filter implementations
CN106872861A (zh) 超特高压输电工程线路直击导线的雷电流反演方法及系统
CN105629168B (zh) 一种基于功率区间的风电机组闪变确定方法
CN109858061A (zh) 面向电压功率灵敏度估计的配电网络等值与化简方法
CN105678640A (zh) 一种考虑变压器饱和影响的交流电网直流电流分布的预测方法
CN102004191A (zh) 消弧线圈并列运行系统中计算线路电容的方法
CN104198979A (zh) 一种用于宽频电能表计量运行误差的实时比对方法
CN104850751B (zh) 一种电流质量评估方法
CN105139286A (zh) 基于无功不平衡量的scada数据有效性评估方法

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

Granted publication date: 20140226

Termination date: 20151026

EXPY Termination of patent right or utility model