前阵子整理旧文件,翻出两年前在MATLAB里搭的一套两端500kV LCC-HVDC输电模型,顺手跑了一遍稳态,直流电压、触发角、电流纹波都还挺争气。想当年搜“MATLAB 高压直流输电模型”,翻来覆去看到的要么是几十万瓦量级的教学demo,要么是PSCAD截图,真正拿Simulink把500kV等级的12脉波LCC-HVDC搭出来、还带完整控制策略的资料,少得可怜。
这篇文章就把我从零搭这套模型的完整思路捋一遍:主电路参数怎么定、6脉波桥怎么组成12脉波阀组、整流侧定电流和逆变侧定熄弧角是怎么配合的、解算器和启动时序怎么设,以及最后怎么判断仿真波形到底是“能跑”还是“可信”。适合两类人:一是准备用MATLAB做直流输电仿真的研究生,二是想快速建立一个HVDC基准模型的电力工程师。
1. 模型定位与主电路参数设计:先回答“搭给谁、参数怎么定”
很多人一上来就打开Simulink拖模块,拖完才发现变压器容量、平波电抗器、线路阻抗全是乱填的,仿出来的波形自己都不信。搭仿真模型的第一原则是:先算参数,再搭电路。尤其是LCC-HVDC这种强非线性系统,参数差10%,暂态特性能差出去一大截。
1.1 选择LCC-HVDC而不选VSC-HVDC的业务原因
LCC-HVDC用的是晶闸管,只能靠交流系统电压自然换相,没法自关断。它跟现在很火的MMC(模块化多电平换流器)思路完全不同。之所以直到今天还要研究LCC-HVDC,是因为在点对点大容量远距离输电里,晶闸管阀的成本、损耗、可靠性仍然有明显优势。我国已经投运的±500kV直流工程,绝大多数走的是LCC路线,两端各装一组12脉波换流器,这已经是教科书级的经典拓扑。
在MATLAB里摸这个模型,核心价值三个:验证控制策略、复现换相失败等典型故障、以及给大电网仿真提供等值接口。不要只把它当成一个“能出波形”的作业,后面做交流系统故障、多馈入直流、新能源基地送出研究时,这套基准模型都能当底层平台用。
1.2 基准参数演算:单极500kV、额定电流1200A
我建的是单极大地回线方式的极级模型,直流额定电压500kV,额定功率600MW。注意这里不谈双极,先把单极跑透,再复制成双极只是工作量问题。
额定直流电流计算很简单:
[ I_{dn} = \frac{P_{dn}}{U_{dn}} = \frac{600\text{MW}}{500\text{kV}} = 1.2\text{kA} ]
12脉波换流器由两个6脉波桥串联组成,每个6脉波桥承担约250kV的直流电压。6脉波桥的理想空载直流电压与换流变压器阀侧线电压有效值的关系是:
[ U_{d0(6p)} = 1.35 \times U_{vll} ]
考虑额定运行时整流侧触发角约15°,换相重叠角造成的压降约5%-7%,再加上直流线路电阻压降,阀侧线电压按200kV左右设计比较合理。这样每个6脉波桥的理想空载电压在280kV附近,留出了足够的调节裕度。
我给这套主电路定的核心参数如下表,可以直接作为初始化脚本的默认值:
| 参数项 | 取值 | 设计依据 |
|---|---|---|
| 直流额定电压 | 500kV | 目标等级 |
| 直流额定电流 | 1200A | 600MW/500kV |
| 换流变压器阀侧线电压 | 200kV | 满足6脉波桥空载电压要求 |
| 换流变压器容量 | 每端约620MVA | 直流功率+无功损耗余量 |
| 短路比SCR | 3.0 | 中等强度交流系统,贴近实际 |
| 直流线路长度 | 200km | 中等距离集中参数建模 |
| 平波电抗器 | 0.6H每极 | 抑制直流电流冲击与谐波 |
| 逆变侧熄弧角 | 17° | 兼顾换相裕度与无功消耗 |
交流系统采用戴维南等值,等效电压源加串联阻抗。短路比3对应的系统阻抗标幺值就是1/3约0.333,折算到500kV交流母线侧后的实际阻抗值要按基准容量和电压换算,这些都在初始化脚本里统一算,不要在Simulink界面里一个个手敲。
1.3 滤波器与无功补偿的初步估算
LCC换流器在运行中会吸收大量无功,经验上约为直流功率的50%-60%。600MW对应需要补偿约300-360Mvar。交流侧配置了11/13次双调谐滤波器、若干高通滤波器,它们既滤特征谐波又提供无功。直流侧因为12脉波的特征谐波是12k次(k=1,2,...),所以直流滤波器主要谐振在720Hz和1440Hz附近。
滤波器参数不需要一开始就精细优化,先按“谐振点对准特征谐波频率、品质因数取30-50”的经验值给一组初值,后面看谐波频谱再微调。重点是别把滤波器直接并联在换流变母线上而不加任何组断路器,那会在启动时造成很大的涌流。
2. 一次系统Simulink搭接:从6脉波桥到12脉波阀组的模块化实现
主电路的搭建顺序很重要,我建议先搭一个6脉波桥,验证触发脉冲相位对了,再升级成12脉波。不要一上来就拖12个桥臂,出错时根本查不过来。
2.1 换流变压器与Universal Bridge阀组建模
换流变压器在Simscape Electrical里用三相双绕组变压器模块,关键设置是漏抗标幺值放在0.1-0.2之间,太低了换相压降不足、谐波变大;太高了直流电压压降过大,还会增大换相失败风险。我们模型中取0.15pu。
两个6脉波桥要实现12脉波效果,关键在于两桥的换流变压器阀侧线电压相位必须错开30°。一般做法是一台变压器用Y-Y接线,另一台用Y-Δ接线。在Simulink里直接改变压器绕组的连接方式即可。接线组别选Yy0和Yd1,这样两个桥的交流侧电压错开30°电角度,直流侧串联后输出电压就是12脉波。
阀组用Universal Bridge模块,Power Electronics Device选Thyristor。晶闸管模型的Ron建议设在0.001Ω量级,Snubber电阻1kΩ、电容0.1μF这个量级起步。Snubber参数太大会引入额外时间常数,导致仿真步长被压得很小,这是性能杀手。
2.2 直流线路、平波电抗器与滤波器摆放
直流线路我直接用集中参数π型模块。200km线路按单位长度参数R=0.012Ω/km、L=0.8mH/km、C=0.013μF/km折算,综合得到串联电阻约2.4Ω,电感约160mH,对地电容约2.6μF。用集中参数模型在工频故障研究里精度够用,但如果要研究行波保护、高次谐波沿线路传播,那必须上分布参数线路模型,后者计算量明显增大,启动阶段也更难收敛。
平波电抗器接在换流器直流出口和直流线路之间,每极0.6H。这个电感的作用一是抑制直流侧谐波电流,二是在换相失败或直流短路时限制电流上升速率。0.6H配1200A电流,电感上存储的能量不小,所以启动和故障切除时要注意过电压问题。
直流滤波器接在平波电抗器靠线路那一侧,避免滤波器与平波电抗器形成不希望的谐振回路。交流滤波器接在换流变压器网侧母线上。
2.3 子系统和初始化脚本的封装思路
不要把所有模块平铺在一个Simulink窗口里,那会让你在查线时崩溃。我把系统拆成了四个子系统:整流站、逆变站、直流线路、控制与测量。整流站里又分换流变压器模块、阀组模块、交流滤波器子系统。每个子系统留出Input/Output端口,端口命名用交流电压、三相电流、触发脉冲、直流电压等明确标签。
所有电气参数不直接写死在模块参数里,而是在一个init_hvdc.m脚本里用变量赋值,然后让Simulink模块引用这些变量。下面是脚本的核心片段:
%% 直流主回路参数 Udn = 500e3; % 额定直流电压, V Pdn = 600e6; % 额定功率, W Idn = Pdn/Udn; % 额定直流电流, A %% 换流变压器 Uvll = 200e3; % 阀侧线电压有效值, V Xt_pu = 0.15; % 漏抗标幺值 S_trafo = 620e6; % 换流变压器容量, VA %% 直流线路 L_line = 200; % 线路长度, km R_line = 0.012*L_line;% 串联电阻, ohm L_lineH = 0.8e-3*L_line; % 串联电感, H C_line = 0.013e-6*L_line; % 对地电容, F %% 平波电抗器 L_smooth = 0.6; % H %% 控制参数 alpha_ref = 15*pi/180; % 整流侧触发角额定值, rad gamma_ref = 17*pi/180; % 逆变侧熄弧角参考值, rad Kp_cc = 0.5; % 定电流控制器比例系数 Ti_cc = 0.02; % 定电流控制器积分时间常数, s这套做法的好处是,参数调整后只要重新运行脚本,整个模型就同步更新,不会出现“变压器容量改了但控制系统基准值没改”这种低级的参数失配问题。
3. 控制系统策略与实现:电流裕度法里到底谁说了算
LCC-HVDC控制系统是整个模型的灵魂。很多初学者在Simulink里搭好了主电路,但控制部分只是简单接了个常数触发角,结果就是模型完全无法应对扰动。要做出一个可信的模型,必须把整流侧定电流、逆变侧定电压/定熄弧角以及低压限流环节完整搭起来。
3.1 整流侧定电流控制
整流侧正常工况下通过调节触发角α来维持直流电流等于参考值。当直流电流偏大时,增大α降低直流电压从而减小电流;电流偏小时则减小α。PI控制器输出就是触发角指令,限幅一般在5°到150°之间。5°最小触发角是为了保证晶闸管有足够的导通时间,150°最大触发角是为了在故障时快速把电流压下来。
电流测量的滤波时间常数很关键,一般取几毫秒到几十毫秒。滤波时间常数太小,测量噪声会放大到触发角里;太大,控制响应变慢,故障时来不及动作。
PI参数我没有用复杂的整定公式,而是从Kp=0.5、Ti=0.02s起步,做直流电流阶跃响应,观察超调量和调节时间,再微调。LCC系统的滞后较大,积分时间常数不能取得太小,否则系统会振荡。这里有一个通用原则:先在纯电流环里把PI调到阶跃超调小于10%、调节时间50ms左右,再接入其他控制环节。
3.2 逆变侧定熄弧角控制与电流裕度
逆变侧运行在整流状态以外,触发角α一般在90°到150°之间,对应熄弧角γ约为15°到18°。定熄弧角控制的目标是维持γ在参考值附近,保证晶闸管退出导通后有足够的恢复阻断能力。γ太小容易换相失败,γ太大无功消耗增加、运行经济性变差。
电流裕度法是两端配合的基本逻辑系统。整流侧电流参考 (I_{ref_rect}) 和逆变侧电流参考 (I_{ref_inv}) 之间保持一个裕度,通常 (I_{ref_inv} = I_{ref_rect} - 0.1 \times I_{dn})。正常时整流侧定电流,逆变侧定γ;当交流电压下降导致整流侧无法维持电流时,直流电流会下滑到逆变侧电流参考以下,此时逆变侧自动切换为定电流控制,防止电流失控。这个切换是控制系统自然发生的,不需要逻辑判断。
实现的时候,把逆变侧电流控制器与定γ控制器都算出来,选两者中触发角更大的那个作为最终触发角指令。这个“取大”逻辑非常关键,很多模型控制紊乱都是因为这里直接相加或者用了错误的切换逻辑。
3.3 VDCOL低压限流环节
VDCOL是LCC控制系统不可省略的一环。它的功能是:直流电压跌落到一定水平后,自动降低两端的电流参考值,避免故障期间持续大电流导致交流电压进一步坍塌,同时为故障切除后电流快速恢复创造条件。
我的VDCOL参数是:直流电压标幺值高于0.7时电流参考为额定值;低于0.4时电流参考限制在0.3pu;中间线性过渡。这个参数在不同文献里略有差异,但逻辑都差不多。在Simulink里用一个二维查表模块就可以实现,输入是直流电压测量值,输出是电流参考的倍率。
3.4 PLL与触发脉冲生成
触发脉冲生成必须基于准确的换相电压相位。实测上我直接对换流变压器阀侧线电压做锁相环,锁相环输出角度与阀侧电压相位对齐后,再叠加触发角延迟。关键坑在于:变压器接线组别不同,阀侧电压相位与网侧不一致,PLL锁错了相位,脉冲就会乱。处理办法简单粗暴:先用示波器看阀侧三相电压,手动对比PLL输出角度是否在对应相电压过零点附近,不对就加固定的相位补偿角。
轻量建模时也可以使用Simscape Electrical自带的六脉冲触发模块,但它的封装逻辑在故障场景下不灵活。我最终还是选择自己写一个触发脉冲生成子系统,用PLL角度加上α指令,经过比较和脉冲分配器生成6路脉冲,再接Universal Bridge的门极输入。这样每个触发角的变化都能实时反映在波形里。
4. 仿真设置与启动过程:解算器、解锁时序和初始化脚本
主电路和控制模型都搭好之后,真正折磨人的是仿真设置。LCC-HVDC是全电力电子系统,晶闸管开关在极短时间里频繁动作,Simulink默认的求解器往往跑不动或者跑出数值泛滥的结果。
4.1 解算器选择与容差设置
我强烈建议用变步长ode23tb解算器。这个求解器在处理刚性系统时比ode45稳定很多,这是电力电子仿真老油条们的共识。相对容差设1e-3,绝对容差设1e-4,够用就好,不要为了“心里踏实”调到1e-6,那会让仿真时间成倍增加。
电力电子器件本身是理想开关模型,存在高频的开关暂态,Simulink为了捕捉每个开关事件会不断压缩步长。如果发现仿真推进速度像蜗牛爬,优先检查两件事:一是Snubber参数是否过大,二是是否误开了连续模式下不必要的高精度器件模型。
4.2 解锁时序:不要一上来就给触发脉冲
LCC换流器启动不能直接把触发角设成15°然后全压启动,那样必炸。我的做法是:0到0.1s内所有晶闸管闭锁,交流系统给换流变压器充电;0.1s后整流侧解锁,但触发角初始设为120°,让直流电压从零缓慢建立;然后通过斜坡方式把触发角从120°过渡到60°以内,等直流电流建立起来后,定电流控制自动接管,触发角平滑过渡到15°左右。逆变侧要等整流侧把直流电压抬升到一定值后再解锁,避免倒灌电流。
这个解锁时序用Step模块加Rate Limiter就能实现,不需要编程。但时序的先后很讲究。我最开始就是整流侧、逆变侧同时解锁,结果逆变侧还没建立反压就收到触发脉冲,晶闸管直接误导通,仿真直接发散。
4.3 控制器的初始值设定
PI控制器的初始输出不要设0,而应该设成对应的稳态触发角,比如整流侧PI输出初始化为约15°对应的弧度值。这样启动瞬间控制系统不会产生很大的反向调节冲击。更重要的是积分器初始值也要对应设定,否则即使输出限幅正确,积分器内部还是会积累很大的偏差。
数据采集方面,建议在整流侧和逆变侧各放一个Scope记录直流电压、直流电流、触发角指令、熄弧角测量值和三相交流电压电流。Simulink的Scope数据建议通过“Log signals”方式记录到工作区,方便用代码做FFT和指标计算。
5. 稳态、动态和故障波形判读:什么样的曲线才算模型可信
模型能跑起来不等于模型是对的。很多人的LCC模型看起来波形正常,但仔细观察会发现电压电流数值与设计值差得离谱。这里给出我用的几个核心判据。
5.1 稳态波形的基本判据
仿真3s左右,系统进入稳态后应看到:
- 整流侧直流电压约500kV,逆变侧直流电压略低,两者之差等于直流线路压降;
- 直流电流稳定在1.2kA附近,纹波在几个百分点内;
- 整流侧触发角稳定在15°附近,逆变侧熄弧角稳定在17°附近;
- 交流侧电压电流呈明显的12脉波特征,即波形叠加有阶梯状谐波,11/13次谐波经滤波后明显削弱;
- 交流电流总谐波畸变率低于5%,但这里取决于滤波器参数,初调阶段不需要追求极致。
我习惯做一个功率平衡校验:整流侧交流输入有功应当等于直流功率加整流站损耗,逆变侧输出有功等于直流功率减逆变站损耗。两端有功差额应该在可控范围内,如果发现直流功率比交流侧输入还高,那一定是测量方向或符号错误。
5.2 直流电流阶跃动态响应的意义
在稳态基础上,把整流侧电流参考从1.0pu阶跃到1.1pu,观察直流电压和电流的响应。合格的响应是:电流在几十毫秒内平滑上升,略微超调后稳定在1.1pu,直流电压短暂下降后恢复。如果电流响应振荡剧烈,首先怀疑PI积分时间常数太小;如果上升极慢,则可能是线路电感和平波电抗器把电流环带宽压得很低,需要调整控制增益。
5.3 逆变侧交流单相接地故障与换相失败
这套模型最值得做的暂态测试就是逆变侧交流母线单相接地故障。故障发生后,换相电压畸变,逆变侧熄弧角裕度消失,会出现换相失败。观察到的现象是:直流电压骤降、直流电流骤升、整流侧定电流控制开始增大触发角,但受电流上升速率限制;随后VDCOL介入把电流参考压低,故障切除后电流恢复。整个过程应该在几百毫秒内结束。
这个故障测试对模型的真实性是很好的验证。一个可信的LCC-HVDC模型必然能复现换相失败,并且能通过VDCOL和电流裕度法实现可控恢复。如果模型在故障发生时直接发散,或者恢复过程极其混乱,那大概率是控制系统逻辑或者参数有问题。
5.4 用FFT辅助判断滤波器效果
Simulink的FFT工具可以对记录的交流电流和直流电流做频谱分析。交流侧应看到11次(550Hz)、13次(650Hz)被滤波器大幅度削弱;直流侧的特征谐波主要集中在720Hz和1440Hz附近。如果频谱中出现了大量的非特征谐波,比如3次、5次,那往往是变压器励磁涌流或者控制不对称引起的,需要检查解锁过程是否对称、PLL相位是否准确。
6. 调模型踩过的坑与提速优化技巧
这部分才是真正省时间的部分。以下问题我全都实际遇到过,有些甚至花了几天时间排查,写在这里帮你绕开。
6.1 仿真不收敛的第一检查顺序
仿真报错“Solver unable to resolve”时,不要慌,按这个顺序查:先查解锁时序,再看Snubber参数,接着检查PLL相位,最后看PI控制器积分溢出。我遇到最多的是第二种——为了省事把Snubber电容设成10μF,结果每个开关事件都引起长时间数值振荡,步长被压到纳秒级,发散只是时间问题。把Snubber电容降回0.1μF,问题立刻消失。
6.2 触发脉冲相位错乱
这个问题非常隐蔽。12脉波阀组中两个换流变压器接线组别不同,它们的阀侧电压相位差正好是30°。触发脉冲生成系统如果只锁了一个桥的电压相位,另一个桥的脉冲就必须额外补偿30°,否则两个6脉波桥的触发时序对不上,直流电压波形会出现严重的非对称凹陷。查这个问题的方法是:分别看两个桥的直流电流,理论上应该完全相等;如果不等,先查触发相位补偿。
6.3 启动阶段直流电压过冲
直流线路电容、滤波器电容在启动时会充电,如果触发角下降太快,产生的电压过冲可能会击穿模型里的理想开关,表现为运行到某时刻突然发散。解决办法是给触发角斜坡加长斜坡时间,比如从闭锁到稳态用0.3s,而不是0.05s。仿真慢一点,但启动过程更真实。
6.4 模型跑太慢的终极手段:均值模型
如果只是研究控制逻辑或者大电网稳定性,不需要关心阀级暂态波形,建议把12脉波换流器替换成均值模型。均值模型把换流器等效为受控电压源和受控电流源,用代数方程描述交流侧与直流侧的关系,不再模拟每个晶闸管的开关过程。Simulink里可以用受控电压源和受控电流源手动搭建,控制系统完全复用。替换后仿真速度能提高5-10倍,有些场景甚至更多。
我自己的习惯是:先搭详细开关模型,用来校核稳态和典型故障镜像;校核通过后另存一份均值模型版本,用在做长时间动态仿真和批量参数扫描时。两份模型的控制系统完全共用,只是主电路不同,保证结果可以互相对标。
6.5 波形看似对但总觉得哪里不对
有一次我调完滤波器,怎么看交流电流波形都整齐,但直流电压总有约3%的纹波。FFT一查,直流侧不仅存在12次谐波,还出现了较强的6次谐波。后来找到原因:两组换流变压器的漏抗标幺值不一致,导致两个6脉波桥的换相压降不同,直流侧6次谐波无法抵消。把两组变压器漏抗改成相同标幺值后,纹波降到一个健康水平。这类问题只有靠FFT才能暴露,肉眼波形很难发现。
还有一次,逆变侧熄弧角读数一直在10°和20°之间剧烈摆动,但触发角指令很平稳。排查了半天,发现是熄弧角测量环节用了未经滤波的瞬时电压信号,测量噪声被PI控制器放大。给测量信号加一个5ms的滤波环节后,熄弧角控制就稳定了。LCC-HVDC里测量回路的滤波和延时建模,对控制系统仿真的真实性影响极大,不要省略。
这套模型后来我经常拿来当基准,换交流系统阻抗、改控制参数、接多馈入直流网络,都从这套500kV基础模型开始改。如果你也要搭,我真心建议先跑通一个6脉波桥,确认触发脉冲、解锁时序、定电流闭环都已经正常,再去折腾12脉波。直接上12脉波不是不行,只是出了问题不好定位,而在电力电子仿真里,问题定位往往比搭模型本身更耗时间。