基于隐式Zbus高斯法的配电网三相不平衡潮流计算程序实现与验证

配电网三相不平衡潮流计算这个方向，做电力系统的人基本都绕不过去。尤其现在分布式光伏、充电桩、农村单相长线路一多，三相不平衡早就不是“偶尔发生”的工况，而是配电网的常态。我见过不少工程师拿单相潮流程序去算三相配网，算出来的节点电压看着还行，但一接上实测数据就完全对不上，问题就出在模型把三相天然的不平衡特性给抹掉了。这些年我用Matlab实现了多套配电网潮流算法，包括前推回代、牛顿法，还有标题里提到的隐式Zbus高斯法，实测下来，隐式Zbus高斯法是“稳定性”和“实现成本”之间权衡得最舒服的一条路，特别是在变压器数量、位置、绕组方式需要灵活设定的场景下，它的优势非常明显。

这篇文章就围绕这套基于隐式Zbus高斯法的三相不平衡潮流计算程序展开，重点讲四件事：为什么配电网潮流必须用三相模型、隐式Zbus高斯法的迭代原理、变压器模型怎么参数化设定（数量、位置、绕组方式），以及如何把IEEE 12、33、36节点这几套经典测试系统跑通并验证结果。适合正在做配电网潮流、分布式电源接入评估、电压质量分析的朋友参考，也适合刚入门想找一套可靠三相潮流代码当起点的研究生。

1. 先搞清楚：配电网的三相不平衡到底是怎么来的，为什么单相潮流不够用

配电网和输电网的物理特性差异非常大。输电网一般是环网结构，三相线路换位，负荷相对对称，用单相模型做等值计算误差不大；配电网则基本是辐射状树形结构，分支多、供电半径长，而且大量负荷是单相接入的。居民负荷、路灯、单相光伏逆变器、空调压缩机，这些都是典型的单相负荷。三条相线上挂的负荷功率根本不可能一样，哪怕总容量设计得比较均衡，实际运行中三相电流也会明显不对称。

这种不平衡会带来什么后果？最直接的是中性点电压偏移。三相四线制系统里，三相负荷越不对称，中性点电流越大，中性点的电位就越偏离零。负荷端的相电压会出现一相偏高、一相偏低的情况，偏低的可能跌破设备正常工作范围，偏高的又可能加速设备绝缘老化。另一条看不见的路径是变压器零序电流：一部分零序分量会通过变压器的连接组别流入接地系统，造成变压器附加发热和中性点过流。

单相潮流为什么算不了这个？因为单相模型把三相负荷给平均化处理了，等于默认三相完全对称。它只能得到一个“等效功率”的电压结果，算不出A相电压偏低、C相偏高的这种真实状态，更算不出中性线电流和零序电压。哪怕你只关心电压合格率，三相不平衡场景下单相法的结果都可能是过乐观的。

另外，配电网还有一个和输电网显著不同的电气特性：R/X比值很大。架空线路和电缆的电阻接近甚至超过电抗，这意味着潮流的收敛特性和有功无功的耦合关系都跟输电网很不一样。牛顿法在多相不平衡配电网里如果初值给得不好，非常容易发散；前推回代法在纯辐射状网络里效率很高、也确实常用，但遇到弱环网和比较复杂的变压器连接方式时，实现起来就比较绕。隐式Zbus高斯法走的是“建导纳阵、形成阻抗阵、反复回代”的路线，不依赖雅可比矩阵，迭代稳定性好，处理辐射状网、少环网的能力都比较自然。这也是我在这个项目里最终没有沿用前推回代、而是把隐式Zbus高斯法作为内核的原因。

2. 隐式Zbus高斯法的核心原理：把非线性负荷问题变成反复求解线性网络

很多第一次接触隐式Zbus高斯法的朋友会问：Zbus不是把网络等值成一个阻抗矩阵嘛，可潮流方程明明是非线性的，一个阻抗矩阵能解决什么事？这个问题问到点子上了。它确实不能一次性解决非线性方程，但可以把非线性问题拆成两步反复交替：线性网络部分交给Zbus，非线性负荷部分每次迭代前单独处理。

配电网里除电源节点外，负荷节点的注入电流和电压是强耦合的。对恒定功率负荷，注入电流表达式是：

I_i = conj(S_i / V_i)

这里的S_i是节点i的负荷复功率，V_i是节点电压。电流取决于电压，电压又由电流决定，这是典型的环比关系，只能迭代逼近。

隐式Zbus高斯法的推导逻辑其实很直接。把电网写成节点电压方程：

Y V = I

两端同时左乘Y的逆：

