基于LTPV理论的MMC低频振荡机理分析与快速预估方法
2026/6/17 16:59:03 网站建设 项目流程

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方法建模,能够自然地将这些周期性的工作点变化纳入考虑,从而得到比传统方法更精确的端口阻抗/导纳模型。更重要的是,通过合理的谐波阶数截断(例如,在本文中主要考虑基频和二次谐波,即h=1),我们可以从庞大的LTPV矩阵模型中,解析地推导出交流侧导纳Y(s)的显式表达式。这个表达式是一个关于复频率s的有理分式函数,其极点直接对应了系统的自然振荡频率。

注意:谐波阶数h的选择是精度与复杂度之间的权衡。h=1(考虑基频和±1次谐波)已能捕获主导的低频振荡模式,且能获得相对简洁的六阶多项式。若追求更高精度(尤其是对高频段),需提高h值,但多项式阶数会急剧上升,解析求解将变得异常困难。在实际工程预估中,h=1通常是一个良好的起点。

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、子模块数量N(C_eq = C_m/N)、调制比m_d等几个关键参数代入公式,即可直接计算出这三个潜在的振荡频率点,而无需进行任何复杂的矩阵运算或特征值求解。

3. MMC固有振荡模式的解析推导与物理诠释

本节将深入推导过程,并详细解释每个振荡模式背后的物理含义。我们从一个简化的单相等效电路出发,逐步构建理解。

3.1 关键假设与模型简化

