3.3 快速解耦法与直流潮流
3.3 快速解耦法与直流潮流
本节摘要：高压输电网的两个物理事实——电抗远大于电阻、相角差通常很小——让潮流方程的雅可比矩阵可以大幅瘦身。快速解耦法利用 P-θ、Q-V 的天然弱耦合，把一个大型修正方程拆成两个小方程并令雅可比常数化；直流潮流更进一步，把问题线性化到一次矩阵乘法就能出解。本节讲清两者的近似条件、速度收益与适用边界，并给出用直流潮流做 N-1 快速扫描的典型用法。
物理事实换计算红利
先确认两个"物理事实"在数量上是否站得住。220 kV 及以上架空线路的电抗与电阻之比普遍在 5 到 15 之间（回顾 1.3 节：x₁ ≈ 0.4 Ω/km、r₁ ≈ 0.1 Ω/km 以下）；正常运行时相邻节点相角差多数在 10 度以内（0.17 弧度以下），于是 sinθ ≈ θ、cosθ ≈ 1 的近似误差在百分位以下。把这两条代入 3.2 节的雅可比四分块，可以验证：有功失配对相角的敏感度 H 远大于交叉项 N；无功失配对电压的敏感度 L 远大于交叉项 J。换句话说，有功主要跟着相角走、无功主要跟着电压走——P-θ 与 Q-V 天然解耦。
这不是数学家的白日梦，而是电网自己的秉性：有功流动由相角差驱动（可以想象成旋转的"落差"），无功流动由电压幅值差驱动（电压"水位差"）。工程口语"有功管频率、无功管电压"说的就是同一件事的运行侧面。
快速解耦法：两步瘦身
在上述物理事实下做三步简化：丢掉雅可比交叉块 N 与 J（P 方程只留对 θ 的导数，Q 方程只留对 U 的导数）；H、L 块再取 U ≈ 1 的进一步近似，使其只依赖导纳矩阵的虚部；简化后的矩阵用 B'、B' 记（它们由支路电抗倒数直接组装），并在迭代中固定不变。得到迭代格式：
B′·Δθ(k) = ΔP(k)/U(k) B″·ΔU(k) = ΔQ(k)/U(k)
两个小方程交替迭代。收益有三：每个修正方程规模减半；B′、B″ 常数化，三角分解一次完成后全程复用，每轮只剩前代回代；雅可比不再每轮重算。对比 N-R 每轮都要组装并分解雅可比的开销，解耦法单轮成本大幅下降，虽然迭代轮数多一些（五到十轮对三到五轮），总时间通常仍占优。这就是它长期作为生产环境主力算法的原因。
适用边界与 N-R 相反相成：高压输电网（x/r 大）如鱼得水；配电网与电缆网 x/r 比小，解耦假设松动，可能出现收敛慢或不收敛——此时老实回到 N-R。记住这个分工：输网用解耦、配网用牛拉，多数商业软件会自动切换。
直流潮流：线性化的极限形态
把近似推到底——U 全取 1、忽略电阻与无功、sinθ ≈ θ——潮流方程塌缩成一个线性系统：
P = B·θ （B 为节点电纳矩阵） 支路潮流 Pij = (θi − θj)/xij
形式与直流电路的欧姆定律加基尔霍夫定律完全同构（功率对应电流、相角对应电位、电抗对应电阻），故得名"直流潮流"。它丢掉了电压幅值与无功的全部信息，但换来极致速度：一次稀疏矩阵求解，万节点系统毫秒级出解。
它的高频用武之地有两个。其一，N-1 安全校核的批量扫描：逐条开断支路重解线性系统，快速筛出可能过载的严重故障集，再对幸存者上精确潮流复核（第 6 章 6.3 节的流程里它是第一道粗筛）。其二，电力市场的阻塞定价与调度计划校核：市场出清要在几秒内评估成千上万种潮流场景，只有线性模型跑得动。一句要刻在脑子里的边界话：直流潮流只能回答"线路堵不堵、相角偏多少"，永远不要用它回答"电压够不够、无功缺不缺"。
演练：三节点直流潮流手算
沿用 3.2 节的三节点布局（平衡节点 1，P₂ = 0.5 注入，P₃ = −0.9 注出），电抗 x12 = 0.24、x23 = 0.105。列线性方程：节点 2、3 的相角为未知（θ₁ = 0 参考）：
节点2：0.5 = (θ₂−θ₁)/0.24 + (θ₂−θ₃)/0.105 节点3：−0.9 = (θ₃−θ₂)/0.105
由第二条 θ₃ = θ₂ − 0.0945，代回第一条：0.5 = θ₂/0.24 + 0.0945/0.105，解得 θ₂ ≈ 0.0985 弧度（5.6 度）、θ₃ ≈ 0.0040 弧度。支路潮流：P12 = (θ₁−θ₂)/0.24 ≈ −0.41（节点 1 向 2 送 0.41），P23 = (θ₂−θ₃)/0.105 = 0.9？验证一下功率平衡：节点 2 收到 0.41、发出 0.5，恰好向 3 送 0.9；节点 3 收到 0.9 与其负荷吻合。账平了。对照 3.2 节 N-R 结果（θ₂ 收敛于 0.041 附近，因 PV 电压设定与电阻影响），直流潮流的相角偏大一些——这正是忽略电压与电阻的系统性偏差，但它给出的"功率往哪边流、哪条线重"的图景完全正确。粗与细，各司其职。
再引入一个市场与安全分析共用的工具概念：功率转移分布因子 PTDF。它回答"节点 A 向节点 B 转移一兆瓦交易，某条线路的潮流变化多少"——由直流潮流的矩阵一次算出全网的转移系数表。有了 PTDF，N-1 粗筛可以不用逐个开断重解，而是用开断分布因子（LODF）直接由基态潮流推算事故后潮流，速度再上一个台阶；市场里评估双边交易对阻塞的影响也用它。PTDF 是"直流潮流线性化"这棵树上结出的最有工程味道的果子——线性模型的所有系数都可以预先算好、反复使用，这是非线性模型给不了的奢侈。
追问与辨析
直流潮流误差有多大？相角误差通常在一到五度以内、支路有功误差在一成上下，重载或高电阻网络更大——作为筛选用足够，作为结算与决策用则不行。市场结算用直流潮流的准确性争议在业界长期存在：阻塞收益按直流潮流分摊可能与实际物理分布有偏差，所以现货市场普遍要求"出清用直流、结算校核用交流"。直流潮流为什么不考虑电压？不是忘了，是刻意的：它把电压幅值钉死在标幺一，丢掉全部无功信息换取线性——代价与收益都极端，适合回答的问题自然受限。理解一个近似模型的正确姿势，是清楚它"看不见什么"。
本节要点回顾
- 解耦的物理依据：x/r 大、相角差小，使有功随相角、无功随电压，交叉项可弃。
- 快速解耦三收益：方程减半、雅可比常数化、分解结果全程复用；边界是 x/r 比小的配电网。
- 直流潮流把潮流压成线性系统 P = Bθ，速度极快、只答有功分布，用于 N-1 粗筛与市场校核。
- 工具匹配问题：精度要求高用 N-R，输网生产用解耦，海量场景扫描用直流。
求解器的工具箱齐了。下面两节转身面向结果：算完之后，看损耗、看电压。