弱电网这个话题,在我手头这几个项目里反复出现。去年做一个光伏电站的并网友好性评估,业主反馈说轻载工况下总有那么几台逆变器莫名其妙地跳闸,现场录波一看,电流波形带着明显的低频包络,频谱图上次同步频段出现了一个不小的峰。排查来排查去,最后定位到逆变器输出阻抗和电网阻抗在特定频段发生了“负阻耦合”,也就是所谓振荡。当时团队内部就靠LCL-VSC阻抗建模把问题给钉死了。这类问题如果不从阻抗的角度去理解,你只会看到波形异常,根本不知道它为什么异常、什么时候会恢复正常,更别说做针对性的阻尼改造。
这次我把这套方法完整梳理了一遍,涵盖阻抗建模的核心思路、次超同步谐振的机理、Nyquist判据验证流程,以及Simulink仿真的完整实现细节,适合正在做并网变流器控制、微电网稳定性研究,或者是做产品并网认证的工程师参考。项目标题里“次+超同步谐振”这几个字,很多人理解成单纯的谐波问题,其实它和LCL滤波器谐振、电流环控制带宽、PLL动态特性都有关联,不是一路简单的“滤波器失稳”。下面我从问题背景开始,把整条思路串起来讲。
1. 项目背景与核心问题:弱电网下LCL-VSC为何容易失稳
1.1 LCL-VSC并网系统的典型结构
VSC就是电压源型换流器,通常叫做电压源型并网逆变器。LCL-VSC指的是逆变器交流侧通过LCL滤波器并网,即逆变器侧电感L1、滤波电容Cf、网侧电感L2构成三阶滤波网络。这个拓扑在光伏、储能、风电变流器里非常常见,因为LCL对高频开关纹波的衰减能力远强于单L滤波器,用相对较小的电感量就能满足并网谐波要求。
LCL滤波器带来滤波性能的同时,也引入了自身的谐振峰。这个谐振频率等于1/2π乘以sqrt((L1+L2)/(L1·L2·Cf)),一般在开关频率的1/10到1/5之间。谐振峰如果不加阻尼,在微弱扰动下就会被激发,和电网阻抗形成串联或者并联谐振回路,产生振荡。很多做控制的人只把LCL当滤波器看待,忽略了它作为能量交换网络的本质,这是后来出现谐振问题的根源。
1.2 弱电网的定义与影响
弱电网这个概念,工程上常用短路比SCR来衡量。SCR等于并网点短路容量除以设备额定容量,数值越小,说明电网相对越弱。SCR在3以下就属于弱电网,低于2就属于极弱电网。弱电网的最直接表现是电网等效阻抗不能忽略,而且这个阻抗主要是感性的,有线路阻抗Zg=Rg+jωLg。
阻抗的存在改变了系统的闭环极点位置。强电网下逆变器看到的是接近理想电压源的母线,控制设计几乎不用考虑电网阻抗的相位影响。弱电网下电网阻抗和LCL滤波器的网侧电感串联,等效改变了LCL的参数;电网阻抗上的电压降还会影响PCC电压的相位和幅值,PLL动态和电流环的交叉耦合会按频率段被放大。简单说,强电网下你精心设计的相位裕量在弱电网下会被电网阻抗“吃掉”一部分,吃完了系统就振荡。
1.3 次超同步谐振的现象与危害
次同步和超同步谐振这两个词在老一代的电力系统教材里和发电机轴系扭振绑定在一起,但在电力电子化电力系统里,含义有所扩展。在并网逆变器的语境下,次同步振荡指的是振荡频率低于工频50Hz,通常在几赫兹到四十几赫兹范围;超同步振荡是高于工频但低于LCL谐振频率范围内的振荡,典型在50到200Hz之间。
实际录波表现是什么样呢?逆变器输出电流的幅值有低频包络,好像有人在周期性地“踩刹车、松刹车”,有功功率也跟着波动,PCC电压相位抖动。频谱分析后能看到工频两侧对称地出现旁瓣:比如主振荡频率是35Hz的次同步分量,那在65Hz附近往往也有一个对应的超同步分量。理论上这对应dq旋转坐标系下的一对共轭模态,这也是标题里“次+超同步”并列出现的原因。危害是长期振荡会触发保护跳闸、损伤滤波电容、缩短IGBT寿命,甚至会引发电网侧原本不存在的轴系扭振。
2. 阻抗建模核心原理:把逆变器“翻译”成阻抗
2.1 为什么选择阻抗建模
阻抗建模的核心思想是:把逆变器整体等效成一个从PCC看进去的谐频阻抗Zv(s),把电网等效为理想电压源串联阻抗Zg(s),然后研究这两个阻抗的匹配问题,系统稳定性由Zg/Zv的回比矩阵来判断。这个思路的优势在于它把控制参数的细节隐藏在阻抗表达式里,可以把一个几千行代码的控制程序浓缩成一个频率响应。
传统状态空间模型在小扰动下也能判稳,借助特征值分析,能精确指示哪些模态危险。但问题在于状态空间模型是白盒模型,参数一改就要重新求特征值,无法从外部实测数据里直接获得,而且物理意义不够直接。阻抗模型可以来自理论推导,也可以来自仿真扫描,甚至可以直接用阻抗分析仪在真实装置上测,三种来源的阻抗谱放在同一个Nyquist判据框架里比较,非常灵活。
这里要敲一下黑板:阻抗建模适用范围是线性时不变系统的小信号稳定性分析,所有非线性环节(死区、饱和、PWM比较器)都要先在运行点附近做小信号线性化。这意味着分析结果是局部稳定性,不能指望它预测大扰动后的暂态行为,但在设计阶段用来筛选临界参数已经完全够用。
2.2 d-q坐标系中的阻抗模型与控制环节转化
三相对称系统里,最自然的做法是Park变换把交流量变成直流量,然后在d-q坐标系下建模,这样稳态工作点就是常数,小信号模型也容易展开。VSC的数学模型是三阶主电路(L1和L2各一阶电流状态、Cf一阶电压状态)加上控制环节动态。控制环节包含电流环PI、PLL、PWM延时、有时还有功率外环,这些全部要线性化后并入阻抗模型。
具体推导时,从VSC交流侧的频域关系出发,定义PCC电压vpcc和并网电流ig在d-q轴上的扰动分量,然后列写从电流指令到PCC电压的传递函数矩阵。经过整理,阻抗Zv(s)是2×2矩阵,形式为:
Zv(s)=Zdd(s),Zdq(s),Zqd(s),Zqq(s)
矩阵里的每一个元素都受到电流环增益、PLL带宽、LCL参数的影响。控制参数在对角元素和交叉元素中都存在,这决定了d-q轴间存在耦合,不能简单当成三个独立的单相系统来看。
推导过程里最容易忽略的是PLL。PLL的q轴电压被调节到零,这个调节过程把PCC电压相位扰动引入到电流控制的参考坐标系中,等效在阻抗模型里附加了一条负反馈回路。PLL带宽越高、电网阻抗越大,这条回路在主频附近的增益就越高,最终在阻抗模型上表现为中低频段出现负电阻特性。这个负电阻是次同步振荡的源头之一,理论推导里必须显式建模PLL,这一点很多论文一笔带过,实操时却被反复验证。
2.3 从阻抗矩阵到正负序阻抗
2×2的阻抗矩阵要用广义Nyquist判据,看特征值轨迹,操作起来麻烦,工程界通常采用对称分量变换把它解耦成正负序分量。利用三相系统的对称性,正序阻抗和负序阻抗可以由d-q阻抗矩阵组合得到:
Zp(s)=Zdd(s)-jZdq(s) Zn(s)=Zdd(s)+jZdq(s)
这样原本的2×2矩阵就转化成了两个标量阻抗,在正序和负序回路里分别使用经典Nyquist判据。实际工程中,电网阻抗在正序和负序下差别不大(尤其是感性线路),因此重点分析正序回比Lp(s)=Zg_p(s)/Zp(s)就能覆盖主要稳定问题。
这个变换不仅简化了判稳过程,也让“次同步/超同步”的物理图像清晰起来。静止坐标系下一个频率为f的扰动分量,在同步旋转坐标系里出现在f-f1的频率(f1为工频50Hz)。于是次同步频率f<f1对应dq坐标系里的负频率,超同步频率f>f1对应正频率。也就是说,次同步谐振和超同步谐振本质上是dq坐标系里一对共轭模态在静止坐标系里的投影,而正负序阻抗恰好能统一表达这两类频段的耦合程度。
2.4 电网侧阻抗模型与系统开环回比
电网侧在弱电网工况下等效为电阻Rg串电感Lg,这是最基本的一阶形式。做稳定性分析时还会考虑线路对地电容,但对中低压分布式并网而言,对地电容影响较小,通常忽略,只保留感性阻抗。引入式子:
Zg(s)=Rg+sLg
其幅频特性随频率升高而增加,相频特性固定在接近90度。逆变器阻抗Zv(s)的幅频特性则受控于闭环控制,通常在低频段呈容性或阻性,在某些频段可能出现负实部,相频特性在高频段由于PWM延时而滞后。把这两个阻抗做比值L(s)=Zg(s)/Zv(s),L(s)的Nyquist轨迹是否绕过临界点(-1,j0),就刻画了整个互联系统的相对稳定性。
这里有个隐含前提:VSC侧本身和电网侧各自都必须是稳定的,然后才谈得上互联后的回比判据。如果逆变器单机在理想电压源下就不稳定,那是控制器设计优先级的问题,阻抗判据不能掩盖这个事实。所以在做互联稳定性分析之前,我总会在Simulink里把电网短路容量设置到极大值,先验证逆变器自身能不能稳定并网。
3. Simulink实现阻抗扫描的完整实操流程
3.1 主电路模型搭建要点
我用的是Simulink自带的Simscape Electrical库(旧版叫SimPowerSystems)。主电路包括:三相可编程电压源用于模拟理想电网(幅值、频率可调)、RLC串联支路用来设置线路阻抗Lg和Rg、LCL滤波器(L1、L2、Cf按设计参数填)、三相两电平IGBT桥、直流侧用电容+直流电压源替代(做阻抗分析时,直流电压动态其实会影响中低频阻抗,但先固定电压可以简化模型)。
有几个模型细节值得单独说明。第一,IGBT桥的导通压降和开关损耗在Simscape的默认模型里是理想开关,不影响小信号阻抗,但如果你要复现真实损耗特性,得换用器件级模型,这会明显拖慢仿真速度。第二,LCL滤波电容的串联等效电阻ESR不要设成0,一个几十毫欧甚至几毫欧的ESR就能给谐振峰带来一点天然阻尼,设成0会导致仿真里的谐振峰比实际情况更尖锐,容易让扫频结果在LCL谐振频率附近出现异常跳变。第三,三相可编程电压源的内阻要设成极小值,否则它会和线路阻抗串联,等效改变Lg参数。
我建议在模型里使用字母命名的常量参数,比如Lg=VariableName,这样后面做参数扫描时可以直接用Simulink的Parameter Sweep,不需要手动改模型里的Mag参数。这个习惯看着简单,实际能帮你在大批量扫短路比时节省大量时间。
3.2 控制器与PLL的实现细节
控制器我用离散域实现,采样周期Ts取50微秒,对应开关频率10kHz下的PWM周期。电流内环用典型PI调节器,输出经过限幅后进入SVPWM调制。PLL我用的是同步坐标系的SRF-PLL,PI参数按经典设计法整定:自然角频率取2π×20到2π×40 rad/s,阻尼比取0.707。
控制器的模块化要注意,尽量把PLL的dq变换和电流环的dq变换分开成两个子系统,因为PLL给出的θ角是控制器参考角度,电流环的Park变换必须用这个θ,不要用直接测量PCC电压相位得到的角度,否则就相当于给PLL动态开了“顺风耳”,把它的负电阻特性人为避开了,阻抗模型和实际物理系统就对不上。
PWM调制器在Simscape里用PWM Generator模块生成12路IGBT门极脉冲。这里有个坑,PWM模块默认模式是“内载波”,载波频率和采样周期必须一致,否则离散控制信号经过连续PWM模块会产生无法预期的相位延迟。实测下来,如果载波频率10kHz、控制频率也是10kHz,那么从电流采样到PWM输出之间的延迟大约是1.5个采样周期,这对超同步频段的阻抗相位影响很大,需要在后续数据处理时给模型结果加上等效延迟校正,或者干脆在控制器里显式加入一拍延时建模。
3.3 PCC扰动注入与扫频数据处理
阻抗扫描的核心方法叫“频率注入法”。在PCC点并联一个受控电流源,让它注入频率为fp的余弦谐波电流,因为并联电流源的等效阻抗极大,不会改变主电路拓扑,注入的功率也极小,系统仍工作在原平衡点附近。然后测量PCC电压和并网电流在频率fp处的分量,两者的比值就是需要的一个频点阻抗值,即Zp(jωp)=Vp(jωp)/Ip(jωp)。
具体到Simulink里的操作流程是这样的:
- 在PCC节点并联一个三相受控电流源,激励信号由MATLAB Function或Signal Builder生成,每相幅值相同、相位相差120度,频率为fp。
- 扰动幅值一般取额定电流的5%到8%。取大了会激发非线性效应,扫描出的阻抗幅值失真;取小了扰动信号被FFT本底噪声淹没,相频曲线会打滑。
- 每次运行先让系统在纯基波工况下稳定运行0.3秒,然后投入扰动,继续运行0.2秒。仿真时长至少覆盖5个扰动周期,比如扫描20Hz时需要0.25秒。
- 仿真结束后,把PCC电压和并网电流数据导入MATLAB工作空间,对扰动投入后的最后0.1秒数据做加窗FFT。提取fp频率处的幅值和相位,计算得到阻抗的一个频点。
假设电流环带宽设计在500Hz,你关心的频段是5Hz到1000Hz,扫频点按对数均匀分布取30个点左右。每个频点仿真0.4秒,30个点总仿真时长12秒,单核跑Simulink大概需要几分钟到十几分钟,视步长设置而定。仿真步长建议固定为1微秒或2微秒,过大的步长会在高频段引入数字阻尼,让超同步频段的谐振峰被人为压低。
3.4 模型验证:扫描结果与解析模型对拍
阻抗扫描做完,不能直接拿去画Nyquist,先要和解析模型对拍一下。解析模型就是我前面说的小信号推导表达式,在MATLAB里用频响函数计算一组频点,把扫频数据和解析数据叠加在同一张Bode图上。
这一步能帮你快速发现建模错误。我在第一次搭建模型时,扫频数据和解析模型在100Hz以上完全对不上,相位差到60度,排查后才发现PWM模块的载波频率设成了20kHz,而控制器采样周期还是50微秒,相当于控制频率只有载波的1/2,导致本该有的PWM延时相位丢失。把载波频率改回10kHz后,两条曲线就基本贴合了。
对拍验收的工程标准是幅值误差不超过3%、相位误差不超过5度,这基本能保证后续Nyquist判据分析足够可靠。如果扫描结果和解析模型在某些频段系统性偏离,优先检查控制器内部有没有限制幅的饱和环节在扰动幅值下被触发,以及PLL输入端的电压测量低通滤波有没有额外相位滞后。
4. 基于Nyquist判据的稳定性分析与裕度评估
4.1 判稳逻辑与使用前提
Nyquist判据本身是自动控制理论里最基本的频域判据:闭环系统特征方程1+L(s)=0的右半平面极点个数,等于开环传递函数L(s)的右半平面极点数,加上L(s)的Nyquist轨迹绕(-1,j0)点的圈数。只要绕点圈数等于L(s)自身的右半平面极点数,系统就是稳定的。
应用到阻抗比分析上,将与电网阻抗Zg和变换器阻抗Zv相关的回比定为L(s)=Zg(s)/Zv(s),当系统稳定时需要考察该回比的Nyquist轨迹是否穿越临界点。
使用前必须确认两点:第一,L(s)本身在右半平面无极点,这对应VSC和电网各自稳定;第二,L(s)的Nyquist轨迹必须是一整条闭环曲线,所以从离散扫频点重建连续频率特性时,必须把负频率部分也补全,通常利用共轭对称性L(-jω)=conj(L(jω)),把正频率半轴关于实轴翻转就得到了负频率轨迹,组合起来得到完整的闭环曲线。
很多教材会直接说“只要Nyquist轨迹不环绕(-1,j0)点就稳定”,那是对L(s)右半平面无极点的情况而言。如果你不校验前提,抓着一个有右半平面极点的L(s)套用这句口诀,结论会恰好相反。我在几次评审汇报里见过这种低级错误,印象很深。
4.2 基于扫频数据的回比计算与Nyquist图绘制
扫频数据的频率点是离散的,不能直接把离散点画成连线就去判断环绕状态。实际操作分三个步骤:
第一步,对每个频率点计算回比数值L(jω)=Zg(jω)/Zv(jω)。Zg是解析式,带入当前频率即可;Zv用扫频数据的幅值和相位合成复数。
第二步,在离散点之间做插值。高频段幅值和相位变化平缓,用线性插值精度足够;但在LCL谐振频率附近,阻抗幅相变化剧烈,线性插值会低估曲线的弯曲程度,稳妥的做法是在对数频率轴上做样条插值,至少每个十倍频程插出200个点。
第三步,把正频率和负频率的轨迹拼接,在MATLAB里用nyquist函数或者plot(real(L), imag(L))画出曲线。
需要注意,回比L(s)的单位不是无穷小量,它是阻抗比,无量纲。画图时要同时画单位圆(半径为1的圆)和临界点(-1,0)。工程判断我只关心临界点附近的那段轨迹,通常在穿越频率附近。穿越频率是指L(s)的相频特性穿过-180度的频率点,在该点处看幅值是否大于1;或者看幅频特性穿过0dB的频率点,在该点处看相位是否比-180度更负。
4.3 从曲线形态到稳定裕度评估
读取稳定裕度是Nyquist判据应用的落点。相位裕度PM按定义是幅值等于1的那一点,对应频率处的相位与-180度之差;幅值裕度GM是相位等于-180度的那一点,对应频率处的幅值与单位圆幅值1之比。
工程上对弱电网并网逆变器,我给自己定的验收线是相位裕度至少45度,幅值裕度不低于6dB。达不到就说明系统离失稳边界太近,电网阻抗批次差异或者温度漂移都可能把它推过去。
更敏感的表征量是谐振频率处的阻尼。在Nyquist图上,轨迹距离临界点最近的位置,对应系统最小阻尼模态,这个距离换算出来就是阻尼比,可以按ζ≈(1-M)/M的近似公式做工程估算,M是关键频率处L(s)到临界点的距离。阻尼比低于0.05时,即使系统临界稳定,振荡衰减也极其缓慢,实用上应视为不稳定,因为控制器抗扰能力不足,会使振荡持续数十个周期。
除了数值裕度,曲线的接近方向也值得观察。如果Nyquist轨迹是从第三象限以几乎平行于实轴的方式逼近临界点,说明谐振主要由纯电抗失配驱动;如果是从第二象限某个方向逼近,说明负阻尼特征明显,大概率是控制带宽和PLL共同作用的结果。这两种情况的治理手段完全不一样,前者要改滤波器参数或加无源阻尼,后者要调控制带宽或加虚拟阻尼。
5. 次超同步谐振的模拟复现与参数敏感性分析
5.1 构建典型的失稳场景
在Simulink里复现次超同步谐振,只需要把关键参数推到边界附近。我通常这样设置参数:
- 电网短路比SCR=1.5,对应Lg约12mH(按380V/10kVA基准折算);
- PLL带宽设到100Hz以上,模拟现场调锁相环只图响应快、忽视稳定性裕度的典型场景;
- 电流环带宽设到1000Hz左右,高带宽会降低相角裕度;
- LCL的阻尼系数设为0,即不进行任何有源阻尼处理。
在这个参数组合下,并网电流会在大约0.5秒后逐渐出现包络振荡,波形上看就像电流幅值按一个低频周期缓慢涨落。对电流录波做FFT,频谱在主频率50Hz两侧出现明显边带:主峰下面是30Hz附近的次同步分量和70Hz附近的超同步分量。
这个仿真现象和现场录波的特征是一致的。值得关注的是,35Hz分量和65Hz分量在数学上是同一对d-q模态在静止坐标系的镜像频率,所以你要治理次同步振荡而只盯着35Hz做滤波,就会漏掉65Hz的那个“影子模态”,治标不治本。
5.2 PLL带宽与电流环带宽对稳定性的影响
把PLL带宽从30Hz逐步增加到150Hz,每次扫一组阻抗并检查Nyquist裕度,能得到一条很清晰的规律:PLL带宽越高,中频段阻抗模型的等效负电阻效应越强,次同步频段的相位裕度下降越快。到PLL带宽超过100Hz、SCR低于2时,次同步振荡基本就躲不掉了。
电流环带宽的影响要按频段区分。电流环闭环带宽决定了电流指令到电流实际值的跟踪能力,也决定了电流环在中等频率段的相位滞后。带宽越高,动态响应越好,但在超同步频段会引入更大的相位滞后。有意思的是,提升电流环带宽对次同步低频段的阻抗影响较小,它主要是改变了LCL谐振频率附件的阻尼特性。两个带宽对应不同的失稳频段,实际现场波形可以先看主振荡频率落在哪个范围,再反推是哪一个环路的带宽超标。
参数影响可以简化为一张规律表:
| 控制参数 | 主要影响的失稳频段 | 失效方向 |
|---|---|---|
| PLL带宽过高 | 次同步频段 | 相位裕度不足,等效负电阻 |
| 电流环带宽过高 | 超同步频段、LCL谐振附近 | 相位滞后增大,谐振峰抬高 |
| 直流电压环带宽过高 | 次同步至中频 | 直流动态与交流侧功率耦合加剧 |
5.3 短路比SCR与失稳模式的逻辑关系
把SCR从10降到2、再到1.2,扫描Nyquist轨迹会发现一个渐变过程。SCR较高时,系统始终有足够的相位裕度,即使控制带宽设得偏大,也只是阻尼比减小,波形衰减变慢。SCR降到3附近,超同步频段先出现失稳迹象,因为电网感抗变大后与LCL网侧电感、滤波电容形成新的谐振网络,谐振峰向低频移动并靠近电流环带宽范围。SCR继续降到2以下,次同步频段成为主导风险,这时PLL的负电阻效应叠加到感性电网阻抗上,振荡频率在20到40Hz之间,波形上的低频包络非常醒目。
这里产生一个重要工程结论。弱电网下的失稳频率不像教科书谐振那样由单一LC决定,而是由电网阻抗和逆变器阻抗共同决定的“穿越频率”。电网阻抗变化,失稳频率就会漂移,所以现场同一批逆变器在不同并网点、不同季节(线路长度变化)会出现不同振荡频率。做稳定性设计时,光在某一个短路比点上验证远远不够,必须在SCR=1.2到10的范围里做全覆盖扫描,我个人的习惯是每0.5倍距设一个协议点,确保系统在最恶劣情况仍然有可接受的裕度。
6. 稳定性改善措施与阻尼设计
6.1 无源阻尼:简单直接但有代价
最简单粗暴的做法是在LCL滤波电容支路串联一个电阻Rf,这个电阻加上电容的阻抗特性,等效在谐振频率附近引入一个并联阻尼。Rf的取值一般按照谐振频率处电容阻抗的1/3到1倍来选择,比如Cf=10μF、谐振频率800Hz时,电容阻抗约20Ω,Rf取5到10Ω。
无源阻尼的问题在于损耗。基波电流流过电容支路产生的损耗虽然不大,但谐振频率附近的谐波能量会全部消耗在Rf上,导致滤波电容发热,而且Rf的引入会削弱LCL高频衰减斜率,EMI性能下降。实际量产产品里现在很少用纯无源阻尼,但在实验室验证控制方案时我经常用,因为它能把“阻尼问题”和“控制问题”解耦,方便先确认主电路本身的稳定性,再切换成有源阻尼验证控制算法的等效性。
6.2 电容电流反馈有源阻尼的实现与参数整定
工程上使用最广泛的是电容电流反馈有源阻尼,做法是把LCL滤波电容电流ic测量回来后,乘以反馈系数Kd,叠加到电流环的输出指令上。电流环输出的电压指令v在叠加Kd·ic后,等效在电容支路并联了一个大小为1/Kd的虚拟电阻,从而抑制谐振峰。
Kd的整定不能图省事拍脑袋。第一步用公式估算初始值:按期望阻尼比0.707设计,Kd≈2·ζ·ωr·L1,其中ωr是LCL谐振角频率。以L1=3mH、谐振频率2kHz为例,ωr约12560rad/s,算得Kd约53。第二步把它拿到仿真的Bode图上看实际阻尼效果,微调。如果Kd偏小,谐振峰的抑制不明显;如果Kd过大,电流环的高频相位滞后会加剧,换来的可能是新的超同步振荡。我见过的团队经常把Kd往大拧,理由是“阻尼越大越稳”,结果把高频段的相位裕度全部牺牲掉,出现了新的振荡模态。
实现这个反馈要注意采样噪声。电容电流是开关纹波很大的信号,直接乘Kd后叠加到电压指令上,会引入大量PWM噪声,导致电流波形毛刺很大。我通常在采样后加一阶低通滤波器,截止频率设为LCL谐振频率的10倍,这样既保留阻尼频段的增益,又滤掉高频开关噪声。
6.3 虚拟阻抗与PLL参数协调的综合调优
当系统已经出现次同步振荡、无源阻尼和电容电流反馈都无法根治时,就要考虑在控制环里注入“虚拟阻抗”。做法是在电流内环的d-q坐标系中,对电流指令的d轴和q轴分别附加一项比例增益Rvir,该虚拟电阻串入输出阻抗的等效回路,增加低频段的阻尼。
虚拟阻抗参数要避免虚增和PLL动态的耦合。虚拟阻抗的本质是在电流反馈回路中修改阻抗模型的实部,理想情况是让Zv(s)在次同步频段的实部从负变正,阻断负电阻能量供给条件。但Rvir取值过大,会导致PCC电压小扰动下电流指令被强行推离稳态,电压支撑能力变差。所以虚拟阻抗设计完成后,一定要回到Nyquist幅相曲线核查次同步频段的相位有没有被重新推回-180度。
PLL参数也是整套调优里避不开的一环。既然PLL带宽过高会诱发次同步问题,有人在现场直接把PLL带宽压到20Hz,次同步振荡是压住了,但PLL响应太慢,电网电压相位跳变时逆变器跟不上,暂态下可能出现电流尖峰甚至过流。我的经验是分两步走:先用阻抗分析判断失稳频段,再用协调整定办法把PLL带宽和电流环带宽放到同一张稳定域图上做折中,而不是单独压低PLL。
7. 常见问题与调试经验实录
7.1 仿真步长、FFT窗口与频谱泄漏
Simulink里做阻抗扫描最大的坑来自仿真步长和FFT窗口的配合。仿真步长太大,高频段的谐波响应被数字算法“平滑”掉,扫描出的超同步频段阻抗幅值会偏低;步长太小时仿真时间成倍增加,我见过有人跑一个频点要一个多小时,最终放弃扫频。
我的经验是:固定步长取PWM周期的1/50到1/100。PWM周期100微秒时,步长取2微秒到1微秒就够。FFT分析时,窗口长度必须包含整数个扰动周期,最好同时包含整数个基波周期,这样能避免基波泄漏到相邻频率分量。实际做的时候我取0.2秒数据,包含10个20Hz扰动周期和10个50Hz工频周期,加汉宁窗后频谱泄漏几乎可以忽略。
7.2 阻抗矩阵非对称与高频段数据异常
理论上对称工况下Zdd≈Zqq、Zdq≈-Zqd,但实际扫频时会发现高频段这两组数据并不完全对称。原因在于三相电路在Simulink里按ABC三相建模,离散化后的PWM调制在每个采样周期内的不对称性会造成细微差异,而扫频信号本身也是三相对称的正序注入,在高频段受PWM离散特性影响,dq元素间会有十几度的相位差。
处理这个问题的办法是接受非对称性,不强行平均。做正负序转换时,先用扫频的Zdd和Zdq原始数据分别计算Zp和Zn,再单独看正序回比的Nyquist轨迹。如果非对称程度超过可接受范围,检查一下三相电路拓扑在模型中是否严格对称,比如L2三相电感参数保持一致,PWM载波是否三相对称。三相稍微不对称的LC参数在数学上会产生三倍频耦合,混进扫频结果里。
7.3 扫频时间太长与批量扫频加速
单频点逐个注入扫描30个频点,需要十几分钟仿真时间,做参数扫描研究时(比如扫SCR从1.2到10的19个点),总时长就不太有性价比了。提速办法有两个。
其一,采用叠加扰动方式,把多个不同频率的余弦信号叠加到同一个扰动源里注入,一次仿真得到多个频点的响应。叠加信号时注意各频率的幅值要给不同的随机初始相位,避免所有频率峰值同时叠加造成扰动源输出饱和。一次注入10个频点,仿真时间可以减少约一半。
其二,充分利用Simulink的FastRestart功能。参数扫描时开启快速重启,可以避免每次参数变化都重新初始化模型,整批扫描时间可以从几小时压缩到几十分钟。这个功能我用了很久才发现它能和ParameterSweep配合,省下的时间够喝好几杯咖啡。
7.4 一个容易被忽视的环节:测量参考点
在Simulink里测PCC电压和并网电流时,测量模块的位置也会影响扫频结果。电压测量必须放在线路阻抗之后、LCL滤波器之前,这个节点才是真正意义上的PCC。很多人图省事直接测理想电压源端电压,结果注入扰动后测得的电压里没有电网阻抗上的压降动态,阻抗谱就少了一段重要相位信息。与之类似,并网电流的测量点要取网侧电感L2的电流,而不是逆变器侧电感L1的电流,两者相差一个电容支路的分流,在LCL谐振频率附近差别非常明显。
还有牵涉到控制回路的参考点:控制环里的电压反馈也必须从PCC测,不能从电容电压测,这两个点在LCL滤波器中频段差异显著,选错一个点,整个Nyquist曲线的形状就变了,而且这种错误非常隐蔽,波形正常但裕度计算结果完全不对。
说到最后,我想把调试阻抗模型时最常踩的一个坑再强调一遍。很多人拿到扫频结果的Nyquist曲线,看到轨迹贴近临界点就急着调控制器参数,折腾一圈没效果,回头才发现在模型里PLL的采样延时没有建模,控制器的数字离散化和实际代码对不上。阻抗建模这个工具本身非常好用,但它严格依赖你对控制环节细节的忠实还原,任何一个二阶低通滤波器、一拍延时、限幅环节的建模偏差,都会被Nyquist曲线放大成“稳定性恶化”的假象。我个人的经验是,每次做扫频结果和解析模型对拍时,如果高频段相位误差大于5度,先停下来检查离散化建模,而不是急着调阻尼,这条路走过几次之后你会形成条件反射。
这套方法后续还能往两个方向扩展。一是多机并联场景下的阻抗聚合与交互稳定性分析,每台逆变器不再是独立的阻抗,而是通过电网阻抗互相耦合,这时候阻抗比要从单输入单输出扩展成多变量回比矩阵,判稳逻辑还是Nyquist,但计算量上了一个台阶。二是把阻抗模型用于控制器参数自动优化,把相位裕度等指标设为目标函数,在参数空间内搜索最优解。这些方向都是在当前项目方法上的自然延伸,核心仍然是准确建模、诚实对拍、谨慎判稳。