三相有功无功电流检测方法
1. p-q检测法
(1) 三相系统中有功无功定义
自从赤木泰文老哥提出了三相瞬时无功功率理论,检测界迎来了春天。不同于单相系统中有功无功的定义,三相系统的有功无功定义应转化到αβ平面上。下面我将简要分析一下。
[
e
α
e
β
]
=
C
32
[
e
a
e
b
e
c
]
\left[\begin{matrix} e_\alpha\\e_\beta\end{matrix}\right] =C_{32}\left[\begin{matrix} e_a\\e_b\\e_c\end{matrix}\right]
[eαeβ]=C32⎣⎡eaebec⎦⎤
[
i
α
i
β
]
=
C
32
[
i
a
i
b
i
c
]
\left[\begin{matrix} i_\alpha\\ i_\beta \end{matrix} \right] =C_{32} \left[\begin{matrix} i_a\\i_b\\i_c \end{matrix}\right]
[iαiβ]=C32⎣⎡iaibic⎦⎤
式中:
C
32
=
2
3
[
1
−
1
2
−
1
2
0
3
2
−
3
2
]
C_{32}= \sqrt {\frac 23 } \left [\begin{matrix} 1&-\frac12&-\frac12\\0&{\frac{ \sqrt{3}}{2}}&{-\frac{ \sqrt{3}}{2}}\end{matrix}\right]
C32=32[10−2123−21−23]
在
α
β
\alpha\beta
αβ平面上,将
e
α
,
e
β
e_\alpha,e_\beta
eα,eβ合成为电压矢量
e
e
e,将
i
α
,
i
β
i_\alpha,i_\beta
iα,iβ合成为电流矢量
i
i
i,即
e
=
e
α
+
e
β
=
e
∠
ϕ
e
i
=
e
α
+
i
β
=
e
∠
ϕ
i
e=e_\alpha+e_\beta=e\angle\phi_e \\ i=e_\alpha+i_\beta=e\angle\phi_i
e=eα+eβ=e∠ϕei=eα+iβ=e∠ϕi
其矢量分解图如图1所示:
其中
ϕ
=
ϕ
e
−
ϕ
i
\phi=\phi_e-\phi_i
ϕ=ϕe−ϕi,满足一下表达式:
i
p
=
i
cos
ϕ
i
q
=
i
sin
ϕ
p
=
e
i
p
q
=
e
i
q
C
p
q
=
[
e
α
e
β
e
β
−
e
α
]
[
p
q
]
=
[
e
α
e
β
e
β
−
e
α
]
[
i
α
i
β
]
=
C
p
q
[
i
α
i
β
]
i_p=i\cos{\phi} \\ i_q=i\sin{\phi} \\ p=ei_p \\ q=ei_q \\ C_{pq}=\left[\begin{matrix} e_\alpha&e_\beta\\ e_\beta&-e_\alpha \end{matrix}\right] \\ \left[\begin{matrix} p\\ q \end{matrix}\right]=\left[\begin{matrix} e_\alpha&e_\beta\\ e_\beta&-e_\alpha \end{matrix}\right] \left[\begin{matrix} i_\alpha\\ i_\beta \end{matrix}\right] =C_{pq}\left[\begin{matrix} i_\alpha\\ i_\beta \end{matrix}\right]
ip=icosϕiq=isinϕp=eipq=eiqCpq=[eαeβeβ−eα][pq]=[eαeβeβ−eα][iαiβ]=Cpq[iαiβ]
这里的
p
p
p和
q
q
q就是三相电路的瞬时有功功率和无功功率。
(2) p-q检测原理
假定电压三相对称不含谐波,而电网电流三相对称含谐波,波形表达式如下:
{
e
a
=
2
E
sin
w
t
e
b
=
2
E
sin
(
w
t
−
2
π
3
)
e
c
=
2
E
sin
(
w
t
+
2
π
3
)
{
i
a
=
2
I
n
∑
n
=
1
n
sin
(
n
w
t
+
ϕ
n
)
i
b
=
2
I
n
∑
n
=
1
n
sin
[
n
(
w
t
−
2
π
3
)
+
ϕ
n
]
i
c
=
2
I
n
∑
n
=
1
n
sin
[
n
(
w
t
+
2
π
3
)
+
ϕ
n
]
\begin{cases} e_a=\sqrt{2}E\sin{wt} \\ e_b=\sqrt2E\sin({wt-{\frac{2\pi}{3}}}) \\ e_c=\sqrt2E\sin({wt+{\frac{2\pi}{3}}}) \end{cases} \\ \begin{cases} i_a=\sqrt{2}I_n \sum_{n=1}^n \sin(nwt+\phi_n ) \\ i_b=\sqrt{2}I_n \sum_{n=1}^n \sin[n(wt-\frac{2\pi}{3})+\phi_n ] \\ i_c=\sqrt{2}I_n \sum_{n=1}^n \sin[n(wt+\frac{2\pi}{3})+\phi_n ]\\ \end{cases}
⎩⎪⎨⎪⎧ea=2Esinwteb=2Esin(wt−32π)ec=2Esin(wt+32π)⎩⎪⎨⎪⎧ia=2In∑n=1nsin(nwt+ϕn)ib=2In∑n=1nsin[n(wt−32π)+ϕn]ic=2In∑n=1nsin[n(wt+32π)+ϕn]
由
[
p
q
]
=
[
e
α
e
β
e
β
−
e
α
]
[
i
α
i
β
]
\left[\begin{matrix}p\\q\end{matrix}\right]=\left[\begin{matrix}e_\alpha&e_\beta\\e_\beta&-e_\alpha\end{matrix}\right]\left[\begin{matrix}i_\alpha\\i_\beta\end{matrix}\right]
[pq]=[eαeβeβ−eα][iαiβ]可求出
p
p
p 和
q
q
q,其含有直流分量和交流分量,并非一个正弦波。经过低通滤波器可获取其直流分量,也就是有功功率
p
‾
\overline p
p和无功功率
q
‾
\overline q
q。其表达式如下:
[
p
‾
q
‾
]
=
3
E
I
1
[
cos
(
−
ϕ
1
)
sin
(
−
ϕ
1
)
]
C
23
=
2
3
[
1
0
−
1
2
3
2
−
1
2
−
3
2
]
[
i
a
f
i
b
f
i
c
f
]
=
C
23
C
p
q
−
1
=
2
I
1
[
s
i
n
(
w
t
+
ϕ
1
)
s
i
n
(
w
t
−
2
π
3
+
ϕ
1
)
s
i
n
(
w
t
+
2
π
3
+
ϕ
1
)
]
\left[\begin{matrix} \overline p\\ \overline q \end{matrix}\right] =3EI_1 \left[\begin{matrix} \cos(-\phi_1) \\ \sin(-\phi_1) \end{matrix}\right]\\ C_{23}=\sqrt\frac{2}{3} \left[\begin{matrix} 1&0\\ -\frac{1}{2} & \frac{\sqrt3}{2}\\ -\frac{1}{2} & -\frac{\sqrt3}{2}\\ \end{matrix}\right] \\ \left[\begin{matrix} i_{af}\\ i_{bf}\\ i_{cf} \end{matrix}\right]=C_{23}C_{pq}^{-1} =\sqrt2I_1 \left[\begin{matrix} sin(wt+\phi_1)\\ sin({wt-{\frac{2\pi}{3}}+\phi_1}) \\ sin({wt+{\frac{2\pi}{3}}+\phi_1}) \end{matrix}\right]
[pq]=3EI1[cos(−ϕ1)sin(−ϕ1)]C23=32⎣⎢⎡1−21−21023−23⎦⎥⎤⎣⎡iafibficf⎦⎤=C23Cpq−1=2I1⎣⎡sin(wt+ϕ1)sin(wt−32π+ϕ1)sin(wt+32π+ϕ1)⎦⎤第三个表达式就是电网电流的基波分量。此方法当三相不对称或电压波形畸变时不能够准确检测出有功电流和无功电流 。
2. d-q检测法
(1) Clark变换和Park变换
众所周知,交流电机三相对称的静止绕组 A 、B 、C ,通以三相平衡的正弦电流时,所产生的合成磁动势是旋转磁动势F,它在空间呈正弦分布,以同步转速
ω
\omega
ω (即电流的角频率)顺着 A-B-C 的相序旋转。abc轴等效模型如图所示。两相绕组
α
\alpha
α 和
β
\beta
β ,它们在空间
9
0
∘
90^\circ
90∘,通以时间上互差
9
0
∘
90^\circ
90∘的两相平衡交流电流,也产生旋转磁动势 F 。 当两个旋转磁动势大小和转速都相等时,即认为图b的两相绕组与图a的三相绕组等效。
α
β
\alpha\beta
αβ 轴等效模型如图所示。两个匝数相等且互相垂直的绕组 d 和 q,其中分别通以直流电流 id 和iq,产生合成磁动势 F ,其位置相对于绕组来说是固定的。如果让包含两个绕组在内的整个铁心以同步转速旋转,则磁动势 F 自然也随之旋转起来,成为旋转磁动势。
在三相静止绕组A、B、C和两相静止绕组
α
β
\alpha\beta
αβ之间的变换,或称三相静止坐标系和两相静止坐标系间的变换,简称 3/2 变换或Clark变换。下面是等功率变换,即变换前后功率不变,如果要求等幅值变换,只需将
2
3
\sqrt\frac{2}{3}
32修改为
2
3
\frac{2}{3}
32即可。
下图是
α
β
\alpha\beta
αβ轴和
d
q
dq
dq轴的关系。但
α
β
\alpha\beta
αβ轴是静止的,
α
\alpha
α 轴与
d
d
d 轴的夹角
θ
\theta
θ 随时间而变化,因此 is 在
α
β
\alpha\beta
αβ轴上的分量的长短也随时间变化,相当于绕组交流磁动势的瞬时值。由图可见,
i
α
i_\alpha
iα ,
i
β
i_\beta
iβ 和
i
d
i_d
id、
i
q
i_q
iq 之间存在下列关系
[
i
d
i
q
]
=
[
cos
θ
sin
θ
−
sin
θ
cos
θ
]
[
i
α
i
β
]
\left[\begin{matrix} i_d\\ i_q \end{matrix}\right]= \left[\begin{matrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{matrix}\right] \left[\begin{matrix} i_\alpha\\ i_\beta \end{matrix}\right]
[idiq]=[cosθ−sinθsinθcosθ][iαiβ]
此处特别注意,该变换(又称Park变换)中
cos
θ
\cos\theta
cosθ,
sin
θ
\sin\theta
sinθ以及正负号,是不固定的,他们和d轴q轴的选择有关,所以使用前请注意其dq轴的定义。
(2) dq检测法
dq检测法的核心就是Park变换和Park反变换。Park变换分为等量变换和等功率变换,其中等量变换前后电压电流矢量在变换的前后是相等的。并且abc-dq的变换矩阵等于变换矩阵abc-αβ和变换矩阵αβ-dq的乘积,其变换矩阵为
等功率变换前后三相功率不变,其变换矩阵为
下面以等量变换为例,说明一下dq检测法。假设三相电压分别是
U
s
i
n
(
ω
t
+
ϕ
)
Usin(\omega t+\phi)
Usin(ωt+ϕ),
U
s
i
n
(
ω
t
−
12
0
∘
+
ϕ
)
Usin(\omega t-120^\circ+\phi)
Usin(ωt−120∘+ϕ),
U
s
i
n
(
ω
t
+
12
0
∘
+
ϕ
)
Usin(\omega t+120^\circ+\phi)
Usin(ωt+120∘+ϕ)。三相电压经过dq变换可得其有功分量和无功分量。
[
U
d
U
q
0
]
=
[
U
cos
ϕ
U
sin
ϕ
0
]
\left[\begin{matrix} U_d\\ U_q \\0\end{matrix}\right]=\left[\begin{matrix} U\cos\phi\\ U\sin\phi \\0\end{matrix}\right]
⎣⎡UdUq0⎦⎤=⎣⎡UcosϕUsinϕ0⎦⎤
三相电流的变换与三相电压一致。三相的有功功率
p
=
u
d
i
d
+
u
q
i
q
p=u_di_d+u_qi_q
p=udid+uqiq。d-q检测法在电压畸变的情况下仍能完成有效检测,更加先进,但是dq轴的选择和等幅值等功率选择直接影响了Park变换的矩阵(可以判断一下MATLAB中的Park变换矩阵),这一点请务必注意。simulink的dq轴相对于书上滞后了90度。所以simulink的q轴是书上的d轴,而其d轴则是书上的q轴取反,而且是等幅值变换,并且
θ
\theta
θ的定义也变了,在仿真时一定要小心。如下图。