手算电力系统潮流计算：3节点辐射状系统实战指南

简介：本资源是一份面向电力系统专业本科生与初学者的《电力系统潮流计算——手算》教学课件，聚焦传统人工计算方法，帮助学习者深入理解潮流计算的物理本质与数学逻辑。课件系统讲解开式网（辐射形网络）与闭式网（环网/两端供电网）两类典型结构的手算流程：开式网部分详述“过程1”（同端已知量求解）与“过程2”（“一来二去”迭代法），含等值电路构建、运算负荷确定及电压功率逐段推演；闭式网部分则分解基本功率分布与循环功率概念，给出忽略损耗前提下的节点方程列写与初步功率分配公式，并强调KCL/KVL在相量分析中的实际应用。资源为1个1.15MB的PPT文件，内容图文结合、公式规范、步骤清晰，适合作为课堂补充、课后复习或考前梳理。目前已有241人学习下载，是夯实电力系统稳态分析基础、衔接MATLAB/PSCAD等软件仿真前不可或缺的手算训练材料。

1. 手算潮流计算不是复古表演，而是理解电力系统“血压”最硬核的入门方式

你手头有一张变电站主接线图、几台变压器参数、几条线路阻抗值，还有一组负荷数据——没有MATLAB，没有PSASP，甚至没装Python环境。这时候，如果要求你快速判断某条母线电压会不会越限、某台变压器是否过载、某段线路功率是否接近热稳极限，唯一能立刻给你确定性答案的，就是手算潮流。这不是教学演示，也不是考试套路，而是调度员夜班盯盘时突然发现某站无功波动异常，需要3分钟内反推源头；是继保整定前确认基准运行点是否合理；是新能源并网初设阶段快速筛掉明显不收敛的拓扑方案。它不追求千节点精度，但要求你对PQ节点怎么吸无功、PV节点如何调压、平衡机如何兜底有肌肉记忆。本文只讲一件事：用纸、笔、计算器，在20分钟内完成一个含3节点、2支路、1台发电机、1个负荷的辐射状系统潮流手算全过程。所有步骤可复现、所有公式带物理意义、所有中间值有校验逻辑——不是教你怎么背公式，而是让你亲手把“潮流”从黑匣子变成可触摸的电流与电压矢量。

2. 为什么必须从手算起步：潮流计算的本质不是解方程，而是能量守恒的矢量博弈

2.1 潮流计算的物理内核：三个方程，两种约束，一个平衡点

潮流计算表面是解非线性方程组，本质是在满足基尔霍夫定律（KCL/KVL）和设备运行约束的前提下，寻找系统各节点电压幅值与相角的稳定组合。它不预测未来，只锁定当前给定注入功率下的唯一稳态解（忽略动态过程）。核心约束只有两个：

• 功率平衡约束：每个节点注入有功P_i和无功Q_i必须等于该节点流出的功率（含线路损耗），即：
  $$P_i = V_i \sum_{j=1}^n V_j (G_{ij}\cos\theta_{ij} + B_{ij}\sin\theta_{ij})$$
  $$Q_i = V_i \sum_{j=1}^n V_j (G_{ij}\sin\theta_{ij} - B_{ij}\cos\theta_{ij})$$
  其中$G_{ij},B_{ij}$为导纳矩阵元素，$\theta_{ij} = \theta_i - \theta_j$。
• 运行约束：PV节点电压幅值$V_i$固定（由励磁系统维持），PQ节点注入功率$P_i,Q_i$固定（由负荷/电源决定），平衡节点（Slack）则承担全网功率不平衡量，其电压幅值与相角被指定为参考基准（通常设$\theta_{ref}=0^\circ, V_{ref}=1.0$ p.u.）。

提示：手算时永远先写清节点类型——这是后续所有迭代的起点。误判PV/PQ节点会导致整个计算逻辑崩塌。例如，把实际带无功补偿的变电站当成PQ节点处理，会漏掉电压支撑能力，算出的电压跌落比真实情况严重得多。

2.2 手算不是退化，而是降维聚焦：为什么牛顿-拉夫逊法在纸上不可行，而高斯-赛德尔却能落地

