CN104200032B - 广义负荷建模中基于时段性的横向时间轴聚类方法 - Google Patents
广义负荷建模中基于时段性的横向时间轴聚类方法 Download PDFInfo
- Publication number
- CN104200032B CN104200032B CN201410453541.XA CN201410453541A CN104200032B CN 104200032 B CN104200032 B CN 104200032B CN 201410453541 A CN201410453541 A CN 201410453541A CN 104200032 B CN104200032 B CN 104200032B
- Authority
- CN
- China
- Prior art keywords
- data
- cluster
- power
- horizontal time
- period
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 239000013598 vector Substances 0.000 claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims abstract description 25
- 230000005611 electricity Effects 0.000 claims abstract description 24
- 239000000203 mixture Substances 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 5
- 230000008859 change Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 230000007935 neutral effect Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 12
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000004422 calculation algorithm Methods 0.000 description 16
- 238000004364 calculation method Methods 0.000 description 8
- 230000010354 integration Effects 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 4
- 206010008190 Cerebrovascular accident Diseases 0.000 description 3
- 241001269238 Data Species 0.000 description 3
- 208000006011 Stroke Diseases 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000002372 labelling Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 206010068052 Mosaicism Diseases 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 238000013481 data capture Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 210000003765 sex chromosome Anatomy 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- 108091092195 Intron Proteins 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000010189 synthetic method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了广义负荷建模中基于时段性的横向时间轴聚类方法,获取全年风电与负荷组成的根母线数据;数据处理,将处理后的全部数据首尾相连,横列组成横向连续数据,按横向时间单元THN分为M段,对全部数据进行横向聚类;基于横向聚类结果,待分析样本数据源被划分为q个横向类,每类由各自聚类中心代表;特征向量交叉匹配。通过特征向量交叉匹配判断样本所属类别,利用广义负荷建模建立精确模型并检验聚类策略有效性,仿真表明经过聚类分析后的广义负荷建模,可在满足精确性的基础上便于模型走向实用化,有利于提高电力系统仿真的精确性与实效性。
Description
技术领域
本发明涉及广义负荷建模中基于时段性的横向时间轴聚类方法。
背景技术
作为一种间歇性能源,风电以其随机性和波动性给电力系统安全稳定带来了较大冲击,也给广义负荷建模带来了巨大挑战。随着风电容量增加,广义负荷节点时而呈现电源特性,时而呈现负荷特性,不同特性对应不同模型,其对电力系统仿真计算会产生质的改变,因此分析考虑风电不确定性的广义负荷建模对于电力系统分析十分重要。
基于量测的负荷建模中,时变性是阻碍其走向应用的最大障碍,然而由于风电接入所带来的不确定性更增加了原有负荷建模问题的难度。研究表明,分类与综合是解决负荷时变性问题的有效途径。为此,负荷建模工作欲从研究阶段走向实际应用,不可避免的面临负荷特性分类与综合。对于风电与负荷组成的广义负荷特性,由于风电随机性与负荷时变性的交互影响而呈现与各自不同的复杂特性且该特性有一定地域特征,因此仅根据单日峰谷时段简单分类过于粗略,如能利用固有的自然规律和人类社会周期性特征,则分类结果会更为合理、有效。基于此,有必要更有可能寻求一种客观实用的聚类方法便于精确建模和现场应用。
以往聚类方法在传统负荷建模场景下较好的解决了时变性问题,但随着风电渗透率增加,广义负荷特性的不确定性加剧,因此现有聚类方法难以满足该场景下的分类需求。其中,文献[2]首次提出了分类综合问题并利用同类全部数据进行模型综合;文献[3]分别按照时段和季节分类并用综合的方法对分类结果进行了验证,但该方法分类较主观且未体现日差异性;文献[4]采用KOHONEN神经网络法以标准电压激励下的负荷模型和负荷有功运行水平为特征向量进行分类;文献[5]将多元统计分析中的系统聚类法引入负荷建模领域分析时变规律;文献[6]基于随机过程相关性理论,利用实测样本间相关系数进行直接分类;文献[7]以变电站负荷构成成分为特征向量,分别采用模糊等价关系和模糊C均值算法进行分类;文献[8]以各负荷节点-变电站的不同类型负荷比例为特征向量,基于模糊等价关系的传递闭包法对实测数据进行了模糊分类;文献[9]对建模样本输入输出数据进行分析,建立山峰密度函数自适应确定聚类数和聚类中心;文献[10]以实测响应空间分类方法为基础,提出动特性直接综合法。上述文献通过合理的分类综合方法能较好的解决时变性问题,但部分文献采用的聚类方法由于需人为设定聚类数、聚类中心等,无法排除主观因素,新形势下不具有普遍适用性;另外由于分析对象样本较小,聚类策略相对简单,仅需确定合适的聚类方法和特征向量进行聚类即可划分成界限明显的类别,而面对全年(或更长时间)的风电与负荷大样本数据时,简单聚类策略无法将全部样本数据进行合理分类;除此,限于传统应用场景,上述文献均未考虑时段连续性对模型的影响,由于人类活动与自然规律都是长期周期性的渐变过程,因此考虑训练样本在时间上的连续性,会使得模型更加精确、完备。由此文献[11]考虑了连续时间信息,但该方法以日为最小分析间隔,未考虑时段相似性与日间差异性规律,对于广义负荷建模而言,该尺度较大,无法给出确定时段的合理模型。因此有必要寻求更合理的聚类策略以满足该场景下的分类需求与现场应用。
[1]张旭,梁军,贠志皓,等.考虑风电接入不确定性的广义负荷建模[J].电力系统自动化.数字优先出版
ZHANG Xu,LIANG Jun,YUN Zhihao,et al.Generalized Load Modeling andApplication Considering Uncertainty of Wind Power Integration[J].Automationof Electric Power Systems.
[2]章健.电力系统负荷建模方法的研究[D].北京:华北电力大学,1997.
ZHANG Jian.Studies on modeling methodology of electric loads[D].Beijing:North China Electric Power University,1997.
[3]张伶俐,周文.面向综合的电力负荷动特性建模[J].中国电机工程学报,1999,19(9):36-40.
ZHANG Lingli,ZHOU Wen.The Synthesis of Dynamic Load Characteristics[J].Proceedings of the CSEE,1999,19(9):36-40.
[4]张红斌,贺仁睦,刘应梅.基于KOHONEN神经网络的电力系统负荷动特性聚类与综合[J].中国电机工程学报,2003,23(5):1-5.
ZHANG Hongbin,HE Renmu,LIU Yingmei.The Characteristics Clustering andSynthesis of Electric Dynamic Loads Based on KOHONEN Neural Network[J].Proceedings of the CSEE,2003,23(5):1-5.
[5]石景海,贺仁睦.基于量测的负荷建模——分类算法[J].中国电机工程学报,2004,24(2):78-82.
SHI Jinghai,HE Renmu.Measurement based load modeling—sortingalgorithm[J].Proceedings of the CSEE,2004,24(78-82):116-126.
[6]李欣然,林舜江,刘杨华,等.基于实测响应空间的负荷动特性分类原理与方法[J].中国电机工程学报,2006,26(8):39-44.
LI Xinran,LIN Shunjiang,LIU Yanghua,et al.A New Classification Methodfor Aggregate Load Dynamic Characteristics Based on Field Measured Response[J].Proceedings of the CSEE,2006,26(8):39-44.
[7]李培强,李欣然,陈辉华,等.基于模糊聚类的电力负荷特性的分类与综合[J].中国电机工程学报,2006,25(24):73-78.
LI Peiqiang,LI Xinran,CHEN Huihua,et al.The CharacteristicsClassification and Synthesis of Power Load Based on Fuzzy Clustering[J].Proceedings of the CSEE,2006,25(24):73-78.
[8]黄梅,贺仁睦,杨少兵.模糊聚类在负荷实测建模中的应用[J].电网技术,2006,30(14):49-52.
HUANG Mei,HE Renmu,YANG Shaobing.Application of Fuzzy Clustering inMeasurement-Based Load Modeling[J].Power System Technology,2006,30(14):49-52.
[9]李培强,李欣然,陈辉华,等.基于减法聚类的模糊神经网络负荷建模[J].电工技术学报,2006,21(9):2-6.
LI Peiqiang,LI Xinran,CHEN Huihua,et al.Fuzzy Neural Network LoadModeling Based on Subtractive Clustering[J].Transactions of ChinaElectrotechnical Society,2006,21(9):2-6.
[10]林舜江,李欣然,李培强,等.基于实测响应空间的负荷动特性直接综合方法[J].中国电机工程学报,2007,26(21):36-42.
LIN Shunjiang,LI Xinran,LI Peiqiang,et al.A Novel Direct Method forAggregate Load Dynamic Characteristics Based on Field Measured Response[J].Proceedings of the CSEE,2010(7):80-83.
[11]蒋小亮,蒋传文,彭明鸿,等.基于时间连续性及季节周期性的风速短期组合预测方法[J].电力系统自动化,2010(15):75-79.
JIANG Xiaoliang,JIANG Chuanwan,PENG Minghong,et al.A Short-termCombination Wind Speed Forecasting Method Considering Seasonal Periodicityand Time-continuity[J].Automation of Electric Power Systems,2010(15):75-79.
[12]Q/GDW 392-2009风电场接入电网技术规定.北京:国家电网公司.2009.
发明内容
为解决现有技术存在的不足,本发明公开了广义负荷建模中基于时段性的横向时间轴聚类方法,利用AP算法和考虑时段特性的横向时间轴聚类策略,可将大样本实测数据进行合理划分,利用欧氏距离进行特征向量交叉匹配以判断待分析样本所属类别,通过广义负荷建模方法对各类中全部样本数据建模,以确定待分析样本具体模型并检验聚类策略有效性。
为实现上述目的,本发明的具体方案如下:
广义负荷建模中基于时段性的横向时间轴聚类方法,包括以下步骤:
步骤一:获取全年风电与负荷组成的根母线数据;
步骤二:数据处理,对于风场输出为负的情况,即风机吸收功率,剔除该类数据;对由于数据漏检造成的数据缺失,根据临近数据点进行补充;
步骤三:将步骤二中处理后的全部数据首尾相连,横列组成横向连续数据,按横向时间单元THN分为M段,对全部数据进行横向聚类;
步骤四:基于横向聚类结果,待分析样本数据源被划分为q个横向类,每类由各自聚类中心代表;
步骤五:特征向量交叉匹配,当需建立某时段负荷模型时,通过该时段样本特征向量与全部类别聚类中心进行匹配,该样本属于相似度高的历史样本类时,待研究样本直接采用该类所属的类别的负荷模型。
所述步骤一中根母线数据为风场实测有功运行数据与变电站110kV侧出线负荷功率数据。
所述步骤二根据临近数据点进行补充时,采用三次样条插值法。
所述步骤三中对全部数据进行横向聚类采用AP聚类。
所述步骤三中由于横向时间单元是由多个最小时间间隔T组成,利用THN内相邻最小时间间隔T的变化率序列表征横向数据波动趋势,辅以实际有功功率统计量,构成特征向量如式(14):
式中,p为横向时间单元序列号;b为THN中最小时间间隔数;Wmaxp和Wminp分别为第p个横向时间单元内的最大有功功率和最小有功功率;和Ws2p分别为第p个横向时间单元内平均有功功率和方差,γp1,γp2,…,γpb分别为THN内b个最小时间间隔波动率序列。
最小时间间隔波动率序列的通式γpi,其计算式如式(15):
式中,i为最小时间间隔序列号;p为横向时间单元序列号;γpi为第p个横向时间单元内第i个最小时间间隔波动率;Pw为负荷与风电组成的根母线有功功率序列,j为风功率序列标号;THN为横向时间单元;tmin为最小备选时间;X1和X2分别为最小时间间隔的始样本点和末样本点;T为最小时间间隔。
所述步骤五中,将各功率段数据以时间为基准,匹配原始一一对应的有功功率与电压值,利用Levenberg-Marquardt神经网络法学习并提取各段节点特征;最终合并各段模型,形成如式(17)的统一模型结构:
式中,psm表征电源特性,plm表征负荷特性;psm(vsm)、plm(vlm)为各段下提取的电源特性和负荷特性关系表达式;vsm、vlm为各段中根母线电压;以电源特性为例,psm(psm∈ps·[0.1m-1.1,0.1m-1)|psm<0)表示在psm<0的条件下,根母线有功出力落在ps·[0.1m-1.1,0.1m-1)这段的概率,在此概率约束下其出力随电压变化的特征关系为psm(vsm),负荷特性类似。m为分段标识;Ps为基准功率,用于数据归一化处理,该值应大于全部功率数据绝对值最大值,在此基础上根据实际数据合理选取;εsm、εlm分别为电源特性与负荷特性分段范围裕度,以保证分段限值为整数;pmin、pmax分别为功率最小值和功率最大值。
所述最小时间间隔应取满足原则情况下的最大值,见式(8)(9)。
式中,t为备选时间间隔,tmin和tmax分别为最小备选时间和最大备选时间,实际分别取采样间隔和全部研究时间,以分钟计;int{·}为取整函数;γti为备选时间间隔t下,采样间隔序列i的功率波动率,Pw为负荷与风电组成的根母线有功功率序列,{·}内为风功率序列标号。σ为变化率阈值。
本发明的有益效果:
本申请提出基于时段相似性的横向时间轴聚类策略。首先,就聚类方法而言,本申请引入的AP算法因特殊的消息传递机制与竞争机制,不必事先确定聚类数和聚类中心,聚类质量高且对大数据聚类效果明显,恰好满足本申请对客观聚类、历史大样本数据的分析需求。就聚类策略而言,本申请以实测响应空间本身的时段特性为依据提出横向时间轴聚类策略,基于时段连续性将研究样本划分为聚类时间单元,由于人类习惯等存在日内时段规律,通过对连续时间序列组成的时间单元进行聚类分析(本申请中称为横向时间轴聚类),可自动划分为符合日内时段特征规律的横向类。由此,利用AP算法和横向时间轴聚类,足以实现考虑风电随机性与负荷时变性的客观分类;利用欧氏距离进行特征向量交叉匹配以判断待分析样本所属类别;利用文献[1]广义负荷建模方法对各类中全部样本数据建模,以确定待分析样本具体模型并检验聚类策略有效性。仿真表明,该聚类方法合理、有效,可在大数据样本基础上构建更为贴近实用的精确模型,为风电接入后的仿真分析和调度控制提供模型基础和辅助参考。
1.本申请将高质量的、排除人为主观因素的、适用于大样本数据的AP算法引入电力系统广义特性聚类,通过实测全年样本空间的直观聚类结果与特性综合说明了该算法的有效性。
2.提出了横向时间轴聚类策略。通过横向时间轴聚类可见,由于存在较强的不确定性,即使对于连续两天,时段特性也不尽相同,因此主观分为高峰低谷时段的划分方法不适用于风电接入的不确定性场景;利用本申请提出的考虑日内波动规律的横向时间轴聚类方法自动划分横向数据,客观合理,更利于该场景下精确建模。
3.利用AP算法和考虑时段特性的横向时间轴聚类策略,可将大样本实测数据进行合理划分,通过特征向量交叉匹配判断样本所属类别,利用广义负荷建模建立精确模型并检验聚类策略有效性。仿真结果表明,经过聚类分析后的广义负荷建模,可在满足精确性的基础上便于模型走向实用化,有利于提高电力系统仿真的精确性与实效性,因此本申请方法对于考虑风电接入不确定性广义负荷建模有一定作用,可为仿真计算与调度控制提供辅助参考。
附图说明
图1(a)AP算法消息传递机制传递responsibility;
图1(b)AP算法消息传递机制传递availability;
图2根母线节点组成示意图;
图3横向时间轴聚类结果示意图;
图4类别3自描述拟合图;
图5(a)类别1待检验时段拟合效果;
图5(b)类别3待检验时段拟合效果。
具体实施方式:
下面结合附图对本发明进行详细说明:
AP聚类算法
传统常用的贝叶斯分类法和距离分类法模型简单,难以反映时段性与负荷运行水平间错综复杂的关系,另外其面向大样本数据的聚类困难。近年来,负荷特性分类与综合中常用K-means聚类算法、模糊聚类算法、神经网络法等,但大部分需要人为设定如聚类数、聚类中心等,无法排除主观因素,由此本申请引入AP聚类算法。Brendan J.Frey等提出的仿射传播聚类(Affinity Propagation Clustering,AP聚类)是一种有效的聚类方法。该方法为无监督聚类方法,不必事先确定聚类数和聚类中心,聚类质量高,其对大数据聚类效果尤为明显,恰好满足本申请的大数据样本点、客观聚类需求。
与传统K-means聚类不同,AP算法将每个初始样本点均看作候选聚类中心,从而消除了由于初始聚类中心随机选择不当对聚类结果造成的影响。聚类过程中通过样本点间的消息机制传递信息从而确定该样本点的聚类中心或是否其本身为聚类中心,其中消息传递机制包括两种:responsibility和availability,如图1(a)-1(b)。
其中,r(i,k)反映k点作为i点聚类中心的合适程度,该值越大表明候选聚类中心k越可能成为真正聚类中心;a(i,k)表示从候选聚类中心k发送到样本点i的消息,反映i点选择k点为其聚类中心的合适程度,该值越大表明样本点越可能成为以k作为聚类中心的类别中。基于AP算法的具体计算步骤如下:
Step1:变量初始化。确定样本点数N、迭代次数M、欲聚类数据的n个特征向量(如时段、季节因素和负荷运行水平等),按式(1)计算N个样本点Xi=(xi1,xi2,…xin),(i=1,2,…N)的相似度矩阵S作为输入量,其中对参考度P(即对角线元素S(i,i),该值越大则聚类数越多)赋以相同的值(一般为S均值),初始化r和a阵为0,设置适当的阻尼系数λ(控制迭代速度,一般取0.5~1)。
S(i,k)=-||(xi1-xk1)2+(xi2-xk2)2+…(xin-xkn)2|| (1)
Step2:计算各样本点间r(i,k)和a(i,k),如式(2)-(5),其中r(k,k)和a(k,k)分别为自responsibility与自availability,二者值越大,则k成为聚类中心可能性越大。
Step3:根据设定的阻尼系数λ循环迭代r和a,如式(6)-(7)。
ri=(1-λ)·ri+λ·ri-1 (6)
ai=(1-λ)·ai+λ·ai-1 (7)
Step4:对于k点,直至寻得r(i,k)+a(i,k)=max{r(i,j)+a(i,j)}(j=1,2,…,N),则k为样本点i的聚类中心。当全部聚类结果较上次变化不大或达到迭代循环上限M,则聚类结束;否则返回Step2。
式中,i,j,k均为样本点序号,此处为便于说明公式,分别以i,j,k表示;r(i,k)为样本点Xi和Xk间的responsibility;S(i,k)为样本点Xi和Xk相似度,计算方法见式(1);a(i,j)表示样本点Xi和Xj间的availability;S(i,j)为样本点Xi和Xj间相似度;r(k,k)为样本点Xk的自responsibility;P(k)为样本点k的参考度;a(k,j)为表示样本点Xk和Xj间的availability;S(k,j)为样本点Xk和Xj相似度;ri表示第i次迭代的responsibility;ri-1表示第i-1次迭代的responsibility;λ为阻尼系数,用以控制迭代速度,一般取0.5-1;ai为第i次迭代的availability;ai-1为第i-1次迭代的availability。
由此该算法可根据固有的时段特征自动挖掘数据的内在规律,与本申请中排除人为因素、客观合理的聚类需求相吻合。
横向时间轴聚类策略
为应对风电与负荷组成根母线功率的不确定性,首先根据时段和负荷运行水平等特征对数据样本进行聚类分析,将广义负荷特性相似的聚为一类,特性相差较大的为异类,不同类别建立不同模型。如此有针对性的选用模型,可在一定程度上克服风电随机性与负荷时变性带来的不确定性,提高模型精确性。本申请横向时间轴聚类,通过横向客观聚类,体现风电和负荷组成的跟母线数据的时段特性。
将广义负荷特性相似的聚为一类,特性相差较大的为异类,该定义是所有负荷建模文献中对特性聚类的公认定义。对相似性的具体定义不同则引申出不同的聚类方法。若必须清晰描述,不妨以特征向量样本点在空间的距离描述,若某两点间空间距离较小则为特性相似,空间距离大则特性相差较大。距离小或大根据一定的阈值判断,该阈值的选取视具体问题选取。
横向时间轴聚类策略相关概念,为便于理解,首先对横向时间轴聚类策略中有关概念作必要说明。
根母线广义负荷特性:本申请研究对象,即图中的根母线节点,相当于从系统侧看入、与该母线相连的风电和邻近基础负荷等的综合等值,其外特性取决于该时刻负荷与风功率的相对大小,它会随负荷时变性和风功率波动性而变化。从系统角度看,这种广义负荷特性相比于传统单纯负荷,有质的变化,对系统的影响也大不相同,因此研究这类节点负荷特性有重要意义。
所谓横向时间轴,指一日(或更长)的时间范围内,连续时间序列组成的坐标轴。分析的时间窗口视具体对象而定。本申请中,以日为横向分析时间尺度,则横向时间轴的时间窗口Td为一天,其由Nu个最小时间间隔T组成,如此分析横向时间轴数据,有助于总结单日不同时段波动规律,从而整体把握日内时段相似性与差异性。
最小时间间隔T,即能保持该时段内根母线功率近似不变的最大时段,为此该值的选取应满足两个原则:一是最小时间间隔内变化率不超过规定的变化率阈值,则认为该时段内功率近似不变,可将其直接作为一个整体单位进行聚类;二是为保证聚类数适中、便于进一步分析,最小时间间隔不能无限小,如此应取满足原则一情况下的最大值,见式(8)(9)。
式中,t为备选时间间隔,tmin和tmax分别为最小备选时间和最大备选时间,实际分别取采样间隔和全部研究时间,以分钟计,如对于采样间隔为5min的全年数据,tmin和tmax分别为5和525600;int{·}为取整函数;γti为备选时间间隔t下,采样间隔序列i的功率波动率,Pw为负荷与风电组成的根母线有功功率序列,{·}内为风功率序列标号。σ为变化率阈值。
时段连续性与聚类时间单元:由于人类活动和自然规律都是长期的渐变过程,如相似时段、周期性等。因此在聚类分析时将时段连续性纳入考虑范畴,可使模型更加完备、健壮,同时避免由于研究对象过多,造成聚类后所建模型不便应用。为此,实际分析前,先将全部研究数据按横向时间单元THN分为若干等长的连续时段。
所谓横向时间单元THN,即由若干最小时间间隔T构成的一段横向连续时间,目的是通过横向时间聚类提取时段特征,横向时间单元THN应介于最小时间间隔T与横向时间序列时间窗口Td之间,显然若THN过长,则分析结果较粗略,不可靠;THN过短,则会导致聚类单元过多而不便于具体应用,因此应兼顾可靠性与应用性合理选取。在此基础上,以聚类时间单元组成的对象进行聚类分析,可增强聚类分析后所建负荷模型的简洁性和实用性。
横向时间轴聚类
将全年数据进行横向时段划分,目的在于通过对时段特性的挖掘,将待研究时段的特征向量与各类聚类中心迅速匹配,找到该时段所属类别并应用该段模型,快速简便。由于横向时间轴聚类着眼于分析较小的数据单元,如此需先对现场数据进行数据处理以获得更可靠合理的分析数据,再构建反映时段特性的特征向量进行聚类,最终通过特征向量交叉匹配找到其所属类别。
数据处理:对于分析较小数据单元间(如分钟级、小时级数据)特性异同的横向时间轴聚类,每个采样点的数据准确性都会对聚类分析结果产生质的影响,因此需对现场采样数据进行处理。对于风场输出为负的情况,即风机吸收功率,一般是由于风机停机造成,而停机对分析风电特性显然意义不大,可剔除该类数据;另外,现场常会遇到由于数据漏检等造成的数据缺失问题,需根据临近数据点进行合理补充。常用插值进行缺失数据值推测,本申请采用三次样条插值,该函数不但可实现平滑接近真实值,还能保证其局部特性全部为正,符合风电特性规律。
特征向量:为提取时段特性,需考虑时段连续性进行横向聚类分析,由此需获取能够反映横向时间单元THN时段特征的特征向量。将全部数据首尾相连,横列组成横向连续数据,按横向时间单元THN分为M段,对全部数据进行横向聚类,如图3,如此打破单日数据单独聚类的限制,将全部数据统一考虑,聚类更合理。由于横向时间单元是由多个最小时间间隔T组成,利用THN内相邻最小时间间隔T的变化率序列表征横向数据波动趋势,辅以实际有功功率统计量,构成特征向量如式(14)。
式中,p为横向时间单元序列号;b为THN中最小时间间隔数;Wmaxp和Wminp分别为第p个横向时间单元内的最大有功功率和最小有功功率;和Ws2p分别为第p个横向时间单元内平均有功功率和方差。γp1,γp2,…,γpb分别为THN内b个最小时间间隔波动率序列,以γpi为例,其计算式如式(15)。
式中,i为最小时间间隔序列号;p为横向时间单元序列号;γpi为第p个横向时间单元内第i个最小时间间隔波动率;Pw为负荷与风电组成的根母线有功功率序列,j为风功率序列标号;THN为横向时间单元;tmin为最小备选时间;X1和X2分别为最小时间间隔的始样本点和末样本点;T为最小时间间隔。如此通过横向时间轴聚类自适应化为q类(如图3),有助于根据单日不同时段波动规律,整体把握日内时段相似性与差异性,从而提取相应时段特征。
特征向量交叉匹配:基于横向聚类结果,待分析样本数据源被划分为q个横向类,如图3,每类可由各自聚类中心代表。当需建立某时段负荷模型时,只需利用该时段波动数据形成式(14)的特征向量,通过式(16)欧氏距离判断该时段样本与历史样本中各类聚类中心特征向量的相似度,判断该样本是属于相似度高(欧式距离小)的历史样本类。
式中,ρ(A,B)为序列A=[a(1),a(2),…,a(n)]与序列B=[b(1),b(2),…,b(n)]的欧式距离。
区别于传统依据特定时段一一对应高峰低谷时段的方法,本申请通过样本特征向量与全部类别聚类中心匹配,打破时段顺序限制,故称之为交叉匹配。
广义负荷建模
利用横向时间轴分析全年数据,可获得考虑时段相似性的分类数据用以精确建模,但如何得到每类数据的广义负荷模型,需借助于负荷特性综合的方法。由于风电接入带来了节点特性的不确定性,包括功率流向改变与节点的不确定性变化,而传统建模方法因不具备随机特征描述能力而难以应用于不确定场景下的分析计算,由此采用文献[1]中引入概率信息的广义负荷学习与建模新方法。
该方法针对风电接入后节点功率流向的改变,将节点特性以有功功率为参考,划分为电源特性与负荷特性;针对节电特性的不确定性变化,基于历史实测有功功率数据对样本空间进行自适应分段并统计其概率分布,各有功功率段下的概率即为其概率标识;将各功率段数据以时间为基准,匹配原始一一对应的有功功率与电压值,利用Levenberg-Marquardt神经网络法学习并提取各段节点特征;最终合并各段模型,形成如式(17)的统一模型结构。
式中,psm表征电源特性,plm表征负荷特性;psm(vsm)、plm(vlm)为各段下提取的电源特性和负荷特性关系表达式;vsm、vlm为各段中根母线电压;以电源特性为例,psm(psm∈ps·[0.1m-1.1,0.1m-1)|psm<0)表示在psm<0的条件下,根母线有功出力落在ps·[0.1m-1.1,0.1m-1)这段的概率,在此概率约束下其出力随电压变化的特征关系为psm(vsm),负荷特性类似。m为分段标识;Ps为基准功率,用于数据归一化处理,该值应大于全部功率数据绝对值最大值,在此基础上根据实际数据合理选取;εsm、εlm分别为电源特性与负荷特性分段范围裕度,以保证分段限值为整数;pmin、pmax分别为功率最小值和功率最大值。
该方法通过引入概率信息,可从节点特性角度进行考虑风电接入后随机特征的各场景下潮流计算、稳定性计算以及系统风险评估等,是对传统建模方法在应用场景上的扩展和延伸。如此通过广义负荷特性的综合,不仅可以得到精确的广义负荷模型,也可以检验负荷特性聚类的正确性。
实例仿真
以山东某风场实测有功运行数据与该地某变电站110kV侧出线负荷功率数据验证本申请方法正确性,以2011年数据为训练样本,以2012年数据为测试样本,采样间隔均为5min。为获得完整的数据样本进行聚类与建模,首先应通过合理计算以获得风场无功功率样本和根母线电压样本,数据获取方法与文献[1]类似。根据实际比例,以New England-39节点标准系统为例,检验本申请所提方法的正确性。
利用风场获得的有功功率样本定功率因数得到风场无功功率,将风场有功功率与负荷有功功率叠加,数据获取方法具体为:将风场无功功率与负荷无功功率叠加,将风场与负荷组成的根母线节点看作PQ节点,如图2所示,根据实际比例,将根母线节点功率数据加入New England-39节点标准系统,进行潮流计算,最终获得根母线节点所有有功功率、无功功率和电压样本。
最小时间间隔T选取:本申请聚类方法的可靠性与应用性取决于最小时间间隔,因此其选取的合理与否尤为重要。该参数据式(8)和(9)合理选取,分别以5、10、15、20、25为备选时间间隔,相应时间内波动率最大值计算如表1所示。
表1备选时间间隔内最大波动率
对于备选时间间隔5min,由于采样间隔tmin为5min,据式(9)得其波动率为0。根据文献[12],本申请中波动率阈值σ为0.3,故据式(8)合理选择最小时间间隔为15min。
横向时间轴聚类:考虑聚类结果可靠性与模型的实用性因素,横向时间单元THN折中取为4h,则单日被分为6个连续时段。以2011年数据研究样本,构造特征向量进行横向时间轴聚类,限于篇幅仅呈现15天聚类结果如表2。
表2横向时间轴聚类结果
依据本申请方法,2011年根母线前三月数据总共自动划分为4类,聚类中心分别为225II、265VI、283II和215V,其中阿拉伯数字对应所属类别聚类中心的单日序号,罗马数字对应聚类中心的单日时段序号,其所在位置对应实际单日和时段。据表中数据分析可得,即使相邻两天,时段划分也不完全一致,如第212、213天,因此以往文献中直接按高峰、低谷主观划分时段的方法,在应对本申请中由于风电接入带来不确定性的场景时,会因无法排除人为因素而导致最终分析结果难以可靠,有必要利用本申请方法打破传统时段观念,考虑日间差异性与相似性逐一分析单日特性,为进一步分析计算提供更客观、合理的数据支持。
广义负荷建模
对利用横向时间轴划分的数据类进行广义负荷建模,不但可获取相应模型,还可通过综合的办法进行聚类结果的验证。分别对4类数据建立独立的广义负荷模型,为保证样本数据充足,需利用各类中全部实测响应样本数据,限于篇幅本申请仅呈现类别3拟合效果以证明该方法的描述能力,如表3和图4。
表3类别3概率分布与拟合误差表
由于打破了传统时间概念,数据集中于某一区域波动而不随时间平缓过渡,即使这种情况下,本申请选用的建模方法追踪拟合效果仍较好,仅在个别端点处略有偏差,如此可为本申请提出的聚类方法提供可靠的检验手段,也为时段仿真分析提供模型基础。
特征向量交叉匹配与横向时间轴聚类策略验证
为验证聚类结果正确性,取2012年某时段进行验证,通过构建特征向量,分别与2011年所建立横向聚类模型的聚类中心特征向量进行欧式距离计算,结果如表4。
表4待验证样本与4类聚类中心欧式距离
该时段与类别1聚类中心的欧式距离最小为3.5582,与其余各类欧式距离相差较大,因此属于类别1。将待验证时段数据进行分别利用类别1和类别3模型进行拟合,如图5(a)-5(b)。
显然,该模型利用类别1拟合效果较好,而类别3拟合误差较大,通过特性综合的方法再次证明了本申请提出聚类方法的正确性。如此利用同类数据建立广义负荷模型,对数据分析更细致,结果也更加精确。本申请方法能够通过对考虑时段性的聚类分析,自动区分不同时段,可在大数据样本基础上构建更为贴近实用的精确模型。
Claims (7)
1.广义负荷建模中基于时段性的横向时间轴聚类方法,所述横向时间轴,指一日或更长的时间范围内,连续时间序列组成的坐标轴,其特征是,包括以下步骤:
步骤一:获取全年风电与负荷组成的根母线数据;
步骤二:数据处理,对于风场输出为负的情况,即风机吸收功率,剔除该类数据;对由于数据漏检造成的数据缺失,根据临近数据点进行补充;
步骤三:将步骤二中处理后的全部数据首尾相连,横列组成横向连续数据,按横向时间单元THN分为M段,对全部数据进行横向聚类;
步骤四:基于横向聚类结果,待分析样本数据源被划分为q个横向类,每类由各自聚类中心代表;
步骤五:特征向量交叉匹配,当需建立某时段负荷模型时,通过该时段样本特征向量与全部类别聚类中心进行匹配,该样本属于相似度高的历史样本类时,待研究样本直接采用该类所属的类别的负荷模型;
横向聚类即为“横向时间轴聚类”,具体为:将全年数据进行横向时段划分,将待研究时段的特征向量与各类聚类中心迅速匹配,找到该时段所属类别并应用该段模型,由于横向时间轴聚类着眼于分析较小的数据单元,如此需先对现场数据进行数据处理以获得更可靠合理的分析数据,再构建反映时段特性的特征向量进行聚类,最终通过特征向量交叉匹配找到其所属类别;
所述步骤三中由于横向时间单元是由多个最小时间间隔T组成,利用THN内相邻最小时间间隔T的变化率序列表征横向数据波动趋势,辅以实际有功功率统计量,构成特征向量如式(14):
式中,p为横向时间单元序列号;b为THN中最小时间间隔数;Wmaxp和Wminp分别为第p个横向时间单元内的最大有功功率和最小有功功率;和Ws2p分别为第p个横向时间单元内平均有功功率和方差,γp1,γp2,…,γpb分别为THN内b个最小时间间隔波动率序列。
2.如权利要求1所述的广义负荷建模中基于时段性的横向时间轴聚类方法,其特征是,所述步骤一中根母线数据为风场实测有功运行数据与变电站110kV侧出线负荷功率数据。
3.如权利要求1所述的广义负荷建模中基于时段性的横向时间轴聚类方法,其特征是,所述步骤二根据临近数据点进行补充时,采用三次样条插值法。
4.如权利要求1所述的广义负荷建模中基于时段性的横向时间轴聚类方法,其特征是,最小时间间隔波动率序列的通式γpi,其计算式如式(15):
式中,i为最小时间间隔序列号;p为横向时间单元序列号;γpi为第p个横向时间单元内第i个最小时间间隔波动率;Pw为负荷与风电组成的根母线有功功率序列,j为风功率序列标号;THN为横向时间单元;tmin为最小备选时间;X1和X2分别为最小时间间隔的始样本点和末样本点;T为最小时间间隔。
5.如权利要求1所述的广义负荷建模中基于时段性的横向时间轴聚类方法,其特征是,所述步骤五中,将各功率段数据以时间为基准,匹配原始一一对应的有功功率与电压值,利用Levenberg-Marquardt神经网络法学习并提取各段节点特征;最终合并各段模型,形成如式(17)的统一模型结构:
式中,psm表征电源特性,plm表征负荷特性;psm(vsm)、plm(vlm)为各段下提取的电源特性和负荷特性关系表达式;vsm、vlm为各段中根母线电压;以电源特性为例,psm(psm∈ps·[0.1m-1.1,0.1m-1)|psm<0)表示在psm<0的条件下,根母线有功出力落在ps·[0.1m-1.1,0.1m-1)这段的概率,在此概率约束下其出力随电压变化的特征关系为psm(vsm),负荷特性类似;m为分段标识;Ps为基准功率,用于数据归一化处理,该值应大于全部功率数据绝对值最大值,在此基础上根据实际数据合理选取;εsm、εlm分别为电源特性与负荷特性分段范围裕度,以保证分段限值为整数;pmin、pmax分别为功率最小值和功率最大值。
6.如权利要求1所述的广义负荷建模中基于时段性的横向时间轴聚类方法,其特征是,所述最小时间间隔应取满足原则情况下的最大值,见式(8)(9):
式中,t为备选时间间隔,tmin和tmax分别为最小备选时间和最大备选时间,实际分别取采样间隔和全部研究时间,以分钟计;int{·}为取整函数;γti为备选时间间隔t下,采样间隔序列i的功率波动率,Pw为负荷与风电组成的根母线有功功率序列,{·}内为风功率序列标号,σ为变化率阈值。
7.如权利要求1所述的广义负荷建模中基于时段性的横向时间轴聚类方法,其特征是,所述步骤三中对全部数据进行横向聚类采用AP聚类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410453541.XA CN104200032B (zh) | 2014-09-05 | 2014-09-05 | 广义负荷建模中基于时段性的横向时间轴聚类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410453541.XA CN104200032B (zh) | 2014-09-05 | 2014-09-05 | 广义负荷建模中基于时段性的横向时间轴聚类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104200032A CN104200032A (zh) | 2014-12-10 |
CN104200032B true CN104200032B (zh) | 2017-04-05 |
Family
ID=52085325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410453541.XA Expired - Fee Related CN104200032B (zh) | 2014-09-05 | 2014-09-05 | 广义负荷建模中基于时段性的横向时间轴聚类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104200032B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005708B (zh) * | 2015-08-13 | 2017-10-24 | 山东大学 | 一种基于ap聚类算法的广义负荷特性聚类方法 |
CN105069236B (zh) * | 2015-08-13 | 2018-11-13 | 山东大学 | 考虑风电场节点空间相关性的广义负荷联合概率建模方法 |
CN105119278A (zh) * | 2015-08-26 | 2015-12-02 | 河海大学 | 一种特殊负荷建模方法 |
CN105138849B (zh) * | 2015-09-07 | 2018-04-10 | 山东大学 | 一种基于ap聚类的无功电压控制分区方法 |
CN106897728B (zh) * | 2015-12-21 | 2019-12-17 | 腾讯科技(深圳)有限公司 | 基于业务监控系统的样本选择方法、装置和系统 |
TWI628935B (zh) * | 2016-01-29 | 2018-07-01 | 中華電信股份有限公司 | Request traffic grouping method |
CN105721086B (zh) * | 2016-03-11 | 2018-05-01 | 重庆科技学院 | 基于ukfnn的无线信道场景识别方法 |
CN107219809B (zh) * | 2016-11-24 | 2019-07-23 | 浙江浙能中煤舟山煤电有限责任公司 | 一次风系统中一次风机失速的预防方法 |
CN108804563B (zh) * | 2018-05-22 | 2021-11-19 | 创新先进技术有限公司 | 一种数据标注方法、装置以及设备 |
CN111523230B (zh) * | 2020-04-22 | 2023-05-26 | 国网能源研究院有限公司 | 一种风光负荷复合典型场景的自适应聚类方法 |
CN113627674A (zh) * | 2021-08-12 | 2021-11-09 | 中国华能集团清洁能源技术研究院有限公司 | 一种分布式光伏电站出力预测方法、装置及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298707A (zh) * | 2011-08-24 | 2011-12-28 | 辽宁力迅风电控制系统有限公司 | 基于连续时间段聚类与支持向量机建模的风电功率预测方法 |
CN103177188A (zh) * | 2013-04-02 | 2013-06-26 | 东南大学 | 一种基于特征映射的电力系统负荷动特性分类方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7565335B2 (en) * | 2006-03-15 | 2009-07-21 | Microsoft Corporation | Transform for outlier detection in extract, transfer, load environment |
-
2014
- 2014-09-05 CN CN201410453541.XA patent/CN104200032B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298707A (zh) * | 2011-08-24 | 2011-12-28 | 辽宁力迅风电控制系统有限公司 | 基于连续时间段聚类与支持向量机建模的风电功率预测方法 |
CN103177188A (zh) * | 2013-04-02 | 2013-06-26 | 东南大学 | 一种基于特征映射的电力系统负荷动特性分类方法 |
Non-Patent Citations (2)
Title |
---|
聚类分析在负荷模型分类研究中的应用;郑晓雨;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20090215;第C042-238/10-11、23、31页 * |
面向综合的电力负荷动特性建模;张伶俐 等;《中国电机工程学报》;19990930;第19卷(第9期);第38页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104200032A (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104200032B (zh) | 广义负荷建模中基于时段性的横向时间轴聚类方法 | |
CN111612651B (zh) | 一种基于长短期记忆网络的异常电量数据检测方法 | |
CN104319807B (zh) | 一种基于Copula函数获取多风电场容量可信度的方法 | |
CN108022001A (zh) | 基于pca和分位数回归森林的短期负荷概率密度预测方法 | |
CN105005708B (zh) | 一种基于ap聚类算法的广义负荷特性聚类方法 | |
Fang et al. | Application of gray relational analysis to k-means clustering for dynamic equivalent modeling of wind farm | |
CN111753893A (zh) | 一种基于聚类和深度学习的风电机组功率集群预测方法 | |
CN107679687A (zh) | 一种光伏出力建模方法以及发电系统可靠性评估方法 | |
CN106529719A (zh) | 基于粒子群优化算法风速融合的风电功率预测的方法 | |
CN112232561B (zh) | 基于约束并行lstm分位数回归的电力负荷概率预测方法 | |
CN104408562A (zh) | 一种基于bp神经网络的光伏系统发电效率综合评估方法 | |
CN107394809A (zh) | 基于风速时间周期特征的风电并网系统风险评估方法 | |
CN107145707A (zh) | 一种计及光伏出力不确定性和全寿命周期成本的配电网变压器规划方法 | |
CN105005878B (zh) | 一种坚强智能电网的综合评估方法 | |
CN107292502A (zh) | 一种配电网可靠性评估方法 | |
CN109165819A (zh) | 一种基于改进AdaBoost.M1-SVM的有源配电网可靠性快速评估方法 | |
CN106611243A (zh) | 一种基于garch模型的风速预测残差修正方法 | |
CN114357670A (zh) | 一种基于bls和自编码器的配电网用电数据异常预警方法 | |
CN110570016A (zh) | 一种计及多阶段性能的抗灾型骨干网架优化方法 | |
CN106768933A (zh) | 一种基于遗传算法的风电机组叶片故障诊断方法 | |
CN103632314B (zh) | 基于概率统计的风电场广义节点特征建模方法 | |
CN109711707A (zh) | 一种船舶动力装置综合状态评估方法 | |
CN115660182A (zh) | 基于最大期望样本加权神经网络模型的光伏出力预测方法 | |
Jin et al. | Wind and photovoltaic power time series data aggregation method based on an ensemble clustering and Markov chain | |
CN104200106B (zh) | 广义负荷建模中基于季节性的纵向时间轴聚类方法 |
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: 20170405 |
|
CF01 | Termination of patent right due to non-payment of annual fee |