CN105139094A - 基于离散粒子群和最小二乘的电网覆冰定量长期预报方法 - Google Patents

基于离散粒子群和最小二乘的电网覆冰定量长期预报方法 Download PDF

Info

Publication number
CN105139094A
CN105139094A CN201510579268.XA CN201510579268A CN105139094A CN 105139094 A CN105139094 A CN 105139094A CN 201510579268 A CN201510579268 A CN 201510579268A CN 105139094 A CN105139094 A CN 105139094A
Authority
CN
China
Prior art keywords
icing
formula
population
particle
forecasting
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
CN201510579268.XA
Other languages
English (en)
Other versions
CN105139094B (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.)
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd, Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510579268.XA priority Critical patent/CN105139094B/zh
Publication of CN105139094A publication Critical patent/CN105139094A/zh
Application granted granted Critical
Publication of CN105139094B publication Critical patent/CN105139094B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明介绍了一种输电线路覆冰长期定量预报方法,该方法为下述步骤:1.收集历史环流指数数据与覆冰资料;2.计算环流指数与平均覆冰日数的相关系数,初步选择覆冰预报因子;3.利用离散粒子群算法,得出可预报性最强的覆冰预报因子;4.以选出的预报因子为基础利用历史数据建立最小二乘回归预测模型;5.利用所建模型对冬季电网覆冰发生程度进行定量预报。本方法的有益效果是:1、可提前一个月对未来一个季度(冬季)的电网覆冰进行定量预测;2、本发明可操作性强;3、预报输电线路电网覆冰的准确率高;4、解决了电网覆冰长期定量预测的难题。根据预测结论,可科学合理地规划电网应对覆冰所需的人力、物资,并进行相关的经济部署,实现电网覆冰的提前应对,减少电网覆冰所造成的损失。

Description

