4. 电力系统稳态分析与潮流计算基础

潮流计算是电网调度运行与方式安排中最基础、使用最频繁的计算工具。无论是编制年度运行方式、安排检修计划,还是校核断面输送能力、评估新能源接入影响,最终都要落到一组稳态潮流方程上。本节从节点分类出发,逐步建立功率方程,再介绍牛顿-拉夫逊法与PQ分解法的求解思路,并说明它们在方式安排中的实际用法。

一、节点分类:把电网抽象成四类已知量

潮流计算的第一步是把实际电网抽象为节点和支路。每个节点有四个状态量:电压幅值 U、电压相角 θ、注入有功 P、注入无功 Q。四个量中必须已知两个、求解两个,否则方程无法定解。按已知量的不同,节点分为三类:

  • PQ节点:已知 P、Q,待求 U、θ。绝大多数变电站母线、负荷节点属于此类。负荷功率由预测或实测给出,视为已知注入。
  • PV节点:已知 P、U,待求 Q、θ。装有自动电压调节器的发电厂高压母线属于此类,有功由调度计划给定,电压由励磁系统维持。
  • 平衡节点(Vθ节点):已知 U、θ(通常取 θ=0 作为相角参考),待求 P、Q。全系统只设一个,一般选容量最大、调频能力最强的主调电厂母线,它承担全网功率不平衡量。

设系统共 n 个节点,其中 PQ 节点 m 个、PV 节点 n−m−1 个、平衡节点 1 个。待求量为 m 个电压幅值和 n−1 个相角,共 n+m−1 个未知量,与方程数一致。这个计数关系是判断潮流问题是否可解的基本依据,也是方式安排中确定"哪些量必须给定"的依据。

二、功率方程:节点注入的约束关系

对任一节点 i,注入电流与节点电压通过节点导纳矩阵 Y 联系:Ii = ΣYijUj。节点复功率为 Si = UiIi*,展开为实部与虚部即得功率方程:

P_i = U_i · Σ_j U_j (G_ij·cosθ_ij + B_ij·sinθ_ij)
Q_i = U_i · Σ_j U_j (G_ij·sinθ_ij − B_ij·cosθ_ij)
其中 θ_ij = θ_i − θ_j,G_ij、B_ij 为导纳矩阵元素

这是一组非线性代数方程,Pi、Qi 与 U、θ 之间是三角函数与乘积关系,无法直接解析求解,必须迭代。方式安排中常见的"某断面潮流是否越限""某机组出力调整后电压如何变化",本质上都是在这组方程上做扰动分析。

三、牛顿-拉夫逊法:以雅可比矩阵为核心的迭代

牛顿-拉夫逊法把功率方程写成偏差形式 f(x)=0,在初值 x(0) 处泰勒展开并保留一阶项,得到修正方程 J·Δx = −f(x),其中 J 为雅可比矩阵。对 PQ 节点,偏差量为 ΔP、ΔQ;对 PV 节点,偏差量只有 ΔP,ΔQ 由电压约束替代。雅可比矩阵元素由功率方程对 U、θ 求偏导得到,例如:

∂P_i/∂θ_j = U_i·U_j·(G_ij·sinθ_ij − B_ij·cosθ_ij)   (i≠j)
∂P_i/∂U_j = U_i·(G_ij·cosθ_ij + B_ij·sinθ_ij)        (i≠j)
∂Q_i/∂θ_j = −U_i·U_j·(G_ij·cosθ_ij + B_ij·sinθ_ij)   (i≠j)
∂Q_i/∂U_j = U_i·(G_ij·sinθ_ij − B_ij·cosθ_ij)        (i≠j)

迭代流程为:给定初值(通常取 U=1.0 p.u.、θ=0)→ 计算 ΔP、ΔQ → 形成雅可比矩阵 → 解线性方程得 ΔU、Δθ → 更新电压 → 判断最大偏差是否小于收敛精度(一般取 10−5~10−4 p.u.)。牛顿法具有平方收敛特性,通常 3~5 次迭代即可收敛,但每次迭代都要重新形成并分解雅可比矩阵,单次计算量大。

