4. 潮流计算原理与ETAP实现
潮流计算是电力系统分析中最基础、使用频率最高的计算任务。它要回答的核心问题是：在给定网络拓扑、线路与变压器参数、发电机出力（PV节点电压幅值与有功）以及负荷功率（PQ节点有功与无功）的条件下，各节点的电压幅值与相角、各支路的功率分布和网络损耗分别是多少。ETAP的潮流计算模块（Load Flow）正是围绕这一非线性方程组求解展开的，理解其算法原理，才能正确设置参数、判断收敛并合理解读结果。
潮流方程与节点分类
对任一节点 i，注入功率可写为电压与导纳矩阵的二次型：
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 为节点导纳矩阵的实部与虚部。由于方程对 V 和 θ 是非线性的，必须迭代求解。节点按已知量分为三类：平衡节点（Slack）已知 V 和 θ，承担全网功率平衡，通常选主电源母线；PV节点已知 P 和 V，未知 Q 和 θ，对应带自动电压调节的发电机母线；PQ节点已知 P 和 Q，未知 V 和 θ，对应负荷母线。未知量总数等于方程数，这是潮流可解的前提。
高斯-赛德尔法与牛顿-拉夫逊法
高斯-赛德尔法把节点方程改写为 V_i 的显式形式，逐节点代入最新电压值反复扫描。它实现简单、内存占用小，但收敛速度是线性的，迭代次数随系统规模迅速增加，且对重载、弱环网可能不收敛，现代工程中多用于教学演示或作为牛顿法的初值生成手段。
牛顿-拉夫逊法将非线性方程在当前点线性化，构造雅可比矩阵 J，求解修正量 Δx：
[Δθ; ΔV/V] = −J⁻¹ · [ΔP; ΔQ]每次迭代重新计算不平衡量 ΔP、ΔQ 和雅可比矩阵，再更新电压。它具有平方收敛特性，通常 3～5 次迭代即可达到 10⁻⁵ 量级的功率不平衡，是ETAP默认且最常用的算法。代价是每次迭代需形成并求解雅可比矩阵，规模大时计算量较高，因此ETAP还提供快速解耦法（Newton-Raphson Fast Decoupled），利用高压电网中 R≪X 的特点，将 P-θ 与 Q-V 解耦并固定雅可比近似，牺牲少量收敛性换取每轮更快的计算速度，适合大规模输电网的批量扫描。
ETAP潮流计算设置要点
在ETAP中执行潮流前，需在 Load Flow Study Case 中确认以下设置：
- 计算方法（Method）：默认 Newton-Raphson；若网络含大量辐射状馈线或需要快速批量计算，可切换为 Fast Decoupled 或 Adaptive Newton。
- 收敛容差（Precision）：默认 0.0001 p.u.（即 0.01%），一般无需修改；若系统接近临界状态可放宽至 0.001 以观察趋势，但正式报告应保持默认。
- 最大迭代次数（Max Iterations）：默认 50，重载或含长距离弱联系线路时可提高到 100～200。
- 变压器分接头与电压调整：需选择是否考虑 LTC 自动调压、是否按实际分接头位置计算。
- 负荷模型：恒功率、恒电流、恒阻抗或 ZIP 组合，直接影响重载下的电压结果。
- 发电机无功限值：PV节点 Q 越限后是否自动转为 PQ 节点，这是实际运行中常见的处理方式。
收敛调整与常见问题
当ETAP提示“System did not converge”时，通常不是软件问题，而是系统本身接近或超过可行解边界。可按以下顺序排查：
- 检查孤岛与失电区域：未连接的母线或断开的开关会形成无源孤岛，导致方程无解。用ETAP的“Island Detection”确认所有母线均被电源覆盖。
- 检查平衡节点设置：每个电气岛必须且只能有一个 Slack 节点，否则功率无法平衡。
- 检查参数合理性：线路阻抗为零、变压器变比异常、负荷功率过大都会导致雅可比矩阵奇异。
- 调整初值：将母线电压初值设为 1.0 p.u.、相角设为 0°，或启用“Flat Start”。
- 放宽迭代限制：将最大迭代次数提高到 200，观察不平衡量是否单调下降。若持续振荡，说明系统确实无解。
- 切换算法：牛顿法不收敛时可尝试 Adaptive Newton 或 Fast Decoupled，后者对初值不敏感。
结果解读与判断依据
潮流收敛后，ETAP会在 Load Flow Report 中给出母线电压、支路潮流、损耗和告警。解读时重点关注：
| 指标 | 典型判据 | 异常含义 |
|---|---|---|
| 母线电压（p.u.） | 0.95～1.05 | 低于 0.95 可能欠压，需调压或增容 |
| 支路载流率（%） | ≤ 100% | 超过 100% 表示过载，需校核导线或调整运行方式 |
| 网络总损耗 | 与历史值对比 | 异常升高可能提示参数错误或环流 |
| 平衡节点出力 | 与总负荷+损耗匹配 | 偏差大说明存在未计入的电源或负荷 |
| 迭代次数 | 3～8 次 | 超过 15 次说明系统接近临界，需关注 |
以某 5 母线示例系统为例（参数为教学假设值）：基准容量 100 MVA，母线 1 为平衡节点 V=1.05∠0°，母线 2、3 为 PV 节点，母线 4、5 为 PQ 节点。采用牛顿-拉夫逊法，容差 0.0001 p.u.，经 4 次迭代收敛。结果显示母线 5 电压为 0.962 p.u.，线路 3-5 载流率 87%，网络总损耗 2.3 MW。此时电压虽在允许范围内，但已接近下限，若负荷再增长 10%，母线 5 可能跌破 0.95 p.u.，需要提前规划无功补偿或调整变压器分接头。
掌握算法原理与ETAP设置之间的对应关系，是正确使用潮流计算的关键。牛顿-拉夫逊法的高效收敛依赖于合理的初值和完整的雅可比矩阵，而ETAP中的每一项设置——容差、迭代次数、负荷模型、发电机限值——都会直接影响最终结果的物理意义。只有把算法逻辑与工程判据结合起来，才能从一堆数字中读出系统真实的运行状态。