第四章 复杂电力系统潮流的计算机算法 《电力系统稳态分析》.

Presentation on theme: "第四章 复杂电力系统潮流的计算机算法 《电力系统稳态分析》."— Presentation transcript:

1 第四章 复杂电力系统潮流的计算机算法 《电力系统稳态分析》

2 第四章 复杂电力系统潮流的计算机算法 基本要求：本章着重介绍运用电子计算机计算电力系统潮流分布的方法。它是复杂电力系统稳态和暂态运行的基础。
第四章 复杂电力系统潮流的计算机算法 基本要求：本章着重介绍运用电子计算机计算电力系统潮流分布的方法。它是复杂电力系统稳态和暂态运行的基础。 运用计算机计算的步骤，一般包括建立数学模型，确定解算方法，制定框图和编制程序，本章着重前两步。

3 本章知识点： 1、节点导纳矩阵，节点导纳矩阵各元素的物理意义，如何由节点导纳矩阵形成节点阻抗矩阵，节点阻抗矩阵各元素的物理意义，导纳矩阵与阻抗矩阵的对称性和稀疏性； 2、网络节点分类，数学模型中已知条件和待求量； 3、牛顿－拉夫逊迭代法原理，牛顿－拉夫逊迭代法直角坐标形式的功率误差方程和电压误差方程，牛顿－拉夫逊迭代法极坐标形式的雅可比矩阵与修正方程，两种修正方程的不同点，牛顿－拉夫逊迭代法两种坐标系潮流计算求解步骤；

4 4、高斯－赛德尔法潮流原理，非线性节点电压方程的高斯－赛德尔迭代形式，PV节点向PQ节点转化的原因和方法；
5、P－Q分解法潮流计算， P－Q分解法与牛顿－拉夫逊的关系，由牛顿－拉夫逊法导出P－Q分解法用到了几个近似条件，各近似条件的物理意义， P－Q分解法的修正方程式， P－Q分解法与牛顿－拉夫逊的迭代次数与解题速度， P－Q分解法分解法潮流计算求解步骤。

5 4－1 电力网络方程 电力网络方程指将网络的有关参数和变量及其相互关系归纳起来组成的，反映网络特性的数学方程式组。如节点电压方程、回路电流方程，割集电压方程。相应有： （1）节点导纳矩阵 （2）节点阻抗矩阵 （3）回路阻抗矩阵

6

7 一、节点电压方程 ～ 电力网 网络元件：恒定参数 发电机：电压源或电流源 负荷：恒定阻抗 代数方程

8 一、节点电压方程 ～ 以母线电压作为待求量 注意： 零电位是不编号的 1 2 3 4 负荷用阻抗表示 1 2 3 4 E1 E4
电力系统结线图 注意： 零电位是不编号的 负荷用阻抗表示 1 2 3 4 E1 E4 电力系统等值网络

9 一、节点电压方程 1、节点导纳方程 电压源变为电流源 1 2 4 I1 y24 3 I4 y10 y12 y20 y23 y34 y40
以零电位作为参考，根据基尔霍夫电流定律

10 一、节点电压方程 1、节点导纳方程

11 一、节点电压方程 1、节点导纳方程 其中

12 一、节点电压方程 1、节点导纳方程 n 个独立节点的网络，n 个节点方程

13 一、节点电压方程 1、节点导纳方程 n 个独立节点的网络，n 个节点方程

14 一、节点电压方程 1、节点导纳方程 Y 节点导纳矩阵 Yii 节点i的自导纳 Yij 节点i、j间的互导纳
n 个独立节点的网络，n 个节点方程 Y 节点导纳矩阵 Yii 节点i的自导纳 Yij 节点i、j间的互导纳

15 一、节点电压方程 1、节点导纳方程 Y 矩阵元素的物理意义

16 一、节点电压方程 1、节点导纳方程 Y 矩阵元素的物理意义 自导纳
Ykk：当网络中除节点k以外所有节点都接地时，从节点k注入网络的电流同施加于节点k的电压之比 Ykk：节点k以外的所有节点都接地时节点k对地的总导纳

