CN104135033B - 新型光伏并网逆变器电压型控制方法 - Google Patents

新型光伏并网逆变器电压型控制方法 Download PDF

Info

Publication number
CN104135033B
CN104135033B CN201410421837.3A CN201410421837A CN104135033B CN 104135033 B CN104135033 B CN 104135033B CN 201410421837 A CN201410421837 A CN 201410421837A CN 104135033 B CN104135033 B CN 104135033B
Authority
CN
China
Prior art keywords
inverter
voltage
grid
control method
power
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
Application number
CN201410421837.3A
Other languages
English (en)
Other versions
CN104135033A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201410421837.3A priority Critical patent/CN104135033B/zh
Publication of CN104135033A publication Critical patent/CN104135033A/zh
Application granted granted Critical
Publication of CN104135033B publication Critical patent/CN104135033B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers

Landscapes

  • Inverter Devices (AREA)
  • Control Of Electrical Variables (AREA)

Abstract

新型光伏并网逆变器电压型控制方法,涉及一种并网逆变器电压型控制方法。为了解决现有电压控制型逆变器并网功率易受电网波动影响,无法直流母线电压稳定的问题。本方法包括如下步骤:步骤一:将闭环反馈控制方法引入到下垂控制中,增加功率环,经拉普拉斯变换后,获得下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0,步骤二:将下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0分别叠加到新型光伏并网逆变器电压型的并网下垂方程,并简化后获得方程,并将简化后获得方程中的功率环替换为直流母线电压环,完成光伏并网逆变器对直流母线电压的控制。本发明主要应用在并网逆变器上。

Description

