微机距离保护阻抗算法

1、微机距离保护的阻抗算法和特性分析摘要分析了常用微机距离保护所采用的阻抗算法原理和动作特性, 以实现距离保护可靠切除区内相间故障和单相接地故障, 而区外故障不误动的功能。关键词微分方程算法阻抗特性前言随着大规模集成电路技术的飞速发展, 微型计算机保护已得到了普遍的应用。在电力系统常规保护中, 距离保护遇到的问题最多, 因此, 在计算机保护的发展过程中, 计算机距离保护吸引了很多人的注意。计算机继电保护是用数学运算方法实现故障量的测量、分析和判断, 而运算的基础是若干个离散的、量化了的数字采样序列i k, uk , 因此微机保护的一个基本问题是寻找适当的离散运算方法, 使运算结果的精度能满足工程要

2、求, 而计算耗时又尽可能短。近10 多a 来, 国内外的继电保护工作者作了大量的研究, 提出了许多适合于计算机保护的计算方法, 如导数算法、采样积分算法、傅氏算法和微分方程算法等。1微机距离保护的算法在现行南京电力自动化设备总厂生产的11, 15型以及四方公司生产的CSL100 系列微机线路距离保护大多采用微分方程算法。它是假设输电线路由电阻和电感组成, 不同故障情况下建立的微分方程如下:1. 1相间短路时此时, 短路点的电压为零, 则有:u = iR + Ldi / dt 或u = L ( R i/ L + di / dt)写成离散形式为:uk = L ( Ri k/ L + ( ik +

3、1 - ik- 1 ) / 2T s)因对输电线路, R / L = 􀀁为常数, 故得L = uk / ( 􀀁ik + ( ik+ 1 - i k- 1) / 2T s)R = ( uk - L ( ik + 1 - ik- 1 ) / 2T s) / i k或R = uk / ( ik + 1/ 􀀁( i k+ 1 - ik- 1) / 2Ts)根据X = 􀀂L 即可算出电抗值。事实上, 电感L与短路距离成正比, 用电感值作距离量, 还可以不受系统频率变化的影响。1. 2短路点经过渡电阻短路时电力系统中短路点实际上经常是

4、有过渡电阻的, 为了克服短路点的过渡电阻给阻抗继电器的测量带来误差, 常用单相接地时的微分方程:u = L d( i + K L3I 0 ) / dt + R( i + K r 3I 0) + uf式中K L = ( L 0 - L 1) / 3L 1K r = ( R 0 - R1 ) / 3R 1uf 为短路点电压写成离散形式时为:Uk = L ( ( ik+ 1 - ik - 1 + 3K L( I 0k+ 1 -I 0k - 1) ) / 2T s + 􀀁( ik + K r 3I 0k ) ) + ufk ( 1)令Dk = ( ik+ 1 - ik - 1 + 3

5、K L( I 0k + 1 - I 0k- 1 ) ) / 2Ts+ 􀀁( ik + K r 3I 0k)Dk 式中各量均为测量值及常数。故DK 为可计算出的系数。计算L 值需要知道Ufk, Ufk 是短路电压, 无法测得。因相对来说, 零序网络是变化不大的, 此时如假定网络结构已知, 则存在下面的关系:uf = 3I 0fR f = 3I 0fR f/ k f0式中Rf 为短路点过渡电阻; kf0 = I 0 / I 0f 为零序网络的零序电流分配系数如果假定短路点两侧零序网络阻抗角相同, 则k f0 为实常数。3I 0 为流过继电器的零序电流, 是可测量的量。此外, 如再

6、假定在2 3 个采样时间间隔内过渡电阻R f 值保持不变, 则在2 个采样时刻根据( 1)式, 可写出下列方程组Uk = L Dk + I 0k3R f/ kf0 ( 2)Uk+ 1 = L Dk + 1 + I 0k+ 13R f/ kf0 ( 3)联解上述方程组可得:L = ( Uk I 0k+ 1 - Uk+ 1I 0k ) / ( Dk I 0k+ 1 - Uk+ 1I 0k )本算法是在上述假定条件下实现的, 因此计算结果存在一定的误差。当采用较完善的滤波方法时,可变为正弦模型下的微分方程算法, 仍可保持良好的克服过渡电阻的优点, 保证计算精度。2. 1多边形方向阻抗特性多边形方向阻

