Simulink中含DFIG风机的IEEE39节点系统暂态仿真建模

IEEE 39节点系统在电力系统仿真圈子里有多经典，不用我多说了。它就是New England电网的简化模型，10台发电机、39个母线节点，涵盖了火电、水电、负荷、励磁系统和调速器，几乎是每个做电力系统稳定性和暂态研究的同学都绕不开的标准测试系统。不过网上能下载到的IEEE39模型，绝大多数是纯同步机版本，传动系统、PSS、励磁那一套齐全，但你要想研究风电接入后对系统频率、暂态稳定、短路电流的影响，直接把风机模块塞进这个模型里，中间会遇到一堆预料不到的坑。

我最近做的这个项目，就是在IEEE39节点标准模型基础上，加入了带完整控制器的风机模块，用Simulink把整个系统搭通、跑出暂态结果。这个模型加完后，可以拿来做三相对称短路、风机脱网、负荷突变、风电渗透率变化对系统频率和功角稳定性的影响研究。这篇就把我从零搭建到调通的全过程记下来，包括风机选型、并网点选择、潮流初始化、故障仿真和常见的报错处理。适合正在做电力系统新能源接入研究、或者被导师安排去搭含风电的IEEE39模型的研究生参考。

1. 模型到底是什么，以及为什么值得做

1.1 一句话说清这个模型的组成

IEEE39节点系统，又叫New England测试系统，是1979年就提出来的经典算例。它由10台发电机、39个节点、46条支路组成，其中39号节点是一台等值的聚合机组，用来模拟外部大电网，其余9台机分别代表New England地区的不同电厂。系统基准频率60Hz，基准容量100MVA。这个系统最大的特点就是：结构不算特别复杂，但足够模拟出机电暂态现象，比如低频振荡、暂态失稳、功角摆动，所以几十年来一直是学术论文里的常客。

我做的版本，是在这套原始模型基础上，选择性地在某个负荷母线附近增加了一个风电场模块。风机用的是双馈感应发电机（DFIG）模型，包含风速输入、风力机、传动轴、感应发电机、转子侧变流器、网侧变流器、撬棒保护和桨距角控制。整个风电场先通过一台690V/34.5kV的箱变升压，再经线路接入39节点系统的某个母线。这样既能保留原系统的主体特性，又能真实模拟风电场并网后的动态过程。

1.2 为什么说带风机模块的版本很稀缺

你可能觉得，IEEE39模型网上一搜一大把，加个风机不是顺手的事吗？还真不是。原版IEEE39模型里的数据，不管是阻抗、导纳，还是发电机参数，都是标幺值形式给出的。而Simulink的Simscape Electrical元件库要求的是实际值和si单位。这就导致你需要把整张线路表手工转换，一个电容标幺值算错，后面潮流就全偏了。更麻烦的是，原系统里根本没有预留风电接口，你在19号节点并一套风机上去，系统的潮流分布会立刻变化，发电机出力、平衡节点功率全都要重新调。

另外，很多流传的所谓“含风机IEEE39模型”，其实是在某个节点上挂了个等效的电压源或者一个PQ负载，并没有真正建模风力机的气动特性、变流器控制和故障穿越逻辑。这种简化模型做潮流计算问题不大，但你要做暂态稳定分析、看故障后风机的有功恢复速度、研究低电压穿越对系统的影响，它就不够用了。所以我这次做的是带完整DFIG控制器的版本，虽然搭建过程麻烦一点，但出来的结果能发论文、能支撑课题研究，这价值就完全不一样。

1.3 这个模型适合谁用

如果你是电力系统方向的研究生，正在做风电并网稳定性、频率响应或者新能源接入比例相关课题，这个模型可以直接拿来做对比实验。如果你的方向是风机控制，你也可以把DFIG的转子侧变流器控制算法拆出来单独调参，再放回全系统验证。对于做电网规划的工程师来说，这个模型还能用来做风电选址的初步筛选——在不同母线上接入风机，看系统的小干扰稳定裕度和短路容量变化趋势。

2. 整体设计思路：风机往哪加、用什么风机、怎么接

2.1 风机类型选择：DFIG还是PMSG

接入模型前，第一个要拍板的问题就是选什么风机。现在主流的风机类型就两种：双馈感应发电机（DFIG）和永磁同步发电机（PMSG）。对这个项目来说，我最终选了DFIG，原因有三个。

