CN101272121A - 一种风力发电机组的最大功率点跟踪方法 - Google Patents

一种风力发电机组的最大功率点跟踪方法 Download PDF

Info

Publication number
CN101272121A
CN101272121A CNA2008101060451A CN200810106045A CN101272121A CN 101272121 A CN101272121 A CN 101272121A CN A2008101060451 A CNA2008101060451 A CN A2008101060451A CN 200810106045 A CN200810106045 A CN 200810106045A CN 101272121 A CN101272121 A CN 101272121A
Authority
CN
China
Prior art keywords
wind
speed
conversion system
energy conversion
speed ratio
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
Application number
CNA2008101060451A
Other languages
English (en)
Other versions
CN101272121B (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.)
Institute of Electrical Engineering of CAS
Original Assignee
Institute of Electrical Engineering of CAS
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 Institute of Electrical Engineering of CAS filed Critical Institute of Electrical Engineering of CAS
Priority to CN2008101060451A priority Critical patent/CN101272121B/zh
Publication of CN101272121A publication Critical patent/CN101272121A/zh
Application granted granted Critical
Publication of CN101272121B publication Critical patent/CN101272121B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Eletrric Generators (AREA)
  • Wind Motors (AREA)

Abstract

一种风力发电机组的最大功率点跟踪方法,计算风能利用系数Cp、叶尖速比λ和
Figure 200810106045.1_AB_0
,当
Figure 200810106045.1_AB_1
,实际工作点在最大风能利用系数点的左侧,风力机的叶尖速比小于最佳叶尖速比,此时增加叶片转速;
Figure 200810106045.1_AB_2
,实际工作点在最大风能利用系数点的右侧,风力机的叶尖速比大于最佳叶尖速比,此时减小叶片转速;,实际工作点在最大风能利用系数点,风力机的实际叶尖速比等于最佳叶尖速比,风力机运行在最佳状态。将此最佳转速与风力机的实际转速相减,它们的差值经过PI电流控制转换成控制电流id *,id *与风力发电机组的实际电流id进行求差运算,差值输入到PWM调制,同时0与风力发电机组实际电流iq的差值也输入到PWM调制,调制结果控制驱动电路,实现最大功率点跟踪。

Description

