1 引言
Bergen和Hill首先提出了网络结构保持模型[1]（Network Structure Preserving Model），随后文[2]、[3]又分别在此基础上做了改进。由于建模考虑了网络结构，涉及到负荷节点，使得模型中引入了潮流约束方程。 这样，结构保持模型被表示成一组微分代数方程DAE。该模型虽然能更加准确地描述电力系统，但将传统的稳定性判定方法应用于其上时却出现了较大的困难。
将传统Controlling UEP推广至结构保持模型，除EEAC方法[4]本身所具有的特点使其不受DAE影响之外，文[5]~[13]分别利用局部分析方法及奇异摄动方法还做出了许多卓有成效的工作。其中，局部分析方法虽然理论上很完备，但由于奇异面的存在而不能获得全局的结论，而且每个代数子流形都对应于不同的ODE，使得它难以应用。奇异摄动方法,即所谓的边界层方程BLE，通过它可得到对应于DAE的BLE。在某些范围内，BLE可以近似地表示DAE的动态特性，而且对基于BLE的电力系统模型运用Controlling UEP方法能够保证持续故障轨迹与故障后系统稳定域边界相撞。虽然将DAE转化为BLE可以获得全局结论，但这并没有得到证明，况且并非任何情况下DAE都可以等价的转化为BLE。
由上可见，Controlling UEP应用于结构保持模型较经典模型要困难得多，而且准确的负荷建模一直是个难题[14]。文[15]的理论分析部分已经讨论了结构保持模型中的奇异性质以及传统Controlling UEP理论应用于该模型的难点，并就此对外部跳变行为和故障后DAE系统的稳定性进行了分析。在此基础上，本文将提出基于结构保持模型的稳定性判定方法。
2 现有方法概述及问题
由于代数流形的约束、跳变行为的存在等原因，使得传统的直接法难以应用于结构保持模型中。文[5]~[13]从两方面提出了解决方法，其分别是应用隐函数定理将DAE在代数流形上等价转化为ODE，本文称之为局部分析方法；另一种是应用奇异摄动方法，将DAE转化为BLE，本文称之为奇异摄动方法。这两种方法的共同目标是相同的，即将DAE化为传统的Controlling UEP能够应用的一组常微分方程。但它们分别存在各自的缺点。
对于局部分析方法，主要存在以下问题：
（1） 不能提出全局的结论[5,16]；
（2） 没有考虑外部跳变行为。
虽然存在以上问题，但局部分析方法在确定代数子流形上是准确的，理论上也是完备的。
S. Sastry等[12,17]首先将奇异摄动方法引入结构保持模型，随后[7,9~11,13]都在该方面做了深入的研究。对于奇异摄动方法，主要存在以下问题：
（1）DAE不能在任何情况下都等价地转化为BLE。
（2）摄动参数e 的符号确定。同样的e®0，但e <0和e <0对应的BLE轨迹完全不同。文中符号选择的原则是保证BLE与原DAE系统对应的平衡点具有相同稳定属性。
（3）对于相同的持续故障轨迹，不同的e可能使Controlling UEP不同。那么当e®0时，Controlling UEP是否趋于同一个。
文[7]、[9]建议通过奇异摄动方法将DAE系统和能量函数方法联系起来，从而获得一个全局的结论。而且BLE存在与e >0无关的能量函数容易构造[7,9]。但因BLE系统并不和DAE系统等价，它们的轨迹在某些时候也不尽相同，现有判定方法的获得仍需通过启发式的推导。
3 扩展Controlling UEP(ECUEP)方法
3.1 方法提出
为了克服上述方法存在的问题并结合它们的优点, 本节将在文[15]中定理3的基础上提出新的暂态稳定判定方法，并给出仿真算例。
定义1 BLE系统是以ODE形式表示的，给定e >0, 通过持续故障轨迹Ff(t) 都能唯一地得到Controlling UEP，记为
为方便起见，将其称作Ff(t)的e-Controlling UEP。
引理1 Ff(t)为给定的持续故障轨迹，初值Ff(tf)ÎD0。当故障持续发展时，则存在：
定义2 Ff(t)为持续故障轨迹，初值Ff(tf)ÎD0。若存在序列ek®0+，使Ff(t)主导的e-Controlling UEP都对应于同一个