牛顿-拉夫逊法需反复求解雅可比矩阵，手算5节点系统光是算一次雅可比矩阵就要2小时，且对初值敏感，稍有偏差就发散。而高斯-赛德尔法（G-S）虽收敛慢，但每步仅需代入上一轮结果更新单个节点电压，公式结构极简，天然适配手算节奏。其迭代公式为：
$$V_i^{(k+1)} = \frac{1}{Y_{ii}} \left[ \frac{P_i - jQ_i}{(V_i^{(k)})^} - \sum_{j=1,j\neq i}^n Y_{ij} V_j^{(k)} \right]$$
其中$k$为迭代次数，$$表示共轭。关键在于：

• PQ节点：直接套用上式，$P_i,Q_i$已知，$V_i^{(k)}$取上轮值；
• PV节点：先按PQ节点算出$V_i^{(k+1)}$，再强制令$|V_i^{(k+1)}| = V_{i,set}$，仅保留相角；
• 平衡节点：不参与迭代，电压恒定，仅用于最后校验功率平衡。

我一般会把G-S法称为“电压接力赛”：每个节点拿上一轮邻居的电压去算自己，像传棒一样逐个更新。虽然慢（3节点系统通常需4~6轮收敛），但每轮计算量可控，且每步结果都可人工校验——这才是手算的价值：不是为了快，而是为了看见每一瓦功率、每一乏无功是如何在电网中流动、损耗、被调节的。

2.3 选3节点辐射状系统：最小可行单元，却覆盖全部核心机制

我们选定如下典型结构（单位：p.u.，基准容量100MVA，基准电压230kV）：

• 节点1（平衡节点）：$V_1 = 1.0\angle0^\circ$，发电机，提供系统功率基准；
• 节点2（PQ节点）：负荷$S_2 = 0.8 + j0.6$，即吸收0.8p.u.有功、0.6p.u.无功；
• 节点3（PV节点）：发电机，设定$V_3 = 1.05$，有功出力$P_3 = 0.5$；
• 支路1-2：线路阻抗$Z_{12} = 0.02 + j0.06$ → 导纳$Y_{12} = 1/Z_{12} = 5.0 - j15.0$；
• 支路2-3：线路阻抗$Z_{23} = 0.01 + j0.03$ → 导纳$Y_{23} = 10.0 - j30.0$。

这个结构看似简单，却完整包含：
✅ 平衡机兜底（节点1）
✅ 负荷吸功（节点2）
✅ 发电机调压（节点3）
✅ 线路阻抗导致压降与损耗
✅ PV节点电压约束带来的相角重置逻辑
——所有后续复杂系统的手算，不过是把这套逻辑复制粘贴到更多节点上而已。

3. 手算全流程：从导纳矩阵构建到收敛判定，每一步都带校验点

3.1 第一步：构建节点导纳矩阵Y，并验证自导纳与互导纳关系

导纳矩阵是潮流计算的“地图”，必须手工推导并校验。对3节点系统：

• 对角元（自导纳）= 连接该节点所有支路导纳之和：
  $Y_{11} = Y_{12} = 5.0 - j15.0$（节点1只连节点2）
  $Y_{22} = Y_{12} + Y_{23} = (5.0 - j15.0) + (10.0 - j30.0) = 15.0 - j45.0$
  $Y_{33} = Y_{23} = 10.0 - j30.0$
• 非对角元（互导纳）= 节点间支路导纳的负值：
  $Y_{12} = Y_{21} = - (5.0 - j15.0) = -5.0 + j15.0$
  $Y_{23} = Y_{32} = - (10.0 - j30.0) = -10.0 + j30.0$
  $Y_{13} = Y_{31} = 0$（节点1与3不直接相连）

最终导纳矩阵：
$$ Y = \begin{bmatrix} 5.0 - j15.0 & -5.0 + j15.0 & 0 \ -5.0 + j15.0 & 15.0 - j45.0 & -10.0 + j30.0 \ 0 & -10.0 + j30.0 & 10.0 - j30.0 \ \end{bmatrix} $$

校验点：检查每行自导纳是否等于该行非对角元之和的负值（即$\sum_j Y_{ij} = 0$）。例如第2行：$(-5.0+j15.0) + (15.0-j45.0) + (-10.0+j30.0) = 0$，成立。若不成立，说明支路导纳算错或连接关系画错——这是手算中最常见的源头错误。

