CN103777076A - 三相四线制***任意次谐波分量和无功电流检测方法 - Google Patents

三相四线制***任意次谐波分量和无功电流检测方法 Download PDF

Info

Publication number
CN103777076A
CN103777076A CN201410037919.8A CN201410037919A CN103777076A CN 103777076 A CN103777076 A CN 103777076A CN 201410037919 A CN201410037919 A CN 201410037919A CN 103777076 A CN103777076 A CN 103777076A
Authority
CN
China
Prior art keywords
omega
phase
sequence
theta
sin
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
CN201410037919.8A
Other languages
English (en)
Other versions
CN103777076B (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.)
Tsinghua University
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
North China Electric Power University
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Tsinghua University
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
North China Electric Power University
Electric Power Research Institute of State Grid Jiangsu Electric Power 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 Tsinghua University, State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, North China Electric Power University, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN201410037919.8A priority Critical patent/CN103777076B/zh
Publication of CN103777076A publication Critical patent/CN103777076A/zh
Application granted granted Critical
Publication of CN103777076B publication Critical patent/CN103777076B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Measurement Of Current Or Voltage (AREA)

Abstract

本发明公开了一种三相四线制***任意次谐波分量和无功电流的检测方法,其特征在于:包括以下步骤:步骤(1):进行锁相,获取基波正序电压的相位θ+(ω0t);步骤(2):进行n次谐波正序分量的检测;步骤(3):进行n次谐波负序分量的检测;步骤(4):进行n次谐波零序分量的检测;步骤(5):进行无功电流的检测。本发明提供的一种三相四线制***任意次谐波和无功电流的检测方法,基于n次谐波的坐标变换及低通滤波技术(LPF),并考虑了零序电流检测的特殊性;采用一种闭环的锁相方法,保证锁相的准确性,从而为坐标变换提供可靠的相位,该方法物理概念清晰,实现方式简单,适用范围广,检测准确性较高且实用性强。

Description

