Simulink双馈风力发电机并网故障仿真与LVRT分析

我自己做双馈风力发电机并网故障仿真也踩了不少坑，从最开始拿Simulink里自带的Demo模型改参数，到后来自己搭完整的DFIG并网系统，中间折腾了很久。这篇就把整个研究过程中最核心的东西整理出来，包括建模思路、故障注入方法、波形分析要点，以及我踩过的各种坑，希望对做同类课题的人有帮助。

1. 项目背景与核心需求解析

1.1 为什么选“并网故障”这个切入点

双馈风力发电机（DFIG，Doubly Fed Induction Generator）是目前陆上风电场装机量最大的机型，核心特点是转子侧通过背靠背变流器接入电网，而定子直接并网。这个拓扑结构决定了它对电网故障非常敏感——电网电压骤降时，定子磁链不能突变，会产生直流分量和负序分量，进而在转子回路感应出很大的电动势，导致转子过流、直流母线过压。搞不定这个问题，风机就得脱网，脱网又反过来加剧电网故障，这就是前些年风电场大规模脱网事故的根源。

所以“并网故障研究”这个题目，本质上研究的是三件事：一是DFIG在电网故障下的电磁暂态特性，二是故障对变流器和直流母线的影响机理，三是故障穿越策略的验证。用Simulink做这件事，是因为它把风力机、发电机、变流器、电网和控制系统都放在同一个仿真环境里，能完整看到机电暂态和电磁暂态的耦合过程，这个优势是单纯用Matlab脚本算解析公式比不了的。

1.2 研究内容边界与目标设定

我给自己定的研究范围是这样的：只研究对称和不对称故障下的并网特性，故障点设在PCC（公共连接点），故障类型覆盖三相对称短路、两相短路和单相接地短路，故障持续时间按并网导则的低电压穿越要求设置在0.625秒到3秒之间可调，重点是看故障期间和故障切除后的动态过程。

目标不是把控制策略做得多花哨，而是先把故障特性吃透：记录转子电流峰值、直流母线电压波动范围、有功无功功率的暂态过程、电磁转矩振荡情况，用这些数据去验证理论分析的结论，比如为什么对称故障下转子过流最严重，为什么单相故障下直流母线二次谐波会变大，这些有了仿真结果支撑，后面写论文或者做工程方案才有说服力。

2. Simulink中DFIG并网系统的建模思路

2.1 整体拓扑结构与建模路径选择

DFIG并网系统在Simulink里有三种建模路径，这个选择直接影响后面所有工作。第一种是用Specialized Power Systems库里的现成DFIG模型，在Simulink的powerlib里直接拖出来用，优点是快，缺点是内部结构不透明，很多参数被封装死了，做故障分析时想看转子磁链的细节就抓瞎。第二种是自己从零搭建，包括异步电机模型、两电平变流器、PWM调制、矢量控制系统全部手写，工作量大，但对机理的理解最透彻。第三种是混合方式——用库里的DFIG电磁模型，但自己搭控制部分和故障注入逻辑，这是我在实际操作中最终选定的方案。

自己的体会是：除非是研究变流器拓扑本身，否则没必要从零写电机模型。Simulink自带的DFIG模型已经考虑了定转子绕组动态、磁链饱和效应（在饱和选项里启用后），做故障暂态分析精度完全够。重点精力应该放在控制策略的修改和故障场景的搭建上，这部分决定了研究成果的差异化。

整个系统的拓扑结构包括以下部分，每个部分在仿真模型里都是一个独立的子系统：

• 风力机模型：输入风速，输出机械转矩，需设置桨距角控制参数
• 齿轮箱模型：通常用刚性轴简化或用两质量块模型模拟轴系扭振
• DFIG电机本体：定子直接并网，转子经滑环接入变流器
• 转子侧变流器（RSC）：控制DFIG的有功无功解耦
• 网侧变流器（GSC）：维持直流母线电压稳定
• 直流母线电容：滤除开关纹波，稳定母线电压
• 网侧滤波电感：抑制PWM谐波，连接GSC与电网
• 撬棒保护电路（Crowbar）：故障时短接转子回路，保护RSC
• 并网开关与故障注入模块：控制并网与故障时序

2.2 控制策略：转子侧与网侧变流器怎么配合

