CN110007183A - 逆变器在线开路故障检测方法 - Google Patents
逆变器在线开路故障检测方法 Download PDFInfo
- Publication number
- CN110007183A CN110007183A CN201910281011.4A CN201910281011A CN110007183A CN 110007183 A CN110007183 A CN 110007183A CN 201910281011 A CN201910281011 A CN 201910281011A CN 110007183 A CN110007183 A CN 110007183A
- Authority
- CN
- China
- Prior art keywords
- phase
- inverter
- current
- open
- circuit fault
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 52
- 238000005070 sampling Methods 0.000 claims abstract description 37
- 238000012545 processing Methods 0.000 claims abstract description 15
- 230000004888 barrier function Effects 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 6
- 230000004807 localization Effects 0.000 abstract description 4
- 239000003990 capacitor Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16566—Circuits and arrangements for comparing voltage or current with one or several thresholds and for indicating the result not covered by subgroups G01R19/16504, G01R19/16528, G01R19/16533
- G01R19/16576—Circuits and arrangements for comparing voltage or current with one or several thresholds and for indicating the result not covered by subgroups G01R19/16504, G01R19/16528, G01R19/16533 comparing DC or AC voltage with one threshold
- G01R19/1658—AC voltage or recurrent signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Inverter Devices (AREA)
Abstract
本发明公开了一种逆变器在线开路故障检测方法,首先采集逆变器的两相电流,然后计算得到采样时刻时三相电流的幅值和相位,再据此计算得到采样时刻时三相电流绝对值之和的标幺值,将此标幺值与预设阈值进行比较,判断是否发生开路故障,如果发生再根据三相电流进行故障定位,其具体方法为对三相电流进行标幺化处理,然后对标幺化电流进行过零置“一”处理,计算采样时刻的故障定位标识,将故障定位标识与预设阈值进行比较,实现开路故障相与故障管的定位检测。本发明通过对检测得到的负载电流进行在线处理进行开路故障检测,减小故障检测时间,实现在线故障检测。
Description
技术领域
本发明属于逆变器故障检测技术领域,更为具体地讲,涉及一种逆变器在线开路故障检测方法。
背景技术
逆变器广泛应用于电网、电机系统和电源等领域,但是系统因半导体功率器件发生故障后所带来的隐患不可忽视,特别是对于不易察觉的逆变器开路故障。目前,通过检测逆变器输出电流,实现快速准确的逆变器故障检测,而且逆变器故障检测和容错控制的关键方法和技术,可以推广到电网逆变器、电机逆变器和电源逆变器故障检测容错中,得到了国内外专家学者的重视。
基于电流的逆变器故障检测方法,通过对已测电流量经过坐标变换、信号处理或者模式识别等技术进行检测以及定位故障。基于负载电流的故障检测方法可以利用电流闭环系统中的传感器,无需额外的电流传感器,但是电流量却是一个状态量,而且容易受到系统运行情况的影响,不能直接及时地将故障信息反映出来,使得其极易受到外界变化的影响,导致错报和漏报故障,抗干扰性能较差,检测可靠性较低。除此之外,基于电流的现有诸多技术在故障检测时间上往往受限于电流基波周期,检测时间通常以电流基波周期为参考,检测时间较长。
发明内容
本发明的目的在于克服现有技术的不足,提供一种逆变器在线开路故障检测方法,通过对检测得到的负载电流进行在线处理进行开路故障检测,减小故障检测时间,实现在线故障检测。
为实现上述发明目的,本发明逆变器在线开路故障检测方法包括以下步骤:
S1:采用传感器采集逆变器的两相电流,分别记为ia(k)、ib(k),k表示传感器的采样时刻;
S2:记第三相电流为ic(k),根据以下公式计算采样时刻k时三相电流的幅值Im和相位ωt:
其中,ω表示电流角频率,t表示时间,与采样时刻k对应;
S3:采用步骤S2计算得到的采样时刻k时三相电流的幅值Im和相位ωt计算得到采样时刻k时三相电流绝对值之和ρ(k)的标幺值fρ(k):
其中,n是满足式的正整数;
S4:根据需要设定检测阈值s1,若fρ(k)>s1,则说明采样时刻k时逆变器未发生开路故障,不作任何操作,若fρ(k)≤s1,则逆变器开路故障发生,进入步骤S5进行故障定位;
S5:根据以下公式对采样时刻k的三相电流进行标幺化处理,得到标幺化三相电流iph(k):
其中,h=a,b,c,当h=a,j=0,当h=b,j=1,当h=c,j=2;
S6:根据需要设定阈值s2,采用以下公式对标幺化电流iph(k)进行过零置“一”处理,得到标幺化电流iph(k)对应的标志参数σh(k):
S7:采用以下公式计算采样时刻k的故障定位标识δh(k):
其中,L表示一个电流基波周期的采样点数;
S8:根据需要设定检测阈值s3,若δh(k)≤-s3,则逆变器h相发生上管开路故障;若|δh(k)|<s3,则逆变器h相未发生开路故障;若δh(k)≥s3,则逆变器h相发生下管开路故障。
本发明逆变器在线开路故障检测方法,首先采集逆变器的两相电流,然后计算得到采样时刻时三相电流的幅值和相位,再据此计算得到采样时刻时三相电流绝对值之和的标幺值,将此标幺值与预设阈值进行比较,判断是否发生开路故障,如果发生再根据三相电流进行故障定位,其具体方法为对三相电流进行标幺化处理,然后对标幺化电流进行过零置“一”处理,计算采样时刻的故障定位标识,将故障定位标识与预设阈值进行比较,实现开路故障相与故障管的定位检测。
本发明具有以下有益效果:
(1)本发明针对现有技术在故障检测时间上往往受限于电流基波周期的问题,利用检测得到的负载电流设置新的变量,使得检测时间仅仅与1/6的电流基波周期相关,减小故障检测时间,提高故障检测的速度,改善故障检测性能;
(2)本发明能够准确定位故障相和故障管,为系统进一步的故障处理提供有效的依据;
(3)本发明不受负载变化的影响,逆变器负载可以是电网、电机、电源等,相比于其它在线检测方法有更好的灵敏度和准确性,非常适用于逆变器快速开路故障在线检测;
(4)本发明适用于三相逆变器单管和多管开路故障的在线监测。
附图说明
图1是本发明逆变器在线开路故障检测方法的具体实施方式流程图;
图2是本实施例中三相逆变器系统结构图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明逆变器在线开路故障检测方法的具体实施方式流程图。如图1所示,本发明逆变器在线开路故障检测方法的具体步骤包括:
S101:采集两相电流:
采用传感器采集逆变器的两相电流,分别记为ia(k)、ib(k),k表示传感器的采样时刻。
S102:计算三相电流参数:
记第三相电流为ic(k),由于三相电流ia、ib、ic满足等式ia+ib+ic=0。那么根据以下公式计算采样时刻k时三相电流的幅值Im和相位ωt:
其中,ω表示电流角频率,t表示时间,与采样时刻k对应。
S103:计算三相电流绝对值之和的标幺值:
记三相电流绝对值之和ρ=|ia|+|ib|+|ic|,则有:
其中,n是满足式的正整数。
记T1为1/6个电流基波周期,则三相电流绝对值之和ρ在每个周期T1内的平均值<ρ>可以采用以下公式计算:
三相电流绝对值之和ρ可以采用以下公式进行标幺化处理,得到三相电流绝对值之和ρ的标幺值fρ:
显然,采用步骤S102计算得到的采样时刻k时三相电流的幅值Im和相位ωt即可计算得到采样时刻k时三相电流绝对值之和ρ(k)的标幺值fρ(k):
在逆变器正常工作情况下,标幺值fρ的取值范围为可近似认为fρ∈(0.9069,1.0472)。当逆变器单桥臂开路故障发生后,标幺值fρ将会脱离区间(0.9069,1.0472)。
S104:故障检测:
根据需要设定检测阈值s1,若fρ(k)>s1,则说明采样时刻k时逆变器未发生开路故障,不作任何操作,若fρ≤s1,则逆变器开路故障发生,进入步骤S105进行故障定位。
在实际应用时,综合考虑正常与故障情况下标幺值fρ的变化,并权衡检测方法的快速性与抗干扰性,避免故障误判的发生,适当选取阈值s1。根据标幺值fρ的取值范围可知,阈值s1应小于0.9069。经过实验发现,为了兼顾可靠性与检测速度,阈值s1可在0.5附近进行取值。由于变量fρ采用了标准化处理,使得其不依赖于负载,对控制系统参数特征不敏感,故阈值s1不会因故障的不同而改变。
S105:三相电流标幺化处理:
根据以下公式对采样时刻k的三相电流进行标幺化处理,得到标幺化三相电流iph(k):
其中,h=a,b,c,当h=a,j=0,当h=b,j=1,当h=c,j=2。
可以得到标幺化三相电流ipa(k)、ipb(k)、ipc(k)的计算公式如下:
S106:标幺化电流过零置“一”处理:
根据需要设定阈值s2,采用以下公式对标幺化电流iph(k)进行过零置“一”处理,得到标幺化电流iph(k)对应的标志参数σh(k):
阈值s2的理论值为0,但是在实际应用中通常取标幺化电流iph(k)的最大值的2%附近,即s2≈0.01。
S107:计算故障定位标识:
采用以下公式计算采样时刻k的故障定位标识δh(k):
其中,L表示一个电流基波周期的采样点数。
也就是说,本发明对于采样时刻k-L+1至采样时刻k中的每一个采样时刻k′,先根据对应电流基波周期[k′-L+1,k′]内的每个采样时刻k″标幺化电流对应的标志参数计算其参数
然后将采样时刻k-L+1至采样时刻k中的每一个采样时刻k′对应的求和后再平均,从而得到故障定位标识δh(k)。
S107:故障定位:
根据需要设定检测阈值s3,若δh(k)≤-s3,则逆变器h相发生上管开路故障;;若|δh(k)|<s3,则逆变器h相未发生开路故障;若δh(k)≥s3,则逆变器h相发生下管开路故障。从而实现开路故障相与故障管的定位检测。
阈值s3的取值应小于1,实际应用时为了兼顾可靠性与检测速度,可以0.5附近进行取值。
为了更好地说明本发明的技术效果,采用一个具体实例对本发明进行详细说明。图2是本实施例中三相逆变器系统结构图。如图2所示,本实施例中三相逆变器包括由功率MOS管T1-T6构成的三相全桥电路,以及三相负载。图2中N为星形联结的三相负载中点,G为直流侧的接地点,L1-L3为电感,C1-C3为电容,共同构成滤波电路,C4和C5为直流侧的等值分压电容,O为分压电容的中点。
为了便于构建开路故障检测表,在故障检测时,若fρ(k)>s1,令故障检测标志Fρ=0,若fρ≤s1,令故障检测变量标志Fρ=1。在故障定位时,若δh(k)≤-s3,令故障定位标志Δh=1,若|δh(k)|<s3,令故障定位标志Δh=0,若δh(k)≥s3,令故障定位标志Δh=-1。表1是本实施例中三相逆变器的开路故障检测表。
故障功率管 | F<sub>ρ</sub> | Δ<sub>a</sub> | Δ<sub>b</sub> | Δ<sub>c</sub> |
正常 | 0 | 0 | 0 | 0 |
T1 | 1 | 1 | 0 | 0 |
T2 | 1 | -1 | 0 | 0 |
T3 | 1 | 0 | 1 | 0 |
T4 | 1 | 0 | -1 | 0 |
T5 | 1 | 0 | 0 | 1 |
T6 | 1 | 0 | 0 | -1 |
T1&T4 | 1 | 1 | -1 | 0 |
T1&T6 | 1 | 1 | 0 | -1 |
T2&T3 | 1 | -1 | 1 | 0 |
T2&T5 | 1 | -1 | 0 | 1 |
T3&T6 | 1 | 0 | 1 | -1 |
T4&T5 | 1 | 0 | -1 | 1 |
表1
根据表1所示的开路故障检测表,可以实现开路故障相与故障管的定位检测。表1中,“&”表示“与”,即故障同时发生。比如:若检测结果出现Fρ=1,Δa=0,Δb=-1,Δc=0,则T4发生开路故障;如果Fρ=1,Δa=1,Δb=-1,Δc=0,则T1和T4同时开路故障。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种逆变器在线开路故障检测方法,其特征在于,包括以下步骤:
S1:采用传感器采集逆变器的两相电流,分别记为ia(k)、ib(k),k表示传感器的采样时刻;
S2:记第三相电流为ic(k),根据以下公式计算采样时刻k时三相电流的幅值Im和相位ωt:
其中,ω表示电流角频率,t表示时间;
S3:采用步骤S2计算得到的采样时刻k时三相电流的幅值Im和相位ωt计算得到采样时刻k时三相电流绝对值之和ρ(k)的标幺值fρ(k):
其中,n是满足式的正整数;
S4:根据需要设定检测阈值s1,若fρ(k)>s1,则说明采样时刻k时逆变器未发生开路故障,不作任何操作,若fρ≤s1,则逆变器开路故障发生,进入步骤S104进行故障定位;
S5:根据以下公式对采样时刻k的三相电流进行标幺化处理,得到标幺化三相电流iph(k):
其中,h=a,b,c,当h=a,j=0,当h=b,j=1,当h=c,j=2;
S6:根据需要设定阈值s2,采用以下公式对标幺化电流iph(k)进行过零置“一”处理,得到标幺化电流iph(k)对应的标志参数σh(k):
S7:采用以下公式计算采样时刻k的故障定位标识δh(k):
其中,L表示一个电流基波周期的采样点数;
S8:根据需要设定检测阈值s3,若δh(k)≤-s3,则逆变器h相发生上管开路故障;若|δh(k)|<s3,则逆变器h相未发生开路故障;若δh(k)≥s3,则逆变器h相发生下管开路故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910281011.4A CN110007183B (zh) | 2019-04-09 | 2019-04-09 | 逆变器在线开路故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910281011.4A CN110007183B (zh) | 2019-04-09 | 2019-04-09 | 逆变器在线开路故障检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110007183A true CN110007183A (zh) | 2019-07-12 |
CN110007183B CN110007183B (zh) | 2020-07-31 |
Family
ID=67170535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910281011.4A Expired - Fee Related CN110007183B (zh) | 2019-04-09 | 2019-04-09 | 逆变器在线开路故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110007183B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111381188A (zh) * | 2020-03-18 | 2020-07-07 | 华中科技大学 | 一种两电平三相电压源逆变器桥臂开路故障诊断方法 |
CN113030785A (zh) * | 2021-03-17 | 2021-06-25 | 青岛大学 | 一种电机系统开路故障检测方法 |
CN113156337A (zh) * | 2020-09-02 | 2021-07-23 | 河北工业大学 | 用于vienna整流器单管开路故障在线辨识方法、装置及存储介质 |
CN115864804A (zh) * | 2022-12-13 | 2023-03-28 | 国网湖北省电力有限公司随州供电公司 | 一种T-type三相四线整流器外功率管开路故障诊断与容错控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8362732B2 (en) * | 2010-02-02 | 2013-01-29 | GM Global Technology Operations LLC | Motor phase winding fault detection method and apparatus |
CN104597368A (zh) * | 2015-01-15 | 2015-05-06 | 电子科技大学 | 一种基于电流的三相逆变器开路故障检测方法 |
CN106597195A (zh) * | 2016-12-06 | 2017-04-26 | 中南大学 | 一种机车逆变器功率管开路故障在线诊断方法 |
CN108490309A (zh) * | 2018-02-12 | 2018-09-04 | 中国电力科学研究院有限公司 | 含逆变器类型电源配电网的故障电流迭代分析方法及装置 |
CN109375040A (zh) * | 2018-12-13 | 2019-02-22 | 国网冀北电力有限公司电力科学研究院 | 一种h桥电容器组的状态监测方法、设备及系统 |
-
2019
- 2019-04-09 CN CN201910281011.4A patent/CN110007183B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8362732B2 (en) * | 2010-02-02 | 2013-01-29 | GM Global Technology Operations LLC | Motor phase winding fault detection method and apparatus |
CN104597368A (zh) * | 2015-01-15 | 2015-05-06 | 电子科技大学 | 一种基于电流的三相逆变器开路故障检测方法 |
CN106597195A (zh) * | 2016-12-06 | 2017-04-26 | 中南大学 | 一种机车逆变器功率管开路故障在线诊断方法 |
CN108490309A (zh) * | 2018-02-12 | 2018-09-04 | 中国电力科学研究院有限公司 | 含逆变器类型电源配电网的故障电流迭代分析方法及装置 |
CN109375040A (zh) * | 2018-12-13 | 2019-02-22 | 国网冀北电力有限公司电力科学研究院 | 一种h桥电容器组的状态监测方法、设备及系统 |
Non-Patent Citations (1)
Title |
---|
戴鹏: "基于标幺化均值的PWM逆变器开路故障诊断", 《变压器世界》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111381188A (zh) * | 2020-03-18 | 2020-07-07 | 华中科技大学 | 一种两电平三相电压源逆变器桥臂开路故障诊断方法 |
CN113156337A (zh) * | 2020-09-02 | 2021-07-23 | 河北工业大学 | 用于vienna整流器单管开路故障在线辨识方法、装置及存储介质 |
CN113156337B (zh) * | 2020-09-02 | 2022-06-17 | 河北工业大学 | 用于vienna整流器单管开路故障在线辨识方法、装置及存储介质 |
CN113030785A (zh) * | 2021-03-17 | 2021-06-25 | 青岛大学 | 一种电机系统开路故障检测方法 |
CN113030785B (zh) * | 2021-03-17 | 2023-08-29 | 青岛大学 | 一种电机系统开路故障检测方法 |
CN115864804A (zh) * | 2022-12-13 | 2023-03-28 | 国网湖北省电力有限公司随州供电公司 | 一种T-type三相四线整流器外功率管开路故障诊断与容错控制方法 |
CN115864804B (zh) * | 2022-12-13 | 2023-11-10 | 国网湖北省电力有限公司随州供电公司 | 一种T-type三相四线整流器外功率管开路故障诊断与容错控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110007183B (zh) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110007183A (zh) | 逆变器在线开路故障检测方法 | |
CN104730410B (zh) | 一种基于电压电流向量的配电线路断线监测方法及装置 | |
CN103823150B (zh) | 多传感器联合的汽轮发电机转子匝间短路故障诊断方法 | |
CN106872808B (zh) | 一种相序自适应的三相电压锁相环算法 | |
CN109830972B (zh) | 一种新能源场站振荡源快速识别系统及方法 | |
CN105891660B (zh) | 一种发电机定子绕组匝间短路故障的检测方法 | |
CN101344549A (zh) | 一种三相电路缺相检测的方法及其装置 | |
CN106597217A (zh) | 一种mmc‑hvdc交流侧不对称故障诊断方法 | |
CN103872679A (zh) | 微弱信号下的电网戴维南等值模型参数辨识方法 | |
CN107045081B (zh) | 一种保护装置中实时频率的测量方法 | |
CN104808112B (zh) | 基于区段瞬时功率的配网线路故障区段定位方法 | |
CN104316888A (zh) | 一种基于sv采样信号的局放监测内同步基准修正方法 | |
CN104410044B (zh) | 基于峭度和偏度的变压器励磁涌流识别方法 | |
CN109188181B (zh) | 电网电压互感器二次回路中性线虚接判断方法 | |
CN104601083B (zh) | 一种电压不平衡条件下的异步电机负序抑制器及实现方法 | |
CN102841244A (zh) | 电网电压骤变的快速检测方法 | |
CN105203865B (zh) | 分布式电源并网逆变器及电容的在线测试系统的工作方法 | |
CN105067960B (zh) | 一种基于大数据的配电网故障定位系统及其方法 | |
CN104062555B (zh) | 配电线路高阻接地故障特征谐波的辨识方法 | |
Badoud et al. | Bond graph algorithms for fault detection and isolation in wind energy conversion | |
CN102914687A (zh) | 一种精确计算电压或电流有效值的方法 | |
CN103424643A (zh) | 接地故障检测装置 | |
CN117148256A (zh) | 一种变电站带负荷的校验方法、装置、设备和存储介质 | |
CN103454505B (zh) | 基于谐波检测的断相检测方法及电路 | |
CN107769203B (zh) | 一种用于电能质量治理装置的相序自适应方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20200731 |
|
CF01 | Termination of patent right due to non-payment of annual fee |