三相四线制***任意次谐波分量和无功电流检测方法
技术领域
本发明涉及一种三相四线制***任意次谐波和无功电流的检测方法,可以在三相电压不对称非正弦的情况下,检测出三相四线制***中的任意次谐波的正序、负序、零序分量及无功电流,该方法同样适用于三相三线制***的检测,属于电力技术应用领域。
背景技术
由于非线性负载和感应电机的大量应用,工业用户普遍存在谐波污染严重、功率因数较低等问题。三相四线制有源滤波器(APF)可以进行用于谐波治理和无功补偿,而谐波及无功电流检测的准确性是APF能够进行谐波和无功精确补偿的前提和关键。目前,基于瞬时无功功率理论的p-q法和ip-iq法应用最为普遍,效果也较好。然而,p-q法在***电压非正弦情况下存在较大的检测误差,ip-iq法不能检测电流中的无功分量且仅适用于三相三线制***。
发明内容
本发明所要解决的技术问题是,提供一种可以在***电压不对称、非正弦的情况下,检测三相电流中存在的任意次谐波的正序、负序、零序分量及无功电流,该方法基于n次谐波的坐标变换及低通滤波技术(LPF),并考虑了零序电流检测的特殊性;该方法采用一种闭环的锁相方法,可以在三相电压不对称、非正弦的情况下,保证锁相的准确性,从而为坐标变换提供可靠的相位,该方法物理概念清晰,实现方式简单,适用范围广,检测准确性较高且实用性强。
为解决上述技术问题,本发明采用的技术方案为:
三相四线制***任意次谐波分量和无功电流检测方法,包括以下步骤:
步骤(1):通过锁相环获取电网电压基波正序电压分量的相位θ+(ω0t):三相电网电压ua、ub、uc经过CLARK变换和PARK变换后,得到d轴电压分量ud和q轴电压分量uq;基波正序电压分量在d轴和q轴上表现为直流量,而基波负序电压分量和其它次谐波分量在d轴和q轴上表现为交流量;uq经过低通滤波器(LPF)滤除交流量,得到直流量
Figure BDA0000462208470000011
与目标值0的误差信号err经过PI控制器输出相位的修正量θerr;θerr与θ0之和即为锁相环的输出θ*,θ*作为PARK变换的参考相位,
Figure BDA0000462208470000021
的大小反映相位误差的大小,若
Figure BDA0000462208470000022
PI控制器将不断地修正θerr的值,达到稳态后
Figure BDA0000462208470000023
等于0,锁相环的输出θ*等于电网电压基波正序电压分量的相位θ+(ω0t)。
所述步骤(1)具体包括如下步骤:
步骤(1-1):电网电压的合成矢量沿A-B-C方向以额定角速度旋转;从t=0时刻开始旋转的角度θ0即为基波正序分量的相位,即
Figure BDA0000462208470000026
式中:uα和uβ分别为电网电压合成矢量在α轴和β轴上的投影;
步骤(1-2):三相电网电压ua、ub、uc经过CLARK变换矩阵为C32的CLARK变换,即得到电网电压合成矢量在α轴和β轴上的投影uα和uβ,如公式(1-2):
u α u β = C 32 u a u b u c - - - ( 1 - 2 )
其中: C 32 = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 - - - ( 1 - 3 )
uα和uβ分别为电网电压合成矢量
Figure BDA00004622084700000211
在α轴和β轴上的投影,C32为CLARK变换矩阵;
步骤(1-3),假设电网电压的基波正序分量为
Figure BDA00004622084700000212
则:
u a 1 + u b 1 + u c 1 + = U + sin ( θ + ) U + sin ( θ + - 2 π / 3 ) U + sin ( θ + + 2 π / 3 ) - - - ( 1 - 4 )
式中:U为电网基波正序相电压的最大值;
步骤(1-4),电网电压基波正序分量
Figure BDA00004622084700000214
经过变换矩阵为C32的CLARK变换得到的
Figure BDA00004622084700000215
Figure BDA00004622084700000216
再经过PARK变换矩阵为Cdq的PARK变换,得到
Figure BDA00004622084700000217
Figure BDA0000462208470000031
如公式(1-5)、公式(1-6):
u α + u β + = C 32 u a 1 + u b 1 + u c 1 + - - - ( 1 - 5 )
u d + u q + = C dq u α + u β + - - - ( 1 - 6 )
其中, C dq = sin θ * - cos θ * cos θ * sin θ * - - - ( 1 - 7 )
θ*为锁相环输出的相位,
Figure BDA0000462208470000035
Figure BDA0000462208470000036
为在dq坐标系中电网电压基波正序电压分量d轴和q轴直流量;
将公式(1-4)和公式(1-5)代入公式(1-6)得:
u d + u q + = U + cos ( θ + - θ * ) sin ( θ + - θ * ) - - - ( 1 - 8 )
由公式(1-8)可以看出,若锁相环得到的相位θ*与电网电压基波正序分量的相位θ+相等,则
Figure BDA0000462208470000038
采用闭环的控制方式,通过控制
Figure BDA0000462208470000039
对公式(1-1)得到的θ0进行修正,获得基波正序电压的准确相位θ+,所述θ+为一个时间参数函数,即为θ+(ω0t)。
步骤(2):将abc静止坐标转换为以角速度nω0沿a-b-c方向旋转的dq坐标系,通过正序dq变换矩阵、负序dq变换矩阵、傅立叶变换和基波正序dq变换矩阵进行n次谐波正序分量、负序分量、零序分量或者无功电流的检测。
1、进行n次谐波正序分量的检测:
三相电流ia、ib、ic经过正序dq变换矩阵为
Figure BDA00004622084700000310
的n次谐波的正序dq变换,得到d轴和q轴上的正序电流分量idn+、iqn+,如式(1):
i dn + i qn + = T abc - dq n + i a i b i c - - - ( 1 )
其中, T abc - dq n + = 2 3 sin n ω 0 t sin ( n ω 0 t - 2 π / 3 ) sin ( n ω 0 t + 2 π / 3 ) cos n ω 0 t cos ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) - - - ( 2 )
idn+和iqn+经过低通滤波器滤除交流分量得到相应的直流量
Figure BDA0000462208470000042
Figure BDA0000462208470000044
再经过正序dq反变换矩阵为
Figure BDA0000462208470000045
的n次谐波的正序dq反变换,即得到n次谐波的正序电流分量ian+、ibn+、icn+,如式(3):
i an + i bn + i cn + = T dq - abc n + i ‾ dn + i ‾ qn + - - - ( 3 )
其中, T dq - abc n + = 2 3 sin n ω 0 t cos n ω 0 t sin ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) sin ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) - - - ( 4 ) .
2、进行n次谐波负序分量的检测:
三相电流ia、ib、ic经过负序dq变换矩阵为的n次谐波的负序dq变换,得到d轴和q轴上的负序电流分量idn-和iqn-,如式(5):
i dn - i qn - = T abc - dq n - i a i b i c - - - ( 5 )
其中, T abc - dq n - = 2 3 sin n ω 0 t sin ( n ω 0 t + 2 π / 3 ) sin ( n ω 0 t - 2 π / 3 ) cos n ω 0 t cos ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) - - - ( 6 )
idn-和iqn-经过低通滤波器滤除交流分量得到相应的直流量
Figure BDA00004622084700000411
Figure BDA00004622084700000412
Figure BDA00004622084700000413
和iqn-再经过负序dq反变换矩阵为
Figure BDA00004622084700000414
的n次谐波的负序dq反变换,即得到n次谐波的负序电流分量ian-、ibn-、icn-,如式(7):
i an - i bn - i cn - = T dq - abc n - i ‾ dn - i ‾ qn - - - - ( 7 )
其中, T dq - abc n - = 2 3 sin n ω 0 t cos n ω 0 t sin ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) sin ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) - - - ( 8 ) .
3、进行n次谐波零序分量的检测:abc三相电流ia、ib、ic的正序、负序分别对称、零序分量相等,即(ia+ib+ic)/3=ia0=ib0=ic0,等于abc三相中任意相上所有次谐波零序分量之和;令待检测信号x(t)表示待检测信号零序电流分量,即x(t)=ia0,检测信号x(t)与参考正弦信号sinnω0t和余弦信号cosnω0t作乘法,得到xd(t)和xq(t),如以下式(9)、式(10)所示:
x d ( t ) = x ( t ) · sin n ω 0 t = x ( t ) · e jn ω 0 t - e - jn ω 0 t 2 - - - ( 9 )
x q ( t ) = x ( t ) · cos n ω 0 t = x ( t ) · e jn ω 0 t - e - jn ω 0 t 2 - - - ( 10 )
对式(9)和式(10)做傅立叶变换,将xd(t)和xq(t)变换为频域信号xd(ω)和xd(ω),如式(11)和式(12):
X d ( ω ) = 1 2 j ( - X ( ω + n ω 0 ) + X ( ω - n ω 0 ) ) - - - ( 11 )
X q ( ω ) = 1 2 ( X ( ω + n ω 0 ) + X ( ω - n ω 0 ) ) - - - ( 12 )
由上述式(11)和式(12)可知:零序电流与正弦和余弦信号相乘,相当于将零序电流的频谱向左和向右平移nω0,且幅值减少为原来的一半,则零序电流的n次谐波分量转换为直流量和2n次分量,其余次谐波分量仍为交流量;再用低通滤波器滤除xd(t)和xq(t)中的交流量,可以得到xd(t)和xq(t)中的直流量
Figure BDA0000462208470000055
Figure BDA0000462208470000056
Figure BDA0000462208470000057
与相应的正弦信号sinnω0t和余弦信号cosnω0t相乘,即得到n次谐波的零序分量ian0、ibn0、icn0,如式(13):
i an 0 = i bn 0 = i cn 0 = 2 x ‾ d ( t ) · sin n ω 0 t + 2 x ‾ q ( t ) · cos n ω 0 t - - - ( 13 ) .
4、进行无功电流的检测:三相电流ia、ib、ic经过基波正序dq变换矩阵为
Figure BDA0000462208470000059
的基波正序dq变换,得到d轴和q轴上的电流id1+和iq1+,如式(14):
i d 1 + i q 1 + = T abc - dq 1 + i a i b i c - - - ( 14 )
其中, T abc - dq 1 + = 2 3 sin θ + sin ( θ + - 2 π / 3 ) sin ( θ + + 2 π / 3 ) cos θ + cos ( θ + - 2 π / 3 ) cos ( θ + + 2 π / 3 ) - - - ( 15 )
dq坐标系的d轴与基波正序电压的合成矢量重合,因此d轴分量对应有功分量,q轴分量对应无功分量;iq1+经过低通滤波器滤除交流量,得到q轴的直流分量即为无功分量
Figure BDA0000462208470000063
Figure BDA0000462208470000064
Figure BDA0000462208470000065
再经过基波正序dq反变换矩阵为基波正序dq反变换即得到无功电流ia_rec、ib_rec、ic_rec,见式(16):
i a _ rec i b _ rec i c _ rec = T dq - abc 1 + i ‾ d 1 + i ‾ q 1 + - - - ( 16 )
其中, T dq - abc 1 + = 2 3 sin θ + cos θ + sin ( θ + - 2 π / 3 ) cos ( θ + - 2 π / 3 ) sin ( θ + + 2 π / 3 ) cos ( θ + + 2 π / 3 ) - - - ( 17 ) .
n次谐波零序分量检测方法还能用于单相***中任意次谐波的检测。
与现有技术相比,本发明提供的三相四线制***任意次谐波和无功电流的检测方法,基于n次谐波的坐标变换及低通滤波技术(LPF),并考虑了零序电流检测的特殊性;该方法采用一种闭环的锁相方法,可以在三相电压不对称、非正弦的情况下,保证锁相的准确性,从而为坐标变换提供可靠的相位,该方法物理概念清晰,实现方式简单,适用范围广,检测准确性较高且实用性强。
附图说明
图1为本发明的锁相环控制的结构示意图;
图2为本发明电网电压的合成矢量
Figure BDA0000462208470000069
沿A-B-C方向以额定角速度旋转示意图;
图3为本发明的正序、负序、零序分量检测原理图;
图4为本发明的无功电流检测原理图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
三相四线制***任意次谐波分量和无功电流检测方法,包括以下步骤:
步骤(1):如图1所示,通过锁相环获取电网电压基波正序电压分量的相位θ+(ω0t):三相电网电压ua、ub、uc经过CLARK变换和PARK变换后,得到d轴电压分量ud和q轴电压分量uq;基波正序电压分量在d轴和q轴上表现为直流量,而基波负序电压分量和其它次谐波分量在d轴和q轴上表现为交流量;uq经过低通滤波器(LPF)滤除交流量,得到q轴直流量
Figure BDA0000462208470000071
与目标值0的误差信号err经过PI控制器输出相位的修正量θerr;θerr与θ0之和即为锁相环的输出θ*,θ*作为PARK变换的参考相位,
Figure BDA0000462208470000072
的大小反映相位误差的大小,若
Figure BDA0000462208470000073
PI控制器将不断地修正θerr的值,达到稳态后当
Figure BDA0000462208470000074
等于0,锁相环的输出θ*等于电网电压基波正序电压分量的相位θ+(ω0t)。
所述步骤(1)具体包括如下步骤:
步骤(1-1):如图2所示,当电网电压三相对称且无畸变时,即电网电压中仅含有基波正序分量,电网电压的合成矢量
Figure BDA0000462208470000075
沿A-B-C方向以额定角速度旋转;
Figure BDA0000462208470000076
从t=0时刻开始旋转的角度θ0即为基波正序分量的相位,即
Figure BDA0000462208470000077
式中:uα和uβ分别为电网电压合成矢量
Figure BDA0000462208470000078
在α轴和β轴上的投影;
步骤(1-2):三相电网电压ua、ub、uc经过CLARK变换矩阵为C32的CLARK变换,即得到电网电压合成矢量
Figure BDA0000462208470000079
在α轴和β轴上的投影uα和uβ,如公式(2):
u α u β = C 32 u a u b u c - - - ( 1 - 2 )
其中: C 32 = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 - - - ( 1 - 3 )
uα和uβ分别为电网电压合成矢量
Figure BDA0000462208470000082
在α轴和β轴上的投影,C32为CLARK变换矩阵;
步骤(1-3),假设电网电压的基波正序分量为
Figure BDA0000462208470000083
则:
u a 1 + u b 1 + u c 1 + = U + sin ( θ + ) U + sin ( θ + - 2 π / 3 ) U + sin ( θ + + 2 π / 3 ) - - - ( 1 - 4 )
式中:U+为电网基波正序相电压的最大值。
步骤(1-4),电网电压基波正序分量
Figure BDA0000462208470000085
经过变换矩阵为C32的CLARK变换得到的
Figure BDA0000462208470000086
Figure BDA0000462208470000087
再经过PARK变换矩阵为Cdq的PARK变换,得到
Figure BDA0000462208470000088
Figure BDA0000462208470000089
如公式(1-5)、公式(1-6):
u α + u β + = C 32 u a 1 + u b 1 + u c 1 + - - - ( 1 - 5 )
u d + u q + = C dq u α + u β + - - - ( 1 - 6 )
其中, C dq = sin θ * - cos θ * cos θ * sin θ * - - - ( 1 - 7 )
θ*为锁相环输出的相位;
将公式(1-4)和公式(1-5)代入公式(6)得:
u d + u q + = U + cos ( θ + - θ * ) sin ( θ + - θ * ) - - - ( 1 - 8 )
由公式(1-8)可以看出,若锁相环得到的相位θ*与电网电压基波正序分量的相位θ+相等,则
Figure BDA00004622084700000814
采用闭环的控制方式,通过控制
Figure BDA00004622084700000815
对公式(1-1)得到的θ0进行修正,获得基波正序电压的准确相位θ+(ω0t)。
步骤(2):将abc静止坐标转换为以角速度nω0沿a-b-c方向旋转的dq坐标系,通过正序dq变换矩阵、负序dq变换矩阵、傅立叶变换和基波正序dq变换矩阵进行n次谐波正序分量、负序分量、零序分量或者无功电流的检测。
如图3和图4所示,步骤(2)n次谐波正序分量、负序分量、零序分量或者无功电流的检测具体步骤如下。
1、进行n次谐波正序分量的检测:
三相电流ia、ib、ic经过正序dq变换矩阵为
Figure BDA0000462208470000091
的n次谐波的正序dq变换,得到d轴和q轴上的正序电流分量idn+、iqn+,如式(1):
i dn + i qn + = T abc - dq n + i a i b i c - - - ( 1 )
其中, T abc - dq n + = 2 3 sin n ω 0 t sin ( n ω 0 t - 2 π / 3 ) sin ( n ω 0 t + 2 π / 3 ) cos n ω 0 t cos ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) - - - ( 2 )
idn+和iqn+经过低通滤波器滤除交流分量得到相应的直流量
Figure BDA0000462208470000094
Figure BDA0000462208470000095
Figure BDA0000462208470000096
Figure BDA0000462208470000097
再经过正序dq反变换矩阵为
Figure BDA0000462208470000098
的n次谐波的正序dq反变换,即得到n次谐波的正序电流分量ian+、ibn+、icn+,如式(3):
i an + i bn + i cn + = T dq - abc n + i ‾ dn + i ‾ qn + - - - ( 3 )
其中, T dq - abc n + = 2 3 sin n ω 0 t cos n ω 0 t sin ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) sin ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) - - - ( 4 ) .
2、进行n次谐波负序分量的检测:
三相电流ia、ib、ic经过负序dq变换矩阵为
Figure BDA00004622084700000911
的n次谐波的负序dq变换,得到d轴和q轴上的负序电流分量idn-和iqn-,如式(5):
i dn - i qn - = T abc - dq n - i a i b i c - - - ( 5 )
其中, T abc - dq n - = 2 3 sin n ω 0 t sin ( n ω 0 t + 2 π / 3 ) sin ( n ω 0 t - 2 π / 3 ) cos n ω 0 t cos ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) - - - ( 6 )
idn-和iqn-经过低通滤波器滤除交流分量得到相应的直流量
Figure BDA0000462208470000101
Figure BDA0000462208470000102
Figure BDA0000462208470000103
和iqn-再经过负序dq反变换矩阵为
Figure BDA0000462208470000104
的n次谐波的负序dq反变换,即得到n次谐波的负序电流分量ian-、ibn-、icn-,如式(7):
i an - i bn - i cn - = T dq - abc n - i ‾ dn - i ‾ qn - - - - ( 7 )
其中, T dq - abc n - = 2 3 sin n ω 0 t cos n ω 0 t sin ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) sin ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) - - - ( 8 ) .
3、进行n次谐波零序分量的检测:
abc三相电流ia、ib、ic的正序、负序分别对称、零序分量相等,即(ia+ib+ic)/3=ia0=ib0=ic0,等于abc三相中任意相上所有次谐波零序分量之和;令待检测信号x(t)表示待检测信号零序电流分量,即x(t)=ia0,检测信号x(t)与参考正弦信号sinnω0t和余弦信号cosnω0t作乘法,得到xd(t)和xq(t),如以下式(9)、式(10)所示:
x d ( t ) = x ( t ) · sin n ω 0 t = x ( t ) · e jn ω 0 t - e - jn ω 0 t 2 - - - ( 9 )
x q ( t ) = x ( t ) · cos n ω 0 t = x ( t ) · e jn ω 0 t - e - jn ω 0 t 2 - - - ( 10 )
对式(9)和式(10)做傅立叶变换,将xd(t)和xq(t)变换为频域信号xd(ω)和xd(ω),,如式(11)和式(12):
X d ( ω ) = 1 2 j ( - X ( ω + n ω 0 ) + X ( ω - n ω 0 ) ) - - - ( 11 )
X q ( ω ) = 1 2 ( X ( ω + n ω 0 ) + X ( ω - n ω 0 ) ) - - - ( 12 )
由上述式(11)和式(12)可知:零序电流与正弦和余弦信号相乘,相当于将零序电流的频谱向左和向右平移nω0,且幅值减少为原来的一半,则零序电流的n次谐波分量转换为直流量和2n次分量,其余次谐波分量仍为交流量;再用低通滤波器滤除xd(t)和xq(t)中的交流量,可以得到xd(t)和xq(t)中的直流量
Figure BDA0000462208470000111
Figure BDA0000462208470000112
Figure BDA0000462208470000113
与相应的正弦信号sinnω0t和余弦信号cosnω0t相乘,即得到n次谐波的零序分量ian0、ibn0、icn0,如式(13):
i an 0 = i bn 0 = i cn 0 = 2 x ‾ d ( t ) · sin n ω 0 t + 2 x ‾ q ( t ) · cos n ω 0 t - - - ( 13 ) .
4、进行无功电流的检测:三相电流ia、ib、ic经过基波正序dq变换矩阵为的基波正序dq变换,得到d轴和q轴上的电流id1+和iq1+,如式(14):
i d 1 + i q 1 + = T abc - dq 1 + i a i b i c - - - ( 14 )
其中, T abc - dq 1 + = 2 3 sin θ + sin ( θ + - 2 π / 3 ) sin ( θ + + 2 π / 3 ) cos θ + cos ( θ + - 2 π / 3 ) cos ( θ + + 2 π / 3 ) - - - ( 15 )
dq坐标系的d轴与基波正序电压的合成矢量重合,因此d轴分量对应有功分量,q轴分量对应无功分量;iq1+经过低通滤波器滤除交流量,得到q轴的直流分量即为无功分量
Figure BDA0000462208470000118
Figure BDA0000462208470000119
Figure BDA00004622084700001110
再经过基波正序dq反变换矩阵为
Figure BDA00004622084700001111
基波正序dq反变换即得到无功电流ia_rec、ib_rec、ic_rec,见式(16):
i a _ rec i b _ rec i c _ rec = T dq - abc 1 + i ‾ d 1 + i ‾ q 1 + - - - ( 16 )
其中, T dq - abc 1 + = 2 3 sin θ + cos θ + sin ( θ + - 2 π / 3 ) cos ( θ + - 2 π / 3 ) sin ( θ + + 2 π / 3 ) cos ( θ + + 2 π / 3 ) - - - ( 17 ) .
n次谐波零序分量检测方法还能用于单相***中任意次谐波的检测。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (6)