一种风力发电机组的最大功率点跟踪方法
技术领域
本发明涉及风力发电机组的最大功率点跟踪方法。
背景技术
目前风电在电网中所占的比例越来越大。现有的风力发电机组大多都采用恒速运行方式,风力机的转速不能随风速的变化而变化,无法达到最优的风能利用效率。变速运行风力发电机组以其能最大限度的捕获风能的优点倍受国内外学者的青睐,文献1(D.S.Zinger,E.Muljadi.Annualized wind energy improvement using variable speeds[J],IEEE Trans.on Ind.Appl,1997,33:1444-1447.)证明变速风力发电系统比恒速系统效率要高38%,装机容量与日俱增。从国外安装运行的变速风力发电机组的运行状况看,目前国外企业生产的变速风力发电机组已经初步具备了最大功率点跟踪功能。国内在这方面的研究虽起步较晚,但发展较快,已有越来越多的学者开始了这个领域的研究,并取得了一定的成果,只是大部分研究还都处于理论研究和仿真验证阶段。
变速风力发电系统首要解决的问题是最大功率点跟踪(MPPT)技术。在过去的几十年里,对最大功率点跟踪提出过很多种控制算法,如爬山法,最佳功率(转矩)-转速曲线跟踪法,模糊控制法,自适应控制法。爬山法首先根据风速计算输出功率,然后对风力机的转速指令值ω*以一定的转速扰动值Δω进行扰动,并观测功率的变化,如果功率增加则转速扰动方向不变,否则转速反向。但是从理论上讲爬山法不可能追踪到真正的最大功率点,实际工作点将在最大功率点附近振荡,并且从现有文献上看,这种控制算法还没有经过大功率实验验证,也没有经过实际风力发电机组的现场验证。
最佳功率(转矩)一转速曲线跟踪法是根据风力机的空气动力学特性和风力发电机组机械与电气模型,可以推导出风力发电机组输出最佳电磁功率(最佳转矩)与风力机转速之间的关系,便可得到最佳功率(转矩)一转速曲线。控制风力发电机组按照这条曲线运行,实现最大风能捕获。但对风力机本身的参数有较强的依赖性,参数误差会引起较大的能量损失。对于大容量大转动惯量的风力机而言,这种控制方法的动态响应效果也不够好。模糊控制对模糊规则和隶属度函数的依赖性较强,而模糊规则和隶属度函数的选取至今还没有形成统一的科学方法。自适应控制具有不依赖系统参数和模型的优点,同时也具有良好的稳定性,只是模型较为复杂,实现难度较大。
发明内容
本发明克服了现有技术计算速度慢,对风力机本身参数依赖性强的缺点,提出一种新的最大功率点跟踪方法。本发明能够有效地跟踪最大功率点,并且实现简单。
风力机是把风能转化为机械能的一种能量转换装置,在转化过程中有一定的能量损失,存在一个能量转换系数Cp C p = P m P w 其中Pm为单位时间内叶轮吸收且转换的机械能,即风力机的机械输出功率,Pw为单位时间内通过叶轮扫掠面的风能,即风力机的输入功率。风能利用系数Cp与风速v,叶片转速nw,叶片直径D和桨叶节距角β均有关系,是叶尖速比λ和桨叶节距角β的函数。叶尖速比λ是风轮叶尖速度与风速之比:
λ = R w ω w v = π R w n w 30 v
ωw为风力机叶片旋转角速度,nw为叶片的转速,v为风速,Rw为风力机的叶轮半径。在桨叶节距角β不变时,风能利用系数只与叶尖速比λ有关。对于一台确定的风力机,在一定风速和桨距角时,存在一个最大风能转换系数Cpmax,此时风力机输出功率最大。因此,在某一风速下,调节风力机转速,使其运行在最佳叶尖速比条件下,就可以达到最大功率点跟踪的目的。
本发明采用增量电导法计算风力机叶片的最佳转速,判断风力机运行区域。增量电导法的基本思想是当 dC p d&lambda; > 0 时,λ<λopt即nw<nopt,风力机叶片的实际转速小于最佳转速,此时应增加叶片的转速;当 dC p d&lambda; < 0 时,λ>λopt,即nw>nopt,风力机叶片的实际转速大于最佳转速,此时应减小叶片的转速;当 dC p d&lambda; = 0 时,λ=λopt即nw=nopt,叶片的实际转速等于最佳转速,其中λopt,为对应于最大风能利用系数下的叶尖速比,nopt为对应于最大风能利用系数下的叶片转速。通过测量和比较,预估出最大风能转换系数Cpmax
本发明方法步骤如下:
首先根据风速和风力发电机组当前速度,运用增量电导法计算出风力机叶片的最佳转速。具体计算过程如下:在某一时刻利用采样得到风力机的机械输出功率Pm和风力机的输入功率Pw,利用公式 C p = P m P w 计算出风能利用系数Cp,并测量出此时风力机叶片的转速nw,利用公式 &lambda; = R w &omega; w v = &pi; R w n w 30 v 计算出此时刻的叶尖速比λ。利用Cp和λ计算出
Figure A20081010604500048
然后根据
Figure A20081010604500051
的值来判断风力机运行区域。若 dC p d&lambda; > 0 , 则说明风力机实际工作点在最大风能利用系数点的左侧,这时应增加叶片转速;若 dC p d&lambda; < 0 , 则风力机实际工作点在最大风能利用系数点的右侧,这时应减小叶片转速; dC p d&lambda; = 0 , 风力机实际工作点在最大风能利用系数点,此时的转速为最佳转速。计算出最佳转速后,将此最佳转速与风力机的实际转速相减,差值经过PI电流控制转换成控制电流id *,id *与风力发电机组的实际电流id进行求差运算,差值输入到PWM调制,同时0与风力发电机组实际电流iq的差值也输入到PWM调制,调制结果控制驱动电路,从而实现了最大功率点跟踪。
附图说明
图1风力机特性曲线;
图2风能利用系数与叶尖速比的关系曲线图;
图3增量电导法流程图;
图4具备最大功率点跟踪功能的风力发电机组。
具体实施方式:
风力机输出功率与转速之间的关系曲线如图1所示,不同风速情况下均存在一个转速使得风力机输出功率最大,即当风力机叶片转速为该特定值时就能实现最大功率点跟踪。如图2所示,对于一台确定的风力机,在一定风速和桨距角时,存在一个最大风能转换系数Cpmax,此时风力机输出功率最大。因此,在某一风速下,调节风力机叶片转速,使其运行在最佳叶尖速比条件下,就可以达到最大功率点跟踪的目的。
本发明首先利用采样得到风力机的机械输出功率Pm和风力机的输入功率Pw,利用公式 C p = P m P w 计算出风能利用系数Cp,并根据风速和风力发电机组的当前转速nw,利用公式 &lambda; = R w &omega; w v = &pi; R w n w 30 v 计算出此时刻的叶尖速比λ。利用Cp和λ计算出
Figure A20081010604500057
然后根据
Figure A20081010604500058
的值来判断风力机运行区域,再调节风力机叶片转速使其运行在最大风能利用系数点,从而达到最大风能捕获的目的。调节叶片转速的依据是: dC p d&lambda; > 0 , 风力机实际工作点在最大风能利用系数点的左侧,风力机的叶尖速比小于最佳叶尖速比,这时应增加叶片转速; dC p d&lambda; < 0 , 风力机实际工作点在最大风能利用系数点的右侧,风力机的叶尖速比大于最佳叶尖速比,这时应减小叶片转速; dC p d&lambda; = 0 , 风力机实际工作点在最大风能利用系数点,风力机的实际叶尖速比等于最佳叶尖速比,风力机运行在最佳状态。但是由于导数
Figure A20081010604500062
无法直接求出,所以工程上一般是采用近似代替
dC p d&lambda; &ap; &Delta; C p &Delta;&lambda; = C p - C pb &lambda; - &lambda; b
式中Cpb为风力机扰动前的风能利用系数,λb为风力机扰动前的叶尖速比。由此方法计算出叶片的最佳转速。得到此最佳转速后,与风力机的实际转速相减,它们的差值经过PI电流控制转换成一个控制电流id *,id *与风力发电机组的实际电流id进行求差运算,差值输入到PWM调制,同时0与风力发电机组实际电流iq的差值也输入到PWM调制,调制结果控制驱动电路,从而实现了最佳转速的跟踪和调整,即实现了最大功率点跟踪。
下面结合图3和图4描述本发明具体的最大功率点跟踪过程。
1、首先利用采样得到风力机的机械输出功率Pm和风力机的输入功率Pw,利用公式 C p = P m P w 计算出风能利用系数Cp,并根据此时风速和风力机叶片转速nw,利用公式 &lambda; = R w &omega; w v = &pi; R w n w 30 v 计算出此时刻的叶尖速比λ。
2、然后利用dCp=Cp-Cpb,dλ=λ-λb计算出dCp和dλ,由于dλ是分母,所以要首先判断dλ是否为0,
3、如果叶尖速比λ没有变化即dλ=0,且风能利用系数Cp也没有发生变化即dCp=0,则说明不需要进行调整;如果叶尖速比λ没有变化即dλ=0,但是风能利用系数发生了变化即dCp≠0,这时就应该根据dCp的值来进行叶片转速的调整从而改变叶尖速比。具体的是如果dCp>0则应该增加风力机叶片转速,从而使叶尖速比增大;如果dCp<0则应该减小风力机叶片转速,从而使叶尖速比减小。
4、假如dλ≠0再根据
Figure A20081010604500066
的值进行调整,根据附图1风能利用系数与叶尖速比的关系曲线可以知道,只有在最大风能利用系数点处才有 dC p d&lambda; = 0 . 所以如果 dC p d&lambda; = 0 , 说明风力机运行在最佳状态,不需要调整;如果 dC p d&lambda; > 0 , 说明风力机实际工作点在最大风能利用系数点的左侧,这时应增加叶片转速从而增加叶尖速比;如果 dC p d&lambda; < 0 , 说明风力机实际工作点在最大风能利用系数点的右侧,这时应减小叶片转速从而减小叶尖速比。
5、判断完这个过程后,把该时刻的风能利用系数Cp和叶尖速比λ赋给Cpb和λb为下一次循环提供数据。循环运算,直到找到 dC p d&lambda; = 0 的最佳转速。
6、确定出最佳转速后,将此最佳转速w*与风力机的实际转速w求差,差值Δw通过PI电流控制技术转换成一个控制电流id *。风力发电机的定子电流ia,ib,ic通过d,q变换算法得到电流id,iq。id *与id作差得到差值电流isd *,0与iq作差得到差值电流isq *。isd *和isq *一起输入到PWM调制,经过调制后,输出控制参量来驱动控制电路,使风力机在最佳转速状态下运行,从而实现了风力发电机组的最大功率点跟踪。

