CN105425079A - 基于实时平均功率频域分析的洗衣机运行非侵入辨识方法 - Google Patents

基于实时平均功率频域分析的洗衣机运行非侵入辨识方法 Download PDF

Info

Publication number
CN105425079A
CN105425079A CN201510946076.8A CN201510946076A CN105425079A CN 105425079 A CN105425079 A CN 105425079A CN 201510946076 A CN201510946076 A CN 201510946076A CN 105425079 A CN105425079 A CN 105425079A
Authority
CN
China
Prior art keywords
real
sequence
frequency
washing machine
average 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.)
Granted
Application number
CN201510946076.8A
Other languages
English (en)
Other versions
CN105425079B (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.)
Wuxi Hengtong Electric Appliance Co., Ltd.
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201510946076.8A priority Critical patent/CN105425079B/zh
Publication of CN105425079A publication Critical patent/CN105425079A/zh
Application granted granted Critical
Publication of CN105425079B publication Critical patent/CN105425079B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,包括:在一定采样频率范围内,对总电源进线的电压和电流进行采样,形成电压采样序列u和电流采样序列i;在一定计算时间窗口内,计算总电源进线处的实时平均功率序列P;选取实时平均功率序列P的一组样本量进行快速傅里叶变换,得到该实时平均功率序列P的幅值结果序列和频率结果序列;根据实时平均功率序列P的幅值结果序列和频率结果序列判断是否有洗衣机启动。本发明提供了一种全新的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,该方法能够准确感知洗衣机的运行,为实现洗衣机的非侵入辨识提供了技术支撑。

Description

