CN105353268A - 一种用于输电线路分布式行波故障判断及定位方法 - Google Patents

一种用于输电线路分布式行波故障判断及定位方法 Download PDF

Info

Publication number
CN105353268A
CN105353268A CN201510650313.6A CN201510650313A CN105353268A CN 105353268 A CN105353268 A CN 105353268A CN 201510650313 A CN201510650313 A CN 201510650313A CN 105353268 A CN105353268 A CN 105353268A
Authority
CN
China
Prior art keywords
fault
traveling wave
line
zero
phase
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
CN201510650313.6A
Other languages
English (en)
Other versions
CN105353268B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510650313.6A priority Critical patent/CN105353268B/zh
Publication of CN105353268A publication Critical patent/CN105353268A/zh
Application granted granted Critical
Publication of CN105353268B publication Critical patent/CN105353268B/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
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • 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
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/085Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead

Abstract

本发明公开了一种用于输电线路分布式行波故障判断及定位方法,涉及输电线路故障定位技术领域。该方法故障定位步骤如下:S1、设在变电站M、N的母线上以及架空输电线上装有n组检测装置;S2、检测并读取每组检测装置上的电流信息;S3、对三相故障电流行波进行相模变换;S4、对故障点故障发生时刻的电压进行过零识别;S5、根据S4的结果选取小波变换或二分SVD算法;S6、根据故障行波波头信息,确定故障区间、计算行波波速和故障距离。本发明解决了系统单独使用小波变换算法时,不能在电压过零点附近发生故障的情况下对故障距离进行准确定位以及系统单独使用二分SVD算法时数据处理效率低,运算量大等问题。

Description