转子侧变流器用定子磁链定向矢量控制，这个选择不是随便做的。定子磁链定向的好处是，稳态时定子磁链矢量基本恒定，转子电流的d轴分量控制有功，q轴分量控制无功，解耦关系简单清晰。故障发生时定子磁链会剧烈变化，定向会失效，但这正是需要观察的暂态过程。

转子侧变流器的控制结构是典型的双闭环：外环是功率环，有功功率参考值由最大功率跟踪（MPPT）给定，无功功率参考值通常设置为0或按电网调度要求给定；内环是转子电流环，用PI调节器跟踪外环给出的电流参考值。电流环的带宽要远高于功率环，一般取电流环带宽200-500rad/s，功率环带宽10-30rad/s，这个比例关系我调PI参数时反复验证过。

网侧变流器相对简单，用电网电压定向矢量控制，d轴控制直流母线电压，q轴控制网侧无功，正常情况下维持单位功率因数运行。直流母线电压环的响应速度要适中，太慢会导致故障时母线过压严重，太快又会让网侧电流出现较大冲击，这里有个系统性的折中。

需要非常注意的是电流内环的解耦补偿项。在Simulink里实现时，转子电流d-q轴之间存在交叉耦合项，分别是转差频率与转子电感的乘积项。如果不做解耦补偿，电流环的动态性能会明显变差，故障时尤其明显，会看到d轴电流变化时q轴电流跟着抖动。我最初省略了这项，做三相对称故障仿真时转子电流波形出现了不该有的振荡，补上解耦项之后波形才合理。

2.3 参数设置的关键细节

DFIG的电机参数我直接用了一台2MW级风机的典型参数，定子额定电压690V，额定频率50Hz，定转子匝数比0.38左右，定子电阻0.0043pu、漏感0.13pu，转子电阻0.0055pu、漏感0.16pu，励磁电感3.2pu。这些参数在Simulink的预设DFIG模型里能找到接近的版本，关键是拿到参数后要先做空载和负载验证，确认模型能稳态运行再开始故障仿真。

直流母线额定电压我设置为1150V，这是基于690V交流侧相电压峰值再留出约15%的裕量算出来的。母线电容取0.05F，这是2MW功率等级下比较常见的取值。网侧滤波电抗器取0.15pu，这个值要综合考虑滤波效果和网侧变流器的电流控制带宽，太大会限制电流响应速度，太小则谐波大。

变流器开关频率设置到2kHz，实际工程中这个等级的变流器一般就是1.5-3kHz。开关频率决定了PWM谐波的主要频段，也直接关系到仿真速度和步长选择，后面讲到仿真性能时还会提到。

3. 并网故障仿真的实现步骤

3.1 故障场景设计与故障时序控制

故障仿真的第一步是设计故障场景。我在模型里用Three-Phase Fault模块做故障注入，放在升压变压器与无穷大电网之间，模拟输电线路出口处发生短路故障。

故障时序控制是这里面的关键技巧。我做了两个方案：一个是在Three-Phase Fault模块里直接设置切换时间（transition times），用时间数组控制合闸和分闸时刻，这种方法适合固定时序的批量仿真；另一个是用Signal Builder或Step模块驱动故障控制端，配合断路器模块实现更复杂的时序，比如故障期间改变故障阻抗模拟电弧特性。实际做研究时我主要用第一个方案，简单可靠、重复性好。

一组典型的故障仿真时序是这样的：0到0.6秒正常运行，风和变流器控制都已进入稳态；0.6秒时刻投入故障，设置故障持续时间100ms（对应并网导则对近端故障的要求）；0.7秒故障切除；继续仿真到1.5秒左右观察恢复过程。仿真总时长1.5秒到2秒就够，再长计算量增大且信息量不再增加。

3.2 故障类型与故障阻抗的参数化设计

为了让研究更有系统性，我把故障类型和故障阻抗设计成参数化，通过Matlab脚本批量修改并运行模型。这里给出我在simulink中设置故障模块参数时使用的核心脚本，以及关键参数含义：

% 故障参数定义 faultType = 1; % 1-三相短路，2-两相短路，3-单相接地 faultR = 0.01; % 故障电阻/pu，越小故障越严重 faultX = 0.1; % 故障电抗/pu，电网等效阻抗的一部分 % 修改Three-Phase Fault模块参数 blockPath = 'DFIG_Grid_Fault/Three-Phase Fault'; set_param([blockPath '/Fault Resistance'], 'Value', num2str(faultR)); set_param([blockPath '/Fault Reactance'], 'Value', num2str(faultX)); % 根据故障类型设置三相故障切换方式 % 1-三相，2-A相接地，3-AB相短路，4-AB相接地 set_param([blockPath '/Fault Type'], 'Value', num2str(faultType));

