在做并网变流器稳定性评估时，LCL-VSC阻抗建模这一个动作，基本决定了后面所有判断的成败。弱电网下的次同步和超同步谐振，归根结底是变流器输出阻抗与电网阻抗在某一频段出现了负阻尼交互，而Nyquist判据则是这套体系的最终裁判。这篇文章把我用Simulink从零搭起整套阻抗建模与稳定性验证流程的完整过程整理出来，重点拆解dq域小信号推导、PLL与电流环对阻抗频段的分区影响、电网阻抗扫描方法以及Nyquist曲线判读的实操要点。内容不涉及太复杂的数学推导，但会讲清楚每一步背后的物理意义，适合正在做新能源并网稳定性分析、变流器控制参数调优，或者准备相关课题仿真的工程师和研究生参考。
1. 先说清楚：弱电网下LCL-VSC为什么容易“抖”
1.1 VSC并网系统的“地基”：LCL滤波与它的双面性
三相电压源型变流器（VSC）通过LCL滤波器接入电网时，电感与电容构成的二阶网络会引入一个谐振峰。这个谐振峰在理想强电网下往往不是问题，因为电网阻抗几乎为零，阻尼相对充足。但在弱电网条件下，电网阻抗从零变成不可忽略的感抗，原本被抑制的谐振问题会被重新激发，甚至衍生出新的振荡频段。
我看到不少刚接触这个课题的同学把LCL谐振理解成“一个固定的频率点”，这是最大的误区。LCL谐振峰不是固定的，它与电网阻抗、控制带宽、运行功率都有关。比如典型的参数组合：L1=2mH，L2=1mH，Cf=20μF，无电网阻抗时的等效谐振频率大约为1.38kHz。这个频率听起来离工频50Hz很远，似乎不会参与次同步震荡。但当电网侧感抗Lg增大到10mH时，等效串联谐振频率会从约1.38kHz一路降到866Hz附近，而电容与网侧电感形成的并联支路谐振点也会向低频偏移。更重要的是，控制环路（尤其PLL和电流环）在低频段会引入负阻特性，一旦负阻区间与系统谐振点重叠，振荡就出现了。
所以研究LCL-VSC的稳定性，第一步必须把LCL滤波器和电网阻抗当做一个整体来看。LCL不是孤立元件，它是整个电流闭环的一部分。忽略电网阻抗谈LCL传函，等于只看地基不看地基下面是什么土。
1.2 弱电网的真实面貌：不是电压低，而是阻抗大
很多人第一反应“弱电网就是电压偏低”，实际工程中的弱电网定义更倾向于基于短路比（SCR）来描述。短路比是指电网在该接入点的短路容量与变流器额定容量之比，SCR越低说明电网越“弱”，变流器接入后对电网电压的影响越大。常见分级是SCR大于10算强电网，3到10之间属于中等强度，低于3就进入了弱电网区间，某些极端场景甚至低于1.5。
弱电网在电气特性上最直接的表现，就是并网点的戴维南等效阻抗显著增大。绝大多数分析场合可以把这个阻抗近似为纯感抗Lg，因为电网线路的高压侧R/X比值通常很小，电阻分量带来的阻尼可以暂时忽略。但要提醒一句，纯感抗假设会让稳定性分析偏保守，尤其是Nyquist曲线靠近-1点时，电网电阻的微小阻尼可能决定系统是否真的失稳。我做仿真时习惯在弱电网模拟里加一个很小的Rg，比如0.1倍感抗值，用来对标实际工况。
弱电网带来的核心问题是系统机电阻尼被削弱。把电网阻抗画进回路，相当于在L2之后又串联了一个Lg。这个新增的电感既会拉低LCL的固有谐振频率，也会改变电流环在并网点的开环增益相位。如果你从控制角度看，电网感抗本身不产生负阻尼，但它会在某一频段把变流器输出阻抗的相位“拉”到危险区间，让本来稳定的系统失去裕度。
1.3 次+超同步振荡：从阻抗视角看谐振问题
标题里提到的“次同步+超同步谐振”，在实际波形上体现为电流电压出现低于50Hz或高于50Hz（但通常低于几百Hz）的持续振荡分量。比如次同步振荡典型落在5Hz到45Hz，超同步振荡落在55Hz到200Hz左右。这类现象在风电、光伏经柔直或长线路并网场景里已经报道了很多次。
从阻抗视角去理解这个问题非常直观。变流器在dq旋转坐标系下可以等效成一个诺顿模型，即一个受控电流源并联一个输出阻抗Zvsc。电网侧则是电压源串联阻抗Zg。系统稳定与否，取决于Zg和Zvsc在频域上的交互关系。如果在某个频率下Zg的幅值等于Zvsc的幅值，且两者相位差达到180°，那在这个频率点就形成了负阻振荡条件，对应时域上幅值不断增长或持续振荡的谐振信号。
为什么低频段特别危险？因为PLL和电流环在低频段的动态响应直接影响变流器等效阻抗的相位。PLL带宽通常设计在10到50Hz，正好落在次同步频段，PLL在追踪电压相位时引入的动态扰动会让VSC在这个频段呈现出负电阻特性。电流环带宽一般在300到1000Hz，其闭环动态会影响到超同步频段的阻抗幅相特性。这就是“次同步找PLL、超同步找电流环”这一经验说法的来源，后面第5节我会用扫参结果验证这个规律。
2. 阻抗建模：一步步把VSC“翻译”成阻抗表达式
2.1 为什么选择阻抗法而不是特征值法
做稳定性分析有几条路：状态空间特征值法、时域仿真法、阻抗频率法。状态空间法精度高，但需要把变流器主电路、控制环路、PLL全部写成大矩阵状态方程，系统一复杂，矩阵规模迅速膨胀，参数物理意义也不直观，尤其当你想分析“电网阻抗从5mH逐步增大到15mH时系统什么时候失稳”这类参数扫描问题时，每改动一个参数都要重新推导和验证状态矩阵，工作量很大。
阻抗法的优势在于“解耦”和“模块化”。变流器归变流器，电网归电网，两边各自建模（或者用仿真扫频测出来），最后通过阻抗比或者回比矩阵来判断稳定性。电网部分就是一个电感Zg=jωLg，极其简单。变流器部分可以先在理想强电网下标定好输出阻抗模型，再叠加电网阻抗的影响。这样参数扫描时只需要重画曲线，不用重新推导整个状态方程。
阻抗法的另一个好处是物理直觉清晰。看到阻抗Bode图上某段相位跌到-90°以下（等效负电阻），你能马上判断哪个频段有阻尼风险，然后针对性调整PLL带宽或电流环参数，这是状态空间法很难做到的直观性。
2.2 dq旋转坐标系下的VSC闭环阻抗推导
LCL-VSC的阻抗建模主流做法是在dq旋转坐标系下进行小信号线性化。为什么用dq而不是直接用abc？因为在三相平衡系统中，基波正序分量在dq坐标系下变成直流量，可以在静态工作点附近做线性化，得到线性时不变（LTI）模型，进而直接套用传递函数和Nyquist判据。如果在abc坐标下，50Hz基波本身是交变量，小信号模型会是线性时变的，处理起来极其麻烦。
具体推导时，先把VSC控制部分拆成几个环节：电流环PI调节器、PWM惯性环节（通常等效为一阶延迟，时间常数约1.5倍开关周期）、滤波电路（LCL）、PLL动态。电流环在dq下的典型结构是：测量电感电流或并网电流→与指令值比较→经过PI控制器→叠加电压前馈和解耦项→生成调制波→PWM输出桥臂电压。
小信号扰动分析时，在每个环节求偏导数。电流环产生的阻抗效应最直观：电流闭环使得VSC对外呈现一个受控电流源特性，其输出阻抗等于电流环开环增益除以闭环增益后的等效表达式，公式形式大致为Zout = (sL_total + G_i(s)·H_PWM) / (1 + G_loop(s))，频率低于电流环带宽时，闭环增益很大，输出阻抗很小，表现为理想电流源；频率远高于带宽时，输出阻抗变大，趋近于无源滤波电感特性。这就是为什么电流环带宽决定了VSC在高频段的表现。
PLL的引入会更复杂，因为PLL输出的相位扰动会同时影响park变换矩阵和逆变器输出电压的角度参考。当PLL带宽附近出现扰动时，VSC在正序阻抗中会表现出一个明显的负电阻凹陷，这个凹陷的深度与PLL带宽、电流环增益、有功功率大小都有关系。实际建模时，很多人把PLL和电流环解耦处理，先推导电流环阻抗，再把PLL的影响以附加阻抗项叠加进去，这种简化在PLL带宽远低于电流环带宽时精度足够。
2.3 PLL和电流环对阻抗频段的分区影响
把一张完整的VSC输出阻抗Bode图摊开来看，能很明显地看到两个“特殊频段”。低频段（跟PLL带宽同量级）受PLL影响最深，阻抗相位会从接近-90°或0°的地方往负方向偏移，甚至跌到-180°以下，这个负相位区间就是次同步振荡的温床。中频段（电流环带宽附近）受电流环PI参数和滤波电感影响最大，相位曲线可能出现尖峰或急剧转折，对应超同步谐振风险。高频段则主要是LCL电容和电感网络的固有特性，几乎不受控制器参数影响，在这个频段VSC更像一个无源LCL网络。
做阻抗建模时有三个细节值得注意。第一，解耦项对阻抗有影响，dq电流控制中通常加入ωL解耦来消除交叉耦合，这个解耦项在小信号模型中会改变阻抗矩阵的交叉项，对高频段影响较大。第二，数字控制延时不可忽略，PWM更新和采样保持带来的约1.5拍延时，会让阻抗相位在高频段显著滞后，LCL谐振峰附近的阻尼特性因此变差，这也是为什么纯连续域建模经常高估稳定裕度。第三，负序阻抗和正序阻抗并不相同，在dq对称控制下系统确实对称，但PLL动态本身不具有负序响应，因此负序阻抗主要表现为无源特性，而正序阻抗包含了PLL带来的负阻效应。
这部分推导做完以后，你手里会有一个(2×2)的dq阻抗矩阵，其中对角项表示dd和qq自阻抗，反对角项表示d轴与q轴之间的耦合阻抗。由于系统对称，多数情况下可以解耦成正序阻抗和负序阻抗两个标量，直接用Nyquist判据判断稳定性。后面第3节的判据判读就是建立在这个解耦基础之上的。
3. Nyquist判据验证稳定性：判据怎么用才不出错
3.1 从“两个阻抗”到“一个判据”的推导逻辑
弱电网下VSC并网系统可以简化为两个阻抗的串联回路。变流器侧用一个诺顿等效（理想电流源并联Zvsc），电网侧用一个戴维南等效（理想电压源串联Zg）。并网点的电压和电流由这两部分共同决定，而系统的闭环特征方程可以整理成1 + Zg(s)/Zvsc(s) = 0的形式。
这里Zg(s)/Zvsc(s)就是回比函数L(s)，也是Nyquist判据要分析的对象。如果L(s)的Nyquist轨迹穿越或包围了复平面上的(-1, j0)点，则说明特征方程存在右半平面根，系统在相应频段失稳。幅值条件|Zg|=|Zvsc|对应回比幅值为1（0dB），这正是两个阻抗幅值交叉的频率点；相位条件则是Zg与Zvsc的相位差达到180°，即回比相位为-180°。
用Nyquist判据之前有一个前提必须确认：Zvsc本身不能含有右半平面极点。如果变流器内部电流环或PLL设计不当导致输出阻抗存在右半平面极点，那么直接套用“不包围-1点就稳定”会得到错误结论。我在实际项目里先画出Zvsc的Bode图，观察有无明显的相位异常翻转，再用状态空间特征值法交叉验证一次，基本能排除这个风险。对于论文或课题研究，建议在代码里直接检查阻抗传递函数的极点实部。
3.2 Nyquist曲线在次/超同步频段的判读方法
拿到Nyquist曲线图以后，很多人盯着整条曲线看半天不知道重点在哪。我的经验是只关注两个区域：低频端（对应次同步频段）和中频段（对应超同步频段）。曲线从右侧接近原点的部分通常平缓，一旦开始向(-1, j0)点靠拢，就要小心了。
判断次同步风险的步骤是这样：把频率范围限制在1Hz到100Hz，观察Zg与Zvsc的幅值交叉点落在哪个频率。如果交叉频率在PLL带宽附近，并且该处相位差接近-180°，Nyquist曲线就会从左侧绕过甚至穿过-1点。超同步风险则要看100Hz到500Hz的交叉点，这个频段的相位主要由电流环延时和LCL谐振峰共同决定，一旦数字控制延时把相位拖到-180°以下，即使幅值有余量，曲线也会甩向危险区域。
我在项目里习惯把Nyquist曲线按频率着色（低频段暖色、高频段冷色），这样“哪个频段在逼近-1点”一目了然。这个方法用MATLAB自带的scatter函数加颜色渐变就能实现，比单纯黑线直观得多，写论文配图也更有说服力。
3.3 稳定裕度的量化：交叉频率与相位裕度
Nyquist判据不只是给出一个“稳定/失稳”的二元结论，它还能告诉你离失稳边界还有多远。最常用的量化指标就是相位裕度。在回比的幅值交叉频率f_c处，计算相位与-180°之间的差值，这个差值就是该频段的相位裕度。一般认为大于30°比较安全，10°到30°需要谨慎评估，低于10°则属于高失稳风险区间。
同时关注幅值裕度也很有价值，尤其是在相位已经接近-180°但幅值还没到0dB的情况。幅值裕度指的是回比相位达到-180°的频率处，幅值距离0dB还有多少余量。如果幅值裕度也在3dB以内，那个工况基本离振荡不远了。我通常把相位裕度和幅值裕度两个指标联合使用，分别标注在Bode图上，配合Nyquist曲线一起说明，审稿人和评审专家看到这种逻辑都比较认可。
还有一个工程上的小习惯：当Lg变化时，把不同Lg对应的Nyquist图画在一张图上，观察曲线族逐渐逼近-1点的过程，可以非常直观地给出“临界短路比”。这个参数对实际项目很有价值，因为它直接告诉业主方“你这个场站的接入强度低于某个值就会振荡”。
4. Simulink仿真实现全过程
4.1 主电路搭建：直流源-VSC-LCL-弱电网串联结构
Simulink仿真模型的骨架分五个部分：直流侧、三相VSC桥、LCL滤波器、并网公共耦合点、弱电网模拟。直流侧直接用一个理想直流电压源，电压设定750V，后面接一个足够大的直流电容（我常用3000μF），用来稳定直流母线并滤除开关纹波。VSC桥用Simscape Electrical里的Universal Bridge模块，选IGBT/Diode类型，导通压降和开关损耗按器件手册填写，仿真步长控制在1μs左右，否则开关瞬态容易不收敛。
LCL滤波器参数我建议先用理论公式预估，再在仿真里微调。L1=2mH，L2=1mH，Cf=20μF这组参数是我常用的起点，对应约1.38kHz的谐振频率。没有阻尼措施时，这个谐振峰在强电网下会因为开关数字延时出现轻微振荡，弱电网下则可能被激发成明显振荡。为了先排除滤波器本身的干扰，我建议第一版仿真先加一个小阻尼电阻（比如串在Cf支路，取值3Ω），等稳定性分析框架跑通了再改成有源阻尼，这样能隔离变量。
弱电网模拟部分用三相电压源模块（Three-Phase Source）串联可调电感Lg实现。Lg从0.5mH到20mH之间扫描，就覆盖了SCR从约17到不足1的完整区间。为了模拟更真实的电网，可以再串联一个小电阻Rg，取值为0.1ωLg。注意电压源内部阻抗一定要设成RL串联形式，而不是直接调“内部阻抗为零”，否则并联的Lg实际上会被短路掉，这是一个很容易踩的坑。
4.2 控制环路实现：双环+PLL的离散化要点
控制部分按工程惯例用双闭环，电流内环+电压/功率外环。如果课题只关注阻抗稳定性，外环可以先省略，直接给电流指令Id_ref和Iq_ref，比如Id_ref对应有功功率50kW，Iq_ref设0或按无功需求设定。这个简化能让你专心观察电流环与电网阻抗的交互，避免外环动态干扰。
PLL用三相SRF-PLL，就是经典的park变换+PI调节+角度积分。PLL带宽我起始设35Hz，这个值在弱电网研究中偏中间档，既有PLL动态又不是过于激进。电流环PI参数按典型I型系统整定，电流环带宽设在500Hz左右，Kp和Ki对应关系为Kp=带宽×L_total，Ki=带宽×R_total。由于L_total在弱电网下会变大，这个整定公式算出的Kp实际是偏保守的，仿真里可以在此基础上乘一个1.2到1.5的增益系数。
离散化是仿真模型和实际工程一致性的关键。控制器采样频率与开关频率对齐，都设10kHz，PWM用单边调制或对称规则采样。PI控制器用离散传递函数实现，积分器采用后向欧拉法离散，解耦项要使用当前的电感值。如果你用了变参数Lg扫描，解耦项里的L值也要跟着变，否则解耦失效会在阻抗相位上留下明显误差。
4.3 弱电网模拟方法：从固定感抗到短路比自动扫描
固定Lg的模型只是最基本的弱电网模拟，真正高效的做法是用一个MATLAB脚本循环改Lg参数值，每改一次跑一个工况，记录该工况下的频域数据和Nyquist曲线，最后汇总成一个“稳定边界扫描图”。
具体操作上，把Lg定义成模型里的Simulink参数（比如用Data Dictionary或者模型工作区变量）。脚本里用sim函数循环调用模型，每次更新Lg、把时域波形保存到workspace。扫描Lg时注意步长策略：在临界稳定区域附近（目测Nyquist曲线开始贴近-1点的区间）加密扫描步长，比如从每2mH一步加密到每0.5mH一步，这样能准确捕捉临界失稳点。整个批量扫描我一般在30分钟内跑完20个工况，如果用并行工具箱还能再压缩到15分钟。
如果你需要扫“短路比”而不是Lg本身，可以把Lg换算成SCR再作为变量。380V、50kW系统下，基准阻抗约2.888Ω，对应纯感抗约9.19mH时SCR=1。所以5mH对应SCR约1.84，10mH对应SCR约0.92，这组换算关系我经常用。要注意SCR数值和Lg不是线性关系，报告里写“临界SCR”时必须附带说明Lg的对应值，避免被读者误解。
4.4 阻抗扫描与Nyquist绘制的实操流程
完成主电路和控制环路的建模仿真后，最关键的一步是获取VSC输出阻抗频率响应。这个环节理论上可以用解析传函直接画Bode，但为了验证建模是否准确，必须做扫频验证。我用的是Simulink Control Design提供的频率响应估计工具，在电流指令Id_ref处注入扫频扰动，扫描范围设0.5Hz到5000Hz，覆盖次同步到高频LCL谐振的完整区间。
扫频工具设置里有几个参数直接影响结果精度：扰动幅值设为额定电流的2%到5%，太大会激发非线性谐波，太小则信噪比不足；扫频方式选正弦扫频，每个频率点至少保留8个完整周期再记录数据，低频点（比如5Hz以下）需要更长收敛时间；线性化工作点由稳态仿真自动获取，注意在扫频前先让模型运行到100ms稳态，避免初始暂态污染估计结果。
扫频结束后会得到VSC输出阻抗的幅频和相频曲线，把理论解析建模的曲线叠上去做对比。如果二者偏差小于3dB和10°，说明理论模型可信，此时基于理论阻抗模型绘制的Nyquist曲线可以放心使用。如果偏差大，优先检查PLL带宽参数是否在理论模型中正确映射、数字控制延时环节是否考虑完整，而不是怀疑扫描工具本身。
最后一步是Nyquist曲线绘制。把每个频率点计算得到Zg/Zvsc的复数值存入向量，用Nyquist函数绘制，频率范围从0.1Hz到5000Hz。我会在同一张图上叠加不同Lg对应的曲线族，并在原图标出(-1, j0)的位置和距离最近的曲线点对应的频率，方便直接读出危险频段。这部分完成后，整个“阻抗建模+稳定性验证”的闭环就算跑通了。
5. 关键参数与谐振规律的扫参结论
5.1 电网感抗从2mH扫到20mH：谐振点怎么迁移
用第4节的仿真模型做批量扫描，第一个看的是电网感抗Lg对系统稳定性的影响。Lg从2mH开始增加时，系统还很稳定，Nyquist曲线远离-1点，时域波形在扰动后3到5个工频周期内衰减完毕。Lg增加到5mH时，Nyquist曲线的低频端开始向-1点靠近，相位裕度从约45°降到25°，仍在安全区间。
Lg继续增加到8到10mH，情况明显变化。此时在80Hz附近（PLL带宽边缘）Nyquist曲线与实轴的交点距离(-1, j0)已经非常近，相位裕度掉到10°以下，时域波形开始出现明显的低频包络波动，频率在几赫兹到十几赫兹之间。这个振荡在频谱上对应次同步分量。继续增大到15mH以上，Nyquist曲线直接包围了-1点，时域波形发散，保护逻辑触发停机。这组扫参结果说明，纯感抗Lg增大是诱发次同步谐振的最直接因素，临界点大约落在SCR等于2附近（对应Lg约9mH），与工程里“SCR低于3要重点做稳定性复核”的经验高度吻合。
值得注意的一个细节是，超同步频段并未随Lg增大出现同样的恶化趋势，说明电网感抗主要作用于次同步区间。如果项目测得的振荡频率在100Hz以上，优先怀疑方向就不是电网感抗，而应该是电流环参数或LCL谐振峰。这个“频段-参数”对应关系能在故障排查时帮你快速锁定问题源头。
5.2 电流环带宽对超同步谐振的影响
固定其他参数，把电流环比例系数Kp从基础值往上调，观察阻抗Bode图的中频段变化。Kp增大意味着电流环带宽从500Hz提升到800Hz甚至1kHz，VSC输出阻抗的相位在中频段会有明显变化，具体表现为阻抗相位低谷向高频方向移动且凹陷加深，对应Nyquist图在超同步频段（大约150到300Hz）的曲线鼓包变大。
这正是超同步谐振的控制面。把Kp调到初始值的1.5倍以上，配合弱电网下的LCL参数，Nyquist曲线在约200Hz处与-1点的距离急剧缩小。时域验证时，在电流指令加入一个200Hz左右的小扰动分量，振荡不会衰减而是维持甚至放大。这说明电流环带宽并非越高越好，在弱电网下过度追求动态响应会牺牲超同步频段的阻尼。
想改善超同步阻尼，常规做法是加入有源阻尼算法。最简单的实现是电容电流反馈有源阻尼，把Cf电容电流引入电流环反馈路径，等效为在电容支路并联一个虚拟电阻。在Simulink里用一个高通滤波器提取电容电流的高频分量，乘以增益后叠加到调制波。我的调试经验是这个增益不能过大，否则高频噪声被放大，建议从0.01起步逐级增加，观察Nyquist曲线在超同步段的收缩效果。
5.3 PLL带宽对次同步谐振的直接影响
PLL带宽的扫参结果非常“干净”：PLL带宽从15Hz逐步提高到80Hz，次同步频段的阻抗相位凹陷越来越深，Nyquist曲线在10到50Hz区间的包围趋势显著增强。PLL带宽取80Hz时，即使Lg只有5mH（原本是稳定工况），系统也会在约35Hz处出现负阻尼振荡。这组对照实验把“PLL带宽越大越容易次同步振荡”这个结论变得非常直观。
PLL带宽为什么影响如此显著？因为PLL在追踪电压相位时，本质上是一个非线性反馈路径。当电网阻抗增大时，并网点电压受注入电流波动影响，PLL感知到的相角扰动被放大。PLL带宽越高，这种扰动被更“积极”地引入参考相位，电流环输出也随之波动，最终在VSC输出阻抗上表现为低频负电阻。想让系统在弱电网下稳定，PLL带宽应保持在30Hz以下，必要时可以加入基于电网阻抗的PLL自适应整定方案。
实际项目中还有一个折中方案：把PLL带宽降低的同时，加一个相位前馈补偿。在PLL输出角速度上叠加一个正比于无功功率波动的高通补偿量，可以在不明显牺牲动态响应的情况下削掉一部分负阻凹陷。不过这个补偿的系数整定依赖具体的运行点，仿真验证时要多跑几个功率点，确认对超同步段没有副作用。
5.4 参数表速查：稳定/失稳边界
参数扫完后，我习惯把结果整理成一张速查表，方便后续复现和汇报使用。以下数据基于50kW样机仿真，具体数值会随参数不同而漂移，但规律具有参考价值。
| 工况描述 | Lg (mH) | 电流环带宽 (Hz) | PLL带宽 (Hz) | Nyquist裕度评估 | 时域结果 |
|---|---|---|---|---|---|
| 强电网基准 | 1 | 500 | 35 | 裕度充足，相位裕度>45° | 稳定，扰动快速收敛 |
| 中弱电网 | 5 | 500 | 35 | 低频段裕度下降，约25° | 稳定，余量一般 |
| 弱电网临界 | 9 | 500 | 35 | 约80Hz处裕度不足10° | 出现次同步包络波动 |
| 极端弱电网 | 15 | 500 | 35 | Nyquist包围-1点 | 发散，保护动作 |
| 弱电网+高PLL带宽 | 5 | 500 | 80 | 约35Hz处裕度为负 | 次同步持续振荡 |
| 弱电网+高电流环增益 | 9 | 900 | 35 | 约200Hz处裕度骤降 | 超同步振荡风险 |
这张表最核心的结论可以浓缩成两句话：电网感抗增大是次同步振荡的总开关，而PLL带宽决定次同步频段负阻尼的深度；电流环带宽在中高频段影响超同步阻尼，弱电网下不能盲目追求高带宽。
6. 常见问题与排查记录
6.1 仿真不收敛或高频振荡发散
最常遇到的是模型刚启动就报“Simulink cannot solve”或者波形直接飞掉。第一个排查点是步长。整流逆变模型里开关频率10kHz，仿真步长超过5μs基本必发散，我一般固定步长1μs，每个工频周期20ms就是两万步，50kW模型跑一个0.5s工况大约3到5分钟，可以接受。
第二个排查点是LCL阻尼。如果模型不带阻尼电阻也不含有源阻尼，LCL谐振峰在数字控制延时下必然振荡。先加一个串联在电容支路的3Ω小电阻，确认基础稳定后再换有源阻尼，是最高效的排查顺序。
还有一个容易被忽略的问题：电网模拟电压源内部阻抗为0时，串联的Lg被旁路。排查方法是给电压源设置一个极小内部电阻（如1mΩ），或者把Lg模型改成“电网电压源+串联RL支路”的独立结构。我用的是后者，物理上更严谨。
6.2 阻抗扫频结果噪声大或者曲线毛刺多
扫频数据不平滑，首先检查扰动幅值。电流环指令中加入的正弦扰动太小，FFT提取的分量被开关纹波淹没；太大则进入调制非线性区，产生谐波干扰。经验值是额定电流幅值的3%左右，不要超过5%。
其次是收敛时间。低频扫频点（低于5Hz）需要更长观测窗口，如果每个频点只采样6个周期，FFT栅栏效应会带来明显误差。建议低频点至少采样12个周期，高频点可以适当放宽到8个周期，保证每个频点的估计置信度。
如果曲线在LCL谐振频率附近出现尖刺，还要检查一下扫频扰动是否与谐振峰发生交叠。我在扫频设置里把频点间隔加密到谐振频率附近的1/10倍频程，然后用平滑滤波（MATLAB的sgolayfilt）做后处理。平滑时要小心不要把真实的相位低谷抹掉，先看一次原始曲线再做平滑。
6.3 Nyquist判据结果与时域仿真对不上
这个坑很值得写出来。理论计算显示系统稳定，但时域波形却有衰减到一半又反弹的振荡，或者理论判据失稳而时域居然还能维持稳定。排查顺序建议这样：
先确认Zvsc有没有右半平面极点。如果控制参数设计不当，VSC输出阻抗本身就不稳定，这时直接用“Nyquist轨迹不包围-1点”判据会失效。我在第3.1节已经强调过，这里再提醒一次。
再检查阻抗模型是否覆盖了所有动态环节。很多简化建模只考虑了电流环+PLL，却忽略了电压外环、直流电压动态、PWM限幅饱和等环节。直流母线电容在小功率等级下影响不大，但大功率变换器直流侧纹波会耦合到输出阻抗里。仿真和理论对不上的时候，把理论模型换成扫描实测数据再算一次Nyquist，基本能定位偏差来源。
还有一个经验：时域扰动试验的初始扰动幅值如果太大，系统非线性阻尼会掩盖线性失稳特征，导致看起来“还能稳定”。做失稳验证时扰动幅值应该控制在额定电流的5%到10%，并且观察至少300ms以上，让振荡充分发展或衰减稳定后再下结论。
6.4 批量扫描耗时太长怎么办
弱电网扫参场景通常需要跑几十个工况，每个工况如果都要完整跑0.5s，累积时间确实可观。用并行工具箱把Lg和PLL带宽的组合分配给不同worker，能省一半时间。另一个思路是只仿真故障区间，比如先扫一位粗略的波形特征，找到临界点后再按小步长扫描，省掉大量稳定工况的冗余计算。
仿真加速还有一个取巧办法：把PWM载波频率在仿真中降到5kHz，同时把开关模型的导通电阻稍微调大，用阻尼损耗换取步长放宽到2μs。这样扫描精度略降但趋势判断完全够用。等锁定临界区域后再用精确模型确认最终临界点，这是工程仿真的常见策略。
7. 个人实操体会与一个小技巧
这套LCL-VSC阻抗建模与Nyquist判据验证流程，我已经在不同参数的样机上复现过多次，最大的体会是“建模的简化假设决定分析的边界”。你把数字控制延时当成理想环节忽略掉，超同步段相位就偏乐观；你把电网电阻完全忽略不计，次同步段判断又偏保守。做稳定性分析要想不被模型细节反噬，最好的习惯是理论、扫频、时域三线并进，任何两条对不上就先停一停，不要急着改参数。
最后分享一个能明显提升效率的小技巧：在Simulink里把PLL带宽、电流环Kp、Lg这三个参数用“模型参数对象”保存，扫参脚本只改对象属性，不直接改模型内部变量。这样做批量扫描时，模型各模块的标签值和实际参数始终一致，不会因为手动改数导致某个模块忘了同步更新。另外保存数据时把每个工况的Nyquist原始复数据一并存入mat文件，而不是只存图片，后续复现论文图或者追溯分析结论都需要这些原始数据。这个习惯帮我少踩了很多“跑完发现记录的曲线修复不回来”的坑。
