2. 交流输电网络数学模型
潮流计算的第一步,是把一张由线路、变压器、并联补偿设备组成的物理电网,翻译成计算机能够求解的代数方程组。这个翻译过程的核心产物就是节点导纳矩阵 Ybus。理解 Ybus 的每一个元素从哪里来、代表什么,是后续牛顿-拉夫逊法、PQ 分解法能够正确收敛的前提。本节从元件级等值电路出发,逐层搭建到全网矩阵。
一、线路的π型等值电路
输电线路在潮流计算中普遍采用π型等值电路:把线路的串联阻抗集中在中间支路,把线路对地电容平均分成两半,分别挂在两端节点上。对于一条长度为 l 的线路,已知单位长度参数 r、x、b(b 为对地电纳),则:
Z = (r + jx) · l # 串联阻抗
Y_c = jb · l # 全线路对地电纳
Y/2 = Y_c / 2 # 每端对地电纳
以一条 220 kV 线路为例,假设 r=0.08 Ω/km、x=0.40 Ω/km、b=2.8×10⁻⁶ S/km、l=100 km,则 Z=8+j40 Ω,Y_c=j2.8×10⁻⁴ S,每端对地电纳为 j1.4×10⁻⁴ S。π型电路的两端对地电纳在 Ybus 中直接叠加到对应节点的自导纳上,而串联支路则同时影响两个节点的自导纳和它们之间的互导纳。
需要特别说明的是,π型电路是集中参数近似。当线路超过约 250 km 时,分布参数效应明显,应改用双曲函数修正或等值π型参数,否则潮流结果会出现可察觉的偏差。
二、节点导纳矩阵的形成规则
设网络有 n 个节点,Ybus 是一个 n×n 的复数方阵。其元素按以下规则累加:
- 互导纳 Yij(i≠j):等于连接节点 i 与节点 j 之间所有支路导纳之和的负值。若两节点间无直接支路,则 Yij=0。
- 自导纳 Yii:等于与节点 i 相连的所有支路导纳之和,再加上节点 i 上挂接的对地导纳(线路π型对地电纳、并联电容器、电抗器、变压器励磁支路等)。
用一句话概括:互导纳看“连不连”,自导纳看“连多少、对地挂多少”。由于输电网络中一个节点通常只与少数几个节点相连,Ybus 是高度稀疏的矩阵,这也是稀疏矩阵技术在潮流计算中被广泛使用的原因。
下面这张图给出了从元件等值电路到 Ybus 元素的映射关系,以及变压器非标准变比支路的处理方式:
三、变压器非标准变比的处理
实际变压器的高低压侧额定电压比,往往不等于系统基准电压比。例如一台 220/121 kV 的变压器,若系统基准为 220/110 kV,则其非标准变比 k = (220/121)/(220/110) = 110/121 ≈ 0.909。这个 k 不能忽略,否则潮流结果中会出现虚假的电压偏差。
处理方法是把变压器等值为“理想变压器 + 漏抗”串联支路。设漏抗为 ZT,变比 k 折算到节点 p 侧,则对 Ybus 的贡献为:
Y_pp = 1/Z_T
Y_pq = -1/(k·Z_T)
Y_qp = -1/(k·Z_T)
Y_qq = 1/(k²·Z_T)
注意 Ypq 与 Yqp 相等,矩阵仍然对称;但 Ypp 与 Yqq 不再相等,这是非标准变比带来的直接后果。若 k=1,则退化为普通串联支路。工程上还有一种做法是把非标准变比“归算”到某一侧,等效为在节点上附加一个对地导纳,两种做法在数学上等价,但前者更直观、更不易出错。
四、对地电容与并联补偿设备建模
线路π型电路的两端对地电纳、变压器励磁支路、并联电容器和并联电抗器,都属于对地并联支路。它们只影响所在节点的自导纳,不影响互导纳。具体地:
- 并联电容器:注入无功 QC(容性),其等值导纳为 jQC/U²,直接加到 Ymm 上。
- 并联电抗器:吸收无功 QL(感性),等值导纳为 −jQL/U²,同样加到 Ymm 上。
- 线路对地电纳:Yc/2 分别加到两端节点的自导纳。
需要强调的是,并联补偿设备在潮流计算中有两种处理方式:一是作为恒定导纳并入 Ybus,其注入无功随电压平方变化;二是作为恒定无功注入放在节点功率方程中。前者适用于电容器组、电抗器等无调节设备,后者适用于可调无功补偿或 SVC 等。选择哪种方式,取决于设备实际运行特性,不能随意混用。
五、一个可手算的示例
设两节点系统:节点1与节点2之间有一条线路,Z=0.02+j0.06(标幺值),线路总对地电纳 Yc=j0.10;节点2上挂一台并联电容器,容量对应导纳 j0.05。则:
Y_12 = Y_21 = -1/(0.02+j0.06) = -5 + j15
Y_11 = 1/(0.02+j0.06) + j0.05 = 5 - j15 + j0.05 = 5 - j14.95
Y_22 = 1/(0.02+j0.06) + j0.05 + j0.05 = 5 - j14.90
可以看到,节点2因为多挂了并联电容器,其自导纳的虚部比节点1更“正”一些,这正是容性补偿在矩阵中的体现。若把该电容器改为电抗器,则虚部会向负方向变化。通过观察 Ybus 对角元虚部的变化,可以快速判断某节点上并联设备的性质,这也是调试潮流程序时一个实用的检查手段。
掌握以上四类元件的建模规则后,任何规模的输电网络都可以按“逐支路累加”的方式生成 Ybus。下一阶段的任务,就是在这个矩阵基础上建立节点功率方程,并选择合适的方法求解。