弱电网下LCL-VSC次/超同步谐振的阻抗建模与Nyquist判据分析

前段时间做并网逆变器稳定性分析时，我又撞上了那个绕不过去的组合：弱电网下面，带LCL滤波器的VSC系统，在次同步和超同步频段冒出谐振隐患。用阻抗建模把系统拆开，再用Nyquist判据验一遍，稳定裕度不足的问题就清清楚楚摆在眼前。这类问题在新能源并网领域几乎天天见，但很多人卡在不知道阻抗模型怎么建、仿真怎么搭、稳定裕度怎么定量算。

这篇围绕“弱电网+LCL-VSC+次/超同步谐振+Nyquist判据+Simulink仿真”这条主线，把我做阻抗建模和稳定性分析的完整思路、参数依据和踩坑经验整理出来。适合正在做并网变流器稳定性课题的研究生，也适合做光伏、储能、风电变流器工程的工程师。无论你手里是做仿真验证还是现场分析，这套流程都能直接拿过来参考。

1. 弱电网环境下LCL-VSC系统的“难”与“病”

1.1 弱电网到底弱在哪

工程上描述电网强度最喜欢用的就是短路比SCR。SCR的含义很简单：公共连接点（PCC）处的短路容量与变流器额定容量之比。SCR在20以上可以视为极强电网，3到10属于中等强度，低于3就进入弱电网区间，低于2基本是极弱电网。弱电网的核心特征是电网等效阻抗（以电感分量为主）不再能忽略，而这个阻抗还会随着运行方式改变——远端线路投切、邻近负荷接入、变压器分接头变化，都会让Lg在很大范围内变化。

写仿真参数之前，我习惯先给一个具体量级的计算。以380V、50kW的三相并网VSC为例，额定电流大概76A。当电网电感Lg取5mH时，工频50Hz下的感抗Xg=2π×50×0.005≈1.57Ω，PCC处的三相短路电流约为380/(√3×1.57)≈140A，对应短路容量约92kVA，SCR≈1.84。这已经是典型的弱电网工况。同一个平台，Lg取2mH时，SCR约4.6，算是弱电网边缘。理解了这层对应关系，你设置仿真参数时才有底气，不会随手填个0.1mH却自称弱电网。

弱电网给控制带来的本质变化，是PCC点的电压不再由电网单方面决定，而是由“电网电压源串联Lg”与“变流器注入电流”共同决定。注入电流在Lg上产生压降，这个压降又会反馈到PLL的输入和电流控制的参考坐标系里。系统内部多了一条看不见的反馈通道，原本在强电网下可以忽略的耦合效应，在弱电网下会被放大成稳定性问题。

1.2 次同步与超同步谐振从哪冒出来

先给谐振一个直观图像。次同步谐振，指的是振荡频率低于工频50Hz，工程上常见的是5到45Hz这个区间；超同步谐振，指的是高于工频的频率段，通常讨论55到150Hz，严重时能到几百Hz。名字听起来复杂，其实在示波器上非常好认：有功功率和PCC电压出现低频包络波动，波形像被一个慢速信号调幅。

谐振的根源，在于VSC控制系统在特定频段呈现出“负电阻”特性。普通电阻是消耗能量的，负电阻等于向外注入能量。当电网阻抗与这个负电阻在某个频率点相遇，整个回路在该频率下的总阻尼为负，扰动不但不会衰减，反而会越涨越大，直到触发保护或者波形严重畸变。

LCL滤波器在这里起了推波助澜的作用。LCL本身在谐振频率处有一个固有谐振峰，典型设计把这个峰放在开关频率的1/10到1/5附近。弱电网的Lg串进来之后，等效谐振点会向低频移动，整个系统的阻尼特性进一步恶化。这就不难理解为什么“弱电网+LCL+VSC”这三个词凑在一起时，稳定性问题会格外突出，因为每个环节都在给次/超同步频段的谐振创造条件。

2. 阻抗建模：把吃力的时域分析变成频域一条曲线

2.1 建模核心思路：变流器也是“阻抗”