四、PQ分解法:利用电网物理特性的简化

高压输电网有两个显著特点:支路电抗远大于电阻(X≫R),且节点电压相角差很小(θij≈0)。由此可得两个近似:有功主要与相角相关,无功主要与电压幅值相关。PQ分解法据此把修正方程拆成两个低阶方程交替求解:

ΔP/U = B' · (U·Δθ)
ΔQ/U = B'' · ΔU

其中 B'、B'' 为常数矩阵,仅由网络电抗参数构成,迭代过程中不变,只需分解一次。每次迭代交替求解有功—相角、无功—电压两个子问题,单次计算量大幅下降。代价是收敛速度由平方降为线性,迭代次数增至 8~15 次,但总耗时通常仍明显低于牛顿法,因此成为在线潮流和方式计算的主力算法。

潮流计算节点分类与求解方法对照 节点分类(已知量 → 待求量) PQ 节点 已知 P、Q 待求 U、θ 负荷母线、变电站 PV 节点 已知 P、U 待求 Q、θ 调压发电厂母线 平衡节点 Vθ 已知 U、θ(θ=0) 待求 P、Q 主调电厂,全系统唯一 求解方法对比 牛顿-拉夫逊法 修正方程:J·Δx = −f(x) 雅可比矩阵 J 每次迭代更新 收敛特性:平方收敛 迭代次数:3~5 次 单次计算量大,适合离线精确计算 适用:方式编制、详细校核、小系统 收敛精度:10⁻⁵~10⁻⁴ p.u. PQ 分解法 ΔP/U = B'·(U·Δθ) ΔQ/U = B''·ΔU B'、B'' 常数矩阵,仅分解一次 收敛特性:线性收敛 迭代次数:8~15 次 适用:在线潮流、大系统、实时分析 前提:X≫R,θ_ij 很小

五、在方式安排中的应用

方式安排中,潮流计算不是孤立的一次计算,而是围绕"给定运行条件—求解稳态分布—校核约束"的闭环。具体做法是:

  1. 确定计算边界:根据方式类型(夏大、冬小、检修方式等)给定各节点负荷、各机组出力计划、平衡节点电压、PV节点电压设定值。这些量就是潮流计算的输入。
  2. 选择算法:离线方式编制通常用牛顿法,收敛可靠、精度高;在线安全分析、实时调度辅助决策多用PQ分解法,速度快、可反复计算。
  3. 求解并检查收敛性:若迭代不收敛,常见原因是初值不合理、网络参数错误或方式本身无解(如无功严重不足导致电压崩溃)。此时应检查平衡节点选取、PV节点电压设定和负荷水平。
  4. 校核约束:重点看三类量——线路和变压器是否过载(S ≤ Smax)、节点电压是否越限(如 220 kV 母线 0.95~1.05 p.u.)、平衡节点出力是否超出计划范围。
  5. 调整与迭代:若越限,通过调整机组出力、投切无功补偿、改变变压器分接头或调整运行方式,重新计算,直到所有约束满足。
方式安排中的实用判断:当某断面潮流接近限额时,可先做灵敏度分析——计算该断面有功对相关机组出力的灵敏度系数,再按灵敏度大小排序调整,比反复试算效率高得多。PQ分解法中的 B' 矩阵本身就包含了网络电抗信息,可用于快速估算调整效果。
常见误区:把平衡节点选在容量小的电厂,会导致该节点出力剧烈波动甚至超出机组调节范围,计算结果失去参考意义;PV节点电压设定值若超出实际励磁调节能力,潮流结果虽然收敛,但方式不可执行。节点分类和参数设定必须与实际调度控制方式一致。

掌握节点分类、功率方程、牛顿法与PQ分解法的原理和差异,是读懂潮流计算结果、判断方式是否可行的前提。实际工作中,调度员不需要手算雅可比矩阵,但必须清楚每个输入量的来源和含义,才能在结果异常时快速定位问题,做出正确的方式调整决策。