电力系统潮流+风光+水电（IEEE33）【Matlab代码实现】

 👨‍🎓个人主页：研学社的博客   

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

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

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

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

目录

💥1 概述

📚2 运行结果

 2.1 案例

 2. 2 结果展现

🎉3 参考文献

🌈4 Matlab代码实现

---

​

💥1 概述

信息可视化的应用为电力行业带来了巨大的希望，但其潜力尚未被可视化社区充分挖掘。先前关于电力系统可视化的工作仅限于在地理布局之上描绘原始或处理过的信息。很少有人致力于可视化电网的物理特性，这最终决定了电力基础设施的状况和稳定性。电力系统的发展还处于数据可视化阶段，主要以图形的形式表示数据。挖掘数据间潜在规律和关系的可视化技术还相当缺乏。在海量电力数据中，只有合理安排和去除噪声，才能为可视化提供可靠的数据源，为电力系统管理者提供更有价值的策略服务，科学可视化是一门多学科的技术，它融合了计算机图形学、数据挖掘、人机界面理论，并与应用科学相结合，为工程领域有价值的数据提供直观的展示方式，为操作者探索和研究提供手段。物理现象。近年来可视化技术在电力系统中的应用越来越广泛，提供的展示手段越来越多，视觉效果也越来越好，但这些都不能满足电力系统调度员的高要求。如何让可视化技术更加实用，是一个亟待解决的难题。本文提出了几种可视化实际研究的方法。基于已在电力系统中部署的应用，分析了如何提高电力系统科学可视化的实用性，介绍了如何利用可视化技术在电力系统中进行高级应用[3]。

电力行业对可视化的应用长期以来一直没有得到充分利用。AREVA的能源管理系统（EMS）和PowerWorld的模拟器代表了两个最广泛使用的工具，在行业内提供可视化支持。然而，它们的主要功能主要是在地理地图等地理布局的基础上用颜色或图标来描述收集或处理的信息。这种以演示为导向的可视化方法，虽然很流行，但却不能利用可视化本身自然获得的分析优势。

然而，自 2000 年以来，许多电力系统控制中心采用了更先进的可视化技术，尤其是用于广域显示。例如，现在使用颜色轮廓来显示大区域的电压幅度变化很常见。动画箭头可用于显示线路流向，而动态大小的饼图用于指示接近或超过其限制或停止服务的输电线路和变压器。图1提供了使用这些技术的一个示例。

​

📚2 运行结果

 2.1 案例

【节点数可以换成其他的】

以某地电网为例对本文模型进行可视化分析。网络如图4.1所示，对IEEE30节点进行建模，电网线路、水力发电机、风机、光伏和负荷（包括住宅区、工厂、第三产业）也被纳入其中。网络参数表写在附录L1矩阵中；住宅区、工厂、第三产业负荷数据；光伏、风机、水力发电机组出力数据；有功功率和无功功率相关性见附录。

                   ​

                                           图4.1 网络结构图

 2. 2 结果展现

针对各机组发电功率和各类负荷，在Matlab2018a平台上进行仿真，得到居民、工厂、学校负荷，风机、光伏、风能出力如图4.2.1所示。

​

                         图4.2.1 各种负荷需求及各类发电机组出力

然后以节点7为松弛节点，得到各点电压，如图4.2.2所示：

下面是24小时的可视化情况（只选取其中几个时刻，就不24张图都截取了）：

部分代码：

% Power line impedance
Zb=diag([0.00782+0.00212i, 0.01564+0.00424i, 0.01173+0.00318i, 0.01173+0.00318i, 0.013685+0.00371i, 0.01173+0.00318i...
    0.01173+0.00318i, 0.01173+0.00318i, 0.00782+0.00212i, 0.00782+0.00212i, 0.01564+0.00424i, 0.01173+0.00318i...
    0.01173+0.00318i, 0.013685+0.00371i, 0.01173+0.00318i, 0.01173+0.00318i, 0.01173+0.00318i, 0.00782+0.00212i...
    0.00782+0.00212i, 0.01564+0.00424i, 0.01173+0.00318i,0.01173+0.00318i, 0.013685+0.00371i, 0.01173+0.00318i...
    0.01173+0.00318i, 0.01173+0.00318i, 0.00782+0.00212i, 0.00782+0.00212i, 0.01564+0.00424i, 0.01173+0.00318i...
    0.01173+0.00318i, 0.013685+0.00371i, 0.01173+0.00318i]);

% Total power line admittance
Ysh=zeros(node_num,1);

% Slack bus
V0=5;
V0_vect=V0*ones(node_num,1);

% Admittance defined loads
Ycar=zeros(node_num,1);

%% Network Structure

% line aggregated admittance
Y_totNodes= zeros(node_num,1) ;
%%  Generation and Loads

hours= 1:1:24;
num_hours=length(hours);

% Load and generation profiles during 24 hours of simulation (step time =
% 1hour)
%%% Residential

peak_factor=0.12;

res = peak_factor*[0.4,0.25,0.2,0.25,0.40,  0.5,  0.65,0.7,0.65,0.65,...
    0.6,0.6,0.6,0.6,0.65,0.7,  0.8, 0.9, 1, 0.95, 0.9, 0.8, 0.65, 0.55];

% Industrial
ind=[0.2,0.2,0.2,0.2,0.2, 0.4, 0.6, 0.9,0.9,0.9,0.9,0.9, 0.7, ...
     0.9,0.9,0.9,0.9, 0.6, 0.4, 0.2,0.2,0.2,0.2,0.2];

 % Tertiary
ter=[0.35,0.3,0.3,0.35,0.4, 0.5, 0.6, 0.7, 0.8, 0.78, 0.76, 0.74, 0.72, 0.7,... 
     0.68, 0.66, 0.64, 0.52, 0.45,0.43,0.42,0.4,0.4,0.4];

% Photovoltaics
PV=[    0,    0,    0,   0,   0,   0, 0.4, 0.7, 0.88, 0.95, 0.988,  1, ... 
    0.988, 0.95, 0.88, 0.7, 0.4,   0,    0,   0,    0,  0,  0,  0];

% Wind Turbine
WT = -0.1 + rand(1, num_hours)*1.1; 

% Hydrauelectric power plant
hydr = ones(1,num_hours);

% Base amplitudes
S_base=1000; % kVA, Base apparent power
V_base=10; % kV, Base voltage
I_base=S_base/(sqrt(3)*V_base); % A, Base current

% Thermal limit
Ilim=60; % A
Ilim_pu=Ilim/I_base*ones(node_num,1);

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]梁晓伟,朱琳艳,胡吕龙,王凯,丁建顺.电力系统的潮流计算和可视化技术研究[J].粘接,2022,49(11):161-164.

[2]李响,邢文涛.基于SVG的可视化技术在电力系统中的应用研究[J].黑龙江电力,2022,44(05):433-435+452.DOI:10.13625/j.cnki.hljep.2022.05.012.

🌈4 Matlab代码实现