基于LTPV理论的MMC低频振荡机理分析与快速预估方法

发布时间:2026/6/3 18:14:01

基于LTPV理论的MMC低频振荡机理分析与快速预估方法 1. 项目概述从“黑盒”到“白盒”快速透视MMC的固有振荡在高压直流输电和新能源并网领域模块化多电平变换器因其模块化结构、低谐波输出和高电压等级等优势已成为核心的功率变换装置。然而随着MMC的大规模应用一个长期困扰工程师的“幽灵”问题——低频振荡——也日益凸显。这种振荡并非源于外部的电网扰动或复杂的控制环路而是深植于MMC拓扑结构本身的“胎里带”特性。传统的稳定性分析无论是基于状态空间模型的特征值分析还是基于DQ坐标系的阻抗扫描在面对MMC这种高维、强耦合的复杂系统时往往显得力不从心。特征值分析难以给出振荡频率的解析表达式而常规的阻抗建模则可能因为模型简化而丢失关键信息导致对振荡机理的理解停留在“知其然不知其所以然”的层面。本文要探讨的正是一种能够快速、直观地“透视”MMC内在低频振荡模式的方法。其核心在于绕开传统建模的复杂性直接基于线性时周期变量理论推导出MMC交流侧导纳的解析表达式。这个表达式是一个六阶多项式分式其分母的根直接对应了系统固有的振荡模式。这就像是为MMC这个复杂的“黑盒”系统建立了一个清晰的“白盒”等效电路模型。我们不仅能快速计算出振荡频率还能通过合成福斯特型等效电路直观地看到这些振荡模式对应的物理LC谐振回路从根本上理解MMC为何会表现出多个谐振峰以及桥臂电容在其中扮演的关键角色。对于从事MMC系统设计、稳定性分析和控制器研发的工程师而言掌握这种方法意味着在面对潜在的振荡风险时不再需要依赖耗时费力的全阶仿真或复杂的特征值计算而是能够通过几个关键参数如桥臂电感、子模块电容、调制比的简单计算快速预估系统的稳定边界为控制器的参数整定和系统优化提供直接的、定量的理论指导。2. 核心思路拆解为何LTPV建模是解开MMC振荡之谜的钥匙要理解本文方法的精妙之处首先需要厘清MMC低频振荡分析的难点所在以及传统方法的局限性。MMC的复杂性源于其分布式储能子模块电容和复杂的内部环流路径。传统的状态空间建模方法虽然系统性强但会得到一个维度极高的状态矩阵。从这个高维矩阵中提取出主导的振荡模式并进行物理意义上的解读是一项极其困难的任务。而基于DQ坐标系的阻抗建模方法在处理非对称工况或单相系统时存在局限且其模型通常难以获得简洁的解析形式更多是用于数值计算和频域扫描验证。2.1 线性时周期变量理论的核心优势线性时周期变量理论为解决上述问题提供了一条新路径。其核心思想在于承认电力电子变换器在稳态工作时其内部状态变量如电容电压、电感电流是周期时变的而非传统小信号线性化所假设的恒定直流工作点。LTPV建模通过将周期时变量在多个谐波频率上进行线性化构建了一个能够精确描述系统在基频及其旁频如±ω0, ±2ω0处动态行为的模型。对于MMC而言其桥臂电流、子模块电容电压均含有显著的二倍频纹波这正是典型的LTPV系统特征。采用LTPV方法建模能够自然地将这些周期性的工作点变化纳入考虑从而得到比传统方法更精确的端口阻抗/导纳模型。更重要的是通过合理的谐波阶数截断例如在本文中主要考虑基频和二次谐波即h1我们可以从庞大的LTPV矩阵模型中解析地推导出交流侧导纳Y(s)的显式表达式。这个表达式是一个关于复频率s的有理分式函数其极点直接对应了系统的自然振荡频率。注意谐波阶数h的选择是精度与复杂度之间的权衡。h1考虑基频和±1次谐波已能捕获主导的低频振荡模式且能获得相对简洁的六阶多项式。若追求更高精度尤其是对高频段需提高h值但多项式阶数会急剧上升解析求解将变得异常困难。在实际工程预估中h1通常是一个良好的起点。2.2 从导纳到振荡模式多项式求根的物理意义得到导纳的解析表达式Y(s) Num(s)/Den(s)后关键的一步是分析其零极点。分母多项式Den(s)的根即极点对应系统导纳的峰值阻抗的谷值也就是系统容易发生谐振的频率点。分子多项式Num(s)的根即零点对应导纳的谷值阻抗的峰值。本文的突破性贡献在于通过对模型进行合理的简化忽略电阻、假设q轴调制量为零、忽略环流抑制控制产生的调制量最终将Den(s)简化成了一个关于s²的三次方程。这使得我们可以直接使用卡尔丹公式解析地求出三个极点sd,1², sd,2², sd,3²。对其开方并考虑虚部即可得到三个固有的振荡频率。其中sd,1对应的模式最为关键其频率约为1/(2π√(L_arm * C_eq))这是一个典型的LC谐振频率直接由桥臂电感和等效子模块电容决定。另外两个模式则与调制比和基频相关。这种方法的巨大优势在于“快速”和“直观”。工程师只需将MMC的桥臂电感L_arm、子模块电容C_m、子模块数量NC_eq C_m/N、调制比m_d等几个关键参数代入公式即可直接计算出这三个潜在的振荡频率点而无需进行任何复杂的矩阵运算或特征值求解。3. MMC固有振荡模式的解析推导与物理诠释本节将深入推导过程并详细解释每个振荡模式背后的物理含义。我们从一个简化的单相等效电路出发逐步构建理解。3.1 关键假设与模型简化为了使推导可行且结果清晰我们基于工程实际做了以下关键假设这些假设在大多数MMC稳定运行工况下是合理的忽略桥臂电阻在高功率MMC中桥臂电感的感抗通常远大于电阻因此忽略R_arm对振荡频率的影响是合理的它主要影响模式的阻尼。忽略环流调制量假设环流抑制控制效果良好其产生的二次谐波调制量m2很小对交流侧谐波影响可忽略。假设q轴调制量为零MMC通常运行在高功率因数下mq近似为0。若不为零可用复数md jmq进行推广。谐波阶数截断在LTPV模型中仅考虑最高谐波阶数h1即模型矩阵的阶数为3。这牺牲了部分高频精度但极大降低了模型复杂度使解析推导成为可能。基于上述假设通过对LTPV模型进行矩阵运算和求逆最终可以得到MMC在“1”序列正序下的单相导纳Y11(s)的近似解析表达式其分母Den(s)可表示为Den(s) (4L_arm C_eq s² 1) * [64L_arm² C_eq² s⁴ (128ω0² L_arm² C_eq² 16L_arm C_eq m_d² 32L_arm C_eq)s² (64ω0⁴ L_arm² C_eq² 16ω0² L_arm C_eq m_d² - 32ω0² L_arm C_eq m_d⁴ - 4m_d² 2)]令Den(s) 0即可解出三个极点s_d,1², s_d,2², s_d,3²。其解析解如原文公式(18)所示。3.2 三个振荡模式的物理内涵主导模式 (sd,1)其对应的频率为f_dom ≈ 1 / (2π √(4 L_arm C_eq))。注意这里的等效电容是C_eq但电感是4L_arm这里需要仔细审视。从表达式s_d,1² -1/(4L_arm C_eq)来看谐振角频率ω_res 1/√(4L_arm C_eq)。这揭示了该模式本质上是桥臂电感L_arm与等效子模块电容C_eq构成的谐振。为什么是4倍这源于MMC的上下桥臂结构。在交流侧看进去上下桥臂电感串联但在特定振荡模式下电流路径可能使得等效电感表现为单个桥臂电感的倍数。该模式是最显著、通常也是频率最低的振荡模式是稳定性分析中需要首要关注的。与调制相关的模式 (sd,2 和 sd,3)这两个模式的频率表达式更为复杂同时包含基频ω0、调制比md、以及L_arm和C_eq。它们反映了MMC内部由于调制过程引入的耦合 dynamics。当调制比md变化时这两个频率会随之移动。在满调制md1和典型参数下它们的频率通常高于主导模式且更接近基频。实操心得在实际工程中当观察到MMC在非基频如几十到几百赫兹存在振荡时首先应计算这个主导模式频率。如果振荡频率与之吻合那么很可能是拓扑固有的LC谐振被激发。此时增大子模块电容C_m即减小C_eq或减小桥臂电感L_arm是移动该谐振频率、规避风险区域最直接的手段。但需注意增大电容会影响动态响应速度和成本减小电感则可能影响故障电流限流能力和谐波性能需要折中考虑。3.3 福斯特等效电路将数学映射为物理元件为了更直观地理解我们可以将推导出的阻抗函数Z_MMC(s) 1/Y11(s)进行部分分式展开并综合成福斯特I型等效电路。该电路由一个串联电感和三组并联的LC谐振支路构成如下图所示此为概念示意L_Z1 ------o---o---o--- 交流端口 | | | C_Z1 L_Z2 L_Z3 | | C_Z2 C_Z3 | | --- ---其中L_Z1 L_arm/2而其他电感电容值L_Z2, C_Z2, L_Z3, C_Z3则由前面计算出的零极点s_n,i和s_d,i通过公式(22)计算得出。这个电路模型清晰地揭示了MMC交流端口阻抗呈现多个谐振峰/谷的根源正是这三组并联的LC谐振支路在频率响应上产生了多个极点和零点。当子模块电容C_m非常大时C_Z1, C_Z2, C_Z3都趋近于开路整个等效电路退化成一个简单的电感L_arm/2这与传统的两电平VSC仅带L滤波器的输出阻抗特性一致。这从电路原理上解释了为何增大电容可以抑制MMC特有的低频振荡。4. 控制策略对振荡模式的影响机理分析MMC的固有振荡模式是其拓扑的“底色”而施加在其上的各种控制策略则会像“滤镜”一样改变这些模式的阻尼甚至频率。理解这种影响对于控制器参数设计和稳定性保障至关重要。4.1 单环构网型控制可能恶化阻尼在一些追求简单可靠的实际工程中MMC会采用单环的构网型控制例如仅采用有功功率-频率下垂控制。在这种控制架构下系统的动态主要由MMC的固有导纳Y_int和很小的下垂系数D_p决定。本文的推导表明随着下垂增益D_p的增大系统主导振荡模式对应sd,1的根轨迹会向复平面的右半平面移动。这意味着该模式的阻尼在减小甚至可能变为负阻尼引发失稳振荡。时域仿真也证实在相同的功率阶跃扰动下更大的D_p会导致更剧烈的、衰减更慢的功率振荡且振荡频率与固有模式频率吻合。关键洞见这提醒我们在为构网型MMC设计下垂控制时不能简单照搬两电平VSC或同步发电机的经验。必须考虑MMC自身固有振荡频率对控制带宽的限制。下垂系数D_p的选择必须确保系统的特征频率远离这些固有谐振点并留有足够的稳定裕度。通常这要求D_p的值非常小例如小于0.005 p.u.。4.2 电流矢量控制提供有源阻尼相比之下内环的电流矢量控制通常包含比例环节对固有振荡模式的影响是积极的。电流环的比例增益k_p在阻抗模型中等效于在原有桥臂电阻R_arm上串联了一个虚拟电阻R_v ≈ k_p。从修改后的振荡模式表达式s_d,i ≈ - (R_arm R_v)/(2L_arm) ± j√( ... )可以看出虚拟电阻R_v直接增加了振荡模式的阻尼。这使得电流控制环能够有效抑制拓扑本身固有的弱阻尼振荡。根轨迹分析显示随着k_p增大主导模式的极点向左半平面移动阻尼增强。时域仿真中更大的电流环比例增益确实带来了更平稳、无超调的电流动态响应。参数设计建议在工程中电流环比例增益的整定不仅需要考虑动态响应速度还应将其对固有振荡模式的阻尼效果作为一个积极因素来考量。在保证系统相位裕度的前提下适当增大k_p有助于提升系统对低频扰动的鲁棒性。4.3 多环控制如PQ控制引入新模态并复杂化分析当在电流内环之外再加入外层的功率控制环如PQ控制时情况变得更加复杂。外环的积分器会引入新的、频率更低的控制模式。从根轨迹可以看到随着外环PI参数增大MMC原有的固有振荡模式会向左移动阻尼改善但积分作用会引入一个接近原点的新极点这个模式可能随着增益增大而变得不稳定。在导纳伯德图上这表现为在原有谐振峰的基础上在基频附近如50Hz出现一个新的谷值。这意味着在多环控制下系统的稳定性分析不能再仅仅关注拓扑的固有模式控制引入的新模式可能成为主导失稳因素。此时基于阻抗的稳定性分析如奈奎斯特判据或基于详细模型的特征值分析变得更为必要。本文提供的快速估计方法可以作为初筛定位由拓扑引起的问题但完整分析必须结合具体的控制回路。5. 仿真与实验验证流程及关键细节理论推导需要经过严格的验证。本文通过离线仿真和控制器硬件在环实验完整地验证了所提模型的准确性和方法的实用性。5.1 阻抗扫描与模型对比验证验证的第一步是获取MMC的“真实”阻抗特性。通常采用频域扫描法在MMC交流端口注入一个幅值很小、频率可变的谐波电压扰动u_h测量对应的谐波电流响应i_h则该频率下的导纳Y i_h / u_h。遍历感兴趣的频率段即可得到导纳的频响曲线。仿真设置要点模型选择在仿真软件如MATLAB/Simulink中应使用包含所有开关细节的详细模型或经过验证的聚合模型以确保结果可信。扰动注入扰动幅值要足够小如额定电压的0.1%-1%以确保系统工作在线性区间。通常注入正序扰动进行分析。对比基准将扫描得到的导纳曲线图8、图9中的圆点与本文推导的解析模型计算曲线图8、图9中的蓝色线进行对比。如图中所示两者在主要谐振峰/谷的频率和幅值上吻合良好验证了模型的正确性。特别是模型成功预测了三个主要的谐振点。5.2 参数影响验证电容的关键作用为了验证理论分析中关于电容影响的结论可以进行参数化仿真。保持其他参数不变逐步增大子模块电容C_m例如增至5倍、10倍重新进行阻抗扫描。预期现象与结果如图11所示随着C_m增大最低频的主导谐振峰对应sd,1的频率向低频移动符合f ∝ 1/√(C_eq)的规律。其他谐振峰/谷的幅值显著减小。整个导纳曲线的形状逐渐趋近于一条斜率为-20dB/dec的直线这正是纯电感L_arm/2的阻抗特性。这直观地证明了增大子模块电容可以抑制MMC特有的多谐振特性使其阻抗特性逼近简单的两电平VSC。5.3 控制器硬件在环实验验证离线仿真之后采用控制器硬件在环进行实验是更接近实际工程验证的关键一步。如图20所示使用实时仿真器如OPAL-RT运行MMC的主电路详细模型而将实际的控制器代码运行在独立的硬件控制器中。CHIL实验关键步骤模型部署与同步将MMC电路模型编译并部署到实时仿真器FPGA核心确保仿真步长足够小如40μs以准确模拟开关细节。控制器与仿真器通过高速IO板卡进行信号交互。谐波注入与测量通过仿真器的模拟输出通道将编程好的多频次谐波电压扰动信号叠加到电网电压上注入MMC端口。同时通过模拟输入通道高速采样MMC的端口电流。数据处理对采集到的时域电压、电流数据进行快速傅里叶变换提取出各次谐波分量的幅值和相位计算得到每个频率点的导纳。结果对比将CHIL实验测量得到的导纳数据点图23中的圆圈与理论计算曲线图23中的实线绘制在同一张伯德图上。如图23所示两者高度重合这强有力地证明了本文推导的阻抗模型不仅在理论上是正确的在接近实际的硬件在环环境中也是精确的。5.4 时域动态验证最后通过时域仿真验证振荡模式在动态过程中的表现。例如对开环控制的MMC施加一个很小的正序电压扰动如0.01 p.u.。现象分析如图24所示扰动后d轴电流id会出现衰减振荡。对该振荡电流进行FFT分析会发现其主要谐波分量频率并非同步频率50Hz而是对应于导纳曲线上某个谐振峰频率的镜像频率例如若导纳在43Hz有峰值时域中可能在7Hz或93Hz出现振荡这是由于DQ变换的频移效应。这个频率与理论计算出的某个sd,i频率相符从而在时域证明了该振荡模式的存在。6. 工程应用指南、常见问题与排查思路掌握了快速估计方法后如何在工程实践中应用它来指导设计和解决问题以下是基于个人经验的总结。6.1 设计阶段参数选择与稳定性预评估在新项目设计或旧项目改造时可以遵循以下流程计算固有频率根据初步选定的L_arm, C_m, N, m_d等参数利用公式(17)(18)快速计算三个固有振荡频率f1, f2, f3。评估控制带宽检查你设计的控制器电流环、功率环、电压环的带宽是否与这些固有频率过于接近。通常要求控制带宽至少远离谐振频率一定范围例如低于最低谐振频率的1/3或高于最高谐振频率的3倍以避免相互作用引发谐振。参数调整如果发现冲突优先考虑调整子模块电容C_m。增大C_m可以降低所有谐振频率使其远离控制带宽但会增加成本和体积。也可以考虑调整桥臂电感L_arm但其设计受限于谐波抑制和故障电流限制调整空间较小。阻尼设计如果频率无法避开则需要专门设计阻尼策略。电流环的比例增益本身提供阻尼。对于构网型控制可以评估引入虚拟阻抗在控制输出前增加低通滤波或电阻性项来增加特定频段的阻尼。6.2 调试与运行阶段振荡问题排查清单当MMC系统在实际运行中出现低频振荡时可以按以下步骤排查排查步骤操作与观察可能原因与应对措施1. 现象记录记录振荡波形电流、电压、频率、幅值、发生工况功率水平、电网条件。建立问题基线。2. 计算固有频率根据实际设备参数用本文方法计算f1, f2, f3。对比振荡频率判断是否与拓扑固有模式匹配。3. 检查控制模式系统处于何种控制模式开环、电流控制、构网控制等开环/单环构网振荡很可能源于固有模式阻尼不足。电流控制检查比例增益是否过小。多环控制外环参数可能激发新模态。4. 扫描阻抗如果条件允许进行在线或离线阻抗扫描获取实际的阻抗曲线。确认谐振峰位置和幅值与理论计算对比。寻找阻抗比电网阻抗/MMC阻抗不满足奈奎斯特判据的区域。5. 参数敏感性测试在安全范围内微调关键参数如电流环Kp、下垂系数Dp观察振荡变化。若增大Kp振荡减弱说明阻尼不足。若改变Dp影响显著说明构网控制与固有模式耦合。6. 检查外部因素分析电网背景谐波、附近其他电力电子设备如风机、光伏逆变器的阻抗特性。可能是与电网或其他设备阻抗交互引发的次同步/超同步振荡需进行系统级阻抗扫描分析。6.3 常见误区与注意事项忽略谐波阶数影响本文解析解基于h1的简化。对于更高频率的振荡分析如500Hz此简化可能误差较大。此时应使用更高阶的LTPV模型进行数值计算或直接采用阻抗扫描法。混淆不同坐标系下的频率注意本文推导的振荡频率是在相域abc或αβ的。当在DQ坐标系下进行特征值分析或观察时域波形时由于坐标变换一个相域频率f会在DQ域表现为f - f0和f f0两个频率分量镜像频率。图24的FFT结果就清晰展示了这一点。过度依赖解析解解析公式提供了清晰的物理洞察和快速估算能力但在最终控制器参数确定和稳定性验证时仍建议结合详细的时域仿真和基于精确模型的频域分析如奈奎斯特判据进行综合判断。电容电压均衡的影响本文模型基于平均模型假设子模块电容电压均衡良好。在实际中如果电容电压均衡动态较慢或性能不佳可能会引入额外的低频动态影响端口阻抗特性。在分析极低频振荡时需考虑这一点。最后我想分享的一点体会是对于MMC这类复杂系统拥有一个像本文这样既能揭示物理本质、又便于快速计算的简化模型就如同在迷雾中拥有了一张简明的导航图。它可能无法标注出每一处细节但足以指引你避开最主要的险滩直达问题的核心。在实际工作中我习惯于在项目初期就用这个方法快速评估一遍系统潜在的振荡风险这常常能帮助我在设计阶段就做出更优的参数选择避免后期调试陷入被动。

相关新闻