高等电力系统分析课件.ppt

《高等电力系统分析课件.ppt》由会员分享，可在线阅读，更多相关《高等电力系统分析课件.ppt（119页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、高等电力系统分析课程简介n 电网络分析基础电网络分析基础n 潮流算法及其扩展潮流算法及其扩展n 电力系统状态估计电力系统状态估计n 电力系统静态安全分析电力系统静态安全分析第一部分：系统分析篇1 电网络分析基础基础知识节点导纳矩阵电力网络方程求解方法第一部分：系统分析篇（续）2 潮流算法及其扩展潮流计算的数学模型潮流计算的经典算法保留非线性的潮流算法最小化潮流潮流计算中的自动调整最优潮流交直流潮流与含FACTS元件的系统潮流第一部分：系统分析篇（续）3 电力系统状态估计电力系统状态估计的基本概念最小二乘估计不良数据检测第一部分：系统分析篇（续）4 电力系统静态安全分析概述电力系统静态等值电力系

2、统预想事故选择1 电网络分析基础1.1 基础知识：电力网络的概念1 电力网络的概念电力网络的概念电力网络是指将输电配电线路、变压器等电气元件按一定形式连接而成的一个整体，达到输送和分配电能的目的。两个要素：电气元件及其连接方式电气元件及其连接方式。元件特性约束元件特性约束欧姆定律欧姆定律网络拓扑约束网络拓扑约束基尔霍夫定律基尔霍夫定律基尔霍夫电流定律（基尔霍夫电流定律（KCLKCL）基尔霍夫电压定律（基尔霍夫电压定律（KVLKVL）1.1 基础知识：电力网络的描述方法2 电力网络的描述方法电力网络的描述方法基于基尔霍夫电流定律：节点电压方程基于基尔霍夫电压定律：回路电流方程用节点电压方程描述电

3、力用节点电压方程描述电力网络的一个例子网络的一个例子按节点电压整理后得到：n左式中，左端是由各节点流出的电流，右端是向各节点注入的电流。n左式可以表示为规范的形式 以基尔霍夫电流定律列出节点方程：前述式子表示为规范形式如下：可以看出，其中的元素如下：n左式中，即为相应节点间的自导纳及互导纳。其余节点间互导纳为零。n上式为电力网络的节点方程节点方程。n在求出节点电压后，就可以求出各支路电流，从而使网络变量得以求解。n节点方程反映了各节点电压与注入电流间的关系。在此例中，除节点4、5外，其余节点注入电流均为0。重写规范形式如下：一般情况下，如果电力网络有n个节点，则有节点方程：式中：式中：nY是导

4、纳矩阵，对角元是节点i的自导纳，非对角元是节点间的互导纳。n分别是节点注入电流列向量及节点电压列向量1.1 基础知识：电力网络的关联矩阵描述例如，对上例所示的网络接线图，其节点-支路关联矩A为 对于N+1个节点b条支路的图，定义一个矩阵(行号对应节点号，列号对应支路号),矩阵中第i行第j列元素定义为 节点支路关联矩阵节点支路关联矩阵支路：支路：1 2 3 4 5 61 2 3 4 5 6节点节点节点节点节点大地作为参考节点节点关联矩阵网络拓扑结构对上图的节点，列KCL方程并写成矩阵形式为此方程组的系数矩阵就是该图的关联矩阵推广到一般情况：推广到一般情况：将将b b个支路电流写成支路电流向量，则

5、基尔个支路电流写成支路电流向量，则基尔霍夫电流定律的关联矩阵形式为霍夫电流定律的关联矩阵形式为 KCL KCL的关联矩阵形式的关联矩阵形式 KVL KVL的关联矩阵形式的关联矩阵形式此方程的系数矩阵等于图的关联矩阵 的转置选上图为例，用节点电压之差表示支路电压，并写成矩阵形式：推广到一般情况：推广到一般情况：设网络有设网络有b b条支路，条支路，n n个节点，第个节点，第n n号节点为号节点为参考节点，支路电压和节点电压向量分别记作参考节点，支路电压和节点电压向量分别记作：则节点电压与支路电压的关系即KVL：T21buuuL=UT1,21-=nnnnuuuLnUUUA=nT第k条广义支路的方程