阻抗建模的基本思想，是把整个并网变流器系统看成一个大“阻抗Zinv(s)”，电网看成“电压源串联Zg(s)”。PCC处电压和注入电流之间的关系，通过变流器的阻抗体现。这样一来，时域里复杂的非线性控制问题，可以转换到频域里，用一条幅频相频曲线来表征。

这个思路之所以好使，是因为它把一个“变流器+滤波+控制”的闭环系统等价成了无源元件网络。得到的阻抗曲线可以直接和电网阻抗做比值分析，借用奈奎斯特判据或者Bode图判断闭环稳定性。相比在时域里观察波形是否发散，这种方法能定量告诉你哪个频段的稳定裕度最薄弱，也能指导参数调整方向。

建模时要注意一个容易混淆的点：VSC的阻抗并不是单纯把LCL滤波器参数代入拉普拉斯变换那么简单，它必须包含控制环路的影响。电流环PI参数、PLL带宽、甚至功率外环的带宽，都会以频率函数的形式出现在阻抗表达式里。我见过不少初学者绕组组时只搭了LCL加一个理想电流源，然后抱怨阻抗曲线和实测对不上——因为控制环节的贡献全被漏掉了。

2.2 LCL滤波器和VSC控制环节的阻抗特性

LCL滤波器作为无源网络，它的阻抗特性是明确的。逆变器侧电感L1、滤波电容Cf、网侧电感L2组成一个三阶网络，在谐振频率fres处出现一个阻抗低谷（或者说谐振峰）。fres的计算公式是：

fres = (1/2π) × √[(L1+L2)/(L1×L2×Cf)]

举个例子：L1=0.5mH，L2=0.15mH，Cf=30μF，则fres≈2.7kHz。这个频率远高于基波，但远低于开关频率（10kHz），是LCL的典型设计区间。如果这个谐振峰不加阻尼，系统在这个频段附近的相位会急剧变化，对稳定性非常不利。

VSC控制环节的阻抗贡献要复杂得多。电流环在dq坐标系下实现对电流的快速跟踪，它本身的带宽决定了变流器呈现的“电感等效值”会随频率变化。PLL更是重头戏，它在弱电网下引入的相角扰动会让dq轴分量之间产生耦合。数学上表现为阻抗矩阵的非对角元素不为零——也就是说，d轴扰动会在q轴产生响应，反之亦然。这种耦合是次/超同步频段负电阻现象的主要来源之一。

建模时常用的框架是dq阻抗矩阵，写成：

[Vd] [Zdd Zdq] [Id] [Vq] = [Zqd Zqq] × [Iq]

其中Zdd、Zqq是自阻抗，Zdq、Zqd是耦合阻抗。稳定性分析时不能只看对角项，耦合项在弱电网下往往才是决定稳定性的关键。若想更精细化，还可以把阻抗表达成正负序的形式，直接对应次同步和超同步频段的物理含义，方便观察各频段的阻尼水平。

2.3 VSC阻抗分布解析式的大致构成

如果要用解析式表达VSC阻抗，可以按“电机侧”和“控制侧”两部分来拆。电机侧包含L1、L2、Cf的串联-并联组合，控制侧包含电流环传递函数、PLL传递函数以及由它们带来的附加阻抗项。

一个简化的VSC序阻抗表达式可以写成：

Zinv(s) = ZLCL(s) + Gc(s) + Gpll(s) × H_coupling(s)

其中Gc(s)是电流环的等效阻抗贡献，Gpll(s)是PLL的贡献，H_coupling(s)则体现了耦合通道。具体表达式在不同坐标系下形式差异较大，但结构可以这样理解：无源滤波器提供基础阻抗，控制环节在基频附近会修正这个阻抗甚至改变其符号。弱电网下，PLL在低频段的响应会让修正项的虚部变得很危险。

实际工程中，我更推荐的做法是先用数学推导写出阻抗解析式，再用Simulink扫频结果校验。两边的曲线如果能重合，说明模型正确；如果不重合，优先检查PLL的线性化是否做对、电流环在dq坐标系下的交叉解耦是否考虑完整。解析式最大的意义不是拿来手算，而是让你在做参数扫描时有明确的目标函数。

3. 次同步与超同步谐振的机理和参数敏感性

3.1 负电阻是怎么产生的