V = Z I

这里的关键技巧是：Z矩阵只需要形成一次并做一次因子分解。因为网络的拓扑和参数在整个迭代过程中是不变的，变化的只有节点注入电流I。所以每次迭代只需要做两件事：按当前电压更新负荷注入电流I(V^(k))，然后回代求解一次线性系统得到新的电压V^(k+1)。

叠代格式可以写成：

V^(k+1) = Zbus × I(V^(k)) + V_slack

其中V_slack是源母线（一般是变电站母线或平衡节点）的已知电压，Zbus反映的是源母线之外各部分网络对注入电流的响应。每次迭代后判断电压更新量的最大值是否小于收敛容差：

max |V^(k+1) - V^(k)| < ε

如果是就退出，否则继续下一次迭代。

这个思路跟牛顿法最大的区别在于不需要雅可比矩阵的更新和因子分解。牛顿法每一轮迭代都要重新计算雅可比矩阵、做LU分解，在有几百个三相节点、每个节点有三相电压和相角的系统里，计算量和代码复杂度都不小。隐式Zbus高斯法在迭代格式上更像“反复修正电压直到负荷电流自洽”，只要网络是线性的——变压器、线路导纳、电容组这些都不随运行点变化——响应矩阵始终保持不变，所以迭代非常稳。

当然，它也有短板。最大的问题是收敛速度通常比牛顿法慢，尤其在网络规模变大、负荷特别重的情况下，可能迭代几十次甚至上百次。再有就是对非常接近电压崩溃边界的运行点，恒定功率负荷的电流模型会导致迭代不收敛，这是电力系统潮流计算共有的难点。但就配电网常规工况来说，隐式Zbus高斯法的表现相当能打。

如果系统规模很大，别直接对完整的三相导纳矩阵求逆，用稀疏分解再回代会更高效。完整的Zbus显式矩阵在节点很多时会非常占内存——一个三相N节点系统，Zbus维度是3N×3N，节点数上千之后就不是闹着玩的了。我的代码里对中小规模算例直接做了inv，主要是为了简化阅读，工程化的版本建议改成对Ybus做LU分解，然后用“分解一次、反复回代”的模式。

3. 变压器建模：数量、位置、绕组方式在代码里怎么落地

标题里特意标了“可设定变压器数量、位置、绕组方式”，这句话听起来简单，实际实现时工作量不小。很多现成的配电网潮流程序把变压器当成固定环节，预先算好参数塞进数据文件，想改一台变压器的位置或者换一种连接组别，得手动改半天。而一个通用的三相配电网潮流程序，变压器必须是数据驱动的参数化模型。

3.1 变压器在潮流模型里的基本等值结构

配电变压器在主电路中起到的是电压变换和电气隔离的作用。从潮流计算的角度看，一台双绕组变压器可以等值为“理想变压器 + 串联阻抗”的串联组合，其中串联阻抗就是短路阻抗（漏抗为主，电阻占一小部分）。理想变压器负责电压变比和相位移，串联阻抗则产生压降和损耗。有的模型还会在电源侧并联一个励磁支路，但绝大多数配电网潮流计算会把励磁电流忽略掉，因为励磁导纳对潮流分布的影响不到0.5%，对结果基本没影响。

在相分量模型里，变压器需要表达成一个3×3的序导纳矩阵或者相导纳矩阵。三相变压器的不同连接组别，对三个序分量的阻抗响应是完全不同的，尤其零序阻抗差异巨大。下面这个表是设计代码时必须搞清楚的核心差异：

为什么零序这么关键？因为三相不平衡的本质就是零序分量不为零。正序和负序分量在变压器里几乎都能顺利穿越，唯独零序分量对连接组别高度敏感。Dyn11是目前国内配电变压器用得最多的连接组别，低压侧星形带中性线，高压侧三角形。这种结构的好处在于：低压侧的不平衡电流可以在高压侧三角形绕组内部形成零序环流，不至于把零序分量强行注入中压网络，同时低压侧中性点可以做直接接地，让单相负荷获得稳定参考电位。相比之下，Yyn0连接组别的变压器，低压侧虽然也有中性线，但由于高压侧没有三角形绕组提供零序环流回路，零序磁通只能走铁心和油箱等路径，零序阻抗大，中性点电压偏移问题就突出得多。这就是为什么同样的三相不平衡负荷，挂在Dyn11变压器下和挂在Yyn0变压器下，运行表现会差很远。

3.2 代码里如何参数化变压器

我在程序里给变压器定义了一个独立的数据结构，每条记录包含以下字段：