为了使推导可行且结果清晰,我们基于工程实际做了以下关键假设,这些假设在大多数MMC稳定运行工况下是合理的:

  1. 忽略桥臂电阻:在高功率MMC中,桥臂电感的感抗通常远大于电阻,因此忽略R_arm对振荡频率的影响是合理的,它主要影响模式的阻尼。
  2. 忽略环流调制量:假设环流抑制控制效果良好,其产生的二次谐波调制量m2很小,对交流侧谐波影响可忽略。
  3. 假设q轴调制量为零:MMC通常运行在高功率因数下,mq近似为0。若不为零,可用复数md + jmq进行推广。
  4. 谐波阶数截断:在LTPV模型中,仅考虑最高谐波阶数h=1,即模型矩阵的阶数为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 三个振荡模式的物理内涵

  1. 主导模式 (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的上下桥臂结构。在交流侧看进去,上下桥臂电感串联,但在特定振荡模式下,电流路径可能使得等效电感表现为单个桥臂电感的倍数。该模式是最显著、通常也是频率最低的振荡模式,是稳定性分析中需要首要关注的。

  2. 与调制相关的模式 (sd,2 和 sd,3):这两个模式的频率表达式更为复杂,同时包含基频ω0、调制比md、以及L_arm和C_eq。它们反映了MMC内部由于调制过程引入的耦合 dynamics。当调制比md变化时,这两个频率会随之移动。在满调制(md=1)和典型参数下,它们的频率通常高于主导模式,且更接近基频。

实操心得:在实际工程中,当观察到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增大:

  1. 最低频的主导谐振峰(对应sd,1)的频率向低频移动,符合f ∝ 1/√(C_eq)的规律。
  2. 其他谐振峰/谷的幅值显著减小。
  3. 整个导纳曲线的形状逐渐趋近于一条斜率为-20dB/dec的直线,这正是纯电感(L_arm/2)的阻抗特性。这直观地证明了增大子模块电容可以抑制MMC特有的多谐振特性,使其阻抗特性逼近简单的两电平VSC

5.3 控制器硬件在环实验验证

离线仿真之后,采用控制器硬件在环进行实验,是更接近实际工程验证的关键一步。如图20所示,使用实时仿真器(如OPAL-RT)运行MMC的主电路详细模型,而将实际的控制器代码运行在独立的硬件控制器中。

CHIL实验关键步骤

  1. 模型部署与同步:将MMC电路模型编译并部署到实时仿真器FPGA核心,确保仿真步长足够小(如40μs),以准确模拟开关细节。控制器与仿真器通过高速IO板卡进行信号交互。
  2. 谐波注入与测量:通过仿真器的模拟输出通道,将编程好的多频次谐波电压扰动信号叠加到电网电压上,注入MMC端口。同时,通过模拟输入通道高速采样MMC的端口电流。
  3. 数据处理:对采集到的时域电压、电流数据进行快速傅里叶变换,提取出各次谐波分量的幅值和相位,计算得到每个频率点的导纳。
  4. 结果对比:将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 设计阶段:参数选择与稳定性预评估

在新项目设计或旧项目改造时,可以遵循以下流程:

  1. 计算固有频率:根据初步选定的L_arm, C_m, N, m_d等参数,利用公式(17)(18)快速计算三个固有振荡频率f1, f2, f3。
  2. 评估控制带宽:检查你设计的控制器(电流环、功率环、电压环)的带宽是否与这些固有频率过于接近。通常要求控制带宽至少远离谐振频率一定范围(例如,低于最低谐振频率的1/3,或高于最高谐振频率的3倍),以避免相互作用引发谐振。
  3. 参数调整:如果发现冲突,优先考虑调整子模块电容C_m。增大C_m可以降低所有谐振频率,使其远离控制带宽,但会增加成本和体积。也可以考虑调整桥臂电感L_arm,但其设计受限于谐波抑制和故障电流限制,调整空间较小。
  4. 阻尼设计:如果频率无法避开,则需要专门设计阻尼策略。电流环的比例增益本身提供阻尼。对于构网型控制,可以评估引入虚拟阻抗(在控制输出前增加低通滤波或电阻性项)来增加特定频段的阻尼。

6.2 调试与运行阶段:振荡问题排查清单

当MMC系统在实际运行中出现低频振荡时,可以按以下步骤排查:

排查步骤操作与观察可能原因与应对措施
1. 现象记录记录振荡波形(电流、电压)、频率、幅值、发生工况(功率水平、电网条件)。建立问题基线。
2. 计算固有频率根据实际设备参数,用本文方法计算f1, f2, f3。对比振荡频率,判断是否与拓扑固有模式匹配。
3. 检查控制模式系统处于何种控制模式?(开环、电流控制、构网控制等)开环/单环构网:振荡很可能源于固有模式阻尼不足。电流控制:检查比例增益是否过小。多环控制:外环参数可能激发新模态。
4. 扫描阻抗如果条件允许,进行在线或离线阻抗扫描,获取实际的阻抗曲线。确认谐振峰位置和幅值,与理论计算对比。寻找阻抗比(电网阻抗/MMC阻抗)不满足奈奎斯特判据的区域。
5. 参数敏感性测试在安全范围内,微调关键参数(如电流环Kp、下垂系数Dp),观察振荡变化。若增大Kp振荡减弱,说明阻尼不足。若改变Dp影响显著,说明构网控制与固有模式耦合。
6. 检查外部因素分析电网背景谐波、附近其他电力电子设备(如风机、光伏逆变器)的阻抗特性。可能是与电网或其他设备阻抗交互引发的次同步/超同步振荡,需进行系统级阻抗扫描分析。

6.3 常见误区与注意事项

  1. 忽略谐波阶数影响:本文解析解基于h=1的简化。对于更高频率的振荡分析(如>500Hz),此简化可能误差较大。此时应使用更高阶的LTPV模型进行数值计算,或直接采用阻抗扫描法。
  2. 混淆不同坐标系下的频率:注意,本文推导的振荡频率是在相域(abc或αβ)的。当在DQ坐标系下进行特征值分析或观察时域波形时,由于坐标变换,一个相域频率f会在DQ域表现为f - f0和f + f0两个频率分量(镜像频率)。图24的FFT结果就清晰展示了这一点。
  3. 过度依赖解析解:解析公式提供了清晰的物理洞察和快速估算能力,但在最终控制器参数确定和稳定性验证时,仍建议结合详细的时域仿真和基于精确模型的频域分析(如奈奎斯特判据)进行综合判断。
  4. 电容电压均衡的影响:本文模型基于平均模型,假设子模块电容电压均衡良好。在实际中,如果电容电压均衡动态较慢或性能不佳,可能会引入额外的低频动态,影响端口阻抗特性。在分析极低频振荡时需考虑这一点。

最后,我想分享的一点体会是,对于MMC这类复杂系统,拥有一个像本文这样既能揭示物理本质、又便于快速计算的简化模型,就如同在迷雾中拥有了一张简明的导航图。它可能无法标注出每一处细节,但足以指引你避开最主要的险滩,直达问题的核心。在实际工作中,我习惯于在项目初期就用这个方法快速评估一遍系统潜在的振荡风险,这常常能帮助我在设计阶段就做出更优的参数选择,避免后期调试陷入被动。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询