基于实时平均功率频域分析的洗衣机运行非侵入辨识方法
技术领域
本发明属于智能用电技术领域,尤其涉及一种基于实时平均功率频域分析的洗衣机运行非侵入辨识方法。
背景技术
居民电力负荷监测分解技术是一门新兴的智能电网基础支撑技术,与目前智能电表仅量测用户总功率不同,它以监测并分解出居民户内所有电器的启动时间、工作状态、能耗情况为目标,从而实现更加可靠、精确的电能量管理。电力负荷监测分解技术使用户的电费清单像电话费清单一样,各类家用电器的用电量一目了然,从而使用户及时了解自己的用电情况,为合理分配各个电器的用电时间及相应的用电量提供参考,最终能够有效减少电费支出和电能浪费。Google统计数据显示,如果家庭用户能够及时了解住宅电器的详细用电信息,就能使每月电费开支下降5%~15%。如果全美国有一半家庭每个月节省这么多开支,减少的碳排放量相当于减少800万辆汽车的使用。对于工业用户而言,其负荷投切安排一般是比较固定的,只需分时计量即可,对负荷分解的需求较少,本项目的主要研究对象是住宅用电负荷。
目前,居民电力负荷监测分解技术主要分为侵入式监测分解(IntrusiveLoadMonitoringanddecomposition,ILMD)和非侵入式监测分解(Non-intrusiveLoadMonitoringanddecomposition,NILMD)两大类:
(1)侵入式负荷监测分解技术(ILMD):侵入式负荷监测将带有数字通信功能的传感器安装在每个电器与电网的接口,可以准确监测每个负荷的运行状态和功率消耗。但大量安装监测传感器造成建设和维护的成本较高,最重要的是侵入式负荷监测需要进入居民家中进行安装调试,容易造成用户抵制心理。
(2)非侵入式负荷监测分解技术(NILMD):仅在用户入口处安装一个传感器,通过采集和分析入口总电流、电压等信息来判断户内每个或每类电器的用电功率和工作状态(例如,空调具有制冷、制热、待机等不同工作状态),从而得出居民的用电规律。和侵入式负荷分解相比,由于只需要安装一个监测传感器,非侵入负荷分解方案的建设成本和后期维护难度都大幅降低;另外,传感器安装位置可以选择在用户电表箱处,完全不会侵入居民户内进行施工。可以认为,NILMD以分解算法代替ILMD系统的传感器网络,具有简单、经济、可靠、数据完整和易于迅速推广应用等优势,有望发展成为高级量测体系(AMI)中新一代核心技术(成熟后,NILMD算法也可以融合到智能电表的芯片内),支持需求侧管理、定制电力等智能用电的高级功能,也适用于临时性的负荷用电细节监测与调查。
洗衣机利用电动机以及相应的机械结构带动衣物上下左右翻转,由于家庭衣物量有限,家用洗衣机的稳态有功功率一般在200W至1000W之间。洗衣机洗衣过程基本由设备自动控制,由同一电机控制洗衣桶以一定频率(约为0.05Hz)实现正反转交替运行。
综上所述,NILMD技术已经逐渐成为一个研究热点,相关技术的突破和产业化对全社会的节能减排具有重要意义。目前,NILMD技术的研究还停留在理论研究阶段,多态电器尤其是洗衣机的分解辨识方法等关键技术还有待突破。
因此,亟待解决上述问题。
发明内容
发明目的:本发明的目的是提供一种能够准确感知洗衣机的运行的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法。
技术方案:为实现以上目的,本发明公开了一种基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,其特征在于:该辨识方法包括如下步骤:
(1)在一定采样频率范围内,对总电源进线的电压和电流进行采样,形成电压采样序列u和电流采样序列i;
(2)在一定计算时间窗口内,计算总电源进线处的实时平均功率序列P;
(3)选取实时平均功率序列P的一组样本量进行快速傅里叶变换,得到该实时平均功率序列P的幅频特性;
(4)根据步骤(3)中实时平均功率序列P的幅频特性判断是否有洗衣机启动,如未启动,则返回步骤(3)。
其中,所述步骤(1)中的采样频率范围为f=0.5kHz~2kHz。
优选的,所述步骤(2)中所述实时平均功率序列P的计算公式为
P = 1 m N Σ k k + m N u * i
其中,m为实时平均功率序列P的计算时间窗口所含工频周期数目,取m=5个工频周期,k为采样点编号,N为一个工频周期包含的采样点数目。
进一步,所述步骤(3)中所述实时平均功率序列P的幅频特性包括幅值结果序列和频率结果序列,其计算公式为
Fp[j]=FFT(P),fp[j]=(j-1)*50/(m*Ns)
其中,j=1,2,…,Ns,FFT表示快速傅里叶变换运算,Fp[j]为幅值结果序列,fp[j]为频率结果序列,Ns为参与本次快速傅里叶变换的P的样本数量,Ns取600~1200个。
其中,所述步骤(4)中判断方法为:扫描幅频特性中的幅值结果序列和频率结果序列,如果在0.01Hz~0.1Hz频率范围内存在幅值大于50W的频谱分量,即可判断当次计算时间窗口内有洗衣机启动,否则返回步骤(3)对实时平均功率序列P的下一组样本量进行快速傅里叶变换。
有益效果:与现有技术相比,本发明具有以下显著优点:因多状态居民电器尤其是洗衣机的负荷识别一直是非侵入辨识的难点,且现有的洗衣机运行非侵入辨识方法主要围绕洗衣机的启功电流、有功功率、无功功率等负荷特性进行探索研究,但实质性的能准确辨识出洗衣机运行状态的方法尚且没有。本发明提供了的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,解决了上述难题,能够准确感知洗衣机的运行,为实现洗衣机的非侵入辨识提供了技术支撑。
附图说明
图1为本发明基于实时平均功率频域分析的洗衣机运行非侵入辨识方法的流程示意图;
图2为基于实时平均功率频域分析的洗衣机运行非侵入辨识方法中实时平均功率的计算结果图;
图3为基于实时平均功率频域分析的洗衣机运行非侵入辨识方法中实时平均功率幅频特性图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1、图2和图3所示,本发明公开了一种基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,具体的流程步骤如下:
(1)取采样频率f=1kHz,对总电源进线的电压和电流进行采样,形成电压采样序列u和电流采样序列i。
(2)取计算时间窗口m=5个工频周期,计算总电源进线处的实时平均功率序列P,P的计算公式为式中,k为采样点编号,N为一个工频周期包含的采样点数目;
从图2为洗衣机一次完整洗涤过程对应的实时平均功率的计算结果图,可以看出洗衣机的功率波动性大,存在较多功率波动“尖峰”;此外,该洗衣机的最大有功功率为200W左右,功率波动一次的时间间隔A为20s左右,波动频率约为0.05Hz,该频率对应洗衣机的正反转频率。
(3)选取实时平均功率序列P的一组样本量进行快速傅里叶变换,得到实时平均功率序列P的幅频特性,具体计算公式为Fp[j]=FFT(P),fp[j]=(j-1)*50/(m*Ns),其中j=1,2,…,Ns,其中FFT表示快速傅里叶变换运算,Fp[j]为快速傅里叶变换的幅值结果序列,fp[j]为快速傅里叶变换的频率结果序列,Ns为参与本次快速傅里叶变换的P的样本数量,算例中Ns取1000个。
(4)根据步骤(3)的实时平均功率序列P的幅频特性判断是否有洗衣机启动。扫描幅频特性中的幅值结果序列Fp[j]和频率结果序列fp[j],如果在0.01Hz~0.1Hz频率范围内存在幅值大于50W的频谱分量,即可判断本次处理的时间窗口内有洗衣机启动,否则返回步骤(3)对实时平均功率序列P的下一组样本量进行快速傅里叶变换。
如图3所示的幅频特性图分析可知,在0.04Hz~0.05Hz频率范围内存在功率幅值约为70W的频谱分量,故可判断本次处理的时间窗口内有洗衣机启动。