负电阻这个说法听起来违反直觉，但在含控制器的变流器系统中确实是真实存在的。它来源于控制环路的“超前响应”效应。当一个频率较低的扰动电流流过VSC时，如果PLL和电流环对这个扰动的响应延迟足够小，VSC对外表现的电压-电流关系就会呈现类似“电流增加时电压也增加”的异相特征——这正是负电阻的定义。

再具体一点：当PCC电压出现次同步扰动时，PLL检测到相位变化，会调整旋转坐标系的角度。坐标系的偏移会改变电流参考值的方向，控制器随之改变输出电压。如果这个电压变化和扰动电流的变化方向一致，VSC就在向系统注入能量。能量注入意味着振荡的阻尼为负，在频域上对应阻抗实部为负。

这个现象在SCR越低的电网中越明显。因为弱电网下，同样的扰动电流会在Lg上产生更大的电压波动，PLL感知到的相位扰动更大，控制器的负阻尼效应也就更突出。这就是为什么强电网下几乎见不到次同步振荡，而弱电网下却很常见。

3.2 PLL和LCL参数的影响

PLL带宽是影响次同步谐振的第一敏感参数。PLL带宽取得高，锁相速度快，但在弱电网下会把电网阻抗的相位扰动无放大引入控制环，负阻尼加剧。降低PLL带宽，负电阻的频段范围缩窄，次同步振荡的阻尼得到改善。代价是动态响应变慢，在电网电压跌落时的暂态性能会变差。工程上弱电网工况的PLL带宽通常控制在5到20Hz，具体取值要看整体稳定裕度仿真结果。

电流环PI参数的影响同样显著。电流环带宽升高，变流器对扰动的响应更快，但是在次同步频段可能带来相位滞后不足的问题，反而降低该频段的阻尼。适当的做法是让电流环带宽与PLL带宽拉开距离，避免两个控制环在目标频段发生交互。

LCL的参数影响主要体现在高频段。谐振频率设计得离控制带宽太近，会在Nyquist曲线的高频段形成一个向（-1,j0）点逼近的环，削弱幅值裕度。这种情况下通常需要给LCL加阻尼，有源阻尼和无源阻尼各有取舍。这里我不展开对比，但可以给一个经验值：LCL谐振频率最好落在开关频率的1/10到1/6，且远离电流环带宽3倍以上。

3.3 次同步与超同步的对照

把次同步和超同步放在一起看，会更有整体感。两者本质上都源于VSC控制环路与电网阻抗的频率耦合，区别在于作用路径不同。

表格里列出的频率范围不是固定的，实际系统可能因为参数不同偏移很大。做研究时不要死记频段，要看阻抗曲线在哪一段变差。我做次同步分析时，关注的是1到50Hz这个区间；做超同步时，关注的是LCL谐振频率附近以及开关频率之外可能的高频耦合。两种谐振在Simulink里都能通过注入扫频直观看到，后文会讲具体方法。

4. Nyquist判据验证稳定性的完整流程

4.1 阻抗比判据的原理

用一个前向通道的观点来看并网系统：电网阻抗Zg(s)和变流器阻抗Zinv(s)串联构成闭环。闭环特征方程可以等价为1+Zg(s)/Zinv(s)=0，因此开环传递函数可以取L(s)=Zg(s)/Zinv(s)。这时候奈奎斯特判据可以直接使用：只要L(s)的奈奎斯特曲线不包围（-1,j0）点，系统闭环稳定。

在MATLAB里做这件事非常简单。你只需要在命令行里完成下面几步：

% 假设已经通过扫频得到频率向量 freq、逆变器阻抗频率响应 Zinv 和电网阻抗 Zg % 计算开环传递函数 L = Zg / Zinv L = Zg ./ Zinv; % 绘制奈奎斯特曲线 figure; nyquist(L); grid on;

Nyquist图既可以基于解析传递函数，也可以基于扫频数据。基于扫频数据的优势是它不需要完整的状态空间模型，直接从仿真结果就能评估稳定性。这在工程上非常实用，尤其是当你没有精确的解析表达式时。