6、可以表示成(k=1,b)b条支路的支路方程矩阵形式是(省略了复变量s)：简写为+-)(sUsk)(sIsk)(sIk)(sZk+-)(sUk)(),(sIsUkk)a()b()()()()()()()()()()(sUsIsZsIsZsUsIsIsZsUSkSkkkkSkSkkkk+-=+-=+-=SbSSSbSSbbbbUUUIIIZZZIIIZZZUUUMMLMOMMLLMLMOMMLLM212121212121000000000000SSUZIZIU+-=如何表示支路特性约束如何表示支路特性约束欧姆定律欧姆定律2023/1/6高等电力网络分析18其中U U、I I-支路电压向量与支路电流

7、向量-支路源电压与支路源电流量-支路阻抗矩阵与支路导纳矩阵若矩阵Z Z存在逆矩阵 ，令 并乘在 两端，得 SSIYUYUI+-=2023/1/6高等电力网络分析19令(称节点导纳矩阵)节点电压方程简化为AIAI0 0移项后得 节点电压方程SSIYUYUI+-=0 0)(SSIYUYUAAI+-=UUA=nT0 0=+-SSnTAIAYUUAYASSnAIAYUUAYA-=TSSSnAIAYUI-=TAYAY=nSnnnIUY=矩阵A反映了网络的拓扑约束，Y反映了网络的支路特性约束，所以节点导纳矩阵集中了网络两种约束的全部信息。边界条件 如何表示整个网络如何表示整个网络节点电压方程节点电压方程2

8、023/1/6高等电力网络分析20nsnnUIZ=若网络参数用阻抗形式表示，则节点网络方程有如下形式：2023/1/6高等电力网络分析21 关联矢量的引入关联矢量的引入一般串联支路2023/1/6高等电力网络分析22广义关联矢量和变压器/移相器支路的数学描述1.2 节点导纳矩阵：物理意义节点导纳阵反映了电力网络的参数及接线情况节点导纳阵节点电压方程 的推导过程SnnnIUY=TAYAY=n1.2 节点导纳矩阵：物理意义由导纳矩阵所构成的节点方程式是电力网络广泛应用的一种数学模型。SnnnIUY=n节点节点导纳物理意义：导纳物理意义：n如果在节点i加一单位电压，而把其余节点全部接地则上述节点方程

9、式成为n节点自导纳Yii 节点i加单位电压，其它节点接地时，节点i向电网注入的电流。n节点互导纳Yji节点i加单位电压，其它节点接地时，节点j向电网注入的电流。n特点：n当不含移相器时，导纳阵为对称矩阵n导纳矩阵为稀疏矩阵n出线数24条，每行非对角元中仅有24个非零元n例如，节点数分别10，1000的两个网络，平均出线为3n前者非零元40个，占总数40。n后者非零元4000个，占总数0.4。n计算时充分利用对称及稀疏性1.2 节点导纳矩阵：导纳矩阵的特点1.2 节点导纳矩阵：导纳矩阵的形成1.1.矩阵计算形成：节点矩阵计算形成：节点-支路关联矩阵支路关联矩阵节点导纳矩阵：TAYAY=n矩阵A为

10、节点-支路关联矩阵，Y为支路原始导纳阵。例：有以下三节点网络n导纳矩阵有如下形式，现考虑如何求其中各元素1.2 节点导纳矩阵：导纳矩阵的形成2.2.按支路逐条形成：关联矢量按支路逐条形成：关联矢量从图中可以看出：n形成导纳阵第一列元素 Y11，Y21，Y31。应在节点1加单位电压，节点2、3接地。从图中可以看出：n形成导纳阵第二列元素 Y12，Y22，Y32。应在节点2加单位电压，节点1、3接地。从图中可以看出：n形成导纳阵第三列元素 Y13，Y23，Y33。应在节点3加单位电压，节点1、2接地。n最后，得到该网络的导纳矩阵推广到一般情况：nk=i时，上式说明，当网络中除节点i以外所有节点都接

11、地时，从节点i注入网络的电流同施加于节点i的电压之比，即节点自导纳Yii。节点i加单位电压，其它节点接地时，节点i向电网注入的电流。n自导纳Yii是节点i以外的所有节点都接地时节点i对地的总导纳。显然，应等于与节点i相接的各支路导纳之和。得则令nki时，上式说明，当网络中除节点k以外所有节点都接地时，从节点i注入网络的电流同施加于节点k的电压之比，即节点互导纳Yik。节点k加单位电压，其它节点接地时，节点i向电网注入的电流。n此时节点i的电流实际上是自网络流出并进入地中的电流，所以互导纳Yik应等于节点i，k间的支路导纳的负值。Y以地为参考点的节点导纳矩阵ANxb阶节点支路关联矩阵MlA的第l

12、个列矢量按支路扫描，累加每条支路对导纳矩阵的贡献，最后就得到Y矩阵。对互感支路，应将互感支路组成一组，共同考虑它们对节点导纳矩阵的贡献。1.2 节点导纳矩阵：导纳矩阵的形成按支路逐条形成：关联矢量按支路逐条形成：关联矢量=i p j q 1、支路的移去和添加2、节点合并pqP导纳矩阵中相应的行列相加，网络方程降低一阶1.2 节点导纳矩阵：导纳矩阵的修正3、节点消去消去节点p，只需对Y阵中和p有支路相连的节点之间的元素进行修正，其他节点之间的元素不需要修正。4、节点电压给定的情况展开得：5、变压器变比发生变化的情况 自己思考（略）6、一条支路导纳参数发生变化的情况 自己思考（略）7、移去和添加带

13、互感支路的情况添加一条和原网络中支路k有互感的连支支路l时，可分两步进行修正：1）将支路k移出；2）将支路l和k成组追加进去。1.3 电力网络方程求解方法：高斯消去法常用方法有高斯消去法和因子表法常用方法有高斯消去法和因子表法1高斯消去法设有n阶线性方程组a11x1+a12x2+a1nxn=b1a21x1+a22x2+a2nxn=b2 .(1)an1x1+an2x2+annxn=bn 或缩记为：AX=B (2)按列消元按行回代的算法按列消元按行回代的算法按列消元按行回代的算法按列消元按行回代的算法增广A阵求解的具体步骤如下：求解的具体步骤如下：第一步：按列消去。消去第1列第1行规格化：得消去第

14、1列下三角元素：则 变成一般地，消去第k列：第k行规格化：消去第k列下三角元素：则 变成最后可得：写成方程组形式：它与原方程 同解第二步：按行回代第n行 将结果代入第n-1行，得一般地，将 代入第i个方程，得例：按列消元按行回代的高斯消去法例：按列消元按行回代的高斯消去法由原方程写出增广矩阵第1列规格化第1列消去第2列规格化第2列消去第3列规格化第3列消去第4列规格化原方程改写成：回代：设有n阶线性方程组a11x1+a12x2+a1nxn=b1a21x1+a22x2+a2nxn=b2 .(1)an1x1+an2x2+annxn=bn或缩记为：AX=B (2)在实际计算中，经常遇到这种情况：对于

15、方程组需要多次求解，每次仅改变其常数项B，而系数矩阵A是不变的。这时，为了提高计算速度，可以利用因子表求解。1.3 电力网络方程求解方法：因子表法因子表法的基本概念因子表法的基本概念因子表可以理解为高斯消去法解线性方程组的过程中对常数项B全部运算的一种记录表格。高斯消去法分为消去过程和回代过程。回代过程的运算由对系数矩阵进行消去运算后得到的上三角矩阵元素确定，公式：为了对常数项进行消去运算（又叫前代过程），还必须记录消去过程运算所需要的运算因子。消去过程中的运算又分为规格化运算和消去运算，以按列消去过程为例，公式：规格化：（i=1,2,n）将上式中的运算因子 及 逐行放在下三角部分，和消去过程

16、得到的上三角矩阵元素合在一起，就得到了因子表。消去：(k=1,2,i-1)（ij）(ji)因子表中下三角部分的元素就是系数矩阵在消去过程中曾用以进行运算的元素，因此只要把它们保留在原来的位置，并把对角元素取倒数就可以得到因子表的下三角部分。而因子表中上三角部分的元素就是系数矩阵在消去过程完成后的结果。下三角上三角阵对角阵记 即 则因子分解迭代格式：n用因子表法求解线性方程组对于方程组，需要多次求解，每次仅改变其常数项B而系数矩阵A是不变的情况，应首先对其系数矩阵A进行消去运算，形成因子表。有了因子表，就可以对不同的常数项B求解。这时，可以直接应用因子表中的元素。消去（i=k+1,n）回代例：用

17、因子表法求解下述方程组例：用因子表法求解下述方程组解：对照前例，形成因子表解方程消去第1列：规格化：消去：得消去第2列：规格化：消去：得得消去第3列：消去第4列：得：原方程变为：逐行回代，得：1、由于电力网络结构的特点，每个节点仅与35个节点相连，因此描述网络结构的矩阵是稀疏矩阵。n*m的矩阵，非零元个，稀疏度等于/n*m 如果系统有N=500个节点，平均每个节点与5条支路相连，则 稀疏度=5*500/（500*500）=1%2、计算中，我们仅关心一部分的变量：稀疏矢量。3、与稀疏矩阵和稀疏矢量相关的运算中，零元素不参与存储和计算排零存储排零存储和排零计算排零计算1.4 电力网络求解的稀疏技术

18、：引入原因特点特点：排零存储，即只存储其中的非零元和有关的检索信息。要求：节省内存方便地检索和存取考虑网络结构变化时能方便地对存储的信息加以修改稀疏矢量：存储矢量中的非零元值和相应的下标稀疏矩阵：考虑稀疏结构和所采用的算法1.4 电力网络求解的稀疏技术：稀疏存储VAa11a12a14a21a22a23a33a42a43a44IA1112223444JA1241233234散居格式a1313 常用存储方式常用存储方式散居格式按行（列）存储格式三角检索存储格式链表存储格式例：VAa11a12a14a21a22a23a33a42a43a44JA1241233234IA1478按行存储格式VAa11a

19、12a13a14a21a22a23a33a42a43a44JA12341233234IA1589修改后三角检索存储格式Ua12a14a23JU243IU1344La21a42a43IL244JL1234Da11a22a33a44Ua12a13a14a23JU2343IU1455La21a42a43IL244JL1234Da11a22a33a44修改后链表存储格式VAa11a12a14a21a22a23a33a42a43a44JA1241233234LINK23056009100IA1478VAa11a12a14a21a22a23a33a42a43a44a13JA12412332343LINK2

20、110560091003IA1478修改后小结小结存储单元检索修改散居格式3工作量大易按行存储格式 2+n按行检索难三角检索存储格式2+n按行（列）检索难链表格式3+n按行检索易1、稀疏矩阵的因子分解、稀疏矩阵的因子分解采用高斯消去法进行计算1）按行规格化2）消去运算1.4 电力网络求解的稀疏技术：排零计算采用三角检索存储格式时例U7-3JU24IU1333L5-2IL24JL1233D2456U3.5-1.57.5JU244IU1344L5-2IL24JL1233D2-13.556U3.5-1.5-0.555JU244IU1344L5-2IL24JL1233D2-13.554.889U3.5

21、-1.5-0.555JU244IU1344L5-2IL24JL1233D2-13.554.8893=p2 2、利用稀疏矩阵因子表求解稀疏线性代数方程组、利用稀疏矩阵因子表求解稀疏线性代数方程组1）前代过程计算流程除法运算回代运算计算流程采用三角检索存储格式时采用按列存储格式Uu12u24u34u15u25u45IU123124JU12247例U3.5-1.5-0.555IU244IU1344L5-2IL24JL1233D2-13.554.889A图有向A图赋权有向A图A图：和矩阵A有相同拓扑结构的网络图有向A图：对给定A图及节点编号，规定边的正方向由小号节点指向大号节点赋权有向A图：在有向A图

22、中，将A的非对角非零元的值赋给互边，将A的对角元素的值赋给自边1 1、基本定义和术语、基本定义和术语1.4 电力网络求解的稀疏技术：基于图论因子图有向因子图赋权有向因子图因子图：和因子表矩阵U有相同拓扑结构的网络图有向因子图：在因子图上规定边的正方向由小号节点指向大号节点形成赋权有向因子图：在有向因子图中，将U的非对角非零元的值赋给互边，将对角线矩阵D的元素的值赋给自边规格化在赋权有向A图上，相当于对于节点p发出的所有互边的边权加以修正，新的边权等于原边权除以节点p的自边边权。消去运算对角元修正：在赋权有向A图上，就是对节点p发出的边的收点上的自边边权进行修正。非对角元修正：在赋权有向A图上，

23、就是对节点p发出的边的中任取两边，其收点所夹的边的边权应减少的数量是p点发出的两条边的边权与p点自边边权的乘积。2 2、因子分解过程的图论描述、因子分解过程的图论描述算法流程在赋权有向A图上按节点号由小到大的顺序（例如对节点p）执行下面的操作：（1）对节点p发出的互边将其边权除以节点p的自边边权；（2）对节点p发出的互边的收点，将该点上的自边边权减去该互边边权平方乘以节点p的自边边权；（3）对节点p发出的所有互边，这些互边两两之间所夹得互边边权应减去两条相夹边边权与节点p的自边边权三者乘积。操作前被节点对之间无边的情况应视为有一条零权值边。例2-1-1-12-1-12-1-1-14节点1规格化

24、消去节点2规格化消去节点3规格化消去3 3、前代回代过程的图论描述、前代回代过程的图论描述计算流程（1）将独立矢量b的非零元赋值为赋权有向因子图上的点位e；（2）扫描i从1到n-1，从公式修正节点i发出的边的收端节点j的点位（3）对所有节点，用公式对点位规格化（4）扫描j从n到2，对所有指向节点j的边的发端节点i，用公式修正其点位例：在下面的赋权有向因子图上进行前代和回代。已知独立矢量为：赋权有向因子图和独立矢量点位前代点位点位1：为零不用计算点位点位2：点位点位3：规格化回代节点4节点3节点2结果结果已知：对称矩阵A的互边的边数是b，自边的边数是n，有向因子图上的互边边数是 。消去运算的乘法

25、次数：节点p发出的边有 条，则乘法次数总乘法次数为 规格化中的乘法次数：总乘法次数1.4 电力网络求解的稀疏技术：计算代价分析相关概念相关概念稀疏独立矢量，一个给定的只有少量非零元的独立矢量。稀疏解矢量，一个只有少数元素待求的解矢量，其余元素我们不关心道路树，在有向因子图上，从每个节点发出的边中取收点号最小的边作为树边，这样得到的道路树。点的路，在道路树上该点沿道路树到树根所经过的路径，它是道路树的一个子集。点集的路集，是该点集中所有点的路的并集。1.4 电力网络求解的稀疏技术：稀疏矢量技术一个例子一个例子(a)有向因子图(b)道路树(c)点1的路(d)点集1，4，8的路集定理1：在有向因子图

26、上，前代运算只在稀疏矢量中非零元点集的路集上进行定理2：路集上任一点的前代运算必须在路集上比该点编号小且其道路经过该点的点的前代完成之后才能进行，而路集中分支点以下的几点路先做哪个没有关系。定理3：在有向因子图上，回代运算只在稀疏解矢量中待解元素的点集的路集上进行定理4：路集上任一点的回代运算必须在路集上比该点编号大且其道路经过该点的点的回代完成之后才能进行，而路集中分支点以上的几点路先做哪个没有关系。道路集的形成道路集的形成例1，按行存储，上三角矩阵中该行第一个非零元的列号12345678910Uu17u1,12u25u29u36u39u47u4,10u59u5,12JU7125969710

27、912IU1357911131517191112131415161718192021Uu69u6,10u7,10u7,12u89u8,11u9,10u9,12u10,11u10,12u11,12JU91010129111012111212IU2122 寻找节点p的道路 寻找点集的路集点集G中点的路集P 计算代价分析计算代价分析不用稀疏矢量法稀疏矢量法前代运算F(PF个节点发出的边数）除法nPF(b中非零元点集的路集中的节点数)回代运算B(x中非零元点集的路集中的节点发出的边数)总计2 +nF+PF+B1、注入元素的多少与消去节点的顺序或节点编号有关1.4 电力网络求解的稀疏技术：节点优化编号所

28、谓节点优化编号，就是寻找一种使注入元素数目最少的节点编号方式。2、三类节点编号优化方法 静态按最少出线支路数编号静态优化法编号之前，首先统计电力网络各节点的出线支路数，然后按出线支路数少的节点顺序编号，当有n个节点的出线支路数相同时，则可以按任意次序对这n个节点编号。依据：在导纳矩阵中，出线支路数最少的节点所对应得行中非零元素也最少，因此在消去过程中产生注入元素的可能性也最小。缺点：未考虑节点消去过程中，每消去一个节点，与该节点相连的各节点的出线支路数将发生变化。动态地按最少出线支路数编号半动态优化法针对静态优化法的缺点，在每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中出线支路

29、数最少的一个节点进行编号。缺点：只能使消去过程中出现新支路的可能性减少，但并不一定保证在消去这些节点时出现的新支路最少。动态按增加出线数最少编号动态优化法针对上述缺点，采用按消去节点后增加出线数最少的原则编号。具体做法：首先，根据星网变换原理，按下式分别统计消去网络节点时增加的出线数，选其中增加出线数最少的被消节点编为第1节点。如果与节点k相连的节点数为 ，则网形网络的支路数为 ，原有支路数为 ，则新增支路数从网络消去该节点，相应修改其余节点的出线数目。然后重复以上过程，一直到编完为止。缺点：工作量大。例：静态优化法统计各节点出线支路数节点ABCDEFGH出线数11322231F点总是编在C和

30、G之前，故CG两点之间出现新支路难以避免。半动态优化法节点ABCDEFGH被编节点节点号各节点出线数的变化情况(1)111(1)3332222(1)22(1)2(1)222222(1)33332(1)1111(1)AEDBHGFC12345678没有出现新支路。结论结论：对于树形网络来说，半动态优化法永远只编出线为1的节点，因此这种任意性不会影响优化结果。动态优化法被消节点ABCDEFGH出现新出线数00311130A:1被消节点BCDEFGH出现新出线数0310130重复。工作量比半动态优化法大得多，但对于树形网络，效果和半动态优化一样。E:2推广：推广：网络类型静态半动态动态树形可能出现新

31、支路无新支路无新支路简单环网（m边形）m-3新支路复杂环网无明确的一般性结论小小 结结1、稀疏技术包括稀疏矩阵技术和稀疏矢量技术，是电网计算中使用最为广泛的计算技术。2、稀疏技术的关键在于排零存储和排零计算。稀疏矩阵技术充分开发网络矩阵的稀疏结构，减少和稀疏矩阵有关的计算量。稀疏矢量技术充分开发矢量的稀疏性，在前代回代计算中只进行和稀疏矢量中非零元有关的计算，省略了不必要的计算，以进一步提高求解网络方程的计算速度。3、稀疏矩阵技术和稀疏矢量技术可用图的方法来描述。赋权有向A图包含了矩阵A的所有信息，赋权有向因子图包含了矩阵A的因子表矩阵的所有信息。图上因子分解形象说明了稀疏矩阵技术中排零存储和

32、排零计算的实质。4、节点优化编号对稀疏技术性能的提高至关重要。半动态节点优化编号简单有效，得到最为广泛的应用，可大大减少因子分解过程中注入元的数量。在对大规模互联电力系统进行统一分析时，分块计算是一种提高计算速度的有效处理手段。电力系统本身所具有的分层分区结构也特别适合分块计算的应用。根据协调变量的不同，网络分块计算主要分为两类：支路切割法支路切割法：通过切割原网络中的某些支路把原网络分解；节点撕裂法节点撕裂法：将原网络的部分节点撕裂开，将网络分解1.5 大型电力网络的分块计算1 1、节点分裂法、节点分裂法在该网络中选择部分节点，把这些节点撕裂，则把原网络可以分解成几个小的独立子网络，这些节点

33、称为分裂点，用下标t表示。1.5 大型电力网络的分块计算：节点分裂法若分裂点电压 已知，则每个子网络的节点电压可以用下式计算：求分裂点电压 消去个子网络所对应的网络，只保留分裂点t相对应的部分，有分裂节点的电压带有各子系统相互之间的协调信息，也称协调变量。节点分裂法的物理解释 计算分裂点的等值导纳矩阵等值导纳亦可如下表示 计算节点等值注入电流综上所述，等值后网络如下图所示计算分裂节点电压分别计算两个子网络的节点电压1 1、常规支路切割法、常规支路切割法在给定的电力网络中选择部分网络，将这些支路切开，如果此时网络能变成几个相互独立的子网络，就把这些支路称为切割支路。1.5 大型电力网络的分块计算

34、：支路切割法 是关键变量，即协调变量消去法求解2 2、广义支路切割法、广义支路切割法问题问题：原网络有一个接地点，现分成k个子网后，有k-1个子网没有接地点，节点导纳矩阵奇异，怎么办？令消去电压 ，得计算切割支路电流这样，并不是物理电网支路L上的电流，可理解为一个广义电流广义支路切割法。3 3、常规支路切割法的物理解释、常规支路切割法的物理解释简记为原网络子网络对既包含撕裂节点，又包含切割支路的网络，写成：或简写成：1.5 大型电力网络的分块计算：统一算法协调变量求解：求解子系统电压网络分块解法的并行计算特性分析1）由于网络方程的系数矩阵是加边块对角结构，所以，不论是求子网络的等值电路（节点分

35、裂法的诺顿等值和支路切割法的戴维南等值）还是求解子网络内部的电量，都可以分别独立地进行，子网络相互之间不受影响；2）在求解网络的诺顿等值或戴维南等值时，并不需要求子网络内部变量，而是求边界协调变量；3）求解协调变量时方程的阶次很低，计算速度很快；4）网络方程的分块解法用相对少的计算代价计算出协调变量，而后用并行计算机计算每个子网络。1.5 大型电力网络的分块计算：实际应用大规模电网的分解协调计算和并行计算大规模电网的分解协调计算和并行计算子系统协调级利用稀疏矩阵技术提高各子系统的计算速度：利用稀疏矩阵技术提高各子系统的计算速度：节点导纳矩阵的因子表各项展开，得所以大规模电网分块计算的实际应用大规模电网分块计算的实际应用应用算法1）多计算机集群的并行计算 多台并行计算机处于同一地点，可以采用高速通讯技术，用同步迭代来实现。2）多控制中心的分解协调计算 各个控制中心计算机处于异地，通信瓶颈问题突出，通常只能用异步迭代实现。应用领域1）短路电流计算采用支路切割法进行分块并行计算，不需迭代。2）暂态稳定计算联立求解DAE3）其他方面多区域经济分配计算，多区域发电计划安排，联络线潮流控制，区域经济交换功率控制等方面