CN111368713A - 一种基于同步压缩小波变换的车网系统谐波时频分析方法 - Google Patents

一种基于同步压缩小波变换的车网系统谐波时频分析方法 Download PDF

Info

Publication number
CN111368713A
CN111368713A CN202010135939.4A CN202010135939A CN111368713A CN 111368713 A CN111368713 A CN 111368713A CN 202010135939 A CN202010135939 A CN 202010135939A CN 111368713 A CN111368713 A CN 111368713A
Authority
CN
China
Prior art keywords
frequency
time
wavelet
voltage
transformation
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
CN202010135939.4A
Other languages
English (en)
Other versions
CN111368713B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202010135939.4A priority Critical patent/CN111368713B/zh
Publication of CN111368713A publication Critical patent/CN111368713A/zh
Application granted granted Critical
Publication of CN111368713B publication Critical patent/CN111368713B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08Feature extraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/148Wavelet transforms
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Algebra (AREA)
  • Signal Processing (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Power Engineering (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

本发明公开了一种基于同步压缩小波变换的车网系统谐波时频分析方法,在等效牵引网和CRH5型车整流器级联运行下,通过设置不同开关频率和牵引网等效参数,在考虑数字控制电路所产生的延时下,观测网侧电流电压信号,并利用同步压缩小波变换对该信号进行分析;本发明分析结果具有良好的时频分辨率,可以有效辨识系统发生谐波不稳定问题时,谐波放大的幅值和频率;本发明实施方便、适应性较高。

Description

一种基于同步压缩小波变换的车网系统谐波时频分析方法
技术领域
本发明属于高速铁路非平稳数据的时频分析处理技术领域,具体涉及一种基于同步压缩小波变换的车网系统谐波时频分析方法。
背景技术
近年来,随着HXD和CRH系列电力机车的广泛使用,牵引系统中引入了大量的电力电子设备。由于电力电子设备的非线性特性,当机车投入到牵引网后可能会和牵引网发生交互耦合作用而给牵引系统稳定性带来新的挑战。此外,不同车型中变流器开关频率相差较大,所以采用数字控制器而产生的延时也不尽相同。
为了研究动车组变流器开关频率对车网系统的影响规律,Assefa H Y等在文献1中指出了开关频率会对车网系统稳定性产生的影响。但是该文献仅仅阐述了该现象,并没有详细分析不同开关频率下的车网电气波形的时频谱特征。文献2,3指出,延时将会导致变流器等效导纳在高频段产生负阻尼效应,进而降低系统稳定性。文献4指出,由于动车组牵引功率较大,所以变流器开关频率设置的较低,这将会产生严重的控制延时。文献5通过对谐波不稳定波形进行傅里叶变换,指出该现象是一种特征次谐波放大现象。此外,针对非平稳信号,小波变换是一种常用的时频分析工具,通过构造合适的小波基可以有效分析非周期信号和含有奇异性的信号。但是小波分析结果的分辨率往往十分依赖与小波基的选取,而小波基函数的选择也没有固定规律,往往全凭经验进行选择。综合上述研究方法:
1)研究中对车网系统中所出现的谐波不稳定问题采用傅里叶变换和小波变换的分析方法,时频分辨率受限于窗函数或者小波基的选取,采用传统的时频分析方法无法达到理想的分析效果。
2)车网系统稳定性问题是一个强耦合系统,对变流器开关频率和牵引网参数变化较为敏感。
因此,有必要采用高时频分辨率的方法引入谐波不稳定问题的分析中。同步压缩小波变换不过分依赖于小波基的选择,可以准确分析非平稳信号的幅值和频率随时间变换的规律,有效判断变流器开关频率和牵引网参数对车网系统稳定性的影响。
参考文献:
[1]:Assefa H Y,Danielsen S,Molinas M.Impact of PWM switching onmodeling of low frequency power oscillation in electrical rail vehicle[C].13th European Conference on Power Electronics and Applications.IEEE,2009:1-9.
[2]:Harnefors L,Finger R,Wang X,et al.VSC input-admittance modelingand analysis above the Nyquist frequency for passivity-based stabilityassessment[J].IEEE Transactions on Industrial Electronics,2017,64(8):6362-6370.
[3]:Harnefors L,Wang X,Yepes A G,et al.Passivity-based stabilityassessment of grid-connected VSCs—An overview[J].IEEE Journal of emergingand selected topics in Power Electronics,2015,4(1):116-125.
[4]:Mollerstedt E,Bernhardsson B.Out of control because of harmonics-an analysis of the harmonic response of an inverter locomotive[J].IEEEControl Systems Magazine,2000,20(4):70-81.
[5]:陶海东,胡海涛,朱晓娟,et al.车网耦合下的牵引供电系统谐振不稳定机理分析[J].中国电机工程学报,2019,39(8):2315-2324.
发明内容
本发明的目的是:在车网系统发生谐波不稳定现象时进行时频分析,分析动车组变流器开关频率和牵引网参数对车网系统稳定性的影响,为解决车网系统谐波不稳定问题提供新思路。为此,本发明提供了一种基于同步压缩小波变换的车网系统谐波时频分析方法。
本发明的一种基于同步压缩小波变换的车网系统谐波时频分析方法,在等效牵引网和CRH5型车整流器级联运行下,通过设置不同开关频率和牵引网等效参数,在考虑数字控制电路所产生的延时下,观测网侧电流电压信号,并利用同步压缩小波变换对该信号进行分析,具体包括以下步骤:
步骤1:给定变流器开关频率的初始值,计算当PWM调制电路为单倍频采样时数字控制器向车网系统中引入的延时;
步骤2:在仿真过程中加入控制电路中所产生的延时,观察并采集牵引网电压和电流;
步骤3:通过同步压缩小波变换对牵引网电压和电流信号进行分析,提取其时频变化规律;
步骤4:提取并追踪同步压缩变换结果的时频脊,对信号进行模态提取。
其中,步骤3具体过程为:
步骤a:选取解析的小波基,对牵引网电压和电流进行小波变换,将得到的结果记作WTf(a,b),过程如下:
Figure BDA0002397324930000031
其中f(t)为待分析的信号,a为小波基ψa,b(t)中的尺度因子,b为其平移因子;
步骤b:将得到的小波变换结果WTf(a,b)对平移因子b求偏导数,得到
Figure BDA0002397324930000032
步骤c:将得到的偏导数结果除以对信号进行小波变换和i2π的乘积,得到信号瞬时频率值:
Figure BDA0002397324930000033
步骤d:将得到的瞬时频率值在小波变换所得到的时频区域进行重新分配,以提升时频分辨率:
Figure BDA0002397324930000034
其中Tf为同步压缩小波变换矩阵,ak是小波变换系数矩阵WTf进行压缩变换的计算区间,可根据式ω(ak,b)-ωl|≤Δω/2得到;压缩变换后时频分布的结果分布在中心频率ωl的附近。
与现有技术相比,本发明的有益效果是:
1、本发明通过同步压缩小波变换可以精确分析当车网系统发生谐波不稳定时,牵引网谐振频率和幅值随时间变化的规律。
2、本发明实施方便,不需要严格选择小波基对谐波不稳定现象进行分析,该时频分析方法的适应性较高。
附图说明
图1为本发明的车网系统等效框图。
图2为本发明的车网系统等效电路。
图3为本发明车载整流器开关频率设置为1250Hz时的牵引网电压和电流波形。
图4为本发明车载整流器开关频率设置为800Hz时的牵引网电压和电流波形。
图5为本发明车载整流器开关频率设置为500Hz时的牵引网电压和电流波形。
图6为开关频率设置为1250Hz时牵引网电压的时频分析图。
图7为开关频率设置为800Hz时牵引网电压的时频分析图。
图8为开关频率设置为500Hz时牵引网电压的时频分析图。
图9为开关频率设置为1250Hz时牵引网电压的模态分解图。
图10为开关频率设置为800Hz时牵引网电压的模态分解图。
图11为开关频率设置为500Hz时牵引网电压的模态分解图。
图12为半实物平台下开关频率设置为1250Hz时牵引网电压和电流的示波器图。
图13为半实物平台下开关频率设置为800Hz时牵引网电压和电流的示波器图。
图14为半实物平台下开关频率设置为500Hz时牵引网电压和电流的示波器图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。
本实例以牵引网和CRH5型动车组级联系统为例。图1为牵引供电系统AT供电方式下,车网系统等效模型。牵引供电系统由牵引变电所和牵引网组成,电力机车由车载整流器、逆变器和电机组成。图2为车网系统等效电路图,其中Us为牵引网电压,Ls、Rs、Cs分别是牵引网等效电感、电阻和电容。Lg为车载变压器等效电感,Cd为直流侧滤波电容,Rd为机车等效负载。简化后的车载整流器的控制环路由采样环节、控制电路、PWM调制电路和延时电路所组成。其中延时主要是在PWM调制过程中所产生的。当采用单倍频调制时,延时时长和开关频率大小相同。
本实施例的同步压缩小波变换车网系统谐波不稳定时频谱分析由以下各步组成。
1、通过在Matlab/Simulink中设置不同车载变流器开关频率,并考虑数字控制器所产生的延时,观察并采集牵引网电压、电流的波形;
图3、4、5分别为开关频率为1250Hz、800Hz、500Hz时的牵引网电压和电流波形。当开关频率为1250Hz时,在第5秒考虑控制电路所产生的延时,之后,牵引网上电压和电流出现谐波放大现象。在实际运行中,当牵引网电压超过29kV时,会触发机车保护装置使机车闭锁,机车停止运行。当开关频率为800Hz时,在加入延时之后,车网系统并没有出现谐波不稳定现象,且于0.15秒之后牵引网电压和电流恢复到正常运行状态。当开关频率为500Hz时,在加入延时之后,车网系统无法稳定,电压和电流失去同步。
2、利用同步压缩小波变换不同开关频率下牵引网电压时频分布图;
图6、7、8分别是对图3、4、5中5~5.1秒间电压波形进行小波变换和同步压缩小波变换后的结果。从图6可以明显看出,同步压缩小波变换的分析结果相较于小波变换具有更好的时频分辨率。小波变换无法精确分析谐波频率以及谐波幅值,但是同步压缩小波变换却可以做到精确定位。由同步压缩小波变换分析结果可知,在第5秒加入延时后,牵引网电压出现了461Hz谐波电压幅值放大现象,这一结果也和图3相吻合。经表Ι中数据可知,461Hz也正好为牵引网并联谐振频率。故当发生谐波不稳定时,牵引网固有谐振频率下电压的幅值会被放大。
图7、8是当开关频率为800Hz、500Hz时牵引网电压波形的时频分析结果。可以看出,此时,车网系统电压并没有出现特定次谐波放大现象。
3、提取并追踪对谐波不稳定波形同步压缩变换结果的时频脊,对信号进行模态提取;
为了进一步分析发生谐波不稳定时牵引网电压的波形,对图3中电压波形进行了模态提取,分析该波形的成分组成,以及各个成分的频率和幅值随时间变化的规律。
图9为当开关频率为1250Hz时电压波形模态分解结果,其中S(t)为5~5.1秒内所截取的谐波不稳定电压信号,IMF为信号S(t)的固有模态成分。可以发现,当发生谐波不稳定时,波形是由50Hz基波和幅值不断增大的谐波成分组成。图10、11分别为当开关频率为800Hz、500Hz时电压波形模态分解结果。可以发现,当开关频率为800Hz时,电压波形是由50Hz基波和幅值在较小范围内变化的谐波成分组成,并没有出现因为谐波放大而出现的不稳定现象。而当开关频率降低为500Hz时,电压波形却是由50Hz基波和无明显变化规律的谐波成分组成。此时对比图5可以发现,由于开关频率过低,导致当考虑了延时之后,牵引网电压和电流信号失去同步。
4、在半实物平台上进行试验,观察开关频率对车网系统系统稳定性的影响;
为了验证开关频率对车网系统稳定性的影响,搭建硬件在环车网系统半实物实验平台。其中控制器装载控制电路,仿真器装载主电路,在PC端运行StarSim软件使二者形成闭环系统。图12、13、14为当开关频率设置为1250、800、500Hz时,从示波器截取牵引网电压和电流的波形。可以发现半实物平台上的实验结果和Matlab\Simulink仿真结果一致。
车网级联系统参数如表Ι所示。
表Ι车网系统模型参数
Figure BDA0002397324930000051
利用同步压缩小波变换对牵引网电压和电流信号进行分析,分析结果具有良好的时频分辨率,可以有效辨识系统发生谐波不稳定问题时,谐波放大的幅值和频率。

