CN103093079B - 基于仿射算术的辐射型配电网不确定潮流分析方法 - Google Patents
基于仿射算术的辐射型配电网不确定潮流分析方法 Download PDFInfo
- Publication number
- CN103093079B CN103093079B CN201210558626.5A CN201210558626A CN103093079B CN 103093079 B CN103093079 B CN 103093079B CN 201210558626 A CN201210558626 A CN 201210558626A CN 103093079 B CN103093079 B CN 103093079B
- Authority
- CN
- China
- Prior art keywords
- node
- epsiv
- affine
- voltage
- sigma
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种基于仿射算术的辐射型配电网不确定潮流分析方法,其步骤包括:遍历整个网络,确定节点个数并为节点编号;将每个节点的负荷波动范围转换为仿射量;给网络中节点的负荷及电压赋初值;运用仿射算术的运算法则,从网络的末端开始不断前推回代,直至满足收敛精度的要求,得到各节点电压的波动范围。根据本发明的不确定潮流分析方法,可以有效地缓解前人方法的保守性,能够得到更加小的电压波动范围并使其包含所有的电压可能解。
Description
技术领域
本发明属于电力系统稳态分析领域,涉及一种电力系统不确定潮流分析方法,更具体地,涉及一种基于仿射算术的辐射型配电网不确定潮流分析方法。
背景技术
潮流分析是根据给定电力网络结构及运行条件分析整个网络的潮流分布,计算结果中最基本的电气量是系统各母线的复电压的稳态值,其内容属于电力系统稳态分析。
不确定潮流是针对电力系统内负荷和发电机的出力不能精确知道,但知道其一定包含在某个给定区间范围内的情况,通过分析得到的潮流分布结果是一个区间解。由于方法的局限性,在这个区间解里不仅包含所有符合负荷与发电机区间约束的可能解,还会包含一些实际上不符合所有约束条件的解。
前推回代类潮流分析方法是求解辐射型配电网潮流的有效方法,以其简单、灵活、方便等优点,在配电网潮流计算中获得了广泛的应用。前人将区间算术引入到前推回代分析方法中,提出基于区间算术的辐射型配电网不确定潮流分析方法,但该方法的缺陷在于分析结果往往过于保守,使用区间算术得到的区间解的范围往往会比实际中可能出现的范围大很多,这个问题在数次前推回代的长计算链中表现得尤其突出。
发明内容
技术问题:本发明的目的是提供一种能够得到更精确的电压波动范围并包含所有潮流分布结果的基于仿射算术的辐射型配电网不确定潮流分析方法。
技术方案:本发明的基于仿射算术的辐射型配电网不确定潮流分析方法,包括以下步骤:
1)选取平衡节点为根节点,采用深度优先法则或相关近似法则遍历整个网络,获取网络中的节点个数N并根据网络节点被访问的先后顺序依次为节点编号;
2)根据步骤1)得到的节点个数N,针对每个节点,用一个包含2N个噪声元的仿射量表示节点负荷的波动范围,从而将每个节点的负荷波动范围转换为仿射量;
3)将步骤2)中得到的仿射量及各个节点的额定电压,作为负荷及电压的初值赋给除平衡节点外的各个节点;
4)运用仿射算术的运算法则,采用基于电流的U-I-PARS型前推回代法,从网络的末端开始不断前推回代,直至满足收敛精度的要求,得到各节点电压的波动范围,基于电流的U-I-PARS型前推回代法中,在前推过程中修正电压,在回代过程中只修正电流。
本发明的步骤1)中,若采用相关近似法则遍历整个网络,则必须使节点编号满足:网络中任意两个处于同一条支路上的节点,编号较小的节点位于另一个节点的上游。
本发明的步骤2)中,将任一节点负荷的波动范围转化为一个仿射量的具体方法为:
21)将噪声元εi从1到2N编号,并将步骤1)中得到的节点编号与噪声元序号一一对应,即编号为k的节点的有功P与无功Q对应的噪声元序号分别为2k-1和2k;
22)确定网络中的任一节点负荷波动范围的仿射量表达式为: 其中表示仿射量,x0为不确定量的波动范围的中间值,εi为噪声元,xi为噪声元系数,k为当前节点的编号,N为网络中的节点个数;
23)将网络中的每个节点负荷的波动范围,转化为一个仿射量,其中节点k的具体转化方法为:将负荷波动范围的中间值作为x0的值,其为一个复数,将节点有功P及无功Q偏离各自中间值的最大值分别作为各自相对应的噪声元系数赋给x2k-1及x2k,剩下的与该节点没有对应关系的噪声元系数,均赋初值0。
本发明的步骤3)中,平衡节点的电压为给定电压,在分析过程中保持不变。
本发明中,步骤4)的具体流程为:
41)采用仿射算术的运算法则进行第一次前推回代:
首先进行回代:根据每个节点的电压以及负荷转化后的仿射量,运用公式
ILk=(Pk+jQk)*/Uk
,计算得到每个节点的负荷电流ILk,其中,Pk+jQk为该节点负荷转化后的仿射量,Uk为该节点的电压;
从配电网各条支路的末端开始,根据基尔霍夫电流定律,从网络中各条支路的末端依次向首端回代,直至根节点,从而得到任意相连的两节点间电流的仿射量;
然后进行前推:根据回代过程中得到的每个节点的电压Uk,以及回代过程中求得的各支路电流的仿射量,运用公式
Ub=Ua-IabZab
从配电网的根节点开始向各条支路末端前推,依次求出各个节点的新的电压,其中Ua、Ub分别表示相连的两个节点a和b的电压,Iab表示连接这两个节点的线路电流,Zab表示该线路的阻抗;
42)采用仿射算术的运算法则进行后续的前推回代:
将前一次前推回代得到的各节点的电压赋给各个节点,作为新的节点电压,其为一个仿射量,然后按步骤41)中的流程完成新一次的前推回代计算,其中,由于节点电压不再是一个确定量,而是一个仿射量,因此在采用仿射运算法则,计算得到每个节点的负荷电流ILk时,产生4个新的噪声元及其系数,将其添加到原有仿射量的后面;
43)进行收敛精度判断,如果所得结果满足收敛精度的要求,则将最后一次前推回代得到的各节点电压的仿射量作为最终的分析结果,并将其转化为区间形式的电压波动范围输出,否则,回到步骤42)。
本发明步骤4)中的仿射运算法则具体如下:
对于给定的两个包含复数系数的仿射量
其中,x0,y0为复数,xi,yi(i=1,3,5......2n-1)为实数,xi,yi(i=2,4,6......2n)为虚数;
其运算法则为:
加法:
减法:
乘法:
其中,对任一复数Z=a+jb,real(Z)表示该复数的实部a,imag(Z)表示该复数的虚部jb;
除法:
其中,对任一实数a,|a|表示其绝对值,对任一虚数jb,|jb|表示j·|b|;
本发明的步骤43)中,进行收敛精度判断的具体方法为:
a)将求得的各个节点电压的仿射量转换为区间形式,其具体方法为:对于一个仿射量 将其对应的区间形式表示为
b)将每个节点区间形式的电压与上一次前推回代得到的区间形式的电压相比较,求出每个节点前后两次区间下限的差值和上限的差值,取所有这些差值中的最大值与事先设置的收敛精度相比较,若小于收敛精度,即
则符合收敛精度要求,否则,不符合收敛精度要求,其中表示编号k的节点第m次迭代后的电压仿射量转换为区间形式后的下限,表示编号k的节点第m次迭代后的电压仿射量转换为区间形式后的上限。
有益效果:本发明与现有技术相比,具有以下优点:
前人提出的基于区间算术的辐射型配电网不确定潮流分析方法,由于未考虑网络中各个不确定量之间的关联及依赖关系,经过一个接一个的区间运算后得到的结果区间往往比实际范围大的多,即分析结果过于保守。本发明提供的基于仿射算术的辐射型配电网不确定潮流分析方法,用噪声元εi记录各个不确定量之间的依赖关系,通过仿射算术得到最终的分析结果,从而可以有效地缓解分析结果的保守性,即剔除原有保守性结果中一些不可能存在的解,得到波动范围更加小的分析结果,并使其包含所有的符合约束条件的可能解。
附图说明
图1是本发明的方法流程示意图。
图2是编号后的IEEE33节点系统结构图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体介绍。
附图1是本发明的方法流程示意图,介绍了本发明的方法基本步骤。附图2是编号后的IEEE33节点系统,给定了网络的基准电压大小、平衡节点的电压及相角大小、所有非平衡节点的负荷波动范围(在本例中给定为上下波动±10%),以下以该系统为例说明本发明方法的具体实现。
1)选取平衡节点为根节点,采用深度优先法则或相关近似法则遍历整个网络,获取网络中的节点个数N并根据所述网络节点被访问的先后顺序依次为节点编号。以图2所示网络为例,其为一个按照相关规则遍历及编号完毕的网络,共有33个节点,其中,节点1的电压及相角已知,故其为平衡节点,将其选为根节点。同时,所有的节点编号满足要求:网络中任意两个处于同一条支路上的节点,编号较小的节点位于另一个节点的上游。
2)根据步骤1)得到的节点个数N,针对每个节点,用一个包含2N个噪声元的仿射量表示节点负荷的波动范围,从而将每个节点的负荷波动范围转换为仿射量。由于系统中每个节点负荷的不确定性包括有功功率P和无功功率Q两种波动,所以针对图2所示的33节点系统,N=33,系统中共有66个不确定源,即转化后的仿射形式应有66个噪声元。以网络中节点16为例,其对应的噪声元序号为31和33,其负荷波动范围的中间值标幺化后的值为0.006+j0.002,由于其有功功率及无功功率波动均为±10%,即有功功率波动±0.0006,无功功率波动±0.0002,则该节点负荷波动范围的仿射量可以表示为
3)将步骤2)中得到的仿射量及各个节点的额定电压,作为负荷及电压的初值赋给除平衡节点外的各个节点。所述的额定电压即所给网络的基准电压,如图2所示系统的基准电压为12.66KV,其标幺化后为1。
4)运用仿射算术的运算法则,采用基于电流的U-I-PARS型前推回代法,从网络的末端开始不断前推回代,直至满足收敛精度的要求,得到各节点电压的波动范围。
41)采用仿射算术的运算法则进行第一次前推回代:
回代过程:
如图2所示网络,节点23,24,25的负荷波动范围转换为仿射量后分别为:
第一次回代时,各个节点的电压均为额定电压,标幺化后为1,根据公式
ILk=(Pk+jQk)*/Uk
计算得到节点23,24,25的负荷电流分别为: 根据基尔霍夫电流定律,可以求出节点24,25之间线路上的电流 节点23,24之间线路上的电流 节点22,23之间线路上的电流
依此类推,从网络中每条支路的末端依次向首端回代,直至根节点,从而得到任意相连的两节点间电流的仿射量。
前推过程:
如图2所示系统,以节点23,24为例,第一次前推进行到节点23时,其电压仿射量 其中省略号部分是余下56个噪声元系数与噪声元的乘积。根据公式
Ub=Ua-IabZab
代入节点23,24之间线路的阻抗以及回代过程求出的流过该线路的电流,计算得到节点24的电压仿射量 其中省略号部分是余下56个噪声元系数与噪声元的乘积。
依此类推,从配电网的根节点开始向各条支路末端前推,可以求出各个节点的新的电压仿射量。
42)采用仿射算术的运算法则进行后续的前推回代:
后续的前推回代与第一次前推回代的区别在于,各节点电压的初值不再是额定电压,而是前一次前推回代得到的各节点电压仿射量,因此在采用仿射运算法则,运用公式计算得到每个节点的负荷电流ILk时,将产生4个新的噪声元及其系数。后续前推回代的其他步骤及流程与第一次前推回代相同,此处不再赘述。
43)进行收敛精度判断,如果所得结果满足收敛精度的要求,则将最后一次前推回代得到的各节点电压的仿射量作为最终的分析结果,并将其转化为区间形式的电压波动范围输出,否则,回到步骤42)。仍以图2所示系统为例,为简化说明,假设第二次前推回代后,系统中某个节点的电压仿射量为 而第一次前推回代得到的该节点的电压仿射量为 将电压仿射量转换为区间量,则第二次前推回代后的电压区间[0.89+j0.93,0.91+j0.97],第一次前推回代后的电压区间[0.89+j0.94,0.93+j0.98],下限与下限之差的模值为||-j0.01||=0.01,上限与上限之差的模为||-0.02-j0.01||=0.02236,依此方法,求出全部32个非平衡节点上限之差及下限之差的模值,取其最大值与收敛精度ξ相比较,若大于ξ,则回到步骤42),若小于ξ,则前推回代结束,将最终的区间形式的电压波动范围输出,运算结束。
Claims (5)
1.一种基于仿射算术的辐射型配电网不确定潮流分析方法,其特征在于,该方法包括如下步骤:
1)选取平衡节点为根节点,采用深度优先法则或相关近似法则遍历整个网络,获取网络中的节点个数N并根据所述网络节点被访问的先后顺序依次为节点编号;
2)根据所述步骤1)得到的节点个数N,针对每个节点,用一个包含2N个噪声元的仿射量表示节点负荷的波动范围,从而将每个节点的负荷波动范围转换为仿射量;
3)将所述步骤2)中得到的仿射量及各个节点的额定电压,作为负荷及电压的初值赋给除平衡节点外的各个节点;
4)运用仿射算术的运算法则,采用基于电流的U-I-PARS型前推回代法,从网络的末端开始不断前推回代,直至满足收敛精度的要求,得到各节点电压的波动范围,所述基于电流的U-I-PARS型前推回代法中,在前推过程中修正电压,在回代过程中只修正电流;
所述步骤1)中,若采用相关近似法则遍历整个网络,则必须使节点编号满足:网络中任意两个处于同一条支路上的节点,编号较小的节点位于另一个节点的上游;
所述步骤2)中,将任一节点负荷的波动范围转化为一个仿射量的具体方法为:
21)将噪声元εi从1到2N编号,并将所述步骤1)中得到的节点编号与噪声元序号一一对应,即编号为k的节点的有功P与无功Q对应的噪声元序号分别为2k-1和2k;
22)确定网络中的任一节点负荷波动范围的仿射量表达式为: 其中表示仿射量,x0为不确定量的波动范围的中间值,εi为噪声元,xi为噪声元系数,k为当前节点的编号,N为网络中的节点个数;
23)将网络中的每个节点负荷的波动范围,转化为一个仿射量,其中节点k的具体转化方法为:将负荷波动范围的中间值作为x0的值,其为一个复数,将节点有功P及无功Q偏离各自中间值的最大值分别作为各自相对应的噪声元系数赋给x2k-1及x2k,剩下的与该节点没有对应关系的噪声元系数,均赋初值0。
2.根据权利要求1所述的基于仿射算术的辐射型配电网不确定潮流分析方法,其特征在于,所述步骤3)中,平衡节点的电压为给定电压,在分析过程中保持不变。
3.根据权利要求1所述的基于仿射算术的辐射型配电网不确定潮流分析方法,其特征在于,所述步骤4)的具体流程为:
41)采用仿射算术的运算法则进行第一次前推回代:
首先进行回代:根据每个节点的电压以及负荷转化后的仿射量,运用公式
ILk=(Pk+jQk)*/Uk
计算得到每个节点的负荷电流ILk,其中,Pk+jQk为该节点负荷转化后的仿射量,Uk为该节点的电压,*表示共轭复数,即对任一复数Z,Z*表示该复数的共轭复数;
从配电网各条支路的末端开始,根据基尔霍夫电流定律,从网络中各条支路的末端依次向首端回代,直至根节点,从而得到任意相连的两节点间电流的仿射量;
然后进行前推:根据回代过程中得到的每个节点的电压Uk,以及回代过程中求得的各支路电流的仿射量,运用公式
Ub=Ua-IabZab
从配电网的根节点开始向各条支路末端前推,依次求出各个节点的新的电压,其中Ua、Ub分别表示相连的两个节点a和b的电压,Iab表示连接这两个节点的线路电流,Zab表示该线路的阻抗;
42)采用仿射算术的运算法则进行后续的前推回代:
将前一次前推回代得到的各节点的电压赋给各个节点,作为新的节点电压,其为一个仿射量,然后按所述步骤41)中的流程完成新一次的前推回代计算,其中,由于节点电压不再是一个确定量,而是一个仿射量,因此在采用仿射运算法则,计算得到每个节点的负荷电流ILk时,产生4个新的噪声元及其系数,将其添加到原有仿射量的后面;
43)进行收敛精度判断,如果所得结果满足收敛精度的要求,则将最后一次前推回代得到的各节点电压的仿射量作为最终的分析结果,并将其转化为区间形式的电压波动范围输出,否则,回到步骤42)。
4.根据权利要求3所述的基于仿射算术的辐射型配电网不确定潮流分析方法,
其特征在于,所述步骤4)中的仿射运算法则具体如下:
对于给定的两个包含复数系数的仿射量
其中,x0,y0为复数,xi,yi(i=1,3,5......2n-1)为实数,xi,yi(i=2,4,6......2n)为虚数;
其运算法则为:
加法:
减法:
乘法:
其中,对任一复数Z=a+jb,real(Z)表示该复数的实部a,imag(Z)表示该复数的虚部jb;
除法:
其中,对任一实数a,|a|表示其绝对值,对任一虚数jb,|jb|表示j·|b|。
5.根据权利要求3所述的基于仿射算术的辐射型配电网不确定潮流分析方法,其特征在于,所述步骤43)中,进行收敛精度判断的具体方法为:
a)将求得的各个节点电压的仿射量转换为区间形式,其具体方法为:对于一个仿射量 将其对应的区间形式表示为
b)将每个节点区间形式的电压与上一次前推回代得到的区间形式的电压相比较,求出每个节点前后两次区间下限的差值和上限的差值,取所有这些差值中的最大值与事先设置的收敛精度相比较,若小于收敛精度,即
则符合收敛精度要求,否则,不符合收敛精度要求,其中表示编号k的节点第m次迭代后的电压仿射量转换为区间形式后的下限,表示编号k的节点第m次迭代后的电压仿射量转换为区间形式后的上限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210558626.5A CN103093079B (zh) | 2012-12-20 | 2012-12-20 | 基于仿射算术的辐射型配电网不确定潮流分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210558626.5A CN103093079B (zh) | 2012-12-20 | 2012-12-20 | 基于仿射算术的辐射型配电网不确定潮流分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103093079A CN103093079A (zh) | 2013-05-08 |
CN103093079B true CN103093079B (zh) | 2016-03-16 |
Family
ID=48205638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210558626.5A Expired - Fee Related CN103093079B (zh) | 2012-12-20 | 2012-12-20 | 基于仿射算术的辐射型配电网不确定潮流分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103093079B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104104081B (zh) * | 2014-07-30 | 2016-03-16 | 东南大学 | 一种基于优化方法的非迭代不确定潮流分析方法 |
CN104795811B (zh) * | 2015-03-17 | 2017-07-18 | 河海大学 | 一种电力系统区间状态估计方法 |
CN107482633B (zh) * | 2017-08-22 | 2020-03-31 | 东南大学 | 一种适用于辐射状配电网的非迭代区间潮流算法 |
CN111181166B (zh) * | 2020-02-17 | 2022-10-28 | 福州大学 | 一种预测校正的不确定性仿射潮流方法 |
CN113300368B (zh) * | 2021-06-16 | 2022-07-05 | 东南大学 | 一种计及车网互动的配网不确定潮流分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615213A (zh) * | 2009-07-21 | 2009-12-30 | 清华大学 | 基于扩展不确定度的电力系统状态估计结果评价方法 |
CN101752903A (zh) * | 2009-11-27 | 2010-06-23 | 清华大学 | 一种时序递进式电力调度方法 |
-
2012
- 2012-12-20 CN CN201210558626.5A patent/CN103093079B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615213A (zh) * | 2009-07-21 | 2009-12-30 | 清华大学 | 基于扩展不确定度的电力系统状态估计结果评价方法 |
CN101752903A (zh) * | 2009-11-27 | 2010-06-23 | 清华大学 | 一种时序递进式电力调度方法 |
Non-Patent Citations (4)
Title |
---|
An affine arithmetic-based methodology for reliable power flow analysis in the presence of data uncertainty;Alfredo Vaccaro等;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20100531;第25卷(第2期);624-632页 * |
一种改进的辐射状配电网潮流计算方法;刘蔚等;《电力系统保护与控制》;20090801;第37卷(第15期);1-4页 * |
基于区间和仿射运算的不确定潮流算法;丁涛等;《电力系统自动化》;20120710;第36卷(第13期);51-55页 * |
复区间潮流保守性问题的解决方案;王守相等;《电力系统自动化》;20051010;第29卷(第19期);25-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103093079A (zh) | 2013-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103093079B (zh) | 基于仿射算术的辐射型配电网不确定潮流分析方法 | |
CN101635457B (zh) | 一种基于状态估计残差参数灵敏度的电网参数估计方法 | |
CN103323688B (zh) | 一种电力系统的谐波状态估计方法及设备 | |
CN103236691B (zh) | 基于复仿射数学理论的三相不平衡潮流计算方法 | |
CN101651343B (zh) | 基于混合动态仿真的电力系统模型校核方法及系统 | |
CN101964525B (zh) | 一种支持大规模电流量测的配电网状态估计方法 | |
CN103646152A (zh) | 一种基于矩阵指数的电力系统电磁暂态仿真方法 | |
CN102590685B (zh) | 一种配电网电流匹配状态估计方法 | |
CN109494724B (zh) | 基于lu分解的大电网戴维南等值参数在线辨识方法 | |
CN103106256A (zh) | 一种基于Markov链的正交化插值的GM(1,1)模型预测方法 | |
CN107681685A (zh) | 一种考虑光伏非线性相关性的电力系统概率潮流计算方法 | |
CN105048468A (zh) | 基于分布式计算的输配电网一体化电压稳定评估方法 | |
CN101976837A (zh) | 一种电力系统潮流无解调整方法 | |
CN103559042A (zh) | 一种建立基于iec61970/61968的配电网线损cim模型的方法 | |
CN104113061A (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN104882884A (zh) | 基于马尔可夫链-蒙特卡罗法的系统谐波概率评估方法 | |
CN107482633A (zh) | 一种适用于辐射状配电网的非迭代区间潮流算法 | |
CN103474989A (zh) | 一种基于灵敏度分析的网络重构方法 | |
CN113937764A (zh) | 一种低压配电网高频量测数据处理与拓扑识别方法 | |
CN103824234A (zh) | 一种基于分块和层级结构的配电系统可靠性评估方法 | |
CN107257130A (zh) | 基于区域量测解耦的低压配电网损耗计算方法 | |
CN106921161A (zh) | 电力系统薄弱点识别方法及薄弱点电压分布特性分析方法 | |
CN104537233B (zh) | 一种基于核密度估计的配电网伪量测生成方法 | |
CN111654392A (zh) | 基于互信息的低压配电网拓扑识别方法及系统 | |
CN107332239B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160316 Termination date: 20201220 |
|
CF01 | Termination of patent right due to non-payment of annual fee |