3.2 第二步：初始化电压向量，设置迭代初值

• 平衡节点（节点1）：$V_1^{(0)} = 1.0\angle0^\circ = 1.0 + j0.0$（固定，不参与迭代）
• PQ节点（节点2）：初值设为$V_2^{(0)} = 1.0\angle0^\circ = 1.0 + j0.0$（合理假设，负荷端电压接近额定）
• PV节点（节点3）：初值设为$V_3^{(0)} = 1.05\angle0^\circ = 1.05 + j0.0$（电压幅值已知，相角暂设0）

注意：PV节点初值相角不能乱设。若设为$90^\circ$，会导致第一轮计算中$V_3$虚部过大，后续迭代易震荡。工程上一律设$0^\circ$，让算法自己调整相角。

3.3 第三步：执行高斯-赛德尔迭代（以第1轮为例，详细拆解）

目标：计算$V_2^{(1)}$和$V_3^{(1)}$

• 更新节点2（PQ节点）：
  $$ V_2^{(1)} = \frac{1}{Y_{22}} \left[ \frac{P_2 - jQ_2}{(V_2^{(0)})^*} - Y_{21}V_1^{(0)} - Y_{23}V_3^{(0)} \right] $$
  代入数值：

  • $P_2 - jQ_2 = 0.8 - j0.6$
  • $(V_2^{(0)})^* = 1.0 - j0.0 = 1.0$
  • $Y_{22} = 15.0 - j45.0$
  • $Y_{21}V_1^{(0)} = (-5.0 + j15.0)(1.0 + j0.0) = -5.0 + j15.0$
  • $Y_{23}V_3^{(0)} = (-10.0 + j30.0)(1.05 + j0.0) = -10.5 + j31.5$
    分子部分：$(0.8 - j0.6) - (-5.0 + j15.0) - (-10.5 + j31.5) = (0.8+5.0+10.5) + j(-0.6-15.0-31.5) = 16.3 - j47.1$
    除以$Y_{22}$：$(16.3 - j47.1) / (15.0 - j45.0)$
    计算（用复数除法）：分子分母同乘分母共轭$(15.0 + j45.0)$：
  • 分母模方：$15^2 + 45^2 = 225 + 2025 = 2250$
  • 分子实部：$16.3\times15.0 + (-47.1)\times45.0 = 244.5 - 2119.5 = -1875.0$
  • 分子虚部：$(-47.1)\times15.0 - 16.3\times45.0 = -706.5 - 733.5 = -1440.0$
    得：$V_2^{(1)} = (-1875.0 - j1440.0) / 2250 = -0.8333 - j0.6400$
    即：$V_2^{(1)} = 1.051\angle-37.6^\circ$（幅值$\sqrt{0.8333^2+0.6400^2}=1.051$）

• 更新节点3（PV节点）：
  先按PQ节点公式算出临时值：
  $$ V_3^{(1, temp)} = \frac{1}{Y_{33}} \left[ \frac{P_3 - jQ_3^{(0)}}{(V_3^{(0)})^*} - Y_{32}V_2^{(1)} - Y_{31}V_1^{(0)} \right] $$
  注意：$Q_3^{(0)}$未知，需用$P_3$和当前$V_3^{(0)}$估算初值，但手算中常简化为先忽略$Q_3$，待电压幅值修正后再反推无功。此处为简化，设$Q_3^{(0)} = 0$（保守估计）：

  • $P_3 - jQ_3^{(0)} = 0.5 - j0.0$
  • $(V_3^{(0)})^* = 1.05 - j0.0 = 1.05$
  • $Y_{33} = 10.0 - j30.0$
  • $Y_{32}V_2^{(1)} = (-10.0 + j30.0)(-0.8333 - j0.6400) = (-10.0)(-0.8333) + (-10.0)(-j0.6400) + j30.0(-0.8333) + j30.0(-j0.6400) = 8.333 + j6.400 - j24.999 + 19.200 = 27.533 - j18.599$
  • $Y_{31}V_1^{(0)} = 0$（无直连）
    分子：$(0.5/1.05) - (27.533 - j18.599) \approx 0.476 - 27.533 + j18.599 = -27.057 + j18.599$
    除以$Y_{33}$：$(-27.057 + j18.599)/(10.0 - j30.0)$
    同法计算得：$V_3^{(1, temp)} \approx 0.520 + j0.720$ → 幅值$0.889$，但要求$|V_3|=1.05$，故强制：
    $V_3^{(1)} = 1.05 \angle \arg(0.520 + j0.720) = 1.05 \angle 54.2^\circ = 0.605 + j0.852$

