包含UPFC和风电场的分接输电线路差动保护方案（Matlab代码实现）

  💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

 ⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、Simulink仿真实现

---

 ⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

在现代电力系统中，统一潮流控制器（UPFC）和风电场的应用越来越广泛，而在包含这些设备的分接输电线路中，传统的差动保护方案可能面临挑战。下面是一个包含UPFC和风电场的分接输电线路差动保护方案的设计思路：

 1. 差动保护的基本原理
差动保护是通过比较输电线路两端的电流来检测故障。如果两端电流之差超过设定值，则认为线路内部发生了故障。该方法对线路外部故障具有很强的选择性。

2. UPFC的影响
UPFC（Unified Power Flow Controller）是一种灵活交流输电系统（FACTS）设备，它可以调节输电线路的电压、阻抗和相位角，从而实现对电力潮流的控制。由于UPFC可以显著改变输电线路的电气特性，这可能会影响差动保护的正确动作。因此，在设计差动保护方案时，必须考虑UPFC的影响。

 3. 风电场的影响
风电场接入输电线路后，风电机组的发电特性（如间歇性和波动性）会影响线路的电流分布。此外，风电场的功率电子设备也会产生谐波，影响差动保护的动作判断。

4. 差动保护方案设计
针对包含UPFC和风电场的分接输电线路，差动保护方案的设计可以参考以下步骤：

4.1 数据采集和同步
- 使用GPS同步装置（如Phasor Measurement Unit，PMU）获取输电线路两端的电流和电压信号，确保数据的时钟同步。
- 在UPFC和风电场接入点也安装PMU，实时监测电流、电压和功率等参数。

4.2 信号处理和滤波
- 对采集到的电流信号进行滤波，去除风电场功率电子设备产生的谐波。
- 采用自适应滤波器，根据实时工况调整滤波参数，以提高差动保护的抗干扰能力。

4.3 差动保护算法
- 采用改进的比率制动特性（restraining characteristic），以适应UPFC对电流的调节作用。
- 引入暂态分量分析，识别并区分线路内部故障和外部故障。可以采用小波变换等先进的信号处理技术，提取故障特征。

 4.4 保护整定与逻辑
- 设置合理的保护整定值，确保在正常运行和轻微故障情况下差动保护不误动。
- 在差动保护动作逻辑中，增加对UPFC和风电场状态的检测和判断。如果UPFC或风电场发生异常，应及时调整保护整定值或触发备用保护。

4.5 系统测试和仿真
- 使用电力系统仿真软件（如PSCAD/EMTDC、MATLAB/Simulink）进行仿真测试，验证差动保护方案在各种故障情况下的性能。
- 进行现场测试，验证保护方案在实际工况下的可靠性和准确性。

 5. 结论
包含UPFC和风电场的分接输电线路差动保护方案需要综合考虑各种影响因素，采用先进的信号处理和保护算法，确保在复杂工况下能够准确检测故障，提高电力系统的安全性和可靠性。

这个方案旨在提供一种框架，具体的实现细节需要结合实际电力系统的参数和工况进行调整和优化。

📚2 运行结果

部分代码：