第一，DFIG在学术研究里的地位太高了。大量风电并网文献用的都是双馈风机，你要是用PMSG，后面和别人的研究结果做对比时会多很多需要解释的变量。第二，DFIG的变流器容量小，只有转子侧30%左右，这意味着它的建模仿真速度比PMSG快不少。PMSG是全功率变流器，开关器件多，仿真步长限制更严格，放在39节点这样的大系统里，一次暂态仿真可能要跑几个小时，太拖进度了。第三，DFIG的动态特性比PMSG复杂，它和电网之间有直接的电磁耦合路径，故障的时候转子侧会产生很大的暂态电流，需要撬棒电阻保护，这个过程本身就是一个很有意思的研究点。

当然，这并不意味着PMSG不值得做。如果你研究的重点是低电压穿越和故障隔离策略，PMSG的全功率变流器架构反而更适合，因为直流母线把电网和发电机完全隔开了。我个人的建议是：先用DFIG跑通整个系统和仿真流程，后面如果想研究PMSG，直接把风机模块替换掉，电网部分不需要任何改动。

2.2 并网节点的选择逻辑

风机并网点的选择是整个项目里最需要动脑子的一件事。很多人上来就把风机接到29号节点，理由是这个节点负荷比较重、靠近系统中心区域。但实际做下来，我觉得要考虑三个因素。

一是节点的短路容量。风电接入后，系统短路容量会被稀释，如果选在原本比较弱的节点，风机一接入就可能进入电压稳定问题区。二是该节点的有功和无功潮流裕度。你可以在做潮流计算前，先跑一遍基础模型的load flow，看哪些节点的无功接近上限，这种节点尽量避开。三是和研究目标的匹配度。如果你要研究风电对区域间振荡的影响，那应该把风机接在能看到明显模态变化的位置，而不是随便找个负荷点。

我最后选在了18号节点附近。原因是这个区域原本负荷集中，线路阻抗比较大，电压支撑能力相对弱，风机接入后动态特性更容易观察出来。同时这个母线离39号平衡节点有一段电气距离，故障后风电与同步机之间的功角摆动非常明显，很适合做暂态对比分析。

2.3 风机接入后的系统结构和控制策略

风机不是直接挂上去就完了。DFIG系统的结构是这样的：风力机把风的动能转换成机械能，通过齿轮箱和传动轴带动感应发电机的转子。转子侧变流器（RSC）负责控制电磁转矩和无功功率，网侧变流器（GSC）维持直流母线电压稳定，同时可以在电网侧提供无功支持。桨距角控制器在大风速时调节叶片角度，限制转子转速不要超速。

把这个结构映射到Simulink里，大概就是6个核心模块：风速模块、风力机模块、传动系统、发电机模块、RSC控制、GSC控制。RSC控制是最核心的，它采用定子磁链定向矢量控制，内部是dq旋转坐标系下的双闭环——外环控制有功功率和无功功率，内环控制转子电流。GSC则是电网电压定向，外环盯直流母线电压，内环盯网侧电流。这套控制器的参数，也就是PI增益和积分限幅，很大程度上决定了风机故障后的动态响应质量。

我具体在建模时，参考了MATLAB官方文档中DFIG平均模型的思路，但保留了撬棒保护机制。因为IEEE39是一个多机系统，单台风机和多机之间的交互本身就复杂，如果再把撬棒逻辑去掉，故障后的暂态曲线会明显失真的。

3. 实操记录：从参数整理到仿真结果完整流程

3.1 准备IEEE39基础模型

这个项目的第一步，是最基础但也是最枯燥的：把IEEE39的参数吃透，整整齐齐地铺到Simulink里去。我用的数据是标准算例里的那份，系统包含39个节点、46条支路、10台机组和19个负荷。线路参数表、变压器变比表、发电机动态参数表，三张表必须一一对应着来。

我在搭基础模型时踩了个大坑：一开始图省事，直接引用别人转好的Simulink模型，结果发现里面的负荷模型全是恒阻抗的。这意味着电网电压一降，负荷功率也跟着降，系统实际阻尼被高估了，后面做短路分析的时候，曲线特别“完美”，反而失真。后来我自己重新搭了一遍，把负荷改成恒功率模型，也就是用三相动态负载模块，把有功和无功设成常数，这样电压下降时负荷功率基本维持，更接近真实电网特性。这里强烈建议大家，尽量自己动手搭无源网络部分，别偷懒。

