电力系统静态安全分析

1、1 2 n概述 nWard 等值 n支路开断模拟 n直流法 n补偿法 n灵敏度分析法 n发电机开断模拟 n预想事故的自动选择 3 n随着系统总容量的增加，网络的不断扩大，系统出现故障的 可能性也日趋增加。最终导致用户供电中断。最终导致用户供电中断。 n为保证供电持续性，要求系统安全可靠。 n可靠性：在互连系统规划设计方面，当出现故障，系统保 证对负荷持续供电的能力。是一个长时间的概念。 n安全性：在互连系统的运行方面，当出现故障，保证对负 荷持续供电的能力。是时变的或瞬时性问题。 n目前的安全分析， 大部分采用确定性方法，用潮流和稳定程 序对最严重的事故情况进行大量运算。 4 n安全分析的目的

2、：提高系统安全性。 n必须从系统规划、系统调度操作、系统维修等 方面统一考虑，最终体现在系统运行状态上。 n电力系统运行状态用四种状态来描述： n安全正常状态 n不安全正常状态 n紧急状态 n恢复状态 5 n对安全的解释，在实用中更确切地用正常供电情况 下，是否能保持潮流及电压模值等在允许的范围以 内表示。 等式的约束形式：g(x)=0. 式中：x为系统运行的状态量。可以认为是功率平 衡。 6 n在具有合格电能质量的条件下，有关设备的运行状态在具有合格电能质量的条件下，有关设备的运行状态 应处于其运行限值以内，即没有过负荷。应处于其运行限值以内，即没有过负荷。 即：即：UiminUiUimax

3、 PkminPiPkmax QkminQiQkmax 也可写成：h(x) 0 n综上所述：电力系统正常运行时应同时满足等式和不 等式两种约束条件。这时处于运行的正常状态。 7 n正常状态正常状态的电力系统可分为安全正常状态与不安全 正常状态。 n已处于正常状态的电力系统，在承受一个合理的预想事故集 （contingency set）的扰动之后，如果仍不违反等约束及不 等约束，则该系统处于安全正常状态安全正常状态。 n如果运行在正常状态下的电力系统，在承受规定预想事故集 的扰动过程中，只要有一个预想事故使得系统不满足运行不 等式约束条件，就称该系统处于不安全正常状态不安全正常状态。 n预防控制预

4、防控制：使系统从不安全正常状态从不安全正常状态转变到安安 全正常状态全正常状态的控制手段。 8 n电力系统安全分析就是应用预想事故分析的方法来 预见知道系统是否存在隐患，即处于不安全正常状 态，采取相应的措施使之恢复到安全正常状态。 n静态安全分析：用来判断在发生预想事故后系统是 否会发生过负荷或电压越限等。 n暂态安全分析：判断系统是否会失稳。 9 n紧急状态紧急状态：运行在只满足等式约束条件但不满足 不等式的状态。 n持久性的紧急状态持久性的紧急状态：没有失去稳定性质，可 通过校正控制校正控制使之回到安全状态。 n稳定性的紧急状态稳定性的紧急状态：可能失去稳定的紧急状 态。通过紧急控制紧急

5、控制到恢复状态。 n紧急控制一般包括甩负荷，切机，解列控制。 10 n系统经紧急控制后回到恢复状态恢复状态，此时系统可能 不满足等式约束，而满足不等式约束，或一部分满 足约束，另一部分不满足。 n对处于恢复状态的系统，一般通过恢复控制恢复控制使之 进入正常状态。 n恢复控制一般有启动备用机组，重新并列系统等。 11 12 n包括SCADA、安全监控及其它调度管理与计划的 功能系统。 n基础：SCADA、状态估计、安全分析 n运行控制：自动发电控制、负荷控制、电压控 制、调度员培训仿真等。 n电能管理：发电计划、经济调度、负荷预测、 电能交易评估、运行规划等 13 n 应用等值方法可以大大缩小问