17 一、节点电压方程 1、节点导纳方程 Y 矩阵元素的物理意义 互导纳
Yki：当网络中除节点k以外所有节点都接地时，从节点i注入网络的电流同施加于节点k的电压之比 节点i的电流实际上是自网络流出并进入地中的电流，所以Yki应等于节点k、i之间导纳的负值

18 一、节点电压方程 1、节点导纳矩阵Y 节点导纳矩阵中自导纳和互导纳的确定 3 I4 y10 y12 y20 y23 y34 y40 y30

19 一、节点电压方程 1、节点导纳矩阵Y 节点导纳矩阵中自导纳和互导纳的确定 1 2 3 4 y10 y12 y20 y23 y24 y34
－ y10 y12 y20 y23 y24 y34 y40 ＋ y30

20 一、节点电压方程 1、节点导纳矩阵Y 节点导纳矩阵中自导纳和互导纳的确定 1 2 3 4 y10 y12 y20 y23 y24 y34
－ y10 y12 y20 y23 y24 y34 y40 ＋ y30

21 一、节点电压方程 1、节点导纳矩阵Y 节点导纳矩阵中自导纳和互导纳的确定 1 2 3 4 y10 y12 y20 y23 y24 y34
－ y10 y12 y20 y23 y24 y34 y40 ＋ y30

22 一、节点电压方程 1、节点导纳矩阵Y 节点导纳矩阵中自导纳和互导纳的确定 1 2 3 4 y10 y12 y20 y23 y24 y34
－ y10 y12 y20 y23 y24 y34 y40 ＋ y30

23 一、节点电压方程 节点导纳矩阵Y 的特点 直观易得 稀疏矩阵 对称矩阵

24 一、节点电压方程 2、节点阻抗矩阵 Z 矩阵元素的物理意义

25 一、节点电压方程 2、节点阻抗矩阵 Z = Y -1 节点阻抗矩阵 Zii 节点i的自阻抗或输入阻抗 Yij 节点i、j间的互阻抗或转移阻抗

26 一、节点电压方程 2、节点阻抗矩阵 Z 矩阵元素的物理意义

27 一、节点电压方程 2、节点阻抗矩阵 Z 矩阵元素的物理意义
在节点 k 单独注入电流，所有其它节点的注入电流都等于 0 时，在节点 k 产生的电压同注入电流之比 从节点 k 向整个网络看进去的对地总阻抗

28 一、节点电压方程 2、节点阻抗矩阵 Z 矩阵元素的物理意义互阻抗
在节点 k 单独注入电流，所有其它节点的注入电流都等于 0 时，在节点 i 产生的电压同注入电流之比

29 一、节点电压方程 2、节点阻抗矩阵 节点阻抗矩阵中自阻抗和互阻抗的确定 1 2 3 4 I1 I4 z10 z12 z20 z23 z24

30 一、节点电压方程 2、节点阻抗矩阵 节点阻抗矩阵中自阻抗和互阻抗的确定 1 2 3 4 z10 z12 z20 z23 z24 z34

31 一、节点电压方程 2、节点阻抗矩阵 节点阻抗矩阵中自阻抗和互阻抗的确定 1 2 3 4 z10 z12 z20 z23 z24 z34

32 一、节点电压方程 2、节点阻抗矩阵 节点阻抗矩阵中自阻抗和互阻抗的确定 1 2 3 4 z10 z12 z20 z23 z24 z34

33 一、节点电压方程 2、节点阻抗矩阵 节点阻抗矩阵中自阻抗和互阻抗的确定 1 2 3 4 z10 z12 z20 z23 z24 z34

34 一、节点电压方程 2、节点阻抗矩阵 Z 矩阵的特点 复杂难求（Y－1，支路追加法） 满矩阵

35 二、回路电流方程 回路阻抗矩阵 1 2 3 4 z'10 z12 z20 z23 z24 z34 z'40 z30 + -

36 二、回路电流方程 回路阻抗矩阵

37 二、回路电流方程 回路阻抗矩阵 m 个独立回路的网络，m个节点方程

38 二、回路电流方程 回路阻抗矩阵 m 个独立回路的网络，m 个节点方程