• 所在支路：连接的首末节点编号
• 额定容量（kVA）
• 电压变比（高压侧/低压侧额定电压）
• 短路阻抗百分比和短路损耗
• 连接组别（比如"Dyn11"、"Yyn0"）
• 分接头位置或固定变比修正系数
• 所属馈线/子系统标识

之所以把“数量、位置、绕组方式”单独拎出来说，因为这些参数在传统算例里通常是被写死在代码里的，而在这个程序里全部变成可配置项。多台变压器串联或并联、变压器放在馈线首端还是中段、用Dyn11还是Yyn0，只需要修改数据输入，不需要动迭代引擎的代码。这与真实配电网的规划业务非常贴合——电网改造、新装配变、调整供电范围，对应的就是增删变压器记录、修改连接节点或改绕组方式。

变压器支路加入三相导纳矩阵时，需要根据连接组别构造对应的导纳子块。以常见的Dyn11为例，高压侧三角形绕组、低压侧星形绕组，向量组别的相位移是-30度（不同国标定义可能带正负号，需要注意统一惯例）。在Matlab里按相分量展开时，实际上是把理想变压器的变化关系结合漏阻抗，形成一条三相等值支路。若一次侧相电压向量为Va、Vb、Vc，二次侧为va、vb、vc，理想变化比是k，则等值关系可以写成：

va = k × Va（考虑相位移角度）

漏阻抗Zt折算到某一侧后，串联在理想变压器之后。最终支路导纳矩阵Yt的定义满足：

[I_primary; I_secondary] = Yt × [V_primary; V_secondary]

这部分展开矩阵手工推算容易出错，建议用符号推导或者矩阵分块计算，我代码里写了一个专门的build_transformer_matrix_of_connection函数，按不同连接组别返回6×6的导纳矩阵，再装配到全局相导纳矩阵对应的位置。

4. Matlab代码实现：从数据输入到迭代求解的整体框架

这套代码我拆成了五个模块：数据读取、导纳阵构建、源母线初始化、迭代求解、结果输出。每个模块独立成函数，后续要把前推回代法、改进算法或分布式电源模型接进来，都不需要大动干戈。

4.1 数据结构设计：用什么方式装下三相系统

做三相潮流，第一件容易翻车的事是用单相潮流的数据结构硬套。单相潮流里每个节点只有一个电压，而三相潮流每个节点有三个电压分量，还可能有中性点。我的做法是每个节点用一个结构体（或表格行）存储三相电压的实部和虚部，比如V_phase = [Va_re Va_im; Vb_re Vb_im; Vc_re Vc_im]，尺寸为3×2或直接使用复向量3×1。所有节点拼起来后，节点电压整体是一个3N×1的复向量，导纳矩阵Ybus的维度也是3N×3N。

负荷数据同样要按相存储。一个节点既可能挂三相对称负荷，也可能只挂A相单相负荷，甚至A-B相间负荷。程序里负荷的输入接口设置为三列复功率数组，每一列对应A、B、C三相，缺相的填0。前处理要做一道校验：把所有负荷功率逐相累加，看看每相的总功率是否合理，这能在源头上防止三相数据填错导致结果离谱。

4.2 导纳矩阵构建和稀疏处理

导纳矩阵的构建是从零开始的：先初始化全零矩阵，然后逐条遍历线路支路和变压器支路，把支路导纳子块叠加到全局矩阵的对应位置。线路在三相模型中要考虑相间互感和电容，配电线路距离短，多数情况下可以忽略对地电容，但相间阻抗矩阵不能忽略。带中线或保护线的四线制线路，如需精确计算中性点电压偏移，可以把中性线也作为独立节点建模，矩阵维度会进一步增加。

Matlab里有个非常重要的性能细节：用稀疏矩阵构建导纳阵。先估计非零元素数量，用spalloc分配内存，再填写元素。这个方法在30节点及以下的测试系统上体验不出差异，但一旦扩展到几百节点，效果天差地别。全矩阵方式在1000个三相节点时（矩阵维度3000×3000）已经明显吃力，而稀疏矩阵几乎毫无压力。

4.3 主迭代循环：核心代码逻辑

整个求解过程的主循环大概长这样：

% 初始化电压，源母线设为额定电压，其余节点设为1.0 p.u. V = ones(3*N, 1); V(slack_idx) = 1.0; for iter = 1:max_iter % 根据当前电压计算负荷注入电流 I_load = zeros(3*N, 1); nonzero_load = find(load_active); % 只处理有负荷的节点 I_load(nonzero_load) = conj(S_load(nonzero_load) ./ V(nonzero_load)); % 源母线注入电流单独处理（通常为松弛节点，不需要显式给出） I_load(slack_idx) = 0; % 用Zbus回代求新电压 V_new = Zbus * I_load; V_new(slack_idx) = V_slack; % 收敛判断：电压最大偏差 err = max(abs(V_new - V)); V = V_new; if err < tol converged = true; break; end end

