CN104037768B - H桥级联statcom自适应死区补偿方法 - Google Patents

H桥级联statcom自适应死区补偿方法 Download PDF

Info

Publication number
CN104037768B
CN104037768B CN201410293425.6A CN201410293425A CN104037768B CN 104037768 B CN104037768 B CN 104037768B CN 201410293425 A CN201410293425 A CN 201410293425A CN 104037768 B CN104037768 B CN 104037768B
Authority
CN
China
Prior art keywords
rsqb
lsqb
statcom
bridge cascade
voltage
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.)
Expired - Fee Related
Application number
CN201410293425.6A
Other languages
English (en)
Other versions
CN104037768A (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.)
HARBIN TONGWEI ELECTRIC CO Ltd
Original Assignee
HARBIN TONGWEI ELECTRIC CO Ltd
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 HARBIN TONGWEI ELECTRIC CO Ltd filed Critical HARBIN TONGWEI ELECTRIC CO Ltd
Priority to CN201410293425.6A priority Critical patent/CN104037768B/zh
Publication of CN104037768A publication Critical patent/CN104037768A/zh
Application granted granted Critical
Publication of CN104037768B publication Critical patent/CN104037768B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/10Flexible AC transmission systems [FACTS]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Feedback Control In General (AREA)

Abstract

H桥级联STATCOM自适应死区补偿方法,涉及一种H桥级联STATCOM自适应死区补偿方法。它是为了解决缺少直接针对H桥级联STATCOM进行死区补偿的问题。本发明根据H桥级联STATCOM的电路拓扑结构,得出三相坐标系下H桥级联STATCOM的离散数学模型;并将其转换成离散状态下的传递函数形式;采用离散扰动观测器,对死区效应引起的STATCOM输出电压与参考电压的差值进行实时的在线观测,并将观测到的差值作为补偿量,引入到H桥级联STATCOM控制系统中,实现对死区的自适应补偿。本发明适用于一种H桥级联STATCOM自适应死区补偿方法。

Description