搭建无源网络的具体流程是：先用Three-Phase V-I Measurement模块连接每个母线，母线之间按线路阻抗参数串RLC支路。这里需要把标幺值转换成实际值，比如一条138kV输电线路的标幺阻抗是0.001pu（基准100MVA），换算成实际欧姆就是0.001×138²/100 ≈ 0.19Ω。变压器部分用Three-Phase Transformer（Two Windings），注意设置饱和特性为不启用，因为机电暂态仿真里磁饱和不是重点。发电机部分用同步发电机模块，配上IEEE Type 1励磁系统、调速器，再单独加一个Power System Stabilizer（PSS）模块，这样模型才具备完整的动态响应能力。

3.2 将风机模块接入系统

风机模块的位置，我放在了18号母线附近，通过一条短线路接入。风力机参数我用的是Simulink里面Simscape Electrical的Wind Turbine模块。这个模块需要给定一个基准风速，以及基本的空气动力学参数。DFIG发电机的转子和定子参数，可以从典型1.5MW风机数据里找到，我直接使用了软件自带的样例参数。这个样例参数虽然是一个具体的风机型号，但因为它是标幺值建模的，放在100MVA基准的39节点系统里，相当于一台容量较小的机组。可别小看这一步，单纯“1.5MW”这个值很容易让人忽略系统的基准容量，直接混进去结果全错。

为了让风电场的容量对系统有可见的影响，我把单台DFIG基准容量调到了50MVA，也就是相当于一个50MW的风电场，通过一台升压变压器接入系统。你可能觉得50MW在39节点系统里不算大，实际上IEEE39的总负荷约6000MW，50MW占比不到1%，动态波形上影响确实有限。所以我又加了一台同样的DFIG并联，组成一个100MW的风电场。100MW大约占系统总容量的1.6%，单机故障扰动下已经能在角速度曲线上看到明显影响，又不会让系统直接失去稳定，作为研究风电平稳接入影响的比例刚刚好。

3.3 潮流初始化与关键仿真参数设置

这一节是整个项目里最枯燥但也是决定性的一步。Simulink里的机电暂态模型要求初值必须满足潮流方程，否则仿真开始就是“炸”的。我的做法是借助Powergui里面的Load Flow工具来做潮流初始化。

Load Flow工具里，需要对每个母线设置节点类型：平衡节点（Node 39）、PQ节点、PV节点。基准电压设置成实际值，比如18号母线是230kV。同步发电机要设定出力Pm和端电压Vt，负荷节点的有功无功直接设成原始数据里的标幺值乘上100MVA。风电场在这个阶段要设置成PQ节点，有功出力和无功出力按目标值设定，我这里设的是100MW、0MVar。

潮流计算跑完以后，重点看两个指标：一是各母线电压是否落在0.95到1.05pu之间，二是平衡节点39号机组的出力是否合理。我第一次跑的时候，34号母线的电压跑到1.12pu，直接超出限制。查了半天，发现是变压器的分接头位置没设对。IEEE39原始数据里的变压器变比不是标准的1:1，需要按照输电数据表的抽头位置设置。把变比校正后，电压就都恢复正常了。

接下来设置仿真参数。Simulink里仿真模式选择“Continuous”还是“Discrete”，我建议视风机模型而定。如果你用的是详细开关模型，建议用离散模式，步长设置在50微秒左右，因为PWM开关频率一般很高，你必须要捕捉到开关动作。如果用的是平均模型，可以放宽到100微秒，求解器选择ode23tb，相对容差设为1e-3。这个设置比系统默认值靠谱得多，否则容易因为刚性方程问题导致“数值不稳定”报错。

3.4 实际仿真场景设计：短路故障与风机动态响应

模型搭建完成、初始化通过之后，就可以开始做正事了。我最常用的一个场景是三相对称短路故障。故障点设在18号母线附近，也就是风电场并网点不远的线路上。这样做的目的是：让风机近距离体验一次电网电压骤降，观察它在低电压期间的表现。

故障时序大概是这样：0秒开始运行，系统稳定在额定工况。1.0秒时故障发生，用Three-Phase Fault模块模拟三相对地短路，故障持续时间设为0.1秒，1.1秒故障切除。仿真总时长设为5秒，足够捕捉故障后系统的暂态恢复全过程。

仿真跑完后，我最关注的曲线有四条：风电场出口电压曲线、有功功率输出曲线、风机转子转速曲线，以及同步发电机的功角曲线。电压曲线能看到故障期间电压跌到约0.3pu，故障切除后电压能否快速恢复到0.95pu以上；有功功率曲线能看出风机在故障期间是否按控制策略降低了输出、切除后多久恢复；转子转速曲线能看出传动轴是否出现过度扭振；功角曲线则反映整条系统在风电接入后的暂态稳定裕度变化。