39 二、回路电流方程 回路阻抗矩阵 ZL 回路阻抗矩阵 IL 回路电流列相量；（习惯取顺时针的电流流向为正）
m 个独立回路的网络，m 个节点方程 ZL 回路阻抗矩阵 IL 回路电流列相量；（习惯取顺时针的电流流向为正） EL 回路电压源电势的列相量，与IL方向一致为正。

40 二、回路电流方程 回路阻抗矩阵 ZL 矩阵元素的物理意义 Zii：自阻抗，环绕回路i所有支路阻抗的总和；
Zij：互阻抗，回路i和回路j共有的阻抗，其中Zij＝Zji，如回路j、i无共有阻抗，则Zij＝Zji＝0

41 二、回路电流方程 回路阻抗矩阵 ZL 矩阵的特点 对称矩阵 稀疏矩阵

42 三、节点导纳矩阵 Y 矩阵的修改 不同的运行状态，（如不同结线方式下的运行状况、变压器的投切或变比的调整等）
改变一个支路的参数或它的投切只影响该支路两端节点的自导纳和它们之间的互导纳，因此仅需对原有的矩阵作某些修改。

43 三、节点导纳矩阵 Y 矩阵的修改 不同的运行状态，（如不同结线方式下的运行状况、变压器的投切或变比的调整等） 电力网

44 三、节点导纳矩阵 Y 矩阵的修改 电力网

45 三、节点导纳矩阵 Y 矩阵的修改 （1）从原网络引出一条支路增加一个节点 Y 增加一行一列（n＋1）×（n＋1） 电力网 yik i k

46 三、节点导纳矩阵 Y 矩阵的修改 （2）在原有网络节点i、j之间增加一条支路 Y 阶次不变 yij i j 电力网

47 三、节点导纳矩阵 Y 矩阵的修改 （3）在原有网络的节点i、j之间切除一条支路 Y 阶次不变 yij 电力网 i j

48 三、节点导纳矩阵 Y 矩阵的修改 （4）在原有网络的节点i、j之间的导纳由yij改变为y'ij i j -yij 电力网 y'ij

49 三、节点导纳矩阵 Y 矩阵的修改 （5）在原有网络的节点i、j之间变压器的变比由k*改变为k*' ZⅠ ZⅡ i j k*:1 ZT ZⅠ
yT/k*

50 三、节点导纳矩阵 Y 矩阵的修改 （5）在原有网络的节点i、j之间变压器的变比由k*改变为k*'

51 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程 等值电源功率 G 1 2 等值负荷功率 （a）简单系统

52 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程 G 1 2 y10 y20 y12 （b）简单系统的等值网络

53 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程 1 2 y10 y20 y12 — （c）注入功率和注入电流

54 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程 1 2 y10 y20 y12 — （c）注入功率和注入电流

55 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程

56 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程

57 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程

58 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 1、功率方程 决定功率大小的是相对相位角或相对功率角 有功、无功功率损耗为：

59 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 2、变量的分类 除网络参数外，共有十二个变量
（1）负荷消耗的有功、无功功率－PL1、PL2、QL1、QL2。取决于用户，不可控变量或扰动变量，用列向量d表示。 （2）电源发出的有功、无功功率－PG1、PG2、QG1、QG2。控制变量，用列向量μ表示。 （3）母线或节点电压的大小和相位角－U1、U2、δ1、δ2。状态变量或受控变量，U→Q， δ →P，用列向量x表示。

60 ? 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 2、变量的分类
对于n个节点，变量数增为6n，其中d、μ、x各2n个。 将上述变量进行分类后，只要已知或给定扰动变量和控制变量，就可运用功率方程式解出状态变量U，δ。 ? 但是当δ1 、δ2 变化同样大小时，功率的数值不变，从而不可能求出绝对相位角，相应的功率损耗也不能确定。

61 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 2、变量的分类
为克服上述困难，在一个具有n个节点的系统中，对变量的给定稍作调整： （1）只给定（n-1)对控制变量PGi、QGi，余下一对控制变量PGs、QGs待定，以使系统功率保持平衡； （2）给定一对δs、Us，其中； PLi、QLi均为已知。 求解（n-1)对状态变量及一对待定的控制变量