这里有个特别容易忽略的地方：Fault Resistance的单位不是实际欧姆，而是标幺值。Simulink里Three-Phase Fault模块的故障阻抗参数是按pu值设定的，基准阻抗等于电网额定相电压的平方除以额定容量。如果拿实际欧姆值直接填进去，故障严重程度会完全偏掉。我自己就犯过这个错误，本来想模拟一个5欧姆的短路电阻，填进去发现电压跌得几乎为零，因为换算成pu后阻抗太小了。正确的做法是先算出基准阻抗，再换算成pu值填进去。

3.3 故障严重程度的量化设置

研究并网故障不能只做“短路”一种工况，工程上更关心不同程度的电压跌落。我设计了四种典型工况：

• 轻度故障：电压跌落至额定值的80%，故障阻抗较大，用于观察控制系统的调节能力
• 中度故障：电压跌落至额定值的50%，这是最常见的低电压穿越考核点之一
• 严重故障：电压跌落至额定值的20%，接近变流器耐受极限
• 极端故障：机端残压接近于零（三相金属性短路），此时必须依赖撬棒保护

通过调整故障阻抗的数值来控制电压跌落深度。具体做法是：先大致估算故障阻抗，仿真后看PCC电压幅值，再反向微调阻抗来逼近目标跌落深度。这个迭代过程通常三四次就能收敛，比直接套理论公式更实用，因为线路阻抗、变压器漏抗都会影响实际跌落深度。

4. 故障波形分析：从仿真结果到物理规律

4.1 三相短路故障下的电磁暂态特征

三相对称短路是研究DFIG故障特性的基准工况，因为它的分析最清晰，也最能暴露DFIG的核心薄弱点。我做了电压跌落至20%的仿真，记录到的现象与理论预期高度吻合。

故障发生的瞬间，定子电流没有突变，因为定子磁链不能突变，但定子磁链中会感应出直流分量和负序分量。这个直流分量在转子回路中感应出很大的转差频率电动势，导致转子电流剧烈升高。仿真中看到转子电流峰值达到额定值的3.5倍左右，持续大约80ms后开始衰减，衰减时间常数由转子回路阻抗决定。

直流母线电压的变化也很有特征：故障瞬间先是快速上升，因为转子侧能量回灌到直流母线，然后被网侧变流器调节拉回来，但故障期间母线电压仍然存在明显的波动。这个波动幅度与故障严重程度正相关，也直接考验直流母线电容和网侧变流器的调节能力。我实测电压跌落至20%时，母线峰值达到额定值的1.3倍以上，这个数据用来校核母线电容选型和撬棒动作阈值非常有参考价值。

电磁转矩在故障期间出现大幅振荡，甚至出现负值——也就是电机短时间变成电动状态，这会给齿轮箱和轴系带来很大的机械冲击。两个质量块模型的仿真结果显示，轴系扭转角在故障切除后还有一个持续的振荡过程，频率大概在1-2Hz，这就是轴系扭振模式。

4.2 不对称故障与对称故障的差异化分析

不对称故障比对称故障更常见，工程上单相接地短路占了所有输电线路故障的八成以上。我在仿真中对单相接地和两相短路做了对比，发现不对称故障的核心特征是负序分量的出现。

定子电压中含有负序分量时，定子磁链也会产生负序旋转磁动势，它相对转子的转速是两倍同步速（因为转子本身已经在以接近同步速旋转，负序磁场则以同步速反转，相对转速接近两倍同步速），所以在转子回路感应出约100Hz的电动势。这个100Hz分量叠加在原本的转差频率分量上，导致转子电流出现明显的二次谐波畸变。

直流母线电压也会出现100Hz的脉动，因为转子侧瞬时功率中存在两倍频分量。如果用FFT分析母线电压波形，能看到100Hz分量占主导，这正是不对称故障的标志性特征。单相接地故障时直流母线电压的波动幅度甚至比两相短路更明显，这一点容易被忽略，因为我一开始以为两相短路更严重，直到看到仿真数据才发现单相故障时负序电压标幺值更高。