6、题的计算规模，系统中某些不可观 察部分也通过等值方法来处理。 n 电力系统按计算要求分研究系统和外部系统。前者要求详细计算， 后者可用等值计算来取代。 n 研究系统可分为边界系统和内部系统。 n 边界系统是指内部系统与外部系统相联系的边界点（或边界 母线）。 n 内部系统与边界系统的联络支路称为联络线。 n 外部等值方法必须保证，当研究系统内运行条件发生变化，其等 值网络分析结果应与未简化前由全系统计算分析的结果相近。 14 n互联系统可用划分成 研究系统ST和外部系 统E 两部分。 n某些文献把研究系统 分成边界系统B和内 部系统I。如右。 n还有一种，把内部系 统称为研究系统，而 边界母线

7、归并在外部 系统中。 n一般WARD等值用前 种，REI等值用后一 种。 互联系统的第一种划分 互联系统的第二种划分 15 n互联系统可用下列一组线性方程组表示 n如将电网节点分为三类：以子集I表示内部系统节点集合，子 集B为边界节点集合，子集E为外部系统节点集合。式（1）可写成 （1） （2） Y UI . 0 0 EE EEEB BB BEBBBI IBII II UI YY YYYUI YY UI . . . 16 或写成 n消去外部节点后YBB受到修正，亦即边界节点的自导纳与互导纳改 变。 （3） （4） 1 BEEE DYY EQEQ EQ YUI . 11 BEB BEEEBBBE

8、EEEBBI IBII II . UIY YIYY Y YY YY UI . . n外部系统的节点注入电流IE通过分配矩阵D被分配到边界节点上， 分配矩阵D为 17 n对线性系统来说式(3)、(4)是一个严格的等值。只要IE不变，在任何IB、 II下，由（3）求得的UB、UI都与未等值网一致。 n但在实际应用中，需要注入功率来代替注入电流，即 n则（3）可写成 （5） * . .* *1 . ()diag S IUS U 1 1 BE BEEE BE B BBBEEEEBBI IBII I I I . SS Y Y UU UYY Y YY YY U S U * . . . .* . . 18

9、（6） n则式(5)可写成 * . * ()0 0() B I diag diag U E U 1 E BBBEEEB EEQ I I diag S SUY Y U EY U U S * . * . * . . . * () 基本情况下外部系统注入 功率分配到边界节点上的 注入功率增量 19 n如果系统是在某一基本运行方式下进行等值，则外部系统注入 功率分配到边界节点上的注入功率增量值为 （7） 1 E BBEEE E diag S UY Y U * . * . () n 等值是不严格的： n 由于外部系统注入功率在边界节点上的分配与U*B有关。等值 后的边界注入功率式（6）与运行方式有关。

10、n 在非基本运行方式下，由于外部节点电压UE不同于基本情况， 而（7）却引入了基本情况下的UE ，也有误差。 20 n（1）选取一种有代表性的基本运行方式，通过潮流计算确定 全网络各节点的复电压。 n（2）选取内部系统的范围和确定边界节点，然后对下列矩阵 进行高斯消元。 目的：消去外部系统，保留边界节点，得到仅含边界节点的 外部等值导纳阵。 YBB YBEYEE-1YEB BBBE EBEE YY YY n(3)根据式（7）计算出分配到边界节点上的注入功率增量， 并将其加到边界节点原有注入上，得到边界节点的等值注入PiEQ、 QiEQ。 21 ij ijijijijjiioiji EQ i b

