微机继电保护的算法

单击此处编辑母版标题样式,*,单击此处编辑母版文本样式,*,*,第3章 微机继电保护的算法,第3章 微机继电保护的算法,第3章 微机继电保护的算法,3.1,概述,定义,根据模数转换器提供的输入电气量的采样数据进行分析、运算和判断，以实现各种继电保护功能的方法称为算法,分类,根据采样值计算出保护需要的量值，求电压、电流、再计算阻抗，然后和定值比较,直接模拟模拟型保护判据，判断故障是否在区内。,评价指标,精度和速度,3.2 正弦函数模型算法,假设被采样的电压、电流信号都是纯正弦量时，可以利用正弦函数的一系列特性，从若干个采样值中计算出电压、电流的幅值、相位以及功率和测量阻抗的量值,。,角频率,I,电流有效值,Ts,采样间隔,电流初相角,3.2.1 两点乘积算法,设,u,1,，u,2,是相差,90,o,的两个采样值，采样时刻分别为,n,1,，n,2,，,则:,直接计算线路电阻和电抗，将电压和电流写成复数形式,所以：,电抗和电阻及电压与电流间的夹角,：,3.2.2,三采样值乘积算法,三采样值积算法是利用三个连续的等时间间隔的采样值，通过适当的组合消去 项以求出采样的幅值和相位的方法。,为分析方便，设电压的初相位为0，电流滞后电压的角度为 ,则,在 时刻的采样值为,在 时刻的采样值为,对,以上采样值进行处理得：,更进一步的处理得：,当采样频率为600,Hz,，,时，则式（,3-17,）、（,3-18,）、（,3-19,）可简化为,：,或写成有效值,U,和,I，,有：,进一步推导：,三采样值简化算法与二采样值算法比较：,优点：三采样值算法只需等待60度的时间，而二采样值积算法则需等待90度的时间，所以三采样值积算法延时稍短一些，速度较快。,缺点：三采样值算法要用较多的乘除法。,3.2.3 导数算法（一次微分算法）,导数算法是利用输入正弦量在某一个时刻的采样值及在该时刻采样值的导数，即可算出有效值和相位的算法，如图,3-3,所示。,设正弦电压,u,、,电流,i,在,t,1,时刻的值为:,则,t,1,时刻的导数为:,由上述四式可以求出:,如何知道该点的导数值呢？,可利用相邻的采样数据近似计算在,t1,时刻电流、电压的导数，即：,可见利用导数法计算电流、电压的幅值、相位等电量时，只需要使用连续的三个采样值；算法的数据窗较短，仅为二个采样间隔，算式和乘积法相似也不复杂，其快速性较好。,3.2.4 二次微分算法,这种算法是为了消除导数算法受直流分量影响的缺点提出的。该算法在导数算法的基础上作了修正，采用一阶导数值和二阶导数值，代替导数算法的用采样值和一阶导数值的方法。,设正弦电压,u,、,电流,i,为:,则对应的一阶、二阶导数为,:,整理上式，可得:,3.2.5 半周积分算法,半周积分算法的依据是一个正弦量在任意半周期内绝对值的积分为一个常数,S,，,并且积分值,S,和积分的起始点初相角无关，如图,3-4,所示。,在上图中，画有断面线的两块面积显然是相等的。,即：,即正弦量半周期绝对值的积分比于幅值,Um，,从而半周期积分算法可用下式表示:,半周积分算法也有一定的滤波作用，因为在半波积分过程中，谐波中的正负半周相消，剩余未被抵消的部分占总和的比重就减少了。但由于它不能全部滤除谐波分量，因仍要求加入滤波环节。,3.3 傅里叶算法,傅立叶级数：,设,x(t),是一个周期,为,T,的时间函数（信号），则可以把它写成：,根据三角函数的正交性，可得基波分量的系数,写成复数形式,X,的有效值和相位,适于微机计算离散化需要，,a,1,b,1,的积分可以用梯形法,则求得,N,基波信号一周采样的点数，一共使用,N1,个采样值,X,k,第,k,点采样值,X,0,，,X,k,首末点采样值,对于基波工频，当,N12，,即30,o,一个采样点时,附注说明：,1.,X(t),是周期函数,，,求,a,1,，,b,1,可以使用任意一段,X(t)，,也就是该正弦函数取不同初相角。,2.,随着所取,X(t)“,段”的不同，相当于起点位置的不同、或者初相角的不同，,a,1,，,b,1,取得不同的值。换句话说，,a,1,，,b,1,是起点位置的函数。若设起点是,t,1,，,则,求得任意次谐波的振幅和相位的算法,：,以上算法在计算一个频率分量时，需要,N,次乘法和2（,N-1）,次加法，计算量相当大。因此，为了提高计算速度，常采用递推的傅氏算法。,递推的傅氏算法：,以上两式的运算只需要2次乘法和4次加减法，且与选取无关，极大地减少了运算量，因此具有广泛的使用价值。,3.4 递推最小二乘算法,最小二乘算法是广泛应用于数据处理和自动控制等领域中的一种经典方法，是误差理论内容。,其基本思想是：将输入待求量与一个已知的预设函数进行拟合，从而使得待求函数与预设函数尽可能逼近。其总方差或者最小均方差最小，从而近似的求出待求函数（如正弦函数的幅值和相角等）,。,对于已知函数的拟合,在线路保护中，采用最小二乘法对基频分量参数进行计算时，输入信号的预设模型可选择为：,用实、虚部表示有,：,而实际故障信号,y(t)，,可视为由预设信号,x(t),与附加的随机噪声信号,w,(t),共同组成，即:,可见：对于上式所示的输入信号而言，待确定的模型参数为：,非周期分量的初值,X,0,衰减时间常数,T,d,及基频分量的实部,X,R,和虚部,X,I,在实际应用中，为了简化计算，参数,T,d,通常可作为事先给定。,这样，待确定的参数可简化为只包括,X,0,、X,R,和,X,I,将式（,3-53,）用离散采样值形式表示时有：,递推型的最小二乘法,在微机保护中，采样数据是按采样频率逐个提供给计算机的，在得到新的采样数据后，若希望利用新的采样数据来改进原有的参数估计结果，以提高估计精度。,所以提出了递推型的最小二乘法。,递推最小二乘法算法的突出优点：,具有可变的数据窗,它的数据窗长度将随着采样值的增多而自动延长，参数的估计精度也随之逐步提高。,算法计算简便、收敛速度快、收敛过程稳定,对非周期分量和各种高频分量具有良好的滤波能力，在实际使用中，无需再附加其他的数字滤波器。,因此，递推最小二乘法算法可较好满足不同场合对参数计算精度和计算速度的不同要求。,3.5,解微分方程算法-,R-L,模型算法,基本原理,忽略线路分布电容，则输电线路等效为集中参数,R-L,模型。当短路发生时，有：,其中,R,L,是未知数，电压电流是可测量的,i R,L,u,在两个不同时刻分别测量,u、i,和,di,/,dt,，,就可以得到两个独立的方程，即,联立求解上述两个方程可求得两个未知数,R,和,L.,利用连续3次（,n-2,n-1,n）,的,采样值，并用采样值的差分代替微分，有：,由于 ，所以有：,3.5.1,差,分法,可解,得：,3.5.2,积分法,在两个不同的时间段内对（,3-57,）积分可得到两个独立的方程,：,在处理各项积分值时，可用梯形积分法来计算，得,可解得：,3.6,移相与滤序算法,已知电量 ，欲将其移相 角，则 的算法如下:,设第,n,个和第,n-k,个采样值分别为,:,则：,令：,代入上式，则有,：,式,中：,当取,k=1,时，,有 ，则:,式中，为将原,x(t),移相 角的电量。,3.6.2,滤序算法,在微机继电保护中经常使用序分量作保护的判据。,例如：,零序电流保护,发电机负序保护,距离保护的振荡闭锁功能,滤序器的基本表达式,：,1.,数据窗为,13.32,ms,和滤序器表达式,2.数据窗为6.66,ms,的滤序器表达式,3.数据窗为3.33,ms,的滤序器表达式,3.7 继电器特性算法,继电器特性算法是已知某种继电器的动作特性、不经电流、电压幅值和相位的中间计算环节，利用采样值直接得到要求的继电器特性方程的算法。,举例：,方向阻抗继电器算法,偏移特性阻抗继电器,相电流突变量算法,3.7.1 方向阻抗继电器算法,当方向阻抗继电器的整定阻抗为,Z,zd,，,按相位比较方法构成时的动作方程为：,3.7.2 偏移特性阻抗继电器,偏移特性阻抗继电器的正、反向整定阻抗分别为,Z,1,、Z,2,时，动作方程为:,与方向性阻抗继电器的算法推导方法相同，同样可写出偏移特性阻抗继电器的算法:,3.7.3 相电流突变量算法,微机保护中常采用相电流突变量作为启动元件,.,为消除因电网频率波动引起的不平衡电流，相电流突变量按下式计算,