一种基于最小二乘法的虚拟电阻型有源滤波控制策略的制作方法

本发明涉及电力滤波器控制技术领域，具体涉及一种基于最小二乘法的虚拟电阻型有源滤波控制策略。

背景技术：

虚拟电阻型有源电力滤波器在吸收各次谐波电流时，既增加用户成本，又难以量化其对电网谐波治理工作的贡献，相关部门也就无法对其作出补偿。因此，尽管可降低系统电压畸变率，但仍难以调动用户以虚拟电阻型apf这种方式治理谐波的积极性。针对此问题，本专利基于虚拟电阻型有源电力滤波器提出一种新的谐波治理策略，虚拟电阻型有源电力滤波器根据并网点各次谐波电压及apf输出的各次谐波电流，计算apf吸收的各次谐波功率，并采用一定算法，实时调节等效谐波电阻，使apf吸收的各次谐波功率均达到最大，并以谐波功率的排放或吸收来界定用户的污染者或治理者角色，对污染或治理行为进行惩罚或奖励。在此策略下，如果配合有效机制，则主动参与谐波治理的用户将依据其吸收的最大谐波功率可获得最大的激励。本专利把apf吸收的谐波功率为最大值相应的各次谐波电阻的阻值定义为“最佳阻值”。

技术实现要素：

本发明的目的是为了解决现有技术中的上述缺陷，提供一种基于最小二乘法的虚拟电阻型有源滤波控制策略，通过控制逆变器输出的谐波电流，实现一个等效的虚拟谐波电阻，其阻值等于系统的等效阻抗的模值，从而实现谐波的最大功率吸收。

本发明的目的可以通过采取如下技术方案达到：

一种基于最小二乘法的虚拟电阻型有源滤波控制策略，所述虚拟电阻型有源滤波控制策略包括下列步骤：

基于谐波电压源与谐波电流源共同作用的简单电力系统，采用小扰动法对系统的等效阻抗进行观测，控制有源滤波器电流进行m次谐波小扰动，分别记录第k次谐波小扰动下，apf的网点第h次谐波电压向量的实部uih(re)k和虚部uih(im)k，第h次谐波电流向量的实部iih(re)k和虚部iih(im)k,其中k＝1,2，……，m，形成超定方程组ax＝b，

其中：

式中，eh(re)为系统h次谐波等效电势向量的实部，eh(im)为系统h次谐波等效电势向量的虚部，rh为系统h次谐波等效阻抗的实部，xh为系统h次谐波等效阻抗的虚部；

求解超定方程ax＝b的最小二乘解即正规方程组：(a^ta)x＝a^tb的解，求得系统各次谐波阻抗：zh＝rh+jxh；

让逆变器虚拟谐波电阻rih匹配系统谐波阻抗模值|zh|，即根据系统谐波电压求出逆变器的输出谐波电流指令，控制逆变器输出相应的谐波电流，从而实现最大化谐波功率吸收。

进一步地，所述的正规方程组(a^ta)x＝a^tb用对称矩阵的三角分解法求解，其过程如下：

记g＝a^ta，则g是对称矩阵，由三角分解g＝ldl^t，其中l是下三角矩阵，d是对角矩阵，正规方程可化为：ldl^tx＝a^tb，正规方程的求解可分为以下三个步骤：

1)解下三角方程组：lz＝a^tb

2)解对角方程组：dy＝z

3)解上三角方程组：l^tx＝y

其中z＝dl^tx，y＝l^tx，

经过以上步骤，即可求解出

本发明相对于现有技术具有如下的优点及效果：

本发明能快速确定各次虚拟谐波电阻的最佳阻值，从而减小调节虚拟谐波电阻过程中造成的谐波功率波动。传统的导纳调节法为基于一维搜索的自动调节法，但此方法搜索速度较慢，在调节过程中产生长时间的谐波功率波动，影响谐波功率的吸收效果。而采用本专利所提的最小二乘法计算系统等效阻抗则能解决上述问题。所述最小二乘法求系统等效阻抗，是指在确定虚拟谐波电导gh时，在逆变器原谐波参考电流的基础上叠加不同幅值相位的若干次小扰动，再用不同小扰动下的频谱形成超定方程，从而得出最小二乘解下的系统等效阻抗zh，gh即为系统等效阻抗模值|zh|的倒数。这样，就可以通过若干次的小扰动就能快速确定各次虚拟谐波电阻的最佳阻值。

