电力系统暂态稳定分析的时域仿真法.doc

引言 暂态能量函数法是基于一个古典的力学概念发展而来的，该概念中指出：“对于一个自由的(无外力作用的)动态系统，若系统的总能量V (V (X)>0，X为系统状态量)随时间变化率恒为负，则系统总能量不断减少，直至最终达到一个最小值，即平衡状态，则此系统是稳定的”。 图9-1 滚球系统稳定原理 图9-1所示的滚球系统在无扰动时，球位于稳定平衡点(stable equilibrium point，SEP)；受扰后，小球在扰动结束时位于高度h处 (以SEP为参考点)，并具有速度v。该质量为m的小球，总能量V由动能及势能mgh (g为重力加速度)的和组成，即 若小球与壁有摩擦力，则受扰后能量在摩擦力作用下逐步减少；设小球所在容器的壁高为H (以SEP为参考点)，当小球位于壁沿上，且速度为零时(即处于不稳定平衡状态)，相应的势能为mgH，称此位置为不稳定平衡点(unstable equilibrium point，UEP)，相应的势能为系统临界能量，即 根据运动原理，我们知道，若忽略容器壁摩擦，在扰动结束时小球总能量V大于临界能量时，则小球最终将滚出容器，而失去稳定性；反之，则小球将在摩擦力作用下，能量逐步减少，最终静止于SEP。而为临界状态，显然可根据判别稳定裕度。 对于一个实际系统要解决两个关键问题：一是对于一个实际系统如何构造(定义)一个合理的暂态能量函数，它的大小应能正确地反映系统失去稳定的严重性；二是如何确定和系统临界稳定相对应的函数值，即临界能量，从而可通过对扰动结束时暂态能量函数值 (即上例中的)和临界值(即上例中的mgH)的比较来判别稳定性或确定稳定域。这种判别稳定的方法统称为暂态能量函数法(transient energy function，TEF法)。它的特点是从能量的观点来判别稳定性，而不是根据系统运动的轨迹(物理量随时间变化的曲线)来判别稳定性，从而计算量少，速度快。 直接法的最大优点：(1)能计及非线性，适应较大系统。(2)计算速度快，不必逐步积分求摇摆曲线，而是通过能量判据来判别稳定。(3)能给出稳定度。 直接法的两个较大缺点：(1)模型较简单。目前真正实用的软件采用发电机二阶经典模型，恒定阻抗负荷，尚不能计及励磁系统对稳定的作用。(2)分析结果容易偏于保守。这是因为李雅普诺夫直接法相应的稳定准则是充分条件，而不是必要条件。此外在系统很大，或受到一系列扰动(如重合闸过程)时，直接法的速度、精度较差，故目前仅用于判别第一摇摆稳定性。 单机无穷大系统的直接法暂态稳定分析 本节介绍在最简单的单机无穷大系统中如何构造暂态能量函数、确定系统临界能量，并进行暂态稳定分析。 对于图9-2中系统，若发电机采用经典二阶模型，忽略原动机及调速器动态，忽略励磁系统动态，则系统完整的标幺值数学模型为 （9-1) 式中，为转子角速度和同步速的偏差；为转子角；，为机械功率；为电磁功率，为发电机内电势及无穷大系统电压间的系统总电抗 (设电阻为零)；其中E和U为常数；M为发电机惯性时间常数(即)。 设图9-2中系统在稳态时，功角特性为；在t=0时，线路上受到三相故障扰动，功角特性变为，此时发电机加速，转子角增加，直到时，切除故障线路，功角特性变为。 图9-2 单机无穷大系统直接法分析 要求研究的问题是：如何用直接法判别故障切除后系统的第一摇摆稳定性。 显然对故障后的系统，稳定平衡点为S (对应)，不稳定平衡点为U (对应)，在这二点上，均有电磁功率平衡，即。 下面定义系统的暂态能量函数，设系统动能为(注意为与同步速之偏差，故稳态时=0) （9-2) 显然可以将式(9-1)的加速方程的二边对积分而求得故障切除时的动能，即 （9-3) 若定义系统的势能为以故障切除后系统稳定平衡点S为参考点的减速面积(反映系统吸收动能的性能)，则故障切除时的系统势能为 （9-4) 从而系统在扰动结束时总暂态能量V为 （9-5) 若将系统处于不稳定平衡点U (转子角)时，系统以S点为参考点的势能作为临界能量(此值相当于滚球系统的=mgH)，则 （9-6) 和滚球系统相似，可作稳定判别如下。 当，即图9-2中面积(A+B)<面积(B+C)，或者说面积A<面积C时，则系统第一摆稳定；反之若，则系统不稳定；时系统为临界状态。显然这和等面积准则完全一致，是一个准确的稳定判据。这里假定系统有足够的阻尼，若第一摆稳定，则以后作衰减振荡，趋于S点。 下面对此作一个简单讨论。 (1) 从上面分析可知，最关键问题是如何构造(定义)一个反映系统稳定性的暂态能量函数，以及如何正确确定系统的临