% ------------------------------------------------------------------------------ % 1. Profilaxia do caso de simulacao % ------------------------------------------------------------------------------ close all; fclose all; clear all; % ------------------------------------------------------------------------------ % 2. Leitura do caso de simulacao % ------------------------------------------------------------------------------ ABC005_RED; % O sinal original tem 32 amostras por ciclo tempo = matriz(:,1); % Sinal temporal no arquivo IAL = matriz(:,2); % Corrente da fase A no terminal local IBL = matriz(:,3); % Corrente da fase B no terminal local ICL = matriz(:,4); % Corrente da fase C no terminal local IAR = matriz(:,5); % Corrente da fase A no terminal remoto IBR = matriz(:,6); % Corrente da fase B no terminal remoto ICR = matriz(:,7); % Corrente da fase C no terminal remoto % ------------------------------------------------------------------------------ % 3. Configuracao do rele de protecao % ------------------------------------------------------------------------------ % 3.1 Configuracoes sobre o sistema e sobre a amostragem f = 60; % Frequencia do sistema de potencia na = 16; % Numero de amostras por ciclo fa = na*f; % Frequencia de amostragem Ta = 1/fa; % Periodo de amostragem ciclosbuff = 4; % Numero de ciclos armazenados no buffer tambuffer = ciclosbuff*na; % Numero de amostras no buffer iaL_ied = zeros(1,tambuffer); % Tamanho do buffer para corrente na fase A do terminal local ibL_ied = zeros(1,tambuffer); % Tamanho do buffer para corrente na fase B do terminal local icL_ied = zeros(1,tambuffer); % Tamanho do buffer para corrente na fase C do terminal local iaR_ied = zeros(1,tambuffer); % Tamanho do buffer para corrente na fase A do terminal remoto ibR_ied = zeros(1,tambuffer); % Tamanho do buffer para corrente na fase B do terminal remoto icR_ied = zeros(1,tambuffer); % Tamanho do buffer para corrente na fase C do terminal remoto % ------------------------------------------------------------------------------ % 3.2 Especificacao do filtro de entrada do IED % a) Dados do filtro Butterworth fp = 90; % Frequencia maxima da banda de passagem, em [Hz] hc = 3; % Harmonica que se deseja eliminar fs = hc*f;% Frequencia da banda de rejeicao, em [Hz] Amax = 3; % Atenuacao fora da banda de passagem, [dB] Amin = 32; % Atenuacao fora da banda de passagem, [dB] % b) Ordem de um filtro Butterworth % [f_order, wc] = buttord(2*pi*fp/(pi*fa), 2*pi*60*hc/(pi*fa), Amin, Amax); %[f_order, wc] = buttord(2*pi*fp/(2*pi*fa), 2*pi*60*hc/(2*pi*fa), 0.1, Aten); % c) Cria o filtro % [num, den] = butter(f_order, wc); % d) Cria a funcao de transferencia % filtro = tf(num,den); % ------------------------------------------------------------------------------ % 3.3 Configuracoes sobre as funcoes de protecao % a) Sobrecorrente curva = []; % Familia e tipo de curva escolhida (IEEE ou IEC) Ipk = []; % Corrente de pickup (neste caso pode ser em termos de valores primarios Dt = []; % Delta de tempo para coordenacao entre as protecoes % ------------------------------------------------------------------------------ % 4. Filtragem analogica e reamostragem do sinal % a) Filtragem do sinal % IALf = filter(num, den, IAL); % Filtragem do sinal de corrente da fase A do terminal local % IBLf = filter(num, den, IBL); % Filtragem do sinal de corrente da fase B do terminal local % ICLf = filter(num, den, ICL); % Filtragem do sinal de corrente da fase C do terminal local % IARf = filter(num, den, IAR); % Filtragem do sinal de corrente da fase A do terminal remoto % IBRf = filter(num, den, IBR); % Filtragem do sinal de corrente da fase B do terminal remoto % ICRf = filter(num, den, ICR); % Filtragem do sinal de corrente da fase C do terminal remoto IALf = filtro_analogico(0, IAL, tempo, 2*pi*fp, 2*pi*fs, Amin, Amax); % Filtragem do sinal de corrente da fase A do terminal local IBLf = filtro_analogico(0, IBL, tempo, 2*pi*fp, 2*pi*fs, Amin, Amax); % Filtragem do sinal de corrente da fase B do terminal local ICLf = filtro_analogico(0, ICL, tempo, 2*pi*fp, 2*pi*fs, Amin, Amax); % Filtragem do sinal de corrente da fase C do terminal local IARf = filtro_analogico(0, IAR, tempo, 2*pi*fp, 2*pi*fs, Amin, Amax); % Filtragem do sinal de corrente da fase A do terminal remoto IBRf = filtro_analogico(0, IBR, tempo, 2*pi*fp, 2*pi*fs, Amin, Amax); % Filtragem do sinal de corrente da fase B do terminal remoto ICRf = filtro_analogico(0, ICR, tempo, 2*pi*fp, 2*pi*fs, Amin, Amax); % Filtragem do sinal de corrente da fase C do terminal remoto % b) Reamostragem do sinal (como o sinal original possui 32 amostras por ciclo, % basta fazer a decimacao convencional, caso contrario seria necessario o resample % com alguma tecnica de PDS, do tipo "zero padding") cont = 1; for aux=1:length(IALf) if mod(aux,2) tempor(cont) = (aux-1)*Ta; IALfr(cont) = IALf(aux); IBLfr(cont) = IBLf(aux); ICLfr(cont) = ICLf(aux); IARfr(cont) = IARf(aux); IBRfr(cont) = IBRf(aux); ICRfr(cont) = ICRf(aux); cont = cont + 1; end end % ------------------------------------------------------------------------------ % 4. Processamento da protecao (na vida real eh um loop infinito do tipo ] % 'while (1)' % ------------------------------------------------------------------------------ tam = 1; posbuffer = 1; while tam<=length(tempor) % ---------------------------------------------------------------------------- % 4.1 Armazenagem das amostras no buffer do IED (verificando se chegou no % final do tamanho dele (caso contrario tem que comecar a sobrescrever as % amostras mais antigas % ---------------------------------------------------------------------------- if posbuffer>tambuffer posbuffer = 1; end iaL_ied(posbuffer) = IALfr(tam); % Buffer de corrente da fase A no terminal Local ibL_ied(posbuffer) = IBLfr(tam); % Buffer de corrente da fase B no terminal Local icL_ied(posbuffer) = ICLfr(tam); % Buffer de corrente da fase C no terminal Local iaR_ied(posbuffer) = IARfr(tam); % Buffer de corrente da fase A no terminal Remoto ibR_ied(posbuffer) = IBRfr(tam); % Buffer de corrente da fase B no terminal Remoto icR_ied(posbuffer) = ICRfr(tam); % Buffer de corrente da fase C no terminal Remoto % ---------------------------------------------------------------------------- % 4.2 Monta o vetor de correntes para cálculo de Fourier % ---------------------------------------------------------------------------- if posbuffer