CN107562992B - 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 - Google Patents
一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 Download PDFInfo
- Publication number
- CN107562992B CN107562992B CN201710609495.1A CN201710609495A CN107562992B CN 107562992 B CN107562992 B CN 107562992B CN 201710609495 A CN201710609495 A CN 201710609495A CN 107562992 B CN107562992 B CN 107562992B
- Authority
- CN
- China
- Prior art keywords
- maximum power
- photovoltaic array
- svm
- photovoltaic
- voltage value
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000004422 calculation algorithm Methods 0.000 title claims description 21
- 238000005286 illumination Methods 0.000 claims abstract description 31
- 238000004088 simulation Methods 0.000 claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 238000012549 training Methods 0.000 claims abstract description 10
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000012795 verification Methods 0.000 claims abstract description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 claims 1
- 238000012706 support-vector machine Methods 0.000 description 20
- 238000010248 power generation Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 7
- 238000003491 array Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process 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
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
Landscapes
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种光伏列阵最大功率跟踪方法,包括以下步骤:S10,建立光伏列阵最大功率跟踪预测模型,对确定规模的光伏列阵进行大量仿真验证,获得不同阴影和温度下光伏列阵的最大功率值及对应电压值;S20,将步骤S10各光伏电池的光照、温度、该条件下光伏列阵的最大功率值及其对应的电压值数据作为训练样本和测试样本,建立和训练SVM最大功率跟踪模型,仿真测试;S30,根据步骤S20得到的SVM最大功率跟踪模型预测出待跟踪条件下最大功率点对应电压值的范围,通过粒子群算法对电压值进行寻优,获得最优电压值及其对应的全局最大功率点。所述方法能够在较短时间内,精确获得任意阴影下光伏列阵的最大功率点,有效提高了光伏电源的输出效率。
Description
技术领域
本发明涉及电力系统仿真领域,具体涉及一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法。
背景技术
随着化石能源的逐渐减少和新能源技术的不断创新,光伏发电和风力发电得到了长足的发展。其中光伏发电由于光能的节能环保,随地使用及成本低,分布式光伏发电逐渐走进千家万户,人们在自家楼顶安装光伏发电装置,通过光伏发电满足自家的用电需求,从而大大的减少了通过燃烧化石能源获得的电力,既节能又环保。但在光伏发电的使用过程中,光伏发电是否能够最大限度的输出功率一直是人们关心的话题,即出现了光伏列阵的最大功率跟踪(MPPT)问题。MPPT是指在光伏列阵温度或者光照发生变化时,MPPT装置能够快速的获取当前条件下的最大功率点。
光伏列阵的最大功率跟踪难点在于,由于外界条件的快速变化,当前的MPPT跟踪技术无法快速的获得相应条件下的最大功率点。其原因在于:在光伏发电的早期,光伏电站往往建立在空旷,阳光充足的大西北地区,几乎整个电站的光伏列阵所受到的光照和温度条件是一样的,当时人们提出了扰动法,恒定电压跟踪法等方法均能够快速的获得光伏列阵最大功率点。但随着光伏发电的快速发展,分布式光伏发电逐渐走入城市,走入社区住宅,由于树的遮挡、云层变化等,使光伏列阵受到的光照不一致,从而产生了局部阴影下光伏列阵的最大功率点跟踪问题。在局部阴影下,光伏列阵各处所受到的光照和温度是不一样的,此时光伏列阵的P-U曲线呈现多极值性。传统的MPPT跟踪方法无法获得全局的极值点,无法获得光伏列阵在局部阴影下的最大功率点,导致光伏列阵无法最大功率输出,造成了光资源的浪费。
因此,局部阴影下的光伏列阵最大功率点跟踪是当前研究的重点,相关学者提出了一些改进型功率跟踪方法。但这些算法不是无法适用于并联的光伏列阵,就是不能保证在任意情况下得到全局最大功率点,也有些算法未对粒子的初始位置进行讨论,容易陷入局部最优值。
基于以上分析,本发明提出了一种基于SVM和粒子群算法相结合的最大功率跟踪方法,能够有效的获得局部阴影下光伏列阵的最大功率点。
发明内容
本发明的目的是针对上述现有技术的不足,提供了一种基于SVM对电压值进行初步预测,再采用粒子群算法进行寻优的光伏列阵最大功率跟踪方法。
本发明的目的可以通过如下技术方案实现:
一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,所述方法包括以下步骤:
步骤S10,建立光伏列阵最大功率跟踪预测模型,对确定规模的光伏列阵进行大量仿真验证,获得不同光照和温度下光伏列阵的最大功率值及其对应的电压值;
步骤S20,将步骤S10各光电池的光照、温度、该光照温度条件下光伏列阵的最大功率值及其对应的电压值数据作为训练样本和测试样本,建立和训练SVM最大功率跟踪模型,仿真测试;
步骤S30,根据步骤S20得到的SVM最大功率跟踪模型预测出待跟踪条件下最大功率点对应电压值的范围,通过粒子群算法对电压值进行寻优,获得最优电压值及其对应的全局最大功率点。
进一步地,所述建立光伏列阵最大功率跟踪预测模型包括以下步骤:
步骤S101,首先确定光伏列阵的规模,即构成光伏列阵的光伏电池数量,并搭建光伏列阵仿真平台;
步骤S102,分析影响光伏列阵输出功率的因素,从而建立相应的光伏列阵最大功率跟踪预测模型。
进一步地,步骤S101中,通过对光伏列阵设置不同的光照和温度情况,进行仿真,获得该条件下的最大功率点及其对应的光伏列阵电压,记录每个光电池所受光照、温度,及整个光伏列阵最大输出功率值和输出该功率值时光伏列阵两端的电压值;步骤S102中,所述影响光伏列阵输出功率的因素包括各光电池的光照和温度情况、光伏列阵的最大功率点及其对应的光伏列阵两端的电压值,以各影响因素作为输入量,建立和训练SVM最大功率跟踪模型。
进一步地,步骤S20中,在建立和训练SVM最大功率跟踪模型前,首先要对训练样本的数据进行归一化处理。
进一步地,步骤S30中,所述利用SVM最大功率跟踪模型预测待跟踪条件下最大功率点对应电压值的范围时,采用核函数代替点积运算来求解对偶问题和非线性映射到高维特征空间的问题,同时采用粒子群算法,通过个体极值和群体极值更新个体位置,通过多次迭代获得最优值。
进一步地,所述核函数包括:
线型核函数:K(x,xi)=xxi (1)
D阶多项式核函数:K(x,xi)=(xxi+1)d (2)
Sigmoid核函数K(x,xi)=tanh(k(xxi))+θ (4)
其中,K(x,xi)表示核函数,x表示自变量,xi表示样本向量,即输入量,d表示多项式核函数的阶数,σ表示径向基核函数参数,θ表示Sigmoid核函数参数1,k(xxi)表示Sigmoid核函数参数2。
1、本发明利用支持向量机与粒子群算法相结合的方法,有效地解决了局部阴影下光伏列阵的最大功率跟踪问题,支持向量机以采用结构风险最小化为准则,在最小化样本点误差的同时,最小化结构风险,提高了模型的泛化能力,且没有数据维数的限制。本发明通过输入各光电池的光照和温度,预测出光伏最大功率点对应电压值的大概范围,然后通过粒子群算法进行寻优,最终获得最大功率点对应的确切电压,从而实现最大功率跟踪。
2、本发明采用支持向量机和粒子群相结合的方法,有效的解决了局部阴影下光伏列阵最大功率寻优难的问题,为局部阴影下光伏列阵的最大功率跟踪提供了一种新的有效方法。
附图说明
图1为本发明实施例一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法的流程图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例:
本实施例提供了一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,所述方法的流程图如图1所示,包括以下步骤:
步骤S10,建立光伏列阵最大功率跟踪预测模型,对确定规模的光伏列阵进行大量仿真验证,获得不同光照和温度下光伏列阵的最大功率值及其对应的电压值;
所述光伏列阵最大功率跟踪预测模型包括光伏电池、电压测量表、功率测量表、输入模块,其中输入模块包括光照强度输入模块和温度输入模块,对于光伏列阵最大功率预测模型,光伏电池的建模最为重要,本实施例采用文献《孙航,杜海江,季迎旭,等.适用不同尺度光伏阵列的数值建模方法[J].电力系统自动化,2014,38(16):35-40.》中的光伏电池数学模型,其中光伏电池的I-U特性为:
其中:
I、U分别为流过光伏电池的电流及光伏电池两端的电压;Isc、Uoc、Im、Um分别为光伏电池的短路电流、开路电压、最大功率点电流和最大功率点电压;
在完成光伏电池建模的基础上,确立需要测试的光伏列阵规模,本实施例的光伏列阵规模为{5×4},表示4列并联的光伏串,每列光伏串由5个光伏电池串联而成。通过输入模块,设置光伏列阵中各光电池的光照和温度,然后通过建立的光伏列阵最大功率跟踪预测模型进行仿真,获得该条件下的P-U输出曲线,其中P为光伏列阵输出功率,U为光伏列阵两端的电压,记录此次仿真中每个光电池的光照和温度情况,最大功率Pmax及其对应的电压值Uc。然后更改各光电池的光照和温度,重新仿真计算,再次记录相关数据,多次重复仿真与记录,直到得到足够多的数据。
对于以上的光照和温度设置,为了能够较好的模拟实际情况,温度设置范围为0℃~50℃,光照设置范围为0W/m2~1500W/m2。
步骤S20,将步骤S10各光电池的光照、温度、该光照温度条件下光伏列阵的最大功率值及其对应的电压值数据作为训练样本和测试样本,建立和训练SVM最大功率跟踪模型,仿真测试;
根据以上数据进行SVM最大功率跟踪模型的建立和训练,其基本基本原理如下:
设训练样本集S为:
S={(x1,y1),(x2,y2),…,(xk,yk)} (8)
其中k为样本长度,yi(i=1,2,…,k)的值为1或-1,分别表示类别一、类别二,本实施例为待预测的电压值;
设分类超平面y(x)为:
y(x)=ωx+b (9)
其中ω为分类超平面法向量,b为分类超平面截距,若公式(4)能够将样本正确的划分为两类,即可等效为:
对公式(10)求解可得最优的ω*和b*,并得到最终的最优分类超平面:
ybest(x)=sgn[ω*x+b*] (11)
其中,C为惩罚因子,表示错分样本的惩罚程度,在求解最优的ω*、b*及将原输入空间映射到高维的特征空间过程中,均需要进行点积计算,导致计算量增加。为了减小计算,引入核函数K(xi,xj)代替点积运算,则公式(2)变更为:
对于本实施例,x为各光电池光照、温度等输入量,预测量y为最大功率点对应的电压值。
粒子群算法的基本原理为:
设粒子总数为N,假设第i个粒子的位置为Ui,第i个粒子搜索到的历史最优位置为Ubest,i,整个粒子群搜索到的最优位置为Ubestg,假设第i个粒子的位置变化率为vi,则速度及位置更新公式为:
vi(t+1)=w×vi(t)+d1×rand×[Ubest,i(t)-Ui(t)]+d2×rand×[Ubestg(t)-Ui(t)](14)
Ui(t+1)=Ui(t)+vi(t+1)…1≤i≤N (15)
其中t为迭代次数;d1和d2称为学习因子,为正常数;rand为[0,1]的随机数;w为惯性权重。
步骤S30,根据步骤S20得到的SVM最大功率跟踪模型预测出待跟踪条件下最大功率点对应电压值的范围,通过粒子群算法对电压值进行寻优,获得最优电压值及其对应的全局最大功率点。
通过SVM预测出最大功率点对应电压值的大概范围后,运用粒子群算法在很小的范围内进行寻优获得最终的电压值,再将光伏列阵的两端电压值设置为寻优电压值,进行电路仿真即可得到待寻优的最大功率点。
以上所述,仅为本发明专利较佳的实施例,但本发明专利的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明专利所公开的范围内,根据本发明专利的技术方案及其发明专利构思加以等同替换或改变,都属于本发明专利的保护范围。
Claims (6)
1.一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,其特征在于,所述方法包括以下步骤:
步骤S10,建立光伏列阵最大功率跟踪预测模型,对确定规模的光伏列阵进行大量仿真验证,获得不同光照和温度下光伏列阵的最大功率值及其对应的电压值;所述光伏列阵最大功率跟踪预测模型输入为光照和温度,输出为最大功率点及其对应的光伏阵列电压;
步骤S20,将步骤S10各光电池的光照、温度、该光照温度条件下光伏列阵的最大功率值及其对应的电压值数据作为训练样本和测试样本,建立和训练SVM最大功率跟踪模型,仿真测试;所述SVM最大功率跟踪模型输入为步骤S10各光电池的光照、温度、该光照温度条件下光伏列阵的最大功率值及其对应的电压值,输出为最大功率点对应电压值的大概范围;
步骤S30,根据步骤S20得到的SVM最大功率跟踪模型预测出待跟踪条件下最大功率点对应电压值的范围,通过粒子群算法对电压值进行寻优,获得最优电压值及其对应的全局最大功率点。
2.根据权利要求1所述的一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,其特征在于,所述建立光伏列阵最大功率跟踪预测模型包括以下步骤:
步骤S101,首先确定光伏列阵的规模,即构成光伏列阵的光伏电池数量,并搭建光伏列阵仿真平台;
步骤S102,分析影响光伏列阵输出功率的因素,从而建立相应的光伏列阵最大功率跟踪预测模型。
3.根据权利要求2所述的一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,其特征在于:步骤S101中,通过对光伏列阵设置不同的光照和温度情况,进行仿真,获得该条件下的最大功率点及其对应的光伏列阵电压,记录每个光电池所受光照、温度,及整个光伏列阵最大输出功率值和输出该功率值时光伏列阵两端的电压值;步骤S102中,所述影响光伏列阵输出功率的因素包括各光电池的光照和温度情况、光伏列阵的最大功率点及其对应的光伏列阵两端的电压值,以各影响因素作为输入量,建立和训练SVM最大功率跟踪模型。
4.根据权利要求1所述的一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,其特征在于:步骤S20中,在建立和训练SVM最大功率跟踪模型前,首先要对训练样本的数据进行归一化处理。
5.根据权利要求1所述的一种基于SVM和粒子群算法的光伏列阵最大功率跟踪方法,其特征在于:步骤S30中,所述利用SVM最大功率跟踪模型预测待跟踪条件下最大功率点对应电压值的范围时,采用核函数代替点积运算来求解对偶问题和非线性映射到高维特征空间的问题,同时采用粒子群算法,通过个体极值和群体极值更新个体位置,通过多次迭代获得最优值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710609495.1A CN107562992B (zh) | 2017-07-25 | 2017-07-25 | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710609495.1A CN107562992B (zh) | 2017-07-25 | 2017-07-25 | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107562992A CN107562992A (zh) | 2018-01-09 |
CN107562992B true CN107562992B (zh) | 2020-04-28 |
Family
ID=60974596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710609495.1A Expired - Fee Related CN107562992B (zh) | 2017-07-25 | 2017-07-25 | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107562992B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109815618B (zh) * | 2019-02-18 | 2022-09-13 | 南昌航空大学 | 基于物理模型与粒子群算法的遮荫下光伏发电跟踪方法 |
CN110471487A (zh) * | 2019-07-19 | 2019-11-19 | 湖南工业大学 | 一种基于自适应粒子群算法在光伏阵列多峰值系统mppt的控制方法 |
CN112947665A (zh) * | 2021-02-26 | 2021-06-11 | 大连海事大学 | 光伏阵列在动态阴影遮挡条件下的最大功率跟踪方法 |
CN114510111B (zh) * | 2021-12-29 | 2023-09-12 | 北京华能新锐控制技术有限公司 | 部分遮阳光伏阵列的全局mppt控制方法与装置 |
CN115951755A (zh) * | 2023-02-06 | 2023-04-11 | 广芯微电子(广州)股份有限公司 | 基于pso-lstm的光伏最大功率点追踪方法及装置 |
CN116896116B (zh) * | 2023-09-11 | 2023-11-28 | 广州德姆达光电科技有限公司 | 基于人工智能的太阳能并网调控方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992248A (zh) * | 2015-07-07 | 2015-10-21 | 中山大学 | 一种微电网光伏电站发电量组合预测方法 |
CN106503850A (zh) * | 2016-11-01 | 2017-03-15 | 上海电机学院 | 一种用于风电功率预测的svm预测模型 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104753461B (zh) * | 2015-04-10 | 2017-04-12 | 福州大学 | 基于粒子群优化支持向量机的光伏发电阵列故障诊断与分类方法 |
CN106952183A (zh) * | 2017-04-24 | 2017-07-14 | 武汉大学 | 一种基于粒子群优化最小二乘支持向量机的短期负荷预测方法 |
-
2017
- 2017-07-25 CN CN201710609495.1A patent/CN107562992B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992248A (zh) * | 2015-07-07 | 2015-10-21 | 中山大学 | 一种微电网光伏电站发电量组合预测方法 |
CN106503850A (zh) * | 2016-11-01 | 2017-03-15 | 上海电机学院 | 一种用于风电功率预测的svm预测模型 |
Also Published As
Publication number | Publication date |
---|---|
CN107562992A (zh) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107562992B (zh) | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 | |
Wang et al. | Photovoltaic cell parameter estimation based on improved equilibrium optimizer algorithm | |
Jiang et al. | Parameter estimation of solar cells and modules using an improved adaptive differential evolution algorithm | |
Xu et al. | Parameter estimation of photovoltaic modules using a hybrid flower pollination algorithm | |
US20230035108A1 (en) | Photovoltaic cell parameter identification method based on improved equilibrium optimizer algorithm | |
Mellit et al. | Performance prediction of 20 kWp grid-connected photovoltaic plant at Trieste (Italy) using artificial neural network | |
CN110556820B (zh) | 用于确定能量系统操作场景的方法和设备 | |
Michael et al. | Short-term solar irradiance forecasting based on a novel Bayesian optimized deep Long Short-Term Memory neural network | |
Patro et al. | Mathematical modeling framework of a PV model using novel differential evolution algorithm | |
CN103489038A (zh) | 基于lm-bp神经网络的光伏超短期功率预测方法 | |
Lun et al. | An I–V model based on time warp invariant echo state network for photovoltaic array with shaded solar cells | |
Piliougine et al. | Multilayer perceptron applied to the estimation of the influence of the solar spectral distribution on thin-film photovoltaic modules | |
Sogabe et al. | Optimization of decentralized renewable energy system by weather forecasting and deep machine learning techniques | |
Ibrahim et al. | A hybrid wind driven-based fruit fly optimization algorithm for identifying the parameters of a double-diode photovoltaic cell model considering degradation effects | |
Omer et al. | Comparison of decision tree based ensemble methods for prediction of photovoltaic maximum current | |
CN103293950B (zh) | 一种基于lssvm的光伏最大功率点跟踪的控制方法 | |
CN111985678A (zh) | 一种光伏功率短期预测方法 | |
CN116632823A (zh) | 基于功率转换模型与多层感知机的短期光伏功率预测方法 | |
Ncir et al. | An intelligent improvement based on a novel configuration of artificial neural network model to track the maximum power point of a photovoltaic panel | |
Choulli et al. | Hybrid optimization based on the analytical approach and the particle swarm optimization algorithm (Ana-PSO) for the extraction of single and double diode models parameters | |
Dong et al. | Simultaneous operating temperature and output power prediction method for photovoltaic modules | |
Khan et al. | Output power prediction of a photovoltaic module through artificial neural network | |
Mazzeo et al. | A user-friendly and accurate machine learning tool for the evaluation of the worldwide yearly photovoltaic electricity production | |
Kamthania et al. | Determination of efficiency of hybrid photovoltaic thermal air collectors using artificial neural network approach for different PV technology | |
Omer et al. | Adaptive boosting and bootstrapped aggregation based ensemble machine learning methods for photovoltaic systems output current prediction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20200428 |
|
CF01 | Termination of patent right due to non-payment of annual fee |