3. 电力系统网络方程与潮流基础

大扰动功角稳定评估的起点，是把物理电网翻译成一组可计算的代数方程。无论后续使用暂态稳定时域仿真、能量函数法还是直接法，第一步都要建立网络方程，并在此基础上求解潮流，得到扰动前的运行点。本节围绕三个核心对象展开：节点导纳矩阵 Ybus、节点功率平衡方程，以及潮流结果如何作为稳定评估的初始条件接口。

一、节点导纳矩阵：网络的代数骨架

设网络有 n 个节点，节点电压列向量为 V，注入电流列向量为 I，则网络方程写成 I = Ybus · V。其中 Ybus 是 n×n 复对称矩阵，其元素由支路参数直接组装：

• 对角元 Yii 为节点 i 的自导纳，等于接于该节点的所有支路导纳与对地导纳之和；
• 非对角元 Yij 为节点 i、j 之间的互导纳，等于两节点间支路导纳的负值，无支路相连时为 0。

以一条 π 型等值线路为例，支路串联阻抗 z = r + jx，导纳 y = 1/z，两端对地并联电纳分别为 jb/2。则该支路对 Ybus 的贡献为：Yii += y + jb/2，Yjj += y + jb/2，Yij -= y，Yji -= y。变压器支路则需按变比 k 折算，非标准变比侧的自导纳与互导纳都要乘以 1/k 或 1/k² 因子。

需要特别强调：Ybus 只反映网络拓扑与元件参数，与发电机、负荷的注入功率无关。因此当网络结构不变、仅调整发电出力或负荷水平时，Ybus 无需重新组装，这是潮流与稳定计算可以共享同一网络模型的基础。

二、节点功率平衡方程：从电流到功率

实际电力系统中，节点注入通常以功率而非电流给出。将 I = Ybus·V 两边左乘 V 的共轭并取节点形式，可得节点 i 的复功率注入：

Si = Pi + jQi = Vi · Σ(j=1..n) conj(Yij · Vj)

展开为实部与虚部，即得到有功、无功功率平衡方程：

Pi = Vi · Σ Vj · (Gij·cosθij + Bij·sinθij)
Qi = Vi · Σ Vj · (Gij·sinθij − Bij·cosθij)

其中 θij = θi − θj，Gij、Bij 分别为 Yij 的电导与电纳。这两个方程是潮流计算的核心约束：每个节点的注入功率必须等于网络从该节点吸收的功率。

按已知量不同，节点分为三类：PQ 节点给定 P、Q，待求 V、θ；PV 节点给定 P、V，待求 Q、θ；平衡节点给定 V、θ，待求 P、Q，承担全网功率缺额。未知量总数与方程数相等，潮流问题才有确定解。

三、潮流计算：牛顿-拉夫逊法的实现要点

工程上最常用的是牛顿-拉夫逊法。设待求变量为 θ（除平衡节点外）和 V（PQ 节点），不平衡量 ΔP、ΔQ 由给定值与计算值之差构成。迭代方程为：

[ΔP]   [H  N] [Δθ]
[ΔQ] = [J  L] [ΔV/V]

雅可比矩阵 J 的元素由功率方程对 θ、V 求偏导得到，例如 Hij = ∂Pi/∂θj、Nij = ∂Pi/∂Vj。每次迭代解一次线性方程组，更新 θ、V，直到最大不平衡量小于收敛容差（典型取 1e-5 p.u.）。

下面给出一个 3 节点示例的组装与求解思路，参数为合理假设值：节点 1 为平衡节点 V=1.05∠0°；节点 2 为 PV 节点 P=0.5、V=1.02；节点 3 为 PQ 节点 P=−0.8、Q=−0.3。线路阻抗均取 z=0.02+j0.06，对地电纳忽略。按前述规则组装 Ybus 后，用牛顿法迭代 3～5 次即可收敛，得到各节点电压幅值与相角。

四、潮流结果与稳定评估的接口

潮流解是稳定评估的初始运行点，接口关系体现在三方面：

• 电压与相角：各节点 Vi、θi 直接作为暂态仿真 t=0 时刻的状态初值；
• 发电机内电势：由 Eq' = V + jXd'·I 反推，得到各发电机暂态电抗后的电动势与功角初值 δ0；
• 负荷等值：潮流中的恒功率负荷在稳定计算中可转为恒阻抗或恒电流模型，其等值导纳并入 Ybus，形成故障前网络。

因此，潮流不收敛或收敛到异常解，会直接导致稳定评估的初始功角、初始电磁功率失真，进而影响极限切除时间与稳定裕度的判断。工程实践中，通常先确认潮流收敛且运行点合理，再进入稳定计算环节。

五、工程实现中的关键检查点

把上述流程落到代码或仿真工具时，建议按以下顺序核对，避免“能算但算错”：

完成上述检查后，潮流解即可作为稳定评估的输入。需要再次明确：Ybus 是网络的结构性描述，功率平衡方程是节点注入的约束，潮流计算是求解这两者联立方程的过程，而稳定评估接口则负责把潮流解转换为发电机内电势、功角初值和负荷等值模型。四者环环相扣，任一环节出错都会传导到最终的功角稳定判断上。