第4节 · 重点难点解析

Deep Dive

重点难点深度解析

逐层展开：概念 → 原理 → 推导 → 应用，建议按顺序阅读。

01节点功率方程的建立与 2n 变量匹配▼

概念：潮流计算的对象是「稳态、对称、三相对称运行」的电网。把每个节点看成电流的汇集点，用节点电压和节点导纳矩阵就能完整描述网络。

原理：由基尔霍夫电流定律（KCL）在节点 i 处写成 I_i = Σ_j Y_ij V_j，而节点的注入复功率为 S_i = V_i I_i^*，代入并对电流取共轭即得节点功率方程。

节点功率方程（复数形式，标幺值） S_i = V_i Σ_j=1ⁿ Y_ij^* V_j^*　　(i = 1, 2, …, n)

推导：把 V_i = V_i∠θ_i、Y_ij = G_ij + jB_ij、θ_ij = θ_i − θ_j 代入并分出实部与虚部，就得到 2n 个实变量方程。

实变量形式（极坐标） P_i = V_i Σ_j V_j ( G_ij cosθ_ij + B_ij sinθ_ij )
Q_i = V_i Σ_j V_j ( G_ij sinθ_ij − B_ij cosθ_ij )

变量匹配：全网 2n 个实变量（V_i、θ_i）。PQ 节点给 P_i、Q_i（2 个条件），PV 节点给 P_i、V_i（2 个条件），平衡节点给 V_s、θ_s（2 个条件）。因此若共 m 个 PQ、r 个 PV，则待求量 = 2m + r，方程数也恰为 2m + r，方程组封闭可解。

易错点：① 忘记平衡节点要给定 θ（相角参考），否则全网相角无参考、方程奇异；② 把平衡节点的 P_s、Q_s 也当作已知量输入——它们是计算结果，用来吸收网损与功率不平衡；③ 标幺制下同一天然功率基准必须统一。

应用：节点功率方程是所有数值解法的共同出发点。程序实现时先用「支路表 + 变压器表」形成 Y 阵，再按节点类型装配方程。

02牛顿-拉夫逊法：修正方程、雅可比矩阵与二次收敛▼

概念：牛顿-拉夫逊法把非线性方程组在初值附近线性化，逐次求解修正量，直到失配量小于收敛判据。它是潮流计算中收敛最快、最可靠的方法。

原理：设待求量列向量为 x（θ 与 V），失配量函数 f(x) = 0。在初值 x⁽⁰⁾ 处作泰勒展开并忽略二阶及以上项：

泰勒展开与修正方程 f(x⁽⁰⁾ + Δx) ≈ f(x⁽⁰⁾) + f′(x⁽⁰⁾) Δx = 0
⟹ Δf = J Δx，其中 J = ∂f/∂x 为雅可比矩阵

推导（潮流的具体形式）：把方程按「有功失配—无功失配」与「相角修正—电压修正」分块：

极坐标牛顿法修正方程 [ ΔP ; ΔQ ] = [ H N ; M L ] · [ Δθ ; ΔV/V ]
其中 ΔP_i = P_i^sp − P_i(V, θ)，ΔQ_i = Q_i^sp − Q_i(V, θ)

雅可比矩阵各元素（θ_ij = θ_i − θ_j）：

雅可比矩阵元素（极坐标，i ≠ j） H_ij = V_i V_j ( G_ij sinθ_ij − B_ij cosθ_ij )　　N_ij = V_i ( G_ij cosθ_ij + B_ij sinθ_ij )
M_ij = −V_i V_j ( G_ij cosθ_ij + B_ij sinθ_ij )　　L_ij = −V_i ( G_ij sinθ_ij − B_ij cosθ_ij )

对角元素 H_ii = −Q_i − V_i²B_ii　　N_ii = P_i + V_i²G_ii
M_ii = P_i − V_i²G_ii　　L_ii = −Q_i − V_i²B_ii

迭代流程：① 平启动 V_i = 1、θ_i = 0（PV 节点取给定值）；② 算 P、Q 与失配量；③ 判 max|ΔP, ΔQ| < ε；④ 形成并解雅可比矩阵得修正量；⑤ 更新 θ_i ← θ_i + Δθ_i，V_i ← V_i(1 + ΔV_i/V_i)；⑥ 回到 ②。

二次收敛：在解附近每迭代一次失配量约按平方衰减（例如 10⁻¹ → 10⁻² → 10⁻⁴ → 10⁻⁸）。代价是每次迭代都要重新形成并分解雅可比矩阵。

易错点：① 修正量用的是 ΔV_i/V_i 而不是 ΔV_i，对应 N、L 列的取值；② 只有 PQ 节点才列 ΔQ 方程、才有 V 修正列，PV 节点的 V 是固定值；③ 雅可比矩阵的元素依赖当前电压，必须每次迭代重新计算。可用数值差分（有限差分）核对解析雅可比，这是查错的有效手段。

应用：牛顿法配合稀疏技术可求解上万节点系统；实际工程软件中常用「牛顿法 + 稀疏因子表 + 节点编号优化」组合。

03高斯-赛德尔法为什么慢：线性收敛与 PV 节点处理▼