基于离散粒子群和最小二乘的电网覆冰定量长期预报方法
技术领域
本发明属于输配电技术领域,尤其涉及一种电网覆冰长期预报方法。
背景技术
电网覆冰是严重危害输电线路安全稳定运行的自然灾害之一,自1954年第1次输电线路冰害事故以来,大面积冰害事故在全国时有发生,给国民经济带来了巨大的损失。特别是2008年初,我国南方地区遭遇了一场罕见的大面积、长时间严重的雨雪冰冻灾害,由于抗冰资源准备不足,造成电网企业直接经济损失达100多亿元,工厂、医院和居民区停电,京广电气化铁路停运,对社会稳定和人民生产与生活构成了严重威胁。而准确的电网覆冰长期预测是电力企业在人力、物资科学规划、经济部署等的前提条件,具有重大意义。
目前电网覆冰长期预测方法只是针对电网覆冰进行了程度预报,而同一程度的覆冰情况也会有重有轻,难以精细化指导抗冰决策,另外程度预报方法需要包含人为经验因素,难以利用计算机进行程序化处理,工作量大,因此,亟需一种能程序化处理的电网覆冰定量长期预报方法。
发明内容
本发明要解决的技术问题是:针对目前电网覆冰定量长期预测现有技术不足,提供一种基于离散粒子群和最小二乘回归的电网覆冰定量长期预报方法,可提前一个月对未来一个季度(冬季)的电网覆冰进行定量预测。
本发明的解决方案是:所提供的这种基于离散粒子群和最小二乘回归的电网覆冰定量长期预报方法,系考虑74项环流指数对冬季电网覆冰的影响,利用离散粒子群算法选择可预报性最强的冬季覆冰预报因子,以选出的预报因子为基础利用历史数据建立最小二乘回归模型,得出冬季电网平均覆冰日数的预报结论。该方法,包括如下步骤:
计算每个月的74项环流指数与预报区域的历年平均覆冰日数的相关系数,并求取每项环流指数与平均覆冰日数的最大相关系数;
取所得最大相关系数绝对值最大的前m个因子,作为覆冰预报的初步因子;
利用离散粒子群算法,得到粒子群的最优位置,在初步因子中选出可预报性最强的覆冰预报因子;
以选出的可预报性最强的覆冰预报因子为基础,建立预测模型:
yt=b0+b1Xt1+b2Xt2+…+bNXtN+ut
式中yt为平均覆冰日数;b0,b1,…,bN为模型参数;Xt1,…XtN为所选的预报因子;ut为随机误差;N为所选预报因子的个数;
利用最小二乘方法求取模型参数,并将当前预报因子的值带入预测模型求解得到未来冬季电网的平均覆冰日数。
所述m取值为8~25。离散粒子群算法设定粒子总数为10~30。
本发明的有益效果是:
1、使用该方法可以在11月份预报未来12月至次年2月的冬季电网整体覆冰程度;
2、本发明流程清晰、可操作性强;
3、预报准确率高(准确度大于90%)、实用性强;
4、解决了电网覆冰长期定量预测的难题。根据预测结论,可科学合理地规划电网应对覆冰所需的人力、物资,并进行相关的经济部署,实现电网覆冰的提前应对,减少电网覆冰所造成的损失。
附图说明
图1为本发明方法流程图。
具体实施方式
如图1所示,本发明的基于离散粒子群和最小二乘回归的电网覆冰定量长期预报方法具体包括下述的步骤:
1.收集历史环流指数数据与覆冰资料。
收集自1951年以来每个月的包括亚洲极涡、西太平洋副热带高压等74项环流指数数据;通过气象部门收集历年平均覆冰日数。
2.计算环流指数与平均覆冰日数的相关系数,初步选择覆冰预报因子。
利用相关系数计算公式计算3-10月份的74项环流指数与平均覆冰日数的相关系数,相关系数计算的公式如式(1)所示:
R X j = n Σ i = 1 n x i j d i - Σ i = 1 n x i j · Σ i = 1 n d i n Σ i = 1 n x i j 2 - ( Σ i = 1 n x i j ) 2 · n Σ i = 1 n d i 2 - ( Σ i - 1 n d i ) 2 - - - ( 1 )
式中,为第j月74项环流指数X与平均覆冰日数的相关系数;xij为74项环流指数的值,i为年份序号,j为月份序号,n为历史年份数量;di为第i年的平均覆冰日数。
利用式(2)计算每项环流指数X与平均覆冰日数的最大相关系数,
R X m a x = m a x ( R X 3 , R X 4 , ... , R X 10 ) - - - ( 2 )
式中RXmax为环流指数X与平均覆冰日数的最大相关系数,max为所对应的月份。
取所得RXmax绝对值最大的前m个因子,作为覆冰预报的初步因子。
3.利用离散粒子群算法,得出可预报性最强的覆冰预报因子。
(1)初始化粒子群的位置矩阵:式(3)为粒子群的初始位置矩阵Ax,Ax中第i行、第j列的元素aij取值如式(4)所示。
A x = a 11 ... a 1 m . . . ... . . . a s 1 ... a s m - - - ( 3 )
a i j = 1 , r a n d ( 0 , 1 ) &GreaterEqual; 0.5 0 , r a n d ( 0 , 1 ) < 0.5 - - - ( 4 )
式(3)中,s为粒子群中粒子的个数,m为初步因子个数,即Ax的每一行对应一个粒子,每一列对应一个初步因子;式(4)中rand(0,1)为按均匀分布函数在区间[0,1]中随机抽取一个数。
(2)初始化粒子群的速度矩阵:粒子群的初始速度矩阵为式(5)所示,Vx中第i行、第j列的元素vij取值如式(6)所示。
V x = v 11 ... v 1 m . . . ... . . . v s 1 ... v s m - - - ( 5 )
vij=rand(0,1)(6)
(3)计算每个粒子的适应度的值,适应度值采用式(7)进行计算:
f x ( A x ( k ) ) = &Sigma; i ( e i ( k ) 1 - h i ( k ) ) 2 + N ( k ) 10000 - - - ( 7 )
式中fx(Ax(k))为第k个粒子的适应度;ei(k)为利用第k个粒子预测得到第i个样本的残差;hi(k)为第k个粒子预测得到第i个样本到样本中心点的广义距离;N(k)为第k个粒子所选取的因子的个数,即Ax第k行元素取1的个数。
(4)获得初始粒子群的最优位置和最优适应度,其最优适应度的计算如式(8)所示,
gfxmax=max(fx(Ax(k))),k=1,2,…s(8)
最优适应度取值所对应的粒子Ax(k)即为粒子群的最优位置,且为防止算法陷入局部最优,若某个粒子的位置等于种群的最优位置,将该粒子的位置重新初始化赋值。
(5)速度更新:采用式(9)对粒子群的速度矩阵进行更新,
Vxnew=w*Vx+c1*R1.*(Axbest-Ax)+c2*R2.*(Axgbest-Ax)(9)
式中Vxnew为更新后的速度矩阵;w为惯性因子,本文取0.8;c1、c2为加速因子,本文取0.7;R1、R2为与Ax同大小、元素为[0,1]的随机矩阵;Axbest为每个粒子的最优位置的组合矩阵;Axgbest为粒子群的最优位置矩阵,每行元素相同,均为粒子群的最优位置。
(6)位置更新:采用式(10)对粒子群的位置矩阵进行更新,
Axnew=round(Ax+Vxnew)(10)
式中Axnew为粒子群更新后的位置矩阵;round为一函数,即当矩阵中的元素小于0.5时取0,否则取1;
(7)判断是否达到迭代次数,或所有粒子取值均为最优粒子。若是,则优化结束,得到粒子群的最优位置,即选出的因子为可预报性最强的覆冰预报因子。否则,返回步骤(3),重新计算。
4.以选出的预报因子为基础利用历史数据建立最小二乘回归预测模型。
(1)建模:利用所选因子建立多元线性回归模型,如式(11)所示。
yt=b0+b1Xt1+b2Xt2+…+bNXtN+ut(11)
式中yt是输出,即为平均覆冰日数;b0,b1,…,bN为模型参数;Xt1,…XtN为所选的预报因子;ut为随机误差;N为所选预报因子的个数。
(2)求取模型参数:利用最小二乘方法求取模型参数,如式(12)所示
b=(XTX)-1XTY(12)
式中 b = b 0 b 1 . . . b N , X = 1 x 11 ... x 1 N 1 x 21 ... x 2 N . . . . . . . . . 1 x n 1 ... x n N , y = y 1 y 2 . . . y n ; n为样本数。
5.利用所建模型对冬季电网覆冰发生程度进行定量预报。
将当前预报因子的值带入式(11)求解得到未来冬季电网的平均覆冰日数。
下面以一个实例对上述方法进行说明。
2013年基于离散粒子群和最小二乘回归的湖南电网冬季覆冰定量长期预报分析:
1.收集历史环流指数数据与覆冰资料。
收集自1951年1月至2013年10月以来每月的74项环流指数,包括亚洲极涡、西太平洋副热带高压等指数。通过气象部门收集湖南历年平均覆冰日数。
平均覆冰日数为:湖南各个覆冰监测站点的覆冰日数相累加除以监测站数。
2.计算环流指数与平均覆冰日数的相关系数,初步选择覆冰预报因子。
计算3-10月的74项环流指数与湖南历年平均覆冰日数的相关系数和每项环流指数X与平均覆冰日数的最大相关系数。取所得RXmax绝对值最大的前17个因子,作为覆冰预报的初步因子。
3.利用离散粒子群算法,得出可预报性最强的覆冰预报因子。
设定粒子种群规模即粒子总数为15,迭代次数为100,进行离散粒子群初始化计算,迭代更新粒子的速度和位置,最终优化得到的可预报性最强的覆冰预报因子。
4.以选出的预报因子为基础利用历史数据建立最小二乘回归预测模型。
以所选出的可预报性最强的因子进行湖南电网冬季覆冰预测建模,利用最小二乘方法求取模型参数b。
5.利用所建模型对冬季电网覆冰发生程度进行定量预报。
将所选因子对应的值带入所建的回归预测模型,计算得到2013年冬季(即2013年12月~2014年2月)湖南电网平均覆冰日数4.4天,实际为4.7天,准确度大于90%。
本发明中基于离散粒子群和最小二乘回归预测的算法,相比传统的最小二乘回归预测算法,更适用于电网冬季覆冰长期预测。由于采用离散粒子群对覆冰预报因子进行了优化选择,提高了回归算法的预测精度。

