课程名称：微机继电保护
实验项目：电力系统继电保护仿真实验
实验地点：电力系统仿真实验室
专业班级：电气
学生姓名：指导教师：
本科实验报告
1200学号：0000000000 000000
000000
2015年 12 月 2 日
一、实验背景
微机继电保护指的是以数字式计算机(包括微型机)为基础而构成的继电保护。众所周知，传统的继电器是由硬件实现的，直接将模拟信号引入保护装置，实现幅值、相位、比率的判断，从而实现保护功能。而微机保护则是由硬件和软件共同实现，将模拟信号转换为数字信号，经过某种运算求出电流、电压的幅值、相位、比值等，并与整定值进行比较，以决定是否发出跳闸命令。
继电保护的种类很多，按保护对象分有元件保护、线路保护等；按保护原理分有差动保护、距离保护和电压、电流保护等。然而，不管哪一类保护的算法，其核心问题归根结底不外乎是算出可表征被保护对象运行特点的物理量，如电压、电流等的有效值和相位以及视在阻抗等，或者算出它们的序分量、或基波分量、或某次谐波分量的大小和相位等。有了这些基本电气量的计算值，就可以很容易地构成各种不同原理的保护。基本上可以说，只要找出任何能够区分正常与短路的特征量，微机保护就可以予以实现。
由此，微机保护算法就成为了电力系统微机保护研究的重点,微机保护不同功能的实现,主要依靠其软件算法来完成。微机保护的其中一个基本问题便是寻找适当的算法,对采集的电气量进行运算，得到跳闸信号，实现微机保护的功能。微机保护算法众多，但各种算法间存在着差异，对微机保护算法的综合性能进行分析,确定特定场合下如何合理的进行选择,并在此基础上对其进行补偿与改进,对进一步提高微机保护的选择性、速动性、灵敏性和可靠性,满足电网安全稳定运行的要求具有现实指导意义。
目前已提出的算法有很多种，本次实验将着重讨论基本电气量的算法，主要介绍突变量电流算法、半周期积分算法、傅里叶级数算法。
二、实验目的
1.了解目前电力系统微机保护的研究现状、发展前景以及一些电力系统微机保护装置。 2.具体分析几种典型的微机保护算法的基本原理。
3.针对线路保护的保护原理和保护配置，选择典型的电力系统模型，在MATLAB软件搭建仿真模型，对微机保护算法进行程序编写。
4.对仿真结果进行总结分析。
三、实验内容
1、采用MATLAB软件搭建电力系统仿真模型 2、采用MATLAB软件编写突变量电流算法 3、采用MATLAB软件编写半周积分算法 4、采用MATLAB软件编写傅里叶级数算法算法
四、实验步骤
1.突变量电流算法、半周积分算法、傅里叶级数算法简介 1.1突变量电流算法
继电保护装置的启动元件用于反应电力系统中的扰动或故障。微机保护装置中的启动元件是由软件来实现的。它的工作原理目前一般采用反映两相电流差的突变量，其公式为
?Iab?iabn?iab(n?N)?iab(n?N)?iab(n?2N)?Ibc?ibcn?ibc(n?N)?ibc(n?N)?ibc(n?2N)（1） ?Ica?ican?ica(n?N)?ica(n?N)?ica(n?2N)其中
iabn?ian?ibn
公式中
N—一个工频周期内的采样点数
ibcn?ibn?icn（2） ican?icn?ianian、ibn、icn—当前时刻的采样值
iab(n?N)、ibc(n?N)、ica(n?N)—一周前对应时刻的采样值 iab(n?2N)、ibc(n?2N)、ica(n?2N)—两周前对应时刻的采样值
以?Iab为例，正常运行时ian、ia(n?N)、ia(n?2N)的值近似相等，所以?Iab?0，启动元件不动作，如图1所示。
图1 系统正常运行时采样值比较
电力系统正常运行但频率发生变化偏离50Hz时，则ian、ia(n?N)、ia(n?2N)的值将不相等，。
这是因为采样时按时间间隔进行的，频率变化时，ian和ia(n?N)两采样值将不是相差一个周期的采样值，于是ian-ia(n?N)、ia(n?N)-ia(n?2N)将出现差值，且差值接近相等。此时?Iab仍然为零或很小。
系统发生故障时，由于故障电流增大，于是ian将增大，ia(n?N)为故障前电流，故ian-ia(n?N)反映出由于故障电流产生的突变量电流，ia(n?N)-ia(n?2N)仍接近为零，从而?Iab反映了故障电流的突变量，如图2所示。
图2 故障后电流的突变
1.2半周积分法
半周积分算法的依据是一个正弦量在任意半个周期内绝对值的积分为一个常数S，即
T2S=?0T22Isin(?t+?)dt（3）
??02Isin(?t)dt?22?I
积分值S与积分起点的初相角?无关，因为画有断面线的两块面积显然是相等的，如图3所示。式（3）的积分可以用梯形法则近似求出：
N??211S??i0??ik?iN?Ts（4）
?222?k?1????式中
ik：第k次采样值
N：每个周期的采样点数
i0：k?0时的采样值
iN：k?N时的采样值
22Ts：采样间隔
图4所示，只要采样率足够高，用梯形法则近似积分的误差可以做到很小。
S?????S??N20k
图3半周期积分法原理示意图图4用梯形法近似半周期积分示意图
求出S值后，应用式（3）即可求得有效值
I=S?1.3傅里叶级数算法
?22
傅里叶级数算法（简称傅氏算法）的基本思路来自傅里叶级数，算法本身具有滤波作用。它假定被采样的模拟信号是一个周期性的时间函数，除基波外还含有不衰减的直流分量和各次谐波，可表示为
x(t)??Xnsin(n?1t??n)??[(Xnsin?n)cosn?1t?(Xnsin?n)sinn?1t]n?0n?0????[bncosn?1t?ansinn?1t]n?0?(n?0,1,2....)（5）
bn分别为直流、式中an、基波和各次谐波的正弦项和余弦相得振幅，其中bn?Xnsin?n、
an?Xncos?n。
由于各次谐波的相位可能是任意值的，所以，把它们分解成有任意振幅的正弦项和余弦项之和。a1、b1分别为基波分量的正、余弦项的振幅，b0为直流分量的值。
根据傅氏级数的原理，可以求出a1、b1分别为
2a1??x(t)sin(?1t)dt（6）
T0T