💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
📋📋📋本文目录如下:🎁🎁🎁
目录
🌈4 Simulink仿真、Matlab代码、数据、pdf文档
💥1 概述
"SimPowerSystems模型的1000兆瓦HVDC-MMC互联。本SPS模拟通过使用一个合并MMC模型进行了优化。直流电链已成为国家之间电力交换以及可再生能源(水电厂、近岸风电和太阳能发电场)向电网输电的首选方式。目前有几个这样的系统正在运行,例如INELFE项目(法国-西班牙互联)、Dolwin1项目(近海风电HVDC链),或者计划中的未来项目,比如北部通道(加拿大-美国1090兆瓦直流链)。
我们的示例将说明目前广泛使用的一种典型VSC拓扑结构:模块化多电平变流器(MMC)技术。在这个示例中,MMC变流器使用一个合并模型来模拟每个臂上的36个功率模块。通过这个合并模型,控制系统动态、变流器谐波和循环电流现象都得到了很好的表征。然而,由于只有一个虚拟电容器来代表臂上的36个电容器,该模型假定所有功率模块的电容器电压都是平衡的。与使用每个单独功率模块的两个开关器件和一个电容器的详细模型相比,这个合并模型运行速度要快得多。这个合并模型也非常适合实时仿真。
1.1 简化电网
电网使用 400 kV、50 Hz 等效电压和两个断路器为转换器 1 供电.
1.2 转换器 1
转换器 1 使用 6 个半桥 MMC 模块实现,每个模块代表 36 个电源模块。此自定义 SPS 模块使用开关功能模型,其中仅使用一个等效模块来表示所有电源模块。控制信号是一个二元向量 [Nin, Nbl],其中 Nin 表示插入模块的数量,Nbl 表示处于阻塞状态的模块数量。输出 Vc(电容器电压)只有一个元素,并给出电容器模块的平均值。
1.3 直流电路
在直流电路子系统中,将会发现电缆的简化模型,以及使用直流电源和理想开关建模的第二个转换器。您还将找到一个开关,用于在电缆上应用故障。
1.4 控制器
控制器子系统包含操作互连所需的各种控制系统。它包括以下子系统:有功和无功功率调节器-直流电压调节器-带前馈的dq电流调节器-PLL和测量子系统-用于控制转换器1(Converter1 )半桥MMC的PWM发电机。在文件夹中还可以找到对转换器1的各种操作模式进行编程的排序器区域。
1.5 示波器和测量
该子系统包含用于在仿真期间观察多个信号的所有示波器。功率和基波电压计算也在子系统中进行。
基于模块化多电平换流器(Modular Multilevel Converter, MMC)技术和电压源换流器(Voltage Source Converter, VSC)的高压直流(High Voltage Direct Current, HVDC)系统,是现代电力传输领域的一项先进技术,尤其适用于远距离、大容量的电力传输以及电网互联。这种技术结合了MMC与VSC的优势,提供了更高效率、更灵活可控性及更低环境影响的解决方案。下面简要介绍其工作原理和特点:
MMC技术基础
MMC是一种先进的电力电子变换器技术,主要特点是将传统的两电平或三电平换流器结构分解为多个子模块级联而成,每个子模块都能独立调节输出电压,从而实现多电平输出。这样做的优势包括但不限于:
- 输出电压质量高:由于采用多电平输出,可以生成接近正弦波的电压,减少谐波污染。
- 损耗低:子模块间的电压差小,降低了开关损耗。
- 模块化设计:易于维护和扩展,单个子模块故障不影响整个系统的运行。
- 灵活性高:可适应不同的电压等级和功率需求。
VSC技术介绍
VSC是一种通过直流侧电容提供能量进行交流-直流或直流-交流转换的换流器。与传统的基于晶闸管的电流源换流器(CSC)不同,VSC具有以下优点:
- 双向功率控制:能够独立控制有功功率和无功功率,适用于电网支撑和可再生能源接入。
- 快速响应:对电网故障能迅速做出反应,增强系统的稳定性。
- 占地面积小:相比传统换流站,VSC系统更加紧凑。
结合MMC与VSC的HVDC系统
将MMC技术应用于VSC型HVDC系统中,形成MMC-VSC HVDC系统,进一步增强了HVDC系统的性能:
- 更高的效率和功率密度:MMC的低损耗特性和VSC的高效控制策略相结合,提高了整体系统的能效。
- 更好的电网兼容性:输出电压质量高,易于与各种电网连接,特别适合接入弱电网或可再生能源发电系统。
- 更强的可控性:VSC的独立功率控制能力加上MMC的精细电压调节,使得系统能在各种操作条件下维持稳定运行,包括有功功率传输、无功补偿以及黑启动等复杂应用场景。
- 环境友好:减少了电磁干扰和 audible noise,降低了对周围环境的影响。
综上所述,基于MMC和VSC技术的HVDC系统在提高电力传输效率、增强系统灵活性和稳定性方面展现出显著优势,是未来智能电网建设和大规模可再生能源并网的关键技术之一。
📚2 运行结果
可视化matlab代码:
close all
%
load CompareResultsData
%
h1=figure;
set(h1,'Name','DC');
ScreenS=get(0,'Screensize');
set(h1,'Position',[ScreenS(3)*0.01 ScreenS(4)*0.52 ScreenS(3)*0.32 ScreenS(4)*0.37])
subplot(2,1,1)
plot(results.time3,results.signals(1).values*1e-3,'b', ...
results.time3,results.signals(11).values*1e-3,'r')
ylabel('(kV)')
grid
title('Vdc Conv1')
axis([0 10 0 700])
legend('SwD','Agg')
subplot(2,1,2)
plot(results.time3,results.signals(2).values,'b', ...
results.time3,results.signals(12).values,'r')
ylabel('(A)')
xlabel('(s)')
grid
title('Idc Conv1')
legend('SwD','Agg')
axis([0 10 -2000 1000])
%
h1=figure;
set(h1,'Name','PQ');
ScreenS=get(0,'Screensize');
set(h1,'Position',[ScreenS(3)*0.34 ScreenS(4)*0.52 ScreenS(3)*0.32 ScreenS(4)*0.37])
subplot(2,1,1)
plot(results.time2,results.signals(3).values,'b', ...
results.time2,results.signals(13).values,'r')
ylabel('(MW)')
grid
title('Active Power')
axis([0 10 -500 1500])
legend('SwD','Agg','Location','southeast')
subplot(2,1,2)
plot(results.time2,results.signals(4).values,'b', ...
results.time2,results.signals(14).values,'r')
ylabel('(Mvar)')
xlabel('(s)')
grid
title('Reactive Power')
axis([0 10 -300 300])
legend('SwD','Agg','Location','southeast')
%
%
h1=figure;
set(h1,'Name','Iprim');
ScreenS=get(0,'Screensize');
set(h1,'Position',[ScreenS(3)*0.67 ScreenS(4)*0.52 ScreenS(3)*0.32 ScreenS(4)*0.37])
plot(results.time1(1:25000),results.signals(5).values(1:25000),'b')
hold on
plot(results.time1(1:25000),results.signals(15).values(1:25000),'r')
ylabel('(A)')
xlabel('(s)')
grid
title('Primary current, phase A')
axis([0.35 0.45 -250 250])
text(0.38,225,'(during capacitors energization)')
legend('SwD','Agg','Location','southeast')
%
h1=figure;
set(h1,'Name','Vcap');
ScreenS=get(0,'Screensize');
set(h1,'Position',[ScreenS(3)*0.01 ScreenS(4)*0.06 ScreenS(3)*0.32 ScreenS(4)*0.37])
subplot(2,2,1)
plot(results.time1(225000:250000),results.signals(6).values)
ylabel('(V)')
axis([9.9 9.98 16000 20000])
grid
title('Vcap 4 modules; SwD')
subplot(2,2,2)
plot(results.time1(225000:250000),results.signals(16).values,'r')
axis([9.9 9.98 16000 20000])
ylabel('(V)')
grid
title('Vcap average: Agg')
subplot(2,2,[3:4])
plot(results.time1(225000:250000),results.signals(7).values,'b', ...
results.time1(225000:250000),results.signals(17).values,'r')
ylabel('(V)')
axis([9.9 9.98 16000 20000])
xlabel('(s)')
grid
title('Capacitors average voltage, upper arm phA')
legend('SwD','Agg','Location','northeast')
%
h1=figure;
set(h1,'Name','Iarm');
ScreenS=get(0,'Screensize');
set(h1,'Position',[ScreenS(3)*0.34 ScreenS(4)*0.06 ScreenS(3)*0.32 ScreenS(4)*0.37])
plot(results.time1(225000:250000),results.signals(8).values,'b', ...
results.time1(225000:250000),results.signals(18).values,'r')
ylabel('(A)')
xlabel('(s)')
grid
title('Iarm, upper arm, phase A')
axis([9.9 10.0 -2000 2000])
legend('SwD','Agg','Location','southeast')
%
h1=figure;
set(h1,'Name','Iprim');
ScreenS=get(0,'Screensize');
set(h1,'Position',[ScreenS(3)*0.67 ScreenS(4)*0.06 ScreenS(3)*0.32 ScreenS(4)*0.37])
plot(results.time1(225000:250000),results.signals(5).values(225000:250000),'b', ...
results.time1(225000:250000),results.signals(15).values(225000:250000),'r')
ylabel('(A)')
xlabel('(s)')
grid
title('Primary current phase A')
axis([9.9 10.0 -2500 2500])
text(9.935,2300,'(in steady-state)')
legend('SwD','Agg','Location','southeast')
%
🎉3 参考文献
文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。