Claims (1)

1、一种风力发电机组的最大功率点跟踪方法,其特征在于:
(1)首先利用采样得到风力机的机械输出功率Pm和风力机的输入功率Pw,根据公式 C p = P m P w 计算出风能利用系数Cp,并根据此时风速v和风力机叶片转速nw,利用公式 &lambda; = R w &omega; w v = &pi; R w n w 30 v 计算出此时刻的叶尖速比λ,其中公式中的Rw为风力机的叶轮半径;
(2)然后利用dCp=Cp-Cpb,dλ=λ-λb计算出dCp和dλ,判断dλ是否为0,式中Cpb为前一时刻的风能利用系数,λb为前一时刻的叶尖速比,;
如果dλ=0,且dCp=0,不需要调整风力机叶片转速;如果dλ=0,dCp>0则增加风力机叶片转速,从而使叶尖速比增大;如果dλ=0,dCp<0,则减小风力机叶片转速,使叶尖速比减小;
(3)如dλ≠0, dC p d&lambda; = 0 , 说明风力机运行在最佳状态,不需要调整;如果dλ≠0, dC p d&lambda; > 0 , 说明风力机实际工作点在最大风能利用系数点的左侧,此时应增加叶片转速从而增加叶尖速比;如果dλ≠0, dC p d&lambda; < 0 说明风力机实际工作点在最大风能利用系数点的右侧,此时应减小叶片转速从而减小叶尖速比。
(4)判断完此过程后,把该时刻的风能利用系数Cp和叶尖速比λ赋给Cpb和λb为下一次循环提供数据;按照步骤(1)至(3)循环运算,直到找到 dC p d&lambda; = 0 的最佳转速;
(5)确定出最佳转速后,将此最佳转速w*与风力机的实际转速w求差,差值Δw通过PI电流控制技术转换成一个控制电流id *;风力发电机定子三相电流ia,ib,ic通过d,q变换算法得到d,q坐标系下的直轴电流分量id和交轴电流分量iq;id *与id作差得到差值电流isd *,0与iq作差得到差值电流isq *;isd *和isq *一起输入到PWM调制,经过调制后,输出控制参量来驱动控制电路,使风力发电机组在最佳转速状态下运行,实现了最大功率点跟踪。
CN2008101060451A 2008-05-07 2008-05-07 一种风力发电机组的最大功率点跟踪方法 Expired - Fee Related CN101272121B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101060451A CN101272121B (zh) 2008-05-07 2008-05-07 一种风力发电机组的最大功率点跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101060451A CN101272121B (zh) 2008-05-07 2008-05-07 一种风力发电机组的最大功率点跟踪方法