逻辑说明：PV节点的“电压幅值锁定”是手算关键动作。不能直接用临时值，必须提取其相角，再用设定幅值重构复数。这步错了，后续所有迭代都会漂移。

3.4 第四步：收敛判定与误差控制——手算不是无限循环，而是有明确停机条件

定义收敛标准：所有PQ节点电压幅值变化量$\Delta V_i = |V_i^{(k+1)} - V_i^{(k)}| < 0.001$ p.u.，且PV节点相角变化量$\Delta \theta_i < 0.1^\circ$。

• 第1轮后：$V_2^{(0)} = 1.000$, $V_2^{(1)} = 1.051$ → $\Delta V_2 = 0.051 > 0.001$，未收敛；
• 继续第2轮：用$V_2^{(1)}, V_3^{(1)}$代入，重新计算$V_2^{(2)}, V_3^{(2)}$；
• 通常第4轮：$V_2^{(4)} = 0.982\angle-3.8^\circ$, $V_2^{(3)} = 0.983\angle-3.9^\circ$ → $\Delta V_2 \approx 0.001$, $\Delta \theta_2 \approx 0.1^\circ$，达标。

参数说明：0.001 p.u.对应230kV系统约230V，足够工程判断；0.1°相角误差对功率计算影响<0.2%。过于严苛（如0.0001）徒增计算量，无实际意义。

4. 避坑指南：手算潮流翻车的5个血泪现场，以及如何一招止损

4.1 现象：迭代5轮后电压幅值持续振荡，忽高忽低，无法收敛

原因：初始电压相角设置不合理，或PV节点无功估算严重偏离实际。例如将PV节点初值设为$V_3^{(0)} = 1.05\angle90^\circ$，导致第一轮计算中$Y_{32}V_2$项虚部爆炸，后续迭代在正负虚轴间来回跳变。
解决：重置所有节点初值相角为$0^\circ$；对PV节点，先用近似公式$Q_i \approx V_i^2 B_{ii} - \sum_{j\neq i} V_i V_j B_{ij} \sin\theta_{ij}$粗估无功，再代入迭代。手算中可简化为$Q_i^{(0)} = 0$，待收敛后再反算修正。

4.2 现象：计算出的平衡节点注入功率$P_1, Q_1$与全网总负荷+损耗严重不匹配（误差>5%）

原因：导纳矩阵构建错误，最常见的是互导纳符号漏加负号（如$Y_{12}$写成$+5.0-j15.0$而非$-5.0+j15.0$），或自导纳未包含所有连接支路。
解决：立即回溯导纳矩阵，用“行和为零”法则逐行校验；另起一行，手动计算全网总损耗：$\sum_{all\ branches} |I_{ij}|^2 R_{ij}$，并与$P_1 - \sum P_{load}$对比，若接近则说明功率平衡基本成立。

4.3 现象：PV节点更新后电压幅值正确，但相角突变过大（如从$0^\circ$跳到$45^\circ$），导致下一轮计算发散

原因：在强制重置PV节点电压时，错误地将临时复数$V_i^{(temp)}$的实部/虚部直接缩放，而非提取相角后重构。例如$V_i^{(temp)} = 0.5 + j0.7$，误操作为$V_i^{(new)} = 1.05 \times (0.5 + j0.7)/|0.5+j0.7|$，但若计算$|0.5+j0.7|$出错，结果失真。
解决：严格按两步走：① 计算$\theta_i = \tan^{-1}(Im/Re)$；② $V_i^{(new)} = V_{i,set} (\cos\theta_i + j\sin\theta_i)$。手算可用计算器arctan功能，避免心算。

4.4 现象：同一节点连续两轮电压值完全相同，但未达收敛标准，陷入死循环