H桥级联STATCOM自适应死区补偿方法
技术领域
本发明涉及一种H桥级联STATCOM自适应死区补偿方法。
背景技术
在工业应用中,随着电网中非线性负载的广泛使用,所产生的无功和谐波电流对电网造成严重冲击,为了保证较好的电能质量,需要对6kV和10kV电网进行无功补偿和谐波消除,以提高电网的可靠性与稳定性。H桥级联静止同步补偿器(STATCOM)利用多电平变换技术,将开关器件电压应力降低,通过级联,使整个装置可以输出高电压,并具有损耗低、响应快、储能元件体积小和输出电流谐波含量低等优点,成为动态无功补偿装置发展的重要方向,在中高压大容量场合得到了迅速的发展。
但由于H桥功率单元每个桥臂的功率管具有互补的驱动信号,为避免开关器件的直通现象,必须在互补的两路信号间插入死区。然而死区时间的设置和开关器件的非理想特性,必然造成了输出电压与参考电压之间的误差,导致功率单元输出电流的谐波含量增大,还会出现零电流箝位现象。
然而,目前现有的死区补偿方法,大多只针对非H桥级联STATCOM,鲜有直接针对H桥级联STATCOM这一新兴结构提出的方法。
发明内容
本发明是为了解决缺少直接针对H桥级联STATCOM进行死区补偿的问题,从而提供了一种H桥级联STATCOM自适应死区补偿方法。
H桥级联STATCOM自适应死区补偿方法,该方法包括以下步骤:
步骤一、根据H桥级联STATCOM的电路拓扑结构,得出三相坐标系下H桥级联STATCOM的离散数学模型;
步骤二、根据步骤一得到的H桥级联STATCOM的离散数学模型,得到离散状态下的传递函数形式;
步骤三、采用离散扰动观测器,对死区效应引起的STATCOM输出电压与参考电压的差值进行实时的在线观测,并将观测到的差值作为补偿量,引入到H桥级联STATCOM控制系统中,实现对死区的自适应补偿。
本发明的有益效果是:采用本发明对实际的H桥级联STATCOM的补偿性能进行实验验证。搭建H桥级联STATCOM实验平台,并完成网侧电压us=400V时和H桥单元直流侧电压参考值Udc=200V时的低压实验。具体实验参数为:每相N=12个H桥单元,f=50Hz,单元直流侧电容C=5600μF,连接电抗器L=10mH,等效开关频率f=1kHz,死区时间Td=6μs。控制器中的DSP选择TI公司的TMS320F28335芯片,负责无功电流检测,参考指令电流计算和直流电压控制;FPGA选择Altera公司的CycloneⅢ系列的EP3C25芯片,负责产生36路PWM脉冲控制信号,并通过光纤将脉冲触发信号送到每个H桥单元。测得的实验波形均为A相电流波形。
图3和图4为H桥级联STATCOM输出的A相补偿电流波形。其中,图3为未进行死区补偿时的电流波形,可以看出电流波形有明显畸变,出现了明显的零点电流箝位现象,并且在电流峰值处出现了较大的畸变;图4为采用死区补偿后的电流波形,电流波形正弦度得到了明显改善,畸变减小,有效的消除了零点电流箝位现象,且峰值处的波形平滑。谐波畸变THD由未补偿时的2.64%,下降到1.53%。
图5和图6为STATCOM输出的A相电压波形。其中,图5为未进行死区补偿时的输出电压波形,可以看出电压含有很大的谐波,并且在电压峰值处,阶梯波形严重畸变;图6为采用死区补偿后的输出电压波形,电压所含谐波降低,在峰值处,波形平滑,畸变减小,波形的正弦度得到了明显改善。
附图说明
图1是H桥级联STATCOM的主电路拓扑框图;
图2是本发明自适应死区补偿方法框图;
图3是未采用本发明提出的自适应死区补偿方法时的H桥级联STATCOM的输出电流波形;
图4是采用本发明提出的自适应死区补偿方法时的H桥级联STATCOM的输出电流波形;
图5是未采用本发明提出的自适应死区补偿方法时的H桥级联STATCOM的输出电压波形;
图6是采用本发明提出的自适应死区补偿方法时的H桥级联STATCOM的输出电压波形。
具体实施方式
具体实施方式一:下面结合图1和图2说明本实施方式,本实施方式所述的H桥级联STATCOM自适应死区补偿方法,该方法包括以下步骤:
步骤一、根据H桥级联STATCOM的电路拓扑结构,得出三相坐标系下H桥级联STATCOM的离散数学模型;
步骤二、根据步骤一得到的H桥级联STATCOM的离散数学模型,得到离散状态下的传递函数形式;
步骤三、采用离散扰动观测器,对死区效应引起的STATCOM输出电压与参考电压的差值进行实时的在线观测,并将观测到的差值作为补偿量,引入到H桥级联STATCOM控制系统中,实现对死区的自适应补偿。
具体实施方式二:本实施方式对具体实施方式一所述的H桥级联STATCOM自适应死区补偿方法作进一步限定,本实施方式中,步骤一中的三相坐标系下的离散数学模型为:
L di a dt = u sa - u a - Ri a L di b dt = u sb - u b - Ri b L di c dt = u sc - u c - Ri c ; - - - ( 1 )
其中,usa、usb和usc为网侧三相电压,ua、ub和uc为STATCOM输出三相电压,isa、isb和isc为网侧三相电流,ia、ib和ic为STATCOM补偿三相电流,ila、ilb和ilc为负载三相电流,Udc为直流侧电容电压参考值,C为直流侧电容,L为连接电抗器,R为等效损耗电阻;
设采样周期为Ts,将式(1)离散化并整理得:
u a ( k ) = u sa ( k ) - L T s [ i a ( k + 1 ) - i a ( k ) ] - Ri a ( k ) u b ( k ) = u sb ( k ) - L T s [ i b ( k + 1 ) - i b ( k ) ] - Ri b ( k ) u c ( k ) = u sc ( k ) - L T s [ i c ( k + 1 ) - i c ( k ) ] - Ri c ( k ) ; - - - ( 2 )
根据式(2)得出,根据k时刻采样的i(k),us(k)和k+1时刻的参考指令电流预测值,计算出k采样时刻STATCOM参考输出的指令电压u(k);根据数字系统的延时性,计算出指令电压的时刻要滞后采样时刻一拍,k时刻的采样值计算得到k+1时刻的指令电压,得到式(3):
u a * ( k + 1 ) = u sa ( k ) - L T s [ i a ( k + 1 ) - i a ( k ) ] - Ri a ( k ) u b * ( k + 1 ) = u sb ( k ) - L T s [ i b ( k + 1 ) - i b ( k ) ] - Ri b ( k ) u c * ( k + 1 ) = u sc ( k ) - L T s [ i c ( k + 1 ) - i c ( k ) ] - Ri c ( k ) ; - - - ( 3 )
令i(k+1)=i*(k),则在下一采样时刻达到电压指令值,实现无差拍控制,上标*表示指令值;电阻在模型建立的过程中忽略不计,得到式(4):
u a * ( k + 1 ) = u sa ( k ) - L T s [ i a * ( k ) - i a ( k ) ] u b * ( k + 1 ) = u sb ( k ) - L T s [ i b * ( k ) - i b ( k ) ] u c * ( k + 1 ) = u sc ( k ) - L T s [ i c * ( k ) - i c ( k ) ] . - - - ( 4 )
具体实施方式三:本实施方式对具体实施方式一所述的H桥级联STATCOM自适应死区补偿方法作进一步限定,本实施方式中,步骤二中,得到离散状态下的传递函数为:
G ( z ) = Z [ H 0 ( s ) · G ( s ) ] = Z [ H 0 ( s ) · ( I x ( s ) V x ( s ) - E x ( s ) ) ] = Z [ H 0 ( s ) · ( - 1 Ls ) ] = - T s L ( z - 1 ) ; - - - ( 5 )
其中是零阶保持器连续状态下的传递函数,Ts为采样周期,(x=a,b,c)。
上述实施方式中,零阶保持器的作用是在信号传递过程中,把采样信号转换成连续信号。
具体实施方式四:本实施方式对具体实施方式一所述的H桥级联STATCOM自适应死区补偿方法作进一步限定,本实施方式中,步骤三中,结合式(4)、式(5)与离散扰动观测器的工作原理,并重新定义各个变量,将观测到的差值作为补偿量,引入到H桥级联STATCOM控制系统中,同时考虑死区效应的影响,由式(4)可得:
U a * - Δ U ^ a = U sa - L T s [ i a * - i a ] - Δ U a U b * - Δ U ^ b = U sb - L T s [ i b * - i b ] - Δ U b U c * - Δ U ^ c = U sc - L T s [ i c * - i c ] - Δ U c - - - ( 6 )
如果观测到的差值与实际的差值相等,就可以消除死区效应对输出电压及电流的影响。
根据自适应补偿原理,利用(7)式就可以得到利用离散扰动观测器观测到的补偿量,
Δ U ^ a = Q ( z ) [ i a L n ( z - 1 ) - T s - U a * * ] Δ U ^ b = Q ( z ) [ i b L n ( z - 1 ) - T s - U b * * ] Δ U ^ c = Q ( z ) [ i c L n ( z - 1 ) - T s - U c * * ] ; - - - ( 7 )
公式(6)和公式(7)中ix,Ux△Ux和Usx(x=a,b,c)分别为H桥级联STATCOM三相参考电流,输出电流,输出电压,参考电压,扰动电压,观测的扰动电压,死区补偿后的参考电压和网侧电压;Ln为标称模型中的标称电感;
最后使得观测到的差值与实际的差值相等,进而消除死区效应对输出电压及电流的影响。