Claims (6)

1.一种基于离散粒子群和最小二乘的电网覆冰定量长期预报方法,其特征在于包括:
计算每个月的74项环流指数与预报区域的历年平均覆冰日数的相关系数,并求取每项环流指数与平均覆冰日数的最大相关系数;
取所得最大相关系数绝对值最大的前m个因子,作为覆冰预报的初步因子;
利用离散粒子群算法,得到粒子群的最优位置,在初步因子中选出可预报性最强的覆冰预报因子;
以选出的可预报性最强的覆冰预报因子为基础,建立预测模型:
yt=b0+b1Xt1+b2Xt2+…+bNXtN+ut
式中yt为平均覆冰日数;b0,b1,…,bN为模型参数;Xt1,…XtN为所选的预报因子;ut为随机误差;N为所选预报因子的个数;
利用最小二乘方法求取模型参数,并将当前预报因子的值带入预测模型求解得到未来冬季电网的平均覆冰日数。
2.根据权利要求1所述的基于离散粒子群和最小二乘的电网覆冰定量长期预报方法,其特征在于所述m取值为8~25。
3.根据权利要求1所述的基于离散粒子群和最小二乘的电网覆冰定量长期预报方法,其特征在于利用下式计算3-10月份的74项环流指数与平均覆冰日数的相关系数:
R X j = n &Sigma; i = 1 n x i j d i - &Sigma; i = 1 n x i j &CenterDot; &Sigma; i = 1 n d i n &Sigma; i = 1 n x i j 2 - ( &Sigma; i = 1 n x i j ) 2 &CenterDot; n &Sigma; i = 1 n d i 2 - ( &Sigma; i - 1 n d i ) 2
式中,为第j月74项环流指数X与平均覆冰日数的相关系数;xij为74项环流指数的值,i为年份序号,j为月份序号,n为历史年份数量;di为第i年的平均覆冰日数;
利用式下式计算每项环流指数X与平均覆冰日数的最大相关系数,
R X m a x = m a x ( R X 3 , R X 4 , ... , R X 10 ) .
4.根据权利要求1所述的基于离散粒子群和最小二乘的电网覆冰定量长期预报方法,其特征在于所述利用离散粒子群算法包括:
(1)初始化粒子群的位置矩阵:式(3)为粒子群的初始位置矩阵Ax,Ax中第i行、第j列的元素aij取值如式(4);
A x = a 11 ... a 1 m . . . ... . . . a s 1 ... a s m - - - ( 1 )
a i j = 1 , r a n d ( 0 , 1 ) &GreaterEqual; 0.5 0 , r a n d ( 0 , 1 ) < 0.5 - - - ( 2 )
式(1)中,s为粒子群中粒子的个数,m为初步因子个数;式(2)中rand(0,1)为按均匀分布函数在区间[0,1]中随机抽取一个数;
(2)初始化粒子群的速度矩阵:粒子群的初始速度矩阵为式(5)所示,Vx中第i行、第j列的元素vij取值如式(6)所示;
V x = v 11 ... v 1 m . . . ... . . . v s 1 ... v s m - - - ( 3 )
vij=rand(0,1)(4)
(3)计算每个粒子的适应度的值,适应度值采用式(7)进行计算:
f x ( A x ( k ) ) = &Sigma; i ( e i ( k ) 1 - h i ( k ) ) 2 + N ( k ) 10000 - - - ( 5 )
式中fx(Ax(k))为第k个粒子的适应度;ei(k)为利用第k个粒子预测得到第i个样本的残差;hi(k)为第k个粒子预测得到第i个样本到样本中心点的广义距离;N(k)为第k个粒子所选取的因子的个数;
(4)获得初始粒子群的最优位置和最优适应度,其最优适应度的计算如式(8)所示,
gfxmax=max(fx(Ax(k))),k=1,2,…s(6)最优适应度取值所对应的粒子Ax(k)即为粒子群的最优位置,若某个粒子的位置等于种群的最优位置,将该粒子的位置重新初始化赋值;
(5)速度更新:采用式(9)对粒子群的速度矩阵进行更新,
Vxnew=w*Vx+c1*R1.*(Axbest-Ax)+c2*R2.*(Axgbest-Ax)(7)
式中Vxnew为更新后的速度矩阵;w为惯性因子;c1、c2为加速因子;R1、R2为与Ax同大小、元素为[0,1]的随机矩阵;Axbest为每个粒子的最优位置的组合矩阵;Axgbest为粒子群的最优位置矩阵,每行元素相同,均为粒子群的最优位置;
(6)位置更新:采用式(10)对粒子群的位置矩阵进行更新,
Axnew=round(Ax+Vxnew)(8)
式中Axnew为粒子群更新后的位置矩阵;round为一函数,即当矩阵中的元素小于0.5时取0,否则取1;
(7)判断是否达到迭代次数,或所有粒子取值均为最优粒子;若是,则优化结束,得到粒子群的最优位置,即选出的因子为可预报性最强的覆冰预报因子;否则,返回步骤(3),重新计算。
5.根据权利要求4所述的基于离散粒子群和最小二乘的电网覆冰定量长期预报方法,其特征在于设定粒子总数为10~30。
6.根据权利要求1所述的基于离散粒子群和最小二乘的电网覆冰定量长期预报方法,其特征在于利用最小二乘方法求取模型参数,如式(12)所示
b=(XTX)-1XTY(9)
式中 b = b 0 b 1 . . . b N , X = 1 x 11 ... x 1 N 1 x 21 ... x 2 N . . . . . . . . . 1 x n 1 ... x n N , y = y 1 y 2 . . . y n ; n为样本数。
CN201510579268.XA 2015-09-11 2015-09-11 基于离散粒子群和最小二乘的电网覆冰定量长期预报方法 Active CN105139094B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510579268.XA CN105139094B (zh) 2015-09-11 2015-09-11 基于离散粒子群和最小二乘的电网覆冰定量长期预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510579268.XA CN105139094B (zh) 2015-09-11 2015-09-11 基于离散粒子群和最小二乘的电网覆冰定量长期预报方法