Claims (5)

1.一种基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,其特征在于:该辨识方法包括如下步骤:
(1)在一定采样频率范围内,对总电源进线的电压和电流进行采样,形成电压采样序列u和电流采样序列i;
(2)在一定计算时间窗口内,计算总电源进线处的实时平均功率序列P;
(3)选取实时平均功率序列P的一组样本量进行快速傅里叶变换,得到该实时平均功率序列P的幅频特性;
(4)根据步骤(3)中实时平均功率序列P的幅频特性判断是否有洗衣机启动,如未启动,则返回步骤(3)。
2.根据权利要求1所述的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,其特征在于:所述步骤(1)中的采样频率范围为f=0.5kHz~2kHz。
3.根据权利要求1所述的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,其特征在于:所述步骤(2)中所述实时平均功率序列P的计算公式为
P = 1 m N Σ k k + m N u * i
其中,m为实时平均功率序列P的计算时间窗口所含工频周期数目,取m=5个工频周期,k为采样点编号,N为一个工频周期包含的采样点数目。
4.根据权利要求1所述的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,其特征在于:所述步骤(3)中所述实时平均功率序列P的幅频特性包括幅值结果序列和频率结果序列,其计算公式为
Fp[j]=FFT(P),fp[j]=(j-1)*50/(m*Ns)
其中,j=1,2,…,Ns,FFT表示快速傅里叶变换运算,Fp[j]为幅值结果序列,fp[j]为频率结果序列,Ns为参与本次快速傅里叶变换的P的样本数量,Ns取600~1200个。
5.根据权利要求1所述的基于实时平均功率频域分析的洗衣机运行非侵入辨识方法,其特征在于:所述步骤(4)中判断方法为:扫描幅频特性中的幅值结果序列和频率结果序列,如果在0.01Hz~0.1Hz频率范围内存在幅值大于50W的频谱分量,即可判断当次计算时间窗口内有洗衣机启动,否则返回步骤(3)对实时平均功率序列P的下一组样本量进行快速傅里叶变换。
CN201510946076.8A 2015-12-17 2015-12-17 基于实时平均功率频域分析的洗衣机运行非侵入辨识方法 Expired - Fee Related CN105425079B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510946076.8A CN105425079B (zh) 2015-12-17 2015-12-17 基于实时平均功率频域分析的洗衣机运行非侵入辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510946076.8A CN105425079B (zh) 2015-12-17 2015-12-17 基于实时平均功率频域分析的洗衣机运行非侵入辨识方法

Publications (2)

Publication Number Publication Date
CN105425079A true CN105425079A (zh) 2016-03-23
CN105425079B CN105425079B (zh) 2018-07-13

Family

ID=55503411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510946076.8A Expired - Fee Related CN105425079B (zh) 2015-12-17 2015-12-17 基于实时平均功率频域分析的洗衣机运行非侵入辨识方法

Country Status (1)