为了更直观地对比，我把三类故障下的关键指标整理成了表格：

表格里的数据是我针对典型工况仿真得到的，不同参数下数值会有差异，但数量级和相对关系有普遍参考价值。做深入研究时建议用FFT工具把转子电流和定子电流做频域分解，负序分量随时间的变化过程很值得单独分析。

这个对比说明一个关键问题：做故障穿越研究时，只做三相对称故障是不够的。实际电网故障大多是不对称的，而不对称故障下的控制策略需要额外处理负序分量，这和对称故障下只要应付正序电压跌落完全不同。

4.3 低电压穿越（LVRT）过程的分析方法

并网故障研究落到工程应用上，最重要的就是低电压穿越能力评估。按照风电场并网导则的要求，并网点电压跌落到20%额定电压时，风电机组必须不脱网连续运行625ms，之后电压恢复过程中也要能保持并网。

在Simulink里做LVRT仿真时，我通常按上述故障场景设置0.6秒开始故障，持续0.625秒，到1.225秒恢复，然后观察以下关键指标：故障期间是否触发变流器过流保护、直流母线电压是否超过保护阈值、故障切除后有功功率能否快速恢复到故障前水平、无功功率在故障期间是否按导则要求发出无功支撑电网电压。

故障期间发出无功支撑这一点，我在控制策略里做了特殊处理：检测到电压跌落后，立刻把转子侧无功电流参考值切换为电压支撑模式，按照“每跌落1%电压，发出2%-6%额定无功电流”的斜率给定。这是现代风电场并网导则的明确要求，仿真结果也显示这种策略能明显改善故障点电压恢复速度。

撬棒保护动作与恢复是LVRT仿真里最考验建模细节的部分。我采用的策略是：当转子电流超过1.5倍额定值或直流母线电压超过1.2倍额定值时，投入撬棒短接转子回路，同时封锁RSC脉冲；当转子电流回落到额定值以下且持续20ms后，切出撬棒并重新解锁RSC。这个滞环逻辑如果不做，撬棒会在故障期间反复投切，给机械开关带来巨大压力。

5. 常见问题与排查技巧实录

5.1 仿真发散：代数环与数值稳定性

DFIG并网模型最大的梦魇就是仿真发散，尤其是在故障瞬间。故障引起状态变量剧烈变化，如果模型中存在代数环或数值刚性，仿真就很容易在故障时刻报错并终止。

排查代数环的一个实用方法：在Simulink模型菜单中执行“Diagnostics > Solver > Algebraic Loop”检查，如果检测到代数环，系统会给出提示。我的模型中曾经在功率计算环节出现过代数环——因为我用了电压电流的瞬时值直接算功率，而功率值又被用去生成电流参考值，形成了瞬时反馈回路。解决方法是加一个内存块或者一阶惯性环节切断代数环，或者改用有功无功计算模块内部带滤波的方式。

代数环的处理原则：能避免就避免，不能避免就加小惯性。加惯性环节时要选时间常数远小于控制周期，否则会影响控制动态响应。

5.2 固定步长与可变步长的拉锯战

变步长求解器（如ode23tb）在正常工况下效率很高，但故障瞬间数值刚性突增，步长会被压缩到极小，甚至出现“步长太小而无法满足精度要求”的报错。固定步长求解器（如ode3、ode4）不存在这个问题，但为了让PWM调制和解算精度足够，步长必须足够小，通常要取开关周期的1/50到1/100，也就是2kHz开关频率下步长取2e-6到5e-6秒，这会带来巨量的计算负担。

我最终采用的是折中方案：正常稳态运行用变步长求解器，设置最大步长限制为1e-4秒，保证计算效率；故障时刻前后改用固定步长，步长设为5e-6秒，观察故障瞬间的快速暂态过程。这个切换不是自动的，需要在仿真前预估故障时段，手动把仿真参数改掉。虽然笨一点，但能兼顾效率和精度。

如果不想手动切换，可以在模型里把故障时刻附近的PWM发生器跳过，直接用等效平均模型代替变流器，这样步长可以放宽到10e-5秒量级，仿真速度快很多。代价是看不到开关纹波，但故障暂态的宏观趋势完全正确。

5.3 “坏变量存储”和MATLAB版本兼容问题

