发明名称 一种基于线性滤波器的三相四线制谐波电流检测方法
摘要 本发明公开一种基于线性滤波器的三相四线制谐波电流检测方法,针对三相四线制系统的结构特点,按照自适应陷波器(ANF)的思想,依据最小方差原则与梯度下降方法,经旋转变换得到二维线性正弦滤波器,分别从实测电流中分离各相基波电流,再利用改进的瞬时对称分量计算方法,得到对称的基波正序分量,与原实测电流运算后抽取出负载电流的谐波分量及基波负序分量及零序分量。与现有技术相比,本发明提供的基于线性滤波器的三相四线制谐波电流检测方法,方法简单,运算量小,无需进行复杂的旋转变换和三角函数的计算,实时性好,检测精度高,能够在不对称非线性负载和非理想的电压情况下检测三相四线制系统中的谐波电流、无功电流、中线电流。
申请公布号 CN102253266B 申请公布日期 2013.06.05
申请号 CN201110109158.9 申请日期 2011.04.28
申请人 西华大学 发明人 邱晓初
分类号 G01R19/00(2006.01)I 主分类号 G01R19/00(2006.01)I
代理机构 北京集佳知识产权代理有限公司 11227 代理人 曹志霞;李赞坚
主权项 1.一种基于线性滤波器的三相四线制谐波电流检测方法,其特征在于,包括如下步骤: (a)利用三相四线制的结构特点,将其各相看成三个独立的单相,不平衡三相负载电流i<sub>a</sub>,i<sub>b</sub>,i<sub>c</sub>,各使用一个二维线性滤波器抽取出基波电流;设电网的基波频率为f,其角频率为ω,滤波器输入的电力信号为i(t),则有: <img file="FDA00002778494200011.GIF" wi="447" he="64" /><img file="FDA00002778494200012.GIF" wi="1888" he="93" />将状态变量选择为x(t)=[x<sub>1</sub>(t) x<sub>2</sub>(t)]<sup>T</sup>:<img file="FDA00002778494200013.GIF" wi="294" he="53" /><img file="FDA00002778494200014.GIF" wi="695" he="62" />则i(t)的估计值<img file="FDA00002778494200015.GIF" wi="93" he="82" />可表示为:<img file="FDA00002778494200016.GIF" wi="1127" he="70" />信号的估计误差等于:<img file="FDA00002778494200017.GIF" wi="1361" he="76" />根据最小方差思想,定义代价函数为:<img file="FDA00002778494200018.GIF" wi="955" he="67" />按照梯度下降法,有: <img file="FDA00002778494200019.GIF" wi="1954" he="111" />矩阵Ν调节x(t)收敛到最优值的速度以及算法的稳定性,令μ<sub>1</sub>,μ<sub>2</sub>>0,选择: <img file="FDA000027784942000110.GIF" wi="1925" he="139" />由此,可以得到周期系数线性动力系统:<img file="FDA000027784942000111.GIF" wi="867" he="60" />其中:<img file="FDA000027784942000112.GIF" wi="434" he="65" />x=[x<sub>1</sub>(t) x<sub>2</sub>(t)]<sup>T</sup><img file="FDA000027784942000113.GIF" wi="975" he="170" /><img file="FDA000027784942000114.GIF" wi="359" he="166" />作旋转变换,令 <img file="FDA00002778494200021.GIF" wi="1866" he="148" />则由式(1)和式(3)知,i(t)的估计值<img file="FDA00002778494200022.GIF" wi="82" he="73" />为:<img file="FDA00002778494200023.GIF" wi="1867" he="66" />将式(9)代入式(8),有: <img file="FDA00002778494200024.GIF" wi="1065" he="160" /><img file="FDA00002778494200025.GIF" wi="1694" he="161" />令μ<sub>1</sub>=μ<sub>2</sub>=μ,η=ω,得到二维线性滤波算法: <img file="FDA00002778494200026.GIF" wi="1865" he="149" />其中:y(t)=[y<sub>1</sub>(t) y<sub>2</sub>(t)]<sup>T</sup>为状态变量;μ,η为滤波器的参数,均大于零的实常数;y<sub>1</sub>(t)为滤波器输出,它表示输入信号中所含的角频率为ω的正弦分量,也即输入信号的基波分量,y<sub>2</sub>(t)为与y<sub>1</sub>(t)正交且超前的同频信号; 设初始条件y<sub>1</sub>(0)=0,y<sub>2</sub>(0)=0,取i(t)为输入,y<sub>1</sub>(t)为输出,则二维线性滤波器的传递函数为 <img file="FDA00002778494200027.GIF" wi="1858" he="128" />令s=jω,有: <img file="FDA00002778494200028.GIF" wi="1831" he="139" />则其幅频特性和相频特性分别为: <img file="FDA00002778494200031.GIF" wi="1831" he="149" /><img file="FDA00002778494200032.GIF" wi="1830" he="131" />当取μ=0.5,η=ω=314rad/s时,有: <img file="FDA00002778494200033.GIF" wi="1831" he="135" /><img file="FDA00002778494200034.GIF" wi="1831" he="119" />(b)不平衡的基波电流i<sub>af</sub>,i<sub>bf</sub>,i<sub>cf</sub>通过瞬时对称分量计算,得到对称的基波正序电流分量i<sub>afp</sub>,i<sub>bfp</sub>,i<sub>cfp</sub>;设不平衡三相负载电流瞬时值为: <img file="FDA00002778494200035.GIF" wi="1757" he="244" />其中:i<sub>a</sub>(t),i<sub>b</sub>(t),i<sub>c</sub>(t)分别为三相电流瞬时值;I<sub>am</sub>,I<sub>bm</sub>,I<sub>cm</sub>分别为三相电流的幅值;<img file="FDA00002778494200036.GIF" wi="191" he="44" />分别为三相电流的初相位;设与三相电流i<sub>a</sub>(t),i<sub>b</sub>(t),i<sub>c</sub>(t)所对应的旋转相量分别为<img file="FDA00002778494200037.GIF" wi="180" he="63" />有:<img file="FDA00002778494200038.GIF" wi="1842" he="257" />由式(20)知,相量<img file="FDA00002778494200039.GIF" wi="170" he="69" />的虚部系数即为三相电流瞬时值,只要能确定出其实部,就可以确定这些旋转相量;设a相电流的瞬时值可表示为:<img file="FDA00002778494200041.GIF" wi="770" he="59" /><img file="2.GIF" wi="403" he="24" /><img file="FDA00002778494200044.GIF" wi="650" he="58" />也即: <img file="FDA00002778494200045.GIF" wi="577" he="61" /><img file="FDA00002778494200046.GIF" wi="1704" he="65" />故有: <img file="FDA00002778494200047.GIF" wi="533" he="57" /><img file="FDA00002778494200048.GIF" wi="1591" he="122" />又因为: <img file="FDA00002778494200049.GIF" wi="390" he="57" /><img file="FDA000027784942000410.GIF" wi="631" he="57" /><img file="3.GIF" wi="403" he="19" /><img file="FDA000027784942000413.GIF" wi="637" he="57" />整理后可得: <img file="FDA000027784942000414.GIF" wi="1817" he="122" />同理可得: <img file="FDA000027784942000415.GIF" wi="1825" he="135" /><img file="FDA000027784942000416.GIF" wi="1816" he="124" />将式(25)~(27)代入式(20)中,构成瞬时旋转相量,则不平衡三相负载电流正序分量值为: <img file="FDA00002778494200051.GIF" wi="1793" he="291" />其中:Im{}为求复数的虚部运算;a=e<sup>j2π/3</sup>,a<sup>2</sup>=e<sup>-j2π/3</sup>,1+a+a<sup>2</sup>=0; (c)与原负载电流运算后,即可分离出谐波电流分量i<sub>ah</sub>,i<sub>bh</sub>,i<sub>ch</sub>。 
地址 610000 四川省成都市金牛区土桥金周路999号