附图说明

图1是谐波电压源与谐波电流源共同作用的简单电力系统；

图2是虚拟电阻型有源电力滤波器工作原理；

图3是虚拟电阻电流生成示意图；

图4是最小二乘法求虚拟电阻电流流程图；

图5是等效负载电阻突然变大时的仿真波形。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例

本实施例基于虚拟电阻型有源电力滤波器提出一种新的谐波治理策略，虚拟电阻型有源电力滤波器根据并网点各次谐波电压及apf输出的各次谐波电流，计算apf吸收的各次谐波功率，并采用一定算法，实时调节等效谐波电阻，使apf吸收的各次谐波功率均达到最大，并以谐波功率的排放或吸收来界定用户的污染者或治理者角色，对污染或治理行为进行惩罚或奖励。在此策略下，如果配合有效机制，则主动参与谐波治理的用户将依据其吸收的最大谐波功率可获得最大的激励。本发明实施例中把apf吸收的谐波功率为最大值相应的各次谐波电阻的阻值定义为“最佳阻值”。

虚拟电阻型apf同时采集并网点电压和输出电流，并对其进行谐波检测，然后分别计算各相所吸收的各次谐波功率，并对三相同次谐波功率求平均值，作为apf吸收的第h次谐波功率ph。根据ph对各次电导值gh作出调节。虚拟电阻型apf在运行过程中，不要求像一般补偿电流型的apf那样对负荷电流快速的响应，而且电力系统对于每一次gh的变化需要一定的响应时间，因此，各次谐波功率计算环节和各次导纳调节环节可按设定的时间间隔进行，在各次谐波功率更新之前，gh保留原来的值不变。

传统的导纳调节法为基于一维搜索的自动调节法，但此方法搜索速度较慢，在调节过程中产生长时间的谐波功率波动，影响谐波功率的吸收效果。而采用本专利所提的最小二乘法计算系统等效阻抗则能解决上述问题。所述最小二乘法求系统等效阻抗，是指在确定虚拟谐波电导gh时，在逆变器原谐波参考电流的基础上叠加不同幅值相位的若干次小扰动，再用不同小扰动下的频谱形成超定方程，从而得出最小二乘解下的系统等效阻抗zh，gh即为系统等效阻抗模值|zh|的倒数。这样，就可以通过若干次的小扰动就能快速确定各次谐波电阻的最佳阻值。

一、图1是谐波电压源与谐波电流源共同作用的简单电力系统，采用小扰动对系统的等效阻抗进行观测，可列出下列方程：

zh＝rh+jxh(5)

uih(re)＝iih(re)rh-iih(im)xh+eh(re)(6)

uih(im)＝iih(im)rh+iih(re)xh+eh(im)(7)

在上述方程中，为apf的网点h次谐波电压向量，uih(re)为apf的网点h次谐波电压向量的实部，uih(im)为apf的网点h次谐波电压向量的虚部，为系统h次谐波等效电势向量，eh(re)为系统h次谐波等效电势向量的实部，eh(im)为系统h次谐波等效电势向量的虚部，为apf入网h次谐波电流向量，iih(re)为apf入网h次谐波电流向量的实部，iih(im)为apf入网h次谐波电流向量的虚部，zh为系统h次谐波等效阻抗，rh为系统h次谐波等效阻抗的实部，xh为系统h次谐波等效阻抗的虚部。其中，eh(re)，eh(im)，rh，xh是待估计的恒定参数，uih(re)，uih(im)，iih(re)，iih(im)是观测量。

根据上述方程，可以得到关于m次小扰动时的矩阵方程：

ax＝b(8)

其中，

式中，uih(re)k为第k次小扰动时通过fft得到的apf入网h次谐波电流向量的实部，iih(im)k为第k次小扰动时通过fft得到的apf入网h次谐波电流向量的虚部，uih(re)k为第k次小扰动时apf的网点h次谐波电压向量的实部，uih(im)k为第k次小扰动时apf的网点h次谐波电压向量的虚部。

