4. 节点导纳矩阵的形成

节点导纳矩阵 Y_bus 是潮流计算的核心数据结构，它把网络拓扑、支路参数和变压器变比全部压缩进一个 n×n 的复数矩阵中。对多电压等级电网而言，节点编号方式、支路追加顺序以及变压器非标准变比的处理方式，都会直接影响矩阵的数值结构和后续求解效率。本节围绕这四个关键点展开。

一、多电压等级网络的节点编号

多电压等级电网中，节点分布在 500 kV、220 kV、110 kV、35 kV 等不同电压层级，节点编号不是简单的顺序排列，而是需要兼顾物理分区与计算效率。工程上常用的编号策略是：先按电压等级分区，再在区内按厂站顺序编号。例如将 500 kV 节点编为 1～20，220 kV 节点编为 21～80，110 kV 节点编为 81～200，低压节点依次后延。这样编号的好处是：同一电压等级的节点在矩阵中聚集，非零元素呈现明显的块状分布，便于分块处理和并行计算。

编号时还需注意两点：一是平衡节点（松弛节点）通常编为 1 号，便于在迭代中优先处理；二是PV 节点与 PQ 节点应尽量分开编号，使雅可比矩阵中对应子块结构清晰。编号一旦确定，全网节点数 n 即确定，Y_bus 为 n×n 复对称矩阵。

二、支路追加法形成导纳矩阵

支路追加法是最直观的 Y_bus 形成方法：从空矩阵开始，逐条支路把其导纳贡献叠加到对应位置。对一条连接节点 i、j 的支路，设其串联导纳为 y_ij，对地并联导纳为 y_i0、y_j0，则追加规则为：

Y[i][i] += y_ij + y_i0
Y[j][j] += y_ij + y_j0
Y[i][j] -= y_ij
Y[j][i] -= y_ij

若支路为变压器支路，则不能直接套用上式，需要先做非标准变比处理（见下节）。支路追加法的优点是逻辑简单、易于编程，且天然保证矩阵对称性；缺点是每次追加都要访问矩阵元素，对大规模网络效率一般。实际程序中常配合稀疏存储（如 CSR 格式）使用。

三、变压器非标准变比的处理

实际变压器变比往往不等于理想变比，例如 220/121 kV 而非 220/110 kV，此时变比 k = 220/121 ≈ 1.818，而标准变比 k_N = 220/110 = 2.0，非标准变比 k* = k/k_N ≈ 0.909。处理方式是把变压器等效为理想变压器 + 串联阻抗，并将非标准变比折算到一侧。

设变压器阻抗归算到 j 侧为 z_T，导纳 y_T = 1/z_T，非标准变比 k* 折算到 i 侧，则追加到 Y_bus 的四个元素为：

Y[i][i] += y_T / (k*)^2
Y[j][j] += y_T
Y[i][j] -= y_T / k*
Y[j][i] -= y_T / k*

注意此时 Y[i][j] 与 Y[j][i] 仍相等，矩阵保持对称。若把 k* 折算到 j 侧，则公式中 k* 的位置互换。工程上还常把非标准变比分解为“标准变比 + 附加理想变压器”，附加变比部分单独处理，便于与分接头调节联动。

四、导纳矩阵的稀疏性与对称性

Y_bus 有两个重要性质。第一是对称性：只要支路模型是互易的（线路、变压器均满足），就有 Y_ij = Y_ji。这意味着只需存储上三角或下三角部分，存储量减半。第二是稀疏性：一个节点通常只与 2～5 个相邻节点有支路连接，因此每行非零元素很少。对 n 个节点的网络，非零元素总数约为 2b + n（b 为支路数），而矩阵总元素为 n²。当 n = 1000、b = 1500 时，非零元素约 4000 个，仅占 0.4%。

稀疏性带来两个直接好处：一是存储效率高，采用 CSR（压缩稀疏行）格式只需存储非零值和列索引；二是求解效率高，稀疏 LU 分解或前代回代的计算量远小于稠密矩阵。但要注意，节点编号顺序会显著影响稀疏分解产生的注入元素数量。若编号混乱，原本稀疏的矩阵在分解后可能产生大量填充，导致效率下降。因此工程上常采用最小度排序、AMD 排序等优化编号策略。

五、一个 5 节点算例的追加过程

以图中 5 节点网络为例：节点 1、2 为 500 kV，节点 3、4 为 220 kV，节点 5 为 110 kV。支路包括 1-2 线路、1-3 变压器、2-4 变压器、3-4 线路、3-5 线路、4-5 线路。设线路串联导纳均为 y = 1/(0.02 + j0.06) ≈ 5 − j15（标幺值，示例假设），变压器 1-3 非标准变比 k* = 0.909，变压器 2-4 非标准变比 k* = 1.0（标准变比）。

按支路追加法，先处理线路 1-2：Y[1][1] += y，Y[2][2] += y，Y[1][2] -= y，Y[2][1] -= y。再处理变压器 1-3：Y[1][1] += y_T/k*²，Y[3][3] += y_T，Y[1][3] -= y_T/k*，Y[3][1] -= y_T/k*。依次追加全部支路后，得到完整的 5×5 矩阵。最终矩阵中，节点 1 与 2、3 有连接，节点 5 只与 3、4 有连接，因此第 5 行只有 3 个非零元素，稀疏特征明显。

完成追加后，可验证：Y[1][3] 与 Y[3][1] 是否相等（对称性），Y[5][1]、Y[5][2] 是否为零（稀疏性）。若两者均满足，说明矩阵形成正确，可进入下一步潮流迭代。若发现 Y[1][3] ≠ Y[3][1]，应检查变压器非标准变比是否只折算了一侧，或是否误用了非互易模型。

总结来说，多电压等级网络 Y_bus 的形成可以概括为：合理编号 → 逐支路追加 → 变压器非标准变比折算 → 利用对称性与稀疏性存储和求解。编号决定稀疏结构，追加法决定填充顺序，非标准变比决定变压器支路的四个元素，而对称性与稀疏性则决定后续计算能否高效进行。四者环环相扣，缺一不可。