基于轨迹灵敏度分析的双馈风力发电系统暂态稳定评估研究学术报告
一、 研究作者、机构及发表信息
本研究由来自印度理工学院的 Upasana Buragohain、来自印度理工学院的 Nilanjan Senroy 以及来自印度理工学院的 Abdul Saleem Mir 共同完成。研究成果以题为《Transient Stability Assessment of a DFIG-WTG System Using Trajectory Sensitivity Analysis》的学术论文形式，发表于 IEEE Transactions on Power Systems 期刊。该期刊是电力系统领域的顶级期刊，论文于2023年12月21日在线发表，并被收录在2024年7月出版的第39卷第4期中。
二、 研究学术背景
本研究的科学领域属于电力系统稳定分析与可再生能源并网技术，具体聚焦于双馈感应发电机风力发电系统的暂态稳定性问题。随着风电在可再生能源电力系统中的占比日益增加，双馈风力发电机因其高效率、宽速度运行范围、仅需部分功率变流器、能提供无功支持等诸多优点而成为主流技术。然而，与基于同步发电机机械转子动态的传统暂态稳定性不同，DFIG-WTG系统的暂态稳定性主要受其背靠背变流器中控制器动态（尤其是锁相环）以及物理设备的共同影响，这使得其稳定性分析更为复杂。准确评估其暂态稳定边界（如临界切除时间和关键参数值）对于保障高比例风电接入下电网的安全可靠运行至关重要。
传统上，用于评估动态系统暂态稳定性的方法包括等面积准则、时域仿真、能量函数法等。近年来，轨迹灵敏度分析作为一种强有力的工具，已被广泛用于解决电力系统的各种问题，如经济负荷分配、预防性控制方案设计、电压稳定指标推导等。然而，从暂态稳定性的角度来看，TSA在风力发电系统，尤其是DFIG-WTG系统稳定性边界表征方面的应用探索仍然非常有限。现有文献在分析WTG系统时，通常未深入涉及系统内部装置和控制动态的详细影响。此外，基于等面积准则或能量函数的分析方法对于高阶系统存在困难，常需进行保守的系统近似，可能导致分析不准确且计算复杂，且通常不直接提供计算参数临界值的公式。
为此，本研究旨在填补这一空白，提出一种基于轨迹灵敏度分析理论的新颖方法，用于评估并网DFIG-WTG系统的暂态稳定性。其具体目标包括：1) 估算基于DFIG的风力发电系统的临界切除时间（Critical Clearing Time， CCT）和参数临界值；2) 识别对系统稳定性影响最大的参数。该方法强调无需对系统进行任何近似或保守假设，能够准确捕捉系统动态，并提供了直接估算可能改变系统稳定性的参数临界值的手段。
三、 研究详细工作流程
本研究的工作流程高度依赖于数学建模、推导和数值仿真，主要包含以下几个核心步骤：
1. 理论框架建立：轨迹灵敏度分析模型推导 研究首先为非线性时不变动力系统建立了一般性的微分-代数方程表示。通过将该方程组对特定参数求导，推导出了描述状态变量轨迹对参数敏感性的状态灵敏度方程组。该灵敏度方程组与原系统的微分-代数方程联立求解。为了量化稳定性裕度，研究定义了基于状态轨迹灵敏度向量L2范数的指标 (n{\alpha})。该范数值越大，表明系统越接近失稳边界。进一步，定义了其倒数的指标 (\sigma{\alpha} = 1 / \max(n{\alpha}))，在稳定边界处，(\sigma{\alpha}) 趋近于零。基于此，通过两次不同参数值下的仿真得到两个 ((\alpha, \sigma{\alpha})) 数据点，进行线性外推至 (\sigma{\alpha} = 0) 处，即可估算出该参数的临界值 (\alpha_{critical})。这一理论框架是后续所有应用的核心。
2. 研究对象建模：DFIG-WTG系统的降阶模型 为了将轨迹灵敏度分析理论应用于DFIG-WTG系统，并聚焦于并网同步稳定性问题，研究采用了降阶模型作为分析对象。鉴于DFIG-WTG完整模型阶数过高，直接应用TSA会导致计算量剧增（状态数n + 参数数α * n）。研究利用了系统在特定控制器增益和装置参数下存在的宽时间尺度分离特性，选择了一个最近开发的、主要放大锁相环同步动态的降阶模型。该模型仅包含两个状态变量（与PLL的PI控制器中间变量和锁相角相关），以及一组代数方程，其微分-代数方程具体形式由文中公式(8)-(10)给出。该模型极大地简化了计算，同时保留了评估电网同步暂态稳定性所需的核心动态。研究考虑的灵敏度分析参数集包括：锁相环比例增益 (K{p}^{pll})、积分增益 (K{i}^{pll})、电网侧变流器q轴电流 (I_{qg}) 以及电网等效电抗 (Xg)。研究的重点参数是 (I{qg})，因为研究表明相较于其他输入，风速变化对 (I_{qg}) 的影响更为显著。
3. 方法应用与案例研究：单机与多机系统分析 研究通过一系列详尽的数值仿真案例来验证所提方法的有效性和普适性。所有仿真均在MATLAB环境中进行，通过数值求解增广的微分-代数方程（包含原系统动态和灵敏度动态）来获取状态轨迹及其灵敏度。 * 单DFIG-WTG系统CCT估算：首先，将故障切除时间 (t{cl}) 视为参数α，应用上述理论框架。对降阶模型进行两次不同(t{cl})（如0.2秒和0.3秒）的仿真，计算对应的状态轨迹对(t{cl})的灵敏度范数最大值(n{t{cl}})及其倒数(\sigma{t{cl}})。在(\sigma{t{cl}}-t{cl})平面上对两点进行线性外推，与横轴（(\sigma{t{cl}}=0)）的交点即被估算为临界切除时间 (t{critical})。论文估算结果为0.5秒，并通过时域仿真（设置不同故障切除时间，观察系统状态响应是否失稳）验证了这一结果，证实当(t{cl}<0.5)秒时系统稳定，(t{cl}=0.51)秒时失去同步。 * 单DFIG-WTG系统参数临界性与排序：在固定故障条件下（电网电压暂降至0.7 p.u.，持续0.06秒），计算状态轨迹对各个参数（(K{p}^{pll}, K{i}^{pll}, I{qg}, Xg)）的灵敏度，并分别计算每个参数对应的灵敏度范数最大值 (max(n{\alpha}))。将这些最大值按升序排列，即可得到参数对系统暂态稳定影响程度的“关键性”排序。研究发现，(I_{qg}) 是最关键的参数，其次是 (X_g)。由于(Xg)（电网等效电抗）通常无法直接调节，因此调整 (I{qg}) 被认定为最有效提升系统暂态稳定边界的手段。 * 单DFIG-WTG系统参数临界值估算：针对排序中发现的关键参数，采用与估算CCT相同的外推法估算其临界值。例如，固定故障切除时间和其他参数，改变 (I{qg})，进行两次仿真得到两个((I{qg}, \sigma{I{qg}}))点，外推至(\sigma{I{qg}}=0)得到 (I{qg}) 的临界值（估算为2.2051 p.u.）。随后通过时域仿真验证，当 (I{qg}) 超过此临界值时系统失步。研究同样演示了对 (Xg) 和 (K{i}^{pll}) 临界值的估算与验证过程，并在不同 (K{p}^{pll}) 设定下重复了 (K{i}^{pll}) 临界值的估算，展示了方法的鲁棒性。 * 多DFIG-WTG系统扩展应用：为了证明方法的通用性，研究将其推广至包含2台和6台DFIG-WTG的串接系统。对于多机系统，建立了相应的微分-代数方程模型（文中公式(13)-(35)），每台机器贡献两个状态变量。研究以2机系统为例，详细展示了CCT的估算流程（估算值为0.41秒）并通过时域仿真验证（稳定极限为0.39秒）。对于6机系统，也成功应用该方法估算出其CCT（0.29秒）并进行了验证。此外，在2机系统中，同样进行了参数关键性排序，发现线路电感 (L{g1}) 的影响最大，其次是两台发电机的 (I{qg1}) 和 (I_{qg2})，这一结论与单机系统关注电流参数的趋势一致。
四、 主要研究结果
本研究在各个流程中均获得了明确且相互印证的结果： 1. 理论框架可行性：成功建立了将轨迹灵敏度分析应用于非线性动力系统暂态稳定边界表征的通用数学框架，并通过定义灵敏度范数 (n{\alpha}) 及其倒数 (\sigma{\alpha})，将灵敏度信息转化为可用于量化稳定裕度和外推临界值的实用指标。 2. CCT准确估算：对于单机、2机和6机DFIG-WTG系统，基于降阶模型和TSA方法估算出的CCT值均得到了后续时域仿真结果的验证。例如，单机系统估算CCT为0.5秒，仿真验证稳定极限略低于0.51秒；2机系统估算CCT为0.41秒，验证稳定极限为0.39秒；6机系统估算CCT为0.29秒，验证稳定极限为0.29秒。这些结果高度吻合，证明了该方法在估算不同规模DFIG系统CCT方面的准确性。 3. 参数影响力排序：通过计算各参数对应的 (max(n{\alpha}))，成功对影响系统暂态稳定的参数进行了定量排序。在单机系统中，关键性顺序为 (I{qg} > Xg > …)；在2机系统中，则为 (L{g1} > I{qg1} > I{qg2} > …)。这一排序结果直观地指出了对系统稳定性最敏感、因而在控制和设计中最需关注的参数。 4. 参数临界值确定：利用外推法，成功估算了关键参数（如 (I_{qg}, Xg, K{i}^{pll})）的临界值。例如，在特定工况下，(I{qg}) 的临界值被估算为2.2051 p.u.，仿真表明超过此值系统失步。对 (K{i}^{pll}) 在不同 (K_{p}^{pll}) 设置下的临界值估算（79.5 rad/s² 和 88.4 rad/s²）及验证，进一步表明该方法在不同运行条件下均能有效工作，具有良好的鲁棒性。 5. 多机系统适用性：研究将所提方法从单机系统无缝推广至多机系统（2机和6机），并成功完成了CCT估算和参数排序，证明了该方法对于风电场等包含多个DFIG单元的复杂系统同样有效。
这些结果之间存在清晰的逻辑递进关系：首先，基于降阶模型和TSA理论框架，可以对系统进行仿真分析；仿真的直接结果是状态轨迹及其灵敏度数据；通过对灵敏度数据进行处理（计算范数、求倒数、外推），得到了CCT和参数临界值等量化结果；这些量化结果又通过独立的时域仿真得到了反向验证，从而形成一个“方法提出-应用验证-结论确认”的完整闭环。每个步骤的结果都为下一步提供了基础或验证，共同支撑了研究结论的可靠性和方法的有效性。
五、 研究结论与价值意义
本研究提出并验证了一种基于轨迹灵敏度分析来评估DFIG-WTG系统暂态稳定性的创新方法论。研究得出结论，该方法能够： 1. 准确估算单机及多机DFIG-WTG系统的临界切除时间。 2. 有效确定系统及控制参数（如PLL控制器增益、电网侧电流参考值、电网阻抗等）的临界值，为系统安全运行范围的划定提供了直接工具。 3. 识别并排序对系统暂态稳定边界影响最大的关键参数，为控制器优化设计和稳定增强措施指明了重点方向。
本研究的科学价值在于，将轨迹灵敏度分析这一成熟工具创造性地应用于以电力电子变流器为主导的DFIG-WTG系统暂态稳定性这一复杂问题上，提供了一套无需系统近似、能精确捕捉系统动态的量化分析框架。该方法克服了传统能量函数法等在高阶系统应用中面临的困难。
其应用价值显著：对于电力系统规划者和运行人员，该方法提供了一种评估含高比例DFIG风电的电力系统暂态稳定水平、确定安全运行边界（如故障清除时间限值、控制器参数安全范围）的有效工具。对于风电设备制造商和控制系统设计师，参数排序结果有助于优化控制器设计，通过调整最关键参数（如 (I_{qg})）来最大程度地提升机组乃至整个风电场的暂态稳定性能。
六、 研究亮点
七、 其他有价值内容
研究在附录或案例研究部分，提供了所分析测试系统的详细标称参数（如表I、III、IV），包括PLL控制器增益、线路参数、电流参考值等，这使得其他研究者可以复现其研究结果。同时，论文通过多个图表（如图5-12，图13-18）直观地展示了CCT和参数临界值的外推估算过程、时域仿真验证对比以及参数排序结果，增强了论述的可读性和说服力。此外，论文对多机系统进行了通用化建模推导（第k台机器），为该方法应用于任意数量DFIG串联的风电场提供了理论公式基础。