1. 串补并网前先搞清DFIG风电场次同步振荡的物理源头1.1 一个让我印象深刻的仿真事故两年前我接手一个风电并网评估项目业主很直白地提了一句线路要加串补仿真里你们得帮我看看会不会出事。当时会议室里大多数人点头应下但说实话真正能讲清楚出事是指什么的人不多。直到我在PSCAD里把双馈风机风电场经串补并网的模型跑起来看到有功功率波形在故障后开始等幅甚至增幅振荡量出来是28Hz左右的次同步分量那一刻才明白这个出事就是标题里那个绕不开的词——次同步振荡/谐振。所谓串补就是在交流输电线路上串联电容器用电容的容抗去抵消线路感抗相当于把电气距离缩短能显著提高输送能力。对风电场来说动辄上百公里的220kV或330kV送出线路若要提高输送容量串补几乎是成本最低的工程手段之一。但它的副作用也很明确电感L和串联电容C组成了一个谐振回路在低于工频的某一段频率上会呈现低阻抗甚至负阻尼一旦和双馈风机的变流器控制产生交互就可能激发持续的次同步振荡。1.2 三条容易混淆的振荡路径SSR、SSCI与SSTI很多刚接触这个领域的工程师会把次同步振荡等同于传统火电机组的次同步谐振SSR这是个需要纠正的认知。电力系统里的次同步问题大致走三条路径次同步谐振SSR汽轮发电机组轴系与串补电容之间发生机电耦合涉及到汽轮机多段转子轴系的扭振模态频率通常在5~45Hz本质是电气谐振频率与机械轴系固有扭振频率接近时发生的能量交换。这是传统火电串补研究的经典问题。次同步控制相互作用SSCIDFIG或直驱风机的变流器控制系统与串补电气谐振回路相互作用。关键点在于机械轴系基本不参与轴系扭振模态不是必要条件仅靠变流器控制器的负阻尼特性就能让次同步电流增长。次同步扭振相互作用SSTI与HVDC、SVC、STATCOM等电力电子装置的控制系统相关火电机组轴系和电力电子装置之间发生的扭振交互。双馈风机风电场经串补并网时最常见、也最让仿真人员头疼的是SSCI。它不需要轴系模态参与只要有合适的串补度、风电场规模和转子侧控制参数组合次同步电流就可能增长而且风速越低越容易出问题。这一点和火电SSR有很大不同也给建模带来了不少麻烦。1.3 谐振频率怎么算一个简洁却关键的公式要判断串补线路会把次同步振荡的中心频率放在哪个位置先要计算电气谐振频率。设输电线路本身感抗为X_L工频下串补容抗为X_C串补度k定义为k X_C / X_L对应的电气谐振频率为f_er f_0 × √k这个公式推导并不复杂。LC串联谐振时谐振角频率ω 1/√(LC)而工频下X_C/X_L 1/(ω_0²LC)代入即可得到上述形式。举一个实际算例。一条220kV送出线路长度200km单位长度感抗约为0.4Ω/km则X_L 80Ω。如果串补度取30%那么X_C 24Ω谐振频率 f_er 50×√0.3 ≈ 27.4Hz。这个频率恰好在次同步频段内也是转子侧变流器控制带宽可能覆盖到的范围所以风险很高。如果串补度提高到50%谐振频率约35.4Hz风险区间也随之移动。仿真的第一步就是先把这条频率线算出来心里有个底再动手搭模型。2. 从零搭一套能复现次同步振荡的DFIG风电场串补并网模型2.1 建模工具选型PSCAD仍是主力但别指望默认配置做这类问题我接触过PSCAD/EMTDC、MATLAB/Simulink和DigSILENT PowerFactory三种工具各有侧重。简单说工具优势短板PSCAD/EMTDC电磁暂态精度高元件库完整风电并网研究社区常用次同步频段细节捕捉准确风电场聚合和控制系统搭建需要手工接线操作繁琐MATLAB/Simulink控制系统搭建灵活自带DFIG和风机库小信号分析方便电磁暂态细节受限长线路分布参数建模不如PSCAD直观PowerFactory机电暂态与电磁暂态混合特征值分析做得好适合大规模系统次同步SSCI的电磁暂态细节需要额外配置默认风机模型不够细如果目标是复现SSCI并且观测波形细节我的首选是PSCAD。但要注意PSCAD自带的通用DFIG模型很多版本默认采用异步电机加背靠背变流器控制方式未必严格对应某款商用风机需要自己修改转子侧电流内环、功率外环和锁相环的参数才能把SSCI风险真实暴露出来。不要图省事直接用库里的demo那往往是稳定工况什么振荡都跑不出来。2.2 风机模型怎么建单机等值、聚合与撬棒保护的取舍建模第一步是确定风电场等效方式。一个50MW的风电场按5MW单机算就有10台风机按1.5MW算就有30多台每台都详细建模既不现实也没必要。工程上通常有以下几种等效路径单台等值机把整个风电场等效为一台容量等于风电场总容量的DFIG。优点是模型规模小仿真速度快缺点是忽略了机组间风速差异和链路上电压分布差异。对小规模、均匀布机的情况尚可。分群等值根据风速、机型、集电线路阻抗把风电场分成若干群每群用一台等值机代表。这种做法的精度比单台等值好能展示部分机组先进入低电压穿越、部分机组仍在正常控制的情况。详细多机模型每台机组单独建模。只适合做机理分析或验证等值结果工程仿真一般不采用。做次同步振荡仿真时我建议至少用分群等值纯单机等值会把某些风况下的振荡风险抹平。另外撬棒crowbar保护不能删掉。故障期间转子过流撬棒动作会短接转子侧变流器风机从可控状态变为感应电机不可控状态系统的次同步阻尼特性会突变。很多仿真里SSCI在故障切除后出现恰恰就是因为撬棒退回后转子侧变流器恢复控制负阻尼重新占据主导。2.3 串补线路建模别忽略MOV和线路模型的选择串补电容不是简单一个电容元件就能万事大吉。真实工程中串补站包含电容器组、氧化锌避雷器MOV、火花间隙和旁路开关。MOV的作用是限制故障或过电压时电容器两端电压而MOV导通与否直接影响串补电容在故障过程中的有效容抗进而影响振荡激发。仿真中如果完全忽略MOV三相短路期间电容承受的电压会异常偏高后续振荡特性也会失真。PSCAD里有现成的串补模型含MOV建议直接用。线路本体模型也有讲究。次同步频率通常在10~45Hz对应的波长比工频长长线路的分布参数特性仍然有影响。对300km以内的220kV线路用贝杰龙Bergeron模型或带集中参数的π型模型都能得到可接受的精度。但如果线路更长或者做特高压场景建议用频率相关模型。2.4 一套可以直接上手的模型参数参考为了让讨论落地我给出一个用于复现SSCI的等值模型参数组合这套组合我在几个项目里验证过出振荡的概率很高适合用来做机理研究或调参数测试参数项数值风电场容量100 MW单台等值机定子电压690 V机端/ 经箱变升压至35kV35kV集电线路10km电缆架空混合只算等值阻抗主变35kV/220kV100MVA漏抗12%220kV送出线路200kmr0.05Ω/kmx0.4Ω/km串补度30%无穷大系统220kV短路容量5000MVA转子侧电流环PIKp0.6Ki8标幺值基准直流母线电容0.01F按等值机折算这组参数下谐振频率约27Hz如果给线路端施加一个瞬间单相接地故障并切除故障清除后风电场有功功率往往会出现明显的次同步振荡包络。具体激发效果还和风速、电流环增益有关后面单独说。3. 复现与诊断三件套时域波形、频率扫描、特征值分析怎么配合3.1 什么样的扰动能把次同步振荡诱出来模型搭好以后不是随便设置一个负载突变就能看到振荡。次同步振荡需要系统先经历一个冲击让次同步频率的电流分量在谐振回路里被激发出来。如果扰动太弱负阻尼增长需要很长时间仿真时间被拉得很长波形上看不明显。如果扰动过强又可能直接触发过流保护一堆动作波形不可读。实务中我常用的扰动方式有三种风电场出口母线三相金属性短路持续80ms后切除。这是最暴力的方式能激发所有频率分量后续观察次同步衰减或发散大多数情况下能快速看到结果。单相接地故障故障相选在220kV出线靠近串补电容的位置持续100ms。这种方式更接近真实线路单瞬故障次同步分量以正序为主波形比较干净。风速阶跃或者桨距角阶跃比如风速从7m/s突然升到8m/s。这种方式激发较弱适合验证系统在较恶劣运行点下会不会自行失稳。我建议先跑一个三相短路故障确认系统能跑出次同步波动再换更温和的扰动做对比。如果三相短路下都没有任何次同步痕迹说明控制参数或串补度设置不对先回去检查。3.2 时域波形和FFT读数的小坑观察变量通常选风电场出口有功功率、线路A相电流、DFIG转子转速。SSCI发生时有功功率呈典型的次同步包络振荡频率肉眼就能估个大概但精确频率必须靠FFT。做FFT时有几个细节容易踩坑。第一数据窗不要太短。次同步频率在10~45Hz之间如果想分辨到1Hz至少需要1秒的时窗如果只取0.2秒的故障后数据频谱上峰值会被抹平甚至看到一堆假谱线。我在项目中习惯取故障清除后2~4秒的稳态振荡段窗函数选汉宁窗减少频谱泄漏。第二如果只看电流信号基波分量会占据绝对主导FFT图上基波峰值会压得次同步峰看不见。办法是先用带通滤波器滤掉基波或者干脆用有功功率、瞬时功率做FFT——功率信号中次同步分量相对更明显。第三读数时注意区分次同步频率f_ss和它在波形上呈现的调制频率。很多教材里说次同步电流表现为工频与次同步频率的差频调制实际功率振荡频率就是f_ss别被某些文献绕晕。3.3 阻抗频率扫描一条曲线看出负阻尼区间时域仿只能告诉你振荡发生了要回答为什么会在某个频率发生阻抗频率扫描是最直观的手段。方法理论基础不复杂在风电场并网点向系统侧注入一系列不同频率的小幅正序电流分别测量对应频率的电压响应得到该频率下的系统等值阻抗Z(f) U(f)/I(f)。如果阻抗实部R(f)在某频段为负说明该频率下系统呈现负电阻特性扰动会在这个频段增长这就是SSCI的温床。在PSCAD里实现这个扫描有两个思路。一是用自带的小信号分析插件设定起始频率、截止频率和谐波扫描步长运行时FSC自动注入二是在模型里手动放置一个受控电流源给不同频率的幅值很小的正弦指令通过测量电压电流的幅值相位计算阻抗。手动方式的好处是完全自己可控缺点是操作繁琐且容易让数值振荡污染结果。早期项目为了验证手动扫描正确性我对比过两种方式扫描结果基本一致但手动方式在高频段40Hz以上对仿真步长更敏感需要把主步长缩到20μs以下否则阻抗实部会出现虚假波动。扫描得到R(f)曲线后判断标准就一句话哪个频段R(f)0哪个频段就存在负阻尼。把曲线和时域FFT峰值放到同一张图上两者能对上是很有说服力的证据。3.4 特征值分析通过状态矩阵把会不会振荡变成数字阻抗扫描和时域仿真都属于看到现象要彻底搞清稳定裕度还需要特征值分析。思路是把DFIG、变压器、串补线路、无穷大系统和控制器线性化写成状态空间形式dΔx/dt A·Δx然后计算状态矩阵A的特征值。特征值实部为正对应不稳定模态实部为负则表示该模态衰减虚部对应振荡频率。对于次同步模态我们关心的是在10~45Hz范围内是否存在实部大于0的特征值。对完整风电场模型做详细线性化建模工作量很大但针对SSCI分析业内常用简化方法风电场用单机等值忽略机械动态和桨距控制保留转子运动方程、定转子磁链动态、变流器电流内环和锁相环动态就可以得到维数可控的状态方程。这块可以用MATLAB的Simulink线性化工具箱从非线性模型里提取也可以自己根据DFIG五阶电磁模型手动推导。如果不想碰矩阵用PowerFactory直接对详细模型做特征值计算也可以节约不少时间。我个人的建议是先通过阻抗扫描锁定问题频率再用特征值验证两者结合足够指导后续抑制参数设计了。4. 把模型当试验台串补度、风速、控制参数与振荡强度的关系4.1 串补度扫描谐振频率和负阻尼强度如何随k移动串补度是影响SSCI最直接的系统级参数。还是按前面那个200km线路模型我做了从20%到60%的扫描结果趋势非常清楚串补度k理论谐振频率时域振荡表现阻抗最小R20%22.4Hz振荡衰减约4~5s后平息正值30%27.4Hz等幅振荡持续不衰接近040%31.6Hz增幅发散约2s内明显增长明显负值50%35.4Hz强烈发散保护动作负值增大60%38.7Hz强烈发散波形难以持续负值进一步增大低串补度时谐振频率较低在这个频段转子侧变流器和锁相环的相位滞后还不太严重系统阻尼尚可。串补度提高后谐振频率接近30~40Hz正好落入电流内环带宽附近控制器的负电阻效应凸显稳定裕度急剧下降。这也解释了为什么工程上串补度往往控制在20%~30%但即便30%也不能掉以轻心——只要控制参数不太理想30%照样发散发散。4.2 风速与运行点低风速为什么反而更危险如果只考虑风机结构和串补线路直觉会认为满发时电流大、振荡能量足问题应该更严重。但造成SSCI的机制不是能量大小而是控制回路的等效电阻正负。在低风速时DFIG的转差率更大转子电流频率相对工频偏移更大转子侧变流器在次同步频率下的等效阻抗实部更容易变为负值。所以实际仿真中常看到低风速工况下SSCI风险更高比如风速6~8m/s比12m/s更容易出现振荡发散。用一组对照仿真的数据来说在串补度40%和统一控制参数下风速从12m/s降到8m/s振荡增长时间常数从0.5s缩短到0.2s左右波形包络线更陡。这也说明在做SSCI稳定性评估时不能只算额定工况至少要把低风速、中等风速、高风速三个点都跑一遍。如果只看额定风速很可能得出系统稳定的错误结论。4.3 电流内环带宽一个参数改写全局阻尼在所有控制参数里对SSCI影响最显著的是转子侧变流器电流内环的比例增益Kp。电流内环带宽越高对指令跟踪越快但在次同步频段控制响应会呈现更明显的相位滞后和等效负电阻。我做过一组参数试验在其他条件不变的条件下把Kp从0.4调到1.2Kp0.4时30%串补度下系统基本能稳定振荡缓慢衰减。Kp0.6时等幅振荡。Kp0.8时发散趋势明显约3s内电流波动越过保护限制。Kp1.2时故障后不到1s就出现大幅振荡甚至触发变频器保护停机逻辑。这个趋势说明SSCI风险高的机组控制参数设置往往是为了提高动态响应速度而调高了Kp却在次同步频段埋下了负阻尼隐患。实际风机厂商的控制器参数通常不透明仿真时需要保守考虑扫一遍Kp的变化范围再判断系统能否保持稳定。5. 抑制措施在仿真里怎么落地以及我踩过的几个坑5.1 附加次同步阻尼控制SSDC的建模要点一旦仿真确认存在次同步振荡风险下一步就是验证抑制措施。目前工程上最常用的是在DFIG转子侧变流器和网侧变流器的控制系统中附加次同步阻尼控制思路和火电机组的附加励磁阻尼控制类似从有功或者电流信号里提取次同步频段的振荡分量经过带通滤波、相位补偿、增益放大之后作为附加控制信号叠加到无功或功率指令上用控制器去抵消负阻尼。在PSCAD里实现这个附加控制通常包括这几步构造一个带通滤波器中心频率可以按阻抗扫描得到的谐振频率设置比如28Hz左右带宽控制在5~10Hz既保证能提取振荡分量又不至于把工频耦合进来。相位补偿环节采用一阶或二阶超前/滞后环节把提取出来的次同步信号相位校正到能提供正阻尼的方向。相位补偿参数可以通过特征值分析来辅助设计通常是让特征值实部往负方向移动。增益限幅增益太小抑制效果不足太大则可能引入新的高频不稳定。限幅环节用来约束附加信号幅值防止抑制器自身过调制。叠加位置常见做法是附加在转子侧变流器的无功电流指令或功率外环输出上也有方案直接叠加到转子电压指令。不同叠加位置效果差别明显建议逐一试验。我实测下来附加阻尼控制在抑制SSCI方面非常有效调整得当可以将振荡从发散转化为衰减。但它对频率适应性有要求如果串补度切换导致谐振频率偏移超过5Hz带通滤波器中心频率可能失配需要设计频率自适应算法或拓宽滤波带宽。5.2 STATCOM/SVG并联补偿由系统侧提供正阻尼除了改造风机控制本身在风电场并网点加STATCOM也是常见的工程手段。STATCOM动态无功调节能力强可以通过电压控制改变并网点的等值阻抗特性从而在次同步频段提供正阻尼。仿真中建模相对简单用PSCAD标准STATCOM模型容量一般取风电场容量10%~30%下垂系数和电压环带宽按厂家默认值开始跑观察振荡幅值和收敛速度的变化。我在一个工程项目中试过用20%容量STATCOM抑制30%串补度下的SSCI效果显著未经补偿时故障后约2s振荡发散到保护动作投入STATCOM后振荡在1~2s内明显衰减。不过STATCOM对低风速工况的改善幅度不一致需要在多个风速点反复测试确认没有削弱抑制效果的风速区间。5.3 改串补度与保护配合从源头避免谐振点串补度本身是可以和前期的风机控制参数一起优化的。如果串补度从40%降到25%理论谐振频率从31.6Hz降到25Hz可能就避开了控制负阻尼最强的频段。仿真中做串补度灵敏度分析时可以顺手确认是否存在一个安全窗口。如果系统极限制约无法调低串补度那就要靠串补电容保护本身比如MOV的触发特性和旁路间隙的动作时延通过仿真优化让故障后串补电容快速退出避免持续谐振。这部分需要和线路保护定值配合不再是单纯的控制器问题。5.4 我踩过的几个坑写在这里省得大家再踩这些经验是我在模型调试和报告编写过程中遇到的最典型问题每一条都付出过实际时间成本。第一仿真步长不能想当然。PSCAD里如果只关注次同步频段主步长50μs通常够用但一旦启用了详细变流器开关模型或者做频率扫描到40Hz以上50μs会出现明显的数字化畸变导致阻抗实部在某些频点误判为负值。建议扫描类场景把步长缩到20μs或10μs多花点算力换可信度。第二串补电容的初始充电不能忽视。直接给串补电容投上再启动暂态过程电容初始电压为零会引起巨大冲击电流掩盖真实的次同步振荡。建议仿真开始阶段先让电容旁路系统稳定后再通过定时开关投入串补等暂态过程结束后再施加故障扰动。第三锁相环不能简化成纯理想同步。有人为了让仿真稳定把PLL的带宽调得很低甚至直接用无穷大系统角度替代这会人为“抹掉”一部分SSCI负阻尼。DFIG并网SSCI恰恰与PLL在次同步频段的相位响应强相关简化PLL等于自废武功。我建议用详细PLL模型带宽按真实产品的100~200Hz范围估算。第四等值机台数太少会掩盖风电场内振荡模式差异。如果仿真目标是证明系统整体稳定单机等值够用如果目标是定位“哪几台风机更容易参与振荡”单机等值完全不行至少分两到三群按机端风速不同分组否则结论可能把风险给平均掉了。第五记录数据时间太短。很多次同步振荡在故障后几百毫秒内看不出趋势要看到等幅还是发散至少观察2~3秒。数据记录文件里如果没有预留足够仿真时间重跑一次的成本可不小。说到这其实已经把我从零搭建双馈风机风电场经串补并网次同步振荡仿真模型的主要思路、工具选型、参数设置、诊断方法和抑制措施都串起来了。这套流程我从最初只会套demo到现在能在两三天内完成一个小型风电场的SSCI风险初评中间确实靠一堆“跑不出来”的报错和“发不出来”的振荡喂出来。如果你正在搭这类模型先把谐振频率算清楚再把电流环增益和风速扫一遍看到波形再谈抑制这条路基本不会走偏。