4. 电力系统潮流计算数学模型

潮流计算是电力系统稳态分析的核心工具，其本质是在给定网络拓扑、支路参数和部分节点运行条件的前提下，求解全网各节点的电压幅值与相角，进而得到各支路的功率分布和网络损耗。要建立可求解的数学模型，必须回答两个基本问题：每个节点已知什么、待求什么，以及节点功率与电压之间满足怎样的方程关系。前者对应节点分类，后者对应功率方程。

一、节点功率方程的基本形式

设网络共有 n 个节点，节点 i 的注入电流与节点电压通过节点导纳矩阵 Y 关联：

İ_i = Σ_j Y_ij · V̇_j    (j = 1, 2, …, n)

节点复功率为 Ṡ_i = V̇_i · İ_i*，展开后得到节点注入功率与电压的关系：

Ṡ_i = P_i + jQ_i = V̇_i · (Σ_j Y_ij · V̇_j)*

将 Y_ij = G_ij + jB_ij 代入，并分别取实部与虚部，即得直角坐标下的功率方程：

P_i = e_i · Σ_j (G_ij·e_j − B_ij·f_j) + f_i · Σ_j (G_ij·f_j + B_ij·e_j)
Q_i = f_i · Σ_j (G_ij·e_j − B_ij·f_j) − e_i · Σ_j (G_ij·f_j + B_ij·e_j)

其中 V̇_i = e_i + jf_i。若改用极坐标 V̇_i = V_i∠θ_i，则功率方程写为：

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)

式中 θ_ij = θ_i − θ_j。两种坐标形式在数学上等价，但数值特性不同：直角坐标下功率方程是电压实部、虚部的二次多项式，雅可比矩阵元素计算简单，但需额外处理电压幅值约束；极坐标下 P、Q 方程对 θ 和 V 的偏导数具有更规整的稀疏结构，且 PV 节点的电压幅值可直接作为已知量固定，因此工程中极坐标牛顿–拉夫逊法应用更为广泛。

二、节点分类的数学描述

每个节点有四个运行变量：P_i、Q_i、V_i、θ_i。由于每个节点可列出两个功率方程（有功、无功），因此必须已知其中两个变量，才能求解另外两个。根据已知量的不同，节点分为三类：

PQ 节点是最常见的类型。负荷节点的有功、无功由负荷预测或实测给出，发电功率固定的节点也可归入此类。其数学描述为：

P_i = P_i^sp （给定值）
Q_i = Q_i^sp （给定值）
待求：V_i、θ_i

PV 节点对应具有自动电压调节器的发电机组。机组有功出力由调度计划确定，机端电压由励磁系统维持为设定值，而无功出力取决于系统需求，是待求量。其数学描述为：

P_i = P_i^sp （给定值）
V_i = V_i^sp （给定值）
待求：Q_i、θ_i

需要注意，PV 节点的无功出力必须满足机组容量约束 Q_i^min ≤ Q_i ≤ Q_i^max。若迭代过程中 Q_i 越限，则该节点应转换为 PQ 节点，将 Q_i 固定为限值，V_i 转为待求量，这是潮流计算中重要的节点类型转换逻辑。

平衡节点全系统通常只设一个。其电压幅值和相角均给定，有功、无功均待求。平衡节点承担两个职能：一是为全网相角提供参考（θ_s = 0），使潮流方程有唯一解；二是平衡全网功率损耗和未知的功率差额，其注入功率由计算结果反推得到。选择平衡节点时，一般取主调频厂母线或与外部系统相连的联络母线。

三、约束条件与变量匹配

潮流方程本身只描述网络物理规律，实际运行还必须满足一系列不等式约束。这些约束决定了求解结果的可行域，也是节点类型转换的触发条件。

• 发电机有功约束：P_i^min ≤ P_i ≤ P_i^max，通常由调度计划给定，潮流计算中作为已知量处理。
• 发电机无功约束：Q_i^min ≤ Q_i ≤ Q_i^max，PV 节点迭代中若 Q_i 越限，节点转为 PQ 节点。
• 电压幅值约束：V_i^min ≤ V_i ≤ V_i^max，一般取 0.95～1.05 p.u.，部分节点可放宽至 0.90～1.10 p.u.。
• 支路传输功率约束：S_ij ≤ S_ij^max，用于校验线路或变压器是否过载。
• 平衡节点功率约束：平衡节点的 P_s、Q_s 也需在合理范围内，否则说明系统功率平衡存在较大缺口。

变量匹配的核心逻辑可以概括为：每个节点四个变量中已知两个、待求两个，全网 2n 个方程对应 2n 个未知量。但平衡节点的存在使这一计数需要修正——平衡节点的 V、θ 已知，P、Q 待求，相当于额外提供了两个已知量，同时其功率方程不参与迭代求解，而是用于事后计算。因此实际参与迭代的方程数为 2(n−1)，未知量也为 2(n−1)，方程组恰好闭合。

四、坐标形式的选择与比较

直角坐标与极坐标潮流方程在工程实现中各有适用场景。直角坐标下，电压实部 e 和虚部 f 的初值容易设定（如 e=1.0、f=0.0），功率方程是二次多项式，雅可比矩阵元素不含三角函数，计算速度快；但 PV 节点的电压幅值约束 V_i² = e_i² + f_i² 需要作为额外方程处理，增加了方程维数。极坐标下，PV 节点的 V_i 直接作为已知量，方程结构更简洁，且 P–θ、Q–V 之间的耦合关系与电力系统物理直觉一致，便于分析电压稳定性和无功平衡，因此成为大多数商业潮流程序的首选。