先说个事。我最初接触VSG虚拟同步发电机序阻抗建模不是因为看论文觉得有趣而是有一次在弱电网环境里做并网测试装置在某个低短路比工况下出现了持续振荡。示波器上那一串正弦波叠着次同步分量怎么看怎么像风电场“次同步振荡”的缩小版。当时第一反应是查锁相环参数但问题在于VSG根本没有传统意义上的PLL同步机制完全靠功率外环的转子运动方程。那时候手边只有dq阻抗的模型分析起来总觉得隔着一层。后来才下了决心把VSG的序阻抗模型从头到尾推一遍。这篇文章就记录我在这个过程中积累的东西从序阻抗模型到底能干什么、怎么从控制环路上拆VSG、谐波线性化推导的每一步怎么走到仿真扫频怎么验证、参数变化怎么移动阻抗轨迹再到我踩过的坑。内容偏硬核但我会把思路和细节都说清楚适合正在做并网逆变器小信号建模、或者被弱网振荡问题缠住没法下手的同行参考。1. 为什么要盯上VSG的序阻抗一次弱电网振荡的复盘1.1 现场现象与分析思路的转变那次振荡发生在短路比差不多3.0左右的电网条件下VSG并网发电输出大概60%额定功率。振荡频率约23Hz介于典型的次同步频段。现场一开始怀疑是电网侧的背景谐波激发的但测了PCC电压谐波5次、7次、11次含量都很正常。接着怀疑VSG内部的功率环参数但时域仿真里无论怎么调虚拟惯量振荡依然存在只是频率在移动。后来我把问题从“时域波形”拉回“频域阻抗”来想如果把逆变器看成一台等效电压源串联输出阻抗把电网看成理想电压源串联电网阻抗那么振荡的必要条件就是两者阻抗之和在某个频率处趋近于零也就是阻抗匹配。要让这个判断落地就必须有准确的逆变器输出阻抗模型而VSG的输出阻抗并不像传统跟网型逆变器那样只由电流环和LCL滤波器决定功率同步环、无功环、电流环、滤波环节全部耦合在一条链路上。当时手里最顺手的是旋转坐标系下的dq阻抗但拿到现场判稳和不对称分析时用起来并不直观于是才开始正经做序阻抗建模。1.2 序阻抗模型的本质从频域“黑箱”看逆变器序阻抗模型说白了是这样一件事在稳定的基波工作点上向逆变器端口注入一个小幅值的正序或负序电压扰动然后测量对应频率下的电流响应把电压扰动和电流响应的复数比定义为该频率下的正序或负序阻抗。如果扰动频率落在基波附近的边频带那么由于控制回路的非线性频率搬移电流响应会同时出现在扰动频率和镜像频率上这就产生了正序与负序之间的交叉耦合描述起来就需要一个2x2的序阻抗矩阵频率分量正序电压激励产生的电流负序电压激励产生的电流正序响应Z_ppZ_pn负序响应Z_npZ_nn大多数三相对称系统里正序和负序基本解耦Z_pn和Z_np量值很小但并非严格为零。VSG内部用旋转坐标系实现控制坐标变换和功率计算都会把正序扰动映射到直流分量再搬移回来只要控制回路存在不对称或谐波耦合耦合项就有存在的空间。对于做稳定性分析的人而言这个矩阵可以直接放进广义Nyquist判据也可以和电网阻抗一起做阻抗比分析不用再费劲把时域非线性模型线性化在某个稳态点上。1.3 与传统同步发电机建模的差别传统同步发电机的阻抗建模有成熟的理论比如Park方程和派克变换能推导出同步电机的d轴、q轴运算电抗然后通过坐标变换转成序分量。VSG虽然名字里有“同步发电机”但它的电磁关系、控制回路和物理电机完全不同。VSG的“转子运动方程”是写在控制器里的程序虚拟惯量和阻尼系数可以任意设定无功环、电流环、电压前馈、数字采样延迟都在同一条信号链路上。所以VSG的序阻抗无法直接从电机理论里抄必须老老实实从控制环路的传递函数出发把每个环节的频域效应叠加进去。这也是这个课题真正有意思的地方它把同步发电机的外特性与控制器的频域响应拼在一起形成一个既像电机又像电力电子装置的混合体。后续所有推导都围绕这条混合链路展开。2. 从控制回路入手拆解VSG的频域画像2.1 VSG的四个控制环路完整的VSG控制可以拆成四个层次功率同步环、无功-电压环、电流内环、滤波与调制环节。功率同步环负责把有功功率误差映射成虚拟功角无功-电压环负责把无功功率误差映射成电压幅值参考两者一起决定VSG端口电压的幅值和相位电流内环负责在电压指令和实际输出电压之间建立快速跟踪关系最后经过SPWM或SVPWM调制叠加数字控制延时后作用在逆变桥输出端再经过L或LCL滤波器连接到电网。对阻抗建模来说四个环节的作用频段是分开的。功率同步环的时间常数通常在几十到几百毫秒量级对应几Hz到几十Hz的低频段电流内环带宽通常设计在1kHz左右对应中高频段无功环也工作在低频但它的影响主要体现在电压幅值路径上与有功环的相角路径不完全一样。这种频段分离意味着全频段阻抗曲线实际上是由不同控制环分段“接管”的分析时不能一锅端。2.2 功率同步环的等效转子运动方程VSG的功率同步环在数学上就是一个带惯量的频率调节器[ J \frac{d\omega}{dt} P_m - P_e - D_p (\omega - \omega_g), \quad \theta \int \omega dt ]其中J是虚拟惯量D_p是阻尼系数P_m是机械功率参考P_e是电磁功率也就是VSG输出到电网的有功ω_g是电网额定角频率。这个方程在频域里引入了一个固有积分环节1/s所以有功扰动会在中低频段产生相位动态体现为阻抗轨迹中的一个谐振峰或负阻区。与同步发电机类似J越大这个动态的固有频率越低D_p越大阻尼越强。在小信号分析时把功率写成电压幅值和功角的函数然后在稳态工作点附近线性化就能得到扰动功角δ与有功功率扰动ΔP之间的二阶传递函数。这个传递函数并不会直接等于输出阻抗但它会通过电压指令路径改变逆变器端口的等效电压源扰动从而影响阻抗。2.3 电流内环和滤波器的作用边界电流内环的作用是把电压指令转换成电感电流指令的快速跟踪常见PI调节器加解耦项和电网电压前馈。在阻抗模型里电流内环的闭环带宽决定了一个重要边界低于电流环带宽时VSG端口的电压特性主要由功率外环决定表现为低阻抗电压源特性高于电流环带宽时电流环来不及跟踪输出阻抗逐渐被滤波电感主导表现为感性阻抗。这个边界在扫频曲线上通常表现为一个相位翻转点。滤波器结构也需要区分。L型滤波器的输出阻抗是一条干净的感性曲线LCL滤波器在谐振频率附近会出现反谐振峰如果不加有源阻尼阻抗曲线会在谐振点附近变得很怪。大部分VSG研究为了方便分析和控制设计用L滤波器但工程装置往往用LCL。建模时一定要明确滤波器结构否则扫频数据对不上理论。3. 谐波线性化下的序阻抗推导过程3.1 从三相对称扰动到正负序分量的映射推导核心方法是谐波线性化。思路是在三相静止坐标系下假设系统原来处于三相对称基波稳态然后向A相端口注入一个频率为f_p的小信号电压扰动。这个扰动的三相对称形式决定了它是正序还是负序正序扰动的三相相位依次滞后120°负序扰动则是依次超前120°。在频域上一个频率为f_p的正序电压扰动加到交流端会在静止坐标系中表现为两个频率分量一个是f_p本身另一个是f_p - 2f1f1是基波频率。后者是三相坐标变换时产生的镜像频率来源于旋转坐标系之间的频率搬移。举个具体例子如果基波50Hz注入80Hz正序扰动那么在A相里你会同时看到80Hz和-20Hz两个频率分量如果注入负序扰动则会出现80Hz和180Hz两个分量。这个频率配对规则是序阻抗建模最关键的一步做错这一步后面全乱。3.2 从时域扰动方程到频域阻抗表达式推导路径大致分三步。第一步把A相电压和电流写成基波稳态分量加小信号扰动分量的形式代入功率计算模块得到有功、无功的小信号表达式。第二步把有功、无功扰动通过功率同步环和无功环的传递函数映射为功角扰动和电压幅值扰动再合成出电压指令扰动。第三步把电压指令扰动送到电流内环和调制延时环节最终写成端口电压扰动与电流扰动的关系式。把上面的关系整理成矩阵形式就得到了VSG在频域的序阻抗模型典型写法是[ \begin{bmatrix} \Delta V_p(s) \ \Delta V_n(s) \end{bmatrix}\begin{bmatrix} Z_{pp}(s) Z_{pn}(s) \ Z_{np}(s) Z_{nn}(s) \end{bmatrix} \begin{bmatrix} \Delta I_p(s) \ \Delta I_n(s) \end{bmatrix} ]按理说推导到最后你会得到一个相对冗长的传递函数表达式里面包含VSC的基波工作点电压幅值、功率角、所有PI参数、滤波器电感、虚拟惯量、阻尼系数和数字延时。表达式长不等于复杂关键是每个环节在哪个频段起作用心里要有数。如果只是要一个能用的工程模型不需要把每个系数展开到极点用符号工具算完后留着化简即可。3.3 推导中做过的简化与边界条件建模过程中我做过几个重要简化这里说清楚免得你照抄出问题。第一控制器和电网都假设三相对称忽略锁相环VSG本来也没有PLL所以负序耦合项可以忽略Z_pn和Z_np近似为零阻抗矩阵退化成两个独立的标量阻抗Z_p和Z_n。这个假设在电网轻度不对称时仍然成立但严重不对称故障分析时必须回头考虑耦合项。第二调制环节用一阶惯性环节e^{-1.5Ts s}近似Ts是采样周期。这个延时影响高频段阻抗特别是超过电流环带宽之后相位延迟会迅速累积。如果想准确预测1kHz以上频段延时模型不能省。第三直流母线电压假设恒定没有把母线电压动态纳入小信号模型。工程上如果直流侧电容很大、电压环带宽很低这个简化能接受但母线上有剧烈负载波动时低频阻抗曲线会有额外偏差。3.4 为什么表达式对不上的时候别急着怀疑推导我在推导过程中遇到过一个典型现象理论表达式在低频段和中频段跟仿真扫频吻合得很好但到了2kHz附近总是多出一个相位偏移。刚开始怀疑是电流环PI参数标错反复检查发现是数字控制延时模型里少算了半个采样周期。很多文献直接用单拍延时e^{-Ts s}但实际SVPWM的调制更新方式会引入1.5拍等效延时这部分差别只有高频段才明显。所以对不上数据时先检查“模型是否漏了环节”而不是急着改参数。这是建模调试的一般规律我在后面的章节里还会展开。4. 仿真扫频验证从理论表达式到可测量数据4.1 扫频注入方案设计模型推完了必须用仿真扫频验证否则无法确认推导没出错。我的验证方案是在PCC处串联一个理想电压源作为扰动注入源输出角频率为2πf_p的微小电压信号。注意扰动电压幅值要控制在额定电压的1%到2%左右太大会激发非线性环节限幅、死区太小又会被FFT数值噪声淹没。扫频范围一般从1Hz到2000Hz低频段步长取2Hz高频段步长取10Hz分三段扫描这样既能覆盖次同步频段的关键动态又不必让仿真时间爆炸。每扫一个频点都让系统先进入稳态再开启扰动源持续运行足够的时间窗口。以50Hz基波为例一个频率点的仿真时长至少要包含20个周期的扰动信号这样FFT提取幅度和相位才有足够的频率分辨率。如果为了省时间把数据窗缩短提取出来的阻抗曲线会带上严重的频谱泄漏尤其在非基波整数倍频率点上会抖动得很厉害。4.2 从PCC数据中提取正负序阻抗扫频完成后拿到的是PCC电压和并网电流的时域波形。提取正序阻抗的步骤是对A相电压和A相电流做FFT取出扰动频率f_p处的复数分量V_p(f_p)和I_p(f_p)然后按阻抗定义计算[ Z_p(f_p) -\frac{V_p(f_p)}{I_p(f_p)} ]这里有个符号约定要小心不同文献里输出阻抗的正方向定义不一样有人用端口电压与注入电流的比值有人用含负号的形式。如果你拿自己的计算结果和论文曲线对比先确认对方用的哪种符号约定否则幅度一致但相位差180°很容易以为自己推错了。负序阻抗的提取思路类似但注入的是负序扰动并在负序旋转坐标系下面提取响应。实际操作时扰动源的三相相位关系要设置成A超前B超前C的负序排列FFT提取频率还是选f_p本身因为负序扰动在A相中映射出的频率分量同样是f_p和f_p2f1这两个。对于Z_pn这类耦合项需要在正序扰动注入时同时提取镜像频率f_p-2f1处的电流响应再与扰动电压做比。4.3 理论模型与扫频结果匹配的经验我做完后的匹配结果大概是这样的1Hz到100Hz频段理论阻抗幅值相位和扫频结果基本重合误差在1%以内100Hz到500Hz相位误差开始增大但仍在5%以内500Hz以上如果不考虑数字延时相位误差会迅速扩大到20%。考虑到次同步振荡通常发生在5Hz到50Hz之间这个精度已经能支撑稳定性分析。如果目标是分析高频谐振比如并联电容补偿引起的数百赫兹振荡就必须把延时模型做进理论推导并检查LCL滤波器的寄生参数。另外还有一个经验不要只看Bode图的对数幅频曲线来验证相位曲线更敏感。理论与扫频的幅值匹配而相位差个十几度的情况往往说明某个控制环路的相位贡献没算对。优先检查功率同步环、无功环反馈路径上的滤波时间常数以及电流环的前馈项有没有被忽略。4.4 扫频测量的适用边界扫频验证本质上是在某个稳态工作点附近做小信号激励因此它只能反映该工作点下的局部动态。如果VSG运行点从半载变到满载功率环的增益和功角都会变阻抗曲线自然跟着移动。所以做验证时要记录完整的运行条件输出有功、无功、电网电压幅值、电网频率、直流母线电压。缺少任何一个别人的实验或仿真都很难复现你的曲线。如果扰动幅值超过限幅阈值或者装置进入了限流模式那么扫频数据就不属于小信号范畴了拿它和线性模型对比没有意义。在工程测量中我会先用小扰动检查谐波畸变率确认畸变率低于设定值再保存数据如果畸变率过高减小扰动幅值重新扫。5. VSG参数如何搬动阻抗轨迹一个参数敏感性实验5.1 虚拟惯量J低频段的“动态搬移手”扫频和模型都跑通之后我做了两组参数敏感性实验。第一组是改变虚拟惯量J从0.5 kg·m²折算成标幺值约2s逐步增加到4s。阻抗曲线的变化很直观在1Hz到20Hz频段阻抗幅值曲线上的谐振峰整体向低频方向移动相位曲线穿越-180°的频率点也跟着降低。这其实就是转子运动方程的固有频率在变J越大等效机电振荡频率越低。这个现象的工程含义是如果你想抑制某个特定频率附近的振荡不能孤立地调J要看振荡频率落在阻抗轨迹的哪个区域。如果振荡发生在20Hz附近增大J可能把阻抗谐振峰移到10Hz反而给目标频率留出更大的负阻区。参数整定必须结合阻抗曲线逐点看不能只靠直觉。5.2 阻尼系数D中频段谐振峰的“吸收器”第二组实验是改阻尼系数D_p。D_p对阻抗轨迹的影响比J更线性增大D_p中低频段的相位曲线更加平滑谐振峰被明显压制阻抗幅值的尖峰变钝。原因是D_p直接出现在转子运动方程特征多项式的阻尼项里它增大后等效增加了功率同步环自身的阻尼整个电压源特性的动态响应更“收敛”。但D_p也不是越大越好。从阻抗判稳角度看过大的D_p会改变中频段阻抗相位使相位曲线在更宽的频段内接近-180°如果该频段恰好与电网感性阻抗的相位区间重叠可能出现新的稳定性风险。在我的实验里D_p从20增加到80时50Hz附近的相位余量确实下降了这提醒我在参数整定时要把D_p和J放在一起看而不是单独追求时域阶跃响应的阻尼比。5.3 电流环带宽与采样延时高频段的主宰者电流环带宽对阻抗的影响集中在200Hz以上。带宽越高电流环对电压指令的跟随越快逆变器在高频段保持“电压源低阻抗”特性的频带越宽带宽低的装置在高频段提前切换到滤波电感的感性阻抗相位曲线提前翻转。数字控制延时则主要在高频段叠加负相位使得相位曲线在接近Nyquist频率时急速下拉。这两者共同决定VSG在高频段与电网阻抗交互的方式。如果你想扩展高频段的稳定边界可以尝试提高电流环带宽、减小采样延时、采用双采样或预测控制但在工程装置上这些都会受限于开关频率和控制板性能。阻抗模型的价值就在这里先把边界算清楚再决定值不值得花钱升级硬件。5.4 无功-电压环与虚拟阻抗的附加影响无功-电压环的作用在低频段与有功环耦合在一起。它的时间常数通常比有功环更慢所以主要影响1Hz以下极低频段的阻抗特性。如果你做的是次同步频段分析10Hz到50Hz无功环的贡献相对有限但分析超低频功率振荡或者孤岛模式下的动态时无功环的参数就必须进模型。虚拟阻抗技术是目前VSG抑制振荡的常用手段它的本质是在电压指令中人为串联一个阻抗项改变端口的等效输出阻抗。在序阻抗模型里这相当于在Z_p和Z_n上直接叠加一个附加阻抗。由于虚拟阻抗的参数可以独立于物理滤波器设置它给阻抗轨迹整形提供了很大自由度但也容易造成“只压谐振峰、不解决负阻区”的表面稳定实际抗扰动能力更差。我建议用阻抗曲线逐频段检查虚拟阻抗的作用而不是只看时域波形是否衰减。6. 建模路上最容易翻车的五个坑6.1 功率计算滤波环节的隐性极点VSG的功率外环通常需要从瞬时功率计算中提取有功和无功的直流分量。常见的做法是使用低通滤波器或移动平均窗这个滤波环节会在功率反馈路径上引入一个低频极点。如果建模时图省事把功率计算近似成瞬时值那么在几十Hz频段的理论阻抗曲线会与实测产生明显偏差尤其在相位上表现为额外的滞后。对应办法是把功率滤波器的传递函数明确写在功率外环反馈链路中。如果用的是二阶低通注意阻尼比和自然角频率都要进模型如果用的是移动平均窗需要把窗函数的s域等效比如一阶惯性近似代入。仿真扫频数据如果和理论值在滤波器截止频率附近对不上优先检查这里。6.2 数字控制延时的半拍误差前面提过数字控制的等效延时经常被简化成e^{-Ts s}但实际系统从采样到PWM更新往往经历一拍半延迟一拍用于ADC采样和计算半拍用于PWM调制波更新。对于f_sw10kHz的装置1.5拍大约150微秒在1kHz处就带来54°的相位滞后。这在高频阻抗匹配分析中已经是不可忽略的误差。我的建议是写模型时直接把延时写成e^{-1.5Ts s}并且根据你所用控制器的PWM更新机制确认是不是1.5。有的控制器采用对称规则采样等效延时就是1拍有的采用异步更新可能接近2拍。延时定不准高频段的对数幅频曲线几乎不可能匹配。6.3 正负序频率配对的张冠李戴我在推导和扫频时多次被频率配对问题绊倒。正序扰动f_p和它的镜像频率f_p-2f1经常被搞混尤其是在手动写FFT提取脚本时。一个典型的错误是扫了正序80Hz扰动却从80Hz处提取电流又从60Hz处提取电压去算耦合阻抗实际镜像频率应该是50Hz-20Hz-20Hz取绝对值后对应20Hz的分量。这种张冠李戴会让Z_pn的幅值和相位完全失真。建议把频率配对规则写在代码注释里正序扰动在主频率f_p和f_p-2f1处产生响应负序扰动在主频率f_p和f_p2f1处产生响应。每个提取函数都标注是“主频率响应”还是“镜像频率响应”并用仿真波形人工核对一次FFT峰值位置再批量处理。6.4 工作点漂移导致阻抗曲线“失灵”阻抗模型和扫频数据都依赖工作点。VSG的功角、输出无功、电网电压幅值变化后功率外环的小信号增益会改变阻抗曲线会整体平移。如果扫频过程中电网电压或负载有缓慢波动FFT提取出的“稳态”分量实际混入了非小信号成分曲线就会变得毛糙。工程上有两个办法一是尽量让仿真和实验在电网模拟器稳压模式下进行并等系统完全进入稳态再启动扫频二是对每个频率点重复两次实验确认两次提取的阻抗曲线重复性良好。如果两次曲线对不上先检查工作点是否漂移再谈模型精度。6.5 限幅与非线性环节在小信号以外的“鬼影”VSG的限幅环节电压幅值限幅、电流限幅、角度变化率限幅在正常小信号扰动下不会激活所以线性模型可以忽略它们。但如果你为了加速仿真把扰动幅值设到5%以上限幅可能间歇性动作产生与扰动频率无关的谐波分量。此时FFT提取出的该频率阻抗含有大量非线性分量和线性模型对不上是正常的不是模型错误。如果你关心大扰动后的稳定性线性阻抗模型本身就不适用需要通过时域仿真或描述函数法另做分析。阻抗模型的正确使用场景是小扰动下的模式分析和参数整定这一点需要时刻记住。7. 从序阻抗到稳定性判断怎么做工程决策7.1 阻抗比判据与振荡模式定位拿到VSG输出序阻抗Z_VSG和电网阻抗Z_g之后判断稳定性的经典方法是看阻抗比Z_g/Z_VSG的Nyquist曲线是否绕过(-1,0)点。对于单输入单输出系统忽略负序耦合这个判据等价于检查Z_g和Z_VSG的幅频与相频曲线是否在某个频率处相交且相位差达到180°。扫频数据画出来以后直接找两条幅频曲线的交点再查交点处相位差比看Nyquist图更直观。我之前遇到的那次23Hz振荡用阻抗模型定位后发现正是VSG功率同步环在23Hz附近产生了一个负阻区而电网阻抗在该频率恰好呈感性两者相位叠加超过180°。定位到这一点之后整定参数就很有针对性增大阻尼D_p把负阻区压下去同时调整虚拟惯量J把谐振峰移出目标频段。7.2 负序阻抗在不对称工况下的作用三相电网发生单相接地、两相短路等不对称故障时负序阻抗就开始起作用。VSG的负序阻抗曲线和正序阻抗通常不完全相同因为无功-电压环和功率同步环对负序分量的响应增益不一样。如果你用对称的正序阻抗模型去分析不对称故障下的电流峰值会明显低估或高估导致过流保护整定错误。完整的做法是把2x2序阻抗矩阵全部纳入网络方程结合电网的序网络建立复合序网分析。工程中如果只需要粗略判断可以先用正负序阻抗标幺值的对比确定最不利工况再在时域里做精确验证。7.3 阻抗轨迹与参数整定的闭环序阻抗模型最有价值的应用是把“稳定性问题”翻译成“阻抗轨迹整形问题”。传统参数整定靠时域试错调J和D_p看到波形稳定就收工。有了阻抗模型后你可以直接把参数变化对阻抗轨迹的影响画成一组曲线目标明确让阻抗轨迹在所有可能电网阻抗条件下满足相位裕度要求。建议的做法是建立电网阻抗从弱到强的扫描矩阵对每个电网阻抗值画出阻抗比Nyquist曲线计算出最小相位裕度然后在这个裕度面里搜索参数范围。这个工作用脚本跑完全自动化比起时域试错效率高出几个量级。7.4 扩展到多机并联与复杂网络单机阻抗模型建立之后多机并联的稳定性分析只是在网络层面做阻抗叠加。如果各台VSG参数相同可以先做单机等值如果参数差异大就需要保留每台机的序阻抗模型形成节点阻抗矩阵做特征值分析。这个方向对新能源场站的并网振荡分析特别实用。另外把VSG阻抗模型和输电线路、无功补偿装置、变压器饱和特性的阻抗放在同一个网络里求解可以复现很多现场复现不了的振荡模式。这也是我后续打算继续深挖的方向——把序阻抗建模推向实际场站的数字化镜像让每次现场振荡都能在仿真环境里被准确复现、快速定位。回到最开始那次23Hz振荡后来我用阻抗模型定位问题后调整了阻尼系数和虚拟惯量再上现场测试振荡消失波形干净了不少。回头想最大的收获不是“调对参数”而是建立了一套可以复用的分析方法从物理现象到控制链路从控制链路到阻抗模型从阻抗模型到参数整定每一步都有清晰的频域视角支撑。如果你也在做VSG相关的工作建议先从仿真扫频把你的阻抗基线测出来再回头推理论模型两条腿走路能少走很多弯路。
