我最近在 MATLAB/Simulink 里重新搭了一套三相 STATCOM 无功补偿仿真模型从主电路搭到控制环调了一周才把动态响应和电压波形调到满意。STATCOM全称静止同步补偿器是现在配电网、风电光伏场站里非常常见的动态无功补偿装置。它跟老式电容器组和 SVC 最大的不同是能连续、快速、双向地输出容性或感性无功不管系统电压怎么晃都能把母线电压和功率因数稳住。用 Simulink 做这套仿真主要解决三件事验证 STATCOM 控制策略行不行、观察各种工况下的暂态波形、把所有 PI 参数和主电路参数预先调出一个可靠初值免得日后去现场瞎试。这篇文章适合电气工程学生、刚接触无功补偿的工程师也适合用 STATCOM 做毕业设计或项目预研的人。我会把整体思路、建模细节、参数计算和调参坑一次性聊透。1. 为什么选 STATCOM 做无功补偿原理与工程价值1.1 STATCOM 与传统补偿装置的本质区别传统无功补偿主要靠电容器、电抗器和 SVC。电容器只能发无功不能收无功想调压只能分组投切动作慢还容易产生涌流。SVC 用晶闸管控制电抗器和电容器组响应速度能到几十毫秒但它本质上是被动器件组合输出无功和母线电压强相关电压跌得越厉害SVC 能提供的无功反而越少这跟调压需求正好相反。STATCOM 不一样。它本质是一个受控电压源换流器通过 IGBT 桥把直流电容电压转换成与系统同步、幅值和相位都可控的三相交流电压再经过连接电抗器并网。改变输出电压和电网电压之间的幅值差、相位差就能控制流过电抗器的电流大小和方向。形象点说STATCOM 是一个“有源”的电流源能在系统电压跌的时候继续顶住无功输出短时过载能力也很强。这就是为什么在电压波动大、负载冲击强的场合越来越多工程直接上 STATCOM。1.2 无功补偿到底解决什么问题电网里大量感性负荷比如异步电机、变压器、电弧炉、光伏逆变器滤波器都会消耗无功。无功本身不做功但会占线路容量、增大电压降落、拉低功率因数。功率因数低了不仅供电部门考核设备容量利用率也下降。STATCOM 的核心任务就是快速检测系统当前的无功缺额输出对应的容性无功或感性无功把并网点功率因数拉回 0.95 以上同时稳定母线电压。这里有个容易误解的地方无功补偿不是把无功补到零才是最好。实际工程中常常按调度指令维持母线电压尤其弱电网末端需要系统性地补偿压降。所以 STATCOM 外环通常有两种模式无功功率模式给定 Qref电压模式用母线电压偏差生成无功指令。在 Simulink 建模时我建议把两种模式都做进去用 Switch 切换方便对比不同工况。后面我讲的参数整定也以这两种模式为基础。1.3 三相不平衡也是 STATCOM 的用武之地传统补偿装置大多只能分相投切但遇到单相焊机、电气化铁路牵引这类严重不平衡负载各相无功需求不一样光补三相总无功无法解决负序电流问题。STATCOM 的拓扑决定了它可以输出三相不等的无功电流只要控制策略配合比如采用分相控制或正负序分解控制就能把三相电压不平衡度压下来。不过负序补偿会让控制器复杂不少入门阶段还是先从三相平衡工况入手把 dq 坐标系下的控制跑通再考虑负序甚至谐波补偿不然很容易被一堆坐标变换绕晕。2. MATLAB/Simulink 仿真建模的整体思路2.1 建模选型Simscape Electrical 还是手写方程Simulink 里做电力电子仿真主要有三条路。第一条是用 Simscape Electrical 里的 Specialized Power Systems 库也就是以前 SimPowerSystems 那套里面有现成的三相电源、变压器、Universal Bridge、RLC 负荷适合快速搭主电路。第二条是 Simscape Electrical 的新架构 Foundation 库建模方式更现代但很多控制接口和老教程对不上排查问题反而麻烦。第三条是纯 Simulink 数学模块手写状态方程把主电路全部用微分方程表达优点是控制算法验证起来轻量缺点是器件细节被简化看不到开关纹波也无法验证 PWM 死区这些工程细节。我的建议很直接入门和项目复现优先选 Specialized Power Systems。原因就一句话元件成熟、资料多、模型稳。Universal Bridge 拖出来就能用可以选平均模型或详细开关模型控制部分能方便地拿到三相电流电压信号。新架构再美观也不如先把功能跑通重要。等你把控制思路彻底理解再迁移到其他建模环境也不迟。2.2 系统参数设定电网、容量、负荷场景建模第一步是定参数。我以一个常见的 10kV 配电系统为例STATCOM 容量设计为 ±10Mvar。先列一张表后面对着仿真调整参数数值说明系统线电压10kV有效值50Hz系统短路容量500MVA决定系统阻抗仿真中简化为串联阻抗连接电抗器3mH/相限制开关纹波配合解耦控制电抗器内阻0.1Ω模拟损耗不能设成 0直流侧电容20mF支撑直流电压动态过程用直流侧额定电压18kV两电平 SPWM 下需足够高开关频率2kHz常见两电平 IGBT 选择负载无功8Mvar 到 -2Mvar模拟感性负载变化连接电抗器不是随便选的。它一方面滤除高频开关纹波另一方面参与和电网的有功无功交换。电感太小电流纹波大、并网容易振荡电感太大无功输出范围受限制。工程上常按基波阻抗 5%~15% 的短路阻抗来取。10kV、10Mvar 系统额定电流约 577A单相电抗 3mH 在 50Hz 下阻抗约 0.942Ω基波压降约 0.942×577≈544V约为相电压 5774V 的 9.4%在合理范围。直流侧电压也要算。10kV 系统相电压峰值是 10k/√3×√2≈8.16kV。三相两电平 VSC 用 SPWM 时输出相电压基波幅值最大约 Vdc/2因此 Vdc 至少要到 16.32kV考虑电抗压降和调节裕量取 18kV。如果后续用 SVPWM直流电压利用率能提高约 15%可以再压低一些但入门阶段别省这个裕量。2.3 主电路拓扑两电平 VSC 加连接电抗器加直流电容主电路拓扑就用经典三相两电平电压源换流器直流侧一个大电容三相桥臂各两个 IGBT 带反并联二极管交流侧通过连接电抗器并入 10kV 母线。为观测效果在公共连接点装电压电流测量模块。Simulink 中的搭建方式是从 Simscape Electrical Specialized Power Systems 库拖出 Universal Bridge选择 Three-arm bridge 和 IGBT/Diodes电源用 Three-Phase Source 或 Three-Phase Programmable Voltage Source后者方便模拟电压跌落负荷用 Three-Phase Parallel RLC Load无功功率设为感性正值。注意三相三线制系统不一定要中性点但电力电子仿真通常需要给求解器一个参考地。我会把电源中性点设为 Yg 接地或者在测量点附近加一个接地阻抗避免出现矩阵奇异。入门阶段不要为了追求“不接地”的完美工程模型把自己卡在数值报错里。3. 控制策略与关键模块实现3.1 dq 坐标变换和瞬时无功理论STATCOM 控制的基础是把三相电流电压从 abc 静止坐标变换到 dq 同步旋转坐标。为什么要旋转变换因为理想情况下三相对称交流量变换过去后就变成了直流量PI 控制器对直流量的跟踪正好没有稳态误差。等幅值 Park 变换的写法是ud 2/3 × [ua×cosθ ub×cos(θ-120°) uc×cos(θ120°)] uq 2/3 × [ua×sinθ ub×sin(θ-120°) uc×sin(θ120°)]相位 θ 来自锁相环 PLL锁定 A 相电网电压相位。Simulink 里有现成的 PLL 模块也可以自己搭一个三相锁相环。我建议先用内置模块等控制跑通再拆开看内部结构。在 dq 坐标系下如果让 d 轴定向到电网电压矢量方向uq0那么有功功率 P1.5×ud×id无功功率 Q-1.5×ud×iq。这意味着控制 iq 就能独立控制无功控制 id 就能控制有功。STATCOM 虽然主要发无功但开关损耗和电抗器损耗需要从电网吸收一点有功来支撑直流电容所以 id 环不能省它承担直流电压稳定任务。3.2 外环电压和无功环、内环电流环结构控制系统用经典电压外环电流内环这是基于 dq 解耦的串级控制。外环负责生成电流指令内环负责快速跟踪电流指令。外环这样接如果工作在无功模式Qref 与实测 Q 的偏差进 PI输出 iq_ref如果工作在电压模式Vref 与 PCC 电压实测值的偏差进 PI输出 Qref再经过一个无功偏差 PI 得到 iq_ref。直流电压环用 Vdc_ref 与实测 Vdc 的偏差进 PI输出 id_ref用来补充有功损耗。内环考虑耦合典型控制方程是Vd_ref -PI_id(id_ref - id) ωL×iq ud Vq_ref -PI_iq(iq_ref - iq) - ωL×id uq其中 ωL×iq 和 -ωL×id 是解耦前馈项ud 和 uq 是 PCC 电压的 dq 分量。很多初学者直接拿 PI 输出当调制波不做解耦结果动态过程中 id 和 iq 互相拉扯波形越调越乱。前馈补偿加上后两个电流环基本独立PI 参数才好整定。3.3 电流环和外环 PI 参数整定电流环带宽一般取开关频率的 1/10 到 1/5。开关频率 2kHz 时带宽取 200~400Hz。内环 PI 可以按典型Ⅰ型系统来初步设计Kp L/τKi Kp×R/L。其中 τ 是电流环等效时间常数取 1msL3mHR0.1Ω则 Kp≈3Ki≈100。这个估算只是起点Simulink 里可以直接拧 PID 参数观察效果。我习惯先把 Kp 设小一点跑通比如 Kp1、Ki50再加增益防止一开始就振荡。外环带宽必须低于内环一般取内环的 1/5 到 1/10也就是 20~80Hz。如果外环太快内环跟不上整个系统会像两个小孩抢方向盘来回甩。直流电压环同理带宽要低于电流环。整定顺序永远是从内到外先锁住电流环再调电压环最后调无功外环。3.4 PWM 调制与开关信号生成调制方式建议先用三相 SPWM思路最直观把三相调制波归一化到 ±1和 2kHz 三角载波比较生成各桥臂上管驱动信号下管取反加死区。Simulink 里可以直接用 Relational Operator也可以使用 Universal Bridge 自带的 PWM Generator。仿真分两种模型平均模型和详细开关模型。平均模型用受控电压源替代开关桥仿真速度快适合验证控制逻辑但看不到开关纹波。详细开关模型用真实 IGBT 和 PWM 脉冲能看 THD 和死区影响但速度慢。我的习惯是先平均模型跑通所有工况确认控制没毛病再切成详细模型分析谐波。PWM Generator 选 Carrier-based PWM调制方式选 SPWM频率 2kHz。死区时间设成 3~5us避免桥臂直通。纯仿真阶段死区对波形影响不大但如果你把这个模型当作工程样机的前期验证死区必须写上。3.5 直流电压稳定和启动冲击抑制STATCOM 上电时直流电容是空的直接给控制指令会导致电容充电电流巨大仿真里表现为直流电压过冲、电流爆表实际工程中可能直接烧 IGBT。仿真时我用两段式启动先不投入 PWM用一个很低电压的直流源给电容预充电到 90%再切到正常 PWM 控制同时把 id 环输出限幅设成安全范围。如果不想额外搭预充电回路可以把直流电压参考值设成斜坡上升比如 0.5s 内从 0 升到 18kV同时把 PI 输出限幅到 ±2kA。这个斜坡时间不能太短否则直流电容充电电流仍然可能超过额定值。电容越大需要的斜坡时间越长但动态波动越小这是工程折中。直流电容容量粗略可按动态能量需求估算如果允许直流电压波动 5%某工况下需要吸收的能量变化 ΔE 除以允许电压变化 0.9kV再除以 2×Vdc可以反推电容大小。比如 ΔE 按 10Mvar 系统一个工频周期的能量波动粗估算下来约 1mF 起步。我选了 20mF 是很保守的值动态响应好但成本高。实际仿真时可以从 10mF 开始看波形不够再加。4. 完整仿真搭建实操流程4.1 新建模型和参数脚本打开 MATLAB新建一个 Simulink 模型再建一个 MATLAB 脚本statcom_params.m把所有参数写在里面方便统一管理和改参。下面是脚本的核心内容%% STATCOM 系统参数 Vll 10e3; % 线电压有效值 V freq 50; % 电网频率 Hz omega 2*pi*freq; % 角频率 rad/s %% 主电路参数 Lc 3e-3; % 连接电抗 mH/相 Rc 0.1; % 电抗器内阻 Ohm Cdc 20e-3; % 直流电容 F Vdc_ref 18e3; % 直流电压参考 V %% STATCOM 额定值 S_rated 10e6; % 额定容量 VA V_phase_peak Vll/sqrt(3)*sqrt(2); % 相电压峰值 %% 控制参数 fsw 2000; % 开关频率 Hz Ts 1e-5; % 仿真采样时间 s Kp_i 3; % 电流环比例 Ki_i 100; % 电流环积分 Kp_v 0.5; % 电压环比例 Ki_v 20; % 电压环积分在模型里所有模块的参数直接填变量名比如直流电容写成Cdc。这么做的好处是改一次脚本整个模型更新不用逐个点模块对话框。4.2 搭建主电路并添加测量从 Library Browser 调出以下模块Three-Phase Source线电压 10kVYg 接地Three-Phase Series RLC Branch作为系统短路阻抗阻抗值按短路容量 500MVA 折算Three-Phase V-I Measurement放在 PCC 点Three-Phase Parallel RLC Load设置有功和无功初始值Universal Bridge选 Three-armIGBT/Diodes直流电容和接地信号顺序是电源出来经过系统阻抗到 PCC 点一部分接负载一部分经连接电抗器进 Universal Bridge 的交流侧。Universal Bridge 的直流正负端接电容电容负极接地。注意 V-I Measurement 要同时输出三相电压和电流用于控制。4.3 搭建控制回路并处理信号控制回路可以分成三块取 PCC 三相电压进 PLL得到 θ。用 abc 到 dq0 变换模块把三相电压和三相电流都变换到 dq。Simulink 里用abc_to_dq0 Transformation要输入 θ。按照前文的外环内环结构搭 PI 调节器、限幅、前馈解耦和逆变换生成三相调制波送入 PWM Generator。这里强烈建议用标幺值。把电压除以额定相电压峰值电流除以额定电流峰值信号都落在 1 附近PI 参数好调得多。直接用 18kV、577A 这种大数值也能调但积分增益和比例增益量级差很多稍微碰一下就发散。控制环里所有信号都走标幺值最后在输出侧还原成实际值。4.4 配置仿真参数并运行仿真器配置我放在 Powergui 模块里。拖一个 Powergui 进模型选择离散仿真。采样时间 Ts 设 1e-5 秒也就是一个 2kHz 开关周期分成 50 步足够看清 PWM 细节。求解器用离散模式仿真时间先设 0.5s观察启动和阶跃过程足够。如果跑详细开关模型太慢可以先在 Universal Bridge 里选择平均值模型也就是 Average-Model Based VSC再用详细模型做一次最终验证。平均模型仿真速度能快十几倍适合反复调 PI 参数。5. 仿真结果分析与波形解读5.1 无功功率跟随效果搭好模型后第一件事看无功跟随。用 Scope 观察 Q 和 Qref。典型工况是0~0.2s 负载为 8Mvar 感性STATCOM 投入后 PCC 无功接近 00.2s 时 Qref 从 0 跳到 5Mvar观察响应时间。正常情况下 Q 应该在一个到两个工频周期内跟上超调小于 10%。如果响应太慢增大外环无功环 PI如果超调大甚至振荡检查电流环带宽是否足够、解耦项是否加对、PI 限幅是不是太宽。调参一定要一次只动一个环节否则波形乱了根本分不清是哪个环引起的。5.2 母线电压支撑与动态响应这是 STATCOM 的核心价值测试。把外环切到电压模式Vref 设为 1.0pu。在 0.3s 时用 Programmable Source 把系统电压跌到 0.85pu持续 0.2s再看 PCC 电压波形。理想情况下 STATCOM 会迅速输出无功把 PCC 电压拉回到接近 0.95pu。不过要清醒认识到 STATCOM 不是万能的。电压提升量和系统短路容量直接相关ΔV ≈ Q / S_SC。以 500MVA 短路容量、STATCOM 输出 10Mvar 计算电压提升只有 2%。要把 0.15pu 的跌落都补回来需要很大无功容量或在更弱的电网里。仿真的意义就在于让你直观看到补偿能力的边界。5.3 并网电流 THD 与滤波效果详细开关模型跑完用 Powergui 的 FFT Analysis 查看 STATCOM 并网电流 THD。两电平 2kHz SPWM 的 THD 可能到 10% 左右电流波形有明显的高频锯齿。如果想降到 5% 以下可以在交流侧加一个小滤波电容组成 LC 支路或者提高开关频率但这都会增加成本和损耗。实际工程里两电平 STATCOM 很少单机做很大容量很多采用级联 H 桥或多电平拓扑谐波性能天然更好。入门阶段不用纠结在滤波上先把控制策略跑通理解 THD 来源后面再升级拓扑。6. 常见问题与排查技巧实录6.1 启动时直流电压跌落或过冲现象仿真一启动直流电压从 18kV 掉到 10kV之后慢慢回升甚至冲顶。常见原因是直流电容初始电压没设置。双击电容模块把 Initial voltage 设成 Vdc_ref或者用状态初始化器设定。另一种原因是启动时 id 环饱和没及时给电容充能。把 Vdc_ref 设成斜坡上升同时给 id_ref 输出限幅就能把冲击压下来。6.2 电流环发散或振荡排查顺序很固定。先看 PLL 是否对准观察 uq 是否接近 0如果 uq 恒定不为 0说明 dq 定向有偏差需要补偿角度。再看电感单位3mH 在 Simulink 里就是 0.003如果误填成 3解耦系数差出千倍系统必炸。再有调制波限幅很重要SPWM 的调制波必须限幅在 ±0.95 以下超过 1 进入过调制区电流波形会明显畸变。6.3 PLL 失锁与相序错误PLL 输出频率不是 314或者 dq 分量剧烈波动先别怀疑 PLL 模型检查三相信号相序是不是 A-B-C。V-I Measurement 出来的信号顺序可能是 A、B、C也可能按母线排布不同Bus Selector 里选错一路就会导致相位乱跳。建议用 Scope 先看三相电压波形确认顺序再进 PLL。6.4 仿真速度慢和不收敛电力电子仿真跑不动多半是步长太小或开关事件太频繁。先确认 Powergui 是离散模式采样时间不要小于 1e-5。再确认 Universal Bridge 的 Snubber 参数保持默认不要设成 0否则数值求解可能出现奇异。最后Scope 不要挂太多把波形用 To Workspace 记录到工作区仿真结束后再统一画图速度提升非常明显。这里整理一张速查表方便你对照现象可能原因处理办法直流电压塌掉电容初始电压没设电容初始电压设为 Vdc_ref电流发散PLL 定向不准或解耦项错误检查 uq 是否≈0核对电感单位调制波饱和限幅缺失或 PI 输出过大调制波限幅 ±0.95降低外环增益仿真太慢步长过小、平均模型未用离散采样 1e-5先用平均模型波形有高频锯齿两电平开关纹波加 LC 滤波提高开关频率最后再分享两个我自己的实操心得。第一一定要先跑平均模型再跑详细开关模型先标幺化再调 PI不然主电路数值问题和控制问题叠在一起根本无从排查。第二所有参数都放进 MATLAB 脚本里统一管理每次仿真前先运行一次脚本。我因为改完脚本忘了运行用旧电感和新 PI 调了半天波形越调越怪最后才发现是参数没刷新。STATCOM 这套系统看着复杂但只要把 dq 变换、功率平衡、串级 PI 这条主线捋清楚后面再做 SVG、储能 PCS 甚至柔性直流换流器思路都是相通的。