11、gUUggUP )sincos()()( 000020 ij ioijiijijijijji EQ i bbUgbUUQ )()()sincos( 200000 n 式中： gij+jbij为与边界节点i相连的联络线或等值支路导纳。 ij0表示边界节点i和相邻节点j之间的电压相角差。 gio+jbio为支路i侧的对地支路导纳；ji表示节点j与i相邻。 n 在已知基本运行方式下的内部与边界节点i电压模值与相角Ui0,i0 后，则PiEQ,QiEQ的另外表达方式为: n 优点：在实时情况下，外部系统运行状态变化不知，而内部和边 界节点复电压和联络线潮流，可以随时由状态估计器提供。 22 23 n

12、用等值法求解潮流时，迭代次数可能过多或完全用等值法求解潮流时，迭代次数可能过多或完全 不收敛。不收敛。 n 等值网的潮流计算可能收敛在不可行解上。等值网的潮流计算可能收敛在不可行解上。 n 潮流计算结果可能误差太大。潮流计算结果可能误差太大。 这是由于求取等值是在基本运行方式下进行的， 而在系统实时情况下，由于运行方式变化会导致 外部系统实际注入变化和参数发生变化，因此造 成潮流计算的误差。这种现象在无功功率方这种现象在无功功率方 面表现得更为突出。面表现得更为突出。 24 n 并联支路的处理 n 等值后的并联支路，代表了从边界节点看出去的外部网络对地电容 和补偿并联支路。 n 因为外部网络的

13、串联阻抗值较小，所以外部系统的并联支路有集聚 于边界节点的趋势。 n 因此：等值时尽量不用并联支路，而通过求边界的等值注入来计及 影响。考虑并联支路聚集效应。 n 等值在边界的并联支路，产生错误的并联支路响应模型。如：边界 节点电压微小变化，导致并联支路无功功率显著增加。 n 而实际外部系统某些节点电压，通常受邻近的PV节点支援，而边 界节点电压的改变，对这些节点电压的影响很小。 25 n 保留外部系统的部分PV节点 n 在等值时，如果外部系统中含有PV节点，则内部系统中发生事故开 断时，应保持外部PV节点对内部系统的无功支援。否则，等值网潮 流解算结果差。 做法： n 进行外部等值时，保留无

14、功出力大，且与内部系统电气距离小 的PV节点。 26 n 非基本运行方式下WARD等值校正： n先以内部系统实时数据作状态估计，求出边界节点的电压模值与电 压相角； n 然后以所有边界节点作为平衡节点，对基本运行方式下的外部等值 系统（由边界节点及保留的外部系统节点组成）作潮流计算。 n对保留的PV节点：有功注入为0，电压模值为给定值，相角取边界节点 相角平均值。 n 潮流计算求得的边界注入用于校正基本运行方式下的注入。 n 如果校正后注入进行状态估计时，与内部信息有较大残差，可修改 边界节点电压模值与相角，重复计算23次。 27 nYEQ的稀疏性取决于消去范围的大小。 n如1000节点150

15、0条支路系统等值成200节点。其中100个是边界节点。 等值后矩阵等值支路有： 10099/24950条 是原来的三倍。 28 n对图所示的电力系统，各条 支路的导纳和节点注入电流 在图上标出。若将系统节点 划分为内部系统节点集I 5，边界系统节点集B3， 4，外部系统节点集E 1,2，对该系统进行WARD 等值。 解：首先写出网络方程： 29 n则有 30 n求边界等值导纳矩阵： 31 n得到边界等值导纳矩阵： n求边界等值注入电流： 1 BE BBEEE IIY YI . 等值后网络方程如下： 电力系统静态安全分析是根据系统中可能发电力系统静态安全分析是根据系统中可能发 生的扰动来评定系统

16、安全性的。预想事故通常生的扰动来评定系统安全性的。预想事故通常 包括包括支路开断支路开断与与发电机开断发电机开断两类。两类。 支路开断模拟支路开断模拟就是对基本运行状态的电力系就是对基本运行状态的电力系 统，通过支路开断的计算（开断潮流）分析来统，通过支路开断的计算（开断潮流）分析来 校核其安全性。校核其安全性。 开断潮流开断潮流：指网络中的元件开断并退出运行：指网络中的元件开断并退出运行 后的潮流。开断潮流是以开断前的潮流作为初后的潮流。开断潮流是以开断前的潮流作为初 值进行计算的。值进行计算的。 三支路开断模拟三支路开断模拟 32电力系统静态安全分析 支路开断模拟或开断潮流作为在线应用时，

17、对支路开断模拟或开断潮流作为在线应用时，对 计算速度要求很高，所以需要快速求解。计算速度要求很高，所以需要快速求解。 常用的计算方法有：常用的计算方法有： 直流法直流法、 补偿法、补偿法、 灵敏度分析法灵敏度分析法， 这些方法各具特色，现分别介绍如下。这些方法各具特色，现分别介绍如下。 三支路开断模拟三支路开断模拟 33电力系统静态安全分析 1 1 直流法直流法 直流潮流模型把非线性电力系统潮流直流潮流模型把非线性电力系统潮流 问题简化为线性电路问题，从而使分析计问题简化为线性电路问题，从而使分析计 算非常方便。直流潮流模型的缺点是精确算非常方便。直流潮流模型的缺点是精确 度差，只能校验过负荷

18、，不能校验电压越度差，只能校验过负荷，不能校验电压越 界的情况。但直流潮流模型计算快，适合界的情况。但直流潮流模型计算快，适合 处理断线分析，而且便于形成用线性规划处理断线分析，而且便于形成用线性规划 求解的优化问题因此，得到了广泛的应求解的优化问题因此，得到了广泛的应 用。用。 第三节第三节 支路开断模拟支路开断模拟 34电力系统静态安全分析 第三节第三节 支路开断模拟支路开断模拟 XP 1 0 XB ijijijijij PBx 0 PB 直流潮流数学模型直流潮流数学模型 写成另一种形式写成另一种形式 其中 35电力系统静态安全分析 直流潮流的断线模型直流潮流的断线模型 应用直流潮流模型求

19、解输电系统的状应用直流潮流模型求解输电系统的状 态和支路有功潮流非常简单。而且，由于态和支路有功潮流非常简单。而且，由于 模型是线性的，故可以模型是线性的，故可以快速进行追加和开快速进行追加和开 断线路后的潮流计算断线路后的潮流计算。 原理：原网络直流潮流公式：原理：原网络直流潮流公式： 当支路（或追加）开断后，而注入功率当支路（或追加）开断后，而注入功率 P P没有变化时，直流潮流公式为：没有变化时，直流潮流公式为： 第三节第三节 支路开断模拟支路开断模拟 ()XX P XP XP 36电力系统静态安全分析 1 1）阻抗矩阵的变化）阻抗矩阵的变化 设原输电系统网络的节点阻抗矩阵为设原输电系统

20、网络的节点阻抗矩阵为 ， 支路支路 两端的节点为两端的节点为 、。这里的支路是、。这里的支路是 指两节点间各线路的并联，线路是支路中指两节点间各线路的并联，线路是支路中 的一个元件。当支路的一个元件。当支路 增加一条电抗为增加一条电抗为 的线路（称追加线路）时，形成新的网络。的线路（称追加线路）时，形成新的网络。 第三节第三节 支路开断模拟支路开断模拟 X i j k k k x X 37电力系统静态安全分析 应用应用支路追加原理支路追加原理，新网络的节点阻抗矩阵为，新网络的节点阻抗矩阵为 (3-64)(3-64) 式中式中 (3-65)(3-65) (3-66)(3-66) 第三节第三节 支

21、路开断模拟支路开断模拟 T LLLL XXX XX T LkLLkkk XXeXxe Xe 0 1 1 0 k i e j TT kkkkk XXXe e Xxe Xe 38电力系统静态安全分析 式式(3-66)(3-66)可以简写为可以简写为 (3-67)(3-67) 式中：式中： (3-68)(3-68) 由式由式(3-67)(3-67)可知节点阻抗矩阵的修正可知节点阻抗矩阵的修正 量为量为 (3-70)(3-70) 第三节第三节 支路开断模拟支路开断模拟 XeXeXX T kkk 1 T kkkk xe Xe T kkk XXXXe e X 2 2）状态量的变化）状态量的变化 39电力系

22、统静态安全分析 在节点注入功率不变的情况下，可以直接在节点注入功率不变的情况下，可以直接 得到追加线路得到追加线路 后状态向量的增量后状态向量的增量 (3-71)(3-71) 3) 3)追加线路后的状态向量追加线路后的状态向量 第三节第三节 支路开断模拟支路开断模拟 k TT kkkkkk XPXe e XPXe e T kkk Xe e 40电力系统静态安全分析 当网络断开支路当网络断开支路 时只要将时只要将 换为换为 ， 以上公式同样适用。必须指出，当网络开以上公式同样适用。必须指出，当网络开 断支路断支路 使系统解列时，新的阻抗矩阵使系统解列时，新的阻抗矩阵 不存在，这时式不存在，这时式

