CN104779642A - 一种双馈风力发电机组的频率与阻尼综合控制方法 - Google Patents
一种双馈风力发电机组的频率与阻尼综合控制方法 Download PDFInfo
- Publication number
- CN104779642A CN104779642A CN201510182041.1A CN201510182041A CN104779642A CN 104779642 A CN104779642 A CN 104779642A CN 201510182041 A CN201510182041 A CN 201510182041A CN 104779642 A CN104779642 A CN 104779642A
- Authority
- CN
- China
- Prior art keywords
- frequency
- wind power
- control
- generation unit
- power generation
- 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
- 238000013016 damping Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 12
- 238000010248 power generation Methods 0.000 title abstract 10
- 230000008859 change Effects 0.000 claims description 13
- 238000011217 control strategy Methods 0.000 claims description 7
- 230000010355 oscillation Effects 0.000 abstract description 8
- 230000004044 response Effects 0.000 abstract description 8
- 230000001105 regulatory effect Effects 0.000 abstract description 2
- 230000001629 suppression Effects 0.000 abstract 1
- 230000001276 controlling effect Effects 0.000 description 5
- 230000001360 synchronised effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group 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
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Classifications
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种双馈风力发电机组的频率与阻尼综合控制方法。当双馈风力发电机组正常运行时,采用最大功率跟踪控制策略,当系统频率出现变化后,通过频率控制器,为风电机组参与系统调频提供快速有功控制指令,模拟风电机组的惯性频率响应。它将现有风电机组频率控制中比例微分控制改进为带有一阶惯性的比例微分控制,使得风电机组频率控制始终表现为正阻尼特性,保证了风电机组兼有频率调节和抑制系统功率振荡两种控制功能,为控制器广泛应用奠定基础。在电网出现频率变化后,频率控制快速调节风电机组的电磁功率,参与电网调频,若系统出现区域间低频振荡,在该控制策略下,风电机组能够抑制系统功率振荡。
Description
技术领域
本发明涉及一种发电机组的频率与阻尼综合控制方法,尤其是一种双馈风力发电机组的频率与阻尼综合控制方法。
背景技术
电网的频率调节能力和阻尼大小对电力系统的安全稳定运行有着重要影响,而对于风电装机比重较高的区域电网,这一问题将变得更加突出。目前,针对变速风电机组频率调节问题多采用比例微分(proportional derivative,PD)控制器,将df/dt微分控制环节和频率下垂特性的比例环节相结合,共同叠加到风电机组的最大功率追踪控制上。研究表明df/dt控制环节可在电网频率变化瞬间快速调节风电机组的电磁功率,释放或者吸收转子动能,提供惯性响应,而下垂控制可使风电机组参与电网的一次调频。然而,变速风电机组采用PD控制器产生附加注入功率进而实现频率调节的控制方案对电网的功率振荡产生了负阻尼影响,从而限制了控制器的推广应用。因此,通过改进现有变速恒频风电机组的频率控制策略,使其兼有频率调节和抑制系统功率振荡的能力,将对风电渗透率较高的区域电网的安全运行具有重要意义。
发明内容
本发明要解决的技术问题是提供一种双馈风力发电机组的频率与阻尼综合控制方法。
本发明所采取的技术方案是:
一种双馈风力发电机组的频率与阻尼综合控制方法,包括以下步骤:
a:判断电网是否稳定运行,如果是,采用最大功率跟踪控制策略对双馈风力发电机组进行有功控制,如果否,转向步骤b;
b:电网频率变化后,调节所述双馈风力发电机组的输出功率,有功增量ΔPf为:
其中,Kdf、Kpf分别为微分控制系数和比例控制系数,Δf为电网频率偏差,TD为惯性环节时间常数,s为复频率。
采用上述技术方案所产生的有益效果在于:它将原有的微分控制环节改变为带有一阶惯性滞后的微分控制环节,通过合理的参数配置,使双馈风力发电机组的频率控制始终表现为正阻尼特性,使风力发电机在参与电网频率调节的动态过程中具备增加系统阻尼的特性,消除了微分环节给电网引入负阻尼而造成的不利影响。它可使所述风电机组在电网频率变化或功率振荡时,通过参与电网的有功功率调节,使风电机组兼有频率调节和抑制系统功率振荡的能力,提高了电力系统的暂态稳定性。
附图说明
图1是本发明的流程图
图2是含双馈风电场的仿真拓扑结构图;
图3是含双馈风电场的系统等值电路图;
图4是本发明频率控制原理框图;
图5是本发明实施例1中风电机组电磁功率动态响应对比曲线;
图6是本发明实施例1中系统频率动态响应对比曲线;
图7是是本发明实施例1中发电机转子角动态响应对比曲线。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1:
如图1所示,一种双馈风力发电机组的频率与阻尼综合控制方法,包括以下步骤:
a:判断电网是否稳定运行,如果是,采用最大功率跟踪控制策略对双馈风力发电机组进行有功控制,如果否,转向步骤b;
b:电网频率变化后,调节所述双馈风力发电机组的输出功率,有功增量ΔPf为:
其中,Kdf、Kpf分别为微分控制系数和比例控制系数,Δf为电网频率偏差,TD为惯性环节时间常数,s为复频率。
为了说明所提改进策略的有效性,将双馈风电机组连接于图2所示的四机两区域电网B5母线处。假定风速为11m/s,双馈风电机组运行在最大功率跟踪运行状态。图5-7分别为系统负荷L1增加后,风电机组电磁功率动态响应特性曲线、系统频率动态响应特性曲线和同步发电机G3转子角动态响应特性曲线。
从图6可以看出,当系统负荷L1增加后,风电机组参加系统频率控制相比不参加系统频率控制时,系统最大频率偏差出现时刻延迟约2秒,最大频率偏差也明显减小,说明风电机组改善系统频率调节能力。在图7中,双馈风电机组动态注入有功功率后,发电机G3转子角摆动曲线衰减速度显著增加,风电机组表现为正阻尼特性,避免了风电调频引起的功率振荡,有效增强了系统动态稳定性。
本发明可使双馈风电机组根据电网频率偏差信号,快速调节风电机组的有功输出,利用风电机组独立、灵活的功率调节特性,在电网出现频率变化后,参与电网的频率调节,提高系统频率动态稳定性,并且可避免调频引起的系统功率振荡,适于风电场有功频率控制器的推广应用。
频率与阻尼控制原理分析如下:
参看图3,假定同步发电机G1采用q轴暂态电动势E′恒定的二阶经典模型且认为机械功率Pm恒定,其运动方程可表示为
式中,δ、TG、PG、Pm分别为G1的功角、惯性时间常数、电磁功率和机械功率;ωs、ωe分别为G1和系统的电角速度;DG为发电机阻尼系数。
对式(2)求小扰动量,则其小扰动运动方程为
TGp2Δδ+DGpΔδ+ΔPG=0 (3)
式中,ΔPG为G1有功变化量;p为微分因子。
根据相关知识,线路B1-B2和B2-B3输送的有功可表示为:
由系统功率平衡关系可得:
PΣ=PG+Pf (6)
式中,Pf为风电机组的有功输出。
为简化分析,假定线路母线电压恒定,因此线路上功率的变化仅取决于母线电压相角的变化。双馈风电机组转子侧变流器的无功电流控制环采用定交流电压控制,可认为风电场并网点电压UG恒定,对式(4)、(5)求小扰动量下线性方程:
式中,
参看图4,双馈风电机组在所提改进频率控制策略下,有功输出变化量ΔPf包含比例控制变化分量ΔPpf和微分控制变化分量ΔPdf两部分。为简化分析,以下分两种情形分析双馈风电机组在所提改进频率控制策略下是如何为发电机提供正阻尼。
①ΔPf仅包含比例控制变化分量ΔPpf,则
ΔPf=-KpfΔf=-KppΔδ (9)
式中,
将式(7)、(8)、(9)代入式(6)可得:
ΔPG=KKppΔδ+KK2Δδ (10)
式中,K=K1/(K1+K2)。
将式(10)代入式(3)可得发电机小扰动运动方程为
TGp2Δδ+(DG+KKp)pΔδ+KK2Δδ=0 (11)
由式(11)可知,Kp>0时,风电机组频率控制器中比例环节等效增大了同步发电机的阻尼系数,有利于阻尼发电机功角振荡。
②ΔPf仅包含微分控制变化分量ΔPdf,设发电机转子角以频率ωosc振荡,即s=jωosc,则ΔPf可表示为
式中
将式(12)、(7)、(8)代入式(6)可得
ΔPG=KK′dp2Δδ+KK′ppΔδ+KK2Δδ (13)
将(13)式代入(3)式得发电机小扰动运动方程为
(TG+KK′d)p2Δδ+(DG+KK′p)pΔδ+KK2Δδ=0 (14)
由式(14)可知,所提频率与阻尼综合控制器通过微分与一阶惯性环节共同模拟风电机组的惯性频率响应,控制环节中的微分环节在等效增加同步发电机惯性时间常数的同时,也使得同步发电机阻尼系数增加为DG+KK'p。因此,通过对惯性环节控制参数的设计可以实现风电机组在频率控制过程中,同时增强区域电网惯性和阻尼的控制目标,分析如下:
求式(14)所示微分方程的特征根,得其实部λ为
改进频率控制器若要增加系统阻尼,则λ应随Kdf的增加而减小,即此外还应保证λ随TD的增加而减小,即综上,解得惯性环节时间常数TD的取值范围为
所提频率与阻尼综合控制器中一阶惯性环节控制参数TD在上述取值范围内,即可实现前面所述双馈风电机组增加系统惯性和阻尼的控制目标。
本发明可以控制风电机组向系统动态注入的电磁功率。系统出现大幅频率波动时,风电机组可为电网提供功率支持,参与系统调频,并且该控制方案,可以有效增加电网阻尼,能够避免频率波动或其他扰动引起的区域间电网的功率振荡。在电网出现频率变化后,频率控制快速调节风电机组的电磁功率,参与电网调频,若系统出现区域间低频振荡,在该控制策略下,风电机组能够抑制系统功率振荡。
Claims (1)
1.一种双馈风力发电机组的频率与阻尼综合控制方法,其特征在于:包括以下步骤:
a:判断电网是否稳定运行,如果是,采用最大功率跟踪控制策略对双馈风力发电机组进行有功控制,如果否,转向步骤b;
b:电网频率变化后,调节所述双馈风力发电机组的输出功率,有功增量ΔPf为:
其中,Kdf、Kpf分别为微分控制系数和比例控制系数,Δf为电网频率偏差,TD为惯性环节时间常数,s为复频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510182041.1A CN104779642B (zh) | 2015-04-17 | 2015-04-17 | 一种双馈风力发电机组的频率与阻尼综合控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510182041.1A CN104779642B (zh) | 2015-04-17 | 2015-04-17 | 一种双馈风力发电机组的频率与阻尼综合控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104779642A true CN104779642A (zh) | 2015-07-15 |
CN104779642B CN104779642B (zh) | 2017-04-19 |
Family
ID=53620948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510182041.1A Expired - Fee Related CN104779642B (zh) | 2015-04-17 | 2015-04-17 | 一种双馈风力发电机组的频率与阻尼综合控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104779642B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105226723A (zh) * | 2015-11-19 | 2016-01-06 | 国家电网公司 | 一种双馈风电机组基于风功率跟踪自适应调节的虚拟惯性控制方法 |
CN105449699A (zh) * | 2016-01-11 | 2016-03-30 | 东北电力大学 | 双馈感应风电机组非线性分数阶自抗扰阻尼控制方法 |
CN106230021A (zh) * | 2016-08-16 | 2016-12-14 | 华北电力大学(保定) | 含可控惯性风电的区域互联电网的暂态功角稳定控制方法 |
CN106356903A (zh) * | 2016-09-21 | 2017-01-25 | 华北电力大学 | 一种基于风电机组和同步发电机控制频率的方法及装置 |
CN106684923A (zh) * | 2017-03-27 | 2017-05-17 | 云南电网有限责任公司 | 一种针对双馈风机的阻尼增强型控制方法 |
CN106786810A (zh) * | 2017-02-09 | 2017-05-31 | 国网浙江省电力公司电力科学研究院 | 一种频率下垂控制双馈风机的暂态有功支撑能力配置方法 |
CN108471120A (zh) * | 2018-03-06 | 2018-08-31 | 浙江运达风电股份有限公司 | 一种风电机组高电压穿越主动支撑优化控制方法 |
CN109687480A (zh) * | 2018-12-27 | 2019-04-26 | 国网江苏省电力有限公司电力科学研究院 | 抑制电力系统低频振荡的方法、系统及计算机存储介质 |
CN110098624A (zh) * | 2019-05-17 | 2019-08-06 | 华中科技大学 | 一种优化系统频率动态的双馈风机附加控制方法及系统 |
CN111478310A (zh) * | 2020-05-11 | 2020-07-31 | 山西大学 | 一种基于变下垂系数的直流配网虚拟惯性控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2254217A2 (en) * | 2006-12-22 | 2010-11-24 | Ingeteam Energy, S.A. | Reactive Power Control for Operating a Wind Farm |
CN102545246A (zh) * | 2012-01-18 | 2012-07-04 | 华北电力大学(保定) | 永磁直驱风力发电机组的阻尼控制方法 |
CN103441524A (zh) * | 2013-08-01 | 2013-12-11 | 国家电网公司 | 一种基于动态功率备用的变速风电机组频率控制方法 |
CN104065320A (zh) * | 2014-06-17 | 2014-09-24 | 合肥学院 | 双馈风力发电机组无速度传感器控制方法 |
-
2015
- 2015-04-17 CN CN201510182041.1A patent/CN104779642B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2254217A2 (en) * | 2006-12-22 | 2010-11-24 | Ingeteam Energy, S.A. | Reactive Power Control for Operating a Wind Farm |
CN102545246A (zh) * | 2012-01-18 | 2012-07-04 | 华北电力大学(保定) | 永磁直驱风力发电机组的阻尼控制方法 |
CN103441524A (zh) * | 2013-08-01 | 2013-12-11 | 国家电网公司 | 一种基于动态功率备用的变速风电机组频率控制方法 |
CN104065320A (zh) * | 2014-06-17 | 2014-09-24 | 合肥学院 | 双馈风力发电机组无速度传感器控制方法 |
Non-Patent Citations (4)
Title |
---|
唐西胜 等: "风力发电的调频技术研究综述", 《中国电机工程学报》 * |
曹张洁: "双馈感应风电机组参与系统一次调频的控制策略研究", 《工程科技Ⅱ辑》 * |
朱仲邃: "不完全微分PID算法在纯滞后系统中的应用", 《微计算机信息》 * |
杨列銮 等: "双馈异步风电机组模拟惯量相应控制技术综述", 《中国电力》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105226723A (zh) * | 2015-11-19 | 2016-01-06 | 国家电网公司 | 一种双馈风电机组基于风功率跟踪自适应调节的虚拟惯性控制方法 |
CN105226723B (zh) * | 2015-11-19 | 2018-02-09 | 国家电网公司 | 一种双馈风电机组基于风功率跟踪自适应调节的虚拟惯性控制方法 |
CN105449699A (zh) * | 2016-01-11 | 2016-03-30 | 东北电力大学 | 双馈感应风电机组非线性分数阶自抗扰阻尼控制方法 |
CN106230021B (zh) * | 2016-08-16 | 2018-12-04 | 华北电力大学(保定) | 含可控惯性风电的区域互联电网的暂态功角稳定控制方法 |
CN106230021A (zh) * | 2016-08-16 | 2016-12-14 | 华北电力大学(保定) | 含可控惯性风电的区域互联电网的暂态功角稳定控制方法 |
CN106356903A (zh) * | 2016-09-21 | 2017-01-25 | 华北电力大学 | 一种基于风电机组和同步发电机控制频率的方法及装置 |
CN106356903B (zh) * | 2016-09-21 | 2019-02-01 | 华北电力大学 | 一种基于风电机组和同步发电机控制频率的方法及装置 |
CN106786810A (zh) * | 2017-02-09 | 2017-05-31 | 国网浙江省电力公司电力科学研究院 | 一种频率下垂控制双馈风机的暂态有功支撑能力配置方法 |
CN106786810B (zh) * | 2017-02-09 | 2019-05-28 | 国网浙江省电力公司电力科学研究院 | 一种频率下垂控制双馈风机的暂态有功支撑能力配置方法 |
CN106684923A (zh) * | 2017-03-27 | 2017-05-17 | 云南电网有限责任公司 | 一种针对双馈风机的阻尼增强型控制方法 |
CN106684923B (zh) * | 2017-03-27 | 2018-02-27 | 云南电网有限责任公司 | 一种针对双馈风机的阻尼增强型控制方法 |
CN108471120A (zh) * | 2018-03-06 | 2018-08-31 | 浙江运达风电股份有限公司 | 一种风电机组高电压穿越主动支撑优化控制方法 |
CN109687480A (zh) * | 2018-12-27 | 2019-04-26 | 国网江苏省电力有限公司电力科学研究院 | 抑制电力系统低频振荡的方法、系统及计算机存储介质 |
CN110098624A (zh) * | 2019-05-17 | 2019-08-06 | 华中科技大学 | 一种优化系统频率动态的双馈风机附加控制方法及系统 |
CN110098624B (zh) * | 2019-05-17 | 2020-10-16 | 华中科技大学 | 一种优化系统频率动态的双馈风机附加控制方法及系统 |
CN111478310A (zh) * | 2020-05-11 | 2020-07-31 | 山西大学 | 一种基于变下垂系数的直流配网虚拟惯性控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104779642B (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104779642A (zh) | 一种双馈风力发电机组的频率与阻尼综合控制方法 | |
CN102074967B (zh) | 一种具有并网特性的储能型风电场控制方法 | |
CN110970933A (zh) | 基于主动支撑控制的光储联合发电系统虚拟惯量补偿方法 | |
CN107453410B (zh) | 负荷扰动的双馈风机参与风柴微网调频控制方法 | |
Yazdi et al. | Analytical modeling and inertia estimation of VSG-controlled Type 4 WTGs: Power system frequency response investigation | |
Zhang et al. | Comparison of inertia control methods for DFIG-based wind turbines | |
CN106208158B (zh) | 微网中多微源并列运行的惯量匹配方法 | |
CN108429289A (zh) | 一种基于虚拟同步发电机的控制方法及系统 | |
Chen et al. | On optimizing the transient load of variable-speed wind energy conversion system during the MPP tracking process | |
Zhang et al. | Control of DFIG-based wind farms for power network frequency support | |
Lao et al. | Innovated inertia control of DFIG with dynamic rotor speed recovery | |
Zeni et al. | Generic models of wind turbine generators for advanced applications in a VSC-based offshore HVDC network | |
Mi et al. | Small signal stability analysis of PMSG-VSG and optimal design for control parameters | |
CN105048511A (zh) | 一种含可控惯性风力发电机的发电系统惯性综合控制方法 | |
Colson et al. | Load-following for wind turbines with permanent magnet synchronous generators | |
Fadaeinedjad et al. | Flicker contribution of a wind turbine in a stand-alone wind diesel system | |
Mahmud et al. | Investigation of critical factors affecting dynamic stability of wind generation systems with permanent magnet synchronous generators | |
Han et al. | Supplementary power control of PMSG-based wind farms for system dynamic stability | |
Guo et al. | Exploration of the relationship between inertia enhancement and DC-link capacitance for grid-connected converters | |
Charles et al. | Field oriented control of Doubly Fed Induction Generator in wind power system | |
Dharmawardena et al. | Modeling variable speed wind turbine for power system dynamic studies | |
Ren et al. | A novel inertia control strategy for dfig wind turbines based on slip feedback | |
Wu et al. | A review of frequency regulation of DFIG-based wind farms | |
Conte et al. | Wind farm contribution to frequency regulation by distributed adaptive kinetic energy control | |
Ahuja et al. | Large scale wind energy conversion system using DFIG in sub-synchronous and super-synchronous speed ranges |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate 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: 20170419 |