62 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 控制变量 2、变量的分类 得出的解应满足如下约束条件：
取决于一系列的技术经济因素

63 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 2、变量的分类 得出的解应满足如下约束条件： 节点状态变量 扰动变量

64 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 3、节点的分类

65 4－2 功率方程及其迭代解法 一、功率方程和变量、节点的分类 3、节点的分类
(1) PQ节点：PLi、QLi；PGi、QGi，即相应的Pi、Qi给定，待求Ui、δi。如按给定有功、无功发电的发电厂母线和没有其他电源的变电所母线 (2) PU节点： PLi、 PGi ，从而Pi给定； QLi 、Ui给定。即相应的Pi、Ui给定，待求QGi、δi。如有一定无功储备电源变电所母线（很少，甚至没有）。 (3) 平衡节点： 一般只有一个。设s节点为平衡节点，则： PLs、QLs ；Us 、 δs 给定， Us ＝1.0， δs ＝0。待求PGs、QGs。

66 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程）

67 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程）

68 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程）

69 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程） 若式中的aij对于Yij、xi对应Ui，yi对应

70 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程）
此时可用迭代法求解。如设节点1为平衡节点，其余为PQ节点，则有：

71 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程）
此时可用迭代法求解。如设节点1为平衡节点，其余为PQ节点，则有：

72 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程）
此时可用迭代法求解。如设节点1为平衡节点，其余为PQ节点，则有： 计算步骤为：

73 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程） 对各类节点的计算和处理
由于节点的类型不同，已知条件和求解对象不同，约束条件不同，在计算过程中的处理不同。 （1）PQ节点：按标准迭代式直接迭代； （2）PV节点：已知的式Pp和Up，求解的是Qp，δp；按标准迭代式算出Up (k)， δp (k)后，首先修正: 然后修正

74 4－2 功率方程及其迭代解法 二、高斯－赛德尔迭代法（既可解线性，也可解非线性方程） 对各类节点的计算和处理
检查无功是否越限，如越限，取限值,此时：PV→PQ

75 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程） 原理： 按泰勒级数展开，并略去高次项

76 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程） 原理：

77 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程）

78 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程）

79 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程）

80 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程）

81 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程）

82 例题：如图所示，母线1为平衡节点，δ1＝0，U1＝1. 0，母线2为PV节点，U2＝0
例题：如图所示，母线1为平衡节点，δ1＝0，U1＝1.0，母线2为PV节点，U2＝0.95，P2＝PG2－PL2＝4－2＝2，母线3为PQ节点， P3＝－PL3＝－4.0 ， Q3＝－QL3＝－1.5 。试写出此系统的功率方程。 ～

83 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程）
（1）将xi(0)代入，算出△f，J中各元素，代入上式方程组，解出△xi(0)； （2）修正xi(1)＝ xi(0)＋ △xi(0) ，算出△f，J中各元素，代入上式方程组，解出△ xi(1) ；

84 4－2 功率方程及其迭代解法 三、牛顿－拉夫逊迭代法（常用于解非线性方程） 计算步骤：
（1）将xi(0)代入，算出△f，J中各元素，代入上式方程组，解出△xi(0)； （2）修正xi(1)＝ xi(0)＋ △xi(0) ，算出△f，J中各元素，代入上式方程组，解出△ xi(1) ； 注意：xi的初值要选得接近其精确值，否则将不迭代。

85 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 节点电压用直角坐标表示：

86 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 首先对网络中各节点作如下约定：
（1）网络中共有n个节点，编号为1，2，3，…，n； （2）网络中（m－1）个PQ节点，一个平衡节点，编号为1，2，…，m，其中1≤s≤m为平衡节点； （3）n－m个PV节点，编号为m+1,m+2,…，n.

87 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 (m-1)个PQ节点＋(n-m)个PV节点，共n-1个 (m-1)个PQ节点

88 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 相应的：

89 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 用直角坐标表示的修正方程 2(m－1) PQ节点 PV节点 2(n－m)

