CN102998535A - 一种基于极大似然估计理论的***谐波阻抗计算方法 - Google Patents
一种基于极大似然估计理论的***谐波阻抗计算方法 Download PDFInfo
- Publication number
- CN102998535A CN102998535A CN2012105285592A CN201210528559A CN102998535A CN 102998535 A CN102998535 A CN 102998535A CN 2012105285592 A CN2012105285592 A CN 2012105285592A CN 201210528559 A CN201210528559 A CN 201210528559A CN 102998535 A CN102998535 A CN 102998535A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- sigma
- harmonic
- maximum likelihood
- likelihood estimation
- 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
Images
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了电力技术领域的一种基于极大似然估计理论的***谐波阻抗估计方法。其技术方案是,采集公共连接点的母线电压瞬时值和用户接入***的电流瞬时值,建立谐波电压相量和谐波电流相量的关系;在定义复协方差的基础上,推导得到一元复正态分布的概率密度函数,进而得到极大似然估计函数,建立***谐波阻抗估计的复数域极大似然估计理论;利用极值理论求解极大似然估计函数,最终得到***谐波阻抗估计值。本发明的有益效果是,基于极大似然估计理论的***等值谐波阻抗计算方法可以比较精确地计算***等值谐波阻抗,这对进一步解决谐波污染问题、提高电能质量管理水平具有重要意义。
Description
技术领域
本发明属于谐波阻抗计算方法设计领域,尤其涉及一种基于极大似然估计理论的***谐波阻抗计算方法。
背景技术
在传统能源和特高压交直流输电网络的基础上,大量非线性负荷设备接入电网,产生了严重的电力***谐波问题。谐波管理和谐波治理是电力***谐波问题的两个主要内容。在谐波管理方面,目前基于国标的谐波管理制度造成了谐波源用户缺乏谐波治理的主动性,定量确定谐波源污染责任是解决这一问题的关键。在谐波治理方面,安装滤波器是谐波治理的一个主要方法,滤波器在安装之前必须进行滤波器设计。无论是定量确定谐波源的污染责任还是滤波器设计,都需要准确估计***谐波阻抗。
投切电容器是传统的谐波阻抗测量方法,该方法可以方便地改变电网的运行方式,创造出阻抗测量的条件,但是会对电网产生干扰。直接利用电压和电流的监测值计算***谐波阻抗的方法可以避免对电网产生干扰,目前该类方法的实现途径是利用测量数据在实数域内对谐波阻抗进行回归。采用实数域进行回归分析的缺点是谐波阻抗估计的误差不可控,其根本原因是回归分析的过程中将电压和电流相量的实部和虚部分别进行计算,导致了两个方面的结果:一方面,回归方程中被回归的参数是变量而不是常量;另一方面,回归方程中自变量和因变量之间的线性关系得不到保证。由于相量在数学上对应于复数,因此,提出复数域上的基于极大似然估计理论的***谐波阻抗计算方法不仅可以避免对非常量进行回归,还可避免回归方程中自变量和因变量之间严格线性关系的破坏,从而提高计算准确性。
发明内容
针对上述背景技术中提到的实数域回归方法在***谐波阻抗计算方面的不足,本发明提出了一种基于极大似然估计理论的***谐波阻抗计算方法。
一种基于极大似然估计理论的***谐波阻抗计算方法,其特征在于,具体包括以下步骤:
步骤1:采集公共连接点的母线电压瞬时值和用户接入***的电流瞬时值,通过傅里叶变换得到谐波电压和谐波电流相量数据序列;并根据电路理论建立谐波电压相量和谐波电流相量关系;
步骤2:定义复协方差的基础上,推导得到一元复正态分布的概率密度函数;
步骤3:基于一元复正态分布的概率密度函数,推导得到极大似然估计函数,从而建立***谐波阻抗估计的复数域极大似然估计理论;
步骤4:利用极值理论求解极大似然估计函数,最终得到***谐波阻抗估计值。
步骤1中,谐波电压和谐波电流相量数据序列
根据电路理论,谐波电压相量和谐波电流相量关系为:
其中,表示公共连接点处接入的某个用户D接入***的第h次谐波电流相量值,表示公共连接点处第h次谐波电压相量值,为背景谐波电压相量,ZhX为***谐波阻抗,ε是测量误差项,ε1,ε2,…,εn独立同分布且服从均值和方差分别为0与σ2的复正态随机分布;n为电压相量和谐波电流相量数据序列的个数。
步骤2中,推导得到一元复正态分布的概率密度函数的过程为:
设复随机变量Z=X+iY的实部X和虚部Y都服从正态分布,则称复随机变量Z为复正态随机变量。由p个复正态随机变量构成的向量ξ′=(Z1,Z2,…,Zp)称为多元复正态随机变量。据此,由p元复正态随机变量的实部和虚部构成的随机变量η′=(X1,Y1,…,Xp,Yp)是服从多元正态分布的实随机变量。
记η′=(x1,y1,…,xp,yp),ξ′=(z1,z2,…,zp),则有下式:
其中,Eη表示多元正态分布随机变量η的期望,∑η表示多元随机变量η的协方差,f(·)表示随机变量的概率密度函数。
其中,
多元正态分布随机变量η的协方差阵∑η为:
设定:
则有:
结合 和 得到多元复正态分布的概率密度函数:
则,当p=1时:一元复正态分布的概率密度函数为: 其中
其中,z为一元复随机变量,σ2是复随机变量z的方差,exp为指数运算,E表示期望。
步骤3中,基于一元复正态分布的概率密度函数,推导得到极大似然估计函数,从而建立***谐波阻抗估计的复数域极大似然估计理论的过程为:
根据谐波电压相量和谐波电流相量关系:
步骤3中,利用最优化理论的极值定理求解极大似然估计函数的过程为:
将极大似然估计函数变换为对数形式:
得到对数似然函数取最大时谐波阻抗ZhX的估计结果:
本发明的有益效果是,基于极大似然估计理论的***等值谐波阻抗计算方法可以比较精确地计算***等值谐波阻抗,这对进一步解决谐波污染问题、提高电能质量管理水平具有重要意义。
附图说明
图1是本发明提供的一种基于极大似然估计理论的***谐波阻抗计算方法的公共连接点处网络示意图;
图2是本发明提供的一种基于极大似然估计理论的***谐波阻抗计算方法的IEEE14节点标准测试***;
图3是本发明提供的一种基于极大似然估计理论的***谐波阻抗计算方法的IEEE14节点标准测试***仿真开始0.02s谐波源HL1接入的5次谐波电流曲线;其中,(a)是5次谐波电流的实部;(b)是5次谐波电流的虚部;
具体实施方式
下面结合附图,对优选的实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
图2是本发明提供的一种基于极大似然估计理论的***谐波阻抗计算方法的IEEE14节点标准测试***。该测试***由2台发电机、3台同步调相机、14条母线、15条输电线路和3台变压器组成。
选择母线11为关注母线,HL1、HL2及L3为母线处接入的三个负荷,将谐波源HL1选择为如图1所示的用户D,HL2为负荷侧其余非线性负荷,L3为线性负荷,同时母线13处接有谐波源负荷HS,以5次谐波为例,计算除用户D以外其余部分***等效的谐波阻抗。
设定谐波源HL2的5次谐波电流幅值为10.00A,初始相角为-74.25°;***侧谐波源HS的5次谐波电流幅值为114.80A,初始相角为-76.56°。设定谐波源HL1的5次谐波电流基准幅值为20A,并且幅值在0.01~1.1倍基准幅值间随机波动,相角在-180°~180°间随机波动。设定采样频率6.4kHz,仿真时间为1.4s。然后执行以下步骤:
步骤1:采集母线11处的母线电压瞬时值和谐波源HL1接入***的电流瞬时值,通过傅里叶变换得到谐波电压和谐波电流相量数据序列:
表示母线11处接入的谐波源HL1(即图1示意的用户D)第5次谐波电流相量值,表示母线11处第5次谐波电压相量值,为背景谐波电压相量,Z5X为除谐波源HL1的5次***谐波阻抗。εk是测量误差项,ε1,ε2,…,ε1440独立同分布且服从均值和方差分别为0与σ2的复正态随机分布;
步骤2:利用步骤1得到的谐波电压和谐波电流相量数据序列,建立极大似然函数:
步骤3:求解极大似然函数,得到***谐波阻抗的计算公式:
利用上述公式,计算得到***等值谐波阻抗为3.0240∠72.7616°。基于极大似然估计理论的***谐波阻抗估计方法计算结果与真实值的对比如下表所示,符号“-”表示此值不需计算。
表1***等值谐波阻抗计算结果对比
从表1可以看出,基于极大似然估计理论的***等值谐波阻抗计算方法可以比较精确地计算***等值谐波阻抗,提高了计算精度。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应该涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (6)
1.一种基于极大似然估计理论的***谐波阻抗计算方法,其特征在于,具体包括以下步骤:
步骤1:采集公共连接点的母线电压瞬时值和用户接入***的电流瞬时值,通过傅里叶变换得到谐波电压和谐波电流相量数据序列;并根据电路理论建立谐波电压相量和谐波电流相量的关系;
步骤2:定义复协方差的基础上,推导得到一元复正态分布的概率密度函数;
步骤3:基于一元复正态分布的概率密度函数,推导得到极大似然估计函数,从而建立***谐波阻抗估计的复数域极大似然估计理论;
步骤4:利用极值理论求解极大似然估计函数,最终得到***谐波阻抗估计值。
4.根据权利要求1所述的一种基于极大似然估计理论的***谐波阻抗计算方法,其特征在于,所述一元复正态分布的概率密度函数为:
其中,z为一元复随机变量,σ2是复随机变量z的方差,exp为指数运算,E表示期望。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210528559.2A CN102998535B (zh) | 2012-12-10 | 2012-12-10 | 一种基于极大似然估计理论的***谐波阻抗计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210528559.2A CN102998535B (zh) | 2012-12-10 | 2012-12-10 | 一种基于极大似然估计理论的***谐波阻抗计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102998535A true CN102998535A (zh) | 2013-03-27 |
CN102998535B CN102998535B (zh) | 2014-09-03 |
Family
ID=47927360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210528559.2A Expired - Fee Related CN102998535B (zh) | 2012-12-10 | 2012-12-10 | 一种基于极大似然估计理论的***谐波阻抗计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102998535B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218534A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN103809023A (zh) * | 2014-01-26 | 2014-05-21 | 西南交通大学 | 基于二分搜索的电网同步谐波相量测量方法 |
CN103995958A (zh) * | 2014-04-23 | 2014-08-20 | 广东易事特电源股份有限公司 | 一种谐波阻抗特征函数构建方法 |
CN105588982A (zh) * | 2014-10-27 | 2016-05-18 | 国家电网公司 | 一种电力***谐波阻抗的测量计算方法 |
CN105606900A (zh) * | 2016-03-18 | 2016-05-25 | 华南理工大学 | 一种基于方波信号的单相谐波阻抗测量方法 |
CN106154040A (zh) * | 2015-03-24 | 2016-11-23 | 国家电网公司 | 一种用于计算电网单点等值时等效谐波阻抗方法 |
CN106300331A (zh) * | 2015-06-04 | 2017-01-04 | 中国电力科学研究院 | 一种精确计算配电网支路阻抗的方法 |
CN106651150A (zh) * | 2016-12-01 | 2017-05-10 | 全球能源互联网研究院 | 一种评估电网谐波影响的方法及装置 |
CN107576853A (zh) * | 2017-07-12 | 2018-01-12 | 国网上海市电力公司 | 基于典型相关性分析的配电网谐波阻抗计算方法 |
CN107576851A (zh) * | 2017-07-26 | 2018-01-12 | 华南理工大学 | 基于旋转电流相量的***谐波阻抗测量方法 |
CN108802491A (zh) * | 2018-07-06 | 2018-11-13 | 四川大学 | 基于谐波源独立性的谐波发射水平估计求解新方法 |
CN110007147A (zh) * | 2018-11-29 | 2019-07-12 | 广西电网有限责任公司电力科学研究院 | 一种单母线多用户的谐波责任计算方法 |
CN110456159A (zh) * | 2019-08-21 | 2019-11-15 | 四川大学 | 基于修正独立随机矢量的***侧谐波阻抗估计方法及*** |
CN110474325A (zh) * | 2019-08-06 | 2019-11-19 | 南方电网科学研究院有限责任公司 | 基于扇形边界的交流***谐波阻抗的边界构建方法及装置 |
CN110518577A (zh) * | 2019-08-06 | 2019-11-29 | 南方电网科学研究院有限责任公司 | 交流***谐波阻抗的多边形边界构建方法及装置 |
CN111077372A (zh) * | 2019-12-24 | 2020-04-28 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种交直流输电***谐波阻抗测定方法及*** |
CN115032458A (zh) * | 2022-06-02 | 2022-09-09 | 北京妙微科技有限公司 | 一种基于电网***的谐波阻抗估计方法及计算机存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09189729A (ja) * | 1996-01-09 | 1997-07-22 | Nissin Electric Co Ltd | 配電系統の高調波インピーダンス測定方法 |
JP2001004685A (ja) * | 1999-06-22 | 2001-01-12 | Chubu Electric Power Co Inc | 電力系統の高調波特性測定方法 |
CN101706532A (zh) * | 2009-11-25 | 2010-05-12 | 国网电力科学研究院武汉南瑞有限责任公司 | 一种谐波阻抗测量方法及测量装置 |
CN102749521A (zh) * | 2012-07-18 | 2012-10-24 | 华北电力大学(保定) | 一种电力***谐波阻抗计算方法 |
-
2012
- 2012-12-10 CN CN201210528559.2A patent/CN102998535B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09189729A (ja) * | 1996-01-09 | 1997-07-22 | Nissin Electric Co Ltd | 配電系統の高調波インピーダンス測定方法 |
JP2001004685A (ja) * | 1999-06-22 | 2001-01-12 | Chubu Electric Power Co Inc | 電力系統の高調波特性測定方法 |
CN101706532A (zh) * | 2009-11-25 | 2010-05-12 | 国网电力科学研究院武汉南瑞有限责任公司 | 一种谐波阻抗测量方法及测量装置 |
CN102749521A (zh) * | 2012-07-18 | 2012-10-24 | 华北电力大学(保定) | 一种电力***谐波阻抗计算方法 |
Non-Patent Citations (4)
Title |
---|
SHUSEN LI ET AL.: "Dual Synchronization Incremental Method for Harmonic Impedance Measurement", 《HARMONICS AND QUALITY OF POWER(ICHQP),2012 IEEE 15TH INTERNATIONAL CONFERENCE ON》 * |
T.GOURAUD ET AL.: "A Maximum Likelihood Approach To Harmonics Measurements In Power Systems", 《COMPUTERS IN POWER ELECTRONICS,1994,IEEE 4TH WORKSHOP ON》 * |
吕洋等: "投切电容器时电网谐波阻抗测量方法", 《高压电技术》 * |
郭福全: "探究电力***谐波阻抗测量", 《民营科技》 * |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218534A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN103218534B (zh) * | 2013-04-22 | 2017-02-08 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN103809023A (zh) * | 2014-01-26 | 2014-05-21 | 西南交通大学 | 基于二分搜索的电网同步谐波相量测量方法 |
CN103809023B (zh) * | 2014-01-26 | 2016-08-24 | 西南交通大学 | 基于二分搜索的电网同步谐波相量测量方法 |
CN103995958A (zh) * | 2014-04-23 | 2014-08-20 | 广东易事特电源股份有限公司 | 一种谐波阻抗特征函数构建方法 |
CN103995958B (zh) * | 2014-04-23 | 2017-02-15 | 广东易事特电源股份有限公司 | 一种谐波阻抗特征函数构建方法 |
CN105588982A (zh) * | 2014-10-27 | 2016-05-18 | 国家电网公司 | 一种电力***谐波阻抗的测量计算方法 |
CN106154040B (zh) * | 2015-03-24 | 2019-09-06 | 国家电网公司 | 一种用于计算电网单点等值时等效谐波阻抗方法 |
CN106154040A (zh) * | 2015-03-24 | 2016-11-23 | 国家电网公司 | 一种用于计算电网单点等值时等效谐波阻抗方法 |
CN106300331A (zh) * | 2015-06-04 | 2017-01-04 | 中国电力科学研究院 | 一种精确计算配电网支路阻抗的方法 |
CN106300331B (zh) * | 2015-06-04 | 2019-02-22 | 中国电力科学研究院 | 一种精确计算配电网支路阻抗的方法 |
CN105606900B (zh) * | 2016-03-18 | 2019-01-18 | 华南理工大学 | 一种基于方波信号的单相谐波阻抗测量方法 |
CN105606900A (zh) * | 2016-03-18 | 2016-05-25 | 华南理工大学 | 一种基于方波信号的单相谐波阻抗测量方法 |
CN106651150B (zh) * | 2016-12-01 | 2021-02-09 | 全球能源互联网研究院有限公司 | 一种评估电网谐波影响的方法及装置 |
CN106651150A (zh) * | 2016-12-01 | 2017-05-10 | 全球能源互联网研究院 | 一种评估电网谐波影响的方法及装置 |
CN107576853A (zh) * | 2017-07-12 | 2018-01-12 | 国网上海市电力公司 | 基于典型相关性分析的配电网谐波阻抗计算方法 |
CN107576853B (zh) * | 2017-07-12 | 2019-11-05 | 国网上海市电力公司 | 基于典型相关性分析的配电网谐波阻抗计算方法 |
CN107576851B (zh) * | 2017-07-26 | 2019-05-14 | 华南理工大学 | 基于旋转电流相量的***谐波阻抗测量方法 |
CN107576851A (zh) * | 2017-07-26 | 2018-01-12 | 华南理工大学 | 基于旋转电流相量的***谐波阻抗测量方法 |
CN108802491B (zh) * | 2018-07-06 | 2019-07-19 | 四川大学 | 基于谐波源独立性的谐波发射水平估计求解新方法 |
CN108802491A (zh) * | 2018-07-06 | 2018-11-13 | 四川大学 | 基于谐波源独立性的谐波发射水平估计求解新方法 |
CN110007147A (zh) * | 2018-11-29 | 2019-07-12 | 广西电网有限责任公司电力科学研究院 | 一种单母线多用户的谐波责任计算方法 |
CN110518577B (zh) * | 2019-08-06 | 2020-12-18 | 南方电网科学研究院有限责任公司 | 交流***谐波阻抗的多边形边界构建方法及装置 |
CN110474325A (zh) * | 2019-08-06 | 2019-11-19 | 南方电网科学研究院有限责任公司 | 基于扇形边界的交流***谐波阻抗的边界构建方法及装置 |
CN110518577A (zh) * | 2019-08-06 | 2019-11-29 | 南方电网科学研究院有限责任公司 | 交流***谐波阻抗的多边形边界构建方法及装置 |
CN110474325B (zh) * | 2019-08-06 | 2020-12-18 | 南方电网科学研究院有限责任公司 | 基于扇形边界的交流***谐波阻抗的边界构建方法及装置 |
CN110456159A (zh) * | 2019-08-21 | 2019-11-15 | 四川大学 | 基于修正独立随机矢量的***侧谐波阻抗估计方法及*** |
CN110456159B (zh) * | 2019-08-21 | 2020-03-17 | 四川大学 | 基于修正独立随机矢量的***侧谐波阻抗估计方法及*** |
CN111077372A (zh) * | 2019-12-24 | 2020-04-28 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种交直流输电***谐波阻抗测定方法及*** |
CN115032458A (zh) * | 2022-06-02 | 2022-09-09 | 北京妙微科技有限公司 | 一种基于电网***的谐波阻抗估计方法及计算机存储介质 |
CN115032458B (zh) * | 2022-06-02 | 2023-08-11 | 北京妙微科技有限公司 | 一种基于电网***的谐波阻抗估计方法及计算机存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102998535B (zh) | 2014-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102998535A (zh) | 一种基于极大似然估计理论的***谐波阻抗计算方法 | |
Jain et al. | Fast harmonic estimation of stationary and time-varying signals using EA-AWNN | |
CN102749521B (zh) | 一种电力***谐波阻抗计算方法 | |
Jo et al. | On improving distortion power quality index in distributed power grids | |
CN106682407B (zh) | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 | |
CN105548718A (zh) | 一种基于混合整体最小二乘法的***谐波阻抗计算方法 | |
CN106154040B (zh) | 一种用于计算电网单点等值时等效谐波阻抗方法 | |
Gerber et al. | Energy and power quality measurement for electrical distribution in AC and DC microgrid buildings | |
Brunoro et al. | Modeling of loads dependent on harmonic voltages | |
Guo et al. | Global-sensitivity-based theoretical analysis and fast prediction of traveling waves with respect to fault resistance on HVDC transmission lines | |
CN107064744A (zh) | 一种谐波源定位方法 | |
CN108763802A (zh) | 计及电力传输线的电力电子换流***电磁干扰估算方法 | |
CN107359612A (zh) | 一种电能质量对配电网能耗影响的综合评估方法 | |
CN105548809A (zh) | 一种电力***谐波责任的确定方法 | |
CN109541305A (zh) | 一种谐波责任划分模型及谐波责任计算方法 | |
Sheng et al. | Online measurement-based estimation of dynamic system state matrix in ambient conditions | |
Zhang et al. | A non-intrusive identification method of harmonic source loads for industrial users | |
CN102280877B (zh) | 一种多量测断面的电力***不良支路参数辨识方法 | |
Kwon et al. | The modeling and harmonic coupling analysis of multiple-parallel connected inverter using harmonic state space (HSS) | |
CN109494743A (zh) | 基于矩阵重构的电力***电网谐波状态的估计方法 | |
Daisy et al. | A robust transient and sustainable faults location approach for AC microgrid based on voltage and current difference measurements | |
CN104201708A (zh) | 基于改进矩阵束算法的多馈入直流***附加阻尼控制器的布点方法 | |
CN104240036B (zh) | 一种基于临界***等效阻抗的暂态电压稳定量化评估方法 | |
Farzanehrafat et al. | Review of power quality state estimation | |
Zhang et al. | Research on power system harmonic state estimation |
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: 20140903 Termination date: 20191210 |
|
CF01 | Termination of patent right due to non-payment of annual fee |