高压直流输电建模与换相失败仿真：三种模型思路解析

我把当时做课题整理的三套模型思路、信号链路和换相失败分析全部梳理出来了内容比较多但每一步都是实际跑过的可以直接对照复现。1. 先捋清楚三种模型解决的是哪三类问题做高压直流输电课题的人大概率都绕不开三件事搭模型、理信号、分析换相失败。我当初接这个课题时最头疼的不是MATLAB/Simulink操作而是“到底该搭哪种模型”。网上资料不少但要么只给模型不给说明要么给详细模型但完全不讲平均值模型怎么来的。这里把我的思路完整讲一遍。先说结论三种模型对应的研究目标完全不一样。第一种详细模型是标准CIGRE HVDC Benchmark模型结构完整包含晶闸管级别的换流器、换流变压器、直流输电线路、交流系统戴维南等效以及完整的控制系统定电流控制、定熄弧角控制、VDCOL低压限流环节等。这套模型适合研究换相失败的物理机理、阀级暂态过程以及各种控制策略对故障响应的影响。我做过仿真在详细模型下能直接看到每个阀的导通和关断波形触发脉冲异常时能逐一对应到某个阀上。第二种详细模型是在Benchmark基础上的“工程化改造版”。核心换流器结构不变但把交流系统的戴维南等效替换成具体的电网等值模型或者加入交流滤波器、无功补偿装置、调相机等外围设备。有的版本还会把平波电抗器参数、直流线路分布参数改成目标工程的真实值。这种模型适合做“带工程背景”的研究比如评估受端系统短路比变化对换相失败的影响或者验证调相机对逆变站换相电压的支撑效果。第三种是平均值模型Average-Value Model也是很多人容易忽略但实际很好用的一类。平均值模型不模拟晶闸管的开关过程而是用受控电压源和受控电流源基于换流器的稳态方程来等效直流侧和交流侧的电气关系。仿真速度比详细模型快一到两个数量级适合研究长时间动态过程比如连续故障下的恢复特性、多直流馈入系统之间的相互作用或者控制参数整定。选哪种模型本质上是“精度”和“速度”的权衡。阀级机理研究、换相失败瞬间的波形分析必须用详细模型做系统级稳定性分析、控制策略多方案对比、批量仿真用平均值模型会轻松很多。2. 信号全链路从交流母线到直流线路每个量从哪来、到哪去我第一次搭模型时最懵的不是模块怎么连而是Scope里那几十个信号到底分别代表什么、怎么串起来的。这套模型的信号链路其实可以拆成四个层次。2.1 交流侧信号母线电压、电流与锁相环整流侧和逆变侧的交流母线电压Uac和电流Iac是所有控制的基础。Simulink里一般通过三相电压电流测量模块Three-Phase V-I Measurement采集母线信号然后经过坐标变换得到d轴和q轴分量。这里的核心用途有三个一是给锁相环PLL提供相位参考保证触发脉冲与交流电压同步二是计算有功功率和无功功率判断系统的功率传输状态三是作为换相电压的“代名词”因为换相过程中晶闸管换相依赖的就是交流母线电压电压幅值一掉换相失败的隐患立刻出现。我在实际建模时特别留意了PLL的输出。很多初学者把PLL看成“无关紧要的配件”其实PLL的动态响应直接影响触发角的测量精度。Benchmark模型里PLL参数已经调好但如果改动了交流系统的等效阻抗建议重新检查PLL的带宽否则在故障瞬间会出现相位跳变干扰后续的判断。2.2 换流器层信号触发角、熄弧角、换相角换流器层最重要的信号是触发延迟角α、熄弧角γ和换相角μ。这三个量直接决定了直流系统的运行状态。触发角α是控制系统输出的核心指令整流侧一般在5°到20°之间调节逆变侧通常工作在100°到140°之间等效于超前角β在40°到80°。熄弧角γ是逆变侧特有的监视量它表示晶闸管关断后承受反向电压的时间对应的电角度。基准模型里逆变侧配置了定γ控制目标值通常设在15°到20°之间留出安全裕度。说句实在话我开始做模型时对“额定工况γ15°”没有直观感觉后来对比波形才发现γ降到10°以下时阀的关断裕度就非常紧张了触发脉冲稍微偏移一点就会换相失败。如果你用详细模型可以直接在示波器里看到每半个工频周期的γ值抖动那是判断系统健康状态最直观的手段。2.3 直流侧信号直流电压、直流电流与输送功率直流侧的标定量主要是直流电压Ud和直流电流Id。整流侧是定电流控制正常情况下电流指令由上级调度给定整流器的触发角跟随指令调节逆变侧则承担定电压或者定熄弧角的任务。这样安排的原因很实际如果两侧都控制电流系统就没有稳定的电压参考点如果都控制电压电流又失去限制。所以整流侧管电流、逆变侧管电压/熄弧角是这条链路的“宪法”。直流功率Pdc就是Ud和Id的乘积。在平均值模型里功率传输关系可以直接用这个公式做受控源的参考不需要再经过阀级事件模拟。直流线路本身还配有平波电抗器它的作用是抑制直流电流的谐波和故障时的电流冲击。信号流里监视平波电抗器两端电流的变化率往往能提前捕捉到换相失败的前兆——电流瞬间上升时说明逆变侧可能已经发生了阀误导通。2.4 控制信号的闭环流向把四个层次的信号串起来完整闭环是这样的测量环节采集交流母线电压、直流电压、直流电流和熄弧角测量值送入控制器控制器经过PI调节、限幅、VDCOL低压限流等环节输出触发角指令α触发脉冲生成环节根据α和PLL相位产生六路或十二路触发脉冲脉冲驱动晶闸管导通改变换流器输出电压和电流直流线路把这些电气量传到对侧形成反馈。这个闭环里有几条反馈回路要注意电流反馈回路决定了整流器α的增减方向γ反馈回路决定了逆变器触发角的前移后退VDCOL环节则是根据直流电压跌落程度压低电流指令是换相失败抑制的第一道防线。信号流理顺后你再去看模型里的每个模块就明白“为什么这里要加一个常数”“为什么那里要取一个测量值”了。3. 手推核心公式平均直流电压、换相角、熄弧角的来龙去脉模型里那些参数不是拍脑袋填进去的很多Simulink模块的系数都来自换流器的稳态方程。哪怕你用平均值模型也得明白背后的数学关系否则故障工况下算出来的结果对不上你根本不知道错在哪。3.1 平均直流电压公式先看六脉动换流桥。设换流变压器阀侧线电压有效值为U2不考虑换相电抗时空载直流电压为Ud0 (3√2 / π) × U2 ≈ 1.35 × U2这个公式的来源是交流线电压整流后的平均值。六脉动桥在一个工频周期内有六个波头每60°有一个换相点计算正弦波在60°区间内的平均值再乘上对应的系数就得到1.35。实际运行中触发角α不为零晶闸管延迟导通直流电压变成Ud Ud0 × cosα一个12脉动换流器相当于两个六脉动桥串联阀侧绕组分别采用星形和三角形接法相位差30°所以总空载电压翻倍这也是为什么12脉动模型里直流侧电压能到500kV量级。3.2 换相电抗带来的电压损失晶闸管换相不是瞬时完成的。从一相向另一相转移电流需要时间这就是换相角μ。换相期间换流变压器漏抗上的电压降会造成直流电压损失精确的表达式是Ud Ud0 × cosα − (3ωLc / π) × Id其中ω是系统角频率Lc是每相的换相电抗主要是换流变压器漏抗折算到阀侧的等效值Id是直流电流。注意这个式的第二项直流电流越大电压损失越大。这正是故障时直流电流飙升会进一步拖垮直流电压的原因。做平均值模型时受控电压源的输出表达式就是这个公式不需要阀级开关直接用代数方程算Ud这是平均值模型速度快的关键。3.3 换相角μ的求解换相角μ由下面的经典方程确定cosα − cos(α μ) (2ωLc / (√2 × U2)) × Id这个式子的物理含义是换相期间电感中存储的电流变化量必须由交流电压“推动”完成。交流电压越高、漏抗越小换相越快μ就越小直流电流越大需要转移的电流越多换相时间越长μ就越大。3.4 熄弧角γ换相失败最直接的判据逆变侧触发角α通常在90°到180°之间为了方便分析定义超前角β 180° − α。逆变侧阀在关断后还有一段时间承受反向电压这个关断裕度对应的电角度就是熄弧角γ 180° − α − μ β − μ由这个式子可以清楚看到γ减小的途径有三种β减小触发角前移、μ增大换相时间变长。而μ增大又由直流电流增大或交流电压跌落引起。所以交流系统故障→母线电压下降→换相电压不足→μ增大→γ跌到晶闸管恢复时间对应的最小熄弧角γmin以下→换相失败这就是一条完整的事故链。晶闸管恢复阻断能力需要几百微秒对应工频50Hz大约是7°到10°。所以工程上通常取γmin为7°到10°。记住这个数字后面做故障仿真判断时直接用。4. 换相失败的边界条件与分类为什么总是逆变侧出问题整流侧和逆变侧都可能发生换相失败但绝大多数情况下出问题的是逆变侧。原因并不神秘整流侧换相时关断的阀马上进入自然导通的方向很快承受反向电压阻断逆变侧换相电压来自交流系统阀在关断后需要在一定时间内保持反向偏置而这个时间窗口本来就很紧张。4.1 逆变侧换相的脆弱窗口正常运行中逆变侧γ维持在15°到20°。晶闸管自身恢复时间对应的γmin约7°到10°所以正常裕度只有5°到13°的电角度折算成时间大约是1到2毫秒。交流系统随便一个扰动电压小幅跌落γ就可能冲破这条线。我经常用这样的类比解释给同组的师弟听逆变侧的换相过程好比两扇门交替开闭门必须关上并锁住一小段时间才能接受下一次开门如果锁的时间不够γ太小下一次开门时门还自动弹开了阀重新导通这就是换相失败。4.2 从波形上识别换相失败详细模型下换相失败最典型的表现是本该关断的阀在关断后不久又重新导通导致直流电压瞬间跌落甚至出现负压尖峰直流电流快速上升。动手做过仿真的同学一定见过这种波形——直流电压出现一个极窄的凹坑或者反向脉冲随后整条Ud曲线像被“砍了一刀”往下掉。故障模拟时判断换相失败是否发生不能只看一条曲线。我通常同时观察三个量γ波形是否跌破γmin、阀电流波形是否有不该导通的阀带了电流、直流电压波形是否跌落超过额定值的50%。三者同时确认才能下“换相失败”的结论。4.3 换相失败的常见分类按照诱因分最典型的三类第一类是交流系统故障引发。这是最常见的尤其单相接地故障。故障造成换相电压幅值跌落和相位偏移直接压缩关断屏段。仿真时故障点通常设在逆变侧交流母线上故障电阻取0.1Ω到1Ω的小值持续0.05到0.1秒。第二类是直流系统故障引发。比如直流线路对地短路或者平波电抗器参数异常导致直流电流振荡使γ发生剧烈波动。第三类是控制系统扰动或触发系统异常。换流变压器分接头调节不当、锁相环失步、触发脉冲丢失等也会诱发换相失败。这类问题在控制器设计不到位时特别容易出现。按失败次数和范围又可以分成单次换相失败故障后γ短暂跌破随后恢复和连续换相失败一个阀刚恢复另一个阀又失败形成连锁单桥换相失败12脉动机组中一个六脉动桥出问题和双桥同时换相失败两个六脉动桥同时出问题后果严重得多。4.4 换相失败对系统的影响换相失败的直接后果是直流功率传输中断。从交流系统角度看逆变侧换相失败期间直流电流得不到正常转移相当于受端交流系统突然失去了大量有功输入同时换流器还会吸收大量无功造成母线电压进一步下跌。如果受端系统较弱短路比SCR低极可能引发电压稳定问题甚至功率振荡。在仿真里你能看到换相失败后直流功率回落到接近零恢复过程又可能出现有功过冲整个动态过程持续几百毫秒到数秒完全取决于控制器的调节能力。5. 抑制与改进措施的落地路径控制参数怎么调、拓扑怎么选这部分是课题里最实用的一块也是我投入时间最多的部分。抑制换相失败可以从控制策略、无功支撑、换流器拓扑三个层面层层递进。5.1 控制层面的第一道防线最基本的是“最小熄弧角限制”。不管定γ控制怎么调都要在控制器里硬性限制γ最小值不低于一个阈值。我在基准模型里把γ限制设为10°到12°故障瞬间控制器的输出不会让触发脉冲提前到γ低于这个限制的程度。第二道防线是VDCOL低压限流。它的核心思想是直流电压跌落后暂时压低直流电流指令帮助系统快速恢复。VDCOL的典型参数包括电压启动值如0.4pu、电压结束值0.9pu、电流最小限制0.3pu。一条典型的VDCOL曲线电压低于0.4pu时电流指令降到0.3pu电压从0.4升到0.9期间电流指令线性回升到1.0pu。这条曲线的作用原理在于换相失败恢复期间直流电流越小换相所需时间越短γ越高因此压低电流就是在给换相“让路”。我也试过不接VDCOL直接跑故障换相失败后系统恢复明显变慢甚至出现二次换相失败所以VDCOL不是可有可无的功能。第三道更高级的是预测型触发控制。原理是检测交流电压跌落速度和幅值预测γ是否会跌破阈值在换相失败发生前把触发脉冲提前或延迟主动避开危险区。这类策略需要控制器有较快的响应速度和可靠的电压检测环节我在Simulink里用M语言写的S-Function实现过一版比纯Simulink模块搭建更灵活调参也方便。5.2 系统层面的无功支撑方案换相失败期间逆变器需要大量无功而交流系统电压又在下跌如果没有额外的无功源恶性循环就开始了。增加调相机是目前工程上用得较多的方案。调相机的旋转惯量和快速励磁调节可以在几十毫秒内提供无功支撑抑制换相母线电压的跌落深度。仿真里给逆变侧交流母线接一个调相机模块对比故障期间母线电压和γ波形差异非常明显。STATCOM和SVC也能起到类似作用。STATCOM响应速度更快没有转动惯量但容量和造价是制约因素。给交流母线加装并联电容器也能提高换相电压不过电容器属于被动补偿电压越低补偿能力越弱效果不如动态补偿。对比数据放在下面措施响应速度无功输出特性实现难度适用场景并联电容器慢电压平方关系电压越低越弱低强系统、正常运行工况调相机较快强励磁短时无功裕度大高弱受端系统、提升SCRSTATCOM快电压跌落时仍能输出额定无功中动态电压支撑、多直流馈入SVC较快受电压影响但优于电容器中常规动态补偿我建议的研究路径是先用并联电容器补充基础无功再在关键母线配置STATCOM对比不同短路比下的换相失败概率仿真结果差异能直接支撑你的结论。5.3 换流器拓扑层面的根治方向换相失败的本质是晶闸管依赖交流电压换相。如果换相电源不完全依赖交流系统电压问题就从根本上得到改善。电容换相换流器CCC的思路是在换流变压器和晶闸管之间串联电容利用电容电压辅助换相。这样即使交流母线电压跌落了电容积存的能量也能维持换相电动势。我在平均值模型里改过这种结构仿真效果确实比传统晶闸管换流器好得多但串联电容在故障时承受的电压应力很大阀和电容的配合设计要精确。VSC-HVDC则彻底绕开了换相失败问题。电压源换流器使用IGBT等全控器件采用PWM或最近电平逼近调制不存在晶闸管换相过程。但VSC的成本比LCC高损耗也更大适合对电能质量要求高的场合。混合直流输电系统的常见方案是整流侧LCC、逆变侧VSC兼顾经济性和换相失败免疫能力这也是一种非常实用的折中方案。5.4 我推荐的仿真验证顺序一次完整的换相失败抑制研究大致可以这样安排先在详细模型上复现单相接地故障引起的换相失败记录未加任何改进措施时的γ、Ud、Id波形然后逐步加VDCOL优化、预测触发控制看γ的最低点被拉到什么水平再加调相机或者STATCOM观察换相失败是否从“发生”变成“不发生”或者从“连续失败”变成“单次失败后快速恢复”。6. 仿真验证的经验补充故障设置与信号观测里的几个坑最后这部分是我踩过坑后总结出来的实操要点直接影响到仿真结果的可信度。6.1 故障点设置与故障时刻的选择换相失败仿真最常见的问题是故障时刻选得不合适导致结果不具备代表性。交流系统电压是正弦波故障发生时刻决定了故障初相角也就决定了换相电压跌落的严重程度。我建议的做法是让故障起始时刻覆盖一个完整工频周期的多个相位点对比最严重的工况。Simulink里用Step模块配合三相故障模块Three-Phase Fault可以设定起始时刻但要注意故障电阻的过渡过程太小的故障电阻在仿真中容易引发数值振荡。接地故障电阻我习惯取0.01Ω到1Ω之间既贴近工程实际又不会让仿真发散。故障持续时间的设定也有讲究。很多研究用0.05秒相当于两个多工频周期足够引发换相失败并观察到后续恢复过程。如果研究连续换相失败把故障时间延长到0.1秒以上可以观察换相失败后的控制保护动作和二次失败过程。6.2 仿真步长与算法选择详细模型在阀级开关切换时会产生很陡的暂态分量变步长求解器很容易把步长压到非常小仿真跑得极慢。我的实际经验是详细模型使用固定步长步长取20到50微秒求解器选ode4经典四阶龙格库塔或者ode3平均值模型没有开关暂态用变步长ode15s即可仿真速度非常理想。如果你发现详细模型仿真速度慢到无法接受优先检查是不是离散化参数设置或者开关缓冲电路Snubber参数太敏感导致步长细化。缓冲电阻电容的参数值会影响数值稳定性这是Simulink LCC模型里最容易被忽略的地方。6.3 信号的观测与导出判断换相失败我习惯用最直接的方法在逆变侧每相阀电流波形上叠加γ信号波形观察阀电流波形中是否出现“不该出现”的电流脉冲。同步观察直流电压Ud波形和直流电流Id波形记录换相失败过程中的数值变化规律。数据导出的操作建议是使用“To Workspace”模块把γ、α、Ud、Id、交流母线电压幅值等关键信号保存到MATLAB工作区再写一个简单的后处理脚本直接计算γ的最小值、换相失败持续时间和直流功率跌落的深度。这样批量仿真时不用每次都盯着Scope看波形直接对比数组数据就能得出结论。6.4 平均值模型的校准问题平均值模型虽然快但也有一个隐藏的坑它是基于稳态方程计算Ud和Id的在故障瞬间的暂态响应需要通过外加惯性环节来模拟。如果不加惯性环节平均值模型对故障的响应比详细模型“乐观”很多γ的最低点可能偏大导致你低估换相失败风险。我的做法是在平均值模型的触发角输出和受控源之间加一阶惯性环节时间常数取2到5毫秒模拟晶闸管触发控制和换流过程的固有响应延迟。校准的方法是拿平均值模型跑同一个故障和详细模型对比γ波形。两者γ最低点偏差在2°以内说明平均值模型参数设置合理偏差太大就调整惯性时间常数和VDCOL响应速度。做完整套模型后我的体会是高压直流输电建模难不在Simulink操作而在对换流过程本身的理解。信号链路本质上就是交直流能量的变换与反馈控制把公式推清楚了、把信号流向理明白了后续加控制策略也好、改拓扑也罢都不会跑偏。