CN105260952A - 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 - Google Patents
基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 Download PDFInfo
- Publication number
- CN105260952A CN105260952A CN201510808313.4A CN201510808313A CN105260952A CN 105260952 A CN105260952 A CN 105260952A CN 201510808313 A CN201510808313 A CN 201510808313A CN 105260952 A CN105260952 A CN 105260952A
- Authority
- CN
- China
- Prior art keywords
- markov chain
- photovoltaic
- lambda
- state
- chain model
- 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
Links
Landscapes
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法,用于由N个发电单元组成的集中式光伏电站,所述各发电单元包括以下四个部件:光伏阵列、直流汇流箱、光伏逆变器、变压器,最后汇入交流电网;包括以下步骤:建立各元件马尔科夫链模型,状态抽样,运行过程模拟,判断算法收敛。本发明具有良好的稳定性,在大型光伏电站的可靠性评估中有着良好的正确性和优越性。
Description
技术领域
本发明涉及一种光伏电站可靠性评估方法,尤其是一种基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法,属于光伏发电技术领域。
背景技术
近年来,太阳能光伏发电技术持续快速发展,并网光伏发电的装机容量迅速增长。2013年,全球太阳能光伏新增容量第一次超过了风电新增容量,预计到2020年,我国累计光伏发电装机容量将突破20GW。然而,由于国内光伏项目长期存在抢装现象,故障多效率低已经成为不可忽视的问题,同时电站整体的长期可靠性也令人担忧,因此针对不同规模的光伏电站的可靠性评估变得尤为重要。目前国内对于光伏并网发电系统的并网研究较多,对于从工程可靠性角度出发的研究较少。对于光伏电站可靠性的评估有很多问题需要解决,其中选取合适且高效的评估方法是主要问题。
典型的集中式光伏电站由N(N>0)个发电单元组成,每个发电单元的结构如图1所示,包括以下五部件:1)光伏阵列。光伏阵列光伏组件通过串并联组成;2)直流汇流箱。减少了太阳能光伏阵列与光伏逆变器之间的连接线;3)光伏逆变器。可同时实现逆变、单位功率因数控制和最大功率跟踪控制等功能;4)变压器。光伏电站一般采用双分裂式变压器实现并网,可以保证两低压绕组所接入的逆变器电流独立汇入电网且互不影响;5)交流电网。
目前在光伏发电系统常用的可靠性评估方法有主要有解析法和模拟法。解析法,如条件概率方法将故障发生的情况以概率的方式表达出来,这种方法的不足是不能够描述光伏电站的时序特性,从而使结果产生比较大的偏差。并且随着所评估系统的规模增长,解析法很难处理故障状态的大量随机因素组合,容易出现具有“维数灾”的状态数目。模拟法主要是序贯蒙特卡洛法,这种方法既能描述光伏发电的时序出力,又能求解频率和持续时间等可靠性指标。其不足表现在计算效率低,维数高,并且具有静态特性,难以从高维的概率分布函数中抽样。
发明内容
本发明要解决的技术问题是提供一种基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法。
本发明采用下述技术方案:
一种基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法,用于由N个发电单元组成的集中式光伏电站,所述各发电单元包括以下部件:光伏阵列、直流汇流箱、光伏逆变器、变压器;其特征在于:
包括以下步骤:
步骤1:建立各元件马尔科夫链模型:筛选与光伏电站可靠性相关的元件为光伏阵列、直流汇流箱、光伏逆变器、变压器;所述直流汇流箱、光伏逆变器和变压器采用第一类马尔科夫链模型;所述光伏组件采用第二类马尔科夫链模型;
随机变量的序列{X(0),X(1),…,X(t)}表示元件在时刻0至时刻t的状态,下一时刻状态X(t+1)由当前时刻状态X(t)的条件分布P(X(t+1)|X(t))决定,与{X(0),X(1),…,X(t-1)}无关;
所述第一类马尔科夫链模型为:
所述直流汇流箱、光伏逆变器和变压器均有两个相互独立的状态:故障停运和正常运行;所述直流汇流箱、光伏逆变器和变压器的状态相互独立;t时刻状态X(t)为:
第一类马尔科夫链模型中元件的一步转移矩阵A1为:
其中,λ为元件的故障率,μ为元件的修复率;
所述第一类马尔科夫链模型的稳态概率分布p*为:
所述第二类马尔科夫链模型为:
当前时刻状态X(t)为:
第二类马尔科夫链模型中元件的一步转移矩阵A2为:
λ4,λ3,λ2,λ1分别为状态间的转化率,μ为元件的修复率;
所述第二类马尔科夫链的稳态分布概率为:
步骤2:状态抽样:一个包含了l个元件的系统,其状态可用向量S=(S1,S2,…,Sl)来表示:
对于采用第一类马尔科夫链模型的元件,其状态确定方法为:
其中,u为一个均匀分布在[0,1]区间的随机数,p0和p1分别为故障状态和正常状态下第i个元件下次抽样状态改变概率,由所述马尔科夫链模型确定。
对于采用第二类马尔科夫链模型的元件,其状态确定方法为:
其中,u为一个均匀分布在[0,1]区间的随机数,p0、p1、p2、p3和p4分别为故障状态和正常状态下第i个元件下次抽样状态改变概率,由所述马尔科夫链模型确定;
步骤3:利用步骤1中的马尔科夫链作为光伏电站的元件状态样本,重复步骤2的状态抽样过程,进行光伏电站的运行状态模拟,从而统计出所需数据,进行大型光伏电站的可靠性指标评估与分析;
步骤4:判断样本是否足够,如果是,输出可靠性评估指标,否则,转向步骤2;
所述判断样本是否足够的判据为:
其中,f为所求的可靠性评估指标,V(f)表示其估计值,为E(f)的估计值,E(f)为可靠性评估指标的期望,β为方差系数,其表达是为:
其中为的估计值。
采用上述技术方案所产生的有益效果在于:
1、本发明采用马尔科夫链建立了元件状态的先验分布,从而大大减小了评估时间,降低了运算的维数。
2、本发明对光伏阵列建立多状态划分模型,使得评估结果更加准确。
3、本发明建立了评估过程中产生的随机序列元素之间的联系,使评估更符合现实光伏电站运行情况。
4、本发明可以适应元件故障数据的动态变化,实现了对可靠性进行动态评估。
附图说明
图1是本发明中光伏电站发电单元的原理框图;
图2是本发明的流程图;
图3是本发明中第一类马尔科夫链的状态转移图;
图4是本发明中第二类马尔科夫链的状态转移图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
典型的集中式光伏电站由N个发电单元组成,每个发电单元的结构如图1所示,包括以下五个部件:光伏阵列、直流汇流箱、光伏逆变器、变压器、交流电网。
如图2所示,一种基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法,用于由N个发电单元组成的集中式光伏电站,所述各发电单元包括以下四个主要部件:光伏阵列、直流汇流箱、光伏逆变器、变压器;其特征在于:
包括以下步骤:
步骤1:建立各元件马尔科夫链模型:筛选与光伏电站可靠性相关的元件为光伏阵列、直流汇流箱、光伏逆变器、变压器;所述直流汇流箱、光伏逆变器和变压器采用第一类马尔科夫链模型;所述光伏组件采用第二类马尔科夫链模型;
随机变量的序列{X(0),X(1),…,X(t)}表示元件在时刻0至时刻t的状态,下一时刻状态X(t+1)由当前时刻状态X(t)的条件分布P(X(t+1)|X(t))决定,与{X(0),X(1),…,X(t-1)}无关,所生成的序列即为马尔科夫链。
所述第一类马尔科夫链模型为:
所述直流汇流箱、光伏逆变器和变压器均有两个相互独立的状态:故障停运和正常运行;所述直流汇流箱、光伏逆变器和变压器的状态相互独立;t时刻状态X(t)为:
第一类马尔科夫链模型中元件的一步转移矩阵A1为:
其中,λ为元件的故障率,μ为元件的修复率;
所述第一类马尔科夫链模型的稳态概率分布p*为:
所述第二类马尔科夫链模型为:
当前时刻状态X(t)为:
第二类马尔科夫链模型中元件的一步转移矩阵A2为:
λ4,λ3,λ2,λ1为图4的转化率,μ为元件的修复率;光伏阵列工作在可用度100%状态的情况包括可用度为95%-100%的所有情况,以此类推。而光伏阵列工作在可用度80%状态的情况则包括可用度为0%-80%的所有情况。
求解马尔科夫链稳态概率分布:对于第二类马尔科夫链,基于C-K方程有如下表达式。
pm(t)表示t时刻处于状态m的概率。在初始状态时,假设元件正常工作,即:
对稳态分布概率有
根据式(5)-(7)便可以计算出稳态分布概率p*。
步骤2:状态抽样:一个包含了l个元件的系统,其状态可用向量S=(S1,S2,…,Sl)来表示:
对于采用第一类马尔科夫链模型的元件,其状态确定方法为:
其中,u为一个均匀分布在[0,1]区间的随机数,p0和p1分别为故障状态和正常状态下第i个元件下次抽样状态改变概率,由所述马尔科夫链模型确定。
对于采用第二类马尔科夫链模型的元件,其状态确定方法为:
其中,u为一个均匀分布在[0,1]区间的随机数,p0、p1、p2、p3和p4分别为故障状态和正常状态下第i个元件下次抽样状态改变概率,由所述马尔科夫链模型确定;
步骤3:利用步骤1中的马尔科夫链作为光伏电站的元件状态样本,重复步骤2的状态抽样过程,进行光伏电站的运行状态模拟,从而统计出所需数据,进行大型光伏电站的可靠性指标评估与分析;
步骤4:判断样本是否足够,如果是,输出可靠性评估指标,否则,转向步骤2;
所述判断样本是否足够的判据为:
其中,f为所求的可靠性评估指标,V(f)表示其估计值,为E(f)的估计值,E(f)为可靠性评估指标的期望,β为方差系数,其表达是为:
其中为的估计值。
本发明首先对于光伏电站的结构进行划分,筛选出与光伏电站可靠性相关的元件,然后通过对每种类型元件的重复抽样,建立以平稳分布与其概率分布相同的马尔科夫链,从而得到元件的状态样本。基于这些元件的状态样本进行蒙特卡洛仿真,可以实现光伏电站运行情况的动态蒙特卡洛模拟,同时由于考虑到了随机序列之间的相关性,使模拟过程更加接近于光伏电站的真实运行状况,最后根据模拟的结果进行对于光伏电站可靠性的评估。此算法具有良好的稳定性,收敛速度比传统的蒙特卡罗方法快,在大型光伏电站的可靠性评估中有着良好的正确性和优越性。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (1)
1.一种基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法,用于由N个发电单元组成的集中式光伏电站,所述各发电单元包括以下部件:光伏阵列、直流汇流箱、光伏逆变器、变压器;其特征在于:
包括以下步骤:
步骤1:建立各元件马尔科夫链模型:筛选与光伏电站可靠性相关的元件为光伏阵列、直流汇流箱、光伏逆变器、变压器;所述直流汇流箱、光伏逆变器和变压器采用第一类马尔科夫链模型;所述光伏组件采用第二类马尔科夫链模型;
随机变量的序列{X(0),X(1),…,X(t)}表示元件在时刻0至时刻t的状态,下一时刻状态X(t+1)由当前时刻状态X(t)的条件分布P(X(t+1)|X(t))决定,与{X(0),X(1),…,X(t-1)}无关;
所述第一类马尔科夫链模型为:
所述直流汇流箱、光伏逆变器和变压器均有两个相互独立的状态:故障停运和正常运行;所述直流汇流箱、光伏逆变器和变压器的状态相互独立;t时刻状态X(t)为:
第一类马尔科夫链模型中元件的一步转移矩阵A1为:
所述第一类马尔科夫链模型的稳态概率分布p*为:
所述第二类马尔科夫链模型为:
当前时刻状态X(t)为:
第二类马尔科夫链模型中元件的一步转移矩阵A2为:
λ4,λ3,λ2,λ1分别为状态间的转化率,μ为元件的修复率;
所述第二类马尔科夫链的稳态分布概率为:
步骤2:状态抽样:一个包含了l个元件的系统,其状态可用向量S=(S1,S2,…,Sl)来表示:
对于采用第一类马尔科夫链模型的元件,其状态确定方法为:
其中,u为一个均匀分布在[0,1]区间的随机数,p0和p1分别为故障状态和正常状态下第i个元件下次抽样状态改变概率,由所述马尔科夫链模型确定。
对于采用第二类马尔科夫链模型的元件,其状态确定方法为:
其中,u为一个均匀分布在[0,1]区间的随机数,p0、p1、p2、p3和p4分别为故障状态和正常状态下第i个元件下次抽样状态改变概率,由所述马尔科夫链模型确定;
步骤3:利用步骤1中的马尔科夫链作为光伏电站的元件状态样本,重复步骤2的状态抽样过程,进行光伏电站的运行状态模拟,从而统计出所需数据,进行大型光伏电站的可靠性指标评估与分析;
步骤4:判断样本是否足够,如果是,输出可靠性评估指标,否则,转向步骤2;
所述判断样本是否足够的判据为:
其中,f为所求的可靠性评估指标,V(f)表示其估计值,为E(f)的估计值,E(f)为可靠性评估指标的期望,β为方差系数,其表达是为:
其中为的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510808313.4A CN105260952B (zh) | 2015-11-19 | 2015-11-19 | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510808313.4A CN105260952B (zh) | 2015-11-19 | 2015-11-19 | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105260952A true CN105260952A (zh) | 2016-01-20 |
CN105260952B CN105260952B (zh) | 2018-12-14 |
Family
ID=55100627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510808313.4A Expired - Fee Related CN105260952B (zh) | 2015-11-19 | 2015-11-19 | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260952B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106980918A (zh) * | 2017-02-16 | 2017-07-25 | 广西电网有限责任公司电力科学研究院 | 一种发输电系统可靠性评估系统 |
CN106992513A (zh) * | 2017-02-16 | 2017-07-28 | 广西电网有限责任公司电力科学研究院 | 一种发输电系统可靠性评估方法 |
CN107067130A (zh) * | 2016-12-12 | 2017-08-18 | 浙江大学 | 一种基于电动汽车马尔可夫充电需求分析模型的快速充电站容量规划方法 |
CN108122049A (zh) * | 2017-12-21 | 2018-06-05 | 郓城金河热电有限责任公司 | 基于光伏组件故障率分析的光伏阵列运行维护方法 |
CN108537413A (zh) * | 2018-03-19 | 2018-09-14 | 国网天津市电力公司 | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 |
CN109062868A (zh) * | 2018-06-27 | 2018-12-21 | 北京航空航天大学 | 一种一般分布下的多态系统状态稳态概率求解方法 |
CN112100821A (zh) * | 2020-08-26 | 2020-12-18 | 西北工业大学 | 光伏电池的稳健性优化设计方法 |
CN113610124A (zh) * | 2021-07-23 | 2021-11-05 | 中国地质大学(武汉) | 基于马尔可夫链蒙特卡洛的人手轨迹生成方法和系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091238A (zh) * | 2014-07-11 | 2014-10-08 | 国家电网公司 | 一种恶劣天气下用户用电安全风险演化分析评估方法 |
CN104218620A (zh) * | 2014-09-26 | 2014-12-17 | 国家电网公司 | 基于伪序贯蒙特卡洛模拟的主动配电网可靠性分析方法 |
CN104239046A (zh) * | 2014-09-05 | 2014-12-24 | 河海大学 | 基于隐马尔科夫模型和多目标进化算法的软件自适应方法 |
-
2015
- 2015-11-19 CN CN201510808313.4A patent/CN105260952B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091238A (zh) * | 2014-07-11 | 2014-10-08 | 国家电网公司 | 一种恶劣天气下用户用电安全风险演化分析评估方法 |
CN104239046A (zh) * | 2014-09-05 | 2014-12-24 | 河海大学 | 基于隐马尔科夫模型和多目标进化算法的软件自适应方法 |
CN104218620A (zh) * | 2014-09-26 | 2014-12-17 | 国家电网公司 | 基于伪序贯蒙特卡洛模拟的主动配电网可靠性分析方法 |
Non-Patent Citations (4)
Title |
---|
M THERISTIS 等: "Markovian Reliability Analysis of Standalone Photovoltaic Systems Incorporating Repairs", 《IEEE JOURNAL OF PHOTOVOLTAICS》 * |
RADWA SAYED 等: "Modeling of Photovoltaic Based Power Stations for Reliability Studies Using Markov Chains", 《INTERNATIONAL CONFERENCE ON RENEWABLE ENERGY RESEARCH AND APPLICATIONS》 * |
戴志辉 等: "基于动态故障树与蒙特卡罗仿真的保护系统动态可靠性评估", 《中国电机工程学报》 * |
王震 等: "分布式光伏发电系统的可靠性模型及指标体系", 《电力系统自动化》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107067130B (zh) * | 2016-12-12 | 2021-05-07 | 浙江大学 | 一种基于电动汽车马尔可夫充电需求分析模型的快速充电站容量规划方法 |
CN107067130A (zh) * | 2016-12-12 | 2017-08-18 | 浙江大学 | 一种基于电动汽车马尔可夫充电需求分析模型的快速充电站容量规划方法 |
CN106992513A (zh) * | 2017-02-16 | 2017-07-28 | 广西电网有限责任公司电力科学研究院 | 一种发输电系统可靠性评估方法 |
CN106980918A (zh) * | 2017-02-16 | 2017-07-25 | 广西电网有限责任公司电力科学研究院 | 一种发输电系统可靠性评估系统 |
CN108122049A (zh) * | 2017-12-21 | 2018-06-05 | 郓城金河热电有限责任公司 | 基于光伏组件故障率分析的光伏阵列运行维护方法 |
CN108537413A (zh) * | 2018-03-19 | 2018-09-14 | 国网天津市电力公司 | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 |
CN108537413B (zh) * | 2018-03-19 | 2021-09-10 | 国网天津市电力公司 | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 |
CN109062868B (zh) * | 2018-06-27 | 2021-09-28 | 北京航空航天大学 | 一种一般分布下的多态系统状态稳态概率求解方法 |
CN109062868A (zh) * | 2018-06-27 | 2018-12-21 | 北京航空航天大学 | 一种一般分布下的多态系统状态稳态概率求解方法 |
CN112100821A (zh) * | 2020-08-26 | 2020-12-18 | 西北工业大学 | 光伏电池的稳健性优化设计方法 |
US20220067258A1 (en) * | 2020-08-26 | 2022-03-03 | Northwestern Polytechnical University | Robust optimal design method for photovoltaic cells |
US11630936B2 (en) * | 2020-08-26 | 2023-04-18 | Northwestern Polytechnical University | Robust optimal design method for photovoltaic cells |
CN113610124A (zh) * | 2021-07-23 | 2021-11-05 | 中国地质大学(武汉) | 基于马尔可夫链蒙特卡洛的人手轨迹生成方法和系统 |
CN113610124B (zh) * | 2021-07-23 | 2024-04-19 | 中国地质大学(武汉) | 基于马尔可夫链蒙特卡洛的人手轨迹生成方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105260952B (zh) | 2018-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105260952A (zh) | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 | |
Basaran et al. | Energy management for on‐grid and off‐grid wind/PV and battery hybrid systems | |
Li et al. | A multi-state model for the reliability assessment of a distributed generation system via universal generating function | |
Gao et al. | Adequacy assessment of generating systems containing wind power considering wind speed correlation | |
CN104504607A (zh) | 一种基于模糊聚类算法的光伏电站故障诊断方法 | |
CN104037790B (zh) | 一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法 | |
CN107465211B (zh) | 孤岛微电网的分布式固定时间协调控制方法 | |
CN104009484B (zh) | 一种基于潮流计算的电网降损评估方法 | |
CN110601204B (zh) | 基于随机变量状态时序模拟的光伏并网系统概率潮流分析方法 | |
Talei et al. | Smart campus microgrid: Advantages and the main architectural components | |
CN109636027A (zh) | 一种基于蒙特卡罗法的多能互补系统供能可靠性评估方法 | |
CN102570451A (zh) | 一种输电网静态可靠性评估方法 | |
CN110363334A (zh) | 基于灰色神经网络模型的光伏并网的电网线损预测方法 | |
CN103293950B (zh) | 一种基于lssvm的光伏最大功率点跟踪的控制方法 | |
Ferdowsi et al. | New monitoring approach for distribution systems | |
CN105552880A (zh) | 基于状态枚举法的电力系统典型故障集确定方法 | |
CN108110789A (zh) | 一种间歇性可再生能源分层分区并网规划方法 | |
Durgadevi et al. | Availability Capacity Evaluation and Reliability Assessment of Integrated Systems Using Metaheuristic Algorithm. | |
Mandal et al. | Assessment of Resiliency by incorporating DGs in Power Network using Graph Theoretic Approach | |
Beyza et al. | Evaluation of Reliability and Robustness of Electric Power Systems with Renewable Energies | |
CN104283208A (zh) | 大型电网概率性可用输电能力的分解协调计算方法 | |
CN103326356A (zh) | 考虑人工调度知识的电网关键断面极限传输容量计算方法 | |
Liu et al. | Reliability assessment of micro-grid based on Markov model | |
Dhua et al. | Optimization of Generation Capacity at the Incoming Microgrid in an interconnected Microgrid System using ANN | |
CN103280841A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181214 Termination date: 20191119 |