4. 潮流计算基础与电压稳定分析的衔接

电压稳定分析并不是脱离潮流计算独立存在的理论,它本质上是在潮流方程解的结构上做文章。要理解“电压崩溃点”为什么对应雅可比矩阵奇异,必须先弄清楚牛顿-拉夫逊潮流在算什么、PV节点为什么会转成PQ节点、连续潮流又是如何绕过临界点继续追踪解的。本节把这几件事串成一条线。

一、牛顿-拉夫逊潮流:把非线性方程线性化迭代

潮流计算的核心是节点功率平衡方程。对每个节点,注入功率由电压幅值和相角决定,写成极坐标形式:

P_i = V_i · Σ V_j (G_ij cosθ_ij + B_ij sinθ_ij)
Q_i = V_i · Σ V_j (G_ij sinθ_ij − B_ij cosθ_ij)
θ_ij = θ_i − θ_j

其中 G、B 来自节点导纳矩阵 Y=G+jB。未知量是各节点的电压幅值 V 和相角 θ。牛顿-拉夫逊法把功率不平衡量 ΔP、ΔQ 与电压修正量 Δθ、ΔV 用雅可比矩阵 J 联系起来:

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

每次迭代求解该线性方程组,得到修正量后更新电压,直到最大不平衡量小于收敛精度(如 1e-5 p.u.)。这里有一个关键点:雅可比矩阵 J 不是常数,它随电压状态变化。正常运行时 J 非奇异,迭代收敛;当系统接近电压崩溃点时,J 逐渐趋于奇异,迭代次数急剧增加甚至发散。这正是电压稳定分析与潮流计算衔接的第一个接口。

二、PV节点与PQ节点:类型转换的物理含义

潮流计算中节点分三类:平衡节点(Slack)给定 V、θ;PV节点给定 P、V;PQ节点给定 P、Q。PV节点通常对应接有调压发电机的母线,发电机通过调节无功维持电压。但发电机无功出力有上限 Qmax 和下限 Qmin。

当某PV节点所需无功超过 Qmax 时,发电机无法继续维持电压,该节点电压将随负荷增加而下降,此时它实际上失去了电压支撑能力,必须转换为PQ节点,Q固定为 Qmax。反之,若无功低于 Qmin,则转为 Q 固定为 Qmin 的PQ节点。这个转换不是数值技巧,而是物理约束的体现。

PV→PQ转换是电压稳定分析中最重要的非线性事件之一。转换发生后,系统的无功支撑减少,雅可比矩阵结构改变,往往使系统更接近临界点。因此在实际分析中,需要监测哪些PV节点在负荷增长过程中先触及无功限值。

转换的判断逻辑可以概括为:每轮潮流收敛后检查各PV节点无功出力,若越限则将其改为PQ节点并重新求解,直到没有新的越限节点为止。这个过程可能反复进行,因为一个节点转换后会引起其他节点无功重新分配。

三、雅可比矩阵奇异性与电压崩溃点

电压崩溃点(也称临界点、鞍结分岔点)的数学标志是:潮流方程对状态变量的雅可比矩阵奇异,即 det(J)=0。此时系统对负荷增长的灵敏度趋于无穷,微小的负荷增加会导致电压大幅下降。

从物理上理解,J 奇异意味着系统失去了维持电压稳定的能力。在 P-V 曲线上,这对应曲线的鼻点(nose point)。鼻点上半支是稳定运行区,下半支是不稳定区。实际系统运行在上半支,当负荷增加到鼻点时,系统到达电压稳定极限。

需要特别注意的是:J 奇异并不等同于潮流不收敛。在接近奇异时,牛顿法可能因数值条件恶化而不收敛,但真正的临界点是 J 行列式为零的那个精确状态。因此,仅靠常规潮流是否收敛来判断电压稳定极限是不够精确的,这也是连续潮流法被提出的原因。

P-V曲线与雅可比矩阵奇异性关系 负荷 P (p.u.) 电压 V 鼻点 / 临界点 det(J)=0 稳定运行区 J 非奇异,牛顿法收敛 不稳定区 J 奇异后解分支 负荷逐步增加 → 连续潮流:预测-校正 沿切线方向预测,再校正回曲线 PV→PQ 节点转换 Q 越限后失去电压支撑

四、连续潮流法初步:绕过临界点追踪完整曲线

常规牛顿潮流在接近鼻点时,雅可比矩阵接近奇异,迭代难以收敛,无法直接算出临界点。连续潮流法通过引入一个连续参数 λ(负荷增长因子),把潮流方程改写为:

F(V, θ, λ) = 0

并增加一个参数化方程,使方程组在临界点处仍然非奇异。基本步骤是:

  1. 预测步:在当前解处计算切向量,沿切线方向前进一个步长,得到预测解。
  2. 校正步:以预测解为初值,用牛顿法求解扩展潮流方程,得到曲线上的新解。
  3. 步长控制:根据曲率自适应调整步长,在平坦区用大步长,在鼻点附近用小步长。

这样就能从稳定支连续追踪到鼻点,再进入不稳定支,完整绘制 P-V 曲线。鼻点处 λ 达到最大值 λmax,对应的负荷水平就是电压稳定极限。实际分析中,λmax 与当前运行点 λ0 的差值,就是系统的电压稳定裕度。

连续潮流法给出的 λmax 是判断电压稳定极限的直接指标。工程上常要求 λmax/λ0 大于某一阈值(如 1.1~1.2),以保证足够的稳定裕度。这个裕度比单纯看潮流是否收敛更有意义。

五、衔接要点与实操判断

把上述内容串起来,电压稳定分析与潮流计算的衔接可以归纳为三个层次:

层次 潮流计算中的表现 电压稳定含义
正常运行 J 非奇异,牛顿法快速收敛 系统有足够电压支撑
PV→PQ转换 节点类型改变,J 结构变化 无功支撑减少,稳定裕度下降
接近临界点 J 接近奇异,迭代次数增加 电压稳定裕度不足
临界点 det(J)=0,常规潮流不收敛 电压崩溃点,λ=λmax

在实际分析中,建议按以下顺序操作:先建立网络模型和节点分类,运行基础潮流确认收敛;然后逐步增加负荷或传输功率,每步检查PV节点无功是否越限,越限则转换节点类型并重新求解;当常规潮流开始不收敛时,改用连续潮流法继续追踪,直到找到 λmax。最后计算稳定裕度,判断是否满足运行要求。

需要强调的是,雅可比矩阵奇异是电压崩溃的数学标志,但实际系统中电压崩溃往往还伴随无功不足、OLTC动作、发电机励磁限制等动态因素。潮流层面的分析给出的是静态电压稳定极限,是进一步动态分析的基础。理解牛顿法、PV/PQ转换、J奇异性和连续潮流之间的逻辑关系,才能真正把潮流计算工具用于电压稳定评估。