超定方程租ax＝b的最小二乘解的为正规方程组：(a^ta)x＝a^tb的解。用对称矩阵的三角分解法求解正规方程组(a^ta)x＝a^tb，记g＝a^ta，则g是对称矩阵，由三角分解g＝ldl^t，其中l是下三角矩阵，d是对角矩阵，正规方程可化为：ldl^tx＝a^tb。正规方程的求解可分为以下三个步骤：

1)解下三角方程组：lz＝a^tb

2)解对角方程组：dy＝z

3)解上三角方程组：l^tx＝y

其中z＝dl^tx，y＝l^tx。

经过以上步骤，即可求解出

二、匹配虚拟电阻吸收最大谐波功率

由电路原理可知，当负载电阻时，电阻负载吸收的有功功率到达最大值。因此，通过控制逆变器的输出谐波电流来等效一个虚拟电阻，可达到最大化吸收谐波有功功率的功能。

三、有源滤波器控制算法

图2是虚拟电阻型有源电力滤波器工作原理的示意图。图中，i^*c1为虚拟电阻参考电流，i^*c2为直流电压控制环输出的参考电流，i^*c为总参考电流，ic为实际输出电流；sa，sb，sc为开关管的驱动信号。从图中可以看出，虚拟电阻型有源电力滤波器与传统的谐波电流源型有源电力滤波器工作原理的主要区别在于，各自参考电流的生成方法不同。传统的有源电力滤波器的参考电流一般由直流侧的电压外环调节器输出值和负载的谐波、无功电流合成，而虚拟电阻型有源电力滤波器相当于一个谐波电阻，其参考电流应包含上述的直流电压外环调节输出值以及与公共耦合点各次谐波电压成比例的各次谐波电流，即图2中的虚拟电阻电流。因此，各次谐波电压与各次谐波电流的比值就是等效的各次谐波电阻。

图3是虚拟电阻电流的生成过程。图中，谐波检测环节用于提取公共耦合点的各次谐波电压，相位检测环节用于实时检测电压的相位信息，应用该相位信息把各次谐波电压转换为瞬时值，然后把各次谐波电压分别与各次谐波电导gh相乘，得到各次参考电流，通过对各次电流进行求和得到虚拟电阻电流参考值。图中所示的谐波电导gh为谐波电阻的倒数，通过调节各次谐波电导gh，就可等效地调节各次谐波电阻。而gh的计算是在逆变器输出谐波电流定期的小扰动下，用由谐波检测得到各次电压和电流谐波的实部、虚部求取最小二乘解所得。

图4是最小二乘法求虚拟电阻电流流程图，其步骤为：

(1)初始化，

(2)令k＝1

(3)在原参考谐波电流的基础上叠加不同幅值和相位小扰动

(4)用fft计算出逆变器的网点电压和入网电流的频谱

(5)令k＝k+1；如果k大于m则进入(6)，否则返回(3)；

(6)用uih(re)k、uih(im)k、iih(re)k、iih(im)k形成超定方程组ax＝b，并用对称矩阵的三角分解法求解正规方程组，得出系统的等效谐波阻抗zh＝rh+jxh；

(7)用fft计算电压频谱根据得到谐波参考电流频谱，再把频域参考电流转换为时域值；

(8)计算并存储功率最优值；

(9)计算实时功率phrt；

(10)计算abs(ph-phrt)，若大于阈值δ，则返回(2)，否则返回(7)。四、为了验证上文所述的最佳阻值自动调节方法和更具普遍性，本专利基于pscad/emtdc对图2进行仿真，其中，rs为0.00001ω，系统电感ls为0.00637h，即对于5次谐波相应的xs为10ω，等效负荷分别取r1为20ω，r2为10ω，假设三相平衡，取5次谐波为讨论对象，其中a相谐波电流为a相系统谐波电压为可以计算得到，当等效负载为r1时，理论最佳阻值为r^*f01＝8.945ω，相应电导值为g^*01＝0.112s；当等效负载为r2，理论最佳阻值为r^*f02＝7.0714ω，相应电导值为g^*02＝0.141s。

图5为等效负载从r2突变为r1时谐波功率变化曲线。在第一阶段0-1s，最小二乘法求得gh的值为0.1426，pfh稳定在1.141kw附近，在1s处，等效负载电阻突变后，最小二乘法求得gh的值为0.1187s，pfh最终稳定在1.62kw附近。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。