新型光伏并网逆变器电压型控制方法
技术领域
本发明涉及一种并网逆变器电压型控制方法。
背景技术
并网有功功率通过逆变器输出电压幅值调节,并网无功功率通过逆变器相位来调节,因此,电压控制型逆变器并网时的输出功率方程与逆变器并联下垂控制类似。传统下垂控制忽略了线路阻抗中的阻性成分,认为传输线路以感性为主。而低压电力传输线以阻性为主,因此,可以得到逆变器的电压型并网下垂方程
U=U0-kp(P-P0)(1)
f=f0+kq(Q-Q0)(2)
式中,U为逆变器输出电压幅值参考值,f为逆变器输出频率参考值,f0为逆变器额定输出频率,U0为逆变器额定输出电压幅值,kp为逆变器输出有功功率的下垂系数,kq为逆变器输出无功功率的下垂系数,P为逆变器输出的有功功率,Q为逆变器输出的无功功率,P0为逆变器输出的额定有功功率,Q0为逆变器输出的额定无功功率。
当下垂方程不变时,无功功率与电网频率呈线性关系,若电网频率发生波动,无功功率将随之发生偏移。另外,因为传输线路阻抗上存在压降,并联系统的电压幅值不可能是处处一致的。逆变器输出端电压不会与电网强制相等,而是与电网电压、阻抗大小等多种因素有关。由于阻抗压降的作用,即使电网为额定值,输出的有功功率也不等于所期望的,而电网电压在额定点附近的波动会进一步导致有功功率的偏移。因此,要想克服电网频率波动的影响(保持输出无功功率稳定)并希望有功功率能稳定控制在期望值,则必须改变下垂控制方程。
在两级式光伏并网逆变器中,后级逆变器下垂控制的P0应始终等于Boost变换器输出的光伏阵列最大功率PPV,而Q0始终为0。然而光伏阵列的输出功率PPV一直在变化,若后级逆变器输出功率不变,会造成前后级输出功率的不平衡,引起系统的不稳定,触发直流母线电压的过压/欠压保护。因此,此时应平移P-U下垂曲线,使前、后两级输出功率保持平衡。
发明内容
本发明为了解决现有电压控制型逆变器并网功率易受电网波动影响,无法直流母线电压稳定的问题,本发明提供了一种新型光伏并网逆变器电压型控制方法。
新型光伏并网逆变器电压型控制方法,它包括下垂控制方法,所述的下垂控制方法包括如下步骤:
步骤一:将闭环反馈控制方法引入到下垂控制中,增加功率环,经拉普拉斯变换后,获得下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0
步骤二:将下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0分别叠加到新型光伏并网逆变器电压型的并网下垂方程,并简化后获得方程,
U = U 0 + ( k pp + k ip s ) ( P PV - P ) - - - ( 7 ) ,
f = f 0 + ( k pq + k iq s ) Q - - - ( 8 ) ,
将公式(7)和(8)中的功率环替换为直流母线电压环,获得
U = U 0 - ( k pp + k ip s ) ( U DCref - U DC ) - - - ( 9 ) ,
f = f 0 + ( k pq + k iq s ) Q - - - ( 8 ) ,
完成光伏并网逆变器对直流母线电压的控制;
其中,UDC为直流母线电压,UDCref为直流母线电压给定值,U为逆变器输出电压幅值参考值,f为逆变器输出频率参考值,U0为逆变器额定输出电压幅值,P为逆变器输出的有功功率,Q为逆变器输出的无功功率,kpp和kip分别为P-U下垂控制环PI调节器的比例、积分系数,kpq和kiq分别为Q-f下垂控制环PI调节器的比例、积分系数,s为复变量,PPV为光伏阵列最大功率。
原理分析:
通过改进的下垂功率控制环,逆变器输出的有功功率P,最终能够稳定在给定值PPV,而与电网电压和线路阻抗无关;逆变器输出的无功功率Q能稳定在给定值0,进一步去除了电网频率的影响。
两级式光伏逆变器存在维持直流母线电压稳定的问题,并网运行时,升压变换器用来实现MPPT控制,所以稳定直流母线电压需要通过后级逆变器来完成。通常,下垂方程仅有功率环,不具有控制直流母线电压的功能,因此本发明提出将直流母线电压环加入下垂控制中,以满足两级式光伏逆变器的控制要求。
本发明提出一种新型的光伏逆变器电压型并网控制方法,该方法解决了电压控制型逆变器并网功率易受电网波动影响问题,增强了光伏逆变器对功率的控制能力,满足了始终以最大有功功率并网和维持直流母线电压稳定的要求。
附图说明
图1为下垂控制方法中公式(7)所对应的改进的P-U控制环的原理框图;其中,E表示电网电压有效值,Z表示并网阻抗,δ表示电网电压相角,fgrid表示电网电压频率;
图2为下垂控制方法中公式(8)所对应的改进的Q-f控制环的原理框图;
图3为现有技术中光伏逆变器电压型并网控制方法的原理示意图;
图4为光伏功率波动实验波形;
图5为电网电压波动时改进下垂控制实验波形中,电网电压升高时暂态波形;
图6为电网电压波动时改进下垂控制实验波形中,电网电压降低时暂态波形;
图7为电网频率波动时改进下垂控制实验波形中,电网频率升高时实验波形;
图8为电网频率波动时改进下垂控制实验波形中,电网频率降低时实验波形。
具体实施方式
具体实施方式一:参见图1和2说明本实施方式,本实施方式所述的新型光伏并网逆变器电压型控制方法,
步骤一:将闭环反馈控制方法引入到下垂控制中,增加功率环,经拉普拉斯变换后,获得下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0
步骤二:将下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0分别叠加到新型光伏并网逆变器电压型的并网下垂方程,并简化后获得方程,
U = U 0 + ( k pp + k ip s ) ( P PV - P ) - - - ( 7 ) ,
f = f 0 + ( k pq + k iq s ) Q - - - ( 8 ) ,
将公式(7)和(8)中的功率环替换为直流母线电压环,获得
U = U 0 - ( k pp + k ip s ) ( U DCref - U DC ) - - - ( 9 ) ,
f = f 0 + ( k pq + k iq s ) Q - - - ( 8 ) ,
完成光伏并网逆变器对直流母线电压的控制;
其中,UDC为直流母线电压,UDCref为直流母线电压给定值,U为逆变器输出电压幅值参考值,f为逆变器输出频率参考值,U0为逆变器额定输出电压幅值,P为逆变器输出的有功功率,Q为逆变器输出的无功功率,kpp和kip分别为P-U下垂控制环PI调节器的比例、积分系数,kpq和kiq分别为Q-f下垂控制环PI调节器的比例、积分系数,s为复变量,PPV为光伏阵列最大功率。
本实施方式中,本发明所述的新型光伏并网逆变器电压型控制方法是针对低压电力传输线,所述低压电力传输线以阻性为主;
本发明在具体的应用过程中,用本发明所述的“下垂控制方法”替换现有技术中光伏逆变器电压型并网控制方法中的“下垂控制”,现有技术中光伏逆变器电压型并网控制方法的原理图具体参见图3;
且替换后的系统控制过程如下:前级升压电路采用MPPT控制,UPV和IPV分别为光伏输出电压和电流,经过最大功率点跟踪(MPPT)算法产生指令电压UPVref,减去光伏板输出电压后,经过光伏电压控制(PI)后,产生升压电路的开关管驱动信号。后级逆变器采用下垂控制,通过逆变器输出电压uAC和电流iAC计算出有功功率P和无功功率Q,据下垂方程不断微调频率f和幅值U,将重新生成的正弦信号uref作为电压环的给定,与逆变器输出电压uAC和电感电流iL进行电压电流双环控制,产生最终调制波,进行PWM调制,产生逆变桥开关管的驱动信号。
原理分析:
通过改进的下垂功率控制环,逆变器输出的有功功率P,最终能够稳定在给定值PPV,而与电网电压和线路阻抗无关;逆变器输出的无功功率Q能稳定在给定值0,进一步去除了电网频率的影响。
两级式光伏逆变器存在维持直流母线电压稳定的问题,并网运行时,升压变换器用来实现MPPT控制,所以稳定直流母线电压需要通过后级逆变器来完成。通常,下垂方程仅有功率环,不具有控制直流母线电压的功能,因此本发明提出将直流母线电压环加入下垂控制中,以满足两级式光伏逆变器的控制要求。
下面解释式(9)和(8)的物理意义:光伏逆变器前后级输出功率的不一致会导致直流母线电压的变化。若后级输出功率小于光伏输出功率,则多余的能量储存在直流母线电容CDC中,引起UDC持续上升;反之若后级输出功率大于光伏输出的功率,则不足的能量从直流母线电容CDC中取出,引起UDC持续下降;逆变器输出无功功率被控制为零,有功功率的大小则可在直流母线电压UDC上体现出来,所以式(9)和(8)能起到控制逆变器输出最大有功功率的作用,同时维持直流母线电压UDC稳定。
具体实施方式二:参见图1和2说明本实施方式,本实施方式与具体实施方式一所述的新型光伏并网逆变器电压型控制方法的区别在于,所述新型光伏并网逆变器的电压型并网下垂方程为:
U=U0-kp(P-P0)(1),
f=f0+kq(Q-Q0)(2),
其中,kp为逆变器输出有功功率的下垂系数,kq为逆变器输出无功功率的下垂系数,Q0为逆变器输出的额定无功功率,P0为逆变器输出的额定有功功率。
具体实施方式三:参见图1和2说明本实施方式,本实施方式与具体实施方式一所述的新型光伏并网逆变器电压型控制方法的区别在于,
所述下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0中,
Δ U 0 = ( k pp + k ip s ) ( P PV - P ) - - - ( 3 ) ,
Δf 0 = ( k pq + k iq s ) ( Q - 0 ) - - - ( 4 ) .
具体实施方式四:参见图1和2说明本实施方式,本实施方式与具体实施方式一所述的新型光伏并网逆变器电压型控制方法的区别在于,
将下垂曲线U0的平移量△U0和f0的平移量△f0分别叠加到新型光伏并网逆变器电压型的并网下垂方程中,获得方程:
U = U 0 + ( k pp + k ip s ) ( P PV - P ) - k p ( P - P PV ) - - - ( 5 ) ,
f = f 0 + ( k pq + k iq s ) ( Q - 0 ) + k q ( Q - 0 ) - - - ( 6 ) ,
然后将上述两个方程进行简化处理,获得公式(7)和(8)。
为了验证本发明方法的实用性,采用两级式光伏并网逆变器拓扑结构,设计了一台额定功率1kW的实验样机。电网电压220VAC,工频50Hz,开关频率为16.6kHz。
正弦波uinv为逆变器输出电压,正弦波iinv为逆变器输出电流,正弦波iload为负载电流,正弦波ugrid为电网电压,正弦波igrid为并网电流,由图4所示的光伏功率波动实验波形可见,当光伏板输出功率变化时,所提出的改进下垂控制策略能够有效地进行功率控制,使逆变器始终跟随给定功率输出,充分利用可再生能源。
图5为电网电压波动时改进下垂控制实验波形中,电网电压升高时暂态波形;图6为电网电压波动时改进下垂控制实验波形中,电网电压降低时暂态波形;由图5和6可见,所提出的控制策略能够克服电网电压扰动引起的功率偏差,使逆变器始终按照给定有功功率输出。
图7为电网频率波动时改进下垂控制实验波形中,电网频率升高时实验波形;图8为电网频率波动时改进下垂控制实验波形中,电网频率降低时实验波形;图7和8可见,所提出的控制策略能够有效克服电网频率波动的影响,保证逆变器输出的无功功率始终稳定在给定值0,使逆变器以单位功率因数并网,符合光伏逆变器并网的要求。