这个版本的代码在概念上是最朴素的隐式Zbus高斯法，适合教学和验证算例。实际工程中可能需要做两点改进：一是把直接求逆改成LU分解，二是为恒功率负荷加阻尼系数，防止某些重负荷场景下出现振荡。

阻尼迭代是实践中非常好用的技巧。标准的迭代格式是直接用新电压替换旧电压，但在负荷重或系统弱的情况下，可能出现电压在两个值之间反复横跳。此时改用：

V^(k+1) = α × V_zbus + (1 - α) × V^(k)

其中α取值0.3~0.7，迭代就会稳定下来。代价是收敛速度降低，但稳定压倒速度。这个技巧我在多个系统上验证过，是解决不收敛问题的最简单手段。

4.4 变压器参数配置示例

看一下参数文件的组织思路，比文字描述直观得多：

% list of transformers % bus_from, bus_to, rating_kva, vratio, zpercent, conn, tap tfm_data = [ 1, 2, 2000, 10/0.4, 4.5, 'Dyn11', 1.0; 4, 7, 800, 10/0.4, 4.0, 'Yyn0', 1.0; 9, 10, 500, 10/0.4, 4.0, 'Dyn11', 0.975; ];

每一行就是一台配电变压器，数量随意增减，位置就是bus_from和bus_to，绕组方式写在conn列。分接头位置通过tap列修正，例如9-10号节点变压器的tap=0.975，相当于低压侧电压比额定值高2.5%，用来补偿线路长距离压降。在程序内部，这个参数会转换为理想变压器的电压变比修正，从而影响潮流结果。

强调一点：分接头设置不是瞎调的。它的本质是牺牲某一侧的电压水平换取另一侧的电压改善。线路末端电压偏低，适当上调变压器低压侧电压是常规手段，但上调太多可能造成首端用户过电压，这是个需要结合全网电压约束来反复试算的工程问题。

5. IEEE 12、33、36节点算例：数据组织、结果验证与收敛性排查

5.1 三个测试系统的特点和差异

IEEE标准算例的典型价值在于提供了可复现的基准结果，很多学术论文都用它做验证。三个系统各有侧重：

• IEEE 12节点：节点规模小，结构简单，适合作为代码调试的第一站。我一般在开发新功能时先用它跑通，输出结果肉眼检查，确认没有低级错误后再上复杂度高的系统。
• IEEE 33节点：这是配电网潮流领域最经典的测试系统之一，采用12.66kV基准电压，含33个节点、32条支路，系统总负荷约3.7MW。大量文献都报告过它的电压分布曲线，非常方便做横向对比。
• IEEE 36节点：规模和复杂度进一步提升，数据通常在配套的测试系统文档中给出，包含更多分支和变压器环节，适合考核程序在处理多变压器场景下的表现。

这三个系统串起来正好是“调试→验证→压力测试”的完整路径。如果程序在12节点系统上结果正确、在33节点系统上和基准结果吻合、再在36节点系统上收敛顺利，那这个程序的基本可靠度就相当高了。

5.2 怎么判断算出来的结果是对的

验证潮流结果的正确性，不能只看“迭代收敛了”就完事，而要从电气物理规律上交叉检查。我常用的检查手段有几个：

第一，功率平衡检查。把源母线输出的总功率算出来，减去全网负荷和网损，应该严格等于零。三相总功率平衡能对上的话，网络模型大概率没有建错。

第二，电压分布趋势检查。辐射状配电网沿馈线方向的电压应该呈阶梯式下降，除非中间有变压器升压或分接头补偿。尤其是重负荷节点附近，电压必然偏低。

第三，和已知基准结果对比。IEEE 33节点系统在标准负荷条件下，典型最小电压出现在节点18附近，大概在0.90~0.92p.u.之间。如果算出来末端电压反而升高，或者最低电压出现在馈线首端，一定哪里出了问题。

第四，三相电压的相位关系检查。正常工况下，三相电压应该大致保持120度相位差，只有严重不平衡或中性点偏移时才会出现明显畸形。如果三相互补关系完全错乱，需要回头检查负荷数据是否填对了相别。

下表列一个IEEE 33节点系统在典型负荷状态下的结果示例，方便对比：