90 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 用直角坐标表示的修正方程

91 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 用直角坐标表示的修正方程

92 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 以极坐标表示的另一种修正方程式为 2(m－1) PQ节点 PV节点
2(n－m) PV节点 2(m－1) 2(n－m)

93 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 用极坐标表示的修正方程式为

94 4-3牛顿－拉夫逊迭代法潮流计算 一、潮流计算时的修正方程式 雅可比矩阵的特点：
（1）雅可比矩阵各元素均是节点电压相量的函数，在迭代过程中，各元素的值将随着节点电压相量的变化而变化。因此，在迭代过程中要不断重新计算雅可比矩阵各元素的值； （2）雅可比矩阵各非对角元素均与Yij＝Gij＋jBij有关，当Yij＝0，这些非对角元素也为0，将雅可比矩阵进行分块，每块矩阵元素均为2×2阶子阵，分块矩阵与节点导纳矩阵有相同的稀疏性结构； （3）非对称矩阵。

95 4-3牛顿－拉夫逊迭代法潮流计算 二、潮流计算基本步骤

96 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式 (n-1)×(m-1) (n-1)×(n-1) (m-1)×(m-1)

97 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式 1、对修正方程式的第一步简化
高压网络中，各元件的X>>R，δ→P，相应的J≈0；U →Q，N ≈0。

98 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式 2、对修正方程式的第二步简化
高压网络中，各元件的X>>R，使Gij<<Bij，再加上系统稳定性的要求，即| δi－ δj|< | δi－ δj|max， | δi－ δj|max＝（10 ° ～20°）。 3、对修正方程式的第三步简化

99 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式

100 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式

101 4-4 P－Q分解法潮流计算 一、潮流计算时的修正方程式 缩写为

102 4－2 网络方程的解法 高斯消去法 带有节点电流移置的星网变换

103 4－3 节点阻抗矩阵 用支路追加法形成 Z 矩阵 1 2 4 3

104 4－3 节点阻抗矩阵 用支路追加法形成 Z 矩阵 第一条支路必须是接地支路 后追加的支路必须至少有一个端点与已出现的节点相接

105 4－3 节点阻抗矩阵 用支路追加法形成 Z 矩阵 k 电力网 电力网 ziq q i zkm m 追加树枝 追加连枝

106 4－3 节点阻抗矩阵 追加树枝 电力网 ziq q i 追加树枝

107 4－3 节点阻抗矩阵 追加树枝 原有阻抗矩阵不变 m 电力网 ziq q i 从节点 m 单独注入电流 Im 追加树枝

108 4－3 节点阻抗矩阵 追加树枝 m 电力网 从节点 q 单独注入电流 Iq ziq q i 追加树枝

109 4－3 节点阻抗矩阵 追加树枝 m 电力网 ziq q i 追加树枝

110 4－3 节点阻抗矩阵 追加连枝 Ik Ik- Ikm 电力网 电力网 k k Ikm zkm m m Im Im+ Ikm Z 矩阵阶数不变

111 4－3 节点阻抗矩阵 追加连枝

112 4－3 节点阻抗矩阵 追加连枝

113 4－3 节点阻抗矩阵 追加连枝

114 4－3 节点阻抗矩阵 追加连枝 如果m点接地

115 4－3 节点阻抗矩阵 追加树枝 m 电力网 ziq q i 追加树枝

116 4－3 节点阻抗矩阵 追加连枝 如果将k、m两点短接，经过修改后，第k行（列）和第m行（列）的对应元素完全相同。只要将原来这两个节点的注入电流并到其中的一个节点，另外一个节点即可取消，并消去阻抗矩阵中对应的行和列，使矩阵降低一阶。

117 4－3 节点阻抗矩阵 追加变压器树枝 电力网 1:K zkq q k 原有阻抗矩阵不变

118 4－3 节点阻抗矩阵 追加变压器树枝 1:K ziq q m k 从节点 m 单独注入电流 Im

119 4－3 节点阻抗矩阵 追加变压器树枝 1:K ziq q m i 从节点 q 单独注入电流 Iq

