弱电网下LCL-VSC阻抗建模与稳定性分析：次超同步谐振与Nyquist验证

做并网逆变器稳定性分析这几年最让我头疼的问题之一就是明明控制器参数是按强电网调的一接到弱电网就出幺蛾子要么波形开始抖要么干脆谐振保护跳闸。后来被“阻抗分析”这套思路点醒之后才发现自己以前看稳定性的方式太片面了。最近刚好把“弱电网下LCL-VSC阻抗建模及稳定性分析次/超同步谐振Nyquist判据验证”这套流程在Simulink里完整跑通趁着印象还热把整个思路和实操过程整理出来给正在跟这个方向死磕的朋友一个参考。这套方法解决的是并网变流器的“环境适应”问题把变流器看成电网眼里的一个“阻抗”把电网也看成一个“阻抗”两者串联稳定性就转化成一个非常直观的判据问题。文章会覆盖阻抗模型的数学推导思路、dq轴小信号建模的要点、Nyquist判据的具体应用流程以及如何在Simulink环境里一步步搭建模型并验证理论结果。适合正在做新能源并网控制、微电网稳定性分析、或者被次同步/超同步振荡困扰的工程师和研究生参考。1. 为什么弱电网下LCL-VSC容易“犯病”1.1 弱电网的“脾气”短路比SCR才是关键先明确一个背景概念。电网的“强弱”在工程上通常用短路比Short Circuit RatioSCR来衡量定义为并网点短路容量与变流器额定容量之比。SCR大于5算强电网2到3之间是弱电网小于2就是极弱电网。SCR越低意味着电网等效阻抗主要是感性越大并网点电压越容易受功率波动影响。之前做过一个2MW的并网变流器模型电网短路容量10MVASCR刚好等于5调试时电流环参数随便整定都能稳定。后来把电网阻抗调大让SCR降到2.5左右同样一套参数逆变器输出电流开始出现明显的低频振荡频率大概在28Hz附近振荡这个频率既不属于基波也不是LCL固有谐振频率当时排查了很久才发现是阻抗交互导致的次同步振荡。这里要强调一个容易混淆的概念传统意义上的次同步谐振SSR最早是在汽轮发电机组串联补偿电容场景里发现的表现为低于同步转速的机械-电气耦合振荡。而变流器接入弱电网产生的次同步/超同步振荡本质上是变流器控制系统动态特性PLL、电流环、电压环与电网阻抗相互作用形成的电气振荡频率通常落在5Hz到基频之间次同步以及基频到几百赫兹之间超同步机械系统不参与这个过程。理解这个区别后续建模时才不会走偏。1.2 LCL滤波器的“双刃剑”效应LCL滤波器相比单L滤波器高频衰减能力更好同样的开关纹波抑制效果下总电感量可以更小成本更低这是它成为并网变流器标配的原因。但代价是引入了谐振峰而且这个谐振峰并不是固定不变的——电网阻抗变化会直接改变LCL的等效结构。从电路拓扑看LCL包含逆变器侧电感L1、滤波电容Cf、网侧电感L2。如果电网侧还有等效阻抗Lg那么L2和Lg实际上是串联的谐振频率会变成类似f_res 1/(2π × sqrt((L1 L2 Lg) × L1×Cf/(L1L2Lg)))的表达式Lg越大谐振峰越低。谐振峰一旦靠近控制带宽甚至控制带宽之内电流环的相位裕度就会被吃掉系统就开始振荡。这也是为什么很多工程师发现“参数明明没变换个并网点就不稳定”的根本原因。所以做弱电网下的稳定性分析本质上回答一个问题在电网阻抗从零变化到某一个临界值的过程中系统哪一刻会从稳定走向振荡这个临界点在哪里。阻抗分析法就是用来精确回答这个问题的手段。2. 阻抗建模把变流器和电网都“翻译”成阻抗2.1 阻抗分析的基本思路阻抗分析法的核心思想并不复杂。把变流器看成一个戴维南等效或诺顿等效电路输出阻抗为Zvsc(s)电网侧等效为一个理想电压源串联阻抗Zg(s)。两者在并网点连接那么并网点电流可以写成两个源分别作用叠加的形式。系统的稳定性取决于闭环传递函数的极点位置而这个闭环传递函数的分母是1 Zg(s)/Zvsc(s)或者等价的1 Yvsc(s)×Zg(s)其中Yvsc是变流器导纳。根据反馈控制理论如果Zvsc(s)本身稳定变流器空载或接理想电网时稳定Zg(s)也稳定电网本身无源那么整个互联系统的稳定性就由阻抗比Zg(s)/Zvsc(s)的Nyquist轨迹是否绕过(-1, j0)点来决定。这就是应用Nyquist判据的前提条件。这个思路最大的好处是“解耦”变流器内部再怎么复杂外部看到的只是一个阻抗特性电网也只是一个阻抗。两者之间的交互问题就变成了两个阻抗之间的匹配问题不需要把变流器内部的每个控制环都拉出来和电网阻抗联立求解。对于实际工程中“换个电网环境系统是否稳定”这类问题这种方法尤其高效。需要说明的是这个判断在单相系统或者三相对称系统中形式相对简单。对于三相VSC通常在dq同步旋转坐标系下建模得到的阻抗是一个2×2的矩阵也就是Zdd(s)、Zdq(s)、Zqd(s)、Zqq(s)四个元素。这时候阻抗比变成了一个矩阵Nyquist判据需要用广义版本比如回比矩阵的特征值轨迹或者用广义Nyquist判据处理。后面讲到实现细节时再展开。2.2 dq轴小信号建模的完整步骤对三相VSC做阻抗建模我习惯按下面四个步骤推进每个步骤都有容易忽略的坑。第一步是确定主电路拓扑和控制结构。常见配置是三相两电平VSC输出接LCL滤波器控制采用dq同步旋转坐标系下的双闭环电压外环或功率外环电流内环PLL锁相负责提供dq变换的相位参考。阻抗建模必须包含PLL的动态因为PLL的相位扰动会导致dq变换坐标系发生偏移这一项在小信号建模中会产生额外的耦合项是dq阻抗矩阵非对角元素的主要来源之一。第二步是建立时域小信号方程。在dq坐标系下电感电流和电容电压的状态方程需要线性化。具体操作是把电压电流变量写成稳态值加小扰动比如i_d I_d Î_d代入方程后去掉直流稳态项和高阶小量保留一阶扰动项。值得注意的是dq坐标系中存在转速项jωL的耦合这一项在扰动方程中会自然出现对应阻抗矩阵中的交叉耦合部分。第三步是建立控制系统的小信号模型。电流环PI调节器输出的是调制波经过PWM调制后在VSC交流侧产生平均电压平均模型下扰动量和控制量之间的关系需要从电流参考到实际电流的整个环路传递函数中提取。这里的关键是要把PLL的相位扰动Δθ纳入坐标变换矩阵的线性化中变换矩阵的小信号增量会带来额外的传递路径。第四步是整理成阻抗矩阵形式。把所有的扰动方程联立消去中间变量占空比扰动、PLL相位扰动等最终整理成[V̂_d; V̂_q] [Zdd(s), Zdq(s); Zqd(s), Zqq(s)] × [Î_d; Î_q]这个矩阵就是变流器在dq坐标系下的输出阻抗模型。推导过程手算非常繁琐我一般先用符号计算工具Matlab Symbolic或Maple把代数化简部分做了然后手工检查关键项的物理意义是否合理比如稳态工作点是否出现在系数中PLL带宽参数以什么形式出现在耦合项里。这些检查能有效发现推导过程中的低级错误。2.3 关键影响因素PLL、电流环、电网阻抗阻抗模型的表达式虽然复杂但几个关键参数的影响规律是清晰的掌握这些规律能让你在仿真结果异常时快速定位原因。PLL带宽对阻抗特性的影响最为显著。PLL带宽越高它对相位扰动的响应越快这会使得Zdq和Zqd这一对耦合项的幅值在低频段低于PLL带宽显著增大。直观的理解是PLL在努力跟踪电压相位而电压相位又受电流扰动影响这就相当于在d轴和q轴之间建立了一条额外的反馈通路。在弱电网下电网阻抗大并网点电压的相位受电流影响更大PLL的动态就被进一步激活次同步频段的负电阻特性经常由此产生。电流环带宽决定了阻抗在中频段通常100Hz到1kHz的形状。电流环PI参数中的比例系数直接贡献该频段的阻尼积分系数影响低频段增益。实践中有个经验电流环带宽一般设计在开关频率的1/10到1/5在这个范围里阻抗模型的形状相对可控过高容易放大高频噪声过低则动态响应太慢。电网阻抗的影响在阻抗图上表现为Zg(s)的幅值随频率增加而增加感性它与Zvsc(s)的交点频率决定了潜在振荡频率。两个阻抗幅值相等的频率点附近如果相位差接近180度那么就可能形成负阻尼振荡。这就是为什么“阻抗交点分析法”在工程上很流行——虽然严格的稳定性判据要看Nyquist轨迹但阻抗幅频曲线的交点频率和相位关系往往能给出非常直观的预警信号。3. Nyquist判据验证从理论到可执行的判断流程3.1 用阻抗比替代开环传递函数前面说过系统闭环传递函数的分母是1 Zg(s)/Zvsc(s)所以定义回比函数T(s) Zg(s)/Zvsc(s)它等价于一个开环传递函数系统的闭环稳定性就可以用Nyquist判据来判断T(s)的轨迹是否绕过(-1, j0)点。这里有一个工程上常用的简化路径先判断变流器本身接理想电网是否稳定再判断电网阻抗本身是否无源实际电网满足然后画T(s)的Nyquist图。如果T(s)在右半平面没有极点并且Nyquist曲线不包围(-1, j0)点那么互联系统稳定。如果曲线穿过或者包围该点则系统不稳定穿越时的频率对应振荡频率。对于三相VSC由于阻抗是2×2矩阵T(s) Zg(s) × [Zvsc(s)]^(-1)也是一个矩阵。有两种方式使用Nyquist判据一种是数值计算回比矩阵的特征值λ1(s)和λ2(s)分别对两个特征值应用Nyquist判据这就是广义Nyquist判据GNC另一种是如果耦合项较小比如Zdq和Zqd远小于Zdd和Zqq可以近似按SISO处理分别分析d轴和q轴通道的阻抗比。第一种严谨第二种直观实际工程中可以先做第二种快速排查再用第一种精确定判定。3.2 理论模型的Nyquist验证方法在Simulink里做仿真验证时我把整个流程分成了三步。第一步用解析阻抗模型画Nyquist图。把推导出的Zvsc(s)解析表达式和电网阻抗Zg表达式导入Matlab直接计算T(s)的频率响应用nyquist函数画图。重点观察曲线是否经过(-1, j0)以及穿越频率对应的数值。第二步做扫频仿真获取“实测”阻抗。在Simulink中搭建VSC的详细开关模型或者平均模型在并网点注入特定频率的小幅扰动电压在旋转坐标系下分别注入d轴和q轴扰动测量对应的电流响应用傅里叶分析提取该频率下电流响应中的同频分量计算阻抗矩阵在该频点的值。扫频范围从几Hz到几千Hz就可以得到一组离散频率点上的阻抗测量值。第三步对比理论模型和扫频结果。把扫频得到的阻抗矩阵画成Bode图叠加到理论模型的Bode图上。如果两者吻合良好说明建模推导正确如果偏差较大首先怀疑变流器控制模型的实现细节比如PLL结构、滤波器参数是否一致需要反复修正。实际做下来谐振频段附近扫频法的误差较大因为谐振处信号幅值变化剧烈且噪声影响明显。建议在谐振附近加密扫频点并且每个频点用不同的扰动幅值重复测量取平均能有效提升测量精度。4. Simulink仿真环境搭建全流程4.1 仿真模型的整体架构Simulink建模需要先规划顶层架构。我把模型分为四个模块电网模块三相电压源串联阻抗LgRg、LCL滤波器模块、VSC功率电路模块、控制器模块含PLL、电流环、坐标变换、SVPWM或SPWM发生。功率电路可以直接用Simulink Simscape/SimPowerSystems库里的Universal Bridge也可以自己用受控电压源搭平均模型。两种方式各有优劣开关模型更真实能看到谐波和开关纹波但仿真速度慢平均模型速度快适合扫频分析和控制参数调试但看不到开关频率附近的响应。我的做法是两种模型都搭调试参数用平均模型最终验证用开关模型。电网阻抗参数需要考虑弱电网场景通常Lg取值为0.5mH到5mH不等。以我最近做的380V/50Hz、额定功率50kW的样机模型为例L1取1.2mHCf取20uFL2取0.6mH电网阻抗Lg从0.5mH调到3mH对应SCR从大约6降到1.5左右此处按10kV侧短路容量计算会需要折算实际仿真中直接用电压源串联阻抗等效即可。4.2 关键控制环节的设置要点控制器是整个阻抗特性的“灵魂”搭建时需要格外注意以下几点。坐标变换必须统一基准。我习惯用余弦对齐d轴与电网电压矢量对齐的Park变换PLL输出角度用于所有dq变换。要特别注意Simulink中三角函数的角度单位弧度和初相设置否则很容易出现d/q轴分量明显不对的情况。电流环PI参数的设计需要一个基准工况。以L1L21.8mH额定电流峰值105A为例电流环带宽按800Hz设计取Kp (L1L2)×2π×800 ≈ 9.0Ki Kp × (R/L) 数量级实际计算可得约200。这些参数在不同功率等级下要重新计算不要直接套用别人的数值。PLL采用经典SRF-PLL结构PI参数按带宽50Hz设计。在弱电网下PLL带宽不宜太高否则会让阻抗模型在低频段表现出明显的负电阻特性。这里有个直观的经验PLL带宽每提高一倍Zdq耦合项的幅值在低频段大概增加6dB带来的相位偏移足以让Nyquist曲线右半部分变形。4.3 扫频法测量阻抗的具体实现扫频测量是整个验证环节中最容易出问题的一步这里给出一个我反复调试后确定可行的操作流程。首先确定扫频点序列。我通常用对数分布从5Hz到2000Hz取30到40个频点谐振频段根据理论模型预估附近加密比如额外增加10个点。每个频点的总仿真时长至少保证20个扰动周期以便后续FFT有足够的频率分辨率。其次是扰动注入方式。在并网点电压参考上叠加一个该频率的微小正弦扰动幅值一般为额定电压幅值的1%到2%。扰动太大容易激发非线性谐波太小则信噪比低。分别在d轴和q轴注入扰动量采集对应轴的电流响应然后利用FFT提取该频率分量的幅值和相位通过两个轴两两组合的计算关系求出阻抗矩阵的四个元素。实际计算时d轴注入可以求出Zdd和Zqdd轴电流响应除以d轴电压扰动q轴响应除以d轴扰动q轴注入同理。这个测量过程对仿真的数值精度要求比较高建议将求解器设为变步长相对容差改到1e-5。另外在测量扰动频率处的响应时要用同步解调而不是简单FFT简单FFT在场量下精度不够。同步解调的做法是把采集到的电流信号乘以该频率的正弦和余弦参考信号再用一阶低通滤波器取出直流分量得到幅值和相位的正交分解。这个方法在实测中也能用仿真里效果更有保证。5. 常见问题与排查技巧实录5.1 理论曲线和扫频实测对不上的典型原因我在仿真中碰到过三次理论模型和扫频结果对不上的情况排查过程很有代表性。第一次对不上低频段偏差大一看PLL模型原来是Simulink里PLL的输出角速度直接用连续积分器实现而我在理论推导时用了简化的一阶惯性近似。修正方法在理论推导中把PLL的传递函数替换成与Simulink实现完全一致的二阶表达式含低通滤波器项同时确保PLL线性化时叠加的稳态工作点转速是ω0而不是0。这类“实现细节不一致”导致的偏差最隐蔽排查时一定要逐模块对照理论假设和实际实现。第二次对不上超同步频段偏差大经过排查发现是控制器中数字延时的处理问题。Simulink中如果使用离散PI控制器PWM更新会有计算延时通常相当于一拍开关周期。而理论建模时如果忽略了这个延时高频段的阻抗相位会差很多。修正方法在建模时加入e^(-1.5Ts×s)的等效延时环节或者干脆把控制器的采样周期和PWM载波周期在参数设置里明确对应让仿真和理论两边都“说同一个口径”。第三次对不上问题出在扫频测量方法上。我在某些频点取的仿真时长不够FFT结果包含了启动暂态导致测量值抖动很大。后来把每个频点抛弃前5个扰动周期后再采样并增加重复次数取平均波形就平滑下来了。这个问题的解法在实测中同样适用建议优先排查。5.2 仿真不收敛和振荡误判的处理经验弱电网模型仿真容易出现数值不收敛的情况特别是开关模型在电网阻抗较小时更容易出现高频数值振荡。我最常用的手段包括在开关器件旁边并联一个小阻值RC缓冲电路通常是0.1Ω、1uF级别或者把LCL滤波器的电阻分量略微加大一点比如L1和L2的等效串联电阻从5mΩ调整为10mΩ这会在阻尼上没有实质影响却能大幅提升仿真收敛性。还有一个容易忽略的坑同时监控输出功率和并网点电压波形。如果已经判定系统不稳定仿真波形看起来可能非常“杂乱”不能简单地当成数值问题处理。我一般先切到平均模型做同样的测试如果平均模型也振荡说明是真实的控制失稳如果平均模型稳定而开关模型振荡才需要怀疑数值问题。总之先做线性化稳定性分析再用确定性仿真复现第三步才做非线性时域仿真验证这个顺序能大幅减少误判概率。5.3 参数快速调优的经验导则当Nyquist判据显示系统临界接近(-1, j0)时通常意味着阻抗幅频曲线在振荡频点附近的相位裕度不足。此时有几个方向可以调整。增加变流器等效阻尼是最直接的手段。在电流环中加入高通阻尼项有源阻尼或者在滤波电容支路串联电阻增大无源阻尼都能有效抬升阻抗相位。参数调整时先用阻抗扫频观察曲线变化找到峰值相移对应的频段再定向调整对应控制环参数。改变LCL谐振频率躲开电网阻抗谐振区也是常用手段。比如增大Cf会把谐振频率拉低这可能让阻抗交点频率从次同步段移到更安全的位置。不过需要注意改主电路参数往往是牵一发动全身会同时影响纹波抑制、损耗以及控制带宽限制要权衡考虑。对于弱电网应用最实用的还是采用自适应控制根据并网点电压谐波畸变率、短路比在线估算结果实时调节电流环带宽和PLL带宽。这个方向实现不算复杂但需要可靠性验证。实际调试时先把固定参数在某一个弱电网工况下调到稳定再加自适应外层能少走很多弯路。6. 实操心得与进阶扩展全套流程跑通之后我个人的体会是阻抗建模和Nyquist验证真正解决的不只是“这个系统稳不稳定”的问题而是让你对每个控制环节在频域中的角色有了具体的感知。以前调PI参数时是靠经验摸现在是先看阻抗曲线再调参数每一次修改都有明确的方向。模型搭好之后还建议做一个批量扫描实验固定控制参数网格化扫描电网阻抗Lg和电阻Rg记录每个工况下系统是否稳定把稳定边界画在阻抗平面上。这张图对工程交付非常有说服力直接标注出“安全运行区域”和“失稳风险区域”方便现场运行人员快速判断某个并网点是否适合接入。另外再分享一个后续很值得做的扩展方向把阻抗模型用于并网逆变器的并网适应性评估或者振荡抑制策略设计。如果某处现场已经出现了次同步振荡可以用实测阻抗曲线快速定位Zg与Zvsc的交点频率基于阻抗重塑思路设计陷波器或者附加阻尼控制。这套流程在风电变流器、光伏逆变器以及储能PCS的弱电网应用中都有很强的迁移性掌握了建模-判据-验证这条主线就能应付大多数并网稳定性分析任务。