CN105004926A - 一种交流电相位频率幅值跟踪重构的方法 - Google Patents
一种交流电相位频率幅值跟踪重构的方法 Download PDFInfo
- Publication number
- CN105004926A CN105004926A CN201510448078.4A CN201510448078A CN105004926A CN 105004926 A CN105004926 A CN 105004926A CN 201510448078 A CN201510448078 A CN 201510448078A CN 105004926 A CN105004926 A CN 105004926A
- Authority
- CN
- China
- Prior art keywords
- phase
- alternating current
- amplitude
- frequency
- phase alternating
- 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
- Particle Accelerators (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明公开了一种交流电相位频率幅值跟踪重构方法,按照该重构方法,对采集的交流电实际值分析出交流电的频率、相位、相序关系,其后利用二阶广义积分器、频率和相位信息得到各相实时的估计幅值,再次估计幅值代入提出的算法,减小估计的频率和相位误差。按照本发明实现的交流电相位频率幅值跟踪重构的方法,无需提供前馈参考较频率,能够有效分析任何位置角频率的交流电相位、频率和幅值,并且对于单相或三相电压、电流信号均可实现分析,而且能够实现有效、快速的实时重构。
Description
技术领域
本发明属于交流电相位检测与重构的技术领域,更具体地,涉及一种交流电相位频率幅值跟踪重构的方法。
背景技术
光伏、风力等新能源并网发电快速发展,分布式发电成为热点研究对象。软件锁相环是分布式发电系统中的关键技术,用以检测污染的电网电压的频率和相位,为并网发电提供支持。受污染的电网不仅电压不平衡,而且含有多次谐波分量。为越来越多的学者对其进行了研究,提出了诸多锁相环的控制与设计方案。例如基于低通滤波锁相环、基于加权最小二乘法估计的锁相环、基于扩展卡尔曼滤波器的锁相环、基于空间矢量滤波器的锁相环等等。他们都存在以下问题:锁相的精度不高,响应慢,对系统频率变化和三相电压不平衡较敏感。随着研究的深入,许多锁相环被提出,从而使锁相环的性能得到不断的改善和提升。基于单同步坐标系软件锁相环,当三相电网电压不平衡时,得到的电压正序分量的相位存在较大2次谐波。基于对称分量法的单同步坐标系软件锁相环,对频率适应性相对不足。为了提高锁相环的锁相精度和响应,提出了基于双同步坐标系的解耦软件锁相环。该方法能在电压平衡和不平衡条件下准确地获取电网电压的相位、幅值和频率信息,具有较好的电网适应性。
由于单相电网电压信号不存在构成同步旋转坐标系的两相正交电压信号,上述基于同步旋转坐标系的三相电网锁相环控制策略一般不适用于单相电网锁相环的控制。单相电网电压的锁相环方案有:基于延迟法构造虚拟两相的单相锁相环,存在周期性的延迟,动态响应慢。基于P-Q理论的单相锁相环,在电压跌落、相位突变时,响应慢。基于Park反变换的单相锁相环,滤波器设计困难,电压不平衡时响应慢。
基于双二阶广义积分器的检测算法在许多方面都用应用。基于双二阶广义积分器的三相软件锁相环能在电压平衡和不平衡条件下分析三相电网电压的相位、幅值和频率信息,响应快。而且基于双二阶广义积分器的单相软件锁相环也能准确快速的获取单相交流电的信息。但是,一方面,以上锁相环技术均需要提供被测对象的基频和三相电压的相序;另一方面,当交流电谐波分量较大时,实际的相位曲线不是标准的正弦曲线,角频率不为恒定值,而以上算法得到的角频率均为恒定值,这样相位误差就会变大,在变频调速、未知电网和锁相环相位跟踪等复杂的应用情况下,上述算法并不能实际解决问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于二阶广义积分器交流电相位频率跟踪重构的算法,其目的在于无需给定参考频率,可以对任意未知幅值、频率的交流电分析,从而得出各自的幅值、频率和相位,具备良好的跟踪和重构效果。
为实现上述目的,按照本发明的一个方面,提供了一种交流电相位频率幅值跟踪重构方法,其特征在于,该方法包括如下步骤:
(1)实时采集交流电,判断其是单相或是三相交流电,若是单相,则进入步骤(2),若是三相交流电,则进入步骤(3);
(2)实现单相交流电的重构;其包括如下子步骤:
(2-1)连续采集N组单相交流电实时数据,并将其输入二阶广义积分器求解出两个正交的输出量,利用上述两个正交的输出量求出实时交流电幅值;
(2-2)求解其中θ为单相交流电实际相位角,为单相交流电跟踪相位角,其中其中X、Y是分别与所述两个正交的输出量相关的系数;
(2-3)将所述步骤(2-2)获得的作为PI控制器的输入,进行调节,进而输出所述单相交流电的角频率,所述角频率积分之后获得所述单相交流电的相位角;
(2-4)由上述步骤获得的所述幅值、角频率和相位角,实现所述单相交流电的重构;
(3)实现三相交流电的重构;其包括如下子步骤:
(3-1)连续采集N组三相交流电实时数据,计算得出其中θ为三相交流电实际相位角,为三相交流电跟踪相位角,其中M、N是与三相电压幅值和三相电压信号相关的表达式;
(3-2)将所述步骤(3-1)获得的作为PI控制器的输入,进行调节,进而输出所述三相交流电的角频率,所述角频率积分之后获得所述三相交流电的相位角;
(3-3)将所述步骤(3-2)的输出输入二阶广义积分器,求解得出三组两个正交的输出量,计算得出各相实时幅值;
(3-4)由所述步骤(3-2)、(3-3)获得的所述幅值、角频率和相位角,实现所述单相交流电的重构。
进一步地,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为其中t和l分别为求出的所述两正交的输出量。
进一步地,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为反Park变换法,即
总体而言,通过本发明所构思的以上技术方案与现有技术相比,由于对采集的三相交流电实际值分析出三相交流电的频率、相位、相序关系,其次利用二阶广义积分器、频率和相位信息得到各相实时的估计幅值,再次估计幅值代入提出的算法,减小估计的频率和相位误差,能够取得下列有益效果:
(1)巧妙地应用了正弦差公式及二阶广义积分器,无需提供前馈参考角频率,能有效分析任意未知角频率的交流电的相位、频率和幅值;
(2)对于单相或三相电压、电流信号,均可实现信号分析;
(3)本发明能适应含有谐波和噪声的交流电信号,并分析出准确的交流电的相位、频率和幅值;
(4)本发明实用有效、计算量小,在0.2s内即计算出交流电的频率、相位和幅值,并且实现交流电的实时重构。
附图说明
图1是按照本发明实现的针对三相交流电锁相环原理框图;
图2是按照本发明实现的针对单相交流电锁相环原理框图;
图3是按照本发明实现的交流电跟踪重构方法流程图;
图4是发电机实际控制绕组三相电流的观测图;
图5是按照本发明实现的三相锁相环计算得出的相位和频率的观测图;
图6是按照本发明的方法重构的B相电流与实际滤波后的三相电流的观测图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
将三相交流电压用正弦函数表示,同理也可以用余弦函数分析求解。三相交流电压表达式如下所示:
va+vb+vc=0 (4)
A、B、C分别表示三相交流电压幅值,当三相电压幅值相等时,满足A=B=C;当三相交流电压不平衡,幅值不相等时。忽略电压幅值恒为零情况,式(1)-(3)也可表示为:
由(5)-(7)可得:
其中为a相电压实际相位角;为a相电压算法跟踪相位角和角频率。当未知三相交流电的相序时,可使用正序锁相环算法。如果计算出来的角频率为负,则三相交流电为负序;如果计算出来的角频率为正,则三相交流电为正序。
以a相电压作为输入信号为例,单相交流电的锁相算法等式为(13)-(14)。系统达到稳态时,输入信号v和输出信号v′、qv′角频率相等幅值相等|v|=|v′|=|qv′|,v′与qv′正交,v′相位角与输入信号v一致且超前qv′相位角90°。
利用改进后的二阶广义积分器,有两种方法求输入信号的幅值。以a相电压作为输入信号为例,第一种是直接求解方法式(15)所示。
第二种是反Park变换方法求解,如式(16)所示。
将作为一个PI控制器的输入,进行调节。PI控制器的输出为交流电的角频率,角频率积分之后得到交流电的相位角。
综上所诉,可以得出的三相交流电锁相环原理框图,如图1所示;单相交流电锁相环原理框图,如图2所示。
基于上述分析,如图3所示的流程图,按照本发明实施例的交流电频率、相位和幅值跟踪与重构方法包括如下步骤:
(1)根据交流电的实际状况,判断交流电是否单相,是则顺序执行如下步骤(A1)~(A4);否则说明交流电是三相交流电,顺序执行如下步骤(B1)~(B4)。
步骤(A1)~(A4)分别为:
(A1)连续采集N组单相交流电实时数据D1~DN,其中,第k组电机运行数据Dk,利用二阶广义积分器求解得出两个正交的输出量。求解利用式(15)或式(16)计算得出实时交流电幅值。
(A2)利用式(13)和式(14)计算得出
(A3)作为一个PI控制器的输入,进行调节。PI控制器的输出为交流电的角频率,角频率积分之后得到交流电的相位角。
(A4)利用求解得出的交流电幅值,角频率和相位角,实现交流电的重构。
步骤(B1)~(B4)分别为:
(B1)连续采集N组三相交流电实时数据D1~DN,其中,第k组电机运行数据Dk=(va,vb,vc),利用式(9)到式(12)计算得出
(B2)作为一个PI控制器的输入,进行调节。PI控制器的输出为交流电的角频率,如果计算得出的角频率为负数,则三相交流电的相许未负;反之为正。角频率积分之后得到交流电的相位角。
(B3)利用二阶广义积分器,将三相交流电求解得出三组两个正交的输出量。再利用式(15)或式(16)计算得出交流电各相实时幅值。
(B4)利用求解得出的三相交流电各相幅值,角频率和相位角,实现三相交流电的重构。
为使本领域技术人员更好地理解本发明,下面结合具体实例对本发明的交流电频率、相位和幅值跟踪与重构方法进行详细说明。
利用发电实验平台,使用提出的方法跟踪和重构发电机控制绕组三相电流。发电机控制绕组三相电流作为研究对象,发电机的同步转速为750RPM。发电机的转速从670RPM匀加速到860RPM,控制绕组的给定电流频率从负序的5.33Hz匀减为0,再从正序的0匀加速到7.33Hz,采样周期为0.1ms,计算时间4s。
采样的三相电流如图4所示,A相(1通道所示),B相(2通道所示),C相(3通道所示)。实际控制绕组三相电流是变频器直接输出的,不仅含有测量噪声,而且含有高次谐波分量。滤波后的三相电流作为提出的算法的输入量。经过上述方法,利用式(9)到式(12)计算得出PI控制器的输入量。进行调节计算得出三相电流的角频率,角频率积分之后得到三相电流的相位角,如图5所示,三相锁相环计算得出的相位(3通道所示)、未滤波频率(2通道所示)和滤波后的频率(4通道所示)。利用二阶广义积分器和式(15)或式(16)计算得出交流电各相实时幅值。利用求解得出的三相交流电各相幅值,角频率和相位角,实现三相交流电的重构,如图6所示,三相锁相环算法重构的B相电流(4通道所示)与实际滤波后的三相电流,A相(1通道所示),B相(2通道所示),C相(3通道所示)。按照本发明实现的重构方法的三相锁相环运行0.2s后,可实时跟踪电网频率、相位和电压幅值。未滤波的频率波动在±0.3Hz,而滤波后的频率波动在±0.05Hz,相位跟踪误差很小。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种交流电相位频率幅值跟踪重构方法,其特征在于,该方法包括如下步骤:
(1)实时采集交流电,判断其是单相或是三相交流电,若是单相,则进入步骤(2),若是三相交流电,则进入步骤(3);
(2)实现单相交流电的重构;其包括如下子步骤:
(2-1)连续采集N组单相交流电实时数据,并将其输入二阶广义积分器求解出两个正交的输出量,利用上述两个正交的输出量求出实时交流电幅值;
(2-2)求解其中θ为单相交流电实际相位角,为单相交流电跟踪相位角,其中其中X、Y是分别与所述两个正交的输出量相关的系数;
(2-3)将所述步骤(2-2)获得的作为PI控制器的输入,进行调节,进而输出所述单相交流电的角频率,所述角频率积分之后获得所述单相交流电的相位角;
(2-4)由上述步骤获得的所述幅值、角频率和相位角,实现所述单相交流电的重构;
(3)实现三相交流电的重构;其包括如下子步骤:
(3-1)连续采集N组三相交流电实时数据,计算得出其中θ为三相交流电实际相位角,为三相交流电跟踪相位角,其中M、N是与三相电压幅值和三相电压信号相关的表达式;
(3-2)将所述步骤(3-1)获得的作为PI控制器的输入,进行调节,进而输出所述三相交流电的角频率,所述角频率积分之后获得所述三相交流电的相位角;
(3-3)将所述步骤(3-2)的输出输入二阶广义积分器,求解得出三组两个正交的输出量,计算得出各相实时幅值;
(3-4)由所述步骤(3-2)、(3-3)获得的所述幅值、角频率和相位角,实现所述单相交流电的重构。
2.如权利要求1所述的交流电相位频率幅值跟踪重构方法,其特征在于,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为其中t和l分别为求出的所述两正交的输出量。
3.如权利要求1所述的交流电相位频率幅值跟踪重构方法,其特征在于,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为反Park变换法,即
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510448078.4A CN105004926B (zh) | 2015-07-27 | 2015-07-27 | 一种交流电相位频率幅值跟踪重构的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510448078.4A CN105004926B (zh) | 2015-07-27 | 2015-07-27 | 一种交流电相位频率幅值跟踪重构的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105004926A true CN105004926A (zh) | 2015-10-28 |
CN105004926B CN105004926B (zh) | 2018-01-26 |
Family
ID=54377669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510448078.4A Expired - Fee Related CN105004926B (zh) | 2015-07-27 | 2015-07-27 | 一种交流电相位频率幅值跟踪重构的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105004926B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108226588A (zh) * | 2017-11-01 | 2018-06-29 | 中国矿业大学(北京) | 一种适用于单相和三相电力系统的软件锁相方法 |
CN108306569A (zh) * | 2018-02-02 | 2018-07-20 | 华中科技大学 | 基于广义积分器的永磁同步电机无速度控制方法及系统 |
CN110635484A (zh) * | 2019-11-13 | 2019-12-31 | 平顶山学院 | 一种单相有源滤波器指令电流的生成方法、装置及设备 |
CN110907030A (zh) * | 2019-10-14 | 2020-03-24 | 中国电力科学研究院有限公司 | 一种用于测量三相电容器装置噪声的方法及系统 |
CN112325933A (zh) * | 2020-10-29 | 2021-02-05 | 贵州电网有限责任公司 | 一种基于边缘计算的多能源系统采集及预处理方法和终端 |
CN113224993A (zh) * | 2021-05-25 | 2021-08-06 | 中国神华能源股份有限公司哈尔乌素露天煤矿 | 矿用电动轮重启方法和系统 |
CN113447715A (zh) * | 2021-07-30 | 2021-09-28 | 燕山大学 | 一种2次谐波式电网同步锁频环方法 |
CN115795253A (zh) * | 2021-09-09 | 2023-03-14 | 武汉市聚芯微电子有限责任公司 | 一种频率跟踪方法、装置、设备以及计算机存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102185319A (zh) * | 2011-05-10 | 2011-09-14 | 重庆新世杰电气股份有限公司 | 一种有源电力滤波器的谐波检测锁相方法和装置 |
JP5255930B2 (ja) * | 2008-07-04 | 2013-08-07 | 新電元工業株式会社 | 位相検出器 |
CN103986458A (zh) * | 2014-04-16 | 2014-08-13 | 国家电网公司 | 一种基于重复控制的微电网单相并网锁相环控制方法 |
CN104158540A (zh) * | 2014-08-27 | 2014-11-19 | 河北大学 | 一种三相数字锁相环及锁相方法 |
CN104410408A (zh) * | 2014-12-25 | 2015-03-11 | 国家电网公司 | 一种单相锁相方法及装置 |
CN104600696A (zh) * | 2015-01-05 | 2015-05-06 | 上海电力学院 | 一种用于不平衡条件下的软件锁相环实现方法 |
CN104786873A (zh) * | 2015-03-25 | 2015-07-22 | 西南交通大学 | 一种基于锁相环的电压重构方法 |
-
2015
- 2015-07-27 CN CN201510448078.4A patent/CN105004926B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5255930B2 (ja) * | 2008-07-04 | 2013-08-07 | 新電元工業株式会社 | 位相検出器 |
CN102185319A (zh) * | 2011-05-10 | 2011-09-14 | 重庆新世杰电气股份有限公司 | 一种有源电力滤波器的谐波检测锁相方法和装置 |
CN103986458A (zh) * | 2014-04-16 | 2014-08-13 | 国家电网公司 | 一种基于重复控制的微电网单相并网锁相环控制方法 |
CN104158540A (zh) * | 2014-08-27 | 2014-11-19 | 河北大学 | 一种三相数字锁相环及锁相方法 |
CN104410408A (zh) * | 2014-12-25 | 2015-03-11 | 国家电网公司 | 一种单相锁相方法及装置 |
CN104600696A (zh) * | 2015-01-05 | 2015-05-06 | 上海电力学院 | 一种用于不平衡条件下的软件锁相环实现方法 |
CN104786873A (zh) * | 2015-03-25 | 2015-07-22 | 西南交通大学 | 一种基于锁相环的电压重构方法 |
Non-Patent Citations (1)
Title |
---|
王鹿军: "分布式发电中三相三电平并网逆变器的若干关键技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108226588A (zh) * | 2017-11-01 | 2018-06-29 | 中国矿业大学(北京) | 一种适用于单相和三相电力系统的软件锁相方法 |
CN108306569A (zh) * | 2018-02-02 | 2018-07-20 | 华中科技大学 | 基于广义积分器的永磁同步电机无速度控制方法及系统 |
CN110907030A (zh) * | 2019-10-14 | 2020-03-24 | 中国电力科学研究院有限公司 | 一种用于测量三相电容器装置噪声的方法及系统 |
CN110907030B (zh) * | 2019-10-14 | 2022-10-04 | 中国电力科学研究院有限公司 | 一种用于测量三相电容器装置噪声的方法及系统 |
CN110635484A (zh) * | 2019-11-13 | 2019-12-31 | 平顶山学院 | 一种单相有源滤波器指令电流的生成方法、装置及设备 |
CN112325933A (zh) * | 2020-10-29 | 2021-02-05 | 贵州电网有限责任公司 | 一种基于边缘计算的多能源系统采集及预处理方法和终端 |
CN113224993A (zh) * | 2021-05-25 | 2021-08-06 | 中国神华能源股份有限公司哈尔乌素露天煤矿 | 矿用电动轮重启方法和系统 |
CN113447715A (zh) * | 2021-07-30 | 2021-09-28 | 燕山大学 | 一种2次谐波式电网同步锁频环方法 |
CN113447715B (zh) * | 2021-07-30 | 2022-04-15 | 燕山大学 | 一种2次谐波式电网同步锁频环方法 |
CN115795253A (zh) * | 2021-09-09 | 2023-03-14 | 武汉市聚芯微电子有限责任公司 | 一种频率跟踪方法、装置、设备以及计算机存储介质 |
CN115795253B (zh) * | 2021-09-09 | 2023-10-03 | 武汉市聚芯微电子有限责任公司 | 一种频率跟踪方法、装置、设备以及计算机存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105004926B (zh) | 2018-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105004926A (zh) | 一种交流电相位频率幅值跟踪重构的方法 | |
CN106501574B (zh) | 一种有源电力滤波器谐波电流检测方法 | |
Qasim et al. | Optimal current harmonic extractor based on unified ADALINEs for shunt active power filters | |
Rani et al. | A three phase PLL with a dynamic feed forward frequency estimator for synchronization of grid connected converters under wide frequency variations | |
CN107706929B (zh) | 基于最小方差滤波的自适应锁相环方法及系统 | |
Yazdani et al. | A real-time extraction of harmonic and reactive current in a nonlinear load for grid-connected converters | |
CN102401858A (zh) | 一种电网电压基波分量及谐波分量的检测方法 | |
CN110557118B (zh) | 一种锁相装置及锁相方法 | |
CN104184148A (zh) | 一种同步旋转坐标系中谐波电流分次控制方法 | |
CN106410858A (zh) | 一种基于双dq坐标变换的软件数字锁相方法 | |
CN105425011A (zh) | 一种适用于单相畸变电网的非线性幅相检测方法 | |
CN109066678A (zh) | 一种改进的apf系统谐波电流检测方法 | |
CN104143837B (zh) | 具有参数自适应特性的逆变器无交流电压传感器控制方法 | |
CN102868163A (zh) | 一种自适应并网变流器三相软锁相环 | |
CN106787650B (zh) | 电网不平衡下有功二次脉动抑制的直接电流指令计算方法 | |
Nwobu et al. | Grid voltage synchronization for unbalanced voltages using the energy operator | |
Mekhilef et al. | Single-phase hybrid active power filter with adaptive notch filter for harmonic current estimation | |
CN113804949A (zh) | 一种适用于宽变频交流发电系统的幅值确定方法及装置 | |
Shi et al. | Adaptive quadrant filter based phase locked loop system | |
CN111200287B (zh) | 一种指定谐波电流注入装置及其谐波电流的给定方法 | |
Lakshmanan et al. | Improved Phase Lock Loop system for grid connected renewable energy sources under distorted conditions | |
Muñoz-Aguilar et al. | A sensor-less sliding mode control scheme for a stand-alone wound rotor synchronous generator under unbalanced load conditions | |
Guo et al. | Improved park inverse transform algorithm for positive and negative sequence separation of grid voltage under unbalanced grid conditions | |
Popa et al. | Improving the Shunt Active Power Filter Control Methods under Distorted and Unbalanced Grid Voltages | |
CN103683320A (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: 20180126 Termination date: 20200727 |
|
CF01 | Termination of patent right due to non-payment of annual fee |