做VSG虚拟同步发电机并网稳定性分析的时候别一头扎进状态空间法就出不来。这几年在新能源并网领域基于序阻抗的频域分析方法热度很高尤其当你的VSG控制器里带着电压电流双闭环电压外环、电流内环时状态空间模型阶数高得吓人参数建模稍微出点偏差特征值轨迹就会完全不同排查起来特别痛苦。而序阻抗扫频提供了一条更工程化的路径在时域仿真模型里注入扰动、扫频跑完得到机组在正序和负序下的阻抗频率特性再结合电网阻抗去判断系统稳定与否。这篇文章就是我基于 VSG 序阻抗扫频、时域阻抗稳定性分析以及建模仿真的一整套完整记录既适合刚上手VSG仿真搭建的人也适合想用频域阻抗法替代纯时域试凑的工程师。1. VSG稳定性分析为什么绕不开序阻抗1.1 状态空间法在带双闭环VSG面前的吃力表现说到并网逆变器的稳定性分析很多人第一反应是用状态空间法在dq坐标系下拉小信号模型然后求特征值。这个方法在单环路的PQ逆变器上确实好用系统阶数不高参数提取也成熟。可一旦换到VSG上情形就完全不一样了。VSG的控制链路本身就比普通逆变器长虚拟调速器、虚拟励磁调节、功率计算、电压外环、电流内环再加上锁相环和滤波电路。每一环都会引入状态变量整个系统模型动辄十几阶甚至二十几阶。模型阶数上去了参数精度就成了硬约束。任何一个PI增益差个百分之二十特征值轨迹都可能从稳定区跑到不稳定区但你很难判断到底是哪个环节的参数偏差带来了这个结果。另一个问题是弱电网下的动态交互。状态空间法里电网阻抗一般作为固定参数带入模型可实际运行中电网阻抗是随时变化的。想要扫描电网阻抗对稳定性的影响相当于要做大量参数扫描每一次都要重新线性化、重新算特征值工程效率很低。VSG的设计初衷恰恰是为了提升弱电网适应能力电网阻抗变化时的稳定性是最核心的关注点状态空间法在这个场景下就显得笨重。1.2 序阻抗的“端口思维”不看内部只看外部响应序阻抗方法的思路完全不同。它不关心变换器内部到底有几个环、系统有多少阶而是把整个VSG当作一个端口网络。从外部看进去在不同频率下施加电压扰动测量电流响应或者施加电流扰动、测量电压响应就能得到这个端口网络在该频率下的等效阻抗。把正序、负序分开测就得到了正序阻抗Z_p(s)和负序阻抗Z_n(s)。这种做法有点像“黑盒测试”你不需要精确知道盒子里每个元件的参数只要在端口上施加一组已知扰动记录响应就能推断出盒子对外的频域行为。对大电网来说VSG可以等效成一个阻抗源电网也可以等效成一个阻抗两者串联在一起稳定性自然就变成了阻抗比或回差矩阵的奈奎斯特判据。对带电压电流双闭环的VSG来说序阻抗方法还有一个额外优势双闭环参数对阻抗特性的影响是连续的、可观测的。你可以在扫频结果上直接看到电压环带宽附近的阻抗峰、电流环带宽附近的相位变化这种直观性远好过对着几十个特征值猜。这也是我在实际项目中坚决转向扫频法的原因。1.3 VSG的阻抗特性不能直接拿普通逆变器结论套用很多人初期会犯一个错误把普通PQ逆变器的序阻抗结论直接挪到VSG上。实际上两者阻抗形态差异很大。普通逆变器的中频段阻抗主要由锁相环主导而VSG有功率外环引入的是同步转矩和阻尼转矩这相当于在阻抗模型中附加了虚拟惯量环节和阻尼环节。惯量时间常数主要影响阻抗低频段的幅相特性阻尼系数影响振荡峰附近的动态。电压电流双闭环又进一步改变中高频段阻抗形态。所以VSG的序阻抗不是简单的RL串联模型也不是单一PI控制器的输出阻抗而是功率环、电压环、电流环共同塑造的、带明显频段特征的曲线。在正序和负序上也有差异正序阻抗会明显受到功率外环和锁相环影响负序则主要经过电压电流双闭环和LC滤波器控制路径不同阻抗特性自然不同。负序阻抗里锁相环的影响很小主要体现双闭环和滤波器的频域特征。做稳定性分析时如果不区分正负序很容易得出错误结论这一点后面我会细说。2. 从控制结构到序阻抗双闭环参数怎么“写进”阻抗曲线2.1 典型的带双闭环VSG控制结构拆解VSG仿真搭建现在基本是标准套路直流电压源接三相逆变桥经LC滤波器并入电网。控制部分分成三层功率外环、电压外环、电流内环。功率外环采集有功无功经过虚拟调速器和虚拟励磁器输出参考电压的幅值和相位电压外环把LC电容电压控制在给定值附近输出电感电流参考电流内环快速跟踪这个参考输出调制波。还有一个可选的是锁相环用于并网同步但在纯VSG模式下很多方案会弱化锁相环作用甚至直接采用功率同步。我在Simulink里搭建双闭环时全部采用PI控制器电压环输出限幅和电流环输出限幅都必须加否则启动瞬间积分饱和会带来很大的超调。电流内环一般是dq轴解耦加电网电压前馈电压外环也设计成dq双通道。电压外环输出的限幅直接影响电流内环参考值所以在阻抗扫频时扰动幅值不能太大否则控制进入非线性区测量出来的阻抗就不再是小信号意义下的结果了。下面是一组我在50kW VSG仿真模型里实际使用过的控制参数可以作为搭建参考控制环节参数项数值说明有功频率环虚拟惯量H0.5 s模拟转子惯量有功频率环阻尼系数D20模拟阻尼绕组效应无功电压环无功积分系数200电压参考调节速度虚拟励磁电压调节系数0.02无功-电压下垂系数电压外环Kp / Ki0.5 / 100带宽约1kHz电流内环Kp / Ki5 / 500带宽约8kHz2.2 双闭环带宽设置如何改变阻抗曲线形状阻抗扫频结果和双闭环参数之间的对应关系可以总结成几个经验规律。电压环带宽决定阻抗曲线在低频到中频段的抬升程度。电压环PI比例系数增大电压环带宽提高VSG在中频段的正序阻抗幅值会明显抬高相位滞后减小与电网阻抗的交点频率会向高频移动。电流内环带宽通常设置为电压环的5到10倍它决定阻抗曲线在更高频段的形状。电流环带宽越高高频段阻抗相位越接近90度呈现感性有利于抑制高频谐振。反过来如果电流环带宽不足阻抗在高频段会出现明显相位跌落甚至呈现容性特征。在高比例新能源系统里容性阻抗和电网感性阻抗一旦发生幅值交叠很容易形成串联谐振条件。我在扫频对比中把电流环穿越频率从2kHz降到800Hz阻抗曲线在高频段的相位从接近90度降到接近-10度原本稳定的系统直接变成了不稳定。这个试验让我深刻意识到双闭环带宽不只是时域的跟踪性能指标它直接决定了频域下的阻抗形态。控制参数方向低频段阻抗中频段阻抗高频段阻抗虚拟惯量增大幅值升高惯性增强振荡峰位置向低频移动几乎无影响阻尼系数增大低频增益略降峰值回落振荡抑制明显几乎无影响电压环比例增大低频抬高明显中频相位滞后减小轻微影响电流环比例增大几乎无影响过渡段变平缓高频相位接近90度2.3 正序负序分解与交叉耦合问题扫频时不能只测单相阻抗正序和负序要分开处理。实际三相系统中注入正序小扰动会同时引起正序和负序电流响应只是因为三相不对称或控制不对称正负序之间存在耦合。理想情况下dq对称控制中负序阻抗只和负序电流响应有关但VSG功率外环会在正负序之间引入耦合路径特别是功率计算环节会把正序电压与负序电流的乘积项带进来。因此在数据处理时正序阻抗Z_p V_p / I_p负序阻抗Z_n V_n / I_n同时还要计算交叉耦合项Z_pn V_p / I_n和Z_np V_n / I_p。做稳定性判断时如果系统三相基本对称、控制对称耦合项很小可以用简化单序阻抗判据如果电网存在严重不平衡耦合项就不能忽略需要采用2x2回差矩阵行列式判据。实际项目里我的处理方法是先扫出正负序各自的自阻抗再评估耦合项幅值。如果耦合项比自阻抗小20dB以上可以忽略如果只差10dB左右就必须考虑扩展判据。这一点在弱电网不对称故障场景下特别重要只看正序阻抗一定会漏掉负序通道的隐患。3. 时域阻抗扫频的完整操作链3.1 仿真模型搭建与测量参考面选取仿真平台我建议直接用Matlab/Simulink加Simscape Electrical。主电路用离散求解器步长取1e-6秒控制部分用离散PI采样率要和PWM载波频率匹配。扫频时电网模型用理想电压源串联RL阻抗在并网点也就是PCC和逆变器输出之间串联扰动注入模块。测量点选取是新手最容易忽视的点。计算并网阻抗时电压要取PCC处的相电压电流取逆变器出口电流两者最好在同一时刻采样避免因数据传输延迟引入相位误差。如果模型里加了滤波器阻抗参考面不同结果差异很大。输出滤波电容两端的电压和电网侧电压不是一回事。我做扫频时统一把测量基准放在LC滤波器电容电压和滤波电感电流上这样控制闭环的作用范围才和测量范围一致。参数设置如下模型模块参数数值备注直流源电压800V稳压源逆变桥开关频率10kHz离散PWMLC滤波器电感0.5mH滤波电感LC滤波器电容50uF滤波电容电网阻抗电阻/电感0.1Ω/1.2mH弱电网工况控制采样周期离散步长1e-6s求解步长3.2 扰动注入方式单频逐点扫频对比Chirp与PRBS时域阻抗测量的核心是“注入扰动、采集响应、提取频域数据”。扰动信号的选择直接决定扫频精度和耗时。单频逐点扫频是最稳妥的办法。在PCC电压参考或电流参考上叠加一个指定频率的正弦小扰动等系统达到稳态后采集一到两个周期的数据提取该频率下的电压电流幅值和相位算出阻抗。优点是基本没有频谱泄漏问题精度最高缺点是每个频点都要等待系统进入稳态50Hz到3000Hz范围内取一百个频点仿真时间就可能很长。Chirp信号扫频是把一个线性调频信号一次性注入一次仿真就能覆盖全频段适合快速摸底。但Chirp注入的扰动能量分散在各频率上幅值大小、扫频速率、傅里叶变换边界的处理都会影响精度不建议在需要高精度的参数整定验证中使用。PRBS伪随机序列激励在电力电子小信号辨识里也比较常见频谱近似平坦且能量集中但PRBS的幅值选择要非常小心幅值过大会把非线性环节激出来幅值过小信噪比又不够。我在双闭环VSG扫频中首选单频逐点法只有在需要快速观察全频段趋势时才用Chirp粗扫然后对敏感频段再做单频点精扫。扰动方式精度速度适用场景主要注意点单频正弦逐点高慢最终结果验证每个频点需等待稳态Chirp线性扫频中快趋势粗扫需抑制频谱泄漏PRBS伪随机中高快在线辨识幅值要反复试凑3.3 FFT窗口、频谱泄漏与Goertzel算法细节数据后处理是基本功。拿到时域波形后要取稳态段的多个周波做FFT。首要原则是每个频点的采样数据长度尽量包含整数个扰动周期否则就会产生频谱泄漏。实际做法是生成扰动频率表时先确定采样窗口时长是扰动周期的整数倍同时满足和基波周期的最小公倍数关系这样FFT后就能直接读出准确的幅值和相位。更精细的做法是用Goertzel算法在指定频率上做离散傅里叶变换。Goertzel比整段FFT更适合逐点阻抗计算因为它只计算关心的频点计算量小且精度高。具体做法是对电压和电流信号分别乘以cos和sin参考信号在窗口长度内积分得到复振幅。公式可以写成V_k (2/N) * Σ v[n] * e^(-j2πk*n/N)然后在同一频点下计算Z(k) V_k / I_k。为了抑制基波和谐波干扰加窗或者差分都可以但加窗后幅值要乘以窗函数带宽的校正系数。我在工程上为了避免麻烦干脆用整数周期截断加Goertzel既省事又准确。特别提醒扫频结果要做三相合成。把abc三相电压电流通过对称分量法变换为正序、负序和零序分量。由于VSG通常没有中线零序分量很小可以忽略只关注正序和负序。合成时务必使用同一个扰动频率的数据不能把不同频段FFT的结果混在一起用否则正序负序相位基准会乱掉。4. 从扫频数据到稳定性结论4.1 扫频数据的平滑与模型验证单频逐点扫完几十个频点后需要对阻抗曲线做平滑处理因为时域仿真总有少量数值噪声。平滑可以用局部多项式拟合也可以用Matlab的smoothdata函数。人工检查时重点看几个特征频率附近工频附近、功率外环带宽附近、电压环带宽附近、LC谐振峰附近。这些区域的阻抗变化率大往往是稳定性判据最敏感的频段。我习惯把扫频数据和解析小信号模型画在同一张图上对比。如果模型参数准确两条曲线在主要特征频段应该吻合。如果不吻合优先检查测量参考面是否正确、锁相环动态是否包含在内、电网阻抗是否被错误叠加。这种相互验证是阻抗法建模里最有效的纠错手段比单纯相信任何一方的结果都可靠。4.2 阻抗比判据与奈奎斯特判据的工程化操作拿到VSG的序阻抗Z_vsg和电网阻抗Z_g之后稳定性判断通常看Z_g / Z_vsg的奈奎斯特曲线或波特图。原理上两个阻抗串联回路的回差函数是1 Z_g / Z_vsg如果这个回差函数在右半平面没有零点系统就是稳定的。工程上用奈奎斯特曲线去判断就是看Z_g / Z_vsg的奈奎斯特曲线是否包围-1点。实际操作中我更常用波特图判据当Z_vsg和Z_g的幅值曲线相交时看该频率处的相位差。如果相位差小于180度说明交点处相位裕度为正系统稳定如果大于180度就存在振荡风险。三相系统要分别检查正序和负序因为负序阻抗对应振荡模式与正序不同只看正序一定会漏掉负序通道的问题。幅值交点处的相位差稳定性结论裕度情况小于180度稳定有正相位裕度等于180度临界状态相位裕度为零波形会有持续振荡大于180度不稳定负相位裕度发散振荡无幅值交点通常稳定无需额外分析4.3 一个具体算例改变阻尼系数后的稳定性变化我举一个自己实际跑过的例子。一台50kW VSGLC滤波器电感0.5mH、电容50uF电流环带宽约8kHz电压环带宽约1kHz虚拟惯量H取0.5秒阻尼系数D取20。扫频范围10Hz到1000Hz电网短路比在5左右。这个工况下正序阻抗和电网阻抗大约在320Hz处发生幅值交点。该频率处VSG阻抗相位约-85度电网阻抗相位约90度相位差175度稳定裕度只有5度系统勉强稳定实测时域波形上已经看到衰减非常慢的振荡。把阻尼系数D从20提高到50后同样频率处VSG阻抗相位变成了-60度相位差变成150度裕度增加到30度时域波形上振荡明显收敛。这个例子说明阻尼系数D对中频段阻抗相位有明显提升作用稳定性设计与VSG阻尼参数直接相关。如果不做阻抗扫频仅靠在时域波形上观察很难直观看到这种参数变化带来的稳定裕度差异因为时域波形上的振荡衰减快慢本身就不容易量化。5. 扫频实测中绕不开的坑与工程建议5.1 扰动幅值选择的“两难”困境扰动幅值太小FFT结果淹没在开关纹波和数值噪声里阻抗曲线毛刺很多幅值太大VSG的限幅环节和非线性环节被触发测出来的阻抗已经偏离小信号意义。我的做法是先在一个频点上做幅值扫描从额定电压的0.5%逐步加到3%观察该频点阻抗读数是否稳定然后选取稳定区间内的最大幅值作为整个扫频的注入幅值。通常在弱电网下取1%到2%比较合适。这个“两难”在带双闭环的VSG上更明显因为电压外环和电流内环的限幅幅值不同相同的电压扰动下电流内环可能先进入饱和。所以每次调整控制参数后我建议把注入幅值验证重做一遍不能一套幅值走到底。5.2 锁相环、直流母线动态对扫频结果的污染带锁相环的VSG在扫低频段时PLL动态会与功率外环产生耦合导致正序低频阻抗出现明显的负阻特性这是真实存在的物理现象不能当成数值误差去掉。如果仿真模型没有锁相环而是用功率同步方式则要注意功率同步环的低频特性可能导致阻抗在某个频段出现负电阻。这个负阻现象恰好可以解释VSG低频振荡的来源。直流母线也不能做理想化处理。母线电容偏小或直流电压控制环带宽和扫频段重叠时会在阻抗曲线上额外汇出一个尖峰干扰稳定性判断。建议扫频前把直流电压波动控制在1%以内或者直接加足够大的母线电容把直流动态从扫频频段里屏蔽掉。我自己最开始就是因为在理想直流源下扫频没有考虑母线动态后来换成带有限功率直流源就多了个尖峰折腾了半天才发现是直流动态混进来了。5.3 批量扫参与并行计算的经验仿真扫频最烦的是逐频点逐参数重复跑。我的做法是把扫频频点和参数组合写成一个批处理脚本在Matlab里用parfor并行跑。每个频点用不同随机扰动初始相位取两次平均既能提高信噪比又能把多核算力用完。参数扫描时先跑低频段粗扫锁定振荡频段后再细化。这套流程下来一条完整阻抗曲线大概50个频点在主流台式机上大约需要10到20分钟工程上完全能接受。还有一个小技巧扫频时把每个频点的电压电流时域波形也存一份方便后面排查异常数据。有一次阻抗曲线在某个频点突然跳变我翻出对应的时域波形发现是那个频点的扰动和PWM载波谐波发生了拍频调整扰动注入时刻后问题就消失了。这种小问题只看最终阻抗数据是发现不了的。最后再分享一点个人体会。阻抗扫频看似只是把时域仿真数据拿去做FFT但真正用起来之后会发现它把小信号稳定性分析从“推导解析模型、计算特征值”的象牙塔里拉了出来变成一套可以反复验证、误差可控的工程工具。尤其对带电压电流双闭环的VSG这种高阶控制系统扫频法几乎是目前性价比最高的稳定分析手段。你不需要把每个控制环的动态方程都推得清清楚楚只要保证扰动设计合理、数据后处理规范就能得到可信的序阻抗特性进而判断系统稳定裕度。下一步我打算把这套流程做成半自动化的扫频分析脚本顺手把正负序耦合项的完整判据也实现进去等后面有了更多实测数据再写一篇更细的案例对比和大家分享。