Publications (2)

Publication Number Publication Date
CN101272121A true CN101272121A (zh) 2008-09-24
CN101272121B CN101272121B (zh) 2012-05-30

Family

ID=40005863

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101060451A Expired - Fee Related CN101272121B (zh) 2008-05-07 2008-05-07 一种风力发电机组的最大功率点跟踪方法

Country Status (1)

Country Link
CN (1) CN101272121B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102242689A (zh) * 2011-06-24 2011-11-16 南京理工大学 基于风力发电的最大功率点跟踪控制的改进爬山算法
CN102307038A (zh) * 2011-09-13 2012-01-04 辽宁力迅风电控制系统有限公司 离网型风力发电系统的最大功率跟踪方法及其控制装置
CN102352814A (zh) * 2011-06-21 2012-02-15 北京交通大学 一种大型直驱风电机组的最大功率跟踪方法
CN102392785A (zh) * 2011-09-22 2012-03-28 南车株洲电力机车研究所有限公司 一种风力机最大功率跟踪曲线的线性化方法
CN102434391A (zh) * 2011-12-27 2012-05-02 南京理工大学 一种基于起始转速调整的改进最大功率点跟踪控制方法
CN103701378A (zh) * 2013-12-27 2014-04-02 内蒙古久和能源科技有限公司 一种风力发电机输出功率控制方法
CN109756162A (zh) * 2019-03-26 2019-05-14 无锡商业职业技术学院 一种永磁直驱风电系统模糊mppt控制方法及其控制装置
CN111022254A (zh) * 2019-12-25 2020-04-17 金陵科技学院 奇异摄动风力发电模型最大功率点跟踪的时滞控制方法
CN111963372A (zh) * 2020-09-01 2020-11-20 北京石油化工学院 一种风力发电机最佳转速跟踪控制方法
CN112343770A (zh) * 2020-11-17 2021-02-09 北京石油化工学院 一种基于观测器的风力发电机最佳转速有限时间跟踪控制方法
CN113323806A (zh) * 2021-06-11 2021-08-31 国电南京自动化股份有限公司 一种风电机组主控制系统的最大风能跟踪控制方法
CN114548611A (zh) * 2022-04-27 2022-05-27 东方电气风电股份有限公司 一种风力发电机组最优增益参数的搜索方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100409562C (zh) * 2006-11-07 2008-08-06 合肥工业大学 风力发电最大功率点跟踪控制方法

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102352814A (zh) * 2011-06-21 2012-02-15 北京交通大学 一种大型直驱风电机组的最大功率跟踪方法
CN102242689A (zh) * 2011-06-24 2011-11-16 南京理工大学 基于风力发电的最大功率点跟踪控制的改进爬山算法
CN102307038A (zh) * 2011-09-13 2012-01-04 辽宁力迅风电控制系统有限公司 离网型风力发电系统的最大功率跟踪方法及其控制装置
CN102307038B (zh) * 2011-09-13 2013-09-18 辽宁力迅风电控制系统有限公司 离网型风力发电系统的最大功率跟踪方法及其控制装置
CN102392785B (zh) * 2011-09-22 2014-05-21 南车株洲电力机车研究所有限公司 一种风力机最大功率跟踪曲线的线性化方法
CN102392785A (zh) * 2011-09-22 2012-03-28 南车株洲电力机车研究所有限公司 一种风力机最大功率跟踪曲线的线性化方法
CN102434391A (zh) * 2011-12-27 2012-05-02 南京理工大学 一种基于起始转速调整的改进最大功率点跟踪控制方法
CN103701378B (zh) * 2013-12-27 2016-11-02 内蒙古久和能源装备有限公司 一种风力发电机输出功率控制方法
CN103701378A (zh) * 2013-12-27 2014-04-02 内蒙古久和能源科技有限公司 一种风力发电机输出功率控制方法
CN109756162A (zh) * 2019-03-26 2019-05-14 无锡商业职业技术学院 一种永磁直驱风电系统模糊mppt控制方法及其控制装置
CN111022254A (zh) * 2019-12-25 2020-04-17 金陵科技学院 奇异摄动风力发电模型最大功率点跟踪的时滞控制方法
CN111963372A (zh) * 2020-09-01 2020-11-20 北京石油化工学院 一种风力发电机最佳转速跟踪控制方法
CN111963372B (zh) * 2020-09-01 2021-12-17 北京石油化工学院 一种风力发电机最佳转速跟踪控制方法
CN112343770A (zh) * 2020-11-17 2021-02-09 北京石油化工学院 一种基于观测器的风力发电机最佳转速有限时间跟踪控制方法
CN112343770B (zh) * 2020-11-17 2021-09-24 北京石油化工学院 一种基于观测器的风力发电机最佳转速有限时间跟踪控制方法
CN113323806A (zh) * 2021-06-11 2021-08-31 国电南京自动化股份有限公司 一种风电机组主控制系统的最大风能跟踪控制方法
CN113323806B (zh) * 2021-06-11 2022-08-19 国电南京自动化股份有限公司 一种风电机组主控制系统的最大风能跟踪控制方法
CN114548611A (zh) * 2022-04-27 2022-05-27 东方电气风电股份有限公司 一种风力发电机组最优增益参数的搜索方法
CN114548611B (zh) * 2022-04-27 2022-07-19 东方电气风电股份有限公司 一种风力发电机组最优增益参数的搜索方法