23、(3-68)(3-68)中的中的 为无穷大，为无穷大， 因此，应用直流潮流模型可以方便地找出因此，应用直流潮流模型可以方便地找出 网络中那些开断后引起系统解列的线路，网络中那些开断后引起系统解列的线路， 对于这些线路不能直接进行断线分析。对于这些线路不能直接进行断线分析。 第三节第三节 支路开断模拟支路开断模拟 k k x k x k X k 41电力系统静态安全分析 42 n补偿法：将支路开断视为该支路未被断补偿法：将支路开断视为该支路未被断 开，而在其两端节点处引入某一待求的开，而在其两端节点处引入某一待求的 补偿电流，以此来模拟支路开断的影响补偿电流，以此来模拟支路开断的影响。 n特点：

24、不必修改导纳矩阵，可以用原特点：不必修改导纳矩阵，可以用原 来的因子表来解算网络的状态。来的因子表来解算网络的状态。 43 图图3-6 3-6 电力网络发生支路变化时的等效电路电力网络发生支路变化时的等效电路 N 1 I 1 I 2 I 2 I iij II jij II n I n I j I ij I ij Z 1 U 2 U j U n U i U 44 n用原网络的因子表对 进行消去回代运算所求出的节点电压向 量，就是待求的发生支路开断后的节点电压向量 ，且 n当网络节点i、j之间发生支路开断，可以等效地认为在i、j节点间并 联了一个追加的支路阻抗Zij，其数值等于被断开支路阻抗的负值