Claims (4)

1.光伏并网逆变器电压型控制方法,其特征在于,它包括下垂控制方法,所述的下垂控制方法包括如下步骤:
步骤一:将闭环反馈控制方法引入到下垂控制中,增加功率环,经拉普拉斯变换后,获得下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0
步骤二:将下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0分别叠加到新型光伏并网逆变器电压型的并网下垂方程,并简化后获得方程,
U = U 0 + ( k p p + k i p s ) ( P P V - P ) - - - ( 7 ) ,
f = f 0 + ( k p q + k i q s ) Q - - - ( 8 ) ,
将公式(7)和(8)中的功率环替换为直流母线电压环,获得
U = U 0 - ( k p p + k i p s ) ( U D C r e f - U D C ) - - - ( 9 ) ,
f = f 0 + ( k p q + k i q s ) Q - - - ( 8 ) ,
完成光伏并网逆变器对直流母线电压的控制;
其中,UDC为直流母线电压,UDCref为直流母线电压给定值,U为逆变器输出电压幅值参考值,f为逆变器输出频率参考值,U0为逆变器额定输出电压幅值,P为逆变器输出的有功功率,Q为逆变器输出的无功功率,kpp和kip分别为P-U下垂控制环PI调节器的比例、积分系数,kpq和kiq分别为Q-f下垂控制环PI调节器的比例、积分系数,s为复变量,PPV为光伏阵列最大功率。
2.根据权利要求1所述的光伏并网逆变器电压型控制方法,其特征在于,所述新型光伏并网逆变器的电压型并网下垂方程为:
U=U0-kp(P-P0)(1),
f=f0+kq(Q-Q0)(2),
其中,kp为逆变器输出有功功率的下垂系数,kq为逆变器输出无功功率的下垂系数,Q0为逆变器输出的额定无功功率,P0为逆变器输出的额定有功功率。
3.根据权利要求1所述的光伏并网逆变器电压型控制方法,其特征在于,所述下垂曲线U0的平移量△U0和逆变器额定输出频率f0的平移量△f0中,
ΔU 0 = ( k p p + k i p s ) ( P P V - P ) - - - ( 3 ) ,
Δf 0 = ( k p q + k i q s ) ( Q - 0 ) - - - ( 4 ) .
4.根据权利要求2所述的光伏并网逆变器电压型控制方法,其特征在于,将下垂曲线U0的平移量△U0和f0的平移量△f0分别叠加到新型光伏并网逆变器电压型的并网下垂方程中,获得方程:
U = U 0 + ( k p p + k i p s ) ( P P V - P ) - k p ( P - P P V ) - - - ( 5 ) ,
f = f 0 + ( k p q + k i q s ) ( Q - 0 ) + k q ( Q - 0 ) - - - ( 6 ) ,
然后将上述方程(5)和方程(6)进行简化处理,获得公式(7)和(8)。
CN201410421837.3A 2014-08-25 2014-08-25 新型光伏并网逆变器电压型控制方法 Expired - Fee Related CN104135033B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410421837.3A CN104135033B (zh) 2014-08-25 2014-08-25 新型光伏并网逆变器电压型控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410421837.3A CN104135033B (zh) 2014-08-25 2014-08-25 新型光伏并网逆变器电压型控制方法