Also Published As

Publication number Publication date
CN101272121B (zh) 2012-05-30

Similar Documents

Publication Publication Date Title
CN101272121B (zh) 一种风力发电机组的最大功率点跟踪方法
CN101603503B (zh) 一种定桨距风力机的内模控制方法
Marques et al. A survey on variable-speed wind turbine system
CN102332727B (zh) 一种利用直流侧飞轮储能单元平滑永磁直驱风力发电系统输出有功功率的方法
CN101581272B (zh) 定桨距变速风力发电机组在失速区的功率控制方法
CN101769232A (zh) 一种定桨距变速风力发电机组的全风速功率控制方法
CN101404476A (zh) 并网变速恒频风力发电机组运行控制方法
Bystryk et al. Small wind turbine power control in intermittent wind gusts
Cho et al. Development and experimental verification of counter-rotating dual rotor/dual generator wind turbine: Generating, yawing and furling
CN104795835A (zh) 一种双馈风电变流器开关频率的控制方法及系统
Du et al. MPPT control strategy of DFIG-based wind turbines using double steps hill climb searching algorithm
CN102156044B (zh) 适用于直驱型风电机组测试的风力机模拟器选型方法
Shao et al. Application research of maximum wind-energy tracing controller based adaptive control strategy in WECS
CN105257475A (zh) 一种失速型风力发电机组的控制方法
CN107882680B (zh) 一种用于风力发电机组的转速控制方法
CN108757306A (zh) 一种风力发电机控制方法、装置以及风力发电机
CN103375332B (zh) 变速变桨风力发电机组最优阻力矩动态优化方法
Zhen et al. Study on control strategy of maximum power capture for DFIG in wind turbine system
CN102013762A (zh) 风力发电机组的变速恒频装置
CN109882355A (zh) 一种提升中低风速下风力发电效率的动态匹配方法及系统
Barakati Wind turbine systems: history, structure, and dynamic model
Xiao et al. VSCF wind turbine control strategy for maximum power generation
Gavgani et al. Exploiting bidirectional power flow control to capture wind gust power in small and medium wind turbines
Yang et al. Hill-Climbing Algorithm for the Wind Turbine Yaw System
CN104863793A (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
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Xu Honghua

Inventor after: Li Jianlin

Inventor after: Zhao Bin

Inventor after: Ma Rui

Inventor before: Li Jianlin

Inventor before: Zhao Bin

Inventor before: Ma Rui

Inventor before: Xu Honghua

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: LI JIANLIN ZHAO BIN MA RUI XU HONGHUA TO: XU HONGHUA LI JIANLIN ZHAO BIN MA RUI

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

Termination date: 20180507

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