25、。 . I 0 . ( ) I 12 . , T iij jijn II IIIIII 0 12 ( ) . , T ijn IIIIII 1 . UYI I . U . I . ij I n这时流入原网络的注入电流将由 变成 n目前关键问题在于要求出追加支路Zij上通过的电流 ，从而求 得 。 45 图3-6 46 n对于线性网络，可以应用迭加原理把图3-6（a）分成两个网络即 图3-6（b）和3-6（c）。这时待求的节点电压 也可看成两个部 分 . U 01 UUU . . ( )( ) 0 U . ( ) 010 UYI . ( )( ) n式中： 相当于没有追加支路情况下的各节点电压，

26、这个向量可 以用原网络的因子表求出，即: 47 n若假定 ，则有 。于是由（1）就可求出当 为 单位电流时，网络各节点电压U(ij)，即 n 是向原网络注入电流向量 时求出的，其值为 . (1) U . (1) I 111 . ()() UYI 1 00 0110 ji T ijij ji T ij ij ij II I I I I ( ) ( ) , , , 1 ij I (1)( )ij II ( )1( )ijij UY I . ij I (3) (2) (1) 48 n这个电源的等值内阻抗ZT可以用其它节点的注入电流为零，仅在i、j 点分别通入正、负单位电流后，在i、j点产生的电压差来