%% PLOTTING THE FAULT VOLTAGE AND FAULT CURRENT AT DIFFERENT BUSES
% Calulation of Operating energy and Restraining energy using Wavelet
​
clear
close all
Ts=0.00005;
sim('final_upfc_400kv_wavelet_windfarm_q')
I1 = x5.signals.values;
I1a = I1(:,1);
I1b = I1(:,2);
I1c = I1(:,3);
I4 = x.signals.values;
I4a = I4(:,1);
I4b = I4(:,2);
I4c = I4(:,3);
I7 = x9.signals.values;
I7a = I7(:,1);
I7b = I7(:,2);
I7c = I7(:,3);
​
%% Wavelet of Phase A current at Bus 1
[C,L] = wavedec(I1a,3,'db4');
cA3 = appcoef(C,L,'db4',3);
cD3 = detcoef(C,L,3);
IB1A = wrcoef('a',C,L,'db4',3);
% t = 0:0.0001:1;
sample = 0:1:10000;
figure(1)
subplot(2,1,1)
plot(sample,I1a,'r','linewidth',2.0)
grid on;
title('Sending end phase A current')
xlabel('
Samples
')
ylabel('
Magnitude 
of
 Fault 
Current
at
 phase A
')
subplot(2,1,2)
plot(sample,IB1A,'
linewidth
',2.0)
grid on;
title('
Reconstructed Waveform 
of
 sending 
end
 phase A 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
')
​
%% Calculation of Wavelet at phase B
[C2,L2] = wavedec(I1b,3,'
db4
');
cA3b = appcoef(C2,L2,'
db4
',3);
cD3b = detcoef(C2,L2,3);
IB1B = wrcoef('
a
',C2,L2,'
db4
',3);
% t = 0:0.0001:1;
sample = 0:1:10000;
figure(2)
subplot(2,1,1)
plot(sample,I1b,'
r
','
linewidth
',2.0)
grid on;
title('
Sending 
end
 phase B 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
at
 phase B
')
subplot(2,1,2)
plot(sample,IB1B,'
linewidth
',2.0)
grid on;
title('
Reconstructed Waveform 
of
 sending 
end
 phase B 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
')
​
%% Calculation of Wavelet at phase C
[C4,L4] = wavedec(I1c,3,'
db4
');
cA3c = appcoef(C4,L4,'
db4
',3);
cD3c = detcoef(C4,L4,3);
IB1C = wrcoef('
a
',C4,L4,'
db4
',3);
% t = 0:0.0001:1;
sample = 0:1:10000;
figure(3)
subplot(2,1,1)
plot(sample,I1c,'
r
','
linewidth
',2.0)
grid on;
title('
Sending 
end
 phase C 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
at
 phase C
')
subplot(2,1,2)
plot(sample,IB1C,'
linewidth
',2.0)
grid on;
title('
Reconstructed Waveform 
of
 sending 
end
 phase C 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
')
​
%% Wavelet of phase A current at Bus 4
[C1,L1] = wavedec(I4a,3,'
db4
');
cAR3 = appcoef(C1,L1,'
db4
',3);
cDR3 = detcoef(C1,L1,3);
IB4A = wrcoef('
a
',C1,L1,'
db4
',3);
% t = 0:0.0001:1;
sample = 0:1:10000;
figure(4)
subplot(2,1,1)
plot(sample,I4a,'
r
','
linewidth
',2.0)
grid on;
title('
Receiving 
end
 phase A 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
 Fault 
Current
at
 phase A
')
subplot(2,1,2)
plot(sample,IB4A,'
linewidth
',2.0)
grid on;
title('
Reconstructed Waveform 
of
 receiving 
end
 phase A 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
')
​
%% Wavelet of phase B current at Bus 4
[C3,L3] = wavedec(I4b,3,'
db4
');
cAR3b = appcoef(C3,L3,'
db4
',3);
cDR3b = detcoef(C3,L3,3);
IB4B = wrcoef('
a
',C3,L3,'
db4
',3);
% t = 0:0.0001:1;
sample = 0:1:10000;
figure(5)
subplot(2,1,1)
plot(sample,I4b,'
r
','
linewidth
',2.0)
grid on;
title('
Receiving 
end
 phase B 
current
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
 Fault 
Current
at
 phase B
')
subplot(2,1,2)
plot(sample,IB4B,'
linewidth
',2.0)
grid on;
title('
Reconstructed Waveform 
of
 phase B 
current
at
 receiving 
end
')
xlabel('
Samples
')
ylabel('
Magnitude 
of
Current
')

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]罗涛.基于故障分量的T接输电线路电流差动保护的研究[D].天津大学,2012.​

[2]范栋琛,宋亮亮,崔玉,等.基于数据挖掘的含UPFC输电线路纵差保护方案[J].电子器件, 2021(004):044.

[3]顾威,李兴源,魏巍.基于UPFC的风电场稳定性动态仿真研究[J].电力系统保护与控制, 2010(11):5.DOI:10.3969/j.issn.1674-3415.2010.11.015.

🌈4 Matlab代码、Simulink仿真实现

资料获取，更多粉丝福利，MATLAB|Simulink|Python资源获取