微机继电保护实验报告

1、.本科实验报告课程名称： 微机继电保护 实验项目： 电力系统继电保护仿真实验 实验地点： 电力系统仿真实验室 专业班级： 电气1200 学号：0000000000 学生姓名： 000000 指导教师： 000000 2015年 12 月 2 日一、实验背景微机继电保护指的是以数字式计算机(包括微型机)为基础而构成的继电保护。众所周知，传统的继电器是由硬件实现的，直接将模拟信号引入保护装置，实现幅值、相位、比率的判断，从而实现保护功能。而微机保护则是由硬件和软件共同实现，将模拟信号转换为数字信号，经过某种运算求出电流、电压的幅值、相位、比值等，并与整定值进行比较，以决定是否发出跳闸命令。继电保护

2、的种类很多，按保护对象分有元件保护、线路保护等；按保护原理分有差动保护、距离保护和电压、电流保护等。然而，不管哪一类保护的算法，其核心问题归根结底不外乎是算出可表征被保护对象运行特点的物理量，如电压、电流等的有效值和相位以及视在阻抗等，或者算出它们的序分量、或基波分量、或某次谐波分量的大小和相位等。有了这些基本电气量的计算值，就可以很容易地构成各种不同原理的保护。基本上可以说，只要找出任何能够区分正常与短路的特征量，微机保护就可以予以实现。由此，微机保护算法就成为了电力系统微机保护研究的重点,微机保护不同功能的实现,主要依靠其软件算法来完成。微机保护的其中一个基本问题便是寻找适当的算法,对采集

3、的电气量进行运算，得到跳闸信号，实现微机保护的功能。微机保护算法众多，但各种算法间存在着差异，对微机保护算法的综合性能进行分析,确定特定场合下如何合理的进行选择,并在此基础上对其进行补偿与改进,对进一步提高微机保护的选择性、速动性、灵敏性和可靠性,满足电网安全稳定运行的要求具有现实指导意义。目前已提出的算法有很多种，本次实验将着重讨论基本电气量的算法，主要介绍突变量电流算法、半周期积分算法、傅里叶级数算法。二、实验目的1. 了解目前电力系统微机保护的研究现状、发展前景以及一些电力系统微机保护装置。2. 具体分析几种典型的微机保护算法的基本原理。3. 针对线路保护的保护原理和保护配置，选择典型的

4、电力系统模型，在MATLAB软件搭建仿真模型，对微机保护算法进行程序编写。4. 对仿真结果进行总结分析。三、实验内容1、采用MATLAB软件搭建电力系统仿真模型2、采用MATLAB软件编写突变量电流算法3、采用MATLAB软件编写半周积分算法4、采用MATLAB软件编写傅里叶级数算法算法四、实验步骤1.突变量电流算法、半周积分算法、傅里叶级数算法简介1.1突变量电流算法继电保护装置的启动元件用于反应电力系统中的扰动或故障。微机保护装置中的启动元件是由软件来实现的。它的工作原理目前一般采用反映两相电流差的突变量，其公式为 （1）其中 （2）公式中N一个工频周期内的采样点数、当前时刻的采样值、一周

5、前对应时刻的采样值、两周前对应时刻的采样值以为例，正常运行时、的值近似相等，所以，启动元件不动作，如图1所示。 图1 系统正常运行时采样值比较电力系统正常运行但频率发生变化偏离50Hz时，则、的值将不相等，。这是因为采样时按时间间隔进行的，频率变化时，和两采样值将不是相差一个周期的采样值，于是-、-将出现差值，且差值接近相等。此时仍然为零或很小。系统发生故障时，由于故障电流增大，于是将增大，为故障前电流，故-反映出由于故障电流产生的突变量电流，-仍接近为零，从而反映了故障电流的突变量，如图2所示。图2 故障后电流的突变1.2半周积分法半周积分算法的依据是一个正弦量在任意半个周期内绝对值的积分为