Claims (2)

1.一种基于同步压缩小波变换的车网系统谐波时频分析方法,其特征在于,在等效牵引网和CRH5型车整流器级联运行下,通过设置不同开关频率和牵引网等效参数,在考虑数字控制电路所产生的延时下,观测网侧电流电压信号,并利用同步压缩小波变换对该信号进行分析,具体包括以下步骤:
步骤1:给定变流器开关频率的初始值,计算当PWM调制电路为单倍频采样时数字控制器向车网系统中引入的延时;
步骤2:在仿真过程中加入控制电路中所产生的延时,观察并采集牵引网电压和电流;
步骤3:通过同步压缩小波变换对牵引网电压和电流信号进行分析,提取其时频变化规律;
步骤4:提取并追踪同步压缩变换结果的时频脊,对信号进行模态提取。
2.根据权利要求1所述的一种基于同步压缩小波变换的车网系统谐波时频分析方法,其特征在于,所述步骤3具体过程为:
步骤a:选取解析的小波基,对牵引网电压和电流进行小波变换,将得到的结果记作WTf(a,b),过程如下:
Figure FDA0002397324920000011
其中f(t)为待分析的信号,a为小波基ψa,b(t)中的尺度因子,b为其平移因子;
步骤b:将得到的小波变换结果WTf(a,b)对平移因子b求偏导数,得到
Figure FDA0002397324920000012
步骤c:将得到的偏导数结果除以对信号进行小波变换和i2π的乘积,得到信号瞬时频率值:
Figure FDA0002397324920000013
步骤d:将得到的瞬时频率值在小波变换所得到的时频区域进行重新分配,以提升时频分辨率:
Figure FDA0002397324920000014
其中Tf为同步压缩小波变换矩阵,ak是小波变换系数矩阵WTf进行压缩变换的计算区间,可根据式|ω(ak,b)-ωl|≤Δω/2得到;压缩变换后时频分布的结果分布在中心频率ωl的附近。
CN202010135939.4A 2020-03-02 2020-03-02 一种基于同步压缩小波变换的车网系统谐波时频分析方法 Expired - Fee Related CN111368713B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010135939.4A CN111368713B (zh) 2020-03-02 2020-03-02 一种基于同步压缩小波变换的车网系统谐波时频分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010135939.4A CN111368713B (zh) 2020-03-02 2020-03-02 一种基于同步压缩小波变换的车网系统谐波时频分析方法