一种用于输电线路分布式行波故障判断及定位方法
技术领域
[0001] 本发明涉及输电线路故障定位技术领域。
背景技术
[0002] 电力是国民经济的命脉,输电线路肩负着输送电能的重任,是电力系统的重要组 成部分,提高电力系统的安全稳定运行水平、降低运行维护成本,提高经济效益,向用户提 供高质量的电能是电力工作的基本任务。由于输电线路分布区域广,通常架设在无人看管 的野外,途径山区、丘陵、江河等多种地形,加之气候环境复杂,会遇上雷雨、覆冰、强风等气 象条件,是电力系统最容易发生故障的电力设备之一。快速、可靠、准确地进行故障定位,及 时发现和处理绝缘隐患,可大大提高电网的可靠性和愈合能力。
[0003] 目前输电线路故障定位主要由行波法和故障分析法两种,行波法由于定位原理简 单在电力系统中得到大量的应用。但是由于输电线路故障行波的非线性衰减、线路噪声等 问题,使得传统的故障电流行波测距系统在应用中还存在一定的问题。近几年先后有发明 人提出行波故障定位方法,如专利201410222008. 2、201410041394. 5和201210496668. 0分 别提出了不同的行波故障定位方法,总体上看都是通过利用小波变换和反射波来对故障行 波进行定位。但是由于故障行波的衰减特性、过零特性以及实际装置采样率的限制,如果单 独使用小波变换算法会导致某些时刻测距失效。以上专利的测距方法存在以下几点不足:
[0004] 1.在某些故障区域会利用反射波进行测距,由于实际装置采样率的限制以及线路 自身特点,当故障点距离母线非常近时会产生死区;
[0005] 2.在某些故障区域会利用反射波进行测距,由于行波在线路上传播具有衰减特 性,当故障点距离检测装置比较远时,检测装置可能检测不到反射波。
[0006] 3.当电压过零点附近发生故障时,单独使用小波变换会导致测距失败。
发明内容
[0007] 本发明的目的是针对现有测距技术中的不足提出了一种基于二分SVD理论与小 波变换结合的改进分布式行波故障判断及定位方法,以解决系统中在任意时刻发生故障下 都能可靠的对线路故障进行定位。本发明系统由沿线安装在架空输电线路上的若干组(每 组3个,对应A、B、C三相)故障检测装置组成,如图1所示。
[0008] 本发明的技术方案是一种用于输电线路分布式行波故障判断及定位方法:
[0009] 步骤1 :检测并读取每个检测点的三相电流行波数据;
[0010] 步骤2 :三相线路相模变换
[0011] 采用相模变换矩阵T对三相电流行波进行解耦,消除相间的耦合影响,获得解耦 后的3个线模行波:
Figure CN105353268AD00051
[0014] 矩阵I^#iA、iB、ic三相电流行波变换为ipi。、%、、模分量,其中i。、%、\为 线模分量,ic为零模分量,若发生接地故障,则存在零模分量;
[0015] 步骤3 :电压过零点附近故障的识别;
[0016] 步骤3. 1 :计算3个线模行波的能量,选取3个线模行波中能量最大的一个;
[0017] 步骤3. 2 :对选取的线模行波进行小波变换,并取出21尺度的小波变换结果W(n);
[0018] 步骤3. 3 :找到W (η)的最大值
[0019] 步骤3. 4 :在W(n)中寻找小于KwW_的点数N,其中Kw为调整系数需根据实际情况 进行设定;
[0020] 步骤3. 5 :若^^《表明发生了电压过零点附近故障,其中n为21尺度小波变换结 η 果的总长度,Κ为阈值,根据实际情况设定;
[0021] 步骤4 :根据步骤3的结果,若发生电压过零点附近故障则,选取二分SVD算法进 行故障测距,否则使用小波变换进行故障测距;
[0022] 步骤5 :故障位置初步判定;
[0023] 若故障发生在i和i+Ι检测点之间(1 < i〈n),根据行波极性特征有
[0024] P1= p2 =…=ρι+1= Ρι+2 =…=pn (3)
[0025] 步骤6:故障点精确定位
[0026] 设故障发生在i和i+Ι检测点之间(1 < i〈n),检测点i和i+Ι第一次检测到故障 行波波头的时间分别为h,t1+1,检测点i和i+Ι离母线Μ端的距离分别为LdP L 1+1,则行波 波速和故障点位置可以通过以下公式计算得出:
Figure CN105353268AD00052
[0029] 式中v表示行波波速,L表示任意相邻两个检测点间的距离,Λ t表示与行波传播 方向一致上的任意相邻两个检测点第一次检测到故障行波的时间差,X表示母线Μ到故障 点的距离,所述小波变换采用三次Β样条小波基。
[0030] 所述二分SVD变换如下:
[0031] 设离散数字信号Χ= [χ(1),χ(2),···,χ(Ν)],可通过X构造Hankel矩阵如下式:
Figure CN105353268AD00061
[0033] 通过SVD分解可得Μ = qwpf +σ2ω2ν『,σ ^ σ 2为SVD分解后两个奇异值;由于 A中的第二行仅比第一行滞后一个数据,因此有〇1>>〇2;由σ i得到的分量信号六1反映了 原始信号的近似成分,σ 2得到的分量信号D 映了原始信号的细节成分;若将近似信号 4再进行SVD分解,不断重复上述过程,就可对原始信号分解成一系列的近似信号和细节信 号;细节信号反映了行波的奇异特性。
[0034] 本发明解决了系统单独使用小波变换算法时,不能在电压过零点附近发生故障的 情况下对故障距离进行准确定位以及系统单独使用二分SVD算法时数据处理效率低,运算 量大等问题。同时改进分布式行波测距方法解决了实际工程中由于反射波的利用而导致存 在死区或捕捉不到二次反射波的问题。本方法具有实时性强、故障定位可靠等特点。
附图说明
[0035] 图1是改进双端故障测距安装位置示意图;
[0036] 图2是故障测距仿真模型图;
[0037] 图3是二分SVD分解过程示意图:
[0038] 图4是四点分布式定位线路图:
[0039] 图5是一般故障时小波变换(左)和二分SVD(右)检测结果图;
[0040] 图6是电压过零点附近故障是小波变换(左)和二分SVD(右)检测结果图。
具体实施方式
[0041] 下面结合附图,给出本发明的较佳实施例,并予以详细描述。
[0042] 如图4所示,以A相接地故障为例,输电线路总长为120km,上面装有4组检测装 置,每组检测装置相隔30km。故障分别发生在合闸角为0° (电压过零)和10° (电压未 过零)时刻,采样频率为1MHz。
[0043] 步骤一:检测并读取4组检测装置的故障电流信息。
[0044] 步骤二:利用T矩阵对故障数据进行相模变换。
[0045] 步骤三:选取21尺度的小波对故障数据进行过零检测,根据线路参数取K w= 0. 07, K = 0. 05。不同故障点的识别结果如表1所示。
[0046] 步骤四:根据电压过零点附近故障识别结果选取三次B样条小波变换或二分SVD 算法对行波的波头进行识别。图5为电压过零点附近故障是小波变换和二分SVD检测结果, 图6为一般故障时小波变换和二分SVD检测结果。两图结果表明,在任意时刻任意位置故 障本系统都能很好的检测到行波波头(行波模极大值处)。
[0047] 步骤五:利用检测到的行波波头,根据S6中的距离计算公式,即可求得故障点的 位置。计算结果与误差如表2所示。
[0048] 结果表明,任意时刻任意点处发生故障,本方法都能很好的进行故障定位,同时由 于计算中未利用反射波,消除了死区的影响。
[0049] 表1电压过零点附近故障识别结果
Figure CN105353268AD00071
[0051] 表2合闸角不同时单相接地短路行波故障测距结果
Figure CN105353268AD00072