Simulink有时候会弹出“Unable to save坏变量”这类错误，尤其是仿真中途强制停止后再次运行时。这通常是MATLAB工作区里有旧版本的变量与模型引用的变量不匹配造成的。排查方法是依次执行clear all、close all、clc清理工作区，删掉模型缓存文件夹（slprj），重新加载模型。

还有一类问题是不同MATLAB版本对模型文件的兼容性。我用MATLAB R2021a建的模型，在R2019b上打开时部分Specialized Power Systems模块会报错。如果要在多个版本之间协作，建议保存模型时指定为较早的版本格式（如R2019b格式），或者导出为slx后统一在目标环境重新验证。

5.4 仿真结果“不像论文里那样平滑”的原因

很多初学者第一次做DFIG故障仿真，出来的波形毛刺很大，甚至看不出规律，就开始怀疑模型错了。其实这大概率是PWM开关纹波叠加在工频波形上导致的，是正常的。

但如果波形毛刺呈现出明显的锯齿状，且与开关频率不同步，我建议先检查以下项：PWM采样方式是否用了过采样而非自然采样；电流测量模块之前是否加了抗混叠滤波；求解器步长是否偏大导致PWM边沿采样失真。我遇到最多的是第三个原因，把步长从默认的自动改成固定小步长后，波形立刻干净了。

5.5 互感与变压器连接方式的细节坑

做并网故障仿真时变压器接线方式会影响零序分量的流通路径。DFIG并网系统中常用的升压变压器是Dyn11连接，原边三角形绕组为零序电流提供了流通回路，副边星形接地给电网故障电流提供了返回路径。如果变压器连接组别设置错误，单相接地故障的零序网络就完全不对，仿真出的故障电流和电压分布都会失真。

设置变压器模型时要注意输入的是连接组别编号（如Dyn11），Simulink里要正确选择绕组连接方式和接地方式。我见过不少人直接用理想变压器代替三绕组变压器模型，这样做稳态仿真没问题，但单相接地故障时因为缺少励磁支路和漏抗分布，仿真结果与实际情况偏差很大。推荐使用Simscape Electrical里的三绕组变压器模型，并严格按照实际铭牌参数填写。

6. 研究结果的工程落地价值与个人经验

这套Simulink双馈风力发电机并网故障模型，不仅能用在校验理论分析上，也能直接支撑工程层面的工作。风电场并网评估、变流器保护定值整定、低电压穿越改造方案验证，都能在这个模型基础上扩展。比如我自己就用这套模型给一个变流器厂家验证过撬棒电阻参数，他们把不同阻值带入我的模型，对比直流母线电压和转子电流峰值，几十组工况跑完，很快就圈定了阻值的合理区间，比用物理样机反复试验省钱省时得多。

再有就是模型的数据后处理，我习惯用MATLAB脚本把Simulink的输出数据自动批量处理，直接生成故障特性报告。核心逻辑很简单：把To Workspace模块记录的波形数据在仿真结束后自动绘制成标准化格式的曲线图，同时提取峰值、恢复时间、振荡频率等指标。这个流程做成脚本后，几十组故障工况的对比分析一晚上就能跑完，第二天起来直接看汇总表格。

个人经验里还有一点忠告：别迷信模型库里的默认参数。仿真之前花半天时间把所有电机的物理参数、控制器的PI参数、保护阈值通读一遍并做合理性分析，后面能少走很多弯路。有一次我做低电压穿越仿真，发现故障切除后有功功率恢复特别慢，排查了两天才发现是转子侧有功外环的PI参数是从别的模型复制过来的，带宽低了一个数量级，根本不是故障特性的问题。

故障注入深度和机组运行点（风速）之间的耦合关系也值得关注。风速不同时，DFIG在故障发生前的转速和转差率不一样，转差率又直接影响故障后转子暂态电动势的大小，所以同样的故障深度，在不同风速下观测到的冲击水平可能差别很大。做系统性研究时建议加入风况这个因素，结果会更完整。

后面准备在这个模型基础上继续扩展的方向，一个是把电网侧换成含多机系统的微电网架构，研究DFIG故障对相邻机组的影响；另一个是把故障类型从简单的短路拓展为电压骤升和频率扰动，看看DFIG在这些非典型故障下的行为。希望这篇内容能给做同类课题的朋友一些实质性帮助。