7、抗特性如图1。角度取值:a. 为防止在保护区末端经过渡电阻短路时可能出现的超范围动作, 一般􀀁可取710°。b. 考虑到经过过渡电阻短路时, 由过渡电阻引起的附加测量阻抗, 始端故障时比末端故障时小, 所以􀀁1 < 90°,通常取60°。c. 为保证出口经过渡电阻短路时能可靠动作,􀀁2 通常取15°。d. 为保证被保护线路发生金属性短路故障时能可靠动作, 􀀁3 同样可取15°。其动作判据为:A : Rm X DZctan( 90°+ 􀀁3

8、)B: X m R DZtan􀀁2C: X m X DZ - Rmtan􀀁D: R m R DZ + X mctan􀀁1整个阻抗元件的动作逻辑方程为: Z = A BCD图1 只有2 个参数X 和R 可以整定, 􀀁, 􀀁1 , 􀀁2, 􀀁2. 2四边形方向阻抗特性四边形方向阻抗特性如图2。其动作判据为:R mtan􀀁2 X m X DZX mctan( 90°+ 􀀁3 ) R m R DZ + X mctan𙦙

9、7;1􀀁1 , 􀀁2, 􀀁3 都是预先整定的参数。因此, ctan􀀁1 ,tan􀀁2, ctan􀀁3 都是常数。2. 3阻抗特性的偏移当采用四边形或多边形阻抗元件时, 基本能保证可靠切除区内相间故障和单相接地故障。为了避免PT 在线路侧而故障为出口三相短路时, 距离保护拒动, 阻抗动作特性在原四边形或多边形特性的基础上加上一个包括座标原点的小矩形特性。并采用记忆特性来计算短路阻抗值。一般情况下, 实现偏移特性的小矩形的X , R 取值如表1 所示。表1X, R 取值X 取值当X DZ &

10、lt; 1 时, 取X DZn/ 2n 指距离n 段当X DZ 1 时, 取保0. 5R 取值取R DZ / 4 与X 偏移量之小者 三 、数字滤波数字滤波器不同于模拟滤波器，它不是一种纯硬件构成的滤波器，而是由软件编程去实现，改变算法或某些系数即可改变滤波性能，即滤波器的幅频特性和相频特性。 以差分滤波为例做简单介绍。差分滤波器输出信号的差分方程形式为 (81)式中，x(n)、y(n)分别是滤波器在采样时刻n(或n)的输入与输出；x(n-k)是n时刻以前第k个采样时刻的输入，k1。对式(8-1)进行变换，可得传递函数H(z) （82）将 代入式(8-2)中，即得差分滤波器的幅频特性

11、和相频特性分别为式(8-3)及式(8-4)  （83） （84）由式(8-3)可知，设需滤除谐波次数为m，差分步长为k(k次采样)，则此时=m1=m·21，应使 =0。令则有 ； （85）当N（即s和1）取值已定时，采用不同的l和k值，便可滤除m次谐波。四、正弦函数模型算法1半周积分算法半周积分算法的依据是 （86）即正弦函数半周积分与其幅值成正比。式(8-6)的积分可以用梯形法则近似求出： （87）式中第K次采样值; N一周期T内的采样点数; k0时的采样值; kN/2时的采样值。求出积分值S后，应用式(8-6)可求得幅值。 2导数算法 导数算法是利用正弦函数的导数为余弦

12、函数这一特点求出采样值的幅值和相位的一种算法。设 则 （88）很容易得出 （89） （810）和 （811）根据式(8-8)，我们也可推导出 （812） （813）式(8-9)式(8-13)中，u、i对应tk 时为uk 、ik，均为已知数，而对应tk-1和tk+1的u、i为uk-1、uk+1、ik-1、ik+1，也为已知数，此时 （814） （815） （816） （817）导数算法最大的优点是它的“数据窗”即算法所需要的相邻采样数据是三个，即计算速度快。导数算法的缺点是当采样频率较低时，计算误差较大。 五两采样值积算法两采样值积算法是利用2个采样值以推算出正弦曲线波形，即用采样值的乘积来计算

