MMC七电平整流器MATLAB建模与故障分析-育师

1. MMC七电平整流器模型概述

模块化多电平换流器（Modular Multilevel Converter, MMC）作为高压直流输电（HVDC）领域的核心技术，近年来在电力电子领域获得了广泛应用。七电平拓扑结构因其在输出波形质量和器件数量之间的良好平衡，成为工程实践中的常见选择。

在实际工程应用中，MMC七电平整流器的主要优势体现在三个方面：首先，通过多电平输出特性，能够显著降低输出电压的谐波含量；其次，模块化设计便于维护和扩展；最后，采用分布式电容设计，单个子模块故障不会导致系统崩溃。

注意：虽然现代MATLAB版本（2020b之后）提供了更完善的MMC建模工具，但在2019及以下版本中，我们需要手动搭建模型并配置参数，这对理解MMC工作原理非常有帮助。

2. MATLAB环境准备与工具配置

2.1 软件版本要求与工具箱检查

本实验基于MATLAB 2019b版本进行，兼容2016a至2019b的各版本。关键依赖是SimPowerSystems工具箱（现更名为Simscape Electrical），可通过以下命令验证是否安装：

>> ver('simscape') >> ver('simscape electrical')

如果未显示相关工具箱信息，需要通过MATLAB的"附加功能"菜单进行安装。安装完成后，建议运行以下测试脚本验证功能完整性：

% 测试SimPowerSystems基础功能 try powergui; vsrc = simscape.powersys.electrical.sources.acVoltageSource; disp('SimPowerSystems工具箱功能正常'); catch error('SimPowerSystems工具箱存在问题，请检查安装'); end

2.2 基础参数计算与设置

在搭建模型前，需要确定几个核心参数：

系统额定电压：10kV（线电压）
频率：50Hz
子模块电容：根据经验公式 C = P/(6NfΔV) 计算
- 其中P为额定功率，N为电平数，f为频率，ΔV为允许电压波动

以10MVA系统为例，七电平结构(N=7)的电容计算过程如下：

P = 10e6; % 额定功率10MW N = 7; % 七电平 f = 50; % 频率50Hz deltaV = 0.1; % 允许10%电压波动 C = P/(6*N*f*deltaV*1e3) % 电容值计算(单位mF)

计算结果约为476μF，实际工程中可选择470μF标准电容。

3. 主电路建模与参数配置

3.1 三相电源模块设置

在Simulink中新建模型后，从Simscape Electrical库中找到"Three-Phase Source"模块。关键参数设置包括：

线电压：10kV
频率：50Hz
内阻抗：根据系统短路容量计算
相角：0, -120, 120度

配置示例代码：

% 创建三相电源参数结构体 sourceParams.Voltage = 10e3; % 线电压10kV sourceParams.Frequency = 50; % 50Hz sourceParams.PhaseAngle = [0 -120 120]; % 三相相位 sourceParams.ShortCircuitLevel = 100e6; % 短路容量100MVA

3.2 MMC子模块设计与级联

七电平MMC需要6个子模块/相，采用半桥结构。每个子模块包含：

两个IGBT模块
一个直流电容
并联电阻（用于均压）

子模块参数设置要点：

IGBT开关频率：通常1-2kHz
电容值：按前文计算的470μF
电阻值：R = 1/(2πfC)≈6.8kΩ

在Simulink中，可以使用"Half-Bridge Converter"模块为基础构建子模块：

% 子模块参数配置 submodule.SwitchDevice = 'IGBT/Diodes'; % 开关器件类型 submodule.Capacitance = 470e-6; % 电容值 submodule.Rs = 1e-3; % 开关电阻 submodule.Ron = 1e-3; % 导通电阻

4. 控制策略实现

4.1 载波移相调制(CPS-SPWM)

七电平系统需要6个相位互差30度的三角载波。调制波生成代码如下：