判稳的时候有一点要提醒：奈奎斯特曲线显示的是频率从-∞到+∞的完整映射。有时我们只看正频率部分，容易被误导。尤其是VSC系统在dq坐标系下的阻抗曲线并不关于实轴对称，负频率部分对判稳结果很有影响，所以务必保持MATLAB默认的完整绘制范围。

4.2 从Nyquist曲线上读出稳定裕度

Nyquist曲线能告诉你的不只是“稳定/不稳定”这个结论，它更是一个诊断工具。曲线距离（-1,j0）点越近，说明闭环特征根越接近虚轴，系统稳定裕度越小。工程上通常要求奈奎斯特曲线距离（-1,j0）点有足够的增益裕度和相位裕度。增益裕度一般要求3dB以上，相位裕度一般要求30°以上。

实操中我习惯把Nyquist曲线按频段分颜色展示。例如次同步频段用实线、超同步频段用虚线，这样一旦曲线逼近（-1,j0）点，立刻能看出是哪个频段在“惹事”，进而回溯到对应的控制参数。这个思维路径比单纯调参碰运气要高效得多。

如果你做的是参数优化，可以把Nyquist曲线随参数变化的趋势叠加在一张图上。比如把Lg从强电网数值逐渐增加到极弱电网数值，看奈奎斯特曲线如何一步步靠近（-1,j0）点。这个渐变过程能直观地告诉你：系统是在哪个电网强度下失稳的，主导因素是哪一个频段，为后续加阻尼措施提供依据。

5. Simulink仿真建模与阻抗扫描验证

5.1 搭建主电路与控制系统

Simulink仿真模型的核心结构分为四块：三相电网与弱电网阻抗、LCL滤波器、VSC主电路、控制系统。电网用Three-Phase Source加串联R-L来实现弱电网效果，Lg按前面计算的参数设置。LCL滤波器直接使用三个串联/并联的RLC支路搭建，注意L1、L2的寄生电阻不要设为零，否则仿真器容易出现数值奇异，而且免费寄生参数会让理论分析产生虚假的谐振峰。

VSC主电路使用Universal Bridge或者Three-Level Converter，根据你研究的对象选择两电平还是三电平拓扑。控制系统的结构从外到内是：功率参考（P/Q）→ 电流环 → PWM调制。PLL从PCC电压提取相位，用于坐标变换。这部分在Simulink里实现时，最容易被忽略的是坐标变换的初始角度对齐。如果初始相位设置不对，启动瞬间会有一个很大的冲击电流，叠加上弱电网的大Lg，可能导致缓启动机制误动作。

一个比较稳妥的做法是先让PLL在开环条件下锁定一段时间，锁稳后再切入控制系统；或者给电网电压一个软启动斜坡，让控制器在电压逐渐建立的过程中完成同步。仿真步长方面，如果开关频率10kHz，PWM生成建议用连续或定步长为1e-6的仿真器，既保证精度又能接受仿真时间成本。纯离散模型步长可以放宽到2e-6，但要注意延时补偿是否考虑进去。

5.2 基于扰动注入的阻抗扫描

阻抗扫描的核心思路，是在PCC处叠加一个小幅值的正弦扰动电压，测量对应的变流器电流响应，做频域傅里叶分析后计算Z(f)=V_perturb(f)/I_react(f)。这个过程要扫若干频率点，每个频率点做一次仿真，最后把所有点拼接成阻抗频响曲线。

实验操作中，扰动源可以是串联电压源，也可以是并联电流源，推荐使用串联电压源形式，因为它更符合“电网戴维南等效”的思路，物理含义清楚。扰动幅值取额定电压的1%到3%，频率点从1Hz扫到2kHz，对数均匀分布30到60个点即可。每个扫描点需要让仿真运行足够时间，比如1到2秒，让扰动稳态建起来。

阻抗提取时的关键技巧，是正确处理基波分量的干扰。FFT输出里，工频分量的幅值远大于扰动分量，如果直接拿FFT结果计算阻抗，50Hz附近的数据基本会被污染。我惯用的做法是先对时域信号做带通滤波，只保留当前注入频率附近窄带（±0.5Hz内）的能量，再用窄带FFT提取幅值和相位。MATLAB里可以用spectrogram或者自编的Goertzel算法实现。