Country Link
CN (1) CN105425079B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105911342A (zh) * 2016-04-15 2016-08-31 东南大学 基于功率持续时间特性的电饭煲非侵入辨识方法
CN106532719A (zh) * 2016-11-28 2017-03-22 国网江苏省电力公司苏州供电公司 基于电流二次谐波与无功的非变频空调非侵入式辨识方法
CN107247201A (zh) * 2017-06-09 2017-10-13 东南大学 基于功率及时间特性的电烤箱非侵入辨识方法
CN107356827A (zh) * 2017-06-09 2017-11-17 东南大学 一种基于有功功率波动性的洗衣机运行非侵入辨识方法
CN107390043A (zh) * 2017-06-09 2017-11-24 东南大学 一种基于有功功率波动性和周期性的洗衣机非侵入辨识方法
CN107478926A (zh) * 2017-07-14 2017-12-15 宁波三星医疗电气股份有限公司 一种用电监视装置及方法
CN107561354A (zh) * 2017-08-17 2018-01-09 江苏智臻能源科技有限公司 基于功率持续时间特性的多电器非侵入辨识方法
CN108267637A (zh) * 2017-12-29 2018-07-10 国网天津市电力公司电力科学研究院 一种基于混合频谱特性的洗衣机非侵入式运行检测方法
CN111983303A (zh) * 2020-07-09 2020-11-24 张子毅 一种基于占空比分解的洗衣机非侵入电量估算方法
CN112034245A (zh) * 2020-11-03 2020-12-04 天津滨电电力工程有限公司 一种用于嵌入式终端的家用洗衣机辨识方法
WO2021109523A1 (zh) * 2019-12-02 2021-06-10 江苏智臻能源科技有限公司 一种基于功率谐波特征的ptc型电采暖非侵入辨识方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738582A (zh) * 2008-11-14 2010-06-16 中国科学院沈阳自动化研究所 非侵入式电机效率在线检测装置及其控制方法
CN103217603A (zh) * 2013-03-22 2013-07-24 重庆大学 一种非侵入式家用电器用电状况在线监测的识别方法
CN103998942A (zh) * 2010-12-13 2014-08-20 美国弗劳恩霍夫股份公司 用于非侵入式负载监测的方法和系统
EP2801792A2 (en) * 2013-05-08 2014-11-12 Samsung Electronics Co., Ltd. Non-intrusive load monitoring apparatus and method
CN104483575A (zh) * 2014-12-22 2015-04-01 天津求实智源科技有限公司 用于非侵入式电力监测的自适应负荷事件检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738582A (zh) * 2008-11-14 2010-06-16 中国科学院沈阳自动化研究所 非侵入式电机效率在线检测装置及其控制方法
CN103998942A (zh) * 2010-12-13 2014-08-20 美国弗劳恩霍夫股份公司 用于非侵入式负载监测的方法和系统
CN103217603A (zh) * 2013-03-22 2013-07-24 重庆大学 一种非侵入式家用电器用电状况在线监测的识别方法
EP2801792A2 (en) * 2013-05-08 2014-11-12 Samsung Electronics Co., Ltd. Non-intrusive load monitoring apparatus and method
CN104483575A (zh) * 2014-12-22 2015-04-01 天津求实智源科技有限公司 用于非侵入式电力监测的自适应负荷事件检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陆明琼: "一种基于实时平均功率和实时电流的非变频空调辨识判据", 《百度文库HTTPS://WENKU.BAIDU.COM/VIEW/BE835FA8B90D6C85ED3AC62E.HTML? FROM》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105911342A (zh) * 2016-04-15 2016-08-31 东南大学 基于功率持续时间特性的电饭煲非侵入辨识方法
CN105911342B (zh) * 2016-04-15 2019-03-05 东南大学 基于功率持续时间特性的电饭煲非侵入辨识方法
CN106532719A (zh) * 2016-11-28 2017-03-22 国网江苏省电力公司苏州供电公司 基于电流二次谐波与无功的非变频空调非侵入式辨识方法
CN106532719B (zh) * 2016-11-28 2019-06-18 国网江苏省电力公司苏州供电公司 基于电流二次谐波与无功的非变频空调非侵入式辨识方法
CN107390043A (zh) * 2017-06-09 2017-11-24 东南大学 一种基于有功功率波动性和周期性的洗衣机非侵入辨识方法
CN107356827B (zh) * 2017-06-09 2019-10-11 东南大学 一种基于有功功率波动性的洗衣机运行非侵入辨识方法
CN107356827A (zh) * 2017-06-09 2017-11-17 东南大学 一种基于有功功率波动性的洗衣机运行非侵入辨识方法
CN107247201A (zh) * 2017-06-09 2017-10-13 东南大学 基于功率及时间特性的电烤箱非侵入辨识方法
CN107390043B (zh) * 2017-06-09 2019-07-30 东南大学 基于有功功率波动性和周期性的洗衣机非侵入辨识方法
CN107478926A (zh) * 2017-07-14 2017-12-15 宁波三星医疗电气股份有限公司 一种用电监视装置及方法
CN107478926B (zh) * 2017-07-14 2020-10-09 宁波三星医疗电气股份有限公司 一种用电监视装置及方法
CN107561354B (zh) * 2017-08-17 2019-09-10 江苏智臻能源科技有限公司 基于功率持续时间特性的多电器非侵入辨识方法
CN107561354A (zh) * 2017-08-17 2018-01-09 江苏智臻能源科技有限公司 基于功率持续时间特性的多电器非侵入辨识方法
CN108267637A (zh) * 2017-12-29 2018-07-10 国网天津市电力公司电力科学研究院 一种基于混合频谱特性的洗衣机非侵入式运行检测方法
WO2021109523A1 (zh) * 2019-12-02 2021-06-10 江苏智臻能源科技有限公司 一种基于功率谐波特征的ptc型电采暖非侵入辨识方法
CN111983303A (zh) * 2020-07-09 2020-11-24 张子毅 一种基于占空比分解的洗衣机非侵入电量估算方法
CN111983303B (zh) * 2020-07-09 2023-09-29 张子毅 一种基于占空比分解的洗衣机非侵入电量估算方法
CN112034245A (zh) * 2020-11-03 2020-12-04 天津滨电电力工程有限公司 一种用于嵌入式终端的家用洗衣机辨识方法

