一种牵引变电所供电臂的谐波电流检测方法
【技术领域】
[0001] 本发明涉及一种牵引变电所供电臂的谐波电流检测方法，属于牵引供电系统的谐 波电流检测技术领域。
【背景技术】
[0002] 在电气化铁道牵引供电系统中，由于电力机车具有输出功率大、速度及负载状 况变化频繁的特点，导致牵引电网功率因数低、谐波含量高、负序电流大，这不但产生了 电能损耗，而且对公共电网及铁路沿线设备带来了严重的危害。作为一种有效地电能 质量控制装置，有源电能质量调节器（APQC)得到了广泛的关注，其工作原理为：通过对 牵引供电系统谐波电流的实时检测，控制逆变器产生一个大小相等、方向反相的谐波电 流并注入到牵引供电系统中，从而实现牵引供电系统谐波污染的抑制。APQC谐波抑制 性能的好坏主要依赖于谐波电流检测的实时性和精确性，因此谐波电流检测方法至关 重要。目前，研究人员提出了多种方法来检测谐波电流。在众多谐波电流检测方法中， Fryze-Buchholz-Dpenbrock(FBD)法由于其结构简单、实时性好，且不限于三相系统，得到 了广泛的应用。
[0003] FBD法根据负载电流与参考电压获取等效电导，从等效电导中分离出与基波有功 电流相关的线性等效电导，计算获取基波有功电流，从负载电流中消去基波有功电流即得 到广义谐波电流。
[0004] 传统的FBD法见参考文献"电气化铁路中谐波、无功、负序电流的实时检测方 法"（孙卓，姜新建，朱东起，电力系统自动化，2003, 27 (15) : 53-57.)。该方法将FBD法 应用于牵引供电系统中，并验证了其可行性。然而，传统的FBD法检测效果与低通滤波器 (LPF)的性能有很大关系。低通滤波器较小的截止频率能够获得很好地稳态检测精度，但其 动态响应速度较慢；较大的截止频率能够获得较快的动态响应速度，但其稳态检测精度就 会变差。因此，传统的FBD法在电流检测过程中存在动态响应速度和稳态检测精度之间的 固有矛盾。
【发明内容】
[0005] 本发明的发明目的就是提供一种牵引变电所供电臂的谐波电流检测方法，该方法 具有检测精度高、动态响应速度快的特点。
[0006] 本发明实现其发明目的所采用的技术方案是，一种牵引变电所供电臂的谐波电流 检测方法，其步骤如下：
[0007] A、信号采样
[0008] 对牵引变电所上、下行供电臂上的负载电流进行采样分别得到当前时刻n的上、 下行供电臂的负载电流离散值ia(n)、％(!〇;通过锁相环处理获得与牵引变电所上、下行 供电臂实际电压同相位的、当前时刻n的上、下行供电臂的参考电压离散值ua (n)、U{! (n); 其中下标a、0分别为上、下行供电臂的标号；
[0009] B、计算等效电导
[0010] 根据当前时刻n的上、下行供电臂的负载电流离散值ia (n)、％ (n)和参考电压离 散值ua (n)、up (n)，计算当前时刻n牵引变电所供电臂的等效电导Gp(n)，
[0011] Gp(n) =ia (n)ua (n)+ip (n)up (n)
[0012] C、获取线性等效电导
[0013] 将电流值为1的直流信号作为参考输入信号，通过最小均方算法的自适应滤波器 得到线性等效电导G(n)，G(n) =w(n) *1，其中，w(n)为自适应滤波器在当前时刻n的权系 数，其初始值为零；
[0014] D、计算基波有功电流
[0015] 根据步骤C1中获取的线性等效电导G(n)，计算当前时刻n的牵引变电所上、下行 供电臂上的基波有功电流ipla (n)，iplp (n)，
[0016]ipia(n) =G(n)ua (n)
[0017]ipip(n) =G(n)up (n)
[0018] E、获取谐波电流
[0019] 将步骤A中获取的牵引变电所上、下行供电臂上的负载电流ia(n)、i^n)分别 减去基波有功电流ipla (n)，iplf!(n)，得到牵引变电所上、下行供电臂上的广义谐波电流 ica (n)、icP (n)，
[0020] ica (n) =ia (n)-ipla (n)
[0021] icP (n) =ip (n)-ipl{! (n)
[0022] F、自适应滤波器权系数更新
[0023]FI、计算自适应滤波器的残差信号
[0024] 将步骤B中获取的当前时刻n牵引变电所供电臂的等效电导G>)减去C1步的线 性等效电导G(n)，得到自适应滤波器的当前时刻n的残差信号e(n)，e(n) =Gp(n)-G(n);
[0025] F2、计算自适应滤波器的步长因子
[0026] 对残差信号e(n)用牵引变电所供电臂等效电导Gp(n)的绝对值进行归一化处理， 得到当前时刻n的归一化残差信号s (n)，s (n)=e (n)/ |Gp (n)|;
[0027] 迭代计算出当前时刻n的归一化残差信号s(n)与前一时刻n-1的归一化残差信 号s (n-1)的自相关估计p (n)，
[0028] p (n) = 0 p (n-1) + (1-0 ) s (n) s (n-1)
[0029] 其中，P为自相关估计遗忘因子，其取值范围为0.800~0.999 ;p(n)的初始值为 〇 ;
[0030]根据自相关估计p (n)的平方p2 (n)，迭代计算出自适应滤波器在当前时刻n的步 长因子y(n)，
[0031] y (n) = a y (n-l) + yp2(n)
[0032] 其中，a为步长遗忘因子，其取值范围为0.800~0.999 ;y为步长调整参数，其 取值范围为〇? 001~〇? 050;y(n)的初始值为0;
[0033] F3、步长因子的修正
[0034]对自适应滤波器在当前时刻n的步长因子y (n)进行修正，得到修正后的步长因 子，，(/?)，
[0035]
[0036] 其中，ynax为步长因子y(n)的上限值，其取值范围为0. 1~1，y_为步长因子 y(n)的下限值，其取值范围为0? 001~0? 010;
[0037] F4、自适应滤波器权系数更新
[0038] 计算下一时刻n+1的自适应滤波器的权系数w(n+1)，
[0039] vr(/? +1} =win)+/Li{n)e{n)
[0040] G、令n=n+l，重复A至F步的操作，即可实时检测出牵引变电所上、下行供电臂 上的广义谐波电流。
[0041] 与现有技术相比，本发明的有益效果是：
[0042] 本发明采用自适应滤波器来替代传统FBD法原有的低通滤波器，通过采用变步长 的最小均方（LMS)算法控制的自适应滤波器，根据误差信号的自相关估计调整LMS算法步 长大小，在动态响应阶段，误差信号较大，根据变步长算法将产生较大的步长因子，加快了 检测方法的动态响应速度；在稳态响应阶段，误差信号较小，根据变步长算法将产生较小的 步长因子，提高了稳态检测精度。因此，本发明克服了传统的FBD法在电流检测过程中存在 的动态响应速度和稳态检测精度之间的固有矛盾，能够同时获得较快的动态响应速度和较 高的稳态检测精度，表现出了比传统的FBD法更好的检测效果。
[0043] 下面结合附图和【具体实施方式】对本发明进行详细说明。
【附图说明】
[0044] 图1是仿真实验中的负载电流波形。
[0045] 图2是传统FBD法在低通滤波器截止频率f= 20Hz时检测出的基波有功电流。
[0046] 图3是传统FBD法在低通滤波器截止频率f= 80Hz时检测出的基波有功电流。
[0047]图4是本发明检测出的基波有功电流。
【具体实施方式】
[0048] 本发明的一种【具体实施方式】是：一种牵引变电所供电臂的谐波电流检测方法，其 步骤如下：
[0049] A、信号采样
[0050] 对牵引变电所上、下行供电臂上的负载电流进行采样分别得到当前时刻n的上、 下行供电臂的负载电流离散值ia(n)、％(!〇;通过锁相环处理获得与牵引变电所上、下行 供电臂实际电压同相位的、当前时刻n的上、下行供电臂的参考电压离散值ua (n)、U{! (n); 其中下标a、0分别为上、下行供电臂的标号；
[0051] B、计算等效电导
[0052]根据当前时刻n的上、下行供电臂的负载电流离散值ia (n)、ip (n)和参考电压离 散值ua (n)、up (n)，计算当前时刻n牵引变电所供电臂的等效电导Gp(n)，
[0053]Gp(n) =ia (n)ua (n)+ip (n)up (n)
[0054] C、获取线性等效电导
[0055] 将电流值为1的直流信号作为参考输入信号，通过最小均方算法的自适应滤波器 得到线性等效电导G(n)，G(n) =w(n) *1，其中，w(n)为自适应滤波器在当前时刻n的权系 数，其初始值为零；
[0056] D、计算基波有功电流
[0057] 根据步骤C1中获取的线性等效电导G(n)，计算当前时刻n的牵引变电所上、下行 供电臂上的基波有功电流ipla (n)，iplp (n)，
[0058]ipia(n) =G(n)ua (n)
[0059]ipip(n) =G(n)u