Claims (2)

1. 一种用于输电线路分布式行波故障判断及定位方法,该方法包括 步骤1 :检测并读取每个检测点的三相电流行波数据; 步骤2 :三相线路相模变换 采用相模变换矩阵T对三相电流行波进行解耦,消除相间的耦合影响,获得解耦后的3 个线模行波:
Figure CN105353268AC00021
矩阵1'将^、^、、三相电流行波变换为1。、1、%、\模分量,其中1。、%、\为线模 分量,i。为零模分量,若发生接地故障,则存在零模分量; 步骤3 :电压过零点附近故障的识别; 步骤3. 1 :计算3个线模行波的能量,选取3个线模行波中能量最大的一个; 步骤3. 2 :对选取的线模行波进行小波变换,并取出21尺度的小波变换结果W(n); 步骤3. 3 :找到W(n)的最大值 步骤3. 4 :在W(η)中寻找小于KwWniax的点数N,其中K w为调整系数需根据实际情况进行 设定; 步骤3. 5 :若
Figure CN105353268AC00022
表明发生了电压过零点附近故障,其中η为21尺度小波变换结果 的总长度,K为阈值,根据实际情况设定; 步骤4 :根据步骤3的结果,若发生电压过零点附近故障则,选取二分SVD算法进行故 障测距,否则使用小波变换进行故障测距; 步骤5 :故障位置初步判定; 若故障发生在i和i+Ι检测点之间(I < i〈n),根据行波极性特征有 P1= P 2=…=P 夕 P 1+1= P 1+2=…=P n (3) 步骤6 :故障点精确定位 设故障发生在i和i+Ι检测点之间(I < i〈n),检测点i和i+Ι第一次检测到故障行波 波头的时间分别为h,t1+1,检测点i和i+Ι离母线M端的距离分别为LjP L 1+1,则行波波速 和故障点位置可以通过以下公式计算得出:
Figure CN105353268AC00023
式中V表示行波波速,L表示任意相邻两个检测点间的距离,Λ t表示与行波传播方向 一致上的任意相邻两个检测点第一次检测到故障行波的时间差,X表示母线M到故障点的 距离,所述小波变换采用三次B样条小波基。
2.如权利要求1所述的一种用于输电线路分布式行波故障判断及定位方法,其特征在 于所述步骤4中二分SVD变换如下: 设离散数字信号X = [X (I),X (2),…,X (N)],可通过X构造 Hankel矩阵如下式:
Figure CN105353268AC00031
通过SVD分解可得
Figure CN105353268AC00032
σ JP σ 2为SVD分解后两个奇异值;由于A中 的第二行仅比第一行滞后一个数据,因此有O1Ao2;由〇 i得到的分量信号A1反映了原始 信号的近似成分,σ 2得到的分量信号D 映了原始信号的细节成分;若将近似信号A 1再 进行SVD分解,不断重复上述过程,就可对原始信号分解成一系列的近似信号和细节信号; 细节信号反映了行波的奇异特性。
CN201510650313.6A 2015-10-10 2015-10-10 一种用于输电线路分布式行波故障判断及定位方法 Expired - Fee Related CN105353268B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510650313.6A CN105353268B (zh) 2015-10-10 2015-10-10 一种用于输电线路分布式行波故障判断及定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510650313.6A CN105353268B (zh) 2015-10-10 2015-10-10 一种用于输电线路分布式行波故障判断及定位方法

Publications (2)

Publication Number Publication Date
CN105353268A true CN105353268A (zh) 2016-02-24
CN105353268B CN105353268B (zh) 2018-04-06

Family

ID=55329270

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510650313.6A Expired - Fee Related CN105353268B (zh) 2015-10-10 2015-10-10 一种用于输电线路分布式行波故障判断及定位方法

Country Status (1)