Also Published As

Publication number Publication date
CN105425079B (zh) 2018-07-13

Similar Documents

Publication Publication Date Title
CN105425079A (zh) 基于实时平均功率频域分析的洗衣机运行非侵入辨识方法
CN108152630B (zh) 一种基于混合判据的非侵入式ih电饭煲运行辨识方法
CN105911342B (zh) 基于功率持续时间特性的电饭煲非侵入辨识方法
CN106093630B (zh) 一种非侵入式家用电器辨识方法
CN106532719B (zh) 基于电流二次谐波与无功的非变频空调非侵入式辨识方法
CN110412347B (zh) 一种基于非侵入式负荷监测的窃电行为识别方法与装置
JP2017011991A (ja) 負荷識別データプロセッサを備えた電気ユーティリティメータ
CN107247201A (zh) 基于功率及时间特性的电烤箱非侵入辨识方法
CN105629065A (zh) 基于平均有功无功综合频谱分析的微波炉非侵入辨识方法
Liu et al. Non-intrusive energy use monitoring for a group of electrical appliances
CN108595376B (zh) 一种区分定频空调与冲击钻的非侵入负荷辨识方法
CN103116098A (zh) 一种基于互相关系数的家用电器运行状态识别方法
CN109470957B (zh) 一种基于弗雷歇距离的非侵入式电器识别方法
CN106771593A (zh) 基于混合判据的非侵入式电磁炉启动辨识方法
CN106340884A (zh) 基于混合功率判据的非侵入式非变频冰箱启动辨识方法
Zhang et al. An improved non-intrusive load monitoring method for recognition of electric vehicle battery charging load
CN105652118A (zh) 一种基于负荷瞬时能量特征的电网电能负荷监测方法
CN105182250B (zh) 一种非侵入式分布式电源出力监测装置
CN108152574A (zh) 一种基于混合判据的微波炉运行的非侵入辨识方法
CN107390043A (zh) 一种基于有功功率波动性和周期性的洗衣机非侵入辨识方法
CN109449919A (zh) 一种非侵入式用电负荷分析方法及装置
CN111025013A (zh) 一种基于功率谐波特征的ptc型电采暖非侵入辨识方法
CN103018611B (zh) 一种基于电流分解的非侵入式负载监测方法及系统
Hoyo-Montaño et al. Overview of Non-Intrusive Load Monitoring: A way to energy wise consumption
CN108267637A (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
CB03 Change of inventor or designer information

Inventor after: Zhou Gan

Inventor after: Wu Bin

Inventor after: Qin Donglei

Inventor after: Fu Wang

Inventor after: Fu Meng

Inventor after: Li Yongkun

Inventor after: Yao Xun

Inventor after: Zhang Liang

Inventor after: Zhang Xu

Inventor after: Qin Chengming

Inventor after: Gu Wei

Inventor before: Zhou Gan

Inventor before: Fu Wang

Inventor before: Fu Meng

Inventor before: Li Yongkun

Inventor before: Yao Xun

Inventor before: Zhang Liang

Inventor before: Zhang Xu

Inventor before: Qin Chengming

Inventor before: Gu Wei

CB03 Change of inventor or designer information
TA01 Transfer of patent application right

Effective date of registration: 20170616

Address after: Four pailou Nanjing Xuanwu District of Jiangsu Province, No. 2 210096

Applicant after: Southeast University

Applicant after: Wuxi Hengtong Electric Appliance Co., Ltd.

Address before: Four pailou Nanjing Xuanwu District of Jiangsu Province, No. 2 210096

Applicant before: Southeast University

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180713

Termination date: 20191217

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