120 4－3 节点阻抗矩阵 追加变压器连枝 1:K ziq q k 电力网 m

121 4－4 节点编号顺序的优化 消去时增加新支路最少的节点应该优先编号 简化为按节点的连接支路数k（接地支路除外）最少进行编号
4－4 节点编号顺序的优化 消去时增加新支路最少的节点应该优先编号 简化为按节点的连接支路数k（接地支路除外）最少进行编号 静态优化编号 动态优化编号

122 第五章 P-Q分解法 P-Q分解法是牛顿-拉夫逊法潮流计算的一种简化方法。
牛顿-拉夫逊法的缺点：牛顿-拉夫逊法的雅可比矩阵在每一次迭代过程中都有变化，需要重新形成和求解，这占据了计算的大部分时间，成为牛顿-拉夫逊法计算速度不能提高的主要原因。 P-Q分解法利用了电力系统的一些特有的运行特性，对牛顿-拉夫逊法做了简化，以改进和提高计算速度。

123 牛顿-拉夫逊法简化形成P-Q分解法的过程
牛顿-拉夫逊法修正方程展开为： 根据电力系统的运行特性进行简化： 考虑到电力系统中有功功率分布主要受节点电压相角的影响，无功功率分布主要受节点电压幅值的影响，所以可以近似的忽略电压幅值变化对有功功率和电压相位变化对无功功率分布的影响，即：

124 根据电力系统的正常运行条件还可作下列假设：
电力系统正常运行时线路两端的电压相位角一般变化不大（不超过10~20度）； 电力系统中一般架空线路的电抗远大于电阻； 节点无功功率相应的导纳Q/U*U远小于该节点的自导纳的虚部。 用算式表示如下：

125 由以上假设，可得到雅可比矩阵的表达式为：
修正方程式为： U为节点电压有效值的对角矩阵，B为电纳矩阵（由节点导纳矩阵中各元素的虚部构成）

126 根据不同的节点还要做一些改变： 在有功功率部分，要除去与有功功率和电压相位关系较小的因素，如不包含各输电线路和变压器支路等值Π型电路的对地电纳。 在无功功率部分，PV节点要做相应的处理。 则修正方程表示为： 一般，由于以上原因，B’和B’’是不相同的，但都是对称的常数矩阵 。

127 P-Q分解法的特点： 以一个n-1阶和一个n-m-1阶线性方程组代替原有的2n-m-1阶线性方程组；
修正方程的系数矩阵B’和B”为对称常数矩阵，且在迭代过程中保持不变； P-Q分解法具有线性收敛特性，与牛顿-拉夫逊法相比，当收敛到同样的精度时需要的迭代次数较多； P-Q分解法一般只适用于110KV及以上电网的计算。因为35KV及以下电压等级的线路r/x比值很大，不满足上述简化条件，可能出现迭代计算不收敛的情况。

128 第六章 直流法潮流计算 直流法的特点：简单、计算工作量小、没有收敛性问题，易于快速地处理投入或断开线路等操作。广泛应用于电力系统规划、静态安全分析以及牛顿-拉夫逊法潮流的初值计算等需要大量计算或运行条件不十分理想的场合。 直流法的适用范围：110KV以上的超高压线路。 直流法的经常处理的问题：处理开断问题，例如，在电力系统规划和电力系统静态安全分析时，需要进行一种所谓N-1校核计算，即对于某一种运行方式要逐一开断系统中的线路或变压器，检查是否存在支路过载情况。

129 直流法计算潮流的过程 电力网中每条支路i-j中通过的有功功率为： 根据电力系统的实际条件可做如下假设：
实际电力系统中输电线路（或变压器）的电阻远小于其电抗，对地电导可忽略不计 在正常运行时线路两端相位差很少超过20° 节点电压值的偏移很少超过10%，且对有功功率分布影响不大

130 用式子表示： 从而可得： 各节点的注入功率为与该节点相连各支路功率之和：

131 令B0表示正常运行时电力网节点导纳矩阵的负数，则所有节点注入功率可用矩阵表示为：
解方程求出各节点的相角后，可利用前面的式子求出各支路的有功潮流。 直流法称呼的说明。