6、一个常数S，即 （3）积分值S与积分起点的初相角无关，因为画有断面线的两块面积显然是相等的，如图3所示。式（3）的积分可以用梯形法则近似求出： （4）式中：第次采样值：每个周期的采样点数：时的采样值：时的采样值：采样间隔图4所示，只要采样率足够高，用梯形法则近似积分的误差可以做到很小。图3 半周期积分法原理示意图 图4 用梯形法近似半周期积分示意图求出S值后，应用式（3）即可求得有效值1.3傅里叶级数算法傅里叶级数算法（简称傅氏算法）的基本思路来自傅里叶级数，算法本身具有滤波作用。它假定被采样的模拟信号是一个周期性的时间函数，除基波外还含有不衰减的直流分量和各次谐波，可表示为 （5）式中、分别

7、为直流、基波和各次谐波的正弦项和余弦相得振幅，其中、。由于各次谐波的相位可能是任意值的，所以，把它们分解成有任意振幅的正弦项和余弦项之和。、分别为基波分量的正、余弦项的振幅，为直流分量的值。根据傅氏级数的原理，可以求出、分别为 （6） （7）由积分过程可以知道，基波分量正、余弦项的振幅、已经消除了直流分量和整次谐波分量的影响。于是中的基波分量为 （8）合并正弦、余弦项，可写为 （9）式中：基波分量的有效值：时基波分量的相角将用角公式展开可得： （10） （11）用复数表示为 （12）因此，可根据、，求出有效值和相角为 （13） （14）用微机处理时，式（13）和式（14）的积分可以用梯形法则求

8、得： （15） （16）式中：基波信号的一周期采样点数；：第次采样点数；、：分别为和时的采样值。2.使用MATLAB软件搭建电力系统仿真模型2.1MATLAB/Simulink仿真软件在电力系统中的应用介绍MATLAB/Simulink软件是由美国Math Works 公司开发的著名的动态仿真系统，它是MATLAB的一个附加组件，为用户提供了一个建模与仿真的工作平台。它能够实现动态系统建模与仿真的模块集成，而且可以根据设计和使用的要求对系统进行优化，提高建模与仿真的效率。MATLAB/Simulink软件提供了多个学科的仿真系统工具箱，和一些常用工具箱模块，用户可以根据需要方便地选用合适的工具

9、箱进行系统的建模与仿真分析。对于电力系统而言就有专门的工具箱模块库SimPowerSystem供用户使用，其功能强大，包含的电气元件种类多，处理函数模块丰富，为电力系统的仿真与研究提供了很大的便利，是电气工程专业必不可少的研发工具。在SimPowerSystem模块库中，包括10类模块库，即电源元件库（Electrical Sources）、线路元件库（Element）、电力电子元件库（Power Electronics）、电机元件库（Machines）、连接器元件库（Connctors）、电路测量模块元件库（Measurements）、附加元件库（Extras）、演示教程（Demos）、电力

10、图形用户分析界面（Powergui）、电力系统元件库（Powerlip-modles）。正是这些丰富的模块库使得电力系统的仿真变得方便、快捷，并且科学精确。本实验选用MATLAB/Simulink软件作为电力系统微机保护仿真的平台。2.2使用MATLAB软件搭建电力系统故障暂态仿真模型图5给出了一个实际简单电力系统模型，以此来建立电力系统暂态仿真模型。G图5 系统模型打开MATLAB后，点击Simulink库，进入Simulink库界面后，点击“File”弹出菜单栏，选择“New”“Model”，打开建立新模型窗口，将建立电力系统暂态仿真模型所需的仿真模块由模块库中拷贝后粘贴到新模型窗口内，或

11、按住所需模块拖至新模型窗口，再按图5所示的系统模型连接关系进行连接。在建立电力系统仿真模型时主要用到以下仿真模块如表1所示。MATLAB2014Ra中SimPowerSystems库在Simscape库目录下。表1 搭建系统所需模块模块名路径三相等值系统模块SimPowerSystems /Specialized Technology /Electrical Sources三相双绕组变压器模块SimPowerSystems /Specialized Technology /Elements / Three-Phase Transformer(Two Windings)三相分布参数线路模块Sim