实测下来，加入风机模块后最明显的变化是故障后的频率最低点比纯同步机系统稍微低了一些。这说明风电替代了部分同步机容量之后，系统的等效转动惯量下降了。如果你研究的是高渗透率风电系统，这个现象就是最直接的观察对象。

4. 我踩过的坑：问题排查与细节建议

4.1 常见问题速查表

我把搭建过程中遇到的各种报错和现象整理成了一张速查表，方便你对照排查：

4.2 数值稳定性与仿真速度的平衡

这是Simulink多机系统建模最折磨人的一点。IEEE39本身是一个刚性系统，发电机的转子运动方程时间常数在秒级，定子电磁暂态在毫秒级，而风机变流器的PWM开关又在微秒级。三种时间尺度叠加在一起，对求解器是个非常残酷的考验。

我的做法是“三层隔离”：先用平均模型的风机跑通整个系统的控制逻辑，也就是把变流器等效成受控电压源；当所有波形都合理、初始化也稳定后，再针对某个特定场景切换成详细开关模型，跑一次精确结果。这样做的好处是，大部分流程性工作用平均模型完成，省80%以上的仿真时间，最后再用详细模型出精确数据。很多初学者一上来就要求“模型越精细越好”，结果搞了三天连初始化都过不去，就是没掌握这个分层的思路。

另外，求解器的相对容差不要设得太紧，1e-3就行。我试过1e-4，误差是小了，但仿真时间翻了四五倍，而波形上的差异肉眼基本看不出来。做工程仿真，要清楚地知道精确度够用就好，系统参数的精确度误差早就超过那一点数值容差了。

4.3 关于“仿真结果可信吗”这件事

模型跑出来的结果，到底能不能信，这个问题值得认真说。我有三个验证手段。第一，把风机模块的出力设成0，也就是风电场不发有功，此时系统的动态响应应该和纯同步机IEEE39基本一致。如果对不上，说明风机模块在零出力状态下还和系统存在异常交互，这时不要急着往下做。第二，对比风机并网后的稳态潮流，电压幅值和相角的分布要和原始模型相差不大，如果某些母线电压明显偏离，那说明风机接入点选择不当或等值容量过大。第三，做一次小扰动验证，比如在1.5秒时给某台同步机增加5%的机械功率，观察系统能否恢复稳定，以及振荡频率是否在0.3到2.0Hz的典型区间内。

这三步验证做完，模型的可靠性才算真正立住了。此后得到的故障暂态曲线、频率动态曲线、功角摇摆曲线，才敢写进论文或者支撑决策。

5. 给后续研究者的扩展思路

模型搭好以后，后面能做的事情其实非常多。我目前自己已经在尝试的一个方向是虚拟惯量控制。因为前面说了，风机接入后系统等效惯量下降，频率跌落变快，这时候可以通过修改DFIG的转子侧变流器外环控制，在频率变化过快时释放或吸收转子动能，模拟同步机的惯量响应。你可以用这个模型先做一组对比：不加虚拟惯量，频率最低点是多少；加入虚拟惯量后，频率最低点能提升多少，要不要牺牲一部分机械疲劳寿命来做交换。

还有一个方向是低电压穿越策略研究。故障期间DFIG转子电流会飙升，如果不投入撬棒电路，变流器可能过流损坏。你可以在这个包含详细电网结构的系统里，对比不同撬棒阻值、不同投入时序下的转子电流峰值和系统恢复速度，这类结论对风电场的实际运行很有参考价值。再多说一句，你还可以把风电渗透率做成一个参数，从5%到30%逐个试，看系统临界切除时间怎么变化，这是分析新能源极限渗透率的一种经典做法，很多论文的精华就在这个扫描结果里。

我自己做这个项目时最深的体会是：仿真模型的精度并不是越高越好，核心是你要回答什么问题。如果你只是想看风电场对系统频率的影响，平均模型完全够用。非要用详细开关模型，那是自找麻烦。如果你是做风机变流器硬件控制的，那又另当别论了。在这个基础上，还有一个很实用的建议：小容量、单机无穷大系统永远是最好的调试场地。先把风机控制调明白，再往IEEE39全系统里移植，别一上来就用全系统挑战自己。这样一步步来，你会少踩很多坑，也更能在论文里讲清楚每个参数选择的来龙去脉。