1.三相四线制***任意次谐波分量和无功电流检测方法,其特征在于:包括以下步骤:
步骤(1):通过锁相环获取电网电压基波正序电压分量的相位θ+(ω0t):三相电网电压ua、ub、uc经过CLARK变换和PARK变换后,得到d轴电压分量ud和q轴电压分量uq;基波正序电压分量在d轴和q轴上表现为直流量,而基波负序电压分量和其它次谐波分量在d轴和q轴上表现为交流量;uq经过低通滤波器滤除交流量,得到q轴直流量
Figure FDA0000462208460000011
与目标值0的误差信号err经过PI控制器输出相位的修正量θerr;θerr与θ0之和即为锁相环的输出θ*,θ*作为PARK变换的参考相位,的大小反映相位误差的大小,若
Figure FDA0000462208460000013
PI控制器将不断地修正θerr的值,达到稳态后
Figure FDA0000462208460000014
等于0,锁相环的输出θ*等于电网电压基波正序电压分量的相位θ+(ω0t);
步骤(2):将abc静止坐标转换为以角速度nω0沿a-b-c方向旋转的dq坐标系,通过正序dq变换矩阵、负序dq变换矩阵、傅立叶变换和基波正序dq变换矩阵进行n次谐波正序分量、负序分量、零序分量或者无功电流的检测。
2.根据权利要求1所述的三相四线制***任意次谐波分量和无功电流检测方法,其特征在于:所述步骤(1)具体包括如下步骤:
步骤(1-1):电网电压的合成矢量
Figure FDA0000462208460000015
沿a-b-c方向以额定角速度旋转;
Figure FDA0000462208460000016
从t=0时刻开始旋转的角度θ0即为基波正序分量的相位,即
Figure FDA0000462208460000017
式中:uα和uβ分别为电网电压合成矢量
Figure FDA0000462208460000018
在α轴和β轴上的投影;
步骤(1-2):三相电网电压ua、ub、uc经过CLARK变换矩阵为C32的CLARK变换,即得到电网电压合成矢量
Figure FDA0000462208460000019
在α轴和β轴上的投影uα和uβ,如公式(1-2):
u α u β = C 32 u a u b u c - - - ( 1 - 2 )
其中: C 32 = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 - - - ( 1 - 3 )
uα和uβ分别为电网电压合成矢量
Figure FDA0000462208460000022
在α轴和β轴上的投影,C32为CLARK变换矩阵;
步骤(1-3),假设电网电压的基波正序分量为
Figure FDA0000462208460000023
则:
u a 1 + u b 1 + u c 1 + = U + sin ( θ + ) U + sin ( θ + - 2 π / 3 ) U + sin ( θ + + 2 π / 3 ) - - - ( 1 - 4 )
式中:U+为电网基波正序相电压的最大值;
步骤(1-4),电网电压基波正序分量
Figure FDA0000462208460000025
经过变换矩阵为C32的CLARK变换得到的
Figure FDA0000462208460000026
Figure FDA0000462208460000027
再经过PARK变换矩阵为Cdq的PARK变换,得到
Figure FDA0000462208460000028
如公式(1-5)、公式(1-6):
u α + u β + = C 32 u a 1 + u b 1 + u c 1 + - - - ( 1 - 5 )
u d + u q + = C dq u α + u β + - - - ( 1 - 6 )
其中, C dq = sin θ * - cos θ * cos θ * sin θ * - - - ( 1 - 7 )
θ*为锁相环输出的相位,
Figure FDA00004622084600000213
为在dq坐标系中电网电压基波正序电压分量d轴和q轴直流量;
将公式(1-4)和公式(1-5)代入公式(1-6)得:
u d + u q + = U + cos ( θ + - θ * ) sin ( θ + - θ * ) - - - ( 1 - 8 )
由公式(1-8)可以看出,若锁相环得到的相位θ*与电网电压基波正序分量的相位θ+相等,则
Figure FDA00004622084600000216
采用闭环的控制方式,通过控制
Figure FDA00004622084600000217
对公式(1-1)得到的θ0进行修正,获得基波正序电压的准确相位θ+,所述θ+为一个时间参数函数,为θ+(ω0t)。
3.根据权利要求1所述的三相四线制***任意次谐波分量和无功电流检测方法,其特征在于,所述n次谐波正序分量检测方法包括以下步骤:
三相电流ia、ib、ic经过正序dq变换矩阵为
Figure FDA0000462208460000031
的n次谐波的正序dq变换,得到d轴和q轴上的正序电流分量idn+、iqn+,如式(1):
i dn + i qn + = T abc - dq n + i a i b i c - - - ( 1 )
其中, T abc - dq n + = 2 3 sin n ω 0 t sin ( n ω 0 t - 2 π / 3 ) sin ( n ω 0 t + 2 π / 3 ) cos n ω 0 t cos ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) - - - ( 2 )
idn+和iqn+经过低通滤波器滤除交流分量得到相应的直流量
Figure FDA0000462208460000034
Figure FDA0000462208460000036
Figure FDA0000462208460000037
再经过正序dq反变换矩阵为的n次谐波的正序dq反变换,即得到n次谐波的正序电流分量ian+、ibn+、icn+,如式(3):
i an + i bn + i cn + = T dq - abc n + i ‾ dn + i ‾ qn + - - - ( 3 )
其中, T dq - abc n + = 2 3 sin n ω 0 t cos n ω 0 t sin ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) sin ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) - - - ( 4 ) .
4.根据权利要求1所述的三相四线制***任意次谐波分量和无功电流检测方法,其特征在于,所述n次谐波负序分量的检测包括以下步骤:
三相电流ia、ib、ic经过负序dq变换矩阵为
Figure FDA00004622084600000311
的n次谐波的负序dq变换,得到d轴和q轴上的负序电流分量idn-和iqn-,如式(5):
i dn - i qn - = T abc - dq n - i a i b i c - - - ( 5 )
其中, T abc - dq n - = 2 3 sin n ω 0 t sin ( n ω 0 t + 2 π / 3 ) sin ( n ω 0 t - 2 π / 3 ) cos n ω 0 t cos ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) - - - ( 6 )
idn-和iqn-经过低通滤波器滤除交流分量得到相应的直流量
Figure FDA0000462208460000041
Figure FDA0000462208460000042
Figure FDA0000462208460000043
和iqn-再经过负序dq反变换矩阵为
Figure FDA0000462208460000044
的n次谐波的负序dq反变换,即得到n次谐波的负序电流分量ian-、ibn-、icn-,如式(7):
i an - i bn - i cn - = T dq - abc n - i ‾ dn - i ‾ qn - - - - ( 7 )
其中, T dq - abc n - = 2 3 sin n ω 0 t cos n ω 0 t sin ( n ω 0 t + 2 π / 3 ) cos ( n ω 0 t + 2 π / 3 ) sin ( n ω 0 t - 2 π / 3 ) cos ( n ω 0 t - 2 π / 3 ) - - - ( 8 ) .
5.根据权利要求1所述的三相四线制***任意次谐波分量和无功电流检测方法,其特征在于,所述n次谐波零序分量的检测包括以下步骤:abc三相电流ia、ib、ic的正序、负序分别对称、零序分量相等,即(ia+ib+ic)/3=ia0=ib0=ic0,等于abc三相中任意相上所有次谐波零序分量之和;令待检测信号x(t)表示待检测信号零序电流分量,即x(t)=ia0,检测信号x(t)与参考正弦信号sinnω0t和余弦信号cosnω0t作乘法,得到xd(t)和xq(t),如以下式(9)、式(10)所示:
x d ( t ) = x ( t ) · sin n ω 0 t = x ( t ) · e jn ω 0 t - e - jn ω 0 t 2 - - - ( 9 )
x q ( t ) = x ( t ) · cos n ω 0 t = x ( t ) · e jn ω 0 t - e - jn ω 0 t 2 - - - ( 10 )
对式(9)和式(10)做傅立叶变换,将xd(t)和xq(t)变换为频域信号xd(ω)和xq(ω),如式(11)和式(12):
X d ( ω ) = 1 2 j ( - X ( ω + n ω 0 ) + X ( ω - n ω 0 ) ) - - - ( 11 )
X q ( ω ) = 1 2 ( X ( ω + n ω 0 ) + X ( ω - n ω 0 ) ) - - - ( 12 )
由上述式(11)和式(12)可知:零序电流与正弦和余弦信号相乘,相当于将零序电流的频谱向左和向右平移nω0,且幅值减少为原来的一半,则零序电流的n次谐波分量转换为直流量和2n次分量,其余次谐波分量仍为交流量;再用低通滤波器滤除xd(t)和xq(t)中的交流量,可以得到xd(t)和xq(t)中的直流量
Figure FDA0000462208460000051
Figure FDA0000462208460000053
与相应的正弦信号sinnω0t和余弦信号cosnω0t相乘,即得到n次谐波的零序分量ian0、ibn0、icn0,如式(13):
i an 0 = i bn 0 = i cn 0 = 2 x ‾ d ( t ) · sin n ω 0 t + 2 x ‾ q ( t ) · cos n ω 0 t - - - ( 13 ) .
6.根据权利要求1所述的三相四线制***任意次谐波分量和无功电流检测方法,其特征在于,所述无功电流的检测方法包括以下步骤:三相电流ia、ib、ic经过基波正序dq变换矩阵为的基波正序dq变换,得到d轴和q轴上的电流id1+和iq1+,如式(14):
i d 1 + i q 1 + = T abc - dq 1 + i a i b i c - - - ( 14 )
其中, T abc - dq 1 + = 2 3 sin θ + sin ( θ + - 2 π / 3 ) sin ( θ + + 2 π / 3 ) cos θ + cos ( θ + - 2 π / 3 ) cos ( θ + + 2 π / 3 ) - - - ( 15 )
dq坐标系的d轴与基波正序电压的合成矢量重合,因此d轴分量对应有功分量,q轴分量对应无功分量;iq1+经过低通滤波器滤除交流量,得到q轴的直流分量即为无功分量
Figure FDA0000462208460000058
Figure FDA0000462208460000059
Figure FDA00004622084600000510
再经过基波正序dq反变换矩阵为
Figure FDA00004622084600000511
基波正序dq反变换即得到无功电流ia_rec、ib_rec、ic_rec,见式(16):
i a _ rec i b _ rec i c _ rec = T dq - abc 1 + i ‾ d 1 + i ‾ q 1 + - - - ( 16 )
其中, T dq - abc 1 + = 2 3 sin θ + cos θ + sin ( θ + - 2 π / 3 ) cos ( θ + - 2 π / 3 ) sin ( θ + + 2 π / 3 ) cos ( θ + + 2 π / 3 ) - - - ( 17 ) .
CN201410037919.8A 2014-01-26 2014-01-26 三相四线制***任意次谐波分量和无功电流检测方法 Active CN103777076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410037919.8A CN103777076B (zh) 2014-01-26 2014-01-26 三相四线制***任意次谐波分量和无功电流检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410037919.8A CN103777076B (zh) 2014-01-26 2014-01-26 三相四线制***任意次谐波分量和无功电流检测方法