5.3 理论曲线与仿真结果的交叉验证

拿到扫频阻抗数据之后，第一步不是直接画Nyquist图，而是先把理论和仿真阻抗曲线叠在一起对比。如果两者在次同步和超同步频段的幅值和相位都吻合，说明阻抗建模是正确的；如果偏差明显，优先怀疑PLL线性化模型里的相位偏差没写对，或者LCL参数在仿真中引入了额外的串联电感。

对比完阻抗曲线，再画Nyquist判据曲线。此时有两种数据来源：一是用解析传递函数直接画，二是用扫频数据插值后画。两个来源的Nyquist图如果趋势一致，说明稳定性结论可信。若两者稳定裕度差异大，通常是扫频点数不够或者FFT提取误差偏大，需要加密频点或者延长仿真时间。

我习惯再做一步时域验证。在阻抗分析显示稳定的参数下，给电流参考或者电网电压一个阶跃扰动，观察系统响应是否收敛；在显示不稳定的参数下，观察扰动之后波形是否出现次同步或超同步振荡。这一步闭环验证做完，整个“理论-仿真-判定”链条才算完整。

6. 仿真与工程应用中的常见问题与排查心得

6.1 常见问题速查表

下面这张表是我个人调试过程中总结的高频问题，每一条都是真实踩过的坑。我尽量把现象和检查方向写清楚，方便你对照排查。

表格里每一条背后都值得再展开聊，这里拣两个最典型的说透。第一行那个问题，是FFT边窗泄漏导致的。你注入1Hz扰动，FFT结果里却可能在1.2Hz出现能量，这是窗函数主瓣宽度有限造成的。解决办法是让仿真运行时间至少包含20个扰动周期，同时用Flattop或者Kaiser窗替代矩形窗，虽然幅值精度会打折扣，但频率定位准确得多。

第二行的耦合项问题，在dq阻抗法里特别隐蔽。有些人只推导了Zdd和Zqq，却把Zdq和Zqd置零简化处理。强电网条件下这种简化误差不大，但弱电网下耦合项正是次同步谐振的关键路径，丢掉它们之后得到的Nyquist曲线会偏乐观，可能导致实际系统已失稳而理论判定稳定。所以建模阶段一定把PLL耦合项完整写出来。

6.2 参数整定的三条建议

第一，PLL带宽和电流环带宽要做“频段隔离”。我的经验是PLL带宽设为电流环带宽的十分之一以下，比如电流环带宽800Hz，PLL带宽最多取到几十Hz，弱电网下取10到20Hz比较稳妥。这样PLL只在低频段起作用，电流环负责中高频段，两个环的交互影响可以降到最低。

第二，LCL阻尼不要完全依赖无源阻尼电阻。虽然无源阻尼在工程上最简单可靠，但阻尼电阻会引入持续损耗。研究阶段可以先做有源阻尼（比如电容电流反馈阻尼）验证稳定性，工程落地时再根据成本和散热决定最终方案。关键是先把阻抗模型里的阻尼环节写对，别在建模时把阻尼忽略掉，否则Nyquist分析的结论会和实际差很远。

第三，参数敏感性分析一定要做。至少对Lg、PLL带宽、电流环Kp三个变量做一组扫描，把每组参数下的Nyquist曲线最小距离（稳定裕度）记录下来。这个最小距离可以用min(abs(L+1))计算，就是一个数值指标，适合做参数寻优的目标函数。这样一来，稳定性分析就从“看完曲线给结论”升级成了“量化指标+自动寻优”，做课题和工程报告都会顺畅很多。

我个人实际做这套分析的时候，最深的体会是：阻抗建模的推导过程虽然繁琐，但它给你的回报是清晰的可解释性。时域波形发散只能告诉你“有问题”，阻抗曲线却能告诉你“问题在哪个频段、跟哪个参数强相关”。稍微花点时间把模型和仿真打通，后续所有参数优化和故障分析都能站在同一个平台上展开，不需要每次从头试错。这套流程我反复用过多次，从并网逆变器到储能变流器都能迁移。希望这份整理能让你少走几条弯路，把有限的精力放在真正需要判断力的分析上。