Claims (2)

1.H桥级联STATCOM自适应死区补偿方法,其特征在于:该方法包括以下步骤:
步骤一、根据H桥级联STATCOM的电路拓扑结构,得出三相坐标系下H桥级联STATCOM的离散数学模型;
步骤二、根据步骤一得到的H桥级联STATCOM的离散数学模型,得到离散状态下的传递函数形式;
步骤三、采用离散扰动观测器,对死区效应引起的STATCOM输出电压与参考电压的差值进行实时的在线观测,并将观测到的差值作为补偿量,引入到H桥级联STATCOM控制系统中,实现对死区的自适应补偿;
步骤一中的三相坐标系下的离散数学模型为:
L di a d t = u s a - u a - Ri a L di b d t = u s b - u b - Ri b L di c d t = u s c - u c - Ri c ; - - - ( 1 )
其中,usa、usb和usc为网侧三相电压,ua、ub和uc为STATCOM输出三相电压,ia、ib和ic为STATCOM补偿三相电流,L为连接电抗器,R为等效损耗电阻;
设采样周期为Ts,将式(1)离散化并整理得:
{ u a ( k ) = u s a ( k ) - L T s [ i a ( k + 1 ) - i a ( k ) ] - Ri a ( k ) u b ( k ) = u s b ( k ) - L T s [ i b ( k + 1 ) - i b ( k ) ] - Ri b ( k ) u c ( k ) = u s c ( k ) - L T s [ i c ( k + 1 ) - i c ( k ) ] - Ri c ( k ) ; - - - ( 2 )
根据式(2)得出,根据k时刻采样的i(k),us(k)和k+1时刻的参考指令电流预测值i(k+1),计算出k采样时刻STATCOM参考输出的指令电压u(k);根据数字系统的延时性,计算出指令电压的时刻要滞后采样时刻一拍,k时刻的采样值计算得到k+1时刻的指令电压,得到式(3):
{ u a * ( k + 1 ) = u s a ( k ) - L T s [ i a ( k + 1 ) - i a ( k ) ] - Ri a ( k ) u b * ( k + 1 ) = u s b ( k ) - L T s [ i b ( k + 1 ) - i b ( k ) ] - Ri b ( k ) u c * ( k + 1 ) = u s c ( k ) - L T s [ i c ( k + 1 ) - i c ( k ) ] - Ri c ( k ) ; - - - ( 3 )
令i(k+1)=i*(k),则在下一采样时刻达到电压指令值,实现无差拍控制,上标*表示指令值;电阻在模型建立的过程中忽略不计,得到式(4):
{ u a * ( k + 1 ) = u s a ( k ) - L T s [ i a * ( k ) - i a ( k ) ] u b * ( k + 1 ) = u s b ( k ) - L T s [ i b * ( k ) - i b ( k ) ] u c * ( k + 1 ) = u s c ( k ) - L T s [ i c * ( k ) - i c ( k ) ] ; - - - ( 4 )
步骤二中,得到离散状态下的传递函数为:
G ( z ) = Z [ H 0 ( s ) g G ( s ) ] = Z [ H 0 ( s ) g ( I x ( s ) V x ( s ) - E x ( s ) ) ] = Z [ H 0 ( s ) g ( - 1 L s ) ] = - T s L ( z - 1 ) ; - - - ( 5 )
其中是零阶保持器连续状态下的传递函数,Ts为采样周期,x=a,b,c。
2.根据权利要求1所述的H桥级联STATCOM自适应死区补偿方法,其特征在于:步骤三中,结合式(4)、式(5)与离散扰动观测器的工作原理,并重新定义各个变量,将观测到的差值作为补偿量,引入到H桥级联STATCOM控制系统中,同时考虑死区效应的影响,由式(4)可得:
U a * - Δ U ^ a = U s a - L T s [ i a * - i a ] - ΔU a U b * - Δ U ^ b = U s b - L T s [ i b * - i b ] - ΔU b U c * - Δ U ^ c = U s c - L T s [ i c * - i c ] - ΔU c - - - ( 6 )
如果观测到的差值与实际的差值相等,就可以消除死区效应对输出电压及电流的影响,
根据自适应补偿原理,利用(7)式就可以得到利用离散扰动观测器观测到的补偿量,
Δ U ^ a = Q ( z ) [ i a L n ( z - 1 ) - T s - U a * * ] Δ U ^ b = Q ( z ) [ i b L n ( z - 1 ) - T s - U b * * ] Δ U ^ c = Q ( z ) [ i c L n ( z - 1 ) - T s - U c * * ] ; - - - ( 7 )
公式(6)和公式(7)中ix,Ux△Ux和Usx分别为H桥级联STATCOM三相参考电流,输出电流,输出电压,参考电压,扰动电压,观测的扰动电压,死区补偿后的参考电压和网侧电压,x=a,b,c;Ln为标称模型中的标称电感;
最后使得观测到的差值与实际的差值相等,进而消除死区效应对输出电压及电流的影响。
CN201410293425.6A 2014-06-25 2014-06-25 H桥级联statcom自适应死区补偿方法 Expired - Fee Related CN104037768B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410293425.6A CN104037768B (zh) 2014-06-25 2014-06-25 H桥级联statcom自适应死区补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410293425.6A CN104037768B (zh) 2014-06-25 2014-06-25 H桥级联statcom自适应死区补偿方法

Publications (2)

Publication Number Publication Date
CN104037768A CN104037768A (zh) 2014-09-10
CN104037768B true CN104037768B (zh) 2016-03-16

Family

ID=51468432

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410293425.6A Expired - Fee Related CN104037768B (zh) 2014-06-25 2014-06-25 H桥级联statcom自适应死区补偿方法

Country Status (1)

Country Link
CN (1) CN104037768B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104201921B (zh) * 2014-09-28 2017-10-03 南京冠亚电源设备有限公司 一种逆变系统中死区补偿控制方法
CN110137965B (zh) * 2019-05-17 2020-03-13 清华大学 高压级联h桥式svg直流脉动自适应滤波方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2854892Y (zh) * 2005-12-12 2007-01-03 天津理工大学 基于dsp芯片的静止同步补偿器非线性控制装置
CN103580026A (zh) * 2013-11-11 2014-02-12 东南大学 一种具有谐波优化调节器的三角形h桥级联statcom
CN203553952U (zh) * 2013-10-29 2014-04-16 山东华天电气有限公司 集中控制的多模块静止同步补偿器

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2854892Y (zh) * 2005-12-12 2007-01-03 天津理工大学 基于dsp芯片的静止同步补偿器非线性控制装置
CN203553952U (zh) * 2013-10-29 2014-04-16 山东华天电气有限公司 集中控制的多模块静止同步补偿器
CN103580026A (zh) * 2013-11-11 2014-02-12 东南大学 一种具有谐波优化调节器的三角形h桥级联statcom

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
混合级联多电平拓扑及其在STATCOM中的应用研究;许赟;《中国博士学位论文全文数据库》;20090531;第1-86页 *

Also Published As

Publication number Publication date
CN104037768A (zh) 2014-09-10

Similar Documents

Publication Publication Date Title
CN111245010B (zh) 基于llcl型三相并网逆变器的双闭环控制方法
CN103311949B (zh) 大功率光伏逆变器的低电压穿越控制方法
CN102075108A (zh) 一种带有lcl滤波器的并网逆变器的电容电流前馈控制方法
CN107069732A (zh) 基于最小电流误差模型预测的有源滤波器谐波电流补偿方法
CN107959289A (zh) 一种有源电力滤波器改进型指定次谐波电流补偿控制方法
CN105162162A (zh) 一种基于相位扰动的分布式并网逆变器孤岛检测系统
CN103972922B (zh) 基于改进型准谐振控制加重复控制的光伏并网控制方法
CN204886199U (zh) 一种适应不平衡和非线性负载的电力电子变压器
CN103904922A (zh) 一种基于虚拟磁链定向用于电压型整流器的控制方法
CN105763094A (zh) 一种基于电压前馈和复合电流控制的逆变器控制方法
CN103441502B (zh) 并联单相h桥级联型有源电力滤波器控制装置及其方法
CN102340257A (zh) 一种lcl滤波并网逆变器的双电流环控制方法
CN104143932A (zh) 双矢量模型预测直接功率控制方法及装置
CN110429603A (zh) 六开关七电平有源电力滤波器及补偿方法
CN110350534A (zh) 一种基于七电平换流器的有源电力滤波系统及其补偿方法
CN104104251B (zh) 一种基于ssr-kdf的并网逆变器的鲁棒控制方法
CN110380416A (zh) 一种储能变流器有源阻尼谐振控制方法
CN105490297A (zh) 基于双逆变器群协调控制的微电网供电电压和电网电流谐波同步补偿方法
CN107302317B (zh) 三相五电平逆变器漏电流抑制的载波实现方法
CN105244882A (zh) 三相三线制有源电力滤波器对电网谐波电流的补偿方法
CN104037768B (zh) H桥级联statcom自适应死区补偿方法
CN102437753A (zh) 一种lcl滤波的三相pwm整流器三环控制方法
CN103457267B (zh) 三相并联有源电力滤波器的空间矢量脉宽调制控制方法
CN102510062B (zh) 有源电力滤波器的改进空间矢量单周控制方法及装置
CN106208066A (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
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Yang Xu

Inventor after: Liu Yuguang

Inventor after: He Chongfei

Inventor before: Yu Yong

Inventor before: Xu Rong

Inventor before: Xu Dianguo

Inventor before: He Chongfei

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: YU YONG XU RONG XU DIANGUO HE CHONGFEI TO: YANG XU LIU YUGUANG HE CHONGFEI

C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Yang Xu

Inventor after: Yu Yong

Inventor after: Xu Rong

Inventor after: Liu Yuguang

Inventor after: He Chongfei

Inventor before: Yang Xu

Inventor before: Liu Yuguang

Inventor before: He Chongfei

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: YANG XU LIU YUGUANG HE CHONGFEI TO: YANG XU YU YONG XU RONG LIU YUGUANG HE CHONGFEI

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

Termination date: 20200625

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