12、PowerSystems /Specialized Technology /Elements /Distributed Parameter Line三相串联负载模块SimPowerSystems /Specialized Technology /Elements / Three-Phase Series RLC Load三相故障模块SimPowerSystems /Specialized Technology /Elements /3-Phase Fault三相电压电流测量模块SimPowerSystems /Specialized Technology /Measurements /Thre

13、e-Phase V-I Measurement电力图形用户分析界面SimPowerSystems /Specialized Technology /Powergui示波器Simulink/Sinks/ Scope多路分配器模块Simulink/ Commonly Used Blocks/ Demux将以上模块由模块库浏览器中拖放到“建立新模型”窗口中。新模型的缺省名为untitled，可保存为其他文件名。本例命名为shortcircuit.mdl。仿真模型图如下图所示。图6 仿真模型做出三相接地短路、BC两相短路三相接地短路IaIbIcV-IBC相短路IaIbIcV-I3、使用MATLAB软件

14、编写算法程序3.1准备仿真数据本节利用上一节建立的电力系统仿真模型在Simulink环境中进行电力系统故障暂态仿真所产生的数据对微机保护算法进行综合仿真，用以分析微机保护保护算法的性能。算法综合仿真的内容包括：（1）准备仿真数据。由于设置了示波器参数，使用的数据可以在MATLAB电力系统仿真运行后自动存入MATLAB主界面的工作区，可以命令行窗口中调用数据。（注意：此时数据为模型系统文件运行结束后示波器的数据。）（2）微机保护保护算法编程。（3）调用电力系统故障暂态数据，对故障数据进行保护算法计算，并以图形方式显示保护算法计算的结果。3.2突变量电流法突变量电流法算法程序：%突变量算法程序Ts

15、=0.0;Te=0.2;N=200; %仿真起始(start)、结束(end)时间，一周期采样点数f=50; %仿真频率T=1/f; %一个周期时长ts=T/N; %一个周期内采样间隔ia=Ia.signals.values; %A相电流采样值ib=Ib.signals.values; %B相电流采样值iab=ia-ib; %计算得iabn=1:1:2001; %n可从1连续取至2001t=(n-1)*ts; %使计算值与时间相对应Delta_ia=; %定义突变量ia为一个数组Delta_iab_n=; %定义突变量iab_n为一个数组Delta_Iab=; %定义突变量Iab为一个数组fo