27、表示。由式 （3）求得 后，便可求得 n应用等效发电机原理，如果把图3-7所示电路上的i、j节点间的整个系 统看成是Zij的等效电源，其空载电压就是 . ij I 00 . . ( )( ) ij EUU ( )ij U . )( . )(ij j ij iT UUZ ijTij ZZZ (4) (5) (6) ZT亦即是从i、j节点看进去的输入阻抗 令： 49 n通过等值电路图3-7可见，利用式51、52可求出 0000( )( )( )( ) . . ( )( ) ijij ij ijij ij ijij UUUU I Z ZUU 11111 010 ijijij ijijij ij ij

28、 III I UY IYIY IU UUUUU . ( )( )( )( )( ) ( )( )( )( ) . ij I 由上式就可求得支路开断后的节点电压向量。 (8) 求得 之后，由式48即可求得 (9) (7) 50 n对快速解耦潮流算法的修正方程式，可以看成是以B及B”作为 “导纳矩阵”的节点方程式，其注入电流分别为P/U及Q/U， 而待求量为及U。这样就可完全套用以上计算过程进行迭代 计算。 (10) n此时，图3-7中的追加支路Zij为 Zij=1/Bij n当开断元件是变压器时，式49中的I(ij)应改写为 I(ij)=0,nT,0,.,1,0,0T 式中：nT为在i侧的非标准

29、变比。 51 n这时式51、52、54也可改写成 00 0 00 . . ( )( ) . ( )( ) . ( )( ) () Tij ij TTij ijTijij ijTij En UU Zn UU In UUZ ZZZ 其中 n在式10中实际上只表示了开断不接地支路，而输电线或非标准变压 器的开断还应考虑接地支路的同时开断。实际应用时忽略接地支路 的影响。 52 n直流法快速，但只能进行有功潮流，没有考虑电压无功问题。 n补偿法要反复迭代。否则结果误差大，特别是电压无功误差 大。 n断线分析的灵敏度法：以节点功率增量模拟断线的影响。 ij ijijijijjii BGUUPsincos

30、 ij ijijijijjii BGUUQcossin (i=1,2,3,n) 节点功率方程： （略） 53 直流法求发电机的开断潮流：直流法求发电机的开断潮流： 当某节点发电机开断，而其它节点注入功率当某节点发电机开断，而其它节点注入功率P P没有没有 变化时，直流潮流公式为：变化时，直流潮流公式为： 以上方法实际上是假定失去的发电机功率全部由平衡以上方法实际上是假定失去的发电机功率全部由平衡 节点的发电机来平衡，而实际并非如此。节点的发电机来平衡，而实际并非如此。 X PP() X P XP 54电力系统静态安全分析 55 n计及电力系统频率特性的静态频率特性法 n发电机开断时，内部系统、

31、外部系统、联络线功率均作 相应调整。 n发电机的频率响应特性FRC (Frequency Response Characteristics)及边界节点上的等值频率响应特性是求 解这些功率变化的依据。 56 n发电机开断模拟的数学模型，通常整个变化过程分为四个时段。 n时段1：电磁暂态过程。系统的暂态潮流是按网络阻抗与机组暂 态电抗来分布，由于系统电磁储能容量很小，暂态过程在数毫 秒内即被阻尼。 n时段2：机械暂态过程。发电机的反应过程决定于机组的惯性， 有功出力的变化是由发电机旋转部分的转动惯量来决定。 n时段3：调速器动作过程。发电机间功率分配的变化是由FRC特 性来决定的。 n时段4：自动

32、发电控制。在一个控制区域内的发电机按自动发电 控制装置（AGC）的整定值进行调节。 57 n对于在线发电机开断模拟，快速反应的时段1、2不予考虑。 n时段3的行为是静态安全分析所需研究的部分，此时各台发电机的 功率变化可以用它的FRC与系统的FRC之间的比例关系来确定。 时间：几秒几十秒。 n时段4中AGC的作用：是通过二次频率调整来消除静态频率偏差， 此外还可控制联络线的功率来调整互联系统间的静态频率偏差。 58 n当系统中有发电机开断时，全系统的静态有功响应是根据调速系 统一次调节所达到的稳定状态来确定。 59 n发电机组在设定的运行点PGi0处的特性可用调差系数RGi来表示 其倒数为 n

33、式中：KGi “发电机组的FRC”，表示了当系统出现f的频率变化 时，机组i将给出怎样的有功功率增量响应。 Gi Gi f R P Gi Gi R K 1 Hz/MW MW/Hz 60 n负荷的频率响应特性：系统频率的变化引起负荷功率的变化。频率 变化幅度不大时，此特性可认为是线性的，以KLi表示。即当频率 下降时，母线i的负荷需量作线性降低。 n母线i的总有功功率变化（总响应），将受KGi和KLi的影响，并用母 线FRC：Ki 表示： n对一个具有n母线的电力系统，在每一母线都接有发电机及负荷的 假定下，系统的总频率特性将是每一处母线频率特性的总和，当节 点数为n时，系统的FRC为 11 n

34、n i si ii P KK f () GiLii iGiLi PPP KKK ff (90) (91) 61 n设在母线k处开断一台机组，而使系统丧失 的有功功率。在此 以后，由于整个系统频率降低了f，以至于系统中所有母线， 都将按各自的Ki作出响应； n假定丧失的 将全部被其它母线的有功增量和母线k的负荷有功 增量所补偿。因此，当系统丧失 的有功功率时，其它任一 母线i的净有功功率变化为： l G P l G P 1 12 - = 1 ()() ()/ () () ll iGiG iln G jLk j j k l iG knGk Pf PPf P P P PPf f Pf P PPPPP

35、 f l G P (92) 62 注意 92是由 的关系得出 LkkGk PPP ()ik l i iG l sGk K PP KK (93) n式92可写成： 上式子表示其它母线之净有功功率变化 63 n母线K失去的净有功功率则为： () l ll GLk kGkLkGLkG l G PP PPPPPP P 11 nn l GjLkjGk jj j k PPPPP (94) 1 1 1 1 () () n jGkLk jll Sk kGGn l SGk jGk j PPP fKK PPP KK PP f (95) n因为，已经假定 n所以式94有 64 n若以向量形式来表示所有节点的功率增量

36、方程式，则可以定义一个向 量H，其中的元素 hi=0 (当ik时） hi=Ks (当i = k时） l G l sG P KK PKH (96) l G s P K PH(97) n从而，式93,95可以写成 n对于大型电力系统来说，因为KsKGl，于是上式可写成 65 n在实际电力系统中，由于KLiKGi，于是有 n当各节点的注入功率变化Pi可以仅用发电机的出力变化亦即PGi 来表示，则式97可以写成 n i GiGss Gii KKK KK 1 HKP G Gs l G G K P (98) 66 n用解耦潮流算法进行交流潮流计算，可以求得发电机k开断后系统中 较精确的潮流分布。有功功率变

37、化关系写成 UBUPdiagdiag UBUAdiagdiag PA l G l sG P KK KHA (99) (100) (100) n将式96代入100的左边，可得 则式99可写成 取 式中：B是直流潮流矩阵；U是扰动前的电压模值向量。 67 n当电力系统进行外部等值时，由于外部系统的发电机也承担了有功功 率的调节任务，因此必须求出外部系统的等值FRC。 n若将全系统的节点分为三类，E为外部系统节点，B为边界节点，I为 内部系统节点，则式100可分解为以下形式 I B E IIIB BIBBBE EBEE I B E AA AAA AA P P P 0 0 I B IIIB BIBB

38、I B AA AA P P * n消去外部系统部分后，则有 68 n上式中：A*BB为外部系统等值后，相应于边界节点的系数矩阵；P*B 则为等值后边界节点的有功功率注入增量。 *1 *1 BBBBBEEEEB BBBEEEE AAAAA PPAAP 1 ll GG BBBEEEEll sGsG PP KKKK KKA A K * (101) (102) n若把式101中的PB*、PB及PE分别用边界节点的等值FRC（即 KB*）、边界节点的FRC（即KB）及外部节点的FRC（即KE）与 PGl/(Ks-KGl)的乘积表示，则式101可写成 其中 69 n于是得 (103) (104) n上面两

39、式表示外部节点FRC（即KE）与边界节点上的等值FRC（即 KB*）间的关系。KE是按ABEAEE-1的关系分配到边界节点上的。令 JE=AEE-1KE，则式104可写成 式中：JE可以由AEE三角分解后通过前代回代来求出。 n由式100可得 *1* BBEEEBEBB KKKAAKK EEEBEB KAAK 1* EBEB JAK * EEEEEE EBEBBE diagdiag diagdiag UBUA UBUA 70 n于是式103中的ABEAEE-1可写成 n代入式103得 n在实时情况下，外部系统电压模值UE是未知的。可取UE=U01p.u.， 于是上式可写成 n在上式中只用到B的

40、有关元素及边界节点电压。因此求出等值FRC， 亦即边界节点响应的总和。 1 11 BEEEBBEEEE diagdiag A AUBBU 1 *1 BBBBEEEEE diagdiag KKUBBUK 1 * 0 1 BBBBEEEE diag U KKUBBK 71 72 n预想事故的自动选择预想事故的自动选择ACS (Automatic Contingency ACS (Automatic Contingency Selection)Selection)：就是在实时条件下利用电力系统实时信息，：就是在实时条件下利用电力系统实时信息， 自动选出那些会引起支路潮流过载、电压违限等危及系自动选出

41、那些会引起支路潮流过载、电压违限等危及系 统安全运行的预想事故，并用行为指标来表示它对系统统安全运行的预想事故，并用行为指标来表示它对系统 造成的危害严重程度，按其顺序排队给出一览表。造成的危害严重程度，按其顺序排队给出一览表。 n预想事故自动选择需要一种快速的开断模拟算法，并在预想事故自动选择需要一种快速的开断模拟算法，并在 精度上能够满足排队的要求，亦即是能够剔除不起作用精度上能够满足排队的要求，亦即是能够剔除不起作用 的预想事故，并将起作用的预想事故按严重程度排队。的预想事故，并将起作用的预想事故按严重程度排队。 73 1）确定预想事故集 2）预想事故选择（ACS） 运用开断模拟算法，计

42、算预想事故的行为指标，并预 想事故的自动排序进行。 3）确定需要详细分析的。 4）对选出的预想事故用交流潮流进行校验。 74 n表征预想事故严重程度的两种行为指标(PI Performance Index): n(1)有功功率行为指标：是一种用来衡量线路有功功率过负荷程度的计 算方式，表示式为 n式中：p为有功功率权因子；Pl为线路l中的有功潮流；Plmax为线路l的 有功潮流限值；为有功功率过负荷的线路集合。 2 max l l pp P P PI 75 n(2)无功功率行为指标：是用来衡量电压与无功功率违限程度的计算公 式，表示式为 n式中：Ui为节点i的电压模值；Uilim为节点i的电压

43、模值限值；u为电压 权因子；Qi为节点i的无功注入；Qilim为节点i的无功注入限值；q为无 功功率权因子；为电压模值超过上、下限的节点集合； 为无功超 过上、下限的节点集合。 limlim limlim iiii vquq ii UUQQ PI UQ 76 n在上列两式中，、均只限于违限的线路或节点；p、u、q 的值则取决于系统的运行经验和在不同违限情况下有关线路的重要程 度；当权因子取为零时，即认为该线路违限并不重要而排除在集合之 外。 n 在研究线路行为指标时，可将所有线路，不论是否过负荷，都参加PI 计算；也可以只对支路潮流增加的线路参加PI计算。 77 n是指对一览表中的需要进行详细交流潮流计算的预想事故数进行选择 的判据 n常用的两种终止判据。 （1）只分析预想事故表