微机继电保护算法.pptVIP

类似地可求得输入信号中的余弦分量为 设法消去Kc、Ks，即可滤除直流分量。令：K=Kc+jKs，则 其中： 从而得到补偿系数Kc、Ks： 由于： 进一步简化补偿系数Kc、Ks，得到： 消除衰减直流分量对傅氏算法影响的校正值IR1j和II1j为: 傅氏修正算法2的MATLAB仿真 clear T=0.02;t1=0.02;N=24;Ts=T/N; t=0:Ts:4*T; y=(exp(-t/t1)-cos(2*pi*50*t)) ; %输入信号 figure; subplot(121); plot(t,y,k); xlabel(t/ms); ylabel(y(t)) ; title(输入信号) a=1; i=1:N; bs=sin(2*pi*i /N) ;bc=cos(2*pi*i/N) ; ys=filter(bs,a,y) ; yc=filter(bc,a,y) ; ym=2*abs(complex(ys,yc))./N; [m n]=size(y); y1=y(1:N); y2=y(2:N+1); r=sum(y1)/sum(y2); Ks=(1-r)*sin(2*pi/N)/(1+r^2-2*r*cos(2*pi/N)); Kc=(1-r)*(cos(2*pi/N)-r)/(1+r^2-2*r*cos(2*pi/N)); b=ones(1,N); y0=filter(b,a,y); Ks=[zeros(1,N) ones(1,n-N)*Ks]; Kc=[zeros(1,N) ones(1,n-N)*Kc]; ys1=ys-Ks.*y0; yc1=yc-Kc.*y0; ym1=2*abs(complex(ys1,yc1))./N; subplot(122); plot(t,ym,:k,t,ym1,k); xlabel(t/ms);ylabel(ym,ym1); legend(ym基本算法,ym1改进算法) title(傅立叶算法计算结果) 傅氏修正算法2的效果 8. 傅氏算法的特点和选择 全周傅氏算法 1）利用一个周期内的全部采样值来进行计算，数据窗为一个周期。 2）稳定性好，能滤除所有整次谐波分量，但所需数据窗较长，若计及微机保护判断和保护出口的延时，一般至少需要1～1.5个周期，响应速度较慢。 半周傅氏算法： 1）仅利用半个周期的采样值来进行计算。 2）响应速度快，但滤波能力相对较弱，计算精度比全周傅氏算法稍差（对消除直流分量和偶次谐波的效果比全周傅氏算法有所削弱）。 傅氏算法的选择 各种保护算法中，准确性和快速性是一对互相矛盾的量，选择算法时，既要注意计算精度又要兼顾计算速度。 电力系统的电气信号是以基波及奇次谐波为主，而没有或仅有很小比例的偶次谐波。 响应速度较慢时可考虑半周傅氏算法。 四. 微机保护中的阻抗算法 用于计算阻抗或R、X。如线路距离保护计算测量阻抗。 1. 基于正弦信号计算电阻和电抗的算法 当U、I均为正弦信号时，有： 若取 fs=600Hz，?Ts=300，则上 式可改写为： 同理： 输电线路集中参数模型 fs=600Hz时，对正弦信号，uk-6=－uk，ik-6=－ik，故 求解上述方程得： 应用该算法计算阻抗时，需配合使用一个基波分量的 带通滤波器，以保证正弦信号。 2. 用傅氏算法的实部、虚部计算阻抗 设基波分量的电压、电流的实部、虚部分别为UR、UI，IR、 II，则： 则阻抗为： 故： 上式中的UR、UI，IR、 II，可由傅氏算法求出，并应根据保护类型分别选取。 对反应相间故障的阻抗继电器： UR＝ (UA－UB)R、 (UB－UC)R 、(UC－UA)R UI ＝ (UA－UB)I、 (UB－UC)I、(UC－UA)I IR＝(IA－IB)R、 (IB－IC)R 、(IC－IA)R II ＝(IA－IB)I、 (IB－IC)I 、(IC－IA)I 对反应接地故障的阻抗继电器： UR＝UAR、UBR 、UC；UI ＝UAI、UBI、UCI IR＝(IA+3KI0)R、 (IB+3KI0)R 、 (IC+3KI0)R； II ＝ (IA+3KI0)I、 (IB+3KI0)I 、 (IC+3KI0)I 3. 由电压电流的微分方程计算阻抗 线路发生短路故障时，测量端电流、电压满足以下微分方程： 测量电压u、电流 i 取值与故障类型有关。例如，AB相间短路，u＝ua－ub、i＝ia－ib，单相接地短路，u用故障相电压，i 则为故障相电流加零序补偿电流。 针对两个不同时刻建立如下方程： (此方程有两种求解方法－短数 据窗法与长数据窗法) 1）短数据窗法计算R、L 由图得： 电流、电压的采样值 求解上述二元一次方程组，可得： 短数据窗法特点： 1）t1、t2两个时刻隔一个