《测量阻抗的计算》由会员分享，可在线阅读，更多相关《测量阻抗的计算（44页珍藏版）》请在金锄头文库上搜索。
1、硬件原理框图数据采集系统开关量输入系统电源系统uiCPU主系统开关量输出系统接口系统接口系统跳合闸输出键盘显示QF1直接计算法v微机保护通过硬件电路采集到保护安装处的三相电压和三相电流的离散值，经过（数字滤波）、算法后，得到各相电压、电流的相量值，进一步代入阻抗的计算式，求出测量阻抗，然后与整定阻抗进行比较，判断故障的范围。v下面以单相系统的电压、电流为例，讨论计算测量阻抗的方法。2单相系统采集的数据v设数据采集系统采集到的电压、电流分别为：v对采样数据的处理运算，包括数字滤波和算法，现分别讨论如下。3数字滤波作用：距离保护反映的是工频信号，短路电压电流中既含有工频信号，又含有衰减非周期分量、
2、谐波分量和高频行波信号。滤波的作用就是滤除非周期分量、 谐波分量和高频行波信号，保留工频信号。数字滤波应与继电保护算法配合工作，根据算法的需要进行配置。例如当采用傅氏算法时，由于傅氏算法本身具有滤除谐波的能力，所以就不需要再配置滤除谐波的数字滤波器，只需要滤除非周期分量即可；解微分方程算法中，只需要滤除高频分量。4v在微机保护中，常用的数字滤波器为FIR滤波器，其差分方程为: 其中x(n)为输入序列，y(n)为输出序列，h(n)为滤波器的单位冲击响应。该方程表明，数字滤波器的输出等于其输入与滤波器单位冲击响应的卷积。 FIR数字滤波器5交叉相乘后相加x(n)h(n)h(1)h(K)x(1)x(
3、K)y(K)y(n)输入序列输出序列单位冲击响应（滤波器一旦设计完成，该部分是不变的常数）FIR数字滤波器的实现6交叉相乘后相加x(n)h(n)h(1)h(K)x(1)x(K)y(K)y(n)y(K+1)x(K+1)x(2)FIR数字滤波器的实现7交叉相乘后相加x(n)h(n)h(1)h(K)x(1)x(nK1)y(n)x(2)x(n)y(n)N项FIR数字滤波器的实现x(nK2)8在FIR的一般表达式：令则公式变为： 这种滤波器称为差分滤波器。 K为滤波器的阶数。（K阶差分滤波）差分滤波器9K原始序列滤波结果两者相减两者相减差分滤波器的实现10差分滤波器的原理v滤波的原理： 设输入信号包括基
4、波和m次谐波 若KTs正好等于m次谐波周期的整倍数，则相减以后就能够消除m次谐波。11差分滤波器的原理12f/f1A0m2m3m幅频特性中为0的点，就是能够滤除的谐波次数差分滤波器的特性13差分滤波器的特点（1）能滤除直流分量。（2）当m=1时，滤去直流，基波及各次谐波输出为零，可做保护增量元件。（3）运算量小。（4）数据窗长（缺点）。14NTsnTs（nN）Tsy(n)x(n)增量元件15测量阻抗的计算方法 在微机保护中，测量阻抗的计算方法有多种，此处介绍最常用的两种，即： 傅氏算法 解微分方程算法首先讨论傅氏算法。16傅氏算法 由数学中傅氏级数的理论可知，如果x(t) 是一个周期函数，且满
5、足一定的条件（对电力系统的信号来说，该条件都满足），则x(t) 可以展开为傅立叶级数，即用图形表示为17tttttttX(t)X0X1X2X3X4X518XmsXmcXmm傅氏算法19vXm 和 m 可以由 Xms 和 Xmc 求出v求 Xms 和 Xmc的方法 : 根据傅立叶级数的理论， Xms 和 Xmc 可以表示为:傅氏算法20v在微机保护中，输入函数x(t)的确切表达式是不知道的，只知道它在一系列离散点上的采样值，所以并不能直接求出积分值，这时可以应用积分的近似计算法来近似地计算上述积分。v由数学知识，常用的近似积分算法有矩形积分法和梯形积分法两种。v下面以矩形积分法为例，进行讨论傅氏
6、算法21sin mttx(t)tnTstx(t) sin mtData WindowConcept of rectangle method矩形积分法图解22sin mttx(t)tnTsx(t) sin mttXms will be different for different nTs, So do Xmc数据窗移动23v由图可以得到面积的计算公式24v当 m=1，可以算出基波，这时公式变为 ：25若采样频率为 600Hz, 即 N=12, 则公式变为26271020t (ms)举例28n012345678A105000866710000866750000-5000-8667A3033330
7、-3333033330-33330A501000-17322000-173210000-10001732Samp.0933369358667693593330-9333-6935n91011121314151617A1-10000-8667-50000500086671000086675000A333330-3333033330-333303333A5-20001732-100001000-17322000-17321000Samp.-8667-6935-933309333693586676935933329For n=12，粉红框为数据窗3031实际值应该是 10000, 产生误差的原因是
8、0.8667,0.3333等小数的截断误差.32For n=13:3334The actual number should be 10000,phase angle should be 30, the error comes from round off of 0.8667,0.3333 etc.35For n=14:3637The actual number should be 10000,phase angle should be 60, the error comes from round off of 0.8667,0.3333 etc.The results of X1s、X1c a
9、nd 1 will change with n, but X1 will not change.38Conclusions The result of X1s, X1c and 1 will change with the number n X1 will not change with the number n. 1就是计算窗口的起始点（或结束点）对应的角度。39距离保护的实现方法v应用傅氏算法，可以计算出电压、电流的实部和虚部（或幅值和相位），进一步可以计算出直角坐标形式或极坐标形式表示的测量阻抗。40解微分方程算法基本原理 方法一式中 、 、 都是可测量、可计算的，未知 、 两个方程 ，两次采样时刻为 、uiR1L141 用D表示 则用计算机处理时，D可用差分近似处理，方法是 和 分别为两个相邻的采样瞬间的中间值如图所示：解微分方程算法42此时有法二：R-L模型的积分算法分别在两个不同的时段内积分，而得： 解微分方程算法43 为积分时间长度， 、 则为两个不同的积分起始时刻其余各项积分用梯形法则近似求得联立求解得两个未知数 、当 取得足够大，则兼有一定的滤波作用，从而可抑制高频分量，但它所需数据窗要相当长。 、 求出后，测量阻抗容易得出。解微分方程算法44
