CN102879084A - 汽轮发电机组低频振动非稳态在线预警方法 - Google Patents

汽轮发电机组低频振动非稳态在线预警方法 Download PDF

Info

Publication number
CN102879084A
CN102879084A CN2012103885544A CN201210388554A CN102879084A CN 102879084 A CN102879084 A CN 102879084A CN 2012103885544 A CN2012103885544 A CN 2012103885544A CN 201210388554 A CN201210388554 A CN 201210388554A CN 102879084 A CN102879084 A CN 102879084A
Authority
CN
China
Prior art keywords
low
frequency
frequency vibration
vibration amplitude
sequence
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
CN2012103885544A
Other languages
English (en)
Other versions
CN102879084B (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201210388554.4A priority Critical patent/CN102879084B/zh
Publication of CN102879084A publication Critical patent/CN102879084A/zh
Application granted granted Critical
Publication of CN102879084B publication Critical patent/CN102879084B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Turbines (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了旋转机械振动状态监测与故障诊断技术领域中的一种汽轮发电机组低频振动非稳态在线预警方法。包括:设定起始时刻、终止时刻和步进长度;从起始时刻到终止时刻,获取每个步进长度所处时刻的低频振动幅值序列,并由此得到从起始时刻到终止时刻的低频振动幅值序列;分别计算每个频率序号下的低频振动幅值的最大值和每个频率序号下的低频振动幅值的参数;根据每个频率序号下的低频振动幅值的最大值和每个频率序号下的低频振动幅值的参数,判断当前时刻汽轮发电机组轴系转子一侧的低频振动是否处于非稳态。本发明实现了转子低频振动是否处于非稳态的自动判定,提高了大型汽轮发电机组低频振动非稳态实时预警分析工作的效率和准确度。

Description

汽轮发电机组低频振动非稳态在线预警方法
技术领域
本发明属于旋转机械振动状态监测与故障诊断技术领域,尤其涉及一种汽轮发电机组低频振动非稳态在线预警方法。
背景技术
汽轮发电机组轴系剧烈的低频振动可导致转子与定子在密封、轴瓦、叶顶等部位发生动静碰摩、轴瓦乌金脱胎,并产生大幅交变应力,其形成的力学疲劳环境有可能诱发转子裂纹,对机组安全性以及可靠性有潜在危害,必须消除和避免机组轴系低频振动故障。大型汽轮发电机组转子的低频振动故障,已经成为影响大型机组安全运行的重大问题。因此,一旦发现低频振动,应立即采取措施。
产生低频振动的原因非常复杂,与机组的设计、安装、工作环境、操作运行、设备状况等很多因素有关。
判断机组转子低频振动是否处于非稳态,通常由具有一定现场运行经验及专业知识技能的专业人员完成,由此带来分析结果客观性较差、分析过程耗费时间资源、人力过高等问题,并且无法做到低频振动非稳态的实时自动在线监测、分析及判别。因此,提出一种大型汽轮发电机组低频振动非稳态快速预警方法就显得十分重要。
发明内容
本发明的目的在于,提出一种汽轮发电机组低频振动非稳态在线预警方法,用以解决现有的汽轮发电机组低频振动非稳态判定过程工作效率和准确度低下的问题。
为了实现上述目的,本发明提供的技术方案是,一种汽轮发电机组低频振动非稳态在线预警方法,其特征在于,所述方法包括:
步骤1:设定起始时刻TM、终止时刻TN和步进长度t;
步骤2:在时刻TM,获取该时刻低频振动幅值序列
Figure BDA00002249655700021
其中,i为频率序号,i=1,2,...,l,l为设定值;
步骤3:存储该时刻低频振动幅值序列
Figure BDA00002249655700022
判断是否达到终止时刻TN,如果已经达到终止时刻TN,则执行步骤4;否则,令TM=TM+t,返回步骤2;
步骤4:将每个时刻存储的低频振动幅值序列
Figure BDA00002249655700023
中的数据按照存储时间的先后顺序排序,同一时刻的低频振动幅值序列
Figure BDA00002249655700024
中的数据按照下标i由小到大的顺序排序,得到从起始时刻TM到终止时刻TM的低频振动幅值序列
Figure BDA00002249655700025
其中, j = 1,2 , . . . , T N - T M t ;
步骤5:分别计算低频振动幅值序列
Figure BDA00002249655700027
的每个频率序号i下的低频振动幅值的最大值
Figure BDA00002249655700028
和低频振动幅值序列
Figure BDA00002249655700029
的每个频率序号i下的低频振动幅值的参数
Figure BDA000022496557000210
步骤6:根据低频振动幅值序列
Figure BDA000022496557000211
的每个频率序号i下的低频振动幅值的最大值
Figure BDA000022496557000212
和低频振动幅值序列
Figure BDA000022496557000213
的每个频率序号i下的低频振动幅值的参数
Figure BDA000022496557000214
判断当前时刻汽轮发电机组轴系转子一侧的低频振动是否处于非稳态。
所述步骤2包括:
步骤21:采集机组转子一侧的轴相对振动信号、转子的转速信号以及键相信号;
步骤22:根据采集的信号,利用快速傅立叶变换频谱分析方法,计算该时刻从低频到高频的振动频率所对应的振动幅值序列;
步骤23:从所述振动幅值序列中,截取所有小于机组工作转速频率的振动频率所对应的振动幅值,得到该时刻低频振动幅值序列
Figure BDA00002249655700031
其中,i为频率序号,i=1,2,...,l,l为设定值。
所述计算低频振动幅值序列
Figure BDA00002249655700032
的每个频率序号i下的低频振动幅值的最大值
Figure BDA00002249655700033
利用公式
Figure BDA00002249655700034
其中,i=1,2,...,l,l为设定值,j=1,2,...,m, m = T N - T M t .
所述计算低频振动幅值序列
Figure BDA00002249655700036
的每个频率序号i下的低频振动幅值的参数
Figure BDA00002249655700037
利用公式 β i lf = ( S i + 0.5 - m ( m - 1 ) / 4 ) / m ( 2 m 2 + 3 m - 5 ) / 72 ; 其中,Si是频率序号为i的低频振动幅值序列
Figure BDA00002249655700039
的逆序数,i=1,2,...,l,l为设定值,j=1,2,...,m,
Figure BDA000022496557000310
逆序数是指一个数据序列中逆序的总数;逆序是指在一个数据序列中,一对数的前后位置与大小顺序相反,即前面的数大于后面的数。
所述步骤6具体是,当存在频率序号i下的低频振动幅值的参数
Figure BDA000022496557000311
大于第一设定值
Figure BDA000022496557000312
并且存在频率序号i下的低频振动幅值的最大值大于第二设定值
Figure BDA000022496557000314
时,则当前时刻汽轮发电机组轴系转子一侧的低频振动处于非稳态;否则,当前时刻汽轮发电机组轴系转子一侧的低频振动不处于非稳态。
所述第一设定值
Figure BDA000022496557000315
所述第二设定值
Figure BDA00002249655700041
本发明提供的方法,对机组运行中转子振动数据进行实时自动在线监测、分析及判别,判定转子低频振动是否处于非稳态,提高了大型汽轮发电机组低频振动非稳态实时预警分析工作的效率和准确度。
附图说明
图1是汽轮发电机组低频振动非稳态在线预警方法流程图;
图2是汽轮发电机组低频振动非稳态在线预警示意图。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
实施例
图1是汽轮发电机组低频振动非稳态在线预警方法流程图。如图1所示,本发明提供的一种汽轮发电机组低频振动非稳态在线预警方法包括:
步骤1:设定起始时刻TM=0秒、终止时刻TN=10秒和步进长度t=0.1秒。
步骤2:在时刻TM,获取该时刻低频振动幅值序列
Figure BDA00002249655700042
其中,i为频率序号,i=1,2,...,l,l为设定值。
获取该时刻低频振动幅值序列具体包括:
步骤21:采集机组转子一侧的轴相对振动信号、转子的转速信号以及键相信号。
图2是汽轮发电机组低频振动非稳态在线预警示意图。如图2所示,转子轴相对振动数据、转子的转速信号以及键相信号可以从配置汽轮发电机组的监视仪表(TSI)获得,机组功率数据信号可以从配置汽轮发电机组的分布式控制系统(DCS)获得。本实施例中,转子轴相对振动数据、转子的转速信号以及键相信号是从配置汽轮发电机组的监视仪表(TSI)获得,机组功率数据信号是从配置汽轮发电机组的分布式控制系统(DCS)获得。数据采集卡插入工业用微型计算机(IPC)提供的插槽内。根据数据采集卡的要求,数据采集调理设备处理来自汽轮发电机组监视仪表(TSI)的轴相对振动信号、转子的转速信号、键相信号,经过处理后的轴相对振动信号、转子的转速信号、键相信号输入IPC内的高速振动数据采集卡。振动数据采集卡每一通道技术参数为50ks/s,24bit。
根据本发明提供的方法设计具体的汽轮发电机组低频振动非稳态实时预警程序,将程序安装在工业用微型计算机(IPC)内。汽轮发电机组低频振动非稳态实时预警程序中的一次诊断循环过程,包括诊断方法中涉及的数据实时采集计算存储、实时判别、低频振动参量实时计算及低频振动非稳态判定等一系列计算分析环节。
利用汽轮发电机组低频振动非稳态实时预警程序监测分析高压转子的低频振动是否处于非稳态。首先,工业用微型计算机(IPC)中的分析程序通过采用高速振动数据采集卡实时采集汽轮发电机组高压转子A侧支持轴承附近测得的轴相对振动数据、转子的转速信号以及键相信号。振动数据采集卡每一通道技术参数为50ks/s,24bit。
步骤22:根据采集的信号,利用快速傅立叶变换频谱分析方法,计算该时刻从低频到高频的振动频率所对应的振动幅值序列。
针对机组高压转子一侧的轴相对振动数据,利用FFT(快速傅立叶变换)频谱分析方法,计算得到当前时刻从低频到高频的不同振动频率所对应的振动幅值数据序列(幅值单位为μm,即微米)。
步骤23:从所述振动幅值序列中,截取得到频率小于机组工作转速对应频率(50Hz)的当前时刻低频振动幅值序列
Figure BDA00002249655700061
(i=1,2,3,...,l)。设定振动数据采集频率及采集数据量,可使得低频振动幅值序列数据个数l=98。i为不同频率对应的频率号,这些不同的频率都是小于机组工作转速对应频率(50Hz)的频率。
事实上,上述获取该时刻低频振动幅值序列
Figure BDA00002249655700062
的过程已经是现有技术,在本发明中不再赘述。
步骤3:存储该时刻低频振动幅值序列
Figure BDA00002249655700063
判断是否达到终止时刻TN,如果已经达到终止时刻TN,则执行步骤4;否则,令TM=TM+t,返回步骤2。
步骤4:将每个时刻存储的低频振动幅值序列
Figure BDA00002249655700064
中的数据按照存储时间的先后顺序排序,同一时刻的低频振动幅值序列
Figure BDA00002249655700065
中的数据按照下标i由小到大的顺序排序,得到从起始时刻TM到终止时刻TN的低频振动幅值序列
Figure BDA00002249655700066
其中, j = 1,2 , . . . , T N - T M t = 10 - 0 0.1 = 100 .
步骤5:分别计算低频振动幅值序列
Figure BDA00002249655700068
的每个频率序号i下的低频振动幅值的最大值
Figure BDA00002249655700069
和低频振动幅值序列
Figure BDA000022496557000610
的每个频率序号i下的低频振动幅值的参数
计算低频振动幅值序列
Figure BDA000022496557000612
的每个频率序号i下的低频振动幅值的最大值
Figure BDA000022496557000613
利用公式
Figure BDA000022496557000614
其中,i=1,2,...,98,j=1,2,...,100。
计算低频振动幅值序列
Figure BDA000022496557000615
的每个频率序号i下的低频振动幅值的参数
Figure BDA000022496557000616
利用公式 β i lf = ( S i + 0.5 - m ( m - 1 ) / 4 ) / m ( 2 m 2 + 3 m - 5 ) / 72 ; 其中,Si是频率序号为i的低频振动幅值序列
Figure BDA000022496557000618
的逆序数,i=1,2,...,98,j=1,2,...,100。逆序数是指一个数据序列中逆序的总数;逆序是指在一个数据序列中,一对数的前后位置与大小顺序相反,即前面的数大于后面的数。
步骤6:根据低频振动幅值序列
Figure BDA00002249655700071
的每个频率序号i下的低频振动幅值的最大值
Figure BDA00002249655700072
和低频振动幅值序列
Figure BDA00002249655700073
的每个频率序号i下的低频振动幅值的参数
Figure BDA00002249655700074
判断当前时刻汽轮发电机组轴系转子一侧的低频振动是否处于非稳态。具体是:当存在某个频率序号i(i=1,2,...,98),使得在该频率号下的低频振动幅值的参数
Figure BDA00002249655700075
大于第一设定值
Figure BDA00002249655700076
并且在该频率号下的低频振动幅值的最大值
Figure BDA00002249655700077
大于第二设定值
Figure BDA00002249655700078
时,则当前时刻汽轮发电机组轴系转子一侧的低频振动处于非稳态;否则,当前时刻汽轮发电机组轴系转子一侧的低频振动不处于非稳态。
比如,假设第一设定值
Figure BDA00002249655700079
第二设定值
Figure BDA000022496557000710
(微米),如果存在频率序号i=48,该频率序号下的低频振动幅值的参数
Figure BDA000022496557000711
该频率序号下的低频振动幅值的最大值
Figure BDA000022496557000712
使得在该频率序号下的低频振动幅值的参数
Figure BDA000022496557000713
和在该频率序号下的低频振动幅值的最大值
Figure BDA000022496557000714
满足条件
Figure BDA000022496557000715
Figure BDA000022496557000716
则当前时刻汽轮发电机组轴系转子一侧的低频振动处于非稳态。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (7)

1.一种汽轮发电机组低频振动非稳态在线预警方法,其特征在于,所述方法包括:
步骤1:设定起始时刻TM、终止时刻TN和步进长度t;
步骤2:在时刻TM,获取该时刻低频振动幅值序列
Figure FDA00002249655600011
其中,i为频率序号,i=1,2,...,l,l为设定值;
步骤3:存储该时刻低频振动幅值序列判断是否达到终止时刻TN,如果已经达到终止时刻TN,则执行步骤4;否则,令TM=TM+t,返回步骤2;
步骤4:将每个时刻存储的低频振动幅值序列
Figure FDA00002249655600013
中的数据按照存储时间的先后顺序排序,同一时刻的低频振动幅值序列
Figure FDA00002249655600014
中的数据按照下标i由小到大的顺序排序,得到从起始时刻TM到终止时刻TN的低频振动幅值序列
Figure FDA00002249655600015
其中, j = 1,2 , . . . , T N - T M t ;
步骤5:分别计算低频振动幅值序列
Figure FDA00002249655600017
的每个频率序号i下的低频振动幅值的最大值
Figure FDA00002249655600018
和低频振动幅值序列
Figure FDA00002249655600019
的每个频率序号i下的低频振动幅值的参数
Figure FDA000022496556000110
步骤6:根据低频振动幅值序列
Figure FDA000022496556000111
的每个频率序号i下的低频振动幅值的最大值
Figure FDA000022496556000112
和低频振动幅值序列
Figure FDA000022496556000113
的每个频率序号i下的低频振动幅值的参数
Figure FDA000022496556000114
判断当前时刻汽轮发电机组轴系转子一侧的低频振动是否处于非稳态。
2.根据权利要求1所述的方法,其特征在于,所述步骤2包括:
步骤21:采集机组转子一侧的轴相对振动信号、转子的转速信号以及键相信号;
步骤22:根据采集的信号,利用快速傅立叶变换频谱分析方法,计算该时刻从低频到高频的振动频率所对应的振动幅值序列;
步骤23:从所述振动幅值序列中,截取所有小于机组工作转速频率的振动频率所对应的振动幅值,得到该时刻低频振动幅值序列
Figure FDA00002249655600021
其中,i为频率序号,i=1,2,...,l,l为设定值。
3.根据权利要求2所述的方法,其特征在于,所述计算低频振动幅值序列
Figure FDA00002249655600022
的每个频率序号i下的低频振动幅值的最大值
Figure FDA00002249655600023
利用公式
Figure FDA00002249655600024
其中,i=1,2,...,l,l为设定值,j=1,2,...,m,
Figure FDA00002249655600025
4.根据权利要求3所述的方法,其特征在于,所述计算低频振动幅值序列
Figure FDA00002249655600026
的每个频率序号i下的低频振动幅值的参数
Figure FDA00002249655600027
利用公式 β i lf = ( S i + 0.5 - m ( m - 1 ) / 4 ) / m ( 2 m 2 + 3 m - 5 ) / 72 ; 其中,Si是频率序号为i的低频振动幅值序列
Figure FDA00002249655600029
的逆序数,i=1,2,...,l,l为设定值,j=1,2,...,m,
Figure FDA000022496556000210
逆序数是指一个数据序列中逆序的总数;逆序是指在一个数据序列中,一对数的前后位置与大小顺序相反,即前面的数大于后面的数。
5.根据权利要求4所述的方法,其特征在于,所述步骤6具体是,当存在频率序号i下的低频振动幅值的参数
Figure FDA000022496556000211
大于第一设定值
Figure FDA000022496556000212
并且存在频率序号i下的低频振动幅值的最大值
Figure FDA000022496556000213
大于第二设定值
Figure FDA000022496556000214
时,则当前时刻汽轮发电机组轴系转子一侧的低频振动处于非稳态;否则,当前时刻汽轮发电机组轴系转子一侧的低频振动不处于非稳态。
6.根据权利要求5所述的方法,其特征在于,所述第一设定值
7.根据权利要求5所述的方法,其特征在于,所述第二设定值
Figure FDA00002249655600032
CN201210388554.4A 2012-10-12 2012-10-12 汽轮发电机组低频振动非稳态在线预警方法 Expired - Fee Related CN102879084B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210388554.4A CN102879084B (zh) 2012-10-12 2012-10-12 汽轮发电机组低频振动非稳态在线预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210388554.4A CN102879084B (zh) 2012-10-12 2012-10-12 汽轮发电机组低频振动非稳态在线预警方法

Publications (2)

Publication Number Publication Date
CN102879084A true CN102879084A (zh) 2013-01-16
CN102879084B CN102879084B (zh) 2014-08-20

Family

ID=47480497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210388554.4A Expired - Fee Related CN102879084B (zh) 2012-10-12 2012-10-12 汽轮发电机组低频振动非稳态在线预警方法

Country Status (1)

Country Link
CN (1) CN102879084B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103335708A (zh) * 2013-06-13 2013-10-02 华北电力大学 汽轮发电机组低频振动实时预警方法
CN108885473A (zh) * 2016-03-30 2018-11-23 深圳市大疆创新科技有限公司 用于控制电机的方法和系统
CN114542995A (zh) * 2022-03-03 2022-05-27 西安热工研究院有限公司 一种管道低频振动非稳态在线预警诊断系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072764A (zh) * 2010-11-24 2011-05-25 华北电力大学 汽轮发电机组低频振动单峰实时分析方法
CN102087140A (zh) * 2010-11-24 2011-06-08 华北电力大学 汽轮发电机组低频振动主峰频率平稳性分析方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072764A (zh) * 2010-11-24 2011-05-25 华北电力大学 汽轮发电机组低频振动单峰实时分析方法
CN102087140A (zh) * 2010-11-24 2011-06-08 华北电力大学 汽轮发电机组低频振动主峰频率平稳性分析方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103335708A (zh) * 2013-06-13 2013-10-02 华北电力大学 汽轮发电机组低频振动实时预警方法
CN103335708B (zh) * 2013-06-13 2014-12-03 华北电力大学 汽轮发电机组低频振动实时预警方法
CN108885473A (zh) * 2016-03-30 2018-11-23 深圳市大疆创新科技有限公司 用于控制电机的方法和系统
US11108352B2 (en) 2016-03-30 2021-08-31 SZ DJI Technology Co., Ltd. Method and system for controlling a motor
CN114542995A (zh) * 2022-03-03 2022-05-27 西安热工研究院有限公司 一种管道低频振动非稳态在线预警诊断系统
CN114542995B (zh) * 2022-03-03 2024-03-15 西安热工研究院有限公司 一种管道低频振动非稳态在线预警诊断系统

Also Published As

Publication number Publication date
CN102879084B (zh) 2014-08-20

Similar Documents

Publication Publication Date Title
CN101738293B (zh) 汽轮发电机组转子原始质量不平衡故障实时诊断方法
CN101451898B (zh) 汽轮发电机组转子热弯曲不平衡故障实时诊断方法
CN105569932A (zh) 风电机组动态不平衡在线检测和故障辨识方法及系统
CN102095491B (zh) 汽轮发电机组低频振动突变性实时分析方法
CN102072764B (zh) 汽轮发电机组低频振动单峰实时分析方法
CN101532911A (zh) 大型汽轮发电机组转子裂纹故障实时诊断方法
CN102087139B (zh) 汽轮发电机组低频振动频率成份实时分析方法
CN102175306B (zh) 汽轮发电机组油膜振荡故障实时辨识方法
CN102087140B (zh) 汽轮发电机组低频振动主峰频率平稳性分析方法
CN104568435A (zh) 一种基于时频能量谱提取风机状态信号调制成分的方法
CN102879084B (zh) 汽轮发电机组低频振动非稳态在线预警方法
CN102096771B (zh) 汽轮发电机组轴承瓦温变化趋势实时分析方法
CN102175409B (zh) 汽轮发电机组油膜涡动故障实时辨识方法
CN102175408B (zh) 汽轮发电机组轴承座刚度实时辨识方法
CN102175307B (zh) 汽轮发电机组低频振动谱阵实时定量分析方法
CN102692303B (zh) 汽轮发电机组汽流激振故障高效识别方法
CN102879085B (zh) 汽轮发电机组低频振动非稳态实时预警方法
CN102012263B (zh) 汽轮机组转子振动同相分量平稳性实时辨识方法
CN102095564B (zh) 汽轮发电机组波动型碰摩故障实时辨识方法
CN102865921B (zh) 汽轮发电机组低频振动非稳态快速预警方法
CN102022143B (zh) 汽轮发电机组低频振动平稳性实时分析方法
CN102103037A (zh) 汽轮发电机组低频振动与功率递增相关性分析方法
CN102645336B (zh) 汽轮发电机组汽流激振故障实时辨识方法
CN102680243B (zh) 汽轮发电机组汽流激振故障在线判别方法
CN102183349B (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
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: 20140820

Termination date: 20151012

EXPY Termination of patent right or utility model