Country Link
CN (1) CN105353268B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093702A (zh) * 2016-06-06 2016-11-09 南京理工大学 一种考虑多点故障的高压输电线路行波测距方法
CN106455041A (zh) * 2016-12-14 2017-02-22 国家电网公司 基于输电线路的电力无线传感器网络时间同步方法及其系统
CN106771861A (zh) * 2016-12-08 2017-05-31 中国矿业大学 基于广域行波能量和时间差的复杂电网故障定位方法
CN106959401A (zh) * 2017-03-24 2017-07-18 国网山东省电力公司 一种分布式行波测距方法
CN108196166A (zh) * 2018-01-11 2018-06-22 武汉三相电力科技有限公司 基于故障波形起始点到达时差的双端行波测距方法
CN108646145A (zh) * 2018-07-26 2018-10-12 南方电网科学研究院有限责任公司 一种输电线路闪络塔定位方法及系统
CN109142966A (zh) * 2018-07-13 2019-01-04 南方电网科学研究院有限责任公司 基于线路实测数据的故障测距方法、装置、设备及介质
CN109387744A (zh) * 2018-12-17 2019-02-26 国网山东省电力公司电力科学研究院 基于奇异值分解的配网线路故障点定位方法及装置
CN109901021A (zh) * 2019-04-08 2019-06-18 珠海妙微科技有限公司 基于二维Hankel矩阵多尺度SVD变换法
CN110161376A (zh) * 2019-06-24 2019-08-23 四川电安智能科技有限公司 一种行波故障时刻提取算法
CN110426592A (zh) * 2019-08-16 2019-11-08 南京国电南自电网自动化有限公司 基于行波时间差的架空和电缆混合线路接地故障定位方法
CN110703012A (zh) * 2019-09-26 2020-01-17 国电南瑞科技股份有限公司 一种输电线路分布式故障诊断方法
CN110895298A (zh) * 2019-12-17 2020-03-20 云南拓洲科技有限公司 一种10kv和35kv架空线路暂态故障监测指示系统
CN111766470A (zh) * 2020-06-24 2020-10-13 湖南大学 高压直流输电线路的故障定位方法、系统及直流输电线路

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101187687A (zh) * 2007-12-21 2008-05-28 清华大学 基于暂态行波的高阻接地故障检测方法
US20090085574A1 (en) * 2007-09-28 2009-04-02 Yan Pan Method for locating phase to ground faults in DC distribution systems
CN101871990A (zh) * 2010-05-24 2010-10-27 北京交通大学 架空输电线路故障定位系统
CN101923139A (zh) * 2010-04-19 2010-12-22 昆明理工大学 一种输电线路单端行波故障测距智能方法
CN102520315A (zh) * 2011-12-05 2012-06-27 西南交通大学 基于行波多尺度信息的输电线路故障单端定位方法
US20130245969A1 (en) * 2012-03-15 2013-09-19 Applied Materials, Inc. Method to detect wafer arcing in semiconductor manufacturing equipment

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090085574A1 (en) * 2007-09-28 2009-04-02 Yan Pan Method for locating phase to ground faults in DC distribution systems
CN101187687A (zh) * 2007-12-21 2008-05-28 清华大学 基于暂态行波的高阻接地故障检测方法
CN101923139A (zh) * 2010-04-19 2010-12-22 昆明理工大学 一种输电线路单端行波故障测距智能方法
CN101871990A (zh) * 2010-05-24 2010-10-27 北京交通大学 架空输电线路故障定位系统
CN102520315A (zh) * 2011-12-05 2012-06-27 西南交通大学 基于行波多尺度信息的输电线路故障单端定位方法
US20130245969A1 (en) * 2012-03-15 2013-09-19 Applied Materials, Inc. Method to detect wafer arcing in semiconductor manufacturing equipment

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘亚东: "输电线路分布式故障测距理论与关键技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093702B (zh) * 2016-06-06 2019-04-12 南京理工大学 一种考虑多点故障的高压输电线路行波测距方法
CN106093702A (zh) * 2016-06-06 2016-11-09 南京理工大学 一种考虑多点故障的高压输电线路行波测距方法
CN106771861B (zh) * 2016-12-08 2019-01-25 中国矿业大学 基于广域行波能量和时间差的复杂电网故障定位方法
CN106771861A (zh) * 2016-12-08 2017-05-31 中国矿业大学 基于广域行波能量和时间差的复杂电网故障定位方法
CN106455041A (zh) * 2016-12-14 2017-02-22 国家电网公司 基于输电线路的电力无线传感器网络时间同步方法及其系统
CN106959401B (zh) * 2017-03-24 2019-11-15 国网山东省电力公司 一种分布式行波测距方法
CN106959401A (zh) * 2017-03-24 2017-07-18 国网山东省电力公司 一种分布式行波测距方法
CN108196166A (zh) * 2018-01-11 2018-06-22 武汉三相电力科技有限公司 基于故障波形起始点到达时差的双端行波测距方法
CN108196166B (zh) * 2018-01-11 2020-02-21 武汉三相电力科技有限公司 基于故障波形起始点到达时差的双端行波测距方法
CN109142966A (zh) * 2018-07-13 2019-01-04 南方电网科学研究院有限责任公司 基于线路实测数据的故障测距方法、装置、设备及介质
CN108646145B (zh) * 2018-07-26 2021-03-19 南方电网科学研究院有限责任公司 一种输电线路闪络塔定位方法及系统
CN108646145A (zh) * 2018-07-26 2018-10-12 南方电网科学研究院有限责任公司 一种输电线路闪络塔定位方法及系统
CN109387744A (zh) * 2018-12-17 2019-02-26 国网山东省电力公司电力科学研究院 基于奇异值分解的配网线路故障点定位方法及装置
CN109901021B (zh) * 2019-04-08 2021-07-16 珠海妙微科技有限公司 基于二维Hankel矩阵多尺度SVD变换法
CN109901021A (zh) * 2019-04-08 2019-06-18 珠海妙微科技有限公司 基于二维Hankel矩阵多尺度SVD变换法
CN110161376B (zh) * 2019-06-24 2021-05-28 四川电安智能科技有限公司 一种行波故障时刻提取算法
CN110161376A (zh) * 2019-06-24 2019-08-23 四川电安智能科技有限公司 一种行波故障时刻提取算法
CN110426592A (zh) * 2019-08-16 2019-11-08 南京国电南自电网自动化有限公司 基于行波时间差的架空和电缆混合线路接地故障定位方法
CN110703012A (zh) * 2019-09-26 2020-01-17 国电南瑞科技股份有限公司 一种输电线路分布式故障诊断方法
CN110703012B (zh) * 2019-09-26 2021-10-15 国电南瑞科技股份有限公司 一种输电线路分布式故障诊断方法
CN110895298A (zh) * 2019-12-17 2020-03-20 云南拓洲科技有限公司 一种10kv和35kv架空线路暂态故障监测指示系统
CN111766470A (zh) * 2020-06-24 2020-10-13 湖南大学 高压直流输电线路的故障定位方法、系统及直流输电线路