Publications (2)

Publication Number Publication Date
CN104135033A CN104135033A (zh) 2014-11-05
CN104135033B true CN104135033B (zh) 2016-04-06

Family

ID=51807609

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410421837.3A Expired - Fee Related CN104135033B (zh) 2014-08-25 2014-08-25 新型光伏并网逆变器电压型控制方法

Country Status (1)

Country Link
CN (1) CN104135033B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105634309B (zh) * 2014-11-06 2018-06-22 台达电子工业股份有限公司 一种用于逆变系统的控制方法及控制装置
CN104377735A (zh) * 2014-12-08 2015-02-25 国家电网公司 一种分布式电源并网逆变器的下垂控制算法
CN104601029B (zh) * 2015-01-14 2016-10-19 重庆大学 逆变器并联控制系统下垂控制方法
CN104935002B (zh) * 2015-05-29 2017-11-03 国家电网公司 一种分布式光伏电源并网发电过电压调节方法
CN104868502B (zh) * 2015-06-16 2018-07-17 北京亿利智慧能源科技有限公司 一种适用于微电网的分布式光伏逆变器控制方法
CN105720607B (zh) * 2016-04-14 2018-04-24 哈尔滨工业大学 具有弹性有功功率分配的双下垂控制的逆变器并联控制方法
CN106249803B (zh) * 2016-09-30 2018-09-04 上海空间电源研究所 基于mppt的太阳能电池功率控制电路
CN106300431A (zh) * 2016-11-07 2017-01-04 燕山大学 一种微电网功率下垂控制方法
CN106771879B (zh) * 2017-01-18 2020-08-28 阳光电源股份有限公司 一种电网状态检测方法、装置及静态切换开关柜
CN107508463B (zh) * 2017-08-14 2019-11-05 华为技术有限公司 光伏直流变换器输出限压方法和装置
CN107947240A (zh) * 2017-12-21 2018-04-20 燕山大学 一种下垂控制并网逆变器的孤岛检测方法
CN108306311B (zh) * 2018-02-09 2021-12-07 南京工程学院 直流负载系统分区间响应电网调频需求的控制系统及方法
CN110277793B (zh) * 2018-03-16 2023-01-17 特变电工新疆新能源股份有限公司 一种两级式光伏逆变器主动功率平衡控制方法
CN111211704A (zh) * 2020-02-17 2020-05-29 中车青岛四方车辆研究所有限公司 不同辅助逆变器并联的控制方法及系统、辅助电源系统
CN111446735B (zh) * 2020-04-07 2022-01-28 清华大学 用于光伏中压集散系统的控制系统及方法
CN112165125B (zh) * 2020-10-09 2022-10-14 国电南瑞科技股份有限公司 一种惯量反下垂控制方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110285130A1 (en) * 2009-01-30 2011-11-24 Jan Thisted Power System Frequency Inertia for Wind Turbines
CN102611138B (zh) * 2012-03-20 2013-11-06 湖南大学 一种无延时的单相光伏并网功率调节方法
CN103944182B (zh) * 2014-04-01 2016-01-13 浙江大学 一种并网状态下储能逆变器pq控制的方法

