硬件原理框图   直接计算法 微机保护通过硬件电路采集到保护安装处的三相电压和三相电流的离散值，经过（数字滤波）、算法后，得到各相电压、电流的相量值，进一步代入阻抗的计算式，求出测量阻抗，然后与整定阻抗进行比较，判断故障的范围。 下面以单相系统的电压、电流为例，讨论计算测量阻抗的方法。 单相系统采集的数据 设数据采集系统采集到的电压、电流分别为：    对采样数据的处理运算，包括数字滤波和算法，现分别讨论如下。 数字滤波 作用：距离保护反映的是工频信号，短路电压电流中既含有工频信号，又含有衰减非周期分量、谐波分量和高频行波信号。滤波的作用就是滤除非周期分量、 谐波分量和高频行波信号，保留工频信号。 数字滤波应与继电保护算法配合工作，根据算法的需要进行配置。例如当采用傅氏算法时，由于傅氏算法本身具有滤除谐波的能力，所以就不需要再配置滤除谐波的数字滤波器，只需要滤除非周期分量即可；解微分方程算法中，只需要滤除高频分量。 FIR数字滤波器 在微机保护中，常用的数字滤波器为FIR滤波器，其差分方程为:       其中x(n)为输入序列，y(n)为输出序列，h(n)为滤波器的单位冲击响应。该方程表明，数字滤波器的输出等于其输入与滤波器单位冲击响应的卷积。     在FIR的一般表达式：    令 则公式变为：       这种滤波器称为差分滤波器。     K为滤波器的阶数。（K阶差分滤波） 差分滤波器的原理 滤波的原理：    设输入信号包括基波和m次谐波       若KTs正好等于m次谐波周期的整倍数，则相减以后就能够消除m次谐波。  差分滤波器的原理              差分滤波器的特点  （1）能滤除直流分量。 （2）当m=1时，滤去直流，基波及各次谐波输出为零，可做保护增量元件。 （3）运算量小。 （4）数据窗长（缺点）。 测量阻抗的计算方法     在微机保护中，测量阻抗的计算方法有多种，此处介绍最常用的两种，即：        傅氏算法        解微分方程算法 首先讨论傅氏算法。 傅氏算法    由数学中傅氏级数的理论可知，如果x(t) 是一个周期函数，且满足一定的条件（对电力系统的信号来说，该条件都满足），则x(t) 可以展开为傅立叶级数，即  傅氏算法 傅氏算法 Xm 和 αm 可以由 Xms  和 Xmc 求出 求 Xms  和 Xmc的方法 :    根据傅立叶级数的理论， Xms  和 Xmc 可以表示为: 傅氏算法 在微机保护中，输入函数x(t)的确切表达式是不知道的，只知道它在一系列离散点上的采样值，所以并不能直接求出积分值，这时可以应用积分的近似计算法来近似地计算上述积分。 由数学知识，常用的近似积分算法有矩形积分法和梯形积分法两种。 下面以矩形积分法为例，进行讨论 由图可以得到面积的计算公式       当  m=1，可以算出基波，这时公式变为 ：  若采样频率为 600Hz, 即 N=12, 则公式变为     Conclusions       The result of X1s, X1c and α1 will change with the number n  X1 will not change with the number n.  α 1就是计算窗口的起始点（或结束点）对应的角度。   距离保护的实现方法 应用傅氏算法，可以计算出电压、电流的实部和虚部（或幅值和相位），进一步可以计算出直角坐标形式或极坐标形式表示的测量阻抗。 解微分方程算法 基本原理  方法一   式中     、   、    都是可测量、可计算的，未知      、      两个方程  ，两次采样时刻为     、 解微分方程算法                                         用D表示        则    用计算机处理时，D可用差分 近似处理，方法是     和    分别 为两个相邻的采样瞬间的中间 值如图所示： 解微分方程算法 此时有      法二：R-L模型的积分算法 分别在两个不同的时段内积分，而得：                                                                                                                               解微分方程算法      为积分时间长度，   、    则为两个不同的积分起始时刻     其余各项积分用梯形法则近似求得 联立求解得两个未知数      、 当     取得足够大，则兼有一定的滤波作用，从而可抑制高频分量，但它所需数据窗要相当长。            、         求出后，测量阻抗容易得出。 实际值应该是 10000, 产生误差的原因是 0.8667,0.3333等小数的截断误差. For n=13: The actual number should be 10000,phase angle should be 30, the error comes from round off of 0.8667,0.3333 etc. For n=14: The 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 and α1 will change with n, but X1 will not change. u i R1 L1 * * 数据采集系统 开关量 输入系统 电源 系统 u i CPU 主系统 开关量 输出系统 人机 接口系统 通信 接口系统 跳合闸输出 键盘 显示 QF 交叉相乘后相加 x(n) h(n) h(1) h(K) x(1) x(K) y(K) y(n) 输入序列 输出序列 单位冲击响应（滤波器一旦设计完成，该部分是不变的常数） FIR数字滤波器的实现 交叉相乘后相加 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数字滤波器的实现 交叉相乘后相加 x(n) h(n) h(1) h(K) x(1) x(n－K＋1) y(n) x(2) x(n) y(n) N项 FIR数字滤波器的实现 x(n－K＋2) 差分滤波器 K 原始序列 滤波结果 两者相减 两者相减 差分滤波器的实现 f/f1 A 0 m 2m 3m 幅频特性中为0的点，就是能够滤除的谐波次数 差分滤波器的特性 NTs nTs （n－N）Ts y(n) x(n) 增量元件 用图形表示为 t t t t t t t X(t) X0 X1 X2 X3 X4 X5 Xms Xmc Xm αm sin mωt t x(t) t nTs t x(t) sin mωt Data Window Concept of rectangle method 矩形积分法图解 sin mωt t x(t) t nTs x(t) sin mωt t Xms will be different for different nTs, So do Xmc 数据窗移
测量阻抗的计算.ppt