基于蒙特卡罗法评估智能电网可靠性研究（Matlab代码实现）

最新推荐文章于 2024-09-21 07:18:39 发布

天天程序猿

最新推荐文章于 2024-09-21 07:18:39 发布

阅读量190

点赞数

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/2301_77414277/article/details/131354777

版权

文章探讨了DG（分布式电源）接入对配电系统可靠性指标的影响，包括发电机组的可靠性、接入容量评估、保护设备协调以及系统指标的调整。通过蒙特卡罗模拟方法，研究了智能电网应用如何改善配电网的可靠性，以IEEE34试验馈线为例，重点测量了自动重合闸器和DG安装对SAIDI、SAIFI、CAIDI和EUE等指标的影响。

摘要由CSDN通过智能技术生成

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章讲解

💥1 概述

配电系统可靠性指标[13]的确立是进行可靠性评估的关键，而 DG 大量接入配电系统，将对可靠性评估指标体系产生如下影响：

（1）需将发电机组的可靠性指标引入配电系统。传统配电系统中的元件主要包括主馈线、分支线、配变及联络开关等，评价各元件可靠性的指标主要有故障率、维修率、平均运行时间及平均修复时间等。DG 接入配电系统后，系统中增加了发电元件，发电机组的可靠性指标就应引入配电系统。故大量接入 DG 的配电系统可靠性评估指标也应包括发电机组的故障率、平均无故障可用小时、启动可靠度等。

（2）需对 DG 的接入容量进行评估。由于 DG 接入配电网后将引发其运行及协调问题，所以需对其接入的容量进行评估，并用安全性指标来评估 DG 的容量可靠度，使系统在引入 DG 后能达到预定的可靠性水平。

（3）需协调系统与 DG 的保护设备及其控制方案。DG 接入系统后，将导致短路电流增大、电压因功率输出与负荷不协调而产生波动，从而干扰电压控制、减少或增加网络损耗等，如果系统与 DG 的保护设备及其控制方案设置不当或不协调，将会削弱系统的供电可靠性和安全性，为此，需引入衡量 DG保护控制对配电系统影响的指标。

（4）需在系统指标中引入评估 DG 影响的指标。传统可靠性评估的系统指标主要有系统停电频率、停电持续时间和平均供电可用率等指标，这些指标依据负荷点指标的加权均值求得，因此 DG 对配电网可靠性的影响将被弱化，或系统指标体现不出 DG 接入的不利影响[14]，故有必要在系统可靠性指标中引入合理评价 DG 影响的指标。

本文采用的系统为 IEEE 34 试验馈线。目的是测量自动重合闸器在可靠性指标SAIDI、SAIFI、CAIDI和EUE（ARs）、分布式发电机（DGs）上的安装。通过蒙特卡罗模拟方法来检验智能电网应用在提高配电网可靠性方面的影响。

📚2 运行结果

部分代码：

switch_834_lines = [20 21 22 31 32];
switch_832_lines = [17 18 19 30 23 24 25 switch_834_lines];
switch_834_elements = [19 20 21 30 31];
switch_832_elements = [switch_834_elements 18 22 23 24 25 29 32];
interaptionXcustomer_number = 0;
interaptionXcustomer_duration = 0;
unservedEnergy = 0;
cur_t = 0;
while cur_t < maxT
% find next failure time
min_t = (duration + 1)*24*365;
fail_line_number = 0;
for i = 1:N
clear indx;
indx = find(downT{i} > cur_t);
if (~isempty(indx))
indx = indx(1);
if (downT{i}(indx) < min_t)
min_t = downT{i}(indx);
fail_line_number = i;
n = indx;
cur_outage_time = upT{i}(indx) - downT{i}(indx);
end
end
end
cur_t = upT{fail_line_number}(n);
in_switch_834 = ~isempty(find(switch_834_lines == fail_line_number));
customer_834 = sum(users(switch_834_elements));
power_834 = sum(power(switch_834_elements));
in_switch_832 = ~isempty(find(switch_832_elements == fail_line_number));
customer_832 = sum(users(switch_832_elements));
power_832 = sum(power(switch_832_elements));