Publications (2)

Publication Number Publication Date
CN103777076A true CN103777076A (zh) 2014-05-07
CN103777076B CN103777076B (zh) 2016-11-16

Family

ID=50569554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410037919.8A Active CN103777076B (zh) 2014-01-26 2014-01-26 三相四线制***任意次谐波分量和无功电流检测方法

Country Status (1)

Country Link
CN (1) CN103777076B (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104184148A (zh) * 2014-08-14 2014-12-03 国家电网公司 一种同步旋转坐标系中谐波电流分次控制方法
CN104181390A (zh) * 2014-08-18 2014-12-03 信元瑞电气有限公司 一种基于三相四线制零序电流分离的谐波检测方法
CN104199506A (zh) * 2014-08-25 2014-12-10 哈尔滨同为电气股份有限公司 一种应用于无功补偿控制***中的设备三相电流反馈单元
CN104280593A (zh) * 2014-10-21 2015-01-14 国家电网公司 基于瞬时无功功率理论的谐波及无功电流快速检测方法
CN104502705A (zh) * 2014-12-04 2015-04-08 西安理工大学 适用于电网电压畸变且不平衡的无锁相环旋转矢量检测法
CN104655909A (zh) * 2015-03-16 2015-05-27 福建农林大学 一种具有对称分量法计算功能的电路
CN106027038A (zh) * 2016-05-13 2016-10-12 电子科技大学 一种基于延时信号消除法的改进型三相锁相环技术
CN106959385A (zh) * 2017-04-21 2017-07-18 燕山大学 电网电压不平衡时基于二倍频锁相的谐波电流检测法
CN107677884A (zh) * 2017-08-01 2018-02-09 中国电力科学研究院 一种用于对三相三线制***的电能进行计量的方法及***
CN108448579A (zh) * 2018-03-12 2018-08-24 天津大学 一种适用于含谐波电压的单相***容性无功电流提取方法
US10215784B1 (en) 2017-12-05 2019-02-26 Industrial Technology Research Institute Measuring apparatus including phase locked loop and measuring method thereof
CN110168384A (zh) * 2017-09-08 2019-08-23 深圳欣锐科技股份有限公司 一种谐波检测方法及有源电力滤波器
CN110308326A (zh) * 2019-07-15 2019-10-08 国网山西省电力公司电力科学研究院 一种可提高开环测相算法抗噪性能的方法
CN112098721A (zh) * 2020-08-13 2020-12-18 闽南理工学院 基于状态空间模型的谐波检测方法
CN112600235A (zh) * 2020-12-08 2021-04-02 广东电网有限责任公司阳江供电局 一种柔性直流换流器等效阻抗的优化控制方法和装置
CN112698095A (zh) * 2021-01-22 2021-04-23 中南大学 一种基波和任意次谐波的高精度检测方法
CN113964848A (zh) * 2021-12-20 2022-01-21 昆山帝森华途工业物联网科技有限公司 一种apf控制***的无功电流提取模块

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05264605A (ja) * 1992-03-23 1993-10-12 Fuji Electric Co Ltd 三相電流または電圧の正相・逆相成分検出回路
CN101777912A (zh) * 2010-01-15 2010-07-14 浙江大学 一种非固定采样频率的软件锁相环实现方法
CN101893651A (zh) * 2010-06-18 2010-11-24 上海理工大学 供电***的正序、负序、无功与谐波电流检测方法
CN101900751A (zh) * 2009-06-01 2010-12-01 雷富特节能科技(上海)有限公司 能量回馈单元中实时电压矢量的检测方法
CN102290815A (zh) * 2011-08-22 2011-12-21 湖北三环发展股份有限公司 一种基于耦合电感的有源电力滤波装置
CN102522749A (zh) * 2011-11-10 2012-06-27 中冶华天工程技术有限公司 H桥级联型有源滤波器及其控制方法
CN102590646A (zh) * 2012-02-15 2012-07-18 上海大学 一种基于相序滤波器的正序、负序与谐波电流检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05264605A (ja) * 1992-03-23 1993-10-12 Fuji Electric Co Ltd 三相電流または電圧の正相・逆相成分検出回路
CN101900751A (zh) * 2009-06-01 2010-12-01 雷富特节能科技(上海)有限公司 能量回馈单元中实时电压矢量的检测方法
CN101777912A (zh) * 2010-01-15 2010-07-14 浙江大学 一种非固定采样频率的软件锁相环实现方法
CN101893651A (zh) * 2010-06-18 2010-11-24 上海理工大学 供电***的正序、负序、无功与谐波电流检测方法
CN102290815A (zh) * 2011-08-22 2011-12-21 湖北三环发展股份有限公司 一种基于耦合电感的有源电力滤波装置
CN102522749A (zh) * 2011-11-10 2012-06-27 中冶华天工程技术有限公司 H桥级联型有源滤波器及其控制方法
CN102590646A (zh) * 2012-02-15 2012-07-18 上海大学 一种基于相序滤波器的正序、负序与谐波电流检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨柳 等: "三相四线制***任意次谐波电流的检测新方法", 《中国电机工程学报》 *
杨浩 等: "三相电压不对称***的谐波和无功电流检测", 《计算机仿真》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104184148B (zh) * 2014-08-14 2017-01-25 国家电网公司 一种同步旋转坐标系中谐波电流分次控制方法
CN104184148A (zh) * 2014-08-14 2014-12-03 国家电网公司 一种同步旋转坐标系中谐波电流分次控制方法
CN104181390A (zh) * 2014-08-18 2014-12-03 信元瑞电气有限公司 一种基于三相四线制零序电流分离的谐波检测方法
CN104199506A (zh) * 2014-08-25 2014-12-10 哈尔滨同为电气股份有限公司 一种应用于无功补偿控制***中的设备三相电流反馈单元
CN104280593A (zh) * 2014-10-21 2015-01-14 国家电网公司 基于瞬时无功功率理论的谐波及无功电流快速检测方法
CN104502705B (zh) * 2014-12-04 2017-09-29 西安理工大学 适用于电网电压畸变且不平衡的无锁相环旋转矢量检测法
CN104502705A (zh) * 2014-12-04 2015-04-08 西安理工大学 适用于电网电压畸变且不平衡的无锁相环旋转矢量检测法
CN104655909A (zh) * 2015-03-16 2015-05-27 福建农林大学 一种具有对称分量法计算功能的电路
CN104655909B (zh) * 2015-03-16 2017-12-08 福建农林大学 一种具有对称分量法计算功能的电路
CN106027038A (zh) * 2016-05-13 2016-10-12 电子科技大学 一种基于延时信号消除法的改进型三相锁相环技术
CN106027038B (zh) * 2016-05-13 2019-03-05 电子科技大学 一种基于延时信号消除法的改进型三相锁相环算法
CN106959385A (zh) * 2017-04-21 2017-07-18 燕山大学 电网电压不平衡时基于二倍频锁相的谐波电流检测法
CN107677884A (zh) * 2017-08-01 2018-02-09 中国电力科学研究院 一种用于对三相三线制***的电能进行计量的方法及***
CN110168384A (zh) * 2017-09-08 2019-08-23 深圳欣锐科技股份有限公司 一种谐波检测方法及有源电力滤波器
US10215784B1 (en) 2017-12-05 2019-02-26 Industrial Technology Research Institute Measuring apparatus including phase locked loop and measuring method thereof
CN108448579A (zh) * 2018-03-12 2018-08-24 天津大学 一种适用于含谐波电压的单相***容性无功电流提取方法
CN108448579B (zh) * 2018-03-12 2021-06-08 天津大学 一种适用于含谐波电压的单相***容性无功电流提取方法
CN110308326A (zh) * 2019-07-15 2019-10-08 国网山西省电力公司电力科学研究院 一种可提高开环测相算法抗噪性能的方法
CN112098721A (zh) * 2020-08-13 2020-12-18 闽南理工学院 基于状态空间模型的谐波检测方法
CN112098721B (zh) * 2020-08-13 2023-05-26 闽南理工学院 基于状态空间模型的谐波检测方法
CN112600235A (zh) * 2020-12-08 2021-04-02 广东电网有限责任公司阳江供电局 一种柔性直流换流器等效阻抗的优化控制方法和装置
CN112698095A (zh) * 2021-01-22 2021-04-23 中南大学 一种基波和任意次谐波的高精度检测方法
CN113964848A (zh) * 2021-12-20 2022-01-21 昆山帝森华途工业物联网科技有限公司 一种apf控制***的无功电流提取模块
CN113964848B (zh) * 2021-12-20 2022-05-31 昆山帝森华途工业物联网科技有限公司 一种apf控制***的无功电流提取模块

Also Published As

Publication number Publication date
CN103777076B (zh) 2016-11-16

Similar Documents

Publication Publication Date Title
CN103777076A (zh) 三相四线制***任意次谐波分量和无功电流检测方法
CN101893652B (zh) 一种基于电压矢量空间变换的谐波和无功电流检测方法
CN101893651B (zh) 供电***的正序、负序、无功与谐波电流检测方法
Zaveri et al. Comparison of control strategies for DSTATCOM in three-phase, four-wire distribution system for power quality improvement under various source voltage and load conditions
CN102735938A (zh) 一种电网电压基波正序相角的快速检测方法
CN103269086B (zh) 光伏并网逆变器低电压穿越控制的正负序分量分离方法
CN104659813B (zh) 一种快速谐波环流抑制的多逆变器并联控制方法
CN102590646A (zh) 一种基于相序滤波器的正序、负序与谐波电流检测方法
CN105823921A (zh) 一种基于瞬时空间电压矢量定向的补偿电流检测方法
CN106532749B (zh) 一种微电网不平衡功率和谐波电压补偿***及其应用
CN113300381B (zh) 不平衡工况下链式statcom的控制方法
CN104201679B (zh) 抑制微网中电流谐波与三相不平衡的电流型逆变控制策略
CN103399202A (zh) 可变速抽水蓄能机组控制***的锁相方法
CN105896554A (zh) 一种三相四线制svg装置输出限幅方法
CN104502705A (zh) 适用于电网电压畸变且不平衡的无锁相环旋转矢量检测法
CN103267897A (zh) 一种基于反Park变换的三相锁相环
CN102385003A (zh) 检测电力***中电压信号或电流信号频率和相位的方法
CN103472301A (zh) 一种电网电压正序、负序分量的提取方法和***
CN105021872A (zh) 一种电网电压畸变不对称状态下电流各分量有功成分和无功成分的检测方法
Suru et al. The synchronous fundamental dq frame theory implementation and adaptation for the active filtering
CN106026094B (zh) 一种应用于有源滤波器的不平衡控制算法
CN111208340A (zh) 一种基于傅里叶变换的单相基波无功电流精确检测方法
CN103235184B (zh) 一种双馈风力发电变流器电网电压矢量角度检测算法
He et al. Research on digital phase locked method in PWM rectifier
CN111200287B (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
C14 Grant of patent or utility model
GR01 Patent grant