Publications (2)

Publication Number Publication Date
CN111368713A true CN111368713A (zh) 2020-07-03
CN111368713B CN111368713B (zh) 2022-06-28

Family

ID=71206656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010135939.4A Expired - Fee Related CN111368713B (zh) 2020-03-02 2020-03-02 一种基于同步压缩小波变换的车网系统谐波时频分析方法

Country Status (1)

Country Link
CN (1) CN111368713B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070009165A1 (en) * 2005-07-08 2007-01-11 Kumar B A Modular, low cost, memory efficient, input resolution independent, frame-synchronous, video compression system using multi stage wavelet analysis and temporal signature analysis with a highly optimized hardware implementation
CN103217213A (zh) * 2013-02-21 2013-07-24 北京工业大学 基于响应信号时频联合分布特征的模态参数辨识方法
US20140089241A1 (en) * 1999-02-01 2014-03-27 Steven M. Hoffberg System and method for intermachine markup language communications
CN104820786A (zh) * 2015-05-13 2015-08-05 西安交通大学 一种瞬时加权同步挤压小波双谱分析方法
CN105699779A (zh) * 2015-12-31 2016-06-22 北京交通大学 一种级联h桥型牵引网阻抗测试谐波发生器及测试方法
CN107608935A (zh) * 2017-08-30 2018-01-19 西安交通大学 基于时间重排压缩变换的冲击类信号时频分析与重构方法
CN108009347A (zh) * 2017-11-30 2018-05-08 南京理工大学 基于同步压缩联合改进广义s变换的时频分析方法
CN108732421A (zh) * 2018-06-08 2018-11-02 中国铁路总公司 高速列车动态响应信号的瞬时频率的获取方法及装置
US20190053023A1 (en) * 2014-08-11 2019-02-14 Neeraj Jhanji Compressed sensing with machine-to-machine communication
CN110244202A (zh) * 2019-06-21 2019-09-17 三峡大学 基于同步压缩小波变换域变压器局部放电去噪方法
CN110617964A (zh) * 2019-07-29 2019-12-27 中国铁道科学研究院集团有限公司城市轨道交通中心 用于滚动轴承故障诊断的同步压缩变换阶比分析法
CN110737023A (zh) * 2019-10-17 2020-01-31 北京矿冶科技集团有限公司 一种矿用微震监测信号的处理方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140089241A1 (en) * 1999-02-01 2014-03-27 Steven M. Hoffberg System and method for intermachine markup language communications
US20070009165A1 (en) * 2005-07-08 2007-01-11 Kumar B A Modular, low cost, memory efficient, input resolution independent, frame-synchronous, video compression system using multi stage wavelet analysis and temporal signature analysis with a highly optimized hardware implementation
CN103217213A (zh) * 2013-02-21 2013-07-24 北京工业大学 基于响应信号时频联合分布特征的模态参数辨识方法
US20190053023A1 (en) * 2014-08-11 2019-02-14 Neeraj Jhanji Compressed sensing with machine-to-machine communication
CN104820786A (zh) * 2015-05-13 2015-08-05 西安交通大学 一种瞬时加权同步挤压小波双谱分析方法
CN105699779A (zh) * 2015-12-31 2016-06-22 北京交通大学 一种级联h桥型牵引网阻抗测试谐波发生器及测试方法
CN107608935A (zh) * 2017-08-30 2018-01-19 西安交通大学 基于时间重排压缩变换的冲击类信号时频分析与重构方法
CN108009347A (zh) * 2017-11-30 2018-05-08 南京理工大学 基于同步压缩联合改进广义s变换的时频分析方法
CN108732421A (zh) * 2018-06-08 2018-11-02 中国铁路总公司 高速列车动态响应信号的瞬时频率的获取方法及装置
CN110244202A (zh) * 2019-06-21 2019-09-17 三峡大学 基于同步压缩小波变换域变压器局部放电去噪方法
CN110617964A (zh) * 2019-07-29 2019-12-27 中国铁道科学研究院集团有限公司城市轨道交通中心 用于滚动轴承故障诊断的同步压缩变换阶比分析法
CN110737023A (zh) * 2019-10-17 2020-01-31 北京矿冶科技集团有限公司 一种矿用微震监测信号的处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘爽等: "基于滑膜控制的牵引网网压低频振荡抑制方法", 《电网技术》, vol. 42, no. 9, 27 July 2018 (2018-07-27), pages 2999 - 3006 *
张闯等: "基于同步压缩变换的临界折射纵波波速精确测量", 《化学仪表学报》, vol. 40, no. 7, 15 July 2019 (2019-07-15), pages 97 - 105 *