fs = 10e3; % 采样频率10kHz t = 0:1/fs:0.1; % 时间向量 m = 0.9; % 调制比 f0 = 50; % 基波频率 % 生成调制波 modWave = m*sin(2*pi*f0*t); % 生成6个移相载波 carriers = zeros(6,length(t)); for k = 1:6 phaseShift = (k-1)*pi/3; % 60度间隔 carriers(k,:) = sawtooth(2*pi*1e3*t + phaseShift, 0.5); end

4.2 电压平衡控制

子模块电容电压平衡是MMC稳定运行的关键。采用排序算法实现电压平衡：

function [gateSignals] = balanceControl(capVoltages, armCurrent) [~, idx] = sort(capVoltages); if armCurrent > 0 % 充电时投入电压最低的子模块 gateSignals = zeros(size(capVoltages)); gateSignals(idx(1)) = 1; else % 放电时投入电压最高的子模块 gateSignals = zeros(size(capVoltages)); gateSignals(idx(end)) = 1; end end

5. 单相接地故障设置与分析

5.1 故障模块参数配置

使用"Three-Phase Fault"模块设置单相接地故障，关键参数包括：

故障类型：A-G（A相接地）
故障电阻：0.1Ω
故障电感：1mH
起始时间：0.05s
持续时间：0.1s

配置代码示例：

fault.Resistance = 0.1; % 故障电阻 fault.Inductance = 1e-3; % 故障电感 fault.InitialStatus = 0; % 初始状态(正常) fault.SwitchingTimes = [0.05 0.15]; % 切换时间 fault.Phases = [1 0 0]; % A相故障

5.2 故障特征分析

当发生单相接地故障时，系统会出现以下特征：

故障相电压降低（接近零）
非故障相电压升高（√3倍）
零序电流出现
直流侧出现二倍频纹波

可通过以下代码提取故障特征：

% 分析故障波形 [Vabc, Iabc] = extractScopeData(simout); V0 = (Vabc(:,1) + Vabc(:,2) + Vabc(:,3))/3; % 零序电压 I0 = (Iabc(:,1) + Iabc(:,2) + Iabc(:,3))/3; % 零序电流 % 计算THD thdA = thd(Vabc(:,1), fs);

6. 仿真结果分析与验证

6.1 正常工况波形

七电平MMC在正常运行时，典型特征包括：

输出电压具有7个离散电平
线电压THD<5%
子模块电容电压波动<10%

使用FFT分析输出波形质量：

[P1,f1] = periodogram(Vab, [], [], fs, 'power'); THD = sqrt(sum(P1(2:end)))/P1(1); % 计算总谐波畸变率

6.2 故障工况保护策略

针对单相接地故障，建议采取以下保护措施：

快速检测算法（基于零序分量）
故障后0.1s内闭锁PWM
直流断路器动作
重启时的预充电流程

保护逻辑实现示例：

function [blockSignal] = protectionLogic(V0, I0, threshold) persistent faultCounter; if isempty(faultCounter) faultCounter = 0; end if abs(V0) > threshold || abs(I0) > threshold faultCounter = faultCounter + 1; else faultCounter = max(0, faultCounter - 1); end blockSignal = faultCounter > 5; % 连续5个周期超过阈值才动作 end

7. 工程实践中的经验分享

在实际项目应用中，有几个容易忽视但至关重要的细节：

参数初始化问题：
- 子模块电容需要预充电至额定电压
- 使用以下代码实现软启动：
```
for k = 1:100 set_param('MMC_Model/Submodule'+string(k), 'InitialVoltage', k/100*Vdc); end
```
仿真步长选择：
- 开关频率1kHz时，建议步长≤10μs
- 使用变步长求解器时，设置最大步长为50μs

实时监测技巧：

% 在仿真过程中实时显示关键参数 add_exec_event_listener('MMC_Model', 'PostOutputs', @(src,evt)disp(... ['Time: ' num2str(evt.Time) ' | CapVolt: ' num2str(mean(capVoltages))]));