改进距离继电器中功率摆动阻塞和解阻塞功能的新方法（Matlab代码实现）

👨‍🎓个人主页：研学社的博客  

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章讲解

---

💥1 概述

电力系统中的功率摆动条件使距离继电器在某些情况下采取错误的断路器动作。在这种情况下，功率摆幅阻塞和功率摆幅去阻塞是距离继电器的两个重要功能，用于区分功率摆幅和故障，从而保持系统的稳定性和电源的连续性。该文提出一种基于泰勒级数展开的样本估计和误差计算方法，以提高继电器的功率摆幅检测性能。该方法的性能与文献中的四种常规方法和一种新方法进行了比较。分析在两个系统中进行：具有平行线的单机无限总线系统和68总线新英格兰测试-纽约电力系统，用于故障，功率摆动和功率摆动情况下的故障。仿真和性能分析分别使用EMTDC/PSCAD和MATLAB进行。结果表明，当考虑正确操作的性能、算法响应时间、输出清晰度和处理器程序执行所经过的时间时，所提方法总体上具有最佳性能。研究发现，采用所提方法提高了正确运行下的距离继电器性能和采取功率摆动阻塞和去阻塞动作的速度。

文献来源：

📚2 运行结果

 

部分代码：

load('current_fault');time1=time;current1=current;
load('current_swing');time2=time;current2=current;
load('current_fault_during_swing');time3=time;current3=current;

N=20;%number of samples in one cycle

for p=1:length(time1)
    DI1(p)=0;
    if p>=2*N-1
        for q=0:N-1
            DI1(p)=DI1(p)+(current1(-q+p-3)-3*current1(-q+p-2)+3*current1(-q+p-1)-current1(-q+p)+mean(current1(-q+p-N+1:-q+p)))^2;           
        end            
    end
end

%fprintf('Minimum (DI) in fault case = %g\n',min(DI1));
fprintf('Maximum (DI) in fault case = %g\n',max(DI1));

for p=1:length(time2)
    DI2(p)=0;
    if p>=2*N-1
        for q=0:N-1
            DI2(p)=DI2(p)+(current2(-q+p-3)-3*current2(-q+p-2)+3*current2(-q+p-1)-current2(-q+p)+mean(current2(-q+p-N+1:-q+p)))^2;           
        end            
    end
end

%fprintf('Minimum (DI) in power swing case = %g\n',min(DI2));
fprintf('Maximum (DI) in power swing case = %g\n',max(DI2));

for p=1:length(time3)
    DI3(p)=0;
    if p>=2*N-1
        for q=0:N-1
            DI3(p)=DI3(p)+(current3(-q+p-3)-3*current3(-q+p-2)+3*current3(-q+p-1)-current3(-q+p)+mean(current3(-q+p-N+1:-q+p)))^2;           
        end            
    end
end

%fprintf('Minimum (DI) in fault during power swing case = %g\n',min(DI3));
fprintf('Maximum (DI) in fault during power swing case = %g',max(DI3));

subplot(1,3,1);plot(time1,DI1);xlim([0.6,0.75]);ylim([0,110]);grid on;xlabel('Time (s)');ylabel('DI');
title('Fault');
subplot(1,3,2);plot(time2,DI2);xlim([2,3]);ylim([0,110]);grid on;xlabel('Time (s)');ylabel('DI');
title('Power Swing');
subplot(1,3,3);plot(time3,DI3);xlim([2,2.3]);ylim([0,110]);grid on;xlabel('Time (s)');ylabel('DI');
title('Fault during Power Swing');
saveas(gcf, '../results/testSMIB.png');

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]I. G. Tekdemir and B. Alboyaci, "A novel approach for improvement of power swing blocking and deblocking functions in distance relays," 2017 IEEE Power & Energy Society General Meeting, Chicago, IL, USA, 2017, pp. 1-1, doi: 10.1109/PESGM.2017.8273926.

🌈4 Matlab代码、数据、文章讲解