Also Published As

Publication number Publication date
CN111368713B (zh) 2022-06-28

Similar Documents

Publication Publication Date Title
Hu et al. A novel controlled frequency band impedance measurement approach for single-phase railway traction power system
CN108875710A (zh) 基于能量阈值算法的电梯门运行速度估计方法
CN105528741B (zh) 一种基于多信号特征融合的断路器状态识别方法
CN105866571B (zh) 一种基于高频谐波补偿迭代emd的暂态电能质量信号分析方法
CN107036709B (zh) 一种变电站噪声匹配分离方法
CN106484962B (zh) 一种谐振式无线输电系统稳态特性的符号分析方法
CN103675544A (zh) 基于优化算法的电力系统故障信号检测与波形识别方法
Abdelsalam et al. A generalized approach for power quality disturbances recognition based on Kalman filter
CN106501602A (zh) 一种基于滑窗频谱分离的基波参数测量方法
Wang et al. Analysis of partial discharge signal using the Hilbert-Huang transform
CN111368713B (zh) 一种基于同步压缩小波变换的车网系统谐波时频分析方法
CN112861328B (zh) 一种基于随机响应信号的发电机阻尼评估装置及方法
CN108090270B (zh) 一种基于形态学滤波和盲源分离的暂态振荡参数识别方法
Borisov et al. A computationally efficient RDFT-based reference signal generator for active compensators
CN106645948A (zh) 一种基于小波变换的自适应电网信号基波检测方法
Kumar et al. S-Transform based detection of multiple and multistage power quality disturbances
Duan et al. Differential mode conducted EMI prediction in three phase SiC inverters
Asman et al. Fault signature analysis for medium voltage overhead line based on new characteristic indices
Li et al. Analysis of harmonic resonance for locomotive and traction network interacted system considering the frequency-domain passivity properties of the digitally controlled converter
Chen et al. Extraction of common-mode impedance of an induction motor by using all-phase FFT with intermediate-frequency filtering
CN106771850A (zh) 基于级联双稳态系统的配电网故障选线方法
Ferrigno et al. Measurement of Passive $ R $, $ L $, and $ C $ Components Under Nonsinusoidal Conditions: The Solution of Some Case Studies
CN110082642A (zh) 基于全相位微分滤波的电网工况故障时刻检测方法及装置
Yu et al. Research on modal parameter identification method of power transformer winding
Ghaemi et al. Voltage notch indices determination using wavelet transform

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: 20220628

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