3. 节点导纳矩阵与网络方程

在电力系统运行方式计算中，无论采用牛顿-拉夫逊法、PQ分解法还是直流潮流，最终都要落到同一个数学核心上：节点导纳矩阵 Y_bus 与由它构成的网络方程 I = Y_busU。理解 Y_bus 的形成规则、稀疏结构和修正方法，是读懂潮流程序、排查不收敛问题、手工校核小系统的基础。

一、节点导纳矩阵的形成规则

设网络有 n 个节点，节点电压列向量为 U = [U₁, U₂, …, U_n]^T，节点注入电流列向量为 I = [I₁, I₂, …, I_n]^T，则网络方程为：

I = Y_bus · U
Y_bus 为 n×n 复数方阵，元素 Y_ij = G_ij + jB_ij

Y_bus 的元素按以下三条规则直接由支路参数填写，不需要列写完整的 KCL 方程组再消元：

• 对角元素（自导纳）Y_ii：等于与节点 i 直接相连的所有支路导纳之和，再加上节点 i 的对地支路（并联电容、电抗器、接地支路）导纳。即 Y_ii = Σ y_ij + y_i0，其中 y_ij 为连接 i、j 的支路导纳，y_i0 为节点 i 的对地导纳。
• 非对角元素（互导纳）Y_ij（i≠j）：等于节点 i 与节点 j 之间支路导纳的负值，即 Y_ij = −y_ij。若 i、j 之间无直接支路，则 Y_ij = 0。
• 对称性：对于不含移相器的普通支路，y_ij = y_ji，因此 Y_bus 是对称矩阵。只有变压器含非标准变比或移相环节时，才会出现 Y_ij ≠ Y_ji 的非对称情况。

以一条 π 型等值线路为例：节点 i、j 之间线路串联阻抗为 z = r + jx，两端各有对地充电电容容纳 b_c/2。则支路导纳 y = 1/z，节点 i 的自导纳为 Y_ii = y + jb_c/2，节点 j 同理，互导纳 Y_ij = Y_ji = −y。可见对地支路只进对角元，串联支路同时进对角元和非对角元，这是填写 Y_bus 时最容易出错的地方。

二、稀疏性特点

实际电网中，一个节点平均只与 2～4 个节点直接相连，而系统节点数可达数百至数万。因此 Y_bus 中绝大多数非对角元素为零，是一个高度稀疏的矩阵。以 IEEE 14 节点系统为例，14×14 = 196 个元素中非零元素仅约 40 余个，稀疏度超过 75%；节点数越大，稀疏度越高。

稀疏性带来两个直接后果：一是存储上只需记录非零元素及其行列位置（如按行压缩存储 CSR），二是求解时采用稀疏三角分解和节点编号优化（如最小度法、AMD 排序），可把计算量从 O(n³) 量级大幅降低。这也是潮流程序能处理上万节点系统的根本原因。需要强调的是，稀疏性由网络拓扑决定，与运行方式无关：只要支路不投退，零元素位置就不变；运行方式变化只改变非零元素的数值。

三、网络方程与节点类型

网络方程 I = Y_busU 本身是线性的，但节点注入电流 I 与节点功率、电压的关系是非线性的：S_i = U_i·I_i^*。因此潮流计算中通常把方程改写为功率形式：

P_i + jQ_i = U_i · Σ (Y_ij · U_j)*   (j = 1…n)

为求解，需按节点给定量的不同划分三类节点：PQ 节点给定 P、Q，待求 U、θ；PV 节点给定 P、U，待求 Q、θ；平衡节点给定 U、θ，待求 P、Q，负责平衡全网功率损耗。n 个节点中，平衡节点 1 个，其余为 PQ 或 PV，未知量个数与方程个数相等，方程组可解。

四、不同接线方式下的矩阵修正

运行方式计算的核心工作之一，就是根据开关状态、变压器分接头、线路投退等改变 Y_bus。修正遵循“支路增删只影响该支路两端节点的相关元素”这一局部性原则：

• 投入一条支路 i–j：Y_ii += y，Y_jj += y，Y_ij −= y，Y_ji −= y。
• 退出一条支路 i–j：按上述量取负号累加，即 Y_ii −= y，Y_ij += y。
• 节点对地支路投切：只修改该节点的 Y_ii，不影响任何非对角元。
• 变压器变比改变：若变比 k 接在 i 侧，则 Y_ii 增加 y/k²，Y_ij、Y_ji 变为 −y/k，Y_jj 增加 y，矩阵不再对称。
• 节点分裂或合并：相当于删除或新增节点，需重新排列矩阵行列，通常配合节点编号优化一起处理。

下面用一张图概括 Y_bus 的填写规则、稀疏结构与接线修正之间的关系。

掌握上述规则后，面对一个具体运行方式，可以按“先列支路表、再逐条累加、最后检查对角元”的顺序手工填写小系统 Y_bus，也可以据此核对程序输出的矩阵是否正确。当潮流不收敛时，首先检查 Y_bus 中是否有误填的对地导纳或漏掉的互导纳，往往能快速定位问题。