13、电流、电压、阻抗的幅值和相角等电气参数的方法，属于正弦曲线拟合法。这种算法的特点是计算的判定时间较短。设有正弦电压、电流波形在任意二个连续采样时刻tk、tk+1（tk）进行采样，并设被采样电流滞后电压的相位角为，则tk和tk1时刻的采样值分别表示为式(8-18)和式(8-19)。 （818） （819）式中，TS为两采样值的时间间隔，即TStk+1tk 。 由式(8-18)和式(8-19)，取两采样值乘积，则有 （820） （821） （822） （823）式(8-20)和式(8-21)相加，得 （824）式(8-22)和(8-23)相加，得 （825）将式(8-25)乘以cosTS再与式(8

14、-24)相减，可消去tk项，得 （826）同理，由式(8-22)与式(8-23)相减消去tk项，得 （827）在式(8-26)中，如用同一电压的采样值相乘，或用同一电流的采样值相乘，则q 0°，此时可得 （828） （829）由于TS、sinTS、cosTS均为常数，只要送入时间间隔TS的两次采样值，便可按式(8-28)和式(8-29)计算出Um、Im 。 以式(8-29)去除式(8-26)和式(8-27)还可得测量阻抗中的电阻和电抗分量，即 （830） （831） 由式(8-28)和式(8-29)也可求出阻抗的模值 （832）由式(8-30)和式(8-31)还可求出U、I之间的相角

15、差， （833）若取TS900 ，则式(8-28)式(8-33)可进一步化简，进而大大减少了计算机的运算时间。六、三采样值积算法 三采样值积算法是利用三个连续的等时间间隔TS的采样值中两两相乘，通过适当的组合消去t项以求出u、i的幅值和其它电气参数。设在tk+1 后再隔一个TS为时刻tk+2 ，此时的u、i采样值为 （834） （835）上式两采样值相乘，得 （836）上式与式(8-20)相加，得 显然，将式(8-37)和式(8-21)经适当组合以消去tk项，得 若要Ts30o ，上式简化为用Im代替Um（或Um代替Im ），并取0o ，则有 （840） （841）由式(8-39)和式(8-4

16、1)可得 （842）由式(8-27)和式(8-41)，并考虑到，得 （843）由式(8-40)和式(8-41)得 （844）由式(8-42)和式(8-43)得 （845）三采样值积算法的数据窗是2Ts。从精确角度看，如果输入信号波形是纯正弦的，这种算法没有误差，因为算法的基础是考虑了采样值在正弦信号中的实际值。七、傅里叶算法（傅氏算法） 1. 全周波傅里叶算法 根据傅里叶级数，我们将待分析的周期函数电流信号i(t)表示为可用和分别乘式(8-46)两边，然后在t0到t0T积分，得到 （847） （848）每工频周期T采样N次，对式(8-47)和式(8-48)用梯形法数值积分来代替，则得 （849

17、） （850）式中k、ik第k采样及第k个采样值电流n次谐波幅值（最大值）和相位（余弦函数的初相）分别为 （851） （852）写成复数形式有 对于基波分量，若每周采样12点（N=12），则式(8-49)和式(8-50)可简化为 （853） （854）在微机保护的实际编程中，为尽量避免采用费时的乘法指令，在准确度容许的情况下，为了获得对采样结果分析计算的快速性，可用（11/8）近似代替上两式中的，而后1/2和1/8采用较省时的移位指令来实现。全周波傅里叶算法本身具有滤波作用，在计算基频分量时，能抑制恒定直流和消除各整数次谐波，但对衰减的直流分量将造成基频（或其它倍频）分量计算结果的误差。另外用

18、近似数值计算代替积也会导致一定的误差。算法的数据窗为一个工频周期，属于长数据窗类型，响应时间较长。八、解微分方程算法解微分方程算法是假定保护线路分布电容可以忽略，故障点到保护安装处的线路段可用一电阻和电感串联电路，即RL串联模型来表示，于是下述微分方程成立 （865）式中R、L1 分别为故障点至保护安装处线路段的正序电阻和电感，u、i 分别为保护安装处的电压和电流。 1差分法为解得R1和Ll必须有两个方程式。一种方法是取采样时刻tk-1和tk的两个采样值，则有 (867) （868）将 ，代入上两式并联立求解，将得到 （869） （870）其中，Ts为采样间隔。2积分法 用分段积分法对式(8-65)在两段采样时刻tk-2至tk-1和tk-1至tk分别进行积分，得到 （87