Also Published As

Publication number Publication date
CN104135033A (zh) 2014-11-05

Similar Documents

Publication Publication Date Title
CN104135033B (zh) 新型光伏并网逆变器电压型控制方法
Lakshmi et al. Decoupled control of grid connected photovoltaic system using fractional order controller
Huka et al. A comprehensive LVRT strategy of two-stage photovoltaic systems under balanced and unbalanced faults
Maknouninejad et al. Analysis and control of PV inverters operating in VAR mode at night
US9490724B2 (en) Methods for controlling electrical inverters and electrical inverters and systems using the same
CN104158169B (zh) 一种光伏直流微网母线电压控制方法
US9948211B2 (en) System and method for controlling the operating area of an inverter coupled to an alternative energy source
CN103490420A (zh) 一种基于自抗扰技术的微电网并网控制方法
CN103929054A (zh) 一种用于光伏并网逆变器的软启动方法
Singh et al. IRPT based control of a 50 kW grid interfaced solar photovoltaic power generating system with power quality improvement
CN105720607B (zh) 具有弹性有功功率分配的双下垂控制的逆变器并联控制方法
Wei et al. Grid-connected PV system based on the series Z-source inverter
Kumar et al. Dynamic operation of ac micro-grid integrated pv-bess using 3-phase97-level chbmli
Modi et al. A maximum correntropy criteria based adaptive algorithm for an improved power quality SPV system
Razi et al. A new battery charger for plug-in hybrid electric vehicle application using back to back converter in a utility connected micro-grid
Méndez et al. Multifunctional grid-connected photovoltaic-system controlled by sliding mode
CN103280835A (zh) 三相并网光伏逆变器的发电状态控制方法
Dehghan et al. Grid connected DG systems based on Z-source NPC inverters
CN101969274A (zh) 一种母线电压稳定控制装置
Sharma et al. Supercapacitor utilization for power smoothening and stability improvement of ahybrid energy system in a weak grid environment
Prompinit et al. Ramp rate consideration of a BESS using active power control for PV generation
Isen Modelling and simulation of hysteresis current controlled single-phase grid-connected inverter
Leong et al. Deadbeat-based PI controller for stand-alone single-phase voltage source inverter using battery cell as primary sources
Rezayi et al. PLL-Less active and reactive power controller for three-phase grid-connected power converters
Lei et al. A universal droop-based control strategy for single-phase two-stage PV inverters

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: 20160406

Termination date: 20160825

CF01 Termination of patent right due to non-payment of annual fee