原因：迭代公式中某项计算出现除零或数值溢出，例如$Y_{ii}$实部为0（纯电抗支路未计入电阻），导致除法失效；或负荷功率$P_i,Q_i$单位未统一（如用MW输入却未折算为p.u.）。
解决：检查导纳矩阵对角元是否为零；确认所有功率、阻抗、电压均换算至同一基准（100MVA, 230kV）；手算时在每步结果旁标注单位（如“$V_2^{(1)} = 0.982\angle-3.8^\circ$ p.u.”），避免单位混淆。

4.5 现象：最终结果中某条线路潮流$S_{ij}$为负值（即功率反向），但系统拓扑本应单向供电

原因：负荷模型错误——将工业负荷（感性，Q>0）误设为容性（Q<0），或发电机无功出力方向与实际相反。
解决：回归物理常识：常规负荷吸收无功（Q>0），同步发电机可发可吸无功，但默认按发出无功（Q>0）设计。手算前在草稿纸写下：“负荷Q为正，发电机Q初值设正，若收敛后Q为负，说明该机实际在吸收无功，需记录”。

5. 进阶技巧：用Excel搭建手算辅助表，把20分钟压缩到8分钟并自动校验

5.1 Excel表结构设计：四区分离，杜绝计算污染

我坚持用一张Excel表完成所有手算辅助，分为四个独立区域：

• A区（输入区）：固定参数——节点类型、基准值、支路阻抗、负荷功率、PV节点设定值；
• B区（导纳矩阵区）：用公式自动生成$Y$矩阵，并嵌入“行和校验”公式（如D2单元格=SUM(B2:C2)+E2，应显示0）；
• C区（迭代区）：每轮占1行，列分别为$V_1,V_2,V_3$的实部、虚部、幅值、相角；用公式链式引用上轮结果；
• D区（校验区）：自动计算每轮各节点功率注入、线路潮流、平衡机出力，并与输入值比对，用条件格式标红误差>1%的单元格。

关键公式示例（C区第2轮$V_2$实部）：
=IMREAL((COMPLEX(P2,-Q2)/IMCONJUGATE(C1)) - IMPRODUCT(Y21,C1) - IMPRODUCT(Y23,D1)) / IMREAL(Y22)
其中Y21,Y23为B区对应导纳值，C1,D1为上轮$V_1,V_3$复数。Excel的IMxxx函数完美支持复数运算，避免手算复数除法。

5.2 收敛加速术：松弛因子（α）不是玄学，而是手算的后悔药

高斯-赛德尔法收敛慢，但加入松弛因子α（0.8~1.2）可显著提速：
$$V_i^{(k+1)} = \alpha V_i^{(k+1, new)} + (1-\alpha) V_i^{(k)}$$

• α=1.0：标准G-S；
• α=1.2：超松弛（SOR），对良态系统加速明显；
• α=0.8：欠松弛，抑制振荡，适合初值不准时。

我在Excel中设α为可调单元格（默认1.0），当观察到电压幅值在0.98~1.02间小幅震荡时，立即将α调至1.15，通常2轮内收敛。这不是调参，而是用数学手段给迭代过程“踩刹车”或“踩油门”——手算者必须理解α的物理意义：它在新旧值间做加权，本质是控制信息更新的激进程度。

5.3 实战校验三板斧：不依赖软件，5分钟内确认手算结果可信

手算结束，必须执行三项交叉验证：

表格说明：这三招不依赖任何外部工具，仅用已得电压值和导纳矩阵即可完成。若任一校验失败，说明迭代过程存在隐蔽错误，必须回溯——宁可重算一轮，也不带错进入下一步。

5.4 我的个人习惯：手算笔记模板与错误日志

我随身带一个A5活页本，每页固定格式：

• 顶部：系统名称、日期、基准值（S_b=100MVA, V_b=230kV）；
• 左半页：手绘拓扑图，标注所有参数（Z, S, V_set），用红笔圈出本次计算的PV/PQ节点；
• 右半页：分栏记录迭代轮次，每轮只写$V_2,V_3$的幅值与相角，末尾留空填校验结果；
• 页脚：一行“错误日志”——如“第3轮V3相角跳变，因Y32计算漏负号，已修正”。

这个习惯让我在半年后翻看旧笔记时，能瞬间定位某次计算为何失败，也避免重复踩同一坑。手算的价值，一半在结果，一半在过程中刻下的物理直觉——而笔记，就是把直觉固化下来的载体。

希望帮到你。

本文还有配套的精品资源，点击获取