概念：高斯-赛德尔法是最早用于潮流的迭代法，逐节点反复代入，直到各节点电压变化量足够小。

迭代格式（逐节点，即时更新） V_i^(k+1) = (1/Y_ii) [ ( S_i / V_i^*(k) ) − Σ_j≠i Y_ij V_j ]

原理：每个节点 i 的方程对 V_i 是「局部线性」的，代入其他节点最新电压即可显式求出 V_i。因为逐点更新，第一次更新的结果立即被后续节点使用，所以收敛比雅可比迭代快一些，但本质仍是线性收敛。

PV 节点处理：PV 节点 Q_i 未知，不能直接使用上式。做法是：① 按给定 P_i 和「上一轮算得的 Q_i」代入迭代式得到电压；② 把结果电压的幅值强制校正为给定值 V_i^sp（保留其相角）；③ 用校正后的电压反算该节点的 Q_i，供下一轮使用。这就是「无功修正与电压幅值校正」。

PV 节点无功增量（第 k 轮） ΔQ_i^(k) ≈ V_i^sp Σ_j V_j^(k) ( G_ij sinθ_ij^(k) − B_ij cosθ_ij^(k) )

收敛特性：设误差为 e^(k)，则有 e^(k+1) ≈ c·e^(k)（c 为常数），即每轮残差按固定比例缩小——线性收敛。本页工具二的 3 节点算例中，残差从 1.50×10⁻¹ 降到 10⁻⁵ 用了 9 次迭代（取 ε = 10⁻⁸ 则需 15 次；初值改为 0.95 时需 11 次），而牛顿法只需 4～5 次，对比非常直观。

易错点：① S_i 是节点注入（发电机 − 负荷），负荷节点应取负号；② 平衡节点的 V、θ 始终不变，不参与更新；③ 收敛判据既可用电压变化量，也可用功率失配量，两种要区分清楚。

应用：赛德尔法程序简单、内存小，适合教学与小型系统；现代工程计算已基本改用牛顿法与 PQ 分解法。

04PQ 分解法：简化假设的边界▼

概念：高压电网中 P—θ 强耦合、Q—V 强耦合，而 P—V、Q—θ 耦合较弱。PQ 分解法正是利用这一物理特性把牛顿法的雅可比矩阵解耦。

简化假设（三条）：① 支路电阻远小于电抗（r ≪ x）；② 相邻节点电压相角差很小，cosθ_ij ≈ 1、sinθ_ij ≈ 0；③ 电压幅值接近 1 且变化量小。

解耦后的修正方程 ΔP/V = B′ Δθ　　（P—θ 迭代，B′ 由 −1/x 组成，忽略电阻与对地支路）
ΔQ/V = B″ ΔV　　（Q—V 迭代，B″ 由 −B_ij 组成，保留对地支路）

推导思路：把雅可比子块 N、M 置零（解耦），再把 H、L 中的 G 项忽略并用 V 归一化，即得常数矩阵 B′、B″。两式交替求解（P—θ 一轮、Q—V 一轮算作一次迭代）。

特点对比：迭代次数一般为牛顿法的 2～4 倍，但每次迭代不需要重新形成和分解矩阵，单次计算量小、内存省，总体速度常优于牛顿法；对 r/x 较大的配电网或重载条件，简化假设不成立，可能收敛困难。

易错点：① B′ 与 B″ 不是同一个矩阵（是否含电阻、是否含对地支路不同）；② 得出的是 Δθ 与 ΔV/V，更新时不要漏乘 V；③ PV 节点不参与 Q—V 迭代。

06结果分析：电压剖面、网损率与稳定裕度▼

概念：潮流算完只是拿到数据，工程上真正要做的是「读结果」：判断电压是否越限、线路是否过载、网损是否合理、无功是否越限、系统是否有裕度。

关键指标：

指标定义电压偏差率 = (V_i − V_N) / V_N × 100%　（判断是否越限）
支路负载率 = S_ij / S_ij^max × 100%　（判断是否过载）
网损率 = ΔP_Σ / P_负荷Σ × 100%　（衡量经济性）
电压稳定裕度 ≈ (V_临界 − V_当前) / V_临界

平衡节点功率的意义：由全网功率平衡，平衡节点注入的有功等于「全网负荷 + 网损 − 其他发电机有功」，因此它自动吸收了网损。网损又随电压水平变化：抬高电压水平可减小电流、降低线路 I²R 损耗，这是「升压降损」的原理。

无功越限与 PV→PQ：若某 PV 节点算出的 Q_Gi 超出 [Q_Gi^min, Q_Gi^max]，说明发电机已无法维持给定电压，应把该节点改为 PQ 节点（Q 固定为限值），重新计算潮流；可能引起电压进一步下降，需迭代判断。

易错点：① 只看电压幅值不看相角差，忽略了线路输送能力；② 把「网损率」与「线损率」混为一谈（后者按电量统计）；③ 忘记检查平衡节点出力是否超出机组容量；④ 认为电压越下限只加电容器就能解决——还要看是否有功不足与线路压降过大。

校核习惯：拿到结果先做三项闭合校验——全网有功平衡、无功平衡、以及各支路损耗之和与平衡节点功率增量一致。本页工具二就自动给出有功平衡校验值。