Publications (2)

Publication Number Publication Date
CN105139094A true CN105139094A (zh) 2015-12-09
CN105139094B CN105139094B (zh) 2016-09-07

Family

ID=54724435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510579268.XA Active CN105139094B (zh) 2015-09-11 2015-09-11 基于离散粒子群和最小二乘的电网覆冰定量长期预报方法

Country Status (1)

Country Link
CN (1) CN105139094B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107092982A (zh) * 2017-04-11 2017-08-25 北京国网富达科技发展有限责任公司 一种电网覆冰预报方法及装置
CN111104754A (zh) * 2019-12-31 2020-05-05 国家电网有限公司 基于粒子群算法的覆冰参数计算方法
CN112114384A (zh) * 2020-08-27 2020-12-22 中国南方电网有限责任公司超高压输电公司检修试验中心 一种输电线覆冰发生概率预报方法
CN116864124A (zh) * 2023-07-06 2023-10-10 深圳市绿航星际太空科技研究院 一种基于大数据的营养健康预测系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012002702A (ja) * 2010-06-17 2012-01-05 Central Res Inst Of Electric Power Ind 送電線部分模型の弾性支持方法及び弾性支持装置
CN103473610A (zh) * 2013-09-06 2013-12-25 国家电网公司 基于数字高程模型的电网覆冰增长趋势预测方法
CN104281883A (zh) * 2014-09-16 2015-01-14 国家电网公司 一种配电网架空线敷冰风险指数预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012002702A (ja) * 2010-06-17 2012-01-05 Central Res Inst Of Electric Power Ind 送電線部分模型の弾性支持方法及び弾性支持装置
CN103473610A (zh) * 2013-09-06 2013-12-25 国家电网公司 基于数字高程模型的电网覆冰增长趋势预测方法
CN104281883A (zh) * 2014-09-16 2015-01-14 国家电网公司 一种配电网架空线敷冰风险指数预测方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107092982A (zh) * 2017-04-11 2017-08-25 北京国网富达科技发展有限责任公司 一种电网覆冰预报方法及装置
CN111104754A (zh) * 2019-12-31 2020-05-05 国家电网有限公司 基于粒子群算法的覆冰参数计算方法
CN111104754B (zh) * 2019-12-31 2023-06-02 国家电网有限公司 基于粒子群算法的覆冰参数计算方法
CN112114384A (zh) * 2020-08-27 2020-12-22 中国南方电网有限责任公司超高压输电公司检修试验中心 一种输电线覆冰发生概率预报方法
CN116864124A (zh) * 2023-07-06 2023-10-10 深圳市绿航星际太空科技研究院 一种基于大数据的营养健康预测系统

