PSCAD下张北柔直四端MMC建模与故障仿真实战

做张北柔直工程仿真这件事我从搭第一个MMC桥臂到把四端模型完整跑通前后折腾了近两个月。那会儿周围不少人一听“±500kV、四座换流站、几千个子模块”第一反应都是这玩意在PSCAD里能跑得动等我真的建完才发现张北柔直的魅力根本不在把每个IGBT都搬进仿真而在于四端拓扑、MMC建模和多站控制这三件事需要你在同一个模型里协调起来。这篇文章就把我从零搭建、调通并用它做故障分析的全过程记录下来给正准备用PSCAD复现柔直工程、或者正在啃MMC多端直流仿真的朋友做个参照。1. 张北柔直工程的四端网架与仿真还原边界张北可再生能源柔性直流电网试验示范工程2020年正式投运是世界上第一个真正意义上的四端环形柔性直流电网。它由张北、康保、丰宁、北京延庆四座换流站构成直流电压等级±500kV总换流容量9000MW直流线路总长超过600公里。很多文献管它叫“四端配电网”严格按电压等级和应用场景来说这是输电网的范畴但四端结构在潮流控制、故障隔离等方面跟未来主动配电网的思路高度相似这个叫法也就被大家叫开了。1.1 四座换流站各自扮演什么角色我刚开始建模时最头痛的不是MMC模型本身而是搞不清四座站到底谁该定直流电压、谁该送功率。后来把每座站的角色拆开看思路一下就顺了。换流站方向定位额定容量接入对象张北站送端3000MW张北地区风电场、光伏电站康保站送端1500MW康保地区风电场、光伏丰宁站功率支撑1500MW丰宁抽水蓄能电站北京站受端3000MW北京电网负荷中心四座站里张北、康保是典型的新能源汇集站只管把风、光出力变成直流功率送出去丰宁接的是抽水蓄能既能从直流电网取电抽水也能放水发电向电网送功率属于天然的功率调节节点北京站是受端整个直流电网最终靠它把功率馈入交流主网。这个分工直接决定了后面控制模式的分配送端站不能承担定直流电压的任务而受端和抽蓄站必须参与电压支撑否则整个四端网络没有电压参考点仿真必崩。1.2 环形拓扑带来的三个“多端特有”问题为什么张北柔直不敢做成简单的两端直流因为环形结构最大的价值是可靠性——任意一条直流线路故障功率可以从另一条路径绕行不会导致整个系统停运。但环形也带来了三个必须用仿真去死磕的问题第一是潮流分布不再天然固定。环形直流电网里各条线路的功率不仅取决于换流站指令还受线路电阻影响不是想让哪条线送多少就送多少。第二是直流电压全网统一任何一点功率不平衡都会直接体现在电压升降上站间交互影响远比两端系统复杂。第三是故障隔离难度大一处直流故障如果不快速切除其余三站都会向故障点馈入电流必须靠换流站控制和直流断路器联合完成故障穿越。这三个问题恰恰是四端仿真最有价值的地方。两端系统里不存在“功率绕行”更不存在“多条路径同时馈流”只有把模型搭成四端这些动态才能真正看出来。1.3 仿真建模必须做哪些“失真”简化工程是双极结构每个换流站有正负两极、极线和金属回线每个桥臂的子模块数在数百到上千量级。如果完全按工程原样建模四座站、每站6个桥臂、每个桥臂500个子模块整个模型就是上万个详细开关器件。在PSCAD里步长开到微秒级这种规模基本属于不可解。所以我在搭建时做了三层简化每个层次保留的动态不同第一把双极系统降为单极等效模型。研究电压平衡和极间故障时才需要恢复双极研究站间功率控制和一般性短路故障单极完全够用。第二把每个桥臂的子模块数从工程量级降到50到200个。子模块数减少后NLM调制和均压算法仍然保留直流电压纹波、桥臂环流、电容电压波动这些关键动态都还在只是开关频率特征和实际工程有一定偏差。第三用戴维南等效模型替代详细IGBT开关模型每个子模块等效为一个受控电压源串联电阻仿真速度能快一到两个数量级。一句话总结我的建模原则模型精度跟着研究问题走而不是越大越好。研究控制策略用平均模型研究故障暂态和子模块均压用戴维南等效模型只有研究器件级损耗和开关应力才需要上详细模型。2. PSCAD中MMC的模型选型、桥臂参数与调制实现MMC是四端柔直的核心也是PSCAD仿真里最吃功夫的部分。很多新手一上来就照着论文里的公式堆参数结果要么仿真慢到跑不动要么波形全是振荡。这里我把自己在PSCAD里折腾MMC的经验捋一遍。2.1 PSCAD中四种主流MMC建模方式对比PSCAD里建MMC主要有四类思路适用场景完全不同模型类型保留的动态相对速度适用场景详细开关模型IGBT开关过程、电容充放电、均压逻辑全部保留很慢几十个模块的小规模验证开关函数/平均值模型换流器外特性、控制器响应很快系统级长过程、机电暂态戴维南等效模型电容电压波动、子模块投切等效中等大规模MMC工程暂态与控制研究PSCAD库自带MMC组件可配置子模块数、桥臂参数中等偏快工程场景快速复现我做张北柔直四端模型时选的是戴维南等效路线。它的核心思想不难理解每个子模块在工作时只看两个状态——投入或切除投入时电容接入桥臂切除时电容被旁路。把IGBT和二极管组合看成可变电阻再结合电容的离散化模型每个子模块就能化简成一个受控电压源串联一个等效电阻。这样做的好处在整个四端网络里非常明显单个桥臂的所有子模块可以合并成一个可变电压源和可变电阻串联六个桥臂最后只解6个状态方程而不是上万个。我在PSCAD里用自定义模块实现这个逻辑桥臂电流每步更新一次子模块电容电压每步更新一次比详细模型快得多。平均值模型我也试过速度快到几乎不占资源但代价是看不到电容电压波动和桥臂环流故障电流的上升形态也不够真实。所以我的建议是做启动过程和功率阶跃用平均值模型做故障和保护研究必须换戴维南等效模型。2.2 桥臂电感和子模块电容的工程估算方法确定桥臂电感和子模块电容是MMC建模的第一步也是最容易拍脑袋的一步。这两个参数直接决定环流大小、电压纹波和故障电流上升率不能随便填。桥臂电感的主要作用是抑制桥臂间的二倍频环流、限制直流短路电流上升速度同时滤除部分谐波。工程上习惯按标幺值取0.1到0.2pu我仿真时一般取0.15pu作为初值。折算到实际值时可以先用这个思路估算阻抗基准Z_B U_N^2 / S_N桥臂感抗X_L 0.15 Z_B再除以角频率ω就是电感值。比如一个1000MW、±500kV等级的换流站折算到单桥臂电抗几十毫亨左右具体数值会随电压容量基准改变。子模块电容的估算可以用储能时间常数法。公式是C ≈ (2 × H × S_N) / (N × U_C^2)其中H是桥臂储能时间常数常规柔直工程大致取30到40msS_N是换流站额定容量N是每个桥臂的子模块数U_C是子模块额定电压。我举个简便算例S_N1000MW、N200、U_C2.5kV、H0.035s代入得到C约56mF。这个数量级和实际柔直工程的子模块电容参数是吻合的。参数选完后要回去看波形验证稳态下子模块电容电压波动率一般控制在±5%到±10%超过10%说明电容偏小或桥臂电感偏小二倍频环流若没有明显被抑制则要调大桥臂电抗或检查环流抑制控制器。2.3 NLM最近电平逼近与子模块均压排序的实现张北柔直的子模块数量是数百到上千量级这种规模用载波移相PWM不太现实工程和仿真里基本都用最近电平逼近调制NLM。核心逻辑很直接每相桥臂需要的电平数由参考电压除以子模块平均电压后取整得到即n round(u_ref / U_C)然后根据桥臂电流方向决定投入哪些子模块。这里最关键的是子模块电容均压排序。PSCAD是Fortran脚本环境对数组循环处理能力有限如果每个仿真步长都做全排序模型会慢得让人崩溃。我在自定义模块里写的是经典排序均压逻辑! 桥臂电流方向判据 if (i_arm 0.0) then ! 电流正向电容被充电优先投入电压最低的n个子模块 call sort_vc(vc_array, sort_index, N_SM) do k 1, n on_signal(sort_index(k)) 1 end do else ! 电流反向电容放电优先投入电压最高的n个子模块 call sort_vc(vc_array, sort_index, N_SM) do k 1, n on_signal(sort_index(N_SM - k 1)) 1 end do end if这段逻辑的物理意义很清晰电流给电容充电时让电压最低的子模块多充电尽快抬平电压电流让电容放电时让电压最高的子模块多放电把电压降下来。这样一轮排序下来所有子模块电压都会被往平均值方向拉。实操中要注意两点一是没必要每个仿真步长都排序可以每隔几十微秒排序一次中间保持开关状态不变速度提升非常明显二是桥臂子模块数降到50到100个时均压效果仍然不错但降到二三十个以下时NLM量化误差变大波形毛刺会明显增多建议不要低于这个下限。2.4 二倍频环流的机制与抑制控制MMC的一个经典问题是桥臂间会流过二倍频环流。原因是三相桥臂的瞬时总电压不完全均衡各相桥臂之间产生流通路径环流虽然不经过直流侧但会增大桥臂电流有效值提高损耗和电容电压波动。抑制方法是在每个桥臂的调制波上叠加一个二倍频负序补偿电压。PSCAD里的实现方式是先计算每相上、下桥臂电流之和的一半得到环流分量再用旋转角速度为2ω的dq坐标变换把二倍频环流变成直流量通过PI控制器输出补偿电压。我把这部分做成独立的环流抑制控制器接在每相调制波输出端。实际调节时环流抑制PI参数不需要特别精确因为它的带宽比主控制回路低只要能把二倍频分量压到额定电流的5%以下就算合格。判断依据就看波形桥臂电流里二倍频分量明显减小且不影响直流电压和输出功率。3. 四端直流电网的多站协同控制不要让定电压站“打架”四端柔直和单端MMC最大的区别就是控制模式必须分配好。谁定直流电压、谁定功率、谁做备用支撑这个分工一旦弄错仿真模型轻则电压发散重则启动就过流跳闸。3.1 各站控制角色如何分配我参考张北柔直的工程定位按如下方案分配控制模式张北站定有功功率P、定无功功率Q作为送端功率源不参与直流电压调节。康保站定有功功率P、定无功功率Q同属新能源汇集站。丰宁站采用下垂控制或定直流电压作为功率调节和电压支撑节点。因为抽蓄机组可以快速增减出力这个站天然适合充当电压平衡角色。北京站定直流电压Udc、定无功Q是整个直流电网的主电压站负责吸收送端功率与受端负荷间的差额。实际工程里不会只有一个站承担电压调节任务那样可靠性太低所以我通常在丰宁站再叠加一个直流电压裕度控制或下垂控制作为后备。二选一的话仿真里我会固定用“一个主电压站一个下垂站”的组合既能保证稳态电压精度又能分摊功率波动。3.2 PSCAD中外环控制器的模式切换与抗饱和PSCAD里MMC的控制器常规定义为内环电流控制加外环功率/电压控制。内环d轴电流id控制有功通道q轴电流iq控制无功通道外环则根据站的定位选择不同模式。比如北京站的有功外环给定的是直流电压参考值电压偏差经PI调节得到id_ref张北站的有功外环给定的则是功率参考值功率偏差经PI得到id_ref。最容易踩坑的是模式切换瞬间的PI输出跳变。很多同学做一个“从定功率切换为定直流电压”的逻辑结果切换瞬间电流指令直接跳到限幅值系统过流。问题不在控制结构而是PI控制器没有做抗积分饱和和初值跟踪。我在PSCAD里给每个PI模块设计了跟踪模式当该外环不参与控制时PI输出始终跟随当前id_ref当模式切换发生切换前PI输出已经被顶到和当前指令一致就不会产生阶跃跳变。此外限幅必须和直流电压过压保护逻辑联动。四端系统功率分布突变时定电压站要能吸收大量过剩功率id的负向限幅如果设得太小直流电压会迅速上升然后触发过压保护跳闸看起来像保护误动实际是控制器限幅不够。3.3 模拟新能源功率波动与潮流反转的仿真方法张北柔直一个非常重要的运行场景就是新能源出力的剧烈波动。风和光从满发到零发可能只需要几分钟在电磁暂态仿真里可以把它压缩到秒级。我的做法是给张北站和康保站的有功功率参考值P_ref接一个可控信号源初始给满功率1秒后以斜坡降到02秒时改为负值模拟功率从受端反送回来。这时观察两件事一是直流电压是否保持在合理范围二是各站电流方向是否发生反转。功率反转在PSCAD里的典型表现是直流电流方向反向换流器触发角等效地从整流状态切到逆变状态站的有功功率符号翻转。如果定电压主站的PI参数太慢直流电压会出现一个明显的暂升尖峰时间大概在几十到几百毫秒。这个尖峰如果超过保护定值就会触发过压保护所以在仿真里做功率反转场景时一定要把直流电压主站的带宽和限幅设计得比功率站动态慢半拍避免站间振荡。3.4 控制参数整定的工程经验PI参数不需要在仿真里靠试错试半天用一阶惯性整定法就能得到合理的初值。内环电流PI可以取Kp L/τKi R/τ其中L、R是桥臂到换流器端口的等效电感和电阻τ是期望的电流环时间常数一般取0.5到2ms。外环比内环慢一个数量级电压环或功率环的带宽取20到50Hz即可。多端系统里直流电压外环要尤其注意直流网络本身有较大的等效电容电压环带宽过高会导致站间功率振荡。我做四端仿真时的经验是主电压站的外环带宽要低于送端功率站电流环带宽的一半以上才能保证直流电压平稳。4. 直流线路故障仿真短路电流每个阶段怎么变化搞清了稳态控制故障工况才是柔直仿真最刺激的部分。张北柔直用架空线作为直流输电线路故障率比电缆高得多所以保护配合的研究在整个工程里分量很重。4.1 直流线路故障场景怎么在PSCAD里设置PSCAD里直流线路我用Bergeron模型或相位域模型搭建架空线路参数按实际杆塔结构设置如果不做精密的行波保护研究用集中参数模型也够用。故障投入用Fault模块支持设置故障时刻、故障持续时间、过渡电阻以及故障类型。对四端直流电网我一般重点做三类故障单极接地故障最常发生故障电流相对较小对保护灵敏度要求高。双极短路故障最严重短路电流上升率极高主要考核直流断路器和换流器过流能力。经高阻接地故障考验保护能否识别过渡电阻我常设置10到100欧姆用来测试故障检测的灵敏度区间。设置故障前还要检查换流站主变压器的接线方式。如果变压器绕组接线导致零序回路不通单极接地故障的动态会跟实际工程差很多我就因为这个问题浪费过好几天——波形看起来像那么回事实际物理过程是错的。4.2 短路电流的发展阶段与MMC响应曲线判读MMC直流侧发生双极短路后短路电流的发展大致可以分成三个阶段。第一个阶段是子模块电容放电阶段。故障瞬间所有投入的子模块电容通过故障点快速放电直流电流在几毫秒内急剧上升这个阶段的电流上升率主要由桥臂电感和故障线路电感决定子模块电容越大放电越剧烈。波形上表现为直流电压瞬间塌陷电流直线爬升。第二个阶段是交流系统馈入阶段。随着子模块电容电压跌落换流器闭锁条件逐渐逼近交流系统开始通过续流二极管向直流故障点馈入电流这个阶段的电流由交流侧阻抗决定上升速度变慢但幅度可能继续增大。第三个阶段是换流器闭锁后的衰减阶段。闭锁后子模块停止放电半桥子模块情况下故障电流只能通过续流二极管自然衰减衰减速度取决于直流回路的时间常数全桥子模块情况下闭锁后电容反向接入会产生反电动势主动压制故障电流电流下降速度明显加快。在PSCAD波形上判断故障特性和保护配合有效性的关键就是找电流变化率的拐点。如果拐点在闭锁指令之前就出现说明故障电流已经开始受限如果拐点在闭锁之后才出现说明保护动作慢了。做故障分析时不要只看电流峰值更要把这个时间轴捋清楚。4.3 直流断路器与换流器故障穿越的配合逻辑张北柔直工程最亮眼的装备之一就是高压直流断路器。工程资料显示采用了兼具快速分断和大电流开断能力的高压直流断路器开断时间在毫秒量级。仿真里我在每条直流线路两端都配置了断路器组件用自定义控制逻辑驱动。保护逻辑我在PSCAD里按如下时序实现检测到直流电流超过过流阈值或di/dt超过变化率阈值判定为故障经过一个设定的延时比如2到3毫秒给故障线路两端的直流断路器发跳闸信号断路器在几百微秒到几毫秒内完成分断故障电流转移到避雷器耗能支路直到电流衰减到零。这里有个特别重要的细节四端环形结构决定了一条直流线路跳闸后功率会从另一条路径绕行所以不能只看故障线断路器能不能断开还要看非故障线路是否过载。我在仿真里经常设置线路1故障后观察线路2电流是否在允许范围内如果超过1.2倍额定值就要考虑系统级联切机或者调整控制策略这在两端直流仿真里完全不需要考虑。4.4 保护定值不能拍脑袋的两个原因保护定值的校验是故障仿真最后的落脚点。我的体会是直流保护阈值设置一定要用仿真波形反推不能直接写“1.5倍额定电流”。原因有两个。第一故障电流上升率极高用单纯的过流判据会存在动作延时故障电流可能轻松冲到好几倍额定值等到电流越限再动作断路器分断压力巨大。通常要配合di/dt快速判据在电流还处于上升初期就给跳闸信号。第二经高阻接地故障的短路电流可能很低过流判据灵敏度不足必须叠加低电压判据或行波判据。我在仿真里经常做一组故障电阻扫描从0欧姆逐步做到100欧姆看保护在哪个点失效这样才能确定定值的边界。5. 从搭模型到调通的完整流程与踩坑清单这一节专门写给那些想在PSCAD里完整复现张北柔直四端模型的人我把自己几个星期的调试过程压缩成一套可复用的流程和经验。5.1 建议的模型搭建顺序先单站后双端再四端不要一上来就搭四端全模型那是自找苦吃。我建议严格按下面的顺序分阶段搭建和验证。阶段一单站MMC接入交流电网。验证目标是解锁不发生过流直流侧电压能稳定建立有功功率可调稳态波形正常。这个阶段重点解决MMC内部参数、调制和均压的问题是整个工程的基石。阶段二双端背靠背柔直。一端定直流电压一端定功率验证直流电流传输、功率反转、启动流程。这个阶段重点解决定电压站和定功率站的协调以及PLL在弱交流系统下的稳定性。阶段三三端或四端扩展。逐步把张北、康保、丰宁、北京四站都接入先把各站设为定功率只让一个站定电压看看系统能不能稳定再引入下垂站和裕度控制测试功率阶跃和负荷变化。最后再把直流线路故障和断路器逻辑加进去。每个阶段转移前都要把当前阶段的波形指标固化下来比如直流电压波动小于5%、功率跟踪误差小于2%、第一次故障切除后非故障站能够继续运行。指标不过就回头检查不要在未达标的基础上继续往上叠复杂度。5.2 仿真步长怎么设置才能兼顾速度和稳定性PSCAD仿真步长是影响MMC模型成败的关键参数之一。步长太大数值振荡明显NLM的电平切换点不准确步长太小模型跑一步要好几分钟根本做不了参数扫描。我常用的配置如下表可以按研究目的直接套用仿真目的主步长说明平均值模型系统级动态20~50微秒速度快适合长时间功率阶跃戴维南等效MMC暂态仿真10~20微秒兼顾MMC动态和整网暂态详细开关模型故障暂态1~5微秒需要捕捉IGBT和故障细节小步长MMC子网0.1~1微秒配合PSCAD小步长子系统使用我实际跑四端戴维南模型时主步长用20微秒桥臂子模块数降到100左右一个秒级故障仿真大概需要跑几十分钟在可接受范围内。如果发现电压波形有高频毛刺先不急着缩步长检查是不是均压排序频率不够或线路模型太复杂。PSCAD里还有一个技巧是把MMC阀组放进小步长子网主步长可以放宽到50微秒小步长部分用几百纳秒到2微秒这样不仅算得快数值稳定性也更好。5.3 那些常见的报错、怪波形和调试方向我调试四端模型时碰到过不少问题印象最深的是这几类。直流电压建立不起来。如果四座站全部设成定功率系统没有任何一个站主动建立直流电压仿真必然发散。检查思路是确认至少有一个站处于定直流电压模式或V/f模式并且这个站先于其他功率站解锁。启动瞬间过流。很多情况下是PLL没有锁定就被强制解锁或者PI控制器初值没对齐。我后来把解锁条件设成“直流电压达到0.9pu以上才允许解锁送端站”过流问题基本消失。高频数值振荡。波形出现密密麻麻的细小毛刺通常是步长和桥臂电抗、电容参数不匹配导致的数值问题。优先检查桥臂电感是否太小再尝试加一个很小的阻尼电阻最后才考虑缩步长否则仿真速度损失太大。断路器跳开后直流电压仍然持续上升。这是环形网特有的问题故障清除后送端功率没有及时降下来功率积在直流网络里导致电压越限。解决方向是功率站和直流断路器联动跳闸信号直接给送端站发降功率指令这叫故障闭锁联动实际工程里也一样。5.4 模型验证拿什么跟真实工程对标很多同学跑完模型都不知道自己搭得对不对。张北柔直工程公开数据有限但足够用来做稳态验证±500kV直流电压、9000MW总换流容量、四座站各自容量、线路总长等。把这些指标作为系统稳态的对照基准看看空载和满载时各站功率分配是否一致。故障暂态方面可以跟学术论文里已经发表的张北柔直仿真参数和波形特征做定性对比故障电流上升率、断路器动作时间、电容电压跌落幅度、故障清除后的恢复时间。我在文章末尾参考了几篇公开论文的参数表和波形数据用来校准自己的保护阈值和故障电阻范围。注意不要直接搬运完整模型代码参数和波形特征借鉴是可以的。另外一个很有用的验证手段是场景互测比如让张北站功率从满发降为零观察系统能否完成潮流反转再断开一条直流线路看其余三站能否继续保持运行。这些场景如果都能跑通模型的工程可信度就有保证了。我实际调试中最喜欢做的一个实验是断开丰宁站到北京站之间的某条直流线路然后观察功率从剩下的路径绕行时各站直流电压短暂凹陷后又恢复的过程。那一刻你会直观地理解“环形柔直电网”到底是靠什么在支撑可靠性也更能体会张北柔直工程里那些看似复杂的技术方案每一个都是被真实物理需求逼出来的。仿真做完了再回头去看这个工程很多当时觉得晦涩的术语都变成了屏幕上某条波形的某个拐点这就是做工程仿真最上头的部分。如果后面你想继续深入我建议顺着这个模型往两个方向扩展一是给张北、康保站接入完整的风电场和光伏动态模型做新能源出力的随机波动影响分析二是把丰宁抽蓄的机组模型加进去研究抽蓄作为电压支撑节点在功率缺额时的快速补偿能力。先把手头的四端基本流程跑通再往上面叠加复杂模型你会比直接啃大模型轻松得多。