注意这里的数值是假设负荷三相基本平衡但略有偏差时的情况。三相不平衡程度越大，各相之间的电压差异就越明显，这也是检验程序是否真正实现三相建模的试金石：如果你把A、B、C三相负荷输入一组完全一样的数据，三相电压也应该完全相等，否则就说明程序里存在相序处理bug。

5.3 常见的收敛问题与定位思路

迭代不收敛是使用这类程序时最常碰到的问题，我把实际调试中积累的排查顺序说一下。

优先检查导纳矩阵是否奇异。导致奇异最常见的原因是线路数据重复、变压器支路连接错误、或者某些节点变成了孤立节点。用rank函数检查矩阵是否满秩（实操中对角占优的网络这块不至于出问题，但数据少填一行就可能产生断网）。

接着检查源母线设置。隐式Zbus高斯法要求至少有一个松弛节点，否则求解就不具备唯一性。如果代码把源母线也当成普通节点处理，V_new = Zbus * I_load之后整条网络的电压会产生整体漂移，表现为迭代一直不收敛或者结果极其离谱。

然后是负荷功率的极性。潮流里约定流入网络的负荷功率为正，源母线输出功率也为正。如果把负荷功率符号搞反，相当于把负荷当成发电机，电压结果可能不降反升。这个低级错误实际出过的概率非常高，因为很多配电数据表里的功率方向习惯跟课本不完全一致。

最后再考虑阻尼和阈值调整。如果前三种情况都排除了，把阻尼系数从1.0降到0.5试一次，很多顽固的振荡问题可以解决。收敛容差没必要设到1e-10，通常1e-6已经完全满足工程精度，迭代次数在弱网情况下能差出好几倍。

6. 一版跑通之后，这套程序还能往哪些方向扩展

当隐式Zbus高斯法的内核和变压器参数化模型都稳定之后，程序的扩展性会显得非常舒服。至少这几个方向我已经自己实测过，效果都还不错。

第一，接入分布式光伏和储能模型。分布式电源在三相配电网潮流里一般当作PQ节点（功率已知为负值）或PV节点（电压幅值限定）来处理。对PV节点，要在每次迭代中修正无功功率，使电压幅值维持在设定值附近。这个逻辑在隐式Zbus高斯法中并不难加，只要在电流注入计算前额外增加一个无功修正环节。不过要提醒的是，配电网里PV节点数量不宜太多，修正逻辑处理不好会引发振荡，这也是很多光伏并网分析文件里把逆变器切成PQ模式的原因。

第二，弱环网处理。前推回代法对付纯辐射状网很舒服，但如果网络里存在合环开关或者用户自备线路构成闭环，就需要额外的回路补偿计算。隐式Zbus高斯法没有这个问题，因为Zbus矩阵本身就包含了环网拓扑，形成环网时只是导纳矩阵的装配规律变化，迭代格式完全不用改。这也是我向朋友推荐这套算法的重要原因。

第三，负荷动态模型。本文一直默认负荷是恒定功率，但实际配电网里恒阻抗负荷和恒电流负荷也大量存在，比如电热设备近似恒阻抗、城市照明近似恒电流。比较好的做法是允许每一类负荷设定百分比，例如总负荷的40%为恒功率、30%为恒电流、30%为恒阻抗，然后在每次迭代时把混合负荷转化为等效注入电流。这个功能对算冬季空调负荷、夜间照明负荷特别有用。

第四，中性线显式建模。当前版本的默认设置相当于假定中性点电位足够稳定，但真要分析严重的单相接地或者三相严重失衡导致的中性点电压偏移，最好把中性线作为独立的导纳支路建进去。代价是矩阵维度增加比较多，但从结果对比来看，中性线建模对低压用户侧电压质量的分析价值是很大的。

第五，运行点连续扫描。配电网负荷日变化曲线通常按小时或15分钟一条点，如果想算全日电压曲线，可以把潮流程序包在循环里，逐点计算。隐式Zbus高斯法每步迭代主要开销就是一次回代，几百个点时长的跨度完全在可接受范围内。配合这套代码的变压器参数化设定，甚至可以做“不同分接头设置下全年电压合格率”这种规划级分析。

我个人在实际工程中体会最深的一点是：三相不平衡潮流程序最容易出错的地方往往不是数学内核，而是三相数据组织和变压器连接方式的处理。一套干净的数据结构、一套参数化变压器配置、再加上一个不依赖良好初值的稳定迭代内核，这套代码就能真正脱离论文Demo的范畴，变成日常分析工具箱里的一员。如果你正卡在配电网三相潮流分析的起步阶段，沿着隐式Zbus高斯法这条线往下走，是一条性价比很高的路线。