Also Published As

Publication number Publication date
CN105353268B (zh) 2018-04-06

Similar Documents

Publication Publication Date Title
CN105353268A (zh) 一种用于输电线路分布式行波故障判断及定位方法
CN102520315B (zh) 基于行波多尺度信息的输电线路故障单端定位方法
Suonan et al. A novel fault-location method for HVDC transmission lines
CN101299538B (zh) 一种电缆-架空混合线路故障行波测距方法
CN202794433U (zh) 基于云计算平台的配电网故障行波定位装置
CN102074937B (zh) 模型自适应输电线路故障测距方法
CN102914726A (zh) 一种用于同杆双回线路的故障定位方法
CN106019079B (zh) 一种同塔双回直流线路新型双端故障测距方法
CN108375713A (zh) 一种新型电网故障行波定位方法和系统
CN108693444B (zh) 一种基于电磁信号时间反演的输电线路故障点定位方法
CN103364693A (zh) 一种基于区域数据的输电线路行波故障测距方法
CN107632236A (zh) 一种基于对端母线反射波识别的单出线输电线路单端故障测距方法
CN106771861B (zh) 基于广域行波能量和时间差的复杂电网故障定位方法
CN110907755A (zh) 一种输电线路在线监测故障识别方法
CN104865495A (zh) 电力线路故障行波定位方法及系统
CN104820168A (zh) 一种基于波形差异度和雷击故障样本库的雷击故障判别方法
CN109387744A (zh) 基于奇异值分解的配网线路故障点定位方法及装置
CN103823156A (zh) 一种带容错功能的输电线路分布式故障定位方法
CN201993433U (zh) 高压直流输电接地极线路行波测距系统
CN103777115A (zh) 基于故障暂态与稳态信号波速差的输电线路单端定位方法
He et al. Natural frequency‐based protection scheme for voltage source converter‐based high‐voltage direct current transmission lines
CN103197204A (zh) 多端线路故障定位的混合型方法
CN105548820A (zh) 基于暂态行波模极大值法的电缆线路故障识别与定位方法
RU2532760C1 (ru) Способ определения места повреждения разветвленной линии электропередачи
CN104198890A (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
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: 20180406

Termination date: 20201010