16、r n=401:1:2001 %n从401连续取至2001Delta_ia(n)=abs(ia(n)-ia(n-N); %计算得数组ia内的值Delta_iab_n(n)=abs(iab(n)-iab(n-N); %计算得数组iab_n内的值Delta_Iab(n)=abs(Delta_iab_n(n)- Delta_iab_n(n-N); %计算得数组Iab内的值endplot(t, Delta_ia,b-,t, Delta_iab_n,g-,t, Delta_Iab,r-.) ; legend(Delta_ia, Delta_iab_n, Delta_Iab) %画图在运行三相接地短路故障

17、仿真结束后，在MATLAB中运行该程序，得到图形如下。图38 突变量算法得到的、曲线程序从0.04s开始计算，0.04s0.06s内计算值为0，0.06s故障发生后，突变量急剧变化，几乎为故障电流的两倍，足以触动继电保护的启动元件。请同学们自行做出两相相间短路、单相接地短路的突变量变化图。AB相间短路BC相间短路A相短路B相短路C相短路三相短路在下表中填入各短路类型下0.07s时的、值。表2 各短路类型下0.07s时的、值计算值短路类型 AB相短路83.51191.9249.8BC相短路85.69161.3219.1A相短路85.88196.8254.6B相短路85.88196.8254.6三

18、相短路398.3191.9249.83.3半周积分法m=floor(n-st)/16)-1; 一周期采样点数为64，每隔16个点计算的半周积分法算法程序：%半周积分法Ts=0.0;Te=0.2;N=128; %仿真起始(start)、结束(end)时间，一个周期内采样点数f=50; %仿真频率T=1/f; %一个周期时长ts=T/N; %仿真步长、采样间隔ia=Ia.signals.values; %A相电流采样值n=floor(N*0.2/T+1); %计算得0.2s时刻所对应的Ia序列值个数st=(0.05-T/2)/T*N+1; %起始计算点m=(n-st-N/2)/16+1; %每隔1

19、6个点向前计算得半周积分，从0.05s开始计算至0.2s需要计算的次数A=; %定义半周积分值A为一个数组for j=0:1:(m-1); %使j从0取至（m-1）for i=(st+j*16):1:(st+j*16+N/2) %使i取半个周期的采样点数N/2 S=0; %使S初始值为0 S=S+abs(Ia.signals.values(i)*ts; %用矩形法计算半个周期内积分值 i=i+1; %使i加1直至循环结束 end A(j+1)=S; %循环所得为半周积分值A中第j+1个元素 j=j+1; %使j加1直至循环结束endplot(st+N/2-1):16:(n-1)*ts,A,*-

20、) %画图输入程序，按下Enter键后，得到图形如图39。在MATLAB命令行窗口输入A，按下Enter键后，得到半周积分法计算的积分值。图39 半周积分算法得到的积分值曲线程序从0.05s开始计算，每次向前计算半个周期（即0.01s），可见0.05s0.06s计算得半周积分值波动，0.06s故障发生后，半周积分量向上攀升，在0.06s0.12s内震荡剧烈，可以触动继电保护的启动元件，0.12s 故障消失后，半周积分值骤然减小，在0.12s0.2s内震荡减小。在电力图形用户分析界面（Powergui）模块的离散系统仿真模型中通过改变采样时间改变采样周期，获得不同采样周期下的A值，按要求填入下表

21、。并按表2要求做出两相相间短路、单相接地短路以及三相短路的半周积分量变化图。表2 不同采样周期下短路前后半周积分法计算的A值一周期采样点数64（每隔8个点计算半周积分）128（每隔16个点计算半周积分）256（每隔32个点计算半周积分）短路发生时刻0.06s前（最小值）0.06s后（最大值）0.06s前（最小值）0.06s后（最大值）0.06s前（最小值）0.06s后（最大值）面积AB相短路0.022360.14070.0014250.082440.00072330.04166A相短路0.022360.026270.0014250.045920.00072330.02297三相短路0.0223

22、60.17730.0014250.088370.00072330.04401AB相间短路64点128点256点A相接地短路64点128点256点三相短路64点128点256点3.4傅里叶算法一周期采样点数为200的傅里叶算法程序：t=Ia.time; N=200;f=50; %仿真频率T=1/f; %一个周期时长ts=T/N; %一个周期内采样间隔i=1:N;%输入信号为Iasubplot(221); %画2行2列图，此图为第1个plot(t,Ia.signals.values);legend(原始数据);xlabel(t/s);ylabel(A);%画出输入信号为Ia%计算基波电流幅值hs(

23、i)=sin(2*pi*i/N); %傅里叶滤波系数hc(i)=cos(2*pi*i/N);ys=filter(hs,1,Ia.signals.values); %正弦幅值yc=filter(hc,1,Ia.signals.values); %余弦幅值ym=2*sqrt(ys.2+yc.2)/N;subplot(222); %在2行2列图中此图为第2个plot(t,ym);legend(基波幅值);xlabel(t/s);ylabel(A); %画出输入信号Ia基波幅值图%计算3次谐波电流幅值hs3(i)=sin(3*2*pi*i/N);%傅里叶滤波系数hc3(i)=cos(3*2*pi*i/N);ys3=filter(hs3,1,Ia.signals.values);%正弦幅值yc3=filter(hc3,1,Ia.signals.values);%余弦幅值ym3=2*sqrt(ys3.2+yc3.2)/N;subplot(223); %在2行2列图中此图为第3个plot(t,ym3);legend(3次谐波幅值);xlabel(t/s);ylabel(A); %画出输入信号Ia三次谐波幅值图%计算5次谐波电流幅值hs5(i)=sin(5*2*pi*i/N);%傅里叶滤波系数hc5(i)=cos(5*2*pi*i/N);ys5=filt