做LCC-HVDC仿真这件事最难的不是MATLAB操作而是搞清楚“这个参数是怎么来的”。我最近完整搭了一套500kV直流电压等级、两侧交流系统345kV的LCC-HVDC输电仿真系统功率等级按1500MW设计控制侧采用经典的整流站定电流、逆变站定熄弧角结构再配上VDCOL低压限流环节。整个过程从换流变压器变比计算开始到滤波器组参数推导再到故障恢复波形验证踩了不少坑也把很多原理层面的东西真正串起来了。这篇文章就把这套500kV LCC-HVDC仿真模型的搭建过程、参数计算和控制思路完整写出来给正在做高压直流输电课程设计、毕业设计或者刚接触直流输电仿真的同学做个参考。1. 为什么是500kV LCC-HVDC加345kV交流侧先搞清方案的工程基线1.1 这个参数组合背后对应着什么样的实际工程500kV直流电压等级在电网里是非常经典的一档。国内很多早期高压直流工程都落在±500kV送电容量通常在一千到三千兆瓦这个区间。两侧交流系统用345kV这个电压等级在国际上非常常见尤其是区域电网之间的互联场景。把直流侧额定电压定为500kV、交流侧额定电压定为345kV本质上是在模拟一个“区域电网对区域电网”的LCC-HVDC输电通道。LCC是Line Commutated Converter的缩写也就是电网换相换流器核心器件是晶闸管。它靠交流系统自身的电压完成换相所以对交流系统强度有依赖这也是为什么在做这类仿真时交流系统短路比、无功补偿配置、滤波器设计会那么重要。相比之下VSC-HVDC用的是可关断器件但LCC在超高电压、超大容量场景下依然有损耗低、技术成熟的优势现在很多实际工程仍然围绕LCC展开。如果你要做的是毕业设计或课程设计我建议先把这套系统理解成三块整流站负责把交流功率送到直流线路逆变站负责把直流功率送进受端交流电网中间是直流输电线路本身。两侧交流系统电压相同都是345kV但整流侧和逆变侧的换流变压器变比、触发角工作点、无功补偿容量会有差异这些正是仿真的价值所在。1.2 整流站、逆变站、直流线路三大块的拓扑骨架完整模型里整流站交流母线进来之后先经过换流变压器变压器二次侧接两组6脉动换流桥一组通过星型/星型变压器供电另一组通过星型/三角型变压器供电两组桥在直流侧串联就构成12脉动换流器。12脉动的好处是交流侧特征谐波变成了12n±1次也就是11、13、23、25这些次数相比6脉动的5、7、11、13…滤波压力小很多。直流侧从换流桥出来之后要接平波电抗器平波电抗器的作用是抑制直流谐波电流、防止电流断续一般每站取0.5到1亨。然后才是直流输电线路。逆变侧的结构和整流侧对称同样有换流变压器、两组6脉动桥、平波电抗器但控制目标完全不同。在我这个模型里直流线路用的是集中参数RLC模型线路长度按500km设置回路总电阻约8Ω电感约1H电容约2μF。这个参数不是某个具体工程的精确值而是课程设计级别的合理简化。如果你要研究线路波过程就得换分布参数模型但做稳态和故障恢复分析集中参数足够。1.3 从额定功率到额定电流的参数推算系统额定功率我取1500MW。直流电压500kV所以额定直流电流就是Id_N 1500MW / 500kV 3000A这个电流水平对晶闸管阀组来说非常正常也为后面计算换流变压器容量、滤波器容量提供了基准。接下来是整流侧触发角。LCC正常运行时空载整流侧触发角取15°左右逆变侧熄弧角取17°左右。为什么要这么取触发角太小会导致直流电压过高、无功消耗特性变差触发角太大又会压缩逆变侧的调节裕度。15°和17°是工程上很成熟的折中选择。再往里算就需要用到12脉动换流器的理想空载直流电压公式Udo (3√2 / π) × Uvll其中Uvll是换流变压器阀侧线电压有效值系数约等于1.35。整流侧直流电压方程是Ud Udo × cos(α) - (3/π) × Xc × Id后面那一项叫换相压降是重叠角引起的等效电压损失。取换相压降约为Udo的5%整流侧额定电压500kVcos15°约0.966于是Udo_rect ≈ 500 / (0.966 - 0.05) ≈ 546kV阀侧线电压Uvll_rect Udo_rect / 1.35 ≈ 404kV所以整流侧换流变压器变比大约是345kV比404kV折算成标幺变比约为0.854。这里有个很容易忽视的点整流侧和逆变侧的阀侧电压并不一样。因为直流线路有压降逆变侧直流电压会低于500kV假设线路压降24kV逆变侧电压约476kV再按cos17°和5%换相压降反推逆变侧Udo约525kV阀侧线电压约389kV变比约0.887。两侧变比不同怎么办实际工程中换流变压器都配有载分接头就是用来补偿直流线路压降和两侧电压差异的。仿真模型里我直接把两侧变压器变比设成各自的计算值运行效果更稳定也更容易复现。2. 换流变压器与换流桥搭建一次系统模型的正确顺序2.1 阀侧电压怎么算触发角、熄弧角与Udo的关系很多新手上来就是“电源接变压器变压器接整流桥完事”但根本不解释为什么变压器二次侧选这个电压。实际上阀侧电压完全是由功率等级、直流电压、触发角和熄弧角共同决定的。整流侧计算里已经展示了这个逻辑500kV直流母线电压要维持住触发角15°决定了cos项只有0.966换相压降还要吃掉5%于是Udo必须比500kV高出约9%才能保证满载时有足够的“电压头”去克服损耗。逆变侧逻辑完全相同只不过熄弧角17°对应的cos值是0.956比整流侧略低一点再加上直流线路压降的因素所以逆变侧阀侧电压会略低。这个计算过程请务必自己按计算器推一遍不要只抄结果。因为一旦你改了功率等级或直流电压所有阀侧电压、变比、滤波器容量都要跟着变只改一个数必然出问题。我在模型里用的变压器参数如下表位置交流侧电压(kV)阀侧线电压(kV)额定容量(MVA)短路阻抗整流侧345404750×20.15 p.u.逆变侧345389750×20.15 p.u.每站两组换流变压器对应两组6脉动桥每组容量750MVA合起来正好和1500MW的输送容量匹配。短路阻抗取0.15p.u.这是LCC换流变压器的典型值太小会加大换相压降和短路电流太大会限制直流电压调节能力。2.2 12脉动换流桥在MATLAB里的组装方式MATLAB里搭12脉动换流桥有两条路。第一条是用现成的例子改。MATLAB的Simscape Electrical库里有基于CIGRE HVDC Benchmark的示例模型里面就是标准的12脉动LCC整流—逆变系统。直接打开这个模型把交流电压改成345kV、直流电压改成500kV、功率改成1500MW然后按前面算的参数替换变压器和滤波器这是最快的方式。适合赶时间、或者第一次接触想先看整体结构的同学。第二条是从零搭用Universal Bridge模块把器件类型选成Thyristor然后搭两组桥。注意两组桥必须分别接到星型和三角型接法的换流变压器上。这个30°的相位差异会直接体现在触发脉冲上Y/Y组和Y/△组的触发角要错开30°电角度。我自己的建议是如果目标是写论文或者做毕设第二条路更值得走。因为用现成模型改你对触发脉冲、换相过程、控制信号的产生机制基本没有体感评审老师随便问一句“逆变侧怎么决定触发时刻”就答不上来。从零搭一遍哪怕只搭6脉动再扩展成12脉动都比直接改demo更有收获。触发脉冲这块MATLAB里最实用的做法是用PLL锁住换流母线电压相位再根据alpha角的控制信号给每支晶闸管生成对应的触发脉冲。触发角alpha的参考点在相电压过零点之后滞后alpha度。6脉动桥的六个晶闸管依次间隔60°12脉动桥因为有两组间隔变成30°。2.3 平波电抗器与直流线路的取值逻辑平波电抗器是LCC直流侧最容易忽略但影响很大的元件。它串在换流器和直流线路之间主要作用是抑制直流侧谐波电流、防止小电流工况下的电流断续。电感取大了动态响应变慢故障后电流恢复时间拉长取小了谐波抑制不住逆变侧还容易换相失败。我在模型里每站取0.6H两侧加起来1.2H配合线路电感1H直流回路总电感大约2.2H。这个数量级在LCC里是合理区间直流电流的时间常数大概在几十毫秒级别恢复过程波形不会太夸张。直流线路用RLC串联支路模拟参数前面说了8Ω电阻、1H电感、2μF电容。电容虽然只是集中参数的近似但它在直流短路故障后的暂态里作用不小不能完全省掉。还有一个小坑如果你用Universal Bridge直接搭晶闸管的缓冲电路参数默认值有时候会导致仿真速度极慢。如果发现仿真跑不动把Powergui里的求解器换成离散模式、步长取50微秒性能会好很多。高压直流系统的电磁暂态时间常数跨度很大连续求解器在这种模型上经常卡死。3. 控制系统才是LCC仿真的灵魂定电流、定熄弧角与VDCOL3.1 整流侧定电流控制器的结构与PI整定LCC直流输电为什么整流侧要用定电流控制因为晶闸管承受过流能力很弱故障时若电流不加限制器件很快就会损坏。通过快速调节触发角让直流电流始终稳定在参考值附近是保护设备的第一道防线。控制结构是典型的单闭环直流电流实测值Id_meas和参考值Id_ref比较误差进PI调节器输出就是触发角alpha命令。alpha命令经过限幅环节后送到触发脉冲生成器。PI参数怎么整定换流器本身的响应速度很快从触发角变化到直流电流变化大约有1到3个毫秒的延迟整个对象近似成一阶惯性加纯延迟。我用的是工程整定法先把PI的积分时间放在0.02秒附近比例增益从10开始往上加观察直流电流阶跃响应直到它出现一次小幅超调再回调一点。实际模型里我最后用Kp大约20、Ki大约500稳定性和响应速度比较均衡。有一个很多人会犯的错误整流侧alpha限幅直接设成0°到90°。这听起来合理因为整流运行确实不会超过90°但在故障恢复和启动阶段整流侧可能会短暂进入逆变运行给后面的恢复过程帮倒忙。所以我把alpha下限设为5°、上限设为150°。这个上限给控制留出了更大的回旋空间实际运行中很少触发但关键时候能避免控制饱和。3.2 逆变侧定熄弧角控制与电流裕度切换逆变侧的控制目标比整流侧复杂一层。正常工况下逆变侧需要维持足够大的熄弧角γ防止晶闸管换相失败。γ太小关断后的晶闸管来不及恢复阻断能力下一次换相还没开始就被重新导通。工程经验是γ取17°到19°我模型里定的是17°。定熄弧角控制的实现也不难实时估算逆变侧晶闸管的熄弧角和17°参考值比较误差进PI输出调节逆变侧触发角。逆变侧触发角正常运行在140°到145°附近对应超前角β大约35°到40°。这中间包含重叠角20多度加上γ裕度17度正好符合物理约束。比定熄弧角更重要的一个细节是电流裕度控制。直流输电必须保证整流侧和逆变侧不会同时抢功率。标准做法是整流侧电流参考值设在Id_ref逆变侧电流参考值设在Id_ref减去电流裕度ΔI工程上一般取0.1p.u.也就是0.1倍额定电流。正常情况下直流电流由整流侧决定逆变侧只负责定γ。一旦整流侧电压太低、无法维持电流直流电流就会下降当它跌到逆变侧参考值以下逆变侧就会自动切到定电流控制模式把直流电流稳定在略低的水平。这个机制保证了系统在任何情况下都有且只有一个站主导电流控制。如果你做仿真时发现逆变侧经常失控大概率就是没做这个电流裕度切换。我建议哪怕课程设计也要加上只做定γ控制的模型在直流短路和启动工况下会有很多不真实的发散现象。3.3 VDCOL低压限流环节的实现方式VDCOL全称Voltage Dependent Current Order Limiter中文叫低压限流环节。它的作用一句话概括直流电压低的时候主动降低电流参考值。为什么需要它直流侧发生短路或者逆变侧换相失败时直流电压会迅速跌落。如果此时整流侧还按额定电流给定参考值直流电流会冲得很高晶闸管承受很大的热应力故障后系统还容易在一连串换相失败里反复折腾。VDCOL通过监测直流电压把它变为电流参考值的限制信号故障期间主动降压运行等直流电压恢复了再逐步把电流参考值抬回去这样恢复过程会平滑得多。在MATLAB里实现VDCOL很简单用Lookup Table模块就够了。我用的分段线性关系如下直流电压(p.u.)电流参考值限制(p.u.)0.00.30.40.30.91.01.21.0也就是说直流电压高于0.9p.u.时电流参考值不受限低于0.4p.u.时被钳位到0.3p.u.中间线性过渡。这个参数不是固定的不同的系统强度可以调节。交流系统比较弱的话低压段的电流限制要取更低比如0.1p.u.不然恢复期间逆变侧还是会反复换相失败。VDCOL放置的位置要注意整流侧和逆变侧都有电压测量但电流参考值限制应该作用在主控制级的整流侧电流定值上。想简单一点就在直流母线测电压、在整流侧控制器的Id_ref输入端加一个limiter模块逻辑完全一致。4. 交流侧滤波与无功补偿容易被忽视的隐形工程量4.1 特征谐波次数12脉动换流器到底产生什么LCC换流器在换相过程中会在交流侧产生特征谐波电流。12脉动换流器的特征谐波是12n±1次也就是11、13、23、25次谐波电流幅值近似按1/n基波电流递减。如果不加滤波这些谐波电流流过交流系统阻抗会在电压波形上产生严重的畸变对附近设备有直接危害。直流侧的特征谐波是12n次主要是12次和24次。但由于已经有了平波电抗器和直流滤波器直流侧谐波一般不会成为主要矛盾。很多简化模型干脆不做直流滤波器只靠平波电抗器扛着对仿真结果影响不大。交流侧就不行滤波器必须做。在MATLAB里验证谐波的方法是把仿真数据导到Powergui的FFT分析工具选好基波频率50Hz、窗口时间看各次谐波占比和THD。正常情况下接入完整滤波器组后换流母线的电压THD应该压到5%以下实测一般在2%到3%附近如果看到THD畸高先怀疑滤波器参数再怀疑触发脉冲对称性。4.2 滤波器组配置与参数计算实例滤波器设计要从无功平衡算起。LCC换流器在运行中会吸收大量无功数量大约是有功功率的50%到60%。模型里1500MW的直流功率每站换流器大约需要吸收850到900Mvar无功。这些无功必须由交流侧的并联电容器和滤波器组提供全指望交流电网撑的话换流母线电压根本稳不住。我用的是经典配置两组单调谐滤波器分别谐振在11次和13次一组二阶高通滤波器覆盖24次及以上再配一组并联电容。容量分配大致如下支路类型谐振次数无功容量(Mvar)单调谐滤波器11250单调谐滤波器13250二阶高通滤波器24200并联电容器-200合计-900每站合计900Mvar正好覆盖换流器的无功需求并留有一点裕量。单调谐滤波器参数怎么算以11次、250Mvar这支为例。母线线电压345kV相电压约199kV每相承担的无功约83Mvar。先算基波容抗X_C1 U_ph² / Q_ph 199² / 83 ≈ 477Ω对应的基波电容C 1 / (2π × 50 × 477) ≈ 6.7μF因为谐振在11次11次下的容抗是基波容抗的1/11也就是X_C11 477/11 ≈ 43.4Ω。谐振时感抗等于容抗所以L 43.4 / (2π × 550) ≈ 12.6mH电阻由品质因数决定Q值取30到50之间我取40R 2π × 550 × 0.0126 / 40 ≈ 1.1Ω把这个算例里得到的R、L、C数值填进Three-Phase Series RLC Branch就完成了一个单调谐滤波支路。13次和高通支路同理。我用表格整理一下具体参数支路R(Ω)L(mH)C(μF)11次单调谐1.112.66.713次单调谐0.99.26.724次高通较高阻尼约3.57.0并联电容00约5.3并联电容支路的容量按200Mvar一相约67Mvar折算相电压199kV下容抗约591ΩC约5.4μF。注意并联电容上一般不串电感和电阻直接接母线即可。4.3 无功功率差额与交流电压稳定性换流器吸收的无功和滤波器提供的无功不可能时刻相等这个差额直接表现为换流母线电压的升高或降低。系统比较强短路比SCR高的时候母线电压变化不大系统弱的时候一点点无功差额就会引起明显的电压偏移反过来又影响触发角和熄弧角形成恶性循环。我在仿真里两侧交流系统都用短路比衡量取SCR等于5。这个强度属于中等偏强比很多实际弱网条件好一些但足以看出无功补偿的作用。你可以在MATLAB里做这样一个实验把滤波器组全部断开跑稳态然后观察换流母线电压的幅值——大概率会跌到0.9p.u.以下这就是LCC对交流系统强度敏感的直接证据。动手做LCC模型的时候我给一个实际的建议先搭好一次系统和控制后别急着做故障仿真第一步就是校核无功平衡。在整流侧和逆变侧的交流母线各放一个功率测量模块运行到稳态后看无功方向。正常情况下交流系统只提供很小的无功无功主要是滤波器给、换流器吃。如果测量结果显示交流系统输送了大量无功说明你的滤波器和并联电容容量不够母线电压会塌。5. 仿真调试实录从启动振荡到波形稳定的排查过程5.1 稳态启动的正确操作方式LCC-HVDC仿真最容易翻车的就是启动阶段。晶闸管换流器在未触发状态下相当于开路一旦控制系统不合拍第一个周波就能把电压和电流炸到飞起。正确的启动顺序是先不触发任何晶闸管让交流系统单独运行换流变压器带上滤波器组一起做稳态初始化跑0.2秒左右。然后给整流侧一个很小的电流参考值比如0.1p.u.让控制慢慢建立直流电流。等直流侧电流稳定了再把电流参考值斜坡上升到1.0p.u.斜坡可以取0.2到0.5秒。逆变侧在整个过程中保持定熄弧角控制它会自动配合整流侧的电流建立。MATLAB里实现斜坡参考值用Ramp模块就行比直接阶跃温柔得多。阶跃启动虽然也能到稳态但直流电流大概率会出现一个大幅过冲运气不好直接触发保护逻辑或者晶闸管err。做仿真追求的是波形干净没必要用阶跃给自己添堵。启动阶段的仿真步长也需要注意。我用的是离散模式50微秒步长这个精度对50Hz系统来说可以接受。如果你发现波形有毛刺或者控制发散先把步长改成20微秒对比一下看是不是数值问题造成的。很多“控制不稳定”其实只是求解器精度不够尤其在晶闸管换相发生的时刻步长太大会把换相过程抹掉。5.2 三种常见波形异常的原因定位我在调试这套模型时遇到三类问题应该也是每个LCC仿真初学者都会撞上的。第一种是交流母线电压畸变严重THD非常高。排查方向非常明确滤波器没接对。可能是电容电感参数填错也可能是两组单调谐支路的谐振点偏移了。遇到这个问题先把滤波器支路逐个断开看FFT里哪个次数的谐波幅值异常基本就能定位是哪一组的锅。第二种是直流电流低频振荡看起来像是控制不稳定但又不是高频发散。这种多半是PI参数过激进或者VDCOL的恢复斜率和PI参数产生了耦合。我的处理办法是把电流环PI的Ki降一半再试同时把VDCOL从0.4到0.9p.u.的过渡斜率调缓一点。两者都倾向于增加阻尼虽然响应会慢一些但系统稳定性好得多。第三种是逆变侧频繁换相失败。这个最棘手因为换相失败的影响因素很多。先看逆变侧熄弧角裕度够不够17°左右是比较合理的再看逆变侧触发角限幅如果上限不够高比如只到150°正常逆变运行在140°附近没问题但故障后调节余量不足就会失控。再看交流系统强度SCR太低的话母线电压在故障后会大幅下降逆变器根本没法正常换相。最后也要检查触发脉冲与换相电压的参考点是否对齐这个误差哪怕几度都很致命。5.3 直流短路与恢复过程的仿真验证模型调稳之后我做了两个故障验证。第一个是直流线路短路。在直流线路中点设置一个接地故障0.1秒后切除。故障期间直流电压几乎跌到零直流电流由于线路电感和平波电抗器的作用会先冲高然后被VDCOL压回去。故障切除后整流侧电流参考值随着电压恢复逐步回升系统在大约0.3到0.5秒内恢复到额定运行。这个波形就是非常典型的LCC故障恢复曲线值得好好观察。第二个是逆变侧交流母线三相短路100ms。这个故障会直接导致逆变侧换相失败直流电流在故障瞬间迅速上升直流电压大幅下降。由于整流侧还维持着电流控制加上VDCOL的限制直流电流不会失控。故障切除后逆变侧重新建立换相系统同样需要几百毫秒才能完全恢复。做完故障仿真一定要看VDCOL环节的电流参考值波形它在故障期间应该明显被压低。如果电压都跌到0.2p.u.了你发现参考值纹丝不动说明VDCOL接线有问题恢复过程一定会过冲。我在实际调试中发现故障恢复波形对控制器的参数非常敏感。同样一套滤波器、同一个故障位置把整流侧PI的Ki从500改成800恢复时间可能缩短一半但超调量会显著增加甚至出现二次换相失败。所以别迷信某一组参数多跑几组扫描用控制器的响应速度和稳定性做权衡才是做仿真该有的习惯。最后再分享一个小技巧把换流变压器的短路阻抗、换流站的Udo、阀侧线电压、滤波器无功容量打印成一张表挂在桌面旁边。LCC的系统所有参数都是环环相扣的改任何一个都会牵连出三四个关联量。我因为懒得重新推算吃过好几次亏记不清哪次是滤波器补偿不够还是触发角工作点偏了后来老老实实把每一步的计算格式都固定下来调试快多了。这套500kV/345kV的模型跑通之后你再回头做800kV或者交流侧500kV的方案其实就是重新走一遍参数推算流程整体框架完全一样。