微机保护实验报告

微机保护实验报告

The Standardization Office was revised on the afternoon of December 13, 2020

微机继电保护实验报告

项目名称：微机距离保护算法（1）姓名：陈发敏

学号：K03134163

班级：K0313416

实验时间：

实验地点：实验楼五楼

实验成绩：

一、 实验目的

1．熟悉MATLAB 桌面和命令窗口；

2．通过编写滤波程序、阻抗计算程序以及距离保护动作判据程序，了解微机保护工作原理。

3．定性分析各种算法的优缺点。

二、 实验内容

1、用“load ”函数导入短路电流数据和短路电压数据，对其进行滤波处理，要求滤除直流分量和二次谐波分量。注意观察数据的特征，数据第一列为时间，第二列为A 相值，第三列为B 相值，第三列为C 相值。观察滤波前后的波形。

2、编写微机保护算法程序，包括短路阻抗算法和动作判据算法（判据为相间距离保护判据），阻抗继电器的动作特性采用方向圆特性。并利用该程序对步骤1处理后的数据进行计算，观察保护的动作情况。距离保护的整定值为:Z set =+ 。

三、 实验模型及程序

1、 绘制滤波前后的电流、电压波形，并进行对比分析；

电流波形滤波前，短路瞬间电流幅值变大，到短路后的稳态呈曲线变化；经过滤波后，短路后的稳态比较平稳。

电压波形滤波前，短路瞬间电压幅值急剧变小；经过滤波后，短路后的稳态比较平稳，且短路后电压波形变化没有电流波形变化大。

4

4

4

4

4

4

5

5

5

5

5

5

2、 设计编写保护算法程序,绘制阻抗幅值变化的波形，并分析保护的动作情况。

由阻抗幅值变化的波形和保护的动作情况可知：左图的B 相的阻抗值太低，所以致使B 相动作有明显的变化。

附MATLAB 程序如下：

%实验3程序 clc; clear;

%电压电流数据导入

a=load('H:\To be completed\微机保护\jibao3_4\'); %导入电压量 b=load('H:\To be completed\微机保护\jibao3_4\'); %导入电流量 t=a(:,1)'; UA=a(:,2)'; UB=a(:,3)'; UC=a(:,4)'; IA=b(:,2)'; IB=b(:,3)'; IC=b(:,4)'; Ts=t(1,2)-t(1,1); N=Ts; m=size(t); %滤波处理 %%电流滤波 IIA=zeros(1,m(2)); IIB=zeros(1,m(2)); IIC=zeros(1,m(2)); for jj=101:m(2);

IIA(jj)=(IA(jj)-IA(jj-100))/2; IIB(jj)=(IB(jj)-IB(jj-100))/2; IIC(jj)=(IC(jj)-IC(jj-100))/2; end

subplot(3,1,1); plot(t,IIA,'r') title('电流滤波') subplot(3,1,2);

plot(t,IIB,'g')

subplot(3,1,3);

plot(t,IIC,'b')

figure

subplot(3,1,1);

plot(t,IA)

title('电流波形')

subplot(3,1,2);

plot(t,IB)

subplot(3,1,3);

plot(t,IC)

%%电压滤波

UUA=zeros(1,m(2));

UUB=zeros(1,m(2));

UUC=zeros(1,m(2));

for jj=101:m(2);

UUA(jj)=(UA(jj)-UA(jj-100))/2;

UUB(jj)=(UB(jj)-UB(jj-100))/2;

UUC(jj)=(UC(jj)-UC(jj-100))/2;

end

figure

subplot(3,1,1);

plot(t,UUA,'r')

title('电压滤波')

subplot(3,1,2);

plot(t,UUB,'g')

subplot(3,1,3);

plot(t,UUC,'b')

%利用两点乘积算法计算

%电压

USA=zeros(1,m(2));

USB=zeros(1,m(2));

USC=zeros(1,m(2));

for jj=N/4+1:m(2)

USA(jj)=sqrt((UUA(jj)*UUA(jj)+UUA(jj-N/4)*UUA(jj-N/4))/2); USB(jj)=sqrt((UUB(jj)*UUB(jj)+UUB(jj-N/4)*UUB(jj-N/4))/2); USC(jj)=sqrt((UUC(jj)*UUC(jj)+UUC(jj-N/4)*UUC(jj-N/4))/2); end

% %电流

for jj=N/4+1:m(2)

ISA(jj)=sqrt((IIA(jj)*IIA(jj)+IIA(jj-N/4)*IIA(jj-N/4))/2);

ISB(jj)=sqrt((IIB(jj)*IIB(jj)+IIB(jj-N/4)*IIB(jj-N/4))/2);

ISC(jj)=sqrt((IIC(jj)*IIC(jj)+IIC(jj-N/4)*IIC(jj-N/4))/2);

end