Also Published As

Publication number Publication date
CN105139094B (zh) 2016-09-07

Similar Documents

Publication Publication Date Title
Han et al. A PV power interval forecasting based on seasonal model and nonparametric estimation algorithm
CN110070226A (zh) 基于卷积神经网络与元学习的光伏功率预测方法及系统
CN105139094B (zh) 基于离散粒子群和最小二乘的电网覆冰定量长期预报方法
CN104123682A (zh) 一种基于气象影响因素的配网故障风险评估方法
CN103198235B (zh) 基于风功率纵向时刻概率分布的风功率预测值预评估方法
CN105069527A (zh) 一种基于数据挖掘技术的台区合理线损预测方法
CN105117988A (zh) 一种电力系统中缺失数据插补的方法
CN109086928A (zh) 基于saga-fcm-lssvm模型的光伏电站实时功率预测方法
CN104463701A (zh) 一种配电系统与电动汽车充电网络的协调规划方法
CN104376371B (zh) 一种基于拓扑的配网分层负荷预测方法
CN104600713A (zh) 含风力/光伏发电配电网日前无功调度的生成装置及方法
CN105160423A (zh) 一种基于马尔科夫残差修正的光伏发电预测方法
CN103441535A (zh) 基于sced的日前发电计划光伏发电接纳能力分析方法
CN105095668B (zh) 基于亚洲极涡因子的电网覆冰长期预报方法
CN103927460A (zh) 一种基于rbf的风电场短期风速预测方法
CN105226648A (zh) 一种基于大数据的分布式电源配电网规划方法
CN102968747A (zh) 一种风电场典型日出力曲线的确定方法
CN105512774A (zh) 接入分布式光伏发电的配电网检修计划优化方法
CN103166248B (zh) 一种独立风-柴-储微电网系统容量的工程配置方法
CN104899689A (zh) 基于dbcc优化算法与熵权理论的配电网故障恢复方法
CN104915727B (zh) 多维度同构异质bp神经网络光功率超短期预测方法
CN105529714A (zh) 基于正态分布组合特征的快速概率潮流计算方法
CN101709959A (zh) 一种输电线覆冰厚度的预测方法
CN109242337B (zh) 一种面向省级电网的变电站巡检机器人资源调配方法
CN105117801A (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
C14 Grant of patent or utility model
GR01 Patent grant